Population pharmacokinetics/pharmacodynamics of anesthetics

In this article we review how population pharmacokinetic/pharmacodynamic (PD) modeling has evolved in the specialty of anesthesiology, how anesthesiology benefited from the mixed-effects approach, and which features of modeling need careful attention. Key articles from the anesthesiology literature are selected to discuss the modeling of typical anesthesiological PD end points, such as level of consciousness and analgesia, interactions between hypnotics and analgesics, estimation with poor and sometimes rich data sets from populations of various sizes, covariate detection, covariances between random effects, and Bayesian forecasting.     

Population pharmacokinetics and pharmacodynamics of oral etoposide

AIMS: To study the population pharmacokinetics and pharmacodynamics of oral etoposide in patients with solid tumours. METHODS: A prospective, open label, cross-over, bioavailability study was performed in 50 adult patients with miscellaneous, advanced stage solid tumours, who were receiving oral (100 mg capsules) etoposide for 14 days and i.v. (50 mg) etoposide on day 1 or day 7 in randomised order during the first cycle treatment. Total and unbound etoposide concentration were assayed by h.p.l.c. Population PK parameters estimation was done by using the P-Pharm software (Simed). Haematological toxicity and tumour response were the main pharmacodynamic endpoints. RESULTS: Mean clearance was 1.14 l h(-1) (CV 25%). Creatinine clearance was the only covariable to significantly reduce clearance variability (residual CV 18%). (CL = 0.74 + 0.0057 CLCR; r(2) = 0.32). Mean bioavailability was 45% (CV 22%) and mean protein binding 91.5% (CV 5%). Exposure to free, pharmacologically active etoposide (free AUC p.o.) was highly variable (mean value 2.8 mg l(-1) h; CV 64%; range 0.4-9.5). It decreased with increased creatinine clearance and increased with age which accounted for 9% of the CV. Mean free AUC p.o. was the best predictor of neutropenia. Free AUC50 (exposure producing a 50% reduction in absolute neutrophil count) was 1.80 mg l(-1) h. In patients with lung cancer, the free AUC p.o. was higher in the two patients with responsive tumour (5.9 mg l(-1) h) than in patients with stable (2.1 mg l-1 h) or progressive disease (2.3 mg l-1 h) (P = 0.01). CONCLUSIONS: Exposure to free etoposide during prolonged oral treatment is highly variable and is the main determinant of pharmacodynamic effects. The population PK model based on creatinine clearance is poorly predictive of exposure. Therapeutic drug monitoring would be necessary for dose individualization or to study the relationship between exposure and antitumour effect.     

Population pharmacokinetics and pharmacodynamics of meropenem in pediatric patients

Meropenem is a highly potent carbapenem antibiotic against gram-positive and gram-negative bacteria. Meropenem plasma concentration data from 99 pediatric patients (aged 0.08-17.3 years) were used to develop a population pharmacokinetic model. Pharmacokinetic analysis was performed using NONMEM with exponential interindividual variability and combinational residual error model. A 2-compartment model was found to fit the data best. Creatinine clearance and body weight were the most significant covariates explaining variabilities in meropenem pharmacokinetics among pediatric patients. The validated final model was used to predict meropenem plasma concentrations in 37 pediatric meningitis patients, receiving 40 mg/kg meropenem, who had minimum inhibitory concentration values of the causative pathogens and outcome available. Since the causative pathogens in all patients were eradicated, no break points for required exposure could be found. The microbiological outcomes indicate that the current clinical dosage regimen provides sufficient drug exposure to eradicate the pathogens commonly involved in pediatric meningitis.     

Population pharmacokinetics and pharmacodynamics of oral levodopa in parkinsonian patients

Objective: The population pharmacokinetics and pharmacodynamics of a standardised oral test dose of levodopa have been determined in patients with mild to severe Parkinson’s disease using parametric, non-linear mixed effect modelling with the program NONMEM. Levodopa plasma concentration data and motor effect behaviour (tapping times) were obtained from 46 patients, for whom a total of 970 observations were available (approximately 21 pharmacokinetic and pharmacodynamic observations per patient). The pharmacokinetic-pharmacodynamic model used was a one-compartment first-order absorption model linked to the sigmoid E_Max representation of the Hill equation via an equilibration rate-constant, ke₀. The model was also tested via a reduction in the number of pharmacokinetic and pharmacodynamic data points to a total of four to eight per patient. Results: In the final regression models the Hoehn and Yahr (HY) status of the patient and duration of disease (DUR) were found to be important determinants of the pharmacodynamic parameters for levodopa. The pharmacokinetic parameters were not significantly affected by any covariates. A test group of 16 additional parkinsonian patients was used to evaluate the predictive performance of the population parameters. The predictive performance of the pharmacokinetic-pharmacodynamic modelling using the full and reduced data sets was evaluated in NONMEM using posthoc, Bayesian forecasting. Statistically insignificant bias existed among predicted and observed levodopa concentrations, whereas the pharmacodynamic model underpredicted the observed tapping times. There was little difference in the pharmacokinetic-pharmacodynamic predictive performance among results for the full and the reduced data sets. Conclusion: In a clinical setting knowledge of the population pharmacokinetic and pharmacodynamic parameters for oral levodopa may prove useful in estimating the duration of the drug’s beneficial motor activity in patients with mild to severe Parkinson’s disease (Hoehn and Yahr status I–IV). Received: 13 November 1995/Accepted in revised form: 13 February 1996     

Population pharmacokinetics and pharmacodynamics of warfarin in healthy young adults

The population pharmacokinetics and pharmacological response — prothrombin complex activity and factor VII activity — were studied in a group of 48 normal, healthy young volunteers. Population parameter estimates were obtained using a standard two-stage method, a nonlinear mixed effect model (NONMEM) and a two-stage Bayesian method (EM algorithm). A modified sigmoid-Imax model was used to relate the concentration of s-warfarin to the rate of clotting factor synthesis. The three methods produced similar estimates of the population pharmacokinetic parameters, although the standard two-stage method overestimated the contribution of the pharmacokinetic parameters to the interindividual variability. It was not possible to partition the interindividual variability in response between the pharmacodynamic parameters with the NONMEM procedure: the estimates obtained from the EM algorithm were generally in good agreement with those obtained using the standard two-stage approach. The variability in the warfarin concentration contributed at most only 40% of the observed variability in the pharmacological response, and then only for times greater than 96 h after the dose. Most of the variability in the pharmacodynamics was due to interindividual differences in the clotting factor degradation rate constant and C_50,s, the s-warfarin concentration causing a 50% decrease in synthesis rate.     

Population pharmacokinetics and pharmacodynamics for treatment optimization in clinical oncology

Population pharmacokinetic and pharmacodynamic analysis is an important tool to support optimal treatment in clinical oncology. The population approach is suitable to explain variability between patients and to establish relationships between drug exposure and a relevant pharmacodynamic parameter. This can facilitate the selection of dosing schedules, the development of strategies for dose individualization and the application of therapeutic drug monitoring of anticancer agents. This review discusses the role of population pharmacokinetics and pharmacodynamics in clinical oncology to enhance the efficiency of drug development and to support the development of safe and effective dosing regimens for optimal treatment of cancer patients. An overview of published population studies of investigational anticancer agents and established treatment regimens is presented.     

Population pharmacokinetics/pharmacodynamics modeling: parametric and nonparametric methods

As clinicians acquire experience with the clinical and pharmacokinetic behavior of a drug, it is usually optimal to record this experience in the form of a population pharmacokinetic model, and then to relate the behavior of the model to the clinical effects of the drug or to a linked pharmacodynamic model. The role of population modeling is thus to describe and record clinical experience with the behavior of a drug in a certain group or population of patients or subjects.     

Population pharmacokinetics and pharmacodynamics of zidovudine in HIV-infected infants and children

The purpose of this study was to assess the population pharmacokinetics (PK) and pharmacodynamics (PD) of zidovudine (ZDV) in infants and children. This evaluation includes 394 subjects who participated in Pediatric AIDS Clinical Trials Group (PACTG) Study 152 and received either ZDV alone or in combination with didanosine. The most significant PK covariate was age, with infants < 2 years of age having reduced size-adjusted clearance. ZDV exposure was weakly related to maximal reduction in immune complex-dissociated (ICD) p24 antigen but not to reduction at 6 months. Mild chronic anemia occurred in 7.6% of subjects with ZDV average concentration < 1.3 microM (350 ng/mL) versus in 23.4% subjects with higher ZDV concentrations (p < 0.001). There was a direct linear relationship between hemoglobin and ZDV levels. It was concluded that ZDV oral clearance is reduced in infants compared to older children. This lower clearance leads to higher ZDV exposure in infants and contributes to increased hematologic toxicity.     

Population pharmacokinetics and pharmacodynamics of mitomycin during intraoperative hyperthermic intraperitoneal chemotherapy

BACKGROUND: During recent years, cytoreductive surgery combined with hyperthermic intraperitoneal chemotherapy (HIPEC) with mitomycin has been used for various malignancies. OBJECTIVE: To characterise the population pharmacokinetics and pharmacodynamics of mitomycin during HIPEC. METHODS: Forty-seven patients received mitomycin 35 mg/m2 intraperitoneally as a perfusion over 90 minutes. Mitomycin concentrations were determined in both the peritoneal perfusate and plasma. The observed concentration-time profiles were used to develop a population pharmacokinetic model using nonlinear mixed-effect modelling (NONMEM). The area under the plasma concentration-time curve (AUC) was related to the haematological toxicity. RESULTS: Concentration-time profiles of mitomycin in perfusate and plasma were adequately described with one- and two-compartment models, respectively. The average volume of distribution of the perfusate compartment (V1) and rate constant from the perfusate to the systemic circulation (k12) were 4.5 +/- 1.1L and 0.014 +/- 0.003 min(-1), respectively (mean +/- SD, n = 47). The average volume of distribution of the central plasma compartment (V2), clearance from the central compartment (CL) and volume of distribution of the peripheral plasma compartment (V3) were 28 +/- 16L, 0.55 +/- 0.18 L/min and 36 +/- 8L, respectively. The relationship between the AUC in plasma and degree of leucopenia was described with a sigmoidal maximum-effect (Emax) model. CONCLUSIONS: The pharmacokinetics of mitomycin during HIPEC could be fitted successfully to a multicompartment model. Relationships between plasma exposure and haematological toxicity were quantified. The developed pharmacokinetic-pharmacodynamic model can be used to simulate different dosage schemes in order to optimise mitomycin administration during HIPEC.     

Population pharmacokinetics and pharmacodynamics of rivaroxaban in patients with acute coronary syndromes

AIMS

The aim of this analysis was to use a population approach to facilitate the understanding of the pharmacokinetics and pharmacodynamics of rivaroxaban in patients with acute coronary syndrome (ACS) and to evaluate the influence of patient covariates on the exposure of rivaroxaban in patients with ACS.

METHODS

A population pharmacokinetic model was developed using pharmacokinetic samples from 2290 patients in Anti-Xa Therapy to Lower Cardiovascular Events in Addition to Standard Therapy in Subjects with Acute Coronary Syndrome Thrombolysis in Myocardial Infarction 46. The relationship between pharmacokinetics and the primary pharmacodynamic end point, prothrombin time, was evaluated.

RESULTS

The pharmacokinetics of rivaroxaban in patients with ACS was adequately described by an oral one-compartment model. The estimated absorption rate, apparent clearance and volume of distribution were 1.24 h⁻¹ (interindividual variability, 139%), 6.48 l h⁻¹ (31%) and 57.9 l (10%), respectively. Simulations indicate that the influences of renal function, age and bodyweight on exposure in ACS patients are consistent with the findings in previous Phase 1 studies. Rivaroxaban plasma concentrations exhibit a close-to-linear relationship with prothrombin time in the ACS population, with little interindividual variability. The estimated pharmacokinetic and pharmacodynamic parameters for the ACS patients were comparable to those for venous thromboembolism prevention, deep vein thrombosis and atrial fibrillation patients.

CONCLUSIONS

The similarity in pharmacokinetics/pharmacodynamics of rivaroxaban among different patient populations and the low interindividual variability in the exposure–prothrombin time relationship indicate that the anticoagulant effect of rivaroxaban is highly predictable and consistent across all the patient populations studied.     

Antibody pharmacokinetics and pharmacodynamics

The U.S. Food and Drug administration (FDA) has approved several polyclonal antibody preparations and at least 18 monoclonal antibody preparations (antibodies, antibody fragments, antibody fusion proteins, etc.). These drugs, which may be considered as a diverse group of therapeutic proteins, are associated with several interesting pharmacokinetic characteristics. Saturable binding with target antigen may influence antibody disposition, potentially leading to nonlinear distribution and elimination. Independent of antigen interaction, concentration-dependent elimination may be expected for IgG antibodies, due to the influence of the Brambell receptor, FcRn, which protects IgG from catabolism. Antibody administration may induce the development of an endogenous antibody response, which may alter the pharmacokinetics of the therapeutic antibody. Additionally, the pharmacodynamics of antibodies are also complex; these drugs may be used for a wide array of therapeutic applications, and effects may be achieved by a variety of mechanisms. This article provides an overview of many of the complexities associated with antibody pharmacokinetics and pharmacodynamics.     

Population pharmacokinetics and pharmacodynamics of bivalirudin in young healthy Chinese volunteers

Aim:

To investigate the population pharmacokinetics (PK) and pharmacodynamics (PD) of bivalirudin, a synthetic bivalent direct thrombin inhibitor, in young healthy Chinese subjects.

Methods:

Thirty-six young healthy volunteers were randomly assigned into 4 groups received bivalirudin 0.5 mg/kg, 0.75 mg/kg, and 1.05 mg/kg intravenous bolus, 0.75 mg/kg intravenous bolus followed by 1.75 mg/kg intravenous infusion per hour for 4 h. Blood samples were collected to measure bivalirudin plasma concentration and activated clotting time (ACT). Population PK-PD analysis was performed using the nonlinear mixed-effects model software NONMEM. The final models were validated with bootstrap and prediction-corrected visual predictive check (pcVPC) approaches.

Results:

The final PK model was a two-compartment model without covariates. The typical PK population values of clearance (CL), apparent distribution volume of the central-compartment (V₁), inter-compartmental clearance (Q) and apparent distribution volume of the peripheral compartment (V₂) were 0.323 L·h^-1·kg^-1, 0.086 L/kg, 0.0957 L·h^-1·kg^-1, and 0.0554 L/kg, respectively. The inter-individual variabilities of these parameters were 14.8%, 24.2%, fixed to 0% and 15.6%, respectively. The final PK-PD model was a sigmoid E_max model without the Hill coefficient. In this model, a covariate, red blood cell count (RBC^*), had a significant effect on the EC₅₀ value. The typical PD population values of maximum effect (E_max), EC₅₀, baseline ACT value (E₀) and the coefficient of RBC^* on EC₅₀ were 318 s, 2.44 mg/L, 134 s and 1.70, respectively. The inter-individual variabilities of E_max, EC₅₀, and E₀ were 6.80%, 46.4%, and 4.10%, respectively.

Conclusion:

Population PK-PD models of bivalirudin in healthy young Chinese subjects have been developed, which may provide a reference for future use of bivalirudin in China.     

Clinical pharmacokinetics and pharmacodynamics and pharmacodynamics of artemether-lumefantrine.

The combination of artemether and lumefantrine (benflumetol) is a new and very well tolerated oral antimalarial drug effective even against multidrug-resistant falciparum malaria. The artemether component is absorbed rapidly and biotransformed to dihydroartemisinin, and both are eliminated with terminal half-lives of around 1 hour. These are very active antimalarials which give a rapid reduction in parasite biomass and consequent rapid resolution of symptoms. The lumefantrine component is absorbed variably in malaria, and is eliminated more slowly (half-life of 3 to 6 days). Absorption is very dependent on coadministration with fat, and so improves markedly with recovery from malaria. Thus artemether clears most of the infection, and the lumefantrine concentrations that remain at the end of the 3- to 5-day treatment course are responsible for eliminating the residual 100 to 10 000 parasites. The area under the curve of plasma lumefantrine concentrations versus time, or its correlate the plasma concentration on day 7. has proved an important determinant of therapeutic response. Characterisation of these pharmacokinetic-pharmacodynamic relationships provided the basis for dosage optimisation, an approach that could be applied to other antimalarial drugs.     

Population pharmacokinetics and pharmacodynamics of escitalopram in overdose and the effect of activated charcoal

AIMS

To describe the pharmacokinetics and pharmacodynamics (PKPD) of escitalopram in overdose and its effect on QT prolongation, including the effectiveness of single dose activated charcoal (SDAC).

METHODS

The data set included 78 escitalopram overdose events (median dose, 140 mg [10–560 mg]). SDAC was administered 1.0 to 2.6 h after 12 overdoses (15%). A fully Bayesian analysis was undertaken in WinBUGS 1.4.3, first for a population pharmacokinetic (PK) analysis followed by a PKPD analysis. The developed PKPD model was used to predict the probability of having an abnormal QT as a surrogate for torsade de pointes.

RESULTS

A one compartment model with first order input and first-order elimination described the PK data, including uncertainty in dose and a baseline concentration for patients taking escitalopram therapeutically. SDAC reduced the fraction absorbed by 31% and reduced the individual predicted area under the curve adjusted for dose (AUC_i/dose). The absolute QT interval was related to the observed heart rate with an estimated individual heart rate correction factor (α = 0.35). The heart rate corrected QT interval (QT_c) was linearly dependent on predicted escitalopram concentration [slope = 87 ms/(mg l^–1)], using a hypothetical effect-compartment (half-life of effect-delay, 1.0h). Administration of SDAC significantly reduced QT prolongation and was shown to reduce the risk of having an abnormal QT by approximately 35% for escitalopram doses above 200 mg.

CONCLUSIONS

There was a dose-related lengthening of the QT interval that lagged the increase in drug concentration. SDAC resulted in a moderate reduction in fraction of escitalopram absorbed and reduced the risk of the QT interval being abnormal.     

Population pharmacokinetics and pharmacodynamics of enoxaparin in unstable angina and non-ST-segment elevation myocardial infarction

AIMS: A major concern with any antithrombotic therapy is an increase in the risk of haemorrhage. The aim of this study was to analyse population pharmacokinetics and pharmacokinetic/pharmacodynamic (PK/PD) relationships for enoxaparin in patients with unstable angina (UA) and non-ST-segment elevation myocardial infarction (NSTEMI), which may help predict risk of haemorrhage. METHODS: Anti-factor Xa (anti-Xa) activity was measured as marker of enoxaparin concentration in 448 patients receiving the drug as a single 30-mg intravenous bolus followed by 1.0 or 1.25 mg kg(-1) subcutaneously twice a day. A population pharmacokinetic analysis was conducted and individual estimates of enoxaparin clearance and area under the curve were tested as prognostic factors for the occurrence of haemorrhagic episodes. RESULTS: Basic population PK parameters were an enoxaparin clearance of 0.733 l h(-1)[95% confidence interval (CI) 0.698, 0.738], a distribution volume of 5.24 l (95% CI 4.20, 6.28) and an elimination half-life of 5.0 h. Enoxaparin clearance was significantly related to patient weight and creatinine clearance, and was the only independent predictor of experiencing both all (10.7%, P = 0.0013) and major (2.2%, P = 0.0004) haemorrhagic events. A creatinine clearance of 30 ml min(-1) was associated with a decrease in enoxaparin clearance of 27% compared with that in a patient with a median creatinine clearance of 88 ml min-1, and was related to a 1.5- and 3.8-fold increase in the risk of 'all' and 'major' haemorrhagic episodes, respectively. CONCLUSIONS: Enoxaparin clearance depends on body weight, and, therefore, weight-adjusted dosing is recommended to minimize interpatient variability in drug exposure and the risk of haemorrhage. The importance of an increased risk of haemorrhage with decreasing renal function must be weighed against the benefit of treatment with enoxaparin in patients with UA and NSTEMI.     

Macrolides: pharmacokinetics and pharmacodynamics

Three pharmacokinetic/pharmacodynamic parameters--(i) the peak concentration to the minimum inhibitory concentration ratio (C(max)/MIC); (ii) the area under the concentration-time curve to MIC ratio (AUC(24h)/MIC); and (iii) the time the concentration exceeds the MIC (T>MIC)--are important predictors of the clinical efficacy of antibiotics. For antibiotics with pronounced concentration-dependent killing, such as the fluoroquinolones or the aminoglycosides, C(max)/MIC and AUC(24)/MIC are the main factors that establish efficacy. Antibiotics with a weak, or no, concentration dependency generally have their efficacy linked to T>MIC, and these include the beta-lactams and the conventional macrolides. Antibiotics with weak concentration-dependent effects, but with prolonged persistent effects, such as tetracyclines and azithromycin, have their activity mostly related to the AUC(24)/MIC. By applying these concepts to current antibiotics, and also to the development of novel agents, it is possible to optimise their dosages and administration schedules. This will maximise therapeutic efficacy, may prevent or delay the emergence of bacterial resistance to antibiotics, and can certainly minimise side-effects.     

Heparin pharmacokinetics and pharmacodynamics.

Heparin was discovered approximately 75 years ago and has been used extensively for the last 50 years to treat thromboembolic disorders. An endogenous glycosaminoglycan, heparin is found largely in the liver, lung and intestine. It is available for exogenous administration both as unfractionated and low molecular weight heparin. Unfractionated heparin is a heterogenous mixture of polysaccharide chains of varying length resulting in a range of molecular weights from 3000 to 30,000D while low molecular weight heparin ranges from 3000 to 6000D. Heparin produces its antithrombotic effect by binding to antithrombin III and this complex then binds to thrombin. In order to accomplish this a total of 18 to 22 monosaccharide units is necessary including a specific pentasaccharide binding site for antithrombin III. After either subcutaneous or intravenous injection heparin is distributed primarily within the intravascular space. A short distribution phase is seen which is thought to correspond to endothelial cell binding and internalisation. The disposition curve for unfractionated heparin has a unique concave-convex shape which is the result of combined saturable and nonsaturable elimination mechanisms. The nonsaturable elimination mechanism is renal and is the primary route of elimination for low molecular weight heparins. For this reason, the concave-convex pattern is not seen with low molecular weight preparations. Both forms of heparin are useful antithrombotic agents; however, the correlation between the antithrombotic effect and an in vitro laboratory test for either type still needs further clarification.     

Population pharmacokinetics and pharmacodynamics of dalfampridine-ER in healthy volunteers and in patients with multiple sclerosis

Abstract

Objective:

Using data pooled from several studies of dalfampridine extended release (ER), a population pharmacokinetic model was developed for the purposes of characterizing the population pharmacokinetics and pharmacodynamics of dalfampridine-ER with respect to variability in pharmacokinetics, covariates affecting the pharmacokinetics, and whether the current therapeutic dosage represents an optimal dosage. Studies were conducted in healthy volunteers and multiple sclerosis (MS) patients over the course of development and registration of dalfampridine extended release tablets (dalfampridine-ER [Ampyra]; prolonged-, modified- or sustained-release fampridine [Fampyra] in some countries).