Alkaline pH-dependent differential unfolding characteristics of mesophilic and thermophilic homologs of dimeric serine hydroxymethyltransferase

Environmental variables such as pH can significantly influence the folding and stability of a protein molecule. In the present investigation, we compared the alkaline pH-induced unfolding of two homologous serine hydroxymethyltransferase from mesophilic Bacillus subtilis (bsSHMT) and thermophilic Bacillus stearothermophilus (bstSHMT) using various biophysical techniques. The thermophilic enzyme bstSHMT was found to be more resistant to alkaline denaturation compared to its mesophilic counterpart, bsSHMT. Unfolding studies using domain-swapped chimera, constructed by swapping the C-terminal domain of these two wild-type proteins, revealed that C-terminal domain plays a pivotal role in the folding, stability and subunit interaction of these proteins. Primary amino acid sequence analysis of the proteins showed that bsSHMT has six unconserved lysine residues in C-terminal domain, which are absent in bstSHMT. Chemical modification of lysine side chains resulted in stabilization of monomers, only in case of bsSHMT. Moreover, comparison between homology model of bsSHMT with the crystal structure of bstSHMT revealed that a small stretch of 11 amino acids at the end of C-terminal domain was found protruding outside the molecule as a flexible coiled structure in bsSHMT. Taken together these findings suggest that possibly the presence of these non-identical lysine moieties and a small extension of C-terminal domain may be responsible for low stability of bsSHMT under alkaline pH condition.     

The C-terminal domain of dimeric serine hydroxymethyltransferase plays a key role in stabilization of the quaternary structure and cooperative unfolding of protein: domain swapping studies with enzymes having high sequence identity

The serine hydroxymethyltransferase from Bacillus subtilis (bsSHMT) and B. stearothermophilus (bstSHMT) are both homodimers and share approximately 77% sequence identity; however, they show very different thermal stabilities and unfolding pathways. For investigating the role of N- and C-terminal domains in stability and unfolding of dimeric SHMTs, we have swapped the structural domains between bs- and bstSHMT and generated the two novel chimeric proteins bsbstc and bstbsc, respectively. The chimeras had secondary structure, tyrosine, and pyridoxal-5'-phosphate microenvironment similar to that of the wild-type proteins. The chimeras showed enzymatic activity slightly higher than that of the wild-type proteins. Interestingly, the guanidium chloride (GdmCl)-induced unfolding showed that unlike the wild-type bsSHMT, which undergoes dissociation of native dimer into monomers at low guanidium chloride (GdmCl) concentration, resulting in a non-cooperative unfolding of enzyme, its chimera bsbstc, having the C-terminal domain of bstSHMT was resistant to low GdmCl concentration and showed a GdmCl-induced cooperative unfolding from native dimer to unfolded monomer. In contrast, the wild-type dimeric bstSHMT was resistant to low GdmCl concentration and showed a GdmCl-induced cooperative unfolding, whereas its chimera bstbsc, having the C- terminal domain of bsSHMT, showed dissociation of native dimer into monomer at low GdmCl concentration and a GdmCl-induced non-cooperative unfolding. These results clearly demonstrate that the C-terminal domain of dimeric SHMT plays a vital role in stabilization of the oligomeric structure of the native enzyme hence modulating its unfolding pathway.     

Structure and function of L-lactate dehydrogenases from thermophilic and mesophilic bacteria. I) Isolation and characterization of lactate dehydrogenases from thermophilic and mesophilic bacilli.

Lactate dehydrogenases from thermophilic bacilli (Bacillus stearothermophilus, Bacillus caldotenax) and from mesophilic bacilli (Bacillus X1, Bacillus subtilis) have been isolated by a two-step purification procedure. Only one type (LDH-P4) composed of four identical subunits (Mr 34 000 or 36 000) was found in each bacillus. The tetrameric enzymes were characterized with respect to thermostability, pH and temperature dependence of the pyruvate reduction and the L-lactate oxidation, substrate specificity, saturation kinetics (Km values of pyruvate, lactate, NAD, NADH), pyruvate and oxamate inhibition, and activation by fructose bisphosphate. The thermophilic and mesophilic enzymes differ characteristically in these parameters. Preliminary structural data (amino acid composition, comparative N-terminal sequence analysis) show the expected close phylogenetic relationship (high degree of sequence homology), but also typical differences between thermophilic and mesophilic dehydrogenases, a suitable basis for further comparative studies.     

Exploiting the Link between Protein Rigidity and Thermostability for Data‐Driven Protein Engineering

Understanding and exploiting the relationship between microscopic structure and macroscopic stability is important for developing strategies to improve protein stability at high temperatures. The thermostability of proteins has been repeatedly linked to an enhanced structural rigidity of the folded native state. In the current study, the rigidity of protein structures from mesophilic and thermophilic organisms along a thermal unfolding trajectory is directly probed. In order to perform this, protein structures were modeled as constraint networks, and the rigidity in these networks was quantified using the Floppy Inclusion and Rigid Substructure Topography (FIRST) method. During the thermal unfolding, a phase transition was observed that defines the rigidity percolation threshold and corresponds to the folded‐unfolded transition in protein folding. Using concepts from percolation theory and network science, a higher phase transition temperature was observed for ca. two‐thirds of the proteins from thermophilic organisms compared to their mesophilic counterparts, when applied to a data set of 20 pairs of homologues. From both the analysis of the microstructure of the constraint networks and monitoring the macroscopic behavior during the thermal unfolding, direct evidence was found for the “corresponding states” concept, which states that mesophilic and thermophilic enzymes are in corresponding states of similar flexibility at their respective optimal temperature. Finally, the current approach facilitated the identification of structural features from which a destabilization of the structure originates upon thermal unfolding. These predictions show a good agreement with the experimental data. Therefore, the information might be exploited in data‐driven protein engineering by pointing to residues that should be varied to obtain a protein with higher thermostability.     

Three-dimensional structures of thermophilic beta-1,4-xylanases from Chaetomium thermophilum and Nonomuraea flexuosa. Comparison of twelve xylanases in relation to their thermal stability.

The crystal structures of thermophilic xylanases from Chaetomium thermophilum and Nonomuraea flexuosa were determined at 1.75 and 2.1 A resolution, respectively. Both enzymes have the overall fold typical to family 11 xylanases with two highly twisted beta-sheets forming a large cleft. The comparison of 12 crystal structures of family 11 xylanases from both mesophilic and thermophilic organisms showed that the structures of different xylanases are very similar. The sequence identity differences correlated well with the structural differences. Several minor modifications appeared to be responsible for the increased thermal stability of family 11 xylanases: (a) higher Thr : Ser ratio (b) increased number of charged residues, especially Arg, resulting in enhanced polar interactions, and (c) improved stabilization of secondary structures involved the higher number of residues in the beta-strands and stabilization of the alpha-helix region. Some members of family 11 xylanases have a unique strategy to improve their stability, such as a higher number of ion pairs or aromatic residues on protein surface, a more compact structure, a tighter packing, and insertions at some regions resulting in enhanced interactions.     

Structural basis of the thermostability of monomeric malate synthase from a thermophilic Bacillus.

Malate synthases from a thermophilic Bacillus and Escherichia coli have been isolated in a high state of purity. Molecular weights of these two proteins determined in the native state and after denaturation in sodium dodecyl sulfate-mercaptoethanol show that the enzymes are monomeric. This conclusion is supported, for the thermophile enzyme, by the result of an electrophoretic analysis of that protein after treatment with dimethylsuberimidate and denaturation. The thermophilic Bacillus malate synthase is considerably more thermostable than its mesophilic counterparts from E. coli, Bacillus licheniformis, and Pseudomonas indigofera. It is, however, markedly labilized by an increase in the ionic strength of the medium brought about by the addition of 0.2 M potassium chloride or in pH above 9. Increased ionic strength has little effect on the thermostability of the mesophilic bacterial malate synthases. These observations provide strong support for the idea that monomeric proteins in thermophiles owe their unusual heat stability to the presence of salt bridges in their tertiary structure.     

Thermal stability of proteins does not correlate with the energy of intramolecular interactions

Small monomeric proteins from mesophilic and thermophilic organisms were studied. They have close structural and physical and chemical properties but vary in thermal stability. A thermodynamic analysis of heat unfolding was made and integral enthalpy of unfolding (DeltaH(unf)), heat capacity of hydration (DeltaC(p)(hyd)) and enthalpy of hydration (DeltaH(hyd)) and of the buried surface area (DeltaASA) of nonpolar and polar groups as well as the enthalpy of disruption of intramolecular interaction (DeltaH(int) in gas phase) at 298 K were determined. The absence of correlation between protein thermostability and energetic components suggests that regulatory mechanism of protein thermal stabilization has entropic nature.     

Comparing Residue Clusters from Thermophilic and Mesophilic Enzymes Reveals Adaptive Mechanisms

Understanding how proteins adapt to function at high temperatures is important for deciphering the energetics that dictate protein stability and folding. While multiple principles important for thermostability have been identified, we lack a unified understanding of how internal protein structural and chemical environment determine qualitative or quantitative impact of evolutionary mutations. In this work we compare equivalent clusters of spatially neighboring residues between paired thermophilic and mesophilic homologues to evaluate adaptations under the selective pressure of high temperature. We find the residue clusters in thermophilic enzymes generally display improved atomic packing compared to mesophilic enzymes, in agreement with previous research. Unlike residue clusters from mesophilic enzymes, however, thermophilic residue clusters do not have significant cavities. In addition, anchor residues found in many clusters are highly conserved with respect to atomic packing between both thermophilic and mesophilic enzymes. Thus the improvements in atomic packing observed in thermophilic homologues are not derived from these anchor residues but from neighboring positions, which may serve to expand optimized protein core regions.     

Characterization of pyridoxal 5'- phosphate-binding domain and folding intermediate of Bacillus subtilis serine hydroxymethyltransferase: an autonomous folding domain

The pyridoxal-5'-phosphate-binding domain (PLPbd) of bsSHMT (Bacillus subtilis serine hydroxymethyltransferase) was cloned and over-expressed in Escherichia coli. The recombinant protein was solublized, refolded and purified from inclusion bodies by rapid mixing followed by ion exchange chromatography. Structural and functional studies suggested the native form of the domain, which obtained as a monomer and had similar secondary and tertiary structural properties as when present in the bsSHMT. The domain also binds to the PLP however with slightly lesser affinity than the native enzyme. GdmCl (guanidium chloride)-induced equilibrium unfolding of the recombinant PLP-binding domain showed a single monophasic transition which corresponds with the second phase transition of the GdmCl-induced unfolding of bsSHMT. The results indicate that PLPbd of bsSHMT is an independent domain, which attains its tertiary structure before the dimerization of partially folded monomer and behaves as a single cooperative unfolding unit under equilibrium conditions.     

Role of entropy in protein thermostability: folding kinetics of a hyperthermophilic cold shock protein at high temperatures using 19F NMR

We used (19)F NMR to extend the temperature range accessible to detailed kinetic and equilibrium studies of a hyperthermophilic protein. Employing an optimized incorporation strategy, the small cold shock protein from the bacterium Thermotoga maritima (TmCsp) was labeled with 5-fluorotryptophan. Although chaotropically induced unfolding transitions revealed a significant decrease in the stabilization free energy upon fluorine labeling, the protein's kinetic folding mechanism is conserved. Temperature- and guanidinium chloride-dependent equilibrium unfolding transitions monitored by (19)F NMR agree well with the results from optical spectroscopy, and provide a stringent test of the two-state folding character of TmCsp. Folding and unfolding rate constants at high temperatures were determined from the (19)F NMR spectra close to the midpoint of thermal unfolding by global line shape analysis. In combination with results from stopped-flow experiments at lower temperatures, they show that the folding rate constant of TmCsp and its temperature dependence closely resemble those of its mesophilic homologue from Bacillus subtilis, BsCspB. However, the unfolding rate constant of TmCsp is two orders of magnitude lower over the entire temperature range that was investigated. Consequently, the difference in conformational stability between the two proteins is solely due to the unfolding rate constant over a wide temperature range. A thermodynamic analysis points to an important role of entropic factors in the stabilization of TmCsp relative to its mesophilic homologues.     

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: Sources, Uses, and Molecular Mechanisms for Thermostability

Enzymes synthesized by hyperthermophiles (bacteria and archaea with optimal growth temperatures of >80°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.     

Similarity and difference in the unfolding of thermophilic and mesophilic cold shock proteins studied by molecular dynamics simulations

Molecular dynamics simulations were performed to unfold a homologous pair of thermophilic and mesophilic cold shock proteins at high temperatures. The two proteins differ in just 11 of 66 residues and have very similar structures with a closed five-stranded antiparallel beta-barrel. A long flexible loop connects the N-terminal side of the barrel, formed by three strands (beta1-beta3), with the C-terminal side, formed by two strands (beta4-beta5). The two proteins were found to follow the same unfolding pathway, but with the thermophilic protein showing much slower unfolding. Unfolding started with the melting of C-terminal strands, leading to exposure of the hydrophobic core. Subsequent melting of beta3 and the beta-hairpin formed by the first two strands then resulted in unfolding of the whole protein. The slower unfolding of the thermophilic protein could be attributed to ion pair formation of Arg-3 with Glu-46, Glu-21, and the C-terminal. These ion pairs were also found to be important for the difference in folding stability between the pair of proteins. Thus electrostatic interactions appear to play similar roles in the difference in folding stability and kinetics between the pair of proteins.     

Bme585 I [5'-CCCGC(4/6)-3'], a new isoschizomer of restriction endonuclease Fau I, isolated from a strain of Bacillus mesentericus

Bme585 I is a new member of the restriction endonuclease type IIS family. It was partially purified from the heterothrophic, mesophilic bacterial strain Bacillus mesentericus 585 by ammonium sulphate precipitation and phosphocellulose column chromatography. Bme585 I is a monomeric protein with a molecular mass of 62 kD. The enzyme is active over a broad pH range from 7.0 to 8.8, has a temperature optimum of 37 degrees C and tolerance of NaCl in reaction buffer from 0 to 400 mM. Bme585 I recognizes the asymmetric sequence 5'-CCCGC(4/6)-3' and is therefore an isoschizomer of restriction endonuclease Fau I.     

Stabilization of native protein fold by intein-mediated covalent cyclization

A mutant version of the N-terminal domain of Escherichia coli DnaB helicase was used as a model system to assess the stabilization against unfolding gained by covalent cyclization. Cyclization was achieved in vivo by formation of an amide bond between the N and C termini with the help of a split mini-intein. Linear and circular proteins were constructed to be identical in amino acid sequence. Mutagenesis of Phe102 to Glu rendered the protein monomeric even at high concentration. A difference in free energy of unfolding, DeltaDeltaG, between circular and linear protein of 2.3(+/-0.5) kcal mol(-1) was measured at 10 degrees C by circular dichroism. A theoretical estimate of the difference in conformational entropy of linear and circular random chains in a three-dimensional cubic lattice model predicted DeltaDeltaG=2.3 kcal mol(-1), suggesting that stabilization by protein cyclization is driven by the reduced conformational entropy of the unfolded state. Amide-proton exchange rates measured by NMR spectroscopy and mass spectrometry showed a uniform, approximately tenfold decrease of the exchange rates of the most slowly exchanging amide protons, demonstrating that cyclization globally decreases the unfolding rate of the protein. The amide proton exchange was found to follow EX1 kinetics at near-neutral pH, in agreement with an unusually slow refolding rate of less than 4 min(-1) measured by stopped-flow circular dichroism. The linear and circular proteins differed more in their unfolding than in their folding rates. Global unfolding of the N-terminal domain of E.coli DnaB is thus promoted strongly by spatial separation of the N and C termini, whereas their proximity is much less important for folding.     

The N-terminal region is crucial for the thermostability of the G-domain of Bacillus stearothermophilus EF-Tu

Bacterial elongation factor Tu (EF-Tu) is a model monomeric G protein composed of three covalently linked domains. Previously, we evaluated the contributions of individual domains to the thermostability of EF-Tu from the thermophilic bacterium Bacillus stearothermophilus. We showed that domain 1 (G-domain) sets up the basal level of thermostability for the whole protein. Here we chose to locate the thermostability determinants distinguishing the thermophilic domain 1 from a mesophilic domain 1. By an approach of systematically swapping protein regions differing between G-domains from mesophilic Bacillus subtilis and thermophilic B. stearothermophilus, we demonstrate that a small portion of the protein, the N-terminal 12 amino acid residues, plays a key role in the thermostability of this domain. We suggest that the thermostabilizing effect of the N-terminal region could be mediated by stabilizing the functionally important effector region. Finally, we demonstrate that the effect of the N-terminal region is significant also for the thermostability of the full-length EF-Tu.     

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.     

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.     

Replacement of staphylococcal nuclease hydrophobic core residues with those from thermophilic homologues indicates packing is improved in some thermostable proteins

The importance of tight hydrophobic core packing in stabilizing proteins found in thermophilic organisms has been vigorously disputed. Here, portions of the cores found in three thermophilic homologues were transplanted into the core of staphylococcal nuclease, a protein of modest stability. Packing of the core was evaluated by comparing interaction energy of the three mutants to the comprehensive mutant library built up previously at these same sites in staphylococcal nuclease. It was found that the interaction energy of one thermophilic sequence is extraordinarily favorable and the interaction energies of other two transplanted thermophilic sequences are good, comparable to the interaction energies of mutant cores based on cores found in mesophilic homologues. As expected when transferring just a portion of the core sequence, the mutant proteins were destabilized overall relative to wild-type staphylococcal nuclease. The overall conclusion is that improvement of packing interactions is a mechanism to confer stability employed in some proteins from thermophiles, but not all.     

Discovery of Microorganisms and Enzymes Involved in High-Solids Decomposition of Rice Straw Using Metagenomic Analyses

High-solids incubations were performed to enrich for microbial communities and enzymes that decompose rice straw under mesophilic (35°C) and thermophilic (55°C) conditions. Thermophilic enrichments yielded a community that was 7.5 times more metabolically active on rice straw than mesophilic enrichments. Extracted xylanase and endoglucanse activities were also 2.6 and 13.4 times greater, respectively, for thermophilic enrichments. Metagenome sequencing was performed on enriched communities to determine community composition and mine for genes encoding lignocellulolytic enzymes. Proteobacteria were found to dominate the mesophilic community while Actinobacteria were most abundant in the thermophilic community. Analysis of protein family representation in each metagenome indicated that cellobiohydrolases containing carbohydrate binding module 2 (CBM2) were significantly overrepresented in the thermophilic community. Micromonospora, a member of Actinobacteria, primarily housed these genes in the thermophilic community. In light of these findings, Micromonospora and other closely related Actinobacteria genera appear to be promising sources of thermophilic lignocellulolytic enzymes for rice straw deconstruction under high-solids conditions. Furthermore, these discoveries warrant future research to determine if exoglucanases with CBM2 represent thermostable enzymes tolerant to the process conditions expected to be encountered during industrial biofuel production.