An integrated approach for thermal stabilization of a mesophilic adenylate kinase

Thermally stable proteins are desirable for research and industrial purposes, but redesigning proteins for higher thermal stability can be challenging. A number of different techniques have been used to improve the thermal stability of proteins, but the extents of stability enhancement were sometimes unpredictable and not significant. Here, we systematically tested the effects of multiple stabilization techniques including a bioinformatic method and structure‐guided mutagenesis on a single protein, thereby providing an integrated approach to protein thermal stabilization. Using a mesophilic adenylate kinase (AK) as a model, we identified stabilizing mutations based on various stabilization techniques, and generated a series of AK variants by introducing mutations both individually and collectively. The redesigned proteins displayed a range of increased thermal stabilities, the most stable of which was comparable to a naturally evolved thermophilic homologue with more than a 25° increase in its thermal denaturation midpoint. We also solved crystal structures of three representative variants including the most stable variant, to confirm the structural basis for their increased stabilities. These results provide a unique opportunity for systematically analyzing the effectiveness and additivity of various stabilization mechanisms, and they represent a useful approach for improving protein stability by integrating the reduction of local structural entropy and the optimization of global noncovalent interactions such as hydrophobic contact and ion pairs. Proteins 2014; 82:1947–1959. © 2014 Wiley Periodicals, Inc.     

Solution structure and backbone dynamics of the K18G/R82E Alicyclobacillus acidocaldarius thioredoxin mutant: a molecular analysis of its reduced thermal stability

No general strategy for thermostability has been yet established, because the extra stability of thermophiles appears to be the sum of different cumulative stabilizing interactions. In addition, the increase of conformational rigidity observed in many thermophilic proteins, which in some cases disappears when mesophilic and thermophilic proteins are compared at their respective physiological temperatures, suggests that evolutionary adaptation tends to maintain corresponding states with respect to conformational flexibility. In this study, we accomplished a structural analysis of the K18G/R82E Alicyclobacillus acidocaldarius thioredoxin (BacTrx) mutant, which has reduced heat resistance with respect to the thermostable wild-type. Furthermore, we have also achieved a detailed study, carried out at 25, 45, and 65 degrees C, of the backbone dynamics of both the BacTrx and its K18G/R82E mutant. Our findings clearly indicate that the insertion of the two mutations causes a loss of energetically favorable long-range interactions and renders the secondary structure elements of the double mutants more similar to those of the mesophilic Escherichia coli thioredoxin. Moreover, protein dynamics analysis shows that at room temperature the BacTrx, as well as the double mutant, are globally as rigid as the mesophilic thioredoxins; differently, at 65 degrees C, which is in the optimal growth temperature range of A. acidocaldarius, the wild-type retains its rigidity while the double mutant is characterized by a large increase of the amplitude of the internal motions. Finally, our research interestingly shows that fast motions on the pico- to nanosecond time scale are not detrimental to protein stability and provide an entropic stabilization of the native state. This study further confirms that protein thermostability is reached through diverse stabilizing interactions, which have the key role to maintain the structural folding stable and functional at the working temperature.     

PROTS: a fragment based protein thermo-stability potential

Designing proteins with enhanced thermo-stability has been a main focus of protein engineering because of its theoretical and practical significance. Despite extensive studies in the past years, a general strategy for stabilizing proteins still remains elusive. Thus effective and robust computational algorithms for designing thermo-stable proteins are in critical demand. Here we report PROTS, a sequential and structural four-residue fragment based protein thermo-stability potential. PROTS is derived from a nonredundant representative collection of thousands of thermophilic and mesophilic protein structures and a large set of point mutations with experimentally determined changes of melting temperatures. To the best of our knowledge, PROTS is the first protein stability predictor based on integrated analysis and mining of these two types of data. Besides conventional cross validation and blind testing, we introduce hypothetical reverse mutations as a means of testing the robustness of protein thermo-stability predictors. In all tests, PROTS demonstrates the ability to reliably predict mutation induced thermo-stability changes as well as classify thermophilic and mesophilic proteins. In addition, this white-box predictor allows easy interpretation of the factors that influence mutation induced protein stability changes at the residue level.     

Using motif-based methods in multiple genome analyses: a case study comparing orthologous mesophilic and thermophilic proteins

Protein motifs represent highly conserved regions within protein families and are generally accepted to describe critical regions required for protein stability and/or function. In this comprehensive analysis, we present a robust, unique approach to identify and compare corresponding mesophilic and thermophilic sequence motifs between all orthologous proteins within 44 microbial genomes. Motif similarity is determined through global sequence alignment of mesophilic and thermophilic motif pairs, which are identified by a greedy algorithm. Our results reveal only modest correlation between motif and overall sequence similarity, highlighting the rationale of motif-based approaches in comprehensive multigenome comparisons. Conserved mutations reflect previously suggested physiochemical principles for conferring thermostability. Additionally, comparisons between corresponding mesophilic and thermophilic motif pairs provide key biochemical insights related to thermostability and can be used to test the evolutionary robustness of individual structural comparisons. We demonstrate the ability of our unique approach to provide key insights in two examples: the TATA-box binding protein and glutamate dehydrogenase families. In the latter example, conserved mutations hint at novel origins leading to structural stability differences within the hexamer structures. Additionally, we present amino acid composition data and average protein length comparisons for all 44 microbial genomes.     

Crystal structure of histidine-containing phosphocarrier protein from <Emphasis Type="Italic">Thermoanaerobacter tengcongensis</Emphasis> MB4 and the implications for thermostability

Protein thermostability is an inherent characteristic of proteins from thermophilic microorganisms, and therefore enables these organisms to survive at extreme temperatures. Although it is well-known that thermostable proteins are critical for the growth of thermophilic organisms, the structural basis of protein thermostability is not yet fully understood. The histidine-containing phosphocarrier (HPr) protein, a phosphate shuttle protein in the phosphoenolpyruvate-dependent sugar transport system (PTS) of bacterial species, is an ideal model for investigating protein thermostability with respect to its small size and deficiency in disulphide bonds or cofactors. In this study, the HPr protein from Thermoanaerobacter tengcongensis (TtHPr) is cloned and purified. Crystal structure with good quality has been determined at 2.3 Å resolution, which provides a firm foundation for exploring the thermostable mechanism. However, it shows that the crystal structure is conserved and no clue can be obtained from this single structure. Furthermore, detailed comparison of sequence and structure with the homologs from meso- or thermophilic bacteria shows no obvious rule for thermostability, but the extra salt-bridge existing only in thermophilic bacteria might be a better explanation for thermostability of HPr. Thus, mutations are performed to interrupt the salt-bridge in HPrs in thermophilic bacteria. Using site-directed mutations and the circular dichroism method, thermostability is evaluated, and the mutational variations are shown to have a faster denaturing rate than for wild-type viruses, indicating that mutations cause instability in the HPrs. Understanding the higher-temperature resistance of thermophilic and hyperthermophilic proteins is essential to studies on protein folding and stability, and is critical in engineering efficient enzymes that can work at a high temperature.     

Phenotypic Interactions among Mutations in a Thermus thermophilus 16S rRNA Gene Detected with Genetic Selections and Experimental Evolution

During protein synthesis, the ribosome undergoes conformational transitions between functional states, requiring communication between distant structural elements of the ribosome. Despite advances in ribosome structural biology, identifying the protein and rRNA residues governing these transitions remains a significant challenge. Such residues can potentially be identified genetically, given the predicted deleterious effects of mutations stabilizing the ribosome in discrete conformations and the expected ameliorating effects of second-site compensatory mutations. In this study, we employed genetic selections and experimental evolution to identify interacting mutations in the ribosome of the thermophilic bacterium Thermus thermophilus. By direct genetic selections, we identified mutations in 16S rRNA conferring a streptomycin dependence phenotype and from these derived second-site suppressor mutations relieving dependence. Using experimental evolution of streptomycin-independent pseudorevertants, we identified additional compensating mutations. Similar mutations could be evolved from slow-growing streptomycin-resistant mutants. While some mutations arose close to the site of the original mutation in the three-dimensional structure of the 30S ribosomal subunit and probably act directly by compensating for local structural distortions, the locations of others are consistent with long-range communication between specific structural elements within the ribosome.     

Comparison of energy components of proteins from thermophilic and mesophilic organisms

In order to infer the energetic determinants of thermophilic proteins, molecular mechanics calculations were applied to five proteins from thermophilic eubacteria and their mesophilic homologs. The energy function includes a hydration term as well as the electrostatic contribution from the solvent in addition to the usual conformational energy terms. We calculated energy values for three different states of each protein: the native, near-native, and unfolded structures. The energy difference and its components between pairs of these states were compared. The hypothetical near-native structures have almost the same backbone conformation as the native structure but with largely distorted side-chain packing, thus enabling us to extract the energy components important for stabilizing the native backbone topology itself, irrespective of structural details. It was found that the sum of the electrostatic and hydration energies, although of large positive values, were consistently lower for the thermophilic proteins than for their mesophilic counterparts. This trend was observed in the energy difference not only between the native and unfolded structures, but also between the near-native and unfolded structures. In contrast, the energy components regarding side-chain packing did not show any clear tendency. These results suggest that the thermophilic proteins are stabilized so that the precise packing of the native structure does not significantly affect the stability. Implications of this conclusion are also discussed.     

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

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

Electrostatic contributions to the stability of hyperthermophilic proteins.

Electrostatic contributions to the folding free energy of several hyperthermophilic proteins and their mesophilic homologs are calculated. In all the cases studied, electrostatic interactions are more favorable in the hyperthermophilic proteins. The electrostatic free energy is found not to be correlated with the number of ionizable amino acid residues, ion pairs or ion pair networks in a protein, but rather depends on the location of these groups within the protein structure. Moreover, due to the large free energy cost associated with burying charged groups, buried ion pairs are found to be destabilizing unless they undergo favorable interactions with additional polar groups, including other ion pairs. The latter case involves the formation of stabilizing ion pair networks as is observed in a number of proteins. Ion pairs located on the protein surface also provide stabilizing interactions in a number of cases. Taken together, our results suggest that many hyperthermophilic proteins enhance electrostatic interactions through the optimum placement of charged amino acid residues within the protein structure, although different design strategies are used in different cases. Other physical mechanisms are also likely to contribute, however optimizing electrostatic interactions offers a simple means of enhancing stability without disrupting the core residues characteristic of different protein families.     

Structural differences between mesophilic, moderately thermophilic and extremely thermophilic protein subunits: results of a comprehensive survey

BACKGROUND: Proteins from thermophilic organisms usually show high intrinsic thermal stability but have structures that are very similar to their mesophilic homologues. From prevous studies it is difficult to draw general conclusions about the structural features underlying the increased thermal stability of thermophilic proteins. RESULTS: In order to reveal the general evolutionary strategy for changing the heat stability of proteins, a non-redundant data set was compiled comprising all high-quality structures of thermophilic proteins and their mesophilic homologues from the Protein Data Bank. The selection (quality) criteria were met by 64 mesophilic and 29 thermophilic protein subunits, representing 25 protein families. From the atomic coordinates, 13 structural parameters were calculated, compared and evaluated using statistical methods. This study is distinguished from earlier ones by the strict quality control of the structures used and the size of the data set. CONCLUSIONS: Different protein families adapt to higher temperatures by different sets of structural devices. Regarding the structural parameters, the only generally observed rule is an increase in the number of ion pairs with increasing growth temperature. Other parameters show just a trend, whereas the number of hydrogen bonds and the polarity of buried surfaces exhibit no clear-cut tendency to change with growth temperature. Proteins from extreme thermophiles are stabilized in different ways to moderately thermophilic ones. The preferences of these two groups are different with regards to the number of ion pairs, the number of cavities, the polarity of exposed surface and the secondary structural composition.     

Strategies for structural proteomics of prokaryotes: Quantifying the advantages of studying orthologous proteins and of using both NMR and X-ray crystallography approaches

Only about half of non-membrane-bound proteins encoded by either bacterial or archaeal genomes are soluble when expressed in Escherichia coli (Yee et al., Proc Natl Acad Sci USA 2002;99:1825-1830; Christendat et al., Prog Biophys Mol Biol 200;73:339-345). This property limits genome-scale functional and structural proteomics studies, which depend on having a recombinant, soluble version of each protein. An emerging strategy to increase the probability of deriving a soluble derivative of a protein is to study different sequence homologues of the same protein, including representatives from thermophilic organisms, based on the assumption that the stability of these proteins will facilitate structural analysis. To estimate the relative merits of this strategy, we compared the recombinant expression, solubility, and suitability for structural analysis by NMR and/or X-ray crystallography for 68 pairs of homologous proteins from E. coli and Thermotoga maritima. A sample suitable for structural studies was obtained for 62 of the 68 pairs of homologs under standardized growth and purification procedures. Fourteen (eight E. coli and six T. maritima proteins) samples generated NMR spectra of a quality suitable for structure determination and 30 (14 E. coli and 16 T. maritima proteins) samples formed crystals. Only three (one E. coli and two T. maritima proteins) samples both crystallized and had excellent NMR properties. The conclusions from this work are: (1) The inclusion of even a single ortholog of a target protein increases the number of samples for structural studies almost twofold; (2) there was no clear advantage to the use of thermophilic proteins to generate samples for structural studies; and (3) for the small proteins analyzed here, the use of both NMR and crystallography approaches almost doubled the number of samples for structural studies.     

Structural basis for altering the stability of homologous RNAs from a mesophilic and a thermophilic bacterium

Tertiary RNA structures from thermophilic bacteria generally are more stable than their mesophilic homologs. To understand the structural basis of the increase in stability, we investigated equilibrium folding of the specificity domain (S-domain) of RNase P RNA from a mesophilic (Escherichia coli) and a thermophilic (Thermus thermophilus) bacterium. Equilibrium folding of both S-domains is described by a minimal, three-state folding scheme, U-to-I-to-N. In the I-to-N transition of the thermophilic S-domain, more structure forms and protections are stronger against T1 nuclease and hydroxyl radical reactions. Phylogenetic comparison in the context of the native structure reveals that among 39 nucleotide differences between these S-domains, 12 likely contribute to higher stability. These residues participate in extensive networks of hydrogen bonding, stacking, and metal ion coordination throughout the molecule. The thermophilic S-domain achieves higher stability by mutating strategic base pairs to G-C, decreasing surface accessibility of the native state, and increasing the amount of structure formation in the native folding transition. An E. coli S-domain mutant containing these 12 nt has the same stability and folding cooperativity as the T. thermophilus S-domain. E. coli S-domain mutants containing a subset of 4 or 6 nt have the same stability as the T. thermophilus S-domain but the same folding cooperativity as the E. coli S-domain. These results show that increasing stability can be accomplished by mutations within a local structure, but increasing folding cooperativity needs concerted changes among multiple structural units.     

Hydrophobic environment is a key factor for the stability of thermophilic proteins

The stability of thermophilic proteins has been viewed from different perspectives and there is yet no unified principle to understand this stability. It would be valuable to reveal the most important interactions for designing thermostable proteins for such applications as industrial protein engineering. In this work, we have systematically analyzed the importance of various interactions by computing different parameters such as surrounding hydrophobicity, inter‐residue interactions, ion‐pairs and hydrogen bonds. The importance of each interaction has been determined by its predicted relative contribution in thermophiles versus the same contribution in mesophilic homologues based on a dataset of 373 protein families. We predict that hydrophobic environment is the major factor for the stability of thermophilic proteins and found that 80% of thermophilic proteins analyzed showed higher hydrophobicity than their mesophilic counterparts. Ion pairs, hydrogen bonds, and interaction energy are also important and favored in 68%, 50%, and 62% of thermophilic proteins, respectively. Interestingly, thermophilic proteins with decreased hydrophobic environments display a greater number of hydrogen bonds and/or ion pairs. The systematic elimination of mesophilic proteins based on surrounding hydrophobicity, interaction energy, and ion pairs/hydrogen bonds, led to correctly identifying 95% of the thermophilic proteins in our analyses. Our analysis was also applied to another, more refined set of 102 thermophilic–mesophilic pairs, which again identified hydrophobicity as a dominant property in 71% of the thermophilic proteins. Further, the notion of surrounding hydrophobicity, which characterizes the hydrophobic behavior of residues in a protein environment, has been applied to the three‐dimensional structures of elongation factor‐Tu proteins and we found that the thermophilic proteins are enriched with a hydrophobic environment. The results obtained in this work highlight the importance of hydrophobicity as the dominating characteristic in the stability of thermophilic proteins, and we anticipate this will be useful in our attempts to engineering thermostable proteins. © Proteins 2013. © 2012 Wiley Periodicals, Inc.     

Distance-dependent statistical potentials for discriminating thermophilic and mesophilic proteins

Identification of the characteristic structural patterns responsible for protein thermostability is theoretically important and practically useful but largely remains an open problem. These patterns may be revealed through comparative study on thermophilic and mesophilic proteins that have distinct thermostability. In this study, we constructed several distance-dependant potentials from thermophilic and mesophilic proteins. These potentials were then used to evaluate the structural difference between thermophilic and mesophilic proteins. We found that using the subtraction or division of the potentials derived from thermophilic and mesophilic proteins can dramatically increase the discriminatory ability. This approach revealed that the ability to distinct the subtle structural features responsible for protein thermostability may be effectively enhanced through rationally designed comparative study.     

Identification of stabilizing point mutations through mutagenesis of destabilized protein libraries

Although there have been recent transformative advances in the area of protein structure prediction, prediction of point mutations that improve protein stability remains challenging. It is possible to construct and screen large mutant libraries for improved activity or ligand binding. However, reliable screens for mutants that improve protein stability do not yet exist, especially for proteins that are well folded and relatively stable. Here, we demonstrate that incorporation of a single, specific, destabilizing mutation termed parent inactivating mutation into each member of a single-site saturation mutagenesis library, followed by screening for suppressors, allows for robust and accurate identification of stabilizing mutations. We carried out fluorescence-activated cell sorting of such a yeast surface display, saturation suppressor library of the bacterial toxin CcdB, followed by deep sequencing of sorted populations. We found that multiple stabilizing mutations could be identified after a single round of sorting. In addition, multiple libraries with different parent inactivating mutations could be pooled and simultaneously screened to further enhance the accuracy of identification of stabilizing mutations. Finally, we show that individual stabilizing mutations could be combined to result in a multi-mutant that demonstrated an increase in thermal melting temperature of about 20 °C, and that displayed enhanced tolerance to high temperature exposure. We conclude that as this method is robust and employs small library sizes, it can be readily extended to other display and screening formats to rapidly isolate stabilized protein mutants.     

Comparison of crystal structure interactions and thermodynamics for stabilizing mutations in the Tetrahymena ribozyme

Although general mechanisms of RNA folding and catalysis have been elucidated, little is known about how ribozymes achieve structural stability at high temperature. A previous in vitro evolution experiment identified a small number of mutations that significantly increase the thermostability of the tertiary structure of the Tetrahymena ribozyme. Because we also determined the crystal structure of this thermostable ribozyme, we have for the first time the opportunity to compare the structural interactions and thermodynamic contributions of individual nucleotides in a ribozyme. We investigated the contribution of five mutations to thermostability by using temperature gradient gel electrophoresis. Unlike the case with several well-studied proteins, the effects of individual mutations on thermostability of this RNA were highly context dependent. The three most important mutations for thermostability were actually destabilizing in the wild-type background. A269G and A304G contributed to stability only when present as a pair, consistent with their proximity in the ribozyme structure. In an evolutionary context, this work supports and extends the idea that one advantage of protein enzyme systems over an RNA world is the ability of proteins to accumulate stabilizing single-site mutations, whereas RNA may often require much rarer double mutations to improve the stability of both its tertiary and secondary structures.     

Structural adaptation to low temperatures--analysis of the subunit interface of oligomeric psychrophilic enzymes

Enzymes from psychrophiles show higher catalytic efficiency in the 0-20 degrees C temperature range and often lower thermostability in comparison with meso/thermophilic homologs. Physical and chemical characterization of these enzymes is currently underway in order to understand the molecular basis of cold adaptation. Psychrophilic enzymes are often characterized by higher flexibility, which allows for better interaction with substrates, and by a lower activation energy requirement in comparison with meso/thermophilic counterparts. In their tertiary structure, psychrophilic enzymes present fewer stabilizing interactions, longer and more hydrophilic loops, higher glycine content, and lower proline and arginine content. In this study, a comparative analysis of the structural characteristics of the interfaces between oligomeric psychrophilic enzyme subunits was carried out. Crystallographic structures of oligomeric psychrophilic enzymes, and their meso/thermophilic homologs belonging to five different protein families, were retrieved from the Protein Data Bank. The following structural parameters were calculated: overall and core interface area, characterization of polar/apolar contributions to the interface, hydrophobic contact area, quantity of ion pairs and hydrogen bonds between monomers, internal area and total volume of non-solvent-exposed cavities at the interface, and average packing of interface residues. These properties were compared to those of meso/thermophilic enzymes. The results were analyzed using Student's t-test. The most significant differences between psychrophilic and mesophilic proteins were found in the number of ion pairs and hydrogen bonds, and in the apolarity of their subunit interface. Interestingly, the number of ion pairs at the interface shows an opposite adaptation to those occurring at the monomer core and surface.     

Conformational changes in redox pairs of protein structures

Disulfides are conventionally viewed as structurally stabilizing elements in proteins but emerging evidence suggests two disulfide subproteomes exist. One group mediates the well known role of structural stabilization. A second redox‐active group are best known for their catalytic functions but are increasingly being recognized for their roles in regulation of protein function. Redox‐active disulfides are, by their very nature, more susceptible to reduction than structural disulfides; and conversely, the Cys pairs that form them are more susceptible to oxidation. In this study, we searched for potentially redox‐active Cys Pairs by scanning the Protein Data Bank for structures of proteins in alternate redox states. The PDB contains over 1134 unique redox pairs of proteins, many of which exhibit conformational differences between alternate redox states. Several classes of structural changes were observed, proteins that exhibit: disulfide oxidation following expulsion of metals such as zinc; major reorganisation of the polypeptide backbone in association with disulfide redox‐activity; order/disorder transitions; and changes in quaternary structure. Based on evidence gathered supporting disulfide redox activity, we propose disulfides present in alternate redox states are likely to have physiologically relevant redox activity.