Combining ancestral sequence reconstruction with protein design to identify an interface hotspot in a key metabolic enzyme complex.

It is important to identify hotspot residues that determine protein-protein interactions in interfaces of macromolecular complexes. We have applied a combination of ancestral sequence reconstruction and protein design to identify hotspots within imidazole glycerol phosphate synthase (ImGPS). ImGPS is a key metabolic enzyme complex, which links histidine and de novo purine biosynthesis and consists of the cyclase subunit HisF and the glutaminase subunit HisH. Initial fluorescence titration experiments showed that HisH from Zymomonas mobilis (zmHisH) binds with high affinity to the reconstructed HisF from the last universal common ancestor (LUCA-HisF) but not to HisF from Pyrobaculum arsenaticum (paHisF), which differ by 103 residues. Subsequent titration experiments with a reconstructed evolutionary intermediate linking LUCA-HisF and paHisF and inspection of the subunit interface of a contemporary ImGPS allowed us to narrow down the differences crucial for zmHisH binding to nine amino acids of HisF. Homology modeling and in silico mutagenesis studies suggested that at most two of these nine HisF residues are crucial for zmHisH binding. These computational results were verified by experimental site-directed mutagenesis, which finally enabled us to pinpoint a single amino acid residue in HisF that is decisive for high-affinity binding of zmHisH. Our work shows that the identification of protein interface hotspots can be very efficient when reconstructed proteins with different binding properties are included in the analysis. Proteins 2017; 85:312-321. © 2016 Wiley Periodicals, Inc.

The starch-binding domain family CBM41-An in silico analysis of evolutionary relationships.

Within the CAZy database, there are 81 carbohydrate-binding module (CBM) families. A CBM represents a non-catalytic domain in a modular arrangement of glycoside hydrolases (GHs). The present in silico study has been focused on starch-binding domains from the family CBM41 that are usually part of pullulanases from the α-amylase family GH13. Currently there are more than 1,600 sequences classified in the family CBM41, almost exclusively from Bacteria, and so a study was undertaken in an effort to divide the members into relevant groups (subfamilies) and also to contribute to the evolutionary picture of family CBM41. The CBM41 members adopt a ß-sandwich fold (∼100 residues) with one carbohydrate-binding site formed by the side-chains of three aromatic residues that interact with carbohydrate. The family CBM41 can be divided into two basic subdivisions, distinguished from each other by a characteristic sequence pattern or motif of the three essential aromatics as follows: (i) "W-W-∼10aa-W" (the so-called Streptococcus/Klebsiella-type); and (ii) "W-W-∼30aa-W" (Thermotoga-type). Based on our bioinformatics analysis it is clear that the first and second positions of the motif can be occupied by aromatic residues (Phe, Tyr, His) other than tryptophan, resulting in the existence of six different carbohydrate-binding CBM41 groups, that reflect mostly differences in taxonomy, but which should retain the ability to bind an α-glucan. In addition, three more groups have been proposed that, although lacking the crucial aromatic motif, could possibly employ other residues from remaining parts of their sequence for binding carbohydrate. Proteins 2017; 85:1480-1492. © 2017 Wiley Periodicals, Inc.

Microbial communities from 20 different hydrogen-producing reactors studied by 454 pyrosequencing.

To provide new insight into the dark fermentation process, a multi-lateral study was performed to study the microbiology of 20 different lab-scale bioreactors operated in four different countries (Brazil, Chile, Mexico, and Uruguay). Samples (29) were collected from bioreactors with different configurations, operation conditions, and performances. The microbial communities were analyzed using 16S rRNA genes 454 pyrosequencing. The results showed notably uneven communities with a high predominance of a particular genus. The phylum Firmicutes predominated in most of the samples, but the phyla Thermotogae or Proteobacteria dominated in a few samples. Genera from three physiological groups were detected: high-yield hydrogen producers (Clostridium, Kosmotoga, Enterobacter), fermenters with low-hydrogen yield (mostly from Veillonelaceae), and competitors (Lactobacillus). Inocula, reactor configurations, and substrates influence the microbial communities. This is the first joint effort that evaluates hydrogen-producing reactors and operational conditions from different countries and contributes to understand the dark fermentation process.

The structural basis of substrate promiscuity in UDP-hexose 4-epimerase from the hyperthermophilic Eubacterium Thermotoga maritima.

UDP-galactose 4-epimerase (GalE) catalyzes the interconversion of UDP-glucose (UDP-Glc) and UDP-galactose (UDP-Gal), which is a pivotal step in the Leloir pathway for d-galactose metabolism. Although GalE is widely distributed in prokaryotes and eukaryotes, little information is available regarding hyperthermophilic GalE. We overexpressed the TM0509 gene, encoding a putative GalE from Thermotoga maritima (TMGalE), in Escherichia coli and characterized the encoded protein. To further investigate the molecular basis of this enzyme's catalytic function, we determined the crystal structures of TMGalE and TMGalE bound to UDP-Glc at resolutions of 1.9 Å and 2.0 Å, respectively. The enzyme was determined to be a homodimer with a molecular mass of 70 kDa. The enzyme could reversibly catalyze the epimerization of UDP-GalNAc/UDP-GlcNAc as well as UDP-Gal/UDP-Glc at elevated temperatures, with an apparent optimal temperature and pH of 80 °C and 7.0, respectively. Our data showed that TM0509 is a UDP-galactosugar 4-epimerase involved in d-galactose metabolism; consequently, this study provides the first detailed characterization of a hyperthermophilic GalE. Moreover, the promiscuous substrate specificity of TMGalE, which is more similar to human GalE than E. coli GalE, supports the notion that TMGalE might exhibit the earliest form of sugar-epimerizing enzymes in the evolution of galactose metabolism.

Tagaturonate-fructuronate epimerase UxaE, a novel enzyme in the hexuronate catabolic network in Thermotoga maritima.

Thermotoga maritima is a marine hyperthermophilic microorganism that degrades a wide range of simple and complex carbohydrates including pectin and produces fermentative hydrogen at high yield. Galacturonate and glucuronate, two abundant hexuronic acids in pectin and xylan, respectively, are catabolized via committed metabolic pathways to supply carbon and energy for a variety of microorganisms. By a combination of bioinformatics and experimental techniques we identified a novel enzyme family (named UxaE) catalysing a previously unknown reaction in the hexuronic acid catabolic pathway, epimerization of tagaturonate to fructuronate. The enzymatic activity of the purified recombinant tagaturonate epimerase from T. maritima was directly confirmed and kinetically characterized. Its function was also confirmed by genetic complementation of the growth of the Escherichia coli uxaB knockout mutant strain on galacturonate. An inferred novel galacturonate to mannonate catabolic pathway in T. maritima was reconstituted in vitro using a mixture of recombinant purified enzymes UxaE, UxaC and UxuB. Members of the newly identified UxaE family were identified in ~50 phylogenetically diverse heterotrophic bacteria from aquatic and soil environments. The genomic context of respective genes and reconstruction of associated pathways suggest that UxaE enzymatic and biological function remains conserved in all of these species.

On the chimeric nature, thermophilic origin, and phylogenetic placement of the Thermotogales.

Since publication of the first Thermotogales genome, Thermotoga maritima strain MSB8, single- and multi-gene analyses have disagreed on the phylogenetic position of this order of Bacteria. Here we present the genome sequences of 4 additional members of the Thermotogales (Tt. petrophila, Tt. lettingae, Thermosipho melanesiensis, and Fervidobacterium nodosum) and a comprehensive comparative analysis including the original T. maritima genome. While ribosomal protein genes strongly place Thermotogales as a sister group to Aquificales, the majority of genes with sufficient phylogenetic signal show affinities to Archaea and Firmicutes, especially Clostridia. Indeed, on the basis of the majority of genes in their genomes (including genes that are also found in Aquificales), Thermotogales should be considered members of the Firmicutes. This result highlights the conflict between the taxonomic goal of assigning every species to a unique position in an inclusive Linnaean hierarchy and the evolutionary goal of understanding phylogenesis in the presence of pervasive horizontal gene transfer (HGT) within prokaryotes. Amino acid compositions of reconstructed ancestral sequences from 423 gene families suggest an origin of this gene pool even more thermophilic than extant members of this order, followed by adaptation to lower growth temperatures within the Thermotogales.

An oxygen reduction chain in the hyperthermophilic anaerobe Thermotoga maritima highlights horizontal gene transfer between Thermococcales and Thermotogales.

The hyperthermophile Thermotoga maritima, although strictly anaerobic, is able to grow in the presence of low amounts of O(2). Here, we show that this bacterium consumes O(2) via a three-partner chain involving an NADH oxidoreductase (NRO), a rubredoxin (Rd) and a flavo-diiron protein (FprA) (locus tags: TM_0754, TM_0659 and TM_0755, respectively). In vitro experiments showed that the NADH-dependent O(2) consumption rate was 881.9 (± 106.7) mol O(2) consumed min(-1) per mol of FprA at 37°C and that water was the main end-product of the reaction. We propose that this O(2) reduction chain plays a central role in the O(2) tolerance of T. maritima. Phylogenetic analyses suggest that the genes coding for these three components were acquired by an ancestor of Thermotogales from an ancestor of Thermococcales via a single gene transfer. This event likely also involved two ROS scavenging enzymes (neelaredoxin and rubrerythrin) that are encoded by genes clustered with those coding for FprA, NRO and Rd in the ancestor of Thermococcales. Such genomic organization would have provided the ancestor of Thermotogales with a complete set of enzymes dedicated to O(2)-toxicity defence. Beside Thermotogales and Thermococcales, horizontal gene transfers have played a major role in disseminating these enzymes within the hyperthermophilic anaerobic prokaryotic communities, allowing them to cope with fluctuating oxidative conditions that exist in situ.

Evaluating the potential of indigenous methanogenic consortium for enhanced oil and gas recovery from high temperature depleted oil reservoir.

In past years, lots of research has been focused on the indigenous bacteria and their mechanisms, which help in enhanced oil recovery. Most of the oil wells in Indian subcontinent have temperature higher than 60 °C. Also, the role of methanogenic consortia from high temperature petroleum reservoir for enhanced oil recovery (EOR) has not been explored much. Hence, in the present study methanogens isolated from thermophilic oil wells (70 °C) were evaluated for enhanced oil recovery. Methane gas is produced by methanogens, which helps in oil recovery from depleted oil wells through reservoir re-pressurization and also can be recovered from reservoir along with crude oil as alternative energy source. Therefore, in this study indigenous methanogenic consortium (TERIL146) was enriched from high temperature oil reservoir showing (12 mmol/l) gas production along with other metabolites. Sequencing analysis revealed the presence of Methanothermobacter sp., Thermoanaerobacter sp., Gelria sp. and Thermotoga sp. in the consortium. Furthermore, the developed indigenous consortium TERIL146 showed 8.3% incremental oil recovery in sandpack assay. The present study demonstrates successful recovery of both oil and energy (gas) by the developed indigenous methanogenic consortium TERIL146 for potential application in thermophilic depleted oil wells of Indian subcontinent.

Evidence for extensive gene flow and Thermotoga subpopulations in subsurface and marine environments.

Oil reservoirs represent a nutrient-rich ecological niche of the deep biosphere. Although most oil reservoirs are occupied by microbial populations, when and how the microbes colonized these environments remains unanswered. To address this question, we compared 11 genomes of Thermotoga maritima-like hyperthermophilic bacteria from two environment types: subsurface oil reservoirs in the North Sea and Japan, and marine sites located in the Kuril Islands, Italy and the Azores. We complemented our genomes with Thermotoga DNA from publicly available subsurface metagenomes from North America and Australia. Our analysis revealed complex non-bifurcating evolutionary history of the isolates' genomes, suggesting high amounts of gene flow across all sampled locations, a conjecture supported by numerous recombination events. Genomes from the same type of environment tend to be more similar, and have exchanged more genes with each other than with geographically close isolates from different types of environments. Hence, Thermotoga populations of oil reservoirs do not appear isolated, a requirement of the 'burial and isolation' hypothesis, under which reservoir bacteria are descendants of the isolated communities buried with sediments that over time became oil reservoirs. Instead, our analysis supports a more complex view, where bacteria from subsurface and marine populations have been continuously migrating into the oil reservoirs and influencing their genetic composition. The Thermotoga spp. in the oil reservoirs in the North Sea and Japan probably entered the reservoirs shortly after they were formed. An Australian oil reservoir, on the other hand, was likely colonized very recently, perhaps during human reservoir development.

Were the original eubacteria thermophiles?

Thermotoga maritima is one of the more unusual eubacteria: It is highly thermophilic, growing at temperatures higher than any other eubacterium; its cell wall appears to have a unique structure and its lipids a unique composition; and the organism is surrounded by a loose-fitting sheath of unknown function. Its phenotypic uniqueness is matched by its phylogenetic position; Thermotoga maritima represents the deepest known branching in the eubacterial line of descent, as measured by ribosomal RNA sequence comparisons. T. maritima also represents the most slowly evolving of eubacterial lineages. The fact that the two deepest branchings in the eubacterial line of descent (the other, the green non-sulfur bacteria and relatives, i.e. Chloroflexus, Thermomicrobium, etc.) are both basically thermophilic and slowly evolving, strongly suggests that all eubacteria have ultimately arisen from a thermophilic ancestor.

Estimation of phylogeny and invariant sites under the general Markov model of nucleotide sequence evolution.

The models of nucleotide substitution used by most maximum likelihood-based methods assume that the evolutionary process is stationary, reversible, and homogeneous. We present an extension of the Barry and Hartigan model, which can be used to estimate parameters by maximum likelihood (ML) when the data contain invariant sites and there are violations of the assumptions of stationarity, reversibility, and homogeneity. Unlike most ML methods for estimating invariant sites, we estimate the nucleotide composition of invariant sites separately from that of variable sites. We analyze a bacterial data set where problems due to lack of stationarity and homogeneity have been previously well noted and use the parametric bootstrap to show that the data are consistent with our general Markov model. We also show that estimates of invariant sites obtained using our method are fairly accurate when applied to data simulated under the general Markov model.

Chromosome evolution in the Thermotogales: large-scale inversions and strain diversification of CRISPR sequences.

In the present study, the chromosomes of two members of the Thermotogales were compared. A whole-genome alignment of Thermotoga maritima MSB8 and Thermotoga neapolitana NS-E has revealed numerous large-scale DNA rearrangements, most of which are associated with CRISPR DNA repeats and/or tRNA genes. These DNA rearrangements do not include the putative origin of DNA replication but move within the same replichore, i.e., the same replicating half of the chromosome (delimited by the replication origin and terminus). Based on cumulative GC skew analysis, both the T. maritima and T. neapolitana lineages contain one or two major inverted DNA segments. Also, based on PCR amplification and sequence analysis of the DNA joints that are associated with the major rearrangements, the overall chromosome architecture was found to be conserved at most DNA joints for other strains of T. neapolitana. Taken together, the results from this analysis suggest that the observed chromosomal rearrangements in the Thermotogales likely occurred by successive inversions after their divergence from a common ancestor and before strain diversification. Finally, sequence analysis shows that size polymorphisms in the DNA joints associated with CRISPRs can be explained by expansion and possibly contraction of the DNA repeat and spacer unit, providing a tool for discerning the relatedness of strains from different geographic locations.

The Thermotoga maritima phenotype is impacted by syntrophic interaction with Methanococcus jannaschii in hyperthermophilic coculture.

Significant growth phase-dependent differences were noted in the transcriptome of the hyperthermophilic bacterium Thermotoga maritima when it was cocultured with the hyperthermophilic archaeon Methanococcus jannaschii. For the mid-log-to-early-stationary-phase transition of a T. maritima monoculture, 24 genes (1.3% of the genome) were differentially expressed twofold or more. In contrast, methanogenic coculture gave rise to 292 genes differentially expressed in T. maritima at this level (15.5% of the genome) for the same growth phase transition. Interspecies H2 transfer resulted in three- to fivefold-higher T. maritima cell densities than in the monoculture, with concomitant formation of exopolysaccharide (EPS)-based cell aggregates. Differential expression of specific sigma factors and genes related to the ppGpp-dependent stringent response suggests involvement in the transition into stationary phase and aggregate formation. Cell aggregation was growth phase dependent, such that it was most prominent during mid-log phase and decayed as cells entered stationary phase. The reduction in cell aggregation was coincidental with down-regulation of genes encoding EPS-forming glycosyltranferases and up-regulation of genes encoding beta-specific glycosyl hydrolases; the latter were presumably involved in hydrolysis of beta-linked EPS to release cells from aggregates. Detachment of aggregates may facilitate colonization of new locations in natural environments where T. maritima coexists with other organisms. Taken together, these results demonstrate that syntrophic interactions can impact the transcriptome of heterotrophs in methanogenic coculture, and this factor should be considered in examining the microbial ecology in anaerobic environments.

Gene transfer and genome plasticity in Thermotoga maritima, a model hyperthermophilic species.

The genome sequence of the hyperthermophilic bacterium Thermotoga maritima MSB8 presents evidence for lateral gene transfer events between bacterial and archaeal species. To estimate the extent of genomic diversity across the order Thermotogales, a comparative genomic hybridization study was initiated to compare nine Thermotoga strains to the sequenced T. maritima MSB8. Many differences could be associated with substrate utilization patterns, which are most likely a reflection of the environmental niche that these individual species occupy. A detailed analysis of some of the predicted variable regions demonstrates many examples of the deletion/insertion of complete cassettes of genes and of gene rearrangements and insertions of DNA within genes, with the C or N terminus being retained. Although the mechanism for gene transfer in this lineage remains to be elucidated, this analysis suggests possible associations with repetitive elements and highlights the possible benefits of rampant genetic exchange to these species.

A highly conserved plasmid from the extreme thermophile Thermotoga maritima MC24 is a member of a family of plasmids distributed worldwide.

We have screened Thermotoga strains, isolated from hydrothermal vents near the Kuril Islands, for the presence of plasmid DNA. The miniplasmid pMC24 was isolated from the extreme thermophilic eubacteria Thermotoga maritima and sequenced, showing it to be a plasmid of 846 bp. It was found, from a search of the databases, to be closely related to the previously described Thermotoga miniplasmid pRQ7, isolated from a strain found on the Azore Islands, and was distinguished by only two point mutations. These changes resulted in two consecutive frameshifts altering a region encoding 9 amino acids in the Rep-coding region. We have also shown that pMC24, as with pRQ7, is negatively supercoiled. It seems that negatively supercoiled miniplasmids related to pRQ7 are spread worldwide and strongly maintained among Thermotoga strains.

Petrotoga mexicana sp. nov., a novel thermophilic, anaerobic and xylanolytic bacterium isolated from an oil-producing well in the Gulf of Mexico.

A novel anaerobic, thermophilic, xylanolytic, motile rod-shaped bacterium with a sheath-like outer structure (toga) was isolated from a Mexican oil well in the Gulf of Mexico. Strain MET12T was a Gram-negative bacterium, reducing elemental sulfur, thiosulfate and sulfite to hydrogen sulfide. Its optimum growth conditions were 55 degrees C, pH 6.6, 3% NaCl and 0.15% MgCl2.6H2O. The DNA G+C content was 36.1 mol%. Phylogenetically, strain MET12T was related to members of genus Petrotoga, with similarities to Petrotoga mobilis, Petrotoga sibirica, Petrotoga miotherma and Petrotoga olearia varying from 97.6 to 98.8%. However DNA-DNA relatedness values between these species and strain MET12T were lower than 70%. As strain MET12T (=DSM 14811T=CIP 107371T) was genomically and phenotypically different from existing Petrotoga species, it is proposed as the type strain of a novel species, Petrotoga mexicana sp. nov.

Evidence for lateral gene transfer between Archaea and bacteria from genome sequence of Thermotoga maritima.

The 1,860,725-base-pair genome of Thermotoga maritima MSB8 contains 1,877 predicted coding regions, 1,014 (54%) of which have functional assignments and 863 (46%) of which are of unknown function. Genome analysis reveals numerous pathways involved in degradation of sugars and plant polysaccharides, and 108 genes that have orthologues only in the genomes of other thermophilic Eubacteria and Archaea. Of the Eubacteria sequenced to date, T. maritima has the highest percentage (24%) of genes that are most similar to archaeal genes. Eighty-one archaeal-like genes are clustered in 15 regions of the T. maritima genome that range in size from 4 to 20 kilobases. Conservation of gene order between T. maritima and Archaea in many of the clustered regions suggests that lateral gene transfer may have occurred between thermophilic Eubacteria and Archaea.

Glucose-to-fructose conversion at high temperatures with xylose (glucose) isomerases from Streptomyces murinus and two hyperthermophilic Thermotoga species.

The conversion of glucose to fructose at elevated temperatures, as catalyzed by soluble and immobilized xylose (glucose) isomerases from the hyperthermophiles Thermotoga maritima (TMGI) and Thermotoga neapolitana 5068 (TNGI) and from the mesophile Streptomyces murinus (SMGI), was examined. At pH 7.0 in the presence of Mg(2+), the temperature optima for the three soluble enzymes were 85 degrees C (SMGI), 95 degrees to 100 degrees C (TNGI), and >100 degrees C (TMGI). Under certain conditions, soluble forms of the three enzymes exhibited an unusual, multiphasic inactivation behavior in which the decay rate slowed considerably after an initial rapid decline. However, the inactivation of the enzymes covalently immobilized to glass beads, monophasic in most cases, was characterized by a first-order decay rate intermediate between those of the initial rapid and slower phases for the soluble enzymes. Enzyme productivities for the three immobilized GIs were determined experimentally in the presence of Mg(2+). The highest productivities measured were 750 and 760 kg fructose per kilogram SMGI at 60 degrees C and 70 degrees C, respectively. The highest productivity for both TMGI and TNGI in the presence of Mg(2+) occurred at 70 degrees C, pH 7.0, with approximately 230 and 200 kg fructose per kilogram enzyme for TNGI and TMGI, respectively. At 80 degrees C and in the presence of Mg(2+), productivities for the three enzymes ranged from 31 to 273. A simple mathematical model, which accounted for thermal effects on kinetics, glucose-fructose equilibrium, and enzyme inactivation, was used to examine the potential for high-fructose corn syrup (HFCS) production at 80 degrees C and above using TNGI and SMGI under optimal conditions, which included the presence of both Co(2+) and Mg(2+). In the presence of both cations, these enzymes showed the potential to catalyze glucose-to-fructose conversion at 80 degrees C with estimated lifetime productivities on the order of 2000 kg fructose per kilogram enzyme, a value competitive with enzymes currently used at 55 degrees to 65 degrees C, but with the additional advantage of higher fructose concentrations. At 90 degrees C, the estimated productivity for SMGI dropped to 200, whereas, for TNGI, lifetime productivities on the order of 1000 were estimated. Assuming that the most favorable biocatalytic and thermostability features of these enzymes can be captured in immobilized form and the chemical lability of substrates and products can be minimized, HFCS production at high temperatures could be used to achieve higher fructose concentrations as well as create alternative processing strategies.