Topological Methods for Exploring Low-density States in Biomolecular Folding Pathways

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 uninteresting 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. It provides us a systematic tool to explore the low density intermediate states in complex biomolecular folding systems.

23 pages, 6 figures

Cited by in corpus (35)

• Representability of algebraic topology for biomolecules in machine learning based scoring and virtual screening
• Persistent homology analysis of protein structure, flexibility and folding
• A review of mathematical representations of biomolecules
• Hierarchical Nystrom Methods for Constructing Markov State Models for Conformational Dynamics
• Categorified Reeb Graphs
• Statistical Analysis and Parameter Selection for Mapper
• Topological portraits of multiscale coordination dynamics
• Computational topology for configuration spaces of hard disks
• Topological data analysis and machine learning
• Protein folding analysis using features obtained by persistent homology
• Persistent homology analysis of ion aggregation and hydrogen-bonding network
• Strong Equivalence of the Interleaving and Functional Distortion Metrics for Reeb Graphs
• On the choice of weight functions for linear representations of persistence diagrams
• Persistent homology analysis of osmolyte molecular aggregation and their hydrogen-bonding networks
• A Note on Learning Rare Events in Molecular Dynamics using LSTM and Transformer
• A review of geometric, topological and graph theory apparatuses for the modeling and analysis of biomolecular data
• Persistent Homology for The Quantitative Prediction of Fullerene Stability
• Generalized Persistence Algorithm for Decomposing Multi-parameter Persistence Modules
• Energy-Based Clustering: Fast and Robust Clustering of Data with Known Likelihood Functions
• Weighted persistent homology for biomolecular data analysis
• Multiscale persistent functions for biomolecular structure characterization
• Topological fingerprints reveal protein-ligand binding mechanism
• A Primer on Persistent Homology of Finite Metric Spaces
• A quantitative structure comparison with persistent similarity
• Classification of COVID-19 via Homology of CT-SCAN
• Weighted persistent homology for osmolyte molecular aggregation and hydrogen-bonding network analysis
• Persistent spectral graph
• Evolutionary de Rham-Hodge method
• Maximum Margin Clustering for State Decomposition of Metastable Systems
• The de Rham-Hodge analysis and modeling of biomolecules
• Multiresolution topological simplification
• Analysis and prediction of protein folding energy changes upon mutation by element specific persistent homology
• Atlasing of Assembly Landscapes using Distance Geometry and Graph Rigidity
• Efficient Atlasing and Search of Configuration Spaces of Point-Sets Constrained by Distance Intervals
• Mathematics at the eve of a historic transition in biology