SUMMARY
Aptamers, short RNA or DNA molecules that bind distinct targets with high affinity and specificity, can be identified using High Throughput Systematic Evolution of Ligands by Exponential Enrichment (HT-SELEX). But scalable analytic tools for understanding sequence-function relationships from diverse HT-SELEX data are not available. Here, we present AptaTRACE, a computational approach that leverages the experimental design of the HT-SELEX protocol, RNA secondary structure, and the potential presence of many secondary motifs to identify sequence-structure motifs that show a signature of selection. We apply AptaTRACE to identify nine motifs in C-C chemokine receptor type 7 targeted by aptamers in an in vitro cell-SELEX experiment. We experimentally validate two aptamers whose binding required both sequence and structural features. AptaTRACE can identify low-abundance motifs, and we show through simulations that because of this it could lower HT-SELEX cost and time by reducing the number of selection cycles required. AptaTRACE is available for download at www.ncbi.nlm.nih.gov/CBBresearch/Przytycka/index.cgi#aptatools.