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

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

Improvement of an antibody neutralizing the anthrax toxin by simultaneous mutagenesis of its six hypervariable loops

The enhancement of antibody affinity by mutagenesis targeting only complementarity determining regions has the advantage of respecting the framework regions, which are important for tolerance if clinical use is envisaged. Here, starting from a Fab (antigen-binding fragment; 35PA₈₃) capable of neutralizing the lethal toxin of anthrax and having an affinity of 3.4 nM for its antigen, a phage-displayed library of variants where all six complementarity determining regions (73 positions) were targeted for mutagenesis was built. This library contained 5 × 10⁸ variants, and each variant carried four mutations on average. The library was first panned with adsorbed antigen and washes of increasing stringency. It was then screened in parallel with either small concentrations of soluble biotinylated antigen or adsorbed antigen and long elution times in the presence of soluble antigen. The stringencies of both selections were pushed as far as possible. Compared with 35PA₈₃, the best selected clone had an affinity enhanced 19-fold, to 180 pM, and its 50% inhibitory concentration was decreased by 40%. The results of the two selection methods were compared, and the generality of these methods was considered.     

In-vitro selection of highly stabilized protein variants with optimized surface

Thermostable proteins are of prime importance in protein science, but it has remained difficult to develop general strategies for stabilizing a protein. Site-directed mutagenesis based on comparisons with thermophilic homologs is rarely successful because the sequence differences are too numerous and dominated by neutral mutations. Here we used a method of directed evolution to increase the stability of a mesophilic protein, the cold shock protein Bs-CspB from Bacillus subtilis. It differs from its thermophilic counterpart Bc-Csp from Bacillus caldolyticus at 12 surface-exposed positions. To elucidate the stabilizing potential of exposed amino acid residues, six of these variant positions were randomized by saturation mutagenesis, the corresponding library of sequences was inserted into the gene-3-protein of the filamentous phage fd, and stabilized variants were selected by the Proside technique. Proside links the increased protease resistance of stabilized protein variants with the infectivity of the phage. Many strongly stabilized variants of Bs-CspB were identified in two selections, one in the presence of a denaturant and the other at elevated temperature. Several of them are significantly more stable than the naturally thermostable homolog Bc-Csp, and the best variant reaches Tm-Csp (the homolog from the hyperthermophile Thermotoga maritima) in stability. Remarkably, this variant differs from Tm-Csp at five and from Bc-Csp at all six randomized positions. This indicates that proteins can be strongly stabilized by many different sets of surface mutations, and Proside selects them efficiently from large libraries. The course of the selection could be directed by the conditions. In an ionic denaturant non-polar surface interactions were optimized, whereas at elevated temperature variants with improved electrostatics were selected, pointing to two different strategies for stabilization at protein surfaces.     

The Association of the Vanin-1 N131S Variant with Blood Pressure Is Mediated by Endoplasmic Reticulum-Associated Degradation and Loss of Function

High blood pressure (BP) is the most common cardiovascular risk factor worldwide and a major contributor to heart disease and stroke. We previously discovered a BP-associated missense SNP (single nucleotide polymorphism)–rs2272996–in the gene encoding vanin-1, a glycosylphosphatidylinositol (GPI)-anchored membrane pantetheinase. In the present study, we first replicated the association of rs2272996 and BP traits with a total sample size of nearly 30,000 individuals from the Continental Origins and Genetic Epidemiology Network (COGENT) of African Americans (P = 0.01). This association was further validated using patient plasma samples; we observed that the N131S mutation is associated with significantly lower plasma vanin-1 protein levels. We observed that the N131S vanin-1 is subjected to rapid endoplasmic reticulum-associated degradation (ERAD) as the underlying mechanism for its reduction. Using HEK293 cells stably expressing vanin-1 variants, we showed that N131S vanin-1 was degraded significantly faster than wild type (WT) vanin-1. Consequently, there were only minimal quantities of variant vanin-1 present on the plasma membrane and greatly reduced pantetheinase activity. Application of MG-132, a proteasome inhibitor, resulted in accumulation of ubiquitinated variant protein. A further experiment demonstrated that atenolol and diltiazem, two current drugs for treating hypertension, reduce the vanin-1 protein level. Our study provides strong biological evidence for the association of the identified SNP with BP and suggests that vanin-1 misfolding and degradation are the underlying molecular mechanism.     

A <Emphasis Type="Italic">Paenibacillus</Emphasis> sp. dextranase mutant pool with improved thermostability and activity

Random mutagenesis was used to create a library of chimeric dextranase (dex1) genes. A plate-screening protocol was developed with improved thermostability as a selection criterion. The mutant library was screened for active dextranase variants by observing clearing zones on dextran-blue agar plates at 50°C after exposure to 68°C for 2 h, a temperature regime at which wild-type activity was abolished. A number of potentially improved variants were identified by this strategy, five of which were further characterised. DNA sequencing revealed ten nucleotide substitutions, ranging from one to four per variant. Thermal inactivation studies showed reduced (2.9-fold) thermostability for one variant and similar thermostability for a second variant, but confirmed improved thermostability for three mutants with 2.3- (28.9 min) to 6.9-fold (86.6 min) increases in half-lives at 62°C compared to that of the wild-type enzyme (12.6 min). Using a 10-min assay, apparent temperature optima of the variants were similar to that of the wild type (T _opt 60°C). However, one of these variants had increased enzyme activity. Therefore, the first-generation dextranase mutant pool obtained in this study has sufficient molecular diversity for further improvements in both thermostability and activity through recombination (gene shuffling).     

Affinity maturation of antiHER2 monoclonal antibody MIL5 using an epitope-specific synthetic phage library by computational design

Increased affinities mainly equal to improved biological efficacy in many cases. By now, display methods including phage library are widely exploited to obtain higher affinity antibodies. Traditional library methods mainly focus on complementary determining region mutagenesis, in which extensive screening of variant combinations as well as large library capacity is required to find higher affinity clones. In this study, based on the modeling 3D complex structure of antigen (HER2)–antibody (MIL5) complex, the key residues of contact surface were predicted and identified to guide the synthetic phage library design. Then, epitope-specific site-directed mutagenesis phage library comprised of MIL5_scFv mutants was constructed, from which a higher affinity single chain antibody (M5scFv_ph) was screened out. Following experimental results showed that the novel antibody M5scFv_ph retained superimposed epitope to the parent antibody MIL5_scFv, and possessed similar tumor growth inhibitory activity in vivo on ovarian carcinoma xenografts.     

The Fitness Effects of Random Mutations in Single-Stranded DNA and RNA Bacteriophages

Mutational fitness effects can be measured with relatively high accuracy in viruses due to their small genome size, which facilitates full-length sequencing and genetic manipulation. Previous work has shown that animal and plant RNA viruses are very sensitive to mutation. Here, we characterize mutational fitness effects in single-stranded (ss) DNA and ssRNA bacterial viruses. First, we performed a mutation-accumulation experiment in which we subjected three ssDNA (ΦX174, G4, F1) and three ssRNA phages (Qβ, MS2, and SP) to plaque-to-plaque transfers and chemical mutagenesis. Genome sequencing and growth assays indicated that the average fitness effect of the accumulated mutations was similar in the two groups. Second, we used site-directed mutagenesis to obtain 45 clones of ΦX174 and 42 clones of Qβ carrying random single-nucleotide substitutions and assayed them for fitness. In ΦX174, 20% of such mutations were lethal, whereas viable ones reduced fitness by 13% on average. In Qβ, these figures were 29% and 10%, respectively. It seems therefore that high mutational sensitivity is a general property of viruses with small genomes, including those infecting animals, plants, and bacteria. Mutational fitness effects are important for understanding processes of fitness decline, but also of neutral evolution and adaptation. As such, these findings can contribute to explain the evolution of ssDNA and ssRNA viruses.     

Protein Evolution via Amino Acid and Codon Elimination

Background

Global residue-specific amino acid mutagenesis can provide important biological insight and generate proteins with altered properties, but at the risk of protein misfolding. Further, targeted libraries are usually restricted to a handful of amino acids because there is an exponential correlation between the number of residues randomized and the size of the resulting ensemble. Using GFP as the model protein, we present a strategy, termed protein evolution via amino acid and codon elimination, through which simplified, native-like polypeptides encoded by a reduced genetic code were obtained via screening of reduced-size ensembles.

Methodology/Principal Findings

The strategy involves combining a sequential mutagenesis scheme to reduce library size with structurally stabilizing mutations, chaperone complementation, and reduced temperature of gene expression. In six steps, we eliminated a common buried residue, Phe, from the green fluorescent protein (GFP), while retaining activity. A GFP variant containing 11 Phe residues was used as starting scaffold to generate 10 separate variants in which each Phe was replaced individually (in one construct two adjacent Phe residues were changed simultaneously), while retaining varying levels of activity. Combination of these substitutions to generate a Phe-free variant of GFP abolished fluorescence. Combinatorial re-introduction of five Phe residues, based on the activities of the respective single amino acid replacements, was sufficient to restore GFP activity. Successive rounds of mutagenesis generated active GFP variants containing, three, two, and zero Phe residues. These GFPs all displayed progenitor-like fluorescence spectra, temperature-sensitive folding, a reduced structural stability and, for the least stable variants, a reduced steady state abundance.

Conclusions/Significance

The results provide strategies for the design of novel GFP reporters. The described approach offers a means to enable engineering of active proteins that lack certain amino acids, a key step towards expanding the functional repertoire of uniquely labeled proteins in synthetic biology.     

Unsupervised Inference of Protein Fitness Landscape from Deep Mutational Scan

The recent technological advances underlying the screening of large combinatorial libraries in high-throughput mutational scans deepen our understanding of adaptive protein evolution and boost its applications in protein design. Nevertheless, the large number of possible genotypes requires suitable computational methods for data analysis, the prediction of mutational effects, and the generation of optimized sequences. We describe a computational method that, trained on sequencing samples from multiple rounds of a screening experiment, provides a model of the genotype–fitness relationship. We tested the method on five large-scale mutational scans, yielding accurate predictions of the mutational effects on fitness. The inferred fitness landscape is robust to experimental and sampling noise and exhibits high generalization power in terms of broader sequence space exploration and higher fitness variant predictions. We investigate the role of epistasis and show that the inferred model provides structural information about the 3D contacts in the molecular fold.     

Deleterious- and Disease-Allele Prevalence in Healthy Individuals: Insights from Current Predictions,Mutation Databases,and Population-Scale Resequencing

We have assessed the numbers of potentially deleterious variants in the genomes of apparently healthy humans by using (1) low-coverage whole-genome sequence data from 179 individuals in the 1000 Genomes Pilot Project and (2) current predictions and databases of deleterious variants. Each individual carried 281–515 missense substitutions, 40–85 of which were homozygous, predicted to be highly damaging. They also carried 40–110 variants classified by the Human Gene Mutation Database (HGMD) as disease-causing mutations (DMs), 3–24 variants in the homozygous state, and many polymorphisms putatively associated with disease. Whereas many of these DMs are likely to represent disease-allele-annotation errors, between 0 and 8 DMs (0–1 homozygous) per individual are predicted to be highly damaging, and some of them provide information of medical relevance. These analyses emphasize the need for improved annotation of disease alleles both in mutation databases and in the primary literature; some HGMD mutation data have been recategorized on the basis of the present findings, an iterative process that is both necessary and ongoing. Our estimates of deleterious-allele numbers are likely to be subject to both overcounting and undercounting. However, our current best mean estimates of ∼400 damaging variants and ∼2 bona fide disease mutations per individual are likely to increase rather than decrease as sequencing studies ascertain rare variants more effectively and as additional disease alleles are discovered.     

Repeated Random Mutagenesis of alpha-Amylase from Bacillus licheniformis for Improved pH Performance

The alpha-amylases activity was improved by random mutagenesis and screening. A region comprising residues from the position 34-281 was randomly mutated in B. licheniformis alpha-amylase (AmyL), and the library with mutations ranging from low, medium, and high frequencies was generated. The library was screened using an effective liquid-phase screening method to isolate mutants with an altered pH profile. The sequencing of improved variants indicated 2-5 amino acid changes. Among them, mutant TP8H5 showed an altered pH profile as compared with that of wild type. The sequencing of variant TP8H5 indicated 2 amino acid changes, Ile157Ser and Trp193Arg, which were located in the solvent accessible flexible loop region in domain B.     

Deletional Protein Engineering Based on Stable Fold

Diversification of protein sequence-structure space is a major concern in protein engineering. Deletion mutagenesis can generate a protein sequence-structure space different from substitution mutagenesis mediated space, but it has not been widely used in protein engineering compared to substitution mutagenesis, because it causes a relatively huge range of structural perturbations of target proteins which often inactivates the proteins. In this study, we demonstrate that, using green fluorescent protein (GFP) as a model system, the drawback of the deletional protein engineering can be overcome by employing the protein structure with high stability. The systematic dissection of N-terminal, C-terminal and internal sequences of GFPs with two different stabilities showed that GFP with high stability (s-GFP), was more tolerant to the elimination of amino acids compared to a GFP with normal stability (n-GFP). The deletion studies of s-GFP enabled us to achieve three interesting variants viz. s-DL4, s-N14, and s-C225, which could not been obtained from n-GFP. The deletion of 191–196 loop sequences led to the variant s-DL4 that was expressed predominantly as insoluble form but mostly active. The s-N14 and s-C225 are the variants without the amino acid residues involving secondary structures around N- and C-terminals of GFP fold respectively, exhibiting comparable biophysical properties of the n-GFP. Structural analysis of the variants through computational modeling study gave a few structural insights that can explain the spectral properties of the variants. Our study suggests that the protein sequence-structure space of deletion mutants can be more efficiently explored by employing the protein structure with higher stability.     

Enhancing the Thermal Tolerance and Gastric Performance of a Microbial Phytase for Use as a Phosphate-Mobilizing Monogastric-Feed Supplement

The inclusion of phytase in monogastric animal feed has the benefit of hydrolyzing indigestible plant phytate (myo-inositol 1,2,3,4,5,6-hexakis dihydrogen phosphate) to provide poultry and swine with dietary phosphorus. An ideal phytase supplement should have a high temperature tolerance, allowing it to survive the feed pelleting process, a high specific activity at low pHs, and adequate gastric performance. For this study, the performance of a bacterial phytase was optimized by the use of gene site saturation mutagenesis technology. Beginning with the appA gene from Escherichia coli, a library of clones incorporating all 19 possible amino acid changes and 32 possible codon variations in 431 residues of the sequence was generated and screened for mutants exhibiting improved thermal tolerance. Fourteen single site variants were discovered that retained as much as 10 times the residual activity of the wild-type enzyme after a heated incubation regimen. The addition of eight individual mutations into a single construct (Phy9X) resulted in a protein of maximal fitness, i.e., a highly active phytase with no loss of activity after heating at 62°C for 1 h and 27% of its initial activity after 10 min at 85°C, which was a significant improvement over the appA parental phytase. Phy9X also showed a 3.5-fold enhancement in gastric stability.     

High-throughput assay and engineering of self-cleaving ribozymes by sequencing

Self-cleaving ribozymes are found in all domains of life and are believed to play important roles in biology. Additionally, self-cleaving ribozymes have been the subject of extensive engineering efforts for applications in synthetic biology. These studies often involve laborious assays of multiple individual variants that are either designed rationally or discovered through selection or screening. However, these assays provide only a limited view of the large sequence space relevant to the ribozyme function. Here, we report a strategy that allows quantitative characterization of greater than 1000 ribozyme variants in a single experiment. We generated a library of predefined ribozyme variants that were converted to DNA and analyzed by high-throughput sequencing. By counting the number of cleaved and uncleaved reads of every variant in the library, we obtained a complete activity profile of the ribozyme pool which was used to both analyze and engineer allosteric ribozymes.     

Structural dynamics shape the fitness window of alanine:glyoxylate aminotransferase

The conformational landscape of a protein is constantly expanded by genetic variations that have a minimal impact on the function(s) while causing subtle effects on protein structure. The wider the conformational space sampled by these variants, the higher the probabilities to adapt to changes in environmental conditions. However, the probability that a single mutation may result in a pathogenic phenotype also increases. Here we present a paradigmatic example of how protein evolution balances structural stability and dynamics to maximize protein adaptability and preserve protein fitness. We took advantage of known genetic variations of human alanine:glyoxylate aminotransferase (AGT1), which is present as a common major allelic form (AGT‐Ma) and a minor polymorphic form (AGT‐Mi) expressed in 20% of Caucasian population. By integrating crystallographic studies and molecular dynamics simulations, we show that AGT‐Ma is endowed with structurally unstable (frustrated) regions, which become disordered in AGT‐Mi. An in‐depth biochemical characterization of variants from an anticonsensus library, encompassing the frustrated regions, correlates this plasticity to a fitness window defined by AGT‐Ma and AGT‐Mi. Finally, co‐immunoprecipitation analysis suggests that structural frustration in AGT1 could favor additional functions related to protein–protein interactions. These results expand our understanding of protein structural evolution by establishing that naturally occurring genetic variations tip the balance between stability and frustration to maximize the ensemble of conformations falling within a well‐defined fitness window, thus expanding the adaptability potential of the protein.     

Engineered ribosomal RNA operon copy-number variants of E. coli reveal the evolutionary trade-offs shaping rRNA operon number

Ribosomal RNA (rrn) operons, characteristically present in several copies in bacterial genomes (7 in E. coli), play a central role in cellular physiology. We investigated the factors determining the optimal number of rrn operons in E. coli by constructing isogenic variants with 5–10 operons. We found that the total RNA and protein content, as well as the size of the cells reflected the number of rrn operons. While growth parameters showed only minor differences, competition experiments revealed a clear pattern: 7–8 copies were optimal under conditions of fluctuating, occasionally rich nutrient influx and lower numbers were favored in stable, nutrient-limited environments. We found that the advantages of quick adjustment to nutrient availability, rapid growth and economic regulation of ribosome number all contribute to the selection of the optimal rrn operon number. Our results suggest that the wt rrn operon number of E. coli reflects the natural, ‘feast and famine’ life-style of the bacterium, however, different copy numbers might be beneficial under different environmental conditions. Understanding the impact of the copy number of rrn operons on the fitness of the cell is an important step towards the creation of functional and robust genomes, the ultimate goal of synthetic biology.     

Drivers of Hsp104 potentiation revealed by scanning mutagenesis of the middle domain

Hsp104, a yeast protein disaggregase, can be potentiated via numerous missense mutations at disparate locations throughout the coiled‐coil middle domain (MD). Potentiated Hsp104 variants can counter the toxicity and misfolding of TDP‐43, FUS, and α‐synuclein, proteins which are implicated in neurodegenerative disorders. However, potentiated MD variants typically exhibit off‐target toxicity. Further, it has remained confounding how numerous degenerate mutations confer potentiation, hampering engineering of therapeutic Hsp104 variants. Here, we sought to comprehensively define the key drivers of Hsp104 potentiation. Using scanning mutagenesis, we iteratively studied the effects of modulation at each position in the Hsp104 MD. Screening this library to identify enhanced variants reveals that missense mutations at 26% of positions in the MD yield variants that counter FUS toxicity. Modulation of the helix 2–helix 3/4 MD interface potentiates Hsp104, whereas mutations in the analogous helix 1–2 interface do not. Surprisingly, we find that there is a higher likelihood of enhancing Hsp104 activity against human disease substrates than impairing Hsp104 native function. We find that single mutations can broadly destabilize the MD structure and lead to functional potentiation, suggesting this may be a common mechanism conferring Hsp104 potentiation. Using this approach, we have demonstrated that modulation of the MD can yield engineered variants with decreased off‐target effects.     

Mutagenesis of mNeptune Red-Shifts Emission Spectrum to 681-685 nm

GFP-like fluorescent proteins with diverse emission wavelengths have been developed through mutagenesis, offering many possible choices in cellular and tissue imaging, such as multi-targets imaging, deep tissue imaging that require longer emission wavelength. Here, we utilized a combined approach of random mutation and structure-based rational design to develop new NIR fluorescent proteins on the basis of a far-red fluorescent protein, mNeptune (Ex/Em: 600/650 nm). We created a number of new monomeric NIR fluorescent proteins with the emission range of 681–685 nm, which exhibit the largest Stocks shifts (77–80 nm) compared to other fluorescent proteins. Among them, mNeptune681 and mNeptune684 exhibit more than 30 nm redshift in emission relative to mNeptune, owing to the major role of the extensive hydrogen-bond network around the chromophore and contributions of individual mutations to the observed redshift. Furthermore, the two variants still maintain monomeric state in solution, which is a trait crucial for their use as protein tags. In conclusion, our results suggest that there is untapped potential for developing fluorescent proteins with desired properties.     

User-friendly algorithms for estimating completeness and diversity in randomized protein-encoding libraries

Directed evolution of proteins depends on the production of molecular diversity by random mutagenesis. While a number of methods have been developed for introducing this diversity, the best ways to sample it are not always clear. Here we used simple statistics to analyse completeness and diversity in randomized libraries generated by oligonucleotide-directed mutagenesis, error-prone polymerase chain reaction (epPCR) and in vitro recombination of highly homologous sequences. For oligonucleotide-directed mutagenesis, we derive equations to estimate how complete a given library is expected to be and also to predict the size of library required to give a fixed probability of being 100% complete. We describe the statistical bases for computer programs which estimate the number of distinct variants represented in epPCR and shuffled libraries, dubbed PEDEL and DRIVeR, respectively. These programs allow the user to calculate (rather than guess) the diversity represented in a given library and also provide empirical guidelines for maximizing this diversity. PEDEL and DRIVeR are available at www.bio.cam.ac.uk/ approximately blackburn/stats.html.