Molecular dynamics study of oxygenation reactions catalysed by the enzyme p-hydroxybenzoate hydroxylase

Combined quantum mechanical and molecular mechanical (QM/MM) methods have been used to model the rate-determining step of the title reaction for three different substrates, i.e., p-oxybenzoate dianion 1, p-hydroxybenzoate anion 2 and p-oxybenzoic acid anion 3. The system studied includes the enzyme, the cofactor, the substrate, and crystal water (ca. 7000 atoms total, ca. 100 QM atoms). Molecular dynamics simulations at the AM1/GROMOS level confirm that the dianion mechanism is favoured over both monoanion mechanisms, with a computed free energy barrier of 12 kcal mol−1. Geometry optimisations of model systems (QM) and of the complete enzyme (QM/MM) provide additional information on the intrinsic barriers and transition states.

[1]  M. Lah,et al.  The mobile flavin of 4-OH benzoate hydroxylase. , 1994, Science.

[2]  M. Stratford,et al.  Absorption spectra of radicals of substrates for p-hydroxybenzoate hydroxylase following electrophilic attack of the .OH radical in the 3 position. , 1987, Journal of Biological Chemistry.

[3]  G Vriend,et al.  Crystal structure of the p-hydroxybenzoate hydroxylase-substrate complex refined at 1.9 A resolution. Analysis of the enzyme-substrate and enzyme-product complexes. , 1989, Journal of molecular biology.

[4]  B. Entsch,et al.  Substrate recognition by "password" in p-hydroxybenzoate hydroxylase. , 1999, Biochemistry.

[5]  W. Thiel,et al.  Semiempirical treatment of electrostatic potentials and partial charges in combined quantum mechanical and molecular mechanical approaches , 1996, J. Comput. Chem..

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

[7]  B. Entsch,et al.  Catalytic function of tyrosine residues in para-hydroxybenzoate hydroxylase as determined by the study of site-directed mutants. , 1991, The Journal of biological chemistry.

[8]  P. Chapman,et al.  Fluoride elimination from substrates in hydroxylation reactions catalyzed by p-hydroxybenzoate hydroxylase. , 1980, The Journal of biological chemistry.

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

[10]  V Massey,et al.  Flavin-oxygen derivatives involved in hydroxylation by p-hydroxybenzoate hydroxylase. , 1976, The Journal of biological chemistry.

[11]  George P. Ford,et al.  New approach to the rapid semiempirical calculation of molecular electrostatic potentials based on the am1 wave function: Comparison with ab initio hf/6‐31g* results , 1993, J. Comput. Chem..

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

[13]  W. Goddard,et al.  Charge equilibration for molecular dynamics simulations , 1991 .

[14]  J. Vervoort,et al.  Frontier orbital study on the 4-hydroxybenzoate-3-hydroxylase-dependent activity with benzoate derivatives. , 1992, European journal of biochemistry.

[15]  W G Hol,et al.  Crystal structure of p-hydroxybenzoate hydroxylase complexed with its reaction product 3,4-dihydroxybenzoate. , 1988, Journal of Molecular Biology.

[16]  J. Vervoort,et al.  Modelling flavin and substrate substituent effects on the activation barrier and rate of oxygen transfer by p‐hydroxybenzoate hydroxylase , 2000, FEBS letters.

[17]  W. Hol,et al.  Crystal structures of wild-type p-hydroxybenzoate hydroxylase complexed with 4-aminobenzoate,2,4-dihydroxybenzoate, and 2-hydroxy-4-aminobenzoate and of the Tyr222Ala mutant complexed with 2-hydroxy-4-aminobenzoate. Evidence for a proton channel and a new binding mode of the flavin ring. , 1994, Biochemistry.

[18]  B. Entsch,et al.  Purification, properties, and oxygen reactivity of p-hydroxybenzoate hydroxylase from Pseudomonas aeruginosa. , 1989, Biochimica et biophysica acta.

[19]  M. Ortiz-Maldonado,et al.  Use of free energy relationships to probe the individual steps of hydroxylation of p-hydroxybenzoate hydroxylase: studies with a series of 8-substituted flavins. , 1999, Biochemistry.

[20]  H. Schreuder,et al.  Crystal structures of mutant Pseudomonas aeruginosa p-hydroxybenzoate hydroxylases: the Tyr201Phe, Tyr385Phe, and Asn300Asp variants. , 1995, Biochemistry.

[21]  M. Ludwig,et al.  pH-dependent structural changes in the active site of p-hydroxybenzoate hydroxylase point to the importance of proton and water movements during catalysis. , 1996, Biochemistry.

[22]  W. V. van Berkel,et al.  Role of Tyr201 and Tyr385 in substrate activation by p-hydroxybenzoate hydroxylase from Pseudomonas fluorescens. , 1993, European journal of biochemistry.

[23]  B. Vojnovic,et al.  Absorption spectra of radical forms of 2,4-dihydroxybenzoic acid, a substrate for p-hydroxybenzoate hydroxylase. , 1991, Journal of Biological Chemistry.

[24]  Eamonn F. Healy,et al.  Development and use of quantum mechanical molecular models. 76. AM1: a new general purpose quantum mechanical molecular model , 1985 .

[25]  A. Mulholland,et al.  CORRELATION OF CALCULATED ACTIVATION ENERGIES WITH EXPERIMENTAL RATE CONSTANTS FOR AN ENZYME CATALYZED AROMATIC HYDROXYLATION , 1998 .

[26]  W. F. V. Gunsteren,et al.  Moleküldynamik‐Computersimulationen; Methodik, Anwendungen und Perspektiven in der Chemie , 1990 .