Thermodynamic stability of archaeal histones

Histone variants in archaea - An undiscovered country

Exchanging core histones in the nucleosome for paralogous variants can have important functional ramifications. Many of these variants, and their physiological roles, have been characterized in exquisite detail in model eukaryotes, including humans. In comparison, our knowledge of histone biology in archaea remains rudimentary. This is true in particular for our knowledge of histone variants. Many archaea encode several histone genes that differ in sequence, but do these paralogs make distinct, adaptive contributions to genome organization and regulation in a manner comparable to eukaryotes? Below, we review what we know about histone variants in archaea at the level of structure, regulation, and evolution. In all areas, our knowledge pales when compared to the wealth of insight that has been gathered for eukaryotes. Recent findings, however, provide tantalizing glimpses into a rich and largely undiscovered country that is at times familiar and eukaryote-like and at times strange and uniquely archaeal. We sketch a preliminary roadmap for further exploration of this country; an undertaking that may ultimately shed light not only on chromatin biology in archaea but also on the origin of histone-based chromatin in eukaryotes.

A conserved π-helix plays a key role in thermoadaptation of catalysis in the glycoside hydrolase family 4

Here, we characterize the role of a π-helix in the molecular mechanisms underlying thermoadaptation in the glycoside hydrolase family 4 (GH4). The interspersed π-helix present in a subgroup is evolutionarily related to a conserved α-helix in other orthologs by a single residue insertion/deletion event. The insertional residue, Phe407, in a hyperthermophilic α-glucuronidase, makes specific interactions across the inter-subunit interface. In order to establish the sequence-structure-stability implications of the π-helix, the wild-type and the deletion variant (Δ407) were characterized. The variant showed a significant lowering of melting temperature and optimum temperature for the highest activity. Crystal structures of the proteins show a transformation of the π-helix to a continuous α-helix in the variant, identical to that in orthologs lacking this insertion. Thermodynamic parameters were determined from stability curves representing the temperature dependence of unfolding free energy. Though the proteins display maximum stabilities at similar temperatures, a higher melting temperature in the wild-type is achieved by a combination of higher enthalpy and lower heat capacity of unfolding. Comparisons of the structural changes, and the activity and thermodynamic profiles allow us to infer that specific non-covalent interactions, and the existence of residual structure in the unfolded state, are crucial determinants of its thermostability. These features permit the enzyme to balance the preservation of structure at a higher temperature with the thermodynamic stability required for optimum catalysis.

Protein thermostability engineering

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

Energetics of oligomeric protein folding and association

In nature, proteins most often exist as complexes, with many of these consisting of identical subunits. Understanding of the energetics governing the folding and misfolding of such homooligomeric proteins is central to understanding their function and misfunction, in disease or biotechnology. Much progress has been made in defining the mechanisms and thermodynamics of homooligomeric protein folding. In this review, we outline models as well as calorimetric and spectroscopic methods for characterizing oligomer folding, and describe extensive results obtained for diverse proteins, ranging from dimers to octamers and higher order aggregates. To our knowledge, this area has not been reviewed comprehensively in years, and the collective progress is impressive. The results provide evolutionary insights into the development of subunit interfaces, mechanisms of oligomer folding, and contributions of oligomerization to protein stability, function and regulation. Thermodynamic analyses have also proven valuable for understanding protein misfolding and aggregation mechanisms, suggesting new therapeutic avenues. Successful recent designs of novel, functional proteins demonstrate increased understanding of oligomer folding. Further rigorous analyses using multiple experimental and computational approaches are still required, however, to achieve consistent and accurate prediction of oligomer folding energetics. Modeling the energetics remains challenging but is a promising avenue for future advances.

The impact of solubility and electrostatics on fibril formation by the H3 and H4 histones

The goal of this study was to examine fibril formation by the heterodimeric eukaryotic histones (H2A-H2B and H3-H4) and homodimeric archaeal histones (hMfB and hPyA1). The histone fold dimerization motif is an obligatorily domain-swapped structure comprised of two fused helix:β-loop:helix motifs. Domain swapping has been proposed as a mechanism for the evolution of protein oligomers as well as a means to form precursors in the formation of amyloid-like fibrils. Despite sharing a common fold, the eukaryotic histones of the core nucleosome and archaeal histones fold by kinetic mechanisms of differing complexity with transient population of partially folded monomeric and/or dimeric species. No relationship was apparent between fibrillation propensity and equilibrium stability or population of kinetic intermediates. Only H3 and H4, as isolated monomers and as a heterodimer, readily formed fibrils at room temperature, and this propensity correlates with the significantly lower solubility of these polypeptides. The fibrils were characterized by ThT fluorescence, FTIR, and far-UV CD spectroscopies and electron microscopy. The helical histone fold comprises the protease-resistant core of the fibrils, with little or no protease protection of the poorly structured N-terminal tails. The highly charged tails inhibit fibrillation through electrostatic repulsion. Kinetic studies indicate that H3 and H4 form a co-fibril, with simultaneous incorporation of both histones. The potential impact of H3 and H4 fibrillation on the cytotoxicity of extracellular histones and α-synuclein-mediated neurotoxicity and fibrillation is considered.

Experimental evidence for the role of domain swapping in the evolution of the histone fold

The histone fold forms the fundamental endoskeleton of the protein core of the nucleosome and is also found in several transcription factors. We have investigated the evolutionary origins of this ubiquitous protein motif, which is found soluble exclusively as an antiparallel (handshake motif) dimer. We introduced a three amino acid insertion into the middle of a homodimeric archaeal histone fold motif. The engineered molecule was found to be a soluble and stable monomer with properties consistent with a four-helix-bundle protein. The experimental evidence presented here support the hypothesis that the handshake association motif characteristic of present-day histone dimers is the evolutionary product of domain swapping between two four-helix bundle domains, each of which derived from the tandem duplication of a primitive helix-strand-helix unit.

Thermal adaptation of viruses and bacteria

A previously established multiscale population genetics model posits that fitness can be inferred from the physical properties of proteins under the physiological assumption that a loss of stability by any protein confers the lethal phenotype to an organism. Here, we develop this model further by positing that replication rate (fitness) of a bacterial or viral strain directly depends on the copy number of folded proteins, which determine its replication rate. Using this model, and both numerical and analytical approaches, we studied the adaptation process of bacteria and viruses at varied environmental temperatures. We found that a broad distribution of protein stabilities observed in the model and in experiment is the key determinant of thermal response for viruses and bacteria. Our results explain most of the earlier experimental observations: the striking asymmetry of thermal response curves; the absence of evolutionary tradeoff, which was expected but not found in experiments; correlation between denaturation temperature for several protein families and the optimal growth temperature of their carrier organisms; and proximity of bacterial or viral optimal growth temperatures to their evolutionary temperatures. Our theory quantitatively and with high accuracy described thermal response curves for 35 bacterial species using, for each species, only two adjustable parameters-the number of rate-determining genes and the energy barrier for metabolic reactions.

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.

Photorhabdus adhesion modification protein (Pam) binds extracellular polysaccharide and alters bacterial attachment

BACKGROUND

Photorhabdus are Gram-negative nematode-symbiotic and insect-pathogenic bacteria. The species Photorhabdus asymbiotica is able to infect humans as well as insects. We investigated the secreted proteome of a clinical isolate of P. asymbiotica at different temperatures in order to identify proteins relevant to the infection of the two different hosts.

RESULTS

A comparison of the proteins secreted by a clinical isolate of P. asymbiotica at simulated insect (28 degrees C) and human (37 degrees C) temperatures led to the identification of a small and highly abundant protein, designated Pam, that is only secreted at the lower temperature. The pam gene is present in all Photorhabdus strains tested and shows a high level of conservation across the whole genus, suggesting it is both ancestral to the genus and probably important to the biology of the bacterium. The Pam protein shows limited sequence similarity to the 13.6 kDa component of a binary toxin of Bacillus thuringiensis. Nevertheless, injection or feeding of heterologously produced Pam showed no insecticidal activity to either Galleria mellonella or Manduca sexta larvae. In bacterial colonies, Pam is associated with an extracellular polysaccharide (EPS)-like matrix, and modifies the ability of wild-type cells to attach to an artificial surface. Interestingly, Surface Plasmon Resonance (SPR) binding studies revealed that the Pam protein itself has adhesive properties. Although Pam is produced throughout insect infection, genetic knockout does not affect either insect virulence or the ability of P. luminescens to form a symbiotic association with its host nematode, Heterorhabditis bacteriophora.

CONCLUSIONS

We studied a highly abundant protein, Pam, which is secreted in a temperature-dependent manner in P. asymbiotica. Our findings indicate that Pam plays an important role in enhancing surface attachment in insect blood. Its association with exopolysaccharide suggests it may exert its effect through mediation of EPS properties. Despite its abundance and conservation in the genus, we find no evidence for a role of Pam in either virulence or symbiosis.

Effective approach for calculations of absolute stability of proteins using focused dielectric constants

The ability to predict the absolute stability of proteins based on their corresponding sequence and structure is a problem of great fundamental and practical importance. In this work, we report an extensive, refinement and validation of our recent approach (Roca et al., FEBS Lett 2007;581:2065-2071) for predicting absolute values of protein stability DeltaG(fold). This approach employs the semimacroscopic protein dipole Langevin dipole method in its linear response approximation version (PDLD/S-LRA) while using the best fitted values of the dielectric constants epsilon'(p) and epsilon'(eff) for the self energy and charge-charge interactions, respectively. The method is validated on a diverse set of 45 proteins. It is found that the best fitted values of both dielectric constants are around 40. However, the self energy of internal residues and the charge-charge interactions of Lys have to be treated with care, using a somewhat lower values of epsilon'(p) and epsilon'(eff). The predictions of DeltaG(fold) reported here, have an average error of only 1.8 kcal/mole compared to the observed values, making our method very promising for estimating protein stability. It also provides valuable insight into the complex electrostatic phenomena taking place in folded proteins.

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.

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.

Unique fluorophores in the dimeric archaeal histones hMfB and hPyA1 reveal the impact of nonnative structure in a monomeric kinetic intermediate

Homodimeric archaeal histones and heterodimeric eukaryotic histones share a conserved structure but fold through different kinetic mechanisms, with a correlation between faster folding/association rates and the population of kinetic intermediates. Wild-type hMfB (from Methanothermus fervidus) has no intrinsic fluorophores; Met35, which is Tyr in hyperthermophilic archaeal histones such as hPyA1 (from Pyrococcus strain GB-3A), was mutated to Tyr and Trp. Two Tyr-to-Trp mutants of hPyA1 were also characterized. All fluorophores were introduced into the long, central alpha-helix of the histone fold. Far-UV circular dichroism (CD) indicated that the fluorophores did not significantly alter the helical content of the histones. The equilibrium unfolding transitions of the histone variants were two-state, reversible processes, with DeltaG degrees (H2O) values within 1 kcal/mol of the wild-type dimers. The hPyA1 Trp variants fold by two-state kinetic mechanisms like wild-type hPyA1, but with increased folding and unfolding rates, suggesting that the mutated residues (Tyr-32 and Tyr-36) contribute to transition state structure. Like wild-type hMfB, M35Y and M35W hMfB fold by a three-state mechanism, with a stopped-flow CD burst-phase monomeric intermediate. The M35 mutants populate monomeric intermediates with increased secondary structure and stability but exhibit decreased folding rates; this suggests that nonnative interactions occur from burial of the hydrophobic Tyr and Trp residues in this kinetic intermediate. These results implicate the long central helix as a key component of the structure in the kinetic monomeric intermediates of hMfB as well as the dimerization transition state in the folding of hPyA1.

Thermodynamic stability and folding of proteins from hyperthermophilic organisms

Life grows almost everywhere on earth, including in extreme environments and under harsh conditions. Organisms adapted to high temperatures are called thermophiles (growth temperature 45-75 degrees C) and hyperthermophiles (growth temperature >or= 80 degrees C). Proteins from such organisms usually show extreme thermal stability, despite having folded structures very similar to their mesostable counterparts. Here, we summarize the current data on thermodynamic and kinetic folding/unfolding behaviors of proteins from hyperthermophilic microorganisms. In contrast to thermostable proteins, rather few (i.e. less than 20) hyperthermostable proteins have been thoroughly characterized in terms of their in vitro folding processes and their thermodynamic stability profiles. Examples that will be discussed include co-chaperonin proteins, iron-sulfur-cluster proteins, and DNA-binding proteins from hyperthermophilic bacteria (i.e. Aquifex and Theromotoga) and archea (e.g. Pyrococcus, Thermococcus, Methanothermus and Sulfolobus). Despite the small set of studied systems, it is clear that super-slow protein unfolding is a dominant strategy to allow these proteins to function at extreme temperatures.

In vivo characterization of thermal stabilities of Aeropyrum pernix cellular components by differential scanning calorimetry

Revival studies of Aeropyrum pernix show that the viability of cells and cell recovery after heat treatment depends on the temperature of treatment. Differential scanning calorimetry (DSC) is used to analyze the relative thermal stabilities of cellular components of A. pernix and to identify the cellular components responsible for the observed lag phase and reduced maximum growth following a heat treatment. DSC thermograms show 5 visible endothermic transitions with 2 major transitions. DSC analysis of isolated crude ribosomes aids the assignment of the 2 major peaks observed in whole-cell thermograms to denaturation of ribosomal structures. A comparison of partial and immediate full rescan thermograms of A. pernix whole cells indicates that both major peaks represent irreversible thermal transitions. A DNA peak is also identified in the whole-cell thermogram by comparison with the optical data of isolated pure DNA. DNA melting is shown to be irreversible in dilute solution, whereas it is partially reversible in whole cells, owing at least in part, to restricted volume effects. In contrast to mesophilic organisms, hyperthermophilic A. pernix ribosomes are more thermally stable than DNA, but in both organisms, irreversible changes leading to cell death occur owing to ribosomal denaturation.

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.

Thermal adaptation of the yeast mitochondrial Hsp70 system is regulated by the reversible unfolding of its nucleotide exchange factor

The Hsp70 protein switches during its functional cycle from an ADP-bound state with a high affinity for substrates to a low-affinity, ATP-bound state, with concomitant release of the client protein. The rate of the chaperone cycle is regulated by co-chaperones such as nucleotide exchange factors that significantly accelerate the ADP/ATP exchange. Mge1p, a mitochondrial matrix protein with homology to bacterial GrpE, serves as the nucleotide exchange factor of mitochondrial Hsp70. Here, we analyze the influence of temperature on the structure and functional properties of Mge1p from the yeast Saccharomyces cerevisiae. Mge1p is a dimer in solution that undergoes a reversible thermal transition at heat-shock temperatures, i.e. above 37 degrees C, that involves protein unfolding and dimer dissociation. The thermally denatured protein is unable to interact stably with mitochondrial Hsp70, and therefore is unable to regulate its ATPase and chaperone cycle. Crosslinking of wild-type mitochondria reveals that Mge1p undergoes the same dimer to monomer temperature-dependent shift, and that the nucleotide exchange factor does not associate with its Hsp70 partner at stress temperatures (i.e. > or =45 degrees C). Once the stress conditions disappear, Mge1p refolds and recovers both structure and functional properties. Therefore, Mge1p can act as a thermosensor for the mitochondrial Hsp70 system, regulating the nucleotide exchange rates under heat shock, as has been described for two bacterial GrpE proteins. The thermosensor activity is conserved in the GrpE-like nucleotide exchange factors although, as discussed here, it is achieved through a different structural mechanism.

The Physics of the Interactions Governing Folding and Association of Proteins

The review discusses the molecular origins of the forces and free energies that determine several things about proteins, and how experiment and theory reveal this information. The first subject is the stability of the folded, native structures. The second is the range of molecular mechanisms by which proteins find their way to those folded structures in laboratory environments. The third is the much more complex problem of how folding occurs in the cellular environment. This topic includes a discussion of crowding and of the roles of chaperone molecules. The review concludes with a discussion of protein aggregation and fibril formation and of misfolding and therapies associated with it.

Additivity in Both Thermodynamic Stability and Thermal Transition Temperature for Rubredoxin Chimeras via Hybrid Native Partitioning

Given any operational definition of pairwise interaction, the set of residues that differ between two structurally homologous proteins can be uniquely partitioned into subsets of clusters for which no such interactions occur between clusters. Although hybrid protein sequences that preserve such clustering are consistent with tertiary structures composed of only parental native-like interactions, the stability of such predicted structures will depend upon the physical robustness of the assumed interaction potential. A simple distance cutoff criterion was applied to the most thermostable protein known to predict such a seven-residue cluster in the metal binding site region of Pyrococcus furiosus rubredoxin and a mesophile homolog. Both conformational stability and thermal transition temperature measurements demonstrate that 39% of the differential stability arises from these seven residues.

Structural features of thermozymes

Enzymes synthesized by thermophiles and hyperthermophiles are known as thermozymes. These enzymes are typically thermostable, or resistant to irreversible inactivation at high temperatures, and thermophilic, i.e. optimally active at elevated temperatures between 60 and 125 degrees C. Enzyme thermostability encompasses thermodynamic stability and kinetic stability. Thermodynamic stability is defined by the enzyme's free energy of stabilization (deltaG(stab)) and by its melting temperature (Tm). An enzyme's kinetic stability is often expressed as its halflife (t1/2) at defined temperature. DeltaG(stab) of thermophilic proteins is 5-20 kcal/mol higher than that of mesophilic proteins. The thermostability mechanisms for thermozymes are varied and depend on the enzyme; nevertheless, some common features can be identified as contributing to stability. These features include more interactions (i.e. hydrogen bonds, electrostatic interactions, hydrophobic interactions, disulfide bonds, metal binding) than in less stable enzymes and superior conformational structure (i.e. more rigid, higher packing efficiency, reduced entropy of unfolding, conformational strain release and stability of alpha-helix). Understanding of the stabilizing features will greatly facilitate reengineering of some of the mesozymes to more stable thermozymes.

The importance of being dimeric

Why are there so many dimeric proteins and enzymes? While for heterodimers a functional explanation seems quite reasonable, the case of homodimers is more puzzling. The number of homodimers found in all living organisms is rapidly increasing. A thorough inspection of the structural data from the available literature and stability (measured from denaturation-renaturation experiments) allows one to suggest that homodimers can be divided into three main types according to their mass and the presence of a (relatively) stable monomeric intermediate in the folding-unfolding pathway. Among other explanations, we propose that an essential advantage for a protein being dimeric may be the proper and rapid assembly in the cellular milieu.

A thermodynamic study of mesophilic, thermophilic, and hyperthermophilic L-arabinose isomerases: the effects of divalent metal ions on protein stability at elevated temperatures

To gain insight into the structural stability of homologous homo-tetrameric l-arabinose isomerases (AI), we have examined the isothermal guanidine hydrochloride (GdnHCl)-induced unfolding of AIs from mesophilic Bacillus halodurans (BHAI), thermophilic Geobacillus stearothermophilus (GSAI), and hyperthermophilic Thermotoga maritima (TMAI) using circular dichroism spectroscopy. The GdnHCl-induced unfolding of the AIs can be well described by a two-state reaction between native tetramers and unfolded monomers, which directly confirms the validity of the linear extrapolation method to obtain the intrinsic stabilities of these proteins. The resulting unfolding free energy (DeltaGU) values of the AIs as a function of temperature were fit to the Gibbs-Helmholtz equation to determine their thermodynamic parameters based on a two-state mechanism. Compared with the stability curves of BHAI in the presence and absence of Mn2+, those of holo GSAI and TMAI were more broadened than those of the apo enzymes at all temperatures, indicating increased melting temperatures (Tm) due to decreased heat capacity (DeltaGp). Moreover, the extent of difference in DeltaCp between the apo and holo thermophilic AIs is larger than that of BHAI. From these studies, we suggest that the metal dependence of the thermophilic AIs, resulting in the reduced DeltaCp, may play a significant role in structural stability compared to their mesophilic analogues, and that the extent of metal dependence of AI stability seems to be highly correlated to oligomerization.

Experiment-guided thermodynamic simulations on reversible two-state proteins: implications for protein thermostability

Here, we perform protein thermodynamic simulations within a set of boundary conditions, effectively blanketing the experimental data. The thermodynamic parameters, melting temperature (TG), enthalpy change at the melting temperature (DeltaHG) and heat capacity change (DeltaCp) were systematically varied over the experimentally observed ranges for small single domain reversible two-state proteins. Parameter sets that satisfy the Gibbs-Helmholtz equation and yield a temperature of maximal stability (TS) around room temperature were selected. The results were divided into three categories by arbitrarily chosen TG ranges. The TG ranges in these categories correspond to typical values of the melting temperatures observed for the majority of the proteins from mesophilic, thermophilic and hyperthermophilic organisms. As expected, DeltaCp values tend to be high in mesophiles and low in hyperthermophiles. An increase in TG is accompanied by an up-shift and broadening of the protein stability curves, however, with a large scatter. Furthermore, the simulations reveal that the average DeltaHG increases with TG up to approximately 360 K and becomes constant thereafter. DeltaCp decreases with TG with different rates before and after approximately 360 K. This provides further justification for the separate grouping of proteins into thermophiles and hyperthermophiles to assess their thermodynamic differences. This analysis of the Gibbs-Helmholtz equation has allowed us to study the interdependence of the thermodynamic parameters TG, DeltaHG and DeltaCp and their derivatives in a more rigorous way than possible by the limited experimental protein thermodynamics data available in the literature. The results provide new insights into protein thermostability and suggest potential strategies for its manipulation.

Stability and folding mechanism of mesophilic, thermophilic and hyperthermophilic archael histones: the importance of folding intermediates

The equilibrium stabilities to guanidinium chloride (GdmCl)-induced denaturation and kinetic folding mechanisms have been characterized for three archael histones: hFoB from the mesophile Methanobacterium formicicum; hMfB from the thermophile Methanothermus fervidus; and hPyA1 from the hyperthermophile Pyrococcus strain GB-3a. These histones are homodimers of 67 to 69 residues per monomer. The equilibrium unfolding transitions, as measured by far-UV circular dichroism (CD) are highly reversible, two-state processes. The mesophilic hFoB is very unstable and requires approximately 1 M trimethyl-amine-N-oxide (TMAO) to completely populate the native state. The thermophilic histones are more stable, with deltaG degrees (H2O) values of 14 and 16 kcal mol(-1) for hMfB and hPyA1, respectively. The kinetic folding of hFoB and hPyA1 are two-state processes, with no detectable transient kinetic intermediates. For hMfB, there is significant development of CD signal in the stopped-flow dead time, indicative of the formation of a monomeric intermediate, which then folds/associates in a single, second-order step to form the native dimer. While the equilibrium stability to chemical denaturation correlates very well with host growth temperature, there is no simple relationship between folding rates and stability for the archael histones. In the absence of denaturant, the log of the unfolding rates correlate with equilibrium stability. The folding/association of the moderately stable hMfB is the most rapid, with a rate constant in the absence of GdmCl of 3 x 10(6) M(-1) s(-1), compared to 9 x 10(5) M(-1) s(-1) for the more stable hPyA1. It appears that the formation of the hMfB burst-phase monomeric ensemble serves to enhance folding efficiency, rather than act as a kinetic trap. The folding mechanism of the archael histones is compared to the folding of other intertwined, segment-swapped, alpha-helical, DNA-binding dimers (ISSADD), including the eukaryotic heterodimeric histones, which fold more rapidly. The importance of monomeric and dimeric kinetic intermediates in accelerating ISSADD folding reactions is discussed.

An electrostatic basis for the stability of thermophilic proteins

Two factors provide key contributions to the stability of thermophilic proteins relative to their mesophilic homologues: electrostatic interactions of charged residues in the folded state and the dielectric response of the folded protein. The dielectric response for proteins in a "thermophilic series" globally modulates the thermal stability of its members, with the calculated dielectric constant for the protein increasing from mesophiles to hyperthermophiles. This variability results from differences in the distribution of charged residues on the surface of the protein, in agreement with structural and genetic observations. Furthermore, the contribution of electrostatic interactions to the stability of the folded state is more favorable for thermophilic proteins than for their mesophilic homologues. This leads to the conclusion that electrostatic interactions play an important role in determining the stability of proteins at high temperatures. The interplay between electrostatic interactions and dielectric response also provides further rationalization for the enhanced stability of thermophilic proteins with respect to cold-denaturation. Taken together, the distribution of charged residues and their fluctuations have been shown to be factors in modulating protein stability over the entire range of biologically relevant temperatures.

Recombinant expression and characterization of an extremely hyperthermophilic archaeal histone from Pyrococcus horikoshii OT3

A histone-like gene, PHS051 from hyperthermophilic archaeon Pyrococcus horikoshii OT3 strain, was cloned, sequenced, and expressed in Escherichia coli. The recombinant histone, HPhA, encodes a protein of 70 amino acids with a molecular weight of 7868Da. Amino acid sequence analysis of HPhA showed high homology with other archaeal histones and eukaryal core histones. The HPhA was purified to homogeneity by heat precipitation and affinity chromatography. Gel electrophoresis mobility shift assays demonstrate that the purified HPhA has high affinity to DNA. The complex of the HPhA and DNA allows DNA to be protected from cleavage by the restriction enzyme TaqI at 65 degrees C. Circular dichroism spectra reveal that the conformation of the recombinant histone HPhA becomes looser when temperatures increase from 25 to 90 degrees C. The HPhA has inherited a remarkable thermostability especially in the presence of 1M KCl and retained DNA binding activity at extreme temperature, which is consistent with our previous report about its structure stability analyzed by X-ray crystallography.

Phenomenological similarities between protein denaturation and small-molecule dissolution: Insights into the mechanism driving the thermal resistance of globular proteins

This article shows that the stability profiles of thermophilic proteins are significantly displaced toward higher temperatures as compared to those of mesophilic proteins. A similar trend characterizes the aqueous transfer of N-alkyl amides. In fact, as a general feature of transfer processes, liquid dissolution profiles are centered at temperatures higher than those of solid ones. This behavior is governed by packing contributions. A partition of the unfolding thermodynamics based on the analysis of phenomenological temperatures common to dissolution and unfolding phenomena provides a clue to understanding the mechanism of thermal stabilization. In fact, the position of stability profiles along the temperature axis does not appear to depend on solvation of internal residues. Instead, it is notably affected by solidlike components, whose progressive decrease appears to drive the heat denaturation temperature increase of most thermostable proteins. As a corollary, it is shown that there are actually two limiting mechanisms of thermal stabilization.

Role of residual structure in the unfolded state of a thermophilic protein

Ribonucleases H from the thermophilic bacterium Thermus thermophilus and the mesophile Escherichia coli demonstrate a dramatic and surprising difference in their change in heat capacity upon unfolding (DeltaCp degrees ). The lower DeltaCp degrees of the thermophilic protein directly contributes to its higher thermal denaturation temperature (Tm). We propose that this DeltaCp degrees difference originates from residual structure in the unfolded state of the thermophilic protein; we verify this hypothesis by using a mutagenic approach. Residual structure in the unfolded state may provide a mechanism for balancing a high Tm with the optimal thermodynamic stability for a protein's function. Structure in the unfolded state is shown to differentially affect the thermodynamic profiles of thermophilic and mesophilic proteins.

Structure based hyperthermostability of archaeal histone HPhA from Pyrococcus horikoshii

The histone protein HPhA from the hyperthermophilic archaeon Pyrococcus horikoshii, shows hyperthermostability, as required for optimal growth of the organism at 95 degrees C. The structure of recombinant P.horikoshii HPhA has been determined to 2.3A resolution by molecular replacement, and refined to R(work) and R(free) values of 20.7% and 27.3%, respectively. The HPhA monomer structure is characterized by the histone fold and assembles into a homodimer like other archaeal histones. There are four anions found in the dimer structure, giving rise to clues as to where DNA might bind. A detailed comparison of four known structures of archaeal histones, which evolve and exist under different temperatures, shows that the thermophilic archaeal histone HPhA has a larger hydrophobic contact area, an increased number of hydrogen bonds and a reduced solvent-accessible area. We also observe a unique network of tyrosine residues located at the crossover point of the two HPhA monomers, which locks them together and stabilizes the dimer. These factors together account for the increased thermal stability.

Toward the physical basis of thermophilic proteins: linking of enriched polar interactions and reduced heat capacity of unfolding

The enrichment of salt bridges and hydrogen bonding in thermophilic proteins has long been recognized. Another tendency, featuring lower heat capacity of unfolding (DeltaC(p)) than found in mesophilic proteins, is emerging from the recent literature. Here we present a simple electrostatic model to illustrate that formation of a salt-bridge or hydrogen-bonding network around an ionized group in the folded state leads to increased folding stability and decreased DeltaC(p). We thus suggest that the reduced DeltaC(p) of thermophilic proteins could partly be attributed to enriched polar interactions. A reduced DeltaC(p) might serve as an indicator for the contribution of polar interactions to folding stability.

Archaeal histone tetramerization determines DNA affinity and the direction of DNA supercoiling

DNA binding and the topology of DNA have been determined in complexes formed by >20 archaeal histone variants and archaeal histone dimer fusions with residue replacements at sites responsible for histone fold dimer:dimer interactions. Almost all of these variants have decreased affinity for DNA. They have also lost the flexibility of the wild type archaeal histones to wrap DNA into a negative or positive supercoil depending on the salt environment; they wrap DNA into positive supercoils under all salt conditions. The histone folds of the archaeal histones, HMfA and HMfB, from Methanothermus fervidus are almost identical, but (HMfA)(2) and (HMfB)(2) homodimers assemble into tetramers with sequence-dependent differences in DNA affinity. By construction and mutagenesis of HMfA+HMfB and HMfB+HMfA histone dimer fusions, the structure formed at the histone dimer:dimer interface within an archaeal histone tetramer has been shown to determine this difference in DNA affinity. Therefore, by regulating the assembly of different archaeal histone dimers into tetramers that have different sequence affinities, the assembly of archaeal histone-DNA complexes could be localized and used to regulate gene expression.

Both DNA and histone fold sequences contribute to archaeal nucleosome stability

The roles and interdependence of DNA sequence and archaeal histone fold structure in determining archaeal nucleosome stability and positioning have been determined and quantitated. The presence of four tandem copies of TTTAAAGCCG in the polylinker region of pLITMUS28 resulted in a DNA molecule with increased affinity (DeltaDeltaG of approximately 700 cal mol(-1)) for the archaeal histone HMfB relative to the polylinker sequence, and the dominant, quantitative contribution of the helical repeats of the dinucleotide TA to this increased affinity has been established. The rotational and translational positioning of archaeal nucleosomes assembled on the (TTTAAAGCCG)(4) sequence and on DNA molecules selectively incorporated into archaeal nucleosomes by HMfB have been determined. Alternating A/T- and G/C-rich regions were located where the minor and major grooves, respectively, sequentially faced the archaeal nucleosome core, and identical positioning results were obtained using HMfA, a closely related archaeal histone also from Methanothermus fervidus. However, HMfA did not have similarly high affinities for the HMfB-selected DNA molecules, and domain-swap experiments have shown that this difference in affinity is determined by residue differences in the C-terminal region of alpha-helix 3 of the histone fold, a region that is not expected to directly interact with DNA. Rather this region is thought to participate in forming the histone dimer:dimer interface at the center of an archaeal nucleosome histone tetramer core. If differences in this interface do result in archaeal histone cores with different sequence preferences, then the assembly of alternative archaeal nucleosome tetramer cores could provide an unanticipated and novel structural mechanism to regulate gene expression.

Comparative analyses of the conformational stability of a hyperthermophilic protein and its mesophilic counterpart

Comparison of the conformational stability of an O(6)-methylguanine-DNA methyltransferase (MGMT) from the hyperthermophilic archaeon Thermococcus kodakaraensis strain KOD1 (Tk-MGMT), and its mesophilic counterpart C-terminal Ada protein from Escherichia coli (Ec-AdaC) was performed in order to obtain information about the relationship between thermal stability and other factors, such as thermodynamic parameters, thermodynamic stability and other unfolding conditions. Tk-MGMT unfolded at Tm = 98.6 degrees C, which was 54.8 degrees C higher than the unfolding temperature of Ec-AdaC. The maximum free energy (DeltaG(max)) of the proteins were different; the value of Tk-MGMT (42.9 kJ.mol-1 at 29.5 degrees C) was 2.6 times higher than that of Ec-AdaC (16.6 kJ.mol-1 at 7.4 degrees C). The high conformational stability of Tk-MGMT was attributed to a 1.6-fold higher enthalpy value than that of Ec-AdaC. In addition, the DeltaG(max) temperature of Tk-MGMT was considerably higher (by 22.1 degrees C). The apparent heat capacity of denaturation (DeltaC(p)) of Tk-MGMT was 0.7-fold lower than that of Ec-AdaC. These three synergistic effects, increasing DeltaGmax, shifted DeltaG vs. temperature curve, and low DeltaC(p), give Tk-MGMT its thermal stability. Unfolding profiles of the two proteins, tested with four alcohols and three denaturants, showed that Tk-MGMT possessed higher stability than Ec-AdaC in all conditions studied. These results indicate that the high stability of Tk-MGMT gives resistance to chemical unfolding, in addition to thermal unfolding.

Thermodynamics of trimer-of-hairpins formation by the SIV gp41 envelope protein

The gp41 envelope protein mediates the entry of primate immunodeficiency viruses into target cells by promoting the fusion of viral and cellular membranes. The structure of the gp41 ectodomain core represents a trimer of identical helical hairpins in which a central trimeric coiled-coil made up of three amino-terminal helices is wrapped in an outer layer of three antiparallel carboxyl-terminal helices. Triggering formation of this fusion-active gp41 conformation appears to cause close membrane apposition and thus overcome the activation energy barrier for lipid bilayer fusion. We present a detailed description of the folding thermodynamics of the simian immunodeficiency virus (SIV) gp41 core by using a recombinant trimeric N34(L6)C28 model. Differential scanning calorimetry and spectroscopic experiments on denaturant-induced and thermal unfolding indicate that the free energy of association of three 68 residue N34(L6)C28 peptides to a trimer-of-hairpins is 76 kJ mol(-1) at pH 7.0 and 25 degrees C in physiological buffer. The associated enthalpy change, Delta H(unf), is 177 kJ mol(-1), while the entropy of unfolding, Delta S(unf), is 0.32 kJ K(-1) mol(-1). The temperature of maximal stability is close to 20 degrees C. The unfolding heat capacity increment is approximately 9 kJ K(-1) mol(-1) (approximately 45 J K(-1) mol residue(-1)), which is lower than expected for unfolding of the trimer to an extended and fully hydrated polypeptide chain. Replacement by isoleucine of the polar residues Thr582 or Thr586 buried in the N-terminal trimeric coiled-coil interface leads to very strong stabilization of the trimer-of-hairpins, 30-35 kJ mol(-1). Single-point mutations in the central coiled-coil thus strongly stabilize the gp41 core structure. These thermodynamic characteristics may be important for the refolding of the gp41 envelope protein into its fusion-active conformation during membrane fusion.

Hyperthermophilic enzymes: sources, uses, and molecular mechanisms for thermostability

Enzymes synthesized by hyperthermophiles (bacteria and archaea with optimal growth temperatures of > 80 degrees C), also called hyperthermophilic enzymes, are typically thermostable (i.e., resistant to irreversible inactivation at high temperatures) and are optimally active at high temperatures. These enzymes share the same catalytic mechanisms with their mesophilic counterparts. When cloned and expressed in mesophilic hosts, hyperthermophilic enzymes usually retain their thermal properties, indicating that these properties are genetically encoded. Sequence alignments, amino acid content comparisons, crystal structure comparisons, and mutagenesis experiments indicate that hyperthermophilic enzymes are, indeed, very similar to their mesophilic homologues. No single mechanism is responsible for the remarkable stability of hyperthermophilic enzymes. Increased thermostability must be found, instead, in a small number of highly specific alterations that often do not obey any obvious traffic rules. After briefly discussing the diversity of hyperthermophilic organisms, this review concentrates on the remarkable thermostability of their enzymes. The biochemical and molecular properties of hyperthermophilic enzymes are described. Mechanisms responsible for protein inactivation are reviewed. The molecular mechanisms involved in protein thermostabilization are discussed, including ion pairs, hydrogen bonds, hydrophobic interactions, disulfide bridges, packing, decrease of the entropy of unfolding, and intersubunit interactions. Finally, current uses and potential applications of thermophilic and hyperthermophilic enzymes as research reagents and as catalysts for industrial processes are described.

Observation of the noncovalent assembly and disassembly pathways of the chaperone complex MtGimC by mass spectrometry

We have analyzed a newly described archaeal GimC/prefoldin homologue, termed MtGimC, by using nanoflow electrospray coupled with time-of-flight MS. The molecular weight of the complex from Methanobacterium thermoautotrophicum corresponds to a well-defined hexamer of two alpha subunits and four beta subunits. Dissociation of the complex within the gas phase reveals a quaternary arrangement of two central subunits, both alpha, and four peripheral beta subunits. By constructing a thermally controlled nanoflow device, we have monitored the thermal stability of the complex by MS. The results of these experiments demonstrate that a significant proportion of the MtGimC hexamer remains intact under low-salt conditions at elevated temperatures. This finding is supported by data from CD spectroscopy, which show that at physiological salt concentrations, the complex remains stable at temperatures above 65 degrees C. Mass spectrometric methods were developed to monitor in real time the assembly of the MtGimC hexamer from its component subunits. By using this methodology, the mass spectra recorded throughout the time course of the experiment showed the absence of any significantly populated intermediates, demonstrating that the assembly process is highly cooperative. Taken together, these data show that the complex is stable under the elevated temperatures that are appropriate for its hyperthermophile host and demonstrate that the assembly pathway leads exclusively to the hexamer, which is likely to be a structural unit in vivo.

Crystal structures of recombinant histones HMfA and HMfB from the hyperthermophilic archaeon Methanothermus fervidus

The hyperthermophilic archaeon Methanothermus fervidus contains two small basic proteins, HMfA (68 amino acid residues) and HMfB (69 residues) that share a common ancestry with the eukaryal nucleosome core histones H2A, H2B, H3, and H4. HMfA and HMfB have sequences that differ at 11 locations, they have different structural stabilities, and the complexes that they form with DNA have different electrophoretic mobilities. Here, crystal structures are documented for recombinant (r) HMfA at a resolution of 1.55 A refined to a crystallographic R-value of 19.8 % (tetragonal form) and at 1.48 A refined to a R-value of 18.8 % (orthorhombic form), and for rHMfB at 1.9 A refined to a R-value of 18.0 %. The rHMfA and rHMfB monomers have structures that are just histone folds in which a long central alpha-helix (alpha2; 29 residues) is separated from shorter N-terminal (alpha1; 11 residues) and C-terminal (alpha3; 10 residues) alpha-helices by two loops (L1 and L2; both 6 residues). Within L1 and L2, three adjacent residues are in extended (beta) conformation. rHMfA and rHMfB assemble into homodimers, with the alpha2 helices anti-parallel aligned and crossing at an angle of close to 35 degrees, and with hydrogen bonds formed between the extended, parallel regions of L1 and L2 resulting in short beta-ladders. Dimerization creates a novel N-terminal structure that contains four proline residues, two from each monomer. As prolines are present at these positions in all archaeal histone sequences, this proline-tetrad structure is likely to be a common feature of all archaeal histone dimers. Almost all residues that participate in monomer-monomer interactions are conserved in HMfA and HMfB, consistent with the ability of these monomers to form both homodimers and (HMfA+HMfB) heterodimers. Differences in side-chain interactions that result from non-conservative residue differences in HMfA and HMfB are identified, and the structure of a (rHMfA)(2)-DNA complex is presented based on the structures documented here and modeled by homology to histone-DNA interactions in the eukaryal nucleosome.

Mutational analysis of archaeal histone-DNA interactions

Site-specific mutagenesis of the hmfB gene cloned from the archaeon Methanothermus fervidus, followed by expression in Escherichia coli, has been used to generate approximately 90 recombinant (r) variants of the archaeal histone HMfB. The abilities of these variants to form stable archaeal nucleosome-containing complexes with linear pBR322 DNA, and with an 89 bp restriction fragment of this DNA have been determined. Variants that failed to form such complexes, based on negative gel-shift assays, had substitutions at the N terminus or within the alpha1, L1 and L2 regions of the rHMfB histone fold, at sites predicted to be homologous to eucaryal histone fold residues that contact the DNA in the eucaryal nucleosome. Variants that failed to give gel shifts were further assayed for their abilities to facilitate ligase-catalyzed circularization of a linear 88 bp DNA molecule, and to reduce the ellipticity of a DNA solution at 275 nm (theta(275)). Consistent with cooperative but independent sites of DNA binding, a combination of three residue substitutions, one each in alpha1, L1 and L2, was required to generate a rHMfB variant with no detectable DNA binding based on gel shift, circularization and theta(275) reduction assays.