The Molecular Determinants of Thermoadaptation: Methanococcales as a Case Study.

Previous reports have shown that environmental temperature impacts proteome evolution in Bacteria and Archaea. However, it is unknown whether thermoadaptation mainly occurs via the sequential accumulation of substitutions, massive horizontal gene transfers, or both. Measuring the real contribution of amino acid substitution to thermoadaptation is challenging, because of confounding environmental and genetic factors (e.g., pH, salinity, genomic G + C content) that also affect proteome evolution. Here, using Methanococcales, a major archaeal lineage, as a study model, we show that optimal growth temperature is the major factor affecting variations in amino acid frequencies of proteomes. By combining phylogenomic and ancestral sequence reconstruction approaches, we disclose a sequential substitutional scheme in which lysine plays a central role by fine tuning the pool of arginine, serine, threonine, glutamine, and asparagine, whose frequencies are strongly correlated with optimal growth temperature. Finally, we show that colonization to new thermal niches is not associated with high amounts of horizontal gene transfers. Altogether, although the acquisition of a few key proteins through horizontal gene transfer may have favored thermoadaptation in Methanococcales, our findings support sequential amino acid substitutions as the main factor driving thermoadaptation.

Complete genome sequence of Methanofervidicoccus sp. A16, a thermophilic methanogen isolated from Mid Cayman Rise hydrothermal vent.

Methanofervidicoccus sp. A16 is a novel thermophilic and obligate hydrogenotrophic methanogen isolated from a deep-sea hydrothermal vent chimney sample at the Mid Cayman spreading center, Caribbean Sea. Here we report the complete genome of strain A16, which has one circular chromosome of 1,485,358 bp with a mean G+C content of 35.01 mol%. The complete genome harbors 1442 predicted protein-encoding genes. Genes involved in hydrogenotrophic methane production and N2 fixation were identified in this genome. This study expands our knowledge of methanogenesis at high temperatures and the involvement of these microorganisms in the carbon and nitrogen cycles of deep-sea hydrothermal environments.

Temporal and spatial impact of Spartina alterniflora invasion on methanogens community in Chongming Island, China.

Methane production by methanogens in wetland is recognized as a significant contributor to global warming. Spartina alterniflora (S. alterniflora), which is an invasion plant in China's wetland, was reported to have enormous effects on methane production. But studies on shifts in the methanogen community in response to S. alterniflora invasion at temporal and spatial scales in the initial invasion years are rare. Sediments derived from the invasive species S. alterniflora and the native species Phragmites australis (P. australis) in pairwise sites and an invasion chronosequence patch (4 years) were analyzed to investigate the abundance and community structure of methanogens using quantitative real-time PCR (qPCR) and Denaturing gradient gel electrophoresis (DGGE) cloning of the methyl-coenzyme M reductase A (mcrA) gene. For the pairwise sites, the abundance of methanogens in S. alterniflora soils was lower than that of P. australis soils. For the chronosequence patch, the abundance and diversity of methanogens was highest in the soil subjected to two years invasion, in which we detected some rare groups including Methanocellales and Methanococcales. These results indicated a priming effect at the initial invasion stages of S. alterniflora for microorganisms in the soil, which was also supported by the diverse root exudates. The shifts of methanogen communities after S. alterniflora invasion were due to changes in pH, salinity and sulfate. The results indicate that root exudates from S. alterniflora have a priming effect on methanogens in the initial years after invasion, and the predominate methylotrophic groups (Methanosarcinales) may adapt to the availability of diverse substrates and reflects the potential for high methane production after invasion by S. alterniflora.

Characterization of the Cation Binding Sites in the NCKX2 Na+/Ca2+-K+ Exchanger.

NCKX1-5 are proteins involved in K+-dependent Na+/Ca2+ exchange in various signal tissues. Here we present a homology model of NCKX2 based on the crystal structure of the NCX_Mj transporter found in Methanoccocus jannaschii. Molecular dynamics simulations were performed on the resultant wild-type NCKX2 model and two mutants (D548N and D575N) loaded with either four Na+ ions or one Ca2+ ion and one K+ ion, in line with the experimentally observed transport stoichiometry. The selectivity of the active site in wild-type NCKX2 for Na+, K+, and Li+ and the electrostatic interactions of the positive Na+ ions in the negatively charged active site of wild-type NCKX2 and the two mutants were evaluated from free energy perturbation calculations. For validation of the homology model, our computational results were compared to available experimental data obtained from numerous prior functional studies. The NCKX2 homology model is in good agreement with the discussed experimental data and provides valuable insights into the structure of the active site, which is lined with acidic and polar residues. The binding of the potassium and calcium ions is accomplished via Asp 575 and 548, respectively. Mutation of these residues to Asn alters the functionality of NCKX2 because of the elimination of the favorable carboxylate-cation interactions. The knowledge obtained from the NCKX2 model can be transferred to other isoforms of the NCKX family: newly discovered pathological mutations in NCKX4 and NCKX5 affect residues that are involved in ion binding and/or transport according to our homology model.

Abundance and potential metabolic activity of methanogens in well-aerated forest and grassland soils of an alpine region.

Although methanogens were recently discovered to occur in aerated soils, alpine regions have not been extensively studied for their presence so far. Here, the abundance of archaea and the methanogenic guilds Methanosarcinales, Methanococcales, Methanobacteriales, Methanomicrobiales and Methanocella spp. was studied at 16 coniferous forest and 14 grassland sites located at the montane and subalpine belts of the Northern Limestone Alps (calcareous) and the Austrian Central Alps (siliceous) using quantitative real-time PCR. Abundance of archaea, methanogens and the methanogenic potentials were significantly higher in grasslands than in forests. Furthermore, methanogenic potentials of calcareous soils were higher due to pH. Methanococcales, Methanomicrobiales and Methanocella spp. were detected in all collected samples, which indicates that they are autochthonous, while Methanobacteriales were absent from 4 out of 16 forest soils. Methanosarcinales were absent from 10 out of 16 forest soils and 2 out of 14 grassland soils. Nevertheless, together with Methanococcales they represented the majority of the 16S rRNA gene copies quantified from the grassland soils. Contrarily, forest soils were clearly dominated by Methanococcales. Our results indicate a higher diversity of methanogens in well-aerated soils than previously believed and that pH mainly influences their abundances and activities.

Long-Lasting Gene Conversion Shapes the Convergent Evolution of the Critical Methanogenesis Genes.

Methanogenesis and its key small-molecule methyltransferase Mtr complex are poorly understood despite their pivotal role in Earth's global carbon cycle. Mtr complex is encoded by a conserved mtrEDCBAFGH operon in most methanogens. Here we report that two discrete lineages, Methanococcales and Methanomicrobiales, have a noncanonical mtr operon carrying two copies of mtrA resulting from an ancient duplication. Compared to mtrA-1, mtrA-2 acquires a distinct transmembrane domain through domain shuffling and gene fusion. However, the nontransmembrane domains (MtrA domain) of mtrA-1 and mtrA-2 are homogenized by gene conversion events lasting throughout the long history of these extant methanogens (over 2410 million years). Furthermore, we identified a possible recruitment of ancient nonmethanogenic methyltransferase genes to establish the methanogenesis pathway. These results not only provide novel evolutionary insight into the methanogenesis pathway and methyltransferase superfamily but also suggest an unanticipated long-lasting effect of gene conversion on gene evolution in a convergent pattern.

Solution nuclear magnetic resonance structure of the GATase subunit and structural basis of the interaction between GATase and ATPPase subunits in a two-subunit-type GMPS from Methanocaldococcus jannaschii.

The solution structure of the monomeric glutamine amidotransferase (GATase) subunit of the Methanocaldococcus janaschii (Mj) guanosine monophosphate synthetase (GMPS) has been determined using high-resolution nuclear magnetic resonance methods. Gel filtration chromatography and ¹5N backbone relaxation studies have shown that the Mj GATase subunit is present in solution as a 21 kDa (188-residue) monomer. The ensemble of 20 lowest-energy structures showed root-mean-square deviations of 0.35 ± 0.06 Å for backbone atoms and 0.8 ± 0.06 Å for all heavy atoms. Furthermore, 99.4% of the backbone dihedral angles are present in the allowed region of the Ramachandran map, indicating the stereochemical quality of the structure. The core of the tertiary structure of the GATase is composed of a seven-stranded mixed ß-sheet that is fenced by five α-helices. The Mj GATase is similar in structure to the Pyrococcus horikoshi (Ph) GATase subunit. Nuclear magnetic resonance (NMR) chemical shift perturbations and changes in line width were monitored to identify residues on GATase that were responsible for interaction with magnesium and the ATPPase subunit, respectively. These interaction studies showed that a common surface exists for the metal ion binding as well as for the protein-protein interaction. The dissociation constant for the GATase-Mg(2+) interaction has been found to be ∼1 mM, which implies that interaction is very weak and falls in the fast chemical exchange regime. The GATase-ATPPase interaction, on the other hand, falls in the intermediate chemical exchange regime on the NMR time scale. The implication of this interaction in terms of the regulation of the GATase activity of holo GMPS is discussed.

Promoter independent abortive transcription assays unravel functional interactions between TFIIB and RNA polymerase.

TFIIB-like general transcription factors are required for transcription initiation by all eukaryotic and archaeal RNA polymerases (RNAPs). TFIIB facilitates both recruitment and post-recruitment steps of initiation; in particular, TFIIB stimulates abortive initiation. X-ray crystallography of TFIIB-RNAP II complexes shows that the TFIIB linker region penetrates the RNAP active center, yet the impact of this arrangement on RNAP activity and underlying mechanisms remains elusive. Promoter-independent abortive initiation assays exploit the intrinsic ability of RNAP enzymes to initiate transcription from nicked DNA templates and record the formation of the first phosphodiester bonds. These assays can be used to measure the effect of transcription factors such as TFIIB and RNAP mutations on abortive transcription.

Directed evolution of a model primordial enzyme provides insights into the development of the genetic code.

The contemporary proteinogenic repertoire contains 20 amino acids with diverse functional groups and side chain geometries. Primordial proteins, in contrast, were presumably constructed from a subset of these building blocks. Subsequent expansion of the proteinogenic alphabet would have enhanced their capabilities, fostering the metabolic prowess and organismal fitness of early living systems. While the addition of amino acids bearing innovative functional groups directly enhances the chemical repertoire of proteomes, the inclusion of chemically redundant monomers is difficult to rationalize. Here, we studied how a simplified chorismate mutase evolves upon expanding its amino acid alphabet from nine to potentially 20 letters. Continuous evolution provided an enhanced enzyme variant that has only two point mutations, both of which extend the alphabet and jointly improve protein stability by >4 kcal/mol and catalytic activity tenfold. The same, seemingly innocuous substitutions (Ile→Thr, Leu→Val) occurred in several independent evolutionary trajectories. The increase in fitness they confer indicates that building blocks with very similar side chain structures are highly beneficial for fine-tuning protein structure and function.

Using E. coli-based cell-free protein synthesis to evaluate the kinetic performance of an orthogonal tRNA and aminoacyl-tRNA synthetase pair.

Even though the orthogonal tRNA and aminoacyl-tRNA synthetase pairs derived from the archaeon Methanocaldococcus jannaschii have been used for many years for site-specific incorporation of non-natural amino acids (nnAAs) in Escherichia coli, their kinetic parameters have not been evaluated. Here we use a cell-free protein synthesis (CFPS) system to control the concentrations of the orthogonal components in order to evaluate their performance while supporting synthesis of modified proteins (i.e. proteins with nnAAs). Titration experiments and estimates of turnover numbers suggest that the orthogonal synthetase is a very slow catalyst when compared to the native E. coli synthetases. The estimated k(cat) for the orthogonal synthetase specific to the nnAA p-propargyloxyphenylalanine (pPaF) is 5.4 × 10(-5) s(-1). Thus, this catalyst may be the limiting factor for nnAA incorporation when using this approach. These titration experiments also resulted in the highest reported cell-free accumulation of two different modified proteins (450 ± 20 µg/ml CAT109pAzF and 428±2µg/ml sfGFP23pPaF) using the standard KC6 cell extract and either the PANOx SP or the inexpensive Glu NMP cell-free recipe.

Insights into dynamics of mobile genetic elements in hyperthermophilic environments from five new Thermococcus plasmids.

Mobilome of hyperthermophilic archaea dwelling in deep-sea hydrothermal vents is poorly characterized. To gain insight into genetic diversity and dynamics of mobile genetic elements in these environments we have sequenced five new plasmids from different Thermococcus strains that have been isolated from geographically remote hydrothermal vents. The plasmids were ascribed to two subfamilies, pTN2-like and pEXT9a-like. Gene content and phylogenetic analyses illuminated a robust connection between pTN2-like plasmids and Pyrococcus abyssi virus 1 (PAV1), with roughly half of the viral genome being composed of genes that have homologues in plasmids. Unexpectedly, pEXT9a-like plasmids were found to be closely related to the previously sequenced plasmid pMETVU01 from Methanocaldococcus vulcanius M7. Our data suggests that the latter observation is most compatible with an unprecedented horizontal transfer of a pEXT9a-like plasmid from Thermococcales to Methanococcales. Gene content analysis revealed that thermococcal plasmids encode Hfq-like proteins and toxin-antitoxin (TA) systems of two different families, VapBC and RelBE. Notably, although abundant in archaeal genomes, to our knowledge, TA and hfq-like genes have not been previously found in archaeal plasmids or viruses. Finally, the plasmids described here might prove to be useful in developing new genetic tools for hyperthermophiles.

Expression, purification and biochemical characterization of Methanocaldococcus jannaschii DNA ligase.

We describe the biochemical characterization of Methanocaldococcus jannaschii (M. jannaschii) DNA ligase and its potential application in single nucleotide polymorphism (SNP) genotyping. The recombinant M. jannaschii DNA ligase is an ATP-dependent ligase. The ligase activity was dependent on metal ions of Mg(2+) and Mn(2+). The optimal concentrations of ATP cofactor and Mg(2+) ion were 0.01-2 and 10 mM, respectively. The optimal pH value for DNA ligation was 8.5. High concentrations of NaCl inhibited DNA ligation. The effects of mismatches on joining short oligonucleotides by M. jannaschii DNA ligase were fully characterized. The mismatches at the first position 5' to the nick inhibited ligation more than those at the first position 3' to the nick. The mismatches at other positions 5' to the nick (3rd to 7th sites) exhibited less inhibition on ligation. However, the introduction of a C/C mismatch at the third position 5' to the nick could completely inhibit the ligation of the terminal-mismatched nick of an oligonucleotide duplex by M. jannaschii DNA ligase. Therefore, introducing an additional mismatch at the third position 5' to the SNP site is a more effective approach in genotyping by M. jannaschii DNA ligase.

Structurally conserved Nop56/58 N-terminal domain facilitates archaeal box C/D ribonucleoprotein-guided methyltransferase activity.

Box C/D RNA-protein complexes (RNPs) guide the 2'-O-methylation of nucleotides in both archaeal and eukaryotic ribosomal RNAs. The archaeal box C/D and C'/D' RNP subcomplexes are each assembled with three sRNP core proteins. The archaeal Nop56/58 core protein mediates crucial protein-protein interactions required for both sRNP assembly and the methyltransferase reaction by bridging the L7Ae and fibrillarin core proteins. The interaction of Methanocaldococcus jannaschii (Mj) Nop56/58 with the methyltransferase fibrillarin has been investigated using site-directed mutagenesis of specific amino acids in the N-terminal domain of Nop56/58 that interacts with fibrillarin. Extensive mutagenesis revealed an unusually strong Nop56/58-fibrillarin interaction. Only deletion of the NTD itself prevented dimerization with fibrillarin. The extreme stability of the Nop56/58-fibrillarin heterodimer was confirmed in both chemical and thermal denaturation analyses. However, mutations that did not affect Nop56/58 binding to fibrillarin or sRNP assembly nevertheless disrupted sRNP-guided nucleotide modification, revealing a role for Nop56/58 in methyltransferase activity. This conclusion was supported with the cross-linking of Nop56/58 to the target RNA substrate. The Mj Nop56/58 NTD was further characterized by solving its three-dimensional crystal structure to a resolution of 1.7 Å. Despite low primary sequence conservation among the archaeal Nop56/58 homologs, the overall structure of the archaeal NTD domain is very well conserved. In conclusion, the archaeal Nop56/58 NTD exhibits a conserved domain structure whose exceptionally stable interaction with fibrillarin plays a role in both RNP assembly and methyltransferase activity.

Identification of the enzyme responsible for N1-methylation of pseudouridine 54 in archaeal tRNAs.

tRNAs from all three kingdoms of life contain a variety of modified nucleotides required for their stability, proper folding, and accurate decoding. One prominent example is the eponymous ribothymidine (rT) modification at position 54 in the T-arm of eukaryotic and bacterial tRNAs. In contrast, in most archaea this position is occupied by another hypermodified nucleotide: the isosteric N1-methylated pseudouridine. While the enzyme catalyzing pseudouridine formation at this position is known, the pseudouridine N1-specific methyltransferase responsible for this modification has not yet been experimentally identified. Here, we present biochemical and genetic evidence that the two homologous proteins, Mja_1640 (COG 1901, Pfam DUF358) and Hvo_1989 (Pfam DUF358) from Methanocaldococcus jannaschii and Haloferax volcanii, respectively, are representatives of the methyltransferase responsible for this modification. However, the in-frame deletion of the pseudouridine N1-methyltransferase gene in H. volcanii did not result in a discernable phenotype in line with similar observations for knockouts of other T-arm methylating enzymes.

The archaeal COG1901/DUF358 SPOUT-methyltransferase members, together with pseudouridine synthase Pus10, catalyze the formation of 1-methylpseudouridine at position 54 of tRNA.

The methylation of pseudouridine (Ψ) at position 54 of tRNA, producing m(1)Ψ, is a hallmark of many archaeal species, but the specific methylase involved in the formation of this modification had yet to be characterized. A comparative genomics analysis had previously identified COG1901 (DUF358), part of the SPOUT superfamily, as a candidate for this missing methylase family. To test this prediction, the COG1901 encoding gene, HVO_1989, was deleted from the Haloferax volcanii genome. Analyses of modified base contents indicated that while m(1)Ψ was present in tRNA extracted from the wild-type strain, it was absent from tRNA extracted from the mutant strain. Expression of the gene encoding COG1901 from Halobacterium sp. NRC-1, VNG1980C, complemented the m(1)Ψ minus phenotype of the ΔHVO_1989 strain. This in vivo validation was extended with in vitro tests. Using the COG1901 recombinant enzyme from Methanocaldococcus jannaschii (Mj1640), purified enzyme Pus10 from M. jannaschii and full-size tRNA transcripts or TΨ-arm (17-mer) fragments as substrates, the sequential pathway of m(1)Ψ54 formation in Archaea was reconstituted. The methylation reaction is AdoMet dependent. The efficiency of the methylase reaction depended on the identity of the residue at position 55 of the TΨ-loop. The presence of Ψ55 allowed the efficient conversion of Ψ54 to m(1)Ψ54, whereas in the presence of C55, the reaction was rather inefficient and no methylation reaction occurred if a purine was present at this position. These results led to renaming the Archaeal COG1901 members as TrmY proteins.

Nucleotide correlation based measure for identifying origin of replication in genomic sequences.

Computational prediction of the origin of replication is a challenging problem and of immense interest to biologists. Several methods have been proposed for identifying the replicon site for various classes of organisms. However, these methods have limited applicability since the replication mechanism is different in different organisms. We propose a correlation measure and show that it is correctly able to predict the origin of replication in most of the bacterial genomes. When applied to Methanocaldococcus jannaschii, Plasmodium falciparum apicoplast and Nicotiana tabacum plastid, this correlation based method is able to correctly predict the origin of replication whereas the generally used GC skew measure fails. Thus, this correlation based measure is a novel and promising tool for predicting the origin of replication in a wide class of organisms. This could have important implications in not only gaining a deeper understanding of the replication machinery in higher organisms, but also for drug discovery.

Structure and activity of the Cas3 HD nuclease MJ0384, an effector enzyme of the CRISPR interference.

Clustered regularly interspaced short palindromic repeats (CRISPRs) and Cas proteins represent an adaptive microbial immunity system against viruses and plasmids. Cas3 proteins have been proposed to play a key role in the CRISPR mechanism through the direct cleavage of invasive DNA. Here, we show that the Cas3 HD domain protein MJ0384 from Methanocaldococcus jannaschii cleaves endonucleolytically and exonucleolytically (3'-5') single-stranded DNAs and RNAs, as well as 3'-flaps, splayed arms, and R-loops. The degradation of branched DNA substrates by MJ0384 is stimulated by the Cas3 helicase MJ0383 and ATP. The crystal structure of MJ0384 revealed the active site with two bound metal cations and together with site-directed mutagenesis suggested a catalytic mechanism. Our studies suggest that the Cas3 HD nucleases working together with the Cas3 helicases can completely degrade invasive DNAs through the combination of endo- and exonuclease activities.

An N-terminal protein degradation tag enables robust selection of highly active enzymes.

Degradation tags are short peptide sequences that target proteins for destruction by housekeeping proteases. We previously utilized the C-terminal SsrA tag in directed evolution experiments to decrease the intracellular lifetime of a growth-limiting enzyme and thereby facilitate selection of highly active variants. In this study, we examine the N-terminal RepA tag as an alternative degradation signal for laboratory evolution. Although RepA proved to be less effective than SsrA at lowering protein concentrations in the cell, its N-terminal location dramatically reduced the occurrence of truncation and frameshift artifacts in selection experiments. We exploited this improvement to evolve a topologically redesigned chorismate mutase that is intrinsically disordered but already highly active for the conversion of chorismate to prephenate. After three rounds of mutagenesis and high-stringency selection, a robust and more nativelike variant was obtained that exhibited a catalytic efficiency (k(cat)/K(M) = 84000 M(-1) s(-1)) comparable to that of a natural dimeric chorismate mutase. Because of concomitant increases in catalyst yield, the level of intracellular prephenate production increased approximately 30-fold overall over the course of evolution.

Metalloproteins and the pyrite-based origin of life: a critical assessment.

We critically examine the proposal by Wächtershäuser (Prokaryotes 1:275-283, 2006a, Philos Trans R Soc Lond B Biol Sci 361: 787-1808, 2006b) that putative transition metal binding sites in protein components of the translation machinery of hyperthermophiles provide evidence of a direct relationship with the FeS clusters of pyrite and thus indicate an autotrophic origin of life in volcanic environments. Analysis of completely sequenced cellular genomes of Bacteria, Archaea and Eucarya does not support the suggestion by Wächtershäuser (Prokaryotes 1:275-283, 2006a, Philos Trans R Soc Lond B Biol Sci 361: 787-1808, 2006b) that aminoacyl-tRNA synthetases and ribosomal proteins bear sequence signatures typical of strong covalent metal bonding whose absence in mesophilic species reveals a process of adaptation towards less extreme environments.

Importance of single molecular determinants in the fidelity of expanded genetic codes.

The site-selective encoding of noncanonical amino acids (NAAs) is a powerful technique for the installation of novel chemical functional groups in proteins. This is often achieved by recoding a stop codon and requires two additional components: an evolved aminoacyl tRNA synthetase (AARS) and a cognate tRNA. Analysis of the most successful AARSs reveals common characteristics. The highest fidelity NAA systems derived from the Methanocaldococcus jannaschii tyrosyl AARS feature specific mutations to two residues reported to interact with the hydroxyl group of the substrate tyrosine. We demonstrate that the restoration of just one of these determinants for amino acid specificity results in the loss of fidelity as the evolved AARSs become noticeably promiscuous. These results offer a partial explanation of a recently retracted strategy for the synthesis of glycoproteins. Similarly, we reinvestigated a tryptophanyl AARS reported to allow the site-selective incorporation of 5-hydroxy tryptophan within mammalian cells. In multiple experiments, the enzyme displayed elements of promiscuity despite its previous characterization as a high fidelity enzyme. Given the many similarities of the TyrRSs and TrpRSs reevaluated here, our findings can be largely combined, and in doing so they reinforce the long-established central dogma regarding the molecular basis by which these enzymes contribute to the fidelity of translation. Thus, our view is that the central claims of fidelity reported in several NAA systems remain unproven and unprecedented.