Reduction kinetics of 3-hydroxybenzoate 6-hydroxylase from Rhodococcus jostii RHA1

Mechanism of the Multistep Catalytic Cycle of 6-Hydroxynicotinate 3-Monooxygenase Revealed by Global Kinetic Analysis

The class A flavoenzyme 6-hydroxynicotinate 3-monooxygenase (NicC) catalyzes a rare decarboxylative hydroxylation reaction in the degradation of nicotinate by aerobic bacteria. While the structure and critical residues involved in catalysis have been reported, the mechanism of this multistep enzyme has yet to be determined. A kinetic understanding of the NicC mechanism would enable comparison to other phenolic hydroxylases and illuminate its bioengineering potential for remediation of N-heterocyclic aromatic compounds. Toward these goals, transient state kinetic analyses by stopped-flow spectrophotometry were utilized to follow rapid changes in flavoenzyme absorbance spectra during all three stages of NicC catalysis: (1) 6-HNA binding; (2) NADH binding and FAD reduction; and (3) O₂ binding with C4a-adduct formation, substrate hydroxylation, and FAD regeneration. Global kinetic simulations by numeric integration were used to supplement analytical fitting of time-resolved data and establish a kinetic mechanism. Results indicate that 6-HNA binding is a two-step process that substantially increases the affinity of NicC for NADH and enables the formation of a charge-transfer-complex intermediate to enhance the rate of flavin reduction. Singular value decomposition of the time-resolved spectra during the reaction of the substrate-bound, reduced enzyme with dioxygen provides evidence for the involvement of C4a-hydroperoxy-flavin and C4a-hydroxy-flavin intermediates in NicC catalysis. Global analysis of the full kinetic mechanism suggests that steady-state catalytic turnover is partially limited by substrate hydroxylation and C4a-hydroxy-flavin dehydration to regenerate the flavoenzyme. Insights gleaned from the kinetic model and determined microscopic rate constants provide a fundamental basis for understanding NicC's substrate specificity and reactivity.

Phenolic hydroxylases

Many flavin-dependent phenolic hydroxylases (monooxygenases) have been extensively investigated. Their crystal structures and reaction mechanisms are well understood. These enzymes belong to groups A and D of the flavin-dependent monooxygenases and can be classified as single-component and two-component flavin-dependent monooxygenases. The insertion of molecular oxygen into the substrates catalyzed by these enzymes is beneficial for modifying the biological properties of phenolic compounds and their derivatives. This chapter provides an in-depth discussion of the structural features of single-component and two-component flavin-dependent phenolic hydroxylases. The reaction mechanisms of selected enzymes, including 3-hydroxy-benzoate 4-hydroxylase (PHBH) and 3-hydroxy-benzoate 6-hydroxylase as representatives of single-component enzymes and 3-hydroxyphenylacetate 4-hydroxylase (HPAH) as a representative of two-component enzymes, are discussed in detail. This chapter comprises the following four main parts: general reaction, structures, reaction mechanisms, and enzyme engineering for biocatalytic applications. Enzymes belonging to the same group catalyze similar reactions but have different unique structural features to control their reactivity to substrates and the formation and stabilization of C4a-hydroperoxyflavin. Protein engineering has been employed to improve the ability to use these enzymes to synthesize valuable compounds. A thorough understanding of the structural and mechanistic features controlling enzyme reactivity is useful for enzyme redesign and enzyme engineering for future biocatalytic applications.

Mechanisms of ligand binding

Many processes in chemistry and biology involve interactions of a ligand with its molecular target. Interest in the mechanism governing such interactions has dominated theoretical and experimental analysis for over a century. The interpretation of molecular recognition has evolved from a simple rigid body association of the ligand with its target to appreciation of the key role played by conformational transitions. Two conceptually distinct descriptions have had a profound impact on our understanding of mechanisms of ligand binding. The first description, referred to as induced fit, assumes that conformational changes follow the initial binding step to optimize the complex between the ligand and its target. The second description, referred to as conformational selection, assumes that the free target exists in multiple conformations in equilibrium and that the ligand selects the optimal one for binding. Both descriptions can be merged into more complex reaction schemes that better describe the functional repertoire of macromolecular systems. This review deals with basic mechanisms of ligand binding, with special emphasis on induced fit, conformational selection, and their mathematical foundations to provide rigorous context for the analysis and interpretation of experimental data. We show that conformational selection is a surprisingly versatile mechanism that includes induced fit as a mathematical special case and even captures kinetic properties of more complex reaction schemes. These features make conformational selection a dominant mechanism of molecular recognition in biology, consistent with the rich conformational landscape accessible to biological macromolecules being unraveled by structural biology.

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 ¹⁹F NMR in conjunction with fluorinated substrate analogs to directly measure pK_a values and to monitor protein events in hydroxylase active sites. We found that the single-component monooxygenase 3-hydroxybenzoate 6-hydroxylase (3HB6H) depresses the pK_a of the bound substrate analog 4-fluoro-3-hydroxybenzoate (4F3HB) by 1.6 pH units, consistent with previously proposed mechanisms. ¹⁹F NMR was applied anaerobically to the two-component monooxygenase 4-hydroxyphenylacetate 3-hydroxylase (HPAH), revealing depression of the pK_a of 3-fluoro-4-hydroxyphenylacetate by 2.5 pH units upon binding to the C₂ component of HPAH. ¹⁹F NMR also revealed a pK_a 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-C₂, respectively, are consistent with pK_a tuning by one or more H-bonding interactions. Our implementation of ¹⁹F NMR in anaerobic samples is applicable to other two-component flavin-dependent hydroxylases and promises to expand our understanding of their catalytic mechanisms.

3-Hydroxybenzoate 6-Hydroxylase from Rhodococcus jostii RHA1 Contains a Phosphatidylinositol Cofactor

3-Hydroxybenzoate 6-hydroxylase (3HB6H, EC 1.13.14.26) is a FAD-dependent monooxygenase involved in the catabolism of aromatic compounds in soil microorganisms. 3HB6H is unique among flavoprotein hydroxylases in that it harbors a phospholipid ligand. The purified protein obtained from expressing the gene encoding 3HB6H from Rhodococcus jostii RHA1 in the host Escherichia coli contains a mixture of phosphatidylglycerol and phosphatidylethanolamine, which are the major constituents of E. coli's cytoplasmic membrane. Here, we purified 3HB6H (RjHB6H) produced in the host R. jostii RHA#2 by employing a newly developed actinomycete expression system. Biochemical and biophysical analysis revealed that Rj3HB6H possesses similar catalytic and structural features as 3HB6H, but now contains phosphatidylinositol, which is a specific constituent of actinomycete membranes. Native mass spectrometry suggests that the lipid cofactor stabilizes monomer-monomer contact. Lipid analysis of 3HB6H from Pseudomonas alcaligenes NCIMB 9867 (Pa3HB6H) produced in E. coli supports the conclusion that 3HB6H enzymes have an intrinsic ability to bind phospholipids with different specificity, reflecting the membrane composition of their bacterial host.

Reaction of pyranose dehydrogenase from Agaricus meleagris with its carbohydrate substrates

Monomeric Agaricus meleagris pyranose dehydrogenase (AmPDH) belongs to the glucose-methanol-choline family of oxidoreductases. An FAD cofactor is covalently tethered to His103 of the enzyme. AmPDH can double oxidize various mono- and oligosaccharides at different positions (C1 to C4). To study the structure/function relationship of selected active-site residues of AmPDH pertaining to substrate (carbohydrate) turnover in more detail, several active-site variants were generated, heterologously expressed in Pichia pastoris, and characterized by biochemical, biophysical and computational means. The crystal structure of AmPDH shows two active-site histidines, both of which could take on the role as the catalytic base in the reductive half-reaction. Steady-state kinetics revealed that His512 is the only catalytic base because H512A showed a reduction in (kcat /KM )glucose by a factor of 10(5) , whereas this catalytic efficiency was reduced by two or three orders of magnitude for His556 variants (H556A, H556N). This was further corroborated by transient-state kinetics, where a comparable decrease in the reductive rate constant was observed for H556A, whereas the rate constant for the oxidative half-reaction (using benzoquinone as substrate) was increased for H556A compared to recombinant wild-type AmPDH. Steady-state kinetics furthermore indicated that Gln392, Tyr510, Val511 and His556 are important for the catalytic efficiency of PDH. Molecular dynamics (MD) simulations and free energy calculations were used to predict d-glucose oxidation sites, which were validated by GC-MS measurements. These simulations also suggest that van der Waals interactions are the main driving force for substrate recognition and binding.

Flavin dependent monooxygenases

Flavin-dependent monooxygenases catalyze a wide variety of chemo-, regio- and enantioselective oxygenation reactions. As such, they are involved in key biological processes ranging from catabolism, detoxification and biosynthesis, to light emission and axon guidance. Based on fold and function, flavin-dependent monooxygenases can be distributed into eight groups. Groups A and B comprise enzymes that rely on NAD(P)H as external electron donor. Groups C-F are two-protein systems, composed of a monooxygenase and a flavin reductase. Groups G and H comprise internal monooxygenases that reduce the flavin cofactor through substrate oxidation. Recently, many new flavin-dependent monooxygenases have been discovered. In addition to posing basic enzymological questions, these proteins attract attention of pharmaceutical and fine-chemical industries, given their importance as regio- and enantioselective biocatalysts. In this review we present an update of the classification of flavin-dependent monooxygenases and summarize the latest advances in our understanding of their catalytic and structural properties.

The reaction kinetics of 3-hydroxybenzoate 6-hydroxylase from Rhodococcus jostii RHA1 provide an understanding of the para-hydroxylation enzyme catalytic cycle

3-Hydroxybenzoate 6-hydroxylase (3HB6H) from Rhodococcus jostii RHA1 is an NADH-specific flavoprotein monooxygenase that catalyzes the para-hydroxylation of 3-hydroxybenzoate (3HB) to form 2,5-dihydroxybenzoate (2,5-DHB). Based on results from stopped-flow spectrophotometry, the reduced enzyme-3HB complex reacts with oxygen to form a C4a-peroxy flavin with a rate constant of 1.13 ± 0.01 × 10(6) m(-1) s(-1) (pH 8.0, 4 °C). This intermediate is subsequently protonated to form a C4a-hydroperoxyflavin with a rate constant of 96 ± 3 s(-1). This step shows a solvent kinetic isotope effect of 1.7. Based on rapid-quench measurements, the hydroxylation occurs with a rate constant of 36 ± 2 s(-1). 3HB6H does not exhibit substrate inhibition on the flavin oxidation step, a common characteristic found in most ortho-hydroxylation enzymes. The apparent kcat at saturating concentrations of 3HB, NADH, and oxygen is 6.49 ± 0.02 s(-1). Pre-steady state and steady-state kinetic data were used to construct the catalytic cycle of the reaction. The data indicate that the steps of product release (11.7 s(-1)) and hydroxylation (36 ± 2 s(-1)) partially control the overall turnover.

Essential role of conformational selection in ligand binding

Two competing and mutually exclusive mechanisms of ligand recognition - conformational selection and induced fit - have dominated our interpretation of ligand binding in biological macromolecules for almost six decades. Conformational selection posits the pre-existence of multiple conformations of the macromolecule from which the ligand selects the optimal one. Induced fit, on the other hand, postulates the existence of conformational rearrangements of the original conformation into an optimal one that are induced by binding of the ligand. In the former case, conformational transitions precede the binding event; in the latter, conformational changes follow the binding step. Kineticists have used a facile criterion to distinguish between the two mechanisms based on the dependence of the rate of relaxation to equilibrium, kobs, on the ligand concentration, [L]. A value of kobs decreasing hyperbolically with [L] has been seen as diagnostic of conformational selection, while a value of kobs increasing hyperbolically with [L] has been considered diagnostic of induced fit. However, this simple conclusion is only valid under the rather unrealistic assumption of conformational transitions being much slower than binding and dissociation events. In general, induced fit only produces values of kobs that increase with [L] but conformational selection is more versatile and is associated with values of kobs that increase with, decrease with or are independent of [L]. The richer repertoire of kinetic properties of conformational selection applies to kinetic mechanisms with single or multiple saturable relaxations and explains the behavior of nearly all experimental systems reported in the literature thus far. Conformational selection is always sufficient and often necessary to account for the relaxation kinetics of ligand binding to a biological macromolecule and is therefore an essential component of any binding mechanism. On the other hand, induced fit is never necessary and only sufficient in a few cases. Therefore, the long assumed importance and preponderance of induced fit as a mechanism of ligand binding should be reconsidered.

Conformational selection is a dominant mechanism of ligand binding

Molecular recognition in biological macromolecules is achieved by binding interactions coupled to conformational transitions that precede or follow the binding step, two limiting mechanisms known as conformational selection and induced fit, respectively. Sorting out the contribution of these mechanisms to any binding interaction remains a challenging task of general interest in biochemistry. Here we show that conformational selection is associated with a vast repertoire of kinetic behaviors, can never be disproved a priori as a mechanism of ligand binding, and is sufficient to explain the relaxation kinetics documented experimentally for a large number of systems. On the other hand, induced fit features a narrow spectrum of kinetic behaviors and can be disproved in many cases in which conformational selection offers the only possible explanation. This conclusion offers a paradigm shift in the analysis of relaxation kinetics, with conformational selection acquiring preeminence as a mechanism of ligand binding. The dominant role of conformational selection supports the emerging structural view of the macromolecule as a conformational ensemble from which the ligand selects the initial optimal fit to produce a biological response.