Dataset Information

Optimized Magnesium Force Field Parameters for Biomolecular Simulations with Accurate Solvation, Ion-Binding, and Water-Exchange Properties.

ABSTRACT: Magnesium ions play an essential role in many vital processes. To correctly describe their interactions in molecular dynamics simulations, an accurate parametrization is crucial. Despite the importance and considerable scientific effort, current force fields based on the commonly used 12-6 Lennard-Jones interaction potential fail to reproduce a variety of experimental solution properties. In particular, no parametrization exists so far that simultaneously reproduces the solvation free energy and the distance to the water oxygens in the first hydration shell. Moreover, current Mg²⁺ force fields significantly underestimate the rate of water exchange leading to unrealistically slow exchange kinetics. In order to make progress in the development of improved models, we systematically optimize the Mg²⁺ parameters in combination with the TIP3P water model in a much larger parameter space than previously done. The results show that a long-ranged interaction potential and modified Lorentz-Berthelot combination rules allow us to accurately reproduce multiple experimental properties including the solvation free energy, the distances to the oxygens of the first hydration shell, the hydration number, the activity coefficient derivative in MgCl₂ solutions, the self-diffusion coefficient, and the binding affinity to the phosphate oxygen of RNA. Matching this broad range of thermodynamic properties, we present two sets of optimal parameters: MicroMg yields water exchange on the microsecond timescale in agreement with experiments. NanoMg yields water exchange on the nanosecond timescale facilitating the direct observation of ion-binding events. As shown for the example of the add A-riboswitch, the optimized parameters correctly reproduce the structure of specifically bound ions and permit the de novo prediction of Mg²⁺-binding sites in biomolecular simulations.

SUBMITTER: Grotz KK

PROVIDER: S-EPMC8047801 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:Alkali (Li(+), Na(+), K(+), Rb(+), and Cs(+)) and halide (F(-), Cl(-), Br(-), and I(-)) ions play an important role in many biological phenomena, roles that range from stabilization of biomolecular structure, to influence on biomolecular dynamics, to key physiological influence on homeostasis and signaling. To properly model ionic interaction and stability in atomistic simulations of biomolecular structure, dynamics, folding, catalysis, and function, an accurate model or representation of the monovalent ions is critically necessary. A good model needs to simultaneously reproduce many properties of ions, including their structure, dynamics, solvation, and moreover both the interactions of these ions with each other in the crystal and in solution and the interactions of ions with other molecules. At present, the best force fields for biomolecules employ a simple additive, nonpolarizable, and pairwise potential for atomic interaction. In this work, we describe our efforts to build better models of the monovalent ions within the pairwise Coulombic and 6-12 Lennard-Jones framework, where the models are tuned to balance crystal and solution properties in Ewald simulations with specific choices of well-known water models. Although it has been clearly demonstrated that truly accurate treatments of ions will require inclusion of nonadditivity and polarizability (particularly with the anions) and ultimately even a quantum mechanical treatment, our goal was to simply push the limits of the additive treatments to see if a balanced model could be created. The applied methodology is general and can be extended to other ions and to polarizable force-field models. Our starting point centered on observations from long simulations of biomolecules in salt solution with the AMBER force fields where salt crystals formed well below their solubility limit. The likely cause of the artifact in the AMBER parameters relates to the naive mixing of the Smith and Dang chloride parameters with AMBER-adapted Aqvist cation parameters. To provide a more appropriate balance, we reoptimized the parameters of the Lennard-Jones potential for the ions and specific choices of water models. To validate and optimize the parameters, we calculated hydration free energies of the solvated ions and also lattice energies (LE) and lattice constants (LC) of alkali halide salt crystals. This is the first effort that systematically scans across the Lennard-Jones space (well depth and radius) while balancing ion properties like LE and LC across all pair combinations of the alkali ions and halide ions. The optimization across the entire monovalent series avoids systematic deviations. The ion parameters developed, optimized, and characterized were targeted for use with some of the most commonly used rigid and nonpolarizable water models, specifically TIP3P, TIP4P EW, and SPC/E. In addition to well reproducing the solution and crystal properties, the new ion parameters well reproduce binding energies of the ions to water and the radii of the first hydration shells.

Dataset Information

Optimized Magnesium Force Field Parameters for Biomolecular Simulations with Accurate Solvation, Ion-Binding, and Water-Exchange Properties.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets