Mathematical Modeling of the Dynamic Mechanical Behavior of Neighboring Sarcomeres in Actin Stress Fibers

Actin stress fibers (SFs) in live cells consist of series of dynamic individual sarcomeric units. Within a group of consecutive SF sarcomeres, individual sarcomeres can spontaneously shorten or lengthen without changing the overall length of this group, but the underlying mechanism is unclear. We used a computational model to test our hypothesis that this dynamic behavior is inherent to the heterogeneous mechanical properties of the sarcomeres and the cytoplasmic viscosity. Each sarcomere was modeled as a discrete element consisting of an elastic spring, a viscous dashpot and an active contractile unit all connected in parallel, and experiences forces as a result of actin filament elastic stiffness, myosin II contractility, internal viscoelasticity, or cytoplasmic drag. When all four types of forces are considered, the simulated dynamic behavior closely resembles the experimental observations, which include a low-frequency fluctuation in individual sarcomere length and compensatory lengthening and shortening of adjacent sarcomeres. Our results suggest that heterogeneous stiffness and viscoelasticity of actin fibers, heterogeneous myosin II contractility, and the cytoplasmic drag are sufficient to cause spontaneous fluctuations in SF sarcomere length. Our results shed new light to the dynamic behavior of SF and help design experiments to further our understanding of SF dynamics.