Divalent cations promote huntingtin fibril formation on endoplasmic reticulum derived and model membranes

Huntington's Disease (HD) is caused by an abnormal expansion of the polyglutamine (polyQ) domain within the first exon of the huntingtin protein (htt). This expansion promotes disease-related htt aggregation into amyloid fibrils and the formation of proteinaceous inclusion bodies within neurons. Fibril formation is a complex heterogenous process involving an array of aggregate species such as oligomers, protofibrils, and fibrils. In HD, structural abnormalities of membranes of several organelles develop. In particular, the accumulation of htt fibrils near the endoplasmic reticulum (ER) impinges upon the membrane, resulting in ER damage, altered dynamics, and leakage of Ca2+. Here, the aggregation of htt at a bilayer interface assembled from ER-derived liposomes was investigated, and fibril formation directly on these membranes was enhanced. Based on these observations, simplified model systems were used to investigate mechanisms associated with htt aggregation on ER membranes. As the ER-derived liposome fractions contained residual Ca2+, the role of divalent cations was also investigated. In the absence of lipids, divalent cations had minimal impact on htt structure and aggregation. However, the presence of Ca2+ or Mg2+ played a key role in promoting fibril formation on lipid membranes despite reduced htt insertion into and association with lipid interfaces, suggesting that the ability of divalent cations to promote fibril formation on membranes is mediated by induced changes to the lipid membrane physicochemical properties. With enhanced concentrations of intracellular calcium being a hallmark of HD, the ability of divalent cations to influence htt aggregation at lipid membranes may play a role in aggregation events that lead to organelle abnormalities associated with disease.

Exploring the Properties of Curved Lipid Membranes: Comparative Analysis of Atomistic and Coarse-Grained Force Fields

Curvature emerges as a fundamental membrane characteristic crucial for diverse biological processes, including vesicle formation, cell signaling, and membrane trafficking. Increasingly valuable insights into atomistic details governing curvature-dependent membrane properties are provided by computer simulations. Nevertheless, the underlying force field models are conventionally calibrated and tested in relation to experimentally derived parameters of planar bilayers, thereby leaving uncertainties concerning their consistency in reproducing curved lipid systems. In this study we compare the depiction of buckled phosphatidylcholine (POPC) and POPC-cholesterol membranes by four popular force field models. Aside from agreement with respect to general trends in curvature dependence of a number of parameters, we observe a few qualitative differences. Among the most prominent ones is the difference between atomistic and coarse grained force fields in their representation of relative compressibility of the polar headgroup region and hydrophobic lipid core. Through a number of downstream effects, this discrepancy can influence the way in which curvature modulates the behavior of membrane bound proteins depending on the adopted simulation model.

Lateral diffusion of ions near membrane surface

Biological membranes isolate living cells from their environment, while allowing selective molecular transport between the inner and outer realms. For example, Na+ and K+ permeability through ionic channels contributes to neural conduction. Whether the ionic currents arise directly from cations in the bulk, or from the interface, is currently unclear. There are only scant results concerning lateral diffusion of ions on aquated membrane surfaces (and strong belief that this occurs through binding to a diffusing lipid). We performed classical molecular dynamics (MD) simulations of monovalent ions, Na+, K+, and Cl-, near the surface of the zwitterionic palmitoyl-oleoyl-phosphatidylcholine (POPC) membrane. Realistic force-fields for lipids (Amber's Lipid17 and Lipid21) and water (TIP4P-Ew) are tested for the mass and charge densities and the electrostatic potential across the membrane. These calculations reveal that the chloride can bind to the choline moiety through an intervening water molecule by forming a CH⋯OH hydrogen bond, while cations bind to both the phosphatic and carbonyl oxygens of phosphatidylcholine moieties. Upon transitioning from the bulk to the interface, a cation sheds some of its hydration water, which are replaced by headgroup atoms. Notably, an interfacial cation can bind 1-4 headgroup atoms, which is a key to understanding its surface hopping mechanism. We find that cation binding to three headgroup atoms immobilizes it, while binding to four energizes it. Consequently, the lateral cation diffusion rate is only 15-25 times slower than in the bulk, and 4-5 times faster than lipid self-diffusion. K+ diffusion is notably more anomalous than Na+, switching from sub- to super-diffusion after about 2 ns.

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.

Photoswitchable Endocytosis of Biomolecular Condensates in Giant Vesicles

Interactions between membranes and biomolecular condensates can give rise to complex phenomena such as wetting transitions, mutual remodeling, and endocytosis. In this study, light-triggered manipulation of condensate engulfment is demonstrated using giant vesicles containing photoswitchable lipids. UV irradiation increases the membrane area, which can be stored in nanotubes. When in contact with a condensate droplet, the UV light triggers rapid condensate endocytosis, which can be reverted by blue light. The affinity of the protein-rich condensates to the membrane and the reversibility of the engulfment processes is quantified from confocal microscopy images. The degree of photo-induced engulfment, whether partial or complete, depends on the vesicle excess area and the relative sizes of vesicles and condensates. Theoretical estimates suggest that utilizing the light-induced excess area to increase the vesicle-condensate adhesion interface is energetically more favorable than the energy gain from folding the membrane into invaginations and tubes. The overall findings demonstrate that membrane-condensate interactions can be easily and quickly modulated via light, providing a versatile system for building platforms to control cellular events and design intelligent drug delivery systems for cell repair.

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.

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.

Lipid and cholesterols modulate the dynamics of SARS-CoV-2 viral ion channel ORF3a and its pathogenic variants

SARS-CoV-2 accessory protein, ORF3a is a putative ion channel which immensely contributes to viral pathogenicity by modulating host immune responses and virus-host interactions. Relatively high expression of ORF3a in diseased individuals and implication with inflammasome activation, apoptosis and autophagy inhibition, ratifies as an effective target for developing vaccines and therapeutics. Herein, we present the elusive dynamics of ORF3a-dimeric state using all-atoms molecular dynamics (MD) simulations at μ-seconds scale in a heterogeneous lipid-mimetic system in multiple replicates. Additionally, we also explore the effect of non-synonymous pathogenic mutations on ORF3a ion channel activity and viral pathogenicity in different SARS-CoV-2 variants using various structure-based protein stability (ΔΔG) tools and computational saturation mutagenesis. Our study ascertains the role of phosphatidylcholines and cholesterol in modulating the structure of ORF3a, which perturbs the size and flexibility of the polar cavity that allows permeation of large cations. Discrete trend in ion channel pore radius and area per lipid arises the premise that presence of lipids might also affect the overall conformation of ORF3a. MD structural-ensembles, in some replicates rationalize the crucial role of TM2 in maintaining the native structure of ORF3a. We also infer that loss of structural stability primarily grounds for pathogenicity in more than half of the pathogenic variants of ORF3a. Overall, the effect of mutation on alteration of ion permeability of ORF3a, proposed in this study brings mechanistic insights into variant consequences on viral membrane proteins of SARS-CoV-2, which can be utilized for the development of novel therapeutics to treat COVID-19 and other coronavirus diseases.

Accurately Predicting Protein pKa Values Using Nonequilibrium Alchemy

The stability, solubility, and function of a protein depend on both its net charge and the protonation states of its individual residues. pKa is a measure of the tendency for a given residue to (de)protonate at a specific pH. Although pKa values can be resolved experimentally, theory and computation provide a compelling alternative. To this end, we assess the applicability of a nonequilibrium (NEQ) alchemical free energy method to the problem of pKa prediction. On a data set of 144 residues that span 13 proteins, we report an average unsigned error of 0.77 ± 0.09, 0.69 ± 0.09, and 0.52 ± 0.04 pK for aspartate, glutamate, and lysine, respectively. This is comparable to current state-of-the-art predictors and the accuracy recently reached using free energy perturbation methods (e.g., FEP+). Moreover, we demonstrate that our open-source, pmx-based approach can accurately resolve the pKa values of coupled residues and observe a substantial performance disparity associated with the lysine partial charges in Amber14SB/Amber99SB*-ILDN, for which an underused fix already exists.

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.

Amyloid-β Tetramers and Divalent Cations at the Membrane/Water Interface: Simple Models Support a Functional Role

Charge polarization at the membrane interface is a fundamental process in biology. Despite the lower concentration compared to the abundant monovalent ions, the relative abundance of divalent cations (Ca2+, Mg2+, Zn2+, Fe2+, Cu2+) in particular spaces, such as the neuron synapse, raised many questions on the possible effects of free multivalent ions and of the required protection of membranes by the eventual defects caused by the free forms of the cations. In this work, we first applied a recent realistic model of divalent cations to a well-investigated model of a polar lipid bilayer, di-myristoyl phosphatidyl choline (DMPC). The full atomistic model allows a fairly good description of changes in the hydration of charged and polar groups upon the association of cations to lipid atoms. The lipid-bound configurations were analyzed in detail. In parallel, amyloid-β 1-42 (Aβ42) peptides assembled into tetramers were modeled at the surface of the same bilayer. Two of the protein tetramers' models were loaded with four Cu2+ ions, the latter bound as in DMPC-free Aβ42 oligomers. The two Cu-bound models differ in the binding topology: one with each Cu ion binding each of the monomers in the tetramer; one with pairs of Cu ions linking two monomers into dimers, forming tetramers as dimers of dimers. The models here described provide hints on the possible role of Cu ions in synaptic plasticity and of Aβ42 oligomers in storing the same ions away from lipids. The release of structurally disordered peptides in the synapse can be a mechanism to recover ion homeostasis and lipid membranes from changes in the divalent cation concentration.

Effect of cholesterol on the ion-membrane interaction: Zeta potential and dynamic light scattering study

Cholesterol in a bio-membrane plays a significant role in many cellular event and is known to regulate the functional activity of protein and ion channel. In this study we report a significant effect of cholesterol on the ion-membrane interaction. We prepare large unilamellar vesicles, composed of zwitterionic lipid DOPC and anionic lipid DOPG with different cholesterol concentration. Electrostatics of anionic membranes containing cholesterol in the presence of NaCl has systematically been explored using dynamic light scattering and zeta potential. Negative zeta potential of the membrane decreases its negative value with increasing ion concentration for all cholesterol concentrations. However, zeta potential itself decreases with increasing cholesterol content even in the absence of monovalent ions. Electrostatic behaviour of the membrane is determined from well-known Gouy Chapmann model. Negative surface charge density of the membrane decreases with increasing cholesterol content. Binding constant, estimated from the electrostatic double layer theory, is found to increase significantly in the presence of cholesterol. Comparison of electrostatic parameters of the membrane in the presence and absence of cholesterol suggests that cholesterol significantly alter the electrostatic behaviour of the membrane.

Peptide-Membrane Binding: Effects of the Amino Acid Sequence

An understanding of how the amino acid sequence affects the interaction of peptides with lipid membranes remains mostly unknown. This type of knowledge is required to rationalize membrane-induced toxicity of amyloid peptides and to design peptides that can interact with lipid bilayers. Here, we perform a systematic study of how variations in the sequence of the amphipathic Ac-(FKFE)₂-NH₂ peptide affect its interaction with zwitterionic lipid bilayers using extensive all-atom molecular dynamics simulations in explicit solvent. Our results show that peptides with a net positive charge bind more frequently to the lipid bilayer than neutral or negatively charged sequences. Moreover, neutral amphipathic peptides made with the same numbers of phenylalanine (F), lysine (K), and glutamic (E) amino acids at different positions in the sequence differ significantly in their frequency of binding to the membrane. We find that peptides bind with a higher frequency to the membrane if their positive lysine side chains are more exposed to the solvent, which occurs if they are located at the extremity (as opposed to the middle) of the sequence. Non-polar residues play an important role in accounting for the adsorption of peptides onto the membrane. In particular, peptides made with less hydrophobic non-polar residues (e.g., valine and alanine) are significantly less adsorbed to the membrane compared to peptides made with phenylalanine. We also find that sequences where phenylalanine residues are located at the extremities of the peptide have a higher tendency to be adsorbed.

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.

Ionic Strength and Solution Composition Dictate the Adsorption of Cell-Penetrating Peptides onto Phosphatidylcholine Membranes

Adsorption of arginine-rich positively charged peptides onto neutral zwitterionic phosphocholine (PC) bilayers is a key step in the translocation of those potent cell-penetrating peptides into the cell interior. In the past, we have shown both theoretically and experimentally that polyarginines adsorb to the neutral PC-supported lipid bilayers in contrast to polylysines. However, comparing our results with previous studies showed that the results often do not match even at the qualitative level. The adsorption of arginine-rich peptides onto 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) may qualitatively depend on the actual experimental conditions where binding experiments have been performed. In this work, we systematically studied the adsorption of R₉ and K₉ peptides onto the POPC bilayer, aided by molecular dynamics (MD) simulations and fluorescence cross-correlation spectroscopy (FCCS) experiments. Using MD simulations, we tested a series of increasing peptide concentrations, in parallel with increasing Na⁺ and Ca²⁺ salt concentrations, showing that the apparent strength of adsorption of R₉ decreases upon the increase of peptide or salt concentration in the system. The key result from the simulations is that the salt concentrations used experimentally can alter the picture of peptide adsorption qualitatively. Using FCCS experiments with fluorescently labeled R₉ and K₉, we first demonstrated that the binding of R₉ to POPC is tighter by almost 2 orders of magnitude compared to that of K₉. Finally, upon the addition of an excess of either Na⁺ or Ca²⁺ ions with R₉, the total fluorescence correlation signal is lost, which implies the unbinding of R₉ from the PC bilayer, in agreement with our predictions from MD simulations.

The role of ion-lipid interactions and lipid packing in transient defects caused by phenolic compounds

The transient disruption of membranes for the passive permeation of ions or small molecules is a complex process relevant to understanding physiological processes and biotechnology applications. Phenolic compounds are widely studied for their antioxidant and antimicrobial properties, and some of these activities are based on the interactions of the phenolic compound with membranes. Ions are ubiquitous in cells and are known to alter the structure of phospholipid bilayers. Yet, ion-lipid interactions are usually ignored when studying the membrane-altering properties of phenolic compounds. This study aims to assess the role of Ca2+ ions on the membrane-disrupting activity of two phenolic acids and to highlight the role of local changes in lipid packing in forming transient defects or pores. Results from tethered bilayer lipid membrane electrical impedance spectroscopy experiments showed that Ca2+ significantly reduces membrane disruption by caffeic acid methyl ester and caffeic acid. As phenolic acids are known metal chelators, we used UV-vis and fluorescence spectroscopy to exclude the possibility that Ca2+ interferes with membrane disruption by binding to the phenolic compound and subsequently preventing membrane binding. Molecular dynamics simulations showed that Ca2+ but not caffeic acid methyl ester or caffeic acid increases lipid packing in POPC bilayers. The combined data confirm that Ca2+ reduces the membrane-disrupting activity of the phenolic compounds, and that Ca2+-induced changes to lipid packing govern this effect. We discuss our data in the context of ion-induced pores and transient defects and how lipid packing affects membrane disruption by small molecules.

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.

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.

Small ionic radii limit time step in Martini 3 molecular dynamics simulations

Among other improvements, the Martini~3 coarse-grained force field provides a more accurate description of the solvation of protein pockets and channels through the consistent use of various bead types and sizes. Here, we show that the representation of Na+ and Cl- ions as ``tiny' (TQ5) beads limits the accessible time step to 25~fs. By contrast, with Martini~2, time steps of 30-40~fs were possible for lipid bilayer systems without proteins. This limitation is relevant for systems that require long equilibration times. We derive a quantitative kinetic model of time-integration instabilities in molecular dynamics as a function of time step, ion concentration and mass, system size, and simulation time. We demonstrate that ion-water interactions are the main source of instability at physiological conditions, followed closely by ion-ion interactions. We show that increasing the ionic masses makes it possible to use time steps up to 40~fs with minimal impact on static equilibrium properties and on dynamical quantities such as lipid and ion diffusion coefficients. Increasing the size of the bead representing the ions (and thus changing their hydration) also permits longer time steps. For a soluble protein, we find that increasing the mass of tiny beads also on the protein permits simulations with 30-fs time step. The use of larger time steps in Martini~3 results in a more efficient exploration of configuration space. The kinetic model of MD simulation crashes can be used to determine the maximum allowed time step upfront for an efficient use of resources and whenever sampling efficiency is critical.

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.

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.

Changes in the Local Conformational States Caused by Simple Na+ and K+ Ions in Polyelectrolyte Simulations: Comparison of Seven Force Fields with and without NBFIX and ECC Corrections

Electrostatic interactions have a determining role in the conformational and dynamic behavior of polyelectrolyte molecules. In this study, anionic polyelectrolyte molecules, poly(glutamic acid) (PGA) and poly(aspartic acid) (PASA), in a water solution with the most commonly used K+ or Na+ counterions, were investigated using atomistic molecular dynamics (MD) simulations. We performed a comparison of seven popular force fields, namely AMBER99SB-ILDN, AMBER14SB, AMBER-FB15, CHARMM22*, CHARMM27, CHARMM36m and OPLS-AA/L, both with their native parameters and using two common corrections for overbinding of ions, the non-bonded fix (NBFIX), and electronic continuum corrections (ECC). These corrections were originally introduced to correct for the often-reported problem concerning the overbinding of ions to the charged groups of polyelectrolytes. In this work, a comparison of the simulation results with existing experimental data revealed several differences between the investigated force fields. The data from these simulations and comparisons with previous experimental data were then used to determine the limitations and strengths of these force fields in the context of the structural and dynamic properties of anionic polyamino acids. Physical properties, such as molecular sizes, local structure, and dynamics, were studied using two types of common counterions, namely potassium and sodium. The results show that, in some cases, both the macroion size and dynamics depend strongly on the models (parameters) for the counterions due to strong overbinding of the ions and charged side chain groups. The local structures and dynamics are more sensitive to dihedral angle parameterization, resulting in a preference for defined monomer conformations and the type of correction used. We also provide recommendations based on the results.

Calcium Ion Binding at the Lipid-Water Interface Alters the Ion Permeability of Phospholipid Bilayers

Calcium ions (Ca²⁺) play a fundamental role in membrane-associated physiological processes. Ca²⁺ can also significantly modulate the physicochemical properties of phospholipid bilayers, but whether this occurs at physiologically relevant concentrations is difficult to determine because of the uncertainty in the reported affinity of Ca²⁺ for phospholipid bilayers. In this article, we determine the apparent affinity of Ca²⁺ for zwitterionic phospholipid bilayers using tethered bilayer lipid membranes (tBLMs) used in conjunction with swept-frequency electrical impedance spectroscopy (EIS). We report that Ca²⁺ binds to phospholipid bilayers at physiologically relevant concentrations and modulates membrane permeability. We present direct experimental evidence that this effect is governed by specific interactions with select lipid headgroup moieties, which is supported by data from molecular dynamics (MD) simulations. This is the first reported use of tBLM/EIS to estimate cation-membrane affinity. Combined with MD simulations, this technique provides a novel methodology to elucidate the molecular details of cation-membrane interactions at the water-phospholipid interface.

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.

Interplay of Interfacial Viscosity, Specific-Ion, and Impurity Adsorption Determines Zeta Potentials of Phospholipid Membranes

Ion-specific induced changes of the ζ-potential of phospholipid vesicles are commonly used to quantify the affinity of different ions to the lipid interface. The negative ζ-potential of zwitterionic net-neutral phospholipid vesicles in neat water, which changes sign and increases in solutions of NaCl or KCl, is a phenomenon consistently observed in experiments but not fully understood theoretically. Using atomistic molecular dynamics simulations in the presence of applied electric fields which drive electroosmotic flows, in combination with an electrostatic continuum model based on the modified Poisson-Boltzmann and Helmholtz-Smoluchowski equations, we study the electrokinetic and electrostatic properties as well as the specific ion affinities to the phospholipid-water interface, in order to resolve these puzzling observations. Our modified continuum equations account for the dielectric profile at the lipid-water interface, ion-specific interactions between ions and the lipid-water interface, and the interfacial viscosity profile, which are all extracted from our atomistic simulations and rather accurately predict ion-density and electrostatic-potential distributions as well as ζ-potentials in comparison with our atomistic simulations. Our continuum model can explain experimental ζ-potentials only when we assume minute amounts of surface-active anionic impurities in the aqueous solution. In fact, the amount of impurities needed to explain the experimental data increases linearly with the salt concentration, suggesting that surface-active species, which might be already present in the lab water or lipid samples, could further be introduced through the added salt.

Binding Mechanisms of Amyloid-like Peptides to Lipid Bilayers and Effects of Divalent Cations

In several neurodegenerative diseases, cell toxicity can emerge from damage produced by amyloid aggregates to lipid membranes. The details accounting for this damage are poorly understood including how individual amyloid peptides interact with phospholipid membranes before aggregation. Here, we use all-atom molecular dynamics simulations to investigate the molecular mechanisms accounting for amyloid-membrane interactions and the role played by calcium ions in this interaction. Model peptides known to self-assemble into amyloid fibrils and bilayer made from zwitterionic and anionic lipids are used in this study. We find that both electrostatic and hydrophobic interactions contribute to peptide-bilayer binding. In particular, the attraction of peptides to lipid bilayers is dominated by electrostatic interactions between positive residues and negative phosphate moieties of lipid head groups. This attraction is stronger for anionic bilayers than for zwitterionic ones. Hydrophobicity drives the burial of nonpolar residues into the interior of the bilayer producing strong binding in our simulations. Moreover, we observe that the attraction of peptides to the bilayer is significantly reduced in the presence of calcium ions. This is due to the binding of calcium ions to negative phosphate moieties of lipid head groups, which leaves phospholipid bilayers with a net positive charge. Strong binding of the peptide to the membrane occurs less frequently in the presence of calcium ions and involves the formation of a "Ca²⁺ bridge".

Transmembrane potential of physiologically relevant model membranes: Effects of membrane asymmetry

Transmembrane potential difference (V_m) plays important roles in regulating various biological processes. At the macro level, V_m can be experimentally measured or calculated using the Nernst or Goldman-Hodgkin-Katz equation. However, the atomic details responsible for its generation and impact on protein and lipid dynamics still need to be further elucidated. In this work, we performed a series of all-atom molecular dynamics (MD) simulations of symmetric model membranes of various lipid compositions and cation contents to evaluate the relationship between membrane asymmetry and V_m. Specifically, we studied the impact of the asymmetric distribution of POPS (1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-l-serine), PIP2 (phosphatidylinositol 4,5-bisphosphate), as well as Na⁺ and K⁺ on V_m using atomically detailed MD simulations of symmetric model membranes. The results suggest that, for an asymmetric POPC-POPC/POPS bilayer in the presence of NaCl, the presence of the monovalent anionic lipid POPS in the inner leaflet polarizes the membrane (ΔV_m < 0). Intriguingly, replacing a third of the POPS lipids by the polyvalent anionic signaling lipid PIP2 counteracts this effect, resulting in a smaller negative membrane potential. We also found that replacing Na⁺ ions in the inner region by K⁺ depolarizes the membrane (ΔV_m > 0). These divergent effects arise from variations in the strength of cation-lipid interactions and are correlated with changes in lipid chain order and head-group orientation.

Budding and Fission of Nanovesicles Induced by Membrane Adsorption of Small Solutes

Membrane budding and fission are essential cellular processes that produce new membrane compartments during cell and organelle division, for intracellular vesicle trafficking as well as during endo- and exocytosis. Such morphological transformations have also been observed for giant lipid vesicles with a size of many micrometers. Here, we report budding and fission processes of lipid nanovesicles with a size below 50 nm. We use coarse-grained molecular dynamics simulations, by which we can visualize the morphological transformations of individual vesicles. The budding and fission processes are induced by low concentrations of small solutes that absorb onto the outer leaflets of the vesicle membranes. In addition to the solute concentration, we identify the solvent conditions as a second key parameter for these processes. For good solvent conditions, the budding of a nanovesicle can be controlled by reducing the vesicle volume for constant solute concentration or by increasing the solute concentration for constant vesicle volume. After the budding process is completed, the budded vesicle consists of two membrane subcompartments which are connected by a closed membrane neck. The budding process is reversible as we demonstrate explicitly by reopening the closed neck. For poor solvent conditions, on the other hand, we observe two unexpected morphological transformations of nanovesicles. Close to the binodal line, at which the aqueous solution undergoes phase separation, the vesicle exhibits recurrent shape changes with closed and open membrane necks, reminiscent of flickering fusion pores (kiss-and-run) as observed for synaptic vesicles. As we approach the binodal line even closer, the recurrent shape changes are truncated by the fission of the membrane neck which leads to the division of the nanovesicle into two daughter vesicles. In this way, our simulations reveal a nanoscale mechanism for the budding and fission of nanovesicles, a mechanism that arises from the interplay between membrane elasticity and solute-mediated membrane adhesion.

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.

A practical guide to biologically relevant molecular simulations with charge scaling for electronic polarization

Molecular simulations can elucidate atomistic-level mechanisms of key biological processes, which are often hardly accessible to experiment. However, the results of the simulations can only be as trustworthy as the underlying simulation model. In many of these processes, interactions between charged moieties play a critical role. Current empirical force fields tend to overestimate such interactions, often in a dramatic way, when polyvalent ions are involved. The source of this shortcoming is the missing electronic polarization in these models. Given the importance of such biomolecular systems, there is great interest in fixing this deficiency in a computationally inexpensive way without employing explicitly polarizable force fields. Here, we review the electronic continuum correction approach, which accounts for electronic polarization in a mean-field way, focusing on its charge scaling variant. We show that by pragmatically scaling only the charged molecular groups, we qualitatively improve the charge-charge interactions without extra computational costs and benefit from decades of force field development on biomolecular force fields.

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.

Role of Ion-Phospholipid Interactions in Zwitterionic Phospholipid Bilayer Ion Permeation

Despite the central role of Na⁺ and K⁺ in physiological processes, it is still unclear whether they interact or alter the physical properties of simple zwitterionic phospholipid bilayers at physiologically relevant concentrations. Here we report a difference in membrane permeability between Na⁺ and K⁺, as measured with electrical impedance spectroscopy and tethered bilayer lipid membranes. We reveal that the differences in membrane permeability originate from distinct ion coordination by carbonyl oxygens at the phospholipid-water interface, altering the propensity for bilayer pore formation. Molecular dynamics simulations showed differences in the coordination of Na⁺ and K⁺ at the phospholipid-water interface of zwitterionic phospholipid bilayers. The ability of Na⁺ to conscript more phospholipids with a greater number of coordinating interactions causes a higher localized energy barrier for pore formation. These results provide evidence that ion-specific interactions at the phospholipid-water interface can modulate the physical properties of zwitterionic phospholipid bilayers.

Influence of Calcium Binding on Conformations and Motions of Anionic Polyamino Acids. Effect of Side Chain Length

Investigation of the effect of CaCl2 salt on conformations of two anionic poly(amino acids) with different side chain lengths, poly-(α-l glutamic acid) (PGA) and poly-(α-l aspartic acid) (PASA), was performed by atomistic molecular dynamics (MD) simulations. The simulations were performed using both unbiased MD and the Hamiltonian replica exchange (HRE) method. The results show that at low CaCl2 concentration adsorption of Ca2+ ions lead to a significant chain size reduction for both PGA and PASA. With the increase in concentration, the chains sizes partially recover due to electrostatic repulsion between the adsorbed Ca2+ ions. Here, the side chain length becomes important. Due to the longer side chain and its ability to distance the charged groups with adsorbed ions from both each other and the backbone, PGA remains longer in the collapsed state as the CaCl2 concentration is increased. The analysis of the distribution of the mineral ions suggests that both poly(amino acids) should induce the formation of mineral with the same structure of the crystal cell.

Charge-dependent interactions of monomeric and filamentous actin with lipid bilayers

The cytoskeletal protein actin polymerizes into filaments that are essential for the mechanical stability of mammalian cells. In vitro experiments showed that direct interactions between actin filaments and lipid bilayers are possible and that the net charge of the bilayer as well as the presence of divalent ions in the buffer play an important role. In vivo, colocalization of actin filaments and divalent ions are suppressed, and cells rely on linker proteins to connect the plasma membrane to the actin network. Little is known, however, about why this is the case and what microscopic interactions are important. A deeper understanding is highly beneficial, first, to obtain understanding in the biological design of cells and, second, as a possible basis for the building of artificial cortices for the stabilization of synthetic cells. Here, we report the results of coarse-grained molecular dynamics simulations of monomeric and filamentous actin in the vicinity of differently charged lipid bilayers. We observe that charges on the lipid head groups strongly determine the ability of actin to adsorb to the bilayer. The inclusion of divalent ions leads to a reversal of the binding affinity. Our in silico results are validated experimentally by reconstitution assays with actin on lipid bilayer membranes and provide a molecular-level understanding of the actin-membrane interaction.

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.

Concurrent Compression of Phospholipid Membranes by Calcium and Cholesterol

Regulation of cell metabolism, membrane fusion, association of proteins with cellular membranes, and cellular signaling altogether would not be possible without Ca²⁺ ions. The distribution of calcium within the cell is uneven with the negatively charged inner leaflet of the plasma membrane being one of the primary targets of its accumulation. Therefore, we decided to map the influence of Ca²⁺ on the properties of lipid bilayers closely resembling natural lipid membranes. We combined fluorescence spectroscopy (analysis of time-resolved emission spectra of Laurdan probe and derived parameters: integrated relaxation time related to local lipid mobility, and total emission shift reflecting membrane polarity and hydration) with molecular dynamics simulations to determine the effect of the increasing CaCl₂ concentration on model lipid membranes containing POPC, POPS, and cholesterol. On top of that, the impact of calcium on the plasma membranes isolated from HEK293 cells was investigated using the steady-state fluorescence of Laurdan. We found that calcium increases rigidity of all the model lipid membranes used, elevates their thickness, increases lipid packing and ordering, and impedes the local lipid mobility. All these effects were to a great extent similar to those elicited by cholesterol. However, the changes of the membrane properties induced by calcium and cholesterol seem largely independent from each other. At sufficiently high concentrations of calcium or cholesterol, the steric effects hindered a further alteration of membrane organization, i.e., the compressibility limit of membrane structures was reached. We found no indication for mutual interaction between Ca²⁺ and cholesterol, nor competition of Ca²⁺ ions and hydroxyl groups of cholesterol for binding to phospholipids. Fluorescence measurements indicated that Ca²⁺ adsorption decreases mobility within the carbonyl region of model bilayers more efficiently than monovalent ions do (Ca²⁺ ≫ Li⁺ > Na⁺ > K⁺ > Cs⁺). The effects of calcium ions were to a great extent mitigated in the plasma membranes isolated from HEK293 cells when compared to the model lipid membranes. Noticeably, the plasma membranes showed remarkably higher resistance toward rigidification induced by calcium ions even when compared with the model membranes containing cholesterol.

Competing for the same space: protons and alkali ions at the interface of phospholipid bilayers

Maintaining gradients of solvated protons and alkali metal ions such as Na⁺ and K⁺ across membranes is critical for cellular function. Over the last few decades, both the interactions of protons and alkali metal ions with phospholipid membranes have been studied extensively and the reported interactions of these ions with phospholipid headgroups are very similar, yet few studies have investigated the potential interdependence between proton and alkali metal ion binding at the water-lipid interface. In this short review, we discuss the similarities between the proton-membrane and alkali ion-membrane interactions. Such interactions include cation attraction to the phosphate and carbonyl oxygens of the phospholipid headgroups that form strong lipid-ion and lipid-ion-water complexes. We also propose potential mechanisms that may modulate the affinities of these cationic species to the water-phospholipid interfacial oxygen moieties. This review aims to highlight the potential interdependence between protons and alkali metal ions at the membrane surface and encourage a more nuanced understanding of the complex nature of these biologically relevant processes.

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.

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.