(Almost) everything you learned about pharmacokinetics was (somewhat) wrong!

Pharmacokinetic-pharmacodynamic relationship of anesthetic drugs: from modeling to clinical use

Anesthesia is a combination of unconsciousness, amnesia, and analgesia, expressed in sleeping patients by limited reaction to noxious stimulations. It is achieved by several classes of drugs, acting mainly on central nervous system. Compared to other therapeutic families, the anesthetic drugs, administered by intravenous or pulmonary route, are quickly distributed in the blood and induce in a few minutes effects that are fully reversible within minutes or hours. These effects change in parallel with the concentration of the drug, and the concentration time course of the drug follows with a reasonable precision mathematical models based on the Fick principle. Therefore, understanding concentration time course allows adjusting the dosing delivery scheme in order to control the effects.   The purpose of this short review is to describe the basis of pharmacokinetics and modeling, the concentration-effects relationship, and drug interactions modeling to offer to anesthesiologists and non-anesthesiologists an overview of the rules to follow to optimize anesthetic drug delivery.

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.

Time-dependent cardioprotection with calcium antagonism and experimental studies with clevidipine in ischemic-reperfused pig hearts: part I

Clevidipine is a new ultrashort-acting dihydropyridine calcium antagonist developed for blood pressure regulation during cardiac surgery. When given locally to the ischemic and reperfused myocardium, clevidipine exerts a cardioprotective effect that varies depending on the timing of administration during ischemia. The current study explored the pharmacokinetics of clevidipine when administered locally into the coronary circulation. Pentobarbital-anesthetized pigs were randomly divided into four groups, of which three were subjected to myocardial ischemia through ligation of the left anterior descending coronary artery (LAD) for 15, 35, or 45 minutes (n = 4 in each group). The fourth group (n = 4) was not subjected to LAD occlusion and acted as nonischemic controls. Clevidipine (0.3 nmol/kg per minute) was infused over a period of 5 minutes through a catheter distal to the LAD ligation in the ischemic pigs, or at the corresponding site of the nonligated LAD in the nonischemic pigs. Following release of the LAD ligation, or in the nonligated group following the drug infusion, the pigs were subjected to 60 minutes of reperfusion. Simultaneous blood samples were obtained for analysis of clevidipine from the femoral artery and the coronary vein during reperfusion and during drug infusion in the nonischemic pigs. Blood samples for estimating the in vitro hydrolysis rate both in whole blood and in plasma were also obtained. In nonischemic hearts, clevidipine reached a steady state level in the coronary venous blood of about 30 nM during the infusion. The concentration declined almost to the detection limit (1 nM) within 3 minutes of the end of infusion. The mean blood clearance of clevidipine was calculated to be 0.17 l/min per kilogram, and the estimated half-life was 0.5 minute. In animals subjected to different periods of ischemia, very low levels of clevidipine were detected in the coronary venous blood only during the first 2 minutes of reperfusion. There were no detectable levels in the arterial blood at any time. Blood concentration profiles of clevidipine did not differ with the length of myocardial ischemia. The in vitro half-life in pig blood was 13 minutes, and the corresponding half-life in plasma was 111 minutes. At a dose known to exert cardioprotection when given as an intracoronary injection, the systemic concentration of clevidipine does not reach pharmacologically active levels. Clevidipine has an ultrashort blood half-life, and ischemic duration of up to 45 minutes does not seem to change the cardiac metabolism of this drug. Thus, it represents a pharmacological tool well suited for the study of time windows in cardioprotection. Moreover, considering the possibility of exerting myocardioprotection without any systemic effects, it could be an interesting compound to test in a clinical setting.

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.

Context-sensitive half-time and anesthesia: how does theory match reality?

The pharmacokinetic model, which describes context-sensitive decrement times, has received validation through direct measurement. Further validation has come from the ability of targeted infusion schemes to obtain reliable and stable concentrations of a variety of intravenous drugs. A complete understanding of the concept will enable the appropriate selection of drugs and an appreciation of the strengths and limitations of current drug delivery systems.