The close interaction between the 3D interphase DNA structure and gene expression has made chromatin folding a rapidly developing field of study. Several previous reports have been continuously challenged by the progress of research. For example, in the solenoid model, chromatin is folded into a 30-nm fibre, which is assembled into higher-order structures
[153]. However, this report has not been substantiated when studied using electron microscopy tomography, which has shown highly distorted chromatin polymers
[154]. We have discussed the role of TADs and the genomic contact and loop extrusion hypothesis based on CTCF and cohesin. It was previously thought that condensin compacts chromosomes by randomly bridging the DNA segment or supercoiling or passively pushing to sites of convergent transcription by RNA polymerase
[155]. However, direct observation of loop extrusion in a single molecule has been found in condensin motor activity when quantum-dot-tagged yeast condensin was translocated along double-tethered DNA curtains
[156]. The loop extrusion in the naked DNA was reported to be faster than that in DNA polymerase. While the loop extrusion of naked DNA translocation proceeded at a speed of 0.5–2 kb/s, those of DNA and RNA polymerases proceeded at approximately 1 kb/s and 1 kb/min, respectively
[157][158]. The SMC complex uses ATP hydrolysis to perform the loop extrusion at a rate of 0.1–2 s
−1 [155]. It has also been reported that cohesin- or condensin-binding factors possibly reduce the rate of chromatin loop extrusion. The major factors are the 11-nm nucleosome, RNA polymerase, protein complexes, and DNA structures
[155]. However, the SMC complex can avoid these obstacles through nontopological binding, involving intermittent interactions with Nipbl1 and Pds5 proteins that alter the extrusion dynamics, where Nipbl1 possibly acts as a ‘dynamic safety belt’ for the cohesin protein
[155]. A recent study has explored these aspects and reported the presence of TAD-like clusters even after cohesin knockout
[159][160]. These results suggest that large cooperation of the architectural regulatory proteins, as well as the interplay of supercoiling, molecular binding, phase separation, crowding effect, and loop extrusion events, is needed. There are also questions regarding the functional units of chromatin and their hierarchy of folding, and the inner functional units of working TADs at the single-cell level. To understand all these intricate events, researchers have applied mathematical rules (stochastic, self-returning event) and studied a folding algorithm that can replicate experimental observations
[161]. The most common type of chromatin interaction in the genome is that of the promoter and enhancer for transcriptional regulation and heterogeneous packing, which disperses local DNA accessibility and allows transcription and nuclear transport. From a polymer physics point of view, there is an apparent conflict between these two chromatin properties. It has been reported that chromatin resembles a fractal globule, which is a self-similar polymer in a collapsed state
[162][163]. Although the fractal globule model observes high contact frequencies, it does not explain the spatial heterogeneity of chromatin packaging, which a 1D polymer cannot provide
[161]. Therefore, researchers have provided a self-returning random walk (SRRW) mathematical model to address the contact-structure paradox
[161]. It provides a non-branching topology of the 10-nm chromatin fiber and generates tree-like topological domains connected to an open chromatin backbone
[161]. The SRRW model, presenting a new picture of genome organization, supports the hypothesis that local DNA density plays a critical role as a transcriptional regulator; the chromatin folds into a variety of minimally entangled hierarchical structures across the length from nanometers to micrometers without the necessity of a 30-nm fiber
[161]. This model also explains the structure–function relationship of the interphase DNA with higher-order folding and a substantial reduction in dimension during genomic landscape exploration
[161]. The model also predicts that the topological domains in single cells contain random-tree structures, where tree domains are regarded as nanoclusters and loops on a kilobase-to-megabase scale, serving as building blocks for large packaging domains. These tree domains are called ‘3D forests’ within the chromosome territory
[161]. The size of a tree domain is positively correlated with the size of the genome, with considerable depression
[161]. There is also a positive correlation between the tree domain size and packaging density, suggesting a size-dependent domain activity, where the nanodomain of the peak radius is approximately 70 nm
[161]. Additionally, the model predicts a correlation between local DNA density and domain size, supporting the view that small domains are more active than large domains
[161]. The first-order genome of a double helix DNA evolves to adopt a ‘virtual tree data structure’ for higher-order genome organization
[161]. This tree-like topological domain is connected by an open functional backbone segment, which facilitates the proper organization of genomic contacts, package-based regulation of transcription, transport and accommodation of nuclear proteins, and transition between the interphase and mitosis
[161].
A computational string and binder (SBS) model was proposed in polymer physics to understand the mechanism of chromosome compartmentalization, pattern formation, and chromatin folding
[164]. According to the SBS model, chromatin folding can be driven thermodynamically by homotypic interactions between DNA sites that share compatible chromatin marks
[164][165][166][167][168]. The chromatin filament acts as a self-avoiding walk string of beads, where specific beads function as binding sites for a cognate diffusing binder that can bridge them to allow folding
[164]. The different binding sites can selectively interact with their cognate binder, and these binding activities can be computationally investigated by molecular dynamics simulations using Langevin dynamics with classical interaction potential
[164][168]. This model explains chromatin folding thermodynamically by homotypic interactions between DNA sites sharing cognate chromatin
[166][167][168][169][170]. This interaction takes place via protein binding to multiple sites, inducing phase separation of chromatin sub-compartments
[171]. The association between chromatin sites and the nuclear lamina and speckles can also be inferred using the SBS model, with the help of the bridging protein transcription factor YY1, RNA polymerase II, and Polycomb repressive complex 1
[164][172].