Structure-driven homology pairing of chromatin fibers: the role of electrostatics and protein-induced bridging

Deciphering principles of nucleosome interactions and impact of cancer-associated mutations from comprehensive interaction network analysis

Nucleosomes represent hubs in chromatin organization and gene regulation and interact with a plethora of chromatin factors through different modes. In addition, alterations in histone proteins such as cancer mutations and post-translational modifications have profound effects on histone/nucleosome interactions. To elucidate the principles of histone interactions and the effects of those alterations, we developed histone interactomes for comprehensive mapping of histone-histone interactions (HHIs), histone-DNA interactions (HDIs), histone-partner interactions (HPIs) and DNA-partner interactions (DPIs) of 37 organisms, which contains a total of 3808 HPIs from 2544 binding proteins and 339 HHIs, 100 HDIs and 142 DPIs across 110 histone variants. With the developed networks, we explored histone interactions at different levels of granularities (protein-, domain- and residue-level) and performed systematic analysis on histone interactions at a large scale. Our analyses have characterized the preferred binding hotspots on both nucleosomal/linker DNA and histone octamer and unraveled diverse binding modes between nucleosome and different classes of binding partners. Last, to understand the impact of histone cancer-associated mutations on histone/nucleosome interactions, we complied one comprehensive cancer mutation dataset including 7940 cancer-associated histone mutations and further mapped those mutations onto 419,125 histone interactions at the residue level. Our quantitative analyses point to histone cancer-associated mutations' strongly disruptive effects on HHIs, HDIs and HPIs. We have further predicted 57 recurrent histone cancer mutations that have large effects on histone/nucleosome interactions and may have driver status in oncogenesis.

Spatial-Temporal Genome Regulation in Stress-Response and Cell-Fate Change

Complex functioning of the genome in the cell nucleus is controlled at different levels: (a) the DNA base sequence containing all relevant inherited information; (b) epigenetic pathways consisting of protein interactions and feedback loops; (c) the genome architecture and organization activating or suppressing genetic interactions between different parts of the genome. Most research so far has shed light on the puzzle pieces at these levels. This article, however, attempts an integrative approach to genome expression regulation incorporating these different layers. Under environmental stress or during cell development, differentiation towards specialized cell types, or to dysfunctional tumor, the cell nucleus seems to react as a whole through coordinated changes at all levels of control. This implies the need for a framework in which biological, chemical, and physical manifestations can serve as a basis for a coherent theory of gene self-organization. An international symposium held at the Biomedical Research and Study Center in Riga, Latvia, on 25 July 2022 addressed novel aspects of the abovementioned topic. The present article reviews the most recent results and conclusions of the state-of-the-art research in this multidisciplinary field of science, which were delivered and discussed by scholars at the Riga symposium.

Networks and Islands of Genome Nano-architecture and Their Potential Relevance for Radiation Biology : (A Hypothesis and Experimental Verification Hints)

The cell nucleus is a complex biological system in which simultaneous reactions and functions take place to keep the cell as an individualized, specialized system running well. The cell nucleus contains chromatin packed in various degrees of density and separated in volumes of chromosome territories and subchromosomal domains. Between the chromatin, however, there is enough "free" space for floating RNA, proteins, enzymes, ATPs, ions, water molecules, etc. which are trafficking by super- and supra-diffusion to the interaction points where they are required. It seems that this trafficking works somehow automatically and drives the system perfectly. After exposure to ionizing radiation causing DNA damage from single base damage up to chromatin double-strand breaks, the whole system "cell nucleus" responds, and repair processes are starting to recover the fully functional and intact system. In molecular biology, many individual epigenetic pathways of DNA damage response or repair of single and double-strand breaks are described. How these responses are embedded into the response of the system as a whole is often out of the focus of consideration. In this article, we want to follow the hypothesis of chromatin architecture's impact on epigenetic pathways and vice versa. Based on the assumption that chromatin acts like an "aperiodic solid state within a limited volume," functionally determined networks and local topologies ("islands") can be defined that drive the appropriate repair process at a given damage site. Experimental results of investigations of the chromatin nano-architecture and DNA repair clusters obtained by means of single-molecule localization microscopy offer hints and perspectives that may contribute to verifying the hypothesis.

Observation of nucleic acid and protein correlation in chromatin of HeLa nuclei using small-angle neutron scattering with D_{2}O-H_{2}O contrast variation

The small-angle neutron scattering (SANS) on HeLa nuclei demonstrates the bifractal nature of the chromatin structural organization. The border line between two fractal structures is detected as a crossover point at Q_{c}≈4×10^{-2}nm^{-1} in the momentum transfer dependence Q^{-D}. The use of contrast variation (D_{2}O-H_{2}O) in SANS measurements reveals clear similarity in the large scale structural organizations of nucleic acids (NA) and proteins. Both NA and protein structures have a mass fractal arrangement with the fractal dimension of D≈2.5 at scales smaller than 150 nm down to 20 nm. Both NA and proteins show a logarithmic fractal behavior with D≈3 at scales larger than 150 nm up to 6000 nm. The combined analysis of the SANS and atomic force microscopy data allows one to conclude that chromatin and its constitutes (DNA and proteins) are characterized as soft, densely packed, logarithmic fractals on the large scale and as rigid, loosely packed, mass fractals on the smaller scale. The comparison of the partial cross sections from NA and proteins with one from chromatin as a whole demonstrates spatial correlation of two chromatin's components in the range up to 900 nm. Thus chromatin in HeLa nuclei is built as the unified structure of the NA and proteins entwined through each other. Correlation between two components is lost upon scale increases toward 6000 nm. The structural features at the large scale, probably, provide nuclei with the flexibility and chromatin-free space to build supercorrelations on the distance of 10^{3} nm resembling cycle cell activity, such as an appearance of nucleoli and a DNA replication.

Heterochromatin Networks: Topology, Dynamics, and Function (a Working Hypothesis)

Open systems can only exist by self-organization as pulsing structures exchanging matter and energy with the outer world. This review is an attempt to reveal the organizational principles of the heterochromatin supra-intra-chromosomal network in terms of nonlinear thermodynamics. The accessibility of the linear information of the genetic code is regulated by constitutive heterochromatin (CHR) creating the positional information in a system of coordinates. These features include scale-free splitting-fusing of CHR with the boundary constraints of the nucleolus and nuclear envelope. The analysis of both the literature and our own data suggests a radial-concentric network as the main structural organization principle of CHR regulating transcriptional pulsing. The dynamic CHR network is likely created together with nucleolus-associated chromatin domains, while the alveoli of this network, including springy splicing speckles, are the pulsing transcription hubs. CHR contributes to this regulation due to the silencing position variegation effect, stickiness, and flexible rigidity determined by the positioning of nucleosomes. The whole system acts in concert with the elastic nuclear actomyosin network which also emerges by self-organization during the transcriptional pulsing process. We hypothesize that the the transcriptional pulsing, in turn, adjusts its frequency/amplitudes specified by topologically associating domains to the replication timing code that determines epigenetic differentiation memory.

Nucleosome-induced homology recognition in chromatin

One of the least understood properties of chromatin is the ability of its similar regions to recognize each other through weak interactions. Theories based on electrostatic interactions between helical macromolecules suggest that the ability to recognize sequence homology is an innate property of the non-ideal helical structure of DNA. However, this theory does not account for the nucleosomal packing of DNA. Can homologous DNA sequences recognize each other while wrapped up in the nucleosomes? Can structural homology arise at the level of nucleosome arrays? Here, we present a theoretical model for the recognition potential well between chromatin fibres sliding against each other. This well is different from the one predicted for bare DNA; the minima in energy do not correspond to literal juxtaposition, but are shifted by approximately half the nucleosome repeat length. The presence of this potential well suggests that nucleosome positioning may induce mutual sequence recognition between chromatin fibres and facilitate the formation of chromatin nanodomains. This has implications for nucleosome arrays enclosed between CTCF-cohesin boundaries, which may form stiffer stem-like structures instead of flexible entropically favourable loops. We also consider switches between chromatin states, e.g. through acetylation/deacetylation of histones, and discuss nucleosome-induced recognition as a precursory stage of genetic recombination.

Phylostratic Shift of Whole-Genome Duplications in Normal Mammalian Tissues towards Unicellularity Is Driven by Developmental Bivalent Genes and Reveals a Link to Cancer

Tumours were recently revealed to undergo a phylostratic and phenotypic shift to unicellularity. As well, aggressive tumours are characterized by an increased proportion of polyploid cells. In order to investigate a possible shared causation of these two features, we performed a comparative phylostratigraphic analysis of ploidy-related genes, obtained from transcriptomic data for polyploid and diploid human and mouse tissues using pairwise cross-species transcriptome comparison and principal component analysis. Our results indicate that polyploidy shifts the evolutionary age balance of the expressed genes from the late metazoan phylostrata towards the upregulation of unicellular and early metazoan phylostrata. The up-regulation of unicellular metabolic and drug-resistance pathways and the downregulation of pathways related to circadian clock were identified. This evolutionary shift was associated with the enrichment of ploidy with bivalent genes (p < 10⁻¹⁶). The protein interactome of activated bivalent genes revealed the increase of the connectivity of unicellulars and (early) multicellulars, while circadian regulators were depressed. The mutual polyploidy-c-MYC-bivalent genes-associated protein network was organized by gene-hubs engaged in both embryonic development and metastatic cancer including driver (proto)-oncogenes of viral origin. Our data suggest that, in cancer, the atavistic shift goes hand-in-hand with polyploidy and is driven by epigenetic mechanisms impinging on development-related bivalent genes.

Amino Acid Sequence of Oligopeptide Causes Marked Difference in DNA Compaction and Transcription

Compaction of T4 phage DNA (166 kbp) by short oligopeptide octamers composed of two types of amino acids, four cationic lysine (K), and four polar nonionic serine (S) having different sequence order was studied by single-molecule fluorescent microscopy. We found that efficient DNA compaction by oligopeptide octamers depends on the geometrical match between phosphate groups of DNA and cationic amines. The amino acid sequence order in octamers dramatically affects the mechanism of DNA compaction, which changes from a discrete all-or-nothing coil-globule transition induced by a less efficient (K4S4) octamer to a continuous compaction transition induced by a (KS)4 octamer with a stronger DNA-binding character. This difference in the DNA compaction mechanism dramatically changes the packaging density, and the morphology of T4 DNA condensates: DNA is folded into ordered toroidal or rod morphologies during all-or-nothing compaction, whereas disordered DNA condensates are formed as a result of the continuous DNA compaction. Furthermore, the difference in DNA compaction mechanism has a certain effect on the inhibition scenario of the DNA transcription activity, which is gradual for the continuous DNA compaction and abrupt for the all-or-nothing DNA collapse.

Potential Role of Phase Separation of Repetitive DNA in Chromosomal Organization

The basic principles of chromosomal organization in eukaryotic cells remain elusive. Current mainstream research efforts largely concentrate on searching for critical packaging proteins involved in organizing chromosomes. I have taken a different perspective, by considering the role of genomic information in chromatins. In particular, I put forward the concept that repetitive DNA elements are key chromosomal packaging modules, and their intrinsic property of homology-based interaction can drive chromatin folding. Many repetitive DNA families have high copy numbers and clustered distribution patterns in the linear genomes. These features may facilitate the interactions among members in the same repeat families. In this paper, the potential liquid–liquid phase transition of repetitive DNAs that is induced by their extensive interaction in chromosomes will be considered. I propose that the interaction among repetitive DNAs may lead to phase separation of interacting repetitive DNAs from bulk chromatins. Phase separation of repetitive DNA may provide a physical mechanism that drives rapid massive changes of chromosomal conformation.

Concentration-Driven Fascinating Vesicle-Fibril Transition Employing Merocyanine 540 and 1-Octyl-3-methylimidazolium Chloride

In this article, anionic lipophilic dye merocyanine 540(MC540) and cationic surface-active ionic liquid (SAIL) 1-octyl-3-methylimidazolium chloride (C₈mimCl) are employed to construct highly ordered fibrillar and vesicular aggregates exploiting an ionic self-assembly (ISA) strategy. It is noteworthy that the concentration of the counterions has exquisite control over the morphology, in which lowering the concentration of both the building blocks in a stoichiometric ratio of 1:1 provides a vesicle to fibril transition. Here, we have reported the concentration-controlled fibril-vesicle transition utilizing the emerging fluorescence lifetime imaging microscopy (FLIM) technique. Furthermore, we have detected this morphological transformation by means of other microscopic techniques such as field emission scanning electron microscopy (FESEM), transmission electron microscopy (TEM), and cryogenic-transmission electron microscopy (cryo-TEM) to gain additional support. Besides, multiwavelength FLIM (MW-FLIM) and atomic force microscopy (AFM) techniques assist us in knowing the microheterogeneity and the height profile of the vesicles, respectively. We have replaced the SAIL, C₈mimCl, by an analogous traditional surfactant, n-octyltrimethylammonium bromide (OTAB), and it provides a discernible change in morphology similar to that of C₈mimCl, whereas 1-octanol is unable to exhibit any structural aggregation and thus reveals the importance of electrostatic interaction in supramolecular aggregate formation. However, the SAILs having the same imidazolium headgroup with different chain lengths other than C₈mimCl are unable to display any structural transition and determine the importance of the correct chain length for efficient packing of the counterions to form a specific self-assembly. Therefore, this study reveals the synergistic interplay of electrostatic, hydrophobic, and π-π stacking interactions to construct the self-assembly and their concentration-dependent morphological transition.

New Evidence for the Theory of Chromosome Organization by Repetitive Elements (CORE)

Repetitive DNA elements were proposed to coordinate chromatin folding and interaction in chromosomes by their intrinsic homology-based clustering ability. A recent analysis of the data sets from chromosome-conformation-capture experiments confirms the spatial clustering of DNA repeats of the same family in the nuclear space, and thus provides strong new support for the CORE theory.

The R-Operon: A Model of Repetitive DNA-Organized Transcriptional Compartmentation of Eukaryotic Chromosomes for Coordinated Gene Expression

In eukaryotic genomes, it is essential to coordinate the activity of genes that function together to fulfill the same biological processes. Genomic organization likely plays a key role in coordinating transcription of different genes. However, little is known about how co-regulated genes are organized in the cell nucleus and how the chromosomal organization facilitates the co-regulation of different genes. I propose that eukaryotic genomes are organized into repeat assembly (RA)-based structural domains (“R-operons”) in the nuclear space. R-operons result from the interaction of homologous DNA repeats. In an R-operon, genes in different loci of the linear genome are brought into spatial vicinity and co-regulated by the same pool of transcription factors. This type of large-scale chromosomal organization may provide a mechanism for functional compartmentation of chromosomes to facilitate the transcriptional coordination of gene expression.

A computational study of dsDNA pairs and vibrational resonance in separating water

This article investigates the relationship between molecular sequence and dependent interacting behavior of molecular segment pairs and secondly, sequence dependent, vibrational resonance in surrounding normal saline, protein-free water. The development of a molecular model to explore these systems phenomena, the results of several nanoscale molecular dynamics simulations, and analysis of behavior of interacting ΦX174 double-stranded DNA segment pair models in various configurations are presented. Fourier analysis revealed intriguing vibration frequencies within the solvent plane between the segments, while subsequent frequency domain transformation of the time domain waveforms revealed statistically significant resonating harmonic signals in the THz range.

Identification of peculiar and common effects of histone modifications on transcription

Histone modifications (HMs) play an important role in controlling eukaryotic gene expression and next generation sequencing (NGS) has greatly advanced the research on this topic with generating many high-resolution maps for HMs. Here, we use these maps to analyze the relationship between HMs and transcription. By incorporating various segments of genes into analysis without restricting the scope only in the promoter region, we have collected more comprehensive data and captured some details of this process. A position effect of gene regions has been revealed and it can even inverse the property of some HMs from activating to repressing genes such as the cases of H3K4me3, H3K36me3 and H3K14ac. Especially H3K36me3, its dual character on gene transcription makes it able to serve as a criterion to distinguish high and low expressed genes. We also study the general property of different HMs based on the comprehensive data. Using exploratory factor analysis (EFA), we have extracted 4 latent structures underlying the HMs, which are able to represent their activating and repressing effects concisely. These 4 factors have fine properties in the aspects of distinguishing high and low expressed genes, predicting transcription level and identifying genes with unique attributes such as stable RNA generating genes found to have a close relationship with lifespan of organisms here. In summary, while the position effect associated HM peculiarities demonstrate some details of the complex HM regulation network divergently, the common factors catch the nature of the network convergently. This deepens our understanding on the HM-transcription relationship.

The 4D nucleome: Evidence for a dynamic nuclear landscape based on co-aligned active and inactive nuclear compartments

Recent methodological advancements in microscopy and DNA sequencing-based methods provide unprecedented new insights into the spatio-temporal relationships between chromatin and nuclear machineries. We discuss a model of the underlying functional nuclear organization derived mostly from electron and super-resolved fluorescence microscopy studies. It is based on two spatially co-aligned, active and inactive nuclear compartments (ANC and INC). The INC comprises the compact, transcriptionally inactive core of chromatin domain clusters (CDCs). The ANC is formed by the transcriptionally active periphery of CDCs, called the perichromatin region (PR), and the interchromatin compartment (IC). The IC is connected to nuclear pores and serves nuclear import and export functions. The ANC is the major site of RNA synthesis. It is highly enriched in epigenetic marks for transcriptionally competent chromatin and RNA Polymerase II. Marks for silent chromatin are enriched in the INC. Multi-scale cross-correlation spectroscopy suggests that nuclear architecture resembles a random obstacle network for diffusing proteins. An increased dwell time of proteins and protein complexes within the ANC may help to limit genome scanning by factors or factor complexes to DNA exposed within the ANC.

Thermodynamic model of heterochromatin formation through epigenetic regulation

Gene regulation in eukaryotes requires the segregation of silenced genomic regions into densely packed heterochromatin, leaving the active genes in euchromatin regions more accessible. We introduce a model that connects the presence of epigenetically inherited histone marks, methylation at histone 3 lysine-9, to the physical compaction of chromatin fibers via the binding of heterochromatin protein 1 (HP1). Our model demonstrates some of the key physical features that are necessary to explain experimental observations. In particular, we demonstrate that strong cooperative interactions among the HP1 proteins are necessary to see the phase segregation of heterochromatin and euchromatin regions. We also explore how the cell can use the concentration of HP1 to control condensation and under what circumstances there is a threshold of methylation over which the fibers will compact. Finally, we consider how different potential in vivo fiber structures as well as the flexibility of the histone 3 tail can affect the bridging of HP1. Many of the observations that we make about the HP1 system are guided by general thermodynamics principles and thus could play a role in other DNA organizational processes such as the binding of linker histones.

Kinetics of polymer looping with macromolecular crowding: effects of volume fraction and crowder size

The looping of polymers such as DNA is a fundamental process in the molecular biology of living cells, whose interior is characterised by a high degree of molecular crowding. We here investigate in detail the looping dynamics of flexible polymer chains in the presence of different degrees of crowding. From the analysis of the looping-unlooping rates and the looping probabilities of the chain ends we show that the presence of small crowders typically slows down the chain dynamics but larger crowders may in fact facilitate the looping. We rationalise these non-trivial and often counterintuitive effects of the crowder size on the looping kinetics in terms of an effective solution viscosity and standard excluded volume. It is shown that for small crowders the effect of an increased viscosity dominates, while for big crowders we argue that confinement effects (caging) prevail. The tradeoff between both trends can thus result in the impediment or facilitation of polymer looping, depending on the crowder size. We also examine how the crowding volume fraction, chain length, and the attraction strength of the contact groups of the polymer chain affect the looping kinetics and hairpin formation dynamics. Our results are relevant for DNA looping in the absence and presence of protein mediation, DNA hairpin formation, RNA folding, and the folding of polypeptide chains under biologically relevant high-crowding conditions.

Electrostatic effect of H1-histone protein binding on nucleosome repeat length

Within a simple biophysical model we describe the effect of electrostatic binding of H1 histone proteins on the nucleosome repeat length in chromatin. The length of wrapped DNA optimizes its binding energy to the histone core and the elastic energy penalty of DNA wrapping. The magnitude of the effect predicted from our model is in agreement with the systematic experimental data on the linear variation of nucleosome repeat lengths with H1/nucleosome ratio (Woodcock C L et al 2006 Chromos. Res. 14 17-25). We compare our model to the data for different cell types and organisms, with a widely varying ratio of bound H1 histones per nucleosome. We underline the importance of this non-specific histone-DNA charge-balance mechanism in regulating the positioning of nucleosomes and the degree of compaction of chromatin fibers in eukaryotic cells.

Regulation of the nucleosome repeat length in vivo by the DNA sequence, protein concentrations and long-range interactions

The nucleosome repeat length (NRL) is an integral chromatin property important for its biological functions. Recent experiments revealed several conflicting trends of the NRL dependence on the concentrations of histones and other architectural chromatin proteins, both in vitro and in vivo, but a systematic theoretical description of NRL as a function of DNA sequence and epigenetic determinants is currently lacking. To address this problem, we have performed an integrative biophysical and bioinformatics analysis in species ranging from yeast to frog to mouse where NRL was studied as a function of various parameters. We show that in simple eukaryotes such as yeast, a lower limit for the NRL value exists, determined by internucleosome interactions and remodeler action. For higher eukaryotes, also the upper limit exists since NRL is an increasing but saturating function of the linker histone concentration. Counterintuitively, smaller H1 variants or non-histone architectural proteins can initiate larger effects on the NRL due to entropic reasons. Furthermore, we demonstrate that different regimes of the NRL dependence on histone concentrations exist depending on whether DNA sequence-specific effects dominate over boundary effects or vice versa. We consider several classes of genomic regions with apparently different regimes of the NRL variation. As one extreme, our analysis reveals that the period of oscillations of the nucleosome density around bound RNA polymerase coincides with the period of oscillations of positioning sites of the corresponding DNA sequence. At another extreme, we show that although mouse major satellite repeats intrinsically encode well-defined nucleosome preferences, they have no unique nucleosome arrangement and can undergo a switch between two distinct types of nucleosome positioning.