Improved Accuracy in RNA-Protein Rigid Body Docking by Incorporating Force Field for Molecular Dynamics Simulation into the Scoring Function

Drug Repositioning Using Computer-aided Drug Design (CADD)

Drug repositioning is a method of using authorized drugs for other unusually complex diseases. Compared to new drug development, this method is fast, low in cost, and effective. Through the use of outstanding bioinformatics tools, such as computer-aided drug design (CADD), computer strategies play a vital role in the re-transformation of drugs. The use of CADD's special strategy for target-based drug reuse is the most promising method, and its realization rate is high. In this review article, we have particularly focused on understanding the various technologies of CADD and the use of computer-aided drug design for target-based drug reuse, taking COVID-19 and cancer as examples. Finally, it is concluded that CADD technology is accelerating the development of repurposed drugs due to its many advantages, and there are many facts to prove that the new ligand-targeting strategy is a beneficial method and that it will gain momentum with the development of technology.

Challenges in structural modeling of RNA-protein interactions

In the past few years, the number of RNA-binding proteins (RBP) and RNA-RBP interactions has increased significantly. Here, we review recent developments in the methodology for protein-RNA and protein-protein complex structure modeling with deep learning and co-evolution, as well as discuss the challenges and opportunities for building a reliable approach for protein-RNA complex structure modelling. Protein Data bank (PDB) and Cross-linking immunoprecipitation (CLIP) data could be combined together and used to infer 2D geometry of protein-RNA interactions by deep learning.

Smart Approach for the Design of Highly Selective Aptamer-Based Biosensors

Aptamers are chemically synthesized single-stranded DNA or RNA oligonucleotides widely used nowadays in sensors and nanoscale devices as highly sensitive biorecognition elements. With proper design, aptamers are able to bind to a specific target molecule with high selectivity. To date, the systematic evolution of ligands by exponential enrichment (SELEX) process is employed to isolate aptamers. Nevertheless, this method requires complex and time-consuming procedures. In silico methods comprising machine learning models have been recently proposed to reduce the time and cost of aptamer design. In this work, we present a new in silico approach allowing the generation of highly sensitive and selective RNA aptamers towards a specific target, here represented by ammonium dissolved in water. By using machine learning and bioinformatics tools, a rational design of aptamers is demonstrated. This “smart” SELEX method is experimentally proved by choosing the best five aptamer candidates obtained from the design process and applying them as functional elements in an electrochemical sensor to detect, as the target molecule, ammonium at different concentrations. We observed that the use of five different aptamers leads to a significant difference in the sensor’s response. This can be explained by considering the aptamers’ conformational change due to their interaction with the target molecule. We studied these conformational changes using a molecular dynamics simulation and suggested a possible explanation of the experimental observations. Finally, electrochemical measurements exposing the same sensors to different molecules were used to confirm the high selectivity of the designed aptamers. The proposed in silico SELEX approach can potentially reduce the cost and the time needed to identify the aptamers and potentially be applied to any target molecule.

Statistical potentials for RNA-protein interactions optimized by CMA-ES

Characterizing RNA-protein interactions remains an important endeavor, complicated by the difficulty in obtaining the relevant structures. Evaluating model structures via statistical potentials is in principle straight-forward and effective. However, given the relatively small size of the existing learning set of RNA-protein complexes optimization of such potentials continues to be problematic. Notably, interaction-based statistical potentials have problems in addressing large RNA-protein complexes. In this study, we adopted a novel strategy with covariance matrix adaptation (CMA-ES) to calculate statistical potentials, successfully identifying native docking poses.

De novo ssRNA Aptamers against the SARS-CoV-2 Main Protease: In Silico Design and Molecular Dynamics Simulation

Herein, we have generated ssRNA aptamers to inhibit SARS-CoV-2 M^pro, a protease necessary for the SARS-CoV-2 coronavirus replication. Because there is no aptamer 3D structure currently available in the databanks for this protein, first, we modeled an ssRNA aptamer using an entropic fragment-based strategy. We refined the initial sequence and 3D structure by using two sequential approaches, consisting of an elitist genetic algorithm and an RNA inverse process. We identified three specific aptamers against SARS-CoV-2 M^pro, called MApta^pro, MApta^pro-IR1, and MApta^pro-IR2, with similar 3D conformations and that fall in the dimerization region of the SARS-CoV-2 M^pro necessary for the enzymatic activity. Through the molecular dynamic simulation and binding free energy calculation, the interaction between the MApta^pro-IR1 aptamer and the SARS-CoV-2 M^pro enzyme resulted in the strongest and the highest stable complex; therefore, the ssRNA MApta^pro-IR1 aptamer was selected as the best potential candidate for the inhibition of SARS-CoV-2 M^pro and a perspective therapeutic drug for the COVID-19 disease.

Cyclodipeptide Synthases of the NYH Subfamily Recognize tRNA Using an α-Helix Enriched with Positive Residues

Cyclodipeptide synthases (CDPSs) perform nonribosomal protein synthesis using two aminoacyl-tRNA substrates to produce cyclodipeptides. At present, there are no structural details of the CDPS:tRNA interaction available. Using AlbC, a CDPS that produces cyclo(l-Phe-l-Phe), the interaction between AlbC and its Phe-tRNA substrate was investigated. Simulations of models of the AlbC:tRNA complex, proposed by rigid-body docking or homology modeling, demonstrated that interactions with residues of an AlbC α-helix, α4, significantly contribute to the free energy of binding of AlbC to tRNA. Individual residue contributions to the tRNA binding free energy of the discovered binding mode explain well the available biochemical data, and the results of in vivo assay experiments performed in this work and guided by simulations. In molecular dynamics simulations, the phenylalanyl group predominantly occupied the two positions observed in the experimental structure of AlbC in the dipeptide intermediate state, suggesting that tRNAs of the first and second substrates interact with AlbC in a similar manner. Overall, given the high degree of sequence and structural similarity among the members of the CDPS NYH protein subfamily, the mechanism of the protein:tRNA interaction is expected to be pertinent to a wide range of proteins interacting with tRNA.

P3DOCK: a protein-RNA docking webserver based on template-based and template-free docking

MOTIVATION

The main function of protein-RNA interaction is to regulate the expression of genes. Therefore, studying protein-RNA interactions is of great significance. The information of three-dimensional (3D) structures reveals that atomic interactions are particularly important. The calculation method for modeling a 3D structure of a complex mainly includes two strategies: free docking and template-based docking. These two methods are complementary in protein-protein docking. Therefore, integrating these two methods may improve the prediction accuracy.

RESULTS

In this article, we compare the difference between the free docking and the template-based algorithm. Then we show the complementarity of these two methods. Based on the analysis of the calculation results, the transition point is confirmed and used to integrate two docking algorithms to develop P3DOCK. P3DOCK holds the advantages of both algorithms. The results of the three docking benchmarks show that P3DOCK is better than those two non-hybrid docking algorithms. The success rate of P3DOCK is also higher (3-20%) than state-of-the-art hybrid and non-hybrid methods. Finally, the hierarchical clustering algorithm is utilized to cluster the P3DOCK's decoys. The clustering algorithm improves the success rate of P3DOCK. For ease of use, we provide a P3DOCK webserver, which can be accessed at www.rnabinding.com/P3DOCK/P3DOCK.html. An integrated protein-RNA docking benchmark can be downloaded from http://rnabinding.com/P3DOCK/benchmark.html.

AVAILABILITY AND IMPLEMENTATION

www.rnabinding.com/P3DOCK/P3DOCK.html.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Protein-ensemble-RNA docking by efficient consideration of protein flexibility through homology models

MOTIVATION

Given the importance of protein-ribonucleic acid (RNA) interactions in many biological processes, a variety of docking algorithms have been developed to predict the complex structure from individual protein and RNA partners in the past decade. However, due to the impact of molecular flexibility, the performance of current methods has hit a bottleneck in realistic unbound docking. Pushing the limit, we have proposed a protein-ensemble-RNA docking strategy to explicitly consider the protein flexibility in protein-RNA docking through an ensemble of multiple protein structures, which is referred to as MPRDock. Instead of taking conformations from MD simulations or experimental structures, we obtained the multiple structures of a protein by building models from its homologous templates in the Protein Data Bank (PDB).

RESULTS

Our approach can not only avoid the reliability issue of structures from MD simulations but also circumvent the limited number of experimental structures for a target protein in the PDB. Tested on 68 unbound-bound and 18 unbound-unbound protein-RNA complexes, our MPRDock/DITScorePR considerably improved the docking performance and achieved a significantly higher success rate than single-protein rigid docking whether pseudo-unbound templates are included or not. Similar improvements were also observed when combining our ensemble docking strategy with other scoring functions. The present homology model-based ensemble docking approach will have a general application in molecular docking for other interactions.

AVAILABILITY AND IMPLEMENTATION

http://huanglab.phys.hust.edu.cn/mprdock/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

In silico design and validation of high-affinity RNA aptamers targeting epithelial cellular adhesion molecule dimers

Nucleic acid aptamers hold great promise for therapeutic applications due to their favorable intrinsic properties, as well as high-throughput experimental selection techniques. Despite the utility of the systematic evolution of ligands by the exponential enrichment (SELEX) method for aptamer determination, complementary in silico aptamer design is highly sought after to facilitate virtual screening and increased understanding of important nucleic acid-protein interactions. Here, with a combined experimental and theoretical approach, we have developed two optimal epithelial cellular adhesion molecule (EpCAM) aptamers. Our structure-based in silico method first predicts their binding modes and then optimizes them for EpCAM with molecular dynamics simulations, docking, and free energy calculations. Our isothermal titration calorimetry experiments further confirm that the EpCAM aptamers indeed exhibit enhanced affinity over a previously patented nanomolar aptamer, EP23. Moreover, our study suggests that EP23 and the de novo designed aptamers primarily bind to EpCAM dimers (and not monomers, as hypothesized in previous published works), suggesting a paradigm for developing EpCAM-targeted therapies.

PredPRBA: Prediction of Protein-RNA Binding Affinity Using Gradient Boosted Regression Trees

Protein-RNA interactions play essential roles in many biological aspects. Quantifying the binding affinity of protein-RNA complexes is helpful to the understanding of protein-RNA recognition mechanisms and identification of strong binding partners. Due to experimentally measured protein-RNA binding affinity data available is still limited to date, there is a pressing demand for accurate and reliable computational approaches. In this paper, we propose a computational approach, PredPRBA, which can effectively predict protein-RNA binding affinity using gradient boosted regression trees. We build a dataset of protein-RNA binding affinity that includes 103 protein-RNA complex structures manually collected from related literature. Then, we generate 37 kinds of sequence and structural features and explore the relationship between the features and protein-RNA binding affinity. We find that the binding affinity mainly depends on the structure of RNA molecules. According to the type of RNA associated with proteins composed of the protein-RNA complex, we split the 103 protein-RNA complexes into six categories. For each category, we build a gradient boosted regression tree (GBRT) model based on the generated features. We perform a comprehensive evaluation for the proposed method on the binding affinity dataset using leave-one-out cross-validation. We show that PredPRBA achieves correlations ranging from 0.723 to 0.897 among six categories, which is significantly better than other typical regression methods and the pioneer protein-RNA binding affinity predictor SPOT-Seq-RNA. In addition, a user-friendly web server has been developed to predict the binding affinity of protein-RNA complexes. The PredPRBA webserver is freely available at http://PredPRBA.denglab.org/.

Long non-coding RNAs in hematological malignancies: translating basic techniques into diagnostic and therapeutic strategies

Recent studies have revealed that non-coding regions comprise the vast majority of the human genome and long non-coding RNAs (lncRNAs) are a diverse class of non-coding RNAs that has been implicated in a variety of biological processes. Abnormal expression of lncRNAs has also been linked to different human diseases including cancers, yet the regulatory mechanisms and functional effects of lncRNAs are still ambiguous, and the molecular details also need to be confirmed. Unlike protein-coding gene, it is much more challenging to unravel the roles of lncRNAs owing to their unique and complex features such as functional diversity and low conservation among species, which greatly hamper their experimental characterization. In this review, we summarize and discuss both conventional and advanced approaches for the identification and functional characterization of lncRNAs related to hematological malignancies. In particular, the utility and advancement of clustered regularly interspaced short palindromic repeats (CRISPR)-Cas system as gene-editing tools are envisioned to facilitate the molecular dissection of lncRNAs via different knock-in/out strategies. Besides experimental considerations specific to lncRNAs, the roles of lncRNAs in the pathogenesis and progression of leukemia are also highlighted in the review. We expect that these insights may ultimately lead to clinical applications including development of biomarkers and novel therapeutic approaches targeting lncRNAs.

Assessing the performance of MM/PBSA and MM/GBSA methods. 8. Predicting binding free energies and poses of protein-RNA complexes

Molecular docking provides a computationally efficient way to predict the atomic structural details of protein-RNA interactions (PRI), but accurate prediction of the three-dimensional structures and binding affinities for PRI is still notoriously difficult, partly due to the unreliability of the existing scoring functions for PRI. MM/PBSA and MM/GBSA are more theoretically rigorous than most scoring functions for protein-RNA docking, but their prediction performance for protein-RNA systems remains unclear. Here, we systemically evaluated the capability of MM/PBSA and MM/GBSA to predict the binding affinities and recognize the near-native binding structures for protein-RNA systems with different solvent models and interior dielectric constants (ε_in). For predicting the binding affinities, the predictions given by MM/GBSA based on the minimized structures in explicit solvent and the GB^GBn1 model with ε_in = 2 yielded the highest correlation with the experimental data. Moreover, the MM/GBSA calculations based on the minimized structures in implicit solvent and the GB^GBn1 model distinguished the near-native binding structures within the top 10 decoys for 117 out of the 148 protein-RNA systems (79.1%). This performance is better than all docking scoring functions studied here. Therefore, the MM/GBSA rescoring is an efficient way to improve the prediction capability of scoring functions for protein-RNA systems.

Bioinformatics Tools and Benchmarks for Computational Docking and 3D Structure Prediction of RNA-Protein Complexes

RNA-protein (RNP) interactions play essential roles in many biological processes, such as regulation of co-transcriptional and post-transcriptional gene expression, RNA splicing, transport, storage and stabilization, as well as protein synthesis. An increasing number of RNP structures would aid in a better understanding of these processes. However, due to the technical difficulties associated with experimental determination of macromolecular structures by high-resolution methods, studies on RNP recognition and complex formation present significant challenges. As an alternative, computational prediction of RNP interactions can be carried out. Structural models obtained by theoretical predictive methods are, in general, less reliable compared to models based on experimental measurements but they can be sufficiently accurate to be used as a basis for to formulating functional hypotheses. In this article, we present an overview of computational methods for 3D structure prediction of RNP complexes. We discuss currently available methods for macromolecular docking and for scoring 3D structural models of RNP complexes in particular. Additionally, we also review benchmarks that have been developed to assess the accuracy of these methods.

Protein-RNA Docking Using ICM

Protein-RNA interactions play an important role in many biological processes. Computational methods such as docking have been developed to complement existing biophysical and structural biology techniques. Computational prediction of protein-RNA complex structures includes two steps: generating candidate structures from the individual protein and RNA parts and scoring the generated poses to pick out the correct one. In this work, we considered three recently developed data sets of protein-RNA complexes to evaluate and improve the performance of the FFT-based rigid-body docking algorithm implemented in the ICM package. An electrostatic term describing interactions between negatively charged phosphate groups and positively charged protein residues was added to the energy function used during the docking step to take into account the greater role that electrostatic interactions play in protein-RNA complexes. Next, the docking results were used to optimize a scoring function including van der Waals, electrostatic, and solvation terms. This optimization yielded a much smaller weight for the solvation term indicating that solvation energy may be less important for the scoring of protein-RNA structures. Rescoring of the generated poses with the new scoring function led to much higher success rates, while pose clustering by contact fingerprints produced further improvements, achieving a success rate of 0.66 for the top 100 structures.

The design of novel inhibitors for treating cancer by targeting CDC25B through disruption of CDC25B-CDK2/Cyclin A interaction using computational approaches

Cell division cycle 25B is a key cell cycle regulator and widely considered as potent clinical drug target for cancers. This research focused on identifying potential compounds in theory which are able to disrupt transient interactions between CDC25B and its CDK2/Cyclin A substrate.By using the method of ZDOCK and RDOCK, the most optimized 3D structure of CDK2/Cyclin A in complex with CDC25B was constructed and validated using two methods: 1) the superimposition of proteins; 2) analysis of the hydrogen bond distances of Arg 488(N1)-Asp 206(OD1), Arg 492(NE)-Asp 206(OD1), Arg 492(N1)-Asp 206(OD2) and Tyr 497(NE)-Asp 210(OD1). A series of new compounds was gained through searching the fragment database derived from ZINC based on the known inhibitor-compound 7 by the means of "replace fragment" technique. The compounds acquired via meeting the requirements of the absorption, distribution, metabolism, and excretion (ADME) predictions. Finally, 12 compounds with better binding affinity were identified. The comp#1, as a representative, was selected to be synthesized and assayed for their CDC25B inhibitory activities. The comp#1 exhibited mild inhibitory activities against human CDC25B with IC50 values at about 39.02 μM. Molecular Dynamic (MD) simulation revealed that the new inhibitor-comp#1 had favorable conformations for binding to CDC25B and disturbing the interactions between CDC25B and CDK2/Cyclin A.

In silico approaches to RNA aptamer design

RNA aptamers are ribonucleic acids that bind to specific target molecules. An RNA aptamer for a disease-related protein has great potential for development into a new drug. However, huge time and cost investments are required to develop an RNA aptamer into a pharmaceutical. Recently, SELEX combined with high-throughput sequencers (i.e., HT-SELEX) has been widely used to select candidate RNA aptamers that bind to a target protein with high affinity and specificity. After candidate selection, further optimizations such as shortening and modifying candidate sequences are performed. In these steps, in silico approaches are expected to reduce the time and cost associated with aptamer drug development. In this article, we review existing in silico approaches to RNA aptamer development, including a method for ranking the candidates of RNA aptamers from HT-SELEX data, clustering a huge number of aptamer sequences, and finding motifs amidst a set of significant RNA aptamers. It is expected that further studies in addition to these methods will be utilized for in silico RNA aptamer design, permitting a minimal number of experiments to be performed through the utilization of sophisticated computational methods.

A non-redundant protein-RNA docking benchmark version 2.0

We present an updated version of the protein-RNA docking benchmark, which we first published four years back. The non-redundant protein-RNA docking benchmark version 2.0 consists of 126 test cases, a threefold increase in number compared to its previous version. The present version consists of 21 unbound-unbound cases, of which, in 12 cases, the unbound RNAs are taken from another complex. It also consists of 95 unbound-bound cases where only the protein is available in the unbound state. Besides, we introduce 10 new bound-unbound cases where only the RNA is found in the unbound state. Based on the degree of conformational change of the interface residues upon complex formation the benchmark is classified into 72 rigid-body cases, 25 semiflexible cases and 19 full flexible cases. It also covers a wide range of conformational flexibility including small side chain movement to large domain swapping in protein structures as well as flipping and restacking in RNA bases. This benchmark should provide the docking community with more test cases for evaluating rigid-body as well as flexible docking algorithms. Besides, it will also facilitate the development of new algorithms that require large number of training set. The protein-RNA docking benchmark version 2.0 can be freely downloaded from http://www.csb.iitkgp.ernet.in/applications/PRDBv2. Proteins 2017; 85:256-267. © 2016 Wiley Periodicals, Inc.

Structural and energy determinants in protein-RNA docking

Deciphering the structural and energetic determinants of protein-RNA interactions harbors the potential to understand key cell processes at molecular level, such as gene expression and regulation. With this purpose, computational methods like docking aim to complement current biophysical and structural biology efforts. However, the few reported docking algorithms for protein-RNA interactions show limited predictive success rates, mainly due to incomplete sampling of the conformational space of both the protein and the RNA molecules, as well as to the difficulties of the scoring function in identifying the correct docking models. Here, we have tested the predictive value of a variety of knowledge-based and energetic scoring functions on a recently published protein-RNA docking benchmark and developed a scoring function able to efficiently discriminate docking decoys. We first performed docking calculations with the bound conformation, which allowed us to analyze the problem in optimal conditions. We found that geometry-based terms and electrostatics were the most important scoring terms, while binding propensities and desolvation were much less relevant for the scoring of protein-RNA models. This is in contrast with what we observed for protein-protein docking. The results also showed an interesting dependence of the predictive rates on the flexibility of the protein molecule, which arises from the observed higher positive charge of flexible interfaces and provides hints for future development of more efficient protein-RNA docking methods.