Dataset Information

A valence bond model for aqueous Cu(II) and Zn(II) ions in the AMOEBA polarizable force field.

ABSTRACT: A general molecular mechanics (MM) model for treating aqueous Cu(2+) and Zn(2+) ions was developed based on valence bond (VB) theory and incorporated into the atomic multipole optimized energetics for biomolecular applications (AMOEBA) polarizable force field. Parameters were obtained by fitting MM energies to that computed by ab initio methods for gas-phase tetra- and hexa-aqua metal complexes. Molecular dynamics (MD) simulations using the proposed AMOEBA-VB model were performed for each transition metal ion in aqueous solution, and solvent coordination was evaluated. Results show that the AMOEBA-VB model generates the correct square-planar geometry for gas-phase tetra-aqua Cu(2+) complex and improves the accuracy of MM model energetics for a number of ligation geometries when compared to quantum mechanical (QM) computations. On the other hand, both AMOEBA and AMOEBA-VB generate results for Zn(2+)-water complexes in good agreement with QM calculations. Analyses of the MD trajectories revealed a six-coordination first solvation shell for both Cu(2+) and Zn(2+) ions in aqueous solution, with ligation geometries falling in the range reported by previous studies.

SUBMITTER: Xiang JY

PROVIDER: S-EPMC3584243 | biostudies-literature | 2013 Apr

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

A valence bond model for aqueous Cu(II) and Zn(II) ions in the AMOEBA polarizable force field.

Xiang Jin Yu JY Ponder Jay W JW

Journal of computational chemistry 20121205 9

A general molecular mechanics (MM) model for treating aqueous Cu(2+) and Zn(2+) ions was developed based on valence bond (VB) theory and incorporated into the atomic multipole optimized energetics for biomolecular applications (AMOEBA) polarizable force field. Parameters were obtained by fitting MM energies to that computed by ab initio methods for gas-phase tetra- and hexa-aqua metal complexes. Molecular dynamics (MD) simulations using the proposed AMOEBA-VB model were performed for each transi ...[more]

PMID: 23212979

Similar Datasets

Project description:The hydration free energy, structure, and dynamics of the zinc divalent cation are studied using a polarizable force field in molecular dynamics simulations. Parameters for the Zn(2+) are derived from gas-phase ab initio calculation of Zn(2+)-water dimer. The Thole-based dipole polarization is adjusted based on the Constrained Space Orbital Variations (CSOV) calculation while the Symmetry Adapted Perturbation Theory (SAPT) approach is also discussed. The vdW parameters of Zn(2+) have been obtained by comparing the AMOEBA Zn(2+)-water dimerization energy with results from several theory levels and basis sets over a range of distances. Molecular dynamics simulations of Zn(2+) solvation in bulk water are subsequently performed with the polarizable force field. The calculated first-shell water coordination number, water residence time and free energy of hydration are consistent with experimental and previous theoretical values. The study is supplemented with extensive Reduced Variational Space (RVS) and Electron Localization Function (ELF) computations in order to unravel the nature of the bonding in Zn(2+)(H(2)O)(n) (n=1,6) complexes and to analyze the charge transfer contribution to the complexes. Results show that the importance of charge transfer decreases as the size of Zn-water cluster grows due to anticooperativity and to changes in the nature of the metal-ligand bonds. Induction could be dominated by polarization when the system approaches condensed-phase and the covelant effects are eliminated from the Zn(II)-water interaction. To construct an "effective" classical polarizable potential for Zn(2+) in bulk water, one should therefore avoid over-fitting to the ab initio charge transfer energy of Zn(2+)-water dimer. Indeed, in order to avoid overestimation of condensed-phase many-body effects, which is crucial to the transferability of polarizable molecular dynamics, charge transfer should not be included within the classical polarization contribution and should preferably be either incorporated in to the pairwise van der Waals contribution or treated explicitly.

Dataset Information

A valence bond model for aqueous Cu(II) and Zn(II) ions in the AMOEBA polarizable force field.

Publications

A valence bond model for aqueous Cu(II) and Zn(II) ions in the AMOEBA polarizable force field.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets