Some factors affecting oxygen uptake by red blood cells in the pulmonary capillaries

Numerical analysis of oxygen uptake processes by red blood cells in stopped-flow measurements: Effects of cell shape, membrane permeability and unstirred layer

The transport process of oxygen and other gas species across red blood cell (RBC) membrane is of great importance for better understanding the critical biological functions of RBCs, and the stopped-flow experiments have often been employed for such investigations. In previous stopped-flow analyses, the RBC had usually been represented by a spherical capsule based on the RBC volume, and an assumed unstirred layer (USL) thickness had been used to determine the membrane permeability. In this research, unlike these previous studies, we simulate the oxygen uptake process with different RBC shapes (shperical, ellipsoidal and biconcave) and examine the effects of USL thickness and membrane permeability over broad ranges based on literature values. Our results show that the excess membrane area can greatly improve the oxygen transport efficiency, and a same uptake half-time can be obtained using different combinations of USL thickness and membrane permeability. These findings raise concerns on the reliability and uncertainty for the results and conclusions in previous studies, and also call for more complete numerical models, for example, with the fluid flow and cell deformation considered, and more in-depth investigations on the oxygen transport processes.

Oxygen transport across tank-treading red blood cell: Individual and joint roles of flow convection and oxygen-hemoglobin reaction

Gas, especially oxygen, transport in the microcirculation is a complex phenomenon, however, of critical importance for maintaining normal biological functions, and the cytoplasm fluid in red blood cells (RBCs) is the major vehicle for transporting oxygen from lungs to tissues via the circulatory system. Existing theoretical and numerical studies have neglected the cytoplasm convection effect by treating RBCs as rigid particles undergoing a constant translation velocity. As a consequence, the influence and mechanism of the cytoplasm flow on oxygen transport are still not clear in microcirculation research. In this study, we consider a tank-treading capsule in shear flow, which is generated with two parallel plates moving in opposite directions: the top plate of a higher oxygen pressure (P_O2) representing the RBC core in the central region of a microvessel and the bottom plate of a lower P_O2 representing the microvessel wall. Numerical simulations are conducted to investigate the individual and combined effects of cytoplasm convection and oxygen-hemoglobin (O₂-Hb) reaction on the oxygen transport efficiency across the tank-treading capsule, and different P_O2 situations and shear rates are also tested. Due to the lower oxygen diffusivity in cytoplasm, the presence of the capsule reduces the oxygen transfer flux across the gap by 7.34 % in the pure diffusion system where the flow convection and O₂-Hb reaction are both neglected. Including the flow convection or the O₂-Hb reaction has little influence on the oxygen flux; however, when they act together as in real microcirculation situations, the enhancement in oxygen transport could be significant, especially in the low P_O2 and high shear rate situations. In particular, with the respective P_O2 at 60 and 30 mmHg on the top and bottom plates and a 400 s^-1 shear rate, the oxygen flux reduction is only 0.02 %, suggesting that the cytoplasm convection can improve the oxygen transport across RBCs considerably. The simulation results are scrutinized to explore the underlying mechanism for the enhancement, and a new nondimensional parameter is introduced to characterize the importance of cytoplasm convection in oxygen transport. These simulation results, discussion and analysis could be helpful for a better understanding of the complex oxygen transport process and therefor valuable for relevant studies.

Reproducing the Hemoglobin Saturation Profile, a Marker of the Blood Oxygenation Level Dependent (BOLD) fMRI Effect, at the Microscopic Level

The advent of functional MRI in the mid-1990s has catalyzed progress pertaining to scientific discoveries in neuroscience. With the prospect of elucidating the physiological aspect of the Blood Oxygenation Level Dependent (BOLD) effect we present a computational capillary-tissue system capable of mapping venous hemoglobin saturation— a marker of the BOLD hemodynamic response. Free and facilitated diffusion and convection for hemoglobin and oxygen are considered in the radial and axial directions. Hemoglobin reaction kinetics are governed by the oxyhemoglobin dissociation curve. Brain activation, mimicked by dynamic transitions in cerebral blood velocity (CBv) and oxidative metabolism (CMR_O2), is simulated by normalized changes in m = (ΔCBv/CBv)/(ΔCMR_O2/CMR_O2) of values 2, 3 and 4. Venous hemoglobin saturation profiles and peak oxygenation results, for m = 2, based upon a 50% and a 25% increase in CBv and CMR_O2, respectively, lie within physiological limits exhibiting excellent correlation with the BOLD signal, for short-duration stimuli. Our analysis suggests basal CBv and CMR_O2 values of 0.6 mm/s and 200 μmol/100g/min. Coupled CBv and CMR_O2 responses, for m = 3 and m = 4, overestimate peak hemoglobin saturation, confirming the system’s responsiveness to changes in hematocrit, CBv and CMR_O2. Finally, factoring in neurovascular effects, we show that no initial dip will be observed unless there is a time delay in the onset of increased CBv relative to CMR_O2.

Theoretical analysis of the determinants of lung oxygen diffusing capacity

The process of pulmonary oxygen uptake is analyzed to obtain an explicit equation for lung oxygen diffusing capacity in terms of hematocrit and pulmonary capillary diameter. An axisymmetric model with discrete cylindrical erythrocytes is used to represent radial diffusion of oxygen from alveoli through the alveolar-capillary membrane into pulmonary capillaries, through the plasma, and into erythrocytes. Analysis of unsteady diffusion due to the passage of the erythrocytes shows that transport of oxygen through the alveolar-capillary membrane occurs mainly in the regions adjacent to erythrocytes, and that oxygen transport through regions adjacent to plasma gaps can be neglected. The model leads to an explicit formula for diffusing capacity as a function of geometric and oxygen transport parameters. For normal hematocrit and a capillary diameter of 6.75 μm, the predicted diffusing capacity is 102 ml O₂ min⁻¹ mmHg⁻¹. This value is 30-40% lower than values estimated previously by the morphometric method, which considers the total membrane area and the specific uptake rate of erythrocytes. Diffusing capacity is shown to increase with increasing hematocrit and decrease with increasing capillary diameter and increasing thickness of the membrane. Simulations of pulmonary oxygen uptake in humans under conditions of exercise or hypoxia based show closer agreement with experimental data than previous models, but still overestimate oxygen uptake. The remaining discrepancy may reflect effects of heterogeneity of perfusion and ventilation in the lung.

A theoretical comparison of macroscopic and microscopic modeling of singlet oxygen during Photofrin and HPPH mediated-PDT

Mathematic models were developed to simulate the complex dynamic process of photodynamic therapy (PDT). Macroscopic or microscopic modeling of singlet oxygen (1O2) is particularly of interest because it is the major cytotoxic agent causing biological effects during PDT. Our previously introduced macroscopic PDT model incorporates the diffusion equation for the light propagation in tissue and the macroscopic kinetic equations for the production of the 1O2. The distance-dependent distribution of 3O2 and reacted 1O2 can be numerically calculated using finite-element method (FEM). We recently improved the model to include microscopic kinetic equations of oxygen diffusion from uniformly distributed blood vessels and within tissue. In the model, the cylindrical blood capillary has radius in the range of 2-5 μm and a mean length of 300 μm, and supplies oxygen into tissue. The blood vessel network is assumed to form a 2-D square grid perpendicular to a linear light source. The spacing of the grid is 60 μm. Oxygen can also diffuse along the radius and the longitudinal axial of the cylinder within tissue. The oxygen depletion during Photofrin-PDT can be simulated using both macroscopic and microscopic approaches. The comparison of the simulation results have reasonable agreements when velocity of blood flow is reduced during PDT.

Some factors affecting pulmonary oxygen transport

Equations governing oxygen transport from alveolar gas to red blood cells flowing through pulmonary capillaries are written down. Some analytical predictions are made on factors affecting the rate at which this process takes place. Numerical simulations are carried out to investigate the effect of red blood cell shape, capillary dimensions, haematocrit and choice of oxygen dissociation curve on pulmonary oxygen transport. These factors all have an effect on pulmonary transport, with the effect being much more marked for simulations with low oxygen levels, typical of those seen in some subjects with respiratory disease.

Diffusing capacity reexamined: relative roles of diffusion and chemical reaction in red cell uptake of O2, CO, CO2, and NO

This paper presents an analytical expression for the diffusing capacity (Θt) of the red blood cell (RBC) for any reactive gas in terms of size and shape of the RBC, thickness of the unstirred plasma layer surrounding the RBC, diffusivities and solubilities of the gas in RBC and boundary layer, hematocrit, and the slope of the dissociation curve. The expression for Θthas been derived by spatial averaging of the fundamental convection-diffusion-reaction equation for O2in the RBC and has been generalized to all cell shapes and for other reactive gases such as CO, NO, and CO2. The effects of size and shape of the RBC, thickness of the unstirred plasma layer, hemoglobin concentration, and hematocrit on Θthave been analyzed, and the analytically obtained expression for Θthas been validated by comparison with different sets of existing experimental data for O2and CO2. Our results indicate that the discoidal shape of the human RBC with average dimensions of 1.6-μm thickness and 8-μm diameter is close to optimal design for O2uptake and that the true reaction velocity in the RBC is suppressed significantly by the mass transfer resistance in the surrounding unstirred layer. In vitro measurements using rapid-mixing technique, which measures Θtin the presence of artificially created large boundary layers, substantially underpredicts the in vivo diffusing capacity of the RBC in the diffusion-controlled regime. Depending on the conditions in the RBC, uptake of less reactive gases (such as CO) undergoes transition from reaction-limited to diffusion-limited regime. For a constant set of morphological parameters, the theoretical expression for Θtpredicts that Θt,NO> Θt,CO2> Θt,O2> Θt,CO.

Global BOLD MRI changes to ventilatory challenges in congenital central hypoventilation syndrome

We evaluated global blood oxygen level dependent (BOLD) signal changes in gray and white matter in 14 congenital central hypoventilation syndrome (CCHS) and 14 control subjects. One baseline image series with room air and three series with 30 s room air followed by 120 s hypercapnia (5% CO2/95% O2), hypoxia (15% O2/85% N2) or hyperoxia (100% O2) were collected. Hypercapnia and hyperoxia raised, and hypoxia lowered gray and white matter global signal in both groups, with smaller changes in white matter. Signal changes in CCHS cases were lower than control subjects for hypercapnia in gray and white matter, slightly more-enhanced in hypoxia, and, except for initial transient responses, were nearly comparable during hyperoxia. Initial signal rate or pattern changes emerged in all three challenges in gray or white matter in control, but not CCHS cases. Neural or vascular mechanisms mediate perfusion differently in CCHS; the aberrant initial transient responses may reflect deficiencies in rapidly-varying physiologic interactions in the syndrome.