The renal cortex drives renal function. Hypoxia/reoxygenation are primary factors in ischemia-reperfusion (IR) injuries, but renal oxygenation per se is complex and awaits full elucidation. Few mathematical models address this issue: none captures cortical tissue heterogeneity. Using agent-based modeling, we develop the first model of cortical oxygenation at the cell-tissue level (RCM), based on first principles and careful bibliographical analysis. Entirely parameterized with Rat data, RCM is a morphometrically equivalent 2D-slice of cortical tissue, featuring peritubular capillaries (PTC), tubules and interstitium. It implements hemoglobin/O 2 binding-release, oxygen diffusion, and consumption, as well as capillary and tubular flows. Inputs are renal blood flow RBF and PO 2 feeds; output is average tissue PO 2 (tPO 2 ). After verification and sensitivity analysis, RCM was validated at steady-state (tPO 2 37.7 ± 2.2 vs. 36.9 ± 6 mmHg) and under transients (ischemic oxygen half-time: 4.5 ± 2.5 vs. 2.3 ± 0.5 s in situ). Simulations confirm that PO 2 is largely independent of RBF, except at low values. They suggest that, at least in the proximal tubule, the luminal flow dominantly contributes to oxygen delivery, while the contribution of capillaries increases under partial ischemia. Before addressing IR-induced injuries, upcoming developments include ATP production, adaptation to minutes-hours scale, and segmental and regional specification. Int. J. Mol. Sci. 2019, 20, 6246 2 of 35heterogeneities [19]. This complexity pushed physiologists to resort to mathematical modeling and computer simulation (MS), as a complement to the experimental approaches [10,[20][21][22][23].Renal oxygen distribution started to be addressed under the MS angle in the late 90's [12,20]. The AV oxygen shunt has been long hinted and its mathematical descriptions are based on careful anatomic determinations. However, its physiological relevance and quantification remain debated [23][24][25][26]. Finally, oxygen distribution within the medulla, with its countercurrent system and complex organization, has been modeled with attention and a great level of histological detail [27][28][29][30].Conversely, oxygen distribution in the cortex has been addressed only lately, in the mid-2010's. Likely, the urine concentration and NaCl reabsorption role of the medulla have somewhat eclipsed the cortex [30,31]. Further, since the cortex is highly perfused and exhibits PO 2 levels similar to other organs' (30-60 mmHg), its oxygenation has been considered unlimited. Finally, the cortical apparent disorder (cortical tubular labyrinth, entangled with the peritubular capillary network), is less amenable to geometric treatment than the highly organized medullary regions [32].Recently though, a few modeling studies of oxygen use and distribution in the renal cortex have been proposed. These MS studies follow either one of two conceptual lines: one focuses on nephron segmentation and epithelial polarity and transport, while the other ignores epithelial "details" b...