The transition of human estrogen sulfotransferase from generalist to specialist using directed enzyme evolution

Broad specificity is believed to be a property of primordial enzymes that diverged during natural protein evolution to produce highly specific and efficient enzymes. Human estrogen sulfotransferase (SULT1E1) is a broad-specificity enzyme that detoxifies a variety of chemicals, including estrogens, by the transfer of sulfate. To study the molecular basis for the broad specificity of this enzyme and to investigate the process of SULT1E1 specialization, we have adopted a directed enzyme evolution approach. Using two iterative rounds of evolution, we generated SULT1E1 mutants with increased thermostability and narrower specificity from the broadly specific wild-type enzyme. To identify mutants with enhanced specificity, we developed an unbiased screening assay to assess sulfate transfer to three different acceptors in parallel. Such an assay enabled the isolation of SULT1E1 mutants with enhanced or wild-type activity toward an estrogen acceptor and significantly reduced activity for phenol or coumarin type of acceptors, leading to up to 3 orders of magnitude increase in specificity. We found that mutations conferring novel specificity are located in the vicinity of the active site and thus may play a direct role in reshaping the acceptor-binding site. Finally, such mutations resulted in reduced SULT1E1 thermostability, revealing a trade-off between SULT1E1 thermostability and acquisition of novel function.     

Directed evolution of the thermostable xylanase from Thermomyces lanuginosus

The thermostability of the endo-beta-1,4-xylanase from Thermomyces lanuginosus (xynA) was improved by directed evolution using error-prone PCR. Transformants expressing the variant xylanases were first selected on 0.4% Remazol Brilliant Blue-xylan and then exposed to 80 degrees C. Whereas the wild type XynA lost 90% activity after 10 min at 80 degrees C, five mutants displayed both higher stabilities and activities than XynA. Four mutants were subjected to further mutagenesis to improve the stability and activity of the xylanase. Subsequent screening revealed three mutants with enhanced thermostability. Mutant 2B7-10 retained 71% of its activity after treatment at 80 degrees C for 60 min and had a half-life of 215 min at 70 degrees C, which is higher than that attained by XynA. Sequence analysis of second generation mutants revealed that mutations were not concentrated in any particular region of the protein and exhibited much variation. The best mutant obtained from this study was variant 2B7-10, which had a single substitution (Y58F) in beta-sheet A of the protein, which is the hydrophilic, solvent-accessible outer surface of the enzyme. Most of the mutants obtained in this study displayed a compromise between stability and activity, the only exception being mutant 2B7-10. This variant showed increased activity and thermostability.     

Improved thermostability of Clostridium thermocellum endoglucanase Cel8A by using consensus-guided mutagenesis

The use of thermostable cellulases is advantageous for the breakdown of lignocellulosic biomass toward the commercial production of biofuels. Previously, we have demonstrated the engineering of an enhanced thermostable family 8 cellulosomal endoglucanase (EC 3.2.1.4), Cel8A, from Clostridium thermocellum, using random error-prone PCR and a combination of three beneficial mutations, dominated by an intriguing serine-to-glycine substitution (M. Anbar, R. Lamed, E. A. Bayer, ChemCatChem 2:997-1003, 2010). In the present study, we used a bioinformatics-based approach involving sequence alignment of homologous family 8 glycoside hydrolases to create a library of consensus mutations in which residues of the catalytic module are replaced at specific positions with the most prevalent amino acids in the family. One of the mutants (G283P) displayed a higher thermal stability than the wild-type enzyme. Introducing this mutation into the previously engineered Cel8A triple mutant resulted in an optimized enzyme, increasing the half-life of activity by 14-fold at 85°C. Remarkably, no loss of catalytic activity was observed compared to that of the wild-type endoglucanase. The structural changes were simulated by molecular dynamics analysis, and specific regions were identified that contributed to the observed thermostability. Intriguingly, most of the proteins used for sequence alignment in determining the consensus residues were derived from mesophilic bacteria, with optimal temperatures well below that of C. thermocellum Cel8A.     

Simultaneous enhancement of thermostability and catalytic activity of phospholipase A(1) by evolutionary molecular engineering

The thermal stability and catalytic activity of phospholipase A(1) from Serratia sp. strain MK1 were improved by evolutionary molecular engineering. Two thermostable mutants were isolated after sequential rounds of error-prone PCR performed to introduce random mutations and filter-based screening of the resultant mutant library; we determined that these mutants had six (mutant TA3) and seven (mutant TA13) amino acid substitutions. Different types of substitutions were found in the two mutants, and these substitutions resulted in an increase in nonpolar residues (mutant TA3) or in differences between side chains for polar or charged residues (mutant TA13). The wild-type and mutant enzymes were purified, and the effect of temperature on the stability and catalytic activity of the enzymes was investigated. The melting temperatures of the TA3 and TA13 enzymes were increased by 7 and 11 degrees C, respectively, compared with the melting temperature of the wild-type enzyme. Thus, we found that evolutionary molecular engineering was an effective and efficient approach for increasing thermostability without compromising enzyme activity.     

蛋白质定向进化的研究进展

定向进化是改造蛋白质分子的一种有效的新策略。主要是在实验室里模拟自然进化过程,通过由易错PCR、致突变菌株诱变等方法对编码蛋白质的基因进行随机诱变,由DNA改组、随机引导重组和交错延伸等方法进行突变基因体外重组,设计高通量筛选方法来选出需要的突变株。它不仅可快速产生工业上有用的新酶,而且对研究蛋白质的结构与功能的关系具有非常重要的意义。     

Conversion of a cyclodextrin glucanotransferase into an alpha-amylase: assessment of directed evolution strategies

Glycoside hydrolase family 13 (GH13) members have evolved to possess various distinct reaction specificities despite the overall structural similarity. In this study we investigated the evolutionary input required to effeciently interchange these specificities and also compared the effectiveness of laboratory evolution techniques applied, i.e., error-prone PCR and saturation mutagenesis. Conversion of our model enzyme, cyclodextrin glucanotransferase (CGTase), into an alpha-amylase like hydrolytic enzyme by saturation mutagenesis close to the catalytic core yielded a triple mutant (A231V/F260W/F184Q) with the highest hydrolytic rate ever recorded for a CGTase, similar to that of a highly active alpha-amylase, while cyclodextrin production was virtually abolished. Screening of a much larger, error-prone PCR generated library yielded far less effective mutants. Our results demonstrate that it requires only three mutations to change CGTase reaction specificity into that of another GH13 enzyme. This suggests that GH13 members may have diversified by introduction of a limited number of mutations to the common ancestor, and that interconversion of reaction specificites may prove easier than previously thought.     

Thermosensitive mutants of the MPTP and hPTP1B protein tyrosine phosphatases: isolation and structural analysis.

A PCR-based random mutagenesis procedure was employed to identify several thermosensitive mutants of the MPTP enzyme, the murine homologue of the human T-cell PTPase and rat PTP-S enzymes. Four mutants with varying degrees of thermosensitivity were characterized for their thermostability and refolding properties following incubation at the nonpermissive temperature. Structure analysis of these mutations based on the hPTP1B co-ordinate structure demonstrates a clear relationship between the position of each mutated amino acid relative to the catalytic cysteine residue and their thermostability. Introduction of two of these mutations in the related enzyme hPTP1B suggests that the structural defects and the resulting thermosensitivity of these mutations may represent an intrinsic property of all PTPase catalytic domains.     

Simultaneous Enhancement of Thermostability and Catalytic Activity of Phospholipase A1 by Evolutionary Molecular Engineering

The thermal stability and catalytic activity of phospholipase A₁ from Serratia sp. strain MK1 were improved by evolutionary molecular engineering. Two thermostable mutants were isolated after sequential rounds of error-prone PCR performed to introduce random mutations and filter-based screening of the resultant mutant library; we determined that these mutants had six (mutant TA3) and seven (mutant TA13) amino acid substitutions. Different types of substitutions were found in the two mutants, and these substitutions resulted in an increase in nonploar residues (mutant TA3) or in differences between side chains for polar or charged residues (mutant TA13). The wild-type and mutant enzymes were purified, and the effect of temperature on the stability and catalytic activity of the enzymes was investigated. The melting temperatures of the TA3 and TA13 enzymes were increased by 7 and 11°C, respectively, compared with the melting temperature of the wild-type enzyme. Thus, we found that evolutionary molecular engineering was an effective and efficient approach for increasing thermostability without compromising enzyme activity.     

Directed evolution of alpha-aspartyl dipeptidase from Salmonella typhimurium.

Model-free approaches (error-prone PCR to introduce random mutations, DNA shuffling to combine positive mutations, and screening of the resultant mutant libraries) have been used to enhance the catalytic activity and thermostability of alpha-aspartyl dipeptidase from Salmonella typhimurium, which is uniquely able to hydrolyze Asp-X dipeptides (where X is any amino acid) and one tripeptide (Asp-Gly-Gly). Under double selective pressures of activity and thermostability, through two rounds of error-prone PCR and three sequential generations of DNA shuffling, coupled with screening, a mutant pepEM3074 with approximately 47-fold increased enzyme activity compared with its wild-type parent was obtained. Moreover, the stability of pepEM3074 is increased significantly. Three amino acid substitutions (Asn89His, Gln153Glu, and Leu205Arg), two of them are near the active site and substrate binding pocket, were identified by sequencing the genes encoding this evolved enzyme. The mechanism of the enhancement of activity and stability was analyzed in this paper.     

Gene-specific random mutagenesis of Escherichia coli in vivo: isolation of temperature-sensitive mutations in the acyl carrier protein of fatty acid synthesis

Acyl carrier proteins (ACPs) are very small acidic proteins that play a key role in fatty acid and complex lipid synthesis. Moreover, recent data indicate that the acyl carrier protein of Escherichia coli has a large protein interaction network that extends beyond lipid synthesis. Despite extensive efforts over many years, no temperature-sensitive mutants with mutations in the structural gene (acpP) that encodes ACP have been isolated. We report the isolation of three such mutants by a new approach that utilizes error-prone PCR mutagenesis, overlap extension PCR, and phage lambda Red-mediated homologous recombination and that should be generally applicable. These mutants plus other experiments demonstrate that ACP function is essential for the growth of E. coli. Each of the mutants was efficiently modified with the phosphopantetheinyl moiety essential for the function of ACP in lipid synthesis, and thus lack of function at the nonpermissive temperature cannot be attributed to a lack of prosthetic group attachment. All of the mutant proteins were largely stable at the nonpermissive temperature except the A68T/N73D mutant protein. Fatty acid synthesis in strains that carried the D38V or A68T/N73D mutations was inhibited upon a shift to the nonpermissive temperature and in the latter case declined to a small percentage of the rate of the wild-type strain.     

Low-fidelity Pyrococcus furiosus DNA polymerase mutants useful in error-prone PCR

Random mutagenesis constitutes an important approach for identifying critical regions of proteins, studying structure-function relations and developing novel proteins with desired properties. Perhaps, the most popular method is the error-prone PCR, in which mistakes are introduced into a gene, and hence a protein, during DNA polymerase-catalysed amplification cycles. Unfortunately, the relatively high fidelities of the thermostable DNA polymerases commonly used for PCR result in too few mistakes in the amplified DNA for efficient mutagenesis. In this paper, we describe mutants of the family B DNA polymerase from Pyrococcus furiosus (Pfu-Pol), with superb performance in error-prone PCR. The key amino acid changes occur in a short loop linking two long α-helices that comprise the ‘fingers’ sub-domain of the protein. This region is responsible for binding the incoming dNTPs and ensuring that only correct bases are inserted opposite the complementary base in the template strand. Mutations in the short loop, when combined with an additional mutation that abolishes the 3′–5′ proof-reading exonuclease activity, convert the extremely accurate wild-type polymerase into a variant with low fidelity. The mutant Pfu-Pols can be applied in error-prone PCR, under exactly the same conditions used for standard, high-fidelity PCR with the wild-type enzyme. Large quantities of amplified product, with a high frequency of nearly indiscriminate mutations, are produced. It is anticipated that the Pfu-Pol variants will be extremely useful for the randomization of gene, and hence protein, sequences.     

Exploring the thermostable properties of halohydrin dehalogenase from Agrobacterium radiobacter AD1 by a combinatorial directed evolution strategy

As a crucial factor for biocatalysts, protein thermostability often arises from a combination of factors that are often difficult to rationalize. In this work, the thermostable nature of halohydrin dehalogenase from Agrobacterium radiobacter AD1 (HheC) was systematically explored using a combinatorial directed evolution approach. For this, a mutagenesis library of HheC mutants was first constructed using error-prone PCR with low mutagenesis frequency. After screening approximately 2000 colonies, six mutants with eight mutation sites were obtained. Those mutation sites were subsequently combined by adopting several rounds of iterative saturation mutagenesis (ISM) approach. After four rounds of saturation mutagenesis, one best mutant ISM-4 with a 3400-fold improvement in half-life (t _1/2) inactivation at 65 °C, 18 °C increase in apparent T _m value, and 20 °C increase in optimum temperature was obtained, compared to wild-type HheC. To the best of our knowledge, the mutant represents the most thermostable HheC variant reported up to now. Moreover, the mutant was as active as wild-type enzyme for the substrate 1,3-dichloro-2-propanol, and they remained most enantioselectivity of wild-type enzyme in the kinetic resolution of rac-2-chloro-1-phenolethanol, exhibiting a great potential for industrial applications. Our structural investigation highlights that surface loop regions are hot spots for modulating the thermostability of HheC, the residues located at these regions contribute to the thermostability of HheC in a cooperative way, and protein rigidity and oligomeric interface connections contribute to the thermostability of HheC. All of these essential factors could be used for further design of an even more thermostable HheC, which, in turn, could greatly facilitate the application of the enzyme as a biocatalyst.

     

Thermostable Bacillus subtilis lipases: in vitro evolution and structural insight

In vitro evolution methods are now being routinely used to identify protein variants with novel and enhanced properties that are difficult to achieve using rational design. However, one of the limitations is in screening for beneficial mutants through several generations due to the occurrence of neutral/negative mutations occurring in the background of positive ones. While evolving a lipase in vitro from mesophilic Bacillus subtilis to generate thermostable variants, we have designed protocols that combine stringent three-tier testing, sequencing and stability assessments on the protein at the end of each generation. This strategy resulted in a total of six stabilizing mutations in just two generations with three mutations per generation. Each of the six mutants when evaluated individually contributed additively to thermostability. A combination of all of them resulted in the best variant that shows a remarkable 15 °C shift in melting temperature and a millionfold decrease in the thermal inactivation rate with only a marginal increase of 3 kcal mol⁻¹ in free energy of stabilization. Notably, in addition to the dramatic shift in optimum temperature by 20 °C, the activity has increased two- to fivefold in the temperature range 25-65 °C. High-resolution crystal structures of three of the mutants, each with 5° increments in melting temperature, reveal the structural basis of these mutations in attaining higher thermostability. The structures highlight the importance of water-mediated ionic networks on the protein surface in imparting thermostability. Saturation mutagenesis at each of the six positions did not result in enhanced thermostability in almost all the cases, confirming the crucial role played by each mutation as revealed through the structural study. Overall, our study presents an efficient strategy that can be employed in directed evolution approaches employed for obtaining improved properties of proteins.     

Improving tolerance of Candida antarctica lipase B towards irreversible thermal inactivation through directed evolution

To expand the functionality of lipase B from Candida antarctica (CALB) we have used directed evolution to create CALB mutants with improved resistance towards irreversible thermal inactivation. Two mutants, 23G5 and 195F1, were generated with over a 20-fold increase in half-life at 70 degrees C compared with the wild-type CALB (WT-CALB). The increase in half-life was attributed to a lower propensity of the mutants to aggregate in the unfolded state and to an improved refolding. The first generation mutant, 23G5, obtained by error-prone PCR, had two amino acid mutations, V210I and A281E. The second generation mutant, 195F1, derived from 23G5 by error-prone PCR, had one additional mutation, V221D. Amino acid substitutions at positions 221 and 281 were determined to be critical for lipase stability, while the residue at position 210 had only a marginal effect. The catalytic efficiency of the mutants with p-nitrophenyl butyrate and 6,8-difluoro-4-methylumbelliferyl octanoate was also found to be superior to that of WT-CALB.     

Identifying functionally important mutations from phenotypically diverse sequence data

Here we present a simple statistical method to determine the phenotypic contribution of a single mutation from libraries of mutants with diverse phenotypes in which each mutant contains a multitude of mutations. The central premise of this method is that, given M phenotypic classes, mutations that do not affect the phenotype should partition among the M classes according to a multinomial distribution. Deviations from this distribution are indicative of a link between specific mutations and phenotypes. We suggest that this method will aid the engineering of functional nucleic acids, proteins, and other biomolecules by uncovering target sites for rational mutagenesis. As a proof of the principle, we show how the method can be used to deduce the individual effects of mutations in a set of 69 P(L)-lambda promoter variants. Each of these promoters was generated by error-prone PCR and incorporated numerous mutations. The activity of the promoters was assayed using flow cytometry to measure the fluorescence of a green fluorescent protein reporter gene. Our analysis of the sequences of these mutants revealed seven positions having a statistically significant correlation with promoter activity. Using site-directed mutagenesis, we constructed point mutations for several sites, both statistically significant and insignificant, and combinations of these sites. Our results show that the statistical method correctly elucidated the phenotypic manifestations of these mutations. We suggest that this method may be useful for expediting directed evolution experiments by allowing both desired and undesired mutations to be identified and incorporated between rounds of mutagenesis.     

Application of Rigidity Theory to the Thermostabilization of Lipase A from Bacillus subtilis

Protein thermostability is a crucial factor for biotechnological enzyme applications. Protein engineering studies aimed at improving thermostability have successfully applied both directed evolution and rational design. However, for rational approaches, the major challenge remains the prediction of mutation sites and optimal amino acid substitutions. Recently, we showed that such mutation sites can be identified as structural weak spots by rigidity theory-based thermal unfolding simulations of proteins. Here, we describe and validate a unique, ensemble-based, yet highly efficient strategy to predict optimal amino acid substitutions at structural weak spots for improving a protein’s thermostability. For this, we exploit the fact that in the majority of cases an increased structural rigidity of the folded state has been found as the cause for thermostability. When applied prospectively to lipase A from Bacillus subtilis, we achieved both a high success rate (25% over all experimentally tested mutations, which raises to 60% if small-to-large residue mutations and mutations in the active site are excluded) in predicting significantly thermostabilized lipase variants and a remarkably large increase in those variants’ thermostability (up to 6.6°C) based on single amino acid mutations. When considering negative controls in addition and evaluating the performance of our approach as a binary classifier, the accuracy is 63% and increases to 83% if small-to-large residue mutations and mutations in the active site are excluded. The gain in precision (predictive value for increased thermostability) over random classification is 1.6-fold (2.4-fold). Furthermore, an increase in thermostability predicted by our approach significantly points to increased experimental thermostability (p < 0.05). These results suggest that our strategy is a valuable complement to existing methods for rational protein design aimed at improving thermostability.     

Change of Bacillus cereus flavonoid O-triglucosyltransferase into flavonoid O-monoglucosyltransferase by error-prone polymerase chain reaction

The attachment of sugar to flavonoids enhances their solubility. Glycosylation is performed primarily by uridine diphosphate-dependent glycosyltransferases (UGTs). The UGT from Bacillus cereus, BcGT-1 transferred three glucose molecules into kaempferol. The structural analysis of BcGT-1 showed that its substrate binding site is wider than that of flavonoid monoglucosyltransferase of plant. In order to create monoglucosyltransferase from BcGT-1, error-prone polymerase chain reaction (PCR) was performed. We analyzed 150 clones. Among them, two mutants generated only kaempferol O-monoglucoside, albeit with reduced reactivity. Unexpectedly, the two mutants harbored mutations in the amino acids located outside of the active sites. Based on the modeled structure of BcGT-1, it was proposed that the local change in the secondary structure of BcGT-1 caused the alteration of triglucosyltransferase into monoglucosyltransferase.     

Directed evolution of acyl-CoA:diacylglycerol acyltransferase: Development and characterization of Brassica napus DGAT1 mutagenized libraries

Metabolic flux to triacylglycerol (TAG) may be limited by the level of acyl-CoA:diacylglycerol acyltransferase (DGAT, EC 2.3.1.20) activity. In some species, this enzyme also appears to play a role in the channeling of specific fatty acyl moieties into TAG. The objective of this work is to implement a directed evolution approach to enhance the catalytic efficiency of type-1 DGAT from Brassica napus (BnDGAT1). We generated randomly mutagenized libraries of BnDGAT1 in a yeast expression vector using error-prone PCR. The mutagenized libraries were used to transform a Saccharomyces cerevisiae strain devoid of neutral lipid biosynthesis and analyzed using a high-throughput screening (HTS) system. The HTS, recently developed for this purpose, consisted of a positive selection of clones expressing active DGAT mutants followed by quantification of DGAT activity by fluorescence detection of TAG in yeast cells. The initial results indicated that the positive selection system efficiently eliminated DGAT mutants lacking enzyme activity. Screening of 1528 selected mutants revealed that some DGAT clones had enhanced ability to synthesize TAG in yeast. This was confirmed by analysis of individual clones that could carry mutations resulting in an increased catalytic efficiency. The directed evolution approach could lead to the development of an improved plant DGAT1 for increasing seed oil content in oleaginous crops.     

Thermostability of In Vitro Evolved Bacillus subtilis Lipase A: A Network and Dynamics Perspective

Proteins in thermophilic organisms remain stable and function optimally at high temperatures. Owing to their important applicability in many industrial processes, such thermostable proteins have been studied extensively, and several structural factors attributed to their enhanced stability. How these factors render the emergent property of thermostability to proteins, even in situations where no significant changes occur in their three-dimensional structures in comparison to their mesophilic counter-parts, has remained an intriguing question. In this study we treat Lipase A from Bacillus subtilis and its six thermostable mutants in a unified manner and address the problem with a combined complex network-based analysis and molecular dynamic studies to find commonality in their properties. The Protein Contact Networks (PCN) of the wild-type and six mutant Lipase A structures developed at a mesoscopic scale were analyzed at global network and local node (residue) level using network parameters and community structure analysis. The comparative PCN analysis of all proteins pointed towards important role of specific residues in the enhanced thermostability. Network analysis results were corroborated with finer-scale molecular dynamics simulations at both room and high temperatures. Our results show that this combined approach at two scales can uncover small but important changes in the local conformations that add up to stabilize the protein structure in thermostable mutants, even when overall conformation differences among them are negligible. Our analysis not only supports the experimentally determined stabilizing factors, but also unveils the important role of contacts, distributed throughout the protein, that lead to thermostability. We propose that this combined mesoscopic-network and fine-grained molecular dynamics approach is a convenient and useful scheme not only to study allosteric changes leading to protein stability in the face of negligible over-all conformational changes due to mutations, but also in other molecular networks where change in function does not accompany significant change in the network structure.     

A multi-factors rational design strategy for enhancing the thermostability of Escherichia coli AppA phytase

Despite recent advances in our understanding of the importance of protein surface properties for protein thermostability,there are seldom studies on multi-factors rational design strategy, so a more scientific, simple and effective rational strategy is urgent for protein engineering. Here, we first attempted to use a three-factors rational design strategy combining three common structural features, protein flexibility, protein surface, and salt bridges. Escherichia coli AppA phytase was used as a model enzyme to improve its thermostability. Moreover, the structure and enzyme features of the thermostable mutants designed by our strategy were analyzed roundly. For the single mutants, two (Q206E and Y311K), in five exhibited thermostable property with a higher success rate of prediction (40 %). For the multiple mutants, the themostable sites were combined with another site, I427L, we obtained by directed evolution, Q206E/I427L, Y311K/I427L, and Q206E/Y311K/I427L, all exhibited thermostable property. The Y311K/I427L doubled thermostability (61.7 %, and was compared to 30.97 % after being heated at 80 °C for 10 min) and catalytic efficiency (4.46 was compared to 2.37) improved more than the wild-type AppA phytase almost without hampering catalytic activity. These multi-factors of rational design strategy can be applied practically as a thermostabilization strategy instead of the conventional single-factor approach.