Structural basis for the inhibition of PRC2 by active transcription histone posttranslational modifications.

Polycomb repressive complex 2 (PRC2) is an epigenetic regulator essential for embryonic development and maintenance of cell identity that trimethylates histone H3 at lysine 27 (H3K27me3) leading to gene silencing. PRC2 is regulated by association with protein cofactors and crosstalk with histone posttranslational modifications. Trimethylated histone H3 K4 (H3K4me3) and K36 (H3K36me3) localize to sites of active transcription where H3K27me3 is absent and inhibit PRC2 activity through unknown mechanisms. Using cryo-electron microscopy we reveal that histone H3 tails modified with H3K36me3 engage poorly with the PRC2 active site and preclude its effective interaction with chromatin, while the H3K4me3 modification binds to the allosteric site in the EED subunit, acting as an antagonist that competes with allosteric activators required for the spreading of the H3K27me3 repressive mark. Thus, the location along the H3 tail of the H3K4me3 and H3K36me3 modifications allow them to target two essential requirements for efficient trimethylation of histone H3K27. We further show that the JARID2 cofactor modulates PRC2 activity in the presence of these histone modifications.

Structural basis for inactivation of PRC2 by G-quadruplex RNA.

Polycomb repressive complex 2 (PRC2) silences genes through trimethylation of histone H3K27. PRC2 associates with numerous precursor messenger RNAs (pre-mRNAs) and long noncoding RNAs (lncRNAs) with a binding preference for G-quadruplex RNA. In this work, we present a 3.3-Å-resolution cryo-electron microscopy structure of PRC2 bound to a G-quadruplex RNA. Notably, RNA mediates the dimerization of PRC2 by binding both protomers and inducing a protein interface composed of two copies of the catalytic subunit EZH2, thereby blocking nucleosome DNA interaction and histone H3 tail accessibility. Furthermore, an RNA-binding loop of EZH2 facilitates the handoff between RNA and DNA, another activity implicated in PRC2 regulation by RNA. We identified a gain-of-function mutation in this loop that activates PRC2 in zebrafish. Our results reveal mechanisms for RNA-mediated regulation of a chromatin-modifying enzyme.

Structural basis for inactivation of PRC2 by G-quadruplex RNA.

The histone methyltransferase PRC2 (Polycomb Repressive Complex 2) silences genes via successively attaching three methyl groups to lysine 27 of histone H3. PRC2 associates with numerous pre-mRNA and lncRNA transcripts with a binding preference for G-quadruplex RNA. Here, we present a 3.3Å-resolution cryo-EM structure of PRC2 bound to a G-quadruplex RNA. Notably, RNA mediates the dimerization of PRC2 by binding both protomers and inducing a protein interface comprised of two copies of the catalytic subunit EZH2, which limits nucleosome DNA interaction and occludes H3 tail accessibility to the active site. Our results reveal an unexpected mechanism for RNA-mediated inactivation of a chromatin-modifying enzyme. Furthermore, the flexible loop of EZH2 that helps stabilize RNA binding also facilitates the handoff between RNA and DNA, an activity implicated in PRC2 regulation by RNA. One-Sentence Summary: Cryo-EM structure of RNA-bound PRC2 dimer elucidates an unexpected mechanism of PRC2 inhibition by RNA.

DNMT1 inhibition by pUG-fold quadruplex RNA.

Aberrant DNA methylation is one of the earliest hallmarks of cancer. DNMT1 is responsible for methylating newly replicated DNA, but the precise regulation of DNMT1 to ensure faithful DNA methylation remains poorly understood. A link between RNA and chromatin-associated proteins has recently emerged, and several studies have shown that DNMT1 can be regulated by a variety of RNAs. In this study, we have confirmed that human DNMT1 indeed interacts with multiple RNAs, including its own nuclear mRNA. Unexpectedly, we found that DNMT1 exhibits a strong and specific affinity for GU-rich RNAs that form a pUG-fold, a noncanonical G-quadruplex. We find that pUG-fold-capable RNAs inhibit DNMT1 activity by inhibiting binding of hemimethylated DNA, and we additionally provide evidence for multiple RNA binding modes with DNMT1. Together, our data indicate that a human chromatin-associated protein binds to and is regulated by pUG-fold RNA.

Structural insights into the interactions of Polycomb Repressive Complex 2 with chromatin.

Polycomb repressive complexes are a family of chromatin modifier enzymes which are critical for regulating gene expression and maintaining cell-type identity. The reversible chemical modifications of histone H3 and H2A by the Polycomb proteins are central to its ability to function as a gene silencer. PRC2 is both a reader and writer of the tri-methylation of histone H3 lysine 27 (H3K27me3) which serves as a marker for transcription repression, and heterochromatin boundaries. Over the last few years, several studies have provided key insights into the mechanisms regulating the recruitment and activation of PRC2 at Polycomb target genes. In this review, we highlight the recent structural studies which have elucidated the roles played by Polycomb cofactor proteins in mediating crosstalk between histone post-translational modifications and the recruitment of PRC2 and the stimulation of PRC2 methyltransferase activity.

JARID2 and AEBP2 regulate PRC2 in the presence of H2AK119ub1 and other histone modifications.

Polycomb repressive complexes 1 and 2 (PRC1 and PRC2) cooperate to determine cell identity by epigenetic gene expression regulation. However, the mechanism of PRC2 recruitment by means of recognition of PRC1-mediated H2AK119ub1 remains poorly understood. Our PRC2 cryo-electron microscopy structure with cofactors JARID2 and AEBP2 bound to a H2AK119ub1-containing nucleosome reveals a bridge helix in EZH2 that connects the SET domain, H3 tail, and nucleosomal DNA. JARID2 and AEBP2 each interact with one ubiquitin and the H2A-H2B surface. JARID2 stimulates PRC2 through interactions with both the polycomb protein EED and the H2AK119-ubiquitin, whereas AEBP2 has an additional scaffolding role. The presence of these cofactors partially overcomes the inhibitory effect that H3K4me3 and H3K36me3 exert on core PRC2 (in the absence of cofactors). Our results support a key role for JARID2 and AEBP2 in the cross-talk between histone modifications and PRC2 activity.

Recent Structural Insights into Polycomb Repressive Complex 2 Regulation and Substrate Binding.

Polycomb group proteins are transcriptional repressors controlling gene expression patterns and maintaining cell type identity. The chemical modifications of histones and DNA caused by the regulated activity of chromatin-modifying enzymes such as Polycomb help establish and maintain such expression patterns. Polycomb repressive complex 2 (PRC2) is the only known methyltransferase specific for histone H3 lysine 27 (H3K27) and catalyzes its trimethylation leading to the repressive H3K27me3 mark. Structural biology has made important contributions to our understanding of the molecular mechanisms that ensure the spatiotemporal regulation of PRC2 activity and the establishment of inactive chromatin domains marked by H3K27me3. In this review, we discuss the recent structural studies that have advanced our understanding of PRC2 function, in particular the roles of intersubunit interactions in complex assembly and the regulation of methyltransferase activity, as well as the mechanism of local H3K27me3 spreading leading to repressive domains.

Cryo-EM structures of PRC2 simultaneously engaged with two functionally distinct nucleosomes.

Epigenetic regulation is mediated by protein complexes that couple recognition of chromatin marks to activity or recruitment of chromatin-modifying enzymes. Polycomb repressive complex 2 (PRC2), a gene silencer that methylates lysine 27 of histone H3, is stimulated upon recognition of its own catalytic product and has been shown to be more active on dinucleosomes than H3 tails or single nucleosomes. These properties probably facilitate local H3K27me2/3 spreading, causing heterochromatin formation and gene repression. Here, cryo-EM reconstructions of human PRC2 bound to bifunctional dinucleosomes show how a single PRC2, via interactions with nucleosomal DNA, positions the H3 tails of the activating and substrate nucleosome to interact with the EED subunit and the SET domain of EZH2, respectively. We show how the geometry of the PRC2-DNA interactions allows PRC2 to accommodate varying lengths of the linker DNA between nucleosomes. Our structures illustrate how an epigenetic regulator engages with a complex chromatin substrate.

Structures of human PRC2 with its cofactors AEBP2 and JARID2.

Transcriptionally repressive histone H3 lysine 27 methylation by Polycomb repressive complex 2 (PRC2) is essential for cellular differentiation and development. Here we report cryo-electron microscopy structures of human PRC2 in a basal state and two distinct active states while in complex with its cofactors JARID2 and AEBP2. Both cofactors mimic the binding of histone H3 tails. JARID2, methylated by PRC2, mimics a methylated H3 tail to stimulate PRC2 activity, whereas AEBP2 interacts with the RBAP48 subunit, mimicking an unmodified H3 tail. SUZ12 interacts with all other subunits within the assembly and thus contributes to the stability of the complex. Our analysis defines the complete architecture of a functionally relevant PRC2 and provides a structural framework to understand its regulation by cofactors, histone tails, and RNA.

Entropy in molecular recognition by proteins.

Molecular recognition by proteins is fundamental to molecular biology. Dissection of the thermodynamic energy terms governing protein-ligand interactions has proven difficult, with determination of entropic contributions being particularly elusive. NMR relaxation measurements have suggested that changes in protein conformational entropy can be quantitatively obtained through a dynamical proxy, but the generality of this relationship has not been shown. Twenty-eight protein-ligand complexes are used to show a quantitative relationship between measures of fast side-chain motion and the underlying conformational entropy. We find that the contribution of conformational entropy can range from favorable to unfavorable, which demonstrates the potential of this thermodynamic variable to modulate protein-ligand interactions. For about one-quarter of these complexes, the absence of conformational entropy would render the resulting affinity biologically meaningless. The dynamical proxy for conformational entropy or "entropy meter" also allows for refinement of the contributions of solvent entropy and the loss in rotational-translational entropy accompanying formation of high-affinity complexes. Furthermore, structure-based application of the approach can also provide insight into long-lived specific water-protein interactions that escape the generic treatments of solvent entropy based simply on changes in accessible surface area. These results provide a comprehensive and unified view of the general role of entropy in high-affinity molecular recognition by proteins.

On the relationship between NMR-derived amide order parameters and protein backbone entropy changes.

Molecular dynamics simulations are used to analyze the relationship between NMR-derived squared generalized order parameters of amide NH groups and backbone entropy. Amide order parameters (O(2) NH ) are largely determined by the secondary structure and average values appear unrelated to the overall flexibility of the protein. However, analysis of the more flexible subset (O(2) NH < 0.8) shows that these report both on the local flexibility of the protein and on a different component of the conformational entropy than that reported by the side chain methyl axis order parameters, O(2) axis . A calibration curve for backbone entropy vs. O(2) NH is developed, which accounts for both correlations between amide group motions of different residues, and correlations between backbone and side chain motions. This calibration curve can be used with experimental values of O(2) NH changes obtained by NMR relaxation measurements to extract backbone entropy changes, for example, upon ligand binding. In conjunction with our previous calibration for side chain entropy derived from measured O(2) axis values this provides a prescription for determination of the total protein conformational entropy changes from NMR relaxation measurements.

A sharp thermal transition of fast aromatic-ring dynamics in ubiquitin.

Aromatic amino acid side chains have a rich role within proteins and are often central to their structure and function. Suitable isotopic-labelling strategies enable studies of sub-nanosecond aromatic-ring dynamics using solution NMR relaxation methods. Surprisingly, it was found that the three aromatic side chains in human ubiquitin show a sharp thermal dynamical transition at approximately 312 K. Hydrostatic pressure has little effect on the low-temperature behavior, but somewhat decreases the amplitude of motion in the high-temperature regime. Therefore, below the transition temperature, ring motion is largely librational. Above this temperature, a complete ring-rotation process that is fully consistent with a continuous diffusion not requiring the transient creation of a large activated free volume occurs. Molecular dynamics simulations qualitatively corroborate this view and reinforce the notion that the dynamical character of the protein interior has much more liquid-alkane-like properties than previously appreciated.

Dynamic and thermodynamic response of the Ras protein Cdc42Hs upon association with the effector domain of PAK3.

Human cell division cycle protein 42 (Cdc42Hs) is a small, Rho-type guanosine triphosphatase involved in multiple cellular processes through its interactions with downstream effectors. The binding domain of one such effector, the actin cytoskeleton-regulating p21-activated kinase 3, is known as PBD46. Nitrogen-15 backbone and carbon-13 methyl NMR relaxation was measured to investigate the dynamical changes in activated GMPPCP·Cdc42Hs upon PBD46 binding. Changes in internal motion of the Cdc42Hs, as revealed by methyl axis order parameters, were observed not only near the Cdc42Hs-PBD46 interface but also in remote sites on the Cdc42Hs molecule. The binding-induced changes in side-chain dynamics propagate along the long axis of Cdc42Hs away from the site of PBD46 binding with sharp distance dependence. Overall, the binding of the PBD46 effector domain on the dynamics of methyl-bearing side chains of Cdc42Hs results in a modest rigidification, which is estimated to correspond to an unfavorable change in conformational entropy of approximately -10kcalmol(-1) at 298K. A cluster of methyl probes closest to the nucleotide-binding pocket of Cdc42Hs becomes more rigid upon binding of PBD46 and is proposed to slow the catalytic hydrolysis of the γ phosphate moiety. An additional cluster of methyl probes surrounding the guanine ring becomes more flexible on binding of PBD46, presumably facilitating nucleotide exchange mediated by a guanosine exchange factor. In addition, the Rho insert helix, which is located at a site remote from the PBD46 binding interface, shows a significant dynamic response to PBD46 binding.

Banding of NMR-derived methyl order parameters: implications for protein dynamics.

Our understanding of protein folding, stability, and function has begun to more explicitly incorporate dynamical aspects. Nuclear magnetic resonance has emerged as a powerful experimental method for obtaining comprehensive site-resolved insight into protein motion. It has been observed that methyl-group motion tends to cluster into three "classes" when expressed in terms of the popular Lipari-Szabo model-free squared generalized order parameter. Here the origins of the three classes or bands in the distribution of order parameters are examined. As a first step, a Bayesian based approach, which makes no a priori assumption about the existence or number of bands, is developed to detect the banding of Oaxis2 values derived either from NMR experiments or molecular dynamics simulations. The analysis is applied to seven proteins with extensive molecular dynamics simulations of these proteins in explicit water to examine the relationship between O2 and fine details of the motion of methyl bearing side chains. All of the proteins studied display banding, with some subtle differences. We propose a very simple yet plausible physical mechanism for banding. Finally, our Bayesian method is used to analyze the measured distributions of methyl group motions in the catabolite activating protein and several of its mutants in various liganded states and discuss the functional implications of the observed banding to protein dynamics and function.

Microscopic insights into the NMR relaxation-based protein conformational entropy meter.

Conformational entropy is a potentially important thermodynamic parameter contributing to protein function. Quantitative measures of conformational entropy are necessary for an understanding of its role but have been difficult to obtain. An empirical method that utilizes changes in conformational dynamics as a proxy for changes in conformational entropy has recently been introduced. Here we probe the microscopic origins of the link between conformational dynamics and conformational entropy using molecular dynamics simulations. Simulation of seven proteins gave an excellent correlation with measures of side-chain motion derived from NMR relaxation. The simulations show that the motion of methyl-bearing side chains are sufficiently coupled to that of other side chains to serve as excellent reporters of the overall side-chain conformational entropy. These results tend to validate the use of experimentally accessible measures of methyl motion--the NMR-derived generalized order parameters--as a proxy from which to derive changes in protein conformational entropy.

A 13C labeling strategy reveals a range of aromatic side chain motion in calmodulin.

NMR relaxation experiments often require site-specific isotopic enrichment schemes in order to allow for quantitative interpretation. Here we describe a new labeling scheme for site-specific (13)C-(1)H enrichment of a single ortho position of aromatic amino acid side chains in an otherwise perdeuterated background by employing a combination of [4-(13)C]erythrose and deuterated pyruvate during growth on deuterium oxide. This labeling scheme largely eliminates undesired contributions to (13)C relaxation and greatly simplifies the fitting of relaxation data using the Lipari-Szabo model-free formalism. This approach is illustrated with calcium-saturated vertebrate calmodulin and oxidized flavodoxin from Cyanobacterium anabaena . Analysis of (13)C relaxation in the aromatic groups of calcium-saturated calmodulin indicates a wide range of motion in the subnanosecond time regime.

Coupled motion in proteins revealed by pressure perturbation.

The cooperative nature of protein substructure and internal motion is a critical aspect of their functional competence about which little is known experimentally. NMR relaxation is used here to monitor the effects of high pressure on fast internal motion in the protein ubiquitin. In contrast to the main chain, the motions of the methyl-bearing side chains have a large and variable pressure dependence. Within the core, this pressure sensitivity correlates with the magnitude of motion at ambient pressure. Spatial clustering of the dynamic response to applied hydrostatic pressure is also seen, indicating localized cooperativity of motion on the sub-nanosecond time scale and suggesting regions of variable compressibility. These and other features indicate that the native ensemble contains a significant fraction of members with characteristics ascribed to the recently postulated "dry molten globule". The accompanying variable side-chain conformational entropy helps complete our view of the thermodynamic architecture underlying protein stability, folding, and function.

Optimized linear prediction for radial sampled multidimensional NMR experiments.

Radial sampling in multidimensional NMR experiments offers greatly decreased acquisition times while also providing an avenue for increased sensitivity. Digital resolution remains a concern and depends strongly upon the extent of sampling of individual radial angles. Truncated time domain data leads to spurious peaks (artifacts) upon FT and 2D FT. Linear prediction is commonly employed to improve resolution in Cartesian sampled NMR experiments. Here, we adapt the linear prediction method to radial sampling. Significantly more accurate estimates of linear prediction coefficients are obtained by combining quadrature frequency components from the multiple angle spectra. This approach results in significant improvement in both resolution and removal of spurious peaks as compared to traditional linear prediction methods applied to radial sampled data. The 'averaging linear prediction' (ALP) method is demonstrated as a general tool for resolution improvement in multidimensional radial sampled experiments.