How protein stability and new functions trade off

Thermal stabilization of the deglycating enzyme Amadoriase I by rational design

Amadoriases are a class of FAD-dependent enzymes that are found in fungi, yeast and bacteria and that are able to hydrolyze glycated amino acids, cleaving the sugar moiety from the amino acidic portion. So far, engineered Amadoriases have mostly found practical application in the measurement of the concentration of glycated albumin in blood samples. However, these engineered forms of Amadoriases show relatively low absolute activity and stability levels, which affect their conditions of use. Therefore, enzyme stabilization is desirable prior to function-altering molecular engineering. In this work, we describe a rational design strategy based on a computational screening method to evaluate a library of potentially stabilizing disulfide bonds. Our approach allowed the identification of two thermostable Amadoriase I mutants (SS03 and SS17) featuring a significantly higher T₅₀ (55.3 °C and 60.6 °C, respectively) compared to the wild-type enzyme (52.4 °C). Moreover, SS17 shows clear hyperstabilization, with residual activity up to 95 °C, whereas the wild-type enzyme is fully inactive at 55 °C. Our computational screening method can therefore be considered as a promising approach to expedite the design of thermostable enzymes.

FoldX as Protein Engineering Tool: Better Than Random Based Approaches?

Improving protein stability is an important goal for basic research as well as for clinical and industrial applications but no commonly accepted and widely used strategy for efficient engineering is known. Beside random approaches like error prone PCR or physical techniques to stabilize proteins, e.g. by immobilization, in silico approaches are gaining more attention to apply target-oriented mutagenesis. In this review different algorithms for the prediction of beneficial mutation sites to enhance protein stability are summarized and the advantages and disadvantages of FoldX are highlighted. The question whether the prediction of mutation sites by the algorithm FoldX is more accurate than random based approaches is addressed.

Directed evolution to improve protein folding in vivo

Recently, several innovative approaches have been developed that allow one to directly screen or select for improved protein folding in the cellular context. These methods have the potential of not just leading to a better understanding of the in vivo folding process, they may also allow for improved production of proteins of biotechnological interest.

Principles of Protein Stability and Their Application in Computational Design

Proteins are increasingly used in basic and applied biomedical research. Many proteins, however, are only marginally stable and can be expressed in limited amounts, thus hampering research and applications. Research has revealed the thermodynamic, cellular, and evolutionary principles and mechanisms that underlie marginal stability. With this growing understanding, computational stability design methods have advanced over the past two decades starting from methods that selectively addressed only some aspects of marginal stability. Current methods are more general and, by combining phylogenetic analysis with atomistic design, have shown drastic improvements in solubility, thermal stability, and aggregation resistance while maintaining the protein's primary molecular activity. Stability design is opening the way to rational engineering of improved enzymes, therapeutics, and vaccines and to the application of protein design methodology to large proteins and molecular activities that have proven challenging in the past.

Enhancement of the thermal and alkaline pH stability of Escherichia coli lysine decarboxylase for efficient cadaverine production

OBJECTIVE

To enhance the thermal and alkaline pH stability of the lysine decarboxylase from Escherichia coli (CadA) by engineering the decameric interface and explore its potential for industrial applications.

RESULTS

The mutant T88S was designed for improved structural stability by computational analysis. The optimal pH and temperature of T88S were 7.0 and 55 °C (5.5 and 50 °C for wild-type). T88S showed higher thermostability with a 2.9-fold increase in the half-life at 70 °C (from 11 to 32 min) and increased melting temperature (from 76 to 78 °C). Additionally, the specific activity and pH stability (residual activity after 10 h incubation) of T88S at pH 8.0 were increased to 164 U/mg and 78% (58 U/mg and 57% for wild-type). The productivity of cadaverine with T88S (284 g L-lysine L^-1 and 5 g DCW L^-1) was 40 g L^-1 h^-1, in contrast to 28 g L^-1 h^-1 with wild-type.

CONCLUSION

The mutant T88S showed high thermostability, pH stability, and activity at alkaline pH, indicating that this mutant is a promising biocatalyst for industrial production of cadaverine.

Identification of a hot-spot to enhance Candida rugosa lipase thermostability by rational design methods

Lipase is one of the most widely used classes of enzymes in biotechnological applications and organic chemistry. Candida rugosa lipases (CRL) can catalyze hydrolysis, esterification and transesterification with high regio-, stereo- and enantio-selectivity. However, thermal inactivation above 45 °C limits CRL's applications. Studies on improving the thermal stability of CRL are often limited by its slow-growing eukaryotic expression host, which is not suitable for large-scale screening. Identification of thermally stable mutants by rational design, regarded as an efficient substitution of experimental efforts, would provide a method for site-directed improvement of CRL. In this study, mutation-induced stability changes in CRL Lip1 were predicted by three rational design methods. Followed by conservative analyses and functional region exclusion, five mutants of a hot-spot, Asp457Phe, Asp457Trp, Asp457Met, Asp457Leu, and Asp457Tyr, were identified and prepared for enzymatic characterization. These five mutants increased the apparent melting temperature of Lip1 from 7.4 °C to 9.3 °C, with the most thermostable mutant, Asp457Phe, exhibiting a 5.5-fold longer half-life at 50 °C and a 10 °C increase in optimum temperature. Furthermore, pH stability of Lip1 was also enhanced due to the introduction of Asp457Phe mutation. The study demonstrates that thermally stable mutants of CRL could be identified with limited experimental efforts using rational design methods.

The thermostability of Candida rugosa lipase expressed in a eukaryotic host is enhanced with limited experimental effort based on rational design methods.

Adaptations of an RNA virus to increasing thermal stress

Environments can change in incremental fashions, where a shift from one state to another occurs over multiple organismal generations. The rate of the environmental change is expected to influence how and how well populations adapt to the final environmental state. We used a model system, the lytic RNA bacteriophage Φ6, to investigate this question empirically. We evolved viruses for thermostability by exposing them to heat shocks that increased to a maximum temperature at different rates. We observed increases in the ability of many heat-shocked populations to survive high temperature heat shocks. On their first exposure to the highest temperature, populations that experienced a gradual increase in temperature had higher average survival than populations that experienced a rapid temperature increase. However, at the end of the experiment, neither the survival of populations at the highest temperature nor the number of mutations per population varied significantly according to the rate of thermal change. We also evaluated mutations from the endpoint populations for their effects on viral thermostability and growth. As expected, some mutations did increase viral thermostability. However, other mutations decreased thermostability but increased growth rate, suggesting that benefits of an increased replication rate may have sometimes outweighed the benefits of enhanced thermostability. Our study highlights the importance of considering the effects of multiple selective pressures, even in environments where a single factor changes.

Improvement in the Thermostability of a β-Amino Acid Converting ω-Transaminase by Using FoldX

ω-Transaminases (ω-TAs) are important biocatalysts for the synthesis of active, chiral pharmaceutical ingredients containing amino groups, such as β-amino acids, which are important in peptidomimetics and as building blocks for drugs. However, the application of ω-TAs is limited by the availability and stability of enzymes with high conversion rates. One strategy for the synthesis and optical resolution of β-phenylalanine and other important aromatic β-amino acids is biotransformation by utilizing an ω-transaminase from Variovorax paradoxus. We designed variants of this ω-TA to gain higher process stability on the basis of predictions calculated by using the FoldX software. We herein report the first thermostabilization of a nonthermostable S-selective ω-TA by FoldX-guided site-directed mutagenesis. The melting point (T_m ) of our best-performing mutant was increased to 59.3 °C, an increase of 4.0 °C relative to the T_m value of the wild-type enzyme, whereas the mutant fully retained its specific activity.

Evolutionary Analysis of Pectin Lyases of the Genus Colletotrichum

Pectin lyases (PNLs) are important enzymes that are involved in plant cell wall degradation during the infection process. Colletotrichum is a diverse genus of fungi, which allows the study of the evolution of PNLs and their possible role in pathogen-host interactions and lifestyle adaptations. The phylogenetic reconstruction of PNLs from Colletotrichum and analysis of selection pressures showed the formation of protein lineages by groups of species with different selection pressures and specific patterns. The analysis of positive selection at individual sites using different methods allowed for the identification of three codons with evidence of positive selection in the oligosaccharide-binding region and two codons on the antiparallel sheet, which may influence the interaction with the substrate. Seven codons on the surface of the protein, mainly in the peripheral helices of the PNLs, could have an important function in evasion of plant defenses, as has been proposed in other enzymes. According to our results, it is possible that events of genetic duplication occurred in ancestral lines, followed by episodes of genetic diversification and gene loss, probably influenced by differences in the composition of the host cell wall. Additionally, different patterns of evolution in Colletotrichum appear to be molded by a strong purifying selection and positive selection episodes that forged the observed evolutionary patterns, possibly influenced by host interaction or substrate specificity. This work represents a starting point for the study of sites that may be important for evasion of plant defenses and biotechnological purposes.

ΔFlucs: Brighter Photinus pyralis firefly luciferases identified by surveying consecutive single amino acid deletion mutations in a thermostable variant

The bright bioluminescence catalyzed by Photinus pyralis firefly luciferase (Fluc) enables a vast array of life science research such as bio imaging in live animals and sensitive in vitro diagnostics. The effectiveness of such applications is improved using engineered enzymes that to date have been constructed using amino acid substitutions. We describe ΔFlucs: consecutive single amino acid deletion mutants within six loop structures of the bright and thermostable ×11 Fluc. Deletion mutations are a promising avenue to explore new sequence and functional space and isolate novel mutant phenotypes. However, this method is often overlooked and to date there have been no surveys of the effects of consecutive single amino acid deletions in Fluc. We constructed a large semi-rational ΔFluc library and isolated significantly brighter enzymes after finding ×11 Fluc activity was largely tolerant to deletions. Targeting an "omega-loop" motif (T352-G360) significantly enhanced activity, altered kinetics, reduced Km for D-luciferin_, altered emission colors, and altered substrate specificity for redshifted analog DL-infraluciferin. Experimental and in silico analyses suggested remodeling of the Ω-loop impacts on active site hydrophobicity to increase light yields. This work demonstrates the further potential of deletion mutations, which can generate useful Fluc mutants and broaden the palette of the biomedical and biotechnological bioluminescence enzyme toolbox.

Role of an N-terminal extension in stability and catalytic activity of a hyperthermostable α/β hydrolase fold esterase

The carbohydrate esterase family 7 (CE7) enzymes catalyze the deacetylation of acetyl esters of a broad range of alcohols and is unique in its activity towards cephalosporin C. The CE7 fold contains a conserved N-terminal extension that distinguishes it from the canonical α/β hydrolase fold. The hexameric quaternary structure indicates that the N-terminus may affect activity and specificity by controlling access of substrates to the buried active sites via an entrance tunnel. In this context, we characterized the catalytic parameters, conformation and thermal stability of two truncation variants lacking four and ten residues of the N-terminal region of the hyperthermostable Thermotoga maritima CE7 acetyl esterase (TmAcE). The truncations did not affect the secondary structure or the fold but modulated the oligomerization dynamics. A modest increase was observed in substrate specificity for acetylated xylose compared with acetylated glucose. A drastic reduction of ~30-40°C in the optimum temperature for activity of the variants indicated lower thermal stability. The loss of hyperthermostability appears to be an indirect effect associated with an increase in the conformational flexibility of an otherwise rigid neighboring loop containing a catalytic triad residue. The results suggest that the N-terminal extension was evolutionarily selected to preserve the stability of the enzyme.

SSADH deficiency in an Italian family: a novel ALDH5A1 gene mutation affecting the succinic semialdehyde substrate binding site

SSADH deficiency (SSADHD) is a rare autosomal recessively inherited metabolic disorder. It is associated with mutations of ALDH5A1 gene, coding for the homotetrameric enzyme SSADH. This enzyme is involved in γ-aminobutyric acid (GABA) catabolism, since it oxidizes succinic semialdehyde (SSA) to succinate. Mutations in ALDH5A1 gene result in the abnormal accumulation of γ-hydroxybutyrate (GHB), which is pathognomonic of SSADHD. In the present report, diagnosis of SSADHD in a three-month-old female was achieved by detection of high levels of GHB in urine. Sequence analysis of ALDH5A1 gene showed that the patient was a compound heterozygote for c.1226G > A (p.G409D) and the novel missense mutation, c.1498G > C (p.V500 L). By ALDH5A1 gene expression in transiently transfected HEK293 cells and enzyme activity assays, we demonstrate that the p.V500 L mutation, despite being conservative, produces complete loss of enzyme activity. In silico protein modelling analysis and evaluation of tetramer destabilizing energies suggest that structural impairment and partial occlusion of the access channel to the active site affect enzyme activity. These findings add further knowledge on the missense mutations associated with SSADHD and the molecular mechanisms underlying the loss of the enzyme activity.

Enzyme Efficiency but Not Thermostability Drives Cefotaxime Resistance Evolution in TEM-1 β-Lactamase

A leading intellectual challenge in evolutionary genetics is to identify the specific phenotypes that drive adaptation. Enzymes offer a particularly promising opportunity to pursue this question, because many enzymes' contributions to organismal fitness depend on a comparatively small number of experimentally accessible properties. Moreover, on first principles the demands of enzyme thermostability stand in opposition to the demands of catalytic activity. This observation, coupled with the fact that enzymes are only marginally thermostable, motivates the widely held hypothesis that mutations conferring functional improvement require compensatory mutations to restore thermostability. Here, we explicitly test this hypothesis for the first time, using four missense mutations in TEM-1 β-lactamase that jointly increase cefotaxime Minimum Inhibitory Concentration (MIC) ∼1500-fold. First, we report enzymatic efficiency (kcat/KM) and thermostability (Tm, and thence ΔG of folding) for all combinations of these mutations. Next, we fit a quantitative model that predicts MIC as a function of kcat/KM and ΔG. While kcat/KM explains ∼54% of the variance in cefotaxime MIC (∼92% after log transformation), ΔG does not improve explanatory power of the model. We also find that cefotaxime MIC rises more slowly in kcat/KM than predicted. Several explanations for these discrepancies are suggested. Finally, we demonstrate substantial sign epistasis in MIC and kcat/KM, and antagonistic pleiotropy between phenotypes, in spite of near numerical additivity in the system. Thus constraints on selectively accessible trajectories, as well as limitations in our ability to explain such constraints in terms of underlying mechanisms are observed in a comparatively "well-behaved" system.

Prokaryotic and Highly-Repetitive WD40 Proteins: A Systematic Study

As an ancient protein family, the WD40 repeat proteins often play essential roles in fundamental cellular processes in eukaryotes. Although investigations of eukaryotic WD40 proteins have been frequently reported, prokaryotic ones remain largely uncharacterized. In this paper, we report a systematic analysis of prokaryotic WD40 proteins and detailed comparisons with eukaryotic ones. About 4,000 prokaryotic WD40 proteins have been identified, accounting for 6.5% of all WD40s. While their abundances are less than 0.1% in most prokaryotes, they are enriched in certain species from Cyanobacteria and Planctomycetes, and participate in various functions such as prokaryotic signal transduction and nutrient synthesis. Comparisons show that a higher proportion of prokaryotic WD40s tend to contain multiple WD40 domains and a large number of hydrogen bond networks. The observation that prokaryotic WD40 proteins tend to show high internal sequence identity suggests that a substantial proportion of them (~20%) should be formed by recent or young repeat duplication events. Further studies demonstrate that the very young WD40 proteins, i.e., Highly-Repetitive WD40s, should be of higher stability. Our results have presented a catalogue of prokaryotic WD40 proteins, and have shed light on their evolutionary origins.

Adaptive HIV-1 evolutionary trajectories are constrained by protein stability

Despite the use of combination antiretroviral drugs for the treatment of HIV-1 infection, the emergence of drug resistance remains a problem. Resistance may be conferred either by a single mutation or a concerted set of mutations. The involvement of multiple mutations can arise due to interactions between sites in the amino acid sequence as a consequence of the need to maintain protein structure. To better understand the nature of such epistatic interactions, we reconstructed the ancestral sequences of HIV-1’s Pol protein, and traced the evolutionary trajectories leading to mutations associated with drug resistance. Using contemporary and ancestral sequences we modelled the effects of mutations (i.e. amino acid replacements) on protein structure to understand the functional effects of residue changes. Although the majority of resistance-associated sequences tend to destabilise the protein structure, we find there is a general tendency for protein stability to decrease across HIV-1’s evolutionary history. That a similar pattern is observed in the non-drug resistance lineages indicates that non-resistant mutations, for example, associated with escape from the immune response, also impacts on protein stability. Maintenance of optimal protein structure therefore represents a major constraining factor to the evolution of HIV-1.

Computational tools help improve protein stability but with a solubility tradeoff

Accurately predicting changes in protein stability upon amino acid substitution is a much sought after goal. Destabilizing mutations are often implicated in disease, whereas stabilizing mutations are of great value for industrial and therapeutic biotechnology. Increasing protein stability is an especially challenging task, with random substitution yielding stabilizing mutations in only ∼2% of cases. To overcome this bottleneck, computational tools that aim to predict the effect of mutations have been developed; however, achieving accuracy and consistency remains challenging. Here, we combined 11 freely available tools into a meta-predictor (meieringlab.uwaterloo.ca/stabilitypredict/). Validation against ∼600 experimental mutations indicated that our meta-predictor has improved performance over any of the individual tools. The meta-predictor was then used to recommend 10 mutations in a previously designed protein of moderate thermodynamic stability, ThreeFoil. Experimental characterization showed that four mutations increased protein stability and could be amplified through ThreeFoil's structural symmetry to yield several multiple mutants with >2-kcal/mol stabilization. By avoiding residues within functional ties, we could maintain ThreeFoil's glycan-binding capacity. Despite successfully achieving substantial stabilization, however, almost all mutations decreased protein solubility, the most common cause of protein design failure. Examination of the 600-mutation data set revealed that stabilizing mutations on the protein surface tend to increase hydrophobicity and that the individual tools favor this approach to gain stability. Thus, whereas currently available tools can increase protein stability and combining them into a meta-predictor yields enhanced reliability, improvements to the potentials/force fields underlying these tools are needed to avoid gaining protein stability at the cost of solubility.

Ancestral Haloalkane Dehalogenases Show Robustness and Unique Substrate Specificity

Ancestral sequence reconstruction (ASR) represents a powerful approach for empirical testing structure-function relationships of diverse proteins. We employed ASR to predict sequences of five ancestral haloalkane dehalogenases (HLDs) from the HLD-II subfamily. Genes encoding the inferred ancestral sequences were synthesized and expressed in Escherichia coli, and the resurrected ancestral enzymes (AncHLD1-5) were experimentally characterized. Strikingly, the ancestral HLDs exhibited significantly enhanced thermodynamic stability compared to extant enzymes (ΔT_m up to 24 °C), as well as higher specific activities with preference for short multi-substituted halogenated substrates. Moreover, multivariate statistical analysis revealed a shift in the substrate specificity profiles of AncHLD1 and AncHLD2. This is extremely difficult to achieve by rational protein engineering. The study highlights that ASR is an efficient approach for the development of novel biocatalysts and robust templates for directed evolution.

Exploring the past and the future of protein evolution with ancestral sequence reconstruction: the 'retro' approach to protein engineering

A central goal in molecular evolution is to understand the ways in which genes and proteins evolve in response to changing environments. In the absence of intact DNA from fossils, ancestral sequence reconstruction (ASR) can be used to infer the evolutionary precursors of extant proteins. To date, ancestral proteins belonging to eubacteria, archaea, yeast and vertebrates have been inferred that have been hypothesized to date from between several million to over 3 billion years ago. ASR has yielded insights into the early history of life on Earth and the evolution of proteins and macromolecular complexes. Recently, however, ASR has developed from a tool for testing hypotheses about protein evolution to a useful means for designing novel proteins. The strength of this approach lies in the ability to infer ancestral sequences encoding proteins that have desirable properties compared with contemporary forms, particularly thermostability and broad substrate range, making them good starting points for laboratory evolution. Developments in technologies for DNA sequencing and synthesis and computational phylogenetic analysis have led to an escalation in the number of ancient proteins resurrected in the last decade and greatly facilitated the use of ASR in the burgeoning field of synthetic biology. However, the primary challenge of ASR remains in accurately inferring ancestral states, despite the uncertainty arising from evolutionary models, incomplete sequences and limited phylogenetic trees. This review will focus, firstly, on the use of ASR to uncover links between sequence and phenotype and, secondly, on the practical application of ASR in protein engineering.

Thermal stability and kinetic constants for 129 variants of a family 1 glycoside hydrolase reveal that enzyme activity and stability can be separately designed

Accurate modeling of enzyme activity and stability is an important goal of the protein engineering community. However, studies seeking to evaluate current progress are limited by small data sets of quantitative kinetic constants and thermal stability measurements. Here, we report quantitative measurements of soluble protein expression in E. coli, thermal stability, and Michaelis-Menten constants (k_cat, K_M, and k_cat/K_M) for 129 designed mutants of a glycoside hydrolase. Statistical analyses reveal that functional T_m is independent of k_cat, K_M, and k_cat/K_M in this system, illustrating that an individual mutation can modulate these functional parameters independently. In addition, this data set is used to evaluate computational predictions of protein stability using the established Rosetta and FoldX algorithms. Predictions for both are found to correlate only weakly with experimental measurements, suggesting improvements are needed in the underlying algorithms.

Mistranslation can enhance fitness through purging of deleterious mutations

Phenotypic mutations are amino acid changes caused by mistranslation. How phenotypic mutations affect the adaptive evolution of new protein functions is unknown. Here we evolve the antibiotic resistance protein TEM-1 towards resistance on the antibiotic cefotaxime in an Escherichia coli strain with a high mistranslation rate. TEM-1 populations evolved in such strains endow host cells with a general growth advantage, not only on cefotaxime but also on several other antibiotics that ancestral TEM-1 had been unable to deactivate. High-throughput sequencing of TEM-1 populations shows that this advantage is associated with a lower incidence of weakly deleterious genotypic mutations. Our observations show that mistranslation is not just a source of noise that delays adaptive evolution. It could even facilitate adaptive evolution by exacerbating the effects of deleterious mutations and leading to their more efficient purging. The ubiquity of mistranslation and its effects render mistranslation an important factor in adaptive protein evolution.

Mistranslation results in amino acid changes in proteins known as phenotypic mutations and these occur at a much higher rate than DNA mutations. Here, the authors show that mistranslation can increase the response to directional selection by exacerbating the fitness effects of deleterious DNA mutations.

Phenotypic heterogeneity promotes adaptive evolution

Genetically identical cells frequently display substantial heterogeneity in gene expression, cellular morphology and physiology. It has been suggested that by rapidly generating a subpopulation with novel phenotypic traits, phenotypic heterogeneity (or plasticity) accelerates the rate of adaptive evolution in populations facing extreme environmental challenges. This issue is important as cell-to-cell phenotypic heterogeneity may initiate key steps in microbial evolution of drug resistance and cancer progression. Here, we study how stochastic transitions between cellular states influence evolutionary adaptation to a stressful environment in yeast Saccharomyces cerevisiae. We developed inducible synthetic gene circuits that generate varying degrees of expression stochasticity of an antifungal resistance gene. We initiated laboratory evolutionary experiments with genotypes carrying different versions of the genetic circuit by exposing the corresponding populations to gradually increasing antifungal stress. Phenotypic heterogeneity altered the evolutionary dynamics by transforming the adaptive landscape that relates genotype to fitness. Specifically, it enhanced the adaptive value of beneficial mutations through synergism between cell-to-cell variability and genetic variation. Our work demonstrates that phenotypic heterogeneity is an evolving trait when populations face a chronic selection pressure. It shapes evolutionary trajectories at the genomic level and facilitates evolutionary rescue from a deteriorating environmental stress.

Phenotypic heterogeneity of genetically identical cells can generate nonheritable variation in a population. Is this heterogeneity favorable for microbes? In a changing environment, the answer is a definite yes. While scholars have argued that stochastically generated variation precedes genetic changes and thereby facilitate the evolution of complex traits, this idea has remained disputed, not least because of the shortage of experimental studies. We address this long-standing and controversial issue by integrating synthetic biology, laboratory experimental evolution, and genomic analyses. We explicitly tested the mechanisms whereby phenotypic heterogeneity may promote evolvability. Our work demonstrates that phenotypic heterogeneity facilitates evolutionary rescue from deteriorating environmental stress by generating individuals with exceptionally high fitness. Remarkably, elevated phenotypic heterogeneity evolves as a direct response to stress and thereby it promotes evolution of rare combinations of mutations. These results indicate that phenotypic heterogeneity might have an important role in the evolution of key innovations.

Improved thermostability of an acidic xylanase from Aspergillus sulphureus by combined disulphide bridge introduction and proline residue substitution

As a feed additive, xylanase has been widely applied in the feed of monogastric animals, which contains multiple plant polysaccharides. However, during feed manufacture, the high pelleting temperatures challenge wild-type xylanases. The aim of this study was to improve the thermostability of Aspergillus sulphureus acidic xylanase. According to the predicted protein structure, a series of disulphide bridges and proline substitutions were created in the xylanase by PCR, and the mutants were expressed in Pichia pastoris. Enzyme properties were evaluated following chromatographic purification. All the recombinant enzymes showed optima at pH 3.0 and 50 °C or 55 °C and better resistance to some chemicals except for CuSO₄. The specific activity of the xylanase was decreased by introduction of the mutations. Compared to the wild-type enzyme, a combined mutant, T53C-T142C/T46P, with a disulphide bond at 53–142 and a proline substitution at 46, showed a 22-fold increase of half-life at 60 °C. In a 10-L fermentor, the maximal xylanase activity of T53C-T142C/T46P reached 1,684 U/mL. It was suggested that the T53C-T142C/T46P mutant xylanase had excellent thermostability characteristics and could be a prospective additive in feed manufacture.

Efficient affinity maturation of antibody variable domains requires co-selection of compensatory mutations to maintain thermodynamic stability

The ability of antibodies to accumulate affinity-enhancing mutations in their complementarity-determining regions (CDRs) without compromising thermodynamic stability is critical to their natural function. However, it is unclear if affinity mutations in the hypervariable CDRs generally impact antibody stability and to what extent additional compensatory mutations are required to maintain stability during affinity maturation. Here we have experimentally and computationally evaluated the functional contributions of mutations acquired by a human variable (V_H) domain that was evolved using strong selections for enhanced stability and affinity for the Alzheimer’s Aβ42 peptide. Interestingly, half of the key affinity mutations in the CDRs were destabilizing. Moreover, the destabilizing effects of these mutations were compensated for by a subset of the affinity mutations that were also stabilizing. Our findings demonstrate that the accumulation of both affinity and stability mutations is necessary to maintain thermodynamic stability during extensive mutagenesis and affinity maturation in vitro, which is similar to findings for natural antibodies that are subjected to somatic hypermutation in vivo. These findings for diverse antibodies and antibody fragments specific for unrelated antigens suggest that the formation of the antigen-binding site is generally a destabilizing process and that co-enrichment for compensatory mutations is critical for maintaining thermodynamic stability.

A Gradient of Sitewise Diversity Promotes Evolutionary Fitness for Binder Discovery in a Three-Helix Bundle Protein Scaffold

Engineered proteins provide clinically and industrially impactful molecules and utility within fundamental research, yet inefficiencies in discovering lead variants with new desired functionality, while maintaining stability, hinder progress. Improved function, which can result from a few strategic mutations, is fundamentally separate from discovering novel function, which often requires large leaps in sequence space. While a highly diverse combinatorial library covering immense sequence space would empower protein discovery, the ability to sample only a minor subset of sequence space and the typical destabilization of random mutations preclude this strategy. A balance must be reached. At library scale, compounding several destabilizing mutations renders many variants unable to properly fold and devoid of function. Broadly searching sequence space while reducing the level of destabilization may enhance evolution. We exemplify this balance with affibody, a three-helix bundle protein scaffold. Using natural ligand data sets, stability and structural computations, and deep sequencing of thousands of binding variants, a protein library was designed on a sitewise basis with a gradient of mutational levels across 29% of the protein. In direct competition of biased and uniform libraries, both with 1 × 10⁹ variants, for discovery of 6 × 10⁴ ligands (5 × 10³ clusters) toward seven targets, biased amino acid frequency increased ligand discovery 13 ± 3-fold. Evolutionarily favorable amino acids, both globally and site-specifically, are further elucidated. The sitewise amino acid bias aids evolutionary discovery by reducing the level of mutant destabilization as evidenced by a midpoint of denaturation (62 ± 4 °C) 15 °C higher than that of unbiased mutants (47 ± 11 °C; p < 0.001). Sitewise diversification, identified by high-throughput evolution and rational library design, improves discovery efficiency.

Trade-offs between enzyme fitness and solubility illuminated by deep mutational scanning

Proteins are marginally stable, and an understanding of the sequence determinants for improved protein solubility is highly desired. For enzymes, it is well known that many mutations that increase protein solubility decrease catalytic activity. These competing effects frustrate efforts to design and engineer stable, active enzymes without laborious high-throughput activity screens. To address the trade-off between enzyme solubility and activity, we performed deep mutational scanning using two different screens/selections that purport to gauge protein solubility for two full-length enzymes. We assayed a TEM-1 beta-lactamase variant and levoglucosan kinase (LGK) using yeast surface display (YSD) screening and a twin-arginine translocation pathway selection. We then compared these scans with published experimental fitness landscapes. Results from the YSD screen could explain 37% of the variance in the fitness landscapes for one enzyme. Five percent to 10% of all single missense mutations improve solubility, matching theoretical predictions of global protein stability. For a given solubility-enhancing mutation, the probability that it would retain wild-type fitness was correlated with evolutionary conservation and distance to active site, and anticorrelated with contact number. Hybrid classification models were developed that could predict solubility-enhancing mutations that maintain wild-type fitness with an accuracy of 90%. The downside of using such classification models is the removal of rare mutations that improve both fitness and solubility. To reveal the biophysical basis of enhanced protein solubility and function, we determined the crystallographic structure of one such LGK mutant. Beyond fundamental insights into trade-offs between stability and activity, these results have potential biotechnological applications.

Evolutionary interplay between structure, energy and epistasis in the coat protein of the ϕX174 phage family

Viral capsids are structurally constrained by interactions among the amino acids (AAs) of their constituent proteins. Therefore, epistasis is expected to evolve among physically interacting sites and to influence the rates of substitution. To study the evolution of epistasis, we focused on the major structural protein of the ϕX174 phage family by first reconstructing the ancestral protein sequences of 18 species using a Bayesian statistical framework. The inferred ancestral reconstruction differed at eight AAs, for a total of 256 possible ancestral haplotypes. For each ancestral haplotype and the extant species, we estimated, in silico, the distribution of free energies and epistasis of the capsid structure. We found that free energy has not significantly increased but epistasis has. We decomposed epistasis up to fifth order and found that higher-order epistasis sometimes compensates pairwise interactions making the free energy seem additive. The dN/dS ratio is low, suggesting strong purifying selection, and that structure is under stabilizing selection. We synthesized phages carrying ancestral haplotypes of the coat protein gene and measured their fitness experimentally. Our findings indicate that stabilizing mutations can have higher fitness, and that fitness optima do not necessarily coincide with energy minima.

Constrained evolution of a bispecific enzyme: lessons for biocatalyst design

Analysis of the natural evolution of bispecificity in triazine hydrolase highlights the importance of epistasis in protein engineering and evolution.

Ultrahigh-throughput-directed enzyme evolution by absorbance-activated droplet sorting (AADS)

Ultrahigh-throughput screening, in which members of enzyme libraries compartmentalized in water-in-oil emulsion droplets are assayed, has emerged as a powerful format for directed evolution and functional metagenomics but is currently limited to fluorescence readouts. Here we describe a highly efficient microfluidic absorbance-activated droplet sorter (AADS) that extends the range of assays amenable to this approach. Using this module, microdroplets can be sorted based on absorbance readout at rates of up to 300 droplets per second (i.e., >1 million droplets per hour). To validate this device, we implemented a miniaturized coupled assay for NAD⁺-dependent amino acid dehydrogenases. The detection limit (10 μM in a coupled assay producing a formazan dye) enables accurate kinetic readouts sensitive enough to detect a minimum of 1,300 turnovers per enzyme molecule, expressed in a single cell, and released by lysis within a droplet. Sorting experiments showed that the AADS successfully enriched active variants up to 2,800-fold from an overwhelming majority of inactive ones at ∼100 Hz. To demonstrate the utility of this module for protein engineering, two rounds of directed evolution were performed to improve the activity of phenylalanine dehydrogenase toward its native substrate. Fourteen hits showed increased activity (improved >4.5-fold in lysate; k_cat increased >2.7-fold), soluble protein expression levels (up 60%), and thermostability (T_m, 12 °C higher). The AADS module makes the most widely used optical detection format amenable to screens of unprecedented size, paving the way for the implementation of chromogenic assays in droplet microfluidics workflows.

Circumventing the stability-function trade-off in an engineered FN3 domain

The favorable biophysical attributes of non-antibody scaffolds make them attractive alternatives to monoclonal antibodies. However, due to the well-known stability-function trade-off, these gains tend to be marginal after functional selection. A notable example is the fibronectin Type III (FN3) domain, FNfn10, which has been previously evolved to bind lysozyme with 1 pM affinity (FNfn10-α-lys), but suffers from poor thermodynamic and kinetic stability. To explore this stability-function compromise further, we grafted the lysozyme-binding loops from FNfn10-α-lys onto our previously engineered, ultra-stable FN3 scaffold, FN3con. The resulting variant (FN3con-α-lys) bound lysozyme with a markedly reduced affinity, but retained high levels of thermal stability. The crystal structure of FNfn10-α-lys in complex with lysozyme revealed unanticipated interactions at the protein-protein interface involving framework residues of FNfn10-α-lys, thus explaining the failure to transfer binding via loop grafting. Utilizing this structural information, we redesigned FN3con-α-lys and restored picomolar binding affinity to lysozyme, while maintaining thermodynamic stability (with a thermal melting temperature 2-fold higher than that of FNfn10-α-lys). FN3con therefore provides an exceptional window of stability to tolerate deleterious mutations, resulting in a substantial advantage for functional design. This study emphasizes the utility of consensus design for the generation of highly stable scaffolds for downstream protein engineering studies.

Bridging the physical scales in evolutionary biology: from protein sequence space to fitness of organisms and populations

Bridging the gap between the molecular properties of proteins and organismal/population fitness is essential for understanding evolutionary processes. This task requires the integration of the several physical scales of biological organization, each defined by a distinct set of mechanisms and constraints, into a single unifying model. The molecular scale is dominated by the constraints imposed by the physico-chemical properties of proteins and their substrates, which give rise to trade-offs and epistatic (non-additive) effects of mutations. At the systems scale, biological networks modulate protein expression and can either buffer or enhance the fitness effects of mutations. The population scale is influenced by the mutational input, selection regimes, and stochastic changes affecting the size and structure of populations, which eventually determine the evolutionary fate of mutations. Here, we summarize the recent advances in theory, computer simulations, and experiments that advance our understanding of the links between various physical scales in biology.

A Distal Disulfide Bridge in OXA-1 β-Lactamase Stabilizes the Catalytic Center and Alters the Dynamics of the Specificity Determining Ω Loop

Widespread antibiotic resistance, particularly when mediated by broad-spectrum β-lactamases, has major implications for public health. Substitutions in the active site often allow broad-spectrum enzymes to accommodate diverse types of β-lactams. Substitutions observed outside the active site are thought to compensate for the loss of thermal stability. The OXA-1 clade of class D β-lactamases contains a pair of conserved cysteines located outside the active site that forms a disulfide bond in the periplasm. Here, the effect of the distal disulfide bond on the structure and dynamics of OXA-1 was investigated via 4 μs molecular dynamics simulations. The results reveal that the disulfide promotes the preorganized orientation of the catalytic residues and affects the conformation of the functionally important Ω loop. Furthermore, principal component analysis reveals differences in the global dynamics between the oxidized and reduced forms, especially in the motions involving the Ω loop. A dynamical network analysis indicates that, in the oxidized form, in addition to its role in ligand binding, the KTG family motif is a central hub of the global dynamics. As activity of OXA-1 has been measured only in the reduced form, we suggest that accurate assessment of its functional profile would require oxidative conditions mimicking periplasm.

The Molecular Chaperone DnaK Is a Source of Mutational Robustness

Molecular chaperones, also known as heat-shock proteins, refold misfolded proteins and help other proteins reach their native conformation. Thanks to these abilities, some chaperones, such as the Hsp90 protein or the chaperonin GroEL, can buffer the deleterious phenotypic effects of mutations that alter protein structure and function. Hsp70 chaperones use a chaperoning mechanism different from that of Hsp90 and GroEL, and it is not known whether they can also buffer mutations. Here, we show that they can. To this end, we performed a mutation accumulation experiment in Escherichia coli, followed by whole-genome resequencing. Overexpression of the Hsp70 chaperone DnaK helps cells cope with mutational load and completely avoid the extinctions we observe in lineages evolving without chaperone overproduction. Additionally, our sequence data show that DnaK overexpression increases mutational robustness, the tolerance of its clients to nonsynonymous nucleotide substitutions. We also show that this elevated mutational buffering translates into differences in evolutionary rates on intermediate and long evolutionary time scales. Specifically, we studied the evolutionary rates of DnaK clients using the genomes of E. coli, Salmonella enterica, and 83 other gamma-proteobacteria. We find that clients that interact strongly with DnaK evolve faster than weakly interacting clients. Our results imply that all three major chaperone classes can buffer mutations and affect protein evolution. They illustrate how an individual protein like a chaperone can have a disproportionate effect on the evolution of a proteome.

Smoothing a rugged protein folding landscape by sequence-based redesign

The rugged folding landscapes of functional proteins puts them at risk of misfolding and aggregation. Serine protease inhibitors, or serpins, are paradigms for this delicate balance between function and misfolding. Serpins exist in a metastable state that undergoes a major conformational change in order to inhibit proteases. However, conformational labiality of the native serpin fold renders them susceptible to misfolding, which underlies misfolding diseases such as α₁-antitrypsin deficiency. To investigate how serpins balance function and folding, we used consensus design to create conserpin, a synthetic serpin that folds reversibly, is functional, thermostable, and polymerization resistant. Characterization of its structure, folding and dynamics suggest that consensus design has remodeled the folding landscape to reconcile competing requirements for stability and function. This approach may offer general benefits for engineering functional proteins that have risky folding landscapes, including the removal of aggregation-prone intermediates, and modifying scaffolds for use as protein therapeutics.

Chromosomal rearrangements and protein globularity changes in Mycobacterium tuberculosis isolates from cerebrospinal fluid

Background

Meningitis is a major cause of mortality in tuberculosis (TB). It is not clear what factors promote central nervous system invasion and pathology but it has been reported that certain strains of Mycobacterium tuberculosis (Mtb) might have genetic traits associated with neurotropism.

Methods

In this study, we generated whole genome sequences of eight clinical strains of Mtb that were isolated from the cerebrospinal fluid (CSF) of patients presenting with tuberculous meningitis (TBM) in Malaysia, and compared them to the genomes of H37Rv and other respiratory Mtb genomes either downloaded from public databases or extracted from local sputum isolates. We aimed to find genomic features that might be distinctly different between CSF-derived and respiratory Mtb.

Results

Genome-wide comparisons revealed rearrangements (translocations, inversions, insertions and deletions) and non-synonymous SNPs in our CSF-derived strains that were not observed in the respiratory Mtb genomes used for comparison. These rearranged segments were rich in genes for PE (proline-glutamate)/PPE (proline-proline-glutamate), transcriptional and membrane proteins. Similarly, most of the ns SNPs common in CSF strains were noted in genes encoding PE/PPE proteins. Protein globularity differences were observed among mycobacteria from CSF and respiratory sources and in proteins previously reported to be associated with TB meningitis. Transcription factors and other transcription regulators featured prominently in these proteins. Homologs of proteins associated with Streptococcus pneumoniae meningitis and Neisseria meningitidis virulence were identified in neuropathogenic as well as respiratory mycobacterial spp. examined in this study.

Discussion

The occurrence of in silico genetic differences in CSF-derived but not respiratory Mtb suggests their possible involvement in the pathogenesis of TBM. However, overall findings in this comparative analysis support the postulation that TB meningeal infection is more likely to be related to the expression of multiple virulence factors on interaction with host defences than to CNS tropism associated with specific genetic traits.

Distinct Structural Elements Govern the Folding, Stability, and Catalysis in the Outer Membrane Enzyme PagP

The outer membrane enzyme PagP is indispensable for lipid A palmitoylation in Gram-negative bacteria and has been implicated in resistance to host immune defenses. PagP possesses an unusual structure for an integral membrane protein, with a highly dynamic barrel domain that is tilted with respect to the membrane normal. In addition, it contains an N-terminal amphipathic helix. Recent functional and structural studies have shown that these molecular factors are critical for PagP to carry out its function in the challenging environment of the bacterial outer membrane. However, the precise contributions of the N-helix to folding and stability and residues that can influence catalytic rates remain to be addressed. Here, we identify a sequence-dependent stabilizing role for the N-terminal helix of PagP in the measured thermodynamic stability of the barrel. Using chimeric barrel sequences, we show that the Escherichia coli PagP N-terminal helix confers 2-fold greater stability to the Salmonella typhimurium barrel. Further, we find that the W78F substitution in S. typhimurium causes a nearly 20-fold increase in the specific activity in vitro for the phospholipase reaction, compared to that of E. coli PagP. Here, phenylalanine serves as a key regulator of catalysis, possibly by increasing the reaction rate. Through coevolution analysis, we detect an interaction network between seemingly unrelated segments of this membrane protein. Exchanging the structural and functional features between homologous PagP enzymes from E. coli and S. typhimurium has provided us with an understanding of the molecular factors governing PagP stability and function.

Conformational Tinkering Drives Evolution of a Promiscuous Activity through Indirect Mutational Effects

How remote mutations can lead to changes in enzyme function at a molecular level is a central question in evolutionary biochemistry and biophysics. Here, we combine laboratory evolution with biochemical, structural, genetic, and computational analysis to dissect the molecular basis for the functional optimization of phosphotriesterase activity in a bacterial lactonase (AiiA) from the metallo-β-lactamase (MBL) superfamily. We show that a 1000-fold increase in phosphotriesterase activity is caused by a more favorable catalytic binding position of the paraoxon substrate in the evolved enzyme that resulted from conformational tinkering of the active site through peripheral mutations. A nonmutated active site residue, Phe68, was displaced by ∼3 Å through the indirect effects of two second-shell trajectory mutations, allowing molecular interactions between the residue and paraoxon. Comparative mutational scanning, i.e., examining the effects of alanine mutagenesis on different genetic backgrounds, revealed significant changes in the functional roles of Phe68 and other nonmutated active site residues caused by the indirect effects of trajectory mutations. Our work provides a quantitative measurement of the impact of second-shell mutations on the catalytic contributions of nonmutated residues and unveils the underlying intramolecular network of strong epistatic mutational relationships between active site residues and more remote residues. Defining these long-range conformational and functional epistatic relationships has allowed us to better understand the subtle, but cumulatively significant, role of second- and third-shell mutations in evolution.

Automated Structure- and Sequence-Based Design of Proteins for High Bacterial Expression and Stability

Upon heterologous overexpression, many proteins misfold or aggregate, thus resulting in low functional yields. Human acetylcholinesterase (hAChE), an enzyme mediating synaptic transmission, is a typical case of a human protein that necessitates mammalian systems to obtain functional expression. We developed a computational strategy and designed an AChE variant bearing 51 mutations that improved core packing, surface polarity, and backbone rigidity. This variant expressed at ∼2,000-fold higher levels in E. coli compared to wild-type hAChE and exhibited 20°C higher thermostability with no change in enzymatic properties or in the active-site configuration as determined by crystallography. To demonstrate broad utility, we similarly designed four other human and bacterial proteins. Testing at most three designs per protein, we obtained enhanced stability and/or higher yields of soluble and active protein in E. coli. Our algorithm requires only a 3D structure and several dozen sequences of naturally occurring homologs, and is available at http://pross.weizmann.ac.il.

•

A new computational method is used to stabilize five recalcitrant proteins

•

Designed variants show higher expression and stability with unmodified function

•

A designed human acetylcholinesterase variant expresses solubly in bacteria

•

The method is fully automated and implemented on a webserver

Simple approach for ranking structure determining residues

Mutating residues has been a common task in order to study structural properties of the protein of interest. Here, we propose and validate a simple method that allows the identification of structural determinants; i.e., residues essential for preservation of the stability of global structure, regardless of the protein topology. This method evaluates all of the residues in a 3D structure of a given globular protein by ranking them according to their connectivity and movement restrictions without topology constraints. Our results matched up with sequence-based predictors that look up for intrinsically disordered segments, suggesting that protein disorder can also be described with the proposed methodology.

Consensus protein design

A popular and successful strategy in semi-rational design of protein stability is the use of evolutionary information encapsulated in homologous protein sequences. Consensus design is based on the hypothesis that at a given position, the respective consensus amino acid contributes more than average to the stability of the protein than non-conserved amino acids. Here, we review the consensus design approach, its theoretical underpinnings, successes, limitations and challenges, as well as providing a detailed guide to its application in protein engineering.

Causes of molecular convergence and parallelism in protein evolution

To what extent is the convergent evolution of protein function attributable to convergent or parallel changes at the amino acid level? The mutations that contribute to adaptive protein evolution may represent a biased subset of all possible beneficial mutations owing to mutation bias and/or variation in the magnitude of deleterious pleiotropy. A key finding is that the fitness effects of amino acid mutations are often conditional on genetic background. This context dependence (epistasis) can reduce the probability of convergence and parallelism because it reduces the number of possible mutations that are unconditionally acceptable in divergent genetic backgrounds. Here, I review factors that influence the probability of replicated evolution at the molecular level.

Strong Selection Significantly Increases Epistatic Interactions in the Long-Term Evolution of a Protein

Epistatic interactions between residues determine a protein’s adaptability and shape its evolutionary trajectory. When a protein experiences a changed environment, it is under strong selection to find a peak in the new fitness landscape. It has been shown that strong selection increases epistatic interactions as well as the ruggedness of the fitness landscape, but little is known about how the epistatic interactions change under selection in the long-term evolution of a protein. Here we analyze the evolution of epistasis in the protease of the human immunodeficiency virus type 1 (HIV-1) using protease sequences collected for almost a decade from both treated and untreated patients, to understand how epistasis changes and how those changes impact the long-term evolvability of a protein. We use an information-theoretic proxy for epistasis that quantifies the co-variation between sites, and show that positive information is a necessary (but not sufficient) condition that detects epistasis in most cases. We analyze the “fossils” of the evolutionary trajectories of the protein contained in the sequence data, and show that epistasis continues to enrich under strong selection, but not for proteins whose environment is unchanged. The increase in epistasis compensates for the information loss due to sequence variability brought about by treatment, and facilitates adaptation in the increasingly rugged fitness landscape of treatment. While epistasis is thought to enhance evolvability via valley-crossing early-on in adaptation, it can hinder adaptation later when the landscape has turned rugged. However, we find no evidence that the HIV-1 protease has reached its potential for evolution after 9 years of adapting to a drug environment that itself is constantly changing. We suggest that the mechanism of encoding new information into pairwise interactions is central to protein evolution not just in HIV-1 protease, but for any protein adapting to a changing environment.

Evolution is often viewed as a process that occurs “mutation by mutation”, suggesting that the effect of each mutation is independent of that of others. However, in reality the effect of a mutation often depends on the context of other mutations, a dependence known as “epistasis”. Even though epistasis can constrain protein evolution, it is actually very common. Such interactions are particularly pervasive in proteins that evolve resistance to a drug via mutations that create defects, and that must be repaired with compensatory mutations. We study how epistasis between protein residues evolves over time in a new and changing environment, and compare these findings to protein evolution in a constant environment. We analyze the sequences of the human immunodeficiency virus type 1 (HIV-1) protease enzyme collected over a period of 9 years from patients treated with anti-viral drugs (as well as from patients that went untreated), and find that epistasis between residues continues to increase as more potent anti-viral drugs enter the market, while epistasis is unchanging in the proteins exposed to a constant environment. Yet, the proteins adapting to the changing landscape do not appear to be constrained by the epistatic interactions and continue to manage to evade new drugs.

Optimization of Conformational Dynamics in an Epistatic Evolutionary Trajectory

The understanding of protein evolution depends on the ability to relate the impact of mutations on molecular traits to organismal fitness. Biological activity and robustness have been regarded as important features in shaping protein evolutionary landscapes. Conformational dynamics, which is essential for protein function, has received little attention in the context of evolutionary analyses. Here we employ NMR spectroscopy, the chief experimental tool to describe protein dynamics at atomic level in solution at room temperature, to study the intrinsic dynamic features of a metallo- Β: -lactamase enzyme and three variants identified during a directed evolution experiment that led to an expanded substrate profile. We show that conformational dynamics in the catalytically relevant microsecond to millisecond timescale is optimized along the favored evolutionary trajectory. In addition, we observe that the effects of mutations on dynamics are epistatic. Mutation Gly262Ser introduces slow dynamics on several residues that surround the active site when introduced in the wild-type enzyme. Mutation Asn70Ser removes the slow dynamics observed for few residues of the wild-type enzyme, but increases the number of residues that undergo slow dynamics when introduced in the Gly262Ser mutant. These effects on dynamics correlate with the epistatic interaction between these two mutations on the bacterial phenotype. These findings indicate that conformational dynamics is an evolvable trait, and that proteins endowed with more dynamic active sites also display a larger potential for promoting evolution.

Why reinvent the wheel? Building new proteins based on ready-made parts

We protein engineers are ambivalent about evolution: on the one hand, evolution inspires us with myriad examples of biomolecular binders, sensors, and catalysts; on the other hand, these examples are seldom well-adapted to the engineering tasks we have in mind. Protein engineers have therefore modified natural proteins by point substitutions and fragment exchanges in an effort to generate new functions. A counterpoint to such design efforts, which is being pursued now with greater success, is to completely eschew the starting materials provided by nature and to design new protein functions from scratch by using de novo molecular modeling and design. While important progress has been made in both directions, some areas of protein design are still beyond reach. To this end, we advocate a synthesis of these two strategies: by using design calculations to both recombine and optimize fragments from natural proteins, we can build stable and as of yet un-sampled structures, thereby granting access to an expanded repertoire of conformations and desired functions. We propose that future methods that combine phylogenetic analysis, structure and sequence bioinformatics, and atomistic modeling may well succeed where any one of these approaches has failed on its own.

The Evolution of New Catalytic Mechanisms for Xenobiotic Hydrolysis in Bacterial Metalloenzymes

An increasing number of bacterial metalloenzymes have been shown to catalyse the breakdown of xenobiotics in the environment, while others exhibit a variety of promiscuous xenobiotic-degrading activities. Several different evolutionary processes have allowed these enzymes to gain or enhance xenobiotic-degrading activity. In this review, we have surveyed the range of xenobiotic-degrading metalloenzymes, and discuss the molecular and catalytic basis for the development of new activities. We also highlight how our increased understanding of the natural evolution of xenobiotic-degrading metalloenzymes can be been applied to laboratory enzyme design.