Chemical affinities are responsible for driving active matter systems out of equilibrium. At the nano-scale, molecular machines interact with the surrounding environment and are subjected to external forces. The mechano-chemical coupling which arises naturally in these systems reveals a complex interplay between chemical and mechanical degrees of freedom with strong impact on their active mechanism. By considering various models far from equilibrium, we show that the tuning of applied forces give rise to a nonlinear response that causes a non-monotonic behaviour in the machines' activity. Our findings have implications in understanding, designing, and triggering such processes by controlled application of external fields, including the collective dynamics of larger non-equilibrium systems where the total dissipation and performance might be affected by internal and inter-particle interactions.
论文ID : 2405.14011标题 : Nonlinear response theory of molecular machines作者 : Michalis Chatzittofi, Jaime Agudo-Canalejo, Ramin Golestanian分类 : cond-mat.soft cond-mat.stat-mech physics.bio-ph发表时间 : 2024年12月31日 (arXiv v2)论文链接 : https://arxiv.org/abs/2405.14011 机构 : Max Planck Institute for Dynamics and Self-Organization, University College London, University of Oxford化学亲和力负责驱动活性物质系统偏离平衡态。在纳米尺度上,分子机器与周围环境相互作用并受到外力作用。这些系统中自然产生的机械-化学耦合揭示了化学和机械自由度之间的复杂相互作用,对其活性机制产生强烈影响。通过考虑各种远离平衡的模型,作者证明了施加力的调节会产生非线性响应,导致机器活性的非单调行为。这些发现对理解、设计和通过受控外场应用触发此类过程具有重要意义,包括更大非平衡系统的集体动力学,其中总耗散和性能可能受到内部和粒子间相互作用的影响。
该研究要解决的核心问题是:如何理解外力对分子机器内部化学反应动力学的非线性影响机制 。传统的线性响应理论无法充分描述在外场作用下分子机器的复杂行为,特别是机械-化学耦合产生的非单调响应现象。
生物学意义 :ATP水解酶、分子马达等生物分子机器是生命活动的基础,理解其在外力作用下的行为对生物物理学具有根本意义技术应用 :超声辐射、电磁场等外部操控方法在实验中已被广泛使用,需要理论指导系统设计 :为设计和优化人工分子机器提供理论基础线性响应理论的不足 :传统理论假设外力对系统的影响是线性的,无法解释实验中观察到的非单调行为缺乏机械-化学耦合描述 :现有模型往往将活性简化为"活性力",忽略了内部化学过程与外部机械自由度的耦合热力学一致性问题 :许多模型未能正确处理涨落-耗散关系和局部详细平衡条件作者的研究动机是建立一个热力学一致的非线性响应理论框架 ,能够:
准确描述外力如何通过机械-化学耦合影响内部反应动力学 预测和解释实验中观察到的非单调行为 为分子机器的设计和控制提供理论指导 建立了分子机器的非线性响应理论框架 :提出了三个典型模型(酶、转子、微游泳器),系统研究了外力对内部化学反应流的非线性影响揭示了机械-化学耦合的多种机制 :酶:外力改变能垒高度 转子:外扭矩改变化学驱动力 微游泳器:外场影响局部详细平衡 预测了非单调响应行为 :理论预测在中等强度外力下反应活性达到最大值,与实验观察一致提供了热力学一致的数学框架 :正确处理了乘性噪声、涨落-耗散关系和熵产生研究外力作用下分子机器内部化学反应流 J J J 和扩散系数 D D D 的响应行为,其中:
输入 :外力强度、化学亲和力、系统几何参数输出 :稳态反应流和扩散系数约束 :热力学一致性、局部详细平衡基本设置 :酶由两个亚基组成,具有内部反应坐标 ϕ \phi ϕ 和空间坐标 x 1 , x 2 x_1, x_2 x 1 , x 2
势能函数 :
U ( x 1 , x 2 , ϕ ) = V ~ ( ϕ ) − f ϕ + k 2 ( x 1 − x 2 − L ( ϕ ) ) 2 − F x 1 U(x_1, x_2, \phi) = \tilde{V}(\phi) - f\phi + \frac{k}{2}(x_1-x_2-L(\phi))^2 - Fx_1 U ( x 1 , x 2 , ϕ ) = V ~ ( ϕ ) − f ϕ + 2 k ( x 1 − x 2 − L ( ϕ ) ) 2 − F x 1
其中:
V ~ ( ϕ ) \tilde{V}(\phi) V ~ ( ϕ ) :化学反应的保守部分f f f :化学驱动力,Δ μ = 2 π f \Delta\mu = 2\pi f Δ μ = 2 π f L ( ϕ ) = L ( 0 ) + ℓ cos ( ϕ + δ ) L(\phi) = L^{(0)} + \ell\cos(\phi + \delta) L ( ϕ ) = L ( 0 ) + ℓ cos ( ϕ + δ ) :构象依赖的休息长度F F F :外加力约化动力学 :通过绝热消除快变量,得到单一相位 ϕ \phi ϕ 的随机动力学方程:
ϕ ˙ = − M ( ϕ ) ∂ ϕ V ( ϕ ) + k B T 2 ∂ ϕ M ( ϕ ) + 2 k B T M ( ϕ ) ξ \dot{\phi} = -M(\phi)\partial_\phi V(\phi) + \frac{k_BT}{2}\partial_\phi M(\phi) + \sqrt{2k_BT M(\phi)}\xi ϕ ˙ = − M ( ϕ ) ∂ ϕ V ( ϕ ) + 2 k B T ∂ ϕ M ( ϕ ) + 2 k B TM ( ϕ ) ξ
有效势和迁移率 :
V ( ϕ ) = − f ϕ − V ~ ( ϕ ) − F L ( ϕ ) V(\phi) = -f\phi - \tilde{V}(\phi) - FL(\phi) V ( ϕ ) = − f ϕ − V ~ ( ϕ ) − F L ( ϕ ) M ( ϕ ) = m ϕ 1 + m ϕ ( ∂ ϕ L ) 2 / m 1 M(\phi) = \frac{m_\phi}{1 + m_\phi(\partial_\phi L)^2/m_1} M ( ϕ ) = 1 + m ϕ ( ∂ ϕ L ) 2 / m 1 m ϕ
势能函数 :
U r o t ( θ , ϕ ) = V ~ ( ϕ ) − f ϕ − k θ cos ( ϕ − n θ ) − τ θ U_{rot}(\theta, \phi) = \tilde{V}(\phi) - f\phi - k_\theta\cos(\phi - n\theta) - \tau\theta U ro t ( θ , ϕ ) = V ~ ( ϕ ) − f ϕ − k θ cos ( ϕ − n θ ) − τ θ
约化动力学 :
V r o t ( ϕ ) = − ϕ ( f + τ / n ) − v cos ( ϕ + arcsin ( f / v ) ) V_{rot}(\phi) = -\phi(f + \tau/n) - v\cos(\phi + \arcsin(f/v)) V ro t ( ϕ ) = − ϕ ( f + τ / n ) − v cos ( ϕ + arcsin ( f / v ))
关键特征:外扭矩 τ \tau τ 直接修改有效驱动力。
四态马尔可夫模型 :A → B → C → D → A
转移速率 :考虑外电场 E E E 作用下的修正
k β α = k 0 , β α exp ( θ β α W β α / k B T ) k_{\beta\alpha} = k_{0,\beta\alpha}\exp(\theta_{\beta\alpha}W_{\beta\alpha}/k_B T) k β α = k 0 , β α exp ( θ β α W β α / k B T )
其中 W β α W_{\beta\alpha} W β α 是外场做功,θ β α = 1 − θ α β \theta_{\beta\alpha} = 1 - \theta_{\alpha\beta} θ β α = 1 − θ α β 保证局部详细平衡。
乘性噪声的正确处理 :通过绝热消除导出的乘性噪声需要额外的虚假漂移项来保证热力学一致性多尺度分析 :利用时间尺度分离(m 1 k ≫ m ϕ Δ μ m_1k \gg m_\phi\Delta\mu m 1 k ≫ m ϕ Δ μ )进行维数约化几何相位效应 :相位差 δ \delta δ 控制最优外力的方向和大小电荷中性约束 :微游泳器模型保持总电荷中性,外场不产生净漂移但影响内部动力学酶模型 :m ϕ ℓ 2 / m 1 = 1 m_\phi\ell^2/m_1 = 1 m ϕ ℓ 2 / m 1 = 1 ,k B T / E b a = 1 k_BT/E_{ba} = 1 k B T / E ba = 1 转子模型 :标准参数设置,考虑不同的 n n n 值微游泳器 :α L = 2.1 / 6 \alpha_L = 2.1/6 α L = 2.1/6 ,α S = 2.2 / 6 \alpha_S = 2.2/6 α S = 2.2/6 ,θ β α = 1 / 2 \theta_{\beta\alpha} = 1/2 θ β α = 1/2 反应流 :
J = lim t → ∞ ⟨ ϕ ⟩ t 2 π t J = \lim_{t\to\infty} \frac{\langle\phi\rangle_t}{2\pi t} J = lim t → ∞ 2 π t ⟨ ϕ ⟩ t
扩散系数 :
D = lim t → ∞ ⟨ ϕ 2 ⟩ t − ⟨ ϕ ⟩ t 2 ( 2 π ) 2 ⋅ 2 t D = \lim_{t\to\infty} \frac{\langle\phi^2\rangle_t - \langle\phi\rangle_t^2}{(2\pi)^2 \cdot 2t} D = lim t → ∞ ( 2 π ) 2 ⋅ 2 t ⟨ ϕ 2 ⟩ t − ⟨ ϕ ⟩ t 2
乘性噪声情况 :使用解析表达式 (Eqs. 12-15)加性噪声情况 :标准Kramers理论离散马尔可夫过程 :主方程分析非单调行为 :反应流 J ( F ) J(F) J ( F ) 在中等力值处达到最大值相位依赖性 :最优力的方向由相位差 δ \delta δ 决定巨扩散现象 :扩散系数在特定力值处出现双峰结构关键发现 :当 F = − v / ℓ F = -v/\ell F = − v / ℓ (δ = arcsin ( f / v ) \delta = \arcsin(f/v) δ = arcsin ( f / v ) 时),势能中的振荡项被完全抵消,达到全局最优。
双向性 :负扭矩可使反应逆转(如ATP合成代替水解)失速扭矩 :τ = − n f \tau = -nf τ = − n f 时反应停止线性区域 :大扭矩时呈现线性响应电荷分布效应 :偶极和四极模型显示不同的最优场方向中性约束 :总电荷为零但仍能调节活性游泳速度调控 :反应流直接关联自推进速度酶模型:最大增强可达 J ( F o p t ) / J ( 0 ) ∼ 1.2 − 1.5 J(F_{opt})/J(0) \sim 1.2-1.5 J ( F o pt ) / J ( 0 ) ∼ 1.2 − 1.5 转子模型:可实现完全逆转,J ( τ ) / J ( 0 ) J(\tau)/J(0) J ( τ ) / J ( 0 ) 从 − 10 -10 − 10 到 + 10 +10 + 10 微游泳器:增强幅度依赖于电荷分布和几何参数 论文引用了多项实验研究 28 报告的非单调行为,与理论预测高度一致。
活性物质理论 :自推进粒子、集体运动分子马达 :ATP合酶、肌球蛋白、驱动蛋白随机热力学 :涨落定理、熵产生、效率优化非平衡统计力学 :马尔可夫过程、主方程理论Seifert等的随机热力学 :提供了理论基础Golestanian的微游泳器理论 :扩展到外场效应Jülicher的分子马达理论 :从运动学扩展到动力学统一框架 :涵盖多种分子机器类型非线性效应 :超越线性响应理论热力学一致性 :正确处理涨落和耗散普遍性 :所有三种模型都展现非单调响应行为机制多样性 :外力通过不同机制影响内部动力学可控性 :外场可用于优化分子机器性能模型简化 :使用最小模型,实际系统更复杂参数依赖 :结果对模型参数较为敏感实验验证 :需要更多定量实验验证集体效应 :多粒子系统的协同行为实时控制 :动态调节外场的策略人工设计 :基于理论的分子机器设计理论创新 :建立了分子机器非线性响应的完整理论框架数学严谨 :正确处理了乘性噪声和热力学一致性物理洞察 :揭示了机械-化学耦合的深层机制实用价值 :为实验设计和应用提供指导模型局限 :过于简化的几何和相互作用参数拟合 :缺乏与具体实验系统的定量对比动态响应 :未考虑时变外场的瞬态行为学术贡献 :为分子机器理论提供新的分析工具应用前景 :在生物物理和纳米技术领域具有重要价值可复现性 :理论模型清晰,易于重现和扩展酶工程 :优化酶活性和选择性分子马达设计 :人工分子机器开发药物递送 :控制微游泳器的运动生物传感 :基于构象变化的检测论文引用了39篇重要文献,涵盖了活性物质、分子马达、随机热力学等相关领域的核心工作,体现了该研究的学术深度和广度。
总体评价 :这是一篇高质量的理论物理论文,在分子机器的非线性响应理论方面做出了重要贡献。论文的数学推导严谨,物理图像清晰,预测结果与实验观察一致,为该领域的进一步发展奠定了重要基础。