 Research
 Open Access
 Published:
Estimates of the permeability of extracellular pathways through the astrocyte endfoot sheath
Fluids and Barriers of the CNS volume 20, Article number: 20 (2023)
Abstract
Background
Astrocyte endfoot processes are believed to cover all microvessels in the brain cortex and may play a significant role in fluid and substance transport into and out of the brain parenchyma. Detailed fluid mechanical models of diffusive and advective transport in the brain are promising tools to investigate theories of transport.
Methods
We derive theoretical estimates of astrocyte endfoot sheath permeability for advective and diffusive transport and its variation in microvascular networks from mouse brain cortex. The networks are based on recently published experimental data and generated endfoot patterns are based on Voronoi tessellations of the perivascular surface. We estimate corrections for projection errors in previously published data.
Results
We provide structuralfunctional relationships between vessel radius and resistance that can be directly used in flow and transport simulations. We estimate endfoot sheath filtration coefficients in the range \(L_p=2\times 10^{11}\,\hbox {m}\,\hbox {Pa}^{1}\,\hbox {s}^{1}\) to \(3\times 10^{10} \,\hbox {m}\,\hbox {Pa}^{1}\,\hbox {s}^{1}\), diffusion membrane coefficients for small solutes in the range \(C_M= 5 \times 10^{2}\,\hbox {m}^{1}\) to \(6\times 10^{3}\,\hbox {m}^{1}\), and gap area fractions in the range 0.2–0.6%, based on a interendfoot gap width of 20 nm.
Conclusions
The astrocyte endfoot sheath surrounding microvessels forms a secondary barrier to extracellular transport, separating the extracellular space of the parenchyma and the perivascular space outside the endothelial layer. The filtration and membrane diffusion coefficients of the endfoot sheath are estimated to be an order of magnitude lower than those of the extracellular matrix while being two orders of magnitude higher than those of the vessel wall.
Background
Astrocyte endfoot processes have been reported to cover virtually all of the microvasculature in brain gray matter [1,2,3,4,5,6]. The endfoot processes overlap [3] and form a sheath that constitutes the outer boundary of the perivascular space (PVS). Exchange of fluid across the endfoot sheath is vital to maintain homeostasis of the central nervous system [7], and a key component of the proposed glymphatic theory [8]. The extracellular transport pathway, between individual endfoot processes, and the associated permeabilities of the endfoot sheath are relevant for the interpretation of transport phenomena observed for cerebrospinal fluid (CSF), and passively transported substances that are believed to not enter astrocytes in large quantities (such as many MRI contrast agents) and are used in the analysis of flow and transport processes into, out of and within the brain parenchyma. We omit here the discussion of intracellular pathways (see [8,9,10,11] for proposed roles and scientific debate) and perivascular pathways (see e.g. [6]).
The astrocyte endfoot sheath enclosing the microvessels in brain tissue can be viewed as the surface of a tube tiled by individual endfoot processes, cf. [6, Fig.2]. Voronoi tessellations have been successfully used to describe the geometric configuration of cell populations and cell dynamics for decades [12]. Voronoi tessellations appear if cells are grown radially from a center point at constant speed until collision with a neighbor cell growing at the same speed, a simulation process used by [13] to construct virtual astrocyte endfoot processes. However, Voronoi tessellations of a point set can also be more directly constructed as the dual graph of a Delaunay triangulation of the point set. Motivated by the recent work of Wang et al. [6], in which the authors visualized endfoot process gaps in mouse brain resembling Voronoi tessellations, we herein propose their use to generate artificial endfoot patterns. An exemplary realization of such a pattern is shown in Fig. 1 (cf. [6, Fig.2]). To estimate the permeability of the generated cell patterns, the parameterized surface model representation has to be combined with a corresponding crosssectional model representation. A schematic crosssectional cut through a capillary in Fig. 2 introduces the considered perivascular structures and parameterization of the interendfoot gaps.
Few quantitive experimental studies have been published on the geometrical configuration of astrocyte endfeet. Mathiisen et al. [3] conducted an ultrastructural analysis of astrocyte endfeet in capillary vessels in rats, reporting on gap sizes, thickness, and coverage. Individual endfoot processes are separated from neighboring endfeet by gaps of 20 nm on average [3]. Moreover, neighboring endfeet are overlapping [3] and are regularly connected at gap junctions as narrow as 5 nm as described in an early ultrastructural analysis by Brightman and Reese [14]. McCaslin and coworkers [5] report average endfoot densities and average endfoot sheath thickness for capillaries and larger arterial and venous vessels using twophoton microscopy acquired invivo in mice. Recently, Wang et al. [6] analyzed variations in astrocytic endfoot sizes along the vascular tree in mouse brain cortex and hippocampus exvivo using confocal microscopy and their data demonstrates significant differences in endfoot sizes between venous and arterial vessels.
Permeability estimates of the endfoot sheath and its variance in microvascular networks are crucial parameters for computational models of diffusive and advective transport in brain tissue [15, 16]. Previous estimates of the permeability of perivascular and interstitial compartments have been obtained in a number of works [9, 16,17,18]. Asgari et al. [9] estimated the resistance of astrocyte interendfoot gaps based on an idealized geometrical configuration. In [18], this estimate was extended to obtain a brainwide resistance between the periarterial/perivenous compartments and the extracellular space.
There are two shortcomings of the previous analyses. Firstly, the available data has only been partially combined into permeability estimates. As permeability is a crucial material property, we here aim to provide a derivation and resulting estimates. Secondly, the variation of permeability values within cortical microvascular networks has not been estimated. For example, Mathiisen et al. [3] estimated the (interendfoot) gap area fraction based on crosssectional data. Wang et al. [6] reported variations of the average endfoot vessel coverage area with vessel type, and estimated resulting water flux variations, but no resulting gap area fractions. That means the data cannot be directly compared and mapping variation onto a microvascular network requires additional data or model assumptions.
In this work, we will focus on the estimation of the extracellular endfoot sheath permeability and its variability incorporating data on endfoot sheath ultrastructure, endfoot surface area, and variations with vessel diameters. To this end, we estimate parameter distribution in microvascular networks extracted from mouse brain [19]. Based on and parameterized by values from published experimental data [3, 6], we propose a theoretical model based on random tessellations of the endfoot sheath. The model provides estimates for the permeability of the astrocyte endfoot sheath around microvessels to transmembrane transport of fluids and transported substances. By use of this model, we can connect and compare the data obtained by Mathiisen et al. [3] (gap area fraction) and Wang et al. [6] (area variation) and discuss both permeability variations within a microvascular network and networkaveraged quantities.
Methods
Theoretical model of endfoot sheath cell areas and gaps
To generate artificial endfoot sheath coverage patterns, we sequentially sampled uniformly distributed random points^{Footnote 1} on the endfoot sheath surface until we arrived at a prescribed target point density \(\varrho\) (points per area). We did not insert points that are closer than a distance of \(0.3 \sqrt{1/\varrho }\) to any already inserted point avoiding very small endfoot areas. Next, we computed a Voronoi tessellation of the generated point set.^{Footnote 2} The tessellations consist of polygonal faces. Each polygon represents the visible surface.^{Footnote 3} covered by an endfoot process and the polygon edges (also called bisector edges) mark the location of endfootendfoot gaps. The surface is represented by a rectangle of width \(2\pi r_{\textrm{o}}\) (where \(r_{\textrm{o}}\) is the endfoot sheath radius) and height L such that the total area divided by the number of polygons equals the desired mean endfoot area, \({A = 1/\varrho }\).^{Footnote 4} Based on the reported image data by Wang et al. [6], we assumed that the reported vessel diameters include the endfoot sheath.
Assuming a constant gap width \(d_{\textrm{g}}\), we computed \(\phi _{\textrm{g}}\), the area fraction of the surface occupied by intercellular gaps (i.e. the surface available for transmembrane exchange via the extracellular pathway). For this, we multiplied the total edge length \(l_\Sigma\) in the Voronoi tessellation with the intercellular gap width and divide by the total surface area, \(\phi _{\textrm{g}}:= d_{\textrm{g}} l_\Sigma /(L 2\pi r_{\textrm{o}})\). Since the estimated gap area fractions are below 1%, we neglect the influence of considering finitesized gaps on the endfoot area.
For comparison with previously published data, we also computed the average number of gaps counted in crosssectional cuts through the vessel as ’gaps/ring’\(:= l^\perp _\Sigma /L\), where \(l^\perp _\Sigma\) is the total length of the edges after projecting each edge in axial vessel direction. A related number is \(\phi ^\perp _{\textrm{g}}:= d_{\textrm{g}} l^\perp _\Sigma /(L 2\pi r_{\textrm{o}})\) corresponding to the (reduced) gap area fraction obtained when using ’gaps/ring’ as the basis for its estimation (as for example used in [18]).
Combination with experimental data
Wang et al. [6] describe how the area covered by single endfoot process varies along the vascular tree for vessels of different diameters. They report endfoot areas for all analyzed vessels [6, Fig.2] and endfoot areas separate for arterial vessel and venous vessel for vessels with radius \(r_{\textrm{o}} > 7.5\,\upmu \hbox {m}\) [6, Fig.4]. To extract the data shown by Wang et al. [6], we used the opensource image analysis tool WebPlotDigitizer [21]. We extracted the linear regression curve (vessel average) of endfoot area as a function of vessel diameter, and all individual data points and regression trends of the data were classified into arteries and veins. To obtain an estimate over the whole range of vessels, separated into arterial and venous vessels, we constructed functions to fit well the entire range of diameters reported by Wang et al. [6]. The linear regression trends of [6] and our approximation overlaid are shown in Fig. 3.
The areas measured by [6] correspond to plane projections of the endfoot area resulting from the analysis of 2D images rather than 3D reconstructions of the endfoot sheath, cf. [6, Fig.2]. The projection into the image plane underestimates the actual endfoot area by introducing two sources of error: (1) orthogonal projection distorts the vessel surface, and (2) half of the vessel surface is not visible in the projection. Both effects are stronger for smaller vessels where endfeet typically wrap around the vessel. We first quantified these errors based on the generated Voronoi tessellations and virtual projection as described in more detail in Appendix A. We then found a unique mapping between measured and corrected areas, which allows us to correct the projection error, see Appendix A. The second graph in Fig. 3 shows the diameterarea relationship after correction predicted by the model. This corrected diameterarea relationship is used as the basis for all parameter estimates in this work.
Permeability for diffusive transport of passive tracers
As proposed previously, e.g. [9, 18], we conceptually model the endfoot sheath as a porous medium. Since we here only consider the extracellular pathway, a tracer will only diffuse through the interendfoot gaps (pore space) and cannot enter the endfoot processes themselves (solid skeleton). Therefore, the diffusive transport across the endfoot sheath will be diminished by its gap area fraction, \(\phi _{\textrm{g}}\).
The endfoot processes are known to partially overlap [3], cf. Fig. 2 and the diffusive flux over the endfoot sheath is inversely proportional to the gap length (not endfoot sheath thickness), \(l_{\textrm{g}}\), which has been estimated for capillaries as \(l_{\textrm{g}} \approx 0.45\,\upmu \hbox {m}\) by Mathiisen and coworkers [3]. The corresponding endfoot sheath thickness, \(h_{\textrm{ES}}\), is reported to be between 0.02 and \(0.3\,\upmu \hbox {m}\) [3] in capillaries for chemically fixated tissue, while McCaslin et al. [5] report \(h_{\textrm{ES}} \approx 1.0\,\upmu \hbox {m}\) for mouse cortex capillaries invivo and even larger \(h_{\textrm{ES}}\) for arterial and venous vessels. For geometrical reasons, \(l_{\textrm{g}} {\ge } h_{\textrm{ES}}\) (cf. Figure 2). In the absence of quantitative invivo data, we used \(l_{\textrm{g}} = 1.5 h_{\textrm{ES}}\) in the following, based on the measurements of [3] in capillaries using \(h_{\textrm{ES}} = 0.3\,\upmu \hbox {m}\) and \(l_{\textrm{g}} = 0.45\,\upmu \hbox {m}\) . For different values of \(h_{\textrm{ES}}\), this corresponds to simple linear extrapolation. (In the section “Microvascular networks” below, we estimate \(h_{\textrm{ES}}\) based on the vessel radius.)
Due to obstructions in the endfoot gap channel (larger proteins, fibers, or gap junctions [14]) the effective diffusivity may be reduced by a factor \(\alpha\). Nicholson and Hrabětová [22] demonstrate that in the extracellular space of the parenchyma where the intercellular space width is between 20 and 100 nm, the measured ratio of effective to free diffusivity is usually smaller than can be explained by the tortuosity of the pore space alone. For molecules with a hydrodynamic diameter that is onetenth or more of the gap width, \(\alpha\) needs to model sizedependent steric exclusion and restricted diffusion effects [23], for instance, with the Renkin equation [24, 25]. In particular, \(\alpha\) and therefore diffusive permeability is zero for molecules much larger than the gap width. For small molecules, we approximate \(\alpha\) with 1.
The diffusive flux \(F_D\) \([\text {M}\text {T}^{1}\text {L}^{2}]\) through the surface of a vessel segment can be computed as
where \(\Delta c\) is the concentration drop across the endfoot sheath, the gap area fraction \(\phi _{\textrm{g}}\) is given by (5), D is the binary diffusion coefficient in aqueous solution, and
is the diffusion membrane coefficient (in units of \(\hbox {m}^{1}\)). Given a surface area S (for example for a tubular segment of length \(l_{\textrm{v}}\) and radius \(r_{\textrm{v}}\), \(S = 2\pi r_{\textrm{v}} l_{\textrm{v}}\)) with constant \(C_M\), we can compute the amount of a tracer crossing the endfoot sheath per unit time as \(F_D S\). The product \(S C_M D\) is sometimes called permeabilitysurface product or diffusion capacity [26, Ch. 10], in particular when referring to the surface integral of \(C_M D\) in a larger tissue portion. In this work, we normalize the diffusion capacity by the (free) binary diffusion coefficient D.
Permeability for fluid flow
Interendfoot gaps can be imagined as narrow sheetlike channels [3, 14]. Based on the electron microscopy study by Mathiisen and coworkers [3], we assume an interendfoot gap width of \(d_{\textrm{g}} = 20\,\hbox {nm}\). Viscous flow continuum theory can be applied when describing liquid water flow through these channels [27]. We remark that the complex interface region (of \(\approx 1\,\hbox {nm}\) thickness) between the endfoot’s lipid bilayer membrane and the bulk fluid adds some uncertainty to the effective gap width in addition to the uncertainty of width measurements and spatial variations. Therefore, we argue that it is sufficient to approximate the hydraulic transmissibility by a simple parallel plate flow model as \(t_{\textrm{g}} = d_{\textrm{g}}^3 (12 \mu l_{\textrm{g}})^{1}\), cf. [18], where \(\mu\) denotes the dynamic fluid viscosity. This model assumes noslip conditions on the endfoot cell membrane surface. With the transmissibility of the parallel plate model, the flow rate Q \([\text {L}^3\text {T}^{1}]\) through the surface of a tubular segment of length \(l_{\textrm{v}}\) and radius \(r_{\textrm{v}}\) (and thus lateral surface area \(S = 2\pi r_{\textrm{v}} l_{\textrm{v}}\)) can be computed as
where \(\Delta p\) is the effective pressure drop across the endfoot sheath (between extracellular space (ECS) and PVS) and the gap area fraction \(\phi _{\textrm{g}}\) is given by (5) in terms of the vessel radius, \(r_{\textrm{v}}\), and
is the filtration coefficient (in units of \(\hbox {m}\,\hbox {Pa}^{1}\,\hbox {s}^{1}\)). The product \(S L_p\) is also called filtration capacity [26, Ch. 11], in particular when referring to the surface integral of \(L_p\) in a larger tissue portion.
Microvascular networks
We analyzed two microvascular networks (MVN) extracted from the mouse brain cortex in [19] as published in [28]. The raw vessel polylines extracted from segmented voxel images by [19] is smoothed by a DouglasPeucker algorithm [29] using the local vessel radius as tolerance. For vessel classification (arteries and veins), blood pressure values (p) in every vessel segment are simulated with a finite volume method as described in [15] (but neglecting filtration across the bloodbrainbarrier). The boundary conditions are based on estimations computed by Schmid et al. [28, 30]. We solve a modified Poiseuilletype flow using the invivo apparent viscosity relation proposed in [31] scaled to mouse red blood cells (using an average volume of 55 fL). The opensource software DuMu\(^{\textrm{x}}\) [32] was used as a finite volume solver with dunefoamgrid [33] for the network representation.
Using the computed pressure maps, vessel segments were classified as arterial vessels if its pressure exceeds the average pressure of all segments with \(r_{\textrm{v}} \le 4.5\,\upmu \hbox {m}\) (vessel radius excluding endfoot sheath), and as venous vessels otherwise. The networks and the obtained pressure distribution are shown in Appendix B (Fig. 9).
Since the network data is associated with vessel lumen radius data, \(r_{\textrm{v}}\), excluding the endfoot sheath and other perivascular structures but the tiling model is formulated in terms of the total outer radius (\(r_{\textrm{o}}\)) of the astrocyte endfoot sheath, we require a model of how these radii are related. Based on data reported in [5], we assumed a thickness of the endfoot sheath, \(h_{\textrm{ES}}\), of 1 to \(2.5\,\upmu \hbox {m}\). Additionally, we chose the relation \(h_{\textrm{ES}} = 1 + 0.15 (r_{\textrm{v}}  3)\) modeling a linear increase with increasing vessel lumen radius. Moreover, larger vessels with \(r_{\textrm{v}} \le 3\,\upmu \hbox {m}\) are assumed to be sheathed by smooth muscle cells or ensheathing pericytes [34] located in between the endothelial layer and the endfoot sheath. Based on [34, Fig.3], we estimated the smooth muscle cell layer thickness to be approximately equal to \(h_{\textrm{ES}}\). This means ca. \(1\,\upmu \hbox {m}\) for a precapillary arteriole with \(r_{\textrm{v}} = 3\,\upmu \hbox {m}\) and ca. \(2\,\upmu \hbox {m}\) for a penetrating vessel with \(r_{\textrm{v}} = 10\,\upmu \hbox {m}\). Finally, we added the thickness of the endothelial cell layer and basement membrane with \(0.4\,\upmu \hbox {m}\) [35] for all vessels. In summary, \(r_{\textrm{o}} = 2h_{\textrm{ES}}(r_{\textrm{v}}) + 0.4 = 1.3 r_{\textrm{v}} +1.5\) for \(r_{\textrm{v}} \ge 3\,\upmu \hbox {m}\) and \(r_{\textrm{o}} =h_{\textrm{ES}}(r_{\textrm{v}}) + 0.4 = 1.15 r_{\textrm{v}} + 0.95\) otherwise. For the network analysis, the networks are split into 6 vertically stacked analysis layers (layer 0 being closest to the pial surface and layer 5 being closest to the white matter) of \(200\,\upmu \hbox {m}\) thickness (\(100\,\upmu \hbox {m}\) for layer 5). Average values (\(r_{\textrm{v}}\), \(r_{\textrm{o}}\), \(L_p\), \(C_M\)) have been computed as surfaceareaweighted arithmetic averages of all vessels contained in the respective analysis layer.
Results
Astrocyte endfoot area distribution
Endfoot area distribution and resulting gap area fraction predicted by the model for 200 realizations with \(r_{\textrm{o}} =2.9\,\upmu \hbox {m}^{2}\) (capillary) and \(r_{\textrm{o}} =15.0\,\upmu \hbox {m}^{2}\) (venule and arteriole), with (corrected) mean endfoot area A shown in Fig. 3, are reported in Fig. 4. The resulting endfoot area distribution is wellmodeled by a Gamma distribution^{Footnote 5} (with mean \(50\,\upmu \hbox {m}^{2}\) (capillary), \(110\,\upmu \hbox {m}^{2}\) (venule) and \(490\,\upmu \hbox {m}^{2}\) (arteriole), respectively). The resulting gap area fraction distribution is wellapproximated by a normal distribution with mean gap area fractions of 0.0056 (capillary), 0.0038 (venule), and 0.0018 (arteriole). For comparison with [3, 18], we also report the resulting number of gaps counted per vessel crosssection (’gaps/ring’) on average, which is lowest in the capillary (3.2), highest in venules (11.5), and intermediate in arterioles (5.4).
Additionally, the model was run for 50 different diameters with modelcorrected endfoot area A (shown in Fig. 3). For each diameter, we generated 20 samples (a total of \(n=1000\) samples). The resulting data including mean, and 5th, 25th, 75th, and 95th percentile are reported in Fig. 5 for both venous and arterial vessels.
For comparison with previously published data, we computed based on the diameterarea relations that on average, small vessels (\(r_v <4.5\,\upmu \hbox {m}\), average not weighted by radius prevalence in a network) have endfoot density of ca. \(1.9 \times 10^{4}\) endfeet per \(\hbox {mm}^{2}\) surface area. Larger venous vessels (\(r_v >4.5\,\upmu \hbox {m}\)) show ca. \(1.0\times 10^{4}\) endfeet/\(\hbox {mm}^{2}\) and larger arterial vessels (\(r_v >4.5\,\upmu \hbox {m}\)) show the lowest density of ca. \(0.4\times 10^{4}\) endfeet/\(\hbox {mm}^{2}\).
Gap area fraction for different vessels
Using the same \(n=1000\) samples as for the data in Fig. 5, in combination with the gap width and length reported by Mathiisen et al. [3], we computed the resulting gap area fraction \(\phi _{\textrm{g}}\) for each realization. The results are shown in Fig. 6. For small diameters (capillaries), the gap area fraction for venous vessels, \(\phi _{\textrm{g}}^V\), and arterial vessels, \(\phi _{\textrm{g}}^A\), are similar, while for increasing vessel diameters \(\phi _{\textrm{g}}^V < \phi _{\textrm{g}}^A\).
In summary, the model predicts mean \(\phi _{\textrm{g}}^A\) and \(\phi _{\textrm{g}}^V\) given mean endfoot area, A, vessel radius including endfoot sheath, \(r_{\textrm{o}}\), and gap with, \(d_{\textrm{g}}\). Using a gap width of \(d_{\textrm{g}} = 20\,\hbox {nm}\) [3, 14], we obtain the empirical relations,
where the radius \(r_{\textrm{o}}\) is inserted in \(\upmu \hbox {m}\), cf. Fig. 6.
Based on a constant gap length, \(l_g\), the gap area fraction for the smallest capillary vessels is about 3 times higher than that of \(15\,\upmu \hbox {m}\) radius arterial vessel, and about 2 times higher than that of venous vessel of the same caliber. The smaller increase in endfoot size reported by [6] for venous vessels in comparison with arterial vessels leads to effectively higher gap area fractions in venous vessels with increasing vessel radius.
A major source of uncertainty are the values of \(d_{\textrm{g}}\) and \(l_{\textrm{g}}\) for which there are currently quantitive exvivo [3, 14] but no invivo measurement data available, see Discussion.
Permeability for diffusive transport of passive tracers
The resulting estimates for \(C_M\) using an obstruction of \(\alpha =1\) (i.e. no obstructions; smaller \(\alpha\) values would decrease the \(C_M\) estimates) and constant \(l_{\textrm{g}} = 1\,\upmu \hbox {m}\) are shown in Fig. 6. Also taking into consideration the variation of \(h_{\textrm{ES}}\) in the mouse microvascular networks (see Methods), estimates range between \(C_M \approx 500\,\hbox {m}^{1}\) for the largest arterioles and \(C_M \approx 6000\,\hbox {m}^{1}\) for the smallest capillaries. The estimated values are based on an interendfoot gap width of 20 nm. They depend linearly on \(\alpha\), \(d_g\), and \(l_g^{1}\) (see Methods) and must be reviewed by using (2) when better measurement data becomes available or when the results are to be transferred to different species, see “Discussion”.
Permeability for fluid flow
The resulting estimates for \(L_p\) using a viscosity of \(\mu = 0.7 \times10 ^ {3}\,\hbox {Pa}\,\hbox {s}\) (water at \(37^{\circ }\hbox {C}\), larger assumed viscosity values would increase the \(L_p\) estimates) and constant \(l_{\textrm{g}} = 1\,\upmu \hbox {m}\) are given in Fig. 6. Also taking into consideration the variation of \(h_{\textrm{ES}}\) in the microvascular networks (see Methods), estimates range between \(L_p \approx 2\times 10^{11}\,\hbox {m}\,\hbox {Pa}^{1}\,\hbox {s}^{1}\) for the largest arterioles and \(L_p \approx 3\times 10^{10}\,\hbox {m} \,\hbox {Pa}^{1}\,\hbox {s}^{1}\) for the smallest capillaries, cf. Fig 7. The estimated values are based on an interendfoot gap width of 20 nm. They are proportional to \(d_g^3 l_g^{1}\) (see “Methods”) and must be reviewed by using (4) when better measurement data becomes available or when the results are to be transferred to different species. Due to the cubic dependence, the value of \(d_{\textrm{g}}\) is a major uncertainty, see Discussion.
Microvascular networks
We analyzed two microvascular networks from the mouse brain cortex [19], labeled MVN1 and MVN2. Volumes and surfaces computed by assuming cylindrical segments with estimated outer radius \(r_{\textrm{o}}\) (including endfoot sheath) and vessel lumen radius \(r_{\textrm{v}}\) (from [28]) are given in Table 1. MVN2 has a 18% larger vessel volume fraction (2.8% and 3.4% including endfoot sheath) and an 10% larger surfacetovolume ratio than MVN1 (\(1.7\times 10^{4}\,\hbox {m}^{2}\,\hbox {m}^{3}\) and \(1.9\times 10^{4}\,\hbox {m}^{2}\,\hbox {m}^{3}\)). In both networks, the surface area of small vessels (\(r_v < 3\,\upmu \hbox {m}\)) exceeds the area of the larger vessels by a factor 10 or more.
For each vertical depth analysis layer (see “Methods”), the filtration and diffusion capacity as well as averaged filtration and diffusion membrane coefficients are shown in Fig 7. While the filtration coefficient \(L_p\) (and similarly the diffusion membrane coefficient, \(C_M\)) in individual segments differs by a factor 10 between the largest arteriole segments and the smallest capillaries, the layeraveraged coefficient only varies by a factor 2. The surfaceweighted averages over the entire network are found to be \(L_p = 2.1\times 10^{10}\,\hbox {m}\,\hbox {Pa}^{1}\,\hbox {s}^{1}\) and \(C_M = 4.4 \times 10^{3}\,\hbox {m}^{1}\), with only 1% difference between MVN1 and MVN2. The total filtration capacity per tissue volume is found to be \(3.4\times 10^{6}\,\hbox {m}\,\hbox {Pa}^{1}\,\hbox {s}^{1}\,\hbox {m}^{1}\) (MVN1), \(4.1\times 10^{6}\,\hbox {m}\,\hbox {Pa}^{1} \,\hbox {s}^{1}\,\hbox {m}^{1}\) (MVN2). The total diffusion capacity per tissue volume is found to be \(7.2 \times 10^{7}\,\hbox {m}^{2}\) (MVN1), \(8.5 \times 10^{7}\,\hbox {m}^{2}\) (MVN2). The filtration and diffusion capacity peaks in layer 2, where it is between (10–14% larger than in layer 1, 3, and 4 (0 and 5 are excluded from this comparison due to differences in layer size and occupied volume). This coincides with the lowest average vessel radius, cf. Fig 7.
Discussion
Endfoot gap area fraction
In [3], the average number of endfoot gaps per capillary crosssection is reported as 2.5 (2.3 to 2.9 in 3 different animals). For the modeled capillary with \(r_{\textrm{o}} =2.9\,\upmu \hbox {m}\), the predicted number (’gaps/ring’ in Fig. 4) of 3.2 is only slightly larger (20%). However, without correction of the projection error (Appendix A), capillary endfoot size is estimated in [6] at only \(25\,\upmu \hbox {m}^{2}\) corresponding to about 4.0 ’gaps/ring’ (simulated with our model). Hence, the correction by the model allows us to partially resolve an apparent mismatch between the data reported by [3] and [6]. The comparison may be further affected by the different measurement methods employed by [3] and [6], measurement errors, and the quality of the area correction computed by our model. Finally, there might be interspecies variations between rats and mice.
By extrapolating ’gaps/ring’ and the gap width of \(d=20\,\hbox {nm}\) to all of the surface, the authors of [3] conclude that about 0.3% of the endfoot sheath surface is comprised of gaps—a number also used by [18] to estimate endfoot sheath permeability. We note that this computation effectively assumes that gaps run parallel to the longitudinal vessel axis. Under this assumption, we compute for capillaries, a reduced gap area fraction \(\phi ^\perp _{\textrm{g}} = 0.35,\%\) matching well the value obtained in [3]. However, the corresponding actual gap area fraction estimated by our surface tessellation model is \(\phi _{\textrm{g}} = 0.56,\%\), cf. Fig. 4, and therefore significantly larger. The latter value can be considered an improved estimate correcting for yet another projection error.
For a simple estimate of gap area fraction, we may assume a regular hexagonal tiling for \(A=50\,\upmu \hbox {m}^{2}\) and \(d=20\,\hbox {nm}\) corresponds to a value \(\phi _{\textrm{g}} = 0.0052\)^{Footnote 6} (to be compared with 0.0056) for the capillary and \(\phi _{\textrm{g}} = 0.0035\) (compared with 0.0038) for the vein of the same caliber analyzed in Fig. 4. However, regular tiling falls short of providing a model for individual endfoot size variability.
With respect to the variation with vessel type, we remark that assuming constant gap width for all vessels results in a linear correlation between filtration and diffusion membrane coefficient (both quantities depend linearly on the gap area fraction). Therefore, differences in permeability result from variations in the gap area fraction rather than individual gap anatomy. Such a correlation is, for example, also observed for the endothelium of different capillary types [26, Ch. 10.6].
Interendfoot gap width
There is significant uncertainty regarding both gap width invivo and general astrocyte endfoot coverage. The estimates in this work consider a continuous coverage with astrocyte endfeet of all microvessels [5]. Firstly, the actual coverage may be reduced with, for example, astrocyte bodies or microglia substituting endfoot processes on the vessel surface. If the intercellular gap size is not significantly altered, the provided estimates by our model still hold. Secondly, both Mathiisen et al. [3] and Wang et al. [6] worked with chemically fixated tissue. Korogod et al. [39] compared cryogenic and chemical fixation techniques, and report significant differences in the resulting endfeet cavity fraction (37% vs. 4%). At cavity fractions this large, the astrocyte endfoot sheath would be irrelevant in terms of a proposed barrier function. This result is contrasted by invivo observation of continuous coverage [2, 5, 40]. Additionally, Kubotera et al. [41] observed that after laser ablation astrocytes restore the endfoot coverage of microvessels invivo. (Mills and coworkers call this tendency to recover blood vessels after disruptions endfoot plasticity [42].) On the other end, the effective gap area fraction is reduced, if interendfoot gap junctions (2–3 nm [14]) are found to be present invivo with significant density (neglected in this work). It is also reduced for molecules whose hydrodynamic radius is a significant fraction of the gap width (modeled by the parameter \(\alpha\)).
Apart from structural uncertainty, astrocytes are known to change their volume under varying conditions [43]. Changes in cell sizes and changes in the radius of the endfoot sheath could alter its hydrodynamic properties—a potential regulatory mechanism of fluid flow and substance transport [6, 44]. Moreover, vessel diameters are highly dynamic and can dilate up to 30–40% of the vessel diameter [45] which leads to mechanical deformation of the astrocyte endfoot sheath observed invivo [46]. Xie and coworkers [47] found an increase in the parenchymal extracellular volume fraction from ca. 0.14 when mice were awake to ca. 0.23 when they were sleeping. To illustrate the dependence of the estimated permeabilities on the gap width \(d_{\textrm{g}}\), we invite the reader to imagine changes of similar magnitude in the interendfoot gap width \(d_g\). Since \(L_p \propto d_{\textrm{g}}^3\), a 60% increase of \(d_{\textrm{g}}\) would lead to an increase of the \(L_p\) by a factor 4 (a decrease of \(d_{\textrm{g}}\) by 40% would lead to a decrease of \(L_p\) by a factor 5). Since \(C_M \propto d_{\textrm{g}}\), a 60% of \(d_{\textrm{g}}\) would lead to an increase of \(C_M\) by a factor 1.6 (a decrease of \(d_{\textrm{g}}\) by 40% would lead to a decrease of \(C_M\) by a factor 1.7). However, we are not aware of quantitative data describing how the endfoot gap width \(d_{\textrm{g}}\) or endfoot sheath thickness \(h_{\textrm{ES}}\) changes with such alterations. If such data became available, (2) and (4) allow estimating the effect of alterations on the gap area fraction (and \(C_M\), \(L_p\)).
Summarizing to the best of our knowledge, direct evidence for full coverage (or its absence), a precise interendfoot gap width quantification invivo and its variation, as well as quantitative data on temporal dynamics are still lacking. The fluid flow rate \(Q \propto d_{\textrm{g}}^3\) and the diffusive flux \(F_D \propto d_{\textrm{g}}\) strongly depend on the assumed gap width \(d_{\textrm{g}}\) resulting in a large (maybe the largest) source of uncertainty for the estimation of extracellular transport across the endfoot sheath.
Permeability for diffusive transport of passive tracers
The estimated permeability of the endfoot sheath can be compared to adjacent barriers and tissues. The vessel wall is assumed to be virtually impermeable to many molecules. In diseased tissue, for example, neurodegenerative diseases such as multiple sclerosis (MS) or glioma higher permeability has been observed in lesion tissue. For example, MRI contrast agents such as gadobutrol (\(D \approx 3.5 \times 10^{10}\hbox {m}^{2}\,\hbox {s}^{1}\) [15]) can leak out of blood vessels in MS lesions or glioma tissue. In [15], \(C_M D \approx 1 \times 10^{7}\,\hbox {ms}^{1}\) has been estimated for Gadobutrol leakage across the vessel wall in MS lesions, corresponding to \(C_M \approx 300\,\hbox {m}^{1}\) which is an order of magnitude smaller than average values obtained in this work for the endfoot sheath, cf. Fig 7. For skeletal muscle microvascular walls and small hydrophilic molecules, \(C_M\) has been estimated at \(C_M \approx 100\) to \(200\,\hbox {m}^{1}\) [23, Fig.2]. The brain cortex microvascular vessel walls are commonly assumed to be orders of magnitude less permeable than in skeletal muscle. This signifies that the endothelial layer is a much less permeable barrier than the astrocyte endfoot sheath, where we estimated the lowest \(C_M\) for large penetrating arterioles with \(C_M \le 1000\,\hbox {m}^{1}\) and values up to \(C_M \le 6000\,\hbox {m}^{1}\) for capillaries.
The estimated permeability can also be compared to that of the extracellular space (ECS). To this end, we consider a \(1\,\upmu \hbox {m}\) thick slab of ECS. With a porosity of 0.2 and tortuosity factor of 0.35 [48\(1/\lambda ^2\)], we obtain \(C_M = 7\times 10^{4}\,\hbox {m}^{1}\). Hence, the endfoot sheath is more than an order of magnitude less permeable than the ECS given a slab of comparable thickness. Therefore, the endfoot sheath could locally act as a barrier. It could also promote the compartmentalization of substances, depending on whether lowpermeability perivascular pathways parallel to the vessel exist in the vicinity.
We remark that a lower permeability does not necessarily mean slower transport. Moreover, it does not necessarily mean that the layer is limiting along a given pathway. The magnitude of diffusive transport depends on the concentration drop \(\Delta c\), cf. (1) as well as on the distance traveled within a given medium along a specified pathway. As vessel structures constitute thin tubular sources (in an infiltration scenario) or sinks (in a clearance scenario), the magnitude of the concentration gradient can be much larger in the vicinity of the vessels and quickly decays with distance. When, for instance, considering a pathway from the PVS into the ECS, solutes travel a much larger distance in ECS (at least an order of magnitude larger) than across the endfoot sheath. Along this specific pathway, the effective permeability may be approximately equal for the ECS and the extracellular pathway through endfoot sheath. With these considerations, the question as to whether the effect of a comparatively low permeability of the endfoot sheath is significant for any given scenario goes beyond the scope of the present work.
Permeability for fluid flow
In the brain cortex microvasculature, the filtration coefficient of the vessel wall is thought to be very low. Kimura and coworker [49, Tab.3] measured \(L_p=2.8\) to \(4.1\times 10^{12}\,\hbox {m}\,\hbox {Pa}^{1}\,\hbox {s}^{1}\) in single rat brain arterioles. Fraser and Dallas [50] report \(L_p = 2\times 10^{13}\,\hbox {m}\,\hbox {Pa}^{1}\,\hbox {s}^{1}\) in frog brain microvessels. A \(1\,\upmu \hbox {m}\) slab of ECS corresponds to a \(L_p\) value of approximately 0.5 to \(3\times 10^{8}\,\hbox {m}\,\hbox {Pa}^{1}\,\hbox {s}^{1}\) [51,52,53]^{Footnote 7} or larger.^{Footnote 8} We estimated the lowest \(L_p\) values for large arterioles with \(L_p \approx 2\times 10^{11}\,\hbox {m}\,\hbox {Pa}^{1}\,\hbox {s}^{1}\), and the largest values for capillaries, \(L_p \approx 2\times 10^{10}\,\hbox {m}\,\hbox {Pa}^{1}\,\hbox {s}^{1}\).
Hence, similar to the results for diffusion, the endfoot sheath filtration coefficient is one order of magnitude larger than that of the vessel wall. On the other hand, it is two orders of magnitude smaller than a slab of ECS of similar thickness, making the astrocyte endfoot sheath a limiting component for the extracellular PVSECS exchange of fluids.
Using the same parallel plate model as for \(L_p\), (4) and assuming small molecules (with a hydrodynamic molecular diameter much smaller than \(d_{\textrm{g}}\)^{Footnote 9}), we can estimate Péclet numbers for transport through the gaps as \(\text {Pe}_{\textrm{g}} = \frac{\Delta p d_{\textrm{g}}^2}{12\mu \alpha D}\) (using \(l_{\textrm{g}}\) as characteristic length). The Péclet number estimates the importance of advective transport over diffusive transport across the end foot sheath. Given the difference in permeability to that of the vessel wall, across which pressure drops on the order of 1000 Pa may occur due to filtration across the capillary endothelium (estimated for example with the models presented in [15, 56]), we assume maximum pressure drop (\(\Delta p\)) across the endfoot sheath of 100 Pa.^{Footnote 10} Since typical binary diffusion coefficients for small molecules in aqueous solution are in the order of \(1\times 10^{9}\hbox {m}^{2}\,\hbox {s}^{1}\), transport of small solutes across the endfoot sheath is dominated by diffusion (\(\text {Pe}_{\textrm{g}} \lessapprox 0.003 \ll 1\)).
In [9], Asgari and coworkers estimate an interendfootgap permeability of \(L_p = 1.35\times 10^{10}\,\hbox {m}\,\hbox {Pa}^{1}\,\hbox {s}^{1}\) (based on a capillary circumference \(15.7\,\upmu \hbox {m}\), endfoot thickness \(h=1\,\upmu \hbox {m}\) and a parallel plate model, and conversion based an endfoot area of \(78\,\upmu \hbox {m}^{2}\)). This matches well with the average value estimated in arterioles in this work but is lower by about half what we estimate for capillaries. The difference can be explained by the assumption of Asgari that the assumed representative endfoot fully wraps around the vessel (‘gaps/ring‘ is 1) which results in a lower gap area fraction.
In [18], Vinje and colleagues estimate comparable parameters for human brain tissue. In particular, the authors estimated the endfoot sheath resistance (inverse permeability) around arterial and venous vessels (excluding capillaries). The reported resistances correspond to \(L_p = 2\times 10^{10}\,\hbox {m}\,\hbox {Pa}^{1}\,\hbox {s}^{1}\) for arterial and \(L_p = 3\times 10^{10}\,\hbox {m}\,\hbox {Pa}^{1}\,\hbox {s}^{1}\) for venous vessels. The numbers are, in part, based on the gap area fraction estimate provided in [3] for capillaries in rats. As explained above, this number (based on the quantity ’gaps/ring’) results in an underestimation of \(\phi _{\textrm{g}}\) and therefore \(L_p\) of ca. 35% and the suggested higher values would be \(L_p \approx 3\times 10^{10}\,\hbox {m}\,\hbox {Pa}^{1} \,\hbox {s}^{1}\) for arterial and \(L_p \approx 4.5\times 10^{10}\,\hbox {m}\,\hbox {Pa}^{1}\,\hbox {s}^{1}\) for venous vessels. However, we used the endfoot thickness distribution estimated by [5] based on invivo mouse brain data, whereas a constant size straight channel model with \(l_g = h_{\textrm{ES}} = 1\,\upmu \hbox {m}\) is used by [18]. Thus, our resulting permeability for arterioles and venules are approximately half the values of [18], respectively, cf. Fig 7.
Microvascular networks
For the two considered microvascular networks, we find that the filtration and the diffusion capacity are largest at about 40% of cortical depth. This layer also shows the smallest average vessel diameters, cf. Fig 7, and a significant peak in neuron density [57]. A high endfoot density per surface area as in the capillaries, cf. [5], corresponds to a higher permeability of the endfoot sheath due to an increase in the gap area fraction. The average filtration and membrane diffusion coefficients are dominated by the average values for capillaries and appear to be independent of depth. Hence, the increased filtration capacity at 40% seems to be a result of an increased surface area rather than an increased endfoot sheath permeability. This matches with the observation that vessel density is largest in this cortical layer [19, 45, 58].
To the best of our knowledge, the variability of endfoot sizes in the endfoot sheath has not been analyzed using microvascular networks comprising all vessels in a given tissue portion before. Based on the distribution of penetrating arterioles and venous from the macaque cortex [59], Vinje et al. [18] estimate the surface permeability product of the human brain (using an approximate human brain volume of \(V= 1\,\hbox {L}\)). If normalized by the sample volume to eliminate the effect of spatial scale, their estimate corresponds to the volumespecific quantities \(2.2\times 10^{7}\,\hbox {Pa}^{1}\,\hbox {s}^{1}\) for arterioles and \(2.0\times 10^{7}\,\hbox {Pa}^{1}\,\hbox {s}^{1}\)^{Footnote 11} for venules, while we obtain \(6.6\times 10^{8}\,\hbox {Pa}^{1}\,\hbox {s}^{1}\) (MVN1), \(6.6\times 10^{8}\,\hbox {Pa}^{1}\,\hbox {s}^{1}\) (MVN2) for arterioles, and \(9.5\times 10^{8}\,\hbox {Pa}^{1}\,\hbox {s}^{1}\) (MVN1), \(9.1\times 10^{8}\,\hbox {Pa}^{1}\,\hbox {s}^{1}\) (MVN2) for venules. The difference is expected since we estimated lower endfoot sheath permeability.
Although not directly significant for the permeability of the endfoot sheath but relevant for propositions about its main function, we additionally provide cell density estimates resulting from the analysis of the microvascular networks in combination with astrocyte endfoot areas. The assumed diameterarea relations mean that on average, small vessels (\(r_v < 3.0\,\upmu \hbox {m}\)) show an average endfoot density of \(2\times 10^{4}\) endfeet/\(\hbox {mm}^{2}\) surface area. Larger venous vessels (\(r_v > 3.0\,\upmu \hbox {m}\)) show \(1\times 10^{4}\) endfeet/\(\hbox {mm}^{2}\) and larger arterial vessels (\(r_v > 3.0\,\upmu \hbox {m}\)) show the lowest density of \(0.4\times 10^{4}\) endfeet/\(\hbox {mm}^{2}\). McCaslin and coworkers [5] find \(1\times 10^{4}\) endfeet/\(\hbox {mm}^{2}\) for capillaries, \(0.4\times 10^{4}\) endfeet/\(\hbox {mm}^{2}\) for venules, \(0.3\times 10^{4}\) endfeet/\(\hbox {mm}^{2}\) for arterioles invivo in mouse cortex.
Using the endfoot sheath surface areas in Table 1 and our density estimates, we compute about 170,000 (MVN1), 245 000 (MVN2) endfeet around small vessels, 2600 (MVN1), 4600 (MVN2) endfeet around larger venous vessels, and 1300 (MVN1), 3000 (MVN2) endfeet around larger arterial vessels. 97% of endfoot processes are therefore expected to be around capillaries. Using an estimate of astrocyte densities in the mouse cortex (\(20,000 \pm 13,000\) cells/mm\(^3\) [57]) this means the domain of the analyzed networks contains about \(11,000 \pm 7000\) (MVN1) and \(14,000 \pm 9000\) (MVN2) astrocytes with 16 (MVN1) and 19 (MVN2) endfoot processes per astrocyte on average.
Finally, we want to stress that with regard to the prediction of transport across or in parallel to the endfoot sheath, in addition to the presented permeability parameters, a dynamic model for pressure and concentration around vessel networks on the \(\upmu \hbox {m}\) to mm scale (mesoscale) is needed. Concerning implication for macroscale transport models (organscale), we remark that the integral values reported here for the microvascular networks may be used as a starting point to estimate parameters for tissue transport models based on homogenization or mixture theory. However, one should be aware that effective filtration and diffusion capacity on the macroscale generally depend on the local mesoscale pressure and concentration distributions which is an unresolved issue of such models [56] in the context of tissue perfusion simulations.
Relevance in light of the glymphatic theory
Cerebrospinal fluid (CSF) flow through perivascular spaces is a crucial component of the recently proposed glymphatic theory [60,61,62]. Pial perivascular CSF flow has been observed and quantified in [63]. Furthermore, intake of various tracers (Dextran, Gadobutrol) into the parenchyma has been reported to be modified by sleep and disease in both mice and humans [47, 64, 65]. Crucial to determining the mechanisms involved in the intake is to determine the type and magnitude of fluid flow and molecular transport along the different pathways: perivascular, intracellular, and extracellular; and the resistance of barriers between these compartments and the resistance of efflux pathways. Therefore, the herein presented permeability estimates for the astrocyte endfoot sheath being a component of all conceived pathways, provide a starting point for estimating diffusive and advective fluxes outside of the microvasculature.
In [6], the authors estimate the effect of varying astrocyte endfoot gap density on transmembrane CSF flux based on (at least) three assumptions^{Footnote 12}: (1) there is a fluidfilled connected perivascular space (PVS) from descending arterioles all the way down to capillaries; (2) there is a net CSF flow within the PVS from the cortical surface into the capillary bed driven by axial pressure gradients in the PVS; (3) water transport across the endfoot sheath (or transport through intracellular pathways) does not affect the pressure distribution in the PVS, i.e. the exchange is small in comparison to the perivascular flow rates. In a theoretical analysis based on these assumptions, the authors conclude that varying endfoot gap fractions help “maintaining perivascularinterstitial flux through the cortical depth” [6]. The arteriole endfoot sheath has a lower surfacespecific permeability (\(L_p\)) but experiences a larger pressure drop (\(\Delta p\)) across the endfoot sheath than the capillary endfoot sheath (given the authors’ assumptions). These competing effects cancel each other out so that the resulting local fluxes across the endfoot sheath are approximately equal in both vessels. However, a flux is a surfacespecific quantity. We want to additionally point out that for a given portion of tissue (as in Fig 7), since there are many more capillaries than arterioles with a much higher total surface area (Table 1), perivascularinterstitial exchange (even with the authors’ assumptions) would happen predominantly around capillaries.^{Footnote 13} Moreover, the latter statement remains true, even if the arteriole endfoot sheath would have the same (higher) \(L_p\) as the capillaries. However, regardless of this remark, the low permeability (high resistance) of the endfoot sheath in comparison with the ECS may lead to slightly enhanced fluid flow parallel to vessels within the PVS (under the premise that a sufficient driving force and a connected pathway exist).
Limitations
The given permeability estimates have been derived within the wellestablished framework of porous media theory. However, the parameter values used as the basis for our permeability estimates are uncertain and due to a lack of quantitive invivo data, it is currently not possible to quantify the parameter uncertainty. The parameter uncertainty has various origins. The published experimental data is sparse and, for instance, the intercellular gap width, \(d_{\textrm{g}}\), or the gap length \(l_{\textrm{g}}\) have not been characterized invivo. Moreover, the uncertainty inherent to certain experimental techniques (for instance, chemical tissue fixation), providing the quantitative basis for the parameter values used, is not wellcharacterized to date. Additionally, parameter values are likely not constants but undergo dynamic changes in living tissue.
We estimated parameter variance within microvascular networks. For instance, the geometry of MVN1 and MVN2 are different, but overall both the endfoot sheath filtration coefficients, \(L_p\), and the diffusion membrane coefficient \(C_M\) of the whole networks vary only by around 1%. The difference between the layers of different cortical depths is more pronounced. Moreover, the estimated parameter values of the smallest and the largest vessels in the mouse cortical brain tissue differ by one order of magnitude. Finally, parameter values may vary between different species.
As discussed above in detail, the estimated diffusion membrane coefficients, \(C_M\), scale linearly in \(d_g\), while the estimated filtration coefficients, \(L_p\) are proportional to \(d_g^3\). Hence, the exact values estimated should be used cautiously. The values have to be reviewed in case new data becomes available and may have to be adapted for different species. For this purpose, we provide the necessary theoretical considerations and equations, see (2) and (4).
Nevertheless, there is value in estimating parameters to the best knowledge we have today. For example, the estimated parameter values may provide a good starting point for parametric studies with hydraulic network models of water transport in cortical brain tissue.
Conclusion
This work shows how a datainformed theoretical model of astrocyte endfoot size distributions (based on Voronoi tessellations) can be used to relate data from various experimental and theoretical works and arrive at estimates for the endfoot sheath permeability and its variation in microvascular networks from mouse brain cortex. We estimated filtration coefficients in the range \(L_p = 2 \times 10^{11}\,\hbox {m}\,\hbox {Pa}^{1}\,\hbox {s}^{1}\) to \(3\times 10^{10}\,\hbox {m}\,\hbox {Pa}^{1}\,\hbox {s}^{1}\) (average \(2.1 \times 10^{10}\,\hbox {m}\,\hbox {Pa}^{1}\,\hbox {s}^{1}\)) and diffusion membrane coefficients in the range \(C_M = 5 \times 10^{2}\,\hbox {m}^{1}\) to \(6\times 10^{3}\,\hbox {m}^{1}\) (average \(4.4\times 10^{3}\,\hbox {m}^{1}\)). This means that the astrocyte endfoot sheath is more than one order of magnitude more permeable than the vessel wall but about two orders of magnitude less permeable than a similarly thick layer of extracellular space. The numbers are complemented by formulas such that they can be adapted in the case that other data becomes available. In particular, we estimated a relation between the interendfoot gap area fraction and the vessel radius given a constant gap width and find values in the range of 0.2–0.6 %. The data is presented with the intent to be useful for detailed modeling studies of transport of substances in the brain cortex including microvascular network architecture. The estimates are based on the assumption of continuous endfoot coverage of cortical microvessels in mice with an approximately constant interendfoot gap width of 20 nm and largest uncertainty for the permeability of the extracellular pathway stems from the absence of direct evidence of continuous endfoot coverage and the precise geometry of interendfoot gaps invivo.
Availability of data and materials
All data generated or analyzed during this study are included in this published article. Any thirdparty data used is publicly available from the cited sources. The source code for the program used to generate Voronoi tessellations and the source code for the program used to compute pressure distributions in the microvascular networks in the current study are available from the corresponding author on request.
Change history
26 July 2023
A Correction to this paper has been published: https://doi.org/10.1186/s1298702300456x
Notes
This is to match the distribution of endfoot areas on single vessels measured by [6]. For example, a regular uniform cell center pattern would lead to uniform endfoot areas instead.
using the opensource software CGAL [20]
The endfoot processes forming the endfoot sheath overlap. The tessellation models the “visible” endfoot sheath surface configuration as seen from outside the vessel, cf. [6, Fig.2].
L is chosen large (here \(L= 20 r_{\textrm{o}}\)) such that a possible bias on cell size due to boundary effects is minimized. The center points are duplicated along the sides so that the generated pattern is periodic and can be mapped onto (wrapped around) a cylinder surface (as shown in Fig. 1).
Wellmodeled based on a twosample KolmogorovSmirnov test for goodness of fit. We remark that it has previously been observed that the cell area in Voronoi tessellations (with uniformlydistributed center point coordinates) may be approximated by a Gamma distribution, e.g. [36,37,38]; Weaire et al. [37] provide an intuitive explanation.
computed as \(\phi _{\textrm{g}} = A^{1}d \root 4 \of {3}\sqrt{2A}\) based on a single regular hexagon with area A.
For this estimate, the reported tissue permeabilities (3.5\(14\times 10^{18}\hbox {m}^{2}\) [51], 0.4\(110\times 10^{18}\hbox {m}^{2}\) (mean: \(16\times 10^{18} \pm 24\times 10^{18}\hbox {m}^{2}\)) [52], \(1020\times 10^{18}\hbox {m}^{2}\) [53]) are divided by the slab thickness of \(1\,\upmu \hbox {m}\) and fluid viscosity of \(0.69 \times 10 ^{3}\,\hbox {Pa}\,\hbox {s}\). We used a density of \(1 \times 10^{3}\,\hbox {kg}\,\hbox {m}^{3}\) and gravitational acceleration of \(9.81\,\hbox {ms}^{2}\) for conversion from hydraulic conductivity.
In [54], the authors estimate tissue permeability based on data from (whole brain) convectionenhanced delivery studies, that are almost three orders of magnitude larger than what is reported from perfusion studies [51] and flow simulations [53]. However, the estimates might be altered in comparison with values for only ECS by effects of tissue deformation during injection [55] or by the inclusion or opening of highly permeable perivascular pathways.
For molecules of size of a similar order of magnitude as \(d_{\textrm{g}}\) and larger, continuum theory is not applicable.
Actual pressure drops might be much smaller depending on how pressure gradients are generated which is currently unknown.
with \(\Sigma SL_A \cong 1/R_{\textrm{IEG}}\) of [18] and \(V =1\,\hbox {L}\).
We mention that all three assumptions have been challenged and the current evidence does not suffice for a resolution of the debate, see e.g. [61] for a review.
Large penetrating arterioles only account for about 1% of the total microvascular surface area
Abbreviations
 CSF:

Cerebrospinal fluid
 ECS:

Extracellular space
 MVN:

Microvascular network
 MRI:

Magnetic resonance imaging
 MS:

Multiple sclerosis
 PVS:

Perivascular space
References
Wolff J. Beiträge zur Ultrastruktur der Kapillaren in der normalen Grosshirnrinde. Z Zellforsch Mikrosk Anat. 1963;60(3):409–31. https://doi.org/10.1007/bf00336616.
Helmchen F, Kleinfeld D. In vivo measurements of blood flow and glial cell function with twophoton laserscanning microscopy. In: Angiogenesis: In Vivo Systems, Part A. Methods in Enzymology, vol. 444, Academic Press; 2008. p. 231–254. Chap. 10. https://doi.org/10.1016/S00766879(08)028103.
Mathiisen TM, Lehre KP, Danbolt NC, Ottersen OP. The perivascular astroglial sheath provides a complete covering of the brain microvessels: an electron microscopic 3d reconstruction. Glia. 2010;58(9):1094–103. https://doi.org/10.1002/glia.20990.
Watanabe K, Takeishi H, Hayakawa T, Sasaki H. Threedimensional organization of the perivascular glial limiting membrane and its relationship with the vasculature: a scanning electron microscope study. Okajimas Folia Anat Jpn. 2010;87(3):109–21. https://doi.org/10.2535/ofaj.87.109.
McCaslin AFH, Chen BR, Radosevich AJ, Cauli B, Hillman EMC. In vivo 3d morphology of astrocyte—vasculature interactions in the somatosensory cortex: implications for neurovascular coupling. J Cerebral Blood Flow Metab. 2010;31(3):795–806. https://doi.org/10.1038/jcbfm.2010.204.
Wang MX, Ray L, Tanaka KF, Iliff JJ, Heys J. Varying perivascular astroglial endfoot dimensions along the vascular tree maintain perivascularinterstitial flux through the cortical mantle. Glia. 2020;69(3):715–28. https://doi.org/10.1002/glia.23923.
Matsumae M, Sato O, Hirayama A, Hayashi N, Takizawa K, Atsumi H, Sorimachi T. Research into the physiology of cerebrospinal fluid reaches a new horizon: intimate exchange between cerebrospinal fluid and interstitial fluid may contribute to maintenance of homeostasis in the central nervous system. Neurol Med Chir. 2016;56(7):416–41.
Iliff JJ, Wang M, Liao Y, Plogg BA, Peng W, Gundersen GA, Benveniste H, Vates GE, Deane R, Goldman SA, Nagelhus EA, Nedergaard M. A paravascular pathway facilitates CSF flow through the brain parenchyma and the clearance of interstitial solutes, including amyloid \(\upbeta\). Sci Transl Med. 2012. https://doi.org/10.1126/scitranslmed.3003748.
Asgari M, de Zélicourt D, Kurtcuoglu V. How astrocyte networks may contribute to cerebral metabolite clearance. Sci Rep. 2015. https://doi.org/10.1038/srep15024.
MacAulay N. Molecular mechanisms of brain water transport. Nat Rev Neurosci. 2021;22(6):326–44. https://doi.org/10.1038/s41583021004548.
MacAulay N. Reply to ‘aquaporin 4 and glymphatic flow have central roles in brain fluid homeostasis’. Nat Rev Neurosci. 2021;22(10):651–2. https://doi.org/10.1038/s4158302100515y.
Honda H. Description of cellular patterns by dirichlet domains: the twodimensional case. J Theor Biol. 1978;72(3):523–43. https://doi.org/10.1016/00225193(78)903156.
Zisis E, Keller D, Kanari L, Arnaudon A, Gevaert M, Delemontex T, Coste B, Foni A, Abdellah M, Calì C, Hess K, Magistretti PJ, Schürmann F, Markram H. Digital reconstruction of the neurogliavascular architecture. Cereb Cortex. 2021;31(12):5686–703. https://doi.org/10.1093/cercor/bhab254.
Brightman MW, Reese TS. Junctions between intimately apposed cell membranes in the vertebrate brain. J Cell Biol. 1969;40(3):648–77. https://doi.org/10.1083/jcb.40.3.648.
Koch T, Flemisch B, Helmig R, Wiest R, Obrist D. A multiscale subvoxel perfusion model to estimate diffusive capillary wall conductivity in multiple sclerosis lesions from perfusion MRI data. Int J Numer Methods Biomed Eng. 2020. https://doi.org/10.1002/cnm.3298.
Tithof J, Boster KA, Bork PA, Nedergaard M, Thomas JH, Kelley DH. A network model of glymphatic flow under different experimentallymotivated parametric scenarios. iScience. 2022;25(5):104258. https://doi.org/10.1016/j.isci.2022.104258.
Faghih MM, Sharp MK. Is bulk flow plausible in perivascular, paravascular and paravenous channels? Fluids Barriers CNS. 2018;15(1):1–10. https://doi.org/10.1186/s1298701801038.
Vinje V, Eklund A, Mardal KA, Rognes ME, Støverud KH. Intracranial pressure elevation alters CSF clearance pathways. Fluids Barriers CNS. 2020. https://doi.org/10.1186/s12987020001891.
Blinder P, Tsai PS, Kaufhold JP, Knutsen PM, Suhl H, Kleinfeld D. The cortical angiome: an interconnected vascular network with noncolumnar patterns of blood flow. Nat Neurosci. 2013;16(7):889–97. https://doi.org/10.1038/nn.3426.
Yvinec M. 2D triangulations. In: CGAL User and Reference Manual, 5.5 edn. CGAL Editorial Board; 2022. https://doc.cgal.org/5.5/Manual/packages.html#PkgTriangulation2.
Rohatgi A. Webplotdigitizer: Version 4.5; 2021. https://automeris.io/WebPlotDigitizer.
Nicholson C, Hrabětová S. Brain extracellular space: the final frontier of neuroscience. Biophys J. 2017;113(10):2133–42. https://doi.org/10.1016/j.bpj.2017.06.052.
Michel CC, Curry FE. Microvascular permeability. Physiol Rev. 1999;79(3):703–61. https://doi.org/10.1152/physrev.1999.79.3.703.
Renkin EM. Filtration, diffusion, and molecular sieving through porous cellulose membranes. J Gen Physiol. 1954;38(2):225.
Beck RE, Schultz JS. Hindrance of solute diffusion within membranes as measured with microporous membranes of known pore geometry. Biochimica et Biophysica Acta (BBA) Biomembranes. 1972;255(1):273–303. https://doi.org/10.1016/00052736(72)900284.
Herring N, Paterson DJ. Levick’s introduction to cardiovascular physiology. 6th ed. Cham: CRC Press; 2018. https://doi.org/10.1201/9781351107754.
Thomas JA, McGaughey AJH, KuterArnebeck O. Pressuredriven water flow through carbon nanotubes: insights from molecular dynamics simulation. Int J Therm Sci. 2010;49(2):281–9. https://doi.org/10.1016/j.ijthermalsci.2009.07.008.
Schmid F. Averaged results of blood flow simulations with discrete RBC tracking for microvascular networks. Cham: Zenodo; 2017. https://doi.org/10.5281/zenodo.269650.
Douglas DH, Peucker TK. Algorithms for the reduction of the number of points required to represent a digitized line or its caricature. Cartographica: the international journal for geographic information and geovisualization. 1973;10(2):112–22. https://doi.org/10.3138/fm576770u75u7727.
Schmid F, Tsai PS, Kleinfeld D, Jenny P, Weber B. Depthdependent flow and pressure characteristics in cortical microvascular networks. PLoS Comput Biol. 2017;13(2):1005392. https://doi.org/10.1371/journal.pcbi.1005392.
Pries AR, Secomb TW, Gessner T, Sperandio MB, Gross JF, Gaehtgens P. Resistance to blood flow in microvessels in vivo. Circ Res. 1994;75(5):904–15. https://doi.org/10.1161/01.res.75.5.904.
Koch T, Gläser D, Weishaupt K, et al. DuMux 3—an opensource simulator for solving flow and transport problems in porous media with a focus on model coupling. Comput Math Appl. 2021;81:423–43. https://doi.org/10.1016/j.camwa.2020.02.012.
Sander O, Koch T, Schröder N, Flemisch B. The dune foamgrid implementation for surface and network grids. Arch Numer Softw. 2017;5:1–2017. https://doi.org/10.11588/ANS.2017.1.28490.
Hill RA, Tong L, Yuan P, Murikinati S, Gupta S, Grutzendler J. Regional blood flow in the normal and ischemic brain is controlled by arteriolar smooth muscle cell contractility and not by capillary pericytes. Neuron. 2015;87(1):95–110. https://doi.org/10.1016/j.neuron.2015.06.001.
Bonney SK, CoelhoSantos V, Huang SF, Takeno M, Kornfeld J, Keller A, Shih AY. Public volume electron microscopy data: an essential resource to study the brain microvasculature. Front Cell Dev Biol. 2022. https://doi.org/10.3389/fcell.2022.849469.
Kiang T. Random fragmentation in two and three dimensions. Z Astrophys. 1966;64:433.
Weaire D, Kermode JP, Wejchert J. On the distribution of cell areas in a voronoi network. Philosophical Magazine B. 1986;53(5):101–5. https://doi.org/10.1080/13642818608240647.
Koufos K, Dettmann CP. Distribution of cell area in bounded poisson voronoi tessellations with application to secure local connectivity. J Stat Phys. 2019;176(5):1296–315. https://doi.org/10.1007/s1095501902343y.
Korogod N, Petersen CC, Knott GW. Ultrastructural analysis of adult mouse neocortex comparing aldehyde perfusion with cryo fixation. eLife. 2015. https://doi.org/10.7554/elife.05793.
HajYasein NN, Vindedal GF, EilertOlsen M, Gundersen GA, Øivind S, Laake P, Klungland A, Thorén AE, Burkhardt JM, Ottersen OP, Nagelhus EA. Glialconditional deletion of aquaporin4 (Aqp4) reduces bloodbrain water uptake and confers barrier function on perivascular astrocyte endfeet. Proc Natl Acad Sci. 2011;108(43):17815–20. https://doi.org/10.1073/pnas.1110655108.
Kubotera H, IkeshimaKataoka H, Hatashita Y, Allegra Mascaro AL, Pavone FS, Inoue T. Astrocytic endfeet recover blood vessels after removal by laser ablation. Sci Rep. 2019;9(1):1263. https://doi.org/10.1038/s41598018374194.
Mills WA, Woo AM, Jiang S, Martin J, Surendran D, Bergstresser M, Kimbrough IF, Eyo UB, Sofroniew MV, Sontheimer H. Astrocyte plasticity in mice ensures continued endfoot coverage of cerebral blood vessels following injury and declines with age. Nat Commun. 2022. https://doi.org/10.1038/s41467022294752.
Florence CM, Baillie LD, Mulligan SJ. Dynamic volume changes in astrocytes are an intrinsic phenomenon mediated by bicarbonate ion flux. PLoS ONE. 2012;7(11):51124. https://doi.org/10.1371/journal.pone.0051124.
Rosic AB, Dukefoss DB, Åbjørsbråten KS, Tang W, Jensen V, Ottersen OP, Enger R, Nagelhus EA. Aquaporin4independent volume dynamics of astroglial endfeet during cortical spreading depression. Glia. 2019;67(6):1113–21. https://doi.org/10.1002/glia.23604.
Schmid F, Barrett MJP, Jenny P, Weber B. Vascular density and distribution in neocortex. Neuroimage. 2019;197:792–805. https://doi.org/10.1016/j.neuroimage.2017.06.046.
Bojarskaite L, Bjørnstad DM, Vallet A, Binder KMG, Cunen C, Heuser K, Kuchta M, Mardal KA, Enger R. Sleep cycledependent vascular dynamics enhance perivascular cerebrospinal fluid flow and solute transport. BioRxiv. 2022. https://doi.org/10.1101/2022.07.14.500017.
Xie L, Kang H, Xu Q, Chen MJ, Liao Y, Thiyagarajan M, O’Donnell J, Christensen DJ, Nicholson C, Iliff JJ, et al. Sleep drives metabolite clearance from the adult brain. Science. 2013;342(6156):373–7. https://doi.org/10.1126/science.1241224.
Syková E, Nicholson C. Diffusion in brain extracellular space. Physiol Rev. 2008;88(4):1277–340. https://doi.org/10.1152/physrev.00027.2007.
Kimura M, Dietrich HH, Huxley VH, Reichner DR, Dacey RG. Measurement of hydraulic conductivity in isolated arterioles of rat brain cortex. Am J PhysiolHeart Circ Physiol. 1993;264(6):1788–97. https://doi.org/10.1152/ajpheart.1993.264.6.h1788.
Fraser PA, Dallas AD, Davies S. Measurement of filtration coefficient in single cerebral microvessels of the frog. J Physiol. 1990;423(1):343–61. https://doi.org/10.1113/jphysiol.1990.sp018026.
Swabb EA, Wei J, Gullino PM. Diffusion and convection in normal and neoplastic tissues. Can Res. 1974;34(10):2814–22.
Franceschini G, Bigoni D, Regitnig P, Holzapfel GA. Brain tissue deforms similarly to filled elastomers and follows consolidation theory. J Mech Phys Solids. 2006;54(12):2592–620. https://doi.org/10.1016/j.jmps.2006.05.004.
Holter KE, Kehlet B, Devor A, Sejnowski TJ, Dale AM, Omholt SW, Ottersen OP, Nagelhus EA, Mardal KA, Pettersen KH. Interstitial solute transport in 3d reconstructed neuropil occurs by diffusion rather than bulk flow. Proc Natl Acad Sci. 2017;114(37):9894–9. https://doi.org/10.1073/pnas.1706942114.
Smith JH, Humphrey JAC. Interstitial transport and transvascular fluid exchange during infusion into brain and tumor tissue. Microvasc Res. 2007;73(1):58–73. https://doi.org/10.1016/j.mvr.2006.07.001.
Støverud KH, Darcis M, Helmig R, Hassanizadeh SM. Modeling concentration distribution and deformation during convectionenhanced drug delivery into brain tissue. Transp Porous Media. 2012;92(1):119–43.
Vidotto E, Koch T, Köppl T, Helmig R, Wohlmuth B. Hybrid models for simulating blood flow in microvascular networks. Multiscale Modeling Simul. 2019;17(3):1076–102. https://doi.org/10.1137/18m1228712.
Keller D, Erö C, Markram H. Cell densities in the mouse brain: a systematic review. Front Neuroanatomy. 2018. https://doi.org/10.3389/fnana.2018.00083.
Ji X, Ferreira T, Friedman B, Liu R, Liechty H, Bas E, Chandrashekar J, Kleinfeld D. Brain microvasculature has a common topology with local differences in geometry that match metabolic load. Neuron. 2021;109(7):1168–118713. https://doi.org/10.1016/j.neuron.2021.02.006.
Adams DL, Piserchia V, Economides JR, Horton JC. Vascular supply of the cerebral cortex is specialized for cell layers but not columns. Cereb Cortex. 2014;25(10):3673–81. https://doi.org/10.1093/cercor/bhu221.
Iliff JJ, Wang M, Liao Y, Plogg BA, Peng W, Gundersen GA, Benveniste H, Vates GE, Deane R, Goldman SA, Nagelhus EA, Nedergaard M. A paravascular pathway facilitates CSF flow through the brain parenchyma and the clearance of interstitial solutes, including amyloid \(\upbeta\). Sci Transl Med. 2012. https://doi.org/10.1126/scitranslmed.3003748.
Hladky SB, Barrand MA. The glymphatic hypothesis: the theory and the evidence. Fluids Barriers CNS. 2022. https://doi.org/10.1186/s1298702100282z.
Bohr T, et al. The glymphatic system: current understanding and modeling. iScience. 2022;25(9):104987. https://doi.org/10.1016/j.isci.2022.104987.
Mestre H, Tithof J, Du T, Song W, Peng W, Sweeney AM, Olveda G, Thomas JH, Nedergaard M, Kelley DH. Flow of cerebrospinal fluid is driven by arterial pulsations and is reduced in hypertension. Nat Commun. 2018;9(1):1–9. https://doi.org/10.1038/s41467018073183.
Eide PK, Valnes LM, Lindstrøm EK, Mardal KA, Ringstad G. Direction and magnitude of cerebrospinal fluid flow vary substantially across central nervous system diseases. Fluids Barriers CNS. 2021;18(1):1–18. https://doi.org/10.1186/s12987021002516.
Eide PK, Vinje V, Pripp AH, Mardal KA, Ringstad G. Sleep deprivation impairs molecular clearance from the human brain. Brain. 2021;144(3):863–74. https://doi.org/10.1093/brain/awaa443.
Acknowledgements
We would like to thank Franca Schmid for her support in interpreting the microvascular network data.
Funding
Open access funding provided by University of Oslo (incl Oslo University Hospital). This project has received funding from the European Union’s Horizon 2020 Research and Innovation programme under the Marie SkłodowskaCurie Actions Grant agreement No 801133
Author information
Authors and Affiliations
Contributions
TK conceived and designed the study, performed the simulations, analyzed the results, and wrote the article with input from the other authors; VV and KAM helped interpret the results and provide context. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The original version of this article was revised: The Result section of the Abstract has been updated.
Appendices
A Underestimation of endfoot area by 2D image analysis
The analysis of varying endfoot area on the surface of vessels with varying diameter, conducted by Wang and coworkers [6], is based on twodimensional image analysis. A twodimensional image of a vessel is taken. (The image looks similar to Fig. 1 (bottom) as shown in [6, Fig.2].) Interendfoot gaps are segmented and the area surrounded by interendfoot gaps is identified as endfoot area. There are two main errors inherent to the methodology:

(1)
[Projection error] The image shows an orthogonal projection of the (actually curved) vessel surface onto the image plane. Due to the projection the actual area is underestimated. The effect is stronger for endfeet that are close to the vessel “boundary” (vessel outline in the 2D image; in 3D there is no such boundary).

(2)
[Visibility error] The image only shows the part of the endfoot on the visible half of the vessel. The area covered by an endfoot wrapping around the vessel is underestimated.
Importantly, the magnitude of these errors depends on the ratio \(\sqrt{A}/r_{\textrm{v}}\), (endfoot dimension over the vessel radius). The effects are therefore more pronounced for capillaries. Consider the following two limit cases:

(a)
\(\sqrt{A} \ll r_{\textrm{v}}\): A small endfoot on a large vessel surface centered in the middle between the vessel boundaries is fully visible. There is only error (1). This error tends to 0 in the limit of infinitely large vessel radius.

(b)
\(\sqrt{A} \gg r_{\textrm{v}}\): Imagine a rectangular endfoot (one side with arbitrary length) fully wrapped around the vessel (other side equal to \(2\pi r_{\textrm{v}}\)). Due to the projection error, the area is underestimated by a factor \(\frac{2 r_{\textrm{v}}}{\pi r_{\textrm{v}}} = \frac{2}{\pi } \approx 0.64\). Additionally, only half of the endfoot is visible (factor 0.5). In combination, the area is underestimated by a factor 0.32, or the endfoot appears (on the image) approximately 3 times smaller than it actually is.
The presented theoretical model based on Voronoi tessellations of the vessel surface, allows us to investigate the error more systematically. To this end, we make one additional assumption: During the image analysis, seeing an image like Fig. 1 (bottom), the scientist is likely aware that counting small polygons close to the vessel boundary (corresponding to incompletely seen endfoot processes) decreases the accuracy of the results. Here, we cannot be sure how many such polygons have been counted. We therefore consider, in the theoretical analysis, varying thresholds based on the position of the centroid of the endfoot polygon. Specifically, we assume that an endfoot polygon is only counted if its centroid (w.r.t. to its visible projected portion) is in the middle \(P\,\%\) between the vessel boundaries. All polygons centered too close to the boundary on any side (somewhere in the \((100P)\,\%\) percent boundary region) will be omitted from the analysis. Since the value of P has a large influence on the results, we tested different numbers, \((100P) \in \lbrace 0, 10, 30, 50, 70 \rbrace\).
Next, we evaluated for different ratios \(\sqrt{A}/r_{\textrm{v}}\) and different values of \((100P)\), the factor of endfoot area underestimation. To this end, we generated patterns like in Fig. 1 (top), wrapped them around the vessel, cf. Fig. 1 (bottom), and projected them onto the image plane. We did so for different angles and for 5 random realizations per data point. The standard deviation of the computed factors was approximately 0.02 for all cases and is omitted in the following figure only showing the mean value. The results are presented in Fig 8. The results show that there is a significant underestimation inherent to the projection, in particular, for small vessels. For a capillary, with \(r_{\textrm{v}} = 3\,\upmu \hbox {m}\) and \(A = 50\,\upmu \hbox {m}\), projection results in an area estimate between \(15\,\upmu \hbox {m}\) and \(35\,\upmu \hbox {m}\). The upper value is obtained, if only polygons are counted that are centered in the middle 30% of the seen vessel section. That means, even if the best is tried to exclude polygons that are clearly cut at the boundary, a significant error remains.
Our analysis suggests a way to correct the error. To simplify the inverse problem of estimating the corrected area from the measured area, we choose \(P=50\) and observe that the least square fit with a cubic polynomial shown in Fig 8 multiplied with the real area, \(A_{\textrm{img}} = A f(\sqrt{A}/r_{\textrm{v}})\), is monotone in the relevant parameter ranges (\(r_{\textrm{v}} \in \left[ 2.5, 40 \right] \,\upmu \hbox {m}\), \(A \in \left[ 10, 500 \right] \,\upmu \hbox {m}^{2}\)) and can therefore uniquely be inverted. Inverting the relationship \(r_{\textrm{v}}\) and \(A_{\textrm{img}}\) obtained by [6] at discrete sampling point with Brent’s root finding algorithm leads to the corrected diameterarea function shown in Fig. 3.
B Pressure distribution in microvascular networks
The vessels in the microvascular networks were classified into arterial and venous vessels based on the pressure in each segment. The pressure was computed based on the network geometries extracted from measurement data [19] and boundary conditions estimated in [28, 30]. The resulting pressures are visualized in Fig. 9.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Koch, T., Vinje, V. & Mardal, KA. Estimates of the permeability of extracellular pathways through the astrocyte endfoot sheath. Fluids Barriers CNS 20, 20 (2023). https://doi.org/10.1186/s12987023004218
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12987023004218
Keywords
 Brain perfusion
 Astrocyte processes
 Endfoot sheath
 Glial cells
 Perivascular space
 Glymphatic theory
 Waste clearance