Charge parameterization of the metal centers in cytochrome c oxidase

Reliable atomic point charges are of key importance for a correct description of the electrostatic interactions when performing classical, force field based simulations. Here, we present a systematic procedure for point charge derivation, based on quantum mechanical methodology suited for the systems at hand. A notable difference to previous procedures is to include an outer region around the actual system of interest. At the cost of increasing the system sizes, here up to 265 atoms, including the surroundings achieves near‐neutrality for the systems as well as structural stability, important factors for reliable charge distributions. In addition, the common problem of converting between CH bonds and CC bonds at the border vanishes. We apply the procedure to the four redox‐active metal centers of cytochrome c oxidase: CuA, haem a, haem a3, and CuB. Several relevant charge and ligand states are considered. Charges for two different force fields, CHARMM and AMBER, are presented. © 2007 Wiley Periodicals, Inc. J Comput Chem, 2008

[1]  H. Michel,et al.  Titration behavior of residues at the entrance of the D-pathway of cytochrome c oxidase from paracoccus denitrificans investigated by continuum electrostatic calculations. , 2005, Biophysical journal.

[2]  J. Fee,et al.  Electron paramagnetic resonance studies of the soluble CuA protein from the cytochrome ba3 of Thermus thermophilus. , 1996, Biophysical journal.

[3]  Ilme Schlichting,et al.  Structure and chemistry of cytochrome P450. , 2005, Chemical reviews.

[4]  C. Wstermeire.  Structure at 2.7 A resolution of the Paracococcus dendrificans two-subunit cytochrome c oxidase coupled with an antibody Fv fragment , 1997 .

[5]  B. Rabenstein,et al.  Calculation of protonation patterns in proteins with structural relaxation and molecular ensembles – application to the photosynthetic reaction center , 1998, European Biophysics Journal.

[6]  G. Hummer,et al.  Water-gated mechanism of proton translocation by cytochrome c oxidase. , 2003, Biochimica et biophysica acta.

[7]  Frank Neese,et al.  The Purple Mixed-Valence CuA Center in Nitrous-oxide Reductase: EPR of the Copper-63-, Copper-65-, and Both Copper-65- and [15N]Histidine-Enriched Enzyme and a Molecular Orbital Interpretation , 1996 .

[8]  A. Stuchebrukhov,et al.  Electrostatic study of the proton pumping mechanism in bovine heart cytochrome C oxidase. , 2004, Journal of the American Chemical Society.

[9]  M. Wikström.  Proton pump coupled to cytochrome c oxidase in mitochondria. , 1977 .

[10]  M. Blomberg,et al.  Methods and models for studying mechanisms of redox-active enzymes , 2005, Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[11]  Ilya Belevich,et al.  Proton-coupled electron transfer drives the proton pump of cytochrome c oxidase , 2006, Nature.

[12]  A. Stuchebrukhov,et al.  Combined DFT and electrostatics study of the proton pumping mechanism in cytochrome c oxidase. , 2006, Biochimica et biophysica acta.

[13]  Sason Shaik,et al.  Theoretical perspective on the structure and mechanism of cytochrome P450 enzymes. , 2005, Chemical reviews.

[14]  P. Siegbahn,et al.  The effect of electron correlation on the metal-ligand interaction in iron pentacarbonyl (Fe(CO)5) , 1985 .

[15]  Wei Zhang,et al.  A point‐charge force field for molecular mechanics simulations of proteins based on condensed‐phase quantum mechanical calculations , 2003, J. Comput. Chem..

[16]  K Schulten,et al.  Oxygen and proton pathways in cytochrome c oxidase , 1998, Proteins.

[17]  T. Straatsma,et al.  Multiplet splittings and other properties from density functional theory: an assessment in iron–porphyrin systems , 2005 .

[18]  Gerrit Groenhof,et al.  GROMACS: Fast, flexible, and free , 2005, J. Comput. Chem..

[19]  H. Michel,et al.  The coupling of electron transfer and proton translocation: electrostatic calculations on Paracoccus denitrificans cytochrome c oxidase. , 1998, Biophysical journal.

[20]  Mikael P. Johansson,et al.  Change in electron and spin density upon electron transfer to haem. , 2002, Biochimica et biophysica acta.

[21]  Jun Ge,et al.  Can the semiempirical PM3 scheme describe iron‐containing bioinorganic molecules? , 2006, J. Comput. Chem..

[22]  R. Gennis,et al.  Mass spectrometric determination of dioxygen bond splitting in the "peroxy" intermediate of cytochrome c oxidase. , 1999, Proceedings of the National Academy of Sciences of the United States of America.

[23]  M. Blomberg,et al.  Transition-metal systems in biochemistry studied by high-accuracy quantum chemical methods. , 2000, Chemical reviews.

[24]  D L Harris,et al.  High-valent intermediates of heme proteins and model compounds. , 2001, Current opinion in chemical biology.

[25]  R.B.M.R Souza,et al.  Cargas Atômicas em Moléculas , 1996 .

[26]  J. Almlöf,et al.  Density functional theoretical study of oxo(porphyrinato)iron(IV) complexes, models of peroxidase compounds I and II , 1994 .

[27]  G. Babcock,et al.  Coordination geometries and vibrational properties of cytochromes alpha and alpha 3 in cytochrome oxidase from Soret excitation Raman spectroscopy. , 1981, Biochemistry.

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

[29]  G H Loew,et al.  Role of the heme active site and protein environment in structure, spectra, and function of the cytochrome p450s. , 2000, Chemical reviews.

[30]  A. Schäfer,et al.  Fully optimized contracted Gaussian basis sets of triple zeta valence quality for atoms Li to Kr , 1994 .

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

[32]  H. Michel,et al.  Dynamic water networks in cytochrome C oxidase from Paracoccus denitrificans investigated by molecular dynamics simulations. , 2004, Biophysical journal.

[33]  S. Ferguson-Miller,et al.  Heme/Copper Terminal Oxidases. , 1996, Chemical reviews.

[34]  S. Yoshikawa,et al.  The low-spin heme of cytochrome c oxidase as the driving element of the proton-pumping process , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[35]  P. C. Hariharan,et al.  The influence of polarization functions on molecular orbital hydrogenation energies , 1973 .

[36]  A. Klamt,et al.  COSMO : a new approach to dielectric screening in solvents with explicit expressions for the screening energy and its gradient , 1993 .

[37]  A. Puustinen,et al.  The calcium binding site in cytochrome aa3 from Paracoccus denitrificans. , 1999, Biochemistry.

[38]  G. Babcock,et al.  Oxygen activation and the conservation of energy in cell respiration , 1992, Nature.

[39]  A. Puustinen,et al.  Gating of proton and water transfer in the respiratory enzyme cytochrome c oxidase. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

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

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

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

[43]  Per E. M. Siegbahn,et al.  The performance of hybrid DFT for mechanisms involving transition metal complexes in enzymes , 2006, JBIC Journal of Biological Inorganic Chemistry.

[44]  M. Wikström,et al.  Coordination of CuB in reduced and CO-liganded states of cytochrome bo3 from Escherichia coli. Is chloride ion a cofactor? , 1999, Biochemistry.

[45]  Michele Parrinello,et al.  Equilibrium Geometries and Electronic Structure of Iron−Porphyrin Complexes: A Density Functional Study , 1997 .

[46]  M. Saraste,et al.  Soluble CuA-binding domain from the Paracoccus cytochrome c oxidase. , 1993, The Journal of biological chemistry.

[47]  D. Hartree.  The Wave Mechanics of an Atom with a Non-Coulomb Central Field. Part I. Theory and Methods , 1928, Mathematical Proceedings of the Cambridge Philosophical Society.

[48]  M. Wikström.  Cytochrome c oxidase: 25 years of the elusive proton pump. , 2004, Biochimica et biophysica acta.

[49]  T. Vangberg,et al.  Molecular structures and electron distributions of higher-valent iron and manganese porphyrins: Density functional theory calculations and some preliminary open-shell coupled-cluster results , 2001 .

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

[51]  M. Wikström,et al.  The K-pathway revisited: a computational study on cytochrome c oxidase. , 2006, Biochimica et biophysica acta.

[52]  S. Tsuzuki,et al.  Comparison of atomic charge distributions obtained from different procedures: basis set and electron correlation effects , 1996 .

[53]  G. Hummer,et al.  Structure and dynamics of a proton shuttle in cytochrome c oxidase , 1998 .

[54]  Holger Gohlke,et al.  The Amber biomolecular simulation programs , 2005, J. Comput. Chem..

[55]  F. Matthias Bickelhaupt,et al.  Voronoi deformation density (VDD) charges: Assessment of the Mulliken, Bader, Hirshfeld, Weinhold, and VDD methods for charge analysis , 2004, J. Comput. Chem..

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

[57]  Mikael P. Johansson,et al.  Spin and charge distribution in iron porphyrin models: a coupled cluster and density-functional study. , 2004, The Journal of chemical physics.

[58]  Mikael P. Johansson,et al.  The spin distribution in low-spin iron porphyrins. , 2002, Journal of the American Chemical Society.

[59]  H. Michel,et al.  Direct evidence for a tyrosine radical in the reaction of cytochrome c oxidase with hydrogen peroxide. , 1999, Biochemistry.

[60]  P. Hohenberg,et al.  Inhomogeneous Electron Gas , 1964 .

[61]  P. G. Gassman,et al.  Electron correlation in tetrapyrroles. Ab initio calculations on porphyrin and the tautomers of chlorin , 1993 .

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

[63]  K. Schulten,et al.  Imaging alpha-hemolysin with molecular dynamics: ionic conductance, osmotic permeability, and the electrostatic potential map. , 2005, Biophysical journal.

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

[65]  F. Neese,et al.  Probing the ground state of the purple mixed valence CuA center in nitrous oxide reductase: a CW ENDOR (X-band) study of the 65Cu, 15N-histidine labeled enzyme and interpretation of hyperfine couplings by molecular orbital calculations , 1998, JBIC Journal of Biological Inorganic Chemistry.

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

[67]  S. Ferguson-Miller,et al.  Energy transduction: proton transfer through the respiratory complexes. , 2006, Annual review of biochemistry.

[68]  H. Gray,et al.  The Cu_A Center of a Soluble Domain from Thermus Cytochrome ba_3. An NMR Investigation of the Paramagnetic Protein , 1996 .

[69]  Hartmut Michel,et al.  Structure at 2.8 Å resolution of cytochrome c oxidase from Paracoccus denitrificans , 1995, Nature.

[70]  C. DeMaso,et al.  Oxygen activation and reduction in respiration: involvement of redox-active tyrosine 244. , 2000, Science.

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

[72]  J. Perdew,et al.  Density-functional approximation for the correlation energy of the inhomogeneous electron gas. , 1986, Physical review. B, Condensed matter.

[73]  Hans W. Horn,et al.  ELECTRONIC STRUCTURE CALCULATIONS ON WORKSTATION COMPUTERS: THE PROGRAM SYSTEM TURBOMOLE , 1989 .

[74]  R. S. Mulliken.  Electronic Population Analysis on LCAO–MO Molecular Wave Functions. I , 1955 .

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

[76]  W. Kohn,et al.  Self-Consistent Equations Including Exchange and Correlation Effects , 1965 .

[77]  H. Michel,et al.  Structure at 2.7 A resolution of the Paracoccus denitrificans two-subunit cytochrome c oxidase complexed with an antibody FV fragment. , 1997, Proceedings of the National Academy of Sciences of the United States of America.

[78]  Jon Baker,et al.  Classical chemical concepts from ab initio SCF calculations , 1985 .

[79]  T. Tomizaki,et al.  Structures of metal sites of oxidized bovine heart cytochrome c oxidase at 2.8 A , 1995, Science.

[80]  Yifan Song,et al.  Calculated proton uptake on anaerobic reduction of cytochrome C oxidase: is the reaction electroneutral? , 2006, Biochemistry.

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

[82]  D. Pérahia,et al.  Internal and interfacial dielectric properties of cytochrome c from molecular dynamics in aqueous solution. , 1995, Proceedings of the National Academy of Sciences of the United States of America.

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

[84]  W. Woodruff,et al.  Time-resolved Raman detection of .nu.(Fe-O) in an early intermediate in the reduction of oxygen by cytochrome oxidase , 1989 .

[85]  V. Fock,et al.  Näherungsmethode zur Lösung des quantenmechanischen Mehrkörperproblems , 1930 .

[86]  David E. Bernholdt,et al.  High performance computational chemistry: An overview of NWChem a distributed parallel application , 2000 .

[87]  Gilda H. Loew,et al.  Structure and Spectra of Ferrous Dioxygen and Reduced Ferrous Dioxygen Model Cytochrome P450 , 1998 .

[88]  Angelo Carotti,et al.  Three‐dimensional model of the human aromatase enzyme and density functional parameterization of the iron‐containing protoporphyrin IX for a molecular dynamics study of heme‐cysteinato cytochromes , 2006, Proteins.