The role of nucleotide opening dynamics in facilitated target search by DNA-repair proteins.

Preserving the genomic integrity stands a fundamental necessity, primarily achieved by the DNA repair proteins through their continuous patrolling on the DNA in search of lesions. However, comprehending how even a single base-pair lesion can be swiftly and specifically recognized amidst millions of base-pair sites remains a formidable challenge. In this study, we employ extensive molecular dynamics simulations using an appropriately tuned model of both protein and DNA to probe the underlying molecular principles. Our findings reveal that the dynamics of a non-canonical base generate an entropic signal that guides the one-dimensional search of a repair protein, thereby facilitating the recognition of the lesion site. The width of the funnel perfectly aligns with the one-dimensional diffusion length of DNA-binding proteins. The generic mechanism provides a physical basis for rapid recognition and specificity of DNA damage sensing and recognition.

The environmentally-regulated interplay between local three-dimensional chromatin organisation and transcription of proVWX in E. coli.

Nucleoid associated proteins (NAPs) maintain the architecture of bacterial chromosomes and regulate gene expression. Thus, their role as transcription factors may involve three-dimensional chromosome re-organisation. While this model is supported by in vitro studies, direct in vivo evidence is lacking. Here, we use RT-qPCR and 3C-qPCR to study the transcriptional and architectural profiles of the H-NS (histone-like nucleoid structuring protein)-regulated, osmoresponsive proVWX operon of Escherichia coli at different osmolarities and provide in vivo evidence for transcription regulation by NAP-mediated chromosome re-modelling in bacteria. By consolidating our in vivo investigations with earlier in vitro and in silico studies that provide mechanistic details of how H-NS re-models DNA in response to osmolarity, we report that activation of proVWX in response to a hyperosmotic shock involves the destabilization of H-NS-mediated bridges anchored between the proVWX downstream and upstream regulatory elements (DRE and URE), and between the DRE and ygaY that lies immediately downstream of proVWX. The re-establishment of these bridges upon adaptation to hyperosmolarity represses the operon. Our results also reveal additional structural features associated with changes in proVWX transcript levels such as the decompaction of local chromatin upstream of the operon, highlighting that further complexity underlies the regulation of this model operon. H-NS and H-NS-like proteins are wide-spread amongst bacteria, suggesting that chromosome re-modelling may be a typical feature of transcriptional control in bacteria.

Superstructure Detection in Nucleosome Distribution Shows Common Pattern within a Chromosome and within the Genome.

Nucleosome positioning plays an important role in crucial biological processes such as replication, transcription, and gene regulation. It has been widely used to predict the genome's function and chromatin organisation. So far, the studies of patterns in nucleosome positioning have been limited to transcription start sites, CTCFs binding sites, and some promoter and loci regions. The genome-wide organisational pattern remains unknown. We have developed a theoretical model to coarse-grain nucleosome positioning data in order to obtain patterns in their distribution. Using hierarchical clustering on the auto-correlation function of this coarse-grained nucleosome positioning data, a genome-wide clustering is obtained for Candida albicans. The clustering shows the existence beyond hetero- and eu-chromatin inside the chromosomes. These non-trivial clusterings correspond to different nucleosome distributions and gene densities governing differential gene expression patterns. Moreover, these distribution patterns inside the chromosome appeared to be conserved throughout the genome and within species. The pipeline of the coarse grain nucleosome positioning sequence to identify underlying genomic organisation used in our study is novel, and the classifications obtained are unique and consistent.

Inter-nucleosomal potentials from nucleosomal positioning data.

No systematic method exists to derive inter-nucleosomal potentials between nucleosomes along a chromosome consistently across a given genome. Such potentials can yield information on nucleosomal ordering, thermal as well as mechanical properties of chromosomes. Thus, indirectly, they shed light on a possible mechanical genomic code along a chromosome. To develop a method yielding effective inter-nucleosomal potentials between nucleosomes, a generalized Lennard-Jones potential for the parameterization is developed based on nucleosomal positioning data. This approach eliminates some of the problems that the underlying nucleosomal positioning data have, rendering the extraction difficult on the individual nucleosomal level. Furthermore, patterns on which to base a classification along a chromosome appear on larger domains, such as hetero- and euchromatin. An intuitive selection strategy for the noisy optimization problem is employed to derive effective exponents for the generalized potential. The method is tested on the Candida albicans genome. Applying k-means clustering based on potential parameters and thermodynamic compressibilities, a genome-wide clustering of nucleosome sequences is obtained for C. albicans. This clustering shows that a chromosome beyond the classical dichotomic categories of hetero- and euchromatin is more feature-rich.

Prediction and comparative analysis of CTCF binding sites based on a first principle approach.

We calculated the patterns for the CCCTC transcription factor (CTCF) binding sites across many genomes on a first principle approach. The validation of the first principle method was done on the human as well as on the mouse genome. The predicted human CTCF binding sites are consistent with the consensus sequence, ChIP-seq data for the K562 cell, nucleosome positions for IMR90 cell as well as the CTCF binding sites in the mouse HOXA gene. The analysis ofHomo sapiens,Mus musculus,Sus scrofa,Capra hircusandDrosophila melanogasterwhole genomes shows: binding sites are organized in cluster-like groups, where two consecutive sites obey a power-law with coefficient ranging from 0.3292 ± 0.0068 to 0.5409 ± 0.0064; the distance between these groups varies from 18.08 ± 0.52 kbp to 42.1 ± 2.0 kbp. The genome ofAedes aegyptidoes not show a power law, but 19.9% of binding sites are 144 ± 4 and 287 ± 5 bp distant of each other. We run negative tests, confirming the under-representation of CTCF binding sites inCaenorhabditis elegans, Plasmodium falciparum andArabidopsis thalianacomplete genomes.

A chromosomal loop anchor mediates bacterial genome organization.

Nucleoprotein complexes play an integral role in genome organization of both eukaryotes and prokaryotes. Apart from their role in locally structuring and compacting DNA, several complexes are known to influence global organization by mediating long-range anchored chromosomal loop formation leading to spatial segregation of large sections of DNA. Such megabase-range interactions are ubiquitous in eukaryotes, but have not been demonstrated in prokaryotes. Here, using a genome-wide sedimentation-based approach, we found that a transcription factor, Rok, forms large nucleoprotein complexes in the bacterium Bacillus subtilis. Using chromosome conformation capture and live-imaging of DNA loci, we show that these complexes robustly interact with each other over large distances. Importantly, these Rok-dependent long-range interactions lead to anchored chromosomal loop formation, thereby spatially isolating large sections of DNA, as previously observed for insulator proteins in eukaryotes.

Telomerase subunit Est2 marks internal sites that are prone to accumulate DNA damage.

BACKGROUND: The main function of telomerase is at the telomeres but under adverse conditions telomerase can bind to internal regions causing deleterious effects as observed in cancer cells. RESULTS: By mapping the global occupancy of the catalytic subunit of telomerase (Est2) in the budding yeast Saccharomyces cerevisiae, we reveal that it binds to multiple guanine-rich genomic loci, which we termed "non-telomeric binding sites" (NTBS). We characterize Est2 binding to NTBS. Contrary to telomeres, Est2 binds to NTBS in G1 and G2 phase independently of Est1 and Est3. The absence of Est1 and Est3 renders telomerase inactive at NTBS. However, upon global DNA damage, Est1 and Est3 join Est2 at NTBS and telomere addition can be observed indicating that Est2 occupancy marks NTBS regions as particular risks for genome stability. CONCLUSIONS: Our results provide a novel model of telomerase regulation in the cell cycle using internal regions as "parking spots" of Est2 but marking them as hotspots for telomere addition.

Topological Analysis of γH2AX and MRE11 Clusters Detected by Localization Microscopy during X-ray-Induced DNA Double-Strand Break Repair.

DNA double-strand breaks (DSBs), known as the most severe damage in chromatin, were induced in breast cancer cells and normal skin fibroblasts by 2 Gy ionizing photon radiation. In response to DSB induction, phosphorylation of the histone variant H2AX to γH2AX was observed in the form of foci visualized by specific antibodies. By means of super-resolution single-molecule localization microscopy (SMLM), it has been recently shown in a first article about these data that these foci can be separated into clusters of about the same size (diameter ~400 nm). The number of clusters increased with the dose applied and decreased with the repair time. It has also been shown that during the repair period, antibody-labeled MRE11 clusters of about half of the γH2AX cluster diameter were formed inside several γH2AX clusters. MRE11 is part of the MRE11-RAD50-NBS1 (MRN) complex, which is known as a DNA strand resection and broken-end bridging component in homologous recombination repair (HRR) and alternative non-homologous end joining (a-NHEJ). This article is a follow-up of the former ones applying novel procedures of mathematics (topology) and similarity measurements on the data set: to obtain a measure for cluster shape and shape similarities, topological quantifications employing persistent homology were calculated and compared. In addition, based on our findings that γH2AX clusters associated with heterochromatin show a high degree of similarity independently of dose and repair time, these earlier published topological analyses and similarity calculations comparing repair foci within individual cells were extended by topological data averaging (2nd-generation heatmaps) over all cells analyzed at a given repair time point; thereby, the two dimensions (0 and 1) expressed by components and holes were studied separately. Finally, these mean value heatmaps were averaged, in addition. For γH2AX clusters, in both normal fibroblast and MCF-7 cancer cell lines, an increased similarity was found at early time points (up to 60 min) after irradiation for both components and holes of clusters. In contrast, for MRE11, the peak in similarity was found at later time points (2 h up to 48 h) after irradiation. In general, the normal fibroblasts showed quicker phosphorylation of H2AX and recruitment of MRE11 to γH2AX clusters compared to breast cancer cells and a shorter time interval of increased similarity for γH2AX clusters. γH2AX foci and randomly distributed MRE11 molecules naturally occurring in non-irradiated control cells did not show any significant topological similarity.

Elucidation of the Clustered Nano-Architecture of Radiation-Induced DNA Damage Sites and Surrounding Chromatin in Cancer Cells: A Single Molecule Localization Microscopy Approach.

In cancer therapy, the application of (fractionated) harsh radiation treatment is state of the art for many types of tumors. However, ionizing radiation is a "double-edged sword"-it can kill the tumor but can also promote the selection of radioresistant tumor cell clones or even initiate carcinogenesis in the normal irradiated tissue. Individualized radiotherapy would reduce these risks and boost the treatment, but its development requires a deep understanding of DNA damage and repair processes and the corresponding control mechanisms. DNA double strand breaks (DSBs) and their repair play a critical role in the cellular response to radiation. In previous years, it has become apparent that, beyond genetic and epigenetic determinants, the structural aspects of damaged chromatin (i.e., not only of DSBs themselves but also of the whole damage-surrounding chromatin domains) form another layer of complex DSB regulation. In the present article, we summarize the application of super-resolution single molecule localization microscopy (SMLM) for investigations of these structural aspects with emphasis on the relationship between the nano-architecture of radiation-induced repair foci (IRIFs), represented here by γH2AX foci, and their chromatin environment. Using irradiated HeLa cell cultures as an example, we show repair-dependent rearrangements of damaged chromatin and analyze the architecture of γH2AX repair clusters according to topological similarities. Although HeLa cells are known to have highly aberrant genomes, the topological similarity of γH2AX was high, indicating a functional, presumptively genome type-independent relevance of structural aspects in DSB repair. Remarkably, nano-scaled chromatin rearrangements during repair depended both on the chromatin domain type and the treatment. Based on these results, we demonstrate how the nano-architecture and topology of IRIFs and chromatin can be determined, point to the methodological relevance of SMLM, and discuss the consequences of the observed phenomena for the DSB repair network regulation or, for instance, radiation treatment outcomes.

Self-organised segregation of bacterial chromosomal origins.

The chromosomal replication origin region (ori) of characterised bacteria is dynamically positioned throughout the cell cycle. In slowly growing Escherichia coli, ori is maintained at mid-cell from birth until its replication, after which newly replicated sister oris move to opposite quarter positions. Here, we provide an explanation for ori positioning based on the self-organisation of the Structural Maintenance of Chromosomes complex, MukBEF, which forms dynamically positioned clusters on the chromosome. We propose that a non-trivial feedback between the self-organising gradient of MukBEF complexes and the oris leads to accurate ori positioning. We find excellent agreement with quantitative experimental measurements and confirm key predictions. Specifically, we show that oris exhibit biased motion towards MukBEF clusters, rather than mid-cell. Our findings suggest that MukBEF and oris act together as a self-organising system in chromosome organisation-segregation and introduces protein self-organisation as an important consideration for future studies of chromosome dynamics.

Processing and Analysis of Hi-C Data on Bacteria.

The study of three-dimensional genome organization has recently gained much attention in the context of novel techniques for detecting genome-wide contacts using next-generation sequencing. These genome-wide chromosome conformation capture-based methods, such as Hi-C, give a deep topological insight into the architecture of the genome inside the cell. This chapter reviews the steps to process next-generation Hi-C sequencing data to generate a final contact probability map. We describe these steps using publicly available Hi-C datasets of different bacteria. We also present strategies to assess the quality of Hi-C datasets.

Deciphering 3D Organization of Chromosomes Using Hi-C Data.

In order to interpret data from Hi-C studies genome-wide contact probability maps need to be translated into models of functional 3D genome organization. Here, we first present an overview of computational methods to analyze contact probability maps in terms of features such as the level and shape of compartmentalization. Next, we describe approaches to modeling 3D genome organization based on Hi-C data.

Using Persistent Homology as a New Approach for Super-Resolution Localization Microscopy Data Analysis and Classification of γH2AX Foci/Clusters.

DNA double strand breaks (DSB) are the most severe damages in chromatin induced by ionizing radiation. In response to such environmentally determined stress situations, cells have developed repair mechanisms. Although many investigations have contributed to a detailed understanding of repair processes, e.g., homologous recombination repair or non-homologous end-joining, the question is not sufficiently answered, how a cell decides to apply a certain repair process at a certain damage site, since all different repair pathways could simultaneously occur in the same cell nucleus. One of the first processes after DSB induction is phosphorylation of the histone variant H2AX to γH2AX in the given surroundings of the damaged locus. Since the spatial organization of chromatin is not random, it may be conclusive that the spatial organization of γH2AX foci is also not random, and rather, contributes to accessibility of special repair proteins to the damaged site, and thus, to the following repair pathway at this given site. The aim of this article is to demonstrate a new approach to analyze repair foci by their topology in order to obtain a cell independent method of categorization. During the last decade, novel super-resolution fluorescence light microscopic techniques have enabled new insights into genome structure and spatial organization on the nano-scale in the order of 10 nm. One of these techniques is single molecule localization microscopy (SMLM) with which the spatial coordinates of single fluorescence molecules can precisely be determined and density and distance distributions can be calculated. This method is an appropriate tool to quantify complex changes of chromatin and to describe repair foci on the single molecule level. Based on the pointillist information obtained by SMLM from specifically labeled heterochromatin and γH2AX foci reflecting the chromatin morphology and repair foci topology, we have developed a new analytical methodology of foci or foci cluster characterization, respectively, by means of persistence homology. This method allows, for the first time, a cell independent comparison of two point distributions (here the point distributions of two γH2AX clusters) with each other of a selected ensample and to give a mathematical measure of their similarity. In order to demonstrate the feasibility of this approach, cells were irradiated by low LET (linear energy transfer) radiation with different doses and the heterochromatin and γH2AX foci were fluorescently labeled by antibodies for SMLM. By means of our new analysis method, we were able to show that the topology of clusters of γH2AX foci can be categorized depending on the distance to heterochromatin. This method opens up new possibilities to categorize spatial organization of point patterns by parameterization of topological similarity.

Robust detection and segmentation of cell nuclei in biomedical images based on a computational topology framework.

The segmentation of cell nuclei is an important step towards the automated analysis of histological images. The presence of a large number of nuclei in whole-slide images necessitates methods that are computationally tractable in addition to being effective. In this work, a method is developed for the robust segmentation of cell nuclei in histological images based on the principles of persistent homology. More specifically, an abstract simplicial homology approach for image segmentation is established. Essentially, the approach deals with the persistence of disconnected sets in the image, thus identifying salient regions that express patterns of persistence. By introducing an image representation based on topological features, the task of segmentation is less dependent on variations of color or texture. This results in a novel approach that generalizes well and provides stable performance. The method conceptualizes regions of interest (cell nuclei) pertinent to their topological features in a successful manner. The time cost of the proposed approach is lower-bounded by an almost linear behavior and upper-bounded by O(n2) in a worst-case scenario. Time complexity matches a quasilinear behavior which is O(n1+ɛ) for ε < 1. Images acquired from histological sections of liver tissue are used as a case study to demonstrate the effectiveness of the approach. The histological landscape consists of hepatocytes and non-parenchymal cells. The accuracy of the proposed methodology is verified against an automated workflow created by the output of a conventional filter bank (validated by experts) and the supervised training of a random forest classifier. The results are obtained on a per-object basis. The proposed workflow successfully detected both hepatocyte and non-parenchymal cell nuclei with an accuracy of 84.6%, and hepatocyte cell nuclei only with an accuracy of 86.2%. A public histological dataset with supplied ground-truth data is also used for evaluating the performance of the proposed approach (accuracy: 94.5%). Further validations are carried out with a publicly available dataset and ground-truth data from the Gland Segmentation in Colon Histology Images Challenge (GlaS) contest. The proposed method is useful for obtaining unsupervised robust initial segmentations that can be further integrated in image/data processing and management pipelines. The development of a fully automated system supporting a human expert provides tangible benefits in the context of clinical decision-making.

A Three-Pronged Attack To Investigate the Electronic Structure of a Family of Ferromagnetic Fe4Ln2 Cyclic Coordination Clusters: A Combined Magnetic Susceptibility, High-Field/High-Frequency Electron Paramagnetic Resonance, and 57Fe Mössbauer Study.

We present the synthesis, structure, magnetic properties, as well as the Mössbauer and electron paramagnetic resonance studies of a ring-shaped [FeIII4LnIII2(Htea)4(µ-N3)4(N3)3(piv)3] (Ln = Y 1, Gd 2, Tb 3, Dy 4, Ho 5, Er, 6) coordination cluster. The Dy, Tb, and Ho analogues show blocking of the magnetization at low temperatures without applied fields. The anisotropy of the 3d ion and the exchange interaction between 3d and 4f ions in Fe4Ln2 complexes are unambiguously determined by high-field/high-frequency electron paramagnetic resonance measurements at low temperature. Ferromagnetic exchange interaction JFe-Ln is found which decreases upon variation of the Ln ions to larger atomic numbers. This dependence is similar to the behavior shown in the effective barrier values of complexes 3-5. Further information about the anisotropy of the Ln3+ ions was gathered with 57Fe Mössbauer spectroscopy, and the combination of these methods provides detailed information regarding the electronic structure of these complexes.

Viral channel forming proteins--How to assemble and depolarize lipid membranes in silico.

Viral channel forming proteins (VCPs) have been discovered in the late 70s and are found in many viruses to date. Usually they are small and have to assemble to form channels which depolarize the lipid membrane of the host cells. Structural information is just about to emerge for just some of them. Thus, computational methods play a pivotal role in generating plausible structures which can be used in the drug development process. In this review the accumulation of structural data is introduced from a historical perspective. Computational performances and their predictive power are reported guided by biological questions such as the assembly, mechanism of function and drug-protein interaction of VCPs. An outlook of how coarse grained simulations can contribute to yet unexplored issues of these proteins is given. This article is part of a Special Issue entitled: Membrane Proteins edited by J.C. Gumbart and Sergei Noskov.

Radiation induced chromatin conformation changes analysed by fluorescent localization microscopy, statistical physics, and graph theory.

It has been well established that the architecture of chromatin in cell nuclei is not random but functionally correlated. Chromatin damage caused by ionizing radiation raises complex repair machineries. This is accompanied by local chromatin rearrangements and structural changes which may for instance improve the accessibility of damaged sites for repair protein complexes. Using stably transfected HeLa cells expressing either green fluorescent protein (GFP) labelled histone H2B or yellow fluorescent protein (YFP) labelled histone H2A, we investigated the positioning of individual histone proteins in cell nuclei by means of high resolution localization microscopy (Spectral Position Determination Microscopy = SPDM). The cells were exposed to ionizing radiation of different doses and aliquots were fixed after different repair times for SPDM imaging. In addition to the repair dependent histone protein pattern, the positioning of antibodies specific for heterochromatin and euchromatin was separately recorded by SPDM. The present paper aims to provide a quantitative description of structural changes of chromatin after irradiation and during repair. It introduces a novel approach to analyse SPDM images by means of statistical physics and graph theory. The method is based on the calculation of the radial distribution functions as well as edge length distributions for graphs defined by a triangulation of the marker positions. The obtained results show that through the cell nucleus the different chromatin re-arrangements as detected by the fluorescent nucleosomal pattern average themselves. In contrast heterochromatic regions alone indicate a relaxation after radiation exposure and re-condensation during repair whereas euchromatin seemed to be unaffected or behave contrarily. SPDM in combination with the analysis techniques applied allows the systematic elucidation of chromatin re-arrangements after irradiation and during repair, if selected sub-regions of nuclei are investigated.

A multiscale approach to simulating the conformational properties of unbound multi-C2H2 zinc finger proteins.

The conformational properties of unbound multi-Cys2 His2 (mC2H2) zinc finger proteins, in which zinc finger domains are connected by flexible linkers, are studied by a multiscale approach. Three methods on different length scales are utilized. First, atomic detail molecular dynamics simulations of one zinc finger and its adjacent flexible linker confirmed that the zinc finger is more rigid than the flexible linker. Second, the end-to-end distance distributions of mC2H2 zinc finger proteins are computed using an efficient atomistic pivoting algorithm, which only takes excluded volume interactions into consideration. The end-to-end distance distribution gradually changes its profile, from left-tailed to right-tailed, as the number of zinc fingers increases. This is explained by using a worm-like chain model. For proteins of a few zinc fingers, an effective bending constraint favors an extended conformation. Only for proteins containing more than nine zinc fingers, is a somewhat compacted conformation preferred. Third, a mesoscale model is modified to study both the local and the global conformational properties of multi-C2H2 zinc finger proteins. Simulations of the CCCTC-binding factor (CTCF), an important mC2H2 zinc finger protein for genome spatial organization, are presented.

The role of loops on the order of eukaryotes and prokaryotes.

The study of the three-dimensional organization of chromatin has recently gained much focus in the context of novel techniques for detecting genome-wide contacts using next-generation sequencing. These chromosome conformation capture-based methods give a deep topological insight into the architecture of the genome inside the nucleus. Several recent studies observe a compartmentalization of chromatin interactions into spatially confined domains. This structural feature of interphase chromosomes is not only supported by conventional studies assessing the interaction data of millions of cells, but also by analysis on the level of a single cell. We first present and examine the different models that have been proposed to elucidate these topological domains in eukaryotes. Then we show that a model which relies on the dynamic formation of loops within domains can account for the experimentally observed contact maps. Interestingly, the topological domain structure is not only found in mammalian genomes, but also in bacterial chromosomes.