Prospective use of population pharmacokinetics/pharmacodynamics in the development of cisatracurium

Location matters: Overlooked ethnic-geographic effect in China and Austria on propofol/cisatracurium sex differences among a population pharmacokinetic/pharmacodynamic (PopPK/PD) covariate analysis in men, women, and one transgender subject

A quick literature search using "sex/gender" vs. the commonly used hypnotic propofol or neuromuscular-blocking agent cisatracurium will reveal numerous contradictory sex difference publications depending on the ethnic-geographic location of where these studies were conducted. We induced anesthesia with cisatracurium besylate (GlaxoSmithKline) 100 μg kg^-1 administered exactly 1 minute following propofol (AstraZeneca) 2 mg kg^-1 . In 20 male and 20 female ethnic Han-Chinese test set patients (Xi'an China), and in similar ethnic white Austrian validation set patients (Graz Austria), we quantified propofol/cisatracurium pharmacodynamic parameters namely propofol onset time, lag time, plasma concentrations (C_p ) at loss-of-behavioral response (LOBR) using bispectral index (BIS); cisatracurium onset time, lag time, and C_p at T₁ % (first twitch of train-of-four) complete twitch suppression using mechanomyography (MMG). Serial arterial blood samples were collected for population pharmacokinetic (PopPK) analysis of all demographic and biological covariates (region, sex, age, weight, and height) versus volumes of distribution and clearances pharmacokinetic parameters. In Chinese women (but not in white women), propofol C_p at LOBR was 33.60% lower than men and cisatracurium C_p at T₁ % complete twitch suppression was 21.49% lower than men, a clear pharmacodynamic assertion. Region and weight were significant PopPK covariates. We demonstrated that sex differences are influenced by ethnic-geographic location as only in Chinese women (but not in white women) propofol C_p at LOBR and cisatracurium C_p at T₁ % complete twitch suppression were lower than in men. When defining sex differences, ethnic-geographic location should be taken into consideration as a predictive factor for optimizing propofol/cisatracurium initial loading recommended dosages.

Pharmacokinetic/Pharmacodynamic Model of CW002, an Investigational Intermediate Neuromuscular Blocking Agent, in Healthy Volunteers

BACKGROUND

CW002 is an investigational nondepolarizing, neuromuscular blocking agent with a rapid onset and intermediate duration of action in animals. This is a single ascending dose, healthy subject study exploring tolerability, pharmacokinetics, and potency.

METHODS

Population pharmacokinetic and pharmacokinetic/pharmacodynamic models were developed using plasma drug concentration data from a previously published dose-response study in 28 healthy subjects receiving single doses of CW002 during sevoflurane anesthesia. Subjects included in the models were from five different dose cohorts (cohorts 3, 4, 5, 6, and 8 receiving 0.04, 0.06, 0.08, 0.10, and 0.14 mg/kg, respectively). Serial arterial plasma concentrations and muscle twitch heights were monitored.

RESULTS

A four-compartment model was fit to the concentration-time data, whereas a transit compartment with a sigmoid Emax model was fit to the pharmacokinetic/pharmacodynamic data. The population pharmacokinetics of CW002 was linear with very low interindividual variability in clearance (10.8%). Simulations were conducted to predict the onset and offset of effect at 2×, 3×, and 4× ED95. The time to 80% block was predicted to be 1.5, 0.8, and 0.7 min for 2×, 3×, and 4× ED95 doses, respectively. The simulated 25 to 75% recovery index was independent of dose.

CONCLUSIONS

CW002 has predictable pharmacokinetics and is likely to have a rapid onset with an intermediate duration of action at 3× ED95. This model provides information to inform critical decisions (e.g., dose, study design) for continued development of CW002.

Role of SLCO1B1, ABCB1, and CHRNA1 gene polymorphisms on the efficacy of rocuronium in Chinese patients

This study explored the role of SLCO1B1, ABCB1, and CHRNA1 gene polymorphisms on the efficacy and duration of action of rocuronium in Chinese patients. Two hundred seven unrelated Chinese patients scheduled for elective surgery were recruited, and 200 completed the study. Their ABCB1, SLCO1B1, and CHRNA1 genotypes were determined. Demographic and clinical non-genetic data also were collected. The SLCO1B1, ABCB1, and CHRNA1 variants did not affect the onset time of rocuronium. Clinical duration and recovery time of rocuronium were prolonged in patients with the ABCB1 rs1128503TT and SLCO1B1 rs2306283 AG and GG genotypes. We demonstrate that the SLCO1B1 and ABCB1 gene variants could affect the pharmacodynamics of rocuronium. The ABCB1 rs1128503 C>T genotype was the most important factor on the efficacy of rocuronium.

Influence of gender on the course of neuromuscular block following a single bolus dose of cisatracurium or rocuronium

BACKGROUND AND OBJECTIVE

There is increasing evidence of gender differences in the pharmacokinetics and pharmacodynamics of aminosteroid neuromuscular blocking agents. Compared to males, females are more susceptible, requiring approximately 30% less rocuronium to achieve the same degree of neuromuscular block. However, little information is available whether this difference is applicable to modern benzylisoquinolines (cisatracurium).

METHODS

In all, 848 patients (423 males, 425 females) undergoing general surgery under total intravenous anaesthesia with muscle relaxation, tracheal intubation and mechanical ventilation were studied. Patients were randomized to receive a single bolus dose of cisatracurium (0.1 mg kg-1, 221 males and 199 females) or rocuronium (0.6 mg kg-1, 202 males and 226 females). The onset time for 95% depression of T1, clinical duration until 25% recovery and recovery index (T1 from 25% to 75%) were determined with an NMT electromyographic module of the Datex-Ohmeda S/5 Anaesthesia Monitor. The data for male and female groups were compared with appropriate statistical tests (unpaired t-test, Mann-Whitney rank sum test and Fisher's exact test).

RESULTS

In both groups (cisatracurium and rocuronium), males were significantly taller (P < 0.001) and heavier (P < 0.001) than females, but the body mass index was comparable. For rocuronium, the onset time was shorter 91.7 +/- 14.3 s vs. 108.0 +/- 14.6 s (P < 0.001) and the clinical duration was increased in females 43.3 +/- 7.8 min vs. 31.3 +/- 5.5 min (P < 0.001). In the cisatracurium group, both onset times (248.9 +/- 60.7 s for males vs. 253.4 +/- 70.9 s for females) and clinical duration (42.6 +/- 6.9 min for males vs. 43.1 +/- 6.9 min for females) were similar. The recovery index was identical for males and females in both groups.

CONCLUSIONS

Females were more sensitive than males to a single bolus dose of rocuronium. Under the study conditions described, the onset time was shorter and the clinical duration was increased in female patients. This suggests that the routine dose of rocuronium should be reduced in females compared to males. On the contrary, we could demonstrate no gender differences in the onset time or clinical duration of cisatracurium.

Population pharmacokinetic-pharmacodynamic modeling of TF-505 using extension of indirect response model by incorporating a circadian rhythm in healthy volunteers

The pharmacokinetic-pharmacodynamic (PK-PD) relationship of the newly developed drug, (-)-(S)-4-[1-[4-[1-(4-isobutylphenyl)butoxy]benzoyl]indolizin-3-yl]butyric acid (TF-505), was characterized via a population approach in early human study. Healthy volunteers were divided into six groups. The groups received four single doses (25, 50, 75 or 100 mg) and 2 multiple doses (12.5 or 25 mg) of TF-505, respectively. Dihydrotestosterone (DHT) data were collected to assess TF-505 pharmacodynamics. Population PK/PD modeling of TF-505 was performed via mixed-effects modeling using the NONMEM software package. The final PK-PD model incorporates a two-compartment PK model and an extended indirect PD model. The population PK parameters were 0.197 h(-1) for the k(a), 0.0678 h(-1) for k(e), 12.5 l for V(c), 0.0645 h(-1) for k(12), 0.0723 h(-1) for k(21). Extension of indirect response model by incorporating a time-dependent periodic function for k(in) takes into account the chronopharmacologic rhythms (I(max): 0.706+/-0.297, IC(50): 1.01+/-1.64 (microg/ml), k(out): 0.221+/-0.0486 (h(-1)), R(m): 20.4+/-8.08 (% h(-1)), R(amp): 5.06+/-3.43 (% h(-1)), T(z): 5.01+/-0.407 (h) (Population mean+/-S.E.)). R(m) is the mean DHT synthesis rate, R(amp) is the amplitude of the DHT synthesis rate, and T(z) is the acrophase time, signifying maximum synthesis rate. The present study represents a successful population PK-PD model using the full data from early human studies. The population parameters thus obtained could provide useful indicators for the determination of dosage regimens in exploratory studies in patient populations.

Challenges in the transition to model-based development

Practitioners of the art and science of pharmacometrics are well aware of the considerable effort required to successfully complete modeling and simulation activities for drug development programs. This is particularly true because of the current, ad hoc implementation wherein modeling and simulation activities are piggybacked onto traditional development programs. This effort, coupled with the failure to explicitly design development programs around modeling and simulation, will continue to be an important obstacle to the successful transition to model-based drug development. Challenges with timely data availability, high data discard rates, delays in completing modeling and simulation activities, and resistance of development teams to the use of modeling and simulation in decision making are all symptoms of an immature process capability for performing modeling and simulation. A process that will fulfill the promise of model-based development will require the development and deployment of three critical elements. The first is the infrastructure--the data definitions and assembly processes that will allow efficient pooling of data across trials and development programs. The second is the process itself--developing guidelines for deciding when and where modeling and simulation should be applied and the criteria for assessing performance and impact. The third element concerns the organization and culture--the establishment of truly integrated, multidisciplinary, and multiorganizational development teams trained in the use of modeling and simulation in decision-making. Creating these capabilities, infrastructure, and incentivizations are critical to realizing the full value of modeling and simulation in drug development.

Cisatracurium Pharmacodynamics in Patients with Oculopharyngeal Muscular Dystrophy

The pharmacodynamics of muscle relaxants in patients with oculopharyngeal muscular dystrophy (OPMD) have never been studied. We designed this study to compare the pharmacodynamics of cisatracurium in OPMD patients versus a control group. Forty patients were enrolled: 20 OPMD patients requiring general anesthesia for cricopharyngeal myotomy and 20 age-matched controls undergoing an operation of similar duration and expected blood loss. Anesthesia was standardized, and both groups received a bolus of cisatracurium 0.1 mg/kg. Onset time, time to 10% T1 recovery, and the intervals 10%-25% and 25%-75% were calculated for both groups. A subgroup analysis was performed in patients with a more severe form of OPMD. Demographic and intraoperative data were similar. Onset time was significantly longer in OPMD patients compared with the control group (4.6 +/- 1.5 min versus 3.4 +/- 1.0 min; P = 0.001). There was no difference in recovery times or indices between groups, regardless of the severity of the disease. In conclusion, there was no difference in the duration of a cisatracurium-induced neuromuscular block between OPMD patients and a control group. A delayed onset of action of the drug may occur.

Gender differences in drug effects: implications for anesthesiologists

BACKGROUND

The gender aspect in pharmacokinetics and pharmacodynamics of anesthetics has attracted little attention. Knowledge of previous work is required to decide if gender-based differences in clinical practice is justified, and to determine the need for research.

METHODS

Basis for this paper was obtained by Medline searches using the key words 'human' and 'gender' or 'sex,' combined with individual drug names. The reference lists of these papers were further checked for other relevant studies.

RESULTS

Females have 20-30% greater sensitivity to the muscle relaxant effects of vecuronium, pancuronium and rocuronium. When rapid onset of or short duration of action is very important, gender-modified dosing may be considered. Males are more sensitive than females to propofol. It may therefore be necessary to decrease the propofol dose by 30-40% in males compared with females in order to achieve similar recovery times. Females are more sensitive than males to opioid receptor agonists, as shown for morphine as well as for a number of kappa (OP2) receptor agonists. On this basis, males will be expected to require 30-40% higher doses of opioid analgesics than females to achieve similar pain relief. On the other hand, females may experience respiratory depression and other adverse effects more easily if they are given the same doses as males.

CONCLUSION

These examples illustrate that gender should be taken into account as a factor that may be predictive for the dosage of several anesthetic drugs. Moreover, there is an obvious need for more research in this area in order to further optimize drug treatment in anesthesia.

Pharmacokinetic/pharmacodynamic studies in drug product development

In the quest of ways for rationalizing and accelerating drug product development, integrated pharmacokinetic/pharmacodynamic (PK/PD) concepts provide a highly promising tool. PK/PD modeling concepts can be applied in all stages of preclinical and clinical drug development, and their benefits are multifold. At the preclinical stage, potential applications might comprise the evaluation of in vivo potency and intrinsic activity, the identification of bio-/surrogate markers, as well as dosage form and regimen selection and optimization. At the clinical stage, analytical PK/PD applications include characterization of the dose-concentration-effect/toxicity relationship, evaluation of food, age and gender effects, drug/drug and drug/disease interactions, tolerance development, and inter- and intraindividual variability in response. Predictive PK/PD applications can also involve extrapolation from preclinical data, simulation of drug responses, as well as clinical trial forecasting. Rigorous implementation of the PK/PD concepts in drug product development provides a rationale, scientifically based framework for efficient decision making regarding the selection of potential drug candidates, for maximum information gain from the performed experiments and studies, and for conducting fewer, more focused clinical trials with improved efficiency and cost effectiveness. Thus, PK/PD concepts are believed to play a pivotal role in streamlining the drug development process of the future.

Pharmacokinetic studies of neuromuscular blocking agents: good clinical research practice (GCRP)

In September 1997, an international consensus conference on standardization of studies of neuromuscular blocking agents was held in Copenhagen, Denmark. Based on the conference, a set of guidelines for good clinical research practice (GCRP) in pharmacokinetic studies of neuromuscular blocking agents is presented. Guidelines include: design of the study; relevant patient groups to investigate; test drug administration, sampling and analysis; pharmacokinetic analysis; pharmacokinetic/pharmacodynamic modeling; population pharmacokinetics; statistics; and presentation of pharmacokinetic data. The guidelines are intended to aid those working in this research area; it is hoped that they will assist researchers, editors of scientific papers, and pharmaceutical companies in improving the quality of pharmacokinetic studies.

Assuming peripheral elimination: its impact on the estimation of pharmacokinetic parameters of muscle relaxants

For anesthetic drugs undergoing nonorgan-based elimination, there is a definite trend towards using pharmacokinetic (PK) models in which elimination can occur from both central (k10) and peripheral compartments (k20). As the latter cannot be assessed directly, assumptions have to be made regarding its value. The primary purpose of this paper is to evaluate the impact of assuming various degrees of peripheral elimination on the estimation of PK parameters. For doing so, an explanatory model is presented where previously published data from our laboratory on three muscle relaxants, i.e., atracurium, doxacurium, and mivacurium, are used for simulations. The mathematical aspects for this explanatory model as well as for two specific applications are detailed. Our simulations show that muscle relaxants having a short elimination half-life are more affected by the presence of peripheral elimination as their distribution phase occupies the major proportion of their total area under the curve. Changes in the exit site dependent PK parameters (Vdss) are also mostly significant when k20 is smaller than k10. Although the physiological processes that determine drug distribution and those affecting peripheral elimination are independent, the two are mathematically tied together in the two-compartment model with both central and peripheral elimination. It follows that, as greater importance is given to k20, the rate of transfer from the central compartment (k12) increases. However, as a result of a proportional increase in the volume of the peripheral compartment, peripheral concentrations remain unchanged whether or not peripheral elimination is assumed. These findings point out the limitations of compartmental analysis when peripheral elimination cannot be measured directly.

Pharmacokinetics and pharmacodynamics of cisatracurium after a short infusion in patients under propofol anesthesia

Fourteen patients, ASA physical status I or II, were recruited to assess the pharmacokinetic-pharmacodynamic relationship of cisatracurium under nitrous oxide/sufentanil/propofol anesthesia. The electromyographic response of the abductor digiti minimi muscle was recorded on train-of-four stimulation of the ulnar nerve. A 0.1-mg/kg dose of cisatracurium was given as an infusion over 5 min. Arterial plasma concentrations of cisatracurium and its major metabolites were measured by using high-performance liquid chromatography. A nontraditional two-compartment pharmacokinetic model with elimination from central and peripheral compartments was used. The elimination rate constant from the peripheral compartment was fixed to the in vitro rate of degradation of cisatracurium in human plasma (0.0237 min(-1)). The mean terminal half-life of cisatracurium was 23.9+/-3.3 min, and its total clearance averaged 3.7+/-0.8 mL x min(-1) x kg(-1). Using this model, the volume of distribution at steady state was significantly increased compared with that obtained when central elimination only was assumed (0.118+/-0.027 vs 0.089+/-0.017 L/kg). The effect-plasma equilibration rate constant was 0.054+/-0.013 min(-1). The 50% effective concentration (153+/-33 ng/mL) was 56% higher than that reported in patients anesthetized with volatile anesthetics, which suggests that, compared with inhaled anesthetics, a cisatracurium neuromuscular block is less enhanced by propofol.

IMPLICATIONS

The drug concentration-effect relationship of the muscle relaxant cisatracurium has been characterized under balanced and isoflurane anesthesia. Because propofol is now widely used as an IV anesthetic, it is important to characterize the biological fate and the concentration-effect relationship of cisatracurium under propofol anesthesia as well.

Population pharmacokinetics in veterinary medicine: potential use for therapeutic drug monitoring and prediction of tissue residues

Population pharmacokinetics can be defined as a study of the basic features of drug disposition in a population, accounting for the influence of diverse pathophysiological factors on pharmacokinetics, and explicitly estimating the magnitude of the interindividual and intraindividual variability. It is used to identify subpopulations of individuals that may present with differences in drug kinetics or in kinetic/dynamic responses. Rooted in procedures used in engineering systems, population pharmacokinetics methods were conceived as a means to determine the pharmacokinetic profile in populations in which a sparse number of samples were obtained per individual, such as those in late stage human clinical trials. This is the situation commonly encountered in all aspects of veterinary medicine. The exploratory nature of this technique allows one to probe relationships between clinical factors (such as age, gender, renal function, etc.) and drug disposition and/or effect. Similarly, the utilization of these techniques in the clinical research phases of drug development optimize the determination of efficacy and safety of drugs. Given the observational nature of most studies published so far, statistical methods to validate the population models are necessary. Simulation studies may be conducted to explore data collection designs that maximize information yield with a minimum expenditure of resources. The breadth of this approach has allowed population studies to be more commonly employed in all areas of drug therapy and clinical research. Finally, in veterinary medicine, there is an additional field in which population studies are potentially ideally suited: the application of this methodology to the study of tissue drug depletion and drug residues in production animals, and the establishment of withdrawal times tailored to the clinical or production conditions of populations or individuals. Such application would provide a major step toward assuring a safe food supply under a wide variety of dose and off-label clinical uses. Population pharmacokinetics is an ideal method for generating data in support of the implementation of flexible labelling policies and extralabel drug use recently approved under AMDUCA (Animal Medicinal Drug Use Clarification Act. 21 CFR Part 530).

Dose proportionality of cisatracurium

The dose proportionality of cisatracurium pharmacokinetics was assessed using a population approach by incorporating the collection of sparse blood samples from patients in clinical trials. Plasma concentration-time data from 131 patients with limited concentration-time data and 38 patients with full sampling were pooled and analyzed using nonlinear mixed-effects modeling (NONMEM). Dose proportionality was assessed using dichotomous parameterization and a linear model. The population pharmacokinetic approach revealed that the pharmacokinetics of cisatracurium are independent of dose between 0.1 mg/kg and 0.4 mg/kg, as was expected based on the importance of Hofmann elimination, a chemical process dependent on pH and temperature.