Molecular mechanisms underlying the extreme mechanical anisotropy of the flaviviral exoribonuclease-resistant RNAs (xrRNAs)

Tick-borne flavivirus exoribonuclease-resistant RNAs contain a double loop structure

Viruses from the Flaviviridae family contain human relevant pathogens that generate subgenomic noncoding RNAs during infection using structured exoribonuclease resistant RNAs (xrRNAs). These xrRNAs block progression of host cell's 5' to 3' exoribonucleases. The structures of several xrRNAs from mosquito-borne and insect-specific flaviviruses reveal a conserved fold in which a ring-like motif encircles the 5' end of the xrRNA. However, the xrRNAs found in tick-borne and no known vector flaviviruses have distinct characteristics, and their 3-D fold was unsolved. Here, we verify the presence of xrRNAs in the encephalitis-causing tick-borne Powassan Virus. We characterize their secondary structure and obtain a mid-resolution map of one of these xrRNAs using cryo-EM, revealing a unique double-loop ring element. Integrating these results with covariation analysis, biochemical data, and existing high-resolution structural information yields a model in which the core of the fold matches the previously solved xrRNA fold, but the expanded double loop ring is remodeled upon encountering the exoribonuclease. These results are representative of a broad class of xrRNAs and reveal a conserved strategy of structure-based exoribonuclease resistance achieved through a unique topology across a viral family of importance to global health.

Identifying the interactions conferring functional mechanical rigidity on RNase-resistant RNA from Zika virus

Some viruses counter host-cell efforts to digest invading viral RNA by using special structures resistant to host RNases, known as exoribonuclease-resistant RNAs (xrRNAs). xrRNAs typically form an unusual fold with the 5'-end threaded through a ring consisting of a multihelix junction closed by a pseudoknot. By using single-molecule force spectroscopy (SMFS), we previously showed that a Zika virus xrRNA is extremely rigid mechanically, withstanding very high forces, and that this mechanical resistance-not simply the knot-like fold topology-is essential for RNase resistance. Here, we have determined which interactions are most important for generating mechanical rigidity in the Zika virus xrRNA, by systematically mutating tertiary contacts. We found that removing any of the tertiary contacts involving the threaded 5' end was sufficient to abrogate mechanical resistance. In contrast, breaking a single pseudoknot base pair was not sufficient to do so: Two broken pairs were needed. This hierarchy of interaction importance for mechanical rigidity was supported by simulations mapping how mechanical tension was distributed within the xrRNA. For all mutants, RNase resistance varied in lock-step with mechanical resistance, confirming the primary role of mechanical rigidity in xrRNA function. This work reveals which interactions are most important for Zika xrRNA function, with implications for targeting the xrRNA therapeutically.

Lifetime of ground conformational state determines the activity of structured RNA

Biomolecules continually sample alternative conformations. Consequently, even the most energetically favored ground conformational state has a finite lifetime. Here, we show that, in addition to the three-dimensional (3D) structure, the lifetime of a ground conformational state determines its biological activity. Using hydrogen-deuterium exchange nuclear magnetic resonance spectroscopy, we found that Zika virus exoribonuclease-resistant RNA (xrRNA) encodes a ground conformational state with a lifetime that is ~105-107 longer than that of canonical base pairs. Mutations that shorten the apparent lifetime of the ground state without affecting its 3D structure decreased exoribonuclease resistance in vitro and impaired virus replication in cells. Additionally, we observed this exceptionally long-lived ground state in xrRNAs from diverse infectious mosquito-borne flaviviruses. These results demonstrate the biological importance of the lifetime of a preorganized ground state and further suggest that elucidating the lifetimes of dominant 3D structures of biomolecules may be crucial for understanding their behaviors and functions.

Identification of phoxim and omethoate using α-hemolysin nanopore and aptamers

Contamination with pesticides has inflicted substantial harm on human health; therefore, developing rapid, ultra-sensitive, and non-labelling simultaneous detection methods for multiple pesticides is necessary. In this study, we demonstrated that α-hemolysin (α-HL) nanopore sensor can detect and discriminate organophosphorus pesticides of phoxim and omethoate in a single nanopore without requiring labels of the probes or purification of the pesticides in real samples. Aptamers specifically recognise and bind pesticides to obtain pesticide-aptamer complexes that produce characteristic current signals while passing through the nanopore. Phoxim and omethoate were accurately distinguished by a portable instrument within minutes, and their detection sensitivity was up to the femtomole level, whether detected alone or simultaneously. The detection limits of phoxim and omethoate were 8.13 × 10-16 M and 4.16 × 10-15 M. The recoveries of phoxim and omethoate from pear, tomato, and cucumber samples were 82.0-107.0 % and 81.9-118.3 % respectively, with coefficient of variable below 8.0 %.

The pseudoknot structure of a viral RNA reveals a conserved mechanism for programmed exoribonuclease resistance

Exoribonuclease-resistant RNAs (xrRNAs) are viral RNA structures that block degradation by cellular 5'-3' exoribonucleases to produce subgenomic viral RNAs during infection. Initially discovered in flaviviruses, xrRNAs have since been identified in wide range of RNA viruses, including those that infect plants. High sequence variability among viral xrRNAs raises questions about the shared molecular features that characterize this functional RNA class. Here, we present the first structure of a plant-virus xrRNA in its active exoribonuclease-resistant conformation. The xrRNA forms a 9 base pair pseudoknot that creates a knot-like topology similar to that of flavivirus xrRNAs, despite lacking sequence similarity. Biophysical assays confirm a compact pseudoknot structure in solution, and functional studies validate its relevance both in vitro and during infection. Our study reveals how viral RNAs achieve a common functional outcome through highly divergent sequences and identifies the knot-like topology as a defining feature of xrRNAs.

Magnesium Regulates RNA Ring Dynamics and Folding in Subgenomic Flaviviral RNA

Mosquito-borne flaviviruses including dengue, Zika, yellow fever, and regional encephalitis produce a large amount of short subgenomic flaviviral RNAs during infection. A segment of these RNAs named as xrRNA1 features a multi-pseudoknot (PK)-associated structure, which resists the host cell enzyme (XRN1) from degrading the viral RNA. We investigate how this long-range RNA PK folds in the presence of counterions, specifically in a mix of monovalent (K+) and divalent (Mg2+) salts at physiological concentrations. In this study, we use extensive explicit solvent molecular dynamics (MD) simulations to characterize the RNA ion environment of the folded RNA conformation, as determined by the crystal structure. This allowed us to identify the precise locations of various coordinated RNA-Mg2+ interactions, including inner-sphere/chelated and outer-sphere coordinated Mg2+. Given that RNA folding involves large-scale conformational changes, making it challenging to explore through classical MD simulations, we investigate the folding mechanism of xrRNA1 using an all-atom structure-based RNA model with a hybrid implicit-explicit treatment of the ion environment via the dynamic counterion condensation model, both with and without physiological Mg2+ concentration. The study reveals potential folding pathways for this xrRNA1, which is consistent with the results obtained from optical tweezer experiments. The equilibrium and free energy simulations both capture a dynamic equilibrium between the ring-open and ring-close states of the RNA, driven by a long-range PK interaction. Free energy calculations reveal that with the addition of Mg2+ ions, the equilibrium shifts more toward the ring-close state. A detailed analysis of the free energy pathways and ion-mediated contact probability map highlights the critical role of Mg2+ in bridging G50 and A33. This Mg2+-mediated connection helps form the long-range PK which in turn controls the transition between the ring-open and ring-close states. The study underscores the critical role of Mg2+ in the RNA folding transition, highlighting specific locations of Mg2+ contributing to the stabilization of long-range PK connections likely to enhance the robustness of Xrn1 resistance of flaviviral xrRNAs.

Structure-based investigation of a DNA aptamer targeting PTK7 reveals an intricate 3D fold guiding functional optimization

DNA aptamers have emerged as novel molecular tools in disease theranostics owing to their high binding affinity and specificity for protein targets, which rely on their ability to fold into distinctive three-dimensional (3D) structures. However, delicate atomic interactions that shape the 3D structures are often ignored when designing and modeling aptamers, leading to inefficient functional optimization. Challenges persist in determining high-resolution aptamer-protein complex structures. Moreover, the experimentally determined 3D structures of DNA molecules with exquisite functions remain scarce. These factors impede our comprehension and optimization of some important DNA aptamers. Here, we performed a streamlined solution NMR-based structural investigation on the 41-nt sgc8c, a prominent DNA aptamer used to target membrane protein tyrosine kinase 7, for cancer theranostics. We show that sgc8c prefolds into an intricate three-way junction (3WJ) structure stabilized by long-range tertiary interactions and extensive base-base stackings. Delineated by NMR chemical shift perturbations, site-directed mutagenesis, and 3D structural information, we identified essential nucleotides constituting the key functional elements of sgc8c that are centralized at the core of 3WJ. Leveraging the well-established structure-function relationship, we efficiently engineered two sgc8c variants by modifying the apical loop and introducing L-DNA base pairs to simultaneously enhance thermostability, biostability, and binding affinity for both protein and cell targets, a feat not previously attained despite extensive efforts. This work showcases a simplified NMR-based approach to comprehend and optimize sgc8c without acquiring the complex structure, and offers principles for the sophisticated structure-function organization of DNA molecules.

Exploring the Molecular Mechanism of Niuxi-Mugua Formula in Treating Coronavirus Disease 2019 via Network Pharmacology, Computational Biology, and Surface Plasmon Resonance Verification

BACKGROUND

In China, Niuxi-Mugua formula (NMF) has been widely used to prevent and treat coronavirus disease 2019 (COVID-19). However, the mechanism of NMF for treating COVID-19 is not yet fully understood.

OBJECTIVE

This study aimed to explore the potential mechanism of NMF for treating COVID- 19 by network pharmacology, computational biology, and surface plasmon resonance (SPR) verification.

MATERIALS AND METHODS

The NMF-compound-target network was constructed to screen the key compounds, and the Molecular Complex Detection (MCODE) tool was used to screen the preliminary key genes. The overlapped genes (OGEs) and the preliminary key genes were further analyzed by enrichment analysis. Then, the correlation analysis of immune signatures and the preliminary key genes was performed. Molecular docking and molecular dynamic (MD) simulation assays were applied to clarify the interactions between key compounds and key genes. Moreover, the SPR interaction experiment was used for further affinity kinetic verification.

RESULTS

Lipid and atherosclerosis, TNF, IL-17, and NF-kappa B signaling pathways were the main pathways of NMF in the treatment of COVID-19. There was a positive correlation between almost the majority of immune signatures and all preliminary key genes. The key compounds and the key genes were screened out, and they were involved in the main pathways of NMF for treating COVID-19. Moreover, the binding affinities of most key compounds binding to key genes were good, and IL1B-Quercetin had the best binding stability. SPR analysis further demonstrated that IL1B-Quercetin showed good binding affinity.

CONCLUSION

Our findings provided theoretical grounds for NMF in the treatment of COVID-19.

Real-time label-free detection of dynamic aptamer-small molecule interactions using a nanopore nucleic acid conformational sensor

Nucleic acids can undergo conformational changes upon binding small molecules. These conformational changes can be exploited to develop new therapeutic strategies through control of gene expression or triggering of cellular responses and can also be used to develop sensors for small molecules such as neurotransmitters. Many analytical approaches can detect dynamic conformational change of nucleic acids, but they need labeling, are expensive, and have limited time resolution. The nanopore approach can provide a conformational snapshot for each nucleic acid molecule detected, but has not been reported to detect dynamic nucleic acid conformational change in response to small -molecule binding. Here we demonstrate a modular, label-free, nucleic acid-docked nanopore capable of revealing time-resolved, small molecule-induced, single nucleic acid molecule conformational transitions with millisecond resolution. By using the dopamine-, serotonin-, and theophylline-binding aptamers as testbeds, we found that these nucleic acids scaffolds can be noncovalently docked inside the MspA protein pore by a cluster of site-specific charged residues. This docking mechanism enables the ion current through the pore to characteristically vary as the aptamer undergoes conformational changes, resulting in a sequence of current fluctuations that report binding and release of single ligand molecules from the aptamer. This nanopore tool can quantify specific ligands such as neurotransmitters, elucidate nucleic acid-ligand interactions, and pinpoint the nucleic acid motifs for ligand binding, showing the potential for small molecule biosensing, drug discovery assayed via RNA and DNA conformational changes, and the design of artificial riboswitch effectors in synthetic biology.

Quantitative characterization of cell physiological state based on dynamical cell mechanics for drug efficacy indication

Cell mechanics is essential to cell development and function, and its dynamics evolution reflects the physiological state of cells. Here, we investigate the dynamical mechanical properties of single cells under various drug conditions, and present two mathematical approaches to quantitatively characterizing the cell physiological state. It is demonstrated that the cellular mechanical properties upon the drug action increase over time and tend to saturate, and can be mathematically characterized by a linear time-invariant dynamical model. It is shown that the transition matrices of dynamical cell systems significantly improve the classification accuracies of the cells under different drug actions. Furthermore, it is revealed that there exists a positive linear correlation between the cytoskeleton density and the cellular mechanical properties, and the physiological state of a cell in terms of its cytoskeleton density can be predicted from its mechanical properties by a linear regression model. This study builds a relationship between the cellular mechanical properties and the cellular physiological state, adding information for evaluating drug efficacy.

3D Structural Analysis of Long Noncoding RNA by Small Angle X-ray Scattering and Computational Modeling

Small angle X-ray scattering (SAXS) has been widely applied as an enabling integrative technique for comprehensive analysis of the structure of biomacromolecules by multiple, complementary techniques in solution. SAXS in combination with computational modeling can be a powerful strategy bridging the secondary and 3D structural analysis of large RNAs, including the long noncoding RNAs (lncRNA). Here, we outline the major procedures and techniques in the combined use of SAXS and computational modeling for 3D structural characterization of a lncRNA, the subgenomic flaviviral RNA from Zika virus. lncRNA production and purification, RNA buffer and sample preparation for SAXS experiments, SAXS data collection and analysis, SAXS-aided RNA 3D structure prediction, and computational modeling are described.

Enhancement of prime editing via xrRNA motif-joined pegRNA

The prime editors (PEs) have shown great promise for precise genome modification. However, their suboptimal efficiencies present a significant technical challenge. Here, by appending a viral exoribonuclease-resistant RNA motif (xrRNA) to the 3'-extended portion of pegRNAs for their increased resistance against degradation, we develop an upgraded PE platform (xrPE) with substantially enhanced editing efficiencies in multiple cell lines. A pan-target average enhancement of up to 3.1-, 4.5- and 2.5-fold in given cell types is observed for base conversions, small deletions, and small insertions, respectively. Additionally, xrPE exhibits comparable edit:indel ratios and similarly minimal off-target editing as the canonical PE3. Of note, parallel comparison of xrPE to the most recently developed epegRNA-based PE system shows their largely equivalent editing performances. Our study establishes a highly adaptable platform of improved PE that shall have broad implications.

Pseudoknot length modulates the folding, conformational dynamics, and robustness of Xrn1 resistance of flaviviral xrRNAs

To understand how RNA dynamics is regulated and connected to its function, we investigate the folding, conformational dynamics and robustness of Xrn1 resistance of a set of flaviviral xrRNAs using SAXS, smFRET and in vitro enzymatic assays. Flaviviral xrRNAs form discrete ring-like 3D structures, in which the length of a conserved long-range pseudoknot (PK2) ranges from 2 bp to 7 bp. We find that xrRNAs' folding, conformational dynamics and Xrn1 resistance are strongly correlated and highly Mg²⁺-dependent, furthermore, the Mg²⁺-dependence is modulated by PK2 length variations. xrRNAs with long PK2 require less Mg²⁺ to stabilize their folding, exhibit reduced conformational dynamics and strong Xrn1 resistance even at low Mg²⁺, and tolerate mutations at key tertiary motifs at high Mg²⁺, which generally are destructive to xrRNAs with short PK2. These results demonstrate an unusual regulatory mechanism of RNA dynamics providing insights into the functions and future biomedical applications of xrRNAs.

Mechanical strength of RNA knot in Zika virus protects against cellular defenses

Unusual knot-like structures recently discovered in viral exoribonuclease-resistant RNAs (xrRNAs) prevent digestion by host RNases to create subgenomic RNAs enhancing infection and pathogenicity. xrRNAs are proposed to prevent digestion through mechanical resistance to unfolding. However, their unfolding force has not been measured, and the factors determining RNase resistance are unclear. Furthermore, how these knots fold remains unknown. Unfolding a Zika virus xrRNA with optical tweezers revealed that it was the most mechanically stable RNA yet observed. The knot formed by threading the 5' end into a three-helix junction before pseudoknot interactions closed a ring around it. The pseudoknot and tertiary contacts stabilizing the threaded 5' end were both required to generate extreme force resistance, whereas removing a 5'-end contact produced a low-force knot lacking RNase resistance. These results indicate mechanical resistance plays a central functional role, with the fraction of molecules forming extremely high-force knots determining the RNase resistance level.