1
|
Li R, Zhu Q, Sun X, Li Z, Liu X. Electrochemical biosensor based on the integration of maple leaf-like gold nanocrystal and truncated aptamer for detection of α-amanitin with high sensitivity, selectivity and rapidity. Food Chem 2024; 453:139639. [PMID: 38759442 DOI: 10.1016/j.foodchem.2024.139639] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/09/2024] [Revised: 04/01/2024] [Accepted: 05/09/2024] [Indexed: 05/19/2024]
Abstract
This study reports the fabrication of three-dimensional gold nanocrystals as sensing material in the presence of l-glutathion and high-performance aptamer with 20 bases of α-amanitin via truncation and optimization of along aptamer. The resulting maple leaf-like gold nanocrystal (ML-Au) exhibits an improved catalytic activity due to more exposed high-index facets. The use of truncated aptamer increases the sensitivity by 15 times and reduces the reaction time by two times compared with those of original aptamer. An α-amanitin electrochemical biosensor constructed by integrating ML-Au nanocrystals with truncated aptamer exhibits high sensitivity, selectivity and rapidity. An increase of the α-amanitin concentration in the range of 1 × 10-14-1 × 10-9 M causes a linear decrease in the amperometric current with a limit of detection of 2.9 × 10-15 M (S/N = 3). The proposed analytical method is satisfactorily used for electrochemical sensing of α-amanitin in urine and wild mushroom samples.
Collapse
Affiliation(s)
- Ruiyi Li
- Key Laboratory of Synthetic and Biological Colloids, Ministry of Education, School of Chemical and Material Engineering, School of Life Sciences and Health Engineering, and School of Food Science and Technology, Jiangnan University, Wuxi 214122, China
| | - Qiyue Zhu
- Key Laboratory of Synthetic and Biological Colloids, Ministry of Education, School of Chemical and Material Engineering, School of Life Sciences and Health Engineering, and School of Food Science and Technology, Jiangnan University, Wuxi 214122, China
| | - Xiulan Sun
- Key Laboratory of Synthetic and Biological Colloids, Ministry of Education, School of Chemical and Material Engineering, School of Life Sciences and Health Engineering, and School of Food Science and Technology, Jiangnan University, Wuxi 214122, China
| | - Zaijun Li
- Key Laboratory of Synthetic and Biological Colloids, Ministry of Education, School of Chemical and Material Engineering, School of Life Sciences and Health Engineering, and School of Food Science and Technology, Jiangnan University, Wuxi 214122, China.
| | - Xiaohao Liu
- Key Laboratory of Synthetic and Biological Colloids, Ministry of Education, School of Chemical and Material Engineering, School of Life Sciences and Health Engineering, and School of Food Science and Technology, Jiangnan University, Wuxi 214122, China
| |
Collapse
|
2
|
Nithin C, Kmiecik S, Błaszczyk R, Nowicka J, Tuszyńska I. Comparative analysis of RNA 3D structure prediction methods: towards enhanced modeling of RNA-ligand interactions. Nucleic Acids Res 2024:gkae541. [PMID: 38917327 DOI: 10.1093/nar/gkae541] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/04/2024] [Revised: 05/23/2024] [Accepted: 06/16/2024] [Indexed: 06/27/2024] Open
Abstract
Accurate RNA structure models are crucial for designing small molecule ligands that modulate their functions. This study assesses six standalone RNA 3D structure prediction methods-DeepFoldRNA, RhoFold, BRiQ, FARFAR2, SimRNA and Vfold2, excluding web-based tools due to intellectual property concerns. We focus on reproducing the RNA structure existing in RNA-small molecule complexes, particularly on the ability to model ligand binding sites. Using a comprehensive set of RNA structures from the PDB, which includes diverse structural elements, we found that machine learning (ML)-based methods effectively predict global RNA folds but are less accurate with local interactions. Conversely, non-ML-based methods demonstrate higher precision in modeling intramolecular interactions, particularly with secondary structure restraints. Importantly, ligand-binding site accuracy can remain sufficiently high for practical use, even if the overall model quality is not optimal. With the recent release of AlphaFold 3, we included this advanced method in our tests. Benchmark subsets containing new structures, not used in the training of the tested ML methods, show that AlphaFold 3's performance was comparable to other ML-based methods, albeit with some challenges in accurately modeling ligand binding sites. This study underscores the importance of enhancing binding site prediction accuracy and the challenges in modeling RNA-ligand interactions accurately.
Collapse
Affiliation(s)
- Chandran Nithin
- Molecure SA, 02-089 Warsaw, Poland
- Laboratory of Computational Biology, Biological and Chemical Research Center, Faculty of Chemistry, University of Warsaw, 02-089 Warsaw, Poland
| | - Sebastian Kmiecik
- Laboratory of Computational Biology, Biological and Chemical Research Center, Faculty of Chemistry, University of Warsaw, 02-089 Warsaw, Poland
| | | | | | | |
Collapse
|
3
|
Bernard C, Postic G, Ghannay S, Tahi F. State-of-the-RNArt: benchmarking current methods for RNA 3D structure prediction. NAR Genom Bioinform 2024; 6:lqae048. [PMID: 38745991 PMCID: PMC11091930 DOI: 10.1093/nargab/lqae048] [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: 02/08/2024] [Revised: 04/05/2024] [Accepted: 05/08/2024] [Indexed: 05/16/2024] Open
Abstract
RNAs are essential molecules involved in numerous biological functions. Understanding RNA functions requires the knowledge of their 3D structures. Computational methods have been developed for over two decades to predict the 3D conformations from RNA sequences. These computational methods have been widely used and are usually categorised as either ab initio or template-based. The performances remain to be improved. Recently, the rise of deep learning has changed the sight of novel approaches. Deep learning methods are promising, but their adaptation to RNA 3D structure prediction remains difficult. In this paper, we give a brief review of the ab initio, template-based and novel deep learning approaches. We highlight the different available tools and provide a benchmark on nine methods using the RNA-Puzzles dataset. We provide an online dashboard that shows the predictions made by benchmarked methods, freely available on the EvryRNA platform: https://evryrna.ibisc.univ-evry.fr/evryrna/state_of_the_rnart/.
Collapse
Affiliation(s)
- Clément Bernard
- Université Paris-Saclay, Univ. Evry, IBISC, 91020 Evry-Courcouronnes, France
- LISN - CNRS/Université Paris-Saclay, 91400 Orsay, France
| | - Guillaume Postic
- Université Paris-Saclay, Univ. Evry, IBISC, 91020 Evry-Courcouronnes, France
| | - Sahar Ghannay
- LISN - CNRS/Université Paris-Saclay, 91400 Orsay, France
| | - Fariza Tahi
- Université Paris-Saclay, Univ. Evry, IBISC, 91020 Evry-Courcouronnes, France
| |
Collapse
|
4
|
Riveros II, Yildirim I. Prediction of 3D RNA Structures from Sequence Using Energy Landscapes of RNA Dimers: Application to RNA Tetraloops. J Chem Theory Comput 2024; 20:4363-4376. [PMID: 38728627 DOI: 10.1021/acs.jctc.4c00189] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/12/2024]
Abstract
Access to the three-dimensional structure of RNA enables an ability to gain a more profound understanding of its biological mechanisms, as well as the ability to design RNA-targeting drugs, which can take advantage of the unique chemical environment imposed by a folded RNA structure. Due to the dynamic and structurally complex properties of RNA, both experimental and traditional computational methods have difficulty in determining RNA's 3D structure. Herein, we introduce TAPERSS (Theoretical Analyses, Prediction, and Evaluation of RNA Structures from Sequence), a physics-based fragment assembly method for predicting 3D RNA structures from sequence. Using a fragment library created using discrete path sampling calculations of RNA dinucleoside monophosphates, TAPERSS can sample the physics-based energy landscapes of any RNA sequence with relatively low computational complexity. We have benchmarked TAPERSS on 21 RNA tetraloops, using a combinatorial algorithm as a proof-of-concept. We show that TAPERSS was successfully able to predict the apo-state structures of all 21 RNA hairpins, with 16 of those structures also having low predicted energies as well. We demonstrate that TAPERSS performs most accurately on GNRA-like tetraloops with mostly stacked loop-nucleotides, while having limited success with more dynamic UNCG and CUYG tetraloops, most likely due to the influence of the RNA force field used to create the fragment library. Moreover, we show that TAPERSS can successfully predict the majority of the experimental non-apo states, highlighting its potential in anticipating biologically significant yet unobserved states. This holds great promise for future applications in drug design and related studies. With discussed improvements and implementation of more efficient sampling algorithms, we believe TAPERSS may serve as a useful tool for a physics-based conformational sampling of large RNA structures.
Collapse
Affiliation(s)
- Ivan Isaac Riveros
- Department of Chemistry and Biochemistry, Florida Atlantic University, Jupiter, Florida 33458, United States
| | - Ilyas Yildirim
- Department of Chemistry and Biochemistry, Florida Atlantic University, Jupiter, Florida 33458, United States
| |
Collapse
|
5
|
Zhang TT, Yan CL, Qiao JX, Yang AS, Liu ML, Kou YX, Li ZH. Demographic dynamics and molecular evolution of the rare and endangered subsect. Gerardianae of Pinus: insights from chloroplast genomes and mitochondrial DNA markers. PLANTA 2024; 259:45. [PMID: 38281265 DOI: 10.1007/s00425-023-04316-8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/11/2023] [Accepted: 12/21/2023] [Indexed: 01/30/2024]
Abstract
MAIN CONCLUSION The divergence of subsect. Gerardianae was likely triggered by the uplift of the Qinghai-Tibetan Plateau and adjacent mountains. Pinus bungeana might have probably experienced expansion since Last Interglacial period. Historical geological and climatic oscillations have profoundly affected patterns of nucleotide variability, evolutionary history, and species divergence in numerous plants of the Northern Hemisphere. However, how long-lived conifers responded to geological and climatic fluctuations in East Asia remain poorly understood. Here, based on paternally inherited chloroplast genomes and maternally inherited mitochondrial DNA markers, we investigated the population demographic history and molecular evolution of subsect. Gerardianae (only including three species, Pinus bungeana, P. gerardiana, and P. squamata) of Pinus. A low level of nucleotide diversity was found in P. bungeana (π was 0.00016 in chloroplast DNA sequences, and 0.00304 in mitochondrial DNAs). The haplotype-based phylogenetic topology and unimodal distributions of demographic analysis suggested that P. bungeana probably originated in the southern Qinling Mountains and experienced rapid population expansion since Last Interglacial period. Phylogenetic analysis revealed that P. gerardiana and P. squamata had closer genetic relationship. The species divergence of subsect. Gerardianae occurred about 27.18 million years ago (Mya) during the middle to late Oligocene, which was significantly associated with the uplift of the Qinghai-Tibetan Plateau and adjacent mountains from the Eocene to the mid-Pliocene. The molecular evolutionary analysis showed that two chloroplast genes (psaI and ycf1) were under positive selection, the genetic lineages of P. bungeana exhibited higher transition and nonsynonymous mutations, which were involved with the strongly environmental adaptation. These findings shed light on the population evolutionary history of white pine species and provide striking insights for comprehension of their species divergence and molecular evolution.
Collapse
Affiliation(s)
- Ting-Ting Zhang
- Key Laboratory of Resource Biology and Biotechnology in Western China (Ministry of Education), College of Life Sciences, Northwest University, Xi'an, 710069, China
| | - Chun-Li Yan
- Key Laboratory of Resource Biology and Biotechnology in Western China (Ministry of Education), College of Life Sciences, Northwest University, Xi'an, 710069, China
| | - Jin-Xia Qiao
- Key Laboratory of Resource Biology and Biotechnology in Western China (Ministry of Education), College of Life Sciences, Northwest University, Xi'an, 710069, China
| | - Ao-Shuang Yang
- Key Laboratory of Resource Biology and Biotechnology in Western China (Ministry of Education), College of Life Sciences, Northwest University, Xi'an, 710069, China
| | - Mi-Li Liu
- Key Laboratory of Resource Biology and Biotechnology in Western China (Ministry of Education), College of Life Sciences, Northwest University, Xi'an, 710069, China
| | - Yi-Xuan Kou
- Laboratory of Subtropical Biodiversity, Jiangxi Agricultural University, Nanchang, 330045, China
| | - Zhong-Hu Li
- Key Laboratory of Resource Biology and Biotechnology in Western China (Ministry of Education), College of Life Sciences, Northwest University, Xi'an, 710069, China.
| |
Collapse
|
6
|
Chen K, Zhou Y, Wang S, Xiong P. RNA tertiary structure modeling with BRiQ potential in CASP15. Proteins 2023; 91:1771-1778. [PMID: 37638558 DOI: 10.1002/prot.26574] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/15/2023] [Revised: 06/22/2023] [Accepted: 08/08/2023] [Indexed: 08/29/2023]
Abstract
We describe the modeling method for RNA tertiary structures employed by team AIchemy_RNA2 in the 15th Critical Assessment of Structure Prediction (CASP15). The method consists of the following steps. Firstly, secondary structure information was derived from various manually-verified sources. With this information, the full length RNA was fragmented into structural modules. The structures of each module were predicted and then assembled into the full structure. To reduce the searching conformational space, an RNA structure was organized into an optimal base folding tree. And to further improve the sampling efficiency, the energy surface was smoothed at high temperatures during the Monte Carlo sampling to make it easier to move across the energy barrier. The statistical potential energy function BRiQ was employed during Monte Carlo energy optimization.
Collapse
Affiliation(s)
- Ke Chen
- University of Science and Technology of China, Hefei, China
- Suzhou Institute for Advanced Research, University of Science and Technology of China, Suzhou, China
| | - Yaoqi Zhou
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen, China
| | - Sheng Wang
- Shanghai Zelixir Biotech Co. Ltd, Shanghai, China
| | - Peng Xiong
- University of Science and Technology of China, Hefei, China
- Suzhou Institute for Advanced Research, University of Science and Technology of China, Suzhou, China
| |
Collapse
|
7
|
Hörberg J, Reymer A. Decoding the dual recognition mechanism of the glucocorticoid receptor for DNA and RNA: sequence versus shape. Sci Rep 2023; 13:16125. [PMID: 37752333 PMCID: PMC10522765 DOI: 10.1038/s41598-023-43244-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/03/2023] [Accepted: 09/21/2023] [Indexed: 09/28/2023] Open
Abstract
Transcription factors (TFs) regulate eukaryotic transcription through selective DNA-binding, can also specifically interact with RNA, which may present another layer of transcriptional control. The mechanisms of the TFs-DNA recognition are often well-characterised, while the details of TFs-RNA complexation are less understood. Here we investigate the dual recognition mechanism of the glucocorticoid receptor (GR), which interacts with similar affinities with consensus DNA and diverse RNA hairpin motifs but discriminates against uniform dsRNA. Using atomic molecular dynamics simulations, we demonstrate that the GR binding to nucleic acids requires a wide and shallow groove pocket. The protein effectively moulds its binding site within DNA major groove, which enables base-specific interactions. Contrary, the GR binding has little effect on the grooves geometry of RNA systems, most notably in uniform dsRNA. Instead, a hairpin motif in RNA yields a wide and shallow major groove pocket, allowing the protein to anchor itself through nonspecific electrostatic contacts with RNA backbone. Addition of a bulge increases RNA hairpin flexibility, which leads to a greater number of GR-RNA contacts and, thus, higher affinity. Thus, the combination of structural motifs defines the GR-RNA selective binding: a recognition mechanism, which may be shared by other zinc finger TFs.
Collapse
Affiliation(s)
- Johanna Hörberg
- Department of Chemistry and Molecular Biology, University of Gothenburg, 405 30, Göteborg, Sweden
| | - Anna Reymer
- Department of Chemistry and Molecular Biology, University of Gothenburg, 405 30, Göteborg, Sweden.
| |
Collapse
|
8
|
Kulkarni M, Thangappan J, Deb I, Wu S. Comparative analysis of RNA secondary structure accuracy on predicted RNA 3D models. PLoS One 2023; 18:e0290907. [PMID: 37656749 PMCID: PMC10473517 DOI: 10.1371/journal.pone.0290907] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/10/2023] [Accepted: 08/18/2023] [Indexed: 09/03/2023] Open
Abstract
RNA structure is conformationally dynamic, and accurate all-atom tertiary (3D) structure modeling of RNA remains challenging with the prevailing tools. Secondary structure (2D) information is the standard prerequisite for most RNA 3D modeling. Despite several 2D and 3D structure prediction tools proposed in recent years, one of the challenges is to choose the best combination for accurate RNA 3D structure prediction. Here, we benchmarked seven small RNA PDB structures (40 to 90 nucleotides) with different topologies to understand the effects of different 2D structure predictions on the accuracy of 3D modeling. The current study explores the blind challenge of 2D to 3D conversions and highlights the performances of de novo RNA 3D modeling from their predicted 2D structure constraints. Our results show that conformational sampling-based methods such as SimRNA and IsRNA1 depend less on 2D accuracy, whereas motif-based methods account for 2D evidence. Our observations illustrate the disparities in available 3D and 2D prediction methods and may further offer insights into developing topology-specific or family-specific RNA structure prediction pipelines.
Collapse
Affiliation(s)
- Mandar Kulkarni
- R&D Center, PharmCADD Co. Ltd., Dong-gu, Busan, Republic of Korea
| | | | - Indrajit Deb
- R&D Center, PharmCADD Co. Ltd., Dong-gu, Busan, Republic of Korea
| | - Sangwook Wu
- R&D Center, PharmCADD Co. Ltd., Dong-gu, Busan, Republic of Korea
- Department of Physics, Pukyong National University, Busan, Republic of Korea
| |
Collapse
|
9
|
Wang X, Yu S, Lou E, Tan YL, Tan ZJ. RNA 3D Structure Prediction: Progress and Perspective. Molecules 2023; 28:5532. [PMID: 37513407 PMCID: PMC10386116 DOI: 10.3390/molecules28145532] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/14/2023] [Revised: 07/05/2023] [Accepted: 07/13/2023] [Indexed: 07/30/2023] Open
Abstract
Ribonucleic acid (RNA) molecules play vital roles in numerous important biological functions such as catalysis and gene regulation. The functions of RNAs are strongly coupled to their structures or proper structure changes, and RNA structure prediction has been paid much attention in the last two decades. Some computational models have been developed to predict RNA three-dimensional (3D) structures in silico, and these models are generally composed of predicting RNA 3D structure ensemble, evaluating near-native RNAs from the structure ensemble, and refining the identified RNAs. In this review, we will make a comprehensive overview of the recent advances in RNA 3D structure modeling, including structure ensemble prediction, evaluation, and refinement. Finally, we will emphasize some insights and perspectives in modeling RNA 3D structures.
Collapse
Affiliation(s)
- Xunxun Wang
- Department of Physics, Key Laboratory of Artificial Micro & Nano-Structures of Ministry of Education, School of Physics and Technology, Wuhan University, Wuhan 430072, China
| | - Shixiong Yu
- Department of Physics, Key Laboratory of Artificial Micro & Nano-Structures of Ministry of Education, School of Physics and Technology, Wuhan University, Wuhan 430072, China
| | - En Lou
- Department of Physics, Key Laboratory of Artificial Micro & Nano-Structures of Ministry of Education, School of Physics and Technology, Wuhan University, Wuhan 430072, China
| | - Ya-Lan Tan
- School of Bioengineering and Health, Wuhan Textile University, Wuhan 430200, China
- Research Center of Nonlinear Science, School of Mathematical and Physical Sciences, Wuhan Textile University, Wuhan 430200, China
| | - Zhi-Jie Tan
- Department of Physics, Key Laboratory of Artificial Micro & Nano-Structures of Ministry of Education, School of Physics and Technology, Wuhan University, Wuhan 430072, China
| |
Collapse
|
10
|
Waqqar S, Lee K, Lawley B, Bilton T, Quiñones-Mateu ME, Bostina M, Burga LN. Directed Evolution of Seneca Valley Virus in Tumorsphere and Monolayer Cell Cultures of a Small-Cell Lung Cancer Model. Cancers (Basel) 2023; 15:cancers15092541. [PMID: 37174006 PMCID: PMC10177334 DOI: 10.3390/cancers15092541] [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: 03/16/2023] [Revised: 04/18/2023] [Accepted: 04/26/2023] [Indexed: 05/15/2023] Open
Abstract
The Seneca Valley virus (SVV) is an oncolytic virus from the picornavirus family, characterized by a 7.3-kilobase RNA genome encoding for all the structural and functional viral proteins. Directed evolution by serial passaging has been employed for oncolytic virus adaptation to increase the killing efficacy towards certain types of tumors. We propagated the SVV in a small-cell lung cancer model under two culture conditions: conventional cell monolayer and tumorspheres, with the latter resembling more closely the cellular structure of the tumor of origin. We observed an increase of the virus-killing efficacy after ten passages in the tumorspheres. Deep sequencing analyses showed genomic changes in two SVV populations comprising 150 single nucleotides variants and 72 amino acid substitutions. Major differences observed in the tumorsphere-passaged virus population, compared to the cell monolayer, were identified in the conserved structural protein VP2 and in the highly variable P2 region, suggesting that the increase in the ability of the SVV to kill cells over time in the tumorspheres is acquired by capsid conservation and positively selecting mutations to counter the host innate immune responses.
Collapse
Affiliation(s)
- Shakeel Waqqar
- Department of Microbiology and Immunology, University of Otago, Dunedin 9016, New Zealand
| | - Kai Lee
- Department of Microbiology and Immunology, University of Otago, Dunedin 9016, New Zealand
| | - Blair Lawley
- Department of Microbiology and Immunology, University of Otago, Dunedin 9016, New Zealand
| | - Timothy Bilton
- Invermay Agricultural Centre, AgResearch, Mosgiel 9092, New Zealand
| | | | - Mihnea Bostina
- Department of Microbiology and Immunology, University of Otago, Dunedin 9016, New Zealand
| | - Laura N Burga
- Department of Microbiology and Immunology, University of Otago, Dunedin 9016, New Zealand
| |
Collapse
|
11
|
Wang X, Tan YL, Yu S, Shi YZ, Tan ZJ. Predicting 3D structures and stabilities for complex RNA pseudoknots in ion solutions. Biophys J 2023; 122:1503-1516. [PMID: 36924021 PMCID: PMC10147842 DOI: 10.1016/j.bpj.2023.03.017] [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: 10/24/2022] [Revised: 03/05/2023] [Accepted: 03/10/2023] [Indexed: 03/17/2023] Open
Abstract
RNA pseudoknots are a kind of important tertiary motif, and the structures and stabilities of pseudoknots are generally critical to the biological functions of RNAs with the motifs. In this work, we have carefully refined our previously developed coarse-grained model with salt effect through involving a new coarse-grained force field and a replica-exchange Monte Carlo algorithm, and employed the model to predict structures and stabilities of complex RNA pseudoknots in ion solutions beyond minimal H-type pseudoknots. Compared with available experimental data, the newly refined model can successfully predict 3D structures from sequences for the complex RNA pseudoknots including SARS-CoV-2 programming-1 ribosomal frameshifting element and Zika virus xrRNA, and can reliably predict the thermal stabilities of RNA pseudoknots with various sequences and lengths over broad ranges of monovalent/divalent salts. In addition, for complex pseudoknots including SARS-CoV-2 frameshifting element, our analyses show that their thermally unfolding pathways are mainly dependent on the relative stabilities of unfolded intermediate states, in analogy to those of minimal H-type pseudoknots.
Collapse
Affiliation(s)
- Xunxun Wang
- Department of Physics and Key Laboratory of Artificial Micro & Nano-structures of Ministry of Education, School of Physics and Technology, Wuhan University, Wuhan, China
| | - Ya-Lan Tan
- Research Center of Nonlinear Science and School of Mathematical and Physical Sciences, Wuhan Textile University, Wuhan, China
| | - Shixiong Yu
- Department of Physics and Key Laboratory of Artificial Micro & Nano-structures of Ministry of Education, School of Physics and Technology, Wuhan University, Wuhan, China
| | - Ya-Zhou Shi
- Research Center of Nonlinear Science and School of Mathematical and Physical Sciences, Wuhan Textile University, Wuhan, China
| | - Zhi-Jie Tan
- Department of Physics and Key Laboratory of Artificial Micro & Nano-structures of Ministry of Education, School of Physics and Technology, Wuhan University, Wuhan, China.
| |
Collapse
|
12
|
Arteaga SJ, Adams MS, Meyer NL, Richardson KE, Hoener S, Znosko BM. Thermodynamic determination of RNA duplex stability in magnesium solutions. Biophys J 2023; 122:565-576. [PMID: 36540026 PMCID: PMC9941723 DOI: 10.1016/j.bpj.2022.12.025] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/08/2022] [Revised: 12/05/2022] [Accepted: 12/16/2022] [Indexed: 12/24/2022] Open
Abstract
The prediction of RNA secondary structure and thermodynamics from sequence relies on free energy minimization and nearest neighbor parameters. Currently, algorithms used to make these predictions are based on parameters from optical melting studies performed in 1 M NaCl. However, many physiological and biochemical buffers containing RNA include much lower concentrations of monovalent cations and the presence of divalent cations. In order to improve these algorithms, thermodynamic data was previously collected for RNA duplexes in solutions containing 71, 121, 221, and 621 mM Na+. From this data, correction factors for free energy (ΔG°37) and melting temperature (Tm) were derived. Despite these newly derived correction factors for sodium, the stabilizing effects of magnesium have been ignored. Here, the same RNA duplexes were melted in solutions containing 0.5, 1.5, 3.0, and 10.0 mM Mg2+ in the absence of monovalent cations. Correction factors for Tm and ΔG°37 were derived to scale the current parameters to a range of magnesium concentrations. The Tm correction factor predicts the melting temperature within 1.2°C, and the ΔG°37 correction factor predicts the free energy within 0.30 kcalmol. These newly derived magnesium correction factors can be incorporated into algorithms that predict RNA secondary structure and stability from sequence.
Collapse
Affiliation(s)
| | - Miranda S Adams
- Department of Chemistry, Saint Louis University, Saint Louis, Missouri
| | - Nicole L Meyer
- Department of Chemistry, Saint Louis University, Saint Louis, Missouri
| | | | - Scott Hoener
- Department of Chemistry, Saint Louis University, Saint Louis, Missouri
| | - Brent M Znosko
- Department of Chemistry, Saint Louis University, Saint Louis, Missouri.
| |
Collapse
|
13
|
Abstract
RNA molecules carry out various cellular functions, and understanding the mechanisms behind their functions requires the knowledge of their 3D structures. Different types of computational methods have been developed to model RNA 3D structures over the past decade. These methods were widely used by researchers although their performance needs to be further improved. Recently, along with these traditional methods, machine-learning techniques have been increasingly applied to RNA 3D structure prediction and show significant improvement in performance. Here we shall give a brief review of the traditional methods and recent related advances in machine-learning approaches for RNA 3D structure prediction.
Collapse
Affiliation(s)
- Xiujuan Ou
- Institute of Biophysics, School of Physics, Huazhong University of Science and Technology, Wuhan 430074, Hubei, China
| | - Yi Zhang
- Institute of Biophysics, School of Physics, Huazhong University of Science and Technology, Wuhan 430074, Hubei, China
| | - Yiduo Xiong
- Institute of Biophysics, School of Physics, Huazhong University of Science and Technology, Wuhan 430074, Hubei, China
| | - Yi Xiao
- Institute of Biophysics, School of Physics, Huazhong University of Science and Technology, Wuhan 430074, Hubei, China
| |
Collapse
|
14
|
Sun D, Sun M, Zhang J, Lin X, Zhang Y, Lin F, Zhang P, Yang C, Song J. Computational tools for aptamer identification and optimization. Trends Analyt Chem 2022. [DOI: 10.1016/j.trac.2022.116767] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
|
15
|
Zhou L, Wang X, Yu S, Tan YL, Tan ZJ. FebRNA: An automated fragment-ensemble-based model for building RNA 3D structures. Biophys J 2022; 121:3381-3392. [PMID: 35978551 PMCID: PMC9515226 DOI: 10.1016/j.bpj.2022.08.017] [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: 04/30/2022] [Revised: 07/19/2022] [Accepted: 08/15/2022] [Indexed: 11/23/2022] Open
Abstract
Knowledge of RNA three-dimensional (3D) structures is critical to understanding the important biological functions of RNAs. Although various structure prediction models have been developed, the high-accuracy predictions of RNA 3D structures are still limited to the RNAs with short lengths or with simple topology. In this work, we proposed a new model, namely FebRNA, for building RNA 3D structures through fragment assembly based on coarse-grained (CG) fragment ensembles. Specifically, FebRNA is composed of four processes: establishing the library of different types of non-redundant CG fragment ensembles regardless of the sequences, building CG 3D structure ensemble through fragment assembly, identifying top-scored CG structures through a specific CG scoring function, and rebuilding the all-atom structures from the top-scored CG ones. Extensive examination against different types of RNA structures indicates that FebRNA consistently gives the reliable predictions on RNA 3D structures, including pseudoknots, three-way junctions, four-way and five-way junctions, and RNAs in the RNA-Puzzles. FebRNA is available on the Web site: https://github.com/Tan-group/FebRNA.
Collapse
Affiliation(s)
- Li Zhou
- Department of Physics and Key Laboratory of Artificial Micro & Nano-structures of Education, School of Physics and Technology, Wuhan University, Wuhan 430072, China
| | - Xunxun Wang
- Department of Physics and Key Laboratory of Artificial Micro & Nano-structures of Education, School of Physics and Technology, Wuhan University, Wuhan 430072, China
| | - Shixiong Yu
- Department of Physics and Key Laboratory of Artificial Micro & Nano-structures of Education, School of Physics and Technology, Wuhan University, Wuhan 430072, China
| | - Ya-Lan Tan
- Research Center of Nonlinear Science, School of Mathematical and Physical Sciences, Wuhan Textile University, Wuhan 430073, China.
| | - Zhi-Jie Tan
- Department of Physics and Key Laboratory of Artificial Micro & Nano-structures of Education, School of Physics and Technology, Wuhan University, Wuhan 430072, China.
| |
Collapse
|
16
|
3dDNA: A Computational Method of Building DNA 3D Structures. Molecules 2022; 27:molecules27185936. [PMID: 36144680 PMCID: PMC9503956 DOI: 10.3390/molecules27185936] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/19/2022] [Revised: 09/04/2022] [Accepted: 09/08/2022] [Indexed: 02/07/2023] Open
Abstract
Considerable progress has been made in the prediction methods of 3D structures of RNAs. In contrast, no such methods are available for DNAs. The determination of 3D structures of the latter is also increasingly needed for understanding their functions and designing new DNA molecules. Since the number of experimental structures of DNA is limited at present, here, we propose a computational and template-based method, 3dDNA, which combines DNA and RNA template libraries to predict DNA 3D structures. It was benchmarked on three test sets with different numbers of chains, and the results show that 3dDNA can predict DNA 3D structures with a mean RMSD of about 2.36 Å for those with one or two chains and fewer than 4 Å with three or more chains.
Collapse
|
17
|
Xu B, Zhu Y, Cao C, Chen H, Jin Q, Li G, Ma J, Yang SL, Zhao J, Zhu J, Ding Y, Fang X, Jin Y, Kwok CK, Ren A, Wan Y, Wang Z, Xue Y, Zhang H, Zhang QC, Zhou Y. Recent advances in RNA structurome. SCIENCE CHINA. LIFE SCIENCES 2022; 65:1285-1324. [PMID: 35717434 PMCID: PMC9206424 DOI: 10.1007/s11427-021-2116-2] [Citation(s) in RCA: 16] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Subscribe] [Scholar Register] [Received: 03/10/2022] [Accepted: 04/01/2022] [Indexed: 12/27/2022]
Abstract
RNA structures are essential to support RNA functions and regulation in various biological processes. Recently, a range of novel technologies have been developed to decode genome-wide RNA structures and novel modes of functionality across a wide range of species. In this review, we summarize key strategies for probing the RNA structurome and discuss the pros and cons of representative technologies. In particular, these new technologies have been applied to dissect the structural landscape of the SARS-CoV-2 RNA genome. We also summarize the functionalities of RNA structures discovered in different regulatory layers-including RNA processing, transport, localization, and mRNA translation-across viruses, bacteria, animals, and plants. We review many versatile RNA structural elements in the context of different physiological and pathological processes (e.g., cell differentiation, stress response, and viral replication). Finally, we discuss future prospects for RNA structural studies to map the RNA structurome at higher resolution and at the single-molecule and single-cell level, and to decipher novel modes of RNA structures and functions for innovative applications.
Collapse
Affiliation(s)
- Bingbing Xu
- MOE Laboratory of Biosystems Homeostasis & Protection, Innovation Center for Cell Signaling Network, College of Life Sciences, Zhejiang University, Hangzhou, 310058, China
| | - Yanda Zhu
- MOE Laboratory of Biosystems Homeostasis & Protection, Innovation Center for Cell Signaling Network, College of Life Sciences, Zhejiang University, Hangzhou, 310058, China
| | - Changchang Cao
- Key Laboratory of RNA Biology, Institute of Biophysics, Chinese Academy of Sciences, Beijing, 100101, China
| | - Hao Chen
- Life Sciences Institute, Zhejiang University, Hangzhou, 310058, China
| | - Qiongli Jin
- State Key Laboratory of Plant Physiology and Biochemistry, College of Life Sciences, Zhejiang University, Hangzhou, 310058, China
| | - Guangnan Li
- State Key Laboratory of Virology, College of Life Sciences, Wuhan University, Wuhan, 430072, China
| | - Junfeng Ma
- Beijing Advanced Innovation Center for Structural Biology, School of Life Sciences, Tsinghua University, Beijing, 100084, China
| | - Siwy Ling Yang
- Stem Cell and Regenerative Biology, Genome Institute of Singapore, A*STAR, Singapore, Singapore
| | - Jieyu Zhao
- Department of Chemistry, and State Key Laboratory of Marine Pollution, City University of Hong Kong, Kowloon Tong, Hong Kong SAR, China
| | - Jianghui Zhu
- MOE Key Laboratory of Bioinformatics, Beijing Advanced Innovation Center for Structural Biology and Frontier Research Center for Biological Structure, Center for Synthetic and Systems Biology, School of Life Sciences, Tsinghua University, Beijing, 100084, China
- Tsinghua-Peking Center for Life Sciences, Beijing, 100084, China
| | - Yiliang Ding
- Department of Cell and Developmental Biology, John Innes Centre, Norwich Research Park, Norwich, NR4 7UH, United Kingdom.
| | - Xianyang Fang
- Beijing Advanced Innovation Center for Structural Biology, School of Life Sciences, Tsinghua University, Beijing, 100084, China.
| | - Yongfeng Jin
- MOE Laboratory of Biosystems Homeostasis & Protection, Innovation Center for Cell Signaling Network, College of Life Sciences, Zhejiang University, Hangzhou, 310058, China.
| | - Chun Kit Kwok
- Department of Chemistry, and State Key Laboratory of Marine Pollution, City University of Hong Kong, Kowloon Tong, Hong Kong SAR, China.
- Shenzhen Research Institute of City University of Hong Kong, Shenzhen, 518057, China.
| | - Aiming Ren
- Life Sciences Institute, Zhejiang University, Hangzhou, 310058, China.
| | - Yue Wan
- Stem Cell and Regenerative Biology, Genome Institute of Singapore, A*STAR, Singapore, Singapore.
| | - Zhiye Wang
- State Key Laboratory of Plant Physiology and Biochemistry, College of Life Sciences, Zhejiang University, Hangzhou, 310058, China.
| | - Yuanchao Xue
- Key Laboratory of RNA Biology, Institute of Biophysics, Chinese Academy of Sciences, Beijing, 100101, China.
- University of Chinese Academy of Sciences, Beijing, 100101, China.
| | - Huakun Zhang
- Key Laboratory of Molecular Epigenetics of the Ministry of Education, Northeast Normal University, Changchun, 130024, China.
| | - Qiangfeng Cliff Zhang
- MOE Key Laboratory of Bioinformatics, Beijing Advanced Innovation Center for Structural Biology and Frontier Research Center for Biological Structure, Center for Synthetic and Systems Biology, School of Life Sciences, Tsinghua University, Beijing, 100084, China.
- Tsinghua-Peking Center for Life Sciences, Beijing, 100084, China.
| | - Yu Zhou
- State Key Laboratory of Virology, College of Life Sciences, Wuhan University, Wuhan, 430072, China.
| |
Collapse
|
18
|
Singh J, Paliwal K, Litfin T, Singh J, Zhou Y. Predicting RNA distance-based contact maps by integrated deep learning on physics-inferred secondary structure and evolutionary-derived mutational coupling. Bioinformatics 2022; 38:3900-3910. [PMID: 35751593 PMCID: PMC9364379 DOI: 10.1093/bioinformatics/btac421] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/24/2021] [Revised: 04/30/2022] [Accepted: 06/28/2022] [Indexed: 12/24/2022] Open
Abstract
MOTIVATION Recently, AlphaFold2 achieved high experimental accuracy for the majority of proteins in Critical Assessment of Structure Prediction (CASP 14). This raises the hope that one day, we may achieve the same feat for RNA structure prediction for those structured RNAs, which is as fundamentally and practically important similar to protein structure prediction. One major factor in the recent advancement of protein structure prediction is the highly accurate prediction of distance-based contact maps of proteins. RESULTS Here, we showed that by integrated deep learning with physics-inferred secondary structures, co-evolutionary information and multiple sequence-alignment sampling, we can achieve RNA contact-map prediction at a level of accuracy similar to that in protein contact-map prediction. More importantly, highly accurate prediction for top L long-range contacts can be assured for those RNAs with a high effective number of homologous sequences (Neff > 50). The initial use of the predicted contact map as distance-based restraints confirmed its usefulness in 3D structure prediction. AVAILABILITY AND IMPLEMENTATION SPOT-RNA-2D is available as a web server at https://sparks-lab.org/server/spot-rna-2d/ and as a standalone program at https://github.com/jaswindersingh2/SPOT-RNA-2D. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
| | | | - Thomas Litfin
- Institute for Glycomics, Griffith University, Parklands Dr. Southport, QLD 4222, Australia
| | - Jaspreet Singh
- Signal Processing Laboratory, School of Engineering and Built Environment, Griffith University, Brisbane, QLD 4111, Australia
| | - Yaoqi Zhou
- To whom correspondence should be addressed. or or
| |
Collapse
|
19
|
Modelling aptamers with nucleic acid mimics (NAM): From sequence to three-dimensional docking. PLoS One 2022; 17:e0264701. [PMID: 35320268 PMCID: PMC8942228 DOI: 10.1371/journal.pone.0264701] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/22/2021] [Accepted: 02/15/2022] [Indexed: 02/07/2023] Open
Abstract
Aptamers are single-stranded oligonucleotides, formerly evolved by Systematic Evolution of Ligands by EXponential enrichment (SELEX), that fold into functional three-dimensional structures. Such conformation is crucial for aptamers' ability to bind to a target with high affinity and specificity. Unnatural nucleotides have been used to develop nucleic acid mimic (NAM) aptamers with increased performance, such as biological stability. Prior knowledge of aptamer-target interactions is critical for applying post-SELEX modifications with unnatural nucleotides since it can affect aptamers' structure and performance. Here, we describe an easy-to-apply in silico workflow using free available software / web servers to predict the tertiary conformation of NAM, DNA and RNA aptamers, as well as the docking with the target molecule. Representative 2'-O-methyl (2'OMe), locked nucleic acid (LNA), DNA and RNA aptamers, with experimental data deposited in Protein Data Bank, were selected to validate the workflow. All aptamers' tertiary structure and docking models were successfully predicted with good structural similarity to the experimental data. Thus, this workflow will boost the development of aptamers, particularly NAM aptamers, by assisting in the rational modification of specific nucleotides and avoiding trial-and-error approaches.
Collapse
|
20
|
Mao K, Wang J, Xiao Y. Length-Dependent Deep Learning Model for RNA Secondary Structure Prediction. MOLECULES (BASEL, SWITZERLAND) 2022; 27:molecules27031030. [PMID: 35164295 PMCID: PMC8838716 DOI: 10.3390/molecules27031030] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 12/09/2021] [Revised: 01/17/2022] [Accepted: 01/29/2022] [Indexed: 11/16/2022]
Abstract
Deep learning methods for RNA secondary structure prediction have shown higher performance than traditional methods, but there is still much room to improve. It is known that the lengths of RNAs are very different, as are their secondary structures. However, the current deep learning methods all use length-independent models, so it is difficult for these models to learn very different secondary structures. Here, we propose a length-dependent model that is obtained by further training the length-independent model for different length ranges of RNAs through transfer learning. 2dRNA, a coupled deep learning neural network for RNA secondary structure prediction, is used to do this. Benchmarking shows that the length-dependent model performs better than the usual length-independent model.
Collapse
|
21
|
3dRNA: 3D structure prediction from linear to circular RNAs. J Mol Biol 2022; 434:167452. [DOI: 10.1016/j.jmb.2022.167452] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/11/2021] [Revised: 01/07/2022] [Accepted: 01/07/2022] [Indexed: 12/30/2022]
|
22
|
3D Modeling of Non-coding RNA Interactions. ADVANCES IN EXPERIMENTAL MEDICINE AND BIOLOGY 2022; 1385:281-317. [DOI: 10.1007/978-3-031-08356-3_11] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
|
23
|
Zhou T, Wang H, Zeng C, Zhao Y. RPocket: an intuitive database of RNA pocket topology information with RNA-ligand data resources. BMC Bioinformatics 2021; 22:428. [PMID: 34496744 PMCID: PMC8424408 DOI: 10.1186/s12859-021-04349-4] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/10/2021] [Accepted: 08/27/2021] [Indexed: 11/13/2022] Open
Abstract
Background RNA regulates a variety of biological functions by interacting with other molecules. The ligand often binds in the RNA pocket to trigger structural changes or functions. Thus, it is essential to explore and visualize the RNA pocket to elucidate the structural and recognition mechanism for the RNA-ligand complex formation. Results In this work, we developed one user-friendly bioinformatics tool, RPocket. This database provides geometrical size, centroid, shape, secondary structure element for RNA pocket, RNA-ligand interaction information, and functional sites. We extracted 240 RNA pockets from 94 non-redundant RNA-ligand complex structures. We developed RPDescriptor to calculate the pocket geometrical property quantitatively. The geometrical information was then subjected to RNA-ligand binding analysis by incorporating the sequence, secondary structure, and geometrical combinations. This new approach takes advantage of both the atom-level precision of the structure and the nucleotide-level tertiary interactions. The results show that the higher-level topological pattern indeed improves the tertiary structure prediction. We also proposed a potential mechanism for RNA-ligand complex formation. The electrostatic interactions are responsible for long-range recognition, while the Van der Waals and hydrophobic contacts for short-range binding and optimization. These interaction pairs can be considered as distance constraints to guide complex structural modeling and drug design. Conclusion RPocket database would facilitate RNA-ligand engineering to regulate the complex formation for biological or medical applications. RPocket is available at http://zhaoserver.com.cn/RPocket/RPocket.html. Supplementary Information The online version contains supplementary material available at 10.1186/s12859-021-04349-4.
Collapse
Affiliation(s)
- Ting Zhou
- Department of Physics, Institute of Biophysics, Central China Normal University, Wuhan, 430079, China
| | - Huiwen Wang
- Department of Physics, Institute of Biophysics, Central China Normal University, Wuhan, 430079, China
| | - Chen Zeng
- Department of Physics, George Washington University, Washington, DC, 20052, USA
| | - Yunjie Zhao
- Department of Physics, Institute of Biophysics, Central China Normal University, Wuhan, 430079, China.
| |
Collapse
|
24
|
Learning the Fastest RNA Folding Path Based on Reinforcement Learning and Monte Carlo Tree Search. Molecules 2021; 26:molecules26154420. [PMID: 34361572 PMCID: PMC8347524 DOI: 10.3390/molecules26154420] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/07/2021] [Revised: 07/17/2021] [Accepted: 07/20/2021] [Indexed: 11/17/2022] Open
Abstract
RNA molecules participate in many important biological processes, and they need to fold into well-defined secondary and tertiary structures to realize their functions. Like the well-known protein folding problem, there is also an RNA folding problem. The folding problem includes two aspects: structure prediction and folding mechanism. Although the former has been widely studied, the latter is still not well understood. Here we present a deep reinforcement learning algorithms 2dRNA-Fold to study the fastest folding paths of RNA secondary structure. 2dRNA-Fold uses a neural network combined with Monte Carlo tree search to select residue pairing step by step according to a given RNA sequence until the final secondary structure is formed. We apply 2dRNA-Fold to several short RNA molecules and one longer RNA 1Y26 and find that their fastest folding paths show some interesting features. 2dRNA-Fold is further trained using a set of RNA molecules from the dataset bpRNA and is used to predict RNA secondary structure. Since in 2dRNA-Fold the scoring to determine next step is based on possible base pairings, the learned or predicted fastest folding path may not agree with the actual folding paths determined by free energy according to physical laws.
Collapse
|