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.     

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]     

Identification of Lysine Succinylation Substrates and the Succinylation Regulatory Enzyme CobB in Escherichia coli

Lysine succinylation is a newly identified protein post-translational modification pathway present in both prokaryotic and eukaryotic cells. However, succinylation substrates and regulatory enzyme(s) remain largely unknown, hindering the biological study of this modification. Here we report the identification of 2,580 bacterial lysine succinylation sites in 670 proteins and 2,803 lysine acetylation (Kac) sites in 782 proteins, representing the first lysine succinylation dataset and the largest Kac dataset in wild-type E. coli. We quantified dynamic changes of the lysine succinylation and Kac substrates in response to high glucose. Our data showed that high-glucose conditions led to more lysine-succinylated proteins and enhanced the abundance of succinyllysine peptides more significantly than Kac peptides, suggesting that glucose has a more profound effect on succinylation than on acetylation. We further identified CobB, a known Sir2-like bacterial lysine deacetylase, as the first prokaryotic desuccinylation enzyme. The identification of bacterial CobB as a bifunctional enzyme with lysine desuccinylation and deacetylation activities suggests that the eukaryotic Kac-regulatory enzymes may have enzymatic activities on various lysine acylations with very different structures. In addition, it is highly likely that lysine succinylation could have unique and more profound regulatory roles in cellular metabolism relative to lysine acetylation under some physiological conditions.Lysine acetylation (Kac)¹ is a dynamic and evolutionarily conserved post-translational modification (PTM) that is known to be involved in the regulation of diverse cellular processes (1–9). The status of this modification is controlled by two groups of enzymes with opposing enzymatic activities, lysine acetyltransferases that add an acetyl group to the lysine (Lys or K) residue, and histone lysine deacetylases (HDACs) that remove the acetyl group (10–16). HDACs are grouped into several categories (17): class I (HDAC1, -2, -3, and -8), class IIA (HDAC4, -5, -7, and -9), class IIB (HDAC6 and -10), class III (Sirt1–7), and class IV (HDAC11). The weak deacetylation activities of some HDACs (e.g. Sirt4–7 and HDAC4, -5, and -7–11), as well as the demonstration of Sirt5 as a desuccinylation and demalonylation enzyme, suggest that some HDAC enzymes have activities that are independent of acetylation (18, 19).For a long period of time, lysine acetylation was considered as a protein modification that was restricted to nuclei (20). The identification of cytosolic Kac substrates and the localization of some HDACs outside nuclei suggest a non-nuclear function of lysine acetylation (13, 21, 22). The first proteomic screening identified hundreds of substrate proteins in cytosolic and mitochondrial fractions and demonstrated high abundance of Kac in mitochondrial proteins and metabolic enzymes (23). This result implies that Kac has diverse non-nuclear roles and can regulate functions of metabolism and mitochondria (23). Since then, we and others have extensively characterized the cellular acetylome (5, 9, 24–26).The lysine succinylation (Ksucc) and lysine malonylation pathways are two PTM pathways that were recently identified and comprehensively validated in both bacterial and mammalian cells, with multiple substrate proteins identified, using HPLC-MS/MS, co-elution of synthetic peptides, isotopic labeling, Western blotting analysis using pan-anti-Ksucc antibodies, and proteomics analysis (18, 27). We also showed that Ksucc is present in core histones (29). In yeast histones, some Ksucc sites are located in regions where histones make close contact with DNA, suggesting that Ksucc sites may be involved in gene regulation by changing the chromatin structure (29). We then found that Sirt5, a member of the class III family of HDACs, can function as a desuccinylation enzyme in vitro and in vivo (18, 19). In a recent study, we revealed that Sirt5 is a key regulatory enzyme of Ksucc and that Ksucc proteins are abundant among a group of mitochondrial enzymes that are predominantly involved in fatty acid metabolism, amino acid degradation, and the tricarboxylic acid cycle (28). Importantly, Ksucc is very dynamic not only in mammalian cells, but also in bacteria (27, 29). These lines of evidence strongly suggest that lysine succinylation is likely an important PTM in the regulation of cellular functions.Although key elements of the Ksucc pathway are being identified in mammalian cells, their counterparts in bacteria remain largely unknown. We and others have used a proteomics approach to identify Kac substrates in bacteria (26, 30, 31, 52). The Sir2-like enzyme CobB is the best-studied protein deacetylase in bacteria (8). CobB was initially identified as an enzyme required for the activation of acetyl-CoA synthetase (8). Recently, CobB was shown to play roles in bacterial energy metabolism (31) and stress response (32). Those studies indicated that Kac is an evolutionarily conserved PTM with a role in energy metabolism in prokaryotes. Nevertheless, dynamic changes of lysine acetylation in bacteria have not been studied. In addition, substrates of lysine succinylation and their regulatory enzymes are not known.In this paper, we report a quantitative proteomic approach based on stable isotope labeling by amino acids in cell culture (SILAC) to identify and quantify changes in bacterial lysine succinylation, as well as lysine acetylation, in response to glucose, a major energy source. Our screening detected 2,580 lysine-succinylated sites in 670 proteins and 2,803 Kac sites in 782 proteins in Escherichia coli. Our quantitative proteomics data show that glucose had a more profound effect on Ksucc than on Kac. In addition, we found that CobB, a known prokaryotic deacetylase, had dual enzymatic activities to catalyze the removal of two structurally different lysine acyl groups, acetyl and succinyl, from the modified lysine residues.     

Algorithms for Finding Small Attractors in Boolean Networks

A Boolean network is a model used to study the interactions between different genes in genetic regulatory networks. In this paper, we present several algorithms using gene ordering and feedback vertex sets to identify singleton attractors and small attractors in Boolean networks. We analyze the average case time complexities of some of the proposed algorithms. For instance, it is shown that the outdegree-based ordering algorithm for finding singleton attractors works in time for , which is much faster than the naive time algorithm, where is the number of genes and is the maximum indegree. We performed extensive computational experiments on these algorithms, which resulted in good agreement with theoretical results. In contrast, we give a simple and complete proof for showing that finding an attractor with the shortest period is NP-hard.[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]     

Proinsulin and the Genetics of Diabetes Mellitus

Insulin plays a central role in the regulation of vertebrate metabolism. The hormone, the post-translational product of a single-chain precursor, is a globular protein containing two chains, A (21 residues) and B (30 residues). Recent advances in human genetics have identified dominant mutations in the insulin gene causing permanent neonatal-onset DM² (1–4). The mutations are predicted to block folding of the precursor in the ER of pancreatic β-cells. Although expression of the wild-type allele would in other circumstances be sufficient to maintain homeostasis, studies of a corresponding mouse model (5–7) suggest that the misfolded variant perturbs wild-type biosynthesis (8, 9). Impaired β-cell secretion is associated with ER stress, distorted organelle architecture, and cell death (10). These findings have renewed interest in insulin biosynthesis (11–13) and the structural basis of disulfide pairing (14–19). Protein evolution is constrained not only by structure and function but also by susceptibility to toxic misfolding.Insulin plays a central role in the regulation of vertebrate metabolism. The hormone, the post-translational product of a single-chain precursor, is a globular protein containing two chains, A (21 residues) and B (30 residues). Recent advances in human genetics have identified dominant mutations in the insulin gene causing permanent neonatal-onset DM² (1–4). The mutations are predicted to block folding of the precursor in the ER of pancreatic β-cells. Although expression of the wild-type allele would in other circumstances be sufficient to maintain homeostasis, studies of a corresponding mouse model (5–7) suggest that the misfolded variant perturbs wild-type biosynthesis (8, 9). Impaired β-cell secretion is associated with ER stress, distorted organelle architecture, and cell death (10). These findings have renewed interest in insulin biosynthesis (11–13) and the structural basis of disulfide pairing (14–19). Protein evolution is constrained not only by structure and function but also by susceptibility to toxic misfolding.     

Deep,Quantitative Coverage of the Lysine Acetylome Using Novel Anti-acetyl-lysine Antibodies and an Optimized Proteomic Workflow

Introduction of antibodies specific for acetylated lysine has significantly improved the detection of endogenous acetylation sites by mass spectrometry. Here, we describe a new, commercially available mixture of anti-lysine acetylation (Kac) antibodies and show its utility for in-depth profiling of the acetylome. Specifically, seven complementary monoclones with high specificity for Kac were combined into a final anti-Kac reagent which results in at least a twofold increase in identification of Kac peptides over a commonly used Kac antibody. We outline optimal antibody usage conditions, effective offline basic reversed phase separation, and use of state-of-the-art LC-MS technology for achieving unprecedented coverage of the acetylome. The methods were applied to quantify acetylation sites in suberoylanilide hydroxamic acid-treated Jurkat cells. Over 10,000 Kac peptides from over 3000 Kac proteins were quantified from a single stable isotope labeling by amino acids in cell culture labeled sample using 7.5 mg of peptide input per state. This constitutes the deepest coverage of acetylation sites in quantitative experiments obtained to-date. The approach was also applied to breast tumor xenograft samples using isobaric mass tag labeling of peptides (iTRAQ4, TMT6 and TMT10-plex reagents) for quantification. Greater than 6700 Kac peptides from over 2300 Kac proteins were quantified using 1 mg of tumor protein per iTRAQ 4-plex channel. The novel reagents and methods we describe here enable quantitative, global acetylome analyses with depth and sensitivity approaching that obtained for other well-studied post-translational modifications such as phosphorylation and ubiquitylation, and should have widespread application in biological and clinical studies employing mass spectrometry-based proteomics.Lysine acetylation (Kac)¹ is a well conserved, reversible post-translational modification (PTM) involved in multiple cellular processes (1). Acetylation is regulated by two classes of enzymes: lysine acetyltransferases (KATs) and histone deacetylases (HDACs) (2–4). This modification was originally identified as a nuclear event on histone proteins and has been long appreciated for its role in epigenetic and DNA-dependent processes. With the help of a growing number of large-scale acetylation studies, it has become evident that lysine acetylation is ubiquitous, also occurring on cytoplasmic and mitochondrial proteins and has a role in signaling, metabolism, and immunity (1, 4–6). Therefore, the examination of lysine acetylation on nonhistone proteins has gained a prominent role in PTM analysis.To date, the identification of large numbers of acetylation sites has been challenging because of the substoichiometric nature of this modification (7, 8). Additionally, global acetylation is generally less abundant than phosphorylation and ubiquitylation (1). The introduction of antibodies specific for lysine acetylation has significantly improved the ability to enrich and identify thousands of sites (9–14). A landmark study by Choudhary et al. used anti-Kac antibodies to globally map 3600 lysine acetylation sites on 1750 proteins, thereby demonstrating the feasibility of profiling the acetylome (10). A more recent study by Lundby et al. investigated the function and distribution of acetylation sites in 16 different rat tissues, and identified, in aggregate, 15,474 acetylation sites from 4541 proteins (12).Although anti-acetyl lysine antibodies have been a breakthrough for globally mapping acetylation sites (9–12), it remains a challenge to identify large numbers of lysine acetylation sites from a single sample, as is now routinely possible for phosphorylation and ubiquitylation (13, 15–18). To improve the depth-of-coverage in acetylation profiling experiments there is a clear need for (1) alternative anti-acetyl lysine antibodies with higher specificity, (2) optimized antibody usage parameters, and (3) robust proteomic workflows that permit low to moderate protein input. In this study, we describe a newly commercialized mixture of anti-Kac antibodies and detail a complete proteomic workflow for achieving unprecedented coverage of the acetylome from a single stable isotope labeling by amino acids in cell culture (SILAC) labeled sample as well as isobaric tags for relative and absolute quantitation (iTRAQ)- and tandem mass tag (TMT)-labeled samples.     

The Wavelet-Based Cluster Analysis for Temporal Gene Expression Data

A variety of high-throughput methods have made it possible to generate detailed temporal expression data for a single gene or large numbers of genes. Common methods for analysis of these large data sets can be problematic. One challenge is the comparison of temporal expression data obtained from different growth conditions where the patterns of expression may be shifted in time. We propose the use of wavelet analysis to transform the data obtained under different growth conditions to permit comparison of expression patterns from experiments that have time shifts or delays. We demonstrate this approach using detailed temporal data for a single bacterial gene obtained under 72 different growth conditions. This general strategy can be applied in the analysis of data sets of thousands of genes under different conditions.[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]     

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]     

Multipattern Consensus Regions in Multiple Aligned Protein Sequences and Their Segmentation

Decomposing a biological sequence into its functional regions is an important prerequisite to understand the molecule. Using the multiple alignments of the sequences, we evaluate a segmentation based on the type of statistical variation pattern from each of the aligned sites. To describe such a more general pattern, we introduce multipattern consensus regions as segmented regions based on conserved as well as interdependent patterns. Thus the proposed consensus region considers patterns that are statistically significant and extends a local neighborhood. To show its relevance in protein sequence analysis, a cancer suppressor gene called p53 is examined. The results show significant associations between the detected regions and tendency of mutations, location on the 3D structure, and cancer hereditable factors that can be inferred from human twin studies.[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]