The computational route from bilayer membranes to vesicle fusion

Direct and cell-mediated EV-ECM interplay

Extracellular vesicles (EV) are a heterogeneous group of lipid particles excreted by cells. They play an important role in regeneration, development, inflammation, and cancer progression, together with the extracellular matrix (ECM), which they constantly interact with. In this review, we discuss direct and indirect interactions of EVs and the ECM and their impact on different physiological processes. The ECM affects the secretion of EVs, and the properties of the ECM and EVs modulate EVs' diffusion and adhesion. On the other hand, EVs can affect the ECM both directly through enzymes and indirectly through the modulation of the ECM synthesis and remodeling by cells. This review emphasizes recently discovered types of EVs bound to the ECM and isolated by enzymatic digestion, including matrix-bound nanovesicles (MBV) and tissue-derived EV (TiEV). In addition to the experimental studies, computer models of the EV-ECM-cell interactions, from all-atom models to quantitative pharmacology models aiming to improve our understanding of the interaction mechanisms, are also considered. STATEMENT OF SIGNIFICANCE: Application of extracellular vesicles in tissue engineering is an actively developing area. Vesicles not only affect cells themselves but also interact with the matrix and change it. The matrix also influences both cells and vesicles. In this review, different possible types of interactions between vesicles, matrix, and cells are discussed. Furthermore, the united EV-ECM system and its regulation through the cellular activity are presented.

Roughness Volumes: An Improved RoughMob Concept for Predicting the Increase of Molecular Mobility upon Coarse-Graining

The reduced number of degrees of freedom in a coarse-grained molecular model compared to its parent atomistic model not only makes it possible to simulate larger systems for longer time scales but also results in an artificial mobility increase. The RoughMob method [Meinel, M. K. and Müller-Plathe, F. J. Chem. Theory Comput. 2020, 16, 1411.] linked the acceleration factor of the dynamics to the loss of geometric information upon coarse-graining. Our hypothesis is that coarse-graining a multiatom molecule or group into a single spherical bead smooths the molecular surface and, thus, leads to reduced intermolecular friction. A key parameter is the molecular roughness difference, which is calculated via a numerical comparison of the molecular surfaces of both the atomistic and coarse-grained models. Augmenting the RoughMob method, we add the concept of the region where the roughness acts. This information is contained in four so-called roughness volumes. For 17 systems of homogeneous hydrocarbon fluids, simple one-bead coarse-grained models are derived by the structure-based iterative Boltzmann inversion. They include 13 different homogeneous aliphatic and aromatic molecules and two different mapping schemes. We present a simple way to correlate the roughness volumes to the acceleration factor. The resulting relation is able to a priori predict the acceleration factors for an extended size and shape range of hydrocarbon molecules, with different mapping schemes and different densities.

Remodeling of Membrane Shape and Topology by Curvature Elasticity and Membrane Tension

Cellular membranes exhibit a fascinating variety of different morphologies, which are continuously remodeled by transformations of membrane shape and topology. This remodeling is essential for important biological processes (cell division, intracellular vesicle trafficking, endocytosis) and can be elucidated in a systematic and quantitative manner using synthetic membrane systems. Here, recent insights obtained from such synthetic systems are reviewed, integrating experimental observations and molecular dynamics simulations with the theory of membrane elasticity. The study starts from the polymorphism of biomembranes as observed for giant vesicles by optical microscopy and small nanovesicles in simulations. This polymorphism reflects the unusual elasticity of fluid membranes and includes the formation of membrane necks or fluid 'worm holes'. The proliferation of membrane necks generates stable multi-spherical shapes, which can form tubules and tubular junctions. Membrane necks are also essential for the remodeling of membrane topology via membrane fission and fusion. Neck fission can be induced by fine-tuning of membrane curvature, which leads to the controlled division of giant vesicles, and by adhesion-induced membrane tension as observed for small nanovesicles. Challenges for future research include the interplay of curvature elasticity and membrane tension during membrane fusion and the localization of fission and fusion processes within intramembrane domains.

Dissipative particle dynamics simulations in colloid and Interface science: a review

Dissipative particle dynamics (DPD) is one of the most efficient mesoscale coarse-grained methodologies for modeling soft matter systems. Here, we comprehensively review the progress in theoretical formulations, parametrization strategies, and applications of DPD over the last two decades. DPD bridges the gap between the microscopic atomistic and macroscopic continuum length and time scales. Numerous efforts have been performed to improve the computational efficiency and to develop advanced versions and modifications of the original DPD framework. The progress in the parametrization techniques that can reproduce the engineering properties of experimental systems attracted a lot of interest from the industrial community longing to use DPD to characterize, help design and optimize the practical products. While there are still areas for improvements, DPD has been efficiently applied to numerous colloidal and interfacial phenomena involving phase separations, self-assembly, and transport in polymeric, surfactant, nanoparticle, and biomolecules systems.

Adhesive contact between cylindrical (Ebola) and spherical (SARS-CoV-2) viral particles and a cell membrane

A critical event during the process of cell infection by a viral particle is attachment, which is driven by adhesive interactions and resisted by bending and tension. The biophysics of this process has been studied extensively, but the additional role of externally applied force or displacement has generally been neglected. In this work, we study the adhesive force-displacement response of viral particles against a cell membrane. We have built two models: one in which the viral particle is cylindrical (say, representative of a filamentous virus such as Ebola) and another in which it is spherical (such as SARS-CoV-2 and Zika). Our interest is in initial adhesion, in which case deformations are small, and the mathematical model for the system can be simplified considerably. The parameters that characterize the process combine into two dimensionless groups that represent normalized membrane bending stiffness and tension. In the limit where bending dominates, for sufficiently large values of normalized bending stiffness, there is no adhesion between viral particles and the cell membrane without applied force. (The zero external force contact width and pull-off force are both zero.) For large values of normalized membrane tension, the adhesion between virus and cell membrane is weak but stable. (The contact width at zero external force has a small value.) Our results for pull-off force and zero force contact width help to quantify conditions that could aid the development of therapies based on denying the virus entry into the cell by blocking its initial adhesion.

The Application of Coarse-Grained Molecular Dynamics to the Evaluation of Liposome Physical Stability

Physical stability is one of critical characteristics of liposome, especially to its clinical application. Vesicle fusion was one of the common physical stability phenomena that occurred during the long storage period. Because vesicle fusion could be easily checked by the change of vesicle size, it was widely applied in the evaluation of liposome physical stability. However, since the method requires the liposome to be placed under certain conditions for long-term observation, a liposome physical stability test usually takes several weeks, which greatly hinders the research efficiency. In this study, to speed up the research efficiency, coarse-grained molecular dynamics was first applied in the study of liposome physical stability. By analyzing the microprocess of vesicle fusion, two parameters including diffusion constant and the total time of the vesicle morphology transition process were employed to study the liposome physical stability. Then, in order to verify the applicability of two parameters, the physical stability of elastic liposomes and conventional liposomes was compared at 3 different temperatures. It was found that the fusion probability and speed of elastic liposomes were higher than those of conventional liposomes. Thus, elastic liposomes showed a worse physical stability compared with that of conventional liposomes, which was consistent with former research. Through this research, a new efficient method based on coarse-grained molecular dynamics was proposed for the study of liposome physical stability.

Polymersomes: Breaking the Glass Ceiling?

Polymer vesicles, also known as polymersomes, have garnered a lot of interest even before the first report of their fabrication in the mid-1990s. These capsules have found applications in areas such as drug delivery, diagnostics and cellular models, and are made via the self-assembly of amphiphilic block copolymers, predominantly with soft, rubbery hydrophobic segments. Comparatively, and despite their remarkable impermeability, glassy polymersomes (GPs) have been less pervasive due to their rigidity, lack of biodegradability and more restricted fabrication strategies. GPs are now becoming more prominent, thanks to their ability to undergo stable shape-change (e.g., into non-spherical morphologies) as a response to a predetermined trigger (e.g., light, solvent). The basics of block copolymer self-assembly with an emphasis on polymersomes and GPs in particular are reviewed here. The principles and advantages of shape transformation of GPs as well as their general usefulness are also discussed, together with some of the challenges and opportunities currently facing this area.

Reciprocal link between cell biomechanics and exocytosis

A cell is able to sense the biomechanical properties of the environment such as the rigidity of the extracellular matrix and adapt its tension via regulation of plasma membrane and underlying actomyosin meshwork properties. The cell's ability to adapt to the changing biomechanical environment is important for cellular homeostasis and also cell dynamics such as cell growth and motility. Membrane trafficking has emerged as an important mechanism to regulate cell biomechanics. In this review, we summarize the current understanding of the role of cell mechanics in exocytosis, and reciprocally, the role of exocytosis in regulating cell mechanics. We also discuss how cell mechanics and membrane trafficking, particularly exocytosis, can work together to regulate cell polarity and motility.

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.

Rate of hemifusion diaphragm dissipation and ability to form three-junction bound HD determined by lipid composition

Though the hemifusion diaphragm (HD) is widely accepted as an intermediate in bilayer membrane fusion, lipid contributions toward HD stability and dynamics are still not fully understood. In this paper, we study large, binary, protein-free HD systems at varying compositions of negative intrinsic curvature (NIC) lipids using molecular dynamics (MD) simulations of a solvent-free coarse-grained lipid model. Under MD, initially created HDs are found to relax to three major end states depending on the composition and lipid intrinsic curvature. Low compositions of NIC lipids or weak intrinsic curvature result in double-bilayer end states, and moderate compositions of moderate to strong NIC lipids result in metastable fusion pores. Importantly, high compositions of moderate NIC lipids result in a metastable HD that persists beyond μs time scales. NIC lipids stabilize the HD by filling the junction core around the HD. Sorting of NIC lipids toward the three-junction region occurs in fused-endpoint systems, but no significant sorting was seen in systems that end in a double bilayer indicating that high line tension at the triple junction drives HD dissipation faster than sorting can enrich that junction enough to lower that line tension. The appearance of three end states dependent on the NIC lipid composition highlights the necessity of NIC lipids for non-leaky fusion.

Carbon Nanotubes Mediate Fusion of Lipid Vesicles

The fusion of lipid membranes is opposed by high energetic barriers. In living organisms, complex protein machineries carry out this biologically essential process. Here we show that membrane-spanning carbon nanotubes (CNTs) can trigger spontaneous fusion of small lipid vesicles. In coarse-grained molecular dynamics simulations, we find that a CNT bridging between two vesicles locally perturbs their lipid structure. Their outer leaflets merge as the CNT pulls lipids out of the membranes, creating an hourglass-shaped fusion intermediate with still intact inner leaflets. As the CNT moves away from the symmetry axis connecting the vesicle centers, the inner leaflets merge, forming a pore that completes fusion. The distinct mechanism of CNT-mediated membrane fusion may be transferable, providing guidance in the development of fusion agents, e.g., for the targeted delivery of drugs or nucleic acids.

The fusion mechanism of small polymersomes formed by rod-coil diblock copolymers

The fusion mechanism of polymersomes self-assembled by rod-coil copolymers is intrinsically different from that of liposomes due to the effect of chain topology on conformational entropy and molecular packing. The influences of membrane tension, coil-block length, rod-block length, mutual compatibility between the solvent and the rod-coil block, and π-π interaction strength on the fusion pathway are explored by dissipative particle dynamics. The fusion process of spontaneously formed polymersomes generally consists of four stages. In the kissing stage, hopping of rod-blocks forms a connection between two vesicles of a one-legged rod-coil copolymer. In the adhesion stage, a stalk is developed by a few link-up rods and then a stretched diaphragm with rods lying parallel to the stretching direction is formed in the hemi-fusion stage. Eventually, a pore is developed and expanded in the fusion stage. If the membrane tension (τ) is adjusted by deflation/inflation of the polymersomes, the hemi-fusion diaphragm disappears. As τ is reduced, multiple stalks take shape and lead to the formation of inverted micelles, which is the rate-determining step and raises the fusion time substantially. As τ is elevated, a neck is developed after the stalk formation. The fusion process is significantly accelerated. τ of spontaneously formed vesicles varies with the coil-block length, rod-block length, solvent quality, and π-π interaction strength. There exists a critical value of τ below which the fusion process cannot be completed and a hemi-fused polymersome is formed. In addition to τ, the anisotropic steric interactions within the rod layers also resist hopping of longer rod-blocks. The coil layers develop a barrier impeding fusion between vesicles with longer coil-blocks. Consequently, lowering the solvent quality for the coil-block or rod-block facilities the fusion process due to the formation of a thinner coil layer.

Lipid nanotechnology

Nanotechnology is a multidisciplinary field that covers a vast and diverse array of devices and machines derived from engineering, physics, materials science, chemistry and biology. These devices have found applications in biomedical sciences, such as targeted drug delivery, bio-imaging, sensing and diagnosis of pathologies at early stages. In these applications, nano-devices typically interface with the plasma membrane of cells. On the other hand, naturally occurring nanostructures in biology have been a source of inspiration for new nanotechnological designs and hybrid nanostructures made of biological and non-biological, organic and inorganic building blocks. Lipids, with their amphiphilicity, diversity of head and tail chemistry, and antifouling properties that block nonspecific binding to lipid-coated surfaces, provide a powerful toolbox for nanotechnology. This review discusses the progress in the emerging field of lipid nanotechnology.

Dynamics of modeled oligonucleosomes and the role of histone variant proteins in nucleosome organization

Elucidation of the structural dynamics of a nucleosome is of primary importance for understanding the molecular mechanisms that control the nucleosomal positioning. The presence of variant histone proteins in the nucleosome core raises the functional diversity of the nucleosomes in gene regulation and has the profound epigenetic consequences of great importance for understanding the fundamental issues like the assembly of variant nucleosomes, chromatin remodeling, histone posttranslational modifications, etc. Here, we report our observation of the dominant mechanisms of relaxation motions of the oligonucleosomes such as dimer, trimer, and tetramer (in the beads on a string model) with conventional core histones and role of variant histone H2A.Z in the chromatin dynamics using normal mode analysis. Analysis of the directionality of the global dynamics of the oligonucleosome reveals (i) the in-planar stretching as well as out-of-planar bending motions as the relaxation mechanisms of the oligonucleosome and (ii) the freedom of the individual nucleosome in expressing the combination of the above-mentioned motions as the global mode of dynamics. The highly dynamic N-termini of H3 and (H2A.Z-H2B) dimer evidence their participation in the transcriptionally active state. The key role of variant H2A.Z histone as a major source of vibrant motions via weaker intra- and intermolecular correlations is emphasized in this chapter.

Vesicles and vesicle fusion: coarse-grained simulations

Biological cells are highly dynamic, and continually move material around their own volume and between their interior and exterior. Much of this transport encapsulates the material inside phospholipid vesicles that shuttle to and from, fusing with, and budding from, other membranes. A feature of vesicles that is crucial for this transport is their ability to fuse to target membranes and release their contents to the distal side. In industry, some personal care products contain vesicles to help transport reagents across the skin, and research on drug formulation shows that packaging active compounds inside vesicles delays their clearance from the blood stream. In this chapter, we survey the biological role and physicochemical properties of phospholipids, and describe progress in coarse-grained simulations of vesicles and vesicle fusion. Because coarse-grained simulations retain only those molecular details that are thought to influence the large-scale processes of interest, they act as a model embodying our current understanding. Comparing the predictions of these models with experiments reveals the importance of the retained microscopic details and also the deficiencies that can suggest missing details, thereby furthering our understanding of the complex dynamic world of vesicles.

An introduction to dissipative particle dynamics

Dissipative particle dynamics (DPD) is a particle-based mesoscopic simulation method, which facilitates the studies of thermodynamic and dynamic properties of soft matter systems at physically interesting length and time scales. In this method, molecule groups are clustered into the dissipative beads, and this coarse-graining procedure is a very important aspect of DPD as it allows significant computational speed-up. In this chapter, we introduce the DPD methodology, including its theoretical foundation and its parameterization. With this simulation technique, we can study complex behaviors of biological systems, such as the formation of vesicles and their fusion and fission processes, and the phase behavior of lipid membranes.

The Martini coarse-grained force field

The Martini force field is a coarse-grained force field suited for molecular dynamics simulations of biomolecular systems. The force field has been parameterized in a systematic way, based on the reproduction of partitioning free energies between polar and apolar phases of a large number of chemical compounds. In this chapter the methodology underlying the force field is presented together with details of its parameterization and limitations. Then currently available topologies are described with a short overview of the key elements of their parameterization. These include the new polarizable Martini water model. A set of three selected ongoing studies using the Martini force field is presented. Finally the latest lines of development are discussed.

Spatial simulations in systems biology: from molecules to cells

Cells are highly organized objects containing millions of molecules. Each biomolecule has a specific shape in order to interact with others in the complex machinery. Spatial dynamics emerge in this system on length and time scales which can not yet be modeled with full atomic detail. This review gives an overview of methods which can be used to simulate the complete cell at least with molecular detail, especially Brownian dynamics simulations. Such simulations require correct implementation of the diffusion-controlled reaction scheme occurring on this level. Implementations and applications of spatial simulations are presented, and finally it is discussed how the atomic level can be included for instance in multi-scale simulation methods.

Energetics of stalk intermediates in membrane fusion are controlled by lipid composition

We have used X-ray diffraction on the rhombohedral phospholipid phase to reconstruct stalk structures in different pure lipids and lipid mixtures with unprecedented resolution, enabling a quantitative analysis of geometry, as well as curvature and hydration energies. Electron density isosurfaces are used to study shape and curvature properties of the bent lipid monolayers. We observe that the stalk structure is highly universal in different lipid systems. The associated curvatures change in a subtle, but systematic fashion upon changes in lipid composition. In addition, we have studied the hydration interaction prior to the transition from the lamellar to the stalk phase. The results indicate that facilitating dehydration is the key to promote stalk formation, which becomes favorable at an approximately constant interbilayer separation of 9.0 ± 0.5 Å for the investigated lipid compositions.

Spatial modeling of vesicle transport and the cytoskeleton: the challenge of hitting the right road

The membrane trafficking machinery provides a transport and sorting system for many cellular proteins. We propose a mechanistic agent-based computer simulation to integrate and test the hypothesis of vesicle transport embedded into a detailed model cell. The method tracks both the number and location of the vesicles. Thus both the stochastic properties due to the low numbers and the spatial aspects are preserved. The underlying molecular interactions that control the vesicle actions are included in a multi-scale manner based on the model of Heinrich and Rapoport (2005). By adding motor proteins we can improve the recycling process of SNAREs and model cell polarization. Our model also predicts that coat molecules should have a high turnover at the compartment membranes, while the turnover of motor proteins has to be slow. The modular structure of the underlying model keeps it tractable despite the overall complexity of the vesicle system. We apply our model to receptor-mediated endocytosis and show how a polarized cytoskeleton structure leads to polarized distributions in the plasma membrane both of SNAREs and the Ste2p receptor in yeast. In addition, we can couple signal transduction and membrane trafficking steps in one simulation, which enables analyzing the effect of receptor-mediated endocytosis on signaling.

VISUALIZING SOFT MATTER: MESOSCOPIC SIMULATIONS OF MEMBRANES, VESICLES AND NANOPARTICLES

Biological membranes have properties and behavior that emerge from the propagation of the molecular characteristics of their components across many scales. Artificial smart materials, such as drug delivery vehicles and nanoparticles, often rely on modifying naturally-occurring soft matter, such as polymers and lipid vesicles, so that they possess useful behavior. Mesoscopic simulations allow in silico experiments to be easily and cheaply performed on complex, soft materials requiring as input only the molecular structure of the constituents at a coarse-grained level. They can therefore act as a guide to experimenters prior to performing costly assays. Additionally, mesoscopic simulations provide the only currently feasible window on the length and time scales relevant to important biophysical processes such as vesicle fusion. We describe here recent work using Dissipative Particle Dynamics simulations to explore the structure and behavior of amphiphilic membranes, the fusion of vesicles, and the interactions between rigid nanoparticles and soft surfaces.

Membrane protein dynamics from femtoseconds to seconds

Membrane proteins play a key role in energy conversion, transport, signal recognition, transduction, and other fundamental biological processes. Despite considerable progress in experimental techniques, the determination of structure and dynamics of membrane proteins still represents a great challenge. Computer simulation methods are becoming an increasingly important tool not only in the interpretation of experiments but also in the prediction of membrane protein dynamics. In the present review, we give a brief introduction to molecular modeling techniques currently used to explore protein dynamics on time scales ranging from femtoseconds to microseconds. We then describe a few recent example applications of these techniques to membrane proteins. In conclusion, we also discuss some of the newest developments in simulation methodology that have the potential to further extend the time scale accessible to explore (membrane) protein dynamics.

A systematically coarse-grained solvent-free model for quantitative phospholipid bilayer simulations

We present an implicit solvent coarse-grained (CG) model for quantitative simulations of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) bilayers. The absence of explicit solvent enables membrane simulations on large length and time scales at moderate computational expense. Despite improved computational efficiency, the model preserves chemical specificity and quantitative accuracy. The bonded and nonbonded interactions together with the effective cohesion mimicking the hydrophobic effect were systematically tuned by matching structural and mechanical properties from experiments and all-atom bilayer simulations, such as saturated area per lipid, radial distribution functions, density and pressure profiles across the bilayer, P(2) order, etc. The CG lipid model is shown to self-assemble into a bilayer starting from a random dispersion. Its line tension and elastic properties, such as bending and stretching modulus, are semiquantitatively consistent with experiments. The effects of (i) reduced molecular friction and (ii) more efficient integration combine to an overall speed-up of 3-4 orders of magnitude compared to all-atom bilayer simulations. Our CG lipid model is especially useful for studies of large-scale phenomena in membranes that nevertheless require a fair description of chemical specificity, e.g., membrane patches interacting with movable and transformable membrane proteins and peptides.

Protein adsorption and desorption on lipid bilayers

The protein surface usually exhibits one or a few charged spots. If a lipid bilayer contains a significant amount of lipids with oppositely charged head groups, protein adsorption on a bilayer may be energetically favourable due to the protein-lipid electrostatic interaction. The specifics of this case are that the lipids are highly mobile and the protein adsorption is accompanied by the redistribution of lipids between the areas covered and not covered by protein. We present a kinetic model illustrating that this effect is especially interesting if the fraction of the surface covered by charged lipids is relatively low. In this situation, with increasing protein coverage, the protein desorption rate constant rapidly increases while the adsorption rate constant drops, so that there is critical fraction of the area covered by protein. Adsorption above this fraction is hindered both kinetically and thermodynamically.

Relaxation dynamics of fluid membranes

We study the effect of membrane viscosity in the dynamics of liquid membranes-possibly with free or internal boundaries-driven by conservative forces (curvature elasticity and line tension) and dragged by the bulk dissipation of the ambient fluid and the friction occurring when the amphiphilic molecules move relative to each other. To this end, we formulate a continuum model which includes a form of the governing equations for a two-dimensional viscous fluid moving on a curved, time-evolving surface. The effect of membrane viscosity has received very limited attention in previous continuum studies of the dynamics of fluid membranes, although recent coarse-grained discrete simulations suggest its importance. By applying our model to the study of vesiculation and membrane fusion in a simplified geometry, we conclude that membrane viscosity plays a dominant role in the relaxation dynamics of fluid membranes of sizes comparable to those found in eukaryotic cells, and is not negligible in many large synthetic systems of current interest.

Lipids on the move: simulations of membrane pores, domains, stalks and curves

In this review we describe the state-of-the-art of computer simulation studies of lipid membranes. We focus on collective lipid-lipid and lipid-protein interactions that trigger deformations of the natural lamellar membrane state, showing that many important biological processes including self-aggregation of membrane components into domains, the formation of non-lamellar phases, and membrane poration and curving, are now amenable to detailed simulation studies.

Insight or illusion? Seeing inside the cell with mesoscopic simulations

The expulsion of material from a cell by fusion of vesicles at the plasma membrane, and the entry of a virus by membrane invagination are complex membrane-associated processes whose control is crucial to cell survival. Our ability to visualize the dynamics of such processes experimentally is limited by spatial resolution and the speed of molecular rearrangements. The increase in computing power of the last few decades enables the construction of computational tools for observing cellular processes in silico. As experiments yield increasing amounts of data on the protein and lipid constituents of the cell, computer simulations parametrized using this data are beginning to allow models of cellular processes to be interrogated in ways unavailable in the laboratory. Mesoscopic simulations retain only those molecular features that are believed to be relevant to the processes of interest. This allows the dynamics of spatially heterogeneous membranes and the crowded cytoplasmic environment to be followed at a modest computational cost. The price for such power is that the atomic detail of the constituents is much lower than in atomistic Molecular Dynamics simulations. We argue that this price is worth paying because mesoscopic simulations can generate new insight into the complex, dynamic life of a cell.

Lipid Exchange during Contact between Oppositely Charged Lipid Bilayers

The mechanistic details of processes occurring in and during contact between lipid bilayers are still poorly understood due to their complexity on the mesoscopic scale. Here, we analyze lipid exchange during contact of oppositely charged lipid bilayers. Specifically, we explore a generic mechanism, where this process occurs via diffusion of individual lipids between the layers. Our estimates indicate that this scenario is feasible on the time scale of conventional experiments and also on the time scale of biochemical processes in cells.

The MARTINI Force Field: Coarse Grained Model for Biomolecular Simulations

We present an improved and extended version of our coarse grained lipid model. The new version, coined the MARTINI force field, is parametrized in a systematic way, based on the reproduction of partitioning free energies between polar and apolar phases of a large number of chemical compounds. To reproduce the free energies of these chemical building blocks, the number of possible interaction levels of the coarse-grained sites has increased compared to those of the previous model. Application of the new model to lipid bilayers shows an improved behavior in terms of the stress profile across the bilayer and the tendency to form pores. An extension of the force field now also allows the simulation of planar (ring) compounds, including sterols. Application to a bilayer/cholesterol system at various concentrations shows the typical cholesterol condensation effect similar to that observed in all atom representations.

Reptational dynamics in dissipative particle dynamics simulations of polymer melts

Understanding the fundamental properties of polymeric liquids remains a challenge in materials science and soft matter physics. Here, we present a simple and computationally efficient criterion for topological constraints, i.e., uncrossability of chains, in polymeric liquids using the dissipative particle dynamics (DPD) method. No new length scales or forces are added. To demonstrate that this approach really prevents chain crossings, we study a melt of linear homopolymers. We show that for short chains the model correctly reproduces Rouse-like dynamics whereas for longer chains the dynamics becomes reptational as the chain length is increased--something that is not attainable using standard DPD or other coarse-grained soft potential methods.

Comparison of Thermodynamic Properties of Coarse-Grained and Atomic-Level Simulation Models

Thermodynamic data are often used to calibrate or test amomic-level (AL) force fields for molecular dynamics (MD) simulations. In contrast, the majority of coarse-grained (CG) force fields do not rely extensively on thermodynamic quantities. Recently, a CG force field for lipids, hydrocarbons, ions, and water, in which approximately four non-hydrogen atoms are mapped onto one interaction site, has been proposed and applied to study various aspects of lipid systems. To date, no extensive investigation of its capability to describe salvation thermodynamics has been undertaken. In the present study, a detailed picture of vaporization, solvation, and phase-partitioning thermodynamics for liquid hydrocarbons and water was obtained at CG and AL resolutions, in order to compare the two types or models and evaluate their ability to describe thermodynamic properties in the temperature range between 263 and 343 K. Both CG and AL models capture the experimental dependence of the thermodynamic properties on the temperature, albeit a systematically weaker dependence is found for the CG model. Moreover, deviations are found for solvation thermodynamics and for the corresponding enthalpy-entropy compensation for the CG model. Particularly water/oil repulsion seems to be overestimated. However, the results suggest that the thermodynamic properties considered should be reproducible by a CG model provided it is reparametrized on the basis of these liquid-phase properties.