DOCSLIB.ORG

Regulatory Motifs of DNA

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
Regulatory Motifs of DNA Lecture 3: PWMs and other models for regulatory elements in DNA -Sequence motifs -Position weight matrices (PWMs) -Learning PWMs combinatorially, or probabilistically -Learning from an alignment -Ab initio learning: Baum-Welch & Gibbs sampling -Visualization of PWMs as sequence logos -Search methods for PWM occurrences -Cis-regulatory modules Regulatory motifs of DNA Measuring gene activity (’expression’) with a microarray 1 Jointly regulated gene sets Cluster of co-expressed genes, apply pattern discovery in regulatory regions 600 basepairs Expression profiles Retrieve Upstream regions Genes with similar expression pattern Pattern over-represented in cluster Find a hidden motif 2 Find a hidden motif Find a hidden motif (cont.) Multiple local alignment! Smith-Waterman?? 3 Motif? • Definition: Motif is a pattern that occurs (unexpectedly) often in a string or in a set of strings • The problem of finding repetitions is similar cgccgagtgacagagacgctaatcagg ctgtgttctcaggatgcgtaccgagtg ggagacagcagcacgaccagcggtggc agagacccttgcagacatcaagctctt tgggaacaagtggagcaccgatgatgt acagccgatcaatgacatttccctaat gcaggattacattgcagtgcccaagga gaagtatgccaagtaatacctccctca cagtg... 4 Longest repeat? Example: A 50 million bases long fragment of human genome. We found (using a suffix-tree) the following repetitive sequence of length 2559, with one occurrence at 28395980, and other at 28401554r ttagggtacatgtgcacaacgtgcaggtttgttacatatgtatacacgtgccatgatggtgtgctgcacccattaactcgtcatttagcgttaggtatatctcc gaatgctatccctcccccctccccccaccccacaacagtccccggtgtgtgatgttccccttcctgtgtccatgtgttctcattgttcaattcccacctatgagt gagaacatgcggtgtttggttttttgtccttgcgaaagtttgctgagaatgatggtttccagcttcatccatatccctacaaaggacatgaactcatcatttttt tatggctgcatagtattccatggtgtatatgtgccacattttcttaacccagtctacccttgttggacatctgggttggttccaagtctttgctattgtgaatag tgccgcaataaacatacgtgtgcatgtgtctttatagcagcatgatttataatcctttgggtatatacccagtaatgggatggctgggtcaaatggtatttcta gttctagatccctgaggaatcaccacactgacttccacaatggttgaactagtttacagtcccagcaacagttcctatttctccacatcctctccagcacctgt tgtttcctgactttttaatgatcgccattctaactggtgtgagatggtatctcattgtggttttgatttgcatttctctgatggccagtgatgatgagcatttttt catgtgttttttggctgcataaatgtcttcttttgagaagtgtctgttcatatccttcgcccacttttgatggggttgtttgtttttttcttgtaaatttgttgga gttcattgtagattctgggtattagccctttgtcagatgagtaggttgcaaaaattttctcccattctgtaggttgcctgttcactctgatggtggtttcttctgc tgtgcagaagctctttagtttaattagatcccatttgtcaattttggcttttgttgccatagcttttggtgttttagacatgaagtccttgcccatgcctatgtcc tgaatggtattgcctaggttttcttctagggtttttatggttttaggtctaacatgtaagtctttaatccatcttgaattaattataaggtgtatattataaggt gtaattataaggtgtataattatatattaattataaggtgtatattaattataaggtgtaaggaagggatccagtttcagctttctacatatggctagccagtt ttccctgcaccatttattaaatagggaatcctttccccattgcttgtttttgtcaggtttgtcaaagatcagatagttgtagatatgcggcattatttctgaggg ctctgttctgttccattggtctatatctctgttttggtaccagtaccatgctgttttggttactgtagccttgtagtatagtttgaagtcaggtagcgtgatggtt ccagctttgttcttttggcttaggattgacttggcaatgtgggctcttttttggttccatatgaactttaaagtagttttttccaattctgtgaagaaattcattg gtagcttgatggggatggcattgaatctataaattaccctgggcagtatggccattttcacaatattgaatcttcctacccatgagcgtgtactgttcttccatt tgtttgtatcctcttttatttcattgagcagtggtttgtagttctccttgaagaggtccttcacatcccttgtaagttggattcctaggtattttattctctttga agcaattgtgaatgggagttcactcatgatttgactctctgtttgtctgttattggtgtataagaatgcttgtgatttttgcacattgattttgtatcctgagac tttgctgaagttgcttatcagcttaaggagattttgggctgagacgatggggttttctagatatacaatcatgtcatctgcaaacagggacaatttgacttcct cttttcctaattgaatacccgttatttccctctcctgcctgattgccctggccagaacttccaacactatgttgaataggagtggtgagagagggcatccctgt cttgtgccagttttcaaagggaatgcttccagtttttgtccattcagtatgatattggctgtgggtttgtcatagatagctcttattattttgagatacatccca tcaatacctaatttattgagagtttttagcatgaagagttcttgaattttgtcaaaggccttttctgcatcttttgagataatcatgtggtttctgtctttggttc tgtttatatgctggagtacgtttattgattttcgtatgttgaaccagccttgcatcccagggatgaagcccacttgatcatggtggataagctttttgatgtgct gctggattcggtttgccagtattttattgaggatttctgcatcgatgttcatcaaggatattggtctaaaattctctttttttgttgtgtctctgtcaggctttgg tatcaggatgatgctggcctcataaaatgagttagg Representations of motifs • pattern – substring – substring with gaps – string in generalized alphabet (e.g., IUPAC) – finite automaton – Hidden Markov Model – binding affinity matrix – cluster of binding affinity matrices – … (= the hidden structure to be learned from data) 5 Types of occurrences • occurrence – exact – approximate – with high probability –… Subsequence motifs with approximate occurrences 6 Motif finding problem • Given: a set of sequences S = x(1),x(2),...,x(n) from alphabet Ȉ, and a motif length m • Find the best motif of length m that occurs in the training data S • Best? Occurs? • Combinatorial approach: motif = sequence of length m – Approximate occurrences defined using Hamming distance (= number of mismatches) • Probabilistic approach: motif = PWM (position weight matrix) of length m – Approximate occurrences defined using the probability distribution Combinatorial approach • d(W,x(i)) = minimum Hamming distance between W and any subsequence of x(i) of length m (i) • d(W,S) = i d(W,x ) 7 Combinatorial approach: a pattern- driven algorithm 1. For all sequences W İȈm of length m do Find d(W,S) 2. Report W* := arg minW (d(W,S)) For DNA, the trivial implementation has running time O(mN4m). Why? Explain! (N = total length of the sequences in S) Combinatorial approach: a sample- driven algorithm 1. For all subsequences W of length m from S do Find d(W,S) 2. Report W* := arg minW(d(W,S)) For DNA, the trivial implementation has running time O(mN2). Why? Explain! (N = total length of sequences in S) 8 Probabilistic approach • Learn Position Weight Matrices (PWMs) from data • Probabilistic model (in fact, a simple HMM) Positionally weighted pattern • Weighted pattern w = (wij) of length m in alphabet Ȉ: _Ȉ| x m matrix of real-valued scores • Also called as: position weight matrix (PWM), position- specific scoring matrix (PSSM), profile(-HMM), motif, ... (Stormo et al 1980, Gribskov et al 1987, Henikoff et al 1990,...) • Public collections of PWMs: TRANSFAC, JASPAR 1 2 3 4 5 6 A 0.3 0.0 0.1 0.2 1.0 0.3 C 0.1 0.8 0.5 0.2 0.0 0.4 G 0.2 0.0 0.4 0.3 0.0 0.0 T 0.4 0.2 0.0 0.3 0.0 0.3 9 Construction of PWM from an alignment • Given: – aligned set of binding sites (= sequences of length m) – background distribution qa for a ɽ Ȉ • PWM construction algorithm: 1. Construct count matrix (Nij): Nij := the number of symbols i ɽ Ȉ on column j of the aligment 2. Add pseudocounts: • Nij := Nij + 1 (= Laplace rule), or • Nij := Nij + ȕqi where ȕ is a scaling parameter that determines the total number of pseudocounts in an alignment column 3. Construct probability matrix (fij)by normalizing the (pseudo)counts: fij := Nij / i Nij 4. Construct log-odds of signal vs background: wij := log(fij / qi) (= weighted pattern w) PWM construction: example 1 2 3 4 5 6 7 CTCACACGTGGG (Nij): A 2 0 5 0 0 0 0 TCACACGTGGGA C 0 5 0 5 0 0 0 ATTAGCACGTTTT G 3 0 0 0 5 1 3 TTAGCACGTTTCGC T 0 0 0 0 0 4 2 CGCTGCACGGGGCC 10 PWM construction: example (cont.) Step 1: (Nij) Normalization without pseudocounts A 2 0 5 0 0 0 0 A 0.4 0 1.0 0 0 0 0 C 0 5 0 5 0 0 0 C 0 1.0 0 1.0 0 0 0 G 3 0 0 0 5 1 3 G 0.6 0 0 0 1.0 0.2 0.6 T 0 0 0 0 0 4 2 T 0 0 0 0 0 0.8 0.4 Nij := Nij + ȕqi where qA = qT = 0.32, qG = qC = 0.18 ȕ = 1 Step 2: (Nij) with pseudocounts A 2.32 0.32 5.32 0.32 0.32 0.32 0.32 C 0.18 5.18 0.18 5.18 0.18 0.18 0.18 G 3.18 0.18 0.18 0.18 5.18 1.18 3.18 T 0.32 0.32 0.32 0.32 0.32 4.32 2.32 Step 4: log-odds wij := log2 (fij / qi ) Step 3: probability matrix (f ) := (N ) / 6 ij ij A ... -2.7 1.5 -2.7 -2.7 ... 0.39 0.05 0.89 0.05 0.05 0.05 0.05 A C 2.3 -2.6 2.3 -2.6 C 0.03 0.87 0.03 0.87 0.03 0.03 0.03 G -2.6 -2.6 -2.6 2.3 G 0.53 0.03 0.03 0.03 0.87 0.20 0.53 T ... -2.7 -2.7 -2.7 -2.7 ... T 0.05 0.05 0.05 0.05 0.05 0.72 0.39 Visualisation of PWMs: sequence logo • PWM probability matrix (faj): faj = probability of symbol a in motif position j f I • Entropy of column j: A Aj j Hj = -a ɽȈ faj log2 faj Ij T C • Information present in column j: G Ij = log2 _Ȉ| - Hj [DNA: Ij 2] • Logo: Column j has total height Ij, symbol a has height faj·Ij 11 Sequence logo: an example A 9 11 49 51 0 1 1 4 C 19 3 0 0 0 45 25 16 G 5 1 2 0 17 0 4 21 T 18 36 0 0 34 5 21 10 Ab initio construction of PWMs • No alignment given, just the training sequences S = x(1), ..., x(N) • Find from S all m-long words • Assume each word comes either from motif or from background (don’t care about overlaps!) • Find the most likely – motif model (PWM), – background model, and – classification of the m-long words of S into motif and background instances 12 HMM for mixture of multinomials PWM ... Ȝ M1 M2 M3 Mm 1-Ȝ Q Background m-1 times PWM by the EM algorithm • Train the mixture model using EM algorithm by iterating – E step – M step • Training data: all m-words from S • For more details, see: – T.L. Bailey & C. Elkan: The value of prior knowledge in discovering motifs with MEME – Other papers on MEME 13 Ab initio motif finding: Gibbs sampling • Another popular probabilistic algorithm for motif (PWM) discovery • Local search algorithm Gibbs sampling: basic idea Current motif = PWM formed by circled substrings 14 Gibbs sampling: basic idea Delete one substring Gibbs sampling: basic idea Try a replacement: Compute its score, and accept the replacement depending on the score. 15 Gibbs sampling: basic idea New current motif Repeat many times! PWM by Gibbs sampling • Goal: Find most probable pattern by sampling from motif probabilities to maximize the ratio of model to background • Given: (1) (N) – Training data S = X , …, X , and background distribution (qa) – motif length m • Notation & algorithm idea: – Model
Load more

Recommended publications
  • Motif Discovery in Sequential Data Kyle L. Jensen
    Motif discovery in sequential data by Kyle L. Jensen Submitted to the Department of Chemical Engineering in partial fulfillment of the requirements for the degree of Doctor of Philosophy at the MASSACHUSETTS INSTITUTE OF TECHNOLOGY May © Kyle L. Jensen, . Author.............................................................. Department of Chemical Engineering April , Certified by . Gregory N. Stephanopoulos Bayer Professor of Chemical Engineering esis Supervisor Accepted by . William M. Deen Chairman, Department Committee on Graduate Students Motif discovery in sequential data by Kyle L. Jensen Submitted to the Department of Chemical Engineering on April , , in partial fulfillment of the requirements for the degree of Doctor of Philosophy Abstract In this thesis, I discuss the application and development of methods for the automated discovery of motifs in sequential data. ese data include DNA sequences, protein sequences, and real– valued sequential data such as protein structures and timeseries of arbitrary dimension. As more genomes are sequenced and annotated, the need for automated, computational methods for analyzing biological data is increasing rapidly. In broad terms, the goal of this thesis is to treat sequential data sets as unknown languages and to develop tools for interpreting an understanding these languages. e first chapter of this thesis is an introduction to the fundamentals of motif discovery, which establishes a common mode of thought and vocabulary for the subsequent chapters. One of the central themes of this work is the use of grammatical models, which are more commonly associated with the field of computational linguistics. In the second chapter, I use grammatical models to design novel antimicrobial peptides (AmPs). AmPs are small proteins used by the innate immune system to combat bacterial infection in multicellular eukaryotes.
    [Show full text]
  • A General Pairwise Interaction Model Provides an Accurate Description of in Vivo Transcription Factor Binding Sites
    A General Pairwise Interaction Model Provides an Accurate Description of In Vivo Transcription Factor Binding Sites Marc Santolini, Thierry Mora, Vincent Hakim* Laboratoire de Physique Statistique, CNRS, Universite´ P. et M. Curie, Universite´ D. Diderot, E´cole Normale Supe´rieure, Paris, France Abstract The identification of transcription factor binding sites (TFBSs) on genomic DNA is of crucial importance for understanding and predicting regulatory elements in gene networks. TFBS motifs are commonly described by Position Weight Matrices (PWMs), in which each DNA base pair contributes independently to the transcription factor (TF) binding. However, this description ignores correlations between nucleotides at different positions, and is generally inaccurate: analysing fly and mouse in vivo ChIPseq data, we show that in most cases the PWM model fails to reproduce the observed statistics of TFBSs. To overcome this issue, we introduce the pairwise interaction model (PIM), a generalization of the PWM model. The model is based on the principle of maximum entropy and explicitly describes pairwise correlations between nucleotides at different positions, while being otherwise as unconstrained as possible. It is mathematically equivalent to considering a TF-DNA binding energy that depends additively on each nucleotide identity at all positions in the TFBS, like the PWM model, but also additively on pairs of nucleotides. We find that the PIM significantly improves over the PWM model, and even provides an optimal description of TFBS statistics within statistical noise. The PIM generalizes previous approaches to interdependent positions: it accounts for co-variation of two or more base pairs, and predicts secondary motifs, while outperforming multiple-motif models consisting of mixtures of PWMs.
    [Show full text]
  • A New Sequence Logo Plot to Highlight Enrichment and Depletion
    bioRxiv preprint doi: https://doi.org/10.1101/226597; this version posted November 29, 2017. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY-NC-ND 4.0 International license. A new sequence logo plot to highlight enrichment and depletion. Kushal K. Dey 1, Dongyue Xie 1, Matthew Stephens 1, 2 1 Department of Statistics, University of Chicago 2 Department of Human Genetics, University of Chicago * Corresponding Email : [email protected] Abstract Background : Sequence logo plots have become a standard graphical tool for visualizing sequence motifs in DNA, RNA or protein sequences. However standard logo plots primarily highlight enrichment of symbols, and may fail to highlight interesting depletions. Current alternatives that try to highlight depletion often produce visually cluttered logos. Results : We introduce a new sequence logo plot, the EDLogo plot, that highlights both enrichment and depletion, while minimizing visual clutter. We provide an easy-to-use and highly customizable R package Logolas to produce a range of logo plots, including EDLogo plots. This software also allows elements in the logo plot to be strings of characters, rather than a single character, extending the range of applications beyond the usual DNA, RNA or protein sequences. We illustrate our methods and software on applications to transcription factor binding site motifs, protein sequence alignments and cancer mutation signature profiles. Conclusion : Our new EDLogo plots, and flexible software implementation, can help data analysts visualize both enrichment and depletion of characters (DNA sequence bases, amino acids, etc) across a wide range of applications.
    [Show full text]
  • Predicting Expression Levels of ​De Novo​Protein Designs in Yeast
    Predicting expression levels of de novo protein designs in yeast through machine learning. ​ ​ Patrick Slogic A thesis submitted in partial fulfilment of the requirements for the degree of: Master of Science in Bioengineering University of Washington 2020 Committee: Eric Klavins Azadeh Yazdan-Shahmorad Program Authorized to Offer Degree: Bioengineering ©Copyright 2020 Patrick Slogic University of Washington Abstract Predicting expression levels of de novo protein designs in yeast through machine learning. ​ ​ Patrick Slogic Chair of the Supervisory Committee: Eric Klavins Department of Electrical and Computer Engineering Protein engineering has unlocked potential for the biomolecular synthesis of amino acid sequences with specified functions. With the many implications in enhancing therapy development, developing accurate functional models of recombinant genetic outcomes has shown great promise in unlocking the potential of novel protein structures. However, the complicated nature of structural determination requires an iterative process of design and testing among an enormous search space for potential designs. Expressed protein sequences through biological platforms are necessary to produce and analyze the functional efficacy of novel designs. In this study, machine learning principles are applied to a yeast surface display protocol with the goal of optimizing the search space of amino acid sequences that will express in yeast transformations. Machine learning networks are compared in predicting expression levels and generating similar sequences to find those likely to be successfully expressed. ACKNOWLEDGEMENTS I would like to thank Dr. Eric Klavins for accepting me as a master’s student in his laboratory and for providing direction and guidance in exploring the topics of machine learning and recombinant genetics throughout my graduate program in bioengineering.
    [Show full text]
  • Low-N Protein Engineering with Data-Efficient Deep Learning
    bioRxiv preprint doi: https://doi.org/10.1101/2020.01.23.917682; this version posted January 24, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY-NC-ND 4.0 International license. Low-N protein engineering with data-efficient deep learning 1† 3† 2† 2 1,4 Surojit Biswas ,​ Grigory Khimulya ,​ Ethan C. Alley ,​ Kevin M. Esvelt ,​ George M. Church ​ ​ ​ ​ ​ 1 Wyss Institute for Biologically Inspired Engineering, Harvard University 2 MIT Media Lab, Massachusetts Institute of Technology 3 Somerville, MA 02143, USA 4 Department of Genetics, Harvard Medical School †These authors contributed equally to this work. *Correspondence to: George M. Church: gchurch{at}genetics.med.harvard.edu ​ Abstract Protein engineering has enormous academic and industrial potential. However, it is limited by the lack of experimental assays that are consistent with the design goal and sufficiently high-throughput to find rare, enhanced variants. Here we introduce a machine learning-guided paradigm that can use as few as 24 functionally assayed mutant sequences to build an accurate virtual fitness landscape and screen ten million sequences via in silico directed evolution. As demonstrated in two highly dissimilar proteins, avGFP and ​ ​ TEM-1 β-lactamase, top candidates from a single round are diverse and as active as engineered mutants ​ ​ obtained from previous multi-year, high-throughput efforts. Because it distills information from both global and local sequence landscapes, our model approximates protein function even before receiving experimental data, and generalizes from only single mutations to propose high-functioning epistatically non-trivial designs.
    [Show full text]
  • Lecture 12. Gene Finding & Rnaseq
    Lecture 12. Gene Finding & RNAseq Michael Schatz March 4, 2020 Applied Comparative Genomics Assignment 4: Due Wed Mar 4 Assignment 5: Due Wed Mar 11 Oxford Nanopore MinION • Thumb drive sized sequencer powered over USB • Contains 512 channels • Four pores per channel, only one pore active at a time • Early access began in 2014 • Officially released in 2015 (the same year I met MiKe!) Assignment 5: Due March 8 Goal: Genome Annotations aatgcatgcggctatgctaatgcatgcggctatgctaagctgggatccgatgacaatgcatgcggctatgctaat gcatgcggctatgcaagctgggatccgatgactatgctaagctgggatccgatgacaatgcatgcggctatgct aatgaatggtcttgggatttaccttggaatgctaagctgggatccgatgacaatgcatgcggctatgctaatgaa tggtcttgggatttaccttggaatatgctaatgcatgcggctatgctaagctgggatccgatgacaatgcatgcg gctatgctaatgcatgcggctatgcaagctgggatccgatgactatgctaagctgcggctatgctaatgcatgcg gctatgctaagctgggatccgatgacaatgcatgcggctatgctaatgcatgcggctatgcaagctgggatcct gcggctatgctaatgaatggtcttgggatttaccttggaatgctaagctgggatccgatgacaatgcatgcggct atgctaatgaatggtcttgggatttaccttggaatatgctaatgcatgcggctatgctaagctgggaatgcatgcg gctatgctaagctgggatccgatgacaatgcatgcggctatgctaatgcatgcggctatgcaagctgggatccg atgactatgctaagctgcggctatgctaatgcatgcggctatgctaagctcatgcggctatgctaagctgggaat gcatgcggctatgctaagctgggatccgatgacaatgcatgcggctatgctaatgcatgcggctatgcaagctg ggatccgatgactatgctaagctgcggctatgctaatgcatgcggctatgctaagctcggctatgctaatgaatg gtcttgggatttaccttggaatgctaagctgggatccgatgacaatgcatgcggctatgctaatgaatggtcttgg gatttaccttggaatatgctaatgcatgcggctatgctaagctgggaatgcatgcggctatgctaagctgggatc cgatgacaatgcatgcggctatgctaatgcatgcggctatgcaagctgggatccgatgactatgctaagctgcg
    [Show full text]
  • Interpretable Convolution Methods for Learning Genomic Sequence Motifs
    bioRxiv preprint doi: https://doi.org/10.1101/411934; this version posted September 8, 2018. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY-NC 4.0 International license. Interpretable Convolution Methods for Learning Genomic Sequence Motifs MS Ploenzke1*, RA Irizarry1,2 1 Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA, USA 2 Department of Biostatistics and Computational Biology, Dana-Farber Cancer Institute, Boston, MA, USA * Corresponding author: ploenzke(at)g.harvard.edu Abstract The first-layer filters employed in convolutional neural networks tend to learn, or extract, spatial features from the data. Within their application to genomic sequence data, these learned features are often visualized and interpreted by converting them to sequence logos; an information-based representation of the consensus nucleotide motif. The process to obtain such motifs, however, is done through post-training procedures which often discard the filter weights themselves and instead rely upon finding those sequences maximally correlated with the given filter. Moreover, the filters collectively learn motifs with high redundancy, often simply shifted representations of the same sequence. We propose a schema to learn sequence motifs directly through weight constraints and transformations such that the individual weights comprising the filter are directly interpretable as either position weight matrices (PWMs) or information gain matrices (IGMs). We additionally leverage regularization to encourage learning highly-representative motifs with low inter-filter redundancy. Through learning PWMs and IGMs directly we present preliminary results showcasing how our method is capable of incorporating previously-annotated database motifs along with learning motifs de novo and then outline a pipeline for how these tools may be used jointly in a data application.
    [Show full text]
  • Specificity and Nonspecificity in RNA–Protein Interactions
    REVIEWS Specificity and nonspecificity in RNA–protein interactions Eckhard Jankowsky1–3 and Michael E. Harris2 Abstract | To fully understand the regulation of gene expression, it is critical to quantitatively define whether and how RNA‑binding proteins (RBPs) discriminate between alternative binding sites in RNAs. Here, we describe new methods that measure protein binding to large numbers of RNA variants and that reveal [Au: ok? ‘measure.. binding patterns’ did not seem correct] the binding patterns they produce, including affinity distributions and free energy landscapes. We discuss how the new methodologies and the associated concepts enable the development of inclusive, quantitative models for RNA–protein interactions that transcend the traditional binary classification of RBPs as either specific or nonspecific. RNA–protein interactions are critical for the regulation RNA–protein interactions in cells or measured RNA– of gene expression1. Research over recent decades has protein associations for large numbers of sequences shown that RNA is invariably bound and often altered in vitro indicated that specificity, or the lack thereof, is by proteins in cells, and that in biological environments considerably more nuanced than suggested by the binary RNAs generally function together with proteins as specific versus nonspecific classification. As descriptions RNA–protein complexes known as ribonucleoproteins of cellular RNA–protein interaction networks move (RNPs)2,3. It has also become clear that cellular RNA– towards systems-level quantitative models4,11,12, and as protein interactions form a highly complex network other lines of research attempt to engineer novel RNA- involving numerous RNAs and proteins4. In addition, a binding proteins (RBPs)13–16, a comprehensive, quantita- multitude of diseases have been linked to misregulation tive view on specificity and nonspecificity is required.
    [Show full text]
  • Ancient Mechanisms for the Evolution of the Bicoid Homeodomain's
    RESEARCH ARTICLE Ancient mechanisms for the evolution of the bicoid homeodomain’s function in fly development Qinwen Liu1†, Pinar Onal2†, Rhea R Datta2†‡, Julia M Rogers3,4, Urs Schmidt-Ott5, Martha L Bulyk3,4,6, Stephen Small2*, Joseph W Thornton1,7* 1Department of Ecology and Evolution, University of Chicago, Chicago, United States; 2Department of Biology, New York University, New York, United States; 3Committee on Higher Degrees in Biophysics, Harvard University, Cambridge, United States; 4Division of Genetics, Department of Medicine, Brigham and Women’s Hospital and Harvard Medical School, Boston, United States; 5Department of Organismal Biology and Anatomy, University of Chicago, Chicago, United States; 6Department of Pathology, Brigham and Women’s Hospital and Harvard Medical School, Boston, United States; 7Department of Human Genetics, University of Chicago, Chicago, United States Abstract The ancient mechanisms that caused developmental gene regulatory networks to diversify among distantly related taxa are not well understood. Here we use ancestral protein reconstruction, biochemical experiments, and developmental assays of transgenic animals carrying *For correspondence: [email protected] (SS); reconstructed ancestral genes to investigate how the transcription factor Bicoid (Bcd) evolved its [email protected] (JWT) central role in anterior-posterior patterning in flies. We show that most of Bcd’s derived functions are attributable to evolutionary changes within its homeodomain (HD) during a phylogenetic †These authors contributed interval >140 million years ago. A single substitution from this period (Q50K) accounts almost equally to this work entirely for the evolution of Bcd’s derived DNA specificity in vitro. In transgenic embryos Present address: ‡Biology expressing the reconstructed ancestral HD, however, Q50K confers activation of only a few of Department, Hamilton College, Bcd’s transcriptional targets and yields a very partial rescue of anterior development.
    [Show full text]
  • Probability and Statistics for Bioinformatics and Genetics Course
    Probability and Statistics for Bioinformatics and Genetics Course Notes Paul Maiste The Johns Hopkins University These notes are meant solely to be used as study notes for students in 550.435 - Bioinformatics and Statistical Genetics, Spring 2006. c Copyright 2006 Paul Maiste, except for various figures and definitions. May 2, 2006 Contents 1 Introduction to Cell Biology and Genetics 5 1.1 The Cell, DNA, and Chromosomes ................................. 5 1.2 Function of the Cell and DNA .................................... 9 1.3 Genetics ................................................ 13 1.4 Genetic Inheritance .......................................... 19 2 Probability and Statistics Review 28 2.1 Introduction to Random Variables ................................. 28 2.2 Discrete Random Variables ..................................... 30 2.2.1 Probability Distribution ................................... 30 2.2.2 Cumulative Distribution Function ............................. 32 2.3 Continuous Random Variables .................................... 33 2.3.1 General Details About a Continuous Random Variable . 33 2.3.2 Properties of the density function f(x) ........................... 35 2.3.3 Cumulative Distribution Function ............................. 36 2.4 Expected Value and Variance .................................... 37 2.5 Percentiles of a Random Variable .................................. 38 2.6 Common Discrete Distributions ................................... 39 2.6.1 Bernoulli Random Variable ................................
    [Show full text]
  • Predicting Variation of DNA Shape Preferences in Protein-DNA Interaction in Cancer Cells with a New Biophysical Model
    G C A T T A C G G C A T genes Article Predicting Variation of DNA Shape Preferences in Protein-DNA Interaction in Cancer Cells with a New Biophysical Model Kirill Batmanov and Junbai Wang * ID Department of Pathology, Oslo University Hospital—Norwegian Radium Hospital, Montebello, 0310 Oslo, Norway; [email protected] * Correspondence: [email protected] Received: 31 July 2017; Accepted: 13 September 2017; Published: 18 September 2017 Abstract: DNA shape readout is an important mechanism of transcription factor target site recognition, in addition to the sequence readout. Several machine learning-based models of transcription factor–DNA interactions, considering DNA shape features, have been developed in recent years. Here, we present a new biophysical model of protein–DNA interactions by integrating the DNA shape properties. It is based on the neighbor dinucleotide dependency model BayesPI2, where new parameters are restricted to a subspace spanned by the dinucleotide form of DNA shape features. This allows a biophysical interpretation of the new parameters as a position-dependent preference towards specific DNA shape features. Using the new model, we explore the variation of DNA shape preferences in several transcription factors across various cancer cell lines and cellular conditions. The results reveal that there are DNA shape variations at FOXA1 (Forkhead Box Protein A1) binding sites in steroid-treated MCF7 cells. The new biophysical model is useful for elucidating the finer details of transcription factor–DNA interaction, as well as for predicting cancer mutation effects in the future. Keywords: transcription factors; DNA shape; protein–DNA interaction 1. Introduction Understanding how transcription factors (TFs) recognize their target DNA binding sites is an important task in the study of gene regulation.
    [Show full text]
  • Regulatory Motif Discovery Using Pwms and the Architecture of Eukaryotic Core Promoters
    Regulatory motif discovery using PWMs and the architecture of eukaryotic core promoters Holger Hartmann München 2012 Dissertation zur Erlangung des Doktorgrades der Fakultät für Chemie und Pharmazie der Ludwig-Maximilians-Universität München Regulatory motif discovery using PWMs and the architecture of eukaryotic core promoters Holger Hartmann aus Bad Friedrichshall, Deutschland 2012 Erklärung: Diese Dissertation wurde im Sinne von §7 der Promotionsordnung vom 28. November 2011 von Herrn Prof. Dr. Patrick Cramer betreut. Eidesstattliche Versicherung: Diese Dissertation wurde eigenständig und ohne unerlaubte Hilfe erarbeitet. München, am 29. März 2012 Holger Hartmann Dissertation eingereicht am 29. März 2012 1. Gutachter: Prof. Dr. Patrick Cramer 2. Gutachter: Dr. Johannes Söding Mündliche Prüfung am 8. Mai 2012 Acknowledgments During the past years of work on this thesis, I was influenced, assisted and motivated by many people. Having the opportunity to thank all of them is a pleasure not to be missed: First, I would like to thank Dr. Johannes Söding for being an amazing supervisor, offering me scientific freedom and constant support at the same time. I would like to thank Prof. Dr. Patrick Cramer for being my doctoral supervisor and the possibility to work in many fruitful collaborations. I am also very grateful to Prof. Dr. Ulrike Gaul, Dr. Dietmar Martin, Prof. Dr. Klaus Förstemann, and Prof. Dr. Roland Beckmann for offering their time as members of my dissertation committee. Furthermore, I would like to thank all students I had the pleasure to supervise (Claudia Gugenmus, Alexej Grjasnow, Lukas Utz, Anja Kiesel, Sebastian Luehr, and Mark Heron) for all input and help during my projects and the knowledge I gained from their analyses.
    [Show full text]