Insights into thermal stability from a comparison of the glutamate dehydrogenases from Pyrococcus furiosus and Thermococcus litoralis

Parallel molecular mechanisms for enzyme temperature adaptation

The mechanisms that underly the adaptation of enzyme activities and stabilities to temperature are fundamental to our understanding of molecular evolution and how enzymes work. Here, we investigate the molecular and evolutionary mechanisms of enzyme temperature adaption, combining deep mechanistic studies with comprehensive sequence analyses of thousands of enzymes. We show that temperature adaptation in ketosteroid isomerase (KSI) arises primarily from one residue change with limited, local epistasis, and we establish the underlying physical mechanisms. This residue change occurs in diverse KSI backgrounds, suggesting parallel adaptation to temperature. We identify residues associated with organismal growth temperature across 1005 diverse bacterial enzyme families, suggesting widespread parallel adaptation to temperature. We assess the residue properties, molecular interactions, and interaction networks that appear to underly temperature adaptation.

Improvement of enzymatic properties of Rhizopus oryzae α-amylase by site-saturation mutagenesis of histidine 286

Optimal pH and ideal functioning temperature for fungal α-amylase can greatly contribute to improving enzyme efficiency in maltose-forming ability. This work aimed to improve the enzymatic properties of Rhizopus oryzae α-amylase by site-saturation mutagenesis of histidine 286. The biochemical properties of selected mutant enzymes were modified to increase their enzymatic efficiencies compared to their wild-type counterparts. For instance, the optimum temperature of mutants H286 L, H286I, H286S and H286 T was increased from 50 °C to 55 °C, while a similar increase was observed for H286 P from 50 °C to 60 °C. The optimum pH of mutants H286 L, H286I and H286D shifted from 5.5 to 5.0, and the optimum pH of mutant H286E shifted from 5.5 to 4.5. The results obtained showed that the mutant H286I showed a 1.5-fold increase in half-life at 55 °C and the mutant H286E showed a 6.43-fold increase in half-life at a pH of 4.5. Furthermore, the ability to form maltose from soluble starch for mutants H286 L and H286 M was significantly improved under the optimum conditions determined in the study. The catalytic mechanism responsible for improved maltose-forming ability was confirmed through molecular docking simulations with maltotriose among wild-type and mutant enzymes. The mutants with improved enzymatic properties that were attained in this work may help in future computer-aided directed evolution of fungal α-amylase.

Multifactorial level of extremostability of proteins: can they be exploited for protein engineering?

Research on extremostable proteins has seen immense growth in the past decade owing to their industrial importance. Basic research of attributes related to extreme-stability requires further exploration. Modern mechanistic approaches to engineer such proteins in vitro will have more impact in industrial biotechnology economy. Developing a priori knowledge about the mechanism behind extreme-stability will nurture better understanding of pathways leading to protein molecular evolution and folding. This review is a vivid compilation about all classes of extremostable proteins and the attributes that lead to myriad of adaptations divulged after an extensive study of 6495 articles belonging to extremostable proteins. Along with detailing on the rationale behind extreme-stability of proteins, emphasis has been put on modern approaches that have been utilized to render proteins extremostable by protein engineering. It was understood that each protein shows different approaches to extreme-stability governed by minute differences in their biophysical properties and the milieu in which they exist. Any general rule has not yet been drawn regarding adaptive mechanisms in extreme environments. This review was further instrumental to understand the drawback of the available 14 stabilizing mutation prediction algorithms. Thus, this review lays the foundation to further explore the biophysical pleiotropy of extreme-stable proteins to deduce a global prediction model for predicting the effect of mutations on protein stability.

Evaluation and directed evolution for thermostability improvement of a GH 13 thermostable α-glucosidase from Thermus thermophilus TC11

Background

Thermal stable α-glucosidases with transglycosylation activity could be applied to the industrial production of oligosaccharides as well as conjugation of sugars to biologically useful materials. Therefore, α-glucosidases isolated from thermophiles have gained attention over the past decade. In this study, the characterization of a highly thermostable α-glucosidase and its thermostability improved mutant from newly isolated strain Thermus thermophilus TC11 were investigated.

Results

The recombinant α-glucosidase (TtAG) from Thermus thermophilus TC11 was expressed in Escherichia coli BL21 (DE3) and purified. The purified enzyme had a molecular mass of 184 kDa and consisted of 59-kDa subunits; it showed hydrolytic activity for pNP-α-d-glucopyranoside (pNPG), sucrose, trehalose, panose, and isomaltooligosaccharides and very low activity for maltose. The highest specific activity of 288.96 U/mg was observed for pNPG at 90 °C and pH 5.0; Pb²⁺ provided a 20 % activity increase. TtAG was stable at 70 °C for more than 7 h and had a half-life of 195 min at 80 °C and 130 min at 90 °C. Transglycosylation activity was also observed with sucrose and trehalose as substrates. TtAG showed differences on substrate specificity, transglycosylation, multimerization, effects of metal ions and optimal pH from other reported Thermus α-glucosidases. One single-substitution TtAG mutant Q10Y with improved thermostability was also obtained from random mutagenesis library. The site-saturation mutagenesis and structural modelling analysis indicated that Q10Y substitution stabilized TtAG structure via additional hydrogen bonding and hydrophobic interactions.

Conclusion

Our findings indicate that TtAG is a highly thermostable and more acidic α-glucosidase distinct from other reported Thermus α-glucosidases. And this work also provides new insights into the catalytic and thermal tolerance mechanisms of α-glucosidases, which may guide molecular engineering of α-glucosidase and other thermostable enzymes for industrial application.

Electronic supplementary material

The online version of this article (doi:10.1186/s12896-015-0197-x) contains supplementary material, which is available to authorized users.

Rigidity versus flexibility: the dilemma of understanding protein thermal stability

The role of fluctuations in protein thermostability has recently received considerable attention. In the current literature a dualistic picture can be found: thermostability seems to be associated with enhanced rigidity of the protein scaffold in parallel with the reduction of flexible parts of the structure. In contradiction to such arguments it has been shown by experimental studies and computer simulation that thermal tolerance of a protein is not necessarily correlated with the suppression of internal fluctuations and mobility. Both concepts, rigidity and flexibility, are derived from mechanical engineering and represent temporally insensitive features describing static properties, neglecting that relative motion at certain time scales is possible in structurally stable regions of a protein. This suggests that a strict separation of rigid and flexible parts of a protein molecule does not describe the reality correctly. In this work the concepts of mobility/flexibility versus rigidity will be critically reconsidered by taking into account molecular dynamics calculations of heat capacity and conformational entropy, salt bridge networks, electrostatic interactions in folded and unfolded states, and the emerging picture of protein thermostability in view of recently developed network theories. Last, but not least, the influence of high temperature on the active site and activity of enzymes will be considered.

Prediction of the determinants of thermal stability by linear discriminant analysis: the case of the glutamate dehydrogenase protein family

Little is known about the determinants of thermal stability in individual protein families. Most of the knowledge on thermostability comes, in fact, from comparative analyses between large, and heterogeneous, sets of thermo- and mesophilic proteins. Here, we present a multivariate statistical approach aimed to detect signature sequences for thermostability in a single protein family. It was applied to the glutamate dehydrogenase (GDH) family, which is a good model for investigating this peculiar process. The structure of GDH consists of six subunits, each of them organized into two domains. Formation of ion-pair networks on the surface of the protein subunits, or increase in the inter-subunit hydrophobic interactions, have been suggested as important factors for explaining stability at high temperatures. However, identification of the amino acid changes that are involved in this process still remains elusive. Our approach consisted of a linear discriminant analysis on a set of GDH sequences from Archaea and Bacteria (33 thermo- and 36 mesophilic GDHs). It led to detection of 3 amino acid clusters as the putative determinants of thermal stability. They were localized at the subunit interface or in close proximity to the binding site of the NAD(P)(+) coenzyme. Analysis within the clusters led to prediction of 8 critical amino acid sites. This approach could have a wide utility, in the ligth of the notion that each protein family seems to adopt its own strategy for achieving thermostability.

Evaluation of sub-cellular distribution of glutamate dehydrogenase (GDH) in Drosophila melanogaster larvae

Glutamate dehydrogenase (GDH) enzyme was conventionally known as a mitochondrial marker. However, subsequently it was reported to be present in the nuclei as well. So far, the nuclear distribution of GDH has been reported in a number of organisms including yeast, rat, cow, chicken. However, the sub-cellular distribution of GDH, illustrated by in situ methods still remains elusive. Here, by assaying the GDH activity and by immuno-blotting using anti-GDH antibody in the fractionated nuclear and cytoplasmic fractions of Drosophila larvae, we demonstrate the cytoplasmic distribution of GDH. This observation was further supported by in situ immunostaining of salivary gland, Malpighian tubules and eye imaginal discs of Drosophila larvae. Collectively, our results demonstrate that in Drosophila larvae, GDH is not found in the nucleus, but is localized exclusively in the cytoplasm.

A role for A-to-I RNA editing in temperature adaptation

A-to-I RNA editing can recode mRNAs, giving organisms the option to express diverse, functionally distinct protein isoforms. Here, we propose that RNA editing is inherently geared for temperature adaptation because it tends to recode to smaller, less stabilizing amino acids. Studies on how editing affects protein function support this idea.

Purification and characterization of a thermostable glutamate dehydrogenase from a thermophilic bacterium isolated from a sterilization drying oven

Glutamate dehydrogenase from axenic bacterial cultures of a new microorganism, called GWE1, isolated from the interior of a sterilization drying oven, was purified by anion-exchange and molecular-exclusion liquid chromatography. The apparent molecular mass of the native enzyme was 250.5 kDa and was shown to be an hexamer with similar subunits of molecular mass 40.5 kDa. For glutamate oxidation, the enzyme showed an optimal pH and temperature of 8.0 and 70 degrees C, respectively. In contrast to other glutamate dehydrogenases isolated from bacteria, the enzyme isolated in this study can use both NAD(+) and NADP(+) as electron acceptors, displaying more affinity for NADP(+) than for NAD(+). No activity was detected with NADH or NADPH, 2-oxoglutarate and ammonia. The enzyme was exceptionally thermostable, maintaining more than 70% of activity after incubating at 100(o)C for more than five hours suggesting being one of the most thermoestable enzymes reported in the family of dehydrogenases.

Enhancing the thermostability of alpha-glucosidase from Thermoanaerobacter tengcongensis MB4 by single proline substitution

Thermostability can be increased by introducing prolines at suitable sites in target proteins. In this study, we compared five thermostable alpha-glucosidases and the moderate thermostable alpha-glucosidase (TtGluA) from Thermoanaerobacter tengcongensis MB4. Based on the amino acid sequence alignment, four sites (Leu152, Asn208, Lys285, and Thr430) of TtGluA were chosen for proline substitution to improve its thermostability. Thermostability of mutants L152P, K285P, and T430P increased evidently, but no thermostability improvement was observed for N208P. Compared to the wild-type enzyme, T(50)(15) of T430P had a rise of 2 degrees C without distinct loss of activity. However, T(50)(15) values of L152P and K285P increased 2 degrees C and 10.5 degrees C, respectively, while retaining activity of only 26.6% and 24.9% of wild-type enzyme. The K(m) of L152P, K285P, T430P and wild-type enzyme was 1.61, 0.32, 1.64, and 1.08 mM, respectively. These indicate that the selected sites are not only important for the thermostability but also related to the substrate binding and catalytic activity of TtGluA. The CD spectra analysis of the improved mutants and wild-type enzyme showed no distinct changes in their secondary structures. Combining analysis of secondary structure prediction and 3D structure modeling, the proline substitution at the three sites stabilized TtGluA possibly by reducing the flexibility of loop and random coil or (and) increasing the hydrophobic effect at these strategic regions with no evident structure change.

Salt bridges in the hyperthermophilic protein Ssh10b are resilient to temperature increases

A double mutant cycle (DMC) approach was employed to estimate the effect of temperature on the contribution of two highly conserved salt bridges to protein stability in the hyperthermophilic protein Ssh10b. The coupling free energy were 2.4 +/- 0.4 kJ/mol at 298 K and 2.2 +/- 0.4 kJ/mol at 353 K for Glu-54/Arg-57, and 6.0 +/- 0.2 kJ/mol at 298 K and 5.9 +/- 0.6 kJ/mol at 353 K for Glu-36/Lys-68. The stability free energy of Ssh10b decrease greatly with increasing temperature, while the direct contribution of these two salt bridges to protein stability remain almost constant, providing evidence supporting the theoretical prediction that salt bridges are extremely resilient to temperature increases and thus are specially suited to improving protein stability at high temperatures. The reason for the difference in coupling free energy between salt bridges Glu-54/Arg-57 and Glu-36/Lys-68 is discussed. Comparing our results with published DMC data for the contribution of salt bridges to stability in other proteins, we found that the energy contribution of a salt bridge formed by two charged residues far apart in the primary sequence is higher than that of those formed between two very close ones. Implications of this finding are useful for engineering proteins with enhanced thermostability.

The backbone structure of the thermophilic Thermoanaerobacter tengcongensis ribose binding protein is essentially identical to its mesophilic E. coli homolog

Background

Comparison of experimentally determined mesophilic and thermophilic homologous protein structures is an important tool for understanding the mechanisms that contribute to thermal stability. Of particular interest are pairs of homologous structures that are structurally very similar, but differ significantly in thermal stability.

Results

We report the X-ray crystal structure of a Thermoanaerobacter tengcongensis ribose binding protein (tteRBP) determined to 1.9 Å resolution. We find that tteRBP is significantly more stable (^appT_m value ~102°C) than the mesophilic Escherichia coli ribose binding protein (ecRBP) (^appT_m value ~56°C). The tteRBP has essentially the identical backbone conformation (0.41 Å RMSD of 235/271 C_α positions and 0.65 Å RMSD of 270/271 C_α positions) as ecRBP. Classification of the amino acid substitutions as a function of structure therefore allows the identification of amino acids which potentially contribute to the observed thermal stability of tteRBP in the absence of large structural heterogeneities.

Conclusion

The near identity of backbone structures of this pair of proteins entails that the significant differences in their thermal stabilities are encoded exclusively by the identity of the amino acid side-chains. Furthermore, the degree of sequence divergence is strongly correlated with structure; with a high degree of conservation in the core progressing to increased diversity in the boundary and surface regions. Different factors that may possibly contribute to thermal stability appear to be differentially encoded in each of these regions of the protein. The tteRBP/ecRBP pair therefore offers an opportunity to dissect contributions to thermal stability by side-chains alone in the absence of large structural differences.

Purification and characterization of trypsin-like enzymes from North Pacific krill (Euphausia pacifica)

Two trypsin-like enzymes (TLEs) were purified from North Pacific krill (Euphausia pacifica) by ammonium sulfate precipitation, ion-exchange and gel-filtration chromatography. The purified enzymes were identified as trypsins by LC-ESI-MS/MS analysis. The relative molecular mass of TLE I and TLE II were 33 and 32.3 kDa, respectively, with isoelectric points of 4.5 and 4.3, respectively. The TLEs showed excellent thermal stable in the crude extract and the purified TLEs were active over a wide pH (6.0-11.0) and temperature (10-70 degrees C) range. Compared with trypsins from other organisms, the purified TLEs had physiological efficiencies of 1.6-6.7-fold. The difference in Arg, Ile and Asp content might explain why E. pacifica TLEs have good thermal stability and physiological efficiency.

Kinetically controlled refolding of a heat-denatured hyperthermostable protein

The thermal denaturation of endo-beta-1,3-glucanase from the hyperthermophilic microorganism Pyrococcus furiosus was studied by calorimetry. The calorimetric profile revealed two transitions at 109 and 144 degrees C, corresponding to protein denaturation and complete unfolding, respectively, as shown by circular dichroism and fluorescence spectroscopy data. Calorimetric studies also showed that the denatured state did not refold to the native state unless the cooling temperature rate was very slow. Furthermore, previously denatured protein samples gave well-resolved denaturation transition peaks and showed enzymatic activity after 3 and 9 months of storage, indicating slow refolding to the native conformation over time.

Amino-acid interactions in psychrophiles, mesophiles, thermophiles, and hyperthermophiles: insights from the quasi-chemical approximation

We investigate the mechanisms used by proteins to maintain thermostability throughout a wide range of temperatures. We use the quasi-chemical approximation to estimate interaction strengths for psychrophiles, mesophiles, thermophiles, and hyperthermophiles. Our results highlight the importance of core packing in thermophilic stability. Although we observed an increase in the number of charged residues, the contribution of salt bridges appears to be relatively modest by comparison. We observed results consistent with a gradual loosening of structure in psychrophiles, including a weakening of almost all types of interactions.

Differences in amino acids composition and coupling patterns between mesophilic and thermophilic proteins

Thermophilic proteins show substantially higher intrinsic thermal stability than their mesophilic counterparts. Amino acid composition is believed to alter the intrinsic stability of proteins. Several investigations and mutagenesis experiment have been carried out to understand the amino acid composition for the thermostability of proteins. This review presents some generalized features of amino acid composition found in thermophilic proteins, including an increase in residue hydrophobicity, a decrease in uncharged polar residues, an increase in charged residues, an increase in aromatic residues, certain amino acid coupling patterns and amino acid preferences for thermophilic proteins. The differences of amino acids composition between thermophilic and mesophilic proteins are related to some properties of amino acids. These features provide guidelines for engineering mesophilic protein to thermophilic protein.

Hyperthermophilic enzymes − stability, activity and implementation strategies for high temperature applications

Current theories agree that there appears to be no unique feature responsible for the remarkable heat stability properties of hyperthermostable proteins. A concerted action of structural, dynamic and other physicochemical attributes are utilized to ensure the delicate balance between stability and functionality of proteins at high temperatures. We have thoroughly screened the literature for hyperthermostable enzymes with optimal temperatures exceeding 100 degrees C that can potentially be employed in multiple biotechnological and industrial applications and to substitute traditionally used, high-cost engineered mesophilic/thermophilic enzymes that operate at lower temperatures. Furthermore, we discuss general methods of enzyme immobilization and suggest specific strategies to improve thermal stability, activity and durability of hyperthermophilic enzymes.

Molecular characterization of cold adaptation based on ortholog protein sequences from Vibrionaceae species

A set of 298 protein families from psychrophilic Vibrio salmonicida was compiled to identify genotypic characteristics that discern it from orthologous sequences from the mesophilic Vibrio/Photobacterium branch of the gamma-Proteobacteria (Vibrionaceae family). In our comparative exploration we employed alignment based bioinformatical and statistical methods. Interesting information was found in the substitution matrices, and the pattern of asymmetries in the amino acid substitution process. Together with the compositional difference, they identified the amino acids Ile, Asn, Ala and Gln as those having the most psycrophilic involvement. Ile and Asn are enhanced whereas Gln and Ala are suppressed. The inflexible Pro residue is also suppressed in loop regions, as expected in a flexible structure. The dataset were also classified and analysed according to the predicted subcellular location, and we made an additional study of 183 intracellular and 65 membrane proteins. Our results revealed that the psychrophilic proteins have similar hydrophobic and charge contributions in the core of the protein as mesophilic proteins, while the solvent-exposed surface area is significantly more hydrophobic. In addition, the psychrophilic intracellular (but not the membrane) proteins are significantly more negatively charged at the surface. Our analysis supports the hypothesis of preference for more flexible amino acids at the molecular surface. Life in cold climate seems to be obtained through many minor structural modifications rather than certain amino acids substitutions.

Analysis of the thermostability determinants of hyperthermophilic esterase EstE1 based on its predicted three-dimensional structure

The three-dimensional (3D) structure of the hyperthermophilic esterase EstE1 was constructed by homology modeling using Archaeoglobus fulgidus esterase as a reference, and the thermostability-structure relationship was analyzed. Our results verified the predicted 3D structure of EstE1 and identified the ion pair networks and hydrophobic interactions that are critical determinants for the thermostability of EstE1.

Stabilization mechanism of the tryptophan synthase alpha-subunit from Thermus thermophilus HB8: X-ray crystallographic analysis and calorimetry

In order to elucidate the thermo-stabilization mechanism of the tryptophan synthase alpha-subunit from the extreme thermophile Thermus thermophilus HB8 (Tt-alpha-subunit), its crystal structure was determined and its stability was examined using DSC. The results were compared to those of other orthologs from mesophilic and hyperthermophilic organisms. The denaturation temperature of the Tt-alpha-subunit was higher than that of the alpha-subunit from S. typhimurium (St-alpha-subunit) but lower than that of the alpha-subunit from P. furiosus (Pf-alpha-subunit). Specific denaturation enthalpy and specific denaturation heat capacity values of the Tt-alpha-subunit were the lowest among the three proteins, suggesting that entropy effects are responsible for the stabilization of the Tt-alpha-subunit. Based on a structural comparison with the St-alpha-subunit, two deletions in loop regions, an increase in the number of ion pairs and a decrease in cavity volume seem to be responsible for the stabilization of the Tt-alpha-subunit. The results of structural comparison suggest that the native structure of the Tt-alpha-subunit is better adapted to an ideally stable structure than that of the St-alpha-subunit, but worse than that of the Pf-alpha-subunit. The results of calorimetry suggest that the residual structure of the Tt-alpha-subunit in the denatured state contributes to the stabilization.

Effective factors in thermostability of thermophilic proteins

Thermostability of proteins in general and especially thermophilic proteins has been subject of a wide variety of studies based on theoretical and experimental investigation. Thermostability seems to be a property obtained through many minor structural modifications rather than certain amino acids substitution. In comparison with its mesophile homologue in a thermostable protein, usually a number of amino acids are exchanged. A wide variety of theoretical studies are based on comparative investigation of thermophilic proteins characteristics with their mesophilic counterparts in order to reveal their sequences, structural differences and consequently, to relate these observed differences to the thermostability properties. In this work we have compared a dataset of thermophilic proteins with their mesophilic homologues and furthermore, a mesophilic proteins dataset was also compared with its mesophilic homologue. This strategy enabled us first, to eliminate noise or background differences from signals and moreover, the important factors which were related to the thermostability were recognized too. Our results reveal that thermophilic and mesophilic proteins have both similar polar and nonpolar contribution to the surface area and compactness. On the other hand, salt bridges and main chain hydrogen bonds show an increase in the majority of thermophilic proteins in comparison to their mesophilic homologues. In addition, in thermophilic proteins hydrophobic residues are significantly more frequent, while polar residues are less. These findings indicate that thermostable proteins through evolution adopt several different strategies to withstand high temperature environments.

Temperature-dependent structural and functional features of a hyperthermostable enzyme using elastic neutron scattering

The dynamic behavior of an endoglucanase from the hyperthermophilic microorganism Pyrococcus furiosus was investigated using elastic neutron scattering. The temperature dependence of the atomic motions was correlated with conformational and functional characteristics of the enzyme. The onset of biological function at temperatures higher than approximately 25 degrees C (the hyperthermostable enzyme is essentially inactive at room temperature) was associated with a dynamical transition in the anharmonic motions domain. This transition from the nonactive to the enzymatically active conformation involved structurally similar conformational substates in the energy landscape. From the mean-square displacement of the protein atoms, the molecular flexibility and the effective force constants were calculated at different temperature zones. The results showed that the activity increases at higher temperatures where the intramolecular bonds are weakened and the overall rigidity of the protein is decreased. Further temperature increase resulted in significantly increased atomic fluctuations featuring heat denaturation of the protein.

Characterization of cyclodextrin glycosyltransferase of the same gene expressed from Bacillus macerans, Bacillus subtilis, and Escherichia coli

The plasmid pHG contains a cyclodextrin glycosyltransferase (CGTase) gene (cgt) derived from Bacillus macerans. Two transformants, Bacillus subtilis (pHG) and Escherichia coli (pHG), were found to produce CGTases with the same primary structure as the enzyme from B. macerans. However, the beta-cyclodextrin coupling activity of the CGTase from E. coli (pHG) was 14-fold higher than that of the enzymes from the other strains. By contrast, no differences in alpha-cyclodextrin coupling activities were observed among these CGTases. CGTase from E. coli (pHG) was found to be less thermostable than the other CGTases. When the CGTase produced by B. subtilis was treated with increasing urea concentrations (10-1000 mM) to promote increasing degrees of protein unfolding, a bell-shaped beta-cyclodextrin coupling activity profile was obtained. Subtle differences in the conformation of the CGTase produced by E. coli are therefore proposed to be responsible for the markedly increased beta-cyclodextrin coupling activity of this enzyme.

The crystal structure of Plasmodium falciparum glutamate dehydrogenase, a putative target for novel antimalarial drugs

Plasmodium falciparum is the main causative agent of tropical malaria, the most severe parasitic disease in the world. Growing resistance of Plasmodia towards available drugs is an increasing problem in countries where malaria is endemic. As Plasmodia are sensitive to oxidative stress, augmenting this in the parasite represents a promising principle for the development of novel antimalarial drugs. The NADP-dependent glutamate dehydrogenase (GDH) of P.falciparum is largely responsible for the production of NADPH in the parasite, which in turn serves as electron source for the antioxidative enzymes glutathione reductase and thioredoxin reductase. As GDH does not occur in the host erythrocyte, GDH is a particularly attractive target for drug therapy. The three-dimensional structure of P.falciparum GDH in the unligated state has been determined by X-ray crystallography to a resolution of 2.7A. Compared to the mammalian enzymes, two amino acid residues are exchanged in the putative active site of the parasite GDH. The most obvious differences between parasite and human GDH are the subunit interfaces of the hexameric proteins. In the parasite protein, several salt-bridges mediate contacts between the subunits whereas in the human enzyme these interactions are mainly of hydrophobic nature. Furthermore, P.falciparum GDH possesses a unique N-terminal extension that does not occur in any other GDH sequence so far studied. These findings might be exploited for the design of peptidomimetics capable of disrupting the oligomeric organisation of the parasite enzyme.

Packing-based difference of structural features between thermophilic and mesophilic proteins

Twenty pairs of thermophilic and mesophilic proteins were compared in terms of residue packing distribution to obtain structural features related to protein thermostability. Based on residue packing concept, structural features of residues such as residue packing distribution, inner/outer position, secondary structure and water solvation were investigated. The statistical tests revealed that higher frequency in well-packed state of residues, lower frequency in exposed state and higher frequency in well-packed state of inner positioned residues, and higher frequency in well-packed state of 3/10 helix residues could be general structural features thermophilic proteins have.

Hyperthermal stability of neuroglobin and cytoglobin

Neuroglobin (Ngb) and cytoglobin (Cygb), recent additions to the globin family, display a hexa-coordinated (bis-histidyl) heme in the absence of external ligands. Although these proteins have the classical globin fold they reveal a very high thermal stability with a melting temperature (Tm) of 100 degrees C for Ngb and 95 degrees C for Cygb. Moreover, flash photolysis experiments at high temperatures reveal that Ngb remains functional at 90 degrees C. Human Ngb may have a disulfide bond in the CD loop region; reduction of the disulfide bond increases the affinity of the iron atom for the distal (E7) histidine, and leads to a 3 degrees C increase in the T(m) for ferrous Ngb. A similar Tm is found for a mutant of human Ngb without cysteines. Apparently, the disulfide bond is not involved directly in protein stability, but may influence the stability indirectly because it modifies the affinity of the distal histidine. Mutation of the distal histidine leads to lower thermal stability, similar to that for other globins. Only globins with a high affinity of the distal histidine show the very high thermal stability, indicating that stable hexa-coordination is necessary for the enhanced thermal stability; the CD loop which contains the cysteines appears as a critical region in the neuroglobin thermal stability, because it may influence the affinity of the distal histidine.

Structure and Function of a Regulated Archaeal Triosephosphate Isomerase Adapted to High Temperature

Triosephophate isomerase (TIM) is a dimeric enzyme in eucarya, bacteria and mesophilic archaea. In hyperthermophilic archaea, however, TIM exists as a tetramer composed of monomers that are about 10% shorter than other eucaryal and bacterial TIM monomers. We report here the crystal structure of TIM from Thermoproteus tenax, a hyperthermophilic archaeon that has an optimum growth temperature of 86 degrees C. The structure was determined from both a hexagonal and an orthorhombic crystal form to resolutions of 2.5A and 2.3A, and refined to R-factors of 19.7% and 21.5%, respectively. In both crystal forms, T.tenax TIM exists as a tetramer of the familiar (betaalpha)(8)-barrel. In solution, however, and unlike other hyperthermophilic TIMs, the T.tenax enzyme exhibits an equilibrium between inactive dimers and active tetramers, which is shifted to the tetramer state through a specific interaction with glycerol-1-phosphate dehydrogenase of T.tenax. This observation is interpreted in physiological terms as a need to reduce the build-up of thermolabile metabolic intermediates that would be susceptible to destruction by heat. A detailed structural comparison with TIMs from organisms with growth optima ranging from 15 degrees C to 100 degrees C emphasizes the importance in hyperthermophilic proteins of the specific location of ionic interactions for thermal stability rather than their numbers, and shows a clear correlation between the reduction of heat-labile, surface-exposed Asn and Gln residues with thermoadaptation. The comparison confirms the increase in charged surface-exposed residues at the expense of polar residues.

Thermodynamics of core hydrophobicity and packing in the hyperthermophile proteins Sac7d and Sso7d

The influence of core hydrophobicity and packing on the structure and stability of the hyperthermophile proteins Sac7d and Sso7d have been studied by calorimetry, circular dichroism, and NMR. Valine 30 is positioned in Sac7d to allow a cavity-filling Val --> Ile substitution which occurs naturally in the homologous more thermostable Sso7d. The cavity-filling mutation in Sac7d has been characterized and compared to the reciprocal Ile --> Val mutation in Sso7d. A detailed analysis of the stability of the proteins was obtained by globally fitting the variation of DSC parameters and circular dichroism intensities as a function of temperature (0-100 degrees C), salt (0-0.3 M), and pH (0-8). A global analysis over such a range of conditions permitted an unusually precise measure of the thermodynamic parameters, as well as the separation of the thermodynamics of the intrinsic unfolding reaction from the linked effects of protonation and chloride binding associated with acid-induced folding. The results indicate differences in the energetics of unfolding Sac7d and Sso7d that would not be apparent from an analysis of DSC data alone using conventional methods. The sign and magnitude of the changes in DeltaG, DeltaH, TDeltaS, and DeltaC(P) of unfolding resulting from core Ile/Val substitutions in the two proteins were consistent with differences in hydrophobicity of Val and Ile and negligible changes in packing (van der Waals) interactions. The benefit of increased hydrophobicity of the core increased with temperature, with maximal effect around 116 degrees C. Increased hydrophobicity of the core achieved not only an increase in the free energy of unfolding, but also a lateral shift of the temperature of maximal stability to higher temperature.

How Oligomerization Contributes to the Thermostability of an Archaeon Protein

To study how oligomerization may contribute to the thermostability of archaeon proteins, we focused on a hexameric protein, protein L-isoaspartyl-O-methyltransferase from Sulfolobus tokodaii (StoPIMT). The crystal structure shows that StoPIMT has a distinctive hexameric structure composed of monomers consisting of two domains: an S-adenosylmethionine-dependent methyltransferase fold domain and a C-terminal alpha-helical domain. The hexameric structure includes three interfacial contact regions: major, minor, and coiled-coil. Several C-terminal deletion mutants were constructed and characterized. The hexameric structure and thermostability were retained when the C-terminal alpha-helical domain (Tyr(206)-Thr(231)) was deleted, suggesting that oligomerization via coiled-coil association using the C-terminal alpha-helical domains did not contribute critically to hexamerization or to the increased thermostability of the protein. Deletion of three additional residues located in the major contact region, Tyr(203)-Asp(204)-Asp(205), led to a significant decrease in hexamer stability and chemico/thermostability. Although replacement of Thr(146) and Asp(204), which form two hydrogen bonds in the interface in the major contact region, with Ala did not affect hexamer formation, these mutations led to a significant decrease in thermostability, suggesting that two residues in the major contact region make significant contributions to the increase in stability of the protein via hexamerization. These results suggest that cooperative hexamerization occurs via interactions of "hot spot" residues and that a couple of interfacial hot spot residues are responsible for enhancing thermostability via oligomerization.

The influence of dipeptide composition on protein thermostability

In this work, the influence of dipeptide composition on protein thermostability was studied. After comparing the normalized dipeptide composition between mesophilic proteins and (hyper)thermophilic proteins, we concluded that when organism optimal growth temperature increased, for archaeal proteins, the compositions of VK, KI, YK, IK, KV, KY, and EV increased significantly and the compositions of DA, AD, TD, DD, DT, HD, DH, DR, and DG decreased significantly; and for bacterial proteins, the compositions of KE, EE, EK, YE, VK, KV, KK, LK, EI, EV, RK, EF, KY, VE, KI, KG, EY, FK, KF, FE, KR, VY, MK, WK, and WE increased significantly and the compositions of WQ, AA, QA, MQ, AW, QW, QQ, RQ, QH, HQ, AD, AQ, WL, QL, HA, and DA decreased significantly. So these characteristic dipeptides are correlative to protein thermostability. At the same time, the influence of single amino acid composition on protein thermostability was also studied for comparison. We found that the influence of single amino acid composition could be deduced from the influence of dipeptide composition. So we thought that the influence of dipeptide composition on protein thermostability is larger than the influence of amino acid composition. The characteristic dipeptides not only describe the dipeptides that influence protein thermostability significantly but also show the relationship among significant single amino acids that influence protein thermostability.

An evolutionary route to xylanase process fitness

Directed evolution technologies were used to selectively improve the stability of an enzyme without compromising its catalytic activity. In particular, this article describes the tandem use of two evolution strategies to evolve a xylanase, rendering it tolerant to temperatures in excess of 90 degrees C. A library of all possible 19 amino acid substitutions at each residue position was generated and screened for activity after a temperature challenge. Nine single amino acid residue changes were identified that enhanced thermostability. All 512 possible combinatorial variants of the nine mutations were then generated and screened for improved thermal tolerance under stringent conditions. The screen yielded eleven variants with substantially improved thermal tolerance. Denaturation temperature transition midpoints were increased from 61 degrees C to as high as 96 degrees C. The use of two evolution strategies in combination enabled the rapid discovery of the enzyme variant with the highest degree of fitness (greater thermal tolerance and activity relative to the wild-type parent).

Crystal structure of an ADP-dependent glucokinase from Pyrococcus furiosus: implications for a sugar-induced conformational change in ADP-dependent kinase

ADP-dependent kinases are used in the modified Embden-Meyerhoff pathway of certain archaea. Our previous study has revealed a mechanism for ADP-dependent phosphoryl transfer by Thermococcus litoralis glucokinase (tlGK), and its evolutionary relationship with ATP-dependent ribokinases and adenosine kinases (PFKB carbohydrate kinase family members). Here, we report the crystal structure of glucokinase from Pyrococcus furiosus (pfGK) in a closed conformation complexed with glucose and AMP at 1.9A resolution. In comparison with the tlGK structure, the pfGK structure shows significant conformational changes in the small domain and a region around the hinge, suggesting glucose-induced domain closing. A part of the large domain next to the hinge is also shifted accompanied with domain closing. In the pfGK structure, glucose binds in a groove between the large and small domains, and the electron density of O1 atoms for both the alpha and beta-anomer configurations was observed. The structural details of the sugar-binding site of ADP-dependent glucokinase were firstly clarified and then site-directed mutagenesis analysis clarified the catalytic residues for ADP-dependent kinase, such as Arg205 and Asp451 of tlGK. Homology search and multiple alignment of amino acid sequences using the information obtained from the structures reveals that eucaryotic hypothetical proteins homologous to ADP-dependent kinases retain the residues for the recognition of a glucose substrate.

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

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

Denaturation studies by fluorescence and quenching of thermophilic protein NAD+-glutamate dehydrogenase from Thermus thermophilus HB8

Fluorescence techniques have been used to study the structural characteristics of many proteins. The thermophilic enzyme NAD-glutamate dehydrogenase from Thermus thermophilus HB8 is found to be a hexameric enzyme. Fluorescence spectra of native and denatured protein and effect of denaturants as urea and guanidine hydrochloride on enzyme activity of thermophilic glutamate dehydrogenase (t-GDH) have been analyzed. Native t-GDH presents the maximum emission at 338 nm. The denaturation process is accompanied by an exposure to the solvent of the tryptophan residues, as manifested by the red shift of the emission maximum. Fluorescence quenching by external quenchers, KI and acrylamide, has also been carried out.

Structural plasticity of thermophilic serine hydroxymethyltransferases

Serine hydroxymethyltransferase (SHMT) catalyzes the reversible cleavage of serine to form glycine and monocarbonic groups, essential in several biosynthetic pathways. The availability of crystallographic structures of SHMT from mesophilic organisms and information produced by the genomic projects prompted the analysis of the adaptation of SHMT to "extreme" environments, such as high temperatures, by exploitation of structural data from thermophilic organisms. The sequences of 10 thermophilic/hyperthermophilic SHMTs were multiply aligned to 53 mesophilic homologs and analyzed by a comparative approach, examining the amino acid compositions and preferred residue exchanges between mesophiles and extremophiles. The structural basis of the observed exchanges was further investigated through the application of homology modeling to the 10 extremophilic SHMTs. The results of this study indicate that, in SHMT, thermal stability can be achieved mainly through three strategies: (i) increased number of charged residues at the protein surface; (ii) increased hydrophobicity of the protein core; and (iii) substitution of thermolabile residues exposed to the solvent. Additional features of the archaeal SHMTs, for which no structural data are available yet, were also investigated to explain their quaternary assemblage and the interaction with modified folates.

Comparative analysis of thermoadaptation within the archaeal glyceraldehyde-3-phosphate dehydrogenases from mesophilic Methanobacterium bryantii and thermophilic Methanothermus fervidus

To gain insight into the molecular determinants of thermoadaptation within the family of archaeal glyceraldehyde-3-phosphate dehydrogenases (GAPDH), a homology-based 3-D model of the mesophilic GAPDH from Methanobacterium bryantii was built and compared with the crystal structure of the thermophilic GAPDH from Methanothermus fervidus. The homotetrameric model of the holoenzyme was initially assembled from identical subunits completed with NADP molecules. The structure was then refined by energy minimization and simulated-annealing procedures. PROCHECK and the 3-D profile method were used to appraise the model reliability. Striking molecular features underlying the difference in stability between the enzymes were deduced from their structural comparison. First, both the increase in hydrophobic contacts and the decrease in accessibility to the protein core were shown to discriminate in favor of the thermophilic enzyme. Besides, but to a lesser degree, the number of ion pairs involved in cooperative clusters appeared to correlate with thermostability. Finally, the decreased stability of the mesophilic enzyme was also predicted to proceed from both the lack of charge-dipole interactions within alpha-helices and the enhanced entropy of unfolding due to an increase in chain flexibility. Thus, archaeal GAPDHs appear to be governed by thermoadaptation rules that differ in some aspects from those previously observed within their eubacterial counterparts.

Comparative structural analysis of psychrophilic and meso- and thermophilic enzymes

Enzymes adapted to cold display structures comparable with those of their meso- and thermophilic homologs but are characterized by a higher catalytic efficiency at low temperatures and by thermolability at moderate temperatures. To identify the structural factors responsible of such features, we undertook a systematic comparative analysis of several structural properties in a data set consisting of 7 cold active enzymes belonging to different structural families and 28 related structures from meso/thermophiles representing most of the structural information now available. Only high-resolution and high-quality structures were considered. Properties were calculated and then compared for each pair of 3D structures displaying different temperatures of adaptation using a temperature-weighting scheme. The significance of the resulting differences was evaluated with a statistical method. Results reveal that each protein family adopts different structural strategies to adapt to low temperatures. However, some common trends are observed: the number of ion pairs, the side-chain contribution to the exposed surface, and the apolar fraction of the buried surface show a consistent decrease with decreasing optimal temperatures.

High thermal stability of 3-isopropylmalate dehydrogenase from Thermus thermophilus resulting from low DeltaC(p) of unfolding

To characterize the thermal stability of 3-isopropylmalate dehydrogenase (IPMDH) from an extreme thermophile, Thermus thermophilus, urea-induced unfolding of the enzyme and of its mesophilic counterpart from Escherichia coli was investigated at various temperatures. The unfolding curves were analyzed with a three-state model for E.coli IPMDH and with a two-state model for T.thermophilus IPMDH, to obtain the free energy change DeltaG degrees of each unfolding process. Other thermodynamic parameters, enthalpy change DeltaH, entropy change DeltaS and heat capacity change DeltaC(p), were derived from the temperature dependence of DeltaG degrees. The main feature of the thermophilic enzyme was its lower dependence of DeltaG degrees on temperature resulting from a low DeltaC(p). The thermophilic IPMDH had a significantly lower DeltaC(p), 1.73 kcal/mol.K, than that of E.coli IPMDH (20.7 kcal/mol.K). The low DeltaC(p) of T.thermophilus IPMDH could not be predicted from its change in solvent-accessible surface area DeltaASA. The results suggested that there is a large structural difference between the unfolded state of T.thermophilus and that of E.coli IPMDH. Another responsible factor for the higher thermal stability of T.thermophilus IPMDH was the increase in the most stable temperature T(s). The DeltaG degrees maximum of T.thermophilus IPMDH was much smaller than that of E.coli IPMDH. The present results clearly demonstrated that a higher melting temperature T(m) is not necessarily accompanied by a higher DeltaG degrees maximum.

Thermal stability of pyrrolidone carboxyl peptidases from the hyperthermophilic Archaeon, Pyrococcus furiosus

The temperature adaptation of pyrrolidone carboxyl peptidase (PCP) from a hyperthermophile, Pyrococcus furiosus (Pf PCP), was characterized in the context of an assembly form of the protein which is a homotetramer at neutral pH. The Pf PCP exhibited maximal catalytic activity at 90-95 degrees C and its activity was higher in the temperature range 30-100 degrees C than its counterpart from the mesophilic Bacillus amyloliquefaciens (BaPCP). Thermal stability was monitored by differential scanning calorimetry (DSC). Two clearly separated peaks appeared on the DSC curves for Pf PCP at alkaline and acidic pH. Using the oxidized Pf PCP and two mutant proteins (Pf C188S and Pf C142/188S), it was found that the peaks on the high and low temperature sides of the DSC curve of Pf PCP were produced by the forms with an intersubunit disulfide bridge between the two subunits and without the bridge, respectively, indicating the stabilization effect of intersubunit disulfide bridges. The denaturation temperature (Td) of Pf PCP with intersubunit disulfide bridges was higher by 53 degrees C at pH 9.0 than that of BaPCP. An analysis of the equilibrium ultracentrifugation patterns showed that the tetrameric Pf C142/188S dissociated into dimers with decreasing pH in the acidic region and became monomer subunits at pH 2.5. The heat denaturation of Pf PCP and its two Cys mutants was highly reversible in the dimeric forms, but completely irreversible in the tetrameric form. The Td of Pf C142/188S decreased as the enzyme became dissociated, but the monomeric form of the protein was still folded at pH 2.5, although BaPCP was completely denatured at acidic pH. These results indicate that subunit interaction plays an important role in stabilizing PCP from P. furiosus in addition to the intrinsic enhanced stability of its monomer.

Kinetic and mechanistic studies of prolyl oligopeptidase from the hyperthermophile Pyrococcus furiosus

Prolyl oligopeptidase (POP) is widely distributed in mammals, where it is implicated in neuropeptide processing. It is also present in some bacteria and archaea. Because POP is found in mesophilic and hyperthermophilic organisms, and is distributed among all three phylogenetic domains, studies of its function and structure could lead to new insights about the evolution of enzyme mechanisms and thermostability. Kinetic studies were conducted on the POP of the hyperthermophilic archaeon Pyrococcus furiosus (Pfu) 85 degrees C in both H(2)O and D(2)O. Pfu POP displayed many similarities to mammalian POPs, however the solvent isotope effect (k(0)/k(1)) was 2.2 at both high and low pH, indicating that general base/acid catalysis is the rate-limiting step. The pH-rate profiles indicated a three-deprotonation process with pK(a) values of 4.3, 7.2, and 9.1. The temperature dependence of these values revealed a heat of ionization of 4.7 kJ/mol for pK(es1) and 22 kJ/mol for pK(es2), suggesting the catalytic involvement of a carboxyl group and an imidazole group, respectively. Temperature dependence of the catalytic rate was assessed at pH 6.0 and 7.6. Entropy values of -119 and -143 Jmol(-1)K(-1) were calculated at the respective pH values, with a corresponding difference in enthalpy of 8.5 kJ/mol. These values suggest that two or three hydrogen bonds are broken during the transition state of the acidic enzyme form, whereas only one or two are broken during the transition state of the basic enzyme form. A model has been constructed for Pfu POP based on the crystal structure of porcine POP and the sequence alignment. The similarities demonstrated for POPs from these two organisms reflect the most highly conserved characteristics of this class of serine protease, whereas the differences between these enzymes highlights the large evolutionary distance between them. Such fundamental information is crucial to our understanding of the function of proteins at high temperature.