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.

Exploring the Application of SiteMap and Site Finder for Focused Cryptic Pocket Identification

CrypToth: Cryptic pocket detection through mixed-solvent molecular dynamics simulations based topological data analysis

Exploration of cryptic pockets using enhanced sampling along normal modes: A case study of KRAS G12D

Divide and Conquer. Pocket-Opening Mixed-Solvent Simulations in the Perspective of Docking Virtual Screening Applications for Drug Discovery

Comprehensive detection and characterization of human druggable pockets through binding site descriptors

Comprehensive detection and characterization of human druggable pockets through novel binding site descriptors

Pocket2Mol: Efficient Molecular Sampling Based on 3D Protein Pockets.

Subpocket-Based Analysis Approach for the Protein Pocket Dynamics

Identification of Cavities on Protein Surface Using Multiple Computational Approaches for Drug Binding Site Prediction

Decipher Fundamental Atomic Interactions to Unify Generative Molecular Docking and Design

VAPPD: Visual Analysis of Protein Pocket Dynamics

Discovery of cryptic allosteric sites using reversed allosteric communication by a combined computational and experimental strategy

SWISH-X, an Expanded Approach to Detect Cryptic Pockets in Proteins and at Protein-Protein Interfaces

SiteFerret: Beyond Simple Pocket Identification in Proteins

Multiscaled Exploration of Coupled Folding and Binding of an Intrinsically Disordered Molecular Recognition Element in Measles Virus Nucleoprotein

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

Which cryptic sites are feasible drug targets?

Pocket Crafter: a 3D generative modeling based workflow for the rapid generation of hit molecules in drug discovery

Inhibition of striatal-enriched protein tyrosine phosphatase by targeting computationally revealed cryptic pockets

Quantum Chemistry in a Pocket: A Multifaceted Tool to Link Structure and Activity

Estimating the Similarity between Protein Pockets