Design of highly stable functional GroEL minichaperones

Design of stable circular permutants of the GroEL chaperone apical domain

Enhancing protein stability holds paramount significance in biotechnology, therapeutics, and the food industry. Circular permutations offer a distinctive avenue for manipulating protein stability while keeping intra-protein interactions intact. Amidst the creation of circular permutants, determining the optimal placement of the new N- and C-termini stands as a pivotal, albeit largely unexplored, endeavor. In this study, we employed PONDR-FIT's predictions of disorder propensity to guide the design of circular permutants for the GroEL apical domain (residues 191-345). Our underlying hypothesis posited that a higher predicted disorder value would correspond to reduced stability in the circular permutants, owing to the increased likelihood of fluctuations in the novel N- and C-termini. To substantiate this hypothesis, we engineered six circular permutants, positioning glycines within the loops as locations for the new N- and C-termini. We demonstrated the validity of our hypothesis along the set of the designed circular permutants, as supported by measurements of melting temperatures by circular dichroism and differential scanning microcalorimetry. Consequently, we propose a novel computational methodology that rationalizes the design of circular permutants with projected stability. Video Abstract.

GroEL—A Versatile Chaperone for Engineering and a Plethora of Applications

Chaperones play a vital role in the life of cells by facilitating the correct folding of other proteins and maintaining them in a functional state, being themselves, as a rule, more stable than the rest of cell proteins. Their functional properties naturally tempt investigators to actively adapt them for biotechnology needs. This review will mostly focus on the applications found for the bacterial chaperonin GroE and its counterparts from other organisms, in biotechnology or for research purposes, both in their engineered or intact versions.

Surface residues and non-additive interactions stabilize a consensus homeodomain protein

Despite the widely reported success of consensus design in producing highly stabilized proteins, little is known about the physical mechanisms underlying this stabilization. Here we explore the potential sources of stabilization by performing a systematic analysis of the 29 substitutions that we previously found to collectively stabilize a consensus homeodomain compared to an extant homeodomain. By separately introducing groups of consensus substitutions that alter or preserve charge state, occur at varying degrees of residue burial, and occur at positions of varying degrees of conservation, we determine the extent to which these three features contribute to the consensus stability enhancement. Surprisingly, we find that the largest total contribution to stability comes from consensus substitutions on the protein surface and that the largest per-substitution contributions come from substitutions that maintain charge state. This finding suggests that although consensus proteins are often enriched in charged residues, consensus stabilization does not result primarily from interactions involving charged residues. Although consensus substitutions at strongly conserved positions also contribute disproportionately to stabilization, significant stabilization is also contributed from substitutions at weakly conserved positions. Furthermore, we find that identical consensus substitutions show larger stabilizing effects when introduced into the consensus background than when introduced into an extant homeodomain, indicating that synergistic, stabilizing interactions among the consensus residues contribute to consensus stability enhancement of the homeodomain. By measuring DNA binding affinity for the same set of variants, we find that although consensus design of the homeodomain increases both affinity and folding stability, it does so using a largely non-overlapping set of substitutions.

Consensus mutagenesis approach improves the thermal stability of system xc - transporter, xCT, and enables cryo-EM analyses

System x_c⁻ is an amino acid antiporter that imports L‐cystine into cells and exports intracellular L‐glutamate, at a 1:1 ratio. As L‐cystine is an essential precursor for glutathione synthesis, system x_c⁻ supports tumor cell growth through glutathione‐based oxidative stress resistance and is considered as a potential therapeutic target for cancer treatment. System x_c⁻ consists of two subunits, the light chain subunit SLC7A11 (xCT) and the heavy chain subunit SLC3A2 (also known as CD98hc or 4F2hc), which are linked by a conserved disulfide bridge. Although the recent structures of another SLC7 member, L‐type amino acid transporter 1 (LAT1) in complex with CD98hc, have provided the structural basis toward understanding the amino acid transport mechanism, the detailed molecular mechanism of xCT remains unknown. To revealthe molecular mechanism, we performed single‐particle analyses of the xCT‐CD98hc complex. As wild‐type xCT‐CD98hc displayed poor stability and could not be purified to homogeneity, we applied a consensus mutagenesis approach to xCT. The consensus mutated construct exhibited increased stability as compared to the wild‐type, and enabled the cryoelectron microscopy (cryo‐EM) map to be obtained at 6.2 Å resolution by single‐particle analysis. The cryo‐EM map revealed sufficient electron density to assign secondary structures. In the xCT structure, the hash and arm domains are well resolved, whereas the bundle domain shows some flexibility. CD98hc is positioned next to the xCT transmembrane domain. This study provides the structural basis of xCT, and our consensus‐based strategy could represent a good choice toward solving unstable protein structures.

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.

Consensus protein engineering on the thermostable histone-like bacterial protein HUs significantly improves stability and DNA binding affinity

Consensus-based protein engineering strategy has been applied to various proteins and it can lead to the design of proteins with enhanced biological performance. Histone-like HUs comprise a protein family with sequence variety within a highly conserved 3D-fold. HU function includes compacting and regulating bacterial DNA in a wide range of biological conditions in bacteria. To explore the possible impact of consensus-based design in the thermodynamic stability of HU proteins, the approach was applied using a dataset of sequences derived from a group of 40 mesostable, thermostable, and hyperthermostable HUs. The consensus-derived HU protein was named HUBest, since it is expected to perform best. The synthetic HU gene was overexpressed in E. coli and the recombinant protein was purified. Subsequently, HUBest was characterized concerning its correct folding and thermodynamic stability, as well as its ability to interact with plasmid DNA. A substantial increase in HUBest stability at high temperatures is observed. HUBest has significantly improved biological performance at ambience temperature, presenting very low K_d values for binding plasmid DNA as indicated from the Gibbs energy profile of HUBest. This K_d may be associated to conformational changes leading to decreased thermodynamic stability and, therefore, higher flexibility at ambient temperature.

Strategies for Increasing Protein Stability

The stability of wild-type proteins is often a hurdle to their practical use in research, industry, and medicine. The route to engineering stability of a protein of interest lies largely with the available data. Where high-resolution structural data is available, rational design, based on fundamental principles of protein chemistry, can improve protein stability. Recent advances in computational biology and the use of nonnatural amino acids have also provided novel rational methods for improving protein stability. Likewise, the explosion of sequence and structural data available in public databases, in combination with improvements in freely available computational tools, has produced accessible phylogenetic approaches. Trawling modern sequence databases can identify the thermostable homologs of a target protein, and evolutionary data can be quickly generated using available phylogenetic tools. Grafting features from those thermostable homologs or ancestors provides stability improvement through a semi-rational approach. Further, molecular techniques such as directed evolution have shown great promise in delivering designer proteins. These strategies are well documented and newly accessible to the molecular biologist, allowing for rapid enhancements of protein stability.

Methionine in a protein hydrophobic core drives tight interactions required for assembly of spider silk

Web spiders connect silk proteins, so-called spidroins, into fibers of extraordinary toughness. The spidroin N-terminal domain (NTD) plays a pivotal role in this process: it polymerizes spidroins through a complex mechanism of dimerization. Here we analyze sequences of spidroin NTDs and find an unusually high content of the amino acid methionine. We simultaneously mutate all methionines present in the hydrophobic core of a spidroin NTD from a nursery web spider’s dragline silk to leucine. The mutated NTD is strongly stabilized and folds at the theoretical speed limit. The structure of the mutant is preserved, yet its ability to dimerize is substantially impaired. We find that side chains of core methionines serve to mobilize the fold, which can thereby access various conformations and adapt the association interface for tight binding. Methionine in a hydrophobic core equips a protein with the capacity to dynamically change shape and thus to optimize its function.

Spider silk is of interest in material science research. Here the authors show that the tight binding of a spider silk protein domain relies on the amino acid methionine, which is abundant in the domain core where it facilitates dynamic shape adaption of the binding interface.

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.

[Prediction of protein thermostability from their primary structure: the current state and development factors]

The construction of proteins and peptides with desired properties, including resistance to high temperatures, as well as optimization of their amino acid composition, is an important and complex task, which attracts much attention in various branches of the basic sciences, and also in biomedicine and biotechnology. This raises the question: what method is more relevant for the at the pilot stage of research in order to estimate the influence of the planned amino acid substitutions on the thermostability of the resultant protein construct? In this brief review we have classified existing basic practical and theoretical approaches used in studies and predicting the thermal stability of native and recombinant polypeptides. Particular attention has been paid to the predictive potential of statistical methods for studying the thermodynamic parameters of the primary protein structure and prospects of their use.

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.

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.

Structural and dynamic properties that govern the stability of an engineered fibronectin type III domain

Consensus protein design is a rapid and reliable technique for the improvement of protein stability, which relies on the use of homologous protein sequences. To enhance the stability of a fibronectin type III (FN3) domain, consensus design was employed using an alignment of 2123 sequences. The resulting FN3 domain, FN3con, has unprecedented stability, with a melting temperature >100°C, a ΔG_D−N of 15.5 kcal mol⁻¹ and a greatly reduced unfolding rate compared with wild-type. To determine the underlying molecular basis for stability, an X-ray crystal structure of FN3con was determined to 2.0 Å and compared with other FN3 domains of varying stabilities. The structure of FN3con reveals significantly increased salt bridge interactions that are cooperatively networked, and a highly optimized hydrophobic core. Molecular dynamics simulations of FN3con and comparison structures show the cooperative power of electrostatic and hydrophobic networks in improving FN3con stability. Taken together, our data reveal that FN3con stability does not result from a single mechanism, but rather the combination of several features and the removal of non-conserved, unfavorable interactions. The large number of sequences employed in this study has most likely enhanced the robustness of the consensus design, which is now possible due to the increased sequence availability in the post-genomic era. These studies increase our knowledge of the molecular mechanisms that govern stability and demonstrate the rising potential for enhancing stability via the consensus method.

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.

Experimental evolution of a green fluorescent protein composed of 19 unique amino acids without tryptophan

At some stage of evolution, genes of organisms may have encoded proteins that were synthesized using fewer than 20 unique amino acids. Similar to evolution of the natural 19-amino-acid proteins GroEL/ES, proteins composed of 19 unique amino acids would have been able to evolve by accumulating beneficial mutations within the 19-amino-acid repertoire encoded in an ancestral genetic code. Because Trp is thought to be the last amino acid included in the canonical 20-amino-acid repertoire, this late stage of protein evolution could be mimicked by experimental evolution of 19-amino-acid proteins without tryptophan (Trp). To further understand the evolution of proteins, we tried to mimic the evolution of a 19-amino-acid protein involving the accumulation of beneficial mutations using directed evolution by random mutagenesis on the whole targeted gene sequence. We created active 19-amino-acid green fluorescent proteins (GFPs) without Trp from a poorly fluorescent 19-amino-acid mutant, S1-W57F, by using directed evolution with two rounds of mutagenesis and selection. The N105I and S205T mutations showed beneficial effects on the S1-W57F mutant. When these two mutations were combined on S1-W57F, we observed an additive effect on the fluorescence intensity. In contrast, these mutations showed no clear improvement individually or in combination on GFPS1, which is the parental GFP mutant composed of 20 amino acids. Our results provide an additional example for the experimental evolution of 19-amino-acid proteins without Trp, and would help understand the mechanisms underlying the evolution of 19-amino-acid proteins. (236 words).

Modulating protein stability - directed evolution strategies for improved protein function

Protein engineering is widely used to generate proteins with novel or enhanced function. However, manipulating protein function in the laboratory can prove laborious, protracted and challenging. Recent developments in the understanding of protein evolutionary dynamics have unveiled the full extent by which the evolution of function is limited by protein stability - a revelation that may be applied to protein engineering on a whole. Thus, strategies that modulate protein stability and reduce its constraining effects may facilitate the engineering of protein function. A combinatorial approach involving the introduction of compensatory mutations and manipulation of the stability threshold by chaperone buffering during directed evolution can improve the functional adaptation of a protein, thereby fostering our ability to attain ever-more ambitious protein functions in the laboratory.

Highly thermostable fungal cellobiohydrolase I (Cel7A) engineered using predictive methods

Building on our previous efforts to generate thermostable chimeric fungal cellobiohydrolase I (CBH I, also known as Cel7A) cellulases by structure-guided recombination, we used FoldX and a 'consensus' sequence approach to identify individual mutations present in the five homologous parent CBH I enzymes which further stabilize the chimeras. Using the FoldX force field, we calculated the effect on ΔG(Folding) of each candidate mutation in a number of CBH I structures and chose those predicted to be stabilizing in multiple structures. With an alignment of 41 CBH I sequences, we also used amino acid frequencies at each candidate position to calculate predicted effects on ΔG(Folding). A combination of mutations chosen using these methods increased the T(50) of the most thermostable chimera by an additional 4.7°C, to yield a CBH I with T(50) of 72.1°C, which is 9.2°C higher than that of the most stable native CBH I, from Talaromyces emersonii. This increased stability resulted in a 10°C increase in the optimal temperature for activity, to 65°C, and a 50% increase in total sugar production from crystalline cellulose at the optimal temperature, compared with native T.emersonii CBH I.

Design of novel FN3 domains with high stability by a consensus sequence approach

The use of consensus design to produce stable proteins has been applied to numerous structures and classes of proteins. Here, we describe the engineering of novel FN3 domains from two different proteins, namely human fibronectin and human tenascin-C, as potential alternative scaffold biotherapeutics. The resulting FN3 domains were found to be robustly expressed in Escherichia coli, soluble and highly stable, with melting temperatures of 89 and 78°C, respectively. X-ray crystallography was used to confirm that the consensus approach led to a structure consistent with the FN3 design despite having only low-sequence identity to natural FN3 domains. The ability of the Tenascin consensus domain to withstand mutations in the loop regions connecting the β-strands was investigated using alanine scanning mutagenesis demonstrating the potential for randomization in these regions. Finally, rational design was used to produce point mutations that significantly increase the stability of one of the consensus domains. Together our data suggest that consensus FN3 domains have potential utility as alternative scaffold therapeutics.

Genomic and structural aspects of protein evolution

It has been known for more than 35 years that, during evolution, new proteins are formed by gene duplications, sequence and structural divergence and, in many cases, gene combinations. The genome projects have produced complete, or almost complete, descriptions of the protein repertoires of over 600 distinct organisms. Analyses of these data have dramatically increased our understanding of the formation of new proteins. At the present time, we can accurately trace the evolutionary relationships of about half the proteins found in most genomes, and it is these proteins that we discuss in the present review. Usually, the units of evolution are protein domains that are duplicated, diverge and form combinations. Small proteins contain one domain, and large proteins contain combinations of two or more domains. Domains descended from a common ancestor are clustered into superfamilies. In most genomes, the net growth of superfamily members means that more than 90% of domains are duplicates. In a section on domain duplications, we discuss the number of currently known superfamilies, their size and distribution, and superfamily expansions related to biological complexity and to specific lineages. In a section on divergence, we describe how sequences and structures diverge, the changes in stability produced by acceptable mutations, and the nature of functional divergence and selection. In a section on domain combinations, we discuss their general nature, the sequential order of domains, how combinations modify function, and the extraordinary variety of the domain combinations found in different genomes. We conclude with a brief note on other forms of protein evolution and speculations of the origins of the duplication, divergence and combination processes.

Directed evolution of a histone acetyltransferase--enhancing thermostability, whilst maintaining catalytic activity and substrate specificity

Histone acetylation plays an integral role in the epigenetic regulation of gene expression. Transcriptional activity reflects the recruitment of opposing classes of enzymes to promoter elements; histone acetyltransferases (EC 2.3.1.48) that deposit acetyl marks at a subset of histone residues and histone deacetylases that remove them. Many histone acetyltransferases are difficult to study in solution because of their limited stability once purified. We have developed a directed evolution protocol that allows the screening of hundreds of histone acetyltransferase mutants for histone acetylating activity, and used this to enhance the thermostability of the human P/CAF histone acetyltransferase. Two rounds of directed evolution significantly stabilized the enzyme without lowering the catalytic efficiency and substrate specificity of the enzyme. Twenty-four variants with higher thermostability were identified. Detailed analysis revealed twelve single amino acid mutants that were found to possess a higher thermostability. The residues affected are scattered over the entire protein structure, and are different from mutations predicted by sequence alignment approaches, suggesting that sequence comparison and directed evolution methods are complementary strategies in engineering increased protein thermostability. The stabilizing mutations are predominately located at surface of the enzyme, suggesting that the protein's surface is important for stability. The directed evolution approach described in the present study is easily adapted to other histone modifying enzymes, requiring only appropriate peptide substrates and antibodies, which are available from commercial suppliers.

Bioinformatic method for protein thermal stabilization by structural entropy optimization

Engineering proteins for higher thermal stability is an important and difficult challenge. We describe a bioinformatic method incorporating sequence alignments to redesign proteins to be more stable through optimization of local structural entropy. Using this method, improved configurational entropy (ICE), we were able to design more stable variants of a mesophilic adenylate kinase with only the sequence information of one psychrophilic homologue. The redesigned proteins display considerable increases in their thermal stabilities while still retaining catalytic activity. ICE does not require a three-dimensional structure or a large number of homologous sequences, indicating a broad applicability of this method. Our results also highlight the importance of entropy in the stability of protein structures.

Intense neutral drifts yield robust and evolvable consensus proteins

What changes occur when a natural protein that had been under low mutation rates is subjected to a neutral drift at high mutational loads, thus generating genetically diverse (polymorphic) gene ensembles that all maintain the protein's original function and structure? To address this question we subjected large populations of TEM-1 beta-lactamase to a prolonged neutral drift, applying high mutation rates and purifying selection to maintain TEM-1's existing penicillinase activity. Purging of deleterious mutations and enrichment of beneficial ones maintained the sequence of these ensembles closer to TEM-1's family consensus and inferred ancestor. In particular, back-to-consensus/ancestor mutations that increase TEM-1's kinetic and thermodynamic stability were enriched. These acted as global suppressors and enabled the tolerance of a broad range of deleterious mutations, thus further increasing the genetic diversity of the drifting populations. The probability of a new function emerging (cefotaxime degradation) was also substantially increased in these ensembles owing to the presence of many gene variants carrying the global suppressors. Our findings indicate the unique features of large, polymorphic neutral ensembles generated under high mutational loads and prompt the speculation that the progenitors of today's proteins may have evolved under high mutational loads. The results also suggest that predictable back-to-consensus/ancestor changes can be used in the laboratory to generate highly diverse and evolvable gene libraries.

Rational design of a novel calcium-binding site adjacent to the ligand-binding site on CD2 increases its CD48 affinity

Electrostatic interactions are important for molecular recognition processes including Ca2+-binding and cell adhesion. To understand these processes, we have successfully introduced a novel Ca2+-binding site into the non-Ca2+-dependent cell adhesion protein CD2 using our criteria that are specifically tailored to the structural and functional properties of the protein environment and charged adhesion surface. This designed site with ligand residues exclusively from the beta-sheets selectively binds to Ca2+ and Ln3+ over other mono- and divalent cations. While Ca2+ and Ln3+ binding specifically alters the local environment of the designed Ca2+-binding site, the designed protein undergoes a significantly smaller conformation change compared with those observed in naturally occurring Ca2+-binding sites that are composed of at least part of the flexible loop and helical regions. In addition, the CD2-CD48-binding affinity increased approximately threefold after protein engineering, suggesting that the cell adhesion of CD2 can be modulated by altering the local electrostatic environment. The study provides site-specific information for regulating cell adhesion within CD2 and gives insight into the structural factors required for Ca2+-modulated biological processes.

Optimized variants of the cold shock protein from in vitro selection: structural basis of their high thermostability

The bacterial cold shock proteins (Csp) are widely used as models for the experimental and computational analysis of protein stability. In a previous study, in vitro evolution was employed to identify strongly stabilizing mutations in Bs-CspB from Bacillus subtilis. The best variant found by this approach contained the mutations M1R, E3K and K65I, which raised the midpoint of thermal unfolding of Bs-CspB from 53.8 degrees C to 83.7 degrees C, and increased the Gibbs free energy of stabilization by 20.9 kJ mol(-1). Another selected variant with the two mutations A46K and S48R was stabilized by 11.1 kJ mol(-1). To elucidate the molecular basis of these stabilizations, we determined the crystal structures of these two Bs-CspB variants. The mutated residues are generally well ordered and provide additional stabilizing interactions, such as charge interactions, additional hydrogen bonds and improved side-chain packing. Several mutations improve the electrostatic interactions, either by the removal of unfavorable charges (E3K) or by compensating their destabilizing interactions (A46K, S48R). The stabilizing mutations are clustered at a contiguous surface area of Bs-CspB, which apparently is critically important for the stability of the beta-barrel structure but not well optimized in the wild-type protein.

The creation of a novel fluorescent protein by guided consensus engineering

Consensus engineering has been used to increase the stability of a number of different proteins, either by creating consensus proteins from scratch or by modifying existing proteins so that their sequences more closely match a consensus sequence. In this paper we describe the first application of consensus engineering to the ab initio creation of a novel fluorescent protein. This was based on the alignment of 31 fluorescent proteins with >62% homology to monomeric Azami green (mAG) protein, and used the sequence of mAG to guide amino acid selection at positions of ambiguity. This consensus green protein is extremely well expressed, monomeric and fluorescent with red shifted absorption and emission characteristics compared to mAG. Although slightly less stable than mAG, it is better expressed and brighter under the excitation conditions typically used in single molecule fluorescence spectroscopy or confocal microscopy. This study illustrates the power of consensus engineering to create stable proteins using the subtle information embedded in the alignment of similar proteins and shows that the benefits of this approach may extend beyond stability.

Point mutations in protein globular domains: contributions from function, stability and misfolding

Several contrasting hypotheses have been formulated about the influence of functional and conformational properties, like stability and avoidance of misfolding, on the evolution of protein globular domains. Selection at functional sites has been suggested to be detrimental to stability or coupled to it. Avoidance of misfolding may be achieved by discarding misfolding-prone sequences or by maintaining a stable native state and thus destabilizing partially or fully unfolded states from which misfolding can take place. We have performed a hierarchical analysis of a large database of point mutations to dissect the relative contributions of function, stability and misfolding in the evolution of natural sequences. We show that at catalytic sites, selection for function overrules selection for stability but find no evidence for an anticorrelation between function and stability. Selection for stability plays a secondary role at binding sites, but is not fully coupled to selection for function. Remarkably, we did not find a selective pressure against misfolding-prone sequences in globular proteins at the level of individual positions. We suggest that such a selection would compromise native-state stability due to a correlation between the stabilities of native and misfolded states. Stabilization of the native state is the most frequent way in which natural proteins avoid misfolding.

Structure-guided consensus approach to create a more thermostable penicillin G acylase

The thermostabilization of penicillin G acylase (PGA) is a difficult problem due to the large size of the protein and its complex maturation process. We developed a data-driven protein design method that requires fewer homologous sequences than the traditional consensus approach and utilizes structural information to limit the number of variants created. Approximately 50% of our 21 single-point mutants were found experimentally to be more thermostable than the wild-type PGA, two had almost threefold longer half-life at 50 degrees C, with very little effect on activity. An analysis of four programs that predict the thermostability conferred by point mutations shows little agreement between the programs and with the experimental data, emphasizing that the chosen stabilizing mutations are very difficult to predict, but that our data-driven design method should prove useful.

Designing thermostable proteins: ancestral mutants of 3-isopropylmalate dehydrogenase designed by using a phylogenetic tree

We have recently developed a new method for designing thermostable proteins using phylogenetic trees of enzymes. In this study, we investigated a method for designing proteins with improved stability using 3-isopropylmalate dehydrogenase (IPMDH) from Thermus thermophilus as a model enzyme. We designed 12 mutant enzymes, each having an ancestral amino acid residue that was present in the common ancestor of Bacteria and Archaea. At least six of the 12 ancestral mutants tested showed thermal stability higher than that of the original enzyme. The results supported the hyperthermophilic universal ancestor hypothesis. The effect of ancestral residues on IPMDHs of several organisms and on the related enzyme isocitrate dehydrogenase was summarised and analysed. The effect of an ancestral residue on thermostability did not depend on the degree of conservation of the residue at the site, suggesting that the stabilisation of these mutant proteins is not related to sequence conservation but to the antiquity of the introduced residues. The results suggest also that this method could be an efficient way of designing mutant enzymes with higher thermostability based only on the primary structure and a phylogenetic tree.

A modified consensus approach to mutagenesis inverts the cofactor specificity of Bacillus stearothermophilus lactate dehydrogenase

Lactate dehydrogenase from Bacillus stearothermophilus is specific for NAD+. There have been several attempts to alter the cofactor specificity of this enzyme, but these have yielded enzymes with relatively low activities that still largely prefer NAD+. A modified consensus approach was used to create a library of phylogenetically preferred amino acids situated near the cofactor binding site, and variants were screened for their ability to utilize NMN+. A triple mutant (Mut31) was discovered that proved to be more catalytically efficient than wild-type. Mut31 was also better at utilizing NAD+ than the wild-type enzyme and was weakly active with NADP+ and NMN+. An analysis of single amino acid substitutions suggested that all three mutations worked in a concerted fashion to yield robust cofactor utilization. When two previously identified amino acid substitutions were introduced into the Mut31 background, the resultant quintuply substituted enzyme not only utilized NADP+ far better than the wild-type enzyme, it actually inverted its preference for NAD+ and NADP+.

Do Proteins Always Benefit from a Stability Increase? Relevant and Residual Stabilisation in a Three-state Protein by Charge Optimisation

The vast majority of our knowledge on protein stability arises from the study of simple two-state models. However, proteins displaying equilibrium intermediates under certain conditions abound and it is unclear whether the energetics of native/intermediate equilibria is well represented in current knowledge. We consider here that the overall conformational stability of three-state proteins is made of a "relevant" term and a "residual" one, corresponding to the free energy differences of the native to intermediate (N-to-I) and intermediate to denatured (I-to-D) equilibria, respectively. The N-to-I free energy difference is considered to be the relevant stability because protein-unfolding intermediates are likely devoid of biological activity. We use surface charge optimisation to first increase the overall (N-to-D) stability of a model three-state protein (apoflavodoxin) and then investigate whether the stabilisation obtained is realised into relevant or into residual stability. Most of the mutations designed from electrostatic calculations or from simple sequence conservation analysis produce large increases in the overall stability of the protein. However, in most cases, this simply leads to similarly large increases of the residual stability. Two mutations, nevertheless, show a different trend and increase the relevant stability of the protein substantially. When all the mutations are mapped onto the structure of the apoflavodoxin thermal-unfolding intermediate (obtained independently by equilibrium phi-analysis and NMR) they cluster perfectly so that the mutations increasing the relevant stability appear in the small unstructured region of the intermediate and the others in the native-like region. This illustrates the need for specific investigation of N-to-I equilibria and the structure of protein intermediates, and indicates that it is possible to rationally stabilise a protein against partial unfolding once the structure of the intermediate conformation is known, even if at low resolution.

Consensus design of repeat proteins

Consensus design is a valuable protein-engineering method that is based on statistical information derived from sequence alignments of homologous proteins. Recently, consensus design was adapted to repeat proteins. We discuss the potential of this novel repeat-based approach for the design of consensus repeat proteins and repeat protein libraries and summarize recent results from such experiments.

Crystal Structure of the Native Chaperonin Complex from Thermus thermophilus Revealed Unexpected Asymmetry at the cis-Cavity

The chaperonins GroEL and GroES are essential mediators of protein folding. GroEL binds nonnative protein, ATP, and GroES, generating a ternary complex in which protein folding occurs within the cavity capped by GroES (cis-cavity). We determined the crystal structure of the native GroEL-GroES-ADP homolog from Thermus thermophilus, with substrate proteins in the cis-cavity, at 2.8 A resolution. Twenty-four in vivo substrate proteins within the cis-cavity were identified from the crystals. The structure around the cis-cavity, which encapsulates substrate proteins, shows significant differences from that observed for the substrate-free Escherichia coli GroEL-GroES complex. The apical domain around the cis-cavity of the Thermus GroEL-GroES complex exhibits a large deviation from the 7-fold symmetry. As a result, the GroEL-GroES interface differs considerably from the previously reported E. coli GroEL-GroES complex, including a previously unknown contact between GroEL and GroES.

Optimization of the Antibody CH3 Domain by Residue Frequency Analysis of IgG Sequences

In an attempt to enhance the overall assembly, yield and half-life of recombinant antibody proteins, we have cloned and expressed several IgG1 C(H)3 domains and examined their folding/refolding characteristics. We utilized a cytoplasmic bacterial expression system with a thioredoxin reductase knock-out strain of BL21(DE3) to produce bovine, murine and human C(H)3. Under identical conditions, expression of bovine C(H)3 resulted consistently in the highest yields of properly folded/oxidized protein. Circular dichroism and fluorescence experiments demonstrate that oxidized bovine and murine C(H)3 have surprisingly similar structures and stabilities, considering the marginal sequence conservation between the two molecules. Residue frequency analysis using a limited data set of 36 unique Fc sequences originating from 19 different mammalian species targeted five specific sites for optimization within bovine C(H)3. Combination of three of these mutants increased the thermal stability of the molecule to 86 degrees C. Comparison of this approach to similar studies using larger sequence databases and/or different selection criteria suggests sequence database design can increase the success rate for identifying residue sites worth optimizing. This optimized C(H)3 domain can be used as a particularly stable platform for functional design and can be grafted into full-length antibody sequences to enhance their thermodynamic parameters and shelf-life.

Design and characterization of a hyperstable p16INK4a that restores Cdk4 binding activity when combined with oncogenic mutations

Cyclin-dependent kinase inhibitor p16(INK4a) is the founding member of the INK4 family of tumor suppressors capable of arresting mammalian cell division. Missense mutations in the p16(INK4a) gene (INK4a/CDKN2A/MTS1) are strongly linked to several types of human cancer. These mutations are evenly distributed throughout this small, ankyrin repeat protein and the majority of them disrupt the native secondary and/or tertiary structure, leading to protein unfolding, aggregation and loss of function. We report here the use of multiple stabilizing substitutions to increase the stability of p16(INK4a) and furthermore, to restore Cdk4 binding activity of several defective, cancer-related mutant proteins. Stabilizing substitutions were predicted using four different techniques. The three most effective substitutions were combined to create a hyperstable p16(INK4a) variant that is 1.4 kcal/mol more stable than wild-type. This engineered construct is monomeric in solution with wild-type-like secondary and tertiary structure and cyclin-dependent kinase 4 binding activity. Interestingly, these hyperstable substitutions, when combined with oncogenic mutations R24P, P81L or V126D, can significantly restore Cdk4 binding activity, despite the divergent features of each destabilizing mutation. Extensive biophysical studies indicate that the hyperstable substitutions enhance the binding activity of mutant p16 through several different mechanisms, including an increased amount of secondary structure and thermostability, reduction in exposed hydrophobic surface(s) and/or a reduced tendency to aggregate. This apparent global suppressor effect suggests that increasing the thermodynamic stability of p16 can be used as a general strategy to restore the biological activity to defective mutants of this important tumor suppressor protein.

Gln277 and Phe554 residues are involved in thermal inactivation of protective antigen of Bacillus anthracis

Protective antigen (PA) is the main component of all the vaccines against anthrax. The currently available vaccines have traces of other proteins that contribute to its reactogenicity. Thus, purified PA is recommended for human vaccination. PA loses its biological activity within 48h at 37 degrees C and its thermolability has been a cause of concern as accidental exposure to higher temperatures during transportation or storage could decrease its efficacy. In the present study, we have used protein engineering approach to increase the thermostability of PA by mutating amino acid residues on the surface as well as the interior of the protein. After screening several mutants, the mutants Gln277Ala and Phe554Ala have been found to be more thermostable than the wild-type PA. Gln277Ala retains approximately 45% and Phe554Ala retains approximately 90% activity, even after incubation at 37 degrees C for 48h while in the same period wild-type PA loses its biological activity completely. It is the first report of increasing thermostability of PA using site-directed mutagenesis. Generation of such mutants could pave the way for better anthrax vaccines with longer shelf life.

Knowledge-based potential functions in protein design

Predicting protein sequences that fold into specific native three-dimensional structures is a problem of great potential complexity. Although the complete solution is ultimately rooted in understanding the physical chemistry underlying the complex interactions between amino acid residues that determine protein stability, recent work shows that empirical information about these first principles is embedded in the statistics of protein sequence and structure databases. This review focuses on the use of 'knowledge-based' potentials derived from these databases in designing proteins. In addition, the data suggest how the study of these empirical potentials might impact our fundamental understanding of the energetic principles of protein structure.

Optimized production of the Diels-Alderase antibody 1E9 as a chimeric Fab

Monoclonal antibody 1E9, which catalyzes the [4+2] cycloaddition between tetrachlorothiophene dioxide and N-ethylmaleimide, has been re-engineered for production as a chimeric human–murine Fab fragment in Escherichia coli. Stabilizing point mutations in the variable regions of the antibody were identified by replacing residues that rarely occur at individual positions in aligned immunoglobulin sequences with their consensus counterparts. By combining favorable substitutions, double (MetH87Thr–GlyL63Ser) and triple (MetH87Thr–GlyL63Ser–PheL95Pro) mutants were created, which can be produced in good yield (4 and 17 mg L–1cell culture, respectively). The triple mutant exhibits a modest fourfold drop in the apparent kcatvalue for the cycloaddition reaction, but the kinetic properties of the double mutant are indistinguishable from those of the parent murine IgG. The availability of recombinant versions of this catalytic antibody will facilitate efforts to determine the origins of its selectivity and catalytic efficiency through mutagenesis.Key words: catalytic antibody, Fab fragment, bacterial production.

The consensus concept for thermostability engineering of proteins: further proof of concept

Previously, we calculated a consensus amino acid sequence from 13 homologous fungal phytases. A synthetic gene was constructed and recombinantly expressed. Surprisingly, consensus phytase-1 was 15-26 degrees C more thermostable than all parent phytases used in its design [Lehmann et al. (2000)Protein Eng., 13, 49-57]. In the present study, inclusion of six further phytase sequences in the amino acid sequence alignment resulted in the replacement of 38 amino acid residues in either one or both of the new consensus phytases-10 and -11. Since consensus phytase-10, again, was 7.4 degrees C more thermostable than consensus phytase-1, the thermostability effects of most of the 38 amino acid substitutions were tested by site-directed mutagenesis. Both stabilizing and destabilizing mutations were identified, but all affected the stability of the enzyme by <3 degrees C. The combination of all stabilizing amino acid exchanges in a multiple mutant of consensus phytase-1 increased the unfolding temperature from 78.0 to 88.5 degrees C. Likewise, back-mutation of four destabilizing amino acids and introduction of an additional stabilizing amino acid in consensus phytase-10 further increased the unfolding temperature from 85.4 to 90.4 degrees C. The thermostabilization achieved is the result of a combination of slight improvements from multiple amino acid exchanges rather than being the effect of a single or of just a few dominating mutations that have been introduced by chance. The present findings support the general validity of the consensus concept for thermostability engineering of proteins.

Engineering proteins for thermostability: the use of sequence alignments versus rational design and directed evolution

With the advent of directed evolution techniques, protein engineering has received a fresh impetus. Engineering proteins for thermostability is a particularly exciting and challenging field, as it is crucial for broadening the industrial use of recombinant proteins. In addition to directed evolution, a variety of partially successful rational concepts for engineering thermostability have been developed in the past. Recent results suggest that amino acid sequence comparisons of mesophilic proteins alone can be used efficiently to engineer thermostable proteins. The potential benefits of the underlying, semirational 'consensus concept' are compared with those of rational design and directed evolution approaches.

Structural genomics: opportunities and challenges

Following the complete genome sequencing of an increasing number of organisms, structural biology is engaging in a systematic approach of high-throughput structure determination called structural genomics to create a complete inventory of protein folds/structures that will help predict functions for all proteins. First results show that structural genomics will be highly effective in finding functional annotations for proteins of unknown function.

Increasing protein stability using a rational approach combining sequence homology and structural alignment: Stabilizing the WW domain

This study shows that a combination of sequence homology and structural information can be used to increase the stability of the WW domain by 2.5 kcal mol(-1) and increase the T(m) by 28 degrees C. Previous homology-based protein design efforts typically investigate positions with low sequence identity, whereas this study focuses on semi-conserved core residues and proximal residues, exploring their role(s) in mediating stabilizing interactions on the basis of structural considerations. The A20R and L30Y mutations allow increased hydrophobic interactions because of complimentary surfaces and an electrostatic interaction with a third residue adjacent to the ligand-binding hydrophobic cluster, increasing stability significantly beyond what additivity would predict for the single mutations. The D34T mutation situated in a pi-turn possibly disengages Asn31, allowing it to make up to three hydrogen bonds with the backbone in strand 1 and loop 2. The synergistic mutations A20R/L30Y in combination with the remotely located mutation D34T add together to create a hYap WW domain that is significantly more stable than any of the protein structures on which the design was based (Pin and FBP28 WW domains).

A thermophilic mini-chaperonin contains a conserved polypeptide-binding surface: combined crystallographic and NMR studies of the GroEL apical domain with implications for substrate interactions

A homologue of the Escherichia coli GroEL apical domain was obtained from thermophilic eubacterium Thermus thermophilus. The domains share 70 % sequence identity (101 out of 145 residues). The thermal stability of the T. thermophilus apical domain (Tm>100 degrees C as evaluated by circular dichroism) is at least 35 degrees C greater than that of the E. coli apical domain (Tm=65 degrees C). The crystal structure of a selenomethione-substituted apical domain from T. thermophilus was determined to a resolution of 1.78 A using multiwavelength-anomalous-diffraction phasing. The structure is similar to that of the E. coli apical domain (root-mean-square deviation 0.45 A based on main-chain atoms). The thermophilic structure contains seven additional salt bridges of which four contain charge-stabilized hydrogen bonds. Only one of the additional salt bridges would face the "Anfinsen cage" in GroEL. High temperatures were exploited to map sites of interactions between the apical domain and molten globules. NMR footprints of apical domain-protein complexes were obtained at elevated temperature using 15N-1H correlation spectra of 15N-labeled apical domain. Footprints employing two polypeptides unrelated in sequence or structure (an insulin monomer and the SRY high-mobility-group box, each partially unfolded at 50 degrees C) are essentially the same and consistent with the peptide-binding surface previously defined in E. coli GroEL and its apical domain-peptide complexes. An additional part of this surface comprising a short N-terminal alpha-helix is observed. The extended footprint rationalizes mutagenesis studies of intact GroEL in which point mutations affecting substrate binding were found outside the "classical" peptide-binding site. Our results demonstrate structural conservation of the apical domain among GroEL homologues and conservation of an extended non-polar surface recognizing diverse polypeptides.

From minichaperone to GroEL 1: information on GroEL-polypeptide interactions from crystal packing of minichaperones

We are reconstructing the mechanism of action of GroEL by a reductionist approach of isolating its minimal fragment that has residual activity (the "minichaperone" core) and then identifying how additional elements of structure confer further activity and function. We report here the 2.0 A resolution crystal structure of the minichaperone GroEL(193-345). The structure provides further clues on the nature of GroEL-polypeptide substrate interactions, because two molecules in the asymmetric unit interact by the binding of one molecule in the active site of its partner, thus mimicking a chaperone-polypeptide substrate complex. The results may explain some experimental observations, including the preference of GroEL for net positive charges (mediated by Glu238 and Glu257) and the key role of Tyr203 in mediating polypeptide binding. The larger binding site identified by these studies forms a continuous surface near the opening of the central cavity of GroEL that can accommodate a wide range of non-native protein conformations that differ in size and in structural and chemical properties.

From minichaperone to GroEL 3: properties of an active single-ring mutant of GroEL

The next step in our reductional analysis of GroEL was to study the activity of an isolated single seven-membered ring of the 14-mer. A known single-ring mutant, GroEL(SR1), contains four point mutations that prevent the formation of double-rings. That heptameric complex is functionally inactive because it is unable to release GroES. We found that the mutation E191G, which is responsible for the temperature sensitive (ts) Escherichia coli allele groEL44 and is located in the hinge region between the intermediate and apical domains of GroEL, appears to function by weakening the binding of GroES, without destabilizing the overall structure of GroEL44 mutant. We introduced, therefore, the mutation E191G into GroEL(SR1) in order to generate a single-ring mutant that may have weaker binding of GroES and hence be active. The new single-ring mutant, GroEL(SR44), was indeed effective in refolding both heat and dithiothreitol-denatured mitochondrial malate dehydrogenase with great efficiency. Further, unlike all smaller constructs of GroEL, the expression of GroEL(SR44) in E. coli that contained no endogenous GroEL restored biological viability, but not as efficiently as does wild-type GroEL. We envisage the notional evolution of the structure and properties of GroEL. The minichaperone core acts as a primitive chaperone by providing a binding surface for denatured states that prevents their self-aggregation. The assembly of seven minichaperones into a ring then enhances substrate binding by introducing avidity. The acquisition of binding sites for ATP then allows the modulation of substrate binding by introducing the allosteric mechanism that causes cycling between strong and weak binding sites. This is accompanied by the acquisition by the heptamer of the binding of GroES, which functions as a lid to the central cavity and competes for peptide binding sites. Finally, dimerization of the heptamer enhances its biological activity.

From minichaperone to GroEL 2: importance of avidity of the multisite ring structure

Structural studies on minichaperones and GroEL imply a continuous ring of binding sites around the neck of GroEL. To investigate the importance of this ring, we constructed an artificial heptameric assembly of minichaperones to mimic their arrangement in GroEL. The heptameric Gp31 co-chaperonin from bacteriophage T4, an analogue of GroES, was used as a scaffold to display the GroEL minichaperones. A fusion protein, MC(7), was generated by replacing a part of the highly mobile loop of Gp31 (residues 23-44) with the sequence of the minichaperone (residues 191-376 of GroEL). The purified recombinant protein assembled into a heptameric ring composed of seven 30.6 kDa subunits. Although single minichaperones (residues 193-335 to 191-376 of GroEL) have certain chaperone activities in vitro and in vivo, they cannot refold heat and dithiothreitol-denatured mitochondrial malate dehydrogenase (mtMDH), a reaction that normally requires GroEL, its co-chaperonin GroES and ATP. But, MC(7) refolded MDH in vitro. The expression of MC(7) complements in vivo two temperature-sensitive Escherichia coli alleles, groEL44 and groEL673, at 43 degrees C. Although MC(7) could not compensate for the complete absence of GroEL in vivo, it enhanced the colony-forming ability of cells containing limiting amounts of wild-type GroEL at 37 degrees C. MC(7 )also reduces aggregate formation and cell death in mammalian cell models of Huntington's disease. The assembly of seven minichaperone subunits on a heptameric ring significantly improves their activity, demonstrating the importance of avidity in GroEL function.