A 3D multi-agent-based model for lumen morphogenesis: the role of the biophysical properties of the extracellular matrix

The correct function of many organs depends on proper lumen morphogenesis, which requires the orchestration of both biological and mechanical aspects. However, how these factors coordinate is not yet fully understood. Here, we focus on the development of a mechanistic model for computationally simulating lumen morphogenesis. In particular, we consider the hydrostatic pressure generated by the cells' fluid secretion as the driving force and the density of the extracellular matrix as regulators of the process. For this purpose, we develop a 3D agent-based-model for lumen morphogenesis that includes cells' fluid secretion and the density of the extracellular matrix. Moreover, this computer-based model considers the variation in the biological behavior of cells in response to the mechanical forces that they sense. Then, we study the formation of the lumen under different-mechanical scenarios and conclude that an increase in the matrix density reduces the lumen volume and hinders lumen morphogenesis. Finally, we show that the model successfully predicts normal lumen morphogenesis when the matrix density is physiological and aberrant multilumen formation when the matrix density is excessive.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s00366-022-01654-1.

Confined Cell Migration and Asymmetric Hydraulic Environments to Evaluate The Metastatic Potential of Cancer Cells

Metastasis, a hallmark of cancer development, is also the leading reason for most cancer-related deaths. Furthermore, cancer cells are highly adaptable to microenvironments and can migrate along pre-existing channel-like tracks of anatomical structures. However, more representative three-dimensional models are required to reproduce the heterogeneity of metastatic cell migration in vivo to further understand the metastasis mechanism and develop novel therapeutic strategies against it. Here, we designed and fabricated different microfluidic-based devices that recreate confined migration and diverse environments with asymmetric hydraulic resistances. Our results show different migratory potential between metastatic and nonmetastatic cancer cells in confined environments. Moreover, although nonmetastatic cells have not been tested against barotaxis due to their low migration capacity, metastatic cells present an enhanced preference to migrate through the lowest resistance path, being sensitive to barotaxis. This device, approaching the study of metastasis capability based on confined cell migration and barotactic cell decisions, may pave the way for the implementation of such technology to determine and screen the metastatic potential of certain cancer cells.

Multiscale modeling of bone tissue mechanobiology

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.

Cell biophysical stimuli in lobopodium formation: a computer based approach

Different cell migration modes have been identified in 3D environments, e.g., modes incorporating lamellopodia or blebs. Recently, a new type of cellular migration has been investigated: lobopodia-based migration, which appears only in three-dimensional matrices under certain conditions. The cell creates a protrusion through which the nucleus slips, dividing the cell into two parts (front and rear) with different hydrostatic pressures. In this work, we elucidate the mechanical conditions that favour this type of migration.One of the hypotheses about this type of migration is that it depends on the mechanical properties of the extracellular matrix. That is, lobopodia-based migration is dependent on whether the extracellular matrix is linearly elastic or non-linearly elastic.To determine whether the mechanical properties of the extracellular matrix are crucial in the choice of cell migration mode and which mechanotransduction mechanism the cell might use, we develop a finite element model. From our simulations, we identify two different possible mechanotransduction mechanisms that could regulate the cell to switch from a lobopodial to a lamellipodial migration mode. The first relies on a differential pressure increase inside the cytoplasm while the cell contracts, and the second relies on a change in the fluid flow direction in non-linearly elastic extracellular matrices but not in linearly elastic matrices. The biphasic nature of the cell has been determined to mediate this mechanism and the different behaviours of cells in linearly elastic and non-linearly elastic matrices.

An Evaluation of Surgical Functional Reconstruction of the Foot Using Kinetic and Kinematic Systems: A Case Report

Most pedobarographic studies of microsurgical foot reconstruction have been retrospective. In the present study, we report the results from a prospective pedobarographic study of a patient after microsurgical reconstruction of her foot with a latissimus dorsi flap and a cutaneous paddle, with a 42-month follow-up period. We describe the foot reconstruction plan and the pedobarographic measurements and analyzed its functional outcome. The goal of the present study was to demonstrate that pedobarography could have a role in the treatment of foot reconstruction from a quantitative perspective. The pedobarographic measurements were recorded after the initial coverage surgery and 2 subsequent foot remodeling procedures. A total of 4 pedobarographic measurements and 2 gait analyses were recorded and compared for both the noninvolved foot and the injured foot. Furthermore, the progress of the reconstructed foot was critically evaluated using this method. Both static and dynamic patterns were compared at subsequent follow-up visits after the foot reconstruction. The values and progression of the foot shape, peak foot pressure (kPa), average foot pressure (kPa), total contact surface (cm²), loading time (%), and step time (ms) were recorded. Initially, the pressure distribution of the reconstructed foot showed higher peak values at nonanatomic locations, revealing a greater ulceration risk. Over time, we found an improvement in the shape and values of these factors in the involved foot. To homogenize the pressure distribution and correct the imbalance between the 2 feet, patient-specific insoles were designed and fabricated. In our patient, pedobarography provided an objective, repeatable, and recordable method for the evaluation of the reconstructed foot. Pedobarography can therefore provide valuable insights into the prevention of pressure ulcers and optimization of rehabilitation.

Computational model of mesenchymal migration in 3D under chemotaxis

Cell chemotaxis is an important characteristic of cellular migration, which takes part in crucial aspects of life and development. In this work, we propose a novel in silico model of mesenchymal 3D migration with competing protrusions under a chemotactic gradient. Based on recent experimental observations, we identify three main stages that can regulate mesenchymal chemotaxis: chemosensing, dendritic protrusion dynamics and cell–matrix interactions. Therefore, each of these features is considered as a different module of the main regulatory computational algorithm. The numerical model was particularized for the case of fibroblast chemotaxis under a PDGF-bb gradient. Fibroblasts migration was simulated embedded in two different 3D matrices – collagen and fibrin – and under several PDGF-bb concentrations. Validation of the model results was provided through qualitative and quantitative comparison with in vitro studies. Our numerical predictions of cell trajectories and speeds were within the measured in vitro ranges in both collagen and fibrin matrices. Although in fibrin, the migration speed of fibroblasts is very low, because fibrin is a stiffer and more entangling matrix. Testing PDGF-bb concentrations, we noticed that an increment of this factor produces a speed increment. At 1 ng mL⁻¹ a speed peak is reached after which the migration speed diminishes again. Moreover, we observed that fibrin exerts a dampening behavior on migration, significantly affecting the migration efficiency.

A cell-regulatory mechanism involving feedback between contraction and tissue formation guides wound healing progression

Wound healing is a process driven by cells. The ability of cells to sense mechanical stimuli from the extracellular matrix that surrounds them is used to regulate the forces that cells exert on the tissue. Stresses exerted by cells play a central role in wound contraction and have been broadly modelled. Traditionally, these stresses are assumed to be dependent on variables such as the extracellular matrix and cell or collagen densities. However, we postulate that cells are able to regulate the healing process through a mechanosensing mechanism regulated by the contraction that they exert. We propose that cells adjust the contraction level to determine the tissue functions regulating all main activities, such as proliferation, differentiation and matrix production. Hence, a closed-regulatory feedback loop is proposed between contraction and tissue formation. The model consists of a system of partial differential equations that simulates the evolution of fibroblasts, myofibroblasts, collagen and a generic growth factor, as well as the deformation of the extracellular matrix. This model is able to predict the wound healing outcome without requiring the addition of phenomenological laws to describe the time-dependent contraction evolution. We have reproduced two in vivo experiments to evaluate the predictive capacity of the model, and we conclude that there is feedback between the level of cell contraction and the tissue regenerated in the wound.

Influence of high-frequency cyclical stimulation on the bone fracture-healing process: mathematical and experimental models

Mechanical stimulation affects the evolution of healthy and fractured bone. However, the effect of applying cyclical mechanical stimuli on bone healing has not yet been fully clarified. The aim of the present study was to determine the influence of a high-frequency and low-magnitude cyclical displacement of the fractured fragments on the bone-healing process. This subject is studied experimentally and computationally for a sheep long bone. On the one hand, the mathematical computational study indicates that mechanical stimulation at high frequencies can stimulate and accelerate the process of chondrogenesis and endochondral ossification and consequently the bony union of the fracture. This is probably achieved by the interstitial fluid flow, which can move nutrients and waste from one place to another in the callus. This movement of fluid modifies the mechanical stimulus on the cells attached to the extracellular matrix. On the other hand, the experimental study was carried out using two sheep groups. In the first group, static fixators were implanted, while, in the second one, identical devices were used, but with an additional vibrator. This vibrator allowed a cyclic displacement with low magnitude and high frequency (LMHF) to be applied to the fractured zone every day; the frequency of stimulation was chosen from mechano-biological model predictions. Analysing the results obtained for the control and stimulated groups, we observed improvements in the bone-healing process in the stimulated group. Therefore, in this study, we show the potential of computer mechano-biological models to guide and define better mechanical conditions for experiments in order to improve bone fracture healing. In fact, both experimental and computational studies indicated improvements in the healing process in the LMHF mechanically stimulated fractures. In both studies, these improvements could be associated with the promotion of endochondral ossification and an increase in the rate of cell proliferation and tissue synthesis.

An interspecies computational study on limb lengthening

Distraction osteogenesis is a surgical technique that produces large volumes of new bone by gradually separating two osteotomized bone segments. A previously proposed mechanical-based model that includes the effect of pre-traction stresses (stress level in the gap tissue before each distraction step) during limb lengthening is used here. In the present work, the spatial and temporal patterns of tissue distribution during distraction osteogenesis in different species (sheep, rabbit) and in the human are compared numerically to predict experimental results. Interspecies differential characteristics such as size, distraction protocol, and rate of distraction, among others, are chosen according to experiments. Tissue distributions and reaction forces are then analysed as indicators of the healing pattern. The results obtained are in agreement with experimental findings regarding both tissue distribution and reaction forces. The ability of the model to qualitatively predict the two animal models and the human healing pattern in distraction osteogenesis indicates its potential in understanding the influence of mechanics in this complex process.

Computational comparison of reamed versus unreamed intramedullary tibial nails

We compared, via a computational model, the biomechanical performance of reamed versus unreamed intramedullary tibial nails to treat fractures in three different locations: proximal, mid-diaphyseal, and distal. Two finite element models were analyzed for the two nail types and the three kinds of fractures. Several biomechanical variables were determined: interfragmentary strains in the fracture site, von Mises stresses in nails and bolts, and strain distributions in the tibia and fibula. Although good mechanical stabilization was achieved in all the simulated fractures, the best results were obtained in the proximal fracture for the unreamed nail and in the mid-diaphyseal and distal fractures for the reamed nail. The interlocking bolts, in general, were subjected to higher stresses in the unreamed tibial nail than in the reamed one; thus the former stabilization technique is more likely to fail due to fatigue.

A damage model for the growth plate: Application to the prediction of slipped capital epiphysis

Despite slipped capital femoral epiphysis (SCFE) being one of the most common disorders of the adolescent hip, its early diagnosis is quite difficult. The main objective of this work is to apply an interface damage model to predict the failure of the bone-growth plate-bone interface. This model allows to evaluate the risk of development of SCFE and to investigate the range of mechanical properties of the physis that may cause slippage of the plate. This paper also studies the influence of different geometrical parameters and body weight of the patient on the development of SCFE. We have demonstrated, thanks to the proposed model, that higher physeal sloping and posterior sloping angles are associated to a higher probability of development of SCFE. In a similar way, increasing body weight results in a more probable slippage.

Computational simulation of fracture healing: Influence of interfragmentary movement on the callus growth

Bone fractures heal through a complex process involving several cellular events. This healing process can serve to study factors that control tissue growth and differentiation from mesenchymal stem cells. The mechanical environment at the fracture site is one of the factors influencing the healing process and controls size and differentiation patterns in the newly formed tissue. Mathematical models can be useful to unravel the complex relation between mechanical environment and tissue formation. In this study, we present a mathematical model that predicts tissue growth and differentiation patterns from local mechanical signals. Our aim was to investigate whether mechanical stimuli, through their influence on stem cell proliferation and chondrocyte hypertrophy, predict characteristic features of callus size and geometry. We found that the model predicted several geometric features of fracture calluses. For instance, callus size was predicted to increase with increasing movement. Also, increases in size were predicted to occur through increase in callus diameter but not callus length. These features agree with experimental observations. In addition, spatial and temporal tissue differentiation patterns were in qualitative agreement with well-known experimental results. We therefore conclude that local mechanical signals can probably explain the shape and size of fracture calluses.

A 3D computational simulation of fracture callus formation: influence of the stiffness of the external fixator

The stiffness of the external fixation highly influences the fracture healing pattern. In this work we study this aspect by means of a finite element model of a simple transverse mid-diaphyseal fracture of an ovine metatarsus fixed with a bilateral external fixator. In order to simulate the regenerative process, a previously developed mechanobiological model of bone fracture healing was implemented in three dimensions. This model is able to simulate tissue differentiation, bone regeneration, and callus growth. A physiological load of 500 N was applied and three different stiffnesses of the external fixator were simulated (2300, 1725, and 1150 N/mm). The interfragmentary strain and load sharing mechanism between bone and the external fixator were compared to those recorded in previous experimental works. The effects of the stiffness on the callus shape and tissue distributions in the fracture site were also analyzed. We predicted that a lower stiffness of the fixator delays fracture healing and causes a larger callus, in correspondence to well-documented clinical observations.

Influence of fracture gap size on the pattern of long bone healing: a computational study

Following fractures, bones restore their original structural integrity through a complex process in which several cellular events are involved. Among other factors, this process is highly influenced by the mechanical environment of the fracture site. In this study, we present a mathematical model to simulate the effect of mechanical stimuli on most of the cellular processes that occur during fracture healing, namely proliferation, migration and differentiation. On the basis of these three processes, the model then simulates the evolution of geometry, distributions of cell types and elastic properties inside a healing fracture. The three processes were implemented in a Finite Element code as a combination of three coupled analysis stages: a biphasic, a diffusion and a thermoelastic step. We tested the mechano-biological regulatory model thus created by simulating the healing patterns of fractures with different gap sizes and different mechanical stimuli. The callus geometry, tissue differentiation patterns and fracture stiffness predicted by the model were similar to experimental observations for every analysed situation.

Finite Element Prediction of Proximal Femoral Fracture Patterns Under Different Loads

The main purpose of this work is to discuss the ability of finite element analyses, together with an appropriate anisotropic fracture criterion, to predict the ultimate load and type of fracture in bones and more specifically in the proximal femur. We show here that the use of a three-dimensional anisotropic criterion provides better results than other well-known isotropic criteria. The criterion parameters and the anisotropic elastic properties were defined in terms of the bone tissue microstructure, quantified by the apparent density and the so-called “fabric tensor”, whose spatial distributions were obtained by means of an anisotropic remodeling model able to capture the main features of the internal structure of long bones. In order to check the validity of the results obtained, they have been compared with those of an experimental work that analyzes different types of fractures induced in the proximal femur by a static overload.