- Department of Mechanical Engineering and Materials Science, Yale University, New Haven, CT, United States
The perturbation response of muscle is important for the versatile, stable and agile control capabilities of animals. Muscle resists being stretched by developing forces in the passive tissues and in the active crossbridges. This review focuses on the active perturbation response of the sarcomere. The active response exhibits typical stress relaxation, and thus approximated by a Maxwell material that has a spring and dashpot arranged in series. The ratio of damping to stiffness in this approximation defines the relaxation timescale for dissipating stresses that are developed in the crossbridges due to external perturbations. Current understanding of sarcomeres suggests that stiffness varies nearly linearly with neural excitation, but not much is known about damping. But if both stiffness and damping have the same functional (linear or not) dependence on neural excitation, then the stress relaxation timescale cannot be varied depending on the demands of the task. This implies an unavoidable and biologically unrealistic trade-off between how freely the crossbridges can yield and dissipate stresses when stretched (injury avoidance in agile motions) vs. how long they can maintain perturbation-induced stresses and behave like a solid material (stiffness maintenance for stability). We hypothesize that muscle circumvents this trade-off by varying damping in a nonlinear manner with neural excitation, unlike stiffness that varies linearly. Testing this hypothesis requires new experimental and mathematical characterization of muscle mechanics, and also identifies new design goals for robotic actuators.
1. Introduction
A muscle develops mechanical forces when neurally or electrically excited, and also when externally perturbed (Rack and Westbury, 1974; Kirsch et al., 1994; Lindstedt and Hoppeler, 2016). The perturbation response ΔFp of passive tissues and the active excitation-dependent perturbation response ΔFa add to the baseline force Fa generated due to neural excitations and to Fp due to passive tissues (Figure 1A). Passive refers to mechanical responses in the absence of neural stimulation while active responses require neural stimulation, and consume metabolic energy. Notably, the active resistance to stretch is one of the first mechanical responses of muscle when stimulated, even before it begins to develop tension (McMahon, 1984).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/344831/frobt-05-00069-HTML/image_m/frobt-05-00069-g001.jpg)
Figure 1. (A) A muscle provides, at every length, a passive tension Fp and an active tension Fa that varies with neural input. Furthermore, an external perturbation to its length elicits additional passive and active contributions, ΔFp and ΔFa, respectively. The tendons attached to the muscle work in series and act to transmit the total tension FT the muscle outputs, i.e., FT = Fp + ΔFp + Fa + ΔFa. (B) The active perturbation response arises from the contractile machinery of sarcomeres, and can be measured by applying length perturbations Δx to an isolated skinned muscle fiber. The perturbation response ΔFa is approximated by a Maxwell-type material that is comprised of a spring and dashpot in series. The spring has stiffness k(u) and dashpot has damping b(u), both of which vary with the neural input u. The ratio of the damping to stiffness determine a stress relaxation timescale τrelax. For a step length perturbation, the active perturbation response ΔFa(t) decays exponentially over the relaxation time.
Active and passive perturbation responses play an essential role in animal motor control because they are faster than any neural response, including the fastest of reflexes (Bizzi et al., 1982; Brown and Loeb, 2000; Dickinson et al., 2000; Hogan and Buerger, 2005; Holmes et al., 2006; Nishikawa et al., 2007; Biewener, 2016; Roberts, 2016). These responses have also been called preflexes or mechanical feedback (Brown and Loeb, 2000; Nishikawa et al., 2007). In robotics as well, the perturbation response of actuators are employed advantageously when appropriately tuned to the task and the environment's mechanical response (Hogan, 1984; Pratt and Williamson, 1995; Hogan and Buerger, 2005; Buerger and Hogan, 2007; Vanderborght et al., 2013). However, current actuator technologies do not yet match the ability of muscle to vary its active perturbation response by several fold, such as stiffness that may vary by over 50× in muscle (Hunter and Lafontaine, 1992; Madden, 2007; Anderson et al., 2012; Vanderborght et al., 2013; Hines et al., 2017). Therefore the mechanical capabilities of skeletal muscle continues to be sought-after by roboticists not only in terms of their force, work and power generation capabilities (Madden, 2007), but also in terms of their active perturbation response (Madden et al., 2004).
In this review, we examine the active perturbation response using the formalism of frequency-dependent dynamic modulus, i.e., the ratio of the force response ΔFa to externally imposed small sinusoidal length perturbations Δx (Figure 1B). Emphasizing qualitative over detailed quantitative comparisons with data, we use approaches from materials science (de Gennes, 1979; Barnes et al., 1989) to simplify the possibly complicated active response in terms of spring-like and dashpot-like elements (Figure 1B). In particular, because the active perturbation response of muscle does not have a single, simple rest length and decays with time, its simplest representation is a Maxwell material with a spring and dashpot in series (Maxwell, 1867; Barnes et al., 1989) whose stiffness k and damping b depend on the neural input u.
We use examples from motor control in section 2 to illustrate the role of active perturbation responses, and present an overview in section 3 of how a muscle's active perturbation response is characterized. We elaborate the Maxwell model in section 4, examine its neural modulation in section 5, and summarize the main conclusions of this review in section 6.
2. Active Perturbation Response in Motor Control
Active perturbation response of muscle depends strongly on the animal's intent and biomechanical context. At one extreme, muscle may behave like a stiff solid that undergoes no appreciable strain (e.g., co-contracted muscles to stiffen a joint). At another extreme it may offer little resistance and yield freely like a fluid (e.g., biceps in throwing). In general, muscles exhibit every behavior in between the extremes depending on the level of neural excitation (Cui et al., 2008; Sponberg and Full, 2008; Farahat and Herr, 2010; Hu et al., 2011; Sawicki et al., 2015). Consider three representative examples with different active perturbation responses of muscles and their analogs in robots.
First, when elastic energy storage and recovery are important, a muscle typically behaves as a stiff strut (Alexander and Bennet-Clark, 1977; Zajac, 1989; Biewener and Roberts, 2000; Roberts, 2016) so that most of the externally imposed strain is elastically stored in and recovered from the tendon (but see George et al., 2013). Use of the tendon's series elasticity to store and recover elastic energy in running is a well-known example (Cavagna et al., 1964). Robotic actuators employ a similar approach using series elastic elements for energy storage and recovery, in addition to protecting the actuator from shock loads (Raibert, 1986; Pratt and Williamson, 1995).
Second, when stabilizing a joint through muscle stiffening or regulating limb stiffness for controlling interactions with the surroundings, the stiffness of muscle is more continuously and appropriately varied to the dynamics and mechanics of the task being performed (Lacquaniti and Maioli, 1987; Lacquaniti et al., 1993; McIntyre et al., 1996; Burdet et al., 2001; Hogan and Buerger, 2005; Cui et al., 2008; Hu et al., 2011). In dynamical contexts, muscle's active perturbation response is generalized to an active impedance (Hogan, 1984; Hogan and Buerger, 2005) that may be approximately understood in terms of spring-like and dashpot-like responses to the external strain, strain-rate and the neural excitation. Varying limb impedance by modulating the impedance of the driving actuators has also been central to controlling interactions in robots (Hogan and Buerger, 2005; Vanderborght et al., 2013).
Third, some tasks involve quick transitions of muscles from a fluid-like to a solid-like response. For example, in throwing, the elbow's joint angular velocity exceeds 5000°/s before the elbow rapidly brakes and stiffens at the end to avoid injuries (Roach et al., 2013). Antagonistic muscles to this movement, such as the biceps brachii have to yield with little resistance like a weak dashpot would, in order to not prematurely decelerate the arm and to avoid injuries when a rapid stretch is imposed upon them (Lindstedt et al., 2001; LaStayo et al., 2003). At the end of the motion, the biceps brachii provide active braking to safely decelerate the arm without themselves suffering damage, thus transitioning into impedance control and ending with high stiffness. Such transitions in the material properties of actuators have not been studied or used in robots.
These examples illustrate that muscle's active perturbation response bridges the gap between two extremes. One extreme is the ability to maintain internal stresses that arise from an external perturbation for prolonged periods of time, lasting several minutes (Rancourt and Hogan, 2001; Loram et al., 2007), so that the muscle may function as a solid spring-like material for stability, elastic energy storage in tendons and so on. The other extreme is the ability to rapidly dissipate perturbation-induced internal stresses in mere tens of milliseconds (Roach et al., 2013) so that the muscle may function as a fluid dashpot-like material that enables agile and rapid motions without suffering damage.
3. Characterization of Muscle's Active Perturbation Response
Current understanding of muscle's mechanical behavior may be encapsulated by simple mathematical models such as Hill-type models of muscles (Zajac, 1989), Huxley-type models of sarcomeres (Walcott, 2014), and more intricate models of non-crossbridge elements such as the winding filament model of titin (Nishikawa et al., 2012; LeMoyne et al., 2014). Such simplified models are essential to elucidate underlying biological principles (Herzog, 2000), and to facilitate intensive computations such as applications of optimal control to study motor behaviors (Todorov, 2004). These mechanical models used to understand and characterize muscle have clear analogs in the centuries-old development of constitutive models of viscoelasticity in materials science (Barnes et al., 1989). An important lesson from materials science (Barnes et al., 1989) and from recent developments in soft (Wyss et al., 2007; Goldenfeld, 2018) and active or biological matter (Mizuno et al., 2007; Marchetti et al., 2013) has been that simple models, although often quantitatively inaccurate, guide experimental design and form the basis for fundamental mechanistic understanding.
Among viscoelastic constitutive models of materials, the Voigt model with a spring and dashpot in parallel, and the Maxwell model with a spring and dashpot in series are the two simplest approximations (Barnes et al., 1989). The Voigt model has been applied to characterize the viscoelastic properties of passive tissues, including that of muscle (Fung, 2013). But it is qualitatively wrong for the active perturbation response of muscle because it implies a single fixed length and the persistence of elastic stresses forever. We therefore consider the Maxwell model (Figure 1B) or a generalization called the standard linear model, which is a combination of a Maxwell body in parallel with a second weak spring representing parallel passive elasticity. These models are not literal representations of microscopic springs and dashpots, but their stiffness and damping are the respective linearized parameterizations of the reversible (elastic) and irreversible (viscous) components of the dynamic response to perturbations. In this sense, they are applicable to passive biological materials (Fung, 2013), as well as to the perturbation response of active biological materials (Deng et al., 2006; Mizuno et al., 2007).
We briefly summarize how muscle generates forces, with special focus on how it resists perturbations, and point the reader to more thorough examinations of the century-old topic of force generation in muscle (Herzog, 2000). Neural excitation in the form of a train of action potentials increases the concentration of freely available intracellular Ca2+ ions, and in turn increases the number of crossbridges formed between interdigitating filaments of actin and myosin (Figure 1B). The crossbridges form transient load-bearing links between actomyosin filaments. By executing a power stroke, each crossbridge contributes approximately 2pN to the active tension Fa. It also acts as a molecular spring with a stiffness of approximately 1pN/nm and contributes to the active perturbation response ΔFa (Finer et al., 1994). The collective behavior of crossbridges is such that they store elastic stresses when the sarcomere is externally perturbed, but slowly dissipate the stresses as crossbridges detach and re-attach elsewhere (Huxley, 1974). The presence of a relaxation timescale for stress dissipation suggests that the active perturbation response resembles a Maxwell viscoelastic material (Palmer et al., 2007).
3.1. Dynamic Modulus and Other Measures of Perturbation Response
The dynamic modulus is the ratio of the active perturbation response ΔFa to externally applied small stretches Δx, and depends on the excitation level (or intracellular [Ca2+]) and the rate (frequency) at which the stretch is applied. The behavior resembles a pure spring if the stresses ΔFa persist without any decay, and resembles a pure dashpot if the stresses decay (exponentially for a linear dashpot). Most materials exhibit both behaviors, depending on the timescales of the applied stretch and of observation. For example, a Maxwell body responds to a sudden stretch with a sharp pure spring-like transient, followed by a slower pure dashpot-like dissipation of stress (Figure 1B). These responses to small perturbations have also been characterized and identified as the short-range elastic component (SREC) (Rack and Westbury, 1974; Campbell, 2010), the complex modulus, or the sinusoidal response (Kawai and Brandt, 1980; Palmer et al., 2007).
The dynamic modulus arising from crossbridge dynamics is not the same as the slopes of the well-known force-length and force-velocity curves of the sarcomere (Zajac, 1989), which are sometimes misinterpreted as stiffness and damping. Although the slope of these respective curves have physical units of stiffness and damping, they do not represent the dynamic perturbation response of the sarcomere. For example, the force-length curve for a sarcomere near its natural length of ≈ 2.2μm has zero slope, which leads to a misinterpretation that an excited sarcomere has no stiffness at its natural length, a provable fallacy (Rack and Westbury, 1974). Seen from the perspective of dynamic systems, the sarcomere may be characterized as a material with a frequency-dependent perturbation response, and the force-length relationship quantifies the zero-frequency stiffness alone (Kawai and Brandt, 1980). It is tempting to set aside these complicated (and complex) perturbation responses of muscle as biological artifacts. However, as we argue here, the frequency-dependent perturbation response of muscle, and its modulation through neural excitation, are central to muscle's utility as a biological actuator.
The perturbation response of muscle has also been extensively used as a window into its microscopic functioning (McMahon, 1984). Among the many insights gleaned on sarcomere function, there remain major open questions such as the molecular origins of force enhancement (Rassier, 2012) and thixotropy (Campbell, 2010). Force enhancement (or depression) is the persistence of additional stresses for several seconds when an active sarcomere is externally stretched (contracted Edman et al., 1982). Thixotropy is a term borrowed from passive materials to refer to the history dependence of a cyclically stretched sarcomere so long as the time elapsed between consecutive stretches is sufficiently small (Campbell and Moss, 2002). These and similar studies have revealed that besides the actomyosin contractile machinery, muscle's force generation and perturbation response are additionally affected by other factors such as the surrounding passive tissues (Roberts, 2016), pressure and viscosity of the intra- and inter-cellular fluid (Baron et al., 2017; Sleboda and Roberts, 2017), and non-crossbridge but calcium-sensitive components like titin (Herzog, 2014; Nishikawa, 2016). Frequency-dependent dynamic modulus is related to the perturbation response protocols used in force enhancement and thixotropy studies, but generalizes to an examination of multiple timescales by adopting established protocols from rheology.
The frequency-dependent dynamic modulus of a sarcomere, myofibril, muscle fiber or a whole muscle is mostly attributable to the crossbridges only when the imposed stretches are sufficiently small so as to not forcibly detach bound crossbridges (McMahon, 1984). A natural concern therefore is whether the dynamic modulus is relevant for animal motor behaviors. Forcibly detaching crossbridges leads to microscopic damage that helps build muscle if the extent of damage is sufficiently small (LaStayo et al., 2003), but excessive damage is injurious. Thus when muscle is highly externally stretched, it is important for the sarcomere itself to remain stiff and enable the softer series elastic elements such as the tendon, the aponeuroses, and other passive tissues to accommodate a majority of the strain. In muscles with short tendons or when forcibly stretched by amounts that cannot be accommodated by tendons, the crossbridges must unbind and dissipate stresses fast enough so that a majority are not forcibly detached. Thus non-injurious perturbation response involves small elastic strains at individual crossbridges although the whole muscle or joint may experience large motions, and large crossbridge strains imply dashpot-like stress dissipation. We therefore argue that the small strain dynamic modulus of the muscle is relevant to non-pathological function, and is indeed known to play a role in a large variety of tasks (Bizzi et al., 1982; Lacquaniti and Maioli, 1987; Lin and Rymer, 1998, 2000, 2001; Loram et al., 2007; Cui et al., 2008; Rancourt and Hogan, 2009; Hu et al., 2011).
4. Active Perturbation Response of a Maxwell Material
A Maxwell-type response is modeled by a series arrangement of a spring with stiffness k(u) and dashpot with damping b(u), and captures the relaxation behavior of the active perturbation response ΔFa (Figure 1B). The dependence on neural input u represents the fact that the number of crossbridges and therefore the active response of a sarcomere varies with neural excitation (Herzog, 2000). Upon imposing a step perturbation to the length, the spring initially assumes all the externally applied strain and ΔFa = ℓk(u) where ℓ is the amplitude of the imposed step. The dashpot dissipates the stored elastic stress over a timescale τrelax = b(u)/k(u) and ΔFa decays exponentially in time. The response resembles an elastic solid against perturbations completed over a duration shorter than the relaxation timescale τrelax and a viscous fluid against longer perturbations.
The perturbation response of a Maxwell material is characterized by its frequency-dependent dynamic modulus K. The imposed length perturbation Δx(t) is distributed between the displacement Δxs of the spring and Δxd of the dashpot so that Δx = Δxs + Δxd. Using the constitutive laws for a spring and dashpot, the active perturbation response is given by
The normalized dynamic modulus K/k(u) is found using the Fourier transform of Equation (1), and depends on the frequency ω of the applied sinusoidal length perturbation and the stress relaxation timescale τrelax according to
where and are the Fourier transforms of Δx(t) and ΔFa(t), respectively. For fast perturbations with ωτrelax ≫ 1, the dynamic modulus K ≈ k(u), and the active perturbation response resembles a pure spring. For slower perturbations with ωτrelax ≪ 1, the dynamic modulus K ≈ τrelaxk(u)ω and decreases linearly with frequency. The prefactor τrelaxk(u) is equal to damping b(u), and the perturbation response resembles a pure dashpot.
The active perturbation response of skinned muscle fibers (Kawai and Brandt, 1980; Palmer et al., 2007) indeed resembles a Maxwell material's active response ΔFa plus a weak passive parallel elastic spring ΔFp (Figure 2A). The dynamic modulus measurement of skinned cardiac muscle fibers using sinusoidal perturbations of different frequencies is shown in Figure 2A, where activated fibers were held at the plateau region of the force-length curve and perturbed sinusoidally. Data from skeletal muscle (Kawai and Brandt, 1980; Miller et al., 2010; Palmer et al., 2013) show a similar response but those experiments did not perturb at high enough frequencies, because of which the dashpot-like response is evident but not the spring-like response.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/344831/frobt-05-00069-HTML/image_m/frobt-05-00069-g002.jpg)
Figure 2. (A) Sinusoidal analysis of skinned mouse cardiac muscle fibers shows an active perturbation response qualitatively similar to a Maxwell-type material for frequencies above 10 Hz. A parallel elastic response dominates at lower frequencies, and the active component is no longer evident. Data replotted from Palmer et al. (2007). (B) plot: The viscoelastic perturbation response for a Maxwell-type material where τrelax is the timescale of stress relaxation, K is the dynamic modulus of the muscle, and ω is the frequency of a sinusoidal perturbation. Against perturbations with periods shorter than τrelax, the response resembles an elastic solid, and it resembles a viscous fluid for slower perturbations. inset: If the relaxation time for the active perturbation response of a sarcomere depends on its neural input u, we hypothesize that the sarcomere may switch between solid and fluid behaves by varying u. For example, the active perturbation response may be fluid-like at low neural input (green) and solid-like at high neural input (red). (C) For the hypothesized dependence (Equation 6), τrelax is increasing with increasing u if d > 1, invariant with u if d = 1, and decreasing with increasing u if d < 1. Each case results in different functional consequences for a sarcomere and ultimately for muscles when coupled with the fact that a sarcomere's active tension Fa increases with u. A sarcomere can vary between a solid behavior with high active tension (mechanically stable against external perturbations, green shading) and a fluid behavior with low active tension (freely yields with minimal resistance, brown shading) only if d > 1.
Cardiac muscles have many physiological differences from skeletal muscles, but both are governed by similar biomechanical principles by virtue of relying on sarcomeres for active force production. Whether the Maxwell body is attributed to crossbridge dynamics, titin, or other active components of the sarcomere, the sinusoidal analysis experiments indicate that a single timescale, namely the stress relaxation timescale τrelax, separates the fluid and solid behaviors of the active perturbation response.
A non-dimensional version of the Maxwell model has no parameters and is obtained by normalizing the dynamic modulus K by the excitation-dependent high-frequency stiffness k(u), and normalizing the frequency ω by the reciprocal of the stress relaxation timescale τrelax (Figure 2B). Therefore if the relaxation timescale τrelax were made excitation-dependent so that the corner of the frequency response shifts right at lower excitation levels, it would also reduce the high frequency stiffness (red vs. green curves in Figure 2B inset). Thus, at the same frequency of perturbation, the active response may transition from a stiff solid to a weakly visous fluid if the excitation is sufficiently reduced.
5. Varying the Relaxation Timescale
The stress relaxation timescale τrelax of the active perturbation response is the ratio of the sarcomere's damping to its stiffness, and may vary with neural excitation depending on how stiffness and damping vary (Figure 2B). While it is known that the stiffness increases linearly with neural input (Rack and Westbury, 1974; Kirsch et al., 1994), it remains unknown how damping may vary. To examine the neural modulation of τrelax, we consider a linear form for the stiffness k(u), but allow a general power-law for how damping b(u) varies with neural input according to
Thus the stress relaxation timescale is given by the ratio of damping to stiffness as
Neural excitation is parameterized by the normalized variable 0 ≤ u ≤ 1, and τ0 = b0/k0 is the relaxation timescale at maximal excitation.
Recall that evidence from motor control (section 2) indicates that the stress relaxation timescale increases with increasing neural drive. Highly activated muscle resembles an elastic solid whose relaxation timescale is much greater than the duration of the experiment or behavior (green shading in Figure 2C), and minimally activated muscle resembles a weakly viscous fluid whose relaxation timescale is much smaller than the fastest perturbation in the experiment or behavior (brown shading in Figure 2C).
Only for d > 1 is the qualitatively correct behavior observed for τrelax as a function of u (blue solid line, left panel of Figure 2C), and if both stiffness and damping scale equally with excitation (d = 1) the timescale becomes an invariant quantity (blue solid line, middle panel Figure 2C). For d < 1 the modulation of τrelax with increasing u makes the active perturbation response weak spring-like for weak neural input and strong dashpot-like for strong neural input (blue solid line, right panel Figure 2C). This case would correspond to large crossbridge strains and frequent forced detachments of crossbridges at both extremes of muscle function, solid or fluid.
In all three cases, the mean active tension Fa would increase with the excitation u (red dashed line, Figure 2C), and therefore the dependence of the relaxation timescale τrelax on the neural input u is equivalent to being dependent on the mean active tension Fa. The functionally desirable behavior is therefore one of a force-dependent slowing of stress relaxation in the active perturbation response.
Hill-type muscle-tendon models that incorporate a series elastic tendon, and treat the force-length and force-velocity curves as the stiffness and damping of the active perturbation response, exhibit a similar frequency response to the Maxwell model (Figures 16, 18, and 19 in Zajac, 1989). The high-pass filtering characteristics of those models arise because the dynamic modulus of the active element is effectively zero and it behaves as a pure damper in response to perturbations. Thus the active response plus the passive elastic spring of the tendon is simply a Maxwell model. Furthermore, when the Hill-type models operate on the ascending limb of the force-length curve, the zero-frequency stiffness of the active response is non-zero and the stress relaxation timescale (corner in the frequency response) becomes excitation-dependent. However, the functional dependence of the stress relaxation timescale on the excitation in the Hill-type models corresponds to the unrealistic d < 1 scenario, where increasing excitation makes the muscle-tendon more fluid-like. This increased fluid-like response may however be relevant in dynamic tasks such as locomotion, where submaximally activated muscle operates primarily on the ascending limb of the force-length curve (Holt and Azizi, 2016).
The active perturbation response of muscle is also modulated through neural feedback circuits, the stretch reflex being the fastest of them (McMahon, 1984). Even those fast reflexes take over 50ms in humans and therefore it can impact the perturbation response only at frequencies 20Hz or below. At these lower frequencies, reflexes have the capability of altering the frequency response, including the effective damping (Lin and Rymer, 1998, 2001). The analysis presented in this paper does not incorporate the effect of neural reflexes, and neither do the single fiber experiments using sinusoidal analysis (Kawai and Brandt, 1980; Palmer et al., 2007, 2013; Miller et al., 2010).
6. Conclusion
We hypothesize that an excitation-dependent increase in the damping associated with the active perturbation response should outstrip the stiffness increase. This is necessary to explain observed muscle behavior and accomplish tasks that require the muscle to mimic phase transitions between a fluid-like and solid-like active response. The microscopic origins of how the stress relaxation timescale may vary is an important avenue for future research. Whether it shares mechanisms with force-enhancement, thixotropy or even the well-known load-dependent changes in the detachment rate of actin-bound myosin remains to be discovered. Our proposal for studying the active perturbation response of sarcomeres has analogies to the applications of rheology in the fields of active and biological matter (Deng et al., 2006; Mizuno et al., 2007; Wyss et al., 2007; Marchetti et al., 2013).
Based on the known uses of muscle in control, and current theories of sarcomere function, we have argued for treating the active perturbation response of the sarcomere as a Maxwell-type material with the elastic and viscous elements in series. We note that our conclusion differs from some authors who express a preference for a Voigt element, although they conclude that both models could fit their experimental data (Ford et al., 1977). A significant reason for the differing viewpoints is because we explicitly separate the active perturbation response from all other force production components, namely the passive perturbation response, and the active and passive mean force generation, while (Ford et al., 1977) have lumped them together.
A renewed examination of the active perturbation response of muscle and its control through neural excitation may provide new design targets for engineered actuators for use in agile animal-like robots. However, significant technological challenges remain in developing actuators that can modulate stiffness and damping. The prevalent use of high-bandwidth feedback control to mimic viscoelasticity creates fragile devices that are sensitive to sensor malfunction, time-delays and noise (Hogan and Buerger, 2005). On the other hand, novel soft actuators (Majidi, 2014; Rus and Tolley, 2015; Hines et al., 2017) such as dielectric (Anderson et al., 2012) or nematic (Hines et al., 2017) elastomers, twisting cable muscles (Haines et al., 2016) or pneumatic actuators (Wehner et al., 2014; Peele et al., 2015) resemble a Voigt material rather than a Maxwell material. Therefore, they suffer the severe limitations imposed by a fixed neutral length and a strong parallel elastic component. As an alternative, variable stiffness and damping capabilities may be achieved through mechanical design of a transmission element (Vanderborght et al., 2013). These designs also have a key limitation, namely that force, stiffness and damping are typically varied through independent control inputs and thereby lead to high control complexity. We do not imply that matching the success of muscles as actuators needs mimicry of its microscopic structure. Rather, developing novel actuators with a Maxwell-type perturbation response that can undergo large changes in the stress relaxation timescale may prove fruitful in mimicking the beneficial principles underlying muscle's versatile use by animals.
Author Contributions
KN and NS contributed equally. MV conceived the paper. KN synthesized the literature on sarcomere mechanics. NS synthesized the motor control literature. KN, NS, and MV wrote and edited the paper together.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
Bertrand Tanner and Samuel Walcott for comments and discussions. Funding support from the Raymond and Beverly Sackler Institute for Biological, Physical and Engineering Sciences at Yale.
References
Alexander, R. M., and Bennet-Clark, H. C. (1977). Storage of elastic strain energy in muscle and other tissues. Nature 265:114.
Anderson, I. A., Gisby, T. A., McKay, T. G., OBrien, B. M., and Calius, E. P. (2012). Multi-functional dielectric elastomer artificial muscles for soft and smart machines. J. Appl. Phys. 112:041101. doi: 10.1063/1.4740023
Barnes, H. A., Hutton, J. F., and Walters, K. (1989). An introduction to rheology, Vol. 3. Amsterdam: Elsevier.
Baron, M., Hosoi, A., Williams, C., and Daniel, T. (2017). “Flow in the lattice of myofilaments,” in Integrative and Comparative Biology, Vol. 57. (Oxford), E200.
Biewener, A. A. (2016). Locomotion as an emergent property of muscle contractile dynamics. J. Exp. Biol. 219, 285–294. doi: 10.1242/jeb.123935
Biewener, A. A., and Roberts, T. J. (2000). Muscle and tendon contributions to force, work, and elastic energy savings: a comparative perspective. Exerc. Sport Sci. Rev. 28, 99–107. Available online at: https://journals.lww.com/acsm-essr/Citation/2000/28030/Muscle_and_Tendon_Contributions_to_Force,_Work,.2.aspx
Bizzi, E., Chapple, W., and Hogan, N. (1982). Mechanical properties of muscles: implications for motor control. Trends Neurosci. 5, 395–398. doi: 10.1016/0166-2236(82)90221-1
Brown, I. E., and Loeb, G. E. (2000). “A reductionist approach to creating and using neuromusculoskeletal models,” in Biomechanics and Neural ontrol of Posture and Movement (New York, NY: Springer), 148–163.
Buerger, S. P., and Hogan, N. (2007). Complementary stability and loop shaping for improved human–robot interaction. IEEE Transactions on Robotics 23, 232–244. doi: 10.1109/TRO.2007.892229
Burdet, E., Osu, R., Franklin, D. W., Milner, T. E., and Kawato, M. (2001). The central nervous system stabilizes unstable dynamics by learning optimal impedance. Nature 414, 446–449. doi: 10.1038/35106566
Campbell, K. S. (2010). “Short-range mechanical properties of skeletal and cardiac muscles,” in Muscle Biophysics, ed D. Rassier (New York, NY: Springer), 223–246.
Campbell, K. S., and Moss, R. L. (2002). History-dependent mechanical properties of permeabilized rat soleus muscle fibers. Biophys. J. 82, 929–943. doi: 10.1016/S0006-3495(02)75454-4
Cavagna, G., Saibene, F., and Margaria, R. (1964). Mechanical work in running. J. Appl. Physiol. 19, 249–256. doi: 10.1152/jappl.1964.19.2.249
Cui, L., Perreault, E. J., Maas, H., and Sandercock, T. G. (2008). Modeling short-range stiffness of feline lower hindlimb muscles. J. Biomechanics 41, 1945–1952. doi: 10.1016/j.jbiomech.2008.03.024
Deng, L., Trepat, X., Butler, J. P., Millet, E., Morgan, K. G., Weitz, D. A., et al. (2006). Fast and slow dynamics of the cytoskeleton. Nat. Materials 5:636–640. doi: 10.1038/nmat1685
Dickinson, M. H., Farley, C. T., Full, R. J., Koehl, M. A., Kram, R., and Lehman, S. (2000). How animals move: an integrative view. Science 288, 100–106. doi: 10.1126/science.288.5463.100
Edman, K. A., Elzinga, G., and Noble, M. (1982). Residual force enhancement after stretch of contracting frog single muscle fibers. J. Gen. Physiol. 80, 769–784. doi: 10.1085/jgp.80.5.769
Farahat, W. A., and Herr, H. M. (2010). Optimal workloop energetics of muscle-actuated systems: an impedance matching view. PLoS Comput. Biol. 6:e1000795. doi: 10.1371/journal.pcbi.1000795
Finer, J. T., Simmons, R. M., and Spudich, J. A. (1994). Single myosin molecule mechanics: piconewton forces and nanometre steps. Nature 368, 113–119. doi: 10.1038/368113a0
Ford, L. E., Huxley, A. F., and Simmons, R. M. (1977). Tension responses to sudden length change in stimulated frog muscle fibres near slack length. J. Physiol. 269, 441–515. doi: 10.1113/jphysiol.1977.sp011911
Fung, Y.-C. (2013). Biomechanics: Mechanical Properties of Living Tissues. New York, NY: Springer Science and Business Media.
George, N. T., Irving, T. C., Williams, C. D., and Daniel, T. L. (2013). The cross-bridge spring: can cool muscles store elastic energy? Science 340, 1217–1220. doi: 10.1126/science.1229573
Goldenfeld, N. (2018). Lectures on Phase Transitions and the Renormalization Group. Boca Raton, FL: CRC Press.
Haines, C. S., Li, N., Spinks, G. M., Aliev, A. E., Di, J., and Baughman, R. H. (2016). New twist on artificial muscles. Proc. Natl. Acad. Sci. U.S.A. 113, 11709–11716. doi: 10.1073/pnas.1605273113
Herzog, W. (2000). Skeletal Muscle Mechanics: From Mechanisms to Function. Hoboken, NJ: John Wiley and Sons.
Herzog, W. (2014). The role of titin in eccentric muscle contraction. J. Exp. Biol. 217, 2825–2833. doi: 10.1242/jeb.099127
Hines, L., Petersen, K., Lum, G. Z., and Sitti, M. (2017). Soft actuators for small-scale robotics. Adv. Mat. 29, 1–43. doi: 10.1002/adma.201603483
Hogan, N. (1984). Adaptive control of mechanical impedance by coactivation of antagonist muscles. IEEE Trans. Automatic Control 29, 681–690. doi: 10.1109/TAC.1984.1103644
Hogan, N., and Buerger, S. P. (2005). Impedance and interaction control. Robot. Autom. Handb. 1, 19-1–19-24.
Holmes, P., Full, R. J., Koditschek, D., and Guckenheimer, J. (2006). The dynamics of legged locomotion: models, analyses, and challenges. SIAM Rev. 48, 207–304. doi: 10.1137/S0036144504445133
Holt, N. C., and Azizi, E. (2016). The effect of activation level on muscle function during locomotion: are optimal lengths and velocities always used? Proc. R. Soc. B 283:20152832. doi: 10.1098/rspb.2015.2832
Hu, X., Murray, W. M., and Perreault, E. J. (2011). Muscle short-range stiffness can be used to estimate the endpoint stiffness of the human arm. J. Neurophysiol. 105, 1633–1641. doi: 10.1152/jn.00537.2010
Hunter, I. W., and Lafontaine, S. (1992). “A comparison of muscle with artificial actuators,” in Solid-State Sensor and Actuator Workshop, 1992. 5th Technical Digest, IEEE (Piscataway, NJ: IEEE), 178–185.
Kawai, M., and Brandt, P. W. (1980). Sinusoidal analysis: a high resolution method for correlating biochemical reactions with physiological processes in activated skeletal muscles of rabbit, frog and crayfish. J. Muscle Res. Cell Motil. 1, 279–303. doi: 10.1007/BF00711932
Kirsch, R. F., Boskov, D., and Rymer, W. Z. (1994). Muscle stiffness during transient and continuous movements of cat muscle: perturbation characteristics and physiological relevance. IEEE Trans. Biomed. Eng. 41, 758–770. doi: 10.1109/10.310091
Lacquaniti, F., Carrozzo, M., and Borghese, N. A. (1993). Time-varying mechanical behavior of multijointed arm in man. J. Neurophysiol. 69, 1443–1464. doi: 10.1152/jn.1993.69.5.1443
Lacquaniti, F., and Maioli, C. (1987). Anticipatory and reflex coactivation of antagonist muscles in catching. Brain Res. 406, 373–378. doi: 10.1016/0006-8993(87)90810-9
LaStayo, P. C., Woolf, J. M., Lewek, M. D., Snyder-Mackler, L., Reich, T., and Lindstedt, S. L. (2003). Eccentric muscle contractions: their contribution to injury, prevention, rehabilitation, and sport. J. Orthopaed. Sports Phys. Ther. 33, 557–571. doi: 10.2519/jospt.2003.33.10.557
LeMoyne, R., Petak, J., Tester, J., and Nishikawa, K. (2014). “Simulation of a computational winding filament model with an exponential spring to represent titin,” in 2014 36th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC) (Piscataway, NJ: IEEE), 836–839.
Lin, D. C., and Rymer, W. Z. (1998). Damping in reflexively active and areflexive lengthening muscle evaluated with inertial loads. J. Neurophysiol. 80, 3369–3372. doi: 10.1152/jn.1998.80.6.3369
Lin, D. C., and Rymer, W. Z. (2000). Damping actions of the neuromuscular system with inertial loads: soleus muscle of the decerebrate cat. J. Neurophysiol. 83, 652–658. doi: 10.1152/jn.2000.83.2.652
Lin, D. C., and Rymer, W. Z. (2001). Damping actions of the neuromuscular system with inertial loads: human flexor pollicis longus muscle. J. Neurophysiol. 85, 1059–1066. doi: 10.1152/jn.2001.85.3.1059
Lindstedt, S. L., and Hoppeler, H. H. (2016). Expanding knowledge of contracting muscle. J. Exp. Biol. 219:134. doi: 10.1242/jeb.135921
Lindstedt, S. L., LaStayo, P. C., and Reich, T. E. (2001). When active muscles lengthen: properties and consequences of eccentric contractions. News Physiol. Sci. 16, 256–261. doi: 10.1152/physiologyonline.2001.16.6.256
Loram, I. D., Maganaris, C. N., and Lakie, M. (2007). The passive, human calf muscles in relation to standing: the short range stiffness lies in the contractile component. J. Physiol. 584, 677–692. doi: 10.1113/jphysiol.2007.140053
Madden, J. D. (2007). Mobile robots: motor challenges and materials solutions. Science 318, 1094–1097. doi: 10.1126/science.1146351
Madden, J. D., Vandesteeg, N. A., Anquetil, P. A., Madden, P. G., Takshi, A., Pytel, R. Z., et al. (2004). Artificial muscle technology: physical principles and naval prospects. IEEE J. Oceanic Eng. 29, 706–728. doi: 10.1109/JOE.2004.833135
Majidi, C. (2014). Soft robotics: a perspective on current trends and prospects for the future. Soft Robot. 1, 5–11. doi: 10.1089/soro.2013.0001
Marchetti, M. C., Joanny, J.-F., Ramaswamy, S., Liverpool, T. B., Prost, J., Rao, M., et al. (2013). Hydrodynamics of soft active matter. Rev. Mod. Phys. 85:1143. doi: 10.1103/RevModPhys.85.1143
Maxwell, J. C. (1867). On the dynamical theory of gases. Philos. Trans. R. Soc. Lond. 157, 49–88. doi: 10.1098/rstl.1867.0004
McIntyre, J., Mussa-Ivaldi, F. A., and Bizzi, E. (1996). The control of stable postures in the multijoint arm. Exp. Brain Res. 110, 248–264. doi: 10.1007/BF00228556
McMahon, T. A. (1984). Muscles, Reflexes, and Locomotion. Princeton, NJ: Princeton University Press.
Miller, M. S., VanBuren, P., LeWinter, M. M., Braddock, J. M., Ades, P. A., Maughan, D. W., et al. (2010). Chronic heart failure decreases cross-bridge kinetics in single skeletal muscle fibres from humans. J. Physiol. 588, 4039–4053. doi: 10.1113/jphysiol.2010.191957
Mizuno, D., Tardin, C., Schmidt, C. F., and MacKintosh, F. C. (2007). Nonequilibrium mechanics of active cytoskeletal networks. Science 315, 370–373. doi: 10.1126/science.1134404
Nishikawa, K. (2016). Eccentric contraction: unraveling mechanisms of force enhancement and energy conservation. J. Exp. Biol. 219, 189–196. doi: 10.1242/jeb.124057
Nishikawa, K., Biewener, A. A., Aerts, P., Ahn, A. N., Chiel, H. J., Daley, M. A., et al. (2007). Neuromechanics: an integrative approach for understanding motor control. Integr. Comp. Biol. 47, 16–54. doi: 10.1093/icb/icm024
Nishikawa, K. C., Monroy, J. A., Uyeno, T. E., Yeo, S. H., Pai, D. K., and Lindstedt, S. L. (2012). Is titin a winding filament? a new twist on muscle contraction. Proc. R. Soc. Lond. Biol. Sci. 279, 981–990. doi: 10.1098/rspb.2011.1304
Palmer, B. M., Suzuki, T., Wang, Y., Barnes, W. D., Miller, M. S., and Maughan, D. W. (2007). Two-state model of acto-myosin attachment-detachment predicts c-process of sinusoidal analysis. Biophys. J. 93, 760–769. doi: 10.1529/biophysj.106.101626
Palmer, B. M., Tanner, B. C., Toth, M. J., and Miller, M. S. (2013). An inverse power-law distribution of molecular bond lifetimes predicts fractional derivative viscoelasticity in biological tissue. Biophys. J. 104, 2540–2552. doi: 10.1016/j.bpj.2013.04.045
Peele, B. N., Wallin, T. J., Zhao, H., and Shepherd, R. F. (2015). 3d printing antagonistic systems of artificial muscle using projection stereolithography. Bioinspir. Biomimet. 10:055003. doi: 10.1088/1748-3190/10/5/055003
Pratt, G. A., and Williamson, M. M. (1995). “Series elastic actuators,” in Proceedings. 1995 IEEE/RSJ International Conference on Intelligent Robots and Systems 95.'Human Robot Interaction and Cooperative Robots', Vol. 1 (Piscataway, NJ: IEEE), 399–406.
Rack, P. M., and Westbury, D. R. (1974). The short range stiffness of active mammalian muscle and its effect on mechanical properties. J. Physiol. 240:331.
Rancourt, D., and Hogan, N. (2001). Stability in force-production tasks. J. Motor Behav. 33, 193–204. doi: 10.1080/00222890109603150
Rancourt, D., and Hogan, N. (2009). “The biomechanics of force production,” in Progress in Motor Control (New York, NY: Springer), 645–661.
Rassier, D. E. (2012). Residual force enhancement in skeletal muscles: one sarcomere after the other. J. Muscle Res. Cell Motil. 33, 155–165. doi: 10.1007/s10974-012-9308-7
Roach, N. T., Venkadesan, M., Rainbow, M. J., and Lieberman, D. E. (2013). Elastic energy storage in the shoulder and the evolution of high-speed throwing in homo. Nature 498, 483–486. doi: 10.1038/nature12267
Roberts, T. J. (2016). Contribution of elastic tissues to the mechanics and energetics of muscle function during movement. J. Exp. Biol. 219, 266–275. doi: 10.1242/jeb.124446
Rus, D., and Tolley, M. T. (2015). Design, fabrication and control of soft robots. Nature 521, 467–475. doi: 10.1038/nature14543
Sawicki, G. S., Robertson, B. D., Azizi, E., and Roberts, T. J. (2015). Timing matters: tuning the mechanics of a muscle–tendon unit by adjusting stimulation phase during cyclic contractions. J. Exp. Biol. 218, 3150–3159. doi: 10.1242/jeb.121673
Sleboda, D. A., and Roberts, T. J. (2017). Incompressible fluid plays a mechanical role in the development of passive muscle tension. Biol. Lett. 13:20160630. doi: 10.1098/rsbl.2016.0630
Sponberg, S., and Full, R. J. (2008). Neuromechanical response of musculo-skeletal structures in cockroaches during rapid running on rough terrain. J. Exp. Biol. 211, 433–446. doi: 10.1242/jeb.012385
Todorov, E. (2004). Optimality principles in sensorimotor control. Nat. Neurosci. 7:907–915. doi: 10.1038/nn1309
Vanderborght, B., Albu-Schäffer, A., Bicchi, A., Burdet, E., Caldwell, D. G., Carloni, R., et al. (2013). Variable impedance actuators: A review. Robot. Auton. Sys. 61, 1601–1614. doi: 10.1016/j.robot.2013.06.009
Walcott, S. (2014). Muscle activation described with a differential equation model for large ensembles of locally coupled molecular motors. Phys. Rev. E 90:042717. doi: 10.1103/PhysRevE.90.042717
Wehner, M., Tolley, M. T., Mengüç, Y., Park, Y.-L., Mozeika, A., Ding, Y., et al. (2014). Pneumatic energy sources for autonomous and wearable soft robotics. Soft Robot. 1, 263–274. doi: 10.1089/soro.2014.0018
Wyss, H. M., Miyazaki, K., Mattsson, J., Hu, Z., Reichman, D. R., and Weitz, D. A. (2007). Strain-rate frequency superposition: a rheological probe of structural relaxation in soft materials. Phys. Rev. Lett. 98:238303. doi: 10.1103/PhysRevLett.98.238303
Keywords: muscle viscoelasticity, sarcomere mechanics, active perturbation response, variable impedance, stress relaxation timescale, dynamic modulus, sinusoidal response
Citation: Nguyen KD, Sharma N and Venkadesan M (2018) Active Viscoelasticity of Sarcomeres. Front. Robot. AI 5:69. doi: 10.3389/frobt.2018.00069
Received: 23 December 2017; Accepted: 25 May 2018;
Published: 14 June 2018.
Edited by:
Monica A. Daley, Royal Veterinary College, United KingdomReviewed by:
Eric D. Tytell, Tufts University, United StatesAllison Arnold, Harvard University, United States
Gregory S. Sawicki, University of North Carolina at Chapel Hill, United States
Copyright © 2018 Nguyen, Sharma and Venkadesan. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Madhusudhan Venkadesan, mv@classicalmechanic.net
†These authors have contributed equally to this work.