1
|
Harrison K, Loundagin L, Hiebert B, Panahifar A, Zhu N, Marchiori D, Arnason T, Swekla K, Pivonka P, Cooper D. Glucocorticoids disrupt longitudinal advance of cortical bone basic multicellular units in the rabbit distal tibia. Bone 2024; 187:117171. [PMID: 38901788 DOI: 10.1016/j.bone.2024.117171] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 02/17/2024] [Revised: 06/13/2024] [Accepted: 06/15/2024] [Indexed: 06/22/2024]
Abstract
Glucocorticoids (GCs) are the leading cause of secondary osteoporosis. The emerging perspective, derived primarily from 2D histological study of trabecular bone, is that GC-induced bone loss arises through the uncoupling of bone formation and resorption at the level of the basic multicellular unit (BMU), which carries out bone remodeling. Here we explore the impact of GCs on cortical bone remodeling in the rabbit model. Based upon the rapid reduction of bone formation and initial elevation of resorption caused by GCs, we hypothesized that the rate of advance (longitudinal erosion rate; LER) of cortical BMUs would be increased. To test this hypothesis we divided 20 female New Zealand White rabbits into four experimental groups: ovariohysterectomy (OVH), glucocorticoid (GC), OVH + GC and SHAM controls (n = 5 animals each). Ten weeks post-surgery (OVH or sham), and two weeks after the initiation of dosing (daily subcutaneous injections of 1.5 mg/kg of methylprednisolone sodium succinate in the GC-treated groups and 1 ml of saline for the others), the right tibiae were scanned in vivo using Synchrotron Radiation (SR) in-line phase contrast micro-CT at the Canadian Light Source. After an additional 2 weeks of dosing, the rabbits were euthanized and ex vivo images were collected using desktop micro-CT. The datasets were co-registered in 3D and LER was calculated as the distance traversed by BMU cutting-cones in the 14-day interval between scans. Counter to our hypothesis, LER was greatly reduced in GC-treated rabbits. Mean LER was lower in GC (4.27 μm/d; p < 0.001) and OVH + GC (4.19 μm/d; p < 0.001), while similar in OVH (40.13 μm/d; p = 0.990), compared to SHAM (40.44 μm/d). This approximately 90 % reduction in LER with GCs was also associated with an overall disruption of BMU progression, with radial expansion of the remodeling space occurring in all directions. This unexpected outcome suggests that GCs do not simply uncouple formation and resorption within cortical BMUs and highlights the value of the time-lapsed 4D approach employed.
Collapse
Affiliation(s)
- Kim Harrison
- Department of Anatomy, Physiology and Pharmacology, College of Medicine, University of Saskatchewan, Saskatoon, Canada
| | - Lindsay Loundagin
- Department of Anatomy, Physiology and Pharmacology, College of Medicine, University of Saskatchewan, Saskatoon, Canada
| | - Beverly Hiebert
- Department of Anatomy, Physiology and Pharmacology, College of Medicine, University of Saskatchewan, Saskatoon, Canada; Max Rady College of Medicine, University of Manitoba, Winnipeg, Canada
| | - Arash Panahifar
- BioMedical Imaging and Therapy Beamline, Canadian Light Source, Saskatoon, Canada; Department of Medical Imaging, College of Medicine, University of Saskatchewan, Saskatoon, Canada
| | - Ning Zhu
- BioMedical Imaging and Therapy Beamline, Canadian Light Source, Saskatoon, Canada; Division of Biomedical Engineering, College of Engineering, University of Saskatchewan, Saskatoon, Canada
| | - Denver Marchiori
- Department of Anatomy, Physiology and Pharmacology, College of Medicine, University of Saskatchewan, Saskatoon, Canada
| | - Terra Arnason
- Medicine Dept of Endocrinology, Faculty of Medicine & Dentistry, University of Alberta, Edmonton, Canada
| | - Kurtis Swekla
- Animal Care and Research Support Office, Office of the Vice President of Research, University of Saskatchewan, Saskatoon, Canada
| | - Peter Pivonka
- School of Mechanical, Medical, and Process Engineering, Queensland University of Technology, Brisbane, Australia
| | - David Cooper
- Department of Anatomy, Physiology and Pharmacology, College of Medicine, University of Saskatchewan, Saskatoon, Canada.
| |
Collapse
|
2
|
Shi V, Morgan EF. Estrogen and estrogen receptors mediate the mechanobiology of bone disease and repair. Bone 2024:117220. [PMID: 39106937 DOI: 10.1016/j.bone.2024.117220] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 06/08/2024] [Revised: 07/28/2024] [Accepted: 07/30/2024] [Indexed: 08/09/2024]
Abstract
It is well understood that the balance of bone formation and resorption is dependent on both mechanical and biochemical factors. In addition to cell-secreted cytokines and growth factors, sex hormones like estrogen are critical to maintaining bone health. Although the direct osteoprotective function of estrogen and estrogen receptors (ERs) has been reported extensively, evidence that estrogen signaling also has a role in mediating the effects of mechanical loading on maintenance of bone mass and healing of bone injuries has more recently emerged. Recent studies have underscored the role of estrogen and ERs in many pathways of bone mechanosensation and mechanotransduction. Estrogen and ERs have been shown to augment integrin-based mechanotransduction as well as canonical Wnt/b-catenin, RhoA/ROCK, and YAP/TAZ pathways. Estrogen and ERs also influence the mechanosensitivity of not only osteocytes but also osteoblasts, osteoclasts, and marrow stromal cells. The current review will highlight these roles of estrogen and ERs in cellular mechanisms underlying bone mechanobiology and discuss their implications for management of osteoporosis and bone fractures. A greater understanding of the mechanisms behind interactions between estrogen and mechanical loading may be crucial to addressing the shortcomings of current hormonal and pharmaceutical therapies. A combined therapy approach including high-impact exercise therapy may mitigate adverse side effects and allow an effective long-term solution for the prevention, treatment, and management of bone fragility in at-risk populations. Furthermore, future implications to novel local delivery mechanisms of hormonal therapy for osteoporosis treatment, as well as the effects on bone health of applications of sex hormone therapy outside of bone disease, will be discussed.
Collapse
Affiliation(s)
- Vivian Shi
- Boston University, Department of Biomedical Engineering, 44 Cummington St, Boston 02215, MA, USA; Center for Multiscale and Translational Mechanobiology, Boston University, 44 Cummington St, Boston 02215, MA, USA
| | - Elise F Morgan
- Boston University, Department of Biomedical Engineering, 44 Cummington St, Boston 02215, MA, USA; Center for Multiscale and Translational Mechanobiology, Boston University, 44 Cummington St, Boston 02215, MA, USA.
| |
Collapse
|
3
|
Sánchez JF, Ramtani S, Boucetta A, Velasco MA, Vaca-González JJ, Duque-Daza CA, Garzón-Alvarado DA. Tumor growth for remodeling process: A 2D approach. J Theor Biol 2024; 585:111781. [PMID: 38432504 DOI: 10.1016/j.jtbi.2024.111781] [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: 11/01/2023] [Revised: 02/07/2024] [Accepted: 02/27/2024] [Indexed: 03/05/2024]
Abstract
This paper aims to present a comprehensive framework for coupling tumor-bone remodeling processes in a 2-dimensional geometry. This is achieved by introducing a bio-inspired damage that represents the growing tumor, which subsequently affects the main populations involved in the remodeling process, namely, osteoclasts, osteoblasts, and bone tissue. The model is constructed using a set of differential equations based on the Komarova's and Ayati's models, modified to incorporate the bio-inspired damage that may result in tumor mass formation. Three distinct models were developed. The first two models are based on the Komarova's governing equations, with one demonstrating an osteolytic behavior and the second one an osteoblastic model. The third model is a variation of Ayati's model, where the bio-inspired damage is induced through the paracrine and autocrine parameters, exhibiting an osteolytic behavior. The obtained results are consistent with existing literature, leading us to believe that our in-silico experiments will serve as a cornerstone for paving the way towards targeted interventions and personalized treatment strategies, ultimately improving the quality of life for those affected by these conditions.
Collapse
Affiliation(s)
| | - Salah Ramtani
- Laboratoire CSPBAT, equipe LBPS, CNRS (UMR 7244), Universit e Sorbonne Paris Nord, France.
| | - Abdelkader Boucetta
- Laboratoire CSPBAT, equipe LBPS, CNRS (UMR 7244), Universit e Sorbonne Paris Nord, France
| | | | - Juan Jairo Vaca-González
- Escuela de Pregrado - Direccion Académica, Universidad Nacional de Colombia, Sede de La Paz, Colombia.
| | | | | |
Collapse
|
4
|
Yang L, Chen H, Yang C, Hu Z, Jiang Z, Meng S, Liu R, Huang L, Yang K. Research progress on the regulatory mechanism of integrin-mediated mechanical stress in cells involved in bone metabolism. J Cell Mol Med 2024; 28:e18183. [PMID: 38506078 PMCID: PMC10951882 DOI: 10.1111/jcmm.18183] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/15/2023] [Revised: 01/14/2024] [Accepted: 02/04/2024] [Indexed: 03/21/2024] Open
Abstract
Mechanical stress is an internal force between various parts of an object that resists external factors and effects that cause an object to deform, and mechanical stress is essential for various tissues that are constantly subjected to mechanical loads to function normally. Integrins are a class of transmembrane heterodimeric glycoprotein receptors that are important target proteins for the action of mechanical stress stimuli on cells and can convert extracellular physical and mechanical signals into intracellular bioelectrical signals, thereby regulating osteogenesis and osteolysis. Integrins play a bidirectional regulatory role in bone metabolism. In this paper, relevant literature published in recent years is reviewed and summarized. The characteristics of integrins and mechanical stress are introduced, as well as the mechanisms underlying responses of integrin to mechanical stress stimulation. The paper focuses on integrin-mediated mechanical stress in different cells involved in bone metabolism and its associated signalling mechanisms. The purpose of this review is to provide a theoretical basis for the application of integrin-mediated mechanical stress to the field of bone tissue repair and regeneration.
Collapse
Affiliation(s)
- Li Yang
- Department of Periodontology, Hospital of StomatologyZunyi Medical UniversityZunyiChina
| | - Hong Chen
- Department of Periodontology, Hospital of StomatologyZunyi Medical UniversityZunyiChina
| | - Chanchan Yang
- Department of Periodontology, Hospital of StomatologyZunyi Medical UniversityZunyiChina
| | - Zhengqi Hu
- Department of Periodontology, Hospital of StomatologyZunyi Medical UniversityZunyiChina
| | - Zhiliang Jiang
- Department of Periodontology, Hospital of StomatologyZunyi Medical UniversityZunyiChina
| | - Shengzi Meng
- Department of Periodontology, Hospital of StomatologyZunyi Medical UniversityZunyiChina
| | | | - Lan Huang
- Department of Periodontology, Hospital of StomatologyZunyi Medical UniversityZunyiChina
| | | |
Collapse
|
5
|
Miller CJ, Pickering E, Martelli S, Dall'Ara E, Delisser P, Pivonka P. Cortical bone adaptation response is region specific, but not peak load dependent: insights from μ CT image analysis and mechanostat simulations of the mouse tibia loading model. Biomech Model Mechanobiol 2024; 23:287-304. [PMID: 37851203 PMCID: PMC10901956 DOI: 10.1007/s10237-023-01775-6] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/06/2023] [Accepted: 09/20/2023] [Indexed: 10/19/2023]
Abstract
The two major aims of the present study were: (i) quantify localised cortical bone adaptation at the surface level using contralateral endpoint imaging data and image analysis techniques, and (ii) investigate whether cortical bone adaptation responses are universal or region specific and dependent on the respective peak load. For this purpose, we re-analyse previously published μ CT data of the mouse tibia loading model that investigated bone adaptation in response to sciatic neurectomy and various peak load magnitudes (F = 0, 2, 4, 6, 8, 10, 12 N). A beam theory-based approach was developed to simulate cortical bone adaptation in different sections of the tibia, using longitudinal strains as the adaptive stimuli. We developed four mechanostat models: universal, surface-based, strain directional-based, and combined surface and strain direction-based. Rates of bone adaptation in these mechanostat models were computed using an optimisation procedure (131,606 total simulations), performed on a single load case (F = 10 N). Subsequently, the models were validated against the remaining six peak loads. Our findings indicate that local bone adaptation responses are quasi-linear and bone region specific. The mechanostat model which accounted for differences in endosteal and periosteal regions and strain directions (i.e. tensile versus compressive) produced the lowest root mean squared error between simulated and experimental data for all loads, with a combined prediction accuracy of 76.6, 55.0 and 80.7% for periosteal, endosteal, and cortical thickness measurements (in the midshaft of the tibia). The largest root mean squared errors were observed in the transitional loads, i.e. F = 2 to 6 N, where inter-animal variability was highest. Finally, while endpoint imaging studies provide great insights into organ level bone adaptation responses, the between animal and loaded versus control limb variability make simulations of local surface-based adaptation responses challenging.
Collapse
Affiliation(s)
- Corey J Miller
- School of Mechanical, Medical and Process Engineering, Queensland University of Technology (QUT), Brisbane, Australia.
| | - Edmund Pickering
- School of Mechanical, Medical and Process Engineering, Queensland University of Technology (QUT), Brisbane, Australia
| | - Saulo Martelli
- School of Mechanical, Medical and Process Engineering, Queensland University of Technology (QUT), Brisbane, Australia
| | - Enrico Dall'Ara
- Department of Oncology and Metabolism and Insigneo Institute for In Silico Medicine, University of Sheffield, Sheffield, UK
| | | | - Peter Pivonka
- School of Mechanical, Medical and Process Engineering, Queensland University of Technology (QUT), Brisbane, Australia.
| |
Collapse
|
6
|
Jafariandehkordi A, Daneshmehr A. Studying the mechanical properties of the mandible and injury prediction under the effect of ossification factors. J Mech Behav Biomed Mater 2023; 148:106209. [PMID: 37918338 DOI: 10.1016/j.jmbbm.2023.106209] [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: 08/27/2023] [Revised: 10/19/2023] [Accepted: 10/23/2023] [Indexed: 11/04/2023]
Abstract
BACKGROUND AND OBJECTIVE It is essential to know the quantitative interactions between biological tissues and external mechanical and chemical stimuli. This assists the physicians to better know the quantitative behavior of the tissue and plan more effective therapy. In the literature, the effect of the chemical and mechanical loading was investigated on the bone biological cell activities and some mechanical features, but a lack of prediction of bone injury under the chemical and mechanical factors was sensed. Therefore, the present study aims to investigate the effects of the application of major chemical factors involved in ossification, including RANKL1 (Receptor Activator of Nuclear Factor Kappa Beta Ligand), PTH2 (Parathyroid Hormone), and OPG3 (Osteoprotegerin) on the mandibular bone biological osteoblast and osteoclast activities and mechanical properties. Moreover, the study assesses the bone injury possibility under uniform mastication pressure applied on the premolar tooth in terms of the mechanostat theory undergoing the effects of the chemical factors. METHODS A 3D geometry of the mandible-tooth assembly was generated from the CT image dataset. The geometry was next purified, solidified, and exported to FEM4 (Finite Element Method) software to be meshed, where boundary conditions and loading were applied. Moreover, the mathematical system of differential equations to model the chemical factor effects on osteoblast and osteoclast activities as well as bone matrix volume fraction and elastic stiffness relations were applied. Next, the values of the equivalent strain were calculated to predict the injury states of the bone. RESULTS The results complied with the literature data. The results showed that RANKL and PTH increased the values of the equivalent strain from 450 με to 11500 με, while OPG reduced that from 450 με to 300 με. CONCLUSIONS Therefore, RANKL and PTH doses of this study put the bone at risk of injury compared to the baseline of no dose applied, while OPG secured the bone from injury.
Collapse
|
7
|
Wang Y, Dong H, Yan Y, Yu J, Wu X, Wang Y, Xue Y, Wang X, Wei X, Li P, Chen W. Biomechanical analysis of a lacunar-canalicular system under different cyclic displacement loading. Comput Methods Biomech Biomed Engin 2023; 26:1806-1821. [PMID: 36377250 DOI: 10.1080/10255842.2022.2145889] [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: 07/14/2022] [Accepted: 11/01/2022] [Indexed: 11/16/2022]
Abstract
The objective of this study is to use the finite element (FE) method to predict the mechanical signals (interstitial fluid velocity, strain, pore pressure, and pore fluid velocity) produced by osteocyte during physiological activities. The model predicts that the amplitude and distribution of the mechanical signals are mainly affected by the loading rate. The magnitude of mechanical signals in the lacunar-canalicular system increases as the amplitude, frequency and amount of direction of load increase. Collagen hillocks can effectively amplify strain signals at the process. The established model can be used for studying the mechanism of bone mechanotransduction at the micro-level.
Collapse
Affiliation(s)
- Yan Wang
- College of Biomedical Engineering, Taiyuan University of Technology, Taiyuan, China
| | - Hao Dong
- College of Biomedical Engineering, Taiyuan University of Technology, Taiyuan, China
| | - Yang Yan
- College of Biomedical Engineering, Taiyuan University of Technology, Taiyuan, China
| | - Jianhao Yu
- College of Biomedical Engineering, Taiyuan University of Technology, Taiyuan, China
| | - Xiaogang Wu
- College of Biomedical Engineering, Taiyuan University of Technology, Taiyuan, China
- Shanxi Provincial Key Laboratory for Repair of Bone and Soft Tissue Injury, Second hospital of Shanxi Medical University, Taiyuan, China
| | - Yanqin Wang
- College of Biomedical Engineering, Taiyuan University of Technology, Taiyuan, China
| | - Yanru Xue
- College of Biomedical Engineering, Taiyuan University of Technology, Taiyuan, China
| | - Xiyu Wang
- College of Biomedical Engineering, Taiyuan University of Technology, Taiyuan, China
| | - Xiaochun Wei
- Shanxi Provincial Key Laboratory for Repair of Bone and Soft Tissue Injury, Second hospital of Shanxi Medical University, Taiyuan, China
| | - Pengcui Li
- Shanxi Provincial Key Laboratory for Repair of Bone and Soft Tissue Injury, Second hospital of Shanxi Medical University, Taiyuan, China
| | - Weiyi Chen
- College of Biomedical Engineering, Taiyuan University of Technology, Taiyuan, China
| |
Collapse
|
8
|
Ait Oumghar I, Barkaoui A, Ghazi AE, Chabrand P. Modeling and simulation of bone cells dynamic behavior under the late effect of breast cancer treatments. Med Eng Phys 2023; 115:103982. [PMID: 37120177 DOI: 10.1016/j.medengphy.2023.103982] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/26/2022] [Revised: 04/14/2023] [Accepted: 04/19/2023] [Indexed: 05/01/2023]
Abstract
Breast Cancer (BC) treatments have been proven to interfere with the health of bones. Chemotherapy and endocrinal treatment regimens such as tamoxifen and aromatase inhibitors are frequently prescribed for women with BC. However, these drugs increase bone resorption and reduce the Bone Mineral Density (BMD), thus increasing the risk of bone fracture. In the current study, a mechanobiological bone remodeling model has been developed by coupling cellular activities, mechanical stimuli, and the effect of breast cancer treatments (chemotherapy, tamoxifen, and aromatase inhibitors). This model algorithm has been programmed and implemented on MATLAB software to simulate different treatment scenarios and their effects on bone remodeling and also predict the evolution of Bone Volume fraction (BV/TV) and the associated Bone Density Loss (BDL) over a period of time. The simulation results, achieved from different combinations of Breast Cancer treatments, allow the researchers to predict the intensity of each combination treatment on BV/TV and BMD. The combination of chemotherapy, tamoxifen, and aromatase inhibitors, followed by the combination of chemotherapy and tamoxifen remain the most harmful regimen. This is because they have a strong ability to induce the bone degradation which is represented by a decrease of 13.55% and 11.55% of the BV/TV value, respectively. These results were compared with the experimental studies and clinical observations which showed good agreement. The proposed model can be used by clinicians and physicians to choose the most appropriate combination of treatments, according to the patient's case.
Collapse
Affiliation(s)
- Imane Ait Oumghar
- Université Internationale de Rabat, LERMA Lab, Rocade Rabat Salé 11100, Rabat-Sala El Jadida, Morocco; Université Aix-Marseille, ISM, 163 av. de Luminy F-13288, Marseille cedex 09, France
| | - Abdelwahed Barkaoui
- Université Internationale de Rabat, LERMA Lab, Rocade Rabat Salé 11100, Rabat-Sala El Jadida, Morocco.
| | - Abdellatif El Ghazi
- Université Internationale de Rabat, TIC Lab, Rocade Rabat Salé 11100, Rabat-Sala El Jadida, Morocco
| | - Patrick Chabrand
- Université Aix-Marseille, ISM, 163 av. de Luminy F-13288, Marseille cedex 09, France
| |
Collapse
|
9
|
Goto S, Zhu Q, Wildsoet CF. Advanced ocular heterotopic bone formation in a guinea pig. VETERINARY RECORD CASE REPORTS 2023. [DOI: 10.1002/vrc2.594] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/17/2023]
Affiliation(s)
- So Goto
- Herbert Wertheim School Optometry and Vision Science University of California Berkeley California USA
- Department of Ophthalmology Osaka University Graduate School of Medicine Suita Osaka Japan
- Department of Ophthalmology National Hospital Organization Tokyo Medical Center Meguro‐ku Tokyo Japan
| | - Qiurong Zhu
- Herbert Wertheim School Optometry and Vision Science University of California Berkeley California USA
- Schepens Eye Research Institute Massachusetts Eye and Ear Harvard Medical School Boston Massachusetts USA
| | - Christine F. Wildsoet
- Herbert Wertheim School Optometry and Vision Science University of California Berkeley California USA
| |
Collapse
|
10
|
Computational modeling for osteogenic potential assessment of physical exercises based on loading-induced mechanobiological environments in cortical bone remodeling. Biomech Model Mechanobiol 2023; 22:281-295. [PMID: 36305993 DOI: 10.1007/s10237-022-01647-5] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/19/2022] [Accepted: 10/08/2022] [Indexed: 11/02/2022]
Abstract
Osteoporosis and disuse can cause bone loss which reduces the weight-bearing strength of long bones. Physical exercise or mechanical loading prevents bone loss as it promotes bone modeling through osteogenesis, i.e., new bone formation. Several studies have observed distinct bone remodeling responses to physical exercises; nevertheless, the underlying mechanism behind such responses is not well established. Loading-induced pore-pressure and fluid motion act as mechanobiological stimuli to bone cells namely osteocytes which further initiate osteoactivities. The shape of loading waveforms also affects the poromechanical environment of bone. Accordingly, the present study hypothesizes that loading waveforms associated with physiological exercises may expose the bone to different mechanobiological stimuli resulting in distinct bone remodeling. A poromechanical finite element model is developed to compute pore-pressure and interstitial fluid velocity in femoral cortical bone tissue (healthy and osteoporotic) subjected to loading waveforms of three physiological exercises namely walking, running, and jumping. The model also computes the mechanobiological stimulus as a function of fluid velocity. The outcomes indicate that pore-pressure and fluid velocity decrease significantly in osteoporotic bone tissue in comparison with healthy tissue. Jumping and running both improve pore-pressure and fluid velocity in healthy and osteoporotic tissues, whereas running significantly enhances mechanobiological stimulus in both the tissues which indicates a possible explanation for distinct bone remodeling to different physical exercises. The present work also suggests that running may be recommended as a potential biomechanical therapeutic to prevent bone loss. Overall, the present work contributes to the area of orthopedic research to develop effective designs of prophylactic exercises to improve bone health.
Collapse
|
11
|
Boaretti D, Marques FC, Ledoux C, Singh A, Kendall JJ, Wehrle E, Kuhn GA, Bansod YD, Schulte FA, Müller R. Trabecular bone remodeling in the aging mouse: A micro-multiphysics agent-based in silico model using single-cell mechanomics. Front Bioeng Biotechnol 2023; 11:1091294. [PMID: 36937760 PMCID: PMC10017748 DOI: 10.3389/fbioe.2023.1091294] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/06/2022] [Accepted: 02/15/2023] [Indexed: 03/06/2023] Open
Abstract
Bone remodeling is regulated by the interaction between different cells and tissues across many spatial and temporal scales. Notably, in silico models are regarded as powerful tools to further understand the signaling pathways that regulate this intricate spatial cellular interplay. To this end, we have established a 3D multiscale micro-multiphysics agent-based (micro-MPA) in silico model of trabecular bone remodeling using longitudinal in vivo data from the sixth caudal vertebra (CV6) of PolgA(D257A/D257A) mice, a mouse model of premature aging. Our in silico model includes a variety of cells as single agents and receptor-ligand kinetics, mechanomics, diffusion and decay of cytokines which regulate the cells' behavior. We highlighted its capabilities by simulating trabecular bone remodeling in the CV6 of five mice over 4 weeks and we evaluated the static and dynamic morphometry of the trabecular bone microarchitecture. Based on the progression of the average trabecular bone volume fraction (BV/TV), we identified a configuration of the model parameters to simulate homeostatic trabecular bone remodeling, here named basal. Crucially, we also produced anabolic, anti-anabolic, catabolic and anti-catabolic responses with an increase or decrease by one standard deviation in the levels of osteoprotegerin (OPG), receptor activator of nuclear factor kB ligand (RANKL), and sclerostin (Scl) produced by the osteocytes. Our results showed that changes in the levels of OPG and RANKL were positively and negatively correlated with the BV/TV values after 4 weeks in comparison to basal levels, respectively. Conversely, changes in Scl levels produced small fluctuations in BV/TV in comparison to the basal state. From these results, Scl was deemed to be the main driver of equilibrium while RANKL and OPG were shown to be involved in changes in bone volume fraction with potential relevance for age-related bone features. Ultimately, this micro-MPA model provides valuable insights into how cells respond to their local mechanical environment and can help to identify critical pathways affected by degenerative conditions and ageing.
Collapse
Affiliation(s)
| | | | - Charles Ledoux
- Institute for Biomechanics, ETH Zurich, Zurich, Switzerland
| | - Amit Singh
- Institute for Biomechanics, ETH Zurich, Zurich, Switzerland
| | | | - Esther Wehrle
- Institute for Biomechanics, ETH Zurich, Zurich, Switzerland
- AO Research Institute Davos, Davos Platz, Switzerland
| | - Gisela A. Kuhn
- Institute for Biomechanics, ETH Zurich, Zurich, Switzerland
| | | | | | - Ralph Müller
- Institute for Biomechanics, ETH Zurich, Zurich, Switzerland
- *Correspondence: Ralph Müller,
| |
Collapse
|
12
|
Ait Oumghar I, Barkaoui A, Chabrand P, Ghazi AE, Jeanneau C, Guenoun D, Pivonka P. Experimental-based mechanobiological modeling of the anabolic and catabolic effects of breast cancer on bone remodeling. Biomech Model Mechanobiol 2022; 21:1841-1856. [PMID: 36001274 DOI: 10.1007/s10237-022-01623-z] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/23/2022] [Accepted: 08/02/2022] [Indexed: 11/02/2022]
Abstract
Bone is a biological tissue characterized by its hierarchical organization. This material has the ability to be continually renewed, which makes it highly adaptative to external loadings. Bone renewing is managed by a dynamic biological process called bone remodeling (BR), where continuous resorption of old bone and formation of new bone permits to change the bone composition and microstructure. Unfortunately, because of several factors, such as age, hormonal imbalance, and a variety of pathologies including cancer metastases, this process can be disturbed leading to various bone diseases. In this study, we have investigated the effect of breast cancer (BC) metastases causing osteolytic bone loss. BC has the ability to affect bone quantity in different ways in each of its primary and secondary stages. Based on a BR mathematical model, we modeled the BC cells' interaction with bone cells to assess their effect on bone volume fraction (BV/TV) evolution during the remodeling process. Some of the parameters used in our model have been determined experimentally using the enzyme-linked immune-sorbent assay (ELISA) and the MTT assay. Our numerical simulations show that primary BC plays a significant role in enhancing bone-forming cells' activity leading to a 6.22% increase in BV/TV over 1 year. On the other hand, secondary BC causes a noticeable decrease in BV/TV reaching 15.74% over 2 years.
Collapse
Affiliation(s)
- Imane Ait Oumghar
- LERMA Lab, Université Internationale de Rabat, Rocade Rabat Salé 11100, Rabat-Sala El Jadida, Morocco.,Université Aix-Marseille, ISM, 163 av. de Luminy, 13288, Marseille Cedex 09, France
| | - Abdelwahed Barkaoui
- LERMA Lab, Université Internationale de Rabat, Rocade Rabat Salé 11100, Rabat-Sala El Jadida, Morocco.
| | - Patrick Chabrand
- Université Aix-Marseille, ISM, 163 av. de Luminy, 13288, Marseille Cedex 09, France
| | - Abdellatif El Ghazi
- TIC Lab, Université Internationale de Rabat, Rocade Rabat Salé 11100, Rabat-Sala El Jadida, Morocco
| | - Charlotte Jeanneau
- Université Aix-Marseille, ISM, 163 av. de Luminy, 13288, Marseille Cedex 09, France
| | - Daphne Guenoun
- Université Aix-Marseille, ISM, 163 av. de Luminy, 13288, Marseille Cedex 09, France
| | - Peter Pivonka
- Biomechanics and Spine Research Group, Queensland University of Technology at the Centre for Children's Health Research, South Brisbane, 4101, QLD, Australia
| |
Collapse
|
13
|
Ledoux C, Boaretti D, Sachan A, Müller R, Collins CJ. Clinical Data for Parametrization of In Silico Bone Models Incorporating Cell-Cytokine Dynamics: A Systematic Review of Literature. Front Bioeng Biotechnol 2022; 10:901720. [PMID: 35910035 PMCID: PMC9335409 DOI: 10.3389/fbioe.2022.901720] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/22/2022] [Accepted: 06/13/2022] [Indexed: 11/13/2022] Open
Abstract
In silico simulations aim to provide fast, inexpensive, and ethical alternatives to years of costly experimentation on animals and humans for studying bone remodeling, its deregulation during osteoporosis and the effect of therapeutics. Within the varied spectrum of in silico modeling techniques, bone cell population dynamics and agent-based multiphysics simulations have recently emerged as useful tools to simulate the effect of specific signaling pathways. In these models, parameters for cell and cytokine behavior are set based on experimental values found in literature; however, their use is currently limited by the lack of clinical in vivo data on cell numbers and their behavior as well as cytokine concentrations, diffusion, decay and reaction rates. Further, the settings used for these parameters vary across research groups, prohibiting effective cross-comparisons. This review summarizes and evaluates the clinical trial literature that can serve as input or validation for in silico models of bone remodeling incorporating cells and cytokine dynamics in post-menopausal women in treatment, and control scenarios. The GRADE system was used to determine the level of confidence in the reported data, and areas lacking in reported measures such as binding site occupancy, reaction rates and cell proliferation, differentiation and apoptosis rates were highlighted as targets for further research. We propose a consensus for the range of values that can be used for the cell and cytokine settings related to the RANKL-RANK-OPG, TGF-β and sclerostin pathways and a Levels of Evidence-based method to estimate parameters missing from clinical trial literature.
Collapse
Affiliation(s)
- Charles Ledoux
- Institute for Biomechanics, ETH Zurich, Zurich, Switzerland
| | | | - Akanksha Sachan
- Institute for Biomechanics, ETH Zurich, Zurich, Switzerland
- Department of Chemical Engineering, Indian Institute of Technology Bombay, Mumbai, India
| | - Ralph Müller
- Institute for Biomechanics, ETH Zurich, Zurich, Switzerland
| | - Caitlyn J. Collins
- Institute for Biomechanics, ETH Zurich, Zurich, Switzerland
- Department for Biomedical Engineering and Mechanics, Virginia Tech, Blacksburg, VI,United States
- *Correspondence: Caitlyn J. Collins,
| |
Collapse
|
14
|
Rayat Pisheh H, Ansari M, Eslami H. How is mechanobiology involved in bone regenerative medicine? Tissue Cell 2022; 76:101821. [DOI: 10.1016/j.tice.2022.101821] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/09/2022] [Revised: 04/27/2022] [Accepted: 05/10/2022] [Indexed: 10/18/2022]
|
15
|
Liu P, Tu J, Wang W, Li Z, Li Y, Yu X, Zhang Z. Effects of Mechanical Stress Stimulation on Function and Expression Mechanism of Osteoblasts. Front Bioeng Biotechnol 2022; 10:830722. [PMID: 35252138 PMCID: PMC8893233 DOI: 10.3389/fbioe.2022.830722] [Citation(s) in RCA: 13] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/07/2021] [Accepted: 01/10/2022] [Indexed: 12/13/2022] Open
Abstract
Osteoclasts and osteoblasts play a major role in bone tissue homeostasis. The homeostasis and integrity of bone tissue are maintained by ensuring a balance between osteoclastic and osteogenic activities. The remodeling of bone tissue is a continuous ongoing process. Osteoclasts mainly play a role in bone resorption, whereas osteoblasts are mainly involved in bone remodeling processes, such as bone cell formation, mineralization, and secretion. These cell types balance and restrict each other to maintain bone tissue metabolism. Bone tissue is very sensitive to mechanical stress stimulation. Unloading and loading of mechanical stress are closely related to the differentiation and formation of osteoclasts and bone resorption function as well as the differentiation and formation of osteoblasts and bone formation function. Consequently, mechanical stress exerts an important influence on the bone microenvironment and bone metabolism. This review focuses on the effects of different forms of mechanical stress stimulation (including gravity, continuously compressive pressure, tensile strain, and fluid shear stress) on osteoclast and osteoblast function and expression mechanism. This article highlights the involvement of osteoclasts and osteoblasts in activating different mechanical transduction pathways and reports changings in their differentiation, formation, and functional mechanism induced by the application of different types of mechanical stress to bone tissue. This review could provide new ideas for further microscopic studies of bone health, disease, and tissue damage reconstruction.
Collapse
Affiliation(s)
- Pan Liu
- School of Clinical Medicine, Chengdu Medical College, Chengdu, China
- The First Affiliated Hospital of Chengdu Medical College, Chengdu, China
| | - Ji Tu
- Spine Labs, St. George & Sutherland Clinical School, University of New South Wales, Sydney, NSW, Australia
| | - Wenzhao Wang
- Department of Orthopedics, West China Hospital of Sichuan University, Chengdu, China
| | - Zheng Li
- People’s Hospital of Jiulongpo District, Chongqing, China
| | - Yao Li
- School of Clinical Medicine, Chengdu Medical College, Chengdu, China
- The First Affiliated Hospital of Chengdu Medical College, Chengdu, China
| | - Xiaoping Yu
- School of Public Health, Chengdu Medical College, Chengdu, China
- Basic Medical College of Chengdu University, Chengdu, China
- *Correspondence: Xiaoping Yu, ; Zhengdong Zhang,
| | - Zhengdong Zhang
- School of Clinical Medicine, Chengdu Medical College, Chengdu, China
- The First Affiliated Hospital of Chengdu Medical College, Chengdu, China
- Department of Orthopedics, The First Affiliated Hospital of Chengdu Medical College, Chengdu, China
- *Correspondence: Xiaoping Yu, ; Zhengdong Zhang,
| |
Collapse
|
16
|
Morin C, Hellmich C, Nejim Z, Avril S. Fiber Rearrangement and Matrix Compression in Soft Tissues: Multiscale Hypoelasticity and Application to Tendon. Front Bioeng Biotechnol 2021; 9:725047. [PMID: 34712652 PMCID: PMC8546211 DOI: 10.3389/fbioe.2021.725047] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/14/2021] [Accepted: 09/06/2021] [Indexed: 11/13/2022] Open
Abstract
It is widely accepted that the nonlinear macroscopic mechanical behavior of soft tissue is governed by fiber straightening and re-orientation. Here, we provide a quantitative assessment of this phenomenon, by means of a continuum micromechanics approach. Given the negligibly small bending stiffness of crimped fibers, the latter are represented through a number of hypoelastic straight fiber phases with different orientations, being embedded into a hypoelastic matrix phase. The corresponding representative volume element (RVE) hosting these phases is subjected to “macroscopic” strain rates, which are downscaled to fiber and matrix strain rates on the one hand, and to fiber spins on the other hand. This gives quantitative access to the fiber decrimping (or straightening) phenomenon under non-affine conditions, i.e. in the case where the fiber orientations cannot be simply linked to the macroscopic strain state. In the case of tendinous tissue, such an RVE relates to the fascicle material with 50 μm characteristic length, made up of crimped collagen bundles and a gel-type matrix in-between. The fascicles themselves act as parallel fibers in a similar matrix at the scale of a tissue-related RVE with 500 μm characteristic length. As evidenced by a sensitivity analysis and confirmed by various mechanical tests, it is the initial crimping angle which drives both the degree of straightening and the shape of the macroscopic stress-strain curve, while the final linear portion of this curve depends almost exclusively on the collagen bundle elasticity. Our model also reveals the mechanical cooperation of the tissue’s key microstructural components: while the fibers carry tensile forces, the matrices undergo hydrostatic pressure.
Collapse
Affiliation(s)
- Claire Morin
- Mines Saint-Etienne, Univ. Lyon, Univ. Jean Monnet, INSERM, U1059 Sainbiose, Centre CIS, Saint-Etienne, France
| | - Christian Hellmich
- Institute for Mechanics of Materials and Structures, TU Wien - Vienna University of Technology, Vienna, Austria
| | - Zeineb Nejim
- Mines Saint-Etienne, Univ. Lyon, Univ. Jean Monnet, INSERM, U1059 Sainbiose, Centre CIS, Saint-Etienne, France
| | - Stéphane Avril
- Mines Saint-Etienne, Univ. Lyon, Univ. Jean Monnet, INSERM, U1059 Sainbiose, Centre CIS, Saint-Etienne, France.,Institute for Mechanics of Materials and Structures, TU Wien - Vienna University of Technology, Vienna, Austria
| |
Collapse
|
17
|
García-Aznar JM, Nasello G, Hervas-Raluy S, Pérez MÁ, Gómez-Benito MJ. Multiscale modeling of bone tissue mechanobiology. Bone 2021; 151:116032. [PMID: 34118446 DOI: 10.1016/j.bone.2021.116032] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 01/13/2021] [Revised: 04/25/2021] [Accepted: 06/02/2021] [Indexed: 02/07/2023]
Abstract
Mechanical environment has a crucial role in our organism at the different levels, ranging from cells to tissues and our own organs. This regulatory role is especially relevant for bones, given their importance as load-transmitting elements that allow the movement of our body as well as the protection of vital organs from load impacts. Therefore bone, as living tissue, is continuously adapting its properties, shape and repairing itself, being the mechanical loads one of the main regulatory stimuli that modulate this adaptive behavior. Here we review some key results of bone mechanobiology from computational models, describing the effect that changes associated to the mechanical environment induce in bone response, implant design and scaffold-driven bone regeneration.
Collapse
Affiliation(s)
- José Manuel García-Aznar
- Multiscale in Mechanical and Biological Engineering, Instituto de Investigación en Ingeniería de Aragón (I3A), Instituto de Investigación Sanitaria Aragón (IIS Aragón), University of Zaragoza, Zaragoza, Spain.
| | - Gabriele Nasello
- Multiscale in Mechanical and Biological Engineering, Instituto de Investigación en Ingeniería de Aragón (I3A), Instituto de Investigación Sanitaria Aragón (IIS Aragón), University of Zaragoza, Zaragoza, Spain; Biomechanics Section, KU Leuven, Leuven, Belgium
| | - Silvia Hervas-Raluy
- Multiscale in Mechanical and Biological Engineering, Instituto de Investigación en Ingeniería de Aragón (I3A), Instituto de Investigación Sanitaria Aragón (IIS Aragón), University of Zaragoza, Zaragoza, Spain
| | - María Ángeles Pérez
- Multiscale in Mechanical and Biological Engineering, Instituto de Investigación en Ingeniería de Aragón (I3A), Instituto de Investigación Sanitaria Aragón (IIS Aragón), University of Zaragoza, Zaragoza, Spain
| | - María José Gómez-Benito
- Multiscale in Mechanical and Biological Engineering, Instituto de Investigación en Ingeniería de Aragón (I3A), Instituto de Investigación Sanitaria Aragón (IIS Aragón), University of Zaragoza, Zaragoza, Spain
| |
Collapse
|
18
|
Larcher I, Scheiner S. Parameter reduction, sensitivity studies, and correlation analyses applied to a mechanobiologically regulated bone cell population model of the bone metabolism. Comput Biol Med 2021; 136:104717. [PMID: 34426166 DOI: 10.1016/j.compbiomed.2021.104717] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/12/2021] [Revised: 07/27/2021] [Accepted: 07/27/2021] [Indexed: 11/30/2022]
Abstract
When striving for reconstructing and predicting bone remodeling processes by means of mathematical models, cell population models have become a popular option. From a conceptual point of view, these models are able to take into account an arbitrary amount of regulatory mechanisms driving the development of bone cells and their activities. However, in most cases, the models include a large number of parameters; and most of those parameters cannot be measured, which certainly compromises the credibility of cell population models. Here, new insights are presented as to the potential improvement of this unsatisfactory situation. In particular, a previously published bone remodeling model was considered, and based on combination and merging of the original parameters, the total number of parameters could be reduced from 28 to 18, without impairing the model's versatility and significance. Furthermore, a comprehensive number of one- and two-variable sensitivity studies were performed, pointing out which parameters (alone and in combination with other parameters) influence the model predictions significantly - for that purpose, the mean squared relative error (MSRE) between simulations based on the original parameters and based on varied parameters was considered as failure measure. It has turned out that the model is significantly more sensitive to parameters which can be considered as phenomenological (such as differentiation, proliferation, and apoptosis rates) than to parameters which are directly related to specific processes (such as dissociation rate constants, and maximum concentrations of the involved factors). Using common correlation measures (such as Pearson, Spearman, and partial ranked correlation coefficients), correlation studies revealed that the correlations between most parameters and the MSRE are weak, while a few parameters exhibited moderate correlations. In conclusion, the results shown in this paper provide valuable insights concerning the design of new experiments allowing for measurement of the parameters which are most influential in the context of bone remodeling simulation.
Collapse
Affiliation(s)
- Isabella Larcher
- Institute for Mechanics of Materials and Structures, (TU Wien) Vienna University of Technology, Karlsplatz 13/202, 1040, Vienna, Austria
| | - Stefan Scheiner
- Institute for Mechanics of Materials and Structures, (TU Wien) Vienna University of Technology, Karlsplatz 13/202, 1040, Vienna, Austria.
| |
Collapse
|
19
|
A General Mechano-Pharmaco-Biological Model for Bone Remodeling Including Cortisol Variation. MATHEMATICS 2021. [DOI: 10.3390/math9121401] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
The process of bone remodeling requires a strict coordination of bone resorption and formation in time and space in order to maintain consistent bone quality and quantity. Bone-resorbing osteoclasts and bone-forming osteoblasts are the two major players in the remodeling process. Their coordination is achieved by generating the appropriate number of osteoblasts since osteoblastic-lineage cells govern the bone mass variation and regulate a corresponding number of osteoclasts. Furthermore, diverse hormones, cytokines and growth factors that strongly link osteoblasts to osteoclasts coordinated these two cell populations. The understanding of this complex remodeling process and predicting its evolution is crucial to manage bone strength under physiologic and pathologic conditions. Several mathematical models have been suggested to clarify this remodeling process, from the earliest purely phenomenological to the latest biomechanical and mechanobiological models. In this current article, a general mathematical model is proposed to fill the gaps identified in former bone remodeling models. The proposed model is the result of combining existing bone remodeling models to present an updated model, which also incorporates several important parameters affecting bone remodeling under various physiologic and pathologic conditions. Furthermore, the proposed model can be extended to include additional parameters in the future. These parameters are divided into four groups according to their origin, whether endogenous or exogenous, and the cell population they affect, whether osteoclasts or osteoblasts. The model also enables easy coupling of biological models to pharmacological and/or mechanical models in the future.
Collapse
|
20
|
Pant A, Paul E, Niebur GL, Vahdati A. Integration of mechanics and biology in computer simulation of bone remodeling. PROGRESS IN BIOPHYSICS AND MOLECULAR BIOLOGY 2021; 164:33-45. [PMID: 33965425 DOI: 10.1016/j.pbiomolbio.2021.05.001] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/16/2021] [Revised: 03/27/2021] [Accepted: 05/03/2021] [Indexed: 12/14/2022]
Abstract
Bone remodeling is a complex physiological process that spans across multiple spatial and temporal scales and is regulated by both mechanical and hormonal cues. An imbalance between bone resorption and bone formation in the process of bone remodeling may lead to various bone pathologies. One powerful and non-invasive approach to gain new insights into mechano-adaptive bone remodeling is computer modeling and simulation. Recent findings in bone physiology and advances in computer modeling have provided a unique opportunity to study the integration of mechanics and biology in bone remodeling. Our objective in this review is to critically appraise recent advances and developments and discuss future research opportunities in computational bone remodeling approaches that enable integration of mechanics and cellular and molecular pathways. Based on the critical appraisal of the relevant recent published literature, we conclude that multiscale in silico integration of personalized bone mechanics and mechanobiology combined with data science and analytics techniques offer the potential to deepen our knowledge of bone remodeling and provide ample opportunities for future research.
Collapse
Affiliation(s)
- Anup Pant
- Multi-disciplinary Mechanics and Modeling Laboratory, Department of Engineering, East Carolina University, Greenville, NC 27858, USA
| | - Elliot Paul
- Multi-disciplinary Mechanics and Modeling Laboratory, Department of Engineering, East Carolina University, Greenville, NC 27858, USA
| | - Glen L Niebur
- Tissue Mechanics Laboratory, Bioengineering Graduate Program, University of Notre Dame, Notre Dame, IN 46556, USA
| | - Ali Vahdati
- Multi-disciplinary Mechanics and Modeling Laboratory, Department of Engineering, East Carolina University, Greenville, NC 27858, USA.
| |
Collapse
|
21
|
Load adaptation through bone remodeling: a mechanobiological model coupled with the finite element method. Biomech Model Mechanobiol 2021; 20:1495-1507. [PMID: 33900492 DOI: 10.1007/s10237-021-01458-0] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/11/2020] [Accepted: 04/05/2021] [Indexed: 10/21/2022]
Abstract
This work proposes a novel tissue-scale mechanobiological model of bone remodeling to study bone's adaptation to distinct loading conditions. The devised algorithm describes the mechanosensitivity of bone and its impact on bone cells' functioning through distinct signaling factors. In this study, remodeling is mechanically ruled by variations of the strain energy density (SED) of bone, which is determined by performing a linear elastostatic analysis combined with the finite element method. Depending on the SED levels and on a set of biological signaling factors ([Formula: see text] parameters), osteoclasts and osteoblasts can be mechanically triggered. To reproduce this phenomenon, this work proposes a new set of [Formula: see text] parameters. The combined response of osteoclasts and osteoblasts will then affect bone's apparent density, which is correlated with other mechanical properties of bone, through a phenomenological law. Thus, this novel model proposes a constant interplay between the mechanical and biological components of the process. The spatiotemporal simulation used to validate this new approach is a benchmark example composed by two distinct phases: (1) pre-orientation and (2) load adaptation. On both of them, bone is able to adapt its morphology according to the loading condition, achieving the required trabecular distribution to withstand the applied loads. Moreover, the equilibrium morphology reflects the orientation of the load. These preliminary results support the new approach proposed in this study.
Collapse
|
22
|
Subject-specific multiscale modeling of aortic valve biomechanics. Biomech Model Mechanobiol 2021; 20:1031-1046. [PMID: 33792805 PMCID: PMC8154826 DOI: 10.1007/s10237-021-01429-5] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/06/2020] [Accepted: 01/28/2021] [Indexed: 11/13/2022]
Abstract
A Finite Element workflow for the multiscale analysis of the aortic valve biomechanics was developed and applied to three physiological anatomies with the aim of describing the aortic valve interstitial cells biomechanical milieu in physiological conditions, capturing the effect of subject-specific and leaflet-specific anatomical features from the organ down to the cell scale. A mixed approach was used to transfer organ-scale information down to the cell-scale. Displacement data from the organ model were used to impose kinematic boundary conditions to the tissue model, while stress data from the latter were used to impose loading boundary conditions to the cell level. Peak of radial leaflet strains was correlated with leaflet extent variability at the organ scale, while circumferential leaflet strains varied over a narrow range of values regardless of leaflet extent. The dependency of leaflet biomechanics on the leaflet-specific anatomy observed at the organ length-scale is reflected, and to some extent emphasized, into the results obtained at the lower length-scales. At the tissue length-scale, the peak diastolic circumferential and radial stresses computed in the fibrosa correlated with the leaflet surface area. At the cell length-scale, the difference between the strains in two main directions, and between the respective relationships with the specific leaflet anatomy, was even more evident; cell strains in the radial direction varied over a relatively wide range (\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$0.36-0.87$$\end{document}0.36-0.87) with a strong correlation with the organ length-scale radial strain (\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$R^{2}= 0.95$$\end{document}R2=0.95); conversely, circumferential cell strains spanned a very narrow range (\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$0.75-0.88$$\end{document}0.75-0.88) showing no correlation with the circumferential strain at the organ level (\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$R^{2}= 0.02$$\end{document}R2=0.02). Within the proposed simulation framework, being able to account for the actual anatomical features of the aortic valve leaflets allowed to gain insight into their effect on the structural mechanics of the leaflets at all length-scales, down to the cell scale.
Collapse
|
23
|
Peyroteo MMA, Belinha J, Natal Jorge RM. A mathematical biomechanical model for bone remodeling integrated with a radial point interpolating meshless method. Comput Biol Med 2020; 129:104170. [PMID: 33352308 DOI: 10.1016/j.compbiomed.2020.104170] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/21/2020] [Revised: 12/03/2020] [Accepted: 12/04/2020] [Indexed: 11/18/2022]
Abstract
Bone remodeling is a highly complex process, in which bone cells interact and regulate bone's apparent density as a response to several external and internal stimuli. In this work, this process is numerically described using a novel 2D biomechanical model. Some of the new features in this model are (i) the mathematical parameters used to determine bone's apparent density and cellular density; (ii) an automatic boundary recognition step to spatially control bone remodeling and (iii) an approach to mimic the mechanical transduction to osteoclasts and osteoblasts. Moreover, this model is combined with a meshless approach - the Radial Point Interpolation Method (RPIM). The use of RPIM is an asset for this application, especially in the construction of the boundary maps. This work studies bone's adaptation to a certain loading regime through bone resorption. The signaling pathways of bone cells are dependent on the level of strain energy density (SED) in bone. So, when SED changes, bone cells' functioning is affected, causing also changes on bone's apparent density. With this model, bone is able to achieve an equilibrium state, optimizing its structure to withstand the applied loads. Results suggest that this model has the potential to provide high quality solutions while being a simpler alternative to more complex bone remodeling models in the literature.
Collapse
Affiliation(s)
- M M A Peyroteo
- INEGI, Institute of Science and Innovation in Mechanical and Industrial Engineering, Rua Dr. Roberto Frias, 400, 4200-465, Porto, Portugal.
| | - J Belinha
- School of Engineering, Polytechnic of Porto (ISEP), Mechanical Engineering Department, Rua Dr. António Bernardino de Almeida, 431, 4200-072, Porto, Portugal.
| | - R M Natal Jorge
- Faculty of Engineering of the University of Porto, Mechanical Engineering Department, FEUP, Rua Dr. Roberto Frias, S/N, 4200-465, Porto, Portugal.
| |
Collapse
|
24
|
Ait Oumghar I, Barkaoui A, Chabrand P. Toward a Mathematical Modeling of Diseases' Impact on Bone Remodeling: Technical Review. Front Bioeng Biotechnol 2020; 8:584198. [PMID: 33224935 PMCID: PMC7667152 DOI: 10.3389/fbioe.2020.584198] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/16/2020] [Accepted: 09/30/2020] [Indexed: 12/18/2022] Open
Abstract
A wide variety of bone diseases have hitherto been discovered, such as osteoporosis, Paget's disease, osteopetrosis, and metastatic bone disease, which are not well defined in terms of changes in biochemical and mechanobiological regulatory factors. Some of these diseases are secondary to other pathologies, including cancer, or to some clinical treatments. To better understand bone behavior and prevent its deterioration, bone biomechanics have been the subject of mathematical modeling that exponentially increased over the last years. These models are becoming increasingly complex. The current paper provides a timely and critical analysis of previously developed bone remodeling mathematical models, particularly those addressing bone diseases. Besides, mechanistic pharmacokinetic/pharmacodynamic (PK/PD) models, which englobe bone disease and its treatment's effect on bone health. Therefore, the review starts by presenting bone remodeling cycle and mathematical models describing this process, followed by introducing some bone diseases and discussing models of pathological mechanisms affecting bone, and concludes with exhibiting the available bone treatment procedures considered in the PK/PD models.
Collapse
Affiliation(s)
- Imane Ait Oumghar
- Laboratoire des Energies Renouvelables et Matériaux Avancés (LERMA), Université Internationale de Rabat, Rabat-Sala El Jadida, Morocco
- Aix Marseille Univ, CNRS, ISM, Inst Movement Sci, Marseille, France
| | - Abdelwahed Barkaoui
- Laboratoire des Energies Renouvelables et Matériaux Avancés (LERMA), Université Internationale de Rabat, Rabat-Sala El Jadida, Morocco
| | - Patrick Chabrand
- Aix Marseille Univ, CNRS, ISM, Inst Movement Sci, Marseille, France
| |
Collapse
|
25
|
Sego TJ, Hsu YT, Chu TM, Tovar A. Modeling Progressive Damage Accumulation in Bone Remodeling Explains the Thermodynamic Basis of Bone Resorption by Overloading. Bull Math Biol 2020; 82:134. [PMID: 33037933 DOI: 10.1007/s11538-020-00808-w] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/04/2020] [Accepted: 09/18/2020] [Indexed: 11/24/2022]
Abstract
Computational modeling of skeletal tissue seeks to predict the structural adaptation of bone in response to mechanical loading. The theory of continuum damage-repair, a mathematical description of structural adaptation based on principles of damage mechanics, continues to be developed and utilized for the prediction of long-term peri-implant outcomes. Despite its technical soundness, CDR does not account for the accumulation of mechanical damage and irreversible deformation. In this work, a nonlinear mathematical model of independent damage accumulation and plastic deformation is developed in terms of the CDR formulation. The proposed model incorporates empirical correlations from uniaxial experiments. Supporting elements of the model are derived, including damage and yielding criteria, corresponding consistency conditions, and the basic, necessary forms for integration during loading. Positivity of mechanical dissipation due to damage is proved, while strain-based, associative plastic flow and linear hardening describe post-yield behavior. Calibration of model parameters to the empirical correlations from which the model was derived is then accomplished. Results of numerical experiments on a point-wise specimen show that damage and plasticity inhibit bone formation by dissipation of energy available to biological processes, leading to material failure that would otherwise be predicted to experience a net gain of bone.
Collapse
Affiliation(s)
- T J Sego
- Department of Intelligent Systems Engineering, Indiana University, Bloomington, IN, USA.
| | - Yung-Ting Hsu
- Department of Periodontics, University of Washington, Seattle, WA, USA
| | - Tien-Min Chu
- Department of Restorative Dentistry, Indiana University, Indianapolis, IN, USA
| | - Andres Tovar
- Department of Mechanical and Energy Engineering, Indiana University-Purdue University Indianapolis, Indianapolis, IN, USA
| |
Collapse
|
26
|
Gubaua JE, Dicati GWO, Mercuri EGF, Pereira JT. Simulation of bone remodeling around a femoral prosthesis using a model that accounts for biological and mechanical interactions. Med Eng Phys 2020; 84:126-135. [PMID: 32977909 DOI: 10.1016/j.medengphy.2020.08.004] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/17/2019] [Revised: 08/10/2020] [Accepted: 08/15/2020] [Indexed: 11/28/2022]
Abstract
The present study focuses on a model for three-dimensional bone remodeling of the human femur that considers cellular dynamics to determine the volume fraction of new bone. The model considers the interaction among responsive osteoblasts, active osteoblasts, and osteoclasts, as well as signaling molecules and parathyroid hormone (PTH). The stimulus of the model has a systemic origin due to the PTH effect, and a local origin due to the action of cytokines, growth factors, and mechanical stimuli near the site of the bone cells. The present work considers that the mechanical stimulus that activates cellular activity is obtained from stresses acting on the bone tissue and the number of daily loading cycles. In addition to simulating the bone modeling process in an intact femur, the numerical model is used to simulate bone adaptation in relation to the stress shielding phenomenon after the implantation of a femoral prosthesis. The results showed that the simulations provide a distribution of bone density that is similar to a radiograph and, in addition, allows the visualization of osteoblast and osteoclast dynamics in bone adaptation response after prosthesis implantation.
Collapse
Affiliation(s)
- José Eduardo Gubaua
- Postgraduate Program in Mechanical Engineering, Federal University of Paraná, Paraná, Brazil.
| | - Gabriela Wessling Oening Dicati
- Postgraduate Program in Mechanical Engineering, Federal University of Paraná, Paraná, Brazil; Mechanical Engineering Department, Federal Technological University of Paraná, Campus Pato Branco, Paraná, Brazil
| | | | - Jucélio Tomás Pereira
- Postgraduate Program in Mechanical Engineering, Federal University of Paraná, Paraná, Brazil; Mechanical Engineering Department, Federal University of Paraná, Paraná, Brazil
| |
Collapse
|
27
|
Zelaya-Lainez L, Kariem H, Nischkauer W, Limbeck A, Hellmich C. "Variances" and "in-variances" in hierarchical porosity and composition, across femoral tissues from cow, horse, ostrich, emu, pig, rabbit, and frog. MATERIALS SCIENCE & ENGINEERING. C, MATERIALS FOR BIOLOGICAL APPLICATIONS 2020; 117:111234. [PMID: 32919621 DOI: 10.1016/j.msec.2020.111234] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/21/2019] [Revised: 06/10/2020] [Accepted: 06/21/2020] [Indexed: 10/24/2022]
Abstract
It is very well known that bone is a hierarchically organized material produced by bone cells residing in the fluid environments filling (larger) vascular pores and (smaller) lacunar pores. The extracellular space consists of hydroxyapatite crystals, collagen type I molecules, and water with non-collageneous organics. It is less known to which extent the associated quantities (mineral, organic, and water concentrations; vascular, lacunar, and extracellular porosities) vary across species, organs, and ages. We here investigate the aforementioned quantities across femoral shaft tissues from cow, horse, emu, frog, ostrich, pig, and rabbit; by means of light microscopy and dehydration-demineralization tests; thereby revealing interesting invariances: The extracellular volume fractions of organic matter turn out to be similar across all tested non-amphibian tissues; as do the extracellular volume fractions of hydroxyapatite across all tested mammals. Hence, the chemical composition of the femoral extracellular bone matrix is remarkably "invariant" across differently aged mammals; while the water content shows significant variations, as does the partitions of water between the different pore spaces. The latter exhibit strikingly varying morphologies as well. This finding adds to the ample "universal patterns" in the sense of evolutionary developmental biology; and it provides interesting design requirements for the development of novel biomimetic tissue engineering solutions.
Collapse
Affiliation(s)
- Luis Zelaya-Lainez
- Institute for Mechanics of Materials and Structures, TU Wien - Vienna University of Technology, Karlsplatz 13/E202, 1040 Vienna, Austria
| | - Hawraa Kariem
- Institute for Mechanics of Materials and Structures, TU Wien - Vienna University of Technology, Karlsplatz 13/E202, 1040 Vienna, Austria
| | - Winfried Nischkauer
- Institute of Chemical Technologies and Analytics, Division of Instrumental Analytical Chemistry, TU Wien - Vienna University of Technology, Getreidemarkt 9/164, 1060 Vienna, Austria
| | - Andreas Limbeck
- Institute of Chemical Technologies and Analytics, Division of Instrumental Analytical Chemistry, TU Wien - Vienna University of Technology, Getreidemarkt 9/164, 1060 Vienna, Austria
| | - Christian Hellmich
- Institute for Mechanics of Materials and Structures, TU Wien - Vienna University of Technology, Karlsplatz 13/E202, 1040 Vienna, Austria.
| |
Collapse
|
28
|
Sree VD, Tepole AB. Computational systems mechanobiology of growth and remodeling: Integration of tissue mechanics and cell regulatory network dynamics. CURRENT OPINION IN BIOMEDICAL ENGINEERING 2020. [DOI: 10.1016/j.cobme.2020.01.002] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
|
29
|
A mechano-chemo-biological model for bone remodeling with a new mechano-chemo-transduction approach. Biomech Model Mechanobiol 2020; 19:2499-2523. [PMID: 32623542 DOI: 10.1007/s10237-020-01353-0] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/04/2019] [Accepted: 05/29/2020] [Indexed: 12/26/2022]
Abstract
Bone remodeling is a fundamental biological process that develops in bone tissue along its whole lifetime. It refers to a continuous bone transformation with new bone formation and old bone resorption that changes the internal microstructure and composition of the tissue. The main objectives of bone remodeling are: repair of the internal microcracks; adaptation of the macroscopic stiffness and strength to the actual changing mechanical demands; and control of the calcium homeostasis. Understanding this process and predicting its evolution is critical to reduce the effects of long-term disuse as happens during periods of reduced mobility. It is also important in the design of bone implants to avoid long-term stress shielding. Many mathematical models have been proposed from the earliest purely phenomenological to the latest that include biological knowledge. However, there still exists a lack of connection between the mechanical driving force and the biochemical and cell processes it triggers. Here, and following previous works that model independently the mechanobiological and biochemical processes in bone remodeling, we present a more complete model, useful for both cortical and trabecular bone, that uses a new mechanotransduction approach based on the effect of strains onto the bonding-unbonding rate of RANK/RANKL/OPG receptor-ligand reactions. We compare the results of this model with previous ones, showing a good agreement in similar conditions. We also apply it to realistic situations such as a femoral bone after implantation of a hip prosthesis, getting similar results to the clinical ones in the final bone density distribution. Finally, we extend this approach to the anisotropic case, getting not only the mean density, but also the directional homogenization of the microstructure. This biochemical approach permits, not only to predict the bone evolution under changes in the mechanical loads, but also, to consider anabolic and catabolic drugs to control bone density, such as those used in osteoporosis.
Collapse
|
30
|
Strain shielding for cemented hip implants. Clin Biomech (Bristol, Avon) 2020; 77:105027. [PMID: 32447179 DOI: 10.1016/j.clinbiomech.2020.105027] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 12/15/2019] [Revised: 04/27/2020] [Accepted: 04/29/2020] [Indexed: 02/07/2023]
Abstract
BACKGROUND Long-term survival of hip implants is of increasing relevance due to the rising life expectancy. The biomechanical effect of strain shielding as a result of implant insertion may lead to bone resorption, thus increasing risk for implant loosening and periprosthetic fractures. Patient-specific quantification of strain shielding could assist orthopedic surgeons in choosing the biomechanically most appropriate prosthesis. METHODS Validated quantitative CT-based finite element models of five femurs in intact and implanted states were considered to propose a systematic algorithm for strain shielding quantification. Three different strain measures were investigated and the most appropriate measure for strain shielding quantification is recommended. It is used to demonstrate a practical femur-specific implant selection among three common designs. FINDINGS Strain shielding measures demonstrated similar trends in all Gruen zones except zone 1, where the volumetric strain measure differed from von-Mises and maximum principal strains. The volumetric strain measure is in better agreement with clinical bone resorption records. It is also consistent with the biological mechanism of bone remodeling so it is recommended for strain shielding quantification. Applying the strain shielding algorithm on three different implants for a specific femur suggests that the collared design is preferable. Such quantitative biomechanical input is valuable for practical patient specific implant selection. INTERPRETATION Volumetric strain should be considered for strain shielding examination. The presented methodology may potentially enable patient-specific pre-operative strain shielding evaluation so to minimize strain shielding. It should be further used in a longitudinal study so to correlate between strain shielding predictions and clinical bone resorption.
Collapse
|
31
|
Chen H, Han Q, Wang C, Liu Y, Chen B, Wang J. Porous Scaffold Design for Additive Manufacturing in Orthopedics: A Review. Front Bioeng Biotechnol 2020; 8:609. [PMID: 32626698 PMCID: PMC7311579 DOI: 10.3389/fbioe.2020.00609] [Citation(s) in RCA: 63] [Impact Index Per Article: 15.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/17/2020] [Accepted: 05/18/2020] [Indexed: 12/15/2022] Open
Abstract
With the increasing application of orthopedic scaffolds, a dramatically increasing number of requirements for scaffolds are precise. The porous structure has been a fundamental design in the bone tissue engineering or orthopedic clinics because of its low Young's modulus, high compressive strength, and abundant cell accommodation space. The porous structure manufactured by additive manufacturing (AM) technology has controllable pore size, pore shape, and porosity. The single unit can be designed and arrayed with AM, which brings controllable pore characteristics and mechanical properties. This paper presents the current status of porous designs in AM technology. The porous structures are stated from the cellular structure and the whole structure. In the aspect of the cellular structure, non-parametric design and parametric design are discussed here according to whether the algorithm generates the structure or not. The non-parametric design comprises the diamond, the body-centered cubic, and the polyhedral structure, etc. The Voronoi, the Triply Periodic Minimal Surface, and other parametric designs are mainly discussed in parametric design. In the discussion of cellular structures, we emphasize the design, and the resulting biomechanical and biological effects caused by designs. In the aspect of the whole structure, the recent experimental researches are reviewed on uniform design, layered gradient design, and layered gradient design based on topological optimization, etc. These parts are summarized because of the development of technology and the demand for mechanics or bone growth. Finally, the challenges faced by the porous designs and prospects of porous structure in orthopedics are proposed in this paper.
Collapse
Affiliation(s)
- Hao Chen
- Department of Orthopedics, Second Hospital of Jilin University, Changchun, China
| | - Qing Han
- Department of Orthopedics, Second Hospital of Jilin University, Changchun, China
| | - Chenyu Wang
- Department of Dermatology, The First Hospital of Jilin University, Changchun, China
| | - Yang Liu
- Department of Orthopedics, Second Hospital of Jilin University, Changchun, China
| | - Bingpeng Chen
- Department of Orthopedics, Second Hospital of Jilin University, Changchun, China
| | - Jincheng Wang
- Department of Orthopedics, Second Hospital of Jilin University, Changchun, China
| |
Collapse
|
32
|
Wang L, Shi Q, Cai Y, Chen Q, Guo X, Li Z. Mechanical–chemical coupled modeling of bone regeneration within a biodegradable polymer scaffold loaded with VEGF. Biomech Model Mechanobiol 2020; 19:2285-2306. [DOI: 10.1007/s10237-020-01339-y] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/28/2019] [Accepted: 05/12/2020] [Indexed: 10/24/2022]
|
33
|
Li Z, Betts D, Kuhn G, Schirmer M, Müller R, Ruffoni D. Mechanical regulation of bone formation and resorption around implants in a mouse model of osteopenic bone. J R Soc Interface 2020; 16:20180667. [PMID: 30890053 DOI: 10.1098/rsif.2018.0667] [Citation(s) in RCA: 22] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
Abstract
Although mechanical stimulation is considered a promising approach to accelerate implant integration, our understanding of load-driven bone formation and resorption around implants is still limited. This lack of knowledge may delay the development of effective loading protocols to prevent implant loosening, especially in osteoporosis. In healthy bone, formation and resorption are mechanoregulated processes. In the intricate context of peri-implant bone regeneration, it is not clear whether bone (re)modelling can still be load-driven. Here, we investigated the mechanical control of peri-implant bone (re)modelling with a well-controlled mechanobiological experiment. We applied cyclic mechanical loading after implant insertion in tail vertebrae of oestrogen depleted mice and we monitored peri-implant bone response by in vivo micro-CT. Experimental data were combined with micro-finite element simulations to estimate local tissue strains in (re)modelling locations. We demonstrated that a substantial increase in bone mass around the implant could be obtained by loading the entire bone. This augmentation could be attributed to a large reduction in bone resorption rather than to an increase in bone formation. We also showed that following implantation, mechanical regulation of bone (re)modelling was transiently lost. Our findings should help to clarify the role of mechanical stimulation on the maintenance of peri-implant bone mass.
Collapse
Affiliation(s)
- Zihui Li
- 1 Institute for Biomechanics, ETH Zurich , Zurich , Switzerland
| | - Duncan Betts
- 1 Institute for Biomechanics, ETH Zurich , Zurich , Switzerland
| | - Gisela Kuhn
- 1 Institute for Biomechanics, ETH Zurich , Zurich , Switzerland
| | | | - Ralph Müller
- 1 Institute for Biomechanics, ETH Zurich , Zurich , Switzerland
| | - Davide Ruffoni
- 1 Institute for Biomechanics, ETH Zurich , Zurich , Switzerland.,3 Mechanics of Biological and Bioinspired Materials Laboratory, Department of Aerospace and Mechanical Engineering, University of Liège , Liège , Belgium
| |
Collapse
|
34
|
Stalidzans E, Zanin M, Tieri P, Castiglione F, Polster A, Scheiner S, Pahle J, Stres B, List M, Baumbach J, Lautizi M, Van Steen K, Schmidt HH. Mechanistic Modeling and Multiscale Applications for Precision Medicine: Theory and Practice. NETWORK AND SYSTEMS MEDICINE 2020. [DOI: 10.1089/nsm.2020.0002] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
Affiliation(s)
- Egils Stalidzans
- Computational Systems Biology Group, University of Latvia, Riga, Latvia
- Latvian Biomedical Reasearch and Study Centre, Riga, Latvia
| | - Massimiliano Zanin
- Centro de Tecnología Biomédica, Universidad Politécnica de Madrid, Pozuelo de Alarcón, Spain
| | - Paolo Tieri
- CNR National Research Council, IAC Institute for Applied Computing, Rome, Italy
| | - Filippo Castiglione
- CNR National Research Council, IAC Institute for Applied Computing, Rome, Italy
| | | | - Stefan Scheiner
- Institute for Mechanics of Materials and Structures, Vienna University of Technology, Vienna, Austria
| | - Jürgen Pahle
- BioQuant, Heidelberg University, Heidelberg, Germany
| | - Blaž Stres
- Department of Animal Science, University of Ljubljana, Ljubljana, Slovenia
- Faculty of Civil and Geodetic Engineering, University of Ljubljana, Ljubljana, Slovenia
- Department of Automation, Biocybernetics and Robotics, Jozef Stefan Institute, Ljubljana, Slovenia
| | - Markus List
- Big Data in BioMedicine Research Group, Chair of Experimental Bioinformatics, TUM School of Weihenstephan, Technical University of Munich, Freising, Germany
| | - Jan Baumbach
- Chair of Experimental Bioinformatics, TUM School of Weihenstephan, Technical University of Munich, Freising, Germany
| | - Manuela Lautizi
- Computational Systems Medicine Research Group, Chair of Experimental Bioinformatics, TUM School of Weihenstephan, Technical University of Munich, Freising, Germany
| | - Kristel Van Steen
- BIO-Systems Genetics, GIGA-R, University of Liège, Liège, Belgium
- BIO3—Systems Medicine, Department of Human Genetics, KU Leuven, Leuven, Belgium
| | - Harald H.H.W. Schmidt
- Department of Pharmacology and Personalised Medicine, Faculty of Health, Medicine and Life Science, Maastricht University, Maastricht, The Netherlands
| |
Collapse
|
35
|
Zhu S, Bennett S, Kuek V, Xiang C, Xu H, Rosen V, Xu J. Endothelial cells produce angiocrine factors to regulate bone and cartilage via versatile mechanisms. Am J Cancer Res 2020; 10:5957-5965. [PMID: 32483430 PMCID: PMC7255007 DOI: 10.7150/thno.45422] [Citation(s) in RCA: 46] [Impact Index Per Article: 11.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/28/2020] [Accepted: 04/14/2020] [Indexed: 02/06/2023] Open
Abstract
Blood vessels are conduits distributed throughout the body, supporting tissue growth and homeostasis by the transport of cells, oxygen and nutrients. Endothelial cells (ECs) form the linings of the blood vessels, and together with pericytes, are essential for organ development and tissue homeostasis through producing paracrine signalling molecules, called angiocrine factors. In the skeletal system, ECs - derived angiocrine factors, combined with bone cells-released angiogenic factors, orchestrate intercellular crosstalk of the bone microenvironment, and the coupling of angiogenesis-to-osteogenesis. Whilst the involvement of angiogenic factors and the blood vessels of the skeleton is relatively well established, the impact of ECs -derived angiocrine factors on bone and cartilage homeostasis is gradually emerging. In this review, we survey ECs - derived angiocrine factors, which are released by endothelial cells of the local microenvironment and by distal organs, and act specifically as regulators of skeletal growth and homeostasis. These may potentially include angiocrine factors with osteogenic property, such as Hedgehog, Notch, WNT, bone morphogenetic protein (BMP), fibroblast growth factor (FGF), insulin-like growth factor (IGF), and platelet-derived growth factor (PDGF). Understanding the versatile mechanisms by which ECs-derived angiocrine factors orchestrate bone and cartilage homeostasis, and pathogenesis, is an important step towards the development of therapeutic potential for skeletal diseases.
Collapse
|
36
|
Martin M, Sansalone V, Cooper DML, Forwood MR, Pivonka P. Assessment of romosozumab efficacy in the treatment of postmenopausal osteoporosis: Results from a mechanistic PK-PD mechanostat model of bone remodeling. Bone 2020; 133:115223. [PMID: 31935526 DOI: 10.1016/j.bone.2020.115223] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 09/03/2019] [Revised: 12/03/2019] [Accepted: 01/03/2020] [Indexed: 01/01/2023]
Abstract
This paper introduces a theoretical framework for the study of the efficacy of romosozumab, a humanized monoclonal antibody targeting sclerostin for the treatment of osteoporosis. We developed a comprehensive mechanistic pharmacokinetic-pharmacodynamic (PK-PD) model of the effect of drug treatment on bone remodeling in postmenopausal osteoporosis (PMO). We utilized a one-compartment PK model to represent subcutaneous injections of romosozumab and subsequent absorption into serum. The PD model is based on a recently-developed bone cell population model describing the bone remodeling process at the tissue scale. The latter accounts for mechanical feedback by incorporating nitric oxide (NO) and sclerostin (Scl) as biochemical feedback molecules. Utilizing a competitive binding model, where Wnt and Scl compete for binding to LRP5/6, allows to regulate anabolic bone remodeling responses. Here, we extended this model with respect to romosozumab binding to sclerostin. For the currently approved monthly injections of 210 mg, the model predicted a 6.59%, 10.38% and 15.25% increase in BMD at the lumbar spine after 6, 12 and 24 months, respectively. These results are in good agreement with the data reported in the literature. Our model is also able to distinguish the bone-site specific drug effects. For instance, at the femoral neck, our model predicts a BMD increase of 3.85% after 12 months of 210 mg injections, which is consistent with literature observations. Finally, our simulations indicate rapid bone loss after treatment discontinuation, indicating that some additional interventions such as use of bisphosphonates are required to maintain bone mass.
Collapse
Affiliation(s)
- Madge Martin
- School of Chemistry, Physics and Mechanical Engineering, Queensland University of Technology, 2 George St, Brisbane, QLD 4000, Australia; Laboratoire Modélisation et Simulation Multi-Echelle (MSME), UMR CNRS 8208, Université Paris-Est Créteil, 61 avenue du Général de Gaulle, Créteil 94010, France.
| | - Vittorio Sansalone
- Laboratoire Modélisation et Simulation Multi-Echelle (MSME), UMR CNRS 8208, Université Paris-Est Créteil, 61 avenue du Général de Gaulle, Créteil 94010, France
| | - David M L Cooper
- Department of Anatomy, Physiology and Pharmacology, University of Saskatchewan, 107 Wiggins Road, Saskatoon, SK, Canada
| | - Mark R Forwood
- School of Medical Science, Griffith University, Gold Goast, QLD 4222, Australia
| | - Peter Pivonka
- School of Chemistry, Physics and Mechanical Engineering, Queensland University of Technology, 2 George St, Brisbane, QLD 4000, Australia
| |
Collapse
|
37
|
Bahia M, Hecke M, Mercuri E, Pinheiro M. A bone remodeling model governed by cellular micromechanics and physiologically based pharmacokinetics. J Mech Behav Biomed Mater 2020; 104:103657. [DOI: 10.1016/j.jmbbm.2020.103657] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/23/2019] [Revised: 01/11/2020] [Accepted: 01/23/2020] [Indexed: 11/29/2022]
|
38
|
Kameo Y, Miya Y, Hayashi M, Nakashima T, Adachi T. In silico experiments of bone remodeling explore metabolic diseases and their drug treatment. SCIENCE ADVANCES 2020; 6:eaax0938. [PMID: 32181336 PMCID: PMC7060067 DOI: 10.1126/sciadv.aax0938] [Citation(s) in RCA: 30] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/18/2019] [Accepted: 12/13/2019] [Indexed: 05/05/2023]
Abstract
Bone structure and function are maintained by well-regulated bone metabolism and remodeling. Although the underlying molecular and cellular mechanisms are now being understood, physiological and pathological states of bone are still difficult to predict due to the complexity of intercellular signaling. We have now developed a novel in silico experimental platform, V-Bone, to integratively explore bone remodeling by linking complex microscopic molecular/cellular interactions to macroscopic tissue/organ adaptations. Mechano-biochemical couplings modeled in V-Bone relate bone adaptation to mechanical loading and reproduce metabolic bone diseases such as osteoporosis and osteopetrosis. V-Bone also enables in silico perturbation on a specific signaling molecule to observe bone metabolic dynamics over time. We also demonstrate that this platform provides a powerful way to predict in silico therapeutic effects of drugs against metabolic bone diseases. We anticipate that these in silico experiments will substantially accelerate research into bone metabolism and remodeling.
Collapse
Affiliation(s)
- Y. Kameo
- Department of Biosystems Science, Institute for Frontier Life and Medical Sciences, Kyoto University, Kyoto, Japan
- Department of Micro Engineering, Graduate School of Engineering, Kyoto University, Kyoto, Japan
- Department of Mammalian Regulatory Network, Graduate School of Biostudies, Kyoto University, Kyoto, Japan
| | - Y. Miya
- Department of Micro Engineering, Graduate School of Engineering, Kyoto University, Kyoto, Japan
| | - M. Hayashi
- Department of Cell Signaling, Graduate School of Medical and Dental Sciences, Tokyo Medical and Dental University, Tokyo, Japan
| | - T. Nakashima
- Department of Cell Signaling, Graduate School of Medical and Dental Sciences, Tokyo Medical and Dental University, Tokyo, Japan
- Core Research for Evolutional Science and Technology, Japan Agency for Medical Research and Development, Tokyo, Japan
| | - T. Adachi
- Department of Biosystems Science, Institute for Frontier Life and Medical Sciences, Kyoto University, Kyoto, Japan
- Department of Micro Engineering, Graduate School of Engineering, Kyoto University, Kyoto, Japan
- Department of Mammalian Regulatory Network, Graduate School of Biostudies, Kyoto University, Kyoto, Japan
- Core Research for Evolutional Science and Technology, Japan Agency for Medical Research and Development, Tokyo, Japan
- Corresponding author.
| |
Collapse
|
39
|
Study of the combined effects of PTH treatment and mechanical loading in postmenopausal osteoporosis using a new mechanistic PK-PD model. Biomech Model Mechanobiol 2020; 19:1765-1780. [PMID: 32100180 DOI: 10.1007/s10237-020-01307-6] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/10/2019] [Accepted: 02/07/2020] [Indexed: 02/02/2023]
Abstract
One of only a few approved and available anabolic treatments for severe osteoporosis is daily injections of PTH (1-34). This drug has a specific dual action which can act either anabolically or catabolically depending on the type of administration, i.e. intermittent or continuous, respectively. In this paper, we present a mechanistic pharmacokinetic-pharmacodynamic model of the action of PTH in postmenopausal osteoporosis. This model accounts for anabolic and catabolic activities in bone remodelling under intermittent and continuous administration of PTH. The model predicts evolution of common bone biomarkers and bone volume fraction (BV/TV) over time. We compared the relative changes in BV/TV resulting from a daily injection of 20 [Formula: see text]g of PTH with experimental data from the literature. Simulation results indicate a site-specific bone gain of 8.66[Formula: see text] (9.4 ± 1.13[Formula: see text]) at the lumbar spine and 3.14[Formula: see text] (2.82 ± 0.72[Formula: see text]) at the femoral neck. Bone gain depends nonlinearly on the administered dose, being, respectively, 0.68[Formula: see text], 3.4[Formula: see text] and 6.16[Formula: see text] for a 10, 20 and 40 [Formula: see text]g PTH dose at the FN over 2 years. Simulations were performed also taking into account a bone mechanical disuse to reproduce elderly frail subjects. The results show that mechanical disuse ablates the effects of PTH and leads to a 1.08% reduction of bone gain at the FN over a 2-year treatment period for the 20 [Formula: see text]g of PTH. The developed model can simulate a range of pathological conditions and treatments in bones including different PTH doses, different mechanical loading environments and combinations. Consequently, the model can be used for testing and generating hypotheses related to synergistic action between PTH treatment and physical activity.
Collapse
|
40
|
Alcântara ACS, Assis I, Prada D, Mehle K, Schwan S, Costa-Paiva L, Skaf MS, Wrobel LC, Sollero P. Patient-Specific Bone Multiscale Modelling, Fracture Simulation and Risk Analysis-A Survey. MATERIALS (BASEL, SWITZERLAND) 2019; 13:E106. [PMID: 31878356 PMCID: PMC6981613 DOI: 10.3390/ma13010106] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 11/13/2019] [Revised: 12/16/2019] [Accepted: 12/17/2019] [Indexed: 12/26/2022]
Abstract
This paper provides a starting point for researchers and practitioners from biology, medicine, physics and engineering who can benefit from an up-to-date literature survey on patient-specific bone fracture modelling, simulation and risk analysis. This survey hints at a framework for devising realistic patient-specific bone fracture simulations. This paper has 18 sections: Section 1 presents the main interested parties; Section 2 explains the organzation of the text; Section 3 motivates further work on patient-specific bone fracture simulation; Section 4 motivates this survey; Section 5 concerns the collection of bibliographical references; Section 6 motivates the physico-mathematical approach to bone fracture; Section 7 presents the modelling of bone as a continuum; Section 8 categorizes the surveyed literature into a continuum mechanics framework; Section 9 concerns the computational modelling of bone geometry; Section 10 concerns the estimation of bone mechanical properties; Section 11 concerns the selection of boundary conditions representative of bone trauma; Section 12 concerns bone fracture simulation; Section 13 presents the multiscale structure of bone; Section 14 concerns the multiscale mathematical modelling of bone; Section 15 concerns the experimental validation of bone fracture simulations; Section 16 concerns bone fracture risk assessment. Lastly, glossaries for symbols, acronyms, and physico-mathematical terms are provided.
Collapse
Affiliation(s)
- Amadeus C. S. Alcântara
- Department of Computational Mechanics, School of Mechanical Engineering, University of Campinas—UNICAMP, Campinas, Sao Paulo 13083-860, Brazil; (A.C.S.A.); (D.P.)
| | - Israel Assis
- Department of Integrated Systems, School of Mechanical Engineering, University of Campinas—UNICAMP, Campinas, Sao Paulo 13083-860, Brazil;
| | - Daniel Prada
- Department of Computational Mechanics, School of Mechanical Engineering, University of Campinas—UNICAMP, Campinas, Sao Paulo 13083-860, Brazil; (A.C.S.A.); (D.P.)
| | - Konrad Mehle
- Department of Engineering and Natural Sciences, University of Applied Sciences Merseburg, 06217 Merseburg, Germany;
| | - Stefan Schwan
- Fraunhofer Institute for Microstructure of Materials and Systems IMWS, 06120 Halle/Saale, Germany;
| | - Lúcia Costa-Paiva
- Department of Obstetrics and Gynecology, School of Medical Sciences, University of Campinas—UNICAMP, Campinas, Sao Paulo 13083-887, Brazil;
| | - Munir S. Skaf
- Institute of Chemistry and Center for Computing in Engineering and Sciences, University of Campinas—UNICAMP, Campinas, Sao Paulo 13083-860, Brazil;
| | - Luiz C. Wrobel
- Institute of Materials and Manufacturing, Brunel University London, Uxbridge UB8 3PH, UK;
- Department of Civil and Environmental Engineering, Pontifical Catholic University of Rio de Janeiro, Rio de Janeiro 22451-900, Brazil
| | - Paulo Sollero
- Department of Computational Mechanics, School of Mechanical Engineering, University of Campinas—UNICAMP, Campinas, Sao Paulo 13083-860, Brazil; (A.C.S.A.); (D.P.)
| |
Collapse
|
41
|
Yu W, Wu X, Cen H, Guo Y, Li C, Wang Y, Qin Y, Chen W. Study on the biomechanical responses of the loaded bone in macroscale and mesoscale by multiscale poroelastic FE analysis. Biomed Eng Online 2019; 18:122. [PMID: 31870380 PMCID: PMC6929473 DOI: 10.1186/s12938-019-0741-3] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/09/2019] [Accepted: 12/10/2019] [Indexed: 11/10/2022] Open
Abstract
BACKGROUND Bone is a hierarchically structured composite material, and different hierarchical levels exhibit diverse material properties and functions. The stress and strain distribution and fluid flow in bone play an important role in the realization of mechanotransduction and bone remodeling. METHODS To investigate the mechanotransduction and fluid behaviors in loaded bone, a multiscale method was developed. Based on poroelastic theory, we established the theoretical and FE model of a segment bone to provide basis for researching more complex bone model. The COMSOL Multiphysics software was used to establish different scales of bone models, and the properties of mechanical and fluid behaviors in each scale were investigated. RESULTS FE results correlated very well with analytical in macroscopic scale, and the results for the mesoscopic models were about less than 2% different compared to that in the macro-mesoscale models, verifying the correctness of the modeling. In macro-mesoscale, results demonstrated that variations in fluid pressure (FP), fluid velocity (FV), von Mises stress (VMS), and maximum principal strain (MPS) in the position of endosteum, periosteum, osteon, and interstitial bone and these variations can be considerable (up to 10, 8, 4 and 3.5 times difference in maximum FP, FV, VMS, and MPS between the highest and the lowest regions, respectively). With the changing of Young's modulus (E) in each osteon lamella, the strain and stress concentration occurred in different positions and given rise to microscale spatial variations in the fluid pressure field. The heterogeneous distribution of lacunar-canalicular permeability (klcp) in each osteon lamella had various influence on the FP and FV, but had little effect on VMS and MPS. CONCLUSION Based on the idealized model presented in this article, the presence of endosteum and periosteum has an important influence on the fluid flow in bone. With the hypothetical parameter values in osteon lamellae, the bone material parameters have effect on the propagation of stress and fluid flow in bone. The model can also incorporate alternative material parameters obtained from different individuals. The suggested method is expected to provide dependable biological information for better understanding the bone mechanotransduction and signal transduction.
Collapse
Affiliation(s)
- WeiLun Yu
- College of Biomedical Engineering, Shanxi Key Lab. of Material Strength, College of Biomedical Engineering & Structural Impact, Taiyuan University of Technology, Taiyuan, 030024, Shanxi, China
| | - XiaoGang Wu
- College of Biomedical Engineering, Shanxi Key Lab. of Material Strength, College of Biomedical Engineering & Structural Impact, Taiyuan University of Technology, Taiyuan, 030024, Shanxi, China.
| | - HaiPeng Cen
- Biological Science and Medical Engineering, Beihang University, Beijing, China
| | - Yuan Guo
- College of Biomedical Engineering, Shanxi Key Lab. of Material Strength, College of Biomedical Engineering & Structural Impact, Taiyuan University of Technology, Taiyuan, 030024, Shanxi, China
| | - ChaoXin Li
- College of Biomedical Engineering, Shanxi Key Lab. of Material Strength, College of Biomedical Engineering & Structural Impact, Taiyuan University of Technology, Taiyuan, 030024, Shanxi, China
| | - YanQin Wang
- College of Biomedical Engineering, Shanxi Key Lab. of Material Strength, College of Biomedical Engineering & Structural Impact, Taiyuan University of Technology, Taiyuan, 030024, Shanxi, China
| | - YiXian Qin
- Orthopaedic Bioengineering Research Laboratory, Department of Biomedical Engineering, Stony Brook University, Stony Brook, NY, USA
| | - WeiYi Chen
- College of Biomedical Engineering, Shanxi Key Lab. of Material Strength, College of Biomedical Engineering & Structural Impact, Taiyuan University of Technology, Taiyuan, 030024, Shanxi, China.
| |
Collapse
|
42
|
Computational model of the dual action of PTH - Application to a rat model of osteoporosis. J Theor Biol 2019; 473:67-79. [PMID: 31009612 DOI: 10.1016/j.jtbi.2019.04.020] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/17/2018] [Revised: 04/18/2019] [Accepted: 04/19/2019] [Indexed: 11/22/2022]
Abstract
This paper presents a pharmacokinetic/pharmacodynamic (PK/PD) model of the action of PTH(1-34) on bone modelling and remodelling, developed for quantitatively investigating the dose- and administration pattern-dependency of the bone tissue response to this drug. Firstly, a PK model of PTH(1-34) was developed, accounting for administration via subcutaneous injections. Subsequently, the PK model was coupled to a (mechanistic) bone cell population model of bone modelling and remodelling, taking into account the effects of PTH(1-34) on the differentiation of lining cells into active osteoblasts, on the apoptosis of active osteoblasts, and on proliferation of osteoblast precursors, as well as on the key regulatory pathways of bone cell activities. Numerical simulations show that the coupled PK/PD model is able to distinguish between continuous and intermittent administration patterns of PTH(1-34), in terms of yielding both catabolic bone responses (if drug administration is carried out continuously) and anabolic bone responses (if drug administration is carried out intermittently). The model also features a non-linear relation between bone gain and drug dose (as known from experiments); doubling the dose from 80 μg/kg/day to 160 μg/kg/day induced a 1.3-fold increase of the bone volume-to-total volume ratio. Furthermore, the model presented in this paper confirmed that bone modelling represents an essential mechanism of the anabolic response of bone to PTH(1-34) administration in rat models, and that the large amount of bone formation observed in such models cannot be explained via remodelling alone.
Collapse
|
43
|
Della Corte A, Giorgio I, Scerrato D. A review of recent developments in mathematical modeling of bone remodeling. Proc Inst Mech Eng H 2019; 234:273-281. [DOI: 10.1177/0954411919857599] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
Abstract
In this article, we summarize the developments in the mathematical modeling of the mechanics of bone and related biological phenomena. We will devote special attention to the results of the last 10–15 years, although we will cover some relevant classical work to better frame the more recent researches. We will propose a division of the literature based on the main aim of the model (mechanical/biomathematical) and the type of biological phenomena considered (stimulus, growth, cell population dynamics). Finally, we will suggest some possible directions for future investigations.
Collapse
Affiliation(s)
- Alessandro Della Corte
- International Research Center for the Mathematics and Mechanics of Complex Systems, University of L’Aquila, L’Aquila, Italy
| | - Ivan Giorgio
- International Research Center for the Mathematics and Mechanics of Complex Systems, University of L’Aquila, L’Aquila, Italy
- Department of Structural and Geotechnical Engineering, SAPIENZA Università di Roma, Rome, Italy
| | - Daria Scerrato
- International Research Center for the Mathematics and Mechanics of Complex Systems, University of L’Aquila, L’Aquila, Italy
| |
Collapse
|
44
|
Mechanobiological osteocyte feedback drives mechanostat regulation of bone in a multiscale computational model. Biomech Model Mechanobiol 2019; 18:1475-1496. [DOI: 10.1007/s10237-019-01158-w] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/29/2018] [Accepted: 04/23/2019] [Indexed: 10/26/2022]
|
45
|
Wang X, Das R, Fernandez J. A surrogate mechanostatistical microstructural model to inform whole hip cortical bone remodelling. INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING 2019; 35:e3183. [PMID: 30701692 DOI: 10.1002/cnm.3183] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/09/2018] [Revised: 01/05/2019] [Accepted: 01/20/2019] [Indexed: 06/09/2023]
Abstract
Multiscale models of bone remodelling which incorporate clinically important features at a wide range of spatial scales are computationally intractable. In addition, there is a large gap in existing bone remodelling literature which deals with Haversian-level cortical bone remodelling behaviour and its effect on whole organ bone strength. In this study, we present an integration between continuum mechanics and surrogate modelling as an efficient circumvention to the problem of intractability, linked with an intricate 3D cortical bone remodelling algorithm which captures realistic microanatomical features at the Haversian level. The surrogate model predicted clear differentiation in long-term bone strength changes in clinically significant cases with extremely high computational efficiency.
Collapse
Affiliation(s)
- Xiaoming Wang
- Auckland Bioengineering Institute, The University of Auckland, Auckland, New Zealand
| | - Raj Das
- Department of Mechanical Engineering, The University of Auckland, Auckland, New Zealand
| | - Justin Fernandez
- Auckland Bioengineering Institute, The University of Auckland, Auckland, New Zealand
- Department of Engineering Science, The University of Auckland, Auckland, New Zealand
| |
Collapse
|
46
|
Burova I, Wall I, Shipley RJ. Mathematical and computational models for bone tissue engineering in bioreactor systems. J Tissue Eng 2019; 10:2041731419827922. [PMID: 30834100 PMCID: PMC6391543 DOI: 10.1177/2041731419827922] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/03/2018] [Accepted: 01/01/2019] [Indexed: 01/13/2023] Open
Abstract
Research into cellular engineered bone grafts offers a promising solution to problems associated with the currently used auto- and allografts. Bioreactor systems can facilitate the development of functional cellular bone grafts by augmenting mass transport through media convection and shear flow-induced mechanical stimulation. Developing successful and reproducible protocols for growing bone tissue in vitro is dependent on tuning the bioreactor operating conditions to the specific cell type and graft design. This process, largely reliant on a trial-and-error approach, is challenging, time-consuming and expensive. Modelling can streamline the process by providing further insight into the effect of the bioreactor environment on the cell culture, and by identifying a beneficial range of operational settings to stimulate tissue production. Models can explore the impact of changing flow speeds, scaffold properties, and nutrient and growth factor concentrations. Aiming to act as an introductory reference for bone tissue engineers looking to direct their experimental work, this article presents a comprehensive framework of mathematical models on various aspects of bioreactor bone cultures and overviews modelling case studies from literature.
Collapse
Affiliation(s)
- Iva Burova
- Department of Mechanical Engineering, University College London (UCL), London, UK
| | - Ivan Wall
- Aston Medical Research Institute and School of Life & Health Sciences, Aston University, Birmingham, UK
- Institute of Tissue Regeneration Engineering (ITREN), Dankook University, Cheonan, Republic of Korea
| | - Rebecca J Shipley
- Department of Mechanical Engineering, University College London (UCL), London, UK
| |
Collapse
|
47
|
Arias CF, Herrero MA, Echeverri LF, Oleaga GE, López JM. Bone remodeling: A tissue-level process emerging from cell-level molecular algorithms. PLoS One 2018; 13:e0204171. [PMID: 30231062 PMCID: PMC6145577 DOI: 10.1371/journal.pone.0204171] [Citation(s) in RCA: 26] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/16/2018] [Accepted: 09/03/2018] [Indexed: 01/08/2023] Open
Abstract
The human skeleton undergoes constant remodeling throughout the lifetime. Processes occurring on microscopic and molecular scales degrade bone and replace it with new, fully functional tissue. Multiple bone remodeling events occur simultaneously, continuously and independently throughout the body, so that the entire skeleton is completely renewed about every ten years.Bone remodeling is performed by groups of cells called Bone Multicellular Units (BMU). BMUs consist of different cell types, some specialized in the resorption of old bone, others encharged with producing new bone to replace the former. These processes are tightly regulated so that the amount of new bone produced is in perfect equilibrium with that of old bone removed, thus maintaining bone microscopic structure.To date, many regulatory molecules involved in bone remodeling have been identified, but the precise mechanism of BMU operation remains to be fully elucidated. Given the complexity of the signaling pathways already known, one may question whether such complexity is an inherent requirement of the process or whether some subset of the multiple constituents could fulfill the essential role, leaving functional redundancy to serve an alternative safety role. We propose in this work a minimal model of BMU function that involves a limited number of signals able to account for fully functional BMU operation. Our main assumptions were i) at any given time, any cell within a BMU can select only one among a limited choice of decisions, i.e. divide, die, migrate or differentiate, ii) this decision is irreversibly determined by depletion of an appropriate internal inhibitor and iii) the dynamics of any such inhibitor are coupled to that of specific external mediators, such as hormones, cytokines, growth factors. It was thus shown that efficient BMU operation manifests as an emergent process, which results from the individual and collective decisions taken by cells within the BMU unit in the absence of any external planning.
Collapse
Affiliation(s)
- Clemente F. Arias
- Grupo Interdisciplinar de Sistemas Complejos (GISC), Universidad Complutense, 28040 Madrid, Spain
| | - Miguel A. Herrero
- Departamento de Análisis Matemático y Matemática Aplicada, Facultad de Matemáticas, Universidad Complutense, 28040 Madrid, Spain
| | - Luis F. Echeverri
- Instituto de Matemáticas, Facultad de Ciencias Exactas y Naturales, Universidad de Antioquia, 53108 Medellín, Colombia
| | - Gerardo E. Oleaga
- Departamento de Análisis Matemático y Matemática Aplicada, Facultad de Matemáticas, Universidad Complutense, 28040 Madrid, Spain
- Instituto de Matemática Interdisciplinar, Facultad de Matemáticas, Universidad Complutense, 28040 Madrid, Spain
| | - José M. López
- Departamento de Morfología y Biología Celular, Universidad de Oviedo, 33006 Oviedo, Asturias, Spain
| |
Collapse
|
48
|
Budyn E, Gaci N, Sanders S, Bensidhoum M, Schmidt E, Cinquin B, Tauc P, Petite H. Human Stem Cell Derived Osteocytes in Bone-on-Chip. ACTA ACUST UNITED AC 2018. [DOI: 10.1557/adv.2018.278] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2022]
|
49
|
Szlazak K, Vass V, Hasslinger P, Jaroszewicz J, Dejaco A, Idaszek J, Scheiner S, Hellmich C, Swieszkowski W. X-ray physics-based CT-to-composition conversion applied to a tissue engineering scaffold, enabling multiscale simulation of its elastic behavior. MATERIALS SCIENCE & ENGINEERING. C, MATERIALS FOR BIOLOGICAL APPLICATIONS 2017; 95:389-396. [PMID: 30573263 DOI: 10.1016/j.msec.2017.11.044] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/20/2017] [Revised: 09/05/2017] [Accepted: 11/29/2017] [Indexed: 12/13/2022]
Abstract
Nowadays, the assessment of the mechanical competence of tissue engineering scaffolds based on computer simulations is a well-accepted technology. Typically, such simulations are performed by means of the Finite Element (FE) method, with the underlying structural model being created based on micro-computed tomography (microCT). Here, this analysis modality is applied to a new, ternary composite, consisting of PHBV, i.e. poly(3-hydroxybutyrate-co-3-hydroxyvalerate), PLGA, i.e. poly(lactic-co-glycolide), as well as of TCP, i.e. tricalcium phosphate hydrate. The studied scaffold structure is made up by fibers of this new composite material, manufactured by means of the rapid prototyping method. The data collected from microCT is utilized for adequately defining the mechanical properties of the FE model. In particular, the three-dimensional field of grey values is interpreted in terms of the underlying field of attenuation coefficients, taking into account the photon energy employed in microCT imaging, eventually allowing for calculation of the three-dimensionally distributed, voxel-specific composition of the studied material. For the sake of keeping the FE simulations as efficient as possible, groups of voxels are combined into one finite element; the grey value of the latter is obtained by volume averaging. Employing a two-step micromechanical homogenization scheme, the experimentally accessible stiffness of the three constituents (PHBV, PLGA, and TCP) is then, finite element by finite element, upscaled to the composition-dependent stiffness of the composite material. The plausibility and adequacy of the FE model is demonstrated by simulating the effects of uniaxial compression on the scaffold structure, in terms of resulting stress and strain fields, highlighting the importance of the fiber junctions (as they are the mechanically most stressed regions), and that neglecting the material heterogeneity would lead to a potentially significant underestimation of stresses and strains. Finally, a comparison is made of the employed analysis modality of microCT data with a previously pursued, simplified analysis strategy, highlighting the conceptual superiority of the former, and pointing out the application limits of the latter.
Collapse
Affiliation(s)
- Karol Szlazak
- Faculty of Materials Science and Engineering, Warsaw University of Technology, Warsaw, Poland
| | - Viktoria Vass
- Institute for Mechanics of Materials and Structures, TU Wien - Vienna University of Technology, Vienna, Austria
| | - Patricia Hasslinger
- Institute for Mechanics of Materials and Structures, TU Wien - Vienna University of Technology, Vienna, Austria
| | - Jakub Jaroszewicz
- Faculty of Materials Science and Engineering, Warsaw University of Technology, Warsaw, Poland
| | - Alexander Dejaco
- Institute for Mechanics of Materials and Structures, TU Wien - Vienna University of Technology, Vienna, Austria
| | - Joanna Idaszek
- Faculty of Materials Science and Engineering, Warsaw University of Technology, Warsaw, Poland
| | - Stefan Scheiner
- Institute for Mechanics of Materials and Structures, TU Wien - Vienna University of Technology, Vienna, Austria.
| | - Christian Hellmich
- Institute for Mechanics of Materials and Structures, TU Wien - Vienna University of Technology, Vienna, Austria
| | - Wojciech Swieszkowski
- Faculty of Materials Science and Engineering, Warsaw University of Technology, Warsaw, Poland
| |
Collapse
|