Design and analysis of ChIP-Seq experiments for nuclear factor I DNA-binding proteins

SUMMARY : Eukaryotic DNA interacts with the nuclear proteins using non-covalent ionic interactions. Proteins can recognize specific nucleotide sequences based on the sterical interactions with the DNA and these specific protein-DNA interactions are the basis for many nuclear processes, e.g. gene transcription, chromosomal replication, and recombination. New technology termed ChIP-Seq has been recently developed for the analysis of protein-DNA interactions on a whole genome scale and it is based on immunoprecipitation of chromatin and high-throughput DNA sequencing procedure. ChIP-Seq is a novel technique with a great potential to replace older techniques for mapping of protein-DNA interactions. In this thesis, we bring some new insights into the ChIP-Seq data analysis. First, we point out to some common and so far unknown artifacts of the method. Sequence tag distribution in the genome does not follow uniform distribution and we have found extreme hot-spots of tag accumulation over specific loci in the human and mouse genomes. These artifactual sequence tags accumulations will create false peaks in every ChIP-Seq dataset and we propose different filtering methods to reduce the number of false positives. Next, we propose random sampling as a powerful analytical tool in the ChIP-Seq data analysis that could be used to infer biological knowledge from the massive ChIP-Seq datasets. We created unbiased random sampling algorithm and we used this methodology to reveal some of the important biological properties of Nuclear Factor I DNA binding proteins. Finally, by analyzing the ChIP-Seq data in detail, we revealed that Nuclear Factor I transcription factors mainly act as activators of transcription, and that they are associated with specific chromatin modifications that are markers of open chromatin. We speculate that NFI factors only interact with the DNA wrapped around the nucleosome. We also found multiple loci that indicate possible chromatin barrier activity of NFI proteins, which could suggest the use of NFI binding sequences as chromatin insulators in biotechnology applications. RESUME : L'ADN des eucaryotes interagit avec les protéines nucléaires par des interactions noncovalentes ioniques. Les protéines peuvent reconnaître les séquences nucléotidiques spécifiques basées sur l'interaction stérique avec l'ADN, et des interactions spécifiques contrôlent de nombreux processus nucléaire, p.ex. transcription du gène, la réplication chromosomique, et la recombinaison. Une nouvelle technologie appelée ChIP-Seq a été récemment développée pour l'analyse des interactions protéine-ADN à l'échelle du génome entier et cette approche est basée sur l'immuno-précipitation de la chromatine et sur la procédure de séquençage de l'ADN à haut débit. La nouvelle approche ChIP-Seq a donc un fort potentiel pour remplacer les anciennes techniques de cartographie des interactions protéine-ADN. Dans cette thèse, nous apportons de nouvelles perspectives dans l'analyse des données ChIP-Seq. Tout d'abord, nous avons identifié des artefacts très communs associés à cette méthode qui étaient jusqu'à présent insoupçonnés. La distribution des séquences dans le génome ne suit pas une distribution uniforme et nous avons constaté des positions extrêmes d'accumulation de séquence à des régions spécifiques, des génomes humains et de la souris. Ces accumulations des séquences artéfactuelles créera de faux pics dans toutes les données ChIP-Seq, et nous proposons différentes méthodes de filtrage pour réduire le nombre de faux positifs. Ensuite, nous proposons un nouvel échantillonnage aléatoire comme un outil puissant d'analyse des données ChIP-Seq, ce qui pourraient augmenter l'acquisition de connaissances biologiques à partir des données ChIP-Seq. Nous avons créé un algorithme d'échantillonnage aléatoire et nous avons utilisé cette méthode pour révéler certaines des propriétés biologiques importantes de protéines liant à l'ADN nommés Facteur Nucléaire I (NFI). Enfin, en analysant en détail les données de ChIP-Seq pour la famille de facteurs de transcription nommés Facteur Nucléaire I, nous avons révélé que ces protéines agissent principalement comme des activateurs de transcription, et qu'elles sont associées à des modifications de la chromatine spécifiques qui sont des marqueurs de la chromatine ouverte. Nous pensons que lés facteurs NFI interagir uniquement avec l'ADN enroulé autour du nucléosome. Nous avons également constaté plusieurs régions génomiques qui indiquent une éventuelle activité de barrière chromatinienne des protéines NFI, ce qui pourrait suggérer l'utilisation de séquences de liaison NFI comme séquences isolatrices dans des applications de la biotechnologie.

Proof of concept for microarray-based detection of DNA-binding oncogenes in cell extracts.

The function of DNA-binding proteins is controlled not just by their abundance, but mainly at the level of their activity in terms of their interactions with DNA and protein targets. Moreover, the affinity of such transcription factors to their target sequences is often controlled by co-factors and/or modifications that are not easily assessed from biological samples. Here, we describe a scalable method for monitoring protein-DNA interactions on a microarray surface. This approach was designed to determine the DNA-binding activity of proteins in crude cell extracts, complementing conventional expression profiling arrays. Enzymatic labeling of DNA enables direct normalization of the protein binding to the microarray, allowing the estimation of relative binding affinities. Using DNA sequences covering a range of affinities, we show that the new microarray-based method yields binding strength estimates similar to low-throughput gel mobility-shift assays. The microarray is also of high sensitivity, as it allows the detection of a rare DNA-binding protein from breast cancer cells, the human tumor suppressor AP-2. This approach thus mediates precise and robust assessment of the activity of DNA-binding proteins and takes present DNA-binding assays to a high throughput level.

DNA-lipid systems: A physical chemistry study

It is well known that the interaction of polyelectrolytes with oppositely charged surfactants leads to an associative phase separation; however, the phase behavior of DNA and oppositely charged surfactants is more strongly associative than observed in other systems. A precipitate is formed with very low amounts of surfactant and DNA. DNA compaction is a general phenomenon in the presence of multivalent ions and positively charged surfaces; because of the high charge density there are strong attractive ion correlation effects. Techniques like phase diagram determinations, fluorescence microscopy, and ellipsometry were used to study these systems. The interaction between DNA and catanionic mixtures (i.e., mixtures of cationic and anionic surfactants) was also investigated. We observed that DNA compacts and adsorbs onto the surface of positively charged vesicles, and that the addition of an anionic surfactant can release DNA back into solution from a compact globular complex between DNA and the cationic surfactant. Finally, DNA interactions with polycations, chitosans with different chain lengths, were studied by fluorescence microscopy, in vivo transfection assays and cryogenic transmission electron microscopy. The general conclusion is that a chitosan effective in promoting compaction is also efficient in transfection.

The Xanthomonas citri effector protein PthA interacts with citrus proteins involved in nuclear transport, protein folding and ubiquitination associated with DNA repair

P>Xanthomonas axonopodis pv. citri utilizes the type III effector protein PthA to modulate host transcription to promote citrus canker. PthA proteins belong to the AvrBs3/PthA family and carry a domain comprising tandem repeats of 34 amino acids that mediates protein-protein and protein-DNA interactions. We show here that variants of PthAs from a single bacterial strain localize to the nucleus of plant cells and form homo- and heterodimers through the association of their repeat regions. We hypothesize that the PthA variants might also interact with distinct host targets. Here, in addition to the interaction with alpha-importin, known to mediate the nuclear import of AvrBs3, we describe new interactions of PthAs with citrus proteins involved in protein folding and K63-linked ubiquitination. PthAs 2 and 3 preferentially interact with a citrus cyclophilin (Cyp) and with TDX, a tetratricopeptide domain-containing thioredoxin. In addition, PthAs 2 and 3, but not 1 and 4, interact with the ubiquitin-conjugating enzyme complex formed by Ubc13 and ubiquitin-conjugating enzyme variant (Uev), required for K63-linked ubiquitination and DNA repair. We show that Cyp, TDX and Uev interact with each other, and that Cyp and Uev localize to the nucleus of plant cells. Furthermore, the citrus Ubc13 and Uev proteins complement the DNA repair phenotype of the yeast Delta ubc13 and Delta mms2/uev1a mutants, strongly indicating that they are also involved in K63-linked ubiquitination and DNA repair. Notably, PthA 2 affects the growth of yeast cells in the presence of a DNA damage agent, suggesting that it inhibits K63-linked ubiquitination required for DNA repair.

Transfection of drug-specific T-cell receptors into hybridoma cells: tools to monitor drug interaction with T-cell receptors and evaluate cross-reactivity to related compounds

In the context of drug hypersensitivity, our group has recently proposed a new model based on the structural features of drugs (pharmacological interaction with immune receptors; p-i concept) to explain their recognition by T cells. According to this concept, even chemically inert drugs can stimulate T cells because certain drugs interact in a direct way with T-cell receptors (TCR) and possibly major histocompatibility complex molecules without the need for metabolism and covalent binding to a carrier. In this study, we investigated whether mouse T-cell hybridomas transfected with drug-specific human TCR can be used as an alternative to drug-specific T-cell clones (TCC). Indeed, they behaved like TCC and, in accordance with the p-i concept, the TCR recognize their specific drugs in a direct, processing-independent, and dose-dependent way. The presence of antigen-presenting cells was a prerequisite for interleukin-2 production by the TCR-transfected cells. The analysis of cross-reactivity confirmed the fine specificity of the TCR and also showed that TCR transfectants might provide a tool to evaluate the potential of new drugs to cause hypersensitivity due to cross-reactivity. Recombining the alpha- and beta-chains of sulfanilamide- and quinolone-specific TCR abrogated drug reactivity, suggesting that both original alpha- and beta-chains were involved in drug binding. The TCR-transfected hybridoma system showed that the recognition of two important classes of drugs (sulfanilamides and quinolones) by TCR occurred according to the p-i concept and provides an interesting tool to study drug-TCR interactions and their biological consequences and to evaluate the cross-reactivity potential of new drugs of the same class.

AN ANALYSIS OF THE DNA DAMAGE INDUCED BY NICKEL COMPOUNDS IN CHINESE HAMSTER OVARY CELLS (CARCINOGENESIS, HETEROCHROMATIN, CYTOTOXICITY, REPLICATION, PROTEIN CROSSLINKING)

The carcinogenic activity of water-insoluble crystalline nickel sulfide requires phagocytosis and lysosome-mediated intracellular dissolution of the particles to yield Ni('2+). This study investigated the extent and nature of the DNA damage in Chinese hamster ovary cells treated with various nickel compounds using the technique of alkaline elution. Crystalline NiS and water-soluble NiCl(,2) induced single strand breaks that were repaired quickly and DNA-protein crosslinks that persisted up to 24 hr after exposure to nickel. The induction of single strand breaks was concentration dependent at both noncytotoxic and lethal amounts of nickel. The induction of DNA-protein crosslinks was concentration dependent but was absent at lethal amounts of nickel. The cytoplasmic and nuclear uptake of nickel was concentration dependent even at the toxic level of nickel. However, the induction of DNA-protein crosslinks by nickel required active cell cycling and occurred predominantly in mid-late S phase of the cell cycle, suggesting that the lethal amounts of nickel inhibited DNA-protein crosslinking by inhibiting active cell cycling. Since the DNA-protein crosslinking induced by nickel was resistant to DNA repair, the nature of this lesion was investigated using various methods of DNA isolation and chromatin fractionation in combination with SDS-polyacrylamide gel electrophoresis. High molecular weight, non-histone chromosomal proteins and possibly histone 1 were preferentially crosslinked to DNA by nickel. The crosslinked proteins were concentrated in a magnesium-insoluble fraction of sonicated chromatin (5% of the total) that was similar to heterochromatin in solubility and protein composition. Alterations in DNA structure and function, brought about by the effect of nickel on protein-DNA interactions, may be related to the carcinogenicity of nickel compounds. ^

Searching for the ideal DNA recognition code of covalent ligands

The combitiatorial approach restriction endonuclease protection selection and amplification REPSA was successfully used to determine ideal DNA interactions sites of covalent ligands. Unlike most other combinatorial methods, REPSA is based on inhibition of enzymatic cleavage by specific ligand-DNA complexes, which enables identification of binding sites of various ligands. However, the inherent nature of this technique posses a problem during selection of binding sites of covalent ligands. By modifying the technique according to the nature of the ligand, we demonstrate the flexibility of REPSA in identifying the preferred binding sites for monocovalent ligands, topoisomerase I and tallimustine, and the bicovalent ligand topoisomerase II. From among the preferred binding sites, we identified the consensus binding sequence of camptothecin induced topoisomerase I cleavage as ‘aGWT/Gc’, and tallimustine consensus sequences as ‘GTTCTA’ and ‘TTTTTTC’. We have shown for the first time that preferential binding of tallimustine occurs at sequences not previously reported. Furthermore, our data indicate that tallimustine is a novel DNA minor groove, guanine-specific alkylating agent. ^ Additionally, we have demonstrated in vivo that sequence-specific covalent DNA-binding small molecules have the ability to regulate transcription by inhibiting RNA polymerase II. Tallimustine, binding to its preferred sequences located in the 5′ untranslated region were an effective impediment for transcribing polymerase II. The ability of covalent binding small molecules to target predetermined DNA sequences located downstream of the promoter suggests a general approach for regulation of gene expression. ^

A structural model for a homeotic protein-extradenticle-DNA complex accounts for the choice of HOX protein in the heterodimer.

The genes of the homeotic complex (HOX) encode DNA binding homeodomain proteins that control developmental fates by differentially regulating the transcription of downstream target genes. Despite their unique in vivo functions, disparate HOX proteins often bind to very similar DNA sequences in vitro. Thus, a critical question is how HOX proteins select the correct sets of target genes in vivo. The homeodomain proteins encoded by the Drosophila extradenticle gene and its mammalian homologues, the pbx genes, contribute to HOX specificity by cooperatively binding to DNA with HOX proteins. For example, the HOX protein labial cooperatively binds with extradenticle protein to a 20-bp oligonucleotide that is sufficient to direct a labial-like expression pattern in Drosophila embryos. Here we have analyzed the protein-DNA interactions that are important for forming the labial-extradenticle-DNA complex. The data suggest a model in which labial and extradenticle, separated by only 4 bp, bind this DNA as a heterodimer in a head-to-tail orientation. We have confirmed several aspects of this model by characterizing extradenticle-HOX binding to mutant oligonucleotides. Most importantly, mutations in base pairs predicted to contact the HOX N-terminal arm resulted in a change in HOX preference in the heterodimer, from labial to Ultrabithorax. These results demonstrate that extradenticle prefers to bind cooperatively with different HOX proteins depending on subtle differences in the heterodimer binding site.

Sequence specificity of triplex DNA formation: Analysis by a combinatorial approach, restriction endonuclease protection selection and amplification.

We have devised a combinatorial method, restriction endonuclease protection selection and amplification (REPSA), to identify consensus ligand binding sequences in DNA. In this technique, cleavage by a type IIS restriction endonuclease (an enzyme that cleaves DNA at a site distal from its recognition sequence) is prevented by a bound ligand while unbound DNA is cleaved. Since the selection step of REPSA is performed in solution under mild conditions, this approach is amenable to the investigation of ligand-DNA complexes that are either insufficiently stable or not readily separable by other methods. Here we report the use of REPSA to identify the consensus duplex DNA sequence recognized by a G/T-rich oligodeoxyribonucleotide under conditions favoring purine-motif triple-helix formation. Analysis of 47 sequences indicated that recognition between 13 bases on the oligonucleotide 3' end and the duplex DNA was sufficient for triplex formation and indicated the possible existence of a new base triplet, G.AT. This information should help identify appropriate target sequences for purine-motif triplex formation and demonstrates the power of REPSA for investigating ligand-DNA interactions.

DNA loops induced by cooperative binding of transcriptional activator proteins and preinitiation complexes.

DNA conformational changes are essential for the assembly of multiprotein complexes that contact several DNA sequence elements. An approach based on atomic force microscopy was chosen to visualize specific protein-DNA interactions occurring on eukaryotic class II nuclear gene promoters. Here we report that binding of the transcription regulatory protein Jun to linearized plasmid DNA containing the consensus AP-1 binding site upstream of a class II gene promoter leads to bending of the DNA template. This binding of Jun was found to be essential for the formation of preinitiation complexes (PICs). The cooperative binding of Jun and PIC led to looping of DNA at the protein binding sites. These loops were not seen in the absence of either PICs, Jun, or the AP-1 binding site, suggesting a direct interaction between DNA-bound Jun homodimers and proteins bound to the core promoter. This direct visualization of functional transcriptional complexes confirms the theoretical predictions for the mode of gene regulation by trans-activating proteins.

GANN: Genetic algorithm neural networks for the detection of conserved combinations of features in DNA

Background: The multitude of motif detection algorithms developed to date have largely focused on the detection of patterns in primary sequence. Since sequence-dependent DNA structure and flexibility may also play a role in protein-DNA interactions, the simultaneous exploration of sequence-and structure-based hypotheses about the composition of binding sites and the ordering of features in a regulatory region should be considered as well. The consideration of structural features requires the development of new detection tools that can deal with data types other than primary sequence. Results: GANN ( available at http://bioinformatics.org.au/gann) is a machine learning tool for the detection of conserved features in DNA. The software suite contains programs to extract different regions of genomic DNA from flat files and convert these sequences to indices that reflect sequence and structural composition or the presence of specific protein binding sites. The machine learning component allows the classification of different types of sequences based on subsamples of these indices, and can identify the best combinations of indices and machine learning architecture for sequence discrimination. Another key feature of GANN is the replicated splitting of data into training and test sets, and the implementation of negative controls. In validation experiments, GANN successfully merged important sequence and structural features to yield good predictive models for synthetic and real regulatory regions. Conclusion: GANN is a flexible tool that can search through large sets of sequence and structural feature combinations to identify those that best characterize a set of sequences.

Discovering properties of new DNA-binding activity of proteins

Protein-DNA interactions are an essential feature in the genetic activities of life, and the ability to predict and manipulate such interactions has applications in a wide range of fields. This Thesis presents the methods of modelling the properties of protein-DNA interactions. In particular, it investigates the methods of visualising and predicting the specificity of DNA-binding Cys2His2 zinc finger interaction. The Cys2His2 zinc finger proteins interact via their individual fingers to base pair subsites on the target DNA. Four key residue positions on the a- helix of the zinc fingers make non-covalent interactions with the DNA with sequence specificity. Mutating these key residues generates combinatorial possibilities that could potentially bind to any DNA segment of interest. Many attempts have been made to predict the binding interaction using structural and chemical information, but with only limited success. The most important contribution of the thesis is that the developed model allows for the binding properties of a given protein-DNA binding to be visualised in relation to other protein-DNA combinations without having to explicitly physically model the specific protein molecule and specific DNA sequence. To prove this, various databases were generated, including a synthetic database which includes all possible combinations of the DNA-binding Cys2His2 zinc finger interactions. NeuroScale, a topographic visualisation technique, is exploited to represent the geometric structures of the protein-DNA interactions by measuring dissimilarity between the data points. In order to verify the effect of visualisation on understanding the binding properties of the DNA-binding Cys2His2 zinc finger interaction, various prediction models are constructed by using both the high dimensional original data and the represented data in low dimensional feature space. Finally, novel data sets are studied through the selected visualisation models based on the experimental DNA-zinc finger protein database. The result of the NeuroScale projection shows that different dissimilarity representations give distinctive structural groupings, but clustering in biologically-interesting ways. This method can be used to forecast the physiochemical properties of the novel proteins which may be beneficial for therapeutic purposes involving genome targeting in general.

PDNAsite:identification of DNA-binding site from protein sequence by incorporating spatial and sequence context

Protein-DNA interactions are involved in many fundamental biological processes essential for cellular function. Most of the existing computational approaches employed only the sequence context of the target residue for its prediction. In the present study, for each target residue, we applied both the spatial context and the sequence context to construct the feature space. Subsequently, Latent Semantic Analysis (LSA) was applied to remove the redundancies in the feature space. Finally, a predictor (PDNAsite) was developed through the integration of the support vector machines (SVM) classifier and ensemble learning. Results on the PDNA-62 and the PDNA-224 datasets demonstrate that features extracted from spatial context provide more information than those from sequence context and the combination of them gives more performance gain. An analysis of the number of binding sites in the spatial context of the target site indicates that the interactions between binding sites next to each other are important for protein-DNA recognition and their binding ability. The comparison between our proposed PDNAsite method and the existing methods indicate that PDNAsite outperforms most of the existing methods and is a useful tool for DNA-binding site identification. A web-server of our predictor (http://hlt.hitsz.edu.cn:8080/PDNAsite/) is made available for free public accessible to the biological research community.

Interaction of Multiwalled Carbon Nanotube Produces Structural Alteration of DNA

Multiwalled carbon nanotube (MWCNT) has been found to produce structural changes in Calf Thymus-DNA (CT-DNA). The interaction or binding of the multi-walled carbon nanotubes (MWCNT) was investigated in order to discover if it brings about any significant changes of the DNA double helix using CD spectra of the CT-DNA at two concentration levels of MWCNT representing an increasing MWCNT/DNA molar ratio. In addition, spectrophotometric titrations between MWCNT and CT-DNA were carried out in order to utilize spectral changes as a means of detecting specific binding modes of either intercalation or degradation of DNA. Interactions of MWCNT induced significant changes in the CD spectra of the B-form of natural DNA. The intensities of the positive CD band at 280 nm decreased significantly. This decrease was found to be concentration-dependent. Following spectrophotometric titrations; specific subtle conformational changes were observed with a molar ratio combination of 2:1 between MWCNT and CT-DNA and these were characterized by a formation constant of the order of 103 M-1 and a negative Gibbs free energy suggesting that MWCNT avidly binds to DNA. Thermodynamic considerations revealed that electrostatic interactions between the DNA base pairs and the MWCNT are taking place accounting for the negative free energy change, positive enthalpy change with a small entropy change. The results obtained in the study of the binding interactions of MWCNT with DNA confirm that a cytogenetic effect of MWCNT with DNA is a possibility in vivo.