A New Pharmacokinetic Model of Propofol for Japanese Patients

Target-controlled infusion (TCI) of propofol is used for general anesthesia. However, the only pharmacokinetic parameter commercially used for Japanese patients is weight, and pharmacokinetic models are based on European physical attributes. Drug metabolism also differs in races. This study aimed to identify optimal continuous doses of propofol for Japanese patients and to create a simulated pharmacokinetic (PK) model. Thirty Japanese patients were enrolled. Patients received a constant infusion of 9 mg/kg/h of propofol. Arterial blood samples were collected and the time course of plasma propofol concentrations was modeled using the nonlinear mixed effects model (NONMEM) three-compartmental PK model. We validated the model by intravenously injecting 10 patients with a TCI driver system programmed with the NONMEM model. Our model’s performance was evaluated using the median prediction error (MDPE), median absolute prediction error (MDAPE), and Wobble. We analyzed 320 blood samples for model building and 160 samples for validating our new model. The calculated parameters for the three-compartmental PK model were volume [V1, 3.58; V2, 13.0 + 0.49 × (Age—64); and V3, 186] and elimination clearance [CL1, 0.77 + (WT—54) × 0.04 + (HT—158) × 0.03; CL2, 0.89 + 0.12 × (Age—64); and CL3, 0.98 × exp ((Age—64)/10)]. The new model improved MDPE, MDAPE, and Wobble values (11.5% ± 43.8%, 14.3% ± 33.0%, and 25.0% ± 21.3%, respectively). We created a new pharmacokinetic model for Japanese patients, which is more accurate than the three existing models applied to Japanese populations. Electronic document is a “live” template. The various components of your paper [title, text, heads, etc.] are already defined on the style sheet, as illustrated by the portions given in this document.


Introduction
Target-controlled infusion (TCI) of propofol is used worldwide to induce and maintain general anesthesia. However, the pharmacokinetic parameters of TCI devices for propofol that are commercially used in Japan include only the patient's weight. Various pharmacokinetic parameters have been published; however, these models are based on European physical attributes, which differ considerably with regard to fat and muscle composition from those of Japanese populations [1] [2]. In addition, drug metabolism differs in races [3]. Using only a European standard of weight may cause bias in the use of a TCI device, which could potentially lead to overdoses in Japanese patients. Furthermore, this bias is likely remarkable, particularly when TCI devices are used for obese patients.
Tachibana et al. [4] reported that the body mass index (BMI) influences blood propofol concentrations, and may lead to a propofol overdose in obese patients when assessing propofol concentrations by using the Marsh model. Lampotang et al. [5] reported that the plasma and effect-site concentrations of propofol predicted using the Marsh model to achieve loss of consciousness vary among Caucasians, Chinese, Africans, and Indians. The present study aimed to identify optimal continuous doses of propofol for Japanese patients and create a simulated PK model using NONMEM version 7.2 software (ICON Clinical Research, LLC, North Wales, PA, USA) [6].

Methods
The protocol of this study was designed in accordance with the principles out-

Recruitment of Patients
This was a prospective observational study. Thirty patients living in Japan from birth were enrolled in this study. All patients were scheduled for elective brain, head, or neck surgery under general anesthesia at Sapporo Medical University Hospital. The exclusion criteria were as follows: age < 20 years or age > 70 years; American Society of Anesthesiologists (ASA) physical status of 3 to 5; allergic to propofol; a daily intake of psychoactive drugs; known or suspected drug or alcohol abuse; and cardiac, hepatic, renal, or neurological impairment; extremely obese or lean body.

Creating a Pharmacokinetic Model for Propofol
We planned to obtain blood samples from 20 patients to create a pharmacokinetic model of propofol. The patients were not premedicated before arriving at the operating room where cardiac parameters were monitored with electrocar-Open Journal of Anesthesiology diography and brain parameters were monitored with the bispectral index (BIS) (VISTA; Covidien, Dublin, Ireland). Blood pressure, oxygen saturation, and end-tidal carbon dioxide concentrations were also monitored. A 20-gauge cannula was inserted into a forearm vein, and a 22-gauge radial artery catheter was also inserted for arterial blood sampling and invasive blood pressure monitoring.

Propofol Assays
Propofol concentrations were measured using a high-performance liquid chromatography (HPLC) system with a fluorescence detector (EICOM, Kyoto, Japan). To confirm the propofol concentration, we used previously described HPLC methods as the basis for the studies we conducted [7].
In brief, 200 µL of a blood plasma sample was mixed with 200 μL of methanol containing 10 μL of 2-tert-butyl-6-methylphenol (0.1 μL/mL) as the internal standard. After mixing, the samples were centrifuged for 10 min at 15,000 × g.
The supernatant was filtered and injected into the high-performance liquid chromatograph for analysis. The propofol concentration was calculated, based on the ratio of the propofol and peaks in the HPLC chromatogram.
The mobile phase comprised 0.1 M phosphate buffer and methanol (3:7). The flow rate was 0.5 mL/min, and the column temperature was 25˚C. The fluorescence detector was set at an excitation of 270 nm and an emission of 310 nm.
Drug separation was achieved with a C18 column (4.6 mm × 30 mm; PP-ODS II; EICOM), which was maintained at 25˚C.

Comparison with Published Pharmacokinetic Models
The Marsh model [8], Schnider model [9], and General Purpose model [10] for patient simulation were compared using TIVA Trainer software (Gutta B.V., Aerdenhout, Netherlands). The simulation program calculated the predicted blood propofol concentrations (Cp pred

Model Validation
We validated our model, as follows. A TCI driver programmed with our new model intravenously injected 10 patients with propofol to achieve a blood concentration (Cp) of 3.0 μg/mL and collected arterial blood samples at the time points described previously. To administer propofol, we used our original program. A syringe pump was connected via RS232C terminal to the computer.
The performance of the model in each patient was assessed using the MDPE, MDAPE, and Wobble. Blood pressure degradation after achieving a Cp of 3.0 μg/mL was measured, and adverse effects such as vascular pain and hypotension were evaluated as the secondary outcomes.  We revised the model by including covariates for age, height, and body weight to optimize the relationship between the observed data, predicted data, and pharmacokinetic data. Table 3 shows that the pharmacokinetic data derived from the new model and the Schnider model were similar. Age was a significant covariate for rapid peripheral volume and rapid peripheral clearance, whereas height and body weight were significant covariates for metabolic clearance.

Results and Discussions
We then evaluated this new model using a TCI driver programmed with the new model to achieve a Cp of 3.0 μg/mL. Table 4 shows the patient demographics, average elapsed time to consciousness, and blood pressure degradation. Figure 2 shows the percentage prediction error for the new model and Table 5 shows the pharmacokinetic performance of intravenously injected propofol. Many studies have used the pharmacokinetic parameters established by Marsh et al. [8], and have assessed these parameters in obese, pediatric, and elderly Open Journal of Anesthesiology   General Purpose (patient) −40.3 ± 30.0 40.7 ± 28.9 11.3 ± 6.9 The data are presented as the mean ± the standard deviation. MDAPE: median absolute prediction error, MDPE: median prediction error. Open Journal of Anesthesiology   Patients with vascular pain (n) 0/10 The data are presented as the mean ± the standard deviation or as the number. Cp: predicted blood propofol concentration, F: female, M: male, MAP: mean arterial pressure. patients. These studies demonstrated that the parameters were clinically acceptable. Li et al. [12] reported that TCI systems using the Marsh model are clinically acceptable for upper abdominal surgery. Some investigators have found that other models are preferable to the Marsh model. Masui et al. [13] recommend the Schnider model rather than the Marsh model. Rigouzzo et al. [14] report that the Schnider model is more useful for children. Glen et al. [15] recommend the White model rather than the Marsh model. However, our results indicated that the Marsh model performed better than the Schnider model. The performance of the General Purpose model was the poorest among the three models in the present study, despite its application in most published data. An MDPE less than 15% and an MDAPE less than 30% indicate that a TCI sys-Open Journal of Anesthesiology tem has a clinically acceptable performance [16]. We took into consideration that the pharmacokinetic parameters may have been biased in published studies because most of their data are based on European physical attributes.
Many factors such as cardiac output and hepatic blood flow influence blood propofol concentrations. In addition, weight and body height, which are affected by ethnic differences, can introduce bias. Lampotang et al. [5] reported that propofol concentrations predicted by the Marsh model varied among ethnicities, and created a pharmacodynamic model of propofol-induced loss of consciousness. Our model outperformed the other three models when applied to a Japanese population and may be useful. Moreover, the costs will be reduced because unnecessary drug administration may be cut when our model is used.
Sex may also affect propofol pharmacokinetics. The present study included more women than men. Vuyk et al. [17] found that sex affects propofol pharmacokinetics at lower concentrations in female elderly patients than in male elderly patients. The age of our patients ranged from 30 years to 70 years. We found comparable sex-related differences in propofol concentrations. However, no significant difference was observed.
Vascular pain on injection is one disadvantage of propofol. Eriksson et al. [18] reported that pH changes on the vascular endothelium may modify propofol-induced pain on injection and that lignocaine mixed with propofol causes less pain. None of our patients experienced vascular pain; however, we may have injected a smaller bolus for induction than that used in the Marsh model. If so, the pH may not have affected the vascular endothelium as much as it does in the Marsh model.
Hypotension after propofol injection is also a concern. Wang et al. [19] found that the mean arterial blood pressure decreased by approximately 18% when achieving an effect-site concentration of 3.0 μg/mL by using a modified Marsh model. We believe that our model achieved a similar performance.
The findings of our study must be evaluated in consideration of some limitations. First, our sampling schedule did not include data for the time after the termination of propofol infusion. These data may be needed for building an exact PK model. Potential complications may be over dose. When the patient's metabolism and excretion are different from usual, propofol administration used by our model may be potency of overdose administration. However, the purpose of this article was to create a new model suitable for Japanese individuals, and to determine how to introduce and maintain anesthesia suitably for Japanese individuals. To that end, we believe that a stable blood concentration is an important serious matter during the introduction and maintenance periods. We believe that the change in the blood concentration from the dosage cancellation depends on the metabolic pathway of an individual, even if we use any models. Therefore, it was not the principal objective in this study. We also believe it can be substituted in the elapsed time to consciousness, which is a factor of pharmacodynamics but not PK.
We could not calculate the pharmacodynamic parameters of propofol because general anesthesia was induced using fentanyl. In addition, fentanyl may affect vascular pain. We postulate that the effect of fentanyl would be small because propofol was injected immediately after the infusion of fentanyl. In addition, we did not have a sufficient number of blood samples to apply our model clinically.
More samples are needed to increase the accuracy of propofol pharmacokinetic parameters.

Conclusion
In conclusion, we demonstrated that the General Purpose model, Schnider model, and Marsh model cannot be directly applied to Japanese patients who are infused with a constant rate of propofol. We created a new pharmacokinetic model that performed better than the other three models for Japanese populations.