1. Introduction
The highly organized living systems are energetically open and far from thermal equilibrium [1]. Its physical phenomena are collective and have strong physical roots [2]. Structures built up by anabolism and store information in the open system [3]. The living matter is heterogeneous, having numerous different electrolytes engulfed by specialized tissues or lipids enveloping isolated aqueous electrolytes in definite structures. The isolating layers control the chemical and physical reactions between the electrolytes and regulate the complex interactions. The fundamental division of electrolytes is between the cytosol (the intracellular electrolyte) and the Extracellular Matrix (ECM). The membrane is a complexly organized multifunctioning part. This double lipid layer regulates the information and ionic exchange between the intra and extracellular reagents, having a vital role in the energy distribution and production of the entire system. The mass of the living object is volume dependent (scaling by 3), while the surface is scaled only by 2. Consequently, we expect an exponent for mass-dependence of
energy exchange (metabolism) as 2/3, the ratio of the cell surface to the cell volume (
). So the expected metabolic power (
) in rest state (Basal Metabolic Power, BMP) dependence vs. mass (M) is expected:
(1)
However, the experiments show a variation of exponents, the
is not shared. When the metabolism is concentrated on surfaces, the
well approaches reality. On the other hand, when it is centered on the energy resources, the exponent is close to
. When the whole mass of the organism is involved in the metabolic energy exchange, the exponent is near to
. In complete demand, the actual body-part (organ or whole-body) needs maximally available energy supply, proportional to its mass, so the scaling exponent is
in this case [4]. In this case, the actual demand decides about the metabolic power and not the geometry. Of course, both the extremes are not ideal for the living object and could not follow evolutional requests. What is optimal? Despite the different exponential power, one feature is strictly common, all the experiments show power-scale (called scaling) in a few orders of magnitudes of the parameters, which is linear in the double logarithmic plot:
(2)
The scaling behavior is the consequence of the self-similarity of the living objects [5] [6]. The fundamental phenomenon behind it is the relative proportional change of the parameters [7]. The fundamental principle was oriented on the changes of the same organism, which has to grow in collective harmony, so the relative growth of parts must have balanced growth [8]. The structure and regulation of biosystems are complex. Various modern approaches have been developed in the last few decades to describe this complexity. The description of statistics of complex systems is far from the normal (Gaussian) distribution. Usually, power-law-tailed distributions (with a general exponent α) are applied:
(3)
There are various phenomena, including social, economic, physical, chemical, and biological, to be described by this function [9] [10] [11] [12] [13]. Despite the somewhat different fields of applications of the power law, it has a common root in complex systems: self-organization. The simplest fingerprint of the self-organized complexity is the self-similar or scale-free structures characterized by a power function. This power-function relation magnifies the
by a constant only, m-dependent
value at any m magnification of x:
(4)
Self-organization explains the evolution of the system [14], expressed in non-linear dynamics [15].
The objective of this present article focuses on analyzing the metabolic alimentation of the healthy tissues in normal conditions and the developing tumors in two different conditions:
1) When the tumor metabolizes as a homeostatic organized unit, the theoretically expected allometric exponent corresponds with the optimal healthy allometry;
2) When the tumor metabolism is not in such an “ideal” optimization of the metabolic supply, its alimentation is suboptimal, using the observed fractal behavior of its angio-structure.
2. Method
Fractal physiology describes the structural and dynamical properties of living organisms and their parts [16] [17], based on physical principles [18]. The self-similar behavior could be described by the normalized relative change of the magnitudes, similarly to the Weber-Fechner law [19] in psychophysics like:
(5)
where α is a constant fitting factor. By integration, we get the (3):
(6)
The self-similar functional relation makes a “scale-invariance” feature due to the independence of the magnification, which is the fundamental behavior of the fractal structures, too [20].
The power function is the central description of the scaling in (1), which bases the allometry of living organisms [21]. The original allometry idea was recognized almost a hundred years ago [22], but the exciting question of the energizing of the life phenomena explained in connection of allometry is a half-century-old knowledge [23]. The connection between the homeostatic energizing level and the basal metabolic rate (B0) as a self-similar function of mass (m) of living objects is [24]:
(7)
where the two parameters are determined experimentally; a is the allometric coefficient, and α is the allometric exponent, and (7) is usually called bioscaling [25]. The usual regression analysis uses the logarithmic transformation of (7):
(8)
which allows high linear accuracy and fits both parameters a and α well [26]. The literature has numerous debates about the theoretical allometric relation based on fractal calculus and the empirical fits based on probability calculus [27]. The B0 of living objects shows allometric scaling to its mass, which refers to the energy supply of the living mass of the volume. The (7) function gives a correct mathematical and biological framework for the complex bio-systems fractal studies [28]. The scaling power function of the mass describes it, and it has been shown valid in a broad category of living structures and processes [29]. The scaling considerations are applied not only in biology but broader, in the complete biosystem as well [30]. The importance of understanding the challenges of the complexity of human medicine was recognized on this basis [31] [32].
In a simple formulation, metabolic processes are surface-dependent, while the mass is proportional to the volume. Therefore, the exponent of their ratio mirrors their dimensionality, and consequently, the exponent is 2/3 [33]. Complex living allometry most likely shows the exponent as 3/4 instead of 2/3 in a broad spectrum of living objects [34], or at least have no linearity in a double-logarithmic plot [35]. However, the large data-mining does not show an overall validity of the 3/4 exponent over 2/3 [36]. The curvature could be size-dependent in developing clusters by their size [37]. The 3/4 exponent could be described as a relation between the three-dimensional surface and the four-dimensional volume [38]. The explanation of the fourth dimension is based on the fractal structure of microcirculation [39], which supplies the energy demand according to a homeostatic equilibrium (B0) in the living complexity. Life in this meaning is “four-dimensional”. Its metabolic exchange processes proceed on fractal surfaces, maximizing the available energy consumption, scaling even the fluctuation of the metabolic power in the universal scaling law as well [40].
The optimization of energy consumption was formulated rigorously by the scaling idea and discussed in a universal frame, even on the energy-consumption subcellular level, including the mitochondria and respiratory complexes [41]. The allometry shows a structural, geometrical constraint for living organisms in homeostatic equilibrium.
The metabolic scaling in cancer development is critical [42]. Contrary to the homeostatic homogeneity of the healthy tissue [43]; the functional heterogeneity of the solid tumor allows an abnormal organ self-possession of multiple cell-types and electrolytes like the Extracellular Matrix (ECM) lymph and blood-transports [44]. The tumor metabolism is based on the blood transport to the tumor. The logarithm of wet-weight of the tumor (
) and the tumor blood-flow (
) have linear dependence [45], which was observed in model xenografts of ovarian cancer, so they have a bioscaling relation:
(9)
where the exponent is close to 3/4.
The allometric scaling supposes three geometrical variables to define the optimization of the circulatory system in living objects:
· The average length of the blood circulatory network (l);
· The surface of the relevant material exchange of the blood circulation system (s);
· The volume of the blood (v).
Furthermore, we suppose that these parameters are represented by the self-similar, self-organized functions of the L value, which is characteristic of a given organ. Hence:
(10)
Using the theoretical fractal explanation, the conditions are:
,
and
, from where:
(11)
where
,
, and
is the characteristic length. The first relation limits the pattern of the circulatory system to the maximum that could be planar, while the second is limited to a maximum, filling up a three-dimensional space. The third exponent
could be calculated because the exponents are not independent. The volume is proportional to the product of the surface and length:
(12)
consequently
(13)
Using these conditions, we obtain from (11):
(14)
Furthermore, the actual volume of the blood is proportional to the actual mass of the given system or organ:
(15)
Considering (14) and (15), now we have:
(16)
The metabolism is a surface-controlled mechanism, so
, consequently:
(17)
If the living structure is geometric in conventional Euclidean meaning, then
and therefore
; consequently, the scaling is
, where
. When at least one of the
,
, which modifies the common simple dimensional approach of the metabolic processes.
3. Results
The allometry gives a possibility to describe the development of the tumor [46]. It is valid for the primary cancer lesions but not always applicable in metastases [47]. We are dealing with primary tumors only. There are two ways of approximating the allometric metabolic rate of a tumor:
1) The theoretical approach accepts that a healthy life has a four-dimensional behavior connected to the highly self-organized, consequently self-similar hierarchic order [48], we fix the exponent to 3/4.
2) The experimental approach assumes that the cell cluster tries to maximize its metabolic rate [49], and this way, it modifies the scaling exponent from the value of 3/4.
Both approaches depend on the environmental conditions of the tumor, mainly on the nourishment of the cells.
3.1. Optimal Alimentation to Maximum Metabolic Rate
Evolution maximized the surface where the nutrients are transferred from the blood to the cells, ensuring the best conditions of the living object, so:
(18)
This task is equivalent to the minimizing of the reciprocal value of the exponent in (16):
(19)
with constraint conditions of:
(20)
(19) can be transformed into
(21)
Hence, considering (20), the minimum condition demands that:
(22)
Substituting (22) into (11), the exponents of the self-similar structures are:
(23)
Consequently, in cases of ideal alimentation, these exponents are the fractal dimensions of the parameters of the network, and while the length is one dimensional, the surface is three, and the volume is four [38]. Because metabolism is a surface-regulated process,
, the scaling exponent of the metabolic rate versus mass using (16) is 3/4:
(24)
and so
in the unit mass
(25)
Primarily the blood stream provides the metabolic supply, so the fractality of the vascular network could be decisional in its allometric evaluation. The condition of (24) maximizes the blood flow energizing all the parts of the volume for their optimum, providing a maximum metabolic rate.
In consequence of (24), the life prefers the large masses as more effective energy-consumers in a unit volume shown in (25). However, on another side this process could lead to the loss of complex information, developing higher instability of the system, arguing that this is a negative tendency manifest the “aging of life’s algorithm as a whole” [50]. The model could be applied by guessing when the energy supply is optimal, so the developed active surface cannot supply the actual demands. Two different sources are possible to create such a situation (1) the length of the supplier system changes (the constructional template differs), or (2) the volume of transport exchange is limited despite the growing demands. Various irregularities originate both challenges could be a symptom of disease, like cancer [51].
3.2. Suboptimal Alimentation for Tumor
The malignancy usually demands a higher energy input from its healthy environment than the available. The tumor supply is suboptimal. The higher energy demand (usually exponential in starting phase [52] ) forces to increase the length of the vessel network. In cancerous clusters, contrary to (22), the vascular fractal dimension (
) of the supplying blood-vessel network (
) is larger than 1, (
) [51]; consequently
in the relation of (11).
could be measured by the box-counting method [51].
According to (11), the actual active surface is evolutionary normal for self-organizing of healthy tissues (
). The extra energy demand of the intensive proliferation changes the exponents of parameters in (11). In this case, the surface of the supply follows the evolution-requested exponent of 3 (
from (22)) in the self-similar conditions, but the requested length changes:
(26)
where
modifies the power of the transport measures, so the fractal organization of the transport lines is different. This type of change could be formed by neoangiogenesis satisfying the higher energy demand in cancerous tissues and could cause abnormalities as inflammation, thrombosis, varicose veins modification of the arteries, etc. The corresponding power-law for the actual metabolic rate at the longer length of vessels, so the suboptimal metabolic rate in this phase (
) from (24) is:
(27)
The apparent “dimension” of the reaction request for volume is
, the dimension increases. According to
the measurable fractal dimension of the blood vessel network is
. In this way the
is measurable by the fractal dimension of the vessel structures [53], for example, with the box-counting method [51]. When
[53],
, and the scaling exponent is
. The Microvessel Fractal Dimension (MFD) (which is equivalent with (
) for renal cell carcinoma ranges between 1.30 - 1.66 [54], and correlates well with the tumor Microvessel Density (MVD) [54]. From (27) we know, when
, the scaling exponent describes a non-fractal-like structure,
.
When the tumor growth is so intensive that the available length of the vessel network cannot deliver appropriate energy, then another possible deviation from the homeostatic self-organization happens. In this case, the volume of the delivered energy remains constant, which limits the energy supply. The tumor-growth turns to sigmoidal this stage [55], usually follows Weibull distribution due to the self-similar development [6]. This could happen in severe hypoxia, low oxygen saturation in blood, anemia, various hematological diseases. In this case, the volume of the supply follows the evolution-requested exponent of 4 (
) (23), in the self-similar conditions in [53], but the requested length and surface is not enough for the proper work, so
and
. The self-similar conditions differ from (22) due to (13):
(28)
Consequently, at fixed four-dimensional volume, the metabolic surface reactions behave by power-law of suboptimal metabolic rate in this phase:
(29)
Here the volume “dimension” of the reaction request is 4, but the actual conditions are worse than optimal. The
again here also is measurable by the fractal dimension of the structures [56], in this case, the fractal dimension of the vessel system is
. For example, measuring the vascular fractal dimension in one disease as
[51], we use
, so the scaling exponent is
. When
[56], the scaling exponent is
. At
, the scaling exponent is the well-known
.
The exponents of the active transport surface in the two suboptimal supplies
(30)
Both these exponents are smaller than the optimal, and the exponent in the second phase of growth is the smallest (Figure 1 ).
![]()
Figure 1. The allometric development with the various exponents. The initial phase of tumor growth is exponential, and the allometric scale follows the phase1 curve, while the intensive development further decreases the exponent, which turns to a sigmoid phase in a tumor-specific time.
Following the idea of “4-dimensionality”, the volume dimension changes in
, while the active surface in
. Both exponents are ≤3/4 because of the length fractal dimension >1. The mechanisms which cause this modification are different. The same exponent could be only in Euclidean non-fractal case, when
. Due to
and
, hence
and
. The Mandelbrot calculated 2.7 for the fractal dimension of the arterial tree of the lung [57], which was supported by experiments later [58], the relevant changes are
and
. Both values are smaller than
.
Both non-optimal situations (defect of the length of transport way or limited transport against the demands) make the tissue under-energized, and the exponent of the power-relation scaling down-regulated. In such a way, measuring the scaling exponent of metabolism and the fractal dimension of the supplying microvessels have a diagnostic value about the actual deviations from normal.
4. Discussion
The optimal alimentation in a healthy system makes the energy distribution balanced, supplying all requirements of the homeostatic state. The exponent 3/4 has a strong predominance on a theoretical and empirical basis [59] in healthy homeostatic basal metabolic activity. The ideal nutrition supply supports ontogenic growth. However, at least at larger sizes, the cancer growth never happens with an optimal nutrition supply; the cells compete intensively for the available energy sources.
The cancer is out of the overall homeostatic balance. The tumor development certainly has a higher energy supply due to its proliferation than its healthy counterpart needs. Due to the extra-large energy demand, the tumor development’s alimentation in most cases is far from optimal, so the tumor is in a permanent energy deficiency. When the oxygen supply is limited, the first attempt to produce more ATP is the massive fermentative use of glucose, a simple and quick production mechanism. The cell extends its ATP production to fermentation by non-mitochondrial respiration, abandoning the more complicated Kerbs-cycle in the mitochondria [60].
While the mitochondrial metabolism is always aerobic, its scaling exponent is nearly
[30]. However, the scaling of metabolic activity is also different in mitochondrial and non-mitochondrial processes [61]. The metabolic power not only depends on the active surface of the transport but also on the transport rate at the same active surface size. Due to the transport modifications at the changed metabolic pathway, the deviation from the 3/4 exponent could be remarkable. The allometric scaling exponent of fermentative processes decreases to nearly 2/3. This last scaling exponent shows that the cell-membrane directly regulates the fermentation, and the surface/volume ratio controls the complete process, which could be anticipated from the direct linear dependence of the lactate production (
) on the glucose-intake (
) with a slope of ≅1 [45]:
(31)
while the bioscaling of the oxygen (
) and glucose (
) intake [45] are even lower than 2/3, −0.570, and −0.523, respectively.
Not only does the malignancy need an intensive extra metabolism. For example, the benthic invertebrates (n = 215) have the lowest average scaling exponent (
, [near to 2/3],
), which metabolizes in an anaerobic way [62]. No regulative factor exists when the cells are entirely independent, and the available alimentation is unlimited (like in most in vitro experiments). The metabolic rate is linearly proportional to the mass, so the exponent is
[30].
The metabolic transformation of the cells [63] is one of the well-recognized hallmarks of malignancy [64] that has an emerging intensive interest in the field of oncology [65], as the core hallmark of cancer [66]. The adaptation of mitochondria in energy-limited conditions is the focus of the research [67]. The tumor forces the development of the angiogenetic processes [68] and overcomes the energy limitations. The vascularity is promoted [69], and the rapid development by intensive proliferation supports the changes of the scaling behavior [70]. Without extra angiogenesis (starting clusters), only the ready-made capacity of the delivery is available, so the tumor has a suboptimal alimentation. In the beginning, its fractal structure was developed, which is similar to the healthy structure, so the four-dimensional scaling remains valid (28).
When the tumor develops, the fractal structure of vascularity changes. Consequently, its fractal dimension changes too. The forced angiogenesis [68] tries to provide a sufficient supply to the hypoxic (insufficiently supported) tumor, and the structure changes rapidly, broadening the scaling exponent in a wide range [71]. The missing supply suppresses the scaling exponent, shown in (27). Still, the angiogenetic pool changes the trend, approaching linearity. The unlimited availability of nutrients for every cell realizes the linearity measured in vitro [48], limited to ~0.9 in vivo by insufficient oxygen transport [72]. However, the angiogenesis is usually not fast enough to supply the faster-growing larger tumors, so the inner part of the tumor becomes necrotic, forming a smaller living mass to supply, easing the energy distribution [68]. The essential message of the cases of insufficient alimentation from the calculations above is that when the fractal dimension of the supplying network grows, the scaling exponent decreases. The four-dimensionality and the allometry with the evolutional optimizing request are not the same approaches: further evolution conditions have a higher than four-dimensional allometric scaling. The tumor mass is a somewhat indefinite parameter because the whole environment of the tumor suffers from suboptimal alimentation. Consequently, we tried to find a more fundamental networking condition parameter published elsewhere [73].
There is a vast number of researches about the vascular development of the tumor progression, calculating the fractal dimension of the vascularity. The in silico modeling of the growing tumor vessel architecture in high-grade gliomas [74] shows that the fractal dimension is less than 1 in the avascular state and growing linearly by time, reaching
at
, by slope approximately
. In a longer time, the development of the fractal dimension drastically changes, follows a less rapid development (slope
) until
. We may assume that the fractal dimension 1.2 characterizes the finally developed vessel structure inside the tumor, followed by neo-angiogenetic processes reaching the tumor-surface, changing the vascular architecture, growing slower to the higher values of the fractal dimension.
In optimal alimentation, the allometric scaling shows exponent 3/4 (24); which supposes the
, so the vascular fractal dimension in this case is
. However
by the growing vessel network, so
in the allometric scaling of tumor-vascularity due to the suboptimal energy supply, which triggers the angiogenesis. Using the results from in silico model-calculations, the internal
growth of the vessels have
, or
,
according to the assumption of suboptimal alimentation by maximal metabolic rate (case 1) or by the metabolic rate forced four-dimensional “optimizing” concept (case 2). When the external angiogenesis is developing, the allometry
changes:
, and
. So, the optimizing
of the suboptimal energy availability in extended angiogenetic cases realizes the allometry, which fits the simple geometrical expectations
well.
The measurements of vascular fractal dimensions in various tumors show a lower scaling exponent than the ideal 3/4, depending on the conditions of the tumor-angiogenesis development. For example, when the epithelial-connective tissue interfaces with a malignant tumor in the oral mucosa, it is
[75], the scaling exponent in suboptimal alimentation situations in cases 1 and 2 are
and
. The last one (optimal distribution of the suboptimally available energy) is near the “conventional” 2/3. Another microscopic evaluation of angio-structures [76] shows lower values of
, like the fractal dimension of the normal and malignant tissues are
and
, respectively [77]; resulting in low
values. In other evaluations, the vascular structure’s dimensionality grows to 1.9, which provides the maximal energy usage of the suboptimal alimentation, and the exponent became as low as
.
It is interesting to see the effect of various anti-tumor treatments on the vascular fractal dimension. The treatment changes the vascularization and suppresses the fractal dimension forms 1.135 1.037, 0.933, 0.982 by Photodynamic Therapy (PDT); Cysteine Proteases Inhibitors (CPI), combined therapy, PDT and CPI [78]; which corresponds in cases when the maximalizing of the energy-supply is equivalent to the allometric exponents of 0.716, 0.741, 0.767 and 0.755, respectively (the optimal distribution of the suboptimal availability would be 0.726, 0.743, 0.763, and 0.753). By treating VEGF165, the fractal dimension increases from 1.65 to 1.69, decreasing the allometric exponent [79] [80]. In matrigel inoculated human umbilical vein endothelial cells (HUVEC) treated by docetaxel, the fractal dimension of the vascular structure has decreased from 1.2 to 1.09, corresponding in case 1
, and in case 2
[81]. The fractal analysis is a successful and rather accurate method for monitoring the efficacy of angiogenic consequences of therapies [82].
(a) (b)
Figure 2. The summary of the structure of calculation. (a) The biophysical considerations (b) The mathematical description.
5. Summary
We had shown that the allometric relation of basal metabolic rate and the tumor mass depends on the fractal dimension of the vascular structure. Due to the desperate need for energy supply and the intensive proliferation of the malignant tumor, cancer does not have an optimal alimentation. Two strategies could distribute the available (not sufficient) energy by the main transport of it, the vascular network:
1) Assuming that the cell cluster tries to maximize its metabolic rate by the surface transports and lowers the scaling exponent from the value of 3/4;
2) Accepting that in the case of a four-dimensional volumetric behavior limits the energy supply. The tumor optimizes the energy distribution in its volume in among these conditions.
The structure of the biophysical considerations and their mathematical steps are summarized in Figure 2 .
The two strategies in consequent phases of tumor growth optimize the available energy by different allometric scalings. The organized optimum of the suboptimal availability of energy gives lowered allometric scaling exponents.
Acknowledgements
This work was supported by the Hungarian National Research Development and Innovation Office PIACI KFI grant: 2019-1.1.1-PIACI-KFI-2019-00011.