Human physiologically based pharmacokinetic model for propofol

Utility of life stage-specific chemical risk assessments based on New Approach Methodologies (NAMs)

The safety assessments for chemicals targeted for use or expected to be exposed to specific life stages, including infancy, childhood, pregnancy and lactation, and geriatrics, need to account for extrapolation of data from healthy adults to these populations to assess their human health risk. However, often adequate and relevant toxicity or pharmacokinetic (PK) data of chemicals in specific life stages are not available. For such chemicals, New Approach Methodologies (NAMs), such as physiologically based pharmacokinetic (PBPK) modeling, biologically based dose response (BBDR) modeling, in vitro to in vivo extrapolation (IVIVE), etc. can be used to understand the variability of exposure and effects of chemicals in specific life stages and assess their associated risk. A life stage specific PBPK model incorporates the physiological and biochemical changes associated with each life stage and simulates their impact on the absorption, distribution, metabolism, and elimination (ADME) of these chemicals. In our review, we summarize the parameterization of life stage models based on New Approach Methodologies (NAMs) and discuss case studies that highlight the utility of a life stage based PBPK modeling for risk assessment. In addition, we discuss the utility of artificial intelligence (AI)/machine learning (ML) and other computational models, such as those based on in vitro data, as tools for estimation of relevant physiological or physicochemical parameters and selection of model. We also discuss existing gaps in the available toxicological datasets and current challenges that need to be overcome to expand the utility of NAMs for life stage-specific chemical risk assessment.

Comparing Perioperative Outcomes of Total Intravenous Anesthesia (TIVA) With Volatile Anesthesia in Patients With Obesity: A Systematic Review

In this systematic review, the perioperative outcomes of total intravenous anesthesia (TIVA) and volatile anesthesia were compared in obese adults (BMI ≥ 30 kg/m²) undergoing elective surgery. The review analyzed data from 12 randomized-controlled trials involving 935 patients, sourced from PubMed/MEDLINE (Medical Literature Analysis and Retrieval System Online), Cochrane, Scopus, and Web of Science databases. The focus was on intraoperative vital signs, emergence time, postoperative nausea and vomiting (PONV), duration of post-anesthesia care unit (PACU) stay, and ICU admission rates. Findings showed that TIVA (using propofol) might reduce PONV, but there were no significant differences in other outcomes compared to volatile anesthesia (with desflurane as the most common agent). The review highlights the need for more research, especially comparing sevoflurane with TIVA, to establish clear clinical guidelines for anesthesia in obese patients.

A Physiologically Based Pharmacokinetic Approach to Recommend an Individual Dose of Tacrolimus in Adult Heart Transplant Recipients

Tacrolimus is the principal immunosuppressive drug which is administered after heart transplantation. Managing tacrolimus therapy is challenging due to a narrow therapeutic index and wide pharmacokinetic (PK) variability. We aimed to establish a physiologically based pharmacokinetic (PBPK) model of tacrolimus in adult heart transplant recipients to optimize dose regimens in clinical practice. A 15-compartment full-PBPK model (Simbiology® Simulator, version 5.8.2) was developed using clinical observations from 115 heart transplant recipients. This study detected 20 genotypes associated with tacrolimus metabolism. CYP3A5*3 (rs776746), CYP3A4*18B (rs2242480), and IL-10 G-1082A (rs1800896) were identified as significant genetic covariates in tacrolimus pharmacokinetics. The PBPK model was evaluated using goodness-of-fit (GOF) and external evaluation. The predicted peak blood concentration (Cmax) and area under the drug concentration-time curve (AUC) were all within a two-fold value of the observations (fold error of 0.68-1.22 for Cmax and 0.72-1.16 for AUC). The patients with the CYP3A5*3/*3 genotype had a 1.60-fold increase in predicted AUC compared to the patients with the CYP3A5*1 allele, and the ratio of the AUC with voriconazole to alone was 5.80 when using the PBPK model. Based on the simulation results, the tacrolimus dosing regimen after heart transplantation was optimized. This is the first PBPK model used to predict the PK of tacrolimus in adult heart transplant recipients, and it can serve as a starting point for research on immunosuppressive drug therapy in heart transplant patients.

Evaluating the trophic transfer of PCBs from fish to humans: Insights from a synergism of environmental monitoring and physiologically-based pharmacokinetic modeling

Accumulation of polychlorinated biphenyls (PCBs) within fish tissues has prompted many states to issue consumption advisories. In Pennsylvania such advisories suggest one meal per month for most game species harvested from Lake Erie; however, these advisories do not account for the emergent properties of regional PCB mixtures, and the downstream accumulation of PCB congeners into human tissues is poorly documented. This study aimed to demonstrate the utility of pairing environmental monitoring with pharmacokinetic modeling for the purpose of estimating dietary PCB exposure in humans. We qualified and quantified the PCB congeners present in the filets of five Lake Erie fish species and used these data to estimate exposure under consumption scenarios that matched or exceeded the advisories. Physiologically-based pharmacokinetic (PBPK) modeling was then employed to predict PCB accumulation within seven tissue compartments of a hypothetical man and woman over 10 years. Twenty-one congeners were detected between the five fish species at concentrations ranging from 56.0 to 411.7 ng/g. Predicted accumulation in human tissues varied based on tissue type, the species consumed, biological sex, and fish-consumption rate. Notably, steady-state concentrations were higher in fatty tissue compartments ("Fat" and "Liver") and across all tissues in women compared to men. This study serves as a preliminary blueprint for generating predictions of site-specific and tissue-specific exposure through the integration of environmental monitoring and pharmacokinetic modeling. Although the details may vary across applications, this simple approach could complement traditional exposure assessments for vulnerable communities in the Great Lakes region that continue to suffer from legacy contamination.

The Impact of Low Cardiac Output on Propofol Pharmacokinetics across Age Groups-An Investigation Using Physiologically Based Pharmacokinetic Modelling

BACKGROUND

pathophysiological changes such as low cardiac output (LCO) impact pharmacokinetics, but its extent may be different throughout pediatrics compared to adults. Physiologically based pharmacokinetic (PBPK) modelling enables further exploration.

METHODS

A validated propofol model was used to simulate the impact of LCO on propofol clearance across age groups using the PBPK platform, Simcyp^® (version 19). The hepatic and renal extraction ratio of propofol was then determined in all age groups. Subsequently, manual infusion dose explorations were conducted under LCO conditions, targeting a 3 µg/mL (80-125%) propofol concentration range.

RESULTS

Both hepatic and renal extraction ratios increased from neonates, infants, children to adolescents and adults. The relative change in clearance following CO reductions increased with age, with the least impact of LCO in neonates. The predicted concentration remained within the 3 µg/mL (80-125%) range under normal CO and LCO (up to 30%) conditions in all age groups. When CO was reduced by 40-50%, a dose reduction of 15% is warranted in neonates, infants and children, and 25% in adolescents and adults.

CONCLUSIONS

PBPK-driven, the impact of reduced CO on propofol clearance is predicted to be age-dependent, and proportionally greater in adults. Consequently, age group-specific dose reductions for propofol are required in LCO conditions.

A "Smart" Biosensor-Enabled Intravascular Catheter and Platform for Dynamic Delivery of Propofol to "Close the Loop" for Total Intravenous Anesthesia

BACKGROUND

Target-controlled infusion anesthesia is used worldwide to provide user-defined, stable, blood concentrations of propofol for sedation and anesthesia. The drug infusion is controlled by a microprocessor that uses population-based pharmacokinetic data and patient biometrics to estimate the required infusion rate to replace losses from the blood compartment due to drug distribution and metabolism. The objective of the research was to develop and validate a method to detect and quantify propofol levels in the blood, to improve the safety of propofol use, and to demonstrate a pathway for regulatory approval for its use in the USA.

METHODS

We conceptualized and prototyped a novel "smart" biosensor-enabled intravenous catheter capable of quantifying propofol at physiologic levels in the blood, in real time. The clinical embodiment of the platform is comprised of a "smart" biosensor-enabled catheter prototype, a signal generation/detection readout display, and a driving electronics software. The biosensor was validated in vitro using a variety of electrochemical methods in both static and flow systems with biofluids, including blood.

RESULTS

We present data demonstrating the experimental detection and quantification of propofol at sub-micromolar concentrations using this biosensor and method. Detection of the drug is rapid and stable with negligible biofouling due to the sensor coating. It shows a linear correlation with mass spectroscopy methods. An intuitive graphical user interface was developed to: (1) detect and quantify the propofol sensor signal, (2) determine the difference between targeted and actual propofol concentration, (3) communicate the variance in real time, and (4) use the output of the controller to drive drug delivery from an in-line syringe pump. The automated delivery and maintenance of propofol levels was demonstrated in a modeled benchtop "patient" applying the known pharmacokinetics of the drug using published algorithms.

CONCLUSIONS

We present a proof-of-concept and in vitro validation of accurate electrochemical quantification of propofol directly from the blood and the design and prototyping of a "smart," indwelling, biosensor-enabled catheter and demonstrate feedback hardware and software architecture permitting accurate measurement of propofol in blood in real time. The controller platform is shown to permit autonomous, "closed-loop" delivery of the drug and maintenance of user-defined propofol levels in a dynamic flow model.

Constructing a control-ready model of EEG signal during general anesthesia in humans

Significant effort toward the automation of general anesthesia has been made in the past decade. One open challenge is in the development of control-ready patient models for closed-loop anesthesia delivery. Standard depth-of-anesthesia tracking does not readily capture inter-individual differences in response to anesthetics, especially those due to age, and does not aim to predict a relationship between a control input (infused anesthetic dose) and system state (commonly, a function of electroencephalography (EEG) signal). In this work, we developed a control-ready patient model for closed-loop propofol-induced anesthesia using data recorded during a clinical study of EEG during general anesthesia in ten healthy volunteers. We used principal component analysis to identify the low-dimensional state-space in which EEG signal evolves during anesthesia delivery. We parameterized the response of the EEG signal to changes in propofol target-site concentration using logistic models. We note that inter-individual differences in anesthetic sensitivity may be captured by varying a constant cofactor of the predicted effect-site concentration. We linked the EEG dose-response with the control input using a pharmacokinetic model. Finally, we present a simple nonlinear model predictive control in silico demonstration of how such a closed-loop system would work.

Desflurane anesthesia compared with total intravenous anesthesia on anesthesia-controlled operating room time in ambulatory surgery following strabotomy: a randomized controlled study

BACKGROUND

Ophthalmic ambulatory surgery is preferred to be performed under general anesthesia either by total intravenous anesthesia (TIVA) or by inhalational anesthesia to increase the patient comfort. However, anesthesia-controlled time (ACT) can cause increased non-operative operating room (OR) time which may adversely affect the ORs efficiency. This study was aimed to compare the ACT of desflurane with that of propofol-remifentanil in strabismus ambulatory surgery.

METHODS

From November 2016 to December 2017, a total of 200 strabismus patients (aged 18-60 years old, and scheduled for elective ambulatory surgery at Zhongshan Ophthalmic Center) were randomly assigned to receive either propofol-based TIVA (group TIVA) or desflurane anesthesia (group DES) for maintenance of anesthesia. The primary outcome was the extubation time. Secondary outcomes included surgical time, anesthetic time, OR exit time, and Phase I and II recovery time. The intraoperative incidences of hypotension, bradycardia and oculocardiac reflex (OCR), and the incidences of any post-operative complications were recorded. Mann-Whitney U test and Chi-square or Fisher exact tests were used to compare the two groups.

RESULTS

We found that the extubation time (5.5 [3.9-7.0] vs. 9.7 [8.5-11.4] min, P < 0.001) and the incidence of prolonged time to extubation (0 vs. 6%, P = 0.029) in the DES group were significantly decreased compared with those in the TIVA group. The patients in the DES group displayed shorter OR exit time as compared with that in the TIVA group (7.3 [5.5-8.7] vs. 10.8 [9.3-12.3] min, P < 0.001). The patients using desflurane exhibited more stable hemodynamics during surgery than the patients using propofol-based TIVA, as demonstrated by lower incidences of hypotension (1% vs. 22%, P < 0.001), bradycardia (2% vs. 13%, P = 0.002), and OCR (17% vs. 44%, P < 0.001).

CONCLUSION

DES enhanced the ophthalmic OR efficiency by reducing the extubation time and OR exit time, and provided more stable hemodynamics intra-operatively than TIVA in patients undergoing strabismus ambulatory surgery.

TRIAL REGISTRATION

ClinicalTrials.gov, No. NCT02922660; https://clinicaltrials.gov/ct2/show/NCT02922660?id=NCT02922660&draw=2&rank=1.

The influence of cardiac output on propofol and fentanyl pharmacokinetics and pharmacodynamics in patients undergoing abdominal aortic surgery

Cardiac output (CO) is expected to affect elimination and distribution of highly extracted and perfusion rate-limited drugs. This work was undertaken to quantify the effect of CO measured by the pulse pressure method on pharmacokinetics and pharmacodynamics of propofol and fentanyl administrated during total intravenous anesthesia (TIVA). The data were obtained from 22 ASA III patients undergoing abdominal aortic surgery. Propofol was administered via target-controlled infusion system (Diprifusor) and fentanyl was administered at a dose of 2-3 µg/kg each time analgesia appeared to be inadequate. Hemodynamic measurements as well as bispectral index were monitored and recorded throughout the surgery. Data analysis was performed by using a non-linear mixed-effect population modeling (NONMEM 7.4 software). Three compartment models that incorporated blood flows as parameters were used to describe propofol and fentanyl pharmacokinetics. The delay of the anesthetic effect, with respect to plasma concentrations, was described using a biophase (effect) compartment. The bispectral index was linked to the propofol and fentanyl effect site concentrations through a synergistic Emax model. An empirical linear model was used to describe CO changes observed during the surgery. Cardiac output was identified as an important predictor of propofol and fentanyl pharmacokinetics. Consequently, it affected the depth of anesthesia and the recovery time after propofol-fentanyl TIVA infusion cessation. The model predicted (not observed) CO values correlated best with measured responses. Patients' age was identified as a covariate affecting the rate of CO changes during the anesthesia leading to age-related difference in individual patient's responses to both drugs.

Neurophysiologic Intraoperative Monitoring for Spine Surgery: A Practical Guide From Past to Present

Intraoperative neuromonitoring was introduced in the second half of the 20th century with the goal of preventing patient morbidity for patients undergoing complex operations of the central and peripheral nervous system. Since its early use for scoliosis surgery, the growth and utilization of IOM techniques expanded dramatically over the past 50 years to include spinal tumor resection and evaluation of cerebral ischemia. The importance of IOM has been broadly acknowledged, and in 1989, the American Academy of Neurology (AAN) released a statement that the use of SSEPs should be standard-of-care during spine surgery. In 2012, both the AAN and the American Clinical Neurophysiology Society (ACNS) recommended that: "Intraoperative monitoring (IOM) using SSEPs and transcranial MEPs be established as an effective means of predicting an increased risk of adverse outcomes, such as paraparesis, paraplegia, and quadriplegia, in spinal surgery." With a multimodal approach that combines SSEPs, MEPs, and sEMG with tEMG and D waves, as appropriate, sensitivity and specificity can be maximized for the diagnosis of reversible insults to the spinal cord, nerve roots, and peripheral nerves. As with most patient safety efforts in the operating room, IOM requires contributions from and communication between a number of different teams. This comprehensive review of neuromonitoring techniques for surgery on the central and peripheral nervous system will highlight the technical, surgical and anesthesia factors required to optimize outcomes. In addition, this review will discuss important trouble shooting measures to be considered when managing ION changes concerning for potential injury.

Bolus pharmacokinetics: moving beyond mass-based dosing to guide drug administration

Despite the common approach of bolus drug dosing using a patient's mass, a more tailored approach would be to use empirically derived pharmacokinetic models. Previously, this could only be possible though the use of computer simulation using programs which are rarely available in clinical practice. Through mathematical manipulations and approximations, a simplified set of equations is demonstrated that can identify a bolus dose required to achieve a specified target effect site concentration. The proposed solution is compared against simulations of a wide variety of pharmacokinetic models. This set of equations provides a near-identical solution to the simulation approach. A boundary condition is established to ensure the derived equations have an acceptable error. This approach may allow for more precise administration of medications with the use of point of care technology and potentially allows for pharmacokinetic dosing in artificial intelligence problems.

Critical view of the effect site modelling of propofol

Target controlled infusion (TCI) of Propofol has been the subject of discussion during its 20 years of use, including the validity of the models that represent the course of the effect, such as: Are the different EEG indexes representative of the effect? Is the reactivity of the EEG index used to build models comparable to each other? What is the real reacting time of each monitor? Is the ke0 influenced by the infusion speed? Is the ke0 or the time to peak effect affected by age? How valid are the current Emax models? Are the induction and wakening simple mirror phenomenon as they are represented in the E max models? This review discusses issues related to the complexity and difficulty in obtaining a representation of the effect, and the lack of agreed definitions to be able to construct representative models of the temporary installation of the effect of Propofol for its use in TCI.

Postoperative Cardiac Troponin I Thresholds Associated With 1-Year Cardiac Mortality After Adult Cardiac Surgery: An Attempt to Link Risk Stratification With Management Stratification in an Observational Study

OBJECTIVES

Cardiac troponin (cTn) concentrations are measured routinely in some centers after cardiac surgery as part of risk stratification, but there are no data on how increased cTn concentrations could change patients' management. The aim of this study was to estimate relevant cTnI thresholds and identify potential interventions (additional monitoring/therapeutic interventions) that could be part of management changes of patients with cTnI greater than relevant thresholds.

DESIGN

Retrospective, single-center, observational study.

SETTING

Bichat-Claude Bernard Hospital, Paris, France, between January 1, 2009, and December 31, 2012.

PARTICIPANTS

Consecutive adult patients undergoing cardiac surgery.

MEASUREMENTS AND MAIN RESULTS

cTnI was measured on the 20th postoperative hour. Causes of death and possible interventions were determined by analysis of individual medical records. cTnI thresholds for 1-year cardiac mortality with a specificity >80% were calculated. For this study, 3,228 procedures were analyzed; 129 deaths occurred (4%), 83 of which (2.6%) were cardiac deaths. Threshold cTnI values were 4.2 µg/L for coronary artery bypass grafting (95% confidence interval [CI] 3.9-4.5) and 10.7 µg/L for non-coronary artery bypass grafting (95% CI 10.0-11.3). In multivariable analysis, the EuroSCORE II (odds ratio 1.1 [95% CI 1.06-1.13]; p < 0.001) and cTnI concentrations greater than the thresholds (odds ratio 5.62 [95% CI 3.37-9.37]; p < 0.001) were associated with significantly increased risk of death. The additive and absolute Net Reclassification Index were 0.288% and 14.1%, respectively, for a logistic model including cTnI and EuroSCORE II (area under the curve C-index 0.82 [95% CI 0.77-0.87]) compared with a model including only EuroSCORE II (area under the curve C-index 0.80 [95% CI 0.75-0.84]). Fifty-three of the 83 patients who experienced cardiac death (64%) had a cTnI concentration greater than the threshold, and an intervention was deemed possible in 47 of those 53 (89%) (mostly patients with mild postoperative cardiac dysfunction). For noncardiac deaths, 28% of patients had a cTnI concentration greater than the threshold and no interventions were deemed possible.

CONCLUSIONS

In an attempt to evolve from risk to management stratification, this study's results identified a subgroup of patients with mild cardiac dysfunction and a cTnI concentration greater than the threshold who could be the target for interventions in future validation studies concerning changes in patient management.

Developing a personalized closed-loop controller of medically-induced coma in a rodent model

OBJECTIVE

Personalized automatic control of medically-induced coma, a critical multi-day therapy in the intensive care unit, could greatly benefit clinical care and further provide a novel scientific tool for investigating how the brain response to anesthetic infusion rate changes during therapy. Personalized control would require real-time tracking of inter- and intra-subject variabilities in the brain response to anesthetic infusion rate while simultaneously delivering the therapy, which has not been achieved. Current control systems for medically-induced coma require a separate offline model fitting experiment to deal with inter-subject variabilities, which would lead to therapy interruption. Removing the need for these offline interruptions could help facilitate clinical feasbility. In addition, current systems do not track intra-subject variabilities. Tracking intra-subject variabilities is essential for studying whether or how the brain response to anesthetic infusion rate changes during therapy. Further, such tracking could enhance control precison and thus help facilitate clinical feasibility.

APPROACH

Here we develop a personalized closed-loop anesthetic delivery (CLAD) system in a rodent model that tracks both inter- and intra-subject variabilities in real time while simultaneously controlling the anesthetic in closed loop. We tested the CLAD in rats by administrating propofol to control the electroencephalogram (EEG) burst suppression. We first examined whether the CLAD can remove the need for offline model fitting interruption. We then used the CLAD as a tool to study whether and how the brain response to anesthetic infusion rate changes as a function of changes in the depth of medically-induced coma. Finally, we studied whether the CLAD can enhance control compared with prior systems by tracking intra-subject variabilities.

MAIN RESULTS

The CLAD precisely controlled the EEG burst suppression in each rat without performing offline model fitting experiments. Further, using the CLAD, we discovered that the brain response to anesthetic infusion rate varied during control, and that these variations correlated with the depth of medically-induced coma in a consistent manner across individual rats. Finally, tracking these variations reduced control bias and error by more than 70% compared with prior systems.

SIGNIFICANCE

This personalized CLAD provides a new tool to study the dynamics of brain response to anesthetic infusion rate and has significant implications for enabling clinically-feasible automatic control of medically-induced coma.

PKQuest: PBPK modeling of highly lipid soluble and extracellular solutes

One of the primary objectives of physiologically based pharmacokinetics (PBPK) is the prediction of a drug’s pharmacokinetics just from knowledge of its physicochemical structure. Unfortunately, at present, the accuracy of this prediction is limited for most drugs because of uncertainty about the drug’s organ/blood partition coefficient (K). However, there are two classes of solutes which are exceptions to this: 1) the highly lipid soluble (HLS) solutes, and 2) the extracellular (ECS) solutes. Since the HLS drugs (eg, volatile anesthetics, propofol, cannabinol) have lipid/water partition coefficients (P_L/W) of 100 or greater, their K is dominated by the tissue fat fraction and one can accurately predict K just from in vitro measurements of P_L/W along with prior anatomic measurements of the fat fraction of the organs in the PBPK model. Since the ECS drugs, such as most antibiotics, cannot penetrate cells, they are not subject to the intracellular binding that complicates the prediction of K for the weak bases and acids. The ECS K is determined primarily by plasma and interstitial albumin binding and can be predicted from in vitro measurements of plasma albumin binding along with prior measurements of interstitial tissue volume and albumin concentrations. This review provides an in depth discussion of the PBPK modeling of these two drug classes along with many specific clinical examples illustrating the good PBPK predictions possible with just zero (volatile anesthetics) or 1 (the clearance) adjustable parameter. The PBPK analysis uses PKQuest, a freely distributed, general purpose pharmacokinetic program. PKQuest is designed so that application to the HLS and ECS solute classes is especially easy. The user only needs to enter the specific parameters that are required to characterize the drug (eg, P_L/W for HLS or plasma albumin binding for ECS) with all the other PBPK parameters (organ blood flow, fat fraction, extracellular volumes, etc.) are set by default.

Target and Tissue Selectivity Prediction by Integrated Mechanistic Pharmacokinetic-Target Binding and Quantitative Structure Activity Modeling

Selectivity is an important attribute of effective and safe drugs, and prediction of in vivo target and tissue selectivity would likely improve drug development success rates. However, a lack of understanding of the underlying (pharmacological) mechanisms and availability of directly applicable predictive methods complicates the prediction of selectivity. We explore the value of combining physiologically based pharmacokinetic (PBPK) modeling with quantitative structure-activity relationship (QSAR) modeling to predict the influence of the target dissociation constant (K _D) and the target dissociation rate constant on target and tissue selectivity. The K _D values of CB1 ligands in the ChEMBL database are predicted by QSAR random forest (RF) modeling for the CB1 receptor and known off-targets (TRPV1, mGlu5, 5-HT1a). Of these CB1 ligands, rimonabant, CP-55940, and Δ⁸-tetrahydrocanabinol, one of the active ingredients of cannabis, were selected for simulations of target occupancy for CB1, TRPV1, mGlu5, and 5-HT1a in three brain regions, to illustrate the principles of the combined PBPK-QSAR modeling. Our combined PBPK and target binding modeling demonstrated that the optimal values of the K _D and k _off for target and tissue selectivity were dependent on target concentration and tissue distribution kinetics. Interestingly, if the target concentration is high and the perfusion of the target site is low, the optimal K _D value is often not the lowest K _D value, suggesting that optimization towards high drug-target affinity can decrease the benefit-risk ratio. The presented integrative structure-pharmacokinetic-pharmacodynamic modeling provides an improved understanding of tissue and target selectivity.

Propofol Breath Monitoring as a Potential Tool to Improve the Prediction of Intraoperative Plasma Concentrations

INTRODUCTION

Monitoring of drug concentrations in breathing gas is routinely being used to individualize drug dosing for the inhalation anesthetics. For intravenous anesthetics however, no decisive evidence in favor of breath concentration monitoring has been presented up until now. At the same time, questions remain with respect to the performance of currently used plasma pharmacokinetic models implemented in target-controlled infusion systems. In this study, we investigate whether breath monitoring of propofol could improve the predictive performance of currently applied, target-controlled infusion models.

METHODS

Based on data from a healthy volunteer study, we developed an addition to the current state-of-the-art pharmacokinetic model for propofol, to accommodate breath concentration measurements. The potential of using this pharmacokinetic (PK) model in a Bayesian forecasting setting was studied using a simulation study. Finally, by introducing bispectral index monitor (BIS) measurements and the accompanying BIS models into our PK model, we investigated the relationship between BIS and predicted breath concentrations.

RESULTS AND DISCUSSION

We show that the current state-of-the-art pharmacokinetic model is easily extended to reliably describe propofol kinetics in exhaled breath. Furthermore, we show that the predictive performance of the a priori model is improved by Bayesian adaptation based on the measured breath concentrations, thereby allowing further treatment individualization and a more stringent control on the targeted plasma concentrations during general anesthesia. Finally, we demonstrated concordance between currently advocated BIS models, relying on predicted effect-site concentrations, and our new approach in which BIS measurements are derived from predicted breath concentrations.

Planning for operating room efficiency and faster anesthesia wake-up time in open major upper abdominal surgery

Reducing anesthesia-controlled time (ACT) may improve operation room (OR) efficiency result from different anesthetic techniques. However, the information about the difference in ACT between desflurane (DES) anesthesia and propofol-based total intravenous anesthesia (TIVA) techniques for open major upper abdominal surgery under general anesthesia (GA) is not available in the literature.This retrospective study uses our hospital database to analyze the ACT of open major upper abdominal surgery without liver resection after either desflurane/fentanyl-based anesthesia or TIVA via target-controlled infusion with fentanyl/propofol from January 2010 to December 2011. The various time intervals including waiting for anesthesia time, anesthesia time, surgical time, extubation time, exit from OR after extubation, total OR time, and postanesthetic care unit (PACU) stay time and percentage of prolonged extubation (≥15 minutes) were compared between these 2 anesthetic techniques.We included data from 343 patients, with 159 patients receiving TIVA and 184 patients receiving DES. The only significant difference is extubation time, TIVA was faster than the DES group (8.5 ± 3.8 vs 9.4 ± 3.7 minutes; P = 0.04). The factors contributed to prolonged extubation were age, gender, body mass index, DES anesthesia, and anesthesia time.In our hospital, propofol-based TIVA by target-controlled infusion provides faster emergence compared with DES anesthesia; however, it did not improve OR efficiency in open major abdominal surgery. Older, male gender, higher body mass index, DES anesthesia, and lengthy anesthesia time were factors that contribute to extubation time.

Comparison of Anesthesia-Controlled Operating Room Time between Propofol-Based Total Intravenous Anesthesia and Desflurane Anesthesia in Open Colorectal Surgery: A Retrospective Study

We conducted a retrospective study to investigate the anesthesia-controlled time and factors that contribute to prolonged extubation in open colorectal surgery. Using our hospital database, demographic data, various time intervals (waiting for anesthesia time, anesthesia time, surgical time, emergence time, exit from operating room after extubation, total operating room time, and post-anesthesia care unit stay time), and incidence of prolonged extubation (≥ 15 mins), were compared between patients who received desflurane/fentanyl-based anesthesia and total intravenous anesthesia via target-controlled infusion with fentanyl/propofol. Logistic regression analyses were performed to assess the association between variables that contributed to prolonged extubation. In conclusion, the anesthesia-controlled time was similar in desflurane anesthesia and propofol-based total intravenous anesthesia for open colorectal surgery in our hospital. Surgical time greater than 210 minutes, as well as age, contributed to prolonged extubation.

Physiologically-based pharmacokinetic models: approaches for enabling personalized medicine

Personalized medicine strives to deliver the 'right drug at the right dose' by considering inter-person variability, one of the causes for therapeutic failure in specialized populations of patients. Physiologically-based pharmacokinetic (PBPK) modeling is a key tool in the advancement of personalized medicine to evaluate complex clinical scenarios, making use of physiological information as well as physicochemical data to simulate various physiological states to predict the distribution of pharmacokinetic responses. The increased dependency on PBPK models to address regulatory questions is aligned with the ability of PBPK models to minimize ethical and technical difficulties associated with pharmacokinetic and toxicology experiments for special patient populations. Subpopulation modeling can be achieved through an iterative and integrative approach using an adopt, adapt, develop, assess, amend, and deliver methodology. PBPK modeling has two valuable applications in personalized medicine: (1) determining the importance of certain subpopulations within a distribution of pharmacokinetic responses for a given drug formulation and (2) establishing the formulation design space needed to attain a targeted drug plasma concentration profile. This review article focuses on model development for physiological differences associated with sex (male vs. female), age (pediatric vs. young adults vs. elderly), disease state (healthy vs. unhealthy), and temporal variation (influence of biological rhythms), connecting them to drug product formulation development within the quality by design framework. Although PBPK modeling has come a long way, there is still a lengthy road before it can be fully accepted by pharmacologists, clinicians, and the broader industry.

Sleep Homeostasis and General Anesthesia: Are Fruit Flies Well Rested after Emergence from Propofol?

BACKGROUND

Shared neurophysiologic features between sleep and anesthetic-induced hypnosis indicate a potential overlap in neuronal circuitry underlying both states. Previous studies in rodents indicate that preexisting sleep debt discharges under propofol anesthesia. The authors explored the hypothesis that propofol anesthesia also dispels sleep pressure in the fruit fly. To the authors' knowledge, this constitutes the first time propofol has been tested in the genetically tractable model, Drosophila melanogaster.

METHODS

Daily sleep was measured in Drosophila by using a standard locomotor activity assay. Propofol was administered by transferring flies onto food containing various doses of propofol or equivalent concentrations of vehicle. High-performance liquid chromatography was used to measure the tissue concentrations of ingested propofol. To determine whether propofol anesthesia substitutes for natural sleep, the flies were subjected to 10-h sleep deprivation (SD), followed by 6-h propofol exposure, and monitored for subsequent sleep.

RESULTS

Oral propofol treatment causes anesthesia in flies as indicated by a dose-dependent reduction in locomotor activity (n = 11 to 41 flies from each group) and increased arousal threshold (n = 79 to 137). Recovery sleep in flies fed propofol after SD was delayed until after flies had emerged from anesthesia (n = 30 to 48). SD was also associated with a significant increase in mortality in propofol-fed flies (n = 44 to 46).

CONCLUSIONS

Together, these data indicate that fruit flies are effectively anesthetized by ingestion of propofol and suggest that homologous molecular and neuronal targets of propofol are conserved in Drosophila. However, behavioral measurements indicate that propofol anesthesia does not satisfy the homeostatic need for sleep and may compromise the restorative properties of sleep.

A Pharmacokinetics-Neural Mass Model (PK-NMM) for the Simulation of EEG Activity during Propofol Anesthesia

Modeling the effects of anesthetic drugs on brain activity is very helpful in understanding anesthesia mechanisms. The aim of this study was to set up a combined model to relate actual drug levels to EEG dynamics and behavioral states during propofol-induced anesthesia. We proposed a new combined theoretical model based on a pharmacokinetics (PK) model and a neural mass model (NMM), which we termed PK-NMM--with the aim of simulating electroencephalogram (EEG) activity during propofol-induced general anesthesia. The PK model was used to derive propofol effect-site drug concentrations (C(eff)) based on the actual drug infusion regimen. The NMM model took C(eff) as the control parameter to produce simulated EEG-like (sEEG) data. For comparison, we used real prefrontal EEG (rEEG) data of nine volunteers undergoing propofol anesthesia from a previous experiment. To see how well the sEEG could describe the dynamic changes of neural activity during anesthesia, the rEEG data and the sEEG data were compared with respect to: power-frequency plots; nonlinear exponent (permutation entropy (PE)); and bispectral SynchFastSlow (SFS) parameters. We found that the PK-NMM model was able to reproduce anesthesia EEG-like signals based on the estimated drug concentration and patients' condition. The frequency spectrum indicated that the frequency power peak of the sEEG moved towards the low frequency band as anesthesia deepened. Different anesthetic states could be differentiated by the PE index. The correlation coefficient of PE was 0.80 ± 0.13 (mean ± standard deviation) between rEEG and sEEG for all subjects. Additionally, SFS could track the depth of anesthesia and the SFS of rEEG and sEEG were highly correlated with a correlation coefficient of 0.77 ± 0.13. The PK-NMM model could simulate EEG activity and might be a useful tool for understanding the action of propofol on brain activity.

Evidence of Different Propofol Pharmacokinetics under Short and Prolonged Infusion Times in Rabbits

Propofol is an anaesthetic widely used in both human beings and animals. However, the characterization of propofol pharmacokinetics (PK) is not well understood when long-term infusions are used. The main objective of this study was to explore the PK behaviour of propofol in a rabbit model during short and prolonged propofol infusions and to develop an internally validated PK model, for propofol dose individualization in the rabbit for future use. Population 1 (P1) was constituted by seven New Zealand rabbits and was used to characterize the PK profile of propofol at short infusions. Animals were anaesthetized with a bolus of 20 mg/kg, followed by an infusion rate of 50 mg/kg/hr of propofol at 1%, which was then maintained for 30 min. A second rabbit population (P2, n = 7) was sedated according to reflexes responses and Index of Consciousness values, for 20 consecutive hours using propofol 2% aiming at characterizing propofol behaviour at long-term infusions. Clinical data and blood samples were collected at specific time-points in both populations. Propofol plasma concentrations were determined by gas chromatography/ion trap mass spectrometry. The NONMEM VII software was used to evaluate the relationships between dose and plasma concentrations. A linear two-compartment model with different central compartment volume and plasma clearance (separately modelled in the two populations) was the one that best described propofol concentrations. The time course of propofol plasma concentrations was well characterized by the PK model developed, which simultaneously accounts for propofol short- and long-term infusions and can be used to optimize future PK studies in rabbits.

Pharmacokinetic-pharmacodynamic modelling in anaesthesia

Anaesthesiologists adjust drug dosing, administration system and kind of drug to the characteristics of the patient. They then observe the expected response and adjust dosing to the specific requirements according to the difference between observed response, expected response and the context of the surgery and the patient. The approach above can be achieved because on one hand quantification technology has made significant advances allowing the anaesthesiologist to measure almost any effect by using noninvasive, continuous measuring systems. On the other the knowledge on the relations between dosing, concentration, biophase dynamics and effect as well as detection of variability sources has been achieved as being the benchmark specialty for pharmacokinetic-pharmacodynamic (PKPD) modelling. The aim of the review is to revisit the most common PKPD models applied in the field of anaesthesia (i.e. effect compartmental, turnover, drug-receptor binding and drug interaction models) through representative examples. The effect compartmental model has been widely used in this field and there are multiple applications and examples. The use of turnover models has been limited mainly to describe respiratory effects. Similarly, cases in which the dissociation process of the drug-receptor complex is slow compared with other processes relevant to the time course of the anaesthetic effect are not frequent in anaesthesia, where in addition to a rapid onset, a fast offset of the response is required. With respect to the characterization of PD drug interactions different response surface models are discussed. Relevant applications that have changed the way modern anaesthesia is practiced are also provided.

Prediction of nomegestrol acetate pharmacokinetics in healthy female adolescents and adults by whole-body physiology-based pharmacokinetic modelling and clinical validation

OBJECTIVE

Nomegestrol acetate (NOMAC), a selective progestogen, and 17β-estradiol (E2), which is identical to endogenous oestrogen, are components of a new monophasic combined oral contraceptive--NOMAC/E2. This study aimed to compare pharmacokinetics (PK) of NOMAC in adolescent and adult women following a single dose of NOMAC/E2.

STUDY DESIGN

Healthy postmenarcheal adolescent (14-17years) and adult (18-50years) women received a single dose of NOMAC/E2 (2.5mg/1.5mg) in this single-centre, open-label, parallel-group Phase 1 study (EudraCT# 2008-002142-38). Blood samples were obtained for PK analysis, and concentrations of NOMAC, E2 and its metabolite estrone (E1) were determined for up to 129h following dosing to obtain PK data. An independent whole-body physiology-based pharmacokinetic (WB-PBPK) simulation model of NOMAC based on an independent Phase 3 dataset was used to scale NOMAC concentration-time plots to adolescents.

RESULTS

Overall, 52 women were screened, of whom 30 (15 adolescents and 15 adults) were enrolled. No statistically significant differences were observed between the adolescent and adult groups for the clinically evaluated NOMAC PK parameters [maximum concentration (Cmax), area under the curve (AUC) and half-life (t1/2)]. The PK of E2 and E1 showed extensive overlap between both age groups. The WB-PBPK model accurately predicted NOMAC AUC and Cmax values in both groups.

CONCLUSIONS

No differences were observed in the clinically evaluated PK parameters for NOMAC between adolescent and adult women after a single dose of NOMAC/E2. The WB-PBPK model accurately predicted NOMAC PK data (EudraCT# 2008-002142-38).

IMPLICATIONS

PK studies in adolescents are challenging because of ethical considerations. The whole-body physiology-based model described here complements classic noncompartmental and population PK approaches. The utility of this method is its ability to expand to adolescent postmenarcheal girls by using virtual postmenarcheal adolescent population data and applying physiological scaling.

Organ-specific microcirculatory mass transport of oxycodone in humans: clinical implications for therapeutic use

OBJECTIVES

To begin to address the problem of heterogeneity of distribution of oxycodone (OC) in humans, we developed an organ-specific microcirculatory capillary-tissue exchange 2-compartment model for studying regional OC mass transport.

MATERIALS AND METHODS

The model was developed in silico. It quantifies OC's organ-specific mass transport rates, clearances and recycling, and it considers the effects of blood flow on OC's convective and diffusive transport.

RESULTS

What is new is the finding that OC undergoes local recycling at the level of organ-specific capillary-tissue exchange units in humans. Results indicate recycled OC occurs in sufficient amounts to function as a reusable source of circulating OC; which has important implications for OC dosing. Results show the brain, which is central to OC effects only receives about 8% of OC delivered to all organs via the microcirculation. This suggests that differential regulation of receptor binding, trafficking, internalization, or desensitization in the brain likely plays a dominant role in OC's central analgesic effects.

DISCUSSION

Organ-specific OC mass transport kinetics provide new information for OC dosing in pain management. The model promotes patient safety in opioid prescribing because it allows predictions to be made about the relative contribution that OC recycling makes to circulating OC levels. The model indicates that pharmacologic modulation of the microcirculation may give way to site-specific delivery of opioids in the future. Our study demonstrates that translation of bench in silico research data into clinical practice, although still challenging, is feasible and can assist in OC dose regimen design for patient safety.

Determination of the most influential sources of variability in tacrolimus trough blood concentrations in adult liver transplant recipients: a bottom-up approach

Tacrolimus, an immunosuppressant drug, presents a narrow therapeutic window and a large pharmacokinetic variability with poor correlation between drug dosing regimen and blood concentration. The objective was to identify predictive factors influencing tacrolimus trough concentrations (C0) using a bottom-up approach. A physiologically based pharmacokinetic (PBPK) model of tacrolimus was proposed, taking into account the body weight, the proportion of fat (P(fat)), hematocrit, lipid fraction of organs, typical intrinsic clearance (CLi(typ)), CYP3A5 genotype of liver donor, plasma unbound fraction of tacrolimus (fu(p)), and concomitant drugs (CYP3A4 inhibitors). For the evaluation of the PBPK model, mean C0 and concentrations 2 h after oral dose of tacrolimus were compared with those from 66 liver transplant recipients included in a multicentric pharmacokinetic study and were found very close. Tacrolimus concentration profiles were simulated in a virtual population defined by a set of covariate values similar to those from the real population. The sensitivity of tacrolimus C0 with respect to each covariate has been tested to identify the most influential ones. With the range of covariate values tested, the impact of each covariate on tacrolimus C0 may be ranked as follows: fu(p), CLi(typ), bioavailability, body weight, hematocrit, CYP3A5 polymorphism, P(fat), and CYP3A4 inhibitory drug-drug interactions. Values for initial dosing regimen of tacrolimus in order to reach a C0 of 10 ng/ml at day 5 (assuming a constant dosing schedule) as a function of CYP3A5 donor genotype and patient's hematocrit and body weight are proposed.

Evaluation of propofol anesthesia in morbidly obese children and adolescents

BACKGROUND

Poor characterization of propofol pharmacokinetics and pharmacodynamics in the morbidly obese (MO) pediatric population poses dosing challenges. This study was conducted to evaluate propofol total intravenous anesthesia (TIVA) in this population.

METHODS

After IRB approval, a prospective study was conducted in 20 MO children and adolescents undergoing laparoscopic surgery under clinically titrated propofol TIVA. Propofol doses/infusion rates, hemodynamic variables, times to induction and emergence, and postoperative occurrence of respiratory adverse events (RAE) were recorded, along with intraoperative blinded Bispectral Index/BIS and postoperative Ramsay sedation scores (RSS). Study subjects completed awareness questionnaires on postoperative days 1 and 3. Propofol concentrations were obtained at predetermined intra- and post-operative time points.

RESULTS

Study subjects ranged 9 - 18 years (age) and 97 - 99.9% (BMI for age percentiles). Average percentage variability of hemodynamic parameters from baseline was ≈ 20%. Patients had consistently below target BIS values (BIS < 40 for >90% of maintenance phase), delayed emergence (25.8 ± 22 minutes), increased somnolence (RSS ≥ 4) in the first 30 minutes of recovery from anesthesia and 30% incidence of postoperative RAE, the odds for which increased by 14% per unit increase in BMI (p ≤ 0.05). Mean propofol concentration was 6.2 mg/L during maintenance and 1.8 mg/L during emergence from anesthesia.

CONCLUSIONS

Our findings indicate clinical overestimation of propofol requirements and highlight the challenges of clinically titrated propofol TIVA in MO adolescents. In this setting, it may be advantageous to titrate propofol to targeted BIS levels until more accurate weight-appropriate dosing regimens are developed, to minimize relative overdosing and its consequences.

Application of a physiologically based pharmacokinetic model to assess propofol hepatic and renal glucuronidation in isolation: utility of in vitro and in vivo data

A physiologically based pharmacokinetic (PBPK) modeling approach was used to assess the prediction accuracy of propofol hepatic and extrahepatic metabolic clearance and to address previously reported underprediction of in vivo clearance based on static in vitro-in vivo extrapolation methods. The predictive capacity of propofol intrinsic clearance data (CLint) obtained in human hepatocytes and liver and kidney microsomes was assessed using the PBPK model developed in MATLAB software. Microsomal data obtained by both substrate depletion and metabolite formation methods and in the presence of 2% bovine serum albumin were considered in the analysis. Incorporation of hepatic and renal in vitro metabolic clearance in the PBPK model resulted in underprediction of propofol clearance regardless of the source of in vitro data; the predicted value did not exceed 35% of the observed clearance. Subsequently, propofol clinical data from three dose levels in intact patients and anhepatic subjects were used for the optimization of hepatic and renal CLint in a simultaneous fitting routine. Optimization process highlighted that renal glucuronidation clearance was underpredicted to a greater extent than liver clearance, requiring empirical scaling factors of 17 and 9, respectively. The use of optimized clearance parameters predicted hepatic and renal extraction ratios within 20% of the observed values, reported in an additional independent clinical study. This study highlights the complexity involved in assessing the contribution of extrahepatic clearance mechanisms and illustrates the application of PBPK modeling, in conjunction with clinical data, to assess prediction of clearance from in vitro data for each tissue individually.

A two-compartment effect site model describes the bispectral index after different rates of propofol infusion

Different estimates of the rate constant for the effect site distribution (k(e0)) of propofol, depending on the rate and duration of administration, have been reported. This analysis aimed at finding a more general pharmacodynamic model that could be used when the rate of administration is changed during the treatment. In a cross-over study, 21 healthy volunteers were randomised to receive a 1 min infusion of 2 mg/kg of propofol at one occasion, and a 1 min infusion of 2 mg/kg of propofol immediately followed by a 29 min infusion of 12 mg kg(-1) h(-1) of propofol at another occasion. Arterial plasma concentrations of propofol were collected up to 4 h after dosing, and BIS was collected before start of infusion and until the subjects were fully awake. The population pharmacokinetic-pharmacodynamic analysis was performed using NONMEM VI. A four-compartment PK model with time-dependent elimination and distribution described the arterial propofol concentrations, and was used as input to the pharmacodynamic model. A standard effect compartment model could not accurately describe the delay in the effects of propofol for both regimens, whereas a two-compartment effect site model significantly improved the predictions. The two-compartment effect site model included a central and a peripheral effect site compartment, possibly representing a distribution within the brain, where the decrease in BIS was linked to the central effect site compartment concentrations through a sigmoidal E(max) model.

PKQuest_Java: free, interactive physiologically based pharmacokinetic software package and tutorial

Background

Physiologically based pharmacokinetics (PBPK) uses a realistic organ model to describe drug kinetics. The blood-tissue exchange of each organ is characterized by its volume, perfusion, metabolism, capillary permeability and blood/tissue partition coefficient. PBPK applications require both sophisticated mathematical modeling software and a reliable complete set of physiological parameters. Currently there are no software packages available that combine ease of use with the versatility that is required of a general PBPK program.

Findings

The program is written in Java and is available for free download at . Included in the download is a detailed tutorial that discusses the pharmacokinetics of 6 solutes (D₂O, amoxicillin, desflurane, propofol, ethanol and thiopental) illustrated using experimental human pharmacokinetic data. The complete PBPK description for each solute is stored in Excel spreadsheets that are included in the download. The main features of the program are: 1) Intuitive and versatile interactive interface; 2) Absolute and semi-logarithmic graphical output; 3) Pre-programmed optimized human parameter data set (but, arbitrary values can be input); 4) Time dependent changes in the PBPK parameters; 5) Non-linear parameter optimization; 6) Unique approach to determine the oral "first pass metabolism" of non-linear solutes (e.g. ethanol); 7) Pulmonary perfusion/ventilation heterogeneity for volatile solutes; 8) Input and output of Excel spreadsheet data; 9) Antecubital vein sampling.

Conclusion

PKQuest_Java is a free, easy to use, interactive PBPK software routine. The user can either directly use the pre-programmed optimized human or rat data set, or enter an arbitrary data set. It is designed so that drugs that are classified as "extracellular" or "highly fat soluble" do not require information about tissue/blood partition coefficients and can be modeled by a minimum of user input parameters. PKQuest_Java, along with the included tutorial, could be used as the basis of an interactive, on-line, pharmacokinetic course.

Whole body physiologically-based pharmacokinetic models: their use in clinical drug development

BACKGROUND

Whole-body physiologically-based pharmacokinetic (WB-PBPK) models mathematically describe an organism as a closed circulatory system consisting of compartments that represent the organs important for compound absorption, distribution, metabolism and elimination.

OBJECTIVES

To review the current state of WB-PBPK model use in the clinical phases of drug development.

METHODS

A qualitative description of the WB-PBPK model structure is included along with a review of the varying methods available for input parameterisation. Current and potential WB-PBPK model application in clinical development is discussed.

CONCLUSIONS

This modelling tool is at present used for small and large molecule drug development primarily as a means to scale pharmacokinetics from animals to humans based on physiology. The pharmaceutical industry is active in employing these models to clinical drug development although the applications in use now are narrow in comparison to the potential. Expanded integration of WB-PBPK models into the drug development process will only be achieved with staff training, managerial will, success stories and regulatory agency openness.

Fospropofol disodium injection for the sedation of patients undergoing colonoscopy

Sedation plays a central role in making colonoscopy tolerable for patients and feasible for the endoscopist to perform. The array of agents used for endoscopic sedation continues to evolve. Fospropofol (FP), a prodrug of propofol with a slower pharmacokinetic profile, is currently under evaluation for use during endoscopic procedures. Preliminary data suggests that FP dosed at 6.5 mg/kg is well tolerated by most patients with perineal paresthesias being the most commonly experienced adverse effect. This article will examine the current literature on the use of FP for the sedation of patients undergoing colonoscopy, highlighting the pharmacokinetics, pharmacodynamics, risks, and common adverse events associated with the novel sedative/hypnotic.

In vitro release of propofol and binding capacity with regard to plasma constituents

PURPOSE

New evidence suggests that the anesthetic effect of parenteral propofol emulsions varies between commercial preparations. We examined and compared different propofol preparations to determine propofol release and binding capacity with regard to plasma lipoproteins and albumin.

METHODS

We created a novel assay consisting of microtiter plates coated with either low-density lipoprotein (LDL) or albumin to analyze propofol binding kinetics. Using high performance liquid chromatography (HPLC), we measured propofol release from the oily phase and the corresponding amount of propofol bound to the plates in a time-dependent manner and at equilibrium conditions attained after 30 min of incubation at 37 degrees C. The concentrations of free propofol in the aqueous phase of different propofol preparations - Diprivan, and the generic formulations Propofol "Fresenius" (1% and 2% propofol) and Propofol-Lipuro - were analyzed using ultracentrifugation or dialysis for phase separation. Finally, we investigated the effect of isolated lipoprotein fractions on propofol release.

RESULTS

Propofol bound to LDL-coated plates with approximately twofold higher affinity than to albumin-coated plates. No significant differences in total propofol release were observed between preparations. Moreover, similar amounts of free propofol were observed in the aqueous phase of all products tested (1% propofol preparations: 18 microg/ml; 2% propofol preparations: 35 microg/ml), except for the medium-chain and long-chain triglyceride (MCT/LCT) preparation studied, in which the concentration of free propofol was lower. Lipoproteins had no effect on propofol release, except for high-density lipoprotein (HDL), which triggered almost 100% release from the oily phase at HDL concentrations above 1000 microg/ml.

CONCLUSIONS

No differences were observed between the binding/release capacity and lipoprotein interactions of any of the propofol preparations tested. We propose that clinical observations of inconsistent propofol activity are related to variations in the lipoprotein profile, enzyme activity or genetic disorders of individual patients, rather than to the propofol preparation itself.

Simulation of propofol anaesthesia for intracranial decompression using brain hypothermia treatment

Background

Although propofol is commonly used for general anaesthesia of normothermic patients in clinical practice, little information is available in the literature regarding the use of propofol anaesthesia for intracranial decompression using brain hypothermia treatment. A novel propofol anaesthesia scheme is proposed that should promote such clinical application and improve understanding of the principles of using propofol anaesthesia for hypothermic intracranial decompression.

Methods

Theoretical analysis was carried out using a previously-developed integrative model of the thermoregulatory, hemodynamic and pharmacokinetic subsystems. Propofol kinetics is described using a framework similar to that of this model and combined with the thermoregulation subsystem through the pharmacodynamic relationship between the blood propofol concentration and the thermoregulatory threshold. A propofol anaesthesia scheme for hypothermic intracranial decompression was simulated using the integrative model.

Results

Compared to the empirical anaesthesia scheme, the proposed anaesthesia scheme can reduce the required propofol dosage by more than 18%.

Conclusion

The integrative model of the thermoregulatory, hemodynamic and pharmacokinetic subsystems is effective in analyzing the use of propofol anaesthesia for hypothermic intracranial decompression. This propofol infusion scheme appears to be more appropriate for clinical application than the empirical one.

Physiological models of body composition and human obesity

Background

The body mass index (BMI) is the standard parameter for predicting body fat fraction and for classifying degrees of obesity. Currently available regression equations between BMI and fat are based on 2 or 3 parameter empirical fits and have not been validated for highly obese subjects. We attempt to develop regression relations that are based on realistic models of body composition changes in obesity. These models, if valid, can then be extrapolated to the high fat fraction of the morbidly obese.

Methods

The analysis was applied to 3 compartment (density and total body water) measurements of body fat. The data was collected at the New York Obesity Research Center, Body Composition Unit, as part of ongoing studies. A total of 1356 subjects were included, with a BMI range of 17 to 50 for males and 17 to 65 for females. The body composition model assumes that obese subjects can be represented by the sum of a standard lean reference subject plus an extra weight that has a constant adipose, bone and muscle fraction.

Results

There is marked age and sex dependence in the relationship between BMI and fat fraction. There was no significant difference among Caucasians, Blacks and Hispanics while Asians had significantly greater fat fraction for the same BMI. A linear relationship between BMI and fat fraction provides a good description for men but overestimates the fat fraction in morbidly obese women for whom a non-linear regression should be used. New regression relations for predicting body fat just from experimental measurements of body density are described that are more accurate then those currently used. From the fits to the experimental BMI and density data, a quantitative description of the bone, adipose and muscle body composition of lean and obese subjects is derived.

Conclusion

Physiologically realistic models of body composition provide both accurate regression relations and new insights about changes in body composition in obesity.

A model to predict patient temperature during cardiac surgery

A core temperature drop after cardiac surgery slows down the patient's recuperation process. In order to minimize the amount of the so-called afterdrop, more knowledge is needed about the impaired thermoregulatory system during anesthesia and the effect of different protocols on temperature distribution. Therefore, a computer model has been developed that describes heat transfer during cardiac surgery. The model consists of three parts: (1) a passive part, which gives a simplified description of the human geometry and the passive heat transfer processes, (2) an active part that takes into account the thermoregulatory system as a function of the amount of anesthesia and (3) submodels, through which it is possible to adjust the boundary conditions. The validity of the new model was tested by comparing the model results to the measurement results of three surgical procedures. A good resemblance was found between simulation results and the experiments. Next, a model application was shown. A parameter study was performed to study the effect of different temperature protocols on afterdrop. It was shown that the effectiveness of forced-air heating is larger than the benefits resulting from increased environmental temperature or usage of a circulating water mattress. Ultimately, the model could be used to develop a monitoring decision system that advises clinicians what temperature protocol will be best for the patient.

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.

Human hepatocytes: isolation, cryopreservation and applications in drug development

The recent developments in the isolation, culturing, and cryopreservation of human hepatocytes, and the application of the cells in drug development are reviewed. Recent advances include the improvement of cryopreservation procedures to allow cell attachment, thereby extending the use of the cells to assays that requires prolong culturing such as enzyme induction studies. Applications of human hepatocytes in drug development include the evaluation of metabolic stability, metabolite profiling and identification, drug-drug interaction potential, and hepatotoxic potential. The use of intact human hepatocytes, because of the complete, undisrupted metabolic pathways and cofactors, allows the development of data more relevant to humans in vivo than tissue fractions such as human liver microsomes. Incorporation of key in vivo factors with the intact hepatocytes in vitro may help predictive human in vivo drug properties. For instance, evaluation of drug metabolism and drug-drug interactions with intact human hepatocytes in 100% human serum may eliminate the need to determine in vivo intracellular concentrations for the extrapolation of in vitro data to in vivo. Co-culturing of hepatocytes and nonhepatic primary cells from other organs in the integrated discrete multiple organ co-culture (IdMOC) may allow the evaluation of multiple organ interactions in drug metabolism and drug toxicity. In conclusion, human hepatocytes represent a critical experimental model for drug development, allowing early evaluation of human drug properties to guide the design and selection of drug candidates with a high probability of clinical success.

Population-based analysis of methadone distribution and metabolism using an age-dependent physiologically based pharmacokinetic model

Limited pharmacokinetic (PK) and pharmacodynamic (PD) data are available to use in methadone dosing recommendations in pediatric patients for either opioid abstinence or analgesia. Considering the extreme inter-individual variability of absorption and metabolism of methadone, population-based PK would be useful to provide insight into the relationship between dose, blood concentrations, and clinical effects of methadone. To address this need, an age-dependent physiologically based pharmacokinetic (PBPK) model has been constructed to systematically study methadone metabolism and PK. The model will facilitate the design of cost-effective studies that will evaluate methadone PK and PD relationships, and may be useful to guide methadone dosing in children. The PBPK model, which includes whole-body multi-organ distribution, plasma protein binding, metabolism, and clearance, is parameterized based on a database of pediatric PK parameters and data collected from clinical experiments. The model is further tailored and verified based on PK data from individual adults, then scaled appropriately to apply to children aged 0-24 months. Based on measured variability in CYP3A enzyme expression levels and plasma orosomucoid (ORM2) concentrations, a Monte-Carlo-based simulation of methadone kinetics in a pediatric population was performed. The simulation predicts extreme variability in plasma concentrations and clearance kinetics for methadone in the pediatric population, based on standard dosing protocols. In addition, it is shown that when doses are designed for individuals based on prior protein expression information, inter-individual variability in methadone kinetics may be greatly reduced.

Human physiologically based pharmacokinetic model for ACE inhibitors: ramipril and ramiprilat

Background

The angiotensin-converting enzyme (ACE) inhibitors have complicated and poorly characterized pharmacokinetics. There are two binding sites per ACE (high affinity "C", lower affinity "N") that have sub-nanomolar affinities and dissociation rates of hours. Most inhibitors are given orally in a prodrug form that is systemically converted to the active form. This paper describes the first human physiologically based pharmacokinetic (PBPK) model of this drug class.

Methods

The model was applied to the experimental data of van Griensven et. al for the pharmacokinetics of ramiprilat and its prodrug ramipril. It describes the time course of the inhibition of the N and C ACE sites in plasma and the different tissues. The model includes: 1) two independent ACE binding sites; 2) non-equilibrium time dependent binding; 3) liver and kidney ramipril intracellular uptake, conversion to ramiprilat and extrusion from the cell; 4) intestinal ramipril absorption. The experimental in vitro ramiprilat/ACE binding kinetics at 4°C and 300 mM NaCl were assumed for most of the PBPK calculations. The model was incorporated into the freely distributed PBPK program PKQuest.

Results

The PBPK model provides an accurate description of the individual variation of the plasma ramipril and ramiprilat and the ramiprilat renal clearance following IV ramiprilat and IV and oral ramipril. Summary of model features: Less than 2% of total body ACE is in plasma; 35% of the oral dose is absorbed; 75% of the ramipril metabolism is hepatic and 25% of this is converted to systemic ramiprilat; 100% of renal ramipril metabolism is converted to systemic ramiprilat. The inhibition was long lasting, with 80% of the C site and 33% of the N site inhibited 24 hours following a 2.5 mg oral ramipril dose. The plasma ACE inhibition determined by the standard assay is significantly less than the true in vivo inhibition because of assay dilution.

Conclusion

If the in vitro plasma binding kinetics of the ACE inhibitor for the two binding sites are known, a unique PBPK model description of the Griensven et. al. experimental data can be obtained.

Pharmacokinetic and Pharmacodynamic Characteristics of??Medications Used for Moderate Sedation

The ability to deliver safe and effective moderate sedation is crucial to the ability to perform invasive procedures. Sedative drugs should have a quick onset of action, provide rapid and clear-headed recovery, and be easy to administer and monitor. A number of drugs have been demonstrated to provide effective sedation for outpatient procedures but since each agent has its own limitations, a thorough knowledge of the available drugs is required to choose the appropriate drug, dose and/or combination regimen for individual patients. Midazolam, propofol, ketamine and sevoflurane are the most frequently used agents, and all have a quick onset of action and rapid recovery. The primary drawback of midazolam is the potential for accumulation of the drug, which can result in prolonged sedation and a hangover effect. The anaesthetics propofol and sevoflurane have recently been used for sedation in procedures of short duration. Although effective, these agents require monitored anaesthesia care. Ketamine is an effective agent, particularly in children, but there is concern regarding emergence reactions. AQUAVAN injection (fospropofol disodium), a phosphorylated prodrug of propofol, is an investigational agent possessing a unique and distinct pharmacokinetic and pharmacodynamic profile. Compared with propofol emulsion, AQUAVAN is associated with a slightly longer time to peak effect and a more prolonged pharmacodynamic effect. Advances in the delivery of sedation, including the development of new sedative agents, have the potential to further improve the provision of moderate sedation for a variety of invasive procedures.