Thermal adaptation analyzed by comparison of protein sequences from mesophilic and extremely thermophilic Methanococcus species

The genome sequence of the extremely thermophilic archaeon Methanococcus jannaschii provides a wealth of data on proteins from a thermophile. In this paper, sequences of 115 proteins from M. jannaschii are compared with their homologs from mesophilic Methanococcus species. Although the growth temperatures of the mesophiles are about 50°C below that of M. jannaschii, their genomic G+C contents are nearly identical. The properties most correlated with the proteins of the thermophile include higher residue volume, higher residue hydrophobicity, more charged amino acids (especially Glu, Arg, and Lys), and fewer uncharged polar residues (Ser, Thr, Asn, and Gln). These are recurring themes, with all trends applying to 83–92% of the proteins for which complete sequences were available. Nearly all of the amino acid replacements most significantly correlated with the temperature change are the same relatively conservative changes observed in all proteins, but in the case of the mesophile/thermophile comparison there is a directional bias. We identify 26 specific pairs of amino acids with a statistically significant (P < 0.01) preferred direction of replacement.

Structural distribution of stability in a thermophilic enzyme

Stability parameters for individual residues in Thermus thermophilus cysteine-free RNase H were determined by native state hydrogen exchange, thus providing a unique comparison of regional thermodynamics between thermophilic and mesophilic homologues. The general distribution of stability in the thermophilic protein is similar to that of its mesophilic homologue, with a proportional increase in stability for almost all residues. As a consequence, the residue-specific stabilities of the two proteins are remarkably similar under conditions where their global stabilities are the same. These results indicate that T. thermophilus RNase H is stabilized in a delocalized fashion, preserving a finely tuned balance of stabilizing interactions throughout the structure. Therefore, although protein stability can be altered by single amino acid substitution, evolution for optimal function may require more subtle and delocalized mechanisms.

Identification of thermophilic species by the amino acid compositions deduced from their genomes

The global amino acid compositions as deduced from the complete genomic sequences of six thermophilic archaea, two thermophilic bacteria, 17 mesophilic bacteria and two eukaryotic species were analysed by hierarchical clustering and principal components analysis. Both methods showed an influence of several factors on amino acid composition. Although GC content has a dominant effect, thermophilic species can be identified by their global amino acid compositions alone. This study presents a careful statistical analysis of factors that affect amino acid composition and also yielded specific features of the average amino acid composition of thermophilic species. Moreover, we introduce the first example of a ‘compositional tree’ of species that takes into account not only homologous proteins, but also proteins unique to particular species. We expect this simple yet novel approach to be a useful additional tool for the study of phylogeny at the genome level.

The thermodynamic origin of the stability of a thermophilic ribozyme

Understanding the mechanism of thermodynamic stability of an RNA structure has significant implications for the function and design of RNA. We investigated the equilibrium folding of a thermophilic ribozyme and its mesophilic homologue by using hydroxyl radical protection, small-angle x-ray scattering, and circular dichroism. Both RNAs require Mg2+ to fold to their native structures that are very similar. The stability is measured as a function of Mg2+ and urea concentrations at different temperatures. The enhanced stability of the thermophilic ribozyme primarily is derived from a tremendous increase in the amount of structure formed in the ultimate folding transition. This increase in structure formation and cooperativity arises because the penultimate and the ultimate folding transitions in the mesophilic ribozyme become linked into a single transition in the folding of the thermophilic ribozyme. Therefore, the starting point, or reference state, for the transition to the native, functional thermophilic ribozyme is significantly less structured. The shift in the reference state, and the resulting increase in folding cooperativity, is likely due to the stabilization of selected native interactions that only form in the ultimate transition. This mechanism of using a less structured intermediate and increased cooperativity to achieve higher functional stability for tertiary RNAs is fundamentally different from that commonly proposed to explain the increased stability of thermophilic proteins.

Crystal structure of a thermostable type B DNA polymerase from Thermococcus gorgonarius

Most known archaeal DNA polymerases belong to the type B family, which also includes the DNA replication polymerases of eukaryotes, but maintain high fidelity at extreme conditions. We describe here the 2.5 Å resolution crystal structure of a DNA polymerase from the Archaea Thermococcus gorgonarius and identify structural features of the fold and the active site that are likely responsible for its thermostable function. Comparison with the mesophilic B type DNA polymerase gp43 of the bacteriophage RB69 highlights thermophilic adaptations, which include the presence of two disulfide bonds and an enhanced electrostatic complementarity at the DNA–protein interface. In contrast to gp43, several loops in the exonuclease and thumb domains are more closely packed; this apparently blocks primer binding to the exonuclease active site. A physiological role of this “closed” conformation is unknown but may represent a polymerase mode, in contrast to an editing mode with an open exonuclease site. This archaeal B DNA polymerase structure provides a starting point for structure-based design of polymerases or ligands with applications in biotechnology and the development of antiviral or anticancer agents.

Adjustment of conformational flexibility is a key event in the thermal adaptation of proteins

3-Isopropylmalate dehydrogenase (IPMDH, E.C. 1.1.1.85) from the thermophilic bacterium Thermus thermophilus HB8 is homologous to IPMDH from the mesophilic Escherichia coli, but has an approximately 17°C higher melting temperature. Its temperature optimum is 22–25°C higher than that of the E. coli enzyme; however, it is hardly active at room temperature. The increased conformational rigidity required to stabilize the thermophilic enzyme against heat denaturation might explain its different temperature-activity profile. Hydrogen/deuterium exchange studies were performed on this thermophilic-mesophilic enzyme pair to compare their conformational flexibilities. It was found that Th. thermophilus IPMDH is significantly more rigid at room temperature than E. coli IPMDH, whereas the enzymes have nearly identical flexibilities under their respective optimal working conditions, suggesting that evolutionary adaptation tends to maintain a “corresponding state” regarding conformational flexibility. These observations confirm that conformational fluctuations necessary for catalytic function are restricted at room temperature in the thermophilic enzyme, suggesting a close relationship between conformational flexibility and enzyme function.

A structural view of evolutionary divergence

Two directed evolution experiments on p-nitrobenzyl esterase yielded one enzyme with a 100-fold increased activity in aqueous-organic solvents and another with a 17°C increase in thermostability. Structures of the wild type and its organophilic and thermophilic counterparts are presented at resolutions of 1.5 Å, 1.6 Å, and 2.0 Å, respectively. These structures identify groups of interacting mutations and demonstrate how directed evolution can traverse complex fitness landscapes. Early-generation mutations stabilize flexible loops not visible in the wild-type structure and set the stage for further beneficial mutations in later generations. The mutations exert their influence on the esterase structure over large distances, in a manner that would be difficult to predict. The loops with the largest structural changes generally are not the sites of mutations. Similarly, none of the seven amino acid substitutions in the organophile are in the active site, even though the enzyme experiences significant changes in the organization of this site. In addition to reduction of surface loop flexibility, thermostability in the evolved esterase results from altered core packing, helix stabilization, and the acquisition of surface salt bridges, in agreement with other comparative studies of mesophilic and thermophilic enzymes. Crystallographic analysis of the wild type and its evolved counterparts reveals networks of mutations that collectively reorganize the active site. Interestingly, the changes that led to diversity within the α/β hydrolase enzyme family and the reorganization seen in this study result from main-chain movements.

Carbonic anhydrase is an ancient enzyme widespread in prokaryotes

Carbonic anhydrases catalyze the reversible hydration of CO2 and are ubiquitous in highly evolved eukaryotes. The recent identification of a third class of carbonic anhydrase (γ class) in a methanoarchaeon and our present finding that the β class also extends into thermophilic species from the Archaea domain led us to initiate a systematic search for these enzymes in metabolically and phylogenetically diverse prokaryotes. Here we show that carbonic anhydrase is widespread in the Archaea and Bacteria domains, and is an ancient enzyme. The occurrence in chemolithoautotrophic species occupying deep branches of the universal phylogenetic tree suggests a role for this enzyme in the proposed autotrophic origin of life. The presence of the β and γ classes in metabolically diverse species spanning the Archaea and Bacteria domains demonstrates that carbonic anhydrases have a far more extensive and fundamental role in prokaryotic biology than previously recognized.

A 55-kDa protein isolated from human cells shows DNA glycosylase activity toward 3,N4-ethenocytosine and the G/T mismatch

Etheno adducts in DNA arise from multiple endogenous and exogenous sources. Of these adducts we have reported that, 1,N6-ethenoadenine (ɛA) and 3,N4-ethenocytosine (ɛC) are removed from DNA by two separate DNA glycosylases. We later confirmed these results by using a gene knockout mouse lacking alkylpurine-DNA-N-glycosylase, which excises ɛA. The present work is directed toward identifying and purifying the human glycosylase activity releasing ɛC. HeLa cells were subjected to multiple steps of column chromatography, including two ɛC-DNA affinity columns, which resulted in >1,000-fold purification. Isolation and renaturation of the protein from SDS/polyacrylamide gel showed that the ɛC activity resides in a 55-kDa polypeptide. This apparent molecular mass is approximately the same as reported for the human G/T mismatch thymine-DNA glycosylase. This latter activity copurified to the final column step and was present in the isolated protein band having ɛC-DNA glycosylase activity. In addition, oligonucleotides containing ɛC⋅G or G/T(U), could compete for ɛC protein binding, further indicating that the ɛC-DNA glycosylase is specific for both types of substrates in recognition. The same substrate specificity for ɛC also was observed in a recombinant G/T mismatch DNA glycosylase from the thermophilic bacterium, Methanobacterium thermoautotrophicum THF.

Conversion of a c type cytochrome to a b type that spontaneously forms in vitro from apo protein and heme: Implications for c type cytochrome biogenesis and folding

Cytochrome c552 from Hydrogenobacter thermophilus, a thermophilic bacterium, has been converted into a b type cytochrome, after mutagenesis of both heme-binding cysteines to alanine and expression in the cytoplasm of Escherichia coli. The b type variant is less stable, with the guanidine hydrochloride unfolding midpoint occurring at a concentration 2 M lower than for the wild-type protein. The reduction potential is 75 mV lower than that of the recombinant wild-type protein. The heme can be removed from the b type variant, thus generating an apo protein that has, according to circular dichroism spectroscopy, an α-helical content different from that of the holo b type protein. The latter is readily reformed in vitro by addition of heme to the apo protein. This reforming suggests that previously observed assembly of cytochrome c552, which has the typical class I cytochrome c fold, in the E. coli cytoplasm is a consequence of spontaneous thioether bond formation after binding of heme to a prefolded polypeptide. These observations have implications for the general problem of c type cytochrome biogenesis.

A double-hexamer archaeal minichromosome maintenance protein is an ATP-dependent DNA helicase

The minichromosome maintenance (MCM) proteins are essential for DNA replication in eukaryotes. Thus far, all eukaryotes have been shown to contain six highly related MCMs that apparently function together in DNA replication. Sequencing of the entire genome of the thermophilic archaeon Methanobacterium thermoautotrophicum has allowed us to identify only a single MCM-like gene (ORF Mt1770). This gene is most similar to MCM4 in eukaryotic cells. Here we have expressed and purified the M. thermoautotrophicum MCM protein. The purified protein forms a complex that has a molecular mass of ≈850 kDa, consistent with formation of a double hexamer. The protein has an ATP-independent DNA-binding activity, a DNA-stimulated ATPase activity that discriminates between single- and double-stranded DNA, and a strand-displacement (helicase) activity that can unwind up to 500 base pairs. The 3′ to 5′ helicase activity requires both ATP hydrolysis and a functional nucleotide-binding site. Moreover, the double hexamer form is the active helicase. It is therefore likely that an MCM complex acts as the replicative DNA helicase in eukaryotes and archaea. The simplified replication machinery in archaea may provide a simplified model for assembly of the machinery required for initiation of eukaryotic DNA replication.

Temperature, template topology, and factor requirements of archaeal transcription

Although Archaea are prokaryotic and resemble Bacteria morphologically, their transcription apparatus is remarkably similar to those of eukaryotic cell nuclei. Because some Archaea exist in environments with temperatures of around 100°C, they are likely to have evolved unique strategies for transcriptional control. Here, we investigate the effects of temperature and DNA template topology in a thermophilic archaeal transcription system. Significantly, and in marked contrast with characterized eucaryal systems, archaeal DNA template topology has negligible effect on transcription levels at physiological temperatures using highly purified polymerase and recombinant transcription factors. Furthermore, archaeal transcription does not require hydrolysis of the β-γ phosphoanhydride bond of ATP. However, at lower temperatures, negatively supercoiled templates are transcribed more highly than those that are positively supercoiled. Notably, the block to transcription on positively supercoiled templates at lowered temperatures is at the level of polymerase binding and promoter opening. These data imply that Archaea do not possess a functional homologue of transcription factor TFIIH, and that for the promoters studied, transcription is mediated by TATA box-binding protein, transcription factor TFB, and RNA polymerase alone. Furthermore, they suggest that the reduction of plasmid linking number by hyperthermophilic Archaea in vivo in response to cold shock is a mechanism to maintain gene expression under these adverse circumstances.

B12-dependent ribonucleotide reductases from deeply rooted eubacteria are structurally related to the aerobic enzyme from Escherichia coli

The ribonucleotide reductases from three ancient eubacteria, the hyperthermophilic Thermotoga maritima (TM), the radioresistant Deinococcus radiodurans (DR), and the thermophilic photosynthetic Chloroflexus aurantiacus, were found to be coenzyme-B12 (class II) enzymes, similar to the earlier described reductases from the archaebacteria Thermoplasma acidophila and Pyrococcus furiosus. Reduction of CDP by the purified TM and DR enzymes requires adenosylcobalamin and DTT. dATP is a positive allosteric effector, but stimulation of the TM enzyme only occurs close to the temperature optimum of 80–90°C. The TM and DR genes were cloned by PCR from peptide sequence information. The TM gene was sequenced completely and expressed in Escherichia coli. The deduced amino acid sequences of the two eubacterial enzymes are homologous to those of the archaebacteria. They can also be aligned to the sequence of the large protein of the aerobic E. coli ribonucleotide reductase that belongs to a different class (class I), which is not dependent on B12. Structure determinations of the E. coli reductase complexed with substrate and allosteric effectors earlier demonstrated a 10-stranded β/α-barrel in the active site. From the conservation of substrate- and effector-binding residues we propose that the B12-dependent class II enzymes contain a similar barrel.

Thermostable and site-specific DNA binding of the gene product ORF56 from the Sulfolobus islandicus plasmid pRN1, a putative archael plasmid copy control protein

There is still a lack of information on the specific characteristics of DNA-binding proteins from hyperthermophiles. Here we report on the product of the gene orf56 from plasmid pRN1 of the acidophilic and thermophilic archaeon Sulfolobus islandicus. orf56 has not been characterised yet but low sequence similarily to several eubacterial plasmid-encoded genes suggests that this 6.5 kDa protein is a sequence-specific DNA-binding protein. The DNA-binding properties of ORF56, expressed in Escherichia coli, have been investigated by EMSA experiments and by fluorescence anisotropy measurements. Recombinant ORF56 binds to double-stranded DNA, specifically to an inverted repeat located within the promoter of orf56. Binding to this site could down-regulate transcription of the orf56 gene and also of the overlapping orf904 gene, encoding the putative initiator protein of plasmid replication. By gel filtration and chemical crosslinking we have shown that ORF56 is a dimeric protein. Stoichiometric fluorescence anisotropy titrations further indicate that ORF56 binds as a tetramer to the inverted repeat of its target binding site. CD spectroscopy points to a significant increase in ordered secondary structure of ORF56 upon binding DNA. ORF56 binds without apparent cooperativity to its target DNA with a dissociation constant in the nanomolar range. Quantitative analysis of binding isotherms performed at various salt concentrations and at different temperatures indicates that approximately seven ions are released upon complex formation and that complex formation is accompanied by a change in heat capacity of –6.2 kJ/mol. Furthermore, recombinant ORF56 proved to be highly thermostable and is able to bind DNA up to 85°C.

Characterization of RNase P from Thermotoga maritima

The protein subunit of RNase P from a thermophilic bacterium, Thermotoga maritima, was overexpressed in and purified from Escherichia coli. The cloned protein was reconstituted with the RNA subunit transcribed in vitro. The temperature optimum of the holoenzyme is near 50°C, with no enzymatic activity at 65°C or above. This finding is in sharp contrast to the optimal growth temperature of T.maritima, which is near 80°C. However, in heterologous reconstitution experiments in vitro with RNase P subunits from other species, we found that the protein subunit from T.maritima was responsible for the comparative thermal stability of such complexes.