Improved nearest-neighbor parameters for the stability of RNA/DNA hybrids under a physiological condition

Strong sequence-dependence in RNA/DNA hybrid strand displacement kinetics

Strand displacement reactions underlie dynamic nucleic acid nanotechnology. The kinetic and thermodynamic features of DNA-based displacement reactions are well understood and well predicted by current computational models. By contrast, understanding of RNA/DNA hybrid strand displacement kinetics is limited, restricting the design of increasingly complex RNA/DNA hybrid reaction networks with more tightly regulated dynamics. Given the importance of RNA as a diagnostic biomarker, and its critical role in intracellular processes, this shortfall is particularly limiting for the development of strand displacement-based therapeutics and diagnostics. Herein, we characterise 22 RNA/DNA hybrid strand displacement systems, alongside 11 DNA/DNA systems, varying a range of common design parameters including toehold length and branch migration domain length. We observe  that differences in stability between RNA-DNA hybrids and DNA-DNA duplexes have large effects on strand displacement rates, with rates for equivalent sequences differing by up to 3 orders of magnitude. Crucially, however, this effect is strongly sequence-dependent, with RNA invaders strongly favoured in a system with RNA strands of high purine content, and disfavoured in a system when the RNA strands have low purine content. These results lay the groundwork for more general design principles, allowing for creation of de novo reaction networks with novel complexity while maintaining predictable reaction kinetics.

A mechanistic study on the tolerance of PAM distal end mismatch by SpCas9

The therapeutic application of CRISPR-Cas9 is limited due to its off-target activity. To have a better understanding of this off-target effect, we focused on its mismatch-prone PAM distal end. The off-target activity of SpCas9 depends directly on the nature of mismatches, which in turn results in deviation of the active site of SpCas9 due to structural instability in the RNA-DNA duplex strand. In order to test the hypothesis, we designed an array of mismatched target sites at the PAM distal end and performed in vitro and cell line-based experiments, which showed a strong correlation for Cas9 activity. We found that target sites having multiple mismatches in the 18th to 15th position upstream of the PAM showed no to little activity. For further mechanistic validation, Molecular Dynamics simulations were performed, which revealed that certain mismatches showed elevated root mean square deviation values that can be attributed to conformational instability within the RNA-DNA duplex. Therefore, for successful prediction of the off-target effect of SpCas9, along with complementation-derived energy, the RNA-DNA duplex stability should be taken into account.

Effects of monovalent and divalent cations on the rheology of entangled DNA

In this paper we investigate the effects of varying cation valency and concentration on the rheology of entangled λDNA solutions. We show that monovalent cations moderately increase the viscoelasticty of the solutions mainly by stabilising linear concatenation of λDNA "monomers" via hybridisation of their sticky ends. On the contrary, divalent cations have a far more complex and dramatic effect on the rheology of the solution and we observe evidence of inter-molecular DNA-DNA bridging by Mg2+. We argue that these results may be interesting in the context of dense solutions of single and double stranded DNA, e.g. in vivo or in biotechnology applications such as DNA origami and DNA hydrogels.

In vitro selection of a trans aptamer complex for target-responsive fluorescence activation

BACKGROUND

Most biological molecular complexes consist of multiple functional domains, yet rationally constructing such multifunctional complexes is challenging. Aptamers, the nucleic acid-based functional molecules, can perform multiple tasks including target recognition, conformational changes, and enzymatic activities, while being chemically synthesizable and tunable, and thus provide a basis for engineering enhanced functionalities through combination of multiple units. However, the conventional approach of simply combining aptamer units in a serial manner is susceptible to undesired crosstalk or interference between the aptamer units and to false interactions with non-target molecules; besides, the approach would require additional mechanisms to separate the units if they are desired to function independently. It is clearly a challenge to develop multi-aptamer complexes that preserve independent functions of each unit while avoiding undesired interference and non-specific interactions.

RESULTS

By directly in vitro selecting a 'trans' aptamer complex, we demonstrate that one aptamer unit ('utility module') can remain hidden or 'inactive' until a target analyte triggers the other unit ('sensing module') and separates the two aptamers. Since the operation of the utility module occurs free from the sensing module, unnecessary crosstalk between the two units can be avoided. Because the utility module is kept inactive until separated from the complex, non-specific interactions of the hidden module with noncognate targets can be naturally prevented. In our demonstration, the sensing module was selected to detect serotonin, a clinically important neurotransmitter, and the target-binding-induced structure-switching of the sensing module reveals and activates the utility module that turns on a fluorescence signal. The aptamer complex exhibited a moderately high affinity and an excellent specificity for serotonin with ∼16-fold discrimination against common neurotransmitter molecules, and displayed strong robustness to perturbations in the design, disallowing nonspecific reactions against various challenges.

SIGNIFICANCE

This work represents the first example of a trans aptamer complex that was in vitro selected de novo. The trans aptamer complex selected by our strategy does not require chemical modifications or immediate optimization processes to function, because the complex is directly selected to perform desired functions. This strategy should be applicable to a wide range of functional nucleic acid moieties, which will open up diverse applications in biosensing and molecular therapeutics.

NMR characterization of RNA binding property of the DEAD-box RNA helicase DDX3X and its implications for helicase activity

The DEAD-box RNA helicase (DDX) plays a central role in many aspects of RNA metabolism by remodeling the defined structure of RNA molecules. While a number of structural studies have revealed the atomistic details of the interaction between DDX and RNA ligands, the molecular mechanism of how this molecule unwinds a structured RNA into an unstructured single-stranded RNA (ssRNA) has largely remained elusive. This is due to challenges in structurally characterizing the unwinding intermediate state and the lack of thermodynamic details underlying this process. In this study, we use solution nuclear magnetic resonance (NMR) spectroscopy to characterize the interaction of human DDX3X, a member of the DDX family, with various RNA ligands. Our results show that the inherent binding affinity of DDX3X for ssRNA is significantly higher than that for structured RNA elements. This preferential binding, accompanied by the formation of a domain-closed conformation in complex with ssRNA, effectively stabilizes the denatured ssRNA state and thus underlies the unwinding activity of DDX3X. Our results provide a thermodynamic and structural basis for the DDX function, whereby DDX can recognize and remodel a distinct set of structured RNAs to participate in a wide range of physiological processes.

Unleashing inkjet-printed nanostructured electrodes and battery-free potentiostat for the DNA-based multiplexed detection of SARS-CoV-2 genes

Following the global COVID-19 pandemic triggered by SARS-CoV-2, the need for rapid, specific and cost-effective point-of-care diagnostic solutions remains paramount. Even though COVID-19 is no longer a public health emergency, the disease still poses a global threat leading to deaths, and it continues to change with the risk of new variants emerging causing a new surge in cases and deaths. Here, we address the urgent need for rapid, cost-effective and point-of-care diagnostic solutions for SARS-CoV-2. We propose a multiplexed DNA-based sensing platform that utilizes inkjet-printed nanostructured gold electrodes and an inkjet-printed battery-free near-field communication (NFC) potentiostat for the simultaneous quantitative detection of two SARS-CoV-2 genes, the ORF1ab and the N gene. The detection strategy based on the formation of an RNA-DNA sandwich structure leads to a highly specific electrochemical output. The inkjet-printed nanostructured gold electrodes providing a large surface area enable efficient binding and increase the sensitivity. The inkjet-printed battery-free NFC potentiostat enables rapid measurements and real-time data analysis via a smartphone application, making the platform accessible and portable. With the advantages of speed (5 min), simplicity, sensitivity (low pM range, ∼450% signal gain) and cost-effectiveness, the proposed platform is a promising alternative for point-of-care diagnostics and high-throughput analysis that complements the COVID-19 diagnostic toolkit.

Coarse-grained modeling of DNA-RNA hybrids

We introduce oxNA, a new model for the simulation of DNA-RNA hybrids that is based on two previously developed coarse-grained models-oxDNA and oxRNA. The model naturally reproduces the physical properties of hybrid duplexes, including their structure, persistence length, and force-extension characteristics. By parameterizing the DNA-RNA hydrogen bonding interaction, we fit the model's thermodynamic properties to experimental data using both average-sequence and sequence-dependent parameters. To demonstrate the model's applicability, we provide three examples of its use-calculating the free energy profiles of hybrid strand displacement reactions, studying the resolution of a short R-loop, and simulating RNA-scaffolded wireframe origami.

A FRET-Based Assay and Computational Tools to Quantify Enzymatic Rates and Explore the Mechanisms of RNA Deadenylases in Heterogeneous Environments

We developed a medium-throughput assay that can measure the time-dependent distribution of RNA products generated as a deadenylase degrades a polyadenosine (poly(A)) RNA tract, thereby providing insight into the mechanism of deadenylation. Importantly, this assay can be performed in both homogeneous and heterogeneous environments without relying on gel electrophoresis of RNA products or coupled enzymatic reactions that indirectly report on the RNA distribution through the detection of freed adenosine monophosphate. In parallel, we have established an open-source, Python-based command-line software package, deadenylationkinetics, that can be used to numerically simulate and/or fit the datasets afforded by our assay with different deadenylation mechanisms to determine the most likely case and estimate the associated rate constants. In this chapter, we detail the implementation of our method and the quantification of poly(A) RNA binding and degradation kinetics in application to a truncated version of CNOT7 from the CCR4-NOT deadenylation complex, which serves as a model deadenylase with enhanced activity.

Application of the weighted histogram method for calculating the thermodynamic parameters of the formation of oligodeoxyribonucleotide duplexes

To date, many derivatives and analogs of nucleic acids (NAs) have been developed. Some of them have found uses in scientific research and biomedical applications. Their effective use is based on the data about their properties. Some of the most important physicochemical properties of oligonucleotides are thermodynamic parameters of the formation of their duplexes with DNA and RNA. These parameters can be calculated only for a few NA derivatives: locked NAs, bridged oligonucleotides, and peptide NAs. Existing predictive approaches are based on an analysis of experimental data and the consequent construction of predictive models. The ongoing pilot studies aimed at devising methods for predicting the properties of NAs by computational modeling techniques are based only on knowledge about the structure of oligonucleotides. In this work, we studied the applicability of the weighted histogram analysis method (WHAM) in combination with umbrella sampling to the calculation of thermodynamic parameters of DNA duplex formation (changes in enthalpy ΔH°, entropy ΔS°, and Gibbs free energy ΔG° 37). A procedure was designed involving WHAM for calculating the hybridization properties of oligodeoxyribonucleotides. Optimal parameters for modeling and calculation of thermodynamic parameters were determined. The feasibility of calculation of ΔH°, ΔS°, and ΔG° 37 was demonstrated using a representative sample of 21 oligonucleotides 4-16 nucleotides long with a GC content of 14-100 %. Error of the calculation of the thermodynamic parameters was 11.4, 12.9, and 11.8 % for ΔH°, ΔS°, and ΔG° 37, respectively, and the melting temperature was predicted with an average error of 5.5 °C. Such high accuracy of computations is comparable with the accuracy of the experimental approach and of other methods for calculating the energy of NA duplex formation. In this paper, the use of WHAM for computation of the energy of DNA duplex formation was systematically investigated for the first time. Our results show that a reliable calculation of the hybridization parameters of new NA derivatives is possible, including derivatives not yet synthesized. This work opens up new horizons for a rational design of constructs based on NAs for solving problems in biomedicine and biotechnology.

In-Cell Stability Prediction of RNA/DNA Hybrid Duplexes for Designing Oligonucleotides Aimed at Therapeutics

In cells, the formation of RNA/DNA hybrid duplexes regulates gene expression and modification. The environment inside cellular organelles is heterogeneously crowded with high concentrations of biomolecules that affect the structure and stability of RNA/DNA hybrid duplexes. However, the detailed environmental effects remain unclear. Therefore, the mechanistic details of the effect of such molecular crowding were investigated at the molecular level by using thermodynamic and nuclear magnetic resonance analyses, revealing structure-dependent destabilization of the duplexes under crowded conditions. The transition from B- to A-like hybrid duplexes due to a change in conformation of the DNA strand guided by purine-pyrimidine asymmetry significantly increased the hydration number, which resulted in greater destabilization by the addition of cosolutes. By quantifying the individual contributions of environmental factors and the bulk structure of the duplex, we developed a set of parameters that predict the stability of hybrid duplexes with conformational dissimilarities under diverse crowding conditions. A comparison of the effects of environmental conditions in living cells and in vitro crowded solutions on hybrid duplex formation using the Förster resonance energy transfer technique established the applicability of our parameters to living cells. Moreover, our derived parameters can be used to estimate the efficiency of transcriptional inhibition, genome editing, and silencing techniques in cells. This supports the usefulness of our parameters for the visualization of cellular mechanisms of gene expression and the development of nucleic acid-based therapeutics targeting different cells.

Unmodificated stepless regulation of CRISPR/Cas12a multi-performance

As CRISPR technology is promoted to more fine-divided molecular biology applications, its inherent performance finds it increasingly difficult to cope with diverse needs in these different fields, and how to more accurately control the performance has become a key issue to develop CRISPR technology to a new stage. Herein, we propose a CRISPR/Cas12a regulation strategy based on the powerful programmability of nucleic acid nanotechnology. Unlike previous difficult and rigid regulation of core components Cas nuclease and crRNA, only a simple switch of different external RNA accessories is required to change the reaction kinetics or thermodynamics, thereby finely and almost steplessly regulating multi-performance of CRISPR/Cas12a including activity, speed, specificity, compatibility, programmability and sensitivity. In particular, the significantly improved specificity is expected to mark advance the accuracy of molecular detection and the safety of gene editing. In addition, this strategy was applied to regulate the delayed activation of Cas12a, overcoming the compatibility problem of the one-pot assay without any physical separation or external stimulation, and demonstrating great potential for fine-grained control of CRISPR. This simple but powerful CRISPR regulation strategy without any component modification has pioneering flexibility and versatility, and will unlock the potential for deeper applications of CRISPR technology in many finely divided fields.

CRISPR-Cas9 Activities with Truncated 16-Nucleotide RNA Guides Are Tuned by Target Duplex Stability Beyond the RNA/DNA Hybrid

CRISPR-Cas9 has been adapted as a readily programmable genome manipulation agent, and continuing technological advances rely on an in-depth mechanistic understanding of Cas9 target discrimination. Cas9 interrogates a target by unwinding the DNA duplex to form an R-loop, where the RNA guide hybridizes with one of the DNA strands. It has been shown that RNA guides shorter than the normal length of 20-nucleotide (-nt) support Cas9 cleavage activity by enabling partial unwinding beyond the RNA/DNA hybrid. To investigate whether DNA segment beyond the RNA/DNA hybrid can impact Cas9 target discrimination with truncated guides, Cas9 double-stranded DNA cleavage rates (kcat) were measured with 16-nt guides on targets with varying sequences at +17 to +20 positions distal to the protospacer-adjacent-motif (PAM). The data reveal a log-linear inverse correlation between kcat and the PAM+(17-20) DNA duplex dissociation free energy (ΔGNN(17-20)0), with sequences having smaller ΔGNN(17-20)0 showing faster cleavage and a higher degree of unwinding. The results indicate that, with a 16-nt guide, "peripheral" DNA sequences beyond the RNA/DNA hybrid contribute to target discrimination by tuning the cleavage reaction transition state through the modulation of PAM-distal unwinding. The finding provides mechanistic insights for the further development of strategies that use RNA guide truncation to enhance Cas9 specificity.

Defining an optimal control for RNAi experiments with adult Schistosoma mansoni

In parasites such as Schistosoma mansoni, gene knockdown by RNA interference (RNAi) has become an indispensable tool for functional gene characterization. To distinguish target-specific RNAi effects versus off-target effects, controls are essential. To date, however, there is still no general agreement about suitable RNAi controls, which limits the comparability between studies. To address this point, we investigated three selected dsRNAs for their suitability as RNAi controls in experiments with adult S. mansoni in vitro. Two dsRNAs were of bacterial origin, the neomycin resistance gene (neoR) and the ampicillin resistance gene (ampR). The third one, the green fluorescent protein gene (gfp), originated from jellyfish. Following dsRNA application, we analyzed physiological parameters like pairing stability, motility, and egg production as well as morphological integrity. Furthermore, using RT-qPCR we evaluated the potential of the used dsRNAs to influence transcript patterns of off-target genes, which had been predicted by si-Fi (siRNA-Finder). At the physiological and morphological levels, we observed no obvious changes in the dsRNA treatment groups compared to an untreated control. However, we detected remarkable differences at the transcript level of gene expression. Amongst the three tested candidates, we suggest dsRNA of the E. coli ampR gene as the most suitable RNAi control.

Effect of Cas9 Protein on the Seed-Target Base Pair of the sgRNA/DNA Hybrid Duplex

CRISPR (clustered regularly interspaced short palindromic repeats)/CRISPR-associated protein (Cas9) has been widely used for gene editing. Not all guide RNAs can cleave the DNA efficiently remains a major challenge to CRISPR/Cas9-mediated genome engineering. Therefore, understanding how the Cas9 complex successfully and efficiently identifies specific functional targets through base-pairing has great implications for such applications. The 10-nt seed sequence at the 3' end of the guide RNA is critical to target recognition and cleavage. Here, through stretching molecular dynamics simulation, we studied the thermodynamics and kinetics of the binding-dissociation process of the seed base and the target DNA base with the Cas9 protein. The results showed that in the presence of Cas9 protein, the enthalpy change and entropy change in binding-dissociation of the seed base with the target are smaller than those without the Cas9 protein. The reduction of entropy penalty upon association with the protein resulted from the pre-organization of the seed base in an A-form helix, and the reduction of enthalpy change was due to the electrostatic attraction of the positively charged channel with the negative target DNA. The binding barrier coming from the entropy loss and the dissociation barrier resulting from the destruction of the base pair in the presence of Cas9 protein were lower than those without protein, which indicates that the seed region is crucial for efficiently searching the correct target by accelerating the binding rate and dissociating fast from the wrong target.

Substrate Specificity and Kinetics of RNA Hydrolysis by SARS-CoV-2 NSP10/14 Exonuclease

Severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2), the virus that causes COVID-19, continues to evolve resistance to vaccines and existing antiviral therapies at an alarming rate, increasing the need for new direct-acting antiviral drugs. Despite significant advances in our fundamental understanding of the kinetics and mechanism of viral RNA replication, there are still open questions regarding how the proofreading exonuclease (NSP10/NSP14 complex) contributes to replication fidelity and resistance to nucleoside analogs. Through single turnover kinetic analysis, we show that the preferred substrate for the exonuclease is double-stranded RNA without any mismatches. Double-stranded RNA containing a 3'-terminal remdesivir was hydrolyzed at a rate similar to a correctly base-paired cognate nucleotide. Surprisingly, single-stranded RNA or duplex RNA containing a 3'-terminal mismatch was hydrolyzed at rates 125- and 45-fold slower, respectively, compared to the correctly base-paired double-stranded RNA. These results define the substrate specificity and rate of removal of remdesivir for the exonuclease and outline rigorous kinetic assays that could help in finding next-generation exonuclease inhibitors or nucleoside analogs that are able to evade excision. These results also raise important questions about the role of the polymerase/exonuclease complex in proofreading during viral replication. Addressing these questions through rigorous kinetic analysis will facilitate the search for desperately needed antiviral drugs to combat COVID-19.

Pattern Recognition of microRNA Expression in Body Fluids Using Nanopore Decoding at Subfemtomolar Concentrations

This paper describes a method for detecting microRNA (miRNA) expression patterns using the nanopore-based DNA computing technology. miRNAs have shown promise as markers for cancer diagnosis due to their cancer type specificity, and therefore simple strategies for miRNA pattern recognition are required. We propose a system for pattern recognition of five types of miRNAs overexpressed in bile duct cancer (BDC). The information of miRNAs from BDC is encoded in diagnostic DNAs (dgDNAs) and decoded electrically by nanopore analysis. With this system, we succeeded in the label-free detection of miRNA expression patterns from the plasma of BDC patients. Moreover, our dgDNA-miRNA complexes can be detected at subfemtomolar concentrations, which is a significant improvement compared to previously reported limits of detection (∼10^-12 M) for similar analytical platforms. Nanopore decoding of dgDNA-encoded information represents a promising tool for simple and early cancer diagnosis.

The stability and number of nucleating interactions determine DNA hybridization rates in the absence of secondary structure

The kinetics of DNA hybridization are fundamental to biological processes and DNA-based technologies. However, the precise physical mechanisms that determine why different DNA sequences hybridize at different rates are not well understood. Secondary structure is one predictable factor that influences hybridization rates but is not sufficient on its own to fully explain the observed sequence-dependent variance. In this context, we measured hybridization rates of 43 different DNA sequences that are not predicted to form secondary structure and present a parsimonious physically justified model to quantify our observations. Accounting only for the combinatorics of complementary nucleating interactions and their sequence-dependent stability, the model achieves good correlation with experiment with only two free parameters. Our results indicate that greater repetition of Watson-Crick pairs increases the number of initial states able to proceed to full hybridization, with the stability of those pairings dictating the likelihood of such progression, thus providing new insight into the physical factors underpinning DNA hybridization rates.

Applicability of the nearest-neighbour model for pseudoknot RNAs

The validity of the nearest-neighbour (NN) model was verified in an RNA pseudoknot (PK) structure. The thermodynamic parameters of the second hairpin stem (S2) region, which separates the PK from a hairpin structure, were monitored using CD and UV melting. Different PKs with identical NN base pairs in the S2 region exhibited similar thermodynamic parameters, highlighting the validity of the NN model in this RNA tertiary structure motif.

Concurrent Gene Insertion, Deletion, and Inversion during the Construction of a Novel Attenuated BoHV-1 Using CRISPR/Cas9 Genome Editing

Bovine herpesvirus type I (BoHV-1) is an important pathogen that causes respiratory disease in bovines. The disease is prevalent worldwide, causing huge economic losses to the cattle industry. Gene-deficient vaccines with immunological markers to distinguish them from wild-type infections have become a mainstream in vaccine research and development. In order to knock out the gE gene BoHV-1, we employed the CRISPR/Cas9 system. Interesting phenomena were observed at the single guide RNA (sgRNA) splicing site, including gene insertion, gene deletion, and the inversion of 5′ and 3′ ends of the sgRNA splicing site. In addition to the deletion of the gE gene, the US9 gene, and the non-coding regions of gE and US9, it was found that the US4 sequence, US6 sequence, and part of the US7 sequence were inserted into the EGFP sgRNA splicing site and the 3′ end of the EGFP sequence was deleted. Similar to the BoHV-1 parent, the BoHV-1 mutants induced high neutralizing antibodies titer levels in mice. In summary, we developed a series of recombinant gE-deletion BoHV-1 samples using the CRISPR/Cas9 gene editing system. The mutant viruses with EGFP⁺ or EGFP⁻ will lay the foundation for research on BoHV-1 and vaccine development in the future.

Calculation of Energy for RNA/RNA and DNA/RNA Duplex Formation by Molecular Dynamics Simulation

Abstract

The development of approaches for predictive calculation of hybridization properties of various nucleic acid (NA) derivatives is the basis for the rational design of the NA-based constructs. Modern advances in computer modeling methods provide the feasibility of these calculations. We have analyzed the possibility of calculating the energy of DNA/RNA and RNA/RNA duplex formation using representative sets of complexes (65 and 75 complexes, respectively). We used the classical molecular dynamics (MD) method, the MMPBSA or MMGBSA approaches to calculate the enthalpy (ΔH°) component, and the quasi-harmonic approximation (Q-Harm) or the normal mode analysis (NMA) methods to calculate the entropy (ΔS°) contribution to the Gibbs energy ($$\Delta G_{{37}}^{^\circ }$$ ) of the NA complex formation. We have found that the MMGBSA method in the analysis of the MD trajectory of only the NA duplex and the empirical linear approximation allow calculation of the enthalpy of formation of the DNA, RNA, and hybrid duplexes of various lengths and GC content with an accuracy of 8.6%. Within each type of complex, the combination of rather efficient MMGBSA and Q-Harm approaches being applied to the trajectory of only the bimolecular complex makes it possible to calculate the $$\Delta G_{{37}}^{^\circ }$$ of the duplex formation with an error value of 10%. The high accuracy of predictive calculation for different types of natural complexes (DNA/RNA, DNA/RNA, and RNA/RNA) indicates the possibility of extending the considered approach to analogs and derivatives of nucleic acids, which gives a fundamental opportunity in the future to perform rational design of new types of NA-targeted sequence-specific compounds.

Transcriptome screening followed by integrated physicochemical and structural analyses for investigating RNA-mediated berberine activity

Non-coding RNAs are regarded as promising targets for the discovery of innovative drugs due to their abundance in the genome and their involvement in many biological processes. Phytochemicals (PCs) are the primary source of ligand-based drugs due to their broad spectrum of biological activities. Since many PCs are heterocyclic and have chemical groups potentially involved in the interaction with nucleic acids, detailed interaction analysis between PCs and RNA is crucial to explore the effect of PCs on RNA functions. In this study, an integrated approach for investigating interactions between PCs and RNAs were demonstrated to verify the RNA-mediated PCs functions by using berberine (BRB) as a model PC. RNA screening of a transcriptome library followed by sequence refinement found minimal RNA motif consisting of a cytosine bulge with U-A and G-U neighbouring base pairs for interaction with BRB. NMR-based structure determination and physicochemical analyses using chemical analogues of BRB demonstrated the importance of electrostatic and stacking interactions for sequence selective interaction and RNA stabilization. The selective interaction with a relatively small RNA motif based on a chemical structure of a planer heterocyclic highlights the biological activities of various PCs mediated by the interactions with particular functional RNAs. In addition, the systematic and quantitative investigations demonstrated in this study could be useful for the development of therapeutic chemicals targeting functional RNAs, based on the PCs, in the future.

Effects of Phosphoryl Guanidine Modification of Phosphate Residues on the Structure and Hybridization of Oligodeoxyribonucleotides

Phosphoryl guanidine oligonucleotides (PGOs) are promising tools for biological research and development of biosensors and therapeutics. We performed structural and hybridization analyses of octa-, deca-, and dodecamers with all phosphate residues modified by 1,3-dimethylimidazolidine-2-imine moieties. Similarity of the B-form double helix between native and modified duplexes was noted. In PGO duplexes, we detected a decrease in the proportion of C2'-endo and an increased proportion of C1'-exo sugar conformations of the modified chain. Applicability of the two-state model to denaturation transition of all studied duplexes was proved for the first time. Sequence-dependent effects of this modification on hybridization properties were observed. The thermal stability of PGO complexes is almost native at 100 mM NaCl and slightly increases with decreasing ionic strength. An increase in water activity and dramatic changes in interaction with cations and in solvation of PGOs and their duplexes were noted, resulting in slight elevation of the melting temperature after an ionic-strength decrease from 1 M NaCl down to deionized water. Decreased binding of sodium ions and decreased water solvation were documented for PGOs and their duplexes. In contrast to DNA, the PGO duplex formation leads to a release of several cations. The water shell is significantly more disordered near PGOs and their complexes. Nevertheless, changes in solvation during the formation of native and PGO complexes are similar and indicate that it is possible to develop models for predictive calculations of the thermodynamic properties of phosphoryl guanidine oligomers. Our results may help devise an approach for the rational design of PGOs as novel improved molecular probes and tools for many modern methods involving oligonucleotides.

Optimized design of antisense oligomers for targeted rRNA depletion

RNA sequencing (RNA-seq) is extensively used to quantify gene expression transcriptome-wide. Although often paired with polyadenylate (poly(A)) selection to enrich for messenger RNA (mRNA), many applications require alternate approaches to counteract the high proportion of ribosomal RNA (rRNA) in total RNA. Recently, digestion using RNaseH and antisense DNA oligomers tiling target rRNAs has emerged as an alternative to commercial rRNA depletion kits. Here, we present a streamlined, more economical RNaseH-mediated rRNA depletion with substantially lower up-front costs, using shorter antisense oligos only sparsely tiled along the target RNA in a 5-min digestion reaction. We introduce a novel Web tool, Oligo-ASST, that simplifies oligo design to target regions with optimal thermodynamic properties, and additionally can generate compact, common oligo pools that simultaneously target divergent RNAs, e.g. across different species. We demonstrate the efficacy of these strategies by generating rRNA-depletion oligos for Xenopus laevis and for zebrafish, which expresses two distinct versions of rRNAs during embryogenesis. The resulting RNA-seq libraries reduce rRNA to <5% of aligned reads, on par with poly(A) selection, and also reveal expression of many non-adenylated RNA species. Oligo-ASST is freely available at https://mtleelab.pitt.edu/oligo to design antisense oligos for any taxon or to target any abundant RNA for depletion.

Stability prediction of canonical and non-canonical structures of nucleic acids in various molecular environments and cells

This review provides the biophysicochemical background and recent advances in stability prediction of canonical and non-canonical structures of nucleic acids in various molecular environments and cells.