Change of carbon source causes dramatic effects in the phospho-proteome of the archaeon Sulfolobus solfataricus

Novel Tyrosine Kinase-Mediated Phosphorylation With Dual Specificity Plays a Key Role in the Modulation of Streptococcus pyogenes Physiology and Virulence

Streptococcus pyogenes (Group A Streptococcus, GAS) genomes do not contain a gene encoding a typical bacterial-type tyrosine kinase (BY-kinase) but contain an orphan gene-encoding protein Tyr-phosphatase (SP-PTP). Hence, the importance of Tyr-phosphorylation is underappreciated and not recognized for its role in GAS pathophysiology and pathogenesis. The fact that SP-PTP dephosphorylates Abl-tyrosine kinase-phosphorylated myelin basic protein (MBP), and SP-STK (S. pyogenes Ser/Thr kinase) also autophosphorylates its Tyr¹⁰¹-residue prompted us to identify a putative tyrosine kinase and Tyr-phosphorylation in GAS. Upon a genome-wide search of kinases possessing a classical Walker motif, we identified a non-canonical tyrosine kinase M5005_Spy_1476, a ∼17 kDa protein (153 aa) (SP-TyK). The purified recombinant SP-TyK autophosphorylated in the presence of ATP. In vitro and in vivo phosphoproteomic analyses revealed two key phosphorylated tyrosine residues located within the catalytic domain of SP-TyK. An isogenic mutant lacking SP-TyK derived from the M1T1 strain showed a retarded growth pattern. It displayed defective cell division and long chains with multiple parallel septa, often resulting in aggregates. Transcriptomic analysis of the mutant revealed 287 differentially expressed genes responsible for GAS pathophysiology and pathogenesis. SP-TyK also phosphorylated GAS CovR, WalR, SP-STP, and SDH/GAPDH proteins with dual specificity targeting their Tyr/Ser/Thr residues as revealed by biochemical and mass-spectrometric-based phosphoproteomic analyses. SP-TyK-phosphorylated CovR bound to PcovR efficiently. The mutant displayed sustained release of IL-6 compared to TNF-α during co-culturing with A549 lung cell lines, attenuation in mice sepsis model, and significantly reduced ability to adhere to and invade A549 lung cells and form biofilms on abiotic surfaces. SP-TyK, thus, plays a critical role in fine-tuning the regulation of key cellular functions essential for GAS pathophysiology and pathogenesis through post-translational modifications and hence, may serve as a promising target for future therapeutic developments.

The biology of thermoacidophilic archaea from the order Sulfolobales

Thermoacidophilic archaea belonging to the order Sulfolobales thrive in extreme biotopes, such as sulfuric hot springs and ore deposits. These microorganisms have been model systems for understanding life in extreme environments, as well as for probing the evolution of both molecular genetic processes and central metabolic pathways. Thermoacidophiles, such as the Sulfolobales, use typical microbial responses to persist in hot acid (e.g. motility, stress response, biofilm formation), albeit with some unusual twists. They also exhibit unique physiological features, including iron and sulfur chemolithoautotrophy, that differentiate them from much of the microbial world. Although first discovered >50 years ago, it was not until recently that genome sequence data and facile genetic tools have been developed for species in the Sulfolobales. These advances have not only opened up ways to further probe novel features of these microbes but also paved the way for their potential biotechnological applications. Discussed here are the nuances of the thermoacidophilic lifestyle of the Sulfolobales, including their evolutionary placement, cell biology, survival strategies, genetic tools, metabolic processes and physiological attributes together with how these characteristics make thermoacidophiles ideal platforms for specialized industrial processes.

Analysis of the phosphorylome of trichoderma reesei cultivated on sugarcane bagasse suggests post-translational regulation of the secreted glycosyl hydrolase Cel7A

•

Phosphorylome of Trichoderma reesei reveals phosphosites in some glycosyl hydrolases.

•

Phosphoserine and phosphothreonine is the major phosphosites identified.

•

Protein Kinase C is the most frequently predicted kinase in phosphorylome.

•

The cellulase Cel7A activity is affected by dephosphorylation.

Trichoderma reesei is one of the major producers of holocellulases. It is known that in T. reesei, protein production patterns can change in a carbon source-dependent manner. Here, we performed a phosphorylome analysis of T. reesei grown in the presence of sugarcane bagasse and glucose as carbon source. In presence of sugarcane bagasse, a total of 114 phosphorylated proteins were identified. Phosphoserine and phosphothreonine corresponded to 89.6% of the phosphosites and 10.4% were related to phosphotyrosine. Among the identified proteins, 65% were singly phosphorylated, 19% were doubly phosphorylated, 12% were triply phosphorylated, and 4% displayed even higher phosphorylation. Seventy-five kinases were predicted to phosphorylate the sites identified in this work, and the most frequently predicted serine/threonine kinase was PKC1. Among phosphorylated proteins, four glycosyl hydrolases were predicted to be secreted. Interestingly, Cel7A activity, the most secreted protein, was reduced to approximately 60% after in vitro dephosphorylation, suggesting that phosphorylation might alter Cel7A structure, substrate affinity, and targeting of the substrate to its carbohydrate-binding domain. These results suggest a novel post-translational regulation of Cel7A.

Comparative Phosphoproteomics of Classical Bordetellae Elucidates the Potential Role of Serine, Threonine and Tyrosine Phosphorylation in Bordetella Biology and Virulence

The Bordetella genus is divided into two groups: classical and non-classical. Bordetella pertussis, Bordetella bronchiseptica and Bordetella parapertussis are known as classical bordetellae, a group of important human pathogens causing whooping cough or whooping cough-like disease and hypothesized to have evolved from environmental non-classical bordetellae. Bordetella infections have increased globally driving the need to better understand these pathogens for the development of new treatments and vaccines. One unexplored component in Bordetella is the role of serine, threonine and tyrosine phosphorylation. Therefore, this study characterized the phosphoproteome of classical bordetellae and examined its potential role in Bordetella biology and virulence. Applying strict identification of localization criteria, this study identified 70 unique phosphorylated proteins in the classical bordetellae group with a high degree of conservation. Phosphorylation was a key regulator of Bordetella metabolism with proteins involved in gluconeogenesis, TCA cycle, amino acid and nucleotide synthesis significantly enriched. Three key virulence pathways were also phosphorylated including type III secretion system, alcaligin synthesis and the BvgAS master transcriptional regulatory system for virulence genes in Bordetella. Seven new phosphosites were identified in BvgA with 6 located in the DNA binding domain. Of the 7, 4 were not present in non-classical bordetellae. This suggests that serine/threonine phosphorylation may play an important role in stabilizing/destabilizing BvgA binding to DNA for fine-tuning of virulence gene expression and that BvgA phosphorylation may be an important factor separating classical from non-classical bordetellae. This study provides the first insight into the phosphoproteome of classical Bordetella species and the role that Ser/Thr/Tyr phosphorylation may play in Bordetella biology and virulence.

Wzb of Vibrio vulnificus represents a new group of low-molecular-weight protein tyrosine phosphatases with a unique insertion in the W-loop

Protein tyrosine phosphorylation regulates the production of capsular polysaccharide, an essential virulence factor of the deadly pathogen Vibrio vulnificus. The process requires the protein tyrosine kinase Wzc and its cognate phosphatase Wzb, both of which are largely uncharacterized. Herein, we report the structures of Wzb of V. vulnificus (VvWzb) in free and ligand-bound forms. VvWzb belongs to the low-molecular-weight protein tyrosine phosphatase (LMWPTP) family. Interestingly, it contains an extra four-residue insertion in the W-loop, distinct from all known LMWPTPs. The W-loop of VvWzb protrudes from the protein body in the free structure, but undergoes significant conformational changes to fold toward the active site upon ligand binding. Deleting the four-residue insertion from the W-loop severely impaired the enzymatic activity of VvWzb, indicating its importance for optimal catalysis. However, mutating individual residues or even substituting the whole insertion with four alanine residues only modestly decreased the enzymatic activity, suggesting that the contribution of the insertion to catalysis is not determined by the sequence specificity. Furthermore, inserting the four residues into Escherichia coli Wzb at the corresponding position enhanced its activity as well, indicating that the four-residue insertion in the W-loop can act as a general activity enhancing element for other LMWPTPs. The novel W-loop type and phylogenetic analysis suggested that VvWzb and its homologs should be classified into a new group of LMWPTPs. Our study sheds new insight into the catalytic mechanism and structural diversity of the LMWPTP family and promotes the understanding of the protein tyrosine phosphorylation system in prokaryotes.

Mitochondrial protein phosphorylation in yeast revisited

Protein phosphorylation is one of the best-known post-translational modifications occurring in all domains of life. In eukaryotes, protein phosphorylation affects all cellular compartments including mitochondria. High-throughput techniques of mass spectrometry combined with cell fractionation and biochemical methods yielded thousands of phospho-sites on hundreds of mitochondrial proteins. We have compiled the information on mitochondrial protein kinases and phosphatases and their substrates in Saccharomyces cerevisiae and provide the current state-of-the-art overview of mitochondrial protein phosphorylation in this model eukaryote. Using several examples, we describe emerging features of the yeast mitochondrial phosphoproteome and present challenges lying ahead in this exciting field.

Origin of the Phosphoprotein Phosphatase (PPP) sequence family in Bacteria: Critical ancestral sequence changes, radiation patterns and substrate binding features

Background

Phosphoprotein phosphatases (PPP) belong to the PPP Sequence family, which in turn belongs to the broader metallophosphoesterase (MPE) superfamily. The relationship between the PPP Sequence family and other members of the MPE superfamily remains unresolved, in particular what transitions took place in an ancestral MPE to ultimately produce the phosphoprotein specific phosphatases (PPPs).

Methods

We use structural and sequence alignment data, phylogenetic tree analysis, sequence signature (Weblogo) analysis, in silico protein-peptide modeling data, and in silico mutagenesis to trace a likely route of evolution from MPEs to the PPP Sequence family. Hidden Markov Model (HMM) based iterative database search strategies were utilized to identify PPP Sequence Family members from numerous bacterial groups.

Results

Using Mre11 as proxy for an ancestral nuclease-like MPE we trace a possible evolutionary route that alters a single active site substrate binding His-residue to yield a new substrate binding accessory, the "2-Arg-Clamp". The 2-Arg-Clamp is not found in MPEs, but is present in all PPP Sequence family members, where the phosphomonesterase reaction predominates. Variation in position of the clamp arginines and a supplemental sequence loop likely provide substrate specificity for each PPP Sequence family group.

Conclusions

Loss of a key substrate binding His-in MPEs opened the path to bind novel substrates and evolution of the 2-Arg-Clamp, a sequence change seen in both bacterial and eukaryotic phosphoprotein phosphatases.General significance: We establish a likely evolutionary route from nuclease-like MPE to PPP Sequence family enzymes, that includes the phosphoprotein phosphatases.

Phosphoproteomic Analysis Reveals Rio1-Related Protein Phosphorylation Changes in Response to UV Irradiation in Sulfolobus islandicus REY15A

DNA damage response (DDR) in eukaryotes is largely regulated by protein phosphorylation. In archaea, many proteins are phosphorylated, however, it is unclear how the cells respond to DNA damage through global protein phosphorylation. We previously found that Δrio1, a Rio1 kinase homolog deletion strain of Sulfolobus islandicus REY15A, was sensitive to UV irradiation. In this study, we showed that Δrio1 grew faster than the wild type. Quantitative phosphoproteomic analysis of the wild type and Δrio1, untreated and irradiated with UV irradiation, revealed 562 phosphorylated sites (with a Ser/Thr/Tyr ratio of 65.3%/23.8%/10.9%) of 333 proteins in total. The phosphorylation levels of 35 sites of 30 proteins changed with >1.3-fold in the wild type strain upon UV irradiation. Interestingly, more than half of the UV-induced changes in the wild type did not occur in the Δrio1 strain, which were mainly associated with proteins synthesis and turnover. In addition, a protein kinase and several transcriptional regulators were differentially phosphorylated after UV treatment, and some of the changes were dependent on Rio1. Finally, many proteins involved in various cellular metabolisms exhibited Riol-related and UV-independent phosphorylation changes. Our results suggest that Rio1 is involved in the regulation of protein recycling and signal transduction in response to UV irradiation, and plays regulatory roles in multiple cellular processes in S. islandicus.

The Phosphatase PP2A Interacts With ArnA and ArnB to Regulate the Oligomeric State and the Stability of the ArnA/B Complex

In the crenarchaeon Sulfolobus acidocaldarius, the archaellum, a type-IV pilus like motility structure, is synthesized in response to nutrient starvation. Synthesis of components of the archaellum is controlled by the archaellum regulatory network (arn). Protein phosphorylation plays an important role in this regulatory network since the deletion of several genes encoding protein kinases and the phosphatase PP2A affected cell motility. Several proteins in the archaellum regulatory network can be phosphorylated, however, details of how phosphorylation levels of different components affect archaellum synthesis are still unknown. To identify proteins interacting with the S. acidocaldarius phosphatases PTP and PP2A, co-immunoprecipitation assays coupled to mass spectrometry analysis were performed. Thirty minutes after growth in nutrient starvation medium, especially a conserved putative ATP/GTP binding protein (Saci_1281), a universal stress protein (Saci_0887) and the archaellum regulators ArnA and ArnB were identified as highly abundant interaction proteins of PP2A. The interaction between ArnA, ArnB, and PP2A was further studied. Previous studies showed that the Forkhead-associated domain containing ArnA interacts with von Willebrand type A domain containing ArnB, and that both proteins could be phosphorylated by the kinase ArnC in vitro. The ArnA/B heterodimer was reconstituted from the purified proteins. In complex with ArnA, phosphorylation of ArnB by the ArnC kinase was strongly stimulated and resulted in formation of (ArnA/B)₂ and higher oligomeric complexes, while association and dephosphorylation by PP2A resulted in dissociation of these ArnA/B complexes.

Post-Translational Modifications Aid Archaeal Survival

Since the pioneering work of Carl Woese, Archaea have fascinated biologists of almost all areas given their unique evolutionary status, wide distribution, high diversity, and ability to grow in special environments. Archaea often thrive in extreme conditions such as high temperature, high/low pH, high salinity, and anoxic ecosystems. All of these are threats to the stability and proper functioning of biological molecules, especially proteins and nucleic acids. Post-translational modifications (PTMs), such as phosphorylation, methylation, acetylation, and glycosylation, are reportedly widespread in Archaea and represent a critical adaptive mechanism to extreme habitats. Here, we summarize our current understanding of the contributions of PTMs to aid in extremophile survival, with a particular focus on the maintenance of genome stability.

Tyrosine Phosphorylation as a Widespread Regulatory Mechanism in Prokaryotes

Phosphorylation events modify bacterial and archaeal proteomes, imparting cells with rapid and reversible responses to specific environmental stimuli or niches. Phosphorylated proteins are generally modified at one or more serine, threonine, or tyrosine residues. Within the last ten years, increasing numbers of global phosphoproteomic surveys of prokaryote species have revealed an abundance of tyrosine-phosphorylated proteins. In some cases, novel phosphorylation-dependent regulatory paradigms for cell division, gene transcription, and protein translation have been identified, suggesting that a wide scope of prokaryotic physiology remains to be characterized. Recent observations of bacterial proteins with putative phosphotyrosine binding pockets or Src homology 2 (SH2)-like domains suggest the presence of phosphotyrosine-dependent protein interaction networks. Here in this minireview, we focus on protein tyrosine phosphorylation, a posttranslational modification once thought to be rare in prokaryotes but which has emerged as an important regulatory facet in microbial biology.

Structure and interactions of the archaeal motility repression module ArnA–ArnB that modulates archaellum gene expression in Sulfolobus acidocaldarius

Phosphorylation-dependent interactions play crucial regulatory roles in all domains of life. Forkhead-associated (FHA) and von Willebrand type A (vWA) domains are involved in several phosphorylation-dependent processes of multiprotein complex assemblies. Although well-studied in eukaryotes and bacteria, the structural and functional contexts of these domains are not yet understood in Archaea. Here, we report the structural base for such an interacting pair of FHA and vWA domain-containing proteins, ArnA and ArnB, in the thermoacidophilic archaeon Sulfolobus acidocaldarius, where they act synergistically and negatively modulate motility. The structure of the FHA domain of ArnA at 1.75 Å resolution revealed that it belongs to the subclass of FHA domains, which recognizes double-pSer/pThr motifs. We also solved the 1.5 Å resolution crystal structure of the ArnB paralog vWA2, disclosing a complex topology comprising the vWA domain, a β-sandwich fold, and a C-terminal helix bundle. We further show that ArnA binds to the C terminus of ArnB, which harbors all the phosphorylation sites identified to date and is important for the function of ArnB in archaellum regulation. We also observed that expression levels of the archaellum components in response to changes in nutrient conditions are independent of changes in ArnA and ArnB levels and that a strong interaction between ArnA and ArnB observed during growth on rich medium sequentially diminishes after nutrient limitation. In summary, our findings unravel the structural features in ArnA and ArnB important for their interaction and functional archaellum expression and reveal how nutrient conditions affect this interaction.

The Carbon Switch at the Level of Pyruvate and Phosphoenolpyruvate in Sulfolobus solfataricus P2

Sulfolobus solfataricus P2 grows on different carbohydrates as well as alcohols, peptides and amino acids. Carbohydrates such as D-glucose or D-galactose are degraded via the modified, branched Entner–Doudoroff (ED) pathway whereas growth on peptides requires the Embden–Meyerhof–Parnas (EMP) pathway for gluconeogenesis. As for most hyperthermophilic Archaea an important control point is established at the level of triosephophate conversion, however, the regulation at the level of pyruvate/phosphoenolpyruvate conversion was not tackled so far. Here we describe the cloning, expression, purification and characterization of the pyruvate kinase (PK, SSO0981) and the phosphoenolpyruvate synthetase (PEPS, SSO0883) of Sul. solfataricus. The PK showed only catabolic activity [catalytic efficiency (PEP): 627.95 mM^-1s^-1, 70°C] with phosphoenolpyruvate as substrate and ADP as phosphate acceptor and was allosterically inhibited by ATP and isocitrate (K_i 0.8 mM). The PEPS was reversible, however, exhibited preferred activity in the gluconeogenic direction [catalytic efficiency (pyruvate): 1.04 mM^-1s^-1, 70°C] and showed some inhibition by AMP and α-ketoglutarate. The gene SSO2829 annotated as PEPS/pyruvate:phosphate dikinase (PPDK) revealed neither PEPS nor PPDK activity. Our studies suggest that the energy charge of the cell as well as the availability of building blocks in the citric acid cycle and the carbon/nitrogen balance plays a major role in the Sul. solfataricus carbon switch. The comparison of regulatory features of well-studied hyperthermophilic Archaea reveals a close link and sophisticated coordination between the respective sugar kinases and the kinetic and regulatory properties of the enzymes at the level of PEP-pyruvate conversion.

Two membrane-bound transcription factors regulate expression of various type-IV-pili surface structures in Sulfolobus acidocaldarius

In Archaea and Bacteria, gene expression is tightly regulated in response to environmental stimuli. In the thermoacidophilic crenarchaeon Sulfolobus acidocaldarius nutrient limitation induces expression of the archaellum, the archaeal motility structure. This expression is orchestrated by a complex hierarchical network of positive and negative regulators-the archaellum regulatory network (arn). The membrane-bound one-component system ArnR and its paralog ArnR1 were recently described as main activators of archaellum expression in S. acidocaldarius. They regulate gene expression of the archaellum operon by targeting the promoter of flaB, encoding the archaellum filament protein. Here we describe a strategy for the isolation and biochemical characterization of these two archaellum regulators. Both regulators are capable of forming oligomers and are phosphorylated by the Ser/Thr kinase ArnC. Apart from binding to pflaB, ArnR but not ArnR1 bound to promoter sequences of aapF and upsX, which encode components of the archaeal adhesive pilus and UV-inducible pili system, demonstrating a regulatory connection between different surface appendages of S. acidocaldarius.

Conserved principles of transcriptional networks controlling metabolic flexibility in archaea

Gene regulation is intimately connected with metabolism, enabling the appropriate timing and tuning of biochemical pathways to substrate availability. In microorganisms, such as archaea and bacteria, transcription factors (TFs) often directly sense external cues such as nutrient substrates, metabolic intermediates, or redox status to regulate gene expression. Intense recent interest has characterized the functions of a large number of such regulatory TFs in archaea, which regulate a diverse array of unique archaeal metabolic capabilities. However, it remains unclear how the co-ordinated activity of the interconnected metabolic and transcription networks produces the dynamic flexibility so frequently observed in archaeal cells as they respond to energy limitation and intermittent substrate availability. In this review, we communicate the current state of the art regarding these archaeal networks and their dynamic properties. We compare the topology of these archaeal networks to those known for bacteria to highlight conserved and unique aspects. We present a new computational model for an exemplar archaeal network, aiming to lay the groundwork toward understanding general principles that unify the dynamic function of integrated metabolic-transcription networks across archaea and bacteria.

Global effect of the lack of inorganic polyphosphate in the extremophilic archaeon Sulfolobus solfataricus: A proteomic approach

Inorganic polyphosphates (polyP) are present in all living cells and several important functions have been described for them. They are involved in the response to stress conditions, such as nutrient depletion, oxidative stress and toxic metals amongst others. A recombinant strain of Sulfolobus solfataricus unable to accumulate polyP was designed by the overexpression of its endogenous ppx gene. The overall impact of the lack of polyP on this S. solfataricus polyP (-) strain was analyzed by using quantitative proteomics (isotope-coded protein label, ICPL). Stress-related proteins, such as peroxiredoxins and heat shock proteins, proteins involved in metabolism and several others were produced at higher levels in the ppx expression strain. The polyP deficient strain showed an increased copper sensitivity and an earlier transcriptional up-regulation of copA gene coding for the P-type copper-exporting ATPase. This implies a complementary function of both copper resistance systems. These results strongly suggests that the lack of polyP makes this hyperthermophilic archaeon more sensitive to toxic conditions, such as an exposure to metals or other harmful stimuli, emphasizing the importance of this inorganic phosphate polymers in the adaptations to live in the environmental conditions in which thermoacidophilic archaea thrive. SIGNIFICANCE: Inorganic polyphosphate (polyP) are ubiquitous molecules with many functions in living organisms. Few studies related to these polymers have been made in archaea. The construction of a polyP deficient recombinant strain of Sulfolobus solfataricus allowed the study of the global changes in the proteome of this thermoacidophilic archaeon in the absence of polyP compared with the wild type strain. The results obtained using quantitative proteomics suggest an important participation of polyP in the oxidative stress response of the cells and as having a possible metabolic role in the cell, as previously described in bacteria. The polyP deficient strain also showed an increased copper sensitivity and an earlier transcriptional up-regulation of copA, implying a complementary role of both copper resistance systems.

Autophosphorylation and Cross-Phosphorylation of Protein Kinases from the Crenarchaeon Sulfolobus islandicus

Protein phosphorylation, one of the most important post-translational modifications, regulates almost every cellular process. Although signal transduction by protein phosphorylation is extensively studied in Eukaryotes and Bacteria, the knowledge of this process in archaea is greatly lagging behind, especially for Ser/Thr/Tyr phosphorylation by eukaryotic-like protein kinases (ePKs). So far, only a few studies on archaeal ePKs have been reported, most of which focused on the phosphorylation activities in vitro, but their physiological functions and interacting network are still largely unknown. In this study, we systematically investigated the autophosphorylation and cross-phosphorylation activities of ePKs from Sulfolobus islandicus REY15A using proteins expressed in Escherichia coli or S. islandicus. In vitro kinase assay showed that 7 out of the 11 putative ePKs have autophosphorylation activity. A protein Ser/Thr phosphatase, SiRe_1009, was able to dephosphorylate various autophosphorylated ePKs, confirming that these proteins are Ser/Thr kinases. Two ePKs, SiRe_2030 and SiRe_2056, homologs of typical eukaryotic PKs involved in peptide synthesis in response to various cellular stresses, exhibit highly efficient phosphorylation activities on both themselves and other ePKs. Overexpression of the protein kinases in vivo revealed that elevated level of either SiRe_1531 or SiRe_2056 inhibited the cell growth of S. islandicus cells. Finally, a phosphorylation network of the protein kinases was proposed and their putative physiological roles were discussed.

Natural Mutagenesis-Enabled Global Proteomic Study of Metabolic and Carbon Source Implications in Mutant Thermoacidophillic Archaeon Sulfolobus solfataricus PBL2025

The thermoacidophilic crenarchaeon Sulfolobus solfataricus has been widely used as a model organism for archaeal systems biology research. Investigation using its spontaneous mutant PBL2025 provides an effective metabolic baseline to study subsequent mutagenesis-induced functional process shifts as well as changes in feedback inhibitions. Here, an untargeted metabolic investigation using quantitative proteomics and metabolomics was performed to correlate changes in S. solfataricus strains P2 against PBL2025 and under both glucose and tryptone. The study is combined with pathway enrichment analysis to identify prominent proteins with differential stoichiometry. Proteome level quantification reveals that over 20% of the observed overlapping proteome is differentially expressed under these conditions. Metabolic-induced differential expressions are observed along the central carbon metabolism, along with 12 other significantly regulated pathways. Current findings suggest that PBL2025 is able to compensate through the induction of carbon metabolism, as well as other anabolic pathways such as Val, Leu and iso-Leu biosynthesis. Studying protein abundance changes after changes in carbon sources also reveals distinct differences in metabolic strategies employed by both strains, whereby a clear down-regulation of carbohydrate and nucleotide metabolism is observed for P2, while a mixed response through down-regulation of energy formation and up-regulation of glycolysis is observed for PBL2025. This study contributes, to date, the most comprehensive network of changes in carbohydrate and amino acid pathways using the complementary systems biology observations at the protein and metabolite levels. Current findings provide a unique insight into molecular processing changes through natural (spontaneous) metabolic rewiring, as well as a systems biology understanding of the metabolic elasticity of thermoacidophiles to environmental carbon source change, potentially guiding more efficient directed mutagenesis in archaea.

Wing phosphorylation is a major functional determinant of the Lrs14-type biofilm and motility regulator AbfR1 in Sulfolobus acidocaldarius

In response to a variety of environmental cues, prokaryotes can switch between a motile and a sessile, biofilm-forming mode of growth. The regulatory mechanisms and signaling pathways underlying this switch are largely unknown in archaea but involve small winged helix-turn-helix DNA-binding proteins of the archaea-specific Lrs14 family. Here, we study the Lrs14 member AbfR1 of Sulfolobus acidocaldarius. Small-angle X-ray scattering data are presented, which are consistent with a model of dimeric AbfR1 in which dimerization occurs via an antiparallel coiled coil as suggested by homology modeling. Furthermore, solution structure data of AbfR1-DNA complexes suggest that upon binding DNA, AbfR1 induces deformations in the DNA. The wing residues tyrosine 84 and serine 87, which are phosphorylated in vivo, are crucial to establish stable protein-DNA contacts and their substitution with a negatively charged glutamate or aspartate residue inhibits formation of a nucleoprotein complex. Furthermore, mutation abrogates the cellular abundance and transcription regulatory function of AbfR1 and thus affects the resulting biofilm and motility phenotype of S. acidocaldarius. This work establishes a novel wHTH DNA-binding mode for Lrs14-like proteins and hints at an important role for protein phosphorylation as a signal transduction mechanism for the control of biofilm formation and motility in archaea.

ArnS, a kinase involved in starvation-induced archaellum expression

Organisms have evolved motility organelles that allow them to move to favourable habitats. Cells integrate environmental stimuli into intracellular signals to motility machineries to direct this migration. Many motility organelles are complex surface appendages that have evolved a tight, hierarchical regulation of expression. In the crenearchaeon Sulfolobus acidocaldarius, biosynthesis of the archaellum is regulated by regulatory network proteins that control expression of archaellum components in a phosphorylation-dependent manner. A major trigger for archaellum expression is nutrient starvation, but although some components are known, the regulatory cascade triggered by starvation is poorly understood. In this work, the starvation-induced Ser/Thr protein kinase ArnS (Saci_1181) which is located proximally to the archaellum operon was identified. Deletion of arnS results in reduced motility, though the archaellum is properly assembled. Therefore, our experimental and modelling results indicate that ArnS plays an essential role in the precisely controlled expression of archaellum components during starvation-induced motility in Sulfolobus acidocaldarius. Furthermore they combined in vivo experiments and mathematical models to describe for the first time in archaea the dynamics of key regulators of archaellum expression.

Abundant Lysine Methylation and N-Terminal Acetylation in Sulfolobus islandicus Revealed by Bottom-Up and Top-Down Proteomics

Protein post-translational methylation has been reported to occur in archaea, including members of the genus Sulfolobus, but has never been characterized on a proteome-wide scale. Among important Sulfolobus proteins carrying such modification are the chromatin proteins that have been described to be methylated on lysine side chains, resembling eukaryotic histones in that aspect. To get more insight into the extent of this modification and its dynamics during the different growth steps of the thermoacidophylic archaeon S. islandicus LAL14/1, we performed a global and deep proteomic analysis using a combination of high-throughput bottom-up and top-down approaches on a single high-resolution mass spectrometer. 1,931 methylation sites on 751 proteins were found by the bottom-up analysis, with methylation sites on 526 proteins monitored throughout three cell culture growth stages: early-exponential, mid-exponential, and stationary. The top-down analysis revealed 3,978 proteoforms arising from 681 proteins, including 292 methylated proteoforms, 85 of which were comprehensively characterized. Methylated proteoforms of the five chromatin proteins (Alba1, Alba2, Cren7, Sul7d1, Sul7d2) were fully characterized by a combination of bottom-up and top-down data. The top-down analysis also revealed an increase of methylation during cell growth for two chromatin proteins, which had not been evidenced by bottom-up. These results shed new light on the ubiquitous lysine methylation throughout the S. islandicus proteome. Furthermore, we found that S. islandicus proteins are frequently acetylated at the N terminus, following the removal of the N-terminal methionine. This study highlights the great value of combining bottom-up and top-down proteomics for obtaining an unprecedented level of accuracy in detecting differentially modified intact proteoforms. The data have been deposited to the ProteomeXchange with identifiers PXD003074 and PXD004179.

Expanding the archaellum regulatory network - the eukaryotic protein kinases ArnC and ArnD influence motility of Sulfolobus acidocaldarius

Expression of the archaellum, the archaeal‐type IV pilus‐like rotating motility structure is upregulated under nutrient limitation. This is controlled by a network of regulators, called the archaellum regulatory network (arn). Several of the components of this network in Sulfolobus acidocaldarius can be phosphorylated, and the deletion of the phosphatase PP2A results in strongly increased motility during starvation, indicating a role for phosphorylation in the regulation of motility. Analysis of the motility of different protein kinase deletion strains revealed that deletion of saci_0965, saci_1181, and saci_1193 resulted in reduced motility, whereas the deletion of saci_1694 resulted in hypermotility. Here ArnC (Saci_1193) and ArnD (Saci_1694) are characterized. Purified ArnC and ArnD phosphorylate serine and threonine residues in the C‐terminus of the repressor ArnB. arnC is upregulated in starvation medium, whereas arnD is constitutively expressed. However, while differences in the expression and levels of flaB were observed in the ΔarnD strain during growth under rich conditions, under nutrient limiting conditions the ΔarnC and ΔarnD strains showed no large differences in the expression levels of the archaellum or of the studied regulators. This suggests that next to the regulation via the archaellum regulatory network additional regulatory mechanisms of expression and/or activity of the archaellum exist.

Protein phosphorylation and its role in archaeal signal transduction

Reversible protein phosphorylation is the main mechanism of signal transduction that enables cells to rapidly respond to environmental changes by controlling the functional properties of proteins in response to external stimuli. However, whereas signal transduction is well studied in Eukaryotes and Bacteria, the knowledge in Archaea is still rather scarce. Archaea are special with regard to protein phosphorylation, due to the fact that the two best studied phyla, the Euryarchaeota and Crenarchaeaota, seem to exhibit fundamental differences in regulatory systems. Euryarchaeota (e.g. halophiles, methanogens, thermophiles), like Bacteria and Eukaryotes, rely on bacterial-type two-component signal transduction systems (phosphorylation on His and Asp), as well as on the protein phosphorylation on Ser, Thr and Tyr by Hanks-type protein kinases. Instead, Crenarchaeota (e.g. acidophiles and (hyper)thermophiles) only depend on Hanks-type protein phosphorylation. In this review, the current knowledge of reversible protein phosphorylation in Archaea is presented. It combines results from identified phosphoproteins, biochemical characterization of protein kinases and protein phosphatases as well as target enzymes and first insights into archaeal signal transduction by biochemical, genetic and polyomic studies.

The authors review the current knowledge about protein phosphorylation in Archaea and its impact on signaling in this organism group.

Role of Protein Phosphorylation in the Regulation of Cell Cycle and DNA-Related Processes in Bacteria

In all living organisms, the phosphorylation of proteins modulates various aspects of their functionalities. In eukaryotes, protein phosphorylation plays a key role in cell signaling, gene expression, and differentiation. Protein phosphorylation is also involved in the global control of DNA replication during the cell cycle, as well as in the mechanisms that cope with stress-induced replication blocks. Similar to eukaryotes, bacteria use Hanks-type kinases and phosphatases for signal transduction, and protein phosphorylation is involved in numerous cellular processes. However, it remains unclear whether protein phosphorylation in bacteria can also regulate the activity of proteins involved in DNA-mediated processes such as DNA replication or repair. Accumulating evidence supported by functional and biochemical studies suggests that phospho-regulatory mechanisms also take place during the bacterial cell cycle. Recent phosphoproteomics and interactomics studies identified numerous phosphoproteins involved in various aspect of DNA metabolism strongly supporting the existence of such level of regulation in bacteria. Similar to eukaryotes, bacterial scaffolding-like proteins emerged as platforms for kinase activation and signaling. This review reports the current knowledge on the phosphorylation of proteins involved in the maintenance of genome integrity and the regulation of cell cycle in bacteria that reveals surprising similarities to eukaryotes.

Dynamic Metabolite Profiling in an Archaeon Connects Transcriptional Regulation to Metabolic Consequences

Previous work demonstrated that the TrmB transcription factor is responsible for regulating the expression of many enzyme-coding genes in the hypersaline-adapted archaeon Halobacterium salinarum via a direct interaction with a cis-regulatory sequence in their promoters. This interaction is abolished in the presence of glucose. Although much is known about the effects of TrmB at the transcriptional level, it remains unclear whether and to what extent changes in mRNA levels directly affect metabolite levels. In order to address this question, here we performed a high-resolution metabolite profiling time course during a change in nutrients using a combination of targeted and untargeted methods in wild-type and ΔtrmB strain backgrounds. We found that TrmB-mediated transcriptional changes resulted in widespread and significant changes to metabolite levels across the metabolic network. Additionally, the pattern of growth complementation using various purines suggests that the mis-regulation of gluconeogenesis in the ΔtrmB mutant strain in the absence of glucose results in low phosphoribosylpyrophosphate (PRPP) levels. We confirmed these low PRPP levels using a quantitative mass spectrometric technique and found that they are associated with a metabolic block in de novo purine synthesis, which is partially responsible for the growth defect of the ΔtrmB mutant strain in the absence of glucose. In conclusion, we show how transcriptional regulation of metabolism affects metabolite levels and ultimately, phenotypes.

A specific proteomic response of Sulfolobus solfataricus P2 to gamma radiations

Sulfolobus solfataricus is an acidophilic hyperthermophilic crenarchaeon living at 80 °C in aerobic conditions. As other thermophilic organisms, S. solfataricus is resistant to gamma irradiation and we studied the response of this microorganism to this ionizing irradiation by monitoring cell growth, DNA integrity and proteome variations. In aerobic conditions, the S. solfataricus genome was fragmented due to the multiple DNA double strand breakages induced by γ-rays and was fully restored within a couple of hours. Comparison of irradiated and unirradiated cell proteomes indicated that only few proteins changed. The proteins identified by mass spectrometry are involved in different cellular pathways including DNA replication, recombination and repair. Interestingly, we observed that some proteins are irradiation dose-specific while others are common to the cell response regardless of the irradiation dose. Most of the proteins highlighted in these conditions seem to act together to allow an efficient cell response to γ-irradiation.

Activation of SsoPK4, an Archaeal eIF2α Kinase Homolog, by Oxidized CoA

The eukaryotic protein kinase (ePK) paradigm provides integral components for signal transduction cascades throughout nature. However, while so-called typical ePKs permeate the Eucarya and Bacteria, atypical ePKs dominate the kinomes of the Archaea. Intriguingly, the catalytic domains of the handful of deduced typical ePKs from the archaeon Sulfolobus solfataricus P2 exhibit significant resemblance to the protein kinases that phosphorylate translation initiation factor 2α (eIF2α) in response to cellular stresses. We cloned and expressed one of these archaeal eIF2α protein kinases, SsoPK4. SsoPK4 exhibited protein-serine/threonine kinase activity toward several proteins, including the S. solfataricus homolog of eIF2α, aIF2α. The activity of SsoPK4 was inhibited in vitro by 3ʹ,5ʹ-cyclic AMP (K_i of ~23 µM) and was activated by oxidized Coenzyme A, an indicator of oxidative stress in the Archaea. Activation enhanced the apparent affinity for protein substrates, K_m, but had little effect on V_max. Autophosphorylation activated SsoPK4 and rendered it insensitive to oxidized Coenzyme A.

Archaeal Clusters of Orthologous Genes (arCOGs): An Update and Application for Analysis of Shared Features between Thermococcales, Methanococcales, and Methanobacteriales

With the continuously accelerating genome sequencing from diverse groups of archaea and bacteria, accurate identification of gene orthology and availability of readily expandable clusters of orthologous genes are essential for the functional annotation of new genomes. We report an update of the collection of archaeal Clusters of Orthologous Genes (arCOGs) to cover, on average, 91% of the protein-coding genes in 168 archaeal genomes. The new arCOGs were constructed using refined algorithms for orthology identification combined with extensive manual curation, including incorporation of the results of several completed and ongoing research projects in archaeal genomics. A new level of classification is introduced, superclusters that unit two or more arCOGs and more completely reflect gene family evolution than individual, disconnected arCOGs. Assessment of the current archaeal genome annotation in public databases indicates that consistent use of arCOGs can significantly improve the annotation quality. In addition to their utility for genome annotation, arCOGs also are a platform for phylogenomic analysis. We explore this aspect of arCOGs by performing a phylogenomic study of the Thermococci that are traditionally viewed as the basal branch of the Euryarchaeota. The results of phylogenomic analysis that involved both comparison of multiple phylogenetic trees and a search for putative derived shared characters by using phyletic patterns extracted from the arCOGs reveal a likely evolutionary relationship between the Thermococci, Methanococci, and Methanobacteria. The arCOGs are expected to be instrumental for a comprehensive phylogenomic study of the archaea.

Insight into the cellular involvement of the two reverse gyrases from the hyperthermophilic archaeon Sulfolobus solfataricus

Background

Reverse gyrases are DNA topoisomerases characterized by their unique DNA positive-supercoiling activity. Sulfolobus solfataricus, like most Crenarchaeota, contains two genes each encoding a reverse gyrase. We showed previously that the two genes are differently regulated according to temperature and that the corresponding purified recombinant reverse gyrases have different enzymatic characteristics. These observations suggest a specialization of functions of the two reverse gyrases. As no mutants of the TopR genes could be obtained in Sulfolobales, we used immunodetection techniques to study the function(s) of these proteins in S. solfataricus in vivo. In particular, we investigated whether one or both reverse gyrases are required for the hyperthermophilic lifestyle.

Results

For the first time the two reverse gyrases of S. solfataricus have been discriminated at the protein level and their respective amounts have been determined in vivo. Actively dividing S. solfataricus cells contain only small amounts of both reverse gyrases, approximately 50 TopR1 and 125 TopR2 molecules per cell at 80°C. S. solfataricus cells are resistant at 45°C for several weeks, but there is neither cell division nor replication initiation; these processes are fully restored upon a return to 80°C. TopR1 is not found after three weeks at 45°C whereas the amount of TopR2 remains constant. Enzymatic assays in vitro indicate that TopR1 is not active at 45°C but that TopR2 exhibits highly positive DNA supercoiling activity at 45°C.

Conclusions

The two reverse gyrases of S. solfataricus are differently regulated, in terms of protein abundance, in vivo at 80°C and 45°C. TopR2 is present both at high and low temperatures and is therefore presumably required whether cells are dividing or not. By contrast, TopR1 is present only at high temperature where the cell division occurs, suggesting that TopR1 is required for controlling DNA topology associated with cell division activity and/or life at high temperature. Our findings in vitro that TopR1 is able to positively supercoil DNA only at high temperature, and TopR2 is active at both temperatures are consistent with them having different functions within the cells.

Agents of change - concepts in Mycobacterium tuberculosis Ser/Thr/Tyr phosphosignalling

The flow of information from the outside to the inside of bacterial cells is largely directed by protein kinases. In addition to histidine/aspartate phosphorelays of two-component response regulators, recent work in Mycobacterium tuberculosis (Mtb) reinforces the idea that phosphorylation on serine (Ser), threonine (Thr) and tyrosine (Tyr) is central to bacterial physiology and pathogenesis, and that the corresponding phosphosystems are highly similar to those in eukaryotes. In this way, eukaryotes are a useful guide to understanding Ser/Thr/Tyr phosphorylation (O-phosphorylation) in prokaryotes such as Mtb. However, as novel functions and components of bacterial O-phosphorylation are identified, distinct differences between pro- and eukaryotic phosphosignalling systems become apparent. The emerging picture of O-phosphorylation in Mtb is complicated, goes beyond the eukaryotic paradigms, and shows the limitations of viewing bacterial phosphosignalling within the confines of the 'eukaryotic-like' model. Here, we summarize recent findings about Ser/Thr and the recently discovered Tyr phosphorylation pathways in Mtb, highlight the similarities and differences between eukaryotic and prokaryotic O-phosphorylation, and pose additional questions about signalling components, pathway organization, and ultimately, the cellular roles of O-phosphorylation in Mtb physiology and pathogenesis.

Carbohydrate metabolism in Archaea: current insights into unusual enzymes and pathways and their regulation

The metabolism of Archaea, the third domain of life, resembles in its complexity those of Bacteria and lower Eukarya. However, this metabolic complexity in Archaea is accompanied by the absence of many "classical" pathways, particularly in central carbohydrate metabolism. Instead, Archaea are characterized by the presence of unique, modified variants of classical pathways such as the Embden-Meyerhof-Parnas (EMP) pathway and the Entner-Doudoroff (ED) pathway. The pentose phosphate pathway is only partly present (if at all), and pentose degradation also significantly differs from that known for bacterial model organisms. These modifications are accompanied by the invention of "new," unusual enzymes which cause fundamental consequences for the underlying regulatory principles, and classical allosteric regulation sites well established in Bacteria and Eukarya are lost. The aim of this review is to present the current understanding of central carbohydrate metabolic pathways and their regulation in Archaea. In order to give an overview of their complexity, pathway modifications are discussed with respect to unusual archaeal biocatalysts, their structural and mechanistic characteristics, and their regulatory properties in comparison to their classic counterparts from Bacteria and Eukarya. Furthermore, an overview focusing on hexose metabolic, i.e., glycolytic as well as gluconeogenic, pathways identified in archaeal model organisms is given. Their energy gain is discussed, and new insights into different levels of regulation that have been observed so far, including the transcript and protein levels (e.g., gene regulation, known transcription regulators, and posttranslational modification via reversible protein phosphorylation), are presented.

Protein Ser/Thr/Tyr phosphorylation in the Archaea

The third domain of life, the Archaea (formerly Archaebacteria), is populated by a physiologically diverse set of microorganisms, many of which reside at the ecological extremes of our global environment. Although ostensibly prokaryotic in morphology, the Archaea share much closer evolutionary ties with the Eukarya than with the superficially more similar Bacteria. Initial genomic, proteomic, and biochemical analyses have revealed the presence of "eukaryotic" protein kinases and phosphatases and an intriguing set of serine-, threonine-, and tyrosine-phosphorylated proteins in the Archaea that may offer new insights into this important regulatory mechanism.

Structure of the CRISPR interference complex CSM reveals key similarities with cascade

The Clustered Regularly Interspaced Palindromic Repeats (CRISPR) system is an adaptive immune system in prokaryotes. Interference complexes encoded by CRISPR-associated (cas) genes utilize small RNAs for homology-directed detection and subsequent degradation of invading genetic elements, and they have been classified into three main types (I-III). Type III complexes share the Cas10 subunit but are subclassifed as type IIIA (CSM) and type IIIB (CMR), depending on their specificity for DNA or RNA targets, respectively. The role of CSM in limiting the spread of conjugative plasmids in Staphylococcus epidermidis was first described in 2008. Here, we report a detailed investigation of the composition and structure of the CSM complex from the archaeon Sulfolobus solfataricus, using a combination of electron microscopy, mass spectrometry, and deep sequencing. This reveals a three-dimensional model for the CSM complex that includes a helical component strikingly reminiscent of the backbone structure of the type I (Cascade) family.

Archaeal signal transduction: impact of protein phosphatase deletions on cell size, motility, and energy metabolism in Sulfolobus acidocaldarius

In this study, the in vitro and in vivo functions of the only two identified protein phosphatases, Saci-PTP and Saci-PP2A, in the crenarchaeal model organism Sulfolobus acidocaldarius were investigated. Biochemical characterization revealed that Saci-PTP is a dual-specific phosphatase (against pSer/pThr and pTyr), whereas Saci-PP2A exhibited specific pSer/pThr activity and inhibition by okadaic acid. Deletion of saci_pp2a resulted in pronounced alterations in growth, cell shape and cell size, which could be partially complemented. Transcriptome analysis of the three strains (Δsaci_ptp, Δsaci_pp2a and the MW001 parental strain) revealed 155 genes that were differentially expressed in the deletion mutants, and showed significant changes in expression of genes encoding the archaella (archaeal motility structure), components of the respiratory chain and transcriptional regulators. Phosphoproteome studies revealed 801 unique phosphoproteins in total, with an increase in identified phosphopeptides in the deletion mutants. Proteins from most functional categories were affected by phosphorylation, including components of the motility system, the respiratory chain, and regulatory proteins. In the saci_pp2a deletion mutant the up-regulation at the transcript level, as well as the observed phosphorylation pattern, resembled starvation stress responses. Hypermotility was also observed in the saci_pp2a deletion mutant. The results highlight the importance of protein phosphorylation in regulating essential cellular processes in the crenarchaeon S. acidocaldarius.

Atypical protein kinases of the RIO family in archaea

Atypical protein kinases of the RIO (right open reading frame) kinase family are found in all three domains of life, emphasizing their essential function. In all archaeal genomes sequenced to date, typically two, but at least one, members of the RIO kinase family have been identified. Although the function of RIO kinases in Archaea remains to be resolved, bioinformatics analysis (e.g. comparison of the phylogenetic distribution and gene neighbourhood analysis, as well as interaction analysis) in combination with the available phosphoproteome study of Sulfolobus solfataricus provided some first hints to the possible function as well as revealed some putative target proteins for RIO kinases in Archaea. This study suggests a possible function of archaeal RIO kinases in RNA and/or DNA binding/processing translation initiation or ribosomal biogenesis resembling the assumed physiological role in yeast.

Protein-DNA binding dynamics predict transcriptional response to nutrients in archaea

Organisms across all three domains of life use gene regulatory networks (GRNs) to integrate varied stimuli into coherent transcriptional responses to environmental pressures. However, inferring GRN topology and regulatory causality remains a central challenge in systems biology. Previous work characterized TrmB as a global metabolic transcription factor in archaeal extremophiles. However, it remains unclear how TrmB dynamically regulates its ∼100 metabolic enzyme-coding gene targets. Using a dynamic perturbation approach, we elucidate the topology of the TrmB metabolic GRN in the model archaeon Halobacterium salinarum. Clustering of dynamic gene expression patterns reveals that TrmB functions alone to regulate central metabolic enzyme-coding genes but cooperates with various regulators to control peripheral metabolic pathways. Using a dynamical model, we predict gene expression patterns for some TrmB-dependent promoters and infer secondary regulators for others. Our data suggest feed-forward gene regulatory topology for cobalamin biosynthesis. In contrast, purine biosynthesis appears to require TrmB-independent regulators. We conclude that TrmB is an important component for mediating metabolic modularity, integrating nutrient status and regulating gene expression dynamics alone and in concert with secondary regulators.

Archaeal aldehyde dehydrogenase ST0064 from Sulfolobus tokodaii, a paralog of non-phosphorylating glyceraldehyde-3-phosphate dehydrogenase, is a succinate semialdehyde dehydrogenase

Aldehyde dehydrogenase ST0064, the closest paralog of previously characterized allosteric non-phosphorylating glyceraldehyde-3-phosphate (GAP) dehydrogenase (GAPN, ST2477) from a thermoacidophilic archaeon, Sulfolobus tokodaii, was expressed heterologously and characterized in detail. ST0064 showed remarkable activity toward succinate semialdehyde (SSA) (Km of 0.0029 mM and kcat of 30.0 s(-1)) with no allosteric regulation. Activity toward GAP was lower (Km of 4.6 mM and kcat of 4.77 s(-1)), and previously predicted succinyl-CoA reductase activity was not detected, suggesting that the enzyme functions practically as succinate semialdehyde dehydrogenase (SSADH). Phylogenetic analysis indicated that archaeal SSADHs and GAPNs are closely related within the aldehyde dehydrogenase superfamily, suggesting that they are of the same origin.

Genome-scale reconstruction and analysis of the metabolic network in the hyperthermophilic archaeon Sulfolobus solfataricus

We describe the reconstruction of a genome-scale metabolic model of the crenarchaeon Sulfolobus solfataricus, a hyperthermoacidophilic microorganism. It grows in terrestrial volcanic hot springs with growth occurring at pH 2–4 (optimum 3.5) and a temperature of 75–80°C (optimum 80°C). The genome of Sulfolobus solfataricus P2 contains 2,992,245 bp on a single circular chromosome and encodes 2,977 proteins and a number of RNAs. The network comprises 718 metabolic and 58 transport/exchange reactions and 705 unique metabolites, based on the annotated genome and available biochemical data. Using the model in conjunction with constraint-based methods, we simulated the metabolic fluxes induced by different environmental and genetic conditions. The predictions were compared to experimental measurements and phenotypes of S. solfataricus. Furthermore, the performance of the network for 35 different carbon sources known for S. solfataricus from the literature was simulated. Comparing the growth on different carbon sources revealed that glycerol is the carbon source with the highest biomass flux per imported carbon atom (75% higher than glucose). Experimental data was also used to fit the model to phenotypic observations. In addition to the commonly known heterotrophic growth of S. solfataricus, the crenarchaeon is also able to grow autotrophically using the hydroxypropionate-hydroxybutyrate cycle for bicarbonate fixation. We integrated this pathway into our model and compared bicarbonate fixation with growth on glucose as sole carbon source. Finally, we tested the robustness of the metabolism with respect to gene deletions using the method of Minimization of Metabolic Adjustment (MOMA), which predicted that 18% of all possible single gene deletions would be lethal for the organism.