Acetyl-methyllysine marks chromatin at active transcription start sites.

Lysine residues in histones and other proteins can be modified by post-translational modifications that encode regulatory information1. Lysine acetylation and methylation are especially important for regulating chromatin and gene expression2-4. Pathways involving these post-translational modifications are targets for clinically approved therapeutics to treat human diseases. Lysine methylation and acetylation are generally assumed to be mutually exclusive at the same residue. Here we report cellular lysine residues that are both methylated and acetylated on the same side chain to form Nε-acetyl-Nε-methyllysine (Kacme). We show that Kacme is found on histone H4 (H4Kacme) across a range of species and across mammalian tissues. Kacme is associated with marks of active chromatin, increased transcriptional initiation and is regulated in response to biological signals. H4Kacme can be installed by enzymatic acetylation of monomethyllysine peptides and is resistant to deacetylation by some HDACs in vitro. Kacme can be bound by chromatin proteins that recognize modified lysine residues, as we demonstrate with the crystal structure of acetyllysine-binding protein BRD2 bound to a histone H4Kacme peptide. These results establish Kacme as a cellular post-translational modification with the potential to encode information distinct from methylation and acetylation alone and demonstrate that Kacme has all the hallmarks of a post-translational modification with fundamental importance to chromatin biology.

Improving the study of RNA dynamics through advances in RNA-seq with metabolic labeling and nucleotide-recoding chemistry.

RNA metabolic labeling using 4-thiouridine (s4U) captures the dynamics of RNA synthesis and decay. The power of this approach is dependent on appropriate quantification of labeled and unlabeled sequencing reads, which can be compromised by the apparent loss of s4U-labeled reads in a process we refer to as dropout. Here we show that s4U-containing transcripts can be selectively lost when RNA samples are handled under sub-optimal conditions, but that this loss can be minimized using an optimized protocol. We demonstrate a second cause of dropout in nucleotide recoding and RNA sequencing (NR-seq) experiments that is computational and downstream of library preparation. NR-seq experiments involve chemically converting s4U from a uridine analog to a cytidine analog and using the apparent T-to-C mutations to identify the populations of newly synthesized RNA. We show that high levels of T-to-C mutations can prevent read alignment with some computational pipelines, but that this bias can be overcome using improved alignment pipelines. Importantly, kinetic parameter estimates are affected by dropout independent of the NR chemistry employed, and all chemistries are practically indistinguishable in bulk, short-read RNA-seq experiments. Dropout is an avoidable problem that can be identified by including unlabeled controls, and mitigated through improved sample handing and read alignment that together improve the robustness and reproducibility of NR-seq experiments.

bakR: uncovering differential RNA synthesis and degradation kinetics transcriptome-wide with Bayesian hierarchical modeling.

Differential expression analysis of RNA sequencing (RNA-seq) data can identify changes in cellular RNA levels, but provides limited information about the kinetic mechanisms underlying such changes. Nucleotide recoding RNA-seq methods (NR-seq; e.g., TimeLapse-seq, SLAM-seq, etc.) address this shortcoming and are widely used approaches to identify changes in RNA synthesis and degradation kinetics. While advanced statistical models implemented in user-friendly software (e.g., DESeq2) have ensured the statistical rigor of differential expression analyses, no such tools that facilitate differential kinetic analysis with NR-seq exist. Here, we report the development of Bayesian analysis of the kinetics of RNA (bakR; https:// github.com/simonlabcode/bakR), an R package to address this need. bakR relies on Bayesian hierarchical modeling of NR-seq data to increase statistical power by sharing information across transcripts. Analyses of simulated data confirmed that bakR implementations of the hierarchical model outperform attempts to analyze differential kinetics with existing models. bakR also uncovers biological signals in real NR-seq data sets and provides improved analyses of existing data sets. This work establishes bakR as an important tool for identifying differential RNA synthesis and degradation kinetics.

Internally controlled RNA sequencing comparisons using nucleoside recoding chemistry.

Quantitative comparisons of RNA levels from different samples can lead to new biological understanding if they are able to distinguish biological variation from variable sample preparation. These challenges are pronounced in comparisons that require complex biochemical manipulations (e.g. isolating polysomes to study translation). Here, we present Transcript Regulation Identified by Labeling with Nucleoside Analogues in Cell Culture (TILAC), an internally controlled approach for quantitative comparisons of RNA content. TILAC uses two metabolic labels, 4-thiouridine (s4U) and 6-thioguanosine (s6G), to differentially label RNAs in cells, allowing experimental and control samples to be pooled prior to downstream biochemical manipulations. TILAC leverages nucleoside recoding chemistry to generate characteristic sequencing signatures for each label and uses statistical modeling to compare the abundance of RNA transcripts between samples. We verified the performance of TILAC in transcriptome-scale experiments involving RNA polymerase II inhibition and heat shock. We then applied TILAC to quantify changes in mRNA association with actively translating ribosomes during sodium arsenite stress and discovered a set of transcripts that are translationally upregulated, including MCM2 and DDX5. TILAC is broadly applicable to uncover differences between samples leading to improved biological insights.

A Dual-Site Inhibitor of CBP/p300 KIX is a Selective and Effective Modulator of Myb.

The protein-protein interaction between the KIX motif of the transcriptional coactivator CBP/p300 and the transcriptional activator Myb is a high-value target due to its established role in certain acute myeloid leukemias (AML) and potential contributions to other cancers. However, the CBP/p300 KIX domain has multiple binding sites, several structural homologues, many binding partners, and substantial conformational plasticity, making it challenging to specifically target using small-molecule inhibitors. Here, we report a picomolar dual-site inhibitor (MybLL-tide) of the Myb-CBP/p300 KIX interaction. MybLL-tide has higher affinity for CBP/p300 KIX than any previously reported compounds while also possessing 5600-fold selectivity for the CBP/p300 KIX domain over other coactivator domains. MybLL-tide blocks the association of CBP and p300 with Myb in the context of the proteome, leading to inhibition of key Myb·KIX-dependent genes in AML cells. These results show that MybLL-tide is an effective, modifiable tool to selectively target the KIX domain and assess transcriptional effects in AML cells and potentially other cancers featuring aberrant Myb behavior. Additionally, the dual-site design has applicability to the other challenging coactivators that bear multiple binding surfaces.