Directed evolution of enzyme stability

Enhancing the efficiency of directed evolution in focused enzyme libraries by the adaptive substituent reordering algorithm

Directed evolution is a broadly successful strategy for protein engineering in the quest to enhance the stereoselectivity, activity, and thermostability of enzymes. To increase the efficiency of directed evolution based on iterative saturation mutagenesis, the adaptive substituent reordering algorithm (ASRA) is introduced here as an alternative to traditional quantitative structure-activity relationship (QSAR) methods for identifying potential protein mutants with desired properties from minimal sampling of focused libraries. The operation of ASRA depends on identifying the underlying regularity of the protein property landscape, allowing it to make predictions without explicit knowledge of the structure-property relationships. In a proof-of-principle study, ASRA identified all or most of the best enantioselective mutants among the synthesized epoxide hydrolase from Aspergillus niger, in the absence of peptide seeds with high E-values. ASRA even revealed a laboratory error from irregularities of the reordered E-value landscape alone.

Biophysical characterization of mutants of Bacillus subtilis lipase evolved for thermostability: factors contributing to increased activity retention

Previously, Lipase A from Bacillus subtilis was subjected to in vitro directed evolution using iterative saturation mutagenesis, with randomization sites chosen on the basis of the highest B-factors available from the crystal structure of the wild-type (WT) enzyme. This provided mutants that, unlike WT enzyme, retained a large part of their activity after heating above 65 °C and cooling down. Here, we subjected the three best mutants along with the WT enzyme to biophysical and biochemical characterization. Combining thermal inactivation profiles, circular dichroism, X-ray structure analyses and NMR experiments revealed that mutations of surface amino acid residues counteract the tendency of Lipase A to undergo precipitation under thermal stress. Reduced precipitation of the unfolding intermediates rather than increased conformational stability of the evolved mutants seems to be responsible for the activity retention.

Thermostabilizing mutations preferentially occur at structural weak spots with a high mutation ratio

We apply Constraint Network Analysis (CNA) to investigate the relationship between structural rigidity and thermostability of five citrate synthase (CS) structures over a temperature range from 37 °C to 100 °C. For the first time, we introduce an ensemble-based variant of CNA and model the temperature-dependence of hydrophobic interactions in the constraint network. A very good correlation between the predicted thermostabilities of CS and optimal growth temperatures of their source organisms (R²=0.88, p=0.017) is obtained, which validates that CNA is able to quantitatively discriminate between less and more thermostable proteins even within a series of orthologs. Structural weak spots on a less thermostable CS, predicted by CNA to be in the top 5% with respect to the frequency of occurrence over an ensemble, have a higher mutation ratio in a more thermostable CS than other sequence positions. Furthermore, highly ranked weak spots that are also highly conserved with respect to the amino acid type found at that sequence position are nevertheless found to be mutated in the more stable CS. As for mechanisms at an atomic level that lead to a reinforcement of weak spots in more stable CS, we observe that the thermophilic CS achieve a higher thermostability by better hydrogen bonding networks whereas hyperthermophilic CS incorporate more hydrophobic contacts to reach the same goal. Overall, these findings suggest that CNA can be applied as a pre-filter in data-driven protein engineering to focus on residues that are highly likely to improve thermostability upon mutation.

Restriction enzyme-free construction of random gene mutagenesis libraries in Escherichia coli

Directed evolution relies on both random and site-directed mutagenesis of individual genes and regulatory elements to create variants with altered activity profiles for engineering applications. Central to these experiments is the construction of large libraries of related variants. However, a number of technical hurdles continue to limit routine construction of random mutagenesis libraries in Escherichia coli, in particular, inefficiencies during digestion and ligation steps. Here, we report a restriction enzyme-free approach to library generation using megaprimers termed MegAnneal. Target DNA is first exponentially amplified using error-prone polymerase chain reaction (PCR) and then linearly amplified with a single 3' primer to generate long, randomly mutated, single-stranded megaprimers. These are annealed to single-stranded dUTP-containing template plasmid and extended with T7 polymerase to create a complementary strand, and the resulting termini are ligated with T4 DNA ligase. Using this approach, we are able to reliably generate libraries of approximately 10⁷ colony-forming units (cfu)/μg DNA/transformation in a single day. We have created MegAnneal libraries based on three different single-chain antibodies and identified variants with enhanced expression and ligand-binding affinity. The key advantages of this approach include facile amplification, restriction enzyme-free library generation, and a significantly reduced risk of mutations outside the targeted region and wild-type contamination as compared with current methods.

In vitro evolved non-aggregating and thermostable lipase: structural and thermodynamic investigation

Rational and in vitro evolutionary approaches to improve either protein stability or aggregation resistance were successful, but empirical rules for simultaneous improvement of both stability and aggregation resistance under denaturing conditions are still to be ascertained. We have created a robust variant of a lipase from Bacillus subtilis named "6B" using multiple rounds of in vitro evolution. T(m) and optimum activity temperature of 6B is 78 °C and 65 °C, respectively, which is ~22 °C and 30 °C higher than that of wild-type lipase. Most significantly, 6B does not aggregate upon heating. Physical basis of remarkable thermostability and non-aggregating behavior of 6B was explored using X-ray crystallography, NMR and differential scanning calorimetry. Our structural investigations highlight the importance of tightening of mobile regions of the molecule such as loops and helix termini to attain higher thermostability. Accordingly, NMR studies suggest a very rigid structure of 6B lipase. Further investigation suggested that reduction/perturbation of the large hydrophobic patches present in the wild-type protein structure, decreased propensity of amino acid sequence for aggregation and absence of aggregation-prone intermediate during thermal unfolding of 6B can account for its resistance to aggregation. Overall, our study suggest that better anchoring of the loops with the rest of the protein molecule through mutations particularly on the sites that perturb/disturb the exposed hydrophobic patches can simultaneously increase protein stability and aggregation resistance.

Thermostabilization of firefly luciferase by in vivo directed evolution

Firefly luciferase is widely used in a number of areas of biotechnology and molecular biology. However, rapid inactivation of wild-type (WT) luciferases at elevated temperatures often hampers their application. A simple non-lethal in vivo screening scheme was used to identify thermostable mutants of luciferase in Escherichia coli colonies. This scheme allowed carrying out each cycle of mutagenesis in a rapid and efficient manner. Four rounds of directed evolution were conducted on a part of the gene coding for amino acid residues 130-390 of Luciola mingrelica luciferase. The resultant mutant designated 4TS had a half-life of 10 h at 42°C, which is 65-fold higher compared with the WT luciferase. Moreover, the mutant 4TS showed a 1.9-fold increase in specific activity, 5.7-fold reduction of K(m) for ATP and a higher-temperature optimum compared with the WT enzyme. 4TS contains eight mutations, four of which are suggested to be mainly responsible for the enhancement of thermostability: R211L, A217V, E356K and S364C. Thus, directed evolution with non-lethal colony screening for in vivo bioluminescence activity proved to be an effective and efficient approach for increasing thermal stability of luciferase while retaining high catalytic activity.

The role of calcium ions in the stability and instability of a thermolysin-like protease

Thermolysin and other secreted broad-specificity proteases, such as subtilisin or alpha-lytic protease, are produced as pre-pro-proteins that stay at least partially unfolded while in the cytosol. After secretion, the pro-proteases fold to their active conformations in a process that includes the autolytic removal of the pro-peptide. We review the life cycle of the thermolysin-like protease from Bacillus stearothermophilus in light of the calcium dependent stability and instability of the N-terminal domain. The protease binds calcium ions in the regions that are involved in the autolytic maturation process. It is generally assumed that the calcium ions contribute to the extreme stability of the protease, but experimental evidence for TLP-ste indicates that at least one of the calcium ions plays a regulatory role. We hypothesize that this calcium ion plays an important role as a switch that modulates the protease between stable and unstable states as appropriate to the biological need.

A rewired green fluorescent protein: folding and function in a nonsequential, noncircular GFP permutant

The sequential order of secondary structural elements in proteins affects the folding and activity to an unknown extent. To test the dependence on sequential connectivity, we reconnected secondary structural elements by their solvent-exposed ends, permuting their sequential order, called "rewiring". This new protein design strategy changes the topology of the backbone without changing the core side chain packing arrangement. While circular and noncircular permutations have been observed in protein structures that are not related by sequence homology, to date no one has attempted to rationally design and construct a protein with a sequence that is noncircularly permuted while conserving three-dimensional structure. Herein, we show that green fluorescent protein can be rewired, still functionally fold, and exhibit wild-type fluorescence excitation and emission spectra.

In vivo protein stabilization based on fragment complementation and a split GFP system

Protein stabilization was achieved through in vivo screening based on the thermodynamic linkage between protein folding and fragment complementation. The split GFP system was found suitable to derive protein variants with enhanced stability due to the correlation between effects of mutations on the stability of the intact chain and the effects of the same mutations on the affinity between fragments of the chain. PGB1 mutants with higher affinity between fragments 1 to 40 and 41 to 56 were obtained by in vivo screening of a library of the 1 to 40 fragments against wild-type 41 to 56 fragments. Colonies were ranked based on the intensity of green fluorescence emerging from assembly and folding of the fused GFP fragments. The DNA from the brightest fluorescent colonies was sequenced, and intact mutant PGB1s corresponding to the top three sequences were expressed, purified, and analyzed for stability toward thermal denaturation. The protein sequence derived from the top fluorescent colony was found to yield a 12 °C increase in the thermal denaturation midpoint and a free energy of stabilization of -8.7 kJ/mol at 25 °C. The stability rank order of the three mutant proteins follows the fluorescence rank order in the split GFP system. The variants are stabilized through increased hydrophobic effect, which raises the free energy of the unfolded more than the folded state; as well as substitutions, which lower the free energy of the folded more than the unfolded state; optimized van der Waals interactions; helix stabilization; improved hydrogen bonding network; and reduced electrostatic repulsion in the folded state.

Site-specific saturation mutagenesis on residues 132 and 450 of Candida rugosa LIP2 enhances catalytic efficiency and alters substrate specificity in various chain lengths of triglycerides and esters

The catalytic versatility of recombinant Candida rugosa LIP2 has been known to have potential applications in industry. In this study, site-specific saturation mutagenesis on residues L132 and G450 of recombinant LIP2 has been employed to investigate the impact of both residues on substrate specificity of LIP2. Point mutations on L132 and G450 were done separately using mutagenic degenerate primer sets containing 32 codons to generate two libraries of mutants in Pichia pastoris . Replacements of amino acid on these mutants were identified as L132A, L132I, G450S, and G450A. In lipase activity assay, L132A and L132I mutants showed a shift of preference from short- to medium-chain triglyceride, whereas G450S and G450A mutants retained preferences as compared to wild-type LIP2. Among mutants, G450A has the highest activity on tributyrin. However, hydrolysis of p-nitrophenyl (p-NP) esters with L132A, L132I, and G450S did not show differences of preferences over medium- to long-chain esters except in G450A, which prefers only medium-chain ester as compared to wild-type LIP2. All mutants showed an enhanced catalytic activity and higher optimal temperature and pH stability as compared to wild-type LIP2.

Increasing the stability of an enzyme toward hostile organic solvents by directed evolution based on iterative saturation mutagenesis using the B-FIT method

Mutants of the lipase from Bacillus subtilis, previously engineered for enhanced thermostability using directed evolution based on the B-FIT method, show significantly increased tolerance to hostile organic solvents.

Iterative saturation mutagenesis accelerates laboratory evolution of enzyme stereoselectivity: rigorous comparison with traditional methods

Efficacy in laboratory evolution of enzymes is currently a pressing issue, making comparative studies of different methods and strategies mandatory. Recent reports indicate that iterative saturation mutagenesis (ISM) provides a means to accelerate directed evolution of stereoselectivity and thermostability, but statistically meaningful comparisons with other methods have not been documented to date. In the present study, the efficacy of ISM has been rigorously tested by applying it to the previously most systematically studied enzyme in directed evolution, the lipase from Pseudomonas aeruginosa as a catalyst in the stereoselective hydrolytic kinetic resolution of a chiral ester. Upon screening only 10,000 transformants, unprecedented enantioselectivity was achieved (E = 594). ISM proves to be considerably more efficient than all previous systematic efforts utilizing error-prone polymerase chain reaction at different mutation rates, saturation mutagenesis at hot spots, and/or DNA shuffling, pronounced positive epistatic effects being the underlying reason.

Enzymes in food processing: a condensed overview on strategies for better biocatalysts

Food and feed is possibly the area where processing anchored in biological agents has the deepest roots. Despite this, process improvement or design and implementation of novel approaches has been consistently performed, and more so in recent years, where significant advances in enzyme engineering and biocatalyst design have fastened the pace of such developments. This paper aims to provide an updated and succinct overview on the applications of enzymes in the food sector, and of progresses made, namely, within the scope of tapping for more efficient biocatalysts, through screening, structural modification, and immobilization of enzymes. Targeted improvements aim at enzymes with enhanced thermal and operational stability, improved specific activity, modification of pH-activity profiles, and increased product specificity, among others. This has been mostly achieved through protein engineering and enzyme immobilization, along with improvements in screening. The latter has been considerably improved due to the implementation of high-throughput techniques, and due to developments in protein expression and microbial cell culture. Expanding screening to relatively unexplored environments (marine, temperature extreme environments) has also contributed to the identification and development of more efficient biocatalysts. Technological aspects are considered, but economic aspects are also briefly addressed.

Distance-dependent statistical potentials for discriminating thermophilic and mesophilic proteins

Identification of the characteristic structural patterns responsible for protein thermostability is theoretically important and practically useful but largely remains an open problem. These patterns may be revealed through comparative study on thermophilic and mesophilic proteins that have distinct thermostability. In this study, we constructed several distance-dependant potentials from thermophilic and mesophilic proteins. These potentials were then used to evaluate the structural difference between thermophilic and mesophilic proteins. We found that using the subtraction or division of the potentials derived from thermophilic and mesophilic proteins can dramatically increase the discriminatory ability. This approach revealed that the ability to distinct the subtle structural features responsible for protein thermostability may be effectively enhanced through rationally designed comparative study.

Miniaturization in biocatalysis

The use of biocatalysts for the production of both consumer goods and building blocks for chemical synthesis is consistently gaining relevance. A significant contribution for recent advances towards further implementation of enzymes and whole cells is related to the developments in miniature reactor technology and insights into flow behavior. Due to the high level of parallelization and reduced requirements of chemicals, intensive screening of biocatalysts and process variables has become more feasible and reproducibility of the bioconversion processes has been substantially improved. The present work aims to provide an overview of the applications of miniaturized reactors in bioconversion processes, considering multi-well plates and microfluidic devices, update information on the engineering characterization of the hardware used, and present perspective developments in this area of research.

Induced allostery in the directed evolution of an enantioselective Baeyer-Villiger monooxygenase

The molecular basis of allosteric effects, known to be caused by an effector docking to an enzyme at a site distal from the binding pocket, has been studied recently by applying directed evolution. Here, we utilize laboratory evolution in a different way, namely to induce allostery by introducing appropriate distal mutations that cause domain movements with concomitant reshaping of the binding pocket in the absence of an effector. To test this concept, the thermostable Baeyer-Villiger monooxygenase, phenylacetone monooxygenase (PAMO), was chosen as the enzyme to be employed in asymmetric Baeyer-Villiger reactions of substrates that are not accepted by the wild type. By using the known X-ray structure of PAMO, a decision was made regarding an appropriate site at which saturation mutagenesis is most likely to generate mutants capable of inducing allostery without any effector compound being present. After screening only 400 transformants, a double mutant was discovered that catalyzes the asymmetric oxidative kinetic resolution of a set of structurally different 2-substituted cyclohexanone derivatives as well as the desymmetrization of three different 4-substituted cyclohexanones, all with high enantioselectivity. Molecular dynamics (MD) simulations and covariance maps unveiled the origin of increased substrate scope as being due to allostery. Large domain movements occur that expose and reshape the binding pocket. This type of focused library production, aimed at inducing significant allosteric effects, is a viable alternative to traditional approaches to "designed" directed evolution that address the binding site directly.

A "molecular evolution" approach for isolation of intrinsically active (MEK-independent) MAP kinases

Mitogen-activated protein (MAP) kinases are a large family of enzymes composed of about four subfamilies, each containing several isoforms and splicing variants. Many MAP kinases are coexpressed in each eukaryotic cell and coactivated in response to various stimuli. It is, therefore, difficult to explore the specific downstream effects of each species of MAPK. Expression of an intrinsically active variant of a MAPK, while other MAPKs are not active, allows for tracking of a specific array of substrates, target genes, and biological/pathological effects corresponding to the expressed molecule. This chapter describes a method for obtaining such intrinsically active MAPKs. Because of the unique mode of MAPK activation, which is absolutely dependent on unconventional phosphorylation (on neighboring Thr + Tyr residues), a rational design of mutations that would render the kinase intrinsically active is currently unfeasible. Our method is based, therefore, on a "Molecular Evolution" approach that uses the power of yeast genetics and is unbiased toward the mutation sites. We describe in detail how to prepare a large population of randomly mutated molecules of the desired MAPK and how to screen this library in a yeast strain lacking the relevant MAPK kinase (MAPKK). The idea is to identify MAPK variants that are fulfilling all MAPK functions and allow growth of this strain - namely, MAPK molecules that function biologically in the complete absence of their upstream activator. We further describe the details of the "plasmid-loss" assay used for distinguishing between true positive and false positive clones. Finally, we report on a new yeast strain lacking four MAPKKs that could serve as a universal target for screening for active MAPK of all subfamilies.

Optimizing protein stability in vivo

Identifying mutations that stabilize proteins is challenging because most substitutions are destabilizing. In addition to being of immense practical utility, the ability to evolve protein stability in vivo may indicate how evolution has formed today's protein sequences. Here we describe a genetic selection that directly links the in vivo stability of proteins to antibiotic resistance. It allows the identification of stabilizing mutations within proteins. The large majority of mutants selected for improved antibiotic resistance are stabilized both thermodynamically and kinetically, indicating that similar principles govern stability in vivo and in vitro. The approach requires no prior structural or functional knowledge and allows selection for stability without a need to maintain function. Mutations that enhance thermodynamic stability of the protein Im7 map overwhelmingly to surface residues involved in binding to colicin E7, showing how the evolutionary pressures that drive Im7-E7 complex formation have compromised the stability of the isolated Im7 protein.

Enhancement of the thermostability and activity of mesophilic Clostridium cellulovorans EngD by in vitro DNA recombination with Clostridium thermocellum CelE

The thermal stability and catalytic activity of endoglucanase (EngD) from mesophilic Clostridium cellulovorans were improved by evolutionary molecular engineering. Thermostable mutants were isolated after staggered extension process (StEP) with celE from thermophilic Clostridium thermocellum performed to conduct family shuffling and overlay screening of the resultant mutant library. The relative activity of the best-evolved clone has been improved of about 2 times higher at 50 degrees C and showed a higher k(cat)/K(m) value than its engD parental clone. We determined that these variants had two amino acid substitutions (L157N, Q158E) and confirmed their effects by substituting these amino acids in the parental gene by site-directed mutagenesis. These substitutions resulted in an increase in hydrophilic or charged residues. Our results demonstrate that in vitro recombination is an effective approach to improve the thermostability and enzymatic activity of a mesophilic enzyme.