Correlating histone acetylation with nucleosome core particle dynamics and function

Unveiling the Conformational Dynamics of the Histone Tails Using Markov State Modeling

Biomolecules predominantly exert their function by altering conformational dynamics. The nucleosome core particle (NCP) is the fundamental unit of chromatin. DNA with ∼146 base pairs wraps around the histone octamer to form a nucleosome. The histone octamer is composed of two copies of each histone protein (H3, H4, H2A, and H2B) with a globular core and disordered N-terminal tails. Epigenetic modifications of the histone N-terminal tails play a critical role in regulating the chromatin structure and biological processes such as transcription and DNA repair. Here, we report all-atom molecular dynamics (MD) simulations of the nucleosome at microsecond time scales to construct Markov state models (MSMs) to elucidate distinct conformations of the histone tails. We employ time-lagged independent component analysis (tICA) to capture their essential slow dynamics, with k-means clustering used to discretize the conformational space. MSMs unveil distinct states and transition probabilities to characterize the dynamics and kinetics of the tails. Next, we focus on the H2B tail, which is one of the least studied tails. We show that acetylation increases secondary structure formation with increased transition rates. These findings will aid in understanding the functional implications of tail conformations for nucleosome stability and gene regulation.

NMR-Based Measurements of Site-Specific Electrostatic Potentials of Histone Tails in Nucleosome Core Particles

Electrostatics play a dominant role in guiding many biological processes. This is especially the case in the context of chromatin, where charge interactions modulate diverse activities such as DNA repair, transcription, replication, condensation, and phase separation. Using NMR experiments quantifying solvent paramagnetic relaxation enhancements of backbone amide and side chain methyl protons in the presence of paramagnetic cosolutes and focusing on the nucleosome core particle (NCP), we report near surface electrostatic potentials of tail residues of each of the four histone components of the NCP. These are all negative, despite sequences comprising a high density of positively charged amino acids, emphasizing the strong contribution of the negatively charged DNA with which the tails interact. Changes in electrostatic potentials of as much as 60 mV between isolated histone tails and tails in the context of the NCP are calculated. Notably, the tail potentials can vary significantly from each other, with enrichment in glycine residues correlating with less negative values, highlighting differences in the interactions with DNA.

Deep Learning of Proteins with Local and Global Regions of Disorder

Although machine learning has transformed protein structure prediction of folded protein ground states with remarkable accuracy, intrinsically disordered proteins and regions (IDPs/IDRs) are defined by diverse and dynamical structural ensembles that are predicted with low confidence by algorithms such as AlphaFold. We present a new machine learning method, IDPForge (Intrinsically Disordered Protein, FOlded and disordered Region GEnerator), that exploits a transformer protein language diffusion model to create all-atom IDP ensembles and IDR disordered ensembles that maintains the folded domains. IDPForge does not require sequence-specific training, back transformations from coarse-grained representations, nor ensemble reweighting, as in general the created IDP/IDR conformational ensembles show good agreement with solution experimental data, and options for biasing with experimental restraints are provided if desired. We envision that IDPForge with these diverse capabilities will facilitate integrative and structural studies for proteins that contain intrinsic disorder.

Connected Chromatin Amplifies Acetylation-Modulated Nucleosome Interactions

Histone acetylation is a key regulatory post-translational modification closely associated with gene transcription. In particular, H4K16 acetylation (H4K16ac) is a crucial gene activation marker that induces an open chromatin configuration. While previous studies have explored the effects of H4K16ac on nucleosome interactions, how this local modification affects higher-order chromatin organization remains unclear. To bridge the chemical modifications of these histone tail lysine residues to global chromatin structure, we utilized a residue-resolution coarse-grained chromatin model and enhanced sampling techniques to simulate charge-neutralization effects of histone acetylation on nucleosome stability, internucleosome interactions, and higher-order chromatin structure. Our simulations reveal that H4K16ac stabilizes a single nucleosome due to the reduced entropic contribution of histone tails during DNA unwrapping. In addition, acetylation modestly weakens internucleosome interactions by diminishing contacts between histone tails, DNA, and nucleosome acidic patches. These weakened interactions are amplified when nucleosomes are connected by linker DNA, where increases in linker DNA entry-exit angles lead to significant chromatin destacking and decompaction, exposing nucleosomes to transcriptional activity. Our findings suggest that the geometric constraint imposed by chromatin DNA plays a critical role in driving chromatin structural reorganization upon post-translational modifications.

Multifunctional histone variants in genome function

Histones are integral components of eukaryotic chromatin that have a pivotal role in the organization and function of the genome. The dynamic regulation of chromatin involves the incorporation of histone variants, which can dramatically alter its structural and functional properties. Contrary to an earlier view that limited individual histone variants to specific genomic functions, new insights have revealed that histone variants exert multifaceted roles involving all aspects of genome function, from governing patterns of gene expression at precise genomic loci to participating in genome replication, repair and maintenance. This conceptual change has led to a new understanding of the intricate interplay between chromatin and DNA-dependent processes and how this connection translates into normal and abnormal cellular functions.

Involvement of Histone Acetyltransferase 1 (HAT1) in the Spermatogenesis of Non-Condensed Nuclear Sperm in Chinese Mitten Crab, Eriocheir sinensis

Chinese mitten crab, Eriocheir sinensis, is a decapod crustacean with a special, non-condensated nucleus in the sperm. Studies have shown that the nuclear compact state of male germ cells during the spermatogenesis is closely related to histone modification. To explore the possible role of histone acetyltransferase 1 (HAT1) in the chromatin organization during the E. sinensis spermatogenesis, we took the testis tissues of both adult and juvenile crabs as the materials of study and analyzed the biological functions of HAT1 by whole transcriptome sequencing and bioinformatics, then further analyzed the expression and distribution of HAT1 using the methods of RT-qRCR, western blotting, and immunofluorescence location. The results showed that HAT1 is an alkaline-unstable hydrophilic protein. It was predicted to interact with a variety of histones and chromosome assembly proteins, including Asf1b, Chaf1b, and Hist1h3f, and is involved in many biological functions pertaining to chromatin dynamics such as chromatin organization, DNA dependent nucleosome assembly, DNA conformational changes, and so on. HAT1 was up-regulated in the adult testes compared to the juvenile (n = 3, P < 0.05). HAT1 was mainly located in the nuclei of male germ cells of E. sinensis. As spermatogenesis proceeded, the expression of HAT1 decreased and even disappeared in the nuclei (n = 3, P < 0.05). HAT1 is an important player in histone acetylation, which facilitates chromatin alteration in a three-dimensional conformation. The expression of HAT1 in different male germ cells might indicate the chromatin dynamics at the diversity stages of spermatogenesis. The high expression of HAT1 at the early stages of E. sinensis spermatogenesis hints the active involvement in chromatin organization, while its progressively reduced expression accompanied by the progression of spermatogenesis suggests a relatively gradual stabilization and stereotyping of chromatin. As for the disappearance of HAT1 in mature sperm with non-condensed nuclei, the reduction in histones targeted by HAT1 or histone acetylation may be an important initiator.

Unveiling the Conformational Dynamics of the Histone Tails Using Markov State Modeling

Biomolecules predominantly exert their function through altering conformational dynamics. The nucleosome core particle (NCP) is the fundamental unit of chromatin. DNA with ~146 base pairs wrap around the histone octamer to form a nucleosome. The histone octamer is comprised of two copies of each histone protein (H3, H4, H2A, and H2B) with a globular core and disordered N-terminal tails. Epigenetic modifications of the histone N-terminal tails play a critical role in the regulation of chromatin structure and biological processes such as transcription and DNA repair. Here, we report all-atomistic molecular dynamics (MD) simulations of the nucleosome at microsecond timescales to construct Markov state models (MSMs) to elucidate distinct conformations of the histone tails. We employ the time-lagged independent component analysis (tICA) to capture their essential slow dynamics, with k-means clustering used to discretize the conformational space. MSMs unveil distinct states and transition probabilities to characterize the dynamics and kinetics of the tails. Next, we focus on the H2B tail, one of the least studied tails. We show that acetylation increases secondary structure formation, with an increase in transition rates. These findings will aid in understanding the functional implications of tail conformations in nucleosome stability and gene regulation.

Histone-acetyl epigenome regulates TGF-β pathway-associated chemoresistance in colorectal cancer

TGF-β signaling pathway has been demonstrated to be closely related to chemoresistance, which is the major cause of recurrence and poor outcome in colorectal cancer (CRC), however, the comprehensive epigenetic landscape that functionally implicates in the regulation of TGF-β pathway-associated chemoresistance has not yet well established in CRC. In our study, chromatin immunoprecipitation sequencing (ChIP-seq) and Western blot were employed to investigate epigenetic modifications for histones in response to TGF-β1 intervene. We found that the activation of the TGF-β pathway was characterized by genome-wide high levels of H3K9ac and H3K18ac. Mechanistically, the activation of the TGF-β signaling pathway leads to the downregulation of the deacetylase HDAC4, resulting in the upregulation of H3K9ac and H3K18ac. Consequently, this cascade induces oxaliplatin chemoresistance in CRC by triggering the anti-apoptotic PI3K/AKT signaling pathway. Our in vivo experiment results confirmed that overexpression of HDAC4 significantly enhances the sensitivity of CRC to oxaliplatin chemotherapy. Moreover, the expression level of HDAC4 was positively correlated with patients' prognosis in CRC. Our data suggest that histone-acetyl modification demonstrates a crucial role in modulating TGF-β pathway-associated chemoresistance in CRC, and HDAC4 would be a biomarker for prognostic prediction and potential therapeutic target for treatment in CRC.

The role of epigenetics in women's reproductive health: the impact of environmental factors

This paper explores the significant role of epigenetics in women's reproductive health, focusing on the impact of environmental factors. It highlights the crucial link between epigenetic modifications-such as DNA methylation and histones post-translational modifications-and reproductive health issues, including infertility and pregnancy complications. The paper reviews the influence of pollutants like PM2.5, heavy metals, and endocrine disruptors on gene expression through epigenetic mechanisms, emphasizing the need for understanding how dietary, lifestyle choices, and exposure to chemicals affect gene expression and reproductive health. Future research directions include deeper investigation into epigenetics in female reproductive health and leveraging gene editing to mitigate epigenetic changes for improving IVF success rates and managing reproductive disorders.

Determining mesoscale chromatin structure parameters from spatially correlated cleavage data using a coarse-grained oligonucleosome model

The three-dimensional structure of chromatin has emerged as an important feature of eukaryotic gene regulation. Recent technological advances in DNA sequencing-based assays have revealed locus- and chromatin state-specific structural patterns at the length scale of a few nucleosomes (~1 kb). However, interpreting these data sets remains challenging. Radiation-induced correlated cleavage of chromatin (RICC-seq) is one such chromatin structure assay that maps DNA-DNA-contacts at base pair resolution by sequencing single-stranded DNA fragments released from irradiated cells. Here, we develop a flexible modeling and simulation framework to enable the interpretation of RICC-seq data in terms of oligonucleosome structure ensembles. Nucleosomes are modeled as rigid bodies with excluded volume and adjustable DNA wrapping, connected by linker DNA modeled as a worm-like chain. We validate the model's parameters against cryo-electron microscopy and sedimentation data. Our results show that RICC-seq is sensitive to nucleosome spacing, nucleosomal DNA wrapping, and the strength of inter-nucleosome interactions. We show that nucleosome repeat lengths consistent with orthogonal assays can be extracted from experimental RICC-seq data using a 1D convolutional neural net trained on RICC-seq signal predicted from simulated ensembles. We thus provide a suite of analysis tools that add quantitative structural interpretability to RICC-seq experiments.

Conformational Dynamics of the Nucleosomal Histone H2B Tails Revealed by Molecular Dynamics Simulations

Epigenetic modifications of histone N-terminal tails play a critical role in regulating the chromatin structure and biological processes such as transcription and DNA repair. One of the key post-translational modifications (PTMs) is the acetylation of lysine residues on histone tails. Epigenetic modifications are ubiquitous in the development of diseases, such as cancer and neurological disorders. Histone H2B tails are critical regulators of nucleosome dynamics, biological processes, and certain diseases. Here, we report all-atomistic molecular dynamics (MD) simulations of the nucleosome to demonstrate that acetylation of the histone tails changes their conformational space and interaction with DNA. We perform simulations of H2B tails, critical regulators of gene regulation, in both the lysine-acetylated (ACK) and unacetylated wild type (WT) states. To explore the effects of salt concentration, we use two different NaCl concentrations to perform simulations at microsecond time scales. Salt can modulate the effects of electrostatic interactions between the DNA phosphate backbone and histone tails. Upon acetylation, H2B tails shift their secondary structure helical propensity. The number of contacts between the DNA and the H2B tail decreases. We characterize the conformational dynamics of the H2B tails by principal component analysis (PCA). The ACK tails become more compact at increased salt concentrations, but conformations from the WT tails display the most contacts with DNA at both salt concentrations. Mainly, H2B acetylation may increase the DNA accessibility for regulatory proteins to bind, which can aid in gene regulation and NCP stability.

Structural dynamics in chromatin unraveling by pioneer transcription factors

Pioneer transcription factors are proteins with a dual function. First, they regulate transcription by binding to nucleosome-free DNA regulatory elements. Second, they bind to DNA while wrapped around histone proteins in the chromatin and mediate chromatin opening. The molecular mechanisms that connect the two functions are yet to be discovered. In recent years, pioneer factors received increased attention mainly because of their crucial role in promoting cell fate transitions that could be used for regenerative therapies. For example, the three factors required to induce pluripotency in somatic cells, Oct4, Sox2, and Klf4 were classified as pioneer factors and studied extensively. With this increased attention, several structures of complexes between pioneer factors and chromatin structural units (nucleosomes) have been resolved experimentally. Furthermore, experimental and computational approaches have been designed to study two unresolved, key scientific questions: First, do pioneer factors induce directly local opening of nucleosomes and chromatin fibers upon binding? And second, how do the unstructured tails of the histones impact the structural dynamics involved in such conformational transitions? Here we review the current knowledge about transcription factor-induced nucleosome dynamics and the role of the histone tails in this process. We discuss what is needed to bridge the gap between the static views obtained from the experimental structures and the key structural dynamic events in chromatin opening. Finally, we propose that integrating nuclear magnetic resonance spectroscopy with molecular dynamics simulations is a powerful approach to studying pioneer factor-mediated dynamics of nucleosomes and perhaps small chromatin fibers using native DNA sequences.

Cryo-EM of the Nucleosome Core Particle Bound to Ran-RCC1 Reveals a Dynamic Complex

Ras-related nuclear protein (Ran) is a member of the Ras superfamily of small guanosine triphosphatases (GTPases) and a regulator of multiple cellular processes. In healthy cells, the GTP-bound form of Ran is concentrated at chromatin, creating a Ran•GTP gradient that provides the driving force for nucleocytoplasmic transport, mitotic spindle assembly, and nuclear envelope formation. The Ran•GTP gradient is maintained by the regulator of chromatin condensation 1 (RCC1), a guanine nucleotide exchange factor that accelerates GDP/GTP exchange in Ran. RCC1 interacts with nucleosomes, which are the fundamental repeating units of eukaryotic chromatin. Here, we present a cryo-EM analysis of a trimeric complex composed of the nucleosome core particle (NCP), RCC1, and Ran. While the contacts between RCC1 and Ran in the complex are preserved compared with a previously determined structure of RCC1-Ran, our study reveals that RCC1 and Ran interact dynamically with the NCP and undergo rocking motions on the nucleosome surface. Furthermore, the switch 1 region of Ran, which plays an important role in mediating conformational changes associated with the substitution of GDP and GTP nucleotides in Ras family members, appears to undergo disorder-order transitions and forms transient contacts with the C-terminal helix of histone H2B. Nucleotide exchange assays performed in the presence and absence of NCPs are not consistent with an active role for nucleosomes in nucleotide exchange, at least in vitro. Instead, the nucleosome stabilizes RCC1 and serves as a hub that concentrates RCC1 and Ran to promote efficient Ran•GDP to Ran•GTP conversion.

Nucleosome conformation dictates the histone code

Histone post-translational modifications (PTMs) play a critical role in chromatin regulation. It has been proposed that these PTMs form localized 'codes' that are read by specialized regions (reader domains) in chromatin-associated proteins (CAPs) to regulate downstream function. Substantial effort has been made to define [CAP: histone PTM] specificities, and thus decipher the histone code and guide epigenetic therapies. However, this has largely been done using the reductive approach of isolated reader domains and histone peptides, which cannot account for any higher-order factors. Here, we show that the [BPTF PHD finger and bromodomain: histone PTM] interaction is dependent on nucleosome context. The tandem reader selectively associates with nucleosomal H3K4me3 and H3K14ac or H3K18ac, a combinatorial engagement that despite being in cis is not predicted by peptides. This in vitro specificity of the BPTF tandem reader for PTM-defined nucleosomes is recapitulated in a cellular context. We propose that regulatable histone tail accessibility and its impact on the binding potential of reader domains necessitates we refine the 'histone code' concept and interrogate it at the nucleosome level.

Reduced SIRT1 and SIRT3 and Lower Antioxidant Capacity of Seminal Plasma Is Associated with Shorter Sperm Telomere Length in Oligospermic Men

Infertility affects millions of couples worldwide and has a profound impact not only on their families, but also on communities. Telomere attrition has been associated with infertility, DNA damage and fragmentation. Oxidative stress has been shown to affect sperm DNA integrity and telomere length. Sirtuins such as SIRT1 and SIRT3 are involved in aging and oxidative stress response. The aim of the present study is to determine the role of SIRT1 and SIRT3 in regulating oxidative stress, telomere shortening, and their association with oligospermia. Therefore, we assessed the protein levels of SIRT1 and SIRT3, total antioxidant capacity (TAC), superoxide dismutase (SOD), malondialdehyde (MDA) and catalase activity (CAT) in the seminal plasma of 272 patients with oligospermia and 251 fertile men. We also measured sperm telomere length (STL) and leukocyte telomere length (LTL) using a standard real-time quantitative PCR assay. Sperm chromatin and protamine deficiency were also measured as per standard methods. Our results for oligospermic patients demonstrate significant reductions in semen parameters, shorter STL and LTL, lower levels of SOD, TAC, CAT, SIRT1 and SIRT3 levels, and also significant protamine deficiency and higher levels of MDA and DNA fragmentation. We conclude that a shorter TL in sperms and leukocytes is associated with increased oxidative stress that also accounts for high levels of DNA fragmentation in sperms. Our results support the hypothesis that various sperm parameters in the state of oligospermia are associated with or caused by reduced levels of SIRT1 and SIRT3 proteins.

Using NMR diffusion data to validate MD models of disordered proteins: Test case of N-terminal tail of histone H4

MD simulations can provide uniquely detailed models of intrinsically disordered proteins (IDPs). However, these models need careful experimental validation. The coefficient of translational diffusion Dtr, measurable by pulsed field gradient NMR, offers a potentially useful piece of experimental information related to the compactness of the IDP's conformational ensemble. Here, we investigate, both experimentally and via the MD modeling, the translational diffusion of a 25-residue N-terminal fragment from histone H4 (N-H4). We found that the predicted values of Dtr, as obtained from mean-square displacement of the peptide in the MD simulations, are largely determined by the viscosity of the MD water (which has been reinvestigated as a part of our study). Beyond that, our analysis of the diffusion data indicates that MD simulations of N-H4 in the TIP4P-Ew water give rise to an overly compact conformational ensemble for this peptide. In contrast, TIP4P-D and OPC simulations produce the ensembles that are consistent with the experimental Dtr result. These observations are supported by the analyses of the 15N spin relaxation rates. We also tested a number of empirical methods to predict Dtr based on IDP's coordinates extracted from the MD snapshots. In particular, we show that the popular approach involving the program HYDROPRO can produce misleading results. This happens because HYDROPRO is not intended to predict the diffusion properties of highly flexible biopolymers such as IDPs. Likewise, recent empirical schemes that exploit the relationship between the small-angle x-ray scattering-informed conformational ensembles of IDPs and the respective experimental Dtr values also prove to be problematic. In this sense, the first-principle calculations of Dtr from the MD simulations, such as demonstrated in this work, should provide a useful benchmark for future efforts in this area.

Conformational and Interaction Landscape of Histone H4 Tails in Nucleosomes Probed by Paramagnetic NMR Spectroscopy

The fundamental repeat unit of chromatin, the nucleosome, consists of approximately 147 base pairs of double-stranded DNA and a histone protein octamer containing two copies each of histones H2A, H2B, H3, and H4. Each histone possesses a dynamically disordered N-terminal tail domain, and it is well-established that the tails of histones H3 and H4 play key roles in chromatin compaction and regulation. Here we investigate the conformational ensemble and interactions of the H4 tail in nucleosomes by means of solution NMR measurements of paramagnetic relaxation enhancements (PREs) in recombinant samples reconstituted with 15N-enriched H4 and nitroxide spin-label tagged H3. The experimental PREs, which report on the proximities of individual H4 tail residues to the different H3 spin-label sites, are interpreted by using microsecond time-scale molecular dynamics simulations of the nucleosome core particle. Collectively, these data enable improved localization of histone H4 tails in nucleosomes and support the notion that H4 tails engage in a fuzzy complex interaction with nucleosomal DNA.

Transcriptional co-activators: emerging roles in signaling pathways and potential therapeutic targets for diseases

Specific cell states in metazoans are established by the symphony of gene expression programs that necessitate intricate synergic interactions between transcription factors and the co-activators. Deregulation of these regulatory molecules is associated with cell state transitions, which in turn is accountable for diverse maladies, including developmental disorders, metabolic disorders, and most significantly, cancer. A decade back most transcription factors, the key enablers of disease development, were historically viewed as 'undruggable'; however, in the intervening years, a wealth of literature validated that they can be targeted indirectly through transcriptional co-activators, their confederates in various physiological and molecular processes. These co-activators, along with transcription factors, have the ability to initiate and modulate transcription of diverse genes necessary for normal physiological functions, whereby, deregulation of such interactions may foster tissue-specific disease phenotype. Hence, it is essential to analyze how these co-activators modulate specific multilateral processes in coordination with other factors. The proposed review attempts to elaborate an in-depth account of the transcription co-activators, their involvement in transcription regulation, and context-specific contributions to pathophysiological conditions. This review also addresses an issue that has not been dealt with in a comprehensive manner and hopes to direct attention towards future research that will encompass patient-friendly therapeutic strategies, where drugs targeting co-activators will have enhanced benefits and reduced side effects. Additional insights into currently available therapeutic interventions and the associated constraints will eventually reveal multitudes of advanced therapeutic targets aiming for disease amelioration and good patient prognosis.

Contributions of histone tail clipping and acetylation in nucleosome transcription by RNA polymerase II

The N-terminal tails of histones protrude from the nucleosome core and are target sites for histone modifications, such as acetylation and methylation. Histone acetylation is considered to enhance transcription in chromatin. However, the contribution of the histone N-terminal tail to the nucleosome transcription by RNA polymerase II (RNAPII) has not been clarified. In the present study, we reconstituted nucleosomes lacking the N-terminal tail of each histone, H2A, H2B, H3 or H4, and performed RNAPII transcription assays. We found that the N-terminal tail of H3, but not H2A, H2B and H4, functions in RNAPII pausing at the SHL(-5) position of the nucleosome. Consistently, the RNAPII transcription assay also revealed that the nucleosome containing N-terminally acetylated H3 drastically alleviates RNAPII pausing at the SHL(-5) position. In addition, the H3 acetylated nucleosome produced increased amounts of the run-off transcript. These results provide important evidence that the H3 N-terminal tail plays a role in RNAPII pausing at the SHL(-5) position of the nucleosome, and its acetylation directly alleviates this nucleosome barrier.