The influence of key residues in the tunnel entrance and the active site on activity and selectivity of toluene-4-monooxygenase

Whole-cell biocatalysis for epoxidation using the non-heme diiron monooxygenase system PmlABCDEF

Non-heme diiron-dependent monooxygenases are versatile oxidative enzymes that catalyze a wide range of oxy-functionalization reactions, including the conversion of CC bonds into epoxides. Despite their potential, their use as biocatalysts is limited due to challenges in expressing their multi-subunit structure in conventional hosts. Consequently, these enzymes exhibit relatively low activity in epoxidation reactions, significantly restricting their practical applications. Nevertheless, pairing appropriate enzymatic systems with carefully chosen expression vectors in suitable host organisms shows great promise. This chapter focuses on the PmlABCDEF monooxygenase system, expressed in Pseudomonas putida KT2440 and utilized as a whole-cell biocatalyst for epoxidation. The methodology details the construction of two distinct expression plasmids and their application in the preparative-scale bio-epoxidation reactions using two different substrate types. This approach is versatile and can be readily adapted for a broad range of epoxide production and other oxygenation-based transformations.

Targeted Mutation of a Non-catalytic Gating Residue Increases the Rate of Pseudomonas aeruginosa d-Arginine Dehydrogenase Catalytic Turnover

Commercial food and l-amino acid industries rely on bioengineered d-amino acid oxidizing enzymes to detect and remove d-amino acid contaminants. However, the bioengineering of enzymes to generate faster biological catalysts has proven difficult as a result of the failure to target specific kinetic steps that limit enzyme turnover, kcat, and the poor understanding of loop dynamics critical for catalysis. Pseudomonas aeruginosa d-arginine dehydrogenase (PaDADH) oxidizes most d-amino acids and is a good candidate for application in the l-amino acid and food industries. The side chain of the loop L2 E246 residue located at the entrance of the PaDADH active site pocket potentially favors the closed active site conformation and secures the substrate upon binding. This study used site-directed mutagenesis, steady-state, and rapid reaction kinetics to generate the glutamine, glycine, and leucine variants and investigate whether increasing the rate of product release could translate to an increased enzyme turnover rate. Upon E246 mutation to glycine, there was an increased rate of d-arginine turnover kcat from 122 to 500 s-1. Likewise, the kcat values increased 2-fold for the glutamine or leucine variants. Thus, we have engineered a faster biocatalyst for industrial applications by selectively increasing the rate of the PaDADH product release.

Heterogeneity in M. tuberculosis β-lactamase inhibition by Sulbactam

For decades, researchers have elucidated essential enzymatic functions on the atomic length scale by tracing atomic positions in real-time. Our work builds on possibilities unleashed by mix-and-inject serial crystallography (MISC) at X-ray free electron laser facilities. In this approach, enzymatic reactions are triggered by mixing substrate or ligand solutions with enzyme microcrystals. Here, we report in atomic detail (between 2.2 and 2.7 Å resolution) by room-temperature, time-resolved crystallography with millisecond time-resolution (with timepoints between 3 ms and 700 ms) how the Mycobacterium tuberculosis enzyme BlaC is inhibited by sulbactam (SUB). Our results reveal ligand binding heterogeneity, ligand gating, cooperativity, induced fit, and conformational selection all from the same set of MISC data, detailing how SUB approaches the catalytic clefts and binds to the enzyme noncovalently before reacting to a trans-enamine. This was made possible in part by the application of singular value decomposition to the MISC data using a program that remains functional even if unit cell parameters change up to 3 Å during the reaction.

Mutation of a distal gating residue modulates NADH binding in NADH:Quinone oxidoreductase from Pseudomonas aeruginosa PAO1

Enzymes require flexible regions to adopt multiple conformations during catalysis. The mobile regions of enzymes include gates that modulate the passage of molecules in and out of the enzyme's active site. The enzyme PA1024 from Pseudomonas aeruginosa PA01 is a recently discovered flavin-dependent NADH:quinone oxidoreductase (NQO, EC 1.6.5.9). Q80 in loop 3 (residues 75-86) of NQO is ∼15 Å away from the flavin and creates a gate that seals the active site through a hydrogen bond with Y261 upon NADH binding. In this study, we mutated Q80 to glycine, leucine, or glutamate to investigate the mechanistic significance of distal residue Q80 in NADH binding in the active site of NQO. The UV-visible absorption spectrum reveals that the mutation of Q80 minimally affects the protein microenvironment surrounding the flavin. The anaerobic reductive half-reaction of the NQO-mutants yields a ≥25-fold increase in the Kd value for NADH compared to the WT enzyme. However, we determined that the kred value was similar in the Q80G, Q80L, and wildtype enzymes and only ∼25% smaller in the Q80E enzyme. Steady-state kinetics with NQO-mutants and NQO-WT at varying concentrations of NADH and 1,4-benzoquinone establish a ≤5-fold decrease in the kcat/KNADH value. Moreover, there is no significant difference in the kcat/KBQ (∼1 × 106 M-1s-1) and kcat (∼24 s-1) values in NQO-mutants and NQO-WT. These results are consistent with the distal residue Q80 being mechanistically essential for NADH binding to NQO with minimal effect on the quinone binding to the enzyme and hydride transfer from NADH to flavin.

Heterogeneity in the M. tuberculosis β-Lactamase Inhibition by Sulbactam

For decades, researchers have been determined to elucidate essential enzymatic functions on the atomic lengths scale by tracing atomic positions in real time. Our work builds on new possibilities unleashed by mix-and-inject serial crystallography (MISC) 1-5 at X-ray free electron laser facilities. In this approach, enzymatic reactions are triggered by mixing substrate or ligand solutions with enzyme microcrystals 6 . Here, we report in atomic detail and with millisecond time-resolution how the Mycobacterium tuberculosis enzyme BlaC is inhibited by sulbactam (SUB). Our results reveal ligand binding heterogeneity, ligand gating 7-9 , cooperativity, induced fit 10,11 and conformational selection 11-13 all from the same set of MISC data, detailing how SUB approaches the catalytic clefts and binds to the enzyme non-covalently before reacting to a trans- enamine. This was made possible in part by the application of the singular value decomposition 14 to the MISC data using a newly developed program that remains functional even if unit cell parameters change during the reaction.

Structural basis for peptide substrate specificities of glycosyltransferase GalNAc-T2

The polypeptide N-acetylgalactosaminyl transferase (GalNAc-T) enzyme family initiates O-linked mucin-type glycosylation. The family constitutes 20 isoenzymes in humans. GalNAc-Ts exhibit both redundancy and finely tuned specificity for a wide range of peptide substrates. In this work, we deciphered the sequence and structural motifs that determine the peptide substrate preferences for the GalNAc-T2 isoform. Our approach involved sampling and characterization of peptide-enzyme conformations obtained from Rosetta Monte Carlo-minimization-based flexible docking. We computationally scanned 19 amino acid residues at positions -1 and +1 of an eight-residue peptide substrate, which comprised a dataset of 361 (19x19) peptides with previously characterized experimental GalNAc-T2 glycosylation efficiencies. The calculations recapitulated experimental specificity data, successfully discriminating between glycosylatable and non-glycosylatable peptides with a probability of 96.5% (ROC-AUC score), a balanced accuracy of 85.5% and a false positive rate of 7.3%. The glycosylatable peptide substrates viz. peptides with proline, serine, threonine, and alanine at the -1 position of the peptide preferentially exhibited cognate sequon-like conformations. The preference for specific residues at the -1 position of the peptide was regulated by enzyme residues R362, K363, Q364, H365 and W331, which modulate the pocket size and specific enzyme-peptide interactions. For the +1 position of the peptide, enzyme residues K281 and K363 formed gating interactions with aromatics and glutamines at the +1 position of the peptide, leading to modes of peptide-binding sub-optimal for catalysis. Overall, our work revealed enzyme features that lead to the finely tuned specificity observed for a broad range of peptide substrates for the GalNAc-T2 enzyme. We anticipate that the key sequence and structural motifs can be extended to analyze specificities of other isoforms of the GalNAc-T family and can be used to guide design of variants with tailored specificity.

BTEX biodegradation by Bacillus amyloliquefaciens subsp. plantarum W1 and its proposed BTEX biodegradation pathways

Benzene, toluene, ethylbenzene and (p-, m- and o-) xylene (BTEX) are classified as main pollutants by several environmental protection agencies. In this study, a non-pathogenic, Gram-positive rod-shape bacterium with an ability to degrade all six BTEX compounds, employed as an individual substrate or as a mixture, was isolated. The bacterial isolate was identified as Bacillus amyloliquefaciens subsp. plantarum strain W1. An overall BTEX biodegradation (as individual substrates) by strain W1 could be ranked as: toluene > benzene, ethylbenzene, p-xylene > m-xylene > o-xylene. When presented in a BTEX mixture, m-xylene and o-xylene biodegradation was slightly improved suggesting an induction effect by other BTEX components. BTEX biodegradation pathways of strain W1 were proposed based on analyses of its metabolic intermediates identified by LC–MS/MS. Detected activity of several putative monooxygenases and dioxygenases suggested the versatility of strain W1. Thus far, this is the first report of biodegradation pathways for all of the six BTEX compounds by a unique bacterium of the genus Bacillus. Moreover, B. amyloliquefaciens subsp. plantarum W1 could be a good candidate for an in situ bioremediation considering its Generally Recognized as Safe (GRAS) status and a possibility to serve as a plant growth-promoting rhizobacterium (PGPR).

Filling the Void: Introducing Aromatic Interactions into Solvent Tunnels To Enhance Lipase Stability in Methanol

An enhanced stability of enzymes in organic solvents is desirable under industrial conditions. The potential of lipases as biocatalysts is mainly limited by their denaturation in polar alcohols. In this study, we focused on selected solvent tunnels in lipase from Geobacillus stearothermophilus T6 to improve its stability in methanol during biodiesel synthesis. Using rational mutagenesis, bulky aromatic residues were incorporated to occupy solvent channels and induce aromatic interactions leading to a better inner core packing. The chemical and structural characteristics of each solvent tunnel were systematically analyzed. Selected residues were replaced with Phe, Tyr, or Trp. Overall, 16 mutants were generated and screened in 60% methanol, from which 3 variants showed an enhanced stability up to 81-fold compared with that of the wild type. All stabilizing mutations were found in the longest tunnel detected in the "closed-lid" X-ray structure. The combination of Phe substitutions in an A187F/L360F double mutant resulted in an increase in unfolding temperature (T_m ) of 7°C in methanol and a 3-fold increase in biodiesel synthesis yield from waste chicken oil. A kinetic analysis with p-nitrophenyl laurate revealed that all mutants displayed lower hydrolysis rates (k _cat), though their stability properties mostly determined the transesterification capability. Seven crystal structures of different variants were solved, disclosing new π-π or CH/π intramolecular interactions and emphasizing the significance of aromatic interactions for improved solvent stability. This rational approach could be implemented for the stabilization of other enzymes in organic solvents.IMPORTANCE Enzymatic synthesis in organic solvents holds increasing industrial opportunities in many fields; however, one major obstacle is the limited stability of biocatalysts in such a denaturing environment. Aromatic interactions play a major role in protein folding and stability, and we were inspired by this to redesign enzyme voids. The rational protein engineering of solvent tunnels of lipase from Geobacillus stearothermophilus is presented here, offering a promising approach to introduce new aromatic interactions within the enzyme core. We discovered that longer tunnels leading from the surface to the enzyme active site were more beneficial targets for mutagenesis for improving lipase stability in methanol during biodiesel biosynthesis. A structural analysis of the variants confirmed the generation of new interactions involving aromatic residues. This work provides insights into stability-driven enzyme design by targeting the solvent channel void.

Mutagenesis and expression of methane monooxygenase to alter regioselectivity with aromatic substrates

Soluble methane monooxygenase (sMMO) from methane-oxidising bacteria can oxygenate more than 100 hydrocarbons and is one of the most catalytically versatile biological oxidation catalysts. Expression of recombinant sMMO has to date not been achieved in Escherichia coli and so an alternative expression system must be used to manipulate it genetically. Here we report substantial improvements to the previously described system for mutagenesis of sMMO and expression of recombinant enzymes in a methanotroph (Methylosinus trichosporium OB3b) expression system. This system has been utilised to make a number of new mutants and to engineer sMMO to increase its catalytic precision with a specific substrate whilst increasing activity by up to 6-fold. These results are the first 'proof-of-principle' experiments illustrating the feasibility of developing sMMO-derived catalysts for diverse applications.

Cavity residue leucine 95 and channel residues glutamine 204, aspartic acid 211, and phenylalanine 269 of toluene o-xylene monooxygenase influence catalysis

Structural analysis of toluene-o-xylene monooxygenase (ToMO) hydroxylase revealed the presence of three hydrophobic cavities, a channel, and a pore leading from the protein surface to the active site. Here, saturation mutagenesis was used to investigate the catalytic roles of alpha-subunit (TouA) second cavity residue L95 and TouA channel residues Q204, D211, and F269. By testing the substrates toluene, phenol, nitrobenzene, and/or naphthalene, these positions were found to influence the catalytic activity of ToMO. Several regiospecific variants were identified from TouA positions Q204, F269, and L95. For example, TouA variant Q204H had the regiospecificity of nitrobenzene changed significantly from 30 to 61 % p-nitrophenol. Interestingly, a combination of mutations at Q204H and A106V altered the regiospecificity of nitrobenzene back to 27 % p-nitrophenol. TouA variants F269Y, F269P, Q204E, and L95D improved the meta-hydroxylating capability of nitrobenzene by producing 87, 85, 82, and 77 % m-nitrophenol, respectively. For naphthalene oxidation, TouA variants F269V, Q204A, Q204S/S222N, and F269T had the regiospecificity changed from 16 to 9, 10, 23, and 25 % 2-naphthol, respectively. Here, two additional TouA residues, S222 and A106, were also identified that may have important roles in catalysis. Most of the isolated variants from D211 remained active, whereas having a hydrophobic residue at this position appeared to diminish the catalytic activity toward naphthalene. The mutational effects on the ToMO regiospecificity described here suggest that it is possible to further fine tune and engineer the reactivity of multicomponent diiron monooxygenases toward different substrates at positions that are relatively distant from the active site.

Substrate tunnels in enzymes: structure-function relationships and computational methodology

In enzymes, the active site is the location where incoming substrates are chemically converted to products. In some enzymes, this site is deeply buried within the core of the protein, and, in order to access the active site, substrates must pass through the body of the protein via a tunnel. In many systems, these tunnels act as filters and have been found to influence both substrate specificity and catalytic mechanism. Identifying and understanding how these tunnels exert such control has been of growing interest over the past several years because of implications in fields such as protein engineering and drug design. This growing interest has spurred the development of several computational methods to identify and analyze tunnels and how ligands migrate through these tunnels. The goal of this review is to outline how tunnels influence substrate specificity and catalytic efficiency in enzymes with buried active sites and to provide a brief summary of the computational tools used to identify and evaluate these tunnels.

Ensemble generation and the influence of protein flexibility on geometric tunnel prediction in cytochrome P450 enzymes

Computational prediction of ligand entry and egress paths in proteins has become an emerging topic in computational biology and has proven useful in fields such as protein engineering and drug design. Geometric tunnel prediction programs, such as Caver3.0 and MolAxis, are computationally efficient methods to identify potential ligand entry and egress routes in proteins. Although many geometric tunnel programs are designed to accommodate a single input structure, the increasingly recognized importance of protein flexibility in tunnel formation and behavior has led to the more widespread use of protein ensembles in tunnel prediction. However, there has not yet been an attempt to directly investigate the influence of ensemble size and composition on geometric tunnel prediction. In this study, we compared tunnels found in a single crystal structure to ensembles of various sizes generated using different methods on both the apo and holo forms of cytochrome P450 enzymes CYP119, CYP2C9, and CYP3A4. Several protein structure clustering methods were tested in an attempt to generate smaller ensembles that were capable of reproducing the data from larger ensembles. Ultimately, we found that by including members from both the apo and holo data sets, we could produce ensembles containing less than 15 members that were comparable to apo or holo ensembles containing over 100 members. Furthermore, we found that, in the absence of either apo or holo crystal structure data, pseudo-apo or -holo ensembles (e.g. adding ligand to apo protein throughout MD simulations) could be used to resemble the structural ensembles of the corresponding apo and holo ensembles, respectively. Our findings not only further highlight the importance of including protein flexibility in geometric tunnel prediction, but also suggest that smaller ensembles can be as capable as larger ensembles at capturing many of the protein motions important for tunnel prediction at a lower computational cost.

Atomic picture of ligand migration in toluene 4-monooxygenase

Computational modeling combined with mutational and activity assays was used to underline the substrate migration pathways in toluene 4-monooxygenase, a member of the important family of bacterial multicomponent monooxygenases (BMMs). In all structurally defined BMM hydroxylases, several hydrophobic cavities in the α-subunit map a preserved path from the protein surface to the diiron active site. Our results confirm the presence of two pathways by which different aromatic molecules can enter/escape the active site. While the substrate is observed to enter from both channels, the more hydrophilic product is withdrawn mainly from the shorter channel ending at residues D285 and E214. The long channel ends in the vicinity of S395, whose variants have been seen to affect activity and specificity. These mutational effects are clearly reproduced and rationalized by the in silico studies. Furthermore, the combined computational and experimental results highlight the importance of residue F269, which is located at the intersection of the two channels.

The key role of a non-active-site residue Met148 on the catalytic efficiency of meta-cleavage product hydrolase BphD

meta-Cleavage product (MCP) hydrolases (EC 3.7.1.9) can catalyze a specific C-C bond fission during the microbial aerobic degradation of aromatics. The previous studies on structure-function relationship of MCP hydrolases mainly focus on the active site residues by site-directed mutagenesis. However, the information about the role of the non-active-site residues is still unclear. In this study, a non-active-site residue Met148 of MCP hydrolase BphD was selected as the mutagenesis site according to the sequence alignments, structure superimpose and the tunnel analysis, which underwent the saturation mutagenesis resulting 19 mutants. The catalytic efficiencies of the mutants on 6-oxo-6-phenylhexa-2,4-dienoic acid (HOPDA) were all decreased compared with the wild-type one except for the M148D mutant. Especially, the M148P mutant exhibited 290-fold lower k cat/K m than that of the wild-type BphD. Transient kinetic analyses of M148P showed the reciprocal relaxation time corresponded to C-C bond cleavage and product release steps (9.6 s(-1)) was 4.08-fold lower than BphD WT (39.2 s(-1)). Tunnel cluster analysis of BphD WT, M148P and M148W demonstrated that only the bulky Trp148 could block tunnel T2 in the BphD WT, but it exhibited slight effects on the catalytic efficiency (0.94-fold of BphD WT). Therefore, product release was not the main reason for the efficiency decrease of M148P. On the other hand, molecular dynamics simulations on the BphD WT and BphD M148P in complex with HOPDA indicated that the dramatic decrease of the catalytic efficiencies of BphD M148P should be due to the unproductive binding of HOPDA. The study demonstrated the catalytic efficiency of MCP hydrolase can be engineered by modification of non-active site residue.

Engineering non-heme mono- and dioxygenases for biocatalysis

Oxygenases are ubiquitous enzymes that catalyze the introduction of one or two oxygen atoms to unreactive chemical compounds. They require reduction equivalents from NADH or NADPH and comprise metal ions, metal ion complexes, or coenzymes in their active site. Thus, for industrial purposes, oxygenases are most commonly employed using whole cell catalysis, to alleviate the need for co-factor regeneration. Biotechnological applications include bioremediation, chiral synthesis, biosensors, fine chemicals, biofuels, pharmaceuticals, food ingredients and polymers. Controlling activity and selectivity of oxygenases is therefore of great importance and of growing interest to the scientific community. This review focuses on protein engineering of non-heme monooxygenases and dioxygenases for generating improved or novel functionalities. Rational mutagenesis based on x-ray structures and sequence alignment, as well as random methods such as directed evolution, have been utilized. It is concluded that knowledge-based protein engineering accompanied with targeted libraries, is most efficient for the design and tuning of biocatalysts towards novel substrates and enhanced catalytic activity while minimizing the screening efforts.

Improving biocatalyst performance by integrating statistical methods into protein engineering

Directed evolution and rational design were used to generate active variants of toluene-4-monooxygenase (T4MO) on 2-phenylethanol (PEA), with the aim of producing hydroxytyrosol, a potent antioxidant. Due to the complexity of the enzymatic system-four proteins encoded by six genes-mutagenesis is labor-intensive and time-consuming. Therefore, the statistical model of Nov and Wein (J. Comput. Biol. 12:247-282) was used to reduce the number of variants produced and evaluated in a lab. From an initial data set of 24 variants, with mutations at nine positions, seven double or triple mutants were identified through statistical analysis. The average activity of these mutants was 4.6-fold higher than the average activity of the initial data set. In an attempt to further improve the enzyme activity to obtain PEA hydroxylation, a second round of statistical analysis was performed. Nine variants were considered, with 3, 4, and 5 point mutations. The average activity of the variants obtained in the second statistical round was 1.6-fold higher than in the first round and 7.3-fold higher than that of the initial data set. The best variant discovered, TmoA I100A E214G D285Q, exhibited an initial oxidation rate of 4.4 ± 0.3 nmol/min/mg protein, which is 190-fold higher than the rate obtained by the wild type. This rate was also 2.6-fold higher than the activity of the wild type on the natural substrate toluene. By considering only 16 preselected mutants (out of ∼13,000 possible combinations), a highly active variant was discovered with minimum time and effort.