Headgroup Structure and Cation Binding in Phosphatidylserine Lipid Bilayers

Elucidating the complex membrane binding of a protein with multiple anchoring domains using extHMMM

Membrane binding is a crucial mechanism for many proteins, but understanding the specific interactions between proteins and membranes remains a challenging endeavor. Coagulation factor Va (FVa) is a large protein whose membrane interactions are complicated due to the presence of multiple anchoring domains that individually can bind to lipid membranes. Using molecular dynamics simulations, we investigate the membrane binding of FVa and identify the key mechanisms that govern its interaction with membranes. Our results reveal that FVa can either adopt an upright or a tilted molecular orientation upon membrane binding. We further find that the domain organization of FVa deviates (sometimes significantly) from its crystallographic reference structure, and that the molecular orientation of the protein matches with domain reorganization to align the C2 domain toward its favored membrane-normal orientation. We identify specific amino acid residues that exhibit contact preference with phosphatidylserine lipids over phosphatidylcholine lipids, and we observe that mostly electrostatic effects contribute to this preference. The observed lipid-binding process and characteristics, specific to FVa or common among other membrane proteins, in concert with domain reorganization and molecular tilt, elucidate the complex membrane binding dynamics of FVa and provide important insights into the molecular mechanisms of protein-membrane interactions. An updated version of the HMMM model, termed extHMMM, is successfully employed for efficiently observing membrane bindings of systems containing the whole FVa molecule.

Revisiting the reaction pathways for phospholipid hydrolysis catalyzed by phospholipase A2 with QM/MM methods

Secreted phospholipase A2 (sPLA2) is a Ca2+-dependent, widely distributed enzyme superfamily in almost all mammalian tissues and bacteria. It is also a critical component of the venom of nearly all snakes, as well as many invertebrate species. In non-venomous contexts, sPLA2 assumes significance in cellular signaling pathways by binding cell membranes and catalyzing ester bond hydrolysis at the sn-2 position of phospholipids. Elevated levels of GIIA sPLA2 have been detected in the synovial fluid of arthritis patients, where it exhibits a pro-inflammatory function. Consequently, identifying sPLA2 inhibitors holds promise for creating highly effective pharmaceutical treatments. Beyond arthritis, the similarities among GIIA sPLA2s offer an opportunity for developing treatments against snakebite envenoming, the deadliest neglected tropical disease. Despite decades of study, the details of PLA2 membrane-binding, substrate-binding, and reaction mechanism remain elusive, demanding a comprehensive understanding of the sPLA2 catalytic mechanism. This study explores two reaction mechanism hypotheses, involving one or two water molecules, and distinct roles for the Ca2+ cofactor. Our research focuses on the human synovial sPLA2 enzyme bound to lipid bilayers of varying phospholipid compositions, and employing adiabatic QM/MM and QM/MM MD umbrella sampling methods to energetically and geometrically characterize the structures found along both reaction pathways. Our studies demonstrate the higher frequency of productive conformations within the single-water pathway, also revealing a lower free energy barrier for hydrolyzing POPC. Furthermore, we observe that the TS of this concerted one-step reaction closely resembles transition state geometries observed in X-ray crystallography complexes featuring high-affinity transition state analogue inhibitors.

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.

Membrane binding and lipid-protein interaction of the C2 domain from coagulation factor V

Anchoring of coagulation factors to anionic regions of the membrane involves the C2 domain as a key player. The rate of enzymatic reactions of the coagulation factors is increased by several orders of magnitude upon membrane binding. However, the precise mechanisms behind the rate acceleration remain unclear, primarily because of a lack of understanding of the conformational dynamics of the C2-containing factors and corresponding complexes. We elucidate the membrane-bound form of the C2 domain from human coagulation factor V (FV-C2) by characterizing its membrane binding the specific lipid-protein interactions. Employing all-atom molecular dynamics simulations and leveraging the highly mobile membrane-mimetic (HMMM) model, we observed spontaneous binding of FV-C2 to a phosphatidylserine (PS)-containing membrane within 2-25 ns across twelve independent simulations. FV-C2 interacted with the membrane through three loops (spikes 1-3), achieving a converged, stable orientation. Multiple HMMM trajectories of the spontaneous membrane binding provided extensive sampling and ample data to examine the membrane-induced effects on the conformational dynamics of C2 as well as specific lipid-protein interactions. Despite existing crystal structures representing presumed "open" and "closed" states of FV-C2, our results revealed a continuous distribution of structures between these states, with the most populated structures differing from both "open" and "closed" states observed in crystal environments. Lastly, we characterized a putative PS-specific binding site formed by K23, Q48, and S78 located in the groove enclosed by spikes 1-3 (PS-specificity pocket), suggesting a different orientation of a bound headgroup moiety compared to previous proposals based upon analysis of static crystal structures.

Elucidating the Membrane Binding Process of a Disordered Protein: Dynamic Interplay of Anionic Lipids and the Polybasic Region

Intrinsically disordered regions of proteins are responsible for many biological processes such as in the case of liver kinase B1 (LKB1)-a serine/threonine kinase relevant for cell proliferation and cell polarity. LKB1 becomes fully activated upon recruitment to the plasma membrane by binding of its disordered C-terminal polybasic motif consisting of eight lysines/arginines to phospholipids. Here, we present extensive molecular dynamics (MD) simulations of the polybasic motif interacting with a model membrane composed of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) and 1-palmitoyl-2-oleyl phosphatidic acid (PA) and cell culture experiments. Protein-membrane binding effects are due to the electrostatic interactions between the polybasic amino acids and PAs. For significant binding, the first three lysines turn out to be dispensable, which was also recapitulated in cell culture using transfected GFP-LKB1 variants. LKB1-membrane binding results in nonmonotonous changes in the structure of the protein as well as the membrane, in particular, accumulation of PAs and reduced thickness at the protein-membrane contact area. The protein-lipid binding turns out to be highly dynamic due to an interplay of PA-PA repulsion and protein-PA attraction. The thermodynamics of this interplay is captured by a statistical fluctuation model, which allows the estimation of both energies. Quantification of the significance of each polar amino acid in the polybasic provides detailed insights into the molecular mechanism of protein-membrane binding of LKB1. These results can likely be transferred to other proteins, which interact by intrinsically disordered polybasic regions with anionic membranes.

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.

Deciphering Lipid Arrangement in Phosphatidylserine/Phosphatidylcholine Mixed Membranes: Simulations and Experiments

Phosphatidylserine (PS) exposure on the plasma membrane is crucial for many cellular processes including apoptotic cell recognition, blood clotting regulation, cellular signaling, and intercellular interactions. In this study, we investigated the arrangement of PS headgroups in mixed PS/phosphatidylcholine (PC) bilayers, serving as a simplified model of the outer leaflets of mammalian cell plasma membranes. Combining atomistic-scale molecular dynamics (MD) simulations with Langmuir monolayer experiments, we unraveled the mutual miscibility of POPC and POPS lipids and the intricate intermolecular interactions inherent to these membranes as well as the disparities in position and orientation of PC and PS headgroups. Our experiments revealed micrometer-scale miscibility at all mole fractions of POPC and POPS, marked by modest deviations from ideal mixing with no apparent microscale phase separation. The MD simulations, meanwhile, demonstrated that these deviations were due to strong electrostatic interactions between like-lipid pairs (POPC-POPC and POPS-POPS), culminating in lateral segregation and nanoscale clustering. Notably, PS headgroups profoundly affect the ordering of the lipid acyl chains, leading to lipid elongation and subtle PS protrusion above the zwitterionic membrane. In addition, PC headgroups are more tilted with respect to the membrane normal, while PS headgroups align at a smaller angle, making them more exposed to the surface of the mixed PC/PS membranes. These findings provide a detailed molecular-level account of the organization of mixed PC/PS membranes, corroborated by experimental data. The insights gained here extend our comprehension of the physiological role of PSs.

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.

Efficient Simulations of Solvent Asymmetry Across Lipid Membranes Using Flat-Bottom Restraints

The routinely employed periodic boundary conditions complicate molecular simulations of physiologically relevant asymmetric lipid membranes together with their distinct solvent environments. Therefore, separating the extracellular fluid from its cytosolic counterpart has often been performed using a costly double-bilayer setup. Here, we demonstrate that the lipid membrane and solvent asymmetry can be efficiently modeled with a single lipid bilayer by applying an inverted flat-bottom potential to ions and other solute molecules, thereby restraining them to only interact with the relevant leaflet. We carefully optimized the parameters of the suggested method so that the results obtained using the flat-bottom and double-bilayer approaches become mutually indistinguishable. Then, we apply the flat-bottom approach to lipid bilayers with various compositions and solvent environments, covering ions and cationic peptides to validate the approach in a realistic use case. We also discuss the possible limitations of the method as well as its computational efficiency and provide a step-by-step guide on how to set up such simulations in a straightforward manner.

The Effect of Cholesterol in SOPC Lipid Bilayers at Low Temperatures

We study the behavior of lipid bilayers composed of SOPC (1-stearoyl-2-oleoyl-sn-glycero-3-phosphocholine) with different concentrations of cholesterol, ranging from 10 mol% to 50 mol% at 273 K. To this end, we carry out extensive atomistic molecular dynamic simulations with the aid of the Slipid force field aiming at computing basic bilayer parameters, as well as thermodynamic properties and structural characteristics. The obtained results are compared to available relevant experimental data and the outcome of atomistic simulations performed on bilayers composed of analogous phospholipids. Our results show a good quantitative, as well as qualitative, agreement with the main trends associated with the concentration increase in cholesterol. Moreover, it comes out that a change in the behavior of the bilayer is brought about at a concentration of about 30 mol% cholesterol. At this very concentration, some of the bilayer properties are found to exhibit a saturation and a significant long-range ordering of the lipid molecules in the membrane shows up.

Lipid Structure Determines the Differential Impact of Single Metal Additions and Binary Mixtures of Manganese, Calcium and Magnesium on Membrane Fluidity and Liposome Size

Unilamellar vesicles of the biologically relevant lipids phosphatidic acid (PA) and phosphatidylserine (PS) with fully saturated (DM-) or partly unsaturated (PO-) acyl side chains were exposed to Ca, Mn and Mg in single metal additions; in equimolar mixtures or by sequential additions of one metal at a time. Laurdan generalized polarization measured the membrane fluidity, while dynamic light scattering reported liposome size changes complemented by zeta potential. All metals induced membrane rigidity and increased liposome sizes across all systems. Mn had the strongest effect overall, but Mg was comparable for DMPS. Lipid side chain architecture was important as GP values for binary mixtures were higher than expected from the sum of values for single additions added to POPS but smaller for DMPS. Sequential additions were predominantly different for Ca:Mg mixtures. Mn induced the strongest increase of liposome size in saturated lipids whereas Ca effects dominated unsaturated matrices. Binary additions induced larger sizes than the sum of single additions for POPS, but much lower changes in DMPA. The order of addition was relevant for PS systems. Thus, lipid structure determines metal effects, but their impact is modulated by other ions. Thus, metal effects may differ with the local lipid architecture and metal concentrations within cells.

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.

Charged Small Molecule Binding to Membranes in MD Simulations Evaluated against NMR Experiments

Interactions of charged molecules with biomembranes regulate many of their biological activities, but their binding affinities to lipid bilayers are difficult to measure experimentally and model theoretically. Classical molecular dynamics (MD) simulations have the potential to capture the complex interactions determining how charged biomolecules interact with membranes, but systematic overbinding of sodium and calcium cations in standard MD simulations raises the question of how accurately force fields capture the interactions between lipid membranes and charged biomolecules. Here, we evaluate the binding of positively charged small molecules, etidocaine, and tetraphenylphosphonium to a phosphatidylcholine (POPC) lipid bilayer using the changes in lipid head-group order parameters. We observed that these molecules behave oppositely to calcium and sodium ions when binding to membranes: (i) their binding affinities are not overestimated by standard force field parameters, (ii) implicit inclusion of electronic polarizability increases their binding affinity, and (iii) they penetrate into the hydrophobic membrane core. Our results can be explained by distinct binding mechanisms of charged small molecules with hydrophobic moieties and monoatomic ions. The binding of the former is driven by hydrophobic effects, while the latter has direct electrostatic interactions with lipids. In addition to elucidating how different kinds of charged biomolecules bind to membranes, we deliver tools for further development of MD simulation parameters and methodology.

The lung surfactant activity probed with molecular dynamics simulations

The surface of pulmonary alveolar subphase is covered with a mixture of lipids and proteins. This lung surfactant plays a crucial role in lung functioning. It shows a complex phase behavior which can be altered by the interaction with third molecules such as drugs or pollutants. For studying multicomponent biological systems, it is of interest to couple experimental approach with computational modelling yielding atomic-scale information. Simple two, three, or four-component model systems showed to be useful for getting more insight in the interaction between lipids, lipids and proteins or lipids and proteins with drugs and impurities. These systems were studied theoretically using molecular dynamic simulations and experimentally by means of the Langmuir technique. A better understanding of the structure and behavior of lung surfactants obtained from this research is relevant for developing new synthetic surfactants for efficient therapies, and may contribute to public health protection.

Lipid headgroup and side chain architecture determine manganese-induced dose dependent membrane rigidification and liposome size increase

Metal ion-membrane interactions have gained appreciable attention over the years resulting in increasing investigations into the mode of action of toxic and essential metals. More work has focused on essential ions like Ca or Mg and toxic metals like Cd and Pb, whereas this study investigates the effects of the abundant essential trace metal manganese with model lipid systems by screening zwitterionic and anionic glycerophospholipids. Despite its essentiality, deleterious impact towards cell survival is known under Mn stress. The fluorescent dyes Laurdan and diphenylhexatriene were used to assess changes in membrane fluidity both in the head group and hydrophobic core region of the membrane, respectively. Mn-rigidified membranes composed of the anionic phospholipids, phosphatidic acid, phosphatidylglycerol, cardiolipin, and phosphatidylserine. Strong binding resulted in large shifts of the phase transition temperature. The increase was in the order phosphatidylserine > phosphatidylglycerol > cardiolipin, and in all cases, saturated analogues > mono-unsaturated forms. Dynamic light scattering measurements revealed that Mn caused extensive aggregation of liposomes composed of saturated analogues of phosphatidic acid and phosphatidylserine, whilst the mono-unsaturated analogue had significant membrane swelling. Increased membrane rigidity may interfere with permeability of ions and small molecules, possibly disrupting cellular homeostasis. Moreover, liposome size changes could indicate fusion, which could also be detrimental to cellular transport. Overall, this study provided further understanding into the effects of Mn with biomembranes, whereby the altered membrane properties are consequential to the proper structural and signalling functions of membrane lipids.

Bacterial Membranes Are More Perturbed by the Asymmetric Versus Symmetric Loading of Amphiphilic Molecules

Characterizing the biophysical properties of bacterial membranes is critical for understanding the protective nature of the microbial envelope, interaction of biological membranes with exogenous materials, and designing new antibacterial agents. Presented here are molecular dynamics simulations for two cationic quaternary ammonium compounds, and the anionic and nonionic form of a fatty acid molecule interacting with a Staphylococcus aureus bacterial inner membrane. The effect of the tested materials on the properties of the model membranes are evaluated with respect to various structural properties such as the lateral pressure profile, lipid tail order parameter, and the bilayer’s electrostatic potential. Conducting asymmetric loading of molecules in only one leaflet, it was observed that anionic and cationic amphiphiles have a large impact on the Staphylococcus aureus membrane’s electrostatic potential and lateral pressure profile as compared to a symmetric distribution. Nonintuitively, we find that the cationic and anionic molecules induce a similar change in the electrostatic potential, which points to the complexity of membrane interfaces, and how asymmetry can induce biophysical consequences. Finally, we link changes in membrane structure to the rate of electroporation for the membranes, and again find a crucial impact of introducing asymmetry to the system. Understanding these physical mechanisms provides critical insights and viable pathways for the rational design of membrane-active molecules, where controlling the localization is key.

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.

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.

Amyloid-β peptide dimers undergo a random coil to β-sheet transition in the aqueous phase but not at the neuronal membrane

Mounting evidence suggests that the neuronal cell membrane is the main site of oligomer-mediated neuronal toxicity of amyloid-β peptides in Alzheimer's disease. To gain a detailed understanding of the mutual interference of amyloid-β oligomers and the neuronal membrane, we carried out microseconds of all-atom molecular dynamics (MD) simulations on the dimerization of amyloid-β (Aβ)42 in the aqueous phase and in the presence of a lipid bilayer mimicking the in vivo composition of neuronal membranes. The dimerization in solution is characterized by a random coil to β-sheet transition that seems on pathway to amyloid aggregation, while the interactions with the neuronal membrane decrease the order of the Aβ42 dimer by attenuating its propensity to form a β-sheet structure. The main lipid interaction partners of Aβ42 are the surface-exposed sugar groups of the gangliosides GM1. As the neurotoxic activity of amyloid oligomers increases with oligomer order, these results suggest that GM1 is neuroprotective against Aβ-mediated toxicity.

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.

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.

Targeted liposomal drug delivery: a nanoscience and biophysical perspective

Liposomes are a unique platform for drug delivery, and a number of liposomal formulations have already been commercialized. Doxil is a representative example, which uses PEGylated liposomes to load doxorubicin for cancer therapy. Its delivery relies on the enhanced permeability and retention (EPR) effect or passive targeting. Drug loading can be achieved using both standard liposomes and also those containing a solid core such as mesoporous silica and poly(lactide-co-glycolide) (PLGA). Developments have also been made on active targeted delivery using bioaffinity ligands such as small molecules, antibodies, peptides and aptamers. Compared to other types of nanoparticles, the surface of liposomes is fluid, allowing dynamic organization of targeting ligands to achieve optimal binding to cell surface receptors. This review article summarizes development of liposomal targeted drug delivery systems, with an emphasis on the biophysical properties of lipids. In both passive and active targeting, the effects of liposome size, charge, fluidity, rigidity, head-group chemistry and PEGylation are discussed along with recent examples. Most of the examples are focused on targeting tumors or cancer cells. Finally, a few examples of commercialized formulations are described, and some future research opportunities are discussed.

In Silico Prediction of the Binding, Folding, Insertion, and Overall Stability of Membrane-Active Peptides

Membrane-active peptides (MAPs) are short-length peptides used for potential biomedical applications in diagnostic imaging of tissues, targeted drug delivery, gene delivery, and antimicrobials and antibiotics. The broad appeal of MAPs is that they are infinitely variable, relatively low cost, and biocompatible. However, experimentally characterizing the specific properties of a MAP or its many variants is a low-resolution and potentially time-consuming endeavor; molecular dynamics (MD) simulations have emerged as an invaluable tool in identifying the biophysical interactions that are fundamental to the function of MAPs. In this chapter, a step-by-step approach to discreetly model the binding, folding, and insertion of a membrane-active peptide to a model lipid bilayer using MD simulations is described. Detailed discussion is devoted to the critical aspects of running these types of simulations: prior knowledge of the system, understanding the strengths and weaknesses of molecular mechanics force fields, proper construction and equilibration of the system, realistically estimating both experimental and computational timescales, and leveraging analysis to make direct comparisons to experimental results as often as possible.

Polycation-Anionic Lipid Membrane Interactions

Natural or synthetic polycations are used as biocides or as drug/gene carriers. Understanding the interactions between these macromolecules and cell membranes at the molecular level is therefore of great importance for the design of effective polymer biocides or biocompatible polycation-based delivery systems. Until now, details of the processes at the interface between polycations and biological systems have not been fully recognized. In this study, we consider the effect of strong polycations with quaternary ammonium groups on the properties of anionic lipid membranes that we use as a model system for protein-free cell membranes. For this purpose, we employed experimental measurements and atomic-scale molecular dynamics (MD) simulations. MD simulations reveal that the polycations are strongly hydrated in the aqueous phase and do not lose the water shell after adsorption at the bilayer surface. As a result of strong hydration, the polymer chains reside at the phospholipid headgroup and do not penetrate to the acyl chain region. The polycation adsorption involves the formation of anionic lipid-rich domains, and the density of anionic lipids in these domains depends on the length of the polycation chain. We observed the accumulation of anionic lipids only in the leaflet interacting with the polymer, which leads to the formation of compositionally asymmetric domains. Asymmetric adsorption of the polycation on only one leaflet of the anionic membrane strongly affects the membrane properties in the polycation-membrane contact areas: (i) anionic lipid accumulates in the region near the adsorbed polymer, (ii) acyl chain ordering and lipid packing are reduced, which results in a decrease in the thickness of the bilayer, and (iii) polycation-anionic membrane interactions are strongly influenced by the presence and concentration of salt. Our results provide an atomic-scale description of the interactions of polycations with anionic lipid bilayers and are fully supported by the experimental data. The outcomes are important for understanding the correlation of the structure of polycations with their activity on biomembranes.

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.

Simulations of Asymmetric Membranes Illustrate Cooperative Leaflet Coupling and Lipid Adaptability

Biological membranes are composed of lipid bilayers that are often asymmetric with regards to the lipid composition and/or aqueous solvent they separate. Studying lipid asymmetry both experimentally and computationally is challenging. Molecular dynamics simulations of lipid bilayers with asymmetry are difficult due to finite system sizes and time scales accessible to simulations. Due to the very slow flip-flop rate for phospholipids, one must first choose how many lipids are on each side of the bilayer, but the resulting bilayer may be unstable (or metastable) due to differing tensile and compressive forces between leaflets. Here we use molecular dynamics simulations to investigate a number of different asymmetric membrane systems, both with atomistic and coarse-grained models. Asymmetries studied include differences in number of lipids, lipid composition (unsaturated and saturated tails and different headgroups), and chemical gradients between the aqueous phases. Extensive analysis of the bilayers' properties such as area per lipid, density, and lateral pressure profiles are used to characterize bilayer asymmetry. We also address how cholesterol (which flip-flops relatively quickly) influences membrane asymmetries. Our results show how each leaflet is influenced by the other and can mitigate the structural changes to the bilayer overall structure. Cholesterol can respond to changes in bilayer asymmetry to alleviate some of the effect on the bilayer structure, but that will alter its leaflet distribution, which in turn affects its chemical potential. Ionic imbalances are shown to have a modest change in bilayer structure, despite large changes in the electrostatic potential. Bilayer asymmetry can also induce a modest electrostatic potential across the membrane. Our results highlight the importance of membrane asymmetry on bilayer properties, the influence of lipid headgroups, tails and cholesterol on asymmetry, and the ability of lipids to adapt to different environments.

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.

Calcium-Lipid Interactions Observed with Isotope-Edited Infrared Spectroscopy

Calcium ions bind to lipid membranes containing anionic lipids; however, characterizing the specific ion-lipid interactions in multicomponent membranes has remained challenging because it requires nonperturbative lipid-specific probes. Here, using a combination of isotope-edited infrared spectroscopy and molecular dynamics simulations, we characterize the effects of a physiologically relevant (2 mM) Ca²⁺ concentration on zwitterionic phosphatidylcholine and anionic phosphatidylserine lipids in mixed lipid membranes. We show that Ca²⁺ alters hydrogen bonding between water and lipid headgroups by forming a coordination complex involving the lipid headgroups and water. These interactions distort interfacial water orientations and prevent hydrogen bonding with lipid ester carbonyls. We demonstrate, experimentally, that these effects are more pronounced for the anionic phosphatidylserine lipids than for zwitterionic phosphatidylcholine lipids in the same membrane.

Functionalization of the Parylene C Surface Enhances the Nucleation of Calcium Phosphate: Combined Experimental and Molecular Dynamics Simulations Approach

Interactions at the solid-body fluid interfaces play a vital role in bone tissue formation at the implant surface. In this study, fully atomistic molecular dynamics (MD) simulations were performed to investigate interactions between the physiological components of body fluids (Ca2+, HPO42-, H2PO4-, Na+, Cl-, and H2O) and functionalized parylene C surface. In comparison to the native parylene C (-Cl surface groups), the introduction of -OH, -CHO, and -COOH surface groups significantly enhances the interactions between body fluid ions and the polymeric surface. The experimentally observed formation of calcium phosphate nanocrystals is discussed in terms of MD simulations of the calcium phosphate clustering. Surface functional groups promote the clustering of calcium and phosphate ions in the following order: -OH > -CHO > -Cl (parent parylene C) ≈ -COO-. This promoting role of surface functional groups is explained as stimulating the number of Ca2+ and HPO42- surface contacts as well as ion chemisorption. The molecular mechanism of calcium phosphate cluster formation at the functionalized parylene C surface is proposed.

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.

Accurate Biomolecular Simulations Account for Electronic Polarization

In this perspective, we discuss where and how accounting for electronic many-body polarization affects the accuracy of classical molecular dynamics simulations of biomolecules. While the effects of electronic polarization are highly pronounced for molecules with an opposite total charge, they are also non-negligible for interactions with overall neutral molecules. For instance, neglecting these effects in important biomolecules like amino acids and phospholipids affects the structure of proteins and membranes having a large impact on interpreting experimental data as well as building coarse grained models. With the combined advances in theory, algorithms and computational power it is currently realistic to perform simulations with explicit polarizable dipoles on systems with relevant sizes and complexity. Alternatively, the effects of electronic polarization can also be included at zero additional computational cost compared to standard fixed-charge force fields using the electronic continuum correction, as was recently demonstrated for several classes of biomolecules.