Parameters for molecular dynamics simulations of iron‐sulfur proteins

Iron‐sulfur proteins involved in electron transfer reactions have finely tuned redox potentials, which allow them to be highly efficient and specific. Factors such as metal center solvent exposure, interaction with charged residues, or hydrogen bonds between the ligand residues and amide backbone groups have all been pointed out to cause such specific redox potentials. Here, we derived parameters compatible with the AMBER force field for the metal centers of iron‐sulfur proteins and applied them in the molecular dynamics simulations of three iron‐sulfur proteins. We used density‐functional theory (DFT) calculations and Seminario's method for the parameterization. Parameter validation was obtained by matching structures and normal frequencies at the quantum mechanics and molecular mechanics levels of theory. Having guaranteed a correct representation of the protein coordination spheres, the amide H‐bonds and the water exposure to the ligands were analyzed. Our results for the pattern of interactions with the metal centers are consistent to those obtained by nuclear magnetic resonance spectroscopy (NMR) experiments and DFT calculations, allowing the application of molecular dynamics to the study of those proteins. © 2013 Wiley Periodicals, Inc.

[1]  W. M. Westler,et al.  NMR Investigations of Clostridium pasteurianum Rubredoxin. Origin of Hyperfine 1H, 2H, 13C, and 15N NMR Chemical Shifts in Iron−Sulfur Proteins As Determined by Comparison of Experimental Data with Hybrid Density Functional Calculations† , 1998 .

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

[3]  A. Wachters,et al.  Gaussian Basis Set for Molecular Wavefunctions Containing Third‐Row Atoms , 1970 .

[4]  J. Dewan,et al.  Models of amide-cysteine hydrogen bonding in rubredoxin: hydrogen bonding between amide and benzenethiolate in [(CH3)3NCH2CONH2]2[Co(SC6H5)4].cntdot.1/2CH3CN and [(CH3)3NCH2CONH2][SC6H5] , 1991 .

[5]  P. E. John,et al.  CALCULATING THE NUMBER OF SPANNING TREES IN A LABELED PLANAR MOLECULAR GRAPH WHOSE INNER DUAL IS A TREE , 1996 .

[6]  L. Sieker,et al.  Crystallographic refinement of rubredoxin at 1 x 2 A degrees resolution. , 1985, Journal of molecular biology.

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

[8]  W. M. Westler,et al.  Correlation between hydrogen bond lengths and reduction potentials in Clostridium pasteurianum rubredoxin. , 2003, Journal of the American Chemical Society.

[9]  Hongbin Li,et al.  Hydrogen bond strength modulates the mechanical strength of ferric-thiolate bonds in rubredoxin. , 2012, Journal of the American Chemical Society.

[10]  Z. Dauter,et al.  Zinc- and iron-rubredoxins from Clostridium pasteurianum at atomic resolution: a high-precision model of a ZnS4 coordination unit in a protein. , 1996, Proceedings of the National Academy of Sciences of the United States of America.

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

[12]  R. Huber,et al.  Crystal structure of desulforedoxin from Desulfovibrio gigas determined at 1.8 A resolution: a novel non-heme iron protein structure. , 1995, Journal of molecular biology.

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

[14]  M. Romão,et al.  Structural studies by X‐ray diffraction on metal substituted desulforedoxin, a rubredoxin‐type protein , 1999, Protein science : a publication of the Protein Society.

[15]  T. Ichiye,et al.  Molecular dynamics simulations of rubredoxin from Clostridium pasteurianum: Changes in structure and electrostatic potential duringredox reactions , 1995, Proteins.

[16]  G. Sheldrick,et al.  Crystal structure determination at 1.4 A resolution of ferredoxin from the green alga Chlorella fusca. , 1999, Structure.

[17]  V. Hornak,et al.  Comparison of multiple Amber force fields and development of improved protein backbone parameters , 2006, Proteins.

[18]  H. Beinert,et al.  Iron-sulfur clusters: nature's modular, multipurpose structures. , 1997, Science.

[19]  W. M. Westler,et al.  Paramagnetic NMR spectroscopy and density functional calculations in the analysis of the geometric and electronic structures of iron-sulfur proteins. , 2005, Inorganic chemistry.

[20]  Angelo Vedani,et al.  A new force field for modeling metalloproteins , 1990 .

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

[22]  W. M. Westler,et al.  Changes in hydrogen-bond strengths explain reduction potentials in 10 rubredoxin variants. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

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

[24]  Masao Kakudo,et al.  X-Ray Analysis of a [2Fe-2S] Ferredoxin from ‘Spirulina platensis. Main Chain Fold and Location of Side Chains at 2.5 Å Resolution , 1981 .

[25]  O. Livnah,et al.  Structural basis for the thermostability of ferredoxin from the cyanobacterium Mastigocladus laminosus. , 2005, Journal of molecular biology.

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

[27]  T. Darden,et al.  A smooth particle mesh Ewald method , 1995 .

[28]  David Baker,et al.  Prediction of structures of zinc‐binding proteins through explicit modeling of metal coordination geometry , 2010, Protein science : a publication of the Protein Society.

[29]  W. M. Westler,et al.  AMPLIFICATION OF ONE-BOND 1H/2H ISOTOPE EFFECTS ON 15N CHEMICAL SHIFTS IN CLOSTRIDIUM PASTEURIANUM RUBREDOXIN BY FERMI-CONTACT EFFECTS THROUGH HYDROGE N BONDS , 1998 .

[30]  H. Beinert.  Iron-sulfur proteins: ancient structures, still full of surprises , 2000, JBIC Journal of Biological Inorganic Chemistry.

[31]  L. Sieker,et al.  The structure of rubredoxin at 1.2 A resolution. , 1979, Journal of molecular biology.

[32]  Joost VandeVondele,et al.  Calculation of redox properties: understanding short- and long-range effects in rubredoxin. , 2007, The journal of physical chemistry. B.