RESUMEN
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.
RESUMEN
Bone defects represent a challenging clinical problem as they can lead to non-union. In silico models are well suited to study bone regeneration under varying conditions by linking both cellular and systems scales. This paper presents an in silico micro-multiphysics agent-based (micro-MPA) model for bone regeneration following an osteotomy. The model includes vasculature, bone, and immune cells, as well as their interaction with the local environment. The model was calibrated by time-lapsed micro-computed tomography data of femoral osteotomies in C57Bl/6J mice, and the differences between predicted bone volume fractions and the longitudinal in vivo measurements were quantitatively evaluated using root mean square error (RMSE). The model performed well in simulating bone regeneration across the osteotomy gap, with no difference (5.5% RMSE, p = 0.68) between the in silico and in vivo groups for the 5-week healing period - from the inflammatory phase to the remodelling phase - in the volume spanning the osteotomy gap. Overall, the proposed micro-MPA model was able to simulate the influence of the local mechanical environment on bone regeneration, and both this environment and cytokine concentrations were found to be key factors in promoting bone regeneration. Further, the validated model matched clinical observations that larger gap sizes correlate with worse healing outcomes and ultimately simulated non-union. This model could help design and guide future experimental studies in bone repair, by identifying which are the most critical in vivo experiments to perform.
RESUMEN
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.
RESUMEN
Postmenopausal osteoporosis is a disease manifesting in degradation of bone mass and microarchitecture, leading to weakening and increased risk of fracture. Clinical trials are an essential tool for evaluating new treatments and may provide further mechanistic understanding of their effects in vivo. However, the histomorphometry from clinical trials is limited to 2D images and reflects single time points. Biochemical markers of bone turnover give global insight into a drug's action, but not the local dynamics of the bone remodeling process and the cells involved. Additionally, comparative trials necessitate separate treatment groups, meaning only aggregated measures can be compared. In this study, in silico modeling based on histomorphometry and pharmacokinetic data was used to assess the effects of treatment versus control on µCT scans of the same biopsy samples over time, matching the changes in bone volume fraction observed in biopsies from denosumab and placebo groups through year 10 of the FREEDOM Extension trial. In the simulation, treatment decreased osteoclast number, which led to a modest increase in trabecular thickness and osteocyte stress shielding. Long-term bone turnover suppression led to increased RANKL production, followed by a small increase in osteoclast number at the end of the 6-month-dosing interval, especially at the end of the Extension study. Lack of treatment led to a significant loss of bone mass and structure. The study's results show how in silico models can generate predictions of denosumab cellular action over a 10-year period, matching static and dynamic morphometric measures assessed in clinical biopsies. The use of in silico models with clinical trial data can be a method to gain further insight into fundamental bone biology and how treatments can perturb this. With rigorous validation, such models could be used for informing the design of clinical trials, such that the number of participants could be reduced to a minimum to show efficacy. © 2021 The Authors. JBMR Plus published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research.
RESUMEN
Local tissue mechanics play a critical role in cell function, but measuring these properties at cellular length scales in living 3D tissues can present considerable challenges. Here we present thermoresponsive, smart material microgels that can be dispersed or injected into tissues and optically assayed to measure residual tissue elasticity after creep over several weeks. We first develop and characterize the sensors, and demonstrate that internal mechanical profiles of live multicellular spheroids can be mapped at high resolutions to reveal broad ranges of rigidity within the tissues, which vary with subtle differences in spheroid aggregation method. We then show that small sites of unexpectedly high rigidity develop in invasive breast cancer spheroids, and in an in vivo mouse model of breast cancer progression. These focal sites of increased intratumoral rigidity suggest new possibilities for how early mechanical cues that drive cancer cells towards invasion might arise within the evolving tumor microenvironment.
