Improving Toluene Dioxygenase Activity for Ester‐Functionalized Substrates through Enzyme Engineering

Rieske dioxygenases are multi-component enzyme systems, naturally found in many soil bacteria, that have been widely applied in the production of fine chemicals, owing to the unique and valuable oxidative dearomatization reactions they catalyze. The range of practical applications for these enzymes in this context has historically been limited, however, due to their limited substrate scope and strict selectivity. In an attempt to overcome these limitations, our research group has employed the tools of enzyme engineering to expand the substrate scope or improve the reactivity of these enzyme systems in specific contexts. Traditionally, enzyme engineering campaigns targeting metalloenzymes have avoided mutations to metal-coordinating residues, based on the assumption that these residues are essential for enzyme activity. Inspired by the success of other recent enzyme engineering reports, our research group investigated the potential to alter or improve the reactivity of Rieske dioxygenases by altering or eliminating iron coordination in the active site of these enzymes. Herein, we report the modification of all three iron-coordinating residues in the active site of toluene dioxygenase both to alternate residues capable of coordinating iron, and to a residue that would eliminate iron coordination. The enzyme variants produced in this way were tested for their activity in the cis-dihydroxylation of a small library of potential aromatic substrates. The results of these studies demonstrated that all three iron-coordinating residues, in their natural state, are essential for enzyme activity in toluene dioxygenase, as the introduction of any mutations at these sites resulted in a complete loss of cis-dihydroxylation activity.

Rieske dioxygenases have a long history of being utilized as green chemical tools in the organic synthesis of high-value compounds, due to their capacity to perform the cis-dihydroxylation of a wide variety of aromatic substrates. The practical utility of these enzymes has been hampered however by steric and electronic constraints on their substrate scopes, resulting in limited reactivity with certain substrate classes. Herein, we report the engineering of a widely used member of the Rieske dioxygenase class of enzymes, toluene dioxygenase (TDO), to produce improved variants with greatly increased activity for the cis-dihydroxylation of benzoates. Through rational mutagenesis and screening, TDO variants with substantially improved activity over the wild-type enzyme were identified. Homology modeling, docking studies, molecular dynamics simulations, and substrate tunnel analysis were applied in an effort to elucidate how the identified mutations resulted in improved activity for this polar substrate class. These analyses revealed modification of the substrate tunnel as the likely cause of the improved activity observed with the best-performing enzyme variants.

Rieske dioxygenases are enzyme systems that have a long history of being applied as chiral, green chemical catalysts in the production of valuable building blocks for organic synthesis, owing to their ability to catalyze the cis‐dihydroxylation of aromatics. The practical utility of these catalysts, however, has been limited by restrictions on their substrate scope and selectivity. Recent studies have demonstrated the potential of modifying the substrate tunnel of oxidase enzymes to modulate the selectivity and activity of these enzymes for specific substrates. Herein, we report the targeted modification of residues lining the substrate tunnel of a representative and widely used Rieske dioxygenase, toluene dioxygenase (TDO). Several enzyme variants generated through modification of the residues lining the substrate tunnel demonstrated substantially improved activity over the wild‐type enzyme for multiple substrates. Homology modeling, docking studies, molecular dynamics simulations, and substrate tunnel analysis were applied in efforts to elucidate how the identified mutations resulted in improved activity. These analyses suggested that new interactions introduced along the substrate tunnel may explain the improved activity observed with the best‐performing enzyme variants.

Comparative analysis of the sequences of enzymes encoded in a variety of prokaryotic and eukaryotic genomes reveals convergence and divergence at several levels. Functional convergence can be inferred when structurally distinct and hence non-homologous enzymes show the ability to catalyze the same biochemical reaction. In contrast, as a result of functional diversification, many structurally similar enzyme molecules act on substantially distinct substrates and catalyze diverse biochemical reactions. Here, we present updates on the ATP-grasp, alkaline phosphatase, cupin, HD hydrolase, and N-terminal nucleophile (Ntn) hydrolase enzyme superfamilies and discuss the patterns of sequence and structural conservation and diversity within these superfamilies. Typically, enzymes within a superfamily possess common sequence motifs and key active site residues, as well as (predicted) reaction mechanisms. These observations suggest that the strained conformation (the entatic state) of the active site, which is responsible for the substrate binding and formation of the transition complex, tends to be conserved within enzyme superfamilies. The subsequent fate of the transition complex is not necessarily conserved and depends on the details of the structures of the enzyme and the substrate. This variability of reaction outcomes limits the ability of sequence analysis to predict the exact enzymatic activities of newly sequenced gene products. Nevertheless, sequence-based (super)family assignments and generic functional predictions, even if imprecise, provide valuable leads for experimental studies and remain the best approach to the functional annotation of uncharacterized proteins from new genomes.

Computational insights into the oxidation of mono- and 1,4 disubstituted arenes by the Toluene Dioxygenase enzymatic complex

An automated, active site-focused, computational method is described for use in predicting the effects of engineered amino acid mutations on enzyme catalytic activity. The method uses structure-based function descriptors (Fuzzy Functional Forms trade mark or FFFs trade mark ) to automatically identify enzyme functional sites in proteins. Three-dimensional sequence profiles are created from the surrounding active site structure. The computationally derived active site profile is used to analyze the effect of each amino acid change by defining three key features: proximity of the change to the active site, degree of amino acid conservation at the position in related proteins, and compatibility of the change with residues observed at that position in similar proteins. The features were analyzed using a data set of individual amino acid mutations occurring at 128 residue positions in 14 different enzymes. The results show that changes at key active site residues and at highly conserved positions are likely to have deleterious effects on the catalytic activity, and that non-conservative mutations at highly conserved residues are even more likely to be deleterious. Interestingly, the study revealed that amino acid substitutions at residues in close contact with the key active site residues are not more likely to have deleterious effects than mutations more distant from the active site. Utilization of the FFF-derived structural information yields a prediction method that is accurate in 79-83% of the test cases. The success of this method across all six EC classes suggests that it can be used generally to predict the effects of mutations and nsSNPs for enzymes. Future applications of the approach include automated, large-scale identification of deleterious nsSNPs in clinical populations and in large sets of disease-associated nsSNPs, and identification of deleterious nsSNPs in drug targets and drug metabolizing enzymes.

Herein we report the development of a new periodate-based reactive assay system for the fluorescent detection of the cis-diol metabolites produced by Rieske dioxygenases. This sensitive and diastereoselective assay system successfully evaluates the substrate scope of Rieske dioxygenases and determines the relative activity of a rationally designed Rieske dioxygenase variant library. The high throughput capacity of the assay system enables rapid and efficient substrate scope investigations and screening of large dioxygenase variant libraries.

L-arginine deiminase (ADI) catalyzes the irreversible hydrolysis of L-arginine to citrulline and ammonia. In a previous report of the structure of apoADI from Pseudomonas aeruginosa, the four residues that form the catalytic motif were identified as Cys406, His278, Asp280, and Asp166. The function of Cys406 in nucleophilic catalysis has been demonstrated by transient kinetic studies. In this study, the structure of the C406A mutant in complex with L-arginine is reported to provide a snapshot of the enzyme.substrate complex. Through the comparison of the structures of apoenzyme and substrate-bound enzyme, a substrate-induced conformational transition, which might play an important role in activity regulation, was discovered. Furthermore, the position of the guanidinium group of the bound substrate relative to the side chains of His278, Asp280, and Asp166 indicated that these residues mediate multiple proton transfers. His278 and Asp280, which are positioned to activate the water nucleophile in the hydrolysis of the S-alkylthiouronium intermediate, were replaced with alanine to stabilize the intermediate for structure determination. The structures determined for the H278A and D280A mutants co-crystallized with L-arginine provide a snapshot of the S-alkylthiouronium adduct formed by attack of Cys406 on the guanidinium carbon of L-arginine followed by the elimination of ammonia. Asp280 and Asp166 engage in ionic interactions with the guanidinium group in the C406A ADI. L-arginine structure and might orient the reaction center and participate in proton transfer. Structure determination of D166A revealed the apoD166A ADI. The collection of structures is interpreted in the context of recent biochemical data to propose a model for ADI substrate recognition and catalysis.

A Conserved Negatively Charged Cluster in the Active Site of Creatine Kinase Is Critical for Enzymatic Activity

Free energy analysis of ω-transaminase reactions to dissect how the enzyme controls the substrate selectivity

Characterization of the Substrate Specificity of PhlD, a Type III Polyketide Synthase from Pseudomonas fluorescens

This work describes a site-directed mutagenesis study of pentaerythritol tetranitrate reductase (PETN reductase) to probe the role of key active site residues in influencing both product enantiopurity and the ratio of CC vs. nitro-group reduction with 2-phenyl-1-nitropropene. Comparative biotransformations of wild type and single/double mutants of PETN reductase with 2-phenyl-1-nitropropene showed that one enzyme scaffold was capable of generating both enantiomeric products with improved enantiopurities by a manipulation of the reaction conditions and/or the presence of a one or two key mutations. These changes located at key active site residues were sufficient to moderately improve product enantiopurity, cause a switch in the major product enantiomer formed and/or promote or eliminate side-product formation. The mutation of substrate-binding residue Y351 to alanine and phenylalanine improved the biocatalytic potential of PETN reductase by the elimination of a competing side reaction. The crystal structures of three mutants at residue Y351 (PDB codes: 3P81, 3P84 and 3P8J) show that only subtle changes in the active site environment may be necessary to generate significantly improved biocatalysts.

Rieske dioxygenases catalyze the initial steps in the hydroxylation of aromatic compounds and are critical for the metabolism of xenobiotic substances. Because substrates do not bind to the mononuclear non-heme FeII center, elementary steps leading to O2 activation and substrate hydroxylation are difficult to delineate, thus making it challenging to rationalize divergent observations on enzyme mechanisms, reactivity, and substrate specificity. Here, we show for nitrobenzene dioxygenase, a Rieske dioxygenase capable of transforming nitroarenes to nitrite and substituted catechols, that unproductive O2 activation with the release of the unreacted substrate and reactive oxygen species represents an important path in the catalytic cycle. Through correlation of O2 uncoupling for a series of substituted nitroaromatic compounds with 18O and 13C kinetic isotope effects of dissolved O2 and aromatic substrates, respectively, we show that O2 uncoupling occurs after the rate-limiting formation of FeIII-(hydro)peroxo species from which substrates are hydroxylated. Substituent effects on the extent of O2 uncoupling suggest that the positioning of the substrate in the active site rather than the susceptibility of the substrate for attack by electrophilic oxygen species is responsible for unproductive O2 uncoupling. The proposed catalytic cycle provides a mechanistic basis for assessing the very different efficiencies of substrate hydroxylation vs unproductive O2 activation and generation of reactive oxygen species in reactions catalyzed by Rieske dioxygenases.

Density functional theory with the B3LYP hybrid functional has been used to study the mechanisms for dioxygen activation by four families of mononuclear non-heme iron enzymes: alpha-ketoacid-dependent dioxygenases, tetrahydrobiopterin-dependent hydroxylases, extradiol dioxygenases, and Rieske dioxygenases. These enzymes have a common active site with a ferrous ion coordinated to two histidines and one carboxylate group (aspartate or glutamate). In contrast to the heme case, this type of weak field environment always leads to a high-spin ground state. With the exception of the Rieske dioxygenases, which have an electron source outside the active site, the dioxygen activation process passes through the formation of a bridging-peroxide species, which then undergoes O-O bond cleavage finally leading to the four electron reduction of O(2). In the case of tetrahydrobiopterin- and alpha-ketoacid-dependent enzymes, the O-O heterolysis yields a high-valent iron-oxo species, which is capable of performing a two-electron oxidation chemistry on various organic substrates. For the other two families of enzymes (extradiol dioxygenases and Rieske dioxygenases) the substrate oxidation and the O-O bond cleavage are found to be coupled. In the extradiol dioxygenases the product of the O-O bond cleavage is a ferric iron with an oxy-substrate with a mixture of radical and anionic character, which is essential for the selectivity of the catechol cleavage.

Rieske oxygenases (ROs) catalyze a large range of oxidative chemistry. We have shown that cis-dihydrodiol-forming Rieske dioxygenases first react with their aromatic substrates via an active site nonheme Fe(III)-superoxide; electron transfer from the Rieske cluster then completes the product-forming reaction. Alternatively, two-electron-reduced Fe(III)-peroxo or hydroxo-Fe(V)-oxo activated oxygen intermediates are possible and may be utilized by other ROs to expand the catalytic range. Here, the reaction of a Rieske monooxygenase, salicylate 5-hydroxylase, that does not form a cis-dihydrodiol is examined. Single-turnover kinetic studies show fast binding of salicylate and O2. Transfer of the Rieske electron required to form the gentisate product occurs through bonds over ∼12 Å and must also be very fast. However, the observed rate constant for this reaction is much slower than expected and sensitive to substrate type. This suggests that initial reaction with salicylate occurs using the same Fe(III)-superoxo-level intermediate as Rieske dioxygenases and that this reaction limits the observed rate of electron transfer. A transient intermediate (λmax = 700 nm) with an electron paramagnetic resonance (EPR) at g = 4.3 is observed after the product is formed in the active site. The use of 17O2 (I = 5/2) results in hyperfine broadening of the g = 4.3 signal, showing that gentisate binds to the mononuclear iron via its C5-OH in the intermediate. The chromophore and EPR signal allow study of product release in the catalytic cycle. Comparison of the kinetics of single- and multiple-turnover reactions shows that re-reduction of the metal centers accelerates product release ∼300-fold, providing insight into the regulatory mechanism of ROs.