Mechanical cell-matrix feedback explains pairwise and collective endothelial cell behavior in vitro

Substrate stiffness heterogeneities disrupt endothelial barrier integrity in a micropillar model of heterogeneous vascular stiffening

Intimal stiffening has been linked with increased vascular permeability and leukocyte transmigration, hallmarks of atherosclerosis. However, recent evidence indicates age-related intimal stiffening is not uniform but rather characterized by increased point-to-point heterogeneity in subendothelial matrix stiffness, the impact of which is much less understood. To investigate the impact of spatially heterogeneous matrix rigidity on endothelial monolayer integrity, we develop a micropillar model to introduce closely-spaced, step-changes in substrate rigidity and compare endothelial monolayer phenotype to rigidity-matched, uniformly stiff and compliant substrates. We found equivalent disruption of adherens junctions within monolayers on step-rigidity and uniformly stiff substrates relative to uniformly compliant substrates. Similarly, monolayers cultured on step-rigidity substrates exhibited equivalent percentages of leukocyte transmigration to monolayers on rigidity-matched, uniformly stiff substrates. Adherens junction tension and focal adhesion density, but not size, increased within monolayers on step-rigidity and uniformly stiff substrates compared to more compliant substrates suggesting that elevated tension is disrupting adherens junction integrity. Leukocyte transmigration frequency and time, focal adhesion size, and focal adhesion density did not differ between stiff and compliant sub-regions of step-rigidity substrates. Overall, our results suggest that endothelial monolayers exposed to mechanically heterogeneous substrates adopt the phenotype associated with the stiffer matrix, indicating that spatial heterogeneities in intimal stiffness observed with age could disrupt endothelial barrier integrity and contribute to atherogenesis.

Data driven mathematical modeling reveals the dynamic mechanism of MSC-induced neovascularization

Coculture of mesenchymal stem cells (MSCs) and vascular endothelial cells (ECs) in vitro leads to the formation of a capillary-like reticular structure by ECs, which has great potential as a better substitute for artificial blood vessels in terms of stability and functionality. To investigate the mechanisms of the early neovascularization induced by MSCs, we analyzed the kinematic features of the motion of ECs and concluded that the dynamic interaction between cells and the extracellular matrix would reveal the capillary-like structure formation. Based on this hypothesis, we proposed a mathematical model to simulate the vascular-like migration pattern of ECs in silico, which was confirmed by in vitro studies. These in vitro studies validated that the dynamic secretion and degradation of collagen I is the critical factor for capillary structure formation. The model proposed based on cell tracking, single cell sequencing, and mathematical simulation provides a better understanding of the neovascularization process induced by MSCs and a possible simple explanation guiding this important cellular behavior.-Yu, Y., Situ, Q., Jia, W., Li, J., Wu, Q., Lei, J. Data driven mathematical modeling reveals the dynamic mechanism of MSC-induced neovascularization.

Soluble VEGFR1 signaling guides vascular patterns into dense branching morphologies

Vascular patterning is a key process during development and disease. The diffusive decoy receptor sVEGFR1 (sFlt1) is a known regulator of endothelial cell behavior, yet the mechanism by which it controls vascular structure is little understood. We propose computational models to shed light on how vascular patterning is guided by self-organized gradients of the VEGF/sVEGFR1 factors. We demonstrate that a diffusive inhibitor can generate structures with a dense branching morphology in models where the activator elicits directed growth. Inadequate presence of the inhibitor leads to compact growth, while excessive production of the inhibitor blocks expansion and stabilizes existing structures. Model predictions were compared with time-resolved experimental data obtained from endothelial sprout kinetics in fibrin gels. In the presence of inhibitory antibodies against VEGFR1 vascular sprout density increases while the speed of sprout expansion remains unchanged. Thus, the rate of secretion and stability of extracellular sVEGFR1 can modulate vascular sprout density.

Capillary network formation from dispersed endothelial cells: Influence of cell traction, cell adhesion, and extracellular matrix rigidity

The formation of a functional vascular network depends on biological, chemical, and physical processes being extremely well coordinated. Among them, the mechanical properties of the extracellular matrix and cell adhesion are fundamental to achieve a functional network of endothelial cells, able to fully cover a required domain. By the use of a Cellular Potts Model and Finite Element Method it is shown that there exists a range of values of endothelial traction forces, cell-cell adhesion, and matrix rigidities where the network can spontaneously be formed, and its properties are characterized. We obtain the analytical relation that the minimum traction force required for cell network formation must obey. This minimum value for the traction force is approximately independent on the considered cell number and cell-cell adhesion. We quantify how these two parameters influence the morphology of the resulting networks (size and number of meshes).

Modelling skin wound healing angiogenesis: A review

The occurrence of wounds is a main health concern in Western society due to their high frequency and treatment cost. During wound healing, the formation of a functional blood vessel network through angiogenesis is an essential process. Angiogenesis allows the reestablishment of the normal blood flow, the sufficient exchange of oxygen and nutrients and the removal of metabolic waste, necessary for cell proliferation and viability. Mathematical and computational models provide new tools to improve the healing process. In fact, over the last thirty years, in silico models have been continuously formulated to describe the effect of several biological and mechanical factors in angiogenesis during wound healing. Additionally, with different levels of complexity, these models allow coupling the human skin structure, to distinct cell types and growth factors, to study extracellular matrix composition and to understand its deformation. This paper discusses how in silico models, which are more economical and less time-consuming comparatively to laboratory methodologies, can help test new strategies to promote/optimize angiogenesis. The continuum, cell-based and hybrid mathematical models of wound healing angiogenesis are reviewed in the present paper, in order to identify possible improvements. Accordingly, the development of higher dimension models incorporating multiscale analysis at molecular, cellular and tissue level remains a challenge that future models should consider.

A local uPAR-plasmin-TGFβ1 positive feedback loop in a qualitative computational model of angiogenic sprouting explains the in vitro effect of fibrinogen variants

In experimental assays of angiogenesis in three-dimensional fibrin matrices, a temporary scaffold formed during wound healing, the type and composition of fibrin impacts the level of sprouting. More sprouts form on high molecular weight (HMW) than on low molecular weight (LMW) fibrin. It is unclear what mechanisms regulate the number and the positions of the vascular-like structures in cell cultures. To address this question, we propose a mechanistic simulation model of endothelial cell migration and fibrin proteolysis by the plasmin system. The model is a hybrid, cell-based and continuum, computational model based on the cellular Potts model and sets of partial-differential equations. Based on the model results, we propose that a positive feedback mechanism between uPAR, plasmin and transforming growth factor β1 (TGFβ1) selects cells in the monolayer for matrix invasion. Invading cells releases TGFβ1 from the extracellular matrix through plasmin-mediated fibrin degradation. The activated TGFβ1 further stimulates fibrin degradation and keeps proteolysis active as the sprout invades the fibrin matrix. The binding capacity for TGFβ1 of LMW is reduced relative to that of HMW. This leads to reduced activation of proteolysis and, consequently, reduced cell ingrowth in LMW fibrin compared to HMW fibrin. Thus our model predicts that endothelial cells in LMW fibrin matrices compared to HMW matrices show reduced sprouting due to a lower bio-availability of TGFβ1.

Therapies for a range of medical conditions, including cancer, wound healing and diabetic retinopathy can benefit from a better control over the growth of blood vessels. The chemical properties of fibrin, the material that forms scabs in wounds and can also occur in large concentrations in tumors, can regulate the degree of blood vessel growth (angiogenesis). Angiogenesis can be mimicked in cell cultures. These allow us to modulate the chemical properties of fibrin and study the effect on angiogenesis. Fibrin occurs in high molecular weight (HMW) and in low molecular weight (LMW) forms. Interestingly, there is more ingrowth of angiogenic-like structures into HMW than in LMW fibrin, but the mechanisms are poorly understood. To get more insight into these, we constructed a computational model. Using the model, we propose and analyse a hypothetical mechanism for sprouting that could explain the differences in endothelial cell sprouting in LMW and HMW fibrin matrices. Our model suggests that cells digest fibrin, thus creating space for ingrowth. At the same time, digestion frees growth factors bound to fibrin, that activates further secretion of digestive enzymes by the cells. We propose that the resulting positive feedback loop spontaneously selects cells in the endothelial monolayer for ingrowth and helps the blood vessel sprout move deeper into the fibrin. This could be a complementary mechanism to lateral-inhibition by Delta-Notch for the selection of leader cells, also called ‘tip cells’. Our model predicts that endothelial cells in LMW fibrin compared to HMW fibrin show reduced sprouting due to a lower bio-availability of TGFβ1.

Angiogenic Factors produced by Hypoxic Cells are a leading driver of Anastomoses in Sprouting Angiogenesis-a computational study

Angiogenesis - the growth of new blood vessels from a pre-existing vasculature - is key in both physiological processes and on several pathological scenarios such as cancer progression or diabetic retinopathy. For the new vascular networks to be functional, it is required that the growing sprouts merge either with an existing functional mature vessel or with another growing sprout. This process is called anastomosis. We present a systematic 2D and 3D computational study of vessel growth in a tissue to address the capability of angiogenic factor gradients to drive anastomosis formation. We consider that these growth factors are produced only by tissue cells in hypoxia, i.e. until nearby vessels merge and become capable of carrying blood and irrigating their vicinity. We demonstrate that this increased production of angiogenic factors by hypoxic cells is able to promote vessel anastomoses events in both 2D and 3D. The simulations also verify that the morphology of these networks has an increased resilience toward variations in the endothelial cell's proliferation and chemotactic response. The distribution of tissue cells and the concentration of the growth factors they produce are the major factors in determining the final morphology of the network.

A phenomenological model for cell and nucleus deformation during cancer metastasis

Cell migration plays an essential role in cancer metastasis. In cancer invasion through confined spaces, cells must undergo extensive deformation, which is a capability related to their metastatic potentials. Here, we simulate the deformation of the cell and nucleus during invasion through a dense, physiological microenvironment by developing a phenomenological computational model. In our work, cells are attracted by a generic emitting source (e.g., a chemokine or stiffness signal), which is treated by using Green’s Fundamental solutions. We use an IMEX integration method where the linear parts and the nonlinear parts are treated by using an Euler backward scheme and an Euler forward method, respectively. We develop the numerical model for an obstacle-induced deformation in 2D or/and 3D. Considering the uncertainty in cell mobility, stochastic processes are incorporated and uncertainties in the input variables are evaluated using Monte Carlo simulations. This quantitative study aims at estimating the likelihood for invasion and the length of the time interval in which the cell invades the tissue through an obstacle. Subsequently, the two-dimensional cell deformation model is applied to simplified cancer metastasis processes to serve as a model for in vivo or in vitro biomedical experiments.

Eigenstrain as a mechanical set-point of cells

Cell contraction regulates how cells sense their mechanical environment. We sought to identify the set-point of cell contraction, also referred to as tensional homeostasis. In this work, bovine aortic endothelial cells (BAECs), cultured on substrates with different stiffness, were characterized using traction force microscopy (TFM). Numerical models were developed to provide insights into the mechanics of cell-substrate interactions. Cell contraction was modeled as eigenstrain which could induce isometric cell contraction without external forces. The predicted traction stresses matched well with TFM measurements. Furthermore, our numerical model provided cell stress and displacement maps for inspecting the fundamental regulating mechanism of cell mechanosensing. We showed that cell spread area, traction force on a substrate, as well as the average stress of a cell were increased in response to a stiffer substrate. However, the cell average strain, which is cell type-specific, was kept at the same level regardless of the substrate stiffness. This indicated that the cell average strain is the tensional homeostasis that each type of cell tries to maintain. Furthermore, cell contraction in terms of eigenstrain was found to be the same for both BAECs and fibroblast cells in different mechanical environments. This implied a potential mechanical set-point across different cell types. Our results suggest that additional measurements of contractility might be useful for monitoring cell mechanosensing as well as dynamic remodeling of the extracellular matrix (ECM). This work could help to advance the understanding of the cell-ECM relationship, leading to better regenerative strategies.

Computational modeling of three-dimensional ECM-rigidity sensing to guide directed cell migration

Filopodia have a key role in sensing both chemical and mechanical cues in surrounding extracellular matrix (ECM). However, quantitative understanding is still missing in the filopodial mechanosensing of local ECM stiffness, resulting from dynamic interactions between filopodia and the surrounding 3D ECM fibers. Here we present a method for characterizing the stiffness of ECM that is sensed by filopodia based on the theory of elasticity and discrete ECM fiber. We have applied this method to a filopodial mechanosensing model for predicting directed cell migration toward stiffer ECM. This model provides us with a distribution of force and displacement as well as their time rate of changes near the tip of a filopodium when it is bound to the surrounding ECM fibers. Aggregating these effects in each local region of 3D ECM, we express the local ECM stiffness sensed by the cell and explain polarity in the cellular durotaxis mechanism.

Traction Force Microscopy for Noninvasive Imaging of Cell Forces

The forces exerted by cells on their surroundings play an integral role in both physiological processes and disease progression. Traction force microscopy is a noninvasive technique that enables the in vitro imaging and quantification of cell forces. Utilizing expertise from a variety of disciplines, recent developments in traction force microscopy are enhancing the study of cell forces in physiologically relevant model systems, and hold promise for further advancing knowledge in mechanobiology. In this chapter, we discuss the methods, capabilities, and limitations of modern approaches for traction force microscopy, and highlight ongoing efforts and challenges underlying future innovations.

Mechanoreciprocity in cell migration

Cell migration is an adaptive process that depends on and responds to physical and molecular triggers. Moving cells sense and respond to tissue mechanics and induce transient or permanent tissue modifications, including extracellular matrix stiffening, compression and deformation, protein unfolding, proteolytic remodelling and jamming transitions. Here we discuss how the bi-directional relationship of cell-tissue interactions (mechanoreciprocity) allows cells to change position and contributes to single-cell and collective movement, structural and molecular tissue organization, and cell fate decisions.

Matrix Stiffness Enhances VEGFR-2 Internalization, Signaling, and Proliferation in Endothelial Cells

Vascular endothelial growth factor (VEGF) can mediate endothelial cell migration, proliferation, and angiogenesis. During cancer progression, VEGF production is often increased to stimulate the growth of new blood vessels to supply growing tumors with the additional oxygen and nutrients they require. Extracellular matrix stiffening also occurs during tumor progression, however, the crosstalk between tumor mechanics and VEGF signaling remains poorly understood. Here, we show that matrix stiffness heightens downstream endothelial cell response to VEGF by altering VEGF receptor-2 (VEGFR-2) internalization, and this effect is influenced by cell confluency. In sub-confluent endothelial monolayers, VEGFR-2 levels, but not VEGFR-2 phosphorylation, are influenced by matrix rigidity. Interestingly, more compliant matrices correlated with increased expression and clustering of VEGFR-2; however, stiffer matrices induced increased VEGFR-2 internalization. These effects are most likely due to actin-mediated contractility, as inhibiting ROCK on stiff substrates increased VEGFR-2 clustering and decreased internalization. Additionally, increasing matrix stiffness elevates ERK 1/2 phosphorylation, resulting in increased cell proliferation. Moreover, cells on stiff matrices generate more actin stress fibers than on compliant substrates, and the addition of VEGF stimulates an increase in fiber formation regardless of stiffness. In contrast, once endothelial cells reached confluency, stiffness-enhanced VEGF signaling was no longer observed. Together, these data show a complex effect of VEGF and matrix mechanics on VEGF-induced signaling, receptor dynamics, and cell proliferation that is mediated by cell confluency.

A model for cell migration in non-isotropic fibrin networks with an application to pancreatic tumor islets

Cell migration, known as an orchestrated movement of cells, is crucially important for wound healing, tumor growth, immune response as well as other biomedical processes. This paper presents a cell-based model to describe cell migration in non-isotropic fibrin networks around pancreatic tumor islets. This migration is determined by the mechanical strain energy density as well as cytokines-driven chemotaxis. Cell displacement is modeled by solving a large system of ordinary stochastic differential equations where the stochastic parts result from random walk. The stochastic differential equations are solved by the use of the classical Euler–Maruyama method. In this paper, the influence of anisotropic stromal extracellular matrix in pancreatic tumor islets on T-lymphocytes migration in different immune systems is investigated. As a result, tumor peripheral stromal extracellular matrix impedes the immune response of T-lymphocytes through changing direction of their migration.

Virtual cardiac monolayers for electrical wave propagation

The complex structure of cardiac tissue is considered to be one of the main determinants of an arrhythmogenic substrate. This study is aimed at developing the first mathematical model to describe the formation of cardiac tissue, using a joint in silico-in vitro approach. First, we performed experiments under various conditions to carefully characterise the morphology of cardiac tissue in a culture of neonatal rat ventricular cells. We considered two cell types, namely, cardiomyocytes and fibroblasts. Next, we proposed a mathematical model, based on the Glazier-Graner-Hogeweg model, which is widely used in tissue growth studies. The resultant tissue morphology was coupled to the detailed electrophysiological Korhonen-Majumder model for neonatal rat ventricular cardiomyocytes, in order to study wave propagation. The simulated waves had the same anisotropy ratio and wavefront complexity as those in the experiment. Thus, we conclude that our approach allows us to reproduce the morphological and physiological properties of cardiac tissue.

Mechanotransduction Dynamics at the Cell-Matrix Interface

The ability of cells to sense and respond to mechanical cues from the surrounding environment has been implicated as a key regulator of cell differentiation, migration, and proliferation. The extracellular matrix (ECM) is an oft-overlooked component of the interface between cells and their surroundings. Cells assemble soluble ECM proteins into insoluble fibrils with unique mechanical properties that can alter the mechanical cues a cell receives. In this study, we construct a model that predicts the dynamics of cellular traction force generation and subsequent assembly of fibrils of the ECM protein fibronectin (FN). FN fibrils are the primary component in primordial ECM and, as such, FN assembly is a critical component in the cellular mechanical response. The model consists of a network of Hookean springs, each representing an extensible domain within an assembling FN fibril. As actomyosin forces stretch the spring network, simulations predict the resulting traction force and FN fibril formation. The model accurately predicts FN fibril morphometry and demonstrates a mechanism by which FN fibril assembly regulates traction force dynamics in response to mechanical stimuli and varying surrounding substrate stiffness.

Cell Contractility Facilitates Alignment of Cells and Tissues to Static Uniaxial Stretch

During animal development and homeostasis, the structure of tissues, including muscles, blood vessels, and connective tissues, adapts to mechanical strains in the extracellular matrix (ECM). These strains originate from the differential growth of tissues or forces due to muscle contraction or gravity. Here we show using a computational model that by amplifying local strain cues, active cell contractility can facilitate and accelerate the reorientation of single cells to static strains. At the collective cell level, the model simulations show that active cell contractility can facilitate the formation of strings along the orientation of stretch. The computational model is based on a hybrid cellular Potts and finite-element simulation framework describing a mechanical cell-substrate feedback, where: 1) cells apply forces on the ECM, such that 2) local strains are generated in the ECM and 3) cells preferentially extend protrusions along the strain orientation. In accordance with experimental observations, simulated cells align and form stringlike structures parallel to static uniaxial stretch. Our model simulations predict that the magnitude of the uniaxial stretch and the strength of the contractile forces regulate a gradual transition between stringlike patterns and vascular networklike patterns. Our simulations also suggest that at high population densities, less cell cohesion promotes string formation.

Mathematical modeling for meshwork formation of endothelial cells in fibrin gels

Vasculogenesis is the earliest process in development for spontaneous formation of a primitive capillary network from endothelial progenitor cells. When human umbilical vein endothelial cells (HUVECs) are cultured on Matrigel, they spontaneously form a network structure which is widely used as an in vitro model of vasculogenesis. Previous studies indicated that chemotaxis or gel deformation was involved in spontaneous pattern formation. In our study, we analyzed the mechanism of vascular pattern formation using a different system, meshwork formation by HUVECs embedded in fibrin gels. Unlike the others, this experimental system resulted in a perfusable endothelial network in vitro. We quantitatively observed the dynamics of endothelial cell protrusion and developed a mathematical model for one-dimensional dynamics. We then extended the one-dimensional model to two-dimensions. The model showed that random searching by endothelial cells was sufficient to generate the observed network structure in fibrin gels.

Talking over the extracellular matrix: How do cells communicate mechanically?

Communication between cells enables them to coordinate their activity and is crucial for the differentiation, development, and function of tissues and multicellular organisms. Cell-cell communication is discussed almost exclusively as having a chemical or electrical origin. Only recently, a new mode of cell communication was elucidated: mechanical communication through the extracellular matrix (ECM). Cells can communicate mechanically by responding either to mechanical deformations generated by their neighbors or to a change in the mechanical properties of the ECM induced by a neighboring cell. This newly resolved mode of communication possesses unique features that complement the cellular ability to receive and share information, and to consequently act in a cooperative way with surrounding cells. Herein, we review several examples of mechanical communication, discuss their unique properties, and comment on the major challenges facing the field.

Engineering biomimetic periosteum with β-TCP scaffolds to promote bone formation in calvarial defects of rats

Background

There is a critical need for the management of large bone defects. The purpose of this study was to engineer a biomimetic periosteum and to combine this with a macroporous β-tricalcium phosphate (β-TCP) scaffold for bone tissue regeneration.

Methods

Rat bone marrow-derived mesenchymal stem cells (rBMSCs) were harvested and cultured in different culture media to form undifferentiated rBMSC sheets (undifferentiated medium (UM)) and osteogenic cell sheets (osteogenic medium (OM)). Simultaneously, rBMSCs were differentiated to induced endothelial-like cells (iECs), and the iECs were further cultured on a UM to form a vascularized cell sheet. At the same time, flow cytometry was used to detect the conversion rates of rBMSCs to iECs. The pre-vascularized cell sheet (iECs/UM) and the osteogenic cell sheet (OM) were stacked together to form a biomimetic periosteum with two distinct layers, which mimicked the fibrous layer and cambium layer of native periosteum. The biomimetic periostea were wrapped onto porous β-TCP scaffolds (BP/β-TCP) and implanted in the calvarial bone defects of rats. As controls, autologous periostea with β-TCP (AP/β-TCP) and β-TCP alone were implanted in the calvarial defects of rats, with a no implantation group as another control. At 2, 4, and 8 weeks post-surgery, implants were retrieved and X-ray, microcomputed tomography (micro-CT), histology, and immunohistochemistry staining analyses were performed.

Results

Flow cytometry results showed that rBMSCs were partially differentiated into iECs with a 35.1% conversion rate in terms of CD31. There were still 20.97% rBMSCs expressing CD90. Scanning electron microscopy (SEM) results indicated that cells from the wrapped cell sheet on the β-TCP scaffold apparently migrated into the pores of the β-TCP scaffold. The histology and immunohistochemistry staining results from in vivo implantation indicated that the BP/β-TCP and AP/β-TCP groups promoted the formation of blood vessels and new bone tissues in the bone defects more than the other two control groups. In addition, micro-CT showed that more new bone tissue formed in the BP/β-TCP and AP/β-TCP groups than the other groups.

Conclusions

Inducing rBMSCs to iECs could be a good strategy to obtain an endothelial cell source for prevascularization. Our findings indicate that the biomimetic periosteum with porous β-TCP scaffold has a similar ability to promote osteogenesis and angiogenesis in vivo compared to the autologous periosteum. This function could result from the double layers of biomimetic periosteum. The prevascularized cell sheet served a mimetic fibrous layer and the osteogenic cell sheet served a cambium layer of native periosteum. The biomimetic periosteum with a porous ceramic scaffold provides a new promising method for bone healing.

Stiffness-Controlled Thermoresponsive Hydrogels for Cell Harvesting with Sustained Mechanical Memory

Most mechanobiological investigations focused on in situ mechanical regulation of cells on stiffness-controlled substrates with few downstream applications, as it is still challenging to harvest and expand mechanically primed cells by enzymatic digestion (e.g., trypsin) without interrupting cellular mechanical memory between passages. This study develops thermoresponsive hydrogels with controllable stiffness to generate mechanically primed cells with intact mechanical memory for augmented wound healing. No significant cellular property alteration of the fibroblasts primed on thermoresponsive hydrogels with varied stiffness has been observed through thermoresponsive harvesting. When reseeding the harvested cells for further evaluation, softer hydrogels are proven to better sustain the mechanical priming effects compared to rigid tissue culture plate, which indicates that both the stiffness-controlled substrate and thermoresponsive harvesting are required to sustain cellular mechanical memory between passages. Moreover, epigenetics analysis reveals that thermoresponsive harvesting could reduce the rearrangement and loss of chromatin proteins compared to that of trypsinization. In vivo wound healing using mechanically primed fibroblasts shows featured epithelium and sebaceous glands, which indicates augmented skin recovery compared with trypsinized fibroblasts. Thus, the thermoresponsive hydrogel-based cell harvesting system offers a powerful tool to investigate mechanobiology between cell passages and produces abundant cells with tailored mechanical priming properties for cell-based applications.

Physical models of collective cell motility: from cell to tissue

In this article, we review physics-based models of collective cell motility. We discuss a range of techniques at different scales, ranging from models that represent cells as simple self-propelled particles to phase field models that can represent a cell's shape and dynamics in great detail. We also extensively review the ways in which cells within a tissue choose their direction, the statistics of cell motion, and some simple examples of how cell-cell signaling can interact with collective cell motility. This review also covers in more detail selected recent works on collective cell motion of small numbers of cells on micropatterns, in wound healing, and the chemotaxis of clusters of cells.

Measurement of dynamic cell-induced 3D displacement fields in vitro for traction force optical coherence microscopy

Traction force microscopy (TFM) is a method used to study the forces exerted by cells as they sense and interact with their environment. Cell forces play a role in processes that take place over a wide range of spatiotemporal scales, and so it is desirable that TFM makes use of imaging modalities that can effectively capture the dynamics associated with these processes. To date, confocal microscopy has been the imaging modality of choice to perform TFM in 3D settings, although multiple factors limit its spatiotemporal coverage. We propose traction force optical coherence microscopy (TF-OCM) as a novel technique that may offer enhanced spatial coverage and temporal sampling compared to current methods used for volumetric TFM studies. Reconstructed volumetric OCM data sets were used to compute time-lapse extracellular matrix deformations resulting from cell forces in 3D culture. These matrix deformations revealed clear differences that can be attributed to the dynamic forces exerted by normal versus contractility-inhibited NIH-3T3 fibroblasts embedded within 3D Matrigel matrices. Our results are the first step toward the realization of 3D TF-OCM, and they highlight the potential use of OCM as a platform for advancing cell mechanics research.

3D hybrid modelling of vascular network formation

We develop an off-lattice, agent-based model to describe vasculogenesis, the de novo formation of blood vessels from endothelial progenitor cells during development. The endothelial cells that comprise our vessel network are viewed as linearly elastic spheres that move in response to the forces they experience. We distinguish two types of endothelial cells: vessel elements are contained within the network and tip cells are located at the ends of vessels. Tip cells move in response to mechanical forces caused by interactions with neighbouring vessel elements and the local tissue environment, chemotactic forces and a persistence force which accounts for their tendency to continue moving in the same direction. Vessel elements are subject to similar mechanical forces but are insensitive to chemotaxis. An angular persistence force representing interactions with the local tissue is introduced to stabilise buckling instabilities caused by cell proliferation. Only vessel elements proliferate, at rates which depend on their degree of stretch: elongated elements have increased rates of proliferation, and compressed elements have reduced rates. Following division, the fate of the new cell depends on the local mechanical environment: the probability of forming a new sprout is increased if the parent vessel is highly compressed and the probability of being incorporated into the parent vessel increased if the parent is stretched. Simulation results reveal that our hybrid model can reproduce the key qualitative features of vasculogenesis. Extensive parameter sensitivity analyses show that significant changes in network size and morphology are induced by varying the chemotactic sensitivity of tip cells, and the sensitivities of the proliferation rate and the sprouting probability to mechanical stretch. Varying the chemotactic sensitivity directly influences the directionality of the networks. The degree of branching, and thereby the density of the networks, is influenced by the sprouting probability. Glyphs that simultaneously depict several network properties are introduced to show how these and other network quantities change over time and also as model parameters vary. We also show how equivalent glyphs constructed from in vivo data could be used to discriminate between normal and tumour vasculature and, in the longer term, for model validation. We conclude that our biomechanical hybrid model can generate vascular networks that are qualitatively similar to those generated from in vitro and in vivo experiments.

Computational Screening of Tip and Stalk Cell Behavior Proposes a Role for Apelin Signaling in Sprout Progression

Angiogenesis involves the formation of new blood vessels by sprouting or splitting of existing blood vessels. During sprouting, a highly motile type of endothelial cell, called the tip cell, migrates from the blood vessels followed by stalk cells, an endothelial cell type that forms the body of the sprout. To get more insight into how tip cells contribute to angiogenesis, we extended an existing computational model of vascular network formation based on the cellular Potts model with tip and stalk differentiation, without making a priori assumptions about the differences between tip cells and stalk cells. To predict potential differences, we looked for parameter values that make tip cells (a) move to the sprout tip, and (b) change the morphology of the angiogenic networks. The screening predicted that if tip cells respond less effectively to an endothelial chemoattractant than stalk cells, they move to the tips of the sprouts, which impacts the morphology of the networks. A comparison of this model prediction with genes expressed differentially in tip and stalk cells revealed that the endothelial chemoattractant Apelin and its receptor APJ may match the model prediction. To test the model prediction we inhibited Apelin signaling in our model and in an in vitro model of angiogenic sprouting, and found that in both cases inhibition of Apelin or of its receptor APJ reduces sprouting. Based on the prediction of the computational model, we propose that the differential expression of Apelin and APJ yields a "self-generated" gradient mechanisms that accelerates the extension of the sprout.

Review on experiment-based two- and three-dimensional models for wound healing

Traumatic and chronic wounds are a considerable medical challenge that affects many populations and their treatment is a monetary and time-consuming burden in an ageing society to the medical systems. Because wounds are very common and their treatment is so costly, approaches to reveal the responses of a specific wound type to different medical procedures and treatments could accelerate healing and reduce patient suffering. The effects of treatments can be forecast using mathematical modelling that has the predictive power to quantify the effects of induced changes to the wound-healing process. Wound healing involves a diverse and complex combination of biophysical and biomechanical processes. We review a wide variety of contemporary approaches of mathematical modelling of gap closure and wound-healing-related processes, such as angiogenesis. We provide examples of the understanding and insights that may be garnered using those models, and how those relate to experimental evidence. Mathematical modelling-based simulations can provide an important visualization tool that can be used for illustrational purposes for physicians, patients and researchers.

Modeling cell shape and dynamics on micropatterns

Adhesive micropatterns have become a standard tool to study cells under defined conditions. Applications range from controlling the differentiation and fate of single cells to guiding the collective migration of cell sheets. In long-term experiments, single cell normalization is challenged by cell division. For all of these setups, mathematical models predicting cell shape and dynamics can guide pattern design. Here we review recent advances in predicting and explaining cell shape, traction forces and dynamics on micropatterns. Starting with contour models as the simplest approach to explain concave cell shapes, we move on to network and continuum descriptions as examples for static models. To describe dynamic processes, cellular Potts, vertex and phase field models can be used. Different types of model are appropriate to address different biological questions and together, they provide a versatile tool box to predict cell behavior on micropatterns.

Cell-cell junctional mechanotransduction in endothelial remodeling

The vasculature is one of the most dynamic tissues that encounter numerous mechanical cues derived from pulsatile blood flow, blood pressure, activity of smooth muscle cells in the vessel wall, and transmigration of immune cells. The inner layer of blood and lymphatic vessels is covered by the endothelium, a monolayer of cells which separates blood from tissue, an important function that it fulfills even under the dynamic circumstances of the vascular microenvironment. In addition, remodeling of the endothelial barrier during angiogenesis and trafficking of immune cells is achieved by specific modulation of cell-cell adhesion structures between the endothelial cells. In recent years, there have been many new discoveries in the field of cellular mechanotransduction which controls the formation and destabilization of the vascular barrier. Force-induced adaptation at endothelial cell-cell adhesion structures is a crucial node in these processes that challenge the vascular barrier. One of the key examples of a force-induced molecular event is the recruitment of vinculin to the VE-cadherin complex upon pulling forces at cell-cell junctions. Here, we highlight recent advances in the current understanding of mechanotransduction responses at, and derived from, endothelial cell-cell junctions. We further discuss their importance for vascular barrier function and remodeling in development, inflammation, and vascular disease.

Extracellular matrix motion and early morphogenesis

For over a century, embryologists who studied cellular motion in early amniotes generally assumed that morphogenetic movement reflected migration relative to a static extracellular matrix (ECM) scaffold. However, as we discuss in this Review, recent investigations reveal that the ECM is also moving during morphogenesis. Time-lapse studies show how convective tissue displacement patterns, as visualized by ECM markers, contribute to morphogenesis and organogenesis. Computational image analysis distinguishes between cell-autonomous (active) displacements and convection caused by large-scale (composite) tissue movements. Modern quantification of large-scale 'total' cellular motion and the accompanying ECM motion in the embryo demonstrates that a dynamic ECM is required for generation of the emergent motion patterns that drive amniote morphogenesis.

Looking Beyond the Genes: The Interplay Between Signaling Pathways and Mechanics in the Shaping and Diversification of Epithelial Tissues

The core of Evo-Devo lies in the intuition that the way tissues grow during embryonic development, the way they sustain their structure and function throughout lifetime, and the way they evolve are closely linked. Epithelial tissues are ubiquitous in metazoans, covering the gut and internal branched organs, as well as the skin and its derivatives (ie, teeth). Here, we discuss in vitro, in vivo, and in silico studies on epithelial tissues to illustrate the conserved, dynamical, and complex aspects of their development. We then explore the implications of the dynamical and nonlinear nature of development on the evolution of their size and shape at the phenotypic and genetic levels. In rare cases, when the interplay between signaling and mechanics is well understood at the cell level, it is becoming clear that the structure of development leads to covariation of characters, an integration which in turn provides some predictable structure to evolutionary changes. We suggest that such nonlinear systems are prone to genetic drift, cryptic genetic variation, and context-dependent mutational effects. We argue that experimental and theoretical studies at the cell level are critical to our understanding of the phenotypic and genetic evolution of epithelial tissues, including carcinomas.

Dynamics of Cell Ensembles on Adhesive Micropatterns: Bridging the Gap between Single Cell Spreading and Collective Cell Migration

The collective dynamics of multicellular systems arise from the interplay of a few fundamental elements: growth, division and apoptosis of single cells; their mechanical and adhesive interactions with neighboring cells and the extracellular matrix; and the tendency of polarized cells to move. Micropatterned substrates are increasingly used to dissect the relative roles of these fundamental processes and to control the resulting dynamics. Here we show that a unifying computational framework based on the cellular Potts model can describe the experimentally observed cell dynamics over all relevant length scales. For single cells, the model correctly predicts the statistical distribution of the orientation of the cell division axis as well as the final organisation of the two daughters on a large range of micropatterns, including those situations in which a stable configuration is not achieved and rotation ensues. Large ensembles migrating in heterogeneous environments form non-adhesive regions of inward-curved arcs like in epithelial bridge formation. Collective migration leads to swirl formation with variations in cell area as observed experimentally. In each case, we also use our model to predict cell dynamics on patterns that have not been studied before.

Mathematical modelling of angiogenesis using continuous cell-based models

In this work, we develop a mathematical formalism based on a 3D in vitro model that is used to simulate the early stages of angiogenesis. The model treats cells as individual entities that are migrating as a result of chemotaxis and durotaxis. The phenotypes used here are endothelial cells that can be distinguished into stalk and tip (leading) cells. The model takes into account the dynamic interaction and interchange between both phenotypes. Next to the cells, the model takes into account several proteins such as vascular endothelial growth factor, delta-like ligand 4, urokinase plasminogen activator and matrix metalloproteinase, which are computed through the solution of a system of reaction–diffusion equations. The method used in the present study is classified into the hybrid approaches. The present study, implemented in three spatial dimensions, demonstrates the feasibility of the approach that is qualitatively confirmed by experimental results.

An in silico study on the role of smooth muscle cell migration in neointimal formation after coronary stenting

Excessive migration and proliferation of smooth muscle cells (SMCs) has been observed as a major factor contributing to the development of in-stent restenosis after coronary stenting. Building upon the results from in vivo experiments, we formulated a hypothesis that the speed of the initial tissue re-growth response is determined by the early migration of SMCs from the injured intima. To test this hypothesis, a cellular Potts model of the stented artery is developed where stent struts were deployed at different depths into the tissue. An extreme scenario with a ruptured internal elastic lamina was also considered to study the role of severe injury in tissue re-growth. Based on the outcomes, we hypothesize that a deeper stent deployment results in on average larger fenestrae in the elastic lamina, allowing easier migration of SMCs into the lumen. The data also suggest that growth of the neointimal lesions owing to SMC proliferation is strongly dependent on the initial number of migrated cells, which form an initial condition for the later phase of the vascular repair. This mechanism could explain the in vivo observation that the initial rate of neointima formation and injury score are strongly correlated.

Cell-ECM Interactions in Tumor Invasion

The cancer cells obtain their invasion potential not only by genetic mutations, but also by changing their cellular biophysical and biomechanical features and adapting to the surrounding microenvironments. The extracellular matrix, as a crucial component of the tumor microenvironment, provides the mechanical support for the tissue, mediates the cell-microenvironment interactions, and plays a key role in cancer cell invasion. The biomechanics of the extracellular matrix, particularly collagen, have been extensively studied in the biomechanics community. Cell migration has also enjoyed much attention from both the experimental and modeling efforts. However, the detailed mechanistic understanding of tumor cell-ECM interactions, especially during cancer invasion, has been unclear. This chapter reviews the recent advances in the studies of ECM biomechanics, cell migration, and cell-ECM interactions in the context of cancer invasion.

Primary Phenomenon in the Network Formation of Endothelial Cells: Effect of Charge

Blood vessels are essential organs that are involved in the supply of nutrients and oxygen and play an important role in regulating the body's internal environment, including pH, body temperature, and water homeostasis. Many studies have examined the formation of networks of endothelial cells. The results of these studies have revealed that vascular endothelial growth factor (VEGF) affects the interactions of these cells and modulates the network structure. Though almost all previous simulation studies have assumed that the chemoattractant VEGF is present before network formation, vascular endothelial cells secrete VEGF only after the cells bind to the substrate. This suggests VEGF is not essential for vasculogenesis especially at the early stage. Using a simple experiment, we find chain-like structures which last quite longer than it is expected, unless the energetically stable cluster should be compact. Using a purely physical model and simulation, we find that the hydrodynamic interaction retard the compaction of clusters and that the chains are stabilized through the effects of charge. The charge at the surface of the cells affect the interparticle potential, and the resulting repulsive forces prevent the chains from folding. The ions surrounding the cells may also be involved in this process.

Tip cell overtaking occurs as a side effect of sprouting in computational models of angiogenesis

BACKGROUND

During angiogenesis, the formation of new blood vessels from existing ones, endothelial cells differentiate into tip and stalk cells, after which one tip cell leads the sprout. More recently, this picture has changed. It has become clear that endothelial cells compete for the tip position during angiogenesis: a phenomenon named tip cell overtaking. The biological function of tip cell overtaking is not yet known. From experimental observations, it is unclear to what extent tip cell overtaking is a side effect of sprouting or to what extent it is regulated through a VEGF-Dll4-Notch signaling network and thus might have a biological function. To address this question, we studied tip cell overtaking in computational models of angiogenic sprouting in absence and in presence of VEGF-Dll4-Notch signaling.

RESULTS

We looked for tip cell overtaking in two existing Cellular Potts models of angiogenesis. In these simulation models angiogenic sprouting-like behavior emerges from a small set of plausible cell behaviors. In the first model, cells aggregate through contact-inhibited chemotaxis. In the second model the endothelial cells assume an elongated shape and aggregate through (non-inhibited) chemotaxis. In both these sprouting models the endothelial cells spontaneously migrate forwards and backwards within sprouts, suggesting that tip cell overtaking might occur as a side effect of sprouting. In accordance with other experimental observations, in our simulations the cells' tendency to occupy the tip position can be regulated when two cell lines with different levels of Vegfr2 expression are contributing to sprouting (mosaic sprouting assay), where cell behavior is regulated by a simple VEGF-Dll4-Notch signaling network.

CONCLUSIONS

Our modeling results suggest that tip cell overtaking can occur spontaneously due to the stochastic motion of cells during sprouting. Thus, tip cell overtaking and sprouting dynamics may be interdependent and should be studied and interpreted in combination. VEGF-Dll4-Notch can regulate the ability of cells to occupy the tip cell position in our simulations. We propose that the function of VEGF-Dll4-Notch signaling might not be to regulate which cell ends up at the tip, but to assure that the cell that randomly ends up at the tip position acquires the tip cell phenotype.

A global sensitivity analysis approach for morphogenesis models

BACKGROUND

Morphogenesis is a developmental process in which cells organize into shapes and patterns. Complex, non-linear and multi-factorial models with images as output are commonly used to study morphogenesis. It is difficult to understand the relation between the uncertainty in the input and the output of such 'black-box' models, giving rise to the need for sensitivity analysis tools. In this paper, we introduce a workflow for a global sensitivity analysis approach to study the impact of single parameters and the interactions between them on the output of morphogenesis models.

RESULTS

To demonstrate the workflow, we used a published, well-studied model of vascular morphogenesis. The parameters of this cellular Potts model (CPM) represent cell properties and behaviors that drive the mechanisms of angiogenic sprouting. The global sensitivity analysis correctly identified the dominant parameters in the model, consistent with previous studies. Additionally, the analysis provided information on the relative impact of single parameters and of interactions between them. This is very relevant because interactions of parameters impede the experimental verification of the predicted effect of single parameters. The parameter interactions, although of low impact, provided also new insights in the mechanisms of in silico sprouting. Finally, the analysis indicated that the model could be reduced by one parameter.

CONCLUSIONS

We propose global sensitivity analysis as an alternative approach to study the mechanisms of morphogenesis. Comparison of the ranking of the impact of the model parameters to knowledge derived from experimental data and from manipulation experiments can help to falsify models and to find the operand mechanisms in morphogenesis. The workflow is applicable to all 'black-box' models, including high-throughput in vitro models in which output measures are affected by a set of experimental perturbations.

The Force at the Tip--Modelling Tension and Proliferation in Sprouting Angiogenesis

Sprouting angiogenesis, where new blood vessels grow from pre-existing ones, is a complex process where biochemical and mechanical signals regulate endothelial cell proliferation and movement. Therefore, a mathematical description of sprouting angiogenesis has to take into consideration biological signals as well as relevant physical processes, in particular the mechanical interplay between adjacent endothelial cells and the extracellular microenvironment. In this work, we introduce the first phase-field continuous model of sprouting angiogenesis capable of predicting sprout morphology as a function of the elastic properties of the tissues and the traction forces exerted by the cells. The model is very compact, only consisting of three coupled partial differential equations, and has the clear advantage of a reduced number of parameters. This model allows us to describe sprout growth as a function of the cell-cell adhesion forces and the traction force exerted by the sprout tip cell. In the absence of proliferation, we observe that the sprout either achieves a maximum length or, when the traction and adhesion are very large, it breaks. Endothelial cell proliferation alters significantly sprout morphology, and we explore how different types of endothelial cell proliferation regulation are able to determine the shape of the growing sprout. The largest region in parameter space with well formed long and straight sprouts is obtained always when the proliferation is triggered by endothelial cell strain and its rate grows with angiogenic factor concentration. We conclude that in this scenario the tip cell has the role of creating a tension in the cells that follow its lead. On those first stalk cells, this tension produces strain and/or empty spaces, inevitably triggering cell proliferation. The new cells occupy the space behind the tip, the tension decreases, and the process restarts. Our results highlight the ability of mathematical models to suggest relevant hypotheses with respect to the role of forces in sprouting, hence underlining the necessary collaboration between modelling and molecular biology techniques to improve the current state-of-the-art.

Sprouting angiogenesis—a process by which new blood vessels grow from existing ones—is an ubiquitous phenomenon in health and disease of higher organisms, playing a crucial role in organogenesis, wound healing, inflammation, as well as on the onset and progression of over 50 different diseases such as cancer, rheumatoid arthritis and diabetes. Mathematical models have the ability to suggest relevant hypotheses with respect to the mechanisms of cell movement and rearrangement within growing vessel sprouts. The inclusion of both biochemical and mechanical processes in a mathematical model of sprouting angiogenesis permits to describe sprout extension as a function of the forces exerted by the cells in the tissue. It also allows to question the regulation of biochemical processes by mechanical forces and vice-versa. In this work we present a compact model of sprouting angiogenesis that includes the mechanical characteristics of the vessel and the tissue. We use this model to suggest the mechanism for the regulation of proliferation within sprout formation. We conclude that the tip cell has the role of creating a tension in the cells that follow its lead. On those first cells of the stalk, this tension produces strain and/or empty spaces, inevitably triggering cell proliferation. The new cells occupy the space behind the tip, the tension decreases, and the process restarts. The modelling strategy used, deemed phase-field, permits to describe the evolution of the shape of different domains in complex systems. It is focused on the movement of the interfaces between the domains, and not on an exhaustive description of the transport properties within each domain. For this reason, it requires a reduced number of parameters, and has been used extensively in modelling other biological phenomena such as tumor growth. The coupling of mechanical and biochemical processes in a compact mathematical model of angiogenesis will enable the study of lumen formation and aneurisms in the near future. Also, this framework will allow the study of the action of flow in vessel remodelling, since local forces can readily be coupled with cell movement to obtain the final vessel morphology.

Multiscale mechanobiology: computational models for integrating molecules to multicellular systems

Mechanical signals exist throughout the biological landscape. Across all scales, these signals, in the form of force, stiffness, and deformations, are generated and processed, resulting in an active mechanobiological circuit that controls many fundamental aspects of life, from protein unfolding and cytoskeletal remodeling to collective cell motions. The multiple scales and complex feedback involved present a challenge for fully understanding the nature of this circuit, particularly in development and disease in which it has been implicated. Computational models that accurately predict and are based on experimental data enable a means to integrate basic principles and explore fine details of mechanosensing and mechanotransduction in and across all levels of biological systems. Here we review recent advances in these models along with supporting and emerging experimental findings.

Emergent systems energy laws for predicting myosin ensemble processivity

In complex systems with stochastic components, systems laws often emerge that describe higher level behavior regardless of lower level component configurations. In this paper, emergent laws for describing mechanochemical systems are investigated for processive myosin-actin motility systems. On the basis of prior experimental evidence that longer processive lifetimes are enabled by larger myosin ensembles, it is hypothesized that emergent scaling laws could coincide with myosin-actin contact probability or system energy consumption. Because processivity is difficult to predict analytically and measure experimentally, agent-based computational techniques are developed to simulate processive myosin ensembles and produce novel processive lifetime measurements. It is demonstrated that only systems energy relationships hold regardless of isoform configurations or ensemble size, and a unified expression for predicting processive lifetime is revealed. The finding of such laws provides insight for how patterns emerge in stochastic mechanochemical systems, while also informing understanding and engineering of complex biological systems.

Three-dimensional numerical model of cell morphology during migration in multi-signaling substrates

Cell Migration associated with cell shape changes are of central importance in many biological processes ranging from morphogenesis to metastatic cancer cells. Cell movement is a result of cyclic changes of cell morphology due to effective forces on cell body, leading to periodic fluctuations of the cell length and cell membrane area. It is well-known that the cell can be guided by different effective stimuli such as mechanotaxis, thermotaxis, chemotaxis and/or electrotaxis. Regulation of intracellular mechanics and cell's physical interaction with its substrate rely on control of cell shape during cell migration. In this notion, it is essential to understand how each natural or external stimulus may affect the cell behavior. Therefore, a three-dimensional (3D) computational model is here developed to analyze a free mode of cell shape changes during migration in a multi-signaling micro-environment. This model is based on previous models that are presented by the same authors to study cell migration with a constant spherical cell shape in a multi-signaling substrates and mechanotaxis effect on cell morphology. Using the finite element discrete methodology, the cell is represented by a group of finite elements. The cell motion is modeled by equilibrium of effective forces on cell body such as traction, protrusion, electrostatic and drag forces, where the cell traction force is a function of the cell internal deformations. To study cell behavior in the presence of different stimuli, the model has been employed in different numerical cases. Our findings, which are qualitatively consistent with well-known related experimental observations, indicate that adding a new stimulus to the cell substrate pushes the cell to migrate more directionally in more elongated form towards the more effective stimuli. For instance, the presence of thermotaxis, chemotaxis and electrotaxis can further move the cell centroid towards the corresponding stimulus, respectively, diminishing the mechanotaxis effect. Besides, the stronger stimulus imposes a greater cell elongation and more cell membrane area. The present model not only provides new insights into cell morphology in a multi-signaling micro-environment but also enables us to investigate in more precise way the cell migration in the presence of different stimuli.

Semi-stochastic cell-level computational modelling of cellular forces: application to contractures in burns and cyclic loading

A phenomenological model is formulated to model cellular forces on extracellular material. The model is capable of modelling both expansion and contractile forces. This work is based on the assumption of linear elasticity, which allows a superposition argument to arrive at fundamental expressions for cellular forces. It is also shown how the cellular forces can be implemented using different strategies, as well as an extension to cellular point sources. Illustrations are given for modelling a (permanent) contraction (e.g. a contracture) of burns and for cyclic loading by the cells.

A coupled model of neovessel growth and matrix mechanics describes and predicts angiogenesis in vitro

During angiogenesis, sprouting microvessels interact with the extracellular matrix (ECM) by degrading and reorganizing the matrix, applying traction forces, and producing deformation. Morphometric features of the resulting microvascular network are affected by the interaction between the matrix and angiogenic microvessels. The objective of this study was to develop a continuous-discrete modeling approach to simulate mechanical interactions between growing neovessels and the deformation of the matrix in vitro. This was accomplished by coupling an existing angiogenesis growth model which uses properties of the ECM to regulate angiogenic growth with the nonlinear finite element software FEBio (www.febio.org). FEBio solves for the deformation and remodeling of the matrix caused by active stress generated by neovessel sprouts, and this deformation was used to update the ECM into the current configuration. After mesh resolution and parameter sensitivity studies, the model was used to accurately predict vascular alignment for various matrix boundary conditions. Alignment primarily arises passively as microvessels convect with the deformation of the matrix, but active alignment along collagen fibrils plays a role as well. Predictions of alignment were most sensitive to the range over which active stresses were applied and the viscoelastic time constant in the material model. The computational framework provides a flexible platform for interpreting in vitro investigations of vessel-matrix interactions, predicting new experiments, and simulating conditions that are outside current experimental capabilities.