Aspects of the mechanism of catalysis of glucose oxidase: A docking, molecular mechanics and quantum chemical study

The complex structure of glucose oxidase (GOX) with the substrate glucose was determined using a docking algorithm and subsequent molecular dynamics simulations. Semiempirical quantum chemical calculations were used to investigate the role of the enzyme and FAD co-enzyme in the catalytic oxidation of glucose. On the basis of a small active site model, substrate binding residues were determined and heats of formation were computed for the enzyme substrate complex and different potential products of the reductive half reaction. The influence of the protein environment on the active site model was estimated with a point charge model using a mixed QM/MM method. Solvent effects were estimated with a continuum model. Possible modes of action are presented in relation to experimental data and discussed with respect to related enzymes. The calculations indicate that the redox reaction of GOX differs from the corresponding reaction of free flavins as a consequence of the protein environment. One of the active site histidines is involved in substrate binding and stabilization of potential intermediates, whereas the second histidine is a proton acceptor. The former one, being conserved in a series of oxidoreductases, is also involved in the stabilization of a C4a-hydroperoxy dihydroflavin in the course of the oxidative half reaction.

[1]  R. Ornstein,et al.  Ab initio quantum mechanics analysis of imidazole C-H...O water hydrogen bonding and a molecular mechanics forcefield correction. , 1997, Journal of biomolecular structure & dynamics.

[2]  P. Kraulis.  A program to produce both detailed and schematic plots of protein structures , 1991 .

[3]  Z. Derewenda,et al.  (His)Cε-H···O=C< Hydrogen Bond in the Active Sites of Serine Hydrolases , 1994 .

[4]  R. Ornstein,et al.  A Theoretical Study of the Structures of Flavin in Different Oxidation and Protonation States , 1996 .

[5]  Donald G. Truhlar,et al.  Improved methods for semiempirical solvation models , 1995, J. Comput. Chem..

[6]  Timothy Clark,et al.  Molecular orbital studies of enzyme mechanisms. II. Catalytic oxidation of alcohols by liver alcohol dehydrogenase , 1993, J. Comput. Chem..

[7]  M. Meyer.  Computational study of imidazole and methyl imidazoles , 1994 .

[8]  Mark S. Gordon,et al.  General atomic and molecular electronic structure system , 1993, J. Comput. Chem..

[9]  George C. Shields,et al.  Ability of the PM3 quantum‐mechanical method to model intermolecular hydrogen bonding between neutral molecules , 1993, J. Comput. Chem..

[10]  Haiyan Liu,et al.  The reaction pathway of the isomerization of d‐xylose catalyzed by the enzyme d‐xylose isomerase: A theoretical study , 1997, Proteins.

[11]  V. Massey.  Activation of molecular oxygen by flavins and flavoproteins. , 1994, The Journal of biological chemistry.

[12]  Tapani A. Pakkanen,et al.  Ab initio models for receptor-ligand interactions in proteins. 3. Model assembly study of the proton transfer in the hydroxylation step of the catalytic mechanism of p-hydroxybenzoate hydroxylase , 1993 .

[13]  Robert J. Woods,et al.  Molecular Mechanical and Molecular Dynamic Simulations of Glycoproteins and Oligosaccharides. 1. GLYCAM_93 Parameter Development , 1995 .

[14]  J. Wouters,et al.  Electronic properties of flavins: Implications on the reactivity and absorption properties of flavoproteins , 1997 .

[15]  D. Nguyen,et al.  Simulation of the enzyme reaction mechanism of malate dehydrogenase. , 1997, Biochemistry.

[16]  P. Kollman,et al.  Chemical reaction mechanisms in vacuo, in solution and in enzyme fields: isomerization catalyzed by triose phosphate isomerase (TIM) , 1996 .

[17]  Peter A. Kollman,et al.  A comparison of the AM1 and PM3 semiempirical models for evaluating model compounds relevant to catalysis by serine proteases , 1991 .

[18]  J. Pople,et al.  Self—Consistent Molecular Orbital Methods. XII. Further Extensions of Gaussian—Type Basis Sets for Use in Molecular Orbital Studies of Organic Molecules , 1972 .

[19]  D Schomburg,et al.  Crystal structure of glucose oxidase from Aspergillus niger refined at 2.3 A resolution. , 1993, Journal of molecular biology.

[20]  G. Shields,et al.  Hydrogen bonding of nucleotide base pairs: Application of the PM3 method , 1994 .

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

[22]  E. Pai,et al.  Absolute stereochemistry of flavins in enzyme-catalyzed reactions. , 1986, Biochemistry.

[23]  M. Meyer.  Density functional study of isoalloxazine and C4a-hydroperoxidihydroisoalloxazine , 1997 .

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

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

[26]  J. Stewart.  Optimization of parameters for semiempirical methods I. Method , 1989 .

[27]  M. Weibel,et al.  The glucose oxidase mechanism. Interpretation of the pH dependence. , 1971, The Journal of biological chemistry.

[28]  A. Bacher,et al.  15N- and 13C-NMR investigations of glucose oxidase from Aspergillus niger. , 1991, European journal of biochemistry.

[29]  F. Mȕller.  Chemistry and Biochemistry of Flavoenzymes: Volume I , 1991 .

[30]  M. Sundaralingam,et al.  C-H...O hydrogen bonding in biology. , 1997, Trends in biochemical sciences.

[31]  V. Moliner,et al.  AM1 AND PM3 TRANSITION STRUCTURE FOR THE HYDRIDE TRANSFER. A MODEL OF REACTION CATALYZED BY DIHYDROFOLATE REDUCTASE , 1995 .

[32]  D. Cavener,et al.  GMC oxidoreductases. A newly defined family of homologous proteins with diverse catalytic activities. , 1992, Journal of molecular biology.

[33]  Bernd Beck,et al.  The natural atomic orbital point charge model for PM3: Multipole moments and molecular electrostatic potentials , 1994, J. Comput. Chem..

[34]  T. Darden,et al.  An ab Initio Quantum Mechanical Model for the Catalytic Mechanism of HIV-1 Protease , 1996 .

[35]  D. Schomburg,et al.  Hydrogen bonding and molecular surface shape complementarity as a basis for protein docking. , 1996, Journal of molecular biology.

[36]  R E Cachau,et al.  Inhibition and catalytic mechanism of HIV-1 aspartic protease. , 1996, Journal of molecular biology.

[37]  J. Voet,et al.  Electrostatic control of enzyme reactions: effect of ionic strength on the pKa of an essential acidic group on glucose oxidase. , 1981, Biochemistry.

[38]  T. Pakkanen,et al.  Model assembly study of the ligand binding by p‐hydroxybenzoate hydroxylase: Correlation between the calculated binding energies and the experimental dissociation constants , 1995, Proteins.

[39]  Brian D. Wladkowski,et al.  Transphosphorylation catalyzed by ribonuclease A: computational study using ab initio effective fragment potentials , 1995 .

[40]  J. Klinman,et al.  Effects of protein glycosylation on catalysis: changes in hydrogen tunneling and enthalpy of activation in the glucose oxidase reaction. , 1997, Biochemistry.

[41]  D. Schomburg,et al.  Semiempirical and ab initio study of closed and open shell derivatives of 10-methylisoalloxazine: a model of flavin redox states , 1996 .

[42]  W. R. Carper,et al.  The hydrolysis of 1,5-gluconolactone: semi-empirical methods and 13C NMR confirmation , 1992 .

[43]  Carol A. Venanzi,et al.  Computational analysis of binding affinity and neural response at the l-alanine receptor , 1995, J. Comput. Aided Mol. Des..

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

[45]  P. Brick,et al.  Crystal structure of cholesterol oxidase complexed with a steroid substrate: implications for flavin adenine dinucleotide dependent alcohol oxidases. , 1993, Biochemistry.