Search

Shear stress induces directed endothelial cell (EC) migration in blood vessels leading to vessel diameter increase and induction of vascular maturation. Other factors, such as EC elongation and interaction between ECs and non-vascular areas are also important. Computational models have previously been used to study collective cell migration. These models can be used to predict EC migration and its effect on vascular remodelling during embryogenesis. We combined live time-lapse imaging of the remodelling vasculature of the quail embryo yolk sac with flow quantification using a combination of micro-Particle Image Velocimetry and computational fluid dynamics. We then used the flow and remodelling data to inform a model of EC migration during remodelling. To obtain the relation between shear stress and velocity in vitro for EC cells, we developed a flow chamber to assess how confluent sheets of ECs migrate in response to shear stress. Using these data as an input, we developed a multiphase, self-propelled particles (SPP) model where individual agents are driven to migrate based on the level of shear stress while maintaining appropriate spatial relationship to nearby agents. These agents elongate, interact with each other, and with avascular agents at each time-step of the model. We compared predicted vascular shape to real vascular shape after 4 hours from our time-lapse movies and performed sensitivity analysis on the various model parameters. Our model shows that shear stress has the largest effect on the remodelling process. Importantly, however, elongation played an especially important part in remodelling. This model provides a powerful tool to study the input of different biological processes on remodelling., Competing Interests: The authors have declared that no competing interests exist.

Actin protrusion dynamics plays an important role in the regulation of three-dimensional (3D) cell migration. Cells form protrusions that adhere to the surrounding extracellular matrix (ECM), mechanically probe the ECM and contract in order to displace the cell body. This results in cell migration that can be directed by the mechanical anisotropy of the ECM. However, the subcellular processes that regulate protrusion dynamics in 3D cell migration are difficult to investigate experimentally and therefore not well understood. Here, we present a computational model of cell migration through a degradable viscoelastic ECM. This model is a 2D representation of 3D cell migration. The cell is modeled as an active deformable object that captures the viscoelastic behavior of the actin cortex and the subcellular processes underlying 3D cell migration. The ECM is regarded as a viscoelastic material, with or without anisotropy due to fibrillar strain stiffening, and modeled by means of the meshless Lagrangian smoothed particle hydrodynamics (SPH) method. ECM degradation is captured by local fluidization of the material and permits cell migration through the ECM. We demonstrate that changes in ECM stiffness and cell strength affect cell migration and are accompanied by changes in number, lifetime and length of protrusions. Interestingly, directly changing the total protrusion number or the average lifetime or length of protrusions does not affect cell migration. A stochastic variability in protrusion lifetime proves to be enough to explain differences in cell migration velocity. Force-dependent adhesion disassembly does not result in faster migration, but can make migration more efficient. We also demonstrate that when a number of simultaneous protrusions is enforced, the optimal number of simultaneous protrusions is one or two, depending on ECM anisotropy. Together, the model provides non-trivial new insights in the role of protrusions in 3D cell migration and can be a valuable contribution to increase the understanding of 3D cell migration mechanics., Competing Interests: The authors have declared that no competing interests exist.

In order to unravel rapid mechano-chemical feedback mechanisms in sprouting angiogenesis, we combine selective plane illumination microscopy (SPIM) and tailored image registration algorithms - further referred to as SPIM-based displacement microscopy - with an in vitro model of angiogenesis. SPIM successfully tackles the problem of imaging large volumes while upholding the spatial resolution required for the analysis of matrix displacements at a subcellular level. Applied to in vitro angiogenic sprouts, this unique methodological combination relates subcellular activity - minute to second time scale growing and retracting of protrusions - of a multicellular systems to the surrounding matrix deformations with an exceptional temporal resolution of 1 minute for a stack with multiple sprouts simultaneously or every 4 seconds for a single sprout, which is 20 times faster than with a conventional confocal setup. Our study reveals collective but non-synchronised, non-continuous activity of adjacent sprouting cells along with correlations between matrix deformations and protrusion dynamics., Competing Interests: The authors have declared that no competing interests exist.

Traction Force Microscopy (TFM) is a widespread method used to recover cellular tractions from the deformation that they cause in their surrounding substrate. Particle Image Velocimetry (PIV) is commonly used to quantify the substrate's deformations, due to its simplicity and efficiency. However, PIV relies on a block-matching scheme that easily underestimates the deformations. This is especially relevant in the case of large, locally non-uniform deformations as those usually found in the vicinity of a cell's adhesions to the substrate. To overcome these limitations, we formulate the calculation of the deformation of the substrate in TFM as a non-rigid image registration process that warps the image of the unstressed material to match the image of the stressed one. In particular, we propose to use a B-spline -based Free Form Deformation (FFD) algorithm that uses a connected deformable mesh to model a wide range of flexible deformations caused by cellular tractions. Our FFD approach is validated in 3D fields using synthetic (simulated) data as well as with experimental data obtained using isolated endothelial cells lying on a deformable, polyacrylamide substrate. Our results show that FFD outperforms PIV providing a deformation field that allows a better recovery of the magnitude and orientation of tractions. Together, these results demonstrate the added value of the FFD algorithm for improving the accuracy of traction recovery.

Although bone has a unique restorative capacity, i.e., it has the potential to heal scarlessly, the conditions for spontaneous bone healing are not always present, leading to a delayed union or a non-union. In this work, we use an integrative in vivo-in silico approach to investigate the occurrence of non-unions, as well as to design possible treatment strategies thereof. The gap size of the domain geometry of a previously published mathematical model was enlarged in order to study the complex interplay of blood vessel formation, oxygen supply, growth factors and cell proliferation on the final healing outcome in large bone defects. The multiscale oxygen model was not only able to capture the essential aspects of in vivo non-unions, it also assisted in understanding the underlying mechanisms of action, i.e., the delayed vascularization of the central callus region resulted in harsh hypoxic conditions, cell death and finally disrupted bone healing. Inspired by the importance of a timely vascularization, as well as by the limited biological potential of the fracture hematoma, the influence of the host environment on the bone healing process in critical size defects was explored further. Moreover, dependent on the host environment, several treatment strategies were designed and tested for effectiveness. A qualitative correspondence between the predicted outcomes of certain treatment strategies and experimental observations was obtained, clearly illustrating the model's potential. In conclusion, the results of this study demonstrate that due to the complex non-linear dynamics of blood vessel formation, oxygen supply, growth factor production and cell proliferation and the interactions thereof with the host environment, an integrative in silico-in vivo approach is a crucial tool to further unravel the occurrence and treatments of challenging critical sized bone defects.

Bioluminescence imaging assays have become a widely integrated technique to quantify effectiveness of cell-based therapies by monitoring fate and survival of transplanted cells. To date these assays are still largely qualitative and often erroneous due to the complexity and dynamics of local micro-environments (niches) in which the cells reside. Here, we report, using a combined experimental and computational approach, on oxygen that besides being a critical niche component responsible for cellular energy metabolism and cell-fate commitment, also serves a primary role in regulating bioluminescent light kinetics. We demonstrate the potential of an oxygen dependent Michaelis-Menten relation in quantifying intrinsic bioluminescence intensities by resolving cell-associated oxygen gradients from bioluminescent light that is emitted from three-dimensional (3D) cell-seeded hydrogels. Furthermore, the experimental and computational data indicate a strong causal relation of oxygen concentration with emitted bioluminescence intensities. Altogether our approach demonstrates the importance of oxygen to evolve towards quantitative bioluminescence and holds great potential for future microscale measurement of oxygen tension in an easily accessible manner.

Adhesion governs to a large extent the mechanical interaction between a cell and its microenvironment. As initial cell spreading is purely adhesion driven, understanding this phenomenon leads to profound insight in both cell adhesion and cell-substrate interaction. It has been found that across a wide variety of cell types, initial spreading behavior universally follows the same power laws. The simplest cell type providing this scaling of the radius of the spreading area with time are modified red blood cells (RBCs), whose elastic responses are well characterized. Using a mechanistic description of the contact interaction between a cell and its substrate in combination with a deformable RBC model, we are now able to investigate in detail the mechanisms behind this universal power law. The presented model suggests that the initial slope of the spreading curve with time results from a purely geometrical effect facilitated mainly by dissipation upon contact. Later on, the spreading rate decreases due to increasing tension and dissipation in the cell's cortex as the cell spreads more and more. To reproduce this observed initial spreading, no irreversible deformations are required. Since the model created in this effort is extensible to more complex cell types and can cope with arbitrarily shaped, smooth mechanical microenvironments of the cells, it can be useful for a wide range of investigations where forces at the cell boundary play a decisive role.

During endochondral ossification, chondrocyte growth and differentiation is controlled by many local signalling pathways. Due to crosstalks and feedback mechanisms, these interwoven pathways display a network like structure. In this study, a large-scale literature based logical model of the growth plate network was developed. The network is able to capture the different states (resting, proliferating and hypertrophic) that chondrocytes go through as they progress within the growth plate. In a first corroboration step, the effect of mutations in various signalling pathways of the growth plate network was investigated.

The healing of a fracture depends largely on the development of a new blood vessel network (angiogenesis) in the callus. During angiogenesis tip cells lead the developing sprout in response to extracellular signals, amongst which vascular endothelial growth factor (VEGF) is critical. In order to ensure a correct development of the vasculature, the balance between stalk and tip cell phenotypes must be tightly controlled, which is primarily achieved by the Dll4-Notch1 signaling pathway. This study presents a novel multiscale model of osteogenesis and sprouting angiogenesis, incorporating lateral inhibition of endothelial cells (further denoted MOSAIC model) through Dll4-Notch1 signaling, and applies it to fracture healing. The MOSAIC model correctly predicted the bone regeneration process and recapitulated many experimentally observed aspects of tip cell selection: the salt and pepper pattern seen for cell fates, an increased tip cell density due to the loss of Dll4 and an excessive number of tip cells in high VEGF environments. When VEGF concentration was even further increased, the MOSAIC model predicted the absence of a vascular network and fracture healing, thereby leading to a non-union, which is a direct consequence of the mutual inhibition of neighboring cells through Dll4-Notch1 signaling. This result was not retrieved for a more phenomenological model that only considers extracellular signals for tip cell migration, which illustrates the importance of implementing the actual signaling pathway rather than phenomenological rules. Finally, the MOSAIC model demonstrated the importance of a proper criterion for tip cell selection and the need for experimental data to further explore this. In conclusion, this study demonstrates that the MOSAIC model creates enhanced capabilities for investigating the influence of molecular mechanisms on angiogenesis and its relation to bone formation in a more mechanistic way and across different time and spatial scales.

Recently developed atrophic non-union models are a good representation of the clinical situation in which many non-unions develop. Based on previous experimental studies with these atrophic non-union models, it was hypothesized that in order to obtain successful fracture healing, blood vessels, growth factors, and (proliferative) precursor cells all need to be present in the callus at the same time. This study uses a combined in vivo-in silico approach to investigate these different aspects (vasculature, growth factors, cell proliferation). The mathematical model, initially developed for the study of normal fracture healing, is able to capture essential aspects of the in vivo atrophic non-union model despite a number of deviations that are mainly due to simplifications in the in silico model. The mathematical model is subsequently used to test possible treatment strategies for atrophic non-unions (i.e. cell transplant at post-osteotomy, week 3). Preliminary in vivo experiments corroborate the numerical predictions. Finally, the mathematical model is applied to explain experimental observations and identify potentially crucial steps in the treatments and can thereby be used to optimize experimental and clinical studies in this area. This study demonstrates the potential of the combined in silico-in vivo approach and its clinical implications for the early treatment of patients with problematic fractures.