A comparative study of various computational approaches in calculating the structure of pyridoxal 5′‐phosphate (PLP)‐dependent β‐lyase protein. The importance of protein environment

Various computational approaches, using molecular mechanics (Amber), semiempirical (AM1), density functional (B3LYP), and various ONIOM methods, have been comparatively investigated for the structure of Escherichia coli NifS CsdB protein. The structure of the entire monomer containing 407 amino acid residues and 579 surrounding water molecules has been optimized. The full geometry optimization in the “active site‐only” approach (including only active site atoms) has been found to give the largest root‐mean‐square (RMS) deviation from the X‐ray structure; a much better agreement has been achieved by keeping the atoms leading to the backbones of some amino acids frozen in their positions in the X‐ray structure. The best agreement has been attained by including the surrounding protein in the calculations using the two‐layer ONIOM (B3LYP:Amber) approach. The results presented in this study conclusively demonstrate the importance of the protein/active‐site interaction on the active‐site structure of the enzyme. The present theoretical study represents the largest system studied at the ONIOM level to date, containing 7992 atoms, including 84 atoms in the QM region and rest in the MM region. © 2005 Wiley Periodicals, Inc. J Comput Chem 26: 443–446, 2005

[1]  Martin Karplus,et al.  Calculation of ground and excited state potential surfaces of conjugated molecules. I. Formulation and parametrization , 1972 .

[2]  Y. Hata,et al.  Structure of external aldimine of Escherichia coli CsdB, an IscS/NifS homolog: implications for its specificity toward selenocysteine. , 2002, Journal of biochemistry.

[3]  Benny G. Johnson,et al.  Linear scaling density functional calculations via the continuous fast multipole method , 1996 .

[4]  K. Morokuma,et al.  Investigation of the S0S1 excitation in bacteriorhodopsin with the ONIOM(MO:MM) hybrid method , 2003 .

[5]  M. Blomberg,et al.  A theoretical study of the mechanism for peptide hydrolysis by thermolysin , 2002, JBIC Journal of Biological Inorganic Chemistry.

[6]  Per E. M. Siegbahn,et al.  A Theoretical Study of the Mechanism for the Reductive Half-Reaction of Pea Seedling Amine Oxidase (PSAO) , 2001 .

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

[8]  Sason Shaik,et al.  Quantum mechanical/molecular mechanical investigation of the mechanism of C-H hydroxylation of camphor by cytochrome P450cam: theory supports a two-state rebound mechanism. , 2004, Journal of the American Chemical Society.

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

[10]  Robert H. White,et al.  Cysteine desulfurase activity indicates a role for NIFS in metallocluster biosynthesis. , 1993, Proceedings of the National Academy of Sciences of the United States of America.

[11]  Michael J. Frisch,et al.  Achieving linear scaling in exchange-correlation density functional quadratures , 1996 .

[12]  R. White,et al.  Mechanism for the desulfurization of L-cysteine catalyzed by the nifS gene product. , 1994, Biochemistry.

[13]  Jian Li,et al.  Quantum chemical studies of intermediates and reaction pathways in selected enzymes and catalytic synthetic systems. , 2004, Chemical reviews.

[14]  K. Morokuma,et al.  Effects of the protein environment on the structure and energetics of active sites of metalloenzymes. ONIOM study of methane monooxygenase and ribonucleotide reductase. , 2002, Journal of the American Chemical Society.

[15]  Victor Guallar,et al.  Cytochrome P450CAM enzymatic catalysis cycle: a quantum mechanics/molecular mechanics study. , 2004, Journal of the American Chemical Society.

[16]  Feliu Maseras,et al.  IMOMM: A new integrated ab initio + molecular mechanics geometry optimization scheme of equilibrium structures and transition states , 1995, J. Comput. Chem..

[17]  Y. Hata,et al.  Structure of a NifS homologue: X-ray structure analysis of CsdB, an Escherichia coli counterpart of mammalian selenocysteine lyase , 2000 .

[18]  Michael J. Frisch,et al.  Achieving Linear Scaling for the Electronic Quantum Coulomb Problem , 1996, Science.

[19]  J. Sudmeier,et al.  Cadmium-113 nuclear magnetic resonance studies of 113Cd(ii)-substituted human carbonic anhydrase B. , 1977, Journal of the American Chemical Society.

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

[21]  K. Morokuma,et al.  A NEW ONIOM IMPLEMENTATION IN GAUSSIAN98. PART I. THE CALCULATION OF ENERGIES, GRADIENTS, VIBRATIONAL FREQUENCIES AND ELECTRIC FIELD DERIVATIVES , 1999 .

[22]  K. Morokuma,et al.  On the application of the IMOMO (integrated molecular orbital + molecular orbital) method , 2000 .

[23]  Sason Shaik,et al.  A predictive pattern of computed barriers for C-h hydroxylation by compound I of cytochrome p450. , 2004, Journal of the American Chemical Society.

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

[25]  C. Lima.  Analysis of the E. coli NifS CsdB protein at 2.0 A reveals the structural basis for perselenide and persulfide intermediate formation. , 2002, Journal of molecular biology.

[26]  T. Koga.  HYPERBOLIC COSINE FUNCTIONS APPLIED TO ATOMIC ROOTHAAN-HARTREE-FOCK WAVE FUNCTIONS , 1998 .

[27]  K. Morokuma,et al.  Does the Active Site of Mammalian Glutathione Peroxidase (GPx) Contain Water Molecules? An ONIOM Study , 2004 .

[28]  A. Warshel,et al.  Structure/function correlations of proteins using MM, QM/MM, and related approaches: methods, concepts, pitfalls, and current progress. , 2003, Advances in protein chemistry.

[29]  R. Huber,et al.  Crystal structure of the cystine C-S lyase from Synechocystis: stabilization of cysteine persulfide for FeS cluster biosynthesis. , 2000, Proceedings of the National Academy of Sciences of the United States of America.

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

[31]  W. Thiel,et al.  Hybrid Models for Combined Quantum Mechanical and Molecular Mechanical Approaches , 1996 .

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