A new framework for computational biomechanical models and 3Rs in musculoskeletal research.

Lead Research Organisation: University of Leeds
Department Name: Sch of Biomedical Sciences

Abstract

Abstracts are not currently available in GtR for all funded research. This is normally because the abstract was not required at the time of proposal submission, but may be because it included sensitive information such as personal details.

Technical Summary

This project has two primary goals: (1) to quantify the experimental input data required for musculoskeletal computer models to deliver accurate biomechanical predictions, and (2) provide quantitative data on the current and future potential of computer simulation approaches as an alternative to animal experimentation in basic science, clinical/veterinary and industrial projects involving biomechanics. To achieve this we will use mastication in rabbits as our model system because it is experimentally tractable and will generate data relevant to basic biological science fields and the health and welfare of rabbits. We will collect a wealth of anatomical, in vivo and in vitro experimental data on the rabbit masticatory system. Biplanar x-ray videography will be synchronized with strain gauges, muscle EMG and sonomicrometry and a bite force transducer to simultaneously record 3D motions, muscle dynamics, bone strains and bite forces during feeding. In vitro physiology experiments will quantify muscle mechanics and key contractile properties. Key anatomical variables from the same rabbits will be quantified using dissection, medical imaging and nanoindentation. Initially this data will be used to test hypotheses relating functional adaptations to processing different food types and to understand how tissue level properties translate to overall skull function. But most importantly this data will allow us to build, drive and validate the most comprehensive musculoskeletal computer models produced to-date. Our initial models will thus provide a best-case scenario for simulation accuracy. We will then incrementally reduce or average-out the resolution of input data given to the model and observe the effects on accuracy. This will tell us how individual input parameters affect accuracy and help the musculoskeletal research community identify which parameters do not need to be measured through experimentation in real animals to achieve specific thresholds of accuracy.

Planned Impact

The main impacts of this research programme are on animal welfare, the general public, and the researchers employed on the grant, together with benefits to the academic community (see Academic Beneficiaries).

APPLIED LINKS WITH THE POTENTIAL TO DELIVER ANIMAL WELFARE BENEFITS THROUGH THE 3Rs
Publication and dissemination of computational models will benefit a very wide range of researchers, potentially leading to significant societal impact in terms of improved animal welfare. Through the development of the most robust, validated computer models ever constructed we will facilitate a step change in our understanding of model accuracy and usefulness, and in doing so provide a marker for future musculoskeletal research. Our models will demonstrate how precise (or how species-specific) input parameters have to be to achieve different thresholds of accuracy in terms of predicting an enormous range of muscle and bone parameters. This will directly guide other researchers utilising musculoskeletal models in a variety of contexts, ranging from basic biologists (e.g. functional anatomists, palaeontologists), to applied researchers (clinicians/veterinarians) interested in understanding healthy and "abnormal" biomechanics, through to bioengineers and roboticists developing medical devices and bioinspired robots.

As a direct consequence of this we believe that our models will demonstrate that computer simulation approaches can contribute significantly to achieving the 3Rs in musculoskeletal research, and in surgical and implant design studies that require biomechanical analyses. Experimental studies of this type are highly invasive and typically cause significant pain and distress to the animals before they are euthanized. The approach we will take herein will mean that the potential for future benefits from 3Rs in this context will be broad. By demonstrating how individual anatomical and physiological parameters (e.g. muscle size, contraction behaviour) impact on model accuracy our data will indicate which variables need to be measured directly from cohorts or species of animal in the future and which do not, given the threshold of accuracy for specific model parameters required. In the future, the application of such models will allow some animal experiments to be replaced and in other cases refined and/or reduced as model simulations will allow research efforts involving animal research to be better designed.

THE GENERAL PUBLIC
Animal anatomy and functional morphology are topics that interest the general public, as testified by the continued success of natural history documentaries and recent documentaries like "Inside Nature's Giants." Medical imaging and computational models are powerful and flexible mediums through which to communicate our understanding of the links between anatomy, mechanics and ecology to the general public. Our work will have a positive impact in this respect and demonstrating technological advances and a wide variety of applications of biological research and will have an impact on the general public as a whole, but in particular school children.

OTHER SPECIFIC IMPACTS
The careers of the PDRAs will be developed considerably by involvement in a multi-disciplinary project that integrates physiology, biomechanics and modelling approaches. They will benefit from working closely with laboratories in different leading institutions. The research will also impact on the training of undergraduates carrying out research projects and postgraduates completing MSc/PhD's within our labs.
 
Description We collected the most exhaustive anatomical, mechanical and physiological data set on the masticatory system of a single cohort of animals from any species to-date. This data set has and will continue to contribute to publications about the masticatory anatomy and function of the rabbit, which is a key model organism on comparative academic research and the bone implant industry. For example, we have shown that the architecture and function of the digastric is mechanistically linked to its unusually low muscle spindle abundance, supporting the hypothesis that muscles with a low spindle abundance appear to primarily operate as energy absorbers or brakes. In other work, we have shown that rabbits dynamically modulate the temporal kinetics of feeding and recruitment of muscles in order effectively process foods of differing material properties.

This data set has also provided a unique basis from which to construct computational models using two separate against which to judge the accuracy of two computer modelling techniques (multi-body dynamics analysis and finite element analysis), and subsequently judge their accuracy and their potential to contribute to the replacement, reduction and refinement (3Rs) of animal use in future musculoskeletal research. At present wWe have fully evaluated the predictive accuracy of a subject-specific multi-body dynamics (MDA) model during chewing on different food types in the in-vivo experiments, and during stimulated maximal incisor and molar biting, by comparing model predictions to measured data from the same rabbit. We have also nearly completed analysedis the effect of where we simplifying the model complexity (e.g. using different methods of muscle modelling) and remove or averaginge out input data (e.g. replacing subject-specific data with either averaged data from our experimental rabbits and/or literature data) into this model (replacing subject-specific data with either averaged data from our experimental rabbits and/or literature data) to assess the level of data required for specific accuracy thresholds. This will inform future studies and potentially allow higher levels of 3Rs in the construction of subject-specific models. In parallel to this, we also evaluateding the accuracy with which these different model iterations match species-averaged data, thereby evaluating the extent to which generic or species-averaged models might enabled greater 3Rs in future studies.

We have measured and analysed the variation of material properties of bone in the rabbit skull. We found the bone properties vary not only in terms of the location in the skull at which they are measured (i.e. front of the skull vs the back), but also in terms of the direction in which they are measured (i.e. they can be measured in one of three axes). For example, the bone properties were found to be lower in one axis, when compared to the other two axes. This will provide a valuable data set that can inform future computational models of the rabbit.

We have constructed a high-resolution finite element (FE) model of the same rabbit used to build the MDA model above. We have are in the process of simulateding chewing cycles and evaluateding the accuracy of this subject-specific model by combining simulated strains in the cranium to those measured during in-vivo chewing experiments and during stimulated maximal incisor and molar biting. We have also in the process of investigated ing the influence of the measured bone properties on the outputs of the FEA model. This has material variation on the FEA outputs. This will be used to provided suggestions for the determine the quantity of the material property input data required (e.g. isotropic vs orthotropic) homogeneous vs subject-specific) to obtained specific levels in accuracy of the outputs and thus provide a quantitative basis for 3R-related evaluation in future experimental design. For example, we have shown that isotropic material properties are suitable when investigating general strain distribution, and it best avoid using material properties measured in certain directions.
Exploitation Route Once published and made freely available, oOur "off the self" computational models can be modified to replicate a range industrial biomedical experiments that observe the cause-effect relationship between the musculoskeletal system and biomechanical loading. The models are directly applicable to experiments that reply on cranial biomechanics as the anatomical system (e.g. dental implant design, biomaterial osseointegration/resorption, effects of tooth loosening and loss, effects of jaw joint disorders). However, Tthe models will be easily modifiable in order to predict muscle/bone adaptations in other anatomical systems and/or species that are commonly used in biomedical experiments. This will supply validated in silico techniques that can initially supplement biomedical experiments, with the real long-term potential to aid industry and academic research in making a significant contribution to reduce, refine and replace animal experiments. The knowledge and understanding of model accuracy that we provide will be directly transferable to other areas of implant and prosthetic design, thus benefiting this industry very widely.

Dental disease is one of the most common welfare issues in domesticated rabbits and our experimental data set and computational models will be directly applicable to veterinary research into its causes and progression. This area is currently controversial, but diet and food processing have been identified as an important factor. In the medium-longer term both our experimental and modelling work should inspire similar approaches in other areas of veterinary and animal welfare research.

We believe that our models demonstrate that computer simulation approaches can contribute significantly to achieving the 3Rs in musculoskeletal research, and in surgical and implant design studies that require biomechanical analyses.. Experimental studies of this type are highly invasive, and typically cause significant pain and distress to the animals before they are euthanized. By demonstrating how individual anatomical and physiological parameters (e.g. muscle size, architecture, contraction behaviour) impact on model accuracy, our data indicates which variables need to be measured directly from cohorts or species of animal in the future and which do not, given the threshold of accuracy for specific model parameters required. Thus future studies will be able to refine and reduce animal use through computer models on the basis of our data.
Sectors Healthcare

Pharmaceuticals and Medical Biotechnology

 
Title ESM Raw data from Skeletal muscle function underpins muscle spindle abundance 
Description Muscle spindle abundance is highly variable within and across species, but we currently lack any clear picture of the mechanistic causes or consequences of this variation. Previous use of spindle abundance as a correlate for muscle function implies a mechanical underpinning to this variation, but these ideas have not been tested. Herein, we use integrated medical imaging and subject-specific musculoskeletal models to investigate the relationship between spindle abundance, muscle architecture and in vivo muscle behaviour in the human locomotor system. These analyses indicate that muscle spindle number is tightly correlated with muscle fascicle length, absolute fascicle length change, velocity of fibre lengthening and active muscle forces during walking. Novel correlations between functional indices and spindle abundance are also recovered, where muscles with a high abundance predominantly function as springs, compared to those with a lower abundance mostly functioning as brakes during walking. These data demonstrate that muscle fibre length, lengthening velocity and fibre force are key physiological signals to the central nervous system and its modulation of locomotion, and that muscle spindle abundance may be tightly correlated to how a muscle generates work. These insights may be combined with neuromechanics and robotic studies of motor control to help further tease apart the functional drivers of muscle spindle composition. 
Type Of Material Database/Collection of data 
Year Produced 2022 
Provided To Others? Yes  
URL https://rs.figshare.com/articles/dataset/ESM_Raw_data_from_Skeletal_muscle_function_underpins_muscle...
 
Title ESM Raw data from Skeletal muscle function underpins muscle spindle abundance 
Description Muscle spindle abundance is highly variable within and across species, but we currently lack any clear picture of the mechanistic causes or consequences of this variation. Previous use of spindle abundance as a correlate for muscle function implies a mechanical underpinning to this variation, but these ideas have not been tested. Herein, we use integrated medical imaging and subject-specific musculoskeletal models to investigate the relationship between spindle abundance, muscle architecture and in vivo muscle behaviour in the human locomotor system. These analyses indicate that muscle spindle number is tightly correlated with muscle fascicle length, absolute fascicle length change, velocity of fibre lengthening and active muscle forces during walking. Novel correlations between functional indices and spindle abundance are also recovered, where muscles with a high abundance predominantly function as springs, compared to those with a lower abundance mostly functioning as brakes during walking. These data demonstrate that muscle fibre length, lengthening velocity and fibre force are key physiological signals to the central nervous system and its modulation of locomotion, and that muscle spindle abundance may be tightly correlated to how a muscle generates work. These insights may be combined with neuromechanics and robotic studies of motor control to help further tease apart the functional drivers of muscle spindle composition. 
Type Of Material Database/Collection of data 
Year Produced 2022 
Provided To Others? Yes  
URL https://rs.figshare.com/articles/dataset/ESM_Raw_data_from_Skeletal_muscle_function_underpins_muscle...
 
Description Collaboration with Dr James Charles 
Organisation Royal Veterinary College (RVC)
Country United Kingdom 
Sector Academic/University 
PI Contribution We have provided physiological data as inputs to a Biomechanical computational model and predictions of muscle performance during locomotion to be compared to the model.
Collaborator Contribution Dr Charles has provided data from a computational model of the predictions of muscle performance during locomotion.
Impact Ongoing collaboration
Start Year 2022
 
Description NC3Rs workshop (University of Liverpool) 
Form Of Engagement Activity A formal working group, expert panel or dialogue
Part Of Official Scheme? No
Geographic Reach Local
Primary Audience Third sector organisations
Results and Impact NC3Rs workshop/ symposium involving talks and networking - attending by University of Liverpool researchers (PIs, postdocs, PhD students). Sparked questions and discussions about our research.
Year(s) Of Engagement Activity 2020