The favorable IFNL3 genotype escapes mRNA decay mediated by AU-rich elements and hepatitis C virus-induced microRNAs.

IFNL3, which encodes interferon-λ3 (IFN-λ3), has received considerable attention in the hepatitis C virus (HCV) field, as many independent genome-wide association studies have identified a strong association between polymorphisms near IFNL3 and clearance of HCV. However, the mechanism underlying this association has remained elusive. In this study, we report the identification of a functional polymorphism (rs4803217) in the 3' untranslated region (UTR) of IFNL3 mRNA that dictated transcript stability. We found that this polymorphism influenced AU-rich element (ARE)-mediated decay (AMD) of IFNL3 mRNA, as well as the binding of HCV-induced microRNAs during infection. Together these pathways mediated robust repression of the unfavorable IFNL3 polymorphism. Our data reveal a previously unknown mechanism by which HCV attenuates the antiviral response and indicate new potential therapeutic targets for HCV treatment.

Structural characterization of a new subclass of panicum mosaic virus-like 3' cap-independent translation enhancer.

Canonical eukaryotic mRNA translation requires 5'cap recognition by initiation factor 4E (eIF4E). In contrast, many positive-strand RNA virus genomes lack a 5'cap and promote translation by non-canonical mechanisms. Among plant viruses, PTEs are a major class of cap-independent translation enhancers located in/near the 3'UTR that recruit eIF4E to greatly enhance viral translation. Previous work proposed a single form of PTE characterized by a Y-shaped secondary structure with two terminal stem-loops (SL1 and SL2) atop a supporting stem containing a large, G-rich asymmetric loop that forms an essential pseudoknot (PK) involving C/U residues located between SL1 and SL2. We found that PTEs with less than three consecutive cytidylates available for PK formation have an upstream stem-loop that forms a kissing loop interaction with the apical loop of SL2, important for formation/stabilization of PK. PKs found in both subclasses of PTE assume a specific conformation with a hyperreactive guanylate (G*) in SHAPE structure probing, previously found critical for binding eIF4E. While PTE PKs were proposed to be formed by Watson-Crick base-pairing, alternative chemical probing and 3D modeling indicate that the Watson-Crick faces of G* and an adjacent guanylate have high solvent accessibilities. Thus, PTE PKs are likely composed primarily of non-canonical interactions.

Dynamic bulge nucleotides in the KSHV PAN ENE triple helix provide a unique binding platform for small molecule ligands.

Cellular and virus-coded long non-coding (lnc) RNAs support multiple roles related to biological and pathological processes. Several lncRNAs sequester their 3' termini to evade cellular degradation machinery, thereby supporting disease progression. An intramolecular triplex involving the lncRNA 3' terminus, the element for nuclear expression (ENE), stabilizes RNA transcripts and promotes persistent function. Therefore, such ENE triplexes, as presented here in Kaposi's sarcoma-associated herpesvirus (KSHV) polyadenylated nuclear (PAN) lncRNA, represent targets for therapeutic development. Towards identifying novel ligands targeting the PAN ENE triplex, we screened a library of immobilized small molecules and identified several triplex-binding chemotypes, the tightest of which exhibits micromolar binding affinity. Combined biophysical, biochemical, and computational strategies localized ligand binding to a platform created near a dinucleotide bulge at the base of the triplex. Crystal structures of apo (3.3 Å) and ligand-soaked (2.5 Å) ENE triplexes, which include a stabilizing basal duplex, indicate significant local structural rearrangements within this dinucleotide bulge. MD simulations and a modified nucleoside analog interference technique corroborate the role of the bulge and the base of the triplex in ligand binding. Together with recently discovered small molecules that reduce nuclear MALAT1 lncRNA levels by engaging its ENE triplex, our data supports the potential of targeting RNA triplexes with small molecules.

Virtual Screening of Hepatitis B Virus Pre-Genomic RNA as a Novel Therapeutic Target.

The global burden imposed by hepatitis B virus (HBV) infection necessitates the discovery and design of novel antiviral drugs to complement existing treatments. One attractive and underexploited therapeutic target is ε, an ~85-nucleotide (nt) cis-acting regulatory stem-loop RNA located at the 3'- and 5'-ends of the pre-genomic RNA (pgRNA). Binding of the 5'-end ε to the viral polymerase protein (P) triggers two early events in HBV replication: pgRNA and P packaging and reverse transcription. Our recent solution nuclear magnetic resonance spectroscopy structure of ε permits structure-informed drug discovery efforts that are currently lacking for P. Here, we employ a virtual screen against ε using a Food and Drug Administration (FDA)-approved compound library, followed by in vitro binding assays. This approach revealed that the anti-hepatitis C virus drug Daclatasvir is a selective ε-targeting ligand. Additional molecular dynamics simulations demonstrated that Daclatasvir targets ε at its flexible 6-nt priming loop (PL) bulge and modulates its dynamics. Given the functional importance of the PL, our work supports the notion that targeting ε dynamics may be an effective anti-HBV therapeutic strategy.

Stealth oxime ether lipid vesicles promote delivery of functional DsiRNA in human lung cancer A549 tumor bearing mouse xenografts.

We previously reported that hydroxylated oxime ether lipids (OELs) efficiently deliver functional Dicer substrate siRNAs (DsiRNAs) in cells. Here, we explored in vivo utility of these OELs, using OEL4 as a prototype and report that surface modification of the OEL4 formulations was essential for their in vivo applications. These surface-modified OEL4 formulations were developed by inclusion of various PEGylated lipids. The vesicle stability and gene knock-down were dependent on the PEG chain length. OEL4 containing DSPE-PEG350 and DSPE-PEG1000 (surprisingly not DSPE2000) promoted gene silencing in cells. In vivo studies demonstrated that OEL4 vesicles formulated using 3 mol% DSPE-PEG350 accumulate in human lung cancer (A549-luc2) xenografts in mice and exhibit a significant increase in tumor to liver ratios. These vesicles also showed a statistically significant reduction of luciferase signal in tumors compared to untreated mice. Taken together, the scalable OEL4:DSPE-PEG350 formulation serves as a novel candidate for delivery of RNAi therapeutics.

Predicting RNA SHAPE scores with deep learning.

Secondary structure prediction approaches rely typically on models of equilibrium free energies that are themselves based on in vitro physical chemistry. Recent transcriptome-wide experiments of in vivo RNA structure based on SHAPE-MaP experiments provide important information that may make it possible to extend current in vitro-based RNA folding models in order to improve the accuracy of computational RNA folding simulations with respect to the experimentally measured in vivo RNA secondary structure. Here we present a machine learning approach that utilizes RNA secondary structure prediction results and nucleotide sequence in order to predict in vivo SHAPE scores. We show that this approach has a higher Pearson correlation coefficient with experimental SHAPE scores than thermodynamic folding. This could be an important step towards augmenting experimental results with computational predictions and help with RNA secondary structure predictions that inherently take in-vivo folding properties into account.

Ribosomal frameshifting in the CCR5 mRNA is regulated by miRNAs and the NMD pathway.

Programmed -1 ribosomal frameshift (-1 PRF) signals redirect translating ribosomes to slip back one base on messenger RNAs. Although well characterized in viruses, how these elements may regulate cellular gene expression is not understood. Here we describe a -1 PRF signal in the human mRNA encoding CCR5, the HIV-1 co-receptor. CCR5 mRNA-mediated -1 PRF is directed by an mRNA pseudoknot, and is stimulated by at least two microRNAs. Mapping the mRNA-miRNA interaction suggests that formation of a triplex RNA structure stimulates -1 PRF. A -1 PRF event on the CCR5 mRNA directs translating ribosomes to a premature termination codon, destabilizing it through the nonsense-mediated mRNA decay pathway. At least one additional mRNA decay pathway is also involved. Functional -1 PRF signals that seem to be regulated by miRNAs are also demonstrated in mRNAs encoding six other cytokine receptors, suggesting a novel mode through which immune responses may be fine-tuned in mammalian cells.

Photoactivation of sulfonated polyplexes enables localized gene silencing by DsiRNA in breast cancer cells.

Translation potential of RNA interference nanotherapeutics remains challenging due to in vivo off-target effects and poor endosomal escape. Here, we developed novel polyplexes for controlled intracellular delivery of dicer substrate siRNA, using a light activation approach. Sulfonated polyethylenimines covalently linked to pyropheophorbide-α for photoactivation and bearing modified amines (sulfo-pyro-PEI) for regulated endosomal escape were investigated. Gene knock-down by the polymer-complexed DsiRNA duplexes (siRNA-NPs) was monitored in breast cancer cells. Surprisingly, sulfo-pyro-PEI/siRNA-NPs failed to downregulate the PLK1 or eGFP proteins. However, photoactivation of these cell associated-polyplexes with a 661-nm laser clearly restored knock-down of both proteins. In contrast, protein down-regulation by non-sulfonated pyro-PEI/siRNA-NPs occurred without any laser treatments, indicating cytoplasmic disposition of DsiRNA followed a common intracellular release mechanism. Therefore, sulfonated pyro-PEI holds potential as a unique trap and release light-controlled delivery platform for on-demand gene silencing bearing minimal off target effects.

RiboSketch: versatile visualization of multi-stranded RNA and DNA secondary structure.

Summary: Creating clear, visually pleasing 2D depictions of RNA and DNA strands and their interactions is important to facilitate and communicate insights related to nucleic acid structure. Here we present RiboSketch, a secondary structure image production application that enables the visualization of multistranded structures via layout algorithms, comprehensive editing capabilities, and a multitude of simulation modes. These interactive features allow RiboSketch to create publication quality diagrams for structures with a wide range of composition, size and complexity. The program may be run in any web browser without the need for installation, or as a standalone Java application. Availability and implementation: https://rnastructure.cancer.gov/ribosketch.

Pre-clinical compartmental pharmacokinetic modeling of 2-[1-hexyloxyethyl]-2-devinyl pyropheophorbide-a (HPPH) as a photosensitizer in rat plasma by validated HPLC method.

A second-generation chlorin-based photosensitizer, 2-[1-hexyloxyethyl]-2-devinyl pyropheophorbide-a (HPPH) has shown tremendous therapeutic potential in clinical trials in the treatment of esophageal cancer. Herein, we have developed and validated a bioanalytical method for estimation of HPPH in rat plasma using High Performance Liquid Chromatography (HPLC) with a photo diode array (PDA) detector. The method was applied for carrying out pharmacokinetic study of HPPH. Further pharmacokinetic modeling was carried out to understand the compartment kinetics of HPPH. The developed method was fully validated as per the United States Food and Drug Administration (US-FDA) guidelines for bioanalytical method validation. The linearity of the method was in the range of 250-8000 ng mL-1, and the plasma recovery was found to be 70%. Pharmacokinetic parameters were evaluated and compared via non-compartment analysis and compartment modeling after the intravenous (i.v.) bolus administration in rats using Phoenix WinNonlin 8.0 (Certara™, USA). From the obtained results, we hypothesize that the HPPH complies with two compartmental pharmacokinetic model. Furthermore, it was observed that HPPH has the rapid distribution from the central compartment to peripheral compartment along with slow elimination from peripheral compartment.

RNA2Drawer: geometrically strict drawing of nucleic acid structures with graphical structure editing and highlighting of complementary subsequences.

RNA structure prediction programs remain imperfect and many substructures are still identified by manual exploration, which is most efficiently conducted within an RNA structure drawing program. However, most nucleic acid structure drawing programs have limited capability for structure modification (i.e., breaking and forming new bonds between bases), often requiring that the structure notation be textually edited. RNA2Drawer was developed to allow for graphical structure editing while maintaining the geometry of a drawing (e.g., ellipsoid loops, stems with evenly stacked base pairs) throughout structural changes and manual adjustments to the layout by the user. In addition, the program allows for annotations such as colouring and circling of bases and drawing of tertiary interactions (e.g., pseudoknots). RNA2Drawer can also draw commonly desired elements such as an optionally flattened outermost loop and assists structure editing by automatically highlighting complementary subsequences, which facilitates the discovery of potentially new and alternative pairings, particularly tertiary pairings over long-distances, which are biologically critical in the genomes of many RNA viruses and cannot be accurately predicted by current structure prediction programs. Additionally, RNA2Drawer outputs drawings either as PNG files, or as PPTX and SVG files, such that every object of a drawing (e.g., bases, bonds) is an individual PPTX or SVG object, allowing for further manipulation in Microsoft PowerPoint or a vector graphics editor such as Adobe Illustrator. PowerPoint is the standard for presentations and is often used to create figures for publications, and RNA2Drawer is the first program to export drawings as PPTX files.

Functionally-interdependent shape-switching nanoparticles with controllable properties.

We introduce a new concept that utilizes cognate nucleic acid nanoparticles which are fully complementary and functionally-interdependent to each other. In the described approach, the physical interaction between sets of designed nanoparticles initiates a rapid isothermal shape change which triggers the activation of multiple functionalities and biological pathways including transcription, energy transfer, functional aptamers and RNA interference. The individual nanoparticles are not active and have controllable kinetics of re-association and fine-tunable chemical and thermodynamic stabilities. Computational algorithms were developed to accurately predict melting temperatures of nanoparticles of various compositions and trace the process of their re-association in silico. Additionally, tunable immunostimulatory properties of described nanoparticles suggest that the particles that do not induce pro-inflammatory cytokines and high levels of interferons can be used as scaffolds to carry therapeutic oligonucleotides, while particles with strong interferon and mild pro-inflammatory cytokine induction may qualify as vaccine adjuvants. The presented concept provides a simple, cost-effective and straightforward model for the development of combinatorial regulation of biological processes in nucleic acid nanotechnology.

Dynamic Behavior of RNA Nanoparticles Analyzed by AFM on a Mica/Air Interface.

RNA is an attractive biopolymer for engineering self-assembling materials suitable for biomedical applications. Previously, programmable hexameric RNA rings were developed for the controlled delivery of up to six different functionalities. To increase the potential for functionalization with little impact on nanoparticle topology, we introduce gaps into the double-stranded regions of the RNA rings. Molecular dynamic simulations are used to assess the dynamic behavior and the changes in the flexibility of novel designs. The changes suggested by simulations, however, cannot be clearly confirmed by the conventional techniques such as nondenaturing polyacrylamide gel electrophoresis (native-PAGE) and dynamic light scattering (DLS). Also, an in vitro analysis in primary cultures of human peripheral blood mononuclear cells does not reveal any discrepancy in the immunological recognition of new assemblies. To address these deficiencies, we introduce a computer-assisted quantification strategy. This strategy is based on an algorithmic atomic force microscopy (AFM)-resolved deformation analysis of the RNA nanoparticles studied on a mica/air interface. We validate this computational method by manual image analysis and fitting it to the simulation-predicted results. The presented nanoparticle modification strategy and subsequent AFM-based analysis are anticipated to provide a broad spectrum approach for the future development of nucleic acid-based nanotechnology.

Design and biological activity of novel stealth polymeric lipid nanoparticles for enhanced delivery of hydrophobic photodynamic therapy drugs.

Advances in in vivo stability and preferential tumor uptake of cancer nanomedicine are warranted for effective chemotherapy. Here, we describe a novel nanoformulation using an unconventional polymeric tubule-forming phospholipid, DC8,9PC. We report that DC8,9PC transitions to stable vesicles (LNPs) in the presence of PEGylated lipid (DSPE-PEG2000); the resulting DC8,9PC:DSPE-PEG2000 LNPs efficiently included a hydrophobic PDT drug, HPPH. Remarkably, these LNPs incorporated unusually high DSPE-PEG2000 concentrations; LNP10-HPPH and LNP20-HPPH (10 & 20 mol% PEGylated lipid, respectively) exhibited >90% serum stability at 37 °C. Increased PEGylation in the LNPs correlated with enhanced tumor accumulation in intravenously injected HT29 tumor mouse xenographs. Colon-26 bearing BALB/c mice, intravenously injected with LNP20-HPPH showed superior PDT efficacy and animal survival (no tumor recurrence up to 100 days) as compared to a formulation currently used in clinical trials. Taken together, we present a simple stealth binary lipid nanosystem with enhanced efficiency of tumor accumulation and superior therapeutic efficacy.

RNA⁻Protein Interactions Prevent Long RNA Duplex Formation: Implications for the Design of RNA-Based Therapeutics.

Cells frequently simultaneously express RNAs and cognate antisense transcripts without necessarily leading to the formation of RNA duplexes. Here, we present a novel transcriptome-wide experimental approach to ascertain the presence of accessible double-stranded RNA structures based on sequencing of RNA fragments longer than 18 nucleotides that were not degraded by single-strand cutting nucleases. We applied this approach to four different cell lines with respect to three different treatments (native cell lysate, removal of proteins, and removal of ribosomal RNA and proteins). We found that long accessible RNA duplexes were largely absent in native cell lysates, while the number of RNA duplexes was dramatically higher when proteins were removed. The majority of RNA duplexes involved ribosomal transcripts. The duplex formation between different non-ribosomal transcripts appears to be largely of a stochastic nature. These results suggest that cells are-via RNA-binding proteins-mostly devoid of long RNA duplexes, leading to low "noise" in the molecular patterns that are utilized by the innate immune system. These findings have implications for the design of RNA interference (RNAi)-based therapeutics by imposing structural constraints on designed RNA complexes that are intended to have specific properties with respect to Dicer cleavage and target gene downregulation.

Combined single molecule experimental and computational approaches for understanding the unfolding pathway of a viral translation enhancer that participates in a conformational switch.

How plus-strand [+]RNA virus genomes transition from translation templates to replication templates is a matter of much speculation. We have previously proposed that, for Turnip crinkle virus, binding of the encoded RNA-dependent RNA polymerase (RdRp) to the 3'UTR of the [+]RNA template promotes a regional wide-spread conformational switch to an alternative structure that disassembles the cap-independent translation enhancer (CITE) in the 3'UTR. The active 3'CITE folds into a tRNA-like T-shaped structure (TSS) that binds to 80S ribosomes and 60S subunits in the P-site. In this Point-of-View, we discuss the history of our research on the TSS and our recent report combining coarse level single molecule force spectroscopy (optical tweezers) with fine-grain computer simulations of this experimental process and biochemical approaches to obtain a detailed understanding of how RdRp binding in the TSS vicinity might lead to an extensive rearrangement of the RNA structure.

Ring Catalog: A resource for designing self-assembling RNA nanostructures.

Designing self-assembling RNA ring structures based on known 3D structural elements connected via linker helices is a challenging task due to the immense number of motif combinations, many of which do not lead to ring-closure. We describe an in silico solution to this design problem by combinatorial assembly of RNA 3-way junctions, bulges, and kissing loops, and tabulating the cases that lead to ring formation. The solutions found are made available in the form of a web-accessible Ring Catalog. As an example of a potential use of this resource, we chose a predicted RNA square structure consisting of five RNA strands and demonstrate experimentally that the self-assembly of those five strands leads to the formation of a square-like complex. This is a demonstration of a novel "design by catalog" approach to RNA nano-structure generation. The URL https://rnajunction.ncifcrf.gov/ringdb can be used to access the resource.

Multistrand Structure Prediction of Nucleic Acid Assemblies and Design of RNA Switches.

RNA is an attractive material for the creation of molecular logic gates that release programmed functionalities only in the presence of specific molecular interaction partners. Here we present HyperFold, a multistrand RNA/DNA structure prediction approach for predicting nucleic acid complexes that can contain pseudoknots. We show that HyperFold also performs competitively compared to other published folding algorithms. We performed a large variety of RNA/DNA hybrid reassociation experiments for different concentrations, DNA toehold lengths, and G+C content and find that the observed tendencies for reassociation correspond well to computational predictions. Importantly, we apply this method to the design and experimental verification of a two-stranded RNA molecular switch that upon binding to a single-stranded RNA toehold disease-marker trigger mRNA changes its conformation releasing an shRNA-like Dicer substrate structure. To demonstrate the concept, connective tissue growth factor (CTGF) mRNA and enhanced green fluorescent protein (eGFP) mRNA were chosen as trigger and target sequences, respectively. In vitro experiments confirm the formation of an RNA switch and demonstrate that the functional unit is being released when the trigger RNA interacts with the switch toehold. The designed RNA switch is shown to be functional in MDA-MB-231 breast cancer cells. Several other switches were also designed and tested. We conclude that this approach has considerable potential because, in principle, it allows the release of an siRNA designed against a gene that differs from the gene that is utilized as a biomarker for a disease state.

The Use of Minimal RNA Toeholds to Trigger the Activation of Multiple Functionalities.

Current work reports the use of single-stranded RNA toeholds of different lengths to promote the reassociation of various RNA-DNA hybrids, which results in activation of multiple split functionalities inside human cells. The process of reassociation is analyzed and followed with a novel computational multistrand secondary structure prediction algorithm and various experiments. All of our previously designed RNA/DNA nanoparticles employed single-stranded DNA toeholds to initiate reassociation. The use of RNA toeholds is advantageous because of the simpler design rules, the shorter toeholds, and the smaller size of the resulting nanoparticles (by up to 120 nucleotides per particle) compared to the same hybrid nanoparticles with single-stranded DNA toeholds. Moreover, the cotranscriptional assemblies result in higher yields for hybrid nanoparticles with ssRNA toeholds.

Co-transcriptional production of RNA-DNA hybrids for simultaneous release of multiple split functionalities.

Control over the simultaneous delivery of different functionalities and their synchronized intracellular activation can greatly benefit the fields of RNA and DNA biomedical nanotechnologies and allow for the production of nanoparticles and various switching devices with controllable functions. We present a system of multiple split functionalities embedded in the cognate pairs of RNA-DNA hybrids which are programmed to recognize each other, re-associate and form a DNA duplex while also releasing the split RNA fragments which upon association regain their original functions. Simultaneous activation of three different functionalities (RNAi, Förster resonance energy transfer and RNA aptamer) confirmed by multiple in vitro and cell culture experiments prove the concept. To automate the design process, a novel computational tool that differentiates between the thermodynamic stabilities of RNA-RNA, RNA-DNA and DNA-DNA duplexes was developed. Moreover, here we demonstrate that besides being easily produced by annealing synthetic RNAs and DNAs, the individual hybrids carrying longer RNAs can be produced by RNA polymerase II-dependent transcription of single-stranded DNA templates.