Salt bridges are ionic interactions that are of great importance in protein recognition. However, their structural description using X-ray crystallography or NMR may be inconclusive. Classical molecular dynamics (MD) used for the interpretation neglects electronic polarization, which results in artifactual overbinding. Here, we resolve the problem via charge scaling, which accounts for electronic polarization in a mean-field way. We study three salt bridges in insulin analogue. New NMR ensembles are generated via NOE-restrained MD using ff19SB and CHARMM36m force fields and the scaled-charge prosECCo75. Tens of μs of unrestrained MD show in a statistically converged manner that ff19SB induces a non-native salt bridge. This behavior is quantified via umbrella sampling of salt bridge dissociation, which indicates a rather high strength of up to 4 and 5 kcal mol-1 for CHARMM36m and ff19SB, respectively. In contrast, prosECCo75 gives a biologically reasonable dissociation barrier of 1 kcal mol-1. Our results indicate that a physically justified description of charge-charge interactions within a nonpolarizable MD framework reliably describes aqueous biomolecular systems.

• Publication type
• Journal Article MeSH

This study employs molecular dynamics (MD) simulations to investigate the adsorption and aggregation behavior of simple polyarginine cell-penetrating peptides (CPPs), specifically modeled as R9 peptides, at zwitterionic phosphocholine POPC membranes under varying ionic strengths of two peptide concentrations and two concentrations of NaCl and CaCl2. The results reveal an intriguing phenomenon of R9 aggregation at the membrane, which is dependent on the ionic strength, indicating a salting-out effect. As the peptide concentration and ionic strength increase, peptide aggregation also increases, with aggregate lifetimes and sizes showing a corresponding rise, accompanied by the total decrease of adsorbed peptides at the membrane surface. Notably, in high ionic strength environments, large R9 aggregates, such as octamers, are also observed occasionally. The salting-out, typically uncommon for short positively charged peptides, is attributed to the unique properties of arginine amino acid, specifically by its side chain containing amphiphilic guanidinium (Gdm+) ion which makes both intermolecular hydrophobic like-charge Gdm+ - Gdm+ and salt-bridge Gdm+ - C-terminus interactions, where the former are increased with the ionic strength, and the latter decreased due to electrostatic screening. The aggregation behavior of R9 peptides at membranes can also be linked to their CPP translocation properties, suggesting that aggregation may aid in translocation across cellular membranes.

• Keywords
• Ionic strength, Molecular dynamics simulations, Peptide aggregation, Phosphocholine lipid bilayers, Polyarginines, Salting-out,
• Publication type
• Journal Article 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

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

Ethylenediaminetetraacetic acid (EDTA) is frequently used in lipid experiments to remove redundant ions, such as Ca2+, from the sample solution. In this work, combining molecular dynamics (MD) simulations and Langmuir monolayer experiments, we show that on top of the expected Ca2+ depletion, EDTA anions themselves bind to phosphatidylcholine (PC) monolayers. This binding, originating from EDTA interaction with choline groups of PC lipids, leads to the adsorption of EDTA anions at the monolayer surface and concentration-dependent changes in surface pressure as measured by monolayer experiments and explained by MD simulations. This surprising observation emphasizes that lipid experiments carried out using EDTA-containing solutions, especially of high concentrations, must be interpreted very carefully due to potential interfering interactions of EDTA with lipids and other biomolecules involved in the experiment, e.g., cationic peptides, that may alter membrane-binding affinities of studied compounds.

• MeSH
• Edetic Acid MeSH
• Phosphatidylcholines * chemistry   MeSH
• Ions MeSH
• Membranes, Artificial * MeSH
• Molecular Dynamics Simulation MeSH
• Publication type
• Journal Article MeSH
• Research Support, Non-U.S. Gov't MeSH
• Names of Substances
• Edetic Acid MeSH
• Phosphatidylcholines * MeSH
• Ions MeSH
• Membranes, Artificial * MeSH

In this work, the influence of membrane curvature on the Ca2+ binding to phospholipid bilayers is investigated by means of molecular dynamics simulations. In particular, we compared Ca2+ binding to flat, elastically buckled, or uniformly bent zwitterionic and anionic phospholipid bilayers. We demonstrate that Ca2+ ions bind preferably to the concave membrane surfaces in both types of bilayers. We also show that the membrane curvature leads to pronounced changes in Ca2+ binding including differences in free ion concentrations, lipid coordination distributions, and the patterns of ion binding to different chemical groups of lipids. Moreover, these effects differ substantially for the concave and convex membrane monolayers. Comparison between force fields with either full or scaled charges indicates that charge scaling results in reduction of the Ca2+ binding to curved phosphatidylserine bilayers, while for phosphatidylcholine membranes, calcium binds only weakly for both force fields.

• MeSH
• Phosphatidylcholines chemistry   MeSH
• Phospholipids * chemistry   MeSH
• Ions MeSH
• Lipid Bilayers * chemistry   MeSH
• Molecular Dynamics Simulation MeSH
• Calcium chemistry   MeSH
• Publication type
• Journal Article MeSH
• Research Support, Non-U.S. Gov't MeSH
• Names of Substances
• Phosphatidylcholines MeSH
• Phospholipids * MeSH
• Ions MeSH
• Lipid Bilayers * MeSH
• Calcium 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

Lipid monolayers provide our lungs and eyes their functionality and serve as proxy systems in biomembrane research. Therefore, lipid monolayers have been studied intensively including using molecular dynamics simulations, which are able to probe their lateral structure and interactions with, e.g., pharmaceuticals or nanoparticles. However, such simulations have struggled in describing the forces at the air-water interface. Particularly, the surface tension of water and long-range van der Waals interactions have been considered critical, but their importance in monolayer simulations has been evaluated only separately. Here, we combine the recent C36/LJ-PME lipid force field that includes long-range van der Waals forces with water models that reproduce experimental surface tensions to elucidate the importance of these contributions in monolayer simulations. Our results suggest that a water model with correct surface tension is necessary to reproduce experimental surface pressure-area isotherms and monolayer phase behavior. The latter includes the liquid expanded and liquid condensed phases, their coexistence, and the opening of pores at the correct area per lipid upon expansion. Despite these improvements of the C36/LJ-PME with certain water models, the standard cutoff-based CHARMM36 lipid model with the 4-point OPC water model still provides the best agreement with experiments. Our results emphasize the importance of using high-quality water models in applications and parameter development in molecular dynamics simulations of biomolecules.

• MeSH
• Lipids MeSH
• Surface Tension MeSH
• Molecular Dynamics Simulation * MeSH
• Water * chemistry   MeSH
• Publication type
• Journal Article MeSH
• Names of Substances
• Lipids MeSH
• Water * 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

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.

• Keywords
• calcium phosphate, functional groups, molecular dynamics, nucleation mechanism, parylene C, polymer surface,
• Publication type
• Journal Article MeSH