Succinylome Analysis Reveals the Involvement of Lysine Succinylation in Metabolism in Pathogenic Mycobacterium tuberculosis

Mycobacterium tuberculosis (Mtb), the causative agent of human tuberculosis, remains one of the most prevalent human pathogens and a major cause of mortality worldwide. Metabolic network is a central mediator and defining feature of the pathogenicity of Mtb. Increasing evidence suggests that lysine succinylation dynamically regulates enzymes in carbon metabolism in both bacteria and human cells; however, its extent and function in Mtb remain unexplored. Here, we performed a global succinylome analysis of the virulent Mtb strain H37Rv by using high accuracy nano-LC-MS/MS in combination with the enrichment of succinylated peptides from digested cell lysates and subsequent peptide identification. In total, 1545 lysine succinylation sites on 626 proteins were identified in this pathogen. The identified succinylated proteins are involved in various biological processes and a large proportion of the succinylation sites are present on proteins in the central metabolism pathway. Site-specific mutations showed that succinylation is a negative regulatory modification on the enzymatic activity of acetyl-CoA synthetase. Molecular dynamics simulations demonstrated that succinylation affects the conformational stability of acetyl-CoA synthetase, which is critical for its enzymatic activity. Further functional studies showed that CobB, a sirtuin-like deacetylase in Mtb, functions as a desuccinylase of acetyl-CoA synthetase in in vitro assays. Together, our findings reveal widespread roles for lysine succinylation in regulating metabolism and diverse processes in Mtb. Our data provide a rich resource for functional analyses of lysine succinylation and facilitate the dissection of metabolic networks in this life-threatening pathogen.Post-translational modifications (PTMs)¹ are complex and fundamental mechanisms modulating diverse protein properties and functions, and have been associated with almost all known cellular pathways and disease processes (1, 2). Among the hundreds of different PTMs, acylations at lysine residues, such as acetylation (3–6), malonylation (7, 8), crotonylation (9, 10), propionylation (11–13), butyrylation (11, 13), and succinylation (7, 14–16) are crucial for functional regulations of many prokaryotic and eukaryotic proteins. Because these lysine PTMs depend on the acyl-CoA metabolic intermediates, such as acetyl-CoA (Ac-CoA), succinyl-CoA, and malonyl-CoA, lysine acylation could provide a mechanism to respond to changes in the energy status of the cell and regulate energy metabolism and the key metabolic pathways in diverse organisms (17, 18).Among these lysine PTMs, lysine succinylation is a highly dynamic and regulated PTM defined as transfer of a succinyl group (-CO-CH₂-CH₂-CO-) to a lysine residue of a protein molecule (8). It was recently identified and comprehensively validated in both bacterial and mammalian cells (8, 14, 16). It was also identified in core histones, suggesting that lysine succinylation may regulate the functions of histones and affect chromatin structure and gene expression (7). Accumulating evidence suggests that lysine succinylation is a widespread and important PTM in both eukaryotes and prokaryotes and regulates diverse cellular processes (16). The system-wide studies involving lysine-succinylated peptide immunoprecipitation and liquid chromatography-mass spectrometry (LC-MS/MS) have been employed to analyze the bacteria (E. coli) (14, 16), yeast (S. cerevisiae), human (HeLa) cells, and mouse embryonic fibroblasts and liver cells (16, 19). These succinylome studies have generated large data sets of lysine-succinylated proteins in both eukaryotes and prokaryotes and demonstrated the diverse cellular functions of this PTM. Notably, lysine succinylation is widespread among diverse mitochondrial metabolic enzymes that are involved in fatty acid metabolism, amino acid degradation, and the tricarboxylic acid cycle (19, 20). Thus, lysine succinylation is reported as a functional PTM with the potential to impact mitochondrial metabolism and coordinate different metabolic pathways in human cells and bacteria (14, 19–22).Mycobacterium tuberculosis (Mtb), the causative agent of tuberculosis (TB), is a major cause of mortality worldwide and claims more human lives annually than any other bacterial pathogen (23). About one third of the world''s population is infected with Mtb, which leads to nearly 1.3 million deaths and 8.6 million new cases of TB in 2012 worldwide (24). Mtb remains a major threat to global health, especially in the developing countries. Emergence of multidrug resistant (MDR) and extensively drug-resistant (XDR) Mtb, and also the emergence of co-infection between TB and HIV have further worsened the situation (25–27). Among bacterial pathogens, Mtb has a distinctive life cycle spanning different environments and developmental stages (28). Especially, Mtb can exist in dormant or active states in the host, leading to asymptomatic latent TB infection or active TB disease (29). To achieve these different physiologic states, Mtb developed a mechanism to sense diverse signals from the host and to coordinately regulate multiple cellular processes and pathways (30, 31). Mtb has evolved its metabolic network to both maintain and propagate its survival as a species within humans (32–35). It is well accepted that metabolic network is a central mediator and defining feature of the pathogenicity of Mtb (23, 36–38). Knowledge of the regulation of metabolic pathways used by Mtb during infection is therefore important for understanding its pathogenicity, and can also guide the development of novel drug therapies (39). On the other hand, increasing evidence suggests that lysine succinylation dynamically regulates enzymes in carbon metabolism in both bacteria and human cells (14, 19–22). It is tempting to speculate that lysine succinylation may play an important regulatory role in metabolic processes in Mtb. However, to the best of our knowledge, no succinylated protein in Mtb has been identified, presenting a major obstacle to understand the regulatory roles of lysine succinylation in this life-threatening pathogen.In order to fill this gap in our knowledge, we have initiated a systematic study of the identities and functional roles of the succinylated protein in Mtb. Because Mtb H37Rv is the first sequenced Mtb strain (40) and has been extensively used for studies in dissecting the roles of individual genes in pathogenesis (41), it was selected as a test case. We analyzed the succinylome of Mtb H37Rv by using high accuracy nano-LC-MS/MS in combination with the enrichment of succinylated peptides from digested cell lysates and subsequent peptide identification. In total, 1545 lysine succinylation sites on 626 proteins were identified in this pathogen. The identified succinylated proteins are involved in various biological processes and render particular enrichment to metabolic process. A large proportion of the succinylation sites are present on proteins in the central metabolism pathway. We further dissected the regulatory role of succinylation on acetyl-CoA synthetase (Acs) via site-specific mutagenesis analysis and molecular dynamics (MD) simulations showed that reversible lysine succinylation could inhibit the activity of Acs. Further functional studies showed that CobB, a sirtuin-like deacetylase in Mtb, functions as a deacetylase and as a desuccinylase of Acs in in vitro assays. Together, our findings provide significant insights into the range of functions regulated by lysine succinylation in Mtb.     

Splitting the BLOSUM Score into Numbers of Biological Significance

Mathematical tools developed in the context of Shannon information theory were used to analyze the meaning of the BLOSUM score, which was split into three components termed as the BLOSUM spectrum (or BLOSpectrum). These relate respectively to the sequence convergence (the stochastic similarity of the two protein sequences), to the background frequency divergence (typicality of the amino acid probability distribution in each sequence), and to the target frequency divergence (compliance of the amino acid variations between the two sequences to the protein model implicit in the BLOCKS database). This treatment sharpens the protein sequence comparison, providing a rationale for the biological significance of the obtained score, and helps to identify weakly related sequences. Moreover, the BLOSpectrum can guide the choice of the most appropriate scoring matrix, tailoring it to the evolutionary divergence associated with the two sequences, or indicate if a compositionally adjusted matrix could perform better.[1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29]     

Analysis of Free Energy Signals Arising from Nucleotide Hybridization Between rRNA and mRNA Sequences during Translation in Eubacteria

A decoding algorithm is tested that mechanistically models the progressive alignments that arise as the mRNA moves past the rRNA tail during translation elongation. Each of these alignments provides an opportunity for hybridization between the single-stranded, -terminal nucleotides of the 16S rRNA and the spatially accessible window of mRNA sequence, from which a free energy value can be calculated. Using this algorithm we show that a periodic, energetic pattern of frequency 1/3 is revealed. This periodic signal exists in the majority of coding regions of eubacterial genes, but not in the non-coding regions encoding the 16S and 23S rRNAs. Signal analysis reveals that the population of coding regions of each bacterial species has a mean phase that is correlated in a statistically significant way with species () content. These results suggest that the periodic signal could function as a synchronization signal for the maintenance of reading frame and that codon usage provides a mechanism for manipulation of signal phase.[1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32]     

A Combined Proteomics and Metabolomics Profiling of Gastric Cardia Cancer Reveals Characteristic Dysregulations in Glucose Metabolism

Gastric cardia cancer (GCC), which occurs at the gastric-esophageal boundary, is one of the most malignant tumors. Despite its high mortality and morbidity, the molecular mechanism of initiation and progression of this disease is largely unknown. In this study, using proteomics and metabolomics approaches, we found that the level of several enzymes and their related metabolic intermediates involved in glucose metabolism were deregulated in GCC. Among these enzymes, two subunits controlling pyruvic acid efflux, lactate dehydrogenase A (LDHA) and pyruvate dehydrogenase B (PDHB), were further analyzed in vitro. Either down-regulation of LDH subunit LDHA or overexpression of PDH subunit PDHB could force pyruvic acid into the Krebs cycle rather than the glycolysis process in AGS gastric cancer cells, which inhibited cell growth and cell migration. Our results reflect an important glucose metabolic signature, especially the dysregulation of pyruvic acid efflux in the development of GCC. Forced transition from glycolysis to the Krebs cycle had an inhibitory effect on GCC progression, providing potential therapeutic targets for this disease.Gastric cardia cancer (GCC),1 which occurs at the gastric-esophageal boundary, is one of the most malignant tumors. Despite the steadily falling incidence of gastric non-cardia cancer in the past two decades (1), the rate of GCC has risen rapidly, establishing gastric cancer as the second major cause of cancer-related deaths throughout the world (2). GCC has become a significant cause of mortality and morbidity both in the west (3–5) and in Asia (6, 7), especially in China (8). Although this cancer has become an important health problem worldwide, the its pathogenesis has not been well characterized (1). Most patients are diagnosed at an advanced stage, contributing to the high mortality rate of the disease.Systematic proteomics analysis has proved to be a powerful approach in a variety of human cancer research, including lung (9), esophagus (10), gastric (11), liver (12), breast (13), and brain cancer (14). Metabolomics, another new bio-omics technology recently introduced into cancer research (15), is the global analysis of the small metabolites produced by normal or pathologic cellular processes. Some metabolic intermediates have been identified as new cancer biomarkers (16).Using proteomics and metabolomics methods in this study, we found that a series of proteins and metabolic intermediates, mainly involved in glucose metabolism, were altered during the development of GCC. The high activity of anaerobic glycolysis and the impairment of aerobic respiration occurring in these cells recapitulated the Warburg effect (17). Further studies using a gastric cancer cell line demonstrated that the predominant anaerobic glycolysis was essential for tumor cells to sustain rapid proliferation, whereas forced transition from anaerobic glycolysis to aerobic respiration inhibited the growth of tumor cells. In conclusion, our study revealed the major metabolic alterations essential for the development of GCC and discovered a biomarker signature of GCC. Such a finding has the potential to improve early diagnosis and prognosis and helps to identify new therapeutic targets.     

Dbf4-Cdc7 Phosphorylation of Mcm2 Is Required for Cell Growth

The Dbf4-Cdc7 kinase (DDK) is required for the activation of the origins of replication, and DDK phosphorylates Mcm2 in vitro. We find that budding yeast Cdc7 alone exists in solution as a weakly active multimer. Dbf4 forms a likely heterodimer with Cdc7, and this species phosphorylates Mcm2 with substantially higher specific activity. Dbf4 alone binds tightly to Mcm2, whereas Cdc7 alone binds weakly to Mcm2, suggesting that Dbf4 recruits Cdc7 to phosphorylate Mcm2. DDK phosphorylates two serine residues of Mcm2 near the N terminus of the protein, Ser-164 and Ser-170. Expression of mcm2-S170A is lethal to yeast cells that lack endogenous MCM2 (mcm2Δ); however, this lethality is rescued in cells harboring the DDK bypass mutant mcm5-bob1. We conclude that DDK phosphorylation of Mcm2 is required for cell growth.The Cdc7 protein kinase is required throughout the yeast S phase to activate origins (1, 2). The S phase cyclin-dependent kinase also activates yeast origins of replication (3–5). It has been proposed that Dbf4 activates Cdc7 kinase in S phase, and that Dbf4 interaction with Cdc7 is essential for Cdc7 kinase activity (6). However, it is not known how Dbf4-Cdc7 (DDK)² acts during S phase to trigger the initiation of DNA replication. DDK has homologs in other eukaryotic species, and the role of Cdc7 in activation of replication origins during S phase may be conserved (7–10).The Mcm2-7 complex functions with Cdc45 and GINS to unwind DNA at a replication fork (11–15). A mutation of MCM5 (mcm5-bob1) bypasses the cellular requirements for DBF4 and CDC7 (16), suggesting a critical physiologic interaction between Dbf4-Cdc7 and Mcm proteins. DDK phosphorylates Mcm2 in vitro with proteins purified from budding yeast (17, 18) or human cells (19). Furthermore, there are mutants of MCM2 that show synthetic lethality with DBF4 mutants (6, 17), suggesting a biologically relevant interaction between DBF4 and MCM2. Nevertheless, the physiologic role of DDK phosphorylation of Mcm2 is a matter of dispute. In human cells, replacement of MCM2 DDK-phosphoacceptor residues with alanines inhibits DNA replication, suggesting that Dbf4-Cdc7 phosphorylation of Mcm2 in humans is important for DNA replication (20). In contrast, mutation of putative DDK phosphorylation sites at the N terminus of Schizosaccharomyces pombe Mcm2 results in viable cells, suggesting that phosphorylation of S. pombe Mcm2 by DDK is not critical for cell growth (10).In budding yeast, Cdc7 is present at high levels in G₁ and S phase, whereas Dbf4 levels peak in S phase (18, 21, 22). Furthermore, budding yeast DDK binds to chromatin during S phase (6), and it has been shown that Dbf4 is required for Cdc7 binding to chromatin in budding yeast (23, 24), fission yeast (25), and Xenopus (9). Human and fission yeast Cdc7 are inert on their own (7, 8), but Dbf4-Cdc7 is active in phosphorylating Mcm proteins in budding yeast (6, 26), fission yeast (7), and human (8, 10). Based on these data, it has been proposed that Dbf4 activates Cdc7 kinase in S phase and that Dbf4 interaction with Cdc7 is essential for Cdc7 kinase activity (6, 9, 18, 21–24). However, a mechanistic analysis of how Dbf4 activates Cdc7 has not yet been accomplished. For example, the multimeric state of the active Dbf4-Cdc7 complex is currently disputed. A heterodimer of fission yeast Cdc7 (Hsk1) in complex with fission yeast Dbf4 (Dfp1) can phosphorylate Mcm2 (7). However, in budding yeast, oligomers of Cdc7 exist in the cell (27), and Dbf4-Cdc7 exists as oligomers of 180 and 300 kDa (27).DDK phosphorylates the N termini of human Mcm2 (19, 20, 28), human Mcm4 (10), budding yeast Mcm4 (26), and fission yeast Mcm6 (10). Although the sequences of the Mcm N termini are poorly conserved, the DDK sites identified in each study have neighboring acidic residues. The residues of budding yeast Mcm2 that are phosphorylated by DDK have not yet been identified.In this study, we find that budding yeast Cdc7 is weakly active as a multimer in phosphorylating Mcm2. However, a low molecular weight form of Dbf4-Cdc7, likely a heterodimer, has a higher specific activity for phosphorylation of Mcm2. Dbf4 or DDK, but not Cdc7, binds tightly to Mcm2, suggesting that Dbf4 recruits Cdc7 to Mcm2. DDK phosphorylates two serine residues of Mcm2, Ser-164 and Ser-170, in an acidic region of the protein. Mutation of Ser-170 is lethal to yeast cells, but this phenotype is rescued by the DDK bypass mutant mcm5-bob1. We conclude that DDK phosphorylation of Ser-170 of Mcm2 is required for budding yeast growth.     

High-definition De Novo Sequencing of Crustacean Hyperglycemic Hormone (CHH)-family Neuropeptides

A complete understanding of the biological functions of large signaling peptides (>4 kDa) requires comprehensive characterization of their amino acid sequences and post-translational modifications, which presents significant analytical challenges. In the past decade, there has been great success with mass spectrometry-based de novo sequencing of small neuropeptides. However, these approaches are less applicable to larger neuropeptides because of the inefficient fragmentation of peptides larger than 4 kDa and their lower endogenous abundance. The conventional proteomics approach focuses on large-scale determination of protein identities via database searching, lacking the ability for in-depth elucidation of individual amino acid residues. Here, we present a multifaceted MS approach for identification and characterization of large crustacean hyperglycemic hormone (CHH)-family neuropeptides, a class of peptide hormones that play central roles in the regulation of many important physiological processes of crustaceans. Six crustacean CHH-family neuropeptides (8–9.5 kDa), including two novel peptides with extensive disulfide linkages and PTMs, were fully sequenced without reference to genomic databases. High-definition de novo sequencing was achieved by a combination of bottom-up, off-line top-down, and on-line top-down tandem MS methods. Statistical evaluation indicated that these methods provided complementary information for sequence interpretation and increased the local identification confidence of each amino acid. Further investigations by MALDI imaging MS mapped the spatial distribution and colocalization patterns of various CHH-family neuropeptides in the neuroendocrine organs, revealing that two CHH-subfamilies are involved in distinct signaling pathways.Neuropeptides and hormones comprise a diverse class of signaling molecules involved in numerous essential physiological processes, including analgesia, reward, food intake, learning and memory (1). Disorders of the neurosecretory and neuroendocrine systems influence many pathological processes. For example, obesity results from failure of energy homeostasis in association with endocrine alterations (2, 3). Previous work from our lab used crustaceans as model organisms found that multiple neuropeptides were implicated in control of food intake, including RFamides, tachykinin related peptides, RYamides, and pyrokinins (4–6).Crustacean hyperglycemic hormone (CHH)¹ family neuropeptides play a central role in energy homeostasis of crustaceans (7–17). Hyperglycemic response of the CHHs was first reported after injection of crude eyestalk extract in crustaceans. Based on their preprohormone organization, the CHH family can be grouped into two sub-families: subfamily-I containing CHH, and subfamily-II containing molt-inhibiting hormone (MIH) and mandibular organ-inhibiting hormone (MOIH). The preprohormones of the subfamily-I have a CHH precursor related peptide (CPRP) that is cleaved off during processing; and preprohormones of the subfamily-II lack the CPRP (9). Uncovering their physiological functions will provide new insights into neuroendocrine regulation of energy homeostasis.Characterization of CHH-family neuropeptides is challenging. They are comprised of more than 70 amino acids and often contain multiple post-translational modifications (PTMs) and complex disulfide bridge connections (7). In addition, physiological concentrations of these peptide hormones are typically below picomolar level, and most crustacean species do not have available genome and proteome databases to assist MS-based sequencing.MS-based neuropeptidomics provides a powerful tool for rapid discovery and analysis of a large number of endogenous peptides from the brain and the central nervous system. Our group and others have greatly expanded the peptidomes of many model organisms (3, 18–33). For example, we have discovered more than 200 neuropeptides with several neuropeptide families consisting of as many as 20–40 members in a simple crustacean model system (5, 6, 25–31, 34). However, a majority of these neuropeptides are small peptides with 5–15 amino acid residues long, leaving a gap of identifying larger signaling peptides from organisms without sequenced genome. The observed lack of larger size peptide hormones can be attributed to the lack of effective de novo sequencing strategies for neuropeptides larger than 4 kDa, which are inherently more difficult to fragment using conventional techniques (34–37). Although classical proteomics studies examine larger proteins, these tools are limited to identification based on database searching with one or more peptides matching without complete amino acid sequence coverage (36, 38).Large populations of neuropeptides from 4–10 kDa exist in the nervous systems of both vertebrates and invertebrates (9, 39, 40). Understanding their functional roles requires sufficient molecular knowledge and a unique analytical approach. Therefore, developing effective and reliable methods for de novo sequencing of large neuropeptides at the individual amino acid residue level is an urgent gap to fill in neurobiology. In this study, we present a multifaceted MS strategy aimed at high-definition de novo sequencing and comprehensive characterization of the CHH-family neuropeptides in crustacean central nervous system. The high-definition de novo sequencing was achieved by a combination of three methods: (1) enzymatic digestion and LC-tandem mass spectrometry (MS/MS) bottom-up analysis to generate detailed sequences of proteolytic peptides; (2) off-line LC fractionation and subsequent top-down MS/MS to obtain high-quality fragmentation maps of intact peptides; and (3) on-line LC coupled to top-down MS/MS to allow rapid sequence analysis of low abundance peptides. Combining the three methods overcomes the limitations of each, and thus offers complementary and high-confidence determination of amino acid residues. We report the complete sequence analysis of six CHH-family neuropeptides including the discovery of two novel peptides. With the accurate molecular information, MALDI imaging and ion mobility MS were conducted for the first time to explore their anatomical distribution and biochemical properties.