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.

ADVANCES IN NICA-DONNAN MODEL ON BINDING HEAVY METAL

Modeling and Mechanism of the Adsorption of Proton and Copper to Natural Bamboo Sawdust Using the NICA–Donnan Model

Metal ions binding onto a lignocellulosic substrate extracted from wheat bran: a NICA-Donnan approach

Kinetics of heavy metal adsorption and desorption in soil: Developing a unified model based on chemical speciation

A General Model For Kinetics Of Heavy Metal Adsorption And Desorption On Soils

Micro-interfacial Mechanism and Model of Metal Ions Adsorption on the Iron ( Hydr ) Oxides and Humic Substances : A Review

A Fixed-K Model for Metal-Humate Binding

Kinetics of Heavy Metal Dissociation from Natural Organic Matter: Roles of the Carboxylic and Phenolic Sites.

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

A General Model for Kinetics of Heavy Metal Adsorption and Desorption on Soils.

Modeling Kinetics of Heavy Metal Release from Field-Contaminated Soils: Roles of Soil Adsorbents and Binding Sites

Predicting Heavy Metal Partition Equilibrium in Soils: Roles of Soil Components and Binding Sites

Nanoscale Insights into the Interaction Mechanism Underlying the Adsorption and Retention of Heavy Metal Ions by Humic Acid

Multi-scale modeling of natural organic matter–heavy metal cations interactions: Aggregation and stabilization mechanisms

Interaction Structure and Affinity of Zwitterionic Amino Acids with Important Metal Cations (Cd2+, Cu2+, Fe3+, Hg2+, Mn2+, Ni2+ and Zn2+) in Aqueous Solution: A Theoretical Study

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

Molecular Dynamics Simulation of the Interaction between Common Metal Ions and Humic Acids

Quantifying Metal Ions Binding Onto Dissolved Organic Matter Using Log-Transformed Absorbance Spectra.

Nature differences of humic acids fractions induced by extracted sequence as explanatory factors for binding characteristics of heavy metals

Interaction Structure and Affinity of Zwitterionic Amino Acids with Important Metal Cations (Cd<SUP>2+</SUP>, Cu<SUP>2+</SUP>, Fe<SUP>3+</SUP>, Hg<SUP>2+</SUP>, Mn<SUP>2+</SUP>, Ni<SUP>2+</SUP> and Zn<SUP>2+</SUP>) in Aqueous Solution: A Theoretical Study

A Qicar Approach for Quantifying Binding Constants for Metal-Ligand Complexes