RESUMEN
K-mer based comparisons have emerged as powerful complements to BLAST-like alignment algorithms, particularly when the sequences being compared lack direct evolutionary relationships. In this chapter, we describe methods to compare k-mer content between groups of long noncoding RNAs (lncRNAs), to identify communities of lncRNAs with related k-mer contents, to identify the enrichment of protein-binding motifs in lncRNAs, and to scan for domains of related k-mer contents in lncRNAs. Our step-by-step instructions are complemented by Python code deposited in Github. Though our chapter focuses on lncRNAs, the methods we describe could be applied to any set of nucleic acid sequences.
Asunto(s)
Biología Computacional/métodos , ARN Largo no Codificante/clasificación , ARN Largo no Codificante/genética , Algoritmos , Análisis por Conglomerados , Motivos de Nucleótidos/genética , Unión ProteicaRESUMEN
The marsupial inactive X chromosome expresses a long noncoding RNA (lncRNA) called Rsx that has been proposed to be the functional analog of eutherian Xist Despite the possibility that Xist and Rsx encode related functions, the two lncRNAs harbor no linear sequence similarity. However, both lncRNAs harbor domains of tandemly repeated sequence. In Xist, these repeat domains are known to be critical for function. Using k-mer based comparison, we show that the repeat domains of Xist and Rsx unexpectedly partition into two major clusters that each harbor substantial levels of nonlinear sequence similarity. Xist Repeats B, C, and D were most similar to each other and to Rsx Repeat 1, whereas Xist Repeats A and E were most similar to each other and to Rsx Repeats 2, 3, and 4. Similarities at the level of k-mers corresponded to domain-specific enrichment of protein-binding motifs. Within individual domains, protein-binding motifs were often enriched to extreme levels. Our data support the hypothesis that Xist and Rsx encode similar functions through different spatial arrangements of functionally analogous protein-binding domains. We propose that the two clusters of repeat domains in Xist and Rsx function in part to cooperatively recruit PRC1 and PRC2 to chromatin. The physical manner in which these domains engage with protein cofactors may be just as critical to the function of the domains as the protein cofactors themselves. The general approaches we outline in this report should prove useful in the study of any set of RNAs.
Asunto(s)
Marsupiales/genética , ARN Largo no Codificante/química , ARN Largo no Codificante/genética , Animales , Análisis por Conglomerados , Humanos , Marsupiales/metabolismo , Proteínas del Grupo Polycomb/metabolismo , Dominios Proteicos , Homología de Secuencia de Ácido Nucleico , Secuencias Repetidas en Tándem , Inactivación del Cromosoma XRESUMEN
The functions of most long non-coding RNAs (lncRNAs) are unknown. In contrast to proteins, lncRNAs with similar functions often lack linear sequence homology; thus, the identification of function in one lncRNA rarely informs the identification of function in others. We developed a sequence comparison method to deconstruct linear sequence relationships in lncRNAs and evaluate similarity based on the abundance of short motifs called k-mers. We found that lncRNAs of related function often had similar k-mer profiles despite lacking linear homology, and that k-mer profiles correlated with protein binding to lncRNAs and with their subcellular localization. Using a novel assay to quantify Xist-like regulatory potential, we directly demonstrated that evolutionarily unrelated lncRNAs can encode similar function through different spatial arrangements of related sequence motifs. K-mer-based classification is a powerful approach to detect recurrent relationships between sequence and function in lncRNAs.
Asunto(s)
Motivos de Nucleótidos , ARN Largo no Codificante/clasificación , ARN Largo no Codificante/genética , Análisis de Secuencia de ARN/métodos , Algoritmos , Animales , Secuencia de Bases , Análisis por Conglomerados , Secuencia Conservada , Bases de Datos Genéticas , Células Hep G2 , Humanos , Células K562 , Ratones , Anotación de Secuencia Molecular , Conformación de Ácido Nucleico , Motivos de Nucleótidos/genética , Canales de Potasio con Entrada de Voltaje/genética , ARN Largo no Codificante/química , Alineación de SecuenciaRESUMEN
In this work, we have examined contributions to the thermodynamics of calmodulin (CaM) binding from the intrinsic propensity for target peptides to adopt an α-helical conformation. CaM target sequences are thought to commonly reside in disordered regions within proteins. Using the ability of TFE to induce α-helical structure as a proxy, the six peptides studied range from having almost no propensity to adopt α-helical structure through to a very high propensity. This despite all six peptides having similar CaM-binding affinities. Our data indicate there is some correlation between the deduced propensities and the thermodynamics of CaM binding. This finding implies that molecular recognition features, such as CaM target sequences, may possess a broad range of propensities to adopt local structure. Given that these peptides bind to CaM with similar affinities, the data suggest that having a higher propensity to adopt α-helical structure does not necessarily result in tighter binding, and that the mechanism of CaM binding is very dependent on the nature of the substrate sequence.