Computational Statistics and Machine Learning, Italian Institute of Technology, Via Enrico Melen 83, 16142 Genoa, Italy.
Atomistic Simulations, Italian Institute of Technology, Via Enrico Melen 83, 16142 Genoa, Italy.
J Chem Theory Comput. 2022 Sep 13;18(9):5195-5202. doi: 10.1021/acs.jctc.2c00393. Epub 2022 Aug 3.
Present-day atomistic simulations generate long trajectories of ever more complex systems. Analyzing these data, discovering metastable states, and uncovering their nature are becoming increasingly challenging. In this paper, we first use the variational approach to conformation dynamics to discover the slowest dynamical modes of the simulations. This allows the different metastable states of the system to be located and organized hierarchically. The physical descriptors that characterize metastable states are discovered by means of a machine learning method. We show in the cases of two proteins, chignolin and bovine pancreatic trypsin inhibitor, how such analysis can be effortlessly performed in a matter of seconds. Another strength of our approach is that it can be applied to the analysis of both unbiased and biased simulations.
目前的原子模拟生成了越来越复杂的系统的长轨迹。分析这些数据、发现亚稳状态并揭示它们的本质变得越来越具有挑战性。在本文中,我们首先使用构象动力学的变分方法来发现模拟的最慢动力学模式。这允许系统的不同亚稳状态被定位并分层组织。通过机器学习方法发现了描述亚稳状态的物理描述符。我们展示了在 chignolin 和牛胰蛋白酶抑制剂这两种蛋白质的情况下,这种分析可以在几秒钟内毫不费力地完成。我们方法的另一个优势是它可以应用于无偏和有偏模拟的分析。