Abstract:Characterization of transient intermediate or transition states is crucial for the description of biomolecular folding pathways, which is, however, difficult in both experiments and computer simulations. Such transient states are typically of low population in simulation samples. Even for simple systems such as RNA hairpins, recently there are mounting debates over the existence of multiple intermediate states. In this paper, we develop a computational approach to explore the relatively low populated transition or intermediate states in biomolecular folding pathways, based on a topological data analysis tool, MAPPER, with simulation data from large-scale distributed computing. The method is inspired by the classical Morse theory in mathematics which characterizes the topology of high-dimensional shapes via some functional level sets. In this paper we exploit a conditional density filter which enables us to focus on the structures on pathways, followed by clustering analysis on its level sets, which helps separate low populated intermediates from high populated folded/unfolded structures. A successful application of this method is given on a motivating example, a RNA hairpin with GCAA tetraloop, where we are able to provide structural evidence from computer simulations on the multiple intermediate states and exhibit different pictures about unfolding and refolding pathways. The method is effective in dealing with high degree of heterogeneity in distribution, capturing structural features in multiple pathways, and being less sensitive to the distance metric than nonlinear dimensionality reduction or geometric embedding methods. The methodology described in this paper admits various implementations or extensions to incorporate more information and adapt to different settings, which thus provides a systematic tool to explore the low-density intermediate states in complex biomolecular folding systems.

Genus Comparisons in the Topological Analysis of RNA Structures

Metrics for RNA Secondary Structure Comparison

Ab initio RNA folding

Influence of Contact Map Topology on RNA Structure Prediction

McGenus: A Monte Carlo algorithm to predict RNA secondary structures with pseudoknots

Comparison of Pseudoknotted RNA Secondary Structures by Topological Centroid Identification and Tree Edit Distance

RNA fold prediction by Monte Carlo in graph space and the statistical mechanics of tertiary interactions

Evaluation of RNA secondary structure prediction for both base-pairing and topology

A Comparative Review of Deep Learning Methods for RNA Tertiary Structure Prediction

Topological methods for exploring low-density states in biomolecular folding pathways

Review of RNA secondary structure prediction algorithms

Improved RNA pseudoknots prediction and classification using a new topological invariant

Statistical-physics approaches to RNA molecules, families and networks

RNA Secondary Structure Prediction Using a Genetic Algorithm with a Selection Method Based on Free Energy Value and Topological Index

ANNaMo: Coarse-grained modelling for folding and assembly of RNA and DNA systems

A Dynamic 3d Graphical Representation for Rna Structure Analysis and Its Application in Non-Coding Rna Classification

Parametric analysis of RNA folding

RNA Folding on the 3D Triangular Lattice

Helix-Based RNA Landscape Partition and Alternative Secondary Structure Determination

Shapes of interacting RNA complexes

A Novel Comparative Sequence Analysis Method for Ncrna Secondary Structure Prediction Without Multiple Sequence Alignment