This paper presents a mathematical model governing the dynamics of a morphogenetic vascular endothelial cell (EC) during angiogenesis, and vascular growth formed by EC. Especially, we adopt a multiparticle system for modeling these cells. This model does not distinguish a tip cell from a stalk cell. A formed vessel is modeled using phase-field equation to prevent capillary expansion with time stepping in particular. Numerical simulation reveals that all cells are moving in the direction of high concentration of vascular endothelial growth factor (VEGF), and that they are mutually repellent in cases in which they are closer than some threshold.
Citation: Takashi Nakazawa, Sohei Tasaki, Kiyohiko Nakai, Takashi Suzuki. Multicellular model of angiogenesis[J]. AIMS Bioengineering, 2022, 9(1): 44-60. doi: 10.3934/bioeng.2022004
This paper presents a mathematical model governing the dynamics of a morphogenetic vascular endothelial cell (EC) during angiogenesis, and vascular growth formed by EC. Especially, we adopt a multiparticle system for modeling these cells. This model does not distinguish a tip cell from a stalk cell. A formed vessel is modeled using phase-field equation to prevent capillary expansion with time stepping in particular. Numerical simulation reveals that all cells are moving in the direction of high concentration of vascular endothelial growth factor (VEGF), and that they are mutually repellent in cases in which they are closer than some threshold.
[1] | Holderfield MT, Hughes CC (2008) Crosstalk between vascular endothelial growth factor, notch, and transforming growth factor-beta in vascular morphogenesis. Circ Res 102: 637-652. https://doi.org/10.1161/CIRCRESAHA.107.167171 |
[2] | De Smet F, Segura I, De Bock K, et al. (2009) Mechanics of vessels branching: filopodia on endothelial tip cells lead the way. Arterioscle Thromb Vasc Biol 29: 639-649. https://doi.org/10.1161/ATVBAHA.109.185165 |
[3] | Armulik A, Abramsson A, Betsholtz C (2005) Endothelial/pericyte interactions. Circ Res 97: 512-523. https://doi.org/10.1161/01.RES.0000182903.16652.d7 |
[4] | Gaengel K, Genove G, Armulik A, et al. (2009) Endothelial-mural cell signaling in vascular development and angiogenesis. Arterioscle Thromb Vasc Biol 29: 630-638. https://doi.org/10.1161/ATVBAHA.107.161521 |
[5] | Hellstrom M, Gerhardt H, Kalen M, et al. (2001) Lack of pericytes engenders endothelial hyperplasia and abnormal vascular morphogenesis. J Cell Biol 153: 543-554. https://doi.org/10.1083/jcb.153.3.543 |
[6] | Lafleur MA, Forsyth PA, Atkinson SJ, et al. (2001) Perivascular cells regulate endothelial membrane type-1 matrix metalloprotainase activity. Biochem Biophy Res 282: 463-473. https://doi.org/10.1006/bbrc.2001.4596 |
[7] | Liu H, Kennard S, Lilly B (2009) NOTCH3 expression is induced in mural cells through an autoregulatory loop that requires endothelial-expressed JAGGED1. Circ Res 104: 466-475. https://doi.org/10.1161/CIRCRESAHA.108.184846 |
[8] | Frieboes HB, Lowengrub JS, Wise SM, et al. (2007) Computer simulation of glioma and morphology. Neuroimage 37: 59-70. https://doi.org/10.1016/j.neuroimage.2007.03.008 |
[9] | Wise SM, Lowengrub JS, Frieboes HB, et al. (2008) Three-dimensional multispecies nonlinear tumor growth – 1: Model and numerical method. J Theor Biol 253: 524-543. https://doi.org/10.1016/j.jtbi.2008.03.027 |
[10] | Frieboes HB, Jin F, Chuang YL, et al. (2010) Three-dimensional multispecies nonlinear tumor growth – 2: tumor invasion and angiogenesis. Neuroimage 264: 1254-1278. https://doi.org/10.1016/j.jtbi.2010.02.036 |
[11] | Lima EA, Oden JT, Almeida RC (2014) A hybrid ten-species phase-field model of tumor growth. Math Models Methods Appl Sci 24: 2569-2599. https://doi.org/10.1142/S0218202514500304 |
[12] | Perfahl H, Bryne HM, Chen T, et al. (2011) Multiscale modelling of vascular tumor growth in 3D: The role of domain size and boundary conditions. PloS One 6: e14790. https://doi.org/10.1371/journal.pone.0014790 |
[13] | Alber MS, Kiskowski MA, Glazierz JA, et al. (2003) On cellular automaton approaches to modelling biological cells. Mathematical Systems Theory in Biology, Communications, Computation, and Finance. New York: Springer 1-39. https://doi.org/10.1007/978-0-387-21696-6_1 |
[14] | McAneney H, O'Rourke SFC (2007) Investigation of various growth mechanisms of solid tumor growth within the linear-quadratic model for radiotherapy. Phys Med Biol 52: 1039-1054. https://doi.org/10.1088/0031-9155/52/4/012 |
[15] | Rejniak KA, Anderson ARA (2010) Hybrid models of tumor growth. Wires SystBiol Med 3: 115-125. https://doi.org/10.1002/wsbm.102 |
[16] | Boondirek A, Lenbury Y (2006) A stochastic model of cancer growth with immune response. J Kor Phy Phys Soc 4: 1652-1666. |
[17] | Gerlee P, Anderson ARA (2008) A hybrid cellular automaton model of clonal evolution in cancer: The emergence of the glycolytic phenotype. J Thero Biol 250: 705-722. https://doi.org/10.1016/j.jtbi.2007.10.038 |
[18] | Arima S, Nishiyama K, Ko T, et al. (2011) Angeogenesis morphogenesis driven by dynamic and heterogeneous collective endothelial cell movement. Development 138: 4763-4776. https://doi.org/10.1242/dev.068023 |
[19] | Takubo N, Yura F, Naemura K, et al. (2019) Cohesive and anisotropic vascular endotherial cell motility driving angiogenic morphogenesis. Scientific Reports 9: 9034. https://doi.org/10.1038/s41598-019-45666-2 |
[20] | Cahn JW, Hilliard JE (1958) Free energy of a nonuniform system. I. Interfacial free generic interface. J Chem Phys 28: 258-267. https://doi.org/10.1063/1.1744102 |
[21] | Furihata D, Matsuo T (2011) Discrete Variational Derivative Method: A structure-preserving numerical method for partial differential equations. Nwe York: CRC Press. https://doi.org/10.1201/b10387 |
[22] | Feng X, Li Y, Xing Y (2016) Analysis of mixed interior penalty discontinuous Galerkin methods for the Cahn–Hilliard equation and the Hele–Shaw flow. SIAM J Numer Anal 54: 825-847. https://doi.org/10.1137/15M1009962 |
[23] | Kovacs M, Larsson S, Mesforsh A (2011) Finite element approximation of the Cahn–Hilliard–Cook equation. SIAM J Numer Anal 49: 2407-2429. https://doi.org/10.1137/110828150 |
[24] | Brenner SC, Sung Li-y, Zhang H, et al. (2013) A Morley finite element method for the displacement obstacle problem of clamped Kirchhoff plate. J Comp Appl Math 254: 31-42. https://doi.org/10.1016/j.cam.2013.02.028 |
[25] | Li M, Guan X, Mao S (2014) New error estimates of the Morley element for the plate bending problems. J Comp Appl Math 263: 405-416. https://doi.org/10.1016/j.cam.2013.12.024 |
[26] | Hecht F (2012) New development in FreeFem++. J Numerical Math 20: 251-266. https://doi.org/10.1515/jnum-2012-0013 |
[27] | Sánchez-Pérez R, Pavan S, Mazzeo R, et al. (2018) Mutation of a bHLH transcription factor allowed almond domestication. Science 364: 1095-1098. https://doi.org/10.1126/science.aav8197 |