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.

Fast topological analysis of 2D and 3D grids of data: Application to the atoms in molecule (AIM) and the electron localization function (ELF)

Molecular Structural Analysis Based on the Topological Extraction of 3D Scalar Fields

Fluid Flow Topology Extraction on 3D Vector Fields

Some Topological Indices of Polar Grid Graph

Numerical computation of critical properties and atomic basins from 3D grid electron densities

Finding Critical Points and Reconstruction of Electron Densities on Grids

Topological Analysis of Molecular Dynamics Simulations using the Euler Characteristic

A new DFT method for atoms and molecules in Cartesian grid

Topological Data Analysis: Concepts, Computation, and Applications in Chemical Engineering

Chemical structure and reactivity by means of quantum chemical topology analysis

A Topological Framework for Local Structure Analysis in Condensed Matter

A Fast and Scalable Computational Topology Framework for the Euler Characteristic

Analyzing the expressions for nanostructures via topological indices

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

Large-scale characterization of chemical bonding and topology in the Materials Project database

Density functional electric response properties of molecules in Cartesian grid

Grid-based state space exploration for molecular binding

A review of geometric, topological and graph theory apparatuses for the modeling and analysis of biomolecular data

Two-grid method on unstructured tetrahedra: Applying computational geometry to staggered solution of coupled flow and mechanics problems

How electrons still guard the space: Electron number distribution functions based on QTAIM∩ELF intersections

Topology-Aware Space-Shared Co-Analysis of Large-Scale Molecular Dynamics Simulations