Towards multi-scale modeling of muscle fibers with sarcomere non-uniformities

Solitary waves in a nonintegrable chain with double-well potentials

We study solitary waves in a one-dimensional lattice of identical masses that are connected in series by nonlinear springs. The potential of each spring is nonconvex, where two disjoint convex regions, phase I and phase II, are separated by a concave, spinodal region. Consequently, the force-strain relation of the spring is nonmonotonous, which gives rise to a bistable behavior. Based on analytical treatment, with some approximations, combined with extensive numerical simulations, we are able to reveal important insights. For example, we find that the solitary-wave solution is indifferent to the energy barrier that separates the two energy wells associated with phase I and phase II, and that the shape of the wave can be described by means of merely two scalar properties of the potential of the springs, namely, the ratio of stiffness in phase II and phase I, and the ratio between the Maxwell's force and corresponding transition strain. The latter ratio provides a useful measure for the significance of the spinodal region. Linear stability of the solitary-wave solution is studied analytically using the Vakhitov-Kolokolov criterion applied to the approximate solutions obtained in the first part. These results are validated by numerical simulations. We find that the solitary-wave solution is stable provided that its velocity is higher than some critical value. It is shown that, practically, the solitary waves are stable for almost the entire range of possible wave velocities. This is also manifested in the interaction between two solitary waves or between a solitary wave and a wall (rigid boundary). Such interaction results in a minor change of height and shape of the solitary wave along with the formation of a trail of small undulations that follow the wave, as expected in a nonintegrable system. Even after a significant number of interactions the changes in the wave height and shape are minor, suggesting that the bistable chain may be a useful platform for delivering information over long distances, even concurrently with additional information (other solitary waves) passing through the chain.

Curvature-Induced Spatial Ordering of Composition in Lipid Membranes

Phase segregation of membranal components, such as proteins, lipids, and cholesterols, leads to the formation of aggregates or domains that are rich in specific constituents. This process is important in the interaction of the cell with its surroundings and in determining the cell's behavior and fate. Motivated by published experiments on curvature-modulated phase separation in lipid membranes, we formulate a mathematical model aiming at studying the spatial ordering of composition in a two-component biomembrane that is subjected to a prescribed (imposed) geometry. Based on this model, we identified key nondimensional quantities that govern the biomembrane response and performed numerical simulations to quantitatively explore their influence. We reproduce published experimental observations and extend them to surfaces with geometric features (imposed geometry) and lipid phases beyond those used in the experiments. In addition, we demonstrate the possibility for curvature-modulated phase separation above the critical temperature and propose a systematic procedure to determine which mechanism, the difference in bending stiffness or difference in spontaneous curvatures of the two phases, dominates the coupling between shape and composition.

Contraction induced muscle injury: towards personalized training and recovery programs

Skeletal muscles can be injured by their own contractions. Such contraction-induced injury, often accompanied by delayed onset of muscle soreness, is a leading cause of the loss of mobility in the rapidly increasing population of elderly people. Unlike other types of muscle injuries which hurt almost exclusively those who are subjected to intensive exercise such as professional athletes and soldiers in training, contraction induced injury is a phenomenon which may be experienced by people of all ages while performing a variety of daily-life activities. Subjects that experience contraction induced injury report on soreness that usually increases in intensity in the first 24 h after the activity, peaks from 24 to 72 h, and then subsides and disappears in a few days. Despite their clinical importance and wide influence, there are almost no studies, clinical, experimental or computational, that quantitatively relate between the extent of contraction induced injury and activity factors, such as number of repetitions, their frequency and magnitude. The lack of such quantitative information is even more emphasized by the fact that contraction induced injury can be used, if moderate and controlled, to improve muscle performance in the long term. Thus, if properly understood and carefully implemented, contraction induced injury can be used for the purpose of personalized training and recovery programs. In this paper, we review experimental, clinical, and theoretical works, attempting towards drawing a more quantitative description of contraction induced injury and related phenomena.

Mechanisms Of Residual Force Enhancement In Skeletal Muscle: Insights From Experiments And Mathematical Models

A skeletal muscle that is stretched while contracting will produce more force at steady state than if it is stretched passively and then stimulated to contract. This phenomenon is known as residual force enhancement and has been widely studied since its description more than sixty years ago. The idea that the mechanical properties of a muscle are governed not just by its present length but also by its length at earlier time points has far reaching implications since muscles stretch and shorten routinely in normal use. In this review, we present the experimental and theoretical advances that have been made toward understanding the mechanisms that underlie residual force enhancement. In the past ten years, experiments and models have focused on essentially three candidate mechanisms for residual force enhancement: (half-) sarcomere inhomogeneity, activity of so-called 'passive' mechanical elements in the sarcomere (titin), and the intrinsic properties of myosin crossbridges. Evidence, both computational and experimental, is accumulating for each of these mechanisms such that a final description of the phenomenon seems attainable in the near future. We conclude that computational models that incorporate more than one putative mechanism may ultimately facilitate reconciliation of the growing number of ideas and experimental data in this field.