Coarse-grain model for lipid bilayer self-assembly and dynamics: Multiparticle collision description of the solvent

Dynamic framework for large-scale modeling of membranes and peripheral proteins

In this chapter, we present a novel computational framework to study the dynamic behavior of extensive membrane systems, potentially in interaction with peripheral proteins, as an alternative to conventional simulation methods. The framework effectively describes the complex dynamics in protein-membrane systems in a mesoscopic particle-based setup. Furthermore, leveraging the hydrodynamic coupling between the membrane and its surrounding solvent, the coarse-grained model grounds its dynamics in macroscopic kinetic properties such as viscosity and diffusion coefficients, marrying the advantages of continuum- and particle-based approaches. We introduce the theoretical background and the parameter-space optimization method in a step-by-step fashion, present the hydrodynamic coupling method in detail, and demonstrate the application of the model at each stage through illuminating examples. We believe this modeling framework to hold great potential for simulating membrane and protein systems at biological spatiotemporal scales, and offer substantial flexibility for further development and parametrization.

Scale-invariance in miniature coarse-grained red blood cells by fluctuation analysis

To accurately represent the morphological and elastic properties of a human red blood cell, Fu et al. [Fu et al., Lennard-Jones type pair-potential method for coarse-grained lipid bilayer membrane simulations in LAMMPS, 2017, 210, 193-203] recently developed a coarse-grained molecular dynamics model with particular detail in the membrane. However, such a model accrues an extremely high computational cost for whole-cell simulation when assuming an appropriate length scaling - that of the bilayer thickness. To date, the model has only simulated "miniature" cells in order to circumvent this, with the a priori assumption that these miniaturised cells correctly represent their full-sized counterparts. The present work assesses the validity of this approach, by testing the scale invariance of the model through simulating cells of various diameters; first qualitatively in their shape evolution, then quantitatively by measuring their bending rigidity through fluctuation analysis. Cells of diameter of at least 0.5 μm were able to form the characteristic biconcave shape of human red blood cells, though smaller cells instead equilibrated to bowl-shaped stomatocytes. Thermal fluctuation analysis showed the bending rigidity to be constant over all cell sizes tested, and consistent between measurements on the whole-cell and on a planar section of bilayer. This is as expected from the theory on both counts. Therefore, we confirm that the evaluated model is a good representation of a full-size RBC when the model diameter is ≥0.5 μm, in terms of the morphological and mechanical properties investigated.

Hydrodynamic coupling for particle-based solvent-free membrane models

The great challenge with biological membrane systems is the wide range of scales involved, from nanometers and picoseconds for individual lipids to the micrometers and beyond millisecond for cellular signaling processes. While solvent-free coarse-grained membrane models are convenient for large-scale simulations and promising to provide insight into slow processes involving membranes, these models usually have unrealistic kinetics. One major obstacle is the lack of an equally convenient way of introducing hydrodynamic coupling without significantly increasing the computational cost of the model. To address this, we introduce a framework based on anisotropic Langevin dynamics, for which major in-plane and out-of-plane hydrodynamic effects are modeled via friction and diffusion tensors from analytical or semi-analytical solutions to Stokes hydrodynamic equations. Using this framework, in conjunction with our recently developed membrane model, we obtain accurate dispersion relations for planar membrane patches, both free-standing and in the vicinity of a wall. We briefly discuss how non-equilibrium dynamics is affected by hydrodynamic interactions. We also measure the surface viscosity of the model membrane and discuss the affecting dissipative mechanisms.

Asymmetric Lipid Membranes under Shear Flows: A Dissipative Particle Dynamics Study

We investigate the phase behavior of the asymmetric lipid membranes under shear flows, using the dissipative particle dynamics simulation. Two cases, the weak and strong shear flows, are considered for the asymmetric lipid microstructures. Three typical asymmetric structures, the membranes, tubes, and vesicle, are included in the phase diagrams, where the effect of two different types of lipid chain length on the formation of asymmetric membranes is evaluated. The dynamic processes are demonstrated for the asymmetric membranes by calculating the average radius of gyration and shape factor. The result indicates that different shear flows will affect the shape of the second type of lipid molecules; the shape of the first type of lipid molecules is more stable than that of the second type of lipid molecules. The mechanical properties are investigated for the asymmetric membranes by analyzing the interface tension. The results reveal an absolute pressure at the junctions of different types of particles under the weak shear flow; the other positions are almost in a state of no pressure; there is almost no pressure inside the asymmetric lipid membrane structure under the strong shear flow. The findings will help us to understand the potential applications of asymmetric lipid microstructures in the biological and medical fields.

Quantification and demonstration of the collective constriction-by-ratchet mechanism in the dynamin molecular motor

Dynamin is a protein that is a central player in endocytosis, a process that mediates the entry of diverse particles into cells, from nutrients to viruses. Dynamin’s primary activity is to use guanosine triphosphate as fuel to constrict and cut membrane tubes. Key quantitative aspects of its function remain yet unclear. In this work, we determine the strength of an individual dynamin motor. Then, by building a detailed computational model resolving individual motors, we demonstrate that dynamin produces sufficient force to tightly constrict a membrane tube when most of its motors are simultaneously cooperating. Hence, we quantitatively validate the prevailing constriction-by-ratchet model for nature’s strongest torque-generating motor: the dynamin “nanomuscle.”

Dynamin oligomerizes into helical filaments on tubular membrane templates and, through constriction, cleaves them in a GTPase-driven way. Structural observations of GTP-dependent cross-bridges between neighboring filament turns have led to the suggestion that dynamin operates as a molecular ratchet motor. However, the proof of such mechanism remains absent. Particularly, it is not known whether a powerful enough stroke is produced and how the motor modules would cooperate in the constriction process. Here, we characterized the dynamin motor modules by single-molecule Förster resonance energy transfer (smFRET) and found strong nucleotide-dependent conformational preferences. Integrating smFRET with molecular dynamics simulations allowed us to estimate the forces generated in a power stroke. Subsequently, the quantitative force data and the measured kinetics of the GTPase cycle were incorporated into a model including both a dynamin filament, with explicit motor cross-bridges, and a realistic deformable membrane template. In our simulations, collective constriction of the membrane by dynamin motor modules, based on the ratchet mechanism, is directly reproduced and analyzed. Functional parallels between the dynamin system and actomyosin in the muscle are seen. Through concerted action of the motors, tight membrane constriction to the hemifission radius can be reached. Our experimental and computational study provides an example of how collective motor action in megadalton molecular assemblies can be approached and explicitly resolved.

Interaction between Bottlebrush Polymers and Phospholipid Membranes in Solutions

In this work, the interactions between bottlebrush polymers and phospholipid membranes were investigated using dissipative particle dynamics simulations. The weak and strong adsorption phenomena between the polymers and membranes were examined by calculating the system parameters. A spring model was introduced to explain the variances in the shape factors and the radius of gyration of the bottlebrush polymers, as well as the order parameters of the phospholipid membrane in the pulling processes. This work provides further understanding for the application of bottlebrush polymers in biological processes.

Mechanochemical enzymes and protein machines as hydrodynamic force dipoles: the active dimer model

Mechanochemically active enzymes change their shapes within every turnover cycle. Therefore, they induce circulating flows in the solvent around them and behave as oscillating hydrodynamic force dipoles. Because of non-equilibrium fluctuating flows collectively generated by the enzymes, mixing in the solution and diffusion of passive particles within it are expected to get enhanced. Here, we investigate the intensity and statistical properties of such force dipoles in the minimal active dimer model of a mechanochemical enzyme. In the framework of this model, novel estimates for hydrodynamic collective effects in solution and in lipid bilayers under rapid rotational diffusion are derived, and available experimental and computational data is examined.

Predicting asymmetric phospholipid microstructures in solutions

Asymmetric phospholipid microstructures, such as asymmetric phospholipid membranes, have potential applications in biological and medicinal processes. Here, we used the dissipative particle dynamics simulation method to predict the asymmetric phospholipid microstructures in aqueous solutions. The asymmetric phospholipid membranes, tubes and vesicles are determined and characterized by the chain density distributions and order parameters. The phase diagrams are constructed to evaluate the effects of the chain length on the asymmetric structure formations at equilibrium states, while the average radius of gyration and shape factors are calculated to analyze the asymmetric structure formations in the non-equilibrium processes. Meanwhile, we predicted the mechanical properties of the asymmetric membranes by analyzing the spatial distributions of the interface tensions and osmotic pressures in solutions.

Large-scale simulation of biomembranes incorporating realistic kinetics into coarse-grained models

Biomembranes are two-dimensional assemblies of phospholipids that are only a few nanometres thick, but form micrometre-sized structures vital to cellular function. Explicit molecular modelling of biologically relevant membrane systems is computationally expensive due to the large number of solvent particles and slow membrane kinetics. Coarse-grained solvent-free membrane models offer efficient sampling but sacrifice realistic kinetics, thereby limiting the ability to predict pathways and mechanisms of membrane processes. Here, we present a framework for integrating coarse-grained membrane models with continuum-based hydrodynamics. This framework facilitates efficient simulation of large biomembrane systems with large timesteps, while achieving realistic equilibrium and non-equilibrium kinetics. It helps to bridge between the nanometer/nanosecond spatiotemporal resolutions of coarse-grained models and biologically relevant time- and lengthscales. As a demonstration, we investigate fluctuations of red blood cells, with varying cytoplasmic viscosities, in 150-milliseconds-long trajectories, and compare kinetic properties against single-cell experimental observations.

Particle-based membrane model for mesoscopic simulation of cellular dynamics

We present a simple and computationally efficient coarse-grained and solvent-free model for simulating lipid bilayer membranes. In order to be used in concert with particle-based reaction-diffusion simulations, the model is purely based on interacting and reacting particles, each representing a coarse patch of a lipid monolayer. Particle interactions include nearest-neighbor bond-stretching and angle-bending and are parameterized so as to reproduce the local membrane mechanics given by the Helfrich energy density over a range of relevant curvatures. In-plane fluidity is implemented with Monte Carlo bond-flipping moves. The physical accuracy of the model is verified by five tests: (i) Power spectrum analysis of equilibrium thermal undulations is used to verify that the particle-based representation correctly captures the dynamics predicted by the continuum model of fluid membranes. (ii) It is verified that the input bending stiffness, against which the potential parameters are optimized, is accurately recovered. (iii) Isothermal area compressibility modulus of the membrane is calculated and is shown to be tunable to reproduce available values for different lipid bilayers, independent of the bending rigidity. (iv) Simulation of two-dimensional shear flow under a gravity force is employed to measure the effective in-plane viscosity of the membrane model and show the possibility of modeling membranes with specified viscosities. (v) Interaction of the bilayer membrane with a spherical nanoparticle is modeled as a test case for large membrane deformations and budding involved in cellular processes such as endocytosis. The results are shown to coincide well with the predicted behavior of continuum models, and the membrane model successfully mimics the expected budding behavior. We expect our model to be of high practical usability for ultra coarse-grained molecular dynamics or particle-based reaction-diffusion simulations of biological systems.

Self-assembly and rheology of dipolar colloids in simple shear studied using multi-particle collision dynamics

Magnetic nanoparticles in a colloidal solution self-assemble in various aligned structures, which has a profound influence on the flow behavior. However, the precise role of the microstructure in the development of the rheological response has not been reliably quantified. We investigate the self-assembly of dipolar colloids in simple shear using hybrid molecular dynamics and multi-particle collision dynamics simulations with explicit coarse-grained hydrodynamics, conduct simulated rheometric studies and apply micromechanical models to produce master curves, showing evidence of the universality of the structural behavior governed by the competition between the bonding (dipolar) and erosive (thermal and/or hydrodynamic) stresses. The simulations display viscosity changes across several orders of magnitude in fair quantitative agreement with various literature sources, substantiating the universality of the approach, which seems to apply generally across vastly different length scales and a broad range of physical systems.

Toward Hydrodynamics with Solvent Free Lipid Models: STRD Martini

Solvent hydrodynamics are incorporated into simulations of the solvent-free Dry Martini model. The solvent hydrodynamics are modeled with the stochastic rotation dynamics (SRD) algorithm, a particle-based method for resolving fluid hydrodynamics. SRD does not require calculation of particle-particle distances in the solvent, and so is scalable to arbitrary volumes of solvent with minimal additional computational overhead. The viscosity of the solvent is easily tuned via parameters of the algorithm to span an order of magnitude in viscosity around the viscosity of water at room temperature. The combination "Stochastic Thermostatted Rotation Dynamics (STRD) with Martini" was implemented in Gromacs v.5.01. Simulations of an SRD/palmitoyloleoylphosphatidylcholine membrane demonstrate that the solvent may be included without reparametrizing the lipid model, with minimal perturbation to the thermodynamics. A recent generalization of Saffman-Delbruck theory to periodic geometries by Camley and Brown indicates that lipid dynamics are contaminated by a finite-size effect in typical molecular dynamics (MD) simulations, and that very large systems are required for quantitative simulation of dynamics. Analysis of lipid translational diffusion in this work shows good agreement with the theory, and with explicitly solvated simulations. This indicates that STRD Martini is a viable approach for quantitative simulation of membrane dynamics and does not require massive computational overhead to model the solvent.

Hydrodynamic collective effects of active proteins in biological membranes

Lipid bilayers forming biological membranes are known to behave as viscous two-dimensional fluids on submicrometer scales; usually they contain a large number of active protein inclusions. Recently, it was shown [A. S. Mikhailov and R. Kapral, Proc. Natl. Acad. Sci. USA 112, E3639 (2015)PNASA60027-842410.1073/pnas.1506825112] that such active proteins should induce nonthermal fluctuating lipid flows leading to diffusion enhancement and chemotaxislike drift for passive inclusions in biomembranes. Here, a detailed analytical and numerical investigation of such effects is performed. The attention is focused on the situations when proteins are concentrated within lipid rafts. We demonstrate that passive particles tend to become attracted by active rafts and are accumulated inside them.

On the dynamics of molecular self-assembly and the structural analysis of bilayer membranes using coarse-grained molecular dynamics simulations

We present a molecular dynamics simulation study of the self-assembly of coarse-grained lipid molecules from unbiased random initial configurations. Our lipid model is based on a well-tried CG polymer model with an additional potential that mimics the hydrophobic properties of lipid tails. We find that several stages of self-organization of lipid clusters are involved in the dynamics of bilayer formation and that the resulting equilibrium structures sensitively depend on the strength of hydrophobic interactions hc of the lipid tails and on temperature T. The obtained stable lipid membranes are quantitatively analyzed with respect to their local structure and their degree of order. At equilibrium, we obtain self-stabilizing bilayer membrane structures that exhibit a bending stiffness κB and compression modulus KC comparable to experimental measurements under physiological conditions. We present a phase diagram of our lipid model which covers a sol-gel transition, a liquid (or gel-like) phase including stable bilayer structures and vesicle formation, as well as a quasi-crystalline phase. We also determine the exact conditions for temperature T and degree of hydrophobicity hc for stable bilayer formation including closed vesicles.

Modeling Nanoparticle Targeting to a Vascular Surface in Shear Flow Through Diffusive Particle Dynamics

Nanoparticles are regarded as promising carriers for targeted drug delivery and imaging probes. A fundamental understanding of the dynamics of polymeric nanoparticle targeting to receptor-coated vascular surfaces is therefore of great importance to enhance the design of nanoparticles toward improving binding ability. Although the effects of particle size and shear flow on the binding of nanoparticles to a vessel wall have been studied at the particulate level, a computational model to investigate the details of the binding process at the molecular level has not been developed. In this research, dissipative particle dynamics simulations are used to study nanoparticles with diameters of several nanometers binding to receptors on vascular surfaces under shear flow. Interestingly, shear flow velocities ranging from 0 to 2000 s(-1) had no effect on the attachment process of nanoparticles very close to the capillary wall. Increased binding energy between the ligands and wall caused a corresponding linear increase in bonding ability. Our simulations also indicated that larger nanoparticles and those of rod shape with a higher aspect ratio have better binding ability than those of smaller size or rounder shape.

Hydrodynamic collective effects of active protein machines in solution and lipid bilayers

The cytoplasm and biomembranes in biological cells contain large numbers of proteins that cyclically change their shapes. They are molecular machines that can function as molecular motors or carry out various other tasks in the cell. Many enzymes also undergo conformational changes within their turnover cycles. We analyze the advection effects that nonthermal fluctuating hydrodynamic flows induced by active proteins have on other passive molecules in solution or membranes. We show that the diffusion constants of passive particles are enhanced substantially. Furthermore, when gradients of active proteins are present, a chemotaxis-like drift of passive particles takes place. In lipid bilayers, the effects are strongly nonlocal, so that active inclusions in the entire membrane contribute to local diffusion enhancement and the drift. All active proteins in a biological cell or in a membrane contribute to such effects and all passive particles, and the proteins themselves, will be subject to them.

Simulating the complex cell design of Trypanosoma brucei and its motility

The flagellate Trypanosoma brucei, which causes the sleeping sickness when infecting a mammalian host, goes through an intricate life cycle. It has a rather complex propulsion mechanism and swims in diverse microenvironments. These continuously exert selective pressure, to which the trypanosome adjusts with its architecture and behavior. As a result, the trypanosome assumes a diversity of complex morphotypes during its life cycle. However, although cell biology has detailed form and function of most of them, experimental data on the dynamic behavior and development of most morphotypes is lacking. Here we show that simulation science can predict intermediate cell designs by conducting specific and controlled modifications of an accurate, nature-inspired cell model, which we developed using information from live cell analyses. The cell models account for several important characteristics of the real trypanosomal morphotypes, such as the geometry and elastic properties of the cell body, and their swimming mechanism using an eukaryotic flagellum. We introduce an elastic network model for the cell body, including bending rigidity and simulate swimming in a fluid environment, using the mesoscale simulation technique called multi-particle collision dynamics. The in silico trypanosome of the bloodstream form displays the characteristic in vivo rotational and translational motility pattern that is crucial for survival and virulence in the vertebrate host. Moreover, our model accurately simulates the trypanosome's tumbling and backward motion. We show that the distinctive course of the attached flagellum around the cell body is one important aspect to produce the observed swimming behavior in a viscous fluid, and also required to reach the maximal swimming velocity. Changing details of the flagellar attachment generates less efficient swimmers. We also simulate different morphotypes that occur during the parasite's development in the tsetse fly, and predict a flagellar course we have not been able to measure in experiments so far.

Typanosoma brucei is a uni-cellular parasite that causes the sleeping sickness, a deadly disease for humans that also occurs in livestock. Injected into the mammalian host by the tsetse fly, the trypanosome travels through the blood stream, where it proliferates, and ultimately can be taken up again by a fly during a bloodmeal. In the tsetse fly, it continues its development with several morphological changes to the cell body plan. During its life cycle, the trypanosome meets different microenvironments, such as the mammalian's bloodstream and the tsetse fly's midgut, proventriculus, foregut, and salivary gland. The cell body of the trypanosome has the shape of a spindle along which an eukaryotic flagellum is attached. We have developed an accurate, in silico model trypanosome using information from live cell analyses. Performing computer simulations, we are able to reproduce all motility patterns of the blood-stream form in typical cell culture medium. Modifying the cell design, we show that the helical course of the flagellar attachment optimizes the trypanosome's swimming speed. We also design trypanosomal morphotypes that occur in the tsetse fly. Simulation science thereby provides an investigative tool to systematically explore the morphologcial diversity during the trypanosome's life cycle even beyond experimental capabilities.

Lattice solution model for order-disorder transitions in membranes and Langmuir monolayers

Lipid monolayers and bilayers have been used as experimental models for the investigation of membrane thermal transitions. The main transition takes place near ambient temperatures for several lipids and reflects the order-disorder transition of lipid hydrocarbonic chains, which is accompanied by a surface density gap. Equivalence between the transitions in the two systems has been argued by several authors. The two-state statistical model adopted by numerous authors for different properties of the membrane, such as permeability, diffusion, and mixture or insertion of cholesterol or protein, is inadequate for the description of charged membranes, since it lacks a proper description of surface density. We propose a lattice solution model which adds interactions with water molecules to lipid-lipid interactions and obtain its thermal properties under a mean-field approach. Density variations, although concomitant with chain order variations, are independent of the latter. The model presents both chain order and gas-liquid transitions, and extends the range of applicability of previous models, yielding Langmuir isotherms in the full range of pressures and areas.

Destruction of cancer cells by laser-induced shock waves: recent developments in experimental treatments and multiscale computer simulations

In this emerging area article we review recent progress in the mechanical destruction of cancer cells using laser-induced shock waves. The pure mechanical damaging and destruction of cancer cells associated with this technique possibly opens up a new route to tumor treatments and the corresponding therapies. At the same time progress in multiscale simulation techniques makes it possible to simulate mechanical properties of soft biological matter such as membranes, cytoskeletal networks and even whole cells and tissue. In this way an interdisciplinary approach to understanding key mechanisms in shock wave interactions with biological matter has become accessible. Mechanical properties of biological materials are also critical for many physiological processes and cover length scales ranging from the atomistic to the macroscopic scale. We argue that the latest developments and progress in experimentation enable the investigation of the shock wave interaction with cancer cells on multiple time- and length-scales. In this way the integrated use of experiment and simulation can address key challenges in this field. The exploration of the biological effects of laser-generated shock waves on a fundamental level constitutes an emerging multidisciplinary research area combining scientific methods from the areas of physics, biology, medicine and computer science.

Coarse-grain simulations of active molecular machines in lipid bilayers

A coarse-grain method for simulations of the dynamics of active protein inclusions in lipid bilayers is described. It combines the previously proposed hybrid simulations of bilayers [M.-J. Huang, R. Kapral, A. S. Mikhailov, and H.-Y. Chen, J. Chem. Phys. 137, 055101 (2012)], based on molecular dynamics for the lipids and multi-particle collision dynamics for the solvent, with an elastic-network description of active proteins. The method is implemented for a model molecular machine which performs active conformational motions induced by ligand binding and its release after reaction. The situation characteristic for peripheral membrane proteins is considered. Statistical investigations of the effects of single active or passive inclusions on the shape of the membrane are carried out. The results show that the peripheral machine produces asymmetric perturbations of the thickness of two leaflets of the membrane. It also produces a local saddle in the midplane height of the bilayer. Analysis of the power spectrum of the fluctuations of the membrane midplane shows that the conformational motion of the machine perturbs these membrane fluctuations. The hydrodynamic lipid flows induced by cyclic conformational changes in the machine are analyzed. It is shown that such flows are long-ranged and should provide an additional important mechanism for interactions between active inclusions in biological membranes.