Abstract:Metal ions exist in almost half of the proteins in the protein databank and they serve as structural, electron-transfer and catalytic elements in the metabolic processes of organisms. Molecular Dynamics (MD) simulation is a powerful tool that provides information about biomolecular systems at the atomic level. Coupled with the growth in computing power, algorithms like the Particle Mesh Ewald (PME) method have become the accepted standard when dealing with long-range interactions in MD simulations. The nonbonded model of metal ions consists of an electrostatic plus 12-6 Lennard Jones (LJ) potential and is used largely because of its speed relative to more accurate models. In previous work we found that ideal parameters do not exist that reproduce several experimental properties for M(II) ions simultaneously using the nonbonded model coupled with the PME method due to the underestimation of metal ion-ligand interactions. Via a consideration of the nature of the nonbonded model, we proposed that the observed error largely arises from overlooking charge-induced dipole interactions. The electrostatic plus 12-6 LJ potential model works reasonably well for neutral systems but does struggle with more highly charged systems. In the present work we designed and parameterized a new nonbonded model for metal ions by adding a 1/r4 term to the 12-6 model. We call it the 12-6-4 LJ-type nonbonded model due to its mathematical construction. Parameters were determined for 16 +2 metal ions for the TIP3P, SPC/E and TIP4PEW water models. The final parameters reproduce the experimental hydration free energies (HFE), ion-oxygen distances (IOD) in the first solvation shell and coordination numbers (CN) accurately for the metal ions investigated. Preliminary tests on MgCl2 at different concentrations in aqueous solution and Mg2+--nucleic acid systems show reasonable results suggesting that the present parameters can work in mixed systems. The 12-6-4 LJ-type nonbonded model is readily adopted into standard force fields like AMBER, CHARMM and OPLS-AA with only a modest computational overhead. The new nonbonded model doesn't consider charge-transfer effects explicitly and, hence, may not suitable for the simulation of systems where charge-transfer effects play a decisive role.

Parameterization and validation of solvation corrected atomic radii.

Variable Atomic Radii For Continuum-Solvent Electrostatics Calculation

Accuracy of Discrete-Continuum Solvation Model for Cations: A Benchmark Study

Development and test of highly accurate endpoint free energy methods. 1: Evaluation of ABCG2 charge model on solvation free energy prediction and optimization of atom radii suitable for more accurate solvation free energy prediction by the PBSA method

Systematic Parameterization of Monovalent Ions Employing the Nonbonded Model

Computer Simulations of Aqua Metal Ions for Accurate Reproduction of Hydration Free Energies and Structures

Parameterization of Highly Charged Metal Ions Using the 12-6-4 LJ-Type Nonbonded Model in Explicit Water

Parameterization of Monovalent Ions for the OPC3, OPC, TIP3P-FB, and TIP4P-FB Water Models

Balancing the Atomic Waistline: Isodensity‐based Scrf Radii for Main‐group Elements and Transition Metals

Experimental Compilation and Computation of Hydration Free Energies for Ionic Solutes

Variable Van Der Waals Radii Derived from a Hybrid Gaussian Charge Distribution Model for Continuum-Solvent Electrostatic Calculations

Model for Aqueous Solvation Based on Class IV Atomic Charges and First Solvation Shell Effects

Crystal lattice properties fully determine short-range interaction parameters for alkali and halide ions

Bridging atomistic simulations and thermodynamic hydration models of aqueous electrolyte solutions

Optimal Dielectric Boundary for Binding Free Energy Estimates in the Implicit Solvent

Taking into Account the Ion-Induced Dipole Interaction in the Nonbonded Model of Ions

Development of Nonbonded Models for Metal Cations Using the Electronic Continuum Correction

Force fields for monovalent and divalent metal cations in TIP3P water based on thermodynamic and kinetic properties

Improved Lennard-Jones Parameters for Accurate Molecular Dynamics Simulations

Rational Design of Nonbonded Point Charge Models for Highly Charged Metal Cations with Lennard-Jones 12-6 Potential.