Radial density distribution of chromatin: evidence that chromatin fibers have solid centers

Implicit Solvent Model for Million-Atom Atomistic Simulations: Insights into the Organization of 30-nm Chromatin Fiber

Molecular dynamics (MD) simulations based on the implicit solvent generalized Born (GB) models can provide significant computational advantages over the traditional explicit solvent simulations. However, the standard GB becomes prohibitively expensive for all-atom simulations of large structures; the model scales poorly, ∼n2, with the number of solute atoms. Here we combine our recently developed optimal point charge approximation (OPCA) with the hierarchical charge partitioning (HCP) approximation to present an ∼n log n multiscale, yet fully atomistic, GB model (GB-HCPO). The HCP approximation exploits the natural organization of biomolecules (atoms, groups, chains, and complexes) to partition the structure into multiple hierarchical levels of components. OPCA approximates the charge distribution for each of these components by a small number of point charges so that the low order multipole moments of these components are optimally reproduced. The approximate charges are then used for computing electrostatic interactions with distant components, while the full set of atomic charges are used for nearby components. We show that GB-HCPO can deliver up to 2 orders of magnitude speedup compared to the standard GB, with minimal impact on its accuracy. For large structures, GB-HCPO can approach the same nominal speed, as in nanoseconds per day, as the highly optimized explicit-solvent simulation based on particle mesh Ewald (PME). The increase in the nominal simulation speed, relative to the standard GB, coupled with substantially faster sampling of conformational space, relative to the explicit solvent, makes GB-HCPO a suitable candidate for MD simulation of large atomistic systems in implicit solvent. As a practical demonstration, we use GB-HCPO simulation to refine a ∼1.16 million atom structure of 30 nm chromatin fiber (40 nucleosomes). The refined structure suggests important details about spatial organization of the linker DNA and the histone tails in the fiber: (1) the linker DNA fills the core region, allowing the H3 histone tails to interact with the linker DNA, which is consistent with experiment; (2) H3 and H4 tails are found mostly in the core of the structure, closer to the helical axis of the fiber, while H2A and H2B are mostly solvent exposed. Potential functional consequences of these findings are discussed. GB-HCPO is implemented in the open source MD software NAB in Amber 2016.

Chromatin structure outside and inside the nucleus

The structure of the 30-nm chromatin fiber has provided, over the years, an important reference in chromatin studies. Originally derived from electron microscopic studies of soluble chromatin fibers released by restriction digestion, the gross structural features of such fragments have been supported by biophysical methods such as low angle X-ray and neutron scattering, sedimentation, light scattering, and electric dichroism. Electron microscopy and sedimentation velocity measurements demonstrated that reconstituted chromatin fibers, prepared from repeating arrays of high affinity nucleosome positioning sequences, retain the same overall features as observed for native chromatin fibers. It had been suggested that the 30 nm fiber might be the form assumed in vivo by transcriptionally silent chromatin, but individual gene or genome-wide studies of chromatin released from nuclei do not reveal any such simple correlation. Furthermore, even though the 30 nm fiber has been thought to represent an intermediate in the hierarchical folding of DNA into chromosomes, most analyses of chromatin folding within the nucleus do not detect any regular extended compact structures. However, there are important exceptions in chicken erythroid cell nuclei as well as in transcribed regions that form extended loops. Localized domains within the nucleus, either at the surface of chromosome domains or constrained as a specialized kind of constitutive heterochromatin by specific DNA binding proteins, may adopt 30 nm fiber-like structures.

Structure and organization of chromatin fiber in the nucleus

Eukaryotic genomes are organized hierarchically into chromatin structures by histones. Despite extensive research for over 30 years, not only the fundamental structure of the 30-nm chromatin fiber is being debated, but the actual existence of such fiber remains hotly contested. In this review, we focus on the most recent progress in elucidating the structure of the 30-nm fiber upon in vitro reconstitution, and its possible organization inside the nucleus. In addition, we discuss the roles of linker histone H1 as well as the importance of specific nucleosome-nucleosome interactions in the formation of the 30-nm fiber. Finally, we discuss the involvement of structural variations and epigenetic mechanisms available for the regulation of this chromatin form.

Cryo-EM study of the chromatin fiber reveals a double helix twisted by tetranucleosomal units

The hierarchical packaging of eukaryotic chromatin plays a central role in transcriptional regulation and other DNA-related biological processes. Here, we report the 11-angstrom-resolution cryogenic electron microscopy (cryo-EM) structures of 30-nanometer chromatin fibers reconstituted in the presence of linker histone H1 and with different nucleosome repeat lengths. The structures show a histone H1-dependent left-handed twist of the repeating tetranucleosomal structural units, within which the four nucleosomes zigzag back and forth with a straight linker DNA. The asymmetric binding and the location of histone H1 in chromatin play a role in the formation of the 30-nanometer fiber. Our results provide mechanistic insights into how nucleosomes compact into higher-order chromatin fibers.

Modeling studies of chromatin fiber structure as a function of DNA linker length

Chromatin fibers encountered in various species and tissues are characterized by different nucleosome repeat lengths (NRLs) of the linker DNA connecting the nucleosomes. While single cellular organisms and rapidly growing cells with high protein production have short NRL ranging from 160 to 189 bp, mature cells usually have longer NRLs ranging between 190 and 220 bp. Recently, various experimental studies have examined the effect of NRL on the internal organization of chromatin fiber. Here, we investigate by mesoscale modeling of oligonucleosomes the folding patterns for different NRL, with and without linker histone (LH), under typical monovalent salt conditions using both one-start solenoid and two-start zigzag starting configurations. We find that short to medium NRL chromatin fibers (173 to 209 bp) with LH condense into zigzag structures and that solenoid-like features are viable only for longer NRLs (226 bp). We suggest that medium NRLs are more advantageous for packing and various levels of chromatin compaction throughout the cell cycle than their shortest and longest brethren; the former (short NRLs) fold into narrow fibers, while the latter (long NRLs) arrays do not easily lead to high packing ratios due to possible linker DNA bending. Moreover, we show that the LH has a small effect on the condensation of short-NRL arrays but has an important condensation effect on medium-NRL arrays, which have linker lengths similar to the LH lengths. Finally, we suggest that the medium-NRL species, with densely packed fiber arrangements, may be advantageous for epigenetic control because their histone tail modifications can have a greater effect compared to other fibers due to their more extensive nucleosome interaction network.

Maintenance of a functional higher order chromatin structure: The role of the nuclear matrix in normal and disease states

The ordered packaging of DNA within the nucleus of somatic cells reflects a dynamic supportive structure that facilitates stable transcription interrupted by intermittent cycles of extreme condensation. This dynamic mode of packing and unpacking chromatin is intimately linked to the ability of the genome to specifically complex with both histones and non-histone proteins. Understanding the underlying mechanism that governs the formation of higher order chromatin structures is a key to understanding how local architecture modulates transcription. In part, the formation of these structures appears to be regulated through genomic looping that is dynamically mediated by attachment to the nuclear scaffold/matrix at S/MARs, i.e., Scaffold/Matrix Attachment Regions. Although the mechanism guiding the formation and use of these higher-ordered structures remains unknown, S/MARs continue to reveal a multitude of roles in development and the pathogenesis of disease.

A variable topology for the 30-nm chromatin fibre

The structure of the 30-nm chromatin fibre is an important determinant of the regulation of eukaryotic transcription. A fundamental issue is whether the stacking of nucleosomes in this fibre is organized as a one-start or two-start helix. We argue that all recent experimental data are compatible with a two-start helix and propose that the topology of the fibre, but not the mode of stacking the nucleosomes, is dependent on the length of the linker DNA. This arrangement conserves nucleosome stacking and thus the external morphology of the fibre, and also ensures that the fibre adopts the highest available packing density.

EM measurements define the dimensions of the "30-nm" chromatin fiber: evidence for a compact, interdigitated structure

Chromatin structure plays a fundamental role in the regulation of nuclear processes such as DNA transcription, replication, recombination, and repair. Despite considerable efforts during three decades, the structure of the 30-nm chromatin fiber remains controversial. To define fiber dimensions accurately, we have produced very long and regularly folded 30-nm fibers from in vitro reconstituted nucleosome arrays containing the linker histone and with increasing nucleosome repeat lengths (10 to 70 bp of linker DNA). EM measurements show that the dimensions of these fully folded fibers do not increase linearly with increasing linker length, a finding that is inconsistent with two-start helix models. Instead, we find that there are two distinct classes of fiber structure, both with unexpectedly high nucleosome density: arrays with 10 to 40 bp of linker DNA all produce fibers with a diameter of 33 nm and 11 nucleosomes per 11 nm, whereas arrays with 50 to 70 bp of linker DNA all produce 44-nm-wide fibers with 15 nucleosomes per 11 nm. Using the physical constraints imposed by these measurements, we have built a model in which tight nucleosome packing is achieved through the interdigitation of nucleosomes from adjacent helical gyres. Importantly, the model closely matches raw image projections of folded chromatin arrays recorded in the solution state by using electron cryo-microscopy.

DNase I digestion reveals alternating asymmetrical protection of the nucleosome by the higher order chromatin structure

DNase I was used to probe the higher order chromatin structure in whole nuclei. The digestion profiles obtained were the result of single-stranded cuts and were independent of pH, type of divalent ion and chromatin repeat length. Furthermore, the protection from digestion of the DNA at the entry/exit points on the nucleosome was found to be caused not by the H1/H5 histone tails, but by the compact structure that these proteins support. In order to resolve symmetry ambiguities, DNase I digestion fragments over several nucleosome repeat lengths were analysed quantitatively and compared with computer simulations using combinations of the experimentally obtained rate constants (some of which were converted to 0 to simulate steric protection from DNase I digestion). A clear picture of precisely defined, alternating, asymmetrically protected nucleosomes emerged. The linker DNA is inside the fibre, while the nucleosomes are positioned above and below a helical path and/or with alternating orientation towards the dyad axis. The dinucleosomal modulation of the digestion patterns comes from alternate protection of cutting sites inside the nucleosome and not from alternating exposure to the enzyme of the linker DNA.

Quantitative analysis of DNase I digestion patterns of oligo- and polynucleosomes

DNase I digestion of unlabelled chromatin has been used to determine the orientation of nucleosomes in the higher-order structure of chromatin. DNA digestion patterns were analysed quantitatively and compared with simulation curves that were generated with the experimentally obtained rate constants for cutting inside nucleosomes and the half-height widths of the bands. The rate constants for cutting at the linker DNA were varied to fit the experimental curves. By comparing the digestion profiles of polynucleosomes, oligonucleosomes and H1/H5-depleted oligonucleosomes we have been able to distinguish between protection caused by H1/H5 histones only and protection caused by the higher-order structure. By the nature of this method the area protected by H1/H5 histones appears symmetrical around the centre of the nucleosomal DNA on the dyad axis (position S[0]), but it can also be interpreted as a superposition of two separate protected areas that are symmetrical around position S[0]. Protection by the higher-order structure shows that the nucleosomes are oriented with their S[0] positions inside the 30 nm fibre and that there is a minimum number of nucleosomes required for this structure to be formed. We have also found that the linker DNA is cut in a continuous (non-periodical) way and that there are considerable amounts of nucleosomes at discrete distances of multiples of ten nucleotides (10a nt) as well as stretches of nucleosomes positioned either at random or at distances of (10a + 5)nt.

Internal structure of the 30 nm chromatin fiber

In the presence of 1.7 mM Mg2+, the diameter of the circular structures produced by small chromatin fragments isolated from chicken erythrocytes remains essentially unchanged when the number of nucleosomes in these fragments increases from 10 to 36. In contrast, the results obtained in unidirectional shadowing experiments show that under the same conditions the height of the chromatin fragments increases with the number of nucleosomes. These observations indicate that the electron microscope images studied in this work correspond to a top view of small chromatin fragments. Rotary-shadowed chromatin fragments show three parts: (a) a contour with a heavy deposition of platinum; (b) an annular zone between the central region and the periphery; and (c) a central hole. The heterogeneous ring generated by the deposition of platinum in the periphery suggests that nucleosomes form a one-start helix (5-7 nucleosomes per turn) that apparently can be left- or right-handed. The annular region (thickness of about 11 nm) shows spokes probably due to flat faces and core DNA of radially oriented nucleosomes. The central hole (8-12 nm) is clearly seen in many images but it is not empty because some deformed fragments show coated material (probably linker DNA) that protrudes from this central depression. We have observed that these structural elements directly detected in short chromatin fragments are also present in long chromatin fibers. This allows us to conclude that these elements are basic structural components of the 30 nm chromatin fiber.

Poly(ADP-ribosyl)ation, DNA strand breaks, chromatin and cancer

Studies on poly(ADP-ribosyl)ation related to chromatin structure and to nuclear functions such as repair, gene expression and replication are reviewed. Poly(ADP-ribosyl)ation might be involved in regulating the activity of nuclear enzymes involved in the metabolism of DNA strand breaks such as ligase II and topoisomerases I and II. In addition, it modifies nuclear proteins participating in gene expression including HMG non-histones, large T antigen, acetylated histone H4 and nuclear matrix proteins. It is speculated that poly(ADP-ribose) can induce free DNA domains by removing histones from specific nucleosomes whose DNA has been damaged. This process is proposed to require specific proteins recognizing lesions on DNA that ultimately attach the damaged site on the nuclear matrix where the repair enzymes are located. The role of poly(ADP-ribosyl)ation in carcinogenesis arises from that inhibitors of this modification potentiate the cytotoxicity of DNA-damaging drugs used in cancer chemotherapy and either enhance or inhibit tumor growth.

Quantitative energy-filtered electron microscopy of biological molecules in ice

The theoretical and experimental bases for quantitative electron microscopy of frozen-hydrated specimens are described, with special considerations of energy filtration to improve the images. The elastic and inelastic scattering from molecules in vacuum and in ice are calculated, and simple methods to approximate scattering are introduced. Multiple scattering calculations are used to describe the scattering from vitreous ice and to predict the characteristics of images of frozen-hydrated molecules as a function of ice thickness and accelerating voltage. Energy filtration is predicted to improve image contrast and signal-to-noise ratio. Experimental values for the inelastic scattering of ice, the energy spectrum of thick ice, and the contrast of biological specimens are determined. The principles of compensation for the contrast transfer function are presented. Tobacco mosaic virus is used to quantify the accuracy of interpreting image intensities to derive the absolute mass, mass per unit length, and internal mass densities of biological molecules. It is shown that compensation for the contrast transfer function is necessary and sufficient to convert the images into accurate representations of molecular density. At a resolution of 2 nm, the radial density reconstructions of tobacco mosaic virus are in quantitative agreement with the atomic model derived from X-ray results.

DNA-bridging by a palindromic alpha-helix

The nucleosomal DNA repeat of 240 base pairs in the chromatin structure of sea urchin sperm is exceptionally long and is accompanied by the presence of a histone H1 molecule larger than is usual in most species of chromatin. I propose how these two features are correlated and how they fit into the solenoidal model for the 300-A-diameter fiber of chromatin. Comparison of the sequence of spermatogenous H1 with other H1 sequences reveals an insert of 55 amino acid residues (residues 122-176). A 37-residue sequence in the insert (residues 140-176) has a palindromic character. I propose that each half of the palindromic sequence constitutes an alpha-helical DNA-binding unit and that the continuous alpha-helix made up of the two halves, by virtue of its palindromic nature, stabilizes the formation of an extra superhelical turn by the long linker DNA between two nucleosome cores. The N-terminal-C-terminal "polarity" of each alpha-helical section of half the palindromic sequence indicates how the arginine/lysine-rich DNA-binding surface of the alpha-helical section is used. The polarity of the H1 insertion sequence supports the so-called "reverse-loop" model or a "figure-eight" model for the path of the DNA within the solenoid structure; i.e., the linker DNA forms a right-handed superhelical turn toward the center of the solenoid structure. This use of a pair of a palindromically related alpha-helical sections has a similarity with the "scissors-grip" model for the interaction of the leucine-zipper proteins with DNA.

Quantitation of molecular densities by cryo-electron microscopy. Determination of the radial density distribution of tobacco mosaic virus

We have determined the absolute mass and radial scattering density distribution of tobacco mosaic virus in the frozen-hydrated state by energy-filtered low-dose bright-field transmission electron microscopy. The absolute magnitude of electron scattering from tobacco mosaic virus in 150 nm of ice was within 3.0% of that predicted, with inelastic scattering accounting for approximately 80% of the scattering contrast. In order to test the accuracy of the radial reconstruction, a computer model of tobacco mosaic virus was built from the atomic co-ordinates assuming uniform solvent density. The validity of the model was confirmed by comparison of X-ray scattering and predictions of the model (R factor = 0.05). First-order corrections for the microscope contrast transfer function were necessary and sufficient for conversion of the cryo-electron microscopy images into accurate representations of the mass density. At 1.9 nm resolution the compensated reconstruction and model had density peaks of similar magnitude at 2.4, 4.2, 6.0 and 7.8 nm radius and a central hole of 2 nm radius. Equatorial Fourier transforms of the corrected electron images were in excellent agreement with predictions of the model (R factor = 0.12). Thus, the uniform solvent approximation was adequate at 1.9 nm resolution to describe quantitatively X-ray scattering in liquid water and electron imaging in vitreous ice. This is the first demonstration that cryo-electron microscopy images can be used to quantitate the absolute mass, mass per unit length and internal density distributions of proteins and nucleic acids.

Small angle x-ray scattering of chromatin. Radius and mass per unit length depend on linker length

Analyses of low angle x-ray scattering from chromatin, isolated by identical procedures but from different species, indicate that fiber diameter and number of nucleosomes per unit length increase with the amount of nucleosome linker DNA. Experiments were conducted at physiological ionic strength to obtain parameters reflecting the structure most likely present in living cells. Guinier analyses were performed on scattering from solutions of soluble chromatin from Necturus maculosus erythrocytes (linker length 48 bp), chicken erythrocytes (linker length 64 bp), and Thyone briareus sperm (linker length 87 bp). The results were extrapolated to infinite dilution to eliminate interparticle contributions to the scattering. Cross-sectional radii of gyration were found to be 10.9 +/- 0.5, 12.1 +/- 0.4, and 15.9 +/- 0.5 nm for Necturus, chicken, and Thyone chromatin, respectively, which are consistent with fiber diameters of 30.8, 34.2, and 45.0 nm. Mass per unit lengths were found to be 6.9 +/- 0.5, 8.3 +/- 0.6, and 11.8 +/- 1.4 nucleosomes per 10 nm for Necturus, chicken, and Thyone chromatin, respectively. The geometrical consequences of the experimental mass per unit lengths and radii of gyration are consistent with a conserved interaction among nucleosomes. Cross-linking agents were found to have little effect on fiber external geometry, but significant effect on internal structure. The absolute values of fiber diameter and mass per unit length, and their dependencies upon linker length agree with the predictions of the double-helical crossed-linker model. A compilation of all published x-ray scattering data from the last decade indicates that the relationship between chromatin structure and linker length is consistent with data obtained by other investigators.

Direct detection of linker DNA bending in defined-length oligomers of chromatin

Linker DNA, which connects between nucleosomes in chromatin, is short and, therefore, may be essentially straight and inflexible. We have carried out hydrodynamic and electron microscopic studies of dinucleosomes--fragments of chromatin containing just two nucleosomes--to test the ability of linker DNA to bend. We find that ionic conditions that stabilize the folding of long chromatin cause linker DNA in dinucleosomes to bend, bringing the two nucleosomes into contact. The results uphold a key prediction of the solenoid model of chromosome folding and suggest a mechanism by which proteins that are separated along the DNA can interact by direct contact.

The diameters of frozen-hydrated chromatin fibers increase with DNA linker length: evidence in support of variable diameter models for chromatin

The diameters of chromatin fibers from Thyone briareus (sea cucumber) sperm (DNA linker length, n = 87 bp) and Necturus maculosus (mudpuppy) erythrocytes (n = 48 bp) were investigated. Soluble fibers were frozen into vitrified aqueous solutions of physiological ionic strength (124 mM), imaged by cryo-EM, and measured interactively using quantitative computer image-processing techniques. Frozen-hydrated Thyone and Necturus fibers had significantly different mean diameters of 43.5 nm (SD = 4.2 nm; SEM = 0.61 nm) and 32.0 nm (SD = 3.0 nm; SEM = 0.36 nm), respectively. Evaluation of previously published EM data shows that the diameters of chromatin from a large number of sources are proportional to linker length. In addition, the inherent variability in fiber diameter suggests a relationship between fiber structure and the heterogeneity of linker length. The cryo-EM data were in quantitative agreement with space-filling double-helical crossed-linker models of Thyone and Necturus chromatin. The data, however, do not support solenoid or twisted-ribbon models for chromatin that specify a constant 30 nm diameter. To reconcile the concept of solenoidal packing with the data, we propose a variable-diameter solid-solenoid model with a fiber diameter that increases with linker length. In principle, each of the variable diameter models for chromatin can be reconciled with local variations in linker length.

A polarized photobleaching study of chromatin reorientation in intact nuclei

Polarized fluorescence recovery after photobleaching (pFRAP) was used to monitor the effects that condensation, i.e. compaction and aggregation, have on the (microseconds and ms) internal dynamics of chromatin in intact nuclei. When divalent cations were present with physiological (approximately 90 mM) monovalent salt the chromatin was found to exist in a compact and aggregated state which was characterized by rotational immobilization over timescales that range from 10 microseconds to 40 milliseconds. This immobilization is attributed to suppression of internal dynamics by intermolecular interactions. When the divalent cations were removed, the compact fibers no longer aggregated and were free to reorient with a characteristic decay time of about 1.2 milliseconds. It is shown that this millisecond relaxation could represent rigid rotation of topologically independent structural domains. Dilution of the monovalent salt induced a gradual change in the structural state of the chromatin that was manifest as a dramatic increase in internal flexibility. At the lowest salt concentration studied (11 mM-monovalent salt) the chromatin reorients in fewer than ten microseconds. These changes in flexibility are continuous with salt concentration, indicating that there are no well-defined endpoints to structural transitions and that the microsecond-millisecond internal dynamics of chromatin are a sensitive measure of structure. Measurements made on nuclei from cells that are either transcriptionally quiescent or active indicate that the dynamics mirrors biological activity.