Nanotechnology and synthetic biology currently constitute one of the most innovative, interdisciplinary fields of research, poised to radically transform society in the 21st century. This paper concerns the synthetic design of ribonucleic acid molecules, using our recent algorithm, , which can determine all RNA sequences whose minimum free energy secondary structure is a user-specified target structure. Using , we design ten cis-cleaving hammerhead ribozymes, all of which are shown to be functional by a cleavage assay. We additionally use to design a functional cis-cleaving hammerhead as a modular unit of a synthetic larger RNA. Analysis of kinetics on this small set of hammerheads suggests that cleavage rate of computationally designed ribozymes may be correlated with positional entropy, ensemble defect, structural flexibility/rigidity and related measures. Artificial ribozymes have been designed in the past either manually or by SELEX (Systematic Evolution of Ligands by Exponential Enrichment); however, this appears to be the first purely computational design and experimental validation of novel functional ribozymes. is available at http://bioinformatics.bc.edu/clotelab/RNAiFold/.

BackgroundAutogenous cis-regulators of ribosomal protein synthesis play a critical role in maintaining the stoichiometry of ribosome components. Structured portions within an mRNA transcript typically interact with specific ribosomal proteins to prevent expression of the entire operon, thus balancing levels of ribosomal proteins across transcriptional units. Three distinct RNA structures from different bacterial phyla have demonstrated interactions with S15 to regulate gene expression; however, these RNAs are distributed across a small fraction of bacterial diversity.ResultsWe used comparative genomics in combination with analysis of existing transcriptomic data to identify three novel putative RNA structures associated with the S15 coding region in microbial genomes. These structures are completely distinct from those previously published and encompass potential regulatory regions including ribosome-binding sites. To validate the biological relevance of our findings, we demonstrate that an example of the Alphaproteobacterial RNA from Rhizobium radiobacter specifically interacts with S15 in vitro, and allows in vivo regulation of gene expression in an E. coli reporter system. In addition, structural probing and nuclease protection assays confirm the predicted secondary structure and indicate nucleotides required for protein interaction.ConclusionsThis work illustrates the importance of integrating comparative genomic and transcriptomic approaches during de novo ncRNA identification and reveals a diversity of distinct natural RNA regulators that support analogous biological functions. Furthermore, this work indicates that many additional uncharacterized RNA regulators likely exist within bacterial genomes and that the plasticity of RNA structure allows unique, and likely independently derived, solutions to the same biological problem.Electronic supplementary materialThe online version of this article (doi:10.1186/1471-2164-15-657) contains supplementary material, which is available to authorized users.

RNA-protein interactions are critical in many biological processes, yet how such interactions affect the evolution of both partners is still unknown. RNA and protein structures are impacted very differently by mechanisms of genomic change. While most protein families are identifiable at the nucleotide level across large phylogenetic distances, RNA families display far less nucleotide similarity and are often only shared by closely related bacterial species. Ribosomal protein S15 has two RNA binding functions. First, it is a ribosomal protein responsible for organizing the rRNA during ribosome assembly. Second, in many bacterial species S15 also interacts with a structured portion of its own transcript to negatively regulate gene expression. While the first interaction is conserved in most bacteria, the second is not. Four distinct mRNA structures interact with S15 to enable regulation, each of which appears to be independently derived in different groups of bacteria. With the goal of understanding how protein-binding specificity may influence the evolution of such RNA regulatory structures, we examine whether examples of these mRNA structures are able to interact with, and regulate in response to, S15 homologs from organisms containing distinct mRNA structures. We find that despite their shared RNA binding function in the rRNA, S15 homologs have distinct RNA recognition profiles. We present a model to explain the specificity patterns observed, and support this model by with further mutagenesis. After analyzing the patterns of conservation for the S15 protein coding sequences, we also identified amino acid changes that alter the binding specificity of an S15 homolog. In this work we demonstrate that homologous RNA-binding proteins have different specificity profiles, and minor changes to amino acid sequences, or to RNA structural motifs, can have large impacts on RNA-protein recognition.

There are several natural examples of distinct RNA structures that interact with the same ligand to regulate the expression of homologous genes in different organisms. One essential question regarding this phenomenon is whether such RNA regulators are the result of convergent or divergent evolution. Are the RNAs derived from some common ancestor and diverged to the point where we cannot identify the similarity, or have multiple solutions to the same biological problem arisen independently? A key variable in assessing these alternatives is how frequently such regulators arise within sequence space. Ribosomal protein S15 is autogenously regulated via an RNA regulator in many bacterial species; four apparently distinct regulators have been functionally validated in different bacterial phyla. Here, we explore how frequently such regulators arise within a partially randomized sequence population. We find many RNAs that interact specifically with ribosomal protein S15 from Geobacillus kaustophilus with biologically relevant dissociation constants. Furthermore, of the six sequences we characterize, four show regulatory activity in an Escherichia coli reporter assay. Subsequent footprinting and mutagenesis analysis indicates that protein binding proximal to regulatory features such as the Shine–Dalgarno sequence is sufficient to enable regulation, suggesting that regulation in response to S15 is relatively easily acquired.

The Burkholderia cepacia complex (Bcc) is a family of closely related bacterial pathogens that are the causative agent of deadly human infections. Virulence in Bcc species has been shown to be controlled by the CepI/CepR quorum sensing (QS) system, which is mediated by an N-acyl L-homoserine lactone (AHL) signal (C8-AHL) and its cognate LuxR-type receptor (CepR). Chemical strategies to block QS in Bcc members would represent an approach to intercept this bacterial communication process and further delineate its role in infection. In the current study, we sought to identify non-native AHLs capable of agonizing or antagonizing CepR, and thereby QS, in a Bcc member. We screened a library of AHL analogs in cell-based reporters for CepR, and identified numerous highly potent CepR agonists and antagonists. These compounds remain active in a Bcc member, B. multivorans, with one agonist 250-fold more potent than the native ligand C8-AHL, and can affect QS-controlled motility. Further, the CepR antagonists prolong C. elegans survival in an infection model. These AHL analogs are the first reported non-native molecules that both directly modulate CepR and impact QS-controlled phenotypes in a Bcc member, and represent valuable chemical tools to assess the role of QS in Bcc infections.