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.

• Publikační typ
• časopisecké články MeSH

Glycosaminoglycans (GAGs) are negatively charged polysaccharides found on cell surfaces, where they regulate transport pathways of foreign molecules toward the cell. The structural and functional diversity of GAGs is largely attributed to varied sulfation patterns along the polymer chains, which makes understanding their molecular recognition mechanisms crucial. Molecular dynamics (MD) simulations, thanks to their unmatched microscopic resolution, have the potential to be a reference tool for exploring the patterns responsible for biologically relevant interactions. However, the capability of molecular dynamics force fields used in biosimulations to accurately capture sulfation-specific interactions is not well established, partly due to the intrinsic properties of GAGs that pose challenges for most experimental techniques. In this work, we evaluate the performance of molecular dynamics force fields for sulfated GAGs by studying ion pairing of Ca2+ to sulfated moieties─N-methylsulfamate and methylsulfate─that resemble N- and O-sulfation found in GAGs, respectively. We tested available nonpolarizable (CHARMM36 and GLYCAM06) and explicitly polarizable (Drude and AMOEBA) force fields, and derived new implicitly polarizable models through charge scaling (prosECCo75 and GLYCAM-ECC75) that are consistent with our developed "charge-scaling" framework. The calcium-sulfamate/sulfate interaction free energy profiles obtained with the tested force fields were compared against reference ab initio molecular dynamics (AIMD) simulations, which serve as a robust alternative to experiments. AIMD simulations indicate that the preferential Ca2+ binding mode to sulfated GAG groups is solvent-shared pairing. Only our scaled-charge models agree satisfactorily with the AIMD data, while all other force fields exhibit poorer agreement, sometimes even qualitatively. Surprisingly, even explicitly polarizable force fields display a notable disagreement with the AIMD data, likely attributed to difficulties in their optimization and possible inherent limitations in depicting high-charge-density ion interactions accurately. Finally, the underperforming force fields lead to unrealistic aggregation of sulfated saccharides, which qualitatively disagrees with our understanding of the soft glycocalyx environment. Our results highlight the importance of accurately treating electronic polarization in MD simulations of sulfated GAGs and caution against over-reliance on currently available models without thorough validation and optimization.

• MeSH
• glykosaminoglykany * chemie   MeSH
• kyseliny sulfonové chemie   MeSH
• simulace molekulární dynamiky * MeSH
• sírany * chemie   MeSH
• statická elektřina * MeSH
• vápník chemie   MeSH
• Publikační typ
• časopisecké články MeSH
• Názvy látek
• glykosaminoglykany * MeSH
• kyseliny sulfonové MeSH
• sírany * MeSH
• sulfamic acid MeSH Prohlížeč
• vápník 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
• peptidy chemie   MeSH
• proteiny chemie   MeSH
• simulace molekulární dynamiky * MeSH
• statická elektřina * MeSH
• voda chemie   MeSH
• Publikační typ
• časopisecké články MeSH
• Názvy látek
• peptidy MeSH
• proteiny MeSH
• voda 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
• fosfatidylcholiny chemie   MeSH
• fosfolipidy * chemie   MeSH
• ionty MeSH
• lipidové dvojvrstvy * chemie   MeSH
• simulace molekulární dynamiky MeSH
• vápník chemie   MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• Názvy látek
• fosfatidylcholiny MeSH
• fosfolipidy * MeSH
• ionty MeSH
• lipidové dvojvrstvy * MeSH
• vápník MeSH

The functional properties of some biological ion channels and membrane transport proteins are proposed to exploit anion-hydrophobic interactions. Here, we investigate a chloride-pumping rhodopsin as an example of a membrane protein known to contain a defined anion binding site composed predominantly of hydrophobic residues. Using molecular dynamics simulations, we explore Cl- binding to this hydrophobic site and compare the dynamics arising when electronic polarization is neglected (CHARMM36 [c36] fixed-charge force field), included implicitly (via the prosECCo force field), or included explicitly (through the polarizable force field, AMOEBA). Free energy landscapes of Cl- moving out of the binding site and into bulk solution demonstrate that the inclusion of polarization results in stronger ion binding and a second metastable binding site in chloride-pumping rhodopsin. Simulations focused on this hydrophobic binding site also indicate longer binding durations and closer ion proximity when polarization is included. Furthermore, simulations reveal that Cl- within this binding site interacts with an adjacent loop to facilitate rebinding events that are not observed when polarization is neglected. These results demonstrate how the inclusion of polarization can influence the behavior of anions within protein binding sites and can yield results comparable with more accurate and computationally demanding methods.

• MeSH
• anionty MeSH
• chloridy * chemie   MeSH
• elektronika MeSH
• rodopsin * MeSH
• simulace molekulární dynamiky MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• Názvy látek
• anionty MeSH
• chloridy * MeSH
• rodopsin * 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.

• Klíčová slova
• carboxyls, counterions, force fields, ions, molecular dynamics, peptides and proteins,
• Publikační typ
• časopisecké články MeSH