The driving mechanisms at the base of the clearance of biological wastes in the brain interstitial space (ISS) are still poorly understood and an actively debated subject. A complete comprehension of the processes that lead to the aggregation of amyloid proteins in such environment, hallmark of the onset and progression of Alzheimer’s disease, is of crucial relevance. Here we employ combined computational fluid dynamics and molecular dynamics techniques to uncover the role of fluid flow and proteins transport in the brain ISS. Our work identifies diffusion as the principal mechanism for amyloid-β proteins clearance, whereas fluid advection may lead transport for larger molecular bodies, like amyloid-β aggregates or extracellular vesicles. We also clearly quantify the impact of large nascent prefibrils on the fluid flowing and shearing. Finally, we show that, even in the irregular brain interstitial space (ISS), hydrodynamic interactions enhance amyloid-β aggregation at all stages of the aggregation pathway. Our results are key to understand the role of fluid flow and solvent-solute interplay on therapeutics like antibodies acting in the brain ISS.

In this work, we investigate how fluid flows impact the aggregation mechanisms of Aβ40 proteins and Aβ16–22 peptides and mechanically perturb their (pre)fibrillar aggregates. We exploit the OPEP coarse-grained model for proteins and the Lattice Boltzmann Molecular Dynamics technique. We show that beyond a critical shear rate, amyloid aggregation speeds up in Couette flow because of the shorter collisions times between aggregates, following a transition from diffusion limited to advection dominated dynamics. We also characterize the mechanical deformation of (pre)fibrillar states due to the fluid flows (Couette and Poiseuille), confirming the capability of (pre)fibrils to form pathological loop-like structures as detected in experiments. Our findings can be of relevance for microfluidic applications and for understanding aggregation in the interstitial brain space.

Langevin and Brownian simulations play a prominent role in computational research, and state of the art integration algorithms provide trajectories with different stability ranges and accuracy in reproducing statistical averages. The practical usability of integrators is an important aspect to allow choosing large time steps while ensuring numerical stability and overall computational efficiency. In this work, different use cases and practical features are selected in order to perform a cumulative comparison of integrators with a focus on evaluating the derived velocity and position autocorrelation functions, a comparison that is often disregarded in the literature. A standard industrial open-source software methodology is suggested to compare systematically the different algorithms.

Macromolecular crowding has profound effects on the mobility of proteins, with strong implications on the rates of intracellular processes. To describe the dynamics of crowded environments, detailed molecular models are needed, capturing the structures and interactions arising in the crowded system. In this work, we present OPEPv7, which is a coarse-grained force field at amino-acid resolution, suited for rigid-body simulations of the structure and dynamics of crowded solutions formed by globular proteins. Using the OPEP protein model as a starting point, we have refined the intermolecular interactions to match the experimentally observed dynamical slowdown caused by crowding. The resulting force field successfully reproduces the diffusion slowdown in homogeneous and heterogeneous protein solutions at different crowding conditions. Coupled with the lattice Boltzmann technique, it allows the study of dynamical phenomena in protein assemblies and opens the way for the in silico rheology of protein solutions.

Mucociliary clearance is the first defense mechanism of the respiratory tract against inhaled particles. This mechanism is based on the collective beating motion of cilia at the surface of epithelial cells. Impaired clearance, either caused by malfunctioning or absent cilia, or mucus defects, is a symptom of many respiratory diseases. Here, by exploiting the lattice Boltzmann particle dynamics technique, we develop a model to simulate the dynamics of multiciliated cells in a two-layer fluid. First, we tuned our model to reproduce the characteristic length- and time-scales of the cilia beating. We then check for the emergence of the metachronal wave as a consequence of hydrodynamic mediated correlations between beating cilia. Finally, we tune the viscosity of the top fluid layer to simulate the mucus flow upon cilia beating, and evaluate the pushing efficiency of a carpet of cilia. With this work, we build a realistic framework that can be used to explore several important physiological aspects of mucociliary clearance.

Polyelectrolytes can electrophoretically be driven through nanopores in order to be detected. The respective translocation events are often very fast and the process needs to be controlled to promote efficient detection. To this end, we attempt to control the translocation dynamics by coating the inner surface of a nanopore. For this, different charge distributions are chosen that result in substantial variations of the pore–polymer interactions. In addition and in view of the existing detection modalities, experimental settings, and nanopore materials, different types of sensors inside the nanopore have been considered to probe the translocation process and its temporal spread. The respective transport of polyelectrolytes through the coated nanopores is modeled through a multi-physics computational scheme that incorporates a mesoscopic/electrokinetic description for the solvent and particle-based scheme for the polymer. This investigation could underline the interplay between sensing modality, nanopore material, and detection accuracy. The electro-osmotic flow and electrophoretic motion in a pore are analyzed together with the polymeric temporal and spatial fluctuations unraveling their correlations and pathways to optimize the translocation speed and dynamics. Accordingly, this work sketches pathways in order to tune the pore–polymer interactions in order to control the translocation dynamics and, in the long run, errors in their measurements.

Abstract Von Willebrand Factor (vWf) is a giant multimeric extracellular blood plasma involved in hemostasis. In this work we present multi-scale simulations of its three-domains fragment A1A2A3. These three domains are essential for the functional regulation of vWf. Namely the A2 domain hosts the site where the protease ADAMTS13 cleavages the multimeric vWf allowing for its length control that prevents thrombotic conditions. The exposure of the cleavage site follows the elongation/unfolding of the domain that is caused by an increased shear stress in blood. By deploying Lattice Boltzmann molecular dynamics simulations based on the OPEP coarse-grained model for proteins, we investigated at molecular level the unfolding of the A2 domain under the action of a perturbing shear flow. We described the structural steps of this unfolding that mainly concerns the beta-strand structures of the domain, and we compared the process occurring under shear with that produced by the action of a directional pulling force, a typical condition of single molecule experiments. We observe, that under the action of shear flow, the competition among the elongational and rotational components of the fluid field leads to a complex behaviour of the domain, where elongated structures can be followed by partially collapsed melted globule structures with a very different degree of exposure of the cleavage site. Our simulations pose the base for the development of a multi-scale in-silico description of vWf dynamics and functionality in physiological conditions, including high resolution details for molecular relevant events, e.g., the binding to platelets and collagen during coagulation or thrombosis.

We present a numerical investigation of the airflow dynamics and particle transport through an averaged human nasal cavity. The effect of particle size and breathing rate on the deposition patterns are explored. The simulations reveal that smaller particles penetrate deeper into the airway, whereas larger particles agglomerate near the anterior portion of the nasal cavity. Increasing the flow rate augmented the penetration of the particles. The complex interplay of the finite particle size and the flow inertia decided the spatial deposition of the particles. The findings from this study demonstrate the efficacy of state-of-art simulation frameworks for targeting respiratory disorders.

We present a new multistage method to study the N-Methyl-D-Aspartate (NMDA) neuroreceptor starting from the reconstruction of its crystallographic structure. Thanks to the combination of Homology Modelling, Molecular Dynamics and Lattice Boltzmann simulations, we analyse the allosteric transition of NDMA upon ligand binding and compute the receptor response to ionic passage across the membrane.

In this paper, we deploy the hybrid Lattice Boltzmann - Particle Dynamics (LBPD) method to investigate the transport properties of blood flow within arterioles and venules. The numerical approach is applied to study the transport of Red Blood Cells (RBC) through plasma, highlighting significant agreement with the experimental data in the seminal work by Fahraeus and Lindqvist. Moreover, the results provide evidence of an interesting hand-shaking between the range of validity of the proposed hybrid approach and the domain of viability of particle methods. A joint inspection of accuracy and computational cost, indicate that LBPD offers an appealing multiscale strategy for the study of blood transport across scales of motion, from macroscopic vessels, down to arterioles and venules, where particle methods can eventually take over.

This review discusses the lattice Boltzmann-particle dynamics (LBPD) multiscale paradigm for the simulation of complex states of flowing matter at the interface between physics, chemistry, and biology. In particular, current large-scale LBPD simulations of biopolymer translocation across cellular membranes, molecular transport in ion channels, and amyloid aggregation in cells are described. Prospects are provided for future LBPD explorations in the direction of cellular organization, the direct simulation of full biological organelles, all the way up to physiological scales of potential relevance to future precision-medicine applications, such as the accurate description of homeostatic processes. It is argued that. with the advent of Exascale computing, the mesoscale physics approach advocated in this review may come to age in the next decade and open up new exciting perspectives for physics-based computational medicine.

We investigate the dynamics of a phase-separating binary fluid, containing colloidal dumbbells anchored to the fluid-fluid interface. Extensive lattice Boltzmann-immersed boundary method simulations reveal that the presence of soft dumbbells can significantly affect the curvature dynamics of the interface between phase-separating fluids, even though the coarsening dynamics is left nearly unchanged. In addition, our results show that the curvature dynamics exhibits distinct non-local effects, which might be exploited for the design of new soft mesoscale materials. We point out that the inspection of the statistical dynamics of the curvature can disclose new insights into local inhomogeneities of the binary fluid configuration, as a function of the volume fraction and aspect ratio of the dumbbells.

Computer simulations of bi-continuous two-phase fluids with interspersed dumbbells show that, unlike rigid colloids, soft dumbbells do not lead to arrested coarsening. However, they significantly alter the curvature dynamics of the fluid-fluid interface, whose probability density distributions are shown to exhibit (i) a universal spontaneous transition (observed even in the absence of colloids) from an initial broad-shape distribution towards a highly localized one and (ii) super-diffusive dynamics with long-range effects. Both features may prove useful for the design of novel families of soft porous materials.

We present a lattice Boltzmann model for charged leaky dielectric multiphase fluids in the context of electrified jet simulations, which are of interest for a number of production technologies including electrospinning. The role of nonlinear rheology on the dynamics of electrified jets is considered by exploiting the Carreau model for pseudoplastic fluids. We report exploratory simulations of charged droplets at rest and under a constant electric field, and we provide results for charged jet formation under electrospinning conditions.

The development of high-efficiency porous catalyst membranes critically depends on our understanding of where the majority of the chemical conversions occur within the porous structure. This requires mapping of chemical reactions and mass transport inside the complex nanoscale architecture of porous catalyst membranes which is a multiscale problem in both the temporal and spatial domains. To address this problem, we developed a multiscale mass transport computational framework based on the lattice Boltzmann method that allows us to account for catalytic reactions at the gas-solid interface by introducing a new boundary condition. In good agreement with experiments, the simulations reveal that most catalytic reactions occur near the gas-flow facing side of the catalyst membrane if chemical reactions are fast compared to mass transport within the porous catalyst membrane.

Recent developments of the lattice Boltzmann method for large-scale haemodynamic applications are presented, with special focus on multiscale aspects, including the self-consistent dynamics of suspended biological bodies and their coupling to surface structures, such as the glycocalyx, in the proximity of endothelium using unstructured grids. The description of such multiscale phenomena, each one treated with a suitable variation of the lattice Boltzmann method, opens up new perspectives for a fundamental understanding of the physical mechanisms underlying cardiovascular pathologies, such as plaque growth and the subsequent development of atherosclerotic diseases.

We review recent advances on the mesoscopic modeling of water-like fluids, based on the lattice Boltzmann (LB) methodology. The main idea is to enrich the basic LB (hydro)-dynamics with angular degrees of freedom responding to suitable directional potentials between water-like molecules. The model is shown to reproduce some microscopic features of liquid water, such as an average number of hydrogen bonds per molecules (HBs) between 3 and 4, as well as a qualitatively correct statistics of the hydrogen bond angle as a function of the temperature. Future developments, based on the coupling the present water-like LB model with the dynamics of suspended bodies, such as biopolymers, may open new angles of attack to the simulation of complex biofluidic problems, such as protein folding and aggregation, and the motion of large biomolecules in complex cellular environments.

We present the results of a computational study of coronary trees obtained from CT acquisition at resolution of 0.35mm x 0.35mm x 0.4mm and presenting significant stenotic plaques. We analyze the cardiovascular implications of stenotic plaques for a sizeable number of patients and show that the standard clinical criterion for surgical or percutaneous intervention, based on the Fractional Flow Reserve (FFR), is well reproduced by simulations in a range of inflow conditions that can be finely controlled. The relevance of the present study is related to the reproducibility of FFR data by simulating the coronary trees at global level via high performance simulation methods together with an independent assessment based on in vitro hemodynamics. The data show that controlling the flow Reynolds number is a viable procedure to account for FFR as heart-cycle time averages and maximal hyperemia, as measured in vivo. The reproducibility of the clinical data with simulation offers a systematic approach to measuring the functional implications of stenotic plaques. © 2014 SPIE.

A hydro-kinetic scheme for water-like fluids, based on a lattice version of the Boltzmann equation, is presented and applied to the popular TIP3P model of liquid water. By proceeding in much larger steps than molecular dynamics, the scheme proves to be very effective in attaining global minima of classical pair potentials with directional and radial interactions, as confirmed by further simulations using the three-dimensional Ben-Naim water potential. The scheme is shown to reproduce the propensity of water to form nearly four hydrogen bonds per molecule, as well as their statistical distribution in the presence of thermal fluctuations, at a linear cost of computational time with the system size. This journal is © the Partner Organisations 2014.