Protein stability and molecular adaptation to extreme conditions.

Proteins, due to the delicate balance of stabilizing and destabilizing interactions, are only marginally stable. Adaptation to extreme environments tends to shift the 'mesophilic' characteristics of proteins to the respective extremes of temperature, hydrostatic pressure, pH and salinity, such that, under the mutual physiological conditions, the molecular properties are similar regarding overall topology, flexibility and solvation. Enhanced intrinsic stability requires only minute local structural changes so that general strategies of stabilization cannot be established. Apart from mutative changes of amino-acid sequences, extrinsic factors (or cellular components) may be involved in 'extremophilic adaptation'. The molecular basis of acidophilic, alkalophilic and barophilic adaptation is still obscure. Mechanisms of enhanced thermal stability involve improved packing density, as well as specific local interactions. In halophiles, water and salt binding of the intrinsically stable protein inventory is accomplished by favoring acidic over basic amino acid residues and decreased hydrophobicity. General limits of viability are: (a) the susceptibility of the covalent structure of the polypeptide chain toward hydrolysis or hydrothermal degradation; (b) the competition of extreme solvent parameters with the weak electrostatic and hydrophobic interactions involved in protein stabilization; (c) perturbations of the folding and assembly of proteins; and (d) 'dislocation' of biochemical pathways due to effects of extreme conditions on the intricate network of metabolic reactions.     

Stability and stabilization of globular proteins in solution

Proteins are multifunctional: their amino acid sequences simultaneously determine folding, function and turnover. Correspondingly, evolution selected for compromises between rigidity (stability) and flexibility (folding/function/degradation), to the result that generally the free energy of stabilization of globular proteins in solution is the equivalent to only a few weak intermolecular interactions. Additional increments may come from extrinsic factors such as ligands or specific compatible solutes. Apart from the enthalpic effects, entropy may play a role by reducing the flexibility (cystine bridges, increased proline content), or by water release from residues buried upon folding and association. Additional quaternary interactions and closer packing are typical characteristics of proteins from thermophiles. In halophiles, protein stability and function are maintained by increased ion binding and glutamic acid content, both allowing the protein inventory to compete for water at high salt. Acidophiles and alkalophiles show neutral intracellular pH; proteins facing the outside extremes of pH possess anomalously high contents in ionizable amino acids. Global comparisons of the amino acid compositions and sequences of proteins from mesophiles and extremophiles did not result in general rules of protein stabilization, even after including complete genome sequences into the search. Obviously, proteins are individuals that optimize internal packing and external solvent interactions by very different mechanisms, each protein in its own way. Strategies deduced from specific ultrastable proteins allow stabilizing point mutations to be predicted.     

Proteins under extreme physical conditions

Life on earth is ubiquitous within the limits from -5 to 110 degrees C for temperature, 0.1 to 120 MPa for hydrostatic pressure, 1.0 to 0.6 for water activity and pH 1 to 12. In general, mutative adaptation of proteins to changing environmental conditions tends to maintain 'corresponding states' regarding overall topology, flexibility and hydration. Due to the minute changes in the free energy of stabilization responsible for enhanced stability, nature provides a wide variety of different adaptative strategies. In the case of thermophilic proteins, improved packing densities are crucial. In halophilic proteins, decreased hydrophobicity and clustered surface charges serve to increase water and salt binding required for solubilization at high salt concentration. In the case of barophiles, high-pressure adaptation is expected to be less important than adaptation to low temperatures governing the deep sea. Nothing is known with respect to the mechanisms underlying psychrophilic and acidophilic/alkalophilic adaptation.     

Review: Protein function at thermal extremes: balancing stability and flexibility.

No organism can survive across the entire temperature range found in the biosphere, and a given species can rarely support active metabolism across more than a few tens of degrees C. Nevertheless, life can be maintained at surprisingly extreme temperatures, from below -50 to over 110 degrees C. That proteins, which are assembled with the same 20 amino acids in all species, can function well at both extremes of this range illustrates the plasticity available in the construction of these macromolecules. In studying proteins from extremophiles, researchers have found no new amino acids, covalent modifications or structural motifs that explain the ability of these molecules to function in such harsh environments. Rather, subtle redistributions of the same intramolecular interactions required for protein stabilization at moderate temperatures are sufficient to maintain structural integrity at hot or cold extremes. The key to protein function, whether in polar seas or hot springs, is the maintenance of an appropriate balance between molecular stability on the one hand and structural flexibility on the other. Stability is needed to ensure the appropriate geometry for ligand binding, as well as to avoid denaturation, while flexibility is necessary to allow catalysis at a metabolically appropriate rate. Comparisons of homologous proteins from organisms spanning a wide range of thermal habitats show that adaptive mutations, as well as stabilizing solutes, maintain a balance between these two attributes, regardless of the temperature at which the protein functions.     

Cold adaptation of enzyme reaction rates

A major issue for organisms living at extreme temperatures is to preserve both stability and activity of their enzymes. Cold-adapted enzymes generally have a reduced thermal stability, to counteract freezing, and show a lower enthalpy and a more negative entropy of activation compared to mesophilic and thermophilic homologues. Such a balance of thermodynamic activation parameters can make the reaction rate decrease more linearly, rather than exponentially, as the temperature is lowered, but the structural basis for rate optimization toward low working temperatures remains unclear. In order to computationally address this problem, it is clear that reaction simulations rather than standard molecular dynamics calculations are needed. We have thus carried out extensive computer simulations of the keto-enol(ate) isomerization steps in differently adapted citrate synthases to explore the structure-function relationships behind catalytic rate adaptation to different temperatures. The calculations reproduce the absolute rates of the psychrophilic and mesophilic enzymes at 300 K, as well as the lower enthalpy and more negative entropy of activation of the cold-adapted enzyme, where the latter simulation result is obtained from high-precision Arrhenius plots. The overall catalytic effect originates from electrostatic stabilization of the transition state and enolate and the reduction of reorganization free energy. The simulations, however, show psychrophilic, mesophilic, and hyperthermophilic citrate synthases to have increasingly stronger electrostatic stabilization of the transition state, while the energetic penalty in terms of internal protein interactions follows the reverse order with the cold-adapted enzyme having the most favorable energy term. The lower activation enthalpy and more negative activation entropy observed for cold-adapted enzymes are found to be associated with a decreased protein stiffness. The origin of this effect is, however, not localized to the active site but to other regions of the protein structure.     

Dynamic arrangement of ion pairs and individual contributions to the thermal stability of the cofactor-binding domain of glutamate dehydrogenase from Thermotoga maritima

The dynamics of a hyperthermophilic protein fragment in a water environment, as studied by performing molecular dynamics (MD) simulations at various temperatures, is compared to the dynamical behavior of a homologous mesophilic protein simulated under identical conditions. The effects on the stability of the spatial arrangement and mobility of the charged residues in solution were quantified by calculating free energy changes upon salt bridge formation in these proteins. Electrostatic free energy terms derived from a thermodynamic cycle were obtained by solving the linearized Poisson-Boltzmann equation for a series of protein conformations generated by MD simulations and placed subsequently in a continuum solvent medium. Our results show that the ion pairs are electrostatically stabilizing in most of the cases, but their individual contributions vary significantly. The greater contribution of the charged residues to the stability of the hyperthermophilic protein as compared with the mesophilic counterpart was evidenced only by the calculations that included conformations sampled at 343 and 373 K. The "dynamic" structure of the hyperthermophilic protein fragment simulated at elevated temperatures reveals an optimum placement of the ionizable residues within the protein structure as well as the role of their cooperative interactions in promoting thermal stability. The thermodynamic properties such as electrostatic free energy differences, configurational entropies, and specific heat capacities calculated in the dynamic context of the protein structure provided new insight into the mechanism of protein thermostabilization.     

Protein stability curves

The stability curve of a protein is defined as the plot of the free energy of unfolding as a function of temperature. For most proteins the change in heat capacity on denaturation, or unfolding, is large but approximately constant. When unfolding is s two-state process, most of the salient features of the stability curves of proteins can be derived from this fact. A number of relations are obtained, including the special features of low-temperature denaturation, the properties of the maximum in stability, and the interrelationships of the characteristic temperatures of the protein. The paper closes with a formula that permits one to calculate small changes in stabilization free energy from changes in the melting temperature of the protein.     

Thermodynamic differences among homologous thermophilic and mesophilic proteins.

Here, we analyze the thermodynamic parameters and their correlations in families containing homologous thermophilic and mesophilic proteins which show reversible two-state folding <--> unfolding transitions between the native and the denatured states. For the proteins in these families, the melting temperatures correlate with the maximal protein stability change (between the native and the denatured states) as well as with the enthalpic and entropic changes at the melting temperature. In contrast, the heat capacity change is uncorrelated with the melting temperature. These and additional results illustrate that higher melting temperatures are largely obtained via an upshift and broadening of the protein stability curves. Both thermophilic and mesophilic proteins are maximally stable around room temperature. However, the maximal stabilities of thermophilic proteins are considerably greater than those of their mesophilic homologues. At the living temperatures of their respective source organisms, homologous thermophilic and mesophilic proteins have similar stabilities. The protein stability at the living temperature of the source organism does not correlate with the living temperature of the protein. We tie thermodynamic observations to microscopics via the hydrophobic effect and a two-state model of the water structure. We conclude that, to achieve higher stability and greater resistance to high and low temperatures, specific interactions, particularly electrostatic, should be engineered into the protein. The effect of these specific interactions is largely reflected in an increased enthalpy change at the melting temperature.     

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.     

Structural stability of binding sites: Consequences for binding affinity and allosteric effects

During the course of biological function, proteins interact with other proteins, ligands, substrates, inhibitors, etc. These interactions occur at precisely defined locations within the protein but their effects are sometimes propagated to distal regions, triggering highly specific responses. These effects can be used as signals directed to activate or inhibit other sites, modulate interactions with other molecules, and/or establish inter‐molecular communication networks. During the past decade, it has become evident that the energy of stabilization of the protein structure is not evenly distributed throughout the molecule and that, under native conditions, proteins lack global cooperativity and are characterized by the occurrence of multiple independent local unfolding events. From a biological point of view, it is important to assess if this uneven distribution reflects specific functional requirements. For example, are binding sites more likely to be found in well structured regions, unstable regions, or mixed regions? In this article, we have addressed these questions by performing a structure‐based thermodynamic stability analysis of non‐structurally homologous proteins for which high resolution structures of their complexes with specific ligands are available. The results of these studies indicate that for all 16 proteins considered, the binding sites have a dual character and are characterized by the presence of regions with very low structural stability and regions with high stability. In many cases the low stability regions are loops that become stable and cover a significant portion of low molecular weight ligands upon binding. For enzymes, catalytic residues are usually, but not always, located in regions with high structural stability. It is shown that this arrangement provides significant advantages for the optimization of binding affinity of small ligands. In allosteric enzymes, low stability regions in the regulatory site are shown to play a crucial role in the transmission of information to the catalytic site. Proteins 2000;41:63–71. © 2000 Wiley‐Liss, Inc.     

Adaptation to high temperatures through macromolecular dynamics by neutron scattering

Work on the relationship between hyperthermophile protein dynamics, stability and activity is reviewed. Neutron spectroscopy has been applied to measure and compare the macromolecular dynamics of various hyperthermophilic and mesophilic proteins, under different conditions. First, molecular dynamics have been analyzed for the hyperthermophile malate dehydrogenase from Methanococcus jannaschii and a mesophilic homologue, the lactate dehydrogenase from Oryctolagus cunniculus (rabbit) muscle. The neutron scattering approach has provided independent measurements of the global flexibility and structural resilience of each protein, and it has been demonstrated that macromolecular dynamics represents one of the molecular mechanisms of thermoadaptation. The resilience was found to be higher for the hyperthermophilic protein, thus ensuring similar flexibilities in both enzymes at their optimal activity temperature. Second, the neutron method has been developed to quantify the average macromolecular flexibility and resilience within the natural crowded environment of the cell, and mean macromolecular motions have been measured in vivo in psychrophile, mesophile, thermophile and hyperthermophile bacteria. The macromolecular resilience in bacteria was found to increase with adaptation to high temperatures, whereas flexibility was maintained within narrow limits, independent of physiological temperature for all cells in their active state. Third, macromolecular motions have been measured in free and immobilized dihydrofolate reductase from Escherichia coli. The immobilized mesophilic enzyme has increased stability and decreased activity, so that its properties are changed to resemble those of a thermophilic enzyme. Quasi-elastic neutron scattering measurements have also been performed to probe the protein motions. Compared to the free enzyme, the average height of the activation free energy barrier to local motions was found to be increased by 0.54 kcal.mol(-1) in the immobilized dihydrofolate reductase, a value that is of the same order as expected from the theoretical rate equation.     

Thermal and Structural Stability of Adsorbed Proteins

Experimental evidence suggests that proteins adsorbed to hydrophobic surfaces at low coverages are stabilized relative to the bulk. For larger coverages, proteins unfold and form β-sheets. We performed computer simulations on model proteins and found that: 1), For weakly adsorbing surfaces, unfolded conformations lose more entropy upon adsorption than folded ones. 2), The melting temperature, both in the bulk and at surfaces, decreases with increasing protein concentration because of favorable interprotein interactions. 3), Proteins in the bulk show large unfolding free energy barriers; this barrier decreases at stronger adsorbing surfaces. We conjecture that typical experimental temperatures appear to be below the bulk melting temperature for a single protein, but above the melting temperature for concentrated protein solutions. Purely thermodynamic factors then explain protein stabilization on adsorption at low concentrations. However, both thermodynamic and kinetic factors are important at higher concentrations. Thus, proteins in the bulk do not denature with increasing concentration due to large kinetic barriers, even though the aggregated state is thermodynamically preferred. However, they readily unfold upon adsorption, with the surface acting as a heterogeneous catalyst. The thermal behavior of proteins adsorbed to hydrophobic surfaces thus appears to follow behavior independent of their chemical specificity.     

Structural investigation of cold activity and regulation of aspartate carbamoyltransferase from the extreme psychrophilic bacterium Moritella profunda

Aspartate carbamoyltransferase (EC 2.1.3.2) is extensively studied as a model for cooperativity and allosteric regulation. The structure of the Escherichia coli enzyme has been thoroughly analyzed by X-ray crystallography, and recently the crystal structures of two hyperthermophilic ATCases of the same structural class have been characterized. We here report the detailed functional and structural investigation of the ATCase from the psychrophilic deep sea bacterium Moritella profunda. Our analysis indicates that the enzyme conforms to the E. coli model in that two allosteric states exist that are influenced by similar homotropic interactions. The heterotropic properties differ in that CTP and UTP inhibit the holoenzyme, but ATP seems to exhibit a dual regulatory pattern, activating the enzyme at low concentrations and inhibiting it in the mM range. The crystal structure of the unliganded M. profunda ATCase shows resemblance to a more extreme T state reported previously for an E. coli ATCase mutant. A detailed molecular analysis reveals potential features of adaptation to cold activity and cold regulation. Moreover, M. profunda ATCase presents similarities with certain mutants of E. coli ATCase altered in their kinetic properties or temperature relationships. Finally, structural and functional comparison of ATCases across the full physiological temperature range agrees with an important, but fundamentally different role for electrostatics in protein adaptation at both extremes, i.e. an increased stability through the formation of ion pairs and ion pair networks at high physiological temperatures, and an increased flexibility through enhanced protein solvation at low temperatures.     

Extreme free energy of stabilization of Taq DNA polymerase

We have examined the chemical denaturations of the Klentaq and Klenow large-fragment domains of the Type 1 DNA polymerases from Thermus aquaticus (Klentaq) and Escherichia coli (Klenow) under identical solution conditions in order to directly compare the stabilization energetics of the two proteins. The high temperature stability of Taq DNA polymerase is common knowledge, and is the basis of its use in the polymerase chain reaction. This study, however, is aimed at understanding the thermodynamic basis for this high-temperature stability. Chemical denaturations with guanidine hydrochloride report a folding free energy (DeltaG) for Klentaq that is over 20 kcal/mol more favorable than that for Klenow under the conditions examined. This difference between the stabilization free energies of a homologous mesophilic-thermophilic protein pair is significantly larger than generally observed. This is due in part to the fact that the stabilization free energy for Klentaq polymerase, at 27.5 kcal/mol, is one of the largest ever determined for a monomeric protein. Large differences in the chemical midpoints of the unfolding (Cm) and the dependences of the unfolding free energy on denaturant concentration in the transition region (m-value) between the two proteins are also observed. Measurements of the sedimentation coefficients of the two proteins in the native and denatured states report that both proteins approximately double in hydrodynamic size upon denaturation, but that Klentaq expands somewhat more than Klenow.     

On the stabilizing action of protein denaturants: acetonitrile effect on stability of lysozyme in aqueous solutions

Stability of hen lysozyme in the presence of acetonitrile (MeCN) at different pH values of the medium was studied by scanning microcalorimetry with a special emphasis on determination of reliable values of the denaturational heat capacity change. It was found that the temperature of denaturation decreases on addition of MeCN. However, the free energy extrapolation showed that below room temperature the thermodynamic stability increases at low concentrations of MeCN in spite of the general destabilizing effect at higher concentrations and temperatures. Charge-induced contribution to this stabilization was shown to be negligible (no pH-dependence was found); therefore, the most probable cause for the phenomenon is an increase of hydrophobic interactions at low temperatures in aqueous solutions containing small amounts of the organic additive. The difference in preferential solvation of native and denatured states of lysozyme was calculated from the stabilization free energy data. It was found that the change in preferential solvation strongly depends on the temperature in the water-rich region. At the higher MeCN content this dependence decreases until, at 0.06 mole fractions of MeCN, the difference in the preferential solvation between native and denatured lysozyme becomes independent of the temperature over a range of 60 K. The importance of taking into account non-ideality of a mixed solution, when analyzing preferential solvation phenomena was emphasized.     

The stability of proteins in extreme environments

Three complete genome sequences of thermophilic bacteria provide a wealth of information challenging current ideas concerning phylogeny and evolution, as well as the determinants of protein stability. Considering known protein structures from extremophiles, it becomes clear that no general conclusions can be drawn regarding adaptive mechanisms to extremes of physical conditions. Proteins are individuals that accumulate increments of stabilization; in thermophiles these come from charge clusters, networks of hydrogen bonds, optimization of packing and hydrophobic interactions, each in its own way. Recent examples indicate ways for the rational design of ultrastable proteins.     

Stability Curve Prediction of Homologous Proteins Using Temperature-Dependent Statistical Potentials

The unraveling and control of protein stability at different temperatures is a fundamental problem in biophysics that is substantially far from being quantitatively and accurately solved, as it requires a precise knowledge of the temperature dependence of amino acid interactions. In this paper we attempt to gain insight into the thermal stability of proteins by designing a tool to predict the full stability curve as a function of the temperature for a set of 45 proteins belonging to 11 homologous families, given their sequence and structure, as well as the melting temperature () and the change in heat capacity () of proteins belonging to the same family. Stability curves constitute a fundamental instrument to analyze in detail the thermal stability and its relation to the thermodynamic stability, and to estimate the enthalpic and entropic contributions to the folding free energy. In summary, our approach for predicting the protein stability curves relies on temperature-dependent statistical potentials derived from three datasets of protein structures with targeted thermal stability properties. Using these potentials, the folding free energies () at three different temperatures were computed for each protein. The Gibbs-Helmholtz equation was then used to predict the protein''s stability curve as the curve that best fits these three points. The results are quite encouraging: the standard deviations between the experimental and predicted ''s, ''s and folding free energies at room temperature () are equal to 13 , 1.3 ) and 4.1 , respectively, in cross-validation. The main sources of error and some further improvements and perspectives are briefly discussed.     

Confinement effects on the thermodynamics of protein folding: Monte Carlo simulations

The effects of chaperonin-like cage-induced confinement on protein stability have been studied for molecules of varying sizes and topologies. Minimalist models based on Gō-like interactions are employed for the proteins, and density-of-states-based Monte Carlo simulations are performed to accurately characterize the thermodynamic transitions. This method permits efficient sampling of conformational space and yields precise estimates of free energy and entropic changes associated with protein folding. We find that confinement-driven stabilization is not only dependent on protein size and cage radius, but also on the specific topology. The choice of the confining potential is also shown to have an effect on the observed stabilization and the scaling behavior of the stabilization with respect to the cage size.     

Low-temperature unfolding of a mutant of phage T4 lysozyme. 1. Equilibrium studies

The mutant protein I3C-C97/C54T of phage T4 lysozyme is free of sulfhydryl groups and has a genetically engineered disulfide bridge between positions 3 and 97 (Perry & Wetzel, 1986). This protein has a maximum stability at 12 degrees C in 3 M guanidinium chloride and undergoes reversible high- and low-temperature melting at 28 and -3 degrees C, respectively, in this medium. The free energy of stabilization of the protein has been studied over a range of temperature that includes both melting transitions. The stability curve fits a constant delta Cp model over the entire range, permitting an unusually complete determination of the thermodynamic parameters of the protein and demonstrating that the low-temperature unfolded form of the protein may be interpreted as an extrapolation with constant delta Cp of the high-temperature unfolded form. The free energy of unfolding is a linear function of guanidinium concentration within experimental error which permits a rough estimate of the stability of the protein at low temperatures and of the differential interaction of the unfolded protein with guanidinium chloride. These equilibrium studies provide a basis for the interpretation of the kinetic studies reported in the following paper.     

Enzymes from psychrophilic organisms

Abstract: Psychrophilic organisms such as micro-organisms and other ectothermic species living in polar, deep- sea or any constantly low temperature environments, produce enzymes adapted to function at low temperature. These enzymes are characterized by a high catalytic efficiency at low and moderate temperatures but are rather thermolabile. Due to their high specific activity and their rapid inactivation at temperatures as low as 30°C, they offer, along with the producing micro-organisms, a great potential in biotechnology. The molecular basis of the adaptation of cold α-amylase, subtilisin, triose phosphate isomerase from Antarctic bacteria and of trypsin from fish living in North Atlantic and in Antarctic sea waters have been studied. The comparison of the 3D structures obtained either by protein modelling or by X-ray crystallography (North Atlantic trypsin) with those of their mesophilic counterparts indicates that the molecular changes tend to increase the flexibility of the structure by a weakening of the intramolecular interactions and by an increase of the interactions with the solvent. For each enzyme, the most appropriate strategy enabling it to accommodate the substrate at a low energy cost is selected. There is a price to pay in terms of thermosensibility because the selective pressure is essentially oriented towards the harmonization of the specific activity with ambient thermal conditions. However, as demonstrated by site-directed mutagenesis experiments carried out on the Antarctic subtilisin, the possibility remains to stabilize the structure of these enzymes without affecting their high catalytic efficiency.