Simulating Monovalent and Divalent Ions in Aqueous Solution Using a Drude Polarizable Force Field.
Ontology highlight
ABSTRACT:
SUBMITTER: Yu H
PROVIDER: S-EPMC2838399 | biostudies-other | 2010
REPOSITORIES: biostudies-other
Publications
Simulating Monovalent and Divalent Ions in Aqueous Solution Using a Drude Polarizable Force Field.
Yu Haibo H Whitfield Troy W TW Harder Edward E Lamoureux Guillaume G Vorobyov Igor I Anisimov Victor M VM Mackerell Alexander D AD Roux Benoît B
Journal of chemical theory and computation 20100101 3
An accurate representation of ion solvation in aqueous solution is critical for meaningful computer simulations of a broad range of physical and biological processes. Polarizable models based on classical Drude oscillators are introduced and parametrized for a large set of monoatomic ions including cations of the alkali metals (Li(+), Na(+), K(+), Rb(+) and Cs(+)) and alkaline earth elements (Mg(2+), Ca(2+), Sr(2+) and Ba(2+)) along with Zn(2+) and halide anions (F(-), Cl(-), Br(-) and I(-)). Th ...[more]
Similar Datasets
Project description:Development of accurate force field parameters for molecular ions in the context of a polarizable energy function based on the classical Drude oscillator is a crucial step toward an accurate polarizable model for modeling and simulations of biological macromolecules. Toward this goal we have undertaken a hierarchical approach in which force field parameter optimization is initially performed for small molecules for which experimental data exists that serve as building blocks of macromolecular systems. Small molecules representative of the ionic moieties of biological macromolecules include the cationic ammonium and methyl substituted ammonium derivatives, imidazolium, guanidinium and methylguanidinium, and the anionic acetate, phenolate, and alkanethiolates. In the present work, parameters for molecular ions in the context of the Drude polarizable force field are optimized and compared to results from the nonpolarizable additive CHARMM general force field (CGenFF). Electrostatic and Lennard-Jones parameters for the model compounds are developed in the context of the polarizable SWM4-NDP water model, with emphasis on assuring that the hydration free energies are consistent with previously reported parameters for atomic ions. The final parameters are shown to be in good agreement with the selected quantum mechanical (QM) and experimental target data. Analysis of the structure of water around the ions reveals substantial differences between the Drude and additive force fields indicating the important role of polarization in dictating the molecular details of aqueous solvation. The presented parameters represent the foundation for the charged functionalities in future generations of the Drude polarizable force field for biological macromolecules as well as for drug-like molecules.
| S-EPMC5975207 | biostudies-literature
Project description:Recently we presented a first-generation all-atom Drude polarizable force field for DNA based on the classical Drude oscillator model, focusing on optimization of key dihedral angles followed by extensive validation of the force field parameters. Presently, we describe the procedure for balancing the electrostatic interactions between ions, water, and DNA as required for development of the Drude force field for DNA. The proper balance of these interactions is shown to impact DNA stability and subtler conformational properties, including the conformational equilibrium between the BI and BII states, and the A and B forms of DNA. The parametrization efforts were simultaneously guided by gas-phase quantum mechanics (QM) data on small model compounds and condensed-phase experimental data on the hydration and osmotic properties of biologically relevant ions and their solutions, as well as theoretical predictions for ionic distribution around DNA oligomer. In addition, fine-tuning of the internal base parameters was performed to obtain the final DNA model. Notably, the Drude model is shown to more accurately reproduce counterion condensation theory predictions of DNA charge neutralization by the condensed ions as compared to the CHARMM36 additive DNA force field, indicating an improved physical description of the forces dictating the ionic solvation of DNA due to the explicit treatment of electronic polarizability. In combination with the polarizable DNA force field, the availability of Drude polarizable parameters for proteins, lipids, and carbohydrates will allow for simulation studies of heterogeneous biological systems.
| S-EPMC4064693 | biostudies-literature
Project description:An empirical all-atom CHARMM polarizable force filed for aldopentofuranoses and methyl-aldopentofuranosides based on the classical Drude oscillator is presented. A single electrostatic model is developed for eight different diastereoisomers of aldopentofuranoses by optimizing the existing electrostatic and bonded parameters as transferred from ethers, alcohols, and hexopyranoses to reproduce quantum mechanical (QM) dipole moments, furanose-water interaction energies and conformational energies. Optimization of selected electrostatic and dihedral parameters was performed to generate a model for methyl-aldopentofuranosides. Accuracy of the model was tested by reproducing experimental data for crystal intramolecular geometries and lattice unit cell parameters, aqueous phase densities, and ring pucker and exocyclic rotamer populations as obtained from NMR experiments. In most cases the model is found to reproduce both QM data and experimental observables in an excellent manner, whereas for the remainder the level of agreement is in the satisfactory regimen. In aqueous phase simulations the monosaccharides have significantly enhanced dipoles as compared to the gas phase. The final model from this study is transferrable for future studies on carbohydrates and can be used with the existing CHARMM Drude polarizable force field for biomolecules.
| S-EPMC4483154 | biostudies-literature
Project description:Explicit treatment of electronic polarizability in empirical force fields (FFs) represents an extension over a traditional additive or pairwise FF and provides a more realistic model of the variations in electronic structure in condensed phase, macromolecular simulations. To facilitate utilization of the polarizable FF based on the classical Drude oscillator model, Drude Prepper has been developed in CHARMM-GUI. Drude Prepper ingests additive CHARMM protein structures file (PSF) and pre-equilibrated coordinates in CHARMM, PDB, or NAMD format, from which the molecular components of the system are identified. These include all residues and patches connecting those residues along with water, ions, and other solute molecules. This information is then used to construct the Drude FF-based PSF using molecular generation capabilities in CHARMM, followed by minimization and equilibration. In addition, inputs are generated for molecular dynamics (MD) simulations using CHARMM, GROMACS, NAMD, and OpenMM. Validation of the Drude Prepper protocol and inputs is performed through conversion and MD simulations of various heterogeneous systems that include proteins, nucleic acids, lipids, polysaccharides, and atomic ions using the aforementioned simulation packages. Stable simulations are obtained in all studied systems, including 5 μs simulation of ubiquitin, verifying the integrity of the generated Drude PSFs. In addition, the ability of the Drude FF to model variations in electronic structure is shown through dipole moment analysis in selected systems. The capabilities and availability of Drude Prepper in CHARMM-GUI is anticipated to greatly facilitate the application of the Drude FF to a range of condensed phase, macromolecular systems.
| S-EPMC8741736 | biostudies-literature
Project description:Small metal ions play critical roles in numerous biological processes. Of particular interest is how metalloenzymes are allosterically regulated by the binding of specific ions. Understanding how ion binding affects these biological processes requires atomic models that accurately treat the microscopic interactions with the protein ligands. Theoretical approaches at different levels of sophistication can contribute to a deeper understanding of these systems, although computational models must strike a balance between accuracy and efficiency in order to enable long molecular dynamics simulations. In this study, we present a systematic effort to optimize the parameters of a polarizable force field based on classical Drude oscillators to accurately represent the interactions between ions (K(+), Na(+), Ca(2+), and Cl(-)) and coordinating amino-acid residues for a set of 30 biologically important proteins. By combining ab initio calculations and experimental thermodynamic data, we derive a polarizable force field that is consistent with a wide range of properties, including the geometries and interaction energies of gas-phase ion/protein-like model compound clusters, and the experimental solvation free-energies of the cations in liquids. The resulting models display significant improvements relative to the fixed-atomic-charge additive CHARMM C36 force field, particularly in their ability to reproduce the many-body electrostatic nonadditivity effects estimated from ab initio calculations. The analysis clarifies the fundamental limitations of the pairwise additivity assumption inherent in classical fixed-charge force fields, and shows its dramatic failures in the case of Ca(2+) binding sites. These optimized polarizable models, amenable to computationally efficient large-scale MD simulations, set a firm foundation and offer a powerful avenue to study the roles of the ions in soluble and membrane transport proteins.
| S-EPMC4516320 | biostudies-literature
Project description:Polarizable force fields, in particular, the Drude polarizable force field (Drude FF), may hold the key to more accurately modeling biomolecules with molecular dynamics simulations by explicitly accounting for atomic polarizability. Previous work has shown promising results in simulating duplex nucleic acids and protein structures with excellent agreement with experimental values. However, benchmarking the Drude polarizable force field with highly flexible, single-stranded structures has yet to be achieved. In this work, the r(GACC) tetranucleotide is simulated over a multimicrosecond time scale, starting with various different initial conformations. Despite the starting conformation, including starting from the expected dominant A-form major conformation, the experimental structural distribution is not matched. In fact, the major NMR conformation is never resampled. Instead, the r(GACC) tetranucleotide becomes stabilized in anomalous structures that are inconsistent with the NMR data and that favor base-pairing and electrostatic interactions over base stacking. These structures are maintained for lengthy time scales (>1 μs) themselves, suggesting a misbalance of forces in the Drude polarizable force field itself. This model system is suggestive of the fact that currently the Drude polarizable force field does not appear to produce the sensitive balance of forces required to accurately model other single-stranded or noncanonical RNA structures.
| S-EPMC10131218 | biostudies-literature
Project description:RNA molecules are highly dynamic and capable of adopting a wide range of complex, folded structures. The factors driving the folding and dynamics of these structures are dependent on a balance of base pairing, hydration, base stacking, ion interactions, and the conformational sampling of the 2'-hydroxyl group in the ribose sugar. The representation of these features is a challenge for empirical force fields used in molecular dynamics simulations. Toward meeting this challenge, the inclusion of explicit electronic polarization is important in accurately modeling RNA structure. In this work, we present a polarizable force field for RNA based on the classical Drude oscillator model, which represents electronic degrees of freedom via negatively charged particles attached to their parent atoms by harmonic springs. Beginning with parametrization against quantum mechanical base stacking interaction energy and conformational energy data, we have extended the Drude-2017 nucleic acid force field to include RNA. The conformational sampling of a range of RNA sequences were used to validate the force field, including canonical A-form RNA duplexes, stem-loops, and complex tertiary folds that bind multiple Mg2+ ions. Overall, the Drude-2017 RNA force field reproduces important properties of these structures, including the conformational sampling of the 2'-hydroxyl and key interactions with Mg2+ ions. © 2018 Wiley Periodicals, Inc.
| S-EPMC6284239 | biostudies-literature
Project description:The CHARMM Drude-2013 polarizable force field (FF) was developed to include the explicit treatment of induced electronic polarizability, resulting in a more accurate description of the electrostatic interactions in molecular dynamics (MD) simulations. While the Drude-2013 protein FF has shown success in improving the folding properties of α-helical peptides and to reproduce experimental observables in simulations up to 1 μs, some limitations were noted regarding the stability of β-sheet structures in simulations longer than 100 ns as well as larger deviations from crystal structures in simulations of a number of proteins compared to the additive CHARMM36 protein FF. The origin of the instability has been identified and appears to be primarily due to overestimated atomic polarizabilities and induced dipole-dipole interactions on the Cβ, Cγ, and Cδ side chain atoms. To resolve this and other issues, a number of aspects of the model were revisited, resulting in Drude-2019 protein FF. Backbone parameters were optimized targeting the conformational properties of the (Ala)5 peptide in solution along with gas phase properties of the alanine dipeptide. Dipeptides that contain N-acetylated and N'-methylamidated termini, excluding Gly, Pro, and Ala, were used as models to optimize the atomic polarizabilities and Thole screening factors on selected Cβ, Cγ, and Cδ carbons by targeting quantum mechanical (QM) dipole moments and molecular polarizabilities. In addition, to obtain better conformational properties, side chain χ1 and χ2 dihedral parameters were optimized targeting QM data for the respective side chain dipeptide conformations as well as Protein Data Bank survey data based on the χ1, χ2 sampling from Hamiltonian replica-exchange MD simulations of (Ala)4-X-(Ala)4 in solution, where X is the amino acid of interest. Further improvements include optimizing nonbonded interactions between charged residues to reproduce QM interaction energies of the charged-protein model compounds and experimental osmotic pressures. Validation of the optimized Drude protein FF includes MD simulations of a collection of peptides and proteins including β-sheet structures, as well as transmembrane ion channels. Results showed that the updated Drude-2019 protein FF yields smaller overall root-mean-square differences of proteins as compared to the additive CHARMM36m and Drude-2013 FFs as well as similar or improved agreement with experimental NMR properties, allowing for long time scale simulation studies of proteins and more complex biomolecular systems in conjunction with the remainder of the Drude polarizable FF.
| S-EPMC7306265 | biostudies-literature
Project description:We present an extension of the CHARMM Drude polarizable force field to enable modeling of polysaccharides containing pyranose and furanose monosaccharides. The new force field parameters encompass 1↔2, 1→3, 1→4, and 1→6 pyranose-furanose linkages, 2→1 and 2→6 furanose-furanose linkages, 2→2, 2→3, and 2→4 furanose-pyranose, and 1↔1, 1→2, 1→3, 1→4, and 1→6 pyranose-pyranose linkages. For the glycosidic linkages, both simple model compounds and the full disaccharides with methylation at the reducing end were used for parameter optimization. The model compounds were chosen to be monomers or glycosidic-linked dimers of tetrahydropyran (THP) and tetrahydrofuran (THF). Target data for optimization included one- and two-dimensional potential energy scans of ω and the Φ/Ψ glycosidic dihedral angles in the model compounds and full disaccharides computed by quantum mechanical (QM) RIMP2/cc-pVQZ single point energies on MP2/6-31G(d) optimized structures. Also included in the target data are extensive sets of QM gas phase monohydrate water-saccharide interactions, dipole moments, and molecular polarizabilities for both model compounds and full disaccharides. The resulting polarizable model is shown to be in good agreement with a range of QM data, offering a significant improvement over the additive CHARMM36 carbohydrate force field, as well as experimental data including crystal structures and conformational properties of disaccharides and a trisaccharide in aqueous solution.
| S-EPMC5997548 | biostudies-literature
Project description:In the present study we report on interactions of and competition between monovalent ions for two DNA sequences in MD simulations. Efforts included the development and validation of parameters for interactions among the first-group monovalent cations, Li(+), Na(+), K(+), and Rb(+), and DNA in the Drude polarizable and additive CHARMM36 force fields (FF). The optimization process targeted gas-phase QM interaction energies of various model compounds with ions and osmotic pressures of bulk electrolyte solutions of chemically relevant ions. The optimized ionic parameters are validated against counterion condensation theory and buffer exchange-atomic emission spectroscopy measurements providing quantitative data on the competitive association of different monovalent ions with DNA. Comparison between experimental and MD simulation results demonstrates that, compared to the additive CHARMM36 model, the Drude FF provides an improved description of the general features of the ionic atmosphere around DNA and leads to closer agreement with experiment on the ionic competition within the ion atmosphere. Results indicate the importance of extended simulation systems on the order of 25 Å beyond the DNA surface to obtain proper convergence of ion distributions.
| S-EPMC4378841 | biostudies-literature