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.

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

Biomolecular Topology: Modelling and Analysis

Contact-Map-Driven Exploration of Heterogeneous Protein-Folding Paths

Exploration of Multi-State Conformational Dynamics and Underlying Global Functional Landscape of Maltose Binding Protein

The Folding Transition State of Protein L is Extensive with Nonnative Interactions (and Not Small and Polarized)

A topological selection of folding pathways from native states of knotted proteins

A WAVELET APPROACH FOR THE ANALYSIS OF FOLDING TRAJECTORY OF PROTEIN TRP-CAGE

Systematically Constructing Kinetic Transition Network in Polypeptide from Top to Down: Trajectory Mapping

RNA Folding Pathways in Stop Motion

Complex Folding Pathways in a Simple Beta-Hairpin

Revealing Multiple Transition Paths During the Unfolding/Folding of Individual Peptides in a Confined Nanopore

Network Approach to Identify the Folding Transition States of Peptides and Proteins.

Minimal models for proteins and RNA: From folding to function

DiffTopo: Fold exploration using coarse grained protein topology representations

Persistent homology analysis of protein structure, flexibility and folding

Folding with a protein's native shortcut network

The Hidden Transition Paths During the Unfolding of Individual Peptides with a Confined Nanopore

Free Energy Landscape and Multiple Folding Pathways of an H-Type RNA Pseudoknot

Energy Landscape And Multiroute Folding Of Topologically Complex Proteins Adenylate Kinase And 2ouf-Knot

Discrete molecular dynamics studies of the folding of a protein-like model

Detecting intermediate protein conformations using algebraic topology