| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Cell, Tumor, and Stem Cell Biology |
Departments of 1 Biological Regulation and 2 Veterinary Resources, Weizmann Institute of Science, Rehovot, Israel
Requests for reprints: Hadassa Degani, Department of Biological Regulation, Weizmann Institute of Science, Rehovot, Israel 76100. Phone: 972-8-9342017; Fax 972-8-936154; E-mail: hadassa.degnani{at}weizmann.ac.il.
| Abstract |
|---|
|
|
|---|
28 mm Hg and, for comparison, in orthotopic MCF7 human breast tumors which exhibited a lower IFP of
14 mm Hg, both implanted in nude mice. The MRI protocol consisted of slow infusion of the contrast agent [gadolinium-diethylenetriaminepentaacetic acid (GdDTPA)] into the blood for
2 hours, sequential acquisition of images before and during the infusion, and measurements of T1 relaxation rates before infusion and after blood and tumor GdDTPA concentration reached a steady state. Image analysis yielded parametric images of steady-state tissue GdDTPA concentration with high values of this concentration outside the tumor boundaries,
1 mmol/L, declining in the tumor periphery to
0.5 mmol/L, and then steeply decreasing to low or null values. The distribution of steady-state tissue GdDTPA concentration reflected the distribution of IFP, showing an increase from the rim inward, with a high IFP plateau inside the tumor. The changes outside the borders of the tumors with high IFP were indicative of convective transport through the interstitium. This work presents a noninvasive method for assessing the spatial distribution of tumor IFP and mapping barriers to drug delivery and transport. (Cancer Res 2006; 66(8): 4159-66) | Introduction |
|---|
|
|
|---|
Thus far, assessment of tumor IFP has been done by methods such as the wick-in-needle and the micropancture technique (12, 14). These methods are invasive and limited to measurements in few locations. It would therefore be highly useful to develop a noninvasive imaging method that would map the spatial distribution of IFP and the net transfer into the tumor interstitium. An early attempt to use magnetic resonance imaging (MRI) for mapping IFP was based on correlating proton relaxation rates with IFP values obtained by the wick-in-needle technique (17). However, the results showed that both T1 and T2 relaxation rates did not correlate with the measured IFP (17). Subsequent dynamic contrast-enhanced MRI studies, using a bolus injection of gadolinium-diethylenetriaminepentaacetic acid (GdDTPA), indicated the presence of disparities between the influx and outflux transcapillary transfer constants in breast tumors (18). This disparity with the outflux exceeding the influx constant increased from the tumor rim to the center. This distribution is in accord with the profile of IFP distribution in tumors as was previously shown (13, 19). Recently, dynamic contrast-enhanced MRI studies of fibrosarcoma mouse model showed that a decline in IFP, induced by thalidomide, was accompanied by increased plasma volume fraction and fractional efflux rate from the interstitial space to the plasma (20).
Here we present an alternative contrast-enhanced MRI method that reveals the distribution of the contrast material due to the net effect of extravasation, diffusion, and convection in ectopic NCI-H460 nonsmall-cell lung cancer tumors implanted in immunodeficient mice, which exhibit high IFP values (
28 mm Hg). For comparison, we also applied this method to investigate orthotopic MCF7 human breast tumors which exhibit a significantly lower IFP (
14 mm Hg). The contrast agent was continuously administered by slow infusion into the blood circulation, raising its blood level to a steady-state concentration. The MRI recordings monitored T1 relaxation rates and signal intensity before the start of the infusion and during the infusion, including at blood and tumor steady-state concentrations. Analysis of the changes in T1 relaxation rates yielded steady-state tissue GdDTPA concentration (mmol/tissue volume) maps of the tumors and their surrounding. These maps reflected inhibition of transfer due to elevated tumor IFP and transfer by convection in the tumor surrounding.
| Materials and Methods |
|---|
|
|
|---|
During the experiments, mice were anesthetized by inhalation of 1% isoflurane (Medeva Pharmaceuticals, Inc., Rochester, NY) in an O2/N2O (3:7) mixture applied through a nose cone. All the protocols were approved by the Weizmann Institute Animal Care and Use Committee.
Measurements of IFP. IFP was measured in H460 tumors (n = 7) 13 days after their implantation and in MCF7 tumors (n = 9)
5 weeks after their implantation using the wick-in-needle apparatus (23). Briefly, a 23-gauge needle with a side hole located at
3 mm from the needle tip was connected to a pressure monitor system (model 295-1 Pressure, Stryker, Kalamazoo, MI) especially designed for measuring tissue fluid pressures. The system was filled with saline. The needle was inserted into a central part of the tumor or into the flank muscle (n = 20) for reference, and 50 µL of 0.9% sodium chloride were injected to ensure fluid communication between the tissue and the pressure monitor system.
Histology. At the end of the MRI experiments, tumors were dissected free from s.c. tissue and cut in the middle in a plane parallel to that of the magnetic resonance images. Tumors were fixed in 2.5% paraformaldehyde for H&E staining or in 4% zinc solution (0.5% zinc chloride and zinc acetate, 0.05% calcium acetate in 0.1 mol/L tris buffer, pH 7.4) for CD31 immunofluorescence staining. Seven-micrometer-thick paraffin-embedded sections were prepared.
For endothelial staining, the sections were deparaffinized in xylene, hydrated in series of graded ethanol, and rinsed in double-distilled water. Sections immersed in cold acetone (20°C, 7 minutes) and rinsed in double-distilled water were then blocked with 20% normal rabbit serum and incubated overnight with primary rat anti-mouse CD31 (platelet/endothelial cell adhesion molecule 1, monoclonal rat anti-mouse, 1:100; PharMingen, San Diego, CA), which is constitutively expressed on the surface of mature endothelial cells. Next, sections were incubated for 1 hour in biotinylated antirat CD31 (1:100; Vector Lab, Burlingame, CA). Antibody distribution was visualized using a fluorescence streptavidin-Cy3 conjugated complex (Jackson ImmunoResearch Laboratories, Inc., Baltimore, MD). Nuclear staining was done using Hoechst solution (Molecular Probes, Eugene, OR) diluted 1:2,000.
Stained sections were examined by fluorescence microscope (E600, Nikon, Toyo, Japan) equipped with Plan Fluor objective connected to a CCD camera (DMX1200F, Nikon). Digital images (3.2 mm2) of all the tumor area were collected and analyzed using the Image-Pro plus 4.1 software. Quantification of endothelial staining was done by measuring the percentage of area with positive staining for CD31.
MRI studies. MRI scans were acquired with a 4.7-T Biospec spectrometer (Bruker Biospin, Rheinstetten, Germany). Fourteen H460 tumors and nine MCF7 tumors were scanned using a protocol that included an initial two-dimensional T2-weighted spin echo sequence with echo time = 68 ms; repetition time = 2,500 ms; 128 x 128 matrix; 1-mm slice thickness; an interslice distance of 1.1; and 3 x 3 cm2 field of view. The tumor region of interest in each slice was traced on the T2-weighted images and this trace was then used for localizing the tumor in the various subsequent images obtained at the same spatial resolution. The size of the tumors was determined from the area of the region of interest and the slice thickness, taking into account the inter-slice distance (24). T1 measurements were then done using two-dimensional sequential inversion recovery snapshot fast low-angle shot imaging with 11 inversion times ranging from 10 ms to 10 seconds; echo time = 3.5 ms; repetition time = 15 ms; flip angle = 10 degrees; and the same matrix size and field of view as the T2-weighted images. Two-dimensional T1-weighted gradient echo images were also scanned using echo time/repetition time = 2.73/35.8 ms; flip angle = 60 degrees; and the same spatial resolution as the T2-weighted images acquiring four scans within 18 seconds. Following these measurements, slow infusion was initiated with 0.05 mol/L GdDTPA solution (gadopentate-dimeglumine, Schering, Berlin, Germany) at a rate of 0.66 mmol/h/kg wt for 2 hours. Sequential images were scanned during the slow infusion using the three-dimensional T1-weighted gradient echo sequence described above. At 90 minutes after the start of infusion, T1 relaxation rates were measured again using the inversion recovery snapshot fast low-angle shot sequence described above.
Separate experiments were done to monitor the GdDTPA enhancement in the carotid arteries and determine the time needed to reach steady-state in the blood during the slow infusion (n = 3). In these experiments, inversion recovery fast low-angle shot sequence was applied using a fixed inversion time of 120 ms; echo time = 3.5 ms; repetition time = 15 ms; flip angle = 10 degrees; 128 x 128 matrix; 1-mm slice thickness; and 3 x 3 cm2 field of view at a temporal resolution of 6 s.
Processing and analysis. T1 relaxation rates were calculated at pixel resolution applying a nonlinear least square fit of the intensity I per pixel at varying inversion times. The curves obtained from measurements before administration of the contrast agent and at steady-state infusion conditions exhibited a single decay time constant according to the following equation: I = Iinf [1 A exp(TI / T1)], with Iinf [maximum intensity at a long inversion time (TI)], A (maximum value 2), and T1 as the free variables in this fit with R2 of the fit ranging from 0.9 to
1. At steady-state concentration of the contrast agent in the tissues, the intracellular and extracellular T1 water relaxation rates differ due to the sole presence of GdDTPA in the extracellular compartment. At a maximum concentration in the extracellular compartment of
1 mmol/L, the upper limit for this difference is
4 s1 whereas the effective intracellular-extracellular water exchange rate is more than an order of magnitude higher, 50 s1 (based on intracellular lifetime of
100 ms and an intracellular to extracellular volume ratio of 4). Hence, water exchange between the intracellular and extracellular compartments is at the fast exchange limit (25) and the T1 relaxation rate at steady state is decaying uniexponentially as was indeed found in the T1 measurements. Furthermore, under this fast exchange condition, tissue GdDTPA concentration (CGd), defined as the amount of GdDTPA in millimoles per tissue volume at GdDTPA steady-state concentration, is obtained from the measured relaxation rates according to the equation CGd = (1/T1ss 1/T10) / r1, where r1 is the water relaxivity of GdDTPA in solution, 4.2 s1 x (mmol/L)1 (26), and 1/T1ss and 1/T10 are the relaxation rates at steady-state concentration and before the infusion, respectively.
Maps of the actual GdDTPA concentration in the extracellular volume fraction of H460 tumors were calculated estimating an extracellular volume fraction of 0.2 (14). The extracellular volume fraction of MCF7 tumors was determined by applying a method based on diffusion MRI (22). Further estimation of IFP in the H460 tumors was obtained by assuming a linear relation between IFP and the calculated GdDTPA concentration in the extracellular volume at steady state using an approximate scale of IFP between 0 mm Hg at the rim and 28 mm Hg at the center.
An attempt was made to analyze the enhancement curves during the first 30 minutes of the infusion using the kinetic model described by Tofts and Berkowitz (27) and a nonlinear least square fit program previously developed in our laboratory (28). The output of this analysis yielded the influx and efflux transcapillary transfer rate constants.
| Results |
|---|
|
|
|---|
The IFP of the tumors was determined by the wick-in-needle method. Attempts were made to measure the pressure close to the center of the tumors. H460 tumors exhibited high IFP values, ranging from 18 to 45 mm Hg with a mean ± SD of 28 ± 8 mm Hg (n = 7). The IFP values of MCF7 tumors were lower and more diverse, ranging from 4 to 32 mm Hg with a mean of 14 ± 10 mm Hg (n = 9). Measurement in the flank muscle opposite to the tumor site and of control mice showed IFP values ranging from 0 to 5 mm Hg.
The analysis of histologic sections of H460 tumors revealed their structural and morphologic features. Sections stained with H&E showed well-delineated, but nonencapsulated, tumors composed of cohesive, densely cellular, and disorganized sheets of pleomorphic cells, exhibiting occasional invasion into the surrounding connective tissue. The tumors predominantly consisted of viable, cellular regions with occasionally small localized regions of necrosis spread homogeneously over the entire tumor section (Fig. 1A and B ). Specific immunostaining of endothelial cells with CD31 monoclonal antibody traced the capillary distribution in the tumors (Fig. 1C-F). The density of the capillaries varied both at the boundaries and in central regions; part of the regions exhibited high density (Fig. 1C and F) and others had low density (Fig. 1D and E). Thus, the heterogeneous distribution of the capillaries was similar in the centers and rims of the tumors with an average CD31 stained area of 1.1 ± 0.4% [7 tumors, in each 4 to 10 fields of 3.2 mm2 (2 x 1.6 mm)]. The histologic features of the orthotopic MCF7 tumors were similar to those previously described (22) and revealed the presence of a larger extracellular volume fraction and a more heterogeneous distribution of the cells than in the H460 tumors. Staining of the capillaries indicated an average stained area similar to that of H460 tumors of 1.3 ± 0.6%.
|
20 minutes after the start of the infusion, the concentration of the contrast agent in the blood circulation reached a steady state (of
0.9 mmol/L); i.e., the rate of infusion was equal to the rate of clearance from the blood through the kidneys into the urine. Following blood steady state, other parts in the mouse body also reached steady-state concentration (Fig. 2B and C).
|
|
70 minutes after blood steady-state concentration was reached; Figs. 4A-C
, 5A-C
, and 6A-C
). The high tissue GdDTPA concentration outside the tumors (in the range 1-2 mmol/L) reflected transport via convection. The steep decrease from the tumor rim towards the center reflected the increase in IFP. Inner regions of part of the H460 tumors were completely void of GdDTPA despite the long duration of infusion, indicating the presence of high IFP and convective flow outwards (Figs. 4C and 5C). This concentration distribution was in contrast to the homogeneous morphology of the tumors as was indicated by histologic sections in planes similar to those of the images (Fig. 1A). It was also in contrast to the distribution of capillaries that did not show large differences in the microvascular density between the rim and the center (Fig. 1).
|
|
|
Assuming a linear relation between GdDTPA concentration at steady state and IFP, and using an average scaling of the range of IFP from 0 to 28 mm Hg (based on the wick-in-needle measurements), yielded estimated IFP maps and profiles as shown in Fig. 7 for three H460 tumors. The linear relation was based on the assumption that the flux of GdDTPA is predominantly determined by the extravasation across the capillary walls and that the effect of transport in the interstitium is negligible.
|
1.0, which, according to the model, estimate the extracellular volume fraction. In contrast, analysis of the dynamics in MCF7 tumors with a lower IFP yielded kin values of (1.1 ± 0.2) x 102 min1 and kep values of (3.0 ± 0.9) x 102 min1, resulting in a median kin/kep of 0.4 ± 0.1, which was similar to the extracellular volume fraction measured by the MRI diffusion method. We therefore conclude that it is important to include in the model of dynamic contrast enhancement terms that will account for pressure gradients, particularly in tumors with high IFP. | Discussion |
|---|
|
|
|---|
28 mm Hg and are therefore suitable as a tumor model in this study. For comparison, we also investigated orthotopic MCF7 human breast tumors that exhibited a lower IFP of
14 mm Hg and higher extracellular volume fraction but a similar distribution of blood vessels. In general, reported IFP values in tumors range between 10 and 50 mm Hg whereas IFP values in normal tissues, including the muscle tissue of nude mice, range between 2 and 0 mm Hg (2931).
In most previous dynamic contrast-enhanced MRI studies of tumors in rodents and humans, the contrast agent was administered by a bolus injection. Our initial attempts to monitor signal enhancement in H460 tumors after a bolus administration of the contrast agent showed that only a small fraction of the tumor pixels (7-25%) were enhanced (data not shown). These results suggested slow transcapillary transfer rates of GdDTPA through the microvascular network of the tumors. To improve detection of enhancement in regions with slow transfer rates, we applied a slow and long infusion protocol. Moreover, extending the infusion time to allow diffusion of GdDTPA in the interstitium (calculated tissue diffusion distance of GdDTPA is
2.4 mm during 90 minutes; refs. 18, 32) and reaching tumor steady-state conditions made it possible to accurately determine T1 relaxation times at this stage. This further enabled us to map the tumor steady-state tissue GdDTPA concentration and estimate the IFP distribution.
In the absence of IFP, GdDTPA transfer from the plasma into the interstitial space and back is due to diffusion through the capillary walls and in the interstitial space in the direction of the concentration gradients. At blood and tumor steady state, the extracellular concentration throughout the tumor should be equal to the plasma concentration. Tissue GdDTPA concentration distribution under this condition would therefore be directly proportional to the extracellular volume fraction. Typically, the distribution of the steady-state tumor tissue GdDTPA concentration showed a steep decline from the rim towards the tumor center. Such a distribution was not found in the histologic sections for the cell density (and hence extracellular volume fraction), as well as for the capillary density. The observed steep decline in the tissue GdDTPA concentration, therefore, reflected an increase in IFP with regions void of GdDTPA exhibiting highest IFP.
By using an approximate homogeneous interstitial volume fraction and scaling IFP according to the wick-in-needle results, we obtained an estimation of IFP distribution in the H460 tumors (Fig. 7). We cannot exclude small local changes in the extracellular volume fraction which could affect the calculated extracellular GdDTPA concentrations and, hence, the calculated IFP values; however, these changes will not markedly modify the general profile across the tumor. This distribution is highly comparable to the IFP distribution measured previously by the wick-in-needle and micropuncture techniques and supports this approximation (13, 19, 33). Furthermore, the model of Baish et al. (34), which described the steady-state coupling between vascular and interstitial flows, predicted uniformly elevated IFP in the central region of tumors and a rapid decline to normal tissue values of 0 mm Hg at the periphery. However, independent measurements of the spatial variations in the extracellular volume fraction and comparison with invasive regional IFP measurements are necessary to verify and standardize this slow infusion MRI method.
We further assumed a linear relation between GdDTPA concentration in the extracellular volume fraction and the IFP. In general, the spatial distribution of GdDTPA concentration in the extracellular volume fraction is a result of the transcapillary extravasation flux and the interstitial transport flux (14, 35). The GdDTPA extravasation flux depends on the diffusion and the hydrostatic water pressure difference between the capillary pressure and IFP [e.g., see Eq. (F) in ref. 35]. The interstitial transport flux depends on the interstitial concentration gradient (diffusion component) and pressure gradient across the tumor interstitium convection component [e.g., see Eq (G) in ref. 35]. We have assumed that in H460 tumors, with viable cells and blood vessels distributed throughout the whole tumor, the extravasation flux is dominant; therefore, under the steady-state condition, when GdDTPA flux is zero, the concentration gradient is linearly related to IFP. However, we cannot exclude the presence of an interstitial convective term which yields at steady state a logarithmic relation between the extracellular volume fraction concentration and IFP, and hence diverts the relation from linearity.
Efforts have previously been made to generally describe fluid flow and tracer transport in solid tumors. Specifically, the pathophysiologic effects of elevated IFP in tumors have been explored by using mathematical models developed by Jain and his group (34, 3638). Their model predicted a distribution of IFP with low IFP outside the tumor in the boundary with the normal tissue, which increases as the distance from the tumor edge towards the center increases depending on the ratio of the vascular to interstitial hydraulic conductivities (35, 36). Thus, for high values of this ratio, a plateau of high IFP is established in a large fraction of the inner tumor volume. Milosevic et al. (29) extended the mathematical modeling to investigate the relationship between elevated IFP and blood flow. Their model predicted elevated tumor IFP in the range of 0 to 56 mm Hg. Recently, a mathematical model was developed which included the microstructure of the vasculature allowing variations of vascular architecture, size, and conductivity (39). This model captured the strong coupling of the capillary and interstitial flow occurring for highly permeable vascular walls and showed that the assumption of uniform IFP is not generally appropriate. Beard and Bassingthwaighte (40) also extended the modeling to include convection and diffusion in three-dimensional complex physiologic geometries. Although these models are very useful to make new predictions, it is necessary to have new experimental methods that will provide data supporting or dismissing the theoretical predictions. Our method opens the way to obtain data for evaluating these theoretical predictions.
We examined here the use of a physiologic two-compartment exchange model that was adapted for a slow infusion protocol by Tofts (41). The mechanisms of transfer did not include the additional contribution of IFP, and the basic assumptions made in formulating this model were the same as those made for the bolus administration of the contrast agent (27, 41). In the MCF7 tumors that exhibited a low IFP, we obtained physiologic reasonable influx and efflux transfer constants as well as extracellular volume fractions. However, neglecting IFP and fluxes due to pressure gradients may have influenced the calculated values of the transcapillary transfer constants, particularly in the central regions. Analysis of the dynamic data of the H460 tumors that exhibited high IFP clearly showed the weakness of this model yielding underestimated efflux transfer constants and unrealistic extracellular volume fraction. This emphasizes the need to extend the physiologic model of dynamic contrast-enhanced MRI and explicitly add the influence of pressure gradients.
In summary, although transcapillary transfer and transport of low molecular weight substances are believed to be driven mainly by passive diffusion, our results indicate that transfer of the common gadolinium-based MRI contrast agents can be pressure-driven in the presence of elevated IFP. We also showed that the distribution of the most common MRI contrast agent under steady-state conditions can serve to map IFP and predict the presence of barriers to drug delivery.
| Acknowledgments |
|---|
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
We thank the Lord David Alliance, CBE, UK and the Estate of Julie Osler, USA.
| Footnotes |
|---|
Received 9/14/05. Revised 1/ 1/06. Accepted 1/31/06.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
T. P.F. Gade, I. M. Buchanan, M. W. Motley, Y. Mazaheri, W. M. Spees, and J. A. Koutcher Imaging Intratumoral Convection: Pressure-Dependent Enhancement in Chemotherapeutic Delivery to Solid Tumors Clin. Cancer Res., January 1, 2009; 15(1): 247 - 255. [Abstract] [Full Text] [PDF] |
||||
![]() |
F. Pastorino, D. Di Paolo, F. Piccardi, B. Nico, D. Ribatti, A. Daga, G. Baio, C. E. Neumaier, C. Brignole, M. Loi, et al. Enhanced Antitumor Efficacy of Clinical-Grade Vasculature-Targeted Liposomal Doxorubicin Clin. Cancer Res., November 15, 2008; 14(22): 7320 - 7329. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Cancer Research | Clinical Cancer Research |
| Cancer Epidemiology Biomarkers & Prevention | Molecular Cancer Therapeutics |
| Molecular Cancer Research | Cancer Prevention Research |
| Cancer Prevention Journals Portal | Cancer Reviews Online |
| Annual Meeting Education Book | Meeting Abstracts Online |