Protein engineering of toluene monooxygenases for synthesis of chiral sulfoxides

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.

RNA-seq transcriptome analysis of a Pseudomonas strain with diversified catalytic properties growth under different culture medium

Biocatalysis is an emerging strategy for the production of enantio-pure organic molecules. However, lacking of commercially available enzymes restricts the widespread application of biocatalysis. In this study, we report a Pseudomonas strain which exhibited versatile oxidation activity to synthesize chiral sulfoxides when growing under M9-toluene medium and reduction activity to synthesize chiral alcohols when on Luria-Bertani (LB) medium, respectively. Further comparative transcriptome analysis on samples from these two cultural conditions has identified 1038 differentially expressed genes (DEG). Gene Ontology (GO) enrichment and KEGG pathways analysis demonstrate significant changes in protein synthesis, energy metabolism, and biosynthesis of metabolites when cells cultured under different conditions. We have identified eight candidate enzymes from this bacterial which may have the potential to be used for synthesis of chiral alcohol and sulfoxide chemicals. This work provides insights into the mechanism of diversity in catalytic properties of this Pseudomonas strain growth with different cultural conditions, as well as candidate enzymes for further biocatalysis of enantiomerically pure molecules and pharmaceuticals.

Production of chiral compounds using immobilized cells as a source of biocatalysts

The importance of chiral compounds in all fields of technology and life sciences is shown. Small chiral molecules are mainly used as building blocks in the synthesis of more complex and functionalized compounds. Nature creates and imposes stereoselectivity by means of enzymes, which are highly efficient biocatalysts. The use of whole cells as a biocatalyst source is a promising strategy for avoiding some drawbacks associated with the use of pure enzymes, especially their high cost. The use of free cells is also challenging, since cell lysis can also occur under the reaction conditions. However, cell immobilization has been employed to increase the catalytic potential of enzymes by extending their lifetimes in organic solvents and non-natural environments. Besides, immobilized cells maintain their biocatalytic performance for several reaction cycles. Considering the above-mentioned arguments, several authors have synthesized different classes of chiral compounds such as alcohols, amines, carboxylic acids, amides, sulfides and lactones by means of immobilized cells. Our aim was to discuss the main aspects of the production of chiral compounds using immobilized cells as a source of biocatalysts, except under fermentation conditions.

The Toluene o-Xylene Monooxygenase Enzymatic Activity for the Biosynthesis of Aromatic Antioxidants

Monocyclic phenols and catechols are important antioxidant compounds for the food and pharmaceutic industries; their production through biotransformation of low-added value starting compounds is of major biotechnological interest. The toluene o-xylene monooxygenase (ToMO) from Pseudomonas sp. OX1 is a bacterial multicomponent monooxygenase (BMM) that is able to hydroxylate a wide array of aromatic compounds and has already proven to be a versatile biochemical tool to produce mono- and dihydroxylated derivatives of aromatic compounds. The molecular determinants of its regioselectivity and substrate specificity have been thoroughly investigated, and a computational strategy has been developed which allows designing mutants able to hydroxylate non-natural substrates of this enzyme to obtain high-added value compounds of commercial interest. In this work, we have investigated the use of recombinant ToMO, expressed in cells of Escherichia coli strain JM109, for the biotransformation of non-natural substrates of this enzyme such as 2-phenoxyethanol, phthalan and 2-indanol to produce six hydroxylated derivatives. The hydroxylated products obtained were identified, isolated and their antioxidant potential was assessed both in vitro, using the DPPH assay, and on the rat cardiomyoblast cell line H9c2. Incubation of H9c2 cells with the hydroxylated compounds obtained from ToMO-catalyzed biotransformation induced a differential protective effect towards a mild oxidative stress induced by the presence of sodium arsenite. The results obtained confirm once again the versatility of the ToMO system for oxyfunctionalization reactions of biotechnological importance. Moreover, the hydroxylated derivatives obtained possess an interesting antioxidant potential that encourages the use of the enzyme for further functionalization reactions and their possible use as scaffolds to design novel bioactive molecules.

Biotechnological production of chiral organic sulfoxides: current state and perspectives

Chiral organic sulfoxides (COSs) are important compounds that act as chiral auxiliaries in numerous asymmetric reactions and as intermediates in chiral drug synthesis. In addition to their optical resolution, stereoselective oxidation of sulfides can be used for COS production. This reaction is facilitated by oxygenases and peroxidases from various microbial resources. To meet the current demand for esomeprazole, a proton pump inhibitor used in the treatment of gastric-acid-related disorders, and the (S)-isomer of an organic sulfoxide compound, omeprazole, a successful biotechnological production method using a Baeyer-Villiger monooxygenase (BVMO), was developed. In this review, we summarize the recent advancements in COS production using biocatalysts, including enzyme identification, protein engineering, and process development.

Enantioselective imidation of sulfides via enzyme-catalyzed intermolecular nitrogen-atom transfer

Engineering enzymes with novel reaction modes promises to expand the applications of biocatalysis in chemical synthesis and will enhance our understanding of how enzymes acquire new functions. The insertion of nitrogen-containing functional groups into unactivated C–H bonds is not catalyzed by known enzymes but was recently demonstrated using engineered variants of cytochrome P450_BM3 (CYP102A1) from Bacillus megaterium. Here, we extend this novel P450-catalyzed reaction to include intermolecular insertion of nitrogen into thioethers to form sulfimides. An examination of the reactivity of different P450_BM3 variants toward a range of substrates demonstrates that electronic properties of the substrates are important in this novel enzyme-catalyzed reaction. Moreover, amino acid substitutions have a large effect on the rate and stereoselectivity of sulfimidation, demonstrating that the protein plays a key role in determining reactivity and selectivity. These results provide a stepping stone for engineering more complex nitrogen-atom-transfer reactions in P450 enzymes and developing a more comprehensive biocatalytic repertoire.

Fitness loss and library size determination in saturation mutagenesis

Saturation mutagenesis is a widely used directed evolution technique, in which a large number of protein variants, each having random amino acids in certain predetermined positions, are screened in order to discover high-fitness variants among them. Several metrics for determining the library size (the number of variants screened) have been suggested in the literature, but none of them incorporates the actual fitness of the variants discovered in the experiment. We present the results of an extensive simulation study, which is based on probabilistic models for protein fitness landscape, and which investigates how the result of a saturation mutagenesis experiment – the fitness of the best variant discovered – varies as a function of the library size. In particular, we study the loss of fitness in the experiment: the difference between the fitness of the best variant discovered, and the fitness of the best variant in variant space. Our results are that the existing criteria for determining the library size are conservative, so smaller libraries are often satisfactory. Reducing the library size can save labor, time, and expenses in the laboratory.

Protein engineering of nirobenzene dioxygenase for enantioselective synthesis of chiral sulfoxides

Nitrobenzene dioxygenase (NBDO) from Comamonas sp. is shown here to perform enantioselective oxidation of aromatic sulfides. Several para-substituted alkyl aryl sulfides were examined and it was found that the activity of the enzyme is dependent on the size of the substrate. Saturation mutagenesis was performed on different residues in the active site in order to improve activity and selectivity. Mutagenesis at position 258 in the α-hydroxylase subunit of NBDO improved both activity and enantioselectivity. Substitutions in position 293 improved the activity on all substrates and had diverse influence on enantioselectivity. Mutagenesis in position 207 provided two interesting variants, V207I and V207A, with opposite enantioselectivities. Furthermore, combining two favorable mutations, N258A and F293H, provided an improved variant with both higher activity (5.20 ± 0.01, 2.12 ± 0.21, 2.64 ± 0.14 and 4.01 ± 0.34 nmol min(-1) mg protein(-1) on thioanisole, ptolyl, Cl-thioanisole and Br-thioanisole, respectively, which is 1.7, 4.6, 7.1 and 26.7-fold compared with wild type) and improved enantioselectivity (e.g. 67% enantiomeric excess for Cl-thioanisole vs. 5% for wild type). Molecular docking and active site volume calculations were used to correlate between the structure of the substrates and the function of the enzymes. The results from this work suggest that the location of pro-chiral sulfides in the active site is coordinated by hydrophobic interactions and by steric considerations, which in turn influences the activity and enantioselectivity of NBDO.

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.

Crystallographic analysis of active site contributions to regiospecificity in the diiron enzyme toluene 4-monooxygenase

Crystal structures of toluene 4-monooxygenase hydroxylase in complex with reaction products and effector protein reveal active site interactions leading to regiospecificity. Complexes with phenolic products yield an asymmetric μ-phenoxo-bridged diiron center and a shift of diiron ligand E231 into a hydrogen bonding position with conserved T201. In contrast, complexes with inhibitors p-NH(2)-benzoate and p-Br-benzoate showed a μ-1,1 coordination of carboxylate oxygen between the iron atoms and only a partial shift in the position of E231. Among active site residues, F176 trapped the aromatic ring of products against a surface of the active site cavity formed by G103, E104 and A107, while F196 positioned the aromatic ring against this surface via a π-stacking interaction. The proximity of G103 and F176 to the para substituent of the substrate aromatic ring and the structure of G103L T4moHD suggest how changes in regiospecificity arise from mutations at G103. Although effector protein binding produced significant shifts in the positions of residues along the outer portion of the active site (T201, N202, and Q228) and in some iron ligands (E231 and E197), surprisingly minor shifts (<1 Å) were produced in F176, F196, and other interior residues of the active site. Likewise, products bound to the diiron center in either the presence or absence of effector protein did not significantly shift the position of the interior residues, suggesting that positioning of the cognate substrates will not be strongly influenced by effector protein binding. Thus, changes in product distributions in the absence of the effector protein are proposed to arise from differences in rates of chemical steps of the reaction relative to motion of substrates within the active site channel of the uncomplexed, less efficient enzyme, while structural changes in diiron ligand geometry associated with cycling between diferrous and diferric states are discussed for their potential contribution to product release.

Study of the TmoS/TmoT two-component system: towards the functional characterization of the family of TodS/TodT like systems

The two‐component system TmoS/TmoT controls the expression of the toluene‐4‐monooxygenase pathway in Pseudomonas mendocina RK1 via modulation of P_tmoX activity. The TmoS/TmoT system belongs to the family of TodS/TodT like proteins. The sensor kinase TmoS is a 108 kDa protein composed of seven different domains. Using isothermal titration calorimetry we show that purified TmoS binds a wide range of aromatic compounds with high affinities. Tightest ligand binding was observed for toluene (K_D = 150 nM), which corresponds to the highest affinity measured between an effector and a sensor kinase. Other compounds with affinities in the nanomolar range include benzene, the 3 xylene isomers, styrene, nitrobenzene or p‐chlorotoluene. We demonstrate that only part of the ligands that bind to TmoS increase protein autophosphorylation in vitro and consequently pathway expression in vivo. These compounds are referred to as agonists. Other TmoS ligands, termed antagonists, failed to increase TmoS autophosphorylation, which resulted in their incapacity to stimulate gene expression in vivo. We also show that TmoS saturated with different agonists differs in their autokinase activities. The effector screening of gene expression showed that promoter activity of P_tmoX and P_todX (controlled by the TodS/TodT system) is mediated by the same set of 22 compounds. The common structural feature of these compounds is the presence of a single aromatic ring. Among these ligands, toluene was the most potent inducer of both promoter activities. Information on the TmoS/TmoT and TodS/TodT system combined with a sequence analysis of family members permits to identify distinct features that define this protein family.

When second best is good enough: another probabilistic look at saturation mutagenesis

We developed new criteria for determining the library size in a saturation mutagenesis experiment. When the number of all possible distinct variants is large, any of the top-performing variants (e.g., any of the top three) is likely to meet the design requirements, so the probability that the library contains at least one of them is a sensible criterion for determining the library size. By using a criterion of this type, one may significantly reduce the library size and thus save costs and labor while minimally compromising the quality of the best variant discovered. We present the probabilistic tools underlying these criteria and use them to compare the efficiencies of four randomization schemes: NNN, which uses all 64 codons; NNB, which uses 48 codons; NNK, which uses 32 codons; and MAX, which assigns equal probabilities to each of the 20 amino acids. MAX was found to be the most efficient randomization scheme and NNN the least efficient. TopLib, a computer program for carrying out the related calculations, is available through a user-friendly Web server.

Whole-cell oxidation of omeprazole sulfide to enantiopure esomeprazole with Lysinibacillus sp. B71

Production of enantiopure esomeprazole by biocatalysis is of great demand by pharmaceutical industry. A Gram-positive bacterium oxidizing omeprazole sulfide 1a (5-methoxy-2-[((4-methoxy-3,5-dimethylpyridin-2-yl)methyl)thio]-1H-benzoimidazole) to (S)-sulfoxide esomeprazole 2a (S)-5-methoxy-2-[(4-methoxy-3,5-dimethylpyridin-2-yl) methylsulfinyl]-3H-benzoimidazole was isolated from soil polluted with elemental sulfur. The strain exhibited the highest identity with the genus Lysinibacillus and catalyzed oxidation of 1a into enantiopure esomeprazole with conversion of 77% in a stirred bioreactor, fed-batch culture. No consecutive oxidation of (S)-sulfoxide to sulfone was observed during whole-cell catalysis. The unique characteristics of the catalyst provide a solid basis for further improvement and development of sustainable green bioprocess.

Tuning the specificity of the recombinant multicomponent toluene o-xylene monooxygenase from Pseudomonas sp. strain OX1 for the biosynthesis of tyrosol from 2-phenylethanol

Biocatalysis is today a standard technology for the industrial production of several chemicals, and the number of biotransformation processes running on a commercial scale is constantly increasing. Among biocatalysts, bacterial multicomponent monooxygenases (BMMs), a diverse group of nonheme diiron enzymes that activate dioxygen, are of primary interest due to their ability to catalyze a variety of complex oxidations, including reactions of mono- and dihydroxylation of phenolic compounds. In recent years, both directed evolution and rational design have been successfully used to identify the molecular determinants responsible for BMM regioselectivity and to improve their activity toward natural and nonnatural substrates. Toluene o-xylene monooxygenase (ToMO) is a BMM isolated from Pseudomonas sp. strain OX1 which hydroxylates a wide spectrum of aromatic compounds. In this work we investigate the use of recombinant ToMO for the biosynthesis in recombinant cells of Escherichia coli strain JM109 of 4-hydroxyphenylethanol (tyrosol), an antioxidant present in olive oil, from 2-phenylethanol, a cheap and commercially available substrate. We initially found that wild-type ToMO is unable to convert 2-phenylethanol to tyrosol. This was explained by using a computational model which analyzed the interactions between ToMO active-site residues and the substrate. We found that residue F176 is the major steric hindrance for the correct positioning of the reaction intermediate leading to tyrosol production into the active site of the enzyme. Several mutants were designed and prepared, and we found that the combination of different mutations at position F176 with mutation E103G allows ToMO to convert up to 50% of 2-phenylethanol into tyrosol in 2 h.

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.

Monooxygenases as biocatalysts: Classification, mechanistic aspects and biotechnological applications

Monooxygenases are enzymes that catalyze the insertion of a single oxygen atom from O(2) into an organic substrate. In order to carry out this type of reaction, these enzymes need to activate molecular oxygen to overcome its spin-forbidden reaction with the organic substrate. In most cases, monooxygenases utilize (in)organic cofactors to transfer electrons to molecular oxygen for its activation. Monooxygenases typically are highly chemo-, regio-, and/or enantioselective, making them attractive biocatalysts. In this review, an exclusive overview of known monooxygenases is presented, based on the type of cofactor that these enzymes require. This includes not only the cytochrome P450 and flavin-dependent monooxygenases, but also enzymes that utilize pterin, metal ions (copper or iron) or no cofactor at all. As most of these monooxygenases require nicotinamide coenzymes as electron donors, also an overview of current methods for coenzyme regeneration is given. This latter overview is of relevance for the biotechnological applications of these oxidative enzymes.

Redesigning dehalogenase access tunnels as a strategy for degrading an anthropogenic substrate

Engineering enzymes to degrade anthropogenic compounds efficiently is challenging. We obtained Rhodococcus rhodochrous haloalkane dehalogenase mutants with up to 32-fold higher activity than wild type toward the toxic, recalcitrant anthropogenic compound 1,2,3-trichloropropane (TCP) using a new strategy. We identified key residues in access tunnels connecting the buried active site with bulk solvent by rational design and randomized them by directed evolution. The most active mutant has large aromatic residues at two out of three randomized positions and two positions modified by site-directed mutagenesis. These changes apparently enhance activity with TCP by decreasing accessibility of the active site for water molecules, thereby promoting activated complex formation. Kinetic analyses confirmed that the mutations improved carbon-halogen bond cleavage and shifted the rate-limiting step to the release of products. Engineering access tunnels by combining computer-assisted protein design with directed evolution may be a valuable strategy for refining catalytic properties of enzymes with buried active sites.

Rapid methods for high-throughput detection of sulfoxides

Enantiopure sulfoxides are prevalent in drugs and are useful chiral auxiliaries in organic synthesis. The biocatalytic enantioselective oxidation of prochiral sulfides is a direct and economical approach for the synthesis of optically pure sulfoxides. The selection of suitable biocatalysts requires rapid and reliable high-throughput screening methods. Here we present four different methods for detecting sulfoxides produced via whole-cell biocatalysis, three of which were exploited for high-throughput screening. Fluorescence detection based on the acid activation of omeprazole was utilized for high-throughput screening of mutant libraries of toluene monooxygenases, but no active variants have been discovered yet. The second method is based on the reduction of sulfoxides to sulfides, with the coupled release and measurement of iodine. The availability of solvent-resistant microtiter plates enabled us to modify the method to a high-throughput format. The third method, selective inhibition of horse liver alcohol dehydrogenase, was used to rapidly screen highly active and/or enantioselective variants at position V106 of toluene ortho-monooxygenase in a saturation mutagenesis library, using methyl-p-tolyl sulfide as the substrate. A success rate of 89% (i.e., 11% false positives) was obtained, and two new mutants were selected. The fourth method is based on the colorimetric detection of adrenochrome, a back-titration procedure which measures the concentration of the periodate-sensitive sulfide. Due to low sensitivity during whole-cell screening, this method was found to be useful only for determining the presence or absence of sulfoxide in the reaction. The methods described in the present work are simple and inexpensive and do not require special equipment.

HotSpot Wizard: a web server for identification of hot spots in protein engineering

HotSpot Wizard is a web server for automatic identification of 'hot spots' for engineering of substrate specificity, activity or enantioselectivity of enzymes and for annotation of protein structures. The web server implements the protein engineering protocol, which targets evolutionarily variable amino acid positions located in the active site or lining the access tunnels. The 'hot spots' for mutagenesis are selected through the integration of structural, functional and evolutionary information obtained from: (i) the databases RCSB PDB, UniProt, PDBSWS, Catalytic Site Atlas and nr NCBI and (ii) the tools CASTp, CAVER, BLAST, CD-HIT, MUSCLE and Rate4Site. The protein structure and e-mail address are the only obligatory inputs for the calculation. In the output, HotSpot Wizard lists annotated residues ordered by estimated mutability. The results of the analysis are mapped on the enzyme structure and visualized in the web browser using Jmol. The HotSpot Wizard server should be useful for protein engineers interested in exploring the structure of their favourite protein and for the design of mutations in site-directed mutagenesis and focused directed evolution experiments. HotSpot Wizard is available at http://loschmidt.chemi.muni.cz/hotspotwizard/.