Protein engineering of toluene 4-monooxygenase of Pseudomonas mendocina KR1 for synthesizing 4-nitrocatechol from nitrobenzene

Regiospecific Oxidation of Chlorobenzene to 4-Chlororesorcinol, Chlorohydroquinone, 3-Chlorocatechol and 4-Chlorocatechol by Engineered Toluene o-Xylene Monooxygenases

Toluene o-xylene monooxygenase (ToMO) was found to oxidize chlorobenzene to form 2-chlorophenol (2-CP, 4%), 3-CP (12%), and 4-CP (84%) with a total product formation rate of 1.2 ± 0.17 nmol/min/mg protein. It was also discovered that ToMO forms 4-chlorocatechol (4-CC) from 3-CP and 4-CP with initial rates of 0.54 ± 0.10 and 0.40 ± 0.04 nmol/min/mg protein, respectively, and chlorohydroquinone (CHQ, 13%), 4-chlororesorcinol (4-CR, 3%), and 3-CC (84%) from 2-CP with an initial product formation rate of 1.1 ± 0.32 nmol/min/mg protein. To increase the oxidation rate and alter the oxidation regiospecificity of chloroaromatics, as well as to study the roles of active site residues L192 and A107 of the alpha hydroxylase fragment of ToMO (TouA), we used the saturation mutagenesis approach of protein engineering. Thirteen TouA variants were isolated, among which some of the best substitutions uncovered here have never been studied before. Specifically, TouA variant L192V was identified which had 1.8-, 1.4-, 2.4-, and 4.8-fold faster hydroxylation activity toward chlorobenzene, 2-CP, 3-CP, and 4-CP, respectively, compared to the native ToMO. The L192V variant also had the regiospecificity of chlorobenzene changed from 4% to 13% 2-CP and produced the novel product 3-CC (4%) from 3-CP. Most of the isolated variants were identified to change the regiospecificity of oxidation. For example, compared to the native ToMO, variants A107T, A107N, and A107M produced 6.3-, 7.0-, and 7.3-fold more 4-CR from 2-CP, respectively, and variants A107G and A107G/L192V produced 3-CC (33 and 39%, respectively) from 3-CP whereas native ToMO did not. IMPORTANCE Chlorobenzene is a commonly used toxic solvent and listed as a priority environmental pollutant by the US Environmental Protection Agency. Here, we report that Escherichia coli TG1 cells expressing toluene o-xylene monooxygenase (ToMO) can successfully oxidize chlorobenzene to form dihydroxy chloroaromatics, which are valuable industrial compounds. ToMO performs this at room temperature in water using only molecular oxygen and a cofactor supplied by the cells. Using protein engineering techniques, we also isolated ToMO variants with enhanced oxidation activity as well as fine-tuned regiospecificities which make direct microbial oxygenations even more attractive. The significance of this work lies in the ability to degrade environmental pollutants while at the same time producing valuable chemicals using environmentally benign biological methods rather than expensive, complex chemical processes.

Nuclear Resonance Vibrational Spectroscopic Definition of the Fe(IV)2 Intermediate Q in Methane Monooxygenase and Its Reactivity

Methanotrophic bacteria utilize the nonheme diiron enzyme soluble methane monooxygenase (sMMO) to convert methane to methanol in the first step of their metabolic cycle under copper-limiting conditions. The structure of the sMMO Fe(IV)₂ intermediate Q responsible for activating the inert C-H bond of methane (BDE = 104 kcal/mol) remains controversial, with recent studies suggesting both "open" and "closed" core geometries for its active site. In this study, we employ nuclear resonance vibrational spectroscopy (NRVS) to probe the geometric and electronic structure of intermediate Q at cryogenic temperatures. These data demonstrate that Q decays rapidly during the NRVS experiment. Combining data from several years of measurements, we derive the NRVS vibrational features of intermediate Q as well as its cryoreduced decay product. A library of 90 open and closed core models of intermediate Q is generated using density functional theory to analyze the NRVS data of Q and its cryoreduced product as well as prior spectroscopic data on Q. Our analysis reveals that a subset of closed core models reproduce these newly acquired NRVS data as well as prior data. The reaction coordinate with methane is also evaluated using both closed and open core models of Q. These studies show that the potent reactivity of Q toward methane resides in the "spectator oxo" of its Fe(IV)₂O₂ core, in contrast to nonheme mononuclear Fe(IV)═O enzyme intermediates that H atoms abstract from weaker C-H bonds.

Transcriptomic analysis of Burkholderia cenocepacia CEIB S5-2 during methyl parathion degradation

Methyl parathion (MP) is a highly toxic organophosphorus pesticide associated with water, soil, and air pollution events. The identification and characterization of microorganisms capable of biodegrading pollutants are an important environmental task for bioremediation of pesticide impacted sites. The strain Burkholderia cenocepacia CEIB S5-2 is a bacterium capable of efficiently hydrolyzing MP and biodegrade p-nitrophenol (PNP), the main MP hydrolysis product. Due to the high PNP toxicity over microbial living forms, the reports on bacterial PNP biodegradation are scarce. According to the genomic data, the MP- and PNP-degrading ability observed in B. cenocepacia CEIB S5-2 is related to the presence of the methyl parathion-degrading gene (mpd) and the gene cluster pnpABA'E1E2FDC, which include the genes implicated in the PNP degradation. In this work, the transcriptomic analysis of the strain in the presence of MP revealed the differential expression of 257 genes, including all genes implicated in the PNP degradation, as well as a set of genes related to the sensing of environmental changes, the response to stress, and the degradation of aromatic compounds, such as translational regulators, membrane transporters, efflux pumps, and oxidative stress response genes. These findings suggest that these genes play an important role in the defense against toxic effects derived from the MP and PNP exposure. Therefore, B. cenocepacia CEIB S5-2 has a great potential for application in pesticide bioremediation approaches due to its biodegradation capabilities and the differential expression of genes for resistance to MP and PNP.

Concerns with computational protein engineering programmes IPRO and OptMAVEn and metabolic pathway engineering programme optStoic

It has become customary in engineering to require a modelling component in research endeavours. In addition, as the code for these models becomes more byzantine in complexity, it is difficult for reviewers and readers to discern their value and understand the underlying code. This opinion piece summarizes the negative experience of the author with the IPRO and OptMAVEn computational protein engineering models as well as problems with the optStoic metabolic pathway model. In our hands, these models often fail to predict reliable ways to engineer proteins and metabolic pathways.

Fungal Unspecific Peroxygenases Oxidize the Majority of Organic EPA Priority Pollutants

Unspecific peroxygenases (UPOs) are secreted fungal enzymes with promiscuity for oxygen transfer and oxidation reactions. Functionally, they represent hybrids of P450 monooxygenases and heme peroxidases; phylogenetically they belong to the family of heme-thiolate peroxidases. Two UPOs from the basidiomycetous fungi Agrocybe aegerita (AaeUPO) and Marasmius rotula (MroUPO) converted 35 out of 40 compounds listed as EPA priority pollutants, including chlorinated benzenes and their derivatives, halogenated biphenyl ethers, nitroaromatic compounds, polycyclic aromatic hydrocarbons (PAHs) and phthalic acid derivatives. These oxygenations and oxidations resulted in diverse products and—if at all—were limited for three reasons: (i) steric hindrance caused by multiple substitutions or bulkiness of the compound as such (e.g., hexachlorobenzene or large PAHs), (ii) strong inactivation of aromatic rings (e.g., nitrobenzene), and (iii) low water solubility (e.g., complex arenes). The general outcome of our study is that UPOs can be considered as extracellular counterparts of intracellular monooxygenases, both with respect to catalyzed reactions and catalytic versatility. Therefore, they should be taken into consideration as a relevant biocatalytic detoxification and biodegradation tool used by fungi when confronted with toxins, xenobiotics and pollutants in their natural environments.

Structure/function correlations over binuclear non-heme iron active sites

Binuclear non-heme iron enzymes activate O2 to perform diverse chemistries. Three different structural mechanisms of O2 binding to a coupled binuclear iron site have been identified utilizing variable-temperature, variable-field magnetic circular dichroism spectroscopy (VTVH MCD). For the μ-OH-bridged Fe(II)2 site in hemerythrin, O2 binds terminally to a five-coordinate Fe(II) center as hydroperoxide with the proton deriving from the μ-OH bridge and the second electron transferring through the resulting μ-oxo superexchange pathway from the second coordinatively saturated Fe(II) center in a proton-coupled electron transfer process. For carboxylate-only-bridged Fe(II)2 sites, O2 binding as a bridged peroxide requires both Fe(II) centers to be coordinatively unsaturated and has good frontier orbital overlap with the two orthogonal O2 π* orbitals to form peroxo-bridged Fe(III)2 intermediates. Alternatively, carboxylate-only-bridged Fe(II)2 sites with only a single open coordination position on an Fe(II) enable the one-electron formation of Fe(III)-O2 (-) or Fe(III)-NO(-) species. Finally, for the peroxo-bridged Fe(III)2 intermediates, further activation is necessary for their reactivities in one-electron reduction and electrophilic aromatic substitution, and a strategy consistent with existing spectral data is discussed.

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.

Molecular-Level Insight into the Differential Oxidase and Oxygenase Reactivities of de Novo Due Ferri Proteins

Using the single-chain due ferri (DFsc) peptide scaffold, the differential oxidase and oxygenase reactivities of two 4A→4G variants, one with two histidines at the diiron center (G4DFsc) and the other with three histidines (3His-G4DFsc(Mut3)), are explored. By controlling the reaction conditions, the active form responsible for 4-aminophenol (4-AP) oxidase activity in both G4DFsc and 3His-G4DFsc(Mut3) is determined to be the substrate-bound biferrous site. Using circular dichroism (CD), magnetic CD (MCD), and variable-temperature, variable-field (VTVH) MCD spectroscopies, 4-AP is found to bind directly to the biferrous sites of the DF proteins. In G4DFsc, 4-AP increases the coordination of the biferrous site, while in 3His-G4DFsc(Mut3), the coordination number remains the same and the substrate likely replaces the additional bound histidine. This substrate binding enables a two-electron process where 4-AP is oxidized to benzoquinone imine and O2 is reduced to H2O2. In contrast, only the biferrous 3His variant is found to be active in the oxygenation of p-anisidine to 4-nitroso-methoxybenzene. From CD, MCD, and VTVH MCD, p-anisidine addition is found to minimally perturb the biferrous centers of both G4DFsc and 3His-G4DFsc(Mut3), indicating that this substrate binds near the biferrous site. In 3His-G4DFsc(Mut3), the coordinative saturation of one iron leads to the two-electron reduction of O2 at the second iron to generate an end-on hydroperoxo-Fe(III) active oxygenating species.

Identification and characterization of a HEPN-MNT family type II toxin-antitoxin in Shewanella oneidensis

Toxin-antitoxin (TA) systems are prevalent in bacteria and archaea. However, related studies in the ecologically and bioelectrochemically important strain Shewanella oneidensis are limited. Here, we show that SO_3166, a member of the higher eukaryotes and prokaryotes nucleotide-binding (HEPN) superfamily, strongly inhibited cell growth in S. oneidensis and Escherichia coli. SO_3165, a putative minimal nucleotidyltransferase (MNT), neutralized the toxicity of SO_3166. Gene SO_3165 lies upstream of SO_3166, and they are co-transcribed. Moreover, the SO_3165 and SO_3166 proteins interact with each other directly in vivo, and antitoxin SO_3165 bound to the promoter of the TA operon and repressed its activity. Finally, the conserved Rx4-6H domain in HEPN family was identified in SO_3166. Mutating either the R or H abolished SO_3166 toxicity, confirming that Rx4-6H domain is critical for SO_3166 activity. Taken together, these results demonstrate that SO_3166 and SO_3165 in S. oneidensis form a typical type II TA pair. This TA pair plays a critical role in regulating bacterial functions because its disruption led to impaired cell motility in S. oneidensis. Thus, we demonstrated for the first time that HEPN-MNT can function as a TA system, thereby providing important insights into the understanding of the function and regulation of HEPNs and MNTs in prokaryotes.

RalR (a DNase) and RalA (a small RNA) form a type I toxin-antitoxin system in Escherichia coli

For toxin/antitoxin (TA) systems, no toxin has been identified that functions by cleaving DNA. Here, we demonstrate that RalR and RalA of the cryptic prophage rac form a type I TA pair in which the antitoxin RNA is a trans-encoded small RNA with 16 nucleotides of complementarity to the toxin mRNA. We suggest the newly discovered antitoxin gene be named ralA for RalR antitoxin. Toxin RalR functions as a non-specific endonuclease that cleaves methylated and unmethylated DNA. The RNA chaperone Hfq is required for RalA antitoxin activity and appears to stabilize RalA. Also, RalR/RalA is beneficial to the Escherichia coli host for responding to the antibiotic fosfomycin. Hence, our results indicate that cryptic prophage genes can be functionally divergent from their active phage counterparts after integration into the host genome.

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.

CAVER 3.0: a tool for the analysis of transport pathways in dynamic protein structures

Tunnels and channels facilitate the transport of small molecules, ions and water solvent in a large variety of proteins. Characteristics of individual transport pathways, including their geometry, physico-chemical properties and dynamics are instrumental for understanding of structure-function relationships of these proteins, for the design of new inhibitors and construction of improved biocatalysts. CAVER is a software tool widely used for the identification and characterization of transport pathways in static macromolecular structures. Herein we present a new version of CAVER enabling automatic analysis of tunnels and channels in large ensembles of protein conformations. CAVER 3.0 implements new algorithms for the calculation and clustering of pathways. A trajectory from a molecular dynamics simulation serves as the typical input, while detailed characteristics and summary statistics of the time evolution of individual pathways are provided in the outputs. To illustrate the capabilities of CAVER 3.0, the tool was applied for the analysis of molecular dynamics simulation of the microbial enzyme haloalkane dehalogenase DhaA. CAVER 3.0 safely identified and reliably estimated the importance of all previously published DhaA tunnels, including the tunnels closed in DhaA crystal structures. Obtained results clearly demonstrate that analysis of molecular dynamics simulation is essential for the estimation of pathway characteristics and elucidation of the structural basis of the tunnel gating. CAVER 3.0 paves the way for the study of important biochemical phenomena in the area of molecular transport, molecular recognition and enzymatic catalysis. The software is freely available as a multiplatform command-line application at http://www.caver.cz.

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.

Bacterial persistence increases as environmental fitness decreases

Since persister cells cause chronic infections and since Escherichia coli toxin MqsR increases persisters, we used protein engineering to increase the toxicity of MqsR to gain insights into persister cell formation. Through two amino acid replacements that increased the stability of MqsR, toxicity and persistence were increased. A whole‐transcriptome study revealed that the MqsR variant increased persistence by repressing genes for acid resistance, multidrug resistance and osmotic resistance. Corroborating these microarray results, deletion of rpoS, as well as the genes that the master stress response regulator RpoS controls, increased persister formation dramatically to the extent that nearly the whole population became persistent. Furthermore, wild‐type cells stressed by prior treatment to acid or hydrogen peroxide increased persistence 12 000‐fold. Whole‐transcriptome analyses of persister cells generated by two different methods (wild‐type cells pretreated with hydrogen peroxide and the rpoS deletion) corroborated the importance of suppressing RpoS in persister cell formation. Therefore, the more toxic MqsR increases persistence by decreasing the ability of the cell to respond to antibiotic stress through its RpoS‐based regulation of acid resistance, multidrug resistance and osmotic resistance systems.

Engineering global regulator Hha of Escherichia coli to control biofilm dispersal

The global transcriptional regulator Hha of Escherichia coli controls biofilm formation and virulence. Previously, we showed that Hha decreases initial biofilm formation; here, we engineered Hha for two goals: to increase biofilm dispersal and to reduce biofilm formation. Using random mutagenesis, Hha variant Hha13D6 (D22V, L40R, V42I and D48A) was obtained that causes nearly complete biofilm dispersal (96%) by increasing apoptosis without affecting initial biofilm formation. Hha13D6 caused cell death probably by the activation of proteases since Hha‐mediated dispersal was dependent on protease HslV. Hha variant Hha24E9 (K62X) was also obtained that decreased biofilm formation by inducing gadW, glpT and phnF but that did not alter biofilm dispersal. Hence, Hha may be engineered to influence both biofilm dispersal and formation.

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.

Controlling biofilm formation, prophage excision and cell death by rewiring global regulator H-NS of Escherichia coli

The global regulator H-NS of Escherichia coli controls genes related to stress response, biofilm formation and virulence by recognizing curved DNA and by silencing acquired genes. Here, we rewired H-NS to control biofilm formation using protein engineering; H-NS variant K57N was obtained that reduces biofilm formation 10-fold compared with wild-type H-NS (wild-type H-NS increases biofilm formation whereas H-NS K57N reduces it). Whole-transcriptome analysis revealed that H-NS K57N represses biofilm formation through its interaction with the nucleoid-associated proteins Cnu and StpA and in the absence of these proteins, H-NS K57N was unable to reduce biofilm formation. Significantly, H-NS K57N enhanced the excision of defective prophage Rac while wild-type H-NS represses excision, and H-NS controlled only Rac excision among the nine resident E. coli K-12 prophages. Rac prophage excision not only led to the change in biofilm formation but also resulted in cell lysis through the expression of toxin HokD. Hence, the H-NS regulatory system may be evolved through a single-amino-acid change in its N-terminal oligomerization domain to control biofilm formation, prophage excision and apoptosis.

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.

Protein engineering of the transcriptional activator FhlA To enhance hydrogen production in Escherichia coli

Escherichia coli produces H(2) from formate via the formate hydrogenlyase (FHL) complex during mixed acid fermentation; the FHL complex consists of formate dehydrogenase H (encoded by fdhF) for forming 2H(+), 2e(-), and CO(2) from formate and hydrogenase 3 (encoded by hycGE) for synthesizing H(2) from 2H(+) and 2e(-). FHL protein production is activated by the sigma(54) transcriptional activator FhlA, which activates transcription of fdhF and the hyc, hyp, and hydN-hypF operons. Here, through random mutagenesis using error-prone PCR over the whole gene, as well as over the fhlA region encoding the first 388 amino acids of the 692-amino-acid protein, we evolved FhlA to increase H(2) production. The amino acid replacements in FhlA133 (Q11H, L14V, Y177F, K245R, M288K, and I342F) increased hydrogen production ninefold, and the replacements in FhlA1157 (M6T, S35T, L113P, S146C, and E363K) increased hydrogen production fourfold. Saturation mutagenesis at the codons corresponding to the amino acid replacements in FhlA133 and at position E363 identified the importance of position L14 and of E363 for the increased activity; FhlA with replacements L14G and E363G increased hydrogen production (fourfold and sixfold, respectively) compared to FhlA. Whole-transcriptome and promoter reporter constructs revealed that the mechanism by which the FhlA133 changes increase hydrogen production is by increasing transcription of all of the genes activated by FhlA (the FHL complex). With FhlA133, transcription of P(fdhF) and P(hyc) is less sensitive to formate regulation, and with FhlA363 (E363G), P(hyc) transcription increases but P(hyp) transcription decreases and hydrogen production is less affected by the repressor HycA.

Pathways and mechanisms for product release in the engineered haloalkane dehalogenases explored using classical and random acceleration molecular dynamics simulations

Eight mutants of the DhaA haloalkane dehalogenase carrying mutations at the residues lining two tunnels, previously observed by protein X-ray crystallography, were constructed and biochemically characterized. The mutants showed distinct catalytic efficiencies with the halogenated substrate 1,2,3-trichloropropane. Release pathways for the two dehalogenation products, 2,3-dichloropropane-1-ol and the chloride ion, and exchange pathways for water molecules, were studied using classical and random acceleration molecular dynamics simulations. Five different pathways, denoted p1, p2a, p2b, p2c, and p3, were identified. The individual pathways showed differing selectivity for the products: the chloride ion releases solely through p1, whereas the alcohol releases through all five pathways. Water molecules play a crucial role for release of both products by breakage of their hydrogen-bonding interactions with the active-site residues and shielding the charged chloride ion during its passage through a hydrophobic tunnel. Exchange of the chloride ions, the alcohol product, and the waters between the buried active site and the bulk solvent can be realized by three different mechanisms: (i) passage through a permanent tunnel, (ii) passage through a transient tunnel, and (iii) migration through a protein matrix. We demonstrate that the accessibility of the pathways and the mechanisms of ligand exchange were modified by mutations. Insertion of bulky aromatic residues in the tunnel corresponding to pathway p1 leads to reduced accessibility to the ligands and a change in mechanism of opening from permanent to transient. We propose that engineering the accessibility of tunnels and the mechanisms of ligand exchange is a powerful strategy for modification of the functional 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.

Improving low-temperature activity of Sulfolobus acidocaldarius 2-keto-3-deoxygluconate aldolase

Sulfolobus acidocaldarius 2-keto-3-deoxygluconate aldolase (SacKdgA) displays optimal activity at 95 degrees C and is studied as a model enzyme for aldol condensation reactions. For application of SacKdgA at lower temperatures, a library of randomly generated mutants was screened for improved synthesis of 2-keto-3-deoxygluconate from pyruvate and glyceraldehyde at the suboptimal temperature of 50 degrees C. The single mutant SacKdgA-V193A displayed a threefold increase in activity compared with wild type SacKdgA. The increased specific activity at 40-60 degrees C of this mutant was observed, not only for the condensation of pyruvate with glyceraldehyde, but also for several unnatural acceptor aldehydes. The optimal temperature for activity of SacKdgA-V193A was lower than for the wild type enzyme, but enzymatic stability of the mutant was similar to that of the wild type, indicating that activity and stability were uncoupled. Valine193 has Van der Waals interactions with Lysine153, which covalently binds the substrate during catalysis. The mutation V193A introduced space close to this essential residue, and the increased activity of the mutant presumably resulted from increased flexibility of Lysine153. The increased activity of SacKdgA-V193A with unaffected stability demonstrates the potential for optimizing extremely thermostable aldolases for synthesis reactions at moderate temperatures.

Reconfiguring the quorum-sensing regulator SdiA of Escherichia coli to control biofilm formation via indole and N-acylhomoserine lactones

SdiA is a homolog of quorum-sensing regulators that detects N-acylhomoserine lactone (AHL) signals from other bacteria. Escherichia coli uses SdiA to reduce its biofilm formation in the presence of both AHLs and its own signal indole. Here we reconfigured SdiA (240 amino acids) to control biofilm formation using protein engineering. Four SdiA variants were obtained with altered biofilm formation, including truncation variants SdiA1E11 (F7L, F59L, Y70C, M94K, and K153X) and SdiA14C3 (W9R, P49T, N87T, frameshift at N96, and L123X), which reduced biofilm formation by 5- to 20-fold compared to wild-type SdiA in the presence of endogenous indole. Whole-transcriptome profiling revealed that wild-type SdiA reduced biofilm formation by repressing genes related to indole synthesis and curli synthesis compared to when no SdiA was expressed, while variant SdiA1E11 induced genes related to indole synthesis in comparison to wild-type SdiA. These results suggested altered indole metabolism, and corroborating the DNA microarray results in regard to indole synthesis, variant SdiA1E11 produced ninefold more indole, which led to reduced swimming motility and cell density. Also, wild-type SdiA decreased curli production and tnaA transcription, while SdiA1E11 increased tnaA transcription (tnaA encodes tryptophanase, which forms indole) compared to wild-type SdiA. Hence, wild-type SdiA decreased biofilm formation by reducing curli production and motility, and SdiA1E11 reduced biofilm formation via indole. Furthermore, an AHL-sensitive variant (SdiA2D10, having four mutations at E31G, Y42F, R116H, and L165Q) increased biofilm formation sevenfold in the presence of N-octanoyl-DL-homoserine lactone and N-(3-oxododecatanoyl)-L-homoserine lactone. Therefore, SdiA can be evolved to increase or decrease biofilm formation, and biofilm formation may be controlled by altering sensors rather than signals.

Protein engineering of toluene monooxygenases for synthesis of chiral sulfoxides

Enantiopure sulfoxides are valuable asymmetric starting materials and are important chiral auxiliaries in organic synthesis. Toluene monooxygenases (TMOs) have been shown previously to catalyze regioselective hydroxylation of substituted benzenes and phenols. Here we show that TMOs are also capable of performing enantioselective oxidation reactions of aromatic sulfides. Mutagenesis of position V106 in the alpha-hydroxylase subunit of toluene ortho-monooxygenase (TOM) of Burkholderia cepacia G4 and the analogous position I100 in toluene 4-monooxygenase (T4MO) of Pseudomonas mendocina KR1 improved both rate and enantioselectivity. Variant TomA3 V106M of TOM oxidized methyl phenyl sulfide to the corresponding sulfoxide at a rate of 3.0 nmol/min/mg protein compared with 1.6 for the wild-type enzyme, and the enantiomeric excess (pro-S) increased from 51% for the wild type to 88% for this mutant. Similarly, T4MO variant TmoA I100G increased the wild-type oxidation rate by 1.7-fold, and the enantiomeric excess rose from 86% to 98% (pro-S). Both wild-type enzymes showed lower activity with methyl para-tolyl sulfide as a substrate, but the improvement in the activity and enantioselectivity of the mutants was more dramatic. For example, T4MO variant TmoA I100G oxidized methyl para-tolyl sulfide 11 times faster than the wild type did and changed the selectivity from 41% pro-R to 77% pro-S. A correlation between regioselectivity and enantioselectivity was shown for TMOs studied in this work. Using in silico homology modeling, it is shown that residue I100 in T4MO aids in steering the substrate into the active site at the end of the long entrance channel. It is further hypothesized that the main function of V106 in TOM is the proper positioning or docking of the substrate with respect to the diiron atoms. The results from this work suggest that when the substrate is not aligned correctly in the active site, the oxidation rate is decreased and enantioselectivity is impaired, resulting in products with both chiral configurations.

Use of Pseudomonas mendocina, or recombinant Escherichia coli cells expressing toluene-4-monooxygenase, and a cell-free tyrosinase for the synthesis of 4-fluorocatechol from fluorobenzene

The transformation of fluorobenzene (FB) by whole cell expressing toluene-4-monooxygenase (T4MO) resulted in the formation of various hydroxylated products. The predominant product was either 4-fluorophenol (4FP) or 4-fluorocatechol (4Fcat) depending on the ratio of biocatalyst to substrate concentration. The transformation of 1 mM FB by whole cells (1.5 mg CDW/ml) gave a 52% yield of 4Fcat as a single product. The yield of 4Fcat was improved 1.6-fold (80%) by adding 10 mM ascorbic acid to the biotransformations. A combination of two biocatalysts (whole cells expressing T4MO and cell free mushroom tyrosinase) also resulted in the transformation of FB (5 mM) to higher concentrations of 4Fcat (1.8 mM) compared to a whole cell biotransformation alone. However, mixed products were formed and the yield of 4Fcat from FB was lower using the two-step (tandem) method (27%) compared to the use of whole cells of P. mendocina KR1 alone (80%).

Directed evolution: an approach to engineer enzymes

Directed evolution is being used increasingly in industrial and academic laboratories to modify and improve commercially important enzymes. Laboratory evolution is thought to make its biggest contribution in explorations of non-natural functions, by allowing us to distinguish the properties nurtured by evolution. In this review we report the significant advances achieved with respect to the methods of biocatalyst improvement and some critical properties and applications of the modified enzymes. The application of directed evolution has been elaborately demonstrated for protein solubility, stability and catalytic efficiency. Modification of certain enzymes for their application in enantioselective catalysis has also been elucidated. By providing a simple and reliable route to enzyme improvement, directed evolution has emerged as a key technology for enzyme engineering and biocatalysis.

Alanine 101 and alanine 110 of the alpha subunit of Pseudomonas stutzeri OX1 toluene-o-xylene monooxygenase influence the regiospecific oxidation of aromatics

Saturation mutagenesis was used to generate 10 mutants of toluene-o-xylene monooxygenase (ToMO) at alpha subunit (TouA) positions A101 and A110: A101G, A101I, A101M, A101VE, A101V, A110G, A110C, A110S, A110P, and A110T; by testing the substrates toluene, o-cresol, m-cresol, p-cresol, phenol, naphthalene, o-methoxyphenol, m-methoxyphenol, p-methoxyphenol, o-xylene, and nitrobenzene, these positions were found to influence the regiospecific oxidation of aromatics. For example, compared to wild-type ToMO, TouA variant A101V produced threefold more 3-methoxycatechol from m-methoxyphenol as well as produced methylhydroquinone from o-cresol whereas wild-type ToMO did not. Similarly, variant A110C synthesized 1.8-fold more o-cresol from toluene and 1.8-fold more 3-methoxycatechol from m-methoxyphenol, and variant A110G synthesized more m-nitrophenol and twofold less p-nitrophenol from nitrobenzene. The A101V and A110C mutations did not affect the rate of reaction with the natural substrate toluene, so the variants had high activity. This is the first report that these or analogous residues influence the catalysis with this class of enzymes. Wild-type ToMO was found to oxidize o-methoxyphenol to methoxyhydroquinone (60%) and 4-methoxyresorcinol (40%), m-methoxyphenol to 4-methoxycatechol (96%) and 3-methoxycatechol (4%), and p-methoxyphenol to 4-methoxycatechol (100%).

Regiospecific oxidation of naphthalene and fluorene by toluene monooxygenases and engineered toluene 4-monooxygenases of Pseudomonas mendocina KR1

The regiospecific oxidation of the polycyclic aromatic hydrocarbons naphthalene and fluorene was examined with Escherichia coli strains expressing wildtype toluene 4-monooxygenase (T4MO) from Pseudomonas mendocina KR1, toluene para-monooxygenase (TpMO) from Ralstonia pickettii PKO1, toluene ortho-monooxygenase (TOM) from Burkholderia cepacia G4, and toluene/ortho-xylene monooxygenase (ToMO) from P. stutzeri OX1. T4MO oxidized toluene (12.1+/-0.8 nmol/min/mg protein at 109 microM), naphthalene (7.7+/-1.5 nmol/min/mg protein at 5 mM), and fluorene (0.68+/-0.04 nmol/min/mg protein at 0.2 mM) faster than the other wildtype enzymes (2-22-fold) and produced a mixture of 1-naphthol (52%) and 2-naphthol (48%) from naphthalene, which was successively transformed to a mixture of 2,3-, 2,7-, 1,7-, and 2,6-dihydroxynaphthalenes (7%, 10%, 20%, and 63%, respectively). TOM and ToMO made 1,7-dihydroxynaphthalene from 1-naphthol, and ToMO made a mixture of 2,3-, 2,6-, 2,7-, and 1,7-dihydroxynaphthalene (26%, 22%, 1%, and 44%, respectively) from 2-naphthol. TOM had no activity on 2-naphthol, and T4MO had no activity on 1-naphthol. To take advantage of the high activity of wildtype T4MO but to increase its regiospecificity on naphthalene, seven engineered enzymes containing mutations in T4MO alpha hydroxylase TmoA were examined; the selectivity for 2-naphthol by T4MO I100A, I100S, and I100G was enhanced to 88-95%, and the selectivity for 1-naphthol was enhanced to 87% and 99% by T4MO I100L and G103S/A107G, respectively, while high oxidation rates were maintained except for G103S/A107G. Therefore, the regiospecificity for naphthalene oxidation was altered to practically pure 1-naphthol or 2-naphthol. All four wildtype monooxygenases were able to oxidize fluorene to different monohydroxylated products; T4MO oxidized fluorene successively to 3-hydroxyfluorene and 3,6-dihydroxyfluorene, which was confirmed by gas chromatography-mass spectrometry and 1H nuclear magnetic resonance analysis. TOM and its variant TomA3 V106A oxidize fluorene to a mixture of 1-, 2-, 3-, and 4-hydroxyfluorene. This is the first report of using enzymes to synthesize 1-, 3-, and 4-hydroxyfluorene, and 3,6-dihydroxyfluorene from fluorene as well as 2-naphthol and 2,6-dihydroxynaphthalene from naphthalene.

Protein engineering of toluene-o-xylene monooxygenase from Pseudomonas stutzeri OX1 for oxidizing nitrobenzene to 3-nitrocatechol, 4-nitrocatechol, and nitrohydroquinone

Toluene-o-xylene monooxygenase (ToMO) from Pseudomonas stutzeri OX1 was found to oxidize nitrobenzene (NB) to form m-nitrophenol (m-NP, 72%) and p-NP (28%) with an initial rate of 0.098 and 0.031 nmol/(min mg protein), respectively. It was also discovered that wild-type ToMO forms 4-nitrocatechol (4-NC) from m-NP and p-NP with an initial rate of 0.15 and 0.0082 nmol/(min mg protein), respectively, and 3-NC (12%) and nitrohydroquinone (NHQ, 88%) from o-NP with an initial rate of 0.11 and 0.8 nmol/(min mg protein), respectively. To increase the oxidation rate and alter the oxidation regiospecificity of nitro aromatics as well as to study the role of the active site residues I100, Q141, T201, and F205 of the alpha hydroxylase fragment of ToMO (TouA), DNA shuffling and saturation mutagenesis were used to generate random mutants. The mutants were initially identified by screening via a rapid agar plate assay and then were further examined by high-performance liquid chromatography (HPLC) and gas chromatography (GC). Several mutants with higher rates of activities and with different regiospecificities were identified; for example, Escherichia coli TG1 cells expressing either TouA mutant M180T/E284G or E214G/D312N/M399V produce 4-NC 4.5- and 20-fold faster than wild-type ToMO (0.037 and 0.16 nmol/min mg protein from p-NP, respectively). TouA mutant A107T/E214A had the regiospecificity of NB changed significantly from 28% to 79% p-NP. From 200 microM NB, TouA variants A101T/M114T, A110T/E392D, M180T/E284G, and E214G/D312N/M399V produce 4-NC whereas wild-type ToMO does not. From m-NP, TouA mutant I100Q produces 4-NC (37%) and NHQ (63%), whereas wild-type ToMO produces only 4-NC (100%). Variant A107T/E214A acts like a para enzyme and forms p-cresol as the major product (93%) from toluene with enhanced activity (2.3-fold), whereas wild-type ToMO forms 32%, 21%, and 47% of o-, m-, and p-cresol, respectively. Hence, the non-specific ToMO was converted into a regiospecific enzyme, which rivals toluene 4-monooxygenase of P. mendocina KR1 and toluene o-monooxygenase of Burkholderia cepacia G4 in its specificity.

Alpha-subunit positions methionine 180 and glutamate 214 of Pseudomonas stutzeri OX1 toluene-o-xylene monooxygenase influence catalysis

Alpha-subunit position M180 of toluene-o-xylene monooxygenase influences the regiospecific oxidation of aromatics (e.g., from o-cresol, M180H forms 3-methylcatechol, methylhydroquinone, and 4-methylresorcinol, whereas the wild type forms only 3-methylcatechol). Position E214 influences the rate of reaction (e.g., E214G increases p-nitrophenol oxidation 15-fold) by controlling substrate entrance and product efflux as a gate residue.

Controlling the regiospecific oxidation of aromatics via active site engineering of toluene para-monooxygenase of Ralstonia pickettii PKO1

A primary goal of protein engineering is to control catalytic activity. Here we show that through mutagenesis of three active site residues, the catalytic activity of a multicomponent monooxygenase is altered so that it hydroxylates all three positions of toluene as well as both positions of naphthalene. Hence, for the first time, an enzyme has been engineered so that its regiospecific oxidation of a substrate can be controlled. Through the A107G mutation in the alpha-subunit of toluene para-monooxygenase, a variant was formed that hydroxylated toluene primarily at the ortho-position while converting naphthalene to 1-naphthol. Conversely, the A107T variant produced >98% p-cresol and p-nitrophenol from toluene and nitrobenzene, respectively, as well as produced 2-naphthol from naphthalene. The mutation I100S/G103S produced a toluene para-monooxygenase variant that formed 75% m-cresol from toluene and 100% m-nitrophenol from nitrobenzene; thus, for the first time a true meta-hydroxylating toluene monooxygenase was created.