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

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

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.     

Comparative Study and Mutational Analysis of Distinctive Structural Elements of Hyperthermophilic Enzymes

Comparison of the three-dimensional structure of hyperthermophilic and mesophilic β-glycosidases shows differences in secondary structure composition. The enzymes from hyperthermophilic archaea have a significantly larger number of β-strands arranged in supernumerary β-sheets compared to mesophilic enzymes from bacteria and other organisms. Amino acid replacements designed to alter the structure of the supernumerary β-strands were introduced by site directed mutagenesis into the sequence encoding the β-glycosidase from Sulfolobus solfataricus. Most of the replacements caused almost complete loss of activity but some yielded enzyme variants whose activities were affected specifically at higher temperatures. Far-UV CD spectra recorded as a function of temperature for both wild type β-glycosidase and mutant V349G, one of the mutants with reduced activity at higher temperatures, were similar, showing that the protein structure of the mutant was stable at the highest temperatures assayed. The properties of mutant V349G show a difference between thermostability (stability of the protein structure at high temperatures) and thermophilicity (optimal activity at high temperatures).     

A comparative infrared spectroscopic study of glycoside hydrolases from extremophilic archaea revealed different molecular mechanisms of adaptation to high temperatures

The identification of the determinants of protein thermal stabilization is often pursued by comparing enzymes from hyperthermophiles with their mesophilic counterparts while direct structural comparisons among proteins and enzymes from hyperthermophiles are rather uncommon. Here, oligomeric beta-glycosidases from the hyperthermophilic archaea Sulfolobus solfataricus (Ss beta-gly), Thermosphaera aggregans (Ta beta-gly), and Pyrococcus furiosus (Pf beta-gly), have been compared. Studies of FTIR spectroscopy and kinetics of thermal inactivation showed that the three enzymes had similar secondary structure composition, but Ss beta-gly and Ta beta-gly (temperatures of melting 98.1 and 98.4 degrees C, respectively) were less stable than Pf beta-gly, which maintained its secondary structure even at 99.5 degrees C. The thermal denaturation of Pf beta-gly, followed in the presence of SDS, suggested that this enzyme is stabilized by hydrophobic interactions. A detailed inspection of the 3D-structures of these enzymes supported the experimental results: Ss beta-gly and Ta beta-gly are stabilized by a combination of ion-pairs networks and intrasubunit S-S bridges while the increased stability of Pf beta-gly resides in a more compact protein core. The different strategies of protein stabilization give experimental support to recent theories on thermophilic adaptation and suggest that different stabilization strategies could have been adopted among archaea.     

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.     

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.     

Structural comparison of tRNA m1A58 methyltransferases revealed different molecular strategies to maintain their oligomeric architecture under extreme conditions

Background

tRNA m¹A58 methyltransferases (TrmI) catalyze the transfer of a methyl group from S-adenosyl-L-methionine to nitrogen 1 of adenine 58 in the T-loop of tRNAs from all three domains of life. The m¹A58 modification has been shown to be essential for cell growth in yeast and for adaptation to high temperatures in thermophilic organisms. These enzymes were shown to be active as tetramers. The crystal structures of five TrmIs from hyperthermophilic archaea and thermophilic or mesophilic bacteria have previously been determined, the optimal growth temperature of these organisms ranging from 37°C to 100°C. All TrmIs are assembled as tetramers formed by dimers of tightly assembled dimers.

Results

In this study, we present a comparative structural analysis of these TrmIs, which highlights factors that allow them to function over a large range of temperature. The monomers of the five enzymes are structurally highly similar, but the inter-monomer contacts differ strongly. Our analysis shows that bacterial enzymes from thermophilic organisms display additional intermolecular ionic interactions across the dimer interfaces, whereas hyperthermophilic enzymes present additional hydrophobic contacts. Moreover, as an alternative to two bidentate ionic interactions that stabilize the tetrameric interface in all other TrmI proteins, the tetramer of the archaeal P. abyssi enzyme is strengthened by four intersubunit disulfide bridges.

Conclusions

The availability of crystal structures of TrmIs from mesophilic, thermophilic or hyperthermophilic organisms allows a detailed analysis of the architecture of this protein family. Our structural comparisons provide insight into the different molecular strategies used to achieve the tetrameric organization in order to maintain the enzyme activity under extreme conditions.     

Structural insights into cold inactivation of tryptophanase and cold adaptation of subtilisin S41

A wide variety of enzymes can undergo a reversible loss of activity at low temperature, a process that is termed cold inactivation. This phenomenon is found in oligomeric enzymes such as tryptophanase (Trpase) and other pyridoxal phosphate dependent enzymes. On the other hand, cold-adapted, or psychrophilic enzymes, isolated from organisms able to thrive in permanently cold environments, have optimal activity at low temperature, which is associated with low thermal stability. Since cold inactivation may be considered "contradictory" to cold adaptation, we have looked into the amino acid sequences and the crystal structures of two families of enzymes, subtilisin and tryptophanase. Two cold adapted subtilisins, S41 and subtilisin-like protease from Vibrio, were compared to a mesophilic and a thermophilic subtilisins, as well as to four PLP-dependent enzymes in order to understand the specific surface residues, specific interactions, or any other molecular features that may be responsible for the differences in their tolerance to cold temperatures. The comparison between the psychrophilic and the mesophilic subtilisins revealed that the cold adapted subtilisins have a high content of acidic residues mainly found on their surface, making it charged. The analysis of the Trpases showed that they have a high content of hydrophobic residues on their surface. Thus, we suggest that the negatively charged residues on the surface of the subtilisins may be responsible for their cold adaptation, whereas the hydrophobic residues on the surface of monomeric Trpase molecules are responsible for the tetrameric assembly, and may account for their cold inactivation and dissociation.     

Ion pairs and the thermotolerance of proteins from hyperthermophiles: a "traffic rule" for hot roads

The proteins from hyperthermophilic organisms maintain their biologically active structure at temperatures that are significantly higher than the denaturation temperatures of their mesophilic counterparts. The fact that there is usually a high degree of sequence and structural homology between these two classes of proteins suggests that the source of this extreme thermal tolerance is hidden in the delicate balance of the non-covalent interactions. Among the large number of factors identified in the literature as being responsible for the thermostability of these proteins, this article focuses on electrostatic interactions. It demonstrates that the optimization of electrostatic interactions by increasing of the number of salt bridges is a driving force for enhancement of the thermotolerance of proteins from hyperthermophilic microorganisms. This feature is less evident in proteins from thermophilic organisms and is absent from mesophile-derived proteins.     

Thermophilic Fungi: Their Physiology and Enzymes

Thermophilic fungi are a small assemblage in mycota that have a minimum temperature of growth at or above 20°C and a maximum temperature of growth extending up to 60 to 62°C. As the only representatives of eukaryotic organisms that can grow at temperatures above 45°C, the thermophilic fungi are valuable experimental systems for investigations of mechanisms that allow growth at moderately high temperature yet limit their growth beyond 60 to 62°C. Although widespread in terrestrial habitats, they have remained underexplored compared to thermophilic species of eubacteria and archaea. However, thermophilic fungi are potential sources of enzymes with scientific and commercial interests. This review, for the first time, compiles information on the physiology and enzymes of thermophilic fungi. Thermophilic fungi can be grown in minimal media with metabolic rates and growth yields comparable to those of mesophilic fungi. Studies of their growth kinetics, respiration, mixed-substrate utilization, nutrient uptake, and protein breakdown rate have provided some basic information not only on thermophilic fungi but also on filamentous fungi in general. Some species have the ability to grow at ambient temperatures if cultures are initiated with germinated spores or mycelial inoculum or if a nutritionally rich medium is used. Thermophilic fungi have a powerful ability to degrade polysaccharide constituents of biomass. The properties of their enzymes show differences not only among species but also among strains of the same species. Their extracellular enzymes display temperature optima for activity that are close to or above the optimum temperature for the growth of organism and, in general, are more heat stable than those of the mesophilic fungi. Some extracellular enzymes from thermophilic fungi are being produced commercially, and a few others have commercial prospects. Genes of thermophilic fungi encoding lipase, protease, xylanase, and cellulase have been cloned and overexpressed in heterologous fungi, and pure crystalline proteins have been obtained for elucidation of the mechanisms of their intrinsic thermostability and catalysis. By contrast, the thermal stability of the few intracellular enzymes that have been purified is comparable to or, in some cases, lower than that of enzymes from the mesophilic fungi. Although rigorous data are lacking, it appears that eukaryotic thermophily involves several mechanisms of stabilization of enzymes or optimization of their activity, with different mechanisms operating for different enzymes.     

Influence of high temperature and ethanol on thermostable lignocellulolytic enzymes

Lignocellulolytic enzymes are among the most costly part in production of bioethanol. Therefore, recycling of enzymes is interesting as a concept for reduction of process costs. However, stability of the enzymes during the process is critical. In this work, focus has been on investigating the influence of temperature and ethanol on enzyme activity and stability in the distillation step, where most enzymes are inactivated due to high temperatures. Two enzyme mixtures, a mesophilic and a thermostable mixture, were exposed to typical process conditions [temperatures from 55 to 65 °C and up to 5 % ethanol (w/v)] followed by specific enzyme activity analyses and SDS-PAGE. The thermostable and mesophilic mixture remained active at up to 65 and 55 °C, respectively. When the enzyme mixtures reached their maximum temperature limit, ethanol had a remarkable influence on enzyme activity, e.g., the more ethanol, the faster the inactivation. The reason could be the hydrophobic interaction of ethanol on the tertiary structure of the enzyme protein. The thermostable mixture was more tolerant to temperature and ethanol and could therefore be a potential candidate for recycling after distillation.     

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.     

"Hot cores" in proteins: Comparative analysis of the apolar contact area in structures from hyper/thermophilic and mesophilic organisms

Background  

A wide variety of stabilizing factors have been invoked so far to elucidate the structural basis of protein thermostability. These include, amongst the others, a higher number of ion-pairs interactions and hydrogen bonds, together with a better packing of hydrophobic residues. It has been frequently observed that packing of hydrophobic side chains is improved in hyperthermophilic proteins, when compared to their mesophilic counterparts. In this work, protein crystal structures from hyper/thermophilic organisms and their mesophilic homologs have been compared, in order to quantify the difference of apolar contact area and to assess the role played by the hydrophobic contacts in the stabilization of the protein core, at high temperatures.     

Proton-Pumping Inorganic Pyrophosphatases in Some Archaea and Other Extremophilic Prokaryotes

Comparative studies between the proton-pumping, membrane-bound inorganic pyrophosphatases (H(+)-PPases) from hyperthermophilic and thermophilic prokaryotes and those from mesophilic organisms can now be performed because of very recent sequence data. Typical overall factors that contribute to protein thermostability are found in H(+)-PPases from extremophiles; nevertheless, putative active site motifs of this class of enzymes may be identical over the whole range of average growth temperatures of the compared prokaryotes. Heterologous expression in yeast of H(+)-PPases from organisms spanning a wide range of thermal habitats has allowed the biochemical comparison among these proteins within the same system, ensuring that differences observed are due to intrinsic characteristics of the proteins and not to their interactions with different cellular environments. On the other hand, the availability of H(+)-PPase sequences from a variety of sources have permitted molecular phylogenetic studies of this class of proton pumps, thus providing information about their general structural and functional properties. A great step forward may be expected when one of the several groups now attempting crystallization and 3D structural determination of H(+)-PPases will be successful.     

Activity, stability and structural studies of lactate dehydrogenases adapted to extreme thermal environments

Lactate dehydrogenase (LDH) catalyzes the conversion of pyruvate to lactate with concomitant oxidation of NADH during the last step in anaerobic glycolysis. In the present study, we present a comparative biochemical and structural analysis of various LDHs adapted to function over a large temperature range. The enzymes were from Champsocephalus gunnari (an Antarctic fish), Deinococcus radiodurans (a mesophilic bacterium) and Thermus thermophilus (a hyperthermophilic bacterium). The thermodynamic activation parameters of these LDHs indicated that temperature adaptation from hot to cold conditions was due to a decrease in the activation enthalpy and an increase in activation entropy. The crystal structures of these LDHs have been solved. Pairwise comparisons at the structural level, between hyperthermophilic versus mesophilic LDHs and mesophilic versus psychrophilic LDHs, have revealed that temperature adaptation is due to a few amino acid substitutions that are localized in critical regions of the enzyme. These substitutions, each having accumulating effects, play a role in either the conformational stability or the local flexibility or in both. Going from hot- to cold-adapted LDHs, the various substitutions have decreased the number of ion pairs, reduced the size of ionic networks, created unfavorable interactions involving charged residues and induced strong local disorder. The analysis of the LDHs adapted to extreme temperatures shed light on how evolutionary processes shift the subtle balance between overall stability and flexibility of an enzyme.     

Subunit interfaces of oligomeric hyperthermophilic enzymes display enhanced compactness

Enzymes from thermophilic and, particularly, from hyperthermophilic organisms are surprisingly stable. Understanding the molecular origin of protein thermostability and thermoactivity attracted the interest of many scientists both for the perspective comprehension of the principles of protein structure and for the possible biotechnological applications through protein engineering. Comparative studies at sequence and structure levels were aimed at detecting significant differences of structural parameters related to protein stability between thermophilic and hyperthermophilic proteins and their mesophilic homologs. In a recent work, we focused attention on structural adaptation occurring at the subunit interface of oligomeric hyper- and thermostable enzymes. A set of structural and chemico-physical parameters were compared to those observed at the corresponding interfaces of homologous mesophilic proteins. Among the most significant variations, a general increase of interface apolarity and packing density in hyperthermophilic enzymes were found. This work was therefore aimed at elucidating whether the increased packing observed is reached also through the reduction of interface cavity number and volume. The results indicate that number of cavities tends to be relatively constant while cavity volume tends to decrease in the hyperthermophilic interfaces. The cavity apolarity increases in thermophiles but, apparently, not in hyperthermophiles. Moreover, interface hot spot residues of the mesophilic interfaces tend to be conserved in the extremophilic counterparts.     

Properties of extracellular proteases from three psychrotolerant Stenotrophomonas maltophilia isolated from Antarctic soil

Three Antarctic psychrotolerant Stenotrophomonas maltophilia were isolated and the characteristics of their extracellular serine proteases were described. The isolates were able to grow at 14 and 34°C, but grew better between 20 and 28°C. The highest protease secretion was reached at 20–24°C. The purified enzyme preparations had maximal activity at 55–60°C and alkaline pH. They showed high pH stability, retaining more than 60% of residual activity after 3 h of incubation at a pH range of 4–12. The thermal stability was slightly lower compared with a commercial mesophilic protease, with 74–79% residual activity after 90 min at 40°C and 50% inactivation at 50°C between 43 and 69 min. These properties suggest that the Antarctic isolates could be adapted to cold by means of synthesising more enzymes with high activity but that the proteases they produce are not truly cold-active, being more similar to mesophilic enzymes.     

Tiny TIM: a small, tetrameric, hyperthermostable triosephosphate isomerase

Comparative structural studies on proteins derived from organisms with growth optima ranging from 15 to 100 degrees C are beginning to shed light on the mechanisms of protein thermoadaptation. One means of sustaining hyperthermostability is for proteins to exist in higher oligomeric forms than their mesophilic homologues. Triosephosphate isomerase (TIM) is one of the most studied enzymes, whose fold represents one of nature's most common protein architectures. Most TIMs are dimers of approximately 250 amino acid residues per monomer. Here, we report the 2.7 A resolution crystal structure of the extremely thermostable TIM from Pyrococcus woesei, a hyperthermophilic archaeon growing optimally at 100 degrees C, representing the first archaeal TIM structure. P. woesei TIM exists as a tetramer comprising monomers of only 228 amino acid residues. Structural comparisons with other less thermostable TIMs show that although the central beta-barrel is largely conserved, severe pruning of several helices and truncation of some loops give rise to a much more compact monomer in the small hyperthermophilic TIM. The classical TIM dimer formation is conserved in P. woesei TIM. The extreme thermostability of PwTIM appears to be achieved by the creation of a compact tetramer where two classical TIM dimers interact via an extensive hydrophobic interface. The tetramer is formed through largely hydrophobic interactions between some of the pruned helical regions. The equivalent helical regions in less thermostable dimeric TIMs represent regions of high average temperature factor. The PwTIM seems to have removed these regions of potential instability in the formation of the tetramer. This study of PwTIM provides further support for the role of higher oligomerisation states in extreme thermal stabilisation.     

Electrostatic contributions to the stability of hyperthermophilic proteins.

Electrostatic contributions to the folding free energy of several hyperthermophilic proteins and their mesophilic homologs are calculated. In all the cases studied, electrostatic interactions are more favorable in the hyperthermophilic proteins. The electrostatic free energy is found not to be correlated with the number of ionizable amino acid residues, ion pairs or ion pair networks in a protein, but rather depends on the location of these groups within the protein structure. Moreover, due to the large free energy cost associated with burying charged groups, buried ion pairs are found to be destabilizing unless they undergo favorable interactions with additional polar groups, including other ion pairs. The latter case involves the formation of stabilizing ion pair networks as is observed in a number of proteins. Ion pairs located on the protein surface also provide stabilizing interactions in a number of cases. Taken together, our results suggest that many hyperthermophilic proteins enhance electrostatic interactions through the optimum placement of charged amino acid residues within the protein structure, although different design strategies are used in different cases. Other physical mechanisms are also likely to contribute, however optimizing electrostatic interactions offers a simple means of enhancing stability without disrupting the core residues characteristic of different protein families.     

Study and design of stability in GH5 cellulases

Thermostable enzymes that hydrolyze lignocellulosic materials provide potential advantages in process configuration and enhancement of production efficiency over their mesophilic counterparts in the bioethanol industry. In this study, the dynamics of β-1,4-endoglucanases (EC: 3.2.1.4) from family 5 of glycoside hydrolases (GH5) were investigated computationally. The conformational flexibility of 12 GH5 cellulases, ranging from psychrophilic to hyperthermophilic, was investigated by molecular dynamics (MD) simulations at elevated temperatures. The results indicated that the protein flexibility and optimum activity temperatures are appreciably correlated. Intra-protein interactions, packing density and solvent accessible area were further examined in crystal structures to investigate factors that are possibly involved in higher rigidity of thermostable cellulases. The MD simulations and the rules learned from analyses of stabilizing factors were used in design of mutations toward the thermostabilization of cellulase C, one of the GH5 endoglucanases. This enzyme was successfully stabilized both chemically and thermally by introduction of a new disulfide cross-link to its highly mobile 56-amino acid subdomain.