A kinematic model coupling stress fiber dynamics with JNK activation in response to matrix stretching

The role of the actin cytoskeleton in regulating mechanotransduction in response to external forces is complex and incompletely understood. Here, we develop a mathematical model coupling the dynamic disassembly and reassembly of actin stress fibers and associated focal adhesions to the activation of c-jun N-terminal kinase (JNK) in cells attached to deformable matrices. The model is based on the assumptions that stress fibers are pre-extended to a preferred level under static conditions and that perturbations from this preferred level destabilize the stress fibers. The subsequent reassembly of fibers upregulates the rate of JNK activation as a result of the formation of new integrin bonds within the associated focal adhesions. Numerical solutions of the model equations predict that different patterns of matrix stretch result in distinct temporal patterns in JNK activation that compare well with published experimental results. In the case of cyclic uniaxial stretching, stretch-induced JNK activation slowly subsides as stress fibers gradually reorient perpendicular to the stretch direction. In contrast, JNK activation is chronically elevated in response to cyclic equibiaxial stretch. A step change in either uniaxial or equibiaxial stretch results in a short, transient upregulation in JNK that quickly returns to the basal level as overly stretched stress fibers disassemble and are replaced by fibers assembled at the preferred level of stretch. In summary, the model describes a mechanism by which the dynamic properties of the actin cytoskeleton allow cells to adapt to applied forces through turnover and reorganization to modulate intracellular signaling.