Nucleus

The spatiotemporal organisation of chromatin in the eukaryotic nucleus is fundamental for genome regulation. Chromatin undergoes rapid remodelling and rearrangements within minutes, altering its diffusion properties. Considering the tight coupling between genome function and nuclear architecture, a key question is how chromatin dynamics adapt to or promote nuclear processes. To elucidate the underlying physical principles, we employed High-resolution Diffusion mapping (Hi-D) to track chromatin motion throughout interphase in live human cells. Our analysis, that considers both diffusive motion and drift generated by active forces, revealed that chromatin dynamics are heterogeneous, with distinct behaviours in different subnuclear zones. Notably, both diffusive and processive contributions to chromatin motion progressively decrease from G1 to G2 phase, with this reduction occurring uniformly across all subzones. This suggests a global mechanism driving the observed decrease in chromatin mobility during cell cycle progression. By combining genetic knockout experiments and polymer modelling, we demonstrate that the doubling of DNA content, rather than cohesin-mediated sister chromatid entrapment, is responsible for the gradual decrease in chromatin motion during the cell cycle in human nuclei. These findings provide new insights into the physical and functional organisation of chromatin and its regulation during cellular proliferation. GRAPHICAL ABSTRACT O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=78 SRC="FIGDIR/small/712877v1_ufig1.gif" ALT="Figure 1"> View larger version (26K): org.highwire.dtl.DTLVardef@75c654org.highwire.dtl.DTLVardef@2fbd3dorg.highwire.dtl.DTLVardef@31025aorg.highwire.dtl.DTLVardef@191808e_HPS_FORMAT_FIGEXP M_FIG C_FIG

Nucleus shape is a sensitive indicator of cell state, influenced by numerous bio-chemical and physiological factors. While prior work has cataloged how perturbations alter nucleus morphology, we address the inverse: inferring underlying molecular changes from nucleus shape alone. We previously developed a mechanical model yielding two nondimensional parameters: flatness index and scale factor, which are surrogate measures for cortical actin tension and nuclear envelope compliance respectively. In this study, we apply these parameters to investigate the dynamics in cellular mechanics during confined migration. We fabricated polydimethylsiloxane (PDMS) microchannels with widths of 3 {micro}m (high confinement) and 10 {micro}m (low confinement) and tracked cells migrating through them. We captured high-frequency 3D nucleus shapes via double fluorescence exclusion microscopy and custom image analysis. Fitting the model and estimating flatness index and scale factor to time-resolved shapes revealed dynamic regulation in 3 {micro}m channels: actin tension decreased and nucleus compliance increased immediately before nucleus entry into the constriction, with rapid restoration to baseline upon exit. No such changes occurred in 10 {micro}m channels, indicating active, confinement-dependent cytoskeletal adaptation. Immunostaining for YAP and lamin-A,C confirmed these model inferences. Our results uncover mechanostasis, active mechanical homeostasis, during confined migration and establish the combination of double fluorescence exclusion microscopy and nondimensional nucleus shape parameters as a powerful, non-invasive tool for single-cell mechanobiology studies.

Understanding how chromatin folds in three dimensions remains challenging because most experimental assays capture low-dimensional projections of an underlying, highly heterogeneous polymer. Here we present an ensemble-based interpretive framework built on the previously introduced Self-Returning Excluded Volume (SR-EV) model, a minimal generator of nucleosome-resolution chromatin conformations based on stochastic return rules and excluded-volume geometry. Despite its simplicity, SR-EV reproduces key experimental signatures across scales: heterogeneous nanoscale packing domains resembling ChromEMT and ChromSTEM observations, sparse and highly variable single-configuration contact patterns analogous to single-cell chromosome conformation capture (Hi-C), and robust ensemble-level contact enrichment consistent with topologically associating domains (TADs). In this framework, Hi-C loop and TAD signatures are interpreted as ensemble-level statistical enrichments rather than invariant features of single-cell conformations. SR-EV is explicitly designed to generate large ensembles of complete three-dimensional chromatin configurations that can be projected consistently onto two-dimensional contact maps and one-dimensional genomic profiles. By introducing architectural-protein effects only through ensemble selection rather than explicit forces, SR-EV supports a separation between intrinsic polymer geometry and regulatory bias and suggests that TAD-like features can emerge as statistical enrichments rather than deterministic three-dimensional structures. Coordination number and probe-based accessibility computed directly from SR-EV provide a unified link between three-dimensional packing, two-dimensional contact maps, and one-dimensional genomic profiles. Together, these results establish SR-EV as a minimal and physically grounded reference framework for interpreting how heterogeneous chromatin ensembles give rise to multimodal experimental observables, while remaining consistent with the fact that chromatin organization is realized in individual cells. SIGNIFICANCEChromatin domains, boundaries, and contact enrichments are often interpreted as fixed structural entities, even though most experimental measurements average over large and heterogeneous cell populations. The SR-EV framework shows that many of these features can be understood as emerging from minimal geometric rules combined with ensemble-level bias, without requiring explicit molecular interactions or deterministic folding mechanisms. By distinguishing single-configuration heterogeneity from ensemble-level statistical organization - including the emergence of packing domains- SR-EV supports an interpretation in which chromatin organization is realized in individual cells but must be analyzed through ensembles. This perspective clarifies the probabilistic nature of genome architecture and provides a tractable reference framework for interpreting multimodal genomic and imaging data.

Quantifying chromatin-state dynamics in living cells remains challenging, in part because most methods require fixation or cell lysis. Here, we benchmark and introduce three simple live-cell image-derived metrics computed from routine DNA staining--the coefficient of variation (CV), 1-Gini, and the Diffuse Signal Index (DSI), introduced here--as fixation-free readouts of chromatin state. Using HL60-derived neutrophils (dHL-60) undergoing NETosis as a model system with a pronounced compact-to-decompact chromatin transition, we show that all three metrics track progressive chromatin reorganization in live-cell trajectories, but differ markedly in sensitivity: DSI provides the strongest trajectory-level discrimination between NETing and non-NETing cells, followed by 1-Gini and CV. Comparison with Tn5-based chromatin accessibility measurements in fixed cells further shows that all three metrics correlate with chromatin accessibility, supporting their biological relevance. Together, our results provide a practical framework for extracting chromatin-state readouts from routine live-cell DNA staining and identify DSI as the most discriminative metric for tracking chromatin reorganization in this benchmark.

The transcription factor Yin Yang 1 (YY1) plays roles in chromatin organization, combining sequence-specific DNA recognition via zinc finger domains with multivalent interactions mediated by intrinsically disordered regions. While YY1 has been implicated in phase separation and enhancer-promoter communication, how its structured and intrinsically disordered regions cooperate to shape DNA-protein assemblies remains unclear. Here, we used single-molecule DNA curtain fluorescence imaging to dissect the molecular basis of YY1-DNA assembly. We found that YY1 induces higher-order assembly in a concentration-dependent manner. At moderately high concentrations, YY1 formed weakly-linked DNA condensates in which dynamic, liquid-like YY1 molecules are scaffolded by relatively immobile, solid-like DNA. At high concentrations, zinc finger-mediated bridging dominates, producing strongly-linked condensates. Distinct domain-deletion mutants selectively impaired either weakly-(soft) or strongly-linked (hard) DNA condensate formation, indicating that these architectures are generated by separate domain-dependent mechanisms. Our findings establish a domain-level framework for DNA condensate formation and highlight how transcription factors can integrate specific and non-specific DNA interactions to control the material state of chromatin and would influence genome regulation.

The PBAF is one of three biochemically distinct BAF chromatin remodelers in humans. We previously proposed the role of ARID2, a PBAF component, as a bonafide tumor suppressor in colorectal cancer (CRC). Here, we validated loss of tumor suppression under conditions of ARID2 deficiency emanating from a marked reduction in PBAF complex assembly resulting from destabilization of PBAF-specific components BRD7, PHF10, and PBRM1. Transcriptome profiling of ARID2 deficient CRC cells revealed perturbation of disease processes, including CRC and neurodegenerative disorders, as well as CRC relevant pathways including Wnt/{beta}-catenin signalling, but transcript levels of PBAF-specific components remained unchanged, confirmed by RT-qPCR and TCGA data analysis. Our study establishes ARID2 as a critical stabilizer of the PBAF complex of relevance to CRC.

Inferring the physical interaction landscape underlying chromatin contact maps remains a central challenge in genome biology. Chromosome conformation capture experiments such as Micro-C provide high-resolution measurements of spatial contacts between genomic loci, yet translating these contact frequencies into an interpretable structural interaction model is a fundamentally underdetermined inverse problem. Here we present a nucleosome-resolution maximum entropy framework that infers effective pairwise interaction parameters directly from experimental Micro-C data. In this approach, chromatin is represented as a heterogeneous nucleosome-linker polymer, and maximum entropy optimization identifies the minimal set of interaction constraints required to reproduce the observed contact statistics. Applied to several human genomic loci, the inferred interaction landscape accurately reconstructs experimental contact maps and remains robust to substantial masking or perturbation of the input data. Forward simulations using the inferred parameters generate structural ensembles that reproduce the observed chromatin organization, demonstrating that the model captures the generative constraints governing chromatin folding rather than merely fitting contact frequencies. Analysis of these ensembles reveals domain boundaries, localized interaction hotspots and emergent chromatin "blobs" that coincide with independently observed architectural features. Together, these results establish a highresolution statistical inference framework that reconstructs locus-scale chromatin structure while providing a physically interpretable interaction landscape linking chromatin contact maps to three-dimensional genome organization.

DNA is the largest biopolymer in nature, and chromatin contact maps are widely interpreted as quantitative readouts of its three-dimensional organization. However, the validity of such interpretations critically depends on how these maps are processed. Here, we identify a previously overlooked but fundamental source of bias in chromatin contact data analysis. We demonstrate that a widely adopted preprocessing convention, namely whole-matrix percentile clipping, systematically distorts sparse contact maps by collapsing their dynamic range. This effect is strongest in near-diagonal interactions, precisely the regime encoding chromatin domains and looping structures, thereby compromising quantitative interpretation while preserving superficial structural features. We show that this distortion represents a sparsity-dependent failure mode of current preprocessing standards and affects the comparability of datasets and computational methods across technologies and sequencing depths. To address this, we introduce a statistically consistent preprocessing framework based on nonzero-percentile clipping and log-space normalization, which preserves the intrinsic dynamic range of observed contacts. Building on this foundation, we present CCUT, a modular deep learning framework for chromatin contact map reconstruction. Under corrected preprocessing, reconstructed maps recover domain organization, contact decay, and scaling behavior consistent with polymer physics. Importantly, we demonstrate quantitative agreement between reconstructed maps and simulated contact patterns derived from a kinetic Monte Carlo loop extrusion model, enabling direct comparison between experimental data and physical models. Together, our results establish preprocessing as a decisive determinant of the physical interpretability of chromatin contact maps and provide a principled framework for robust and comparable analysis across chromatin conformation capture technologies.

Phase separated condensates are recognized as a ubiquitous mechanism of spatial organization in cell biology. Biophysical modeling of condensates provides critical insights into the dynamics and functions of these subcellular structures that are difficult to extract via experiments. Here we present an efficient computational pipeline, CASPULE (Condensate Analysis of Sticker Spacer Polymers Using the LAMMPS Engine), to simulate and analyze the biological condensates made of sticker-spacer polymers. CASPULE implements a unique force field that combines traditional Langevin dynamics with a "detailed balance proof" protocol for single-valent bond formation between stickers. This framework allows us to study the non-trivial biophysics that emerge out of the single-valent sticker interactions coupled with the effect of separation in energetic contribution by stickers and spacers. We provide detailed documentation on how to setup the simulation environment, perform simulations and analyze the results. Through case studies, we highlight the utility and efficacy of our pipeline. Importantly, we provide statistical parameters to characterize the cluster size distribution often observed in biological systems. We envision this tool to be broadly useful in decoding the interplay of kinetics and thermodynamics underlying the formation and function of biological condensates.

The cohesin complex has conserved roles in sister chromatid cohesion, DNA replication, genome organization, and the DNA damage response. We heterologously expressed the human cohesin complex in yeast to probe the behaviour of human cohesin. Human cohesin was unable to complement loss of function mutations in yeast cohesin, either as single subunits or as complexes, including in the context of co-expressing up to 12 human cohesin-associated genes. Heterologous expression of human cohesin in yeast expressing wildtype yeast cohesin resulted in dominant cohesion dysregulation and DNA damage sensitivity phenotypes. We used co-immunoprecipitation to demonstrate that human SMC proteins interact with endogenous yeast cohesin rings creating dominant-negative hybrid complexes that disrupt endogenous cohesin biology.

BackgroundGene expression changes in response to developmental and environmental cues rely on cis-regulatory sequence elements (cREs). BRG1/BRM-Associated Factors (BAF) chromatin remodeling complexes maintain chromatin accessibility at many cREs, enabling binding by transcription factors (TFs). However, cREs exhibit a broad range of sensitivity to loss of BAF function, and the basis of this variability remains unknown. ResultsTo identify the characteristics of BAF-dependent cREs, we mapped chromatin accessibility changes following acute pharmacologic BAF inhibition in GM12878 lymphoblastoid cells. We integrated these results with over 100 TF and histone modification ChIP-seq datasets and used machine learning to identify features that predict chromatin accessibility changes. We found that Activator Protein 1 (AP-1) factors and lymphoid lineage-defining TFs including RUNX3 and PU.1 predicted BAF-dependence. Strikingly, we found that cREs bearing the chromatin signature of "primed" enhancers - enriched for H3K4me1 but lacking H3K27ac - were significantly more sensitive to BAF inhibition than typical active enhancers. As primed enhancers are known to facilitate transcriptional responses to stimuli, we tested the requirement of BAF activity in these responses. Acute BAF inhibition was sufficient to prevent both chromatin and transcriptional responses to interferon gamma and dexamethasone. cREs which normally gained accessibility in response to stimulation failed to do so with BAF inhibition, and these cREs were linked to genes with suppressed transcriptional induction. ConclusionsCollectively, our results demonstrate a requirement for continuous BAF activity to enable stimulus response and suggest that defective signal responsiveness may be a pathogenic mechanism in disease states caused by loss-of-function mutations in BAF subunits.

Bacteria microcompartments (BMCs) are pseudo-organelles comprised of a self-assembling, semi-permeable protein shell, most commonly enclosing components of enzymatic pathways. -Carboxysomes (-CBs) are anabolic BMCs known for their role in sequestering Rubisco, the enzyme responsible for carbon fixation in plants, algae and bacteria, along with an upstream enzyme and an assembly protein. Rubisco has low selectivity for its substrate, CO2, and has a slow enzymatic turnover rate, resulting in an inefficient metabolic pathway. Within the -CB, Rubisco has been observed at a range of concentrations and with either a liquid-like assembly or a pseudo-lattice of polymerized fibrils. The biophysical origins of the fibril ultrastructure organization are unclear; however, it is only observed inside -CBs. Quantitative knowledge of the binding constants and energies for assembly and maintenance of these fibrils is critical for understanding this organization and Rubisco regulation, but quantitative methods for in situ analysis of Rubisco polymerization have been lacking. Here, we present an approach to convert tomography-derived -CB volumes and Rubisco particle positions into polymerization binding curves. We used this procedure to determine the Rubisco polymerization constants, including the nucleus size (n) and equilibrium polymerization constant (Kpol). The adopted modeling approach is consistent with in situ constraints, such as concentration-dependent binding interactions and confinement. This approach offers a powerful tool to evaluate both in vitro and potentially in vivo biomolecular interactions, both of Rubisco and of other proteins and polymers suitable for analysis by cryo-electron tomography. Significance StatementCryogenic electron tomography (cryoET) is a powerful method to resolve structures of proteins in their native environment at subnanometer-level resolution. Because tomography data retains spatial relationships of all particles, it intrinsically contains information about component (e.g., protein) binding interactions. Here, we use Rubisco polymerization in -carboxysomes as a model system to demonstrate that quantitative, biochemical binding analysis is possible with cryoET.

Dominantly inherited missense and gene dosage mutations in SNCA, the -synuclein gene, cause familial forms of Parkinsons disease and dementia. Here we report the structures of -synuclein filaments from the brains of such individuals. Pathogenic mutations A53T and G51D in SNCA give rise to singlets and doublets of the Lewy fold with a left-handed helical twist in the absence of a peptide-like density for island A. By contrast, filaments from the non-pathogenic variant H50Q consist of singlets of the right-handed Lewy fold with a density for island A, like filaments of wild-type -synuclein. The structures of filaments from homozygous mice transgenic for human mutant A53T -synuclein (line M83) are unlike those from human brains. They are more similar to the multiple system atrophy folds than to the Lewy fold of Parkinsons disease, Parkinsons disease dementia and dementia with Lewy bodies

In living systems across developmental and cancer biology, populations of cells on surfaces organize themselves into aggregates that mediate function and disease. Recent experimental studies have identified that such aggregates can have emergent fluid-like properties such as surface tension, yet the physical origin of these properties is not clear. Here, we develop a minimal cell-based model in which cell-cell and cell-substrate interactions are governed by active intermittent attachments. We explain the transition of cells from a dilute population to a dense aggregate, and quantify the emergent material properties underpinning this transition. We use our model to interpret experiments on dewetting of aggregates of MDA-MDB-231 cancer cells and shape fluctuations of surface-associated OVCAR3 cell aggregates. Finally, we show how spatial heterogeneity in attachments governs collective chemotaxis of cell aggregates. Together, these results reveal how active intermittent attachments generate cell aggregates with emergent material properties, with broad implications for development and cancer.

Cells dynamically integrate biochemical and mechanical signals arising from their surrounding microenvironment to regulate morphology and behavior. Mechanical cues like matrix stiffness, surface topography, and other physical perturbations modify biophysical signals. Surface topography, particularly curvature regime acts as any important mediator of mechanotransduction by coordinating cytoskeletal organization, focal adhesion dynamics, and nuclear architecture. Curvature response has been demonstrated at broader length scales and influences nucleus shape change, chromatin organization, and gene regulation, positioning the nucleus as an active mechanosensitive hub. Bone tissue consists of a curvature-rich microenvironment defined by a trabecular architecture at tissue scale and by resorption cavities such as Howships lacunae at cellular scale. While these geometries are essential for homeostasis, their role in pathological context remains poorly understood. Osteosarcoma develops within this mechanically complex multiscale architecture, but how bone-inspired curvature regulates nuclear behavior and signaling in osteosarcoma cells remains unclear. Here, we engineered three-dimensional (3D) concave hemispherical substrates that recapitulate nucleus-scale bone micro-curvature and assessed their effects on human SaOS-2 osteosarcoma cells. In comparison with flat surfaces, concave confinement resulted in pronounced nuclear rounding and softening, accompanied by Lamin A/C reorganization and increased heterochromatin compaction marked by H3K9me3. Curvature-driven nuclear remodeling selectively modulated Hippo pathway main effectors YAP/TAZ without activating NF-{kappa}B mediated canonical inflammatory responses. Furthermore, cells maintained overall viability without elevated pathological DNA damage or apoptotic signaling, suggesting an adaptive, damage-tolerant nuclear response. Overall, these findings indicate nucleus-scale curvature as a critical regulator within the bone microenvironment that governs nuclear modelling and mechanosensitive signaling in osteosarcoma cells. Incorporating physiologically relevant geometry into in vitro models establishes new insight into cancer microenvironment crosstalk and highlights nuclear interior and outer architecture as a key regulator of tumor cell behavior.

Pregnancy results in profound physiological changes driven by dynamic and precisely programmed molecular processes. Maternal peripheral blood is generally the specimen of choice for studying these processes, as it is easily accessible and essential for many aspects of maintaining a healthy pregnancy. Here, we present a high-resolution atlas of the dynamic temporal changes in the transcriptome of maternal peripheral blood in healthy human pregnancy. We generated comprehensive RNA sequencing data in 802 weekly samples from 31 healthy pregnant women from the first trimester until after delivery. Using a strict discovery and replication setup, our longitudinal analysis of gene expression identified 720 genes with robust pregnancy-specific expression patterns. Using weighted graph correlation network analysis, we identified nine pregnancy-associated transcriptional modules that reveal a strong, coordinated enrichment of innate/neutrophil and antiviral immune programs, alongside changes in adaptive immunity (T cell differentiation and signaling), erythropoiesis and hemoglobin metabolism. Cell-type deconvolution revealed that these transcriptomic shifts were accompanied by increased relative neutrophil proportions and reduced naive CD4 and CD8 T cells in pregnancy. We provide a comprehensive characterization of dynamic changes across pregnancy, highlighting maternal blood as a key systemic regulator in healthy gestation. Together, our findings establish a reference atlas of healthy pregnancy, which can be used to identify dysregulated processes and mechanisms in women with pregnancy complications. Graphical abstract O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=168 SRC="FIGDIR/small/715300v1_ufig1.gif" ALT="Figure 1"> View larger version (34K): org.highwire.dtl.DTLVardef@2a4b28org.highwire.dtl.DTLVardef@ac49d9org.highwire.dtl.DTLVardef@12468c8org.highwire.dtl.DTLVardef@15b282f_HPS_FORMAT_FIGEXP M_FIG C_FIG O_LI720 genes showed robust pregnancy specific expression patterns. C_LIO_LICo-expression analysis clustered the genes into nine modules with distinct dynamics. C_LIO_LIEnrichment in pathways involved in innate and neutrophil-mediated immunity, antiviral responses, T cell differentiation and signaling, erythropoiesis and hemoglobin metabolism. C_LIO_LICell-type deconvolution showed increases in neutrophils and decreases in naive CD4 and CD8 T cells. C_LIO_LIThe atlas of detailed longitudinal transcriptional changes provides a baseline reference for healthy pregnancy. C_LIO_LIResults for all genes and protein-protein interaction networks are made available for interactive exploration. C_LI

Spatiotemporal organisation of biological molecules is a key driver of cellular processes, including many post-transcriptional epigenetic processes. The germline-specific germ granules are biomolecular condensates that act as hubs for mRNA and small RNA processing and are core regulators of germline gene expression programming. Germ granules have been studied extensively in C. elegans, and recent developments have led to many subdivisions of the germ granule into specialised compartments. Rapid advancements in microscopy and protein-protein interaction (PPI) screening techniques have produced a large amount of data towards characterising the localisation of proteins to specific granules. However, common methods used to probe PPIs are limited in their ability to robustly detect valid interactions, especially the multivalent and sometimes transient ones observed in granule environments. Here we perform a meta-analysis of granule protein interaction screens. While these experiments generally enrich for proteins matching the profile of granule-associated proteins, we find that when considering screens individually, reproducibility is surprisingly low, highlighting not only the variability inherent in these methods but also the dynamic nature of the PPI networks present in granules. We developed an algorithm to provide a measure of each proteins association with specific granules across various experiments. By further clustering and investigation of the resulting score matrix, we demonstrate the power of this holistic approach to provide deeper insights into germ granule organisation and highlight novel can provide a resource to better inform future investigations into granules and their constituent proteins.

Background/ObjectivesRNA polymerase II is a multifunctional complex that is critical for gene regulation and environmental responses. Its POLR2I subunit in human is associated with various pathologies, including cancer chemoresistance. However, much of our understanding of how POLR2I could function indirectly derives from studies of its homologs in yeasts called Rpb9. Here, we endogenously humanized the rpb9 gene of the fission yeast Schizosaccharomyces pombe to examine the functional capabilities of POLR2I. MethodsWe edited the genomic rpb9 locus in S. pombe so that it encodes the human POLR2I protein, and investigated functional and structural conservation. ResultsWith our humanized yeast system, we find widespread functional complementation by human POLR2I of S. pombe rpb9 roles in yeast growth, chronological aging, and stress responses. We also find that POLR2I complements novel roles for yeast rpb9 in facultative heterochromatin assembly, resistance against the chemotherapy 5-fluorouracil, and resistance against the fungicide thiabendazole. In contrast, we find that POLR2I cannot complement the role of rpb9 in resistance against the transcription elongation inhibitor 6-azauracil (6-AU) in our system. Interestingly, POLR2I could complement 6-AU resistance if ectopically expressed. Lastly, we observe extensive structural homology between Rpb9 and POLR2I proteins. ConclusionsOur study establishes an endogenous cross-species gene complementation strategy that uncovers both conserved and rewired functions of fission yeast rpb9 and its human homolog, POLR2I. In addition to validating conserved roles, we also identified conservation of previously unrecognized roles of rpb9 in heterochromatin formation and chemoresistance.

Intercalation of small molecules between DNA base pairs affects DNA conformation, disrupting essential cellular processes including replication, transcription, and repair. We investigated conformational changes in 18-mer DNA upon intercalation of doxorubicin, SYBR Gold and YOYO-1 using extensive MD simulations. Two main patterns for the intercalation were identified: RISE-type intercalation occurs between adjacent base pairs and extends the DNA helix with decreased twist angles, while OPEN-type intercalation proceeds through base-pair opening without significant DNA extension. Kinetic analysis revealed that association rates for intercalation followed the order: first YO-moiety (mono-intercalation) > SYBR Gold > doxorubicin > YOYO-1 (bis-intercalation). Free energy landscape showed that forces at DNA termini reached up to 117 pN during stretching. Notably, base pairs adjacent to intercalators were protected from strand separation, accompanied by additional helical unwinding. MM-PBSA/GBSA analysis revealed that the driving force for intercalation is the stacking energy, and the binding affinity was highest for minor groove binding. Persistence length decreased with single molecule binding but recovered with two molecules due to their electrostatic repulsion. Mechanical properties of intercalated DNA showed position-dependence, demonstrating that multiple intercalation modes coexist in solution. The heterogeneous nature of intercalation explains why experimental measurements reflect ensemble averages rather than single binding configurations.

BackgroundConventional models of human ribosomal DNA (rDNA) array organization have historically depended on transcription-centric boundaries, partitioning the unit into a [~]13 kb rDNA transcription region and a monolithic [~]31 kb intergenic spacer (IGS). While our previous identification of Duplication Segment Units (DSUs) mapped these arrays based on an intuitive analysis of the microsatellite density landscape of the complete reference human genome, our present deep mining of this landscape has revealed a more accurate rDNA Gene Unit Pattern. Methods & ResultsIn this study, we conducted a deep mining analysis of our previously established microsatellite density landscape of the T2T-CHM13 assembly, focusing specifically on nucleolar organizing regions (NORs). We suggest a more accurate rDNA Gene Unit Pattern containing a (CTTT)n microsatellite aggregation ahead of the rDNA gene and a (CT)n microsatellite aggregation behind the gene, rather than a pattern featuring an IGS region inserted between two rDNA genes. ConclusionsA correct rDNA gene pattern of the human genome probably includes a (CTTT)n microsatellite aggregation ahead of the gene and a (CT)n microsatellite aggregation behind it, which possibly constitute cis- and trans-regulating regions; the (CTTT)n and (CT)n microsatellite aggregations may provide two different local stable DNA structures for regulatory protein binding.