A three-point method for evaluations of AMBER force field parameters: an application to copper-based artificial nucleases

We present the theoretical evaluation of new AMBER force field parameters for 12 copper-based nucleases with bis(2-pyridylmethyl) amine, 2,2′-dipyridylamine, imidazole, N,N-bis(2-benzimidazolylmethyl) amine and their derivative ligands based on first-principles electronic structure calculations at the B3LYP level of theory. A three-point approach was developed to accurately and efficiently evaluate the force field parameters for the copper-based nucleases with the ligands. The protocol of RESP atomic charges has been used to calculate the atomic charge distributions of the studied copper-based nucleases. The evaluated force field parameters and RESP atomic charges have been successfully applied in the testing molecular mechanics calculations and molecular dynamics simulations on the nucleases and the nuclease–DNA complexes, respectively. It has been demonstrated that the developed force field parameters and atomic charges can consistently reproduce molecular geometries and conformations in the available X-ray crystal structures and can reasonably predict the interaction properties of the nucleases with DNA. The developed force field parameters in this work provide an extension of the AMBER force field for its application to computational modeling and simulations of the copper-based artificial nucleases associated with DNA.

[1]  Robert James Deeth,et al.  A test of ligand field molecular mechanics as an efficient alternative to QM/MM for modelling metalloproteins: the structures of oxidised type I copper centres. , 2006, Chemical communications.

[2]  W. A. Ortiz,et al.  Phosphate diester hydrolysis and DNA damage promoted by new cis-aqua/hydroxy copper(II) complexes containing tridentate imidazole-rich ligands. , 2003, Inorganic chemistry.

[3]  A. Bond,et al.  Isomerism in copper(II) chloride complexes of bis(2-pyridylmethyl)amine and N-substituted derivatives: Synthesis and X-ray structural characterisation , 2007 .

[4]  K. Karlin,et al.  Reactivity study of a hydroperoxodicopper(II) complex: hydroxylation, dehydrogenation, and ligand cross-link reactions. , 2006, Inorganic chemistry.

[5]  P. Kollman,et al.  New-generation amber united-atom force field. , 2006, The journal of physical chemistry. B.

[6]  P. Fernandes,et al.  Molecular dynamics simulations of the enzyme Cu, Zn superoxide dismutase. , 2006, The journal of physical chemistry. B.

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

[8]  F. Hampel,et al.  The activation of tertiary carboxamides in metal complexes: an experimental and theoretical study on the methanolysis of acylated bispicolylamine copper(II) complexes. , 2004, Inorganic chemistry.

[9]  M. Greenberg,et al.  A minor groove binding copper-phenanthroline conjugate produces direct strand breaks via beta-elimination of 2-deoxyribonolactone. , 2002, Journal of the American Chemical Society.

[10]  P. Kollman,et al.  An approach to computing electrostatic charges for molecules , 1984 .

[11]  K. Karlin,et al.  Changing selectivity of DNA oxidation from deoxyribose to Guanine by ligand design and a new binuclear copper complex. , 2005, Journal of the American Chemical Society.

[12]  Miaoli Zhu,et al.  Crystal structure and nuclease activity of mono(1,10-phenanthroline) copper complex. , 2003, Journal of inorganic biochemistry.

[13]  R. S. Brown,et al.  Rapid three-step cleavage of RNA and DNA model systems promoted by a dinuclear Cu(II) complex in methanol. energetic origins of the catalytic efficacy. , 2007, Journal of the American Chemical Society.

[14]  Junmei Wang,et al.  Development and testing of a general amber force field , 2004, J. Comput. Chem..

[15]  T. D. Stack,et al.  Structure and spectroscopy of copper-dioxygen complexes. , 2004, Chemical reviews.

[16]  P. Kollman,et al.  A well-behaved electrostatic potential-based method using charge restraints for deriving atomic char , 1993 .

[17]  D. Truhlar,et al.  A new local density functional for main-group thermochemistry, transition metal bonding, thermochemical kinetics, and noncovalent interactions. , 2006, The Journal of chemical physics.

[18]  David S. Goodsell,et al.  Automated docking using a Lamarckian genetic algorithm and an empirical binding free energy function , 1998 .

[19]  F. V. Pamatong,et al.  Double-Strand Cleavage of DNA by a Monofunctional Transition Metal Cleavage Agent , 1996 .

[20]  Peter Comba,et al.  Inorganic and bioinorganic molecular mechanics modeling—the problem of the force field parameterization , 2003 .

[21]  D. Sigman.  Nuclease activity of 1,10-phenanthroline-copper ion , 1986 .

[22]  C. Breneman,et al.  Determining atom‐centered monopoles from molecular electrostatic potentials. The need for high sampling density in formamide conformational analysis , 1990 .

[23]  Tore Brinck,et al.  Force field parameterization of copper(I)-olefin systems from density functional calculations , 1997 .

[24]  P. Kollman,et al.  A Second Generation Force Field for the Simulation of Proteins, Nucleic Acids, and Organic Molecules , 1995 .

[25]  Robert J Deeth,et al.  On the performance of ligand field molecular mechanics for model complexes containing the peroxido-bridged [Cu2O2]2+ center. , 2008, Inorganic chemistry.

[26]  John SantaLucia,et al.  AMBER Force Field Parameters for the Naturally Occurring Modified Nucleosides in RNA. , 2007, Journal of chemical theory and computation.

[27]  Xichen Li,et al.  Evaluation of Amber force field parameters for copper(II) with pyridylmethyl-amine and benzimidazolylmethyl-amine ligands: A quantum chemical study , 2008 .

[28]  M. Halcrow,et al.  A new synthesis of bis(2-{pyrid-2-yl}ethyl)amine (LH) from bis(2-{pyrid-2-yl}ethyl)hydroxylamine (LOH), and the copper-dependent reduction of LOH to LHElectronic supplementary information (ESI) available: Full synthetic procedures and analytical data for the compounds in this study. See http://www.r , 2003 .

[29]  Ning Gan,et al.  Oxidative DNA strand scission induced by a trinuclear copper(II) complex. , 2004, Inorganic chemistry.

[31]  K. Karlin,et al.  Targeted strand scission of DNA substrates by a tricopper(II) coordination complex. , 2002, Journal of the American Chemical Society.

[32]  Donald G Truhlar,et al.  Density functionals for inorganometallic and organometallic chemistry. , 2005, The journal of physical chemistry. A.

[33]  H. Sugiyama,et al.  Mechanism of DNA Strand Scission Induced by (1,10-Phenanthroline)copper Complex: Major Direct DNA Cleavage Is Not through 1‘,2‘-Dehydronucleotide Intermediate nor β-Elimination of Forming Ribonolactone , 2000 .

[34]  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.

[35]  U. Singh,et al.  A NEW FORCE FIELD FOR MOLECULAR MECHANICAL SIMULATION OF NUCLEIC ACIDS AND PROTEINS , 1984 .

[36]  David E Reichert,et al.  Molecular modeling of bifunctional chelate peptide conjugates. 1. Copper and indium parameters for the AMBER force field. , 2001, Inorganic chemistry.

[37]  M. Karplus,et al.  CHARMM: A program for macromolecular energy, minimization, and dynamics calculations , 1983 .

[38]  Hui-Bie F. Xu,et al.  DNA damage by copper(II) complexes: coordination-structural dependence of reactivities. , 1999, Journal of inorganic biochemistry.

[39]  D. Truhlar,et al.  How well can new-generation density functionals describe the energetics of bond-dissociation reactions producing radicals? , 2008, The journal of physical chemistry. A.

[40]  J. Pople,et al.  Self‐Consistent Molecular‐Orbital Methods. IX. An Extended Gaussian‐Type Basis for Molecular‐Orbital Studies of Organic Molecules , 1971 .

[41]  F. V. Pamatong,et al.  Nonrandom Double Strand Cleavage of DNA by a Monofunctional Metal Complex: Mechanistic Studies , 1996 .

[42]  D. Wahnon,et al.  Reactivity of Copper(II) Hydroxides and Copper(II) Alkoxides for Cleaving an Activated Phosphate Diester , 1995 .

[43]  W. Goddard,et al.  Chromophore-in-Protein Modeling of the Structures and Resonance Raman Spectra for Type 1 Copper Proteins , 1998 .

[44]  K. Karlin,et al.  Unsymmetrical dicopper complexes. Direct observation of reversible O2 binding in a copper monooxygenase model system , 1991 .

[45]  Zaida Luthey-Schulten,et al.  Classical force field parameters for the heme prosthetic group of cytochrome c , 2004, J. Comput. Chem..

[46]  NATHAN E. SCHULTZ,et al.  Benchmarking approximate density functional theory for s/d excitation energies in 3d transition metal cations , 2008, J. Comput. Chem..

[47]  K. Karlin,et al.  Efficient and specific strand scission of DNA by a dinuclear copper complex: comparative reactivity of complexes with linked tris(2-pyridylmethyl)amine moieties. , 2002, Journal of the American Chemical Society.

[48]  Kenneth D. Karlin,et al.  Facile amide hydrolysis effected by dinuclear copper complexes , 1993 .

[49]  B. Nair,et al.  Oxidative cleavage of DNA by tridentate copper (II) complex. , 2003, Journal of inorganic biochemistry.

[50]  Peter A. Kollman,et al.  Application of the multimolecule and multiconformational RESP methodology to biopolymers: Charge derivation for DNA, RNA, and proteins , 1995, J. Comput. Chem..

[51]  Alexander D. MacKerell,et al.  All-atom empirical potential for molecular modeling and dynamics studies of proteins. , 1998, The journal of physical chemistry. B.

[52]  Peter Comba,et al.  A new molecular mechanics force field for the oxidized form of blue copper proteins , 2002, J. Comput. Chem..

[53]  J. Bashkin.  Introduction to RNA/DNA Cleavage. , 1998, Chemical reviews.

[54]  W. Goddard,et al.  UFF, a full periodic table force field for molecular mechanics and molecular dynamics simulations , 1992 .

[55]  J. Cowan,et al.  Efficient Inorganic Deoxyribonucleases. Greater than 50-Million-Fold Rate Enhancement in Enzyme-Like DNA Cleavage , 2000 .

[56]  Robert J. Deeth,et al.  The ligand field molecular mechanics model and the stereoelectronic effects of d and s electrons , 2001 .

[57]  K. Karlin,et al.  Selective DNA strand scission with binuclear copper complexes: implications for an active Cu2-O2 species. , 2006, Journal of the American Chemical Society.

[58]  T. Tullius,et al.  Oxidative Strand Scission of Nucleic Acids: Routes Initiated by Hydrogen Abstraction from the Sugar Moiety. , 1998, Chemical reviews.

[59]  Donald G Truhlar,et al.  Databases for transition element bonding: metal-metal bond energies and bond lengths and their use to test hybrid, hybrid meta, and meta density functionals and generalized gradient approximations. , 2005, The journal of physical chemistry. A.

[60]  Yan Zhao,et al.  Exchange-correlation functional with broad accuracy for metallic and nonmetallic compounds, kinetics, and noncovalent interactions. , 2005, The Journal of chemical physics.

[61]  Hans W. Horn,et al.  Fully optimized contracted Gaussian basis sets for atoms Li to Kr , 1992 .

[62]  Christian Buning,et al.  Molecular mechanics modeling of organic backbone of metal‐free and coordinated ligands , 1998 .

[63]  Jun Yu Li,et al.  Electronic effect of different positions of the -NO2 group on the DNA-intercalator of chiral complexes [Ru(bpy)2L]2+ (L =o-npip, m-npip and p-npip). , 2005, Dalton transactions.

[64]  Yan Jin,et al.  DNA cleavage by copper-ATCUN complexes. Factors influencing cleavage mechanism and linearization of dsDNA. , 2005, Journal of the American Chemical Society.

[65]  A. Neves,et al.  Structural, electrochemical and spectroscopic characterization of a new [Cu2L2(μ-Cl)2]2[CuCl4]2 dimer complex , 2004 .

[66]  Liang-Nian Ji,et al.  Cleavage of double-strand DNA by linear and triangular trinuclear copper complexes. , 2006, Journal of inorganic biochemistry.

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

[68]  Zijian Guo,et al.  Design of artificial metallonucleases with oxidative mechanism , 2007 .

[69]  Donald G Truhlar,et al.  PdnCO (n = 1,2): accurate Ab initio bond energies, geometries, and dipole moments and the applicability of density functional theory for fuel cell modeling. , 2006, The journal of physical chemistry. B.

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

[71]  U. Ryde,et al.  Geometry, reduction potential, and reorganization energy of the binuclear Cu(A) site, studied by density functional theory. , 2001, Journal of the American Chemical Society.

[72]  Huai Sun,et al.  Computer simulations of poly(ethylene oxide): force field, pvt diagram and cyclization behaviour , 1997 .

[73]  AKIFUMI ODA,et al.  New AMBER force field parameters of heme iron for cytochrome P450s determined by quantum chemical calculations of simplified models , 2005, J. Comput. Chem..

[74]  C. Burrows,et al.  Mechanisms of DNA cleavage by copper complexes of 3-clip-phen and of its conjugate with a distamycin analogue. , 2000, Nucleic acids research.

[75]  S. Oikawa,et al.  Site-specific DNA damage induced by NADH in the presence of copper(II): role of active oxygen species. , 1996, Biochemistry.

[76]  David M. Golden,et al.  Hydrocarbon Bond Dissociation Energies , 1982 .

[77]  Nicolai Lehnert,et al.  Direct hydrogen-atom abstraction by activated bleomycin: an experimental and computational study. , 2006, Journal of the American Chemical Society.

[78]  D. Truhlar,et al.  Improved description of nuclear magnetic resonance chemical shielding constants using the M06-L meta-generalized-gradient-approximation density functional. , 2008, The journal of physical chemistry. A.

[79]  Krishnan Raghavachari,et al.  Highly correlated systems. Ionization energies of first row transition metals Sc–Zn , 1989 .

[80]  Judith N. Burstyn,et al.  Toward the development of metal-based synthetic nucleases and peptidases: a rationale and progress report in applying the principles of coordination chemistry , 1998 .

[81]  Andrew E. Torda,et al.  The GROMOS biomolecular simulation program package , 1999 .

[82]  Yi-zhi Li,et al.  A trinuclear copper(II) complex of 2,4,6-tris(di-2-pyridylamine)-1,3,5-triazine shows prominent DNA cleavage activity. , 2007, Inorganic chemistry.

[83]  K Schulten,et al.  VMD: visual molecular dynamics. , 1996, Journal of molecular graphics.

[84]  D. Odom,et al.  Recognition and reaction of metallointercalators with DNA. , 1999, Chemical reviews.

[85]  Peter Comba,et al.  A new molecular mechanics force field for the design of oxotechnetium(V) and oxorhenium(V) radiopharmaceuticals , 2006 .

[86]  Donald G Truhlar,et al.  Design of Density Functionals by Combining the Method of Constraint Satisfaction with Parametrization for Thermochemistry, Thermochemical Kinetics, and Noncovalent Interactions. , 2006, Journal of chemical theory and computation.

[87]  T. Sakurai,et al.  Molecular orbital analysis of active site of oxidized azurin: Dependency of electronic properties on molecular structure , 2005 .

[88]  Robert J. Deeth,et al.  Molecular Mechanics for Coordination Complexes: The Impact of Adding d-Electron Stabilization Energies , 1995 .

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

[90]  Yi-zhi Li,et al.  Oxidative DNA cleavage promoted by multinuclear copper complexes: activity dependence on the complex structure. , 2006, Chemistry.