Matrix architecture plays a pivotal role in 3D osteoblast migration: The effect of interstitial fluid flow

Osteoblast migration is a crucial process in bone regeneration, which is strongly regulated by interstitial fluid flow. However, the exact role that such flow exerts on osteoblast migration is still unclear. To deepen the understanding of this phenomenon, we cultured human osteoblasts on 3D microfluidic devices under different fluid flow regimes. Our results show that a slow fluid flow rate by itself is not able to alter the 3D migratory patterns of osteoblasts in collagen-based gels but that at higher fluid flow rates (increased flow velocity) may indirectly influence cell movement by altering the collagen microstructure. In fact, we observed that high fluid flow rates (1 µl/min) are able to alter the collagen matrix architecture and to indirectly modulate the migration pattern. However, when these collagen scaffolds were crosslinked with a chemical crosslinker, specifically, transglutaminase II, we did not find significant alterations in the scaffold architecture or in osteoblast movement. Therefore, our data suggest that high interstitial fluid flow rates can regulate osteoblast migration by means of modifying the orientation of collagen fibers. Together, these results highlight the crucial role of the matrix architecture in 3D osteoblast migration. In addition, we show that interstitial fluid flow in conjunction with the matrix architecture regulates the osteoblast morphology in 3D.

Degradation of extracellular matrix regulates osteoblast migration: A microfluidic-based study

Bone regeneration is strongly dependent on the capacity of cells to move in a 3D microenvironment, where a large cascade of signals is activated. To improve the understanding of this complex process and to advance in the knowledge of the role of each specific signal, it is fundamental to analyze the impact of each factor independently. Microfluidic-based cell culture is an appropriate technology to achieve this objective, because it allows recreating realistic 3D local microenvironments by taking into account the extracellular matrix, cells and chemical gradients in an independent or combined scenario. The main aim of this work is to analyze the impact of extracellular matrix properties and growth factor gradients on 3D osteoblast movement, as well as the role of cell matrix degradation. For that, we used collagen-based hydrogels, with and without crosslinkers, under different chemical gradients, and eventually inhibiting metalloproteinases to tweak matrix degradation. We found that osteoblast's 3D migratory patterns were affected by the hydrogel properties and the PDGF-BB gradient, although the strongest regulatory factor was determined by the ability of cells to remodel the matrix.

Force loading explains spatial sensing of ligands by cells

Cells can sense the density and distribution of extracellular matrix (ECM) molecules by means of individual integrin proteins and larger, integrin-containing adhesion complexes within the cell membrane. This spatial sensing drives cellular activity in a variety of normal and pathological contexts. Previous studies of cells on rigid glass surfaces have shown that spatial sensing of ECM ligands takes place at the nanometre scale, with integrin clustering and subsequent formation of focal adhesions impaired when single integrin-ligand bonds are separated by more than a few tens of nanometres. It has thus been suggested that a crosslinking 'adaptor' protein of this size might connect integrins to the actin cytoskeleton, acting as a molecular ruler that senses ligand spacing directly. Here, we develop gels whose rigidity and nanometre-scale distribution of ECM ligands can be controlled and altered. We find that increasing the spacing between ligands promotes the growth of focal adhesions on low-rigidity substrates, but leads to adhesion collapse on more-rigid substrates. Furthermore, disordering the ligand distribution drastically increases adhesion growth, but reduces the rigidity threshold for adhesion collapse. The growth and collapse of focal adhesions are mirrored by, respectively, the nuclear or cytosolic localization of the transcriptional regulator protein YAP. We explain these findings not through direct sensing of ligand spacing, but by using an expanded computational molecular-clutch model, in which individual integrin-ECM bonds-the molecular clutches-respond to force loading by recruiting extra integrins, up to a maximum value. This generates more clutches, redistributing the overall force among them, and reducing the force loading per clutch. At high rigidity and high ligand spacing, maximum recruitment is reached, preventing further force redistribution and leading to adhesion collapse. Measurements of cellular traction forces and actin flow speeds support our model. Our results provide a general framework for how cells sense spatial and physical information at the nanoscale, precisely tuning the range of conditions at which they form adhesions and activate transcriptional regulation.

A hybrid computational model to explore the topological characteristics of epithelial tissues

Epithelial tissues show a particular topology where cells resemble a polygon-like shape, but some biological processes can alter this tissue topology. During cell proliferation, mitotic cell dilation deforms the tissue and modifies the tissue topology. Additionally, cells are reorganized in the epithelial layer and these rearrangements also alter the polygon distribution. We present here a computer-based hybrid framework focused on the simulation of epithelial layer dynamics that combines discrete and continuum numerical models. In this framework, we consider topological and mechanical aspects of the epithelial tissue. Individual cells in the tissue are simulated by an off-lattice agent-based model, which keeps the information of each cell. In addition, we model the cell-cell interaction forces and the cell cycle. Otherwise, we simulate the passive mechanical behaviour of the cell monolayer using a material that approximates the mechanical properties of the cell. This continuum approach is solved by the finite element method, which uses a dynamic mesh generated by the triangulation of cell polygons. Forces generated by cell-cell interaction in the agent-based model are also applied on the finite element mesh. Cell movement in the agent-based model is driven by the displacements obtained from the deformed finite element mesh of the continuum mechanical approach. We successfully compare the results of our simulations with some experiments about the topology of proliferating epithelial tissues in Drosophila. Our framework is able to model the emergent behaviour of the cell monolayer that is due to local cell-cell interactions, which have a direct influence on the dynamics of the epithelial tissue.

Quantifying 3D chemotaxis in microfluidic-based chips with step gradients of collagen hydrogel concentrations

Cell migration is an essential process involved in crucial stages of tissue formation, regeneration or immune function as well as in pathological processes including tumor development or metastasis. During the last few years, the effect of gradients of soluble molecules on cell migration has been widely studied, and complex systems have been used to analyze cell behavior under simultaneous mechano-chemical stimuli. Most of these chemotactic assays have, however, focused on specific substrates in 2D. The aim of the present work is to develop a novel microfluidic-based chip that allows the long-term chemoattractant effect of growth factors (GFs) on 3D cell migration to be studied, while also providing the possibility to analyze the influence of the interface generated between different adjacent hydrogels. Namely, 1.5, 2, 2.5 and 4 mg ml^-1 concentrations of collagen type I were alternatively combined with 5, 10 or 50 ng ml^-1 concentrations of PDGF and VEGF (as a negative control). To achieve this goal, we have designed a new microfluidic device including three adjacent chambers to introduce hydrogels that allow the generation of a collagen concentration step gradient. This versatile and simple platform was tested by using dermal human fibroblasts embedded in 3D collagen matrices. Images taken over a week were processed to quantify the number of cells in each zone. We found, in terms of cell distribution, that the presence of PDGF, especially in small concentrations, was a strong chemoattractant for dermal human fibroblasts across the gels regardless of their collagen concentration and step gradient direction, whereas the effects of VEGF or collagen step gradient concentrations alone were negligible.

A phenomenological cohesive model for the macroscopic simulation of cell-matrix adhesions

Cell adhesion is crucial for cells to not only physically interact with each other but also sense their microenvironment and respond accordingly. In fact, adherent cells can generate physical forces that are transmitted to the surrounding matrix, regulating the formation of cell-matrix adhesions. The main purpose of this work is to develop a computational model to simulate the dynamics of cell-matrix adhesions through a cohesive formulation within the framework of the finite element method and based on the principles of continuum damage mechanics. This model enables the simulation of the mechanical adhesion between cell and extracellular matrix (ECM) as regulated by local multidirectional forces and thus predicts the onset and growth of the adhesion. In addition, this numerical approach allows the simulation of the cell as a whole, as it models the complete mechanical interaction between cell and ECM. As a result, we can investigate and quantify how different mechanical conditions in the cell (e.g., contractile forces, actin cytoskeletal properties) or in the ECM (e.g., stiffness, external forces) can regulate the dynamics of cell-matrix adhesions.

Phenomenological modelling and simulation of cell clusters in 3D cultures

Cell clustering and aggregation are fundamental processes in the development of several tissues and the progression of many diseases. The formation of these aggregates also has a direct impact on the oxygen concentration in their surroundings due to cellular respiration and poor oxygen diffusion through clusters. In this work, we propose a mathematical model that is capable of simulating cell cluster formation in 3D cultures through combining a particle-based and a finite element approach to recreate complex experimental conditions. Cells are modelled considering cell proliferation, cell death and cell-cell mechanical interactions. Additionally, the oxygen concentration profile is calculated through finite element analysis using a reaction-diffusion model that considers cell oxygen consumption and diffusion through the extracellular matrix and the cell clusters. In our model, the local oxygen concentration in the medium determines both cell proliferation and cell death. Numerical predictions are also compared with experimental data from the literature. The simulation results indicate that our model can predict cell clustering, cluster growth and oxygen distribution in 3D cultures. We conclude that the initial cell distribution, cell death and cell proliferation dynamics determine the size and density of clusters. Moreover, these phenomena are directly affected by the oxygen transport in the 3D culture.

Subject-specific musculoskeletal loading of the tibia: Computational load estimation

The systematic development of subject-specific computer models for the analysis of personalized treatments is currently a reality. In fact, many advances have recently been developed for creating virtual finite element-based models. These models accurately recreate subject-specific geometries and material properties from recent techniques based on quantitative image analysis. However, to determine the subject-specific forces, we need a full gait analysis, typically in combination with an inverse dynamics simulation study. In this work, we aim to determine the subject-specific forces from the computer tomography images used to evaluate bone density. In fact, we propose a methodology that combines these images with bone remodelling simulations and artificial neural networks. To test the capability of this novel technique, we quantify the personalized forces for five subject-specific tibias using our technique and a gait analysis. We compare both results, finding that similar vertical loads are estimated by both methods and that the dominant part of the load can be reliably computed. Therefore, we can conclude that the numerical-based technique proposed in this work has great potential for estimating the main forces that define the mechanical behaviour of subject-specific bone.

Biomechanical assessment and clinical analysis of different intramedullary nailing systems for oblique fractures

The aim of this study is to evaluate the fracture union or non-union for a specific patient that presented oblique fractures in tibia and fibula, using a mechanistic-based bone healing model. Normally, this kind of fractures can be treated through an intramedullary nail using two possible configurations that depends on the mechanical stabilisation: static and dynamic. Both cases are simulated under different fracture geometries in order to understand the effect of the mechanical stabilisation on the fracture healing outcome. The results of both simulations are in good agreement with previous clinical experience. From the results, it is demonstrated that the dynamization of the fracture improves healing in comparison with a static or rigid fixation of the fracture. This work shows the versatility and potential of a mechanistic-based bone healing model to predict the final outcome (union, non-union, delayed union) of realistic 3D fractures where even more than one bone is involved.

Modeling the formation of cell-matrix adhesions on a single 3D matrix fiber

Cell-matrix adhesions are crucial in different biological processes like tissue morphogenesis, cell motility, and extracellular matrix remodeling. These interactions that link cell cytoskeleton and matrix fibers are built through protein clutches, generally known as adhesion complexes. The adhesion formation process has been deeply studied in two-dimensional (2D) cases; however, the knowledge is limited for three-dimensional (3D) cases. In this work, we simulate different local extracellular matrix properties in order to unravel the fundamental mechanisms that regulate the formation of cell-matrix adhesions in 3D. We aim to study the mechanical interaction of these biological structures through a three dimensional discrete approach, reproducing the transmission pattern force between the cytoskeleton and a single extracellular matrix fiber. This numerical model provides a discrete analysis of the proteins involved including spatial distribution, interaction between them, and study of the different phenomena, such as protein clutches unbinding or protein unfolding.

Fibroblast Migration in 3D is Controlled by Haptotaxis in a Non-muscle Myosin II-Dependent Manner

Cell migration in 3D is a key process in many physiological and pathological processes. Although valuable knowledge has been accumulated through analysis of various 2D models, some of these insights are not directly applicable to migration in 3D. In this study, we have confined biomimetic hydrogels within microfluidic platforms in the presence of a chemoattractant (platelet-derived growth factor-BB). We have characterized the migratory responses of human fibroblasts within them, particularly focusing on the role of non-muscle myosin II. Our results indicate a prominent role for myosin II in the integration of chemotactic and haptotactic migratory responses of fibroblasts in 3D confined environments.

Inducing chemotactic and haptotactic cues in microfluidic devices for three-dimensional in vitro assays

Microfluidic devices allow for the production of physiologically relevant cellular microenvironments by including biomimetic hydrogels and generating controlled chemical gradients. During transport, the biomolecules interact in distinct ways with the fibrillar networks: as purely diffusive factors in the soluble fluid or bound to the matrix proteins. These two main mechanisms may regulate distinct cell responses in order to guide their directional migration: caused by the substrate-bound chemoattractant gradient (haptotaxis) or by the gradient established within the soluble fluid (chemotaxis). In this work 3D diffusion experiments, in combination with ELISA assays, are performed using microfluidic platforms in order to quantify the distribution of PDGF-BB and TGF-β1 across collagen and fibrin gels. Furthermore, to gain a deeper understanding of the fundamental processes, the experiments are reproduced by computer simulations based on a reaction-diffusion transport model. This model yields an accurate prediction of the experimental results, confirming that diffusion and binding phenomena are established within the microdevice.

Image analysis for the quantitative comparison of stress fibers and focal adhesions

Actin stress fibers (SFs) detect and transmit forces to the extracellular matrix through focal adhesions (FAs), and molecules in this pathway determine cellular behavior. Here, we designed two different computational tools to quantify actin SFs and the distribution of actin cytoskeletal proteins within a normalized cellular morphology. Moreover, a systematic cell response comparison between the control cells and those with impaired actin cytoskeleton polymerization was performed to demonstrate the reliability of the tools. Indeed, a variety of proteins that were present within the string beginning at the focal adhesions (vinculin) up to the actin SFs contraction (non-muscle myosin II (NMMII)) were analyzed. Finally, the software used allows for the quantification of the SFs based on the relative positions of FAs. Therefore, it provides a better insight into the cell mechanics and broadens the knowledge of the nature of SFs.

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.

A time-dependent phenomenological model for cell mechano-sensing

Adherent cells normally apply forces as a generic means of sensing and responding to the mechanical nature of their surrounding environment. How these forces vary as a function of the extracellular rigidity is critical to understanding the regulatory functions that drive important phenomena such as wound healing or muscle contraction. In recognition of this fact, experiments have been conducted to understand cell rigidity-sensing properties under known conditions of the extracellular environment, opening new possibilities for modeling this active behavior. In this work, we provide a physics-based constitutive model taking into account the main structural components of the cell to reproduce its most significant contractile properties such as the traction forces exerted as a function of time and the extracellular stiffness. This model shows how the interplay between the time-dependent response of the acto-myosin contractile system and the elastic response of the cell components determines the mechano-sensing behavior of single cells.

Dynamic mechanisms of cell rigidity sensing: insights from a computational model of actomyosin networks

Cells modulate themselves in response to the surrounding environment like substrate elasticity, exhibiting structural reorganization driven by the contractility of cytoskeleton. The cytoskeleton is the scaffolding structure of eukaryotic cells, playing a central role in many mechanical and biological functions. It is composed of a network of actins, actin cross-linking proteins (ACPs), and molecular motors. The motors generate contractile forces by sliding couples of actin filaments in a polar fashion, and the contractile response of the cytoskeleton network is known to be modulated also by external stimuli, such as substrate stiffness. This implies an important role of actomyosin contractility in the cell mechano-sensing. However, how cells sense matrix stiffness via the contractility remains an open question. Here, we present a 3-D Brownian dynamics computational model of a cross-linked actin network including the dynamics of molecular motors and ACPs. The mechano-sensing properties of this active network are investigated by evaluating contraction and stress in response to different substrate stiffness. Results demonstrate two mechanisms that act to limit internal stress: (i) In stiff substrates, motors walk until they exert their maximum force, leading to a plateau stress that is independent of substrate stiffness, whereas (ii) in soft substrates, motors walk until they become blocked by other motors or ACPs, leading to submaximal stress levels. Therefore, this study provides new insights into the role of molecular motors in the contraction and rigidity sensing of cells.

A mathematical model for cell differentiation, as an evolutionary and regulated process

We introduce an approach which allows one to introduce the concept of cell plasticity into models for tissue regeneration. In contrast to most of the recent models for tissue regeneration, cell differentiation is considered a gradual process, which evolves in time and which is regulated by an arbitrary number of parameters. In the current approach, cell differentiation is modelled by means of a differentiation state variable. Cells are assumed to differentiate into an arbitrary number of cell types. The differentiation path is considered as reversible, unless differentiation has fully completed. Cell differentiation is incorporated into the partial differential equations (PDEs), which model the tissue regeneration process, by means of an advection term in the differentiation state space. This allows one to consider the differentiation path of cells, which is not possible if a reaction-like term is used for differentiation. The boundary conditions, which should be specified for the general PDEs, are derived from the flux of the fully non-differentiated cells and from the irreversibility of the fully completed differentiation process. An application of the proposed model for peri-implant osseointegration is considered. Numerical results are compared with experimental data. Potential lines of further development of the present approach are proposed.

Numerical simulation of bone remodelling around dental implants

Crestal bone loss can result in the failure of dental implants and can be caused, by among other factors, the development of non-physiological mechanical conditions. Bone remodelling (BR) is the physiological process through which bone adapts itself to the mechanical environment. A previously published mathematical model of BR is used in this work to study the homogenized structural evolution of peri-implant bone. This model is used to study the influence of the diameter and length of a dental implant of pure titanium on its long-term stability. The temporal evolution of porosity and microstructural damage of the peri-implant bone are the variables analysed in this study. The results show that damage and porosity increase as the implant length decreases and, more pronouncedly, as its diameter decreases. The increase in damage and porosity levels is localized, as many other studies confirm, at the implant neck due to the stress concentration that is created in that area. The main conclusion of this study is that in implants with a diameter equal to or greater than 3 mm the damage is under control and there is no mechanical failure of the peri-implant bone in the long term.