Pharmacokinetics and Pharmacodynamics of Cisatracurium After a Short Infusion in Patients Under Propofol Anesthesia

Synergism between rocuronium and cisatracurium: comparison of the Minto and Greco interaction models

Background

This study was conducted to investigate the pharmacodynamic interaction between rocuronium and cisatracurium using the response surface model, which is not subject to the limitations of traditional isobolographic analysis.

Methods

One hundred and twenty patients were randomly allocated to receive one of the fifteen predefined combinations of rocuronium and cisatracurium. To study single drugs, cisatracurium 0.2, 0.15, or 0.1 mg/kg or rocuronium 0.8, 0.6 or 0.4 mg/kg doses were administered alone. To study the pharmacodynamic interaction, drugs were applied in three types of combination ratio, i.e., half dose of each drug alone, 75% of each single dose of rocuronium and 25% of each single dose of cisatracurium, and vice versa. Train-of-four (TOF) ratio and T1% (first twitch of the TOF presented as percentage compared to the initial T1) were used as pharmacodynamic endpoints, and the Greco and Minto models were used as surface interaction models.

Results

The interaction term α of the Greco model for TOF ratio and T1% measurements showed synergism with values of 0.977 and 1.12, respectively. Application of the Minto model resulted in U₅₀ (θ) values (normalized unit of concentration that produces 50% of the maximal effect in the 0 < θ < 1 region) less than 1 for both TOF ratio and T1% measurements, indicating that rocuronium and cisatracurium exhibit synergism.

Conclusions

Response surface modeling of the interaction between rocuronium and cisatracurium, based on considerations of their effects on muscle relaxation as measured by TOF ratio and T1%, indicated that the two drugs show considerable synergism.

Simulation of the kinetics of neuromuscular block: implications for speed of onset

BACKGROUND

The onset time for paralysis varies 3-fold among nondepolarizing muscle relaxants. Possible explanations include: (a) pharmacokinetic differences among drugs and (b) buffering of drug molecules by acetylcholine receptors as they diffuse into the neuromuscular junction. Although some pharmacokinetic models consider buffered diffusion, these models do not account for either the high density of receptors or synapse geometry. Here, I used computer simulations to calculate the kinetics of buffered diffusion. The goal was to determine the conditions under which buffered diffusion could account for differences in onset time among nondepolarizing muscle relaxants.

METHODS

Monte Carlo simulation was used along with a realistic 3-dimensional model of the rat neuromuscular junction. Simulations determined the time dependence of the number of drug-bound receptors. A 1000-fold range of drug potency was examined. In some simulations, the drug concentration outside the junction was changed instantaneously. In other simulations, the concentration changed according to predictions of pharmacokinetic models assuming time-dependent changes in plasma drug concentration. The rate constant for equilibration of drug between plasma and muscle, keo, was varied between 0.15 and 0.6 min(-1). Twitch amplitude was calculated from receptor occupancy assuming a high safety margin for neuromuscular transmission. Some simulations used a synaptic model with an increased nerve-muscle contact width.

RESULTS

Simulations with instantaneous changes in drug concentration at the synapse, indicated that the time to 50% twitch depression (onset time) was 0.1 to 30 seconds and was proportional to drug potency. This corresponds to iontophoretic application of drug to isolated neuromuscular junctions, but is too fast to explain onset times in humans. When pharmacokinetic models were used to calculate the drug concentration outside the synapse, buffered diffusion increased onset times of potent drugs (drugs for which the effective concentration at 50% twitch height is <600 nM). Simulations using keo = 0.6 min(-1) and a model with a 2- to 3-fold wider nerve-muscle contact width indicated that buffered diffusion could account for the differences in clinical onset times among the nondepolarizing muscle relaxants.

CONCLUSION

Monte Carlo simulation provides a biophysically appropriate way to incorporate buffered diffusion into pharmacokinetic modeling. The simulations indicated that buffered diffusion could account for differences in onset time among drugs. However, a better understanding of the geometry of the human neuromuscular junction is needed before the magnitude of the effect of buffered diffusion can be quantified.

MUSCLE RELAXATION CONTROLLED BY AUTOMATED ADMINISTRATION OF CISATRACURIUM

In this paper, the muscle relaxant agent (i.e., cisatracurium) and three clinical control methods (i.e., 13 patients undergoing intermitted bolus control, 15 patients undergoing intensive manual control and 15 patients) undergoing automatic fuzzy logic control (FLC), were used for maintaining depth of muscle relaxation (DOM) during surgery. Cisatracurium, a muscle relaxation drug with long-term effect, low metabolic loading, but long delay time, is widely used in operating rooms and ICUs. Meanwhile, the rules for the FLC were developed from the experimental experience of intensive manual control after learning from 15 patient trials. According to experts' experimental experience, our FLC inputs were chosen from T1% error and trend of T1% which differ from other previous studies on eliminating the effect of time delay from cisatracurium. In individual clinical experimental results, the mean(SD) of the mean T1% error in 13 patients for intermitted bolus control, in 15 patients for intensive manual control, and in 15 patients for automatic control was 8.76(1.46), 1.65(1.67), and 0.48(1.43), respectively. The t test results show that automatic control is not significantly different from intensive manual control. The results show that a simple fuzzy logic controller derived from anesthetists ' clinical trials can provide good accuracy without being affected by the pharmacological time delay problem.

Effect of isoflurane versus propofol-remifentanil anesthesia on neuromuscular blockade and hemodynamic responses by cisatracurium bolus injection

BACKGROUND

Inhalational anesthetics potentiate nondepolarizing muscle relaxants. Cisatracurium is a recently introduced neuromuscular blocker in Korea. We studied the effect of inhalational anesthesia and total intravenous anesthesia (TIVA) on neuromuscular blockades and hemodynamic responses by cisatracurium bolus injection.

METHODS

Forty patients undergoing elective surgery were randomly divided into isoflurane and propofol-remifentanil groups. A bolus dose of cisatracurium of 0.15 mg/kg (3 × ED(95)) was administered after induction and the onset time and clinical duration of action were recorded. The nueromuscular blockade was monitored using train-of-four (TOF) stimulation. Hemodynamic parameters were also recorded.

RESULTS

Onset time was 194.0 ± 39.1 sec in the isoflurane group and 226.5 ± 62.2 sec in the propofol-remifentanil group. Clinical duration of action was 49.2 ± 9.0 min in the isoflurane group and 43.0 ± 9.2 min in the propofol-remifentanil group. Systolic blood pressure (SBP) and diastolic blood pressure (DBP) immediately before intubation decreased in the propofol-remifentanil group. Heart rate (HR), SBP and DBP 1 and 3 min after tracheal intubation increased in the isoflurane group.

CONCLUSIONS

Onset time was similar between isoflurane and propofol-remifentanil anesthesia. Clinical duration of action was significantly longer in isoflurane anesthesia. SBP and DBP immediately before intubation and HR, SBP and DBP 1 and 3 min after tracheal intubation were significantly different between the two groups.

Physicochemical properties of neuromuscular blocking agents and their impact on the pharmacokinetic-pharmacodynamic relationship

BACKGROUND

Among the factors influencing the onset of action of neuromuscular blocking agents (NMBA), the potency (EC50) and the rate of equilibration between blood and the effect compartment (k(e0)) have been highlighted. Although these descriptors are intrinsically influenced by the physicochemical characteristics of the drug, the impact of lipid solubility, molecular weight and protein binding on pharmacokinetic-pharmacodynamic (PK-PD) descriptors has not been established for most NMBA.

METHODS

The octanol/phosphate buffer distribution coefficients (logD) of various NMBA (vecuronium, rocuronium, mivacurium isomers (cis-cis, cis-trans and trans-trans), doxacurium, cisatracurium, atracurium, succinylcholine) were determined. The free fraction for each drug was measured using an ultrafiltration technique. PK-PD descriptors were obtained from selected clinical studies. Correlations between physicochemical parameters (including molecular weight) and PK-PD descriptors were assessed by linear or multiple linear regression.

RESULTS

A wide range of log D (-4.15 for succinylcholine to 0.75 for vecuronium) and free fraction (from 31% for vecuronium to 80% for succinylcholine) is observed for NMBA. Molecular weight combined with either lipid solubility (r2=0.70; P=0.001) or free fraction (r2=0.84; P<0.001) were highly correlated with potency, while for k(e0) a greater degree of correlation was obtained when both lipid solubility and free fraction (r2=0.74; P=0.002) were included.

CONCLUSIONS

The basic characteristics of NMBAs, namely, molecular weight, lipid solubility and protein binding, are strongly associated with the kinetics of the drug response.