The optimization of protein-solvent interactions: thermostability and the role of hydrophobic and electrostatic interactions

Crystal structure of a thermophilic alcohol dehydrogenase substrate complex suggests determinants of substrate specificity and thermostability

The crystal structure of a thermophilic alcohol dehydrogenase (TBAD) from Thermoanaerobacter brockii has been determined in a binary complex with sec-butanol as substrate to a resolution of 3.0 A. Van der Waals interactions of the carbon C1 atom of sec-butanol with atoms in His59, Ala85, Trp110, Asp150, and Leu294 account for the substrate preference of this enzyme for secondary over primary alcohols. A crevice from the surface to the active site provides access for substrates and products. This opening is lined with the hydrophobic residues Ile49, Leu107, Trp110, Tyr267, Leu294 as well as Cys283 and Met285 from another molecule within the tetrameric assembly. This might explain the tolerance of this enzyme toward organic solvents. The zinc ion occupies a position in the active site, which is too remote for direct interaction with the alcohol group. A mechanism is suggested whereby the introduction of NADP would trigger a displacement of the zinc ion to its catalytic site. Features important for the unusually high melting temperature of 98 degrees C are suggested by comparison to the crystal structure of a highly homologous mesophilic alcohol dehydrogenase from Clostridium beijerinckii (CBAD). The thermophilic enzyme has a more hydrophilic exterior, a more hydrophobic interior, a smaller surface area, more prolines, alanines, and fewer serines than CBAD. Furthermore, in the thermophilic enzyme the number of all types of intersubunit interactions in these tetrameric enzymes is increased: more salt bridges, hydrogen bonds, and hydrophobic interactions. All these effects combined can account for the higher melting temperature of the thermophilic enzyme.

Crystal structures of thermostable xylose isomerases from Thermus caldophilus and Thermus thermophilus: possible structural determinants of thermostability

The crystal structures of highly thermostable xylose isomerases from Thermus thermophilus (TthXI) and Thermus caldophilus (TcaXI), both with the optimum reaction temperature of 90 degrees C, have been determined by X-ray crystallography. The model of TcaXI has been refined to an R-factor of 17.8 % for data extending to 2.3 A and that of TthXI to 17.1 % for data extending to 2.2 A. The tetrameric arrangement of subunits characterized by the 222-symmetry and the tertiary fold of each subunit in both TcaXI and TthXI are basically the same as in other xylose isomerases. Each monomer is composed of two domains. Domain I (residues 1 to 321) folds into the (beta/alpha)8-barrel. Domain II (residues 322 to 387), lacking beta-strands, makes extensive contacts with domain I of an adjacent subunit. Each monomer of TcaXI contains ten beta-strands, 15 alpha-helices, and six 310-helices, while that of TthXI contains ten beta-strands, 16 alpha-helices, and five 310-helices. Although the electron density does not indicate the presence of bound metal ions in the present models of both TcaXI and TthXI, the active site residues show the conserved structural features. In order to understand the structural basis for thermostability of these enzymes, their structures have been compared with less thermostable XIs from Arthrobacter B3728 and Actinoplanes missouriensis (AXI and AmiXI), with the optimum reaction temperatures of 80 degrees C and 75 degrees C, respectively. Analyses of various factors that may affect protein thermostability indicate that the possible structural determinants of the enhanced thermostability of TcaXI/TthXI over AXI/AmiXI are (i) an increase in ion pairs and ion-pair networks, (ii) a decrease in the large inter-subunit cavities, (iii) a removal of potential deamidation/isoaspartate formation sites, and (iv) a shortened loop.

Analysis of the dynamics of rhizomucor miehei lipase at different temperatures

The dynamics of Rhizomucor miehei lipase has been studied by molecular dynamics simulations at temperatures ranging from 200-500K. Simulations carried out in periodic boundary conditions and using explicit water molecules were performed for 400 ps at each temperature. Our results indicate that conformational changes and internal motions in the protein are significantly influenced by the temperature increase. With increasing temperature, the number of internal hydrogen bonds decreases, while surface accessibility, radius of gyration and the number of residues in random coil conformation increase. In the temperature range studied, the motions can be described in a low dimensional subspace, whose dimensionality decreases with increasing temperature. Approximately 80% of the total motion is described by the first (i) 80 eigenvectors at T=200K, (ii) 30 eigenvectors at T=300K and (iii) 10 eigenvectors at T=400K. At high temperature, the alpha-helix covering the active site in the native Rhizomucor miehei lipase, the helix at which end the active site is located, and in particular, the loop (Gly35-Lys50) show extensive flexibility.

Molecular dynamics simulations of the hyperthermophilic protein sac7d from Sulfolobus acidocaldarius: contribution of salt bridges to thermostability

Hyperthermophilic proteins often possess an increased number of surface salt bridges compared with their mesophilic homologues. However, salt bridges are generally thought to be of minor importance in protein stability at room temperature. In an effort to understand why this may no longer be true at elevated temperatures, we performed molecular dynamics simulations of the hyperthermophilic protein Sac7d at 300 K, 360 K, and 550 K. The three trajectories are stable on the nanosecond timescale, as evidenced by the analysis of several time-resolved properties. The simulations at 300 K and (to a lesser extent) 360 K are also compatible with nuclear Overhauser effect-derived distances. Raising the temperature from 300 K to 360 K results in a less favourable protein-solvent interaction energy, and a more favourable intraprotein interaction energy. Both effects are almost exclusively electrostatic in nature and dominated by contributions due to charged side-chains. The reduced solvation is due to a loss of spatial and orientational structure of water around charged side-chains, which is a consequence of the increased thermal motion in the solvent. The favourable change in the intraprotein Coulombic interaction energy is essentially due to the tightening of salt bridges. Assuming that charged side-chains are on average more distant from one another in the unfolded state than in the folded state, it follows that salt bridges may contribute to protein stability at elevated temperatures because (i) the solvation free energy of charged side-chains is more adversely affected in the unfolded state than in the folded state by an increase in temperature, and (ii) due to the tightening of salt bridges, unfolding implies a larger unfavourable increase in the intraprotein Coulombic energy at higher temperature. Possible causes for the unexpected stability of the protein at 550 K are also discussed.

Refined crystal structure of a superoxide dismutase from the hyperthermophilic archaeon Sulfolobus acidocaldarius at 2.2 A resolution

The extremely thermostable superoxide dismutase from the hyperthermophilic archaeon Sulfolobus acidocaldarius was crystallized and the three-dimensional structure was determined by X-ray diffraction methods. The enzyme crystallized in the monoclinic spacegroup C2 with the cell dimensions a=168.1 A, b=91.3 A, c=85.7 A, beta=91.4 degrees. The diffraction limit of these crystals was 2.2 A. The crystals were very stable in the X-ray beam and measured diffraction data of a single crystal had a completeness of 99.5 % up to a resolution of 2.2 A. The crystal structure of S. acidocaldarius superoxide dismutase was solved by Patterson search methods using a dimer of Thermus thermophilus superoxide dismutase as a search model. The asymmetric unit accommodates three dimers. Two dimers form a tetramer by using only local symmetries; the third dimer forms a tetramer as well, however, by using the crystallographic 2-fold symmetry. The three-dimensional structure of the S. acidocaldarius dismutase has typical features of tetrameric dismutases. Secondary structure elements as well as residues important for the catalytic activity of the enzyme were found to be highly conserved. The model was refined at a resolution of 2.2 A and yielded a crystallographic R-value of 17.4 % (Rfree=22.3 %). A structural comparison of the two extremely stable tetrameric dismutases from S. acidocaldarius and Aquifex pyrophilus with the less stable enzyme from T. thermophilus and Mycoplasma tuberculosis revealed the structural determinants which are probably responsible for the high intrinsic stability of S. acidocaldarius dismutase. The most obvious factor which may give rise to the extraordinary thermal stability of S. acidocaldarius dismutase (melting temperature of about 125 degreesC) is the increase in intersubunit ion pairs and hydrogen bonds and, more importantly, the significant reduction of solvent-accessible hydrophobic surfaces, as well as an increase in the percentage of buried hydrophobic residues.

The stability of salt bridges at high temperatures: implications for hyperthermophilic proteins

Salt bridges have been proposed to play a crucial role in promoting hyperthermostability in proteins, yet they appear to make little contribution to protein stability at room temperature. The latter point has been rationalized previously on the basis that the association of two charged molecules to form a salt bridge incurs a substantial desolvation penalty, which is seldom completely compensated by favourable interactions within the salt bridge and with the rest of the protein. Here a continuum solvation model is used to investigate how this same argument applies at temperatures more appropriate to hyperthermophiles. The solvation model employed was previously parameterised to reproduce the hydration free energies of neutral and charged amino acid side-chains in the temperature range from 5-100 degreesC. A key result of the previous work was that the hydration free energies of charged side-chains are more adversely affected by increasing temperature than are the hydration free energies of hydrophobic side-chains of identical size and shape (isosteres). As is shown here, a direct consequence of the temperature dependence of the hydration free energies is that at high temperatures the desolvation penalty for formation of a salt bridge is markedly reduced in magnitude. As a result, the argument that relative to hydrophobic isosteres, salt bridges destabilise proteins, may no longer be true at high temperatures. We demonstrate this point first in the setting of a small model system, but then also show that the same argument is likely to carry over to real proteins. We compare three hyperthermophilic proteins with their mesophilic homologues and find that hydration effects preferentially stabilise the hyperthermophiles at high temperatures. When the hydration effects are incorporated into a model for the free energy of folding of the proteins, it is found that in each case, the hyperthermophile is predicted to remain stable to a temperature 20-25 deg.C higher than the corresponding mesophile. Higher thermal stability for the hyperthermophile is obtained even if the mesophile is more stable at room temperature. The results obtained therefore suggest one possible way in which the apparently destabilising effects of salt bridges at room temperature can be reconciled with their increased abundance in hyperthermophilic proteins.

Protein thermostability above 100 degreesC: a key role for ionic interactions

The discovery of hyperthermophilic microorganisms and the analysis of hyperthermostable enzymes has established the fact that multisubunit enzymes can survive for prolonged periods at temperatures above 100 degreesC. We have carried out homology-based modeling and direct structure comparison on the hexameric glutamate dehydrogenases from the hyperthermophiles Pyrococcus furiosus and Thermococcus litoralis whose optimal growth temperatures are 100 degreesC and 88 degreesC, respectively, to determine key stabilizing features. These enzymes, which are 87% homologous, differ 16-fold in thermal stability at 104 degreesC. We observed that an intersubunit ion-pair network was substantially reduced in the less stable enzyme from T. litoralis, and two residues were then altered to restore these interactions. The single mutations both had adverse effects on the thermostability of the protein. However, with both mutations in place, we observed a fourfold improvement of stability at 104 degreesC over the wild-type enzyme. The catalytic properties of the enzymes were unaffected by the mutations. These results suggest that extensive ion-pair networks may provide a general strategy for manipulating enzyme thermostability of multisubunit enzymes. However, this study emphasizes the importance of the exact local environment of a residue in determining its effects on stability.

Ion pairs involved in maintaining a thermostable structure of glutamate dehydrogenase from a hyperthermophilic archaeon

Intersubunit ion pairs are considered to be involved for maintaining a stable structure of the glutamate dehydrogenase (GDH) from hyperthermophiles. In order to demonstrate an effect of intersubunit ion pairs on the structural stability, two kinds of mutation (T138E, Thr at position 138 was replaced by Glu; E158Q, Glu at position 158 was replaced by Gln) which add and remove ion pairs, respectively, were introduced into Pk-gdhA gene encoding GDH from Pyrococcus kodakaraensis KOD1. Addition of one ion pair (Pk-GDHA-T138E) increased the optimum temperature and thermostability. In contrast, Pk-GDH-E158Q showed lower optimum temperature and less thermostability than wild type GDH. Structure analysis of GDHs was performed by circular dichroism (CD) and indicated that all recombinant enzymes (Pk-GDH, Pk-GDH-T138E, Pk-GDH-E158Q) possess different structures from that of natural GDH. Upon heat treatment (60 degrees C, 2 h), the structures of Pk-GDH and Pk-GDH-T138E were converted to another form close to the natural structure. However, no structural conversion by heat treatment was observed in Pk-GDH-E158Q. These results indicate that intersubunit ion pairs play an important role in forming thermostable structure of Pk-GDH.

Proteins from hyperthermophiles: stability and enzymatic catalysis close to the boiling point of water

It has become clear since about a decade ago, that the biosphere contains a variety of microorganisms that can live and grow in extreme environments. Hyperthermophilic microorganisms, present among Archaea and Bacteria, proliferate at temperatures of around 80-100 degrees C. The majority of the genera known to date are of marine origin, however, some of them have been found in continental hot springs and solfataric fields. Metabolic processes and specific biological functions of these organisms are mediated by enzymes and proteins that function optimally under these extreme conditions. We are now only starting to understand the structural, thermodynamic and kinetic basis for function and stability under conditions of high temperature, salt and extremes of pH. Insights gained from the study of such macromolecules help to extend our understanding of protein biochemistry and -biophysics and are becoming increasingly important for the investigation of fundamental problems in structure biology such as protein stability and protein folding. Extreme conditions in the biosphere require either the adaptation of the amino acid sequence of a protein by mutations, the optimization of weak interactions within the protein and at the protein-solvent boundary, the influence of extrinsic factors such as metabolites, cofactors, compatible solutes. Furthermore folding catalysts, known as chaperones, that assist the folding of proteins may be involved or increased protein protein synthesis in order to compensate for destruction by extreme conditions. The comparison of structure and stability of homologous proteins from mesophiles and hyperthermophiles has revealed important determinants of thermal stability of proteins. Rather than being the consequence of one dominant type of interactions or of a general stabilization strategy, it appears that the adaptation to high temperatures reflects a number of subtle interactions, often characteristic for each protein species, that minimize the surface energy and the hydration of apolar surface groups while burying hydrophobic residues and maximizing packing of the core as well as the energy due to charge-charge interactions and hydrogen bonds. In this article, mechanisms of intrinsic stabilization of proteins are reviewed. These mechanisms are found on different levels of structural organization. Among the extrinsic stabilization factors, emphasis is put on archaea chaperonins and their still strongly debated function. It will be shown, that optimization of weak protein-protein and protein-solvent interactions plays a key role in gaining thermostability. The difficulties in correlating suitable optimization criteria with real thermodynamic stability measures are due to experimental difficulties in measuring stabilization energies in large proteins or protein oligomers and will be discussed. Thus small single domain proteins or isolated domains of larger proteins may serve as model systems for large or multidomain proteins which due to the complexity of their thermal unfolding transitions cannot be analyzed by equilibrium thermodynamics. The analysis of the energetics of the thermal unfolding of a small, hyperthermostable DNA binding protein from Sulfolobus has revealed that a high melting temperature is not synonymous with a larger maximum thermodynamic stability. Finally, it is now well documented, that many thermophilic and hyperthermophilic proteins show a statistically increased number of salt bridges and salt bridge networks. However their contribution to thermodynamic and functional stability is still obscure.

Engineering activity and stability of Thermotoga maritima glutamate dehydrogenase. I. Introduction of a six-residue ion-pair network in the hinge region

Comparison of the recently determined three-dimensional structures of several glutamate dehydrogenases allowed for the identification of a five-residue ion-pair network in the hinge region of Pyrococcus furiosus glutamate dehydrogenase (melting temperature 113 degrees C), that is not present in the homologous glutamate dehydrogenase from Thermotoga maritima (melting temperature 93 degrees C). In order to study the role of this ion-pair network, we introduced it into the T. maritima enzyme using a site-directed mutagenesis approach. The resulting T. maritima glutamate dehydrogenases N97D, G376 K and N97D/G376 K as well as the wild-type enzyme were overproduced in Escherichia coli and subsequently purified. Elucidation of the three-dimensional structure of the double mutant N97D/G376 K at 3.0 A, showed that the designed ion-pair interactions were indeed formed. Moreover, because of interactions with an additional charged residue, a six-residue network is present in this double mutant. Melting temperatures of the mutant enzymes N97D, G376 K and N97D/G376 K, as determined by differential scanning calorimetry, did not differ significantly from that of the wild-type enzyme. Identical transition midpoints in guanidinium chloride-induced denaturation experiments were found for the wild-type and all mutant enzymes. Thermal inactivation at 85 degrees C occured more than twofold faster for all mutant enzymes than for the wild-type glutamate dehydrogenase. At temperatures of 65 degrees C and higher, the wild-type and the three mutant enzymes showed identical specific activities. However, at 58 degrees C the specific activity of N97D/G376 K and G376 K was found to be significantly higher than that of the wild-type and N97D enzymes. These results suggest that the engineered ion-pair interactions in the hinge region do not affect the stability towards temperature or guanidinium chloride-induced denaturation but rather affect the specific activity of the enzyme and the temperature at which it functions optimally.

Linkage of protonation and anion binding to the folding of Sac7d

The temperature, pH, and salt dependence of the folding of recombinant Sac7d from the hyperthermophile Sulfolobus acidocaldarius is mapped using multi-dimensional differential scanning calorimetry (DSC) and folding progress surfaces followed by circular dichroism. Linkage relations are derived to explain the observed dependencies, and it is shown that the data can be explained by the linkage of at least two protonation reactions and two anion binding sites to a two-state unfolding process. Circular dichroism spectra indicate that a native-like fold is stabilized at acid pH by anion binding. An apparent binding isotherm surface (folding progress versus pH and salt) is used to obtain intrinsic chloride binding constants as a function of pH for both sites. A saddle is predicted in the folding progress surface (progress versus temperature and pH) at low salt with a minimum near pH 2 and 20 degrees C with approximately 25% of the protein folded. The position of the saddle is sensitive to the intrinsic delta C degrees of unfolding and provides a third measure of delta C degrees independent of that obtained by a Kirchoff plot of DSC data and chemical denaturation. The observed enthalpy of unfolding approaches zero near the saddle making the unfolding largely invisible to DSC under these conditions. The linkage analysis demonstrates that the delta C degrees for unfolding obtained from a Kirchoff plot of DSC data should be distinguished from the intrinsic delta C degrees of unfolding. It is shown that the discrepancy between the free energy of unfolding for Sac7d obtained by DSC and that obtained by chemical denaturation may be explained by the linkage of protonation and anion binding to protein folding. The linkage analysis demonstrates the limitations of using the delta Hcal/ delta Hvh ratio an indication of two-state unfolding.

Dynamics and unfolding pathways of a hyperthermophilic and a mesophilic rubredoxin

Molecular dynamics simulations in solution are performed for a rubredoxin from the hyperthermophilic archaeon Pyrococcus furiosus (RdPf) and one from the mesophilic organism Desulfovibrio vulgaris (RdDv). The two proteins are simulated at four temperatures: 300 K, 373 K, 473 K (two sets), and 500 K; the various simulations extended from 200 ps to 1,020 ps. At room temperature, the two proteins are stable, remain close to the crystal structure, and exhibit similar dynamic behavior; the RMS residue fluctuations are slightly smaller in the hyperthermophilic protein. An analysis of the average energy contributions in the two proteins is made; the results suggest that the intraprotein energy stabilizes RdPf relative to RdDv. At 373 K, the mesophilic protein unfolds rapidly (it begins to unfold at 300 ps), whereas the hyperthermophilic does not unfold over the simulation of 600 ps. This is in accord with the expected stability of the two proteins. At 473 K, where both proteins are expected to be unstable, unfolding behavior is observed within 200 ps and the mesophilic protein unfolds faster than the hyperthermophilic one. At 500 K, both proteins unfold; the hyperthermophilic protein does so faster than the mesophilic protein. The unfolding behavior for the two proteins is found to be very similar. Although the exact order of events differs from one trajectory to another, both proteins unfold first by opening of the loop region to expose the hydrophobic core. This is followed by unzipping of the beta-sheet. The results obtained in the simulation are discussed in terms of the factors involved in flexibility and thermostability.

Optimization of the electrostatic interactions between ionized groups and peptide dipoles in proteins

The three-dimensional optimization of the electrostatic interactions between the charged amino acid residues and the peptide partial charges was studied by Monte-Carlo simulations on a set of 127 nonhomologous protein structures with known atomic coordinates. It was shown that this type of interaction is very well optimized for all proteins in the data set, which suggests that they are a valuable driving force, at least for the native side-chain conformations. Similar to the optimization of the charge-charge interactions (Spassov VZ, Karshikoff AD, Ladenstein R, 1995, Protein Sci 4:1516-1527), the optimization effect was found more pronounced for enzymes than for proteins without enzymatic function. The asymmetry in the interactions of acidic and basic groups with the peptide dipoles was analyzed and a hypothesis was proposed that the properties of peptide dipoles are a factor contributing to the natural selection of the basic amino acids in the chemical composition of proteins.

Crystal structure of glutamate dehydrogenase from the hyperthermophilic eubacterium Thermotoga maritima at 3.0 A resolution

The extremely thermostable glutamate dehydrogenase from the hyperthermophilic bacterium Thermotoga maritima has been crystallized and the three-dimensional structure has been determined by X-ray diffraction methods. Crystals up to a maximum size of 1.2 mm have been grown in 3% polyethylene glycol, 120 mM ammonium acetate and 50 mM bis-tris propane (pH 6.5). The enzyme crystallized in the trigonal space group P3(1)21 with the cell dimensions a = b = 147.3 A, c = 273.6 A. The diffraction limit of these crystals is 3.0 A. Measured diffraction data have a completeness of 94% up to a resolution of 3.0 A and contain 75% of all possible data in the last resolution shell between 3.1 and 3.0 A. The crystal structure of T. maritima glutamate dehydrogenase has been solved by Patterson search methods using the hexameric Pyrococcus furiosus glutamate dehydrogenase as a search model. The crystallographic refinement has been carried out to a maximum resolution of 3.1 A and an crystallographic R-value of 22.5% (Rfree = 29.5%). The three-dimensional structure of the T. maritima enzyme shows typical features of hexameric glutamate dehydrogenases: six subunits are arranged in 32 symmetry. Each subunit consists of two domains connected by a flexible hinge region. Secondary structure elements as well as residues important for the catalytic activity of the enzyme are highly conserved. A structural comparison of the two glutamate dehydrogenases from the hyperthermophiles T. maritima and P. furiosus with the enzyme from the mesophilic bacterium Clostridium symbiosum has revealed that common as well as distinct mechanisms contribute to the thermal stability of these enzymes. The number of intrasubunit ion pairs is increased and the volume of intrasubunit cavities decreased in both thermostable enzymes, whereas striking differences have been observed in the subunit interfaces. In P. furiosus glutamate dehydrogenase the subunit interactions are dominated by ionic interactions realized by large saltbridge networks. However, in T. maritima glutamate dehydrogenase the number of intersubunit ion pairs is reduced and the hydrophobic interactions are increased.

Bacillus thermoamyloliquefaciens KP1071 alpha-glucosidase II is a thermostable M(r) 540,000 homohexameric alpha-glucosidase with both exo-alpha-1,4-glucosidase and oligo-1,6-glucosidase activities

alpha-Glucosidase II of the facultative thermophile Bacillus thermoamyloliquefaciens KP1071 (FERM-P8477; growth over 30-66 degrees C) was purified to a homogeneous state. Its M(r) was estimated as 90000 by SDS/PAGE. However, the enzyme behaved as an active Mr 540000 protein on gel filtration with each of two gels of different matrices as well as on gel electrophoresis under native conditions. The enzyme was not glycosylated. Its isoelectric point was estimated as 5.7. The N-terminal sequence of 20 residues was determined asAla1-Ile-Gln-Pro-Glu-Gln-Asp-Asp-Lys-Thr-Gln-Glu-Asp-Gly- Tyr-Ile-Asp-Ile-Gly-Asn20. The sequence did not resemble those of procaryotic and eucaryotic proteins hitherto reported including the monomeric exo-alpha-1,4-glucosidase and the monomeric oligo-1,6-glucosidase from the same microorganism. The alpha-glucosidase II had no antigenic group shared with the latter two enzymes. Analysis of substrate specificity showed that the alpha-glucosidase II has dual activity towards oligo-1,6-glucosidases and exo-alpha-1,4-glucosidases, but its preference is for non-reducing terminal alpha-1,4 glucosidic bonds in substrates. Kinetic studies proved that both activities are attributed to the same catalytic site. The enzyme was most active at 81 degrees C and pH 7.0. Its half-life at pH 6.8 was 10 min at 81 degrees C, and 5 h at 55 degrees C in 6.4 M urea, 26% ethanol or 2.5% SDS. We suggest that the alpha-glucosidase II is a thermostable, homohexameric enzyme of origin distinct from the exo-alpha-1,4-glucosidase and the oligo-1,6-glucosidase present in the same strain.