Computational tools for calculating alternative muscle force patterns during motion: A comparison of possible solutions

Modeling Human Suboptimal Control: A Review

This review paper provides an overview of the approaches to model neuromuscular control, focusing on methods to identify nonoptimal control strategies typical of populations with neuromuscular disorders or children. Where possible, the authors tightened the description of the methods to the mechanisms behind the underlying biomechanical and physiological rationale. They start by describing the first and most simplified approach, the reductionist approach, which splits the role of the nervous and musculoskeletal systems. Static optimization and dynamic optimization methods and electromyography-based approaches are summarized to highlight their limitations and understand (the need for) their developments over time. Then, the authors look at the more recent stochastic approach, introduced to explore the space of plausible neural solutions, thus implementing the uncontrolled manifold theory, according to which the central nervous system only controls specific motions and tasks to limit energy consumption while allowing for some degree of adaptability to perturbations. Finally, they explore the literature covering the explicit modeling of the coupling between the nervous system (acting as controller) and the musculoskeletal system (the actuator), which may be employed to overcome the split characterizing the reductionist approach.

Neuromuscular control: from a biomechanist's perspective

Understanding neural control of movement necessitates a collaborative approach between many disciplines, including biomechanics, neuroscience, and motor control. Biomechanics grounds us to the laws of physics that our musculoskeletal system must obey. Neuroscience reveals the inner workings of our nervous system that functions to control our body. Motor control investigates the coordinated motor behaviours we display when interacting with our environment. The combined efforts across the many disciplines aimed at understanding human movement has resulted in a rich and rapidly growing body of literature overflowing with theories, models, and experimental paradigms. As a result, gathering knowledge and drawing connections between the overlapping but seemingly disparate fields can be an overwhelming endeavour. This review paper evolved as a need for us to learn of the diverse perspectives underlying current understanding of neuromuscular control. The purpose of our review paper is to integrate ideas from biomechanics, neuroscience, and motor control to better understand how we voluntarily control our muscles. As biomechanists, we approach this paper starting from a biomechanical modelling framework. We first define the theoretical solutions (i.e., muscle activity patterns) that an individual could feasibly use to complete a motor task. The theoretical solutions will be compared to experimental findings and reveal that individuals display structured muscle activity patterns that do not span the entire theoretical solution space. Prevalent neuromuscular control theories will be discussed in length, highlighting optimality, probabilistic principles, and neuromechanical constraints, that may guide individuals to families of muscle activity solutions within what is theoretically possible. Our intention is for this paper to serve as a primer for the neuromuscular control scientific community by introducing and integrating many of the ideas common across disciplines today, as well as inspire future work to improve the representation of neural control in biomechanical models.

Electroplating of HAp-brushite coating on metallic bioimplants with advanced hemocompatibility and osteocompatibility properties

In cases of severe bone tissue injuries, the use of metallic bioimplants is quite widespread due to their high strength, high fracture toughness, hardness, and corrosion resistance. However, they lack adequate biocompatibility and show poor metal-tissue integration during the post-operative phase. To mitigate this drawback, it is beneficial to add a biocompatible polymer layer to ensure a quick growth of cell or tissue over the surface of metallic bioimplant material. Furthermore, this additional layer should possess good adherence with the underlying material and also accompany a rapid bonding between the tissue and the implant material, in order to reduce the recovery time for the patient. Therefore, in this work, we report a novel green electroplating route for growing porous hydroxyapatite-brushite coatings on a stainless steel surface. The malic acid used for the production of hydroxyapatite-brushite coatings has been obtained from an extract of locally available apple fruit (Malus domestica). We demonstrate the effect of electroplating parameters on the structural morphology of the electroplated composite layer via XRD, SEM with EDS, and FTIR characterization techniques and report an optimized set of electroplating parameters that will yield the best composite coating in terms of thickness, adherence to substrate and speed. The hemocompatibility and osteocompatibility studies on the electroplated composites coating show this technology's effectiveness and potential applicability in biomedical applications. Compared to other routes reported in the literature, this electroplating route is quicker and yields better composite coatings with faster bone tissue growth potential.

Effects of kinematic complexity and number of muscles on musculoskeletal model robustness to muscle dysfunction

The robustness of motor outputs to muscle dysfunction has been investigated using musculoskeletal modeling, but with conflicting results owing to differences in model complexity and motor tasks. Our objective was to systematically study how the number of kinematic degrees of freedom, and the number of independent muscle actuators alter the robustness of motor output to muscle dysfunction. We took a detailed musculoskeletal model of the human leg and systematically varied the model complexity to create six models with either 3 or 7 kinematic degrees of freedom and either 14, 26, or 43 muscle actuators. We tested the redundancy of each model by quantifying the reduction in sagittal plane feasible force set area when a single muscle was removed. The robustness of feasible force set area to the loss of any single muscle, i.e. general single muscle loss increased with the number of independent muscles and decreased with the number of kinematic degrees of freedom, with the robust area varying from 1% and 52% of the intact feasible force set area. The maximum sensitivity of the feasible force set to the loss of any single muscle varied from 75% to 26% of the intact feasible force set area as the number of muscles increased. Additionally, the ranges of feasible muscle activation for maximum force production were largely unconstrained in many cases, indicating ample musculoskeletal redundancy even for maximal forces. We propose that ratio of muscles to kinematic degrees of freedom can be used as a rule of thumb for estimating musculoskeletal redundancy in both simulated and real biomechanical systems.

In silico assessment of biomedical products: The conundrum of rare but not so rare events in two case studies

In silico clinical trials, defined as "The use of individualized computer simulation in the development or regulatory evaluation of a medicinal product, medical device, or medical intervention," have been proposed as a possible strategy to reduce the regulatory costs of innovation and the time to market for biomedical products. We review some of the the literature on this topic, focusing in particular on those applications where the current practice is recognized as inadequate, as for example, the detection of unexpected severe adverse events too rare to be detected in a clinical trial, but still likely enough to be of concern. We then describe with more details two case studies, two successful applications of in silico clinical trial approaches, one relative to the University of Virginia/Padova simulator that the Food and Drug Administration has accepted as possible replacement for animal testing in the preclinical assessment of artificial pancreas technologies, and the second, an investigation of the probability of cardiac lead fracture, where a Bayesian network was used to combine in vivo and in silico observations, suggesting a whole new strategy of in silico-augmented clinical trials, to be used to increase the numerosity where recruitment is impossible, or to explore patients' phenotypes that are unlikely to appear in the trial cohort, but are still frequent enough to be of concern.

The Virtual Physiological Human: Ten Years After

Biomedical research and clinical practice are struggling to cope with the growing complexity that the progress of health care involves. The most challenging diseases, those with the largest socioeconomic impact (cardiovascular conditions; musculoskeletal conditions; cancer; metabolic, immunity, and neurodegenerative conditions), are all characterized by a complex genotype-phenotype interaction and by a "systemic" nature that poses a challenge to the traditional reductionist approach. In 2005 a small group of researchers discussed how the vision of computational physiology promoted by the Physiome Project could be translated into clinical practice and formally proposed the term Virtual Physiological Human. Our knowledge about these diseases is fragmentary, as it is associated with molecular and cellular processes on the one hand and with tissue and organ phenotype changes (related to clinical symptoms of disease conditions) on the other. The problem could be solved if we could capture all these fragments of knowledge into predictive models and then compose them into hypermodels that help us tame the complexity that such systemic behavior involves. In 2005 this was simply not possible-the necessary methods and technologies were not available. Now, 10 years later, it seems the right time to reflect on the original vision, the results achieved so far, and what remains to be done.

Muscle parameters estimation based on biplanar radiography

The evaluation of muscle and joint forces in vivo is still a challenge. Musculo-Skeletal (musculo-skeletal) models are used to compute forces based on movement analysis. Most of them are built from a scaled-generic model based on cadaver measurements, which provides a low level of personalization, or from Magnetic Resonance Images, which provide a personalized model in lying position. This study proposed an original two steps method to access a subject-specific musculo-skeletal model in 30 min, which is based solely on biplanar X-Rays. First, the subject-specific 3D geometry of bones and skin envelopes were reconstructed from biplanar X-Rays radiography. Then, 2200 corresponding control points were identified between a reference model and the subject-specific X-Rays model. Finally, the shape of 21 lower limb muscles was estimated using a non-linear transformation between the control points in order to fit the muscle shape of the reference model to the X-Rays model. Twelfth musculo-skeletal models were reconstructed and compared to their reference. The muscle volume was not accurately estimated with a standard deviation (SD) ranging from 10 to 68%. However, this method provided an accurate estimation the muscle line of action with a SD of the length difference lower than 2% and a positioning error lower than 20 mm. The moment arm was also well estimated with SD lower than 15% for most muscle, which was significantly better than scaled-generic model for most muscle. This method open the way to a quick modeling method for gait analysis based on biplanar radiography.

Suboptimal Muscle Synergy Activation Patterns Generalize their Motor Function across Postures

We used a musculoskeletal model to investigate the possible biomechanical and neural bases of using consistent muscle synergy patterns to produce functional motor outputs across different biomechanical conditions, which we define as generalizability. Experimental studies in cats demonstrate that the same muscle synergies are used during reactive postural responses at widely varying configurations, producing similarly-oriented endpoint force vectors with respect to the limb axis. However, whether generalizability across postures arises due to similar biomechanical properties or to neural selection of a particular muscle activation pattern has not been explicitly tested. Here, we used a detailed cat hindlimb model to explore the set of feasible muscle activation patterns that produce experimental synergy force vectors at a target posture, and tested their generalizability by applying them to different test postures. We used three methods to select candidate muscle activation patterns: (1) randomly-selected feasible muscle activation patterns, (2) optimal muscle activation patterns minimizing muscle effort at a given posture, and (3) generalizable muscle activation patterns that explicitly minimized deviations from experimentally-identified synergy force vectors across all postures. Generalizability was measured by the deviation between the simulated force direction of the candidate muscle activation pattern and the experimental synergy force vectors at the test postures. Force angle deviations were the greatest for the randomly selected feasible muscle activation patterns (e.g., >100°), intermediate for effort-wise optimal muscle activation patterns (e.g., ~20°), and smallest for generalizable muscle activation patterns (e.g., <5°). Generalizable muscle activation patterns were suboptimal in terms of effort, often exceeding 50% of the maximum possible effort (cf. ~5% in minimum-effort muscle activation patterns). The feasible muscle activation ranges of individual muscles associated with producing a specific synergy force vector was reduced by ~45% when generalizability requirements were imposed. Muscles recruited in the generalizable muscle activation patterns had less sensitive torque-producing characteristics to changes in postures. We conclude that generalization of function across postures does not arise from limb biomechanics or a single optimality criterion. Muscle synergies may reflect acquired motor solutions globally tuned for generalizability across biomechanical contexts, facilitating rapid motor adaptation.

Feasible muscle activation ranges based on inverse dynamics analyses of human walking

Although it is possible to produce the same movement using an infinite number of different muscle activation patterns owing to musculoskeletal redundancy, the degree to which observed variations in muscle activity can deviate from optimal solutions computed from biomechanical models is not known. Here, we examined the range of biomechanically permitted activation levels in individual muscles during human walking using a detailed musculoskeletal model and experimentally-measured kinetics and kinematics. Feasible muscle activation ranges define the minimum and maximum possible level of each muscle's activation that satisfy inverse dynamics joint torques assuming that all other muscles can vary their activation as needed. During walking, 73% of the muscles had feasible muscle activation ranges that were greater than 95% of the total muscle activation range over more than 95% of the gait cycle, indicating that, individually, most muscles could be fully active or fully inactive while still satisfying inverse dynamics joint torques. Moreover, the shapes of the feasible muscle activation ranges did not resemble previously-reported muscle activation patterns nor optimal solutions, i.e. static optimization and computed muscle control, that are based on the same biomechanical constraints. Our results demonstrate that joint torque requirements from standard inverse dynamics calculations are insufficient to define the activation of individual muscles during walking in healthy individuals. Identifying feasible muscle activation ranges may be an effective way to evaluate the impact of additional biomechanical and/or neural constraints on possible versus actual muscle activity in both normal and impaired movements.

Stochastic modelling of muscle recruitment during activity

Muscle forces can be selected from a space of muscle recruitment strategies that produce stable motion and variable muscle and joint forces. However, current optimization methods provide only a single muscle recruitment strategy. We modelled the spectrum of muscle recruitment strategies while walking. The equilibrium equations at the joints, muscle constraints, static optimization solutions and 15-channel electromyography (EMG) recordings for seven walking cycles were taken from earlier studies. The spectrum of muscle forces was calculated using Bayesian statistics and Markov chain Monte Carlo (MCMC) methods, whereas EMG-driven muscle forces were calculated using EMG-driven modelling. We calculated the differences between the spectrum and EMG-driven muscle force for 1-15 input EMGs, and we identified the muscle strategy that best matched the recorded EMG pattern. The best-fit strategy, static optimization solution and EMG-driven force data were compared using correlation analysis. Possible and plausible muscle forces were defined as within physiological boundaries and within EMG boundaries. Possible muscle and joint forces were calculated by constraining the muscle forces between zero and the peak muscle force. Plausible muscle forces were constrained within six selected EMG boundaries. The spectrum to EMG-driven force difference increased from 40 to 108 N for 1-15 EMG inputs. The best-fit muscle strategy better described the EMG-driven pattern (R (2) = 0.94; RMSE = 19 N) than the static optimization solution (R (2) = 0.38; RMSE = 61 N). Possible forces for 27 of 34 muscles varied between zero and the peak muscle force, inducing a peak hip force of 11.3 body-weights. Plausible muscle forces closely matched the selected EMG patterns; no effect of the EMG constraint was observed on the remaining muscle force ranges. The model can be used to study alternative muscle recruitment strategies in both physiological and pathophysiological neuromotor conditions.

A general-purpose framework to simulate musculoskeletal system of human body: using a motion tracking approach

Computation of muscle force patterns that produce specified movements of muscle-actuated dynamic models is an important and challenging problem. This problem is an undetermined one, and then a proper optimization is required to calculate muscle forces. The purpose of this paper is to develop a general model for calculating all muscle activation and force patterns in an arbitrary human body movement. For this aim, the equations of a multibody system forward dynamics, which is considered for skeletal system of the human body model, is derived using Lagrange-Euler formulation. Next, muscle contraction dynamics is added to this model and forward dynamics of an arbitrary musculoskeletal system is obtained. For optimization purpose, the obtained model is used in computed muscle control algorithm, and a closed-loop system for tracking desired motions is derived. Finally, a popular sport exercise, biceps curl, is simulated by using this algorithm and the validity of the obtained results is evaluated via EMG signals.