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]     

Inflammasome-dependent Caspase-1 Activation in Cervical Epithelial Cells Stimulates Growth of the Intracellular Pathogen Chlamydia trachomatis

Inflammasomes have been extensively characterized in monocytes and macrophages, but not in epithelial cells, which are the preferred host cells for many pathogens. Here we show that cervical epithelial cells express a functional inflammasome. Infection of the cells by Chlamydia trachomatis leads to activation of caspase-1, through a process requiring the NOD-like receptor family member NLRP3 and the inflammasome adaptor protein ASC. Secretion of newly synthesized virulence proteins from the chlamydial vacuole through a type III secretion apparatus results in efflux of K⁺ through glibenclamide-sensitive K⁺ channels, which in turn stimulates production of reactive oxygen species. Elevated levels of reactive oxygen species are responsible for NLRP3-dependent caspase-1 activation in the infected cells. In monocytes and macrophages, caspase-1 is involved in processing and secretion of pro-inflammatory cytokines such as interleukin-1β. However, in epithelial cells, which are not known to secrete large quantities of interleukin-1β, caspase-1 has been shown previously to enhance lipid metabolism. Here we show that, in cervical epithelial cells, caspase-1 activation is required for optimal growth of the intracellular chlamydiae.Chlamydia trachomatis is the most common cause of bacterial sexually transmitted disease in the United States, and it is the leading cause of preventable blindness in the world (1–5). Untreated, C. trachomatis infection in women can cause pelvic inflammatory disease, which can lead to infertility and ectopic pregnancy because of scarring of the ovaries and the Fallopian tubes (6). Infection by the lymphogranuloma venereum (LGV)² strain of C. trachomatis, which has become more common in North America and Europe (7, 8), is characterized by swelling and inflammation of the lymph nodes in the groin (9).Chlamydiae are intracellular pathogens that preferentially infect epithelial mucosa and have a biphasic infection cycle (10). A metabolically inactive form, the elementary body, infects the epithelial host cells through entry vesicles that avoid fusion with host cell lysosomes and develop into a membrane-bound inclusion (11–13). Despite their intravacuolar localization, chlamydiae are still able to acquire nutrients from the host cell and interact with host-cell signaling pathways (13–23). Within a few hours, the elementary bodies differentiate into larger, metabolically active reticulate bodies, which proliferate but are noninfectious. Depending on the strain of C. trachomatis, the reticulate bodies transform back into elementary bodies after 1–3 days and are released into the extracellular medium to infect other cells (11, 24, 25). Chlamydial species possess a type III secretion (T3S) system that secretes bacterial virulence factors into host cell cytosol and may control interactions between the inclusion and host-cell compartments (26).Long before the adaptive immune response is activated, infected epithelial cells produce proinflammatory cytokines and chemokines, including interleukin (IL)-6, IL-8, and granulocyte-macrophage colony-stimulating factor (27), which recruit neutrophils to the site of infection and activate other immune effector cells. However, in many cases the immune system fails to clear the infection, and the chronic release of cytokines becomes a major contributor to the scarring and damage associated with the infection (28–30).The innate immune response during C. trachomatis infection is initiated by chlamydial pathogen-associated molecular patterns, including lipopolysaccharides, which bind to pattern recognition receptors such as Toll-like receptors and cytosolic NOD-like receptors (NLRs), ultimately promoting pro-inflammatory cytokine gene expression and secretion of the cytokine proteins (31–37). However, secretion of the key pro-inflammatory cytokine IL-1β is tightly regulated (38). First, pro-IL-1β is produced following activation of pattern recognition receptor, and the precursor is then cleaved into the mature form by the pro-inflammatory cysteine protease, caspase-1 (also known as interleukin-1 converting enzyme or ICE). The mechanism by which caspase-1 is activated in response to infection or tissue damage was found to be modulated by a macromolecular protein complex termed the “inflammasome,” which consists of an NLR family member, an adaptor protein (apoptosis-associated speck-like protein containing a caspase activation recruitment domain or ASC), and an inactive caspase-1 precursor (pro-caspase-1) (39, 40). Previous studies demonstrated that IL-1β is produced in response to chlamydial infection in dendritic cells, macrophages, and monocytes (41–44). Moreover, C. trachomatis or Chlamydia caviae infection activates caspase-1 in epithelial cells or monocytes (43, 45, 46). However, whether caspase-1 activation during chlamydial infection requires the formation of an inflammasome remains unclear.Previous studies have shown that different pathogens can cause inflammasome-mediated caspase-1 activation in macrophages and monocytes (47). However, epithelial cells lining mucosal surfaces are not only the preferred target for chlamydial infection and other intracellular pathogens but also play an important role in early host immune response to infection by secreting proinflammatory cytokines and chemokines (27). Although epithelial cells are not known to secrete large amounts of IL-1β, inflammasome-dependent caspase-1 activation in epithelial cells is known to contribute to lipid metabolism and membrane regeneration in epithelial cells damaged by the membrane-disrupting toxin, aerolysin (48). As lipids are sorted from the Golgi apparatus to the chlamydial inclusion (13, 15, 49), we therefore investigated whether C. trachomatis induces caspase-1 activation in epithelial cells via the assembly of an inflammasome. We demonstrated that C. trachomatis-induced caspase-1 activation is mediated by an inflammasome containing the NLR member, NLRP3. Several studies have demonstrated the involvement of T3S apparatus in inflammasome-mediated caspase-1 activation by different pathogens in macrophages and monocytes (50–56). Therefore, we further investigated the mechanism by which C. trachomatis triggers the formation of the NLRP3 inflammasome. Our results showed that metabolically active chlamydiae, relying on their T3S apparatus, cause K⁺ efflux, which in turn leads to formation of reactive oxygen species (ROS) and ultimately NLRP3-dependent caspase-1 activation. Epithelial cells do not typically secrete large amounts of IL-1β; instead, caspase-1 activation in cervical epithelial cells contributes to development of the chlamydial inclusion.     

Functions of Manduca sexta Hemolymph Proteinases HP6 and HP8 in Two Innate Immune Pathways

Serine proteinases in insect plasma have been implicated in two types of immune responses; that is, activation of prophenoloxidase (proPO) and activation of cytokine-like proteins. We have identified more than 20 serine proteinases in hemolymph of the tobacco hornworm, Manduca sexta, but functions are known for only a few of them. We report here functions of two additional M. sexta proteinases, hemolymph proteinases 6 and 8 (HP6 and HP8). HP6 and HP8 are each composed of an amino-terminal clip domain and a carboxyl-terminal proteinase domain. HP6 is an apparent ortholog of Drosophila Persephone, whereas HP8 is most similar to Drosophila and Tenebrio spätzle-activating enzymes, all of which activate the Toll pathway. proHP6 and proHP8 are expressed constitutively in fat body and hemocytes and secreted into plasma, where they are activated by proteolytic cleavage in response to infection. To investigate activation and biological activity of HP6 and HP8, we purified recombinant proHP8, proHP6, and mutants of proHP6 in which the catalytic serine was replaced with alanine, and/or the activation site was changed to permit activation by bovine factor Xa. HP6 was found to activate proPO-activating proteinase (proPAP1) in vitro and induce proPO activation in plasma. HP6 was also determined to activate proHP8. Active HP6 or HP8 injected into larvae induced expression of antimicrobial peptides and proteins, including attacin, cecropin, gloverin, moricin, and lysozyme. Our results suggest that proHP6 becomes activated in response to microbial infection and participates in two immune pathways; activation of PAP1, which leads to proPO activation and melanin synthesis, and activation of HP8, which stimulates a Toll-like pathway.Innate immune systems of mammals and arthropods include extracellular serine proteinase cascade pathways, which rapidly amplify responses to infection and stimulate killing of pathogens. These proteinase-driven processes include the complement system of vertebrates (1, 2) and pathways in arthropods involving proteinases containing amino-terminal clip domains (3). Clip domain proteinases function in blood coagulation (4, 5), activation of prophenoloxidase (proPO) that leads to melanin synthesis (6–9), and stimulation of the Toll pathway to promote synthesis of antimicrobial peptides/proteins (AMPs)² secreted into the hemolymph (10, 11).The serine proteinase systems best characterized in arthropods are the horseshoe crab hemolymph coagulation pathway and the cascade leading to activation of the Toll pathway in dorsal-ventral development in Drosophila (12–14). Recent research also has led to better characterization of the proPO activation pathway in Manduca sexta (7, 15, 16) and the Toll-signaling pathway in the Drosophila immune response (17, 18) and to both the proPO and Toll pathways in the beetle Tenebrio molitor (11, 19).In the proPO activation pathway, soluble pattern recognition proteins initially recognize pathogen-associated molecular patterns such as bacterial peptidoglycan or fungal β-1,3-glucan (20–22). This interaction stimulates the sequential activation of a series of serine proteinases in hemolymph, leading to the activation of proPO-activating proteinase (PAP), also known as proPO activating enzyme (7, 23). Activated PAP converts inactive proPO to PO. PO catalyzes the hydroxylation of monophenols to o-diphenols and the oxidation of o-diphenols to quinones that are involved in microbial killing, melanin synthesis, sequestration of parasites or pathogens, and wound healing (24, 25). Other proteins required for proPO activation are clip-domain serine proteinase homologs (SPHs), whose catalytic serine is replaced with glycine and, therefore, lack proteolytic activity (26, 27). Serine proteinase inhibitors, including members of the serpin superfamily, regulate the activation of proPO by inhibiting the activating proteinases (28, 29).Drosophila clip-domain serine proteinases Persephone, Grass, Spirit, and spätzle-processing enzyme (SPE) participate in the activation of Toll pathway, stimulating synthesis of antimicrobial peptides as an innate immune response (18, 30–32). Although genetic evidence indicates that Persephone and Spirit are upstream of SPE in the cascade, the substrate(s) of Persephone and Spirit have not been identified, and which proteinase directly activates SPE is unknown. Neither is it clear whether these enzymes may be related to the melanization pathway, which involves clip-domain proteinases MP2 and MP1 (33).Here we report the functional characterization of M. sexta HP6 and HP8, probable orthologs of Drosophila Persephone and SPE, respectively. We developed methods to activate purified recombinant proHP6 and proHP8 and discovered that HP6 participates in proPO activation by activating proPAP1 and that both HP6 and HP8 function in a pathway that stimulates the synthesis of AMPs in M. sexta.     

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]     

The Kinase Activity of Rip2 Determines Its Stability and Consequently Nod1- and Nod2-mediated Immune Responses

Rip2 (RICK, CARD3) has been identified as a key effector molecule downstream of the pattern recognition receptors, Nod1 and Nod2; however, its mechanism of action remains to be elucidated. In particular, it is unclear whether its kinase activity is required for signaling or for maintaining protein stability. We have investigated the expression level of different retrovirally expressed kinase-dead Rip2 mutants and the role of Rip2 kinase activity in the signaling events that follow Nod1 and Nod2 stimulation. We show that in primary cells expressing kinase-inactive Rip2, protein levels were severely compromised, and stability could not be reconstituted by the addition of a phospho-mimetic mutation in its autophosphorylation site. Consequently, inflammatory cytokine production in response to Nod1 and Nod2 ligands was abrogated both in vitro and in vivo in the absence of Rip2 kinase activity. Our results highlight the central role that Rip2 kinase activity plays in conferring stability to the protein and thus in the preservation of Nod1- and Nod2-mediated innate immune responses.A key step in the initiation of effector immune responses is the recognition of highly conserved molecules expressed by microbial pathogens. The immune system has developed specific receptors that sense these so-called pathogen-associated molecular patterns and initiate appropriate immune responses. One key family of pattern recognition receptors is the Nod-like receptor (NLR)² family (1–3), of which two members, Nod1 and Nod2, have been implicated in the recognition of bacterial peptidoglycan derivatives released into the cytosol upon bacterial infection (4–6). Several studies have shown that Nod1 plays a role in host defense against invasive pathogens such as Helicobacter pylori and Escherichia coli (7, 8), and Nod2 mutations have been associated with a higher incidence of Crohn disease (9, 10), thus highlighting these NLRs as important regulators of inflammatory immune responses.Rip2, also called CARD3, RICK, or CARDIAK, is a serine/threonine kinase, which was implicated in the induction of NF-κB activation and apoptosis (11–13). Rip2 has been described to be critical for responses against Toll-like receptor ligands such as LPS (14, 15), although findings from recent studies did not support this conclusion (16). Rip2 contains a caspase-recruitment domain (CARD), which mediates interaction with other CARD-containing proteins such as Nod1 and Nod2, in addition to an N-terminal kinase domain and an intermediate domain. Nod1 and Nod2 associate with Rip2 upon peptidoglycan ligation (17) leading to downstream signaling events that culminate in NF-κB and mitogen-activated protein kinase activation (15, 18–20). Recent reports have suggested that the mitogen-activated protein kinase kinase kinase family member TAK1 provides the link between Rip2 and NF-κB activation upon Nod1 and Nod2 stimulation (21–23). However, the exact role of Rip2 and in particular its kinase activity in mediating downstream effector activation in NLR signaling still remains unclear. Notably, in vitro investigations have suggested that Rip2 kinase activity may be dispensable for the induction of immune responses initiated by NLR-ligands (21, 24, 25) and that disruption of Rip2 kinase activity is associated with a loss in protein stability (23); however, such studies utilized protein overexpression in cell lines and are yet to be tested in primary cells or in vivo.In the current investigation we sought to elucidate the role of Rip2 kinase activity in transducing inflammatory signals upon NLR stimulation in vitro and in vivo. To this end, we utilized both Rip2 knock-out (15) and Rip2 kinase-dead knock-in mice (24) in addition to Rip2 deficient primary cells that were retrovirally reconstituted with different kinase-inactive mutants. We show here that in the absence of intact kinase activity, Rip2 protein is not stable and that insertion of a phospho-mimetic mutation is not sufficient to restore stability. Moreover, pharmacological abrogation of Rip2 kinase activity in primary cells similarly leads to destabilization of the molecule. As a consequence, signaling downstream of Nod1 and Nod2 and inflammatory cytokine production is impaired both in vivo and in vitro. Our results highlight Rip2 kinase activity as a central regulator of protein stability and consequently innate immune responses triggered by Nod1 and Nod2 ligands.