Toward Atomistic Resolution Structure of Phosphatidylcholine Headgroup and Glycerol Backbone at Different Ambient Conditions

Interaction of MRI Contrast Agent [Gd(DOTA)]- with Lipid Membranes: A Molecular Dynamics Study

Contrast agents are important imaging probes in clinical MRI, allowing the identification of anatomic changes that otherwise would not be possible. Intensive research on the development of new contrast agents is being made to image specific pathological markers or sense local biochemical changes. The most widely used MRI contrast agents are based on gadolinium(III) complexes. Due to their very high charge density, they have low permeability through tight biological barriers such as the blood-brain barrier, hampering their application in the diagnosis of neurological disorders. In this study, we explore the interaction between the widely used contrast agent [Gd(DOTA)]- (Dotarem) and POPC lipid bilayers by means of molecular dynamics simulations. This metal complex is a standard reference where several chemical modifications have been introduced to improve key properties such as bioavailability and targeting. The simulations unveil detailed insights into the agent's interaction with the lipid bilayer, offering perspectives beyond experimental methods. Various properties, including the impact on global and local bilayer properties, were analyzed. As expected, the results indicate a low partition coefficient (KP) and high permeation barrier for this reference compound. Nevertheless, favorable interactions are established with the membrane leading to moderately long residence times. While coordination of one inner-sphere water molecule is maintained for the membrane-associated chelate, the physical-chemical attributes of [Gd(DOTA)]- as a MRI contrast agent are affected. Namely, increases in the rotational correlation times and in the residence time of the inner-sphere water are observed, with the former expected to significantly increase the water proton relaxivity. This work establishes a reference framework for the use of simulations to guide the rational design of new contrast agents with improved relaxivity and bioavailability and for the development of liposome-based formulations for use as imaging probes or theranostic agents.

Evaluating Polarizable Biomembrane Simulations against Experiments

Owing to the increase of available computational capabilities and the potential for providing a more accurate description, polarizable molecular dynamics force fields are gaining popularity in modeling biomolecular systems. It is, however, crucial to evaluate how much precision is truly gained with increasing cost and complexity of the simulation. Here, we leverage the NMRlipids open collaboration and Databank to assess the performance of available polarizable lipid models─the CHARMM-Drude and the AMOEBA-based parameters─against high-fidelity experimental data and compare them to the top-performing nonpolarizable models. While some improvement in the description of ion binding to membranes is observed in the most recent CHARMM-Drude parameters, and the conformational dynamics of AMOEBA-based parameters are excellent, the best nonpolarizable models tend to outperform their polarizable counterparts for each property we explored. The identified shortcomings range from inaccuracies in describing the conformational space of lipids to excessively slow conformational dynamics. Our results provide valuable insights for the further refinement of polarizable lipid force fields and for selecting the best simulation parameters for specific applications.

Insight into the Molecular Mechanism of Surface Interactions of Phosphatidylcholines─Langmuir Monolayer Study Complemented with Molecular Dynamics Simulations

Mutual interactions between components of biological membranes are pivotal for maintaining their proper biophysical properties, such as stability, fluidity, or permeability. The main building blocks of biomembranes are lipids, among which the most important are phospholipids (mainly phosphatidylcholines (PCs)) and sterols (mainly cholesterol). Although there is a plethora of reports on interactions between PCs, as well as between PCs and cholesterol, their molecular mechanism has not yet been fully explained. Therefore, to resolve this issue, we carried out systematic investigations based on the classical Langmuir monolayer technique complemented with molecular dynamics simulations. The studies involved systems containing 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) analogues possessing in the structure one or two polar functional groups similar to those of DPPC. The interactions and rheological properties of binary mixtures of DPPC analogues with 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) and cholesterol were compared with reference systems (DPPC/POPC and DPPC/cholesterol). This pointed to the importance of the ternary amine group in PC/cholesterol interactions, while in PC mixtures, the phosphate group played a key role. In both cases, the esterified glycerol group had an effect on the magnitude of interactions. The obtained results are crucial for establishing structure-property relationships as well as for designing substitutes for natural lipids.

Overlay databank unlocks data-driven analyses of biomolecules for all

Tools based on artificial intelligence (AI) are currently revolutionising many fields, yet their applications are often limited by the lack of suitable training data in programmatically accessible format. Here we propose an effective solution to make data scattered in various locations and formats accessible for data-driven and machine learning applications using the overlay databank format. To demonstrate the practical relevance of such approach, we present the NMRlipids Databank-a community-driven, open-for-all database featuring programmatic access to quality-evaluated atom-resolution molecular dynamics simulations of cellular membranes. Cellular membrane lipid composition is implicated in diseases and controls major biological functions, but membranes are difficult to study experimentally due to their intrinsic disorder and complex phase behaviour. While MD simulations have been useful in understanding membrane systems, they require significant computational resources and often suffer from inaccuracies in model parameters. Here, we demonstrate how programmable interface for flexible implementation of data-driven and machine learning applications, and rapid access to simulation data through a graphical user interface, unlock possibilities beyond current MD simulation and experimental studies to understand cellular membranes. The proposed overlay databank concept can be further applied to other biomolecules, as well as in other fields where similar barriers hinder the AI revolution.

Quantitative Comparison against Experiments Reveals Imperfections in Force Fields' Descriptions of POPC-Cholesterol Interactions

Cholesterol is a central building block in biomembranes, where it induces orientational order, slows diffusion, renders the membrane stiffer, and drives domain formation. Molecular dynamics (MD) simulations have played a crucial role in resolving these effects at the molecular level; yet, it has recently become evident that different MD force fields predict quantitatively different behavior. Although easily neglected, identifying such limitations is increasingly important as the field rapidly progresses toward simulations of complex membranes mimicking the in vivo conditions: pertinent multicomponent simulations must capture accurately the interactions between their fundamental building blocks, such as phospholipids and cholesterol. Here, we define quantitative quality measures for simulations of binary lipid mixtures in membranes against the C-H bond order parameters and lateral diffusion coefficients from NMR spectroscopy as well as the form factors from X-ray scattering. Based on these measures, we perform a systematic evaluation of the ability of commonly used force fields to describe the structure and dynamics of binary mixtures of palmitoyloleoylphosphatidylcholine (POPC) and cholesterol. None of the tested force fields clearly outperforms the others across the tested properties and conditions. Still, the Slipids parameters provide the best overall performance in our tests, especially when dynamic properties are included in the evaluation. The quality evaluation metrics introduced in this work will, particularly, foster future force field development and refinement for multicomponent membranes using automated approaches.

Hydrophobic Matching Dictates over the Linear Rule of Mixtures in Binary Lipid Membranes

Biological membranes feature heterogeneous mixtures of lipids with different head and tail characteristics. Their biophysical properties are dictated by the intimate interaction among different constituent lipids. Previous studies suggest that the membrane area-per-lipid (APL) deviates from the linear rule of mixtures (LRM) for binary lipid membranes, but the underlying mechanism remains elusive. Our molecular dynamics (MD) simulations of binary lipid membranes consisting of lipids with different tail characteristics reveal a competitive mechanism whereby lipids tend to deform each other to minimize the hydrophobic mismatch between their tails. Depending on the relative tail lengths and saturation levels, this may result in an either positive or negative deviation of APL from the LRM. As lipid packing plays an essential role in membrane fusion and peptide-membrane binding, our findings may help guide the selection of lipids for the effective rational design of nanoliposomes and membrane-targeting peptides.

Conditions for the stable adsorption of lipid monolayers to solid surfaces

Lipid monolayers are ubiquitous in biological systems and have multiple roles in biotechnological applications, such as lipid coatings that enhance colloidal stability or prevent surface fouling. Despite the great technological importance of surface-adsorbed lipid monolayers, the connection between their formation and the chemical characteristics of the underlying surfaces has remained poorly understood. Here, we elucidate the conditions required for stable lipid monolayers nonspecifically adsorbed on solid surfaces in aqueous solutions and water/alcohol mixtures. We use a framework that combines the general thermodynamic principles of monolayer adsorption with fully atomistic molecular dynamics simulations. We find that, very universally, the chief descriptor of adsorption free energy is the wetting contact angle of the solvent on the surface. It turns out that monolayers can form and remain thermodynamically stable only on substrates with contact angles above the adsorption contact angle, θads. Our analysis establishes that θads falls into a narrow range of around 60∘-70∘ in aqueous media and is only weakly dependent on the surface chemistry. Moreover, to a good approximation, θads is roughly determined by the ratio between the surface tensions of hydrocarbons and the solvent. Adding small amounts of alcohol to the aqueous medium lowers θads and thereby facilitates monolayer formation on hydrophilic solid surfaces. At the same time, alcohol addition weakens the adsorption strength on hydrophobic surfaces and results in a slowdown of the adsorption kinetics, which can be useful for the preparation of defect-free monolayers.

Structure of POPC Lipid Bilayers in OPLS3e Force Field

It is crucial for molecular dynamics simulations of biomembranes that the force field parameters give a realistic model of the membrane behavior. In this study, we examined the OPLS3e force field for the carbon-hydrogen order parameters S_CH of POPC (1-palmitoyl-2-oleoylphosphatidylcholine) lipid bilayers at varying hydration conditions and ion concentrations. The results show that OPLS3e behaves similarly to the CHARMM36 force field and relatively accurately follows the experimentally measured S_CH for the lipid headgroup, the glycerol backbone, and the acyl tails. Thus, OPLS3e is a good choice for POPC bilayer simulations under many biologically relevant conditions. The exception are systems with an abundancy of ions, as similarly to most other force fields OPLS3e strongly overestimates the membrane-binding of cations, especially Ca²⁺. This leads to undesirable positive charge of the membrane surface and drastically lowers the concentration of Ca²⁺ in the surrounding solvent, which might cause issues in systems sensitive to correct charge distribution profiles across the membrane.

Determining structure and action mechanism of LBF14 by molecular simulation

The recently discovered, membrane-active peptide LBF14 contains several non-proteinogenic amino acids and is able to transform vesicles into tubule networks. The exact membrane interaction mechanism and detailed secondary structure are yet to be determined. We performed molecular dynamics simulations of LBF14 and let it fold de novo into its ensemble of native secondary structures. Histidine protonation state effects on secondary structure were investigated. An MD simulation of the peptide with a lipid bilayer was performed. Simulation results were compared to circular dichroism and electron paramagnetic resonance data of previous studies. LBF14 contains a conserved helical section in an otherwise random structure. Helical stability is influenced by histidine protonation. The peptide localized to the polar layer of the membrane, consistent with experimental results. While the overall secondary structure is unaffected by membrane interaction, Ramachandran plot analysis yielded two distinct peptide conformations during membrane interaction. This conformational change was accompanied by residue repositioning within the membrane. LBF14 only affected the local order in the membrane, and had no measurable effect on pressure. The simulation results are consistent with the previously proposed membrane interaction mechanism of LBF14 and can additionally explain the local interaction mechanism. Communicated by Ramaswamy H. Sarma.

Molecular dynamics simulations and solid-state nuclear magnetic resonance spectroscopy measurements of C-H bond order parameters and effective correlation times in a POPC-GM3 bilayer

Glycolipids such as gangliosides affect the properties of lipid membranes and in extension the interactions between membranes and other biomolecules like proteins. To better understand how the properties of individual lipid molecules can contribute to shape the functional aspects of a membrane, the spatial restriction and dynamics of C-H bond segments can be measured using nuclear magnetic resonance (NMR) spectroscopy. We combine solid-state NMR spectroscopy with all-atom molecular dynamics (MD) simulations to investigate how ganglioside GM3 affects the bilayer structure and dynamics of C-H bond segments. These two methods yield reorientational correlation functions, molecular profiles of C-H bond order parameters |S_CH| and effective correlation times τ_e, which we compare for lipids in POPC bilayers with and without 30 mol% GM3. Our results revealed that all C-H segments of POPC reorient slower in the presence of GM3 and that the defining features of the GM3-POPC bilayer lie in the GM3 headgroup; it gives the bilayer an extended headgroup layer with high order (|S_CH| up to 0.3-0.4) and slow dynamics (τ_e up to 100 ns), a character that may be mechanistically important in ganglioside interactions with other biomolecules.

The Structures of Heterogeneous Membranes and Their Interactions with an Anticancer Peptide: A Molecular Dynamics Study

We conduct molecular dynamics simulations of model heterogeneous membranes and their interactions with a 24-amino acid peptide-NAF-144-67. NAF-144-67 is an anticancer peptide that selectively permeates and kills malignant cells; it does not permeate normal cells. We examine three membranes with different binary mixtures of lipids, DOPC-DOPA, DOPC-DOPS, and DOPC-DOPE, with a single peptide embedded in each as models for the diversity of biological membranes. We illustrate that the peptide organization in the membrane depends on the types of nearby phospholipids and is influenced by the charge and size of the head groups. The present study sheds light on early events of permeation and the mechanisms by which an amphiphilic peptide crosses from an aqueous solution to a hydrophobic membrane. Understanding the translocation mechanism is likely to help the design of new permeants.

Cations Do Not Alter the Membrane Structure of POPC—A Lipid With an Intermediate Area

Combining small-angle neutron scattering (SANS), small-angle X-ray scattering (SAXS), and densitometric measurements, we have studied the interactions of the divalent cations Ca²⁺ and Mg²⁺ with the lipid vesicles prepared of a mixed-chain palmitoyl-oleoyl-phosphatidylcholine (POPC) at 25°C. The structural parameters of the POPC bilayer, such as the bilayer thickness, lateral area, and volume per lipid, displayed no changes upon the ion addition at concentrations up to 30 mM and minor changes at > 30 mM Ca²⁺ and Mg²⁺, while some decrease in the vesicle radius was observed over the entire concentration range studied. This examination allows us to validate the concept of lipid–ion interactions governed by the area per lipid suggested previously and to propose the mixed mode of those interactions that emerge in the POPC vesicles. We speculate that the average area per POPC lipid that corresponds to the cutoff length of lipid–ion interactions generates an equal but opposite impact on ion bridges and separate lipid–ion pairs. As a result of the dynamic equilibrium, the overall structural properties of bilayers are not affected. As the molecular mechanism proposed is affected by the structural properties of a particular lipid, it might help us to understand the fundamentals of processes occurring in complex multicomponent membrane systems.

Molecular Modeling is an Enabling Approach to Complement and Enhance Channelopathy Research

Hundreds of human membrane proteins form channels that transport necessary ions and compounds, including drugs and metabolites, yet details of their normal function or how function is altered by genetic variants to cause diseases are often unknown. Without this knowledge, researchers are less equipped to develop approaches to diagnose and treat channelopathies. High-resolution computational approaches such as molecular modeling enable researchers to investigate channelopathy protein function, facilitate detailed hypothesis generation, and produce data that is difficult to gather experimentally. Molecular modeling can be tailored to each physiologic context that a protein may act within, some of which may currently be difficult or impossible to assay experimentally. Because many genomic variants are observed in channelopathy proteins from high-throughput sequencing studies, methods with mechanistic value are needed to interpret their effects. The eminent field of structural bioinformatics integrates techniques from multiple disciplines including molecular modeling, computational chemistry, biophysics, and biochemistry, to develop mechanistic hypotheses and enhance the information available for understanding function. Molecular modeling and simulation access 3D and time-dependent information, not currently predictable from sequence. Thus, molecular modeling is valuable for increasing the resolution with which the natural function of protein channels can be investigated, and for interpreting how genomic variants alter them to produce physiologic changes that manifest as channelopathies. © 2022 American Physiological Society. Compr Physiol 12:3141-3166, 2022.

Lipid21: Complex Lipid Membrane Simulations with AMBER

We extend the modular AMBER lipid force field to include anionic lipids, polyunsaturated fatty acid (PUFA) lipids, and sphingomyelin, allowing the simulation of realistic cell membrane lipid compositions, including raft-like domains. Head group torsion parameters are revised, resulting in improved agreement with NMR order parameters, and hydrocarbon chain parameters are updated, providing a better match with phase transition temperature. Extensive validation runs (0.9 μs per lipid type) show good agreement with experimental measurements. Furthermore, the simulation of raft-like bilayers demonstrates the perturbing effect of increasing PUFA concentrations on cholesterol molecules. The force field derivation is consistent with the AMBER philosophy, meaning it can be easily mixed with protein, small molecule, nucleic acid, and carbohydrate force fields.

Rotational decoupling between the hydrophilic and hydrophobic regions in lipid membranes

Cells use homeostatic mechanisms to ensure an optimal composition of distinct types of lipids in cellular membranes. The hydrophilic region of biological lipid membranes is mainly composed of several types of phospholipid headgroups that interact with incoming molecules, nanoparticles, and viruses, whereas the hydrophobic region consists of a distribution of acyl chains and sterols affecting membrane fluidity/rigidity related properties and forming an environment for membrane-bound molecules such as transmembrane proteins. A fundamental open question is to what extent the motions of these regions are coupled and, consequently, how strongly the interactions of phospholipid headgroups with other molecules depend on the properties and composition of the membrane hydrophobic core. We combine advanced solid-state nuclear magnetic resonance spectroscopy with high-fidelity molecular dynamics simulations to demonstrate how the rotational dynamics of choline headgroups remain nearly unchanged (slightly faster) with incorporation of cholesterol into a phospholipid membrane, contrasting the well-known extreme slowdown of the other phospholipid segments. Notably, our results suggest a new paradigm in which phospholipid dipole headgroups interact as quasi-freely rotating flexible dipoles at the interface, independent of the properties in the hydrophobic region.

Improved Parameterization of Phosphatidylinositide Lipid Headgroups for the Martini 3 Coarse-Grain Force Field

Phosphoinositides are a family of membrane phospholipids that play crucial roles in membrane regulatory events. As such, these lipids are often a key part of molecular dynamics simulation studies of biological membranes, in particular of those employing coarse-grain models because of the potential long times and sizes of the involved membrane processes. Version 3 of the widely used Martini coarse-grain force field has been recently published, greatly refining many aspects of biomolecular interactions. In order to properly use it for lipid membrane simulations with phosphoinositides, we put forth the Martini 3-specific parameterization of inositol, phosphatidylinositol, and seven physiologically relevant phosphorylated derivatives of phosphatidylinositol. Compared to parameterizations for earlier Martini versions, focus was put on a more accurate reproduction of the behavior seen in both atomistic simulations and experimental studies, including the signaling-relevant phosphoinositide interaction with divalent cations. The models that we develop improve upon the conformational dynamics of phosphoinositides in the Martini force field and provide stable topologies at typical Martini time steps. They are able to reproduce experimentally known protein-binding poses as well as phosphoinositide aggregation tendencies. The latter was tested both in the presence and absence of calcium and included correct behavior of PI(4,5)P₂ calcium-induced clusters, which can be of relevance for regulation.

Mechanistic Understanding from Molecular Dynamics in Pharmaceutical Research 2: Lipid Membrane in Drug Design

We review the use of molecular dynamics (MD) simulation as a drug design tool in the context of the role that the lipid membrane can play in drug action, i.e., the interaction between candidate drug molecules and lipid membranes. In the standard "lock and key" paradigm, only the interaction between the drug and a specific active site of a specific protein is considered; the environment in which the drug acts is, from a biophysical perspective, far more complex than this. The possible mechanisms though which a drug can be designed to tinker with physiological processes are significantly broader than merely fitting to a single active site of a single protein. In this paper, we focus on the role of the lipid membrane, arguably the most important element outside the proteins themselves, as a case study. We discuss work that has been carried out, using MD simulation, concerning the transfection of drugs through membranes that act as biological barriers in the path of the drugs, the behavior of drug molecules within membranes, how their collective behavior can affect the structure and properties of the membrane and, finally, the role lipid membranes, to which the vast majority of drug target proteins are associated, can play in mediating the interaction between drug and target protein. This review paper is the second in a two-part series covering MD simulation as a tool in pharmaceutical research; both are designed as pedagogical review papers aimed at both pharmaceutical scientists interested in exploring how the tool of MD simulation can be applied to their research and computational scientists interested in exploring the possibility of a pharmaceutical context for their research.

Inverse Conformational Selection in Lipid-Protein Binding

Interest in lipid interactions with proteins and other biomolecules is emerging not only in fundamental biochemistry but also in the field of nanobiotechnology where lipids are commonly used, for example, in carriers of mRNA vaccines. The outward-facing components of cellular membranes and lipid nanoparticles, the lipid headgroups, regulate membrane interactions with approaching substances, such as proteins, drugs, RNA, or viruses. Because lipid headgroup conformational ensembles have not been experimentally determined in physiologically relevant conditions, an essential question about their interactions with other biomolecules remains unanswered: Do headgroups exchange between a few rigid structures, or fluctuate freely across a practically continuous spectrum of conformations? Here, we combine solid-state NMR experiments and molecular dynamics simulations from the NMRlipids Project to resolve the conformational ensembles of headgroups of four key lipid types in various biologically relevant conditions. We find that lipid headgroups sample a wide range of overlapping conformations in both neutral and charged cellular membranes, and that differences in the headgroup chemistry manifest only in probability distributions of conformations. Furthermore, the analysis of 894 protein-bound lipid structures from the Protein Data Bank suggests that lipids can bind to proteins in a wide range of conformations, which are not limited by the headgroup chemistry. We propose that lipids can select a suitable headgroup conformation from the wide range available to them to fit the various binding sites in proteins. The proposed inverse conformational selection model will extend also to lipid binding to targets other than proteins, such as drugs, RNA, and viruses.

A Data-Driven Dimensionality Reduction Approach to Compare and Classify Lipid Force Fields

Molecular dynamics simulations of all-atom and coarse-grained lipid bilayer models are increasingly used to obtain useful insights for understanding the structural dynamics of these assemblies. In this context, one crucial point concerns the comparison of the performance and accuracy of classical force fields (FFs), which sometimes remains elusive. To date, the assessments performed on different classical potentials are mostly based on the comparison with experimental observables, which typically regard average properties. However, local differences of the structure and dynamics, which are poorly captured by average measurements, can make a difference, but these are nontrivial to catch. Here, we propose an agnostic way to compare different FFs at different resolutions (atomistic, united-atom, and coarse-grained), by means of a high-dimensional similarity metrics built on the framework of Smooth Overlap of Atomic Position (SOAP). We compare and classify a set of 13 FFs, modeling 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) bilayers. Our SOAP kernel-based metrics allows us to compare, discriminate, and correlate different FFs at different model resolutions in an unbiased, high-dimensional way. This also captures differences between FFs in modeling nonaverage events (originating from local transitions), for example, the liquid-to-gel phase transition in dipalmitoylphosphatidylcholine (DPPC) bilayers, for which our metrics allows us to identify nucleation centers for the phase transition, highlighting some intrinsic resolution limitations in implicit versus explicit solvent FFs.

Effect of magnesium sulfate in oxidized lipid bilayers properties by using molecular dynamics

Magnesium sulfate (MgSO₄) has been used as a protector agent for many diseases related to oxidative stress. The effect of MgSO₄ on the oxidized lipid bilayer has not yet been studied using molecular dynamics calculations. In this work, the effects of oxidation were evaluated by using a POPC membrane model at different concentrations of its aldehyde (-CHO) and hydroperoxide (-OOH) derivatives with and without MgSO₄. Several quantitative and qualitative properties were evaluated, such as membrane thickness, area per lipid, area compressibility modulus, snapshots after simulation finish, density distributions, time evolutions of oxidized group positions, and radial distributions of oxidized group concerning Mg. Results indicate that in the absence of MgSO₄ the mobility of oxidized groups, particularly –CHO, toward the surface interface is high. At a low oxidation level of the bilayer there is an increase in the compressibility modulus as compared to the unoxidized bilayer. MgSO₄, at a low oxidation level, tends to lessen the oxidation effects by lowering the dispersion in the distribution of oxidized species toward the membrane surface and the water region. However, MgSO₄ does not change the trends of decreasing membrane thickness and area compressibility modulus and increasing area per lipid upon oxidation. In this regard, MgSO₄ diminishes the electrostatic long-distance attractive interactions between the oxidized groups and the charged headgroups of the interface, owing to the Mg⁺² and SO₄^-2 screening effects and an electrostatic stabilization of the headgroups, preventing the pore formation, which is well-known to occur in oxidized membranes.

•

MgSO₄ in vitro restores oxidized membranes but its molecular mechanism is unclear.

•

MD simulations of oxidized lipid bilayers were performed with and without of MgSO₄.

•

A restriction in the mobility of oxidized groups is produced by MgSO₄.

•

Mg⁺² and SO₄⁼ produce screening effects on the oxidized membranes.

•

MgSO₄ produce a diminution of electrostatic long-distance attractive interactions.

Molecular Study of Ultrasound-Triggered Release of Fluorescein from Liposomes

Several investigations have suggested that ultrasound triggers the release of drugs encapsulated into liposomes at acoustic pressures low enough to avoid cavitation or high hyperthermia. However, the mechanism leading to this triggered release as well as the adequate composition of the liposome membrane remains unknown. Here, we investigate the ultrasound-triggered release of fluorescein disodium salt encapsulated into liposomes made of 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) or 1,2-distearoylphosphatidyl-ethanolamine (DSPC) lipids with various concentrations of cholesterol (from 0 to 44 mol %). The passive release of encapsulated fluorescein was first characterized. It was observed to be higher when the membrane is in a fluid phase and increased with temperature but decreased upon addition of cholesterol. Next, the release of fluorescein was measured at different acoustic frequencies (0.8, 1.1, and 3.3 MHz) and peak-to-peak pressures (0, 2, 2.5, 5, and 8 MPa). Measurements were performed at temperatures where DOPC and DSPC liposomes were, respectively, in the fluid or gel phase. We found that the release rate did not depend on the ultrasound frequency. For DOPC liposomes, the ultrasound-triggered release of fluorescein decreased with increasing concentration of cholesterol in liposomes, while the behavior was more complex for DSPC liposomes. Overall, the triggered release from DSPC liposomes was up to ten times less than DOPC liposomes. Molecular dynamics simulations performed on a pure DOPC membrane showed that a membrane experiences, under a directional pressure of ±2.4 MPa, various changes in properties such as the area per lipid (APL). An increase in the APL was notably observed when the simulation box was laterally stretched or perpendicularly compressed, which was accompanied by an increase in the number of water molecules crossing the membrane. This suggests that ultrasound most probably enhances the diffusion of encapsulated molecules at small acoustic pressures by increasing the distance between lipids.

A narrative review of urinary phospholipids: from biochemical aspect towards clinical application

As a newly emerged discipline, lipidomic studies have focused on the comprehensive characterization and quantification of lipids in a given biological system, which has remarkably advanced in recent years owing to the rapid development of analytical techniques, especially mass spectrometry. Among diverse lipid classes, phospholipids, which have fundamental roles in the formation of cellular membranes, signaling processes, and bioenergetics have gained momentum in several fields of research. The altered composition, concentration, spatial distribution, and metabolism of phospholipids in cells, tissues, and body fluids have been elucidated in various human diseases such as cancer, inflammation, as well as cardiovascular and metabolic disorders. Among the different kinds of phospholipid sources in the human body, urine has not been extensively investigated in recent years owing to the extremely low concentrations of phospholipids and high levels of salts and other contaminants, which can interfere with precise detection. However, with profound advances and rapid expansion in analytical methods, urinary phospholipids have attracted increasing attention in current biomedical research as urine is an easily available source for the discovery of noninvasive biomarkers. In this review, we provide an overview of urinary phospholipids, including their biochemical aspects and clinical applications, aimed at promoting this field of research.

CHARMM36 Lipid Force Field with Explicit Treatment of Long-Range Dispersion: Parametrization and Validation for Phosphatidylethanolamine, Phosphatidylglycerol, and Ether Lipids

Long-range Lennard-Jones (LJ) interactions have been incorporated into the CHARMM36 (C36) lipid force field (FF) using the LJ particle-mesh Ewald (LJ-PME) method in order to remove the inconsistency of bilayer and monolayer properties arising from the exclusion of long-range dispersion [Yu, Y.; Semi-automated Optimization of the CHARMM36 Lipid Force Field to Include Explicit Treatment of Long-Range Dispersion. J. Chem. Theory Comput. 2021, 10.1021/acs.jctc.0c01326. (preceding article in this issue)]. The new FF is denoted C36/LJ-PME. While the first optimization was based on three phosphatidylcholines (PCs), this work extends the validation and parametrization to more lipids including PC, phosphatidylethanolamine (PE), phosphatidylglycerol (PG), and ether lipids. The agreement with experimental structure data is excellent for PC, PE, and ether lipids. C36/LJ-PME also compares favorably with scattering data of PG bilayers but less so with NMR deuterium order parameters of 1,2-dimyristoyl-sn-glycero-3-phospho-(1'-rac-glycerol) (DMPG) at 303.15 K, indicating a need for future optimization regarding PG-specific parameters. Frequency dependence of NMR T1 spin-lattice relaxation times is well-described by C36/LJ-PME, and the overall agreement with experiment is comparable to C36. Lipid diffusion is slower than C36 due to the added long-range dispersion causing a higher viscosity, although it is still too fast compared to experiment after correction for periodic boundary conditions. When using a 10 Å real-space cutoff, the simulation speed of C36/LJ-PME is roughly equal to C36. While more lipids will be incorporated into the FF in the future, C36/LJ-PME can be readily used for common lipids and extends the capability of the CHARMM FF by supporting monolayers and eliminating the cutoff dependence.

Semi-automated Optimization of the CHARMM36 Lipid Force Field to Include Explicit Treatment of Long-Range Dispersion

The development of the CHARMM lipid force field (FF) can be traced back to the early 1990s with its current version denoted CHARMM36 (C36). The parametrization of C36 utilized high-level quantum mechanical data and free energy calculations of model compounds before parameters were manually adjusted to yield agreement with experimental properties of lipid bilayers. While such manual fine-tuning of FF parameters is based on intuition and trial-and-error, automated methods can identify beneficial modifications of the parameters via their sensitivities and thereby guide the optimization process. This work introduces a semi-automated approach to reparametrize the CHARMM lipid FF with consistent inclusion of long-range dispersion through the Lennard-Jones particle-mesh Ewald (LJ-PME) approach. The optimization method is based on thermodynamic reweighting with regularization with respect to the C36 set. Two independent optimizations with different topology restrictions are presented. Targets of the optimizations are primarily liquid crystalline phase properties of lipid bilayers and the compression isotherm of monolayers. Pair correlation functions between water and lipid functional groups in aqueous solution are also included to address headgroup hydration. While the physics of the reweighting strategy itself is well-understood, applying it to heterogeneous, complex anisotropic systems poses additional challenges. These were overcome through careful selection of target properties and reweighting settings allowing for the successful incorporation of the explicit treatment of long-range dispersion, and we denote the newly optimized lipid force field as C36/LJ-PME. The current implementation of the optimization protocol will facilitate the future development of the CHARMM and related lipid force fields.

Ferreira T, Ollila OHS, Miettinen MS. Using Open Data to Rapidly Benchmark Biomolecular Simulations: Phospholipid Conformational Dynamics

Molecular dynamics (MD) simulations are widely used to monitor time-resolved motions of biomacromolecules, although it often remains unknown how closely the conformational dynamics correspond to those occurring in real life. Here, we used a large set of open-access MD trajectories of phosphatidylcholine (PC) lipid bilayers to benchmark the conformational dynamics in several contemporary MD models (force fields) against nuclear magnetic resonance (NMR) data available in the literature: effective correlation times and spin-lattice relaxation rates. We found none of the tested MD models to fully reproduce the conformational dynamics. That said, the dynamics in CHARMM36 and Slipids are more realistic than in the Amber Lipid14, OPLS-based MacRog, and GROMOS-based Berger force fields, whose sampling of the glycerol backbone conformations is too slow. The performance of CHARMM36 persists when cholesterol is added to the bilayer, and when the hydration level is reduced. However, for conformational dynamics of the PC headgroup, both with and without cholesterol, Slipids provides the most realistic description because CHARMM36 overestimates the relative weight of ∼1 ns processes in the headgroup dynamics. We stress that not a single new simulation was run for the present work. This demonstrates the worth of open-access MD trajectory databanks for the indispensable step of any serious MD study: benchmarking the available force fields. We believe this proof of principle will inspire other novel applications of MD trajectory databanks and thus aid in developing biomolecular MD simulations into a true computational microscope-not only for lipid membranes but for all biomacromolecular systems.

High-performance macromolecular data delivery and visualization for the web

Biomacromolecular structural data make up a vital and crucial scientific resource that has grown not only in terms of its amount but also in its size and complexity. Furthermore, these data are accompanied by large and increasing amounts of experimental data. Additionally, the macromolecular data are enriched with value-added annotations describing their biological, physicochemical and structural properties. Today, the scientific community requires fast and fully interactive web visualization to exploit this complex structural information. This article provides a survey of the available cutting-edge web services that address this challenge. Specifically, it focuses on data-delivery problems, discusses the visualization of a single structure, including experimental data and annotations, and concludes with a focus on the results of molecular-dynamics simulations and the visualization of structural ensembles.

Mechanistic Understanding From Molecular Dynamics Simulation in Pharmaceutical Research 1: Drug Delivery

In this review, we outline the growing role that molecular dynamics simulation is able to play as a design tool in drug delivery. We cover both the pharmaceutical and computational backgrounds, in a pedagogical fashion, as this review is designed to be equally accessible to pharmaceutical researchers interested in what this new computational tool is capable of and experts in molecular modeling who wish to pursue pharmaceutical applications as a context for their research. The field has become too broad for us to concisely describe all work that has been carried out; many comprehensive reviews on subtopics of this area are cited. We discuss the insight molecular dynamics modeling has provided in dissolution and solubility, however, the majority of the discussion is focused on nanomedicine: the development of nanoscale drug delivery vehicles. Here we focus on three areas where molecular dynamics modeling has had a particularly strong impact: (1) behavior in the bloodstream and protective polymer corona, (2) Drug loading and controlled release, and (3) Nanoparticle interaction with both model and biological membranes. We conclude with some thoughts on the role that molecular dynamics simulation can grow to play in the development of new drug delivery systems.

Outer Membranes of Polymyxin-Resistant Acinetobacter baumannii with Phosphoethanolamine-Modified Lipid A and Lipopolysaccharide Loss Display Different Atomic-Scale Interactions with Polymyxins

Resistance to the last-line polymyxins is increasingly reported in multidrug-resistant Gram-negative pathogens, including Acinetobacter baumannii, which develops resistance via either lipid A modification (e.g., with phosphoethanolamine [pEtN]) or even lipopolysaccharide (LPS) loss in the outer membrane (OM). Considering these two different mechanisms, quantitative membrane lipidomics data were utilized to develop three OM models representing polymyxin-susceptible and -resistant A. baumannii strains. Through all-atom molecular simulations with enhanced sampling techniques, the effect of lipid A-pEtN modification and LPS loss on the action of colistin (i.e., polymyxin E) was examined for the first time, with a focus on the dynamics and energetics of colistin penetration into these OMs. Lipid A-pEtN modification improved the OM stability, impeding the penetration of colistin into the OM; this differed from the current literature that lipid A-pEtN modification confers resistance by diminishing the initial interaction with polymyxins. In contrast, the LPS deficiency significantly reduced the negative charges on the OM surface, diminishing the binding of colistin. Moreover, both lipid A-pEtN modification and LPS loss also constituted colistin resistance through disturbing the conformational transitions of the colistin molecule. Collectively, atomic-scale interactions between polymyxins and different bacterial OMs are very different and the findings may facilitate the discovery of new-generation polymyxins against Gram-negative 'superbugs'.

Optimization of Slipids Force Field Parameters Describing Headgroups of Phospholipids

The molecular mechanics force field Slipids developed in a series of works by Jämbeck and Lyubartsev (J. Phys. Chem. B 2012, 116, 3164-3179; J. Chem. Theory Comput. 2012, 8, 2938-2948) generally provides a good description of various lipid bilayer systems. However, it was also found that order parameters of C-H bonds in the glycerol moiety of the phosphatidylcholine headgroup deviate significantly from NMR results. In this work, the dihedral force field parameters have been reparameterized in order to improve the agreement with experiment. For this purpose, we have computed energies for a large amount of lipid headgroup conformations using density functional theory on the B3P86/cc-pvqz level and optimized dihedral angle parameters simultaneously to provide the best fit to the quantum chemical energies. The new parameter set was validated for three lipid bilayer systems against a number of experimental properties including order parameters, area per lipid, scattering form factors, bilayer thickness, area compressibility and lateral diffusion coefficients. In addition, the order parameter dependence on cholesterol content in the POPC bilayer was investigated. It is shown that the new force field significantly improves agreement with the experimental order parameters for the lipid headgroup while keeping good agreement with other experimentally measured properties.

How Do Molecular Dynamics Data Complement Static Structural Data of GPCRs

G protein-coupled receptors (GPCRs) are implicated in nearly every physiological process in the human body and therefore represent an important drug targeting class. Advances in X-ray crystallography and cryo-electron microscopy (cryo-EM) have provided multiple static structures of GPCRs in complex with various signaling partners. However, GPCR functionality is largely determined by their flexibility and ability to transition between distinct structural conformations. Due to this dynamic nature, a static snapshot does not fully explain the complexity of GPCR signal transduction. Molecular dynamics (MD) simulations offer the opportunity to simulate the structural motions of biological processes at atomic resolution. Thus, this technique can incorporate the missing information on protein flexibility into experimentally solved structures. Here, we review the contribution of MD simulations to complement static structural data and to improve our understanding of GPCR physiology and pharmacology, as well as the challenges that still need to be overcome to reach the full potential of this technique.

Interfacial Water Structure at Zwitterionic Membrane/Water Interface: The Importance of Interactions between Water and Lipid Carbonyl Groups

In this work, atomistic molecular dynamics (MD) simulations of palmitoyl-oleoyl-phosphatidylcholine (POPC) bilayer were carried out to investigate the effect of water models on membrane dipole potential, which is primarily associated with the preferential orientation of molecular dipoles at the membrane-water interface. We discovered that the overestimation of the dipole potential by the TIPS3P water model can be effectively reduced by the TIP4P water model. On the one hand, the TIP4P water model decreases the negative contribution of lipid to the dipole potential through influencing the orientation of lipid headgroups. On the other hand, the TIP4P water model reduces the positive contribution of water to the dipole potential by increasing the preference of H-down orientation (the water dipole orients toward the bilayer center). Interestingly, the TIP4P water model affects the orientation of interfacial water molecules more obviously than that of lipid headgroups, leading to the decrease in the dipole potential. Furthermore, the MD results revealed that the water close to the positively charged choline (namely, N-associated water) prefers the H-down orientation while the water around the negatively charged phosphate (namely, P-associated water) favors the H-up orientation, in support of recent experimental and MD studies. However, interfacial water molecules are more strongly influenced by the phosphate groups than by the choline groups, resulting in the net H-up orientation (the water dipole orients toward the bilayer center) in the region of lipid headgroups. In addition, it is intriguing that the preference of H-up orientation decreases when water molecules penetrate more deeply into the lipid bilayer. This is attributed to the counteracting effect of lipid carbonyl groups, and the effect varies with the lipid chains (oleoyl and palmitoyl chains), suggesting the important role of lipid carbonyl groups.

Principal component analysis highlights the influence of temperature, curvature and cholesterol on conformational dynamics of lipids

Membrane lipids are inherently highly dynamic molecules. Currently, it is difficult to probe the structures of individual lipids experimentally at the timescales corresponding to atomic motions, and consequently molecular dynamics simulations are used widely. In our previous work, we have introduced the principal component analysis (PCA) as a convenient framework for comprehensive quantitative description of lipid motions. Here, we present a newly developed open source script, PCAlipids, which automates the analysis and allows us to refine the approach and test its limitations. We use PCAlipids to determine the influence of temperature, cholesterol and curvature on individual lipids, and show that the most prominent lipid tail scissoring motion is strongly affected by these factors and allows tracking of phase transition. Addition of cholesterol affects the conformations and selectively changes the dynamics of lipid molecules, impacting the large-amplitude motions. Introduction of curvature biases the conformational ensembles towards more extended structures. We hope that the developed approach will be useful for understanding the molecular basis of different processes occurring in lipid membrane systems and will stimulate development of complementary experimental techniques probing the conformations of individual lipid molecules.

Dual Resolution Membrane Simulations Using Virtual Sites

All-atomistic (AA) and coarse-grain (CG) simulations have been successfully applied to investigate a broad range of biomolecular processes. However, the accessible time and length scales of AA simulation are limited and the specific molecular details of CG simulation are simplified. Here, we propose a virtual site (VS) based hybrid scheme that can concurrently couple AA and CG resolutions in a single membrane simulation, mitigating the shortcomings of either representation. With some adjustments to make the AA and CG force fields compatible, we demonstrate that lipid bilayer properties are well kept in our hybrid approach. Our VS hybrid method was also applied to simulate a small lipid vesicle, with the inner leaflet and interior solvent represented in AA, and the outer leaflet together with exterior solvent at the CG level. Our multiscale method opens the way to investigate biomembrane properties at increased computational efficiency, in particular applications involving large solvent filled regions.

Improved Cation Binding to Lipid Bilayers with Negatively Charged POPS by Effective Inclusion of Electronic Polarization

Phosphatidylserine (PS) lipids are important signaling molecules and the most common negatively charged lipids in eukaryotic membranes. The signaling can be often regulated by calcium, but its interactions with PS headgroups are not fully understood. Classical molecular dynamics (MD) simulations can potentially give detailed description of lipid-ion interactions, but the results strongly depend on the used force field. Here, we apply the electronic continuum correction (ECC) to the Amber Lipid17 parameters of 1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-l-serine (POPS) lipid to improve its interactions with K⁺, Na⁺, and Ca²⁺ ions. The partial charges of the headgroup, glycerol backbone, and carbonyls of POPS, bearing a unit negative charge, were scaled with a factor of 0.75, derived for monovalent ions, and the Lennard-Jones σ parameters of the same segments were scaled with a factor of 0.89. The resulting ECC-POPS model gives more realistic interactions with Na⁺ and Ca²⁺ cations than the original Amber Lipid17 parameters when validated using headgroup order parameters and the "electrometer concept". In ECC-lipids simulations, populations of complexes of Ca²⁺ cations with more than two PS lipids are negligible, and interactions of Ca²⁺ cations with only carboxylate groups are twice more likely than with only phosphate groups, while interactions with carbonyls almost entirely involve other groups as well. Our results pave the way for more realistic MD simulations of biomolecular systems with anionic membranes, allowing signaling processes involving PS and Ca²⁺ to be elucidated.

Molecular insights into the effect of an apoptotic raft-like bilayer on the conformation and dynamics of calreticulin

Cell surface calreticulin (CRT) can mediate apoptotic cells removal by binding and activating LDL receptor-related protein (LRP1). Phosphatidylserine (PS) lipids in the inner leaflet of the cell membrane are externalized and become exposed in cholesterol (CHOL)-rich membrane raft-like microdomain during apoptosis and co-localized with cell surface CRT. How the apoptotic raft-like membrane microdomain affects the structure and dynamics of CRT, further affecting CRT binding with LRP1 to signal apoptotic-cell clearance, remains unknown. In this study, we investigate the interactions between CRT and raft-like bilayers with or without POPS lipids with molecular dynamics simulations. In addition, the effect of an apoptotic raft bilayer on the binding between CRT and thrombospondin-1 (TSP1), a ligand of CRT on the cell surface to signal focal adhesion disassembly, was also investigated. Results of single CRT interactions with raft-like bilayers show that PS lipids in apoptotic raft-like bilayer increased the interactions between CRT and lipid bilayer, which enhanced the conformational stability and increased dynamical motion of CRT. The microscopic and mesoscopic properties of apoptotic raft-like bilayer were altered by the binding of CRT with lipid bilayer. Results of CRT-TSP1 complex interactions with raft-like bilayers show that the binding free energy between TSP1 and CRT was reduced in apoptotic raft-like bilayer environment. This study provides molecular and structural insight into the effect of an apoptotic raft-like bilayer on the conformation and dynamics of CRT, which could enrich our understanding of CRT-mediated apoptotic-cell clearance and focal adhesion disassembly.

A molecular dynamics study proposing the existence of statistical structural heterogeneity due to chain orientation in the POPC-cholesterol bilayer

We performed molecular dynamics simulations of a lipid bilayer consisting of POPC and cholesterol at temperatures from 283 to 308K and cholesterol concentrations from 0 to 50% mol/mol. The purpose of this study was to look for the existence of structural differences in the region delimited by these parameters and, in particular, in a region where coexistence of liquid disordered and liquid ordered phases has been proposed. Our interest in this range of concentration and temperature responds to the fact that polyene ionophore activity varies considerably along it. Two force fields, CHARMM36 and Slipids, were compared in order to determine the most suitable. Both force fields predict non-monotonic behaviors consistent with the existence of phase transitions. We found the presence of lateral structural heterogeneity, statistical in nature, in some of the bilayers occurring in this range of temperatures and sterol concentrations. This heterogeneity was produced by correlated ordering of the POPC tails and not due to cholesterol enrichment, and lasts for tens of nanoseconds. We relate these observations to the action of polyenes in these membranes.

Sharing Data from Molecular Simulations

Given the need for modern researchers to produce open, reproducible scientific output, the lack of standards and best practices for sharing data and workflows used to produce and analyze molecular dynamics (MD) simulations has become an important issue in the field. There are now multiple well-established packages to perform molecular dynamics simulations, often highly tuned for exploiting specific classes of hardware, each with strong communities surrounding them, but with very limited interoperability/transferability options. Thus, the choice of the software package often dictates the workflow for both simulation production and analysis. The level of detail in documenting the workflows and analysis code varies greatly in published work, hindering reproducibility of the reported results and the ability for other researchers to build on these studies. An increasing number of researchers are motivated to make their data available, but many challenges remain in order to effectively share and reuse simulation data. To discuss these and other issues related to best practices in the field in general, we organized a workshop in November 2018 ( https://bioexcel.eu/events/workshop-on-sharing-data-from-molecular-simulations/ ). Here, we present a brief overview of this workshop and topics discussed. We hope this effort will spark further conversation in the MD community to pave the way toward more open, interoperable, and reproducible outputs coming from research studies using MD simulations.

Lipid Head Group Parameterization for GROMOS 54A8: A Consistent Approach with Protein Force Field Description

Membranes are a crucial component of both bacterial and mammalian cells, being involved in signaling, transport, and compartmentalization. This versatility requires a variety of lipid species to tailor the membrane’s behavior as needed, increasing the complexity of the system. Molecular dynamics simulations have been successfully applied to study model membranes and their interactions with proteins, elucidating some crucial mechanisms at the atomistic detail and thus complementing experimental techniques. An accurate description of the functional interplay of the diverse membrane components crucially depends on the selected parameters that define the adopted force field. A coherent parameterization for lipids and proteins is therefore needed. In this work, we propose and validate new lipid head group parameters for the GROMOS 54A8 force field, making use of recently published parametrizations for key chemical moieties present in lipids. We make use additionally of a new canonical set of partial charges for lipids, chosen to be consistent with the parameterization of soluble molecules such as proteins. We test the derived parameters on five phosphocholine model bilayers, composed of lipid patches four times larger than the ones used in previous studies, and run 500 ns long simulations of each system. Reproduction of experimental data like area per lipid and deuterium order parameters is good and comparable with previous parameterizations, as well as the description of liquid crystal to gel-phase transition. On the other hand, the orientational behavior of the head groups is more realistic for this new parameter set, and this can be crucial in the description of interactions with other polar molecules. For that reason, we tested the interaction of the antimicrobial peptide lactoferricin with two model membranes showing that the new parameters lead to a weaker peptide–membrane binding and give a more realistic outcome in comparing binding to antimicrobial versus mammal membranes.

Accelerating Membrane Simulations with Hydrogen Mass Repartitioning

The time step of atomistic molecular dynamics (MD) simulations is determined by the fastest motions in the system and is typically limited to 2 fs. An increasingly popular approach is to increase the mass of the hydrogen atoms to ∼3 amu and decrease the mass of the parent atom by an equivalent amount. This approach, known as hydrogen-mass repartitioning (HMR), permits time steps up to 4 fs with reasonable simulation stability. While HMR has been applied in many published studies to date, it has not been extensively tested for membrane-containing systems. Here, we compare the results of simulations of a variety of membranes and membrane-protein systems run using a 2 fs time step and a 4 fs time step with HMR. For pure membrane systems, we find almost no difference in structural properties, such as area-per-lipid, electron density profiles, and order parameters, although there are differences in kinetic properties such as the diffusion constant. Conductance through a porin in an applied field, partitioning of a small peptide, hydrogen-bond dynamics, and membrane mixing show very little dependence on HMR and the time step. We also tested a 9 Å cutoff as compared to the standard CHARMM cutoff of 12 Å, finding significant deviations in many properties tested. We conclude that HMR is a valid approach for membrane systems, but a 9 Å cutoff is not.

Cholesterol-Induced Conformational Change in the Sphingomyelin Headgroup

Sphingomyelin (SM) and cholesterol (Cho) are the important lipids for the formation of biologically functional membrane domains, lipid rafts. However, the interaction between Cho and the headgroup of SM remains unclear. In this study, we performed solid-state NMR experiments to reveal the Cho effects on the headgroup conformation using ²H-labeled stearoyl-SM (SSM). Deuterated SSMs at the Cα, Cβ, and Cγ positions of a choline moiety were separately prepared and subjected to NMR measurements to determine the quadrupolar splitting of ²H signals in hydrated SSM unitary and SSM/Cho (1:1) bilayers. Using ²H NMR and ¹³C-³¹P REDOR data, the conformation and orientation of the choline moiety were deduced and compared with those derived from molecular dynamics simulations. In SSM unitary bilayers, three torsional angles in the phosphocholine moiety, P-O-Cα-Cβ, were found to be consecutive +gauche(g)/+g/+g or -g/-g/-g. The orientation and conformation of the SSM headgroup were consistent with the results of our molecular dynamics simulations and the previous results on phosphatidylcholines. The quadrupolar coupling at the α methylene group slightly increased in the presence of Cho, and those at the Cβ and Cγ decreased more significantly, thus suggesting that Cho reduced the gauche conformation at the Cα-Cβ torsion. The conformational ensemble in the presence of Cho may enhance the so-called umbrella effect of the SSM headgroup, resulting in the stabilization of Cho near the SM molecules by concealing the hydrophobic Cho core from interfacial water. We also examined the effect of the chiral centers at the sphingosine chain to the headgroup conformation by determining the enantiomeric excess between the diastereomeric +g/+g/+g and -g/-g/-g conformers using (S)-Cα-deuterated and (R)-Cα-deuterated SSMs. Their ²H NMR measurements showed that the chiral centers induced the slight diastereomeric excess in the SM headgroup conformation.

Computational Modeling of Realistic Cell Membranes

Cell membranes contain a large variety of lipid types and are crowded with proteins, endowing them with the plasticity needed to fulfill their key roles in cell functioning. The compositional complexity of cellular membranes gives rise to a heterogeneous lateral organization, which is still poorly understood. Computational models, in particular molecular dynamics simulations and related techniques, have provided important insight into the organizational principles of cell membranes over the past decades. Now, we are witnessing a transition from simulations of simpler membrane models to multicomponent systems, culminating in realistic models of an increasing variety of cell types and organelles. Here, we review the state of the art in the field of realistic membrane simulations and discuss the current limitations and challenges ahead.

Interaction of Bile Salts With Lipid Bilayers: An Atomistic Molecular Dynamics Study

Bile salts (BS) are biosurfactants crucial for emulsification and intestinal absorption of cholesterol and other hydrophobic compounds such as vitamins and fatty acids. Interaction of BS with lipid bilayers is important for understanding their effects on membranes properties. The latter have relevance in passive diffusion processes through intestinal epithelium such as reabsorption of BS, as well as their degree of toxicity to intestinal flora and their potential applications in drug delivery. In this work, we used molecular dynamics simulations to address at the atomic scale the interactions of cholate, deoxycholate, and chenodeoxycholate, as well as their glycine conjugates with POPC bilayers. In this set of BS, variation of three structural aspects was addressed, namely conjugation with glycine, number and position of hydroxyl substituents, and ionization state. From atomistic simulations, the location and orientation of BS inside the bilayer, and their specific interactions with water and host lipid, such as hydrogen bonding and ion-pair formation, were studied in detail. Membrane properties were also investigated to obtain information on the degree of perturbation induced by the different BS. The results are described and related to a recent experimental study (Coreta-Gomes et al., 2015). Differences in macroscopic membrane partition thermodynamics and translocation kinetics are rationalized in terms of the distinct structures and atomic-scale behavior of the bile salt species. In particular, the faster translocation of cholate is explained by its higher degree of local membrane perturbation. On the other hand, the relatively high partition of the polar glycine conjugates is related to the longer and more flexible side chain, which allows simultaneous efficient solvation of the ionized carboxylate and deep insertion of the ring system.

Free energy analysis of membrane pore formation process in the presence of multiple melittin peptides

Understanding the molecular mechanism underlying pore formation in lipid membranes by antimicrobial peptides is of great importance in biological sciences as well as in drug design applications. Melittin has been widely studied as a pore forming peptide, though the molecular mechanism for pore formation is still illusive. We examined the free energy barrier for the creation of a pore in lipid membranes with and without multiple melittin peptides. It was found that six melittin peptides significantly stabilized a pore, though a small barrier (a few k_BT) for the formation still existed. With five melittin peptides or fewer, the pore formation barrier was much higher, though the established pore was in a local energy minimum. Although seven melittins effectively reduced the free energy barrier, a single melittin peptide left the pore after a long time MD simulation probably because of the overcrowded environment around the bilayer pore. Thus, it is highly selective for the number of melittin peptides to stabilize the membrane pore, as was also suggested by the line tension evaluations. The free energy cost required to insert a single melittin into the membrane is too high to explain the one-by-one insertion mechanism for pore formation, which also supports the collective melittin mechanism for pore formation.

Multiscale Simulations of Biological Membranes: The Challenge To Understand Biological Phenomena in a Living Substance

Biological membranes are tricky to investigate. They are complex in terms of molecular composition and structure, functional over a wide range of time scales, and characterized by nonequilibrium conditions. Because of all of these features, simulations are a great technique to study biomembrane behavior. A significant part of the functional processes in biological membranes takes place at the molecular level; thus computer simulations are the method of choice to explore how their properties emerge from specific molecular features and how the interplay among the numerous molecules gives rise to function over spatial and time scales larger than the molecular ones. In this review, we focus on this broad theme. We discuss the current state-of-the-art of biomembrane simulations that, until now, have largely focused on a rather narrow picture of the complexity of the membranes. Given this, we also discuss the challenges that we should unravel in the foreseeable future. Numerous features such as the actin-cytoskeleton network, the glycocalyx network, and nonequilibrium transport under ATP-driven conditions have so far received very little attention; however, the potential of simulations to solve them would be exceptionally high. A major milestone for this research would be that one day we could say that computer simulations genuinely research biological membranes, not just lipid bilayers.

Influence of Bilayer Size and Number in Multi-Bilayer DOPC Simulations at Full and Low Hydration

Biophysical studies of model cell membranes at full and low hydration are usually carried out using scattering measurements on multi-bilayer systems. Molecular simulations of lipid bilayers aimed at reproducing those experimental conditions are usually conducted using single bilayers with different amounts of water. These simulation conditions may lead to artifacts arising from size effects and self-interactions because of periodic boundary conditions. We have tested the influence of the size and number of bilayers on membrane properties using the Lipid14 force field for lipids in molecular dynamics simulations of 1,2-dioleoyl- sn-glycero-3-phosphocholine bilayers at full hydration (44 water molecules per lipid), low hydration (18 water molecules per lipid), and dehydration (9 water molecules per lipid). A number of additional simulations were conducted with the Slipids force field for comparison. We have found that the average area per lipid (APL), thickness, mass density profiles, and acyl tail order parameters are insensitive to the size and the number of bilayers for all hydration states. The Lipid14 force field can also successfully reproduce the experimentally observed decrease in APL and corresponding increase in bilayer thickness upon dehydration, reflecting the increase in ordering as the system becomes more gel-like. Additionally, decreasing hydration levels were associated with a trend away from normal lateral diffusion and toward more subdiffusive regimes across both force fields. In summary, at least for the Lipid14 force field, the use of a single bilayer with 128 phospholipid molecules provides an adequate representation of multi-bilayer systems at varying levels of hydration.

Curvature sensing by cardiolipin in simulated buckled membranes

Cardiolipin is a non-bilayer phospholipid with a unique dimeric structure. It localizes to negative curvature regions in bacteria and is believed to stabilize respiratory chain complexes in the highly curved mitochondrial membrane. Cardiolipin's localization mechanism remains unresolved, because important aspects such as the structural basis and strength for lipid curvature preferences are difficult to determine, partly due to the lack of efficient simulation methods. Here, we report a computational approach to study curvature preferences of cardiolipin by simulated membrane buckling and quantitative modeling. We combine coarse-grained molecular dynamics with simulated buckling to determine the curvature preferences in three-component bilayer membranes with varying concentrations of cardiolipin, and extract curvature-dependent concentrations and lipid acyl chain order parameter profiles. Cardiolipin shows a strong preference for negative curvatures, with a highly asymmetric chain order parameter profile. The concentration profiles are consistent with an elastic model for lipid curvature sensing that relates lipid segregation to local curvature via the material constants of the bilayers. These computations constitute new steps to unravel the molecular mechanism by which cardiolipin senses curvature in lipid membranes, and the method can be generalized to other lipids and membrane components as well.