Proteomic and transcriptomic analysis of selenium utilization in Methanococcus maripaludis

Methanococcus maripaludis utilizes selenocysteine- (Sec-) containing proteins (selenoproteins), mostly active in the organism's primary energy metabolism, methanogenesis. During selenium depletion, M. maripaludis employs a set of enzymes containing cysteine (Cys) instead of Sec. The genes coding for these Sec-/Cys-containing isoforms were the only genes known of which expression is influenced by the selenium status of the cell. Using proteomics and transcriptomics, approx. 7% and 12%, respectively, of all genes/proteins were found differentially expressed/synthesized in response to the selenium supply. Some of the genes identified involve methanogenesis, nitrogenase functions, and putative transporters. An increase of transcript abundance for putative transporters under selenium depletion indicated the organism's effort to tap into alternative sources of selenium. M. maripaludis is known to utilize selenite and dimethylselenide as selenium sources. To expand this list, a selenium-responsive reporter strain was assessed with nine other, environmentally relevant selenium species. While the effect of some was very similar to that of selenite, others were effectively utilized at lower concentrations. Conversely, selenate and seleno-amino acids were only utilized at unphysiologically high concentrations and two compounds were not utilized at all. To address the role of the selenium-regulated putative transporters, M. maripaludis mutant strains lacking one or two of the putative transporters were tested for the capability to utilize the different selenium species. Of the five putative transporters analyzed by loss-of-function mutagenesis, none appeared to be absolutely required for utilizing any of the selenium species tested, indicating they have redundant and/or overlapping specificities or are not dedicated selenium transporters.

IMPORTANCE

While selenium metabolism in microorganisms has been studied intensively in the past, global gene expression approaches have not been employed so far. Furthermore, the use of different selenium sources, widely environmentally interconvertible via biotic and abiotic processes, was also not extensively studied before. Methanococcus maripaludis JJ is ideally suited for such analyses, thanks to its known selenium usage and available genetic tools. Thus, an overall view on the selenium regulon of M. maripaludis was obtained via transcriptomic and proteomic analyses, which inspired further experimentation. This led to demonstrating the use of selenium sources M. maripaludis was previously not known to employ. Also, an attempt-although so far unsuccessful-was made to pinpoint potential selenium transporter genes, in order to deepen our understanding of trace element utilization in this important model organism.

Microbial electrosynthesis of single cell protein and methane by coupling fast-growing Methanococcus maripaludis with microbial electrolysis cells

Single cell protein (SCP) is a promising alternative protein source, as its production bypasses the disadvantages of animal protein production in industrial agriculture. Coupling a fast-growing hydrogen consuming organism with microbial electrolysis cells (MECs) could be a viable method for SCP production. In this study, a fast-growing and protein-rich methanogen, Methanococcus maripaludis was selected as the primary SCP source. The inoculation of M. maripaludis in MECs triggered cell synthesis with methane production. The doubling time measured was 11.2 h and the specific growth rate was 0.062 1/h. The highest SCP production rate was 13.7 mg/L/h. In the dried biomass, the weight of protein was over 60 %. Amino acid profiling of the harvested biomass demonstrated high abundance of essential amino acids. The electron flux analysis indicated that 31.3 % electrons in the electrochemical systems were directed into SCP synthesis. These results illustrated the potential for SCP production by coupling a fast-growing methanogen with MECs.

Functionally redundant formate dehydrogenases enable formate-dependent growth in Methanococcus maripaludis

Methanogens are essential for the complete remineralization of organic matter in anoxic environments. Most cultured methanogens are hydrogenotrophic, using H2 as an electron donor to reduce CO2 to CH4, but in the absence of H2 many can also use formate. Formate dehydrogenase (Fdh) is essential for formate oxidation, where it transfers electrons for the reduction of coenzyme F420 or to a flavin-based electron bifurcating reaction catalyzed by heterodisulfide reductase (Hdr), the terminal reaction of methanogenesis. Furthermore, methanogens that use formate encode at least two isoforms of Fdh in their genomes, but how these different isoforms participate in methanogenesis is unknown. Using Methanococcus maripaludis, we undertook a biochemical characterization of both Fdh isoforms involved in methanogenesis. Both Fdh1 and Fdh2 interacted with Hdr to catalyze the flavin-based electron bifurcating reaction, and both reduced F420 at similar rates. F420 reduction preceded flavin-based electron bifurcation activity for both enzymes. In a Δfdh1 mutant background, a suppressor mutation was required for Fdh2 activity. Genome sequencing revealed that this mutation resulted in the loss of a specific molybdopterin transferase (moeA), allowing for Fdh2-dependent growth, and the metal content of the proteins suggested that isoforms are dependent on either molybdenum or tungsten for activity. These data suggest that both isoforms of Fdh are functionally redundant, but their activities in vivo may be limited by gene regulation or metal availability under different growth conditions. Together these results expand our understanding of formate oxidation and the role of Fdh in methanogenesis.

A robust genetic toolbox for fine-tuning gene expression in the CO2-Fixing methanogenic archaeon Methanococcus maripaludis

Libraries of well-characterized genetic elements for fine-tuning gene expression are essential for biological and biotechnological research and applications. The fast-growing and genetically tractable methanogen, Methanococcus maripaludis, is a promising host organism for biotechnological conversion of carbon dioxide and renewable hydrogen into fuels and value-added products, as well as fundamental biological studies of archaea. However, the lack of molecular tools for gene expression has hindered its application as a workhorse to fine-tune gene and metabolic pathway expressions. In this study, we developed a genetic toolbox, including libraries of promoters, ribosome binding sites (RBS), and neutral sites for chromosomal integration, to facilitate precise gene expression in M. maripaludis. We generated a promoter library consisting of 81 constitutive promoters with expression strengths spanning a ∼104-fold dynamic range. Importantly, we identified a base composition rule for strong archaeal promoters and successfully remodeled weak promoters, enhancing their activities by up to 120-fold. We also established an RBS library containing 42 diverse RBS sequences with translation strengths covering a ∼100-fold dynamic range. Additionally, we identified eight neutral sites and developed a one-step, Cas9-based marker-less knock-in approach for chromosomal integration. We successfully applied the characterized promoter and RBS elements to significantly improve recombinant protein expression by 41-fold and modulate essential gene expression to generate corresponding physiological changes in M. maripaludis. Therefore, this work establishes a solid foundation for utilizing this autotrophic methanogen as an ideal workhorse for archaeal biology and biotechnological studies and applications.

Growth rate-dependent coordination of catabolism and anabolism in the archaeon Methanococcus maripaludis under phosphate limitation

Catabolic and anabolic processes are finely coordinated in microorganisms to provide optimized fitness under varying environmental conditions. Understanding this coordination and the resulting physiological traits reveals fundamental strategies of microbial acclimation. Here, we characterized the system-level physiology of Methanococcus maripaludis, a niche-specialized methanogenic archaeon, at different dilution rates ranging from 0.09 to 0.003 h-1 in chemostat experiments under phosphate (i.e., anabolic) limitation. Phosphate was supplied as the limiting nutrient, while formate was supplied in excess as the catabolic substrate and carbon source. We observed a decoupling of catabolism and anabolism resulting in lower biomass yield relative to catabolically limited cells at the same dilution rates. In addition, the mass abundance of several coarse-grained proteome sectors (i.e., combined abundance of proteins grouped based on their function) exhibited a linear relationship with growth rate, mostly ribosomes and their biogenesis. Accordingly, cellular RNA content also correlated with growth rate. Although the methanogenesis proteome sector was invariant, the metabolic capacity for methanogenesis, measured as methane production rates immediately after transfer to batch culture, correlated with growth rate suggesting translationally independent regulation that allows cells to only increase catabolic activity under growth-permissible conditions. These observations are in stark contrast to the physiology of M. maripaludis under formate (i.e., catabolic) limitation, where cells keep an invariant proteome including ribosomal content and a high methanogenesis capacity across a wide range of growth rates. Our findings reveal that M. maripaludis employs fundamentally different strategies to coordinate global physiology during anabolic phosphate and catabolic formate limitation.

Proteomic Analysis of Methanococcus voltae Grown in the Presence of Mineral and Nonmineral Sources of Iron and Sulfur

Iron sulfur (Fe-S) proteins are essential and ubiquitous across all domains of life, yet the mechanisms underpinning assimilation of iron (Fe) and sulfur (S) and biogenesis of Fe-S clusters are poorly understood. This is particularly true for anaerobic methanogenic archaea, which are known to employ more Fe-S proteins than other prokaryotes. Here, we utilized a deep proteomics analysis of Methanococcus voltae A3 cultured in the presence of either synthetic pyrite (FeS2) or aqueous forms of ferrous iron and sulfide to elucidate physiological responses to growth on mineral or nonmineral sources of Fe and S. The liquid chromatography-mass spectrometry (LCMS) shotgun proteomics analysis included 77% of the predicted proteome. Through a comparative analysis of intra- and extracellular proteomes, candidate proteins associated with FeS2 reductive dissolution, Fe and S acquisition, and the subsequent transport, trafficking, and storage of Fe and S were identified. The proteomic response shows a large and balanced change, suggesting that M. voltae makes physiological adjustments involving a range of biochemical processes based on the available nutrient source. Among the proteins differentially regulated were members of core methanogenesis, oxidoreductases, membrane proteins putatively involved in transport, Fe-S binding ferredoxin and radical S-adenosylmethionine proteins, ribosomal proteins, and intracellular proteins involved in Fe-S cluster assembly and storage. This work improves our understanding of ancient biogeochemical processes and can support efforts in biomining of minerals. IMPORTANCE Clusters of iron and sulfur are key components of the active sites of enzymes that facilitate microbial conversion of light or electrical energy into chemical bonds. The proteins responsible for transporting iron and sulfur into cells and assembling these elements into metal clusters are not well understood. Using a microorganism that has an unusually high demand for iron and sulfur, we conducted a global investigation of cellular proteins and how they change based on the mineral forms of iron and sulfur. Understanding this process will answer questions about life on early earth and has application in biomining and sustainable sources of energy.

Fine-Tuning Protein Self-Organization by Orthogonal Chemo-Optogenetic Tools

A universal gain-of-function approach for the spatiotemporal control of protein activity is highly desirable when reconstituting biological modules in vitro. Here we used orthogonal translation with a photocaged amino acid to map and elucidate molecular mechanisms in the self-organization of the prokaryotic filamentous cell-division protein (FtsZ) that is highly relevant for the assembly of the division ring in bacteria. We masked a tyrosine residue of FtsZ by site-specific incorporation of a photocaged tyrosine analogue. While the mutant still shows self-assembly into filaments, dynamic self-organization into ring patterns can no longer be observed. UV-mediated uncaging revealed that tyrosine 222 is essential for the regulation of the protein's GTPase activity, self-organization, and treadmilling dynamics. Thus, the light-mediated assembly of functional protein modules appears to be a promising minimal-regulation strategy for building up molecular complexity towards a minimal cell.

Performance of different methanogenic species for the microbial electrosynthesis of methane from carbon dioxide

Microbial electrosynthesis (MES) is a promising technology to convert CO₂ and electricity into the biofuel methane using methanogens. Until now, most investigations on electro-methanogenesis are "proof-of-principle" studies. In this paper, different strains were quantitatively compared in regard to final methane concentration, yields based on CO₂-conversion, productivities as well as Coulombic efficiencies in order to identify suitable organisms for MES. Methanococcus vannielii, Methanococcus maripaludis, Methanolacinia petrolearia, Methanobacterium congolense, and Methanoculleus submarinus were able to produce methane via MES at -700 mV vs. standard hydrogen electrode (SHE). Beside methane also biological H₂ production was detected during MES, which might be due to the involvement of hydrogenases. A direct electron transfer pathway is most likely. Obviously, M. maripaludis is the most resource efficient methane producer in microbial electrosynthesis regarding the methane productivity (8.81 ± 0.51 mmol m^-2 d^-1) and the Coulombic efficiency (58.9 ± 0.8%).

Biological methane production under putative Enceladus-like conditions

The detection of silica-rich dust particles, as an indication for ongoing hydrothermal activity, and the presence of water and organic molecules in the plume of Enceladus, have made Saturn's icy moon a hot spot in the search for potential extraterrestrial life. Methanogenic archaea are among the organisms that could potentially thrive under the predicted conditions on Enceladus, considering that both molecular hydrogen (H2) and methane (CH4) have been detected in the plume. Here we show that a methanogenic archaeon, Methanothermococcus okinawensis, can produce CH4 under physicochemical conditions extrapolated for Enceladus. Up to 72% carbon dioxide to CH4 conversion is reached at 50 bar in the presence of potential inhibitors. Furthermore, kinetic and thermodynamic computations of low-temperature serpentinization indicate that there may be sufficient H2 gas production to serve as a substrate for CH4 production on Enceladus. We conclude that some of the CH4 detected in the plume of Enceladus might, in principle, be produced by methanogens.

Crystallization and biochemical characterization of an archaeal lectin from Methanococcus voltae A3

A lectin from Methanococcus voltae A3 has been cloned, expressed, purified and characterized. The lectin appears to be specific for complex sugars. The protein crystallized in a tetragonal space group, with around 16 subunits in the asymmetric unit. Sequence comparisons indicate the lectin to have a β-prism I fold, with poor homology to lectins of known three-dimensional structure.

A stable genetic polymorphism underpinning microbial syntrophy

Syntrophies are metabolic cooperations, whereby two organisms co-metabolize a substrate in an interdependent manner. Many of the observed natural syntrophic interactions are mandatory in the absence of strong electron acceptors, such that one species in the syntrophy has to assume the role of electron sink for the other. While this presents an ecological setting for syntrophy to be beneficial, the potential genetic drivers of syntrophy remain unknown to date. Here, we show that the syntrophic sulfate-reducing species Desulfovibrio vulgaris displays a stable genetic polymorphism, where only a specific genotype is able to engage in syntrophy with the hydrogenotrophic methanogen Methanococcus maripaludis. This 'syntrophic' genotype is characterized by two genetic alterations, one of which is an in-frame deletion in the gene encoding for the ion-translocating subunit cooK of the membrane-bound COO hydrogenase. We show that this genotype presents a specific physiology, in which reshaping of energy conservation in the lactate oxidation pathway enables it to produce sufficient intermediate hydrogen for sustained M. maripaludis growth and thus, syntrophy. To our knowledge, these findings provide for the first time a genetic basis for syntrophy in nature and bring us closer to the rational engineering of syntrophy in synthetic microbial communities.

Effects of salinity on simultaneous reduction of perchlorate and nitrate in a methane-based membrane biofilm reactor

This study builds upon prior work showing that methane (CH₄) could be utilized as the sole electron donor and carbon source in a membrane biofilm reactor (MBfR) for complete perchlorate (ClO₄^-) and nitrate (NO₃^-) removal. Here, we further investigated the effects of salinity on the simultaneous removal of the two contaminants in the reactor. By testing ClO₄^- and NO₃^- at different salinities, we found that the reactor performance was very sensitive to salinity. While 0.2 % salinity did not significantly affect the hydrogen-based MBfR for ClO₄^- and NO₃^- removals, 1 % salinity completely inhibited ClO₄^- reduction and significantly lowered NO₃^- reduction in the CH₄-based MBfR. In salinity-free conditions, NO₃^- and ClO₄^- removal fluxes were 0.171 g N/m²-day and 0.091 g/m²-day, respectively, but NO₃^- removal fluxes dropped to 0.0085 g N/m²-day and ClO₄^- reduction was completely inhibited when the medium changed to 1 % salinity. Scanning electron microscopy (SEM) showed that the salinity dramatically changed the microbial morphology, which led to the development of wire-like cell structures. Quantitative real-time PCR (qPCR) indicated that the total number of microorganisms and abundances of functional genes significantly declined in the presence of NaCl. The relative abundances of Methylomonas (methanogens) decreased from 31.3 to 5.9 % and Denitratisoma (denitrifiers) decreased from 10.6 to 4.4 % when 1 % salinity was introduced.

Identification of the first transcriptional activator of an archaellum operon in a euryarchaeon

The archaellum is the swimming organelle of the third domain, the Archaea. In the euryarchaeon Methanococcus maripaludis, genes involved in archaella formation, including the three archaellins flaB1, flaB2 and flaB3, are mainly located in the fla operon. Previous studies have shown that transcription of fla genes and expression of Fla proteins are regulated under different growth conditions. In this study, we identify MMP1718 as the first transcriptional activator that directly regulates the fla operon in M. maripaludis. Mutants carrying an in-frame deletion in mmp1718 did not express FlaB2 detected by western blotting. Quantitative reverse transcription PCR analysis of purified RNA from the Δmmp1718 mutant showed that transcription of flaB2 was negligible compared to wildtype cells. In addition, no archaella were observed on the cell surface of the Δmmp1718 mutant. FlaB2 expression and archaellation were restored when the Δmmp1718 mutant was complemented with mmp1718 in trans. Electrophoretic motility shift assay and isothermal titration calorimetry results demonstrated the specific binding of purified MMP1718 to DNA fragments upstream of the fla promoter. Four 6 bp consensus sequences were found immediately upstream of the fla promoter and are considered the putative MMP1718-binding sites. Herein, we designate MMP1718 as EarA, the first euryarchaeal archaellum regulator.

Metabolic processes of Methanococcus maripaludis and potential applications

Methanococcus maripaludis is a rapidly growing, fully sequenced, genetically tractable model organism among hydrogenotrophic methanogens. It has the ability to convert CO2 and H2 into a useful cleaner energy fuel (CH4). In fact, this conversion enhances in the presence of free nitrogen as the sole nitrogen source due to prolonged cell growth. Given the global importance of GHG emissions and climate change, diazotrophy can be attractive for carbon capture and utilization applications from appropriately treated flue gases, where surplus hydrogen is available from renewable electricity sources. In addition, M. maripaludis can be engineered to produce other useful products such as terpenoids, hydrogen, methanol, etc. M. maripaludis with its unique abilities has the potential to be a workhorse like Escherichia coli and S. cerevisiae for fundamental and experimental biotechnology studies. More than 100 experimental studies have explored different specific aspects of the biochemistry and genetics of CO2 and N2 fixation by M. maripaludis. Its genome-scale metabolic model (iMM518) also exists to study genetic perturbations and complex biological interactions. However, a comprehensive review describing its cell structure, metabolic processes, and methanogenesis is still lacking in the literature. This review fills this crucial gap. Specifically, it integrates distributed information from the literature to provide a complete and detailed view for metabolic processes such as acetyl-CoA synthesis, pyruvate synthesis, glycolysis/gluconeogenesis, reductive tricarboxylic acid (RTCA) cycle, non-oxidative pentose phosphate pathway (NOPPP), nitrogen metabolism, amino acid metabolism, and nucleotide biosynthesis. It discusses energy production via methanogenesis and its relation to metabolism. Furthermore, it reviews taxonomy, cell structure, culture/storage conditions, molecular biology tools, genome-scale models, and potential industrial and environmental applications. Through the discussion, it develops new insights and hypotheses from experimental and modeling observations, and identifies opportunities for further research and applications.

Flux measurements and maintenance energy for carbon dioxide utilization by Methanococcus maripaludis

BACKGROUND

The rapidly growing mesophilic methanogen Methanococcus maripaludis S2 has a unique ability to consume both CO2 and N2, the main components of a flue gas, and produce methane with H2 as the electron donor. The existing literature lacks experimental measurements of CO2 and H2 uptake rates and CH4 production rates on M. maripaludis. Furthermore, it lacks estimates of maintenance energies for use with genome-scale models. In this paper, we performed batch culture experiments on M. maripaludis S2 using CO2 as the sole carbon substrate to quantify three key extracellular fluxes (CO2, H2, and CH4) along with specific growth rates. For precise computation of these fluxes from experimental measurements, we developed a systematic process simulation approach. Then, using an existing genome-scale model, we proposed an optimization procedure to estimate maintenance energy parameters: growth associated maintenance (GAM) and non-growth associated maintenance (NGAM).

RESULTS

The measured extracellular fluxes for M. maripaludis showed excellent agreement with in silico predictions from a validated genome-scale model (iMM518) for NGAM = 7.836 mmol/gDCW/h and GAM = 27.14 mmol/gDCW. M. maripaludis achieved a CO2 to CH4 conversion yield of 70-95 % and a growth yield of 3.549 ± 0.149 g DCW/mol CH4 during the exponential phase. The ATP gain of 0.35 molATP/molCH4 for M. maripaludis, computed using NGAM, is in the acceptable range of 0.3-0.7 mol ATP/molCH4 reported for methanogens. Interestingly, the uptake distribution of amino acids, quantified using iMM518, confirmed alanine to be the most preferred amino acids for growth and methanogenesis.

CONCLUSIONS

This is the first study to report experimental gas consumption and production rates for the growth of M. maripaludis on CO2 and H2 in minimal media. A systematic process simulation and optimization procedure was successfully developed to precisely quantify extracellular fluxes along with cell growth and maintenance energy parameters. Our growth yields, ATP gain, and energy parameters fall within acceptable ranges known in the literature for hydrogenotrophic methanogens.

Study of the Binding Energies between Unnatural Amino Acids and Engineered Orthogonal Tyrosyl-tRNA Synthetases

We utilized several computational approaches to evaluate the binding energies of tyrosine (Tyr) and several unnatural Tyr analogs, to several orthogonal aaRSes derived from Methanocaldococcus jannaschii and Escherichia coli tyrosyl-tRNA synthetases. The present study reveals the following: (1) AutoDock Vina and ROSETTA were able to distinguish binding energy differences for individual pairs of favorable and unfavorable aaRS-amino acid complexes, but were unable to cluster together all experimentally verified favorable complexes from unfavorable aaRS-Tyr complexes; (2) MD-MM/PBSA provided the best prediction accuracy in terms of clustering favorable and unfavorable enzyme-substrate complexes, but also required the highest computational cost; and (3) MM/PBSA based on single energy-minimized structures has a significantly lower computational cost compared to MD-MM/PBSA, but still produced sufficiently accurate predictions to cluster aaRS-amino acid interactions. Although amino acid-aaRS binding is just the first step in a complex series of processes to acylate a tRNA with its corresponding amino acid, the difference in binding energy, as shown by MD-MM/PBSA, is important for a mutant orthogonal aaRS to distinguish between a favorable unnatural amino acid (unAA) substrate from unfavorable natural amino acid substrates. Our computational study should assist further designing and engineering of orthogonal aaRSes for the genetic encoding of novel unAAs.

Effects of N-glycosylation site removal in archaellins on the assembly and function of archaella in Methanococcus maripaludis

In Methanococcus maripaludis S2, the swimming organelle, the archaellum, is composed of three archaellins, FlaB1_S2, FlaB2_S2 and FlaB3_S2. All three are modified with an N-linked tetrasaccharide at multiple sites. Disruption of the N-linked glycosylation pathway is known to cause defects in archaella assembly or function. Here, we explored the potential requirement of N-glycosylation of archaellins on archaellation by investigating the effects of eliminating the 4 N-glycosylation sites in the wildtype FlaB2_S2 protein in all possible combinations either by Asn to Glu (N to Q) substitution or Asn to Asp (N to D) substitutions of the N-glycosylation sequon asparagine. The ability of these mutant derivatives to complement a non-archaellated ΔflaB2_S2 strain was examined by electron microscopy (for archaella assembly) and swarm plates (for analysis of swimming). Western blot results showed that all mutated FlaB2_S2 proteins were expressed and of smaller apparent molecular mass compared to wildtype FlaB2_S2, consistent with the loss of glycosylation sites. In the 8 single-site mutant complements, archaella were observed on the surface of Q2, D2 and D4 (numbers after N or Q refer to the 1^st to 4^th glycosylation site). Of the 6 double-site mutation complementations all were archaellated except D1,3. Of the 4 triple-site mutation complements, only D2,3,4 was archaellated. Elimination of all 4 N-glycosylation sites resulted in non-archaellated cells, indicating some minimum amount of archaellin glycosylation was necessary for their incorporation into stable archaella. All complementations that led to a return of archaella also resulted in motile cells with the exception of the D4 version. In addition, a series of FlaB2_S2 scanning deletions each missing 10 amino acids was also generated and tested for their ability to complement the ΔflaB2_S2 strain. While most variants were expressed, none of them restored archaellation, although FlaB2_S2 harbouring a smaller 3-amino acid deletion was able to partially restore archaellation.

A peptide-based method for 13C Metabolic Flux Analysis in microbial communities

The study of intracellular metabolic fluxes and inter-species metabolite exchange for microbial communities is of crucial importance to understand and predict their behaviour. The most authoritative method of measuring intracellular fluxes, ¹³C Metabolic Flux Analysis (¹³C MFA), uses the labeling pattern obtained from metabolites (typically amino acids) during ¹³C labeling experiments to derive intracellular fluxes. However, these metabolite labeling patterns cannot easily be obtained for each of the members of the community. Here we propose a new type of ¹³C MFA that infers fluxes based on peptide labeling, instead of amino acid labeling. The advantage of this method resides in the fact that the peptide sequence can be used to identify the microbial species it originates from and, simultaneously, the peptide labeling can be used to infer intracellular metabolic fluxes. Peptide identity and labeling patterns can be obtained in a high-throughput manner from modern proteomics techniques. We show that, using this method, it is theoretically possible to recover intracellular metabolic fluxes in the same way as through the standard amino acid based ¹³C MFA, and quantify the amount of information lost as a consequence of using peptides instead of amino acids. We show that by using a relatively small number of peptides we can counter this information loss. We computationally tested this method with a well-characterized simple microbial community consisting of two species.

Microbial communities underlie a variety of important biochemical processes ranging from underground cave formation to gold mining or the onset of obesity. Metabolic fluxes describe how carbon and energy flow through the microbial community and therefore provide insights that are rarely captured by other techniques, such as metatranscriptomics or metaproteomics. The most authoritative method to measure fluxes for pure cultures consists of feeding the cells a labeled carbon source and deriving the fluxes from the ensuing metabolite labeling pattern (typically amino acids). Since we cannot easily separate cells of metabolite for each species in a community, this approach is not generally applicable to microbial communities. Here we present a method to derive fluxes from the labeling of peptides, instead of amino acids. This approach has the advantage that peptides can be assigned to each species in a community in a high-throughput fashion through modern proteomic methods. We show that, by using this method, it is theoretically possible to recover the same amount of information as through the standard approach, if enough peptides are used. We computationally tested this method with a well-characterized simple microbial community consisting of two species.

Identification of an additional minor pilin essential for piliation in the archaeon Methanococcus maripaludis

Methanococcus maripaludis is an archaeon with two studied surface appendages, archaella and type IV-like pili. Previously, the major structural pilin was identified as MMP1685 and three additional proteins were designated as minor pilins (EpdA, EpdB and EpdC). All of the proteins are likely processed by the pilin-specific prepilin peptidase EppA. Six other genes were identified earlier as likely encoding pilin proteins processed also by EppA. In this study, each of the six genes (mmp0528, mmp0600, mmp0601, mmp0709, mmp0903 and mmp1283) was deleted and the mutants examined by electron microscopy to determine their essentiality for pili formation. While mRNA transcripts of all genes were detected by RT-PCR, only the deletion of mmp1283 led to nonpiliated cells. This strain could be complemented back to a piliated state by supplying a wildtype copy of the mmp1283 gene in trans. This study adds to the complexity of the type IV pili system in M. maripaludis and raises questions about the functions of the remaining five pilin-like genes and whether M. maripaludis under other growth conditions may be able to assemble additional pili-like structures.

Archaeology of eukaryotic DNA replication

Recent advances in the characterization of the archaeal DNA replication system together with comparative genomic analysis have led to the identification of several previously uncharacterized archaeal proteins involved in replication and currently reveal a nearly complete correspondence between the components of the archaeal and eukaryotic replication machineries. It can be inferred that the archaeal ancestor of eukaryotes and even the last common ancestor of all extant archaea possessed replication machineries that were comparable in complexity to the eukaryotic replication system. The eukaryotic replication system encompasses multiple paralogs of ancestral components such that heteromeric complexes in eukaryotes replace archaeal homomeric complexes, apparently along with subfunctionalization of the eukaryotic complex subunits. In the archaea, parallel, lineage-specific duplications of many genes encoding replication machinery components are detectable as well; most of these archaeal paralogs remain to be functionally characterized. The archaeal replication system shows remarkable plasticity whereby even some essential components such as DNA polymerase and single-stranded DNA-binding protein are displaced by unrelated proteins with analogous activities in some lineages.

Selenomodification of tRNA in archaea requires a bipartite rhodanese enzyme

5-Methylaminomethyl-2-selenouridine (mnm(5)Se(2)U) is found in the first position of the anticodon in certain tRNAs from bacteria, archaea and eukaryotes. This selenonucleoside is formed in Escherichia coli from the corresponding thionucleoside mnm(5)S(2)U by the monomeric enzyme YbbB. This nucleoside is present in the tRNA of Methanococcales, yet the corresponding 2-selenouridine synthase is unknown in archaea and eukaryotes. Here we report that a bipartite ybbB ortholog is present in all members of the Methanococcales. Gene deletions in Methanococcus maripaludis and in vitro activity assays confirm that the two proteins act in trans to form in tRNA a selenonucleoside, presumably mnm(5)Se(2)U. Phylogenetic data suggest a primal origin of seleno-modified tRNAs.

Kinetics of the association/dissociation cycle of an ATP-binding cassette nucleotide-binding domain

Most ATP binding cassette (ABC) proteins are pumps that transport substrates across biological membranes using the energy of ATP hydrolysis. Functional ABC proteins have two nucleotide-binding domains (NBDs) that bind and hydrolyze ATP, but the molecular mechanism of nucleotide hydrolysis is unresolved. This is due in part to the limited kinetic information on NBD association and dissociation. Here, we show dimerization of a catalytically active NBD and follow in real time the association and dissociation of NBDs from the changes in fluorescence emission of a tryptophan strategically located at the center of the dimer interface. Spectroscopic and structural studies demonstrated that the tryptophan can be used as dimerization probe, and we showed that under hydrolysis conditions (millimolar MgATP), not only the dimer dissociation rate increases, but also the dimerization rate. Neither dimer formation or dissociation are clearly favored, and the end result is a dynamic equilibrium where the concentrations of monomer and dimer are very similar. We proposed that based on their variable rates of hydrolysis, the rate-limiting step of the hydrolysis cycle may differ among full-length ABC proteins.

Conformational dynamics of the plug domain of the SecYEG protein-conducting channel

The central pore of the SecYEG preprotein-conducting channel is closed at the periplasmic face of the membrane by a plug domain. To study its conformational dynamics, the plug was labeled site-specifically with an environment-sensitive fluorophore. In the presence of a stable preprotein translocation inter-mediate, the SecY plug showed an enhanced solvent exposure consistent with a displacement from the hydrophobic central pore region. In contrast, binding and insertion of a ribosome-bound nascent membrane protein did not alter the plug conformation. These data indicate different plug dynamics depending on the ligand bound state of the SecYEG channel.

Mechanism of N-methylation by the tRNA m1G37 methyltransferase Trm5

Trm5 is a eukaryal and archaeal tRNA methyltransferase that catalyzes methyl transfer from S-adenosylmethionine (AdoMet) to the N(1) position of G37 directly 3' to the anticodon. While the biological role of m(1)G37 in enhancing translational fidelity is well established, the catalytic mechanism of Trm5 has remained obscure. To address the mechanism of Trm5 and more broadly the mechanism of N-methylation to nucleobases, we examined the pH-activity profile of an archaeal Trm5 enzyme, and performed structure-guided mutational analysis. The data reveal a marked dependence of enzyme-catalyzed methyl transfer on hydrogen ion equilibria: the single-turnover rate constant for methylation increases by one order of magnitude from pH 6.0 to reach a plateau at pH 7.0. This suggests a mechanism involving proton transfer from G37 as the key element in catalysis. Consideration of the kinetic data in light of the Trm5-tRNA-AdoMet ternary cocrystal structure, determined in a precatalytic conformation, suggests that proton transfer is associated with an induced fit rearrangement of the complex that precedes formation of the reactive configuration in the active site. Key roles for the conserved R145 side chain in stabilizing a proposed oxyanion at G37-O(6), and for E185 as a general base to accept the proton from G37-N(1), are suggested based on the mutational analysis.

Environments that induce synthetic microbial ecosystems

Interactions between microbial species are sometimes mediated by the exchange of small molecules, secreted by one species and metabolized by another. Both one-way (commensal) and two-way (mutualistic) interactions may contribute to complex networks of interdependencies. Understanding these interactions constitutes an open challenge in microbial ecology, with applications ranging from the human microbiome to environmental sustainability. In parallel to natural communities, it is possible to explore interactions in artificial microbial ecosystems, e.g. pairs of genetically engineered mutualistic strains. Here we computationally generate artificial microbial ecosystems without re-engineering the microbes themselves, but rather by predicting their growth on appropriately designed media. We use genome-scale stoichiometric models of metabolism to identify media that can sustain growth for a pair of species, but fail to do so for one or both individual species, thereby inducing putative symbiotic interactions. We first tested our approach on two previously studied mutualistic pairs, and on a pair of highly curated model organisms, showing that our algorithms successfully recapitulate known interactions, robustly predict new ones, and provide novel insight on exchanged molecules. We then applied our method to all possible pairs of seven microbial species, and found that it is always possible to identify putative media that induce commensalism or mutualism. Our analysis also suggests that symbiotic interactions may arise more readily through environmental fluctuations than genetic modifications. We envision that our approach will help generate microbe-microbe interaction maps useful for understanding microbial consortia dynamics and evolution, and for exploring the full potential of natural metabolic pathways for metabolic engineering applications.