Search

The SARS-CoV-2 coronavirus, which causes COVID-19, uses a viral surface spike protein for host cell entry and the human cell-surface transmembrane serine protease, TMPRSS2, to process the spike protein. Camostat mesylate, an orally available and clinically used serine protease inhibitor, inhibits TMPRSS2, supporting clinical trials to investigate its use in COVID-19. A one-compartment pharmacokinetic (PK)/pharmacodynamic (PD) model for camostat and the active metabolite FOY-251 was developed, incorporating TMPRSS2 reversible covalent inhibition by FOY-251, and empirical equations linking TMPRSS2 inhibition of SARS-CoV-2 cell entry. The model predicts that 95% inhibition of TMPRSS2 is required for 50% inhibition of viral entry efficiency. For camostat 200 mg dosed four times daily, 90% inhibition of TMPRSS2 is predicted to occur but with only about 40% viral entry inhibition. For 3-fold higher camostat dosing, marginal improvement of viral entry rate inhibition, up to 54%, is predicted. Because respiratory tract viral load may be associated with negative outcome, even modestly reducing viral entry and respiratory tract viral load may reduce disease progression. This modeling also supports medicinal chemistry approaches to enhancing PK/PD and potency of the camostat molecule. IMPORTANCE Strategies to repurpose already-approved drugs for the treatment of COVID-19 has been attractive since the beginning of the pandemic. Camostat mesylate, a serine protease inhibitor approved in Japan for the treatment of acute exacerbations of chronic pancreatitis, inhibits TMPRSS1, a host cell surface serine protease essential for SARS-CoV-2 viral entry. In vitro experiments provided data suggesting that camostat might be effective in the treatment of COVID-19. Multiple clinical trials were planned to test the hypothesis that camostat would be beneficial for treating COVID-19 (for example, clinicaltrials.gov, NCT04353284). The present work used a one-compartment pharmacokinetic (PK)/pharmacodynamic (PD) mathematical model for camostat and the active metabolite FOY-251, incorporating TMPRSS2 reversible covalent inhibition by FOY-251, and empirical equations linking TMPRSS2 inhibition of SARS-CoV-2 cell entry. This work is valuable to guide further development of camostat mesylate and possible medicinal chemistry derivatives for the treatment of COVID-19.

To provide insight into pharmacological treatment of hyperuricemia we developed a semi-mechanistic, dynamical model of uric acid (UA) disposition in human. Our model represents the hyperuricemic state in terms of production of UA (rate, PUA), its renal filtration (glomerular filtration rate, GFR) and proximal tubular reabsorption (fractional excretion coefficient, FE). Model parameters were estimated using data from 9 Phase I studies of xanthine oxidase inhibitors (XOI) allopurinol and febuxostat and a novel uricosuric, the selective UA reabsorption inhibitor lesinurad, approved for use in combination with a XOI. The model was qualified for prediction of the effect of patients' GFR and FE on concentration of UA in serum (sUA) and UA excretion in urine and their response to drug treatment, using data from 2 Phase I and 4 Phase III studies of lesinurad. Percent reduction in sUA from baseline by a XOI is predicted to be independent of GFR, FE or PUA. Uricosurics are more effective in underexcreters of UA or patients with normal GFR. Co-administration of a XOI and an uricosuric agent should be considered for patients with high sUA first in the treatment algorithm of gout before uptitration of XOI. The XOI dose in combination with a uricosuric can be reduced compared to XOI alone for the same target sUA to the degree dependent on patient's GFR and FE. This exposure-response model of UA can be used to rationally select the best drug treatment option to lower elevated sUA in gout patients under differing pathophysiological situations., (© 2018 Astra Zeneca. Physiological Reports published by Wiley Periodicals, Inc. on behalf of The Physiological Society and the American Physiological Society.)

Background: Numerous oncology combination therapies involving modulators of the cancer immune cycle are being developed, yet quantitative simulation models predictive of outcome are lacking. We here present a model-based analysis of tumor size dynamics and immune markers, which integrates experimental data from multiple studies and provides a validated simulation framework predictive of biomarkers and anti-tumor response rates, for untested dosing sequences and schedules of combined radiation (RT) and anti PD-(L)1 therapies., Methods: A quantitative systems pharmacology model, which includes key elements of the cancer immunity cycle and the tumor microenvironment, tumor growth, as well as dose-exposure-target modulation features, was developed to reproduce experimental data of CT26 tumor size dynamics upon administration of RT and/or a pharmacological IO treatment such as an anti-PD-L1 agent. Variability in individual tumor size dynamics was taken into account using a mixed-effects model at the level of tumor-infiltrating T cell influx., Results: The model allowed for a detailed quantitative understanding of the synergistic kinetic effects underlying immune cell interactions as linked to tumor size modulation, under these treatments. The model showed that the ability of T cells to infiltrate tumor tissue is a primary determinant of variability in individual tumor size dynamics and tumor response. The model was further used as an in silico evaluation tool to quantitatively predict, prospectively, untested treatment combination schedules and sequences. We demonstrate that anti-PD-L1 administration prior to, or concurrently with RT reveal further synergistic effects, which, according to the model, may materialize due to more favorable dynamics between RT-induced immuno-modulation and reduced immuno-suppression of T cells through anti-PD-L1., Conclusions: This study provides quantitative mechanistic explanations of the links between RT and anti-tumor immune responses, and describes how optimized combinations and schedules of immunomodulation and radiation may tip the immune balance in favor of the host, sufficiently to lead to tumor shrinkage or rejection.

Sphingosine 1-phosphate (S1P) receptor agonists are associated with cardiovascular effects in humans. This study aims to develop a systems pharmacology model to identify the site of action (i.e., primary hemodynamic response variable) of S1P receptor agonists, and to predict, in a quantitative manner, the cardiovascular effects of novel S1P receptor agonists in vivo. The cardiovascular effects of once-daily fingolimod (0, 0.1, 0.3, 1, 3, and 10 mg/kg) and siponimod (3 and 15 mg/kg) were continuously recorded in spontaneously hypertensive rats and Wistar-Kyoto rats. The results were analyzed using a recently developed systems cardiovascular pharmacology model, i.e. the CVS model; total peripheral resistance and heart rate were identified as the site of action for fingolimod. Next, the CVS model was interfaced with an S1P agonist pharmacokinetic-pharmacodynamic (PKPD) model. This combined model adequately predicted, in a quantitative manner, the cardiovascular effects of siponimod using in vitro binding assays. In conclusion, the combined CVS and S1P agonist PKPD model adequately describes the hemodynamic effects of S1P receptor agonists in rats and constitutes a basis for the prediction, in a strictly quantitative manner, of the cardiovascular effects of novel S1P receptor agonists., (Copyright © 2017 by The American Society for Pharmacology and Experimental Therapeutics.)

Background and Purpose: Previously, a systems pharmacology model was developed characterizing drug effects on the interrelationship between mean arterial pressure (MAP), cardiac output (CO) and total peripheral resistance (TPR). The present investigation aims to (i) extend the previously developed model by parsing CO into heart rate (HR) and stroke volume (SV) and (ii) evaluate if the mechanism of action (MoA) of new compounds can be elucidated using only HR and MAP measurements., Experimental Approach: Cardiovascular effects of eight drugs with diverse MoAs (amiloride, amlodipine, atropine, enalapril, fasudil, hydrochlorothiazide, prazosin and propranolol) were characterized in spontaneously hypertensive rats (SHR) and normotensive Wistar-Kyoto (WKY) rats following single administrations of a range of doses. Rats were instrumented with ascending aortic flow probes and aortic catheters/radiotransmitters for continuous recording of MAP, HR and CO throughout the experiments. Data were analysed in conjunction with independent information on the time course of the drug concentration following a mechanism-based pharmacokinetic-pharmacodynamic modelling approach., Key Results: The extended model, which quantified changes in TPR, HR and SV with negative feedback through MAP, adequately described the cardiovascular effects of the drugs while accounting for circadian variations and handling effects., Conclusions and Implications: A systems pharmacology model characterizing the interrelationship between MAP, CO, HR, SV and TPR was obtained in hypertensive and normotensive rats. This extended model can quantify dynamic changes in the CVS and elucidate the MoA for novel compounds, with one site of action, using only HR and MAP measurements. Whether the model can be applied for compounds with a more complex MoA remains to be established., (© 2014 The British Pharmacological Society.)

A complicating factor in the translational pharmacology of sphingosine 1-phosphate agonists is that they exert their pharmacological effect through their respective phosphate metabolites, which are formed by the enzyme sphingosine kinase (S1PHK). In this investigation, we present a semimechanistic pharmacokinetic model for the interconversion of S1PHK substrates and their respective phosphates in rats and humans with the aim of investigating whether characterization of the rate of phosphorylation in blood platelets constitutes a basis for interspecies scaling using fingolimod as a model compound. Data on the time course of fingolimod and fingolimod-phosphate (fingolimod-P) blood concentrations after intravenous and oral administration of fingolimod and/or fingolimod-P in rats and after oral administration of fingolimod in doses of 0.5, 1.25, and 5 mg once daily in healthy volunteers were analyzed in conjunction with data on the ex vivo interconversion and blood-plasma distribution in rat and human blood, respectively. Integrating the data from the ex vivo and in vivo studies enabled simulation of fingolimod and fingolimod-P concentrations in plasma rather than blood, which are more relevant for characterizing drug effects. Large interspecies differences in the rate of phosphorylation between rats and humans were quantified. In human, phosphorylation of fingolimod in the platelets was four times slower compared with rat, whereas the dephosphorylation rates were comparable in both species. This partly explained the 10-12-fold overprediction of fingolimod-P exposure in human when applying a dose-by-factor approach on the developed rat model. Additionally, differences in presystemic phosphorylation should also be taken into account., (Copyright © 2014 by The American Society for Pharmacology and Experimental Therapeutics.)

Background and Purpose: The homeostatic control of arterial BP is well understood with changes in BP resulting from changes in cardiac output (CO) and/or total peripheral resistance (TPR). A mechanism-based and quantitative analysis of drug effects on this interrelationship could provide a basis for the prediction of drug effects on BP. Hence, we aimed to develop a mechanism-based pharmacokinetic-pharmacodynamic (PKPD) model in rats that could be used to characterize the effects of cardiovascular drugs with different mechanisms of action (MoA) on the interrelationship between BP, CO and TPR., Experimental Approach: The cardiovascular effects of six drugs with diverse MoA, (amlodipine, fasudil, enalapril, propranolol, hydrochlorothiazide and prazosin) were characterized in spontaneously hypertensive rats. The rats were chronically instrumented with ascending aortic flow probes and/or aortic catheters/radiotransmitters for continuous recording of CO and/or BP. Data were analysed in conjunction with independent information on the time course of drug concentration using a mechanism-based PKPD modelling approach., Key Results: By simultaneous analysis of the effects of six different compounds, the dynamics of the interrelationship between BP, CO and TPR were quantified. System-specific parameters could be distinguished from drug-specific parameters indicating that the model developed is drug-independent., Conclusions and Implications: A system-specific model characterizing the interrelationship between BP, CO and TPR was obtained, which can be used to quantify and predict the cardiovascular effects of a drug and to elucidate the MoA for novel compounds. Ultimately, the proposed PKPD model could be used to predict the effects of a particular drug on BP in humans based on preclinical data., (© 2013 The Authors. British Journal of Pharmacology © 2013 The British Pharmacological Society.)

The development of fundamental pharmacokinetics and pharmacodynamics concepts has enabled anesthesiologists to choose and dose anesthetic agents on a rational basis. The application of these concepts to a variety of clinical scenarios and patient populations makes it possible to individualize the dose, thereby decreasing the risk of complications. As more knowledge is gained about the sometimes profound differences in drug response, empirical dosing such as in milligrams per kilogram of total body weight is disappearing from the anesthesia specialty.

Background: Indacaterol is a once-daily long-acting inhaled β2-agonist indicated for maintenance treatment of moderate-to-severe chronic obstructive pulmonary disease (COPD). The large inter-patient and inter-study variability in forced expiratory volume in 1 second (FEV1) with bronchodilators makes determination of optimal doses difficult in conventional dose-ranging studies. We considered alternative methods of analysis., Methods: We utilized a novel modelling approach to provide a robust analysis of the bronchodilatory dose response to indacaterol. This involved pooled analysis of study-level data to characterize the bronchodilatory dose response, and nonlinear mixed-effects analysis of patient-level data to characterize the impact of baseline covariates., Results: The study-level analysis pooled summary statistics for each steady-state visit in 11 placebo-controlled studies. These study-level summaries encompassed data from 7476 patients at indacaterol doses of 18.75-600 μg once daily, and showed that doses of 75 μg and above achieved clinically important improvements in predicted trough FEV1 response. Indacaterol 75 μg achieved 74% of the maximum effect on trough FEV1, and exceeded the midpoint of the 100-140 mL range that represents the minimal clinically important difference (MCID; ≥120 mL vs placebo), with a 90% probability that the mean improvement vs placebo exceeded the MCID. Indacaterol 150 μg achieved 85% of the model-predicted maximum effect on trough FEV1 and was numerically superior to all comparators (99.9% probability of exceeding MCID). Indacaterol 300 μg was the lowest dose that achieved the model-predicted maximum trough response.The patient-level analysis included data from 1835 patients from two dose-ranging studies of indacaterol 18.75-600 μg once daily. This analysis provided a characterization of dose response consistent with the study-level analysis, and demonstrated that disease severity, as captured by baseline FEV1, significantly affects the dose response, indicating that patients with more severe COPD require higher doses to achieve optimal bronchodilation., Conclusions: Comprehensive assessment of the bronchodilatory dose response of indacaterol in COPD patients provided a robust confirmation that 75 μg is the minimum effective dose, and that 150 and 300 μg are expected to provide optimal bronchodilation, particularly in patients with severe disease.

Model-based drug development (MBDD) is a promising approach to improve decision making in drug development. The pharmaceutical industry has made substantial progress from engaging in empirical decision making to increasingly using pharmacometrics (ie, modeling and simulation [M&S]) as a quantitative decision-making tool. Focusing on culture and an organizational structure perspective, this commentary summarizes experiences and vision from industry M&S leaders on implementing MBDD. A culture for MBDD needs to have wide acceptance of MBDD, enhanced decision making with probability-based evidence and transparent rationale, quantitative impact metrics, and a brand that emphasizes cross-disciplinary collaboration and ownership. An organizational structure for MBDD needs to have a dedicated pharmacometrics function, fine balance between quick wins and impact on long-term R&D goals, and collaborative MBDD teams among clinical pharmacologists, statisticians, pharmacometricians, and clinicians. Pharmaceutical companies with these characteristics are prepared to fully embrace and implement MBDD.

In this chapter, drawn largely from the synthesis of material that we first presented in the sixth edition of Miller's Anesthesia, Chap 31 (Stanski and Shafer 2005; used by permission of the publisher), we have defined anesthetic depth as the probability of non-response to stimulation, calibrated against the strength of the stimulus, the difficulty of suppressing the response, and the drug-induced probability of non-responsiveness at defined effect site concentrations. This definition requires measurement of multiple different stimuli and responses at well-defined drug concentrations. There is no one stimulus and response measurement that will capture depth of anesthesia in a clinically or scientifically meaningful manner. The "clinical art" of anesthesia requires calibration of these observations of stimuli and responses (verbal responses, movement, tachycardia) against the dose and concentration of anesthetic drugs used to reduce the probability of response, constantly adjusting the administered dose to achieve the desired anesthetic depth. In our definition of "depth of anesthesia" we define the need for two components to create the anesthetic state: hypnosis created with drugs such as propofol or the inhalational anesthetics and analgesia created with the opioids or nitrous oxide. We demonstrate the scientific evidence that profound degrees of hypnosis in the absence of analgesia will not prevent the hemodynamic responses to profoundly noxious stimuli. Also, profound degrees of analgesia do not guarantee unconsciousness. However, the combination of hypnosis and analgesia suppresses hemodynamic response to noxious stimuli and guarantees unconsciousness.

Background and Objective: Enrichment strategies which select subjects who appear to respond to the drug have been used in drug studies to demonstrate clinical efficacy. We have used clinical trial simulation techniques to examine factors that are relevant in clinical trial design based on enrichment where poor responders are excluded from the double-blind phase of the study., Methods: Simulations were performed for an analgesic trial design involving an open-dose titration phase (enrichment phase) followed by a double-blind, randomized, placebo-controlled maintenance phase. Enrichment was examined by excluding subjects above a predefined pain score (cutoff) from analysis of efficacy for the maintenance phase. Cutoff pain scores ranging from 4 to 7 on a 0 to 10 categorical scale were examined. A database consisting of chronic pain patients who participated in studies with a new formulation of buprenorphine was used to build the simulation model. Since no data were available for the key model variable "correlation between treatment and placebo response", values of 0.25, 0.5, and 0.75 were used for the simulations., Results: A correlation between treatment and placebo effect ranging from 0.75 to 0.25 will cause the likelihood of trial success to vary from 50% to 95%. This model also shows that recruitment efficiency will decrease with the use of lower cutoff pain scores., Conclusion: Prior to using enrichment techniques, investigators must consider the correlation between treatment effect and placebo response to optimize clinical trial design.

Background: To compare the pharmacokinetics of intravenous, oral immediate-release (IR), and oral extended-release (OROS ) formulations of hydromorphone., Methods: In this randomized, six-session, crossover-design study, 12 subjects received hydromorphone 8-mg intravenous, 8-mg IR oral, and 8-, 16-, and 32-mg OROS formulations or placebo orally followed by plasma sampling for hydromorphone determination. Pharmacokinetic analysis was performed using NONMEM. Using the disposition of hydromorphone from the intravenous administration, deconvolution was used to estimate the input rate function (release rate from the gut to the blood) for the IR and OROS formulations. A linear spline was used to describe the drug input rate function., Results: The deconvolution using linear splines described the release characteristics of both the IR and OROS formulations. The mean absolute bioavailability for the 8-mg OROS formulation was significantly larger ( = 0.025) than for the 8-mg IR formulation: 0.24 (SD 0.059) versus 0.19 (SD 0.054), respectively. The bioavailability was the same for the three doses of the OROS formulation. Predicted degree of fluctuation of plasma concentrations would be expected to be 130% and 39% for the IR and OROS 8-mg doses, respectively., Conclusions: The OROS formulation of hydromorphone produced continued release of medication over 24 h, which should allow for once-daily oral dosing. The extended release of hydromorphone will produce less fluctuation of plasma concentrations compared with IR formulations, which should provide for more constant pain control. The in vivo release of hydromorphone from both IR and OROS formulations were adequately described using a linear spline deconvolution approach. The increased bioavailability from the OROS formulation may be related to decreased metabolism by a first-pass effect or enterohepatic recycling of hydromorphone.

Background: The disposition kinetics of hydromorphone generally necessitates oral administration every 4 h of the conventional immediate-release tablet to provide sustained pain relief. This trial examined time course and magnitude of analgesia to experimental pain after administration of sustained-release hydromorphone as compared with that after immediate-release hydromorphone or placebo., Methods: Using a 4 x 4 Latin square double-blind design, 12 subjects were randomized to receive a single dose of 8, 16, and 32 mg sustained-release hydromorphone and placebo. The same subjects had received 8 mg immediate-release hydromorphone before this study. Using an electrical experimental pain paradigm, analgesic effects were assessed for up to 30 h after administration, and venous hydromorphone plasma concentrations were measured at corresponding times., Results: The hydromorphone plasma concentration peaked significantly later (12.0 h [12.0--18.0] vs. 0.8 h [0.8--1.0]; median and interquartile range) but was maintained significantly longer at greater than 50% of peak concentration (22.7 +/- 8.2 h vs. 1.1 +/- 0.7 h; mean +/- SD) after sustained-release than after immediate-release hydromorphone. Similarly, sustained-release hydromorphone produced analgesic effects that peaked significantly later (9.0 h [9.0--12.0] vs. 1.5 h [1.0--2.0]) but were maintained significantly longer at greater than 50% of peak analgesic effect (13.3 +/- 6.3 h vs. 3.6 +/- 1.7 h). A statistically significant linear relation between the hydromorphone plasma concentration and the analgesic effect on painful stimuli existed., Conclusion: A single oral dose of a new sustained-release formulation of hydromorphone provided analgesia to experimental pain beyond 24 h of its administration.

Background: Epidural administration of morphine is a common analgesic technique to manage pain. Morphine spreads from the epidural space to the cerebrospinal fluid and then rostrally, causing side effects mediated by the brain stem. However, data on the rostral spread of morphine-mediated analgesia are sparse. This study examined the rostral spread of analgesic effects on heat and electrical pain after epidural administration of morphine., Methods: In a randomized, double-blinded, placebo-controlled, crossover study, 5 mg morphine or saline placebo were injected into the lumbar epidural space in nine healthy volunteers. Correct needle placement was confirmed with fluoroscopy. Analgesia to experimental nociceptive heat and electrical stimuli was measured at lumbar (L4), thoracic (T10), cervical (C2), and trigeminal (V2) levels before and 2, 5, 10, and 24 h after epidural injection. Plasma samples for assaying morphine concentrations were drawn before and after each analgesic evaluation., Results: Epidural morphine significantly attenuated experimental heat pain at all dermatomes tested compared with saline placebo. Analgesic effects were significant at L4 after 2, 5, and 10 h, at T10 after 5, 10, and 24 h, and at V2 after 10 h. Electrical pain was attenuated at the lumbar and thoracic but not at the cervical dermatome. Analgesic effects were significant at L4 after 2, 5, and 10 h and at T10 after 5 and 10 h. Morphine plasma concentrations were below the detection limit (1 ng/ml) in eight of the nine subjects 10 h after epidural injection., Conclusions: Lumbar epidural injection of morphine attenuated cutaneous heat pain up to the trigeminal dermatome during a 24-h observation period. In a clinical context, this implies that some types of pain may be attenuated up to the supraspinal level after lumbar epidural administration of morphine.

Background: The influence of changes in the physiologic state of a patient on the disposition of fentanyl and alfentanil is poorly understood. The aims of this study were to determine whether physiologic pharmacokinetic models for fentanyl and alfentanil, based on data from rats, could predict plasma concentrations of these opioids in humans and to determine how changes in physiology would influence the predictions of their disposition., Methods: The predictions of the models were tested against plasma concentration data from published pharmacokinetic studies. The influences of changes in body composition, cardiac output, and regional blood flows on the disposition of the opioids were simulated., Results: The models could predict independently measured plasma concentrations of the opioids after short infusions in humans. Simulations then predicted that differences in body composition between men and women would have little influence on the pharmacokinetics of the opioids. Changes in cardiac output would affect drug redistribution, and consequently the early decay of the plasma concentrations, but not markedly influence rates of elimination. Further, the clearance of the opioids would decrease and their volumes of distribution increase with the age of the patient, but this would only marginally affect the early disposition of the drugs. Even large fluctuations in peripheral or hepatic blood flows would have modest effects on arterial plasma concentrations of the opioids, and sudden "postoperative" increases in peripheral blood flows would cause minor secondary plasma concentration peaks., Conclusions: The ability of the physiologic models to predict plasma concentrations of fentanyl and alfentanil in humans was confirmed. When changes in physiologic condition were simulated, effects on the pharmacokinetics of the opioids with possible implications for dosing were obtained only if cardiac output was varied over a wide range.

This study characterizes the pharmacokinetic-pharmacodynamic (PK-PD) relationships of the cardiovascular, EEG, hypnotic and ventilatory effects of the alpha-2 adrenergic agonist dexmedetomidine in rats. Dexmedetomidine was administered by a single rapid infusion (n = 6) and by an infusion regimen of gradually increasing rate (n = 8). HR, mean arterial pressure (MAP) and EEG signals were recorded continuously, as was the time at which the rats woke up spontaneously from drug-induced sleep, a measure of hypnosis. Arterial concentrations of dexmedetomidine and blood gases were determined regularly. A sigmoidal Emax model was used to describe the HR, MAP and EEG concentration-effect relationships, with the EEG effect (activity in 0.5-3.5-Hz frequency band) linked to an effect-site model. The PK of dexmedetomidine could be described by a two-compartment model, with similar PK parameters for both infusion regimens. Plasma protein binding was 84.1[0.7]%. Because of complex cardiovascular homeostatic reflex mechanisms, HR and MAP could only be analyzed during gradually increasing infusions. The maximal decrease in HR was 35(2)%, and the maximal increase in MAP was 37(2)%. For both infusion regimens, similar PD parameters were found for the EEG and the hypnotic measure. These data suggest the absence of active metabolites or tolerance of the EEG and hypnotic effects. Judging on the basis of concentrations of dexmedetomidine (mean (S.E. M.)), HR decrease was the most sensitive response [EC50 of 0.65(0. 09) ng/ml], followed by increase in MAP [EC50 of 2.01(0.14) ng/ml], change in EEG activity [EC50 of 2.24(0.16) ng/ml] and the hypnotic measure [Cwake-up of 2.64(0.10) ng/ml]. Ventilatory effects were minor.

Background: Understanding the influence of physiological variables on thiopental pharmacokinetics would enhance the scientific basis for the clinical usage of this anesthetic., Methods: A physiological pharmacokinetic model for thiopental previously developed in rats was scaled to humans by substituting human values for tissue blood flows, tissue masses, and elimination clearance in place of respective rat values. The model was validated with published serum concentration data from 64 subjects. The model was simulated after intravenous thiopental administration, 250 mg, over 1 min, to predict arterial plasma concentrations under conditions of different cardiac outputs, degrees of obesity, gender, or age., Results: The human pharmacokinetic model is characterized by a steady state volume of distribution of 2.2 l/kg, an elimination clearance of 0.22 l/min, and a terminal half-life of 9 h. Measured thiopental concentrations are predicted with an accuracy of 6 +/- 37% (SD). Greater peak arterial concentrations are predicted in subjects with a low versus a high cardiac output (3.1 and 9.4 l/min), and in subjects who are lean versus obese (56 and 135 kg). Acutely, obesity influences concentrations because it affects cardiac output. Prolonged changes are due to differences in fat mass. Changes with gender and age are relatively minor., Conclusions: The physiological pharmacokinetic model developed in rats predicts thiopental pharmacokinetics in humans. Differences in basal cardiac output may explain much of the variability in early thiopental disposition between subjects.

With increasing human age (20-80 years), the electroencephalogram (EEG) dose requirement for the intravenous anesthetic thiopental decreases approximately 10% per decade of life. The goal of this study was to compare the dose required to attain isoelectric EEG in young (4-5 month) vs. aged (24-25-month) Fischer 344 rats. One second isoelectricity was found to be an endpoint where minimal cardiorespiratory depression occurred. The effects of age, infusion rate, and repeated administration were examined in nine young and nine old rodents. Thiopental dose requirement increased with increasing infusion rates. Repeated administration at two-day intervals did not demonstrate tolerance to thiopental. No difference in thiopental dose requirement was detected in the young vs. elderly rats. In a separate group of five young and five old rats, thiopental plasma, brain, heart, and CSF concentrations were measured when five seconds of EEG isoelectricity was achieved: no consistent differences were noted. The rat may not be an appropriate model to investigate acute age-related anesthetic effects in humans, because cardiovascular changes with age are dissimilar between species.

Background: The effects of anesthetic drugs on electroencephalograms (EEG) have been studied to develop the EEG as a measure of anesthetic depth. Bispectral analysis is a new quantitative technique that measures the consistency of the phase and power relationships and returns a single measure, the bispectral index. The purpose of this study was to compare the performance of the bispectral index, version 1.1, with other spectral analysis EEG measures of drug effect for three commonly used anesthetic drugs., Methods: The EEG waveforms from 31 adults receiving infusions of alfentanil, propofol, or midazolam were analyzed. The time course of spectral edge (SE95), relative power in delta band, and bispectral index were related to the estimated effect-site concentration with use of a sigmoidal Emax model to estimate the potency (IC50) and the plasma effect-site equilibration rate constant (Ke0) for each measure. The performance of the fitting was assessed by the coefficient of correlation between predicted and observed effect., Results: Alfentanil induced a high-amplitude low-frequency EEG response. Propofol induced a biphasic response. At low concentrations, both frequency and amplitude increased. When the concentration increased, the EEG slowed and the amplitude decreased. High concentration produced burst suppression. Midazolam increased EEG frequency and amplitude. Bispectral index, SE95, and delta power yield similar estimates of IC50 and ke0. Except for alfentanil, the performance of the modeling with the bispectral index was as good that with SE95 or delta power., Conclusion: Bispectral analysis can be used as a measure of the EEG effects of anesthetic drugs.

1. Oxycodone is a strong opioid agonist that is currently available in immediate-release (IR) formulations for the treatment of moderate to severe pain. Recently, controlled-release (CR) oxycodone tablets were developed to provide the benefits of twice-a-day dosing to patients treated with oxycodone. The purpose of this investigation was to develop and validate a pharmacokinetic model for CR oxycodone tablets in comparison with IR oxycodone solution. 2. Twenty-four normal male volunteers were enrolled in a single-dose, randomized, analytically blinded, two-way crossover study designed to compare the pharmacokinetics of two 10 mg CR oxycodone tablets with 20 mg IR oxycodone oral solution. Pharmacokinetic models describing the oxycodone plasma concentration vs time profiles of CR tablets and IR solution were derived using NONMEM version IV. The predictive performance of the models was assessed by comparison of predicted oxycodone plasma concentrations with actual oxycodone plasma concentrations observed in a separate group of 21 volunteers who received repeated doses of IR and CR oxycodone for 4 days. 3. The unit impulse disposition function of oxycodone was best described by a one-compartment model. Absorption rate of the IR solution was best described by a mono-exponential model with a lag time, whereas absorption rate of the CR tablet was best described using a bi-exponential model. The absorption profile of the CR tablets was characterized by a rapid absorption component (t1/2abs = 37 min) accounting for 38% of the available dose and a slow absorption phase (t1/2abs = 6.2 h) accounting for 62% of the available dose. Two 10 mg tablets of oral CR oxycodone hydrochloride were 102.7% bioavailable relative to 20 mg of IR oxycodone hydrochloride oral solution. The population model derived after administration of a single dose accurately predicted both the mean and range of oxycodone concentrations observed during 4 days of repeated dosing. The mean prediction error was 2.7% with a coefficient of variation of 54%. 4. The absorption characteristics of CR oxycodone tablets should allow effective plasma concentrations of oxycodone to be reached quickly and for effective concentrations to be maintained for a longer period after dosing compared with the IR oral solution. The CR dosage form has pharmacokinetic characteristics that permit 12 hourly dosing.

Objective: To derive a population pharmacokinetic-pharmacodynamic model that characterizes the distribution of pain relief scores and remedication times observed in patients receiving intramuscular ketorolac for the treatment of moderate to severe postoperative pain., Background: The data analysis approach deals with the complexities of analyzing analgesic trial data: (1) repeated measurements, (2) ordered categorical response variables, and (3) nonrandom censoring because the patients can take a rescue medication if their pain relief is insufficient., Methods: Patients (n = 522) received a single oral or intramuscular administration of placebo or a single intramuscular dose of 10, 30, 60, or 90 mg ketorolac for postoperative pain relief. Pain relief was measured periodically with use of a five-category ordinal scale up to 6 hours after dosing. In this period, 288 patients received additional medication because of insufficient pain relief. Pharmacokinetic data was available for 85 subjects. Models were fitted to the data with the NONMEM program., Results: The pharmacokinetic data was best described by a two-compartment model with first-order absorption. Pain relief was found to be a function of drug concentration (Emax model), time (waxing and waning of placebo effect), and an individual random effect. The drug concentration at half-maximal effect (EC50) and the first-order rate constant (keo) half-life for pain relief were 0.37 mg/L and 24 minutes. The probability of remedication was found to be a function of the observed level of pain relief and was found to increase with time. Monte Carlo simulations showed that adequate pain relief was achieved in 50% of the patients at 41, 27, 23, and 21 minutes after 10, 30, 60, or 90 mg of intramuscular ketorolac. Adequate pain relief was maintained up to 6 hours in 50%, 70%, 78%, and 81% of patients after these four doses. Only 25% of the patients achieved adequate pain relief with placebo., Conclusions: A population pharmacokinetic-pharmacodynamic model for the analgesic efficacy of intramuscular ketorolac was derived. The simulated relationship between dose, time, and percentage of patients with adequate pain relief suggested that 30 mg intramuscular ketorolac was the optimal initial dose for postoperative pain relief.

Background: Thiopental produces a concentration-dependent continuum of effects on the cortical electroencephalogram (EEG) that has been linked to behavioral measures of anesthetic depth. The complexity of the response, however, limits a clear insight into the neurophysiologic actions of thiopental. The current study investigated thiopental actions on cortical EEG and hippocampal electrical activity, to determine whether similar effects occur on both structures and to compare synchronized activity between these structures., Methods: Thiopental was administered intravenously via an implanted catheter in freely moving rats. Arterial blood oxygen/carbon dioxide concentration, thiopental concentrations, and temperature were monitored and controlled. Neocortical EEG was recorded from implanted dural surface electrodes and hippocampal neuron electrical activity was recorded from stereotaxically placed microelectrodes. Pharmacokinetic models were used to determine effect site concentrations., Results: Thiopental produced an increase in EEG frequency and amplitude at low concentrations (15-20 micrograms/ml total plasma, approximately 10 microM unbound), which produced a loss of righting reflex. This was followed by a frequency decrease and burst suppression activity at higher concentrations (50-80 micrograms/ml, approximately 60 microM), which produced a loss of tail pinch and corneal reflexes. Higher concentrations of thiopental ( > 60 micrograms/ml) uncoupled synchronized burst discharges recorded in hippocampus and cortex. Isoelectric EEG activity was associated with concentrations of 70-90 micrograms/ml (approximately 80 microM) and a deep level of anesthesia; motor reflexes were abolished, although cardiovascular reflexes remained. In all frequency bands, similar concentration-EEG effect relationships were observed for cortical and hippocampal signals, only differing in the magnitude of response. A reversed progression of effects was observed on recovery., Conclusions: The results confirm earlier findings in humans and animals and demonstrate that both the hippocampus and neocortex exhibit burst suppression and isoelectric activity during thiopental anesthesia. Thiopental-induced synchronized burst activity was depressed by progressively higher concentrations. The lost synchronization suggests a depression of synaptic coupling between cortical structures contributes to anesthesia.

Purpose: This study investigates the rate and extent of absorption following intramuscular injection of midazolam and diazepam., Methods: Four healthy male volunteers were recruited in this randomized three-way cross-over study. On one occasion each subject received simultaneous im injections of 5 mg midazolam and 10 mg diazepam in separate deltoid muscles. On two other separate occasions each subject received an iv infusion of 7.5 mg midazolam and 30 mg diazepam over five minutes. Frequent arterial blood samples were collected for up to two hours and venous blood samples were collected for up to 24 hours for midazolam and ten days for diazepam. A gas chromatography assay was used to determine the plasma concentrations of midazolam and diazepam. The im absorption profiles were estimated using constrained least-squares deconvolution., Results: There were substantial intersubject variabilities in the estimated pharmacokinetic parameters (volume and clearances) of intravenous midazolam and diazepam. The mean (+/-sd) time to peak plasma concentration (Cmax) was shorter for im midazolam (17.5 +/- 6.5 min) relative to diazepam (33.8 +/- 7.5 min). The mean (+sd) time to peak absorption rate was also shorter for midazolam (9 +/- 2 vs 13.8 +/- 7.5 min). The peak rate of absorption was identical (0.18 mg. min-1) and bioavailability was 1.0 for both drugs., Conclusions: We conclude that midazolam has more rapid absorption than diazepam following im administration.

This paper describes a sensitive HPLC-electrochemical detection analytical method for determining the concentration of the intravenous anesthetic, propofol, in human or rat plasma or serum and a variety of rat tissues. Internal standard and drug are extracted from serum or plasma and other tissues with pentane. 2,6-tert.-Butylmethylphenol is used as internal standard. It includes a novel steam distillation procedure for separating the highly lipophilic propofol from skin and fat. The plasma/serum assay has a precision of 1-4% (C.V.) in the range 10 ng/ml to 1 microgram/ml and permits the assay of assay of 5 ng/ml from 0.1 ml of plasma/serum. The tissue procedure allows the estimation of 50 ng/g in 0.1 g of tissue for most of the major organs with less than 2% (C.V.) precision. This assay was used to measure propofol concentrations in plasma/serum and tissue samples in support of a project to develop a physiological pharmacokinetic model for propofol in the rat.

Background: The goal of this investigation was to characterize the effects of thiopental on cardia output and regional blood flows in the rat. Blood flows influence thiopental pharmacokinetics. Acquisition of these data may ultimately permit evaluation of the contribution of thiopental-induced alterations in regional blood flows to the disposition and hypnotic effect of this drug., Methods: Chronically instrumented unrestrained Wistar rats (n=20) aged 3-4 months received either a dose of thiopental sufficient to induce a brief period of unconsciousness (20 mg.kg(-1)) or a larger dose achieving electroencephalographic burst suppression (45 mg.kg(-1)). Cardiac output and blood flows to 14 tissues were determined at 4 times in each rat for a period of 420 min using injections of radioactive microspheres (expressed as mean +/- SD). Mean arterial pressure, heart rate, and blood gas tensions were determined at all measurement times. Arterial plasma concentrations were sampled at postinfusion times., Results: No important changes in systemic cardiovascular measurements were detected after the smaller dose of thiopental. One minute after the larger dose, cardiac output decreased from baseline (123 +/- 14 to 84 +/- ml.min (-1), P< 0.01), flow to muscle and fat decreased, and muscle and fat resistance increased. At 5 min, compared to baseline, no difference in cardiac output was detected (123 +/- vs. 119 +/- ml.min (-1)), intestinal flows increased and intestinal resistances decreased. Cardiac output was again depressed at 30, 90, and 180 min. Brain blood flow decreased 25 +/- 19 % (P< 0.01) from baseline for the duration of the study., Conclusions: Thiopental acutely decreases cardiac output, and blood flows to muscle and fat tissue. The temporary return of cardiac output to baseline may be related to intestinal vasodilation. These blood flow alterations may influence the pharmacokinetics of thiopental.

Background: In contrast to that of inhalational anesthetics, quantitation of anesthetic depth for intravenous agents has not been well defined. In this study, using rodents, the relationship between the constant plasma thiopental concentrations and the clinical response to multiple nociceptive stimuli were investigated characterizing the anesthetic state from light sedation to deep anesthesia and correlated to the degree of electroencephalogram (EEG) drug effect., Methods: Thirty rats were instrumented with chronically implanted EEG electrodes, arterial and venous catheters. A computer-driven infusion pump was used to rapidly attain and then maintain constant, target plasma thiopental concentrations ranging from 7 to 100 micrograms/ml. Three different target plasma thiopental concentrations were achieved in each rat. Electroencephalographic effects were monitored with aperiodic waveform analysis. The following nociceptive stimuli were applied: (1) unprovoked righting reflex, (2) provoked righting reflex, (3) noise stimulus, (4) tail clamping with an alligator clip, (5) constant tail pressure with an analgesiameter, (6) corneal reflex, and (7) tracheal intubation. For tail clamping, tail pressure, and intubation, either purposeful extremity movement or abdominal muscle contraction response was noted to be present or absent. The clinical responses (present or absent) were modeled using logistic regression to estimate the Cp50, the plasma thiopental concentration with a 50% probability of no response., Results: The following mean Cp50 values (95% confidence interval) were obtained: unprovoked righting reflex, 15.9 (15.1-16.6) micrograms/ml; provoked righting reflex, 21.4 (20.2-22.7) micrograms/ml; noise stimuli, 31.3 (29.7-33.0) micrograms/ml; tail clamp and limb movement, 38.3 (36.1-40.4) micrograms/ml; tail pressure and limb movement, 39.2 (37.1-41.3) micrograms/ml; tail pressure and abdominal muscle contraction, 52.5 (50.0- 55) micrograms/ml; tail clamping and abdominal muscle contraction, 56.1 (50.0-56.2) micrograms/ml; corneal reflex, 60.0 (56.6-63.4) micrograms/ml; and limb movement or muscle abdominal contraction response to intubation, 67.7 (59.2-76.1) micrograms/ml. At an EEG-effect of 9.1 and 2.2 waves/s, there was a 50% chance of limb movement response to tail clamping and tracheal intubation, respectively. There was a poor relationship between the plasma thiopental concentration and the percent increase of either heart rate or mean arterial blood pressure after applying either tail pressure or tail clamp stimuli., Conclusions: A range of nociceptive stimuli and their observed clinical responses can be used to quantitate thiopental anesthetic depth, ranging from light sedation to deep anesthesia (isoelectric EEG and unresponsive to intubation) in the rodent. Clinical response can be mapped to surrogate EEG measures.

Background: The pharmacokinetic profiles of sufentanil available in the literature are conflicting because of methodologic differences. Length of sampling and assay sensitivity are key factors involved in accurately estimating the volumes of distribution, clearances, and elimination phase. The unit disposition function of increasing doses of sufentanil were investigated and the influence of dose administered on the linearity of pharmacokinetics was assessed., Methods: The pharmacokinetics of sufentanil were investigated in 23 patients, aged 14-68 yr, scheduled for surgery with postoperative ventilation. After induction of anesthesia, sufentanil was administered as a short infusion (10-20 min) in doses ranging from 250 micrograms to 1,500 micrograms. Frequent arterial blood samples were gathered during and at the end of infusion, then at specific intervals up to 48 h after infusion. Plasma concentrations of sufentanil were measured by radioimmunoassay (limit of sensitivity 0.02 ng.ml-1). The data were analyzed with the standard two-stage, naive pooled-data and the mixed effect pharmacokinetic approaches., Results: The pharmacokinetics of sufentanil were adequately described by a linear three-compartmental mamillary model with the following parameters, expressed as log mean values with 95% confidence intervals: the central volume of distribution = 14.3 l (13.1-15.41), the rapidly equilibrating volume = 63.1 l (61.9-64.3 l), the slowly equilibrating volume = 261.6 l (260.2-262.9 l), the steady-state distribution volume = 339 l (335-343 l), metabolic clearance = 0.92 l.min-1 (0.84-1.05 l.min-1), rapid distribution clearance = 1.55 l.min-1 (1.34-2.14 l.min-1), slow distribution clearance = 0.33 l.min-1 (0.27-0.49 l.min-1), and elimination half-life = 769 min (690-1011 min). No relation to age, weight, or lean body mass was found for any of the parameters., Conclusions: Sufentanil pharmacokinetics were linear within the dose range studied. Drug detection up to 24 h after dosing was necessary to define the terminal elimination phase. The metabolic clearance approached liver blood flow and a large volume of distribution was identified, consistent with the long terminal elimination half-life. Simulations predicted that plasma sufentanil steady-state concentrations would rapidly decline after termination of an infusion despite the long half-lives.

The purpose of this study was to model pharmacodynamically the reversal of midazolam sedation with flumazenil. Ten human volunteers underwent four different sessions. In session 1, individual midazolam pharmacokinetics and electroencephalographic pharmacodynamics were determined. In sessions 2 and 3, a computer-controlled infusion of midazolam with individual volunteer pharmacokinetic data was administered, targeting a plasma concentration corresponding to a light or deep level of sedation (20% or 80% of the maximal midazolam electroencephalographic effect) for a period of 210 minutes. After obtaining a stable electroencephalographic effect and constant midazolam plasma concentrations, a zero-order infusion of flumazenil was started until complete reversal of midazolam electroencephalographic effect was obtained. The flumazenil infusion was then stopped and the volunteer was allowed to resedate because of the constant midazolam drug effect. The electroencephalographic response was measured during a 180-minute period and analyzed by aperiodic analysis and fast-Fourier transforms. In session 4, a midazolam plasma concentration corresponding to a deep level of sedation was targeted for 210 minutes to examine for the possible development of acute tolerance. No flumazenil was given in session 4. For a light sedation level, with a mean midazolam plasma concentration of 160 +/- 64 ng/ml, the mean half-life of the equilibration rate constant of flumazenil reversal is 5.0 +/- 2.5 minutes, and the mean effect site concentration causing 50% of Emax is 13.7 +/- 5.8 ng/ml. For a deep level of sedation, with a mean midazolam plasma concentration of 551 +/- 196 ng/ml, the mean half-life of the equilibration rate constant is 3.9 +/- 1.5 minutes, and the mean effect site concentration causing 50% of Emax is 20.6 +/- 6.8 ng/ml. This study provides an estimate of the magnitude of the blood/central nervous system equilibration delay for flumazenil antagonism of midazolam sedation and further defines the usefulness of the electroencephalogram as a measure of midazolam pharmacodynamic effect.

The safety, pharmacokinetics, and pharmacodynamics of the investigational partial opioid agonist, mirfentanil, were determined in a dose-escalating, Phase 1 study in healthy male volunteers. Hemodynamic, central nervous system, and respiratory monitoring were used for safety assessment. The electroencephalogram (EEG) was evaluated as a surrogate measure of drug effect. Butorphanol was chosen as the control drug. In the mirfentanil group (n = 8) the dose was increased in sequential subjects from 25 micrograms.kg-1.min-1 for 30 min to 450 micrograms.kg-1.min-1 for 15 min, and in the butorphanol group (n = 10) from 2 micrograms.kg-1.min-1 for 30 min to 25 micrograms.kg-1.min-1 for 15 min. In the mirfentanil group, serious side effects were observed at plasma concentrations more than 2000 ng/mL: heart rates exceeded 130 bpm (n = 2), epileptiform EEG potentials (n = 2), and a convulsion (n = 1). The clearance of mirfentanil was high (5.8-7.2 L/min), and the volume of distribution large (247-348 L). The EEG of the subjects receiving mirfentanil showed no changes typical for opioids. Butorphanol however, caused intermittent slowing in the delta and theta ranges. The results of our study define the upper limit of safe plasma concentrations in future mirfentanil studies.

Background: When given as an intravenous bolus for induction of anesthesia, propofol can decrease postintubation hypertension but can also create moderate to severe postinduction, preintubation hypotension. The addition of fentanyl usually decreases the postintubation hypertension but can increase the propofol-induced preintubation hypotension. The goal of the study was to determine the relation between propofol and fentanyl doses and the hemodynamic changes post-induction, preintubation and postintubation., Methods: Twelve groups of 10 patients, ASA physical status 1 or 2, first received fentanyl 0, 2, or 4 micrograms.kg-1 and then 5 min later received propofol 2.0, 2.5, 3.0, or 3.5 mg.kg-1 as an intravenous bolus for induction of anesthesia. Arterial blood pressure was continuously monitored. The trachea was intubated 4 min after propofol administration., Results: The mean decrease in systolic blood pressure after propofol was 28 mmHg when no fentanyl was given, 53 mmHg after 2 microgram.kg-1 of fentanyl (P < 0.05 vs. no fentanyl), and 50 mmHg after 4 micrograms.kg-1 (P < 0.05 vs. no fentanyl; no statistically significant difference 4 vs. 2 micrograms.kg-1). There was no statistically significant difference in hemodynamic response to intubation relative to propofol dose. Hemodynamic response to intubation was decreased by the administration of fentanyl; the mean increase of systolic blood pressure after intubation was 65 mmHg from preintubation value without fentanyl, 50 mmHg after 2 micrograms.kg-1, and 37 mmHg after 4 micrograms.kg-1 (P < 0.05 for 2 and 4 micrograms.kg-1 vs. no fentanyl and for 4 vs. 2 micrograms.kg-1). Hemodynamic changes postintubation were not statistically different with increasing doses of propofol., Conclusions: Hemodynamic changes after induction with propofol or propofol/fentanyl, pre- or postintubation, are not modified when the propofol dose is increased from 2 to 3.5 mg.kg-1. Maximal hypotension preintubation occurs with a fentanyl dose of 2 micrograms.kg-1, whereas the magnitude of postintubation hypertension is significantly decreased with an increase in the fentanyl dose to 4 micrograms.kg-1.

The objectives of this investigation were to characterize the disposition of fentanyl and alfentanil in 14 tissues in the rat, and to create physiological pharmacokinetic models for these opioids that would be scalable to man. We first created a parametric submodel for the disposition of either drug in each tissue and then assembled these submodels into whole-body models. The disposition of fentanyl and alfentanil in the heart and brain and of fentanyl in the lungs could be described by perfusion-limited 1-compartment models. The disposition of both opioids in all other examined tissues was characterized by 2- or 3-compartment models. From these models, the extraction ratios of the opioids in the various tissues could be calculated, confirming the generally lower extraction of alfentanil as compared to fentanyl. Assembly of the single-tissue models resulted in a wholebody model for fentanyl that accurately described its disposition in the rat. A similar assembly of the tissue models for alfentanil revealed non-first-order elimination kinetics that were not apparent in the blood concentration data. Michaelis-Menten parameters for the hepatic metabolism of alfentanil were determined by iterative optimization of the entire model. The parametric models were finally scaled to describe the disposition of fentanyl and alfentanil in humans.

Objective: We determined the possible benefits of a new opioid, trefentanil, relative to fentanyl and alfentanil using high-resolution pharmacokinetic-pharmacodynamic modeling and computer simulations of clinical dosing scenarios., Methods: First, we determined in nine volunteers the electroencephalographic (EEG) effects and the trefentanil infusion rate that gave maximal EEG changes in 3 to 10 minutes. Then, in a crossover fashion in five volunteers, we compared the pharmacokinetics and EEG pharmacodynamics of trefentanil with fentanyl and alfentanil. Finally, we used computer simulations to predict offset of opioid effects of trefentanil, fentanyl, and alfentanil when given in different dosing schemes., Results: The pharmacokinetic-pharmacodynamic profile of trefentanil was similar to alfentanil, except for a higher elimination clearance. Trefentanil versus alfentanil pharmacokinetic parameters were as follows: Elimination clearance, 0.444 +/- 0.073 versus 0.184 +/- 0.031 L/min; steady-state distribution volume, 37 +/- 7 versus 23 +/- 3 L; and elimination half-life, 127 +/- 24 versus 114 +/- 19 minutes. Trefentanil versus alfentanil pharmacodynamics were as follows: the equilibration half-time between EEG effect and arterial drug concentration, 1.2 +/- 0.5 versus 0.6 +/- 0.4 minutes; and the concentration resulting in 50% of maximal EEG effect, 429 +/- 313 versus 577 +/- 273 ng/ml. The pharmacokinetic-pharmacodynamic profile of fentanyl was significantly different from trefentanil and alfentanil. Simulation of effect compartment concentration decay curves after variable-length infusions predicted more rapid recovery from trefentanil than from alfentanil or fentanyl., Conclusion: We suggest that high-resolution pharmacokinetic-pharmacodynamic studies and computer simulations of clinical dosing scenarios may have significant usefulness in appreciating differences between new and established drugs in early phase I studies.

Physiologically based pharmacokinetic modeling procedures employ anatomical tissue weight, blood flow, and steady tissue/blood partition data, often obtained from different sources, to construct a system of differential equations that predict blood and tissue concentrations. Because the system of equations and the number of variables optimized is considerable, physiologic modeling frequently remains a simulation activity where fits to the data are adjusted by eye rather than with a computer-driven optimization algorithm. We propose a new approach to physiological modeling in which we characterize drug disposition in each tissue separately using constrained numerical deconvolution. This technique takes advantage of the fact that the drug concentration time course, CT(t), in a given tissue can be described as the convolution of an input function with the unit disposition function (UDFT) of the drug in the tissue, (i.e., CT(t) = (Ca(t)QT)*UDFT(t) where Ca(t) is the arterial concentration, Q tau is the tissue blood flow and * is the convolution operator). The obtained tissue until disposition function (UDF) for each tissue describes the theoretical disposition of a unit amount of drug infected into the tissue in the absence of recirculation. From the UDF, a parametric model for the intratissue disposition of each tissue can be postulated. Using as input the product of arterial concentration and blood flow, this submodel is fit separately utilizing standard nonlinear regression programs. In a separate step, the entire body is characterized by reassembly of the individuals submodels. Unlike classical physiologic modeling the fit for a given tissue is not dependent on the estimates obtained for other tissues in the model. Additionally, because this method permits examination of individual UDFs, appropriate submodel selection is driven by relevant information. This paper reports our experience with a piecewise modeling approach for thiopental disposition in the rat.

Background: alpha 2-Adrenergic agonists such as dexmedetomidine can be used to reduce the dose requirement of intravenous and volatile anesthetics. Whereas dexmedetomidine and volatile anesthetics interact pharmacodynamically (reduction of MAC), the mechanism of interaction between dexmedetomidine and intravenous anesthetics is not known., Methods: Fourteen male ASA physical status 1 patients were randomly assigned to serve as control subjects (n = 7) or to be treated with dexmedetomidine (n = 7; 100, 30, and 6 ng.kg-1.min-1 for 10 min, 15 min, and thereafter, respectively). After 35 min, in all patients, thiopental (100 mg/min) was infused until burst suppression appeared in the raw tracing of the electroencephalogram. By using concentrations of thiopental in plasma and the electroencephalogram as a continuous pharmacologic effect measure, the apparent effect site concentrations for thiopental were estimated in both groups. Three-compartment pharmacokinetics were calculated for thiopental., Results: Dexmedetomidine reduced the thiopental dose requirement for electroencephalographic burst suppression by 30%. There was no difference in estimated thiopental effect site concentrations between dexmedetomidine and control patients, suggesting the absence of a major pharmacodynamic interaction. Dexmedetomidine significantly decreased distribution volumes (V2, V3, and Vdss) and distribution clearances (Cl12 and Cl13) of thiopental., Conclusions: The thiopental dose-sparing effect of dexmedetomidine on the electroencephalogram is not the result of a pharmacodynamic interaction but rather can be explained by a dexmedetomidine-induced decrease in thiopental distribution volume and distribution clearances. Dexmedetomidine reduces thiopental distribution, most probably by decreasing cardiac output and regional blood flow.

Background: The terminal pharmacokinetic parameters of sufentanil have, until now, been poorly characterized. This is probably because of the poor sensitivity or unreliability of the assay methods used. Radioimmunoassay (RIA) can be a very helpful assay method for sufentanil. However, before application to key pharmacokinetic studies, it requires adequate validation, e.g., by comparison with a method of proven sensitivity and specificity, such as gas chromatography-mass spectrometry (GC-MS)., Methods: Spiked control plasma samples and 135 plasma samples obtained from five patients receiving intravenous doses of 500 or 750 micrograms sufentanil, as a 10-20-min infusion, were analyzed by an improved, sensitive RIA and capillary GC-MS., Results: Both techniques had comparable limits of quantitation (0.02 ng/ml). Between-day coefficients of variation in the 0.05-10-ng/ml concentration range were 8.5-10.5% for the RIA and less than 10% for the GC-MS method. The patient plasma concentrations determined by RIA (y) and GC-MS (x) showed a good agreement (y = 1.016x + 0.002) and a correlation coefficient of 0.97., Conclusions: The results demonstrate the validity of the improved RIA method for the determination of sufentanil plasma concentrations.

Background: Remifentanil (GI87084B) is a new short-acting opioid with a unique ester structure. Metabolism of remifentanil by ester hydrolysis results in very rapid elimination. The aim of this study was to characterize in detail the pharmacokinetic profile of remifentanil in healthy male volunteers., Methods: Ten healthy adult male volunteers received a zero-order infusion of remifentanil at doses ranging from 1 to 8 micrograms.kg-1.min-1 for 20 min. Frequent arterial blood samples were drawn and analyzed by gas chromatographic mass spectroscopy to determine the remifentanil blood concentrations. The raw pharmacokinetic data were analyzed using three different parametric compartmental modeling methods (traditional two-stage, naive pooled data, and NONMEM). The raw pharmacokinetic data also were analyzed using numeric deconvolution and a nonparametric moment technique. A computer simulation using hte pharmacokinetic parameters of the NONMEM compartmental model was performed to provide a more intuitively meaningful and clinically relevant description of the pharmacokinetics. The simulation estimated the time necessary to achieve a 50% decrease in remifentanil concentration after a variable-length infusion., Results: For each parametric method, a three-compartment mamillary model that accurately describes remifentanil's concentration decay curve was constructed. The NONMEM analysis population pharmacokinetic parameters included a central clearance of 2.8 l/min, a volume of distribution at steady state of 32.8 l, and a terminal half-life of 48 min. The mean results of the nonparametric moment analysis included a clearance of 2.9 l/min, a volume of distribution at steady state of 31.8 l, and a mean residence time of 10.9 min. The computer simulation revealed the strikingly unique pharmacokinetic profile of remifentanil compared to that of the currently available fentanyl family of opioids., Conclusions: Remifentanil is a new, short-acting opioid with promising clinical potential in anesthesiology.

To determine if the pharmacokinetics of alfentanil varied with the type of surgery (site of operation and intensity of surgical stimulation) or with the mode of administration, three groups of patients were given a variable-rate infusion of alfentanil (supplemented with bolus injections) as an adjunct to nitrous oxide-oxygen anaesthesia for breast surgery, lower abdominal surgery, or upper abdominal surgery. A two-compartment pharmacokinetic model adequately characterized the plasma concentration profiles of alfentanil for 31 of 34 patients. Mean distribution half-life (6 min), elimination half-life (87 min), initial volume of distribution (10 litres), steady state volume of distribution (31 litres), total body clearance (308 ml min-1), and distribution clearance (656 ml min-1) did not vary with the type of surgery. The pharmacokinetics could not be determined for three patients undergoing upper abdominal surgery, probably because of nonstationarity. Comparison of our data with those obtained after a single bolus injection or a constant-rate infusion indicates that the pharmacokinetics of alfentanil do not vary with the mode of administration.

Traditionally, physiological pharmacokinetic models assume that arterial blood flow to tissue is the rate-limiting step in the transfer of drug into tissue parenchyma. When this assumption is made the tissue can be described as a well-stirred single compartment. This study presents the tissue washout concentration curves of the two opioid analgesics fentanyl and alfentanil after simultaneous 1-min iv infusions in the rat and explores the feasibility of characterizing their tissue pharmacokinetics, modeling each of the 12 tissues separately, by means of either a one-compartment model or a unit disposition function. The tissue and blood concentrations of the two opioids were measured by gas-liquid chromatography. The well-stirred one-compartment tissue model could reasonably predict the concentration-time course of fentanyl in the heart, pancreas, testes, muscle, and fat, and of alfentanil in the brain and heart only. In most other tissues, the initial uptake of the opioids was considerably lower than predicted by this model. The unit disposition functions of the opioids in each tissue could be estimated by nonparametric numerical deconvolution, using the arterial concentration times tissue blood flow as the input and measured tissue concentrations as the response function. The observed zero-time intercepts of the unit disposition functions were below the theoretical value of one, and were invariably lower for alfentanil than for fentanyl. These findings can be explained by the existence of diffusion barriers within the tissues and they also indicate that alfentanil is less efficiently extracted by the tissue parenchyma than the more lipophilic compound fentanyl. The individual unit disposition functions obtained for fentanyl and alfentanil in 12 rat tissues provide a starting point for the development of models of intratissue kinetics of these opioids. These submodels can then be assembled into full physiological models of drug disposition.