How RNA sequences fold to specific tertiary structures is one of the key problems for understanding their dynamics and functions. Here, we study the folding process of an H-type RNA pseudoknot by performing a large-scale all-atom MD simulation and bias-exchange metadynamics. The folding free energy landscapes are obtained and several folding intermediates are identified. It is suggested that the folding occurs via multiple mechanisms, including a step-wise mechanism starting either from the first helix or the second, and a cooperative mechanism with both helices forming simultaneously. Despite of the multiple mechanism nature, the ensemble folding kinetics estimated from a Markov state model is single-exponential. It is also found that the correlation between folding and binding of metal ions is significant, and the bound ions mediate long-range interactions in the intermediate structures. Non-native interactions are found to be dominant in the unfolded state and also present in some intermediates, possibly hinder the folding process of the RNA.

Glioma is a lethal malignant brain cancer, and many reports have shown that abnormalities in the behavior of water and ion channels play an important role in regulating tumor proliferation, migration, apoptosis, and differentiation. Recently, new studies have suggested that some long noncoding RNAs containing small open reading frames can encode small peptides and form oligomers for water or ion regulation. However, because the peptides are difficult to identify, their functional mechanisms are far from being clearly understood. In this study, we used bioinformatics methods to identify and evaluate lncRNAs, which may encode small transmembrane peptides in gliomas. Combining ab initio homology modeling, molecular dynamics simulations, and free energy calculations, we constructed a predictive model and predicted the oligomer channel activity of peptides by identifying the lncRNA ORFs. We found that one key hub lncRNA, namely, DLEU1, which contains two smORFs (ORF1 and ORF8), encodes small peptides that form pentameric channels. The mechanics of water and ion (Na+ and Cl-) transport through this pentameric channel were simulated. The potential mean force of the H2O molecules along the two ORF-encoded peptide channels indicated that the energy barrier was different between ORF1 and ORF8. The ORF1-encoded peptide pentamer acted as a self-assembled water channel but not as an ion channel, and the ORF8 permeated neither ions nor water. This work provides new methods and theoretical support for further elucidation of the function of lncRNA-encoded small peptides and their role in cancer. Additionally, this study provides a theoretical basis for drug development.