Journal of Applied Mathematics

Volume 2013 (2013), Article ID 396486, 6 pages

http://dx.doi.org/10.1155/2013/396486

Research Article

Analysis of a Model for the Morphological Structure of Renal Arterial Tree: Fractal Structure

Aurora Espinoza-Valdez,¹ Francisco C. Ordaz-Salazar,² Edgardo Ugalde,³ and Ricardo Femat⁴

¹Departamento de Ciencias Computacionales, CUCEI, Universidad de Guadalajara, Avenida Revolución 1500, CP 44430, Guadalajara, JAL, Mexico

²Universidad Politécnica de San Luis Potosí, Urbano Villalón 500, CP 78363, La Ladrillera, SLP, Mexico

³Instituto de Física UASLP, Universidad Autónoma de San Luis Potosí, Avenida Manuel Nava No. 6., Zona Universitaria, CP 78200, SLP, Mexico

⁴Laboratorio para Biodinámica y Sistemas Alineales, División de Matemáticas Aplicadas, IPICYT, Apartado Postal 3-90, CP 78231, Tangamanga, SLP, Mexico

Received 7 March 2013; Accepted 12 June 2013

Academic Editor: Kiwoon Kwon

Copyright © 2013 Aurora Espinoza-Valdez et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

One of the fields of applied mathematics is related to model analysis. Biomedical systems are suitable candidates for this field because of their importance in life sciences including therapeutics. Here we deal with the analysis of a model recently proposed by Espinoza-Valdez et al. (2010) for the kidney vasculature developed via angiogenesis. The graph theory allows one to model quantitatively a vascular arterial tree of the kidney in sense that (1) the vertex represents a vessels bifurcation, whereas (2) each edge stands for a vessel including physiological parameters. The analytical model is based on the two processes of sprouting and splitting angiogeneses, the concentration of the vascular endothelial growth factor (VEGF), and the experimental data measurements of the rat kidneys. The fractal dimension depends on the probability of sprouting angiogenesis in the development of the arterial vascular tree of the kidney, that is, of the distribution of blood vessels in the morphology generated by the analytical model. The fractal dimension might determine whether a suitable renal vascular structure is capable of performing physiological functions under appropriate conditions. The analysis can describe the complex structures of the development vasculature in kidney.

1. Introduction

The arterial structure of organs has been the subject of many studies [1–5]. However, not all systems have similar functions; for the organs the purpose of arterial structure is to provide the blood required in the metabolic process and other specific functions. The arterial structure is highly nonuniform because it is determined by reasons of anatomy and local flow requirements [4]. The kidney is one of the most complicated organs in terms of structure and physiology because it is highly vascularized and constitutes the main organ for maintaining chemical balance in blood [6]. The kidney consists of three trees: arterial, venous, and ureter. The arterial vascular tree of the kidney is structured by the renal artery branches into interlobar arteries, arcuate arteries, and interlobular arteries, which are formed by bifurcation [6]. The development of the arterial vascular tree of the kidney can be formed mainly by angiogenesis [1]. The process of angiogenesis is the formation of new blood vessels from preexisting vessels and consists of two different processes: sprouting and splitting. Sprouting refers to the case in which the new branch literally sprouts to some existing branch. In splitting a branch is divided into two new branches [7]. The factors involved at the vascular development are VEGF, renin-expression, ephrins A and B, platelet-derived growth factor-B (PDGF-B), Ets family (as Ets-1 and TEL), and angiopoietins 1-2. The VEGF is an essential regulatory factor for both processes of angiogenesis in the development of renal arterial vascularization [1]. Then, we only consider VEGF, because the other factors increase the complexity of the analytical model.

Zamir and Phipps in [8] studied the morphological characteristics of the rat kidney and found that the branching rules of these vessels are determined by considerations of the angiological function. Subsequently, Zamir generated a model based on L-systems (Lindenmayer System) and incorporated some physiological laws, that is, random parameters [4]. Zamir's results suggest that the arterial structure of the tree is determined mainly by the growth rules of arterial tree branching [4]. The structural morphologic reconstruction of renal vasculature from microcomputed tomography (micro-CT) images was presented by Nordsletten et. al. The arterial and venous trees of the rat kidney were generated numerically using micro-TC [3]. These morphological data provide a statistical basis from which renal vascular topology can be generated [5]. Graph theory generates branching tree structures incorporating the physiological laws of the renal artery branching through the process of angiogenesis [5]. Many pathological conditions, such as cancers, arteriovenous malformations, and diabetes, induce changes to vessel's morphology or spatial organization. The fractal analysis has been applied to a large variety of healthy or pathological vascular networks [9–14]. Sabo et al. concluded in [15] that the microvessel fractal dimension as a marker of tumor microvascular complexity might provide important pronostic information as well as shed light on the complex interactions between tumor angiogenesis and growth. Moreover, Cross et al. analyzed the fractal dimension of renal angiograms: a normal kidney, a congenitally dysplastic kidney, a kidney with renal artery stenosis, and a kidney with recurrent thromboembolism lesions (see Figure 3 in [16]). The analysis of fractal dimension may allow for characterization of the vascular morphology.

In this paper we extend the results of [5] by analysing mathematically the structure of the arterial vascular tree of the kidney and by incorporating physiological parameters generated from a computational model. In the implementation of the proposed model the branches are represented by the edges of the graph and the ramification point is represented by the nodes [5]. We expect that this representation implemented in computer methodology will help to understand the mechanism involved in biological systems, specifically the development of the renal arterial bifurcation. The analysis provides the mean of the length in each level, and the average width of the arterial vascular tree of the rat kidney is in the range of what is observed in previous experimental and computational models. The fractal dimension depends on the sprouting angiogenesis that appears in the development of the arterial vascular tree of the kidney. The complex structures of the renal arterial tree are consistent with previous fractal analysis of vascular structures [9, 11–14]. These study opens the possibility of a new taxonomy for normal kidneys and for the pathological injuries related to the vascular morphology of the kidneys.

2. Theoretical Basis

The representation of the renal arterial vascularization in the computational model is as follows: the arterial vascular tree of the kidney is defined by labeled and oriented graphs; we can identify two objects: the branches (edges) and the branching point (nodes). That is, the edge represents a blood vessel and the vertex is where the angiological stimulus appears. The length and diameter of the new branches are smaller than their parent branches [5, 6].

The arterial vascular tree of the kidney down to the interlobular arteries can be structured as follows: renal artery, level 0: interlobar arteries, levels 1-2: arcuate arteries, levels 3-4: interlobular arteries, and levels 5–9. We based in the experiments of Tomanek [3] for the morphological structure of the renal arterial tree. In the model, we develop the arterial vascular tree of the rat kidney. Therefore, the model based on the graph theory for the arterial tree is a binary tree of ten levels and has a total of vessels equal to 1023. The physiological parameters, and rules are based on parent vessel branch level, the parameters, and the morphological structure of their child branches, and this process deploys the structure of the model. The parameters included in the edges of the binary tree model are ,  ,  ,  , and   as follows (for details see [5]).

(i)The Type of Angiogenesis: Sprouting or Splitting . Angiogenesis is identified by the variable which can take the values in the set , where denotes sprouting angiogenesis and denotes splitting angiogenesis. Both sprouting and splitting angiogeneeis generate branching tree structures, which can be represented using graph theory. It is conjectured that the sprouting and splitting processes occur with different probability in the renal arterial tree for the morphology of the kidney. In the proposed computational model this assumption is controlled probabilistically. Sprouting angiogenesis is more probably than splitting angiogenesis because it requires only reorganization of existing endothelial cells (not migration), that is, having less energy requirements [1].(ii)Concentration of the VEGF in the Vessel . Both processes, sprouting and splitting angiogeneses, depend on the regulation of in the pre-existing vessel [1]. While in sprouting angiogenesis it is known that endothelial cells promotes the differentiation, migration, proliferation, and assembly, the effect on splitting angiogenesis is still in research. In the model is generated randomly in a uniform way within the range  ng/mL reported in [5, 17].(iii)Length of the Vessel . Here, we focus on determining the value of the lengths in the arterial bifurcations depending of the concentration of the VEGF. Although the effect of over length for the splitting angiogenesis is unknown, in the sprouting angiogenesis we consider that the length of the new vessels depends of the . We obtain the value of the dimensionless length from experimental data; the function which can be approximated for VEGF₁₂₁ is [5] . The   is applied in this model; the length of the vessel depends on and is compared with the length ranges reported in the literature according to the level [3, 5]. The formed by splitting is defined by a contraction factor. In the case that is out of these ranges, the process is repeated finite number of times. If is still out of range, then it is adjusted to the length of its parent branch.(iv)Diameter of the Vessel . In this work, we do not discuss the value of the diameters in the blood vessels. This is due to the fact that we do not have sufficient information of in for sprouting and splitting angiogeneses. Then, in sprouting angiogenesis, the diameter of the sprout is adjusted to a factor of the diameter of its parent branch. In splitting angiogenesis the diameter of every new vessel is adjusted to one-half of the diameter of its parent branch.(v)Angle of Bifurcation . In sprouting angiogenesis the angle of bifurcation is adjusted within the range with respect to its parent branch [1, 18]. In splitting angiogenesis, one branch is adjusted to and the other branch to [18].

The physiological parameters are included in the label of the edge and have the form . The arterial vascular tree of the kidney has vertices with oriented edges in such a way that from each vertex one edge arrives and two edges leave (the orientation symbolizes the circulation of blood flow in arteries of the kidney). The renal arterial vascular tree of the kidney has 10 levels, because we based in the experiments of Tomanek. Therefore, the arterial vascular tree of the kidney has nodes. The set of integers to name every node is .

An algorithm generates step by step the physiological characteristics of every branch and saves it as labeled edges with the format . Just as the tree is generated, the other structure is constructed in which we save the position of every node in . The algorithm calculates the position of every node in function of the physiological characteristics of the branches and the position of nodes of its parent branch. The physiological parameters of the root (the initial condition) are determined: with probability and , where ,  ng/mL, , and they continue in numerical order with the other branches using the information of parent branches. Some characteristics such as and are generated probabilistically in every step, and the others depend on these parameters, their parent branch, and the rules that were enumerated in the previous section.

3. Main Results

3.1. Length and Width

The length of each level into the kidney depends on and ; that is, the length for each level denoted by in the arterial vascular tree of the kidney is analytically determined as follows:

where the parameters are defined as in the previous section. Then, the mean length for each segment is

These data are within the range of rat trees reported in [3]; we compare the results with Nordsletten from depth to .

Moreover, the analytical average width denoted by is the root length (renal artery, ) to the leaves of the tree (interlobular arteries, ), and then

where ;   is the length of the renal artery [5]; is defined in the set and ; is the average of the contraction factor for ; and is the contraction factor for . Then,

These data are within the range of rat trees reported in [3]; while we consider 10 levels [1, 5], Nordsletten considers 11 levels.

However, if , the average with respect to is as follows:

Let and , and then

where

Similarly, , and then

where  mm is the error.

Substituting the values ,  ,  , and  ,

We generate 5000 trees of renal arterial vasculature for different probabilities in and . In Figure 1 the mean ± SD (standard deviation) width of the kidney has an intersection with these data and previous result where the mean width is  mm [3].

Figure 1: The mean ± SD of width of 5000 trees generated for 5 different probabilities of and .

Therefore, the depends on the process of angiogenesis; that is, is inversely proportional to the probability of the .

3.2. Fractal Dimension of Kidney Vascular Tree

The fractal dimension quantifies through dimension the ability of an object to occupy a space. Different methods exist and are based on the different definitions of fractal dimension [9, 11–14, 19].

The analytical fractal dimension can be derived using the following result [19]. Let be a hyperbolic iterative function system (IFS), and let denote its attractor. Suppose that is a similitude of scaling factor for each . If the IFS is totally disconnected or just touching, the attractor has fractal dimension , which is given by the unique solution of

If the IFS is overlapping, then , where is the solution of

The morphology of arterial vascular tree of the kidney, produced by an IFS, is composed of applications of contraction with a factor . Thus, from the subsection of length and width, the contraction factor for the splitting angiogenesis is while for sprouting angiogenesis the average contraction factor is . Since we have , the analytical fractal dimension of the renal arterial tree is

where is generated randomly within the range defined by . In other words, the fractal dimension depends on the number of vessels, the spatial relationships between the vascular components, and the surrounding environment.

3.2.1. Box Counting of the Arterial Vascular Tree of the Kidney

Additionally, the box counting theorem by Barnsley [19] states that, denoting with as the number of boxes of side ,

where the arterial vascular tree of the kidney has fractal dimension . The box counting results allow one to compute the rate of change in complexity with scale as well as a measure of hetreogeneity.

Here, we used the box counting theorem for determining the fractal dimension. Figure 2 shows the average of the fractal dimension in each probability for the arterial vascular tree of the kidney. The fractal dimension decreases with the increase in the probability of sprouting in the development of the arterial vascular tree of the kidney, whereas the mean width of the kidney decreases. Then, the fractal dimension of the arterial vascular tree of the kidney depends on the probability of and , that is, the distribution of blood vessels in the morphology generated by graph theory model.

Figure 2: The mean ± SD of fractal dimension for different probabilities of , where has the probability .

The results suggest that the fractal dimension is inversely proportional to . This behavior is congruent in the sense that, by symmetry structure, the arterial vascular tree is bigger for small values of ; that is, the renal space is covered more efficiently. The complex structures of the vasculature in kidney are consistent with previous studies of vascular structures about their fractal nature [9–16].

Examples derived from graph model using different probabilities in and , mean width and fractal dimension, are shown in Figure 3. The analysis of these responses may allow for characterization of the vascular morphology on the arterial vascular tree of the kidney.

Figure 3: Trees generated with different probabilities in and (where ), mean width and fractal dimension of the kidney.

4. Conclusions

The arterial vascular tree of kidney development was analyzed through angiogenesis. The renal arterial vascular tree of the kidney goes down into the interlobular arteries. We generate 5000 trees of renal arterial vasculature for different probabilities on and . The analytical mean length in each level and the average width of the arterial vascular tree of the kidney have an intersection with previous studies. Then, the graph theory allows the vascular tree to model the vascular growth; that is, it generates the ramification of structures arborescent incorporating physiological laws of arterial branching.

In conclusion, the analytical expression of the fractal dimension depends on the number of vessels, the spatial relationships between the vascular components, and the surrounding environment. The arterial vascular arterial tree of the kidney has a fractal dimension which is inversely proportional to the probability of the occurring sprouting angiogenesis. The fractal dimensions determined for the development of the arterial vascular tree of the kidney by box counting may allow for characterization of the vascular morphology. As a conjecture, the fractal dimension might determine whether a suitable renal vascular structure is capable of performing physiological functions under appropriate conditions (hemodynamics). These studies could be expanded to include those pathologies originating from arterial alterations.

References

R. J. Tomake, Assembly of the Vasculature and Its Regulation, Birkhäuser, Berlin, Germany, 2001.

E. M. Wahl, L. V. Quintas, L. L. Lurie, and M. L. Gargano, “A graph theory analysis of renal glomerular microvascular networks,” Microvascular Research, vol. 67, no. 3, pp. 223–230, 2004.

D. A. Nordsletten, S. Blackett, M. D. Bentley, E. L. Ritman, and N. P. Smith, “Structural morphology of renal vasculature,” American Journal of Physiology, vol. 291, no. 1, pp. H296–H309, 2006.

M. Zamir, “Arterial branching within the confines of fractal L-system formalism,” Journal of General Physiology, vol. 118, no. 3, pp. 267–275, 2001.

A. Espinoza-Valdez, R. Femat, and F. C. Ordaz-Salazar, “A model for renal arterial branching based on graph theory,” Mathematical Biosciences, vol. 225, no. 1, pp. 36–43, 2010.

A. C. Guyton and J. E. Hall, Tratado de Fisiología Médica, McGraw-Hill, 10th edition, 2000, in Spanish.

S. Patan, “Vasculogenesis and angiogenesis as mechanisms of vascular network formation, growth and remodeling,” Journal of Neuro-Oncology, vol. 50, no. 1-2, pp. 1–15, 2000.

M. Zamir and S. Phipps, “Morphometric analysis of the distributing vessels of the kidney,” Canadian Journal of Physiology and Pharmacology, vol. 65, no. 12, pp. 2433–2440, 1987.

Y. Gazit, D. A. Berk, M. Leunig, et al., “Scale-invariant behavior and vascular network formation in normal and tumor tissue,” Physical Review Letters, vol. 75, no. 12, pp. 2428–2431, 1995.

Y. Gazit, J. W. Baish, N. Safabakhsh, M. Leunig, L. T. Baxter, and R. K. Jain, “Fractal characteristics of tumor vascular architecture during tumor growth and regression,” Microcirculation, vol. 4, no. 4, pp. 395–402, 1997.

J. Panico and P. Sterling, “Retinal neurons and vessels are not fractal but space-filling,” Journal of Comparative Neurology, vol. 361, no. 3, pp. 479–490, 1995.

P. G. Vico, S. Kyriacos, O. Heymans, S. Louryan, and L. Cartilier, “Dynamic study of the extraembryonic vascular network of the chick embryo by fractal analysis,” Journal of Theoretical Biology, vol. 195, no. 4, pp. 525–532, 1998.

C. Arlt, H. Schmid-Schönbein, and M. Baumann, “Measuring the fractal dimension of the microvascular network of the chorioallantoic membrane,” Fractals, vol. 11, no. 2, pp. 205–212, 2003.

S. Lorthois and F. Cassot, “Fractal analysis of vascular networks: insights from morphogenesis,” Journal of Theoretical Biology, vol. 262, no. 4, pp. 416–433, 2010.

E. Sabo, A. Boltenko, Y. Sova, A. Stein, S. Kleinhaus, and M. B. Resnick, “Microscopic analysis and significance of vascular architectural complexity in renal cell carcinoma,” Clinical Cancer Research, vol. 7, no. 3, pp. 533–537, 2001.

S. S. Cross, R. D. Start, P. B. Silcocks, A. D. Bull, D. W. K. Cotton, and J. C. E. Underwood, “Quantitation of the renal arterial tree by fractal analysis,” Journal of Pathology, vol. 170, no. 4, pp. 479–484, 1993.

M. N. Nakatsu, R. C. A. Sainson, S. Pérez-del-Pulgar et al., “VEGF₁₂₁ and VEGF₁₆₅ regulate blood vessel diameter through vascular endothelial growth factor receptor 2 in an in vitro angiogenesis model,” Laboratory Investigation, vol. 83, no. 12, pp. 1873–1885, 2003.

E. Gabryś, M. Rybaczuk, and A. Kȩdzia, “Fractal models of circulatory system. Symmetrical and asymmetrical approach comparison,” Chaos, Solitons & Fractals, vol. 24, no. 3, pp. 707–715, 2005.

M. Barnsley, Fractals Everywhere, Academic Press, 1988.