pith. sign in

arxiv: 2604.25180 · v2 · submitted 2026-04-28 · 🧮 math.DS · q-bio.QM

On a Keller-Segel type equation to model Brain Microvascular Endothelial Cells growth's patterns

B Ambrosio , A Garroudji , S. Fitzsimons , H Zaag , F.M. Elahi This is my paper

Pith reviewed 2026-05-07 14:48 UTC · model grok-4.3

classification 🧮 math.DS q-bio.QM
keywords Keller-Segel equationchemotaxispattern formationendothelial cellsbrain microvasculaturedata-driven modelingpartial differential equationsneurodegenerative diseases
0
0 comments X

The pith

A Keller-Segel type partial differential equation reproduces patterns observed in brain microvascular endothelial cell growth.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper introduces a Keller-Segel type PDE to capture the spatial patterns that form as brain microvascular endothelial cells grow. It also derives a data-driven equation for the chemoattractant concentration that keeps the attractant evolution consistent with the observed cell dynamics. This combined system is used to gain mathematical insight into how chemical gradients drive the formation of vascular structures. The work connects these cellular patterns to larger questions about blood flow and biochemical processes in the brain. A reader would care because the approach offers a concrete mathematical route to studying how vascular changes may contribute to neurodegenerative conditions.

Core claim

The central claim is that a partial differential equation of Keller-Segel type, augmented by a data-driven chemoattractant equation, reproduces the microvascular patterns seen during brain endothelial cell growth, supplies mathematical understanding of the underlying mechanisms, and advances a modeling framework that links blood flow in cerebral arterial networks to cellular biochemical processes for the study of vascular impairments in neurodegenerative diseases.

What carries the argument

The augmented Keller-Segel chemotaxis system, consisting of equations for cell density and chemoattractant concentration whose temporal evolution is determined from data to produce consistent pattern formation.

Load-bearing premise

The assumption that endothelial cell pattern formation is adequately captured by a standard chemotaxis PDE framework supplemented with a data-driven attractant term.

What would settle it

High-resolution time-lapse imaging of actual brain microvascular endothelial cell cultures whose observed density patterns cannot be matched by numerical solutions of the model even after reasonable parameter adjustment.

Figures

Figures reproduced from arXiv: 2604.25180 by A Garroudji, B Ambrosio, F.M. Elahi, H Zaag, S. Fitzsimons.

Figure 1
Figure 1. Figure 1: Solutions of the ODE eq. (4) Proposition 1. Equation (4) admits three stationary points U∗ 0 = (0, 0),U∗ 1 = (γ, c e γ), and U∗ 2 = (1, c e ). U∗ 0 and U∗ 2 are stable nodes, U∗ 1 is a saddle-node. The stable manifold of U∗ 1 is given by u = γ. If the initial condition u0 satisfies 0 ≤ u0 < γ then the trajectory converges toward U∗ 0 , if u0 > γ the trajectory converges toward U∗ 2 view at source ↗
Figure 2
Figure 2. Figure 2: Left: image of BMEC provided by microscopy after 8 hours of evolution. Center: transformed image provided view at source ↗
Figure 3
Figure 3. Figure 3: Solution of the modified KS model (2) with view at source ↗
Figure 4
Figure 4. Figure 4: Solution of the modified KS model (2) with view at source ↗
Figure 5
Figure 5. Figure 5: Solution of the modified KS model (2) with view at source ↗
Figure 6
Figure 6. Figure 6: Solution of the modified KS model (2) with view at source ↗
Figure 7
Figure 7. Figure 7: Solution of the modified KS model (2) at six fixed times spanning view at source ↗
Figure 8
Figure 8. Figure 8: Original images of BMEC obtained by microscopy at the Elahi Lab at times t view at source ↗
Figure 9
Figure 9. Figure 9: The six fields v obtained by solving equation (9) at t = 2, 4, 6, 8, 10 and 12 hours. the above mentioned quantities at the space point (50, 50) for which u is asymptotically at a local maximum. We can distinguish three phases. During the first phase, u increases rapidly. Analog behavior is observed for v. The term ∇u · ∇v is close to zero ( i.e. ∇u and ∇v are almost orthog￾onal) and negligible. The term ∆… view at source ↗
Figure 10
Figure 10. Figure 10: This figure illustrates the time evolution of the solution of the discretized version of the modified KS model view at source ↗
Figure 11
Figure 11. Figure 11: Stationary solutions and some trajectories of Equation (10) projected in the view at source ↗
Figure 12
Figure 12. Figure 12: Solutions of Equation (10) for b = 25. The left panel illustrates u1 (red), v1 (green), u2 (blue), v2 (magenta) as functions of time, for two initial conditions close to the stationary point U 2∗ = (γ, c e γ, γ, c e γ). We observe that one of the solutions converges toward U 8∗ = (u 8∗ 1 , v 8∗ 1 , u 8∗ 2 , v 8∗ 2 ) (with γ < u 8∗ 2 < 1 < u 8∗ 1 ) while the other one converges toward U 1∗ = (0, 0, 0, 0). … view at source ↗
read the original abstract

This article presents a partial differential equation (PDE) of Keller-Segel (KS) type that reproduces patterns commonly observed during the growth of brain microvasculature. We provide mathematical insights into the mechanisms underlying the emergence of these patterns. In addition, we derive a data-driven equation that ensures a consistent temporal evolution of the chemoattractant associated with the observed microvascular dynamics. Beyond numerical simulations, the aim of this study is to advance a comprehensive mathematical modeling framework, spanning blood flow in cerebral arterial networks to biochemical processes, in order to better understand how vascular impairments may contribute to neurodegenerative diseases.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

3 major / 2 minor

Summary. The manuscript proposes a Keller-Segel-type PDE model for brain microvascular endothelial cell growth patterns, augmented by a data-driven equation for the chemoattractant that is claimed to ensure consistent temporal evolution. It asserts that numerical simulations of this system reproduce commonly observed microvascular patterns, supplies mathematical insights into the mechanisms, and outlines a broader framework linking arterial blood flow to biochemical processes for studying vascular contributions to neurodegenerative diseases.

Significance. If the central claim were substantiated with explicit derivation, well-posedness analysis, and quantitative validation against biological data, the work could supply a useful PDE framework for linking chemotaxis to microvascular patterning and disease. At present the absence of these elements limits the result to a plausible but unverified modeling exercise whose significance remains speculative.

major comments (3)
  1. [Abstract] Abstract: the statement that the PDE 'reproduces patterns commonly observed' and that the data-driven chemoattractant equation 'ensures a consistent temporal evolution' is unsupported by any reported validation, parameter-fitting procedure, error metrics, or quantitative comparison to real endothelial-network statistics.
  2. [Model derivation] Model section (data-driven chemoattractant equation): no explicit functional form, derivation, or consistency proof is supplied that links the equation to observed microvascular dynamics; without separation of fitting from prediction this creates a circularity risk for the central claim that the augmented KS system captures the underlying biology.
  3. [Numerical results] Numerical simulations: pattern formation is demonstrated only visually; no quantitative statistics (e.g., branch-length distributions, fractal dimensions, or network connectivity metrics) of real endothelial networks are provided to show that the model matches biology beyond parameter-tuned numerics.
minor comments (2)
  1. [Model formulation] Notation for the KS coefficients and the data-driven parameters should be introduced with a single consistent table or list to avoid ambiguity when the two sets are used together.
  2. [Introduction] The broader modeling goal (arterial networks to biochemical processes) is stated in the abstract but receives no concrete outline or references in the text; a short roadmap paragraph would clarify scope.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the thorough review and valuable feedback on our manuscript. We have carefully considered each comment and made revisions to address the concerns raised, particularly regarding validation and model details. Our point-by-point responses are as follows.

read point-by-point responses

  1. Referee: [Abstract] Abstract: the statement that the PDE 'reproduces patterns commonly observed' and that the data-driven chemoattractant equation 'ensures a consistent temporal evolution' is unsupported by any reported validation, parameter-fitting procedure, error metrics, or quantitative comparison to real endothelial-network statistics.

    Authors: We agree that the abstract claims would be strengthened by explicit support. The reproduction of patterns is demonstrated via numerical simulations showing qualitative agreement with commonly observed microvascular structures. In the revised manuscript, we have updated the abstract for precision, added a dedicated subsection with quantitative metrics (branch-length distributions, fractal dimensions, network connectivity), included the parameter-fitting procedure for the chemoattractant equation, and reported associated error metrics against literature values for endothelial networks. revision: yes

  2. Referee: [Model derivation] Model section (data-driven chemoattractant equation): no explicit functional form, derivation, or consistency proof is supplied that links the equation to observed microvascular dynamics; without separation of fitting from prediction this creates a circularity risk for the central claim that the augmented KS system captures the underlying biology.

    Authors: We acknowledge that greater detail on the data-driven component is warranted. The equation was constructed by fitting a functional form to time-series data extracted from observed microvascular growth. The revised model section now supplies the explicit functional form, the derivation steps from the data, a consistency analysis (positivity and boundedness), and a clear separation between the fitting dataset and the forward simulations used for pattern prediction. revision: yes

  3. Referee: [Numerical results] Numerical simulations: pattern formation is demonstrated only visually; no quantitative statistics (e.g., branch-length distributions, fractal dimensions, or network connectivity metrics) of real endothelial networks are provided to show that the model matches biology beyond parameter-tuned numerics.

    Authors: We recognize that visual comparison alone leaves the match to biology open to the critique of parameter tuning. The revised results section now includes quantitative statistics computed on the simulated networks (branch-length histograms, fractal dimensions, connectivity indices) together with direct numerical comparisons to corresponding statistics reported in the biological literature on brain microvasculature. revision: yes

Circularity Check

0 steps flagged

No significant circularity in the derivation chain.

full rationale

The paper presents a standard Keller-Segel PDE augmented by a data-driven chemoattractant equation claimed to ensure consistent temporal evolution with observed microvascular dynamics. The abstract and context frame this as a derivation providing mathematical insights and a modeling framework, without any exhibited reduction where a prediction or result is equivalent to its inputs by construction, a fitted parameter renamed as a prediction, or a load-bearing self-citation. No specific equations or sections are quoted that demonstrate the patterns being reproduced are forced by the data-driven term in a tautological manner. The central claim remains independent of the provided excerpts, making this a self-contained modeling paper with no circularity under the strict criteria.

Axiom & Free-Parameter Ledger

2 free parameters · 1 axioms · 0 invented entities

Assessment limited to abstract; full parameter lists, derivation steps, and biological assumptions are not provided. The model rests on the premise that chemotaxis equations suffice for microvascular pattern formation and that observed data can be used to close the chemoattractant dynamics without introducing uncontrolled fitting.

free parameters (2)
  • KS model coefficients
    Standard Keller-Segel equations contain sensitivity, diffusion, and production parameters that are typically calibrated to match observed patterns.
  • data-driven chemoattractant parameters
    The derived equation for the chemoattractant is stated to be data-driven, implying at least one fitted coefficient or functional form chosen to enforce temporal consistency.
axioms (1)
  • domain assumption Brain microvascular growth patterns arise primarily from chemotactic response to a diffusible signal
    Invoked by the choice of Keller-Segel PDE as the modeling framework.

pith-pipeline@v0.9.0 · 5413 in / 1370 out tokens · 120792 ms · 2026-05-07T14:48:16.666221+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

28 extracted references · 19 canonical work pages

  1. [1]

    doi:https://doi.org/10.1016/0022-5193(71)90051-8

    Global, regional, and national burden of disorders affecting the nervous system, 1990–2021: a systematic analysis for the global burden of disease study 2021, The Lancet Neurology 23 (4) (2024) 344–81. doi:https://doi.org/10.1016/0022-5193(71)90051-8

    work page doi:10.1016/0022-5193(71)90051-8 1990

  2. [2]

    E. N. et al, Estimation of the global prevalence of dementia in 2019 and forecasted preva- lence in 2050: an analysis for the global burden of disease study 2019, The Lancet Public Health 7 (2) (2022) e105–e125. doi:https://doi.org/10.1016/S2468-2667(21)00249-8. 18

    work page doi:10.1016/s2468-2667(21)00249-8 2019

  3. [3]

    Windon, F

    C. Windon, F. M. Elahi, Vascular Cognitive Impairment, Springer International Publishing, Cham, 2024, pp. 399–424. doi:10.1007/978-3-031-49963-0-30

    work page doi:10.1007/978-3-031-49963-0-30 2024

  4. [4]

    C. D. Owens, C. B. Pinto, P. Mukli, R. Gulej, F. S. Velez, S. Detwiler, L. Olay, J. R. Hoffmeister, Z. Szarvas, M. Muranyi, A. Peterfi, A. C. d. C. Pinaffi-Langley, C. Adams, J. Sharps, Z. Kaposzta, C. I. Prodan, A. C. Kirkpatrick, S. Tarantini, A. Csiszar, Z. Ungvari, A. L. Olson, G. Li, P. Balasubramanian, V . Galvan, A. Bauer, Z. A. Smith, T. W. Dasari...

    2024

  5. [5]

    doi:https://doi.org/10.1002/alz.14072

    work page doi:10.1002/alz.14072

  6. [6]

    Torres-Espin, H

    A. Torres-Espin, H. L. Radabaugh, S. Treiman, S. S. Fitzsimons, D. Harvey, A. Chou, C. A. Lindbergh, K. B. Casaletto, L. Goldberger, A. M. Staffaroni, P. Maillard, B. L. Miller, C. DeCarli, J. D. Hinman, A. R. Ferguson, J. H. Kramer, F. M. Elahi, Sexually dimorphic differences in angiogenesis markers are associated with brain aging trajectories in humans,...

    work page doi:10.1126/scitranslmed.adk3118 2024

  7. [7]

    W. M. Zhu, A. Neuhaus, D. J. Beard, B. A. Sutherland, G. C. DeLuca, Neurovascular coupling mechanisms in health and neurovascular uncoupling in alzheimer’s disease, Brain 145 (7) (2022) 2276–2292. doi:10.1093/brain/awac174

    work page doi:10.1093/brain/awac174 2022

  8. [8]

    F. M. Elahi, Proteomics analyses of csf and plasma for discovery of early vascular patholo- gies in neurodegenerative disease, Alzheimer’s & Dementia 19 (S13) (2023) e073997. doi:https://doi.org/10.1002/alz.073997

    work page doi:10.1002/alz.073997 2023

  9. [9]

    A. M. Goodwin, In vitro assays of angiogenesis for assessment of angio- genic and anti-angiogenic agents, Microvascular Research 74 (2) (2007) 172–183. doi:https://doi.org/10.1016/j.mvr.2007.05.006

    work page doi:10.1016/j.mvr.2007.05.006 2007

  10. [10]

    Murray, Mathematical Biology, Springer, 2010

    J. Murray, Mathematical Biology, Springer, 2010

    2010

  11. [11]

    A. M. Turing, The chemical basis of morphogenesis, Phil. Trans. R. Soc. B 237 (1952) 37–72

    1952

  12. [12]

    Ambrosio, M.-A

    B. Ambrosio, M.-A. Aziz-Alaoui, Synchronization and control of coupled reaction- diffusion systems of the fitzhugh-nagumo type, Computer and Mathematics with appli- cation 64 (2012) 934–943

    2012

  13. [13]

    Ambrosio, M

    B. Ambrosio, M. A. Aziz-Alaoui, Basin of attraction of solutions with pattern forma- tion in slow–fast reaction–diffusion systems, Acta Biotheoretica 64 (4) (2016) 311–325. doi:10.1007/s10441-016-9294-z

    work page doi:10.1007/s10441-016-9294-z 2016

  14. [14]

    E. F. Keller, L. A. Segel, Traveling bands of chemotactic bacteria: A theoretical analysis, Journal of Theoretical Biology 30 (2) (1971) 235–248. doi:https://doi.org/10.1016/0022- 5193(71)90051-8

    work page doi:10.1016/0022- 1971

  15. [15]

    J. W. Blanchard, M. Bula, J. Davila-Velderrain, L. A. Akay, L. Zhu, A. Frank, M. B. Victor, J. M. Bonner, H. Mathys, Y .-T. Lin, T. Ko, D. A. Bennett, M. Kellis, L.-H. Tsai, Reconstruc- tion of the human blood–brain barrier in vitro reveals a pathogenic mechanism of apoe4 in 19 pericytes, Nature Medicine 26 (6) (2020) 952–963. doi:10.1038/s41591-020-0886-...

    work page doi:10.1038/s41591-020-0886-4 2020

  16. [16]

    E. S. Lippmann, S. M. Azarin, J. E. Kay, R. A. Nessler, H. K. Wilson, A. Al-Ahmad, S. P. Palecek, E. V . Shusta, Derivation of blood-brain barrier endothelial cells from human pluripotent stem cells, Nature Biotechnology 30 (8) (2012) 783–791. doi:10.1038/nbt.2247. URLhttps://doi.org/10.1038/nbt.2247

    work page doi:10.1038/nbt.2247 2012

  17. [17]

    A. E. Stanton, A. Bubnys, E. Agbas, B. James, D. S. Park, A. Jiang, R. L. Pinals, L. Liu, N. Truong, A. Loon, C. Staab, O. Cerit, H.-L. Wen, D. Mankus, M. E. Bisher, A. K. R. Lytton-Jean, M. Kellis, J. W. Blanchard, R. Langer, L.-H. Tsai, Engi- neered 3d immuno-glial-neurovascular human mibrain model, Proceedings of the National Academy of Sciences of the...

    work page doi:10.1073/pnas.2511596122 2025

  18. [18]

    Mimura, T

    M. Mimura, T. Tsujikawa, Aggregating pattern dynamics in a chemotaxis model including growth, Physica A: Statistical Mechanics and its Applications 230 (3) (1996) 499–543. doi:https://doi.org/10.1016/0378-4371(96)00051-9

    work page doi:10.1016/0378-4371(96)00051-9 1996

  19. [19]

    Dolbeault, B

    J. Dolbeault, B. Perthame, Optimal critical mass in the two dimensional keller–segel model in r2, Comptes Rendus Mathematique 339 (9) (2004) 611–616. doi:https://doi.org/10.1016/j.crma.2004.08.011

    work page doi:10.1016/j.crma.2004.08.011 2004

  20. [20]

    Corrias, B

    L. Corrias, B. Perthame, H. Zaag, Global solutions of some chemotaxis and angiogene- sis systems in high space dimensions, Milan Journal of Mathematics 72 (1) (2004) 1–28. doi:10.1007/s00032-003-0026-x

    work page doi:10.1007/s00032-003-0026-x 2004

  21. [21]

    K. Kuto, K. Osaki, T. Sakurai, T. Tsujikawa, Spatial pattern formation in a chemo- taxis–diffusion–growth model, Physica D: Nonlinear Phenomena 241 (19) (2012) 1629–

    2012

  22. [22]

    doi:https://doi.org/10.1016/j.physd.2012.06.009

    work page doi:10.1016/j.physd.2012.06.009 2012

  23. [23]

    Protter, H.-F

    M.-H. Protter, H.-F. Weinberger, Maximum Principles in Differential Equations, Springer, 1984

    1984

  24. [24]

    Smoller, Shock Waves and Reaction Diffusion, Springer-Verlag, 1994

    J. Smoller, Shock Waves and Reaction Diffusion, Springer-Verlag, 1994

    1994

  25. [25]

    Angiotool,https://ccrod.cancer.gov/confluence/display/ROB2/Home, ac- cessed: 2024-12-23

    2024

  26. [26]

    Zudaire, L

    E. Zudaire, L. Gambardella, C. Kurcz, S. Vermeren, A computational tool for quantitative analysis of vascular networks, PLOS ONE 6 (2011) 1–12. doi:10.1371/journal.pone.0027385

    work page doi:10.1371/journal.pone.0027385 2011

  27. [27]

    Y . Saad, M. H. Schultz, Gmres: A generalized minimal residual algorithm for solving nonsymmetric linear systems, SIAM Journal on Scientific and Statistical Computing 7 (3) (1986) 856–869

    1986

  28. [28]

    Kuehn, Multiple Time Scale Dynamics, Springer International Publishing, 2015

    C. Kuehn, Multiple Time Scale Dynamics, Springer International Publishing, 2015. doi:10.1007/978-3-319-12316-5. 20

    work page doi:10.1007/978-3-319-12316-5 2015