Comparison of cell response to chromatin and DNA damage.

DNA-targeting drugs are widely used for anti-cancer treatment. Many of these drugs cause different types of DNA damage, i.e. alterations in the chemical structure of DNA molecule. However, molecules binding to DNA may also interfere with DNA packing into chromatin. Interestingly, some molecules do not cause any changes in DNA chemical structure but interfere with DNA binding to histones and nucleosome wrapping. This results in histone loss from chromatin and destabilization of nucleosomes, a phenomenon that we call chromatin damage. Although the cellular response to DNA damage is well-studied, the consequences of chromatin damage are not. Moreover, many drugs used to study DNA damage also cause chromatin damage, therefore there is no clarity on which effects are caused by DNA or chromatin damage. In this study, we aimed to clarify this issue. We treated normal and tumor cells with bleomycin, nuclease mimicking drug which cut predominantly nucleosome-free DNA and therefore causes DNA damage in the form of DNA breaks, and CBL0137, which causes chromatin damage without direct DNA damage. We describe similarities and differences between the consequences of DNA and chromatin damage. Both agents were more toxic for tumor than normal cells, but while DNA damage causes senescence in both normal and tumor cells, chromatin damage does not. Both agents activated p53, but chromatin damage leads to the accumulation of higher levels of unmodified p53, which transcriptional activity was similar to or lower than that of p53 activated by DNA damage. Most importantly, we found that while transcriptional changes caused by DNA damage are limited by p53-dependent activation of a small number of p53 targets, chromatin damage activated many folds more genes in p53 independent manner.

Comparison of cell response to chromatin and DNA damage.

DNA-targeting drugs may damage DNA or chromatin. Many anti-cancer drugs damage both, making it difficult to understand their mechanisms of action. Using molecules causing DNA breaks without altering nucleosome structure (bleomycin) or destabilizing nucleosomes without damaging DNA (curaxin), we investigated the consequences of DNA or chromatin damage in normal and tumor cells. As expected, DNA damage caused p53-dependent growth arrest followed by senescence. Chromatin damage caused higher p53 accumulation than DNA damage; however, growth arrest was p53-independent and did not result in senescence. Chromatin damage activated the transcription of multiple genes, including classical p53 targets, in a p53-independent manner. Although these genes were not highly expressed in basal conditions, they had chromatin organization around the transcription start sites (TSS) characteristic of most highly expressed genes and the highest level of paused RNA polymerase. We hypothesized that nucleosomes around the TSS of these genes were the most sensitive to chromatin damage. Therefore, nucleosome loss upon curaxin treatment would enable transcription without the assistance of sequence-specific transcription factors. We confirmed this hypothesis by showing greater nucleosome loss around the TSS of these genes upon curaxin treatment and activation of a p53-specific reporter in p53-null cells by chromatin-damaging agents but not DNA-damaging agents.

A multi-reference poly-conformational method for in silico design, optimization, and repositioning of pharmaceutical compounds illustrated for selected SARS-CoV-2 ligands.

Background: This work presents a novel computational multi-reference poly-conformational algorithm for design, optimization, and repositioning of pharmaceutical compounds. Methods: The algorithm searches for candidates by comparing similarities between conformers of the same compound and identifies target compounds, whose conformers are collectively close to the conformers of each compound in the reference set. Reference compounds may possess highly variable MoAs, which directly, and simultaneously, shape the properties of target candidate compounds. Results: The algorithm functionality has been case study validated in silico, by scoring ChEMBL drugs against FDA-approved reference compounds that either have the highest predicted binding affinity to our chosen SARS-CoV-2 targets or are confirmed to be inhibiting such targets in-vivo. All our top scoring ChEMBL compounds also turned out to be either high-affinity ligands to the chosen targets (as confirmed in separate studies) or show significant efficacy, in-vivo, against those selected targets. In addition to method case study validation, in silico search for new compounds within two virtual libraries from the Enamine database is presented. The library's virtual compounds have been compared to the same set of reference drugs that we used for case study validation: Olaparib, Tadalafil, Ergotamine and Remdesivir. The large reference set of four potential SARS-CoV-2 compounds has been selected, since no drug has been identified to be 100% effective against the virus so far, possibly because each candidate drug was targeting only one, particular MoA. The goal here was to introduce a new methodology for identifying potential candidate(s) that cover multiple MoA-s presented within a set of reference compounds.

Novel Efficient Multistage Lead Optimization Pipeline Experimentally Validated for DYRK1B Selective Inhibitors.

In addition to general challenges in drug discovery such as the identification of lead compounds in time- and cost-effective ways, specific challenges also exist. Particularly, it is necessary to develop pharmacological inhibitors that effectively discriminate between closely related molecular targets. DYRK1B kinase is considered a valuable target for cancer-specific mono- or combination chemotherapy; however, the inhibition of its closely related DYRK1A kinase is not beneficial. Existing inhibitors target both kinases with essentially the same efficiency, and the unavailability of the DYRK1B crystal structure makes the discovery of DYRK1B-specific inhibitors even more challenging. Here, we propose a novel multi-stage compound discovery pipeline aimed at in silico identification of both potent and selective small molecules from a large set of initial candidates. The method uses structure-based docking and ligand-based quantitative structure-activity relationship modeling. This approach allowed us to identify lead and runner-up small-molecule compounds targeting DYRK1B with high efficiency and specificity.

Domain Model of Eukaryotic Genome Organization: From DNA Loops Fixed on the Nuclear Matrix to TADs.

The article reviews the development of ideas on the domain organization of eukaryotic genome, with special attention on the studies of DNA loops anchored to the nuclear matrix and their role in the emergence of the modern model of eukaryotic genome spatial organization. Critical analysis of results demonstrating that topologically associated chromatin domains are structural-functional blocks of the genome supports the notion that these blocks are fundamentally different from domains whose existence was proposed by the domain hypothesis of eukaryotic genome organization formulated in the 1980s. Based on the discussed evidence, it is concluded that the model postulating that eukaryotic genome is built from uniformly organized structural-functional blocks has proven to be untenable.

The twisted path of the 3D genome: where does it lead?

A new era in 3D genome studies began with the development of the so-called 'C-methods', used for the analysis of spatial contacts between distant genomic elements. However, the idea that spatial genome organization, partitioning of the genome into structural/functional units, and the functional compartmentalization of the cell nucleus are important for the implementation of key functions of the genome arose much earlier. In this Opinion article, we briefly overview how the concept of spatial genome organization has changed over recent decades, discuss current views on the 3D genome and cell nucleus organization, and compare the experimental evidence for the inter-relation between gene regulation and the 3D genome.

DNA G-Quadruplexes Contribute to CTCF Recruitment.

G-quadruplex (G4) sites in the human genome frequently colocalize with CCCTC-binding factor (CTCF)-bound sites in CpG islands (CGIs). We aimed to clarify the role of G4s in CTCF positioning. Molecular modeling data suggested direct interactions, so we performed in vitro binding assays with quadruplex-forming sequences from CGIs in the human genome. G4s bound CTCF with Kd values similar to that of the control duplex, while respective i-motifs exhibited no affinity for CTCF. Using ChIP-qPCR assays, we showed that G4-stabilizing ligands enhance CTCF occupancy at a G4-prone site in STAT3 gene. In view of the reportedly increased CTCF affinity for hypomethylated DNA, we next questioned whether G4s also facilitate CTCF recruitment to CGIs via protecting CpG sites from methylation. Bioinformatics analysis of previously published data argued against such a possibility. Finally, we questioned whether G4s facilitate CTCF recruitment by affecting chromatin structure. We showed that three architectural chromatin proteins of the high mobility group colocalize with G4s in the genome and recognize parallel-stranded or mixed-topology G4s in vitro. One of such proteins, HMGN3, contributes to the association between G4s and CTCF according to our bioinformatics analysis. These findings support both direct and indirect roles of G4s in CTCF recruitment.

Co-Regulated Genes and Gene Clusters.

There are many co-regulated genes in eukaryotic cells. The coordinated activation or repression of such genes occurs at specific stages of differentiation, or under the influence of external stimuli. As a rule, co-regulated genes are dispersed in the genome. However, there are also gene clusters, which contain paralogous genes that encode proteins with similar functions. In this aspect, they differ significantly from bacterial operons containing functionally linked genes that are not paralogs. In this review, we discuss the reasons for the existence of gene clusters in vertebrate cells and propose that clustering is necessary to ensure the possibility of selective activation of one of several similar genes.

Treacle and TOPBP1 control replication stress response in the nucleolus.

Replication stress is one of the main sources of genome instability. Although the replication stress response in eukaryotic cells has been extensively studied, almost nothing is known about the replication stress response in nucleoli. Here, we demonstrate that initial replication stress-response factors, such as RPA, TOPBP1, and ATR, are recruited inside the nucleolus in response to drug-induced replication stress. The role of TOPBP1 goes beyond the typical replication stress response; it interacts with the low-complexity nucleolar protein Treacle (also referred to as TCOF1) and forms large Treacle-TOPBP1 foci inside the nucleolus. In response to replication stress, Treacle and TOPBP1 facilitate ATR signaling at stalled replication forks, reinforce ATR-mediated checkpoint activation inside the nucleolus, and promote the recruitment of downstream replication stress response proteins inside the nucleolus without forming nucleolar caps. Characterization of the Treacle-TOPBP1 interaction mode leads us to propose that these factors can form a molecular platform for efficient stress response in the nucleolus.

Suppression of liquid-liquid phase separation by 1,6-hexanediol partially compromises the 3D genome organization in living cells.

Liquid-liquid phase separation (LLPS) contributes to the spatial and functional segregation of molecular processes within the cell nucleus. However, the role played by LLPS in chromatin folding in living cells remains unclear. Here, using stochastic optical reconstruction microscopy (STORM) and Hi-C techniques, we studied the effects of 1,6-hexanediol (1,6-HD)-mediated LLPS disruption/modulation on higher-order chromatin organization in living cells. We found that 1,6-HD treatment caused the enlargement of nucleosome clutches and their more uniform distribution in the nuclear space. At a megabase-scale, chromatin underwent moderate but irreversible perturbations that resulted in the partial mixing of A and B compartments. The removal of 1,6-HD from the culture medium did not allow chromatin to acquire initial configurations, and resulted in more compact repressed chromatin than in untreated cells. 1,6-HD treatment also weakened enhancer-promoter interactions and TAD insulation but did not considerably affect CTCF-dependent loops. Our results suggest that 1,6-HD-sensitive LLPS plays a limited role in chromatin spatial organization by constraining its folding patterns and facilitating compartmentalization at different levels.

LASCA: loop and significant contact annotation pipeline.

Chromatin loops represent one of the major levels of hierarchical folding of the genome. Although the situation is evolving, current methods have various difficulties with the accurate mapping of loops even in mammalian Hi-C data, and most of them fail to identify chromatin loops in animal species with substantially different genome architecture. This paper presents the loop and significant contact annotation (LASCA) pipeline, which uses Weibull distribution-based modeling to effectively identify loops and enhancer-promoter interactions in Hi-C data from evolutionarily distant species: from yeast and worms to mammals. Available at: https://github.com/ArtemLuzhin/LASCA_pipeline .

Genetically Encoded Fluorescent Sensor for Poly-ADP-Ribose.

Poly-(ADP-ribosyl)-ation (PARylation) is a reversible post-translational modification of proteins and DNA that plays an important role in various cellular processes such as DNA damage response, replication, transcription, and cell death. Here we designed a fully genetically encoded fluorescent sensor for poly-(ADP-ribose) (PAR) based on Förster resonance energy transfer (FRET). The WWE domain, which recognizes iso-ADP-ribose internal PAR-specific structural unit, was used as a PAR-targeting module. The sensor consisted of cyan Turquoise2 and yellow Venus fluorescent proteins, each in fusion with the WWE domain of RNF146 E3 ubiquitin ligase protein. This bipartite sensor named sPARroW (sensor for PAR relying on WWE) enabled monitoring of PAR accumulation and depletion in live mammalian cells in response to different stimuli, namely hydrogen peroxide treatment, UV irradiation and hyperthermia.

Chromatin Trapping of Factors Involved in DNA Replication and Repair Underlies Heat-Induced Radio- and Chemosensitization.

Hyperthermia has been used as an adjuvant treatment for radio- and chemotherapy for decades. In addition to its effects on perfusion and oxygenation of cancer tissues, hyperthermia can enhance the efficacy of DNA-damaging treatments such as radiotherapy and chemotherapy. Although it is believed that the adjuvant effects are based on hyperthermia-induced dysfunction of DNA repair systems, the mechanisms of these dysfunctions remain elusive. Here, we propose that elevated temperatures can induce chromatin trapping (c-trapping) of essential factors, particularly those involved in DNA repair, and thus enhance the sensitization of cancer cells to DNA-damaging therapeutics. Using mass spectrometry-based proteomics, we identified proteins that could potentially undergo c-trapping in response to hyperthermia. Functional analyses of several identified factors involved in DNA repair demonstrated that c-trapping could indeed be a mechanism of hyperthermia-induced transient deficiency of DNA repair systems. Based on our proteomics data, we showed for the first time that hyperthermia could inhibit maturation of Okazaki fragments and activate a corresponding poly(ADP-ribose) polymerase-dependent DNA damage response. Together, our data suggest that chromatin trapping of factors involved in DNA repair and replication contributes to heat-induced radio- and chemosensitization.

Studying RNA-DNA interactome by Red-C identifies noncoding RNAs associated with various chromatin types and reveals transcription dynamics.

Non-coding RNAs (ncRNAs) participate in various biological processes, including regulating transcription and sustaining genome 3D organization. Here, we present a method termed Red-C that exploits proximity ligation to identify contacts with the genome for all RNA molecules present in the nucleus. Using Red-C, we uncovered the RNA-DNA interactome of human K562 cells and identified hundreds of ncRNAs enriched in active or repressed chromatin, including previously undescribed RNAs. Analysis of the RNA-DNA interactome also allowed us to trace the kinetics of messenger RNA production. Our data support the model of co-transcriptional intron splicing, but not the hypothesis of the circularization of actively transcribed genes.

Weak interactions in higher-order chromatin organization.

The detailed principles of the hierarchical folding of eukaryotic chromosomes have been revealed during the last two decades. Along with structures composing three-dimensional (3D) genome organization (chromatin compartments, topologically associating domains, chromatin loops, etc.), the molecular mechanisms that are involved in their establishment and maintenance have been characterized. Generally, protein-protein and protein-DNA interactions underlie the spatial genome organization in eukaryotes. However, it is becoming increasingly evident that weak interactions, which exist in biological systems, also contribute to the 3D genome. Here, we provide a snapshot of our current understanding of the role of the weak interactions in the establishment and maintenance of the 3D genome organization. We discuss how weak biological forces, such as entropic forces operating in crowded solutions, electrostatic interactions of the biomolecules, liquid-liquid phase separation, DNA supercoiling, and RNA environment participate in chromosome segregation into structural and functional units and drive intranuclear functional compartmentalization.

The 3D Genome as a Target for Anticancer Therapy.

The role of 3D genome organization in the precise regulation of gene expression is well established. Accordingly, the mechanistic connections between 3D genome alterations and disease development are becoming increasingly apparent. This opinion article provides a snapshot of our current understanding of the 3D genome alterations associated with cancers. We discuss potential connections of the 3D genome and cancer transcriptional addiction phenomenon as well as molecular mechanisms of action of 3D genome-disrupting drugs. Finally, we highlight issues and perspectives raised by the discovery of the first pharmaceutical strongly affecting 3D genome organization.

C-TALE, a new cost-effective method for targeted enrichment of Hi-C/3C-seq libraries.

Studies performed using Hi-C and other high-throughput whole-genome C-methods have demonstrated that 3D organization of eukaryotic genomes is functionally relevant. Unfortunately, ultra-deep sequencing of Hi-C libraries necessary to detect loop structures in large vertebrate genomes remains rather expensive. However, many studies are in fact aimed at determining the fine-scale 3D structure of comparatively small genomic regions up to several Mb in length. Such studies typically focus on the spatial structure of domains of coregulated genes, molecular mechanisms of loop formation, and interrogation of functional significance of GWAS-revealed polymorphisms. Therefore, a handful of molecular techniques based on Hi-C have been developed to address such issues. These techniques commonly rely on in-solution hybridization of Hi-C/3C-seq libraries with pools of biotinylated baits covering the region of interest, followed by deep sequencing of the enriched library. Here, we describe a new protocol of this kind, C-TALE (Chromatin TArget Ligation Enrichment). Preparation of hybridization probes from bacterial artificial chromosomes and an additional round of enrichment make C-TALE a cost-effective alternative to existing many-versus-all C-methods.

Hypoosmotic stress induces R loop formation in nucleoli and ATR/ATM-dependent silencing of nucleolar transcription.

The contribution of nucleoli to the cellular stress response has been discussed for over a decade. Stress-induced inhibition of RNA polymerase I-dependent transcription is hypothesized as a possible effector program in such a response. In this study, we report a new mechanism by which ribosomal DNA transcription can be inhibited in response to cellular stress. Specifically, we demonstrate that mild hypoosmotic stress induces stabilization of R loops in ribosomal genes and thus provokes the nucleoli-specific DNA damage response, which is governed by the ATM- and Rad3-related (ATR) kinase. Activation of ATR in nucleoli strongly depends on Treacle, which is needed for efficient recruitment/retention of TopBP1 in nucleoli. Subsequent ATR-mediated activation of ATM results in repression of nucleolar transcription.

The anti-cancer drugs curaxins target spatial genome organization.

Recently we characterized a class of anti-cancer agents (curaxins) that disturbs DNA/histone interactions within nucleosomes. Here, using a combination of genomic and in vitro approaches, we demonstrate that curaxins strongly affect spatial genome organization and compromise enhancer-promoter communication, which is necessary for the expression of several oncogenes, including MYC. We further show that curaxins selectively inhibit enhancer-regulated transcription of chromatinized templates in cell-free conditions. Genomic studies also suggest that curaxins induce partial depletion of CTCF from its binding sites, which contributes to the observed changes in genome topology. Thus, curaxins can be classified as epigenetic drugs that target the 3D genome organization.

Synthetically Lethal Interactions of ATM, ATR, and DNA-PKcs.

Synthetic lethality occurs when simultaneous perturbations of two genes or molecular processes result in a loss of cell viability. The number of known synthetically lethal interactions is growing steadily. We review here synthetically lethal interactions of ataxia-telangiectasia mutated (ATM), ATM- and Rad3-related (ATR), and DNA-dependent protein kinase catalytic subunit (DNA-PKcs). These kinases are appropriate for synthetic lethal therapies because their genes are frequently mutated in cancer, and specific inhibitors are currently in clinical trials. Understanding synthetically lethal interactions of a particular gene or gene family can facilitate predicting new synthetically lethal interactions, therapy toxicity, and mechanisms of resistance, as well as defining the spectrum of tumors amenable to these therapeutic approaches.