arXiv:2409.06428physics.chem-phcond-mat.stat-mech2024-09中稿 · as part of J被引 9

用谱映射法学习蛋白折叠的慢变量,构建马尔可夫动力学模型。

Spectral Map for Slow Collective Variables, Markovian Dynamics, and Transition State Ensembles

  • 通过最大化转移矩阵谱隙,学习能捕捉慢动态的集体变量。
  • 单个慢变量即可作为物理反应坐标,准确描述蛋白折叠特征。
  • 适用于研究蛋白质折叠机制,尤其适合关注动力学路径的科研人员。

理解复杂分子系统的演化是物理化学的核心问题。为描述其长时动力学——决定系统关键特性的部分——可识别少数慢集体变量(CVs),并将其余快变量视为热噪声,从而将动力学简化为慢变量空间中的自由能景观扩散,实现马尔可夫化。我们先前提出的谱映射方法(Rydzewski, J. Phys. Chem. Lett. 2023, 14, 22, 5216-5220)通过最大化转移矩阵的谱隙来学习慢CVs。本文在此框架中引入若干改进,以蛋白质高维可逆折叠过程为例:提出一种粗粒化马尔可夫转移矩阵的算法,基于动力学对慢变量空间进行分区,并定义过渡态集合;结果表明,谱映射学习的慢变量已接近过阻尼扩散的马尔可夫极限;坐标依赖的扩散系数对构建的自由能景观影响微弱;此外,展示了谱映射用于量化特征重要性,并与常见的结构描述符比较慢变量表现。总体而言,单一谱映射学习的慢变量可作为物理反应坐标,有效捕捉蛋白折叠的本质特征。

原文摘要 · Abstract (English)

Understanding the behavior of complex molecular systems is a fundamental problem in physical chemistry. To describe the long-time dynamics of such systems, which is responsible for their most informative characteristics, we can identify a few slow collective variables (CVs) while treating the remaining fast variables as thermal noise. This enables us to simplify the dynamics and treat it as diffusion in a free-energy landscape spanned by slow CVs, effectively rendering the dynamics Markovian. Our recent statistical learning technique, spectral map [Rydzewski, J. Phys. Chem. Lett. 2023, 14, 22, 5216-5220], explores this strategy to learn slow CVs by maximizing a spectral gap of a transition matrix. In this work, we introduce several advancements into our framework, using a high-dimensional reversible folding process of a protein as an example. We implement an algorithm for coarse-graining Markov transition matrices to partition the reduced space of slow CVs kinetically and use it to define a transition state ensemble. We show that slow CVs learned by spectral map closely approach the Markovian limit for an overdamped diffusion. We demonstrate that coordinate-dependent diffusion coefficients only slightly affect the constructed free-energy landscapes. Finally, we present how spectral map can be used to quantify the importance of features and compare slow CVs with structural descriptors commonly used in protein folding. Overall, we demonstrate that a single slow CV learned by spectral map can be used as a physical reaction coordinate to capture essential characteristics of protein folding.

分子动力学蛋白折叠慢变量马尔可夫模型

Thank you to arXiv for use of its open access interoperability. PaperDance 不是 arXiv 官方产品;中文卡片由大模型生成,请以原文为准。