R. Hsu

Analysis of the effects of oxygen supply and demand on hypoxic fraction in tumors.

The extent of hypoxic regions in a tumor tissue depends on the arrangement, blood flow rate and blood oxygen content of microvessels, and on the tissues oxygen consumption rate. Here, the effects of blood flow rate, blood oxygen content and oxygen consumption on hypoxic fraction are simulated theoretically, for a region whose microvascular geometry was derived from observations of a transplanted mammary andenocarcinoma (R3230AC) in a rat dorsal skin flap preparation. In the control state, arterial PO2 is 100 mmHg, consumption rate is 2.4 cm3 O2/100 g/min, and hypoxic fraction (tissue with PO2 < 3 mmHg) is 30%. Hypoxia is abolished by a reduction in consumption rate of at least 30%, relative to control, or an increase in flow rate by a factor of 4 or more, or an increase in arterial PO2 by a factor of 11 or more. These results suggest that reducing oxygen consumption rate may be more effective than elevating blood flow or oxygen content as a method to reduce tumor hypoxia.

Analysis of oxygen transport to tumor tissue by microvascular networks

We present theoretical simulations of oxygen delivery to tumor tissues by networks of microvessels, based on in vivo observations of vascular geometry and blood flow in the tumor microcirculation. The aim of these studies is to investigate the impact of vascular geometry on the occurrence of tissue hypoxia. The observations were made in the tissue (thickness 200 microns) contained between two glass plates in a dorsal skin flap preparation in the rat. Mammary adenocarcinomas (R3230 AC) were introduced and allowed to grow, and networks of microvessels in the tumors were mapped, providing data on length, geometric orientation, diameter and blood velocity in each segment. Based on these data, simulations were made of a 1 mm x 1 mm region containing five unbranched vascular segments and a 0.25 mm x 0.35 mm region containing 22 segments. Generally, vessels were assumed to lie in the plane midway between the glass plates, at 100 microns depth. Flow rates in the vessels were based on measured velocities and diameters. The assumed rate of oxygen consumption in the tissue was varied over a range of values. Using a Greens function method, partial pressure of oxygen (PO2) was computed at each point in the tissue region. As oxygen consumption is increased, tissue PO2 falls, with hypoxia first appearing at points relatively distant from the nearest blood vessel. The width of the well-oxygenated region is comparable to that predicted by simpler analyses. Cumulative frequency distributions of tissue PO2 were compared with predictions of a Krogh-type model with the same vascular densities, and it was found that the latter approach, which assumes a uniform spacing of vessels, may underestimate the extent of the hypoxic tissue. Our estimates of the maximum consumption rate that can be sustained without tissue hypoxia were substantially lower than those obtained from the Krogh-type model. We conclude that the heterogeneous structure of tumor microcirculation can have a substantial effect on the occurrence of hypoxic micro-regions.

Green's function methods for analysis of oxygen delivery to tissue by microvascular networks.

Delivery of oxygen to tissue is an essential function of the circulatory system. The distance that oxygen can diffuse into oxygen-consuming tissue is small, and so tissue oxygenation is critically dependent on the spatial arrangement of microvessels in tissue. Theoretical methods have been developed to simulate the spatial distribution of oxygen levels in tissue surrounding a network of microvessels. Here, numerical methods based on a Greens function approach are presented, for realistic three-dimensional network geometries derived from observations of skeletal muscle, brain, and tumor tissues. Relative to finite-difference methods, the Greens function approach reduces the number of unknowns in the numerical formulation and allows rapid computations even for complex vascular geometries. Generally, the boundary conditions on the exterior of the computational domain are not known. Imposition of a no-flux boundary condition can lead to exaggerated estimates of the extent of hypoxia in the tissue region. A new version of the method is described that avoids this problem and can be applied to arbitrarily shaped tissue domains.

A model for red blood cell motion in glycocalyx-lined capillaries

The interior surfaces of capillaries are lined with a layer (glycocalyx) of macromolecules bound or absorbed to the endothelium. Here, a theoretical model is used to analyze the effects of the glycocalyx on hematocrit and resistance to blood flow in capillaries. The glycocalyx is represented as a porous layer that resists penetration by red blood cells. Axisymmetric red blood cell shapes are assumed, and effects of cell membrane shear elasticity are included. Lubrication theory is used to compute the flow of plasma around the cell and within the glycocalyx. The effects of the glycocalyx on tube hematocrit (Fahraeus effect) and on flow resistance are predicted as functions of the width and hydraulic resistivity of the layer. A layer of width 1 micron and resistivity 10(8) dyn.s/cm4 leads to a relative apparent viscosity of approximately 10 in a 6-micron capillary at discharge hematocrit 45% and flow velocity of approximately 1 mm/s. This is consistent with experimental observations of increased flow resistance in microvessels in vivo, relative to glass tubes with the same diameters.

Theoretical Simulation of Oxygen Transport to Tumors by Three-Dimensional Networks of Microvessels

Hypoxic cancer cells are resistant to irradiation and to some types of chemotherapy. When these methods are used to treat tumors, presence of small hypoxic regions can compromise the outcome. Therefore, it is important to understand the factors determining the partial pressure of oxygen (pO2) in tumors at the microscopic level. These factors include the oxygen supply to the tumor, which depends on the flow rate and oxygen content of blood and the arrangement of microvessels within the tissue, and the rate of oxygen consumption by the tissue (Vaupel, 1979).

Theoretical Simulation of Oxygen Transport to Brain by Networks of Microvessels: Effects of Oxygen Supply and Demand on Tissue Hypoxia

Objective: Simulations of oxygen delivery by a three‐dimensional network of microvessels in rat cerebral cortex were used to examine how the distribution of partial pressure of oxygen (PO2) in tissue depends on blood flow and oxygen consumption rates.

Blood Flow and Red Blood Cell Deformation in Nonuniform Capillaries: Effects of the Endothelial Surface Layer

Objective: A theoretical model is used to examine the mechanics of red blood cell (RBC) motion in nonuniform capillaries. The model includes effects of the endothelial surface layer (ESL), which is a layer of macromolecules adjacent to the endothelium and which impedes plasma flow.

Angiogenesis: an adaptive dynamic biological patterning problem.

Formation of functionally adequate vascular networks by angiogenesis presents a problem in biological patterning. Generated without predetermined spatial patterns, networks must develop hierarchical tree-like structures for efficient convective transport over large distances, combined with dense space-filling meshes for short diffusion distances to every point in the tissue. Moreover, networks must be capable of restructuring in response to changing functional demands without interruption of blood flow. Here, theoretical simulations based on experimental data are used to demonstrate that this patterning problem can be solved through over-abundant stochastic generation of vessels in response to a growth factor generated in hypoxic tissue regions, in parallel with refinement by structural adaptation and pruning. Essential biological mechanisms for generation of adequate and efficient vascular patterns are identified and impairments in vascular properties resulting from defects in these mechanisms are predicted. The results provide a framework for understanding vascular network formation in normal or pathological conditions and for predicting effects of therapies targeting angiogenesis.

A Green's function method for analysis of oxygen delivery to tissue by microvascular networks

A theoretical model is formulated for analyzing oxygen delivery from an arbitrary network configuration of cylindrical microvessels to a finite region of tissue. In contrast to models based on the classical Krogh cylinder approach, this model requires no a priori assumptions concerning the extent of the tissue region supplied with oxygen by each vessel segment. Steady-state conditions are assumed, and oxygen consumption in the tissue is assumed to be uniform. The nonlinear dissociation characteristics of oxyhemoglobin are taken into account. A computationally efficient Greens function approach is used, in which the tissue oxygen field is expressed in terms of the distribution of source strengths along each segment. The utility of the model is illustrated by analyses of oxygen delivery to a cuboidal tissue region by a single segment and by a six-segment network. It is found that the fractional contribution of the proximal segments to total oxygen delivery increases with decreasing flow rate and metabolic rate.

Analysis of red blood cell motion through cylindrical micropores: effects of cell properties

Filtration through micropores is frequently used to assess red blood cell deformability, but the dependence of pore transit time on cell properties is not well understood. A theoretical model is used to simulate red cell motion through cylindrical micropores with diameters of 3.6, 5, and 6.3 microns, and 11-microns length, at driving pressures of 100-1000 dyn/cm2. Cells are assumed to have axial symmetry and to conserve surface area during deformation. Effects of membrane shear viscosity and elasticity are included, but bending resistance is neglected. A time-dependent lubrication equation describing the motion of the suspending fluid is solved, together with the equations for membrane equilibrium, using a finite difference method. Predicted transit times are consistent with previous experimental observations. Time taken for cells to enter pores represents more than one-half of the transit time. Predicted transit time increases with increasing membrane viscosity and with increasing cell volume. It is relatively insensitive to changes in internal viscosity and to changes in membrane elasticity except in the narrowest pores at low driving pressures. Elevating suspending medium viscosity does not increase sensitivity of transit time to membrane properties. Thus filterability of red cells is sensitively dependent on their resistance to transient deformations, which may be a key determinant of resistance to blood flow in the microcirculation.