Motor Clutch Model

Motor-Clutch Model

Overview

The motor-clutch model describes how a cell converts substrate stiffness into mechanical traction.

In the mechanogenomic virtual cell, this module represents the cell-substrate interface. Myosin motors pull the actin cytoskeleton, while integrin-based molecular clutches transmit force to the extracellular matrix.

The main output of this module is the stiffness-dependent traction:

$$T(E)$$

where E is the Young's modulus of the hydrogel or tissue.

This output is passed to the nuclear mechanics module as the mechanical input that drives nuclear deformation and mechanotranscription.

---

Biological interpretation

Cells do not sense stiffness as an abstract material property. Instead, they generate forces through actomyosin contraction and probe how much the substrate resists deformation.

The motor-clutch system includes:

• Myosin motors, which pull actin filaments.
• Actin retrograde flow, which loads molecular adhesions.
• Integrin clutches, which connect actin to the extracellular matrix.
• An elastic substrate, whose stiffness controls force transmission.

On soft substrates, the substrate deforms easily and little force is transmitted.

On stiff substrates, clutches load more efficiently and cellular traction increases.

At very high stiffness, clutches may load and fail rapidly, producing nonlinear or biphasic behavior.

---

Stiffness mapping

The experimentally controlled stiffness is the Young's modulus, E, given in kPa.

The motor-clutch model uses an effective clutch-scale substrate stiffness, kappa:

$$\kappa = \alpha E$$

Symbol	Meaning
E	Young's modulus of the hydrogel or tissue
kappa	Effective substrate stiffness perceived by clutches
alpha	Coupling parameter between macroscopic stiffness and clutch-scale stiffness

This mapping connects hydrogel or fibrosis-like tissue stiffness to the computational motor-clutch module.

---

Myosin motor force

The total stall force generated by the myosin motor ensemble is:

$$F_{stall} = n_m F_m$$

Symbol	Meaning
n_m	Effective number of myosin motors
F_m	Force generated by each motor
F_stall	Maximum force before actin flow stalls

The actin retrograde flow velocity decreases as substrate force increases:

$$v = v_u \left(1 - \frac{F_{sub}}{F_{stall}}\right)_+$$

Symbol	Meaning
v	Instantaneous actin retrograde velocity
v_u	Unloaded actin velocity
F_sub	Force transmitted to the substrate
(.)_+	Rectification operator; negative values are set to zero

As the transmitted force approaches the stall force, actin flow slows down.

---

Molecular clutches

Each integrin-based molecular clutch is modeled as an elastic spring.

For a bound clutch i, the force is:

$$F_i = k_c x_i$$

Symbol	Meaning
F_i	Force carried by clutch i
k_c	Clutch stiffness
x_i	Extension of clutch i

The total force transmitted to the substrate is:

$$F_{sub} = \sum_{i \in bound} k_c x_i$$

Only bound clutches contribute to force transmission.

---

Clutch binding and unbinding

Unbound clutches attach with a constant binding rate:

$$k_{on}$$

Bound clutches detach with a force-dependent unbinding rate. In this implementation, clutches behave as slip bonds:

$$k_{off} = k_{off}^{0} \exp\left(\frac{k_c |x_i|}{F_b}\right)$$

Symbol	Meaning
k_off^0	Basal clutch unbinding rate
F_b	Characteristic bond rupture force
`k_c	x_i

This means that the probability of clutch detachment increases with force.

---

Displacement partition

When actin moves, displacement is partitioned between clutch extension and substrate deformation.

The clutch extension evolves as:

$$\frac{dx_i}{dt} = v \frac{\kappa}{\kappa + k_c}$$

This term captures the relative compliance between the substrate and the clutch.

If the substrate is soft, much of the displacement is absorbed by the substrate.

If the substrate is stiff, clutches are loaded more efficiently.

---

Traction output

The main output of the motor-clutch module is the mean traction:

$$T(E) = \langle F_{sub} \rangle$$

This value is computed by stochastic simulation of clutch binding, loading, and unbinding events.

The simulation follows these steps:

1. Initialize a population of molecular clutches.
2. Allow unbound clutches to bind with rate k_on.
3. Load bound clutches through actin retrograde flow.
4. Detach bound clutches according to force-dependent k_off.
5. Compute the total substrate force F_sub.
6. Average F_sub after a burn-in period.

---

Expected stiffness response

The motor-clutch model can produce nonlinear stiffness sensing.

Regime	Behavior
Very soft substrate	Low force transmission because the substrate deforms easily
Intermediate stiffness	Efficient clutch loading and increased traction
Very stiff substrate	Rapid clutch loading and possible load-and-fail behavior

In the hepatocyte virtual-cell model, this module provides the first physical step linking hydrogel or tissue stiffness to intracellular mechanical signaling.

---

Role in the mechanogenomic virtual cell

The motor-clutch module converts substrate stiffness into cellular traction:

$$E \rightarrow T(E)$$

This traction is then transmitted to the nucleus through the cytoskeleton and LINC-associated mechanics:

$$T(E) \rightarrow \sigma_{nuc}(E)$$

Thus, the motor-clutch model acts as the mechanical engine of the virtual cell.