Dataset Information

Neuro-musculoskeletal flexible multibody simulation yields a framework for efficient bone failure risk assessment.

ABSTRACT: Fragility fractures are a major socioeconomic problem. A non-invasive, computationally-efficient method for the identification of fracture risk scenarios under the representation of neuro-musculoskeletal dynamics does not exist. We introduce a computational workflow that integrates modally-reduced, quantitative CT-based finite-element models into neuro-musculoskeletal flexible multibody simulation (NfMBS) for early bone fracture risk assessment. Our workflow quantifies the bone strength via the osteogenic stresses and strains that arise due to the physiological-like loading of the bone under the representation of patient-specific neuro-musculoskeletal dynamics. This allows for non-invasive, computationally-efficient dynamic analysis over the enormous parameter space of fracture risk scenarios, while requiring only sparse clinical data. Experimental validation on a fresh human femur specimen together with femur strength computations that were consistent with literature findings provide confidence in the workflow: The simulation of an entire squat took only 38?s CPU-time. Owing to the loss (16% cortical, 33% trabecular) of bone mineral density (BMD), the strain measure that is associated with bone fracture increased by 31.4%; and yielded an elevated risk of a femoral hip fracture. Our novel workflow could offer clinicians with decision-making guidance by enabling the first combined in-silico analysis tool using NfMBS and BMD measurements for optimized bone fracture risk assessment.

SUBMITTER: Geier A

PROVIDER: S-EPMC6503141 | biostudies-literature | 2019 May

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Neuro-musculoskeletal flexible multibody simulation yields a framework for efficient bone failure risk assessment.

Geier Andreas A Kebbach Maeruan M Soodmand Ehsan E Woernle Christoph C Kluess Daniel D Bader Rainer R

Scientific reports 20190506 1

Fragility fractures are a major socioeconomic problem. A non-invasive, computationally-efficient method for the identification of fracture risk scenarios under the representation of neuro-musculoskeletal dynamics does not exist. We introduce a computational workflow that integrates modally-reduced, quantitative CT-based finite-element models into neuro-musculoskeletal flexible multibody simulation (NfMBS) for early bone fracture risk assessment. Our workflow quantifies the bone strength via the ...[more]

PMID: 31061388

Similar Datasets

Project description:Optimal control simulations have shown that both musculoskeletal dynamics and physiological noise are important determinants of movement. However, due to the limited efficiency of available computational tools, deterministic simulations of movement focus on accurately modelling the musculoskeletal system while neglecting physiological noise, and stochastic simulations account for noise while simplifying the dynamics. We took advantage of recent approaches where stochastic optimal control problems are approximated using deterministic optimal control problems, which can be solved efficiently using direct collocation. We were thus able to extend predictions of stochastic optimal control as a theory of motor coordination to include muscle coordination and movement patterns emerging from non-linear musculoskeletal dynamics. In stochastic optimal control simulations of human standing balance, we demonstrated that the inclusion of muscle dynamics can predict muscle co-contraction as minimal effort strategy that complements sensorimotor feedback control in the presence of sensory noise. In simulations of reaching, we demonstrated that nonlinear multi-segment musculoskeletal dynamics enables complex perturbed and unperturbed reach trajectories under a variety of task conditions to be predicted. In both behaviors, we demonstrated how interactions between task constraint, sensory noise, and the intrinsic properties of muscle influence optimal muscle coordination patterns, including muscle co-contraction, and the resulting movement trajectories. Our approach enables a true minimum effort solution to be identified as task constraints, such as movement accuracy, can be explicitly imposed, rather than being approximated using penalty terms in the cost function. Our approximate stochastic optimal control framework predicts complex features, not captured by previous simulation approaches, providing a generalizable and valuable tool to study how musculoskeletal dynamics and physiological noise may alter neural control of movement in both healthy and pathological movements.

Project description:BackgroundSpasticity is an important complication after stroke, especially in the anti-gravity muscles, i.e. lower limb extensors. However the contribution of hyperexcitable muscle spindle reflex loops to gait impairments after stroke is often disputed. In this study a neuro-musculoskeletal model was developed to investigate the contribution of an increased length and velocity feedback and altered reflex modulation patterns to hemiparetic gait deficits.MethodsA musculoskeletal model was extended with a muscle spindle model providing real-time length and velocity feedback of gastrocnemius, soleus, vasti and rectus femoris during a forward dynamic simulation (neural control model). By using a healthy subject's base muscle excitations, in combination with increased feedback gains and altered reflex modulation patterns, the effect on kinematics was simulated. A foot-ground contact model was added to account for the interaction effect between the changed kinematics and the ground. The qualitative effect i.e. the directional effect and the specific gait phases where the effect is present, on the joint kinematics was then compared with hemiparetic gait deviations reported in the literature.ResultsOur results show that increased feedback in combination with altered reflex modulation patterns of soleus, vasti and rectus femoris muscle can contribute to excessive ankle plantarflexion/inadequate dorsiflexion, knee hyperextension/inadequate flexion and increased hip extension/inadequate flexion during dedicated gait cycle phases. Increased feedback of gastrocnemius can also contribute to excessive plantarflexion/inadequate dorsiflexion, however in combination with excessive knee and hip flexion. Increased length/velocity feedback can therefore contribute to two types of gait deviations, which are both in accordance with previously reported gait deviations in hemiparetic patients. Furthermore altered modulation patterns, in particular the reduced suppression of the muscle spindle feedback during swing, can contribute largely to an increased plantarflexion and knee extension during the swing phase and consequently to hampered toe clearance.ConclusionsOur results support the idea that hyperexcitability of length and velocity feedback pathways, especially in combination with altered reflex modulation patterns, can contribute to deviations in hemiparetic gait. Surprisingly, our results showed only subtle temporal differences between length and velocity feedback. Therefore, we cannot attribute the effects seen in kinematics to one specific type of feedback.

Project description:Musculoskeletal geometry and muscle volumes vary widely in the population and are intricately linked to the performance of tasks ranging from walking and running to jumping and sprinting. As an alternative to experimental approaches, where it is difficult to isolate factors and establish causal relationships, simulations can be used to independently vary musculoskeletal geometry and muscle volumes, and develop a fundamental understanding. However, our ability to understand how these parameters affect task performance has been limited due to the high computational cost of modelling the necessary complexity of the musculoskeletal system and solving the requisite multi-dimensional optimization problem. For example, sprinting and running are fundamental to many forms of sport, but past research on the relationships between musculoskeletal geometry, muscle volumes, and running performance has been limited to observational studies, which have not established cause-effect relationships, and simulation studies with simplified representations of musculoskeletal geometry. In this study, we developed a novel musculoskeletal simulator that is differentiable with respect to musculoskeletal geometry and muscle volumes. This simulator enabled us to find the optimal body segment dimensions and optimal distribution of added muscle volume for sprinting and marathon running. Our simulation results replicate experimental observations, such as increased muscle mass in sprinters, as well as a mass in the lower end of the healthy BMI range and a higher leg-length-to-height ratio in marathon runners. The simulations also reveal new relationships, for example showing that hip musculature is vital to both sprinting and marathon running. We found hip flexor and extensor moment arms were maximized to optimize sprint and marathon running performance, and hip muscles the main target when we simulated strength training for sprinters. Our simulation results provide insight to inspire future studies to examine optimal strength training. Our simulator can be extended to other athletic tasks, such as jumping, or to non-athletic applications, such as designing interventions to improve mobility in older adults or individuals with movement disorders.

Project description:PurposeEpistasis, the interaction between two or more genes, is integral to the study of genetics and is present throughout nature. Yet, it is seldom fully explored as most approaches primarily focus on single-locus effects, partly because analyzing all pairwise and higher-order interactions requires significant computational resources. Furthermore, existing methods for epistasis detection only consider a Cartesian (multiplicative) model for interaction terms. This is likely limiting as epistatic interactions can evolve to produce varied relationships between genetic loci, some complex and not linearly separable.MethodsWe present new algorithms for the interaction coefficients for standard regression models for epistasis that permit many varied models for the interaction terms for loci and efficient memory usage. The algorithms are given for two-way and three-way epistasis and may be generalized to higher order epistasis. Statistical tests for the interaction coefficients are also provided. We also present an efficient matrix based algorithm for permutation testing for two-way epistasis. We offer a proof and experimental evidence that methods that look for epistasis only at loci that have main effects may not be justified. Given the computational efficiency of the algorithm, we applied the method to a rat data set and mouse data set, with at least 10,000 loci and 1,000 samples each, using the standard Cartesian model and the XOR model to explore body mass index.ResultsThis study reveals that although many of the loci found to exhibit significant statistical epistasis overlap between models in rats, the pairs are mostly distinct. Further, the XOR model found greater evidence for statistical epistasis in many more pairs of loci in both data sets with almost all significant epistasis in mice identified using XOR. In the rat data set, loci involved in epistasis under the XOR model are enriched for biologically relevant pathways.ConclusionOur results in both species show that many biologically relevant epistatic relationships would have been undetected if only one interaction model was applied, providing evidence that varied interaction models should be implemented to explore epistatic interactions that occur in living systems.

Dataset Information

Neuro-musculoskeletal flexible multibody simulation yields a framework for efficient bone failure risk assessment.

Publications

Neuro-musculoskeletal flexible multibody simulation yields a framework for efficient bone failure risk assessment.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets