Abstract:A modeling framework for the internal conformational dynamics and external mechanical movement of single biological macromolecules in aqueous solution at constant temperature is developed. Both the internal dynamics and external movement are stochastic; the former is represented by a master equation for a set of discrete states, and the latter is described by a continuous Smoluchowski equation. Combining these two equations into one, a comprehensive theory for the Brownian dynamics and statistical thermodynamics of single macromolecules arises. This approach is shown to have wide applications. It is applied to protein-ligand dissociation under external force, unfolding of polyglobular proteins under extension, movement along linear tracks of motor proteins against load, and enzyme catalysis by single fluctuating proteins. As a generalization of the classic polymer theory, the dynamic equation is capable of characterizing a single macromolecule in aqueous solution, in probabilistic terms, (1) its thermodynamic equilibrium with fluctuations, (2) transient relaxation kinetics, and most importantly and novel (3) nonequilibrium steady-state with heat dissipation. A reversibility condition which guarantees an equilibrium solution and its thermodynamic stability is established, an H-theorem like inequality for irreversibility is obtained, and a rule for thermodynamic consistency in chemically pumped nonequilibrium steady-state is given.
What problem does this paper attempt to address?
The core problem that this paper attempts to solve is to develop a theoretical framework that can describe the internal conformational dynamics and external mechanical motion of a single biomacromolecule in a constant - temperature aqueous solution. Specifically, the author aims to combine the discrete - state Markov model and the continuous Brownian motion model to propose a unified Stochastic Macromolecular Mechanics (SM3) theory. This theory needs to satisfy thermodynamic consistency and be able to integrate quantitative measurement data at the macroscopic, mesoscopic, and microscopic levels.
### Main contributions of the paper
1. **Unified theoretical framework**: By combining internal conformational dynamics (represented by the master equation) with external mechanical motion (represented by the Smoluchowski equation), a comprehensive theoretical framework is proposed.
2. **Wide range of applications**: This theory can be applied to processes such as protein - ligand dissociation, the unfolding of multi - globular proteins, the movement of kinesin along the track, and the enzymatic catalysis of a single fluctuating protein.
3. **Nonequilibrium steady - state**: Special emphasis is placed on the study of the nonequilibrium steady - state, which is a new area not covered by classical polymer theory. The irreversible process caused by chemical pumping (such as ATP hydrolysis) is introduced, and its thermodynamic characteristics are explored.
4. **Thermodynamic consistency**: The conditions for ensuring thermodynamic stability, namely the detailed balance condition, are established, and the thermodynamic consistency rules for the nonequilibrium steady - state under chemical pumping are given.
### Summary of mathematical formulas
- **Smoluchowski equation**:
\[
\beta\frac{\partial P(x,t)}{\partial t}=k_BT\frac{\partial^2P(x,t)}{\partial x^2}-\frac{\partial}{\partial x}[(F_{\text{int}}(x)+F_{\text{ext}}(x,t))P(x,t)]
\]
where \(P(x,t)\) is the probability density of the ligand at position \(x\), \(k_B\) is the Boltzmann constant, \(T\) is the temperature, and \(F_{\text{int}}(x)\) and \(F_{\text{ext}}(x,t)\) are the internal and external forces respectively.
- **Master equation**:
\[
\frac{\partial P(n,t)}{\partial t}=(N - n + 1)\lambda_u(x,n - 1)P(n - 1,t)+(n + 1)\lambda_f(x,n + 1)P(n + 1,t)-[n\lambda_f(x,n)+(N - n)\lambda_u(x,n)]P(n,t)
\]
where \(P(n,t)\) is the probability of being in state \(n\) at time \(t\), and \(\lambda_f\) and \(\lambda_u\) are the folding and unfolding rates respectively.
- **Evolution equation of the joint probability distribution**:
\[
\frac{\partial P(x,n,t)}{\partial t}=\left(\frac{k_BT}{\beta}\right)\frac{\partial^2P(x,n,t)}{\partial x^2}+\frac{1}{\beta}\frac{\partial}{\partial x}\left[\left(\frac{\partial U_{\text{int}}(x,n)}{\partial x}+F_{\text{ext}}\right)P(x,n,t)\right]+\sum_{k = 1}^N[\lambda_{kn}(x)P(x,k,t)-\lambda_{nk}(x)P(x,n,
\]