Transit time dispersion in pulmonary and systemic circulation: effects of cardiac output and solute diffusivity

Histogram-based comparison between dynamic and static lung perfused blood volume images using dual energy CT

OBJECTIVE

The purpose of this study was to compare the results of a histogram-based analysis of static and dynamic lung perfused blood volume (LPBV) images.

METHODS

Sixty-five patients (mean age: 61.3 years, 36 male) underwent dynamic and static LPBV for evaluation of pulmonary vascular diseases (n = 11), lung carcinoma (n = 27) or pulmonary thromboembolism (PTE: n = 27). Seven sets of dynamic sequential scans were performed at the pulmonary trunk using dual-energy technique before the static LPBV scan. The image of lung parenchyma that showed the greatest mean attenuation in dynamic series was defined as the peak dynamic LPBV image. The differences and correlations in the mean attenuation, image noise, signal-to-noise ratio (SNR), contrast-to-noise ratio (CNR), histogram skewness and histogram kurtosis were evaluated according to the type of disease in static and dynamic LPBV images.

RESULTS

Static LPBV images showed significantly larger mean attenuation (Rt:24.2, Lt: 24.2), SNR (Rt:2.31, Lt:2.30), and CNR (Rt:2.40, Lt:2.39), and smaller kurtosis values (Rt:1.06, Lt:0.61) values in comparison to dynamic LPBV images (p < 0.001); however, with the exception of kurtosis of the left lung (r = 0.17), these values were well-corrected with that of the dynamic LPBV images in these values (r = 0.4-0.77, p ≤ 0.001) without kurtosis of left lung (r = 0.17) in all patients. The histogram kurtosis of static LPBV image showed a good correlation with that of dynamic LPBV (r = 0.41-0.77, p < 0.05), especially in patients with PTE.

CONCLUSION

In patients with PTE, the static LPBV image valueswere well correlated with the peak dynamic LPBV images which demonstrated pulmonary artery-dominant flow.

Beta blockade increases pulmonary and systemic transit time heterogeneity: evaluation based on indocyanine green kinetics in healthy volunteers

Knowledge of factors influencing the heterogeneity of blood transit times is important in cardiovascular physiology. The aim of the study was to investigate the effect of beta-adrenergic blockade on blood transit time dispersion in awake, anxious volunteers. Recirculatory modelling of the disposition of intravascular markers using parametric forms for transit time distributions, such as the inverse Gaussian distribution, provides the opportunity to estimate the systemic and pulmonary transit time dispersion in vivo. The latter is determined by the flow heterogeneity in the microcirculatory network. Using this approach, we have analysed indocyanine green (ICG) disposition data obtained in four subjects by frequent early arterial blood sampling before and after beta-adrenergic blockade by propranolol. Propranolol decreased cardiac output from 9·3 ± 2·8 l min^-1 to 3·5 ± 0·47 l min^-1 (P<0·05). This reduction was accompanied by a 4·5 ± 0·6-fold and 2·1 ± 0·3-fold increase (P<0·001) in the relative dispersion (dimensionless variance) of blood transit times through the systemic and pulmonary circulation, respectively.

A physiologically based model of hepatic ICG clearance: interplay between sinusoidal uptake and biliary excretion

Although indocyanine green (ICG) has long been used for the assessment of liver function, the respective roles of sinusoidal uptake and canalicular excretion in determining hepatic ICG clearance remain unclear. Here this issue was addressed by incorporating a liver model into a minimal physiological model of ICG disposition that accounts of the early distribution phase after bolus injection. Arterial ICG concentration-time data from awake dogs under control conditions and from the same dogs while anesthetized with 3.5% isoflurane were subjected to population analysis. The results suggest that ICG elimination in dogs is uptake limited since it depends on hepatocellular uptake capacity and on biliary excretion but not on hepatic blood flow. Isoflurane caused a 63% reduction in cardiac output and a 33% decrease in the ICG biliary excretion rate constant (resulting in a 26% reduction in elimination clearance) while leaving unchanged the sinusoidal uptake rate. The terminal slope of the concentration-time curve, K, correlated significantly with elimination clearance. The model could be useful for assessing the functions of sinusoidal and canalicular ICG transporters.

Cardiac output and systemic transit time dispersion as determinants of circulatory mixing time: a simulation study

A new approach to characterize the kinetics of intravascular mixing process is presented. The mixing time, defined as the time required for achieving 95% homogeneity, is calculated by numerical simulations using a circulatory model applied to the intravascular marker indocyanine green (ICG). The results suggest that the mixing time is determined by cardiac output and the relative dispersion of transit time distribution across the systemic circulation, whereby the rate of mixing increases with increasing cardiac output and decreasing transit time dispersion, and vice versa. The estimation of plasma volume from simulated ICG dilution data using the backextrapolation method shows that slow mixing is accompanied by an overestimation of blood volume. This error may be negligible for mixing times of less than approximately 3 min but high in disease states characterized by low cardiac output and/or high transit time dispersion. In view of the role of transit time dispersion as determinant of intravascular mixing, it would be interesting to know more about the effect of disease states on systemic transit time dispersion.

A note on population analysis of dissolution-absorption models using the inverse Gaussian function

Because conventional absorption models often fail to describe plasma concentration-time profiles following oral administration, empirical input functions such as the inverse Gaussian function have been successfully used. The purpose of this note is to extend this model by adding a first-order absorption process and to demonstrate the application of population analysis using maximum likelihood estimation via the EM algorithm (implemented in ADAPT 5). In one example, the analysis of bioavailability data of an extended-release formulation, as well as the mean dissolution times estimated in vivo and in vitro with the use of the inverse Gaussian function, is well in accordance, suggesting that the inverse Gaussian function indeed accounts for the in vivo dissolution process. In the other example, the kinetics of trapidil in patients with liver disease, the absorption/dissolution parameters are characterized by a high interindividual variability. Adding a first-order absorption process to the inverse Gaussian function improved the fit in both cases.

Advanced pharmacokinetic models based on organ clearance, circulatory, and fractal concepts

Three advanced models of pharmacokinetics are described. In the first class are physiologically based pharmacokinetic models based on in vitro data on transport and metabolism. The information is translated as transporter and enzyme activities and their attendant heterogeneities into liver and intestine models. Second are circulatory models based on transit time distribution and plasma concentration time curves. The third are fractal models for nonhomogeneous systems and non-Fickian processes are presented. The usefulness of these pharmacokinetic models, with examples, is compared.

Circulatory transport and capillary-tissue exchange as determinants of the distribution kinetics of inulin and antipyrine in dog

A pharmacokinetic model was developed to estimate physiologically meaningful parameters of distribution kinetics from plasma concentration-time data. The model is based on simultaneously measured disposition curves of drug and vascular marker. Employing residence time distribution theory, a recirculatory model with two subsystems, the pulmonary and systemic circulation, was constructed. In addition to intravascular mixing, the axially distributed model of the systemic circulation accounts for transcapillary transport of solutes, quantified by permeability-surface area product (PS) and diffusional equilibration time. Parameters of ICG, inulin, and antipyrine were estimated from disposition data obtained in awake dogs under control conditions and during an isoproterenol infusion or moderate hypovolemia. Results suggest that distribution kinetics is (1) governed by extravascular diffusion and (2) its dependency on cardiac output decreases with increasing diffusional resistance. Hemorrhage decreased the effective PS of inulin. In conclusion, this novel mechanistic model effectively described both the permeability-limited distribution of inulin into interstitial fluid and the flow-limited distribution of antipyrine into total body water and might be useful for other drugs.

Residence Time Dispersion as a General Measure of Drug Distribution Kinetics: Estimation and Physiological Interpretation

PURPOSE

To evaluate distribution kinetics of drugs by the relative dispersion of disposition residence time and demonstrate its uses, interpretation and limitations.

MATERIALS AND METHODS

The relative dispersion was estimated from drug disposition data of inulin and digoxin fitted by three-exponential functions, and calculated from compartmental parameters published for fentanyl and alfentanil. An interpretation is given in terms of a lumped organs model and the distributional equilibration process in a noneliminating system.

RESULTS

As a measure of the deviation from mono-exponential disposition (one-compartment behavior), the relative dispersion provides information on the distribution kinetics of drugs, i.e., diffusion-limited distribution or slow tissue binding, without assuming a specific structural model. It also defines the total distribution clearance which has a clear physical meaning.

CONCLUSION

The residence time dispersion is a model-independent measure that can be used to characterize the distribution kinetics of drugs and to reveal the influence of disease states. It can be estimated with high precision from drug disposition data.

Mechanistic modeling of digoxin distribution kinetics incorporating slow tissue binding

This study aims to develop a mechanistic pharmacokinetic model that accounts for the kinetics of tissue binding in order to evaluate the effect of slow binding of digoxin to skeletal muscular Na(+)/K(+)-ATPase in humans. The approach is based on a minimal circulatory model with a systemic transit time density function that accounts for vascular mixing, transcapillary permeation and extravascular binding of the drug. The model parameters were estimated using previously published disposition data of digoxin in healthy volunteers and physiological distribution volumes taken from the literature. A time constant of the binding process of 34min was estimated indicating that receptor binding and not permeation clearance is the rate-limiting step of the distribution process. Model simulations suggest that up- or downregulation of sodium pumps, typically observed under physiological or pathophysiological conditions, could be detected with this method. The model allows a quantitative prediction of the effect of changes in skeletal muscular sodium pump activity on plasma levels of digoxin.