Parameters for Molecular Dynamics Simulations of Manganese-Containing Metalloproteins.

A set of geometrical parameters has been determined for single manganese metalloproteins for the AMBER force field, and ultimately to other force fields with a similar philosophy. Twelve (12) models from 9 different single-cluster manganese proteins were optimized and parametrized, using a bonded model approach. Mn-ligand bonds, Mn-ligand angles, and Restrained Electrostatic Potential charges for all the 74 residues in the first metal coordination sphere of each Mn metalloprotein were parametrized. The determined parameters were validated with molecular dynamics simulations and several statistics strategies were used to analyze the results. In addition, to validate the parametrized models, frequency and normal mode calculations were performed and comparisons were obtained for the overall structures both with quantum mechanics and molecular mechanics calculations. Linear and polynomial fittings were performed to estimate Mn-ligand bond force constants for generic manganese centers. Furthermore, averages are proposed for the main Mn-ligand angle interactions of typical manganese coordination centers: axial, square and triangular equatorial planes, and tetrahedral positions, for the different combinations of donor atoms from waters and hard ligands.

[1]  Patricia C Babbitt,et al.  Evolution of enzymatic activities in the enolase superfamily: stereochemically distinct mechanisms in two families of cis,cis-muconate lactonizing enzymes. , 2009, Biochemistry.

[2]  G J Kleywegt,et al.  A re-evaluation of the crystal structure of chloromuconate cycloisomerase. , 1996, Acta crystallographica. Section D, Biological crystallography.

[3]  Michael Dolg,et al.  Energy‐adjusted ab initio pseudopotentials for the first row transition elements , 1987 .

[4]  Glauciete S. Maciel,et al.  Study of charges transferability for use in force fields , 2006 .

[5]  V. Pecoraro,et al.  Recent advances in the understanding of the biological chemistry of manganese. , 1999, Current opinion in chemical biology.

[6]  M. Lawrence,et al.  The crystal structure of pneumococcal surface antigen PsaA reveals a metal-binding site and a novel structure for a putative ABC-type binding protein. , 1998, Structure.

[7]  Kenneth M Merz,et al.  Structural Survey of Zinc Containing Proteins and the Development of the Zinc AMBER Force Field (ZAFF). , 2010, Journal of chemical theory and computation.

[8]  K. Burke.  Perspective on density functional theory. , 2012, The Journal of chemical physics.

[9]  Maria João Ramos,et al.  The carboxylate shift in zinc enzymes: a computational study. , 2007, Journal of the American Chemical Society.

[10]  William L. Jorgensen,et al.  OPLS all‐atom force field for carbohydrates , 1997 .

[11]  H. Michel,et al.  Cytochrome c oxidase: structure and spectroscopy. , 1998, Annual review of biophysics and biomolecular structure.

[12]  Ruibo Wu,et al.  A Transferable Non-bonded Pairwise Force Field to Model Zinc Interactions in Metalloproteins. , 2011, Journal of chemical theory and computation.

[13]  Pedro Alexandrino Fernandes,et al.  Comparative analysis of the performance of commonly available density functionals in the determination of geometrical parameters for zinc complexes , 2009, J. Comput. Chem..

[14]  Ulf Ryde,et al.  Comparison of methods for deriving atomic charges from the electrostatic potential and moments , 1998 .

[15]  Pedro Alexandrino Fernandes,et al.  General performance of density functionals. , 2007, The journal of physical chemistry. A.

[16]  P. Siegbahn,et al.  Theoretical study of the catalytic reaction mechanism of MndD , 2006, JBIC Journal of Biological Inorganic Chemistry.

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

[18]  E. Nakamura,et al.  Mechanistic study of the manganese-catalyzed [2 + 2 + 2] annulation of 1,3-dicarbonyl compounds and terminal alkynes. , 2009, Journal of the American Chemical Society.

[19]  J. Cano,et al.  Two proton-one electron coupled transfer in iron and manganese superoxide dismutases: A density functional study , 2007 .

[20]  Michael Bühl,et al.  Geometries of Second-Row Transition-Metal Complexes from Density-Functional Theory. , 2007, Journal of chemical theory and computation.

[21]  M. Harding,et al.  The geometry of metal-ligand interactions relevant to proteins. , 1999, Acta crystallographica. Section D, Biological crystallography.

[22]  T. N. Bhat,et al.  The Protein Data Bank , 2000, Nucleic Acids Res..

[23]  P. P. Ewald.  Die Berechnung optischer und elektrostatischer Gitterpotentiale , 1921 .

[24]  Peter Comba,et al.  The Directionality of d‐Orbitals and Molecular‐Mechanics Calculations of Octahedral Transition‐Metal Compounds , 1995 .

[25]  William L. Jorgensen,et al.  OPLS potential functions for nucleotide bases. Relative association constants of hydrogen-bonded base pairs in chloroform , 1991 .

[26]  Alexander D. MacKerell.  Empirical force fields for biological macromolecules: Overview and issues , 2004, J. Comput. Chem..

[27]  W. Stallings,et al.  Structure-function in Escherichia coli iron superoxide dismutase: comparisons with the manganese enzyme from Thermus thermophilus. , 1995, Biochemistry.

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

[29]  P A Kollman,et al.  Calculation and prediction of binding free energies for the matrix metalloproteinases. , 2000, Journal of medicinal chemistry.

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

[31]  Rochus Schmid,et al.  Ab initio parametrized MM3 force field for the metal‐organic framework MOF‐5 , 2007, J. Comput. Chem..

[32]  T. Straatsma,et al.  Force-Field Development and Molecular Dynamics of [NiFe] Hydrogenase. , 2012, Journal of chemical theory and computation.

[33]  U. Ryde,et al.  Structure of reduced and oxidized manganese superoxide dismutase: a combined computational and experimental approach. , 2006, The journal of physical chemistry. B.

[34]  Kenneth M. Merz,et al.  Force Field Design for Metalloproteins , 1991 .

[35]  C. Lim,et al.  Empirical force fields for biologically active divalent metal cations in water. , 2006, The journal of physical chemistry. A.

[36]  L. Wackett,et al.  Crystallographic Comparison of Manganese- and Iron-Dependent Homoprotocatechuate 2,3-Dioxygenases , 2004, Journal of bacteriology.

[37]  Maria J Ramos,et al.  Chemically modified tetracyclines as inhibitors of MMP-2 matrix metalloproteinase: a molecular and structural study. , 2012, The journal of physical chemistry. B.

[38]  M. Klein,et al.  Modeling the charge distribution at metal sites in proteins for molecular dynamics simulations. , 2007, Journal of structural biology.

[39]  P. Kollman,et al.  Automatic atom type and bond type perception in molecular mechanical calculations. , 2006, Journal of molecular graphics & modelling.

[40]  Alexander D. MacKerell,et al.  An all-atom empirical energy function for the simulation of nucleic acids , 1995 .

[41]  V. DeRose,et al.  Multifrequency pulsed EPR studies of biologically relevant manganese(II) complexes , 2007, Applied magnetic resonance.

[42]  C. Gui,et al.  Differentiating a ligand's chemical requirements for allosteric interactions from those for protein binding. Phenylalanine inhibition of pyruvate kinase. , 2006, Biochemistry.

[43]  P. Siegbahn.  Quantum chemical studies of manganese centers in biology. , 2002, Current opinion in chemical biology.

[44]  M. Harding,et al.  The geometry of metal-ligand interactions relevant to proteins. II. Angles at the metal atom, additional weak metal-donor interactions. , 2000, Acta crystallographica. Section D, Biological crystallography.

[45]  Malcolm D. Walkinshaw,et al.  Metals in protein structures: a review of their principal features , 2010 .

[46]  P. Kollman,et al.  How well does a restrained electrostatic potential (RESP) model perform in calculating conformational energies of organic and biological molecules? , 2000 .

[47]  G. Borgstahl,et al.  The structure of human mitochondrial manganese superoxide dismutase reveals a novel tetrameric interface of two 4-helix bundles , 1992, Cell.

[48]  J. Koepke,et al.  High resolution crystal structure of Paracoccus denitrificans cytochrome c oxidase: new insights into the active site and the proton transfer pathways. , 2009, Biochimica et biophysica acta.

[49]  Daniel J. Sindhikara,et al.  Apo and nickel-bound forms of the Pyrococcus horikoshii species of the metalloregulatory protein: NikR characterized by molecular dynamics simulations. , 2009, Biochemistry.

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

[51]  Renxiao Wang,et al.  Systematic Derivation of AMBER Force Field Parameters Applicable to Zinc-Containing Systems. , 2010, Journal of chemical theory and computation.

[52]  Min Wei,et al.  Valence Force Field for Layered Double Hydroxide Materials Based on the Parameterization of Octahedrally Coordinated Metal Cations , 2012 .

[53]  David van der Spoel,et al.  The Origin of Layer Structure Artifacts in Simulations of Liquid Water. , 2006 .

[54]  A. Becke,et al.  Density-functional exchange-energy approximation with correct asymptotic behavior. , 1988, Physical review. A, General physics.

[55]  Sérgio F. Sousa,et al.  Effective tailor-made force field parameterization of the several Zn coordination environments in the puzzling FTase enzyme: opening the door to the full understanding of its elusive catalytic mechanism , 2006 .

[56]  M Karplus,et al.  Zinc binding in proteins and solution: A simple but accurate nonbonded representation , 1995, Proteins.

[57]  David van der Spoel,et al.  The Origin of Layer Structure Artifacts in Simulations of Liquid Water. , 2006, Journal of chemical theory and computation.

[58]  R. Friesner,et al.  Evaluation and Reparametrization of the OPLS-AA Force Field for Proteins via Comparison with Accurate Quantum Chemical Calculations on Peptides† , 2001 .

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

[60]  A Wlodawer,et al.  Structure of the catalytic domain of avian sarcoma virus integrase with a bound HIV-1 integrase-targeted inhibitor. , 1998, Proceedings of the National Academy of Sciences of the United States of America.

[61]  Yuzhen Ye,et al.  Solution structure of Mycobacterium tuberculosis NmtR in the apo state: insights into Ni(II)-mediated allostery. , 2012, Biochemistry.

[62]  Lihong Hu,et al.  Reorganization energy for internal electron transfer in multicopper oxidases. , 2011, The journal of physical chemistry. B.

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

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

[65]  M. Wagner,et al.  Copper complexes of mono- and ditopic [(methylthio)methyl]borates: missing links and linked systems en route to copper enzyme models. , 2008, Chemistry.

[66]  A. Becke.  Density-functional thermochemistry. , 1996 .

[67]  M. Ramos,et al.  Metals in proteins: cluster analysis studies , 2011, Journal of molecular modeling.

[68]  G. Petsko,et al.  Mechanism of the reaction catalyzed by mandelate racemase. 2. Crystal structure of mandelate racemase at 2.5-A resolution: identification of the active site and possible catalytic residues. , 1991, Biochemistry.

[69]  W. L. Jorgensen,et al.  Comparison of simple potential functions for simulating liquid water , 1983 .

[70]  Mikael P. Johansson,et al.  Charge parameterization of the metal centers in cytochrome c oxidase , 2008, J. Comput. Chem..

[71]  L. Banci.  Molecular dynamics simulations of metalloproteins. , 2003, Current opinion in chemical biology.

[72]  P Brandt,et al.  Deriving force field parameters for coordination complexes , 2001 .

[73]  G. Ciccotti,et al.  Numerical Integration of the Cartesian Equations of Motion of a System with Constraints: Molecular Dynamics of n-Alkanes , 1977 .

[74]  Stepan Sklenak,et al.  A molecular dynamics exploration of the catalytic mechanism of yeast cytosine deaminase. , 2005, The journal of physical chemistry. B.

[75]  U. Ryde,et al.  Reaction mechanism of manganese superoxide dismutase studied by combined quantum and molecular mechanical calculations and multiconfigurational methods. , 2009, The journal of physical chemistry. B.

[76]  Y. Pang,et al.  Successful molecular dynamics simulation of the zinc-bound farnesyltransferase using the cationic dummy atom approach. , 2000, Protein science : a publication of the Protein Society.

[77]  Thomas E. Cheatham,et al.  Quantum mechanically derived AMBER‐compatible heme parameters for various states of the cytochrome P450 catalytic cycle , 2012, J. Comput. Chem..

[78]  K. Merz,et al.  Insights into the structure and dynamics of the dinuclear zinc beta-lactamase site from Bacteroides fragilis. , 2002, Biochemistry.

[79]  E. W. Czerwinski,et al.  Structural and functional linkages between subunit interfaces in mammalian pyruvate kinase. , 2001, Journal of molecular biology.

[80]  Matthias Rarey,et al.  Modeling of metal interaction geometries for protein–ligand docking , 2007, Proteins.

[81]  Per E. M. Siegbahn,et al.  A Quantum Chemical Approach to the Study of Reaction Mechanisms of Redox-Active Metalloenzymes , 2001 .

[82]  Kenneth M. Merz,et al.  Binding of Bicarbonate to Human Carbonic Anhydrase II: A Continuum of Binding States , 1997 .

[83]  P. Fernandes,et al.  Molecular Dynamics Simulations: Difficulties, Solutions and Strategies for Treating Metalloenzymes , 2010 .

[84]  Ulf Ryde,et al.  Comparison of Methods to Obtain Force-Field Parameters for Metal Sites. , 2011, Journal of chemical theory and computation.