Consensus-Based Engineering of Protein Stability: From Intrabodies to Thermostable Enzymes

Comparative analysis of reconstructed ancestral proteins with their extant counterparts suggests primitive life had an alkaline habitat

To understand the origin and early evolution of life it is crucial to establish characteristics of the primordial environment that facilitated the emergence and evolution of life. One important environmental factor is the pH of the primordial environment. Here, we assessed the pH-dependent thermal stabilities of previously reconstructed ancestral nucleoside diphosphate kinases and ribosomal protein uS8s. The selected proteins were likely to be present in ancient organisms such as the last common ancestor of bacteria and that of archaea. We also assessed the thermal stability of homologous proteins from extant acidophilic, neutralophilic, and alkaliphilic microorganisms as a function of pH. Our results indicate that the reconstructed ancestral proteins are more akin to those of extant alkaliphilic bacteria, which display greater stability under alkaline conditions. These findings suggest that the common ancestors of bacterial and archaeal species thrived in an alkaline environment. Moreover, we demonstrate the reconstruction method employed in this study is a valuable technique for generating alkali-tolerant proteins that can be used in a variety of biotechnological and environmental applications.

MEnTaT: A machine-learning approach for the identification of mutations to increase protein stability

Enhancing protein thermal stability is important for biomedical and industrial applications as well as in the research laboratory. Here, we describe a simple machine-learning method which identifies amino acid substitutions that contribute to thermal stability based on comparison of the amino acid sequences of homologous proteins derived from bacteria that grow at different temperatures. A key feature of the method is that it compares the sequences based not simply on the amino acid identity, but rather on the structural and physicochemical properties of the side chain. The method accurately identified stabilizing substitutions in three well-studied systems and was validated prospectively by experimentally testing predicted stabilizing substitutions in a polyamine oxidase. In each case, the method outperformed the widely used bioinformatic consensus approach. The method can also provide insight into fundamental aspects of protein structure, for example, by identifying how many sequence positions in a given protein are relevant to temperature adaptation.

Engineering functional thermostable proteins using ancestral sequence reconstruction

Natural proteins are often only slightly more stable in the native state than the denatured state, and an increase in environmental temperature can easily shift the balance toward unfolding. Therefore, the engineering of proteins to improve protein stability is an area of intensive research. Thermostable proteins are required to withstand industrial process conditions, for increased shelf-life of protein therapeutics, for developing robust 'biobricks' for synthetic biology applications, and for research purposes (e.g., structure determination). In addition, thermostability buffers the often destabilizing effects of mutations introduced to improve other properties. Rational design approaches to engineering thermostability require structural information, but even with advanced computational methods, it is challenging to predict or parameterize all the relevant structural factors with sufficient precision to anticipate the results of a given mutation. Directed evolution is an alternative when structures are unavailable but requires extensive screening of mutant libraries. Recently, however, bioinspired approaches based on phylogenetic analyses have shown great promise. Leveraging the rapid expansion in sequence data and bioinformatic tools, ancestral sequence reconstruction can generate highly stable folds for novel applications in industrial chemistry, medicine, and synthetic biology. This review provides an overview of the factors important for successful inference of thermostable proteins by ancestral sequence reconstruction and what it can reveal about the determinants of stability in proteins.

Improving antibody thermostability based on statistical analysis of sequence and structural consensus data

Background

The use of Monoclonal Antibodies (MAbs) as therapeutics has been increasing over the past 30 years due to their high specificity and strong affinity towards the target. One of the major challenges towarding their use as drugs is their low thermostability, which impacts both efficacy as well as manufacturing and delivery.

Methods

To aid the design of thermally more stable mutants, consensus sequence-based method has been widely used. These methods typically have a success rate of about 50% with maximum melting temperature increment ranging from 10 to 32 °C. In order to improve the prediction performance, we have developed a new and fast MAbs specific method by adding a 3D structural layer to the consensus sequence method. This is done by analyzing the close-by residue pairs which are conserved in more than eight hundred MAbs’ 3D structures.

Results

Combining consensus sequence and structural residue pair covariance methods, we developed an in-house application for predicting human MAb thermostability to guide protein engineers to design stable molecules. Major advantage of this structural level assessment is in significantly reducing the false positives by almost half from the consensus sequence method alone. This application has shown success in designing MAb engineering panels in multiple biologics programs.

Conclusions

Our data science-based method shows impacts in Mab engineering.

Ancestral sequences of a large promiscuous enzyme family correspond to bridges in sequence space in a network representation

Evolutionary relationships of protein families can be characterized either by networks or by trees. Whereas trees allow for hierarchical grouping and reconstruction of the most likely ancestral sequences, networks lack a time axis but allow for thresholds of pairwise sequence identity to be chosen and, therefore, the clustering of family members with presumably more similar functions. Here, we use the large family of arylsulfatases and phosphonate monoester hydrolases to investigate similarities, strengths and weaknesses in tree and network representations. For varying thresholds of pairwise sequence identity, values of betweenness centrality and clustering coefficients were derived for nodes of the reconstructed ancestors to measure the propensity to act as a bridge in a network. Based on these properties, ancestral protein sequences emerge as bridges in protein sequence networks. Interestingly, many ancestral protein sequences appear close to extant sequences. Therefore, reconstructed ancestor sequences might also be interpreted as yet-to-be-identified homologues. The concept of ancestor reconstruction is compared to consensus sequences, too. It was found that hub sequences in a network, e.g. reconstructed ancestral sequences that are connected to many neighbouring sequences, share closer similarity with derived consensus sequences. Therefore, some reconstructed ancestor sequences can also be interpreted as consensus sequences.

Comprehensive mutagenesis to identify amino acid residues contributing to the difference in thermostability between two originally thermostable ancestral proteins

Further improvement of the thermostability of inherently thermostable proteins is an attractive challenge because more thermostable proteins are industrially more useful and serve as better scaffolds for protein engineering. To establish guidelines that can be applied for the rational design of hyperthermostable proteins, we compared the amino acid sequences of two ancestral nucleoside diphosphate kinases, Arc1 and Bac1, reconstructed in our previous study. Although Bac1 is a thermostable protein whose unfolding temperature is around 100°C, Arc1 is much more thermostable with an unfolding temperature of 114°C. However, only 12 out of 139 amino acids are different between the two sequences. In this study, one or a combination of amino acid(s) in Bac1 was/were substituted by a residue(s) found in Arc1 at the same position(s). The best mutant, which contained three amino acid substitutions (S108D, G116A and L120P substitutions), showed an unfolding temperature more than 10°C higher than that of Bac1. Furthermore, a combination of the other nine amino acid substitutions also led to improved thermostability of Bac1, although the effects of individual substitutions were small. Therefore, not only the sum of the contributions of individual amino acids, but also the synergistic effects of multiple amino acids are deeply involved in the stability of a hyperthermostable protein. Such insights will be helpful for future rational design of hyperthermostable proteins.

Surface salt bridges contribute to the extreme thermal stability of an FN3-like domain from a thermophilic bacterium

This study uses differential scanning calorimetry, X-ray crystallography, and molecular dynamics simulations to investigate the structural basis for the high thermal stability (melting temperature 97.5°C) of a FN3-like protein domain from thermophilic bacteria Thermoanaerobacter tengcongensis (FN3tt). FN3tt adopts a typical FN3 fold with a three-stranded beta sheet packing against a four-stranded beta sheet. We identified three solvent exposed arginine residues (R23, R25, and R72), which stabilize the protein through salt bridge interactions with glutamic acid residues on adjacent strands. Alanine mutation of the three arginine residues reduced melting temperature by up to 22°C. Crystal structures of the wild type (WT) and a thermally destabilized (∆Tm -19.7°C) triple mutant (R23L/R25T/R72I) were found to be nearly identical, suggesting that the destabilization is due to interactions of the arginine residues. Molecular dynamics simulations showed that the salt bridge interactions in the WT were stable and provided a dynamical explanation for the cooperativity observed between R23 and R25 based on calorimetry measurements. In addition, folding free energy changes computed using free energy perturbation molecular dynamics simulations showed high correlation with melting temperature changes. This work is another example of surface salt bridges contributing to the enhanced thermal stability of thermophilic proteins. The molecular dynamics simulation methods employed in this study may be broadly useful for in silico surface charge engineering of proteins.

Mining and Statistical Modeling of Natural and Variant Class IIa Bacteriocins Elucidate Activity and Selectivity Profiles across Species

Class IIa bacteriocin antimicrobial peptides (AMPs) are a compelling alternative to current antimicrobials because of potential specific activity toward antibiotic-resistant bacteria, including vancomycin-resistant enterococci. Engineering of these molecules would be enhanced by a better understanding of AMP sequence-activity relationships to improve efficacy in vivo and limit effects of off-target activity. Toward this goal, we experimentally evaluated 210 natural and variant class IIa bacteriocins for antimicrobial activity against six strains of enterococci. Inhibitory activity was ridge regressed to AMP sequence to predict performance, achieving an area under the curve of 0.70 and demonstrating the potential of statistical models for identifying and designing AMPs. Active AMPs were individually produced and evaluated against eight enterococcus strains and four Listeria strains to elucidate trends in susceptibility. It was determined that the mannose phosphotransferase system (manPTS) sequence is informative of susceptibility to class IIa bacteriocins, yet other factors, such as membrane composition, also contribute strongly to susceptibility. A broadly potent bacteriocin variant (lactocin DT1) from a Lactobacillus ruminis genome was identified as the only variant with inhibitory activity toward all tested strains, while a novel enterocin variant (DT2) from an Enterococcus faecium genome demonstrated specificity toward Listeria strains. Eight AMPs were evaluated for proteolytic stability to trypsin, chymotrypsin, and pepsin, and three C-terminal disulfide-containing variants, including divercin V41, were identified as compelling for future in vivo studies, given their high potency and proteolytic stability.IMPORTANCE Class IIa bacteriocin antimicrobial peptides (AMPs), an alternative to traditional small-molecule antibiotics, are capable of selective activity toward various Gram-positive bacteria, limiting negative side effects associated with broad-spectrum activity. This selective activity is achieved through targeting of the mannose phosphotransferase system (manPTS) of a subset of Gram-positive bacteria, although factors affecting this mechanism are not entirely understood. Peptides identified from genomic data, as well as variants of previously characterized AMPs, can offer insight into how peptide sequence affects activity and selectivity. The experimental methods presented here identify promising potent and selective bacteriocins for further evaluation, highlight the potential of simple computational modeling for prediction of AMP performance, and demonstrate that factors beyond manPTS sequence affect bacterial susceptibility to class IIa bacteriocins.

Solubility, Stability, and Avidity of Recombinant Antibody Fragments Expressed in Microorganisms

Solubility of recombinant proteins (i.e., the extent of soluble versus insoluble expression in heterogeneous hosts) is the first checkpoint criterion for determining recombinant protein quality. However, even soluble proteins often fail to represent functional activity because of the involvement of non-functional, misfolded, soluble aggregates, which compromise recombinant protein quality. Therefore, screening of solubility and folding competence is crucial for improving the quality of recombinant proteins, especially for therapeutic applications. The issue is often highlighted especially in bacterial recombinant hosts, since bacterial cytoplasm does not provide an optimal environment for the folding of target proteins of mammalian origin. Antibody fragments, such as single-chain variable fragment (scFv), single-chain antibody (scAb), and fragment antigen binding (Fab), have been utilized for numerous applications such as diagnostics, research reagents, or therapeutics. Antibody fragments can be efficiently expressed in microorganisms so that they offer several advantages for diagnostic applications such as low cost and high yield. However, scFv and scAb fragments have generally lower stability to thermal stress than full-length antibodies, necessitating a judicious combination of designer antibodies, and bacterial hosts harnessed with robust chaperone function. In this review, we discuss efforts on not only the production of antibodies or antibody fragments in microorganisms but also scFv stabilization via (i) directed evolution of variants with increased stability using display systems, (ii) stabilization of the interface between variable regions of heavy (V_H) and light (V_L) chains through the introduction of a non-native covalent bond between the two chains, (iii) rational engineering of V_H-V_L pair, based on the structure, and (iv) computational approaches. We also review recent advances in stability design, increase in avidity by multimerization, and maintaining the functional competence of chimeric proteins prompted by various types of chaperones.

The use of consensus sequence information to engineer stability and activity in proteins

The goal of protein design is to create proteins that are stable, soluble, and active. Here we focus on one approach to protein design in which sequence information is used to create a "consensus" sequence. Such consensus sequences comprise the most common residue at each position in a multiple sequence alignment (MSA). After describing some general ideas that relate MSA and consensus sequences and presenting a statistical thermodynamic framework that relates consensus and non-consensus sequences to stability, we detail the process of designing a consensus sequence and survey reports of consensus design and characterization from the literature. Many of these consensus proteins retain native biological activities including ligand binding and enzyme activity. Remarkably, in most cases the consensus protein shows significantly higher stability than extant versions of the protein, as measured by thermal or chemical denaturation, consistent with the statistical thermodynamic model. To understand this stability increase, we compare various features of consensus sequences with the extant MSA sequences from which they were derived. Consensus sequences show enrichment in charged residues (most notably glutamate and lysine) and depletion of uncharged polar residues (glutamine, serine, and asparagine). Surprisingly, a survey of stability changes resulting from point substitutions show little correlation with residue frequencies at the corresponding positions within the MSA, suggesting that the high stability of consensus proteins may result from interactions among residue pairs or higher-order clusters. Whatever the source, the large number of reported successes demonstrates that consensus design is a viable route to generating active and in many cases highly stabilized proteins.

3D-Printed Phenacrylate Decarboxylase Flow Reactors for the Chemoenzymatic Synthesis of 4-Hydroxystilbene

Continuous flow systems for chemical synthesis are becoming a major focus in organic chemistry and there is a growing interest in the integration of biocatalysts due to their high regio- and stereoselectivity. Methods established for 3D bioprinting enable the fast and simple production of agarose-based modules for biocatalytic reactors if thermally stable enzymes are available. We report here on the characterization of four different cofactor-free phenacrylate decarboxylase enzymes suitable for the production of 4-vinylphenol and test their applicability for the encapsulation and direct 3D printing of disk-shaped agarose-based modules that can be used for compartmentalized flow microreactors. Using the most active and stable phenacrylate decarboxylase from Enterobacter spec. in a setup with four parallel reactors and a subsequent palladium(II) acetate-catalysed Heck reaction, 4-hydroxystilbene was synthesized from p-coumaric acid with a total yield of 14.7 % on a milligram scale. We believe that, due to the convenient direct immobilization of any thermostable enzyme and straightforward tuning of the reaction sequence by stacking of modules with different catalytic activities, this simple process will facilitate the establishment and use of cascade reactions and will therefore be of great advantage for many research approaches.

Improvement in affinity and thermostability of a fully human antibody against interleukin-17A by yeast-display technology and CDR grafting

Monoclonal antibodies (mAbs) are widely used in many fields due to their high specificity and ability to recognize a broad range of antigens. IL-17A can induce a rapid inflammatory response both alone and synergistically with other proinflammatory cytokines. Accumulating evidence suggests that therapeutic intervention of IL-17A signaling offers an attractive treatment option for autoimmune diseases and cancer. Here, we present a combinatorial approach for optimizing the affinity and thermostability of a novel anti-hIL-17A antibody. From a large naïve phage-displayed library, we isolated the anti-IL-17A mAb 7H9 that can neutralize the effects of recombinant human IL-17A. However, the modest neutralization potency and poor thermostability limit its therapeutic applications. In vitro affinity optimization was then used to generate 8D3 by using yeast-displayed random mutagenesis libraries. This resulted in four key amino acid changes and provided an approximately 15-fold potency increase in a cell-based neutralization assay. Complementarity-determining regions (CDRs) of 8D3 were further grafted onto the stable framework of the huFv 4D5 to improve thermostability. The resulting hybrid antibody 9NT/S has superior stabilization and affinities beyond its original antibody. Human fibrosarcoma cell-based assays and in vivo analyses in mice indicated that the anti-IL-17A antibody 9NT/S efficiently inhibited the secretion of IL-17A-induced proinflammatory cytokines. Therefore, this lead anti-IL-17A mAb might be used as a potential best-in-class candidate for treating IL-17A related diseases.

Design and Production of Bispecific Antibodies

With the current biotherapeutic market dominated by antibody molecules, bispecific antibodies represent a key component of the next-generation of antibody therapy. Bispecific antibodies can target two different antigens at the same time, such as simultaneously binding tumor cell receptors and recruiting cytotoxic immune cells. Structural diversity has been fast-growing in the bispecific antibody field, creating a plethora of novel bispecific antibody scaffolds, which provide great functional variety. Two common formats of bispecific antibodies on the market are the single-chain variable fragment (scFv)-based (no Fc fragment) antibody and the full-length IgG-like asymmetric antibody. Unlike the conventional monoclonal antibodies, great production challenges with respect to the quantity, quality, and stability of bispecific antibodies have hampered their wider clinical application and acceptance. In this review, we focus on these two major bispecific types and describe recent advances in the design, production, and quality of these molecules, which will enable this important class of biologics to reach their therapeutic potential.

Consensus designs and thermal stability determinants of a human glutamate transporter

Human excitatory amino acid transporters (EAATs) take up the neurotransmitter glutamate in the brain and are essential to maintain excitatory neurotransmission. Our understanding of the EAATs’ molecular mechanisms has been hampered by the lack of stability of purified protein samples for biophysical analyses. Here, we present approaches based on consensus mutagenesis to obtain thermostable EAAT1 variants that share up to ~95% amino acid identity with the wild type transporters, and remain natively folded and functional. Structural analyses of EAAT1 and the consensus designs using hydrogen-deuterium exchange linked to mass spectrometry show that small and highly cooperative unfolding events at the inter-subunit interface rate-limit their thermal denaturation, while the transport domain unfolds at a later stage in the unfolding pathway. Our findings provide structural insights into the kinetic stability of human glutamate transporters, and introduce general approaches to extend the lifetime of human membrane proteins for biophysical analyses.

An Atypical Mechanism of Split Intein Molecular Recognition and Folding

Split inteins associate to trigger protein splicing in trans, a post-translational modification in which protein sequences fused to the intein pair are ligated together in a traceless manner. Recently, a family of naturally split inteins has been identified that is split at a noncanonical location in the primary sequence. These atypically split inteins show considerable promise in protein engineering applications; however, the mechanism by which they associate is unclear and must be different from that of previously characterized canonically split inteins due to unique topological restrictions. Here, we use a consensus design strategy to generate an atypical split intein pair (Cat) that has greatly improved activity and is amenable to detailed biochemical and biophysical analysis. Guided by the solution structure of Cat, we show that the association of the fragments involves a disorder-to-order structural transition driven by hydrophobic interactions. This molecular recognition mechanism satisfies the topological constraints of the intein fold and, importantly, ensures that premature chemistry does not occur prior to fragment complementation. Our data lead a common blueprint for split intein complementation in which localized structural rearrangements are used to drive folding and regulate protein-splicing activity.

Engineered Charge Redistribution of Gp2 Proteins through Guided Diversity for Improved PET Imaging of Epidermal Growth Factor Receptor

The Gp2 domain is a protein scaffold for synthetic ligand engineering. However, the native protein function results in a heterogeneous distribution of charge on the conserved surface, which may hinder further development and utility. We aim to modulate charge, without diminishing function, which is challenging in small proteins where each mutation is a significant fraction of protein structure. We constructed rationally guided combinatorial libraries with charge-neutralizing or charge-flipping mutations and sorted them, via yeast display and flow cytometry, for stability and target binding. Deep sequencing of functional variants revealed effective mutations both in clone-dependent contexts and broadly across binders to epidermal growth factor receptor (EGFR), insulin receptor, and immunoglobulin G. Functional mutants averaged 4.3 charge neutralizing mutations per domain while maintaining net negative charge. We evolved an EGFR-targeted Gp2 mutant that reduced charge density by 33%, maintained net charge, and improved charge distribution homogeneity while elevating thermal stability ( Tm = 87 ± 1 °C), improving binding specificity, and maintaining affinity ( Kd = 8.8 ± 0.6 nM). This molecule was conjugated with 1,4,7-triazacyclononane,1-glutaric acid-4,7-acetic acid for 64Cu chelation and evaluated for physiological distribution in mice with xenografted A431 (EGFRhigh) and MDA-MB-435 (EGFRlow) tumors. Excised tissue gamma counting and positron emission tomography/computed tomography imaging revealed good EGFRhigh tumor signal (4.7 ± 0.5%ID/g) at 2 h post-injection and molecular specificity evidenced by low uptake in EGFRlow tumors (0.6 ± 0.1%ID/g, significantly lower than for non-charge-modified Gp2, p = 0.01). These results provide charge mutations for an improved Gp2 framework, validate an effective approach to charge engineering, and advance performance of physiological EGFR targeting for molecular imaging.

On-Demand Production of Flow-Reactor Cartridges by 3D Printing of Thermostable Enzymes

The compartmentalization of chemical reactions is an essential principle of life that provides a major source of innovation for the development of novel approaches in biocatalysis. To implement spatially controlled biotransformations, rapid manufacturing methods are needed for the production of biocatalysts that can be applied in flow systems. Whereas three-dimensional (3D) printing techniques offer high-throughput manufacturing capability, they are usually not compatible with the delicate nature of enzymes, which call for physiological processing parameters. We herein demonstrate the utility of thermostable enzymes in the generation of biocatalytic agarose-based inks for a simple temperature-controlled 3D printing process. As examples we utilized an esterase and an alcohol dehydrogenase from thermophilic organisms as well as a decarboxylase that was thermostabilized by directed protein evolution. We used the resulting 3D-printed parts for a continuous, two-step sequential biotransformation in a fluidic setup.

Evaluation of affibody charge modification identified by synthetic consensus design in molecular PET imaging of epidermal growth factor receptor

Tumor overexpression of epidermal growth factor receptor (EGFR) correlates to therapeutic response in select patient populations. Thus, molecular positron emission tomography (PET) imaging of EGFR could stratify responders versus non-responders. We previously demonstrated effectiveness of a "synthetic consensus" design principle to identify six neutralizing mutations within a 58-amino acid EGFR-targeted affibody domain. Herein, we extend the approach to identify additional neutralized variants that vary net charge from -2 to either -4 or +4 while retaining high affinity (1.6 ± 1.2 nM and 2.5 ± 0.7 nM), specific binding to EGFR, secondary structure, and stability (T_m = 68 °C and 59 °C). We radiolabeled the resultant collection of five charge variants with ⁶⁴Cu and evaluated PET imaging performance in murine models with subcutaneously xenografted EGFR^high and EGFR^low tumors. All variants exhibited good EGFR^high tumor imaging as early as 1 h, with EA35S (+3/-5) achieving 7.7 ± 1.4 %ID/g tumor at 4 h with 1.5 ± 0.3%ID/g EGFR^low tumor, 34 ± 5 tumor:muscle and 12 ± 3 tumor:blood ratios. The positively charged EA62S mutant (+6/-2) exhibited 2.2-3.3-fold higher liver signal than the other variants (p<0.01). The EA68 variant with higher charge density was more stable to human and mouse serum than neutralized variants. In a comparison of radiometal chelators, 1,4,7-triazacyclononane,1-glutaric acid-4,7-acetic acid (NODAGA) exhibited superior physiological specificity to 1,4,7,10-tetraazacyclododecane-1,4,7,10-tetraacetic acid (DOTA). In total, these studies comparatively evaluated a set of EGFR-targeted affibodies varying in net charge and charge density, which revealed functional variations that are useful in engineering an ideal probe for translational studies.

Comprehensive reduction of amino acid set in a protein suggests the importance of prebiotic amino acids for stable proteins

Modern organisms commonly use the same set of 20 genetically coded amino acids for protein synthesis with very few exceptions. However, earlier protein synthesis was plausibly much simpler than modern one and utilized only a limited set of amino acids. Nevertheless, few experimental tests of this issue with arbitrarily chosen amino acid sets had been reported prior to this report. Herein we comprehensively and systematically reduced the size of the amino acid set constituting an ancestral nucleoside kinase that was reconstructed in our previous study. We eventually found that two convergent sequences, each comprised of a 13-amino acid alphabet, folded into soluble, stable and catalytically active structures, even though their stabilities and activities were not as high as those of the parent protein. Notably, many but not all of the reduced-set amino acids coincide with those plausibly abundant in primitive Earth. The inconsistent amino acids appeared to be important for catalytic activity but not for stability. Therefore, our findings suggest that the prebiotically abundant amino acids were used for creating stable protein structures and other amino acids with functional side chains were recruited to achieve efficient catalysis.

Characterization of Reconstructed Ancestral Proteins Suggests a Change in Temperature of the Ancient Biosphere

Understanding the evolution of ancestral life, and especially the ability of some organisms to flourish in the variable environments experienced in Earth’s early biosphere, requires knowledge of the characteristics and the environment of these ancestral organisms. Information about early life and environmental conditions has been obtained from fossil records and geological surveys. Recent advances in phylogenetic analysis, and an increasing number of protein sequences available in public databases, have made it possible to infer ancestral protein sequences possessed by ancient organisms. However, the in silico studies that assess the ancestral base content of ribosomal RNAs, the frequency of each amino acid in ancestral proteins, and estimate the environmental temperatures of ancient organisms, show conflicting results. The characterization of ancestral proteins reconstructed in vitro suggests that ancient organisms had very thermally stable proteins, and therefore were thermophilic or hyperthermophilic. Experimental data supports the idea that only thermophilic ancestors survived the catastrophic increase in temperature of the biosphere that was likely associated with meteorite impacts during the early history of Earth. In addition, by expanding the timescale and including more ancestral proteins for reconstruction, it appears as though the Earth’s surface temperature gradually decreased over time, from Archean to present.

Exploiting sequence and stability information for directing nanobody stability engineering

BACKGROUND

Variable domains of camelid heavy-chain antibodies, commonly named nanobodies, have high biotechnological potential. In view of their broad range of applications in research, diagnostics and therapy, engineering their stability is of particular interest. One important aspect is the improvement of thermostability, because it can have immediate effects on conformational stability, protease resistance and aggregation propensity of the protein.

METHODS

We analyzed the sequences and thermostabilities of 78 purified nanobody binders. From this data, potentially stabilizing amino acid variations were identified and studied experimentally.

RESULTS

Some mutations improved the stability of nanobodies by up to 6.1°C, with an average of 2.3°C across eight modified nanobodies. The stabilizing mechanism involves an improvement of both conformational stability and aggregation behavior, explaining the variable degree of stabilization in individual molecules. In some instances, variations predicted to be stabilizing actually led to thermal destabilization of the proteins. The reasons for this contradiction between prediction and experiment were investigated.

CONCLUSIONS

The results reveal a mutational strategy to improve the biophysical behavior of nanobody binders and indicate a species-specificity of nanobody architecture.

GENERAL SIGNIFICANCE

This study illustrates the potential and limitations of engineering nanobody thermostability by merging sequence information with stability data, an aspect that is becoming increasingly important with the recent development of high-throughput biophysical methods.

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.

Protein Stability: Enhancement and Measurement

This article defines protein stability, emphasizes its importance and surveys the field of protein stabilization, with summary reference to a selection of 2009-2015 publications. One can enhance stability by, in particular, protein engineering strategies and by chemical modification (including conjugation) in solution. General protocols are set out on how to measure a given protein's (1) kinetic thermal stability, and (2) oxidative stability, and (3) how to undertake chemical modification of a protein in solution.

Synthetic and natural consensus design for engineering charge within an affibody targeting epidermal growth factor receptor

Protein ligand charge can impact physiological delivery with charge reduction often benefiting performance. Yet neutralizing mutations can be detrimental to protein function. Herein, three approaches are evaluated to introduce charged-to-neutral mutations of three cations and three anions within an affibody engineered to bind epidermal growth factor receptor. These approaches-combinatorial library sorting or consensus design, based on natural homologs or library-sorted mutants-are used to identify mutations with favorable affinity, stability, and recombinant yield. Consensus design, based on 942 affibody homologs, yielded a mutant of modest function (Kd  = 11 ±4 nM, Tm  = 62°C, and yield = 4.0 ± 0.8 mg/L as compared to 5.3 ± 1.7 nM, 71°C, and 3.5 ± 0.3 mg/L for the parental affibody). Extension of consensus design to 10 additional mutants exhibited varied performance including a substantially improved mutant (Kd  = 6.9 ± 1.4 nM, Tm  = 71°C, and 12.7 ± 0.9 mg/L yield). Sorting a homolog-based combinatorial library of 7 × 10(5) mutants generated a distribution of mutants with lower stability and yield, but did identify one strongly binding variant (Kd  = 1.2 ± 0.3 nM, Tm  = 69°C, and 6.0 ± 0.4 mg/L yield). Synthetic consensus design, based on the amino acid distribution in functional library mutants, yielded higher affinities (P = 0.05) with comparable stabilities and yields. The best of four analyzed clones had Kd  = 1.7 ± 0.5 nM, Tm  = 68°C, and 7.0 ± 0.5 mg/L yield. While all three approaches were effective in creating targeted affibodies with six charged-to-neutral mutations, synthetic consensus design proved to be the most robust. Synthetic consensus design provides a valuable tool for ligand engineering, particularly in the context of charge manipulation. Biotechnol. Bioeng. 2016;113: 1628-1638. © 2016 Wiley Periodicals, Inc.

Gene-specific amplicons from metagenomes as an alternative to directed evolution for enzyme screening: a case study using phenylacetaldehyde reductases

Screening gene‐specific amplicons from metagenomes (S‐GAM) is a highly promising technique for the isolation of genes encoding enzymes for biochemical and industrial applications. From metagenomes, we isolated phenylacetaldehyde reductase (par) genes, which code for an enzyme that catalyzes the production of various Prelog's chiral alcohols. Nearly full‐length par genes were amplified by PCR from metagenomic DNA, the products of which were fused with engineered par sequences at both terminal regions of the expression vector to ensure proper expression and then used to construct Escherichia coli plasmid libraries. Sequence‐ and activity‐based screening of these libraries identified different homologous par genes, Hpar‐001 to ‐036, which shared more than 97% amino acid sequence identity with PAR. Comparative characterization of these active homologs revealed a wide variety of enzymatic properties including activity, substrate specificity, and thermal stability. Moreover, amino acid substitutions in these genes coincided with those of Sar268 and Har1 genes, which were independently engineered by error‐prone PCR to exhibit increased activity in the presence of concentrated 2‐propanol. The comparative data from both approaches suggest that sequence information from homologs isolated from metagenomes is quite useful for enzyme engineering. Furthermore, by examining the GAM‐based sequence dataset derived from soil metagenomes, we easily found amino acid substitutions that increase the thermal stability of PAR/PAR homologs. Thus, GAM‐based approaches can provide not only useful homologous enzymes but also an alternative to directed evolution methodologies.

Design of a Split Intein with Exceptional Protein Splicing Activity

Protein trans-splicing (PTS) by split inteins has found widespread use in chemical biology and biotechnology. Herein, we describe the use of a consensus design approach to engineer a split intein with enhanced stability and activity that make it more robust than any known PTS system. Using batch mutagenesis, we first conduct a detailed analysis of the difference in splicing rates between the Npu (fast) and Ssp (slow) split inteins of the DnaE family and find that most impactful residues lie on the second shell of the protein, directly adjacent to the active site. These residues are then used to generate an alignment of 73 naturally occurring DnaE inteins that are predicted to be fast. The consensus sequence from this alignment (Cfa) demonstrates both rapid protein splicing and unprecedented thermal and chaotropic stability. Moreover, when fused to various proteins including antibody heavy chains, the N-terminal fragment of Cfa exhibits increased expression levels relative to other N-intein fusions. The durability and efficiency of Cfa should improve current intein based technologies and may provide a platform for the development of new protein chemistry techniques.

Protein thermostability engineering

Using structure and sequence based analysis we can engineer proteins to increase their thermal stability.

Epistasis effects of multiple ancestral-consensus amino acid substitutions on the thermal stability of glycerol kinase from Cellulomonas sp. NT3060

Thermostable variants of the Cellulomonas sp. NT3060 glycerol kinase have been constructed by through the introduction of ancestral-consensus mutations. We produced seven mutants, each having an ancestral-consensus amino acid residue that might be present in the common ancestors of both bacteria and of archaea, and that appeared most frequently at the position of 17 glycerol kinase sequences in the multiple sequence alignment. The thermal stabilities of the resulting mutants were assessed by determining their melting temperatures (Tm), which was defined as the temperature at which 50% of the initial catalytic activity is lost after 15 min of incubation, as well as when the half-life of the catalytic activity occurs at a temperature of 60°C (t1/2). Three mutants showed increased stabilities compared to the wild-type protein. We then produced five more mutants with multiple amino acid substitutions. Some of the resulting mutants showed thermal stabilities much greater than those expected given the stabilities of the respective mutants with single mutations. Therefore, the effects of mutations are not always simply additive and some amino acid substitutions, which do not affect or only slightly improve stability when individually introduced into the protein, show substantial stabilizing effects in combination with other mutations.

A general approach to antibody thermostabilization

Antibody engineering to enhance thermostability may enable further application and ease of use of antibodies across a number of different areas. A modified human IgG framework has been developed through a combination of engineering approaches, which can be used to stabilize antibodies of diverse specificity. This is achieved through a combination of complementarity-determining region (CDR)-grafting onto the stable framework, mammalian cell display and in vitro somatic hypermutation (SHM). This approach allows both stabilization and maturation to affinities beyond those of the original antibody, as shown by the stabilization of an anti-HA33 antibody by approximately 10°C and affinity maturation of approximately 300-fold over the original antibody. Specificities of 10 antibodies of diverse origin were successfully transferred to the stable framework through CDR-grafting, with 8 of these successfully stabilized, including the therapeutic antibodies adalimumab, stabilized by 9.9°C, denosumab, stabilized by 7°C, cetuximab stabilized by 6.9°C and to a lesser extent trastuzumab stabilized by 0.8°C. This data suggests that this approach may be broadly useful for improving the biophysical characteristics of antibodies across a number of applications.

Fusion to a highly stable consensus albumin binding domain allows for tunable pharmacokinetics

A number of classes of proteins have been engineered for high stability using consensus sequence design methods. Here we describe the engineering of a novel albumin binding domain (ABD) three-helix bundle protein. The resulting engineered ABD molecule, called ABDCon, is expressed at high levels in the soluble fraction of Escherichia coli and is highly stable, with a melting temperature of 81.5°C. ABDCon binds human, monkey and mouse serum albumins with affinity as high as 61 pM. The solution structure of ABDCon is consistent with the three-helix bundle design and epitope mapping studies enabled a precise definition of the albumin binding interface. Fusion of a 10 kDa scaffold protein to ABDCon results in a long terminal half-life of 60 h in mice and 182 h in cynomolgus monkeys. To explore the link between albumin affinity and in vivo exposure, mutations were designed at the albumin binding interface of ABDCon yielding variants that span an 11 000-fold range in affinity. The PK properties of five such variants were determined in mice in order to demonstrate the tunable nature of serum half-life, exposure and clearance with variations in albumin binding affinity.

The consensus-based approach for gene/enzyme replacement therapies and crystallization strategies: the case of human alanine-glyoxylate aminotransferase

Protein stability is a fundamental issue in biomedical and biotechnological applications of proteins. Among these applications, gene- and enzyme-replacement strategies are promising approaches to treat inherited diseases that may benefit from protein engineering techniques, even though these beneficial effects have been largely unexplored. In the present study we apply a sequence-alignment statistics procedure (consensus-based approach) to improve the activity and stability of the human AGT (alanine-glyoxylate aminotransferase) protein, an enzyme which causes PH1 (primary hyperoxaluria type I) upon mutation. By combining only five consensus mutations, we obtain a variant (AGT-RHEAM) with largely enhanced in vitro thermal and kinetic stability, increased activity, and with no side effects on foldability and peroxisomal targeting in mammalian cells. The structure of AGT-RHEAM reveals changes at the dimer interface and improved electrostatic interactions responsible for increased kinetic stability. Consensus-based variants maintained the overall protein fold, crystallized more easily and improved the expression as soluble proteins in two different systems [AGT and CIPK24 (CBL-interacting serine/threonine-protein kinase) SOS2 (salt-overly-sensitive 2)]. Thus the consensus-based approach also emerges as a simple and generic strategy to increase the crystallization success for hard-to-get protein targets as well as to enhance protein stability and function for biomedical applications.

The lower limits for protein stability and foldability in primary hyperoxaluria type I

Mutational effects on protein stability and foldability are important to understand conformational diseases and protein evolution. In this work, we perform a comprehensive investigation on the energetic basis underlying mutational effects on the stability of human alanine:glyoxylate aminotransferase (AGT). We study twenty two variants whose kinetic stabilities span over eleven orders of magnitude and are classified into two groups: i) ten naturally-occurring variants, including the most common mutations causing primary hyperoxaluria type I (PH1); and ii) twelve consensus variants obtained by sequence-alignment statistics. We show that AGT dimer stability determines denaturation rates, and mutations modulate stability by changes in the effective thermodynamic stability, the aggregation propensity of partially/globally unfolded states and subtle energetic changes in the rate-limiting denaturation step. In combination with our previous expression analyses in eukaryotic cells, we propose the existence of two lower limits for AGT stability, one linked to optimal folding efficiency (close to the major allele stability) and the other setting a minimal efficiency compatible with glyoxylate detoxification in vivo (close to the minor allele stability). These lower limits could explain the high prevalence of misfolding as a disease mechanism in PH1 and support the use of pharmacological ligands aimed to increase AGT stability as therapies for this disease.

Structural basis of the regulatory mechanism of the plant CIPK family of protein kinases controlling ion homeostasis and abiotic stress

Plant cells have developed specific protective molecular machinery against environmental stresses. The family of CBL-interacting protein kinases (CIPK) and their interacting activators, the calcium sensors calcineurin B-like (CBLs), work together to decode calcium signals elicited by stress situations. The molecular basis of biological activation of CIPKs relies on the calcium-dependent interaction of a self-inhibitory NAF motif with a particular CBL, the phosphorylation of the activation loop by upstream kinases, and the subsequent phosphorylation of the CBL by the CIPK. We present the crystal structures of the NAF-truncated and pseudophosphorylated kinase domains of CIPK23 and CIPK24/SOS2. In addition, we provide biochemical data showing that although CIPK23 is intrinsically inactive and requires an external stimulation, CIPK24/SOS2 displays basal activity. This data correlates well with the observed conformation of the respective activation loops: Although the loop of CIPK23 is folded into a well-ordered structure that blocks the active site access to substrates, the loop of CIPK24/SOS2 protrudes out of the active site and allows catalysis. These structures together with biochemical and biophysical data show that CIPK kinase activity necessarily requires the coordinated releases of the activation loop from the active site and of the NAF motif from the nucleotide-binding site. Taken all together, we postulate the basis for a conserved calcium-dependent NAF-mediated regulation of CIPKs and a variable regulation by upstream kinases.

Selection of high-affinity Centyrin FN3 domains from a simple library diversified at a combination of strand and loop positions

Alternative scaffold molecules represent a class of proteins important to the study of protein design and mechanisms of protein-protein interactions, as well as for the development of therapeutic proteins. Here, we describe the generation of a library built upon the framework of a consensus FN3 domain sequence resulting in binding proteins we call Centyrins. This new library employs diversified positions within the C-strand, CD-loop, F-strand and FG-loop of the FN3 domain. CIS display was used to select high-affinity Centyrin variants against three targets; c-MET, murine IL-17A and rat TNFα and scanning mutagenesis studies were used to define the positions of the library most important for target binding. Contributions from both the strand and loop positions were noted, although the pattern was different for each molecule. In addition, an affinity maturation scheme is described that resulted in a significant improvement in the affinity of one selected Centyrin variant. Together, this work provides important data contributing to our understanding of potential FN3 binding interfaces and a new tool for generating high-affinity scaffold molecules.

Antibody variable domain interface and framework sequence requirements for stability and function by high-throughput experiments

Protein structural stability and biological functionality are dictated by the formation of intradomain cores and interdomain interfaces, but the intricate sequence-structure-function interrelationships in the packing of protein cores and interfaces remain difficult to elucidate due to the intractability of enumerating all packing possibilities and assessing the consequences of all the variations. In this work, groups of β strand residues of model antibody variable domains were randomized with saturated mutagenesis and the functional variants were selected for high-throughput sequencing and high-throughput thermal stability measurements. The results show that the sequence preferences of the intradomain hydrophobic core residues are strikingly flexible among hydrophobic residues, implying that these residues are coupled indirectly with antigen binding through energetic stabilization of the protein structures. By contrast, the interdomain interface residues are directly coupled with antigen binding. The interdomain interface should be treated as an integral part of the antigen-binding site.

Experimental evidence for the thermophilicity of ancestral life

Theoretical studies have focused on the environmental temperature of the universal common ancestor of life with conflicting conclusions. Here we provide experimental support for the existence of a thermophilic universal common ancestor. We present the thermal stabilities and catalytic efficiencies of nucleoside diphosphate kinases (NDK), designed using the information contained in predictive phylogenetic trees, that seem to represent the last common ancestors of Archaea and of Bacteria. These enzymes display extreme thermal stabilities, suggesting thermophilic ancestries for Archaea and Bacteria. The results are robust to the uncertainties associated with the sequence predictions and to the tree topologies used to infer the ancestral sequences. Moreover, mutagenesis experiments suggest that the universal ancestor also possessed a very thermostable NDK. Because, as we show, the stability of an NDK is directly related to the environmental temperature of its host organism, our results indicate that the last common ancestor of extant life was a thermophile that flourished at a very high temperature.

Thermodynamic stability of domain III from the envelope protein of flaviviruses and its improvement by molecular design

The Flavivirus genus includes widespread and severe human pathogens like the four serotypes of dengue virus (DENV1 to DENV4), yellow fever virus, Japanese encephalitis virus and West Nile virus. Domain III (ED3) of the viral envelope protein interacts with cell receptors and contains epitopes recognized by virus neutralizing antibodies. Its structural, antigenic and immunogenic properties have been thoroughly studied contrary to its physico-chemical properties. Here, the ED3 domains of the above pathogenic flaviviruses were produced in the periplasm of Escherichia coli. Their thermodynamic stabilities were measured and compared in experiments of unfolding equilibriums, induced with chemicals or heat and monitored through protein fluorescence. A designed ED3 domain, with the consensus sequence of DENV strains from all serotypes, was highly stable. The low stability of the ED3 domain from DENV3 was increased by three changes of residues in the protein core without affecting its reactivity towards DENV-infected human serums. Additional changes showed that the stability of ED3 varied with the DENV3 genotype. The T(m) of ED3 was higher than 69°C for all the tested viruses and reached 86°C for the consensus ED3. The latter, deprived of its disulfide bond by mutations, was predominantly unfolded at 20°C. These results will help better understand and design the properties of ED3 for its use as diagnostic, vaccine or therapeutic tools.

Engineering protein stability

This article defines protein stability, emphasizes its importance and surveys some notable recent publications (2004-2008) in the field of protein stability/stabilization. Knowledge of the factors stabilizing proteins has emerged from denaturation studies and from study of thermophilic (and other extremophilic) proteins. One can enhance stability by protein engineering strategies, the judicious use of solutes and additives, immobilization, and chemical modification in solution. General protocols are set out on how to measure the kinetic thermal stability of a given protein and how to undertake chemical modification of a protein in solution.

Generic approach for the generation of stable humanized single-chain Fv fragments from rabbit monoclonal antibodies

Despite their favorable pharmacokinetic properties, single-chain Fv antibody fragments (scFvs) are not commonly used as therapeutics, mainly due to generally low stabilities and poor production yields. In this work, we describe the identification and optimization of a human scFv scaffold, termed FW1.4, which is suitable for humanization and stabilization of a broad variety of rabbit antibody variable domains. A motif consisting of five structurally relevant framework residues that are highly conserved in rabbit variable domains was introduced into FW1.4 to generate a generically applicable scFv scaffold, termed FW1.4gen. Grafting of complementarity determining regions (CDRs) from 15 different rabbit monoclonal antibodies onto FW1.4 and their derivatives resulted in humanized scFvs with binding affinities in the range from 4.7 x 10(-9) to 1.5 x 10(-11) m. Interestingly, minimalistic grafting of CDRs onto FW1.4gen, without any substitutions in the framework regions, resulted in affinities ranging from 5.7 x 10(-10) to <1.8 x 10(-12) m. When compared with progenitor rabbit scFvs, affinities of most humanized scFvs were similar. Moreover, in contrast to progenitor scFvs, which were difficult to produce, biophysical properties of the humanized scFvs were significantly improved, as exemplified by generally good production yields in a generic refolding process and by apparent melting temperatures between 53 and 86 degrees C. Thus, minimalistic grafting of rabbit CDRs on the FW1.4gen scaffold presents a simple and reproducible approach to humanize and stabilize rabbit variable domains.

Conserved amino acid networks involved in antibody variable domain interactions

Engineered antibodies are a large and growing class of protein therapeutics comprising both marketed products and many molecules in clinical trials in various disease indications. We investigated naturally conserved networks of amino acids that support antibody V(H) and V(L) function, with the goal of generating information to assist in the engineering of robust antibody or antibody-like therapeutics. We generated a large and diverse sequence alignment of V-class Ig-folds, of which V(H) and V(L) domains are family members. To identify conserved amino acid networks, covariations between residues at all possible position pairs were quantified as correlation coefficients (phi-values). We provide rosters of the key conserved amino acid pairs in antibody V(H) and V(L) domains, for reference and use by the antibody research community. The majority of the most strongly conserved amino acid pairs in V(H) and V(L) are at or adjacent to the V(H)-V(L) interface suggesting that the ability to heterodimerize is a constraining feature of antibody evolution. For the V(H) domain, but not the V(L) domain, residue pairs at the variable-constant domain interface (V(H)-C(H)1 interface) are also strongly conserved. The same network of conserved V(H) positions involved in interactions with both the V(L) and C(H)1 domains is found in camelid V(HH) domains, which have evolved to lack interactions with V(L) and C(H)1 domains in their mature structures; however, the amino acids at these positions are different, reflecting their different function. Overall, the data describe naturally occurring amino acid networks in antibody Fv regions that can be referenced when designing antibodies or antibody-like fragments with the goal of improving their biophysical properties.

Inferring stabilizing mutations from protein phylogenies: application to influenza hemagglutinin

One selection pressure shaping sequence evolution is the requirement that a protein fold with sufficient stability to perform its biological functions. We present a conceptual framework that explains how this requirement causes the probability that a particular amino acid mutation is fixed during evolution to depend on its effect on protein stability. We mathematically formalize this framework to develop a Bayesian approach for inferring the stability effects of individual mutations from homologous protein sequences of known phylogeny. This approach is able to predict published experimentally measured mutational stability effects (ΔΔG values) with an accuracy that exceeds both a state-of-the-art physicochemical modeling program and the sequence-based consensus approach. As a further test, we use our phylogenetic inference approach to predict stabilizing mutations to influenza hemagglutinin. We introduce these mutations into a temperature-sensitive influenza virus with a defect in its hemagglutinin gene and experimentally demonstrate that some of the mutations allow the virus to grow at higher temperatures. Our work therefore describes a powerful new approach for predicting stabilizing mutations that can be successfully applied even to large, complex proteins such as hemagglutinin. This approach also makes a mathematical link between phylogenetics and experimentally measurable protein properties, potentially paving the way for more accurate analyses of molecular evolution.

Mutating a protein frequently causes a change in its stability. As scientists, we often care about these changes because we would like to engineer a protein's stability or understand how its stability is impacted by a naturally occurring mutation. Evolution also cares about mutational stability changes, because a basic evolutionary requirement is that proteins remain sufficiently stable to perform their biological functions. Our work is based on the idea that it should be possible to use the fact that evolution selects for stability to infer from related proteins the effects of specific mutations. We show that we can indeed use protein evolutionary histories to computationally predict previously measured mutational stability changes more accurately than methods based on either of the two main existing strategies. We then test whether we can predict mutations that increase the stability of hemagglutinin, an influenza protein whose rapid evolution is partly responsible for the ability of this virus to cause yearly epidemics. We experimentally create viruses carrying predicted stabilizing mutations and find that several do in fact improve the virus's ability to grow at higher temperatures. Our computational approach may therefore be of use in understanding the evolution of this medically important virus.