Peripheral link model as an alternative for pharmacokinetic-pharmacodynamic modeling of drugs having a very short elimination half-life

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.

Control of Muscle Relaxation During Anesthesia: A Novel Approach for Clinical Routine

During general anesthesia drugs are administered to provide hypnosis, ensure analgesia, and skeletal muscle relaxation. In this paper, the main components of a newly developed controller for skeletal muscle relaxation are described. Muscle relaxation is controlled by administration of neuromuscular blocking agents. The degree of relaxation is assessed by supramaximal train-of-four stimulation of the ulnar nerve and measuring the electromyogram response of the adductor pollicis muscle. For closed-loop control purposes, a physiologically based pharmacokinetic and pharmacodynamic model of the neuromuscular blocking agent mivacurium is derived. The model is used to design an observer-based state feedback controller. Contrary to similar automatic systems described in the literature this controller makes use of two different measures obtained in the train-of-four measurement to maintain the desired level of relaxation. The controller is validated in a clinical study comparing the performance of the controller to the performance of the anesthesiologist. As presented, the controller was able to maintain a preselected degree of muscle relaxation with excellent precision while minimizing drug administration. The controller performed at least equally well as the anesthesiologist.

Volume of the effect compartment in simulations of neuromuscular block

Background

The study examines the role of the volume of the effect compartment in simulations of neuromuscular block (NMB) produced by nondepolarizing muscle relaxants.

Methods

The molar amount of the postsynaptic receptors at the motor end plates in muscle was assumed constant; the apparent receptor concentration in the effect compartment is the ratio of this amount and the volume arbitrarily assigned to the effect compartment. The muscle relaxants were postulated to diffuse between the central and the effect compartment and to bind to the postsynaptic receptors. NMB was calculated from the free concentration of the muscle relaxant in the effect compartment.

Results

The simulations suggest that the time profiles of NMB and the derived pharmacokinetic and pharmacodynamic variables are dependent on the apparent receptor concentration in the effect compartment. For small, but not for large, volumes, times to peak submaximal NMB are projected to depend on the magnitude of NMB and on the binding affinities.

Conclusion

An experimental design to estimate the volume of the effect compartment is suggested.

Wirkortäquilibration, Anschlagzeit, „time to peak effect“

Contrary to the situation in "classical" clinical pharmacology, non-steady state phenomena play a fundamental role for clinical pharmacology in anesthesia. Their understanding is of tantamount importance for the safe and efficient application of drugs relevant to anesthesia. Concepts like optimised target-controlled infusion (TCI), effect compartment targeting and the small margin of error tolerable during maintained spontaneous ventilation, force the anesthesiologist to acquire a firm understanding of the difference between the concentration time course at the effect side vs. time course of the plasma concentration. The underlying concepts, their application for the rational use of muscle relaxants, propofol with TCI systems, volatile anaesthetics and opioids will be discussed.

Impact of peripheral elimination on the concentration-effect relationship of remifentanil in anaesthetized dogs

BACKGROUND

This study elucidates the impact of sampling site when estimating pharmacokinetic-pharmacodynamic (PK-PD) parameters of drugs such as remifentanil that undergo tissue extraction in the biophase. The interrelationship between the concentrations of remifentanil predicted for the effect compartment and those measured in arterial, venous, and cerebrospinal fluid were investigated under steady-state conditions.

METHODS

Following induction of anaesthesia with pentobarbital, an arterial cannula (femoral) and two venous catheters (jugular and femoral) were inserted. Electrodes were placed for EEG recording of theta wave activity. Each dog received two consecutive 5-min infusions for the PK-PD study and a bolus followed by a 60-min infusion was started for the steady-state study. Cerebrospinal fluid, arterial and venous blood samples were drawn simultaneously after 30, 40, and 50 min. At the end of the infusion, arterial blood samples were collected for pharmacokinetic analysis.

RESULTS

Remifentanil PK-PD parameters based on theta wave activity were as follows: apparent volume of distribution at steady-state (V(ss)) (231+/-37 ml kg(-1)), total body clearance (Cl) (63+/-16 ml min(-1) kg(-1)), terminal elimination half-life (t(1/2 beta)) (7.71 min), effect compartment concentration at 50% of maximal observed effect (EC(50)) (21+/-13 ng ml(-1)), and equilibration rate constant between plasma and effect compartment (k(e0)) (0.48+/-0.24 min). The mean steady-state cerebrospinal fluid concentration of 236 ng ml(-1) represented 52 and 74% of that in arterial and venous blood, respectively.

CONCLUSIONS

Our study re-emphasizes the importance of a sampling site when performing PK-PD modelling for drugs undergoing elimination from the effect compartment. For a drug undergoing tissue elimination such as remifentanil, venous rather than arterial concentrations will reflect more exactly the effect compartment concentrations, under steady-state conditions.

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.