Structure and binding of the H4 histone tail and the effects of lysine 16 acetylation

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.

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.

Modulation of cellular processes by histone and non-histone protein acetylation

Lysine acetylation is a widespread and versatile protein post-translational modification. Lysine acetyltransferases and lysine deacetylases catalyse the addition or removal, respectively, of acetyl groups at both histone and non-histone targets. In this Review, we discuss several features of acetylation and deacetylation, including their diversity of targets, rapid turnover, exquisite sensitivity to the concentrations of the cofactors acetyl-CoA, acyl-CoA and NAD⁺, and tight interplay with metabolism. Histone acetylation and non-histone protein acetylation influence a myriad of cellular and physiological processes, including transcription, phase separation, autophagy, mitosis, differentiation and neural function. The activity of lysine acetyltransferases and lysine deacetylases can, in turn, be regulated by metabolic states, diet and specific small molecules. Histone acetylation has also recently been shown to mediate cellular memory. These features enable acetylation to integrate the cellular state with transcriptional output and cell-fate decisions.

Is the H4 histone tail intrinsically disordered or intrinsically multifunctional?

The structural versatility of histone tails is one of the key elements in the organisation of chromatin, which allows for the compact storage of genomic information. However, this structural diversity also complicates experimental and computational studies. Here, the potential and free energy landscape for the isolated and bound H4 histone tail are explored. The landscapes exhibit a set of distinct structural ensembles separated by high energy barriers, with little difference between isolated and bound tails. This consistency is a desirable feature that facilitates the formation of transient interactions, which are required for the liquid-like chromatin organisation. The existence of multiple, distinct structures on a multifunnel energy landscape is likely to be associated with multifunctionality, i.e. a set of evolved, distinct functions. Contrasting it with previously reported results for other disordered peptides, this type of landscape may be associated with a conformational selection based binding mechanism. Given the similarity to other systems exhibiting similar multifunnel energy landscapes, the disorder in histone tails might be better described in context of multifunctionality.

Structural Alterations of Histone Proteins in DNA-Damaged Cells Revealed by Synchrotron Radiation Circular Dichroism Spectroscopy: A New Piece of the DNA-Damage-Response Puzzle

Double-strand breaks of DNA may lead to discontinuous DNA and consequent loss of genetic information, which may result in mutations or, ultimately, carcinogenesis. To avoid such potentially serious situations, cells have evolved efficient DNA damage repair systems. It is thought that DNA-repair processes involve drastic alterations of chromatin and histone structures, but detection of these altered structures in DNA-damaged cells remains rare in the literature. Recently, synchrotron radiation circular dichroism (SRCD) spectroscopy, which can provide secondary structural information of proteins in solution, has identified structural alterations of histone proteins induced by DNA damage responses. In this review, these results and experimental procedures are discussed with the aim of facilitating further studies of the chromatin remodeling and DNA damage repair pathways using SRCD spectroscopy.

Chromatin regulation by Histone H4 acetylation at Lysine 16 during cell death and differentiation in the myeloid compartment

Histone H4 acetylation at Lysine 16 (H4K16ac) is a key epigenetic mark involved in gene regulation, DNA repair and chromatin remodeling, and though it is known to be essential for embryonic development, its role during adult life is still poorly understood. Here we show that this lysine is massively hyperacetylated in peripheral neutrophils. Genome-wide mapping of H4K16ac in terminally differentiated blood cells, along with functional experiments, supported a role for this histone post-translational modification in the regulation of cell differentiation and apoptosis in the hematopoietic system. Furthermore, in neutrophils, H4K16ac was enriched at specific DNA repeats. These DNA regions presented an accessible chromatin conformation and were associated with the cleavage sites that generate the 50 kb DNA fragments during the first stages of programmed cell death. Our results thus suggest that H4K16ac plays a dual role in myeloid cells as it not only regulates differentiation and apoptosis, but it also exhibits a non-canonical structural role in poising chromatin for cleavage at an early stage of neutrophil cell death.

Biophysics of Chromatin Dynamics

Nucleosomes and chromatin control eukaryotic genome accessibility and thereby regulate DNA processes, including transcription, replication, and repair. Conformational dynamics within the nucleosome and chromatin structure play a key role in this regulatory function. Structural fluctuations continuously expose internal DNA sequences and nucleosome surfaces, thereby providing transient access for the nuclear machinery. Progress in structural studies of nucleosomes and chromatin has provided detailed insight into local chromatin organization and has set the stage for recent in-depth investigations of the structural dynamics of nucleosomes and chromatin fibers. Here, we discuss the dynamic processes observed in chromatin over different length scales and timescales and review current knowledge about the biophysics of distinct structural transitions.

General Purpose Water Model Can Improve Atomistic Simulations of Intrinsically Disordered Proteins

Unconstrained atomistic simulations of intrinsically disordered proteins and peptides (IDP) remain a challenge: widely used, "general purpose" water models tend to favor overly compact structures relative to experiment. Here we have performed a total of 93 μs of unrestrained MD simulations to explore, in the context of IDPs, a recently developed "general-purpose" 4-point rigid water model OPC, which describes liquid state of water close to experiment. We demonstrate that OPC, together with a popular AMBER force field ff99SB, offers a noticeable improvement over TIP3P in producing more realistic structural ensembles of three common IDPs benchmarks: 55-residue apo N-terminal zinc-binding domain of HIV-1 integrase ("protein IN"), amyloid β-peptide (Aβ₄₂) (residues 1-42), and 26-reside H4 histone tail. As a negative control, computed folding profile of a regular globular miniprotein (CLN025) in OPC water is in appreciably better agreement with experiment than that obtained in TIP3P, which tends to overstabilize the compact native state relative to the extended conformations. We employed Aβ₄₂ peptide to investigate the possible influence of the solvent box size on simulation outcomes. We advocate a cautious approach for simulations of IDPs: we suggest that the solvent box size should be at least four times the radius of gyration of the random coil corresponding to the IDP. The computed free energy landscape of protein IN in OPC resembles a shallow "tub" - conformations with substantially different degrees of compactness that are within 2 k_B T of each other. Conformations with very different secondary structure content coexist within 1 k_B T of the global free energy minimum. States with higher free energy tend to have less secondary structure. Computed low helical content of the protein has virtually no correlation with its degree of compactness, which calls into question the possibility of using the helicity as a metric for assessing performance of water models for IDPs, when the helicity is low. Predicted radius of gyration ( R _g) of H4 histone tail in OPC water falls in-between that of a typical globular protein and a fully denatured protein of the same size; the predicted R _g is consistent with two independent predictions. In contrast, H4 tail in TIP3P water is as compact as the corresponding globular protein. The computed free energy landscape of H4 tail in OPC is relatively flat over a significant range of compactness, which, we argue, is consistent with its biological function as facilitator of internucleosome interactions.

Sample Volume Reduction Using the Schwarzschild Objective for a Circular Dichroism Spectrophotometer and an Application to the Structural Analysis of Lysine-36 Trimethylated Histone H3 Protein

With the increasing interest in scarce proteins, reducing the sample volume for circular dichroism (CD) spectroscopy has become desirable. Demagnification of the incident beam size is required to reduce the sample volume for CD spectroscopy detecting transmitted light passed through the sample. In this study, the beam size was demagnified using a focal mirror, and small-capacity sample cells were developed in an attempt to reduce the sample volume. The original beam size was 6 × 6 mm²; we successfully converged it to a size of 25 × 25 μm² using the Schwarzschild objective (SO). The new sample cell and SO allowed the required sample volume to be reduced to 1/10 (15 → 1.5 μL), when using a 15 μm path length cell. By adopting a smaller sample cell, further sample reduction could be achieved. By using the SO system, the secondary structural contents of the lysine-36 trimethylated histone H3 protein were analyzed. The trimethylation induced the increment of helix structures and decrement of unordered structures. These structural alterations may play a role in regulating cellular function(s), such as DNA damage repair processes.

Inter-nucleosomal communication between histone modifications for nucleosome phasing

Combinatorial effects of epigenetic modifications on transcription activity have been proposed as “histone codes”. However, it is unclear whether there also exist inter-nucleosomal communications among epigenetic modifications at single nucleosome level, and if so, what functional roles they play. Meanwhile, how clear nucleosome patterns, such as nucleosome phasing and depletion, are formed at functional regions remains an intriguing enigma. To address these questions, we developed a Bayesian network model for interactions among different histone modifications across neighboring nucleosomes, based on the framework of dynamic Bayesian network (DBN). From this model, we found that robust inter-nucleosomal interactions exist around transcription start site (TSS), transcription termination sites (TTS) or around CTCF binding sites; and these inter-nucleosomal interactions are often involved in transcription regulation. In addition to these general principles, DBN also uncovered a novel specific epigenetic interaction between H2A.Z and H4K20me1 on neighboring nucleosomes, involved in nucleosome free region (NFR) and nucleosome phasing establishment or maintenance. The level of negative correlation between neighboring H2A.Z and H4K20me1 strongly correlate with the size of NFR and the strength of nucleosome phasing around TSS. Our study revealed inter-nucleosomal communications as important players in signal propagation, chromatin remodeling and transcription regulation.

Nucleosomes are the basic unit of chromatin organization. At a global level, they fold up to form chromatin fibers in higher order structure to control the activation/repression states of chromatins. At a local level, especially around transcriptional starting sites (TSSs), nucleosomes play an important role in regulating gene expression by dynamically positioning to affect the recruitment of RNA polymerase II and transcriptional factors. In particular around actively transcribed TSSs, nucleosomes are regularly positioned to form a typical pattern of nucleosome phasing. As it suggests that the forming of nucleosome phasing is a synergistic behavior across the nucleosomes around TSS, we hypothesize that there exist communications, which is probably some propagations of histone modifications, between neighboring nucleosomes, as nucleosome functions are essentially due to histone modifications. Here, to address the question, we investigated the correlations of histone modifications across neighboring nucleosomes, and revealed a negative correlation between H2A.Z and H4K20me1 across neighboring nucleosomes. It is a development to the well accepted knowledge that H2A.Z and H4K20me1 are positively correlated at genome-wide level. In addition, we revealed a probable contribution of H2A.Z-H4K20me1 anti-correlation in nucleosome phasing around active TSSs, therefore, shedding light on understanding the forming of nucleosome phasing.

The effect of epigenetic modifications on the secondary structures and possible binding positions of the N-terminal tail of histone H3 in the nucleosome: a computational study

The roles of histone tails as substrates for reversible chemical modifications and dynamic cognate surfaces for the binding of regulatory proteins are well established. Despite these crucial roles, experimentally derived knowledge of the structure and possible binding sites of histone tails in chromatin is limited. In this study, we utilized molecular dynamics of isolated histone H3 N-terminal peptides to investigate its structure as a function of post-translational modifications that are known to be associated with defined chromatin states. We observed a structural preference for α-helices in isoforms associated with an inactive chromatin state, while isoforms associated with active chromatin states lacked α-helical content. The physicochemical effect of the post-translational modifications was highlighted by the interaction of arginine side-chains with the phosphorylated serine residues in the inactive isoform. We also showed that the isoforms exhibit different tail lengths, and, using molecular docking of the first 15 N-terminal residues of an H3 isoform, identified potential binding sites between the superhelical gyres on the octamer surface, close to the site of DNA entry/exit in the nucleosome. We discuss the possible functional role of the binding of the H3 tail within the nucleosome on both nucleosome and chromatin structure and stability.

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

Histone tails in nucleosomes play critical roles in regulation of many biological processes, including chromatin compaction, transcription, and DNA repair. Moreover, post-translational modifications, notably lysine acetylation, are crucial to these functions. While the tails have been intensively studied, how the structures and dynamics of tails are impacted by the presence of a nearby bulky DNA lesion is a frontier research area, and how these properties are impacted by tail lysine acetylation remains unexplored. To obtain molecular insight, we have utilized all atom 3 μs molecular dynamics simulations of nucleosome core particles (NCPs) to determine the impact of a nearby DNA lesion, 10S (+)-trans-anti-B[a]P-N²-dG-the major adduct derived from the procarcinogen benzo[a]pyrene-on H2B tail behavior in unacetylated and acetylated states. We similarly studied lesion-free NCPs to investigate the normal properties of the H2B tail in both states. In the lesion-free NCPs, charge neutralization upon lysine acetylation causes release of the tail from the DNA. When the lesion is present, it stably engulfs part of the nearby tail, impairing the interactions between DNA and tail. With the tail in an acetylated state, the lesion still interacts with part of it, although unstably. The lesion's partial entrapment of the tail should hinder the tail from interacting with other nucleosomes, and other proteins such as acetylases, deacetylases, and acetyl-lysine binding proteins, and thus disrupt critical tail-governed processes. Hence, the lesion would impede tail functions modulated by acetylation or deacetylation, causing aberrant chromatin structures and impaired biological transactions such as transcription and DNA repair.

Regulation of Nucleosome Stacking and Chromatin Compaction by the Histone H4 N-Terminal Tail-H2A Acidic Patch Interaction

Chromatin folding and dynamics are critically dependent on nucleosome-nucleosome interactions with important contributions from internucleosome binding of the histone H4 N-terminal tail K16-R23 domain to the surface of the H2A/H2B dimer. The H4 Lys16 plays a pivotal role in this regard. Using in vitro reconstituted 12-mer nucleosome arrays, we have investigated the mechanism of the H4 N-terminal tail in maintaining nucleosome-nucleosome stacking and mediating intra- and inter-array chromatin compaction, with emphasis on the role of K16 and the positive charge region, R17-R23. Analytical ultracentrifugation sedimentation velocity experiments and precipitation assays were employed to analyze effects on chromatin folding and self-association, respectively. Effects on chromatin folding caused by various mutations and modifications at position K16 in the H4 histone were studied. Additionally, using charge-quenching mutations, we characterized the importance of the interaction of the residues within the H4 positive charge region R17-R23 with the H2A acidic patch of the adjacent nucleosome. Furthermore, crosslinking experiments were conducted to establish the proximity of the basic tail region to the acidic patch. Our data indicate that the positive charge and length of the side chain of H4 K16 are important for its access to the adjacent nucleosome in the process of nucleosome-nucleosome stacking and array folding. The location and orientation of the H4 R17-R23 domain on the H2A/H2B dimer surface of the neighboring nucleosome core particle (NCP) in the compacted chromatin fiber were established. The dominance of electrostatic interactions in maintaining intra-array interaction was demonstrated.

Histone Acetylation Regulates Chromatin Accessibility: Role of H4K16 in Inter-nucleosome Interaction

The N-terminal tail of histone H4 is an indispensable mediator for inter-nucleosome interaction, which is required for chromatin fiber condensation. H4K16 acetylation (H4K16Ac) activates gene transcription by influencing both chromatin structure and interplay with nonhistone proteins. To understand the influence of H4K16Ac on inter-nucleosome interaction, we performed a simulation study for the H4 tail in the context of two nucleosomes in neighboring unit cells in the crystal structure. The binding conformation of H4 tail with/without K16Ac was sampled by replica exchange with solute tempering, and the free energy landscape was explored by metadynamics. The results indicate two important features of H4K16: 1) it is the first button to anchor the H4 tail on the adjacent nucleosome; and 2) it is the only acetylation site interacting with the acidic patch. H4K16Ac disrupts the electrostatic interactions of K16, weakens H4 tail-acidic patch binding, and significantly increases H4 tail conformation diversity. Our study suggests that H4K16Ac directly reduces the inter-nucleosome interaction mediated by the H4 tail, which might further encourage the binding of nonhistone proteins on the acidic patch.

Distinct Roles of Histone H3 and H2A Tails in Nucleosome Stability

Nucleosome breathing potentially increases the DNA exposure, which in turn recruits DNA-binding protein and regulates gene transcription. Numerous studies have shown the critical roles of N-terminal tails of histones H3 and H4 in gene expression; however, few studies have focused on the H2A C-terminal tail. Here we present thorough computational studies on a single nucleosome particle showing the linker DNA closing and opening, which is thought to be nucleosome breathing. With our simulation, the H2A C-terminal and H3 N-terminal tails were found to modulate the nucleosome conformation differently. The H2A C-terminal tail regulates nucleosome conformation by binding to linker DNA at different locations, whereas the H3 N-terminal tail regulates linker DNA by binding to it in different patterns. Further MD simulation on tail truncated structures corroborates this analysis. These findings replenish our understanding of the histone tail regulation mechanism on atomic level.

Mechanism for verification of mismatched and homoduplex DNAs by nucleotides-bound MutS analyzed by molecular dynamics simulations

In order to understand how MutS recognizes mismatched DNA and induces the reaction of DNA repair using ATP, the dynamics of the complexes of MutS (bound to the ADP and ATP nucleotides, or not) and DNA (with mismatched and matched base-pairs) were investigated using molecular dynamics simulations. As for DNA, the structure of the base-pairs of the homoduplex DNA which interacted with the DNA recognition site of MutS was intermittently disturbed, indicating that the homoduplex DNA was unstable. As for MutS, the disordered loops in the ATPase domains, which are considered to be necessary for the induction of DNA repair, were close to (away from) the nucleotide-binding sites in the ATPase domains when the nucleotides were (not) bound to MutS. This indicates that the ATPase domains changed their structural stability upon ATP binding using the disordered loop. Conformational analysis by principal component analysis showed that the nucleotide binding changed modes which have structurally solid ATPase domains and the large bending motion of the DNA from higher to lower frequencies. In the MutS-mismatched DNA complex bound to two nucleotides, the bending motion of the DNA at low frequency modes may play a role in triggering the formation of the sliding clamp for the following DNA-repair reaction step. Moreover, MM-PBSA/GBSA showed that the MutS-homoduplex DNA complex bound to two nucleotides was unstable because of the unfavorable interactions between MutS and DNA. This would trigger the ATP hydrolysis or separation of MutS and DNA to continue searching for mismatch base-pairs. Proteins 2016; 84:1287-1303. © 2016 Wiley Periodicals, Inc.

Molecular dynamics simulations demonstrate the regulation of DNA-DNA attraction by H4 histone tail acetylations and mutations

The positively charged N-terminal histone tails play a crucial role in chromatin compaction and are important modulators of DNA transcription, recombination, and repair. The detailed mechanism of the interaction of histone tails with DNA remains elusive. To model the unspecific interaction of histone tails with DNA, all-atom molecular dynamics (MD) simulations were carried out for systems of four DNA 22-mers in the presence of 20 or 16 short fragments of the H4 histone tail (variations of the 16-23 a. a. KRHRKVLR sequence, as well as the unmodified fragment a. a.13-20, GGAKRHRK). This setup with high DNA concentration, explicit presence of DNA-DNA contacts, presence of unstructured cationic peptides (histone tails) and K(+) mimics the conditions of eukaryotic chromatin. A detailed account of the DNA interactions with the histone tail fragments, K(+) and water is presented. Furthermore, DNA structure and dynamics and its interplay with the histone tail fragments binding are analysed. The charged side chains of the lysines and arginines play major roles in the tail-mediated DNA-DNA attraction by forming bridges and by coordinating to the phosphate groups and to the electronegative sites in the minor groove. Binding of all species to DNA is dynamic. The structure of the unmodified fully-charged H4 16-23 a.a. fragment KRHRKVLR is dominated by a stretched conformation. The H4 tail a. a. fragment GGAKRHRK as well as the H4 Lys16 acetylated fragment are highly flexible. The present work allows capturing typical features of the histone tail-counterion-DNA structure, interaction and dynamics.

Cis and trans internucleosomal interactions of H3 and H4 tails in tetranucleosomes

Chromatin structure and the transcriptional state of a gene can be modulated by modifications made on H3 and H4 tails of histones. Elucidating the internucleosomal interactions of these tails is vital to understanding epigenetic regulation. Differentiation between cis (intra-nucleosomal) and trans (inter-nucleosomal) interactions is often difficult with conventional techniques since H3 and H4 tails are flexible. The distinction, however, is important because these interactions model short- and long-range chromatin interactions respectively and have different bearings in biological processes. Combining FCS and PCH analysis, we can decouple the contribution of histone tails to cis and trans effects. A Mg(2+) gradient was employed to facilitate compaction and oligomerization of tetranucleosomes with H3 and/or H4 tail truncations. H3 tails were found to play a multifunctional role and exhibit the ability to partake in both attractive cis and trans interactions simultaneously. H4 tails partake in attractive cis and repulsive trans interactions at low [Mg(2+)]. These interactions are diminished at higher [Mg(2+)]. Simultaneous H3 and H4 tail truncation inhibited array oligomerization but maintained local array compaction at relatively high [Mg(2+)]. The established experimental approach can be extended to study the detailed molecular interactions mediated by epigenetic modification of flexible histone tails.

The role of histone tails in the nucleosome: a computational study

Histone tails play an important role in gene transcription and expression. We present here a systematic computational study of the role of histone tails in the nucleosome, using replica exchange molecular dynamics simulations with an implicit solvent model and different well-established force fields. We performed simulations for all four histone tails, H4, H3, H2A, and H2B, isolated and with inclusion of the nucleosome. The results confirm predictions of previous theoretical studies for the secondary structure of the isolated tails but show a strong dependence on the force field used. In the presence of the entire nucleosome for all force fields, the secondary structure of the histone tails is destabilized. Specific contacts are found between charged lysine and arginine residues and DNA phosphate groups and other binding sites in the minor and major DNA grooves. Using cluster analysis, we found a single dominant configuration of binding to DNA for the H4 and H2A histone tails, whereas H3 and H2B show multiple binding configurations with an equal probability. The leading stabilizing contribution for those binding configurations is the attractive interaction between the positively charged lysine and arginine residues and the negatively charged phosphate groups, and thus the resulting charge neutralization. Finally, we present results of molecular dynamics simulations in explicit solvent to confirm our conclusions. Results from both implicit and explicit solvent models show that large portions of the histone tails are not bound to DNA, supporting the complex role of these tails in gene transcription and expression and making them possible candidates for binding sites of transcription factors, enzymes, and other proteins.

Structural insights of nucleosome and the 30-nm chromatin fiber

The eukaryotic genome is hierarchically packaged into chromatin in the nucleus. The organization and dynamics of 30-nm chromatin fibers, which is typically regarded as the secondary structure of chromatin, play a crucial role in regulating DNA accessibility for gene expression. Here we reviewed some recent progresses on the structural studies on nucleosomes, nucleosome-protein complexes, and chromatin fibers, focusing on the structural insights how the chromatin structure is regulated by different epigenetic regulation factors.

Recognition of the nucleosome by chromatin factors and enzymes

Dynamic expression of the genome requires coordinated binding of chromatin factors and enzymes that carry out genome-templated processes. Until recently, the molecular mechanisms governing how these factors and enzymes recognize and act on the fundamental unit of chromatin, the nucleosome core particle, have remained a mystery. A small, yet growing set of structures of the nucleosome in complex with chromatin factors and enzymes highlights the importance of multivalency in defining nucleosome binding and specificity. Many such interactions include an arginine anchor motif, which targets a unique acidic patch on the nucleosome surface. These emerging paradigms for chromatin recognition will be discussed, focusing on several recent structural breakthroughs.

Chromatin Unfolding by Epigenetic Modifications Explained by Dramatic Impairment of Internucleosome Interactions: A Multiscale Computational Study

Histone tails and their epigenetic modifications play crucial roles in gene expression regulation by altering the architecture of chromatin. However, the structural mechanisms by which histone tails influence the interconversion between active and inactive chromatin remain unknown. Given the technical challenges in obtaining detailed experimental characterizations of the structure of chromatin, multiscale computations offer a promising alternative to model the effect of histone tails on chromatin folding. Here we combine multimicrosecond atomistic molecular dynamics simulations of dinucleosomes and histone tails in explicit solvent and ions, performed with three different state-of-the-art force fields and validated by experimental NMR measurements, with coarse-grained Monte Carlo simulations of 24-nucleosome arrays to describe the conformational landscape of histone tails, their roles in chromatin compaction, and the impact of lysine acetylation, a widespread epigenetic change, on both. We find that while the wild-type tails are highly flexible and disordered, the dramatic increase of secondary-structure order by lysine acetylation unfolds chromatin by decreasing tail availability for crucial fiber-compacting internucleosome interactions. This molecular level description of the effect of histone tails and their charge modifications on chromatin folding explains the sequence sensitivity and underscores the delicate connection between local and global structural and functional effects. Our approach also opens new avenues for multiscale processes of biomolecular complexes.

The chromatin fiber: multiscale problems and approaches

The structure of chromatin, affected by many factors from DNA linker lengths to posttranslational modifications, is crucial to the regulation of eukaryotic cells. Combined experimental and computational methods have led to new insights into its structural and dynamical features, from interactions due to the flexible core histone tails or linker histones to the physical mechanism driving the formation of chromosomal domains. Here we present a perspective of recent advances in chromatin modeling techniques at the atomic, mesoscopic, and chromosomal scales with a view toward developing multiscale computational strategies to integrate such findings. Innovative modeling methods that connect molecular to chromosomal scales are crucial for interpreting experiments and eventually deciphering the complex dynamic organization and function of chromatin in the cell.

Genetically encoding lysine modifications on histone H4

Post-translational modifications of proteins are important modulators of protein function. In order to identify the specific consequences of individual modifications, general methods are required for homogeneous production of modified proteins. The direct installation of modified amino acids by genetic code expansion facilitates the production of such proteins independent of the knowledge and availability of the enzymes naturally responsible for the modification. The production of recombinant histone H4 with genetically encoded modifications has proven notoriously difficult in the past. Here, we present a general strategy to produce histone H4 with acetylation, propionylation, butyrylation, and crotonylation on lysine residues. We produce homogeneous histone H4 containing up to four simultaneous acetylations to analyze the impact of the modifications on chromatin array compaction. Furthermore, we explore the ability of antibodies to discriminate between alternative lysine acylations by incorporating these modifications in recombinant histone H4.

The histone H3 N-terminal tail: a computational analysis of the free energy landscape and kinetics

Extensive molecular dynamics simulations and Markov State models are used to characterize the free energy landscape and kinetics of the histone H3 N-terminal tail, which plays a critical role in regulating chromatin dynamics and gene activity.

Biomimetic Material-Assisted Delivery of Human Embryonic Stem Cell Derivatives for Enhanced In Vivo Survival and Engraftment

The ability of human embryonic stem cells (hESCs) and their derivatives to differentiate and contribute to tissue repair has enormous potential to treat various debilitating diseases. However, improving the in vivo viability and function of the transplanted cells, a key determinant of translating cell-based therapies to the clinic, remains a daunting task. Here, we develop a hybrid biomaterial consisting of hyaluronic acid (HA) grafted with 6-aminocaproic acid moieties (HA-6ACA) to improve cell delivery and their subsequent in vivo function using skeletal muscle as a model system. Our findings show that the biomimetic material-assisted delivery of hESC-derived myogenic progenitor cells into cardiotoxin-injured skeletal muscles of NOD/SCID mice significantly promotes survival and engraftment of transplanted cells in a dose-dependent manner. The donor cells were found to contribute to the regeneration of damaged muscle fibers and to the satellite cell (muscle specific stem cells) compartment. Such biomimetic cell delivery vehicles that are cost-effective and easy-to-synthesize could play a key role in improving the outcomes of other stem cell-based therapies.

Molecular interactions and residues involved in force generation in the T4 viral DNA packaging motor

Many viruses utilize molecular motors to package their genomes into preformed capsids. A striking feature of these motors is their ability to generate large forces to drive DNA translocation against entropic, electrostatic, and bending forces resisting DNA confinement. A model based on recently resolved structures of the bacteriophage T4 motor protein gp17 suggests that this motor generates large forces by undergoing a conformational change from an extended to a compact state. This transition is proposed to be driven by electrostatic interactions between complementarily charged residues across the interface between the N- and C-terminal domains of gp17. Here we use atomistic molecular dynamics simulations to investigate in detail the molecular interactions and residues involved in such a compaction transition of gp17. We find that although electrostatic interactions between charged residues contribute significantly to the overall free energy change of compaction, interactions mediated by the uncharged residues are equally if not more important. We identify five charged residues and six uncharged residues at the interface that play a dominant role in the compaction transition and also reveal salt bridging, van der Waals, and solvent hydrogen-bonding interactions mediated by these residues in stabilizing the compact form of gp17. The formation of a salt bridge between Glu309 and Arg494 is found to be particularly crucial, consistent with experiments showing complete abrogation in packaging upon Glu309Lys mutation. The computed contributions of several other residues are also found to correlate well with single-molecule measurements of impairments in DNA translocation activity caused by site-directed mutations.

A distinct switch in interactions of the histone H4 tail domain upon salt-dependent folding of nucleosome arrays

The core histone tail domains mediate inter-nucleosomal interactions that direct folding and condensation of nucleosome arrays into higher-order chromatin structures. The histone H4 tail domain facilitates inter-array interactions by contacting both the H2A/H2B acidic patch and DNA of neighboring nucleosomes. Likewise, H4 tail-H2A contacts stabilize array folding. However, whether the H4 tail domains stabilize array folding via inter-nucleosomal interactions with the DNA of neighboring nucleosomes remains unclear. We utilized defined oligonucleosome arrays containing a single specialized nucleosome with a photo-inducible cross-linker in the N terminus of the H4 tail to characterize these interactions. We observed that the H4 tail participates exclusively in intra-array interactions with DNA in unfolded arrays. These interactions are diminished during array folding, yet no inter-nucleosome, intra-array H4 tail-DNA contacts are observed in condensed chromatin. However, we document contacts between the N terminus of the H4 tail and H2A. Installation of acetylation mimics known to disrupt H4-H2A surface interactions did not increase observance of H4-DNA inter-nucleosomal interactions. These results suggest the multiple functions of the H4 tail require targeted distinct interactions within condensed chromatin.

Chromatin without the 30-nm fiber: constrained disorder instead of hierarchical folding

Several hierarchical levels of DNA packaging are believed to exist in chromatin, starting from a 10-nm chromatin fiber that is further packed into a 30-nm fiber. Transitions between the 30-nm and 10-nm fibers are thought to be essential for the control of chromatin transcriptional status. However, recent studies demonstrate that in the nuclei, DNA is packed in tightly associated 10-nm fibers that are not compacted into 30-nm fibers. Additionally, the accessibility of DNA in chromatin depends on the local mobility of nucleosomes rather than on decompaction of chromosome regions. These findings argue for reconsidering the hierarchical model of chromatin packaging and some of the basic definitions of chromatin. In particular, chromatin domains should be considered as three-dimensional objects, which may include genomic regions that do not necessarily constitute a continuous domain on the DNA chain.

Epigenetic molecular recognition: a biomolecular modeling perspective

The abnormal regulation of epigenetic protein families is associated with the onset and progression of various human diseases. However, epigenetic processes remain relatively obscure at the molecular level, thus preventing the rational design of chemical therapeutics. An array of robust computational and modeling approaches can complement experiments to shed light on the complex mechanisms of epigenetic molecular recognition and can guide medicinal chemists in designing selective and potent drug molecules. Herein we present a review of studies focused on epigenetic molecular recognition from a biomolecular modeling viewpoint. Although the known epigenetic targets are numerous, this review focuses on the more limited protein families on which computational modeling has been successfully applied. Therefore, we review three main topics: 1) histone deacetylases, 2) histone demethylases, and 3) histone tail dynamics. A brief review of the biological background and biomedical relevance is presented for each topic, followed by a detailed discussion of the computational studies and their relevance.

Nucleosome signalling; an evolving concept

The nucleosome core particle is the first stage of DNA packaging in virtually all eukaryotes. It both organises nuclear DNA and protects it from adventitious binding of transcription factors and the consequent deregulation of gene expression. Both properties are essential to allow the genome expansion characteristic of complex eukaryotes. The nucleosome is a flexible structure in vivo, allowing selective relaxation of its intrinsically inhibitory effects in response to external signals. Structural changes are brought about by dedicated remodelling enzymes and by posttranslational modifications of the core histones. Histone modifications occasionally alter nucleosome structure directly, but their more usual roles are to act as receptors on the nucleosome surface that are recognised by specific protein domains. The bound proteins, in turn, affect nucleosome structure and function. This strategy enormously expands the signalling capacity of the nucleosome and its ability to influence both the initiation and elongation stages of transcription. The enzymes responsible for placing and removing histone modifications, and the modification-binding proteins themselves, are ubiquitous, numerous and conserved amongst eukaryotes. Like the nucleosome, they date back to the earliest eukaryotes and may have played integral and essential roles in eukaryotic evolution. The present properties and epigenetic functions of the nucleosome reflect its evolutionary past and the selective pressures to which it has responded and can be better understood in this context. This article is part of a Special Issue entitled: Molecular mechanisms of histone modification function.

Establishing epigenetic domains via chromatin-bound histone modifiers

The eukaryotic nucleus harbors the DNA genome, which associates with histones and other chromosomal proteins into a complex referred to as chromatin. It provides an additional layer of so-called epigenetic information via histone modifications and DNA methylation on top of the DNA sequence that determines the cell's active gene expression program. The nucleus is devoid of internal organelles separated by membranes. Thus, free diffusive transport of proteins and RNA can occur throughout the space accessible for a given macromolecule. At the same time, chromatin is partitioned into different specialized structures such as nucleoli, chromosome territories, and heterochromatin domains that serve distinct functions. Here, we address the question of how the activity of chromatin-modifying enzymes is confined to chromatin subcompartments. We discuss mechanisms for establishing activity gradients of diffusive chromatin-modifying enzymes that could give rise to distinct chromatin domains within the cell nucleus. Interestingly, such gradients might directly result from immobilization of the enzymes on the flexible chromatin chain. Thus, locus-specific tethering of these enzymes to chromatin could have the potential to establish, maintain, or modulate epigenetic patterns of characteristic domain size.

Clipping of flexible tails of histones H3 and H4 affects the structure and dynamics of the nucleosome

Förster resonance energy transfer was used to monitor the dynamic conformations of mononucleosomes under different chromatin folding conditions to elucidate the role of the flexible N-terminal regions of H3 and H4 histones. The H3 tail was shown to partake in intranucleosomal interactions by restricting the DNA breathing motion and compacting the nucleosome. The H3 tail effects were mostly independent of the ionic strength and valency of the ions. The H4 tail was shown to not greatly affect the nucleosome conformation, but did slightly influence the relative population of the preferred conformation. The role of the H4 tail varied depending on the valency and ionic strength, suggesting that electrostatic forces play a primary role in H4 tail interactions. Interestingly, despite the H4 tail's lack of influence, when H3 and H4 tails were simultaneously clipped, a more dramatic effect was seen than when only H3 or H4 tails were clipped. The combinatorial effect of H3 and H4 tail truncation suggests a potential mechanism by which various combinations of histone tail modifications can be used to control accessibility of DNA-binding proteins to nucleosomal DNA.

The role of the nucleosome acidic patch in modulating higher order chromatin structure

Higher order folding of chromatin fibre is mediated by interactions of the histone H4 N-terminal tail domains with neighbouring nucleosomes. Mechanistically, the H4 tails of one nucleosome bind to the acidic patch region on the surface of adjacent nucleosomes, causing fibre compaction. The functionality of the chromatin fibre can be modified by proteins that interact with the nucleosome. The co-structures of five different proteins with the nucleosome (LANA, IL-33, RCC1, Sir3 and HMGN2) recently have been examined by experimental and computational studies. Interestingly, each of these proteins displays steric, ionic and hydrogen bond complementarity with the acidic patch, and therefore will compete with each other for binding to the nucleosome. We first review the molecular details of each interface, focusing on the key non-covalent interactions that stabilize the protein-acidic patch interactions. We then propose a model in which binding of proteins to the nucleosome disrupts interaction of the H4 tail domains with the acidic patch, preventing the intrinsic chromatin folding pathway and leading to assembly of alternative higher order chromatin structures with unique biological functions.

An advanced coarse-grained nucleosome core particle model for computer simulations of nucleosome-nucleosome interactions under varying ionic conditions

In the eukaryotic cell nucleus, DNA exists as chromatin, a compact but dynamic complex with histone proteins. The first level of DNA organization is the linear array of nucleosome core particles (NCPs). The NCP is a well-defined complex of 147 bp DNA with an octamer of histones. Interactions between NCPs are of paramount importance for higher levels of chromatin compaction. The polyelectrolyte nature of the NCP implies that nucleosome-nucleosome interactions must exhibit a great influence from both the ionic environment as well as the positively charged and highly flexible N-terminal histone tails, protruding out from the NCP. The large size of the system precludes a modelling analysis of chromatin at an all-atom level and calls for coarse-grained approximations. Here, a model of the NCP that include the globular histone core and the flexible histone tails described by one particle per each amino acid and taking into account their net charge is proposed. DNA wrapped around the histone core was approximated at the level of two base pairs represented by one bead (bases and sugar) plus four beads of charged phosphate groups. Computer simulations, using a Langevin thermostat, in a dielectric continuum with explicit monovalent (K(+)), divalent (Mg(2+)) or trivalent (Co(NH(3))(6) (3+)) cations were performed for systems with one or ten NCPs. Increase of the counterion charge results in a switch from repulsive NCP-NCP interaction in the presence of K(+), to partial aggregation with Mg(2+) and to strong mutual attraction of all 10 NCPs in the presence of CoHex(3+). The new model reproduced experimental results and the structure of the NCP-NCP contacts is in agreement with available data. Cation screening, ion-ion correlations and tail bridging contribute to the NCP-NCP attraction and the new NCP model accounts for these interactions.

Secondary structures of the core histone N-terminal tails: their role in regulating chromatin structure

The core histone N-terminal tails dissociate from their binding positions in nucleosomes at moderate salt concentrations, and appear unstructured in the crystal. This suggested that the tails contributed minimally to chromatin structure. However, in vitro studies have shown that the tails were involved in a range of intra- and inter-nucleosomal as well as inter-fibre contacts. The H4 tail, which is essential for chromatin compaction, was shown to contact an adjacent nucleosome in the crystal. Acetylation of H4K16 was shown to abolish the ability of a nucleosome array to fold into a 30 nm fibre. The application of secondary structure prediction software has suggested the presence of extended structured regions in the histone tails. Molecular Dynamics studies have further shown that sections of the H3 and H4 tails assumed α-helical and β-strand content that was enhanced by the presence of DNA, and that post-translational modifications of the tails had a major impact on these structures. Circular dichroism and NMR showed that the H3 and H4 tails exhibited significant α-helical content, that was increased by acetylation of the tail. There is thus strong evidence, both from biophysical and from computational approaches, that the core histones tails, particularly that of H3 and H4, are structured, and that these structures are influenced by post-translational modifications. This chapter reviews studies on the position, binding sites and secondary structures of the core histone tails, and discusses the possible role of the histone tail structures in the regulation of chromatin organization, and its impact on human disease.

Structural dynamics of nucleosomes at single-molecule resolution

The detailed mechanisms of how DNA that is assembled around a histone core can be accessed by DNA-binding proteins for transcription, replication, or repair, remain elusive nearly 40 years after Kornberg's nucleosome model was proposed. Uncovering the structural dynamics of nucleosomes is a crucial step in elucidating the mechanisms regulating genome accessibility. This requires the deconvolution of multiple structural states within an ensemble. Recent advances in single-molecule methods enable unprecedented efficiency in examining subpopulation dynamics. In this review, we summarize studies of nucleosome structure and dynamics from single-molecule approaches and how they advance our understanding of the mechanisms that govern DNA transactions.

Modelling chromatin structure and dynamics: status and prospects

The packaging of genomic DNA into chromatin in the eukaryotic cell nucleus demands extensive compaction. This requires attractive nucleosome-nucleosome interactions to overcome repulsion between the negatively charged DNA segments as well as other constraints. At the same time, DNA must be dynamically accessible to the cellular machinery that operates on it. Recent progress in the experimental characterisation of the higher order structure and dynamics of well-defined chromatin fibres has stimulated the attempts at theoretical description of chromatin and the nucleosome. Here we review the present status of chromatin modelling, with particular emphasis on coarse-grained computer simulation models, the role of electrostatic interactions, and discuss future perspectives in the field.

Toward convergence of experimental studies and theoretical modeling of the chromatin fiber

Understanding the structural organization of eukaryotic chromatin and its control of gene expression represents one of the most fundamental and open challenges in modern biology. Recent experimental advances have revealed important characteristics of chromatin in response to changes in external conditions and histone composition, such as the conformational complexity of linker DNA and histone tail domains upon compact folding of the fiber. In addition, modeling studies based on high-resolution nucleosome models have helped explain the conformational features of chromatin structural elements and their interactions in terms of chromatin fiber models. This minireview discusses recent progress and evidence supporting structural heterogeneity in chromatin fibers, reconciling apparently contradictory fiber models.

Finding a balance: how diverse dosage compensation strategies modify histone h4 to regulate transcription

Dosage compensation balances gene expression levels between the sex chromosomes and autosomes and sex-chromosome-linked gene expression levels between the sexes. Different dosage compensation strategies evolved in different lineages, but all involve changes in chromatin. This paper discusses our current understanding of how modifications of the histone H4 tail, particularly changes in levels of H4 lysine 16 acetylation and H4 lysine 20 methylation, can be used in different contexts to either modulate gene expression levels twofold or to completely inhibit transcription.

Structure and function of active chromatin and DNase I hypersensitive sites

Chromatin is by its very nature a repressive environment which restricts the recruitment of transcription factors and acts as a barrier to polymerases. Therefore the complex process of gene activation must operate at two levels. In the first instance, localized chromatin decondensation and nucleosome displacement is required to make DNA accessible. Second, sequence-specific transcription factors need to recruit chromatin modifiers and remodellers to create a chromatin environment that permits the passage of polymerases. In this review I will discuss the chromatin structural changes that occur at active gene loci and at regulatory elements that exist as DNase I hypersensitive sites.