Structural and genomic correlates of hyperthermostability

Structural and functional analyses of an L-asparaginase from Geobacillus thermopakistaniensis

Genome sequence of Geobacillus thermopakistaniensis contains an open reading frame annotated as a type II L-asparaginase (ASNaseGt). Critical structural analysis disclosed that ASNaseGt might be a type I L-asparaginase. In order to determine whether it is a type I or type II L-asparaginase, we have performed the structural-functional characterization of the recombinant protein as well as analyzed the localization of ASNaseGt in G. thermopakistaniensis. ASNaseGt exhibited optimal activity at 52 °C and pH 9.5. There was a > 3-fold increase in activity in the presence of β-mercaptoethanol. Apparent Vmax and Km values were 2735 U/mg and 0.35 mM, respectively. ASNaseGt displayed high thermostability with >80 % residual activity even after 6 h of incubation at 55 °C. Recombinant ASNaseGt existed in oligomeric form. Addition of β-mercaptoethanol lowered the degree of oligomerization and displayed that tetrameric form was the most active, with a specific activity of 4300 U/mg. Under physiological conditions, ASNaseGt displayed >50 % of the optimal activity. Localization studies in G. thermopakistaniensis revealed that ASNaseGt is a cytosolic protein. Structural and functional characterization, and localization in G. thermopakistaniensis displayed that ASNaseGt is not a type II but a type I L-asparaginase.

Living in trinity of extremes: Genomic and proteomic signatures of halophilic, thermophilic, and pH adaptation

Since nucleic acids and proteins of unicellular prokaryotes are directly exposed to extreme environmental conditions, it is possible to explore the genomic-proteomic compositional determinants of molecular mechanisms of adaptation developed by them in response to harsh environmental conditions. Using a wealth of currently available complete genomes/proteomes we were able to explore signatures of adaptation to three environmental factors, pH, salinity, and temperature, observing major trends in compositions of their nucleic acids and proteins. We derived predictors of thermostability, halophilic, and pH adaptations and complemented them by the principal components analysis. We observed a clear difference between thermophilic and salinity/pH adaptations, whereas latter invoke seemingly overlapping mechanisms. The genome-proteome compositional trade-off reveals an intricate balance between the work of base paring and base stacking in stabilization of coding DNA and r/tRNAs, and, at the same time, universal requirements for the stability and foldability of proteins regardless of the nucleotide biases. Nevertheless, we still found hidden fingerprints of ancient evolutionary connections between the nucleotide and amino acid compositions indicating their emergence, mutual evolution, and adjustment. The evolutionary perspective on the adaptation mechanisms is further studied here by means of the comparative analysis of genomic/proteomic traits of archaeal and bacterial species. The overall picture of genomic/proteomic signals of adaptation obtained here provides a foundation for future engineering and design of functional biomolecules resistant to harsh environments.

An extensive ion-pair/hydrogen-bond network contributes to the thermostability of the MutL ATPase domain from Aquifex aeolicus

Proteins from hyperthermophiles often contain a large number of ionic interactions. Close examination of the previously determined crystal structure of the ATPase domain of MutL from a hyperthermophile, Aquifex aeolicus, revealed that the domain contains a continuous ion-pair/hydrogen-bond network consisting of 11 charged amino acid residues on a β-sheet. Mutations were introduced to disrupt the network, showing that the more extensively the network was disrupted, the greater the thermostability of the protein was decreased. Based on urea denaturation analysis, a thermodynamic parameter, energy for the conformational stability, was evaluated, which indicated that amino acid residues in the network contributed additively to the protein stability. A continuous network rather than a cluster of isolated interactions would pay less entropic penalty upon fixing the side chains to make the same number of ion pairs/hydrogen bonds, which might contribute more favorably to the structural formation of thermostable proteins.

Idiosyncratic genome evolution of the thermophilic cyanobacterium Synechococcus at the limits of phototrophy

Thermophilic microorganisms are expected to have smaller cells and genomes compared with mesophiles, a higher proportion of horizontally acquired genes, and distinct nucleotide and amino acid composition signatures. Here, we took an integrative approach to investigate these apparent correlates of thermophily for Synechococcus A/B cyanobacteria, which include the most heat-tolerant phototrophs on the planet. Phylogenomics confirmed a unique origin of different thermotolerance ecotypes, with low levels of continued gene flow between ecologically divergent but overlapping populations, which has shaped the distribution of phenotypic traits along these geothermal gradients. More thermotolerant strains do have smaller genomes, but genome reduction is associated with a decrease in community richness and metabolic diversity, rather than with cell size. Horizontal gene transfer played only a limited role during Synechococcus evolution, but, the most thermotolerant strains have acquired a Thermus tRNA modification enzyme that may stabilize translation at high temperatures. Although nucleotide base composition was not associated with thermotolerance, we found a general replacement of aspartate with glutamate, as well as a dramatic remodeling of amino acid composition at the highest temperatures that substantially differed from previous predictions. We conclude that Synechococcus A/B genome diversification largely does not conform to the standard view of temperature adaptation. In addition, carbon fixation was more thermolabile than photosynthetic oxygen evolution for the most thermotolerant strains compared with less tolerant lineages. This suggests that increased flow of reducing power generated during the light reactions to an electron sink(s) beyond carbon dioxide has emerged during temperature adaptation of these bacteria.

Investigation of thermal stability characteristic in family A DNA polymerase - A theoretical study

DNA polymerases create complementary DNA strands in living cells and are crucial to genome transmission and maintenance. These enzymes possess similar human right-handed folds which contain thumb, fingers, and palm subdomains and contribute to polymerization activities. These enzymes are classified into seven evolutionary families, A, B, C, D, X, Y, and RT, based on amino acid sequence analysis and biochemical characteristics. Family A DNA polymerases exist in an extended range of organisms including mesophilic, thermophilic, and hyper-thermophilic bacteria, participate in DNA replication and repair, and have a broad application in molecular biology and biotechnology. In this study, we attempted to detect factors that play a role in the thermostability properties of this family member despite their remarkable similarities in structure and function. For this purpose, similarities and differences in amino acid sequences, structure, and dynamics of these enzymes have been inspected. Our results demonstrated that thermophilic and hyper-thermophilic enzymes have more charged, aromatic, and polar residues than mesophilic ones and consequently show further electrostatic and cation-pi interactions. In addition, in thermophilic enzymes, aliphatic residues tend to position in buried states more than mesophilic enzymes. These residues within their aliphatic parts increase hydrophobic core packing and therefore enhance the thermostability of these enzymes. Furthermore, a decrease in thermophilic cavities volumes assists in the protein compactness enhancement. Moreover, molecular dynamic simulation results revealed that increasing temperature impacts mesophilic enzymes further than thermophilic ones that reflect on polar and aliphatic residues surface area and hydrogen bonds changes.

Thermo-L-Asparaginases: From the Role in the Viability of Thermophiles and Hyperthermophiles at High Temperatures to a Molecular Understanding of Their Thermoactivity and Thermostability

L-asparaginase (L-ASNase) is a vital enzyme with a broad range of applications in medicine, food industry, and diagnostics. Among various organisms expressing L-ASNases, thermophiles and hyperthermophiles produce enzymes with superior performances-stable and heat resistant thermo-ASNases. This review is an attempt to take a broader view on the thermo-ASNases. Here we discuss the position of thermo-ASNases in the large family of L-ASNases, their role in the heat-tolerance cellular system of thermophiles and hyperthermophiles, and molecular aspects of their thermoactivity and thermostability. Different types of thermo-ASNases exhibit specific L-asparaginase activity and additional secondary activities. All products of these enzymatic reactions are associated with diverse metabolic pathways and are important for mitigating heat stress. Thermo-ASNases are quite distinct from typical mesophilic L-ASNases based on structural properties, kinetic and activity profiles. Here we attempt to summarize the current understanding of the molecular mechanisms of thermo-ASNases' thermoactivity and thermostability, from amino acid composition to structural-functional relationships. Research of these enzymes has fundamental and biotechnological significance. Thermo-ASNases and their improved variants, cloned and expressed in mesophilic hosts, can form a large pool of enzymes with valuable characteristics for biotechnological application.

Sourcing thermotolerant poly(ethylene terephthalate) hydrolase scaffolds from natural diversity

Enzymatic deconstruction of poly(ethylene terephthalate) (PET) is under intense investigation, given the ability of hydrolase enzymes to depolymerize PET to its constituent monomers near the polymer glass transition temperature. To date, reported PET hydrolases have been sourced from a relatively narrow sequence space. Here, we identify additional PET-active biocatalysts from natural diversity by using bioinformatics and machine learning to mine 74 putative thermotolerant PET hydrolases. We successfully express, purify, and assay 51 enzymes from seven distinct phylogenetic groups; observing PET hydrolysis activity on amorphous PET film from 37 enzymes in reactions spanning pH from 4.5-9.0 and temperatures from 30-70 °C. We conduct PET hydrolysis time-course reactions with the best-performing enzymes, where we observe differences in substrate selectivity as function of PET morphology. We employed X-ray crystallography and AlphaFold to examine the enzyme architectures of all 74 candidates, revealing protein folds and accessory domains not previously associated with PET deconstruction. Overall, this study expands the number and diversity of thermotolerant scaffolds for enzymatic PET deconstruction.

iThermo: A Sequence-Based Model for Identifying Thermophilic Proteins Using a Multi-Feature Fusion Strategy

Thermophilic proteins have important application value in biotechnology and industrial processes. The correct identification of thermophilic proteins provides important information for the application of these proteins in engineering. The identification method of thermophilic proteins based on biochemistry is laborious, time-consuming, and high cost. Therefore, there is an urgent need for a fast and accurate method to identify thermophilic proteins. Considering this urgency, we constructed a reliable benchmark dataset containing 1,368 thermophilic and 1,443 non-thermophilic proteins. A multi-layer perceptron (MLP) model based on a multi-feature fusion strategy was proposed to discriminate thermophilic proteins from non-thermophilic proteins. On independent data set, the proposed model could achieve an accuracy of 96.26%, which demonstrates that the model has a good application prospect. In order to use the model conveniently, a user-friendly software package called iThermo was established and can be freely accessed at http://lin-group.cn/server/iThermo/index.html. The high accuracy of the model and the practicability of the developed software package indicate that this study can accelerate the discovery and engineering application of thermally stable proteins.

Aspartate-phobia of thermophiles as a reaction to deleterious chemical transformations

Prokaryotes growing at high temperatures have a high proportion of charged residues in their proteins to stabilize their 3D structure. By mining 175 disparate bacterial and archaeal proteomes we found that, against the general trend for charged residues, the frequency of aspartic acid residues decreases strongly as natural growth temperature increases. In search of the explanation, we hypothesized that the reason for such unusual correlation is the deleterious consequences of spontaneous chemical transformations of aspartate at high temperatures. Our subsequent statistical analysis supported this hypothesis. This finding reveals that organisms have likely adapted to high temperatures by minimizing the harmful consequences of spontaneous chemical transformations.

The Mutational Robustness of the Genetic Code and Codon Usage in Environmental Context: A Non-Extremophilic Preference?

The genetic code was evolved, to some extent, to minimize the effects of mutations. The effects of mutations depend on the amino acid repertoire, the structure of the genetic code and frequencies of amino acids in proteomes. The amino acid compositions of proteins and corresponding codon usages are still under selection, which allows us to ask what kind of environment the standard genetic code is adapted to. Using simple computational models and comprehensive datasets comprising genomic and environmental data from all three domains of Life, we estimate the expected severity of non-synonymous genomic mutations in proteins, measured by the change in amino acid physicochemical properties. We show that the fidelity in these physicochemical properties is expected to deteriorate with extremophilic codon usages, especially in thermophiles. These findings suggest that the genetic code performs better under non-extremophilic conditions, which not only explains the low substitution rates encountered in halophiles and thermophiles but the revealed relationship between the genetic code and habitat allows us to ponder on earlier phases in the history of Life.

Phylogenetic, functional and structural characterization of a GH10 xylanase active at extreme conditions of temperature and alkalinity

Endoxylanases active under extreme conditions of temperature and alkalinity can replace the use of highly pollutant chemicals in the pulp and paper industry. Searching for enzymes with these properties, we carried out a comprehensive bioinformatics study of the GH10 family. The phylogenetic analysis allowed the construction of a radial cladogram in which protein sequences putatively ascribed as thermophilic and alkaliphilic appeared grouped in a well-defined region of the cladogram, designated TAK Cluster. One among five TAK sequences selected for experimental analysis (Xyn11) showed extraordinary xylanolytic activity under simultaneous conditions of high temperature (90 °C) and alkalinity (pH 10.5). Addition of a carbohydrate binding domain (CBM2) at the C-terminus of the protein sequence further improved the activity of the enzyme at high pH. Xyn11 structure, which has been solved at 1.8 Å resolution by X-ray crystallography, reveals an unusually high number of hydrophobic, ionic and hydrogen bond atomic interactions that could account for the enzyme’s extremophilic nature.

Molecular strategies to enhance stability and catalysis of extremophile-derived α-amylase using computational biology

α-Amylase is the most significant glycoside hydrolase having applications in various industries. It cleaves the α,1-4 glucosidic linkages of polysaccharides like starch, glycogen to yield a small polymer of glucose in α-anomeric configuration. α-Amylase is produced by all the three domains of life but microorganisms are preferred sources for industrial-scale production due to several advantages. Enormous studies and research have been done in this field in the past few decades. Still, it is requisite to work on enzyme stability and catalysis, as it loses its functionality in extreme. As the enzyme loses its structural and catalytic property under extreme environmental conditions, it is mandatory to confer some potential strategies for enhancing enzyme behaviour in such conditions. This limitation of an enzyme can be overcome up to some extent by extremophiles. They serve as an excellent source of α-amylase with outstanding features. This review is an attempt to encapsulate some structure-based strategies for improving enzyme behaviour thereby enabling researchers to selectively amend any of the strategies as per requirement during upstream and downstream processing for higher enzyme yield and stability. Thus, it will provide some cutting-edge strategies for tailoring α-amylase producing organism and enzyme with the help of several computational biology tools.

Threonine-77 Is a Determinant of the Low-Temperature Conditioning of the Most Abundant Isoform of Tropomyosin in Atlantic Salmon

The Atlantic salmon Salmo salar survives below 10 °C. The main skeletal muscle is composed of a single isoform of tropomyosin (classified as Tpm1 α-fast) that is >92% identical to the mammalian homologue. How salmon Tpm1 maintains flexibility is investigated by reversing the only full charge substitution; threonine-77(g) in salmon and lysine in other vertebrates. The mutation (Thr-77 to Lys), which falls within a known destabilizing alanine cluster, (i) yields a useful electrophoretic shift in the absence and presence of an anionic detergent, (ii) increases the T_ms of both cooperative transitions (calorimetry, 0.1 M salt, pH 7) [35 °C (minor) and 44 °C (major); ΔT_m1 = 5 °C, ΔT_m2 = 3.5 °C], (iii) increases the T_m of CN1A (residues 11-127) to 53 °C (ΔT_m = 13 °C), a value similar to that of mammalian CN1A, (iv) markedly reduces the rate of proteolysis at Leu-169, and (v) weakens the affinity of salmon Tpm1 for troponin-Sepharose. Glu-82(e), the interstrand ionic partner of Lys-77(g), is conserved. The change in ionic interactions at this locus is postulated to be "sensed" in actin period 5 (residues 166-207) and likely beyond. Wild type (acetylated) salmon Tmp1 binds more tightly to F-actin at 4 °C than at 22 °C, which is the opposite of the long-known relationship displayed by the mammalian homologue. All of the evidence indicates that the presence of a neutral 77th amino acid destabilizes a sizable portion of salmon Tpm1 that includes the midregion. Threonine-77 is a key factor in rescuing the thin filament from the peril of cold-induced rigidity.

Predicting the optimal growth temperatures of prokaryotes using only genome derived features

MOTIVATION

Optimal growth temperature is a fundamental characteristic of all living organisms. Knowledge of this temperature is central to the study of a prokaryote, the thermal stability and temperature dependent activity of its genes, and the bioprospecting of its genome for thermally adapted proteins. While high throughput sequencing methods have dramatically increased the availability of genomic information, the growth temperatures of the source organisms are often unknown. This limits the study and technological application of these species and their genomes. Here, we present a novel method for the prediction of growth temperatures of prokaryotes using only genomic sequences.

RESULTS

By applying the reverse ecology principle that an organism's genome includes identifiable adaptations to its native environment, we can predict a species' optimal growth temperature with an accuracy of 5.17°C root-mean-square error and a coefficient of determination of 0.835. The accuracy can be further improved for specific taxonomic clades or by excluding psychrophiles. This method provides a valuable tool for the rapid calculation of organism growth temperature when only the genome sequence is known.

AVAILABILITY AND IMPLEMENTATION

Source code, genomes analyzed and features calculated are available at: https://github.com/DavidBSauer/OGT_prediction.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Protein adaptation to high hydrostatic pressure: Computational analysis of the structural proteome

Hydrostatic pressure has a vital role in the biological adaptation of the piezophiles, organisms that live under high hydrostatic pressure. However, the mechanisms by which piezophiles are able to adapt their proteins to high hydrostatic pressure is not well understood. One proposed hypothesis is that the volume changes of unfolding (ΔVTot ) for proteins from piezophiles is distinct from those of nonpiezophilic organisms. Since ΔVTot defines pressure dependence of stability, we performed a comprehensive computational analysis of this property for proteins from piezophilic and nonpiezophilic organisms. In addition, we experimentally measured the ΔVTot of acylphosphatases and thioredoxins belonging to piezophilic and nonpiezophilic organisms. Based on this analysis we concluded that there is no difference in ΔVTot for proteins from piezophilic and nonpiezophilic organisms. Finally, we put forward the hypothesis that increased concentrations of osmolytes can provide a systemic increase in pressure stability of proteins from piezophilic organisms and provide experimental thermodynamic evidence in support of this hypothesis.

High-throughput quantification of protein structural change reveals potential mechanisms of temperature adaptation in Mytilus mussels

Background

Temperature exerts a strong influence on protein evolution: species living in thermally distinct environments often exhibit adaptive differences in protein structure and function. However, previous research on protein temperature adaptation has focused on small numbers of proteins and on proteins adapted to extreme temperatures. Consequently, less is known about the types and quantity of evolutionary change that occurs to proteins when organisms adapt to small shifts in environmental temperature. In this study, these uncertainties were addressed by developing software that enabled comparison of structural changes associated with temperature adaptation (hydrogen bonding, salt bridge formation, and amino acid use) among large numbers of proteins from warm- and cold-adapted species of marine mussels, Mytilus galloprovincialis and Mytilus trossulus, respectively.

Results

Small differences in habitat temperature that characterize the evolutionary history of Mytilus mussels were sufficient to cause protein structural changes consistent with temperature adaptation. Hydrogen bonds and salt bridges that increase stability and protect against heat-induced denaturation were more abundant in proteins from warm-adapted M. galloprovincialis compared with proteins from cold-adapted M. trossulus. These structural changes were related to deviations in the use of polar and charged amino acids that facilitate formation of hydrogen bonds and salt bridges within proteins, respectively. Enzymes, in particular those within antioxidant and cell death pathways, were over-represented among proteins with the most hydrogen bonds and salt bridges in warm-adapted M. galloprovincialis. Unlike extremophile proteins, temperature adaptation in Mytilus proteins did not involve substantial changes in the number of hydrophobic or large volume amino acids, nor in the content of glycine or proline.

Conclusions

Small shifts in organism temperature tolerance, such as that needed to cope with climate warming, may result from structural and functional changes to a small percentage of the proteome. Proteins in which function is dependent on large conformational change, notably enzymes, may be particularly sensitive to temperature perturbation and represent foci for natural selection. Protein temperature adaptation can occur through different types and frequencies of structural change, and adaptive mechanisms used to cope with small shifts in habitat temperature appear different from mechanisms used to retain protein function at temperature extremes.

Misprediction of Structural Disorder in Halophiles

Whereas the concept of intrinsic disorder derives from biophysical observations of the lack of structure of proteins or protein regions under native conditions, many of our respective concepts rest on proteome-scale bioinformatics predictions. It is established that most predictors work reliably on proteins commonly encountered, but it is often neglected that we know very little about their performance on proteins of microorganisms that thrive in environments of extreme temperature, pH, or salt concentration, which may cause adaptive sequence composition bias. To address this issue, we predicted structural disorder for the complete proteomes of different extremophile groups by popular prediction methods and compared them to those of the reference mesophilic group. While significant deviations from mesophiles could be explained by a lack or gain of disordered regions in hyperthermophiles and radiotolerants, respectively, we found systematic overprediction in the case of halophiles. Additionally, examples were collected from the Protein Data Bank (PDB) to demonstrate misprediction and to help understand the underlying biophysical principles, i.e., halophilic proteins maintain a highly acidic and hydrophilic surface to avoid aggregation in high salt conditions. Although sparseness of data on disordered proteins from extremophiles precludes the development of dedicated general predictors, we do formulate recommendations for how to address their disorder with current bioinformatics tools.

Net Evolutionary Loss of Residue Polarity in Drosophilid Protein Cores Indicates Ongoing Optimization of Amino Acid Composition

Amino acid frequencies in proteins may not be at equilibrium. We consider two possible explanations for the nonzero net residue fluxes in drosophilid proteins. First, protein interiors may have a suboptimal residue composition and be under a selective pressure favoring stability, that is, leading to the loss of polar (and the gain of large) amino acids. One would then expect stronger net fluxes on the protein interior than at the exposed sites. Alternatively, if most of the polarity loss occurs at the exposed sites and the selective constraint on amino acid composition at such sites decreases over time, net loss of polarity may be neutral and caused by disproportionally high occurrence of polar residues at exposed, least constrained sites. We estimated net evolutionary fluxes of residue polarity and volume at sites with different solvent accessibility in conserved protein families from 12 species of Drosophila. Net loss of polarity, miniscule in magnitude, but consistent across all lineages, occurred at all sites except the most exposed ones, where net flux of polarity was close to zero or, in membrane proteins, even positive. At the intermediate solvent accessibility the net fluxes of polarity and volume were similar to neutral predictions, whereas much of the polarity loss not attributable to neutral expectations occurred at the buried sites. These observations are consistent with the hypothesis that residue composition in many proteins is structurally suboptimal and continues to evolve toward lower polarity in the protein interior, in particular in proteins with intracellular localization. The magnitude of polarity and volume changes was independent from the protein’s evolutionary age, indicating that the approach to equilibrium has been slow or that no such single equilibrium exists.

Arabidopsis Heat Stress-Induced Proteins Are Enriched in Electrostatically Charged Amino Acids and Intrinsically Disordered Regions

Comparison of the proteins of thermophilic, mesophilic, and psychrophilic prokaryotes has revealed several features characteristic to proteins adapted to high temperatures, which increase their thermostability. These characteristics include a profusion of disulfide bonds, salt bridges, hydrogen bonds, and hydrophobic interactions, and a depletion in intrinsically disordered regions. It is unclear, however, whether such differences can also be observed in eukaryotic proteins or when comparing proteins that are adapted to temperatures that are more subtly different. When an organism is exposed to high temperatures, a subset of its proteins is overexpressed (heat-induced proteins), whereas others are either repressed (heat-repressed proteins) or remain unaffected. Here, we determine the expression levels of all genes in the eukaryotic model system Arabidopsis thaliana at 22 and 37 °C, and compare both the amino acid compositions and levels of intrinsic disorder of heat-induced and heat-repressed proteins. We show that, compared to heat-repressed proteins, heat-induced proteins are enriched in electrostatically charged amino acids and depleted in polar amino acids, mirroring thermophile proteins. However, in contrast with thermophile proteins, heat-induced proteins are enriched in intrinsically disordered regions, and depleted in hydrophobic amino acids. Our results indicate that temperature adaptation at the level of amino acid composition and intrinsic disorder can be observed not only in proteins of thermophilic organisms, but also in eukaryotic heat-induced proteins; the underlying adaptation pathways, however, are similar but not the same.

Molecular improvements in microbial α-amylases for enhanced stability and catalytic efficiency

α-Amylases is one of the most important industrial enzyme which contributes to 25% of the industrial enzyme market. Though it is produced by plant, animals and microbial source, those from microbial source seems to have potential applications due to their stability and economic viability. However a large number of α-amylases from different sources have been detailed in the literature, only few numbers of them could withstand the harsh industrial conditions. Thermo-stability, pH tolerance, calcium independency and oxidant stability and starch hydrolyzing efficiency are the crucial qualities for α-amylase in starch based industries. Microbes can be genetically modified and fine tuning can be done for the production of enzymes with desired characteristics for specific applications. This review focuses on the native and recombinant α-amylases from microorganisms, their heterologous production and the recent molecular strategies which help to improve the properties of this industrial enzyme.

Determinants of thermostability in the cytochrome P450 fold

Cytochromes P450 are found throughout the biosphere in a wide range of environments, serving a multitude of physiological functions. The ubiquity of the P450 fold suggests that it has been co-opted by evolution many times, and likely presents a useful compromise between structural stability and conformational flexibility. The diversity of substrates metabolized and reactions catalyzed by P450s makes them attractive starting materials for use as biocatalysts of commercially useful reactions. However, process conditions impose different requirements on enzymes to those in which they have evolved naturally. Most natural environments are relatively mild, and therefore most P450s have not been selected in Nature for the ability to withstand temperatures above ~40°C, yet industrial processes frequently require extended incubations at much higher temperatures. Thus, there has been considerable interest and effort invested in finding or engineering thermostable P450 systems. Numerous P450s have now been identified in thermophilic organisms and analysis of their structures provides information as to mechanisms by which the P450 fold can be stabilized. In addition, protein engineering, particularly by directed or artificial evolution, has revealed mutations that serve to stabilize particular mesophilic enzymes of interest. Here we review the current understanding of thermostability as it applies to the P450 fold, gleaned from the analysis of P450s characterized from thermophilic organisms and the parallel engineering of mesophilic forms for greater thermostability. We then present a perspective on how this information might be used to design stable P450 enzymes for industrial application. This article is part of a Special Issue entitled: Cytochrome P450 biodiversity and biotechnology, edited by Erika Plettner, Gianfranco Gilardi, Luet Wong, Vlada Urlacher, Jared Goldstone.

Comparative analysis of microbial diversity in two hot springs of Bakreshwar, West Bengal, India

Various aspects of hot springs at Bakreshwar (Lat. 23°52'48″N; Long. 87°22'40″E) in West Bengal, India have been investigated since the middle of 20th century, but comprehending the complete diversity and the complexity of the microbial population therein has been in the continuing process. Some of these microorganisms are found to have immense industrial importance. Microbes generally exist in milieus of varying complexities and diversities. Attempting the usually employed cultivation-based techniques in experimentation with those microbes had confronted various limitations. To overcome these limitations a strategy based on high-throughput sequencing of 16S rRNA gene amplicon analysis was employed for studying the differential diversity and the detailed nature of microbial population of the two hot springs of Bakreshwar (54 °C & 65 °C). Paired-end libraries of amplified V-3 hyper-variable 16S rDNA fragments from sets of samples that varied in their contents, ranging from a single bacterium to highly complex communities were sequenced. The comparison revealed the differential aspects in the two hot spring waters; the samples at 54 °C showed the bacterial phylum Firmicutes (65.85%) and Synergistetes (27.24%) predominating and those from hot spring water at 65 °C showed the abundance of the phyla Firmicutes (96.10%) and Proteobacteria (3.36%). The presence of Archaea in the hot springs could not be ascertained.

Insights into thermoadaptation and the evolution of mesophily from the bacterial phylum Thermotogae

Thermophiles are extremophiles that grow optimally at temperatures >45 °C. To survive and maintain function of their biological molecules, they have a suite of characteristics not found in organisms that grow at moderate temperature (mesophiles). At the cellular level, thermophiles have mechanisms for maintaining their membranes, nucleic acids, and other cellular structures. At the protein level, each of their proteins remains stable and retains activity at temperatures that would denature their mesophilic homologs. Conversely, cellular structures and proteins from thermophiles may not function optimally at moderate temperatures. These differences between thermophiles and mesophiles presumably present a barrier for evolutionary transitioning between the 2 lifestyles. Therefore, studying closely related thermophiles and mesophiles can help us determine how such lifestyle transitions may happen. The bacterial phylum Thermotogae contains hyperthermophiles, thermophiles, mesophiles, and organisms with temperature ranges wide enough to span both thermophilic and mesophilic temperatures. Genomic, proteomic, and physiological differences noted between other bacterial thermophiles and mesophiles are evident within the Thermotogae. We argue that the Thermotogae is an ideal group of organisms for understanding of the response to fluctuating temperature and of long-term evolutionary adaptation to a different growth temperature range.

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.

Mining the entire Protein DataBank for frequent spatially cohesive amino acid patterns

BACKGROUND

The three-dimensional structure of a protein is an essential aspect of its functionality. Despite the large diversity in protein structures and functionality, it is known that there are common patterns and preferences in the contacts between amino acid residues, or between residues and other biomolecules, such as DNA. The discovery and characterization of these patterns is an important research topic within structural biology as it can give fundamental insight into protein structures and can aid in the prediction of unknown structures.

RESULTS

Here we apply an efficient spatial pattern miner to search for sets of amino acids that occur frequently in close spatial proximity in the protein structures of the Protein DataBank. This allowed us to mine for a new class of amino acid patterns, that we term FreSCOs (Frequent Spatially Cohesive Component sets), which feature synergetic combinations. To demonstrate the relevance of these FreSCOs, they were compared in relation to the thermostability of the protein structure and the interaction preferences of DNA-protein complexes. In both cases, the results matched well with prior investigations using more complex methods on smaller data sets.

CONCLUSIONS

The currently characterized protein structures feature a diverse set of frequent amino acid patterns that can be related to the stability of the protein molecular structure and that are independent from protein function or specific conserved domains.

Long-Range PCR Amplification of DNA by DNA Polymerase III Holoenzyme from Thermus thermophilus

DNA replication in bacteria is accomplished by a multicomponent replicase, the DNA polymerase III holoenzyme (pol III HE). The three essential components of the pol III HE are the α polymerase, the β sliding clamp processivity factor, and the DnaX clamp-loader complex. We report here the assembly of the functional holoenzyme from Thermus thermophilus (Tth), an extreme thermophile. The minimal holoenzyme capable of DNA synthesis consists of α, β and DnaX (τ and γ), δ and δ′ components of the clamp-loader complex. The proteins were each cloned and expressed in a native form. Each component of the system was purified extensively. The minimum holoenzyme from these five purified subunits reassembled is sufficient for rapid and processive DNA synthesis. In an isolated form the α polymerase was found to be unstable at temperatures above 65°C. We were able to increase the thermostability of the pol III HE to 98°C by addition and optimization of various buffers and cosolvents. In the optimized buffer system we show that a replicative polymerase apparatus, Tth pol III HE, is capable of rapid amplification of regions of DNA up to 15,000 base pairs in PCR reactions.

Characterization of a novel thermostable N-acylhomoserine lactonase from the thermophilic bacterium Thermaerobacter marianensis

Thermaerobacter marianensis is an extremely thermophilic bacterium, which was isolated from the Mariana Trench, with an optimal growth temperature of approximately 75 °C. N-Acylhomoserine lactone (AHL) is a quorum-sensing signal molecule used by many gram-negative bacteria. Here, we report the identification of an AHL-degrading gene homolog (designated aiiT) in the genome of T. marianensis JCM 10246. AiiT has 59.7%, 21.2%, and 11.2% identity to AhlS from Solibacillus silvestris, AiiA from Bacillus cereus, and AidC from Chryseobacterium sp., respectively. Homologs of aiiT were also found in Thermaerobacter nagasakiensis, T. composti, and T. subterraneus. A purified AiiT-maltose binding fusion showed high AHL-degrading activity against N-hexanoyl-L-homoserine lactone, N-octanoyl-L-homoserine lactone, and N-decanoyl-L-homoserine lactone at temperatures ranging from 40 to 80 °C. HPLC analysis revealed that AiiT functions as an AHL-lactonase that catalyzes AHL ring opening by hydrolyzing lactones. AiiT displayed maximal activity at high temperatures (60-80 °C) and showed higher thermostability than other AHL lactonases.

The fundamental tradeoff in genomes and proteomes of prokaryotes established by the genetic code, codon entropy, and physics of nucleic acids and proteins

Background

Mutations in nucleotide sequences provide a foundation for genetic variability, and selection is the driving force of the evolution and molecular adaptation. Despite considerable progress in the understanding of selective forces and their compositional determinants, the very nature of underlying mutational biases remains unclear.

Results

We explore here a fundamental tradeoff, which analytically describes mutual adjustment of the nucleotide and amino acid compositions and its possible effect on the mutational biases. The tradeoff is determined by the interplay between the genetic code, optimization of the codon entropy, and demands on the structure and stability of nucleic acids and proteins.

Conclusion

The tradeoff is the unifying property of all prokaryotes regardless of the differences in their phylogenies, life styles, and extreme environments. It underlies mutational biases characteristic for genomes with different nucleotide and amino acid compositions, providing foundation for evolution and adaptation.

Reviewers

This article was reviewed by Eugene Koonin, Michael Gromiha, and Alexander Schleiffer.

Electronic supplementary material

The online version of this article (doi:10.1186/s13062-014-0029-2) contains supplementary material, which is available to authorized users.

Stability curve prediction of homologous proteins using temperature-dependent statistical potentials

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

The prediction of protein stability remains one of the key goals of protein science. Despite the significant efforts of the last decades, faster and more accurate stability predictors on the proteomic-wide scale are currently demanded. The determination and control of protein stability are indeed fundamental steps on the path towards de novo design. In this paper we develop a method for predicting the stability curve of proteins. This curve encodes the temperature dependence of the folding free energy (). Its knowledge is important in the study of protein stability since all the thermodynamic parameters characterizing the folding transition can be extracted from it. Our prediction method is based on temperature-dependent mean force potentials and uses the tertiary structure of the target protein as well as the melting temperature () and the heat capacity change () of some other proteins belonging to the same family. From the predicted stability curves, the , the and the at room temperature can be inferred. The predictions obtained are compared with experimental data and show reasonable performances.

A loose domain swapping organization confers a remarkable stability to the dimeric structure of the arginine binding protein from Thermotoga maritima

The arginine binding protein from Thermatoga maritima (TmArgBP), a substrate binding protein (SBP) involved in the ABC system of solute transport, presents a number of remarkable properties. These include an extraordinary stability to temperature and chemical denaturants and the tendency to form multimeric structures, an uncommon feature among SBPs involved in solute transport. Here we report a biophysical and structural characterization of the TmArgBP dimer. Our data indicate that the dimer of the protein is endowed with a remarkable stability since its full dissociation requires high temperature as well as SDS and urea at high concentrations. In order to elucidate the atomic level structural properties of this intriguing protein, we determined the crystallographic structures of the apo and the arginine-bound forms of TmArgBP using MAD and SAD methods, respectively. The comparison of the liganded and unliganded models demonstrates that TmArgBP tertiary structure undergoes a very large structural re-organization upon arginine binding. This transition follows the Venus Fly-trap mechanism, although the entity of the re-organization observed in TmArgBP is larger than that observed in homologous proteins. Intriguingly, TmArgBP dimerizes through the swapping of the C-terminal helix. This dimer is stabilized exclusively by the interactions established by the swapping helix. Therefore, the TmArgBP dimer combines a high level of stability and conformational freedom. The structure of the TmArgBP dimer represents an uncommon example of large tertiary structure variations amplified at quaternary structure level by domain swapping. Although the biological relevance of the dimer needs further assessments, molecular modelling suggests that the two TmArgBP subunits may simultaneously interact with two distinct ABC transporters. Moreover, the present protein structures provide some clues about the determinants of the extraordinary stability of the biomolecule. The availability of an accurate 3D model represents a powerful tool for the design of new TmArgBP suited for biotechnological applications.

Molecular mechanisms of adaptation emerging from the physics and evolution of nucleic acids and proteins

DNA, RNA and proteins are major biological macromolecules that coevolve and adapt to environments as components of one highly interconnected system. We explore here sequence/structure determinants of mechanisms of adaptation of these molecules, links between them, and results of their mutual evolution. We complemented statistical analysis of genomic and proteomic sequences with folding simulations of RNA molecules, unraveling causal relations between compositional and sequence biases reflecting molecular adaptation on DNA, RNA and protein levels. We found many compositional peculiarities related to environmental adaptation and the life style. Specifically, thermal adaptation of protein-coding sequences in Archaea is characterized by a stronger codon bias than in Bacteria. Guanine and cytosine load in the third codon position is important for supporting the aerobic life style, and it is highly pronounced in Bacteria. The third codon position also provides a tradeoff between arginine and lysine, which are favorable for thermal adaptation and aerobicity, respectively. Dinucleotide composition provides stability of nucleic acids via strong base-stacking in ApG dinucleotides. In relation to coevolution of nucleic acids and proteins, thermostability-related demands on the amino acid composition affect the nucleotide content in the second codon position in Archaea.

Uniquely localized intra-molecular amino acid concentrations at the glycolytic enzyme catalytic/active centers of Archaea, Bacteria and Eukaryota are associated with their proposed temporal appearances on earth

The distributions of amino acids at most-conserved sites nearest catalytic/active centers (C/AC) in 4,645 sequences of ten enzymes of the glycolytic Embden-Meyerhof-Parnas pathway in Archaea, Bacteria and Eukaryota are similar to the proposed temporal order of their appearance on Earth. Glycine, isoleucine, leucine, valine, glutamic acid and possibly lysine often described as prebiotic, i.e., existing or occurring before the emergence of life, were localized in positional and conservational defined aggregations in all enzymes of all Domains. The distributions of all 20 biologic amino acids in most-conserved sites nearest their C/ACs were quite different either from distributions in sites less-conserved and further from their C/ACs or from all amino acids regardless of their position or conservation. The major concentrations of glycine, e.g., perhaps the earliest prebiotic amino acid, occupies ≈ 16 % of all the most-conserved sites within a volume of ≈ 7-8 Å radius from their C/ACs and decreases linearly towards the molecule's peripheries. Spatially localized major concentrations of isoleucine, leucine and valine are in the mid-conserved and mid-distant sites from their C/ACs in protein interiors. Lysine and glutamic acid comprise ≈ 25-30 % of all amino acids within an irregular volume bounded by ≈ 24-28 Å radii from their C/ACs at the most-distant least-conserved sites. The unreported characteristics of these amino acids: their spatially and conservationally identified concentrations in Archaea, Bacteria and Eukaryota, suggest some common structural organization of glycolytic enzymes that may be relevant to their evolution and that of other proteins. We discuss our data in relation to enzyme evolution, their reported prebiotic putative temporal appearances on Earth, abundances, biological "cost", neighbor-sequence preferences or "ordering" and some thermodynamic parameters.

Structural differences between thermophilic and mesophilic membrane proteins

The evolutionary adaptations of thermophilic water-soluble proteins required for maintaining stability at high temperature have been extensively investigated. Little is known about the adaptations in membrane proteins, however. Here, we compare many properties of mesophilic and thermophilic membrane protein structures, including side-chain burial, packing, hydrogen bonding, transmembrane kinks, loop lengths, hydrophobicity, and other sequence features. Most of these properties are quite similar between mesophiles and thermophiles although we observe a slight increase in side-chain burial and possibly a slight decrease in the frequency of transmembrane kinks in thermophilic membrane protein structures. The most striking difference is the increased hydrophobicity of thermophilic transmembrane helices, possibly reflecting more stringent hydrophobicity requirements for membrane partitioning at high temperature. In agreement with prior work examining transmembrane sequences, we find that thermophiles have an increase in small residues (Gly, Ala, Ser, and Val) and a strong suppression of Cys. We also find a relative dearth of most strongly polar residues (Asp, Asn, Glu, Gln, and Arg). These results suggest that in thermophiles, there is significant evolutionary pressure to offload destabilizing polar amino acids, to decrease the entropy cost of side chain burial, and to eliminate thermally sensitive amino acids.

Deep transcriptome-sequencing and proteome analysis of the hydrothermal vent annelid Alvinella pompejana identifies the CvP-bias as a robust measure of eukaryotic thermostability

BACKGROUND

Alvinella pompejana is an annelid worm that inhabits deep-sea hydrothermal vent sites in the Pacific Ocean. Living at a depth of approximately 2500 meters, these worms experience extreme environmental conditions, including high temperature and pressure as well as high levels of sulfide and heavy metals. A. pompejana is one of the most thermotolerant metazoans, making this animal a subject of great interest for studies of eukaryotic thermoadaptation.

RESULTS

In order to complement existing EST resources we performed deep sequencing of the A. pompejana transcriptome. We identified several thousand novel protein-coding transcripts, nearly doubling the sequence data for this annelid. We then performed an extensive survey of previously established prokaryotic thermoadaptation measures to search for global signals of thermoadaptation in A. pompejana in comparison with mesophilic eukaryotes. In an orthologous set of 457 proteins, we found that the best indicator of thermoadaptation was the difference in frequency of charged versus polar residues (CvP-bias), which was highest in A. pompejana. CvP-bias robustly distinguished prokaryotic thermophiles from prokaryotic mesophiles, as well as the thermophilic fungus Chaetomium thermophilum from mesophilic eukaryotes. Experimental values for thermophilic proteins supported higher CvP-bias as a measure of thermal stability when compared to their mesophilic orthologs. Proteome-wide mean CvP-bias also correlated with the body temperatures of homeothermic birds and mammals.

CONCLUSIONS

Our work extends the transcriptome resources for A. pompejana and identifies the CvP-bias as a robust and widely applicable measure of eukaryotic thermoadaptation.

DRE-binding transcription factor gene (LlaDREB1b) is regulated by various abiotic stresses in Lepidium latifolium L

We have isolated and in silico characterized a cold-induced LlaDREB1b encoding a putative DRE-binding transcription factor from Lepidium latifolium. Its cDNA (JN214345) sequence (998 bp) consisted of a 642 bp ORF, 168 and 188 bp of 5' and 3' UTR regions, respectively, encoding a protein of 213 aa with deduced molecular mass 23.85 kDa and pI of 4.63. In silico and phylogenetic analysis further suggested that the protein showed features of a typical member of the AP2/EREBP family of DNA-binding proteins. Southern blot analysis indicated that multiple copies of the gene could be present in the genome. Its transcripts were abundant in leaves, petiole and stem, but scarce in roots and could strongly be induced by cold treatment (4 °C), weakly by drought and salt stress, and did not respond to ABA treatment. Thus, LlaDREB1b is a potential candidate for abiotic stress-tolerance engineering in crop plants upon its further functional validation.

Preservation and evolution of organic matter during experimental fossilisation of the hyperthermophilic archaea Methanocaldococcus jannaschii

Identification of the earliest traces of life is made difficult by the scarcity of the preserved microbial remains and by the alteration and potential contamination of the organic matter (OM) content of rocks. These factors can confuse interpretations of the biogenicity and syngenicity of fossilised structures and organic molecules found in ancient rocks. In order to improve our knowledge of the fossilisation processes and their effects at the molecular level, we made a preliminary study of the fate of OM during experimental fossilisation. Changes in the composition and quantity of amino acids, monosaccharides and fatty acids were followed with HPLC, GC and GC-MS analyses during 1 year of silicification of the hyperthermophilic Archaea Methanocaldococcus jannaschii. Although the cells themselves did not fossilise and the accompanying extracellular polymeric substances (EPS) did, our analyses showed that the OM initially present in both cells and EPS was uniformly preserved in the precipitated silica, with amino acids and fatty acids being the best preserved compounds. This study thus completes previous data obtained by electron microscopy investigations of simulated microbial fossilisation and can help better identification and interpretation of microbial biosignatures in both ancient rocks and in recent hydrothermal formations and sediments.

A similarity distance of diversity measure for discriminating mesophilic and thermophilic proteins

The successful prediction of thermophilic proteins is useful for designing stable enzymes that are functional at high temperature. We have used the increment of diversity (ID), a novel amino acid composition-based similarity distance, in a 2-class K-nearest neighbor classifier to classify thermophilic and mesophilic proteins. And the KNN-ID classifier was successfully developed to predict the thermophilic proteins. Instead of extracting features from protein sequences as done previously, our approach was based on a diversity measure of symbol sequences. The similarity distance between each pair of protein sequences was first calculated to quantitatively measure the similarity level of one given sequence and the other. The query protein is then determined using the K-nearest neighbor algorithm. Comparisons with multiple recently published methods showed that the KNN-ID proposed in this study outperforms the other methods. The improved predictive performance indicated it is a simple and effective classifier for discriminating thermophilic and mesophilic proteins. At last, the influence of protein length and protein identity on prediction accuracy was discussed further. The prediction model and dataset used in this article can be freely downloaded from http://wlxy.imu.edu.cn/college/biostation/fuwu/KNN-ID/index.htm .

Ca(2+)-binding reduces conformational flexibility of RC-LH1 core complex from thermophile Thermochromatium tepidum

The light-harvesting complex, LH1, of thermophile purple bacteria Thermochromatium tepidum consists of an array of α- and β-polypeptides which assemble the photoactive bacteriochlorophyll and closely interact with the membrane-lipids. In this study, we investigated the effect of calcium and manganese ions on the protein structure and thermostability of the reaction centre (RC)-LH1/lipid complex. The binding of Ca(2+), but not Mn(2+) is shown to shift the LH1 Q ( y ) absorption maximum from ~889 to 915 nm and to significantly raise the thermostability of the RC-LH1 complex. The ATR-FTIR spectra indicate that interaction of Ca(2+) as monitored by the carboxylates' vibration of aspartate residues, but not Mn(2+) induces changes in the α-helix packing arrangement. The reduced rate of (1)H/(2)H exchange of proteins' amide protons shows that the accessibility to (2)H(2)O is significantly lowered in Ca(2+)-substituted RC-LH1/lipid complexes. In particular, exchange with the associated lipid molecules, is significantly retarded. These results suggest that the thermostability of the RC-LH1 complex is raised by the distinct interaction with calcium cations which reduces the RC-LH1/lipid dynamics, particularly, at the membrane-water interface.

Kinetic and thermodynamic study of cloned thermostable endo-1,4-β-xylanase from Thermotoga petrophila in mesophilic host

The 1,044 bp endo-1,4-β-xylanase gene of a hyperthermophilic Eubacterium, "Thermotoga petrophila RKU 1" (T. petrophila) was amplified, from the genomic DNA of donor bacterium, cloned and expressed in mesophilic host E. coli strain BL21 Codon plus. The extracellular target protein was purified by heat treatment followed by anion and cation exchange column chromatography. The purified enzyme appeared as a single band, corresponding to molecular mass of 40 kDa, upon SDS-PAGE. The pH and temperature profile showed that enzyme was maximally active at 6.0 and 95 °C, respectively against birchwood xylan as a substrate (2,600 U/mg). The enzyme also exhibited marked activity towards beech wood xylan (1,655 U/mg). However minor activity against CMC (61 U/mg) and β-Glucan barley (21 U/mg) was observed. No activity against Avicel, Starch, Laminarin and Whatman filter paper 42 was observed. The K(m), V(max) and K (cat) of the recombinant enzyme were found to be 3.5 mg ml(-1), 2778 μmol mg(-1)min(-1) and 2,137,346.15 s(-1), respectively against birchwood xylan as a substrate. The recombinant enzyme was found very stable and exhibited half life (t(½)) of 54.5 min even at temperature as high as 96 °C, with enthalpy of denaturation (ΔH*(D)), free energy of denaturation (ΔG*(D)) and entropy of denaturation (ΔS*(D)) of 513.23 kJ mol(-1), 104.42 kJ mol(-1) and 1.10 kJ mol(-1)K(-1), respectively at 96 °C. Further the enthalpy (ΔH*), Gibbs free energy (ΔG*) and entropy (ΔS*) for birchwood xylan hydrolysis by recombinant endo-1,4-β-xylanase were calculated at 95 °C as 62.45 kJ mol(-1), 46.18 kJ mol(-1) and 44.2 J mol(-1) K(-1), respectively.