“…Although experimental methods, such as X-ray crystallography, nuclear magnetic resonance spectroscopy, and cryo-electron microscopy, can be used to determine the structures of RNAs including pseudoknots, the structures in Protein Data Bank (PDB; https://www.rcsb.org) are still limited due to the high cost of the experimental measurements (Hajdin et al, 2010;Rose et al, 2011;Shi et al, 2014b;Schlick and Pyle, 2017). To complement the experiments, some computational models/methods (e.g., FARNA, MC-Fold/MC-Sym, Vfold, iFoldRNA, 3dRNA, RNAComposer, SimRNA, oxRNA, HiRE-RNA, and pk3D) have been developed for predicting RNA 3D structures (Cao and Chen, 2005;Ding et al, 2008;Parisien and Major, 2008;Zhang et al, 2009;Das et al, 2010;Popenda et al, 2012;Zhao et al, 2012;He et al, 2013He et al, , 2015He et al, , 2019Kim et al, 2014;Liwo et al, 2014Liwo et al, , 2020Sulc et al, 2014;Cragnolini et al, 2015;Wang et al, 2015a,b;Boniecki et al, 2016;Dawson et al, 2016;Li et al, 2016Li et al, , 2018Tan et al, 2019). Most of these models/methods are primarily designed to predict folded structures and cannot predict the stability of RNAs, especially in ion solutions (Shi et al, 2014b;Dawson et al, 2016;Schlick and Pyle, 2017), whereas the structural stability of RNAs can be very sensitive to ion conditions due to their polyanionic nature (Das et al, 2005;Draper et al, 2005;Chen, 2007, 2011;Qiu et al, 2010;Lipfert et al, 2014;Wang et al, 2018<...>…”