Exploring the RNA-Recognition Mechanism Using Supervised Molecular Dynamics (SuMD) Simulations: Toward a Rational Design for Ribonucleic-Targeting Molecules?

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.

Thermal titration molecular dynamics (TTMD): shedding light on the stability of RNA-small molecule complexes

Ribonucleic acids are gradually becoming relevant players among putative drug targets, thanks to the increasing amount of structural data exploitable for the rational design of selective and potent binders that can modulate their activity. Mainly, this information allows employing different computational techniques for predicting how well would a ribonucleic-targeting agent fit within the active site of its target macromolecule. Due to some intrinsic peculiarities of complexes involving nucleic acids, such as structural plasticity, surface charge distribution, and solvent-mediated interactions, the application of routinely adopted methodologies like molecular docking is challenged by scoring inaccuracies, while more physically rigorous methods such as molecular dynamics require long simulation times which hamper their conformational sampling capabilities. In the present work, we present the first application of Thermal Titration Molecular Dynamics (TTMD), a recently developed method for the qualitative estimation of unbinding kinetics, to characterize RNA-ligand complexes. In this article, we explored its applicability as a post-docking refinement tool on RNA in complex with small molecules, highlighting the capability of this method to identify the native binding mode among a set of decoys across various pharmaceutically relevant test cases.

Fighting Antimicrobial Resistance: Insights on How the Staphylococcus aureus NorA Efflux Pump Recognizes 2-Phenylquinoline Inhibitors by Supervised Molecular Dynamics (SuMD) and Molecular Docking Simulations

The superbug Staphylococcus aureus (S. aureus) exhibits several resistance mechanisms, including efflux pumps, that strongly contribute to antimicrobial resistance. In particular, the NorA efflux pump activity is associated with S. aureus resistance to fluoroquinolone antibiotics (e.g., ciprofloxacin) by promoting their active extrusion from cells. Thus, since efflux pump inhibitors (EPIs) are able to increase antibiotic concentrations in bacteria as well as restore their susceptibility to these agents, they represent a promising strategy to counteract bacterial resistance. Additionally, the very recent release of two NorA efflux pump cryo-electron microscopy (cryo-EM) structures in complex with synthetic antigen-binding fragments (Fabs) represents a real breakthrough in the study of S. aureus antibiotic resistance. In this scenario, supervised molecular dynamics (SuMD) and molecular docking experiments were combined to investigate for the first time the molecular mechanisms driving the interaction between NorA and efflux pump inhibitors (EPIs), with the ultimate goal of elucidating how the NorA efflux pump recognizes its inhibitors. The findings provide insights into the dynamic NorA-EPI intermolecular interactions and lay the groundwork for future drug discovery efforts aimed at the identification of novel molecules to fight antimicrobial resistance.

RNA-protein complexes and force field polarizability

Molecular dynamic (MD) simulations offer a way to study biomolecular interactions and their dynamics at the atomistic level. There are only a few studies of RNA-protein complexes in MD simulations, and here we wanted to study how force fields differ when simulating RNA-protein complexes: 1) argonaute 2 with bound guide RNA and a target RNA, 2) CasPhi-2 bound to CRISPR RNA and 3) Retinoic acid-inducible gene I C268F variant in complex with double-stranded RNA. We tested three non-polarizable force fields: Amber protein force fields ff14SB and ff19SB with RNA force field OL3, and the all-atom OPLS4 force field. Due to the highly charged and polar nature of RNA, we also tested the polarizable AMOEBA force field and the ff19SB and OL3 force fields with a polarizable water model O3P. Our results show that the non-polarizable force fields lead to compact and stable complexes. The polarizability in the force field or in the water model allows significantly more movement from the complex, but in some cases, this results in the disintegration of the complex structure, especially if the protein contains longer loop regions. Thus, one should be cautious when running long-scale simulations with polarizability. As a conclusion, all the tested force fields can be used to simulate RNA-protein complexes and the choice of the optimal force field depends on the studied system and research question.

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.

Geneticin shows selective antiviral activity against SARS-CoV-2 by interfering with programmed -1 ribosomal frameshifting

SARS-CoV-2 is currently causing an unprecedented pandemic. While vaccines are massively deployed, we still lack effective large-scale antiviral therapies. In the quest for antivirals targeting conserved structures, we focused on molecules able to bind viral RNA secondary structures. Aminoglycosides are a class of antibiotics known to interact with the ribosomal RNA of both prokaryotes and eukaryotes and have previously been shown to exert antiviral activities by interacting with viral RNA. Here we show that the aminoglycoside geneticin is endowed with antiviral activity against all tested variants of SARS-CoV-2, in different cell lines and in a respiratory tissue model at non-toxic concentrations. The mechanism of action is an early inhibition of RNA replication and protein expression related to a decrease in the efficiency of the -1 programmed ribosomal frameshift (PRF) signal of SARS-CoV-2. Using in silico modeling, we have identified a potential binding site of geneticin in the pseudoknot of frameshift RNA motif. Moreover, we have selected, through virtual screening, additional RNA binding compounds, interacting with the same site with increased potency.

Investigating RNA-protein recognition mechanisms through supervised molecular dynamics (SuMD) simulations

Ribonucleic acid (RNA) plays a key regulatory role within the cell, cooperating with proteins to control the genome expression and several biological processes. Due to its characteristic structural features, this polymer can mold itself into different three-dimensional structures able to recognize target biomolecules with high affinity and specificity, thereby attracting the interest of drug developers and medicinal chemists. One successful example of the exploitation of RNA's structural and functional peculiarities is represented by aptamers, a class of therapeutic and diagnostic tools that can recognize and tightly bind several pharmaceutically relevant targets, ranging from small molecules to proteins, making use of the available structural and conformational freedom to maximize the complementarity with their interacting counterparts. In this scientific work, we present the first application of Supervised Molecular Dynamics (SuMD), an enhanced sampling Molecular Dynamics-based method for the study of receptor-ligand association processes in the nanoseconds timescale, to the study of recognition pathways between RNA aptamers and proteins, elucidating the main advantages and limitations of the technique while discussing its possible role in the rational design of RNA-based therapeutics.

Geneticin shows selective antiviral activity against SARS-CoV-2 by interfering with programmed -1 ribosomal frameshifting

SARS-CoV-2 is currently causing an unprecedented pandemic. While vaccines are massively deployed, we still lack effective large-scale antiviral therapies. In the quest for antivirals targeting conserved structures, we focused on molecules able to bind viral RNA secondary structures. Aminoglycosides are a class of antibiotics known to interact with the ribosomal RNA of both prokaryotes and eukaryotes and have previously been shown to exert antiviral activities by interacting with viral RNA. Here we show that the aminoglycoside geneticin is endowed with antiviral activity against all tested variants of SARS-CoV-2, in different cell lines and in a respiratory tissue model at non-toxic concentrations. The mechanism of action is an early inhibition of RNA replication and protein expression related to a decrease in the efficiency of the -1 programmed ribosomal frameshift (PRF) signal of SARS-CoV-2. Using in silico modelling, we have identified a potential binding site of geneticin in the pseudoknot of frameshift RNA motif. Moreover, we have selected, through virtual screening, additional RNA binding compounds, interacting with the same site with increased potency.

Protein DEK and DTA Aptamers: Insight Into the Interaction Mechanisms and the Computational Aptamer Design

By blocking the DEK protein, DEK-targeted aptamers (DTAs) can reduce the formation of neutrophil extracellular traps (NETs) to reveal a strong anti-inflammatory efficacy in rheumatoid arthritis. However, the poor stability of DTA has greatly limited its clinical application. Thus, in order to design an aptamer with better stability, DTA was modified by methoxy groups (DTA_OMe) and then the exact DEK–DTA interaction mechanisms were explored through theoretical calculations. The corresponding 2′-OCH3-modified nucleotide force field was established and the molecular dynamics (MD) simulations were performed. It was proved that the 2′-OCH3-modification could definitely enhance the stability of DTA on the premise of comparative affinity. Furthermore, the electrostatic interaction contributed the most to the binding of DEK–DTA, which was the primary interaction to maintain stability, in addition to the non-specific interactions between positively-charged residues (e.g., Lys and Arg) of DEK and the negatively-charged phosphate backbone of aptamers. The H-bond network analysis reminded that eight bases could be mutated to probably enhance the affinity of DTA_OMe. Therein, replacing the 29th base from cytosine to thymine of DTA_OMe was theoretically confirmed to be with the best affinity and even better stability. These research studies imply to be a promising new aptamer design strategy for the treatment of inflammatory arthritis.

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.

Shedding Light on the Molecular Recognition of Sub-Kilodalton Macrocyclic Peptides on Thrombin by Supervised Molecular Dynamics

Macrocycles are attractive structures for drug development due to their favorable structural features, potential in binding to targets with flat featureless surfaces, and their ability to disrupt protein–protein interactions. Moreover, large novel highly diverse libraries of low-molecular-weight macrocycles with therapeutically favorable characteristics have been recently established. Considering the mentioned facts, having a validated, fast, and accurate computational protocol for studying the molecular recognition and binding mode of this interesting new class of macrocyclic peptides deemed to be helpful as well as insightful in the quest of accelerating drug discovery. To that end, the ability of the in-house supervised molecular dynamics protocol called SuMD in the reproduction of the X-ray crystallography final binding state of a macrocyclic non-canonical tetrapeptide—from a novel library of 8,988 sub-kilodalton macrocyclic peptides—in the thrombin active site was successfully validated. A comparable binding mode with the minimum root-mean-square deviation (RMSD) of 1.4 Å at simulation time point 71.6 ns was achieved. This method validation study extended the application domain of the SuMD sampling method for computationally cheap, fast but accurate, and insightful macrocycle–protein molecular recognition studies.

Exploring Ligand Binding to Calcitonin Gene-Related Peptide Receptors

Class B1 G protein-coupled receptors (GPCRs) are important targets for many diseases, including cancer, diabetes, and heart disease. All the approved drugs for this receptor family are peptides that mimic the endogenous activating hormones. An understanding of how agonists bind and activate class B1 GPCRs is fundamental for the development of therapeutic small molecules. We combined supervised molecular dynamics (SuMD) and classic molecular dynamics (cMD) simulations to study the binding of the calcitonin gene-related peptide (CGRP) to the CGRP receptor (CGRPR). We also evaluated the association and dissociation of the antagonist telcagepant from the extracellular domain (ECD) of CGRPR and the water network perturbation upon binding. This study, which represents the first example of dynamic docking of a class B1 GPCR peptide, delivers insights on several aspects of ligand binding to CGRPR, expanding understanding of the role of the ECD and the receptor-activity modifying protein 1 (RAMP1) on agonist selectivity.

Probing the mechanism of peptide binding to REV response element RNA of HIV-1; MD simulations and free energy calculations

Ribonucleic acid (RNA) of HIV-1 contains a 350 nucleotide, highly structured, cis-acting element called RRE (REV-response-element RNA), essential for virus replication. REV is a natural peptide that binds to RRE and transports it from the nucleus to cytoplasm where it is expressed into a new virus. The synthetic peptide known as RSG-1.2 also binds the RRE element and competes with REV. The purpose of study is to rationally design novel peptides such as RSG peptide with improved binding affinity to prevent the transport of HIV-1 RNA and so replication of virus. Herein, we performed MD simulation and free energy calculations to evaluate the interactions and binding free energies of REV (PDB ID: 4PMI) and RSGs peptides (PDB IDs: 1G70 and 1I9F) with RRE. The protein-RNA interactions were analyzed using the MM-PBSA method. Results suggest that REV has more binding free energy -188.41 kcal/mol than two RSG peptides with total binding free energy -120.97 and -141.46 kcal/mol. The ARG and ASN were found to be important residues of REV. In the RRE sequence, the nucleotides 62-67 and 78-84 were found to be important contributors in binding free energy. This study play a major role in elaboration of binding REV and RSG1-2 with RRE element and pave the way for further synthesis of peptide that can bind with RRE element and can be selected as therapeutic agent for HIV.Communicated by Ramaswamy H. Sarma.