Zwitterionic lipid assemblies: molecular dynamics studies of monolayers, bilayers, and vesicles using a new coarse grain force field

Molecular-level insight into the binding of arginine to a zwitterionic Langmuir monolayer

Arginine molecules bind to a DPPC monolayer, altering the interfacial electrostatic potential and the lateral mobility of the lipids, while having little effect on the compression isotherm of the monolayer.

The structure and dynamics of Nano Particles encapsulated by the SDS monolayer collapse at the water/TCE interface

The super-saturated surfactant monolayer collapses with the nanoparticles (NPs) at the water/trichloroethylene (TCE) interface are investigated using molecular dynamics (MD) simulations. The results show that sodium alkyl sulfate (SDS) monolayer collapse is initiated by buckling and followed primarily by budding and the bud encapsulating the NPs and oil molecules. The developed bud detaches from the monolayer into a water phase and forms the swollen micelle emulsion with NPs and oil molecules. We investigate the wavelength of the initial budding and the theoretical description of the budding process. The wavelength of the monolayer increases with bending modulus. The energy barrier of the budding can be easily overcome by thermal fluctuation energy, which indicates that budding process proceeds rapidly.

The cellular membrane as a mediator for small molecule interaction with membrane proteins

The cellular membrane constitutes the first element that encounters a wide variety of molecular species to which a cell might be exposed. Hosting a large number of structurally and functionally diverse proteins associated with this key metabolic compartment, the membrane not only directly controls the traffic of various molecules in and out of the cell, it also participates in such diverse and important processes as signal transduction and chemical processing of incoming molecular species. In this article, we present a number of cases where details of interaction of small molecular species such as drugs with the membrane, which are often experimentally inaccessible, have been studied using advanced molecular simulation techniques. We have selected systems in which partitioning of the small molecule with the membrane constitutes a key step for its final biological function, often binding to and interacting with a protein associated with the membrane. These examples demonstrate that membrane partitioning is not only important for the overall distribution of drugs and other small molecules into different compartments of the body, it may also play a key role in determining the efficiency and the mode of interaction of the drug with its target protein. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

Molecular dynamics simulations of cholesterol-rich membranes using a coarse-grained force field for cyclic alkanes

The architecture of a biological membrane hinges upon the fundamental fact that its properties are determined by more than the sum of its individual components. Studies on model membranes have shown the need to characterize in molecular detail how properties such as thickness, fluidity, and macroscopic bending rigidity are regulated by the interactions between individual molecules in a non-trivial fashion. Simulation-based approaches are invaluable to this purpose but are typically limited to short sampling times and model systems that are often smaller than the required properties. To alleviate both limitations, the use of coarse-grained (CG) models is nowadays an established computational strategy. We here present a new CG force field for cholesterol, which was developed by using measured properties of small molecules, and can be used in combination with our previously developed force field for phospholipids. The new model performs with precision comparable to atomistic force fields in predicting the properties of cholesterol-rich phospholipid bilayers, including area per lipid, bilayer thickness, tail order parameter, increase in bending rigidity, and propensity to form liquid-ordered domains in ternary mixtures. We suggest the use of this model to quantify the impact of cholesterol on macroscopic properties and on microscopic phenomena involving localization and trafficking of lipids and proteins on cellular membranes.

Coarse-grained molecular dynamics study of membrane fusion: Curvature effects on free energy barriers along the stalk mechanism

The effects of membrane curvature on the free energy barrier for membrane fusion have been investigated using coarse-grained molecular dynamics (CG-MD) simulations, assuming that fusion takes place through a stalk intermediate. Free energy barriers were estimated for stalk formation as well as for fusion pore formation using the guiding potential method. Specifically, the three different geometries of two apposed membranes were considered: vesicle-vesicle, vesicle-planar, and planar-planar membranes. The free energy barriers for the resulting fusion were found to depend importantly on the fusing membrane geometries; the lowest barrier was obtained for vesicular membranes. Further, lipid sorting was observed in fusion of the mixed membranes of dimyristoyl phosphatidylcholine and dioleoyl phosphatidylethanolamine (DOPE). Specifically, DOPE molecules were found to assemble around the stalk to support the highly negative curved membrane surface. A consistent result for lipid sorting was observed when a simple continuum model (CM) was used, where the Helfrich energy and mixing entropy of the lipids were taken into account. However, the CM predicts a much higher free energy barrier than found using CG-MD. This discrepancy originates from the conformational changes of lipids, which were not considered in the CM. The results of the CG-MD simulations reveal that a large conformational change in the lipid takes place around the stalk region, which results in a reduction of free energy barriers along the stalk mechanism of membrane fusion.

Multiscale approach for the construction of equilibrated all-atom models of a poly(ethylene glycol)-based hydrogel

A multiscale modeling approach is presented for the efficient construction of an equilibrated all-atom model of a cross-linked poly(ethylene glycol) (PEG)-based hydrogel using the all-atom polymer consistent force field (PCFF). The final equilibrated all-atom model was built with a systematic simulation toolset consisting of three consecutive parts: (1) building a global cross-linked PEG-chain network at experimentally determined cross-link density using an on-lattice Monte Carlo method based on the bond fluctuation model, (2) recovering the local molecular structure of the network by transitioning from the lattice model to an off-lattice coarse-grained (CG) model parameterized from PCFF, followed by equilibration using high performance molecular dynamics methods, and (3) recovering the atomistic structure of the network by reverse mapping from the equilibrated CG structure, hydrating the structure with explicitly represented water, followed by final equilibration using PCFF parameterization. The developed three-stage modeling approach has application to a wide range of other complex macromolecular hydrogel systems, including the integration of peptide, protein, and/or drug molecules as side-chains within the hydrogel network for the incorporation of bioactivity for tissue engineering, regenerative medicine, and drug delivery applications.

Collapse of a lipid-coated nanobubble and subsequent liposome formation

We investigate the collapse of a lipid-coated nanobubble and subsequent formation of a lipid vesicle by coarse grained molecular dynamics simulations. A spherical nanobubble coated with a phospholipid monolayer in water is a model of an aqueous dispersion of phospholipids under negative pressure during sonication. When subjected to a positive pressure, the bubble shape deforms into an irregular spherical shape and the monolayer starts to buckle and fold locally. The local folds grow rapidly in multiple directions and forming a discoidal membrane with folds of various amplitudes. Folds of small amplitude disappear in due course and the membrane develops into a unilamellar vesicle via a bowl shape. Folds with large amplitude develop into a bowl shape and a multivesicular shape forms. The membrane shape due to bubble collapse can be an important factor governing the vesicular shape during sonication.

Molecular Mechanism of Ionic-Liquid-Induced Membrane Disruption: Morphological Changes to Bilayers, Multilayers, and Vesicles

The application of ionic liquids (ILs) in many industrially relevant processes provides an urgent need to better understand their molecular interactions with biological systems. A detailed understanding of the cytotoxicity mechanism of ILs can be helpful in facilitating the molecular design of nontoxic ILs. Using coarse-grained molecular dynamics (MD) simulations, we investigate the effects of imidazolium-based ILs on several lipid bilayer morphologies. Our results demonstrate that the asymmetric insertion of IL cations into one side of a lipid bilayer leaflet enhances the leaflet strain, which upon reaching a critical value triggers a morphological disruption in the bilayer. Consistently, the bending modulus of the bilayer is reduced by 1 to 2 orders of magnitude relative to that of an IL-free planar bilayer prior to the disruption event. Our results suggest that ILs that can easily insert into the lipid bilayer without diffusing across or inducing lipid flip-flop can be more disruptive to a lipid biomembrane.

Permeability across lipid membranes

Molecular permeation through lipid membranes is a fundamental biological process that is important for small neutral molecules and drug molecules. Precise characterization of free energy surface and diffusion coefficients along the permeation pathway is required in order to predict molecular permeability and elucidate the molecular mechanisms of permeation. Several recent technical developments, including improved molecular models and efficient sampling schemes, are illustrated in this review. For larger penetrants, explicit consideration of multiple collective variables, including orientational, conformational degrees of freedom, are required to be considered in addition to the distance from the membrane center along the membrane normal. Although computationally demanding, this method can provide significant insights into the molecular mechanisms of permeation for molecules of medical and pharmaceutical importance. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

Simulation of lipid bilayer self-assembly using all-atom lipid force fields

In this manuscript we expand significantly on our earlier communication by investigating the bilayer self-assembly of eight different types of phospholipids in unbiased molecular dynamics (MD) simulations using three widely used all-atom lipid force fields. Irrespective of the underlying force field, the lipids are shown to spontaneously form stable lamellar bilayer structures within 1 microsecond, the majority of which display properties in satisfactory agreement with the experimental data. The lipids self-assemble via the same general mechanism, though at formation rates that differ both between lipid types, force fields and even repeats on the same lipid/force field combination. In addition to zwitterionic phosphatidylcholine (PC) and phosphatidylethanolamine (PE) lipids, anionic phosphatidylserine (PS) and phosphatidylglycerol (PG) lipids are represented. To our knowledge this is the first time bilayer self-assembly of phospholipids with negatively charged head groups is demonstrated in all-atom MD simulations.

Interaction of peptides with cell membranes: insights from molecular modeling

The investigation of the interaction of peptides with cell membranes is the focus of active research. It can enhance the understanding of basic membrane functions such as membrane transport, fusion, and signaling processes, and it may shed light on potential applications of peptides in biomedicine. In this review, we will present current advances in computational studies on the interaction of different types of peptides with the cell membrane. Depending on the properties of the peptide, membrane, and external environment, the peptide-membrane interaction shows a variety of different forms. Here, on the basis of recent computational progress, we will discuss how different peptides could initiate membrane pores, translocate across the membrane, induce membrane endocytosis, produce membrane curvature, form fibrils on the membrane surface, as well as interact with functional membrane proteins. Finally, we will present a conclusion summarizing recent progress and providing some specific insights into future developments in this field.

Computer simulations of lung surfactant

Lung surfactant lines the gas-exchange interface in the lungs and reduces the surface tension, which is necessary for breathing. Lung surfactant consists mainly of lipids with a small amount of proteins and forms a monolayer at the air-water interface connected to bilayer reservoirs. Lung surfactant function involves transfer of material between the monolayer and bilayers during the breathing cycle. Lipids and proteins are organized laterally in the monolayer; selected species are possibly preferentially transferred to bilayers. The complex 3D structure of lung surfactant and the exact roles of lipid organization and proteins remain important goals for research. We review recent simulation studies on the properties of lipid monolayers, monolayers with phase coexistence, monolayer-bilayer transformations, lipid-protein interactions, and effects of nanoparticles on lung surfactant. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

Molecular mechanisms of ionic liquid cytotoxicity probed by an integrated experimental and computational approach

Ionic liquids (ILs) are salts that remain liquid down to low temperatures, and sometimes well below room temperature. ILs have been called "green solvents" because of their extraordinarily low vapor pressure and excellent solvation power, but ecotoxicology studies have shown that some ILs exhibit greater toxicity than traditional solvents. A fundamental understanding of the molecular mechanisms responsible for IL toxicity remains elusive. Here we show that one mode of IL toxicity on unicellular organisms is driven by swelling of the cell membrane. Cytotoxicity assays, confocal laser scanning microscopy, and molecular simulations reveal that IL cations nucleate morphological defects in the microbial cell membrane at concentrations near the half maximal effective concentration (EC50) of several microorganisms. Cytotoxicity increases with increasing alkyl chain length of the cation due to the ability of the longer alkyl chain to more easily embed in, and ultimately disrupt, the cell membrane.

Force Field Development for Lipid Membrane Simulations

With the rapid development of computer power and wide availability of modelling software computer simulations of realistic models of lipid membranes, including their interactions with various molecular species, polypeptides and membrane proteins have become feasible for many research groups. The crucial issue of the reliability of such simulations is the quality of the force field, and many efforts, especially in the latest several years, have been devoted to parametrization and optimization of the force fields for biomembrane modelling. In this review, we give account of the recent development in this area, covering different classes of force fields, principles of the force field parametrization, comparison of the force fields, and their experimental validation. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

Prediction of Surface and Bulk Partition of Nonionic Surfactants Using Free Energy Calculations

The surface-bulk partition of nonionic surfactants was predicted by calculating chemical potential differences using two simulation boxes representing the surface and the bulk phases separately. A published coarse grained force field was modified and validated for this application. Thermodynamic integration (TI) was applied to compute excess chemical potentials. The high concentration surface was stabilized by applying an external harmonic potential, and the bulk was treated as ideal solution, which was confirmed by simulation using a lattice model at conditions near the critical micelle concentration (∼10(-5) mol/L). Based on the calculated chemical potential differences with precision of ca. 1 kJ/mol, the equilibria of surface-bulk concentration was predicted well in comparison with the experimental data.

All-atom lipid bilayer self-assembly with the AMBER and CHARMM lipid force fields

This communication reports the first example of spontaneous lipid bilayer formation in unbiased all-atom molecular dynamics (MD) simulations. Using two different lipid force fields we show simulations started from random mixtures of lipids and water in which four different types of phospholipids self-assemble into organized bilayers in under 1 microsecond.

Molecular dynamics simulations and neutron reflectivity as an effective approach to characterize biological membranes and related macromolecular assemblies

In combination with other spectroscopy, microscopy, and scattering techniques, neutron reflectivity is a powerful tool to characterize biological systems. Specular reflection of neutrons provides structural information at the nanometer and subnanometer length scales, probing the composition and organization of layered materials. Currently, analysis of neutron reflectivity data involves several simplifying assumptions about the structure of the sample under study, affecting the extraction and interpretation of information from the experimental data. Computer simulations can be used as a source of structural and dynamic data with atomic resolution. We present a novel tool to compare the structural properties determined by neutron reflectivity experiments with those obtained from molecular simulations. This tool allows benchmarking the ability of molecular dynamics simulations to reproduce experimental data, but it also promotes unbiased interpretation of experimentally determined quantities. Two application examples are presented to illustrate the capabilities of the new tool. The first example is the generation of reflectivity profiles for a 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) lipid bilayer from molecular dynamics simulations using data from both atomistic and coarse-grained models, and comparison with experimentally measured data. The second example is the calculation of lipid volume changes with temperature and composition from all atoms simulations of single and mixed 1,2-di-palmitoyl-sn-glycero-3-phosphocholine (DOPC) and 1,2-dihexadecanoyl-sn-glycero-3-phosphocholine (DPPC) bilayers.

CHARMM-GUI Martini Maker for Coarse-Grained Simulations with the Martini Force Field

Coarse-grained simulations are widely used to study large biological systems. Nonetheless, building such simulation systems becomes nontrivial, especially when membranes with various lipid types are involved. Taking advantage of the frameworks in all-atom CHARMM-GUI modules, we have developed CHARMM-GUI Martini Maker for building solution, micelle, bilayer, and vesicle systems as well as systems with randomly distributed lipids using the Martini force field. Martini Maker supports 82 lipid types and different flavors of the Martini force field, including polar and nonpolar Martini, Dry Martini, and ElNeDyn (an elastic network model for proteins). The qualities of the systems generated by Martini Maker are validated by simulations of various examples involving proteins and lipids. We expect Martini Maker to be a useful tool for modeling large, complicated biomolecular systems in a user-friendly way.

Biomembranes in atomistic and coarse-grained simulations

The architecture of biological membranes is tightly coupled to the localization, organization, and function of membrane proteins. The organelle-specific distribution of lipids allows for the formation of functional microdomains (also called rafts) that facilitate the segregation and aggregation of membrane proteins and thus shape their function. Molecular dynamics simulations enable to directly access the formation, structure, and dynamics of membrane microdomains at the molecular scale and the specific interactions among lipids and proteins on timescales from picoseconds to microseconds. This review focuses on the latest developments of biomembrane force fields for both atomistic and coarse-grained molecular dynamics (MD) simulations, and the different levels of coarsening of biomolecular structures. It also briefly introduces scale-bridging methods applicable to biomembrane studies, and highlights selected recent applications.

Three Dimensional Nano "Langmuir Trough" for Lipid Studies

A three-dimensional-phospholipid monolayer with tunable molecular structure was created on the surface of oil nanodroplets from a mixture of phospholipids, oil, and water. This simple nanoemulsion preparation technique generates an in situ prepared membrane model system with controllable molecular surface properties that resembles a lipid droplet. The molecular interfacial structure of such a nanoscopic system composed of hexadecane, 1,2-dihexadecanoyl-sn-glycero-3-phosphocholine (DPPC), and water was determined using vibrational sum frequency scattering and second harmonic scattering techniques. The droplet surface structure of DPPC can be tuned from a tightly packed liquid condensed phase like monolayer to a more dilute one that resembles the liquid condensed/liquid expanded coexistence phase by varying the DPPC/oil/water ratio. The tunability of the chemical structure, the high surface-to-volume ratio, and the small sample volume make this system an ideal model membrane for biochemical research.

All-atom molecular dynamics simulations of an artificial sodium channel in a lipid bilayer: the effect of water solvation/desolvation of the sodium ion

All-atom molecular dynamics is used to investigate the transport of Na(+) across a 1,2-dioleoyl-sn-glycero-3-phosphocholine lipid bilayer facilitated by a diazacrown hydraphile. Specifically, the free energy of Na(+) passing through the bilayer is calculated using the adaptive biasing force method to study the free energy associated with the increase in Na(+) transport in the presence of the hydraphile molecule. The results show that water interaction greatly influences Na(+) transport through the lipid bilayer as water is pulled through the bilayer with Na(+) forming a water channel. The hydraphile causes a reduction in the free energy barrier for the transport of Na(+) through the head group part of the lipid bilayer since it complexes the Na(+) reducing the necessity for water to be complexed and, therefore, dragged through with Na(+), an energetically unfavorable process. The free energy associated with Na(+) being desolvated within the bilayer is significantly decreased in the presence of the hydraphile molecule; the hydraphile increases the number of solvation states of Na(+) that can be adopted, and this increase in the number of available configurations provides an entropic explanation for the success of the hydraphile.

Coarse-grain simulations of skin ceramide NS with newly derived parameters clarify structure of melted phase

Ceramides are lipids that are involved in numerous biologically important structures (e.g., the stratum corneum and ceramide-rich platforms) and processes (e.g., signal transduction and membrane fusion), but their behavior is not fully understood. We report coarse-grain force field parameters for N-lignocerylsphingosine (ceramide NS, also known as ceramide 2) that are consistent with the Martini force field. These parameters were optimized for simulations in the gel phase and validated against atomistic simulations. Coarse-grained simulations with our parameters provide areas per lipid, membrane thicknesses, and electron density profiles that are in good agreement with atomistic simulations. Properties of the simulated membranes are compared with available experimental data. The obtained parameters were used to model the phase behavior of ceramide NS as a function of temperature and hydration. At low water content and above the main phase transition temperature, the bilayer melts into an irregular phase, which may correspond to the unstructured melted-chain phase observed in X-ray diffraction experiments. The developed parameters also reproduce the extended conformation of ceramide, which may occur in the stratum corneum. The parameters presented herein will facilitate studies on important complex functional structures such as the uppermost layer of the skin and ceramide-rich platforms in phospholipid membranes.

Coarse-grained force fields for molecular simulations

Molecular dynamics (MD) simulations at the atomic scale are a powerful tool to study the structure and dynamics of model biological systems. However, because of their high computational cost, the time and length scales of atomistic simulations are limited. Biologically important processes, such as protein folding, ion channel gating, signal transduction, and membrane remodeling, are difficult to investigate using atomistic simulations. Coarse-graining reduces the computational cost of calculations by reducing the number of degrees of freedom in the model, allowing simulations of larger systems for longer times. In the first part of this chapter we review briefly some of the coarse-grained models available for proteins, focusing on the specific scope of each model. Then we describe in more detail the MARTINI coarse-grained force field, and we illustrate how to set up and run a simulation of a membrane protein using the Gromacs software package. We explain step-by-step the preparation of the protein and the membrane, the insertion of the protein in the membrane, the equilibration of the system, the simulation itself, and the analysis of the trajectory.

Computer simulations of phase separation in lipid bilayers and monolayers

Studying phase coexistence in lipid bilayers and monolayers is important for understanding lipid-lipid interactions underlying lateral organization in biological membranes. Computer simulations follow experimental approaches and use model lipid mixtures of simplified composition. Atomistic simulations give detailed information on the specificity of intermolecular interactions, while coarse-grained simulations achieve large time and length scales and provide a bridge towards state-of-the-art experimental techniques. Computer simulations allow characterizing the structure and composition of domains during phase transformations at Angstrom and picosecond resolution, and bring new insights into phase behavior of lipid membranes.

Instability of carbon nanoparticles interacting with lipid bilayers

For nanoparticles penetrating biological tissue, modelling indicates that without external forces, carbon nanoparticles will remain trapped in lipid bilayers.

TCR triggering by pMHC ligands tethered on surfaces via poly(ethylene glycol) depends on polymer length

Antigen recognition by T cells relies on the interaction between T cell receptor (TCR) and peptide-major histocompatibility complex (pMHC) at the interface between the T cell and the antigen presenting cell (APC). The pMHC-TCR interaction is two-dimensional (2D), in that both the ligand and receptor are membrane-anchored and their movement is limited to 2D diffusion. The 2D nature of the interaction is critical for the ability of pMHC ligands to trigger TCR. The exact properties of the 2D pMHC-TCR interaction that enable TCR triggering, however, are not fully understood. Here, we altered the 2D pMHC-TCR interaction by tethering pMHC ligands to a rigid plastic surface with flexible poly(ethylene glycol) (PEG) polymers of different lengths, thereby gradually increasing the ligands' range of motion in the third dimension. We found that pMHC ligands tethered by PEG linkers with long contour length were capable of activating T cells. Shorter PEG linkers, however, triggered TCR more efficiently. Molecular dynamics simulation suggested that shorter PEGs exhibit faster TCR binding on-rates and off-rates. Our findings indicate that TCR signaling can be triggered by surface-tethered pMHC ligands within a defined 3D range of motion, and that fast binding rates lead to higher TCR triggering efficiency. These observations are consistent with a model of TCR triggering that incorporates the dynamic interaction between T cell and antigen-presenting cell.

[70]Fullerenes assist the formation of phospholipid bicelles at low lipid concentrations

The incorporation of neutral [70]fullerenes (C70) led to bicelle formation in a relatively low lipid concentration range from neutral lipid mixtures (DMPC/DHPC). Furthermore, C70 addition resulted in the formation of large bicelles with a radius of ca. 100 nm, in contrast to C70-free bicelles that were formed from anionic lipid mixtures (DMPC/DHPC/DMPG). The stabilization of these bicelles was attributed to C70 incorporation into the membranes.

Interactions of peripheral proteins with model membranes as viewed by molecular dynamics simulations

Many cellular signalling and related events are triggered by the association of peripheral proteins with anionic lipids in the cell membrane (e.g. phosphatidylinositol phosphates or PIPs). This association frequently occurs via lipid-binding modules, e.g. pleckstrin homology (PH), C2 and four-point-one, ezrin, radixin, moesin (FERM) domains, present in peripheral and cytosolic proteins. Multiscale simulation approaches that combine coarse-grained and atomistic MD simulations may now be applied with confidence to investigate the molecular mechanisms of the association of peripheral proteins with model bilayers. Comparisons with experimental data indicate that such simulations can predict specific peripheral protein–lipid interactions. We discuss the application of multiscale MD simulation and related approaches to investigate the association of peripheral proteins which contain PH, C2 or FERM-binding modules with lipid bilayers of differing phospholipid composition, including bilayers containing multiple PIP molecules.

Interleaflet sliding in lipidic bilayers under shear flow: comparison of the gel and fluid phases using reversed non-equilibrium molecular dynamics simulations

The friction between two rubbing surfaces lubricated by water can be diminished if they are coated with phospholipidic bilayers or brushes of polyelectrolytes. In the case of a coating by lipid membranes, the friction is lower when the lipids are in the gel phase rather than in the liquid phase. We investigated the response of fluid or gel bilayers to a mechanical load or under shear using non-equilibrium molecular dynamics simulations (NEMD) to understand whether this difference could come from intermonolayer sliding. The system is composed of a single fully hydrated bilayer of coarse grained phospholipids under a parallel shear with vorticity parallel to the bilayer. In both the liquid and the gel phases, an intermonolayer slip was measured in the velocity profile. In the liquid phase this slip is proportional to the shear stress. In the tilted gel phase of our model the stress is not systematically linear and relaxes differently when the shear is in the direction of the tilt or perpendicular to it. The impact of surface tension (or load) on the friction is different for the liquid and gel phases, but grossly the slip remains of the same order of magnitude.

Free energy analysis along the stalk mechanism of membrane fusion

The free energy profile of the stalk model of membrane fusion has been calculated using coarse-grained molecular dynamics simulations. The proposed method guides the lipid configuration using a guiding wall potential to make the transition from two apposed membranes to a stalk and a fusion pore. The free energy profile is obtained with a thermodynamic integration scheme using the mean force working on the guiding wall as a response of the system. We applied the method to two apposed flat bilayers composed of dioleoyl phosphatidylethanolamine/dioleoyl phosphatidylcholine expanding over the simulation box under the periodic boundary conditions. The two transition states are identified as pre-stalk and pre-pore states. The free energy barrier for the latter is confirmed to be in good agreement with that estimated by the pulling method. The present method provides a practical way to calculate the free energy profile along the stalk mechanism.

The power of coarse graining in biomolecular simulations

Computational modeling of biological systems is challenging because of the multitude of spatial and temporal scales involved. Replacing atomistic detail with lower resolution, coarse grained (CG), beads has opened the way to simulate large-scale biomolecular processes on time scales inaccessible to all-atom models. We provide an overview of some of the more popular CG models used in biomolecular applications to date, focusing on models that retain chemical specificity. A few state-of-the-art examples of protein folding, membrane protein gating and self-assembly, DNA hybridization, and modeling of carbohydrate fibers are used to illustrate the power and diversity of current CG modeling.

A guiding potential method for evaluating the bending rigidity of tensionless lipid membranes from molecular simulation

A new method is proposed to estimate the bending rigidity of lipid membranes from molecular dynamics simulations. An external cylindrical guiding potential is used to impose a sinusoidal deformation to a planar membrane. The bending rigidity is obtained from the mean force acting on the cylinder by calibrating against a discretized Helfrich model that accounts for thermal fluctuations of the membrane surface. The method has been successfully applied to a dimyristoyl phosphatidylcholine bilayer simulated with a coarse-grained model. A well-converged bending rigidity was obtained for the tension-free membrane and showed reasonable agreement with that obtained from the height fluctuation spectrum.