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.

Effect of the intrinsic Josephson coupling on the pancake lattice in layered superconductors

Diatomic molecules of alkali-metal and alkaline-earth-metal atoms: interaction potentials, dipole moments, and polarizabilities

Theoretical Investigation on excited dipole bound states of alkali-containing diatomic anions

Bond-counting potentials -- A classical many-body model of covalent bonding with exact solutions in one dimension

Rational Design of Nonbonded Point Charge Models for Monovalent Ions with Lennard-Jones 12–6 Potential

Spin structure of diatomic van der Waal molecules of alkali atoms

Systematic Parameterization of Monovalent Ions Employing the Nonbonded Model

A new parametrizable model of molecular electronic structure

The diabatic valence-hole concept

Helium cluster ions: coherent charge sharing and the general trimerization trend

Interaction between LiH molecule and Li atom from state-of-the-art electronic structure calculations

Long range intermolecular forces in triatomic systems: connecting the atom-diatom and atom-atom-atom representations

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

Van der Waals molecules consisting of a zinc or cadmium atom interacting with an alkali-metal or alkaline-earth-metal atom

Near-threshold bound states of the dipole-dipole interaction

A simple model for the Pauli Repulsion with possible utility in QM, MM and Chemical Education

Accurate Predictions of Electron Binding Energies of Dipole-Bound Anions via Quantum Monte Carlo Methods

Universal Pairwise Interatomic van der Waals Potentials Based on Quantum Drude Oscillators

Near-threshold structures in inner-shell photoabsorption processes of CH_{4}, NH_{3}, H_{2}O, and HF

Diatomic molecules embedded with Morse potential and under the effect of external fields

A set of molecular models based on quantum mechanical ab initio calculations and thermodynamic data