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

A general molecular mechanics (MM) model for treating aqueous Cu2+ and Zn2+ 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 Cu2+ 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 Zn2+–water complexes in good agreement with QM calculations. Analyses of the MD trajectories revealed a six‐coordination first solvation shell for both Cu2+ and Zn2+ ions in aqueous solution, with ligation geometries falling in the range reported by previous studies. © 2012 Wiley Periodicals, Inc.

[1]  Robert J. Deeth,et al.  Factors influencing Jahn-Teller distortions in six-coordinate copper(II) and low-spin nickel(II) complexes , 1986 .

[2]  Clark R. Landis,et al.  Valence Bond Concepts Applied to the Molecular Mechanics Description of Molecular Shapes. 2. Applications to Hypervalent Molecules of the P-Block , 1993 .

[3]  Robert J Deeth.  Comprehensive molecular mechanics model for oxidized type I copper proteins: active site structures, strain energies, and entatic bulging. , 2007, Inorganic chemistry.

[4]  Daniel Spångberg,et al.  Distorted five-fold coordination of Cu2+(aq) from a Car-Parrinello molecular dynamics simulation. , 2005, Physical chemistry chemical physics : PCCP.

[5]  Peter Comba,et al.  Molecular modeling of inorganic compounds , 1995 .

[6]  Linus Pauling,et al.  A resonating-valence-bond theory of metals and intermetallic compounds , 1949, Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences.

[7]  T. Darden,et al.  Efficient particle-mesh Ewald based approach to fixed and induced dipolar interactions , 2000 .

[8]  Robert J. Deeth,et al.  Molecular modelling for transition metal complexes: Dealing with d-electron effects , 2009 .

[9]  Pengyu Y. Ren,et al.  Temperature and Pressure Dependence of the AMOEBA Water Model , 2004 .

[10]  J. Harvey.  On the accuracy of density functional theory in transition metal chemistry , 2006 .

[11]  Rudi van Eldik,et al.  Hydration and Water Exchange of Zinc(II) Ions. Application of Density Functional Theory , 1997 .

[12]  Pengyu Y. Ren,et al.  Towards accurate solvation dynamics of divalent cations in water using the polarizable amoeba force field: From energetics to structure. , 2006, The Journal of chemical physics.

[13]  Clark R. Landis,et al.  VALENCE BOND CONCEPTS APPLIED TO THE MOLECULAR MECHANICS DESCRIPTION OF MOLECULAR SHAPES. 3. APPLICATIONS TO TRANSITION METAL ALKYLS AND HYDRIDES , 1998 .

[14]  Christian Ochsenfeld,et al.  Linear and sublinear scaling formation of Hartree-Fock-type exchange matrices , 1998 .

[15]  Clark R. Landis,et al.  Valency and Bonding: Author index , 2005 .

[16]  Benjamin Williams-Hubbard,et al.  DommiMOE: An implementation of ligand field molecular mechanics in the molecular operating environment , 2005, J. Comput. Chem..

[17]  Gustavo E. Scuseria,et al.  IMPROVING SELF-CONSISTENT FIELD CONVERGENCE BY VARYING OCCUPATION NUMBERS , 1999 .

[18]  Magnus Sandström,et al.  Neutron diffraction study of aqueous transition metal salt solutions by isomorphic substitution , 1981 .

[19]  R. Car,et al.  First solvation shell of the Cu(II) aqua ion: evidence for fivefold coordination. , 2001, Science.

[20]  A. Becke.  Density-functional thermochemistry. III. The role of exact exchange , 1993 .

[21]  A. E. Carlsson.  ANGULAR FORCES AROUND TRANSITION METALS IN BIOMOLECULES , 1998 .

[22]  Robert J. Deeth,et al.  Analytical derivatives, ?? bonding and d???s mixing in the ligand field molecular mechanics methodElectronic supplementary information (ESI) available: theoretical background. See http://www.rsc.org/suppdata/cp/b2/b203815c/ , 2002 .

[23]  Ingmar Persson,et al.  A computational and experimental study on the Jahn-Teller effect in the hydrated copper (II) ion. Comparisons with hydrated nickel (II) ions in aqueous solution and solid Tutton's salts , 1989 .

[24]  Christian G Hartinger,et al.  Bioorganometallic chemistry--from teaching paradigms to medicinal applications. , 2009, Chemical Society reviews.

[25]  Shridhar R. Gadre,et al.  Ab initio quality one‐electron properties of large molecules: Development and testing of molecular tailoring approach , 2003, J. Comput. Chem..

[26]  Pengyu Y. Ren,et al.  Polarizable Atomic Multipole-based Molecular Mechanics for Organic Molecules. , 2011, Journal of chemical theory and computation.

[27]  Peter Comba,et al.  Modeling of molecular properties , 2011 .

[28]  Parr,et al.  Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density. , 1988, Physical review. B, Condensed matter.

[29]  Peter Comba,et al.  Molecular mechanics and the Jahn-Teller effect , 1994 .

[30]  Gérard Jaouen,et al.  Bioorganometallics: Future Trends in Drug Discovery, Analytical Chemistry, and Catalysis†,‡ , 2011 .

[31]  Pengyu Y. Ren,et al.  Consistent treatment of inter‐ and intramolecular polarization in molecular mechanics calculations , 2002, J. Comput. Chem..

[32]  R. Friesner,et al.  Ab initio quantum chemical and mixed quantum mechanics/molecular mechanics (QM/MM) methods for studying enzymatic catalysis. , 2005, Annual review of physical chemistry.

[33]  L. Helm,et al.  Water exchange on metal ions: experiments and simulations , 1999 .

[34]  Pengyu Y. Ren,et al.  Polarizable Atomic Multipole Water Model for Molecular Mechanics Simulation , 2003 .

[35]  Robert J Deeth,et al.  Molecular modelling of Jahn-Teller distortions in Cu(II)N6 complexes: elongations, compressions and the pathways in between. , 2006, Dalton transactions.

[36]  Thomas S. Hofer,et al.  Coordination and ligand exchange dynamics of solvated metal ions , 2005 .

[37]  Per Linse,et al.  Ewald summation and reaction field methods for potentials with atomic charges, dipoles, and polarizabilities , 2000 .

[38]  Thomas A. Halgren,et al.  The representation of van der Waals (vdW) interactions in molecular mechanics force fields: potential form, combination rules, and vdW parameters , 1992 .

[39]  Donald G Truhlar,et al.  Density functional theory for transition metals and transition metal chemistry. , 2009, Physical chemistry chemical physics : PCCP.

[40]  Richard M. Noyes,et al.  Mechanisms of Inorganic Reactions. A Study of Metal Complexes in Solution. , 1959 .

[41]  Kirk A Peterson,et al.  Systematically convergent basis sets for transition metals. I. All-electron correlation consistent basis sets for the 3d elements Sc-Zn. , 2005, The Journal of chemical physics.

[42]  Linus Pauling,et al.  The Nature of the Interatomic Forces in Metals , 1938 .

[43]  B. Thole.  Molecular polarizabilities calculated with a modified dipole interaction , 1981 .

[44]  L Pauling.  Angles between orthogonal spd bond orbitals with maximum strength. , 1976, Proceedings of the National Academy of Sciences of the United States of America.

[45]  Jean-Philip Piquemal,et al.  Polarizable molecular dynamics simulation of Zn(II) in water using the AMOEBA force field. , 2010, Journal of Chemical Theory and Computation.

[46]  Zilvinas Rinkevicius,et al.  Conformations, structural transitions and visible near-infrared absorption spectra of four-, five- and six-coordinated Cu(II) aqua complexes. , 2009, Physical chemistry chemical physics : PCCP.

[47]  J. Ponder,et al.  Force fields for protein simulations. , 2003, Advances in protein chemistry.

[48]  T. Bredow,et al.  Theory and range of modern semiempirical molecular orbital methods , 2005 .

[49]  H. A. Jahn,et al.  Stability of Polyatomic Molecules in Degenerate Electronic States. I. Orbital Degeneracy , 1937 .

[50]  A. V. van Duin,et al.  Development of the ReaxFF reactive force field for describing transition metal catalyzed reactions, with application to the initial stages of the catalytic formation of carbon nanotubes. , 2005, The journal of physical chemistry. A.

[51]  T. Keyes,et al.  Classical simulations with the POLIR potential describe the vibrational spectroscopy and energetics of hydration: divalent cations, from solvation to coordination complex. , 2011, Journal of the American Chemical Society.

[52]  William W. Porterfield,et al.  Inorganic chemistry : a unified approach , 1984 .

[53]  J. B. Hastings,et al.  Application of the EXAFS method to Jahn—Teller ions: static and dynamic behavior of Cu(H2O)62+ and Ci(H2O)2+6 in aqueous solution , 1981 .

[54]  A E Carlsson,et al.  The functional form of angular forces around transition metal ions in biomolecules. , 2001, Biophysical journal.

[55]  L. Curtiss,et al.  Compact contracted basis sets for third‐row atoms: Ga–Kr , 1990 .

[56]  Alan Grossfield,et al.  Simulation of Ca2+ and Mg2+ solvation using polarizable atomic multipole potential. , 2006, The journal of physical chemistry. B.

[57]  Henry Margenau,et al.  Theory of intermolecular forces , 1969 .

[58]  Kersti Hermansson,et al.  Development and validation of a ReaxFF reactive force field for Cu cation/water interactions and copper metal/metal oxide/metal hydroxide condensed phases. , 2010, The journal of physical chemistry. A.

[59]  Clark R. Landis,et al.  Valency and Bonding: Contents , 2005 .

[60]  T. H. Dunning.  Gaussian basis sets for use in correlated molecular calculations. I. The atoms boron through neon and hydrogen , 1989 .

[61]  Robert J Deeth,et al.  Molecular modelling for coordination compounds: Cu(II)-amine complexes. , 2005, Dalton transactions.

[62]  Pengyu Y. Ren,et al.  Ion solvation thermodynamics from simulation with a polarizable force field. , 2003, Journal of the American Chemical Society.

[63]  C. Jørgensen,et al.  The angular overlap model, an attempt to revive the ligand field approaches , 1965 .

[64]  I. Bersuker,et al.  Modern aspects of the Jahn-Teller effect theory and applications to molecular problems. , 2001, Chemical reviews.

[65]  C. Landis,et al.  Valence bond concepts applied to the molecular mechanics description of molecular shapes. 4. Transition metals with pi-bonds. , 2001, Journal of the American Chemical Society.

[66]  Robert J. Deeth,et al.  Structural and mechanistic insights into the oxy form of tyrosinase from molecular dynamics simulations , 2010, JBIC Journal of Biological Inorganic Chemistry.

[67]  Timothy Clark,et al.  Towards a ‘‘next generation’’ neglect of diatomic differential overlap based semiempirical molecular orbital technique , 2003 .