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.     

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.     

Proteome-wide Analysis of Lysine Acetylation Suggests its Broad Regulatory Scope in Saccharomyces cerevisiae

Post-translational modification of proteins by lysine acetylation plays important regulatory roles in living cells. The budding yeast Saccharomyces cerevisiae is a widely used unicellular eukaryotic model organism in biomedical research. S. cerevisiae contains several evolutionary conserved lysine acetyltransferases and deacetylases. However, only a few dozen acetylation sites in S. cerevisiae are known, presenting a major obstacle for further understanding the regulatory roles of acetylation in this organism. Here we use high resolution mass spectrometry to identify about 4000 lysine acetylation sites in S. cerevisiae. Acetylated proteins are implicated in the regulation of diverse cytoplasmic and nuclear processes including chromatin organization, mitochondrial metabolism, and protein synthesis. Bioinformatic analysis of yeast acetylation sites shows that acetylated lysines are significantly more conserved compared with nonacetylated lysines. A large fraction of the conserved acetylation sites are present on proteins involved in cellular metabolism, protein synthesis, and protein folding. Furthermore, quantification of the Rpd3-regulated acetylation sites identified several previously known, as well as new putative substrates of this deacetylase. Rpd3 deficiency increased acetylation of the SAGA (Spt-Ada-Gcn5-Acetyltransferase) complex subunit Sgf73 on K33. This acetylation site is located within a critical regulatory domain in Sgf73 that interacts with Ubp8 and is involved in the activation of the Ubp8-containing histone H2B deubiquitylase complex. Our data provides the first global survey of acetylation in budding yeast, and suggests a wide-ranging regulatory scope of this modification. The provided dataset may serve as an important resource for the functional analysis of lysine acetylation in eukaryotes.Lysine acetylation is a dynamic and reversible post-translational modification. Acetylation of lysines on their ε-amino group is catalyzed by lysine acetyltransferases (KATs¹, also known as histone acetyltrasferases (HATs)), and reversed by lysine deacetylases (KDACs, also known as histone deacetylases (HDACs)) (1). The enzymatic machinery involved in lysine acetylation is evolutionary conserved in all forms of life (2–4). The role of acetylation has been extensively studied in the regulation of gene expression via modification of histones (5). Acetylation also plays important roles in controlling cellular metabolism (6–10), protein folding (11), and sister chromatid cohesion (12). Furthermore, acetylation has been implicated in regulating the beneficial effects of calorie restriction (13), a low nutrient diet without starvation, and aging. Based on these findings, it is proposed that the functional roles of acetylation in these processes are evolutionary conserved from yeast to mammals.Advancements in mass spectrometry (MS)-based proteomics have greatly facilitated identification of thousands of post-translational modification (PTM) sites in eukaryotic cells (14–18). Proteome-wide mapping of PTM sites can provide important leads for analyzing the functional relevance of individual sites and a systems-wide view of the regulatory scope of post-translational modifications. Also, large-scale PTM datasets are an important resource for the in silico analysis of PTMs, which can broaden the understanding of modification site properties and their evolutionary trajectories.The budding yeast Saccharomyces cerevisiae is a commonly used unicellular eukaryotic model organism. Yeast has been used in many pioneering “-omics” studies, including sequencing of the first eukaryotic genome (19), systems-wide genetic interactions analysis (20, 21), MS-based comprehensive mapping of a eukaryotic proteome (22), and proteome-wide analysis of protein-protein interactions (23, 24). In addition, S. cerevisiae has been extensively used to study the molecular mechanisms of acetylation. Many lysine acetyltransferases and deacetylases were discovered in this organism (2, 25), and their orthologs were subsequently identified in higher eukaryotes. Furthermore, the functional roles of many well-studied acetylation sites on histones are conserved from yeast to mammals. Recent data from human and Drosophila cells show that acetylation is present on many highly conserved metabolic enzymes (26–28). However, only a few dozen yeast acetylation sites are annotated in the Uniprot database. Given the presence of a well-conserved and elaborate acetylation machinery in yeast, we reasoned that many more acetylation sites exist in this organism that remained to be identified.Here we used high resolution mass spectrometry-based proteomics to investigate the scope of acetylation in S. cerevisiae. We identified about 4000 unique acetylation sites in this important model organism. Bioinformatic analysis of yeast acetylation sites and comparison with previously identified human and Drosophila acetylation sites indicates that many acetylation sites are evolutionary conserved. Furthermore, quantitative analysis of the Rpd3-regulated acetylation sites identified several nuclear proteins that showed increased acetylation in rpd3 knockout cells. Our results provide a systems-wide view of acetylation in budding yeast, and a rich dataset for functional analysis of acetylation sites in this organism.     

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]