From Enhanced Sampling to Human-Readable Representations of Protein Dynamics

TL;DR

This paper introduces an automated method to convert complex enhanced sampling trajectories into simple, human-readable representations of protein dynamics, improving interpretability without sacrificing essential information.

Contribution

The authors develop a fully automated framework that transforms enhanced sampling data into interpretable geometric descriptors, facilitating understanding of protein motions without prior knowledge.

Findings

Successfully applied to five proteins, including KRAS and HIV-1 protease.

Consistently identifies biologically relevant domains and motions.

Reproduces known conformational states with low statistical uncertainty.

Abstract

Understanding protein conformational dynamics is essential for elucidating biological function but remains challenging due to the wide range of timescales and the complexity of collective motions. Enhanced sampling methods overcome timescale limitations of conventional molecular dynamics, yet their effectiveness depends on the choice of collective variables (CVs), which are often difficult to define and may lack physical interpretability. In particular, collective variables derived from machine learning or collective vibrational modes can efficiently capture slow dynamics but are not easily mapped onto intuitive structural descriptors. Here, we present a fully automated framework that transforms enhanced sampling trajectories into human-readable representations of protein dynamics. Our approach combines enhanced sampling along CVs derived from frequency-selective anharmonic mode analysis with a post hoc analysis of biased trajectories using weighted dynamic cross-correlation matrices. From these, we identify residue pairs and domains exhibiting correlated and anti-correlated motions, yielding simple domain-domain distances that serve as physically interpretable CVs. We apply this method to five proteins, including KRAS and HIV-1 protease, and show that it consistently identifies biologically relevant domains and motions without prior system-specific knowledge. Projection onto these distances produces free energy surfaces that reproduce known conformational states with low statistical uncertainty while maximizing independent dynamical information. This workflow enables systematic recasting of complex CVs into simple geometric descriptors without loss of essential dynamics. Its generality and automation make it broadly applicable for interpreting enhanced sampling simulations and generating interpretable conformational ensembles for integration with emerging machine learning approaches.