Transvascular exchange of fluid and plasma proteins

A theoretical model of transvascular exchange of fluid and plasma proteins in the microcirculation is developed based on fundamental laws of the fluid mechanics and on phenomenological transport equations of the irreversible thermodynamics. Intravascular axial changes of the pressure, flow and plasma protein concentration are taken into account as well as axial gradients of vascular permeability. Proper nondimensionalization of the resulting equations leads to the identification of dimensionless parameters which combine the transport characteristics of the endothelial wall and the intravascular flow resistance. In the theory, the dependence of the reflection coefficient on the transport coefficients of the vascular wall and on the plasma protein concentration is established. The model is applied to the cat mesentery and the rat intestinal muscle. The numerical simulations indicate that taking into account vascular protein permeability yields considerable differences in the axial distribution of the plasma protein concentration and transvascular fluxes in comparison with the case of protein impermeability of the endothelial wall. The results show that the maximum of the transvascular fluid and plasma protein movement resides at the site of the small venules while a minimum of the exchange occurs at the site of the midcapillaries.     

In vivo mimicking model for solid tumor towards hydromechanics of tissue deformation and creation of necrosis

The present work addresses transvascular and interstitial fluid transport inside a solid tumor surrounded by normal tissue (close to an in vivo mimicking setup). In general, biological tissues behave like a soft porous material and show mechanical behavior towards the fluid motion through the interstitial space. In general, forces like viscous drag that are associated with the fluid flow may compress the tissue material. On the macroscopic level, we try to model the motion of fluids and macromolecules through the interstitial space of solid tumor and the normal tissue layer. The transvascular fluid transport is assumed to be governed by modified Starling’s law. The poroelastohydrodynamics (interstitial hydrodynamics and the deformation of tissue material) inside the tumor and normal tissue regions is modeled using linearized biphasic mixture theory. Correspondingly, the velocity distribution of fluid is coupled to the displacement field of the solid phase (mainly cellular phase and extracellular matrix) in both the normal and tumor tissue regions. The corresponding velocity field is used within the transport reaction equation for fluids and macromolecules through interstitial space to get the overall solute (e.g., nutrients, drug, and other macromolecules) distribution. This study justifies that the presence of the normal tissue layer plays a significant role in delaying/assisting necrosis inside the tumor tissue. It is observed that the exchange process of fluids and macromolecules across the interface of the tumor and normal tissue affects the effectiveness factor corresponding to the tumor tissue.     

A model of unsteady-state transvascular fluid and protein transport in the lung

Models of steady-state fluid and solute transport in the microcirculation are used primarily to characterize filtration and permeability properties of the transport barrier. Important transient relationships, such as the rate of fluid accumulation in the tissue, cannot be predicted with steady-state models. In this paper we present three simple models of unsteady-state fluid and protein exchange between blood plasma and interstitial fluid. The first treats the interstitium as a homogeneous well-mixed compliant compartment, the second includes an interstitial gel, and the third allows for both gel and free fluid in the interstitium. Because we are primarily interested in lung transvascular exchange we used the multiple-pore model and pore sizes described by Harris and Roselli (J. Appl. Physiol.: Respirat . Environ. Exercise Physiol. 50: 1-14, 1981) to characterize the microvascular barrier. However, the unsteady-state transport theory presented here should apply to other organ systems and can be used with different conceptual models of the blood-lymph barrier. For a step increase in microvascular pressure we found good agreement between theoretical and experimental lymph flow and lymph concentrations in the sheep lung when the following parameter ranges were used: base-line interstitial volume, 150-190 ml; interstitial compliance, 7-10 ml/Torr; initial interstitial fluid pressure, -1 Torr; pressure in initial lymphatics, -5 to -6 Torr; and conductivity of the interstitium and lymphatic barrier, 4.25 X 10(-4) ml X s-1 X Torr-1. Based on these values the model predicts 50% of the total change in interstitial water volume occurs in the first 45 min after a step change in microvascular pressure.(ABSTRACT TRUNCATED AT 250 WORDS)     

A mixture theory model of fluid and solute transport in the microvasculature of normal and malignant tissues. I. Theory

In order to better understand the mechanisms governing transport of drugs, nanoparticle-based treatments, and therapeutic biomolecules, and the role of the various physiological parameters, a number of mathematical models have previously been proposed. The limitations of the existing transport models indicate the need for a comprehensive model that includes transport in the vessel lumen, the vessel wall, and the interstitial space and considers the effects of the solute concentration on fluid flow. In this study, a general model to describe the transient distribution of fluid and multiple solutes at the microvascular level was developed using mixture theory. The model captures the experimentally observed dependence of the hydraulic permeability coefficient of the capillary wall on the concentration of solutes present in the capillary wall and the surrounding tissue. Additionally, the model demonstrates that transport phenomena across the capillary wall and in the interstitium are related to the solute concentration as well as the hydrostatic pressure. The model is used in a companion paper to examine fluid and solute transport for the simplified case of an axisymmetric geometry with no solid deformation or interconversion of mass.     

Interstitial hydraulic conductivity in a fibrosarcoma

Convective transport of therapeutic agents in solid tumors can be improved through intratumoral infusion. To optimize the convection, we investigated the dependence of the hydraulic conductivity on tissue deformation induced by interstitial fluid pressure gradient during the infusion. Two experimental systems were used in the investigation: 1) one-dimensional perfusion through tumor slices and 2) intratumoral infusion using a needle. With these systems, we found that the apparent hydraulic conductivity (K(app)) could be altered by several orders of magnitude in fibrosarcomas through changes in perfusion conditions. When the perfusion pressure was less than a threshold level, fluid flow in tissues could not be detected. When the perfusion pressure was increased above the threshold level, K(app) depended on perfusion system and pressure. The maximum variation in K(app) in fibrosarcomas reached 80,260-fold in our experiments. The large variation in K(app) could be explained by perfusion pressure-induced tissue deformation. These experimental data suggest that the hydraulic conductivity is very sensitive to tissue deformation and imply that it is possible to improve intratumoral infusion of therapeutic agents through optimization of infusion conditions.     

Effect of complement depletion on lung fluid balance after thrombin

We examined the effect of complement depletion on lung fluid and protein exchange after thrombin-induced pulmonary thromboembolization. Sheep were prepared with lung lymph fistulas to assess pulmonary transvascular fluid and protein dynamics. Studies were made in three groups: in group I (n = 5) pulmonary thromboembolization (PT) was induced by an iv infusion of thrombin (55.0 +/- 12.9 NIH U/kg); in group II (n = 6) cobra venom factor (CVF) was given ip (94.5 +/- 18.8 U/kg/day) for 2 days to deplete complement, and then thrombin (66.4 +/- 37.0 NIH U/kg) was infused to raise pulmonary vascular resistance to the same level as in group I; in group III (n = 10) left atrial pressure (Pla) was increased by 10-15 Torr in normal animals by inflation of a Foley balloon catheter. In group I, thrombin infusion caused an increase in pulmonary lymph flow (Qlym) with a gradual increase in the lymph-to-plasma protein concentration ratio (L/P). In complement-depleted sheep, thrombin caused a transient increase in Qlym, which was associated with a decrease in L/P. In group I an increase in Pla further increased Qlym but without a change in L/P, indicating an increase in lung vascular permeability to proteins; whereas in the decomplemented-thrombin sheep raising Pla increased Qlym but decreased L/P. Results in the latter group were similar to those obtained in normal animals after left atrial hypertension (group III). Therefore the complement system participates in the increase in lung vascular permeability following thrombin-induced microembolization.     

Thrombin-induced alterations in lung fluid balance in awake sheep

We examined the effect of fibrinolysis depression on thrombin-induced pulmonary microembolism in awake sheep prepared with chronic lung lymph fistulas. Fibrinolysis was depressed by an intravenous infusion (100 mg) of tranexamic acid [trans-4-(Aminomethyl)cyclohexanecarboxylic acid]. Pulmonary microembolism was induced by an intravenous infusion of alpha-thrombin (80 NIH U/kg) in normal (n = 7) and in tranexamic acid-treated (n = 6) sheep. Thrombin immediately increased pulmonary lymph flow (Qlym) in both groups. The increased Qlym was not associated with a change in the lymph-to-plasma protein concentration (L/P) ratio in the control group and with a small decrease in the tranexamic acid-treated group. The increases in Qlym and pulmonary transvascular protein clearance (Qlym X L/P ratio) in the tranexamic acid-treated group were greater and sustained at four- to fivefold above base line for 10 h after the thrombin and remained elevated at twofold above base line even at 24 h. In contrast, Qlym and protein clearance were transiently increased in the control group. The mean pulmonary arterial pressure (Ppa) and pulmonary vascular resistance (PVR) increased after thrombin in tranexamic acid-treated group; the increases in Ppa and PVR in the control group were transient. Protein reflection coefficient as determined by the filtration independent method decreased after thrombin in tranexamic acid-treated sheep (n = 5), indicating an increased vascular permeability to proteins. We conclude that prolongation of microthrombi retention in the pulmonary circulation results in an increased vascular permeability to proteins. Both increased vascular permeability and vascular hydrostatic pressure are important determinants of the increases in Qlym and transvascular protein clearance after thrombin-induced pulmonary microembolism.     

Study of tumor blood perfusion and its variation due to vascular normalization by anti-angiogenic therapy based on 3D angiogenic microvasculature

The coupling of intravascular and interstitial flow is a distinct feature of tumor microcirculation, due to high vessel permeability, low osmotic pressure gradient and absence of functional lymphatic system inside tumors. We have previously studied the tumor microcirculation by using a 2D coupled model. In this paper, we extend it to a 3D case with some new considerations, to investigate tumor blood perfusion on a more realist microvasculature, and the effects of vascular normalization by anti-angiogenic therapy on tumor microenvironment.The model predict the abnormal tumor microcirculation and the resultant hostile microenvironment: (1) in the intra-tumoral vessels, blood flows slowly with almost constant pressure values, haematocrit is much lower which contributes to hypoxia and necrosis formation of the tumor centre; (2) the total transvascular flux is at the same order of magnitude as intravascular flux, the intravasation appears inside of the tumor, the ratio of the total amount of intravasation flux to extravasation flux is about 16% for the present model; (3) the interstitial pressure is uniformly high throughout the tumor and drops precipitously at the periphery, which leads to an extremely slow interstitial flow inside the tumor, and a rapidly rising convective flow oozing out from the tumor margin into the surrounding normal tissue. The investigation of the sensitivity of flows to changes in transport properties of vessels and interstitium as well as the vascular density of the vasculature, gains an insight into how normalization of tumor microenvironment by anti-angiogenic therapies influences the blood perfusion.     

Effects of inspiratory resistance loading on lung fluid balance in awake sheep

Because pulmonary edema has been associated clinically with airway obstruction, we sought to determine whether decreased intrathoracic pressure, created by selective inspiratory obstruction, would affect lung fluid balance. We reasoned that if decreased intrathoracic pressure caused an increase in the transvascular hydrostatic pressure gradient, then lung lymph flow would increase and the lymph-to-plasma protein concentration ratio (L/P) would decrease. We performed experiments in six awake sheep with chronic lung lymph cannulas. After a base-line period, we added an inspiratory load (20 cmH2O) and allowed normal expiration at atmospheric pressure. Inspiratory loading was associated with a 12-cmH2O decrease in mean central airway pressure. Mean left atrial pressure fell 11 cmH2O, and mean pulmonary arterial pressure was unchanged; calculated microvascular pressure decreased 8 cmH2O. The changes that occurred in lung lymph were characteristic of those seen after other causes of increased transvascular hydrostatic gradient, such as increased intravascular pressure. Lung lymph flow increased twice base line, and L/P decreased. We conclude that inspiratory loading is associated with an increase in the pulmonary transvascular hydrostatic gradient, possibly by causing a greater fall in interstitial perimicrovascular pressure than in microvascular pressure.     

Microperfusion Technique to Investigate Regulation of Microvessel Permeability in Rat Mesentery

Experiments to measure the permeability properties of individually perfused microvessels provide a bridge between investigation of molecular and cellular mechanisms regulating vascular permeability in cultured endothelial cell monolayers and the functional exchange properties of whole microvascular beds. A method to cannulate and perfuse venular microvessels of rat mesentery and measure the hydraulic conductivity of the microvessel wall is described. The main equipment needed includes an intravital microscope with a large modified stage that supports micromanipulators to position three different microtools: (1) a beveled glass micropipette to cannulate and perfuse the microvessel; (2) a glass micro-occluder to transiently block perfusion and enable measurement of transvascular water flow movement at a measured hydrostatic pressure, and (3) a blunt glass rod to stabilize the mesenteric tissue at the site of cannulation. The modified Landis micro-occlusion technique uses red cells suspended in the artificial perfusate as markers of transvascular fluid movement, and also enables repeated measurements of these flows as experimental conditions are changed and hydrostatic and colloid osmotic pressure difference across the microvessels are carefully controlled. Measurements of hydraulic conductivity first using a control perfusate, then after re-cannulation of the same microvessel with the test perfusates enable paired comparisons of the microvessel response under these well-controlled conditions. Attempts to extend the method to microvessels in the mesentery of mice with genetic modifications expected to modify vascular permeability were severely limited because of the absence of long straight and unbranched microvessels in the mouse mesentery, but the recent availability of the rats with similar genetic modifications using the CRISPR/Cas9 technology is expected to open new areas of investigation where the methods described herein can be applied.     

A Mathematical Model of the Enhanced Permeability and Retention Effect for Liposome Transport in Solid Tumors

The discovery of the enhanced permeability and retention (EPR) effect has resulted in the development of nanomedicines, including liposome-based formulations of drugs, as cancer therapies. The use of liposomes has resulted in substantial increases in accumulation of drugs in solid tumors; yet, significant improvements in therapeutic efficacy have yet to be achieved. Imaging of the tumor accumulation of liposomes has revealed that this poor or variable performance is in part due to heterogeneous inter-subject and intra-tumoral liposome accumulation, which occurs as a result of an abnormal transport microenvironment. A mathematical model that relates liposome accumulation to the underlying transport properties in solid tumors could provide insight into inter and intra-tumoral variations in the EPR effect. In this paper, we present a theoretical framework to describe liposome transport in solid tumors. The mathematical model is based on biophysical transport equations that describe pressure driven fluid flow across blood vessels and through the tumor interstitium. The model was validated by direct comparison with computed tomography measurements of tumor accumulation of liposomes in three preclinical tumor models. The mathematical model was fit to liposome accumulation curves producing predictions of transport parameters that reflect the tumor microenvironment. Notably, all fits had a high coefficient of determination and predictions of interstitial fluid pressure agreed with previously published independent measurements made in the same tumor type. Furthermore, it was demonstrated that the model attributed inter-subject heterogeneity in liposome accumulation to variations in peak interstitial fluid pressure. These findings highlight the relationship between transvascular and interstitial flow dynamics and variations in the EPR effect. In conclusion, we have presented a theoretical framework that predicts inter-subject and intra-tumoral variations in the EPR effect based on fundamental properties of the tumor microenvironment and forms the basis for transport modeling of liposome drug delivery.     

CVF-induced decomplementation: effect on lung transvascular protein flux after thrombin

We examined the effects of cobra venom factor (CVF) on the changes in pulmonary hemodynamics and transvascular fluid and protein exchange following thrombin-induced pulmonary microembolism. Studies were made in unanesthetized sheep prepared with lung lymph fistulas. The animals received tranexamic acid (100 mg) to suppress fibrinolysis and were then challenged with an intravenous infusion of alpha-thrombin (80 U/kg). Control-thrombin challenged sheep were compared with the CVF-treated sheep challenged with the same thrombin dosage. CVF treatment (187 U X kg-1 X day-1 for 4 days) decreased the total hemolytic complement activity by 45% of control. Thrombin infusion in control sheep increased the mean pulmonary arterial pressure (Ppa), pulmonary vascular resistance (PVR), and lymph protein clearance (pulmonary lymph flow X lymph-to-plasma protein concentration ratio, Clym). Thrombin infusion in CVF-treated sheep produced smaller increments in Ppa, PVR, and Clym. Pulmonary lymph obtained from control-thrombin and CVF-thrombin sheep induced migration of granulocytes obtained from normal unchallenged sheep. The granulocytes obtained from CVF-treated sheep responded relatively less to the migratory and O-2-generating stimuli (i.e., zymosan-treated serum, pulmonary lymph from sheep after thrombin challenge, and plasma from sheep after CVF treatment) compared with normal granulocytes. The attenuation of the thrombin-induced increases in Ppa, PVR, and lung transvascular fluid and protein exchange by CVF treatment may be the result of impaired function of granulocytes.     

A nonlinear biphasic model of flow-controlled infusions in brain: Mass transport analyses

A biphasic nonlinear mathematical model is proposed for the mass transport that occurs during constant flow-rate infusions into brain tissue. The model takes into account geometric and material nonlinearities and a hydraulic conductivity dependent upon strain. The biphasic and convective–diffusive transport equations were implemented in a custom-written code assuming spherical symmetry and using an updated Lagrangian finite element algorithm. Results of the model indicate that the inclusion of these nonlinearities produced modest changes in the interstitial concentration but important variations in drug penetration and bulk concentration. Increased penetration of the drug but smaller bulk concentrations were obtained at smaller strains caused by combination of parameters such as increased Young’s modulus and initial hydraulic conductivity. This indicates that simulations of constant flow-rate infusions under the assumption of infinitesimal deformations or rigidity of the tissue may yield lower bulk concentrations near the infusion cavity and over-predictions of the penetration of the infused agent. The analyses also showed that decrease in the infusion flow rate of a fixed amount of drug results in increased penetration of the infused agent. From the clinical point-of-view, this may promote a safer infusion that delivers the therapeutic range over the desired volume while avoiding damage to the tissue by minimizing deformation and strain.     

Evaluating pulmonary vascular permeability with radiolabeled proteins: an error analysis

Using techniques of mathematical simulation, we compared two methods of evaluating pulmonary vascular permeability, i.e., transvascular protein flux. Both methods calculate a transport rate constant [pulmonary transcapillary escape rate (PTCER)] after making external radiation detection measurements of an intravenously administered radiolabeled protein. With one method, lung tissue time-activity data are acquired by positron emission tomography (PET) and are interpreted with a two-compartment model. With the other method, the time-activity data are acquired with simple detector probes and are interpreted by linear regression after normalizing for various physical factors (slope-intercept or SI method). The results show that significant errors in calculating PTCER can result from using the SI method, because it ignores the effects of back-flux on the tissue time-activity measurements. Both methods produce errors if the data analysis includes activity from vascular volumes not involved in tracer exchange with the extravascular compartment. Significant errors can also occur with the PET method, particularly when permeability is nearly normal, if pulmonary vascular volume changes significantly during the period of data collection. On balance, the PET method appears to be the method of choice for accurately evaluating pulmonary vascular permeability by protein flux measurements, although both methods may be useful in clinical applications.     

Pulmonary microcirculatory responses to leukotrienes B4, C4 and D4 in sheep

The pulmonary microvascular responses to leukotrienes B4, C4, and D4 (total dosage of 4 micrograms/kg i.v.) were examined in acutely-prepared halothane anesthetized and awake sheep prepared with lung lymph fistulas. In anesthetized as well as unanesthetized sheep, LTB4 caused a marked and transient decrease in the circulating leukocyte count. Pulmonary transvascular protein clearance (pulmonary lymph flow X lymph-to-plasma protein concentration ratio) increased transiently in awake sheep, suggesting a small increase in pulmonary vascular permeability. The mean pulmonary artery pressure (Ppa) also increased. In the acutely-prepared sheep, the LTB4-induced pulmonary hemodynamic and lymph flow responses were damped. Leukotriene C4 increased Ppa to a greater extent in awake sheep than in anesthetized sheep, but did not significantly affect the pulmonary lymph flow rate (Qlym) and lymph-to-plasma protein concentration (L/P) ratio in either group. LTD4 increased Ppa and Qlym in both acute and awake sheep; Qlym increased without a significant change in the L/P ratio. The LTD4-induced rise in Ppa occurred in association with an increase in plasma thromboxane B2 (TxB2) concentration. The relatively small increase in Qlym with LTD4 suggests that the increase in the transvascular fluid filtration rate is the result of a rise in the pulmonary capillary hydrostatic pressure. In conclusion, LTB4 induces a marked neutropenia, pulmonary hypertension, and may transiently increase lung vascular permeability. Both LTC4 and LTD4 cause a similar degree of pulmonary hypertension in awake sheep, but had different lymph flow responses which may be due to pulmonary vasoconstriction at different sites, i.e. greater precapillary constriction with LTC4 because Qlym did not change and greater postcapillary constriction with LTD4 because Qlym increased with the same rise in Ppa.     

Evaluation of Starling forces in the equine digit

A pump-perfused extracorporeal digital preparation was used to evaluate blood flow, arterial pressure, venous pressure, isogravimetric capillary filtration coefficient, capillary pressure, and vascular compliance in six normal horses. From these data, pre- and postcapillary resistances and pre- and postcapillary resistance ratios were determined. Vascular and tissue oncotic pressures were estimated from plasma and lymph protein concentrations, respectively. By use of the collected and calculated data, tissue pressure in the digit was calculated using the Starling equation. In the isolated equine digit, isogravimetric capillary pressure averaged 36.7 mmHg, plasma and lymph oncotic pressures averaged aged 19.12 and 6.6 mmHg, respectively, interstitial fluid pressure averaged 25.6 mmHg, and the capillary filtration coefficient averaged 0.0013 ml.min-1.mm-1.100 g-1. Our results indicate that digital capillary pressure in the laterally recumbent horse is much higher than in analogous tissues in other species such as dog and human. However, the potential edemagenic effects of this high digital capillary pressure are opposed by at least two mechanisms: 1) a high tissue pressure and 2) a low microvascular surface area for fluid exchange and/or a low microvascular permeability to filtered fluid.     

Low-density lipoprotein concentration in interstitial fluid from human atherosclerotic lesions. Relation to theories of endothelial damage and lipoprotein binding

Increased endothelial permeability to low-density lipoprotein (LDL) is believed to be an initiating factor for atherosclerotic lesions. Concentrations of LDL, alpha 2-macroglobulin and albumin were measured by immunoassay in interstitial fluid collected from normal intima and atherosclerotic lesions of human aortas. The concentration of LDL in interstitial fluid from normal intima was twice the concentration in the patient's serum. In early proliferative (gelatinous) lesions the amount of interstitial fluid was consistently increased but its LDL concentration varied between 80 and 200% of adjacent normal intima. Highest concentrations of LDL were found in interstitial fluid from more advanced proliferative lesions, but the amount was reduced, suggesting a shift in tissue water. LDL was consistently low in interstitial fluid from fatty streaks comprised of lipid-filled cells, and in four of 12 lesions it was absent although alpha 2-macroglobulin and albumin concentrations were normal. Electrophoretic mobility of LDL, reflecting surface charge, was unchanged or increased in interstitial fluid from normal intima and fatty streaks, but decreased in gelatinous lesions. The ratio of LDL to alpha 2-macroglobulin and albumin in interstitial fluid was higher than in adjacent intact tissue. The results do not support the idea that increased endothelial permeability to LDL initiates atherogenesis.     

Poroelastic analysis of interstitial fluid flow in a single lamellar trabecula subjected to cyclic loading

Trabecula, an anatomical unit of the cancellous bone, is a porous material that consists of a lamellar bone matrix and interstitial fluid in a lacuno-canalicular porosity. The flow of interstitial fluid caused by deformation of the bone matrix is believed to initiate a mechanical response in osteocytes for bone remodeling. In order to clarify the effect of the lamellar structure of the bone matrix—i.e., variations in material properties—on the fluid flow stimuli to osteocytes embedded in trabeculae, we investigated the mechanical behavior of an individual trabecula subjected to cyclic loading based on poroelasticity. We focused on variations in the trabecular permeability and developed an analytical solution containing both transient and steady-state responses for interstitial fluid pressure in a single trabecular model represented by a multilayered two-dimensional poroelastic slab. Based on the obtained solution, we calculated the pressure and seepage velocity of the interstitial fluid in lacuno-canalicular porosity, within the single trabecula, under various permeability distributions. Poroelastic analysis showed that a heterogeneous distribution of permeability produces remarkable variations in the fluid pressure and seepage velocity in the cross section of the individual trabecula, and suggests that fluid flow stimuli to osteocytes are mostly governed by the value of permeability in the neighborhood of the trabecular surfaces if there is no difference in the average permeability in a single trabecula.     

Interstitial Fluid Flow and Drug Delivery in Vascularized Tumors: A Computational Model

Interstitial fluid is a solution that bathes and surrounds the human cells and provides them with nutrients and a way of waste removal. It is generally believed that elevated tumor interstitial fluid pressure (IFP) is partly responsible for the poor penetration and distribution of therapeutic agents in solid tumors, but the complex interplay of extravasation, permeabilities, vascular heterogeneities and diffusive and convective drug transport remains poorly understood. Here we consider–with the help of a theoretical model–the tumor IFP, interstitial fluid flow (IFF) and its impact upon drug delivery within tumor depending on biophysical determinants such as vessel network morphology, permeabilities and diffusive vs. convective transport. We developed a vascular tumor growth model, including vessel co-option, regression, and angiogenesis, that we extend here by the interstitium (represented by a porous medium obeying Darcy''s law) and sources (vessels) and sinks (lymphatics) for IFF. With it we compute the spatial variation of the IFP and IFF and determine its correlation with the vascular network morphology and physiological parameters like vessel wall permeability, tissue conductivity, distribution of lymphatics etc. We find that an increased vascular wall conductivity together with a reduction of lymph function leads to increased tumor IFP, but also that the latter does not necessarily imply a decreased extravasation rate: Generally the IF flow rate is positively correlated with the various conductivities in the system. The IFF field is then used to determine the drug distribution after an injection via a convection diffusion reaction equation for intra- and extracellular concentrations with parameters guided by experimental data for the drug Doxorubicin. We observe that the interplay of convective and diffusive drug transport can lead to quite unexpected effects in the presence of a heterogeneous, compartmentalized vasculature. Finally we discuss various strategies to increase drug exposure time of tumor cells.