Protocatechuate hydroxylase is a novel group A flavoprotein monooxygenase with a unique substrate recognition mechanism.

Para-hydroxybenzoate hydroxylase (PHBH) is a group A flavoprotein monooxygenase that hydroxylates p-hydroxybenzoate to protocatechuate (PCA). Despite intensive studies of Pseudomonas aeruginosa p-hydroxybenzoate hydroxylase (PaPobA), the catalytic reactions of extremely diverse putative PHBH isozymes remain unresolved. We analyzed the phylogenetic relationships of known and predicted PHBHs and identified eight divergent clades. Clade F contains a protein that lacks the critical amino acid residues required for PaPobA to generate PHBH activity. Among proteins in this clade, Xylophilus ampelinus PobA (XaPobA) preferred PCA as a substrate and is the first known natural PCA 5-hydroxylase (PCAH). Crystal structures and kinetic properties revealed similar mechanisms of substrate carboxy group recognition between XaPobA and PaPobA. The unique Ile75, Met72, Val199, Trp201, and Phe385 residues of XaPobA form the bottom of a hydrophobic cavity with a shape that complements the 3-and 4-hydroxy groups of PCA and its binding site configuration. An interaction between the δ-sulfur atom of Met210 and the aromatic ring of PCA is likely to stabilize XaPobA-PCA complexes. The 4-hydroxy group of PCA forms a hydrogen bond with the main chain carbonyl of Thr294. These modes of binding constitute a novel substrate recognition mechanism that PaPobA lacks. This mechanism characterizes XaPobA and sheds light on the diversity of catalytic mechanisms of PobA-type PHBHs and group A flavoprotein monooxygenases.

Genetic redundancy of 4-hydroxybenzoate 3-hydroxylase genes ensures the catabolic safety of Pigmentiphaga sp. H8 in 3-bromo-4-hydroxybenzoate-contaminated habitats.

Genetic redundancy is prevalent in organisms and plays important roles in the evolution of biodiversity and adaptation to environmental perturbation. However, selective advantages of genetic redundancy in overcoming metabolic disturbance due to structural analogues have received little attention. Here, functional divergence of the three 4-hydroxybenzoate 3-hydroxylase (PHBH) genes (phbh1~3) was found in Pigmentiphaga sp. strain H8. The genes phbh1/phbh2 were responsible for 3-bromo-4-hydroxybenzoate (3-Br-4-HB, an anthropogenic pollutant) catabolism, whereas phbh3 was primarily responsible for 4-hydroxybenzoate (4-HB, a natural intermediate of lignin) catabolism. 3-Br-4-HB inhibited 4-HB catabolism by competitively binding PHBH3 and was toxic to strain H8 cells especially at high concentrations. The existence of phbh1/phbh2 not only enabled strain H8 to utilize 3-Br-4-HB but also ensured the catabolic safety of 4-HB. Molecular docking and site-directed mutagenesis analyses revealed that Val199 and Phe384 of PHBH1/PHBH2 were required for the hydroxylation activity towards 3-Br-4-HB. Phylogenetic analysis indicated that phbh1 and phbh2 originated from a common ancestor and evolved specifically in strain H8 to adapt to 3-Br-4-HB-contaminated habitats, whereas phbh3 evolved independently. This study deepens our understanding of selective advantages of genetic redundancy in prokaryote's metabolic robustness and reveals the factors driving the divergent evolution of redundant genes in adaptation to environmental perturbation.

Tuning of pKa values activates substrates in flavin-dependent aromatic hydroxylases.

Hydroxylation of substituted phenols by flavin-dependent monooxygenases is the first step of their biotransformation in various microorganisms. The reaction is thought to proceed via electrophilic aromatic substitution, catalyzed by enzymatic deprotonation of substrate, in single-component hydroxylases that use flavin as a cofactor (group A). However, two-component hydroxylases (group D), which use reduced flavin as a co-substrate, are less amenable to spectroscopic investigation. Herein, we employed 19F NMR in conjunction with fluorinated substrate analogs to directly measure pKa values and to monitor protein events in hydroxylase active sites. We found that the single-component monooxygenase 3-hydroxybenzoate 6-hydroxylase (3HB6H) depresses the pKa of the bound substrate analog 4-fluoro-3-hydroxybenzoate (4F3HB) by 1.6 pH units, consistent with previously proposed mechanisms. 19F NMR was applied anaerobically to the two-component monooxygenase 4-hydroxyphenylacetate 3-hydroxylase (HPAH), revealing depression of the pKa of 3-fluoro-4-hydroxyphenylacetate by 2.5 pH units upon binding to the C2 component of HPAH. 19F NMR also revealed a pKa of 8.7 ± 0.05 that we attributed to an active-site residue involved in deprotonating bound substrate, and assigned to His-120 based on studies of protein variants. Thus, in both types of hydroxylases, we confirmed that binding favors the phenolate form of substrate. The 9 and 14 kJ/mol magnitudes of the effects for 3HB6H and HPAH-C2, respectively, are consistent with pKa tuning by one or more H-bonding interactions. Our implementation of 19F NMR in anaerobic samples is applicable to other two-component flavin-dependent hydroxylases and promises to expand our understanding of their catalytic mechanisms.

Understanding the Molecular Mechanism Underlying the High Catalytic Activity of p-Hydroxybenzoate Hydroxylase Mutants for Producing Gallic Acid.

p-Hydroxybenzoate hydroxylase (PHBH) is a flavoprotein monooxygenase that catalyzes the hydroxylation of p-hydroxybenzoate (p-OHB) to 3,4-dihydroxybenzoate (3,4-DOHB). PHBH can bind to other benzoate derivatives in addition to p-OHB; however, hydroxylation does not occur on 3,4-DOHB. Replacement of Tyr385 with Phe forms a mutant, which enables the production of 3,4,5-trihydroxybenzonate (gallic acid) from 3,4-DOHB, although the catalytic activity of the mutant is quite low. In this study, we report how the L199V/Y385F double mutant exhibits activity for producing gallic acid 4.3-fold higher than that of the Y385F single mutant. This improvement in catalytic activity is primarily due to the suppression of a shunt reaction that wastes reduced nicotinamide adenine dinucleotide phosphate by producing H2O2. To further elucidate the molecular mechanism underlying this higher catalytic activity, we performed molecular dynamics simulations and quantum mechanics/molecular mechanics calculations, in addition to determining the crystal structure of the Y385F·3,4-DOHB complex. The simulations showed that the Y385F mutation facilitates the deprotonation of the 4-hydroxy group of 3,4-DOHB, which is necessary for initiating hydroxylation. Moreover, the L199V mutation in addition to the Y385F mutation allows the OH moiety in the peroxide group of C-(4a)-flavin hydroperoxide to come into the proximity of the C5 atom of 3,4-DOHB. Overall, this study provides a consistent explanation for the change in the catalytic activity of PHBH caused by mutations, which will enable us to better design an enzyme with different activities.

Hydroxyl Radical-Coupled Electron-Transfer Mechanism of Flavin-Dependent Hydroxylases.

Class A flavin-dependent hydroxylases (FdHs) catalyze the hydroxylation of organic compounds in a site- and stereoselective manner. In stark contrast, conventional synthetic routes require environmentally hazardous reagents and give modest yields. Thus, understanding the detailed mechanism of this class of enzymes is essential to their rational manipulation for applications in green chemistry and pharmaceutical production. Both electrophilic substitution and radical intermediate mechanisms have been proposed as interpretations of FdH hydroxylation rates and optical spectra. While radical mechanistic steps are often difficult to examine directly, modern quantum chemistry calculations combined with statistical mechanical approaches can yield detailed mechanistic models providing insights that can be used to differentiate reaction pathways. In the current work, we report quantum mechanical/molecular mechanical (QM/MM) calculations on the fungal TropB enzyme that shows an alternative reaction pathway in which hydroxylation through a hydroxyl radical-coupled electron-transfer mechanism is significantly favored over electrophilic substitution. Furthermore, QM/MM calculations on several modified flavins provide a more consistent interpretation of the experimental trends in the reaction rates seen experimentally for a related enzyme, para-hydroxybenzoate hydroxylase. These calculations should guide future enzyme and substrate design strategies and broaden the scope of biological spin chemistry.

PqsL uses reduced flavin to produce 2-hydroxylaminobenzoylacetate, a preferred PqsBC substrate in alkyl quinolone biosynthesis in Pseudomonas aeruginosa.

Alkyl hydroxyquinoline N-oxides (AQNOs) are antibiotic compounds produced by the opportunistic bacterial pathogen Pseudomonas aeruginosa They are products of the alkyl quinolone (AQ) biosynthetic pathway, which also generates the quorum-sensing molecules 2-heptyl-4(1H)-quinolone (HHQ) and 2-heptyl-3-hydroxy-4(1H)-quinolone (PQS). Although the enzymatic synthesis of HHQ and PQS had been elucidated, the route by which AQNOs are synthesized remained elusive. Here, we report on PqsL, the key enzyme for AQNO production, which structurally resembles class A flavoprotein monooxygenases such as p-hydroxybenzoate 3-hydroxylase (pHBH) and 3-hydroxybenzoate 6-hydroxylase. However, we found that unlike related enzymes, PqsL hydroxylates a primary aromatic amine group, and it does not use NAD(P)H as cosubstrate, but unexpectedly required reduced flavin as electron donor. We also observed that PqsL is active toward 2-aminobenzoylacetate (2-ABA), the central intermediate of the AQ pathway, and forms the unstable compound 2-hydroxylaminobenzoylacetate, which was preferred over 2-ABA as substrate of the downstream enzyme PqsBC. In vitro reconstitution of the PqsL/PqsBC reaction was feasible by using the FAD reductase HpaC, and we noted that the AQ:AQNO ratio is increased in an hpaC-deletion mutant of P. aeruginosa PAO1 compared with the ratio in the WT strain. A structural comparison with pHBH, the model enzyme of class A flavoprotein monooxygenases, revealed that structural features associated with NAD(P)H binding are missing in PqsL. Our study completes the AQNO biosynthetic pathway in P. aeruginosa, indicating that PqsL produces the unstable product 2-hydroxylaminobenzoylacetate from 2-ABA and depends on free reduced flavin as electron donor instead of NAD(P)H.

A peculiar IclR family transcription factor regulates para-hydroxybenzoate catabolism in Streptomyces coelicolor.

In Streptomyces coelicolor, we identified a para-hydroxybenzoate (PHB) hydroxylase, encoded by gene pobA (SCO3084), which is responsible for conversion of PHB into PCA (protocatechuic acid), a substrate of the ß-ketoadipate pathway which yields intermediates of the Krebs cycle. We also found that the transcription of pobA is induced by PHB and is negatively regulated by the product of SCO3209, which we named PobR. The product of this gene is highly unusual in that it is the apparent fusion of two IclR family transcription factors. Bioinformatic analyses, in vivo transcriptional assays, electrophoretic mobility shift assays (EMSAs), DNase I footprinting, and isothermal calorimetry (ITC) were used to elucidate the regulatory mechanism of PobR. We found that PobR loses its high affinity for DNA (i.e., the pobA operator) in the presence of PHB, the inducer of pobA transcription. PHB binds to PobR with a KD of 5.8 µM. Size-exclusion chromatography revealed that PobR is a dimer in the absence of PHB and a monomer in the presence of PHB. The crystal structure of PobR in complex with PHB showed that only one of the two IclR ligand binding domains was occupied, and defined how the N-terminal ligand binding domain engages the effector ligand.

Microbial community structure analysis of a benzoate-degrading halophilic archaeal enrichment.

A benzoate-degrading archaeal enrichment was developed using sediment samples from Rozel Point at Great Salt Lake, UT. The enrichment degraded benzoate as the sole carbon source at salinity ranging from 2.0 to 5.0 M NaCl with highest rate of degradation observed at 4.0 M. The enrichment was also tested for its ability to grow on other aromatic compounds such as 4-hydroxybenzoic acid (4-HBA), gentisic acid, protocatechuic acid (PCA), catechol, benzene and toluene as the sole sources of carbon and energy. Of these, the culture only utilized 4-HBA as the carbon source. To determine the initial steps in benzoate degradation pathway, a survey of ring-oxidizing and ring-cleaving genes was performed using degenerate PCR primers. Results showed the presence of 4-hydroxybenzoate 3-monooxygenase (4-HBMO) and protocatechuate 3, 4-dioxygenase (3,4-PCA) genes suggesting that the archaeal enrichment might degrade benzoate to 4-HBA that is further converted to PCA by 4-HBMO and, thus, formed PCA would undergo ring-cleavage by 3,4-PCA to form intermediates that enter the Krebs cycle. Small subunit rRNA gene-based diversity survey revealed that the enrichment consisted entirely of class Halobacteria members belonging to the genera Halopenitus, Halosarcina, Natronomonas, Halosimplex, Halorubrum, Salinarchaeum and Haloterrigena. Of these, Halopenitus was the dominant group accounting for almost 91 % of the total sequences suggesting their potential role in degrading oxygenated aromatic compounds at extreme salinity.

Hybrid Quantum Mechanics/Molecular Mechanics/Coarse Grained Modeling: A Triple-Resolution Approach for Biomolecular Systems.

We present a hybrid quantum mechanics/molecular mechanics/coarse-grained (QM/MM/CG) multiresolution approach for solvated biomolecular systems. The chemically important active-site region is treated at the QM level. The biomolecular environment is described by an atomistic MM force field, and the solvent is modeled with the CG Martini force field using standard or polarizable (pol-CG) water. Interactions within the QM, MM, and CG regions, and between the QM and MM regions, are treated in the usual manner, whereas the CG-MM and CG-QM interactions are evaluated using the virtual sites approach. The accuracy and efficiency of our implementation is tested for two enzymes, chorismate mutase (CM) and p-hydroxybenzoate hydroxylase (PHBH). In CM, the QM/MM/CG potential energy scans along the reaction coordinate yield reaction energies that are too large, both for the standard and polarizable Martini CG water models, which can be attributed to adverse effects of using large CG water beads. The inclusion of an atomistic MM water layer (10 Å for uncharged CG water and 5 Å for polarizable CG water) around the QM region improves the energy profiles compared to the reference QM/MM calculations. In analogous QM/MM/CG calculations on PHBH, the use of the pol-CG description for the outer water does not affect the stabilization of the highly charged FADHOOH-pOHB transition state compared to the fully atomistic QM/MM calculations. Detailed performance analysis in a glycine-water model system indicates that computation times for QM energy and gradient evaluations at the density functional level are typically reduced by 40-70% for QM/MM/CG relative to fully atomistic QM/MM calculations.

Conserved and non-conserved residues and their role in the structure and function of p-hydroxybenzoate hydroxylase.

In order to elucidate the molecular mechanism of the catalytic reaction and enzyme conformation, we substituted 53 conserved residues identified by aligning 92 p-hydroxybenzoate hydroxylase sequences and 19 non-conserved residues selected from crystallographic studies of Pseudomonas fluorescens NBRC14160 p-hydroxybenzoate hydroxylase with 19 other naturally occurring amino acids, yielding a database of 619 active single mutants. The database contained 365 and 254 active single mutants for 44/53 conserved residues and 19 non-conserved residues, respectively; the data included main activity, sub-activity for NADPH and NADPH reaction specificity. Active mutations were not observed for the G14, Q102, G160, E198, R220, R246, N300, F342 and G387 conserved residues, and only one active mutant was obtained at the G9, G11, G187, D286, Y201, R214 and G295 conserved residues and the S13, E32 and R42 non-conserved residues. Only seven active mutants with higher activity than the wild-type enzyme were observed at conserved residues, and only two were observed at non-conserved residues. The 365 mutants at conserved residues included 64 active mutants with higher NADPH reaction specificity than the wild-type enzyme, and some Y181X single mutants exhibited considerable changes in NADPH reaction specificity. A Y181X/L268G double-mutant database was constructed to computationally analyze the effects of these substitutions on structural conformation and function. These results indicated that some conserved or non-conserved residues are important for structural stability or enzyme function.

Structure and ligand binding properties of the epoxidase component of styrene monooxygenase .

Styrene monooxygenase (SMO) is a two-component flavoprotein monooxygenase that transforms styrene to styrene oxide in the first step of the styrene catabolic and detoxification pathway of Pseudomonas putida S12. The crystal structure of the N-terminally histidine-tagged epoxidase component of this system, NSMOA, determined to 2.3 A resolution, indicates the enzyme exists as a homodimer in which each monomer forms two distinct domains. The overall architecture is most similar to that of p-hydroxybenzoate hydroxylase (PHBH), although there are some significant differences in secondary structure. Structural comparisons suggest that a large cavity open to the surface forms the FAD binding site. At the base of this pocket is another cavity that likely represents the styrene binding site. Flavin binding and redox equilibria are tightly coupled such that reduced FAD binds apo NSMOA approximately 8000 times more tightly than the oxidized coenzyme. Equilibrium fluorescence and isothermal titration calorimetry data using benzene as a substrate analogue indicate that the oxidized flavin and substrate analogue binding equilibria of NSMOA are linked such that the binding affinity of each is increased by 60-fold when the enzyme is saturated with the other. A much weaker approximately 2-fold positive cooperative interaction is observed for the linked binding equilibria of benzene and reduced FAD. The low affinity of the substrate analogue for the reduced FAD complex of NSMOA is consistent with a preferred reaction order in which flavin reduction and reaction with oxygen precede the binding of styrene, identifying the apoenzyme structure as the key catalytic resting state of NSMOA poised to bind reduced FAD and initiate the oxygen reaction.

Toward accurate barriers for enzymatic reactions: QM/MM case study on p-hydroxybenzoate hydroxylase.

The hydroxylation reaction catalyzed by p-hydroxybenzoate hydroxylase has been investigated by quantum mechanical/molecular mechanical (QM/MM) calculations at different levels of QM theory. The solvated enzyme was modeled (approximately 23,000 atoms in total, 49 QM atoms). The geometries of reactant and transition state were optimized for ten representative pathways using semiempirical (AM1) and density functional (B3LYP) methods as QM components. Single-point calculations at B3LYP/MM optimized geometries were performed with local correlation methods [LMP2, LCCSD(T0)] and augmented triple-zeta basis sets. A careful validation of the latter approach with regard to all computational parameters indicates convergence of the QM contribution to the computed barriers to within approximately 1 kcal mol(-1). Comparison with the available experimental data supports this assessment.

Characterization of MobR, the 3-hydroxybenzoate-responsive transcriptional regulator for the 3-hydroxybenzoate hydroxylase gene of Comamonas testosteroni KH122-3s.

Comamonas testosteroni KH122-3s is an aerobic soil bacterium that utilizes 3-hydroxybenzoate as a sole carbon and energy source. In this strain, 3-hydroxybenzoate hydroxylase (MobA) acts on the initial step of the degradation to produce 3,4-dihydroxybenzoate, which is subsequently subjected to the meta-cleavage pathway leading to tricarboxylic acid cycle intermediates. Gene walking analysis of the upstream region of mobA revealed an open reading frame (mobR) that encodes a transcriptional regulator of the MarR family. Here, we report that MobR negatively regulates the expression of mobA, and that the repression is relieved by binding of 3-hydroxybenzoate, the substrate for MobA. A primer extension experiment was performed to determine the transcription start site for mobA and identified it at 83 bp upstream of the mobA start codon, accompanied by a typical sigma70-type promoter. The mobR gene was expressed in Escherichia coli cells and the recombinant product was purified to homogeneity. Gel mobility-shift assays and DNase I footprinting analyses indicated that MobR binds as a homodimer to an imperfect inverted repeat within the mobA-mobR intergenic region, with an apparent dissociation constant of 11.5(+/- 0.5) nM. The operator site is located between the start codon and the promoter region for mobA, suggesting that MobR functions as a transcriptional repressor for mobA expression. The results of effector-binding assays indicated that MobR, but not its isomers 4-hydroxybenzoate and salicylate, is released from the operator site by the addition of 3-hydroxybenzoate. This dissociation process is highly cooperative, with a Hill coefficient of approximately 2. In addition, CD spectroscopic studies demonstrated that MobR adopts two conformational states corresponding to the effector-bound and unbound forms. These results suggest that the MobR dimer possesses at least two effector-binding sites, and that the effector binding to MobR induces an allosteric conformational change required for dissociation of the protein-DNA complex.

Conformational dynamics of the isoalloxazine in substrate-free p-hydroxybenzoate hydroxylase: single-molecule studies.

p-Hydroxybenzoate hydroxylase (PHBH) is a homodimeric enzyme in which each subunit noncovalently binds one molecule of FAD in the active site. PHBH is a model system for how flavoenzymes regulate reactions with oxygen. We report single-molecule fluorescence studies of PHBH in the absence of substrate that provide data consistent with the hypothesis that a critical step in substrate binding is the movement of the isoalloxazine between an "in" conformation and a more exposed or "open" conformation. The isoalloxazine is observed to move between these conformations in the absence of substrate. Studies with the Y222A mutant form of PHBH suggest that the exposed conformation is fluorescent while the in-conformation is quenched. Finally, we note that many of the single-molecule-fluorescence trajectories reveal a conformational heterogeneity, with populations of the enzyme characterized by either fast or slow switching between the in- and open-conformations. Our data also allow us to hypothesize a model in which one flavin in the dimer inhibits the motion of the other.

Properties of p-hydroxybenzoate hydroxylase when stabilized in its open conformation.

p-Hydroxybenzoate hydroxylase is extensively studied as a model for single-component flavoprotein monooxygenases. It catalyzes a reaction in two parts: (1) reduction of the FAD in the enzyme by NADPH in response to binding of p-hydroxybenzoate to the enzyme and (2) oxidation of reduced FAD with oxygen in an environment free from solvent to form a hydroperoxide, which then reacts with p-hydroxybenzoate to form an oxygenated product. These different reactions are coordinated through conformational rearrangements of the protein and the isoalloxazine ring during catalysis. Until recently, it has not been clear how p-hydroxybenzoate gains access to the buried active site. In 2002, a structure of a mutant form of the enzyme without substrate was published that showed an open conformation with solvent access to the active site [Wang, J., Ortiz-Maldonado, M., Entsch, B., Massey, V., Ballou, D., and Gatti, D. L. (2002) Proc. Natl. Acad. Sci. U.S.A. 99, 608-613]. The wild-type enzyme does not form high-resolution crystals without substrate. We hypothesized that the wild-type enzyme without substrate also forms an open conformation for binding p-hydroxybenzoate, but only transiently. To test this idea, we have studied the properties of two different mutant forms of the enzyme that are stabilized in the open conformation. These mutant enzymes bind p-hydroxybenzoate very fast, but with very low affinity, as expected from the open structure. The mutant enzymes are extremely inactive, but are capable of slowly forming small amounts of product by the normal catalytic pathway. The lack of activity results from the failure of the mutants to readily form the out conformation required for flavin reduction by NADPH. The mutants form a large fraction of an abnormal conformation of the reduced enzyme with p-hydroxybenzoate bound. This conformation of the enzyme is unreactive with oxygen. We conclude that transient formation of this open conformation is the mechanism for sequestering p-hydroxybenzoate to initiate catalysis. This overall study emphasizes the role that protein dynamics can play in enzymatic catalysis.

Dynamics involved in catalysis by single-component and two-component flavin-dependent aromatic hydroxylases.

Flavoprotein monooxygenases are involved in a wide variety of biological processes including drug detoxification, biodegradation of aromatic compounds in the environment, biosynthesis of antibiotics and siderophores, and many others. The reactions use NAD(P)H and O2 as co-substrates and insert one atom of oxygen into the substrate. The flavin-dependent monooxygenases utilize a general cycle in which NAD(P)H reduces the flavin, and the reduced flavin reacts with O2 to form a C4a-(hydro)peroxyflavin intermediate, which is the oxygenating agent. This complicated catalytic process has diverse requirements that are difficult to be satisfied by a single site. Two general strategies have evolved to satisfy these requirements. para-Hydroxybenzoate hydroxylase, the paradigm for the single-component flavoprotein monooxygenases, is one of the most thoroughly studied of all enzymes. This enzyme undergoes significant protein and flavin dynamics during catalysis. There is an open conformation that gives access of substrate and product to solvent, and a closed or in conformation for the reaction with oxygen and the hydroxylation to occur. This closed form prevents solvent from destabilizing the hydroperoxyflavin intermediate. Finally, there is an out conformation achieved by movement of the isoalloxazine toward the solvent, which exposes its N5 for hydride delivery from NAD(P)H. The protein coordinates these dynamic events during catalysis. The second strategy uses a reductase to catalyze the reduction of the flavin and an oxygenase that uses the reduced flavin as a substrate to react with oxygen and hydroxylate the organic substrate. These two-component systems must be able to transfer reduced flavin from the reductase to the oxygenase and stabilize a C4a-peroxyflavin until a substrate binds to be hydroxylated, all before flavin oxidation and release of H2O2. Again, protein dynamics are important for catalytic success.

Protein dynamics and electrostatics in the function of p-hydroxybenzoate hydroxylase.

para-Hydroxybenzoate hydroxylase is a flavoprotein monooxygenase that catalyzes a reaction in two parts: reduction of the enzyme cofactor, FAD, by NADPH in response to binding p-hydroxybenzoate to the enzyme, then oxidation of reduced FAD by oxygen to form a hydroperoxide, which oxygenates p-hydroxybenzoate to form 3,4-dihydroxybenzoate. These diverse reactions all occur within a single polypeptide and are achieved through conformational rearrangements of the isoalloxazine ring and protein residues within the protein structure. In this review, we examine the complex dynamic behavior of the protein that enables regulated fast and specific catalysis to occur. Original research papers (principally from the past 15 years) provide the information that is used to develop a comprehensive overview of the catalytic process. Much of this information has come from detailed analysis of many specific mutants of the enzyme using rapid reaction technology, biophysical measurements, and high-resolution structures obtained by X-ray crystallography. We describe how three conformations of the enzyme provide a foundation for the catalytic cycle. One conformation has a closed active site for the conduct of the oxygen reactions, which must occur in the absence of solvent. The second conformation has a partly open active site for exchange of substrate and product, and the third conformation has a closed protein structure with the isoalloxazine ring rotated out to the surface for reaction with NADPH, which binds in a surface cleft. A fundamental feature of the enzyme is a H-bond network that connects the phenolic group of the substrate in the buried active site to the surface of the protein. This network serves to protonate and deprotonate the substrate and product in the active site to promote catalysis and regulate the coordination of conformational states for efficient catalysis.

Oxygen reactions in p-hydroxybenzoate hydroxylase utilize the H-bond network during catalysis.

para-Hydroxybenzoate hydroxylase is a flavoprotein monooxygenase that catalyses a reaction in two parts: reduction of the flavin adenine dinucleotide (FAD) in the enzyme by reduced nicotinamide adenine dinucleotide phosphate (NADPH) in response to binding p-hydroxybenzoate to the enzyme and oxidation of reduced FAD with oxygen to form a hydroperoxide, which then oxygenates p-hydroxybenzoate. These different reactions are coordinated through conformational rearrangements of the protein and isoalloxazine ring during catalysis. Earlier research showed that reduction of FAD occurs when the isoalloxazine of the FAD moves to the surface of the protein to allow hydride transfer from NADPH. This move is coordinated with protein rearrangements that are triggered by deprotonation of buried p-hydroxybenzoate through a H-bond network that leads to the surface of the protein. In this paper, we examine the involvement of this same H-bond network in the oxygen reactions-the initial formation of a flavin-C4a-hydroperoxide from the reaction between oxygen and reduced flavin, the electrophilic attack of the hydroperoxide upon the substrate to form product, and the elimination of water from the flavin-C4a-hydroxide to form oxidized enzyme in association with product release. These reactions were measured through absorbance and fluorescence changes in the FAD during the reactions. Results were collected over a range of pH for the reactions of wild-type enzyme and a series of mutant enzymes with the natural substrate and substrate analogues. We discovered that the rate of formation of the flavin hydroperoxide is not influenced by pH change, which indicates that the proton required for this reaction does not come from the H-bond network. The rate of the hydroxylation reaction increases with pH in a manner consistent with a pK(a) of 7.1. We conclude that the H-bond network abstracts the phenolic proton from p-hydroxybenzoate in the transition state of oxygen transfer. The rate of formation of oxidized enzyme increases with pH in a manner consistent with a pK(a) of 7.1, indicating the involvement of the H-bond network. We conclude that product deprotonation enhances the rate of a specific conformational change required for both product release and the elimination of water from C4a-OH-FAD.

Model studies on p-hydroxybenzoate hydroxylase. The catalytic role of Arg-214 and Tyr-201 in the hydroxylation step.

A model C-(4a)-flavinhydroperoxide (FlHOOH) is described that contains the tricyclic isoalloxazine moiety, the C-4a-hydroperoxide functionality, and a beta-hydroxyethyl group to model the effect of the 2'-OH group of the ribityl side chain of native FADHOOH. The electronic structures of this reduced flavin (H(3)()Fl(red)()), its N1 anion (H(2)()Fl(red)()(-)()), oxidized flavin (HFl(ox)()), and FlHOOH have been fully optimized at the B3LYP/ 6-31+G(d,p) level of theory. This model C-4a-flavinhydroperoxide is used to describe the transition state for the key step in the paradigm aromatic hydroxylase, p-hydroxybenzoate hydroxylase (PHBH): the oxidation of p-hydroxybenzoate (p-OHB). The Tyrosine-201 residue in PHBH is modeled by phenol, and Arginine-214 is modeled by guanidine. Electrophilic aromatic substitution proceeds by an S(N)2-like attack of the aromatic sextet of p-OHB phenolate anion on the distal oxygen of FlHOOH 3. The transition structure for oxygen atom transfer is fully optimized [B3LYP/6-31+G(d,p)] and has a classical activation barrier of 24.9 kcal/mol. These data suggest that the role of the Tyr-201 is to orient the p-OHB substrate and to properly align it for the oxygen transfer step. Although the negatively charged phenolate oxygen does activate the C-3 carbon of p-OHB phenolate anion toward oxidation relative to ortho oxidation of the carboxylate anion, it appears that H-bonding the Tyr-201 residue to this phenolic oxygen stabilizes both the ground state (GS) and the transition state (TS) approximately equally and therefore plays only a minor role, if any, in lowering the activation barrier. Complexation of p-OHB with guanidine has only a modest effect upon the oxidation barriers. When the complex is in the form of a salt-bridge (10a), the barrier is only slightly reduced. When the TSs are placed in THF solvent (COSMO) with full geometry optimization, salt-bridge TS-A is slightly favored (DeltaDeltaE() = 2.3 kcal/mol).