- Research article
- Open Access
- Open Peer Review
Pharmacokinetic and pharmacodynamic modeling of oral mitiglinide on glucose lowering in healthy Chinese volunteers
BMC Pharmacology and Toxicologyvolume 18, Article number: 54 (2017)
Mitiglinide is a widely used agent for diabetic treatment. We established a pharmacokinetic-pharmacodynamic (PK-PD) model to illustrate the relationship between mitiglinide plasma concentration and its glucose lowering effects in healthy volunteers.
The volunteers participated in the test after the administration of a single dose of 10 mg mitiglinide. The drug concentration in Plasma and the values of glucose levels were determined by LC-MS/MS assay and hexokinase method. A PK-PD model was established with a series of equations to describe the relationship between plasma medicine and glucose, and the equations were solved numerically and fitted to the data with the Phoenix NLME software.
The results of the two-compartment model analysis were based on the maximum likelihood criterion and visual inspection of the fittings. The terminal elimination half-life (t 1/2) was 1.69 ± 0.16 h and the CL/F was 7.80 ± 1.84 L/h. The plasma glucose levels began to decline by 0.2 h, and hit its bottom decreasing values of 2.6 mg/L at 0.5 h after administration. The calculated parameter and fitting curve indicated that the model established in our experiment fitted well.
A PK/PD model illustrates that the relationship between mitiglinide concentration in plasma and glucose lowering effect in healthy volunteers was established. The results of our experiment suggested that the model can be used reasonably to predict the relationship between PK and PD in mitiglinide, which could be used in diabetes mellitus dosage control in clinical trials and other fields.
Type 2 diabetes mellitus, which affects approximately 11% of adults in the U.S. and 8% worldwide , is resulted from a progressive decline in insulin action and the reduced response of β-cells to insulin [2, 3]. ATP-sensitive K+ (KATP) channels exists widely in pancreatic cells, cardiac, smooth and skeletal muscles, and central nervous system . Different KATP-channel may exist in different tissues, Kir6.2/SUR1 for β-cells, some neurons cardiac, Kir6.2/SUR2A for skeletal KATP-channels, and Kir6.2/SUR2B for vascular smooth muscle KATP-channels .
Mitiglinide is a rapid acting drug for diabetic drugs, which is selectively active on KATP channel of Kir6.2/SUR1 [6, 7]. Previous study has demonstrated that the preprandial administration of mitiglinide was benefit for postprandial hyperglycemia and glycemic control . A series of sensitive HPLC-MS/MS methods have been used to determine the pharmacokinetics (PK) properties of mitiglinide in rats  and humans. The PK property of 24 volunteers who received oral administration mitiglinide at a single dose from 5 to 20 mg/kg was studied . The maximum plasma concentration (C max ) in volunteers with mitiglinide from 5 to 20 mg/kg was 565.30 and 2625.26 ng/mL respectively, and the terminal half-life (T 1/2 ) was 1.19 and 1.43 h respectively. Pharmacokinetic–pharmacodynamic population model-based studies may suggest an accurate model to provide information for impact assessment of the treatment . The relationship between the pharmacokinetics of mitiglinide and the drug effects has not been studied in humans. The objective of present study is to measure the relationship between plasma drug concentration and its effects on healthy volunteers with oral administration. The experiment result may be used for prediction the therapeutic active of mitiglinide and its side effect profiles in clinic.
Mitiglinide (10 mg/tablet, lot no.14032011, expiration: 29/02/2016) was manufactured by Suzhong Pharma Co. Ltd., Jiangsu, China.
The study was approved by the Ethics Committee of the Affiliated Hospital of Nanjing University of Chinese Medicine (IRB#2014NL-020-02). Eighteen healthy male volunteers participated in this study, with average age of 24 (ranging from 18 to 40) and average weight of 67.5 kg (ranging from 56 to 80 kg). All volunteers went through thorough physical examinations and routine laboratory evaluations, including haematology, serum chemistry, virus and urinalysis. None of the volunteers in this study took any drugs which would have affected this study for at least 3 months beforehand. The limitation of our study lies in the fact that subjects were excluded from health problems including drug or alcohol abuse, and the abnormalities in the experiment. All volunteers were informed of the purpose and possible risks before the test and they could withdraw from the test any time for any reason. Signed informed consent was obtained before any trial-related activities.
Subjects were fasted for 10 h before the experiment. An intravenous cannula was inserted into the forearm vein to collect blood samples at time zero. The test group took the agent with 250 mL of water after baseline sampling. All subjects were accommodated in the care unit in the Affiliated Hospital of Nanjing University of Chinese Medicine during the test. Blood samples were taken at 0, 0.08, 0.17, 0.25, 0.33, 0.5, 0.75, 1, 1.5, 2, 3 h after drug administration, centrifuged immediately for 10 min at 3000 rpm, and stored at −80 °C until analysis.
Analysis of the plasma drug concentrations and blood glucose levels
Plasma mitiglinide concentration was determined by the modification of the methods of Yu and Takanohashi [9, 12]. A 1200 HPLC system (Agilent technologies, Palo Alto, CA, USA) coupled with a triple-quadrupole tandem API 4000 mass spectrometer (AB/MDS-Sciex, Concord, Ontario, Canada) were used for analysis. Instrumental control and data processing was performed by the Analyst software (version 1.4.2). The LC separation was performed on an Agilent Zorbax SB-C18column (150 mm × 2.1 mm I.D., 3.5 μm, Agilent Technologies, Wilmington, DE, USA) with a security guard column (12.5 mm × 2.1 mm I.D., 5 μm, Agilent Zorbax SB-C18, DE, USA). The mobile phase consisted of methanol and deionized water (v:v, 60:40) containing 0.1% formic acid at a flow rate of 0.30 mL/min, the autosampler temperature was set at 15 °C and the column temperature was set at 30 °C. A MS detector with electrospray ionization (ESI) interface was used in our experiment and the positive ion mode was selected for quantitative analysis. Quantitation was detected by the multi-reaction-monitoring (MRM) mode of transitions of m/z 316.2 → 298.2 for Mitiglinide, m/z 318.2 → 120.2 for Nateglinide. The optimized conditions used for the ESI+ source were set as follows: capillary voltage 5.5 KV; turbo heater temperature 600 °C; curtain gas (CUR) 40 psi; collision activation dissociation (CAD) 8 psi; declustering potential (DP) 81 V; collision energy (CE) 23 eV for mitiglinide and 29 eV for IS respectively. Blood glucose concentration was examined by the glucose oxidase method with the use of an autoanalyzer (TosohG8, Japan).
Graphical and all other statistical analyses, including the evaluation outputs, were performed within the Phoenix platform .
Pharmacokinetic/pharmacodynamic model and data analysis
A two-compartment linked with Emax PK-PD model was implemented by maximizing the log-likelihood using the first-order conditional estimation (FOCE) method of the Phoenix NLME software (Pharsight Corp, St. Louis, Missouri, United States of America). Pharmacokinetics and pharmacodynamics values of mitiglinide were modelled sequentially. The models were established as a series of equations, describing the relationship between plasma medicine and glucose. The equations were solved numerically and fitted to the data with Phoenix NLME . Plasma concentration of mitiglinide was used for the PK analysis. The two-compartment model for oral administration was established by the least-squares method, and the pharmacokinetic parameters including V, V2, CL and CL2 were also calculated. The amounts of mitiglinide and plasma concentrations in the central compartment (C) were determined and calculated by the distribution volume for the central compartment (V). The two-compartment model as shown in Fig. 1, and the overall PK model could be described as Eqs. (1) to (11) as below:
A means the amount of mitiglinide, A1 and A2 suggest the amounts of drug in the central and peripheral compartments respectively. The concentrations of mitiglinide in the central and peripheral compartment are represented as C and C2 respectively. The elimination clearances of the drug in the central and peripheral compartments are represented as CL and CL2 respectively. TV indicates the typical value of population mean, and η indicates the inter-individual variation. PK parameters were obtained by simultaneously fitting the plasma concentration data after oral administration of mitiglinide (10 mg) to volunteers using the 2-compartment PK model (Table 1).
An Emax PD model was selected based on previous research on describing the relationship between mitiglinide concentrations and the glucose levels in the plasma [14,15,16]. The Biophase model which consists of compartmental PK model in conjunction with Emax PD model is one of the commonly used PK-PD models. If biophase model could adequately fit both the PK and PD profiles, it would render it possible for further application. The glucose plasma concentration (mg/mL) was the response variables in the experiment, and the model can be described by equations as shown below:
Where E is the level of glucose in the plasma, E0 is the blood glucose maintained at a certain level in absence of the drug. C is the plasma concentration in the central compartment and Ce represents the plasma concentration in the effect compartment. Keo represents the glucose disappearance rate constant. TV represents the typical value of the population mean and η means the inter-individual variation, and the inter-individual variation for parameters with high shrinkage value (>0.5) were not included in the model.
The plasma concentration vs time curve following oral administration of 10 mg of mitiglinide was shown in Fig. 2. Solid line represents the mean concentration of the test coupled with the standard deviation on each time point. Terminal elimination half-life (t 1/2) is 1.69 ± 0.16 h and the CL/F 7.80 ± 1.84 L/h. Figure 2 also shows the glucose profiles after drug administration, the model independent parameters are calculated by using individual approach with Non compartmental analysis, and the model dependent parameters are obtained by using population approach with clearance parameterization. The plasma glucose levels declined after a lag of 0.2 h, and hit its bottom values of 2.6 mg/L by 0.5 h after administration. Mitiglinide declined from 0.2 to 2.2 h after administration. PK parameters were obtained by simultaneously fitting the plasma concentration data after oral administration of mitiglinide to volunteers using the 2-compartment model.
The two-compartment model analyzes the results based on the maximum likelihood criterion and visual inspection of the fittings. The data was calculated by the first-order absorption with a lag time and linear elimination. Non-compartmental and population modeling parameters are listed in Tables 1 and 2. PK/PD was analyzed based on the pharmacokinetic parameters obtained from the test as an input function. Figure 3 illustrates the fitting result of the plasma concentration and plasma glucose levels with population PK-PD model where the solid line represents the best fit of the model.
The diagnosis chart of final model was shown in Fig. 4. A is the scatter plots which represents the conditional weighted residuals (CWRES) and independent variable (IVAR). The CWRES distribution of the final model is narrow, and uniform on both sides of the zero line, which indicates that the model fitted well. B is the final model CWRES on the PRED of the scatter diagram. The value of CWRES ranged between 2 and −2 indicating the minor differences in prediction model. C is the observed value (DV) of individual predicted value (IPRE) scatter diagram and D is the observed concentrations (DV) versus population predicted values (PRED).
Based the results above, the two-compartment PK model of first-order absorption with a lag time and linear elimination conjunction with biophase PD model was developed which could capture both PK and PD profiles of mitiglinide in the present reasonably well.
Mitiglinide could reduce the early-phase insulin release and promote the controlling of the postprandial glucose levels with its ability to fasten insulin secretagogue . The characteristic of mitiglinide of anti-hyperglycemic can take effect more rapidly compared with sulphonylureas. The fast release of insulin within a few minutes was due to the rapid absorption rate of mitiglinide after oral administration . In the 90% confidence internal of geometric mean ratios, mitiglinide reached its peak plasma concentration (C max) and AUClast, which represented the area under the plasma concentration-time curve from time zero to the time of the last measurable concentration, were observed as 0.9694 (0.8120, 1.1573) and 0.8951 (0.8440, 0.9494), respectively . Mitiglinide could bind to the SUR and activate it in the pancreatic β-cells, which was considered as the mechanism of the drug underlying the therapeutic effect [20, 21]. In the present study, mitiglinide was given by oral administration and the short elimination half-life, t 1/2 (h) was 1.69. Therefore, the fast absorption in the first 2 h may be caused by the binding of the free form drug to the SUR. The total clearance (CL/F) was 7.79 mL/h in our study, which is comparable to rat (0.83 mL/h) in previous studies .
PK-PD modeling of glucose lowering effects
The characters of PK behaviors and its influences to the values of plasma glucose levels have been studied in several hypoglycemic drugs [23, 24]. Toshiyuki Takanohashi et al. reported the pk-pd. modeling methodology to describe the rapidity therapeutic effect of mitiglinide by using the receptor-binding-dissociation in rats . However, model of inhaled glucose in healthy humans has not been published previously. Phoenix platform is the commonly used software to investigate PK/PD characters depending on the 2-compartment model [25, 26]. Our study analyzed the effects of mitiglinide on glucose lowering after a single oral dose. The shorter duration action of mitiglinide and the decreasing glucose levels in plasma was described by a two-compartment model followed by the least-squares method. The fitting results were adequate, and the efficacy of the drug has been appropriated proposed by the empirical biexponential equation. Furthermore, the model is flexible for the concentration and efficiency of the drug. The T max (h) of mitiglinide was 0.25 (0.17–0.50) in our study and the glucose decrease reached its maximum.
Although mitiglinide is a commonly used and effective drug for diabetes mellitus in clinic, its side effects and the interaction remain to be studied further. A recent investigation shows that mitiglinide treatment could significantly decrease plasma fibroblast growth factor levels, which would have various effects on pancreatic islets . Our results may illustrate the physiological effects of mitiglinide and its relationship with plasma concentration. Different kinds of indirect response models can be used to explain the pk-pd. property of drug. In our study, we investigated glucose disappearance and mitiglinide concentration levels in plasma by building the pk-pd. modeling. The CWRES distribution of the final model uniformed on both sides of the zero line, and the value ranged between 2 and −2 indicating the minor differences in prediction model. The result of our research could describe the drug’s pharmacological action. The results of our experiment are benefit for study designs and appropriate dosage selection of diabetes mellitus in clinical trials.
A PK/PD model illustrates that the relationship between mitiglinide concentration in plasma and glucose lowering effect in healthy Chinese volunteers was established. The PK data of mitiglinide after oral administration were analyzed by a two-compartment model of first-order absorption with a lag time and linear elimination. The results of our experiment suggested that the model can be used reasonably to predict the relationship between PK and PD in mitiglinide. The established PK-PD model may be used to describe the drug’s pharmacological action. The results of our experiment are flexible for the concentration and efficiency of the drug and may be used for study designs and diabetes mellitus dosage control in clinical trials.
Franz MJ, Boucher JL, Rutten-Ramos S, VanWormer JJ. Lifestyle weight-loss intervention outcomes in overweight and obese adults with type 2 diabetes: a systematic review and meta-analysis of randomized clinical trials. J Acad Nutr Diet. 2015;115:1447–63.
Tahara A, Matsuyama-Yokono A, Shibasaki M. Effects of antidiabetic drugs in high-fat diet and streptozotocin–nicotinamide-induced type 2 diabetic mice. Eur J Pharmacol. 2011;655:108–16.
Guerrero-Berroa E, Schmeidler J, Beeri MS. Neuropathology of type 2 diabetes: a short review on insulin-related mechanisms. Eur Neuropsychopharmacol. 2014;24:1961–6.
Ashcroft FM, Gribble FM. ATP-sensitive K+ channels in health and disease. Diabetologia. 1999;42:903–19.
Lim JG, Lee HY, Yun JE, Kim SP, Park JW, Suh SI, Jang BC, Cho CH, Bae JH, Kim SS, Han J, Park MJ, Song DK. Taurine block of cloned ATP-sensitive Kt channels with different sulfonylurea receptor subunits expressed in Xenopus laevis oocytes. Biochem Pharmacol. 2004;68:901–10.
Ohnota H, Koizumi T, Tsutsumi N, Kobayashi M, Inoue S, Sato F, Novel rapidand short-acting hypoglycemic agent. a calcium (2s)-2-benzyl-3-(cishexahydro-2-isoindolinylcarbonyl) propionate (KAD- 1229) that acts on the sulfonylurea receptor: comparison of effects between KAD-1229 and gliclazide. J Pharmacol Exp Ther. 1994;269:489–95.
Sunaga Y, Gonoi T, Shibasaki T, Ichikawa K, Kusama H, Yano H, Seino S. The effects of mitiglinide (KAD-1229), a new anti-diabetic drug, on ATP-sensitive Kq channels and insulin secretion: comparison with the sulfonylureas and nateglinide. Eur J Pharmacol. 2001;431:119–25.
Assaloni R, Da Ros R, Quagliaro L, Piconi L, Maier A, Zuodar G, Motz E, Ceriello A. Effects of S21403 (mitiglinide) on postprandial generation of oxidative stress and inflammation in type 2 diabetic patients. Diabetologia. 2005;48:1919–24.
Lushan Y, Su Z. Determination of mitiglinide in rat plasma by high-performance liquid chromatography with UV detection. J Chromatogr B AnalytTechnol Biomed Life Sci. 2006;834:204–7.
Zhang J, Cai LJ, Peng WX, Zhu RH, Yang J, Zhang QZ. The effects of food on the pharmacokinetics of mitiglinide tablets in healthy volunteers and a novel mass-spectrometric (UPLC-MS/MS) method for such studies. J Clin Pharm Ther. 2012;37:95–9.
Lemenuel-Diot A, Laveille C, Frey N, Jochemsen R, Mallet A. Mixture modeling for the detection of subpopulations in a pharmacokinetic pharmacodynamic analysis. J Pharmacokinet Pharmacodyn. 2007;34:157–81.
Takanohashi T, Arisaka H, Ubukata K, Hayashi M, Yamada Y. Comparison of the rapidity of onset of the therapeutic effect between mitiglinide by PK/PD analysis in rats. Eur J Drug Metab Pharmacokinet. 2012;37:9–15.
Vu T, Ma P, Xiao JJ, Wang YM, Malik FI, Chow AT. Population Pharmacokinetic-Pharmacodynamic Modeling of Omecamtiv Mecarbil, a Cardiac Myosin Activator, in Healthy Volunteers and Patients With Stable Heart Failure. J Clin Pharmacol. 2015;55:1236–47.
Yun HY, Park HC, Kang W, Kwon KI. Pharmacokinetic and pharmacodynamic modeling of the effects of mitiglinide on insulin secretion and glucose lowering in healthy humans. J Clin Pharm Ther. 2006;31:469–76.
Reddy MB, Yang KH, Rao G, Rayner CR, Nie J, Pamulapati C, Marathe BM, Forrest A, Govorkova EA. Oseltamivir Population Pharmacokinetics in the Ferret: Model Application for Pharmacokinetic/Pharmacodynamic Study Design. PLoS One. 2015;10:e0138069.
Liu D, Lon HK, Dubois DC, Almon RR, Jusko WJ. Population pharmacokinetic–pharmacodynamic–disease progression model for effects of anakinra in Lewis rats with collagen-induced arthritis. J Pharmacokinet Pharmacodyn. 2011;38:769–86.
Guardado-Mendoza R, Prioletta A, Jimenez-Ceja LM, Sosale A, Folli F. The role of nateglinide and repaglinide, derivatives of meglitinide, in the treatment of Type 2 diabetes mellitus. Arch Med Sci. 2013;9:936–43.
Dornhorst A. Insulinotropic meglitinide analogues. Lancet. 2001;358:1709–16.
Jung JA, Kim JR, Kim SR, Kim TE, Lee SY, Ko JW, Huh W. Pharmacokinetics of a fixed-dose combination of mitiglinide and metformin versus concurrent administration of individual formulations in healthy subjects: a randomized, open-label, two-treatment, two-period, two-sequence, single-dose, crossover study. Clin Drug Investig. 2012;32:799–804.
Maruyama I, Tomiyama Y, Maruyama K, Ojima K, Kobayashi K, Kobayashi M, Yamazaki Y, Kojima M, Shibata N. Effects of mitiglinide and sulfonylureas in isolated canine coronary arteries and perfused rat hearts. Eur J Pharmacol. 2001;531:194–200.
Blicklé JF. Meglitinide analogues: a review of clinical data focused on recent trials. Diabetes Metab. 2006;32:113–20.
Zeng Y, Xie X, Duan J, Zhou T, Zhang Y, Yang M, Xu F. Perturbation of mitiglinide metabolism by chronic unpredicted mild stress in rats. Sci Rep. 2014;4:3794.
Li H, Xu J, Fan X. Target-mediated pharmacokinetic/pharmacodynamic model based meta-analysis and dosing regimen optimization of a long-acting release formulation of exenatide in patients with type 2 diabetes mellitus. J PharmacolSci. 2015;127:170–80.
Landersdorfer CB, Jusko WJ. Pharmcokinetic pharmacodynamics modeling of glucose clamp effects of inhaled and subcutaneous insulin in healthy volunteers and diabetic patients. Drug Metab Pharmacokinet. 2010;25:418–29.
Kobuchi S, Ito Y, Hayakawa T, Nishimura A, Shibata N, Takada K, Sakaeda T. Pharmacokinetic-pharmacodynamic (PK-PD) modeling and simulation of 5 fluorouracil for erythropenia in rats. J Pharmacol Toxicol Methods. 2014;70:134–44.
Wang J, Xia S, Xue W, Wang D, Sai Y, Liu L, Liu X. A semi physiologically based pharmacokinetic model characterizing mechanism-based auto-inhibition to predictstereoselective pharmacokinetics of verapamil and its metabolite norverapamil in human. Eur J Pharm Sci. 2013;50:290–302.
Wang B, Yang G, Yang M, Liao Y, Li Z, Boden G, Li L. Mitiglinide treatment may decreases plasma fibroblast growth factor-21 levels in individuals with new-onset T2DM. Cytokine. 2012;57:300–3.
Availability of data and materials
All the data and material of the project are in availability and saved at Affiliated Hospital of Nanjing University of Chinese Medicine.
The authors are grateful for the financial support provided by the National Natural Science Foundation of China (No. 81403174 and 81503300), the Natural Science Foundation of Jiangsu province (No. BK20131035, BK20161610 and BK20151602), Jiangsu Provincial Medical Youth Talent (No. QNRC2016642) and the top six talent project of Jiangsu province 2016(WSN-051).
Ethics approval and consent to participate
The study was approved by the Ethics Committee of the Affiliated Hospital of Nanjing University of Chinese Medicine (IRB#2014NL-020-02). All subjects were informed in oral and written form that they could withdraw from the trial any time for any reason. Consent was obtained in written form before any trial-related activities, and the investigator retained the consent forms.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.