RLDOCK: A New Method for Predicting RNA-Ligand Interactions

Single-step discovery of high-affinity RNA ligands by UltraSelex

Aptamers, nucleic acid ligands targeting specific molecules, have emerged as drug candidates, sensors, imaging tools and nanotechnology building blocks. The predominant method for their discovery, systematic evolution of ligands by exponential enrichment, while successful, is laborious, time-consuming and often results in candidates enriched for unintended criteria. Here we present UltraSelex, a noniterative method that combines biochemical partitioning, high-throughput sequencing and computational signal-to-background rank modeling for discovering RNA aptamers in about 1 day. UltraSelex identified high-affinity RNA aptamers capable of binding a fluorogenic silicon rhodamine dye and two protein targets, the SARS-CoV-2 RNA-dependent RNA polymerase and HIV reverse transcriptase, enabling live-cell RNA imaging and efficient enzyme inhibition, respectively. From the ranked sequences, minimal aptamer motifs could be easily inferred. UltraSelex provides a rapid route to reveal new drug candidates and diagnostic tools.

RNAmigos2: accelerated structure-based RNA virtual screening with deep graph learning

RNAs are a vast reservoir of untapped drug targets. Structure-based virtual screening (VS) identifies candidate molecules by leveraging binding site information, traditionally using molecular docking simulations. However, docking struggles to scale with large compound libraries and RNA targets. Machine learning offers a solution but remains underdeveloped for RNA due to limited data and practical evaluations. We introduce a data-driven VS pipeline tailored for RNA, utilizing coarse-grained 3D modeling, synthetic data augmentation, and RNA-specific self-supervision. Our model achieves a 10,000x speedup over docking while ranking active compounds in the top 2.8% on structurally distinct test sets. It is robust to binding site variations and successfully screens unseen RNA riboswitches in a 20,000-compound in-vitro microarray, with a mean enrichment factor of 2.93 at 1%. This marks the first experimentally validated success of structure-based deep learning for RNA VS.

Optimization of Structure-Guided Development of Chemical Probes for the Pseudoknot RNA of the Frameshift Element in SARS-CoV-2

Targeting the RNA genome of SARS-CoV-2 is a viable option for antiviral drug development. We explored three ligand binding sites of the core pseudoknot RNA of the SARS-CoV-2 frameshift element. We iteratively optimized ligands, based on improved affinities, targeting these binding sites and report on structural and dynamic properties of the three identified binding sites. Available experimental 3D structures of the pseudoknot element were compared to SAXS and NMR data to validate its dominant folding state in solution. In order to experimentally map in silico predicted binding sites, NMR assignments of the majority of nucleobases were achieved by segmental labeling of the pseudoknot RNA and isotope-filtered NMR experiments at 1.2 GHz, demonstrating the value of NMR spectroscopy to supplement modelling and docking data. Optimized ligands with enhanced affinity were shown to specifically inhibit frameshifting without affecting 0-frame translation in cell-free translation assays, establishing the frameshift element as target for drug-like ligands of low molecular weight.

Comprehensive Evaluation of End-Point Free Energy Methods in DNA-Ligand Interaction Predictions

Deoxyribonucleic acid (DNA) serves as a repository of genetic information in cells and is a critical molecular target for various antibiotics and anticancer drugs. A profound understanding of small molecule interaction with DNA is crucial for the rational design of DNA-targeted therapies. While the molecular mechanics/Poisson-Boltzmann surface area (MM/PBSA) and molecular mechanics/generalized Born surface area (MM/GBSA) approaches have been well established for predicting protein-ligand binding, their application to DNA-ligand interactions has been less explored. In this study, we systematically investigated the binding of 13 diverse small molecules to DNA, evaluating the accuracy of DNA-ligand interaction predictions across different solvation approaches, interior dielectric constants (εin), and molecular force fields. Our results demonstrate that MM/PBSA, using energy-minimized structures (the bsc1 force field and εin = 20), provides the best correlation (Rp = -0.742) with experimental binding affinities, surpassing the performance of rDock scoring functions (best Rp = -0.481). Notably, the interior dielectric constant was found to significantly impact DNA-ligand binding free energy predictions, especially for MM/PBSA. Moreover, both MM/PBSA and MM/GBSA predictions (εin = 16 or 20) exhibited superior performance in distinguishing native-like binding modes within the top-10 poses from decoys, compared to the molecular docking tools used in this study. However, the popular docking software PLANTS demonstrates notable efficacy in predicting the top-1 binding pose. Given the considerably higher computational cost of MM/PBSA, MM/GBSA rescoring with higher εin = 16 or 20 is more efficient for recognizing the native-like binding poses for DNA-ligand systems. This study presents the first detailed exploration of end-point free energy calculations in the context of DNA-ligand interactions and offers valuable insights for the application of the MM/PB(GB)SA methods in this domain.

XDock: A General Docking Method for Modeling Protein-Ligand and Nucleic Acid-Ligand Interactions

Molecular docking is an essential computational tool in structure-based drug discovery and the investigation of the molecular mechanisms underlying biological processes. Despite the development of many molecular docking programs for various systems, a universal tool that can accurately dock ligands across multiple system types remains elusive. Meeting the need, we developed XDock, a versatile docking framework built for both protein-ligand and nucleic acid-ligand interactions. XDock efficiently accounts for ligand flexibility by docking multiple conformations of a ligand and flexibly refining the final binding poses. It utilizes a distance geometric method for ligand sampling and leverages our knowledge-based scoring functions for assessing protein-ligand and nucleic acid-ligand interactions. XDock has undergone extensive validations on diverse benchmarks of protein-ligand and nucleic acid-ligand complexes and was compared with six other docking methods, including DOCK 6, AutoDock Vina, PLANTS, LeDock, rDock, and RLDock. In addition, XDock is also computationally efficient and on average can dock a ligand within 1 min.

Harnessing Computational Approaches for RNA-Targeted Drug Discovery

RNA molecules have emerged as promising therapeutic targets due to their diverse functional and regulatory roles within cells. Computational modeling in RNA-targeted drug discovery presents a significant opportunity to expedite the discovery of novel small molecule compounds. However, this field encounters unique challenges compared to protein-targeted drug design, primarily due to limited experimental data availability and current models' inability to adequately address RNA's conformational flexibility during ligand recognition. Despite these challenges, several studies have successfully identified active RNA-targeting compounds using structure-based approaches or quantitative structure-activity relationship (QSAR) models. This review offers an overview of recent advancements in modeling RNA-small molecule interactions, emphasizing practical applications of computational methods in RNA-targeted drug discovery. Additionally, we survey existing databases that catalog nucleic acid-small molecule interactions. As interest in RNA-small molecule interactions grows and curated databases expand, the field anticipates rapid development. Novel computational models are poised to enhance the identification of potent and selective small-molecule modulators for therapeutic needs.

Rationally Designed G-Quadruplex Selective "Turn-On" NIR Fluorescent Probe with Large Stokes Shift for Nucleic Acid Research-Based Applications

Guanine-rich DNA/RNA sequences can form Hoogsteen bonds to adopt noncanonical secondary structures called G-quadruplexes, and these have been associated with diverse cellular processes. There has been considerable research interest in the design of G4-interacting ligands for cellular probing of the G4 structure and understanding its associated biological function. Most of the fluorescent G4 ligands either do not have significant selectivity over other nucleic acid structures, have high Stokes shift, or are not in the near-infrared (NIR) region, which limits its cellular visualization. The current work involves the rational design and synthesis of NIR fluorescent probes comprising a (Z)-1-methyl-2-((3-methylbenzo[d]thiazol-2(3H)-ylidene)methyl)quinolin-1-ium scaffold. Among the designed molecules, 4a exhibited far-red fluorescence (λmax = 680 nm) with large Stokes shift (∼182 nm) upon selective binding to human telomeric G-quadruplexes. The dye 4a does not disturb the conformation and stability of G-quadruplexes, thereby making it suitable for nucleic acid research based applications. Interestingly, 4a showed remarkable selectivity over single- and double-stranded structures in contrast to a commercially available quadruplex binding probe, Thiazole orange (TO). The molecular docking studies indicate that 4a binds at the groove region of the telomeric DNA G-quadruplex through π-π stacking interactions with the quinoline and amine-substituted phenyl ring and with the phosphate backbone through anion-π interactions with the benzothiazole ring. The designed molecule 4a has interesting photophysical properties, cell permeability, and biocompatibility with minimal cytotoxicity. Fluorescence imaging studies in live HeLa cells showed that probe 4a binds to the transient population of the DNA G-quadruplex in the nucleus and RNA quadruplexes in the cytoplasm. In brief, G-quadruplex NIR fluorescent probe 4a with a higher signal/noise ratio has significant potential for cellular imaging studies and thus opens avenues to decipher the biological pathways for better understanding of G-quadruplex biology.

A "pseudo" fluorescence indicator displacement assay based on "light-up" sanguinarine for identifying small ligands targeting HIV-1 RRE RNA

Developing effective methods to identify drugs that can target HIV-1 Rev response element (RRE) RNA and block the interaction between Rev and RRE has practical significance in the treatment of AIDS. Fluorescence indicator displacement (FID) assay was commonly employed to identify small ligands binding to RNA. In this study, the non-fluorescent sanguinarine (Sang) was used as a novel "pseudo" fluorescence indicator to identify small ligands targeting RRE through its fluorescence "light-up", assisted by β-cyclodextrin (β-CD) or nano-SiO2. The fluorescence enhancement effect of β-CD or nano-SiO2 on Sang was initially examined. Subsequently, the Sang-RRE interaction was investigated and validated using circular dichroism, and the binding constants and binding sites of Sang to RRE were obtained; the competition between Sang and Rev was explored and confirmed through molecular docking. Finally, the experimental parameters of the "pseudo" FID assay including reaction time, along with β-CD and nano-SiO2 concentrations, were optimized using a well-known Rev inhibitor (neomycin). Under the optimized experimental conditions, the evaluation of three positive inhibitors and one negative control was performed. These findings presented a novel approach for identifying small ligands targeting HIV-1 RRE RNA using a non-fluorescent ligand in the FID assay and established a promising screening platform.

SPRank─A Knowledge-Based Scoring Function for RNA-Ligand Pose Prediction and Virtual Screening

The growing interest in RNA-targeted drugs underscores the need for computational modeling of interactions between RNA molecules and small compounds. Having a reliable scoring function for RNA-ligand interactions is essential for effective computational drug screening. An ideal scoring function should not only predict the native pose for ligand binding but also rank the affinity of the binding for different ligands. However, existing scoring functions are primarily designed to predict the native binding modes for a given RNA-ligand pair and have not been thoroughly assessed for virtual screening purposes. In this paper, we introduce SPRank, a combination of machine-learning and knowledge-based scoring functions developed through a weighted iterative approach, specifically designed to tackle both binding mode prediction and virtual screening challenges. Our approach incorporates third-party docking software, such as rDock and AutoDock Vina, to sample flexible ligands against an ensemble of RNA structures, capturing the conformational flexibility of both the RNA and the ligand. Through rigorous testing, SPRank demonstrates improved performance compared to the tested scoring functions across four test sets comprising 122, 42, 55, and 71 nucleic acid-ligand complexes. Furthermore, SPRank exhibits improved performance in virtual screening tests targeting the HIV-1 TAR ensemble, which highlights its advantage in drug discovery. These results underscore the advantages of SPRank as a potentially promising tool for the RNA-targeted drug design. The source code of SPRank and the data sets are freely accessible at https://github.com/Vfold-RNA/SPRank.

Can Current Molecular Docking Methods Accurately Predict RNA Inhibitors?

Ribonucleic acids (RNAs), particularly the noncoding RNAs, play key roles in cancer, making them attractive drug targets. While conventional methods such as high throughput screening are resource-intensive, computational methods such as RNA-ligand docking can be used as an alternative. However, currently available docking methods are fine-tuned to perform protein-ligand and protein-protein docking. In this work, we evaluated three commonly used docking methods─AutoDock Vina, HADDOCK, and HDOCK─alongside RLDOCK, which is specifically designed for RNA-ligand docking. Our evaluation was based on several criteria including cognate docking, blind docking, scoring potential, and ranking potential. In cognate docking, only RLDOCK showed a success rate of 70% for the top-scoring docked pose. Despite this, all four docking methods did not achieve an overall success rate exceeding 50% amidst our attempt to refine the top-scoring docked poses using molecular dynamics simulations. Meanwhile, all four docking methods showed poor performance in scoring potential evaluation. Although AutoDock Vina achieved an area under the receiver operating characteristic curve of 0.70, it showed poor performance in terms of Matthews' correlation coefficient, precision, enrichment factors, and normalized enrichment factors at 1, 2, and 5%. These results highlight the growing need for further optimization of docking methods to assess RNA-ligand interactions.

Cracking the Code: Enhancing Molecular Tools for Progress in Nanobiotechnology

Nature continually refines its processes for optimal efficiency, especially within biological systems. This article explores the collaborative efforts of researchers worldwide, aiming to mimic nature's efficiency by developing smarter and more effective nanoscale technologies and biomaterials. Recent advancements highlight progress and prospects in leveraging engineered nucleic acids and proteins for specific tasks, drawing inspiration from natural functions. The focus is developing improved methods for characterizing, understanding, and reprogramming these materials to perform user-defined functions, including personalized therapeutics, targeted drug delivery approaches, engineered scaffolds, and reconfigurable nanodevices. Contributions from academia, government agencies, biotech, and medical settings offer diverse perspectives, promising a comprehensive approach to broad nanobiotechnology objectives. Encompassing topics from mRNA vaccine design to programmable protein-based nanocomputing agents, this work provides insightful perspectives on the trajectory of nanobiotechnology toward a future of enhanced biomimicry and technological innovation.

Advances in machine-learning approaches to RNA-targeted drug design

RNA molecules play multifaceted functional and regulatory roles within cells and have garnered significant attention in recent years as promising therapeutic targets. With remarkable successes achieved by artificial intelligence (AI) in different fields such as computer vision and natural language processing, there is a growing imperative to harness AI's potential in computer-aided drug design (CADD) to discover novel drug compounds that target RNA. Although machine-learning (ML) approaches have been widely adopted in the discovery of small molecules targeting proteins, the application of ML approaches to model interactions between RNA and small molecule is still in its infancy. Compared to protein-targeted drug discovery, the major challenges in ML-based RNA-targeted drug discovery stem from the scarcity of available data resources. With the growing interest and the development of curated databases focusing on interactions between RNA and small molecule, the field anticipates a rapid growth and the opening of a new avenue for disease treatment. In this review, we aim to provide an overview of recent advancements in computationally modeling RNA-small molecule interactions within the context of RNA-targeted drug discovery, with a particular emphasis on methodologies employing ML techniques.

Assessing the performance of MM/PBSA and MM/GBSA methods. 10. Prediction reliability of binding affinities and binding poses for RNA-ligand complexes

Ribonucleic acid (RNA)-ligand interactions play a pivotal role in a wide spectrum of biological processes, ranging from protein biosynthesis to cellular reproduction. This recognition has prompted the broader acceptance of RNA as a viable candidate for drug targets. Delving into the atomic-scale understanding of RNA-ligand interactions holds paramount importance in unraveling intricate molecular mechanisms and further contributing to RNA-based drug discovery. Computational approaches, particularly molecular docking, offer an efficient way of predicting the interactions between RNA and small molecules. However, the accuracy and reliability of these predictions heavily depend on the performance of scoring functions (SFs). In contrast to the majority of SFs used in RNA-ligand docking, the end-point binding free energy calculation methods, such as molecular mechanics/generalized Born surface area (MM/GBSA) and molecular mechanics/Poisson Boltzmann surface area (MM/PBSA), stand as theoretically more rigorous approaches. Yet, the evaluation of their effectiveness in predicting both binding affinities and binding poses within RNA-ligand systems remains unexplored. This study first reported the performance of MM/PBSA and MM/GBSA with diverse solvation models, interior dielectric constants (εin) and force fields in the context of binding affinity prediction for 29 RNA-ligand complexes. MM/GBSA is based on short (5 ns) molecular dynamics (MD) simulations in an explicit solvent with the YIL force field; the GBGBn2 model with higher interior dielectric constant (εin = 12, 16 or 20) yields the best correlation (Rp = -0.513), which outperforms the best correlation (Rp = -0.317, rDock) offered by various docking programs. Then, the efficacy of MM/GBSA in identifying the near-native binding poses from the decoys was assessed based on 56 RNA-ligand complexes. However, it is evident that MM/GBSA has limitations in accurately predicting binding poses for RNA-ligand systems, particularly compared with notably proficient docking programs like rDock and PLANTS. The best top-1 success rate achieved by MM/GBSA rescoring is 39.3%, which falls below the best results given by docking programs (50%, PLNATS). This study represents the first evaluation of MM/PBSA and MM/GBSA for RNA-ligand systems and is expected to provide valuable insights into their successful application to RNA targets.

RmsdXNA: RMSD prediction of nucleic acid-ligand docking poses using machine-learning method

Small molecule drugs can be used to target nucleic acids (NA) to regulate biological processes. Computational modeling methods, such as molecular docking or scoring functions, are commonly employed to facilitate drug design. However, the accuracy of the scoring function in predicting the closest-to-native docking pose is often suboptimal. To overcome this problem, a machine learning model, RmsdXNA, was developed to predict the root-mean-square-deviation (RMSD) of ligand docking poses in NA complexes. The versatility of RmsdXNA has been demonstrated by its successful application to various complexes involving different types of NA receptors and ligands, including metal complexes and short peptides. The predicted RMSD by RmsdXNA was strongly correlated with the actual RMSD of the docked poses. RmsdXNA also outperformed the rDock scoring function in ranking and identifying closest-to-native docking poses across different structural groups and on the testing dataset. Using experimental validated results conducted on polyadenylated nuclear element for nuclear expression triplex, RmsdXNA demonstrated better screening power for the RNA-small molecule complex compared to rDock. Molecular dynamics simulations were subsequently employed to validate the binding of top-scoring ligand candidates selected by RmsdXNA and rDock on MALAT1. The results showed that RmsdXNA has a higher success rate in identifying promising ligands that can bind well to the receptor. The development of an accurate docking score for a NA-ligand complex can aid in drug discovery and development advancements. The code to use RmsdXNA is available at the GitHub repository https://github.com/laiheng001/RmsdXNA.

Frontiers and Challenges of Computing ncRNAs Biogenesis, Function and Modulation

Non-coding RNAs (ncRNAs), generated from nonprotein coding DNA sequences, constitute 98-99% of the human genome. Non-coding RNAs encompass diverse functional classes, including microRNAs, small interfering RNAs, PIWI-interacting RNAs, small nuclear RNAs, small nucleolar RNAs, and long non-coding RNAs. With critical involvement in gene expression and regulation across various biological and physiopathological contexts, such as neuronal disorders, immune responses, cardiovascular diseases, and cancer, non-coding RNAs are emerging as disease biomarkers and therapeutic targets. In this review, after providing an overview of non-coding RNAs' role in cell homeostasis, we illustrate the potential and the challenges of state-of-the-art computational methods exploited to study non-coding RNAs biogenesis, function, and modulation. This can be done by directly targeting them with small molecules or by altering their expression by targeting the cellular engines underlying their biosynthesis. Drawing from applications, also taken from our work, we showcase the significance and role of computer simulations in uncovering fundamental facets of ncRNA mechanisms and modulation. This information may set the basis to advance gene modulation tools and therapeutic strategies to address unmet medical needs.

Comparative Assessment of Pose Prediction Accuracy in RNA-Ligand Docking

Structure-based virtual high-throughput screening is used in early-stage drug discovery. Over the years, docking protocols and scoring functions for protein-ligand complexes have evolved to improve the accuracy in the computation of binding strengths and poses. In the past decade, RNA has also emerged as a target class for new small-molecule drugs. However, most ligand docking programs have been validated and tested for proteins and not RNA. Here, we test the docking power (pose prediction accuracy) of three state-of-the-art docking protocols on 173 RNA-small molecule crystal structures. The programs are AutoDock4 (AD4) and AutoDock Vina (Vina), which were designed for protein targets, and rDock, which was designed for both protein and nucleic acid targets. AD4 performed relatively poorly. For RNA targets for which a crystal structure of a bound ligand used to limit the docking search space is available and for which the goal is to identify new molecules for the same pocket, rDock performs slightly better than Vina, with success rates of 48% and 63%, respectively. However, in the more common type of early-stage drug discovery setting, in which no structure of a ligand-target complex is known and for which a larger search space is defined, rDock performed similarly to Vina, with a low success rate of ∼27%. Vina was found to have bias for ligands with certain physicochemical properties, whereas rDock performs similarly for all ligand properties. Thus, for projects where no ligand-protein structure already exists, Vina and rDock are both applicable. However, the relatively poor performance of all methods relative to protein-target docking illustrates a need for further methods refinement.

How Good Are Current Docking Programs at Nucleic Acid-Ligand Docking? A Comprehensive Evaluation

Nucleic acid (NA)-ligand interactions are of paramount importance in a variety of biological processes, including cellular reproduction and protein biosynthesis, and therefore, NAs have been broadly recognized as potential drug targets. Understanding NA-ligand interactions at the atomic scale is essential for investigating the molecular mechanism and further assisting in NA-targeted drug discovery. Molecular docking is one of the predominant computational approaches for predicting the interactions between NAs and small molecules. Despite the availability of versatile docking programs, their performance profiles for NA-ligand complexes have not been thoroughly characterized. In this study, we first compiled the largest structure-based NA-ligand binding data set to date, containing 800 noncovalent NA-ligand complexes with clearly identified ligands. Based on this extensive data set, eight frequently used docking programs, including six protein-ligand docking programs (LeDock, Surflex-Dock, UCSF Dock6, AutoDock, AutoDock Vina, and PLANTS) and two specific NA-ligand docking programs (rDock and RLDOCK), were systematically evaluated in terms of binding pose and binding affinity predictions. The results demonstrated that some protein-ligand docking programs, specifically PLANTS and LeDock, produced more promising or comparable results compared with the specialized NA-ligand docking programs. Among the programs evaluated, PLANTS, rDock, and LeDock showed the highest performance in binding pose prediction, and their top-1 and best root-mean-square deviation (rmsd) success rates were as follows: PLANTS (35.93 and 76.05%), rDock (27.25 and 72.16%), and LeDock (27.40 and 64.37%). Compared with the moderate level of binding pose prediction, few programs were successful in binding affinity prediction, and the best correlation (Rp = -0.461) was observed with PLANTS. Finally, further comparison with the latest NA-ligand docking program (NLDock) on four well-established data sets revealed that PLANTS and LeDock outperformed NLDock in terms of binding pose prediction on all data sets, demonstrating their significant potential for NA-ligand docking. To the best of our knowledge, this study is the most comprehensive evaluation of popular molecular docking programs for NA-ligand systems.

RPflex: A Coarse-Grained Network Model for RNA Pocket Flexibility Study

RNA regulates various biological processes, such as gene regulation, RNA splicing, and intracellular signal transduction. RNA’s conformational dynamics play crucial roles in performing its diverse functions. Thus, it is essential to explore the flexibility characteristics of RNA, especially pocket flexibility. Here, we propose a computational approach, RPflex, to analyze pocket flexibility using the coarse-grained network model. We first clustered 3154 pockets into 297 groups by similarity calculation based on the coarse-grained lattice model. Then, we introduced the flexibility score to quantify the flexibility by global pocket features. The results show strong correlations between the flexibility scores and root-mean-square fluctuation (RMSF) values, with Pearson correlation coefficients of 0.60, 0.76, and 0.53 in Testing Sets I–III. Considering both flexibility score and network calculations, the Pearson correlation coefficient was increased to 0.71 in flexible pockets on Testing Set IV. The network calculations reveal that the long-range interaction changes contributed most to flexibility. In addition, the hydrogen bonds in the base–base interactions greatly stabilize the RNA structure, while backbone interactions determine RNA folding. The computational analysis of pocket flexibility could facilitate RNA engineering for biological or medical applications.

Comparative Assessment of Docking Programs for Docking and Virtual Screening of Ribosomal Oxazolidinone Antibacterial Agents

Oxazolidinones are a broad-spectrum class of synthetic antibiotics that bind to the 50S ribosomal subunit of Gram-positive and Gram-negative bacteria. Many crystal structures of the ribosomes with oxazolidinone ligands have been reported in the literature, facilitating structure-based design using methods such as molecular docking. It would be of great interest to know in advance how well docking methods can reproduce the correct ligand binding modes and rank these correctly. We examined the performance of five molecular docking programs (AutoDock 4, AutoDock Vina, DOCK 6, rDock, and RLDock) for their ability to model ribosomal-ligand interactions with oxazolidinones. Eleven ribosomal crystal structures with oxazolidinones as the ligands were docked. The accuracy was evaluated by calculating the docked complexes' root-mean-square deviation (RMSD) and the program's internal scoring function. The rankings for each program based on the median RMSD between the native and predicted were DOCK 6 > AD4 > Vina > RDOCK >> RLDOCK. Results demonstrate that the top-performing program, DOCK 6, could accurately replicate the ligand binding in only four of the eleven ribosomes due to the poor electron density of said ribosomal structures. In this study, we have further benchmarked the performance of the DOCK 6 docking algorithm and scoring in improving virtual screening (VS) enrichment using the dataset of 285 oxazolidinone derivatives against oxazolidinone binding sites in the S. aureus ribosome. However, there was no clear trend between the structure and activity of the oxazolidinones in VS. Overall, the docking performance indicates that the RNA pocket's high flexibility does not allow for accurate docking prediction, highlighting the need to validate VS. protocols for ligand-RNA before future use. Later, we developed a re-scoring method incorporating absolute docking scores and molecular descriptors, and the results indicate that the descriptors greatly improve the correlation of docking scores and pMIC values. Morgan fingerprint analysis was also used, suggesting that DOCK 6 underpredicted molecules with tail modifications with acetamide, n-methylacetamide, or n-ethylacetamide and over-predicted molecule derivatives with methylamino bits. Alternatively, a ligand-based approach similar to a field template was taken, indicating that each derivative's tail groups have strong positive and negative electrostatic potential contributing to microbial activity. These results indicate that one should perform VS. campaigns of ribosomal antibiotics with care and that more comprehensive strategies, including molecular dynamics simulations and relative free energy calculations, might be necessary in conjunction with VS. and docking.

Dynamic docking of small molecules targeting RNA CUG repeats causing myotonic dystrophy type 1

Expansion of RNA CUG repeats causes myotonic dystrophy type 1 (DM1). Once transcribed, the expanded CUG repeats strongly attract muscleblind-like 1 (MBNL1) proteins and disturb their functions in cells. Because of its unique structural form, expanded RNA CUG repeats are prospective drug targets, where small molecules can be utilized to target RNA CUG repeats to inhibit MBNL1 binding and ameliorate DM1-associated defects. In this contribution, we developed two physics-based dynamic docking approaches (DynaD and DynaD/Auto) and applied them to nine small molecules known to specifically target RNA CUG repeats. While DynaD uses a distance-based reaction coordinate to study the binding phenomenon, DynaD/Auto combines results of umbrella sampling calculations performed on 1 × 1 UU internal loops and AutoDock calculations to efficiently sample the energy landscape of binding. Predictions are compared with experimental data, displaying a positive correlation with correlation coefficient (R) values of 0.70 and 0.81 for DynaD and DynaD/Auto, respectively. Furthermore, we found that the best correlation was achieved with MM/3D-RISM calculations, highlighting the importance of solvation in binding calculations. Moreover, we detected that DynaD/Auto performed better than DynaD because of the use of prior knowledge about the binding site arising from umbrella sampling calculations. Finally, we developed dendrograms to present how bound states are connected to each other in a binding process. Results are exciting, as DynaD and DynaD/Auto will allow researchers to utilize two novel physics-based and computer-aided drug-design methodologies to perform in silico calculations on drug-like molecules aiming to target complex RNA loops.

Biomotors, viral assembly, and RNA nanobiotechnology: Current achievements and future directions

The International Society of RNA Nanotechnology and Nanomedicine (ISRNN) serves to further the development of a wide variety of functional nucleic acids and other related nanotechnology platforms. To aid in the dissemination of the most recent advancements, a biennial discussion focused on biomotors, viral assembly, and RNA nanobiotechnology has been established where international experts in interdisciplinary fields such as structural biology, biophysical chemistry, nanotechnology, cell and cancer biology, and pharmacology share their latest accomplishments and future perspectives. The results summarized here highlight advancements in our understanding of viral biology and the structure-function relationship of frame-shifting elements in genomic viral RNA, improvements in the predictions of SHAPE analysis of 3D RNA structures, and the understanding of dynamic RNA structures through a variety of experimental and computational means. Additionally, recent advances in the drug delivery, vaccine design, nanopore technologies, biomotor and biomachine development, DNA packaging, RNA nanotechnology, and drug delivery are included in this critical review. We emphasize some of the novel accomplishments, major discussion topics, and present current challenges and perspectives of these emerging fields.

An overview of structural approaches to study therapeutic RNAs

RNAs provide considerable opportunities as therapeutic agent to expand the plethora of classical therapeutic targets, from extracellular and surface proteins to intracellular nucleic acids and its regulators, in a wide range of diseases. RNA versatility can be exploited to recognize cell types, perform cell therapy, and develop new vaccine classes. Therapeutic RNAs (aptamers, antisense nucleotides, siRNA, miRNA, mRNA and CRISPR-Cas9) can modulate or induce protein expression, inhibit molecular interactions, achieve genome editing as well as exon-skipping. A common RNA thread, which makes it very promising for therapeutic applications, is its structure, flexibility, and binding specificity. Moreover, RNA displays peculiar structural plasticity compared to proteins as well as to DNA. Here we summarize the recent advances and applications of therapeutic RNAs, and the experimental and computational methods to analyze their structure, by biophysical techniques (liquid-state NMR, scattering, reactivity, and computational simulations), with a focus on dynamic and flexibility aspects and to binding analysis. This will provide insights on the currently available RNA therapeutic applications and on the best techniques to evaluate its dynamics and reactivity.

SILCS-RNA: Toward a Structure-Based Drug Design Approach for Targeting RNAs with Small Molecules

RNA molecules can act as potential drug targets in different diseases, as their dysregulated expression or misfolding can alter various cellular processes. Noncoding RNAs account for ∼70% of the human genome, and these molecules can have complex tertiary structures that present a great opportunity for targeting by small molecules. In the present study, the site identification by ligand competitive saturation (SILCS) computational approach is extended to target RNA, termed SILCS-RNA. Extensions to the method include an enhanced oscillating excess chemical potential protocol for the grand canonical Monte Carlo calculations and individual simulations of the neutral and charged solutes from which the SILCS functional group affinity maps (FragMaps) are calculated for subsequent binding site identification and docking calculations. The method is developed and evaluated against seven RNA targets and their reported small molecule ligands. SILCS-RNA provides a detailed characterization of the functional group affinity pattern in the small molecule binding sites, recapitulating the types of functional groups present in the ligands. The developed method is also shown to be useful for identification of new potential binding sites and identifying ligand moieties that contribute to binding, granular information that can facilitate ligand design. However, limitations in the method are evident including the ability to map the regions of binding sites occupied by ligand phosphate moieties and to fully account for the wide range of conformational heterogeneity in RNA associated with binding of different small molecules, emphasizing inherent challenges associated with applying computer-aided drug design methods to RNA. While limitations are present, the current study indicates how the SILCS-RNA approach may enhance drug discovery efforts targeting RNAs with small molecules.

Rosavin Ameliorates Hepatic Inflammation and Fibrosis in the NASH Rat Model via Targeting Hepatic Cell Death

BACKGROUND

Non-alcoholic fatty liver disease (NAFLD) represents the most common form of chronic liver disease that urgently needs effective therapy. Rosavin, a major constituent of the Rhodiola Rosea plant of the family Crassulaceae, is believed to exhibit multiple pharmacological effects on diverse diseases. However, its effect on non-alcoholic steatohepatitis (NASH), the progressive form of NAFLD, and the underlying mechanisms are not fully illustrated.

AIM

Investigate the pharmacological activity and potential mechanism of rosavin treatment on NASH management via targeting hepatic cell death-related (HSPD1/TNF/MMP14/ITGB1) mRNAs and their upstream noncoding RNA regulators (miRNA-6881-5P and lnc-SPARCL1-1:2) in NASH rats.

RESULTS

High sucrose high fat (HSHF) diet-induced NASH rats were treated with different concentrations of rosavin (10, 20, and 30 mg/kg/day) for the last four weeks of dietary manipulation. The data revealed that rosavin had the ability to modulate the expression of the hepatic cell death-related RNA panel through the upregulation of both (HSPD1/TNF/MMP14/ITGB1) mRNAs and their epigenetic regulators (miRNA-6881-5P and lnc-SPARCL1-1:2). Moreover, rosavin ameliorated the deterioration in both liver functions and lipid profile, and thereby improved the hepatic inflammation, fibrosis, and apoptosis, as evidenced by the decreased protein levels of IL6, TNF-α, and caspase-3 in liver sections of treated animals compared to the untreated NASH rats.

CONCLUSION

Rosavin has demonstrated a potential ability to attenuate disease progression and inhibit hepatic cell death in the NASH animal model. The produced effect was correlated with upregulation of the hepatic cell death-related (HSPD1, TNF, MMP14, and ITGB1) mRNAs-(miRNA-6881-5P-(lnc-SPARCL1-1:2) RNA panel.

Protein-Based Virtual Screening Tools Applied for RNA-Ligand Docking Identify New Binders of the preQ1-Riboswitch

Targeting RNA with small molecules is an emerging field. While several ligands for different RNA targets are reported, structure-based virtual screenings (VSs) against RNAs are still rare. Here, we elucidated the general capabilities of protein-based docking programs to reproduce native binding modes of small-molecule RNA ligands and to discriminate known binders from decoys by the scoring function. The programs were found to perform similar compared to the RNA-based docking tool rDOCK, and the challenges faced during docking, namely, protomer and tautomer selection, target dynamics, and explicit solvent, do not largely differ from challenges in conventional protein-ligand docking. A prospective VS with the Bacillus subtilis preQ₁-riboswitch aptamer domain performed with FRED, HYBRID, and FlexX followed by microscale thermophoresis assays identified six active compounds out of 23 tested VS hits with potencies between 29.5 nM and 11.0 μM. The hits were selected not solely based on their docking score but for resembling key interactions of the native ligand. Therefore, this study demonstrates the general feasibility to perform structure-based VSs against RNA targets, while at the same time it highlights pitfalls and their potential solutions when executing RNA-ligand docking.

RNA-ligand molecular docking: advances and challenges

With rapid advances in computer algorithms and hardware, fast and accurate virtual screening has led to a drastic acceleration in selecting potent small molecules as drug candidates. Computational modeling of RNA-small molecule interactions has become an indispensable tool for RNA-targeted drug discovery. The current models for RNA-ligand binding have mainly focused on the docking-and-scoring method. Accurate docking and scoring should tackle four crucial problems: (1) conformational flexibility of ligand, (2) conformational flexibility of RNA, (3) efficient sampling of binding sites and binding poses, and (4) accurate scoring of different binding modes. Moreover, compared with the problem of protein-ligand docking, predicting ligand binding to RNA, a negatively charged polymer, is further complicated by additional effects such as metal ion effects. Thermodynamic models based on physics-based and knowledge-based scoring functions have shown highly encouraging success in predicting ligand binding poses and binding affinities. Recently, kinetic models for ligand binding have further suggested that including dissociation kinetics (residence time) in ligand docking would result in improved performance in estimating in vivo drug efficacy. More recently, the rise of deep-learning approaches has led to new tools for predicting RNA-small molecule binding. In this review, we present an overview of the recently developed computational methods for RNA-ligand docking and their advantages and disadvantages.

RLDOCK method for predicting RNA-small molecule binding modes

RNA molecules play critical roles in cellular functions at the level of gene expression and regulation. The intricate 3D structures and the functional roles of RNAs make RNA molecules ideal targets for therapeutic drugs. The rational design of RNA-targeted drug requires accurate modeling of RNA-ligand interactions. Recently a new computational tool, RLDOCK, was developed to predict ligand binding sites and binding poses. Using an iterative multiscale sampling and search algorithm and a energy-based evaluation of ligand poses, the method enables efficient and accurate predictions for RNA-ligand interactions. Here we present a detailed illustration of the computational procedure for the practical implementation of the RLDOCK method. Using Flavin mononucleotide (FMN) docking to F. nucleatum FMN riboswitch as an example, we illustrate the computational protocol for RLDOCK-based prediction of RNA- ligand interactions. The RLDOCK software is freely accessible at http://https://github.com/Vfold-RNA/RLDOCK.

Docking and scoring for nucleic acid-ligand interactions: Principles and current status

Nucleic acid (NA)-ligand interactions have crucial roles in many cellular processes and, thus, are increasingly attracting therapeutic interest in drug discovery. Molecular docking is a valuable tool for studying molecular interactions. However, because NAs differ significantly from proteins in both their physical and chemical properties, traditional docking algorithms and scoring functions for protein-ligand interactions might not be applicable to NA-ligand docking. Therefore, various sampling strategies and scoring functions for NA-ligand interactions have been developed. Here, we review the basic principles and current status of docking algorithms and scoring functions for DNA/RNA-ligand interactions. We also discuss challenges and limitations of current docking and scoring approaches.

NLDock: a Fast Nucleic Acid-Ligand Docking Algorithm for Modeling RNA/DNA-Ligand Complexes

Nucleic acid-ligand interactions play an important role in numerous cellular processes such as gene function expression and regulation. Therefore, nucleic acids such as RNAs have become more and more important drug targets, where the structural determination of nucleic acid-ligand complexes is pivotal for understanding their functions and thus developing therapeutic interventions. Molecular docking has been a useful computational tool in predicting the complex structure between molecules. However, although a number of docking algorithms have been developed for protein-ligand interactions, only a few docking programs were presented for nucleic acid-ligand interactions. Here, we have developed a fast nucleic acid-ligand docking algorithm, named NLDock, by implementing our intrinsic scoring function ITScoreNL for nucleic acid-ligand interactions into a modified version of the MDock program. NLDock was extensively evaluated on four test sets and compared with five other state-of-the-art docking algorithms including AutoDock, DOCK 6, rDock, GOLD, and Glide. It was shown that our NLDock algorithm obtained a significantly better performance than the other docking programs in binding mode predictions and achieved the success rates of 73%, 36%, and 32% on the largest test set of 77 complexes for local rigid-, local flexible-, and global flexible-ligand docking, respectively. In addition, our NLDock approach is also computationally efficient and consumed an average of as short as 0.97 and 2.08 min for a local flexible-ligand docking job and a global flexible-ligand docking job, respectively. These results suggest the good performance of our NLDock in both docking accuracy and computational efficiency.

Thermostability, Tunability, and Tenacity of RNA as Rubbery Anionic Polymeric Materials in Nanotechnology and Nanomedicine-Specific Cancer Targeting with Undetectable Toxicity

RNA nanotechnology is the bottom-up self-assembly of nanometer-scale architectures, resembling LEGOs, composed mainly of RNA. The ideal building material should be (1) versatile and controllable in shape and stoichiometry, (2) spontaneously self-assemble, and (3) thermodynamically, chemically, and enzymatically stable with a long shelf life. RNA building blocks exhibit each of the above. RNA is a polynucleic acid, making it a polymer, and its negative-charge prevents nonspecific binding to negatively charged cell membranes. The thermostability makes it suitable for logic gates, resistive memory, sensor set-ups, and NEM devices. RNA can be designed and manipulated with a level of simplicity of DNA while displaying versatile structure and enzyme activity of proteins. RNA can fold into single-stranded loops or bulges to serve as mounting dovetails for intermolecular or domain interactions without external linking dowels. RNA nanoparticles display rubber- and amoeba-like properties and are stretchable and shrinkable through multiple repeats, leading to enhanced tumor targeting and fast renal excretion to reduce toxicities. It was predicted in 2014 that RNA would be the third milestone in pharmaceutical drug development. The recent approval of several RNA drugs and COVID-19 mRNA vaccines by FDA suggests that this milestone is being realized. Here, we review the unique properties of RNA nanotechnology, summarize its recent advancements, describe its distinct attributes inside or outside the body and discuss potential applications in nanotechnology, medicine, and material science.

Design and Implementation of Synthetic RNA Binders for the Inhibition of miR-21 Biogenesis

Targeting RNAs using small molecules is an emerging field of medicinal chemistry and holds promise for the discovery of efficient tools for chemical biology. MicroRNAs are particularly interesting targets since they are involved in a number of pathologies such as cancers. Indeed, overexpressed microRNAs in cancer are oncogenic and various series of inhibitors of microRNAs biogenesis have been developed in recent years. Here, we describe the structure-based design of new efficient inhibitors of microRNA-21. Starting from a previously identified hit, we performed biochemical studies and molecular docking to design a new series of optimized conjugates of neomycin aminoglycoside with artificial nucleobases and amino acids. Investigation about the mode of action and the site of the interaction of the newly synthesized compounds allowed for the description of structure-activity relationships and the identification of the most important parameters for miR-21 inhibition.