About the Authors:
Long Jian Liu
* E-mail: [email protected]
Affiliation: Department of Physics, University of Windsor, Windsor, Ontario, Canada
Stephen L. Brown
Affiliations Department of Physics, University of Windsor, Windsor, Ontario, Canada, Department of Radiation Oncology, Henry Ford Hospital, Detroit, Michigan, United States of America
James R. Ewing
Affiliations Department of Radiation Oncology, Henry Ford Hospital, Detroit, Michigan, United States of America, Department of Neurology, Henry Ford Hospital, Detroit, Michigan, United States of America
Brigitte D. Ala
Affiliation: Department of Radiology, Windsor Regional Hospital, Windsor, Ontario, Canada
Kenneth M. Schneider
Affiliation: Department of Radiation Oncology, Windsor Regional Hospital, Windsor, Ontario, Canada
Mordechay Schlesinger
† Deceased.
Affiliation: Department of Physics, University of Windsor, Windsor, Ontario, Canada
Introduction
Despite significant strides in the treatment of solid tumors in tissues such as breast and prostate, local control of other cancers such as primary brain tumors and pancreatic cancer remain challenging. The determinants of local tumor control are complex and multifactorial including microscopic spread beyond traditionally accepted clear margins at surgical resection, inherent radio-resistance and chemo-resistance in both the primary and ‘adjuvant’ (i.e postoperative) use of these therapies and a wide range of tumor and host biological characteristics that are factors in the growth and potential control of malignant tumors. As new therapies are developed and medicines become increasingly individualized, a need arises for an effective and early predictor of response. Tumor interstitial fluid pressure (TIFP) is a physiological parameter with demonstrated predictive value for a tumor’s aggressiveness, as well as response to chemotherapy and radiotherapy [1–7]. However, its clinical use has been limited because the technique available for its measurement has been invasive.
Because of its clinical importance, a noninvasive estimate of TIFP, if such could be developed, might have a profound effect on the assessment and treatment of many solid tumors. By using GdDTPA contrast enhanced MRI, Hassid et al. [8] found that an elevated TIFP of human lung tumors in the central region is associated with a decreased concentration of contrast agent. This shows that the elevated TIFP hinders drug delivery, and may be an indication that a noninvasive measurement of TIFP can be developed. However, to date there is no quantitative relation between TIFP and tissue contrast agent concentration. We will examine the flow of tumor exudate in and around its vascular sources and establish a theoretical framework for a noninvasive estimate of TIFP.
We had previously suggested a positive linear relation between TIFP in the central region and fluid velocity at the tumor surface [9]; an index of TIFP can be formed by noninvasively measuring the fluid velocity at the tumor surface. Supporting the importance of this measurement is a finding [7] that primary tumors in metastasis-positive mice displayed higher IFP and fluid velocity than those in metastasis-negative mice. Similarly, the same investigators found that the fluid velocity at the tumor periphery in cervical cancer patients with pelvic lymph node metastases was faster than fluid velocity in cancer patients without lymph node involvement. The fluid velocity at the tumor surface may be used as a biomarker of tumor aggressiveness [7]. It has been proposed [7, 9] that tumor exudate fluid flow velocity has a negative linear relation with the distance from the tumor surface in the region outside the tumor. Studies have shown that TIFP decreases steeply in the periphery [1, 10], and that fluid velocity increases from zero centrally to a maximum at the tumor surface, in accordance with Darcy’s law, which posits a relationship between fluid velocity in a porous medium and the local gradient of interstitial pressure.
The purpose of this paper is to present a practical model that may be applied for a quantitative noninvasive measurement of TIFP. Experiments (see below) show that the fluid flow within the tumor rim is not uniform, even for a round shaped tumor. A theoretical basis for a TIFP noninvasive measurement is formulated for locally one-dimensional geometries often encountered in practice. The variation of interstitial fluid flow in the tumor rim is described and visualized through MRI images. The relationship between TIFP, velocity, maximum exudate distance, and fluid conductivity is derived, and the noninvasive measurement of these variables with DCE-MRI is demonstrated. This work provides a theoretical basis for the very useful empirical relationship proposed by Hompland et al. [7] to relate the change with time in the position of the point of peak contrast outside a tumor.
Materials and Methods
Model description
We divide a tumor into three regions: 1) The central region, which contains a necrotic core. 2) Tumor periphery region. In this region, tumor blood vessels are plentiful. This region is the main source of tumor interstitial fluid. 3) The intermediary region [9].
We assume that there are no (or little) functional lymphatics or drainage vessels in the central region, but a functional drainage system exists at the periphery. If the drainage ability within the tumor (mainly from tumor periphery) is greater than or equal to fluid source (JS≤JLmax), the net fluid flow at the tumor surface is zero. In this case, the interstitial fluid flow is limited to the tumor and TIFP at the surface (and outside it) is equal to the pressure of the surroundings, which is usually near atmospheric pressure in normal tissue. This case is defined as that of an “isolated” tumor [9, 10]. The following case also belongs to an isolated tumor: there are no functional drainage vessels such as lymphatics surrounding the tumor and no way for fluid to flow out. This condition may only exist in a malignant tumor and TIFP might be very high and could even reach the extreme value pm [9], a maximum pressure barrier. An “embedded” tumor is defined as the case in which the tumor is enclosed by normal tissue and JS>JLmax, i.e., fluid exits the tumor surface, Rs, and flows into the normal tissue. Lymphatic vessels (and/or other drainage channels) are numerous and functional in normal tissue, forming a functional drainage system. One defines the maximum position Rm that tumor interstitial fluid can spread to as the place at which the pressure becomes the same as that of normal tissue.
The interstitial fluid pressure p0 within the central region is uniform in the steady state, as both we and others have described [1, 9, 10, 11], and by Darcy’s law the fluid velocity is zero. The boundary of the central region, R0, is an iso-pressure and iso-velocity surface. In the intermediary region, there is almost no tumor fluid source and only a drain for tumor fluid to be absorbed or normalized. This gives the boundary conditions at the interface, Rm. TIFP is the same as the pressure p∞ in the normal tissue or surrounding environment; tumor interstitial fluid velocity is zero there.
Though we assign R0, Rs and Rm three parameters to represent the interfaces of various regions, the interfaces do not necessarily have to be spherical in shape as shown in Fig 1. In practicality, most tumors are irregular in shape and the curvature of a tumor varies locally. Also, tumor vascular vessels and their permeability are not uniform. Therefore, even a spherically shaped tumor may not give a spherically symmetrical distribution of TIFP. Usually, one surface area may have a large curvature, while another may have a small curvature that is close to zero. Similarly, the vasculature is distributed mostly in the periphery of the tumors [10] and is an irregular distribution of fluid sources. To take the fluid source networks generated by the arteriole-venule pairs as locally one-dimensional may be more practical and realistic, especially when the width of the vascularized peripheral region is relatively narrow compared to the overall dimensions of the tumor.
[Figure omitted. See PDF.]
Fig 1. Schematic of fluid flow in the regions close to tumour surface.
A tumour boundary can be chosen so that it is locally one-dimensional (the figure on the left is from Fig 2B-4 in Ref. 7; it shows a T1-weighted MRI image of a TS-415 cervical carcinoma). S represents the base of the fluid flux in the local coordinate system in which a 1D model can be applied.
https://doi.org/10.1371/journal.pone.0140892.g001
We have derived the distributions of TIFP and fluid velocity in a one-dimensional case in the appendix. Fig 2 shows TIFP distribution for various vascular and lymphatic conditions. For a tumor with poor drainage system (e.g. malignant tumor), the solution for TIFP distribution is approximated as:(1)
Correspondingly, the fluid velocity is expressed as:(2)(3)where dm = Rm-Rs, and d0 = Rs-R0, p∞ is the pressure of the surrounding tissue, usually zero relative to the atmosphere, K and K′are the high hydraulic conductivities of interstitium in the normal tissue and tumor respectively. All the variables are measurable. In the case of an “isolated” tumor as we discussed above, no exudate crosses the tumor surface. The following conditions are satisfied: p(R0) = p0, dp(R0)/dr = 0 and p(Rs) = p∞ (dm = 0). Thus, we have(4)(5)(6)
Based on Eqs 3 and/or 6, once we know the hydraulic conductivity, the fluid velocity and the size of the tumor rim, we can estimate TIFP.
[Figure omitted. See PDF.]
Fig 2. Comparison of TIFP distribution as a function of position (depth) at steady state for various β.
Other parameters are the same based on Eq 18 in appendix.
https://doi.org/10.1371/journal.pone.0140892.g002
A contrast agent is commonly used to study interstitial fluid flow. As a tracer, it reflects (visualizes) interstitial fluid flow. We assume that the contrast agent does not affect interstitial fluid flow and instead becomes a part of the fluid. Interstitial fluid flow is relatively steady for a tumor at a specific stage (e.g. the fluid velocity at a specific point today is almost the same as that of tomorrow). However, when the interstitial fluid with contrast agent flows from tumor surface Rs to an arbitrary position r through time t continuously, the fluid velocity becomes the same at that position, u(t) = u(r). For the sake of convenience, we let s = r(t)−Rs which is defined as the distance that the contrast agent travels from the tumor surface. Therefore, we have (for the tracer):(7)(8)
The circle that is composed of the brightest points in a specific direction within the tumor rim at various times can be determined. Therefore, both the u(Rs) and dm can be estimated through curve fitting. The derivative of Eq 8 is(9)
It shows that the fluid velocity decreases exponentially. When t = 0, r(t) = Rs, s(t) = 0, u(t) = u(Rs); t→∞, r(t) = Rm, s(t) = dm, u(t) = 0. It is consistent with the description above.
Noninvasive measurement of tumor interstitial fluid velocity based on DCE-MRI images
Eq 9 indicates a method for measuring the interstitial fluid velocity at tumor surface u(Rs). Tumor surface can be determined by find the circle which is composed of the brightest points in various directions before contrast agent is applied or within a short time after contrast agent is applied. We can then measure the displacement of circle after a certain time t delay. Through curve fitting based on Eq 8, we can find u(Rs).
For a specific tumor, the data can be measured from a chosen direction. After the bolus administration of a contrast agent, an estimate of the peak position versus time can be fitted to estimate the value of dm and u(Rs), as shown in Fig 3. Also, the position of R0 can be aligned based on the signal intensity vs. position, as shown in Fig 3(A). Fluid within R0 is still and DCE-MRI has the least intensity of pixels. Similarly, we may be able to locate the position of Rm. Rs is determined by find the brightest point before contrast agent is applied or within a short time, which is long enough for visualizing the tumor rim, after contrast agent is applied. In fact, Rs may also be determined by curve fitting.
[Figure omitted. See PDF.]
Fig 3. Method to determine dm and u(Rs).
(a) Signal intensity vs. position in a tumor periphery at 15, 80, 110 and 450 seconds respectively. (b) Curve fitting for determining dm and u(Rs) for a tumor with low IFP (lower) and a tumor with high IFP (upper) [7].
https://doi.org/10.1371/journal.pone.0140892.g003
Measurement of TIFV based on phase contrast imaging
Phase contrast imaging is a tool for visualizing fluid flow. Spins that are moving along the direction of magnetic field gradient causes a phase shift proportional to their velocity: ϕ = γMu + ϕ0 = γG(T/2)2u + ϕ0, where γ is the gyromagnetic ratio, M the momentum, G represents the magnetic gradient, T the bipolar gradient duration, and ϕ0 is the background phase effect due to susceptibility of field inhomogeneity [12]. It requires two repetitions of the sequence, the second of which used bipolar gradients of opposite polarity to the first. Measuring two data sets, we can have the velocity: u = Δϕ/(γΔM). By comparing the phase of signals from each location in the two sequences the exact amount of motion induced phase change can be determined to have a map where pixel brightness is proportional to spatial velocity. Regions that are stationary remain black while moving regions are represented as gray to white.
However, the method based on a velocity contrast sequence may assume that the concentration of the contrast agent is spatially and temporally uniform. Actually, the agent concentration may change with time and the gray value is correlated with the agent concentration. Note that even though the gray value is proportional to the fluid velocity, the proportional coefficient may relate to the concentration of the contrast agent. Therefore, even though the fluid velocity is the same, the gray value may vary if the concentration of the contrast agent is different. Usually, the concentration of contrast agent is time dependent. Therefore, the gray value is also time dependent as shown in Tables 1–3. However, in the same slice with a certain time delay, the concentration of the contrast agent is fixed. Thus, the velocity of the interstitial fluid is proportional to the average gray value in a specific direction within the tumor rim. Through post-processing, phase contrast imaging can give maps of fluid velocity streamlines [13]. As we know that the interstitial fluid velocity at the tumor surface is maximal. This way we can find the fluid velocity at tumor surface u(Rs) in a specific direction. Measuring fluid velocity noninvasively is an important step for actualizing TIFP noninvasive measurement.
[Figure omitted. See PDF.]
Table 1. The average gray-value of the rim of tumor 1 in liver 1 in different directions at various times.
https://doi.org/10.1371/journal.pone.0140892.t001
[Figure omitted. See PDF.]
Table 2. The average gray-value of the rim of tumor 2 in liver 1 in different directions at various times.
https://doi.org/10.1371/journal.pone.0140892.t002
[Figure omitted. See PDF.]
Table 3. The average gray-value of the rim of tumor 1 in liver 2 in different directions at various times.
https://doi.org/10.1371/journal.pone.0140892.t003
Results
DCE-MRI images of tumors in human livers were obtained for study and analysis. We can covert the pixels to millimeters. The length in mm per pixel is determined by measuring the known length of scale in pixel as shown in Fig 4. The contrast agent visualizes the rim variation at tumor periphery. Three tumors (tumors 1 and 2 in liver 1, and tumor 1 in liver 2) are measured as shown in Fig 5. The tumor rim (Rs-R0) or Rm-R0 can be measured although the accuracy may be improved by using a proper method and specifying and standardizing the process. The position of R0 can be aligned based on the signal intensity vs. position. Fluid within R0 is still and DCE-MRI has the least intensity of pixels. Similarly, we can locate the position of Rm, where the intensity of the signal is the least again. The gray value—distance (pixel) distribution is measured by using the software “ImageJ”. Although the images are too coarse (~0.8mm/pixel) for measuring the r(t) accurately, Fig 5 clearly shows that this method is practical for determining the positions of R0, r(t) and Rm, which are vital parameters for TIFP noninvasive measurement.
[Figure omitted. See PDF.]
Fig 4. Schematic for determining the unit length in mm per pixel.
https://doi.org/10.1371/journal.pone.0140892.g004
[Figure omitted. See PDF.]
Fig 5. Liver tumors with a contrast enhancing rim visible.
Several measurements of the rim were made along tumor periphery to infer TIFP. A judicious local coordinate system may be chosen, in which a 1D model can be applied. The rim of tumor 1 in liver 1 is approximately 3.2 mm in 4 pixels, ~0.8mm/pixel.
https://doi.org/10.1371/journal.pone.0140892.g005
We examined the time delay MRI images without and with contrast agent in 0s (without contrast agent), 30s, 60s, 90s, 120s, 150s, 5mins, and 10 minutes for tumors 1 and 2 in the same host liver but different slices. Tables 1 and 2 show the average gray value within the tumor rim at different times for tumors 1 and 2 respectively. For tumor 1, the images with 30s delay do not show a clear tumor rim. The gray-value variation cannot be measured, and it is ignored in Table 1.
The MRI images clearly show that the contrast between the rim and the background decreases as time increases from 150s to 10min. The change rate of concentration of contrast agent should be correlated with the fluid velocity. Therefore, the concentration of contrast agent at different time may predict the fluid velocity.
We also measured the average gray value along the rim in specific directions of tumor 1 in liver 2. The results are listed in Table 3.
However, the data in Tables 1 and 2 show that the gray values decrease in the first 60s after contrast agent is applied. The gray value then increases, though some fluctuation occurs for data in Table 1. The change in gray value may indicate the variation of the contrast agent concentration.
For tumor 1 in liver 2, contrast starts from A5. However, the rims could not be seen from A1 to A5 (until the images with 60s delay after contrast agent is applied). Also, the rims cannot be seen with 5 and 10 min delay. From Table 3, we can see that the gray value of the tumor rim increases from 60 s to 120 s, and then it decreases at 150 s.
The experiment data show that the fluid flow within the tumor rim is not uniform, even though the shape of the tumor is fairly round. This is because the leakage and permeability of tumor blood vessels are not uniform. Therefore, a spherical model usually cannot be applied to predict the interstitial fluid velocity distribution and the pressure distribution. In this case, as we have presented, the one-dimensional model under a judicious local coordinate system can be applied for approaching the tumor interstitial fluid pressure in the central region.
TIFP correlates with hydraulic conductivity. Boucher et al. [14] developed a method for estimating the hydraulic conductivity of tumor interstitium in vivo by measuring the radial pressure difference based on a spherical model. However, as we showed above, the conditions for a spherical model are usually not satisfied, even for a spherically shaped tumor. Hompland et al. [7] experimentally showed the linear relation between fluid velocity u(Rs) and TIFP based on linear regression analysis of their experimental data. The slope k can be found for various tumors and contrast agents as shown in Table 4.
[Figure omitted. See PDF.]
Table 4. The slope of u(Rs)~p0 for various tumors and contrast agents.
https://doi.org/10.1371/journal.pone.0140892.t004
The slope is described as for an isolated tumor (dm = 0) and for an embedded tumor (, if the hydraulic conductivities in the periphery region from R0 to Rs within the tumor and the intermediary region from Rs to Rm are approximately the same K = K′). The result shows that the ratio of the hydraulic conductivity to the rim size is nearly the same for different tumors. It is related to the types of contrast agents. For a tumor with a 2 mm rim, the hydraulic conductivity of tumor interstitium is around 1.30x10-6 cm2/mmHgs for Gd-DTPA and 1.50x10-6 cm2/mmHgs for Gadomelitol. We should choose a proper contrast agent that can best represent interstitial fluid flow. Contrast agent Gd-DTPA gives a better result.
Usually, the hydraulic conductivity in normal tissue is steady since the structure of normal tissue is stable. The hydraulic conductivity within a certain type of normal organ (tissue) should be almost always the same. However, the interstitium structure within the tumor may change with the progression of the tumor stage. Thus, the hydraulic conductivity within the tumor may depend on the tumor and its stage. Even so, based on our method, we can estimate Rm-R0, dm and d0. The fluid velocity can also be found. Usually, the hydraulic conductivity K and the fluid pressure p∞ are known in a specific normal organ tissue. For instance, IFPs of normal tissues are usually as follows: -8mmHg in the lungs, -3 or -2 mmHg in subcutaneous tissues, 0 to +2 mmHg in the liver and kidneys, and +6 mmHg in the brain [15]. Therefore, the hydraulic conductivity K′ and TIFP within the tumor can be estimated. Generally, a malignant tumor with a higher TIFP causes a larger dm. Also, a malignant tumor may have a larger hydraulic conductivity. A larger necrotic core may lead to a smaller d0. Therefore, the slope k should not change much. Once we find the fluid velocity u(Rs) and the pressure in the host tissue, we can estimate TIFP.
For the data shown in Fig 3(B), we fit the two cases with the following two functions:
shigh(t) = 0.60(1−e−0.0053t/0.6) for the high TIFP, and slow(t) = 0.39(1−e−0.0023t/0.39) for the low TIFP case. Based on the data in Fig 3(A), we estimate that Rm-R0 = 1.86mm. If dm = 0.60mm, then d0 = 1.26mm; if dm = 0.39mm, then d0 = 1.47mm. The fluid velocity u(Rs) is 5.3x10-3mm/s and TIFP is approximately 40.8 mmHg for the upper curve case; u(Rs) is 2.3x10-3mm/s and TIFP 17.7 mmHg for the lower case in Fig 3(B). It may be possible to directly measure dm. We may approach the results more accurately by combining this with curve fitting.
Boucher et al. [1] applied Baxter and Jain’s mathematical model [10] to fit one set of their experimental readings for subcutaneous mammary adenocarcinoma. Though there are only two explicit parameters (α and Pe, see Eq 8A of Ref. 9), more parameters are needed for determining these two. Also, the tumor radius R (we use Rs) must be known in order to do curve fitting. The two explicit parameters are theoretically given. They are difficult to determine in practical application. However, when Eq 4 is employed to fit the experimental points from Fig 8 of Boucher et al. [1], only three parameters are needed. Most importantly, these three parameters are measurable. This is essential for noninvasive measurement. Curve fit to the data in Fig 6 is obtained, with an estimate of the distance (Rs-R0) of 0.85mm. In this case, we get K′ = 5.53×10−7 cm2/mmHgSec and u(Rs)≈1.33μm/s.
[Figure omitted. See PDF.]
Fig 6. Comparison of the model for an ‘isolated’ tumor based on Eq 4 with the experimental data from Fig 8 in Ref. 1.
Only three parameters are needed and they are measurable.
https://doi.org/10.1371/journal.pone.0140892.g006
Though many models take the hydraulic conductivities inside and outside a tumor as two different constants [10, 11, 16, 17], some take them to be the same [18, 19]. The hydraulic conductivity in the periphery region and that in the intermediary region should not differ from each other significantly. In this case, based on Eq 3, TIFP can be estimated from the following:(10)
Similarly, we can fit the experimental data from Fig 3A (track 5) of Ref. 1 well (for a mammary adenocarcinoma t.i. tumor). We choose these two sets of readings because they distribute completely from tumor central region to the periphery, and the readings are clear and can be digitalized accurately. Also, the readings from Fig 6 (track 4) in Ref. 1 almost represent the average effect of the data in this figure. The results are shown in Fig 7. The results of both Figs 6 and 7 show that our model can estimate the hydraulic conductivity K′ of a tumor or its IFP in a practical manner. Noninvasive tools such as contrast enhanced MRI can be applied to determine the size of the region from R0 to Rm (or Rs) and the highest fluid velocity u(Rs) as well as dm, and then used to estimate the hydraulic conductivity K′ of the tumor interstitium and TIFP.
[Figure omitted. See PDF.]
Fig 7. Comparison of our model for ‘embedded’ tumors with the experimental data from Figs 3A track 5 (lower) and 6 track 4 (upper, which represents approximately the average value of all readings from all tracks) in Ref. 1.
The readings from these two tracks are distributed across all regions of the tumors.
https://doi.org/10.1371/journal.pone.0140892.g007
Assuming that the conductivity within the periphery region is the same as that within the intermediary region, we find K′upper = 8.45x10-7 cm2/mmHg, K′lower = 7.80 x10-7 cm2/mmHg, and uupper(Rs) = 4.16x10-3mm/s, ulower(Rs) = 1.63x10-3 mm/s based on the data in Fig 7 by curve fitting. It shows that our model gives us information regarding the interstitial fluid flow and interstitium in the rim, e.g. the hydraulic conductivity K′, dm, d0 and u(Rs). This provides a method for investigating tumor conditions at various stages.
Discussion
The one-dimensional (1D) model may also be applied to approximate TIFP and fluid flow for a spherically symmetrical tumor when the ratio of the tumor rim to its radius dm/R is small. For example, theoretically, when dm/R0<5%, the relative error of the fluid flow (caused by the increase in surface area) is less than 10% when using the 1D approximation. Usually, the tumor periphery is around 2mm. A 1D model can give a good approximation for a tumor and ΔR~±1mm.
There are some other noninvasive methods for measuring the velocity of interstitial fluid. For instance, Chary and Jain demonstrated the measurement of interstitial convection by fluorescence photobleaching [20]. Mostly, these measurements can only give the average fluid velocity at the tumor periphery. However, as having described above and in Ref. 9, for an isolated tumor, the fluid velocity at tumor surface is approximately twice of the average fluid velocity since . In fact, when the interstitial fluid velocity of both isolated and embedded tumors is linearly distributed as
, the average velocity is . Therefore, the interstitial fluid velocity at the tumor surface may be estimated.
The hydraulic conductivity has a strong effect on TIFP, just as fluid velocity does. Aside from the methods for determining the hydraulic conductivity that we discussed above, we may also be able to determine the hydraulic conductivity indirectly by measuring other parameters such as the glycosaminoglycan (GAG) concentration, CGAG (g/100g tissue), and tissue water content z. It was found that the hydraulic conductivity correlates with GAG concentration as K = 4.6×10−13CGAG−1.202 [21, 22, 23]. The hydraulic conductivity in the normal host tissue is usually known. Thus, we may find the hydraulic conductivity in the tumor by comparing the GAG concentrations. Another method is based on a relationship between the hydraulic conductivity and the water content since K = azb, where a and b are constants that relate to interstitium conditions [21, 24]. Volumetric water content can be measured by using a near-infrared spectroscopy [25].
When there is no other drain (e.g. lymphatic system) and the fluid is enclosed, TIFP will be uniform and will reach an equilibrium state under this extreme condition. Thus, the fluid flow rate in such a region becomes zero (JA = JV). We have (see Eq 17 in appendix)(11)
When SALA = SVLV, the equilibrium value of TIFP is(12)
Tumors have elevated interstitial fluid pressure due to high permeability of tumor blood vessels and poor drainage system. In fact, a poor drainage system may cause more serious problems. When the drainage system malfunctions, a small leakage source may cause high fluid pressure. If there are no venules or other drainage system such as lymphatics, the highest TIFP can reach is pm = pA−σA(πA−πAi). For a tumor, the surface area of blood vessels in a unit volume S/V and their conductivity L are much larger than that in normal tissue [10]. Also, the oncotic pressure is elevated [26]. These are the factors that contribute to an increase in TIFP. As an example, here we cite the numerical values in Ref. 27: pA = 35mmHg, pV = 15mmHg, πA = πV = 28mmHg. For normal tissue: πAi = 0.1 mmHg, πVi = 3mmHg [27], σN = 0.91 [10]. Thus, pm≈9.6 mmHg and p(TIFP)eq = 0.9 mmHg. For a tumor, σT = 0.82 [10], when the osmotic pressure of interstitial fluid πAi and πvi, which should be approximately equal, are both elevated to approximately 20 mmHg [26], we get pm≈28.4 mmHg and p(TIFP)eq = 18.4 mmHg.
High TIFP may collapse or create malfunctioning blood vessels in tumors. We assume that there are some blood vessels with vascular fluid pressure pv1, pv2, …pvm, where pvm is the highest vascular pressure in the tumor. When the interstitial fluid pressure increases to a value that is higher than pv1, there is a positive net pressure difference Δp = p-pv1>0 (Here we ignore the osmotic pressure) between the interstitial fluid and the blood plasma. In this case, blood vessel with vascular pressure pv1 cannot contribute to interstitial fluid (no fluid can come out from this blood vessel). Similarly, when p>pvi (1<i<m), all blood vessels with pv<pvi cannot contribute to interstitial fluid and will collapse or lose their function. The higher the TIFP, the more vessels will lose their function. This process will repeat until TIFP equilibrium peq has been reached. In this case, the amount of fluid leaked out equals the amount of interstitial fluid drained away. Depending on the condition of the tumor, peq can take on any value between 0 and pvm. If peq reaches pvm, no fluid can leak out from blood vessels within the tumor. This extreme condition gives maximal TIFP limit, which should be pVm−σV(πV−πVm).
Conclusion
The absolute value of TIFP can be approached by converting theoretical variables into measurable ones. We have presented a fairly realistic analysis of tumor exudate flow in tissue, with emphasis on the eventual application of a technique to noninvasively estimate TIFP. It is practical to select a specific direction for the approach to measuring TIFP, and it is viable for tumors of various shapes. However, velocity at the tumor surface as a measure of TIFP depends on other properties of the tumor, and particularly on its surround, with the parameters Rm-R0 (or Rs-R0) and tissue conductivity K of primary importance in defining the relation between TIFP and u(Rs). These parameters can be measured. We provide a means of estimating tumor fluid velocity, and show that the ratio of fluid conductivity of a tumor to its rim size is independent of tumors. It shows that TIFP noninvasive measurement can be actualized. Also, this work provides a method for investigating tumor conditions at various stages by measuring the hydraulic conductivity and fluid flow within the tumor rim.
Appendix
Mathematical model derivation
Interstitial fluid flow in a tumor is formed by the sum of the net flow from various abnormally leaky blood vessels. Generally, it is expressed as [10,16]:(13)where Js/V (fluid source) is the volumetric net flow rate for fluid coming out of blood vessels, which is based on Starling’s law [28–30], and JL/V (fluid drain) the volumetric flow rate for fluid being reabsorbed into lymphatics.
Lymphatic vessels are present in most tissues and their function is to drain the interstitial fluid. Even though there are no lymphatics in the brain, there may be a similar system for draining the fluid [31]. The drainage term is usually described as [10]:(14)where p is the interstitial fluid pressure, pL is the pressure in lymphatic vessels. We will use the subscripts ‘L’, ‘A’ and ‘V’ represent the corresponding counterparts of the lymphatic, arteriole and venule vessels. L represents the hydraulic conductivity of the vessels (m2s/kg), S the surface area of the vessels. Usually, arterioles and venules are distributed in pairs. Normally, the circulatory system processes 20 liters of blood on average per day through capillary filtration. Roughly 17 liters of the filtered plasma is reabsorbed into the blood vessels, and the remaining 3 liters in the interstitial fluid returns to the blood through the lymph system [32]. The fluid lost from, and gained by, plasma is balanced in normal tissue.
The microvascular sources of fluid follow Starling’s law [10,16] as follows.(15)where JA is the fluid flow rate for blood vessels. The flow rate of fluid that is reabsorbed by venules is expressed as:(16)
To simplify the situation, we take the arteriole and the venule as a pair and state the net fluid flow rate:(17)where pA (or pv) is the vascular plasma pressure, σA (σv) the osmotic reflection coefficient, π A (or π v) the osmotic pressure of the plasma, and πAi (πvi) the osmotic pressure of interstitial fluid [10], S/V the surface area of the vessels in a unit volume.
In many cases, a tumor boundary can be chosen so that it is locally one-dimensional. In this case, a 1D model can better describe the actual situation and easier for application. Combining Starling’s law with Darcy’s law (see appendix) for 1D case u(r) = -Kdp(r)/dr, we have:(18)where K is hydraulic conductivity of interstitium, (an index that reflects the strength of fluid source absorption), , and
The variables are primed (′) to signify that they relate to the fluid flow within the tumor.
The boundary and continuity conditions are: p(R0) = p0, u(R0) = -K′dp(R0)/dr = 0, p(Rs+) = p(Rs-) = p(Rs), Kdp(Rs+)/dr = K′dp(Rs-)/dr = -u(Rs), p(Rm) = p∞ (usually it is around 0 relative to the atmosphere since it is at the interface of normal tissue), u(Rm) = −Kdp(Rm)/dr = 0.
Under the boundary and continuity conditions, Eq 18 can be solved. Now the theoretical variables have been converted into measureable variables. The solution of Eq 18 becomes(19)where d0 = Rs-R0 and dm = Rm-Rs. Correspondingly, the tumor interstitial fluid velocity is:(20)
Here u(Rs), d0 and dm can be measured by administering a contrast agent that extravasates and then travels at the same velocity as the tumor exudate. Based on Eq 19, we obtain p(Rs) and p0,(21)
When and are small, Eqs 19–21 reduce to Eqs 1–3.
With a loss of functioning lymphatics, the effective surface area of the drainage vessels becomes smaller and hence β becomes smaller. Thus, the value of β reflects tumor aggressiveness. The fluid flow (and the fluid velocity) at the tumor surface is determined by the fluid source and sink within the tumor. When there is no sink, even a slight vascular leakage may cause high TIFP p0, which can reach pm (0<p0<pm) [9]. Many tumors, especially the malignant ones, may not have a functional drainage system within.
Acknowledgments
The work was funded by a grant from Seeds 4 Hope, a program that provides seed money to foster local cancer research projects through the Windsor Essex County Cancer Centre Foundation.
Dr. M. Schlesinger, the principal investigator of the Seeds 4 Hope project which in part supported the TIFP studies reported, tragically passed away days before this article was accepted for publication. Dr. Schlesinger will be dearly missed for his mentorship, mathematical expertise, passion for scientific discovery, and his drive to contribute to the fight against cancer.
Author Contributions
Conceived and designed the experiments: LJL SLB JRE BDA KMS MS. Performed the experiments: LJL SLB JRE BDA KMS MS. Analyzed the data: LJL SLB JRE BDA KMS MS. Contributed reagents/materials/analysis tools: LJL SLB JRE BDA KMS MS. Wrote the paper: LJL SLB JRE BDA KMS MS.
Citation: Liu LJ, Brown SL, Ewing JR, Ala BD, Schneider KM, Schlesinger M (2016) Estimation of Tumor Interstitial Fluid Pressure (TIFP) Noninvasively. PLoS ONE 11(7): e0140892. https://doi.org/10.1371/journal.pone.0140892
1. Boucher Y, Baxter LT, Jain RK. Interstitial pressure gradients in tissue-isolated and subcutaneous s: implications for therapy Cancer Res. 1990; 50: 4478–4484.
2. Lunt SJ, Kalliomaki TM, Brown A, Yang VX, Milosevic M, Hill RP. Interstitial fluid pressure, vascularity and metastasis in ectopic, orthotopic and spontaneous tumors. BMC Cancer. 2008; 8: 2. pmid:18179711
3. Ferretti S, Allegrini PR, Becquet MM, McSheehy PM. Tumor interstitial fluid pressure as an early-response marker for anticancer therapeutics. Neoplasia. 2009; 11: 874–881. pmid:19724681
4. Salnikov AV, Iversen VV, Koisti M, Sundberg C, Johansson L, Stuhr LB, et al. Lowering of tumor interstitial fluid pressure specifically augments efficacy of chemotherapy. FASEB J. 2003; 17: 1756–1758. pmid:12958200
5. Rofstad EK, Ruud EB, Mathiesen B, Galappathi K. Associations between radiocurability and interstitial fluid pressure in human tumor xenografts without hypoxic tissue. Clin Cancer Res. 2010; 16: 936–945. pmid:20103667
6. Yeo SG, Kim JS, Cho MJ, Kim KH and Kim JS. Interstitial fluid pressure as a prognostic factor in cervical cancer following radiation therapy. Clin Cancer Res. 2009; 15: 6201–6207. pmid:19773374
7. Hompland T, Ellingsen C, Ovrebo KM, Rofstad EK. Interstitial fluid pressure and associated lymph node metastasis revealed in tumors by dynamic contrast-enhanced MRI. Cancer Res. 2012; 72: 4899–4908. pmid:23027087
8. Hassid Y, Furman-Haran E, Margalit R, Eilam R, Degani H. Noninvasive magnetic resonance imaging of transport and interstitial fluid pressure in ectopic human lung tumors. Cancer Res. 2006; 66: 4159–4166. pmid:16618737
9. Liu LJ, Brown SL, Ewing JR, Schlesinger M. Phenomenological model of interstitial fluid pressure in a solid tumor. Physics Review E. 2011; 84: 021919–1–021919–9.
10. Baxter LT, Jain RK. Transport of fluid and macromolecules in tumors. I. Role of interstitial pressure and convection. Microvasc Res. 1989; 37: 77–104. pmid:2646512
11. Baxter LT, Jain RK. Transport of fluid and macromolecules in tumors. II. Role of heterogeneous perfusion and lymphatics. Microvasc Res. 1990; 40: 246–263. pmid:2250603
12. Markl M. Velocity encoding and flow imaging, Markl-FlowImaging.pdf. 2005; http://ee-classes.usc.edu/ee591/library/.
13. Walker-Samuel S, Ramasawmy R, Wells J, Siow B, Johnson SP, Pedley RB, et al. Investigating tumor interstitial convection currents using extra-vascular convection (EVAC) MRI. Proceedings of the British Chapter of the International Society for Magnetic Resonance in Medicine. Manchester. 2011.
14. Boucher Y, Brekken C, Netti PA, Baxter LT, Jain RK. Intratumoral infusion of fluid: estimation of hydraulic conductivity and implications for the delivery of therapeutic agents, Br J Cancer. 1998; 78: 1442–1448. pmid:9836476
15. Kurbel S, Flam J. Interstitial hydrostatic pressure: a manual for students. Adv Physiol Ed. 2007; 31: 116–117.
16. Soltani M, Chen P. Numerical modeling of fluid flow in solid tumors. PLoS ONE. 2011; 6: e20344. pmid:21673952
17. Sarntinoranont M, Rooney F, Ferrari M. Interstitial stress and fluid pressure within a growing tumor. Ann Biomed Eng. 2003; 31: 327–335. pmid:12680730
18. Jain RK, Tong RT, Munn LL. Effect of vascular normalization by antiangiogenic therapy on interstitial hypertension, peri edema, and lymphatic metastasis: insights from a mathematical model. Cancer Res. 2007; 67: 2729–2935. pmid:17363594
19. Phipps C, Kohandel M. Mathematical model of the effect of interstitial fluid pressure on angiogenic behavior in solid tumors. Comput Math Methods Med. 2011;
20. Chary SR, Jain RK. Direct measurement of interstitial convection and diffusion of albumin in normal and neoplastic tissue by fluorescence photoleaching, Proc Natl Acad Sci USA. 1989; 86: 5385–5389. pmid:2748592
21. Jain RK. Transport of molecules in the tumor interstitium: A Review, Cancer Res. 1987; 47: 3039–3051. pmid:3555767
22. Levick JR. Flow through interstitium and other fibrous matrices. Q J Exp Physiol. 1987; 72: 409–438. pmid:3321140
23. Swabb EA, Wei J, Gullino PM. Diffusion and convection in normal and neoplastic tissues. Cancer Res. 1974; 34: 2814–2822. pmid:4369924
24. Bert JL, Fatt I. Relation of water transport to water content in swelling biological membranes. In: Blank M, editors. Surface Chemistry of Biological Systems. New York: Plenum Publishing Corp; 1970. pp. 287–294.
25. Ali H, Wang WB, Zevallos M, Alfano RR. Near infrared spectroscopy and imaging to probe differences in water content in normal and cancer human prostate tissues. Technol Cancer Res Treat. 2004; 3: 491–497. pmid:15453814
26. Stohrer M, Boucher Y, Stangassinger M, Jain RK. Oncotic pressure in solid tumors is elevated. Cancer Res. 2000; 60: 4251–4225. pmid:10945638
27. Boron WF. Medical Physiology; A cellular and molecular approach. 1st ed. Philadelphia: Saunders; 2003.
28. Starling EH. On the absorption of fluids from the connective tissue spaces. J Physiol. 1986; 19: 312–326.
29. Curry FE. Mechanics and thermodynamics of transcapillary exchange. In: Renken EM, Michel CC, editors. Handbook of Physiology, Section 2: The Cardiovascular System. Bethesda, MD: Amer Physiol Soc; 1984. pp. 309–374.
30. Iliff JJ, Wang M, Liao Y, Plogg BA, Peng W, Gundersen GA, et al. A paravascular pathway facilitates CSF flow through the brain parenchyma and the clearance of interstitial solutes, including amyloid β. Sci Trans Med 2012; 4:
31. Levick JR, Michel CC. Microvascular fluid exchange and the revised Starling principle. Cardiovasc Res. 2010; 87: 198–210. pmid:20200043
32. Sherwood L. Human physiology: from cells to systems. 8th ed. Boston: Cengage Learning; 2012.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
© 2016 Liu et al. This is an open access article distributed under the terms of the Creative Commons Attribution License: http://creativecommons.org/licenses/by/4.0/ (the “License”), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.
Abstract
Tumor interstitial fluid pressure (TIFP), is a physiological parameter with demonstrated predictive value for a tumor’s aggressiveness, drug delivery, as well as response to treatments such as radiotherapy and chemotherapy. Despite its utility, measurement of TIFP has been limited by the need for invasive procedures. In this work, the theoretical basis for approaching the absolute value of TIFP and the experimental method for noninvasively measuring TIFP are presented. Given specific boundary and continuity conditions, we convert theoretical variables into measurable variables by applying MRI technology. The work shows that TIFP in the central region of the tumor can be estimated by an analysis of the variation of tissue fluid motion in the tumor rim and surrounding tissue. It is determined from three noninvasive measurable parameters: i) an estimate of the velocity of the tumor interstitial fluid at the tumor surface, which is maximal, ii) a measurement of the distance from the tumor surface to where the tumor exudates are absorbed (or normalized), and iii) an estimate of the hydraulic conductivity of the interstitium through which the tumor exudate travels. We experimentally show that the fluid flow within the tumor rim is not uniform, even for a round shaped tumor, and demonstrate the procedures for the noninvasive measurement of TIFP.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer