Nucleosome Histone Tail Conformation and Dynamics: Impacts of Lysine Acetylation and a Nearby Minor Groove Benzo[a]pyrene-Derived Lesion

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.

Dynamic changes of DNA methylation induced by benzo(a)pyrene in cancer

Benzo(a)pyrene (BaP), the earliest and most significant carcinogen among polycyclic aromatic hydrocarbons (PAHs), has been found in foods, tobacco smoke, and automobiles exhaust, etc. Exposure to BaP induced DNA damage directly, or oxidative stress-related damage, resulting in cell apoptosis and carcinogenesis in human respiratory system, digestive system, reproductive system, etc. Moreover, BaP triggered genome-wide epigenetic alterations by methylation, which might cause disturbances in regulation of gene expression, and thereby induced cancer. It has been proved that BaP reduced genome-wide DNA methylation, and activated proto-oncogene by hypomethylation in the promoter region, but silenced tumor suppressor genes by promoter hypermethylation, resulting in cancer initiation and progression. Here we summarized the changes in DNA methylation in BaP exposure, and revealed the methylation of DNA plays a role in cancer development.

Ionic strength modulates excision of uracil by SMUG1 from nucleosome core particles

Ionic strength affects many cellular processes including the packaging of genetic material in eukaryotes. For example, chromatin fibers are compacted in high ionic strength environments as are the minimal unit of packaging in chromatin, nucleosome core particles (NCPs). Furthermore, ionic strength is known to modulate several aspects of NCP dynamics including transient unwrapping of DNA from the histone protein core, nucleosome gaping, and intra- and internucleosomal interactions of the N-terminal histone tails. Changes in NCP structure may also impact interactions of transcriptional, repair, and other cellular machinery with nucleosomal DNA. One repair process, base excision repair (BER), is impacted by NCP structure and may be further influenced by changes in ionic strength. Here we examine the effects of ionic strength on the initiation of BER using biochemical assays. Using a population of NCPs containing uracil (U) at dozens of geometric locations, excision of U by single-strand selective monofunctional uracil DNA glycosylase (SMUG1) is assessed at higher and lower ionic strengths. SMUG1 has increased excision activity in the lower ionic strength conditions. On duplex DNA, however, SMUG1 activity is largely unaffected by ionic strength except at short incubation times, suggesting that changes in SMUG1 activity are likely due to alterations in NCP structure and dynamics. These results allow us to further understand the cellular role of SMUG1 in a changing ionic environment and broadly contribute to the understanding of BER on chromatin and genomic stability.

Histone variants H3.3 and H2A.Z/H3.3 facilitate excision of uracil from nucleosome core particles

At the most fundamental level of chromatin organization, DNA is packaged as nucleosome core particles (NCPs) where DNA is wound around a core of histone proteins. This ubiquitous sequestration of DNA within NCPs presents a significant barrier to many biological processes, including DNA repair. We previously demonstrated that histone variants from the H2A family facilitate excision of uracil (U) lesions by DNA base excision repair (BER) glycosylases. Here, we consider how the histone variant H3.3 and double-variant H2A.Z/H3.3 modulate the BER enzymes uracil DNA glycosylase (UDG) and single-strand selective monofunctional uracil DNA glycosylase (SMUG1). Using an NCP model system with U:G base pairs at a wide variety of geometric positions we generate the global repair profile for both glycosylases. Enhanced excision of U by UDG and SMUG1 is observed with the H3.3 variant. We demonstrate that these H3.3-containing NCPs form two species: (1) octasomes, which contain the full complement of eight histone proteins and (2) hexasomes which are sub-nucleosomal particles that contain six histones. Both the octasome and hexasome species facilitate excision activity of UDG and SMUG1, with the largest impacts observed at sterically-occluded lesion sites and in terminal regions of DNA of the hexasome that do not closely interact with histones. For the double-variant H2A.Z/H3.3 NCPs, which exist as octasomes, the global repair profile reveals that UDG but not SMUG1 has increased U excision activity. The enhanced glycosylase activity reveals potential functions for these histone variants to facilitate BER in packaged DNA and contributes to our understanding of DNA repair in chromatin and its significance regarding mutagenesis and genomic integrity.

A Dynamic View of the Interaction of Histone Tails with Clustered Abasic Sites in a Nucleosome Core Particle

Apurinic/apyrimidinic sites are the most common forms of DNA damage under physiological conditions, yet their structural and dynamical behavior within nucleosome core particles has just begun to be investigated and is dramatically different from that of abasic sites in B-DNA. Clusters of two or more abasic sites are repaired even less efficiently and hence constitute hot spots of high mutagenicity notably due to enhanced double-strand break formation. On the basis of an X-ray structure of a 146 bp DNA wrapped onto a histone core, we investigate the structural behavior of two bistranded abasic sites positioned at mutational hot spots during microsecond-range molecular dynamics simulations. Our simulations allow us to probe interactions of histone tails at clustered abasic site locations, with a definitive assignment of the key residues involved in the NCP-catalyzed formation of DNA-protein cross-linking in line with recent experimental findings, and pave the way for a systematic assessment of the response of histone tails to DNA lesions.

The N-terminal Tails of Histones H2A and H2B Adopt Two Distinct Conformations in the Nucleosome with Contact and Reduced Contact to DNA

The nucleosome comprises two histone dimers of H2A-H2B and one histone tetramer of (H3-H4)₂, wrapped around by ~145 bp of DNA. Detailed core structures of nucleosomes have been established by X-ray and cryo-EM, however, histone tails have not been visualized. Here, we have examined the dynamic structures of the H2A and H2B tails in 145-bp and 193-bp nucleosomes using NMR, and have compared them with those of the H2A and H2B tail peptides unbound and bound to DNA. Whereas the H2A C-tail adopts a single but different conformation in both nucleosomes, the N-tails of H2A and H2B adopt two distinct conformations in each nucleosome. To clarify these conformations, we conducted molecular dynamics (MD) simulations, which suggest that the H2A N-tail can locate stably in either the major or minor grooves of nucleosomal DNA. While the H2B N-tail, which sticks out between two DNA gyres in the nucleosome, was considered to adopt two different orientations, one toward the entry/exit side and one on the opposite side. Then, the H2A N-tail minor groove conformation was obtained in the H2B opposite side and the H2B N-tail interacts with DNA similarly in both sides, though more varied conformations are obtained in the entry/exit side. Collectively, the NMR findings and MD simulations suggest that the minor groove conformer of the H2A N-tail is likely to contact DNA more strongly than the major groove conformer, and the H2A N-tail reduces contact with DNA in the major groove when the H2B N-tail is located in the entry/exit side.

Histone Tail Conformations: A Fuzzy Affair with DNA

The core histone tails are critical in chromatin structure and signaling. Studies over the past several decades have provided a wealth of information on the histone tails and their interaction with chromatin factors. However, the conformation of the histone tails in a chromatin relevant context has remained elusive. Only recently has enough evidence emerged to start to build a structural model of the tails in the context of nucleosomes and nucleosome arrays. Here, we review these studies and propose that the histone tails adopt a high-affinity fuzzy complex with DNA, characterized by robust but dynamic association. Furthermore, we discuss how these DNA-bound conformational ensembles promote distinct chromatin structure and signaling, and that their fuzzy nature is important in transitioning between functional states.

Nucleosome Core Particles Lacking H2B or H3 Tails Are Altered Structurally and Have Differential Base Excision Repair Fingerprints

A recently discovered post-translational modification of histone proteins is the irreversible proteolytic clipping of the histone N-terminal tail domains. This modification is involved in the regulation of various biological processes, including the DNA damage response. In this work, we used chemical footprinting to characterize the structural alterations to nucleosome core particles (NCPs) that result from a lack of a histone H2B or H3 tail. We also examine the influence of these histone tails on excision of the mutagenic lesion 1,N⁶-ethenoadenine (εA) by the repair enzyme alkyladenine DNA glycosylase. We found that the absence of the H2B or H3 tail results in altered DNA periodicity relative to that of native NCPs. We correlated these structural alterations to εA excision by utilizing a global analysis of 21 εA sites in NCPs and unincorporated duplex DNA. In comparison to native NCPs, there is enhanced excision of εA in tailless H2B NCPs in regions that undergo DNA unwrapping. This enhanced excision is not observed for tailless H3 NCPs; rather, excision is inhibited in more static areas of the NCP not prone to unwrapping. Our results support in vivo observations of alkylation damage profiles and the potential role of tail clipping as a mechanism for overcoming physical obstructions caused by packaging in NCPs but also reveal the potential inhibition of repair by tail clipping in some locations. Taken together, these results further our understanding of how base excision repair can be facilitated or diminished by histone tail removal and contribute to our understanding of the underlying mechanism that leads to mutational hot spots.

Breaths, Twists, and Turns of Atomistic Nucleosomes

Gene regulation programs establish cellular identity and rely on dynamic changes in the structural packaging of genomic DNA. The DNA is packaged in chromatin, which is formed from arrays of nucleosomes displaying different degree of compaction and different lengths of inter-nucleosomal linker DNA. The nucleosome represents the repetitive unit of chromatin and is formed by wrapping 145-147 basepairs of DNA around an octamer of histone proteins. Each of the four histones is present twice and has a structured core and intrinsically disordered terminal tails. Chromatin dynamics are triggered by inter- and intra-nucleosome motions that are controlled by the DNA sequence, the interactions between the histone core and the DNA, and the conformations, positions, and DNA interactions of the histone tails. Understanding chromatin dynamics requires studying all these features at the highest possible resolution. For this, molecular dynamics simulations can be used as a powerful complement or alternative to experimental approaches, from which it is often very challenging to characterize the structural features and atomic interactions controlling nucleosome motions. Molecular dynamics simulations can be performed at different resolutions, by coarse graining the molecular system with varying levels of details. Here we review the successes and the remaining challenges of the application of atomic resolution simulations to study the structure and dynamics of nucleosomes and their complexes with interacting partners.

Nucleosomal embedding reshapes the dynamics of abasic sites

Apurinic/apyrimidinic (AP) sites are the most common DNA lesions, which benefit from a most efficient repair by the base excision pathway. The impact of losing a nucleobase on the conformation and dynamics of B-DNA is well characterized. Yet AP sites seem to present an entirely different chemistry in nucleosomal DNA, with lifetimes reduced up to 100-fold, and the much increased formation of covalent DNA-protein cross-links leading to strand breaks, refractory to repair. We report microsecond range, all-atom molecular dynamics simulations that capture the conformational dynamics of AP sites and their tetrahydrofuran analogs at two symmetrical positions within a nucleosome core particle, starting from a recent crystal structure. Different behaviours between the deoxyribo-based and tetrahydrofuran-type abasic sites are evidenced. The two solvent-exposed lesion sites present contrasted extrahelicities, revealing the crucial role of the position of a defect around the histone core. Our all-atom simulations also identify and quantify the frequency of several spontaneous, non-covalent interactions between AP and positively-charged residues from the histones H2A and H2B tails that prefigure DNA-protein cross-links. Such an in silico mapping of DNA-protein cross-links gives important insights for further experimental studies involving mutagenesis and truncation of histone tails to unravel mechanisms of DPCs formation.

Molecular dynamics simulations of DNA-DNA and DNA-protein interactions

The all-atom molecular dynamics method can characterize the molecular-level interactions in DNA and DNA-protein systems with unprecedented resolution. Recent advances in computational technologies have allowed the method to reveal the unbiased behavior of such systems at the microseconds time scale, whereas enhanced sampling approaches have matured enough to characterize the interaction free energy with quantitative precision. Here, we describe recent progress toward increasing the realism of such simulations by refining the accuracy of the molecular dynamics force field, and we highlight recent application of the method to systems of outstanding biological interest.

Comparison of the Base Excision and Direct Reversal Repair Pathways for Correcting 1,N6-Ethenoadenine in Strongly Positioned Nucleosome Core Particles

1,N⁶-ethenoadenine (εA) is a mutagenic lesion and biomarker observed in numerous cancerous tissues. Two pathways are responsible for its repair: base excision repair (BER) and direct reversal repair (DRR). Alkyladenine DNA glycosylase (AAG) is the primary enzyme that excises εA in BER, generating stable intermediates that are processed by downstream enzymes. For DRR, the Fe(II)/α-ketoglutarate-dependent ALKBH2 enzyme repairs εA by direct conversion of εA to A. While the molecular mechanism of each enzyme is well understood on unpackaged duplex DNA, less is known about their actions on packaged DNA. The nucleosome core particle (NCP) forms the minimal packaging unit of DNA in eukaryotic organisms and is composed of 145–147 base pairs wrapped around a core of eight histone proteins. In this work, we investigated the activity of AAG and ALKBH2 on εA lesions globally distributed at positions throughout a strongly positioned NCP. Overall, we examined the repair of εA at 23 unique locations in packaged DNA. We observed a strong correlation between rotational positioning of εA and AAG activity but not ALKBH2 activity. ALKBH2 was more effective than AAG at repairing occluded εA lesions, but only AAG was capable of full repair of any εA in the NCP. However, notable exceptions to these trends were observed, highlighting the complexity of the NCP as a substrate for DNA repair. Modeling of binding of the repair enzymes to NCPs revealed that some of these observations can be explained by steric interference caused by DNA packaging. Specifically, interactions between ALKBH2 and the histone proteins obstruct binding to DNA, which leads to diminished activity. Taken together, these results support in vivo observations of alkylation damage profiles and contribute to our understanding of mutational hotspots.

Variable impact of conformationally distinct DNA lesions on nucleosome structure and dynamics: Implications for nucleotide excision repair

The packaging of DNA in nucleosomes presents a barrier for biological transactions including replication, transcription and repair. However, despite years of research, how the DNA is freed from the histone proteins and thereby allows the molecular machines to access the DNA remains poorly understood. We are interested in global genomic nucleotide excision repair (GG-NER). It is established that the histones are obstacles to this process, and DNA lesions are repaired less efficiently in nucleosomes than in free DNA. In the present study, we utilized molecular dynamics simulations to elucidate the nature of the distortions and dynamics imposed in the nucleosome by a set of three structually different lesions that vary in GG-NER efficiencies in free DNA, and in nucleosomes [Shafirovich, Geacintov, et. al, 2019]. Two of these are bulky lesions derived from metabolic activation of the environmental carcinogen benzo[a]pyrene, the 10R (+)-cis-anti-B[a]P-N²-dG and the stereoisomeric 10S (+)-trans-anti-B[a]P-N²-dG, which respectively adopt base-displaced/intercalated and minor groove-aligned conformations in DNA. The third is a non-bulky lesion, the 5'R-8-cyclo-2'-deoxyguanosine cross-link, produced by reactive oxygen and nitrogen species; cyclopurine lesions are highly mutagenic. These adducts are placed near the dyad axis, and rotationally with the lesion-containing strand facing towards or away from the histones. While each lesion has distinct conformational characteristics that are retained in the nucleosome, a spectrum of structural and dynamic disturbances, from slight to substantial, are displayed that depend on the lesion's structure and position in the nucleosome. We hypothesize that these intrinsic structural and dynamic distinctions provide different signals to initiate the cascade of chromatin-opening processes, including acetylation and other post translational modifications, remodeling by ATP-dependent complexes and spontaneous unwrapping that regulate the rate of access to the lesion; this may translate ultimately into varying GG-NER efficiencies, including repair resistance when signals for access are too weak.

Diiron Dithiolate Complex Induced Helical Structure of Histone and Application in Photochemical Hydrogen Generation

Very-lysine-rich calf thymus histone proteins form disordered structure and hydrophobic interaction-driven aggregates in weakly acidic solution. We reported that the conjugation of diiron dithiolate complex to the lysine residues induced formation of helical conformation and condensed nanoassemblies with a high loading capacity up to 18.7 wt %. The incorporated diiron dithiolate complex showed photocatalytic activity for hydrogen evolution in aqueous solutions, with a turnover number (based on [FeFe] catalyst moiety) up to 359 that was more than 6 times that of the free catalyst. The increase of helical conformation in proteins was well correlated to the increasing enhancement of photocatalytic activity. We demonstrated that the [FeFe]-hydrogenase-mimic biohybrid system based on the photocatalyst-induced protein conformational conversion and reassembly is efficient for hydrogen generation regardless of the relatively large size.

Rotational and translational positions determine the structural and dynamic impact of a single ribonucleotide incorporated in the nucleosome

Ribonucleotides misincorporated by replicative DNA polymerases are by far the most common DNA lesion. The presence of ribonucleotides in DNA is associated with genome instability, causing replication stress, chromosome fragility, gross chromosomal rearrangements, and other mutagenic events. Furthermore, nucleosome and chromatin assembly as well as nucleosome positioning are affected by the presence of ribonucleotides. Notably, nucleosome formation is significantly reduced by a single ribonucleotide. Single ribonucleotides are primarily removed from DNA by the ribonucleotide excision repair (RER) pathway via the RNase H2 enzyme, which incises the DNA backbone on the 5'-side of the ribonucleotide. While the structural implications of a single ribonucleotide in free duplex DNA have been well studied, how a single ribonucleotide embedded in nucleosomal DNA impacts nucleosome structure and dynamics, and the possible consequent impact on RER, have not been explored. We have carried out 3.5 μs molecular dynamics simulations of a single ribonucleotide incorporated at various translational and rotational positions in a nucleosome core particle. We find that the presence of the 2'-OH group on the ribose impacts the local conformation and dynamics of both the ribonucleotide and nearby DNA nucleotides as well as their interactions with histones; the nature of these disturbances depends on the rotational and translational setting, including whether the ribose faces toward or away from the histones. The ribonucleotide's preferred C3'-endo pucker is stabilized by interactions with the histones, and furthermore the ribonucleotide can cause dynamic local duplex disturbance involving an abnormal C3'-endo population of the adjacent deoxyribose pucker, minor groove opening, ruptured Watson-Crick pairing, and duplex unwinding that are governed by translation-dependent histone-nucleotide interactions. Possible effects of these disturbances on RER are considered.

Initiating base excision repair in chromatin

The base excision repair (BER) pathway removes modified nucleobases that can be deleterious to an organism. BER is initiated by a glycosylase, which finds and removes these modified nucleobases. Most of the characterization of glycosylase activity has been conducted in the context of DNA oligomer substrates. However, DNA within eukaryotic organisms exists in a packaged environment with the basic unit of organization being the nucleosome core particle (NCP). The NCP is a complex substrate for repair in which a variety of factors can influence glycosylase activity. In this Review, we focus on the geometric positioning of modified nucleobases in an NCP and the consequences on glycosylase activity and initiating BER.

Synergistic effects of H3 and H4 nucleosome tails on structure and dynamics of a lesion-containing DNA: Binding of a displaced lesion partner base to the H3 tail for GG-NER recognition

How DNA lesions in nucleosomes are recognized for global genome nucleotide excision repair (GG-NER) remains poorly understood, and the roles that histone tails may play remains to be established. Histone H3 and H4 N-terminal tails are of particular interest as their acetylation states are important in regulating nucleosomal functions in transcription, replication and repair. In particular the H3 tail has been the focus of recent attention as a site for the interaction with XPC, the GG-NER lesion recognition factor. Here we have investigated how the structure and dynamics of the DNA lesion cis-B[a]P-dG, derived from the environmental carcinogen benzo[a]pyrene (B[a]P), is impacted by the presence of flanking H3 and H4 tails. This lesion is well-repaired by GG-NER, and adopts a base-displaced/intercalated conformation in which the lesion partner C is displaced into the major groove. We used molecular dynamics simulations to obtain structural and dynamic characterizations for this lesion positioned in nucleosomal DNA so that it is bracketed by the H3 and H4 tails. The H4 tail was studied in unacetylated and acetylated states, while the H3 tail was unacetylated, its state when it binds XPC (Kakumu, Nakanishi et al., 2017). Our results reveal that upon acetylation, the H4 tail is released from the DNA surface; the H3 tail then forms a pocket that induces flipping and capture of the displaced lesion partner base C. This reveals synergistic effects of the behavior of the two tails. We hypothesize that the dual capability of the H3 tail to sense the displaced lesion partner base and to bind XPC could foster recognition of this lesion by XPC for initiation of GG-NER in nucleosomes.

Human OGG1 activity in nucleosomes is facilitated by transient unwrapping of DNA and is influenced by the local histone environment

If unrepaired, damage to genomic DNA can cause mutations and/or be cytotoxic. Single base lesions are repaired via the base excision repair (BER) pathway. The first step in BER is the recognition and removal of the nucleobase lesion by a glycosylase enzyme. For example, human oxoguanine glycosylase 1 (hOGG1) is responsible for removal of the prototypic oxidatively damaged nucleobase, 8-oxo-7,8-dihydroguanine (8-oxoG). To date, most studies of glycosylases have used free duplex DNA substrates. However, cellular DNA is packaged as repeating nucleosome units, with 145 base pair segments of DNA wrapped around histone protein octamers. Previous studies revealed inhibition of hOGG1 at the nucleosome dyad axis and in the absence of chromatin remodelers. In this study, we reveal that even in the absence of chromatin remodelers or external cofactors, hOGG1 can initiate BER at positions off the dyad axis and that this activity is facilitated by spontaneous and transient unwrapping of DNA from the histones. Additionally, we find that solution accessibility as determined by hydroxyl radical footprinting is not fully predictive of glycosylase activity and that histone tails can suppress hOGG1 activity. We therefore suggest that local nuances in the nucleosome environment and histone-DNA interactions can impact glycosylase activity.