- Open Access
Hierarchies in eukaryotic genome organization: Insights from polymer theory and simulations
© Iyer et al; licensee BioMed Central Ltd. 2011
- Received: 7 January 2011
- Accepted: 15 April 2011
- Published: 15 April 2011
Eukaryotic genomes possess an elaborate and dynamic higher-order structure within the limiting confines of the cell nucleus. Knowledge of the physical principles and the molecular machinery that govern the 3D organization of this structure and its regulation are key to understanding the relationship between genome structure and function. Elegant microscopy and chromosome conformation capture techniques supported by analysis based on polymer models are important steps in this direction. Here, we review results from these efforts and provide some additional insights that elucidate the relationship between structure and function at different hierarchical levels of genome organization.
- Persistence Length
- Contour Length
- Chromosome Territory
- Chromatin Fiber
- Linker Histone
The term "genome" refers to the complete linear DNA sequence containing all of the hereditary material possessed by an organism. One of the goals of the human genome project was to determine the sequence of the 3 billion base pairs that constitute human DNA. Despite the wealth of information this tour-de-force of scientific enterprise has generated, it is increasingly becoming clear that the cellular function of the human genome is not merely determined by the linear ordering of its DNA base pairs. In fact, many of the functional aspects of the genome are governed by its three dimensional (3D) structure, which involves meters long DNA packaged into the limiting space of a micrometer sized cell nucleus.
The classic image of double stranded DNA (dsDNA) is that of a naked double helix. However, in eukaryotes, dsDNA seldom occurs in its naked form. Instead, it is packaged through a ubiquitous hierarchical process involving specialized proteins called histones. This packaging serves two purposes. First, it compacts the DNA allowing it to fit into the confines of the cell nucleus, and second, it controls the accessibility of DNA to cellular machinery for transcription, regulation, repair and recombination .
The primary structure of the eukaryotic genome consists of DNA wrapping ~1.7 times around histone octamers comprising of two copies of the four histone proteins H2A, H2B, H3, and H4 [4–7]. The combined histone octamer-DNA complex is called the nucleosome. This nucleosomal organization of DNA is considered to be the primary determinant to accessibility of genetic information. Although the atomic structure of the nucleosome has been resolved through X-ray crystallography [6, 7], the dynamics of nucleosomes remains far from fully understood. However, theoretical modeling and simulations are beginning to provide new insights into nucleosomal dynamics by answering questions related to: (1) how histone post-translational modifications and histone variants affect nucleosome structure and intra/inter-nucleosome interactions and (2) how nucleosomes undergo spontaneous conformation transitions between fully- and partially-wrapped states.
Recently, we have  employed all-atom molecular dynamics (MD) simulations and molecular docking to examine the structure of the H4 histone tail and its interactions with the acidic patch of nucleosomes. The H4 tail was found to exhibit a propensity for an α-helical conformation within a stretch of six residues encompassing the lysine 16 (K16) residue. Interestingly, this α-helical region interacted very strongly with the acidic patch; K16, in particular, mediated strong electrostatic interactions with the negatively-charged residues of the acidic patch. Acetylation of K16 diminished these binding interactions, suggesting a plausible mechanism by which post-translational acetylation of K16 could trigger chromatin unfolding . To examine nucleosomal dynamics beyond the time scales accessible to all-atom simulations, Sharma et al.  and Voltz et al.  have developed coarse-grained models of the nucleosome. The former has helped in the identification of important histone residues, termed "cold sites", that maintain the stability of the histone octamer and the latter has been used to compute long wavelength fluctuations of the nucleosome. Some effort has also been devoted to elucidating the energy barriers and the kinetic rate constants associated with the accessibility of DNA in the nucleosome  and explaining specific features of force-induced nucleosome unraveling observed in single-molecule experiments [12–14].
Beyond the single nucleosome level, contiguous nucleosomes separated by short sections of naked DNA called linker DNAs yield the classic beads-on-a-string structure depicted in Figure 2. The secondary structure of chromatin involves folding of this beads-on-a-string motif into a ~30-nm thick fiber called chromatin. While, in vivo imaging of cell nuclei has yielded little information on the secondary structure of chromatin, and has even brought into question the very existence of a 30-nm fiber [15, 16], electron microscopy of isolated nucleosomal arrays have been successful in imaging the transition from 10-nm beads-on-a-string structures to 30-nm condensed structures with increasing salt concentrations [15, 17]. Regardless of the gaps in our knowledge of the secondary structure of chromatin in vivo, this structure is expected to dictate the accessibility of DNA sequences for interactions with nuclear machinery and is also likely to play a role in the recruitment of histone modifying and remodeling factors to specific regions of the genome. Local packing of nucleosomes within the chromatin fiber could also potentially affect its interactions with distant portions of the fiber or other chromatin fibers through processes like interdigitation .
Despite decades of research, the internal structure of the 30-nm chromatin fiber remains controversial. This controversy arises because of two main reasons: (1) in vivo chromatin is too "messy" to be visualized with even the most advanced microscopy techniques  and (2) in vitro reconstituted nucleosome arrays are too large and flexible to be crystallized, and too compact at physiological conditions for their linkers to be fully resolved through microscopy. These limitations have led researchers to adopt indirect ways of deducing the internal structure of chromatin, using chemical cross-linking and single-molecule pulling techniques, with varying degrees of success. Based on the data, two types of models for the chromatin structure have been proposed that differ mainly in the location and configuration of the linkers. In the one-start solenoid model, linkers exhibit a strongly bent configuration and reside at the fiber interior [17, 19, 20]. In the two-start helix, linkers exhibit a straight or gently bent configuration, and they could reside at the periphery of the fiber in one version of the model  or inside the fiber (close to its axis) in a zigzag manner in another version [22–24].
Computational modeling has played an integral role in providing new insights into chromatin architecture. The main challenge in modeling chromatin is the vast degrees of freedom possessed by even small segments of the fiber (e.g., an array of 50 nucleosomes contains >1 million atoms). Hence, all-atom approaches demand prohibitive amounts of computational resources to converge to equilibrium structures. Further, the larger and the more flexible the system, the less adequate a description by one single equilibrium conformation because of the large variation in possible structures around thermal equilibrium. In other words, the systems become "fuzzy" and can only be described in terms of statistical averages. The advantage of this situation is that atomic-detail resolution may easily be abandoned for vast gains in computational speed in a coarse-grained model. To this end, a major focus of computational modeling has been on developing lower-resolution models or "force fields" of nucleosome arrays that still account for the energetic interactions and constraints between the different components of the arrays. A large effort in computational modeling also lies in developing efficient procedures for "sampling" low-energy conformations of the nucleosome array subject to appropriate force fields.
The E2A model is an extension of the two angle model  that also accounts for the cylindrical shape of the nucleosomes . Simulations based on the E2A model are ideally suited for studying long nucleosomal arrays as they overcome the difficulty of lack of knowledge of actual interaction potentials between nucleosomes by using model parameters obtained from experimental data. In particular, the E2A model has been used for studying the effects of variability in linker lengths and variability in linker histone occupancy . In combination with high-resolution light microscopy experiments these simulations promise a way to capture local chromatin structure in the range of 10 - 40 nm resolution .
The next class of models, in addition to including the features of the two-angle model, include the mechanics and electrostatics of the linker DNAs along with a simple treatment of inter-nucleosome interactions [34–36]. The linker DNA is treated as a discretized wormlike chain, possessing energy terms for stretching and bending of the bead-chain with an additional energy term for the relative twist angle between adjacent beads to account for the twisting rigidity of DNA  (Figure 3b). The linker beads are also assigned suitable effective charges that interact with each other via the Debye-Hückel potential . The inter-nucleosome interactions have thus far been treated as hard spheres, interacting with a hard-core attraction  or as cylinders interacting with Gay-Berne [35, 36] and Zewdie potentials . The nucleosome array conformations and packing ratios obtained using Monte Carlo simulations [35, 36] based on these models generally resemble those obtained by electron microscopy for moderately folded chromatin. These models  have also provided valuable insights into the stretching behavior of chromatin and yielded first estimates of the strength of internucleosome interactions by matching simulated force-extension curves to those obtained experimentally . Importantly, the models have revealed the sensitivity of the chromatin structure to the internucleosome interaction strength and linker lengths . Furthermore, it was observed that nucleosome arrays exhibit a zigzag structure without the linker histone and a more solenoidal structure for linker histone-bound arrays .
Very detailed models of local interactions within oligonucleosomes have been developed by Schlick and coworkers [41–49] (Figure 3c). The earliest model [41, 42] treated the nucleosome as a rigid cylinder (along with a small protrusion to represent the H3 histone tail) with hundreds of "pseudo" charges scattered uniformly on the surface. The magnitude of the charges were optimized to reproduce as closely as possible the electric field in the vicinity of the nucleosome. This approach allowed one to account for the salt-dependence in the inter-nucleosome interactions. The model was later refined by employing an irregular-shaped representation of the nucleosome  that was based on a more recent nucleosome crystal structure with all histone tails fully resolved . The model reproduced the experimentally observed [15, 17] compaction of the arrays with increasing salt concentration and indicated that the arrays maintain a zigzag morphology under monovalent salt conditions. Further, the simulations demonstrated that reduced electrostatic repulsion between the linkers is the main mechanism responsible for the folding of arrays at high salt. These models are also being used to study the dynamics of chromatin arrays, especially under different kinds of forces including torsional stresses .
Recently, this model was further improved by accounting for histone tail flexibility , linker histone binding , and effects of divalent ions . The tails were treated as coarse-grained bead-chains, where each bead represented five amino acid residues. The stretching, bending, and the electrostatic terms in the bead-chain were parametrized using an iterative procedure. The linker histone was coarse-grained as three charged beads rigidly bound at the nucleosome dyad with the magnitude of the charges optimized to reproduce the electric field of the atomic linker histone. Divalent ions were treated phenomenologically in terms of their effect on flexibility and electrostatic screening of the linker DNAs. A configurational-bias Monte Carlo approach was used to sample the tail configurations and translation, rotation and pivot moves were used to sample the global array configurations . This model helped elucidate the role of each histone tail, the linker histone and physiological salt condition in chromatin folding [45, 49]. Specifically, the H4 tails were found to mediate the strongest internucleosome interactions, the H3 tails mediated strong internucleosome interactions and screened electrostatic repulsion between the entering/exiting linkers, and the H2A and H2B tails mediated inter-fiber interactions . The linker histones constricted the linker entry/exit angle to bring alternate nucleosomes together. Divalent ions were also found to facilitate tight packing of nucleosomes by allowing a fraction of the linkers to bend and by strongly screening the linker repulsion at the fiber axis. Moreover, the model, in conjunction with sophisticated cross-linking experiments, confirmed the existence of a heteromorphic fiber containing both zigzag and solenoid conformations in the presence of additional divalent cations and linker histones . This model has also recently been used to reproduce the linker length dependence in the observed chromatin structures .
In summary, mesoscale models at varying levels of sophistication such as those discussed above are proving to be valuable tools for examining chromatin structure. The choice of which model to use is dictated by the amount of detail required and the amount of computational resources available. For instance, the detailed models accounting for histone tail flexibility and nucleosome geometry may provide the most accurate representation of short nucleosome arrays. For long arrays containing hundreds of nucleosomes, these models rapidly become computationally intractable and the intermediate-resolution models like the E2A model  become more suitable.
In the previous section, we discussed the secondary structure of chromatin from a static perspective. In reality, the chromatin fiber within the cell nucleus is present in a dynamic state -it is flexible over lengths much larger than the fiber diameter  and it is constantly being subjected to various kinds of remodeling activities, including histone modifications [54–56], sliding and depletion of nucleosomes [57–64], and incorporation of histone variants . There is strong evidence from light microscopy studies indicating that at the gene locus level the chromatin fiber is organized into loops . Studies on several multigene clusters conceive such loops as instrumental in bringing together distant enhancer and promoter regions crucial for gene activation, regulation and recombination [67–69]. Although the detailed mechanism of looping is not fully understood, there is no question that the "intrinsic" bending rigidity of the chromatin fiber dictates to a considerable degree the loop size-dependent statistical probability of two distant regions of the fiber coming into close proximity to form a loop.
Emerging evidence suggests that the relationship between flexibility and looping probability may be utilized by cells for gene regulation. Specifically, the idea that modulation of flexibility through remodeling processes like acetylation of histones alters the probability of interactions between a remote enhancer and cognate gene by means of looping has been successfully applied to explain gene regulation of the Hoxd gene cluster and the β-globin locus . Strong support for this idea also comes from recent experiments in the Murre lab that show large-scale conformational changes in the IgH-locus during B-Cell development accompanying genome-wide and locus-specific histone modifications and nucleosome depletion events [66, 71, 72]. Thus, it seems that chromatin remodeling events, apart from modulating the local structure of chromatin and DNA accessibility, could lead to changes in the higher-order folding of chromatin through its effects on macroscopic properties of the fiber such as its flexibility.
So far, little effort has been devoted to investigating the intrinsic flexibility of chromatin and the associated looping probability as a function of loop size; especially how the flexibility correlates with external conditions like monovalent/divalent salt concentration and system parameters like nucleosome repeat length, DNA wrapping angle, histone variants, histone modifications and presence/absence of linker histones. Chromatin flexibility as characterized by its persistence length, L p , seems to exhibit large variations, depending on the experimental method used for its determination and the conditions and type of chromatin investigated. For example, single-molecule pulling of nucleosome arrays at low salt  and in vivo looping driven cross-linking/recombination assays [53, 73] measure L p as low as 30-50 nm while analysis of fluorescent markers in the genome of erythrocytes measure L p of 100-200 nm [74, 75]. However, it is important to review these results in light of a recent fundamental study  indicating that standard definitions of persistence length as used in these studies may not describe the local intrinsic flexibility of chromatin.
Aumann et al.  recently examined L p of nucleosome arrays using Monte Carlo simulations of a mesoscale model of chromatin. The persistence length was obtained from the decay in the correlation of the tangent vector representing the local fiber axis. It was found that L p decreased strongly with increasing nucleosome repeat length and increasing entry/exit angle, and that binding of the linker histone led to an increase in L p , consistent with experimental observations of linker histone deficient and inclusive chromatin. A comparison between the magnitude of the bending and elastic rigidity suggested that chromatin is much easier to bend than stretch, leading to the interesting hypothesis that it may be easier for the cells to pack chromatin via tight loops rather than by linear compression of the fiber.
We have recently begun to examine the mechanisms behind the conformational collapse in the IgH - locus observed by Murre and coworkers during B-cell development . We hypothesize that chromatin remodeling events (nucleosome depletion and histone modifications [71, 72]) introduce flexible "hinges" within the chromatin fiber causing it to collapse (Figure 4). By treating the chromatin fiber as a worm-like chain (WLC) with fixed contour length L c and variable persistence length L p , we showed that the compaction, as characterized by the ratio of the final to initial mean square end-to-end distances, is given by the ratio of the final to initial persistence length in the ideal chain limit (see Appendix for complete derivation). Using the initial and final L p observed by Heermann and coworkers , this ratio becomes 1/2, indicating that the gene locus compacts by ≈30%. Interestingly, 3D distance measurements in the IgH-locus during progressive stages of B-cell development show a similar ~30% decrease in locus dimensions in response to remodeling events . To examine the functional implications of this collapse, we used a simple Flory type argument to show that the collapse leads to an 8-fold increase in the number of binary interactions n b within the gene locus. Assuming that the frequency of promoter-enhancer (P-E) interactions is proportional to n b , we conclude that the conformational collapse facilitates transcriptional regulation and/or recombination by allowing for higher probability of P-E interactions. Thus, the chromatin remodeling induced changes in the persistence length offers one possible mechanism for governing the conformational state of the gene locus and consequent modification of the functional state of the cell (Figure 4).
Interestingly, the condensed higher-order structures of chromatin, namely chromosomes, were observed as distinct entities during mitosis as early as the 19th century, much before the primary and secondary structures of chromatin were known to exist. However, disappearance of this condensed structure during interphase and the underlying chromatin organization remained a mystery till the late 20th century . In the beginning two different models were proposed for the interphase chromatin organization: (1) random organization akin to a bowl of spaghetti without any apparent structure and (2) organization in territories that later condense during mitosis to form distinct chromosomes . The earliest clues to deciding between the two models came from experiments shining a laser light onto a specific volume of the cell nucleus and observing the effects of the consequent damage on replication . It was found that only a few chromosomes were affected by the laser light, indicating existence of distinct chromosome territories within the cell nucleus. This conjectural evidence has now been confirmed beyond doubt by advances in imaging techniques like fluorescent in situ hybridization (FISH) [82–84] and chromosome painting [82, 83].
While there is concrete evidence for the existence of territories, the internal architecture of chromatin within chromosome territories, the interaction between territories, and the organizational principles governing their formation remain poorly understood . Results from MC simulations of confined polymer chains suggest that territory formation could arise from simple, non-specific entropic forces and from segregation of long chains attempting to conserve their topological state while undergoing confined Brownian motion [86, 87]. Increasingly, insights from such coarse-grained polymer models coupled with experiments are being used to examine the internal architecture of chromosome territories.
Two predominant themes recur during discussions of higher-order structures of interphase chromosomes: (1) formation of loop structures and (2) confined fractal organization. Interestingly, each of these themes derive their support from different kind of experiments. The loop structure theme is predominantly supported by light microscopy experiments [75, 88–90] and the fractal theme by chromosome conformation capture , small angle neutron scattering  and tracer diffusion  experiments. Although the two themes are not mutually exclusive, it is often possible to classify the polymer models based on them. Here, we follow this classification in discussing the polymer models and their central features.
The proposed loop models show a power-law behavior for geometric distances versus contour length at scales <10 Mbp, with scaling exponents ranging from ν = 1/3-1/2. This suggests two different internal structures at work: (1) equilibrium/globular state at short scales and (2) random looped structure at large scales. The transition between these structures is not clearly understood. In contrast, the confined fractals do not encounter these difficulties as they propose a single size exponent at all scales. Recent advances in chromosome conformation capture  have enabled genome-wide measurements of probability of formation of loops as a function of the size of the loops in the human genome (Figure 6). Comparison of the probability measurements with confined freely jointed chain (FJC) simulations indicate that the human genome may be organized as a fractal globule with size scaling exponent ν = 1/3 . However, we note that the agreement of the probability measurements is observed only in the 1-10 Mbp range (Figure 6).
The simple arguments we used for the persistence length calculations inherently assume that the scaling exponent ν = 1/3 is universal within the human genome. However, simulations investigating properties of compact polymers have pointed out that the chromatin fiber, although in a state that resembles the fractal globular state of the compact polymers with ν = 1/3, shows deviations from these type of polymers . Specifically, the moment ratios calculated from intrachain distance distributions for chromatin are observed to be different from that observed for globular state . Instead, the intrachain distributions are observed to be reasonably approximated by the random loop model . However, the random loop model also shows that the distances become independent of genomic separations, i.e., ν → 0. The anomaly in these results can possibly be explained by the argument that in the presence of loops it is inappropriate to use the intrachain distances to determine the scaling exponent ν. Instead, the radius of gyration, a measure of size for closed structures like rings, may be a more appropriate quantity for determining the exponent.
We end this section with a perspective on combined use of the information at the gene locus level and the chromosomal level to infer qualitative aspects governing the functional state of the cell. In the previous section, we have shown that the decrease in persistence length L p causes a collapse and consequent increase in the concentration (density) of segments in a gene locus, via Eqs. (2) and (3). In contrast, in this section, we have shown using Figure 7 that a decrease in L p is accompanied by a decrease in the linear mass density c (since c and Φ are inversely related). The paradoxical nature of this result indicates that chromatin remodeling events, which lower the effective persistence length of chromatin, could have opposing effects on gene-related activities like transcription. On one hand, locus-wide remodeling could open up regions for transcription and lead to increased frequency of interactions between promoters and enhancers; while on the other hand, genome-wide remodeling events could cause genome folding leading to formation of inaccessible, dense regions and potential suppression of transcription. Thus, nucleosome depletion, in addition to its established role in gene activation, may also lead to repression. Such opposing effects are in fact observed when the nucleosomal content is reduced through H4 depletion in the yeast genome, resulting in increased expression in 15% and reduced expression in 10% of genes .
Understanding genome organization within the cell nucleus and its functional implications are fundamental issues in modern cell biology. In this article, we have reviewed the incisive role played by polymer theory and simulations in providing new insights into the structure of the genome at various levels of organization, ranging from packing of nucleosomes into the 30-nm chromatin fiber to the higher-order folding of chromatin via looping into chromosomes. Despite enormous progress in both experimental and computational fronts, a comprehensive model of genome organization at the chromatin and chromosome level is still lacking. Below we list four main challenges that computational models would have to overcome to tackle this highly intriguing and important problem.
First, genome organization is heterogeneous. At the chromatin level, this heterogeneity arises from variations in the nucleosome repeat length, nucleosome composition (histone modifications, histone variants, and linker histones), and potential entrapment of chromatin in metastable states. Most models of chromatin, including those developed in our group, tend to assume uniform nucleosome positioning and composition. Only recently have studies begun to examine the effects of chromatin heterogeneity; preliminary work in this direction has already yielded some very interesting results. At the chromosome level, the size, shape and location of chromatin loops is not fixed. Thus, a key aspect towards examining heterogeneity in chromosome models would be through more realistic modeling of the interactions across chromatin fibers, e.g., allowing loops to dynamically form and break according to the associated energetics of loop formation in chromosome models, as opposed to "fixing" looping points and loop sizes.
Second, the chromatin fiber and chromosomes are highly dynamic. Nucleosomes are constantly being displaced, modified, dissolved and reformed through various mechanisms. At the higher scale, chromatin loops are continuously being broken and reformed. In addition, a range of nuclear proteins including transcription factors and architectural proteins dynamically bind and dissociate from chromatin . Such binding/dissociation events could affect chromatin structure and in turn may be affected by chromatin structure. Most models currently look at the genome from a static perspective. Incorporation of the above mentioned kinetic features into polymer models coupled with structural heterogeneity would be a crucial step in building more comprehensive models of the genome.
Third, development of accurate, coarse-grained models of chromosomes poses a grand challenge. One potential route is via "multi-scale" approaches, where each level represents coarse-graining of the previous, higher-resolution level. Such an approach generally involves the use of potential of mean forces (PMFs) for treating the interactions between coarse-grained subunits. The PMFs account for effects of the degrees of freedom that are "averaged out" during each coarse-graining step. The challenge is to maintain self-consistency from one level to another, as the PMFs are valid only when there is a clear separation of length and time scales across the levels . When there is no clear separation of time/length scales, one might need to resort to "multi-resolution" approaches. In this approach, different portions of the system are treated at different resolutions, depending on their relative importance to the phenomena being investigated (in the same spirit as the QM/MM method  employed in enzyme catalysis). The challenge here is to identify the portions of the system (or degrees of freedom) important enough to be treated at higher resolution as opposed to those regions that can be treated at lower resolution. Another challenge is developing suitable potentials for linking low- and high-resolution regions.
Fourth, little modeling effort has been invested in determining structure-function relationships in genome organization. A definitive connection between structure and function comes from the observation that gene density is high in largely decondensed euchromatin and low in highly condensed heterochromatin . Further, gene-rich and gene-poor regions are found to be physically separated from each other [105, 106]. Specifically, in human lymphocytes gene-rich regions are positioned in the nuclear interior while the gene-poor regions are positioned towards its periphery . Apart from gene density patterns, the molecular species responsible for nuclear processes and gene expression are connected through a complex and extensively coupled network giving rise to functional coupling between macromolecular structures and compartments [108, 109]. Though many aspects of these studies are been driven by the experiments, polymer theory and simulations could contribute to better understanding of such structure-function relationships.
Effect of persistence length modifications on gene locus conformation
To study the generic effects of changes in the persistence length of the chromatin fiber on the conformation of the gene locus, we have treated the chromatin fiber as a worm-like chain (WLC) with fixed contour length L c . The WLC is the standard model used to describe semi-flexible polymer chains, and it has been successfully used to explain the conformational properties of double-stranded DNA and other biopolymers . Though this model is too idealized to represent in vivo chromatin (as it neglects fiber-fiber interactions, existence of stable loops and heterogeneous flexibility), it can nonetheless provide qualitative insights into the magnitude of the conformational changes.
Ratio of mean-square radius of gyration of chains
In the case of the IgH locus, this translates to an 8-fold increase in n b .
Persistence length as a function of folding index
where ν is a scaling exponent that connects the size of the chain to its number of repeat units. Note that the fractal dimension, d f , is the exponent that connects the mass of the chain M to its size via and that M is also proportional to the number of repeat units M ~ N r ~ (C s )1/ν(from Eq. 5). Hence, ν is the inverse of the fractal dimension ν = 1/d f and determines the organization of the FJC at all length scales.
Since the contour length L c = d/c, where d is the genomic distance in bp and c is the linear mass density in bp/nm, Φ increases with decreasing linear mass density.
The contour length within the FJC framework can be partitioned arbitrarily, i.e., if the contour length of a structure is 100 μm then it can be partitioned as, say, 100 repeat units each of 1 μm length or 1000 repeat units of length 0.1 μm. The flexibility of the contour length partitioning allows for choice of either number of repeat units N r or the persistence length L p as a variable in Eq. (6). However, the contour length L c according to Eq. (6) is determined by the product of the folding index Φ and size C s . This constraint has to be satisfied when L p is chosen as the variable and the contour length in the framework changes with the change in the folding index. Thus Φ can be used like a lens of variable resolution to view genome organization at different length scales.
We thank Prof. Cornelis Murre of the Division of Biological Sciences at UCSD for many delightful and insightful discussions on organization of chromosomes.
- Cook PR: The Organization of Replication and Transcription. Science. 1999, 284 (5421): 1790-1795. 10.1126/science.284.5421.1790.Google Scholar
- Misteli T: Beyond the Sequence: Cellular Organization of Genome Function. Cell. 2007, 128: 787-800. 10.1016/j.cell.2007.01.028.Google Scholar
- Lucas JS, Bossen C, Murre C: Transcription and recombination factories: common features?. Curr Opin Cell Biol. 2010, 23: 1-7.Google Scholar
- Kornberg RD, Lorch Y: Twenty-five years of the nucleosome, fundamental particle of the eukaryote chromosome. Cell. 1999, 98: 285-294. 10.1016/S0092-8674(00)81958-3.Google Scholar
- Kornberg RD: Chromatin Structure: A Repeating Unit of Histones and DNA. Science. 1974, 184: 868-871. 10.1126/science.184.4139.868.Google Scholar
- Luger K, Mäder AW, Richmond RK, Sargent DF, Richmond TJ: Crystal structure of the nucleosome core particel at 2.8 Å resolution. Nature. 1997, 389: 251-260. 10.1038/38444.Google Scholar
- Richmond TJ, Davey CA: The structure of DNA in the nucleosome core. Nature. 2003, 423: 145-150. 10.1038/nature01595.Google Scholar
- Yang D, Arya G: Structure and binding of the H4 histone tail and the effects of lysine 16 acetylation. Phys Chem Chem Phys. 2011, 13: 2911-2921. 10.1039/c0cp01487g.Google Scholar
- Shogren-Knaak M, Ishil H, Sun JM, Pazin MJ, Davie JR, Peterson CL: Histone H4-K16 acetylation controls chromatin structure and protein interactions. Science. 2006, 311 (5762): 844-847. 10.1126/science.1124000.Google Scholar
- Sharma S, Ding F, Dokholyan NV: Multiscale modeling of nucleosome dynamics. Biophys J. 2007, 92 (5): 1457-1470. 10.1529/biophysj.106.094805.Google Scholar
- Voltz K, Trylska J, Tozzini V, Kurkal-Siebert V, Langowski J, Smith J: Coarse-Grained Force Field for the Nucleosome from Self-Consistent Multiscaling. J Comp Chem. 2008, 29 (9): 2006-2008. 10.1002/jcc.20902.Google Scholar
- Wocjan T, Klenin K, Langowski J: Brownian Dynamics Simulation of DNA Unrolling from the Nucleosome. J Phys Chem B. 2009, 113 (9): 2639-2646. 10.1021/jp806137e.Google Scholar
- Kulic IM, Schiessel H: DNA spools under tension. Phys Rev Lett. 2004, 92 (22): 228101-10.1103/PhysRevLett.92.228101.Google Scholar
- Sudhanshu B, Mihardja S, Koslover EF, Mehraeen S, Bustamante C, Spakowitz AJ: Tension-dependent structural deformation alters single-molecule transition kinetics. Proc Natl Acad Sci USA. 2011, 108 (5): 1885-1890. 10.1073/pnas.1010047108.Google Scholar
- Woodcock CL, Ghosh RP: Chromatin Higher-order Structure and Dynamics. Cold Spring Harb Perspect Biol. 2010, 2 (5): 10.1101/cshperspect.a000596. a000596-1:22Google Scholar
- Fussner E, Ching RW, Bazett-Jones DP: Living without 30 nm chromatin fibers. Trends Biochem Sci. 2011, 36: 1-6. 10.1016/j.tibs.2010.09.002.Google Scholar
- Thoma F, Koller T, Klug A: Involvement of histone H1 in the organization of the nucleosome and of the salt-dependent superstructures of chromatin. J Cell Biol. 1979, 83: 403-427. 10.1083/jcb.83.2.403.Google Scholar
- Grigoryev SA: Keeping fingers crossed: heterochromatin spreading through interdigitation of nucleosome arrays. FEBS Lett. 2004, 564: 4-8. 10.1016/S0014-5793(04)00258-3.Google Scholar
- Finch JT, Klug A: Solenoidal model for superstructure in chromatin. Proc Natl Acad Sci USA. 1976, 73 (6): 1897-1901. 10.1073/pnas.73.6.1897.Google Scholar
- Widom J, Klug A: Structure of the 300 Å chromatin filament: X-ray diffraction from oriented samples. Cell. 1985, 43: 207-213. 10.1016/0092-8674(85)90025-X.Google Scholar
- Woodcock CLF, Frado LLY, Rattner JB: The higher-order structure of chromatin: evidence for a helical ribbon arrangement. J Cell Biol. 1984, 99: 42-52. 10.1083/jcb.99.1.42.Google Scholar
- Williams SP, Athey BD, Muglia LJ, Schappe RS, Gough AH, Langmore JP: Chromatin fibers are left-handed double helices with diameter and mass per unit length that depend on linker length. Biophys J. 1986, 49: 233-248. 10.1016/S0006-3495(86)83637-2.Google Scholar
- Dorigo B, Schalch T, Kulangara A, Duda S, Schroeder RR, Richmond TJ: Nucleosome arrays reveal the two-start organization of the chromatin fiber. Science. 2004, 306: 1571-1573. 10.1126/science.1103124.Google Scholar
- Schalch T, Duda S, Sargent D, Richmond TJ: X-ray structure of a tetranucleosome and its implications for the chromatin fibre. Nature. 2005, 436: 138-141. 10.1038/nature03686.Google Scholar
- Woodcock CL, Grigoryev SA, Horowitz RA, Whitaker N: A chromatin folding model that incorporates linker variability generates fibers resembling the native structures. Proc Natl Acad Sci USA. 1993, 90 (19): 9021-9025. 10.1073/pnas.90.19.9021.Google Scholar
- van Holde K, Zlatanova J: Chromatin Higher Order Structure: Chasing a Mirage?. J Biol Chem. 1995, 270 (15): 8373-8376. 10.1074/jbc.270.15.8373.Google Scholar
- Schiessel H, Gelbart WM, Bruinsma R: DNA folding: Structural and mechanical properties of the two-angle model for chromatin. Biophys J. 2001, 80: 1940-1956. 10.1016/S0006-3495(01)76164-4.Google Scholar
- Mergell B, Everaers R, Schiessel H: Nucleosome interactions in chromatin: Fiber stiffening and hairpin formation. Phys Rev E. 2004, 70: 10.1103/PhysRevE.70.011915. 011915-1-9)Google Scholar
- Diesinger PM, Heermann DW: Two-angle model and phase diagram for chromatin. Phys Rev E. 2006, 74: 031904-10.1103/PhysRevE.74.031904.Google Scholar
- Wong H, Victor J, Mozziconacci J: An all-atom model of the chromatin fiber containing linker histones reveals a versatile structure tuned by the nucleosomal repeat length. PloS One. 2007, 2 (9): 10.1371/journal.pone.0000877. e877-1:8Google Scholar
- Elena K, Fuller CJ, Straight AF, Spakowitz AJ: Local geometry and elasticity in compact chromatin structure. Biophys J. 2010, 99: 3941-3950. 10.1016/j.bpj.2010.10.024.Google Scholar
- Diesinger PM, Heerman DW: The Influence of the Cylindrical Shape of the Nucleosomes and H1 Defects on Properties of Chromatin. Biophys J. 2008, 94: 4165-4172. 10.1529/biophysj.107.113902.Google Scholar
- Diesinger PM, Heermann DW: Monte Carlo Simulations indicate that Chromatin Nanostructure is accessible by Light Microscopy. PMC Biophys. 2010, 3 (11): 1-20.Google Scholar
- Katritch KV, Bustamante C, Olson WK: Pulling chromatin fibers: computer simulations of direct physical micromanipulations. J Mol Biol. 2000, 295: 29-40. 10.1006/jmbi.1999.3021.Google Scholar
- Wedemann G, Langowski J: Computer simulation of the 30-nanometer chromatin fiber. Biophys J. 2002, 82 (6): 2847-2859. 10.1016/S0006-3495(02)75627-0.Google Scholar
- Stehr R, Kepper N, Rippe K, Wedemann G: The Effect of Internucleosomal Interaction on Folding of the Chromatin Fiber. Biophys J. 2008, 95: 10.1529/biophysj.107.120543.Google Scholar
- Allison SA, Austin R, Hogan M: Bending and twisting dynamics of short linear DNAs: analysis of the triplet anistropy decay of a 209 base pair fragment by Brownian simulation. J Chem Phys. 1989, 90: 3843-3854. 10.1063/1.455790.Google Scholar
- Stigter D: Interactions of highly charged colloidal cylinders with applications to double-stranded DNA. Biopolymers. 1977, 16: 1435-1448. 10.1002/bip.1977.360160705.Google Scholar
- Cui Y, Bustamante C: Pulling a single chromatin fiber reveals the forces that maintain its higher-order structure. Proc Natl Acad Sci USA. 2000, 97: 127-132. 10.1073/pnas.97.1.127.Google Scholar
- Kepper N, Foethke D, Stehr R, Wedemann G, Rippe K: Nucleosome Geometry and Internucleosomal Interactions Control the Chromatin Fiber Conformation. Biophys J. 2008, 95: 3692-3705. 10.1529/biophysj.107.121079.Google Scholar
- Beard D, Schlick T: Modeling salt-mediated electrostatics of macromolecules: The discrete surface charge optimization algorithm and its application to the nucleosome. Biopolymers. 2001, 58: 106-115. 10.1002/1097-0282(200101)58:1<106::AID-BIP100>3.0.CO;2-#.Google Scholar
- Beard D, Schlick T: Computational modeling predicts the structure and dynamics of chromatin fiber. Structure. 2001, 9: 105-114. 10.1016/S0969-2126(01)00572-X.Google Scholar
- Sun J, Zhang Q, Schlick T: Electrostatic mechanism of nucleosomal array folding revealed by computer simulation. Proc Natl Acad Sci USA. 2005, 102: 8180-8185. 10.1073/pnas.0408867102.Google Scholar
- Arya G, Schlick T: Role of histone tails in chromatin folding revealed by a mesoscopic oligonucleosome model. Proc Natl Acad Sci USA. 2006, 103 (44): 16236-16241. 10.1073/pnas.0604817103.Google Scholar
- Arya G, Zhang Q, Schlick T: Flexible histone tails in a new mesoscopic oligonucleosome model. Biophys J. 2006, 91: 133-150. 10.1529/biophysj.106.083006.Google Scholar
- Arya G, Schlick T: Efficient global biopolymer sampling with end-transfer configurational bias Monte Carlo. J Chem Phys. 2007, 126: 044107-10.1063/1.2428305.Google Scholar
- Arya G, Schlick T: A tale of tails: how histone tails mediate chromatin compaction in different salt and linker histone environments. J Phys Chem A. 2009, 113: 4045-4059. 10.1021/jp810375d.Google Scholar
- Schlick T, Perisić O: Mesoscale simulations of two nucleosome-repeat length oligonucleosomes. Phys Chem Chem Phys. 2009, 11 (45): 10729-10737. 10.1039/b918629h.Google Scholar
- Grigoryev SA, Arya G, Correll S, Woodcock CL, Schlick T: Evidence for heteromorphic chromatin fibers from analysis of nucleosome interactions. Proc Natl Acad Sci USA. 2009, 106 (32): 13317-13322. 10.1073/pnas.0903280106.Google Scholar
- Davey CA, Sargent DF, Luger K, Maeder AW, Richmond TJ: Solvent-mediated interactions in the structure of the nucleosome core particle at 1.9 Åresolution. J Mol Biol. 2002, 319: 1097-1113. 10.1016/S0022-2836(02)00386-8.Google Scholar
- Dobrovolskaia IV, Kenward M, Arya G: Twist propagation in dinucleosome arrays. Biophys J. 2010, 99: 3355-3364. 10.1016/j.bpj.2010.09.055.Google Scholar
- Deal R, Henikoff S: Capturing the dynamic epigenome. Genome Biol. 2010, 11 (10): 218-10.1186/gb-2010-11-10-218.Google Scholar
- Ringrose L, Chabanis S, Angrand PO, Woodroofe C, Stewart AF: Quantitative comparison of DNA looping in vitro and in vivo: chromatin increases effective DNA flexibility at short distances. EMBO J. 1999, 18: 6630-6641. 10.1093/emboj/18.23.6630.Google Scholar
- Bruno M, Flaus A, Stockdale C, Rencurel C, Helder Ferreira, Owen Hughes: Histone H2A/H2B dimer exchange by ATP-dependent chromatin remodeling activities. Mol Cell. 2003, 12: 1599-1606. 10.1016/S1097-2765(03)00499-4.Google Scholar
- Kouzarides T: Chromatin Modifications and Their Function. Cell. 2007, 128: 693-705. 10.1016/j.cell.2007.02.005.Google Scholar
- Yang X, Zaurin R, Beato M, Peterson CL: Swi3p controls SWI/SNF assembly and ATP-dependent H2A-H2B displacement. Nat Struct Mol Biol. 2007, 14: 540-547. 10.1038/nsmb1238.Google Scholar
- Lorch Y, Zhang M, Kornberg RD: Histone Octamer Transfer by a Chromatin-Remodeling Complex. Cell. 1999, 96: 389-392. 10.1016/S0092-8674(00)80551-6.Google Scholar
- Whitehouse I, Flaus A, Cairns BR, White MF, Workman JL, Owen-Hughes T: Nucleosome mobilization catalysed by the yeast SWI/SNF complex. Nature. 1999, 400: 784-787. 10.1038/23506.Google Scholar
- Hamiche A, Sandaltzopoulos R, Gdula DA, Wu C: ATP-Dependent Histone Octamer Sliding Mediated by the Chromatin Remodeling Complex NURF. Cell. 1999, 97: 833-842. 10.1016/S0092-8674(00)80796-5.Google Scholar
- Längst G, Bonte EJ, Corona DFV, Becker PB: Nucleosome Movement by CHRAC and ISWI without disruption or trans-displacement of theHhistone Octamer. Cell. 1999, 97: 843-852.Google Scholar
- Boeger H, Griesenbeck J, Strattan J, Kornberg RD: Nucleosomes unfold completely at a transcriptionally active promoter. Mol Cell. 2003, 11: 1587-1598. 10.1016/S1097-2765(03)00231-4.Google Scholar
- Reinke H, Hörz W: Histones are first hyperacetylated and then lose contact with the activated PHO5 promoter. Mol Cell. 2003, 11: 1599-1607. 10.1016/S1097-2765(03)00186-2.Google Scholar
- Boeger H, Griesenbeck J, Strattan J, Kornberg R: Removal of promoter nucleosomes by disassembly rather than sliding in vivo. Mol Cell. 2004, 14: 667-673. 10.1016/j.molcel.2004.05.013.Google Scholar
- Cairns BR: Chromatin remodeling: insights and intrigue from single-molecule studies. Nat Struct Mol Biol. 2007, 14 (11): 989-996. 10.1038/nsmb1333.Google Scholar
- Mizuguchi G, Shen X, Landry J, Wu WH, Sen S, Wu C: ATP-driven exchange of histone H2AZ variant catalyzed by SWR1 chromatin remodeling complex. Science. 2004, 303: 343-348. 10.1126/science.1090701.Google Scholar
- Jhunjhunwala S, van Zelm MC, Peak MM, Cutchin S, Roy R, van Dongen JJM, Grosveld FG, Knoch TA, Murre C: The 3D Structure of the Immunoglobulin Heavy-Chain Locus: Implications for Long-Range Genomic Interactions. Cell. 2008, 133: 265-279. 10.1016/j.cell.2008.03.024.Google Scholar
- Carter D, Chakalova L, Osborne CS, Dai Y, Fraser P: Long-range chromatin regulatory interactions in vivo. Nat Genet. 2002, 32 (4): 623-626. 10.1038/ng1051.Google Scholar
- Tolhuis B, Palstra R, Splinter E, Grosveld F, de Laat W: Looping and interaction between hypersensitive sites in the active beta-globin locus. Mol Cell. 2002, 10: 1453-1465. 10.1016/S1097-2765(02)00781-5.Google Scholar
- Sayegh C, Jhunjhunwala S, Riblet R, Murre C: Visualization of looping involving the immunoglobulin heavy-chain locus in developing B cells. Gene Dev. 2005, 19: 322-327. 10.1101/gad.1254305.Google Scholar
- Li Q, Barkess G, Qian H: Chromatin looping and the probability of transcription. Trends Genet. 2006, 22 (4): 197-202. 10.1016/j.tig.2006.02.004.Google Scholar
- Heinz S, Benner C, Spann N, Bertolino E, Lin YC, Laslo P, Cheng JX, Murre C, Singh H, Glass CK: Simple Combinations of Lineage-Determining Factors Prime cis-Regulatory Elements Required for Macrophage and B-Cell Identities. Mol Cell. 2010, 38: 576-589. 10.1016/j.molcel.2010.05.004.Google Scholar
- Lin YC, Jhunjhunwala S, Benner C, Heinz S, Welinder E, Mansson R, Sigvardsson M, Hagman J, Espinoza CA, Dutkowski J, Ideker T, Glass CK, Murre C: A global network of transcription factors, involving E2A, EBF1 and Foxo1, that orchestrates B cell fate. Nat Immunol. 2010, 11 (7): 635-643. 10.1038/ni.1891.Google Scholar
- Dekker J, Rippe K, Dekker M, Kleckner N: Capturing Chromosome Conformation. Science. 2002, 295: 1306-1311. 10.1126/science.1067799.Google Scholar
- Bystricky K, Heun P, Gehlen L, Langowski J, Gasser SM: Long-range compaction and flexibility of interphase chromatin in budding yeast analyzed by high-resolution imaging techniques. Proc Natl Acad Sci USA. 2004, 101: 16495-6500. 10.1073/pnas.0402766101.Google Scholar
- Sachs RK, van den Engh G, Trask B, Yokota H, Hearst JE: A random-walk/giant-loop model for inter-phase chromosomes. Proc Natl Acad Sci USA. 1995, 92: 2710-2714. 10.1073/pnas.92.7.2710.Google Scholar
- Hsu HP, Paul W, Binder K: Standard Definitions of Persistence Length Do Not Describe the Local "Intrinsic" Stiffness of Real Polymer Chains. Macromolecules. 2010, 43 (6): 3094-3102. 10.1021/ma902715e.Google Scholar
- Aumann F, Lankas F, Caudron M, Langowski J: Monte Carlo simulation of chromatin stretching. Phys Rev E. 2006, 73: 041927-10.1103/PhysRevE.73.041927.Google Scholar
- Diesinger PM, Heermann DW: Depletion Effects Massively Change Chromatin Properties and Influence Genome Folding. Biophys J. 2009, 97: 2146-2153. 10.1016/j.bpj.2009.06.057.Google Scholar
- Diesinger PM, Kunkel S, Langowski J, Heermann DW: Histone depletion facilitates chromatin loops on the kilobasepair scale. Biophys J. 2010, 99 (9): 2995-3001. 10.1016/j.bpj.2010.08.039.Google Scholar
- Razin SV, Petrov A, Hair A, Vassetzky YS: Chromatin Domains and Territories: Flexibly Rigid. Crit Rev Eukar Gene. 2004, 14 (1-2): 79-88. 10.1615/CritRevEukaryotGeneExpr.v14.i12.50.Google Scholar
- Cremer T, Cremer C, Baumann H, Luedtke EK, Sperling K, Teuber V, Zorn C: Rabl's Model of the Inter-phase Chromosome Arrangement Tested in Chinese Hamster Cells by Premature Chromosome Condensation and Laser-UV-Microbeam Experiments. Hum Genet. 1982, 60: 46-56. 10.1007/BF00281263.Google Scholar
- van der Ploeg M: Cytochemical nucleic acid research during the twentieth century. Eur J Histochem. 2000, 44: 7-42.Google Scholar
- Cremer T, Cremer C: Chromosome territories, nuclear architecture and gene regulation in mammalian cells. Nat Rev Genet. 2001, 2: 292-301. 10.1038/35066075.Google Scholar
- Parada LA, Misteli T: Chromosome positioning in the interphase nucleus. Trends Cell Biol. 2002, 12 (9): 425-432. 10.1016/S0962-8924(02)02351-6.Google Scholar
- Cremer T, Cremer M, Dietzel S, Müller S, Solovei I, Fakan S: Chromosome territories - a functional nuclear landscape. Curr Opin Cell Biol. 2006, 18: 307-316. 10.1016/j.ceb.2006.04.007.Google Scholar
- Rosa A, Everaers R: Structure and Dynamics of Interphase Chromosomes. PLOS Comput Biol. 2008, 4 (8): e1000153-1-10. 10.1371/journal.pcbi.1000153.MathSciNetGoogle Scholar
- Cook PR, Marenduzzo D: Entropic organization of interphase chromosomes. J Cell Biol. 2009, 186 (6): 825-834. 10.1083/jcb.200903083.Google Scholar
- van den Engh G, Sachs RK, Trask BJ: Estimating genomic distance from DNA sequence location in cell nuclei by a random walk model. Science. 1992, 257: 1410-1412. 10.1126/science.1388286.Google Scholar
- Münkel C, Eils R, Dietzel S, Zink D, Mehring C, Wedemann G, Cremer T, Langowski J: Compartmentalization of interphase chromosomes observed in simulation and experiment. J Mol Biol. 1999, 285: 1053-1065.Google Scholar
- Mateos-Langerak J, Bohn M, de Leeuw W, Giromus O, Manders EMM, Verschure PJ, Indemans MH, Gierman HJ, Heermann DW, van Driel R, Goetze S: Spatially confined folding of chromatin in the interphase nucleus. Proc Natl Acad Sci USA. 2009, 106 (10): 3812-3817. 10.1073/pnas.0809501106.Google Scholar
- Lieberman-Aiden E, van Berkum NL, Williams L, Imakaev M, Ragoczy T, Telling A, Amit I, Lajoie BR, Sabo PJ, Dorschner MO, Sandstrom R, Bernstein B, Bender MA, Groudine M, Gnirke A, Stamatoyannopoulos J, Mirny LA, Lander ES, Dekker J: Comprehensive Mapping of Long-Range Interactions Reveals Folding Principles of the Human Genome. Science. 2009, 326 (5950): 289-293. 10.1126/science.1181369.Google Scholar
- Lebedev DV, Filatov MV, Kuklin AI, Islamov AK, Kentzinger E, Pantina R, Toperverg BP, Isaev-Ivanov VV: Fractal nature of chromatin organization in interphase chicken erythrocyte nuclei: DNA structure exhibits biphasic fractal properties. FEBS Lett. 2005, 579: 1465-1468. 10.1016/j.febslet.2005.01.052.Google Scholar
- Bancaud A, Huet S, Daigle N, Mozziconacci J, Ellenberg J: Molecular crowding affects diffusion and binding of nuclear proteins in heterochromatin and reveals. EMBO J. 2009, 28 (24): 3785-3798. 10.1038/emboj.2009.340.Google Scholar
- Hahnfeldt P, Hearst JE, Brenner DJ, Sachs RK, Hlatky LR: Polymer models for interphase chromosomes. Proc Natl Acad Sci USA. 1993, 90 (16): 7854-7858. 10.1073/pnas.90.16.7854.Google Scholar
- Ostashevsky J, Lange C: The 30 nm chromatin fiber as a felxible polymer. J Biomol Struct Dyn. 1994, 11: 813-820.Google Scholar
- Münkel C, Langowski J: Chromosome structure predicted by a polymer model. Phys Rev E. 1998, 57 (5): 5888-5896.Google Scholar
- Bohn M, Heermann DW, van Driel R: A Random Loop Model for Long Polymers. Phys Rev E. 2007, 76: 051805-1-051805-8. 10.1103/PhysRevE.76.051805.Google Scholar
- Bohn M, Heermann DW: Topological interactions between ring polymers: Implications for chromatin loops. J Chem Phys. 2010, 132 (4): 10.1063/1.3302812. 044904-1:6Google Scholar
- Rodley CDM, Bertels F, Jones B, O'Sullivan JM: Global identification of yeast chromosome interactions using Genome conformation capture. Fungal Genet Biol. 2009, 46 (11): 879-886. 10.1016/j.fgb.2009.07.006.Google Scholar
- Bohn M, Heermann DW: Conformational properties of compact polymers. J Chem Phys. 2009, 130: 174901-1-174901-11. 10.1063/1.3126651.Google Scholar
- Wyrick JJ, Holstege FCP, Jennings EG, Causton HC, Shore D, Grunstein M, Lander ES, Young RA: Chromosomal landscape of nucleosome-dependent gene expression and silencing in yeast. Nature. 1999, 402: 418-421. 10.1038/46567.Google Scholar
- Müller-Plathe F: Coarse-Graining in Polymer Simulation: From the Atomistic to the Mesoscopic Scale and Back. Chemphyschem. 2002, 3 (9): 755-769.Google Scholar
- Warshel A, Levitt M: Theoretical studies of enzymic reactions: Dielectric, electrostatic and steric stabilization of the carbonium ion in the reaction of lysozyme. J Mol Biol. 1976, 103: 227-249. 10.1016/0022-2836(76)90311-9.Google Scholar
- Gilbert N, Boyle S, Fiegler H, Woodfine K, Carter NP, Bickmore WA: Chromatin architecture of the human genome: gene-rich domains are enriched in open chromatin fibers. Cell. 2004, 118 (5): 555-566. 10.1016/j.cell.2004.08.011.Google Scholar
- Boutanaev AM, Mikhaylova LM, Nurminsky DI: The Pattern of Chromosome Folding in Interphase Is Outlined by the Linear Gene Density Profile. Mol Cell Biol. 2005, 25 (18): 8379-8386. 10.1128/MCB.25.18.8379-8386.2005.Google Scholar
- Shopland LS, Lynch CR, Peterson KA, Thornton K, Kepper N, von Hase J, Stein S, Vincent S, Molloy KR, Kreth G, Cremer C, Bult CJ, O'Brien TP: Folding and organization of a contiguous chromosome region according to the gene distribution pattern in primary genomic sequence. J Cell Biol. 2006, 174: 27-38. 10.1083/jcb.200603083.Google Scholar
- Boyle S, Gilchrist S, Bridger JM, Mahy NL, Ellis JA, Bickmore WA: The spatial organization of human chromosomes within the nuclei of normal and emerin-mutant cells. Hum Mol Genet. 2001, 10 (3): 211-219. 10.1093/hmg/10.3.211.Google Scholar
- Maniatis T, Reed R: An extensive network of coupling among gene expression machines. Nature. 2002, 416: 499-506. 10.1038/416499a.Google Scholar
- O'Brien TP, Bult CJ, Cremer C, Grunze M, Knowles BB, Langowski J, McNally J, Pederson T, Politz JC, Pombo A, Schmahl G, Spatz JP, van Driel R: Genome Function and Nuclear Architecture: From Gene Expression to Nanoscience. Genome Res. 2003, 13: 1029-1041.Google Scholar
- Lu Y, Weers B, Stellwagen NC: DNA persistence length revisited. Biopolymers. 2002, 61: 261-275. 10.1002/bip.10151.Google Scholar
- Langowski J, Heermann DW: Computational modeling of the chromatin fiber. Semin Cell Dev Biol. 2007, 18 (5): 659-667. 10.1016/j.semcdb.2007.08.011.Google Scholar
- de Gennes PG: Scaling Concepts in Polymer Physics. 1985, Ithaca and London: Cornell University PressGoogle Scholar
- Rubinstein M, Colby RH: Polymer Physics. 2003, New York: Oxford University PressGoogle Scholar
- Langowski J: Polymer chain models of DNA and chromatin. Eur Phys J E Soft Matter. 2006, 19: 241-249. 10.1140/epje/i2005-10067-9.Google Scholar
- Woodcock CL, Dimitrov S: Higher-order structure of chromatin and chromosomes. Curr Opin Genet Dev. 2001, 11 (2): 130-135. 10.1016/S0959-437X(00)00169-6.Google Scholar
- Woodcock CL: Chromatin architecture. Curr Opin Struc Biol. 2006, 16 (2): 213-220. 10.1016/j.sbi.2006.02.005.Google Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.