Protein thermostability in extremophiles

Extremophiles and Extremophilic Behaviour-New Insights and Perspectives

Extremophiles, throughout evolutionary time, have evolved a plethora of unique strategies to overcome hardships associated with the environments they are found in. Modifying their genome, showing a bias towards certain amino acids, redesigning their proteins, and enhancing their membranes and other organelles with specialised chemical compounds are only some of those strategies. Scientists can utilise such attributes of theirs for a plethora of biotechnological and astrobiological applications. Moreover, the rigorous study of such microorganisms regarding their evolution and ecological niche can offer deep insight into science's most paramount inquiries such as how life originated on Earth and whether we are alone in the universe. The intensification of studies involving extremophiles in the future can prove to be highly beneficial for humanity, even potentially ameliorating modern problems such as those related to climate change while also expanding our knowledge about the complex biochemical reactions that ultimately resulted in life as we know it today.

The structure of Mycobacterium thermoresistibile MmpS5 reveals a conserved disulfide bond across mycobacteria

The tuberculosis (TB) emergency has been a pressing health threat for decades. With the emergence of drug-resistant TB and complications from the COVID-19 pandemic, the TB health crisis is more serious than ever. Mycobacterium tuberculosis (Mtb), the causative agent of TB, requires iron for its survival. Thus, Mtb has evolved several mechanisms to acquire iron from the host. Mtb produces two siderophores, mycobactin and carboxymycobactin, which scavenge for host iron. Mtb siderophore-dependent iron acquisition requires the export of apo-siderophores from the cytosol to the host environment and import of iron-bound siderophores. The export of Mtb apo-siderophores across the inner membrane is facilitated by two mycobacterial inner membrane proteins with their cognate periplasmic accessory proteins, designated MmpL4/MmpS4 and MmpL5/MmpS5. Notably, the Mtb MmpL4/MmpS4 and MmpL5/MmpS5 complexes have also been implicated in the efflux of anti-TB drugs. Herein, we solved the crystal structure of M. thermoresistibile MmpS5. The MmpS5 structure reveals a previously uncharacterized, biologically relevant disulfide bond that appears to be conserved across the Mycobacterium MmpS4/S5 homologs, and comparison with structural homologs suggests that MmpS5 may be dimeric.

Structural Determination and Chemical Synthesis of the N-Glycan from the Hyperthermophilic Archaeon Thermococcus kodakarensis

Asparagine-linked protein glycosylations (N-glycosylations) are one of the most abundant post-translational modifications and are essential for various biological phenomena. Herein, we describe the isolation, structural determination, and chemical synthesis of the N-glycan from the hyperthermophilic archaeon Thermococcus kodakarensis. The N-glycan from the organism possesses a unique structure including myo-inositol, which has not been found in previously characterized N-glycans. In this structure, myo-inositol is highly glycosylated and linked with a disaccharide unit through a phosphodiester. The straightforward synthesis of this glycan was accomplished through diastereoselective phosphorylation and phosphodiester construction by S_N 2 coupling. Considering the early divergence of hyperthermophilic organisms in evolution, this study can be expected to open the door to approaching the primitive function of glycan modification at the molecular level.

Charge reversal mutations in mesophilic-thermophilic orthologous protein pairs and their role in enhancing coulombic interaction energy

Proteins from thermophilic organisms are a matter of immense interest for decades because of its application in fields like de-novo protein design, thermostable variants of biocatalysts etc. Previous studies have found several sequence and structural adaptations related to thermal stability, while charge reversal study remains ignored. Here we address whether charge reversal mutations naturally occur in mesophilic-thermophilic/hyperthermophilic orthologous proteins. Do they contribute to thermal stability? Our systematic study on 1550 mesophilic-thermophilic/hyperthermophilic orthologous protein pairs with remarkable structural and topological similarity, shows gain in coulombic interaction energy in thermophilic/hyperthermophilic proteins at short range associated with partially exposed and buried charge reversal mutations, which may enhance thermostability. Our findings call forth its application in future protein engineering studies. Communicated by Ramaswamy H. Sarma.

The Thermostability of Nucleoside Phosphorylases from Prokaryotes. I. The Role of the Primary Structure of the N-terminal fragment of the Protein in the Thermostability of Uridine Phosphorylases

Mutant uridine phosphorylase genes from Shewanella oneidensis MR-1 (S. oneidensis) were constructed by site-directed mutagenesis and strains-producers of the corresponding recombinant (F5I and F5G) proteins were obtained on the basis of Escherichia coli cells. The mutant proteins were purified and their physicochemical and enzymatic properties were studied. It was shown that the N-terminal fragment of uridine phosphorylase plays an important role in the thermal stabilization of the enzyme as a whole. The role of the aminoacid (a.a.) residue phenylalanine (F5) in the formation of thermotolerance of uridine phosphorylases from gamma-proteobacteria was revealed.

Molecular Characterization of a DNA Polymerase from Thermus thermophilus MAT72 Phage vB_Tt72: A Novel Type-A Family Enzyme with Strong Proofreading Activity

We present a structural and functional analysis of the DNA polymerase of thermophilic Thermus thermophilus MAT72 phage vB_Tt72. The enzyme shows low sequence identity (<30%) to the members of the type-A family of DNA polymerases, except for two yet uncharacterized DNA polymerases of T. thermophilus phages: φYS40 (91%) and φTMA (90%). The Tt72 polA gene does not complement the Escherichia colipolA− mutant in replicating polA-dependent plasmid replicons. It encodes a 703-aa protein with a predicted molecular weight of 80,490 and an isoelectric point of 5.49. The enzyme contains a nucleotidyltransferase domain and a 3′-5′ exonuclease domain that is engaged in proofreading. Recombinant enzyme with His-tag at the N-terminus was overproduced in E. coli, subsequently purified by immobilized metal affinity chromatography, and biochemically characterized. The enzyme exists in solution in monomeric form and shows optimum activity at pH 8.5, 25 mM KCl, and 0.5 mM Mg2+. Site-directed analysis proved that highly-conserved residues D15, E17, D78, D180, and D184 in 3′-5′ exonuclease and D384 and D615 in the nucleotidyltransferase domain are critical for the enzyme’s activity. Despite the source of origin, the Tt72 DNA polymerase has not proven to be highly thermoresistant, with a temperature optimum at 55 °C. Above 60 °C, the rapid loss of function follows with no activity > 75 °C. However, during heat treatment (10 min at 75 °C), trehalose, trimethylamine N-oxide, and betaine protected the enzyme against thermal inactivation. A midpoint of thermal denaturation at Tm = 74.6 °C (ΔHcal = 2.05 × 104 cal mol−1) and circular dichroism spectra > 60 °C indicate the enzyme’s moderate thermal stability.

Structural and thermodynamic analysis of factors governing the stability and thermal folding/unfolding of SazCA

Molecular basis of protein stability at different temperatures is a fundamental problem in protein science that is substantially far from being accurately and quantitatively solved as it requires an explicit knowledge of the temperature dependence of folding free energy of amino acid residues. In the present study, we attempted to gain insights into the thermodynamic stability of SazCA and its implications on protein folding/unfolding. We report molecular dynamics simulations of water solvated SazCA in a temperature range of 293-393 K to study the relationship between the thermostability and flexibility. Our structural analysis shows that the protein maintains the highest structural stability at 353 K and the protein conformations are highly flexible at temperatures above 353 K. Larger exposure of hydrophobic surface residues to the solvent medium for conformations beyond 353 K were identified from H-bond analysis. Higher number of secondary structure contents exhibited by SazCA at 353 K corroborated the conformations at 353 K to exhibit the highest thermal stability. The analysis of thermodynamics of protein stability revealed that the conformations that denature at higher melting temperatures tend to have greater maximum thermal stability. Our analysis shows that 353 K conformations have the highest melting temperature, which was found to be close to the experimental optimum temperature. The enhanced protein stability at 353 K due the least value of heat capacity at unfolding suggested an increase in folding. Comparative Gibbs free energy analysis and funnel shaped energy landscape confirmed a transition in folding/unfolding pathway of SazCA at 353 K.

Molecular dynamics simulations identify the regions of compromised thermostability in SazCA

The present study examined the structure and dynamics of the most active and thermostable carbonic anhydrase, SazCA, probed using molecular dynamics simulations. The molecular system was described by widely used biological force-fields (AMBER, CHARMM22, CHARMM36, and OPLS-AA) in conjunction with TIP3P water model. The comparison of molecular dynamics simulation results suggested AMBER to be a suitable choice to describe the structure and dynamics of SazCA. In addition to this, we also addressed the effect of temperature on the stability of SazCA. We performed molecular dynamics simulations at 313, 333, 353, 373, and 393 K to study the relationship between thermostability and flexibility in SazCA. The amino acid residues VAL98, ASN99, GLY100, LYS101, GLU145, and HIS207 were identified as the most flexible residues from root-mean-square fluctuations. The salt bridge analysis showed that ion-pairs ASP113-LYS81, ASP115-LYS81, ASP115-LYS114, GLU144-LYS143, and GLU144-LYS206, were responsible for the compromised thermal stability of SazCA.

Heat-resistant cytosolic malate dehydrogenases (cMDHs) of thermophilic intertidal snails (genus Echinolittorina): protein underpinnings of tolerance to body temperatures reaching 55°C

Snails of the genus Echinolittorina are among the most heat-tolerant animals; they experience average body temperatures near 41-44°C in summer and withstand temperatures up to at least 55°C. Here, we demonstrate that heat stability of function (indexed by the Michaelis-Menten constant of the cofactor NADH, K_M^NADH) and structure (indexed by rate of denaturation) of cytosolic malate dehydrogenases (cMDHs) of two congeners (E. malaccana and E. radiata) exceeds values previously found for orthologs of this protein from less thermophilic species. The ortholog of E. malaccana is more heat stable than that of E. radiata, in keeping with the congeners' thermal environments. Only two inter-congener differences in amino acid sequence in these 332 residue proteins were identified. In both cases (positions 48 and 114), a glycine in the E. malaccana ortholog is replaced by a serine in the E. radiata protein. To explore the relationship between structure and function and to characterize how amino acid substitutions alter stability of different regions of the enzyme, we used molecular dynamics simulation methods. These computational methods allow determination of thermal effects on fine-scale movements of protein components, for example, by estimating the root mean square deviation in atom position over time and the root mean square fluctuation for individual residues. The minor changes in amino acid sequence favor temperature-adaptive change in flexibility of regions in and around the active sites. Interspecific differences in effects of temperature on fine-scale protein movements are consistent with the differences in thermal effects on binding and rates of heat denaturation.

Structural and Dynamics Comparison of Thermostability in Ancient, Modern, and Consensus Elongation Factor Tus

Rationally engineering thermostability in proteins would create enzymes and receptors that function under harsh industrial applications. Several sequence-based approaches can generate thermostable variants of mesophilic proteins. To gain insight into the mechanisms by which proteins become more stable, we use structural and dynamic analyses to compare two popular approaches, ancestral sequence reconstruction (ASR) and the consensus method, used to generate thermostable variants of Elongation Factor Thermo-unstable (EF-Tu). We present crystal structures of ancestral and consensus EF-Tus, accompanied by molecular dynamics simulations aimed at probing the strategies employed to enhance thermostability. All proteins adopt crystal structures similar to extant EF-Tus, revealing no difference in average structure between the methods. Molecular dynamics reveals that ASR-generated sequences retain dynamic properties similar to extant, thermostable EF-Tu from Thermus aquaticus, while consensus EF-Tu dynamics differ from evolution-based sequences. This work highlights the advantage of ASR for engineering thermostability while preserving natural motions in multidomain proteins.

Exploring Covalent Allosteric Inhibition of Antigen 85C from Mycobacterium tuberculosis by Ebselen Derivatives

Previous studies identified ebselen as a potent in vitro and in vivo inhibitor of the Mycobacterium tuberculosis (Mtb) antigen 85 (Ag85) complex, comprising three homologous enzymes required for the biosynthesis of the mycobacterial cell wall. In this study, the Mtb Ag85C enzyme was cocrystallized with azido and adamantyl ebselen derivatives, resulting in two crystallographic structures of 2.01 and 1.30 Å resolution, respectively. Both structures displayed the anticipated covalent modification of the solvent accessible, noncatalytic Cys209 residue forming a selenenylsulfide bond. Continuous difference density for both thiol modifiers allowed for the assessment of interactions that influence ebselen binding and inhibitor orientation that were unobserved in previous Ag85C ebselen structures. The k_inact/K_I values for ebselen, adamantyl ebselen, and azido ebselen support the importance of observed constructive chemical interactions with Arg239 for increased in vitro efficacy toward Ag85C. To better understand the in vitro kinetic properties of these ebselen derivatives, the energetics of specific protein-inhibitor interactions and relative reaction free energies were calculated for ebselen and both derivatives using density functional theory. These studies further support the different in vitro properties of ebselen and two select ebselen derivatives from our previously published ebselen library with respect to kinetics and protein-inhibitor interactions. In both structures, the α9 helix was displaced farther from the enzyme active site than the previous Ag85C ebselen structure, resulting in the restructuring of a connecting loop and imparting a conformational change to residues believed to play a role in substrate binding specific to Ag85C. These notable structural changes directly affect protein stability, reducing the overall melting temperature by up to 14.5 °C, resulting in the unfolding of protein at physiological temperatures. Additionally, this structural rearrangement due to covalent allosteric modification creates a sizable solvent network that encompasses the active site and extends to the modified Cys209 residue. In all, this study outlines factors that influence enzyme inhibition by ebselen and its derivatives while further highlighting the effects of the covalent modification of Cys209 by said inhibitors on the structure and stability of Ag85C. Furthermore, the results suggest a strategy for developing new classes of Ag85 inhibitors with increased specificity and potency.

Differential Effects of Hydrophobic Core Packing Residues for Thermodynamic and Mechanical Stability of a Hyperthermophilic Protein

Proteins from organisms that have adapted to environmental extremes provide attractive systems to explore and determine the origins of protein stability. Improved hydrophobic core packing and decreased loop-length flexibility can increase the thermodynamic stability of proteins from hyperthermophilic organisms. However, their impact on protein mechanical stability is not known. Here, we use protein engineering, biophysical characterization, single-molecule force spectroscopy (SMFS), and molecular dynamics (MD) simulations to measure the effect of altering hydrophobic core packing on the stability of the cold shock protein TmCSP from the hyperthermophilic bacterium Thermotoga maritima. We make two variants of TmCSP in which a mutation is made to reduce the size of aliphatic groups from buried hydrophobic side chains. In the first, a mutation is introduced in a long loop (TmCSP L40A); in the other, the mutation is introduced on the C-terminal β-strand (TmCSP V62A). We use MD simulations to confirm that the mutant TmCSP L40A shows the most significant increase in loop flexibility, and mutant TmCSP V62A shows greater disruption to the core packing. We measure the thermodynamic stability (ΔGD-N) of the mutated proteins and show that there is a more significant reduction for TmCSP L40A (ΔΔG = 63%) than TmCSP V62A (ΔΔG = 47%), as might be expected on the basis of the relative reduction in the size of the side chain. By contrast, SMFS measures the mechanical stability (ΔG*) and shows a greater reduction for TmCSP V62A (ΔΔG* = 8.4%) than TmCSP L40A (ΔΔG* = 2.5%). While the impact on the mechanical stability is subtle, the results demonstrate the power of tuning noncovalent interactions to modulate both the thermodynamic and mechanical stability of a protein. Such understanding and control provide the opportunity to design proteins with optimized thermodynamic and mechanical properties.

An overview of siderophores for iron acquisition in microorganisms living in the extreme

Siderophores are iron-chelating molecules produced by microbes when intracellular iron concentrations are low. Low iron triggers a cascade of gene activation, allowing the cell to survive due to the synthesis of important proteins involved in siderophore synthesis and transport. Generally, siderophores are classified by their functional groups as catecholates, hydroxamates and hydroxycarboxylates. Although other chemical structural modifications and functional groups can be found. The functional groups participate in the iron-chelating process when the ferri-siderophore complex is formed. Classified as acidophiles, alkaliphiles, halophiles, thermophiles, psychrophiles, piezophiles, extremophiles have particular iron requirements depending on the environmental conditions in where they grow. Most of the work done in siderophore production by extremophiles is based in siderophore concentration and/or genomic studies determining the presence of siderophore synthesis and transport genes. Siderophores produced by extremophiles are not well known and more work needs to be done to elucidate chemical structures and their role in microorganism survival and metal cycling in extreme environments.

Extremophilic 50S Ribosomal RNA-Binding Protein L35Ae as a Basis for Engineering of an Alternative Protein Scaffold

Due to their remarkably high structural stability, proteins from extremophiles are particularly useful in numerous biological applications. Their utility as alternative protein scaffolds could be especially valuable in small antibody mimetic engineering. These artificial binding proteins occupy a specific niche between antibodies and low molecular weight substances, paving the way for development of innovative approaches in therapeutics, diagnostics, and reagent use. Here, the 50S ribosomal RNA-binding protein L35Ae from the extremophilic archaea Pyrococcus horikoshii has been probed for its potential to serve as a backbone in alternative scaffold engineering. The recombinant wild type L35Ae has a native-like secondary structure, extreme thermal stability (mid-transition temperature of 90°C) and a moderate resistance to the denaturation by guanidine hydrochloride (half-transition at 2.6 M). Chemical crosslinking and dynamic light scattering data revealed that the wild type L35Ae protein has a propensity for multimerization and aggregation correlating with its non-specific binding to a model cell surface of HEK293 cells, as evidenced by flow cytometry. To suppress these negative features, a 10-amino acid mutant (called L35Ae 10X) was designed, which lacks the interaction with HEK293 cells, is less susceptible to aggregation, and maintains native-like secondary structure and thermal stability. However, L35Ae 10X also shows lowered resistance to guanidine hydrochloride (half-transition at 2.0M) and is more prone to oligomerization. This investigation of an extremophile protein’s scaffolding potential demonstrates that lowered resistance to charged chemical denaturants and increased propensity to multimerization may limit the utility of extremophile proteins as alternative scaffolds.

Protein Thermostability Is Owing to Their Preferences to Non-Polar Smaller Volume Amino Acids, Variations in Residual Physico-Chemical Properties and More Salt-Bridges

Introduction

Protein thermostability is an important field for its evolutionary perspective of mesophilic versus thermophilic relationship and for its industrial/ therapeutic applications.

Methods

Presently, a total 400 (200 thermophilic and 200 mesophilic homologue) proteins were studied utilizing several software/databases to evaluate their amino acid preferences. Randomly selected 50 homologous proteins with available PDB-structure of each group were explored for the understanding of the protein charges, isoelectric-points, hydrophilicity, hydrophobicity, tyrosine phosphorylation and salt-bridge occurrences. These 100 proteins were further probed to generate Ramachandran plot/data for the gross secondary structure prediction in and comparison between the thermophilic and mesophilic proteins.

Results

Present results strongly suggest that nonpolar smaller volume amino acids Ala (χ² = 238.54, p<0.001) and Gly (χ² = 73.35, p<0.001) are highly and Val moderately (χ² = 144.43, p<0.001) occurring in the 85% of thermophilic proteins. Phospho-regulated Tyr and redox-sensitive Cys are also moderately distributed (χ²~20.0, p<0.01) in a larger number of thermophilic proteins. A consistent lower distribution of thermophilicity and discretely higher distribution of hydrophobicity is noticed in a large number of thermophilic versus their mesophilic protein homolog. The mean differences of isoelectric points and charges are found to be significantly less (7.11 vs. 6.39, p<0.05 and 1 vs. -0.6, p<0.01, respectively) in thermophilic proteins compared to their mesophilic counterpart. The possible sites for Tyr phosphorylation are noticed to be 25% higher (p<0.05) in thermophilic proteins. The 60% thermophiles are found with higher number of salt bridges in this study. The average percentage of salt-bridge of thermophiles is found to be higher by 20% than their mesophilic homologue. The GLU-HIS and GLU-LYS salt-bridge dyads are calculated to be significantly higher (p<0.05 and p<0.001, respectively) in thermophilic and GLU-ARG is higher in the mesophilic proteins. The Ramachandran plot/ data suggest a higher abundance of the helix, left-handed helix, sheet, nonplanar peptide and lower occurrence of cis peptide, loop/ turn and outlier in thermophiles. Pearson’s correlation result suggests that the isoelectric points of mesophilic and thermophilic proteins are positively correlated (r = 0.93 and 0.84, respectively; p<0.001) to their corresponding charges. And their hydrophilicity is negatively associated with the corresponding hydrophobicity (r = -0.493, p<0.001 and r = -0.324, p<0.05) suggesting their reciprocal evolvement.

Conclusions

Present results for the first time with this large amount of datasets and multiple contributing factors suggest the greater occurrence of hydrophobicity, salt-bridges and smaller volume nonpolar residues (Gly, Ala and Val) and lesser occurrence of bulky polar residues in the thermophilic proteins. A more stoichiometric relationship amongst these factors minimized the hindrance due to side chain burial and increased compactness and secondary structural stability in thermophilic proteins.

Vibrational entropy differences between mesophile and thermophile proteins and their use in protein engineering

We recently introduced ENCoM, an elastic network atomic contact model, as the first coarse-grained normal mode analysis method that accounts for the nature of amino acids and can predict the effect of mutations on thermostability based on changes vibrational entropy. In this proof-of-concept article, we use pairs of mesophile and thermophile homolog proteins with identical structures to determine if a measure of vibrational entropy based on normal mode analysis can discriminate thermophile from mesophile proteins. We observe that in around 60% of cases, thermophile proteins are more rigid at equivalent temperatures than their mesophile counterpart and this difference can guide the design of proteins to increase their thermostability through series of mutations. We observe that mutations separating thermophile proteins from their mesophile orthologs contribute independently to a decrease in vibrational entropy and discuss the application and implications of this methodology to protein engineering.

Molecular interactions within the halophilic, thermophilic, and mesophilic prokaryotic ribosomal complexes: clues to environmental adaptation

Using the available crystal structures of 50S ribosomal subunits from three prokaryotic species: Escherichia coli (mesophilic), Thermus thermophilus (thermophilic), and Haloarcula marismortui (halophilic), we have analyzed different structural features of ribosomal RNAs (rRNAs), proteins, and of their interfaces. We have correlated these structural features with the environmental adaptation strategies of the corresponding species. While dense intra-rRNA packing is observed in thermophilic, loose intra-rRNA packing is observed in halophilic (both compared to mesophilic). Interestingly, protein-rRNA interfaces of both the extremophiles are densely packed compared to that of the mesophilic. The intersubunit bridge regions are almost devoid of cavities, probably ensuring the proper formation of each bridge (by not allowing any loosely packed region nearby). During rRNA binding, the ribosomal proteins experience some structural transitions. Here, we have analyzed the intrinsically disordered and ordered regions of the ribosomal proteins, which are subjected to such transitions. The intrinsically disordered and disorder-to-order transition sites of the thermophilic and mesophilic ribosomal proteins are simultaneously (i) highly conserved and (ii) slowly evolving compared to rest of the protein structure. Although high conservation is observed at such sites of halophilic ribosomal proteins, but slow rate of evolution is absent. Such differences between thermophilic, mesophilic, and halophilic can be explained from their environmental adaptation strategy. Interestingly, a universal biophysical principle evident by a linear relationship between the free energy of interface formation, interface area, and structural changes of r-proteins during assembly is always maintained, irrespective of the environmental conditions.

Reductive activation and structural rearrangement in superoxide reductase: a combined infrared spectroscopic and computational study

Local and global structural changes that enable reductive activation of superoxide reductase are revealed by a combined approach of infrared difference spectroscopy and computational methods.

Analysis of protein thermostability enhancing factors in industrially important thermus bacteria species

Elucidation of evolutionary factors that enhance protein thermostability is a critical problem and was the focus of this work on Thermus species. Pairs of orthologous sequences of T. scotoductus SA-01 and T. thermophilus HB27, with the largest negative minimum folding energy (MFE) as predicted by the UNAFold algorithm, were statistically analyzed. Favored substitutions of amino acids residues and their properties were determined. Substitutions were analyzed in modeled protein structures to determine their locations and contribution to energy differences using PyMOL and FoldX programs respectively. Dominant trends in amino acid substitutions consistent with differences in thermostability between orthologous sequences were observed. T. thermophilus thermophilic proteins showed an increase in non-polar, tiny, and charged amino acids. An abundance of alanine substituted by serine and threonine, as well as arginine substituted by glutamine and lysine was observed in T. thermophilus HB27. Structural comparison showed that stabilizing mutations occurred on surfaces and loops in protein structures.

A comparison of structural and evolutionary attributes of Escherichia coli and Thermus thermophilus small ribosomal subunits: signatures of thermal adaptation

Here we compare the structural and evolutionary attributes of Thermus thermophilus and Escherichia coli small ribosomal subunits (SSU). Our results indicate that with few exceptions, thermophilic 16S ribosomal RNA (16S rRNA) is densely packed compared to that of mesophilic at most of the analogous spatial regions. In addition, we have located species-specific cavity clusters (SSCCs) in both species. E. coli SSCCs are numerous and larger compared to T. thermophilus SSCCs, which again indicates densely packed thermophilic 16S rRNA. Thermophilic ribosomal proteins (r-proteins) have longer disordered regions than their mesophilic homologs and they experience larger disorder-to-order transitions during SSU-assembly. This is reflected in the predicted higher conformational changes of thermophilic r-proteins compared to their mesophilic homologs during SSU-assembly. This high conformational change of thermophilic r-proteins may help them to associate with the 16S ribosomal RNA with high complementary interfaces, larger interface areas, and denser molecular contacts, compared to those of mesophilic. Thus, thermophilic protein-rRNA interfaces are tightly associated with 16S rRNA than their mesophilic homologs. Densely packed 16S rRNA interior and tight protein-rRNA binding of T. thermophilus (compared to those of E. coli) are likely the signatures of its thermal adaptation. We have found a linear correlation between the free energy of protein-RNA interface formation, interface size, and square of conformational changes, which is followed in both prokaryotic and eukaryotic SSU. Disorder is associated with high protein-RNA interface polarity. We have found an evolutionary tendency to maintain high polarity (thereby disorder) at protein-rRNA interfaces, than that at rest of the protein structures. However, some proteins exhibit exceptions to this general trend.

Enhancing the functional properties of thermophilic enzymes by chemical modification and immobilization

The immobilization of proteins (mostly typically enzymes) onto solid supports is mature technology and has been used successfully to enhance biocatalytic processes in a wide range of industrial applications. However, continued developments in immobilization technology have led to more sophisticated and specialized applications of the process. A combination of targeted chemistries, for both the support and the protein, sometimes in combination with additional chemical and/or genetic engineering, has led to the development of methods for the modification of protein functional properties, for enhancing protein stability and for the recovery of specific proteins from complex mixtures. In particular, the development of effective methods for immobilizing large multi-subunit proteins with multiple covalent linkages (multi-point immobilization) has been effective in stabilizing proteins where subunit dissociation is the initial step in enzyme inactivation. In some instances, multiple benefits are achievable in a single process. Here we comprehensively review the literature pertaining to immobilization and chemical modification of different enzyme classes from thermophiles, with emphasis on the chemistries involved and their implications for modification of the enzyme functional properties. We also highlight the potential for synergies in the combined use of immobilization and other chemical modifications.

Large-scale experimental studies show unexpected amino acid effects on protein expression and solubility in vivo in E. coli

The biochemical and physical factors controlling protein expression level and solubility in vivo remain incompletely characterized. To gain insight into the primary sequence features influencing these outcomes, we performed statistical analyses of results from the high-throughput protein-production pipeline of the Northeast Structural Genomics Consortium. Proteins expressed in E. coli and consistently purified were scored independently for expression and solubility levels. These parameters nonetheless show a very strong positive correlation. We used logistic regressions to determine whether they are systematically influenced by fractional amino acid composition or several bulk sequence parameters including hydrophobicity, sidechain entropy, electrostatic charge, and predicted backbone disorder. Decreasing hydrophobicity correlates with higher expression and solubility levels, but this correlation apparently derives solely from the beneficial effect of three charged amino acids, at least for bacterial proteins. In fact, the three most hydrophobic residues showed very different correlations with solubility level. Leu showed the strongest negative correlation among amino acids, while Ile showed a slightly positive correlation in most data segments. Several other amino acids also had unexpected effects. Notably, Arg correlated with decreased expression and, most surprisingly, solubility of bacterial proteins, an effect only partially attributable to rare codons. However, rare codons did significantly reduce expression despite use of a codon-enhanced strain. Additional analyses suggest that positively but not negatively charged amino acids may reduce translation efficiency in E. coli irrespective of codon usage. While some observed effects may reflect indirect evolutionary correlations, others may reflect basic physicochemical phenomena. We used these results to construct and validate predictors of expression and solubility levels and overall protein usability, and we propose new strategies to be explored for engineering improved protein expression and solubility.

Discrimination of thermophilic and mesophilic proteins

Background

There is a considerable literature on the source of the thermostability of proteins from thermophilic organisms. Understanding the mechanisms for this thermostability would provide insights into proteins generally and permit the design of synthetic hyperstable biocatalysts.

Results

We have systematically tested a large number of sequence and structure derived quantities for their ability to discriminate thermostable proteins from their non-thermostable orthologs using sets of mesophile-thermophile ortholog pairs. Most of the quantities tested correspond to properties previously reported to be associated with thermostability. Many of the structure related properties were derived from the Delaunay tessellation of protein structures.

Conclusions

Carefully selected sequence based indices discriminate better than purely structure based indices. Combined sequence and structure based indices improve performance somewhat further. Based on our analysis, the strongest contributors to thermostability are an increase in ion pairs on the protein surface and a more strongly hydrophobic interior.

Distance-dependent statistical potentials for discriminating thermophilic and mesophilic proteins

Identification of the characteristic structural patterns responsible for protein thermostability is theoretically important and practically useful but largely remains an open problem. These patterns may be revealed through comparative study on thermophilic and mesophilic proteins that have distinct thermostability. In this study, we constructed several distance-dependant potentials from thermophilic and mesophilic proteins. These potentials were then used to evaluate the structural difference between thermophilic and mesophilic proteins. We found that using the subtraction or division of the potentials derived from thermophilic and mesophilic proteins can dramatically increase the discriminatory ability. This approach revealed that the ability to distinct the subtle structural features responsible for protein thermostability may be effectively enhanced through rationally designed comparative study.

Contributions of the C-terminal helix to the structural stability of a hyperthermophilic Fe-superoxide dismutase (TcSOD)

Hyperthermophilic superoxide dismutases (SODs) are of particular interest due to their potential industrial importance and scientific merit in studying the molecular mechanisms of protein folding and stability. Compared to the mesophilic SODs, the hyperthermostable Fe-SODs (TcSOD and ApSOD) have an extended C-terminal helix, which forms an additional ion-pairing network. In this research, the role of the extended C-terminus in the structural stability of TcSOD was studied by investigating the properties of two deletion mutants. The results indicated that the ion-pairing network at the C-terminus had limited contributions to the stability of TcSOD against heat- and GdnHCl-induced inactivation. The intactness of the C-terminal helix had dissimilar impact on the two stages of TcSOD unfolding induced by guanidinium chloride. The mutations slightly decreased the Gibbs free energy of the dissociation of the tetrameric enzymes, while greatly affected the stability of the molten globule-like intermediate. These results suggested that the additional ion-pairing network mainly enhanced the structural stability of TcSOD by stabilizing the monomers.

alpha-Amylase: an ideal representative of thermostable enzymes

The conditions prevailing in the industrial applications in which enzymes are used are rather extreme, especially with respect to temperature and pH. Therefore, there is a continuing demand to improve the stability of enzymes and to meet the requirements set by specific applications. In this respect, thermostable enzymes have been proposed to be industrially relevant. In this review, alpha-amylase, a well-established representative of thermostable enzymes, providing an attractive model for the investigation of the structural basis of thermostability of proteins, has been discussed.

Cloning, expression, purification and crystallization of the Rho transcription termination factor from Thermotoga maritima

Rho is an essential ATP-dependent homohexameric helicase that is found in the vast majority of bacterial species. It is responsible for transcription termination at factor-dependent terminators. Rho binds to a specific region of the newly-synthesised mRNA and translocates along the chain until it reaches and disassembles the transcription complex. Basically, two crystallographic structures of Rho hexamer from Escherichia coli have been reported: an open ring with RNA (or ssDNA) bound to the RNA-binding domain, and a closed ring with the RNA bound to both the RNA-binding domain and the ATP-ase domain. The structure of the protein free from RNA is still unknown, but thermophilic bacteria enable an alternative approach to its characterization as their proteins often crystallize more easily than those of their mesophilic homologs. We report here the heterologous expression in E. coli of full-length Rho from the thermophile Thermotoga maritima, a simple protocol for the purification of its hexameric nucleic acid-free form, and the obtainment of 2.4 A-diffracting crystals.

Enzyme sequence and its relationship to hyperbaric stability of artificial and natural fish lactate dehydrogenases

The cDNAs of lactate dehydrogenase b (LDH-b) from both deep-sea and shallow living fish species, Corphaenoides armatus and Gadus morhua respectively, have been isolated, sequenced and their encoded products overproduced as recombinant enzymes in E. coli. The proteins were characterised in terms of their kinetic and physical properties and their ability to withstand high pressures. Although the two proteins are very similar in terms of their primary structure, only 21 differences at the amino acid level exist between them, the enzyme from the deep-sea species has a significantly increased tolerance to pressure and a higher thermostability. It was possible to investigate whether the changes in the N-terminal or C-terminal regions played a greater role in barophilic adaptation by the construction of two chimeric enzymes by use of a common restriction site within the cDNAs. One of these hybrids was found to have even greater pressure stability than the recombinant enzyme from the deep-living fish species. It was possible to conclude that the major adaptive changes to pressure tolerance must be located in the N-terminal region of the protein. The types of changes that are found and their spatial location within the protein structure are discussed. An analysis of the kinetic parameters of the enzymes suggests that there is clearly a trade off between K_m and k_cat values, which likely reflects the necessity of the deep-sea enzyme to operate at low temperatures.

Protein dynamics and stability: the distribution of atomic fluctuations in thermophilic and mesophilic dihydrofolate reductase derived using elastic incoherent neutron scattering

The temperature dependence of the dynamics of mesophilic and thermophilic dihydrofolate reductase is examined using elastic incoherent neutron scattering. It is demonstrated that the distribution of atomic displacement amplitudes can be derived from the elastic scattering data by assuming a (Weibull) functional form that resembles distributions seen in molecular dynamics simulations. The thermophilic enzyme has a significantly broader distribution than its mesophilic counterpart. Furthermore, although the rate of increase with temperature of the atomic mean-square displacements extracted from the dynamic structure factor is found to be comparable for both enzymes, the amplitudes are found to be slightly larger for the thermophilic enzyme. Therefore, these results imply that the thermophilic enzyme is the more flexible of the two.

Ribulose-1,5-bisphosphate carboxylase/oxygenase from thermophilic cyanobacterium Thermosynechococcus elongatus

Ribulose-1,5-bisphosphate carboxylase/oxygenase (rubisco) can be divided into two branches: the "red-like type" of marine algae and the "green-like type" of cyanobacteria, green algae, and higher plants. We found that the "green-like type" rubisco from the thermophilic cyanobacterium Thermosynechococcus elongatus has an almost 2-fold higher specificity factor compared with rubiscos of mesophilic cyanobacteria, reaching the values of higher plants, and simultaneously revealing an improvement in enzyme thermostability. The difference in the activation energies at the transition stages between the oxygenase and carboxylase reactions for Thermosynechococcus elongatus rubisco is very close to that of Galdieria partita and significantly higher than that of spinach. This is the first characterization of a "green-like type" rubisco from thermophilic organism.

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

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

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

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

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.

Adopting selected hydrogen bonding and ionic interactions from Aspergillus fumigatus phytase structure improves the thermostability of Aspergillus niger PhyA phytase

Although it has been widely used as a feed supplement to reduce manure phosphorus pollution of swine and poultry, Aspergillus niger PhyA phytase is unable to withstand heat inactivation during feed pelleting. Crystal structure comparisons with its close homolog, the thermostable Aspergillus fumigatus phytase (Afp), suggest associations of thermostability with several key residues (E35, S42, R168, and R248) that form a hydrogen bond network in the E35-to-S42 region and ionic interactions between R168 and D161 and between R248 and D244. In this study, loss-of-function mutations (E35A, R168A, and R248A) were introduced singularly or in combination into seven mutants of Afp. All seven mutants displayed decreases in thermostability, with the highest loss (25% [P<0.05]) in the triple mutant (E35A R168A R248A). Subsequently, a set of corresponding substitutions were introduced into nine mutants of PhyA to strengthen the hydrogen bonding and ionic interactions. While four mutants showed improved thermostability, the best response came from the quadruple mutant (A58E P65S Q191R T271R), which retained 20% greater (P<0.05) activity after being heated at 80 degrees C for 10 min and had a 7 degrees C higher melting temperature than that of wild-type PhyA. This study demonstrates the functional importance of the hydrogen bond network and ionic interaction in supporting the high thermostability of Afp and the feasibility of adopting these structural units to improve the thermostability of a homologous PhyA phytase.

Differentiation of mesophilic and thermophilic bacteria with fourier transform infrared spectroscopy

In the present study the characterization and differentiation of mesophilic and thermophilic bacteria were investigated by using Fourier transform infrared (FT-IR) spectroscopy. Our results showed significant differences between the FT-IR spectra of mesophilic and thermophilic bacteria. The protein-to-lipid ratio was significantly higher for thermophiles compared to mesophiles. The absorption intensity of the CH(3) asymmetric stretching vibration was higher in thermophilic bacteria, indicating a change in the composition of the acyl chains. The higher intensity/area observed in the CH(2) symmetric stretching mode at 2857 cm(-1), and the CH(2) bending vibration band at 1452 cm(-1), indicated a higher amount of saturated lipids in thermophilic bacteria. The lipid C=O stretching vibration at 1739 cm(-1), which was observed in the mesophilic group, was not observed clearly in the thermophilic group, indicating a difference in packing that is presumably due to the decreased proportion of unsaturated acyl chains in thermophilic bacteria. In addition, the carbonyl groups become hydrogen bonded and the cellular DNA content was lower in thermophilic bacteria. Moreover, in the 1000-400 cm(-1) frequency region, the spectra of each bacterial species belonging to both the mesophilic and thermophilic bacterial groups, showed characteristic differences that were discriminated via dendrogram using cluster analysis. The current study implies that FT-IR spectroscopy could be successfully applied for the rapid comparison of bacterial groups and species to establish either similarities or discrepancies, as well as to confirm biochemical or physiological characteristics.

Folding and assembly pathways of co-chaperonin proteins 10: Origin of bacterial thermostability

To compare folding/assembly processes of heptameric co-chaperonin proteins 10 (cpn10) from different species and search for the origin of thermostability in hyper-thermostable Aquifex aeolicus cpn10 (Aacpn10), we have studied two bacterial variants-Aacpn10 and Escherichia coli cpn10 (GroES)-and compared the results to data on Homo sapiens cpn10 (hmcpn10). Equilibrium denaturation of GroES by urea, guanidine hydrochloride (GuHCl) and temperature results in coupled heptamer-to-monomer transitions in all cases. This is similar to the behavior of Aacpn10 but differs from hmcpn10 denaturation in urea. Time-resolved experiments reveal that GroES unfolds before heptamer dissociation, whereas refolding/reassembly begins with folding of individual monomers; these assemble in a slower step. The sequential folding/assembly mechanism for GroES is rather similar to that observed for Aacpn10 but contradicts the parallel paths of hmcpn10. We reveal that Aacpn10's stability profile is shifted upwards, broadened, and also moved horizontally to higher temperatures, as compared to that of GroES.

Crystal structure of the bifunctional ATP sulfurylase-APS kinase from the chemolithotrophic thermophile Aquifex aeolicus

The thermophilic chemolithotroph, Aquifex aeolicus, expresses a gene product that exhibits both ATP sulfurylase and adenosine-5'-phosphosulfate (APS) kinase activities. These enzymes are usually segregated on two separate proteins in most bacteria, fungi, and plants. The domain arrangement in the Aquifex enzyme is reminiscent of the fungal ATP sulfurylase, which contains a C-terminal domain that is homologous to APS kinase yet displays no kinase activity. Rather, in the fungal enzyme, the motif serves as a sulfurylase regulatory domain that binds the allosteric effector 3'-phosphoadenosine-5'-phosphosulfate (PAPS), the product of true APS kinase. Therefore, the Aquifex enzyme may represent an ancestral homolog of a primitive bifunctional enzyme, from which the fungal ATP sulfurylase may have evolved. In heterotrophic sulfur-assimilating organisms such as fungi, ATP sulfurylase catalyzes the first committed step in sulfate assimilation to produce APS, which is subsequently metabolized to generate all sulfur-containing biomolecules. In contrast, ATP sulfurylase in sulfur chemolithotrophs catalyzes the reverse reaction to produce ATP and sulfate from APS and pyrophosphate. Here, the 2.3 A resolution X-ray crystal structure of Aquifex ATP sulfurylase-APS kinase bifunctional enzyme is presented. The protein dimerizes through its APS kinase domain and contains ADP bound in all four active sites. Comparison of the Aquifex ATP sulfurylase active site with those from sulfate assimilators reveals similar dispositions of the bound nucleotide and nearby residues. This suggests that minor perturbations are responsible for optimizing the kinetic properties for the physiologically relevant direction. The APS kinase active-site lid adopts two distinct conformations, where one conformation is distorted by crystal contacts. Additionally, a disulfide bond is observed in one ATP-binding P-loop of the APS kinase active site. This linkage accounts for the low kinase activity of the enzyme under oxidizing conditions. The thermal stability of the Aquifex enzyme can be explained by the 43% decreased cavity volume found within the protein core.