A. Fiasconaro,J.J. Mazo,F. Falo
Abstract:In this work we study the assisted translocation of a polymer across a membrane nanopore, inside which a molecular motor exerts a force fuelled by the hydrolysis of ATP molecules. In our model the motor switches to its active state for a fixed amount of time, while it waits for an ATP molecule binding and triggering the impulse, during an exponentially distributed time lapse. The polymer is modelled as a beads-springs chain with both excluded volume and bending contributions, and moves in a stochastic three dimensional environment modelled with a Langevin dynamics at fixed temperature. The resulting dynamics shows a Michaelis-Menten translocation velocity that depends on the chain flexibility. The scaling behavior of the mean translocation time with the polymer length for different bending values is also investigated.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to study the transmembrane transport process of semi - flexible polymers assisted by ATP - driven molecular motors in nanopores. Specifically, the paper focuses on:
1. **Power mechanism**: Explore how the hydrolysis of ATP molecules provides energy for molecular motors, thereby promoting polymers to pass through nanopores.
2. **Kinetic behavior**: Analyze the dynamic behavior of the average transport speed and transport time of polymers under different conditions (such as polymer length, bending stiffness, molecular motor activation rate, etc.).
3. **Michaelis - Menten kinetics**: Verify and explain the Michaelis - Menten relationship between the transport speed and the molecular motor activation rate.
4. **Influencing factors**: Study the influence of the flexibility of polymers and the activity of molecular motors on the transport process, especially how they affect the transport time and speed.
### Formula summary
- **Elastic potential energy**:
\[
V_{\text{el}}(d_i)=\frac{k_e}{2}\sum_{i = 1}^N(d_i - d_0)^2
\]
where \(k_e\) is the elastic parameter, \(d_i=|r_{i + 1}-r_i|\) is the distance between the \(i\)th and \((i + 1)\)th monomers, and \(d_0\) is the equilibrium distance between adjacent monomers.
- **Bending potential energy**:
\[
V_{\text{ben}}(\theta_i)=\frac{k_b}{2}\sum_{i = 1}^N[1-\cos(\theta_i-\theta_0)]
\]
where \(k_b\) is the bending elastic constant, \(\theta_i\) is the angle between the \(i\)th and \((i + 1)\)th monomers, and \(\theta_0 = 0\) is the equilibrium angle.
- **Lennard - Jones potential energy**:
\[
V_{\text{LJ}}(r_{ij})=4\epsilon\sum_{i\neq j}^N\left[\left(\frac{\sigma}{r_{ij}}\right)^{12}-\left(\frac{\sigma}{r_{ij}}\right)^6\right]
\]
For \(r_{ij}\leq2^{1/6}\sigma\), otherwise \(-\epsilon\). Where \(r_{ij}\) is the distance between the \(i\)th and \(j\)th monomers.
- **Equation of motion**:
\[
m\gamma\dot{r}_i=-\nabla_i V_{\text{el}}(d_i)-\nabla_i V_{\text{ben}}(\theta_i)-\nabla_i V_{\text{LJ}}(r_{ij})+F_{\text{drv},i}+F_{\text{sp},i}+\sqrt{2m\gamma k_B T}\xi_i(t)
\]
where \(\xi_i(t)\) is a Gaussian uncorrelated thermal fluctuation, satisfying the statistical properties \(\langle\xi_i,\alpha(t)\rangle = 0\) and \(\langle\xi_i,\alpha(t)\xi_j,\beta(t')\rangle=\delta_{ij}\delta_{\alpha,\beta}\delta(t' - t)\).
- **Driving force**:
\[
F_{\text{drv}}(x)=
\begin{cases}
(F_c,F_c + F_M)&\text{if }x\in[0,L_