<title>Introduction</title> <p>Blood-vessel development is essential for myriad biological phenomena in healthy and diseased individuals, including wound healing and tumor growth ##REF##11001068##[1]##,##REF##16355210##[2]##. Blood vessels form either de novo, via <italic>vasculogenesis</italic> or by sprouting or splitting of existing blood vessels via <italic>angiogenesis</italic>.</p> <p>In vasculogenesis, dispersed <italic>endothelial cells</italic> (<italic>ECs</italic>; the cells lining the inner walls of fully-formed blood vessels) organize into a primary vascular plexus of solid cords which then remodel into a vascular network. ECs elongate parallel to the cords, with final aspect ratios of tens to one. Because the early stages of vasculogenesis depend on a single cell type, vasculogenesis is relatively easy to reproduce in vitro. When cultured in vitro on <italic>Matrigel</italic>, a commercial product mimicking the extracellular matrix (<italic>ECM</italic>; the mixture of proteins, growth-factors and carbohydrates surrounding cells in vivo), even in the absence of other cell types or positional cues, ECs organize into cords which form large-scale, honey-comb like patterns, with cords of ECs surrounding regions devoid of ECs. This network slowly reorganizes, with the size of the polygonal, cell-free <italic>lacunae</italic>, gradually increasing. This observation suggests that ECs have <italic>autonomous</italic> patterning ability, rather than following morphogen pre-patterns.</p> <p>The sprouting or splitting of existing blood vessels during <italic>angiogenesis</italic> is more complex. In the first step of angiogenesis, a vessel dilates and releases plasma proteins that induce a series of changes in EC behavior. The ECs which will form the sprout next detach from each other and from the surrounding smooth-muscle cells, destabilizing the vessel. These detached ECs proliferate, migrate out of the vessel and organize into a sprout. EC proliferation continues in the sprout and is fastest just behind the leading <italic>tip cell</italic>, which is selected using a lateral-inhibition mechanism mediated by Dll4 and Notch1 ##REF##17259973##[3]##. Finally, the sprout forms a lumen, secretes a basal lamina and associates with pericytes that stabilize the sprout to form a mature new vessel ##REF##10742145##[4]##.</p> <p>Two fundamental questions concerning vasculogenesis and angiogenesis and their relation to each other are: 1) Does blood-vessel formation require external patterning cues (<italic>pre-patterns</italic> of morphogens) to define the precise position of the ECs, or can ECs organize into vascular patterns autonomously, with external cues merely initiating and fine-tuning vascular morphogenesis? 2) Do vasculogenesis and angiogenesis require the same or different cell behaviors, molecular signals and biomechanics?</p> <title>Experimental Background</title> <p>Despite the biomedical importance of angiogenesis and vasculogenesis, existing experiments are sufficiently ambiguous that even the fundamental mechanisms guiding patterning are uncertain. Experiments suggest a central role for chemotaxis in both de novo and sprouting blood-vessel growth ##REF##10742145##[4]##–##REF##16355211##[6]##. ECs respond to, and often produce, a wide range of chemoattractants and chemorepellants, including the many isoforms of vascular-endothelial growth factor A (<italic>VEGF-A</italic>) ##REF##16355211##[6]##, the chemokine SDF-1 ##REF##11929756##[7]##,##REF##16840724##[8]##, which ECs secrete ##REF##11929756##[7]##, fibroblast growth factor 2 (<italic>FGF-2</italic>), which induces ECs in developing vessels to secrete VEGF ##REF##9647657##[9]##, Slit-2, which can act either as a chemoattractant or a chemorepellant depending on the receptor to which it binds ##REF##15664389##[10]##, and the chemorepelling semaphorins ##REF##15664389##[10]##.</p> <p>Which of these molecules (if any) govern vascular patterning is still unclear. The Torino Group (e.g., ##REF##12682010##[11]##,##REF##12688968##[12]##) argued that a VEGF-A was the short-range autocrine chemoattractant that their chemotaxis-based blood-vessel-growth model required, since ECs express receptors for VEGF (<italic>VEGFR-2</italic>), chemotax towards sources of VEGF under favorable conditions, and secrete VEGFs. However, experiments suggest that cell-autonomous secretion of VEGF is essential only for vascular maintenance, not for angiogenesis <italic>per se</italic>: mice genetically-engineered to lack the VEGF gene only in their ECs have normal vascular density and patterning, but impaired vascular homeostasis and EC survival ##REF##17719546##[13]##. A plausible, alternative cell-autonomous chemoattractant to guide EC aggregation is the chemokine SDF-1/CXCL12, which ECs both secrete and respond to ##REF##16840724##[8]##.</p> <p>However, based on experiments that suggest that ECs can follow stresses in the ECM (see, e.g., ##REF##7573362##[14]## for review), Manoussaki and Murray ##REF##8953213##[15]##, and Namy et al. ##REF##14969709##[16]## proposed that mechanical interactions rather than, or in addition to, chemical interactions govern vasculogenesis. Further complicating this picture, Szabo and coworkers ##REF##17358734##[17]## showed that non-vascular, glia or muscle cells cultured on rigid, plastic culture dishes in continuously-shaken medium can form linear structures. Such culture conditions should reduce both the formation of chemoattractant gradients or migration along stress lines in the ECM. In the absence of ECM, they hypothesized that cells preferentially move towards elongated structures. Szabo and coworkers ##REF##17358734##[17]## proposed two mechanisms for such cell behavior: cells would align to surrounding cells, or they would mechanotactically follow stress fields in the cytoskeleton of neighboring cells. The molecular mechanisms of such cell behavior remains unclear as is the relevance of these results to ECs.</p> <p>Angiogenesis and vasculogenesis also require a number of local, contact-dependent (<italic>juxtracrine</italic>) signals: Tip-cell selection during angiogenic sprouting depends on Delta-notch signaling ##REF##17259973##[3]##, while Eph receptor-ephrin ligand binding amplifies ECs' response to SDF-1 ##REF##16840724##[8]##. All ECs express vascular-endothelial-cadherin (<italic>VE-cadherin</italic>), a homophilic, trans-membrane cell-adhesion molecule, which appears to play a crucial role in vascular patterning ##REF##10207135##[18]##,##REF##18062955##[19]##. Besides its role in cell-cell adhesion, VE-cadherin has a signaling function that determines how ECs respond to VEGF-A. When ECs bind to other ECs through their VE-cadherin, VEGF-A reduces their motility and proliferation. In the absence of VE-cadherin binding, VEGF-A activates pathways related to actin polymerization and the cell cycle, enhancing cell motility and proliferation in sub-confluent monolayers, and causes preferential extension of pseudopods in directions with higher VEGF-A concentrations ##REF##15071551##[20]##. We hypothesize that VE-cadherin-binding acts locally to prevent extension of pseudopods in the direction of cell-cell contacts for all critical chemoattractants, not only to VEGF-A. VE-cadherin −/− double-knock-out mice develop abnormal vascular networks in the yolk sac ##REF##10207135##[18]##, with ECs forming isolated vascular islands instead of wild-type polygonal vascular networks. These mice also have defective angiogenic sprouting, suggesting that both vasculogenesis and angiogenesis require VE-cadherin. VE cadherin −/− ECs still form strong adhesive junctions, so loss of VE-cadherin-mediated signaling rather than loss of intercellular adhesion seems to be responsible for the knock-out phenotype ##REF##10207135##[18]##.</p> <title>Computational Background</title> <p>A number of models and simulations replicate features of in vitro vascular patterning and can help partially reconstruct minimal sets of behaviors ECs require to self-organize into polygonal, vascular patterns ##REF##12682010##[11]##, ##REF##12688968##[12]##, ##REF##8953213##[15]##–##REF##17358734##[17]##, ##UREF##1##[21]##–##REF##19526066##[23]##. Because of the experiments we discussed above, and others which have demonstrated that sprouting angiogenesis and vasculogenesis both require chemotaxis (see, e.g., ##REF##11929756##[7]##,##REF##16840724##[8]##,##REF##12810700##[24]##), most models of vasculogenesis assume that intercellular signaling occurs via a diffusible chemoattractant. Using continuum models deriving from the fluid-dynamic Burgers' equation, Preziosi and coworkers (called the <italic>Torino Group</italic> in this paper) showed that simulated ECs secreting a chemoattractant that attracts surrounding ECs, could self-organize into polygonal patterns similar to the patterns in EC cultures and in vivo ##REF##12682010##[11]##,##REF##12688968##[12]##,##UREF##2##[25]##,##UREF##3##[26]##. However, their work assumed that endothelial cells accelerate in chemical gradients, which is not plausible in the highly viscous, non-inertial environment of the ECM. Microfluidic evidence indicates that mammalian cells (HL60) rapidly reach a flow-dependent, constant velocity ##REF##15915253##[27]## in chemoattractant gradients rather than continuously accelerating. We have previously suggested that ##REF##16325173##[22]## a linear force-velocity relation is the most appropriate model of ECs' experimental response, with the velocity of ECs proportional to the strength of the gradient of the chemoattractant. However, in simulations of this simple model, isotropic ECs form well-separated rounded clusters instead of networks. We have shown that adding one of a number of mechanisms (including cell adhesion ##UREF##1##[21]## and cell elongation ##REF##16325173##[22]##) to chemotactic aggregation suffices to produce quasi-polygonal networks. Section “<xref ref-type="sec" rid="s2">Results</xref>” discusses these mechanisms in more detail.</p> <p>In the mechanical models of Manoussaki and Murray ##REF##8953213##[15]##, and Namy et al. ##REF##14969709##[16]## ECs pull on the elastic ECM and aggregate by <italic>haptotactically</italic> migrating along the resulting ECM stress lines. Surprisingly, the mathematical form of the chemical and mechanical models is practically identical. Because these mechanical models assume that ECs exert radially-symmetric stresses on the ECM, modeling stress fields and EC haptotaxis or EC secretion and response to a chemoattractant, results in the same cell movement. Since simulations of the two mechanisms are identical, distinguishing between the effects of chemical and mechanical mechanisms will require additional experiments (such experiments are currently underway in the Glazier laboratory (Shirinifard, Alileche and Glazier, <italic>preprint, 2008</italic>)).</p> <p>A separate set of simulations addresses angiogenesis. Many models of sprouting blood-vessel growth introduce blood-vessel-level phenomenology by hand through high-level rules for branching ##UREF##4##[28]##–##UREF##5##[30]##. Attempts to derive blood-vessel sprouting and splitting from the underlying behavior of ECs include Levine and coworkers' ##REF##11315313##[31]## model of the onset of angiogenic sprouting as a reinforced random walk, where the ECs degrade the ECM, which locally enhances EC motility and produces paths of degraded ECM, and Bauer and Jiang's ##REF##17277180##[32]## cell-based model of blood-vessel sprouting along externally generated morphogen gradients, which assumed that branch splitting results from ECM inhomogeneities. Neither model can explain both EC assembly and blood-vessel sprouting.</p> <p>Could the behavior of the individual ECs also explain aspects of blood-vessel sprouting? Because the same genetic machinery regulates both angiogenesis and vasculogenesis ##REF##10742145##[4]##, a common set of mechanisms is plausible. Manoussaki ##UREF##6##[33]## extended her mechanical model of vasculogenesis to describe angiogenesis by adding long-range, chemotactic guidance cues. In her simulations, ECs migrated from an aggregate towards a chemoattractant source and cell-traction-driven migration contracted the sprout into a narrow, vessel-like cord.</p> <p>In this paper we present an alternative chemotaxis-based mechanism that can produce networks both from dispersed ECs and EC clusters without requiring long-range guidance cues. Instead, in our model long-range signals would only steer the self-organized vessels, a more biologically-realistic mechanism. Extending simulations that we have briefly introduced elsewhere ##REF##19526066##[23]##, we show that VE-cadherin-mediated contact inhibition of chemotactic pseudopod projections, in combination with secretion of a diffusing, rapidly decaying chemoattractant by ECs, suffices to reproduce aspects of both de novo and sprouting blood-vessel growth. In our simulations ECs: a) secrete a chemoattractant and b) preferentially extend pseudopods up gradients of the chemoattractant, unless, c) contact inhibition locally prevents chemotactic pseudopod extension. Thus, cell-cell binding suppresses the extension of chemotactic pseudopods, while unbound cell surfaces in contact with the ECM continue to extend pseudopods towards sources of chemoattractant ##REF##12810700##[24]##. We compare two biologically-plausible scenarios for chemotaxis, one in which ECs actively extend and retract pseudopods along chemoattractant gradients, and one in which the pseudopods' retractions are chemotactically neutral. The second scenario suggests a sprouting mechanism where a secreted autocrine factor acts both as a long-range chemoattractant and a local inhibitor of pseudopod sprouting.</p>
<title>Materials and Methods</title> <title>The Glazier–Graner–Hogeweg (GGH) Model</title> <p>The GGH represents biological cells as patches of identical lattice indices  on a square or triangular lattice, where each index uniquely identifies, or labels a single biological cell. Connections (<italic>links</italic>) between neighboring lattice sites of unlike index  represent bonds between apposing cell membranes, where the <italic>bond energy</italic> is , assuming that the types and numbers of adhesive cell-surface proteins determine <italic>J</italic>. A penalty increasing with the cell's deviation from a designated target volume <italic>A</italic>
<sub>target</sub>(<italic>σ</italic>) imposes a <italic>volume constraint</italic> on the simulated ECs. We define the pattern's <italic>effective energy</italic>:where  and  are neighboring lattice sites (up to fourth-order neighbors), <italic>a</italic> is the current area of cell <italic>σ</italic>, <italic>A</italic>
<sub>target</sub>(<italic>σ</italic>) is its target area, <italic>λ</italic> represents a cell's resistance to compression, and the Kronecker delta is <italic>δ</italic>(<italic>x</italic>,<italic>y</italic>) = {1,<italic>x</italic> = <italic>y</italic>; 0,<italic>x</italic>≠<italic>y</italic>. Each lattice site represents an area of 2 µm×2 µm. Since we assume that ECs do not divide or grow during patterning, we set <italic>A</italic>
<sub>target</sub>(<italic>σ</italic>) = 50 lattice sites, corresponding to a cell diameter of about 16 µm, and <italic>λ</italic> = 25 for all cells. The ECs reside in a very thin layer of extracellular fluid, which is a generalized cell without a volume constraint and with <italic>σ</italic> = 0. We assume that the ECs and fluid sit on top of a rigid ECM through which the chemoattractant diffuses, but we do not represent this ECM in the GGH lattice. We also assume that the presence of the fluid does not disturb the chemoattractant distribution in the ECM. Unless we specify otherwise, we use a bond energy <italic>J</italic>(<italic>c</italic>,<italic>c</italic>) = 40 between the ECs, and <italic>J</italic>(<italic>c</italic>,<italic>M</italic>) = 20 between the ECs and the ECM. For these settings the ECs do not adhere without chemotaxis. We define a special, high <italic>cell-border energy J</italic>(<italic>c</italic>,B) = 100 to prevent ECs from adhering to the lattice boundaries. We use fixed boundary conditions.</p> <p>To mimic cytoskeletally-driven pseudopod extensions and retractions, we randomly choose a source lattice site , and attempt to copy its index  into a randomly-chosen neighboring lattice site . For better isotropy we select the source site from the twenty, first- to fourth-nearest neighbors ##REF##9905332##[46]##. During a <italic>Monte Carlo Step</italic> (<italic>MCS</italic>) we carry out <italic>N copy attempts</italic>, with <italic>N</italic> the number of sites in the lattice. We set the experimental time per MCS to 30 s; for this setting the simulated ECs move with nearly their experimental velocity ##REF##16325173##[22]##. We calculate how much the effective energy would change if we performed the copy, and accept the attempt with probability , where <italic>T</italic> defines the <italic>intrinsic cell motility</italic>. All our simulations, except those in ##FIG##10##Figures 11##–
##FIG##12##13##, use <italic>T</italic> = 50.</p> <p>In experiments, cells respond to chemoattractant gradients by executing a more-or-less-strongly biased random walk up or down the gradient, where, over times short enough to allow us to neglect adaptation, the velocity of the drift depends on the gradient strength and the absolute concentration. We therefore define a set of extensions to the basic GGH model which reproduce these empirical behaviors due to preferential extension and retraction of pseudopods up chemoattractant gradients ##REF##12810700##[24]## by including a <italic>chemical effective-energy change</italic> at each copy attempt ##UREF##1##[21]##,##UREF##8##[36]##,where <italic>c</italic> is the concentration of the chemoattractant, which we assume is present everywhere in a layer of ECM <italic>under</italic> the ECs,  is the target site,  the source site, and <italic>s</italic> regulates the saturation of the chemotactic response. Unless we specify otherwise, we set <italic>s</italic> = 0, in which case chemotaxis depends linearly on the chemoattractant gradient only, independent of the chemoattractant concentration. Δ<italic>H</italic>
<sub>chemotaxis</sub> → 0 for large values of <italic>s</italic> and, for <italic>s</italic>≠0, for high chemical concentrations. The <italic>chemotaxis coefficient</italic> is <italic>μ</italic> = <italic>Χ</italic>(<italic>c</italic>,<italic>M</italic>) at cell-ECM interfaces and <italic>μ</italic> = <italic>Χ</italic>(<italic>c</italic>,<italic>c</italic>) at cell-cell interfaces respectively. Setting <italic>Χ</italic>(<italic>c</italic>,<italic>c</italic>) = 0 <italic>Χ</italic>(<italic>c</italic>,<italic>M</italic>) = 500 ensures that chemotactic extensions occur only at cell-ECM interfaces, reflecting VE-cadherin's suppression of pseudopods. Both extending and retracting pseudopods contribute to the chemical effective-energy change. To implement <italic>pseudopod-extension-only chemotaxis</italic> (see ##FIG##10##Figures 11##–
##FIG##12##13##), where only <italic>extending</italic> pseudopods at the cell-ECM interface respond to the chemoattractant, cells experience a chemical effective-energy change only if the source lattice site  belongs to an EC, i.e.,
</p> <p>For a more detailed discussion of chemotaxis in the GGH model see ##UREF##11##[47]##. We solve the partial-differential equation for chemoattractant diffusion and degradation (Eq. 1) numerically using a finite-difference scheme on a lattice matching the GGH lattice. We use 15 diffusion steps per MCS, with Δ<italic>t</italic> = 2 s. For these parameters, the chemoattractant diffuses more rapidly than the ECs, enabling us to ignore advection in the medium as the cells push the fluid.</p> <p>Source code and parameters for the simulations in this paper are available online in ##SUPPL##1##Protocol S1## from the supporting material, and from <ext-link ext-link-type="uri" xlink:href="http://sourceforge.net/projects/tst">http://sourceforge.net/projects/tst</ext-link>. Parameter files for the simulations in this paper are included in ##SUPPL##0##Dataset S1##.</p> <title>Allantois Culture and Immunolabeling</title> <p>We dissected allantoides from mouse embryos at embryonic stages 7.5–8.0. We washed the explants in fresh, cold ePBS and pipetted them into fibronectin-coated (5 mg/ml) Delta-T culture dishes (Bioptechs, Butler, PA) containing high-glucose, phenol-red-free Dulbecco's modified Eagles' medium (<italic>DMEM</italic>) supplemented with 10% fetal bovine serum (<italic>FBS</italic>), 1% penicillin-streptomycin, and 1% L-glutamine (GibcoBRL, Grand Island, NY). We maintained the allantoic explants using standard culture conditions (37°C and 5% CO<sub>2</sub>/95% air atmosphere) in a custom-designed culture chamber for 12–24 hours in the presence of an endothelial-specific marker, CD34 monoclonal antibody (BD PharMingen, San Diego, CA) directly conjugated to Cy3 (Amersham Biosciences). We fixed the allantoides in 3% paraformaldehyde for 20 minutes at room temperature, followed by an ePBS wash. For VE-cadherin antibody perturbations, we added anti-VE-cadherin monoclonal antibody (BD PharMingen, San Diego, CA) at 25 µg/ml to the culture medium.</p> <title>Image Acquisition</title> <p>We observed the cultures with a 10× objective (0.30 N.A.) on an inverted, automated, wide-field, epifluorescence/differential-interference-contrast (<italic>DIC</italic>) microscope (Leica DMIRE2, Leica Microsystems, Germany). We recorded images (608×512 pixel spatial and 12-bit intensity resolution) with a cooled Retiga 1300 camera (QImaging, Burnaby, British Columbia) in 2×2 binned acquisition mode, using 100–300 ms exposures. Image acquisition and microscope settings used software described in ##REF##12067365##[44]##.</p>
<title>Results</title> <p>We modeled endothelial cell behavior at a mesoscopic level using the Glazier-Graner-Hogeweg (GGH) model, also known as the Cellular Potts Model (CPM) ##REF##10046374##[34]##–##UREF##9##[37]##. The GGH is a lattice-based Monte-Carlo approach that describes biological cells as spatially extended patches of identical lattice indices. Intercellular junctions and cell junctions to the ECM determine adhesive (or binding) energies. The GGH algorithm, which we describe in more detail in the section <xref ref-type="sec" rid="s4">Materials and Methods</xref>, models pseudopod protrusions by iteratively displacing cell interfaces, with a preference for displacements which reduce the local <italic>effective energy</italic> of the configuration. Cells reorganize to favor stronger rather than weaker cell-cell and cell-ECM bonds and shorter rather than longer cell boundaries. In addition to interface displacements that reduce the effective energy, active cell motility also allows displacements that increase the effective energy. The likelihood of these active displacements increases with the <italic>cell-motility</italic> parameter <italic>T</italic>. Further constraints regulate cell volumes, surface areas, and chemotaxis. To model chemotaxis, we use the Savill and Hogeweg ##UREF##8##[36]## algorithm that favors extensions and retractions of pseudopods up concentration gradients of a chemoattractant (see Eq. 3 in the section <xref ref-type="sec" rid="s4">Materials and Methods</xref>). In the simplest implementation of chemotaxis in the GGH, cell velocity is proportional to the strength of the chemical gradient, in general agreement with experiments; see, e.g., ##REF##16325173##[22]## (we discuss the details of chemotaxis implementation below in the subsections Sensitivity analysis and A dissipative sprouting mechanism and in <xref ref-type="sec" rid="s4">Materials and Methods</xref>; see especially Eq. 3).</p> <p>The advantage of the GGH over alternative cell-based modeling approaches ##UREF##10##[38]## that represent cells as point particles or fixed-sized spheres or ellipsoids is that we can differentiate between bound and unbound regions of cell membrane. The GGH naturally represents the stochastic, exploratory behavior of migrating cells, modeling it as the biased extension and retraction of pseudopods, instead of a biologically-implausible single force acting on cells' centers of mass as in some cell-based simulations.</p> <p>We described chemoattractant diffusion and degradation macroscopically, using a continuum approximation. In analogy to the Torino Group's continuum model of de novo blood-vessel growth ##REF##12688968##[12]##,##UREF##2##[25]##, ECs secrete a diffusing chemoattractant at a rate <italic>α</italic>, which degrades in the ECM at a rate <italic>ε</italic> (e.g., due to proteolytic enzymes or by binding to ECM components), obeying:where  inside cells and is  in the ECM. Because we wish to compare our simulations to experimental yolk-sac cultures, where the vascular patterns are essentially monolayers, we use a two-dimensional GGH.</p> <p>We set the chemoattractant's secretion rate by cells <italic>α</italic> = 10<sup>−3</sup>
<italic>s</italic>
<sup>−1</sup>, its decay rate <italic>ε</italic> = <italic>α</italic>, and its diffusion constant in ECM to a slow <italic>D</italic> = 10<sup>−13</sup> m<sup>2</sup> s<sup>−1</sup>. These parameter values produce steeper gradients than those for VEGF-A<sub>165</sub>, the chemoattractant which Gamba et al. suggested was responsible for vasculogenesis, which has a diffusion coefficient of <italic>D</italic>∼10<sup>−11</sup> m<sup>2</sup> s<sup>−1</sup>
##REF##12688968##[12]##. The diffusion coefficient of SDF-1/CXCL12 is in the range of 1.7×10<sup>−13</sup> m<sup>2</sup> s<sup>−1</sup>
##REF##17066171##[39]##. However, the phenomena we observe in our simulations hold over a large range of diffusion coefficients.</p> <title>EC Aggregation and Vasculogenesis in the Absence of Contact Inhibition</title> <p>In in vitro cultures of mouse allantois explants, ECs (fluorescently labeled in red) organized into polygonal patterns (##FIG##0##Figure 1A–1C##). When we blocked VE-cadherin receptors with anti-VE-cadherin antibodies, thus preventing VE-cadherin receptors from binding to those on apposing cells, the mouse ECs formed isolated vascular islands (##FIG##0##Figure 1D–1F##). We hypothesize that anti-VE-cadherin's antibody blockage of VE-cadherin signaling prevents contact inhibition of chemotactic motility, sensitizing the endothelial cells to the chemoattractant at cell–cell interfaces.</p> <p>In our corresponding simulations (##FIG##1##Figure 2A–2C## and ##SUPPL##4##Video S1##), we randomly distributed 1,000 ECs, each with an area of ∼200 µm<sup>2</sup> over an area of ≈700 µm×700 µm (333×333 lattice sites, or <italic>pixels</italic>, of 2 µm×2 µm each), which we positioned inside a larger lattice of 1,00 µm×1,00 µm to minimize boundary effects. In this cell-based simulation of the Torino Group's continuum model ##REF##12682010##[11]##,##REF##12688968##[12]##, without endothelial-cell acceleration in chemoattractant gradients our cells form disconnected, vascular islands rather than a vascular network. We would expect this result, because with the more realistic chemotactic response we employ, the Torino Group's model reduces to the classic Keller-Segel equations ##REF##5462335##[40]## of chemotactic aggregation ##UREF##2##[25]##, which, like our simulations, form isolated vascular islands. Apparently, the basic Torino-Group model of chemotactic cell aggregation misses a biological mechanism essential for vasculogenesis. We have previously suggested a number of additional mechanisms, any one of which, together with cell aggregation, suffices to induce vasculogenesis-like patterning. For example, when we gave the ECs the elongated shapes observed in later stages of experiments, neighboring cells aligned with each other, causing cell clusters to elongate and interconnect, creating a vascular network, in a mechanism similar to Szabo's ##REF##17358734##[17]##. These vascular networks remodel gradually, with dynamics resembling those of in vitro vascular networks. The causes of cell elongation in experiments are not clear. ECs could elongate either cell-autonomously (e.g., by remodeling their cytoskeletons), or non-cell-autonomously, by maximizing their contact areas with surrounding cells or by aligning to morphogen gradients in the ECM ##REF##16325173##[22]##. Unless we state otherwise, in this paper we neglect cell-autonomous elongation.</p> <p>Even without strong cell-cell adhesion the ECs can form vascular-like structures in simulations of vasculogenesis if the <italic>diffusion length</italic> of the chemoattractant (the length <italic>L</italic> over which the concentration drops to half its value at the EC membrane) is short enough, because the ECs align with the chemical gradients ##REF##19526066##[23]##. This length scale <italic>L</italic> depends on the diffusion coefficient <italic>D</italic> and the chemoattractant decay rate <italic>ε</italic> as 
##REF##12688968##[12]##.</p> <title>Sprouting Angiogenesis in the Absence of Contact-Inhibition</title> <p>To investigate whether the Torino-Group Model could reproduce sprouting angiogenesis, we started our simulations with rounded clusters of simulated ECs representing a blood vessel's surface after degradation of the ECM, keeping the simulation parameters unchanged from ##FIG##1##Figure 2##. As in vasculogenesis, cell-elongation sufficed to drive angiogenesis-like sprouting (see ##FIG##2##Figure 3A–3C##), where we used a length constraint, see ##REF##16325173##[22]##). EC clusters also produced sprouts for strong cell-cell adhesion (i.e., for values of <italic>J</italic>(<italic>c</italic>,<italic>c</italic>)&lt;10); ##FIG##2##Figure 3D–3F##), via a mechanism similar to the cell-elongation-dependent mechanism for vasculogenesis ##REF##16325173##[22]##. Adhesion-independent sprouting occurred only for a narrow range of very small diffusion constants of the chemoattractant, between <italic>D</italic>&lt;2·10<sup>−14</sup> m<sup>2</sup> s<sup>−1</sup> and <italic>D</italic>&gt;2·10<sup>−14</sup> m<sup>2</sup> s<sup>−</sup> (see ##FIG##2##Figure 3G–3I##). The allowable range of <italic>D</italic> increased for bigger cells ##REF##19526066##[23]##. We also systematically screened for sprouting in the absence of contact-inhibited chemotaxis. We present the results of these screens in the section Sensitivity analysis and in ##SUPPL##2##Figure S1##, but we defer an in-depth study of these phenomena to our future work.</p> <title>Contact-Inhibited Chemotaxis in De Novo Blood Vessel Growth</title> <p>In this paper, we focus on the role of contact-inhibited chemotaxis in sprouting blood-vessel growth. We hypothesize that VE-cadherin's local inhibition of chemotaxis-induced pseudopod extensions at EC-EC boundaries, may be responsible for ECs' self-organization into vascular-like networks.</p> <p>We modeled contact inhibition of chemotaxis in our simulations by suppressing chemotaxis at cell-cell interfaces. Thus, only interfaces between cells and ECM respond to the chemoattractant. ##FIG##1##Figure 2D## and ##SUPPL##5##Video S2## and ##SUPPL##6##Video S3## show typical simulations of de novo blood-vessel growth with contact inhibition. The ECs assemble into a structure resembling a capillary plexus: cords of cells enclose lacunae, which grow slowly. Smaller lacunae shrink and disappear, while larger lacunae subdivide via vessel sprouting as, for example, in the quail yolk sac ##REF##9109485##[41]##.</p> <title>Contact-Inhibited Chemotaxis in Blood Vessel Sprouting</title> <p>To investigate the role of contact-inhibited chemotaxis in blood vessel sprouting, we ran a set of simulations with a large cluster of endothelial cells representing a blood vessel's surface after degradation of the ECM, keeping all simulation parameters the same as those in ##FIG##1##Figure 2D##. The surface of the cluster first roughens, with some cells protruding from the surface, then digitates into a structure reminiscent of a primary vascular plexus (##FIG##3##Figure 4A–4C## and ##SUPPL##7##Video S4## and ##SUPPL##8##Video S5##), the first type of structure to develop in both de novo and sprouting blood-vessel growth ##REF##9109485##[41]##. The sprouting instability requires contact inhibition of chemotaxis. Without it, the clusters remained rounded and compact (##FIG##3##Figure 4D##). Thus our simulations suggest that a process operating at the level of individual cells—chemotaxis with contact inhibition—may drive in vitro blood-vessel growth both sprouting and de novo.</p> <p>What drives blood vessel sprouting in our model? At equilibrium, the chemoattractant has a quasi-Gaussian profile across the cluster. It levels off towards the cluster's center, while its inflection point is at the cluster boundary. Chemotaxis produces a continuous, inward, normal force at the cluster boundary, creating a buckling instability (see, e.g., ##REF##10990656##[42]##); chemotactic forces also compress small initial bumps laterally, producing sprouts. Since contact inhibition of chemotaxis leaves the interior cells insensitive to the chemoattractant, ingressing surface cells easily push them aside. When we omit contact inhibition of motility to mimic anti-VE-cadherin-antibody-treated allantois cultures, the interior cells also feel the inward-directed chemotactic forces and resist displacement (##FIG##3##Figure 4D## and ##SUPPL##9##Video S6##).</p> <p>To explore this idea, we varied the ratio of the chemotactic response at cell–cell interfaces relative to the chemotactic response at cell-ECM interfaces (<italic>Χ</italic>(<italic>c</italic>,<italic>c</italic>)/<italic>Χ</italic>(<italic>c</italic>,<italic>M</italic>)), where <italic>Χ</italic>(<italic>c</italic>,<italic>c</italic>)is the ECs' sensitivity to the chemoattractant at cell-cell interfaces and <italic>Χ</italic>(<italic>c</italic>,<italic>M</italic>) the sensitivity at cell-ECM interfaces (see the section <xref ref-type="sec" rid="s4">Materials and Methods</xref> for details). We looked for sprouting in clusters of 128 cells, each of area ∼200 µm<sup>2</sup>, placed in a 400 µm×400 µm lattice, keeping all other parameters unchanged from their values in ##FIG##3##Figure 4##.</p> <p>We defined the clusters' compactness after 10,000 Monte Carlo Steps (the time unit of the simulation, see the section <xref ref-type="sec" rid="s4">Materials and Methods</xref>, with 1 MCS equivalent to about 30 s) to be <italic>C</italic> = <italic>A</italic>
<sub>cluster</sub>/<italic>A</italic>
<sub>hull</sub>, the ratio between the cluster's area, <italic>A</italic>
<sub>cluster</sub>, and the area of its convex hull (that is the tightest possible “gift wrapping” around the cluster), <italic>A</italic>
<sub>hull</sub>. The compactness <italic>C</italic> = 1 for a perfectly circular cluster, whereas <italic>C</italic> → 0 for highly branched or dispersed clusters of cells.</p> <p>We found a phase transition at (<italic>Χ</italic>(<italic>c</italic>,<italic>c</italic>)/<italic>Χ</italic>(<italic>c</italic>,<italic>M</italic>))≈0.5 separating sprouting from non-sprouting clusters (##FIG##4##Figure 5##), suggesting that the sprouting instability only occurs when the core of the cluster behaves as a fluid: because each cell's volume is nearly conserved (apart from small fluctuations around its target volume), the core cells can only release the pressure the ingressing cells exert on them by moving outwards as sprouts. Our ongoing work characterizes this instability mathematically, proving that the cluster self-organizes into a network structure with fixed cord width (A. Shirinifard and J. A. Glazier, preprint 2008).</p> <p>To validate our model against published EC tracking experiments ##REF##18062955##[19]##, we compared the trajectories of cells in sprouting and non-sprouting clusters. ##FIG##5##Figure 6A–6D## shows the trajectories of ten cells in a sprouting cluster (with contact-inhibition; ##FIG##5##Figure 6A–6B##), and ten cells in a non-sprouting cluster (without contact-inhibition; ##FIG##5##Figure 6C–6D##). In non-sprouting clusters, cells followed random-walk trajectories, while in sprouting clusters, they followed biased random-walk trajectories. To further characterize cell motility, we measured cells' average displacements and velocities over 10 independent simulations of 128 cells each. In sprouting clusters, the cells moved further during a given interval than in non-sprouting clusters. Thus, the cell velocity 
##REF##18062955##[19]## is larger during sprouting if the interval Δ<italic>t</italic> between subsequent cell positions is sufficiently large (here we use Δ<italic>t</italic> = 2.5 h as in Perryn et al. ##REF##18062955##[19]##); for shorter intervals (e.g., 30 s) the cell velocity is highest in non-sprouting clusters (not shown), indicating that ECs in sprouting clusters moved faster, but had a somewhat slower random motility.</p> <p>Our simulations agree with recent experiments tracking ECs in embryonic mouse allantoides ##REF##18062955##[19]## that measured the cell-autonomous motility of ECs cells in allantoides relative to the motility of the surrounding mesothelium in which the ECs reside. Administration of anti-VE-cadherin antibodies reduced both cell-autonomous motion and net displacement of ECs. Thus, our simulations suggest that VE-cadherin's role as a contact-dependent inhibitor of cell motility suffices to explain the reduced cell motility observed in anti-VE-cadherin-treated allantoides cultures.</p> <title>Sensitivity Analysis</title> <p>Contact-inhibited sprouting occurs for a wide range of parameter values. In most of our simulations we set the EC-EC adhesion equal to the EC-ECM adhesion (i.e., <italic>J</italic>(<italic>c</italic>,<italic>c</italic>) = <italic>2J</italic>(<italic>c</italic>,<italic>M</italic>); the factor of 2 arises because we model the ECM as a single large generalized cell), which is equivalent to setting the surface tension of the cluster to zero ##UREF##7##[35]##. Zero surface tension clarifies the role of contact inhibition in sprouting, but real ECs adhere strongly to each other via <italic>adherens junctions</italic>
##REF##10207135##[18]##. In ##FIG##6##Figure 7## and in ##SUPPL##10##Video S7##, ##SUPPL##11##S8##, ##SUPPL##12##S9##, ##SUPPL##13##S10##, ##SUPPL##14##S11##, ##SUPPL##15##S12##, ##SUPPL##16##S13##, ##SUPPL##17##S14##, ##SUPPL##18##S15##, ##SUPPL##19##S16##, ##SUPPL##20##S17##, ##SUPPL##21##S18##, ##SUPPL##22##S19##, ##SUPPL##23##S20##, ##SUPPL##24##S21##, ##SUPPL##25##S22##, we studied the effect of cell-cell adhesion on sprouting in clusters of 128 cells (256 cells in the Videos). For stronger EC-EC adhesion, equivalent to positive surface tension, <italic>J</italic>(<italic>c</italic>,<italic>c</italic>)&lt;<italic>2J</italic>(<italic>c</italic>,<italic>M</italic>), the sprouts are longer and thinner and the network less compact than for zero surface tension. For very weak EC-EC adhesion <italic>J</italic>(<italic>c</italic>,<italic>c</italic>)≫2<italic>J</italic>(<italic>c</italic>,<italic>M</italic>), equivalent to strong negative surface tension, the ECs separate from each other, so contact-inhibition no longer occurs, and the clusters do not sprout. For small negative surface tensions, with values of <italic>J</italic>(<italic>c</italic>,<italic>c</italic>)&gt;2<italic>J</italic>(<italic>c</italic>,<italic>M</italic>), chemotaxis overcomes the negative surface tension, so ECs still touch each other and sprouting occurs as for zero surface tension, producing thickened sprouts and elongated clusters. The insets to ##FIG##6##Figure 7## and ##SUPPL##23##Video S20##, ##SUPPL##24##S21##, ##SUPPL##25##S22## show the results for 50≤<italic>J</italic>(<italic>c</italic>,<italic>c</italic>)≤70.</p> <p>We also investigated how sprouting depends on the chemotactic strength <italic>Χ</italic>(<italic>c</italic>,<italic>M</italic>) (##FIG##7##Figure 8## and ##SUPPL##26##Video S23##, ##SUPPL##27##S24##, ##SUPPL##28##S25##, ##SUPPL##29##S26##, ##SUPPL##30##S27##
##SUPPL##31##S28##, ##SUPPL##32##S29##, ##SUPPL##33##S30##, ##SUPPL##34##S31##, ##SUPPL##35##S32##, ##SUPPL##36##S33##). For <italic>Χ</italic>(<italic>c</italic>,<italic>M</italic>) = 500, most vascular cords are two cells wide (##SUPPL##27##Video S24##), while for For <italic>Χ</italic>(<italic>c</italic>,<italic>M</italic>)&gt;500 the cords become thinner and longer, with cords only one cell wide (##SUPPL##28##Videos S25##, ##SUPPL##29##S26##, ##SUPPL##30##S27##, ##SUPPL##31##S28##, ##SUPPL##32##S29##, ##SUPPL##33##S30##, ##SUPPL##34##S31##, ##SUPPL##35##S32##, ##SUPPL##36##S33##). For higher chemotactic forces, the cells intercalate, moving to the chemical gradients' peak. We have derived the conditions for this <italic>folding instability</italic> in our ongoing work (A. Shirinifard and J. A. Glazier, preprint, 2008). Higher chemotactic strengths increase ruffling of the cluster boundary, reducing the cluster's compactness in the absence of contact inhibition (##FIG##7##Figure 8##).</p> <p>We assumed that ECs extend or retract pseudopods depending on the difference in chemoattractant concentration between the retracted and extended positions, independent of the absolute chemoattractant concentrations. However, at higher chemoattractant concentrations, most chemoattractant receptors will saturate with chemoattractant and become insensitive to chemoattractant levels. To study the effect of saturated chemotactic response ##UREF##1##[21]## on angiogenic sprouting, we varied the saturation parameter s (see Eq. 3 in <xref ref-type="sec" rid="s4">Materials and Methods</xref>) leaving all other parameters unchanged. For <italic>s</italic> = 0, the chemotactic response is linear; for <italic>s</italic>&gt;0, the response to the chemoattractant gradient vanishes at high concentrations (see <xref ref-type="sec" rid="s4">Materials and Methods</xref>). For small positive <italic>s</italic>, the clusters sprout normally (see ##FIG##8##Figure 9## and ##SUPPL##37##Videos S34##, ##SUPPL##38##S35##, ##SUPPL##39##S36##); however, for large <italic>s</italic>, the chemotactic response weakens at the chemoattractant levels present at the edge of the cell cluster; thus cells no longer chemotact towards the cluster's interior and the sprouting instability disappears (##SUPPL##40##Videos S37##, ##SUPPL##41##S38##, ##SUPPL##42##S39##). We could test this prediction experimentally by partially inactivating the ECs' chemoattractant receptors. We observed the same effect when we increased the chemoattractant secretion rate for moderate response saturation (<italic>s</italic> = 0.05; see ##SUPPL##3##Figure S2##, bottom panel) leading to higher overall chemoattractant concentrations. We could test this situation experimentally by overexpressing the chemoattractant in ECs. Since for unsaturated chemotactic response (<italic>s</italic> = 0), multiplying the chemoattractant concentrations is equivalent to multiplying the chemotactic strength (<italic>Χ</italic>(<italic>c</italic>,<italic>M</italic>)) by the same factor, increasing the secretion rate first thins and lengthens the cords by increasing the chemotactic strength, then eventually prevents sprouting as the chemotactic response saturates. This effect is most apparent for <italic>s</italic> = 0.01 (##SUPPL##3##Figure S2##, top panel).</p> <p>In the Torino Group's continuum model, the separation between the cords increases with the diffusion length <italic>L</italic> of the chemoattractant, ##FIG##9##Figure 10## and ##SUPPL##43##Videos S40##, ##SUPPL##44##S41##, ##SUPPL##45##S42##, ##SUPPL##46##S43##, ##SUPPL##47##S44##, ##SUPPL##48##S45##, ##SUPPL##49##S46## show sprouting clusters for a range of diffusion lengths. In agreement with the Torino Group's model, longer diffusion lengths produce thicker cords with larger intercord spaces. The clusters do not sprout well when <italic>L</italic> approaches the EC-cluster diameter. Clusters consisting of 1,024 cells sprout for <italic>D</italic>&gt;3·10<sup>−13</sup> m<sup>2</sup> s<sup>−1</sup> (<italic>L</italic>&gt;17.3 µm), while 128-cell clusters do not (##FIG##9##Figure 10## and ##SUPPL##50##Video S47##, ##SUPPL##51##S48##, ##SUPPL##52##S49##, ##SUPPL##53##S50##, ##SUPPL##54##S51##, ##SUPPL##55##S52##, ##SUPPL##56##S53##). If the diffusion length is shorter than the ECs' diameter, the clusters dissociate: the ECs perform random walks with long persistence lengths, moving up the chemoattractant gradients they leave behind themselves (##SUPPL##43##Video S40## and ##SUPPL##50##S47##).</p> <title>A Dissipative Sprouting Mechanism</title> <p>In our simulations, the trailing edges of the ECs retract actively in response to the chemoattractant and exert an inward-normal, compressive force on the EC cluster. To check if sprouting requires this compressive force, we also simulated a situation in which only <italic>extending</italic> pseudopods at cell-ECM interfaces respond to the chemoattractant, while retraction is chemotactically neutral. Both sprouting-angiogenesis and vasculogenesis occurred, but required higher intrinsic cell motilities (larger values of the parameter <italic>T</italic>). ##FIG##10##Figure 11## shows the motilities required under both assumptions. We looked for sprouting after 5000 MCS (∼40 h) in clusters of 128 cells, each of area ≈200 µm<sup>2</sup>, placed in a 400 µm×400 µm lattice, with all other parameters the same as in ##FIG##3##Figure 4##. For <italic>T</italic>&lt;100, our original chemotaxis assumptions produced sprouts, while no sprouting occurred if pseudopods responded to the chemoattractant only during extension. For 100&lt;<italic>T</italic>&gt;400, both mechanisms produced sprouts. For <italic>T</italic>&gt;400, the ECs broke up into small pieces, a well-characterized, non-biological artifact of the GGH ##UREF##7##[35]##. With extension-only chemotaxis, sprouting was slightly slower than for standard, extension-retraction Savill-Hogeweg ##UREF##8##[36]## chemotaxis, as a plot of the time evolution of the clusters' compactness shows (##FIG##11##Figure 12## and ##SUPPL##57##Video S54##, ##SUPPL##58##S55##, ##SUPPL##59##S56##). However, at long times (<italic>t</italic>&gt;2500 MCS) the compactness of clusters decreased at identical rates for both methods.</p> <p>These results suggest an additional mechanism for blood-vessel sprouting: at the cluster surface, <italic>all</italic> pseudopod extensions increase the effective energy slightly, so the chemoattractant <italic>inhibits</italic> pseudopod extension. A recent experimental study ##REF##17038622##[43]## found that autocrine secretion of the sprouting inhibitor TGF-β1 enhances branching in mammary epithelial tubes. Our model suggests a mechanism by which an autocrine, secreted chemical can act <italic>both</italic> as a chemoattractant and as an inhibitor. The rates of pseudopod extensions and retractions are critical to pattern evolution (##FIG##10##Figure 11##). Cells in growing tips see a shallower gradient than do those in valleys between the tips (see, e.g., ##FIG##3##Figure 4B##), so pseudopod extensions at growing tips are more frequent than in the valleys between tips because they have a lower effective-energy cost. During sprouting, conservation of cell area requires that the cells in the valleys must retract, while those in the tips protrude. In the Savill-Hogeweg algorithm, retraction is energetically favorable, while it is energetically neutral in our pseudopod-extension-only chemotaxis algorithm, making the net change in effective energy positive with a rate depending on the cell motility. The effective-energy change is negative in the Savill-Hogeweg algorithm and thus nearly independent of <italic>T</italic> (##FIG##12##Figure 13##, where <italic>H</italic>
<sub>0</sub> is the initial effective energy).</p>
<title>Discussion</title> <p>We have shown that a single set of cell behaviors, i.e., contact-inhibited chemotaxis to an autocrine, secreted chemoattractant can explain aspects of both de novo and sprouting blood-vessel growth. Our results suggest that branching in aggregates of chemotacting ECs could result from two separate effects of the same mechanism. For low cell motilities <italic>T</italic>, i.e., a low probability for active, dissipative cellular protrusion, the branching resembles a buckling instability (see, e.g., ##REF##10990656##[42]##), in which the surface cells exert a surface-normal force on the cluster's inner core. For larger cell motilities, the shallower chemoattractant gradients at protrusions make the ECs there more likely to extend outward-directed pseudopods than cells in the valleys between the protrusions.</p> <p>While we have adopted the Torino Group's assumption that ECs chemotax in response to gradients of a diffusible, autocrine, secreted chemoattractant ##REF##12688968##[12]##,##UREF##2##[25]##, our simulation also reproduces continuum models that assume that ECs stress the ECM ##REF##8953213##[15]##, which either pulls on the surrounding ECs, provides haptotactic cues for active EC migration ##REF##14969709##[16]##, or both ##UREF##3##[26]##. Because these models assume that ECs exert radially-symmetric stresses on the ECM, the underlying mathematical descriptions of the chemotactic and haptotactic mechanisms are equivalent. In both cases, contact inhibition should still operate and the patterning mechanism we have proposed should still apply, with traction or haptotaxis replacing chemotaxis and the mechanical screening length replacing the diffusion length. Our simulation may also apply to the formation of linear structures by non-vascular, glia or muscle cells cultured on rigid, plastic culture dishes in continuously-shaken medium ##REF##17358734##[17]## in which cells explore their environment using long filopodia, then move towards their neighbors by pulling themselves along bound filopodia. Thus, the combination of cell aggregation and contact-inhibition that drives patterning in our model, could also occur without chemical gradients and even without ECM.</p> <p>Our simulations also allow us to clarify a number of subtleties concerning the interpretation of our own and others' experiments in which blocking VE-cadherin interfered with normal vascular patterning. In our in vitro experiments, anti-VE-cadherin treatment caused ECs to round, in addition to its hypothesized effect on contact inhibition, so our experiments cannot rule out the possibility that the anti-VE-cadherin treatment inhibited vascular patterning because of its effect on EC shape. A further complication is that anti-VE-cadherin treatment could conceivably reduce the adhesion between ECs. As we noted above, In VE-cadherin −/− knock-out mice, ECs still form strong adhesive junctions ##REF##10207135##[18]##, suggesting that VE-cadherin is not required for EC-EC binding.</p> <p>Our simulations show that the contact-inhibition patterning mechanism operates over a wide range of cell-cell adhesions, suggesting that changes in adhesivity are not significant provided that contact-inhibition persists, and independent of cell shape ##REF##19526066##[23]##, suggesting that the shape change is not significant. However, we have also shown that strong cell-cell adhesion plus chemotaxis can produce vascular-like patterns in simulations ##UREF##1##[21]##. Fortunately, the three vascular patterning mechanisms (contact-inhibition, cell-elongation and cell-cell adhesion) have vastly different kinetics ##REF##16325173##[22]##. Thus time-lapse microscopy experiments ##REF##18062955##[19]##,##REF##12067365##[44]## quantifying the kinetics of capillary-plexus development (see, e.g., ##REF##16325173##[22]##), will allow us to definitively distinguissh among these three patterning mechanisms. Already, we can say that adhesion-driven patterning is so slow and requires such strong adhesion that it appears incompatible with the available qualitative data from experiments.</p> <p>To further test if VE-cadherin-mediated, contact-dependent signaling to VEGF-R2 ##REF##15071551##[20]##, rather than VE-cadherin's function as a cell-adhesion molecule is responsible for the effects of anti-VE-cadherin treatment in mouse yolk sacs, we could experimentally block signal transduction from VE-cadherin to VEGFR-2, specifically interfering with VE-cadherin's signaling function, while leaving its role as an adhesion molecule intact. A possible target would be CD148, which phosphorylates VEGFR-2 after VE-cadherin binding ##REF##15071551##[20]##,##REF##12588999##[45]##. Embryonic vascularization and angiogenic sprouting are severely deficient in CD148 −/− knock-out mice ##REF##12588999##[45]##, further supporting our hypothesis that VE-cadherin's contact-dependent intercellular signaling is crucial to vasculogenesis and angiogenesis.</p> <p>Perryn et al. ##REF##18062955##[19]## showed that anti-VE-cadherin treatment reduced sprout extension in murine allantois cultures by 70%, while it reduced cell-autonomous motility along sprout segments by 50%. Based on these results, they postulated that VE-cadherin is required for the motility of ECs along sprouts towards the tip. However, our simulations show that the observed cell slow-down after anti-VE-cadherin administration may be an indirect effect of a reduction of sprouting. Furthermore, our simulations suggest that even substantially reduced cell motility may not prevent patterning, though it does slow it down.</p> <p>In our simulations, branching and pattern formation require only experimentally-observed cell-level mechanisms, instead of the blood-vessel-level phenomenology in some other angiogenesis models ##UREF##4##[28]##–##UREF##5##[30]##. However, by starting with a cluster of endothelial cells, our simulations ignore many events preceding sprout formation, including the release of plasma proteins by the vessel, the breakdown of the basal lamina, the detachment of the ECs from surrounding ECs and smooth muscle cells, and cell proliferation. They also ignore subsequent processes consolidating outgrowth of the sprout, including tip-cell selection, any long-range chemoattractants and chemorepellants that guide the vessel to its target, the formation of new basal lamina, the sprout's association with stabilizing cells including pericytes, lumen formation within the sprout, and flow-induced remodeling of the developed vasculature. The mechanism for sprouting and network formation we have proposed forms a firm basis for future, more complete models of angiogenesis which include basal lamina and pericytes. We are currently studying the formation of directed sprouts with proliferating ECs in response to additional chemoattractants or chemorepellants and analyzing the role of cell elongation during sprouting. We are also studying the effect of additional, cell-cell contact-dependent signaling mechanisms, including delta-notch tip-cell selection ##REF##17259973##[3]## and chemoattractant-response amplifying Eph receptor-ephrin ligand interactions ##REF##16840724##[8]##.</p>
<p><bold>¤a:</bold> Current address: Netherlands Institute for Systems Biology and CWI (Center for Mathematics and Computer Science), Amsterdam, The Netherlands</p> <p><bold>¤b:</bold> Current address: Krumlauf Laboratory, Stowers Institute for Medical Research, Kansas City, Missouri, United States of America</p> <p>Conceived and designed the experiments: RMHM EDP JG. Performed the experiments: RMHM EDP. Analyzed the data: RMHM AS. Contributed reagents/materials/analysis tools: RMHM EDP AS. Wrote the paper: RMHM JG.</p> <p>Blood vessels form either when dispersed <italic>endothelial</italic> cells (the cells lining the inner walls of fully formed blood vessels) organize into a vessel network (<italic>vasculogenesis</italic>), or by sprouting or splitting of existing blood vessels (<italic>angiogenesis</italic>). Although they are closely related biologically, no current model explains both phenomena with a single biophysical mechanism. Most computational models describe sprouting at the level of the blood vessel, ignoring how cell behavior drives branch splitting during sprouting. We present a cell-based, Glazier–Graner–Hogeweg model (also called Cellular Potts Model) simulation of the initial patterning before the vascular cords form lumens, based on plausible behaviors of endothelial cells. The endothelial cells secrete a chemoattractant, which attracts other endothelial cells. As in the classic Keller–Segel model, chemotaxis by itself causes cells to aggregate into isolated clusters. However, including experimentally observed VE-cadherin–mediated contact inhibition of chemotaxis in the simulation causes randomly distributed cells to organize into networks and cell aggregates to sprout, reproducing aspects of both de novo and sprouting blood-vessel growth. We discuss two branching instabilities responsible for our results. Cells at the surfaces of cell clusters attempting to migrate to the centers of the clusters produce a buckling instability. In a model variant that eliminates the surface–normal force, a dissipative mechanism drives sprouting, with the secreted chemical acting both as a chemoattractant and as an inhibitor of pseudopod extension. Both mechanisms would also apply if force transmission through the extracellular matrix rather than chemical signaling mediated cell–cell interactions. The branching instabilities responsible for our results, which result from contact inhibition of chemotaxis, are both generic developmental mechanisms and interesting examples of unusual patterning instabilities.</p> <title>Author Summary</title> <p>A better understanding of the mechanisms by which <italic>endothelial</italic> cells (the cells lining the inner walls of blood vessels) organize into blood vessels is crucial if we need to enhance or suppress blood vessel growth under pathological conditions, including diabetes, wound healing, and tumor growth. During embryonic development, endothelial cells initially self-organize into a network of solid cords via blood vessel growth. The vascular network expands by splitting of existing blood vessels and by sprouting. Using computer simulations, we have captured a small set of biologically plausible cell behaviors that can reproduce the initial self-organization of endothelial cells, the sprouting of existing vessels, and the immediately subsequent remodeling of the resulting networks. In this model, endothelial cells both secrete diffusible <italic>chemoattractants</italic> and move up gradients of those chemicals by extending and retracting small <italic>pseudopods</italic>. By itself, this behavior causes simulated cells to accumulate to aggregate into large, round clusters. We propose that endothelial cells stop extending pseudopods along a given section of cell membrane as soon as the membrane touches the membrane of another endothelial cell (<italic>contact inhibition</italic>). Adding such contact-inhibition to our simulations allows vascular cords to form sprouts under a wide range of conditions.</p>
<title>Supporting Information</title>
<p>We gratefully acknowledge discussions with Stuart A. Newman, András Czirók, and Charles D. Little, and helpful comments of the anonymous referees. The simulation software is based on codes RM developed during work with Paulien Hogeweg at Utrecht University, The Netherlands.</p>
<fig id="pcbi-1000163-g001" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1000163.g001</object-id><label>Figure 1</label><caption><title>Anti-VE-cadherin antibody treatment inhibits de novo blood-vessel growth in mouse allantois cultures.</title><p>Endothelial cells fluorescently labeled in red with endothelium-specific CD34-Cy3 antibody. DIC/fluorescent image overlays. (A–C) Control. (D–F) Anti-VE-cadherin-treated cell cultures.</p></caption></fig> <fig id="pcbi-1000163-g002" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1000163.g002</object-id><label>Figure 2</label><caption><title>Endothelial cell aggregation; simulation initiated with 1000 scattered cells.</title><p>(A) After 10 Monte Carlo steps (<italic>MCS</italic>) (∼5 min). (B) After 1000 MCS (∼8 h). (C) After 10,000 MCS (∼80 h). (D) Contact-inhibited chemotaxis drives formation of vascular networks. Scale bar: 50 lattice sites (≈100 µm). Contour levels (green) indicate ten chemoattractant levels relative to the maximum concentration in the simulation. Grey shading indicates absolute concentration on a saturating scale.</p></caption></fig> <fig id="pcbi-1000163-g003" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1000163.g003</object-id><label>Figure 3</label><caption><title>Sprout formation in the absence of contact inhibition.</title><p>(A–C) Cell-autonomous cell elongation; (A) <italic>Λ</italic> = 22 µm; (B) <italic>Λ</italic> = 24 µm; (C) <italic>Λ</italic> = 32 µm; (D–F) Adhesion-driven sprouting. (D) <italic>J</italic>(<italic>c</italic>,<italic>c</italic>) = 1; (E) <italic>J</italic>(<italic>c</italic>,<italic>c</italic>) = 5; (F) <italic>J</italic>(<italic>c</italic>,<italic>c</italic>) = 10; (G–I) Passive cell elongation at short diffusion lengths; (G) <italic>D</italic> = 1·10<sup>−14</sup> m<sup>2</sup> s<sup>−1</sup>; (H) <italic>D</italic> = 2·10<sup>−14</sup> m<sup>2</sup> s<sup>−1</sup>; (I) <italic>D</italic> = 3·10<sup>−14</sup> m<sup>2</sup> s<sup>−1</sup>.</p></caption></fig> <fig id="pcbi-1000163-g004" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1000163.g004</object-id><label>Figure 4</label><caption><title>Sprouting instability in a simulation initiated with a cluster of endothelial cells.</title><p>(A) After 10 MCS (5 min). (B) After 1,000 MCS (∼8 h). (C) After 10,000 MCS (∼80 h). (D) No sprouting in a simulation without contact inhibition of chemotaxis (<italic>Χ</italic>(<italic>c</italic>,<italic>c</italic>)/<italic>Χ</italic>(<italic>c</italic>,<italic>M</italic>) = 1) at 10,000 MCS (∼80 h). Scale bar: 50 lattice sites (≈100 µm).</p></caption></fig> <fig id="pcbi-1000163-g005" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1000163.g005</object-id><label>Figure 5</label><caption><title>Compactness (<italic>C</italic> = <italic>A</italic>
<sub>cluster</sub>/<italic>A</italic>
<sub>hull</sub>) of 128-cell clusters after 10,000 MCS (∼80 h) as a function of the relative chemotactic response at cell-cell <italic>vs.</italic> cell-ECM interfaces.</title><p>Error bars show standard deviations over ten simulations.</p></caption></fig> <fig id="pcbi-1000163-g006" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1000163.g006</object-id><label>Figure 6</label><caption><title>Cell trajectories of simulated endothelial cells in 128-cell clusters in a contact-inhibited, sprouting cluster (A,B) and in a non-contact-inhibited, non-sprouting cluster (C,D).</title><p>(A,C) Cell trajectories during initial sprouting, indicating the cells' centers of mass at 100 MCS (∼50 min) intervals from 100 to 5,000 MCS (∼1–40 h). (B,D) Cell trajectories after initial sprouting, indicating the cells' centers of mass at 1,000 MCS (∼8 h) intervals from 4,000 to 20,000 MCS (∼30–170 h). Closed circles indicate initial cell positions; open circles indicate final cell positions. Colors identify individual cells; brightness increases from dark (initial positions) to bright (final positions). Outlines of clusters shown at 1000 MCS (∼8 h) intervals (A,C) or 4000 MCS (∼33 h) intervals (B,D). (E) Average displacement of cells from original positions over time in 10 simulations with 128 cells each, in contact-inhibited (solid curves) and non-contact-inhibited simulations (dashed curves). Grey curves indicate standard deviations. (F) Cell velocity 
##REF##18062955##[19]## with Δ<italic>t</italic> = 300 MCS (∼2.5 h) for contact-inhibited (solid curves) and non-contact-inhibited (dashed curves) simulations.</p></caption></fig> <fig id="pcbi-1000163-g007" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1000163.g007</object-id><label>Figure 7</label><caption><title>Compactness (<italic>C</italic> = <italic>A</italic>
<sub>cluster</sub>/<italic>A</italic>
<sub>hull</sub>) of 128-cell clusters on 200×200-pixel lattices (∼400 µm×400 µm) after 5,000 MCS (∼40 h) for standard chemotaxis, as a function of the adhesion between endothelial cells, <italic>J</italic>(<italic>c</italic>,<italic>c</italic>).</title><p>For <italic>J</italic>(<italic>c</italic>,<italic>c</italic>)&lt;40 (i.e., <italic>J</italic>(<italic>c</italic>,<italic>c</italic>)&lt;2 <italic>J</italic>(<italic>c</italic>,<italic>M</italic>) the cells adhere without chemotaxis. Insets: Representative configurations after 5000 MCS (∼40 h).</p></caption></fig> <fig id="pcbi-1000163-g008" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1000163.g008</object-id><label>Figure 8</label><caption><title>Compactness (<italic>C</italic> = <italic>A</italic>
<sub>cluster</sub>/<italic>A</italic>
<sub>hull</sub>) of 128-cell clusters on 200×200-pixel lattices (∼400 µm×400 µm ) after 5000 MCS (∼40 h) for standard chemotaxis as a function of absolute chemotactic strength, <italic>Χ</italic>(<italic>c</italic>,<italic>M</italic>).</title><p>Insets: Representative configurations after 5000 MCS (∼40 h).</p></caption></fig> <fig id="pcbi-1000163-g009" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1000163.g009</object-id><label>Figure 9</label><caption><title>Compactness (<italic>C</italic> = <italic>A</italic>
<sub>cluster</sub>/<italic>A</italic>
<sub>hull</sub>) of 128-cell clusters on 200×200-lattices (∼400 µm×400 µm) after 5000 MCS (∼40 h) for standard chemotaxis as a function of the saturation of the chemotactic response, <italic>s</italic>.</title><p>Insets: representative configurations after 5000 MCS (∼40 h).</p></caption></fig> <fig id="pcbi-1000163-g010" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1000163.g010</object-id><label>Figure 10</label><caption><title>Compactness (<italic>C</italic> = <italic>A</italic>
<sub>cluster</sub>/<italic>A</italic>
<sub>hull</sub>) of 128-cell clusters (solid curve) on 200×200-pixel lattices (∼400 µm×400 µm) and 1,024-cell clusters (dashed-dotted curve) on 400×400-pixel lattices (∼800 µm×800 µm) after 5000 MCS (∼40 h) for Savill-Hogeweg chemotaxis as a function of the chemoattractant diffusion constant <italic>D</italic>.</title><p>Larger diffusion constants have longer diffusion lengths, . Dashed curve shows the compactness of VE-cadherin-inhibited 128-cell clusters. Insets: Representative configurations after 5,000 MCS (∼40 h) of the 128-cell clusters (left panels) and 1,024-cell clusters (right panels; not to scale).</p></caption></fig> <fig id="pcbi-1000163-g011" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1000163.g011</object-id><label>Figure 11</label><caption><title>Compactness (<italic>C</italic> = <italic>A</italic>
<sub>cluster</sub>/<italic>A</italic>
<sub>hull</sub>) of 128-cell clusters on 400×400-pixel lattices (∼400 µm×400 µm) after 5,000 MCS (∼40 h) as a function of the cell motility <italic>T</italic>, for standard Savill–Hogeweg ##UREF##8##[36]##
<italic>extension–retraction</italic> chemotaxis (solid curve), and for <italic>extension-only</italic> chemotaxis (dashed curve).</title><p>Black curves show the mean over 100 simulations for each <italic>T</italic> (with a <italic>T</italic>-increment of 10). Dotted grey curves indicate one standard deviation. Insets: Representative configurations after 5,000 MCS (∼40 h).</p></caption></fig> <fig id="pcbi-1000163-g012" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1000163.g012</object-id><label>Figure 12</label><caption><title>Evolution of the compactness (<italic>C</italic> = <italic>A</italic>
<sub>cluster</sub>/<italic>A</italic>
<sub>hull</sub>) of 256-cell clusters on 500×500-pixel lattices (∼1,000 µm×1,000 µm) <italic>vs.</italic> time for standard Savill-Hogeweg ##UREF##8##[36]##
<italic>extension-retraction</italic> chemotaxis (solid and dashed curves, for <italic>T</italic> = 50 and <italic>T</italic> = 200 respectively), and for <italic>extension-only</italic> chemotaxis (dash-dotted curve, <italic>T</italic> = 200), with only extending pseudopods responding to the chemoattractant.</title><p>Black curves show the mean of 100 simulations. Dotted grey curves mark one standard deviation. Insets: Representative configurations after 1000 (∼8 h) and 5000 MCS (∼40 h). Videos available online.</p></caption></fig> <fig id="pcbi-1000163-g013" position="float"><object-id pub-id-type="doi">10.1371/journal.pcbi.1000163.g013</object-id><label>Figure 13</label><caption><title>Cumulative energy differences for standard Savill–Hogeweg ##UREF##8##[36]##
<italic>extension</italic>–<italic>retraction</italic> chemotaxis (solid and dashed curves, for <italic>T</italic> = 50 and <italic>T</italic> = 200, respectively), and for <italic>extension-only</italic> chemotaxis (dash-dotted curve, <italic>T</italic> = 200) as a function of time.</title><p>Black curves show the mean of 100 simulations. Dotted grey curves mark one standard deviation.</p></caption></fig>
<disp-formula><label>(1)</label></disp-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <disp-formula><label>(2)</label></disp-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <disp-formula><label>(3)</label></disp-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <disp-formula><label>(4)</label></disp-formula>
<supplementary-material content-type="local-data" id="pcbi.1000163.s001"><label>Dataset S1</label><caption><p>Parameter files for the simulations shown in ##FIG##1##Figures 2##, ##FIG##3##4##, ##FIG##10##11##, and ##FIG##11##12##, packed as a tar.gz archive To use, unpack the parameter-file archive and install the Tissue Simulation Toolkit (##SUPPL##1##Protocol S1##). Run the simulations from the command line using the command “vessel [parameter-file]”. Reproduce the other simulations by editing the parameter files using a standard text editor to set the values specified in the text.</p><p>(4 KB ZIP)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pcbi.1000163.s002"><label>Protocol S1</label><caption><p>Tissue Simulation Toolkit v0.1.3. The source code for the software used for the simulations presented in this paper is also available from <ext-link ext-link-type="uri" xlink:href="http://sourceforge.net/projects/tst">http://sourceforge.net/projects/tst</ext-link>. Installation: Unpack and compile according to the instructions given in the INSTALL file The code is written in C++ using the cross-platform (Windows, Mac, or Unix/Linux) library Qt (available from <ext-link ext-link-type="uri" xlink:href="http://www.trolltech.com">www.trolltech.com</ext-link>).</p><p>(332 KB ZIP)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pcbi.1000163.s003"><label>Figure S1</label><caption><p>Compactness (<italic>C</italic> = <italic>A</italic>
<sub>cluster</sub>/<italic>A</italic>
<sub>hull</sub>) of 128-cell clusters after 5000 MCS (∼40 h) as a function of the cell length, in the presence (solid curve) or absence (dashed curve) of contact inhibition. Lengths given in terms of the <italic>target length Λ</italic> as defined in Merks et al. 2006 ##REF##16325173##[22]##. Grey curves show standard deviations over ten simulations.</p><p>(1278 KB EPS)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pcbi.1000163.s004"><label>Figure S2</label><caption><p>Configurations of 128-cell clusters after 5,000 MCS (∼40 h) for increasing chemoattractant secretion rates for low (<italic>s</italic> = 0.01) and high (<italic>s</italic> = 0.05) chemoattractant receptor saturations.</p><p>(3536 KB EPS)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pcbi.1000163.s005"><label>Video S1</label><caption><p>Endothelial cell aggregation without contact-inhibited chemotaxis. Simulation initiated with 1,000 scattered cells. 0 MCS to 20,000 MCS (∼0–170 h), 100 MCS per frame.</p><p>(1902 KB AVI)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pcbi.1000163.s006"><label>Video S2</label><caption><p>Endothelial cell aggregation in with contact-inhibited chemotaxis. Simulation initiated with 1000 scattered cells. 0 MCS to 2000 MCS (∼0–20 h), 10 MCS per frame.</p><p>(4534 KB AVI)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pcbi.1000163.s007"><label>Video S3</label><caption><p>Same simulation as ##SUPPL##5##Video S2##. 0 MCS to 20,000 MCS (∼0–170 h), 100 MCS per frame.</p><p>(2000 KB AVI)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pcbi.1000163.s008"><label>Video S4</label><caption><p>Sprouting instability in a simulation with contact-inhibited chemotaxis, initiated with a cluster of 256 endothelial cells. 0 MCS to 2,000 MCS (∼0–20 h), 10 MCS per frame.</p><p>(4554 KB AVI)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pcbi.1000163.s009"><label>Video S5</label><caption><p>Same simulation as ##SUPPL##7##Video S4##. 0 MCS to 20,000 MCS (∼0–170 h), 100 MCS per frame.</p><p>(1895 KB AVI)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pcbi.1000163.s010"><label>Video S6</label><caption><p>Simulation with non-contact–inhibited chemotaxis, initiated with a cluster of 256 endothelial cells. 0 MCS to 2,000 MCS (∼0–20 h), 10 MCS per frame.</p><p>(1640 KB AVI)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pcbi.1000163.s011"><label>Video S7</label><caption><p>Effect of cell adhesion on the proposed sprouting-angiogenesis mechanism (<italic>cf.</italic>
##FIG##6##Figure 7##). Simulations with contact-inhibited chemotaxis, initiated with a cluster of 256 endothelial cells. 0 MCS to 20,000 MCS (∼0–170 h), 100 MCS per frame. <italic>J</italic>(<italic>c</italic>,<italic>c</italic>) = 5.</p><p>(699 KB AVI)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pcbi.1000163.s012"><label>Video S8</label><caption><p>Effect of cell adhesion on the proposed sprouting-angiogenesis mechanism (<italic>cf.</italic>
##FIG##6##Figure 7##). Simulations with contact-inhibited chemotaxis, initiated with a cluster of 256 endothelial cells. 0 MCS to 20,000 MCS (∼0–170 h), 100 MCS per frame. <italic>J</italic>(<italic>c</italic>,<italic>c</italic>) = 10.</p><p>(685 KB AVI)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pcbi.1000163.s013"><label>Video S9</label><caption><p>Effect of cell adhesion on the proposed sprouting-angiogenesis mechanism (<italic>cf.</italic>
##FIG##6##Figure 7##). Simulations with contact-inhibited chemotaxis, initiated with a cluster of 256 endothelial cells. 0 MCS to 20,000 MCS (∼0–170 h), 100 MCS per frame. <italic>J</italic>(<italic>c</italic>,<italic>c</italic>) = 15.</p><p>(691 KB AVI)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pcbi.1000163.s014"><label>Video S10</label><caption><p>Effect of cell adhesion on the proposed sprouting-angiogenesis mechanism (<italic>cf.</italic>
##FIG##6##Figure 7##). Simulations with contact-inhibited chemotaxis, initiated with a cluster of 256 endothelial cells. 0 MCS to 20,000 MCS (∼0–170 h), 100 MCS per frame. <italic>J</italic>(<italic>c</italic>,<italic>c</italic>) = 20.</p><p>(669 KB AVI)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pcbi.1000163.s015"><label>Video S11</label><caption><p>Effect of cell adhesion on the proposed sprouting-angiogenesis mechanism (<italic>cf.</italic>
##FIG##6##Figure 7##). Simulations with contact-inhibited chemotaxis, initiated with a cluster of 256 endothelial cells. 0 MCS to 20,000 MCS (∼0–170 h), 100 MCS per frame. <italic>J</italic>(<italic>c</italic>,<italic>c</italic>) = 25.</p><p>(677 KB AVI)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pcbi.1000163.s016"><label>Video S12</label><caption><p>Effect of cell adhesion on the proposed sprouting-angiogenesis mechanism (<italic>cf.</italic>
##FIG##6##Figure 7##). Simulations with contact-inhibited chemotaxis, initiated with a cluster of 256 endothelial cells. 0 MCS to 20,000 MCS (∼0–170 h), 100 MCS per frame. <italic>J</italic>(<italic>c</italic>,<italic>c</italic>) = 30.</p><p>(669 KB AVI)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pcbi.1000163.s017"><label>Video S13</label><caption><p>Effect of cell adhesion on the proposed sprouting-angiogenesis mechanism (<italic>cf.</italic>
##FIG##6##Figure 7##). Simulations with contact-inhibited chemotaxis, initiated with a cluster of 256 endothelial cells. 0 MCS to 20,000 MCS (∼0–170 h), 100 MCS per frame. <italic>J</italic>(<italic>c</italic>,<italic>c</italic>) = 35.</p><p>(679 KB AVI)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pcbi.1000163.s018"><label>Video S14</label><caption><p>Effect of cell adhesion on the proposed sprouting-angiogenesis mechanism (<italic>cf.</italic>
##FIG##6##Figure 7##). Simulations with contact-inhibited chemotaxis, initiated with a cluster of 256 endothelial cells. 0 MCS to 20,000 MCS (∼0–170 h), 100 MCS per frame. <italic>J</italic>(<italic>c</italic>,<italic>c</italic>) = 40.</p><p>(648 KB AVI)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pcbi.1000163.s019"><label>Video S15</label><caption><p>Effect of cell adhesion on the proposed sprouting-angiogenesis mechanism (<italic>cf.</italic>
##FIG##6##Figure 7##). Simulations with contact-inhibited chemotaxis, initiated with a cluster of 256 endothelial cells. 0 MCS to 20,000 MCS (∼0–170 h), 100 MCS per frame. <italic>J</italic>(<italic>c</italic>,<italic>c</italic>) = 45.</p><p>(383 KB AVI)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pcbi.1000163.s020"><label>Video S16</label><caption><p>Effect of cell adhesion on the proposed sprouting-angiogenesis mechanism (<italic>cf.</italic>
##FIG##6##Figure 7##). Simulations with contact-inhibited chemotaxis, initiated with a cluster of 256 endothelial cells. 0 MCS to 20,000 MCS (∼0–170 h), 100 MCS per frame. <italic>J</italic>(<italic>c,c</italic>) = 50.</p><p>(557 KB AVI)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pcbi.1000163.s021"><label>Video S17</label><caption><p>Effect of cell adhesion on the proposed sprouting-angiogenesis mechanism (<italic>cf.</italic>
##FIG##6##Figure 7##). Simulations with contact-inhibited chemotaxis, initiated with a cluster of 256 endothelial cells. 0 MCS to 20,000 MCS (∼0–170 h), 100 MCS per frame. <italic>J</italic>(<italic>c</italic>,<italic>c</italic>) = 55.</p><p>(528 KB AVI)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pcbi.1000163.s022"><label>Video S18</label><caption><p>Effect of cell adhesion on the proposed sprouting-angiogenesis mechanism (<italic>cf.</italic>
##FIG##6##Figure 7##). Simulations with contact-inhibited chemotaxis, initiated with a cluster of 256 endothelial cells. 0 MCS to 20,000 MCS (∼0–170 h), 100 MCS per frame. <italic>J</italic>(<italic>c</italic>,<italic>c</italic>) = 60.</p><p>(516 KB AVI)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pcbi.1000163.s023"><label>Video S19</label><caption><p>Effect of cell adhesion on the proposed sprouting-angiogenesis mechanism (<italic>cf.</italic>
##FIG##6##Figure 7##). Simulations with contact-inhibited chemotaxis, initiated with a cluster of 256 endothelial cells. 0 MCS to 20,000 MCS (∼0–170 h), 100 MCS per frame. <italic>J</italic>(<italic>c</italic>,<italic>c</italic>) = 65.</p><p>(457 KB AVI)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pcbi.1000163.s024"><label>Video S20</label><caption><p>Effect of cell adhesion on the proposed sprouting-angiogenesis mechanism (<italic>cf.</italic>
##FIG##6##Figure 7##). Simulations with contact-inhibited chemotaxis, initiated with a cluster of 256 endothelial cells. 0 MCS to 20,000 MCS (∼0–170 h), 100 MCS per frame. <italic>J</italic>(<italic>c,c</italic>) = 70.</p><p>(439 KB AVI)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pcbi.1000163.s025"><label>Video S21</label><caption><p>Effect of cell adhesion on the proposed sprouting-angiogenesis mechanism (<italic>cf.</italic>
##FIG##6##Figure 7##). Simulations with contact-inhibited chemotaxis, initiated with a cluster of 256 endothelial cells. 0 MCS to 20,000 MCS (∼0–170 h), 100 MCS per frame. <italic>J</italic>(<italic>c</italic>,<italic>c</italic>) = 75.</p><p>(411KB AVI)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pcbi.1000163.s026"><label>Video S22</label><caption><p>Effect of cell adhesion on the proposed sprouting-angiogenesis mechanism (<italic>cf.</italic>
##FIG##6##Figure 7##). Simulations with contact-inhibited chemotaxis, initiated with a cluster of 256 endothelial cells. 0 MCS to 20,000 MCS (∼0–170 h), 100 MCS per frame. <italic>J</italic>(<italic>c</italic>,<italic>c</italic>) = 80.</p><p>(383 KB AVI)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pcbi.1000163.s027"><label>Video S23</label><caption><p>Effect of the chemotactic strength on the proposed sprouting-angiogenesis mechanism (<italic>cf.</italic>
##FIG##7##Figure 8##). Simulations with contact-inhibited chemotaxis, initiated with a cluster of 256 endothelial cells. 0 MCS to 20,000 MCS (∼0–170 h), 100 MCS per frame. <italic>Χ</italic>(<italic>c</italic>,<italic>M</italic>) = 0.</p><p>(238 KB AVI)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pcbi.1000163.s028"><label>Video S24</label><caption><p>Effect of the chemotactic strength on the proposed sprouting-angiogenesis mechanism (<italic>cf.</italic>
##FIG##7##Figure 8##). Simulations with contact-inhibited chemotaxis, initiated with a cluster of 256 endothelial cells. 0 MCS to 20,000 MCS (∼0–170 h), 100 MCS per frame. <italic>Χ</italic>(<italic>c</italic>,<italic>M</italic>) = 500.</p><p>(653 KB AVI)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pcbi.1000163.s029"><label>Video S25</label><caption><p>Effect of the chemotactic strength on the proposed sprouting-angiogenesis mechanism (<italic>cf.</italic>
##FIG##7##Figure 8##). Simulations with contact-inhibited chemotaxis, initiated with a cluster of 256 endothelial cells. 0 MCS to 20,000 MCS (∼0–170 h), 100 MCS per frame. <italic>Χ</italic>(<italic>c</italic>,<italic>M</italic>) = 1,000.</p><p>(707 KB AVI)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pcbi.1000163.s030"><label>Video S26</label><caption><p>Effect of the chemotactic strength on the proposed sprouting-angiogenesis mechanism (<italic>cf.</italic>
##FIG##7##Figure 8##). Simulations with contact-inhibited chemotaxis, initiated with a cluster of 256 endothelial cells. 0 MCS to 20,000 MCS (∼0–170 h), 100 MCS per frame. <italic>Χ</italic>(<italic>c</italic>,<italic>M</italic>) = 1,500.</p><p>(714 KB AVI)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pcbi.1000163.s031"><label>Video S27</label><caption><p>Effect of the chemotactic strength on the proposed sprouting-angiogenesis mechanism (<italic>cf.</italic>
##FIG##7##Figure 8##). Simulations with contact-inhibited chemotaxis, initiated with a cluster of 256 endothelial cells. 0 MCS to 20,000 MCS (∼0–170 h), 100 MCS per frame. <italic>Χ</italic>(<italic>c</italic>,<italic>M</italic>) = 2,000.</p><p>(727 KB AVI)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pcbi.1000163.s032"><label>Video S28</label><caption><p>Effect of the chemotactic strength on the proposed sprouting-angiogenesis mechanism (<italic>cf.</italic>
##FIG##7##Figure 8##). Simulations with contact-inhibited chemotaxis, initiated with a cluster of 256 endothelial cells. 0 MCS to 20,000 MCS (∼0–170 h), 100 MCS per frame. <italic>Χ</italic>(<italic>c</italic>,<italic>M</italic>) = 2,500.</p><p>(765 KB AVI)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pcbi.1000163.s033"><label>Video S29</label><caption><p>Effect of the chemotactic strength on the proposed sprouting-angiogenesis mechanism (<italic>cf.</italic>
##FIG##7##Figure 8##). Simulations with contact-inhibited chemotaxis, initiated with a cluster of 256 endothelial cells. 0 MCS to 20,000 MCS (∼0–170 h), 100 MCS per frame. <italic>Χ</italic>(<italic>c</italic>,<italic>M</italic>) = 3,000.</p><p>(797 KB AVI)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pcbi.1000163.s034"><label>Video S30</label><caption><p>Effect of the chemotactic strength on the proposed sprouting-angiogenesis mechanism (<italic>cf.</italic>
##FIG##7##Figure 8##). Simulations with contact-inhibited chemotaxis, initiated with a cluster of 256 endothelial cells. 0 MCS to 20,000 MCS (∼0–170 h), 100 MCS per frame. <italic>Χ</italic>(<italic>c</italic>,<italic>M</italic>) = 3,500.</p><p>(805 KB AVI)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pcbi.1000163.s035"><label>Video S31</label><caption><p>Effect of the chemotactic strength on the proposed sprouting-angiogenesis mechanism (<italic>cf.</italic>
##FIG##7##Figure 8##). Simulations with contact-inhibited chemotaxis, initiated with a cluster of 256 endothelial cells. 0 MCS to 20,000 MCS (∼0–170 h), 100 MCS per frame. <italic>Χ</italic>(<italic>c</italic>,<italic>M</italic>) = 4,000.</p><p>(835 KB AVI)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pcbi.1000163.s036"><label>Video S32</label><caption><p>Effect of the chemotactic strength on the proposed sprouting-angiogenesis mechanism (<italic>cf.</italic>
##FIG##7##Figure 8##). Simulations with contact-inhibited chemotaxis, initiated with a cluster of 256 endothelial cells. 0 MCS to 20,000 MCS (∼0–170 h), 100 MCS per frame. <italic>Χ</italic>(<italic>c</italic>,<italic>M</italic>) = 4,500.</p><p>(853 KB AVI)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pcbi.1000163.s037"><label>Video S33</label><caption><p>Effect of the chemotactic strength on the proposed sprouting-angiogenesis mechanism (<italic>cf.</italic>
##FIG##7##Figure 8##). Simulations with contact-inhibited chemotaxis, initiated with a cluster of 256 endothelial cells. 0 MCS to 20,000 MCS (∼0–170 h), 100 MCS per frame. <italic>Χ</italic>(<italic>c</italic>,<italic>M</italic>) = 5,000.</p><p>(858 KB AVI)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pcbi.1000163.s038"><label>Video S34</label><caption><p>Effect of chemotaxis saturation on the proposed sprouting-angiogenesis mechanism (<italic>cf.</italic>
##FIG##8##Figure 9##). Simulations with contact-inhibited chemotaxis, initiated with a cluster of 256 endothelial cells. 0 MCS to 20,000 MCS (∼0–170 h), 100 MCS per frame. <italic>s</italic> = 0.0.</p><p>(655 KB AVI)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pcbi.1000163.s039"><label>Video S35</label><caption><p>Effect of chemotaxis saturation on the proposed sprouting-angiogenesis mechanism (<italic>cf.</italic>
##FIG##8##Figure 9##). Simulations with contact-inhibited chemotaxis, initiated with a cluster of 256 endothelial cells. 0 MCS to 20,000 MCS (∼0–170 h), 100 MCS per frame. <italic>s</italic> = 0.05.</p><p>(577 KB AVI)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pcbi.1000163.s040"><label>Video S36</label><caption><p>Effect of chemotaxis saturation on the proposed sprouting-angiogenesis mechanism (<italic>cf.</italic>
##FIG##8##Figure 9##). Simulations with contact-inhibited chemotaxis, initiated with a cluster of 256 endothelial cells. 0 MCS to 20,000 MCS (∼0–170 h), 100 MCS per frame. <italic>s</italic> = 0.1.</p><p>(409 KB AVI)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pcbi.1000163.s041"><label>Video S37</label><caption><p>Effect of chemotaxis saturation on the proposed sprouting-angiogenesis mechanism (<italic>cf.</italic>
##FIG##8##Figure 9##). Simulations with contact-inhibited chemotaxis, initiated with a cluster of 256 endothelial cells. 0 MCS to 20,000 MCS (∼0–170 h), 100 MCS per frame. <italic>s</italic> = 0.15.</p><p>(224 KB AVI)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pcbi.1000163.s042"><label>Video S38</label><caption><p>Effect of chemotaxis saturation on the proposed sprouting-angiogenesis mechanism (<italic>cf.</italic> Figure 39). Simulations with contact-inhibited chemotaxis, initiated with a cluster of 256 endothelial cells. 0 MCS to 20,000 MCS (∼0–170 h), 100 MCS per frame. <italic>s</italic> = 0.2.</p><p>(213 KB AVI)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pcbi.1000163.s043"><label>Video S39</label><caption><p>Effect of chemotaxis saturation on the proposed sprouting-angiogenesis mechanism (<italic>cf.</italic>
##FIG##8##Figure 9##). Simulations with contact-inhibited chemotaxis, initiated with a cluster of 256 endothelial cells. 0 MCS to 20,000 MCS (∼0–170 h), 100 MCS per frame. <italic>s</italic> = 0.25.</p><p>(215 KB AVI)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pcbi.1000163.s044"><label>Video S40</label><caption><p>Effect of the diffusion constant <italic>D</italic> on the proposed sprouting-angiogenesis mechanism (<italic>cf.</italic>
##FIG##9##Figure 10##). Simulations with contact-inhibited chemotaxis, initiated with a cluster of 256 endothelial cells. 0 MCS to 20,000 MCS (∼0–170 h), 100 MCS per frame. <italic>D</italic> = 4·10<sup>−14</sup> m<sup>2</sup> s<sup>−1</sup>.</p><p>(1053 KB AVI)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pcbi.1000163.s045"><label>Video S41</label><caption><p>Effect of the diffusion constant <italic>D</italic> on the proposed sprouting-angiogenesis mechanism (<italic>cf.</italic>
##FIG##9##Figure 10##). Simulations with contact-inhibited chemotaxis, initiated with a cluster of 256 endothelial cells. 0 MCS to 20,000 MCS (∼0–170 h), 100 MCS per frame. <italic>D</italic> = 5·10<sup>−14</sup> m<sup>2</sup> s<sup>−1</sup>.</p><p>(717 KB AVI)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pcbi.1000163.s046"><label>Video S42</label><caption><p>Effect of the diffusion constant <italic>D</italic> on the proposed sprouting-angiogenesis mechanism (<italic>cf.</italic>
##FIG##9##Figure 10##). Simulations with contact-inhibited chemotaxis, initiated with a cluster of 256 endothelial cells. 0 MCS to 20,000 MCS (∼0–170 h), 100 MCS per frame. <italic>D</italic> = 1·10<sup>−13</sup> m<sup>2</sup> s<sup>−1</sup>.</p><p>(631 KB AVI)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pcbi.1000163.s047"><label>Video S43</label><caption><p>Effect of the diffusion constant <italic>D</italic> on the proposed sprouting-angiogenesis mechanism (<italic>cf.</italic>
##FIG##9##Figure 10##). Simulations with contact-inhibited chemotaxis, initiated with a cluster of 256 endothelial cells. 0 MCS to 20,000 MCS (∼0–170 h), 100 MCS per frame. <italic>D</italic> = 2·10<sup>−13</sup> m<sup>2</sup> s<sup>−1</sup>.</p><p>(501 KB AVI)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pcbi.1000163.s048"><label>Video S44</label><caption><p>Effect of the diffusion constant <italic>D</italic> on the proposed sprouting-angiogenesis mechanism (<italic>cf.</italic>
##FIG##9##Figure 10##). Simulations with contact-inhibited chemotaxis, initiated with a cluster of 256 endothelial cells. 0 MCS to 20,000 MCS (∼0–170 h), 100 MCS per frame. <italic>D</italic> = 3·10<sup>−13</sup> m<sup>2</sup> s<sup>−1</sup>.</p><p>(418 KB AVI)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pcbi.1000163.s049"><label>Video S45</label><caption><p>Effect of the diffusion constant <italic>D</italic> on the proposed sprouting-angiogenesis mechanism (<italic>cf.</italic>
##FIG##9##Figure 10##). Simulations with contact-inhibited chemotaxis, initiated with a cluster of 256 endothelial cells. 0 MCS to 20,000 MCS (∼0–170 h), 100 MCS per frame. <italic>D</italic> = 4·10<sup>−13</sup> m<sup>2</sup> s<sup>−1</sup>.</p><p>(371 KB AVI)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pcbi.1000163.s050"><label>Video S46</label><caption><p>Effect of the diffusion constant <italic>D</italic> on the proposed sprouting-angiogenesis mechanism (<italic>cf.</italic>
##FIG##9##Figure 10##). Simulations with contact-inhibited chemotaxis, initiated with a cluster of 256 endothelial cells. 0 MCS to 20,000 MCS (∼0–170 h), 100 MCS per frame. <italic>D</italic> = 5·10<sup>−13</sup> m<sup>2</sup> s<sup>−1</sup>.</p><p>(316KB AVI)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pcbi.1000163.s051"><label>Video S47</label><caption><p>Effect of the diffusion constant <italic>D</italic> on the proposed sprouting-angiogenesis mechanism (<italic>cf.</italic>
##FIG##9##Figure 10##). Simulations with contact-inhibited chemotaxis, initiated with a cluster of 1024 endothelial cells. MCS 0 to 20,000 (∼0–170 h), 100 MCS per frame. <italic>D</italic> = 1·10<sup>−14</sup> m<sup>2</sup> s<sup>−1</sup>.</p><p>(1376 KB AVI)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pcbi.1000163.s052"><label>Video S48</label><caption><p>Effect of the diffusion constant <italic>D</italic> on the proposed sprouting-angiogenesis mechanism (<italic>cf.</italic>
##FIG##9##Figure 10##). Simulations with contact-inhibited chemotaxis, initiated with a cluster of 1024 endothelial cells. MCS 0 to 20,000 (∼0–170 h), 100 MCS per frame. <italic>D</italic> = 5·10<sup>−14</sup> m<sup>2</sup> s<sup>−1</sup>.</p><p>(1011 KB AVI)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pcbi.1000163.s053"><label>Video S49</label><caption><p>Effect of the diffusion constant <italic>D</italic> on the proposed sprouting-angiogenesis mechanism (<italic>cf.</italic>
##FIG##9##Figure 10##). Simulations with contact-inhibited chemotaxis, initiated with a cluster of 1024 endothelial cells. MCS 0 to 20,000 (∼0–170 h), 100 MCS per frame. <italic>D</italic> = 1·10<sup>−13</sup> m<sup>2</sup> s<sup>−1</sup>.</p><p>(1003 KB AVI)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pcbi.1000163.s054"><label>Video S50</label><caption><p>Effect of the diffusion constant <italic>D</italic> on the proposed sprouting-angiogenesis mechanism (<italic>cf.</italic>
##FIG##9##Figure 10##). Simulations with contact-inhibited chemotaxis, initiated with a cluster of 1024 endothelial cells. MCS 0 to 20,000 (∼0–170 h), 100 MCS per frame. <italic>D</italic> = 2·10<sup>−13</sup> m<sup>2</sup> s<sup>−1</sup>.</p><p>(994 KB AVI)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pcbi.1000163.s055"><label>Video S51</label><caption><p>Effect of the diffusion constant <italic>D</italic> on the proposed sprouting-angiogenesis mechanism (<italic>cf.</italic>
##FIG##9##Figure 10##). Simulations with contact-inhibited chemotaxis, initiated with a cluster of 1024 endothelial cells. MCS 0 to 20,000 (∼0–170 h), 100 MCS per frame. <italic>D</italic> = 3·10<sup>−13</sup> m<sup>2</sup> s<sup>−1</sup>.</p><p>(976 KB AVI)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pcbi.1000163.s056"><label>Video S52</label><caption><p>Effect of the diffusion constant <italic>D</italic> on the proposed sprouting-angiogenesis mechanism (<italic>cf.</italic>
##FIG##9##Figure 10##). Simulations with contact-inhibited chemotaxis, initiated with a cluster of 1024 endothelial cells. MCS 0 to 20,000 (∼0–170 h), 100 MCS per frame. <italic>D</italic> = 4·10<sup>−13</sup> m<sup>2</sup> s<sup>−1</sup>.</p><p>(978 KB AVI)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pcbi.1000163.s057"><label>Video S53</label><caption><p>Effect of the diffusion constant <italic>D</italic> on the proposed sprouting-angiogenesis mechanism (<italic>cf.</italic>
##FIG##9##Figure 10##). Simulations with contact-inhibited chemotaxis, initiated with a cluster of 1024 endothelial cells. MCS 0 to 20,000 (∼0–170 h), 100 MCS per frame. <italic>D</italic> = 5·10<sup>−13</sup> m<sup>2</sup> s<sup>−1</sup>.</p><p>(982 KB AVI)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pcbi.1000163.s058"><label>Video S54</label><caption><p>Sprouting of a 256-cell cluster on a 500×500-pixel lattice (∼1 mm×1 mm) with standard Savill–Hogeweg, extension-retraction chemotaxis at <italic>T</italic> = 50, as in ##FIG##11##Figure 12##.</p><p>(2366 KB AVI)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pcbi.1000163.s059"><label>Video S55</label><caption><p>Sprouting of a 256-cell cluster on a 500×500-pixel lattice (∼1 mm×1 mm) with standard Savill–Hogeweg, extension–retraction chemotaxis at <italic>T</italic> = 200, as in ##FIG##11##Figure 12##.</p><p>(2372 KB AVI)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pcbi.1000163.s060"><label>Video S56</label><caption><p>Sprouting of a 256-cell cluster on a 500×500-pixel lattice (∼1 mm×1 mm) with extension-only chemotaxis at <italic>T</italic> = 200, as in ##FIG##11##Figure 12##.</p><p>(2381 KB AVI)</p></caption></supplementary-material>
<fn-group><fn fn-type="COI-statement"><p>The authors have declared that no competing interests exist.</p></fn><fn fn-type="financial-disclosure"><p>This work received support from NIH grants NIGMS 1R01 GM077138-01A1 and 5R01 GM076692 (JAG; AS), the Indiana University Biocomplexity Institute (RM), the VIB Department of Plant Systems Biology (RM), the Indiana University AVIDD program (JAG; RM) and an American Heart Association predoctoral fellowship 0410084Z (EDP). None of these funding agencies has had any role in the design and conduct of the study, in the collection, analysis, and interpretation of the data, or in the preparation, review, or approval of the manuscript.</p></fn></fn-group>
<graphic xlink:href="pcbi.1000163.e001.jpg" mimetype="image" position="float"/> <inline-graphic xlink:href="pcbi.1000163.e002.jpg" mimetype="image"/> <inline-graphic xlink:href="pcbi.1000163.e003.jpg" mimetype="image"/> <graphic xlink:href="pcbi.1000163.g001"/> <graphic xlink:href="pcbi.1000163.g002"/> <inline-graphic xlink:href="pcbi.1000163.e004.jpg" mimetype="image"/> <graphic xlink:href="pcbi.1000163.g003"/> <graphic xlink:href="pcbi.1000163.g004"/> <graphic xlink:href="pcbi.1000163.g005"/> <inline-graphic xlink:href="pcbi.1000163.e005.jpg" mimetype="image"/> <inline-graphic xlink:href="pcbi.1000163.e006.jpg" mimetype="image"/> <graphic xlink:href="pcbi.1000163.g006"/> <graphic xlink:href="pcbi.1000163.g007"/> <graphic xlink:href="pcbi.1000163.g008"/> <graphic xlink:href="pcbi.1000163.g009"/> <inline-graphic xlink:href="pcbi.1000163.e007.jpg" mimetype="image"/> <graphic xlink:href="pcbi.1000163.g010"/> <graphic xlink:href="pcbi.1000163.g011"/> <graphic xlink:href="pcbi.1000163.g012"/> <graphic xlink:href="pcbi.1000163.g013"/> <inline-graphic xlink:href="pcbi.1000163.e008.jpg" mimetype="image"/> <inline-graphic xlink:href="pcbi.1000163.e009.jpg" mimetype="image"/> <inline-graphic xlink:href="pcbi.1000163.e010.jpg" mimetype="image"/> <graphic xlink:href="pcbi.1000163.e011.jpg" mimetype="image" position="float"/> <inline-graphic xlink:href="pcbi.1000163.e012.jpg" mimetype="image"/> <inline-graphic xlink:href="pcbi.1000163.e013.jpg" mimetype="image"/> <inline-graphic xlink:href="pcbi.1000163.e014.jpg" mimetype="image"/> <inline-graphic xlink:href="pcbi.1000163.e015.jpg" mimetype="image"/> <inline-graphic xlink:href="pcbi.1000163.e016.jpg" mimetype="image"/> <inline-graphic xlink:href="pcbi.1000163.e017.jpg" mimetype="image"/> <graphic xlink:href="pcbi.1000163.e018.jpg" mimetype="image" position="float"/> <inline-graphic xlink:href="pcbi.1000163.e019.jpg" mimetype="image"/> <inline-graphic xlink:href="pcbi.1000163.e020.jpg" mimetype="image"/> <inline-graphic xlink:href="pcbi.1000163.e021.jpg" mimetype="image"/> <graphic xlink:href="pcbi.1000163.e022.jpg" mimetype="image" position="float"/>
<media xlink:href="pcbi.1000163.s001.zip"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pcbi.1000163.s002.zip"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pcbi.1000163.s003.eps"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pcbi.1000163.s004.eps"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pcbi.1000163.s005.avi"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pcbi.1000163.s006.avi"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pcbi.1000163.s007.avi"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pcbi.1000163.s008.avi"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pcbi.1000163.s009.avi"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pcbi.1000163.s010.avi"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pcbi.1000163.s011.avi"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pcbi.1000163.s012.avi"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pcbi.1000163.s013.avi"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pcbi.1000163.s014.avi"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pcbi.1000163.s015.avi"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pcbi.1000163.s016.avi"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pcbi.1000163.s017.avi"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pcbi.1000163.s018.avi"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pcbi.1000163.s019.avi"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pcbi.1000163.s020.avi"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pcbi.1000163.s021.avi"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pcbi.1000163.s022.avi"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pcbi.1000163.s023.avi"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pcbi.1000163.s024.avi"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pcbi.1000163.s025.avi"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pcbi.1000163.s026.avi"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pcbi.1000163.s027.avi"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pcbi.1000163.s028.avi"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pcbi.1000163.s029.avi"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pcbi.1000163.s030.avi"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pcbi.1000163.s031.avi"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pcbi.1000163.s032.avi"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pcbi.1000163.s033.avi"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pcbi.1000163.s034.avi"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pcbi.1000163.s035.avi"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pcbi.1000163.s036.avi"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pcbi.1000163.s037.avi"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pcbi.1000163.s038.avi"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pcbi.1000163.s039.avi"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pcbi.1000163.s040.avi"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pcbi.1000163.s041.avi"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pcbi.1000163.s042.avi"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pcbi.1000163.s043.avi"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pcbi.1000163.s044.avi"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pcbi.1000163.s045.avi"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pcbi.1000163.s046.avi"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pcbi.1000163.s047.avi"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pcbi.1000163.s048.avi"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pcbi.1000163.s049.avi"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pcbi.1000163.s050.avi"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pcbi.1000163.s051.avi"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pcbi.1000163.s052.avi"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pcbi.1000163.s053.avi"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pcbi.1000163.s054.avi"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pcbi.1000163.s055.avi"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pcbi.1000163.s056.avi"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pcbi.1000163.s057.avi"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pcbi.1000163.s058.avi"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pcbi.1000163.s059.avi"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pcbi.1000163.s060.avi"><caption><p>Click here for additional data file.</p></caption></media>
[{"label": ["5"], "element-citation": ["\n"], "surname": ["Helmlinger", "Endo", "Ferrara", "Hlatky", "Jain"], "given-names": ["G", "M", "N", "L", "RK"], "year": ["2000"], "article-title": ["Growth factors - formation of endothelial cell networks."], "source": ["Nature"], "volume": ["405"], "fpage": ["139"], "lpage": ["141"]}, {"label": ["21"], "element-citation": ["\n"], "surname": ["Merks", "Newman", "Glazier", "Sloot", "Chopard", "Hoekstra"], "given-names": ["RMH", "SA", "JA", "PMA", "B", "AG"], "year": ["2004"], "article-title": ["Cell-oriented modeling of in vitro capillary development."], "source": ["Cellular Automata: 6th International Conference on Cellular Automata for Research and Industry"], "publisher-loc": ["Berlin"], "publisher-name": ["Springer Verlag, Lect. Notes Comput. Sc. 3305,"], "fpage": ["425"], "lpage": ["434"]}, {"label": ["25"], "element-citation": ["\n"], "surname": ["Ambrosi", "Gamba", "Serini"], "given-names": ["D", "A", "G"], "year": ["2004"], "article-title": ["Cell directional persistence and chemotaxis in vascular morphogenesis."], "source": ["B Math Biol"], "volume": ["66"], "fpage": ["1851"], "lpage": ["1873"]}, {"label": ["26"], "element-citation": ["\n"], "surname": ["Tosin", "Ambrosi", "Preziosi"], "given-names": ["A", "D", "L"], "year": ["2006"], "article-title": ["Mechanics and chemotaxis in the morphogenesis of vascular networks."], "source": ["B Math Biol"], "volume": ["68"], "fpage": ["1819"], "lpage": ["1836"]}, {"label": ["28"], "element-citation": ["\n"], "surname": ["Anderson", "Chaplain"], "given-names": ["ARA", "MAJ"], "year": ["1998"], "article-title": ["Continuous and discrete mathematical models of tumor-induced angiogenesis."], "source": ["B Math Biol"], "volume": ["60"], "fpage": ["857"], "lpage": ["899"]}, {"label": ["30"], "element-citation": ["\n"], "surname": ["Sun", "Wheeler", "Obeyesekere", "Patrick"], "given-names": ["SY", "MF", "M", "CW"], "year": ["2005"], "article-title": ["A deterministic model of growth factor-induced angiogenesis."], "source": ["B Math Biol"], "volume": ["67"], "fpage": ["313"], "lpage": ["337"]}, {"label": ["33"], "element-citation": ["\n"], "surname": ["Manoussaki"], "given-names": ["D"], "year": ["2003"], "article-title": ["A mechanochemical model of angiogenesis and vasculogenesis."], "source": ["ESAIM-Mathematical Modelling and Numerical Analysis"], "volume": ["37"], "fpage": ["581"], "lpage": ["599"]}, {"label": ["35"], "element-citation": ["\n"], "surname": ["Glazier", "Graner"], "given-names": ["JA", "F"], "year": ["1993"], "article-title": ["Simulation of the differential adhesion driven rearrangement of biological cells."], "source": ["Phys Rev E"], "volume": ["47"], "fpage": ["2128"], "lpage": ["2154"]}, {"label": ["36"], "element-citation": ["\n"], "surname": ["Savill", "Hogeweg"], "given-names": ["NJ", "P"], "year": ["1997"], "article-title": ["Modelling morphogenesis: from single cells to crawling slugs."], "source": ["J Theor Biol"], "volume": ["184"], "fpage": ["229"], "lpage": ["235"]}, {"label": ["37"], "element-citation": ["\n"], "surname": ["Merks", "Glazier"], "given-names": ["RMH", "JA"], "year": ["2005"], "article-title": ["A cell-centered approach to developmental biology."], "source": ["Phys A"], "volume": ["352"], "fpage": ["113"], "lpage": ["113"]}, {"label": ["38"], "element-citation": ["\n"], "surname": ["Anderson", "Chaplain", "Rejniak"], "given-names": ["ARA", "MAJ", "RA"], "year": ["2007"], "source": ["Single-cell-based models in biology and medicine"], "publisher-loc": ["Basel, Boston and Berlin"], "publisher-name": ["Birkha\u00fcser Verlag"]}, {"label": ["47"], "element-citation": ["\n"], "surname": ["Glazier", "Balter", "Pop\u0142awksi", "Anderson", "Chaplain", "Rejniak"], "given-names": ["JA", "A", "NJ", "ARA", "MAJ", "RA"], "year": ["2007"], "article-title": ["Magnetization to morphogenesis: A history of the GGH model."], "source": ["Single-cell-based models in biology and medicine"], "publisher-loc": ["Basel Boston Berlin"], "publisher-name": ["Birkha\u00fcser, Mathematics and Biosciences in Interaction"], "fpage": ["79"], "lpage": ["106"]}]
PLoS Comput Biol. 2008 Sep 19; 4(9):e1000163