A thermodynamic comparison of HPr proteins from extremophilic organisms

Effect of the Lys62Ala Mutation on the Thermal Stability of BstHPr Protein by Molecular Dynamics

We analyzed the thermal stability of the BstHPr protein through the site-directed point mutation Lys62 replaced by Ala residue using molecular dynamics simulations at five different temperatures: 298, 333, 362, 400, and 450 K, for periods of 1 μs and in triplicate. The results from the mutant thermophilic BstHPrm protein were compared with those of the wild-type thermophilic BstHPr protein and the mesophilic BsHPr protein. Structural and molecular interaction analyses show that proteins lose stability as temperature increases. Mutant and wild-type proteins behave similarly up to 362 K. However, at 400 K the mutant protein shows greater structural instability, losing more buried hydrogen bonds and exposing more of its non-polar residues to the solvent. Therefore, in this study, we confirmed that the salt bridge network of the Glu3-Lys62-Glu36 triad, made up of the Glu3-Lys62 and Glu36-Lys62 ion pairs, provides thermal stability to the thermophilic BstHPr protein.

Molecular Dynamics Simulations of HPr Proteins from a Thermophilic and a Mesophilic Organism: A Comparative Thermal Study

The histidine-containing phosphocarrier (HPr) is a monomeric protein conserved in Gram-positive bacteria, which may be of mesophilic or thermophilic nature. In particular, the HPr protein from the thermophilic organism B. stearothermophilus is a good model system for thermostability studies, since experimental data, such as crystal structure and thermal stability curves, are available. However, its unfolding mechanism at higher temperatures is yet unclear at a molecular level. Therefore, in this work, we researched the thermal stability of this protein using molecular dynamics simulations, subjecting it to five different temperatures during a time span of 1 μs. The analyses of the structural parameters and molecular interactions were compared with those of the mesophilic homologue HPr protein from B. subtilis. Each simulation was run in triplicate using identical conditions for both proteins. The results showed that the two proteins lose stability as the temperature increases, but the mesophilic structure is more affected. We found that the salt bridge network formed by the triad of Glu3-Lys62-Glu36 residues and the salt bridge made up of Asp79-Lys83 ion pair are key factors to keep stable the thermophilic protein, maintaining the hydrophobic core protected and the structure packed. In addition, these molecular interactions neutralize the negative surface charge, acting as "natural molecular staples".

The Molecular Basis for Life in Extreme Environments

Sampling and genomic efforts over the past decade have revealed an enormous quantity and diversity of life in Earth's extreme environments. This new knowledge of life on Earth poses the challenge of understandingits molecular basis in such inhospitable conditions, given that such conditions lead to loss of structure and of function in biomolecules from mesophiles. In this review, we discuss the physicochemical properties of extreme environments. We present the state of recent progress in extreme environmental genomics. We then present an overview of our current understanding of the biomolecular adaptation to extreme conditions. As our current and future understanding of biomolecular structure-function relationships in extremophiles requires methodologies adapted to extremes of pressure, temperature, and chemical composition, advances in instrumentation for probing biophysical properties under extreme conditions are presented. Finally, we briefly discuss possible future directions in extreme biophysics.

Structural evidence for solvent-stabilisation by aspartic acid as a mechanism for halophilic protein stability in high salt concentrations

Halophilic organisms have adapted to survive in high salt environments, where mesophilic organisms would perish. One of the biggest challenges faced by halophilic proteins is the ability to maintain both the structure and function at molar concentrations of salt. A distinct adaptation of halophilic proteins, compared to mesophilic homologues, is the abundance of aspartic acid on the protein surface. Mutagenesis and crystallographic studies of halophilic proteins suggest an important role for solvent interactions with the surface aspartic acid residues. This interaction, between the regions of the acidic protein surface and the solvent, is thought to maintain a hydration layer around the protein at molar salt concentrations thereby allowing halophilic proteins to retain their functional state. Here we present neutron diffraction data of the monomeric zwitterionic form of aspartic acid solutions at physiological pH in 0.25 M and 2.5 M concentration of potassium chloride, to mimic mesophilic and halophilic-like environmental conditions. We have used isotopic substitution in combination with empirical potential structure refinement to extract atomic-scale information from the data. Our study provides structural insights that support the hypothesis that carboxyl groups on acidic residues bind water more tightly under high salt conditions, in support of the residue-ion interaction model of halophilic protein stabilisation. Furthermore our data show that in the presence of high salt the self-association between the zwitterionic form of aspartic acid molecules is reduced, suggesting a possible mechanism through which protein aggregation is prevented.

An Aromatic Cap Seals the Substrate Binding Site in an ECF-Type S Subunit for Riboflavin

ECF transporters are a family of active membrane transporters for essential micronutrients, such as vitamins and trace metals. Found exclusively in archaea and bacteria, these transporters are composed of four subunits: an integral membrane substrate-binding subunit (EcfS), a transmembrane coupling subunit (EcfT), and two ATP-binding cassette ATPases (EcfA and EcfA'). We have characterized the structural basis of substrate binding by the EcfS subunit for riboflavin from Thermotoga maritima, TmRibU. TmRibU binds riboflavin with high affinity, and the protein-substrate complex is exceptionally stable in solution. The crystal structure of riboflavin-bound TmRibU reveals an electronegative binding pocket at the extracellular surface in which the substrate is completely buried. Analysis of the intermolecular contacts indicates that nearly every available substrate hydrogen bond is satisfied. A conserved aromatic residue at the extracellular end of TM5, Tyr130, caps the binding site to generate a substrate-bound, occluded state, and non-conservative mutation of Tyr130 reduces the stability of this conformation. Using a novel fluorescence binding assay, we find that an aromatic residue at this position is essential for high-affinity substrate binding. Comparison with other S subunit structures suggests that TM5 and Loop5-6 contain a dynamic, conserved motif that plays a key role in gating substrate entry and release by S subunits of ECF transporters.

The histidine-phosphocarrier protein of the phosphoenolpyruvate: sugar phosphotransferase system of Bacillus sphaericus self-associates

The phosphotransferase system (PTS) is involved in the use of carbon sources in bacteria. Bacillus sphaericus, a bacterium with the ability to produce insecticidal proteins, is unable to use hexoses and pentoses as the sole carbon source, but it has ptsHI genes encoding the two general proteins of the PTS: enzyme I (EI) and the histidine phosphocarrier (HPr). In this work, we describe the biophysical and structural properties of HPr from B. sphaericus, HPr^bs, and its affinity towards EI of other species to find out whether there is inter-species binding. Conversely to what happens to other members of the HPr family, HPr^bs forms several self-associated species. The conformational stability of the protein is low, and it unfolds irreversibly during heating. The protein binds to the N-terminal domain of EI from Streptomyces coelicolor, EIN^sc, with a higher affinity than that of the natural partner of EIN^sc, HPr^sc. Modelling of the complex between EIN^sc and HPr^bs suggests that binding occurs similarly to that observed in other HPr species. We discuss the functional implications of the oligomeric states of HPr^bs for the glycolytic activity of B. sphaericus, as well as a strategy to inhibit binding between HPr^sc and EIN^sc.

Is it possible to stabilize a thermophilic protein further using sequences and structures of mesophilic proteins: a theoretical case study concerning DgAS

Incorporating structural elements of thermostable homologs can greatly improve the thermostability of a mesophilic protein. Despite the effectiveness of this method, applying it is often hampered. First, it requires alignment of the target mesophilic protein sequence with those of thermophilic homologs, but not every mesophilic protein has a thermophilic homolog. Second, not all favorable features of a thermophilic protein can be incorporated into the structure of a mesophilic protein. Furthermore, even the most stable native protein is not sufficiently stable for industrial applications. Therefore, creating an industrially applicable protein on the basis of the thermophilic protein could prove advantageous. Amylosucrase (AS) can catalyze the synthesis of an amylose-like polysaccharide composed of only α-1,4-linkages using sucrose as the lone energy source. However, industrial development of AS has been hampered owing to its low thermostability. To facilitate potential industrial applications, the aim of the current study was to improve the thermostability of Deinococcus geothermalis amylosucrase (DgAS) further; this is the most stable AS discovered to date. By integrating ideas from mesophilic AS with well-established protein design protocols, three useful design protocols are proposed, and several promising substitutions were identified using these protocols. The successful application of this hybrid design method indicates that it is possible to stabilize a thermostable protein further by incorporating structural elements of less-stable homologs.

The structure of Aquifex aeolicus ribosomal protein S8 reveals a unique subdomain that contributes to an extremely tight association with 16S rRNA

The assembly of ribonucleoprotein complexes occurs under a broad range of conditions, but the principles that promote assembly and allow function at high temperature are poorly understood. The ribosomal protein S8 from Aquifex aeolicus (AS8) is unique in that there is a 41-residue insertion in the consensus S8 sequence. In addition, AS8 exhibits an unusually high affinity for the 16S ribosomal RNA, characterized by a picomolar dissociation constant that is approximately 26,000-fold tighter than the equivalent interaction from Escherichia coli. Deletion analysis demonstrated that binding to the minimal site on helix 21 occurred at the same nanomolar affinity found for other bacterial species. The additional affinity required the presence of a three-helix junction between helices 20, 21, and 22. The crystal structure of AS8 was solved, revealing the helix-loop-helix geometry of the unique AS8 insertion region, while the core of the molecule is conserved with known S8 structures. The AS8 structure was modeled onto the structure of the 30S ribosomal subunit from E. coli, suggesting the possibility that the unique subdomain provides additional backbone and side-chain contacts between the protein and an unpaired base within the three-way junction of helices 20, 21, and 22. Point mutations in the protein insertion subdomain resulted in a significantly reduced RNA binding affinity with respect to wild-type AS8. These results indicate that the AS8-specific subdomain provides additional interactions with the three-way junction that contribute to the extremely tight binding to ribosomal RNA.

How do thermophilic proteins and proteomes withstand high temperature?

We attempt to understand the origin of enhanced stability in thermophilic proteins by analyzing thermodynamic data for 116 proteins, the largest data set achieved to date. We compute changes in entropy and enthalpy at the convergence temperature where different driving forces are maximally decoupled, in contrast to the majority of previous studies that were performed at the melting temperature. We find, on average, that the gain in enthalpy upon folding is lower in thermophiles than in mesophiles, whereas the loss in entropy upon folding is higher in mesophiles than in thermophiles. This implies that entropic stabilization may be responsible for the high melting temperature, and hints at residual structure or compactness of the denatured state in thermophiles. We find a similar trend by analyzing a homologous set of proteins classified based only on the optimum growth temperature of the organisms from which they were extracted. We find that the folding free energy at the temperature of maximal stability is significantly more favorable in thermophiles than in mesophiles, whereas the maximal stability temperature itself is similar between these two classes. Furthermore, we extend the thermodynamic analysis to model the entire proteome. The results explain the high optimal growth temperature in thermophilic organisms and are in excellent quantitative agreement with full thermal growth rate data obtained in a dozen thermophilic and mesophilic organisms.

Crystal structure of histidine-containing phosphocarrier protein from Thermoanaerobacter tengcongensis 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.

Evolution and thermodynamics of the slow unfolding of hyperstable monomeric proteins

BACKGROUND

The unfolding speed of some hyperthermophilic proteins is dramatically lower than that of their mesostable homologs. Ribonuclease HII from the hyperthermophilic archaeon Thermococcus kodakaraensis (Tk-RNase HII) is stabilized by its remarkably slow unfolding rate, whereas RNase HI from the thermophilic bacterium Thermus thermophilus (Tt-RNase HI) unfolds rapidly, comparable with to that of RNase HI from Escherichia coli (Ec-RNase HI).

RESULTS

To clarify whether the difference in the unfolding rate is due to differences in the types of RNase H or differences in proteins from archaea and bacteria, we examined the equilibrium stability and unfolding reaction of RNases HII from the hyperthermophilic bacteria Thermotoga maritima (Tm-RNase HII) and Aquifex aeolicus (Aa-RNase HII) and RNase HI from the hyperthermophilic archaeon Sulfolobus tokodaii (Sto-RNase HI). These proteins from hyperthermophiles are more stable than Ec-RNase HI over all the temperature ranges examined. The observed unfolding speeds of all hyperstable proteins at the different denaturant concentrations studied are much lower than those of Ec-RNase HI, which is in accordance with the familiar slow unfolding of hyperstable proteins. However, the unfolding rate constants of these RNases H in water are dispersed, and the unfolding rate constant of thermophilic archaeal proteins is lower than that of thermophilic bacterial proteins.

CONCLUSIONS

These results suggest that the nature of slow unfolding of thermophilic proteins is determined by the evolutionary history of the organisms involved. The unfolding rate constants in water are related to the amount of buried hydrophobic residues in the tertiary structure.

Stability curves of laboratory evolved thermostable mutants of a Bacillus subtilis lipase

Shape of the protein stability curves changes to achieve higher melting temperature. Broadly, these changes have been classified as upward shift (increased G(s)), rightward shift (increase in T(s)) and flattening of the stability curves (decrease in C(p)). Comparative studies on homologous mesophilic-thermophilic protein pairs highlighted the differential contribution of these three strategies amongst proteins. But unambiguous way of identification of the strategies, which will be preferred for a protein, is still not achieved. We have performed comparative thermodynamic studies using differential scanning calorimeter (DSC) on thermostable variants of a lipase from Bacillus subtilis. These variants are products of 1, 2, 3 and 4 rounds of directed evolution and harbor mutations having definite contribution in thermostability unlike natural thermophilic proteins. We have shown that upward and rightward shift in stability curves are prime strategies in this lipase. Our results along with that from the other study on laboratory evolved xylanase A suggest that optimization of suboptimal thermodynamic parameters is having a dominant influence in selection of thermodynamic strategies for higher thermostability.

Increasing protein stability by improving beta-turns

Our goal was to gain a better understanding of how protein stability can be increased by improving beta-turns. We studied 22 beta-turns in nine proteins with 66-370 residues by replacing other residues with proline and glycine and measuring the stability. These two residues are statistically preferred in some beta-turn positions. We studied: Cold shock protein B (CspB), Histidine-containing phosphocarrier protein, Ubiquitin, Ribonucleases Sa2, Sa3, T1, and HI, Tryptophan synthetase alpha-subunit, and Maltose binding protein. Of the 15 single proline mutations, 11 increased stability (Average = 0.8 +/- 0.3; Range = 0.3-1.5 kcal/mol), and the stabilizing effect of double proline mutants was additive. On the basis of this and our previous work, we conclude that proteins can generally be stabilized by replacing nonproline residues with proline residues at the i + 1 position of Type I and II beta-turns and at the i position in Type II beta-turns. Other turn positions can sometimes be used if the phi angle is near -60 degrees for the residue replaced. It is important that the side chain of the residue replaced is less than 50% buried. Identical substitutions in beta-turns in related proteins give similar results. Proline substitutions increase stability mainly by decreasing the entropy of the denatured state. In contrast, the large, diverse group of proteins considered here had almost no residues in beta-turns that could be replaced by Gly to increase protein stability. Improving beta-turns by substituting Pro residues is a generally useful way of increasing protein stability.

Slow unfolding of monomeric proteins from hyperthermophiles with reversible unfolding

Based on the differences in their optimal growth temperatures microorganisms can be classified into psychrophiles, mesophiles, thermophiles, and hyperthermophiles. Proteins from hyperthermophiles generally exhibit greater stability than those from other organisms. In this review, we collect data about the stability and folding of monomeric proteins from hyperthermophilies with reversible unfolding, from the equilibrium and kinetic aspects. The results indicate that slow unfolding is a general strategy by which proteins from hyperthermophiles adapt to higher temperatures. Hydrophobic interaction is one of the factors in the molecular mechanism of the slow unfolding of proteins from hyperthermophiles.

Proline effect on the thermostability and slow unfolding of a hyperthermophilic protein

Ribonuclease HII from hyperthermophile Thermococcus kodakaraensis (Tk-RNase HII) is a robust monomeric protein under kinetic control, which possesses some proline residues at the N-terminal of alpha-helices. Proline residue at the N-terminal of an alpha-helix is thought to stabilize a protein. In this work, the thermostability and folding kinetics of Tk-RNase HII were measured for mutant proteins in which a proline residue is introduced (Xaa to Pro) or removed (Pro to Ala) at the N-terminal of alpha-helices. In the folding experiments, the mutant proteins examined exhibit little influence on the remarkably slow unfolding of Tk-RNase HII. In contrast, E111P and K199P exhibit some thermostabilization, whereas P46A, P70A and P174A have some thermodestabilization. E111P/K199P and P46A/P70A double mutations cause cumulative changes in stability. We conclude that the proline effect on protein thermostability is observed in a hyperthermophilic protein, but each proline residue at the N-terminal of an alpha-helix slightly contributes to the thermostability. The present results also mean that even a natural hyperthermophilic protein can acquire improved thermostability.

Extreme temperature tolerance of a hyperthermophilic protein coupled to residual structure in the unfolded state

Understanding the mechanisms that dictate protein stability is of large relevance, for instance, to enable design of temperature-tolerant enzymes with high enzymatic activity over a broad temperature interval. In an effort to identify such mechanisms, we have performed a detailed comparative study of the folding thermodynamics and kinetics of the ribosomal protein S16 isolated from a mesophilic (S16(meso)) and hyperthermophilic (S16(thermo)) bacterium by using a variety of biophysical methods. As basis for the study, the 2.0 A X-ray structure of S16(thermo) was solved using single wavelength anomalous dispersion phasing. Thermal unfolding experiments yielded midpoints of 59 and 111 degrees C with associated changes in heat capacity upon unfolding (DeltaC(p)(0)) of 6.4 and 3.3 kJ mol(-1) K(-1), respectively. A strong linear correlation between DeltaC(p)(0) and melting temperature (T(m)) was observed for the wild-type proteins and mutated variants, suggesting that these variables are intimately connected. Stopped-flow fluorescence spectroscopy shows that S16(meso) folds through an apparent two-state model, whereas S16(thermo) folds through a more complex mechanism with a marked curvature in the refolding limb indicating the presence of a folding intermediate. Time-resolved energy transfer between Trp and N-(4,4-difluoro-5,7-dimethyl-4-bora-3a,4a-diaza-s-indacene-3-yl)methyl iodoacetamide of proteins mutated at selected positions shows that the denatured state ensemble of S16(thermo) is more compact relative to S16(meso). Taken together, our results suggest the presence of residual structure in the denatured state ensemble of S16(thermo) that appears to account for the large difference in quantified DeltaC(p)(0) values and, in turn, parts of the observed extreme thermal stability of S16(thermo). These observations may be of general importance in the design of robust enzymes that are highly active over a wide temperature span.

Prediction of protein stability upon point mutations

Prediction of protein stability upon amino acid substitution is a challenging problem and it will be helpful for designing stable mutants. We have developed a thermodynamic database for proteins and mutants (ProTherm), which has more than 20000 thermodynamic data along with sequence and structure information, experimental conditions and literature information. It is freely accessible at http://gibk26.bse.kyutech.ac.jp/jouhou/protherm/protherm.html. Utilizing the database, we have analysed the relationship between amino acid properties and protein stability and developed different methods, such as average assignment method, distance and torsion potentials and decision tree models to discriminate the stabilizing and destabilizing mutants, and to predict the stability change upon mutation. Our method could distinguish the stabilizing and destabilizing mutants with an accuracy of 82 and 85% respectively from amino acid sequence and protein three-dimensional structure. We obtained the correlation of 0.70 and 0.87, between the experimental and predicted stability changes upon mutations, from sequence and structure respectively. Furthermore, we have developed different web servers for discrimination and prediction and they are freely accessible at http://bioinformatics.myweb.hinet.net/iptree.htm and http://cupsat.tu-bs.de/.

Lessons in stability from thermophilic proteins

Studies that compare proteins from thermophilic and mesophilic organisms can provide insights into ability of thermophiles to function at their high habitat temperatures and may provide clues that enable us to better define the forces that stabilize all proteins. Most of the comparative studies have focused on thermal stability and show, as expected, that thermophilic proteins have higher Tm values than their mesophilic counterparts. Although these comparisons are useful, more detailed thermodynamic analyses are required to reach a more complete understanding of the mechanisms thermophilic protein employ to remain folded over a wider range of temperatures. This complete thermodynamic description allows one to generate a stability curve for a protein that defines how the conformational stability (DeltaG) varies with temperature. Here we compare stability curves for many pairs of homologous proteins from thermophilic and mesophilc organisms. Of the basic methods that can be employed to achieve enhanced thermostability, we find that most thermophilic proteins use the simple method that raises the DeltaG at all temperatures as the principal way to increase their Tm. We discuss and compare this thermodynamic method with the possible alternatives. In addition we propose ways that structural alterations and changes to the amino acid sequences might give rise to varied methods used to obtain thermostability.