Indocyanine green kinetics characterize blood volume and flow distribution and their alteration by propranolol

Prediction of Survival After Partial Hepatectomy Using a Physiologically Based Pharmacokinetic Model of Indocyanine Green Liver Function Tests

The evaluation of hepatic function and functional capacity of the liver are essential tasks in hepatology as well as in hepatobiliary surgery. Indocyanine green (ICG) is a widely applied test compound that is used in clinical routine to evaluate hepatic function. Important questions for the functional evaluation with ICG in the context of hepatectomy are how liver disease such as cirrhosis alters ICG elimination, and if postoperative survival can be predicted from preoperative ICG measurements. Within this work a physiologically based pharmacokinetic (PBPK) model of ICG was developed and applied to the prediction of the effects of a liver resection under various degrees of cirrhosis. For the parametrization of the computational model and validation of model predictions a database of ICG pharmacokinetic data was established. The model was applied (i) to study the effect of liver cirrhosis and liver resection on ICG pharmacokinetics; and (ii) to evaluate the model-based prediction of postoperative ICG-R15 (retention ratio 15 min after administration) as a measure for postoperative outcome. Key results are the accurate prediction of changes in ICG pharmacokinetics caused by liver cirrhosis and postoperative changes of ICG-elimination after liver resection, as validated with a wide range of data sets. Based on the PBPK model, individual survival after liver resection could be classified, demonstrating its potential value as a clinical tool.

Advanced Assessment of Hemodynamic Parameters during Donor Care

Detailed evaluation of cardiovascular function by hemodynamic “profiling” during donor care may be necessary to select appropriate treatment for hypotension. Preload, systemic vascular resistance, heart rate, and contractility are all determinants of cardiac output and blood pressure. This review discusses the measurement, physiological limitations, and pharmacological manipulation of those variables during care of selected donors.

Maternal-fetal transfer of indocyanine green across the perfused human placenta

Indocyanine green (ICG) is an FDA-approved near-infrared imaging probe, given also to pregnant women. We aimed to characterize ICG's transplacental transfer using the ex-vivo perfusion model. Placentas were obtained from caesarean deliveries. Cotyledons were cannulated and dually perfused. ICG, 9.6μg/mL and antipyrine (50μg/mL) were added to the maternal circulation in the absence (n=4) or the presence of the organic anion transporting polypeptide (OATPs) inhibitor rifampin (10μg/mL; n=5) or the P-glycoprotein inhibitor valspodar (2μg/mL; n=3). ICG's maternal-to-fetal transfer was evaluated over 180min. The cumulative percent of ICG in the fetal reservoir was minor. When ICG transfer was normalized to that of antipyrine, it was lower in the presence of rifampin (a 41% decrease; p<0.05). Valspodar did not appear to modify the kinetics of ICG. ICG's transplacental transfer is minimal and is probably OATP-mediated. The placenta is an effective protective barrier to ICG's distribution into the fetus.

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.

Optimal back-extrapolation method for estimating plasma volume in humans using the indocyanine green dilution method

BACKGROUND

The indocyanine green dilution method is one of the methods available to estimate plasma volume, although some researchers have questioned the accuracy of this method.

METHODS

We developed a new, physiologically based mathematical model of indocyanine green kinetics that more accurately represents indocyanine green kinetics during the first few minutes postinjection than what is assumed when using the traditional mono-exponential back-extrapolation method. The mathematical model is used to develop an optimal back-extrapolation method for estimating plasma volume based on simulated indocyanine green kinetics obtained from the physiological model.

RESULTS

Results from a clinical study using the indocyanine green dilution method in 36 subjects with type 2 diabetes indicate that the estimated plasma volumes are considerably lower when using the traditional back-extrapolation method than when using the proposed back-extrapolation method (mean (standard deviation) plasma volume = 26.8 (5.4) mL/kg for the traditional method vs 35.1 (7.0) mL/kg for the proposed method). The results obtained using the proposed method are more consistent with previously reported plasma volume values.

CONCLUSIONS

Based on the more physiological representation of indocyanine green kinetics and greater consistency with previously reported plasma volume values, the new back-extrapolation method is proposed for use when estimating plasma volume using the indocyanine green dilution method.

Population pharmacokinetics of cyclosporine in transplant recipients

A number of classical pharmacokinetic studies have been conducted in transplant patients. However, they suffer from some limitations, for example, (1) the study design was limited to intense blood sampling in small groups of patients during a certain posttransplant period, (2) patient factors were evaluated one at a time to identify their association with the pharmacokinetic parameters, and (3) mean pharmacokinetic parameters often cannot be precisely estimated due to large intraindividual variability. Population pharmacokinetics provides a potential means of addressing these limitations and is a powerful tool to evaluate the magnitude and consistency of drug exposure. Population pharmacokinetic studies of cyclosporine focused solely on developing limited sampling strategies and Bayesian estimators to estimate drug exposure, have been summarized before, and are, therefore, not a subject of this review. The major focus of this review is to describe factors (demographic factors, hepatic and gastrointestinal functions, drug-drug interactions, genetic polymorphisms of drug metabolizing enzymes and transporters) that have been identified to contribute to the large portion of observed variability in the pharmacokinetics of cyclosporine in transplant patients. This review summarizes and interprets the conclusions as well as the nonlinear mixed-effects modeling methodologies used in such studies. A highly diversified collection of structural models, variability models, and covariate submodels have been evaluated and validated using internal or external validation methods. This review also highlights areas where additional research is warranted to improve the models since a portion of model variability still remains unexplained.

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.

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.

β-Adrenergic Blockade Affects Initial Drug Distribution Due to Decreased Cardiac Output and Altered Blood Flow Distribution

Beta-adrenergic receptor blockers decrease intravenous anesthetic dose requirements. The present study determined the effect of propranolol on indocyanine green and antipyrine disposition from the moment of rapid intravenous injection. Anti-pyrine is a physiological marker that distributes to a volume as large as total body water in a blood flow-dependent manner and is a pharmacokinetic surrogate for many lipophilic drugs, including intravenous anesthetics. Antipyrine and indocyanine green disposition were determined twice in five healthy adult males in this Institutional Review Board-approved study, once during propranolol infusion. After rapid indocyanine green and antipyrine injection, arterial blood samples were collected frequently for 2 min and less frequently thereafter. Plasma indocyanine green and antipyrine concentrations were measured by high-performance liquid chromatography. Indocyanine green and antipyrine disposition were characterized, using SAAM II, by a recirculatory pharmacokinetic model that describes drug disposition from the moment of injection. Parameters were compared using the paired t test. The disposition of indocyanine green demonstrated that propranolol decreased cardiac output at the expense of the fast peripheral (nonsplanchnic) intravascular circuit. The area under the antipyrine concentration versus time relationship was doubled for at least the first 3 min after injection due to both decreased cardiac output and maintenance of nondistributive blood flow at the expense of a two-thirds reduction of blood flow (intercompartmental clearance) to the rapidly equilibrating (fast, splanchnic) tissue volume. The increase in antipyrine area under the curve due to propranolol-induced alteration of initial antipyrine disposition could explain decreased intravenous anesthetic dose requirements in the presence of beta-adrenergic receptor blockade.

Population pharmacokinetic modeling of oral cyclosporin using NONMEM: comparison of absorption pharmacokinetic models and design of a Bayesian estimator

There have been very few population pharmacokinetic (PopPK) studies and Bayesian forecasting methods dealing with cyclosporin (CsA) so far, probably because of the difficulty of modeling the particular absorption profiles of CsA. The present study was conducted in stable renal transplant patients treated with Neoral and employed the NONMEM program. Its goals were (1) to develop a population pharmacokinetic model for CsA based on an Erlang frequency distribution (which describes asymmetric S-shaped absorption profiles) combined with a 2-compartment model; (2) to compare this model with models combining a time-lag parameter and either a zero-order or first-order rate constant and with a model based on a Weibull distribution; and (3) to develop a PK Bayesian estimator for full AUC estimation based on that "Erlang model." The PopPK model was developed in an index set of 70 patients, and then individual PK parameters and AUC were estimated in 10 other patients using Bayesian estimation. The "Erlang" model best described the data, with mean absorption time (MAT), apparent clearance (CL/F), and apparent volume of the central compartment (Vc/F) of 0.78 hours, 26.3 L/h, and 76 L, respectively (interindividual variability CV = 33, 30, and 48%). Bayesian estimation allowed accurate prediction of systemic exposure using only 3 samples collected at 0, 1, and 3 hours. Regression analysis found no significant difference between the predicted and observed concentrations (10 per patient), and AUC(0-12) were estimated with a nonsignificant bias (0.6 to 8.7%) and good precision (RMSE = 5.3%). In conclusion, the Erlang distribution best described CsA absorption profiles, and a Bayesian estimator developed using this model and a mixed-effect PK modeling program provided accurate estimates of CsA systemic exposure using only 3 blood samples.

A Physiologically Based Approach for the Estimation of Recirculatory Parameters

Indicator dilution studies are used to provide estimates for several physiological parameters such as cardiac output as well as intra- and extravascular volumes. This study introduces a novel technique for the estimation of recirculatory parameters. A mathematical model based on a dispersion-convection partial differential equation (PDE), derived from the fractal geometry of the vascular tree and the hydrodynamics of the blood flow, is used to describe the spatiotemporal profile of tracers in the circulatory system. Initially, the equation is fitted to concentration-time (C,t) data of a tracer to derive the parameter estimates of the model equation; in a subsequent step, these estimates along with appropriate changes of the parameters of the PDE are used to generate the early concentration-time profile of a hypothetical appropriate tracer without recirculation. Thus, the area under the concentration-time curve of the first passage of the tracer is calculated and used for the estimation of various physiological parameters, including cardiac output, miscellaneous partial blood volumes, and the corresponding mean transit times. The procedure was applied successfully to literature data of various tracers from humans and dogs.

In vivo continuous-wave optical breast imaging enhanced with Indocyanine Green

We investigate the uptake of a nontargeted contrast agent by breast tumors using a continuous wave diffuse optical tomography apparatus. The instrument operates in the near-infrared spectral window and employs 16 sources and 16 detectors to collect light in parallel on the surface of the tumor-bearing breast (coronal geometry). In our protocol an extrinsic contrast agent, Indocyanine Green (ICG), was injected by bolus. Three clinical scenarios with three different pathologies were investigated. A two-compartment model was used to analyze the pharmacokinetics of ICG and preprocess the data, and diffuse optical tomography was used for imaging. Localization and delineation of the tumor was achieved in good agreement with a priori information. Moreover, different dynamical features were observed for differing pathologies. The malignant cases exhibited slower rate constants (uptake and outflow) compared to healthy tissue. These results provide further evidence that in vivo pharmacokinetics of ICG in breast tumors may be a useful diagnostic tool for differentiation of benign and malignant pathologies.

Evaluation of the splanchnic circulation with indocyanine green pharmacokinetics in liver transplant patients

Although indocyanine green (ICG) can be used to estimate cardiac output (CO) and blood volume independently, a recirculatory multicompartmental ICG model enables description of these and additional intravascular factors. This model was used to describe the effect of end-stage liver disease (ESLD) on systemic and splanchnic hemodynamics in patients undergoing orthotopic liver transplantation. ICG disposition was determined during the dissection phase in six patients with ESLD undergoing orthotopic liver transplantation and six healthy adult living liver donors. After injecting ICG, plasma concentrations were obtained for approximately 10 to 12 minutes by noninvasive pulse dye densitometry. The recirculatory model characterizes three distinct intravascular circuits: lumped parallel fast (presumably nonsplanchnic circulation) and slow peripheral (splanchnic) circuits and a central circuit (central blood volume). Mean transit time (MTT) in the fast peripheral circuit was not different in patients with ESLD and controls. However, ESLD resulted in a significant decrease in MTT in the central (0.11 +/- 0.028 [SD] v 0.24 +/- 0.094 minutes in controls; P <.001) and slow peripheral circuit (0.67 +/- 0.41 v 1.37 +/- 0.37 minutes in controls; P <.001) because of increased flows to the central and slow peripheral circuits. These findings are consistent with the described hyperdynamic systemic and splanchnic circulations in patients with ESLD. In conclusion, the ICG model is able to derive estimates of not only blood volume and CO, but also splanchnic hemodynamics under different physiological conditions. This model can be a useful tool to evaluate the effect of pharmacological manipulation of splanchnic hemodynamics.