Search

Deficits in the ankle plantarflexor muscles, such as weakness and contracture, occur commonly in conditions such as cerebral palsy, stroke, muscular dystrophy, Charcot-Marie-Tooth disease, and sarcopenia. While these deficits likely contribute to observed gait pathologies, determining cause-effect relationships is difficult due to the often co-occurring biomechanical and neural deficits. To elucidate the effects of weakness and contracture, we systematically introduced isolated deficits into a musculoskeletal model and generated simulations of walking to predict gait adaptations due to these deficits. We trained a planar model containing 9 degrees of freedom and 18 musculotendon actuators to walk using a custom optimization framework through which we imposed simple objectives, such as minimizing cost of transport while avoiding falling and injury, and maintaining head stability. We first generated gaits at prescribed speeds between 0.50 m/s and 2.00 m/s that reproduced experimentally observed kinematic, kinetic, and metabolic trends for walking. We then generated a gait at self-selected walking speed; quantitative comparisons between our simulation and experimental data for joint angles, joint moments, and ground reaction forces showed root-mean-squared errors of less than 1.6 standard deviations and normalized cross-correlations above 0.8 except for knee joint moment trajectories. Finally, we applied mild, moderate, and severe levels of muscle weakness or contracture to either the soleus (SOL) or gastrocnemius (GAS) or both of these major plantarflexors (PF) and retrained the model to walk at a self-selected speed. The model was robust to all deficits, finding a stable gait in all cases. Severe PF weakness caused the model to adopt a slower, "heel-walking" gait. Severe contracture of only SOL or both PF yielded similar results: the model adopted a "toe-walking" gait with excessive hip and knee flexion during stance. These results highlight how plantarflexor weakness and contracture may contribute to observed gait patterns. [ABSTRACT FROM AUTHOR]

This study presents a modelling framework in which information on muscle fiber direction and orientation during contraction is derived from diffusion tensor imaging (DTI) and incorporated in a computational model of the surface electromyographic (EMG) signal. The proposed model makes use of the principle of reciprocity to simultaneously calculate the electric potentials produced at the recording electrode by charges distributed along an arbitrary number of muscle fibers within the muscle, allowing for a computationally efficient evaluation of extracellular motor unit action potentials. The approach is applied to the complex architecture of the first dorsal interosseous (FDI) muscle of the hand to simulate EMG during index finger flexion and abduction. Using diffusion tensor imaging methods, the results show how muscle fiber orientation and curvature in this intrinsic hand muscle change during flexion and abduction. Incorporation of anatomically accurate muscle architecture and other hand tissue morphologies enable the model to capture variations in extracellular action potential waveform shape across the motor unit population and to predict experimentally observed differences in EMG signal features when switching from index finger abduction to flexion. The simulation results illustrate how structural and electrical properties of the tissues comprising the volume conductor, in combination with fiber direction and curvature, shape the detected action potentials. Using the model, the relative contribution of motor units of different sizes located throughout the muscle under both conditions is examined, yielding a prediction of the detection profile of the surface EMG electrode array over the muscle cross-section. [ABSTRACT FROM AUTHOR]

Cytoplasmic flows are an ubiquitous feature of biological systems, in particular in large cells, such as oocytes and eggs in early animal development. Here we show that cytoplasmic flows in starfish oocytes, which can be imaged well with transmission light microscopy, are fully determined by the cortical dynamics during surface contraction waves. We first show that the dynamics of the oocyte surface is highly symmetric around the animal-vegetal axis. We then mathematically solve the Stokes equation for flows inside a deforming sphere using the measured surface displacements as boundary conditions. Our theoretical predictions agree very well with the intracellular flows quantified by particle image velocimetry, proving that during this stage the starfish cytoplasm behaves as a simple Newtonian fluid on the micrometer scale. We calculate the pressure field inside the oocyte and find that its gradient is too small as to explain polar body extrusion, in contrast to earlier suggestions. Myosin II inhibition by blebbistatin confirms this conclusion, because it diminishes cell shape changes and hydrodynamic flow, but does not abolish polar body formation. [ABSTRACT FROM AUTHOR]

Hill-type muscle models are widely used within the field of biomechanics to predict and understand muscle behaviour, and are often essential where muscle forces cannot be directly measured. However, these models have limited accuracy, particularly during cyclic contractions at the submaximal levels of activation that typically occur during locomotion. To address this issue, recent studies have incorporated effects into Hill-type models that are oftentimes neglected, such as size-dependent, history-dependent, and activation-dependent effects. However, the contribution of these effects on muscle performance has yet to be evaluated under common contractile conditions that reflect the range of activations, strains, and strain rates that occur in vivo. The purpose of this study was to develop a modelling framework to evaluate modifications to Hill-type muscle models when they contract in cyclic loops that are typical of locomotor muscle function. Here we present a modelling framework composed of a damped harmonic oscillator in series with a Hill-type muscle actuator that consists of a contractile element and parallel elastic element. The intrinsic force-length and force-velocity properties are described using Bezier curves where we present a system to relate physiological parameters to the control points for these curves. The muscle-oscillator system can be geometrically scaled while preserving dynamic and kinematic similarity to investigate the muscle size effects while controlling for the dynamics of the harmonic oscillator. The model is driven by time-varying muscle activations that cause the muscle to cyclically contract and drive the dynamics of the harmonic oscillator. Thus, this framework provides a platform to test current and future Hill-type model formulations and explore factors affecting muscle performance in muscles of different sizes under a range of cyclic contractile conditions. [ABSTRACT FROM AUTHOR]

Involuntary force variability below 15 Hz arises from, and is influenced by, many factors including descending neural drive, proprioceptive feedback, and mechanical properties of muscles and tendons. However, their potential interactions that give rise to the well-structured spectrum of involuntary force variability are not well understood due to a lack of experimental techniques. Here, we investigated the generation, modulation, and interactions among different sources of force variability using a physiologically-grounded closed-loop simulation of an afferented muscle model. The closed-loop simulation included a musculotendon model, muscle spindle, Golgi tendon organ (GTO), and a tracking controller which enabled target-guided force tracking. We demonstrate that closed-loop control of an afferented musculotendon suffices to replicate and explain surprisingly many cardinal features of involuntary force variability. Specifically, we present 1) a potential origin of low-frequency force variability associated with co-modulation of motor unit firing rates (i.e.,‘common drive’), 2) an in-depth characterization of how proprioceptive feedback pathways suffice to generate 5-12 Hz physiological tremor, and 3) evidence that modulation of those feedback pathways (i.e., presynaptic inhibition of Ia and Ib afferents, and spindle sensitivity via fusimotor drive) influence the full spectrum of force variability. These results highlight the previously underestimated importance of closed-loop neuromechanical interactions in explaining involuntary force variability during voluntary ‘isometric’ force control. Furthermore, these results provide the basis for a unifying theory that relates spinal circuitry to various manifestations of altered involuntary force variability in fatigue, aging and neurological disease. [ABSTRACT FROM AUTHOR]

Current models of sensorimotor control posit that motor commands are generated by combining multiple modules which may consist of internal models, motor primitives or motor synergies. The mechanisms which select modules based on task requirements and modify their output during learning are therefore critical to our understanding of sensorimotor control. Here we develop a novel modular architecture for multi-dimensional tasks in which a set of fixed primitives are each able to compensate for errors in a single direction in the task space. The contribution of the primitives to the motor output is determined by both top-down contextual information and bottom-up error information. We implement this model for a task in which subjects learn to manipulate a dynamic object whose orientation can vary. In the model, visual information regarding the context (the orientation of the object) allows the appropriate primitives to be engaged. This top-down module selection is implemented by a Gaussian function tuned for the visual orientation of the object. Second, each module's contribution adapts across trials in proportion to its ability to decrease the current kinematic error. Specifically, adaptation is implemented by cosine tuning of primitives to the current direction of the error, which we show to be theoretically optimal for reducing error. This error-tuned model makes two novel predictions. First, interference should occur between alternating dynamics only when the kinematic errors associated with each oppose one another. In contrast, dynamics which lead to orthogonal errors should not interfere. Second, kinematic errors alone should be sufficient to engage the appropriate modules, even in the absence of contextual information normally provided by vision. We confirm both these predictions experimentally and show that the model can also account for data from previous experiments. Our results suggest that two interacting processes account for module selection during sensorimotor control and learning. [ABSTRACT FROM AUTHOR]

Muscle spindle proprioceptive receptors play a primary role in encoding the effects of external mechanical perturbations to the body. During externally-imposed stretches of passive, i.e. electrically-quiescent, muscles, the instantaneous firing rates (IFRs) of muscle spindles are associated with characteristics of stretch such as length and velocity. However, even in passive muscle, there are history-dependent transients of muscle spindle firing that are not uniquely related to muscle length and velocity, nor reproduced by current muscle spindle models. These include acceleration-dependent initial bursts, increased dynamic response to stretch velocity if a muscle has been isometric, and rate relaxation, i.e., a decrease in tonic IFR when a muscle is held at a constant length after being stretched. We collected muscle spindle spike trains across a variety of muscle stretch kinematic conditions, including systematic changes in peak length, velocity, and acceleration. We demonstrate that muscle spindle primary afferents in passive muscle fire in direct relationship to muscle force-related variables, rather than length-related variables. Linear combinations of whole muscle-tendon force and the first time derivative of force (dF/dt) predict the entire time course of transient IFRs in muscle spindle Ia afferents during stretch (i.e., lengthening) of passive muscle, including the initial burst, the dynamic response to lengthening, and rate relaxation following lengthening. Similar to acceleration scaling found previously in postural responses to perturbations, initial burst amplitude scaled equally well to initial stretch acceleration or dF/dt, though later transients were only described by dF/dt. The transient increase in dF/dt at the onset of lengthening reflects muscle short-range stiffness due to cross-bridge dynamics. Our work demonstrates a critical role of muscle cross-bridge dynamics in history-dependent muscle spindle IFRs in passive muscle lengthening conditions relevant to the detection and sensorimotor response to mechanical perturbations to the body, and to previously-described history-dependence in perception of limb position. [ABSTRACT FROM AUTHOR]

Protein design remains an important problem in computational structural biology. Current computational protein design methods largely use physics-based methods, which make use of information from a single protein structure. This is despite the fact that multiple structures of many protein folds are now readily available in the PDB. While ensemble protein design methods can use multiple protein structures, they treat each structure independently. Here, we introduce a flexible backbone strategy, FlexiBaL-GP, which learns global protein backbone movements directly from multiple protein structures. FlexiBaL-GP uses the machine learning method of Gaussian Process Latent Variable Models to learn a lower dimensional representation of the protein coordinates that best represent backbone movements. These learned backbone movements are used to explore alternative protein backbones, while engineering a protein within a parallel tempered MCMC framework. Using the human ubiquitin–USP21 complex as a model we demonstrate that our design strategy outperforms current strategies for the interface design task of identifying tight binding ubiquitin variants for USP21. [ABSTRACT FROM AUTHOR]

Repeated exposure to a novel physical environment eventually leads to a mature adaptive response whereby feedforward changes in motor output mirror both the amplitude and temporal structure of the environmental perturbations. However, adaptive responses at the earliest stages of learning have been found to be not only smaller, but systematically less specific in their temporal structure compared to later stages of learning. This observation has spawned a lively debate as to whether the temporal structure of the initial adaptive response is, in fact, stereotyped and non-specific. To settle this debate, we directly measured the adaptive responses to velocity-dependent and position-dependent force-field perturbations (vFFs and pFFs) at the earliest possible stage of motor learning in humans–after just a single-movement exposure. In line with previous work, we found these earliest stage adaptive responses to be more similar than the perturbations that induced them. However, the single-trial adaptive responses for vFF and pFF perturbations were clearly distinct, and the disparity between them reflected the difference between the temporal structure of the perturbations that drove them. Critically, we observed these differences between single-trial adaptive responses when vFF and pFF perturbations were randomly intermingled from one trial to the next within the same block, indicating perturbation response specificity at the single trial level. These findings demonstrate that the initial adaptive responses to physical perturbations are not stereotyped. Instead, the neural plasticity in sensorimotor areas is sensitive to the temporal structure of a movement perturbation even at the earliest stage in learning. This insight has direct implications for the development of computational models of early-stage motor adaptation and the evolution of this adaptive response with continued training. [ABSTRACT FROM AUTHOR]

Muscle fatigue is a temporary decline in the force and power capacity of skeletal muscle resulting from muscle activity. Because control of muscle is realized at the level of the motor unit (MU), it seems important to consider the physiological properties of motor units when attempting to understand and predict muscle fatigue. Therefore, we developed a phenomenological model of motor unit fatigue as a tractable means to predict muscle fatigue for a variety of tasks and to illustrate the individual contractile responses of MUs whose collective action determines the trajectory of changes in muscle force capacity during prolonged activity. An existing MU population model was used to simulate MU firing rates and isometric muscle forces and, to that model, we added fatigue-related changes in MU force, contraction time, and firing rate associated with sustained voluntary contractions. The model accurately estimated endurance times for sustained isometric contractions across a wide range of target levels. In addition, simulations were run for situations that have little experimental precedent to demonstrate the potential utility of the model to predict motor unit fatigue for more complicated, real-world applications. Moreover, the model provided insight, into the complex orchestration of MU force contributions during fatigue, that would be unattainable with current experimental approaches. [ABSTRACT FROM AUTHOR]

Motor adaptation paradigms provide a quantitative method to study short-term modification of motor commands. Despite the growing understanding of the role motion states (e.g., velocity) play in this form of motor learning, there is little information on the relative stability of memories based on these movement characteristics, especially in comparison to the initial adaptation. Here, we trained subjects to make reaching movements perturbed by force patterns dependent upon either limb position or velocity. Following training, subjects were exposed to a series of error-clamp trials to measure the temporal characteristics of the feedforward motor output during the decay of learning. The compensatory force patterns were largely based on the perturbation kinematic (e.g., velocity), but also showed a small contribution from the other motion kinematic (e.g., position). However, the velocity contribution in response to the position-based perturbation decayed at a slower rate than the position contribution to velocity-based training, suggesting a difference in stability. Next, we modified a previous model of motor adaptation to reflect this difference and simulated the behavior for different learning goals. We were interested in the stability of learning when the perturbations were based on different combinations of limb position or velocity that subsequently resulted in biased amounts of motion-based learning. We trained additional subjects on these combined motion-state perturbations and confirmed the predictions of the model. Specifically, we show that (1) there is a significant separation between the observed gain-space trajectories for the learning and decay of adaptation and (2) for combined motion-state perturbations, the gain associated to changes in limb position decayed at a faster rate than the velocity-dependent gain, even when the position-dependent gain at the end of training was significantly greater. Collectively, these results suggest that the state-dependent adaptation associated with movement velocity is relatively more stable than that based on position. [ABSTRACT FROM AUTHOR]

Nutation is an oscillatory movement that plants display during their development. Despite its ubiquity among plants movements, the relation between the observed movement and the underlying biological mechanisms remains unclear. Here we show that the kinematics of the full organ in 3D give a simple picture of plant nutation, where the orientation of the curvature along the main axis of the organ aligns with the direction of maximal differential growth. Within this framework we reexamine the validity of widely used experimental measurements of the apical tip as markers of growth dynamics. We show that though this relation is correct under certain conditions, it does not generally hold, and is not sufficient to uncover the specific role of each mechanism. As an example we re-interpret previously measured experimental observations using our model. [ABSTRACT FROM AUTHOR]

We developed a model of cardiac sarcomere contraction to study the calcium-tension relationship in cardiac muscle. Calcium mediates cardiac contraction through its interactions with troponin (Tn) and subsequently tropomyosin molecules. Experimental studies have shown that a slight increase in intracellular calcium concentration leads to a rapid increase in sarcomeric tension. Though it is widely accepted that the rapid increase is not possible without the concept of cooperativity, the mechanism is debated. We use the hypothesis that there exists a base level of cooperativity intrinsic to the thin filament that is boosted by mechanical tension, i.e. a high level of mechanical tension in the thin filament impedes the unbinding of calcium from Tn. To test these hypotheses, we developed a computational model in which a set of three parameters and inputs of calcium concentration and sarcomere length result in output tension. Tension as simulated appeared in good agreement with experimentally measured tension. Our results support the hypothesis that high tension in the thin filament impedes Tn deactivation by increasing the energy required to detach calcium from the Tn. Given this hypothesis, the model predicted that the areas with highest tension, i.e. closest to the Z-disk end of the single overlap region, show the largest concentration of active Tn’s. [ABSTRACT FROM AUTHOR]

While it is known that musculotendon units adapt to their load environments, there is only a limited understanding of tendon adaptation in vivo. Here we develop a computational model of tendon remodeling based on the premise that mechanical damage and tenocyte-mediated tendon damage and repair processes modify the distribution of its collagen fiber lengths. We explain how these processes enable the tendon to geometrically adapt to its load conditions. Based on known biological processes, mechanical and strain-dependent proteolytic fiber damage are incorporated into our tendon model. Using a stochastic model of fiber repair, it is assumed that mechanically damaged fibers are repaired longer, whereas proteolytically damaged fibers are repaired shorter, relative to their pre-damage length. To study adaptation of tendon properties to applied load, our model musculotendon unit is a simplified three-component Hill-type model of the human Achilles-soleus unit. Our model results demonstrate that the geometric equilibrium state of the Achilles tendon can coincide with minimization of the total metabolic cost of muscle activation. The proposed tendon model independently predicts rates of collagen fiber turnover that are in general agreement with in vivo experimental measurements. While the computational model here only represents a first step in a new approach to understanding the complex process of tendon remodeling in vivo, given these findings, it appears likely that the proposed framework may itself provide a useful theoretical foundation for developing valuable qualitative and quantitative insights into tendon physiology and pathology. [ABSTRACT FROM AUTHOR]

A central task in the analysis of human movement behavior is to determine systematic patterns and differences across experimental conditions, participants and repetitions. This is possible because human movement is highly regular, being constrained by invariance principles. Movement timing and movement path, in particular, are linked through scaling laws. Separating variations of movement timing from the spatial variations of movements is a well-known challenge that is addressed in current approaches only through forms of preprocessing that bias analysis. Here we propose a novel nonlinear mixed-effects model for analyzing temporally continuous signals that contain systematic effects in both timing and path. Identifiability issues of path relative to timing are overcome by using maximum likelihood estimation in which the most likely separation of space and time is chosen given the variation found in data. The model is applied to analyze experimental data of human arm movements in which participants move a hand-held object to a target location while avoiding an obstacle. The model is used to classify movement data according to participant. Comparison to alternative approaches establishes nonlinear mixed-effects models as viable alternatives to conventional analysis frameworks. The model is then combined with a novel factor-analysis model that estimates the low-dimensional subspace within which movements vary when the task demands vary. Our framework enables us to visualize different dimensions of movement variation and to test hypotheses about the effect of obstacle placement and height on the movement path. We demonstrate that the approach can be used to uncover new properties of human movement. [ABSTRACT FROM AUTHOR]

We study inter-trial movement fluctuations exhibited by human participants during the repeated execution of a virtual shuffleboard task. Focusing on skilled performance, theoretical analysis of a previously-developed general model of inter-trial error correction is used to predict the temporal and geometric structure of variability near a goal equivalent manifold (GEM). The theory also predicts that the goal-level error scales linearly with intrinsic body-level noise via the total body-goal sensitivity, a new derived quantity that illustrates how task performance arises from the interaction of active error correction and passive sensitivity properties along the GEM. Linear models estimated from observed fluctuations, together with a novel application of bootstrapping to the estimation of dynamical and correlation properties of the inter-trial dynamics, are used to experimentally confirm all predictions, thus validating our model. In addition, we show that, unlike “static” variability analyses, our dynamical approach yields results that are independent of the coordinates used to measure task execution and, in so doing, provides a new set of task coordinates that are intrinsic to the error-regulation process itself. [ABSTRACT FROM AUTHOR]

Physical performance emerges from complex interactions among many physiological systems that are largely driven by the metabolic energy demanded. Quantifying metabolic demand is an essential step for revealing the many mechanisms of physical performance decrement, but accurate predictive models do not exist. The goal of this study was to investigate if a recently developed model of muscle energetics and force could be extended to reproduce the kinematics, kinetics, and metabolic demand of submaximal effort movement. Upright dynamic knee extension against various levels of ergometer load was simulated. Task energetics were estimated by combining the model of muscle contraction with validated models of lower limb musculotendon paths and segment dynamics. A genetic algorithm was used to compute the muscle excitations that reproduced the movement with the lowest energetic cost, which was determined to be an appropriate criterion for this task. Model predictions of oxygen uptake rate (VO2) were well within experimental variability for the range over which the model parameters were confidently known. The model's accurate estimates of metabolic demand make it useful for assessing the likelihood and severity of physical performance decrement for a given task as well as investigating underlying physiologic mechanisms. [ABSTRACT FROM AUTHOR]

Humans and animals control their walking rhythms to maintain motion in a variable environment. The neural mechanism for controlling rhythm has been investigated in many studies using mechanical and electrical stimulation. However, quantitative evaluation of rhythm variation in response to perturbation at various timings has rarely been investigated. Such a characteristic of rhythm is described by the phase response curve (PRC). Dynamical simulations of human skeletal models with changing walking rhythms (phase reset) described a relation between the effective phase reset on stability and PRC, and phase reset around touch-down was shown to improve stability. A PRC of human walking was estimated by pulling the swing leg, but such perturbations hardly influenced the stance leg, so the relation between the PRC and walking events was difficult to discuss. This research thus examines human response to variations in floor velocity. Such perturbation yields another problem, in that the swing leg is indirectly (and weakly) perturbed, so the precision of PRC decreases. To solve this problem, this research adopts the weighted spike-triggered average (WSTA) method. In the WSTA method, a sequential pulsed perturbation is used for stimulation. This is in contrast with the conventional impulse method, which applies an intermittent impulsive perturbation. The WSTA method can be used to analyze responses to a large number of perturbations for each sequence. In the experiment, perturbations are applied to walking subjects by rapidly accelerating and decelerating a treadmill belt, and measured data are analyzed by the WSTA and impulse methods. The PRC obtained by the WSTA method had clear and stable waveforms with a higher temporal resolution than those obtained by the impulse method. By investigation of the rhythm transition for each phase of walking using the obtained PRC, a rhythm change that extends the touch-down and mid-single support phases is found to occur. [ABSTRACT FROM AUTHOR]

Humans employ a high degree of redundancy in joint actuation, with different combinations of muscle and tendon action providing the same net joint torque. Both the resolution of these redundancies and the energetics of such systems depend on the dynamic properties of muscles and tendons, particularly their force-length relations. Current walking models that use stock parameters when simulating muscle-tendon dynamics tend to significantly overestimate metabolic consumption, perhaps because they do not adequately consider the role of elasticity. As an alternative, we posit that the muscle-tendon morphology of the human leg has evolved to maximize the metabolic efficiency of walking at self-selected speed. We use a data-driven approach to evaluate this hypothesis, utilizing kinematic, kinetic, electromyographic (EMG), and metabolic data taken from five participants walking at self-selected speed. The kinematic and kinetic data are used to estimate muscle-tendon lengths, muscle moment arms, and joint moments while the EMG data are used to estimate muscle activations. For each subject we perform an optimization using prescribed skeletal kinematics, varying the parameters that govern the force-length curve of each tendon as well as the strength and optimal fiber length of each muscle while seeking to simultaneously minimize metabolic cost and maximize agreement with the estimated joint moments. We find that the metabolic cost of transport (MCOT) values of our participants may be correctly matched (on average 0.36±0.02 predicted, 0.35±0.02 measured) with acceptable joint torque fidelity through application of a single constraint to the muscle metabolic budget. The associated optimal muscle-tendon parameter sets allow us to estimate the forces and states of individual muscles, resolving redundancies in joint actuation and lending insight into the potential roles and control objectives of the muscles of the leg throughout the gait cycle. [ABSTRACT FROM AUTHOR]

The laws of physics establish the energetic efficiency of our movements. In some cases, like locomotion, the mechanics of the body dominate in determining the energetically optimal course of action. In other tasks, such as manipulation, energetic costs depend critically upon the variable properties of objects in the environment. Can the brain identify and follow energy-optimal motions when these motions require moving along unfamiliar trajectories? What feedback information is required for such optimal behavior to occur? To answer these questions, we asked participants to move their dominant hand between different positions while holding a virtual mechanical system with complex dynamics (a planar double pendulum). In this task, trajectories of minimum kinetic energy were along curvilinear paths. Our findings demonstrate that participants were capable of finding the energy-optimal paths, but only when provided with veridical visual and haptic information pertaining to the object, lacking which the trajectories were executed along rectilinear paths. [ABSTRACT FROM AUTHOR]