Polarizable molecular dynamics simulation of Zn(II) in water using the AMOEBA force field

Johnny C. Wu, Jean Philip Piquemal, Robin Chaudret, Peter Reinhardt, Pengyu Ren

Research output: Contribution to journalArticle

107 Citations (Scopus)

Abstract

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 Zn2+are derived from gas-phase ab initio calculation of the Zn2+-water dimer. The Thole-based dipole polarization is adjusted on the basis of the constrained space orbital variations (CSOV) calculation, while the symmetry adapted perturbation theory (SAPT) approach is also discussed. The vdW parameters of Zn2+ have been obtained by comparing the AMOEBA Zn2+-water dimerization energy with results from several theory levels and basis sets over a range of distances. Molecular dynamics simulations of Zn2+ 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 Zn2+ (H2O)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 the 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 the condensed phase and the covalent effects are eliminated from the Zn(II)-water interaction. To construct an "effective" classical polarizable potential for Zn2+ in bulk water, one should therefore avoid overfitting to the ab initio charge transfer energy of the Zn 2+-water dimer. Indeed, in order to avoid overestimation of the 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 into the pairwise van der Waals contribution or treated explicitly.

Original languageEnglish
Pages (from-to)2059-2070
Number of pages12
JournalJournal of Chemical Theory and Computation
Volume6
Issue number7
DOIs
Publication statusPublished - Jul 13 2010
Externally publishedYes

Fingerprint

field theory (physics)
Molecular dynamics
molecular dynamics
Water
Computer simulation
water
simulation
Charge transfer
charge transfer
Polarization
Hydration
Dimers
Free energy
hydration
polarization
free energy
dimers
Orbital calculations
Dimerization
Divalent Cations

ASJC Scopus subject areas

  • Computer Science Applications
  • Physical and Theoretical Chemistry

Cite this

Polarizable molecular dynamics simulation of Zn(II) in water using the AMOEBA force field. / Wu, Johnny C.; Piquemal, Jean Philip; Chaudret, Robin; Reinhardt, Peter; Ren, Pengyu.

In: Journal of Chemical Theory and Computation, Vol. 6, No. 7, 13.07.2010, p. 2059-2070.

Research output: Contribution to journalArticle

Wu, Johnny C. ; Piquemal, Jean Philip ; Chaudret, Robin ; Reinhardt, Peter ; Ren, Pengyu. / Polarizable molecular dynamics simulation of Zn(II) in water using the AMOEBA force field. In: Journal of Chemical Theory and Computation. 2010 ; Vol. 6, No. 7. pp. 2059-2070.
@article{069150346bc446d5b3e2042c952239d6,
title = "Polarizable molecular dynamics simulation of Zn(II) in water using the AMOEBA force field",
abstract = "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 Zn2+are derived from gas-phase ab initio calculation of the Zn2+-water dimer. The Thole-based dipole polarization is adjusted on the basis of the constrained space orbital variations (CSOV) calculation, while the symmetry adapted perturbation theory (SAPT) approach is also discussed. The vdW parameters of Zn2+ have been obtained by comparing the AMOEBA Zn2+-water dimerization energy with results from several theory levels and basis sets over a range of distances. Molecular dynamics simulations of Zn2+ 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 Zn2+ (H2O)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 the 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 the condensed phase and the covalent effects are eliminated from the Zn(II)-water interaction. To construct an {"}effective{"} classical polarizable potential for Zn2+ in bulk water, one should therefore avoid overfitting to the ab initio charge transfer energy of the Zn 2+-water dimer. Indeed, in order to avoid overestimation of the 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 into the pairwise van der Waals contribution or treated explicitly.",
author = "Wu, {Johnny C.} and Piquemal, {Jean Philip} and Robin Chaudret and Peter Reinhardt and Pengyu Ren",
year = "2010",
month = "7",
day = "13",
doi = "10.1021/ct100091j",
language = "English",
volume = "6",
pages = "2059--2070",
journal = "Journal of Chemical Theory and Computation",
issn = "1549-9618",
publisher = "American Chemical Society",
number = "7",

}

TY - JOUR

T1 - Polarizable molecular dynamics simulation of Zn(II) in water using the AMOEBA force field

AU - Wu, Johnny C.

AU - Piquemal, Jean Philip

AU - Chaudret, Robin

AU - Reinhardt, Peter

AU - Ren, Pengyu

PY - 2010/7/13

Y1 - 2010/7/13

N2 - 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 Zn2+are derived from gas-phase ab initio calculation of the Zn2+-water dimer. The Thole-based dipole polarization is adjusted on the basis of the constrained space orbital variations (CSOV) calculation, while the symmetry adapted perturbation theory (SAPT) approach is also discussed. The vdW parameters of Zn2+ have been obtained by comparing the AMOEBA Zn2+-water dimerization energy with results from several theory levels and basis sets over a range of distances. Molecular dynamics simulations of Zn2+ 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 Zn2+ (H2O)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 the 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 the condensed phase and the covalent effects are eliminated from the Zn(II)-water interaction. To construct an "effective" classical polarizable potential for Zn2+ in bulk water, one should therefore avoid overfitting to the ab initio charge transfer energy of the Zn 2+-water dimer. Indeed, in order to avoid overestimation of the 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 into the pairwise van der Waals contribution or treated explicitly.

AB - 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 Zn2+are derived from gas-phase ab initio calculation of the Zn2+-water dimer. The Thole-based dipole polarization is adjusted on the basis of the constrained space orbital variations (CSOV) calculation, while the symmetry adapted perturbation theory (SAPT) approach is also discussed. The vdW parameters of Zn2+ have been obtained by comparing the AMOEBA Zn2+-water dimerization energy with results from several theory levels and basis sets over a range of distances. Molecular dynamics simulations of Zn2+ 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 Zn2+ (H2O)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 the 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 the condensed phase and the covalent effects are eliminated from the Zn(II)-water interaction. To construct an "effective" classical polarizable potential for Zn2+ in bulk water, one should therefore avoid overfitting to the ab initio charge transfer energy of the Zn 2+-water dimer. Indeed, in order to avoid overestimation of the 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 into the pairwise van der Waals contribution or treated explicitly.

UR - http://www.scopus.com/inward/record.url?scp=78649507575&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=78649507575&partnerID=8YFLogxK

U2 - 10.1021/ct100091j

DO - 10.1021/ct100091j

M3 - Article

AN - SCOPUS:78649507575

VL - 6

SP - 2059

EP - 2070

JO - Journal of Chemical Theory and Computation

JF - Journal of Chemical Theory and Computation

SN - 1549-9618

IS - 7

ER -