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.

Time-Resolved Graphs of Polymorphic Cycles for H-Bonded Network Identification in Flexible Biomolecules

Learning Effective Molecular Models from Experimental Observables.

A Molecular Dynamics Study of Bond Exchange Reactions in Covalent Adaptable Networks.

Graph theory approaches for molecular dynamics simulations

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

Network theory approach for data evaluation in the dynamic force spectroscopy of biomolecular interactions

Quantification of Cyclic Topology in Polymer Networks Using 3D Nets

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

Quantifying the Kinetic Paths of Flexible Biomolecular Recognition

Conformational Preference of Macrocycles Investigated by Ion-Mobility Mass Spectrometry and Distance Geometry Modeling.

Simulating Nanoscale Functional Motions of Biomolecules

A new graph modelisation for molecule similarity

Identification and Analysis of Transition and Metastable Markov States

Systematically Defining Coarse-Grained Representations of Large Biomolecular Complexes

Automated Event Detection and Activity Monitoring in Long Molecular Dynamics Simulations

Hydrogen bonding heterogeneity correlates with protein folding transition state passage time as revealed by data sonification

Training Neural Network Models Using Molecular Dynamics Simulation Results to Efficiently Predict Cyclic Hexapeptide Structural Ensembles

Graph‐based algorithms to dissect long‐distance water‐mediated H‐bond networks for conformational couplings in GPCRs

Collective Variables for Conformational Polymorphism in Molecular Crystals

Molecular dynamics-based approaches for enhanced sampling of long-time, large-scale conformational changes in biomolecules

Simulate Time-integrated Coarse-grained Molecular Dynamics with Multi-Scale Graph Networks