Computational Pipeline for Reference-Free Comparative Analysis of RNA 3D Structures Applied to SARS-CoV-2 UTR Models.

RNA is a unique biomolecule that is involved in a variety of fundamental biological functions, all of which depend solely on its structure and dynamics. Since the experimental determination of crystal RNA structures is laborious, computational 3D structure prediction methods are experiencing an ongoing and thriving development. Such methods can lead to many models; thus, it is necessary to build comparisons and extract common structural motifs for further medical or biological studies. Here, we introduce a computational pipeline dedicated to reference-free high-throughput comparative analysis of 3D RNA structures. We show its application in the RNA-Puzzles challenge, in which five participating groups attempted to predict the three-dimensional structures of 5'- and 3'-untranslated regions (UTRs) of the SARS-CoV-2 genome. We report the results of this puzzle and discuss the structural motifs obtained from the analysis. All simulated models and tools incorporated into the pipeline are open to scientific and academic use.

MAPIYA contact map server for identification and visualization of molecular interactions in proteins and biological complexes.

Correct identification and effective visualization of interactions in biomolecular structures facilitate understanding of their functions and molecular design. In response to the practical needs of structure-based analysis, we have created a Mapiya web server. The Mapiya integrates four main functionalities: (i) generation of contact maps - intramolecular and intermolecular-for proteins, nucleic acids, and their complexes; (ii) characterization of the interactions physicochemical nature, (iii) interactive visualization of biomolecular conformations with automatic zoom on selected contacts using Molstar and (iv) additional sequence- and structure-based analyses performed with third-party software and in-house algorithms combined into an easy-to-use interface. Thus, Mapiya offers a highly customized analysis of the molecular interactions' in various biological systems. The web server is available at: http://mapiya.lcbio.pl/.

Molecular insights into RNA recognition and gene regulation by the TRIM-NHL protein Mei-P26.

The TRIM-NHL protein Meiotic P26 (Mei-P26) acts as a regulator of cell fate in Drosophila Its activity is critical for ovarian germline stem cell maintenance, differentiation of oocytes, and spermatogenesis. Mei-P26 functions as a post-transcriptional regulator of gene expression; however, the molecular details of how its NHL domain selectively recognizes and regulates its mRNA targets have remained elusive. Here, we present the crystal structure of the Mei-P26 NHL domain at 1.6 Å resolution and identify key amino acids that confer substrate specificity and distinguish Mei-P26 from closely related TRIM-NHL proteins. Furthermore, we identify mRNA targets of Mei-P26 in cultured Drosophila cells and show that Mei-P26 can act as either a repressor or activator of gene expression on different RNA targets. Our work reveals the molecular basis of RNA recognition by Mei-P26 and the fundamental functional differences between otherwise very similar TRIM-NHL proteins.

NCodR: A multi-class support vector machine classification to distinguish non-coding RNAs in Viridiplantae.

Non-coding RNAs (ncRNAs) are major players in the regulation of gene expression. This study analyses seven classes of ncRNAs in plants using sequence and secondary structure-based RNA folding measures. We observe distinct regions in the distribution of AU content along with overlapping regions for different ncRNA classes. Additionally, we find similar averages for minimum folding energy index across various ncRNAs classes except for pre-miRNAs and lncRNAs. Various RNA folding measures show similar trends among the different ncRNA classes except for pre-miRNAs and lncRNAs. We observe different k-mer repeat signatures of length three among various ncRNA classes. However, in pre-miRs and lncRNAs, a diffuse pattern of k-mers is observed. Using these attributes, we train eight different classifiers to discriminate various ncRNA classes in plants. Support vector machines employing radial basis function show the highest accuracy (average F1 of ~96%) in discriminating ncRNAs, and the classifier is implemented as a web server, NCodR.

Constrained peptides mimic a viral suppressor of RNA silencing.

The design of high-affinity, RNA-binding ligands has proven very challenging. This is due to the unique structural properties of RNA, often characterized by polar surfaces and high flexibility. In addition, the frequent lack of well-defined binding pockets complicates the development of small molecule binders. This has triggered the search for alternative scaffolds of intermediate size. Among these, peptide-derived molecules represent appealing entities as they can mimic structural features also present in RNA-binding proteins. However, the application of peptidic RNA-targeting ligands is hampered by a lack of design principles and their inherently low bio-stability. Here, the structure-based design of constrained α-helical peptides derived from the viral suppressor of RNA silencing, TAV2b, is described. We observe that the introduction of two inter-side chain crosslinks provides peptides with increased α-helicity and protease stability. One of these modified peptides (B3) shows high affinity for double-stranded RNA structures including a palindromic siRNA as well as microRNA-21 and its precursor pre-miR-21. Notably, B3 binding to pre-miR-21 inhibits Dicer processing in a biochemical assay. As a further characteristic this peptide also exhibits cellular entry. Our findings show that constrained peptides can efficiently mimic RNA-binding proteins rendering them potentially useful for the design of bioactive RNA-targeting ligands.

RNA and DNA G-quadruplexes bind to human dicer and inhibit its activity.

Guanine (G)-rich single-stranded nucleic acids can adopt G-quadruplex structures. Accumulating evidence indicates that G-quadruplexes serve important regulatory roles in fundamental biological processes such as DNA replication, transcription, and translation, while aberrant G-quadruplex formation is linked to genome instability and cancer. Understanding the biological functions played by G-quadruplexes requires detailed knowledge of their protein interactome. Here, we report that both RNA and DNA G-quadruplexes are bound by human Dicer in vitro. Using in vitro binding assays, mutation studies, and computational modeling we demonstrate that G-quadruplexes can interact with the Platform-PAZ-Connector helix cassette of Dicer, the region responsible for anchoring microRNA precursors (pre-miRNAs). Consequently, we show that G-quadruplexes efficiently and stably inhibit the cleavage of pre-miRNA by Dicer. Our data highlight the potential of human Dicer for binding of G-quadruplexes and allow us to propose a G-quadruplex-driven sequestration mechanism of Dicer regulation.

Genome-wide mapping of SARS-CoV-2 RNA structures identifies therapeutically-relevant elements.

SARS-CoV-2 is a betacoronavirus with a linear single-stranded, positive-sense RNA genome, whose outbreak caused the ongoing COVID-19 pandemic. The ability of coronaviruses to rapidly evolve, adapt, and cross species barriers makes the development of effective and durable therapeutic strategies a challenging and urgent need. As for other RNA viruses, genomic RNA structures are expected to play crucial roles in several steps of the coronavirus replication cycle. Despite this, only a handful of functionally-conserved coronavirus structural RNA elements have been identified to date. Here, we performed RNA structure probing to obtain single-base resolution secondary structure maps of the full SARS-CoV-2 coronavirus genome both in vitro and in living infected cells. Probing data recapitulate the previously described coronavirus RNA elements (5' UTR and s2m), and reveal new structures. Of these, ∼10.2% show significant covariation among SARS-CoV-2 and other coronaviruses, hinting at their functionally-conserved role. Secondary structure-restrained 3D modeling of these segments further allowed for the identification of putative druggable pockets. In addition, we identify a set of single-stranded segments in vivo, showing high sequence conservation, suitable for the development of antisense oligonucleotide therapeutics. Collectively, our work lays the foundation for the development of innovative RNA-targeted therapeutic strategies to fight SARS-related infections.

Modeling of Three-Dimensional RNA Structures Using SimRNA.

The molecules of the ribonucleic acid (RNA) perform a variety of vital roles in all living cells. Their biological function depends on their structure and dynamics, both of which are difficult to experimentally determine but can be theoretically inferred based on the RNA sequence. SimRNA is one of the computational methods for molecular simulations of RNA 3D structure formation. The method is based on a simplified (coarse-grained) representation of nucleotide chains, a statistically derived model of interactions (statistical potential), and the Monte Carlo method as a conformational sampling scheme.The current version of SimRNA (3.22) is able to predict basic topologies of RNA molecules with sizes up to about 50-70 nucleotides, based on their sequences only, and larger molecules if supplied with appropriate distance restraints. The user can specify various types of restraints, including secondary structure, pairwise atom-atom distances, and positions of atoms. SimRNA can be also used for studying systems composed of several chains of RNA. SimRNA is a folding simulations method, thus it allows for examining folding pathways, getting an approximate view of the energy landscapes.

A structure-based model for the prediction of protein-RNA binding affinity.

Protein-RNA recognition is highly affinity-driven and regulates a wide array of cellular functions. In this study, we have curated a binding affinity data set of 40 protein-RNA complexes, for which at least one unbound partner is available in the docking benchmark. The data set covers a wide affinity range of eight orders of magnitude as well as four different structural classes. On average, we find the complexes with single-stranded RNA have the highest affinity, whereas the complexes with the duplex RNA have the lowest. Nevertheless, free energy gain upon binding is the highest for the complexes with ribosomal proteins and the lowest for the complexes with tRNA with an average of -5.7 cal/mol/Å2 in the entire data set. We train regression models to predict the binding affinity from the structural and physicochemical parameters of protein-RNA interfaces. The best fit model with the lowest maximum error is provided with three interface parameters: relative hydrophobicity, conformational change upon binding and relative hydration pattern. This model has been used for predicting the binding affinity on a test data set, generated using mutated structures of yeast aspartyl-tRNA synthetase, for which experimentally determined ΔG values of 40 mutations are available. The predicted ΔGempirical values highly correlate with the experimental observations. The data set provided in this study should be useful for further development of the binding affinity prediction methods. Moreover, the model developed in this study enhances our understanding on the structural basis of protein-RNA binding affinity and provides a platform to engineer protein-RNA interfaces with desired affinity.

QRNAS: software tool for refinement of nucleic acid structures.

BACKGROUND: Computational models of RNA 3D structure often present various inaccuracies caused by simplifications used in structure prediction methods, such as template-based modeling or coarse-grained simulations. To obtain a high-quality model, the preliminary RNA structural model needs to be refined, taking into account atomic interactions. The goal of the refinement is not only to improve the local quality of the model but to bring it globally closer to the true structure. RESULTS: We present QRNAS, a software tool for fine-grained refinement of nucleic acid structures, which is an extension of the AMBER simulation method with additional restraints. QRNAS is capable of handling RNA, DNA, chimeras, and hybrids thereof, and enables modeling of nucleic acids containing modified residues. CONCLUSIONS: We demonstrate the ability of QRNAS to improve the quality of models generated with different methods. QRNAS was able to improve MolProbity scores of NMR structures, as well as of computational models generated in the course of the RNA-Puzzles experiment. The overall geometry improvement may be associated with increased model accuracy, especially on the level of correctly modeled base-pairs, but the systematic improvement of root mean square deviation to the reference structure should not be expected. The method has been integrated into a computational modeling workflow, enabling improved RNA 3D structure prediction.

Computational modeling of RNA 3D structure based on experimental data.

RNA molecules are master regulators of cells. They are involved in a variety of molecular processes: they transmit genetic information, sense cellular signals and communicate responses, and even catalyze chemical reactions. As in the case of proteins, RNA function is dictated by its structure and by its ability to adopt different conformations, which in turn is encoded in the sequence. Experimental determination of high-resolution RNA structures is both laborious and difficult, and therefore the majority of known RNAs remain structurally uncharacterized. To address this problem, predictive computational methods were developed based on the accumulated knowledge of RNA structures determined so far, the physical basis of the RNA folding, and taking into account evolutionary considerations, such as conservation of functionally important motifs. However, all theoretical methods suffer from various limitations, and they are generally unable to accurately predict structures for RNA sequences longer than 100-nt residues unless aided by additional experimental data. In this article, we review experimental methods that can generate data usable by computational methods, as well as computational approaches for RNA structure prediction that can utilize data from experimental analyses. We outline methods and data types that can be potentially useful for RNA 3D structure modeling but are not commonly used by the existing software, suggesting directions for future development.

Dissecting water binding sites at protein-protein interfaces: a lesson from the atomic structures in the Protein Data Bank.

We dissect the protein-protein interfaces into water preservation (WP), water hydration (WH) and water dehydration (WD) sites by comparing the water-mediated hydrogen bonds (H-bond) in the bound and unbound states of the interacting subunits. Upon subunit complexation, if a H-bond between an interface water and a protein polar group is retained, we assign it as WP site; if it is lost, we assign it as WD site and if a new H-bond is created, we assign it as WH site. We find that the density of WD sites is highest followed by WH and WP sites except in antigen and (or) antibody complexes, where the density of WH sites is highest followed by WD and WP sites. Furthermore, we find that WP sites are the most conserved followed by WD and WH sites in all class of complexes except in antigen and (or) antibody complexes, where WD sites are the most conserved followed by WH and WP sites. A significant number of WP and WH sites are involved in water bridges that stabilize the subunit interactions. At WH sites, the residues involved in water bridges are significantly better conserved than the other residues. However, no such difference is observed at WP sites. Interestingly, WD sites are generally replaced with direct H-bonds upon subunit complexation. Significantly, we observe many water-mediated H-bonds remain preserved in spite of large conformational changes upon subunit complexation. These findings have implications in predicting and engineering water binding sites at protein-protein interfaces.

Identification and characterization of differentially expressed Phaseolus vulgaris miRNAs and their targets during mungbean yellow mosaic India virus infection reveals new insight into Phaseolus-MYMIV interaction.

Phaseolus vulgaris is an economically important legume in tropical and subtropical regions of Asia, Africa, Latin-America and parts of USA and Europe. However, its production gets severely affected by mungbean yellow mosaic India virus (MYMIV). We aim to identify and characterize differentially expressed miRNAs during MYMIV-infection in P. vulgaris. A total of 422 miRNAs are identified of which 292 are expressed in both MYMIV-treated and mock-treated samples, 109 are expressed only in MYMIV-treated and 21 are expressed only in mock-treated samples. Selected up- and down-regulated miRNAs are validated by RT-qPCR. 3367 target ORFs are identified for 270 miRNAs. Selected targets are validated by 5' RLM-RACE. Differentially expressed miRNAs regulate transcription factors and are involved in improving stress tolerance to MYMIV. These findings will provide an insight into the role of miRNAs during MYMIV infection in P. vulgaris in particular and during any biotic stress conditions in Leguminosae family in general.

Bioinformatics Tools and Benchmarks for Computational Docking and 3D Structure Prediction of RNA-Protein Complexes.

RNA-protein (RNP) interactions play essential roles in many biological processes, such as regulation of co-transcriptional and post-transcriptional gene expression, RNA splicing, transport, storage and stabilization, as well as protein synthesis. An increasing number of RNP structures would aid in a better understanding of these processes. However, due to the technical difficulties associated with experimental determination of macromolecular structures by high-resolution methods, studies on RNP recognition and complex formation present significant challenges. As an alternative, computational prediction of RNP interactions can be carried out. Structural models obtained by theoretical predictive methods are, in general, less reliable compared to models based on experimental measurements but they can be sufficiently accurate to be used as a basis for to formulating functional hypotheses. In this article, we present an overview of computational methods for 3D structure prediction of RNP complexes. We discuss currently available methods for macromolecular docking and for scoring 3D structural models of RNP complexes in particular. Additionally, we also review benchmarks that have been developed to assess the accuracy of these methods.

Genome-wide identification of miRNAs and lncRNAs in Cajanus cajan.

BACKGROUND: Non-coding RNAs (ncRNAs) are important players in the post transcriptional regulation of gene expression (PTGR). On one hand, microRNAs (miRNAs) are an abundant class of small ncRNAs (~22nt long) that negatively regulate gene expression at the levels of messenger RNAs stability and translation inhibition, on the other hand, long ncRNAs (lncRNAs) are a large and diverse class of transcribed non-protein coding RNA molecules (> 200nt) that play both up-regulatory as well as down-regulatory roles at the transcriptional level. Cajanus cajan, a leguminosae pulse crop grown in tropical and subtropical areas of the world, is a source of high value protein to vegetarians or very poor populations globally. Hence, genome-wide identification of miRNAs and lncRNAs in C. cajan is extremely important to understand their role in PTGR with a possible implication to generate improve variety of crops. RESULTS: We have identified 616 mature miRNAs in C. cajan belonging to 118 families, of which 578 are novel and not reported in MirBase21. A total of 1373 target sequences were identified for 180 miRNAs. Of these, 298 targets were characterized at the protein level. Besides, we have also predicted 3919 lncRNAs. Additionally, we have identified 87 of the predicted lncRNAs to be targeted by 66 miRNAs. CONCLUSIONS: miRNA and lncRNAs in plants are known to control a variety of traits including yield, quality and stress tolerance. Owing to its agricultural importance and medicinal value, the identified miRNA, lncRNA and their targets in C. cajan may be useful for genome editing to improve better quality crop. A thorough understanding of ncRNA-based cellular regulatory networks will aid in the improvement of C. cajan agricultural traits.

A non-redundant protein-RNA docking benchmark version 2.0.

We present an updated version of the protein-RNA docking benchmark, which we first published four years back. The non-redundant protein-RNA docking benchmark version 2.0 consists of 126 test cases, a threefold increase in number compared to its previous version. The present version consists of 21 unbound-unbound cases, of which, in 12 cases, the unbound RNAs are taken from another complex. It also consists of 95 unbound-bound cases where only the protein is available in the unbound state. Besides, we introduce 10 new bound-unbound cases where only the RNA is found in the unbound state. Based on the degree of conformational change of the interface residues upon complex formation the benchmark is classified into 72 rigid-body cases, 25 semiflexible cases and 19 full flexible cases. It also covers a wide range of conformational flexibility including small side chain movement to large domain swapping in protein structures as well as flipping and restacking in RNA bases. This benchmark should provide the docking community with more test cases for evaluating rigid-body as well as flexible docking algorithms. Besides, it will also facilitate the development of new algorithms that require large number of training set. The protein-RNA docking benchmark version 2.0 can be freely downloaded from http://www.csb.iitkgp.ernet.in/applications/PRDBv2. Proteins 2017; 85:256-267. © 2016 Wiley Periodicals, Inc.

Probing binding hot spots at protein-RNA recognition sites.

We use evolutionary conservation derived from structure alignment of polypeptide sequences along with structural and physicochemical attributes of protein-RNA interfaces to probe the binding hot spots at protein-RNA recognition sites. We find that the degree of conservation varies across the RNA binding proteins; some evolve rapidly compared to others. Additionally, irrespective of the structural class of the complexes, residues at the RNA binding sites are evolutionary better conserved than those at the solvent exposed surfaces. For recognitions involving duplex RNA, residues interacting with the major groove are better conserved than those interacting with the minor groove. We identify multi-interface residues participating simultaneously in protein-protein and protein-RNA interfaces in complexes where more than one polypeptide is involved in RNA recognition, and show that they are better conserved compared to any other RNA binding residues. We find that the residues at water preservation site are better conserved than those at hydrated or at dehydrated sites. Finally, we develop a Random Forests model using structural and physicochemical attributes for predicting binding hot spots. The model accurately predicts 80% of the instances of experimental ΔΔG values in a particular class, and provides a stepping-stone towards the engineering of protein-RNA recognition sites with desired affinity.

Targeting Non-Coding RNAs in Plants with the CRISPR-Cas Technology is a Challenge yet Worth Accepting.

Non-coding RNAs (ncRNAs) have emerged as versatile master regulator of biological functions in recent years. MicroRNAs (miRNAs) are small endogenous ncRNAs of 18-24 nucleotides in length that originates from long self-complementary precursors. Besides their direct involvement in developmental processes, plant miRNAs play key roles in gene regulatory networks and varied biological processes. Alternatively, long ncRNAs (lncRNAs) are a large and diverse class of transcribed ncRNAs whose length exceed that of 200 nucleotides. Plant lncRNAs are transcribed by different RNA polymerases, showing diverse structural features. Plant lncRNAs also are important regulators of gene expression in diverse biological processes. There has been a breakthrough in the technology of genome editing, the CRISPR-Cas9 (clustered regulatory interspaced short palindromic repeats/CRISPR-associated protein 9) technology, in the last decade. CRISPR loci are transcribed into ncRNA and eventually form a functional complex with Cas9 and further guide the complex to cleave complementary invading DNA. The CRISPR-Cas technology has been successfully applied in model plants such as Arabidopsis and tobacco and important crops like wheat, maize, and rice. However, all these studies are focused on protein coding genes. Information about targeting non-coding genes is scarce. Hitherto, the CRISPR-Cas technology has been exclusively used in vertebrate systems to engineer miRNA/lncRNAs, but it is still relatively unexplored in plants. While briefing miRNAs, lncRNAs and applications of the CRISPR-Cas technology in human and animals, this review essentially elaborates several strategies to overcome the challenges of applying the CRISPR-Cas technology in editing ncRNAs in plants and the future perspective of this field.

Computational prediction of miRNAs and their targets in Phaseolus vulgaris using simple sequence repeat signatures.

BACKGROUND: MicroRNAs (miRNAs) are endogenous, noncoding, short RNAs directly involved in regulating gene expression at the post-transcriptional level. In spite of immense importance, limited information of P. vulgaris miRNAs and their expression patterns prompted us to identify new miRNAs in P. vulgaris by computational methods. Besides conventional approaches, we have used the simple sequence repeat (SSR) signatures as one of the prediction parameter. Moreover, for all other parameters including normalized Shannon entropy, normalized base pairing index and normalized base-pair distance, instead of taking a fixed cut-off value, we have used 99% probability range derived from the available data. RESULTS: We have identified 208 mature miRNAs in P. vulgaris belonging to 118 families, of which 201 are novel. 97 of the predicted miRNAs in P. vulgaris were validated with the sequencing data obtained from the small RNA sequencing of P. vulgaris. Randomly selected predicted miRNAs were also validated using qRT-PCR. A total of 1305 target sequences were identified for 130 predicted miRNAs. Using 80% sequence identity cut-off, proteins coded by 563 targets were identified. The computational method developed in this study was also validated by predicting 229 miRNAs of A. thaliana and 462 miRNAs of G. max, of which 213 for A. thaliana and 397 for G. max are existing in miRBase 20. CONCLUSIONS: There is no universal SSR that is conserved among all precursors of Viridiplantae, but conserved SSR exists within a miRNA family and is used as a signature in our prediction method. Prediction of known miRNAs of A. thaliana and G. max validates the accuracy of our method. Our findings will contribute to the present knowledge of miRNAs and their targets in P. vulgaris. This computational method can be applied to any species of Viridiplantae for the successful prediction of miRNAs and their targets.