New understandings of thermostable and peizostable enzymes

Biochemical characterization of a xylose-tolerant GH43 β-xylosidase from Geobacillus thermodenitrificans

Hemicelluloses are the second most abundant polysaccharide in plant biomass, in which xylan is the main constituent. Aiming at the total degradation of xylan and the obtention of fermentable sugars, several enzymes acting synergistically are required, especially β-xylosidases. In this study, β-xylosidase from Geobacillus thermodenitrificans (GtXyl) was expressed in E. coli BL21 and characterized. The enzyme GtXyl has been grouped within the family of glycoside hydrolases 43 (GH43). Results showed that GtXyl obtained the highest activity at pH 5.0 and temperature of 60 °C. In the additive's tests, the enzyme remained stable in the presence of metal ions and EDTA, and showed high tolerance to xylose, with a relative activity of 55.4% at 400 mM. The enzyme also presented bifunctional activity of β-xylosidase and α-l-arabinofuranosidase, with the highest activity on the substrate p-nitrophenyl-β-d-xylopyranoside. The specific activity on p-nitrophenyl-β-d-xylopyranoside was 18.33 U mg-1 and catalytic efficiency of 20.21 mM-1 s-1, which is comparable to other β-xylosidases reported in the literature. Putting together, the GtXyl enzyme presented interesting biochemical characteristics that are desirable for the application in the enzymatic hydrolysis of plant biomass, such as activity at higher temperatures, high thermostability and stability to metal ions.

Thermotolerant and protease-resistant GH5 family β-mannanase with CBM1 from Penicillium aculeatum APS1: purification and characterization

In past several years, mannanases has attracted many researchers owing to its extensive industrial applications. The search for novel mannanases with high stability still continues. Present investigation was focused on purification of extracellular β-mannanase from Penicillium aculeatum APS1 and its characterization. APS1 mannanase was purified to homogeneity by chromatography techniques. Protein identification by MALDI-TOF MS/MS revealed that the enzyme belongs to GH family 5 and subfamily 7, and possesses CBM1. The molecular weight was found to be 40.6 kDa. The optimum temperature and pH of APS1 mannanase were 70 °C and 5.5, respectively. APS1 mannanase was found to be highly stable at 50 °C and tolerant at 55-60 °C. The enzyme was very sensitive to Mn+2, Hg+2 and Co+2 metal ions and stimulated by Zn+2. Inhibition of activity by N-bromosuccinimide suggested key role of tryptophan residues for catalytic activity. The purified enzyme was efficient in hydrolysis of locust bean gum, guar gum and konjac gum and kinetic studies revealed highest affinity towards locust bean gum (LBG). APS1 mannanase was found to be protease resistant. Looking at the properties, APS1 mannanase can be a valuable candidate for applications in bioconversion of mannan-rich substrates into value-added products and also in food and feed processing.

Less Unfavorable Salt Bridges on the Enzyme Surface Result in More Organic Cosolvent Resistance

Biocatalysis for the synthesis of fine chemicals is highly attractive but usually requires organic (co-)solvents (OSs). However, native enzymes often have low activity and resistance in OSs and at elevated temperatures. Herein, we report a smart salt bridge design strategy for simultaneously improving OS resistance and thermostability of the model enzyme, Bacillus subtilits Lipase A (BSLA). We combined comprehensive experimental studies of 3450 BSLA variants and molecular dynamics simulations of 36 systems. Iterative recombination of four beneficial substitutions yielded superior resistant variants with up to 7.6-fold (D64K/D144K) improved resistance toward three OSs while exhibiting significant thermostability (thermal resistance up to 137-fold, and half-life up to 3.3-fold). Molecular dynamics simulations revealed that locally refined flexibility and strengthened hydration jointly govern the highly increased resistance in OSs and at 50-100 °C. The salt bridge redesign provides protein engineers with a powerful and likely general approach to design OSs- and/or thermal-resistant lipases and other α/β-hydrolases.

Efficient biomass saccharification using a novel cellobiohydrolase from Clostridium clariflavum for utilization in biofuel industry

The present study describes the cloning of the cellobiohydrolase gene from a thermophilic bacterium Clostridium clariflavum and its expression in Escherichia coli BL21(DE3) utilizing the expression vector pET-21a(+). The optimization of various parameters (pH, temperature, isopropyl β-d-1-thiogalactopyranoside (IPTG) concentration, time of induction) was carried out to obtain the maximum enzyme activity (2.78 ± 0.145 U ml⁻¹) of recombinant enzyme. The maximum expression of recombinant cellobiohydrolase was obtained at pH 6.0 and 70 °C respectively. Enzyme purification was performed by heat treatment and immobilized metal anionic chromatography. The specific activity of the purified enzyme was 57.4 U mg⁻¹ with 35.17% recovery and 3.90 purification fold. Sodium dodecyl sulphate polyacrylamide gel electrophoresis (SDS-PAGE) showed that the molecular weight of cellobiohydrolase was 78 kDa. Among metal ions, Ca²⁺ showed a positive impact on the cellobiohydrolase enzyme with increased activity by 115%. Recombinant purified cellobiohydrolase enzyme remained stable and exhibited 77% and 63% residual activity in comparison to control in the presence of n-butanol and after incubation at 80 °C for 1 h, respectively. Our results indicate that our purified recombinant cellobiohydrolase can be used in the biofuel industry.

Successful expression of a novel cellobiohydrolase enzyme from Clostridium clariflavum with efficient saccharification potential of plant biomass for the biofuel industry.

Coevolution of the Activity and Thermostability of an ϵ-Keto Ester Reductase for Better Synthesis of an (R)-α-Lipoic Acid Precursor

In this work, we have identified a significantly improved variant (S131Y/Q252I) of the natural ϵ-keto ester reductase CpAR2 from Candida parapsilosis for efficiently manufacturing (R)-8-chloro-6-hydroxyoctanoic acid [(R)-ECHO] through co-evolution of activity and thermostability. The activity of the variant CpAR2_S131Y/Q252I towards the ϵ-keto ester ethyl 8-chloro-6-oxooctanoate was improved to 214 U mg^-1 -from 120 U mg^-1 in the case of the wild-type enzyme (CpAR2_WT )-and the half-deactivating temperature (T₅₀ , for 15 min incubation) was simultaneously increased by 2.3 °C in relation to that of CpAR2_WT . Consequently, only 2 g L^-1 of lyophilized E. coli cells harboring CpAR2_S131Y/Q252I and a glucose dehydrogenase (GDH) were required in order to achieve productivity similar to that obtained in our previous work, under optimized reaction conditions (530 g L^-1  d^-1 ). This result demonstrated a more economical and efficient process for the production of the key (R)-α-lipoic acid intermediate ethyl 8-chloro-6-oxooctanoate.

Filling the Void: Introducing Aromatic Interactions into Solvent Tunnels To Enhance Lipase Stability in Methanol

An enhanced stability of enzymes in organic solvents is desirable under industrial conditions. The potential of lipases as biocatalysts is mainly limited by their denaturation in polar alcohols. In this study, we focused on selected solvent tunnels in lipase from Geobacillus stearothermophilus T6 to improve its stability in methanol during biodiesel synthesis. Using rational mutagenesis, bulky aromatic residues were incorporated to occupy solvent channels and induce aromatic interactions leading to a better inner core packing. The chemical and structural characteristics of each solvent tunnel were systematically analyzed. Selected residues were replaced with Phe, Tyr, or Trp. Overall, 16 mutants were generated and screened in 60% methanol, from which 3 variants showed an enhanced stability up to 81-fold compared with that of the wild type. All stabilizing mutations were found in the longest tunnel detected in the "closed-lid" X-ray structure. The combination of Phe substitutions in an A187F/L360F double mutant resulted in an increase in unfolding temperature (T_m ) of 7°C in methanol and a 3-fold increase in biodiesel synthesis yield from waste chicken oil. A kinetic analysis with p-nitrophenyl laurate revealed that all mutants displayed lower hydrolysis rates (k _cat), though their stability properties mostly determined the transesterification capability. Seven crystal structures of different variants were solved, disclosing new π-π or CH/π intramolecular interactions and emphasizing the significance of aromatic interactions for improved solvent stability. This rational approach could be implemented for the stabilization of other enzymes in organic solvents.IMPORTANCE Enzymatic synthesis in organic solvents holds increasing industrial opportunities in many fields; however, one major obstacle is the limited stability of biocatalysts in such a denaturing environment. Aromatic interactions play a major role in protein folding and stability, and we were inspired by this to redesign enzyme voids. The rational protein engineering of solvent tunnels of lipase from Geobacillus stearothermophilus is presented here, offering a promising approach to introduce new aromatic interactions within the enzyme core. We discovered that longer tunnels leading from the surface to the enzyme active site were more beneficial targets for mutagenesis for improving lipase stability in methanol during biodiesel biosynthesis. A structural analysis of the variants confirmed the generation of new interactions involving aromatic residues. This work provides insights into stability-driven enzyme design by targeting the solvent channel void.

Structural signatures of thermal adaptation of bacterial ribosomal RNA, transfer RNA, and messenger RNA

Temperature adaptation of bacterial RNAs is a subject of both fundamental and practical interest because it will allow a better understanding of molecular mechanism of RNA folding with potential industrial application of functional thermophilic or psychrophilic RNAs. Here, we performed a comprehensive study of rRNA, tRNA, and mRNA of more than 200 bacterial species with optimal growth temperatures (OGT) ranging from 4°C to 95°C. We investigated temperature adaptation at primary, secondary and tertiary structure levels. We showed that unlike mRNA, tRNA and rRNA were optimized for their structures at compositional levels with significant tertiary structural features even for their corresponding randomly permutated sequences. tRNA and rRNA are more exposed to solvent but remain structured for hyperthermophiles with nearly OGT-independent fluctuation of solvent accessible surface area within a single RNA chain. mRNA in hyperthermophiles is essentially the same as random sequences without tertiary structures although many mRNA in mesophiles and psychrophiles have well-defined tertiary structures based on their low overall solvent exposure with clear separation of deeply buried from partly exposed bases as in tRNA and rRNA. These results provide new insight into temperature adaptation of different RNAs.

Pressure adaptation of 3-isopropylmalate dehydrogenase from an extremely piezophilic bacterium is attributed to a single amino acid substitution

3-Isopropylmalate dehydrogenase (IPMDH) from the extreme piezophile Shewanella benthica (SbIPMDH) is more pressure-tolerant than that from the atmospheric pressure-adapted Shewanella oneidensis (SoIPMDH). To understand the molecular mechanisms of this pressure tolerance, we analyzed mutated enzymes. The results indicate that only a single mutation at position 266, corresponding to Ala (SbIPMDH) and Ser (SoIPMDH), essentially affects activity under higher-pressure conditions. Structural analyses of SoIPMDH suggests that penetration of three water molecules into the cleft around Ser266 under high-pressure conditions could reduce the activity of the wild-type enzyme; however, no water molecule is observed in the Ala266 mutant.

Discovery of an intermolecular disulfide bond required for the thermostability of a heterodimeric protein from the thermophile Hydrogenobacter thermophilus

Factors that increase protein thermostability are of considerable interest in both scientific and industrial fields. Disulfide bonds are one of such factors that increase thermostability, but are rarely found in intracellular proteins because of the reducing environment of the cytosol. Here, we report the first example of an intermolecular disulfide bond between heteromeric subunits of a novel-type phosphoserine phosphatase from a thermophilic bacterium Hydrogenobacter thermophilus, which contributes to the protein thermostability at the physiological temperature. Comparison of remaining soluble proteins between wild-type and cysteine-deleted mutant using SDS-PAGE revealed that the disulfide bond increases the thermostability of the whole protein by tightly connecting a subunit with low solubility to the partner with higher solubility. Furthermore, it was strongly suggested that the disulfide bond is formed and contributes to the stability in vivo. This finding will open new avenues for the design of proteins with increased thermostability.

Archaeal genome guardians give insights into eukaryotic DNA replication and damage response proteins

As the third domain of life, archaea, like the eukarya and bacteria, must have robust DNA replication and repair complexes to ensure genome fidelity. Archaea moreover display a breadth of unique habitats and characteristics, and structural biologists increasingly appreciate these features. As archaea include extremophiles that can withstand diverse environmental stresses, they provide fundamental systems for understanding enzymes and pathways critical to genome integrity and stress responses. Such archaeal extremophiles provide critical data on the periodic table for life as well as on the biochemical, geochemical, and physical limitations to adaptive strategies allowing organisms to thrive under environmental stress relevant to determining the boundaries for life as we know it. Specifically, archaeal enzyme structures have informed the architecture and mechanisms of key DNA repair proteins and complexes. With added abilities to temperature-trap flexible complexes and reveal core domains of transient and dynamic complexes, these structures provide insights into mechanisms of maintaining genome integrity despite extreme environmental stress. The DNA damage response protein structures noted in this review therefore inform the basis for genome integrity in the face of environmental stress, with implications for all domains of life as well as for biomanufacturing, astrobiology, and medicine.

Engineering a serum-resistant and thermostable vesicular stomatitis virus G glycoprotein for pseudotyping retroviral and lentiviral vectors

Vesicular stomatitis virus G glycoprotein (VSV-G) is the most widely used envelope protein for retroviral and lentiviral vector pseudotyping; however, serum inactivation of VSV-G pseudotyped vectors is a significant challenge for in vivo gene delivery. To address this problem, we conducted directed evolution of VSV-G to increase its resistance to human serum neutralization. After six selection cycles, numerous common mutations were present. Based on their location within VSV-G, we analyzed whether substitutions in several surface exposed residues could endow viral vectors with higher resistance to serum. S162T, T230N, and T368A mutations enhanced serum resistance, and additionally K66T, T368A, and E380K substitutions increased the thermostability of VSV-G pseudotyped retroviral vectors, an advantageous byproduct of the selection strategy. Analysis of a number of combined mutants revealed that VSV-G harboring T230N + T368A or K66T + S162T + T230N + T368A mutations exhibited both higher in vitro resistance to human serum and higher thermostability, as well as enhanced resistance to rabbit and mouse serum. Finally, lentiviral vectors pseudotyped with these variants were more resistant to human serum in a murine model. These serum-resistant and thermostable VSV-G variants may aid the application of retroviral and lentiviral vectors to gene therapy.

Terminal amino acids disturb xylanase thermostability and activity

Protein structure is composed of regular secondary structural elements (α-helix and β-strand) and non-regular region. Unlike the helix and strand, the non-regular region consists of an amino acid defined as a disordered residue (DR). When compared with the effect of the helix and strand, the effect of the DR on enzyme structure and function is elusive. An Aspergillus niger GH10 xylanase (Xyn) was selected as a model molecule of (β/α)(8) because the general structure consists of ~10% enzymes. The Xyn has five N-terminal DRs and one C-terminal DR, respectively, which were deleted to construct three mutants, XynΔN, XynΔC, and XynΔNC. Each mutant was ~2-, 3-, or 4-fold more thermostable and 7-, 4-, or 4-fold more active than the Xyn. The N-terminal deletion decreased the xylanase temperature optimum for activity (T(opt)) 6 °C, but the C-terminal deletion increased its T(opt) 6 °C. The N- and C-terminal deletions had opposing effects on the enzyme T(opt) but had additive effects on its thermostability. The five N-terminal DR deletions had more effect on the enzyme kinetics but less effect on its thermo property than the one C-terminal DR deletion. CD data showed that the terminal DR deletions increased regular secondary structural contents, and hence, led to slow decreased Gibbs free energy changes (ΔG(0)) in the thermal denaturation process, which ultimately enhanced enzyme thermostabilities.

Thermal stability of cytochrome c₅ of pressure-sensitive Shewanella livingstonensis

Cytochrome c₅ of pressure-sensitive Shewanella livingstonensis (SL cytc₅) exhibits lower thermal stability than a highly homologous counterpart of pressure-tolerant Shewanella violacea. This stability difference is due to an enthalpic effect that can be attributed to the amino acid residue at position 50 (Leu or Lys). These cytc₅ proteins are appropriate materials for understanding the protein stability mechanism.

In silico characterization of thermostable lipases

Thermostable lipases are of high priority for industrial applications as they are endowed with the capability of carrying out diversified reactions at elevated temperatures. Extremophiles are their potential source. Sequence and structure annotation of thermostable lipases can elucidate evolution of lipases from their mesophilic counterparts with enhanced thermostability hence better industrial potential. Sequence analysis highlighted the conserved residues in bacterial and fungal thermostable lipases. Higher frequency of AXXXA motif and poly Ala residues in lid domain of thermostable Bacillus lipases were distinguishing characteristics. Comparison of amino acid composition among thermostable and mesostable lipases brought into light the role of neutral, charged and aromatic amino acid residues in enhancement of thermostability. Structural annotation of thermostable lipases with that of mesostable lipases revealed some striking features which are increment of gamma turns in thermostable lipases; being first time reported in our paper, longer beta strands, lesser beta-branched residues in helices, increase in charged-neutral hydrogen bonding pair, hydrophobic-hydrophobic contact and differences in the N-cap and C-cap residues of the α helices. Conclusively, it can be stated that subtle changes in the arrangement of amino acid residues in the tertiary structure of lipases contributes to enhanced thermostability.

Temperature dependence of the flexibility of thermophilic and mesophilic flavoenzymes of the nitroreductase fold

A widely held hypothesis regarding the thermostability of thermophilic proteins states asserts that, at any given temperature, thermophilic proteins are more rigid than their mesophilic counterparts. Many experimental and computational studies have addressed this question with conflicting results. Here, we compare two homologous enzymes, one mesophilic (Escherichia coli FMN-dependent nitroreductase; NTR) and one thermophilic (Thermus thermophilus NADH oxidase; NOX), by multiple molecular dynamics simulations at temperatures from 5 to 100 degrees C. We find that the global rigidity/flexibility of the two proteins, assessed by a variety of metrics, is similar on the time scale of our simulations. However, the thermophilic enzyme retains its native conformation to a much greater degree at high temperature than does the mesophilic enzyme, both globally and within the active site. The simulations identify the helix F-helix G 'arm' as the region with the greatest difference in loss of native contacts between the two proteins with increasing temperature. In particular, a network of electrostatic interactions holds helix F to the body of the protein in the thermophilic protein, and this network is absent in the mesophilic counterpart.

Inhibition of hemorragic snake venom components: old and new approaches

Snake venoms are complex toxin mixtures. Viperidae and Crotalidae venoms, which are hemotoxic, are responsible for most of the envenomations around the world. Administration of antivenins aimed at the neutralization of toxins in humans is prone to potential risks. Neutralization of snake venom toxins has been achieved through different approaches: plant extracts have been utilized in etnomedicine. Direct electric current from low voltage showed neutralizing properties against venom phospholipase A2 and metalloproteases. This mini-review summarizes new achievements in venom key component inhibition. A deeper knowledge of alternative ways to inhibit venom toxins may provide supplemental treatments to serum therapy.

Shape and evolution of thermostable protein structure

Organisms evolved at high temperatures must maintain their proteins' structures in the face of increased thermal disorder. This challenge results in differences in residue utilization and overall structure. Focusing on thermostable/mesostable pairs of homologous structures, we have examined these differences using novel geometric measures: specifically burial depth (distance from the molecular surface to each atom) and travel depth (distance from the convex hull to the molecular surface that avoids the protein interior). These along with common metrics like packing and Wadell Sphericity are used to gain insight into the constraints experienced by thermophiles. Mean travel depth of hyperthermostable proteins is significantly less than that of their mesostable counterparts, indicating smaller, less numerous and less deep pockets. The mean burial depth of hyperthermostable proteins is significantly higher than that of mesostable proteins indicating that they bury more atoms further from the surface. The burial depth can also be tracked on the individual residue level, adding a finer level of detail to the standard exposed surface area analysis. Hyperthermostable proteins for the first time are shown to be more spherical than their mesostable homologues, regardless of when and how they adapted to extreme temperature. Additionally, residue specific burial depth examinations reveal that charged residues stay unburied, most other residues are slightly more buried and Alanine is more significantly buried in hyperthermostable proteins.

Structure of the alkalohyperthermophilic Archaeoglobus fulgidus lipase contains a unique C-terminal domain essential for long-chain substrate binding

Several crystal structures of AFL, a novel lipase from the archaeon Archaeoglobus fulgidus, complexed with various ligands, have been determined at about 1.8 A resolution. This enzyme has optimal activity in the temperature range of 70-90 degrees C and pH 10-11. AFL consists of an N-terminal alpha/beta-hydrolase fold domain, a small lid domain, and a C-terminal beta-barrel domain. The N-terminal catalytic domain consists of a 6-stranded beta-sheet flanked by seven alpha-helices, four on one side and three on the other side. The C-terminal lipid binding domain consists of a beta-sheet of 14 strands and a substrate covering motif on top of the highly hydrophobic substrate binding site. The catalytic triad residues (Ser136, Asp163, and His210) and the residues forming the oxyanion hole (Leu31 and Met137) are in positions similar to those of other lipases. Long-chain lipid is located across the two domains in the AFL-substrate complex. Structural comparison of the catalytic domain of AFL with a homologous lipase from Bacillus subtilis reveals an opposite substrate binding orientation in the two enzymes. AFL has a higher preference toward long-chain substrates whose binding site is provided by a hydrophobic tunnel in the C-terminal domain. The unusually large interacting surface area between the two domains may contribute to thermostability of the enzyme. Two amino acids, Asp61 and Lys101, are identified as hinge residues regulating movement of the lid domain. The hydrogen-bonding pattern associated with these two residues is pH dependent, which may account for the optimal enzyme activity at high pH. Further engineering of this novel lipase with high temperature and alkaline stability will find its use in industrial applications.

Superoxide dismutase from the eukaryotic thermophile Alvinella pompejana: structures, stability, mechanism, and insights into amyotrophic lateral sclerosis

Prokaryotic thermophiles supply stable human protein homologs for structural biology; yet, eukaryotic thermophiles would provide more similar macromolecules plus those missing in microbes. Alvinella pompejana is a deep-sea hydrothermal-vent worm that has been found in temperatures averaging as high as 68 degrees C, with spikes up to 84 degrees C. Here, we used Cu,Zn superoxide dismutase (SOD) to test if this eukaryotic thermophile can provide insights into macromolecular mechanisms and stability by supplying better stable mammalian homologs for structural biology and other biophysical characterizations than those from prokaryotic thermophiles. Identification, cloning, characterization, X-ray scattering (small-angle X-ray scattering, SAXS), and crystal structure determinations show that A. pompejana SOD (ApSOD) is superstable, homologous, and informative. SAXS solution analyses identify the human-like ApSOD dimer. The crystal structure shows the active site at 0.99 A resolution plus anchoring interaction motifs in loops and termini accounting for enhanced stability of ApSOD versus human SOD. Such stabilizing features may reduce movements that promote inappropriate intermolecular interactions, such as amyloid-like filaments found in SOD mutants causing the neurodegenerative disease familial amyotrophic lateral sclerosis or Lou Gehrig's disease. ApSOD further provides the structure of a long-sought SOD product complex at 1.35 A resolution, suggesting a unified inner-sphere mechanism for catalysis involving metal ion movement. Notably, this proposed mechanism resolves apparent paradoxes regarding electron transfer. These results extend knowledge of SOD stability and catalysis and suggest that the eukaryote A. pompejana provides macromolecules highly similar to those from humans, but with enhanced stability more suitable for scientific and medical applications.

A stabilizing alpha/beta-hydrophobic core greatly contributes to hyperthermostability of archaeal [P62A]Ssh10b

The hyperthermophilic Ssh10b from Sulfolobus shibatae is a member of the Sac10b family, which has been postulated to play a role in chromosomal organization in Archaea. Ssh10b is capable of significantly constraining negative DNA supercoils at elevated temperatures. In this study, the solution structure of the dimeric P62A mutant Ssh10b ([P62A]Ssh10b) was determined by multidimensional NMR spectroscopy. The backbone 15N dynamics, H/D exchange with and without the denaturant GdmSCN, and chemical and thermal denaturation experiments were performed to investigate the molecular basis of high thermostability of [P62A]Ssh10b. Data analysis has revealed an alpha/beta-hydrophobic core consisting of two alpha-helices and one beta-sheet which are stabilized by cooperative hydrophobic and hydrogen-bonding interactions. This stabilizing alpha/beta-hydrophobic core of [P62A]Ssh10b exhibiting highly restricted internal motions is composed of residues having highly protected amide protons which exchange with solvent mostly by means of a global unfolding process. The K40N mutation greatly destabilizes the mutant [P62A]Ssh10b because this mutation disturbs the packing of alpha-helix against the beta-sheet reducing the stability of the alpha/beta-hydrophobic core in the mutant protein. In comparison with homologous mesophilic and thermophilic proteins, it can be presumed that the stabilizing alpha/beta-hydrophobic core in the [P62A]Ssh10b structure greatly contributes to the high thermostability of the protein.

Discrimination of mesophilic and thermophilic proteins using machine learning algorithms

Discriminating thermophilic proteins from their mesophilic counterparts is a challenging task and it would help to design stable proteins. In this work, we have systematically analyzed the amino acid compositions of 3075 mesophilic and 1609 thermophilic proteins belonging to 9 and 15 families, respectively. We found that the charged residues Lys, Arg, and Glu as well as the hydrophobic residues, Val and Ile have higher occurrence in thermophiles than mesophiles. Further, we have analyzed the performance of different methods, based on Bayes rules, logistic functions, neural networks, support vector machines, decision trees and so forth for discriminating mesophilic and thermophilic proteins. We found that most of the machine learning techniques discriminate these classes of proteins with similar accuracy. The neural network-based method could discriminate the thermophiles from mesophiles at the five-fold cross-validation accuracy of 89% in a dataset of 4684 proteins. Moreover, this method is tested with 325 mesophiles in Xylella fastidosa and 382 thermophiles in Aquifex aeolicus and it could successfully discriminate them with the accuracy of 91%. These accuracy levels are better than other methods in the literature and we suggest that this method could be effectively used to discriminate mesophilic and thermophilic proteins.

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

Background

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

Results

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

Conclusion

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

Enhancing the thermostability of a novel beta-agarase AgaB through directed evolution

To increase the thermostability of beta-agarase AgaB by directed evolution, the mutant gene libraries were generated by error-prone polymerase chain reaction (PCR) and deoxyribonucleic acid (DNA) shuffling. Mutants with high thermostability were screened by a simple method based on agarase-degrading agar to generate a clear zone on the agar plate. A mutant S2 was obtained through two rounds of error-prone PCR and a single round of DNA shuffling and selection. It has higher thermostability and slightly increased catalytic activity than wild-type AgaB. Melting temperature (T(m)) of S2, as determined by circular dichroism, is 4.6 degrees C higher than that of wild-type AgaB, and the half-life of S2 is 350 min at 40 degrees C, which is 18.4-fold longer than that of the wild-type enzyme. Saturation mutagenesis and hydrophobic cluster analysis indicated that hydrophobic interaction might be the key factor that enhances the enzyme stability.

Comparative study and mutational analysis of distinctive structural elements of hyperthermophilic enzymes

Comparison of the three-dimensional structure of hyperthermophilic and mesophilic beta-glycosidases shows differences in secondary structure composition. The enzymes from hyperthermophilic archaea have a significantly larger number of beta-strands arranged in supernumerary beta-sheets compared to mesophilic enzymes from bacteria and other organisms. Amino acid replacements designed to alter the structure of the supernumerary beta-strands were introduced by site directed mutagenesis into the sequence encoding the beta-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 beta-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).

Post-duplication charge evolution of phosphoglucose isomerases in teleost fishes through weak selection on many amino acid sites

Background

The partitioning of ancestral functions among duplicated genes by neutral evolution, or subfunctionalization, has been considered the primary process for the evolution of novel proteins (neofunctionalization). Nonetheless, how a subfunctionalized protein can evolve into a more adaptive protein is poorly understood, mainly due to the limitations of current analytical methods, which can detect only strong selection for amino acid substitutions involved in adaptive molecular evolution. In this study, we employed a comparative evolutionary approach to this question, focusing on differences in the structural properties of a protein, specifically the electric charge, encoded by fish-specific duplicated phosphoglucose isomerase (Pgi) genes.

Results

Full-length cDNA cloning, RT-PCR based gene expression analyses, and comparative sequence analyses showed that after subfunctionalization with respect to the expression organ of duplicate Pgi genes, the net electric charge of the PGI-1 protein expressed mainly in internal tissues became more negative, and that of PGI-2 expressed mainly in muscular tissues became more positive. The difference in net protein charge was attributable not to specific amino acid sites but to the sum of various amino acid sites located on the surface of the PGI molecule.

Conclusion

This finding suggests that the surface charge evolution of PGI proteins was not driven by strong selection on individual amino acid sites leading to permanent fixation of a particular residue, but rather was driven by weak selection on a large number of amino acid sites and consequently by steady directional and/or purifying selection on the overall structural properties of the protein, which is derived from many modifiable sites. The mode of molecular evolution presented here may be relevant to various cases of adaptive modification in proteins, such as hydrophobic properties, molecular size, and electric charge.

Microbial extremophiles at the limits of life

Prokaryotic extremophiles were the first representatives of life on Earth and they are responsible for the genesis of geological structures during the evolution and creation of all currently known ecosystems. Flexibility of the genome probably allowed life to adapt to a wide spectrum of extreme environments. As a result, modern prokaryotic diversity formed in a framework of physico-chemical factors, and it is composed of: thermophilic, psychrophilic, acidophilic, alkaliphilic, halophilic, barophilic, and radioresistant species. This artificial systematics cannot reflect the multiple actions of different environmental factors since one organism could unite characteristics of several extreme-groups. In this review we show the current status of studies in all fields of extremophiles and summarize the limits of life for different species of microbial extremophiles. We also discuss the finding of extremophiles from unusual places such as soils, and briefly review recent studies of microfossils in meteorites in the context of the significance of microbial extremophiles to Astrobiology.

Thermostable variants of the recombinant xylanase a from Bacillus subtilis produced by directed evolution show reduced heat capacity changes

Directed evolution techniques have been used to improve the thermal stability of the xylanase A from Bacillus subtilis (XylA). Two generations of random mutant libraries generated by error prone PCR coupled with a single generation of DNA shuffling produced a series of mutant proteins with increasing thermostability. The most Thermostable XylA variant from the third generation contained four mutations Q7H, G13R, S22P, and S179C that showed an increase in melting temperature of 20 degrees C. The thermodynamic properties of a representative subset of nine XylA variants showing a range of thermostabilities were measured by thermal denaturation as monitored by the change in the far ultraviolet circular dichroism signal. Analysis of the data from these thermostable variants demonstrated a correlation between the decrease in the heat capacity change (deltaC(p)) with an increase in the midpoint of the transition temperature (T(m)) on transition from the native to the unfolded state. This result could not be interpreted within the context of the changes in accessible surface area of the protein on transition from the native to unfolded states. Since all the mutations are located at the surface of the protein, these results suggest that an explanation of the decrease in deltaC(p) should include effects arising from the protein/solvent interface.

Stabilisation of a (betaalpha)8-barrel protein designed from identical half barrels

It has been suggested that the common (betaalpha)(8)-barrel enzyme fold has evolved by the duplication and fusion of identical (betaalpha)(4)-half barrels, followed by the optimisation of their interface. In our attempts to reconstruct these events in vitro we have previously linked in tandem two copies of the C-terminal half barrel HisF-C of imidazole glycerol phosphate synthase from Thermotoga maritima and subsequently reconstituted in the fusion construct HisF-CC a salt bridge cluster present in wild-type HisF. The resulting recombinant protein HisF-C*C, which was produced in an insoluble form and unfolded with low cooperativity at moderate urea concentrations has now been stabilised and solubilised by a combination of random mutagenesis and selection in vivo. For this purpose, Escherichia coli cells were transformed with a plasmid-based gene library encoding HisF-C*C variants fused to chloramphenicol acetyltransferase (CAT). Stable and soluble variants were identified by the survival of host cells on solid medium containing high concentrations of the antibiotic. The selected HisF-C*C proteins, which were characterised in vitro in the absence of CAT, contained eight different amino acid substitutions. One of the exchanges (Y143C) stabilised HisF-C*C by the formation of an intermolecular disulfide bond. Three of the substitutions (G245R, V248M, L250Q) were located in the long loop connecting the two HisF-C copies, whose subsequent truncation from 13 to 5 residues yielded the stabilised variant HisF-C*C Delta. From the remaining substitutions, Y143H and V234M were most beneficial, and molecular dynamics simulations suggest that they strengthen the interactions between the half barrels by establishing a hydrogen-bonding network and an extensive hydrophobic cluster, respectively. By combining the loop deletion of HisF-C*C Delta with the Y143H and V234M substitutions, the variant HisF-C**C was generated. Recombinant HisF-C**C is produced in soluble form, forms a pure monomer with its tryptophan residues shielded from solvent and unfolds with similar cooperativity as HisF. Our results show that, starting from two identical and fused half barrels, few amino acid exchanges are sufficient to generate a highly stable and compact (betaalpha)(8)-barrel protein with wild-type like structural properties.

Mechanisms for stabilisation and the maintenance of solubility in proteins from thermophiles

Background

The database of protein structures contains representatives from organisms with a range of growth temperatures. Various properties have been studied in a search for the molecular basis of protein adaptation to higher growth temperature. Charged groups have emerged as key distinguishing factors for proteins from thermophiles and mesophiles.

Results

A dataset of 291 thermophile-derived protein structures is compared with mesophile proteins. Calculations of electrostatic interactions support the importance of charges, but indicate that increases in charge contribution to folded state stabilisation do not generally correlate with the numbers of charged groups. Relative propensities of charged groups vary, such as the substitution of glutamic for aspartic acid sidechains. Calculations suggest an energetic basis, with less dehydration for longer sidechains. Most other properties studied show weak or insignificant separation of proteins from moderate thermophiles or hyperthermophiles and mesophiles, including an estimate of the difference in sidechain rotameric entropy upon protein folding. An exception is increased burial of alanine and proline residues and decreased burial of phenylalanine, methionine, tyrosine and tryptophan in hyperthermophile proteins compared to those from mesophiles.

Conclusion

Since an increase in the number of charged groups for hyperthermophile proteins is separable from charged group contribution to folded state stability, we hypothesise that charged group propensity is important in the context of protein solubility and the prevention of aggregation. Accordingly we find some separation between mesophile and hyperthermophile proteins when looking at the largest surface patch that does not contain a charged sidechain. With regard to our observation that aromatic sidechains are less buried in hyperthermophile proteins, further analysis indicates that the placement of some of these groups may facilitate the reduction of folding fluctuations in proteins of the higher growth temperature organisms.

Potential and utilization of thermophiles and thermostable enzymes in biorefining

In today's world, there is an increasing trend towards the use of renewable, cheap and readily available biomass in the production of a wide variety of fine and bulk chemicals in different biorefineries. Biorefineries utilize the activities of microbial cells and their enzymes to convert biomass into target products. Many of these processes require enzymes which are operationally stable at high temperature thus allowing e.g. easy mixing, better substrate solubility, high mass transfer rate, and lowered risk of contamination. Thermophiles have often been proposed as sources of industrially relevant thermostable enzymes. Here we discuss existing and potential applications of thermophiles and thermostable enzymes with focus on conversion of carbohydrate containing raw materials. Their importance in biorefineries is explained using examples of lignocellulose and starch conversions to desired products. Strategies that enhance thermostablity of enzymes both in vivo and in vitro are also assessed. Moreover, this review deals with efforts made on developing vectors for expressing recombinant enzymes in thermophilic hosts.

Purification and characterization of two sugarcane bagasse-absorbable thermophilic xylanases from the mesophilic Cellulomonas flavigena

We report the purification and characterization of two thermophilic xylanases from the mesophilic bacteria Cellulomonas flavigena grown on sugarcane bagasse (SCB) as the only carbon source. Extracellular xylanase activity produced by C. flavigena was found both free in the culture supernatant and associated with residual SCB. To identify some of the molecules responsible for the xylanase activity in the substrate-bound fraction, residual SCB was treated with 3 M guanidine hydrochloride and then with 6 M urea. Further analysis of the eluted material led to the identification of two xylanases Xyl36 (36 kDa) and Xyl53 (53 kDa). The pI for Xyl36 was 5.0, while the pI for Xyl53 was 4.5. Xyl36 had a Km value of 1.95 mg/ml, while Xyl53 had a Km value of 0.78 mg/ml. In addition to SCB, Xyl36 and Xyl53 were also able to bind to insoluble oat spelt xylan and Avicel, as shown by substrate-binding assays. Xyl36 and Xyl53 showed optimal activity at pH 6.5, and at optimal temperature 65 and 55 degrees C, respectively. Xyl36 and Xyl53 retained 24 and 35%, respectively, of their original activity after 8 h of incubation at their optimal temperature. As far as we know, this is the first study on the thermostability properties of purified xylanases from microorganisms belonging to the genus Cellulomonas.

Directed evolution: an approach to engineer enzymes

Directed evolution is being used increasingly in industrial and academic laboratories to modify and improve commercially important enzymes. Laboratory evolution is thought to make its biggest contribution in explorations of non-natural functions, by allowing us to distinguish the properties nurtured by evolution. In this review we report the significant advances achieved with respect to the methods of biocatalyst improvement and some critical properties and applications of the modified enzymes. The application of directed evolution has been elaborately demonstrated for protein solubility, stability and catalytic efficiency. Modification of certain enzymes for their application in enantioselective catalysis has also been elucidated. By providing a simple and reliable route to enzyme improvement, directed evolution has emerged as a key technology for enzyme engineering and biocatalysis.

Comparison of crystal structure interactions and thermodynamics for stabilizing mutations in the Tetrahymena ribozyme

Although general mechanisms of RNA folding and catalysis have been elucidated, little is known about how ribozymes achieve structural stability at high temperature. A previous in vitro evolution experiment identified a small number of mutations that significantly increase the thermostability of the tertiary structure of the Tetrahymena ribozyme. Because we also determined the crystal structure of this thermostable ribozyme, we have for the first time the opportunity to compare the structural interactions and thermodynamic contributions of individual nucleotides in a ribozyme. We investigated the contribution of five mutations to thermostability by using temperature gradient gel electrophoresis. Unlike the case with several well-studied proteins, the effects of individual mutations on thermostability of this RNA were highly context dependent. The three most important mutations for thermostability were actually destabilizing in the wild-type background. A269G and A304G contributed to stability only when present as a pair, consistent with their proximity in the ribozyme structure. In an evolutionary context, this work supports and extends the idea that one advantage of protein enzyme systems over an RNA world is the ability of proteins to accumulate stabilizing single-site mutations, whereas RNA may often require much rarer double mutations to improve the stability of both its tertiary and secondary structures.

On the role of basic residues in adapting the reaction centre-LH1 complex for growth at elevated temperatures in purple bacteria

The purple photosynthetic bacterium Thermochromatium tepidum is a moderate thermophile, with a growth optimum of 48-50 degrees C. The X-ray crystal structure of the reaction centre from this organism has been determined, and compared with that from mesophilic bacteria such as Blastochloris viridis and Rhodobacter sphaeroides (Nogi T et al. (2000) Proc Natl Acad Sci USA 97: 13561-13566). Structural features that could contribute to the enhanced thermal stability of the Thermochromatium tepidum reaction centre were discussed, including three arginine residues exposed at the periplasmic side of the membrane that are not present in reaction centres from mesophilic organisms, and potentially could increase the affinity of the complex for the surrounding membrane. In the present report these arginine residues, plus a histidine identified from an extensive sequence alignment, were engineered into structurally homologous positions in the Rhodobacter sphaeroides reaction centre, and the effect on the thermal stability of the Rhodobacter sphaeroides complex was examined. We find that these residues do not enhance the thermal stability of the reaction centre, as assessed by absorbance spectroscopy of the bacteriochlorin cofactors in membrane-bound reaction centres. Possible roles of these residues in the Thermochromatium tepidum reaction centre are discussed, and it is proposed that they facilitate stronger binding of the reaction centre to the encircling LH1 antenna complex, through ionic interactions with acidic residues at the C-terminal end of the LH1 alpha-polypeptide. Such an interaction could enhance the stability of the so-called 'RC-LH1 core' complex that is formed between the reaction centre and the LH1 antenna, and which represents the minimal functional photosynthetic unit in all known purple photosynthetic bacteria. Stronger bonding interactions between the two complexes could also contribute to an increase in the rigidity of the photosynthetic membrane in Thermochromatium tepidum, in accord with the general finding that the cytoplasmic membrane from thermophilic eubacteria is less fluid than its counterpart in mesophilic bacteria.

Calcium-induced tertiary structure modifications of endo-beta-1,3-glucanase from Pyrococcus furiosus in 7.9 M guanidinium chloride

The family 16 endo-beta-1,3 glucanase from the extremophilic archaeon Pyrococcus furiosus is a laminarinase, which in 7.9 M GdmCl (guanidinium chloride) maintains a significant amount of tertiary structure without any change of secondary structure. The addition of calcium to the enzyme in 7.9 M GdmCl causes significant changes to the near-UV CD and fluorescence spectra, suggesting a notable increase in the tertiary structure which leads to a state comparable, but not identical, to the native state. The capability to interact with calcium in 7.9 M GdmCl with a consistent recovery of native tertiary structure is a unique property of this extremely stable endo-beta-1,3 glucanase. The effect of calcium on the thermodynamic parameters relative to the GdmCl-induced equilibrium unfolding has been analysed by CD and fluorescence spectroscopy. The interaction of calcium with the native form of the enzyme is studied by Fourier-transform infrared spectroscopy in the absorption region of carboxylate groups and by titration in the presence of a chromophoric chelator. A homology-based model of the enzyme is generated and used to predict the putative binding site(s) for calcium and the structural interactions potentially responsible for the unusual stability of this protein, in comparison with other family 16 glycoside hydrolases.

Linear and non-linear pressure dependence of enzyme catalytic parameters

The pressure dependence of enzyme catalytic parameters allows volume changes associated with substrate binding and activation volumes for the chemical steps to be determined. Because catalytic constants are composite parameters, elementary volume change contributions can be calculated from the pressure differentiation of kinetic constants. Linear and non-linear pressure-dependence of single-step enzyme reactions and steady-state catalytic parameters can be observed. Non-linearity can be interpreted either in terms of interdependence between the pressure and other environmental parameters (i.e., temperature, solvent composition, pH), pressure-induced enzyme unfolding, compressibility changes and pressure-induced rate limiting changes. These different situations are illustrated with several examples.