Additional file 7: Table S1. of Tip cell overtaking occurs as a side effect of sprouting in computational models of angiogenesis

Effect of VEGF gradients on the sprout tip occupancy by WT cells. The mean (out of 10 simulations) occupancy of sprout tips by WT cells at the end of a simulation with the contact inhibition model with differential adhesion between tip and stalk cells, in which only tip cells have a chemotactic sensitivity (λ c,VEGF = 0.1) to an external VEGF field, is listed for different VEGF gradient shapes (columns) and for different ratios of WT and Vegfr2 haploids in the spheroid (rows). The columns represent different shapes of the gradients of the external VEGF field ranging from concentration 0 to 1, which was uniformly spread over the grid, or increased from left to right over the grid in a linear, exponential or sigmoidal fashion. The p-values represent the probability that the total number of simulated sprouts were occupied by at least the indicated percentage of WT cells when assuming only random motion (calculated with a binomial distribution, with n the number of sprouts, k the number of sprouts occupied by WT cells, and p the mixing ratio. Table S2. Effect of VEGF gradient on cell trajectory data. Anterograde coordination, retrograde coordination, and the directional motility is listed for cells in the contact inhibition model (average of ten simulations) with differential adhesion between tip and stalk cells,, in which only tip cells have a chemotactic sensitivity (λ c,VEGF = 0.1) to an external VEGF field. The simulations were initialized with a mix of WT cells and Vegfr2 haploids in a 1:1 ratio. The columns represent different shapes of gradients of the external VEGF field ranging from concentration 0 to 1, which was uniformly spread over the grid, or increased from left to right over the grid in a linear, exponential or sigmoidal fashion. Table S3. Parameter values of the contact inhibition model and the cell elongation model. Table S4. Parameter values VEGF-Dll4-Notch signaling model. Dimensional units: decay rates, γ N , γ D , γ S , γ R , γ A are per 30 s (1 MCS = 30 s), production rates β N , β D , β Dc , β R , β RC in RU/(30 s) and affinities k S , k D , k A in RU ∙ 30 s. Here Relative Units (RU) replace concentrations which are unknown. (PDF 209 kb