Molecular basis of CENP-C association with the CENP-A nucleosome at yeast centromeres.

Histone CENP-A-containing nucleosomes play an important role in nucleating kinetochores at centromeres for chromosome segregation. However, the molecular mechanisms by which CENP-A nucleosomes engage with kinetochore proteins are not well understood. Here, we report the finding of a new function for the budding yeast Cse4/CENP-A histone-fold domain interacting with inner kinetochore protein Mif2/CENP-C. Strikingly, we also discovered that AT-rich centromere DNA has an important role for Mif2 recruitment. Mif2 contacts one side of the nucleosome dyad, engaging with both Cse4 residues and AT-rich nucleosomal DNA. Both interactions are directed by a contiguous DNA- and histone-binding domain (DHBD) harboring the conserved CENP-C motif, an AT hook, and RK clusters (clusters enriched for arginine-lysine residues). Human CENP-C has two related DHBDs that bind preferentially to DNA sequences of higher AT content. Our findings suggest that a DNA composition-based mechanism together with residues characteristic for the CENP-A histone variant contribute to the specification of centromere identity.

Structure and Dynamics of Compact Dinucleosomes: Analysis by Electron Microscopy and spFRET.

Formation of compact dinucleosomes (CODIs) occurs after collision between adjacent nucleosomes at active regulatory DNA regions. Although CODIs are likely dynamic structures, their structural heterogeneity and dynamics were not systematically addressed. Here, single-particle Förster resonance energy transfer (spFRET) and electron microscopy were employed to study the structure and dynamics of CODIs. spFRET microscopy in solution and in gel revealed considerable uncoiling of nucleosomal DNA from the histone octamer in a fraction of CODIs, suggesting that at least one of the nucleosomes is destabilized in the presence of the adjacent closely positioned nucleosome. Accordingly, electron microscopy analysis suggests that up to 30 bp of nucleosomal DNA are involved in transient uncoiling/recoiling on the octamer. The more open and dynamic nucleosome structure in CODIs cannot be stabilized by histone chaperone Spt6. The data suggest that proper internucleosomal spacing is an important determinant of chromatin stability and support the possibility that CODIs could be intermediates of chromatin disruption.

Interactions of Nucleosomes with Acidic Patch-Binding Peptides: A Combined Structural Bioinformatics, Molecular Modeling, Fluorescence Polarization, and Single-Molecule FRET Study.

In eukaryotic organisms, genomic DNA associates with histone proteins to form nucleosomes. Nucleosomes provide a basis for genome compaction, epigenetic markup, and mediate interactions of nuclear proteins with their target DNA loci. A negatively charged (acidic) patch located on the H2A-H2B histone dimer is a characteristic feature of the nucleosomal surface. The acidic patch is a common site in the attachment of various chromatin proteins, including viral ones. Acidic patch-binding peptides present perspective compounds that can be used to modulate chromatin functioning by disrupting interactions of nucleosomes with natural proteins or alternatively targeting artificial moieties to the nucleosomes, which may be beneficial for the development of new therapeutics. In this work, we used several computational and experimental techniques to improve our understanding of how peptides may bind to the acidic patch and what are the consequences of their binding. Through extensive analysis of the PDB database, histone sequence analysis, and molecular dynamic simulations, we elucidated common binding patterns and key interactions that stabilize peptide-nucleosome complexes. Through MD simulations and FRET measurements, we characterized changes in nucleosome dynamics conferred by peptide binding. Using fluorescence polarization and gel electrophoresis, we evaluated the affinity and specificity of the LANA1-22 peptide to DNA and nucleosomes. Taken together, our study provides new insights into the different patterns of intermolecular interactions that can be employed by natural and designed peptides to bind to nucleosomes, and the effects of peptide binding on nucleosome dynamics and stability.

Hydroxyl-radical footprinting combined with molecular modeling identifies unique features of DNA conformation and nucleosome positioning.

Nucleosomes are the most abundant protein-DNA complexes in eukaryotes that provide compaction of genomic DNA and are implicated in regulation of transcription, DNA replication and repair. The details of DNA positioning on the nucleosome and the DNA conformation can provide key regulatory signals. Hydroxyl-radical footprinting (HRF) of protein-DNA complexes is a chemical technique that probes nucleosome organization in solution with a high precision unattainable by other methods. In this work we propose an integrative modeling method for constructing high-resolution atomistic models of nucleosomes based on HRF experiments. Our method precisely identifies DNA positioning on nucleosome by combining HRF data for both DNA strands with the pseudo-symmetry constraints. We performed high-resolution HRF for Saccharomyces cerevisiae centromeric nucleosome of unknown structure and characterized it using our integrative modeling approach. Our model provides the basis for further understanding the cooperative engagement and interplay between Cse4p protein and the A-tracts important for centromere function.

Structural analysis of nucleosomal barrier to transcription.

Thousands of human and Drosophila genes are regulated at the level of transcript elongation and nucleosomes are likely targets for this regulation. However, the molecular mechanisms of formation of the nucleosomal barrier to transcribing RNA polymerase II (Pol II) and nucleosome survival during/after transcription remain unknown. Here we show that both DNA-histone interactions and Pol II backtracking contribute to formation of the barrier and that nucleosome survival during transcription likely occurs through allosterically stabilized histone-histone interactions. Structural analysis indicates that after Pol II encounters the barrier, the enzyme backtracks and nucleosomal DNA recoils on the octamer, locking Pol II in the arrested state. DNA is displaced from one of the H2A/H2B dimers that remains associated with the octamer. The data reveal the importance of intranucleosomal DNA-protein and protein-protein interactions during conformational changes in the nucleosome structure on transcription. Mechanisms of nucleosomal barrier formation and nucleosome survival during transcription are proposed.

Analysis of the mechanism of nucleosome survival during transcription.

Maintenance of nucleosomal structure in the cell nuclei is essential for cell viability, regulation of gene expression and normal aging. Our previous data identified a key intermediate (a small intranucleosomal DNA loop, Ø-loop) that is likely required for nucleosome survival during transcription by RNA polymerase II (Pol II) through chromatin, and suggested that strong nucleosomal pausing guarantees efficient nucleosome survival. To evaluate these predictions, we analysed transcription through a nucleosome by different, structurally related RNA polymerases and mutant yeast Pol II having different histone-interacting surfaces that presumably stabilize the Ø-loop. The height of the nucleosomal barrier to transcription and efficiency of nucleosome survival correlate with the net negative charges of the histone-interacting surfaces. Molecular modeling and analysis of Pol II-nucleosome intermediates by DNase I footprinting suggest that efficient Ø-loop formation and nucleosome survival are mediated by electrostatic interactions between the largest subunit of Pol II and core histones.

Structural Perspectives on the Evolutionary Expansion of Unique Protein-Protein Binding Sites.

Structures of protein complexes provide atomistic insights into protein interactions. Human proteins represent a quarter of all structures in the Protein Data Bank; however, available protein complexes cover less than 10% of the human proteome. Although it is theoretically possible to infer interactions in human proteins based on structures of homologous protein complexes, it is still unclear to what extent protein interactions and binding sites are conserved, and whether protein complexes from remotely related species can be used to infer interactions and binding sites. We considered biological units of protein complexes and clustered protein-protein binding sites into similarity groups based on their structure and sequence, which allowed us to identify unique binding sites. We showed that the growth rate of the number of unique binding sites in the Protein Data Bank was much slower than the growth rate of the number of structural complexes. Next, we investigated the evolutionary roots of unique binding sites and identified the major phyletic branches with the largest expansion in the number of novel binding sites. We found that many binding sites could be traced to the universal common ancestor of all cellular organisms, whereas relatively few binding sites emerged at the major evolutionary branching points. We analyzed the physicochemical properties of unique binding sites and found that the most ancient sites were the largest in size, involved many salt bridges, and were the most compact and least planar. In contrast, binding sites that appeared more recently in the evolution of eukaryotes were characterized by a larger fraction of polar and aromatic residues, and were less compact and more planar, possibly due to their more transient nature and roles in signaling processes.

Voltage-gated ion channel modulation by lipids: insights from molecular dynamics simulations.

Cells commonly use lipids to modulate the function of ion channels. The lipid content influences the amplitude of the ionic current and changes the probability of voltage-gated ion channels being in the active or in the resting states. Experimental findings inferred from a variety of techniques and molecular dynamics studies have revealed a direct interaction between the lipid headgroups and the ion channel residues, suggesting an influence on the ion channel function. On the other hand the alteration of the lipids may in principle modify the overall electrostatic environment of the channel, and hence the transmembrane potential, leading to an indirect modulation, i.e. a global effect. Here we have investigated the structural and dynamical properties of the voltage-gated potassium channel Kv1.2 embedded in bilayers with modified upper or lower leaflet compositions corresponding to realistic biological scenarios: the first relates to the effects of sphingomyelinase, an enzyme that modifies the composition of lipids of the outer membrane leaflets, and the second to the effect of the presence of a small fraction of PIP2, a highly negatively charged lipid known to modulate voltage-gated channel function. Our molecular dynamics simulations do not enable to exclude the global effect mechanism in the former case. For the latter, however, it is shown that local interactions between the ion channel and the lipid headgroups are key-elements of the modulation.

Direct prediction of residual dipolar couplings of small molecules in a stretched gel by stochastic molecular dynamics simulations.

Residual dipolar couplings are highly useful NMR parameters for calculating and refining molecular structures, dynamics, and interactions. For some applications, however, it is inevitable that the preferred orientation of a molecule in an alignment medium is calculated a priori. Several methods have been developed to predict molecular orientations and residual dipolar couplings. Being beneficial for macromolecules and selected small-molecule applications, such approaches lack sufficient accuracy for a large number of organic compounds for which the fine structure and eventually the flexibility of all involved molecules have to be considered or are limited to specific, well-studied liquid crystals. We introduce a simplified model for detailed all-atom molecular dynamics calculations with a polymer strand lined up along the principal axis as a new approach to simulate the preferred orientation of small to medium-sized solutes in polymer-based, gel-type alignment media. As is shown by a first example of strychnine in a polystyrene/CDCl3 gel, the simulations potentially enable the accurate prediction of residual dipolar couplings taking into account structural details and dynamic averaging effects of both the polymer and the solute.

Nucleosomes and their complexes in the cryoEM era: Trends and limitations.

Twenty-five years have passed since the appearance of the first atomistic model of the nucleosome structure, and since then the number of new structures has gradually increased. With the advent of cryo-microscopy, the rate of accumulation of models has increased significantly. New structures are emerging with different histone variants and a variety of proteins that bind to nucleosomes. At the moment, there are more than four hundred structures containing nucleosomes in the Protein Data Bank. Many of these structures represent similar complexes, others differ in composition, conformation and quality. In this perspective, we investigate the diversity of known nucleosome structures, analyze data and model quality, variations in histone/DNA content of nucleosomes and spectrum of their interactors. We outline those parts of the nucleosome "structurome" that are already explored and those awaiting further exploration.

H2A-H2B Histone Dimer Plasticity and Its Functional Implications.

The protein core of the nucleosome is composed of an H3-H4 histone tetramer and two H2A-H2B histone dimers. The tetramer organizes the central 60 DNA bp, while H2A-H2B dimers lock the flanking DNA segments. Being positioned at the sides of the nucleosome, H2A-H2B dimers stabilize the overall structure of the nucleosome and modulate its dynamics, such as DNA unwrapping, sliding, etc. Such modulation at the epigenetic level is achieved through post-translational modifications and the incorporation of histone variants. However, the detailed connection between the sequence of H2A-H2B histones and their structure, dynamics and implications for nucleosome functioning remains elusive. In this work, we present a detailed study of H2A-H2B dimer dynamics in the free form and in the context of nucleosomes via atomistic molecular dynamics simulations (based on X. laevis histones). We supplement simulation results by comparative analysis of information in the structural databases. Particularly, we describe a major dynamical mode corresponding to the bending movement of the longest H2A and H2B α-helices. This overall bending dynamics of the H2A-H2B dimer were found to be modulated by its interactions with DNA, H3-H4 tetramer, the presence of DNA twist-defects with nucleosomal DNA and the amino acid sequence of histones. Taken together, our results shed new light on the dynamical mechanisms of nucleosome functioning, such as nucleosome sliding, DNA-unwrapping and their epigenetic modulation.

Electron microscopy analysis of ATP-independent nucleosome unfolding by FACT.

FACT is a histone chaperone that participates in nucleosome removal and reassembly during transcription and replication. We used electron microscopy to study FACT, FACT:Nhp6 and FACT:Nhp6:nucleosome complexes, and found that all complexes adopt broad ranges of configurations, indicating high flexibility. We found unexpectedly that the DNA binding protein Nhp6 also binds to the C-terminal tails of FACT subunits, inducing more open geometries of FACT even in the absence of nucleosomes. Nhp6 therefore supports nucleosome unfolding by altering both the structure of FACT and the properties of nucleosomes. Complexes formed with FACT, Nhp6, and nucleosomes also produced a broad range of structures, revealing a large number of potential intermediates along a proposed unfolding pathway. The data suggest that Nhp6 has multiple roles before and during nucleosome unfolding by FACT, and that the process proceeds through a series of energetically similar intermediate structures, ultimately leading to an extensively unfolded form.

A standardized nomenclature for mammalian histone genes.

Histones have a long history of research in a wide range of species, leaving a legacy of complex nomenclature in the literature. Community-led discussions at the EMBO Workshop on Histone Variants in 2011 resulted in agreement amongst experts on a revised systematic protein nomenclature for histones, which is based on a combination of phylogenetic classification and historical symbol usage. Human and mouse histone gene symbols previously followed a genome-centric system that was not applicable across all vertebrate species and did not reflect the systematic histone protein nomenclature. This prompted a collaboration between histone experts, the Human Genome Organization (HUGO) Gene Nomenclature Committee (HGNC) and Mouse Genomic Nomenclature Committee (MGNC) to revise human and mouse histone gene nomenclature aiming, where possible, to follow the new protein nomenclature whilst conforming to the guidelines for vertebrate gene naming. The updated nomenclature has also been applied to orthologous histone genes in chimpanzee, rhesus macaque, dog, cat, pig, horse and cattle, and can serve as a framework for naming other vertebrate histone genes in the future.

N-Terminal Tails of Histones H2A and H2B Differentially Affect Transcription by RNA Polymerase II In Vitro.

Histone N-terminal tails and their post-translational modifications affect various biological processes, often in a context-specific manner; the underlying mechanisms are poorly studied. Here, the role of individual N-terminal tails of histones H2A/H2B during transcription through chromatin was analyzed in vitro. spFRET data suggest that the tail of histone H2B (but not of histone H2A) affects nucleosome stability. Accordingly, deletion of the H2B tail (amino acids 1-31, but not 1-26) causes a partial relief of the nucleosomal barrier to transcribing RNA polymerase II (Pol II), likely facilitating uncoiling of DNA from the histone octamer during transcription. Taken together, the data suggest that residues 27-31 of histone H2B stabilize DNA-histone interactions at the DNA region localized ~25 bp in the nucleosome and thus interfere with Pol II progression through the region localized 11-15 bp in the nucleosome. This function of histone H2B requires the presence of the histone H2A N-tail that mediates formation of nucleosome-nucleosome dimers; however, nucleosome dimerization per se plays only a minimal role during transcription. Histone chaperone FACT facilitates transcription through all analyzed nucleosome variants, suggesting that H2A/H2B tails minimally interact with FACT during transcription; therefore, an alternative FACT-interacting domain(s) is likely involved in this process.

Structure of an Intranucleosomal DNA Loop That Senses DNA Damage during Transcription.

Transcription through chromatin by RNA polymerase II (Pol II) is accompanied by the formation of small intranucleosomal DNA loops containing the enzyme (i-loops) that are involved in survival of core histones on the DNA and arrest of Pol II during the transcription of damaged DNA. However, the structures of i-loops have not been determined. Here, the structures of the intermediates formed during transcription through a nucleosome containing intact or damaged DNA were studied using biochemical approaches and electron microscopy. After RNA polymerase reaches position +24 from the nucleosomal boundary, the enzyme can backtrack to position +20, where DNA behind the enzyme recoils on the surface of the histone octamer, forming an i-loop that locks Pol II in the arrested state. Since the i-loop is formed more efficiently in the presence of SSBs positioned behind the transcribing enzyme, the loop could play a role in the transcription-coupled repair of DNA damage hidden in the chromatin structure.

Molecular Mechanisms of Oncogenesis through the Lens of Nucleosomes and Histones.

At the cellular level, cancer is the disease of both the genome and the epigenome, and the interplay between genetic mutations and epigenetic states may occur at the level of elementary chromatin units, the nucleosomes. They are formed by a segment of DNA wrapped around an octamer of histone proteins. In this review, we survey various mechanisms of cancer etiology and progression mediated by histones and nucleosomes. In particular, we discuss the effects of mutations in histones, changes in their expression and slicing on epigenetic dysregulation and carcinogenesis. The links between cancer phenotypes and differential expression of histone variants and isoforms are summarized. Finally, we discourse the geometric and steric effects of DNA compaction in nucleosomes on DNA mutation rate, interactions with transcription factors, including pioneer transcription factors, and prospects of cancer cells' genome and epigenome editing.

Histone dynamics mediate DNA unwrapping and sliding in nucleosomes.

Nucleosomes are elementary building blocks of chromatin in eukaryotes. They tightly wrap ∼147 DNA base pairs around an octamer of histone proteins. How nucleosome structural dynamics affect genome functioning is not completely clear. Here we report all-atom molecular dynamics simulations of nucleosome core particles at a timescale of 15 microseconds. At this timescale, functional modes of nucleosome dynamics such as spontaneous nucleosomal DNA breathing, unwrapping, twisting, and sliding were observed. We identified atomistic mechanisms of these processes by analyzing the accompanying structural rearrangements of the histone octamer and histone-DNA contacts. Octamer dynamics and plasticity were found to enable DNA unwrapping and sliding. Through multi-scale modeling, we showed that nucleosomal DNA dynamics contribute to significant conformational variability of the chromatin fiber at the supranucleosomal level. Our study further supports mechanistic coupling between fine details of histone dynamics and chromatin functioning, provides a framework for understanding the effects of various chromatin modifications.

Free energy profiles of amino acid side chain analogs near water-vapor interface obtained via MD simulations.

Solvation of 13 neutral amino acid side chain analogs at water-vapor interface was studied by computing high precision free energy profiles of the molecules across the interface using molecular dynamics (MD) simulations. The SPC water model (Berendsen, H. J. C., Postma, J. P. M., van Gunsteren, W. F., Hermans, P. A. K. J., Dixon, R., Cornell, W., Fox, T., Chipot, C., Pohorille, A. In: Wilkinson, A., Weiner, P. and van Gunsteren, W. F. editors. Intermolecular Forces, 1981, 3, 331) and OPLS-AA (Jorgensen et al., J Am Chem Soc, 1996, 118, 11225) potential parameter sets were used. A rigorous approach for the computation of high precision free energy profiles at water-vapor interface using constraint force technique is implemented. Methodology of obtaining high precision potential of mean force (PMF) profiles free of simulation artifacts in MD simulations is outlined and discussed. The accuracy of the calculations is examined by comparing the hydration free energies of studied solutes obtained from PMF calculations and separately using Bennett acceptance ratio technique by decoupling solvent-solute interactions in the bulk. All molecules exhibit a free energy minimum at interface. No significant desolvation barrier is observed for any of the studied species. Adsorption energies for studied molecules at water-vapor interface are estimated and compared with experimental observations. We find that for modeled neutral compounds pronounced surface influence on solute solvation vanishes already at 6-7 A behind the water surface as the solvation free energy approaches the bulk value. The possibility of force field refinement using adsorption free energies is outlined.

Solvent accessible surface area of amino acid residues in globular proteins: correlation of apparent transfer free energies with experimental hydrophobicity scales.

It is known that the distribution of amino acid residues in globular proteins between surface and interior is in certain correlation with various experimental scales based on partitioning of amino acids or their analogs between water and organic solvents. These scales are often used in various quantitative structure-activity relationship (QSAR) studies as well as for evaluation of stability of proteins. In this work we have analyzed the distribution of residues based on their solvent accessible surface area in more than 8000 protein structures. Using extensive statistical sampling, we have computed residue apparent free energies of transfer between protein interior and surface applying various criteria for classifying residues as exposed or buried. The correlation of these statistical energies with several experimental hydrophobicity scales is discussed. We propose three types of statistical apparent transfer free energy scales and show that each of these scales is in better correlation with one of the experimental hydrophobicity scales (water/vapor, water/cyclohexane, and water/octanol transfer scales). The data are interpreted through the application of theoretical considerations by Finkelstein et al. (Protein Struct. Funct. Genet. 1995, 23, 142) based on random energy model of heteropolymer globules. The deviation of apparent transfer free energies from experimental scales is discussed and analyzed. The variations of amino acid distribution in proteins with the size of protein structure is discussed and the final protein set is chosen to minimize these variations.

Linking chromatin composition and structural dynamics at the nucleosome level.

Nucleosomes are fundamental units of chromatin compaction, which organize ∼200 DNA base pairs using an octamer of histone proteins. Their ubiquitous presence in the cell nucleus since the first eukaryotes compelled the chromatin machinery to coevolve and learn how to exploit various modes of nucleosome dynamics and sense differences in nucleosome composition. Alterations to histone or DNA sequences, post-translational modifications (PTM) of histones, recruitment of chromatin proteins modulate nucleosome dyn amics and provide epigenetic regulation to the DNA processing pathways (transcription, replication, repair, etc.). Our understanding of this complex interplay between nucleosome composition, dynamics and functioning is constantly evolving through new insights and discoveries. In this review, we highlight recent contributions to the field while attempting to organize them in a unified framework.