2025-11-25T05:19:17.586114

Nonlinear response theory of molecular machines

Chatzittofi, Agudo-Canalejo, Golestanian
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.
academic

Nonlinear response theory of molecular machines

Basic Information

  • Paper ID: 2405.14011
  • Title: Nonlinear response theory of molecular machines
  • Authors: Michalis Chatzittofi, Jaime Agudo-Canalejo, Ramin Golestanian
  • Classification: cond-mat.soft cond-mat.stat-mech physics.bio-ph
  • Publication Date: December 31, 2024 (arXiv v2)
  • Paper Link: https://arxiv.org/abs/2405.14011
  • Institutions: Max Planck Institute for Dynamics and Self-Organization, University College London, University of Oxford

Abstract

Chemical affinity drives active matter systems away from equilibrium. At the nanoscale, molecular machines interact with their surroundings and are subjected to external forces. The mechano-chemical coupling naturally arising in these systems reveals complex interactions between chemical and mechanical degrees of freedom, producing strong effects on their active mechanisms. By considering various non-equilibrium models, the authors demonstrate that modulation of applied force generates nonlinear responses, leading to non-monotonic behavior in machine activity. These findings are important for understanding, designing, and triggering such processes through controlled external fields, including collective dynamics of larger non-equilibrium systems where total dissipation and performance may be affected by internal and inter-particle interactions.

Research Background and Motivation

Core Problem

The core problem this research addresses is: How to understand the nonlinear mechanisms by which external forces affect the internal chemical reaction kinetics of molecular machines. Traditional linear response theory cannot adequately describe the complex behavior of molecular machines under external fields, particularly the non-monotonic response phenomena arising from mechano-chemical coupling.

Problem Significance

  1. Biological Significance: Biological molecular machines such as ATP hydrolases and molecular motors are fundamental to life processes; understanding their behavior under external forces is of fundamental importance to biophysics
  2. Technological Applications: External manipulation methods such as ultrasonic radiation and electromagnetic fields have been widely used in experiments and require theoretical guidance
  3. System Design: Provides theoretical foundations for designing and optimizing artificial molecular machines

Limitations of Existing Approaches

  1. Inadequacy of Linear Response Theory: Traditional theory assumes external forces affect systems linearly, unable to explain non-monotonic behavior observed in experiments
  2. Lack of Mechano-Chemical Coupling Description: Existing models often simplify activity as "active forces," ignoring the coupling between internal chemical processes and external mechanical degrees of freedom
  3. Thermodynamic Consistency Issues: Many models fail to properly handle fluctuation-dissipation relations and local detailed balance conditions

Research Motivation

The authors' research motivation is to establish a thermodynamically consistent nonlinear response theory framework that can:

  • Accurately describe how external forces affect internal reaction kinetics through mechano-chemical coupling
  • Predict and explain non-monotonic behavior observed in experiments
  • Provide theoretical guidance for the design and control of molecular machines

Core Contributions

  1. Established a nonlinear response theory framework for molecular machines: Proposed three typical models (enzyme, rotor, microswimmer) and systematically studied the nonlinear effects of external forces on internal chemical reaction flows
  2. Revealed multiple mechanisms of mechano-chemical coupling:
    • Enzyme: external force alters energy barrier height
    • Rotor: external torque modifies chemical driving force
    • Microswimmer: external field affects local detailed balance
  3. Predicted non-monotonic response behavior: Theory predicts that reaction activity reaches maximum at moderate external force strength, consistent with experimental observations
  4. Provided a thermodynamically consistent mathematical framework: Correctly handled multiplicative noise, fluctuation-dissipation relations, and entropy production

Methodology Details

Task Definition

Study the response behavior of internal chemical reaction flow JJ and diffusion coefficient DD of molecular machines under external force, where:

  • Input: external force strength, chemical affinity, system geometric parameters
  • Output: steady-state reaction flow and diffusion coefficient
  • Constraints: thermodynamic consistency, local detailed balance

Model Architecture

1. Enzyme Model

Basic Setup: Enzyme consists of two subunits with internal reaction coordinate ϕ\phi and spatial coordinates x1,x2x_1, x_2

Potential Energy Function: U(x1,x2,ϕ)=V~(ϕ)fϕ+k2(x1x2L(ϕ))2Fx1U(x_1, x_2, \phi) = \tilde{V}(\phi) - f\phi + \frac{k}{2}(x_1-x_2-L(\phi))^2 - Fx_1

where:

  • V~(ϕ)\tilde{V}(\phi): conservative part of chemical reaction
  • ff: chemical driving force, Δμ=2πf\Delta\mu = 2\pi f
  • L(ϕ)=L(0)+cos(ϕ+δ)L(\phi) = L^{(0)} + \ell\cos(\phi + \delta): conformation-dependent rest length
  • FF: applied external force

Reduced Dynamics: Through adiabatic elimination of fast variables, obtain stochastic dynamics equation for single phase ϕ\phi: ϕ˙=M(ϕ)ϕV(ϕ)+kBT2ϕM(ϕ)+2kBTM(ϕ)ξ\dot{\phi} = -M(\phi)\partial_\phi V(\phi) + \frac{k_BT}{2}\partial_\phi M(\phi) + \sqrt{2k_BT M(\phi)}\xi

Effective Potential and Mobility: V(ϕ)=fϕV~(ϕ)FL(ϕ)V(\phi) = -f\phi - \tilde{V}(\phi) - FL(\phi)M(ϕ)=mϕ1+mϕ(ϕL)2/m1M(\phi) = \frac{m_\phi}{1 + m_\phi(\partial_\phi L)^2/m_1}

2. Rotor Model

Potential Energy Function: Urot(θ,ϕ)=V~(ϕ)fϕkθcos(ϕnθ)τθU_{rot}(\theta, \phi) = \tilde{V}(\phi) - f\phi - k_\theta\cos(\phi - n\theta) - \tau\theta

Reduced Dynamics: Vrot(ϕ)=ϕ(f+τ/n)vcos(ϕ+arcsin(f/v))V_{rot}(\phi) = -\phi(f + \tau/n) - v\cos(\phi + \arcsin(f/v))

Key Feature: External torque τ\tau directly modifies the effective driving force.

3. Microswimmer Model

Four-State Markov Model: A → B → C → D → A

Transition Rates: Corrected under external electric field EEkβα=k0,βαexp(θβαWβα/kBT)k_{\beta\alpha} = k_{0,\beta\alpha}\exp(\theta_{\beta\alpha}W_{\beta\alpha}/k_B T)

where WβαW_{\beta\alpha} is work done by external field, θβα=1θαβ\theta_{\beta\alpha} = 1 - \theta_{\alpha\beta} ensures local detailed balance.

Technical Innovations

  1. Correct Treatment of Multiplicative Noise: Multiplicative noise derived through adiabatic elimination requires additional spurious drift terms to ensure thermodynamic consistency
  2. Multiscale Analysis: Utilize time scale separation (m1kmϕΔμm_1k \gg m_\phi\Delta\mu) for dimensional reduction
  3. Geometric Phase Effects: Phase difference δ\delta controls the direction and magnitude of optimal external force
  4. Charge Neutrality Constraint: Microswimmer model maintains total charge neutrality; external field produces no net drift but affects internal dynamics

Experimental Setup

Model Parameters

  • Enzyme Model: mϕ2/m1=1m_\phi\ell^2/m_1 = 1, kBT/Eba=1k_BT/E_{ba} = 1
  • Rotor Model: Standard parameter settings considering different nn values
  • Microswimmer: αL=2.1/6\alpha_L = 2.1/6, αS=2.2/6\alpha_S = 2.2/6, θβα=1/2\theta_{\beta\alpha} = 1/2

Evaluation Metrics

Reaction Flow: J=limtϕt2πtJ = \lim_{t\to\infty} \frac{\langle\phi\rangle_t}{2\pi t}

Diffusion Coefficient: D=limtϕ2tϕt2(2π)22tD = \lim_{t\to\infty} \frac{\langle\phi^2\rangle_t - \langle\phi\rangle_t^2}{(2\pi)^2 \cdot 2t}

Analysis Methods

  • Multiplicative Noise Case: Use analytical expressions (Eqs. 12-15)
  • Additive Noise Case: Standard Kramers theory
  • Discrete Markov Process: Master equation analysis

Experimental Results

Main Results

1. Enzyme Model

  • Non-monotonic Behavior: Reaction flow J(F)J(F) reaches maximum at moderate force values
  • Phase Dependence: Optimal force direction determined by phase difference δ\delta
  • Giant Diffusion Phenomenon: Diffusion coefficient exhibits bimodal structure at specific force values

Key Finding: When F=v/F = -v/\ell (at δ=arcsin(f/v)\delta = \arcsin(f/v)), oscillatory terms in potential are completely canceled, achieving global optimum.

2. Rotor Model

  • Bidirectionality: Negative torque can reverse reaction (e.g., ATP synthesis instead of hydrolysis)
  • Stall Torque: Reaction stops at τ=nf\tau = -nf
  • Linear Regime: Linear response at large torques

3. Microswimmer Model

  • Charge Distribution Effects: Dipole and quadrupole models show different optimal field directions
  • Neutrality Constraint: Zero total charge yet still capable of modulating activity
  • Swimming Speed Control: Reaction flow directly correlates with self-propulsion speed

Quantitative Results

  • Enzyme Model: Maximum enhancement reaches J(Fopt)/J(0)1.21.5J(F_{opt})/J(0) \sim 1.2-1.5
  • Rotor Model: Complete reversal achievable, J(τ)/J(0)J(\tau)/J(0) ranges from 10-10 to +10+10
  • Microswimmer: Enhancement magnitude depends on charge distribution and geometric parameters

Experimental Validation

The paper cites multiple experimental studies 28 reporting non-monotonic behavior, highly consistent with theoretical predictions.

Main Research Directions

  1. Active Matter Theory: Self-propelled particles, collective motion
  2. Molecular Motors: ATP synthase, myosin, kinesin
  3. Stochastic Thermodynamics: Fluctuation theorems, entropy production, efficiency optimization
  4. Non-equilibrium Statistical Mechanics: Markov processes, master equation theory
  • Seifert's Stochastic Thermodynamics: Provides theoretical foundation
  • Golestanian's Microswimmer Theory: Extended to external field effects
  • Jülicher's Molecular Motor Theory: Extended from kinematics to dynamics

Advantages of This Work

  1. Unified Framework: Encompasses multiple types of molecular machines
  2. Nonlinear Effects: Goes beyond linear response theory
  3. Thermodynamic Consistency: Correctly handles fluctuations and dissipation

Conclusions and Discussion

Main Conclusions

  1. Universality: All three models exhibit non-monotonic response behavior
  2. Mechanism Diversity: External forces affect internal dynamics through different mechanisms
  3. Controllability: External fields can be used to optimize molecular machine performance

Limitations

  1. Model Simplification: Uses minimal models; actual systems are more complex
  2. Parameter Dependence: Results are relatively sensitive to model parameters
  3. Experimental Validation: Requires more quantitative experimental verification

Future Directions

  1. Collective Effects: Cooperative behavior in multi-particle systems
  2. Real-time Control: Strategies for dynamically adjusting external fields
  3. Artificial Design: Molecular machine design based on theory

In-Depth Evaluation

Strengths

  1. Theoretical Innovation: Establishes a complete theoretical framework for nonlinear response of molecular machines
  2. Mathematical Rigor: Correctly handles multiplicative noise and thermodynamic consistency
  3. Physical Insight: Reveals deep mechanisms of mechano-chemical coupling
  4. Practical Value: Provides guidance for experimental design and applications

Weaknesses

  1. Model Limitations: Oversimplified geometry and interactions
  2. Parameter Fitting: Lacks quantitative comparison with specific experimental systems
  3. Dynamic Response: Does not consider transient behavior under time-varying external fields

Impact

  1. Academic Contribution: Provides new analytical tools for molecular machine theory
  2. Application Prospects: Significant value in biophysics and nanotechnology
  3. Reproducibility: Clear theoretical models, easy to reproduce and extend

Applicable Scenarios

  1. Enzyme Engineering: Optimize enzyme activity and selectivity
  2. Molecular Motor Design: Artificial molecular machine development
  3. Drug Delivery: Control microswimmer motion
  4. Biosensing: Detection based on conformational changes

References

The paper cites 39 important references covering core works in active matter, molecular motors, stochastic thermodynamics and related fields, reflecting the academic depth and breadth of this research.


Overall Assessment: This is a high-quality theoretical physics paper that makes important contributions to nonlinear response theory of molecular machines. The mathematical derivations are rigorous, the physical picture is clear, the predicted results are consistent with experimental observations, and it establishes an important foundation for further development in this field.