1
|
Heras-Sádaba A, Pérez-Ruiz A, Martins P, Ederra C, de Solórzano CO, Abizanda G, Pons-Villanueva J, Calvo B, Grasa J. Exploring the muscle architecture effect on the mechanical behaviour of mouse rotator cuff muscles. Comput Biol Med 2024; 174:108401. [PMID: 38603897 DOI: 10.1016/j.compbiomed.2024.108401] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/06/2023] [Revised: 02/15/2024] [Accepted: 04/01/2024] [Indexed: 04/13/2024]
Abstract
Incorporating detailed muscle architecture aspects into computational models can enable researchers to gain deeper insights into the complexity of muscle function, movement, and performance. In this study, we employed histological, multiphoton image processing, and finite element method techniques to characterise the mechanical dependency on the architectural behaviour of supraspinatus and infraspinatus mouse muscles. While mechanical tests revealed a stiffer passive behaviour in the supraspinatus muscle, the collagen content was found to be two times higher in the infraspinatus. This effect was unveiled by analysing the alignment of fibres during muscle stretch with the 3D models and the parameters obtained in the fitting. Therefore, a strong dependence of muscle behaviour, both active and passive, was found on fibre orientation rather than collagen content.
Collapse
Affiliation(s)
- A Heras-Sádaba
- Aragón Institute of Engineering Research (i3A), Universidad de Zaragoza, Spain
| | - A Pérez-Ruiz
- Technological Innovation Division, Foundation for Applied Medical Research (FIMA), University of Navarra (UNAV), Spain; Instituto de Investigación Sanitaria de Navarra (IdiSNA), Pamplona, Spain
| | - P Martins
- Aragón Institute of Engineering Research (i3A), Universidad de Zaragoza, Spain
| | - C Ederra
- Technological Innovation Division, Foundation for Applied Medical Research (FIMA), University of Navarra (UNAV), Spain; Instituto de Investigación Sanitaria de Navarra (IdiSNA), Pamplona, Spain
| | - C Ortiz de Solórzano
- Technological Innovation Division, Foundation for Applied Medical Research (FIMA), University of Navarra (UNAV), Spain; Instituto de Investigación Sanitaria de Navarra (IdiSNA), Pamplona, Spain
| | - G Abizanda
- Technological Innovation Division, Foundation for Applied Medical Research (FIMA), University of Navarra (UNAV), Spain; Instituto de Investigación Sanitaria de Navarra (IdiSNA), Pamplona, Spain
| | - J Pons-Villanueva
- Technological Innovation Division, Foundation for Applied Medical Research (FIMA), University of Navarra (UNAV), Spain; Instituto de Investigación Sanitaria de Navarra (IdiSNA), Pamplona, Spain; Orthopedic Department, Clínica Universidad de Navarra, Pamplona, Spain
| | - B Calvo
- Aragón Institute of Engineering Research (i3A), Universidad de Zaragoza, Spain; Centro de Investigación Biomédica en Red en Bioingeniería, Biomateriales y Nanomedicina (CIBER-BBN), Spain
| | - J Grasa
- Aragón Institute of Engineering Research (i3A), Universidad de Zaragoza, Spain; Centro de Investigación Biomédica en Red en Bioingeniería, Biomateriales y Nanomedicina (CIBER-BBN), Spain.
| |
Collapse
|
2
|
Multiscale agent-based modeling of restenosis after percutaneous transluminal angioplasty: Effects of tissue damage and hemodynamics on cellular activity. Comput Biol Med 2022; 147:105753. [DOI: 10.1016/j.compbiomed.2022.105753] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/14/2021] [Revised: 04/13/2022] [Accepted: 05/13/2022] [Indexed: 11/17/2022]
|
3
|
Parikh J, Rumbell T, Butova X, Myachina T, Acero JC, Khamzin S, Solovyova O, Kozloski J, Khokhlova A, Gurev V. Generative adversarial networks for construction of virtual populations of mechanistic models: simulations to study Omecamtiv Mecarbil action. J Pharmacokinet Pharmacodyn 2021; 49:51-64. [PMID: 34716531 PMCID: PMC8837558 DOI: 10.1007/s10928-021-09787-4] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/21/2021] [Accepted: 09/23/2021] [Indexed: 11/30/2022]
Abstract
Biophysical models are increasingly used to gain mechanistic insights by fitting and reproducing experimental and clinical data. The inherent variability in the recorded datasets, however, presents a key challenge. In this study, we present a novel approach, which integrates mechanistic modeling and machine learning to analyze in vitro cardiac mechanics data and solve the inverse problem of model parameter inference. We designed a novel generative adversarial network (GAN) and employed it to construct virtual populations of cardiac ventricular myocyte models in order to study the action of Omecamtiv Mecarbil (OM), a positive cardiac inotrope. Populations of models were calibrated from mechanically unloaded myocyte shortening recordings obtained in experiments on rat myocytes in the presence and absence of OM. The GAN was able to infer model parameters while incorporating prior information about which model parameters OM targets. The generated populations of models reproduced variations in myocyte contraction recorded during in vitro experiments and provided improved understanding of OM’s mechanism of action. Inverse mapping of the experimental data using our approach suggests a novel action of OM, whereby it modifies interactions between myosin and tropomyosin proteins. To validate our approach, the inferred model parameters were used to replicate other in vitro experimental protocols, such as skinned preparations demonstrating an increase in calcium sensitivity and a decrease in the Hill coefficient of the force–calcium (F–Ca) curve under OM action. Our approach thereby facilitated the identification of the mechanistic underpinnings of experimental observations and the exploration of different hypotheses regarding variability in this complex biological system.
Collapse
Affiliation(s)
| | | | - Xenia Butova
- Institute of Immunology and Physiology, Ural Branch of the Russian Academy of Sciences (UB RAS), Yekaterinburg, Russia
| | - Tatiana Myachina
- Institute of Immunology and Physiology, Ural Branch of the Russian Academy of Sciences (UB RAS), Yekaterinburg, Russia
| | - Jorge Corral Acero
- Department of Engineering Science, Institute of Biomedical Engineering, University of Oxford, Oxford, UK
| | - Svyatoslav Khamzin
- Institute of Immunology and Physiology, Ural Branch of the Russian Academy of Sciences (UB RAS), Yekaterinburg, Russia
| | - Olga Solovyova
- Ural Federal University, Yekaterinburg, Russia.,Institute of Immunology and Physiology, Ural Branch of the Russian Academy of Sciences (UB RAS), Yekaterinburg, Russia
| | | | - Anastasia Khokhlova
- Ural Federal University, Yekaterinburg, Russia.,Institute of Immunology and Physiology, Ural Branch of the Russian Academy of Sciences (UB RAS), Yekaterinburg, Russia
| | | |
Collapse
|
4
|
Lafuente-Gracia L, Borgiani E, Nasello G, Geris L. Towards in silico Models of the Inflammatory Response in Bone Fracture Healing. Front Bioeng Biotechnol 2021; 9:703725. [PMID: 34660547 PMCID: PMC8514728 DOI: 10.3389/fbioe.2021.703725] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/30/2021] [Accepted: 09/07/2021] [Indexed: 12/21/2022] Open
Abstract
In silico modeling is a powerful strategy to investigate the biological events occurring at tissue, cellular and subcellular level during bone fracture healing. However, most current models do not consider the impact of the inflammatory response on the later stages of bone repair. Indeed, as initiator of the healing process, this early phase can alter the regenerative outcome: if the inflammatory response is too strongly down- or upregulated, the fracture can result in a non-union. This review covers the fundamental information on fracture healing, in silico modeling and experimental validation. It starts with a description of the biology of fracture healing, paying particular attention to the inflammatory phase and its cellular and subcellular components. We then discuss the current state-of-the-art regarding in silico models of the immune response in different tissues as well as the bone regeneration process at the later stages of fracture healing. Combining the aforementioned biological and computational state-of-the-art, continuous, discrete and hybrid modeling technologies are discussed in light of their suitability to capture adequately the multiscale course of the inflammatory phase and its overall role in the healing outcome. Both in the establishment of models as in their validation step, experimental data is required. Hence, this review provides an overview of the different in vitro and in vivo set-ups that can be used to quantify cell- and tissue-scale properties and provide necessary input for model credibility assessment. In conclusion, this review aims to provide hands-on guidance for scientists interested in building in silico models as an additional tool to investigate the critical role of the inflammatory phase in bone regeneration.
Collapse
Affiliation(s)
- Laura Lafuente-Gracia
- Biomechanics Section, Department of Mechanical Engineering, KU Leuven, Leuven, Belgium.,Prometheus: Division of Skeletal Tissue Engineering, KU Leuven, Leuven, Belgium
| | - Edoardo Borgiani
- Biomechanics Section, Department of Mechanical Engineering, KU Leuven, Leuven, Belgium.,Prometheus: Division of Skeletal Tissue Engineering, KU Leuven, Leuven, Belgium.,Biomechanics Research Unit, GIGA in silico Medicine, University of Liège, Liège, Belgium
| | - Gabriele Nasello
- Biomechanics Section, Department of Mechanical Engineering, KU Leuven, Leuven, Belgium.,Prometheus: Division of Skeletal Tissue Engineering, KU Leuven, Leuven, Belgium.,Skeletal Biology and Engineering Research Center, KU Leuven, Leuven, Belgium
| | - Liesbet Geris
- Biomechanics Section, Department of Mechanical Engineering, KU Leuven, Leuven, Belgium.,Prometheus: Division of Skeletal Tissue Engineering, KU Leuven, Leuven, Belgium.,Biomechanics Research Unit, GIGA in silico Medicine, University of Liège, Liège, Belgium.,Skeletal Biology and Engineering Research Center, KU Leuven, Leuven, Belgium
| |
Collapse
|
5
|
Trube N, Riedel W, Boljen M. How muscle stiffness affects human body model behavior. Biomed Eng Online 2021; 20:53. [PMID: 34078371 PMCID: PMC8170985 DOI: 10.1186/s12938-021-00876-6] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/09/2020] [Accepted: 04/04/2021] [Indexed: 11/10/2022] Open
Abstract
BACKGROUND Active human body models (AHBM) consider musculoskeletal movement and joint stiffness via active muscle truss elements in the finite element (FE) codes in dynamic application. In the latest models, such as THUMS™ Version 5, nearly all human muscle groups are modeled in form of one-dimensional truss elements connecting each joint. While a lot of work has been done to improve the active and passive behavior of this 1D muscle system in the past, the volumetric muscle system of THUMS was modeled in a much more simplified way based on Post Mortem Human Subject (PMHS) test data. The stiffness changing effect of isometric contraction was hardly considered for the volumetric muscle system of whole human body models so far. While previous works considered this aspect for single muscles, the effect of a change in stiffness due to isometric contraction of volumetric muscles on the AHBM behavior and computation time is yet unknown. METHODS In this study, a simplified frontal impact using the THUMS Version 5 AM50 occupant model was simulated. Key parameters to regulate muscle tissue stiffness of solid elements in THUMS were identified for the material model MAT_SIMPLIFIED_FOAM and different stiffness states were predefined for the buttock and thigh. RESULTS During frontal crash, changes in muscle stiffness had an effect on the overall AHBM behavior including expected injury outcome. Changes in muscle stiffness for the thigh and pelvis, as well as for the entire human body model and for strain-rate-dependent stiffness definitions based on literature data had no significant effect on the computation time. DISCUSSION Kinematics, peak impact force and stiffness changes were in general compliance with the literature data. However, different experimental setups had to be considered for comparison, as this topic has not been fully investigated experimentally in automotive applications in the past. Therefore, this study has limitations regarding validation of the frontal impact results. CONCLUSION Variations of default THUMS material model parameters allow an efficient change in stiffness of volumetric muscles for whole AHBM applications. The computation time is unaffected by altering muscle stiffness using the method suggested in this work. Due to a lack of validation data, the results of this work can only be validated with certain limitations. In future works, the default material models of THUMS could be replaced with recently published models to achieve a possibly more biofidelic muscle behavior, which would even allow a functional dependency of the 1D and 3D muscle systems. However, the effect on calculation time and model stability of these models is yet unknown and should be considered in future studies for efficient AHBM applications.
Collapse
Affiliation(s)
- Niclas Trube
- Fraunhofer-Institute for High-Speed Dynamics, Ernst-Mach-Institut, EMI, Ernst-Zermelo-Straße 4, 79104, Freiburg, Germany.
| | - Werner Riedel
- Fraunhofer-Institute for High-Speed Dynamics, Ernst-Mach-Institut, EMI, Ernst-Zermelo-Straße 4, 79104, Freiburg, Germany
| | - Matthias Boljen
- Fraunhofer-Institute for High-Speed Dynamics, Ernst-Mach-Institut, EMI, Ernst-Zermelo-Straße 4, 79104, Freiburg, Germany
| |
Collapse
|
6
|
Characterization of hyperelastic mechanical properties for youth corneal anterior central stroma based on collagen fibril crimping constitutive model. J Mech Behav Biomed Mater 2020; 103:103575. [DOI: 10.1016/j.jmbbm.2019.103575] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/17/2019] [Revised: 11/03/2019] [Accepted: 11/29/2019] [Indexed: 11/19/2022]
|
7
|
Lawson BAJ, Drovandi CC, Cusimano N, Burrage P, Rodriguez B, Burrage K. Unlocking data sets by calibrating populations of models to data density: A study in atrial electrophysiology. SCIENCE ADVANCES 2018; 4:e1701676. [PMID: 29349296 PMCID: PMC5770172 DOI: 10.1126/sciadv.1701676] [Citation(s) in RCA: 47] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/19/2017] [Accepted: 12/08/2017] [Indexed: 05/08/2023]
Abstract
The understanding of complex physical or biological systems nearly always requires a characterization of the variability that underpins these processes. In addition, the data used to calibrate these models may also often exhibit considerable variability. A recent approach to deal with these issues has been to calibrate populations of models (POMs), multiple copies of a single mathematical model but with different parameter values, in response to experimental data. To date, this calibration has been largely limited to selecting models that produce outputs that fall within the ranges of the data set, ignoring any trends that might be present in the data. We present here a novel and general methodology for calibrating POMs to the distributions of a set of measured values in a data set. We demonstrate our technique using a data set from a cardiac electrophysiology study based on the differences in atrial action potential readings between patients exhibiting sinus rhythm (SR) or chronic atrial fibrillation (cAF) and the Courtemanche-Ramirez-Nattel model for human atrial action potentials. Not only does our approach accurately capture the variability inherent in the experimental population, but we also demonstrate how the POMs that it produces may be used to extract additional information from the data used for calibration, including improved identification of the differences underlying stratified data. We also show how our approach allows different hypotheses regarding the variability in complex systems to be quantitatively compared.
Collapse
Affiliation(s)
- Brodie A. J. Lawson
- Australian Research Council Centre of Excellence for Mathematical and Statistical Frontiers, School of Mathematical Sciences, Queensland University of Technology, Brisbane, Queensland, Australia
- Corresponding author.
| | - Christopher C. Drovandi
- Australian Research Council Centre of Excellence for Mathematical and Statistical Frontiers, School of Mathematical Sciences, Queensland University of Technology, Brisbane, Queensland, Australia
| | | | - Pamela Burrage
- Australian Research Council Centre of Excellence for Mathematical and Statistical Frontiers, School of Mathematical Sciences, Queensland University of Technology, Brisbane, Queensland, Australia
| | - Blanca Rodriguez
- Department of Computer Science, University of Oxford, Oxford, UK
| | - Kevin Burrage
- Australian Research Council Centre of Excellence for Mathematical and Statistical Frontiers, School of Mathematical Sciences, Queensland University of Technology, Brisbane, Queensland, Australia
- Department of Computer Science, University of Oxford, Oxford, UK
| |
Collapse
|
8
|
Grasa J, Sierra M, Lauzeral N, Muñoz M, Miana-Mena F, Calvo B. Active behavior of abdominal wall muscles: Experimental results and numerical model formulation. J Mech Behav Biomed Mater 2016; 61:444-454. [DOI: 10.1016/j.jmbbm.2016.04.013] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/07/2016] [Revised: 03/29/2016] [Accepted: 04/06/2016] [Indexed: 10/22/2022]
|
9
|
Muszkiewicz A, Britton OJ, Gemmell P, Passini E, Sánchez C, Zhou X, Carusi A, Quinn TA, Burrage K, Bueno-Orovio A, Rodriguez B. Variability in cardiac electrophysiology: Using experimentally-calibrated populations of models to move beyond the single virtual physiological human paradigm. PROGRESS IN BIOPHYSICS AND MOLECULAR BIOLOGY 2015; 120:115-27. [PMID: 26701222 PMCID: PMC4821179 DOI: 10.1016/j.pbiomolbio.2015.12.002] [Citation(s) in RCA: 104] [Impact Index Per Article: 11.6] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 09/16/2015] [Revised: 11/24/2015] [Accepted: 12/02/2015] [Indexed: 01/13/2023]
Abstract
Physiological variability manifests itself via differences in physiological function between individuals of the same species, and has crucial implications in disease progression and treatment. Despite its importance, physiological variability has traditionally been ignored in experimental and computational investigations due to averaging over samples from multiple individuals. Recently, modelling frameworks have been devised for studying mechanisms underlying physiological variability in cardiac electrophysiology and pro-arrhythmic risk under a variety of conditions and for several animal species as well as human. One such methodology exploits populations of cardiac cell models constrained with experimental data, or experimentally-calibrated populations of models. In this review, we outline the considerations behind constructing an experimentally-calibrated population of models and review the studies that have employed this approach to investigate variability in cardiac electrophysiology in physiological and pathological conditions, as well as under drug action. We also describe the methodology and compare it with alternative approaches for studying variability in cardiac electrophysiology, including cell-specific modelling approaches, sensitivity-analysis based methods, and populations-of-models frameworks that do not consider the experimental calibration step. We conclude with an outlook for the future, predicting the potential of new methodologies for patient-specific modelling extending beyond the single virtual physiological human paradigm.
Collapse
Affiliation(s)
- Anna Muszkiewicz
- Department of Computer Science, University of Oxford, Parks Road, Oxford OX1 3QD, United Kingdom
| | - Oliver J Britton
- Department of Computer Science, University of Oxford, Parks Road, Oxford OX1 3QD, United Kingdom
| | - Philip Gemmell
- Clyde Biosciences Ltd, West Medical Building, University of Glasgow, Glasgow G12 8QQ, United Kingdom
| | - Elisa Passini
- Department of Computer Science, University of Oxford, Parks Road, Oxford OX1 3QD, United Kingdom
| | - Carlos Sánchez
- Center for Computational Medicine in Cardiology (CCMC), Institute of Computational Science, Università della Svizzera italiana, Lugano, Switzerland
| | - Xin Zhou
- Department of Computer Science, University of Oxford, Parks Road, Oxford OX1 3QD, United Kingdom
| | | | - T Alexander Quinn
- Department of Physiology and Biophysics, Dalhousie University, Halifax, Nova Scotia, Canada
| | - Kevin Burrage
- Department of Computer Science, University of Oxford, Parks Road, Oxford OX1 3QD, United Kingdom; Mathematical Sciences, Queensland University of Technology, Queensland 4072, Australia; ACEMS, ARC Centre of Excellence for Mathematical and Statistical Frontiers, Queensland University of Technology, Queensland 4072, Australia
| | - Alfonso Bueno-Orovio
- Department of Computer Science, University of Oxford, Parks Road, Oxford OX1 3QD, United Kingdom
| | - Blanca Rodriguez
- Department of Computer Science, University of Oxford, Parks Road, Oxford OX1 3QD, United Kingdom.
| |
Collapse
|