Rapid evolution of beta-glucuronidase specificity by saturation mutagenesis of an active site loop

Affinity maturation of antibody fragments: A review encompassing the development from random approaches to computational rational optimization

Routinely screened antibody fragments usually require further in vitro maturation to achieve the desired biophysical properties. Blind in vitro strategies can produce improved ligands by introducing random mutations into the original sequences and selecting the resulting clones under more and more stringent conditions. Rational approaches exploit an alternative perspective that aims first at identifying the specific residues potentially involved in the control of biophysical mechanisms, such as affinity or stability, and then to evaluate what mutations could improve those characteristics. The understanding of the antigen-antibody interactions is instrumental to develop this process the reliability of which, consequently, strongly depends on the quality and completeness of the structural information. Recently, methods based on deep learning approaches critically improved the speed and accuracy of model building and are promising tools for accelerating the docking step. Here, we review the features of the available bioinformatic instruments and analyze the reports illustrating the result obtained with their application to optimize antibody fragments, and nanobodies in particular. Finally, the emerging trends and open questions are summarized.

Nonspecific interaction and overlap concentration influence macromolecular crowding effect on glucose oxidase activity

Macromolecular crowding can change kinetics of enzyme catalysis. How interaction between enzymes and neighboring macromolecules contributes to the crowding effect on enzyme catalysis has not been quantitatively revealed. In this study, crowding effects of dextran and poly(ethylene glycol) (PEG) on glucose oxidase (GOx) are studied. Fluorescence resonance energy transfer experiments show the high transfer efficiency and stable interaction between the dextran and GOx. Further fluorescence quenching analysis also proves that the association of the dextran-GOx pair can become stronger than that of the PEG-GOx pair. Dextrans with concentrations above or below their chain overlap concentrations (c*) reduce Michaelis constants (K_m) of GOx catalysis by 90 % or 45 %, respectively, through volume exclusion mechanism, and in the meantime elevate the enzymatic efficiency (k_cat/K_m) by 8-fold or by 3-fold, respectively, which is more dramatic than that found in other enzymes before. Strong association between the enzyme and the dextran results in slow turnover rates (k_cat). Intermediate crowding with weak to moderate affinity to the enzyme below the c* can tune the k_cat higher than in the free state. Catalysis under crowded conditions is a joint effect of the enzyme-crowder nonspecific interaction, volume exclusion and overlap condition of the crowders.

Enhancing Enzyme Activity and Thermostability of Bacillus amyloliquefaciens Chitosanase BaCsn46A Through Saturation Mutagenesis at Ser196

Chitosanase plays an important role in chitooligosaccharides (COS) production. We found that the chitosanase (BaCsn46A) of Bacillus amyloliquefacien was a good candidate for chitosan hydrolysis of COS. In order to further improve the enzyme properties of BaCsn46A, the S196 located near the active center was found to be a critical site impacts on enzyme properties by sequence alignment analysis. Herein, saturation mutation was carried out to study role of 196 site on BaCsn46A catalytic function. Compared with WT, the specific enzyme activity of S196A increased by 118.79%, and the thermostability of S196A was much higher than WT. In addition, we found that the enzyme activity of S196P was 2.41% of that of WT, indicating that the type of amino acid in 196 site could significant affect the catalytic activity and thermostability of BaCsn46A. After molecular docking analysis we found that the increase in hydrogen bonds and decrease in unfavorable bonds interacting with the substrate were the main reason for the change of enzyme properties which is valuable for future studies on Bacillus species chitosanase.

Directed evolution of a cyclodipeptide synthase with new activities via label-free mass spectrometric screening

Directed evolution is a powerful approach to engineer enzymes via iterative creation and screening of variant libraries. However, assay development for high-throughput mutant screening remains challenging, particularly for new catalytic activities. Mass spectrometry (MS) analysis is label-free and well suited for untargeted discovery of new enzyme products but is traditionally limited by slow speed. Here we report an automated workflow for directed evolution of new enzymatic activities via high-throughput library creation and label-free MS screening. For a proof of concept, we chose to engineer a cyclodipeptide synthase (CDPS) that synthesizes diketopiperazine (DKP) compounds with therapeutic potential. In recombinant Escherichia coli, site-saturation mutagenesis (SSM) and error-prone PCR (epPCR) libraries expressing CDPS mutants were automatically created and cultivated on an integrated work cell. Culture supernatants were then robotically processed for matrix-assisted laser desorption/ionization time-of-flight (MALDI-ToF) MS analysis at a rate of 5 s per sample. The resulting mass spectral data were processed via custom computational algorithms, which performed a multivariant analysis of 108 theoretical mass-to-charge (m/z) values of 190 possible DKP molecules within a mass window of 115–373 Da. An F186L CDPS mutant was isolated to produce cyclo(l-Phe–l-Val), which is undetectable in the product profile of the wild-type enzyme. This robotic, label-free MS screening approach may be generally applicable to engineering other enzymes with new activities in high throughput.

A robotic workflow for directed evolution of new enzymatic activities via high-throughput library creation and label-free MS screening.

Trends in lipase engineering for enhanced biocatalysis

Lipases, also known as triacylglycerol hydrolases (E.C.No. 3.1.1.3), are considered as leading biocatalysts in the lipid modification business. With properties like ease of availability, capability to work in heterogeneous media, stability in organic solvents, property of catalyzing at the lipid-water interface and even in nonaqueous conditions, have made them a versatile choice for applications in the food, flavor, detergent, pharmaceutical, leather, textile, cosmetic, and paper industries [1]. The increasing alertness toward sustainable technologies, lesser waste generation and solvent usage and minimization of energy input has brought light toward the production and usage of recombinant/improved lipases. For example, Novozym 435, a broadly used recombinant lipase isolated from Candida antarctica, dominates the lipase industry and has even created a supplier bias in the market. This shows that there is a desperate need for novel, low-cost lipases with better properties. For this, mining of existing extremophilic genomes seems more rewarding. But considering the diversity of industrial requirements such as types of solvents used or carrier systems employed for enzyme immobilization, tailor-designed enzymes are an unrealized pressing priority. Therefore, protein engineering strategies in collaboration with the discovery of new lipases can serve as a vital tool to obtain tailor-made enzymes with specific characteristics.

Temperature-resistant and solvent-tolerant lipases as industrial biocatalysts: Biotechnological approaches and applications

The development of new biocatalytic systems to replace the chemical catalysts, with suitable characteristics in terms of efficiency, stability under high temperature reactions and in the presence of organic solvents, reusability, and eco-friendliness is considered a very important step to move towards the green processes. From this basis, the use of lipase as a catalyst is highly desired for many industrial applications because it offers the reactions in which could be used, stability in harsh conditions, reusability and a greener process. Therefore, the introduction of temperature-resistant and solvent-tolerant lipases have become essential and ideal for industrial applications. Temperature-resistant and solvent-tolerant lipases have been involved in many large-scale applications including biodiesel, detergent, food, pharmaceutical, organic synthesis, biosensing, pulp and paper, textile, animal feed, cosmetics, and leather industry. So, the present review provides a comprehensive overview of the industrial use of lipase. Moreover, special interest in biotechnological and biochemical techniques for enhancing temperature-resistance and solvent-tolerance of lipases to be suitable for the industrial uses.

Main Structural Targets for Engineering Lipase Substrate Specificity

Microbial lipases represent one of the most important groups of biotechnological biocatalysts. However, the high-level production of lipases requires an understanding of the molecular mechanisms of gene expression, folding, and secretion processes. Stable, selective, and productive lipase is essential for modern chemical industries, as most lipases cannot work in different process conditions. However, the screening and isolation of a new lipase with desired and specific properties would be time consuming, and costly, so researchers typically modify an available lipase with a certain potential for minimizing cost. Improving enzyme properties is associated with altering the enzymatic structure by changing one or several amino acids in the protein sequence. This review detailed the main sources, classification, structural properties, and mutagenic approaches, such as rational design (site direct mutagenesis, iterative saturation mutagenesis) and direct evolution (error prone PCR, DNA shuffling), for achieving modification goals. Here, both techniques were reviewed, with different results for lipase engineering, with a particular focus on improving or changing lipase specificity. Changing the amino acid sequences of the binding pocket or lid region of the lipase led to remarkable enzyme substrate specificity and enantioselectivity improvement. Site-directed mutagenesis is one of the appropriate methods to alter the enzyme sequence, as compared to random mutagenesis, such as error-prone PCR. This contribution has summarized and evaluated several experimental studies on modifying the substrate specificity of lipases.

Rapid and Error-Free Site-Directed Mutagenesis by a PCR-Free In Vitro CRISPR/Cas9-Mediated Mutagenic System

The quality and efficiency of any PCR-based mutagenesis technique may not be optimal due to, among other things, amino acid bias, which means that the development of efficient PCR-free methods is desirable. Here, we present a highly efficient in vitro CRISPR/Cas9-mediated mutagenic (ICM) system that allows rapid construction of designed mutants in a PCR-free manner. First, it involves plasmid digestion by utilizing a complex of Cas9 with specific single guide RNA (sgRNA) followed by degradation with T5 exonuclease to generate a 15 nt homologous region. Second, primers containing the desired mutations are annealed to form the double-stranded DNA fragments, which are then ligated into the linearized plasmid. In theory, neither the size of the target plasmid nor the unavailable restriction enzyme site poses any problems that may arise in traditional techniques. In this study, single and multiple site-directed mutagenesis were successfully performed even for a large size plasmid (up to 9.0 kb). Moreover, a PCR-free site-saturation mutagenesis library on single site and two adjacent sites of a green fluorescent protein was also generated with promising results. This demonstrates the great potential of the ICM system for creating high-quality mutant libraries in directed evolution as an alternative to PCR-based saturation mutagenesis, thus facilitating research on synthetic biology.

Controlled Aggregation and Increased Stability of β-Glucuronidase by Cellulose Binding Domain Fusion

Cellulose-binding domains (CBDs) are protein domains with cellulose-binding activity, and some act as leaders in the localization of cellulosomal scaffoldin proteins to the hydrophobic surface of crystalline cellulose. In this study, we found that a CBD fusion enhanced and improved soluble β-glucuronidase (GusA) enzyme properties through the formation of an artificially oligomeric state. First, a soluble CBD fused to the C-terminus of GusA (GusA-CBD) was obtained and characterized. Interestingly, the soluble GusA-CBD showed maximum activity at higher temperatures (65°C) and more acidic pH values (pH 6.0) than free GusA did (60°C and pH 7.5). Moreover, the GusA-CBD enzyme showed higher thermal and pH stabilities than the free GusA enzyme did. Additionally, GusA-CBD showed higher enzymatic activity in the presence of methanol than free GusA did. Evaluation of the protease accessibility of both enzymes revealed that GusA-CBD retained 100% of its activity after 1 h incubation in 0.5 mg/ml protease K, while free GusA completely lost its activity. Simple fusion of CBD as a single domain may be useful for tunable enzyme states to improve enzyme stability in industrial applications.

Ten plus one challenges in diseases of the lysosomal system

The advent of the first effective specific therapy for a lysosomal storage disease (LSDs), Gaucher disease type 1, by Roscoe O. Brady was foundational for development of additional treatments for this group of rare diseases. The past 26years, since the approval of enzyme therapy for Gaucher disease type 1, have witnessed a burgeoning understanding of LSDs at genetic, molecular, biochemical, cell biologic, and clinical levels. Simultaneously, this expansion of knowledge has exposed our incomplete understanding of the individual pathophysiologies of LSDs as well as difficult challenges for improvement in therapy and therapeutic outcomes for afflicted individuals. Here, 10 such challenges/problems representing major impediments, which need to be overcome, to move forward toward the goals of more effective and complete therapies for these devastating diseases.

In vitro Engineering of Novel Bioactivity in the Natural Enzymes

Enzymes catalyze various biochemical functions with high efficiency and specificity. In vitro design of the enzyme leads to novel bioactivity in this natural biomolecule that give answers of some vital questions like crucial residues in binding with substrate, molecular evolution, cofactor specificity etc. Enzyme engineering technology involves directed evolution, rational designing, semi-rational designing, and structure-based designing using chemical modifications. Similarly, combined computational and in vitro evolution approaches together help in artificial designing of novel bioactivity in the natural enzyme. DNA shuffling, error prone PCR and staggered extension process are used to artificially redesign active site of enzyme, which can alter its efficiency and specificity. Modifications of the enzyme can lead to the discovery of new path of molecular evolution, designing of efficient enzymes, locating active sites and crucial residues, shift in substrate, and cofactor specificity. The methods and thermodynamics of in vitro designing of the enzyme are also discussed. Similarly, engineered thermophilic and psychrophilic enzymes attain substrate specificity and activity of mesophilic enzymes that may also be beneficial for industry and therapeutics.

Cytochrome P450 genes in coronary artery diseases: Codon usage analysis reveals genomic GC adaptation

Establishing codon usage biases are imperative for understanding the etiology of coronary artery diseases (CAD) as well as the genetic factors associated with these diseases. The aim of this study was to evaluate the contribution of 18 responsible cytochrome P450 (CYP) genes for the risk of CAD. Effective number of codon (Nc) showed a negative correlation with both GC3 and synonymous codon usage order (SCUO) suggesting an antagonistic relationship between codon usage and Nc of genes. The dinucleotide analysis revealed that CG and TA dinucleotides have the lowest odds ratio in these genes. Principal component analysis showed that GC composition has a profound effect in separating the genes along the first major axis. Our findings revealed that mutational pressure and natural selection could possibly be the major factors responsible for codon bias in these genes. The study not only offers an insight into the mechanisms of genomic GC adaptation, but also illustrates the complexity of CYP genes in CAD.

Improved detection of Escherichia coli and coliform bacteria by multiplex PCR

BACKGROUND

The presence of coliform bacteria is routinely assessed to establish the microbiological safety of water supplies and raw or processed foods. Coliforms are a group of lactose-fermenting Enterobacteriaceae, which most likely acquired the lacZ gene by horizontal transfer and therefore constitute a polyphyletic group. Among this group of bacteria is Escherichia coli, the pathogen that is most frequently associated with foodborne disease outbreaks and is often identified by β-glucuronidase enzymatic activity or by the redundant detection of uidA by PCR. Because a significant fraction of essential E. coli genes are preserved throughout the bacterial kingdom, alternative oligonucleotide primers for specific E. coli detection are not easily identified.

RESULTS

In this manuscript, two strategies were used to design oligonucleotide primers with differing levels of specificity for the simultaneous detection of total coliforms and E. coli by multiplex PCR. A consensus sequence of lacZ and the orphan gene yaiO were chosen as targets for amplification, yielding 234 bp and 115 bp PCR products, respectively.

CONCLUSIONS

The assay designed in this work demonstrated superior detection ability when tested with lab collection and dairy isolated lactose-fermenting strains. While lacZ amplicons were found in a wide range of coliforms, yaiO amplification was highly specific for E. coli. Additionally, yaiO detection is non-redundant with enzymatic methods.

A rapid, efficient, and economical inverse polymerase chain reaction-based method for generating a site saturation mutant library

With the development of deep sequencing methodologies, it has become important to construct site saturation mutant (SSM) libraries in which every nucleotide/codon in a gene is individually randomized. We describe methodologies for the rapid, efficient, and economical construction of such libraries using inverse polymerase chain reaction (PCR). We show that if the degenerate codon is in the middle of the mutagenic primer, there is an inherent PCR bias due to the thermodynamic mismatch penalty, which decreases the proportion of unique mutants. Introducing a nucleotide bias in the primer can alleviate the problem. Alternatively, if the degenerate codon is placed at the 5' end, there is no PCR bias, which results in a higher proportion of unique mutants. This also facilitates detection of deletion mutants resulting from errors during primer synthesis. This method can be used to rapidly generate SSM libraries for any gene or nucleotide sequence, which can subsequently be screened and analyzed by deep sequencing.

Patch cloning method for multiple site-directed and saturation mutagenesis

BACKGROUND

Various DNA manipulation methods have been developed to prepare mutant genes for protein engineering. However, development of more efficient and convenient method is still demanded. Homologous DNA assembly methods, which do not depend on restriction enzymes, have been used as convenient tools for cloning and have been applied to site-directed mutagenesis recently. This study describes an optimized homologous DNA assembly method, termed as multiple patch cloning (MUPAC), for multiple site-directed and saturation mutagenesis.

RESULTS

To demonstrate MUPAC, we introduced five back mutations to a mutant green fluorescent protein (GFPuv) with five deleterious mutations at specific sites and transformed Escherichia coli (E. coli) with the plasmids obtained. We observed that the over 90% of resulting colonies possessed the plasmids containing the reverted GFPuv gene and exhibited fluorescence. We extended the test to introduce up to nine mutations in Moloney Murine Leukemia Virus reverse transcriptase (M-MLV RT) by assembling 11 DNA fragments using MUPAC. Analysis of the cloned plasmid by electrophoresis and DNA sequencing revealed that approximately 30% of colonies had the objective mutant M-MLV RT gene. Furthermore, we also utilized this method to prepare a library of mutant GFPuv genes containing saturation mutations at five specific sites, and we found that MUPAC successfully introduced NNK codons at all five sites, whereas other site remained intact.

CONCLUSIONS

MUPAC could efficiently introduce various mutations at multiple specific sites within a gene. Furthermore, it could facilitate the preparation of experimental gene materials important to molecular and synthetic biology research.

OptZyme: computational enzyme redesign using transition state analogues

OptZyme is a new computational procedure for designing improved enzymatic activity (i.e., kcat or kcat/KM) with a novel substrate. The key concept is to use transition state analogue compounds, which are known for many reactions, as proxies for the typically unknown transition state structures. Mutations that minimize the interaction energy of the enzyme with its transition state analogue, rather than with its substrate, are identified that lower the transition state formation energy barrier. Using Escherichia coli β-glucuronidase as a benchmark system, we confirm that KM correlates (R(2) = 0.960) with the computed interaction energy between the enzyme and the para-nitrophenyl- β, D-glucuronide substrate, kcat/KM correlates (R(2) = 0.864) with the interaction energy of the transition state analogue, 1,5-glucarolactone, and kcat correlates (R(2) = 0.854) with a weighted combination of interaction energies with the substrate and transition state analogue. OptZyme is subsequently used to identify mutants with improved KM, kcat, and kcat/KM for a new substrate, para-nitrophenyl- β, D-galactoside. Differences between the three libraries reveal structural differences that underpin improving KM, kcat, or kcat/KM. Mutants predicted to enhance the activity for para-nitrophenyl- β, D-galactoside directly or indirectly create hydrogen bonds with the altered sugar ring conformation or its substituents, namely H162S, L361G, W549R, and N550S.

Synonymous codon usage in forty staphylococcal phages identifies the factors controlling codon usage variation and the phages suitable for phage therapy

The immergence and dissemination of multidrug-resistant strains of Staphylococcus aureus in recent years have expedited the research on the discovery of novel anti-staphylococcal agents promptly. Bacteriophages have long been showing tremendous potentialities in curing the infections caused by various pathogenic bacteria including S. aureus. Thus far, only a few virulent bacteriophages, which do not carry any toxin-encoding gene but are capable of eradicating staphylococcal infections, were reported. Based on the codon usage analysis of sixteen S. aureus phages, previously three phages were suggested to be useful as the anti-staphylococcal agents. To search for additional S. aureus phages suitable for phage therapy, relative synonymous codon usage bias has been investigated in the protein-coding genes of forty new staphylococcal phages. All phages appeared to carry A and T ending codons. Several factors such as mutational pressure, translational selection and gene length seemed to be responsible for the codon usage variation in the phages. Codon usage indeed varied phage to phage. Of the phages, phages G1, Twort, 66 and Sap-2 may be extremely lytic in nature as majority of their genes possess high translational efficiency, indicating that these phages may be employed in curing staphylococcal infections.

A thermostable β-glucuronidase obtained by directed evolution as a reporter gene in transgenic plants

A β-glucuronidase variant, GUS-TR3337, that was obtained by directed evolution exhibited higher thermostability than the wild-type enzyme, GUS-WT. In this study, the utility of GUS-TR337 as an improved reporter was evaluated. The corresponding gus-tr3337 and gus-wt genes were independently cloned in a plant expression vector and introduced into Arabidopsis thaliana. With 4-MUG as a substrate, plants containing the gus-wt gene showed no detectable β-glucuronidase activity after exposure to 60°C for 10 min, while those hosting the gus-tr3337 gene retained 70% or 50% activity after exposure to 80°C for 10 min or 30 min, respectively. Similarly, in vivo β-glucuronidase activity could be demonstrated by using X-GLUC as a substrate in transgenic Arabidopsis plants hosting the gus-tr3337 gene that were exposed to 80°C for up to 30 min. Thus, the thermostability of GUS-TR3337 can be exploited to distinguish between endogenous and transgenic β-glucuronidase activity, which is a welcome improvement in its use as a reporter.

Industrial scale gene synthesis

The most recent developments in the area of deep DNA sequencing and downstream quantitative and functional analysis are rapidly adding a new dimension to understanding biochemical pathways and metabolic interdependencies. These increasing insights pave the way to designing new strategies that address public needs, including environmental applications and therapeutic inventions, or novel cell factories for sustainable and reconcilable energy or chemicals sources. Adding yet another level is building upon nonnaturally occurring networks and pathways. Recent developments in synthetic biology have created economic and reliable options for designing and synthesizing genes, operons, and eventually complete genomes. Meanwhile, high-throughput design and synthesis of extremely comprehensive DNA sequences have evolved into an enabling technology already indispensable in various life science sectors today. Here, we describe the industrial perspective of modern gene synthesis and its relationship with synthetic biology. Gene synthesis contributed significantly to the emergence of synthetic biology by not only providing the genetic material in high quality and quantity but also enabling its assembly, according to engineering design principles, in a standardized format. Synthetic biology on the other hand, added the need for assembling complex circuits and large complexes, thus fostering the development of appropriate methods and expanding the scope of applications. Synthetic biology has also stimulated interdisciplinary collaboration as well as integration of the broader public by addressing socioeconomic, philosophical, ethical, political, and legal opportunities and concerns. The demand-driven technological achievements of gene synthesis and the implemented processes are exemplified by an industrial setting of large-scale gene synthesis, describing production from order to delivery.

Single mutations outside the active site affect the substrate specificity in a β-glycosidase

A library of random mutants of the β-glycosidase Sfβgly was screened for mutations that affect its specificity for the substrate glycone (β-d-fucoside versus β-d-glucoside). Among mutations selected (T35A, R189G, Y345C, P348L, S358F, S378G, N400D, S424F, F460L, and R474H), eight occurred in the C-terminal half of Sfβgly and only two were at the active site (R189G and N400D). Tryptophan fluorescence spectra and thermal inactivation showed that the selected mutants and wild-type Sfβgly are similarly folded. Enzyme kinetics confirmed that these mutations resulted in broadening or narrowing of the preference for the substrate glycone. Structural modeling and interaction maps revealed contact pathways that connect the sites of the selected mutations through up to three interactions to the active site residues E399, W444, and E187, which are involved in substrate binding and catalysis. Interestingly, independently selected mutations (Y345C, P348L, and R189G; S424F and N400D) were placed on the same contact pathway. Moreover, (k(cat)/K(m) fucoside)/(k(cat)/K(m) glucoside) ratios showed that mutations at intermediate residues of the same contact pathway often had similar effects on substrate specificity. Finally mutations in the same contact pathway caused similar structural disturbance as evidenced by acrylamide quenching of the Sfβgly fluorescence. Based on these data, it is proposed that the effects of the selected mutations were propagated into the active site through groups of interacting residues (contact pathways) changing the Sfβgly substrate specificity.

Applying neutral drift to the directed molecular evolution of a β-glucuronidase into a β-galactosidase: Two different evolutionary pathways lead to the same variant

Background

Directed protein evolution has been used to modify protein activity and research has been carried out to enhance the production of high quality mutant libraries. Many theoretical approaches suggest that allowing a population to undergo neutral selection may be valuable in directed evolution experiments.

Findings

Here we report on an investigation into the value of neutral selection in a classical model system for directed evolution, the conversion of the E. coli β-glucuronidase to a β-galactosidase activity. We find that neutral selection, i.e. selection for retaining glucuronidase activity, can efficiently identify the majority of sites of mutation that have been identified as beneficial for galactosidase activity in previous experiments. Each variant demonstrating increased galactosidase activity identified by our neutral drift experiments contained a mutation at one of four sites, T509, S557, N566 or W529. All of these sites have previously been identified using direct selection for beta galactosidase activity.

Conclusions

Our results are consistent with others that show that a neutral selection approach can be effective in selecting improved variants. However, we interpret our results to show that neutral selection is, in this case, not a more efficient approach than conventional directed evolution approaches. However, the neutral approach is likely to be beneficial when the resulting library can be screened for a range of related activities. More detailed statistical studies to resolve the apparent differences between this system and others are likely to be a fruitful avenue for future research.

Directed in vitro evolution of reporter genes based on semi-rational design and high-throughput screening

Marker genes, such as gusA, lacZ, and gfp, have been applied comprehensively in biological studies. Directed in vitro evolution provides a powerful tool for modifying genes and for studying gene structure, expression, and function. Here, we describe a strategy for directed in vitro evolution of reporter genes based on semi-rational design and high-throughput screening. The protocol involves two processes of DNA shuffling and screening. The first DNA shuffling and screening process involves eight steps: (1) amplifying the target gene by PCR, (2) cutting the product into random fragments with DNase I, (3) purification of 50-100 bp fragments, (4) reassembly of the fragments in a primerless PCR, (5) amplification of the reassembled product by primer PCR, (6) cloning into expression vector, (7) transformation of E. coli by electroporation, and (8) screening the target mutants using a nitrocellulose filter. The second DNA shuffling and screening process also involves the same eight steps, except that degenerate oligonucleotide primers are based on the sequence of the selected mutant.

In silico saturation mutagenesis and docking screening for the analysis of protein-ligand interaction: the Endothelial Protein C Receptor case study

Background

The design of mutants in protein functional regions, such as the ligand binding sites, is a powerful approach to recognize the determinants of specific protein activities in cellular pathways. For an exhaustive analysis of selected positions of protein structure large scale mutagenesis techniques are often employed, with laborious and time consuming experimental set-up. 'In silico' mutagenesis and screening simulation represents a valid alternative to laboratory methods to drive the 'in vivo' testing toward more focused objectives.

Results

We present here a high performance computational procedure for large-scale mutant modelling and subsequent evaluation of the effect on ligand binding affinity. The mutagenesis was performed with a 'saturation' approach, where all 20 natural amino acids were tested in positions involved in ligand binding sites. Each modelled mutant was subjected to molecular docking simulation and stability evaluation. The simulated protein-ligand complexes were screened for their impairment of binding ability based on change of calculated Ki compared to the wild-type.

An example of application to the Endothelial Protein C Receptor residues involved in lipid binding is reported.

Conclusion

The computational pipeline presented in this work is a useful tool for the design of structurally stable mutants with altered affinity for ligand binding, considerably reducing the number of mutants to be experimentally tested. The saturation mutagenesis procedure does not require previous knowledge of functional role of the residues involved and allows extensive exploration of all possible substitutions and their pairwise combinations. Mutants are screened by docking simulation and stability evaluation followed by a rationally driven selection of those presenting the required characteristics. The method can be employed in molecular recognition studies and as a preliminary approach to select models for experimental testing.

Mutagenesis of the phosphate-binding pocket of KDPG aldolase enhances selectivity for hydrophobic substrates

Narrow substrate specificities often limit the use of enzymes in biocatalysis. To further the development of Escherichia coli 2-keto-3-deoxy-6-phosphogluconate (KDPG) aldolase as a biocatalyst, the molecular determinants of substrate specificity were probed by mutagenesis. Our data demonstrate that S184 is located in the substrate-binding pocket and interacts with the phosphate moiety of KDPG, providing biochemical support for the binding model proposed on the basis of crystallographic data. An analysis of the substrate selectivity of the mutant enzymes indicates that alterations to the phosphate-binding site of KDPG aldolase changes the substrate selectivity. We report mutations that enhance catalysis of aldol cleavage of substrates lacking a phosphate moiety and demonstrate that electrophile reactivity correlates with the hydrophobicity of the substituted side chain. These mutations improve the selectivity for unnatural substrates as compared to KDPG by up to 2000-fold. Furthermore, the S184L KDPG aldolase mutant improves the catalytic efficiency for the synthesis of a precursor for nikkomycin by 40-fold, making it a useful biocatalyst for the preparation of fine chemicals.

Better library design: data-driven protein engineering

Data-driven protein engineering is increasingly used as an alternative to rational design and combinatorial engineering because it uses available knowledge to limit library size, while still allowing for the identification of unpredictable substitutions that lead to large effects. Recent advances in computational modeling and bioinformatics, as well as an increasing databank of experiments on functional variants, have led to new strategies to choose particular amino acid residues to vary in order to increase the chances of obtaining a variant protein with the desired property. Strategies for limiting diversity at each position, design of small sub-libraries, and the performance of scouting experiments, have also been developed or even automated, further reducing the library size.

Directed evolution of a beta-galactosidase from Pyrococcus woesei resulting in increased thermostable beta-glucuronidase activity

We performed directed evolution on a chemically synthesized 1,533-bp recombinant beta-galactosidase gene from Pyrococcus woesei. More than 200,000 variant colonies in each round of directed evolution were screened using the pYPX251 vector and host strain Rosetta-Blue (DE3). One shifted beta-galactosidase to beta-glucuronidase mutant, named YG6762, was obtained after four rounds of directed evolution and screening. This mutant had eight mutated amino acid residues. T29A, V213I, L217M, N277H, I387V, R491C, and N496D were key mutations for high beta-glucuronidase activity, while E414D was not essential because the mutation did not lead to a change in beta-glucuronidase activity. The amino acid site 277 was the most essential because mutating H back to N resulted in a 50% decrease in beta-glucuronidase activity at 37 degrees C. We also demonstrated that amino acid 277 was the most essential site, as the mutation from N to H resulted in a 1.5-fold increase in beta-glucuronidase activity at 37 degrees C. Although most single amino acid changes lead to less than a 20% increase in beta-glucuronidase activity, the YG6762 variant, which was mutated at all eight amino acid sites, had a beta-glucuronidase activity that was about five and seven times greater than the wild-type enzyme at 37 and 25 degrees C, respectively.

Antibody-induced oligomerization and activation of an engineered reporter enzyme

Our objective is to produce a protein biosensor (or molecular switch) that is specifically activated in solution by a monoclonal antibody. Many effector-dependent enzymes have evolved in nature, but the introduction of a novel regulatory mechanism into a normally unregulated enzyme poses a difficult design problem. We used site-saturation mutagenesis and screening to generate effector-activated variants of the reporter enzyme beta-glucuronidase (GUS). The specific activity of the purified epitope-tagged GUS variant was increased by up to approximately 500-fold by the addition of an equimolar concentration of a monoclonal antibody. This molecular switch is modular in design, so it can easily be re-engineered for the detection of other peptide-specific antibodies. Such antibody-activated reporters could someday enable point-of-care serological assays for the rapid detection of infectious diseases.

Concurrent mutations in six amino acids in β-glucuronidase improve its thermostability

To achieve a thermostable beta-glucuronidase (GUS) and identify key mutation sites, we applied in vitro directed evolution strategy through DNA shuffling and obtained a highly thermostable mutant GUS gene, gus-tr, after four rounds of DNA shuffling and screening. This variant had mutations in 15 nucleic acid sites, resulting in changes in 12 amino acids (AAs). Using gus-tr as the template, we further performed site-directed mutagenesis to reverse the individual mutation to the wild-type protein. We found that six sites (Q493R, T509A, M532T, N550S, G559S and N566S) present in GUS-TR3337, were the key AAs needed to confer its high thermostability. Of these, Q493R and T509A were not reported previously as important residues for thermostability of GUS. Furthermore, all of these six mutations must be present concurrently to confer the high thermostability. We expressed the gus-tr3337 gene and purified the GUS-TR3337 protein that contained the six AA mutations. Compared with the wild-type protein which lost its activity completely after 10 min at 70 degrees C, the mutant GUS-TR3337 protein retained 75% of its activity when heated at 80 degrees C for 10 min. The GUS-TR3337 exhibited high activity even heated at 100 degrees C for 30 min on nitrocellulose filter. The comparison of molecular models of the mutated and wild-type enzyme revealed the relation of protein function and these structural modifications.

Directed Evolution of Beta-galactosidase from Escherichia coli into Beta-glucuronidase

In vitro directed evolution through DNA shuffling is a powerful molecular tool for creation of new biological phenotypes. E. coli beta-galactosidase and beta-glucuronidase are widely used, and their biological function, catalytic mechanism, and molecular structures are well characterized. We applied an in vitro directed evolution strategy through DNA shuffling and obtained five mutants named YG6764, YG6768, YG6769, YG6770 and YG6771 after two rounds of DNA shuffling and screening, which exhibited more beta-glucuronidase activity than wild-type beta-galactosidase. These variants had mutations at fourteen nucleic acid sites, resulting in changes in ten amino acids: S193N, T266A, Q267R, V411A, D448G, G466A, L527I, M543I, Q626R and Q951R. We expressed and purified those mutant proteins. Compared to the wild-type protein, five mutant proteins exhibited high beta-glucuronidase activity. The comparison of molecular models of the mutated and wildtype enzymes revealed the relationship between protein function and structural modification.

Semi-rational approaches to engineering enzyme activity: combining the benefits of directed evolution and rational design

Many research groups successfully rely on whole-gene random mutagenesis and recombination approaches for the directed evolution of enzymes. Recent advances in enzyme engineering have used a combination of these random methods of directed evolution with elements of rational enzyme modification to successfully by-pass certain limitations of both directed evolution and rational design. Semi-rational approaches that target multiple, specific residues to mutate on the basis of prior structural or functional knowledge create 'smart' libraries that are more likely to yield positive results. Efficient sampling of mutations likely to affect enzyme function has been conducted both experimentally and, on a much greater scale, computationally, with remarkable improvements in substrate selectivity and specificity and in the de novo design of enzyme activities within scaffolds of known structure.

Iterative saturation mutagenesis (ISM) for rapid directed evolution of functional enzymes

Iterative saturation mutagenesis (ISM) is a new and efficient method for the directed evolution of functional enzymes. It reduces the necessary molecular biological work and the screening effort drastically. It is based on a Cartesian view of the protein structure, performing iterative cycles of saturation mutagenesis at rationally chosen sites in an enzyme, a given site being composed of one, two or three amino acid positions. The basis for choosing these sites depends on the nature of the catalytic property to be improved, e.g., enantioselectivity, substrate acceptance or thermostability. In the case of thermostability, sites showing highest B-factors (available from X-ray data) are chosen. The pronounced increase in thermostability of the lipase from Bacillus subtilis (Lip A) as a result of applying ISM is illustrated here.

Rapid expression of functional genomic libraries

Genomic-scale analysis of protein function is currently limited by the ability to rapidly express the enormous diversity of protein targets in their active form. We describe a method to construct transcriptionally active expression templates (ETs) in parallel using a single PCR step wherein the overlap-extension reaction for addition of transcription regulatory elements is separated from the amplification of the full-length product by using a GC-rich single primer. Over 90% of 55 diverse genomic targets were extended with T7 regulatory elements to form ETs in high yield and purity. The unpurified ETs directed protein expression using a cell-free protein synthesis (CFPS) system supplemented with cofactors and metal ions to activate a variety of enzymes. Higher activities were obtained in the modified CFPS reactions compared to standard reaction conditions. Protein purification was avoided because the expressed enzyme activity was significantly greater than the background activity associated with the cell extract. These improvements in the parallel synthesis of linear ETs combined with enhanced in vitro enzyme activation help to make CFPS systems more attractive platforms for high-throughput evaluation of protein function.

High efficiency and throughput system in directed evolution in vitro of reporter gene

In vitro directed evolution, especially with DNA shuffling, is a powerful means in biological studies of protein structure and function, and consequently for industrial applications. Escherichia coli beta-glucuronidase (gusA) gene, a versatile and efficient reporter gene, was the model for studying in vitro directed evolution because of its stability, easy analysis of the enzyme properties and conveniently visible phenotype. We developed a high efficiency, throughput system for in vitro directed evolution using gusA reporter gene as the model. The system consisted mainly of three aspects: a prokaryotic expression vector pYPX251, an easy method for obtaining the mutated gene from DNA shuffling and a suitable selected strategy. The vector pYPX251 carried the moderately strong aacC1 gene promoter and T1T2 transcription terminator that allowed expression in E. coli. Over 10,000 individuals could be selected individually in a 9 cm Petri dish after colonies were absorbed on a nitrocellulose filter. A library, which contained 100,000 individuals was screened by incubating ten filter papers with X-Glu. The polymerase chain reaction products of the gusA gene, the fragments of 50-100 bp, with high mutation rates were purified using a dialysis bag from 10% PAGE after electrophoresis. The possibility of obtaining desirable mutations was increased dramatically as the size of the library expanded. A GUS variant, named GUS-TR, was obtained through this system, which is significantly more resistant to high temperature than the wild type enzyme. GUS-TR maintained its high activity even when the nitrocellulose filter containing the variant colony was heated at 100 degrees C for 30 min.

Diversification and specialization of HIV protease function during in vitro evolution

Our goal is to understand how enzymes adapt to utilize novel substrates. We and others have shown that directed evolution tends to generate enzyme variants with broadened substrate specificity. Broad-specificity enzymes are generally deleterious to living cells, so this observed trend might be an artifact of the most commonly employed high throughput screens. Here, we demonstrate a more natural and effective screening strategy for directed evolution. The gene encoding model enzyme HIV protease was randomly mutated, and the resulting library was expressed in Escherichia coli cells to eliminate cytotoxic broad-specificity variants. The surviving variants were screened for clones with activity against a reporter enzyme. The wild-type human immunodeficiency virus type I protease (HIV PR) is cytotoxic and exhibits no detectable activity in reactions with beta-galactosidase (BGAL). In contrast, the selected variants were nontoxic and exhibited greater activity and specificity against BGAL than did the wild-type HIV PR in reactions with any substrate. A single round of whole gene random mutagenesis and conventional high-throughput screening does not usually effect complete inversions of substrate specificity. This suggests that a combination of positive and purifying selection engenders more rapid adaptation than positive selection alone.

Directed evolution of RuBisCO hypermorphs through genetic selection in engineered E.coli

The Calvin Cycle is the primary conduit for the fixation of carbon dioxide into the biosphere; ribulose 1,5-bisphosphate carboxylase/oxygenase (RuBisCO) catalyzes the rate-limiting fixation step. Our goal is to direct the evolution of RuBisCO variants with improved kinetic and biophysical properties. The Calvin Cycle was partially reconstructed in Escherichia coli; the engineered strain requires the Synechococcus PCC6301 RuBisCO for growth in minimal media supplemented with a pentose. We randomly mutated the gene encoding the large subunit of RuBisCO (rbcL), co-expressed the resulting library with the small subunit (rbcS) and the Synechococcus PCC7492 phosphoribulokinase (prkA), and selected hypermorphic variants. The RuBisCO variants that evolved during three rounds of random mutagenesis and selection were over-expressed, and exhibited 5-fold improvement in specific activity relative to the wild-type enzyme. These results demonstrate a new strategy for the artificial selection of RuBisCO and other non-native metabolic enzymes.