The genetic code as expressed through relationships between mRNA structure and protein function

MulStack: An ensemble learning prediction model of multilabel mRNA subcellular localization

Subcellular localization of mRNA is related to protein synthesis, cell polarity, cell movement and other biological regulation mechanisms. The distribution of mRNAs in subcellulars is similar to that of proteins, and most mRNAs are distributed in multiple subcellulars. Recently, some computational methods have been designed to predict the subcellular localization of mRNA. However, these methods only employed a sin-gle level of mRNA features and did not employ the position encoding of nucleotides in mRNA. In this paper, an ensemble learning prediction model is proposed, named MulStack, which is based on random forest and deep learning for multilabel mRNA subcellular localization. The proposed method employs two levels of mRNA features, including sequence-level and residue-level features, and position encoding is employed for the first time in the field of subcellular localization of mRNA. Random forest is employed to learn mRNA sequence-level feature, deep learning is employed to learn mRNA sequence-level feature and mRNA residue-level combined with position encoding. And the outputs of random forest and deep learning model will be weighted sum as the prediction probability. Compared with existing methods, the results show that MulStack is the best in the localization of the nucleus, cytosol and exosome. In addition, position weight matrices (PWMs) are extracted by convolutional neural networks (CNNs) that can be matched with known RNA binding protein motifs. Gene ontology (GO) enrichment analysis shows biological processes, molecular functions and cellular components of mRNA genes. The prediction web server of MulStack is freely accessible at http://bliulab.net/MulStack.

The effect of pseudoknot base pairing on cotranscriptional structural switching of the fluoride riboswitch

A central question in biology is how RNA sequence changes influence dynamic conformational changes during cotranscriptional folding. Here we investigated this question through the study of transcriptional fluoride riboswitches, non-coding RNAs that sense the fluoride anion through the coordinated folding and rearrangement of a pseudoknotted aptamer domain and a downstream intrinsic terminator expression platform. Using a combination of Escherichia coli RNA polymerase in vitro transcription and cellular gene expression assays, we characterized the function of mesophilic and thermophilic fluoride riboswitch variants. We showed that only variants containing the mesophilic pseudoknot function at 37°C. We next systematically varied the pseudoknot sequence and found that a single wobble base pair is critical for function. Characterizing thermophilic variants at 65°C through Thermus aquaticus RNA polymerase in vitro transcription showed the importance of this wobble pair for function even at elevated temperatures. Finally, we performed all-atom molecular dynamics simulations which supported the experimental findings, visualized the RNA structure switching process, and provided insight into the important role of magnesium ions. Together these studies provide deeper insights into the role of riboswitch sequence in influencing folding and function that will be important for understanding of RNA-based gene regulation and for synthetic biology applications.

The Effect of Pseudoknot Base Pairing on Cotranscriptional Structural Switching of the Fluoride Riboswitch

A central question in biology is how RNA sequence changes influence dynamic conformational changes during cotranscriptional folding. Here we investigated this question through the study of transcriptional fluoride riboswitches, non-coding RNAs that sense the fluoride anion through the coordinated folding and rearrangement of a pseudoknotted aptamer domain and a downstream intrinsic terminator expression platform. Using a combination of E. coli RNA polymerase in vitro transcription and cellular gene expression assays, we characterized the function of mesophilic and thermophilic fluoride riboswitch variants. We showed that only variants containing the mesophilic pseudoknot function at 37 °C. We next systematically varied the pseudoknot sequence and found that a single wobble base pair is critical for function. Characterizing thermophilic variants at 65 °C through Thermus aquaticus RNA polymerase in vitro transcription showed the importance of this wobble pair for function even at elevated temperatures. Finally, we performed all-atom molecular dynamics simulations which supported the experimental findings, visualized the RNA structure switching process, and provided insight into the important role of magnesium ions. Together these studies provide deeper insights into the role of riboswitch sequence in influencing folding and function that will be important for understanding of RNA-based gene regulation and for synthetic biology applications.

Mechanisms governing codon usage bias and the implications for protein expression in the chloroplast of Chlamydomonas reinhardtii

Chloroplasts possess a considerably reduced genome that is decoded via an almost minimal set of tRNAs. These features make an excellent platform for gaining insights into fundamental mechanisms that govern protein expression. Here, we present a comprehensive and revised perspective of the mechanisms that drive codon selection in the chloroplast of Chlamydomonas reinhardtii and the functional consequences for protein expression. In order to extract this information, we applied several codon usage descriptors to genes with different expression levels. We show that highly expressed genes strongly favor translationally optimal codons, while genes with lower functional importance are rather affected by directional mutational bias. We demonstrate that codon optimality can be deduced from codon-anticodon pairing affinity and, for a small number of amino acids (leucine, arginine, serine, and isoleucine), tRNA concentrations. Finally, we review, analyze, and expand on the impact of codon usage on protein yield, secondary structures of mRNA, translation initiation and termination, and amino acid composition of proteins, as well as cotranslational protein folding. The comprehensive analysis of codon choice provides crucial insights into heterologous gene expression in the chloroplast of C. reinhardtii, which may also be applicable to other chloroplast-containing organisms and bacteria.

A synonymous mutation in IGF-1 impacts the transcription and translation process of gene expression

Insulin-like growth factor 1 (IGF-1) is considered to be a crucial gene in the animal development of bone and body size. In this study, a unique synonymous mutation (c.258 A > G) of the IGF-1 gene was modified with an adenine base editor to observe the growth and developmental situation of mutant mice. Significant expression differences and molecular mechanisms among vectors with different alanine synonymous codons were explored. Although modification of a single synonymous codon rarely interferes with animal phenotypes, we observed that the expression and secretion of IGF-1 were different between 8-week-old homozygous (Ho) and wild-type (WT) mice. In addition, the IGF-1 with optimal codon combinations showed a higher expression content than other codon combination modes at both transcription and translation levels and performed proliferation promotion. The gene stability and translation initiation efficiency also changed significantly. Our findings illustrated that the synonymous mutation altered the IGF-1 gene expression in individual mice and suggested that the synonymous mutation affected the IGF-1 expression and biological function through the transcription and translation processes.

Ribosome as a Translocase and Helicase

During protein synthesis, ribosome moves along mRNA to decode one codon after the other. Ribosome translocation is induced by a universally conserved protein, elongation factor G (EF-G) in bacteria and elongation factor 2 (EF-2) in eukaryotes. EF-G-induced translocation results in unwinding of the intramolecular secondary structures of mRNA by three base pairs at a time that renders the translating ribosome a processive helicase. Professor Alexander Sergeevich Spirin has made numerous seminal contributions to understanding the molecular mechanism of translocation. Here, we review Spirin's insights into the ribosomal translocation and recent advances in the field that stemmed from Spirin's pioneering work. We also discuss key remaining challenges in studies of translocase and helicase activities of the ribosome.

IRC3 regulates mitochondrial translation in response to metabolic cues in Saccharomyces cerevisiae

Mitochondrial oxidative phosphorylation (OXPHOS) enzymes are made up of dual genetic origin. Mechanisms regulating the expression of nuclear-encoded OXPHOS subunits in response to metabolic cues (glucose vs. glycerol), is significantly understood while regulation of mitochondrially encoded OXPHOS subunits is poorly defined. Here, we show that IRC3 a DEAD/H box helicase, previously implicated in mitochondrial DNA maintenance, is central to integrating metabolic cues with mitochondrial translation. Irc3 associates with mitochondrial small ribosomal subunit in cells consistent with its role in regulating translation elongation based on Arg8^m reporter system. IRC3 deleted cells retained mitochondrial DNA despite growth defect on glycerol plates. Glucose grown Δirc3ρ⁺ and irc3 temperature-sensitive cells at 37⁰C have reduced translation rates from majority of mRNAs. In contrast, when galactose was the carbon source, reduction in mitochondrial translation was observed predominantly from Cox1 mRNA in Δirc3ρ⁺ but no defect was observed in irc3 temperature-sensitive cells, at 37⁰C. In support, of a model whereby IRC3 responds to metabolic cues to regulate mitochondrial translation, suppressors of Δirc3 isolated for restoration of growth on glycerol media restore mitochondrial protein synthesis differentially in presence of glucose vs. glycerol.

mLoc-mRNA: predicting multiple sub-cellular localization of mRNAs using random forest algorithm coupled with feature selection via elastic net

BACKGROUND

Localization of messenger RNAs (mRNAs) plays a crucial role in the growth and development of cells. Particularly, it plays a major role in regulating spatio-temporal gene expression. The in situ hybridization is a promising experimental technique used to determine the localization of mRNAs but it is costly and laborious. It is also a known fact that a single mRNA can be present in more than one location, whereas the existing computational tools are capable of predicting only a single location for such mRNAs. Thus, the development of high-end computational tool is required for reliable and timely prediction of multiple subcellular locations of mRNAs. Hence, we develop the present computational model to predict the multiple localizations of mRNAs.

RESULTS

The mRNA sequences from 9 different localizations were considered. Each sequence was first transformed to a numeric feature vector of size 5460, based on the k-mer features of sizes 1-6. Out of 5460 k-mer features, 1812 important features were selected by the Elastic Net statistical model. The Random Forest supervised learning algorithm was then employed for predicting the localizations with the selected features. Five-fold cross-validation accuracies of 70.87, 68.32, 68.36, 68.79, 96.46, 73.44, 70.94, 97.42 and 71.77% were obtained for the cytoplasm, cytosol, endoplasmic reticulum, exosome, mitochondrion, nucleus, pseudopodium, posterior and ribosome respectively. With an independent test set, accuracies of 65.33, 73.37, 75.86, 72.99, 94.26, 70.91, 65.53, 93.60 and 73.45% were obtained for the respective localizations. The developed approach also achieved higher accuracies than the existing localization prediction tools.

CONCLUSIONS

This study presents a novel computational tool for predicting the multiple localization of mRNAs. Based on the proposed approach, an online prediction server "mLoc-mRNA" is accessible at http://cabgrid.res.in:8080/mlocmrna/ . The developed approach is believed to supplement the existing tools and techniques for the localization prediction of mRNAs.

SHAPE Directed Discovery of New Functions in Large RNAs

RNA lies upstream of nearly all biology and functions as the central conduit of information exchange in all cells. RNA molecules encode information both in their primary sequences and in complex structures that form when an RNA folds back on itself. From the time of discovery of mRNA in the late 1950s until quite recently, we had only a rudimentary understanding of RNA structure across vast regions of most messenger and noncoding RNAs. This deficit is now rapidly being addressed, especially by selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) chemistry, mutational profiling (MaP), and closely related platform technologies that, collectively, create chemical microscopes for RNA. These technologies make it possible to interrogate RNA structure, quantitatively, at nucleotide resolution, and at large scales, for entire mRNAs, noncoding RNAs, and viral RNA genomes. By applying comprehensive structure probing to diverse problems, we and others are showing that control of biological function mediated by RNA structure is ubiquitous across prokaryotic and eukaryotic organisms.Work over the past decade using SHAPE-based analyses has clarified key principles. First, the method of RNA structure probing matters. SHAPE-MaP, with its direct and one-step readout that probes nearly every nucleotide by reaction at the 2'-hydroxyl, gives a more detailed and accurate readout than alternatives. Second, comprehensive chemical probing is essential. Focusing on fragments of large RNAs or using meta-gene or statistical analyses to compensate for sparse data sets misses critical features and often yields structure models with poor predictive power. Finally, every RNA has its own internal structural personality. There are myriad ways in which RNA structure modulates sequence accessibility, protein binding, translation, splice-site choice, phase separation, and other fundamental biological processes. In essentially every instance where we have applied rigorous and quantitative SHAPE technologies to study RNA structure-function interrelationships, new insights regarding biological regulatory mechanisms have emerged. RNA elements with more complex higher-order structures appear more likely to contain high-information-content clefts and pockets that bind small molecules, broadly informing a vigorous field of RNA-targeted drug discovery.The broad implications of this collective work are twofold. First, it is long past time to abandon depiction of large RNAs as simple noodle-like or gently flowing molecules. Instead, we need to emphasize that nearly all RNAs are punctuated with distinctive internal structures, a subset of which modulate function in profound ways. Second, structure probing should be an integral component of any effort that seeks to understand the functional nexuses and biological roles of large RNAs.

Making ends meet: New functions of mRNA secondary structure

The 5' cap and 3' poly(A) tail of mRNA are known to synergistically regulate mRNA translation and stability. Recent computational and experimental studies revealed that both protein-coding and non-coding RNAs will fold with extensive intramolecular secondary structure, which will result in close distances between the sequence ends. This proximity of the ends is a sequence-independent, universal property of most RNAs. Only low-complexity sequences without guanosines are without secondary structure and exhibit end-to-end distances expected for RNA random coils. The innate proximity of RNA ends might have important biological implications that remain unexplored. In particular, the inherent compactness of mRNA might regulate translation initiation by facilitating the formation of protein complexes that bridge mRNA 5' and 3' ends. Additionally, the proximity of mRNA ends might mediate coupling of 3' deadenylation to 5' end mRNA decay. This article is categorized under: RNA Structure and Dynamics > RNA Structure, Dynamics, and Chemistry RNA Structure and Dynamics > Influence of RNA Structure in Biological Systems Translation > Translation Regulation.

High-resolution modeling of the selection on local mRNA folding strength in coding sequences across the tree of life

BACKGROUND

mRNA can form local secondary structure within the protein-coding sequence, and the strength of this structure is thought to influence gene expression regulation. Previous studies suggest that secondary structure strength may be maintained under selection, but the details of this phenomenon are not well understood.

RESULTS

We perform a comprehensive study of the selection on local mRNA folding strengths considering variation between species across the tree of life. We show for the first time that local folding strength selection tends to follow a conserved characteristic profile in most phyla, with selection for weak folding at the two ends of the coding region and for strong folding elsewhere in the coding sequence, with an additional peak of selection for strong folding located downstream of the start codon. The strength of this pattern varies between species and organism groups, and we highlight contradicting cases. To better understand the underlying evolutionary process, we show that selection strengths in the different regions are strongly correlated, and report four factors which have a clear predictive effect on local mRNA folding selection within the coding sequence in different species.

CONCLUSIONS

The correlations observed between selection for local secondary structure strength in the different regions and with the four genomic and environmental factors suggest that they are shaped by the same evolutionary process throughout the coding sequence, and might be maintained under direct selection related to optimization of gene expression and specifically translation regulation.

The Role of the RNA-RNA Interactome in the Hepatitis C Virus Life Cycle

RNA virus genomes are multifunctional entities endowed with conserved structural elements that control translation, replication and encapsidation, among other processes. The preservation of these structural RNA elements constraints the genomic sequence variability. The hepatitis C virus (HCV) genome is a positive, single-stranded RNA molecule with numerous conserved structural elements that manage different steps during the infection cycle. Their function is ensured by the association of protein factors, but also by the establishment of complex, active, long-range RNA-RNA interaction networks-the so-called HCV RNA interactome. This review describes the RNA genome functions mediated via RNA-RNA contacts, and revisits some canonical ideas regarding the role of functional high-order structures during the HCV infective cycle. By outlining the roles of long-range RNA-RNA interactions from translation to virion budding, and the functional domains involved, this work provides an overview of the HCV genome as a dynamic device that manages the course of viral infection.

Effects of Refolding on Large-Scale RNA Structure

Understanding RNA-mediated functions requires a detailed characterization of the underlying RNA structure. In many cases, structure probing experiments are performed on RNA that has been "refolded" in some way, which may cause the conformation to differ from that of the native RNA. We used SHAPE-MaP (selective 2'-hydroxyl acylation analyzed by primer extension, read out by mutational profiling), to probe the structure of the Dengue virus (DENV) RNA genome after gentle extraction of the native RNA from intact virions (ex virion) and after heat denaturation and refolding. Comparison of mutiple SHAPE-informed structural features revealed that refolded RNA is more highly structured and samples fewer conformations than the ex virion RNA. Regions with similar structural features are generally those with low SHAPE reactivity and low Shannon entropy (lowSS regions), which correspond to elements with high levels of well-determined structure. This high-structure and low-entropy analysis framework, previously shown to make possible discovery of functional RNA structures, is thus now shown to allow de novo identification of structural elements in a refolded RNA that are likely to recapitulate RNA structures in the ex virion RNA state. Regions with less well-defined structures, which occurred more frequently in the more native-like ex virion RNA and may contain RNA switches, are challenging to recapitulate using refolded RNA.

Deletion of exons 17 and 18 in prestin's STAS domain results in loss of function

Cochlear outer hair cells (OHC) express the motor protein, prestin, which is required for sensitivity and frequency selectivity. Because our previous work showed that a calmodulin binding site (CBS) was located in prestin's C-terminal, specifically within the intrinsically disordered region, we sought to delete the IDR to study the functional significance of calcium-dependent, calmodulin binding on OHC function. Although the construct lacking the IDR (∆IDR prestin) demonstrated wildtype-like nonlinear capacitance (NLC) in HEK293T cells, the phenotype in ∆IDR prestin knockins (KI) was similar to that in prestin knockouts: thresholds were elevated, NLC was absent and OHCs were missing from basal regions of the cochlea. Although ∆IDR prestin mRNA was measured, no prestin protein was detected. At the mRNA level, both of prestin's exons 17 and 18 were entirely removed, rather than the smaller region encoding the IDR. Our hybrid exon that contained the targeted deletion (17-18 ∆IDR) failed to splice in vitro and prestin protein lacking exons 17 and 18 aggregated and failed to target the cell membrane. Hence, the absence of prestin protein in ∆IDR KI OHCs may be due to the unexpected splicing of the hybrid 17-18 ∆IDR exon followed by rapid degradation of nonfunctional prestin protein.

Conformational Dynamics of the RNA G-Quadruplex and its Effect on Translation Efficiency

During translation, intracellular mRNA folds co-transcriptionally and must refold following the passage of ribosome. The mRNAs can be entrapped in metastable structures during these folding events. In the present study, we evaluated the conformational dynamics of the kinetically favored, metastable, and hairpin-like structure, which disturbs the thermodynamically favored G-quadruplex structure, and its effect on co-transcriptional translation in prokaryotic cells. We found that nascent mRNA forms a metastable hairpin-like structure during co-transcriptional folding instead of the G-quadruplex structure. When the translation progressed co-transcriptionally before the metastable hairpin-like structure transition to the G-quadruplex, function of the G-quadruplex as a roadblock of the ribosome was sequestered. This suggested that kinetically formed RNA structures had a dominant effect on gene expression in prokaryotes. The results of this study indicate that it is critical to consider the conformational dynamics of RNA-folding to understand the contributions of the mRNA structures in controlling gene expression.

The evolution of RNA structural probing methods: From gels to next-generation sequencing

RNA molecules are important players in all domains of life and the study of the relationship between their multiple flexible states and the associated biological roles has increased in recent years. For several decades, chemical and enzymatic structural probing experiments have been used to determine RNA structure. During this time, there has been a steady improvement in probing reagents and experimental methods, and today the structural biologist community has a large range of tools at its disposal to probe the secondary structure of RNAs in vitro and in cells. Early experiments used radioactive labeling and polyacrylamide gel electrophoresis as read-out methods. This was superseded by capillary electrophoresis, and more recently by next-generation sequencing. Today, powerful structural probing methods can characterize RNA structure on a genome-wide scale. In this review, we will provide an overview of RNA structural probing methodologies from a historical and technical perspective. This article is categorized under: RNA Structure and Dynamics > RNA Structure, Dynamics, and Chemistry RNA Methods > RNA Analyses in vitro and In Silico RNA Methods > RNA Analyses in Cells.

mRNAs and lncRNAs intrinsically form secondary structures with short end-to-end distances

The 5' and 3' termini of RNA play important roles in many cellular processes. Using Förster resonance energy transfer (FRET), we show that mRNAs and lncRNAs have an intrinsic propensity to fold in the absence of proteins into structures in which the 5' end and 3' end are ≤7 nm apart irrespective of mRNA length. Computational estimates suggest that the inherent proximity of the ends is a universal property of most mRNA and lncRNA sequences. Only guanosine-depleted RNA sequences with low sequence complexity are unstructured and exhibit end-to-end distances expected for the random coil conformation of RNA. While the biological implications remain to be explored, short end-to-end distances could facilitate the binding of protein factors that regulate translation initiation by bridging mRNA 5' and 3' ends. Furthermore, our studies provide the basis for measuring, computing and manipulating end-to-end distances and secondary structure in RNA in research and biotechnology.

Chloroplast Translation: Structural and Functional Organization, Operational Control, and Regulation

Chloroplast translation is essential for cellular viability and plant development. Its positioning at the intersection of organellar RNA and protein metabolism makes it a unique point for the regulation of gene expression in response to internal and external cues. Recently obtained high-resolution structures of plastid ribosomes, the development of approaches allowing genome-wide analyses of chloroplast translation (i.e., ribosome profiling), and the discovery of RNA binding proteins involved in the control of translational activity have greatly increased our understanding of the chloroplast translation process and its regulation. In this review, we provide an overview of the current knowledge of the chloroplast translation machinery, its structure, organization, and function. In addition, we summarize the techniques that are currently available to study chloroplast translation and describe how translational activity is controlled and which cis-elements and trans-factors are involved. Finally, we discuss how translational control contributes to the regulation of chloroplast gene expression in response to developmental, environmental, and physiological cues. We also illustrate the commonalities and the differences between the chloroplast and bacterial translation machineries and the mechanisms of protein biosynthesis in these two prokaryotic systems.

Ensemble and single-molecule FRET studies of protein synthesis

Protein synthesis is a complex, multi-step process that involves large conformational changes of the ribosome and protein factors of translation. Over the last decade, Förster resonance energy transfer (FRET) has become instrumental for studying structural rearrangements of the translational apparatus. Here, we discuss the design of ensemble and single-molecule (sm) FRET assays of translation. We describe a number of experimental strategies that can be used to introduce fluorophores into the ribosome, tRNA, mRNA and protein factors of translation. Alternative approaches to tethering of translation components to the microscope slide in smFRET experiments are also reviewed. Finally, we discuss possible challenges in the interpretation of FRET data and ways to address these challenges.

Modeling RNA Secondary Structure with Sequence Comparison and Experimental Mapping Data

Secondary structure prediction is an important problem in RNA bioinformatics because knowledge of structure is critical to understanding the functions of RNA sequences. Significant improvements in prediction accuracy have recently been demonstrated though the incorporation of experimentally obtained structural information, for instance using selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) mapping. However, such mapping data is currently available only for a limited number of RNA sequences. In this article, we present a method for extending the benefit of experimental mapping data in secondary structure prediction to homologous sequences. Specifically, we propose a method for integrating experimental mapping data into a comparative sequence analysis algorithm for secondary structure prediction of multiple homologs, whereby the mapping data benefits not only the prediction for the specific sequence that was mapped but also other homologs. The proposed method is realized by modifying the TurboFold II algorithm for prediction of RNA secondary structures to utilize basepairing probabilities guided by SHAPE experimental data when such data are available. The SHAPE-mapping-guided basepairing probabilities are obtained using the RSample method. Results demonstrate that the SHAPE mapping data for a sequence improves structure prediction accuracy of other homologous sequences beyond the accuracy obtained by sequence comparison alone (TurboFold II). The updated version of TurboFold II is freely available as part of the RNAstructure software package.

Equalizing the information amounts of protein and mRNA by information theory

Based on the Shannon's information communication theory, information amount of the entire length of a polymeric macromolecule can be calculated in bits through adding the entropies of each building block. Proteins, DNA and RNA are such macromolecules. When only the building blocks' variation is considered as the source of entropy, there is seemingly lower information in case of the protein if this approach is applied directly on a protein of specific size and the coding sequence size of the mRNA corresponding to the particular length of the protein. This decrease in the information amount seems contradictory but this apparent conflict is resolved by considering the conformational variations in proteins as a new variable in the calculation and balancing the approximated entropy of the coding part of the mRNA and the protein. Probabilities can change therefore we also assigned hypothetical probabilities to the conformational states, which represent the uneven distribution as the time spent in one conformation, providing the probability of the presence in either or one of the possible conformations. Results that are obtained by using hypothetical probabilities are in line with the experimental values of variations in the conformational-state of protein populations. This equalization approach has further biological relevance that it compensates for the degeneracy in the codon usage during protein translation and it leads to the conclusion that the alphabet size for the protein is rather optimal for the proper protein functioning within the thermodynamic milieu of the cell. The findings were also discussed in relation to the codon bias and have implications in relation to the codon evolution concept. Eventually, this work brings the fields of protein structural studies and molecular protein translation processes together with a novel approach.

Detection of Ligand-Induced Conformational Changes in Oligonucleotides by Second-Harmonic Generation at a Supported Lipid Bilayer Interface

There is a high demand for characterizing oligonucleotide structural changes associated with binding interactions as well as identifying novel binders that modulate their structure and function. In this study, second-harmonic generation (SHG) was used to study RNA and DNA oligonucleotide conformational changes associated with ligand binding. For this purpose, we developed an avidin-based biotin capture surface based on a supported lipid bilayer membrane. The technique was applied to two well-characterized aptamers, both of which undergo conformational changes upon binding either a protein or a small molecule ligand. In both cases, SHG was able to resolve conformational changes in these oligonucleotides sensitively and specifically, in solution and in real time, using nanogram amounts of material. In addition, we developed a competition assay for the oligonucleotides between the specific ligands and known, nonspecific binders, and we demonstrated that intercalators and minor groove binders affect the conformation of the DNA and RNA oligonucleotides in different ways upon binding and subsequently block specific ligand binding in all cases. Our work demonstrates the broad potential of SHG for studying oligonucleotides and their conformational changes upon interaction with ligands. As SHG offers a powerful, high-throughput screening approach, our results here also open an important new avenue for identifying novel chemical probes or sequence-targeted drugs that disrupt or modulate DNA or RNA structure and function.

A small stem-loop structure of the Ebola virus trailer is essential for replication and interacts with heat-shock protein A8

Ebola virus (EBOV) is a single-stranded negative-sense RNA virus belonging to the Filoviridae family. The leader and trailer non-coding regions of the EBOV genome likely regulate its transcription, replication, and progeny genome packaging. We investigated the cis-acting RNA signals involved in RNA–RNA and RNA–protein interactions that regulate replication of eGFP-encoding EBOV minigenomic RNA and identified heat shock cognate protein family A (HSC70) member 8 (HSPA8) as an EBOV trailer-interacting host protein. Mutational analysis of the trailer HSPA8 binding motif revealed that this interaction is essential for EBOV minigenome replication. Selective 2′-hydroxyl acylation analyzed by primer extension analysis of the secondary structure of the EBOV minigenomic RNA indicates formation of a small stem-loop composed of the HSPA8 motif, a 3′ stem-loop (nucleotides 1868–1890) that is similar to a previously identified structure in the replicative intermediate (RI) RNA and a panhandle domain involving a trailer-to-leader interaction. Results of minigenome assays and an EBOV reverse genetic system rescue support a role for both the panhandle domain and HSPA8 motif 1 in virus replication.

Regulatory effects of cotranscriptional RNA structure formation and transitions

RNAs, which play significant roles in many fundamental biological processes of life, fold into sophisticated and precise structures. RNA folding is a dynamic and intricate process, which conformation transition of coding and noncoding RNAs form the primary elements of genetic regulation. The cellular environment contains various intrinsic and extrinsic factors that potentially affect RNA folding in vivo, and experimental and theoretical evidence increasingly indicates that the highly flexible features of the RNA structure are affected by these factors, which include the flanking sequence context, physiochemical conditions, cis RNA-RNA interactions, and RNA interactions with other molecules. Furthermore, distinct RNA structures have been identified that govern almost all steps of biological processes in cells, including transcriptional activation and termination, transcriptional mutagenesis, 5'-capping, splicing, 3'-polyadenylation, mRNA export and localization, and translation. Here, we briefly summarize the dynamic and complex features of RNA folding along with a wide variety of intrinsic and extrinsic factors that affect RNA folding. We then provide several examples to elaborate RNA structure-mediated regulation at the transcriptional and posttranscriptional levels. Finally, we illustrate the regulatory roles of RNA structure and discuss advances pertaining to RNA structure in plants. WIREs RNA 2016, 7:562-574. doi: 10.1002/wrna.1350 For further resources related to this article, please visit the WIREs website.

The Coding Region of the HCV Genome Contains a Network of Regulatory RNA Structures

RNA is a versatile macromolecule that accommodates functional information in primary sequence and secondary and tertiary structure. We use a combination of chemical probing, RNA structure modeling, comparative sequence analysis, and functional assays to examine the role of RNA structure in the hepatitis C virus (HCV) genome. We describe a set of conserved but functionally diverse structural RNA motifs that occur in multiple coding regions of the HCV genome, and we demonstrate that conformational changes in these motifs influence specific stages in the virus' life cycle. Our study shows that these types of structures can pervade a genome, where they play specific mechanistic and regulatory roles, constituting a "code within the code" for controlling biological processes.

LOOP IIId of the HCV IRES is essential for the structural rearrangement of the 40S-HCV IRES complex

As obligatory intracellular parasites, viruses rely on cellular machines to complete their life cycle, and most importantly they recruit the host ribosomes to translate their mRNA. The Hepatitis C viral mRNA initiates translation by directly binding the 40S ribosomal subunit in such a way that the initiation codon is correctly positioned in the P site of the ribosome. Such a property is likely to be central for many viruses, therefore the description of host-pathogen interaction at the molecular level is instrumental to provide new therapeutic targets. In this study, we monitored the 40S ribosomal subunit and the viral RNA structural rearrangement induced upon the formation of the binary complex. We further took advantage of an IRES viral mutant mRNA deficient for translation to identify the interactions necessary to promote translation. Using a combination of structure probing in solution and molecular modeling we establish a whole atom model which appears to be very similar to the one obtained recently by cryoEM. Our model brings new information on the complex, and most importantly reveals some structural rearrangement within the ribosome. This study suggests that the formation of a ‘kissing complex’ between the viral RNA and the 18S ribosomal RNA locks the 40S ribosomal subunit in a conformation proficient for translation.

Selective 2'-hydroxyl acylation analyzed by primer extension and mutational profiling (SHAPE-MaP) for direct, versatile and accurate RNA structure analysis

Selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) chemistries exploit small electrophilic reagents that react with 2'-hydroxyl groups to interrogate RNA structure at single-nucleotide resolution. Mutational profiling (MaP) identifies modified residues by using reverse transcriptase to misread a SHAPE-modified nucleotide and then counting the resulting mutations by massively parallel sequencing. The SHAPE-MaP approach measures the structure of large and transcriptome-wide systems as accurately as can be done for simple model RNAs. This protocol describes the experimental steps, implemented over 3 d, that are required to perform SHAPE probing and to construct multiplexed SHAPE-MaP libraries suitable for deep sequencing. Automated processing of MaP sequencing data is accomplished using two software packages. ShapeMapper converts raw sequencing files into mutational profiles, creates SHAPE reactivity plots and provides useful troubleshooting information. SuperFold uses these data to model RNA secondary structures, identify regions with well-defined structures and visualize probable and alternative helices, often in under 1 d. SHAPE-MaP can be used to make nucleotide-resolution biophysical measurements of individual RNA motifs, rare components of complex RNA ensembles and entire transcriptomes.

Toward a next-generation atlas of RNA secondary structure

RNA structure plays a crucial role in gene maturation, regulation and function. Determining the form and frequency of RNA folds is essential for a better understanding of how RNA exerts its functions. Low-throughput studies have focused on RNA primary sequences and expression levels, but with an emphasis on relatively small numbers of transcripts. However, with the recent advent of high-throughput technologies, it is realistic to begin analyzing RNA secondary structures on a genome-wide scale. Here, we review genome-wide RNA secondary structure profiles as well as advances in computational structure predictions. We further discuss the novel characteristics of RNA secondary structure across messenger RNAs. Probing RNA secondary structure by high-throughput sequencing will enable us to build atlases of RNA secondary structures, an important step in helping us to understand the versatility of RNA functions in diverse cellular processes.

The Landscape of long noncoding RNA classification

Advances in the depth and quality of transcriptome sequencing have revealed many new classes of long noncoding RNAs (lncRNAs). lncRNA classification has mushroomed to accommodate these new findings, even though the real dimensions and complexity of the noncoding transcriptome remain unknown. Although evidence of functionality of specific lncRNAs continues to accumulate, conflicting, confusing, and overlapping terminology has fostered ambiguity and lack of clarity in the field in general. The lack of fundamental conceptual unambiguous classification framework results in a number of challenges in the annotation and interpretation of noncoding transcriptome data. It also might undermine integration of the new genomic methods and datasets in an effort to unravel the function of lncRNA. Here, we review existing lncRNA classifications, nomenclature, and terminology. Then, we describe the conceptual guidelines that have emerged for their classification and functional annotation based on expanding and more comprehensive use of large systems biology-based datasets.

Decoding mechanisms by which silent codon changes influence protein biogenesis and function

SCOPE

Synonymous codon usage has been a focus of investigation since the discovery of the genetic code and its redundancy. The occurrences of synonymous codons vary between species and within genes of the same genome, known as codon usage bias. Today, bioinformatics and experimental data allow us to compose a global view of the mechanisms by which the redundancy of the genetic code contributes to the complexity of biological systems from affecting survival in prokaryotes, to fine tuning the structure and function of proteins in higher eukaryotes. Studies analyzing the consequences of synonymous codon changes in different organisms have revealed that they impact nucleic acid stability, protein levels, structure and function without altering amino acid sequence. As such, synonymous mutations inevitably contribute to the pathogenesis of complex human diseases. Yet, fundamental questions remain unresolved regarding the impact of silent mutations in human disorders. In the present review we describe developments in this area concentrating on mechanisms by which synonymous mutations may affect protein function and human health.

PURPOSE

This synopsis illustrates the significance of synonymous mutations in disease pathogenesis. We review the different steps of gene expression affected by silent mutations, and assess the benefits and possible harmful effects of codon optimization applied in the development of therapeutic biologics.

PHYSIOLOGICAL AND MEDICAL RELEVANCE

Understanding mechanisms by which synonymous mutations contribute to complex diseases such as cancer, neurodegeneration and genetic disorders, including the limitations of codon-optimized biologics, provides insight concerning interpretation of silent variants and future molecular therapies.

IRE mRNA riboregulators use metabolic iron (Fe(2+)) to control mRNA activity and iron chemistry in animals

A family of noncoding RNAs bind Fe(2+) to increase protein synthesis. The structures occur in messenger RNAs encoding animal proteins for iron metabolism. Each mRNA regulatory sequence, ∼30 ribonucleotides long, is called an IRE (Iron Responsive Element), and folds into a bent, A-RNA helix with a terminal loop. Riboregulatory RNAs, like t-RNAs, r-RNAs micro-RNAs, etc. contrast with DNA, since single-stranded RNA can fold into a variety of complex, three-dimensional structures. IRE-RNAs bind two types of proteins: (1) IRPs which are protein repressors, sequence-related to mitochondrial aconitases. (2) eIF-4F, which bind ribosomes and enhances general protein biosynthesis. The competition between IRP and eIF-4F binding to IRE-RNA is controlled by Fe(2+)-induced changes in the IRE-RNA conformation. Mn(2+), which also binds to IRE-RNA in solution, is a convenient experimental proxy for air-sensitive Fe(2+) studies of in vitro protein biosynthesis and protein binding. However, only Fe(2+) has physiological effects on protein biosynthesis directed by IRE-mRNAs. The structures of the IRE-RNA riboregulators is known indirectly from effects of base substitutions on function, from solution NMR of the free RNA, and of X-ray crystallography of the IRE-RNA-IRP repressor complex. However, the inability to date, to crystallize the free IRE-RNA, and the dissociation of the IRE-RNA-IRP complex when metal binds, have hampered direct identification and characterization of the RNA-metal binding sites. The high conservation of the primary sequence in IRE-mRNA control elements has facilitated their identification and analysis of metal-assisted riboregulator function. Expansion of RNA search analyses beyond primary will likely reveal other, metal-dependent families of mRNA riboregulators.

RNA structural analysis by evolving SHAPE chemistry

RNA is central to the flow of biological information. From transcription to splicing, RNA localization, translation, and decay, RNA is intimately involved in regulating every step of the gene expression program, and is thus essential for health and understanding disease. RNA has the unique ability to base-pair with itself and other nucleic acids to form complex structures. Hence the information content in RNA is not simply its linear sequence of bases, but is also encoded in complex folding of RNA molecules. A general chemical functionality that all RNAs have is a 2'-hydroxyl group in the ribose ring, and the reactivity of the 2'-hydroxyl in RNA is gated by local nucleotide flexibility. In other words, the 2'-hydroxyl is reactive at single-stranded and conformationally flexible positions but is unreactive at nucleotides constrained by base-pairing. Recent efforts have been focused on developing reagents that modify RNA as a function of RNA 2' hydroxyl group reactivity. Such RNA structure probing techniques can be read out by primer extension in experiments termed RNA SHAPE (selective 2'- hydroxyl acylation and primer extension). Herein, we describe the efforts devoted to the design and utilization of SHAPE probes for characterizing RNA structure. We also describe current technological advances that are being applied to utilize SHAPE chemistry with deep sequencing to probe many RNAs in parallel. The merging of chemistry with genomics is sure to open the door to genome-wide exploration of RNA structure and function.

RNA secondary structure modeling at consistent high accuracy using differential SHAPE

RNA secondary structure modeling is a challenging problem, and recent successes have raised the standards for accuracy, consistency, and tractability. Large increases in accuracy have been achieved by including data on reactivity toward chemical probes: Incorporation of 1M7 SHAPE reactivity data into an mfold-class algorithm results in median accuracies for base pair prediction that exceed 90%. However, a few RNA structures are modeled with significantly lower accuracy. Here, we show that incorporating differential reactivities from the NMIA and 1M6 reagents--which detect noncanonical and tertiary interactions--into prediction algorithms results in highly accurate secondary structure models for RNAs that were previously shown to be difficult to model. For these RNAs, 93% of accepted canonical base pairs were recovered in SHAPE-directed models. Discrepancies between accepted and modeled structures were small and appear to reflect genuine structural differences. Three-reagent SHAPE-directed modeling scales concisely to structurally complex RNAs to resolve the in-solution secondary structure analysis problem for many classes of RNA.

Insights into RNA structure and function from genome-wide studies

A comprehensive understanding of RNA structure will provide fundamental insights into the cellular function of both coding and non-coding RNAs. Although many RNA structures have been analysed by traditional biophysical and biochemical methods, the low-throughput nature of these approaches has prevented investigation of the vast majority of cellular transcripts. Triggered by advances in sequencing technology, genome-wide approaches for probing the transcriptome are beginning to reveal how RNA structure affects each step of protein expression and RNA stability. In this Review, we discuss the emerging relationships between RNA structure and the regulation of gene expression.

Deciphering the rules by which dynamics of mRNA secondary structure affect translation efficiency in Saccharomyces cerevisiae

Messenger RNA (mRNA) secondary structure decreases the elongation rate, as ribosomes must unwind every structure they encounter during translation. Therefore, the strength of mRNA secondary structure is assumed to be reduced in highly translated mRNAs. However, previous studies in vitro reported a positive correlation between mRNA folding strength and protein abundance. The counterintuitive finding suggests that mRNA secondary structure affects translation efficiency in an undetermined manner. Here, we analyzed the folding behavior of mRNA during translation and its effect on translation efficiency. We simulated translation process based on a novel computational model, taking into account the interactions among ribosomes, codon usage and mRNA secondary structures. We showed that mRNA secondary structure shortens ribosomal distance through the dynamics of folding strength. Notably, when adjacent ribosomes are close, mRNA secondary structures between them disappear, and codon usage determines the elongation rate. More importantly, our results showed that the combined effect of mRNA secondary structure and codon usage in highly translated mRNAs causes a short ribosomal distance in structural regions, which in turn eliminates the structures during translation, leading to a high elongation rate. Together, these findings reveal how the dynamics of mRNA secondary structure coupling with codon usage affect translation efficiency.

Regulation of mRNA transport, localization and translation in the nervous system of mammals (Review)

Post-transcriptional control of mRNA trafficking and metabolism plays a critical role in the actualization and fine tuning of the genetic program of cells, both in development and in differentiated tissues. Cis-acting signals, responsible for post-transcriptional regulation, reside in the RNA message itself, usually in untranslated regions, 5′ or 3′ to the coding sequence, and are recognized by trans-acting factors: RNA-binding proteins (RBPs) and/or non-coding RNAs (ncRNAs). ncRNAs bind short mRNA sequences usually present in the 3′-untranslated (3′-UTR) region of their target messages. RBPs recognize specific nucleotide sequences and/or secondary/tertiary structures. Most RBPs assemble on mRNA at the moment of transcription and shepherd it to its destination, somehow determining its final fate. Regulation of mRNA localization and metabolism has a particularly important role in the nervous system where local translation of pre-localized mRNAs has been implicated in developing axon and dendrite pathfinding, and in synapse formation. Moreover, activity-dependent mRNA trafficking and local translation may underlie long-lasting changes in synaptic efficacy, responsible for learning and memory. This review focuses on the role of RBPs in neuronal development and plasticity, as well as possible connections between ncRNAs and RBPs.

Transcriptome analyses of a Chinese hazelnut species Corylus mandshurica

BACKGROUND

Corylus was renowned for its production of hazelnut and taxol. To understand the local adaptation of Chinese species and speed up breeding efforts in China, we analyzed the leaf transcriptome of Corylus mandshurica, which had a high tolerance to fungal infections and cold.

RESULTS

A total of 12,255,030 clean pair-end reads were generated and then assembled into 37,846 Expressed Sequence Tag (EST) sequences. During functional annotation, 26,565 ESTs were annotated with Gene Ontology (GO) terms using Blast2go and 11,056 ESTs were grouped into the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways using KEGG Automatic Annotation Server (KAAS). We identified 45 ESTs that were homologous to enzymes and transcription factors responsible for taxol synthesis. The most differentiated orthologs between C. mandshurica and a European congener, C. avellana, were enriched in stress tolerance to fungal resistance and cold.

CONCLUSIONS

In this study, we detected a set of genes related to taxol synthesis in a taxol-producing angiosperm species for the first time and found a close relationship between most differentiated genes and different adaptations to fungal infection and cold in C. mandshurica and C. avellana. These findings provided tools to improve our understanding of local adaptation, genetic breeding and taxol production in hazelnut.

The genetic code as expressed through relationships between mRNA structure and protein function

Structured RNA elements within messenger RNA often direct or modulate the cellular production of active proteins. As reviewed here, RNA structures have been discovered that govern nearly every step in protein production: mRNA production and stability; translation initiation, elongation, and termination; protein folding; and cellular localization. Regulatory RNA elements are common within RNAs from every domain of life. This growing body of RNA-mediated mechanisms continues to reveal new ways in which mRNA structure regulates translation. We integrate examples from several different classes of RNA structure-mediated regulation to present a global perspective that suggests that the secondary and tertiary structure of RNA ultimately constitutes an additional level of the genetic code that both guides and regulates protein biosynthesis.