PKQuest: volatile solutes - application to enflurane,nitrous oxide,halothane, methoxyflurane and toluene pharmacokinetics

BACKGROUND: The application of physiologically based pharmacokinetic models (PBPK) to human studies has been limited by the lack of the detailed organ information that is required for this analysis. PKQuest is a new generic PBPK that is designed to avoid this problem by using a set of "standard human" default parameters that are applicable to most solutes. RESULTS: PKQuest is used to model the human pharmacokinetics of the volatile solutes. A "standard human" value for the lipid content of the blood and each organ (klip) was chosen. This set of klip and the oil/water partition coefficient then specifies the organ/blood partition for each organ. Using this approach, the pharmacokinetics of inert volatile solute is completely specified by just 2 parameters: the water/air and oil/water partition coefficients. The model predictions of PKQuest were in good agreement with the experimental data for the inert solutes enflurane and nitrous oxide and the metabolized solutes halothane and toluene. METHODS: The experimental data that was modeled was taken from previous publications. CONCLUSIONS: This approach greatly increases the predictive power of the PBPK. For inert volatile solutes the pharmacokinetics are determined just from the water/air and oil/water partition coefficient. Methoxyflurane cannot be modeled by this PBPK because the arterial and end tidal partial pressures are not equal (as assumed in the PBPK). This inequality results from the "washin-washout" artifact in the large airways that is established for solutes with large water/air partition coefficients.PKQuest and the worked examples are available on the web www.pkquest.com.     

Heterogeneity of human adipose blood flow

Background

The long time pharmacokinetics of highly lipid soluble compounds is dominated by blood-adipose tissue exchange and depends on the magnitude and heterogeneity of adipose blood flow. Because the adipose tissue is an infinite sink at short times (hours), the kinetics must be followed for days in order to determine if the adipose perfusion is heterogeneous. The purpose of this paper is to quantitate human adipose blood flow heterogeneity and determine its importance for human pharmacokinetics.

Methods

The heterogeneity was determined using a physiologically based pharmacokinetic model (PBPK) to describe the 6 day volatile anesthetic data previously published by Yasuda et. al. The analysis uses the freely available software PKQuest and incorporates perfusion-ventilation mismatch and time dependent parameters that varied from the anesthetized to the ambulatory period. This heterogeneous adipose perfusion PBPK model was then tested by applying it to the previously published cannabidiol data of Ohlsson et. al. and the cannabinol data of Johansson et. al.

Results

The volatile anesthetic kinetics at early times have only a weak dependence on adipose blood flow while at long times the pharmacokinetics are dominated by the adipose flow and are independent of muscle blood flow. At least 2 adipose compartments with different perfusion rates (0.074 and 0.014 l/kg/min) were needed to describe the anesthetic data. This heterogeneous adipose PBPK model also provided a good fit to the cannabinol data.

Conclusion

Human adipose blood flow is markedly heterogeneous, varying by at least 5 fold. This heterogeneity significantly influences the long time pharmacokinetics of the volatile anesthetics and tetrahydrocannabinol. In contrast, using this same PBPK model it can be shown that the long time pharmacokinetics of the persistent lipophilic compounds (dioxins, PCBs) do not depend on adipose blood flow. The ability of the same PBPK model to describe both the anesthetic and cannabinol kinetics provides direct qualitative evidence that their kinetics are flow limited and that there is no significant adipose tissue diffusion limitation.     

The use of a physiologically based pharmacokinetic model to evaluate deconvolution measurements of systemic absorption

Background

An unknown input function can be determined by deconvolution using the systemic bolus input function (r) determined using an experimental input of duration ranging from a few seconds to many minutes. The quantitative relation between the duration of the input and the accuracy of r is unknown. Although a large number of deconvolution procedures have been described, these routines are not available in a convenient software package.

Methods

Four deconvolution methods are implemented in a new, user-friendly software program (PKQuest, http://www.pkquest.com). Three of these methods are characterized by input parameters that are adjusted by the user to provide the "best" fit. A new approach is used to determine these parameters, based on the assumption that the input can be approximated by a gamma distribution. Deconvolution methodologies are evaluated using data generated from a physiologically based pharmacokinetic model (PBPK).

Results and Conclusions

The 11-compartment PBPK model is accurately described by either a 2 or 3-exponential function, depending on whether or not there is significant tissue binding. For an accurate estimate of r the first venous sample should be at or before the end of the constant infusion and a long (10 minute) constant infusion is preferable to a bolus injection. For noisy data, a gamma distribution deconvolution provides the best result if the input has the form of a gamma distribution. For other input functions, good results are obtained using deconvolution methods based on modeling the input with either a B-spline or uniform dense set of time points.

     

Influence of time of day on propofol pharmacokinetics and pharmacodynamics in rabbits

This study evaluates the administration time-of-day effects on propofol pharmacokinetics and sedative response in rabbits. Nine rabbits were sedated with 5?mg/kg propofol at three local clock times: 10:00, 16:00, and 22:00?h. Each rabbit served as its own control by being given a single infusion at the three different times of day on three separate occasions. Ten arterial blood samples were collected during each clock-time experiment for propofol assay. A two-compartment model was used to describe propofol pharmacokinetics, and the pedal withdrawal reflex was used as the sedation pharmacodynamic response. The categorical data comprising the presence or absence of pedal withdrawal reflex was described by a logistic model. The typical volume of the central compartment equaled 7.67?L and depended on rabbit body weight. The elimination rate constant depended on drug administration time; it was lowest at 10:00?h, highest at 16:00?h, and intermediate at 22:00?h. Delay of the anesthetic effect, with respect to plasma concentrations, was described by the effect compartment, with the rate constant for the distribution to the effector compartment equal to 0.335?min(-1). Drug concentration had a large effect on the probability of anesthesia. The degree of anesthesia was largest at 10:00?h, lowest at 16:00?h, and intermediate at 22:00?h. In summary, both the pharmacokinetics and pharmacodynamics of propofol in rabbits depended on administration time. The developed population approach may be used to assess chronopharmacokinetics and chronopharmacodynamics of medications in animals and humans.     

Influence of Time of Day on Propofol Pharmacokinetics and Pharmacodynamics in Rabbits

This study evaluates the administration time-of-day effects on propofol pharmacokinetics and sedative response in rabbits. Nine rabbits were sedated with 5?mg/kg propofol at three local clock times: 10:00, 16:00, and 22:00?h. Each rabbit served as its own control by being given a single infusion at the three different times of day on three separate occasions. Ten arterial blood samples were collected during each clock-time experiment for propofol assay. A two-compartment model was used to describe propofol pharmacokinetics, and the pedal withdrawal reflex was used as the sedation pharmacodynamic response. The categorical data comprising the presence or absence of pedal withdrawal reflex was described by a logistic model. The typical volume of the central compartment equaled 7.67?L and depended on rabbit body weight. The elimination rate constant depended on drug administration time; it was lowest at 10:00?h, highest at 16:00?h, and intermediate at 22:00?h. Delay of the anesthetic effect, with respect to plasma concentrations, was described by the effect compartment, with the rate constant for the distribution to the effector compartment equal to 0.335?min^?1. Drug concentration had a large effect on the probability of anesthesia. The degree of anesthesia was largest at 10:00?h, lowest at 16:00?h, and intermediate at 22:00?h. In summary, both the pharmacokinetics and pharmacodynamics of propofol in rabbits depended on administration time. The developed population approach may be used to assess chronopharmacokinetics and chronopharmacodynamics of medications in animals and humans. (Author correspondence: agnbienert@op.pl)     

Use of stable isotopically labeled tracers to measure very low density lipoprotein-triglyceride turnover

Tracer methods for VLDL-TG kinetics vary in their ability to account for the effect of tracer recycling, which can influence the calculation of VLDL-TG fractional catabolic rates (FCRs). We evaluated a novel approach, involving stable isotopically labeled glycerol or palmitate tracers in conjunction with compartmental modeling, for measuring VLDL-TG kinetics in normolipidemic human subjects. When administered as a bolus simultaneously, both tracers provided identical VLDL-TG FCRs when the data were analyzed by a compartmental model that accounted for hepatic lipid tracer recycling, but not by non-compartmental analysis. The model-derived FCR was greater than that determined using a non-compartmental approach, and was 2- to 3-fold higher than that usually reported by using a bolus of radioactive [3H]glycerol. When palmitate tracer was given as a constant infusion, VLDL-TG turnover appeared 5-fold slower, because tracer recycling through hepatic lipid pools could not be resolved with the infusion protocol. We conclude that accounting for tracer recycling, particularly the contribution of hepatic glycerolipid pools, is essential to accurately measure VLDL-TG kinetics, and that bolus injection of stable isotopically labeled glycerol or palmitate tracers in conjunction with compartmental modeling analysis offers a reliable approach for measuring VLDL-TG kinetics.     

Determination of kinetic parameters of apolipoprotein B metabolism using amino acids labeled with stable isotopes.

The use of amino acids labeled with stable isotopes represents a relatively new approach for determining kinetic parameters of apolipoprotein metabolism; thus, several aspects of experimental protocols need to be defined. The aims of the present study were to determine whether a) different amino acid tracers or b) different methods of tracer administration affected apolipoprotein (apo) B kinetic parameters obtained by multicompartmental modeling, and c) to compare very low density lipoprotein (VLDL)-apoB metabolic parameters determined by multicompartmental modeling with those estimated by linear regression or by monoexponential analysis. [1-13C]leucine and [15N]glycine were given either as bolus injections or as primed constant infusions. A bolus of one amino acid was administered simultaneously with a primed constant infusion (8 h) of the other amino acid into four healthy normolipidemic subjects (age 23.0 +/- 1.4 yr; BMI 20.9 +/- 0.9 kg.m-2). VLDL-, intermediate density lipoprotein (IDL)-, and low density lipoprotein (LDL)-apoB enrichments were followed over 110 h. For subsequent analysis these values were converted to tracer/tracee ratios. Using the multicompartmental model, the fractional catabolic rate (FCR) for VLDL-apoB was estimated to be 0.36 +/- 0.09 h-1 after the administration of the tracer as a primed constant infusion and 0.35 +/- 0.07 h-1 when the tracer was administered as a bolus. The values for VLDL-apoB production were 14.6 +/- 6.5 mg.kg-1.d-1 and 14.1 +/- 5.4 mg.kg-1.d-1, respectively. The corresponding values for LDL-apoB were 0.027 +/- 0.016 h-1 (0.026 +/- 0.018 h-1) for the FCR and 10.5 +/- 3.7 mg.kg-1.d-1 (10.4 +/- 3.8 mg.kg-1.d-1) for the production following administration of the tracer as a primed constant infusion and a bolus, respectively. Approximately 47% of VLDL-apoB ultimately reached the LDL fraction via the VLDL-IDL-LDL pathway. Thirty-five percent of LDL-apoB did not originate from this cascade pathway, but was shunted from a rapidly turning over VLDL compartment directly into the LDL fraction. While there was some variation between individuals, VLDL-apoB and LDL-apoB parameters derived from the bolus and the primed constant infusions showed no significant differences and were closely correlated. Metabolic parameters were also independent of the two amino acids tested. Although values for FCRs of VLDL-apoB obtained from linear regression (0.36 +/- 0.19 h-1) or monoexponential analysis (0.50 +/- 0.36 h-1) did not differ significantly from those obtained by the multicompartmental model, there was considerable variation and no significant correlation in a given individual.(ABSTRACT TRUNCATED AT 400 WORDS)     

Tissue distribution and physiologically based pharmacokinetics of antisense phosphorothioate oligonucleotide ISIS 1082 in rat

The aim of this study was to develop a whole body physiologically based model of the pharmacokinetics (PBPK) of the phosphorothioate oligonucleotide (PS-ODN) ISIS 1082 in vivo. Rats were administered an intravenous (i.v.) bolus dose of ISIS 1082 (10 mg/kg plus 3H tracer), and arterial blood and tissues were taken at specific times up to 72 hours. Radioactivity was measured in all samples. The parent compound was determined specifically in blood and tissues at 90 minutes and in liver and kidney also at 24 hours, using capillary gel electrophoresis (CGE). A whole body PBPK model was fitted to the combined blood and tissue radioactivity data using nonlinear regression analysis. CGE analysis indicated that the predominant species in plasma and all tissues is ISIS 1082, together with some n-1 and n-2 metabolites. Total radioactivity primarily reflects these species. The whole body model successfully described temporal events in all tissues. However, to adequately model the experimental data, all tissues had to be partitioned into vascular and extravascular spaces to accommodate the relatively slow distribution of ISIS 1082 out of blood because of a permeability rate limitation. ISIS 1082 distributes extensively into tissues, but the relative affinity varies enormously, being highest for kidney and liver and lowest for muscle and brain. A whole body PBPK model with a permeability rate limited tissue distribution was developed that adequately described events in both blood and tissue for an oligonucleotide. This model has the potential not only to characterize the events in individual tissues throughout the body for such compounds but also to scale across animal species, including human.     

In vivo simulation of human pharmacokinetics in the rabbit

The evaluation of drugs in vivo is often based on experimental models using small animals such as mice, rats and rabbits. However, these models could be improved to correspond more closely to the human situation if the pharmacokinetics of the drugs tested in animals were similar to that observed in humans. The use of a computer-controlled pump allowing an adequate flow of tobramycin and amikacin to be infused into rabbits enabled us to simulate the human pharmacokinetics of these antibiotics in vivo in this study. The function defining the rate of infusion required to perform the simulation of an intravenous bolus was first determined generally and symbolically for linear pharmacokinetic models independently from the number of compartments involved. The practical simulation of a decreasing monoexponential serum profile with a half-life of 2 h (one-compartment model for the human pharmacokinetics of aminoglycosides) was then studied for tobramycin and amikacin on the basis of a two-compartment model in the animal. The kinetics obtained had an apparent elimination half-life of 1.97 and 1.86 h, respectively. Linearity of the semilogarithmic regressions of the profiles obtained was quite sound. Finally, an a posteriori analysis of the pharmacokinetic model and its parameters is proposed on the basis of the results obtained after simulation.     

Dobutamine as selective beta(1)-adrenoceptor agonist in in vivo studies on human thermogenesis and lipid utilization.

The use of dobutamine as selective beta(1)-adrenoceptor agonist in in vivo studies on human thermogenesis and lipid utilization was investigated in 20 men. At 2.5, 5, and 10 microg x kg(-1) x min(-1), dobutamine induced significant increases in energy expenditure, lipid oxidation, and lipolysis. The beta(1)-adrenoceptor antagonist atenolol (bolus: 42.5 microg/kg, infusion: 1.02 microg x kg(-1) x min(-1)) blocked all dobutamine-induced effects on thermogenesis and lipid utilization. All parameters remained at levels comparable to those during saline infusion. The dose of atenolol used did not inhibit beta(2)-adrenoceptor-specific changes in energy expenditure, lipid oxidation, and lipolysis during salbutamol infusion (85 ng x kg(-1) x min(-1)). This indicates that atenolol was specific for beta(1)-adrenoceptors and did not camouflage concomitant beta(2)-adrenoceptor stimulation during dobutamine infusion. Therefore, we conclude that dobutamine can be used as a selective beta(1)-adrenoceptor agonist at dosages 相似文献   

Tissue expression profile of human neonatal Fc receptor (FcRn) in Tg32 transgenic mice

The neonatal Fc receptor (FcRn) is a homeostatic receptor responsible for prolonging immunoglobulin G (IgG) half-life by protecting it from lysosomal degradation and recycling it to systemic circulation. Tissue-specific FcRn expression is a critical parameter in physiologically-based pharmacokinetic (PBPK) modeling for translational pharmacokinetics of Fc-containing biotherapeutics. Using online peptide immuno-affinity chromatography coupled with high resolution mass spectrometry, we established a quantitative FcRn tissue protein expression profile in human FcRn (hFcRn) transgenic mice, Tg32 homozygous and hemizygous strains. The concentration of hFcRn across 14 tissues ranged from 3.5 to 111.2 pmole per gram of tissue. Our hFcRn quantification data from Tg32 mice will enable a more refined PBPK model to improve the accuracy of human PK predictions for Fc-containing biotherapeutics.     

PBPK models in risk assessment--A focus on chloroprene

Mathematical models are increasingly being used to simulate events in the exposure-response continuum, and to support quantitative predictions of risks to human health. Physiologically based pharmacokinetic (PBPK) models address that portion of the continuum from an external chemical exposure to an internal dose at a target site. Essential data needed to develop a PBPK model include values of key physiological parameters (e.g., tissue volumes, blood flow rates) and chemical specific parameters (rate of chemical absorption, distribution, metabolism, and elimination) for the species of interest. PBPK models are commonly used to: (1) predict concentrations of an internal dose over time at a target site following external exposure via different routes and/or durations; (2) predict human internal concentration at a target site based on animal data by accounting for toxicokinetic and physiological differences; and (3) estimate variability in the internal dose within a human population resulting from differences in individual pharmacokinetics. Himmelstein et al. [M.W. Himmelstein, S.C. Carpenter, P.M. Hinderliter, Kinetic modeling of beta-chloroprene metabolism. I. In vitro rates in liver and lung tissue fractions from mice, rats, hamsters, and humans, Toxicol. Sci. 79 (1) (2004) 18-27; M.W. Himmelstein, S.C. Carpenter, M.V. Evans, P.M. Hinderliter, E.M. Kenyon, Kinetic modeling of beta-chloroprene metabolism. II. The application of physiologically based modeling for cancer dose response analysis, Toxicol. Sci. 79 (1) (2004) 28-37] developed a PBPK model for chloroprene (2-chloro-1,3-butadiene; CD) that simulates chloroprene disposition in rats, mice, hamsters, or humans following an inhalation exposure. Values for the CD-PBPK model metabolic parameters were obtained from in vitro studies, and model simulations compared to data from in vivo gas uptake studies in rats, hamsters, and mice. The model estimate for total amount of metabolite in lung correlated better with rodent tumor incidence than did the external dose. Based on this PBPK model analytical approach, Himmelstein et al. [M.W. Himmelstein, S.C. Carpenter, M.V. Evans, P.M. Hinderliter, E.M. Kenyon, Kinetic modeling of beta-chloroprene metabolism. II. The application of physiologically based modeling for cancer dose response analysis, Toxicol. Sci. 79 (1) (2004) 28-37; M.W. Himmelstein, R. Leonard, R. Valentine, Kinetic modeling of beta-chloroprene metabolism: default and physiologically-based modeling approaches for cancer dose response, in: IISRP Symposium on Evaluation of Butadiene & Chloroprene Health Effects, September 21, 2005, TBD--reference in this proceedings issue of Chemical-Biological Interactions] propose that observed species differences in the lung tumor dose-response result from differences in CD metabolic rates. The CD-PBPK model has not yet been submitted to EPA for use in developing the IRIS assessment for chloroprene, but is sufficiently developed to be considered. The process that EPA uses to evaluate PBPK models is discussed, as well as potential applications for the CD-PBPK model in an IRIS assessment.     

Intragastric pH and Pharmacokinetics of Intravenous Ranitidine During Sinusoidal and Constant-Rate Infusions

Six patients with healed duodenal ulcer completed two treatment periods with continuous i.v. infusion ranitidine. A 25-mg i.v. bolus was followed by a constant infusion at 6.25 mg/h or a sinusoidal infusion with infusion rates ranging from 3.125 to 9.375 mg/h. The sinusoidal infusion rate was designed to match the previously observed circadian changes in basal acid secretion. The peak infusion rate occurred at 19:30 h. A pharmacokinetic method was designed to predict the resultant plasma concentrations of ranitidine. Intragastric pH and plasma ranitidine concentration data were fit to a cosine function to evaluate circadian and ultradian rhythms. Plasma concentrations during the sinusoidal infusion exhibited a circadian rhythm according to model predictions. Cosinor analyses of the mean ranitidine plasma concentration data showed a mesor concentration of 237 ng/mL and amplitude of 76 ng/mL (coefficient of determination [CD] = 0.98). The acrophase in plasma concentration occurred at 2223 h, a delay of approximately 2.9 hours from the peak in the infusion rate. The constant-rate infusion resulted in a mean plasma concentration of 222 ± 32 ng/mL. The 24-h mean intragastric pH values for the sinusoidal and constant regimens were 5.4 and 5.1, respectively (p = 0.170). The intragastric pH during the constant-rate infusion exhibited a significant circadian rhythm (CD = 0.52). The minimum pH (bathy-phase) occurred at 2031 h. No circadian rhythm was present during the sinusoidal-rate infusion (CD = 0.08). At the approximate time of the peak basal acid secretion, between 21:00 hours and midnight, the mean pH for the sinusoidal infusion was 5.77 versus 4.5 for the constant-rate infusion (p = 0.112). Sinusoidal infusions or alternate methods of increased doses at the times of peak acid output may improve around-the-clock control of intragastric pH.     

Prediction of Pharmacokinetic Profile of Valsartan in Humans Based on in vitro Uptake‐Transport Data

The aim of this study was to evaluate a physiologically based pharmacokinetic (PBPK) model for predicting PK profiles in humans based on a model refined in rats and humans in vitro uptake‐transport data using valsartan as a probe substrate. Valsartan is eliminated unchanged, mostly through biliary excretion, both in humans and rats. It was, therefore, chosen as model compound to predict in vivo elimination based on in vitro hepatic uptake‐transport data using a fully mechanistic PBPK model. Plated rat and human hepatocytes, and cell lines overexpressing human OATP1B1 and OATP1B3 were used for in vitro uptake experiments. A mechanistic two‐compartment model was used to derive the active and passive transport parameters, namely uptake Michaelis–Menten parameters (V_max and K_m,u) together with passive diffusion (P_dif). These transport parameters were then used as input in a whole body physiologically based pharmacokinetic (PBPK) model. The uptake rate of valsartan was higher for rat hepatocytes (K_m,u=28.4±3.7 μM , V_max=1320±180 pmol/mg/min, and P_dif =1.21±0.42 μl/mg/min) compared to human hepatocytes (K_m,u=44.4±14.6 μM , V_max=304±85 pmol/mg/min, and P_dif=0.724±0.271 μl/mg/min). OATP1B1 and ‐1B3 parameters were correlated to human hepatocyte data, using experimentally established relative activity factors (RAF). Resulting PBPK simulations were compared for plasma‐ (humans and rats) and bile‐ (rats) concentration–time profiles following iv bolus administration of valsartan. Plasma clearances (CL_P) for rats and humans were predicted within twofold relative to predictions based on respective in vitro data. The simulations were extended to simulate the impact of either OATP1B1 or ‐1B3 inhibition on plasma profile. The limited data set indicates that the mechanistic model allowed for accurate evaluation of in vitro transport data; and the resulting hepatic uptake transport kinetic parameters enabled the prediction of in vivo PK profiles and plasma clearances, using PBPK modelling. Moreover, the interspecies difference in elimination rate observed in vivo was correctly reflected in the transport parameters determined in vitro.     

PKQuest: a general physiologically based pharmacokinetic model. Introduction and application to propranolol

Background  

A "physiologically based pharmacokinetic" (PBPK) approach uses a realistic model of the animal to describe the pharmacokinetics. Previous PBPKs have been designed for specific solutes, required specification of a large number of parameters and have not been designed for general use.     

Propofol Pharmacokinetics and Estimation of Fetal Propofol Exposure during Mid-Gestational Fetal Surgery: A Maternal-Fetal Sheep Model

Background

Measuring fetal drug concentrations is extremely difficult in humans. We conducted a study in pregnant sheep to simultaneously describe maternal and fetal concentrations of propofol, a common intravenous anesthetic agent used in humans. Compared to inhalational anesthesia, propofol supplemented anesthesia lowered the dose of desflurane required to provide adequate uterine relaxation during open fetal surgery. This resulted in better intraoperative fetal cardiac outcome. This study describes maternal and fetal propofol pharmacokinetics (PK) using a chronically instrumented maternal-fetal sheep model.

Methods

Fetal and maternal blood samples were simultaneously collected from eight mid-gestational pregnant ewes during general anesthesia with propofol, remifentanil and desflurane. Nonlinear mixed-effects modeling was performed by using NONMEM software. Total body weight, gestational age and hemodynamic parameters were tested in the covariate analysis. The final model was validated by bootstrapping and visual predictive check.

Results

A total of 160 propofol samples were collected. A 2-compartment maternal PK model with a third fetal compartment appropriately described the data. Mean population parameter estimates for maternal propofol clearance and central volume of distribution were 4.17 L/min and 37.7 L, respectively, in a typical ewe with a median heart rate of 135 beats/min. Increase in maternal heart rate significantly correlated with increase in propofol clearance. The estimated population maternal-fetal inter-compartment clearance was 0.0138 L/min and the volume of distribution of propofol in the fetus was 0.144 L. Fetal propofol clearance was found to be almost negligible compared to maternal clearance and could not be robustly estimated.

Conclusions

For the first time, a maternal-fetal PK model of propofol in pregnant ewes was successfully developed. This study narrows the gap in our knowledge in maternal-fetal PK model in human. Our study confirms that maternal heart rate has an important influence on the pharmacokinetics of propofol during pregnancy. Much lower propofol concentration in the fetus compared to maternal concentrations explain limited placental transfer in in-vivo paired model, and less direct fetal cardiac depression we observed earlier with propofol supplemented inhalational anesthesia compared to higher dose inhalational anesthesia in humans and sheep.     

Effects of propofol intravenous injection bolus on the left ventricular function and the myocardial beta-adrenoceptor in rats

Propofol bolus injection has been reported to influence cardiovascular functions. However, the detailed mechanism underlying this action has not been elucidated. This study was designed to investigate the effects of propofol i.v. bolus on the left ventricular function, the myocardial beta-adrenoceptor (beta-AR) binding-site density (Bmax) and Kd (apparent dissociation constant) in a 30-minute period. One hundred and four male Wistar rats were randomly divided into four groups: group C (control group), group I (intralipid group), group P1 (5 mg/kg propofol) and group P2 (10 mg/kg propofol). The results showed a significant downregulation of HR, LVSP, +dp/dtmax and -dp/dtmax in both groups P1 and P2 (especially after bolus injection in 7 min) than those of group C (P < 0.05), whereas no significant difference was found between the P1 and P2 groups (P > 0.05). Likely, Bmax was remarkably upregulated in both groups P1 and P2 (P < 0.05, vs. groups C and I), and there was no significant difference between these two groups (P > 0.05). Of note, the Kd value in group P2 (10 mg/kg propofol) was found dramatically increased in 30 min than that in the low-dose propofol-treated group (group P1) as well as in groups C and I (P < 0.05). In conclusion, these results indicate that intravenous injection of propofol bolus can inhibit the cardiac function partially via upregulation of Bmax and downregulation of the beta-AR affinity at higher-dose injection of propofol bolus.     

Bayesian Population Physiologically-Based Pharmacokinetic (PBPK) Approach for a Physiologically Realistic Characterization of Interindividual Variability in Clinically Relevant Populations

Interindividual variability in anatomical and physiological properties results in significant differences in drug pharmacokinetics. The consideration of such pharmacokinetic variability supports optimal drug efficacy and safety for each single individual, e.g. by identification of individual-specific dosings. One clear objective in clinical drug development is therefore a thorough characterization of the physiological sources of interindividual variability. In this work, we present a Bayesian population physiologically-based pharmacokinetic (PBPK) approach for the mechanistically and physiologically realistic identification of interindividual variability. The consideration of a generic and highly detailed mechanistic PBPK model structure enables the integration of large amounts of prior physiological knowledge, which is then updated with new experimental data in a Bayesian framework. A covariate model integrates known relationships of physiological parameters to age, gender and body height. We further provide a framework for estimation of the a posteriori parameter dependency structure at the population level. The approach is demonstrated considering a cohort of healthy individuals and theophylline as an application example. The variability and co-variability of physiological parameters are specified within the population; respectively. Significant correlations are identified between population parameters and are applied for individual- and population-specific visual predictive checks of the pharmacokinetic behavior, which leads to improved results compared to present population approaches. In the future, the integration of a generic PBPK model into an hierarchical approach allows for extrapolations to other populations or drugs, while the Bayesian paradigm allows for an iterative application of the approach and thereby a continuous updating of physiological knowledge with new data. This will facilitate decision making e.g. from preclinical to clinical development or extrapolation of PK behavior from healthy to clinically significant populations.