Molecular Dynamics (MD) simulations have been widely used to study the dynamics and mechanistic processes of biomolecules, including protein folding , protein-ligand binding , and large-scale functional conformational changes . However, MD simulations generate time sequences of high-dimensional protein conformations each containing the Cartesian coordinates of all protein atoms, posing a significant challenge for the analysis and interpretation of biomolecular MD simulation datasets.