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.     

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]     

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]     

Proteome-wide Substrate Analysis Indicates Substrate Exclusion as a Mechanism to Generate Caspase-7 Versus Caspase-3 Specificity

Caspase-3 and -7 are considered functionally redundant proteases with similar proteolytic specificities. We performed a proteome-wide screen on a mouse macrophage lysate using the N-terminal combined fractional diagonal chromatography technology and identified 46 shared, three caspase-3-specific, and six caspase-7-specific cleavage sites. Further analysis of these cleavage sites and substitution mutation experiments revealed that for certain cleavage sites a lysine at the P5 position contributes to the discrimination between caspase-7 and -3 specificity. One of the caspase-7-specific substrates, the 40 S ribosomal protein S18, was studied in detail. The RPS18-derived P6–P5′ undecapeptide retained complete specificity for caspase-7. The corresponding P6–P1 hexapeptide still displayed caspase-7 preference but lost strict specificity, suggesting that P′ residues are additionally required for caspase-7-specific cleavage. Analysis of truncated peptide mutants revealed that in the case of RPS18 the P4–P1 residues constitute the core cleavage site but that P6, P5, P2′, and P3′ residues critically contribute to caspase-7 specificity. Interestingly, specific cleavage by caspase-7 relies on excluding recognition by caspase-3 and not on increasing binding for caspase-7.Caspases, a family of evolutionarily conserved proteases, mediate apoptosis, inflammation, proliferation, and differentiation by cleaving many cellular substrates (1–3). The apoptotic initiator caspases (caspase-8, -9, and -10) are activated in large signaling platforms and propagate the death signal by cleavage-induced activation of executioner caspase-3 and -7 (4, 5). Most of the cleavage events occurring during apoptosis have been attributed to the proteolytic activity of these two executioner caspases, which can act on several hundreds of proteins (2, 3, 6, 7). The substrate degradomes of the two main executioner caspases have not been determined but their identification is important to gaining greater insight in their cleavage specificity and biological functions.The specificity of caspases was rigorously profiled by using combinatorial tetrapeptide libraries (8), proteome-derived peptide libraries (9), and sets of individual peptide substrates (10, 11). The results of these studies indicate that specificity motifs for caspase-3 and -7 are nearly indistinguishable with the canonical peptide substrate, DEVD, used to monitor the enzymatic activity of both caspase-3 and -7 in biological samples. This overlap in cleavage specificity is manifested in their generation of similar cleavage fragments from a variety of apoptosis-related substrates such as inhibitor of caspase-activated DNase, keratin 18, PARP,1 protein-disulfide isomerase, and Rho kinase I (for reviews, see Refs. 2, 3, and 7). This propagated the view that these two caspases have completely redundant functions during apoptosis. Surprisingly, mice deficient in one of these caspases (as well as mice deficient in both) have distinct phenotypes. Depending on the genetic background of the mice, caspase-3-deficient mice either die before birth (129/SvJ) or develop almost normally (C57BL/6J) (12–14). This suggests that dynamics in the genetic background, such as increased caspase-7 expression, compensate for the functional loss of caspase-3 (15). In the C57BL/6J background, caspase-7 single deficient mice are also viable, whereas caspase-3 and -7 double deficient mice die as embryos, further suggesting redundancy (12–14). However, because caspase-3 and -7 probably arose from gene duplication between the Cephalochordata-Vertebrata diversion (16), they might have acquired different substrate specificities during evolution. Caspase-3 and -7 do exhibit different activities on a few arbitrarily identified natural substrates, including BID, X-linked inhibitor of apoptosis protein, gelsolin, caspase-6, ataxin-7, and co-chaperone p23 (17–20). In addition, caspase-3 generally cleaves more substrates during apoptosis than caspase-7 and therefore appears to be the major executioner caspase. Moreover, a recent report describing caspase-1-dependent activation of caspase-7, but not of caspase-3, in macrophages in response to microbial stimuli supports the idea of a non-redundant function for caspase-7 downstream of caspase-1 (21).Commercially available “caspase-specific” tetrapeptide substrates are widely used for specific caspase detection, but they display substantial promiscuity and cannot be used to monitor individual caspases in cells (22, 23). Detecting proteolysis by measuring the release of C-terminal fluorophores, such as 7-amino-4-methylcoumarin (amc), restricts the specificity of these peptide substrates to non-prime cleavage site residues, which may have hampered the identification of specific cleavage events. To address this limitation, a recently developed proteomics technique, called proteomic identification of protease cleavage sites, was used to map both non-prime and prime preferences for caspase-3 and -7 on a tryptic peptide library (9). However, no clear distinction in peptide recognition motifs between caspase-3 and -7 could be observed (9). Because not all classical caspase cleavage sites are processed (7), structural or post-translational higher order constraints are likely involved in steering the cleavage site selectivity. Peptide-based approaches generally overlook such aspects.We made use of the COFRADIC N-terminal peptide sorting methodology (24–26) to profile proteolytic events of caspase-3 and -7 in a macrophage proteome labeled by triple stable isotope labeling by amino acids in cell culture (SILAC), which allowed direct comparison of peak intensities in peptide MS spectra and consequent quantification of N termini that are equally, preferably, or exclusively generated by the action of caspase-3 or -7 (26, 27). We identified 55 cleavage sites in 48 protein substrates, encompassing mutual, preferred, and unique caspase-3 and -7 cleavage sites.     

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]