Transforming modeling in neurorehabilitation: clinical insights for personalized rehabilitation.

Practicing clinicians in neurorehabilitation continue to lack a systematic evidence base to personalize rehabilitation therapies to individual patients and thereby maximize outcomes. Computational modeling- collecting, analyzing, and modeling neurorehabilitation data- holds great promise. A key question is how can computational modeling contribute to the evidence base for personalized rehabilitation? As representatives of the clinicians and clinician-scientists who attended the 2023 NSF DARE conference at USC, here we offer our perspectives and discussion on this topic. Our overarching thesis is that clinical insight should inform all steps of modeling, from construction to output, in neurorehabilitation and that this process requires close collaboration between researchers and the clinical community. We start with two clinical case examples focused on motor rehabilitation after stroke which provide context to the heterogeneity of neurologic injury, the complexity of post-acute neurologic care, the neuroscience of recovery, and the current state of outcome assessment in rehabilitation clinical care. Do we provide different therapies to these two different patients to maximize outcomes? Asking this question leads to a corollary: how do we build the evidence base to support the use of different therapies for individual patients? We discuss seven points critical to clinical translation of computational modeling research in neurorehabilitation- (i) clinical endpoints, (ii) hypothesis- versus data-driven models, (iii) biological processes, (iv) contextualizing outcome measures, (v) clinical collaboration for device translation, (vi) modeling in the real world and (vii) clinical touchpoints across all stages of research. We conclude with our views on key avenues for future investment (clinical-research collaboration, new educational pathways, interdisciplinary engagement) to enable maximal translational value of computational modeling research in neurorehabilitation.

NSF DARE-Transforming modeling in neurorehabilitation: Four threads for catalyzing progress.

We present an overview of the Conference on Transformative Opportunities for Modeling in Neurorehabilitation held in March 2023. It was supported by the Disability and Rehabilitation Engineering (DARE) program from the National Science Foundation's Engineering Biology and Health Cluster. The conference brought together experts and trainees from around the world to discuss critical questions, challenges, and opportunities at the intersection of computational modeling and neurorehabilitation to understand, optimize, and improve clinical translation of neurorehabilitation. We organized the conference around four key, relevant, and promising Focus Areas for modeling: Adaptation & Plasticity, Personalization, Human-Device Interactions, and Modeling 'In-the-Wild'. We identified four common threads across the Focus Areas that, if addressed, can catalyze progress in the short, medium, and long terms. These were: (i) the need to capture and curate appropriate and useful data necessary to develop, validate, and deploy useful computational models (ii) the need to create multi-scale models that span the personalization spectrum from individuals to populations, and from cellular to behavioral levels (iii) the need for algorithms that extract as much information from available data, while requiring as little data as possible from each client (iv) the insistence on leveraging readily available sensors and data systems to push model-driven treatments from the lab, and into the clinic, home, workplace, and community. The conference archive can be found at (dare2023.usc.edu). These topics are also extended by three perspective papers prepared by trainees and junior faculty, clinician researchers, and federal funding agency representatives who attended the conference.

Dynamical Analyses Show That Professional Archers Exhibit Tighter, Finer and More Fluid Dynamical Control Than Neophytes.

Quantifying the dynamical features of discrete tasks is essential to understanding athletic performance for many sports that are not repetitive or cyclical. We compared three dynamical features of the (i) bow hand, (ii) drawing hand, and (iii) center of mass during a single bow-draw movement between professional and neophyte archers: dispersion (convex hull volume of their phase portraits), persistence (tendency to continue a trend as per Hurst exponents), and regularity (sample entropy). Although differences in the two groups are expected due to their differences in skill, our results demonstrate we can quantify these differences. The center of mass of professional athletes exhibits tighter movements compared to neophyte archers (6.3 < 11.2 convex hull volume), which are nevertheless less persistent (0.82 < 0.86 Hurst exponent) and less regular (0.035 > 0.025 sample entropy). In particular, the movements of the bow hand and center of mass differed more between groups in Hurst exponent analysis, and the drawing hand and center of mass were more different in sample entropy analysis. This suggests tighter neuromuscular control over the more fluid dynamics of the movement that exhibits more active corrections that are more individualized. Our work, therefore, provides proof of principle of how well-established dynamical analysis techniques can be used to quantify the nature and features of neuromuscular expertise for discrete movements in elite athletes.

Force variability is mostly not motor noise: Theoretical implications for motor control.

Variability in muscle force is a hallmark of healthy and pathological human behavior. Predominant theories of sensorimotor control assume 'motor noise' leads to force variability and its 'signal dependence' (variability in muscle force whose amplitude increases with intensity of neural drive). Here, we demonstrate that the two proposed mechanisms for motor noise (i.e. the stochastic nature of motor unit discharge and unfused tetanic contraction) cannot account for the majority of force variability nor for its signal dependence. We do so by considering three previously underappreciated but physiologically important features of a population of motor units: 1) fusion of motor unit twitches, 2) coupling among motoneuron discharge rate, cross-bridge dynamics, and muscle mechanics, and 3) a series-elastic element to account for the aponeurosis and tendon. These results argue strongly against the idea that force variability and the resulting kinematic variability are generated primarily by 'motor noise.' Rather, they underscore the importance of variability arising from properties of control strategies embodied through distributed sensorimotor systems. As such, our study provides a critical path toward developing theories and models of sensorimotor control that provide a physiologically valid and clinically useful understanding of healthy and pathologic force variability.

Temporal control of muscle synergies is linked with alpha-band neural drive.

KEY POINTS: It is theorized that the nervous system controls groups of muscles together as functional units, or 'synergies', resulting in correlated electromyographic (EMG) signals among muscles. However, such correlation does not necessarily imply group-level neural control. Oscillatory synchronization (coherence) among EMG signals implies neural coupling, but it is not clear how this relates to control of muscle synergies. EMG was recorded from seven arm muscles of 10 adult participants rotating an upper limb ergometer, and EMG-EMG coherence, EMG amplitude correlations and their relationship with each other were characterized. A novel method to derive multi-muscle synergies from EMG-EMG coherence is presented and these are compared with classically defined synergies. Coherent alpha-band (8-16 Hz) drive was strongest among muscles whose gross activity levels are well correlated within a given task. The cross-muscle distribution and temporal modulation of coherent alpha-band drive suggests a possible role in the neural coordination/monitoring of synergies. ABSTRACT: During movement, groups of muscles may be controlled together by the nervous system as an adaptable functional entity, or 'synergy'. The rules governing when (or if) this occurs during voluntary behaviour in humans are not well understood, at least in part because synergies are usually defined by correlated patterns of muscle activity without regard for the underlying structure of their neural control. In this study, we investigated the extent to which comodulation of muscle output (i.e. correlation of electromyographic (EMG) amplitudes) implies that muscles share intermuscular neural input (assessed via EMG-EMG coherence analysis). We first examined this relationship among pairs of upper limb muscles engaged in an arm cycling task. We then applied a novel multidimensional EMG-EMG coherence analysis allowing synergies to be characterized on the basis of shared neural drive. We found that alpha-band coherence (8-16 Hz) is related to the degree to which overall muscle activity levels correlate over time. The extension of this coherence analysis to describe the cross-muscle distribution and temporal modulation of alpha-band drive revealed a close match to the temporal and structural features of traditionally defined muscle synergies. Interestingly, the coherence-derived neural drive was inversely associated with, and preceded, changes in EMG amplitudes by ∼200 ms. Our novel characterization of how alpha-band neural drive is dynamically distributed among muscles is a fundamental step forward in understanding the neural origins and correlates of muscle synergies.

Cardinal features of involuntary force variability can arise from the closed-loop control of viscoelastic afferented muscles.

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.

Physiological tremor increases when skeletal muscle is shortened: implications for fusimotor control.

KEY POINTS: In tonic, isometric, plantarflexion contractions, physiological tremor increases as the ankle joint becomes plantarflexed. Modulation of physiological tremor as a function of muscle stretch differs from that of the stretch reflex amplitude. Amplitude of physiological tremor may be altered as a function of reflex pathway gains. Healthy humans likely increase their γ-static fusimotor drive when muscles shorten. Quantification of physiological tremor by manipulation of joint angle may be a useful experimental probe of afferent gains and/or the integrity of automatic fusimotor control. ABSTRACT: The involuntary force fluctuations associated with physiological (as distinct from pathological) tremor are an unavoidable component of human motor control. While the origins of physiological tremor are known to depend on muscle afferentation, it is possible that the mechanical properties of muscle-tendon systems also affect its generation, amplification and maintenance. In this paper, we investigated the dependence of physiological tremor on muscle length in healthy individuals. We measured physiological tremor during tonic, isometric plantarflexion torque at 30% of maximum at three ankle angles. The amplitude of physiological tremor increased as calf muscles shortened in contrast to the stretch reflex whose amplitude decreases as muscle shortens. We used a published closed-loop simulation model of afferented muscle to explore the mechanisms responsible for this behaviour. We demonstrate that changing muscle lengths does not suffice to explain our experimental findings. Rather, the model consistently required the modulation of  Î³-static fusimotor drive to produce increases in physiological tremor with muscle shortening - while successfully replicating the concomitant reduction in stretch reflex amplitude. This need to control γ-static fusimotor drive explicitly as a function of muscle length has important implications. First, it permits the amplitudes of physiological tremor and stretch reflex to be decoupled. Second, it postulates neuromechanical interactions that require length-dependent γ drive modulation to be independent from α drive to the parent muscle. Lastly, it suggests that physiological tremor can be used as a simple, non-invasive measure of the afferent mechanisms underlying healthy motor function, and their disruption in neurological conditions.

Intermuscular coherence reflects functional coordination.

Coherence analysis has the ability to identify the presence of common descending drive shared by motor unit pools and reveals its spectral properties. However, the link between spectral properties of shared neural drive and functional interactions among muscles remains unclear. We assessed shared neural drive between muscles of the thumb and index finger while participants executed two mechanically distinct precision pinch tasks, each requiring distinct functional coordination among muscles. We found that shared neural drive was systematically reduced or enhanced at specific frequencies of interest (~10 and ~40 Hz). While amplitude correlations between surface EMG signals also exhibited changes across tasks, only their coherence has strong physiological underpinnings indicative of neural binding. Our results support the use of intermuscular coherence as a tool to detect when coactivated muscles are members of a functional group or synergy of neural origin. Furthermore, our results demonstrate the advantages of considering neural binding at 10, ~20, and >30 Hz, as indicators of task-dependent neural coordination strategies.NEW & NOTEWORTHY It is often unclear whether correlated activity among muscles reflects their neural binding or simply reflects the constraints defining the task. Using the fact that high-frequency coherence between EMG signals (>6 Hz) is thought to reflect shared neural drive, we demonstrate that coherence analysis can reveal the neural origin of distinct muscle coordination patterns required by different tasks.

Muscle Synergies Heavily Influence the Neural Control of Arm Endpoint Stiffness and Energy Consumption.

Much debate has arisen from research on muscle synergies with respect to both limb impedance control and energy consumption. Studies of limb impedance control in the context of reaching movements and postural tasks have produced divergent findings, and this study explores whether the use of synergies by the central nervous system (CNS) can resolve these findings and also provide insights on mechanisms of energy consumption. In this study, we phrase these debates at the conceptual level of interactions between neural degrees of freedom and tasks constraints. This allows us to examine the ability of experimentally-observed synergies--correlated muscle activations--to control both energy consumption and the stiffness component of limb endpoint impedance. In our nominal 6-muscle planar arm model, muscle synergies and the desired size, shape, and orientation of endpoint stiffness ellipses, are expressed as linear constraints that define the set of feasible muscle activation patterns. Quadratic programming allows us to predict whether and how energy consumption can be minimized throughout the workspace of the limb given those linear constraints. We show that the presence of synergies drastically decreases the ability of the CNS to vary the properties of the endpoint stiffness and can even preclude the ability to minimize energy. Furthermore, the capacity to minimize energy consumption--when available--can be greatly affected by arm posture. Our computational approach helps reconcile divergent findings and conclusions about task-specific regulation of endpoint stiffness and energy consumption in the context of synergies. But more generally, these results provide further evidence that the benefits and disadvantages of muscle synergies go hand-in-hand with the structure of feasible muscle activation patterns afforded by the mechanics of the limb and task constraints. These insights will help design experiments to elucidate the interplay between synergies and the mechanisms of learning, plasticity, versatility and pathology in neuromuscular systems.

On neuromechanical approaches for the study of biological and robotic grasp and manipulation.

Biological and robotic grasp and manipulation are undeniably similar at the level of mechanical task performance. However, their underlying fundamental biological vs. engineering mechanisms are, by definition, dramatically different and can even be antithetical. Even our approach to each is diametrically opposite: inductive science for the study of biological systems vs. engineering synthesis for the design and construction of robotic systems. The past 20 years have seen several conceptual advances in both fields and the quest to unify them. Chief among them is the reluctant recognition that their underlying fundamental mechanisms may actually share limited common ground, while exhibiting many fundamental differences. This recognition is particularly liberating because it allows us to resolve and move beyond multiple paradoxes and contradictions that arose from the initial reasonable assumption of a large common ground. Here, we begin by introducing the perspective of neuromechanics, which emphasizes that real-world behavior emerges from the intimate interactions among the physical structure of the system, the mechanical requirements of a task, the feasible neural control actions to produce it, and the ability of the neuromuscular system to adapt through interactions with the environment. This allows us to articulate a succinct overview of a few salient conceptual paradoxes and contradictions regarding under-determined vs. over-determined mechanics, under- vs. over-actuated control, prescribed vs. emergent function, learning vs. implementation vs. adaptation, prescriptive vs. descriptive synergies, and optimal vs. habitual performance. We conclude by presenting open questions and suggesting directions for future research. We hope this frank and open-minded assessment of the state-of-the-art will encourage and guide these communities to continue to interact and make progress in these important areas at the interface of neuromechanics, neuroscience, rehabilitation and robotics.

Robot-assisted and conventional therapies produce distinct rehabilitative trends in stroke survivors.

BACKGROUND: Comparing the efficacy of alternative therapeutic strategies for the rehabilitation of motor function in chronically impaired individuals is often inconclusive. For example, a recent randomized clinical trial (RCT) compared robot-assisted vs. conventional therapy in 77 patients who had had chronic motor impairment after a cerebrovascular accident. While patients assigned to robotic therapy had greater improvements in the primary outcome measure (change in score on the upper extremity section of the Fugl-Meyer assessment), the absolute difference between therapies was small, which left the clinical relevance in question. METHODS: Here we revisit that study to test whether the multidimensional rehabilitative response of these patients can better distinguish between treatment outcomes. We used principal components analysis to find the correlation of changes across seven outcome measures between the start and end of 8 weeks of therapy. Permutation tests verified the robustness of the principal components found. RESULTS: Each therapy in fact produces different rehabilitative trends of recovery across the clinical, functional, and quality of life domains. A rehabilitative trend is a principal component that quantifies the correlations among changes in outcomes with each therapy. CONCLUSIONS: These findings challenge the traditional emphasis of RCTs on using a single primary outcome measure to compare rehabilitative responses that are naturally multidimensional. This alternative approach to, and interpretation of, the results of RCTs may will lead to more effective therapies targeted for the multidimensional mechanisms of recovery. TRIAL REGISTRATION: ClinicalTrials.gov number NCT00719433 . Registered July 17, 2008.

Innovative evaluation of dexterity in pediatrics.

STUDY DESIGN: Review paper. INTRODUCTION: Hand dexterity is multifaceted and essential to the performance of daily tasks. Timed performance and precision demands are the most common features of quantitative dexterity testing. Measurement concepts such as rate of completion, in-hand manipulation and dynamic force control of instabilities are being integrated into assessment tools for the pediatric population. PURPOSE: To review measurement concepts inherent in pediatric dexterity testing and introduce concepts that are infrequently measured or novel as exemplified with two assessment tools. METHODS: Measurement concepts included in common assessment tools are introduced first. We then describe seldom measured and novel concepts embedded in two instruments; the Functional Dexterity Test (FDT) and the Strength-Dexterity (SD) Test. DISCUSSION: The inclusion of novel yet informative tools and measurement concepts in our assessments could aid our understanding of atypical dexterity, and potentially contribute to the design of targeted therapy programs.

Long term functional outcomes after early childhood pollicization.

STUDY DESIGN: Retrospective Cohort INTRODUCTION: Important outcomes of polliciation to treat thumb hypoplasia/aplasia include strength, function, dexterity, and quality of life. PURPOSE OF THE STUDY: To evaluate outcomes and examine predictors of outcome after early childhood pollicization. METHODS: 8 children (10 hands) were evaluated 3-15 years after surgery. Physical examination, questionnaires, grip and pinch strength, Box and Blocks, 9-hole pegboard, and strength-dexterity (S-D) tests were performed. RESULTS: Pollicized hands had poor strength and performance on functional tests. Six of 10 pollicized hands had normal dexterity scores but less stability in maintaining a steady-state force. Predictors of poorer outcomes included older age at surgery, reduced metacarpophalangeal and interphalangeal range of motion, and radial absence. DISCUSSION: Pollicization resulted in poor strength and overall function, but normal dexterity was often achieved using altered control strategies. CONCLUSIONS: Most children should obtain adequate dexterity despite weakness after pollicization except older or severely involved children. LEVEL OF EVIDENCE: IV.

Decrease in muscle contraction time complements neural maturation in the development of dynamic manipulation.

Developmental improvements in dynamic manipulation abilities are typically attributed to neural maturation, such as myelination of corticospinal pathways, neuronal pruning, and synaptogenesis. However the contributions from changes in the peripheral motor system are less well understood. Here we investigated whether there are developmental changes in muscle activation-contraction dynamics and whether these changes contribute to improvements in dynamic manipulation in humans. We compared pinch strength, dynamic manipulation ability, and contraction time of the first dorsal interosseous muscle in typically developing preadolescent, adolescent, and young adults. Both strength and dynamic manipulation ability increased with age (p < 0.0001 and p < 0.00001, respectively). Surprisingly, adults had a 33% lower muscle contraction time compared with preadolescents (p < 0.01), and contraction time showed a significant (p < 0.005) association with dynamic manipulation abilities. Whereas decreases in muscle contraction time during development have been reported in the animal literature, our finding, to our knowledge, is the first report of this phenomenon in humans and the first finding of its association with manipulation. Consequently, the changes in the muscle contractile properties could be an important complement to neural maturation in the development of dynamic manipulation. These findings have important implications for understanding central and peripheral contributors to deficits in manipulation in atypical development, such as in children with cerebral palsy.

Generalizability of foot placement control strategies during unperturbed and perturbed gait.

Control of foot placement is an essential strategy for maintaining balance during walking. During unperturbed, steady-state walking, foot placement can be accurately described as a linear function of the body's centre of mass (CoM) state at midstance. However, it is uncertain if this mapping from CoM state to foot placement generalizes to larger perturbations that could potentially cause falls. Recovery from these perturbations may require reactive control strategies not observed during unperturbed walking. Here, we used unpredictable changes in treadmill belt speed to assess the generalizability of foot placement mappings identified during unperturbed walking. We found that foot placement mappings generalized poorly from unperturbed to perturbed walking and differed for forward perturbation versus backward perturbation. We also used the singular value decomposition of the mapping matrix to reveal that people were more sensitive to backward versus forward perturbations. Together, these results indicate that a single linear mapping cannot describe the foot placement control during both forward and backward losses of balance induced by treadmill belt speed perturbations. Better characterization of human balance control strategies could improve our understanding of why different neuromotor disorders result in heightened fall risk and inform the design of controllers for balance-assisting devices.

An alpha- to gamma-motoneurone collateral can mitigate velocity-dependent stretch reflexes during voluntary movement: A computational study.

The primary motor cortex does not uniquely or directly produce alpha motoneurone (α-MN) drive to muscles during voluntary movement. Rather, α-MN drive emerges from the synthesis and competition among excitatory and inhibitory inputs from multiple descending tracts, spinal interneurons, sensory inputs, and proprioceptive afferents. One such fundamental input is velocity-dependent stretch reflexes in lengthening muscles, which should be inhibited to enable voluntary movement. It remains an open question, however, the extent to which unmodulated stretch reflexes disrupt voluntary movement, and whether and how they are inhibited in limbs with numerous multi-articular muscles. We used a computational model of a Rhesus Macaque arm to simulate movements with feedforward α-MN commands only, and with added velocity-dependent stretch reflex feedback. We found that velocity-dependent stretch reflex caused movement-specific, typically large and variable disruptions to arm movements. These disruptions were greatly reduced when modulating velocity-dependent stretch reflex feedback (i) as per the commonly proposed (but yet to be clarified) idealized alpha-gamma (α-γ) co-activation or (ii) an alternative α-MN collateral projection to homonymous γ-MNs. We conclude that such α-MN collaterals are a physiologically tenable, but previously unrecognized, propriospinal circuit in the mammalian fusimotor system. These collaterals could still collaborate with α-γ co-activation, and the few skeletofusimotor fibers (ß-MNs) in mammals, to create a flexible fusimotor ecosystem to enable voluntary movement. By locally and automatically regulating the highly nonlinear neuro-musculo-skeletal mechanics of the limb, these collaterals could be a critical low-level enabler of learning, adaptation, and performance via higher-level brainstem, cerebellar and cortical mechanisms.

Developmental improvements in dynamic control of fingertip forces last throughout childhood and into adolescence.

While it is clear that the development of dexterous manipulation in children exhibits dramatic improvements over an extended period, it is difficult to separate musculoskeletal from neural contributors to these important functional gains. This is in part due to the inability of current methods to disambiguate improvements in hand strength from gains in finger dexterity (i.e., the dynamic control of fingertip force vectors at low magnitudes). We adapted our novel instrumentation to evaluate finger dexterity in 130 typically developing children between the ages of 4 and 16 yr. We find that finger dexterity continues to develop well into late adolescence and musculoskeletal growth and strength are poorly correlated with the improvements in dexterity. Importantly, because these behavioral results seem to mirror the known timelines of neuroanatomical development up to adolescence, we speculate that they reflect the functional benefits of such continual neural maturation. This novel perspective now enables the systematic study of the functional roles of specific neuroanatomical structures and their connectivity, maturity, and plasticity. Moreover, the temporal dynamics of the fingertip force vectors shows improvements in stability that provide a novel way to look at the maturation of finger control. From a clinical perspective, our results provide a practical means to chart functional development of dexterous manipulation in typically developing children and could be adapted for clinical use and for use in children with developmental disorders.

Functional inference of complex anatomical tendinous networks at a macroscopic scale via sparse experimentation.

In systems and computational biology, much effort is devoted to functional identification of systems and networks at the molecular-or cellular scale. However, similarly important networks exist at anatomical scales such as the tendon network of human fingers: the complex array of collagen fibers that transmits and distributes muscle forces to finger joints. This network is critical to the versatility of the human hand, and its function has been debated since at least the 16(th) century. Here, we experimentally infer the structure (both topology and parameter values) of this network through sparse interrogation with force inputs. A population of models representing this structure co-evolves in simulation with a population of informative future force inputs via the predator-prey estimation-exploration algorithm. Model fitness depends on their ability to explain experimental data, while the fitness of future force inputs depends on causing maximal functional discrepancy among current models. We validate our approach by inferring two known synthetic Latex networks, and one anatomical tendon network harvested from a cadaver's middle finger. We find that functionally similar but structurally diverse models can exist within a narrow range of the training set and cross-validation errors. For the Latex networks, models with low training set error [<4%] and resembling the known network have the smallest cross-validation errors [∼5%]. The low training set [<4%] and cross validation [<7.2%] errors for models for the cadaveric specimen demonstrate what, to our knowledge, is the first experimental inference of the functional structure of complex anatomical networks. This work expands current bioinformatics inference approaches by demonstrating that sparse, yet informative interrogation of biological specimens holds significant computational advantages in accurate and efficient inference over random testing, or assuming model topology and only inferring parameters values. These findings also hold clues to both our evolutionary history and the development of versatile machines.

Challenges and new approaches to proving the existence of muscle synergies of neural origin.

Muscle coordination studies repeatedly show low-dimensionality of muscle activations for a wide variety of motor tasks. The basis vectors of this low-dimensional subspace, termed muscle synergies, are hypothesized to reflect neurally-established functional muscle groupings that simplify body control. However, the muscle synergy hypothesis has been notoriously difficult to prove or falsify. We use cadaveric experiments and computational models to perform a crucial thought experiment and develop an alternative explanation of how muscle synergies could be observed without the nervous system having controlled muscles in groups. We first show that the biomechanics of the limb constrains musculotendon length changes to a low-dimensional subspace across all possible movement directions. We then show that a modest assumption--that each muscle is independently instructed to resist length change--leads to the result that electromyographic (EMG) synergies will arise without the need to conclude that they are a product of neural coupling among muscles. Finally, we show that there are dimensionality-reducing constraints in the isometric production of force in a variety of directions, but that these constraints are more easily controlled for, suggesting new experimental directions. These counter-examples to current thinking clearly show how experimenters could adequately control for the constraints described here when designing experiments to test for muscle synergies--but, to the best of our knowledge, this has not yet been done.

Generalizability of foot-placement control strategies during unperturbed and perturbed gait.

Control of foot placement is an essential strategy for maintaining balance during walking. During unperturbed, steady-state walking, foot placement can be accurately described as a linear function of the body's center of mass state at midstance. However, it is uncertain if this mapping from center of mass state to foot placement generalizes to larger perturbations that may be more likely to cause falls. These perturbations may cause balance disturbances and generate reactive control strategies not observed during unperturbed walking. Here, we used unpredictable changes in treadmill speed to assess the generalizability of foot placement mappings identified during unperturbed walking. We found that foot placement mappings generalized poorly from unperturbed to perturbed walking and differed for forward versus backward perturbations. We also used singular value decomposition of the mapping matrix to reveal that people were more sensitive to backward versus forward perturbations. Together, these results indicate that control of foot placement during losses of balance differs from the control strategies used during unperturbed walking. Better characterization of human balance control strategies could improve our understanding of why different neuromotor disorders result in heightened fall risk and inform the design of controllers for balance-assisting devices.