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.

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.

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.

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.

Genetic incorporation of unnatural amino acids into proteins in Mycobacterium tuberculosis

New tools are needed to study the intracellular pathogen Mycobacterium tuberculosis (Mtb), the causative agent of tuberculosis (TB), to facilitate new drug discovery and vaccine development. We have developed methodology to genetically incorporate unnatural amino acids into proteins in Mycobacterium smegmatis, BCG and Mtb, grown both extracellularly in culture and inside host cells. Orthogonal mutant tRNA^Tyr/tyrosyl-tRNA synthetase pairs derived from Methanococcus jannaschii and evolved in Escherichia coli incorporate a variety of unnatural amino acids (including photocrosslinking, chemically reactive, heavy atom containing, and immunogenic amino acids) into proteins in response to the amber nonsense codon. By taking advantage of the fidelity and suppression efficiency of the MjtRNA/pIpaRS pair in mycobacteria, we are also able to use p-iodophenylalanine to induce the expression of proteins in mycobacteria both extracellularly in culture and inside of mammalian host cells. This provides a new approach to regulate the expression of reporter genes or mycobacteria endogenous genes of interest. The establishment of the unnatural amino acid expression system in Mtb, an intracellular pathogen, should facilitate studies of TB biology and vaccine development.

Splicing of mRNA mediated by tRNA sequences in mouse cells

tRNA splicing is essential for the formation of tRNAs and therefore for gene expression. A circularly permuted sequence of an amber-suppressor pre-tRNA gene was inserted into the sequence encoding the mouse NEMO protein. We demonstrated that, in mouse cells, the hybrid pre-tRNA/pre-mRNAs can be spliced precisely at the sites of the pre-tRNA intron. This splicing reaction produces functional tRNAs that suppress amber codons as well as translatable mRNAs that sustain the NF-kappaB activation pathway. The RNA molecules extracted from mouse cells were amplified by RT-PCR, and their sequences were determined, confirming the identity of the splice junctions. We then applied the Archaea-express technology, in which an archaeal RNA endonuclease is expressed in mouse cells. We show that both the endogenous eukaryal endonuclease and the archaeal one cleave the hybrid pre-tRNA/pre-mRNAs in the same manner with an additive effect.

Tryptophan auxotrophs were obtained by random transposon insertions in the Methanococcus maripaludis tryptophan operon

Methanococcus maripaludis is an anaerobic, methane-producing archaeon that utilizes H(2) or formate for the reduction of CO(2) to methane. Tryptophan auxotrophs were constructed by in vitro insertions of the Tn5 transposon into the tryptophan operon, followed by transformation into M. maripaludis. This method could serve for rapid insertions into large cloned DNA regions.

Synthesis of a cross-linked branched polymer network in the interior of a protein cage

A goal of biomimetic chemistry is to use the hierarchical architecture inherent in biological systems to guide the synthesis of functional three-dimensional structures. Viruses and other highly symmetrical protein cage architectures provide defined scaffolds to initiate hierarchical structure assembly. Here we demonstrate that a cross-linked branched polymer can be initiated and synthesized within the interior cavity of a protein cage architecture. Creating this polymer network allows for the spatial control of pendant reactive sites and dramatically increases the stability of the cage architecture. This material was generated by the sequential coupling of multifunctional monomers using click chemistry to create a branched cross-linked polymer network. Analysis of polymer growth by mass spectrometry demonstrated that the polymer was initiated at the interior surface of the cage at genetically introduced cysteine reactive sites. The polymer grew as expected to generation 2.5 where it was limited by the size constraints of the cavity. The polymer network was fully cross-linked across protein subunits that make up the cage and extended the thermal stability for the cage to at least 120 degrees C. The introduced reactive centers were shown to be active and their number density increased with increasing generation. This synthetic approach provides a new avenue for creating defined polymer networks, spatially constrained by a biological template.

Archaeal Pus10 proteins can produce both pseudouridine 54 and 55 in tRNA

Pus10, a recently identified pseudouridine (Psi) synthase, does not belong to any of the five commonly identified families of Psi synthases. Pyrococcus furiosus Pus10 has been shown to produce Psi55 in tRNAs. However, in vitro studies have identified another mechanism for tRNA Psi55 production in Archaea, which uses Cbf5 and other core proteins of the H/ACA ribonucleoprotein complex, in a guide RNA-independent manner. Pus10 homologs have been observed in nearly all sequenced archaeal genomes and in some higher eukaryotes, but not in yeast and bacteria. This coincides with the presence of Psi54 in the tRNAs of Archaea and higher eukaryotes and its absence in yeast and bacteria. No tRNA Psi54 synthase has been reported so far. Here, using recombinant Methanocaldococcus jannaschii and P. furiosus Pus10, we show that these proteins can function as synthase for both tRNA Psi54 and Psi55. The two modifications seem to occur independently. Salt concentration dependent variations in these activities of both proteins are observed. The Psi54 synthase activity of M. jannaschii protein is robust, while the same activity of P. furiosus protein is weak. Probable reasons for these differences are discussed. Furthermore, unlike bacterial TruB and yeast Pus4, archaeal Pus10 does not require a U54 x A58 reverse Hoogstein base pair and pyrimidine at position 56 to convert tRNA U55 to Psi55. The homology of eukaryal Pus10 with archaeal Pus10 suggests that the former may also have a tRNA Psi54 synthase activity.

Atomic structure of the KEOPS complex: an ancient protein kinase-containing molecular machine

Kae1 is a universally conserved ATPase and part of the essential gene set in bacteria. In archaea and eukaryotes, Kae1 is embedded within the protein kinase-containing KEOPS complex. Mutation of KEOPS subunits in yeast leads to striking telomere and transcription defects, but the exact biochemical function of KEOPS is not known. As a first step to elucidating its function, we solved the atomic structure of archaea-derived KEOPS complexes involving Kae1, Bud32, Pcc1, and Cgi121 subunits. Our studies suggest that Kae1 is regulated at two levels by the primordial protein kinase Bud32, which is itself regulated by Cgi121. Moreover, Pcc1 appears to function as a dimerization module, perhaps suggesting that KEOPS may be a processive molecular machine. Lastly, as Bud32 lacks the conventional substrate-recognition infrastructure of eukaryotic protein kinases including an activation segment, Bud32 may provide a glimpse of the evolutionary history of the protein kinase family.

Methanogen homoaconitase catalyzes both hydrolyase reactions in coenzyme B biosynthesis

Homoaconitase enzymes catalyze hydrolyase reactions in the alpha-aminoadipate pathway for lysine biosynthesis or the 2-oxosuberate pathway for methanogenic coenzyme B biosynthesis. Despite the homology of this iron-sulfur protein to aconitase, previously studied homoaconitases catalyze only the hydration of cis-homoaconitate to form homoisocitrate rather than the complete isomerization of homocitrate to homoisocitrate. The MJ1003 and MJ1271 proteins from the methanogen Methanocaldococcus jannaschii formed the first homoaconitase shown to catalyze both the dehydration of (R)-homocitrate to form cis-homoaconitate, and its hydration is shown to produce homoisocitrate. This heterotetrameric enzyme also used the analogous longer chain substrates cis-(homo)(2)aconitate, cis-(homo)(3)aconitate, and cis-(homo)(4)aconitate, all with similar specificities. A combination of the homoaconitase with the M. jannaschii homoisocitrate dehydrogenase catalyzed all of the isomerization and oxidative decarboxylation reactions required to form 2-oxoadipate, 2-oxopimelate, and 2-oxosuberate, completing three iterations of the 2-oxoacid elongation pathway. Methanogenic archaeal homoaconitases and fungal homoaconitases evolved in parallel in the aconitase superfamily. The archaeal homoaconitases share a common ancestor with isopropylmalate isomerases, and both enzymes catalyzed the hydration of the minimal substrate maleate to form d-malate. The variation in substrate specificity among these enzymes correlated with the amino acid sequences of a flexible loop in the small subunits.

The 5e motif of eukaryotic signal recognition particle RNA contains a conserved adenosine for the binding of SRP72

The signal recognition particle (SRP) plays a pivotal role in transporting proteins to cell membranes. In higher eukaryotes, SRP consists of an RNA molecule and six proteins. The largest of the SRP proteins, SRP72, was found previously to bind to the SRP RNA. A fragment of human SRP72 (72c') bound effectively to human SRP RNA but only weakly to the similar SRP RNA of the archaeon Methanococcus jannaschii. Chimeras between the human and M. jannaschii SRP RNAs were constructed and used as substrates for 72c'. SRP RNA helical section 5e contained the 72c' binding site. Systematic alteration within 5e revealed that the A240G and A240C changes dramatically reduced the binding of 72c'. Human SRP RNA with a single A240G change was unable to form a complex with full-length human SRP72. Two small RNA fragments, one composed of helical section 5ef, the other of section 5e, competed equally well for the binding of 72c', demonstrating that no other regions of the SRPR RNA were required. The biochemical data completely agreed with the nucleotide conservation pattern observed across the phylogenetic spectrum. Thus, most eukaryotic SRP RNAs are likely to require for function an adenosine within their 5e motifs. The human 5ef RNA was remarkably resistant to ribonucleolytic attack suggesting that the 240-AUC-242 "loop" and its surrounding nucleotides form a peculiar compact structure recognized only by SRP72.

Protein abundance ratios for global studies of prokaryotes

The use of multidimensional capillary HPLC combined with MS/MS has allowed high qualitative and quantitative proteome coverage of prokaryotic organisms. The determination of protein abundance change between two or more conditions has matured to the point that false discovery rates can be very low and for smaller proteomes coverage is sufficiently high to explicitly consider false negative error. Selected aspects of using these methods for global protein abundance assessments are reviewed. These include instrumental issues that influence the reliability of abundance ratios; a comparison of sources of nonlinearity, errors, and data compression in proteomics and spotted cDNA arrays; strengths and weaknesses of spectral counting versus stable isotope metabolic labeling; and a survey of microbiological applications of global abundance analysis at the protein level. Proteomic results for two organisms that have been studied extensively using these methods are reviewed in greater detail. Spectral counting and metabolic labeling data are compared and the utility of proteomics for global gene regulation studies are discussed for the methanogenic Archaeon Methanococcus maripaludis. The oral pathogen Porphyromonas gingivalis is discussed as an example of an organism where a large percentage of the proteome differs in relative abundance between the intracellular and extracellular phenotype.

Differential roles of archaeal box H/ACA proteins in guide RNA-dependent and independent pseudouridine formation

RNA-guided pseudouridine (Psi) synthesis in Archaea and Eukarya requires a four-protein one-RNA containing box H/ACA ribonucleoprotein (RNP) complex. The proteins in the archaeal RNP are aCbf5, aNop10, aGar1 and L7Ae. Pyrococcus aCbf5-aNop10 is suggested to be the minimal catalytic core in this synthesis and the activity is enhanced by L7Ae and aGar1. The protein aCbf5 is homologous to eukaryal Cbf5 (dyskerin, NAP57) as well as to bacterial TruB and eukaryal Pus4; the last two produce YPsi55 in tRNAs in a guide RNA-independent manner. Here, using recombinant Methanocaldococcus jannaschii proteins, we report that aCbf5 and aGar1 together can function as a tRNA Psi55 synthase in a guide RNA-independent manner. This activity is enhanced by aNop10, but not by L7Ae. The aCbf5 alone can also produce Psi55 in tRNAs that contain the canonical 3'-CCA sequence and this activity is stimulated by aGar1. These results suggest that the roles of accessory proteins are different in guide RNA-dependent and independent Psi synthesis by aCbf5. The presence of conserved C (or U) and A at tRNA positions 56 and 58, respectively, which are required for TruB/Pus4 activity, is not essential for aCbf5-mediated Psi55 formation. Conserved A58 in tRNA normally forms a tertiary reverse Hoogstein base pair with an equally conserved U54. This base pair is recognized by TruB. Apparently aCbf5 does not require this base pair to recognize U55 for conversion to Psi55.

Interaction of signal-recognition particle 54 GTPase domain and signal-recognition particle RNA in the free signal-recognition particle

The signal-recognition particle (SRP) is a ubiquitous protein-RNA complex that targets proteins to cellular membranes for insertion or secretion. A key player in SRP-mediated protein targeting is the evolutionarily conserved core consisting of the SRP RNA and the multidomain protein SRP54. Communication between the SRP54 domains is critical for SRP function, where signal sequence binding at the M domain directs receptor binding at the GTPase domain (NG domain). These SRP activities are linked to domain rearrangements, for which the role of SRP RNA is not clear. In free SRP, a direct interaction of the GTPase domain with SRP RNA has been proposed but has never been structurally verified. In this study, we present the crystal structure at 2.5-A resolution of the SRP54-SRP19-SRP RNA complex of Methanococcus jannaschii SRP. The structure reveals an RNA-bound conformation of the SRP54 GTPase domain, in which the domain is spatially well separated from the signal peptide binding site. The association of both the N and G domains with SRP RNA in free SRP provides further structural evidence for the pivotal role of SRP RNA in the regulation of the SRP54 activity.

The biochemistry of methane oxidation

Methanotrophic bacteria oxidize methane to methanol in the first step of their metabolic pathway. Two forms of methane monooxygenase (MMO) enzymes catalyze this reaction: soluble MMO (sMMO) and membrane-bound or particulate MMO (pMMO). pMMO is expressed when copper is available, and its active site is believed to contain copper. Whereas sMMO is well characterized, most aspects of pMMO biochemistry remain unknown and somewhat controversial. This review emphasizes advances in the past two to three years related to pMMO and to copper uptake and copper-dependent regulation in methanotrophs. The pMMO metal centers have been characterized spectroscopically, and the first pMMO crystal structure has been determined. Significant effort has been devoted to improving in vitro pMMO activity. Proteins involved in sMMO regulation and additional copper-regulated proteins have been identified, and the Methylococcus capsulatus (Bath) genome has been sequenced. Finally, methanobactin (mb), a small copper chelator proposed to facilitate copper uptake, has been characterized.

Genetic screen for regulatory mutations in Methanococcus maripaludis and its use in identification of induction-deficient mutants of the euryarchaeal repressor NrpR

NrpR is an euryarchaeal transcriptional repressor of nitrogen assimilation genes. Previous studies with Methanococcus maripaludis demonstrated that NrpR binds to palindromic operator sequences, blocking transcription initiation. The metabolite 2-oxoglutarate, an indicator of cellular nitrogen deficiency, induces transcription by lowering the affinity of NrpR for operator DNA. In this report we build on existing genetic tools for M. maripaludis to develop a screen for change-of-function mutations in a transcriptional regulator and demonstrate the use of an X-Gal (5-bromo-4-chloro-3-indolyl-beta-d-galactopyranoside) screen for strict anaerobes. We use the approach to address the primary structural requirements for the response of NrpR to 2-oxoglutarate. nrpR genes from the mesophilic M. maripaludis and the hyperthermophilic Methanopyrus kandleri were targeted for mutagenesis. M. maripaludis nrpR encodes a protein with two homologous NrpR domains while the M. kandleri nrpR homolog encodes a single NrpR domain. Random point mutagenesis and alanine replacement mutagenesis identified two amino acid residues of M. kandleri NrpR involved in induction of gene expression under nitrogen-deficient conditions and thus in the response to 2-oxoglutarate. Mutagenesis of the corresponding regions in either domain of M. maripaludis NrpR resulted in a similar effect, demonstrating a conserved structure-function relationship between the two repressors. The results indicate that in M. maripaludis, both NrpR domains participate in the 2-oxoglutarate response. The approach used here has wide adaptability to other regulatory systems in methanogenic Archaea and other strict anaerobes.

The plug domain of the SecY protein stabilizes the closed state of the translocation channel and maintains a membrane seal

Proteins are translocated across membranes through a channel that is formed by the prokaryotic SecY or eukaryotic Sec61 complex. The crystal structure of the SecY channel from M. jannaschii revealed a plug domain that appears to seal the channel in its closed state. However, the role of the plug remains unclear, particularly because plug deletion mutants in S. cerevisiae are functional. Here, we demonstrate that plug deletion mutants in E. coli SecY are also functional and even efficiently translocate proteins with defective or missing signal sequences. The crystal structures of equivalent plug deletions in SecY of M. jannaschii show that, although the overall structures are maintained, new plugs are formed. These lack many interactions that normally stabilize the closed channel, explaining why the channels can open for proteins with signal-sequence mutations. Our data show that the plug domain is required to maintain a closed state of the channel and suggest a mechanism for channel gating.

[Crystal structures of mutant ribosomal proteins L1]

Nine mutant forms of ribosomal proteins L1 from the bacterium Thermus thermophilus and the archaeon Methanococcus jannaschii were obtained. Their crystal structures were determined and analyzed. Earlier determined structure of S179C TthL1 was also thoroughly analyzed. Five from ten mutant proteins reveal essential changes of spatial structure caused by surface point mutation. It proves that for correct studies of biological processes by site-directed mutagenesis it is necessary to determine or at least to model spatial structures of mutant proteins. Detailed comparison of mutant L1 structures with that of corresponding wild type proteins reveals that side chain of a mutated amino acid residue tries to locate like the side chain of the original residue in the wild type protein. This observation helps to model the mutant structures.

ATP-induced structural transitions in PAN, the proteasome-regulatory ATPase complex in Archaea

ATP binding to the PAN-ATPase complex in Archaea or the homologous 19 S protease-regulatory complex in eukaryotes induces association with the 20 S proteasome and opening of its substrate entry channel, whereas ATP hydrolysis allows unfolding of globular substrates. To clarify the conformational changes associated with ATP binding and hydrolysis, we used protease sensitivity to monitor the conformations of the PAN ATPase from Methanococcus jannischii. Exhaustive trypsin treatment of PAN generated five distinct fragments, two of which differed when a nucleotide (either ATP, ATP gamma S, or ADP) was bound. Surprisingly, the nucleotide concentrations altering protease sensitivity were much lower (K(a) 20-40 microm) than are required for ATP-dependent protein breakdown by the PAN-20S proteasome complex (K(m) approximately 300-500 microm). Unlike trypsin, proteinase K yielded several fragments that differed in the ATP gamma S and ADP-bound forms, and thus revealed conformational transitions associated with ATP hydrolysis. Mapping the fragments generated by each revealed that nucleotide binding and hydrolysis induce local conformational changes, affecting the Walker A and B nucleotide-binding motif, as well as global changes extending to its carboxyl terminus. The location and overlap of the fragments also suggest that the conformation of the six subunits is not identical, probably because they do not all bind ATP simultaneously. Partial nucleotide occupancy was supported by direct assays, which demonstrated that, at saturating conditions, only four nucleotides are bound to hexameric PAN. Using the protease protection maps, we modeled the conformational changes associated with ATP binding and hydrolysis in PAN based on the x-ray structures of the homologous AAA ATPase, HslU.

Diverse homologues of the archaeal repressor NrpR function similarly in nitrogen regulation

NrpR is a transcriptional repressor of nitrogen assimilation genes that was recently discovered and characterized in the methanogenic archaeon Methanococcus maripaludis. NrpR homologues are widely distributed in Euryarchaeota and present in a few bacterial species. They exist in three different domain configurations: a single ORF encoding one NrpR domain following an N-terminal helix-turn-helix (HTH); a single ORF encoding two NrpR domains fused in tandem following an N-terminal HTH; and two separate ORFs, one with a single domain following an N-terminal HTH and one with a single domain without a HTH. Phylogenetic analysis indicated that the NrpR family forms five distinct groups: the single domain HTH type, the two domains of the double domain HTH type and the two separately encoded domains. To determine the function of diverse NrpR homologues, representative genes in were expressed an Methanococcus maripaludis nrpR deletion mutant. Homologues from species that possess a single gene restored regulated repression, regardless of domain structure. In the case of Methanosarcina acetivorans that contains two genes, both were required. The results show that distantly related NrpR homologues that are present in widely dispersed phyla regulate the expression of nitrogen assimilation genes in a similar fashion.

Functionally distinct genes regulated by hydrogen limitation and growth rate in methanogenic Archaea

The use of molecular hydrogen as electron donor for energy generation is a defining characteristic of the hydrogenotrophic methanogens, an ancient group that dominates the phylum Euryarchaeota. We present here a global study of changes in mRNA abundance in response to hydrogen availability for a hydrogenotrophic methanogen. Cells of Methanococcus maripaludis were grown by using continuous culture to deconvolute the effects of hydrogen limitation and growth rate, and microarray analyses were conducted. Hydrogen limitation markedly increased mRNA levels for genes encoding enzymes of the methanogenic pathway that reduce or oxidize the electron-carrying deazaflavin, coenzyme F(420). F(420)-dependent redox functions in energy-generating metabolism are characteristic of the methanogenic Archaea, and the results show that their regulation is distinct from other redox processes in the cell. Rapid growth increased mRNA levels of the gene for an unusual hydrogenase, the hydrogen-dependent methylenetetrahydromethanopterin dehydrogenase.

Nucleotide-dependent allostery within the ABC transporter ATP-binding cassette: a computational study of the MJ0796 dimer

ATP-binding cassette transporters perform energy-dependent transmembrane solute trafficking in all organisms. These proteins often mediate cellular resistance to therapeutic drugs and are involved in a range of human genetic diseases. Enzymological studies have implicated a helical subdomain within the ATP-binding cassette nucleotide-binding domain in coupling ATP hydrolysis to solute transport in the transmembrane domains. Consistent with this, structural and computational analyses have indicated that the helical subdomain undergoes nucleotide-dependent movement relative to the core of the nucleotide-binding domain fold. Here we use theoretical methods to examine the allosteric nucleotide dependence of helical subdomain transitions to further elucidate its role in interactions between the transmembrane and nucleotide-binding domains. Unrestrained 30-ns molecular dynamics simulations of the ATP-bound, ADP-bound, and apo states of the MJ0796 monomer support the idea that interaction of a conserved glutamine residue with the catalytic metal mediates the rotation of the helical subdomain in response to nucleotide binding and hydrolysis. Simulations of the nucleotide-binding domain dimer revealed that ATP hydrolysis induces a large transition of one helical subdomain, resulting in an asymmetric conformation of the dimer not observed previously. A coarse-grained elastic network analysis supports this finding, revealing the existence of corresponding dynamic modes intrinsic to the contact topology of the protein. The implications of these findings for the coupling of ATP hydrolysis to conformational changes in the transmembrane domains required for solute transport are discussed in light of recent whole transporter structures.

A dynamic programming approach for finding common patterns in RNAs

We developed a dynamic programming approach of computing common sequence structure patterns among two RNAs given their primary sequences and their secondary structures. Common patterns between two RNAs are defined to share the same local sequential and structural properties. The locality is based on the connections of nucleotides given by their phosphodiester and hydrogen bonds. The idea of interpreting secondary structures as chains of structure elements leads us to develop an efficient dynamic programming approach in time O(nm) and space O(nm), where n and m are the lengths of the RNAs. The biological motivation is given by detecting common, local regions of RNAs, although they do not necessarily share global sequential and structural properties. This might happen if RNAs fold into different structures but share a lot of local, stable regions. Here, we illustrate our algorithm on Hepatitis C virus internal ribosome entry sites. Our method is useful for detecting and describing local motifs as well. An implementation in C++ is available and can be obtained by contacting one of the authors.

Metabolic modeling of a mutualistic microbial community

The rate of production of methane in many environments depends upon mutualistic interactions between sulfate-reducing bacteria and methanogens. To enhance our understanding of these relationships, we took advantage of the fully sequenced genomes of Desulfovibrio vulgaris and Methanococcus maripaludis to produce and analyze the first multispecies stoichiometric metabolic model. Model results were compared to data on growth of the co-culture on lactate in the absence of sulfate. The model accurately predicted several ecologically relevant characteristics, including the flux of metabolites and the ratio of D. vulgaris to M. maripaludis cells during growth. In addition, the model and our data suggested that it was possible to eliminate formate as an interspecies electron shuttle, but hydrogen transfer was essential for syntrophic growth. Our work demonstrated that reconstructed metabolic networks and stoichiometric models can serve not only to predict metabolic fluxes and growth phenotypes of single organisms, but also to capture growth parameters and community composition of simple bacterial communities.

Detecting molecular interactions that stabilize, activate and guide ligand-binding of the sodium/proton antiporter MjNhaP1 from Methanococcus jannaschii

Integral membrane proteins are involved in virtually every cellular process. Precisely regulating these machineries would allow controlling many human and vertebrate diseases. Embedded into cellular membranes, membrane proteins establish molecular interactions that sensitively react to environmental changes and to molecular compounds, such as ligands or inhibitors. We applied atomic force microscopy (AFM) to image the Na(+)/H(+) antiporter MjNhaP1 from Methanococcus jannaschii, and single-molecule force spectroscopy (SMFS) to probe molecular interactions that drive the protein structure-function relationship. High-resolution AFM topographs showed the dimeric assembly of MjNhaP1 being reconstituted into a lipid bilayer. SMFS of MjNhaP1 unraveled molecular interactions stabilizing individual structural domains. Transmembrane domains exhibited certain probabilities to unfold individually or cooperatively with other domains resulting in different unfolding pathways. Helices VIII and X established pH sensitive interactions altering significantly upon MjNhaP1 activation, while removal of the ligand (Na(+)) destabilized the entire antiporter except helix VIII. It is assumed that Asp234/235 of helix VIII are involved in the ligand-binding site and that helix X plays a functional role in the activation of the transporter.

RNA-dependent conversion of phosphoserine forms selenocysteine in eukaryotes and archaea

The trace element selenium is found in proteins as selenocysteine (Sec), the 21st amino acid to participate in ribosome-mediated translation. The substrate for ribosomal protein synthesis is selenocysteinyl-tRNA(Sec). Its biosynthesis from seryl-tRNA(Sec) has been established for bacteria, but the mechanism of conversion from Ser-tRNA(Sec) remained unresolved for archaea and eukarya. Here, we provide evidence for a different route present in these domains of life that requires the tRNA(Sec)-dependent conversion of O-phosphoserine (Sep) to Sec. In this two-step pathway, O-phosphoseryl-tRNA(Sec) kinase (PSTK) converts Ser-tRNA(Sec) to Sep-tRNA(Sec). This misacylated tRNA is the obligatory precursor for a Sep-tRNA:Sec-tRNA synthase (SepSecS); this protein was previously annotated as SLA/LP. The human and archaeal SepSecS genes complement in vivo an Escherichia coli Sec synthase (SelA) deletion strain. Furthermore, purified recombinant SepSecS converts Sep-tRNA(Sec) into Sec-tRNA(Sec) in vitro in the presence of sodium selenite and purified recombinant E. coli selenophosphate synthetase (SelD). Phylogenetic arguments suggest that Sec decoding was present in the last universal common ancestor. SepSecS and PSTK coevolved with the archaeal and eukaryotic lineages, but the history of PSTK is marked by several horizontal gene transfer events, including transfer to non-Sec-decoding Cyanobacteria and fungi.

Biochemical and genetic characterization of an early step in a novel pathway for the biosynthesis of aromatic amino acids and p-aminobenzoic acid in the archaeon Methanococcus maripaludis

Methanococcus maripaludis is a strictly anaerobic, methane-producing archaeon and facultative autotroph capable of biosynthesizing all the amino acids and vitamins required for growth. In this work, the novel 6-deoxy-5-ketofructose-1-phosphate (DKFP) pathway for the biosynthesis of aromatic amino acids (AroAAs) and p-aminobenzoic acid (PABA) was demonstrated in M. maripaludis. Moreover, PABA was shown to be derived from an early intermediate in AroAA biosynthesis and not from chorismate. Following metabolic labelling with [U-(13)C]-acetate, the expected enrichments for phenylalanine and arylamine derived from PABA were observed. DKFP pathway activity was reduced following growth with aryl acids, an alternative source of the AroAAs. Lastly, a deletion mutant of aroA', which encodes the first step in the DKFP pathway, required AroAAs and PABA for growth. Complementation of the mutants by an aroA' expression vector restored the wild-type phenotype. In contrast, a deletion of aroB', which encodes the second step in the DKFP pathway, did not require AroAAs or PABA for growth. Presumably, methanococci contain an alternative activity for this step. These results identify the initial reactions of a new pathway for the biosynthesis of PABA in methanococci.

Calcium stiffens archaeal Rad51 recombinase from Methanococcus voltae for homologous recombination

Archaeal RadA or Rad51 recombinases are close homologues of eukaryal Rad51 and DMC1. These and bacterial RecA orthologues play a key role in DNA repair by forming helical nucleoprotein filaments in which a hallmark strand exchange reaction between homologous DNA substrates occurs. Recent studies have discovered the stimulatory role by calcium on human and yeast recombinases. Here we report that the strand exchange activity but not the ATPase activity of an archaeal RadA/Rad51 recombinase from Methanococcus voltae (MvRadA) is also subject to calcium stimulation. Crystallized MvRadA filaments in the presence of CaCl(2) resemble that of the recently reported ATPase active form in the presence of an activating dose of KCl. At the ATPase center, one Ca(2+) ion takes the place of two K(+) ions in the K(+)-bound form. The terminal phosphate of the nonhydrolyzable ATP analogue is in a staggered conformation in the Ca(2+)-bound form. In comparison, an eclipsed conformation was seen in the K(+)-bound form. Despite the changes in the ATPase center, both forms harbor largely ordered L2 regions in essentially identical conformations. These data suggest a unified stimulation mechanism by potassium and calcium because of the existence of a conserved ATPase center promiscuous in binding cations.

Simulations of a Protein Translocation Pore: SecY

SecY is the central channel protein of the SecYEbeta translocon, the structure of which has been determined by X-ray diffraction. Extended (15 ns) MD simulations of the isolated SecY protein in a phospholipid bilayer have been performed to explore the relationship between protein flexibility and the mechanisms of channel gating. In particular, principal components analysis of the simulation trajectory has been used to probe the intrinsic flexibility of the isolated SecY protein in the absence of the gamma-subunit (SecE) clamp. Analysis and visualization of the principal eigenvectors support a "plug and clamshell" model of SecY channel gating. The simulation results also indicate that hydrophobic gating at the central pore ring prevents leakage of water and ions through the channel in the absence of a translocating peptide.

Identification of genes involved in the biosynthesis and attachment of Methanococcus voltae N-linked glycans: insight into N-linked glycosylation pathways in Archaea

N-linked glycosylation is recognized as an important post-translational modification across all three domains of life. However, the understanding of the genetic pathways for the assembly and attachment of N-linked glycans in eukaryotic and bacterial systems far outweighs the knowledge of comparable processes in Archaea. The recent characterization of a novel trisaccharide [beta-ManpNAcA6Thr-(1-4)-beta-GlcpNAc3NAcA-(1-3)-beta-GlcpNAc]N-linked to asparagine residues in Methanococcus voltae flagellin and S-layer proteins affords new opportunities to investigate N-linked glycosylation pathways in Archaea. In this contribution, the insertional inactivation of several candidate genes within the M. voltae genome and their resulting effects on flagellin and S-layer glycosylation are reported. Two of the candidate genes were shown to have effects on flagellin and S-layer protein molecular mass and N-linked glycan structure. Further examination revealed inactivation of either of these two genes also had effects on flagella assembly. These genes, designated agl (archaeal glycosylation) genes, include a glycosyl transferase (aglA) involved in the attachment of the terminal sugar to the glycan and an STT3 oligosaccharyl transferase homologue (aglB) involved in the transfer of the complete glycan to the flagellin and S-layer proteins. These findings document the first experimental evidence for genes involved in any glycosylation process within the domain Archaea.