Specific Sequence Preferences Research Articles

Modelling the evolution of transcription factor binding preferences in complex eukaryotes

Transcription factors (TFs) exert their regulatory action by binding to DNA with specific sequence preferences. However, different TFs can partially share their binding sequences due to their common evolutionary origin. This “redundancy” of binding defines a way of organizing TFs in “motif families” by grouping TFs with similar binding preferences. Since these ultimately define the TF target genes, the motif family organization entails information about the structure of transcriptional regulation as it has been shaped by evolution. Focusing on the human TF repertoire, we show that a one-parameter evolutionary model of the Birth-Death-Innovation type can explain the TF empirical repartition in motif families, and allows to highlight the relevant evolutionary forces at the origin of this organization. Moreover, the model allows to pinpoint few deviations from the neutral scenario it assumes: three over-expanded families (including HOX and FOX genes), a set of “singleton” TFs for which duplication seems to be selected against, and a higher-than-average rate of diversification of the binding preferences of TFs with a Zinc Finger DNA binding domain. Finally, a comparison of the TF motif family organization in different eukaryotic species suggests an increase of redundancy of binding with organism complexity.

High-throughput functional testing of ENCODE segmentation predictions.

The histone modification state of genomic regions is hypothesized to reflect the regulatory activity of the underlying genomic DNA. Based on this hypothesis, the ENCODE Project Consortium measured the status of multiple histone modifications across the genome in several cell types and used these data to segment the genome into regions with different predicted regulatory activities. We measured the cis-regulatory activity of more than 2000 of these predictions in the K562 leukemia cell line. We tested genomic segments predicted to be Enhancers, Weak Enhancers, or Repressed elements in K562 cells, along with other sequences predicted to be Enhancers specific to the H1 human embryonic stem cell line (H1-hESC). Both Enhancer and Weak Enhancer sequences in K562 cells were more active than negative controls, although surprisingly, Weak Enhancer segmentations drove expression higher than did Enhancer segmentations. Lower levels of the covalent histone modifications H3K36me3 and H3K27ac, thought to mark active enhancers and transcribed gene bodies, associate with higher expression and partly explain the higher activity of Weak Enhancers over Enhancer predictions. While DNase I hypersensitivity (HS) is a good predictor of active sequences in our assay, transcription factor (TF) binding models need to be included in order to accurately identify highly expressed sequences. Overall, our results show that a significant fraction (∼26%) of the ENCODE enhancer predictions have regulatory activity, suggesting that histone modification states can reflect the cis-regulatory activity of sequences in the genome, but that specific sequence preferences, such as TF-binding sites, are the causal determinants of cis-regulatory activity.

Understanding the Sequence Preference of Recurrent RNA Building Blocks using Quantum Chemistry: The Intrastrand RNA Dinucleotide Platform.

Folded RNA molecules are shaped by an astonishing variety of highly conserved noncanonical molecular interactions and backbone topologies. The dinucleotide platform is a widespread recurrent RNA modular building submotif formed by the side-by-side pairing of bases from two consecutive nucleotides within a single strand, with highly specific sequence preferences. This unique arrangement of bases is cemented by an intricate network of noncanonical hydrogen bonds and facilitated by a distinctive backbone topology. The present study investigates the gas-phase intrinsic stabilities of the three most common RNA dinucleotide platforms - 5'-GpU-3', ApA, and UpC - via state-of-the-art quantum-chemical (QM) techniques. The mean stability of base-base interactions decreases with sequence in the order GpU > ApA > UpC. Bader's atoms-in-molecules analysis reveals that the N2(G)…O4(U) hydrogen bond of the GpU platform is stronger than the corresponding hydrogen bonds in the other two platforms. The mixed-pucker sugar-phosphate backbone conformation found in most GpU platforms, in which the 5'-ribose sugar (G) is in the C2'-endo form and the 3'-sugar (U) in the C3'-endo form, is intrinsically more stable than the standard A-RNA backbone arrangement, partially as a result of a favorable O2'…O2P intra-platform interaction. Our results thus validate the hypothesis of Lu et al. (Lu Xiang-Jun, et al. Nucleic Acids Res. 2010, 38, 4868-4876), that the superior stability of GpU platforms is partially mediated by the strong O2'…O2P hydrogen bond. In contrast, ApA and especially UpC platform-compatible backbone conformations are rather diverse and do not display any characteristic structural features. The average stabilities of ApA and UpC derived backbone conformers are also lower than those of GpU platforms. Thus, the observed structural and evolutionary patterns of the dinucleotide platforms can be accounted for, to a large extent, by their intrinsic properties as described by modern QM calculations. In contrast, we show that the dinucleotide platform is not properly described in the course of atomistic explicit-solvent simulations. Our work also gives methodological insights into QM calculations of experimental RNA backbone geometries. Such calculations are inherently complicated by rather large data and refinement uncertainties in the available RNA experimental structures, which often preclude reliable energy computations.

Distinct p53 genomic binding patterns in normal and cancer-derived human cells

We report here genome-wide analysis of the tumor suppressor p53 binding sites in normal human cells. 743 high-confidence ChIP-seq peaks representing putative genomic binding sites were identified in normal IMR90 fibroblasts using a reference chromatin sample. More than 40% were located within 2 kb of a transcription start site (TSS), a distribution similar to that documented for individually studied, functional p53 binding sites and, to date, not observed by previous p53 genome-wide studies. Nearly half of the high-confidence binding sites in the IMR90 cells reside in CpG islands, in marked contrast to sites reported in cancer-derived cells. The distinct genomic features of the IMR90 binding sites do not reflect a distinct preference for specific sequences, since the de novo developed p53 motif based on our study is similar to those reported by genome-wide studies of cancer cells. More likely, the different chromatin landscape in normal, compared with cancer-derived cells, influences p53 binding via modulating availability of the sites. We compared the IMR90 ChIPseq peaks to the recently published IMR90 methylome1 and demonstrated that they are enriched at hypomethylated DNA. Our study represents the first genome-wide, de novo mapping of p53 binding sites in normal human cells and reveals that p53 binding sites reside in distinct genomic landscapes in normal and cancer-derived human cells.

Comparative analysis of H2A.Z nucleosome organization in the human and yeast genomes

Eukaryotic DNA is wrapped around a histone protein core to constitute the fundamental repeating units of chromatin, the nucleosomes. The affinity of the histone core for DNA depends on the nucleotide sequence; however, it is unclear to what extent DNA sequence determines nucleosome positioning in vivo, and if the same rules of sequence-directed positioning apply to genomes of varying complexity. Using the data generated by high-throughput DNA sequencing combined with chromatin immunoprecipitation, we have identified positions of nucleosomes containing the H2A.Z histone variant and histone H3 trimethylated at lysine 4 in human CD4(+) T-cells. We find that the 10-bp periodicity observed in nucleosomal sequences in yeast and other organisms is not pronounced in human nucleosomal sequences. This result was confirmed for a broader set of mononucleosomal fragments that were not selected for any specific histone variant or modification. We also find that human H2A.Z nucleosomes protect only approximately 120 bp of DNA from MNase digestion and exhibit specific sequence preferences, suggesting a novel mechanism of nucleosome organization for the H2A.Z variant.

Temperature- and H-NS-dependent regulation of a plasmid-encoded virulence operon expressing Escherichia coli hemolysin.

Proteins H-NS and Hha form a nucleoprotein complex that modulates expression of the thermoregulated hly operon of Escherichia coli. We have been able to identify two H-NS binding sites in the hly regulatory region. One of them partially overlaps the promoter region (site II), and the other is located about 2 kbp upstream (site I). In contrast, Hha protein did not show any preference for specific sequences. In vitro, temperature influences the affinity of H-NS for a DNA fragment containing both binding sites and H-NS-mediated repression of hly operon transcription. Deletion analysis of the hly regulatory region confirms the relevance of site I for thermoregulation of this operon. We present a model to explain the temperature-modulated repression of the hly operon, based on the experiments reported here and other, preexisting data.

The Sso7d DNA-binding protein from Sulfolobus solfataricus has ribonuclease activity

Sso7d is a small, basic, abundant protein from the thermoacidophilic archaeon Sulfolobus solfataricus. Previous research has shown that Sso7d can bind double-stranded DNA without sequence specificity by placing its triple-stranded β-sheet across the minor groove. We previously found RNase activity both in preparations of Sso7d purified from its natural source and in recombinant, purified protein expressed in Escherichia coli. This paper provides conclusive evidence that supports the assignment of RNase activity to Sso7d, shown by the total absence of activity in the single-point mutants E35L and K12L, despite the preservation of their overall structure under the assay conditions. In keeping with our observation that the residues putatively involved in RNase activity and those playing a role in DNA binding are located on different surfaces of the molecule, the activity was not impaired in the presence of DNA. If a small synthetic RNA was used as a substrate, Sso7d attacked both predicted double- and single-stranded RNA stretches, with no evident preference for specific sequences or individual bases. Apparently, the more readily attacked bonds were those intrinsically more unstable.

Localized chemical reactivity in double-stranded DNA associated with the intercalative binding of benzo[e]pyridoindole and benzo[g]pyridoindole triple-helix-stabilizing ligands.

Footprinting with methidiumpropyl-EDTA.FeII has been used to map the binding sites on duplex DNA of two closely related benzopyridoindole derivatives which selectively stabilize triple-helical DNA-oligonucleotide complexes. Both ligands bind to many sites, including certain oligopurine.oligopyrimidine tracts, with a weak preference for some (but not all) sequences rich in A.T base pairs. This indifference to primary sequence, with evidence of binding to the commonly disfavoured (A)n.(T)ntracts, may at least partially explain why the ligands stabilize triplex structures composed of T.A.T pairings. Neither 3-methoxy-7H-8-methyl-11- [(3'amino)propylamino]benzo[e]pyrido[4, 3-b]indole (BePI) nor 3- methoxy-7-[3'-diethylamino)propylamino]-10-methyl-11H- benzo[g]pyrido[4,3-b]indole (BgPI) affect the reaction of dimethyl sulphate or potassium tetrachloropalladinate with the N7 of purines in the major groove, but both enhance the reactivity of purines (mostly adenine residues) towards diethylpyrocarbonate, both proximal and distal to their identified binding sites. With potassium permanganate and osmium tetroxide/pyridine, probes for the accessibility of the 5,6 double bond of pyrimidine residues, BgPI has a more potent effect than BePI and, generally, the reaction with KMnO4 is more pronounced than that with OsO4. BgPI conspicuously potentiates the oxidation of pyrimidines in the triplet sequences 3'-ATA, 3'-GTA and 3'-GCA, whereas BePI enhances the reactivity of OsO4 towards thymine in sequences 3'-ATYR, with no effect on cytosine residues. Thus, despite their structural homology and common lack of specific sequence preferences, the two benzopyridoindole derivatives induce distinct conformational changes in duplex DNA, not just within the sites where footprints can be detected.

Mechanisms of DNA sequence selective alkylation of guanine-N7 positions by nitrogen mustards

Abstract Quantitative determinations were carried out of the relative reaction rates of several nitrogen mustards at various guanine-N7 positions in DNA fragments of known sequence. The findings suggest structural hypotheses of the origins of the reaction selectivities. End-labeled DNA fragments were reacted with nitrogen mustards, and the guanine-N7 alkylation sites were analyzed by gel electrophoresis. Relative reaction intensities were determined by computer analysis of digitized densitometer scans. The differences in reaction intensities at different G's were in part attributable to the effects of nearest neighbor base pairs on the molecular electrostatic potential near the reaction site. Uracil and quinacrine mustards have specific sequence preferences for reaction that differ from other mustards. The nature of the specific sequence preferences were determined and hypotheses are proposed to explain their origin.

Mechanisms of DNA sequence selective alkylation of guanine-N7 positions by nitrogen mustards.

Quantitative determinations were carried out of the relative reaction rates of several nitrogen mustards at various guanine-N7 positions in DNA fragments of known sequence. The findings suggest structural hypotheses of the origins of the reaction selectivities. End-labeled DNA fragments were reacted with nitrogen mustards, and the guanine-N7 alkylation sites were analyzed by gel electrophoresis. Relative reaction intensities were determined by computer analysis of digitized densitometer scans. The differences in reaction intensities at different G's were in part attributable to the effects of nearest neighbor base pairs on the molecular electrostatic potential near the reaction site. Uracil and quinacrine mustards have specific sequence preferences for reaction that differ from other mustards. The nature of the specific sequence preferences were determined and hypotheses are proposed to explain their origin.