Understanding the molecular mechanisms of pore formation is crucial for elucidating fundamental biological processes and developing therapeutic strategies, such as the design of drug delivery systems and antimicrobial agents. Although experimental methods can provide valuable information, they often lack the temporal and spatial resolution necessary to fully capture the dynamic stages of pore formation. In this study, we present two novel collective variables (CVs) designed to characterize membrane pore behavior, particularly its energetics, through molecular dynamics (MD) simulations. The first CV─termed Full-Path─effectively tracks both the nucleation and expansion phases of pore formation. The second CV─called Rapid─is tailored to accurately assess pore expansion in the limit of large pores, providing quick and reliable method for evaluating membrane line tension under various conditions. Our results clearly demonstrate that the line tension predictions from both our CVs are in excellent agreement. Moreover, these predictions align qualitatively with available experimental data. Specifically, they reflect higher line tension of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) membranes containing 1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-l-serine (POPS) lipids compared to pure POPC, the decrease in line tension of POPC vesicles as the 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoglycerol (POPG) content increases, and higher line tension when ionic concentration is increased. Notably, these experimental trends are accurately captured only by the all-atom CHARMM36 and prosECCo75 force fields. In contrast, the all-atom Slipids force field, along with the coarse-grained Martini 2.2, Martini 2.2 polarizable, and Martini 3 models, show varying degrees of agreement with experiments. Our developed CVs can be adapted to various MD simulation engines for studying pore formation, with potential implications in membrane biophysics. They are also applicable to simulations involving external agents, offering an efficient alternative to existing methodologies.

• MeSH
• Cell Membrane * chemistry   metabolism   MeSH
• Phosphatidylcholines chemistry   MeSH
• Lipid Bilayers * chemistry   MeSH
• Porosity MeSH
• Molecular Dynamics Simulation * MeSH
• Thermodynamics * MeSH
• Publication type
• Journal Article MeSH
• Names of Substances
• 1-palmitoyl-2-oleoylphosphatidylcholine MeSH Browser
• Phosphatidylcholines MeSH
• Lipid Bilayers * MeSH

prosECCo75 is an optimized force field effectively incorporating electronic polarization via charge scaling. It aims to enhance the accuracy of nominally nonpolarizable molecular dynamics simulations for interactions in biologically relevant systems involving water, ions, proteins, lipids, and saccharides. Recognizing the inherent limitations of nonpolarizable force fields in precisely modeling electrostatic interactions essential for various biological processes, we mitigate these shortcomings by accounting for electronic polarizability in a physically rigorous mean-field way that does not add to computational costs. With this scaling of (both integer and partial) charges within the CHARMM36 framework, prosECCo75 addresses overbinding artifacts. This improves agreement with experimental ion binding data across a broad spectrum of systems─lipid membranes, proteins (including peptides and amino acids), and saccharides─without compromising their biomolecular structures. prosECCo75 thus emerges as a computationally efficient tool providing enhanced accuracy and broader applicability in simulating the complex interplay of interactions between ions and biomolecules, pivotal for improving our understanding of many biological processes.

• MeSH
• Peptides chemistry   MeSH
• Proteins chemistry   MeSH
• Molecular Dynamics Simulation * MeSH
• Static Electricity * MeSH
• Water chemistry   MeSH
• Publication type
• Journal Article MeSH
• Names of Substances
• Peptides MeSH
• Proteins MeSH
• Water MeSH

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.

• MeSH
• Cell Membrane MeSH
• Membrane Lipids * MeSH
• Molecular Dynamics Simulation MeSH
• Machine Learning MeSH
• Artificial Intelligence * MeSH
• Publication type
• Journal Article MeSH
• Names of Substances
• Membrane Lipids * MeSH

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.

• MeSH
• Lipid Bilayers * chemistry   MeSH
• Solvents MeSH
• Molecular Dynamics Simulation * MeSH
• Publication type
• Journal Article MeSH
• Names of Substances
• Lipid Bilayers * MeSH
• Solvents MeSH

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.

• MeSH
• Cholesterol chemistry   MeSH
• Phosphatidylcholines * chemistry   MeSH
• Lipid Bilayers * chemistry   MeSH
• Molecular Dynamics Simulation MeSH
• Publication type
• Journal Article MeSH
• Names of Substances
• 1-palmitoyl-2-oleoylphosphatidylcholine MeSH Browser
• Cholesterol MeSH
• Phosphatidylcholines * MeSH
• Lipid Bilayers * MeSH

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 R9 and K9 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 Ca2+ salt concentrations, showing that the apparent strength of adsorption of R9 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 R9 and K9, we first demonstrated that the binding of R9 to POPC is tighter by almost 2 orders of magnitude compared to that of K9. Finally, upon the addition of an excess of either Na+ or Ca2+ ions with R9, the total fluorescence correlation signal is lost, which implies the unbinding of R9 from the PC bilayer, in agreement with our predictions from MD simulations.

• MeSH
• Adsorption MeSH
• Arginine MeSH
• Phosphatidylcholines chemistry   MeSH
• Phosphorylcholine MeSH
• Lecithins MeSH
• Lipid Bilayers * chemistry   MeSH
• Osmolar Concentration MeSH
• Cell-Penetrating Peptides * chemistry   MeSH
• Publication type
• Journal Article MeSH
• Research Support, Non-U.S. Gov't MeSH
• Names of Substances
• Arginine MeSH
• Phosphatidylcholines MeSH
• Phosphorylcholine MeSH
• Lecithins MeSH
• Lipid Bilayers * MeSH
• Cell-Penetrating Peptides * MeSH

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.

• Keywords
• carboxyls, counterions, force fields, ions, molecular dynamics, peptides and proteins,
• Publication type
• Journal Article MeSH

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.

• MeSH
• Models, Biological * MeSH
• Phospholipids chemistry   metabolism   MeSH
• Carboxylic Acids chemistry   metabolism   MeSH
• Humans MeSH
• Lipidomics methods   MeSH
• Membrane Lipids chemistry   metabolism   MeSH
• Membranes chemistry   metabolism   physiology   MeSH
• Computer Simulation MeSH
• Animals MeSH
• Check Tag
• Humans MeSH
• Animals MeSH
• Publication type
• Journal Article MeSH
• Research Support, Non-U.S. Gov't MeSH
• Names of Substances
• Phospholipids MeSH
• Carboxylic Acids MeSH
• Membrane Lipids MeSH

Understanding interactions of calcium with lipid membranes at the molecular level is of great importance in light of their involvement in calcium signaling, association of proteins with cellular membranes, and membrane fusion. We quantify these interactions in detail by employing a combination of spectroscopic methods with atomistic molecular dynamics simulations. Namely, time-resolved fluorescent spectroscopy of lipid vesicles and vibrational sum frequency spectroscopy of lipid monolayers are used to characterize local binding sites of calcium in zwitterionic and anionic model lipid assemblies, while dynamic light scattering and zeta potential measurements are employed for macroscopic characterization of lipid vesicles in calcium-containing environments. To gain additional atomic-level information, the experiments are complemented by molecular simulations that utilize an accurate force field for calcium ions with scaled charges effectively accounting for electronic polarization effects. We demonstrate that lipid membranes have substantial calcium-binding capacity, with several types of binding sites present. Significantly, the binding mode depends on calcium concentration with important implications for calcium buffering, synaptic plasticity, and protein-membrane association.

• MeSH
• Cell Membrane metabolism   MeSH
• Phospholipids chemistry   metabolism   MeSH
• Lipid Bilayers chemistry   metabolism   MeSH
• Liposomes chemistry   metabolism   MeSH
• Models, Molecular MeSH
• Molecular Dynamics Simulation MeSH
• Calcium metabolism   MeSH
• Calcium Signaling MeSH
• Binding Sites MeSH
• Publication type
• Journal Article MeSH
• Research Support, Non-U.S. Gov't MeSH
• Research Support, U.S. Gov't, Non-P.H.S. MeSH
• Names of Substances
• Phospholipids MeSH
• Lipid Bilayers MeSH
• Liposomes MeSH
• Calcium MeSH