1
|
Rawat S, Subramaniam K, Subramanian SK, Subbarayan S, Dhanabalan S, Chidambaram SKM, Stalin B, Roy A, Nagaprasad N, Aruna M, Tesfaye JL, Badassa B, Krishnaraj R. Drug Repositioning Using Computer-aided Drug Design (CADD). Curr Pharm Biotechnol 2024; 25:301-312. [PMID: 37605405 DOI: 10.2174/1389201024666230821103601] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/27/2022] [Revised: 03/03/2023] [Accepted: 03/20/2023] [Indexed: 08/23/2023]
Abstract
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.
Collapse
Affiliation(s)
- Sona Rawat
- School of Life Sciences, Jaipur National University, Jaipur-302017, India
| | - Kanmani Subramaniam
- Department of Civil Engineering, KPR Institute of Engineering and Technology, Coimbatore-641407, Tamil Nadu, India
| | - Selva Kumar Subramanian
- Department of Sciences, Amrita School of Engineering, Coimbatore - 641112, Tamil Nadu, India
| | - Saravanan Subbarayan
- Department of Civil Engineering, National Institute of Technology, Trichy-620015, Tamil Nadu, India
| | - Subramanian Dhanabalan
- Department of Mechanical Engineering, M. Kumarasamy College of Engineering, Karur - 639113, Tamil Nadu, India
| | | | - Balasubramaniam Stalin
- Department of Mechanical Engineering, Anna University, Regional Campus Madurai, Madurai - 625 019, Tamil Nadu, India
| | - Arpita Roy
- Department of Biotechnology, School of Engineering & Technology, Sharda University, Greater Noida 201310, India
| | - Nagaraj Nagaprasad
- Department of Mechanical Engineering, ULTRA College of Engineering and Technology, Madurai - 625104, Tamilnadu, India
| | - Mahalingam Aruna
- College of Engineering and Computing, Al Ghurair University, Academic City, Dubai, UAE
| | - Jule Leta Tesfaye
- Dambi Dollo University, College of Natural and Computational Science, Department of Physics, Ethiopia
- Centre for Excellence-Indigenous Knowledge, Innovative Technology Transfer and Entrepreneurship, Dambi Dollo University, Dambi Dollo, Ethiopia
- Ministry of innovation and technology, Ethiopia
| | - Bayissa Badassa
- Department of Mechanical Engineering, Dambi Dollo University, Dambi Dollo, Ethiopia
| | - Ramaswamy Krishnaraj
- Centre for Excellence-Indigenous Knowledge, Innovative Technology Transfer and Entrepreneurship, Dambi Dollo University, Dambi Dollo, Ethiopia
- Ministry of innovation and technology, Ethiopia
- Department of Mechanical Engineering, Dambi Dollo University, Dambi Dollo, Ethiopia
| |
Collapse
|
2
|
Liu X, Duan Y, Hong X, Xie J, Liu S. Challenges in structural modeling of RNA-protein interactions. Curr Opin Struct Biol 2023; 81:102623. [PMID: 37301066 DOI: 10.1016/j.sbi.2023.102623] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/31/2023] [Revised: 05/14/2023] [Accepted: 05/16/2023] [Indexed: 06/12/2023]
Abstract
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.
Collapse
Affiliation(s)
- Xudong Liu
- School of Physics, Huazhong University of Science and Technology, Wuhan, Hubei, 430074, China
| | - Yingtian Duan
- School of Physics, Huazhong University of Science and Technology, Wuhan, Hubei, 430074, China
| | - Xu Hong
- School of Physics, Huazhong University of Science and Technology, Wuhan, Hubei, 430074, China
| | - Juan Xie
- School of Physics, Huazhong University of Science and Technology, Wuhan, Hubei, 430074, China
| | - Shiyong Liu
- School of Physics, Huazhong University of Science and Technology, Wuhan, Hubei, 430074, China.
| |
Collapse
|
3
|
Mollica L, Cupaioli FA, Rossetti G, Chiappori F. An overview of structural approaches to study therapeutic RNAs. Front Mol Biosci 2022; 9:1044126. [PMID: 36387283 PMCID: PMC9649582 DOI: 10.3389/fmolb.2022.1044126] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/14/2022] [Accepted: 10/18/2022] [Indexed: 11/07/2023] Open
Abstract
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.
Collapse
Affiliation(s)
- Luca Mollica
- Department of Medical Biotechnologies and Translational Medicine, L.I.T.A/University of Milan, Milan, Italy
| | | | | | - Federica Chiappori
- National Research Council—Institute for Biomedical Technologies, Milan, Italy
| |
Collapse
|
4
|
Bheemireddy S, Sandhya S, Srinivasan N, Sowdhamini R. Computational tools to study RNA-protein complexes. Front Mol Biosci 2022; 9:954926. [PMID: 36275618 PMCID: PMC9585174 DOI: 10.3389/fmolb.2022.954926] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/27/2022] [Accepted: 09/20/2022] [Indexed: 11/19/2022] Open
Abstract
RNA is the key player in many cellular processes such as signal transduction, replication, transport, cell division, transcription, and translation. These diverse functions are accomplished through interactions of RNA with proteins. However, protein–RNA interactions are still poorly derstood in contrast to protein–protein and protein–DNA interactions. This knowledge gap can be attributed to the limited availability of protein-RNA structures along with the experimental difficulties in studying these complexes. Recent progress in computational resources has expanded the number of tools available for studying protein-RNA interactions at various molecular levels. These include tools for predicting interacting residues from primary sequences, modelling of protein-RNA complexes, predicting hotspots in these complexes and insights into derstanding in the dynamics of their interactions. Each of these tools has its strengths and limitations, which makes it significant to select an optimal approach for the question of interest. Here we present a mini review of computational tools to study different aspects of protein-RNA interactions, with focus on overall application, development of the field and the future perspectives.
Collapse
Affiliation(s)
- Sneha Bheemireddy
- Molecular Biophysics Unit, Indian Institute of Science, Bangalore, India
| | - Sankaran Sandhya
- Department of Biotechnology, Faculty of Life and Allied Health Sciences, M.S. Ramaiah University of Applied Sciences, Bengaluru, India
- *Correspondence: Sankaran Sandhya, ; Ramanathan Sowdhamini,
| | | | - Ramanathan Sowdhamini
- Molecular Biophysics Unit, Indian Institute of Science, Bangalore, India
- National Centre for Biological Sciences, TIFR, GKVK Campus, Bangalore, India
- Institute of Bioinformatics and Applied Biotechnology, Bangalore, India
- *Correspondence: Sankaran Sandhya, ; Ramanathan Sowdhamini,
| |
Collapse
|
5
|
Kawai K, Asanuma Y, Kato T, Karuo Y, Tarui A, Sato K, Omote M. LCP: Simple Representation of Docking Poses for Machine Learning: A Case Study on Xanthine Oxidase Inhibitors. Mol Inform 2021; 41:e2100245. [PMID: 34843171 DOI: 10.1002/minf.202100245] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/08/2021] [Accepted: 11/21/2021] [Indexed: 11/05/2022]
Abstract
In this paper, we propose a simple descriptor called the ligand coordinate profile (LCP) for describing docking poses. The LCP descriptor is generated from the coordinates of the polar hydrogen and heavy atoms of the docked ligand. We hypothesize that the prediction of binding poses can be enhanced through the combination of machine learning methods with the LCP descriptor. Two docking programs were used to predict ligand docking against xanthine oxidase. Four machine learning methods-k-nearest neighbors, random forest, support vector machine, and LightGBM-were used to determine whether machine learning-based models could be used to accurately identify the correct binding poses. Regardless of the machine learning method employed, the LCP descriptor demonstrated improved performance compared to the existing descriptor. The results of the leave-one-pdb-out approach revealed that the influence of the pose descriptor was also significant, as demonstrated through cross-validation. When evaluated using top-N metrics, the machine learning models were generally more effective than the docking programs. In addition, the LCP-based models outperformed those based on the existing descriptor. The results obtained in this study suggest that our proposed binding pose descriptor is effective for improving the docking accuracy of xanthine oxidase inhibitors.
Collapse
Affiliation(s)
- Kentaro Kawai
- Faculty of Pharmaceutical Sciences, Setsunan University, 45-1, Nagaotoge-cho, Hirakata, Osaka, 573-0101, Japan
| | - Yoshitaka Asanuma
- Faculty of Pharmaceutical Sciences, Setsunan University, 45-1, Nagaotoge-cho, Hirakata, Osaka, 573-0101, Japan
| | - Toshiki Kato
- Faculty of Pharmaceutical Sciences, Setsunan University, 45-1, Nagaotoge-cho, Hirakata, Osaka, 573-0101, Japan
| | - Yukiko Karuo
- Faculty of Pharmaceutical Sciences, Setsunan University, 45-1, Nagaotoge-cho, Hirakata, Osaka, 573-0101, Japan
| | - Atsushi Tarui
- Faculty of Pharmaceutical Sciences, Setsunan University, 45-1, Nagaotoge-cho, Hirakata, Osaka, 573-0101, Japan
| | - Kazuyuki Sato
- Faculty of Pharmaceutical Sciences, Setsunan University, 45-1, Nagaotoge-cho, Hirakata, Osaka, 573-0101, Japan
| | - Masaaki Omote
- Faculty of Pharmaceutical Sciences, Setsunan University, 45-1, Nagaotoge-cho, Hirakata, Osaka, 573-0101, Japan
| |
Collapse
|
6
|
Blersch KF, Burchert JP, August SC, Welp L, Neumann P, Köster S, Urlaub H, Ficner R. Structural model of the M7G46 Methyltransferase TrmB in complex with tRNA. RNA Biol 2021; 18:2466-2479. [PMID: 34006170 DOI: 10.1080/15476286.2021.1925477] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022] Open
Abstract
TrmB belongs to the class I S-adenosylmethionine (SAM)-dependent methyltransferases (MTases) and introduces a methyl group to guanine at position 7 (m7G) in tRNA. In tRNAs m7G is most frequently found at position 46 in the variable loop and forms a tertiary base pair with C13 and U22, introducing a positive charge at G46. The TrmB/Trm8 enzyme family is structurally diverse, as TrmB proteins exist in a monomeric, homodimeric, and heterodimeric form. So far, the exact enzymatic mechanism, as well as the tRNA-TrmB crystal structure is not known. Here we present the first crystal structures of B. subtilis TrmB in complex with SAM and SAH. The crystal structures of TrmB apo and in complex with SAM and SAH have been determined by X-ray crystallography to 1.9 Å (apo), 2.5 Å (SAM), and 3.1 Å (SAH). The obtained crystal structures revealed Tyr193 to be important during SAM binding and MTase activity. Applying fluorescence polarization, the dissociation constant Kd of TrmB and tRNAPhe was determined to be 0.12 µM ± 0.002 µM. Luminescence-based methyltransferase activity assays revealed cooperative effects during TrmB catalysis with half-of-the-site reactivity at physiological SAM concentrations. Structural data retrieved from small-angle x-ray scattering (SAXS), mass-spectrometry of cross-linked complexes, and molecular docking experiments led to the determination of the TrmB-tRNAPhe complex structure.
Collapse
Affiliation(s)
- Katharina F Blersch
- Department of Molecular Structural Biology, Institute of Microbiology and Genetics, GZMB, Georg August University Göttingen, Göttingen, Germany
| | - Jan-Philipp Burchert
- Institute for X-Ray Physics, Georg August University Göttingen, Göttingen, Germany
| | | | - Luisa Welp
- Bioanalytical Mass Spectrometry Research Group, Max Planck Institute for Biophysical Chemistry, Göttingen, Germany
| | - Piotr Neumann
- Department of Molecular Structural Biology, Institute of Microbiology and Genetics, GZMB, Georg August University Göttingen, Göttingen, Germany
| | - Sarah Köster
- Institute for X-Ray Physics, Georg August University Göttingen, Göttingen, Germany
| | - Henning Urlaub
- Bioanalytical Mass Spectrometry Research Group, Max Planck Institute for Biophysical Chemistry, Göttingen, Germany.,Bioanalytics Group, Institute for Clinical Chemistry, University Medical Center Göttingen, Göttingen, Germany
| | - Ralf Ficner
- Department of Molecular Structural Biology, Institute of Microbiology and Genetics, GZMB, Georg August University Göttingen, Göttingen, Germany
| |
Collapse
|
7
|
Xia CQ, Pan X, Yang Y, Huang Y, Shen HB. Recent Progresses of Computational Analysis of RNA-Protein Interactions. SYSTEMS MEDICINE 2021. [DOI: 10.1016/b978-0-12-801238-3.11315-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022] Open
|
8
|
Zheng J, Hong X, Xie J, Tong X, Liu S. P3DOCK: a protein-RNA docking webserver based on template-based and template-free docking. Bioinformatics 2020; 36:96-103. [PMID: 31173056 DOI: 10.1093/bioinformatics/btz478] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/02/2019] [Revised: 05/24/2019] [Accepted: 06/04/2019] [Indexed: 01/02/2023] Open
Abstract
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.
Collapse
Affiliation(s)
- Jinfang Zheng
- School of Physics, Huazhong University of Science and Technology, Wuhan, Hubei 430074, China
| | - Xu Hong
- School of Physics, Huazhong University of Science and Technology, Wuhan, Hubei 430074, China
| | - Juan Xie
- School of Physics, Huazhong University of Science and Technology, Wuhan, Hubei 430074, China
| | - Xiaoxue Tong
- School of Physics, Huazhong University of Science and Technology, Wuhan, Hubei 430074, China
| | - Shiyong Liu
- School of Physics, Huazhong University of Science and Technology, Wuhan, Hubei 430074, China
| |
Collapse
|
9
|
He J, Tao H, Huang SY. Protein-ensemble-RNA docking by efficient consideration of protein flexibility through homology models. Bioinformatics 2020; 35:4994-5002. [PMID: 31086984 DOI: 10.1093/bioinformatics/btz388] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/01/2019] [Revised: 04/28/2019] [Accepted: 05/03/2019] [Indexed: 12/18/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Jiahua He
- Institute of Biophysics, School of Physics, Huazhong University of Science and Technology, Wuhan, Hubei, China
| | - Huanyu Tao
- Institute of Biophysics, School of Physics, Huazhong University of Science and Technology, Wuhan, Hubei, China
| | - Sheng-You Huang
- Institute of Biophysics, School of Physics, Huazhong University of Science and Technology, Wuhan, Hubei, China
| |
Collapse
|
10
|
Integrative Structural Biology of Protein-RNA Complexes. Structure 2020; 28:6-28. [DOI: 10.1016/j.str.2019.11.017] [Citation(s) in RCA: 23] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/05/2019] [Revised: 11/17/2019] [Accepted: 11/27/2019] [Indexed: 12/16/2022]
|
11
|
Deng L, Yang W, Liu H. PredPRBA: Prediction of Protein-RNA Binding Affinity Using Gradient Boosted Regression Trees. Front Genet 2019; 10:637. [PMID: 31428122 PMCID: PMC6688581 DOI: 10.3389/fgene.2019.00637] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/03/2019] [Accepted: 06/18/2019] [Indexed: 01/24/2023] Open
Abstract
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/.
Collapse
Affiliation(s)
- Lei Deng
- School of Computer Science and Engineering, Central South University, Changsha, China.,School of Software, Xinjiang University, Urumqi, China
| | - Wenyi Yang
- School of Computer Science and Engineering, Central South University, Changsha, China
| | - Hui Liu
- Lab of Information Management, Changzhou University, Changzhou, China
| |
Collapse
|
12
|
Pal A, Levy Y. Structure, stability and specificity of the binding of ssDNA and ssRNA with proteins. PLoS Comput Biol 2019; 15:e1006768. [PMID: 30933978 PMCID: PMC6467422 DOI: 10.1371/journal.pcbi.1006768] [Citation(s) in RCA: 37] [Impact Index Per Article: 7.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/13/2018] [Revised: 04/16/2019] [Accepted: 01/01/2019] [Indexed: 02/06/2023] Open
Abstract
Recognition of single-stranded DNA (ssDNA) or single-stranded RNA (ssRNA) is important for many fundamental cellular functions. A variety of single-stranded DNA-binding proteins (ssDBPs) and single-stranded RNA-binding proteins (ssRBPs) have evolved that bind ssDNA and ssRNA, respectively, with varying degree of affinities and specificities to form complexes. Structural studies of these complexes provide key insights into their recognition mechanism. However, computational modeling of the specific recognition process and to predict the structure of the complex is challenging, primarily due to the heterogeneity of their binding energy landscape and the greater flexibility of ssDNA or ssRNA compared with double-stranded nucleic acids. Consequently, considerably fewer computational studies have explored interactions between proteins and single-stranded nucleic acids compared with protein interactions with double-stranded nucleic acids. Here, we report a newly developed energy-based coarse-grained model to predict the structure of ssDNA–ssDBP and ssRNA–ssRBP complexes and to assess their sequence-specific interactions and stabilities. We tuned two factors that can modulate specific recognition: base–aromatic stacking strength and the flexibility of the single-stranded nucleic acid. The model was successfully applied to predict the binding conformations of 12 distinct ssDBP and ssRBP structures with their cognate ssDNA and ssRNA partners having various sequences. Estimated binding energies agreed well with the corresponding experimental binding affinities. Bound conformations from the simulation showed a funnel-shaped binding energy distribution where the native-like conformations corresponded to the energy minima. The various ssDNA–protein and ssRNA–protein complexes differed in the balance of electrostatic and aromatic energies. The lower affinity of the ssRNA–ssRBP complexes compared with the ssDNA–ssDBP complexes stems from lower flexibility of ssRNA compared to ssDNA, which results in higher rate constants for the dissociation of the complex (koff) for complexes involving the former. Quantifying bimolecular self-assembly is pivotal to understanding cellular function. In recent years, a large progress has been made in understanding the structure and biophysics of protein-protein interactions. Particularly, various computational tools are available for predicting these structures and to estimate their stability and the driving forces of their formation. The understating of the interactions between proteins and nucleic acids, however, is still limited, presumably due to the involvement of non-specific interactions as well as the high conformational plasticity that may demand an induced-fit mechanism. In particular, the interactions between proteins and single-stranded nucleic acids (i.e., single-stranded DNA and RNA) is very challenging due to their high flexibility. Furthermore, the interface between proteins and single-stranded nucleic acids is often chemically more heterogeneous than the interface between proteins and double-stranded DNA. In this study, we developed a coarse-grained computational model to predict the structure of complexes between proteins and single-stranded nucleic acids. The model was applied to estimate binding affinities and the estimated binding energies agreed well with the corresponding experimental binding affinities. The kinetics of association as well as the specificity of the complexes between proteins and ssDNA are different than those with ssRNA, mostly due to differences in their conformational flexibility.
Collapse
Affiliation(s)
- Arumay Pal
- Department of Structural Biology, Weizmann Institute of Science, Rehovot, Israel
| | - Yaakov Levy
- Department of Structural Biology, Weizmann Institute of Science, Rehovot, Israel
- * E-mail:
| |
Collapse
|
13
|
Kappel K, Das R. Sampling Native-like Structures of RNA-Protein Complexes through Rosetta Folding and Docking. Structure 2018; 27:140-151.e5. [PMID: 30416038 DOI: 10.1016/j.str.2018.10.001] [Citation(s) in RCA: 26] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/07/2018] [Revised: 08/27/2018] [Accepted: 10/05/2018] [Indexed: 10/27/2022]
Abstract
RNA-protein complexes underlie numerous cellular processes including translation, splicing, and posttranscriptional regulation of gene expression. The structures of these complexes are crucial to their functions but often elude high-resolution structure determination. Computational methods are needed that can integrate low-resolution data for RNA-protein complexes while modeling de novo the large conformational changes of RNA components upon complex formation. To address this challenge, we describe RNP-denovo, a Rosetta method to simultaneously fold-and-dock RNA to a protein surface. On a benchmark set of diverse RNA-protein complexes not solvable with prior strategies, RNP-denovo consistently sampled native-like structures with better than nucleotide resolution. We revisited three past blind modeling challenges involving the spliceosome, telomerase, and a methyltransferase-ribosomal RNA complex in which previous methods gave poor results. When coupled with the same sparse FRET, crosslinking, and functional data used previously, RNP-denovo gave models with significantly improved accuracy. These results open a route to modeling global folds of RNA-protein complexes from low-resolution data.
Collapse
Affiliation(s)
- Kalli Kappel
- Biophysics Program, Stanford University, Stanford, CA 94305, USA
| | - Rhiju Das
- Biophysics Program, Stanford University, Stanford, CA 94305, USA; Department of Biochemistry, Stanford University School of Medicine, Stanford, CA 94305, USA; Department of Physics, Stanford University, Stanford, CA 94305, USA.
| |
Collapse
|
14
|
Chen F, Sun H, Wang J, Zhu F, Liu H, Wang Z, Lei T, Li Y, Hou T. Assessing the performance of MM/PBSA and MM/GBSA methods. 8. Predicting binding free energies and poses of protein-RNA complexes. RNA (NEW YORK, N.Y.) 2018; 24:1183-1194. [PMID: 29930024 PMCID: PMC6097651 DOI: 10.1261/rna.065896.118] [Citation(s) in RCA: 67] [Impact Index Per Article: 11.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/29/2018] [Accepted: 06/13/2018] [Indexed: 05/10/2023]
Abstract
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 GBGBn1 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 GBGBn1 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.
Collapse
Affiliation(s)
- Fu Chen
- College of Pharmaceutical Sciences, Zhejiang University, Hangzhou, Zhejiang 310058, China
- State Key Lab of CAD&CG, Zhejiang University, Hangzhou, Zhejiang 310058, China
- College of Life and Environmental Sciences, Shanghai Normal University, Shanghai 200234, China
| | - Huiyong Sun
- College of Pharmaceutical Sciences, Zhejiang University, Hangzhou, Zhejiang 310058, China
| | - Junmei Wang
- Department of Pharmaceutical Sciences, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, USA
| | - Feng Zhu
- College of Pharmaceutical Sciences, Zhejiang University, Hangzhou, Zhejiang 310058, China
| | - Hui Liu
- College of Pharmaceutical Sciences, Zhejiang University, Hangzhou, Zhejiang 310058, China
| | - Zhe Wang
- College of Pharmaceutical Sciences, Zhejiang University, Hangzhou, Zhejiang 310058, China
- State Key Lab of CAD&CG, Zhejiang University, Hangzhou, Zhejiang 310058, China
| | - Tailong Lei
- College of Pharmaceutical Sciences, Zhejiang University, Hangzhou, Zhejiang 310058, China
| | - Youyong Li
- Institute of Functional Nano and Soft Materials (FUNSOM), Soochow University, Suzhou, Jiangsu 215123, China
| | - Tingjun Hou
- College of Pharmaceutical Sciences, Zhejiang University, Hangzhou, Zhejiang 310058, China
- State Key Lab of CAD&CG, Zhejiang University, Hangzhou, Zhejiang 310058, China
| |
Collapse
|
15
|
Nithin C, Ghosh P, Bujnicki JM. Bioinformatics Tools and Benchmarks for Computational Docking and 3D Structure Prediction of RNA-Protein Complexes. Genes (Basel) 2018; 9:genes9090432. [PMID: 30149645 PMCID: PMC6162694 DOI: 10.3390/genes9090432] [Citation(s) in RCA: 26] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/14/2018] [Revised: 07/26/2018] [Accepted: 08/21/2018] [Indexed: 12/29/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Chandran Nithin
- Laboratory of Bioinformatics and Protein Engineering, International Institute of Molecular and Cell Biology in Warsaw, ul. Ks. Trojdena 4, PL-02-109 Warsaw, Poland.
| | - Pritha Ghosh
- Laboratory of Bioinformatics and Protein Engineering, International Institute of Molecular and Cell Biology in Warsaw, ul. Ks. Trojdena 4, PL-02-109 Warsaw, Poland.
| | - Janusz M Bujnicki
- Laboratory of Bioinformatics and Protein Engineering, International Institute of Molecular and Cell Biology in Warsaw, ul. Ks. Trojdena 4, PL-02-109 Warsaw, Poland.
- Bioinformatics Laboratory, Institute of Molecular Biology and Biotechnology, Faculty of Biology, Adam Mickiewicz University, ul. Umultowska 89, PL-61-614 Poznan, Poland.
| |
Collapse
|
16
|
Abstract
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.
Collapse
Affiliation(s)
- Yelena A Arnautova
- Molsoft L.L.C., 11199 Sorrento Valley Road, S209 , San Diego , California 92121 , United States
| | - Ruben Abagyan
- Skaggs School of Pharmacy and Pharmaceutical Sciences , University of California San Diego , La Jolla , California 92093 , United States
| | - Maxim Totrov
- Molsoft L.L.C., 11199 Sorrento Valley Road, S209 , San Diego , California 92121 , United States
| |
Collapse
|
17
|
Structural insights into the stimulation of S. pombe Dnmt2 catalytic efficiency by the tRNA nucleoside queuosine. Sci Rep 2018; 8:8880. [PMID: 29892076 PMCID: PMC5995894 DOI: 10.1038/s41598-018-27118-5] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/22/2018] [Accepted: 05/24/2018] [Indexed: 01/16/2023] Open
Abstract
Dnmt2 methylates cytosine at position 38 of tRNAAsp in a variety of eukaryotic organisms. A correlation between the presence of the hypermodified nucleoside queuosine (Q) at position 34 of tRNAAsp and the Dnmt2 dependent C38 methylation was recently found in vivo for S. pombe and D. discoideum. We demonstrate a direct effect of the Q-modification on the methyltransferase catalytic efficiency in vitro, as Vmax/K0.5 of purified S. pombe Dnmt2 shows an increase for in vitro transcribed tRNAAsp containing Q34 to 6.27 ∗ 10–3 s−1 µM−1 compared to 1.51 ∗ 10–3 s−1 µM−1 for the unmodified substrate. Q34tRNAAsp exhibits an only slightly increased affinity for Dnmt2 in comparison to unmodified G34tRNA. In order to get insight into the structural basis for the Q-dependency, the crystal structure of S. pombe Dnmt2 was determined at 1.7 Å resolution. It closely resembles the known structures of human and E. histolytica Dnmt2, and contains the entire active site loop. The interaction with tRNA was analyzed by means of mass-spectrometry using UV cross-linked Dnmt2-tRNA complex. These cross-link data and computational docking of Dnmt2 and tRNAAsp reveal Q34 positioned adjacent to the S-adenosylmethionine occupying the active site, suggesting that the observed increase of Dnmt2 catalytic efficiency by queuine originates from optimal positioning of the substrate molecules and residues relevant for methyl transfer.
Collapse
|
18
|
A pair-conformation-dependent scoring function for evaluating 3D RNA-protein complex structures. PLoS One 2017; 12:e0174662. [PMID: 28358834 PMCID: PMC5373608 DOI: 10.1371/journal.pone.0174662] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/12/2016] [Accepted: 03/13/2017] [Indexed: 01/04/2023] Open
Abstract
Computational prediction of RNA-protein complex 3D structures includes two basic steps: one is sampling possible structures and another is scoring the sampled structures to pick out the correct one. At present, constructing accurate scoring functions is still not well solved and the performances of the scoring functions usually depend on used benchmarks. Here we propose a pair-conformation-dependent scoring function, 3dRPC-Score, for 3D RNA-protein complex structure prediction by considering the nucleotide-residue pairs having the same energy if their conformations are similar, instead of the distance-only dependence of the most existing scoring functions. Benchmarking shows that 3dRPC-Score has a consistent performance in three test sets.
Collapse
|
19
|
Elucidating Mechanisms of Molecular Recognition Between Human Argonaute and miRNA Using Computational Approaches. Methods Mol Biol 2016. [PMID: 27924488 DOI: 10.1007/978-1-4939-6563-2_18] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/05/2024]
Abstract
MicroRNA (miRNA) and Argonaute (AGO) protein together form the RNA-induced silencing complex (RISC) that plays an essential role in the regulation of gene expression. Elucidating the underlying mechanism of AGO-miRNA recognition is thus of great importance not only for the in-depth understanding of miRNA function but also for inspiring new drugs targeting miRNAs. In this chapter we introduce a combined computational approach of molecular dynamics (MD) simulations, Markov state models (MSMs), and protein-RNA docking to investigate AGO-miRNA recognition. Constructed from MD simulations, MSMs can elucidate the conformational dynamics of AGO at biologically relevant timescales. Protein-RNA docking can then efficiently identify the AGO conformations that are geometrically accessible to miRNA. Using our recent work on human AGO2 as an example, we explain the rationale and the workflow of our method in details. This combined approach holds great promise to complement experiments in unraveling the mechanisms of molecular recognition between large, flexible, and complex biomolecules.
Collapse
|
20
|
Nithin C, Mukherjee S, Bahadur RP. A non-redundant protein-RNA docking benchmark version 2.0. Proteins 2016; 85:256-267. [PMID: 27862282 DOI: 10.1002/prot.25211] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/08/2016] [Revised: 10/27/2016] [Accepted: 11/08/2016] [Indexed: 12/23/2022]
Abstract
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.
Collapse
Affiliation(s)
- Chandran Nithin
- Computational Structural Biology Lab, Department of Biotechnology, Indian Institute of Technology Kharagpur, 721302, India
| | - Sunandan Mukherjee
- Computational Structural Biology Lab, Department of Biotechnology, Indian Institute of Technology Kharagpur, 721302, India
| | - Ranjit Prasad Bahadur
- Computational Structural Biology Lab, Department of Biotechnology, Indian Institute of Technology Kharagpur, 721302, India
| |
Collapse
|
21
|
Pérez-Cano L, Romero-Durana M, Fernández-Recio J. Structural and energy determinants in protein-RNA docking. Methods 2016; 118-119:163-170. [PMID: 27816523 DOI: 10.1016/j.ymeth.2016.11.001] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/20/2016] [Revised: 10/14/2016] [Accepted: 11/01/2016] [Indexed: 01/02/2023] Open
Abstract
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.
Collapse
Affiliation(s)
- Laura Pérez-Cano
- Joint BSC-CRG-IRB Research Program in Computational Biology, Life Sciences Department, Barcelona Supercomputing Center (BSC), Jordi Girona 29, Barcelona 08034, Spain; Center for Neurobehavioral Genetics and Center for Autism Research and Treatment, Semel Institute, David Geffen School of Medicine, University of California, Los Angeles, CA, USA
| | - Miguel Romero-Durana
- Joint BSC-CRG-IRB Research Program in Computational Biology, Life Sciences Department, Barcelona Supercomputing Center (BSC), Jordi Girona 29, Barcelona 08034, Spain
| | - Juan Fernández-Recio
- Joint BSC-CRG-IRB Research Program in Computational Biology, Life Sciences Department, Barcelona Supercomputing Center (BSC), Jordi Girona 29, Barcelona 08034, Spain.
| |
Collapse
|
22
|
Zheng J, Kundrotas PJ, Vakser IA, Liu S. Template-Based Modeling of Protein-RNA Interactions. PLoS Comput Biol 2016; 12:e1005120. [PMID: 27662342 PMCID: PMC5035060 DOI: 10.1371/journal.pcbi.1005120] [Citation(s) in RCA: 25] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/25/2016] [Accepted: 08/25/2016] [Indexed: 12/29/2022] Open
Abstract
Protein-RNA complexes formed by specific recognition between RNA and RNA-binding proteins play an important role in biological processes. More than a thousand of such proteins in human are curated and many novel RNA-binding proteins are to be discovered. Due to limitations of experimental approaches, computational techniques are needed for characterization of protein-RNA interactions. Although much progress has been made, adequate methodologies reliably providing atomic resolution structural details are still lacking. Although protein-RNA free docking approaches proved to be useful, in general, the template-based approaches provide higher quality of predictions. Templates are key to building a high quality model. Sequence/structure relationships were studied based on a representative set of binary protein-RNA complexes from PDB. Several approaches were tested for pairwise target/template alignment. The analysis revealed a transition point between random and correct binding modes. The results showed that structural alignment is better than sequence alignment in identifying good templates, suitable for generating protein-RNA complexes close to the native structure, and outperforms free docking, successfully predicting complexes where the free docking fails, including cases of significant conformational change upon binding. A template-based protein-RNA interaction modeling protocol PRIME was developed and benchmarked on a representative set of complexes. Structures of protein-RNA complexes are important for characterization of biological processes. The number of experimentally determined protein-RNA complexes is limited. Thus modeling of these complexes is important. Reliable structural predictions of proteins and their complexes are provided by comparative modeling, which takes advantage of similar complexes with experimentally determined structures. Thus, in the case of protein-RNA complexes, it is important to determine if similar proteins and RNAs bind in a similar way. We show that, similarly to the earlier published results on protein-protein complexes, such correlation of the protein-RNA binding mode and the monomers similarity indeed exists, and is stronger when the similarity is determined by structure rather than sequence alignment. The data shows clear transition from random to similar binding mode with the increase of the structural similarity of the monomers. On the basis of the results we designed and implemented a predictive tool, which should be useful for the biological community interested in modeling of protein-RNA interactions.
Collapse
Affiliation(s)
- Jinfang Zheng
- School of Physics and Key Laboratory of Molecular Biophysics of the Ministry of Education, Huazhong University of Science and Technology, Wuhan, Hubei, China
| | - Petras J. Kundrotas
- Center for Computational Biology and Department of Molecular Biosciences, The University of Kansas, Lawrence, Kansas, United States of America
| | - Ilya A. Vakser
- Center for Computational Biology and Department of Molecular Biosciences, The University of Kansas, Lawrence, Kansas, United States of America
- * E-mail: (IAV); (SL)
| | - Shiyong Liu
- School of Physics and Key Laboratory of Molecular Biophysics of the Ministry of Education, Huazhong University of Science and Technology, Wuhan, Hubei, China
- * E-mail: (IAV); (SL)
| |
Collapse
|
23
|
Iwakiri J, Hamada M, Asai K, Kameda T. Improved Accuracy in RNA-Protein Rigid Body Docking by Incorporating Force Field for Molecular Dynamics Simulation into the Scoring Function. J Chem Theory Comput 2016; 12:4688-97. [PMID: 27494732 DOI: 10.1021/acs.jctc.6b00254] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
RNA-protein interactions play fundamental roles in many biological processes. To understand these interactions, it is necessary to know the three-dimensional structures of RNA-protein complexes. However, determining the tertiary structure of these complexes is often difficult, suggesting that an accurate rigid body docking for RNA-protein complexes is needed. In general, the rigid body docking process is divided into two steps: generating candidate structures from the individual RNA and protein structures and then narrowing down the candidates. In this study, we focus on the former problem to improve the prediction accuracy in RNA-protein docking. Our method is based on the integration of physicochemical information about RNA into ZDOCK, which is known as one of the most successful computer programs for protein-protein docking. Because recent studies showed the current force field for molecular dynamics simulation of protein and nucleic acids is quite accurate, we modeled the physicochemical information about RNA by force fields such as AMBER and CHARMM. A comprehensive benchmark of RNA-protein docking, using three recently developed data sets, reveals the remarkable prediction accuracy of the proposed method compared with existing programs for docking: the highest success rate is 34.7% for the predicted structure of the RNA-protein complex with the best score and 79.2% for 3,600 predicted ones. Three full atomistic force fields for RNA (AMBER94, AMBER99, and CHARMM22) produced almost the same accurate result, which showed current force fields for nucleic acids are quite accurate. In addition, we found that the electrostatic interaction and the representation of shape complementary between protein and RNA plays the important roles for accurate prediction of the native structures of RNA-protein complexes.
Collapse
Affiliation(s)
- Junichi Iwakiri
- Graduate School of Frontier Sciences, The University of Tokyo , 5-1-5 Kashiwanoha, Kashiwa, Chiba 277-8562, Japan
| | - Michiaki Hamada
- Department of Electrical Engineering and Bioscience, Faculty of Science and Engineering, Waseda University , 55N-06-10, 3-4-1, Okubo, Shinjuku-ku, Tokyo 169-8555, Japan.,Artificial Intelligence Research Center, National Institute of Advanced Industrial Science and Technology (AIST) , 2-4-7 Aomi, Koto-ku, Tokyo 135-0064, Japan
| | - Kiyoshi Asai
- Graduate School of Frontier Sciences, The University of Tokyo , 5-1-5 Kashiwanoha, Kashiwa, Chiba 277-8562, Japan.,Artificial Intelligence Research Center, National Institute of Advanced Industrial Science and Technology (AIST) , 2-4-7 Aomi, Koto-ku, Tokyo 135-0064, Japan
| | - Tomoshi Kameda
- Artificial Intelligence Research Center, National Institute of Advanced Industrial Science and Technology (AIST) , 2-4-7 Aomi, Koto-ku, Tokyo 135-0064, Japan
| |
Collapse
|
24
|
Bender BJ, Cisneros A, Duran AM, Finn JA, Fu D, Lokits AD, Mueller BK, Sangha AK, Sauer MF, Sevy AM, Sliwoski G, Sheehan JH, DiMaio F, Meiler J, Moretti R. Protocols for Molecular Modeling with Rosetta3 and RosettaScripts. Biochemistry 2016; 55:4748-63. [PMID: 27490953 PMCID: PMC5007558 DOI: 10.1021/acs.biochem.6b00444] [Citation(s) in RCA: 144] [Impact Index Per Article: 18.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
Abstract
![]()
Previously, we published an article
providing an overview of the
Rosetta suite of biomacromolecular modeling software and a series
of step-by-step tutorials [Kaufmann, K. W., et al. (2010) Biochemistry 49, 2987–2998]. The overwhelming positive
response to this publication we received motivates us to here share
the next iteration of these tutorials that feature de novo folding, comparative modeling, loop construction, protein docking,
small molecule docking, and protein design. This updated and expanded
set of tutorials is needed, as since 2010 Rosetta has been fully redesigned
into an object-oriented protein modeling program Rosetta3. Notable
improvements include a substantially improved energy function, an
XML-like language termed “RosettaScripts” for flexibly
specifying modeling task, new analysis tools, the addition of the
TopologyBroker to control conformational sampling, and support for
multiple templates in comparative modeling. Rosetta’s ability
to model systems with symmetric proteins, membrane proteins, noncanonical
amino acids, and RNA has also been greatly expanded and improved.
Collapse
Affiliation(s)
- Brian J Bender
- Department of Pharmacology, Vanderbilt University , Nashville, Tennessee 37232-6600, United States.,Center for Structural Biology, Vanderbilt University , Nashville, Tennessee 37240-7917, United States
| | - Alberto Cisneros
- Center for Structural Biology, Vanderbilt University , Nashville, Tennessee 37240-7917, United States.,Chemical and Physical Biology Program, Vanderbilt University , Nashville, Tennessee 37232-0301, United States
| | - Amanda M Duran
- Center for Structural Biology, Vanderbilt University , Nashville, Tennessee 37240-7917, United States.,Department of Chemistry, Vanderbilt University , Nashville, Tennessee 37235, United States
| | - Jessica A Finn
- Center for Structural Biology, Vanderbilt University , Nashville, Tennessee 37240-7917, United States.,Department of Pathology, Microbiology and Immunology, Vanderbilt University , Nashville, Tennessee 37232-2561, United States
| | - Darwin Fu
- Center for Structural Biology, Vanderbilt University , Nashville, Tennessee 37240-7917, United States.,Department of Chemistry, Vanderbilt University , Nashville, Tennessee 37235, United States
| | - Alyssa D Lokits
- Center for Structural Biology, Vanderbilt University , Nashville, Tennessee 37240-7917, United States.,Neuroscience Program, Vanderbilt University , Nashville, Tennessee 37235, United States
| | - Benjamin K Mueller
- Center for Structural Biology, Vanderbilt University , Nashville, Tennessee 37240-7917, United States.,Department of Chemistry, Vanderbilt University , Nashville, Tennessee 37235, United States
| | - Amandeep K Sangha
- Center for Structural Biology, Vanderbilt University , Nashville, Tennessee 37240-7917, United States.,Department of Chemistry, Vanderbilt University , Nashville, Tennessee 37235, United States
| | - Marion F Sauer
- Center for Structural Biology, Vanderbilt University , Nashville, Tennessee 37240-7917, United States.,Chemical and Physical Biology Program, Vanderbilt University , Nashville, Tennessee 37232-0301, United States
| | - Alexander M Sevy
- Center for Structural Biology, Vanderbilt University , Nashville, Tennessee 37240-7917, United States.,Chemical and Physical Biology Program, Vanderbilt University , Nashville, Tennessee 37232-0301, United States
| | - Gregory Sliwoski
- Center for Structural Biology, Vanderbilt University , Nashville, Tennessee 37240-7917, United States.,Department of Chemistry, Vanderbilt University , Nashville, Tennessee 37235, United States
| | - Jonathan H Sheehan
- Center for Structural Biology, Vanderbilt University , Nashville, Tennessee 37240-7917, United States
| | - Frank DiMaio
- Department of Biochemistry, University of Washington , Seattle, Washington 98195, United States
| | - Jens Meiler
- Department of Pharmacology, Vanderbilt University , Nashville, Tennessee 37232-6600, United States.,Center for Structural Biology, Vanderbilt University , Nashville, Tennessee 37240-7917, United States.,Chemical and Physical Biology Program, Vanderbilt University , Nashville, Tennessee 37232-0301, United States.,Department of Chemistry, Vanderbilt University , Nashville, Tennessee 37235, United States.,Department of Pathology, Microbiology and Immunology, Vanderbilt University , Nashville, Tennessee 37232-2561, United States.,Neuroscience Program, Vanderbilt University , Nashville, Tennessee 37235, United States
| | - Rocco Moretti
- Center for Structural Biology, Vanderbilt University , Nashville, Tennessee 37240-7917, United States.,Department of Chemistry, Vanderbilt University , Nashville, Tennessee 37235, United States
| |
Collapse
|
25
|
In Silico Aptamer Docking Studies: From a Retrospective Validation to a Prospective Case Study'TIM3 Aptamers Binding. MOLECULAR THERAPY-NUCLEIC ACIDS 2016; 5:e376. [DOI: 10.1038/mtna.2016.84] [Citation(s) in RCA: 26] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/08/2016] [Accepted: 08/22/2016] [Indexed: 12/25/2022]
|
26
|
Dawson WK, Bujnicki JM. Computational modeling of RNA 3D structures and interactions. Curr Opin Struct Biol 2015; 37:22-8. [PMID: 26689764 DOI: 10.1016/j.sbi.2015.11.007] [Citation(s) in RCA: 39] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/01/2015] [Revised: 11/11/2015] [Accepted: 11/12/2015] [Indexed: 11/25/2022]
Abstract
RNA molecules have key functions in cellular processes beyond being carriers of protein-coding information. These functions are often dependent on the ability to form complex three-dimensional (3D) structures. However, experimental determination of RNA 3D structures is difficult, which has prompted the development of computational methods for structure prediction from sequence. Recent progress in 3D structure modeling of RNA and emerging approaches for predicting RNA interactions with ions, ligands and proteins have been stimulated by successes in protein 3D structure modeling.
Collapse
Affiliation(s)
- Wayne K Dawson
- Laboratory of Bioinformatics and Protein Engineering, International Institute of Molecular and Cell Biology, ul. Ks. Trojdena 4, 02-109 Warsaw, Poland
| | - Janusz M Bujnicki
- Laboratory of Bioinformatics and Protein Engineering, International Institute of Molecular and Cell Biology, ul. Ks. Trojdena 4, 02-109 Warsaw, Poland; Bioinformatics Laboratory, Institute of Molecular Biology and Biotechnology, Faculty of Biology, Adam Mickiewicz University, ul. Umultowska 89, 61-614 Poznan, Poland.
| |
Collapse
|
27
|
Nahalka J, Hrabarova E, Talafova K. Protein-RNA and protein-glycan recognitions in light of amino acid codes. Biochim Biophys Acta Gen Subj 2015; 1850:1942-52. [PMID: 26145579 DOI: 10.1016/j.bbagen.2015.06.013] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/17/2015] [Revised: 06/18/2015] [Accepted: 06/22/2015] [Indexed: 12/13/2022]
Abstract
BACKGROUND RNA-binding proteins, in cooperation with non-coding RNAs, play important roles in post-transcriptional regulation. Non-coding micro-RNAs control information flow from the genome to the glycome by interacting with glycan-synthesis enzymes. Glycan-binding proteins read the cell surface and cytoplasmic glycome and transfer signals back to the nucleus. The profiling of the protein-RNA and protein-glycan interactomes is of significant medicinal importance. SCOPE OF REVIEW This review discusses the state-of-the-art research in the protein-RNA and protein-glycan recognition fields and proposes the application of amino acid codes in profiling and programming the interactomes. MAJOR CONCLUSIONS The deciphered PUF-RNA and PPR-RNA amino acid recognition codes can be explained by the protein-RNA amino acid recognition hypothesis based on the genetic code. The tripartite amino acid code is also involved in protein-glycan interactions. At present, the results indicate that a system of four codons ("gnc", where n=g - guanine, c - cytosine, u - uracil or a - adenine) and four amino acids (G - glycine, A - alanine, V - valine, D - aspartic acid) could be the original genetic code that imprinted "rules" into both recognition processes. GENERAL SIGNIFICANCE Amino acid recognition codes have provocative potential in the profiling and programming of the protein-RNA and protein-glycan interactomes. The profiling and even programming of the interactomes will play significant roles in diagnostics and the development of therapeutic procedures against cancer and neurodegenerative, developmental and other diseases.
Collapse
Affiliation(s)
- Jozef Nahalka
- Institute of Chemistry, Centre for Glycomics, Slovak Academy of Sciences, Dubravska cesta 9, SK-84538 Bratislava, Slovak Republic; Institute of Chemistry, Centre of Excellence for White-green Biotechnology, Slovak Academy of Sciences, Trieda Andreja Hlinku 2, SK-94976 Nitra, Slovak Republic.
| | - Eva Hrabarova
- Institute of Chemistry, Centre for Glycomics, Slovak Academy of Sciences, Dubravska cesta 9, SK-84538 Bratislava, Slovak Republic; Institute of Chemistry, Centre of Excellence for White-green Biotechnology, Slovak Academy of Sciences, Trieda Andreja Hlinku 2, SK-94976 Nitra, Slovak Republic
| | - Klaudia Talafova
- Institute of Chemistry, Centre for Glycomics, Slovak Academy of Sciences, Dubravska cesta 9, SK-84538 Bratislava, Slovak Republic; Institute of Chemistry, Centre of Excellence for White-green Biotechnology, Slovak Academy of Sciences, Trieda Andreja Hlinku 2, SK-94976 Nitra, Slovak Republic
| |
Collapse
|