Numerous local anesthetics have an asymmetric tetrahedron carbon, which confers stereoselective differences between the isomers. The authors attempted to quantify the depressant effect of racemic bupivacaine, levobupivacaine, and ropivacaine on myocardial ventricular conduction and on myocardial contractility.


The authors studied the pharmacokinetics (outflow concentration) and pharmacodynamics (QRS widening) of the three drugs infused in an isolated rabbit heart preparation. All data were fitted simultaneously with use of mixed-effect modeling, thus allowing precise statistical comparison between the three drug parameters. The rate dependence of QRS widening was fitted separately.


Racemic bupivacaine, levobupivacaine, and ropivacaine induced a calculated maximum increase in QRS duration in the ratio 1:0.4:0.3. Css50, the dose which caused half the maximum increase in QRS duration at steady state, was similar for all three drugs (22 micrometer free concentration). A rate dependence of QRS widening was observed, which was in the ratio 1:0.5:0.25 for racemic bupivacaine, levobupivacaine, and ropivacaine, respectively.


In the isolated rabbit heart, racemic bupivacaine, levobupivacaine, and ropivacaine induce an increase in QRS duration in the respective ratio of 1:0.4:0.3, which was rate dependent in approximately the same ratio.


Click on the links below to access all the ArticlePlus for this article.

Please note that ArticlePlus files may launch a viewer application outside of your web browser.

NUMEROUS local anesthetics, including bupivacaine, have an asymmetrically substituted carbon, which confers stereoselective differences between the enantiomers. 1Bupivacaine binding to sodium channels and to serum proteins is stereoselective, 2,3and the levo-(S −)-bupivacaine enantiomer is less cardiotoxic than the dextro-(R +)-enantiomer. 4–7However, new experiments strongly suggest that this stereospecific binding to the sodium channel is less important at the nerve site than at the heart site, 8thus explaining that levobupivacaine may cause nerve block of similar or greater intensity and duration than the racemic mixture. 9–11Similar observations have been made about R  (−)- and S  (+)-mepivacaine. 9These findings indicate that, among local anesthetic stereoisomers, some are safer than their mirror enantiomer and than the usual racemic mixture. Ropivacaine ([S −]1-propyl-2′,6′-pipecoloxylidide) is the only local anesthetic available as a pure enantiomer. Ropivacaine is believed to be safer than bupivacaine. 12,13Ropivacaine, which has one carbon less and is slightly less lipophilic than bupivacaine, is also slightly less potent than bupivacaine. 14QRS widening, observed with lidocaine, 15is related to ventricular conduction velocity slowing. 16Long-acting local anesthetics, such as bupivacaine, also impair ventricular conduction, primarily by blocking voltage-sensitive sodium channels, 1,17and this effect is more pronounced with bupivacaine than with lidocaine because of the rate dependence of the block. 18Moreover, stereospecificity has been observed. 1,4–6,17These effects on Na+channels may cause life-threatening arrhythmias, which are thought to be enhanced by an effect on K+channels; however, stereospecificity may vary depending on the drug and the channel studied. 19–21 

To compare racemic bupivacaine, levobupivacaine, and ropivacaine, we attempted to quantify the depressant effect of these agents on myocardial ventricular conduction. Because bupivacaine toxicity has been shown to be rate dependent, 22we also studied the effect of heart rate on QRS duration changes induced by the three drugs.

Materials and Methods

We studied the effects of local anesthetics on an isolated rabbit heart model with use of a modification of our previously described procedure. 6,15Twenty-one male New Zealand rabbits, weighing 1,550–2,040 g, were studied in a random block design of three groups of seven animals each. This study was approved by our institutional animal care committee. Care of the animals conformed to the recommendation of the Helsinki Declaration and to the guidelines of the European Communities and French laws for animal experiments (accreditation No. 1989/2559 to Dr. Mazoit). Group 1 animals were infused with racemic bupivacaine, group 2 animals were infused with levobupivacaine, and group 3 animals were infused with ropivacaine. The experimenters were blind to the drug used until study completion. Nine control animals also were studied at random to ensure the stability of the preparation. Finally, we incorporated a pilot group of five previously studied rabbit hearts that had been infused with racemic bupivacaine; the results of this group were incorporated in the rate dependence part of the study with use of the interoccasion variability concept (see Statistics).

Drugs and Reagents

The drugs used were the commercial solutions for racemic bupivacaine and ropivacaine (ASTRA France, Paris, France). The solutions were tested for concentration accuracy with use of hydrochloride salts (ASTRA Pain Control, Södertalje, Sweden). Levobupivacaine hydrochloride was a gift from Chiroscience (Cambridge, United Kingdom). Racemic bupivacaine was tested with use of another hydrochloride monohydrate salt, from Sigma (St. Quentin Fallavier, France). After a concentration check, blind stock solutions of the drugs were prepared.

The same buffer, with the following composition, was used throughout the study: 118 mm NaCl, 4.7 mm KCl, 2.5 mm CaCl2, 1.2 mm MgSO4, 1.2 mm KH2PO4, 25 mm NaHCO3, 5.5 mm glucose, and 2.0 mm Na pyruvate. The pH of the perfusate at heart inflow was maintained between 7.37 and 7.42. Reagents for chromatography and salts for buffer were purchased from Prolabo (Paris, France).

Study Procedure

The rabbits were anesthetized with 6 mg/kg pentobarbital intraperitoneally. Tracheotomy was performed, and the animals were ventilated manually. The chest was opened, and, after intravenous heparin injection, the heart was removed and mounted quickly on a nonrecirculating Langendorff apparatus, and the coronary arteries were perfused via  the aorta at a constant flow of 30 ml/min with use of a modified Krebs-Henseleit buffer bubbled with a mixture of 95% oxygen and 5% carbon dioxide at 37°C. The hearts were paced atrially throughout the study with a bipolar electrode at 210 beats/min using a Chronocor IV stimulator (Telectronics, SOREM Presles en Brie, France), which delivered a square pulse of 3.5 mA. We used the following exclusion criteria for the preparation: (1) the presence of aortic valve regurgitation, (2) a rhythm (before pacing) less than 120 beats/min or greater than 170 beats/min, (3) the presence of arrhythmias, and (4) a dP/dt maximum lower than 1,000 mmHg/s. After an 8- to 12-min stabilization period, the drug (racemic bupivacaine, levobupivacaine, or ropivacaine) was infused into the inflow perfusate at 20 μm for 5 min (from T0 to T5) and at 5 μm during 15 min (from T5 to T20) with use of a model 33 Harvard pump (Harvard, Les Ulis, France). The concentrations infused correspond to amounts of 0.6 μm/min and 0.15 μm/min, respectively (at 30 ml/min buffer infusion). The hearts were studied during a total period of 60 min. The outflow perfusate was sampled with use of a fraction collector at frequent intervals up to 60 min. Pharmacodynamic variables (electrocardiography and left ventricular pressure) were recorded at the end of each effluent sampling time. Racemic bupivacaine, levobupivacaine, and ropivacaine were measured with use of gas chromatography. Electrocardiography was measured with use of surface electrodes. Data were recorded on a Gould 8000s chart recorder (Gould, Les Ulis, France). QRS duration of three consecutive beats, recorded at a paper speed of 200 mm/s, were averaged. The rate of pacing was modified between 17 and 19 min and between 50 and 52 min, with use of steps every 10 s from 170 to 350 beats/min, in increments of 20 beats/min (the starting point of the sequence [initial heart rate] and the order of change [increase or decrease in stimulation rate] were chosen at random). This additional procedure was performed to quantify the rate dependence of QRS duration. A preliminary study performed on five rabbit hearts with use of a step increase of racemic bupivacaine at 0, 1.535, 3.07, and 6.14 μm showed that, when the pacing rate was changed, approximate steady state QRS duration was attained after 4 or 5 s. A 10-min stabilization period was observed between changes in dose amount. The results of this preliminary group have been incorporated into the study (see after).

Pharmacokinetics and Pharmacodynamics.

We used the model previously described. 6,15Briefly, linear pharmacokinetics were assumed, and the heart was described by a two-compartment open model (with the assumption of venous equilibrium). If C0is drug concentration in the inflow perfusate, the outflow perfusate concentration (C) can be expressed as a function of time (t) by the following relation:


where k10is the exit rate constant from the central compartment, A1is the amount of drug in the central compartment at time t, and Q is the perfusate flow. Fitting was performed with use of standard equations 23with the procedure ADVAN 3 from the program NONMEM 24(version V, level 1). The volume of the central compartment was not measurable with data from outside the heart. Therefore, we set all volumes to unity. We calculated the three rate constants, k12, k21, and k10, from the central to the peripheral compartment, from the peripheral to the central compartment, and from the central compartment to outflow, respectively.

The increase in QRS duration (E) was fitted simultaneously to the Emax:


where E0is the basal QRS duration, Emaxis the maximum increase in QRS duration, and Ai50is the drug amount in compartment i that produces half Emaxat steady state. 25A special-effect compartment model, as described by Sheiner et al. , 26was also tried. The steady state perfusate concentration that produced half Emax(Css50) was calculated as Css50= Ai50/k10, in case of an effect occurring in the central compartment.


Between-group comparison of QRS duration measured before and at the end of each infusion was performed using the Student t  test with the Bonferroni correction. Results are expressed as the arithmetic mean and SD, except for figure 3, in which the standard error of the mean was used for clarity of the figure.

The data were fitted using the program NONMEM. Its use permits the fitting of mixed-effects models by using two levels of random errors (intraindividual and interindividual variability). By using nested models, it allows testing of the statistical difference between drugs for a specified parameter in a wide range of experiments. 27Extended least squares were used as measure of goodness-of-fit 28(see appendix in Web enhancement).

The concentration–time data first were fitted with use of the first-order method. The parameter estimates obtained at this step were used then (fixed at the value obtained at this step) for the calculation of the effect parameter estimates using the whole data set. The choice between the different pharmacokinetic (PK) models (one or two compartments) and pharmacodynamic (PD) models (effect in the central, peripheral, or special-effect compartment) was made with use of the Akaike criterion. 29Thus, a full model, with all interindividual variability parameters considered to be relevant, was defined (see appendix in Web enhancement). After this full model was defined, the choice between the full model and successive reduced models was made with use of the log-likelihood ratio test. 30To avoid overparametrization, we only considered parameters with an estimated coefficient of variation (CV) of the estimates that was less than 60%. 31,32Intraindividual variability (assay error, model mispecification, and so forth) was modeled using a combined constant CV and additive error model without interaction for pharmacokinetic and an additive error model for pharmacodynamic. Interindividual variability was modeled as θexp(η) (assuming a log-normal distribution), in which θ is the vector of the fixed-effect parameter and η is the vector of interindividual variability, with variance ω2. We used the hybrid method, with the mean η corresponding to pharmacokinetic parameters set to 0. We assumed no covariance either between the elements of ε, the vector of residual error resulting from intraindividual and measurement variability, or between the elements of η and the elements of ε.

The rate dependence of QRS duration was tested with use of the following linear model: where HR is heart rate and Ciis the concentration of drug i in perfusate. The slope was considered linearly dependent on drug concentration and, therefore, was set as the product of an intrinsic parameter for drug i (Slopei) times the concentration of that drug (Ci). Both parameters (Intercept and Slopei) had a fixed effect and a random effect (modeled as exp(η)) component. A constant CV was used to model the residual intraindividual error. Data from a pilot study were incorporated in the fitting procedure, considering interoccasion variability with a different η and ε.


QRS duration was constant throughout the study period in all control hearts. Local anesthetic infusion was followed by rapid QRS widening. Arrhythmias occurred during the rapid infusion phase (3–6 min after infusion initiation) and at the time of discontinuation of drug infusion (table 1). Because of the small number of hearts used, it was impossible to compare the number of arrhythmias that occurred with each drug. In fact, onlyone heart in the racemic bupivacaine group, four hearts in the levobupivacaine group, and four hearts in the ropivacaine group had no arrhythmias (table 1). Four hearts experienced arrhythmias at infusion discontinuation.

The maximum observed increase in QRS was significantly greater with racemic bupivacaine than with the two other drugs at the end of the two infusion phases (i.e. , at T5 and T20) (table 1).

Fitting was adequate with use of the two-compartment open model and the QRS widening effect located in the central compartment, i.e. , with no delay between outflow concentration and effect (table 2and fig. 1). Table 2shows the statistical difference between models. Racemic bupivacaine and levobupivacaine had a similar k21 and k10. Emaxwas significantly different among the three drugs, whereas Css50was similar for the three drugs (tables 2 and 3). Emaxwas more than twice as much for racemic bupivacaine than for levobupivacaine and approximately 4 times as much for racemic bupivacaine than for ropivacaine (table 3and fig. 2). The approximate Emaxratio was 1:0.4:0.25 for racemic bupivacaine, levobupivacaine, and ropivacaine, respectively. QRS widening showed a marked rate dependence, which was linearly related to dose amount, at least for racemic bupivacaine (table 3and fig. 3). The rate dependence of QRS widening (slope of the QRS duration–heart rate relation) was significantly different between the three drugs, with an approximate ratio of 1:0.5:0.2 for racemic bupivacaine, levobupivacaine, and ropivacaine, respectively.


The current study confirms that levobupivacaine and ropivacaine induce a much lower impairment of intraventricular conduction than does racemic bupivacaine. The number of hearts used in each group does not allow statistical comparison of the number and type of arrhythmias. However, either block (defined as the absence of electrical activity during at least one cardiac cycle, despite pacing) or premature ventricular contraction was more frequent in the racemic bupivacaine group (table 1). These arrhythmias are usually described with use of bupivacaine and are associated with decreased intraventricular conduction velocity and reentry phenomenon. 22This result is in accordance with the fact that racemic bupivacaine significantly leads to a greater impairment of intraventricular conduction and to a higher rate dependence. The use of nonlinear mixed-effect modeling allowed us to show that the theoretical maximum effect on QRS duration was significantly greater with levobupivacaine than with ropivacaine, whereas Css50was similar for all three drugs. Therefore, at similar free concentrations in blood, the three drugs are expected to induce ventricular conduction impairment in the approximate ratio (intrinsic activity) of 1:0.4:0.25 for racemic bupivacaine, levobupivacaine, and ropivacaine, respectively.

We used the global mixed-effect modeling technique, rather than the classic two-stage method because the former approach permits more accurate calculation of confidence intervals for parameter estimates than does the latter. 31,32Two assumptions were made for interindividual variability estimation. First, we assumed a log-normal distribution for pharmacokinetic and pharmacodynamic parameters, and, therefore, we modeled interindividual variability as exp(η). Second, we assumed a nonlinear behavior for pharmacodynamic data. Therefore, we used the hybrid method in NONMEM.

Fitting of pharmacokinetic data was adequate with use of the well-stirred, two-compartment model. 6,15,33The three drugs showed significantly different kinetic parameters. However, these differences are relatively minor, as shown in figure 2. At the time of discontinuation of drug administration, myocardial washout was rapid, even in the racemic bupivacaine group. These results are in accordance with our previous studies, which showed that bupivacaine did not accumulate in the myocardium 6,33and that the toxic effect of long-lasting local anesthetics was not the consequence of drug accumulation in tissue.

We used the simple Emaxmodel for QRS widening fitting because the addition of a sigmoid parameter resulted in overparametrization. Fitting was adequate, but a small bias that caused underestimation of the highest QRS values was observed (fig. 2). The Css50(the inflow or outflow steady state concentration that produced half the maximum effect) was similar for all three drugs. We have already shown that the Css50was similar between racemic bupivacaine and lidocaine and between racemic bupivacaine, levobupivacaine, and R  (+)-bupivacaine. The Css50for the three drugs was similar to the value previously reported (43 μm vs.  29–39 μm in the current study and previous studies, respectively). 6,33The calculated Css50(43 μm, i.e. , approximately 14 to 15 μg/ml) needs to be interpreted carefully because we used a protein-free perfusate solution. We may estimate that the approximate free concentration that is necessary to double the basal QRS duration at 210 beats/min was 2.4, 7.2, and 14.4 μg/ml for racemic bupivacaine, levobupivacaine, and ropivacaine. 3These concentrations are in the range of the free concentrations expected to occur during accidental massive intravenous injection for racemic bupivacaine, but they are likely more than those expected during the same complication for levobupivacaine and ropivacaine. 34 

All three drugs showed a marked rate dependence of QRS widening (table 3and fig. 3). This rate dependence was statistically different between the three drugs (table 2). With racemic bupivacaine, the slope of the relation between heart rate and QRS duration was related linearly to the dose within the range of frequencies used (fig. 3, top), and nothing indicated that this phenomenon is different with the other two local anesthetics. Therefore, for the comparison between racemic bupivacaine, levobupivacaine, and ropivacaine, we modeled QRS duration as a linear function of inflow concentration with use of equation 3. Fitting was adequate (fig. 3, bottom). The rate dependence of QRS widening was in the range of 1:0.5:0.2 for racemic bupivacaine, levobupivacaine, and ropivacaine, respectively, which approximates the ratio of Emaxcalculated for these drugs. The dose–effect curve parameters (Emaxand Css50) were calculated at a fixed frequency of 210 beats/min. Because QRS duration was related linearly to drug concentration and heart rate, changes in heart rate might change Emax, with a fixed ratio between drugs. However, it may be reasonably speculated that even drugs that rapidly dissociate from the receptor may increase intrinsic activity at extreme heart rates. In contrast, Css50is not expected to vary with changes in heart rate.

In conclusion, using mixed-effect modeling, we showed that racemic bupivacaine, levobupivacaine, and ropivacaine block intraventricular conduction in the rabbit heart in the respective ratio of 1:0.4:0.25. This impairment of conduction was rate dependent in approximately the same ratio (1:0.5:0.2). This must be interpreted with consideration of the supposed slightly lower nerve block potency of ropivacaine when compared with levobupivacaine and with racemic bupivacaine. 14Bupivacaine toxicity is rare, but its potential life-threatening effect needs to be taken into account. With these conditions, extrapolation of our results to humans leads to the conclusion that levobupivacaine and ropivacaine are safer than racemic bupivacaine. The choice between levobupivacaine and ropivacaine necessitates further investigation, particularly to compare toxicity with nerve-blocking potency.

The authors thank Mrs. Régine le Guen (Kremlin-Bicêtre College of Medicine, Paris-Sud University, France) for her technical assistance and Dr. Genery from ChiroScience (Cambridge, United Kingdom) for the gift of levobupivacaine.


Clarkson CW: Stereoselective block of cardiac sodium channels by RAC109 in single guinea pig ventricular myocytes. Circ Res 1989; 65: 1306–23
Lee-Son S, Wang GK, Concus A, Crill E, Strichartz GR: Stereoselective inhibition of neuronal sodium channels by local anesthetics: Evidence for two sites of action? A nesthesiology 1992, 77: 324–35
Mazoit JX, Cao LS, Samii K: Binding of bupivacaine to serum proteins, isolated albumin and isolated alpha1-acid glycoprotein: Differences between the two enantiomers are partly due to cooperativity. J Pharmacol Exp Ther 1996; 256: 109–15
Vanhoute F, Vereecke J, Verbeke N, Carmeliet E: Stereoselective effects of the enantiomers of bupivacaine on the electrophysiological properties of the guinea-pig papillary muscle. Br J Pharmacol 1991; 163: 1275–81
Denson DD, Behbehani MM, Gregg RV: Enantiomer-specific effect of an intravenously administered arrhythmogenic dose of bupivacaine on neurons of nucleus tractus solitarius and the cardiovascular system in the anesthetized rat. Reg Anesth 1992; 17: 311–6
Mazoit JX, Boïco O, Samii K: Myocardial uptake of bupivacaine, I: Pharmacokinetics and pharmacodynamics of bupivacaine enantiomers in the isolated perfused rabbit heart. Anesth Analg 1993; 77: 477–82
Bardsley H, Gristwood R, Baker H, Watson N, Nimmo W: A comparison of the cardiovascular effects of levobupivacaine and rac-bupivacaine following intravenous administration to healthy volunteers. Br J Clin Pharmacol 1998; 46: 245–9
Nau C, Vogel W, Hempelmann G, Brau ME: Stereoselectivity of bupivacaine in local anesthetic-sensitive ion channels of peripheral nerve. A nesthesiology 1999; 91: 786–95
Aberg G: Toxicological and local anaesthetic effects of optically active isomers of two local anaesthetic compounds. Acta Pharmacol Toxicol (Copenh) 1972; 31: 273–86
Akerman B: Uptake and retention of the enantiomers of a local anaesthetic in isolated nerve in relation to different degrees of blocking of nervous conduction. Acta Pharmacol Toxicol (Copenh) 1973; 32: 225–36
Dyhre H, Lang M, Wallin R, Renck H: The duration of action of bupivacaine, levobupivacaine, ropivacaine and pethidine in peripheral nerve block in the rat. Acta Anaesthesiol Scand 1997; 41: 1346–52
Scott DB, Lee A, Fagan D, Bowler GM, Bloomfield P, Lundh R: Acute toxicity of ropivacaine compared with that of bupivacaine. Anesth Analg 1989; 69: 563–9
Pitkanen M, Feldman HS, Arthur GR, Covino BG: Chronotropic and inotropic effects of ropivacaine, bupivacaine, and lidocaine in the spontaneously beating and electrically paced isolated perfused rabbit heart. Reg Anesth 1992; 17: 183–92
Polley LS, Columb MO, Naughton NN, Wagner DS, van de Ven CJ: Relative analgesic potencies of ropivacaine and bupivacaine for epidural analgesia in labor: Implications for therapeutic indexes. A nesthesiology 1999; 90: 944–50
Mazoit JX, Kantelip JP, Orhant EE, Talmant JM: Myocardial uptake of lignocaine: Pharmacokinetics and pharmacodynamics in the isolated perfused heart of the rabbit. Br J Pharmacol 1990; 101: 843–6
Anderson KP, Walker R, Lux RL, Ershler PR, Menlove R, Williams MR, Krall R, Moddrelle D: Conduction velocity depression and drug-induced ventricular tachyarrhythmias: Effects of lidocaine in the intact canine heart. Circulation 1990; 81: 1024–38
Valenzuela C, Snyders DJ, Bennett PB, Tamargo J, Hondeghem LM: Stereoselective block of cardiac sodium channels by bupivacaine in guinea pig ventricular myocytes. Circulation 1995; 92: 3014–24
Clarkson CW, Hondeghem LM: Evidence for a specific receptor site for lidocaine, quinidine, and bupivacaine associated with cardiac sodium channels in guinea-pig ventricular myocardium. Circ Res 1985; 56: 496–506
Castle NA: Bupivacaine inhibits the transient outward K+ current but not the inward rectifier in rat ventricular myocytes. J Pharmacol Exp Ther 1990; 255: 1038–46
Valenzuella C, Delpon E, Tamkun MM, Tarmago J, Snyders DJ: Stereoselective block of a human cardiac potassium channel (Kv1.5) by bupivacaine enantiomers. Biophys J 1995; 69: 418–27
Valenzuella C, Delpon E, Franqueza L, Gay P, Snyders DJ, Tamargo J: Effects of ropivacaine on a potassium channel (hKv1.5) cloned from human ventricle. A nesthesiology 1997; 86: 718–28
de La Coussaye JE, Brugada J, Allessie MA: Electrophysiologic and arrhythmogenic effects of bupivacaine: A study with high-resolution ventricular epicardial mapping in rabbit hearts. A nesthesiology 1992; 77: 132–41
Gibaldi M, Perrier D: Pharmacokinetics, 2nd edition. New York, Marcel Dekker, 1982
Scheiner LB, Beal SL, Boeckmann A: NONMEM Users Guides 1–6. Regents of the University of California, 1979–1999
Holford NH, Sheiner LB: Understanding the dose-effect relationship: Clinical application of pharmacokinetic-pharmacodynamic models. Clin Pharmacokinet 1981; 6: 429–53
Sheiner LB, Stanski DR, Vozech S, Miller RD, Ham J: Simultaneous modeling of pharmacokinetics and pharmacodynamics: Application to d-tubocurarine. Clin Pharmacol Ther 1979; 25: 358–71
Moine P, Mazoit JX: Streptococcus pneumoniae pneumonia in mice: Optimal amoxicillin dosing predicted from a pharmacokinetic-pharmacodynamic model. J Pharmacol Exp Ther 1999; 291: 1086–92
Sheiner LB, Beal SL: Pharmacokinetic parameter estimates from several least squares procedures: Superiority of extended least squares. J Pharmacokinet Biopharm 1985; 13: 185–201
Yamaoka K, Nakagawa T, Uno T: Application of Akaike’s information criterion (AIC) in the evaluation of linear pharmacokinetic equations. J Pharmacokinet Biopharm 1978; 6: 165–75
Eadie WT, Drijard D, James FE, Roos M, Sadoulet B: Statistical Methods in Experimental Physics. Amsterdam, North-Holland, 1982, pp 471–84
Sheiner LB, Beal SL: Evaluation of methods for estimating population pharmacokinetics parameters, I: Michaelis-Menten model: Routine clinical pharmacokinetic data. J Pharmacokinet Biopharm 1980; 8: 553–71
Sheiner BL, Beal SL: Evaluation of methods for estimating population pharmacokinetic parameters, II: Biexponential model and experimental pharmacokinetic data. J Pharmacokinet Biopharm 1981; 9: 635–51
Mazoit JX, Orhant EE, Boïco O, Kantelip JP, Samii K: Myocardial uptake of bupivacaine, I: Pharmacokinetics and pharmacodynamics of lidocaine and bupivacaine in the isolated perfused rabbit heart. Anesth Analg 1993; 77: 469–76
Feldman HS, Arthur GR, Covino BG: Comparative systemic toxicity of convulsant and supraconvulsant doses of intravenous ropivacaine, bupivacaine, and lidocaine in the conscious dog. Anesth Analg 1989; 69: 794–801