P-hydroxybenzoate Hydroxylase Research Articles

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.

Toward QM/MM Simulation of Enzymatic Reactions with the Drude Oscillator Polarizable Force Field.

The polarization of the environment can influence the results from hybrid quantum mechanical/molecular mechanical (QM/MM) simulations of enzymatic reactions. In this article, we address several technical aspects in the development of polarizable QM/MM embedding using the Drude Oscillator (DO) force field. We propose a stable and converging update of the DO polarization state for geometry optimizations and a suitable treatment of the QM/MM-DO boundary when the QM and MM regions are separated by cutting through a covalent bond. We assess the performance of our approach by computing binding energies and geometries of three selected complexes relevant to biomolecular modeling, namely the water trimer, the N-methylacetamide dimer, and the cationic bis(benzene)sodium sandwich complex. Using a recently published MM-DO force field for proteins, we evaluate the effect of MM polarization on the QM/MM energy profiles of the enzymatic reactions catalyzed by chorismate mutase and by p-hydroxybenzoate hydroxylase. We find that inclusion of MM polarization affects the computed barriers by about 10%.

Amperometric bienzyme screen-printed biosensor for the determination of leucine

Leucine plays an important role in protein synthesis, brain functions, building muscle mass, and helping the body when it undergoes stress. Here, we report a new amperometric bienzyme screen-printed biosensor for the determination of leucine, by coimmobilizing p-hydroxybenzoate hydroxylase (HBH) and leucine dehydrogenase (LDH) on a screen-printed electrode with NADP(+) and p-hydroxybenzoate as the cofactors. The detection principle of the sensor is that LDH catalyzes the specific dehydrogenation of leucine by using NADP(+) as a cofactor. The product, NADPH, triggers the hydroxylation of p-hydroxybenzoate by HBH in the presence of oxygen to produce 3,4-dihydroxybenzoate, which results in a change in electron concentration at the working carbon electrode, which is detected by the potentiostat. The sensor shows a linear detection range between 10 and 600μM with a detection limit of 2μM. The response is reproducible and has a fast measuring time of 5-10s after the addition of a given concentration of leucine.

Crystal Structure of 3-Hydroxybenzoate 6-Hydroxylase Uncovers Lipid-assisted Flavoprotein Strategy for Regioselective Aromatic Hydroxylation

3-Hydroxybenzoate 6-hydroxylase (3HB6H) from Rhodococcus jostii RHA1 is a dimeric flavoprotein that catalyzes the NADH- and oxygen-dependent para-hydroxylation of 3-hydroxybenzoate to 2,5-dihydroxybenzoate. In this study, we report the crystal structure of 3HB6H as expressed in Escherichia coli. The overall fold of 3HB6H is similar to that of p-hydroxybenzoate hydroxylase and other flavoprotein aromatic hydroxylases. Unexpectedly, a lipid ligand is bound to each 3HB6H monomer. Mass spectral analysis identified the ligand as a mixture of phosphatidylglycerol and phosphatidylethanolamine. The fatty acid chains occupy hydrophobic channels that deeply penetrate into the interior of the substrate-binding domain of each subunit, whereas the hydrophilic part is exposed on the protein surface, connecting the dimerization domains via a few interactions. Most remarkably, the terminal part of a phospholipid acyl chain is directly involved in the substrate-binding site. Co-crystallized chloride ion and the crystal structure of the H213S variant with bound 3-hydroxybenzoate provide hints about oxygen activation and substrate hydroxylation. Essential roles are played by His-213 in catalysis and Tyr-105 in substrate binding. This phospholipid-assisted strategy to control regioselective aromatic hydroxylation is of relevance for optimization of flavin-dependent biocatalysts.

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.

Is Renalase a Novel Player in Catecholaminergic Signaling? The Mystery of the Catalytic Activity of an Intriguing New Flavoenzyme

Renalase is a flavoprotein recently discovered in humans, preferentially expressed in the proximal tubules of the kidney and secreted in blood and urine. It is highly conserved in vertebrates, with homologs identified in eukaryotic and prokaryotic organisms. Several genetic, epidemiological, clinical and experimental studies show that renalase plays a role in the modulation of the functions of the cardiovascular system, being particularly active in decreasing the catecholaminergic tone, in lowering blood pressure and in exerting a protective action against myocardial ischemic damage. Deficient renalase synthesis might be the cause of the high occurrence of hypertension and adverse cardiac events in kidney disease patients. Very recently, recombinant human renalase has been structurally and functionally characterized in vitro. Results show that it belongs to the p-hydroxybenzoate hydroxylase structural family of flavoenzymes, contains non-covalently bound FAD with redox features suggestive of a dehydrogenase activity, and is not a catecholamine-degrading enzyme,either through oxidase or NAD(P)H-dependent monooxygenase reactions. The biochemical data now available will hopefully provide the basis for a systematic and rational quest toward the identification of the reaction catalyzed by renalase and of the molecular mechanism of its physiological action, which in turn are expected to favor the development of novel therapeutic tools for the treatment of kidney and cardiovascular diseases.

Expression and bioconversion of recombinant m- and p-hydroxybenzoate hydroxylases from a novel moderate halophile, Chromohalobacter sp.

p-Hydroxybenzoate hydroxylase (pobA) and m-hydroxybenzoate hydroxylase (mobA) genes, from the moderate halophile Chromohalobacter sp. HS-2, were expressed and characterized. Solubilities of overexpressed recombinant MobA and PobA were enhanced by the induction of the heat-shock proteins DnaJ and DnaK. Each MobA and PobA maintained stable activity under high NaCl concentrations. V (max) and K (m) values for MobA with m-hydroxybenzoate were 70 μmol min(-1) mg(-1) protein and 81 μM, respectively. Similarly, those of PobA with p-hydroxybenzoate as substrate were 5 μmol min(-1) mg(-1) protein and 129 μM, respectively. The Escherichia coli expression system, including induction of heat shock proteins, was used to convert hydroxybenzoates into protocatechuate (3,4-dihydroxybenzoate) and revealed that resting cells harboring mobA converted 15 mM m-hydroxybenzoate to 15 mM protocatechuate while those harboring pobA converted 50 mM p-hydroxybenzoate to 35 mM protocatechuate at 30 °C, respectively.

ChemInform Abstract: Catalytic and Structural Features of Flavoprotein Hydroxylases and Epoxidases.

Monooxygenases perform chemo-, regio- and/or enantioselective oxygenations of organic substrates under mild reaction conditions. These properties and the increasing number of representatives along with effective preparation methods place monooxygenases in the focus of industrial biocatalysis. Mechanistic and structural insights reveal reaction sequences and allow turning them into efficient tools for the production of valuable products. Herein we describe two biocatalytically relevant subclasses of flavoprotein monooxygenases with a close evolutionary relation: subclass A represented by p-hydroxybenzoate hydroxylase (PHBH) and subclass E formed by styrene monooxygenases (SMOs). PHBH family members perform highly regioselective hydroxylations on a wide variety of aromatic compounds. The more recently discovered SMOs catalyze a number of stereoselective epoxidation and sulfoxidation reactions. Mechanistic and structural studies expose distinct characteristics, which provide a promising source for future biocatalyst development.

Role of the Somersault Rearrangement in the Oxidation Step for Flavin Monooxygenases (FMO). A Comparison between FMO and Conventional Xenobiotic Oxidation with Hydroperoxides

Model quantum mechanical calculations presented for C-4a-flavin hydroperoxide (FlHOOH) at the B3LYP/6-311+G(d,p) level suggest a new mechanism for flavoprotein monooxygenase (FMO) oxidation involving a concerted homolytic O-O bond cleavage in concert with hydroxyl radical transfer from the flavin hydroperoxide rather than an S(N)2-like displacement by the substrate on the C-4a-hydroperoxide OOH group. Homolytic O-O bond cleavage in a somersault-like rearrangement of hydroperoxide C-4a-flavinhydroperoxide (1) (FLHO-OH → FLHO···HO) produces an internally hydrogen-bonded HO(•) radical intermediate with a classical activation barrier of 27.0 kcal/mol. Model hydroperoxide 1 is used to describe the transition state for the key oxidation step in the paradigm aromatic hydroxylase, p-hydroxybenzoate hydroxylase (PHBH). A comparison of the electron distribution in the transition structures for the PHBH hydroxylation of p-hydroxybenzoic acid (ΔE(‡) = 23.0 kcal/mol) with that of oxidation of trimethylamine (ΔE(‡) = 22.3 kcal/mol) and dimethyl sulfide (ΔE‡ = 14.1 kcal/mol) also suggests a mechanism involving a somersault mechanism in concert with transfer of an HO(•) radical to the nucleophilic heteroatom center with a hydrogen transfer back to the FLH-O residue after the barrier is crossed to produce the final product, FLH-OH. In each case the hydroxylation barrier was less than that of the O-O rearrangement barrier in the absence of a substrate supporting an overall concerted process. All three transition structures bear a resemblance to the TS for the comparable hydroxylation of isobutane (ΔE(‡) = 29.2 kcal/mol) and for simple Fenton oxidation by aqueous iron(III) hydroperoxides. To our surprise the oxidation of N- and S-nucleophiles with conventional oxidants such as alkyl hydroperoxides and peracids also proceeds by HO(•) radical transfer in a manner quite similar to that for tricyclic hydroperoxide 1. Stabilization of the developing oxyradical produced by somersault rearrangement for concerted enzymatic oxidation with tricyclic hydroperoxide 1 results in a reduced overall activation barrier.

Catalytic and Structural Features of Flavoprotein Hydroxylases and Epoxidases

AbstractMonooxygenases perform chemo‐, regio‐ and/or enantioselective oxygenations of organic substrates under mild reaction conditions. These properties and the increasing number of representatives along with effective preparation methods place monooxygenases in the focus of industrial biocatalysis. Mechanistic and structural insights reveal reaction sequences and allow turning them into efficient tools for the production of valuable products. Herein we describe two biocatalytically relevant subclasses of flavoprotein monooxygenases with a close evolutionary relation: subclass A represented by p‐hydroxybenzoate hydroxylase (PHBH) and subclass E formed by styrene monooxygenases (SMOs). PHBH family members perform highly regioselective hydroxylations on a wide variety of aromatic compounds. The more recently discovered SMOs catalyze a number of stereoselective epoxidation and sulfoxidation reactions. Mechanistic and structural studies expose distinct characteristics, which provide a promising source for future biocatalyst development.

Kinetic and spectroscopic characterization of the putative monooxygenase domain of human MICAL-1

MICALs form a conserved multidomain protein family essential for cytoskeletal rearrangements. To complement structural information available, we produced the FAD-containing monooxygenase-like domain of human MICAL-1 (MICAL-MO) in forms differing for the presence and location of a His-tag, which only influences the protein yields. The Km for NADPH of the NADPH oxidase reaction is sensitive to ionic strength and type of ions. The apparent kcat (pH 7) is limited by enzyme reduction by NADPH, which occurs without detectable intermediates, as established by anaerobic rapid reaction experiments. The sensitivity to ionic strength and type of ions and the pH dependence of the steady-state kinetic parameters extend MICAL-MO similarity with enzymes of the p-hydroxybenzoate hydroxylase class at the functional level. The reaction is also sensitive to solvent viscosity, providing a tool to monitor the conformational changes predicted to occur during turnover. Finally, it was confirmed that MICAL-MO promotes actin depolymerization, and it was shown that F-actin, but not G-actin, stimulates NADPH oxidation by increasing kcat and kcat/KNADPH (≈5 and ≈200-fold, respectively) with an apparent Km for actin of 4.7μM, under conditions that stabilize F-actin. The time-course of NADPH oxidation shows substrate recycling, indicating the possible reversibility of MICAL effect.

FAD-Binding Site and NADP Reactivity in Human Renalase: A New Enzyme Involved in Blood Pressure Regulation

Renalase is a recently discovered flavoprotein that regulates blood pressure, regulates sodium and phosphate excretion, and displays cardioprotectant action through a mechanism that is barely understood to date. It has been proposed to act as a catecholamine-degrading enzyme, via either O 2-dependent or NADH-dependent mechanisms. Here we report the renalase crystal structure at 2.5 Å resolution together with new data on its interaction with nicotinamide dinucleotides. Renalase adopts the p-hydroxybenzoate hydroxylase fold topology, comprising a Rossmann-fold-based flavin adenine dinucleotide (FAD)-binding domain and a putative substrate-binding domain, the latter of which contains a five-stranded anti-parallel β-sheet. A large cavity (228 Å 3), facing the flavin ring, presumably represents the active site. Compared to monoamine oxidase or polyamine oxidase, the renalase active site is fully solvent exposed and lacks an ‘aromatic cage’ for binding the substrate amino group. Renalase has an extremely low diaphorase activity, displaying lower k cat but higher k cat/ K m for NADH compared to NADPH. Moreover, its FAD prosthetic group becomes slowly reduced when it is incubated with NADPH under anaerobiosis, and binds NAD + or NADP + with K d values of ca 2 mM. The absence of a recognizable NADP-binding site in the protein structure and its poor affinity for, and poor reactivity towards, NADH and NADPH suggest that these are not physiological ligands of renalase. Although our study does not answer the question on the catalytic activity of renalase, it provides a firm framework for testing hypotheses on the molecular mechanism of its action.

Structure and Mutation Analysis of Archaeal Geranylgeranyl Reductase

The crystal structure of geranylgeranyl reductase (GGR) from Sulfolobus acidocaldarius was determined in order to elucidate the molecular mechanism of the catalytic reaction. The enzyme is a flavoprotein and is involved in saturation of the double bonds on the isoprenoid moiety of archaeal membranes. The structure determined in this study belongs to the p-hydroxybenzoate hydroxylase family in the glutathione reductase superfamily. GGR functions as a monomer and is divided into the FAD-binding, catalytic and C-terminal domains. The catalytic domain has a large cavity surrounded by a characteristic YxWxFPx 7-8GxG motif and by the isoalloxazine ring of an FAD molecule. The cavity holds a lipid molecule, which is probably derived from Escherichia coli cells used for over-expression. One of the two forms of the structure clarifies the presence of an anion pocket holding a pyrophosphate molecule, which might anchor the phosphate head of the natural ligands. Mutational analysis supports the suggestion that the three aromatic residues of the YxWxFPx 7-8GxG motif hold the ligand in the appropriate position for reduction. Cys47, which is widely conserved in GGRs, is located at the si-side of the isoalloxazine ring of FAD and is shown by mutational analysis to be involved in catalysis. The catalytic cycle, including the FAD reducing factor binding site, is proposed on the basis of the detailed analysis of the structure.

Long-Range Electrostatic Effects in QM/MM Studies of Enzymatic Reactions: Application of the Solvated Macromolecule Boundary Potential.

Long-range electrostatic interactions are important in simulations of enzymatic reactions. They can be divided into the effects due to bulk solvent and those due to the electrostatic potential of the outer macromolecule. We study and quantify the importance of these two effects for two test systems by application of the solvated macromolecule boundary potential (SMBP) [J. Chem. Theory Comput. 2009, 5, 3114-3128]. We validate the accuracy of the SMBP for these test systems and present a transferable protocol for determination of optimal SMBP parameters as well as recommended default values for these parameters. Two enzymatic reactions with different characteristics are studied: the intramolecular Claisen rearrangement in chorismate mutase that is associated with little charge transfer and the hydroxylation reaction in p-hydroxybenzoate hydroxylase that corresponds to a formal "OH(+)" transfer and thus involves significant charge transfer. It is found that the effects of the electrostatic potential of the outer macromolecule and of bulk solvent are only important in the latter case, where their neglect causes deviations in the computed barriers on the order of 1-2 kcal/mol, respectively. Even larger deviations on the order of several kilocalories per mole are observed for the reaction energies in p-hydroxybenzoate hydroxylase if the electrostatic potential of the outer macromolecule is neglected.

Insights into Substrate Specificity of Geranylgeranyl Reductases Revealed by the Structure of Digeranylgeranylglycerophospholipid Reductase, an Essential Enzyme in the Biosynthesis of Archaeal Membrane Lipids

Archaeal membrane lipids consist of branched, saturated hydrocarbons distinct from those found in bacteria and eukaryotes. Digeranylgeranylglycerophospholipid reductase (DGGR) catalyzes the hydrogenation process that converts unsaturated 2,3-di- O-geranylgeranylglyceryl phosphate to saturated 2,3-di- O-phytanylglyceryl phosphate as a critical step in the biosynthesis of archaeal membrane lipids. The saturation of hydrocarbon chains confers the ability to resist hydrolysis and oxidation and helps archaea withstand extreme conditions. DGGR is a member of the geranylgeranyl reductase family that is also widely distributed in bacteria and plants, where the family members are involved in the biosynthesis of photosynthetic pigments. We have determined the crystal structure of DGGR from the thermophilic heterotrophic archaea Thermoplasma acidophilum at 1.6 Å resolution, in complex with flavin adenine dinucleotide (FAD) and a bacterial lipid. The DGGR structure can be assigned to the well-studied, p-hydroxybenzoate hydroxylase (PHBH) SCOP superfamily of flavoproteins that include many aromatic hydroxylases and other enzymes with diverse functions. In the DGGR complex, FAD adopts the IN conformation (closed) previously observed in other PHBH flavoproteins. DGGR contains a large substrate-binding site that extends across the entire ligand-binding domain. Electron density corresponding to a bacterial lipid was found within this cavity. The cavity consists of a large opening that tapers down to two, narrow, curved tunnels that closely mimic the shape of the preferred substrate. We identified a sequence motif, PxxYxWxFP, that defines a specificity pocket in the enzyme and precisely aligns the double bond of the geranyl group with respect to the FAD cofactor, thus providing a structural basis for the substrate specificity of geranylgeranyl reductases. DGGR is likely to share a common mechanism with other PHBH enzymes in which FAD switches between two conformations that correspond to the reductive and oxidative half cycles. The structure provides evidence that substrate binding likely involves conformational changes, which are coupled to the two conformational states of the FAD.

Amperometric ATP Biosensors Based on Coimmobilizations of p-Hydroxybenzoate Hydroxylase, Glucose-6-Phosphate Dehydrogenase, and Hexokinase on Clark-Type and Screen-Printed Electrodes

Two types of amperometric trienzyme ATP biosensors were developed based on new combinations of enzymes and electrodes by using the coimmobilizations of p-hydroxybenzoate hydroxylase (HBH), glucose-6-phosphate dehydrogenase (G6PDH), and hexokinase (HEX) on a Clark-type oxygen electrode and on a screen-printed electrode. The principles of the determinations are as follows: HEX transfers the phosphate group from ATP to glucose to form glucose-6-phosphate. G6PDH catalyzes the specific dehydrogenation of glucose-6-phosphate by consuming NADP+. The product, NADPH initiates the irreversible hydroxylation of p-hydroxybenzoate by HBH to consume dissolved oxygen which results in a detectable signal on a Clark-type electrode and generate 3, 4-dihydroxybenzoate which results in a detectable signal on a screen-printed electrode. Both sensors show high-performance characteristics with broad detection ranges, short measuring times and good specificities.

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.

A General Boundary Potential for Hybrid QM/MM Simulations of Solvated Biomolecular Systems

We present a general boundary potential for the efficient and accurate evaluation of electrostatic interactions in hybrid quantum mechanical/molecular mechanical (QM/MM) approaches called solvated macromolecule boundary potential (SMBP), which is designed for QM/MM calculations with any kind of QM method. The SMBP targets QM/MM single-point energy calculations and geometry optimizations. In the SMBP scheme, the outer solvent and macromolecule region is described by a boundary potential obtained with the use of Poisson-Boltzmann calculations (treating the bulk solvent as a dielectric continuum). In the QM calculations, the SMBP is represented by virtual point charges on a surface enclosing the explicitly treated inner region. These charges and their interactions with the QM density are determined through a self-consistent reaction field procedure. The accuracy of the SMBP is evaluated on three diverse test systems: the intramolecular proton transfer of glycine in water, the hydroxylation reaction in p-hydroxybenzoate hydroxylase, and the spin state energy splittings in the pentacoordinated ferric complex of cytochrome P450cam. In the case of solvated glycine, application of the SMBP turns out to be problematic since analogous QM/MM/SMBP and full QM/MM geometry optimizations lead to different close-lying local minima. In both enzymes, the SMBP performs very well and closely reproduces the results from full QM/MM optimizations of these more rigid test systems. Starting from optimized QM/MM/SMBP structures along a reaction path, one can apply the previously implemented generalized solvent boundary potential (GSBP) to sample over MM phase space in QM/MM free energy calculations within the framework of free energy perturbation theory. This reduces the overall computational costs of sampling by 1 order of magnitude while maintaining good accuracy. The combined use of SMBP and GSBP thus allows for efficient QM/MM free energy studies of enzymes.

Structural Basis for Substrate Recognition and Specificity in Aklavinone-11-Hydroxylase from Rhodomycin Biosynthesis

In the biosynthesis of several anthracyclines, aromatic polyketides produced by many Streptomyces species, the aglycone core is modified by a specific flavin adenine dinucleotide (FAD)- and NAD(P)H-dependent aklavinone-11-hydroxylase. Here, we report the crystal structure of a ternary complex of this enzyme from Streptomyces purpurascens, RdmE, with FAD and the substrate aklavinone. The enzyme is built up of three domains, a FAD-binding domain, a domain involved in substrate binding, and a C-terminal thioredoxin-like domain of unknown function. RdmE exhibits structural similarity to aromatic hydroxylases from the p-hydroxybenzoate hydroxylase family, but unlike most other related enzymes, RdmE is a monomer. The substrate is bound in a hydrophobic pocket in the interior of the enzyme, and access to this pocket is provided through a different route than for the isoalloxazine ring of FAD—the backside of the ligand binding cleft. The architecture of the substrate binding pocket and the observed enzyme–aklavinone interactions provide a structural explanation for the specificity of the enzyme for non-glycosylated substrates with C9-R stereochemistry. The isoalloxazine ring of the flavin cofactor is bound in the “out” conformation but can be modeled in the “in” conformation without invoking large conformational changes of the enzyme. This model places the flavin ring in a position suitable for catalysis, almost perpendicular to the tetracyclic ring system of the substrate and with a distance of the C4a carbon atom of the isoalloxazine ring to the C-11 carbon atom of the substrate of 4.8 Å. The structure suggested that a Tyr224-Arg373 pair might be involved in proton abstraction at the C-6 hydroxyl group, thereby increasing the nucleophilicity of the aromatic ring system and facilitating electrophilic attack by the perhydroxy-flavin intermediate. Replacement of Tyr224 by phenylalanine results in inactive enzyme, whereas mutants at position Arg373 retain catalytic activity close to wild-type level. These data establish an essential role of residue Tyr224 in catalysis, possibly in aligning the substrate in a position suitable for catalysis.

Kinetic Mechanism of Ornithine Hydroxylase (PvdA) from Pseudomonas aeruginosa: Substrate Triggering of O2 Addition but Not Flavin Reduction

PvdA catalyzes the hydroxylation of the side chain primary amine of ornithine in the initial step of the biosynthesis of the Pseudomonas aeruginosa siderophore pyoverdin. The reaction requires FAD, NADPH, and O(2). PvdA uses the same cosubstrates as several flavin-dependent hydroxylases that differ one from another in the kinetic mechanisms of their oxidative and reductive half-reactions. Therefore, the mechanism of PvdA was determined by absorption stopped-flow experiments. By contrast to some flavin-dependent hydroxylases (notably, p-hydroxybenzoate hydroxylase), binding of the hydroxylation target is not required to trigger reduction of the flavin by NADPH: the reductive half-reaction is equally facile in the presence and absence of ornithine. Reaction of O(2) with FADH(2) in the oxidative half-reaction is accelerated by ornithine 80-fold, providing a mechanism by which PvdA can ensure coupling of NADPH and ornithine oxidation. In the presence of ornithine, the expected C(4a)-hydroperoxyflavin intermediate with 390 nm absorption accumulates and decays to the C(4a)-hydroxyflavin in a kinetically competent fashion. The slower oxidative half-reaction that occurs in the absence of ornithine involves accumulation of an oxygenated flavin species and two subsequent states that are tentatively assigned as C(4a)-peroxy- and C(4a)-hydroperoxyflavin intermediates and the oxidized flavin. The enzyme generates stoichiometric hydrogen peroxide in lieu of hydroxyornithine. The data suggest that PvdA employs a kinetic mechanism that is a hybrid of those previously documented for other flavin-dependent hydroxylases.