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.

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.

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.

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.

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.

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.

A topological similarity measure for proteins.

We introduce a new measure for assessing similarity among chemical structures, based on well-established computational-topology algorithms. We argue that although the method considers geometry, it is more than a mere geometric similarity measure, as it takes into account, on different geometric scales, the important topological features of the compared structures. We prove that our measure is rigorous and complies with the proper mathematical requirements. We validate the method through comparing different configurations of simple zinc finger proteins and present an application on ligands binding to membrane-proteINS extracted from the Directory of Useful Decoys: Enhanced database and corresponding decoys. This article is part of a Special Issue entitled: Viral membrane proteins - Channels for cellular networking.

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.

DNA double-strand breaks: linking gene expression to chromosome morphology and mobility.

Ionizing radiation can lead to DNA double-strand breaks (DSBs) which belong to the most dangerous forms of damage to the DNA. Cells possess elaborate repair mechanisms and react in a complex manner to the emergence of DSBs. Experiments have shown that gene expression levels in irradiated cells are changed, and thousands of radiation-responsive genes have been identified. On the other hand, recent studies have shown that gene expression is tightly connected to the three-dimensional organization of the genome. In this work, we analyzed the chromatin organization in the cell nuclei before and after exposure to ionizing radiation with an expression-dependent folding model. Our results indicate that the alteration of the chromosome organization on the scale of a complete chromosome is rather limited despite the expression level change of a large number of genes. We further modelled breaks within sub-compartments of the model chromosomes and showed that entropic changes caused by a break lead to increased mobility of the break sites and help to locate break ends further to the periphery of the sub-compartments. We conclude that the changes in the chromatin structure after irradiation are limited to local scales and demonstrate the importance of entropy for the behaviour of break ends.

Depletion of the chromatin looping proteins CTCF and cohesin causes chromatin compaction: insight into chromatin folding by polymer modelling.

Folding of the chromosomal fibre in interphase nuclei is an important element in the regulation of gene expression. For instance, physical contacts between promoters and enhancers are a key element in cell-type-specific transcription. We know remarkably little about the principles that control chromosome folding. Here we explore the view that intrachromosomal interactions, forming a complex pattern of loops, are a key element in chromosome folding. CTCF and cohesin are two abundant looping proteins of interphase chromosomes of higher eukaryotes. To investigate the role of looping in large-scale (supra Mb) folding of human chromosomes, we knocked down the gene that codes for CTCF and the one coding for Rad21, an essential subunit of cohesin. We measured the effect on chromosome folding using systematic 3D fluorescent in situ hybridization (FISH). Results show that chromatin becomes more compact after reducing the concentration of these two looping proteins. The molecular basis for this counter-intuitive behaviour is explored by polymer modelling usingy the Dynamic Loop model (Bohn M, Heermann DW (2010) Diffusion-driven looping provides a consistent framework for chromatin organization. PLoS ONE 5: e12218.). We show that compaction can be explained by selectively decreasing the number of short-range loops, leaving long-range looping unchanged. In support of this model prediction it has recently been shown by others that CTCF and cohesin indeed are responsible primarily for short-range looping. Our results suggest that the local and the overall changes in of chromosome structure are controlled by a delicate balance between short-range and long-range loops, allowing easy switching between, for instance, open and more compact chromatin states.

Chromosome segregation by the Escherichia coli Min system.

The mechanisms underlying chromosome segregation in prokaryotes remain a subject of debate and no unifying view has yet emerged. Given that the initial disentanglement of duplicated chromosomes could be achieved by purely entropic forces, even the requirement of an active prokaryotic segregation machinery has been questioned. Using computer simulations, we show that entropic forces alone are not sufficient to achieve and maintain full separation of chromosomes. This is, however, possible by assuming repeated binding of chromosomes along a gradient of membrane-associated tethering sites toward the poles. We propose that, in Escherichia coli, such a gradient of membrane tethering sites may be provided by the oscillatory Min system, otherwise known for its role in selecting the cell division site. Consistent with this hypothesis, we demonstrate that MinD binds to DNA and tethers it to the membrane in an ATP-dependent manner. Taken together, our combined theoretical and experimental results suggest the existence of a novel mechanism of chromosome segregation based on the Min system, further highlighting the importance of active segregation of chromosomes in prokaryotic cell biology.

A model for Escherichia coli chromosome packaging supports transcription factor-induced DNA domain formation.

What physical mechanism leads to organization of a highly condensed and confined circular chromosome? Computational modeling shows that confinement-induced organization is able to overcome the chromosome's propensity to mix by the formation of topological domains. The experimentally observed high precision of separate subcellular positioning of loci (located on different chromosomal domains) in Escherichia coli naturally emerges as a result of entropic demixing of such chromosomal loops. We propose one possible mechanism for organizing these domains: regulatory control defined by the underlying E. coli gene regulatory network requires the colocalization of transcription factor genes and target genes. Investigating this assumption, we find the DNA chain to self-organize into several topologically distinguishable domains where the interplay between the entropic repulsion of chromosomal loops and their compression due to the confining geometry induces an effective nucleoid filament-type of structure. Thus, we propose that the physical structure of the chromosome is a direct result of regulatory interactions. To reproduce the observed precise ordering of the chromosome, we estimate that the domain sizes are distributed between 10 and 700 kb, in agreement with the size of topological domains identified in the context of DNA supercoiling.

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.

Chromatin folding--from biology to polymer models and back.

There is rapidly growing evidence that folding of the chromatin fibre inside the interphase nucleus has an important role in the regulation of gene expression. In particular, the formation of loops mediated by the interaction between specific regulatory elements, for instance enhancers and promoters, is crucial in gene control. Biochemical studies that were based on the chromosome conformation capture (3C) technology have confirmed that eukaryotic genomes are highly looped. Insight into the underlying principles comes from polymer models that explore the properties of the chromatin fibre inside the nucleus. Recent models indicate that chromatin looping can explain various properties of interphase chromatin, including chromatin compaction and compartmentalisation of chromosomes. Entropic effects have a key role in these models. In this Commentary, we give an overview of the recent conjunction of ideas regarding chromatin looping in the fields of biology and polymer physics. Starting from simple linear polymer models, we explain how specific folding properties emerge upon introducing loops and how this explains a variety of experimental observations. We also discuss different polymer models that describe chromatin folding and compare them to experimental data. Experimentally testing the predictions of such polymer models and their subsequent improvement on the basis of measurements provides a solid framework to begin to understand how our genome is folded and how folding relates to function.

A model for the 3D chromatin architecture of pro and eukaryotes.

How do chromosomes fold and how does this depend on the organism and type of cell? This question has been addressed by a model and a methodology that uses only very basic constituents to capture the relevant features of folding. Key is the dynamic formation of loops within the chromosome. With this and entropy we show that the model is capable to describe the folding of human chromosomes in inter- and metaphase as well as for the Escherichia coli circular chromosomes.

Mitotic chromosome structure.

Mounting evidence is compiling linking the physical organizational structure of chromosomes and the nuclear structure to biological function. At the base of the physical organizational structure of both is the concept of loop formation. This implies that physical proximity within chromosomes is provided for otherwise distal genomic regions and thus hierarchically organizing the chromosomes. Together with entropy many experimental observations can be explained with these two concepts. Among the observations that can be explained are the measured physical extent of the chromosomes, their shape, mechanical behavior, the segregation into territories (chromosomal and territories within chromosomes), the results from chromosome conformation capture experiments, as well as linking gene expression to structural organization.

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.

Spatially confined folding of chromatin in the interphase nucleus.

Genome function in higher eukaryotes involves major changes in the spatial organization of the chromatin fiber. Nevertheless, our understanding of chromatin folding is remarkably limited. Polymer models have been used to describe chromatin folding. However, none of the proposed models gives a satisfactory explanation of experimental data. In particularly, they ignore that each chromosome occupies a confined space, i.e., the chromosome territory. Here, we present a polymer model that is able to describe key properties of chromatin over length scales ranging from 0.5 to 75 Mb. This random loop (RL) model assumes a self-avoiding random walk folding of the polymer backbone and defines a probability P for 2 monomers to interact, creating loops of a broad size range. Model predictions are compared with systematic measurements of chromatin folding of the q-arms of chromosomes 1 and 11. The RL model can explain our observed data and suggests that on the tens-of-megabases length scale P is small, i.e., 10-30 loops per 100 Mb. This is sufficient to enforce folding inside the confined space of a chromosome territory. On the 0.5- to 3-Mb length scale chromatin compaction differs in different subchromosomal domains. This aspect of chromatin structure is incorporated in the RL model by introducing heterogeneity along the fiber contour length due to different local looping probabilities. The RL model creates a quantitative and predictive framework for the identification of nuclear components that are responsible for chromatin-chromatin interactions and determine the 3-dimensional organization of the chromatin fiber.

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.