Simulations Of Bilayers Research Articles

Lipid Packing Defects are Necessary and Sufficient for Membrane Binding of α-Synuclein.

α-Synuclein (αSyn), an intrinsically disordered protein implicated in Parkinson's disease, is potentially thought to initiate aggregation through binding to cellular membranes. Previous studies have suggested that anionic membrane charge is necessary for this binding. However, these studies largely focus on unmodified αSyn, while nearly all αSyn in the body is N-terminally acetylated (NTA). NTA dramatically shifts the narrative by diminishing αSyn's reliance on anionic charge for membrane binding. Instead, we demonstrate that membrane packing defects are the dominant forces driving NTA-αSyn interactions, challenging the long-standing paradigm that anionic membranes are essential for αSyn binding. Using fluorescence microscopy and circular dichroism spectroscopy, we monitored the binding of NTA-αSyn to reconstituted membrane surfaces with different lipid compositions. Phosphatidylcholine and phosphatidylserine concentrations were varied to control surface charge, while phospholipid tail unsaturation and methylation were varied to control lipid packing. All-atom molecular dynamics simulations of lipid bilayers supported the observation that membrane packing defects are necessary for NTA-αSyn binding and that defect-rich membranes are sufficient for NTA-αSyn binding regardless of membrane charge. We further demonstrated that this affinity for membrane defects persisted in reconstituted, cholesterol-containing membranes that mimicked the physiological lipid composition of synaptic vesicles. Increasing phospholipid unsaturation in these mimics led to more membrane packing defects and a corresponding increase in NTA-αSyn binding. Altogether, our results point to a mechanism for the regulation of NTA-αSyn binding in biological membranes that extends beyond phospholipid charge to the structural properties of the lipids themselves.

Split Membrane: A New Model to Accelerate All-Atom MD Simulation of Phospholipid Bilayers.

All-atom molecular dynamics simulations are powerful tools for studying cell membranes and their interactions with proteins and other molecules. However, these processes occur on time scales determined by the diffusion rate of phospholipids, which are challenging to achieve in all-atom models. Here, we present a new all-atom model that accelerates lipid diffusion by splitting phospholipid molecules into head and tail groups. The bilayer structure is maintained by using external lateral potentials, which compensate for the lipid split. This split model enhances lateral lipid diffusion more than ten times, allowing faster and cheaper equilibration of large systems with different phospholipid types. The current model has been tested on membranes containing PSM, POPC, POPS, POPE, POPA, and cholesterol. We have also evaluated the interaction of the split model membranes with the Disheveled DEP domain and amphiphilic helix motif of the transcriptional repressor Opi1 as representative of peripheral proteins as well as the dimeric fragment of the epidermal growth factor receptor transmembrane domain and the Human A2A Adenosine of G protein-coupled receptors as representative of transmembrane proteins. The split model can predict the interaction sites of proteins and their preferred phospholipid type. Thus, the model could be used to identify lipid binding sites and equilibrate large membranes at an affordable computational cost.

Dynamic properties of isotropic DMPC/DHPC bicelles: Insights from solution NMR and MD simulations.

Bicelles, an artificial disk-shaped lipid bilayer, are commonly used for the structural and functional characterization of membrane-bound proteins in an environment similar to that in intracellular membranes. Because the dynamics of the lipids that constitute bicelles exert a significant impact on the structure and function of the inserted proteins, we investigated the mobility of lipid molecules in bicelles composed of DMPC (14:0PC) and DHPC (06:0PC) using solution NMR and MD calculations. 13C R1 relaxation experiments for the acyl groups demonstrated that increasing bicelle sizes limit the rotational diffusion of acyl chain H-C bonds in DMPC. Such a limited local motion around H-C bonds was also predicted in the MD simulations of DMPC bilayers with decreased area per lipid, indicating that the limited mobility of the hydrophobic core in larger bicelles is due to the tighter lipid packing. The downfield shifts of the 13C NMR signals of the acyl groups supported the restricted mobility, corresponding to the conformational changes of the acyl chains from the flexible gauche rotamers to the less mobile trans rotamers with increase in bicelle size. These data suggest that larger bicelles pack lipids more densely, leading to increased trans conformation of the acyl chains and, consequently, less lipid motility, which can dynamically modulate the structure and function of membrane-bound proteins inserted into the bicelles.

Fully atomistic molecular dynamics modeling of photoswitchable azo-PC lipid bilayers: structure, mechanical properties, and drug permeation.

Phospholipid based vesicles called liposomes are commonly used as packaging in advanced drug delivery applications. Stimuli-responsive liposomes have been designed to release their contents under certain conditions, for example through heating or illumination. However, in the case of photosensitive liposomes based on azo-PC, namely phosphatidylcholine lipids with azobenzene incorporated into one of the two lipid tails, the release mechanism is not known. Here we show, using fully-atomistic molecular dynamics simulations of pure azo-PC bilayers, that drug permeation through the bilayer is driven by a light-induced gel-to-liquid lipid phase transition that softens the membrane bending rigidity by an order of magnitude, increases the area per lipid, and decreases the membrane thickness. Furthermore, using phenol as a model drug, we quantified its translocation free energy and its ability to cross the bilayer as a result of a chemical potential gradient induced through a double-bilayer simulation unit cell. The molecular level structural and dynamic information obtained in this study should be of help in designing new azo-PC based liposomes.

Studying the Effects of Dissolved Noble Gases and High Hydrostatic Pressure on the Spherical DOPC Bilayer Using Molecular Dynamic Simulations.

Fine-grained molecular dynamics simulations have been conducted to depict lipid objects enclosed in water and interacting with a series of noble gases dissolved in the medium. The simple point-charge (SPC) water system, featuring a boundary composed of 1,2-Dioleoyl-sn-glycero-3-phosphocholine (DOPC) molecules, maintained stability throughout the simulation under standard conditions. This allowed for the accurate modeling of the effects of hydrostatic pressure at an ambient pressure of 25 bar. The chosen pressure references the 240 m depth of seawater: the horizon frequently used by commercial divers, who comprise the primary patient population of the neurological complication of inert gas narcosis and the consequences of high-pressure neurological syndrome. To quantify and validate the neurological effects of noble gases and discriminate them from high hydrostatic pressure, we reduced the dissolved gas molar concentration to 1.5%, three times smaller than what we previously tested for the planar bilayer (3.5%). The nucleation and growth of xenon, argon and neon nanobubbles proved consistent with the data from the planar bilayer simulations. On the other hand, hyperbaric helium induces only a residual distorting effect on the liposome, with no significant condensed gas fraction observed within the hydrophobic core. The bubbles were distributed over a large volume-both in the bulk solvent and in the lipid phase-thereby causing substantial membrane distortion. This finding serves as evidence of the validity of the multisite distortion hypothesis for the neurological effect of inert gases at high pressure.

Coupling deformation analysis of self-morphing bilayers with mismatch strain

Self-morphing bilayers driven by the inner mismatch strain have inspired numerous actuators and soft robots in the four-dimensional (4D) printing field. Analyzing their coupling deformation is crucial for the design of actuators, but current methods face inefficiencies or accuracy issues. This study provides convenient and efficient approaches, from both deformation theory and model aspects, to analyze the coupling deformations of self-morphing bilayers. First, we derived a simplified energy expression for self-morphing bilayers to reflect their deformation process. This concise form permits direct derivation of some deformation theories for self-morphing bilayers, thereby significantly simplifying the complexity of deformation analysis. Second, we established the deformation simulation model of bilayers by discretizing and minimizing the simplified energy expression. The checkerboard pattern method, a new discrete differential geometric (DDG) approach, was introduced and improved to discretize the derived energy expression. The Levenberg-Marquardt (LM) algorithm was used to establish the deformation simulation model for bilayers by minimizing the discrete energy function. Further, a series of numerical simulations of bending bilayers were conducted. Results validated that the proposed model achieved a faster simulation speed compared to the finite element method (FEM)-based model in our previous work, and was more accurate for finite strain situations. Finally, three bio-inspired bilayers and a bio-inspired gripper (successfully grasped a sphere) were designed using the resulting model, demonstrating its practical application capabilities. These proposed methods hold substantial potential to facilitate the advancement of the design field of 4D printing and actuators.

Ion-Induced PIP2 Clustering with Martini3: Modification of Phosphate-Ion Interactions and Comparison with CHARMM36.

Phosphatidylinositol 4,5-bisphosphate (PIP2) is a critical lipid for cellular signaling. The specific phosphorylation of the inositol ring controls protein binding as well as clustering behavior. Two popular models to describe ion-mediated clustering of PIP2 are Martini3 (M3) and CHARMM36 (C36). Molecular dynamics simulations of PIP2-containing bilayers in solutions of potassium chloride, sodium chloride, and calcium chloride, and at two different resolutions are performed to understand the aggregation and the model parameters that drive it. The average M3 clusters of PIP2 in bilayers of 1-palmitoyl-2-oleoyl-glycero-3-phosphocholine and PIP2 bilayers in the presence of K+, Na+, or Ca2+ contained 2.2, 2.6, and 6.4 times more PIP2 than C36 clusters, respectively. Indeed, the Ca2+-containing systems often formed a single large aggregate. Reparametrization of the M3 ion-phosphate Lennard-Jones interaction energies to reproduce experimental osmotic pressure of sodium dimethyl phosphate (DMP), K[DMP], and Ca[DMP]2 solutions, the same experimental target as C36, yielded comparably sized PIP2 clusters for the two models. Furthermore, C36 and the modified M3 predict similar saturation of the phosphate groups with increasing Ca2+, although the coarse-grained model does not capture the cooperativity between K+ and Ca2+. This characterization of the M3 behavior in the presence of monovalent and divalent ions lays a foundation to study cation/protein/PIP2 clustering.

Simulating asymmetric membranes using P21 periodic boundary conditions.

Molecular dynamics (MD) simulations of symmetric lipid bilayers are now well established, while those of asymmetric ones are considerably less developed. This disjunction arises in part because the surface tensions of leaflets in asymmetric bilayers can differ (unlike those of symmetric ones), and there is no simple way to determine them without assumptions. This chapter describes the use of P21 periodic boundary conditions (PBC), which allow lipids to switch leaflets, to generate asymmetric bilayers under the assumption of equal chemical potentials of lipids in opposing leaflets. A series of examples, ranging from bilayers with one lipid type to those with peptides and proteins, provides a guide for the use of P21 PBC. Critical properties of asymmetric membranes, such as spontaneous curvature, are highly sensitive to differences in the leaflet surface tensions (or differential stress), and equilibration with P21 PBC substantially reduces differential stress of asymmetric bilayers assembled with surface area-based methods. Limitations of the method are discussed. Technically, the nonstandard unit cell is difficult to parallelize and to incorporate restraints. Inherently, the assumption of equal chemical potentials, and therefore the method itself, is not applicable to all target systems. Despite these limitations, it is argued that P21 simulations should be considered when designing equilibration protocols for MD studies of most asymmetric membranes.

Binary bilayer simulations for partitioning within membranes.

Membrane proteins (MPs) often show preference for one phase over the other, which is characterized by the partition coefficient, Kp. The physical mechanisms underlying Kp have been only inferred indirectly from experiments due to the unavailability of detailed structures and compositions of ordered phases. Molecular dynamics (MD) simulations can complement these details and thus, in principle, provide further insights into the partitioning of MPs between two phases. However, the application of MD has remained difficult due to long time scales required for equilibration and large system size for the phase stability, which have not been fully resolved even in free energy simulations. This chapter describes the recently developed binary bilayer simulation method, where the membrane is composed of two laterally attached membrane patches. The binary bilayer system (BBS) is designed to preserve the lateral packing of both phases in a significantly smaller size compared to that required for macroscopic phase separation. These characteristics are advantageous in partitioning simulations, as the length scale for diffusion across the system can be significantly smaller. Hence the BBS can be efficiently employed in both conventional MD and free energy simulations, though sampling in ordered phases remains difficult due to slow diffusion. Development of efficient lipid swapping methods and its combination with the BBS would be a useful approach for partitioning in coexisting phases.

Constructing bilayer and volumetric atrial models at scale.

To enable large in silico trials and personalized model predictions on clinical timescales, it is imperative that models can be constructed quickly and reproducibly. First, we aimed to overcome the challenges of constructing cardiac models at scale through developing a robust, open-source pipeline for bilayer and volumetric atrial models. Second, we aimed to investigate the effects of fibres, fibrosis and model representation on fibrillatory dynamics. To construct bilayer and volumetric models, we extended our previously developed coordinate system to incorporate transmurality, atrial regions and fibres (rule-based or data driven diffusion tensor magnetic resonance imaging (MRI)). We created a cohort of 1000 biatrial bilayer and volumetric models derived from computed tomography (CT) data, as well as models from MRI, and electroanatomical mapping. Fibrillatory dynamics diverged between bilayer and volumetric simulations across the CT cohort (correlation coefficient for phase singularity maps: left atrial (LA) 0.27 ± 0.19, right atrial (RA) 0.41 ± 0.14). Adding fibrotic remodelling stabilized re-entries and reduced the impact of model type (LA: 0.52 ± 0.20, RA: 0.36 ± 0.18). The choice of fibre field has a small effect on paced activation data (less than 12 ms), but a larger effect on fibrillatory dynamics. Overall, we developed an open-source user-friendly pipeline for generating atrial models from imaging or electroanatomical mapping data enabling in silico clinical trials at scale (https://github.com/pcmlab/atrialmtk).

Dynamics behavior of PE and PET oligomers in lipid bilayer simulations

In recent years many investigators have been concerned about the toxicity and potential health hazards of micro- and nanoplastics. However, we are still lacking a good understanding of the methods of their transport into the human body and subsequently within cells. This is especially true at the lower nanometer scale; these particles are potentially more dangerous than their micrometer counterparts due to their easier permeation into cells. In this study we used both unbiased molecular dynamics simulations and steered umbrella sampling simulations to explore the interactions of polyethylene terephthalate (PET) and polyethylene (PE) oligomers in phospholipid bilayers. Our simulations revealed that the bilayers did not represent significant energy barriers to the small oligomers; not only did they readily enter the cell membrane but they also became concentrated into specific parts of the membrane. The larger PET tetramers exhibited a strong aggregation in water but were the least likely to permeate through or into the membranes. It is possible that PE monomers and tetramers can become concentrated into membranes while PET monomers are more likely to pass through or concentrate just inside the membrane surface. Passive transport of microplastics into cells is, however, likely limited to particles of a few nanometers in diameter.

Molecular dynamics study of the mechanical properties of drug loaded model systems: A comparison of a polymersome with a bilayer.

In this work we implement a new methodology to study structural and mechanical properties of systems having spherical and planar symmetries throughout Molecular Dynamics simulations. This methodology is applied here to a drug delivery system based in polymersomes, as an example. The chosen model drug was the local anesthetic prilocaine due to previous parameterization within the used coarse grain scheme. In our approach, mass density profiles (MDPs) are used to obtain key structural parameters of the systems, and pressure profiles are used to estimate the curvature elastic parameters. The calculation of pressure profiles and radial MPDs required the development of specific methods, which were implemented in an in-house built version of the GROMACS 2018 code. The methodology presented in this work is applied to characterize poly(ethylene oxide)-poly(butadiene) polymersomes and bilayers loaded with the model drug prilocaine. Our results show that structural properties of the polymersome membrane could be obtained from bilayer simulations, with significantly lower computational cost compared to whole polymersome simulations, but the bilayer simulations are insufficient to get insights on their mechanical aspects, since the elastic parameters are canceled out for the complete bilayer (as consequence of the symmetry). The simulations of entire polymersomes, although more complex, offer a complementary approach to get insights on the mechanical behavior of the systems.

Understanding the Mechanical Properties of Ultradeformable Liposomes Using Molecular Dynamics Simulations.

Improving drug delivery efficiency to solid tumor sites is a central challenge in anticancer therapeutic research. Our previous experimental study (Guo et al., Nat. Commun. 2018, 9, 130) showed that soft, elastic liposomes had increased uptake and accumulation in cancer cells and tumors in vitro and in vivo respectively, relative to rigid particles. As a first step toward understanding how liposomes' molecular structure and composition modulates their elasticity, we performed all-atom and coarse-grained classical molecular dynamics (MD) simulations of lipid bilayers formed by mixing a long-tailed unsaturated phospholipid with a short-tailed saturated lipid with the same headgroup. The former types of phospholipids considered were 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) and 1,2-dipalmitoleoyl-sn-glycero-3-phosphocholine (termed here DPMPC). The shorter saturated lipids examined were 1,2-diheptanoyl-sn-glycero-3-phosphocholine (DHPC), 1,2-didecanoyl-sn-glycero-3-phosphocholine (DDPC), 1,2-dilauroyl-sn-glycero-3-phosphocholine (DLPC), and 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC). Several lipid concentrations and surface tensions were considered. Our results show that DOPC or DPMPC systems having 25-35 mol % of the shortest lipids DHPC or DDPC are the least rigid, having area compressibility moduli KA that are ∼10% smaller than the values observed in pure DOPC or DPMPC bilayers. These results agree with experimental measurements of the stretching modulus and lysis tension in liposomes with the same compositions. These mixed systems also have lower areas per lipid and form more uneven x-y interfaces with water, the tails of both primary and secondary lipids are more disordered, and the terminal methyl groups in the tails of the long lipid DOPC or DPMPC wriggle more in the vertical direction, compared to pure DOPC or DPMPC bilayers or their mixtures with the longer saturated lipid DLPC or DMPC. These observations confirm our hypothesis that adding increasing concentrations of the short unsaturated lipid DHPC or DDPC to DOPC or DPMPC bilayers alters lipid packing and thus makes the resulting liposomes more elastic and less rigid. No formation of lipid nanodomains was noted in our simulations, and no clear trends were observed in the lateral diffusivities of the lipids as the concentration, type of secondary lipid, and surface tension were varied.

Cellular excitability and ns-pulsed electric fields: Potential involvement of lipid oxidation in the action potential activation

Recent studies showed that nanosecond pulsed electric fields (nsPEFs) can activate voltage-gated ion channels (VGICs) and trigger action potentials (APs) in excitable cells. Under physiological conditions, VGICs’ activation takes place on time scales of the order 10–100 µs. These time scales are considerably longer than the applied pulse duration, thus activation of VGICs by nsPEFs remains puzzling and there is no clear consensus on the mechanisms involved. Here we propose that changes in local electrical properties of the cell membrane due to lipid oxidation might be implicated in AP activation. We first use MD simulations of model lipid bilayers with increasing concentration of primary and secondary lipid oxidation products and demonstrate that oxidation not only increases the bilayer conductance, but also the bilayer capacitance. Equipped with MD-based characterization of electrical properties of oxidized bilayers, we then resort to AP modelling at the cell level with Hodgkin-Huxley-type models. We confirm that a local change in membrane properties, particularly the increase in membrane conductance, due to formation of oxidized membrane lesions can be high enough to trigger an AP, even when no external stimulus is applied. However, excessive accumulation of oxidized lesions (or other conductive defects) can lead to altered cell excitability.

Extension of the iSoLF implicit-solvent coarse-grained model for multicomponent lipid bilayers.

iSoLF is a coarse-grained (CG) model for lipid molecules with the implicit-solvent approximation used in molecular dynamics (MD) simulations of biological membranes. Using the original iSoLF (iSoLFv1), MD simulations of lipid bilayers consisting of either POPC or DPPC and these bilayers, including membrane proteins, can be performed. Here, we improve the original model, explicitly treating the electrostatic interactions between different lipid molecules and adding CG particle types. As a result, the available lipid types increase to 30. To parameterize the potential functions of the new model, we performed all-atom MD simulations of each lipid at three different temperatures using the CHARMM36 force field and the modified TIP3P model. Then, we parameterized both the bonded and non-bonded interactions to fit the area per lipid and the membrane thickness of each lipid bilayer by using the multistate Boltzmann Inversion method. The final model reproduces the area per lipid and the membrane thickness of each lipid bilayer at the three temperatures. We also examined the applicability of the new model, iSoLFv2, to simulate the phase behaviors of mixtures of DOPC and DPPC at different concentrations. The simulation results with iSoLFv2 are consistent with those using Dry Martini and Martini 3, although iSoLFv2 requires much fewer computations. iSoLFv2 has been implemented in the GENESIS MD software and is publicly available.

Cholesterol Biases the Conformational Landscape of the Chemokine Receptor CCR3: A MAS SSNMR-Filtered Molecular Dynamics Study.

Cholesterol directs the pathway of ligand-induced G protein-coupled receptor (GPCR) signal transduction. The GPCR C-C motif chemokine receptor 3 (CCR3) is the principal chemotactic receptor for eosinophils, with roles in cancer metastasis and autoinflammatory conditions. Recently, we discovered a direct correlation between bilayer cholesterol and increased agonist-triggered CCR3 signal transduction. However, the allosteric molecular mechanism escalating ligand affinity and G protein coupling is unknown. To study cholesterol-guided CCR3 conformational selection, we implement comparative, objective measurement of protein architectures by scoring shifts (COMPASS) to grade model structures from molecular dynamics simulations. In this workflow, we scored predicted chemical shifts against 2-dimensional solid-state NMR 13C-13C correlation spectra of U-15N,13C-CCR3 samples prepared with and without cholesterol. Our analysis of trajectory model structures uncovers that cholesterol induces site-specific conformational restraint of extracellular loop (ECL) 2 and conserved motion in transmembrane helices and ECL3 not observed in simulations of bilayers with only phosphatidylcholine lipids. PyLipID analysis implicates direct cholesterol agency in CCR3 conformational selection and dynamics. Residue-residue contact scoring shows that cholesterol biases the conformational selection of the orthosteric pocket involving Y411.39, Y1133.32, and E2877.39. Lastly, we observe contact remodeling in activation pathway residues centered on the initial transmission switch, Na+ pocket, and R3.50 in the DRY motif. Our observations have unique implications for understanding of CCR3 ligand recognition and specificity and provide mechanistic insight into how cholesterol functions as an allosteric regulator of CCR3 signal transduction.

Drude Polarizable Lipid Force Field with Explicit Treatment of Long-Range Dispersion: Parametrization and Validation for Saturated and Monounsaturated Zwitterionic Lipids.

Accurate empirical force fields of lipid molecules are a critical component of molecular dynamics simulation studies aimed at investigating properties of monolayers, bilayers, micelles, vesicles, and liposomes, as well as heterogeneous systems, such as protein-membrane complexes, bacterial cell walls, and more. While the majority of lipid force field-based simulations have been performed using pairwise-additive nonpolarizable models, advances have been made in the development of the polarizable force field based on the classical Drude oscillator model. In the present study, we undertake further optimization of the Drude lipid force field, termed Drude2023, including improved treatment of the phosphate and glycerol linker region of PC and PE headgroups, additional optimization of the alkene group in monounsaturated lipids, and inclusion of long-range Lennard-Jones interactions using the particle-mesh Ewald method. Initial optimization targeted quantum mechanical (QM) data on small model compounds representative of the linker region. Subsequent optimization targeted QM data on larger model compounds, experimental data, and dihedral potentials of mean force from the CHARMM36 additive lipid force field using a parameter reweighting protocol. The use of both experimental and QM target data during the reweighting protocol is shown to produce physically reasonable parameters that reproduce a collection of experimental observables. Target data for optimization included surface area/lipid for DPPC, DSPC, DMPC, and DLPC bilayers and nuclear magnetic resonance (NMR) order parameters for DPPC bilayers. Validation data include prediction of membrane thickness, scattering form factors, electrostatic potential profiles, compressibility moduli, surface area per lipid, water permeability, NMR T1 relaxation times, diffusion constants, and monolayer surface tensions for a variety of saturated and unsaturated lipid mono- and bilayers. Overall, the agreement with experimental data is quite good, though the results are less satisfactory for the NMR T1 relaxation times for carbons near the ester groups. Notable improvements compared to the additive C36 force field were obtained for membrane dipole potentials, lipid diffusion coefficients, and water permeability with the exception of monounsaturated lipid bilayers. It is anticipated that the optimized polarizable Drude2023 force field will help generate more accurate molecular simulations of pure bilayers and heterogeneous systems containing membranes, advancing our understanding of the role of electronic polarization in these systems.

Moiré patterns and inversion boundaries in graphene/hexagonal boron nitride bilayers

In this paper a systematic examination of graphene/hexagonal boron nitride (g/hBN) bilayers is presented, through a recently developed two-dimensional phase field crystal model that incorporates out-of-plane deformations. The system parameters are determined by closely matching the stacking energies and heights of g/hBN bilayers to those obtained from existing quantum-mechanical density functional theory calculations. Out-of-plane deformations are shown to reduce the energies of inversion domain boundaries in hBN, and the coupling between graphene and hBN layers leads to a bilayer defect configuration consisting of an inversion boundary in hBN and a domain wall in graphene. Simulations of twisted bilayers reveal the structure, energy, and elastic properties of the corresponding moir\'e patterns and show a crossover as the misorientation angle between the layers increases from a well-defined hexagonal network of domain boundaries and junctions to smeared-out patterns. The transition occurs when the thickness of domain walls approaches the size of the moir\'e patterns and coincides with the peaks in the average von Mises and volumetric stresses of the bilayer.

Surface viscosities of lipid bilayers determined from equilibrium molecular dynamics simulations

Lipid membrane viscosity is critical to biological function. Bacterial cells grown in different environments alter their lipid composition in order to maintain a specific viscosity, and membrane viscosity has been linked to the rate of cellular respiration. To understand the factors that determine the viscosity of a membrane, we ran equilibrium all-atom simulations of single component lipid bilayers and calculated their viscosities. The viscosity was calculated via a Green-Kubo relation, with the stress-tensor autocorrelation function modeled by a stretched exponential function. By simulating a series of lipids at different temperatures, we establish the dependence of viscosity on several aspects of lipid chemistry, including hydrocarbon chain length, unsaturation, and backbone structure. Sphingomyelin is found to have a remarkably high viscosity, roughly 20 times that of DPPC. Furthermore, we find that inclusion of the entire range of the dispersion interaction increases viscosity by up to 140%. The simulated viscosities are similar to experimental values obtained from the rotational dynamics of small chromophores and from the diffusion of integral membrane proteins but significantly lower than recent measurements based on the deformation of giant vesicles.

Surface viscosities of several lipids determined from equilibrium molecular dynamics simulations.

Lipid membrane viscosity is critical to biological function. Bacterial cells grown in different environments have been shown to alter their lipid composition in order to maintain a viscosity in a narrow range,1 and membrane viscosity has been shown to control the rate of cellular respiration.2 In order to understand the factors which determine membrane viscosity, we ran equilibrium all-atom simulations of single component lipid bilayers and calculated their viscosities. The viscosity was calculated via a Green-Kubo relation, with the stress tensor autocorrelation function fit to a stretched exponential function. By simulating a series of lipids at different temperatures, we establish the dependence of viscosity on several aspects of lipid chemistry, including hydrocarbon chain length, unsaturation and backbone structure. Sphingomyelin is found to have a remarkably high viscosity, roughly twenty times that of DPPC. Including the full range of the dispersion increases the membrane viscosity by up to 140%, alleviating some of the discrepancy between simulated and experimental viscosity measurements, although the simulated viscosities remain significantly lower than recent experimental values reported for a series of lipids.3