PBPK modeling to predict the pharmacokinetics of venlafaxine and its active metabolite in different CYP2D6 genotypes and drug-drug interactions with clarithromycin and paroxetine.

Venlafaxine, a serotonin-norepinephrine reuptake inhibitor (SNRI), is indicated for the treatment of major depressive disorder, social anxiety disorder, generalized anxiety disorder, and panic disorder. Venlafaxine is metabolized to the active metabolite desvenlafaxine mainly by CYP2D6. Genetic polymorphism of CYP2D6 and coadministration with other medications can significantly affect the pharmacokinetics and/or pharmacodynamics of venlafaxine and its active metabolite. This study aimed to establish the PBPK models of venlafaxine and its active metabolite related to CYP2D6 genetic polymorphism and to predict drug-drug interactions (DDIs) with clarithromycin and paroxetine in different CYP2D6 genotypes. Clinical pharmacogenomic data for venlafaxine and desvenlafaxine were collected to build the PBPK model. Physicochemical and absorption, distribution, metabolism, and excretion (ADME) characteristics of respective compounds were obtained from previously reported data, predicted by the PK-Sim® software, or optimized to capture the plasma concentration-time profiles. Model evaluation was performed by comparing the predicted pharmacokinetic parameters and plasma concentration-time profiles to the observed data. Predicted plasma concentration-time profiles of venlafaxine and its active metabolite were visually similar to the observed profiles and all predicted AUC and Cmax values for respective compounds were included in the twofold error range of observed values in non-genotyped populations and different CYP2D6 genotypes. When clarithromycin or clarithromycin plus paroxetine was concomitantly administered, predicted plasma concentration-time profiles of venlafaxine properly captured the observed profiles in two different CYP2D6 genotypes and all predicted DDI ratios for AUC and Cmax were included within the acceptance range. Consequently, the present model successfully captured the pharmacokinetic alterations of venlafaxine and its active metabolite according to CYP2D6 genetic polymorphism as well as the DDIs between venlafaxine and two CYP inhibitors. The present model can be used to predict the pharmacokinetics of venlafaxine and its active metabolite considering different races, ages, coadministered drugs, and CYP2D6 activity of individuals and it can contribute to individualized pharmacotherapy of venlafaxine.

PBPK modeling to predict the pharmacokinetics of pantoprazole in different CYP2C19 genotypes.

Pantoprazole is used to treat gastroesophageal reflux disease (GERD), maintain healing of erosive esophagitis (EE), and control symptoms related to Zollinger-Ellison syndrome (ZES). Pantoprazole is mainly metabolized by cytochrome P450 (CYP) 2C19, converting to 4'-demethyl pantoprazole. CYP2C19 is a genetically polymorphic enzyme, and the genetic polymorphism affects the pharmacokinetics and/or pharmacodynamics of pantoprazole. In this study, we aimed to establish the physiologically based pharmacokinetic (PBPK) model to predict the pharmacokinetics of pantoprazole in populations with various CYP2C19 metabolic activities. A comprehensive investigation of previous reports and drug databases was conducted to collect the clinical pharmacogenomic data, physicochemical data, and disposition properties of pantoprazole, and the collected data were used for model establishment. The model was evaluated by comparing the predicted plasma concentration-time profiles and/or pharmacokinetic parameters (AUC and Cmax) with the clinical observation results. The predicted plasma concentration-time profiles in different CYP2C19 phenotypes properly captured the observed profiles. All fold error values for AUC and Cmax were included in the two-fold range. Consequently, the minimal PBPK model for pantoprazole related to CYP2C19 genetic polymorphism was properly established and it can predict the pharmacokinetics of pantoprazole in different CYP2C19 phenotypes. The present model can broaden the insight into the individualized pharmacotherapy for pantoprazole.

Physiologically based pharmacokinetic (PBPK) modeling of pitavastatin in relation to SLCO1B1 genetic polymorphism.

Pitavastatin, a potent 3-hydroxymethylglutaryl coenzyme A reductase inhibitor, is indicated for the treatment of hypercholesterolemia and mixed dyslipidemia. Hepatic uptake of pitavastatin is predominantly occupied by the organic anion transporting polypeptide 1B1 (OATP1B1) and solute carrier organic anion transporter family member 1B1 (SLCO1B1) gene, which is a polymorphic gene that encodes OATP1B1. SLCO1B1 genetic polymorphism significantly alters the pharmacokinetics of pitavastatin. This study aimed to establish the physiologically based pharmacokinetic (PBPK) model to predict pitavastatin pharmacokinetics according to SLCO1B1 genetic polymorphism. PK-Sim® version 10.0 was used to establish the whole-body PBPK model of pitavastatin. Our pharmacogenomic data and a total of 27 clinical pharmacokinetic data with different dose administration and demographic properties were used to develop and validate the model, respectively. Physicochemical properties and disposition characteristics of pitavastatin were acquired from previously reported data or optimized to capture the plasma concentration-time profiles in different SLCO1B1 diplotypes. Model evaluation was performed by comparing the predicted pharmacokinetic parameters and profiles to the observed data. Predicted plasma concentration-time profiles were visually similar to the observed profiles in the non-genotyped populations and different SLCO1B1 diplotypes. All fold error values for AUC and Cmax were included in the two fold range of observed values. Thus, the PBPK model of pitavastatin in different SLCO1B1 diplotypes was properly established. The present study can be useful to individualize the dose administration strategy of pitavastatin in individuals with various ages, races, and SLCO1B1 diplotypes.

Physiologically based pharmacokinetic (PBPK) modeling to predict the pharmacokinetics of irbesartan in different CYP2C9 genotypes.

Irbesartan, a potent and selective angiotensin II type-1 (AT1) receptor blocker (ARB), is one of the representative medications for the treatment of hypertension. Cytochrome P450 (CYP) 2C9 is primarily involved in the oxidation of irbesartan. CYP2C9 is highly polymorphic, and genetic polymorphism of this enzyme is the leading cause of significant alterations in the pharmacokinetics of irbesartan. This study aimed to establish the physiologically based pharmacokinetic (PBPK) model to predict the pharmacokinetics of irbesartan in different CYP2C9 genotypes. The irbesartan PBPK model was established using the PK-Sim® software. Our previously reported pharmacogenomic data for irbesartan was leveraged in the development of the PBPK model and collected clinical pharmacokinetic data for irbesartan was used for the validation of the model. Physicochemical and ADME properties of irbesartan were obtained from previously reported data, predicted by the modeling software, or optimized to fit the observed plasma concentration-time profiles. Model evaluation was performed by comparing the predicted plasma concentration-time profiles and pharmacokinetic parameters to the observed results. Predicted plasma concentration-time profiles were visually similar to observed profiles. Predicted AUCinf in CYP2C9*1/*3 and CYP2C9*1/*13 genotypes were increased by 1.54- and 1.62-fold compared to CYP2C9*1/*1 genotype, respectively. All fold error values for AUC and Cmax in non-genotyped and CYP2C9 genotyped models were within the two-fold error criterion. We properly established the PBPK model of irbesartan in different CYP2C9 genotypes. It can be used to predict the pharmacokinetics of irbesartan for personalized pharmacotherapy in individuals of various races, ages, and CYP2C9 genotypes.

Effects of CYP2D6 and CYP2C19 genetic polymorphisms and cigarette smoking on the pharmacokinetics of tolperisone.

Tolperisone, a muscle relaxant used for post-stroke spasticity, is metabolized to its main metabolite by CYP2D6 and to a lesser extent by CYP2C19 and CYP1A2. We investigated the effects of CYP2D6 and CYP2C19 genetic polymorphisms and cigarette smoking on tolperisone pharmacokinetics. A 150 mg oral dose of tolperisone was given to 184 healthy Korean subjects and plasma concentrations of tolperisone were measured by liquid chromatography-tandem mass spectrometry (LC-MS/MS). A 3.14-fold significant increase in AUC0-∞ was observed in the CYP2D6*10/*10 group compared with the CYP2D6*wt/*wt group, whereas a 3.59-fold increase in AUC0-∞ was observed in CYP2C19PMs compared to CYP2C19EMs. Smokers had a 38.5% decrease in AUC0-∞ when compared to non-smokers. When these effects were combined, CYP2D6*10/*10-CYP2C19PM-Non-smokers had a 25.9-fold increase in AUC0-∞ compared to CYP2D6*wt/*wt-CYP2C19EM-Smokers. Genetic polymorphisms of CYP2D6 and CYP2C19 and cigarette smoking independently and significantly affected tolperisone pharmacokinetics and these effects combined resulted in a much greater impact on tolperisone pharmacokinetics.

Effects of CYP2C9 and CYP2C19 genetic polymorphisms on the pharmacokinetics and pharmacodynamics of gliclazide in healthy subjects.

Gliclazide metabolism is mediated by genetically polymorphic CYP2C9 and CYP2C19 enzymes. We investigated the effects of CYP2C9 and CYP2C19 genetic polymorphisms on the pharmacokinetics and pharmacodynamics of gliclazide. Twenty-seven Korean healthy volunteers were administered a single oral dose of gliclazide 80 mg. The plasma concentration of gliclazide was quantified for the pharmacokinetic analysis and plasma concentrations of glucose and insulin were measured as pharmacodynamic parameters. The pharmacokinetics of gliclazide showed a significant difference according to the number of defective alleles of combined CYP2C9 and CYP2C19. The two defective alleles group (group 3) and one defective allele group (group 2) showed 2.34- and 1.46-fold higher AUC0-∞ (P < 0.001), and 57.1 and 32.3% lower CL/F (P < 0.001), compared to those of the no defective allele group (group 1), respectively. The CYP2C9IM-CYP2C19IM group had AUC0-∞ increase of 1.49-fold (P < 0.05) and CL/F decrease by 29.9% (P < 0.01), compared with the CYP2C9 Normal Metabolizer (CYP2C9NM)-CYP2C19IM group. The CYP2C9NM-CYP2C19PM group and CYP2C9NM-CYP2C19IM group showed 2.41- and 1.51-fold higher AUC0-∞ (P < 0.001), and 59.6 and 35.4% lower CL/F (P < 0.001), compared to those of the CYP2C9NM-CYP2C19NM group, respectively. The results represented that CYP2C9 and CYP2C19 genetic polymorphisms significantly affected the pharmacokinetics of gliclazide. Although the genetic polymorphism of CYP2C19 had a greater effect on the pharmacokinetics of gliclazide, the genetic polymorphism of CYP2C9 also had a significant effect. On the other hand, plasma glucose and insulin responses to gliclazide were not significantly affected by the CYP2C9-CYP2C19 genotypes, requiring further well-controlled studies with long-term dosing of gliclazide in diabetic patients.

Effects of CYP2D6*10 allele on the pharmacokinetics of tolperisone.

Tolperisone, a muscle relaxant used for post-stroke spasticity, has been reported to have a very wide interindividual pharmacokinetic variability. It is metabolized mainly by CYP2D6 and, to a lesser extent, by CYP2C19 and CYP1A2. CYP2D6 is a highly polymorphic enzyme, and CYP2D6*wt/*wt, CYP2D6*wt/*10 and CYP2D6*10/*10 genotypes constitute more than 90% of the CYP2D6 genotypes in the Korean population. Thus, effects of the CYP2D6*10 on tolperisone pharmacokinetics were investigated in this study to elucidate the reasons for the wide interindividual variability. Oral tolperisone 150 mg was given to sixty-four healthy Koreans, and plasma concentrations of tolperisone were measured by liquid chromatography-tandem mass spectrometry (LC-MS/MS). The CYP2D6*10/*10 and CYP2D6*wt/*10 groups had significantly higher Cmax and lower CL/F values than the CYP2D6*wt/*wt group. The AUCinf of CYP2D6*10/*10 and CYP2D6*wt/*10 groups were 5.18-fold and 2.25-fold higher than the CYP2D6*wt/*wt group, respectively. There were considerable variations in the Cmax and AUC values within each genotype group, and the variations were greater as the activity of CYP2D6 decreased. These results suggest that the genetic polymorphism of CYP2D6 significantly affected tolperisone pharmacokinetics and factor(s) other than CYP2D6 may also have significant effects on the pharmacokinetics of tolperisone.

Effects of CYP2C19 genetic polymorphism on the pharmacokinetics of tolperisone in healthy subjects.

Tolperisone hydrochloride is a centrally-acting muscle relaxant used for relieving spasticities of neurological origin and muscle spasms associated with painful locomotor diseases. It is metabolized to the inactive metabolite mainly by CYP2D6 and, to a lesser extent, by CYP2C19 and CYP1A2. In our previous study, the pharmacokinetics of tolperisone was significantly affected by the genetic polymorphism of CYP2D6, but the wide interindividual variation of tolperisone pharmacokinetics was not explained by genetic polymorphism of CYP2D6 alone. Thus, we studied the effects of CYP2C19 genetic polymorphism on tolperisone pharmacokinetics. Eighty-one subjects with different CYP2C19 genotypes received a single oral dose of 150 mg tolperisone with 240 mL of water, and blood samples were collected up to 12 h after dosing. The plasma concentration of tolperisone was measured by a liquid chromatography-tandem mass spectrometry system. The CYP2C19PM group had significantly higher Cmax and lower CL/F values than the CYP2C19EM and CYP2C19IM groups. The AUCinf of the CYP2C19PM group was 2.86-fold and 3.00-fold higher than the CYP2C19EM and CYP2C19IM groups, respectively. In conclusion, the genetic polymorphism of CYP2C19 significantly affected tolperisone pharmacokinetics.

Physiologically based pharmacokinetic (PBPK) modeling of flurbiprofen in different CYP2C9 genotypes.

The aim of this study was to establish the physiologically based pharmacokinetic (PBPK) model of flurbiprofen related to CYP2C9 genetic polymorphism and describe the pharmacokinetics of flurbiprofen in different CYP2C9 genotypes. PK-Sim® software was used for the model development and validation. A total of 16 clinical pharmacokinetic data for flurbiprofen in different CYP2C9 genotypes, dose regimens, and age groups were used for the PBPK modeling. Turnover number (kcat) of CYP2C9 values were optimized to capture the observed profiles in different CYP2C9 genotypes. In the simulation, predicted fraction metabolized by CYP2C9, fraction excreted to urine, bioavailability, and volume of distribution were similar to previously reported values. Predicted plasma concentration-time profiles in different CYP2C9 genotypes were visually similar to the observed profiles. Predicted AUCinf in CYP2C9*1/*2, CYP2C9*1/*3, and CYP2C9*3/*3 genotypes were 1.44-, 2.05-, and 3.67-fold higher than the CYP2C9*1/*1 genotype. The ranges of fold errors for AUCinf, Cmax, and t1/2 were 0.84-1.00, 0.61-1.22, and 0.74-0.94 in development and 0.59-0.98, 0.52-0.97, and 0.61-1.52 in validation, respectively, which were within the acceptance criterion. Thus, the PBPK model was successfully established and described the pharmacokinetics of flurbiprofen in different CYP2C9 genotypes, dose regimens, and age groups. The present model could guide the decision-making of tailored drug administration strategy by predicting the pharmacokinetics of flurbiprofen in various clinical scenarios.

Physiologically based pharmacokinetic modelling to predict the pharmacokinetics of metoprolol in different CYP2D6 genotypes.

Metoprolol, a selective ß1-adrenoreceptor blocking agent used in the treatment of hypertension, angina, and heart failure, is primarily metabolized by the CYP2D6 enzyme, which catalyzes α-hydroxylation and O-desmethylation. As CYP2D6 is genetically highly polymorphic and the enzymatic activity differs greatly depending on the presence of the mutant allele(s), the pharmacokinetic profile of metoprolol is highly variable depending on the genotype of CYP2D6. The aim of study was to develop the physiologically based pharmacokinetic (PBPK) model of metoprolol related to CYP2D6 genetic polymorphism for personalized therapy with metoprolol. For PBPK modelling, our previous pharmacogenomic data were used. To obtain kinetic parameters (Km, Vmax, and CLint) of each genotype, the recombinant CYP enzyme of each genotype was incubated with metoprolol and metabolic rates were assayed. Based on these data, the PBPK model of metoprolol was developed and validated in different CYP2D6 genotypes using PK-Sim® software. As a result, the input values for various parameters for the PBPK model were presented and the PBPK model successfully described the pharmacokinetics of metoprolol in each genotype group. The simulated values were within the acceptance criterion (99.998% confidence intervals) compared with observed values. The PBPK model developed in this study can be used for personalized pharmacotherapy with metoprolol in individuals of various races, ages, and CYP2D6 genotypes.

Physiologically based pharmacokinetic (PBPK) modeling of piroxicam with regard to CYP2C9 genetic polymorphism.

Piroxicam is a non-steroidal anti-inflammatory drug used to alleviate symptoms of osteoarthritis and rheumatoid arthritis. CYP2C9 genetic polymorphism significantly influences the pharmacokinetics of piroxicam. The objective of this study was to develop and validate the piroxicam physiologically based pharmacokinetic (PBPK) model related to CYP2C9 genetic polymorphism. PK-Sim® version 10.0 was used for the PBPK modeling. The PBPK model was evaluated by predicted and observed plasma concentration-time profiles, fold errors of predicted to observed pharmacokinetic parameters, and a goodness-of-fit plot. The turnover number (kcat) of CYP2C9 was adjusted to capture the pharmacokinetics of piroxicam in different CYP2C9 genotypes. The population PBPK model overall accurately described and predicted the plasma concentration-time profiles in different CYP2C9 genotypes. In our simulations, predicted AUCinf in CYP2C9*1/*2, CYP2C9*1/*3, and CYP2C9*3/*3 genotypes were 1.83-, 2.07-, and 6.43-fold higher than CYP2C9*1/*1 genotype, respectively. All fold error values for AUC, Cmax, and t1/2 were included in the acceptance criterion with the ranges of 0.57-1.59, 0.63-1.39, and 0.65-1.51, respectively. The range of fold error values for predicted versus observed plasma concentrations was 0.11-3.13. 93.9% of fold error values were within the two-fold range. Average fold error, absolute average fold error, and root mean square error were 0.93, 1.27, and 0.72, respectively. Our model accurately captured the pharmacokinetic alterations of piroxicam according to CYP2C9 genetic polymorphism.

Effects of CYP2C9*3 and *13 alleles on the pharmacokinetics and pharmacodynamics of glipizide in healthy Korean subjects.

Glipizide is a second-generation sulfonylurea antidiabetic drug. It is principally metabolized to inactive metabolites by genetically polymorphic CYP2C9 enzyme. In this study, we investigated the effects of CYP2C9*3 and *13 variant alleles on the pharmacokinetics and pharmacodynamics of glipizide. Twenty-four healthy Korean volunteers (11 subjects with CYP2C9*1/*1, 8 subjects with CYP2C9*1/*3, and 5 subjects with CYP2C9*1/*13) were recruited for this study. They were administered a single oral dose of glipizide 5 mg. The plasma concentration of glipizide was quantified for pharmacokinetic analysis and plasma glucose and insulin concentrations were measured as pharmacodynamic parameters. The results represented that CYP2C9*3 and *13 alleles significantly affected the pharmacokinetics of glipizide. In subjects with CYP2C9*1/*3 and CYP2C9*1/*13 genotypes, the mean AUC0-∞ were increased by 44.8% and 58.2%, respectively (both P < 0.001), compared to those of subjects with CYP2C9*1/*1 genotype, while effects of glipizide on plasma glucose and insulin levels were not significantly different between CYP2C9 genotype groups. In conclusion, individuals carrying the defective CYP2C9*3 and CYP2C9*13 alleles have markedly elevated plasma concentrations of glipizide compared with CYP2C9*1/*1 wild-type.

Physiologically based pharmacokinetic (PBPK) modeling of meloxicam in different CYP2C9 genotypes.

Meloxicam, a non-steroidal anti-inflammatory drug, is used for the treatment of rheumatoid arthritis and osteoarthritis. Cytochrome P450 (CYP) 2C9 and CYP3A4 are major and minor enzymes involved in the metabolism of meloxicam. Impaired enzyme activity of CYP2C9 variants increases the plasma exposures of meloxicam and the risk of adverse events. The objective of our study is to develop and validate the physiologically based pharmacokinetic (PBPK) model of meloxicam related to CYP2C9 genetic polymorphism using the PK-Sim® software. In vitro kcat of CYP2C9 was optimized in different CYP2C9 genotypes. The demographic and pharmacokinetic dataset for the development of the PBPK model was extracted from two previous clinical pharmacokinetic studies. Thirty-one clinical datasets, representing different dose regimens and demographic characteristics, were utilized to validate the PBPK model. The shapes of simulated plasma concentration-time profiles in each CYP2C9 genotype were visually similar to observed profiles. The predicted exposures (AUCinf) of meloxicam in CYP2C9*1/*3, CYP2C9*1/*13, and CYP2C9*3/*3 genotypes were increased by 1.77-, 2.91-, and 8.35-fold compared to CYP2C9*1/*1 genotype, respectively. In all datasets for the development and validations, fold errors between predicted and observed pharmacokinetic parameters were within the two-fold error criteria. As a result, the PBPK model was appropriately established and properly described the pharmacokinetics of meloxicam in different CYP2C9 genotypes. This study is expected to contribute to reducing the risk of adverse events of meloxicam through optimization of meloxicam dosing in different CYP2C9 genotypes.

Physiologically based pharmacokinetic modeling of candesartan related to CYP2C9 genetic polymorphism in adult and pediatric patients.

Candesartan cilexetil is an angiotensin II receptor blocker and it is widely used to treat hypertension and heart failure. This drug is a prodrug that rapidly converts to candesartan after oral administration. Candesartan is metabolized by cytochrome P450 2C9 (CYP2C9) enzyme or uridine diphosphate glucurinosyltransferase 1A3, or excreted in an unchanged form through urine, biliary tract and feces. We investigated the effect of genetic polymorphism of CYP2C9 enzyme on drug pharmacokinetics using physiologically based pharmacokinetic (PBPK) modeling. In addition, by introducing the age and ethnicity into the model, we developed a model that can propose an appropriate dosage regimen taking into account the individual characteristics of each patient. To evaluate the suitability of the model, the results of a clinical trial on twenty-two healthy Korean subjects and their CYP2C9 genetic polymorphism data was applied. In this study, PK-Sim® was used to develop the PBPK model of candesartan.

Physiologically based pharmacokinetic (PBPK) modelling of tamsulosin related to CYP2D6*10 allele.

Tamsulosin, a selective [Formula: see text]-adrenoceptor blocker, is commonly used for alleviation of lower urinary tract symptoms related to benign prostatic hyperplasia. Tamsulosin is predominantly metabolized by CYP3A4 and CYP2D6 enzymes, and several studies reported the effects of CYP2D6 genetic polymorphism on the pharmacokinetics of tamsulosin. This study aims to develop and validate the physiologically based pharmacokinetic (PBPK) model of tamsulosin in CYP2D6*wt/*wt, CYP2D6*wt/*10, and CYP2D6*10/*10 genotypes, using Simcyp® simulator. Physicochemical, and formulation properties and data for absorption, distribution, metabolism and excretion were collected from previous publications, predicted in the simulator, or optimized in different CYP2D6 genotypes. The tamsulosin PBPK model in CYP2D6*wt/*wt and CYP2D6*wt/*10 genotypes were developed based on the clinical pharmacokinetic study where a single oral dose of 0.2 mg tamsulosin was administered to 25 healthy Korean male volunteers with CYP2D6*wt/*wt and CYP2D6*wt/*10 genotypes. A previous pharmacokinetic study was used to develop the model in CYP2D6*10/*10 genotype. The developed model was validated using other clinical pharmacokinetic studies not used in development. The predicted exposures via the PBPK model in CYP2D6*wt/*10 and CYP2D6*10/*10 genotype was 1.23- and 1.76-fold higher than CYP2D6*wt/*wt genotype, respectively. The simulation profiles were visually similar to the observed profiles, and fold errors of all development and validation datasets were included within the criteria. Therefore, the tamsulosin PBPK model in different CYP2D6 genotypes with regards to CYP2D6*10 alleles was appropriately established. Our model can contribute to the implementation of personalized pharmacotherapy of patients, appropriately predicting the pharmacokinetics of tamsulosin reflecting their demographic and CYP2D6 genotype characteristics without unnecessary drug exposure.

Physiologically based pharmacokinetic (PBPK) modeling for prediction of celecoxib pharmacokinetics according to CYP2C9 genetic polymorphism.

Celecoxib is a non-steroidal anti-inflammatory drug (NSAID) and a representative selective cyclooxygenase (COX)-2 inhibitor, which is commonly prescribed for osteoarthritis, rheumatoid arthritis, ankylosing spondylitis, acute pain, and primary dysmenorrhea. It is mainly metabolized by CYP2C9 and partly by CYP3A4 after oral administration. Many studies reported that CYP2C9 genetic polymorphism has significant effects on the pharmacokinetics of celecoxib and the occurrence of adverse drug reactions. The aim of this study was to develop a physiologically based pharmacokinetic (PBPK) model of celecoxib according to CYP2C9 genetic polymorphism for personalized pharmacotherapy. Initially, a clinical pharmacokinetic study was conducted where a single dose (200 mg) of celecoxib was administered to 39 healthy Korean subjects with CYP2C9*1/*1 or CYP2C9*1/*3 genotypes to obtain data for PBPK development. Based on the conducted pharmacokinetic study and a previous pharmacokinetic study involving subjects with CYP2C9*1/*13 and CYP2C9*3/*3 genotype, PBPK model for celecoxib was developed. A PBPK model for CYP2C9*1/*1 genotype group was developed and then scaled to other genotype groups (CYP2C9*1/*3, CYP2C9*1/*13 and CYP2C9*3/*3). After model development, model validation was performed with comparison of five pharmacokinetic studies. As a result, the developed PBPK model of celecoxib successfully described the pharmacokinetics of each CYP2C9 genotype group and its predicted values were within the acceptance criterion. Additionally, all the predicted values were within two-fold error range in comparison to the previous pharmacokinetic studies. This study demonstrates the possibility of determining the appropriate dosage of celecoxib for each individual through the PBPK modeling with CYP2C9 genomic information. This approach could contribute to the reduction of adverse drug reactions of celecoxib and enable precision medicine.

Near-Infrared Autofluorescence Imaging May Reduce Temporary Hypoparathyroidism in Patients Undergoing Total Thyroidectomy and Central Neck Dissection.

Background: Near-infrared autofluorescence (NIRAF) imaging is known to reduce the incidence of post-thyroidectomy hypocalcemia. However, there are no studies on how much NIRAF imaging affects the serum parathyroid hormone (PTH) level after surgery. We investigated the changes of the serum PTH level and ionized calcium (iCa.) in patients undergoing total thyroidectomy with central neck dissection (CND). Materials and Methods: This retrospective study with historical control enrolled 542 patients who underwent total thyroidectomy with CND. Patients were divided into two groups: the NIRAF group (261 patients) and the control group (281 patients). PTH and iCa. levels were measured at the hospital stay, 1, 3, and 6 months after surgery. In addition, the number of identified parathyroid glands (PGs), autotransplanted PGs, and the inadvertent resection rate of PGs was evaluated. Results: The incidence of postoperative hypoparathyroidism (PTH <15 pg/mL) was significantly lower in the NIRAF group during the hospitalization (88 patients: 33.7% vs. 131 patients: 46.6%; p = 0.002) and at 1 month postoperatively (23 patients: 8.8% vs. 53 patients: 18.9%; p = 0.001). There was no difference in the permanent hypoparathyroidism rate (6 months after surgery) between the NIRAF group and the control group (4.2% vs. 4.6%; p = 0.816). There was no difference in the incidence of hypocalcemia (iCa. <1.09 mmol/L) (during hospitalization: 6.5% vs. 10.0%; 1 month: 2.3% vs. 2.5%; 3 months: 0.8% vs. 0.7%; 6 months after surgery: 1.1% vs. 1.1%) between the two groups. The number of inadvertently resected PGs was significantly lower in the NIRAF group (18:6.9% vs. 36:12.8%; p = 0.021). Conclusions: These results suggest that NIRAF imaging may reduce temporary hypoparathyroidism and the risk of inadvertent resection of PGs in patients undergoing total thyroidectomy with CND.

ABCB1 c.2677G>T/c.3435C>T diplotype increases the early-phase oral absorption of losartan.

Losartan has been shown to be a substrate of the drug-efflux transporter MDR1, encoded by the ABCB1 gene. ABCB1 c.2677G>T and c.3435C>T variants are known to be associated with reduced expression and function of P-glycoprotein (P-gp). We investigated the effects of ABCB1 diplotype on the pharmacokinetics of losartan. Thirty-eight healthy Korean volunteers with different ABCB1 diplotypes [c.2677G> T and c.3435C>T; carriers of GG/CC (n = 13), GT/CT (n = 12) and TT/TT (n = 13) diplotype] were recruited and administered a single 50 mg oral dose of losartan potassium. Losartan and its active metabolite E-3174 samples in plasma and urine were collected up to 10 and 8 h after drug administration, respectively, and the concentrations of both samples were determined by HPLC method. Significant differences were observed in Cmax of losartan and losartan plus E-3174 (Lo + E) among the three diplotype groups (both P < 0.01). However, the power of the performed test is less than the desired power (0.800). The tmax of losartan and E-3174 in three diplotype groups were also significantly different (both P < 0.01). The AUC values of Lo + E were significantly different among the three diplotype groups until 6 h after losartan administration (P < 0.01). On the contrary, AUC at the periods of 8-10 h and 10 h-infinity of Lo + E were significantly lower in the TT/TT group than in the GG/CC group. Urinary excretion of losartan until 4 h after losartan administration in the TT/TT group was higher than that of the GG/CC group. These results suggest that c.2677G>T/c.3435C>T diplotypes of ABCB1 may significantly increase the early-phase absorption of losartan, but not the total absorption.

Effects of paroxetine on the pharmacokinetics of atomoxetine and its metabolites in different CYP2D6 genotypes.

The aim of this study was to investigate the effects of paroxetine, a potent inhibitor of CYP2D6, on the pharmacokinetics of atomoxetine and its two metabolites, 4-hydroxyatomoxetine and N-desmethylatomoxetine, in different CYP2D6 genotypes. Twenty-six healthy subjects were recruited and divided into CYP2D6*wt/*wt (*wt=*1 or *2, n = 10), CYP2D6*wt/*10 (n = 9), and CYP2D6*10/*10 groups (n = 7). In atomoxetine phase, all subjects received a single oral dose of atomoxetine (20 mg). In paroxetine phase, after administration of a single oral dose of paroxetine (20 mg) for six consecutive days, all subjects received a single oral dose of atomoxetine with paroxetine. Plasma concentrations of atomoxetine and its metabolites were determined up to 24 h after dosing. During atomoxetine phase, there were significant differences in Cmax and AUC0-24 of atomoxetine and N-desmethylatomoxetine among three genotype groups, whereas significant differences were not found in relation to CYP2D6*10 allele after administration of paroxetine. AUC ratios of 4-hydroxyatomoxetine and N-desmethylatomoxetine to atomoxetine were significantly different among three genotype groups during atomoxetine phase (all, P < 0.001), but after paroxetine treatment significant differences were not found. After paroxetine treatment, AUC0-24 of atomoxetine was increased by 2.3-, 1.7-, and 1.3-fold, in CYP2D6*wt/*wt, CYP2D6*wt/*10, and CYP2D6*10/*10 groups in comparison to atomoxetine phase, respectively. AUC ratio of 4-hydroxyatomoxetine to atomoxetine in each group was significantly decreased, whereas AUC ratio of N-desmethylatomoxetine to atomoxetine significantly increased after administration of paroxetine. In conclusion, paroxetine coadministration significantly affected pharmacokinetic parameters of atomoxetine and its two metabolites, 4-hydroxyatomoxetine and N-desmethylatomoxetine. When atomoxetine was administered alone, Cmax, AUC0-24 and CL/F of atomoxetine were significantly different among the three CYP2D6 genotype groups. However, after paroxetine coadministration, no significant differences in these pharmacokinetic parameters were observed among the CYP2D6 genotype groups.

Effects of CYP2D6 genetic polymorphism on the pharmacokinetics of metoclopramide.

Metoclopramide inhibits the central and peripheral D2 receptors and is frequently prescribed in adults and children as an antiemetic or a prokinetic drug to control symptoms of upper gastrointestinal motor disorders. Metoclopramide is predominantly metabolized via N-dealkylation and it is primarily mediated by CYP2D6 which is highly polymorphic. Thus, the effects of CYP2D6 genetic polymorphism on the pharmacokinetics of metoclopramide were evaluated in this study. All volunteers were genotyped for CYP2D6 and divided into four different genotype groups (CYP2D6*wt/*wt [*wt = *1 or *2], CYP2D6*wt/*10, CYP2D6*10/*10, and CYP2D6*5/*10). Each subject received a single oral dose of metoclopramide 10 mg. Plasma concentrations of metoclopramide were measured by using HPLC-UV. Compared to CYP2D6*wt/*wt, AUCinf of CYP2D6*wt/*10, CYP2D6*10/*10, and CYP2D6*5/*10 significantly increased by 1.5-, 2.3-, and 2.5-fold, respectively. Cmax also increased significantly in comparison to CYP2D6*wt/*wt across all genotype groups, with 1.5-, 1.7-, and 1.7-fold increases seen in CYP2D6*wt/*10, CYP2D6*10/*10, and CYP2D6*5/*10 groups, respectively. The CL/F of metoclopramide decreased in CYP2D6 genotype groups with decreased function alleles, as decreases of 37%, 56% and 61% were observed in CYP2D6*wt/10, *10/10, and *5/*10 genotype groups in comparison to the CYP2D6*wt/*wt group. In conclusion, the genetic polymorphisms of CYP2D6 significantly affected metoclopramide pharmacokinetics.