Dataset Information

Automated incorporation of pairwise dependency in transcription factor binding site prediction using dinucleotide weight tensors.

ABSTRACT: Gene regulatory networks are ultimately encoded by the sequence-specific binding of (TFs) to short DNA segments. Although it is customary to represent the binding specificity of a TF by a position-specific weight matrix (PSWM), which assumes each position within a site contributes independently to the overall binding affinity, evidence has been accumulating that there can be significant dependencies between positions. Unfortunately, methodological challenges have so far hindered the development of a practical and generally-accepted extension of the PSWM model. On the one hand, simple models that only consider dependencies between nearest-neighbor positions are easy to use in practice, but fail to account for the distal dependencies that are observed in the data. On the other hand, models that allow for arbitrary dependencies are prone to overfitting, requiring regularization schemes that are difficult to use in practice for non-experts. Here we present a new regulatory motif model, called dinucleotide weight tensor (DWT), that incorporates arbitrary pairwise dependencies between positions in binding sites, rigorously from first principles, and free from tunable parameters. We demonstrate the power of the method on a large set of ChIP-seq data-sets, showing that DWTs outperform both PSWMs and motif models that only incorporate nearest-neighbor dependencies. We also demonstrate that DWTs outperform two previously proposed methods. Finally, we show that DWTs inferred from ChIP-seq data also outperform PSWMs on HT-SELEX data for the same TF, suggesting that DWTs capture inherent biophysical properties of the interactions between the DNA binding domains of TFs and their binding sites. We make a suite of DWT tools available at dwt.unibas.ch, that allow users to automatically perform 'motif finding', i.e. the inference of DWT motifs from a set of sequences, binding site prediction with DWTs, and visualization of DWT 'dilogo' motifs.

SUBMITTER: Omidi S

PROVIDER: S-EPMC5550003 | biostudies-literature | 2017 Jul

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Automated incorporation of pairwise dependency in transcription factor binding site prediction using dinucleotide weight tensors.

Omidi Saeed S Zavolan Mihaela M Pachkov Mikhail M Breda Jeremie J Berger Severin S van Nimwegen Erik E

PLoS computational biology 20170728 7

Gene regulatory networks are ultimately encoded by the sequence-specific binding of (TFs) to short DNA segments. Although it is customary to represent the binding specificity of a TF by a position-specific weight matrix (PSWM), which assumes each position within a site contributes independently to the overall binding affinity, evidence has been accumulating that there can be significant dependencies between positions. Unfortunately, methodological challenges have so far hindered the development ...[more]

PMID: 28753602

Similar Datasets

Project description:Finding where transcription factors (TFs) bind to the DNA is of key importance to decipher gene regulation at a transcriptional level. Classically, computational prediction of TF binding sites (TFBSs) is based on basic position weight matrices (PWMs) which quantitatively score binding motifs based on the observed nucleotide patterns in a set of TFBSs for the corresponding TF. Such models make the strong assumption that each nucleotide participates independently in the corresponding DNA-protein interaction and do not account for flexible length motifs. We introduce transcription factor flexible models (TFFMs) to represent TF binding properties. Based on hidden Markov models, TFFMs are flexible, and can model both position interdependence within TFBSs and variable length motifs within a single dedicated framework. The availability of thousands of experimentally validated DNA-TF interaction sequences from ChIP-seq allows for the generation of models that perform as well as PWMs for stereotypical TFs and can improve performance for TFs with flexible binding characteristics. We present a new graphical representation of the motifs that convey properties of position interdependence. TFFMs have been assessed on ChIP-seq data sets coming from the ENCODE project, revealing that they can perform better than both PWMs and the dinucleotide weight matrix extension in discriminating ChIP-seq from background sequences. Under the assumption that ChIP-seq signal values are correlated with the affinity of the TF-DNA binding, we find that TFFM scores correlate with ChIP-seq peak signals. Moreover, using available TF-DNA affinity measurements for the Max TF, we demonstrate that TFFMs constructed from ChIP-seq data correlate with published experimentally measured DNA-binding affinities. Finally, TFFMs allow for the straightforward computation of an integrated TF occupancy score across a sequence. These results demonstrate the capacity of TFFMs to accurately model DNA-protein interactions, while providing a single unified framework suitable for the next generation of TFBS prediction.

Project description:BackgroundEpitope identification is an essential step toward synthetic vaccine development since epitopes play an important role in activating immune response. Classical experimental approaches are laborious and time-consuming, and therefore computational methods for generating epitope candidates have been actively studied. Most of these methods, however, are based on sophisticated nonlinear techniques for achieving higher predictive performance. The use of these techniques tend to diminish their interpretability with respect to binding potential: that is, they do not provide much insight into binding mechanisms.ResultsWe have developed a novel epitope prediction method named EpicCapo and its variants, EpicCapo(+) and EpicCapo(+REF). Nonapeptides were encoded numerically using a novel peptide-encoding scheme for machine learning algorithms by utilizing 40 amino acid pairwise contact potentials (referred to as AAPPs throughout this paper). The predictive performances of EpicCapo(+) and EpicCapo(+REF) outperformed other state-of-the-art methods without losing interpretability. Interestingly, the most informative AAPPs estimated by our study were those developed by Micheletti and Simons while previous studies utilized two AAPPs developed by Miyazawa & Jernigan and Betancourt & Thirumalai. In addition, we found that all amino acid positions in nonapeptides could effect on performances of the predictive models including non-anchor positions. Finally, EpicCapo(+REF) was applied to identify candidates of promiscuous epitopes. As a result, 67.1% of the predicted nonapeptides epitopes were consistent with preceding studies based on immunological experiments.ConclusionsOur method achieved high performance in testing with benchmark datasets. In addition, our study identified a number of candidates of promiscuous CTL epitopes consistent with previously reported immunological experiments. We speculate that our techniques may be useful in the development of new vaccines. The R implementation of EpicCapo(+REF) is available at http://pirun.ku.ac.th/~fsciiok/EpicCapoREF.zip. Datasets are available at http://pirun.ku.ac.th/~fsciiok/Datasets.zip.

Project description:One of the most straightforward and commonly used chemical modifications of proteins is to react surface amino groups (lysine residues) with activated esters. This chemistry has been used to generate protein-polymer conjugates, many of which are now approved therapeutics. Similar conjugates have also been generated by reacting activated ester atom transfer polymerization initiators with lysine residues to create biomacromolecular initiators for polymerization reactions. The reaction between activated esters and lysine amino groups is rapid and has been consistently described in almost every publication on the topic as a "random reaction". A random reaction implies that every accessible lysine amino group on a protein molecule is equally reactive, and as a result, that the reaction is indiscriminate. Nonetheless, the literature contradicts itself by also suggesting that some lysine amino groups are more reactive than others (as a function of pKa, surface accessibility, temperature, and local environment). If the latter assumption is correct, then the outcome of these reactions cannot be random at all, and we should be able to predict the outcome from the structure of the protein. Predicting the non-random outcome of a reaction between surface lysines and reactive esters could transform the speed at which active bioconjugates can be developed and engineered. Herein, we describe a robust integrated tool that predicts the activated ester reactivity of every lysine in a protein, thereby allowing us to calculate the non-random sequence of reaction as a function of reaction conditions. Specifically, we have predicted the intrinsic reactivity of each lysine in multiple proteins with a bromine-functionalised N-hydroxysuccinimide initiator molecule. We have also shown that the model applied to PEGylation. The rules-based analysis has been coupled together in a single Python program that can bypass tedious trial and error experiments usually needed in protein-polymer conjugate design and synthesis.

Dataset Information

Automated incorporation of pairwise dependency in transcription factor binding site prediction using dinucleotide weight tensors.

Publications

Automated incorporation of pairwise dependency in transcription factor binding site prediction using dinucleotide weight tensors.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets