AJP - Heart Add DOIs to your references at manuscript stage!
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


Am J Physiol Heart Circ Physiol 278: H137-H150, 2000;
0363-6135/00 $5.00
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Web of Science (10)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Audi, S. H.
Right arrow Articles by Dawson, C. A.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Audi, S. H.
Right arrow Articles by Dawson, C. A.
Vol. 278, Issue 1, H137-H150, January 2000

Toluidine blue O and methylene blue as endothelial redox probes in the intact lung

Said H. Audi1,2, Lars E. Olson1,3, Robert D. Bongard3, David L. Roerig4,5,6, Marie L. Schulte3, and Christopher A. Dawson1,2,3,4,5,6

1 Department of Biomedical Engineering, Marquette University, Milwaukee 53201-1881; Departments of 2 Pulmonary and Critical Care Medicine, 3 Physiology, 4 Anesthesiology, and 5 Pharmacology and Toxicology, Medical College of Wisconsin, Milwaukee 53226; and 6 Zablocki Veterans Affairs Medical Center, Milwaukee, Wisconsin 53295


    ABSTRACT
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES
APPENDIX

There is increasing evidence that the redox activities of the pulmonary endothelial surface may have important implications for the function of both lungs and blood. Because of the inherent complexity of intact organs, it can be difficult to study these activities in situ. Given the availability of appropriate indicator probes, the multiple-indicator dilution (MID) method is one approach for dealing with some aspects of this complexity. Therefore, the objectives of the present study were to 1) evaluate the potential utility of two thiazine redox indicators, methylene blue (MB) and toluidine blue O (TBO), as MID electron acceptor probes for in situ pulmonary endothelium and 2) develop a mathematical model of the pulmonary disposition of these indicators as a tool for quantifying their reduction on passage through the lungs. Experiments were carried out using isolated rabbit lungs perfused with physiological salt solution with or without plasma albumin over a range of flow rates. A large fraction of the injected TBO disappeared from the perfusate on passage through the lungs. The reduction of its oxidized, strongly polar, relatively hydrophilic blue form to its colorless, highly lipophilic reduced form was revealed by the presence of the reduced form in the venous effluent when plasma albumin was included in the perfusate. MB was also lost from the perfusate, but the fraction was considerably smaller than for TBO. A distributed-in-space-and-time model was developed to estimate the reduction rate parameter, which was ~29 and 1.0 ml/s for TBO and MB, respectively, and almost flow rate independent for both indicators. The results suggest the utility particularly of TBO as an electron acceptor probe for MID studies of in situ pulmonary endothelium and of the model for quantitative evaluation of the data.

transplasma membrane electron transport; multiple-indicator dilution; mathematical modeling


    INTRODUCTION
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES
APPENDIX

THE PULMONARY ENDOTHELIUM is capable of reducing certain blood-borne electron acceptors as they pass through the pulmonary capillary bed. The mechanisms involved include transplasma membrane electron transport systems that utilize intracellular electron donors to reduce extracellular acceptors (1, 10, 12, 34, 50). These redox systems control the pulmonary disposition of certain redox active drugs (10), and, by analogy with the functions of transplasma membrane transport systems in other cell types, they may be mechanisms by which the endothelial cells influence and/or sense the redox status of the blood (15, 19, 22, 28, 32, 35-37, 46, 47, 51). Studying metabolic processes such as these electron transport systems within an intact functioning organ is complicated by the many factors that can influence the disposition of probes for such processes within an organ (4, 7, 8, 23, 29), and probe disposition can be affected by factors other than the targeted process. The bolus injection multiple-indicator dilution (MID) method has been used for studies of other pulmonary endothelial surface reactions (9, 16-18, 20, 21, 24, 30, 49). A principal feature of this MID method is that it provides the temporal information needed to identify separately the contributions of various factors that may influence the overall disposition of a probe for a particular metabolic process. Identification of useful probes is one step in the development of an MID method for a specific cellular function. The purpose of the present study was to evaluate two thiazine compounds, methylene blue (MB) and toluidine blue O (TBO), as electron acceptor probes for the pulmonary endothelium in intact lungs. Both were previously shown to be electron acceptors for transplasma membrane electron transport in bovine pulmonary arterial endothelial cells grown in cell culture (12, 34, 40). A useful MID probe must fulfill several criteria; some are specific to the particular study conditions, but some are general. An important general criterion is that the probe must be a good enough substrate that the rate of the reaction under study is fast enough to have a significant impact on probe disposition within the 0.5- to 2-s pulmonary capillary transit time but slow enough that the reaction kinetics of interest, rather than only the rate of probe delivery, dominates the probe disposition (4, 23). Therefore, the determination of the kinetics of MB and TBO disposition on passage through the pulmonary circulation was a central focus of the present study. We addressed the problem of quantifying the kinetics of probe disposition on passage through the lungs by representing the hypothesized processes involved in the form of a mathematical model.


    METHODS
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES
APPENDIX

Isolated Lung Preparation

The isolated rabbit lung preparation used has been described previously (3, 10). Each of 22 New Zealand White rabbits [2.51 ± 0.23 (SD) kg body wt] was given chlorpromazine hydrochloride (25 mg/kg im) followed by pentobarbital sodium (15-20 mg/kg iv). A carotid artery catheter was inserted for heparinization (1,200 IU/kg) followed by exsanguination. The chest was opened, and cannulas were placed in the pulmonary artery, left atrium, and trachea. The cannulated lungs were removed and suspended vertically from the trachea, with the arterial and venous cannulas connected to a temperature-controlled (37°C) recirculating perfusion system primed with a physiological salt solution (PSS) containing (in mM) 4.7 KCl, 2.51 CaCl2, 1.19 MgSO4, 2.5 KH2PO4, 118 NaCl, 25 NaHCO3, and 5.5 glucose and, in most cases, either 5% BSA (referred to as PSS-albumin) or 5% dextran (average mol wt 74,200) (referred to as PSS-dextran). The perfusate was pumped (Masterflex roller pump) into the pulmonary artery from a reservoir into which the perfusate drained from the left atrium. The first 200 ml to pass through the lungs were used to clear the lungs of residual blood and were discarded before commencing recirculation. The standard flow rate was 200 ml/min, with the flow rate adjusted temporarily according to the experimental protocol as indicated in Bolus Injections. Pulmonary venous pressure was set equal to pleural pressure between adjustments required for a particular experimental protocol. Pulmonary arterial and left atrial pressures, measured relative to the level of the left atrium, were monitored continuously. The lungs were ventilated with a gas mixture of ~17% O2-5% CO2-balance N2 at a frequency of 11 breaths/min, with end-inspiratory and end-expiratory pressures of ~6 and 1 Torr, respectively, maintained with water overflow valves. The resulting perfusate PO2, PCO2, and pH were 123 ± 4.6 (SD) Torr, 36.2 ± 2.9 Torr, and 7.40 ± 0.05, respectively.

To measure the venous effluent indicator concentration versus time curves, the venous effluent could be diverted into the sample tubes of a Gilson-Escargot fraction collector modified to sample at equal6.67 samples/s. For experiments in which TBO was studied in PSS-albumin-perfused lungs, the venous effluent tubing passed through a previously described photodetector (40), measuring optical absorbance at wavelengths of 490 and 590 nm, before flowing into the fraction collector tubes. Each fraction collector tube contained 20 µl of 0.8 mM potassium ferricyanide solution so that complete oxidation of any reduced MB (MBH) or TBO (TBOH) in the effluent would occur before the absorbance was measured. Thus total MB, MBH + oxidized MB (MB+), total TBO, or TBOH + oxidized TBO (TBO+) concentration in the venous effluent was measured in the fraction collector tubes. The photodetector measured TBO+ rapidly enough after leaving the lung that autoxidation of any TBOH present was small. FITC-dextran, which was also measured by the photodetector, was included in the bolus for the purpose of determining the transport function between the two sampling sites as indicated in Bolus Injections.

Bolus Injections

To measure the disappearance of MB+, TBO+, or TBOH from a bolus during passage through the lungs, the ventilation was halted in end expiration and the venous outflow was diverted into the fraction collector. A 0.9-ml bolus containing reference indicators 125I-labeled human serum albumin (~0.09 µCi/ml) and/or FITC-dextran (1.5 mg/ml, average mol wt 2,000) and test indicators MB+, TBO+, or TBOH (0.1-0.15 nmol/ml) was then introduced into the arterial inflow tubing via an injection loop. Consecutive 2-ml samples were collected at a rate and duration appropriate for the flow rate used. The first sample was used as a background reference sample. Measured quantities of the injectate solution were added to the next three samples. All of these samples emerged in the venous effluent before the appearance of the injected indicators, and they were used as standards for determining the effluent concentrations of 125I, FITC, MB+, or TBO+ as fractions of their injected concentrations. A 1-ml aliquot of each collected sample was used to determine 125I concentration by gamma scintillation counting. The remainder was analyzed spectrophotometrically. The photodetector output was calibrated with known concentrations of FITC, TBO+, or MB+ in the perfusate solution. Eight percent of the TBO signal at 590 nm was measured at 490 nm, with no detectable FITC signal at 590 nm. The reference indicator recoveries based on the standards were not significantly different from 100%. For MB, the perfusate was PSS-dextran. For TBO, both PSS-dextran and PSS-albumin were used.

To prepare the TBOH injectates, 3 ml of the respective perfusate containing 160 µM TBO, 1.5 mg FITC-dextran and 200 µM NADH were added to a test tube having two ports. The mixture was bubbled through one of the ports for 10 min with 5.8% CO2 in N2 to deoxygenate while maintaining the pH at 7.4. Diaphorase (0.2 U in 0.02 ml) was then added to the solution under anoxic conditions. After the disappearance of blue color from the solution, the ports were sealed and low-dead space syringes were filled with the injectate directly from the sealed tube. The injection loop was prefilled with 3 ml of deoxygenated perfusate just before the addition of the 0.9-ml TBOH injectate.

At the end of each experiment, the lungs were removed and the pulmonary arterial and left atrial cannulas were connected directly together. Boluses containing the 125I-labeled albumin and FITC-dextran were injected and fraction collector samples and photodetector output were collected at each flow rate used in the experiment. At least one such bolus also included the respective thiazine compound. The data obtained from these samples were used to measure the part of the bolus transit time and dispersion that was caused by the injection, tubing, and sampling system and to affirm that any separation between the thiazine test indicator concentration curves and the reference indicator concentration curves was caused only by passage through the lungs and not some interaction with the perfusion system connected to the lungs.

To measure the indicator concentration versus time curves over a range of flow rates, the first step was to raise the flow rate to the highest level to be used in the experiment, with the left atrial pressure set equal to atmospheric pressure at the level of the left atrium. The resulting arterial pressure divided by 2 was the value at which the average of the arterial left atrial pressure was set for each of the lower flow rates by adjusting the height of the reservoir. This was done so that the vascular volume was approximately the same at each flow rate. The flow rate was then set at ~200, 400, 600, or 800 ml/min in TBO experiments or 50, 100, 200, or 300 ml/min in MB experiments, and a bolus was injected and samples collected as indicated above. Injections were then made at two or three of the other designated flow rates in similar fashion.

To determine whether dextran in the perfusate might have any influence on TBO or MB disposition in the lungs, after the above-described protocol had been completed in one of the lungs from the TBO and MB groups perfused with PSS-dextran, the perfusate was changed to PSS with neither dextran nor albumin. Then a TBO+ or MB+ bolus (also dextran free) was injected at a flow rate of 800 or 300 ml/min, respectively. Also, in one lung perfused with PSS-dextran the photodetector was included to confirm that effluent MBH and TBOH concentrations were below detectable levels when there was no plasma albumin in the perfusate.

Additional Measurements

The equilibrium binding of TBO+ and TBOH to BSA was measured by ultrafiltration using the Amicon MPS-1 micropartition system with a YM30 membrane at room temperature as previously described (3). The concentration of TBO+ or TBOH was 12 µM in phosphate-buffered (pH 7.4) 0.9% NaCl solution containing 5% BSA. TBOH was prepared by adding 0.5 µmol sodium hydrosulfite/ml of TBO+ solution. TBOH was measured after reoxidation to TBO+. The albumin-bound fractions were 59 and 79% for TBO+ and TBOH, respectively.

The apparent octanol-to-water partition coefficient for TBOH and TBO+ at room temperature was determined as previously described (10) by using a concentration of 10.4 µM TBO+ or TBOH in 0.1 M phosphate buffer (pH 7.4). TBOH was prepared as indicated in Bolus Injections for the ultrafiltration studies. The values of the partition coefficients were 0.95 and 10.3 for TBO+ and TBOH, respectively. The previously reported values for MB+ and MBH were 0.22 and 16, respectively (10).

To determine the TBOH autoxidation rate, TBOH was prepared in a manner similar to the TBOH injectate described above. TBO and NADH (final concentration of 160 µM for both) were added to tubes containing 3 ml Hanks' balanced salt solution and 10 mM HEPES at pH 7.4 with or without 5% BSA or 5% dextran. Each solution was bubbled with N2 for 10 min, and diaphorase (0.2 U) was added. The resulting TBOH solution (0.1 ml) was added to a spectrophotometric cuvette containing 2 ml Hanks'-HEPES buffer with or without 5% BSA or 5% dextran at pH 7.4 and equilibrated with atmospheric PO2. Absorbance at 590 nm was recorded at 5-s intervals for 3 min in a Beckman DU 7400 spectrophotometer at 37°C.


    RESULTS
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES
APPENDIX

Experimental Results

The pressure, flow, and volume data reflecting the conditions under which the injections were carried out are given in Tables 1 and 2. The average pressures and vascular volumes were higher for the TBO than the MB experiments because of the higher flow rate range covered. The injections themselves had no detectable effect on perfusion pressures, suggesting that any effects that these compounds might have on vascular tone were below the level of detectability in these studies.

                              
View this table:
[in this window]
[in a new window]
 
Table 1.   Pressure, flow, and volume data for methylene blue experiments


                              
View this table:
[in this window]
[in a new window]
 
Table 2.   Pressure, flow, and volume data for toluidine blue O experiments

Examples of graphs of indicator concentrations normalized to their respective injected amounts versus time, for the various combinations of PSS-dextran, PSS-albumin, MB+, TBO+, and TBOH and at different flow rates, are shown in Figs. 1-5. With PSS-dextran perfusion, only the oxidized forms were detected in the venous effluent. The presence of TBOH in the effluent after TBO+ injections when the lungs were perfused with PSS-albumin is reflected by the fact that the total TBO (TBOH + TBO+) concentrations were higher than the concentrations of TBO+ alone (Fig. 3). The fractions of MB+ and TBO+ that disappeared from the PSS-dextran and PSS-albumin perfusates on passage through the lungs are summarized as extractions in Fig. 6. The TBO+ extractions were higher than the MB+ extractions even though the flow rate range was also higher for TBO+. The extraction of TBO+ from PSS-albumin was lower than from PSS-dextran, reflecting the effects of the TBO+-BSA binding. When the lungs were perfused with PSS with neither albumin nor dextran, the effluent TBO+ and MB+ were indistinguishable from those when the lungs were perfused with PSS-dextran, indicating that the dextran had no detectable effect on TBO+ or MB+ extraction. When TBOH rather than TBO+ was injected with PSS-albumin as the perfusate, the effluent TBOH concentration curves tended to be displaced to the right in comparison to the reference indicator or to the TBO+ curves after TBO+ injection (Fig. 4). This difference reflects the tissue distribution of the TBOH (see DISCUSSION). When TBOH was injected with PSS-dextran as the perfusate, the TBOH almost completely disappeared from the perfusate (Fig. 5).


View larger version (21K):
[in this window]
[in a new window]
 
Fig. 1.   Concentration (normalized to amount of injected indicator) vs. time curves for oxidized methylene blue (MB+) and reference indicator (CR) obtained at 4 flow rates in 1 lung perfused with albumin-free perfusate [physiological salt solution (PSS)-dextran]. Solid lines are model fits to data.



View larger version (18K):
[in this window]
[in a new window]
 
Fig. 2.   Concentration (normalized to amount of injected indicator) vs. time curves for oxidized toluidine blue O (TBO+) and reference indicator (CR) obtained at 3 flow rates in 1 lung perfused with albumin-free perfusate (PSS-dextran). Solid lines are model fits to data.



View larger version (21K):
[in this window]
[in a new window]
 
Fig. 3.   Concentration (normalized to amount of injected indicator) vs. time curves for TBO+, total TBO [reduced TBO (TBOH) + TBO+], and reference indicator (CR) obtained at 3 flow rates in 1 lung perfused with albumin-containing perfusate (PSS-albumin). Solid lines are model fits to data.



View larger version (18K):
[in this window]
[in a new window]
 
Fig. 4.   Concentration (normalized to amount of injected indicator) vs. time curves for TBOH and reference indicator (CR) obtained at 3 flow rates in 1 lung perfused with albumin-containing perfusate (PSS-albumin). Solid lines are model fits to data.



View larger version (17K):
[in this window]
[in a new window]
 
Fig. 5.   Concentration (normalized to amount of injected indicator) vs. time curves for TBOH and reference indicator (CR) obtained at 3 flow rates in 1 lung perfused with albumin-free perfusate (PSS-dextran). Solid lines are model fits to data.



View larger version (29K):
[in this window]
[in a new window]
 
Fig. 6.   Average extraction (±SE) of MB+ from PSS-dextran (n = 6-8) and of TBO+ from PSS-dextran (n = 8) or PSS-albumin (n = 4) vs. flow rate. Extraction = 1 - ([MB+] or [TBO+])/ CR(t), where concentrations were obtained at time of maximum of CR(t). Dashed lines are isopleths for constant values of -F × ln(1 - extraction) indicated.

Further evaluation of these results was carried out by expressing the hypothesized interactions between these thiazine compounds and the lung tissue in the following model.

Model

The thiazine dyes are assumed to participate in the following associations and reactions within the perfusate and tissue.

Stoichiometric equations.

where B+ is the oxidized thiazine test indicator; BH is the reduced thiazine test indicator; A represents the sites of sequestration of BH within the tissue; E is thiazine reductase; DH is the electron donor; D+ is the oxidized form of DH; P is plasma protein; Z represents the nonspecific B+ tissue binding sites; BHA is BH bound to A; BHP is BH bound to plasma protein; B+P is B+ bound to plasma protein; B+Z is B+ associated nonspecifically with tissue; k1 and k-1 are plasma protein association and dissociation rate constants, respectively, for B+; k2 and k-2 are plasma protein association and dissociation rate constants, respectively, for BH; ko is the autoxidation rate constant (k7[O2]1/2[H+]) for BH within the perfusate; kr1 and kr2 are the thiazine reductase association rate constants k5[DH] and k6[DH], for B+ and B+P, respectively; k3 and k-3 are the nonspecific tissue association and dissociation rate constants, respectively, for B+; and k4 is the tissue sequestration rate constant for BH.

Species balance equations. A single capillary element of the kinetic model is composed of a capillary volume (Qc) and a surrounding tissue volume (Qt). The spatial and temporal variations in the concentrations of the vascular reference indicator and B+ and BH are described by the following species balance equations based on the above stoichiometric equations and the following assumptions. 1) The reference indicator, having perfusate concentration [R], was convected through the pulmonary vascular bed without interacting with the tissue. 2) When the plasma protein (P) was present in the perfusate, the equilibration between P and B+ or BH was rapid relative to the capillary transit time (<OVL><IT>t</IT></OVL>c). 3) The endothelial cell reduction of B+ mediated by E was irreversible. 4) The association of B+ with Z was rapid relative to <OVL><IT>t</IT></OVL>c. 5) The equilibration of BH between Qc and Qt with tissue to perfusate partition coefficient lambda  was rapid relative to <OVL><IT>t</IT></OVL>c. 6) The PO2 and [DH] were constant under the study conditions.

SINGLE CAPILLARY ELEMENT.
<FR><NU>∂[R]</NU><DE>∂<IT>t</IT></DE></FR> + <IT>W</IT> <FR><NU>∂[R]</NU><DE>∂<IT>x</IT></DE></FR> = 0 (1)

<FR><NU>∂[B<SUP>+</SUP>]</NU><DE>∂<IT>t</IT></DE></FR> + <IT>W</IT> <FR><NU>∂[B<SUP>+</SUP>]</NU><DE>∂<IT>x</IT></DE></FR> = [B<SUP>+</SUP>P]<IT>k</IT><SUB>−1</SUB> − <IT>k</IT><SUB>1</SUB>[B<SUP>+</SUP>][P] 

− <FR><NU>E<IT>k</IT><SUB>rl</SUB>[B<SUP>+</SUP>]</NU><DE>Q<SUB>c</SUB></DE></FR> + <IT>k</IT><SUB>o</SUB>[BH] + <FR><NU>(<IT>k</IT><SUB>−3</SUB>B<SUP>+</SUP>Z − <IT>k</IT><SUB>3</SUB>[B<SUP>+</SUP>]Z)</NU><DE>Q<SUB>c</SUB></DE></FR> (2)

<FR><NU>∂[B<SUP>+</SUP>P]</NU><DE>∂<IT>t</IT></DE></FR> + <IT>W</IT> <FR><NU>∂[B<SUP>+</SUP>P]</NU><DE>∂<IT>x</IT></DE></FR> = −[B<SUP>+</SUP>P]<IT>k</IT><SUB>−1</SUB> 

+ <IT>k</IT><SUB>1</SUB>[B<SUP>+</SUP>][P] − <FR><NU>E<IT>k</IT><SUB>r2</SUB>[B<SUP>+</SUP>P]</NU><DE>Q<SUB>c</SUB></DE></FR> + <IT>k</IT><SUB>o</SUB>[BHP] (3)

<FR><NU>∂B<SUP>+</SUP>Z</NU><DE>∂<IT>t</IT></DE></FR> = − <IT>k</IT><SUB>−3</SUB> B<SUP>+</SUP>Z + <IT>k</IT><SUB>3</SUB> [B<SUP>+</SUP>] Z (4)

<FR><NU>∂[BH]</NU><DE>∂<IT>t</IT></DE></FR> + <IT>W</IT> <FENCE><FR><NU>Q<SUB>c</SUB></NU><DE>Q<SUB>c</SUB> + &lgr;Q<SUB>t</SUB></DE></FR></FENCE> <FR><NU>∂[BH]</NU><DE>∂<IT>x</IT></DE></FR> = <FENCE><FR><NU>Q<SUB>c</SUB></NU><DE>Q<SUB>c</SUB> + &lgr;Q<SUB>t</SUB></DE></FR></FENCE>

⋅ ([BHP]<IT>k</IT><SUB>−2</SUB> − <IT>k</IT><SUB>2</SUB>[BH][P]) − <FENCE><FR><NU>Q<SUB>c</SUB></NU><DE>Q<SUB>c</SUB> + &lgr;Q<SUB>t</SUB></DE></FR></FENCE> <IT>k</IT><SUB>o</SUB>[BH] 

+ <FR><NU>E<IT>k</IT><SUB>rl</SUB>[B<SUP>+</SUP>]</NU><DE>(Q<SUB>c</SUB> + &lgr;Q<SUB>t</SUB>)</DE></FR> − <FR><NU>&lgr;Q<SUB>t</SUB></NU><DE>(Q<SUB>c</SUB> + &lgr;Q<SUB>t</SUB>)</DE></FR> (<IT>k</IT><SUB>4</SUB>[A][BH]) (5)

<FR><NU>∂[BHP]</NU><DE>∂<IT>t</IT></DE></FR> + <IT>W</IT> <FR><NU>∂[BHP]</NU><DE>∂<IT>x</IT></DE></FR> = − [BHP]<IT>k</IT><SUB>−2</SUB> 

+ <IT>k</IT><SUB>2</SUB>[BH][P] − <IT>k</IT><SUB>o</SUB> [BHP] + <FR><NU>E<IT>k</IT><SUB>r2</SUB>[B<SUP>+</SUP>P]</NU><DE>Q<SUB>c</SUB></DE></FR> (6)
E = [E]Qc, Z = [Z]Qc, and B+Z = [B+Z]Qc, where [E], [Z], and [B+Z] are the moles of thiazine reductase and free and bound nonspecific tissue binding sites per milliliter of capillary volume, respectively. [R](x,t), [B+](x,t), [BH](x,t), [BHP](x,t), and [B+P](x,t) are the respective vascular concentrations at distance x from the capillary inlet (x = 0) and time t. W is the average linear flow velocity within Qc.

Under the assumption of rapid equilibration between the free and plasma protein-bound B+ and BH

[BHP] = <FR><NU>[BH][P]</NU><DE><IT>K</IT><SUB>2</SUB></DE></FR>
and
[B<SUP>+</SUP>P] = <FR><NU>[B<SUP>+</SUP>][P]</NU><DE><IT>K</IT><SUB>1</SUB></DE></FR> (7)
where K1 = k-1/k1 and K2 = k-2/k2 are the plasma protein equilibrium dissociation constants for B+ and BH, respectively.

Under the assumption of rapid equilibration between B+ and Z
<IT>k</IT><SUB>−3</SUB> B<SUP>+</SUP>Z = <IT>k</IT><SUB>3</SUB> [B<SUP>+</SUP>]Z
and
B<SUP>+</SUP>Z = <FR><NU>[B<SUP>+</SUP>]</NU><DE>Q<SUB>f</SUB></DE></FR> (8)
where Qf = k-3/Zk3 has units of volume and acts as a virtual volume of distribution for B+.

Adding Eqs. 2-4 and 5-6 and substituting Eqs. 7 and 8 results in
<FR><NU>∂[<A><AC>B</AC><AC>˜</AC></A><SUP>+</SUP>]</NU><DE>∂<IT>t</IT></DE></FR> + <IT>W</IT> <FENCE><FR><NU>Q<SUB>c</SUB></NU><DE>Q<SUB>c</SUB> + Q<SUB>F</SUB></DE></FR></FENCE> <FR><NU>∂[<A><AC>B</AC><AC>˜</AC></A><SUP>+</SUP>]</NU><DE>∂<IT>x</IT></DE></FR> 

= − <FR><NU>[<A><AC>B</AC><AC>˜</AC></A><SUP>+</SUP>] <IT>K</IT><SUB>red</SUB></NU><DE>Q<SUB>c</SUB> + Q<SUB>F</SUB></DE></FR> + <FR><NU>Q<SUB>c</SUB> [<A><AC>BH</AC><AC>∼</AC></A>] <IT>k</IT><SUB>o</SUB></NU><DE>Q<SUB>c</SUB> + Q<SUB>F</SUB></DE></FR> (9)

<FR><NU>∂[<A><AC>BH</AC><AC>∼</AC></A>]</NU><DE>∂<IT>t</IT></DE></FR> + <IT>W</IT> <FENCE><FR><NU>Q<SUB>c</SUB></NU><DE>Q<SUB>c</SUB> + Q<SUB>T</SUB></DE></FR></FENCE> <FR><NU>∂[<A><AC>BH</AC><AC>∼</AC></A>]</NU><DE>∂<IT>x</IT></DE></FR>

= <FR><NU>[<A><AC>B</AC><AC>˜</AC></A><SUP>+</SUP>] <IT>K</IT><SUB>red</SUB></NU><DE>Q<SUB>c</SUB> + Q<SUB>T</SUB></DE></FR> − <FR><NU>Q<SUB>c</SUB> [<A><AC>BH</AC><AC>∼</AC></A>] <IT>k</IT><SUB>o</SUB></NU><DE>Q<SUB>c</SUB> + Q<SUB>T</SUB></DE></FR> − <FR><NU>[<A><AC>BH</AC><AC>∼</AC></A>] <IT>K</IT><SUB>seq</SUB></NU><DE>Q<SUB>c</SUB> + Q<SUB>T</SUB></DE></FR> (10)
where [<A><AC>BH</AC><AC>∼</AC></A>] = [BH]alpha 2, [&Btilde;+] = [B+]alpha 1, alpha 1 = 1 + [P]/K1, and alpha 2 = 1 + [P]/K2. The model parameters are QF = Qf/alpha 1 (ml), which is the measure of the magnitude of the rapidly equilibrating nonspecific interactions of B+ with the tissue; QT = lambda Qt/alpha 2 (ml), which is the measure of the magnitude of the rapidly equilibrating partitioning of the BH between tissue and perfusate; Kred = (Ekr1 + Ekr2 [P]/K1)/alpha 1 (ml/s), which is the measure of the B+ reduction rate; Kseq = [A]k4QT (ml/s), which is the measure of the BH sequestration rate within the tissue; and ko (s-1), which is the rate of autoxidation of BH within the perfusate.

Whole organ. To construct an organ model from the single capillary element model, the distribution of pulmonary capillary transit times [hc(t)] needs to be taken into account (4, 29). Previously, we (5) estimated that for normal rabbit lungs in this perfusion system, the pulmonary capillary mean transit time (<OVL><IT>t</IT></OVL>c) was ~44% of the total vascular mean transit time ( <OVL><IT>t</IT></OVL> ), the relative dispersion of hc(t) (RDc = sigma c/<OVL><IT>t</IT></OVL>c) was ~0.9, and the skewness coefficient of hc(t) (m3c/sigma 3c) was ~2, where m3c and sigma 3c are the third central moment and standard deviation of hc(t), respectively. For the present analysis, we used these values to approximate hc(t) using a shifted random walk function, which is a probability density function whose functional values are determined by these three moments as previously described (2).

The organ reference indicator outflow concentration vs. time curve [CR(t) = (q/F)hc(t)*hn(t), where * is the convolution operator, q is the mass of the injected indicator, F is the total flow through the organ, and hn(t) is the noncapillary (arteries, veins, connecting tubing and the injection system) transit time distribution] was obtained. The hn(t) was also represented by a shifted random walk function whose parameters were specified by iteratively convolving a trial hn(t) with hc(t) until the optimal least-squares fit to CR(t) was obtained (2).

To estimate model parameters from the data from each injection, Eqs. 9 and 10 were solved numerically for the appropriate boundary conditions (given in Estimation of Model Parameters for Specific Experiments) using the finite-difference method at each iteration of a Levenberg-Marquardt optimization routine (4, 31). The solution is for a single capillary element with Cin(t) as the capillary input concentration curve. The model solution for a single capillary having the maximum capillary transit time also provides the output for all capillary transit times between the minimum and maximum capillary transit times (4). To provide the whole organ output for vascular reference indicator [CR(t)], and test indicator [C(t)], the outputs for all transit times are summed, each weighted according to hc(t) (4).

Autoxidation of BH in arteries, veins, and connecting tubing. Because BH autoxidation can occur outside the capillary region within the conducting arteries and veins and the tubing connecting the lungs to the injection and sampling sites, this autoxidation was also addressed in the model. In this part of the system, Eqs. 9 and 10 reduce to
<FR><NU>∂[<A><AC>BH</AC><AC>∼</AC></A>]</NU><DE>∂<IT>t</IT></DE></FR> + <IT>W</IT> <FR><NU>∂[<A><AC>BH</AC><AC>∼</AC></A>]</NU><DE>∂<IT>x</IT></DE></FR> = − <IT>k</IT><SUB>o</SUB> [<A><AC>BH</AC><AC>∼</AC></A>] (11)

<FR><NU>∂[<A><AC>B</AC><AC>˜</AC></A><SUP>+</SUP>]</NU><DE>∂<IT>t</IT></DE></FR> + <IT>W</IT> <FR><NU>∂[<A><AC>B</AC><AC>˜</AC></A><SUP>+</SUP>]</NU><DE>∂<IT>x</IT></DE></FR> = <IT>k</IT><SUB>o</SUB> [<A><AC>BH</AC><AC>∼</AC></A>] (12)
which can be analytically solved for [<A><AC>BH</AC><AC>˜</AC></A>](&xgr;,<IT>t</IT>) and [&Btilde;+](xi ,t), where x = xi  is the inlet to either the capillaries or to the photodetector
[<A><AC>BH</AC><AC>∼</AC></A>] (&xgr;,<IT>t</IT>) = [<A><AC>BH</AC><AC>∼</AC></A>]<SUB>i</SUB>(0,<IT>t</IT> − <OVL><IT>t</IT></OVL><SUB>o</SUB>)<IT>e</IT><SUP>−<OVL><IT>t</IT></OVL><SUB>o</SUB><IT>k</IT><SUB>o</SUB></SUP> for <IT>t</IT> > <OVL><IT>t</IT></OVL><SUB>o</SUB> (13a)

[<A><AC>B</AC><AC>˜</AC></A><SUP>+</SUP>](&xgr;,<IT>t</IT>) = [<A><AC>B</AC><AC>˜</AC></A><SUP>+</SUP>]<SUB>i</SUB>(0,<IT>t</IT> − <OVL><IT>t</IT></OVL><SUB>o</SUB>) 

+ [<A><AC>BH</AC><AC>∼</AC></A>]<SUB>i</SUB>(0,<IT>t</IT> − <OVL><IT>t</IT></OVL><SUB>o</SUB>) (1 − <IT>e</IT><SUP>−<OVL><IT>t</IT></OVL><SUB>o</SUB><IT>k</IT><SUB>o</SUB></SUP>) for <IT>t</IT> > <OVL><IT>t</IT></OVL><SUB>o</SUB>  (13b)

where [&Btilde;+]i(0,t) and [<A><AC>BH</AC><AC>∼</AC></A>]<SUB>i</SUB>(0,<IT>t</IT>) are [&Btilde;](x,t) and [<A><AC>BH</AC><AC>∼</AC></A>](x,t), respectively, at either the injection site or the capillary outlet, and <OVL><IT>t</IT></OVL>o is the mean transit time either between the injection site and the capillary inlet or between the capillary outlet and the sampling site. The values of ko for TBOH were 2.73 × 10-2 and 2.45 × 10-2 s-1 in PSS-albumin and PSS-dextran, respectively, estimated from the data in Fig. 7 as indicated in the APPENDIX.


View larger version (17K):
[in this window]
[in a new window]
 
Fig. 7.   TBO+ vs. time in 5% dextran and 5% albumin solutions when respective deoxygenated TBOH solutions were rapidly equilibrated with atmospheric oxygen.

Estimation of Model Parameters for Specific Experiments

TBO+ injections in lungs perfused with PSS-albumin. In the lungs perfused with PSS-albumin, the outflow concentrations of TBO+ and total TBO (TBOH + TBO+) were measured at different sampling sites as described in METHODS. The first step in the model interpretation of these data was to reconstruct the TBO+ curve measured by the photodetector as it would have appeared had it been measured at the fraction collector sampling site where the total TBO was measured. This was carried out as follows. The photodetector concentration curves acquired at 25 Hz for TBO+ and FITC were converted to discretely sampled curves with the same sampling time interval as the sample collector time interval using a infinite impulse response filter and decimator. The transfer function from the photodetector sampling site to the fraction collector was obtained by numerically deconvolving (13) the photodetector FITC tubing curve and the 125I fraction collector tubing curve. The resulting transfer function was then convolved with the discretized photodetector TBO+ curve to estimate the TBO+ curve at the sample collector site.

For a given bolus injection, the model was fit to the measured total TBO and the reconstructed TBO+ data as follows. Equations 9 and 10 were first solved with the initial (t = 0) conditions [&Btilde;+](x,0) = [<A><AC>BH</AC><AC>∼</AC></A>(x,0) = 0 and boundary conditions [<A><AC>BH</AC><AC>∼</AC></A>](0,t) = 0 and [&Btilde;+](0,t) = Cin(t)q/F. The solution of Eqs. 9 and 10 under these conditions provides the concentrations of TBO+ and total TBO at the capillary outlet. To account for any autoxidation of TBOH that might occur between capillary outlet and photodetector sampling site, the capillary outlet was subjected to Eq. 13, a and b, for <OVL><IT>t</IT></OVL>o, where in this case, [&Btilde;+]i(0,t) and [<A><AC>BH</AC><AC>∼</AC></A>]<SUB>i</SUB>(0,<IT>t</IT>) are the concentration versus time curves for [&Btilde;+] and [<A><AC>BH</AC><AC>∼</AC></A>], respectively, at the capillary outlet, and <OVL><IT>t</IT></OVL>o is the mean transit time between the capillary outlet and the photodetector. For these data, the identifiable model parameters are QF, QT, Kred, and Kseq, with ko fixed at the value estimated from the data in Fig. 7.

TBO+ and MB+ injections in lungs perfused with PSS-dextran. For TBO+ and MB+ injections in lungs perfused with PSS-dextran, because [P] = 0, alpha 1 = 1, QF = Qf, and Kred = Ekr1 and Eqs. 9 and 10 simplify to
<FR><NU>∂[B<SUP>+</SUP>]</NU><DE>∂<IT>t</IT></DE></FR> + <IT>W</IT> <FENCE><FR><NU>Q<SUB>c</SUB></NU><DE>Q<SUB>c</SUB> + Q<SUB>F</SUB></DE></FR></FENCE> <FR><NU>∂[B<SUP>+</SUP>]</NU><DE>∂<IT>x</IT></DE></FR> = − <FR><NU>[B<SUP>+</SUP>] <IT>K</IT><SUB>red</SUB></NU><DE>Q<SUB>c</SUB> + Q<SUB>F</SUB></DE></FR> (14)
The model fit to the TBO+ or MB+ data from each injection was obtained by solving Eq. 14 with the initial condition [B+](x,0) = 0 and boundary condition [B+](0,t) = Cin(t)q/F. For these data, the identifiable model parameters are QF and Kred.

TBOH injections in lungs perfused with PSS-albumin or PSS-dextran. For the TBOH injection data, the model was fit to the total TBO data from each injection by solving Eqs. 9 and 10 with the initial conditions [&Btilde;+](x,0) = [<A><AC>BH</AC><AC>∼</AC></A>](x,0) = 0 and boundary conditions
[<A><AC>BH</AC><AC>∼</AC></A>](0,<IT>t</IT>) = <FR><NU>C<SUB>in</SUB>(<IT>t</IT>) <IT>q</IT></NU><DE>F</DE></FR> <IT>e</IT><SUP>−<IT>k</IT><SUB>o</SUB><OVL><IT>t</IT></OVL><SUB>o</SUB></SUP>
and
[<A><AC>B</AC><AC>˜</AC></A><SUP>+</SUP>](0,<IT>t</IT>) = <FR><NU>C<SUB>in</SUB>(<IT>t</IT>) <IT>q</IT></NU><DE>F</DE></FR> (1 − <IT>e</IT><SUP>−<IT>k</IT><SUB>o</SUB><OVL><IT>t</IT></OVL><SUB>o</SUB></SUP>)
The above boundary conditions account for the autoxidation of the TBOH during transit time <OVL><IT>t</IT></OVL>o upstream from the capillary inlet. For the PSS-albumin data, the identifiable model parameters are QT and Kseq, with QF and Kred fixed to the mean values estimated from the TBO+ injections given in Table 3. For the PSS-dextran data, the identifiable model parameter is QT, with QF and Kred fixed to the mean values estimated from the TBO+ injections given in Table 4 and Kseq fixed to the value given in Table 5.

                              
View this table:
[in this window]
[in a new window]
 
Table 3.   Model parameters for TBO+ injections in lungs perfused with PSS-albumin


                              
View this table:
[in this window]
[in a new window]
 
Table 4.   Model parameters for TBO+ injections in lungs perfused with PSS-dextran


                              
View this table:
[in this window]
[in a new window]
 
Table 5.   Model parameters for TBOH injections in lungs perfused with PSS-albumin

Model Results

Examples of the model fits to the data are included in Figs. 1-5, and Tables 3-6 summarize the model parameter values. Examples of the sensitivity functions theta iS(t) (6) obtained for the optimized parameters are given in Figs. 8 and 9 for TBO+ and total TBO after TBO+ injections into the PSS-albumin perfusate and for TBOH after TBOH injections into the PSS-albumin perfusate, respectively. For the ith model parameter, theta i, Si(t) = partial C(t)/partial theta i, where C(t) is the calculated test indicator effluent concentration. Si(t) was approximated by the change in C(t) resulting from a 1% change in theta i divided by the change in theta i (6). When Si(t) is multiplied by the value of the parameter estimate, theta i, the relative amplitude of this function for each respective parameter provides an indication of the relative contribution of the parameter to the model fit to the data at a given time. Comparison of the respective shapes of their respective theta iSi(t) reveals the degree to which any pair of parameters is correlated. Interpretation with respect to the individual parameters is discussed below.

                              
View this table:
[in this window]
[in a new window]
 
Table 6.   Model parameters for MB+ injections in lungs perfused with PSS-dextran



View larger version (29K):
[in this window]
[in a new window]
 
Fig. 8.   Sensitivity functions [theta S(t)] for 4 parameters, Kred, QF, QT, and Kseq, from model fit to TBO+ and total TBO (TBOH + TBO+) concentration data obtained after injection of TBO+ into a lung perfused with PSS-albumin. Reference indicator concentration curves [CR(t)] are also provided for timing perspective. See text for definition of model parameters.



View larger version (14K):
[in this window]
[in a new window]
 
Fig. 9.   Sensitivity functions [theta S(t)] for the 2 parameters QT and Kseq from model fit to TBOH concentration data obtained after injection of TBOH into a lung perfused with PSS-albumin.


    DISCUSSION
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES
APPENDIX

Both TBO+ and MB+ were extracted from the PSS-dextran perfusate on passage through the lungs. However, the TBO+ extractions were considerably greater than those of MB+. The TBO+ extractions were >70% even at flow rates above those reported for the normal resting pulmonary flow rate for rabbits in this size range (276-540 ml/min; Ref. 38). This is in contrast to MB+, for which the extraction was only ~15% at only 300 ml/min. The role of reduction in the extraction of MB+ during its passage through the lungs was demonstrated previously (10). The fact that TBO+ was also reduced during passage through the lungs was revealed only when albumin was present in the perfusate. This is because the highly lipophilic TBOH has such a high solubility in the lung tissue that, with no protein in the perfusate, the effluent contained little trace of any TBOH that might have been formed. The presence of TBOH in the effluent when the lungs were perfused with the albumin-containing perfusate resulted from its affinity for the albumin, which offset its affinity for the tissue. The present study does not identify the location of the TBO+ reduction. However, a previous study, with a TBO-containing polymer that could not enter the cells, established the ability of endothelial cells to reduce TBO+ on the cell surface (12). Therefore, it appears likely that this surface reduction accounts for at least part of the disappearance of the relatively hydrophilic TBO+ within the lungs.

Another key observation was that the extractions of both MB+ and TBO+ were inversely proportional to flow rate to the extent that variations in the reduction rate parameter, Kred, with flow rate were small over the fairly wide range of flow rates studied. The implication of this observation is that the reduction rate, rather than the rate of convective supply of the indicators, dominates the extraction. One might visualize the reduction reaction as a barrier that must be traversed for TBO or MB to enter the tissue. In the terminology commonly used in the indicator dilution field, their extractions tend toward "barrier-limited" rather than "flow-limited" (23) behavior. The isopleths in Fig. 6 are for the most parsimonious nested version of the model developed herein, commonly referred to as the Crone-Renkin model (7), for which the solution for Kred is Kred = -F × ln(1 - extraction). For this model, barrier-limited extraction as a function of flow rate produces curves parallel to these isopleths. Thus the extent to which the TBO+ and MB+ extraction data parallel these isopleths may be thought of as a preliminary indication of the extent to which the extractions correspond to barrier- or reaction-limited behavior. The alternative of flow-limited behavior would produce horizontal lines in Fig. 6. As discussed below, the inclusion of additional detail in the model can explain most of the deviation from the Crone-Renkin model prediction observable in Fig. 6.

The model developed herein to interpret the test indicator probe disposition represents a simplistic view of the overall processing of the indicators by the lung. Steps have been left out of the stoichiometric equations that could be added if warranted by experimental manipulations of the relevant variables in future studies. The attempt was to accommodate aspects of the tissue disposition of TBO and MB predicted by previous studies under the present range of experimental conditions. Even so, the model is fairly complex, having several parameters that were either input, i.e., obtained from sources separate from the data from an individual bolus, or estimated by fitting the model to the bolus data. The fitted parameters include terms that are not separately identifiable without a more complex set of experimental conditions. Furthermore, all of the model parameters do not contribute equally to the fits to the data from each type of experiment performed or from the same type of experiment but at different flow rates. In fact, the effects of some phenomena represented in the model are not detectable in all experimental conditions. To help put this in perspective, the sensitivity functions for the model fits to the TBO+ and total TBO data from the PSS-albumin experiments are shown in Fig. 8. The sensitivity functions reveal the extent, and the time epoch, to which the optimized model parameters make their contributions to the fit of the model to the data (2, 6). Thus the dominant role of Kred in the fits to the TBO+ data, especially near the peak of the reference indicator curve can be seen in Fig. 8. QF has most of its influence at early and late times with respect to the reference indicator passage and little influence near the peak of the reference indicator. The very different shapes of the sensitivity functions for these two parameters also imply their relatively independent contributions to the model fit. This is expressed quantitatively for the examples in Fig. 8 by the fact that the correlations (3, 26) between QF and Kred were <0.34 at all flow rates studied. The physical counterpart of QF is not certain. It may include a volume into which at least some TBO+ diffuses before its reduction, possibly inside the cell. However, a compound that binds to plasma protein may likewise have an affinity for molecules on the cell surface, presumably including plasma albumin itself associated with the endothelial surface or glycocalyx (2, 43, 44).

The contribution of QT to the fit to the TBO+ data is small and barely detectable on the same scale as theta S(t) for Kred and QF at the high end of the range of flow rates studied. On the other hand, when the model is fit to the total TBO data, the contribution of QT is substantial and the impact of Kseq, which does not contribute to the TBO+ fit, is evident.

Kred, QF, and QT are the dominant parameters in the PSS-albumin experiment. The sensitivity functions for the MB+ and TBO+ experiments in the PSS-dextran-perfused lungs are not shown because they are qualitatively similar to those for TBO+ in PSS-albumin sensitivity functions, but without QT. Thus, for the PSS-dextran experiments, wherein the returning flux of TBOH is insignificant, Kred and QF are the controlling parameters.

The sensitivity functions for the TBOH injected into PSS-albumin reveal the importance of QT and Kseq (Fig. 9). TBOH is a lipophilic amine compound, and except for the small effects due to autoxidation, the model for its tissue disposition is the same as the general model developed for lipophilic amines in a previous study (3). The reasons for variations in the model parameters QF, QT, and Kseq for the MB+, TBO+, and TBOH injections with flow rate have been discussed previously (2, 3). In essence, the phenomenon occurs because two or more reversible processes having different dissociation rate constants, which are not sufficiently different for the processes to be distinguishable in the data obtained at a given flow rate, become lumped into one process in an MID model. When the flow rate is changed the relative contribution of each process to the data also changes, and the lumped parameter(s) representing these processes will be flow rate dependent. One conclusion with regard to the present study is that the reversible processes represented by QF and QT include distributed dissociation rate constants not explicitly represented in the model or evident in the data from a single flow rate (3). Likewise, the term sequestration applied to Kseq means that the dissociation rates are long relative to the capillary transit time, not that the processes are necessarily irreversible.

The overall tissue disposition of TBO and MB apparently involves the reduction, and nonspecific binding of the oxidized forms, and the partitioning of the reduced forms between tissue and perfusate volumes as represented in the model. However, it has also been demonstrated in endothelial cells in culture that TBOH can be reoxidized and sequestered within the cells (33). The sequestration apparently occurs because the oxidized forms of the thiazines are strongly polar, cationic (14), relatively hydrophilic compounds and thus relatively membrane impermeant. Thus the oxidized forms cannot so readily escape their intracellular site(s) of oxidation, and partitioning in accordance with membrane potentials of the sequestering intracellular organelle(s) may play a role as well (27). This sequestration contributes relatively little to the disposition of the TBO within the time frame of the bolus transient in normal lungs, as evidenced by the small value of Kseq (<3 ml/s) in comparison to Kred (>17 ml/s). However, preliminary results suggest that the sequestration rate may be increased in the lungs of rats adapted to hyperoxia, to the extent that it can have a substantial influence even during the bolus transient (39). Thus Kseq may take on greater significance in conditions in which intracellular reactions that reoxidize the reduced thiazine compounds may be affected.

General aspects of the sensitivity to the input parameters controlling bolus dispersion and capillary transit times have been discussed at length previously (4, 29). The main reason for the larger values of Kred obtained from the present model, reported in Tables 3, 4, and 6, in comparison to the Crone-Renkin model predictions, which can be read in Fig. 6, is the effect of capillary perfusion heterogeneity not included in the Crone-Renkin model. One of the important reasons for parameterizing the data using the particular modeling approach described herein is the potential for preventing a change in perfusion heterogeneity from being aliased by a change in tissue parameters (4). In the present study, the capillary transport function was assumed to be typical of the normal rabbit lung (5). However, in future studies of abnormal conditions, in which that assumption would not necessarily hold, inclusion of a flow-limited indicator such as labeled water in the bolus would provide the necessary information to account for a change in perfusion heterogeneity as discussed previously (4, 42).

Another input parameter was the TBOH autoxidation rate (ko), which was included to evaluate the assumption that the TBO+ data measured on-line can represent the TBO+ exiting the lung capillaries without contamination by autoxidation of effluent TBOH. The autoxidation turned out to be a minimal contributor to the on-line signal as can been seen in Fig. 10, in which a model simulation of the effluent TBO+, the effluent TBOH, and the TBO+ resulting from autoxidation of effluent TBOH is shown for the lowest flow rate in the range studied (i.e., the flow rate providing the longest time for autoxidation to occur).


View larger version (16K):
[in this window]
[in a new window]
 
Fig. 10.   Simulated effluent concentration vs. time curves after injection of TBO+ in a lung perfused with PSS-albumin. Curves are for TBO+ that passed through lungs without being reduced, TBOH generated by reduction of TBO+, and TBO+ resulting from autoxidation within perfusate of TBOH formed in lungs.

As indicated above, it was the presence of the albumin in the perfusate that revealed the reduction of TBO+. Because even without albumin binding the MB+ extraction was low, we did not pursue the consequences of including albumin in the perfusate for MB uptake. On the other hand, TBO uptake was high enough to provide an adequate window for examining the influence of albumin binding. As expected from the fact that TBO+ associates with plasma albumin, the presence of albumin decreased the rate of TBO+ disappearance from the perfusate. However, the decrease was not as great as predicted from the TBO+-BSA binding measured by ultrafiltration, under the assumption that only the free TBO participates in the uptake process. Under the assumption that the BSA-bound TBO+ is not reducible and using the average values of Kred from the PSS-albumin and PSS-dextran experiments, the predicted albumin binding would be only 39%, instead of the 59% obtained by ultrafiltration. This difference between the two is a common observation for the extraction of protein-bound substrates in the liver (45, 48) and other organs (25, 41). Several explanations have been proposed, including the possibility that albumin-bound ligand or substrate can interact directly with the cell surface process. Thus, for the TBO interaction with the pulmonary endothelium, it is conceivable that binding to BSA slows but does not prevent TBO reduction. This would be consistent with the observation that TBO covalently bound to an acrylamide polymer is reduced on the surface of pulmonary arterial endothelial cells in culture (12). It would also be consistent with the notion that endogenous electron acceptors such as the quinonoid forms of antioxidants such as tocopherylquinone and ubiquinone associated with plasma proteins such as low-density lipoproteins might be reducible without dissociating (11).

A similar analysis can be carried out on the data from the TBOH injections, from which albumin binding can be predicted by comparing the QT obtained in the PSS-dextran (average value 483 ml) and PSS-albumin experiments (Table 5; Ref. 3). The TBOH-BSA binding predicted on this basis is ~91%, which is higher than the 79% obtained by ultrafiltration, but the prediction is fairly sensitive to accuracy in the TBOH concentration data from the PSS-dextran experiment, in which the signal is quite small (Fig. 5).

This study was carried out in anticipation of the use of these thiazine compounds to determine the effects of various interventions directed at determining the influence of physiological and/or pathophysiological stresses on these redox functions of the pulmonary endothelium within intact lungs. The results suggest that both MB and TBO have potential utility in this regard. However, the low Kred for MB suggests that MB would be useful mainly for detecting substantial increases in Kred, because any decrease would be barely detectable at normal flow rates. TBO may be more versatile because of the higher extraction. For studies in isolated lungs, the presence of albumin in the perfusate would be required to evaluate the role of the sequestration process in the overall disposition of TBO+. The use of TBOH itself as test indicator might also be useful in this regard. The model should provide a tool for quantitative comparisons in such experiments, and it can also be used to provide a sense of the robustness of the quantification under a given set of experimental conditions.


    APPENDIX
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES
APPENDIX

Estimation of TBOH Autoxidation Rate From Test Tube Experiments

For autoxidation of BH (and for BHP in PSS-albumin) within the perfusate solutions in the test tube, the temporal variation in the concentrations of the oxidized (B+) and reduced (BH) forms is described by the following
<FR><NU>d[BH]</NU><DE>d<IT>t</IT></DE></FR> = [BHP]<IT>k</IT><SUB>−2</SUB> − <IT>k</IT><SUB>2</SUB>[BH][P] − <IT>k</IT><SUB>o</SUB>[BH] (A1)

<FR><NU>d[BHP]</NU><DE>d<IT>t</IT></DE></FR> = − [BHP]<IT>k</IT><SUB>−2</SUB> + <IT>k</IT><SUB>2</SUB>[BH][P] − <IT>k</IT><SUB>o</SUB>[BHP] (A2)

<FR><NU>d[B<SUP>+</SUP>]</NU><DE>d<IT>t</IT></DE></FR> = [B<SUP>+</SUP>P]<IT>k</IT><SUB>−1</SUB> − <IT>k</IT><SUB>1</SUB>[B<SUP>+</SUP>][P] + <IT>k</IT><SUB>o</SUB>[BH] (A3)

<FR><NU>d[B<SUP>+</SUP>P]</NU><DE>d<IT>t</IT></DE></FR> = − [B<SUP>+</SUP>P]<IT>k</IT><SUB>−1</SUB> + <IT>k</IT><SUB>1</SUB>[B<SUP>+</SUP>][P] + <IT>k</IT><SUB>o</SUB>[BHP] (A4)
Under the assumption that when albumin is present the equilibration is rapid between the free and albumin-bound species, Eqs. A1-A4 reduce to
<FR><NU>d[<A><AC>BH</AC><AC>∼</AC></A>]</NU><DE>d<IT>t</IT></DE></FR> = −<IT>k</IT><SUB>o</SUB>[<A><AC>BH</AC><AC>∼</AC></A>] (A5)

<FR><NU>d[<A><AC>B</AC><AC>˜</AC></A><SUP>+</SUP>]</NU><DE>d<IT>t</IT></DE></FR> = <IT>k</IT><SUB>o</SUB>[<A><AC>BH</AC><AC>∼</AC></A>] (A6)
For the initial conditions (t = 0) in the test tube experiments [<A><AC>BH</AC><AC>∼</AC></A>](0) = [<A><AC>BH</AC><AC>∼</AC></A><SUB>0</SUB>] and [&Btilde;](0) = [<A><AC>B</AC><AC>˜</AC></A><SUP>+</SUP><SUB>0</SUB>], the solution for [&Btilde;+](t) is
[<A><AC>B</AC><AC>˜</AC></A><SUP>+</SUP>](<IT>t</IT>) = [<A><AC>B</AC><AC>˜</AC></A><SUP>+</SUP><SUB>0</SUB>] + [<A><AC>BH</AC><AC>∼</AC></A><SUB>0</SUB>] (1 − <IT>e</IT><SUP>−<IT>k</IT><SUB>o</SUB><IT>t</IT></SUP>) (A7)
To estimate a value for ko, Eq. A7 was fit to the test tube [&Btilde;+](t) autoxidation data shown in Fig. 7.


    ACKNOWLEDGEMENTS

This study was supported by National Heart, Lung, and Blood Institute Grant NH-24349, the Falk Trust, the Whitaker Foundation, and the Department of Veterans Affairs.


    FOOTNOTES

The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. §1734 solely to indicate this fact.

Address for reprint requests and other correspondence: C. A. Dawson, Research Service 151, Zablocki VA Med. Ctr., 5000 W. National Ave., Milwaukee, WI 53295-1000 (E-mail: dawsonc{at}vms.csd.mu.edu).

Received 24 May 1999; accepted in final form 20 July 1999.


    REFERENCES
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES
APPENDIX

1.   Al-mehdi, A., G. Zhao, C. Dodia, K. Tozawa, K. Costa, V. Muzykantov, C. Ross, F. Blecha, M. Dinauer, and A. B. Fisher. Endothelial NADPH oxidase as the source of oxidants in lungs exposed to ischemia or high K+. Circ. Res. 83: 730-737, 1998[Abstract/Free Full Text].

2.   Audi, S. H., C. A. Dawson, J. H. Linehan, G. S. Krenz, S. B. Ahlf, and D. L. Roerig. An interpretation of 14C-urea and 14C-primidone extraction in isolated rabbit lungs. Ann. Biomed. Eng. 24: 337-351, 1996[Web of Science][Medline].

3.   Audi, S. H., C. A. Dawson, J. H. Linehan, G. S. Krenz, S. B. Ahlf, and D. L. Roerig. Pulmonary disposition of lipophilic amine compounds in the isolated perfused rabbit lung. J. Appl. Physiol. 84: 516-530, 1998[Abstract/Free Full Text].

4.   Audi, S. H., J. H. Linehan, G. S. Krenz, and C. A. Dawson. Accounting for the heterogeneity of capillary transit times in modeling multiple indicator dilution data. Ann. Biomed. Eng. 26: 914-930, 1998[Web of Science][Medline].

5.   Audi, S. H., J. H. Linehan, G. S. Krenz, C. A. Dawson, S. B. Ahlf, and D. L. Roerig. Estimation of the pulmonary transport function in isolated rabbit lungs. J. Appl. Physiol. 78: 1004-1014, 1995[Abstract/Free Full Text].

6.   Bassingthwaighte, J. B., and M. Chaloupka. Sensitivity functions in the estimation of parameters of cellular exchange. Federation Proc. 43: 181-184, 1984[Medline].

7.   Bassingthwaighte, J. B., and C. A. Goresky. Modeling in the analysis of solute and water exchange in the microvasculature. In: Handbook of Physiology. The Cardiovascular System. Microcirculation. Bethesda, MD: Am. Physiol. Soc, 1984, sect. 2, vol. IV, pt. 1, chapt. 13, p. 549-626.

8.   Bassingthwaighte, J. B., C. A. Goresky, and J. H. Linehan. Modeling in the analysis of the processes of uptake and metabolism in the whole organ. In: Whole Organ Approaches to Cellular Metabolism, edited by J. B. Bassingthwaighte, C. A. Goresky, and J. H. Linehan. New York: Springer, 1998, p. 3-27.

9.   Bassingthwaighte, J. B., K. Kroll, L. M. Schwartz, G. M. Raymond, and R. B. King. Strategies for uncovering the kinetics of nucleoside transport and metabolism in capillary endothelial cells. In: Whole Organ Approaches to Cellular Metabolism, edited by J. B. Bassingthwaighte, C. A. Goresky, and J. H. Linehan. New York: Springer, 1998, p. 163-188.

10.   Bongard, R. D., G. S. Krenz, J. H. Linehan, D. L. Roerig, M. P. Merker, J. L. Widell, and C. A. Dawson. Reduction and accumulation of methylene blue by the lung. J. Appl. Physiol. 77: 1480-1491, 1994[Abstract/Free Full Text].

11.   Bongard, R. D., M. P. Merker, J. M. Daum, and C. A. Dawson. Quinone reduction by endothelial cells: potential mechanism for regulating redox status of low density lipoproteins (LDL) (Abstract). FASEB J. 13: A185, 1999.

12.   Bongard, R. D., M. P. Merker, R. Shundo, Y. Okamoto, D. L. Roerig, J. H. Linehan, and C. A. Dawson. Reduction of thiazine dyes by bovine pulmonary arterial endothelial cells in culture. Am. J. Physiol. Lung Cell. Mol. Physiol. 269: L78-L84, 1995[Abstract/Free Full Text].

13.   Bronikowski, T. A., C. A. Dawson, and J. H. Linehan. Model-free deconvolution techniques for estimating vascular transport functions. Int. J. Biomed. Comput. 14: 411-429, 1983[Medline].

14.   Clark, W. M., B. Cohen, and H. D. Gibbs. Studies on oxidization and reduction. VIII. Methylene blue. Public Health Rep. 40: 1131-1201, 1925.

15.   Crane, F. L., I. L. Sun, R. Barr, and H. Low. Electron and proton transport across plasma membrane. J. Bioenerg. Biomembr. 23: 773-803, 1991[Web of Science][Medline].

16.   Dawson, C. A., C. W. Christensen, D. A. Rickaby, J. H. Linehan, and M. R. Johnston. Lung damage and pulmonary uptake of serotonin in intact dogs. J. Appl. Physiol. 58: 1761-1766, 1985[Abstract/Free Full Text].

17.   Dawson, C. A., and J. H. Linehan. Biogenic amines. In: Lung Biology in Health and Disease, edited by D. Massaro. New York: Dekker, 1988, vol. 41, p. 1091-1139. (Lung Cell Biology Ser.)

18.   Dawson, C. A., D. L. Roerig, and J. H. Linehan. Evaluation of endothelial injury in the human lung. Clin. Chest Med. 10: 13-24, 1989[Web of Science][Medline].

19.   DeLeo, F. R., and M. T. Quinn. Assembly of the phagocyte NADPH oxidase: molecular interaction of oxidase proteins. J. Leukoc. Biol. 60: 677-691, 1996[Abstract].

20.   Dupuis, J., C. A. Goresky, and A. Fournier. Pulmonary clearance of circulating endothelin-1 in dogs in vivo: exclusive role of ETB receptors. J. Appl. Physiol. 81: 1510-1515, 1996[Abstract/Free Full Text].

21.   Gillis, C. N., L. H. Cronau, S. Mandel, and G. L. Hammond. Indicator dilution measurement of 5-hydroxytryptamine clearance by human lung. J. Appl. Physiol. 46: 1178-1183, 1979[Abstract/Free Full Text].

22.   Giulivi, C., and E. Cadenas. Extracellular activation of fluorinated aziridinylbenzoquinone in HT29 cells EPR studies. Chem. Biol. Interact. 113: 191-204, 1998[Web of Science][Medline].

23.   Goresky, C. A., W. H. Ziegler, and G. G. Bach. Capillary exchange modeling: barrier-limited and flow-limited distribution. Circ. Res. 27: 739-764, 1970[Abstract/Free Full Text].

24.   Harris, T. R. The transport of small molecules across the microvascular barrier as measure of permeability and functioning exchange area in the normal and acutely injured lung. In: Whole Organ Approaches to Cellular Metabolism, edited by J. B. Bassingthwaighte, C. A. Goresky, and J. H. Linehan. New York: Springer, 1998, p. 439-454.

25.   Hutter, J. F., H. M. Piper, and P. G. Spieckermann. Kinetic analysis of myocardial fatty acid oxidation suggesting an albumin receptor mediated uptake process. J. Mol. Cell. Cardiol. 16: 219-226, 1984[Web of Science][Medline].

26.   Jacquez, J. A., and T. Perry. Parameter estimation: local identifiability of parameters. Am. J. Physiol. Endocrinol. Metab. 258: E727-E736, 1990[Abstract/Free Full Text].

27.   Johnson, L. V., M. L. Walsh, and L. B. Chen. Localization of mitochondria in living cells with rhodamine 123. Proc. Natl. Acad. Sci. USA 77: 990-994, 1980[Abstract/Free Full Text].

28.   Kaul, N., J. Choi, and H. J. Forman. Transmembrane redox signaling activates NF-kappa B in macrophages. Free Radic. Biol. Med. 24: 202-207, 1998[Web of Science][Medline].

29.   King, R. B., G. M. Raymond, and J. B. Bassingthwaighte. Modeling blood flow heterogeneity. Ann. Biomed. Eng. 24: 352-372, 1996[Web of Science][Medline].

30.   Linehan, J. H., S. H. Audi, and C. A. Dawson. The uptake and metabolism of substrates by endothelium in the lung. In: Whole Organ Approaches to Cellular Metabolism, edited by J. B. Bassingthwaighte, C. A. Goresky, and J. H. Linehan. New York: Springer, 1998, p. 427-437.

31.   Marquardt, D. W. An algorithm for least-squares estimation of nonlinear parameters. J. Soc. Ind. Appl. Math. 11: 431-441, 1963.

32.   May, J. M., Z.-C. Qu, and R. R. Whitesell. Ascorbic acid recycling enhances the antioxidant reserve of human erythrocytes. Biochemistry 34: 12721-12728, 1995[Medline].

33.   Merker, M. P., R. D. Bongard, J. H. Linehan, Y. Okamoto, D. Yvprachticky, B. M. Brantmeier, D. L. Roerig, and C. A. Dawson. Pulmonary endothelial thiazine uptake: separation of cell surface reduction from intracellular reoxidation. Am. J. Physiol. Lung Cell. Mol. Physiol. 272: L673-L680, 1997[Abstract/Free Full Text].

34.   Merker, M. P., L. E. Olson, R. D. Bongard, M. K. Patel, J. H. Linehan, and C. A. Dawson. Ascorbate-mediated transplasma membrane electron transport in pulmonary arterial endothelial cells. Am. J. Physiol. Lung Cell. Mol. Physiol. 274: L685-L693, 1998[Abstract/Free Full Text].

35.   Mohazzab-H, K. M., P. M. Kaminski, R. Agarwal, and M. S. Wolin. Potential role of a membrane-bound NADH oxidoreductase in nitric oxide release and arterial relaxation to nitroprusside. Circ. Res. 84: 220-228, 1999[Abstract/Free Full Text].

36.   Morré, D. J., L.-Y. Wu, and D. M. Morré. Response of a cell surface NADH oxidase to the antitumor sulfonylurea N-(4-methylphenylsulfonyl)-N'-(4-chlorophenylurea) (LY181984) modulated by redox. Biochim. Biophys. Acta 1369: 185-192, 1998[Medline].

37.   Navarro, F., P. Navas, J. R. Burgess, R. I. Bello, R. De Cabo, A. Arroyo, and J. M. Villalba. Vitamin E and selenium deficiency induces expression of the ubiquinone-dependent antioxidant system at the plasma membrane. FASEB J. 12: 1665-1673, 1998[Abstract/Free Full Text].

38.   Neutze, J. M., F. Wyler, and A. M. Rudolph. Use of radioactive microspheres to assess distribution of cardiac output in rabbits. Am. J. Physiol. 215: 486-495, 1968.

39.   Olson, L. E., S. H. Audi, N. Lin, W. Lin, and C. A. Dawson. Adaptation to hyperoxia results in pulmonary sequestration of toluidine blue O (TBO+) (Abstract). FASEB J. 13: A184, 1999.

40.   Olson, L. E., M. P. Merker, R. D. Bongard, B. M. Brantmeier, S. H. Audi, J. H. Linehan, and C. A. Dawson. Kinetics of plasma membrane electron transport in a pulmonary endothelial cell-column. Ann. Biomed. Eng. 26: 117-127, 1998[Web of Science][Medline].

41.   Pardridge, W. M., and E. M. Landaw. Tracer kinetic model of blood brain barrier transport of plasma protein-bound ligands. Empiric testing of the free hormone hypothesis. J. Clin. Invest. 74: 746-752, 1984.

42.   Roerig, D. L., S. H. Audi, J. H. Linehan, G. S. Krenz, S. B. Ahlf, W. Lin, and C. A. Dawson. Detection of changes in lung tissue properties with multiple indicator dilution. J. Appl. Physiol. 86: 1866-1880, 1999[Abstract/Free Full Text].

43.   Schneeberger, E. E., R. D. Lynch, and B. A. Neary. Interaction of native and chemically modified albumin with pulmonary microvascular endothelium. Am. J. Physiol. Lung Cell. Mol. Physiol. 258: L89-L98, 1990[Abstract/Free Full Text].

44.   Schnitzer, J. E., W. W. Carley, and G. E. Palade. Specific albumin binding to microvascular endothelium in culture. Am. J. Physiol. Heart Circ. Physiol. 254: H425-H437, 1988[Abstract/Free Full Text].

45.   Schwab, A. J., and C. A. Goresky. Hepatic uptake of protein-bound ligands: effect of an unstirred Disse space. Am. J. Physiol. Gastrointest. Liver Physiol. 270: G869-G880, 1996[Abstract/Free Full Text].

46.   Thorstensen, K., and I. Romslo. Uptake of iron from transferin by isolated rat hepatocytes. A redox-mediated plasma membrane process. J. Biol. Chem. 263: 8844-8850, 1988[Abstract/Free Full Text].

47.   Villalba, J. M., F. Navarro, F. Cordoba, A. Serrano, A. Arroyo, F. L. Crane, and P. Navas. Coenzyme Q reductase from liver plasma membrane: purification and role in trans-plasma-membrane electron transport. Proc. Natl. Acad. Sci. USA 92: 4887-4891, 1995[Abstract/Free Full Text].

48.   Weisiger, R. A., S. Pond, and L. Bass. Hepatic uptake of protein-bound ligands: extended sinusoidal perfusion model. Am. J. Physiol. Gastrointest. Liver Physiol. 261: G872-G884, 1991[Abstract/Free Full Text].

49.   Wiedemann, H. P., M. A. Matthay, and C. N. Gillis. Pulmonary endothelial cell injury and altered lung metabolic function. Clin. Chest Med. 11: 723-736, 1990[Web of Science][Medline].

50.   Zulueta, J. J., F.-S. Yu, I. A. Hertig, and V. J. Thannickal. Release of hydrogen peroxide in response to hypoxia-reoxygenation: role of an NAD(P)H oxidase-like enzyme in endothelial cell plasma membrane. Am. J. Respir. Cell Mol. Biol. 12: 41-49, 1995[Abstract].

51.   Zurbriggen, R., and J. L. Dreyer. An NADH-diaphorase is located at the cell plasma membrane in a mouse neuroblastoma cell line NB41A3. Biochem. Biophys. Acta 1183: 513-520, 1994[Medline].


Am J Physiol Heart Circ Physiol 278(1):H137-H150



This article has been cited by other articles:


Home page
J. Appl. Physiol.Home page
S. H. Audi, M. P. Merker, G. S. Krenz, T. Ahuja, D. L. Roerig, and R. D. Bongard
Coenzyme Q1 redox metabolism during passage through the rat pulmonary circulation and the effect of hyperoxia
J Appl Physiol, October 1, 2008; 105(4): 1114 - 1126.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Physiol. Lung Cell. Mol. Physiol.Home page
S. H. Audi, R. D. Bongard, G. S. Krenz, D. A. Rickaby, S. T. Haworth, J. Eisenhauer, D. L. Roerig, and M. P. Merker
Effect of chronic hyperoxic exposure on duroquinone reduction in adult rat lungs
Am J Physiol Lung Cell Mol Physiol, November 1, 2005; 289(5): L788 - L797.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Physiol. Lung Cell. Mol. Physiol.Home page
S. H. Audi, R. D. Bongard, C. A. Dawson, D. Siegel, D. L. Roerig, and M. P. Merker
Duroquinone reduction during passage through the pulmonary circulation
Am J Physiol Lung Cell Mol Physiol, November 1, 2003; 285(5): L1116 - L1131.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Physiol. Lung Cell. Mol. Physiol.Home page
S. H. Audi, C. A. Dawson, S. B. Ahlf, and D. L. Roerig
Oxygen dependency of monoamine oxidase activity in the intact lung
Am J Physiol Lung Cell Mol Physiol, October 1, 2001; 281(4): L969 - L981.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Physiol. Lung Cell. Mol. Physiol.Home page
S. H. Audi, R. D. Bongard, Y. Okamoto, M. P. Merker, D. L. Roerig, and C. A. Dawson
Pulmonary reduction of an intravascular redox polymer
Am J Physiol Lung Cell Mol Physiol, June 1, 2001; 280(6): L1290 - L1299.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Web of Science (10)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Audi, S. H.
Right arrow Articles by Dawson, C. A.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Audi, S. H.
Right arrow Articles by Dawson, C. A.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
Visit Other APS Journals Online