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.

Reconstruction of the free energy in the metastable region using the path ensemble.

By quenching into the metastable region of the three-dimensional Ising model, we investigate the paths that the magnetization (energy) takes as a function of time. We accumulate the magnetization (energy) paths into time-dependent distributions from which we reconstruct the free energy as a function of the magnetic field, temperature, and system size. From the reconstructed free energy, we obtain the free-energy barrier that is associated with the transition from a metastable state to the stable equilibrium state. Although mean-field theory predicts a sharp transition between the metastable and the unstable region where the free-energy barrier is zero, the results for the nearest-neighbor Ising model show that the free-energy barrier does not go to zero.