Ketolysis drives CD8+ T cell effector function through effects on histone acetylation.

Environmental nutrient availability influences T cell metabolism, impacting T cell function and shaping immune outcomes. Here, we identified ketone bodies (KBs)-including ß-hydroxybutyrate (ßOHB) and acetoacetate (AcAc)-as essential fuels supporting CD8+ T cell metabolism and effector function. ßOHB directly increased CD8+ T effector (Teff) cell cytokine production and cytolytic activity, and KB oxidation (ketolysis) was required for Teff cell responses to bacterial infection and tumor challenge. CD8+ Teff cells preferentially used KBs over glucose to fuel the tricarboxylic acid (TCA) cycle in vitro and in vivo. KBs directly boosted the respiratory capacity and TCA cycle-dependent metabolic pathways that fuel CD8+ T cell function. Mechanistically, ßOHB was a major substrate for acetyl-CoA production in CD8+ T cells and regulated effector responses through effects on histone acetylation. Together, our results identify cell-intrinsic ketolysis as a metabolic and epigenetic driver of optimal CD8+ T cell effector responses.

Examining the Roles of H3K4 Methylation States with Systematically Characterized Antibodies.

Histone post-translational modifications (PTMs) are important genomic regulators often studied by chromatin immunoprecipitation (ChIP), whereby their locations and relative abundance are inferred by antibody capture of nucleosomes and associated DNA. However, the specificity of antibodies within these experiments has not been systematically studied. Here, we use histone peptide arrays and internally calibrated ChIP (ICeChIP) to characterize 52 commercial antibodies purported to distinguish the H3K4 methylforms (me1, me2, and me3, with each ascribed distinct biological functions). We find that many widely used antibodies poorly distinguish the methylforms and that high- and low-specificity reagents can yield dramatically different biological interpretations, resulting in substantial divergence from the literature for numerous H3K4 methylform paradigms. Using ICeChIP, we also discern quantitative relationships between enhancer H3K4 methylation and promoter transcriptional output and can measure global PTM abundance changes. Our results illustrate how poor antibody specificity contributes to the "reproducibility crisis," demonstrating the need for rigorous, platform-appropriate validation.

DNA strand asymmetry generated by CpG hemimethylation has opposing effects on CTCF binding.

CpG methylation generally occurs on both DNA strands and is essential for mammalian development and differentiation. Until recently, hemimethylation, in which only one strand is methylated, was considered to be simply a transitory state generated during DNA synthesis. The discovery that a subset of CCCTC-binding factor (CTCF) binding sites is heritably hemimethylated suggests that hemimethylation might have an unknown biological function. Here we show that the binding of CTCF is profoundly altered by which DNA strand is methylated and by the specific CTCF binding motif. CpG methylation on the motif strand can inhibit CTCF binding by up to 7-fold, whereas methylation on the opposite strand can stimulate binding by up to 4-fold. Thus, hemimethylation can alter binding by up to 28-fold in a strand-specific manner. The mechanism for sensing methylation on the opposite strand requires two critical residues, V454 and S364, within CTCF zinc fingers 7 and 4. Similar to methylation, CpG hydroxymethylation on the motif strand can inhibit CTCF binding by up to 4-fold. However, hydroxymethylation on the opposite strand removes the stimulatory effect. Strand-specific methylation states may therefore provide a mechanism to explain the transient and dynamic nature of CTCF-mediated chromatin interactions.

An Interactive Database for the Assessment of Histone Antibody Specificity.

Access to high-quality antibodies is a necessity for the study of histones and their posttranslational modifications (PTMs). Here we debut the Histone Antibody Specificity Database (http://www.histoneantibodies.com), an online and expanding resource cataloging the behavior of widely used, commercially available histone antibodies by peptide microarray. This interactive web portal provides a critical resource to the biological research community that routinely uses these antibodies as detection reagents for a wide range of applications.

A physical basis for quantitative ChIP-sequencing.

ChIP followed by next-generation sequencing (ChIP-Seq) is a key technique for mapping the distribution of histone posttranslational modifications (PTMs) and chromatin-associated factors across genomes. There is a perceived challenge to define a quantitative scale for ChIP-Seq data, and as such, several approaches making use of exogenous additives, or "spike-ins," have recently been developed. Herein, we report on the development of a quantitative, physical model defining ChIP-Seq. The quantitative scale on which ChIP-Seq results should be compared emerges from the model. To test the model and demonstrate the quantitative scale, we examine the impacts of an EZH2 inhibitor through the lens of ChIP-Seq. We report a significant increase in immunoprecipitation of presumed off-target histone PTMs after inhibitor treatment, a trend predicted by the model but contrary to spike-in-based indications. Our work also identifies a sensitivity issue in spike-in normalization that has not been considered in the literature, placing limitations on its utility and trustworthiness. We call our new approach the sans-spike-in method for quantitative ChIP-sequencing (siQ-ChIP). A number of changes in community practice of ChIP-Seq, data reporting, and analysis are motivated by this work.

Molecular mechanisms of fentanyl mediated ß-arrestin biased signaling.

The development of novel analgesics with improved safety profiles to combat the opioid epidemic represents a central question to G protein coupled receptor structural biology and pharmacology: What chemical features dictate G protein or ß-arrestin signaling? Here we use adaptively biased molecular dynamics simulations to determine how fentanyl, a potent ß-arrestin biased agonist, binds the µ-opioid receptor (µOR). The resulting fentanyl-bound pose provides rational insight into a wealth of historical structure-activity-relationship on its chemical scaffold. Following an in-silico derived hypothesis we found that fentanyl and the synthetic opioid peptide DAMGO require M153 to induce ß-arrestin coupling, while M153 was dispensable for G protein coupling. We propose and validate an activation mechanism where the n-aniline ring of fentanyl mediates µOR ß-arrestin through a novel M153 "microswitch" by synthesizing fentanyl-based derivatives that exhibit complete, clinically desirable, G protein biased coupling. Together, these results provide molecular insight into fentanyl mediated ß-arrestin biased signaling and a rational framework for further optimization of fentanyl-based analgesics with improved safety profiles.

Chromatin structure and its chemical modifications regulate the ubiquitin ligase substrate selectivity of UHRF1.

Mitotic inheritance of DNA methylation patterns is facilitated by UHRF1, a DNA- and histone-binding E3 ubiquitin ligase that helps recruit the maintenance DNA methyltransferase DNMT1 to replicating chromatin. The DNA methylation maintenance function of UHRF1 is dependent on its ability to bind chromatin, where it facilitates monoubiquitination of histone H3 at lysines 18 and 23, a docking site for DNMT1. Because of technical limitations, this model of UHRF1-dependent DNA methylation inheritance has been constructed largely based on genetics and biochemical observations querying methylated DNA oligonucleotides, synthetic histone peptides, and heterogeneous chromatin extracted from cells. Here, we construct semisynthetic mononucleosomes harboring defined histone and DNA modifications and perform rigorous analysis of UHRF1 binding and enzymatic activity with these reagents. We show that multivalent engagement of nucleosomal linker DNA and dimethylated lysine 9 on histone H3 directs UHRF1 ubiquitin ligase activity toward histone substrates. Notably, we reveal a molecular switch, stimulated by recognition of hemimethylated DNA, which redirects UHRF1 ubiquitin ligase activity away from histones in favor of robust autoubiquitination. Our studies support a noncompetitive model for UHRF1 and DNMT1 chromatin recruitment to replicating chromatin and define a role for hemimethylated linker DNA as a regulator of UHRF1 ubiquitin ligase substrate selectivity.

Multivalent histone engagement by the linked tandem Tudor and PHD domains of UHRF1 is required for the epigenetic inheritance of DNA methylation.

Histone post-translational modifications regulate chromatin structure and function largely through interactions with effector proteins that often contain multiple histone-binding domains. While significant progress has been made characterizing individual effector domains, the role of paired domains and how they function in a combinatorial fashion within chromatin are poorly defined. Here we show that the linked tandem Tudor and plant homeodomain (PHD) of UHRF1 (ubiquitin-like PHD and RING finger domain-containing protein 1) operates as a functional unit in cells, providing a defined combinatorial readout of a heterochromatin signature within a single histone H3 tail that is essential for UHRF1-directed epigenetic inheritance of DNA methylation. These findings provide critical support for the "histone code" hypothesis, demonstrating that multivalent histone engagement plays a key role in driving a fundamental downstream biological event in chromatin.

The finger loop of the SRA domain in the E3 ligase UHRF1 is a regulator of ubiquitin targeting and is required for the maintenance of DNA methylation.

The Su(var)3-9, enhancer of zeste, and trithorax (SET) and really interesting new gene (RING) finger-associated (SRA) protein domain is conserved across bacteria and eukaryota and coordinates extrahelical or "flipped" DNA bases. A functional SRA domain is required for ubiquitin-like with PHD and RING finger domains 1 (UHRF1) E3 ubiquitin ligase activity toward histone H3, a mechanism for recruiting the DNA methylation maintenance enzyme DNA methyltransferase 1 (DNMT1). The SRA domain supports UHRF1 oncogenic activity in colon cancer cells, highlighting that UHRF1 SRA antagonism could be a cancer therapeutic strategy. Here we used molecular dynamics simulations, DNA binding assays, in vitro ubiquitination reactions, and DNA methylation analysis to identify the SRA finger loop as a regulator of UHRF1 ubiquitin targeting and DNA methylation maintenance. A chimeric UHRF1 (finger swap) with diminished E3 ligase activity toward nucleosomal histones, despite tighter binding to unmodified or asymmetric or symmetrically methylated DNA, uncouples DNA affinity from regulation of E3 ligase activity. Our model suggests that SRA domains sample DNA bases through flipping in the presence or absence of a cytosine modification and that specific interactions of the SRA finger loop with DNA are required for downstream host protein function. Our findings provide insight into allosteric regulation of UHRF1 E3 ligase activity, suggesting that UHRF1's SRA finger loop regulates its conformation and function.

How fluxional reactants limit the accuracy/efficiency of infrequent metadynamics.

In an infrequent metadynamics (iMetaD) simulation, a well-tempered metadynamics bias accumulates in the reactant basin, accelerating escapes to the product state. Like the earlier hyperdynamics strategy, iMetaD enables estimates of the unbiased escape rates. However, iMetaD applies the bias to visited locations in a collective variable (CV) space, not to the more specific visited locations in a full configuration space as done in hyperdynamics. This difference makes rate estimates from iMetaD sensitive to the choice of CVs, to parameters that control the bias deposition rate, and to the preparation of the initial state within the reactant basin. This paper uses an extremely simple discrete state model to illustrate complications that can arise in systems that exhibit fluxional transitions between sub-basins of the reactant state. Specifically, we show how the reactant-to-product escape time and relaxation times within the reactant basin(s) impose bounds on the admissible parameter choices for an iMetaD calculation. Predictions from the discrete state model are validated by iMetaD simulations on a corresponding two-dimensional potential energy surface.

Comparative biochemical analysis of UHRF proteins reveals molecular mechanisms that uncouple UHRF2 from DNA methylation maintenance.

UHRF1 is a histone- and DNA-binding E3 ubiquitin ligase that functions with DNMT1 to maintain mammalian DNA methylation. UHRF1 facilitates DNMT1 recruitment to replicating chromatin through a coordinated mechanism involving histone and DNA recognition and histone ubiquitination. UHRF2 shares structural homology with UHRF1, but surprisingly lacks functional redundancy to facilitate DNA methylation maintenance. Molecular mechanisms uncoupling UHRF2 from DNA methylation maintenance are poorly defined. Through comprehensive and comparative biochemical analysis of recombinant human UHRF1 and UHRF2 reader and writer activities, we reveal conserved modes of histone PTM recognition but divergent DNA binding properties. While UHRF1 and UHRF2 diverge in their affinities toward hemi-methylated DNA, we surprisingly show that both hemi-methylated and hemi-hydroxymethylated DNA oligonucleotides stimulate UHRF2 ubiquitin ligase activity toward histone H3 peptide substrates. This is the first example of an E3 ligase allosterically regulated by DNA hydroxymethylation. However, UHRF2 is not a productive histone E3 ligase toward purified mononucleosomes, suggesting UHRF2 has an intra-domain architecture distinct from UHRF1 that is conformationally constrained when bound to chromatin. Collectively, our studies reveal that uncoupling of UHRF2 from the DNA methylation maintenance program is linked to differences in the molecular readout of chromatin signatures that connect UHRF1 to ubiquitination of histone H3.

Quantitative Characterization of Bivalent Probes for a Dual Bromodomain Protein, Transcription Initiation Factor TFIID Subunit 1.

Multivalent binding is an efficient means to enhance the affinity and specificity of chemical probes targeting multidomain proteins in order to study their function and role in disease. While the theory of multivalent binding is straightforward, physical and structural characterization of bivalent binding encounters multiple technical difficulties. We present a case study where a combination of experimental techniques and computational simulations was used to comprehensively characterize the binding and structure-affinity relationships for a series of Bromosporine-based bivalent bromodomain ligands with a bivalent protein, Transcription Initiation Factor TFIID subunit 1 (TAF1). Experimental techniques-Isothermal Titration Calorimetry, X-ray Crystallography, Circular Dichroism, Size Exclusion Chromatography-Multi-Angle Light Scattering, and Surface Plasmon Resonance-were used to determine structures, binding affinities, and kinetics of monovalent ligands and bivalent ligands with varying linker lengths. The experimental data for monomeric ligands were fed into explicit computational simulations, in which both ligand and protein species were present in a broad range of concentrations, and in up to a 100 s time regime, to match experimental conditions. These simulations provided accurate estimates for apparent affinities (in good agreement with experimental data), individual dissociation microconstants and other microscopic details for each type of protein-ligand complex. We conclude that the expected efficiency of bivalent ligands in a cellular context is difficult to estimate by a single technique in vitro, due to higher order associations favored at the concentrations used, and other complicating processes. Rather, a combination of structural, biophysical, and computational approaches should be utilized to estimate and characterize multivalent interactions.

A cellular chemical probe targeting the chromodomains of Polycomb repressive complex 1.

We report the design and characterization of UNC3866, a potent antagonist of the methyllysine (Kme) reading function of the Polycomb CBX and CDY families of chromodomains. Polycomb CBX proteins regulate gene expression by targeting Polycomb repressive complex 1 (PRC1) to sites of H3K27me3 via their chromodomains. UNC3866 binds the chromodomains of CBX4 and CBX7 most potently, with a K(d) of ∼100 nM for each, and is 6- to 18-fold selective as compared to seven other CBX and CDY chromodomains while being highly selective over >250 other protein targets. X-ray crystallography revealed that UNC3866's interactions with the CBX chromodomains closely mimic those of the methylated H3 tail. UNC4195, a biotinylated derivative of UNC3866, was used to demonstrate that UNC3866 engages intact PRC1 and that EED incorporation into PRC1 is isoform dependent in PC3 prostate cancer cells. Finally, UNC3866 inhibits PC3 cell proliferation, consistent with the known ability of CBX7 overexpression to confer a growth advantage, whereas UNC4219, a methylated negative control compound, has negligible effects.

The structural basis of the dominant negative phenotype of the Gαi1ß1γ2 G203A/A326S heterotrimer.

AIM: Dominant negative mutant G proteins have provided critical insight into the mechanisms of G protein-coupled receptor (GPCR) signaling, but the mechanisms underlying the dominant negative characteristics are not completely understood. The aim of this study was to determine the structure of the dominant negative Gαi1ß1γ2 G203A/A326S complex (Gi-DN) and to reveal the structural basis of the mutation-induced phenotype of Gαi1ß1γ2. METHODS: The three subunits of the Gi-DN complex were co-expressed with a baculovirus expression system. The Gi-DN heterotrimer was purified, and the structure of its complex with GDP was determined through X-ray crystallography. RESULTS: The Gi-DN heterotrimer structure revealed a dual mechanism underlying the dominant negative characteristics. The mutations weakened the hydrogen bonding network between GDP/GTP and the binding pocket residues, and increased the interactions in the Gα-Gßγ interface. Concomitantly, the Gi-DN heterotrimer adopted a conformation, in which the C-terminus of Gαi and the N-termini of both the Gß and Gγ subunits were more similar to the GPCR-bound state compared with the wild type complex. From these structural observations, two additional mutations (T48F and D272F) were designed that completely abolish the GDP binding of the Gi-DN heterotrimer. CONCLUSION: Overall, the results suggest that the mutations impede guanine nucleotide binding and Gα-Gßγ protein dissociation and favor the formation of the G protein/GPCR complex, thus blocking signal propagation. In addition, the structure provides a rationale for the design of other mutations that cause dominant negative effects in the G protein, as exemplified by the T48F and D272F mutations.

A fast, open source implementation of adaptive biasing potentials uncovers a ligand design strategy for the chromatin regulator BRD4.

In this communication we introduce an efficient implementation of adaptive biasing that greatly improves the speed of free energy computation in molecular dynamics simulations. We investigated the use of accelerated simulations to inform on compound design using a recently reported and clinically relevant inhibitor of the chromatin regulator BRD4 (bromodomain-containing protein 4). Benchmarking on our local compute cluster, our implementation achieves up to 2.5 times more force calls per day than plumed2. Results of five 1 µs-long simulations are presented, which reveal a conformational switch in the BRD4 inhibitor between a binding competent and incompetent state. Stabilization of the switch led to a -3 kcal/mol improvement of absolute binding free energy. These studies suggest an unexplored ligand design principle and offer new actionable hypotheses for medicinal chemistry efforts against this druggable epigenetic target class.

µ-tempered metadynamics: Artifact independent convergence times for wide hills.

Recent analysis of well-tempered metadynamics (WTmetaD) showed that it converges without mollification artifacts in the bias potential. Here, we explore how metadynamics heals mollification artifacts, how healing impacts convergence time, and whether alternative temperings may be used to improve efficiency. We introduce "µ-tempered" metadynamics as a simple tempering scheme, inspired by a related mollified adaptive biasing potential, that results in artifact independent convergence of the free energy estimate. We use a toy model to examine the role of artifacts in WTmetaD and solvated alanine dipeptide to compare the well-tempered and µ-tempered frameworks demonstrating fast convergence for hill widths as large as 60(∘) for µTmetaD.

UHRF1 ubiquitin ligase activity supports the maintenance of low-density CpG methylation.

The RING E3 ubiquitin ligase UHRF1 is an established cofactor for DNA methylation inheritance. Nucleosomal engagement through histone and DNA interactions directs UHRF1 ubiquitin ligase activity toward lysines on histone H3 tails, creating binding sites for DNMT1 through ubiquitin interacting motifs (UIM1 and UIM2). Here, we profile contributions of UHRF1 and DNMT1 to genome-wide DNA methylation inheritance and dissect specific roles for ubiquitin signaling in this process. We reveal DNA methylation maintenance at low-density CpGs is vulnerable to disruption of UHRF1 ubiquitin ligase activity and DNMT1 ubiquitin reading activity through UIM1. Hypomethylation of low-density CpGs in this manner induces formation of partially methylated domains (PMD), a methylation signature observed across human cancers. Furthermore, disrupting DNMT1 UIM2 function abolishes DNA methylation maintenance. Collectively, we show DNMT1-dependent DNA methylation inheritance is a ubiquitin-regulated process and suggest a disrupted UHRF1-DNMT1 ubiquitin signaling axis contributes to the development of PMDs in human cancers.

Molecular dynamics simulations provide insights into ULK-101 potency and selectivity toward autophagic kinases ULK1/2.

Kinase domains are highly conserved within protein kinases in both sequence and structure. Many factors, including phosphorylation, amino acid substitutions or mutations, and small molecule inhibitor binding, influence conformations of the kinase domain and enzymatic activity. The serine/threonine kinases ULK1 and ULK2 are highly conserved with N- and C-terminal domains, phosphate-binding P-loops, αC-helix, regulatory and catalytic spines, and activation loop DFG and APE motifs. Here, we performed molecular dynamics (MD) simulations to understand better the potency and selectivity of the ULK1/2 small molecule inhibitor, ULK-101. We observed stable bound states for ULK-101 to the adenosine triphosphate (ATP)-binding site of ULK2, coordinated by hydrogen bonding with the hinge backbone and the catalytic lysine sidechain. Notably, ULK-101 occupies a hydrophobic pocket associated with the N-terminus of the αC-helix. Large movements in the P-loop are also associated with ULK-101 inhibitor binding and exit from ULK2. Our data further suggests that ULK-101 could induce a folded P-loop conformation and hydrophobic pocket reflected in its nanomolar potency and kinome selectivity.