Lysine catabolism reprograms tumour immunity through histone crotonylation.

Cancer cells rewire metabolism to favour the generation of specialized metabolites that support tumour growth and reshape the tumour microenvironment1,2. Lysine functions as a biosynthetic molecule, energy source and antioxidant3-5, but little is known about its pathological role in cancer. Here we show that glioblastoma stem cells (GSCs) reprogram lysine catabolism through the upregulation of lysine transporter SLC7A2 and crotonyl-coenzyme A (crotonyl-CoA)-producing enzyme glutaryl-CoA dehydrogenase (GCDH) with downregulation of the crotonyl-CoA hydratase enoyl-CoA hydratase short chain 1 (ECHS1), leading to accumulation of intracellular crotonyl-CoA and histone H4 lysine crotonylation. A reduction in histone lysine crotonylation by either genetic manipulation or lysine restriction impaired tumour growth. In the nucleus, GCDH interacts with the crotonyltransferase CBP to promote histone lysine crotonylation. Loss of histone lysine crotonylation promotes immunogenic cytosolic double-stranded RNA (dsRNA) and dsDNA generation through enhanced H3K27ac, which stimulates the RNA sensor MDA5 and DNA sensor cyclic GMP-AMP synthase (cGAS) to boost type I interferon signalling, leading to compromised GSC tumorigenic potential and elevated CD8+ T cell infiltration. A lysine-restricted diet synergized with MYC inhibition or anti-PD-1 therapy to slow tumour growth. Collectively, GSCs co-opt lysine uptake and degradation to shunt the production of crotonyl-CoA, remodelling the chromatin landscape to evade interferon-induced intrinsic effects on GSC maintenance and extrinsic effects on immune response.

Lysine L-lactylation is the dominant lactylation isomer induced by glycolysis.

Lysine L-lactylation (Kl-la) is a novel protein posttranslational modification (PTM) driven by L-lactate. This PTM has three isomers: Kl-la, N-ε-(carboxyethyl)-lysine (Kce) and D-lactyl-lysine (Kd-la), which are often confused in the context of the Warburg effect and nuclear presence. Here we introduce two methods to differentiate these isomers: a chemical derivatization and high-performance liquid chromatography analysis for efficient separation, and isomer-specific antibodies for high-selectivity identification. We demonstrated that Kl-la is the primary lactylation isomer on histones and dynamically regulated by glycolysis, not Kd-la or Kce, which are observed when the glyoxalase system was incomplete. The study also reveals that lactyl-coenzyme A, a precursor in L-lactylation, correlates positively with Kl-la levels. This work not only provides a methodology for distinguishing other PTM isomers, but also highlights Kl-la as the primary responder to glycolysis and the Warburg effect.

Quantifying Turnover Dynamics of Selenoproteome by Isotopic Perturbation.

Selenium, as an essential trace element of life, is closely related to human health and is required to produce selenoproteins, a family of important functional proteins in many living organisms. All selenoproteins contain a special amino acid, selenocysteine, which often serves as their active-site residue, and the expression and activity of selenoproteins are fine-tuned. However, the turnover dynamics of selenoproteome has never been systematically investigated, especially in a site-specific manner for selenocysteines. In the current work, we developed a chemical proteomic strategy named "SElenoprotein Turnover Rate by Isotope Perturbation (SETRIP)" to quantitatively monitor the turnover dynamics of selenoproteins at the proteomic level. The kinetic rates and half-lives of nine selenoproteins were accurately measured by combining Na274SeO3 metabolic labeling with pulse-chase chemoproteomics. The half-lives of selenoproteins were measured to range from 6 to 32 h with the housekeeping selenoprotein glutathione peroxidases (GPX4) showing a faster turnover rate, implying that the hierarchy regulation also exists in the turnover of selenoproteins in addition to expression and activity. Our study generated a global portrait of dynamic changes in the selenoproteome and provided important clues to study the roles of selenium in biology.

CIMAGE2.0: An Expanded Tool for Quantitative Analysis of Activity-Based Protein Profiling (ABPP) Data.

Activity-based protein profiling (ABPP) is a powerful chemical proteomic method for studying protein activity, modifications, and interactions in a high-throughput manner. In ABPP experiments, accurate quantification is crucial to determine the extent of probe labeling at the level of either target proteins or specific amino acid side chains. CIMAGE has been developed as an in-house quantification software specifically designed for ABPP data analysis that incorporates (1) a relaxed peak extraction algorithm and (2) stringent post-quantification checks for efficient and accurate quantification. It also can generate table and image data for users to conveniently visualize their results. Here we provide a retrospective introduction of the software and describe our recent upgrade efforts to enable (1) interfacing with different database search engines as input, (2) triplex quantification of ABPP data by reductive dimethylation, and (3) envelope checking for chemical elements with special isotopic distributions. We show that the updated CIMAGE can maintain its ability to quantify ABPP data with dramatic depth and high accuracy, and it also has similar quantification performance in benchmarked SILAC data as compared with MaxQuant. We believe that CIMAGE2.0 will continue to serve as a powerful analytical tool for ABPP studies.

Quantitative Profiling of Protein Carbonylations in Ferroptosis by an Aniline-Derived Probe.

Ferroptosis is a regulated form of necrotic cell death implicated in carcinogenesis and neurodegeneration that is driven by phospholipid peroxidation. Lipid-derived electrophiles (LDEs) generated during this process can covalently modify proteins ("carbonylation") and affect their functions. Here we report the development of a quantitative chemoproteomic method to profile carbonylations in ferroptosis by an aniline-derived probe. Using the method, we established a global portrait of protein carbonylations in ferroptosis with >400 endogenously modified proteins and for the first time, identified >20 residue sites with endogenous LDE modifications in ferroptotic cells. Specifically, we discovered and validated a novel cysteine site of modification on voltage-dependent anion-selective channel protein 2 (VDAC2) that might play an important role in sensitizing LDE signals and mediating ferroptosis. Our results will contribute to the understanding of ferroptotic signaling and pathogenesis and provide potential biomarkers for ferroptosis detection.

A Dimethyl-Labeling-Based Strategy for Site-Specifically Quantitative Chemical Proteomics.

Activity-based protein profiling (ABPP) has emerged as a powerful functional chemoproteomic strategy which enables global profiling of proteome reactivity toward bioactive small molecules in complex biological and/or pathological processes. To quantify the degree of reactivity in a site-specific manner, an isotopic tandem orthogonal proteolysis (isoTOP)-ABPP strategy has been developed; however, the high cost and long workflow associated with the synthesis of isotopically labeled cleavable tags limit its wide use. Herein, we combined reductive dimethyl labeling with TOP-ABPP to develop a fast, affordable, and efficient method, termed "rdTOP-ABPP", for quantitative chemical proteomics with site-specific precision and triplex quantification. The rdTOP-ABPP method shows high accuracy and precision, good reproducibility, and better capacity for site identification and quantification and is highly compatible with many commercially available cleavable tags. We demonstrated the power of rdTOP-ABPP by profiling the target of (1 S,3 R)-RSL3, a canonical inducer for cell ferroptosis, and provided the first global portrait of its proteome reactivity in a quantitative and site-specific manner.

Identification of 113 new histone marks by CHiMA, a tailored database search strategy.

Shotgun proteomics has been widely used to identify histone marks. Conventional database search methods rely on the "target-decoy" strategy to calculate the false discovery rate (FDR) and distinguish true peptide-spectrum matches (PSMs) from false ones. This strategy has a caveat of inaccurate FDR caused by the small data size of histone marks. To address this challenge, we developed a tailored database search strategy, named "Comprehensive Histone Mark Analysis (CHiMA)." Instead of target-decoy-based FDR, this method uses "50% matched fragment ions" as the key criterion to identify high-confidence PSMs. CHiMA identified twice as many histone modification sites as the conventional method in benchmark datasets. Reanalysis of our previous proteomics data using CHiMA led to the identification of 113 new histone marks for four types of lysine acylations, almost doubling the number of previously reported marks. This tool not only offers a valuable approach for identifying histone modifications but also greatly expands the repertoire of histone marks.

An accelerated and optimized algorithm of selenium-encoded isotopic signature targeted profiling for global selenoproteome analysis.

Selenoproteins play crucial roles including protection and recovery from oxidative stress in organisms. Direct profiling of selenoproteins in proteomes is challenging due to their extremely low abundance. We have developed a computational algorithm termed selenium-encoded isotopic signature targeted profiling (SESTAR) to increase the sensitivity of detecting selenoproteins in complex proteomic samples. In this chapter, we briefly described the basic algorithm of SESTAR. We then introduced SESTAR++, an updated version of SESTAR, with accelerated computation speed and lowered false positive rate. We also provided a detailed workflow to apply SESTAR++ to proteomic profiling of selenoproteins, including the instruction of running the software and implementing it in a targeted profiling mode.

Class I histone deacetylases (HDAC1-3) are histone lysine delactylases.

Lysine L-lactylation [K(L-la)] is a newly discovered histone mark stimulated under conditions of high glycolysis, such as the Warburg effect. K(L-la) is associated with functions that are different from the widely studied histone acetylation. While K(L-la) can be introduced by the acetyltransferase p300, histone delactylases enzymes remained unknown. Here, we report the systematic evaluation of zinc- and nicotinamide adenine dinucleotide­dependent histone deacetylases (HDACs) for their ability to cleave ε-N-L-lactyllysine marks. Our screens identified HDAC1­3 and SIRT1­3 as delactylases in vitro. HDAC1­3 show robust activity toward not only K(L-la) but also K(D-la) and diverse short-chain acyl modifications. We further confirmed the de-L-lactylase activity of HDACs 1 and 3 in cells. Together, these data suggest that histone lactylation is installed and removed by regulatory enzymes as opposed to spontaneous chemical reactivity. Our results therefore represent an important step toward full characterization of this pathway's regulatory elements.

Structural insights into p300 regulation and acetylation-dependent genome organisation.

Histone modifications are deposited by chromatin modifying enzymes and read out by proteins that recognize the modified state. BRD4-NUT is an oncogenic fusion protein of the acetyl lysine reader BRD4 that binds to the acetylase p300 and enables formation of long-range intra- and interchromosomal interactions. We here examine how acetylation reading and writing enable formation of such interactions. We show that NUT contains an acidic transcriptional activation domain that binds to the TAZ2 domain of p300. We use NMR to investigate the structure of the complex and found that the TAZ2 domain has an autoinhibitory role for p300. NUT-TAZ2 interaction or mutations found in cancer that interfere with autoinhibition by TAZ2 allosterically activate p300. p300 activation results in a self-organizing, acetylation-dependent feed-forward reaction that enables long-range interactions by bromodomain multivalent acetyl-lysine binding. We discuss the implications for chromatin organisation, gene regulation and dysregulation in disease.

Histone lysine methacrylation is a dynamic post-translational modification regulated by HAT1 and SIRT2.

Histone lysine crotonylation is a posttranslational modification with demonstrated functions in transcriptional regulation. Here we report the discovery of a new type of histone posttranslational modification, lysine methacrylation (Kmea), corresponding to a structural isomer of crotonyllysine. We validate the identity of this modification using diverse chemical approaches and further confirm the occurrence of this type of histone mark by pan specific and site-specific anti-methacryllysine antibodies. In total, we identify 27 Kmea modified histone sites in HeLa cells using affinity enrichment with a pan Kmea antibody and mass spectrometry. Subsequent biochemical studies show that histone Kmea is a dynamic mark, which is controlled by HAT1 as a methacryltransferase and SIRT2 as a de-methacrylase. Altogether, these investigations uncover a new type of enzyme-catalyzed histone modification and suggest that methacrylyl-CoA generating metabolism is part of a growing number of epigenome-associated metabolic pathways.

A Quantitative Chemoproteomic Platform to Monitor Selenocysteine Reactivity within a Complex Proteome.

Mammalian selenocysteine (Sec)-containing proteins, selenoproteins, are important to (patho)physiological processes, including redox homeostasis. Sec residues have been recalcitrant to mass spectrometry-based chemoproteomic methods that enrich for reactive cysteine (Cys) residues with electrophilic chemical probes, despite confirmed reactivity of Sec with these electrophiles. Highly abundant Cys peptides likely suppress low-abundant Sec peptides. By exploiting the decreased pKa of Sec relative to Cys, we have developed a chemoproteomic platform that relies on low pH (pH 5.75) electrophile labeling, reducing Cys reactivity and enhancing identification of Sec-containing peptides across mouse tissues and cell lines. The utility of this Sec-profiling platform is underscored by evaluation of the selectivity of auranofin, an inhibitor of the selenoprotein, thioredoxin reductase, against both reactive Cys- and Sec-containing proteins. Platform limitations pertain to the non-physiological low-pH conditions that could perturb protein structure and function. Future work necessitates the discovery of Sec-selective electrophiles that function at physiological pH.

Quantitative Profiling of Protein O-GlcNAcylation Sites by an Isotope-Tagged Cleavable Linker.

Large-scale quantification of protein O-linked ß- N-acetylglucosamine (O-GlcNAc) modification in a site-specific manner remains a key challenge in studying O-GlcNAc biology. Herein, we developed an isotope-tagged cleavable linker (isoTCL) strategy, which enabled isotopic labeling of O-GlcNAc through bioorthogonal conjugation of affinity tags. We demonstrated the application of the isoTCL in mapping and quantification of O-GlcNAcylation sites in HeLa cells. Furthermore, we investigated the O-GlcNAcylation sensitivity to the sugar donor by quantifying the levels of modification under different concentrations of the O-GlcNAc labeling probe in a site-specific manner. In addition, we applied isoTCL to compare the O-GlcNAcylation stoichiometry levels of more than 100 modification sites between placenta samples from male and female mice and confirmed site-specifically that female placenta has a higher O-GlcNAcylation than its male counterpart. The isoTCL platform provides a powerful tool for quantitative profiling of O-GlcNAc modification.

Selenium-Encoded Isotopic Signature Targeted Profiling.

Selenium (Se), as an essential trace element, plays crucial roles in many organisms including humans. The biological functions of selenium are mainly mediated by selenoproteins, a unique class of selenium-containing proteins in which selenium is inserted in the form of selenocysteine. Due to their low abundance and uneven tissue distribution, detection of selenoproteins within proteomes is very challenging, and therefore functional studies of these proteins are limited. In this study, we developed a computational method, named as selenium-encoded isotopic signature targeted profiling (SESTAR), which utilizes the distinct natural isotopic distribution of selenium to assist detection of trace selenium-containing signals from shotgun-proteomic data. SESTAR can detect femtomole quantities of synthetic selenopeptides in a benchmark test and dramatically improved detection of native selenoproteins from tissue proteomes in a targeted profiling mode. By applying SESTAR to screen publicly available datasets from Human Proteome Map, we provide a comprehensive picture of selenoprotein distributions in human primary hematopoietic cells and tissues. We further demonstrated that SESTAR can aid chemical-proteomic strategies to identify additional selenoprotein targets of RSL3, a canonical inducer of cell ferroptosis. We believe SESTAR not only serves as a powerful tool for global profiling of native selenoproteomes, but can also work seamlessly with chemical-proteomic profiling strategies to enhance identification of target proteins, post-translational modifications, or protein-protein interactions.

RNA-seq analysis of transcriptome and glucosinolate metabolism in seeds and sprouts of broccoli (Brassica oleracea var. italic).

BACKGROUND: Broccoli (Brassica oleracea var. italica), a member of Cruciferae, is an important vegetable containing high concentration of various nutritive and functional molecules especially the anticarcinogenic glucosinolates. The sprouts of broccoli contain 10-100 times higher level of glucoraphanin, the main contributor of the anticarcinogenesis, than the edible florets. Despite the broccoli sprouts' functional importance, currently available genetic and genomic tools for their studies are very limited, which greatly restricts the development of this functionally important vegetable. RESULTS: A total of ∼85 million 251 bp reads were obtained. After de novo assembly and searching the assembled transcripts against the Arabidopsis thaliana and NCBI nr databases, 19,441 top-hit transcripts were clustered as unigenes with an average length of 2,133 bp. These unigenes were classified according to their putative functional categories. Cluster analysis of total unigenes with similar expression patterns and differentially expressed unigenes among different tissues, as well as transcription factor analysis were performed. We identified 25 putative glucosinolate metabolism genes sharing 62.04-89.72% nucleotide sequence identity with the Arabidopsis orthologs. This established a broccoli glucosinolate metabolic pathway with high colinearity to Arabidopsis. Many of the biosynthetic and degradation genes showed higher expression after germination than in seeds; especially the expression of the myrosinase TGG2 was 20-130 times higher. These results along with the previous reports about these genes' studies in Arabidopsis and the glucosinolate concentration in broccoli sprouts indicate the breakdown products of glucosinolates may play important roles in the stage of broccoli seed germination and sprout development. CONCLUSION: Our study provides the largest genetic resource of broccoli to date. These data will pave the way for further studies and genetic engineering of broccoli sprouts and will also provide new insight into the genomic research of this species and its relatives.