Time course of enzyme induction in humans: effect of pentobarbital on nortriptyline metabolism

Pharmacokinetics during therapeutic hypothermia in neonates: from pathophysiology to translational knowledge and physiologically-based pharmacokinetic (PBPK) modeling

INTRODUCTION

Perinatal asphyxia (PA) still causes significant morbidity and mortality. Therapeutic hypothermia (TH) is the only effective therapy for neonates with moderate to severe hypoxic-ischemic encephalopathy after PA. These neonates need additional pharmacotherapy, and both PA and TH may impact physiology and, consequently, pharmacokinetics (PK) and pharmacodynamics (PD).

AREAS COVERED

This review provides an overview of the available knowledge in PubMed (until November 2022) on the pathophysiology of neonates with PA/TH. In vivo pig models for this setting enable distinguishing the effect of PA versus TH on PK and translating this effect to human neonates. Available asphyxia pig models and methodological considerations are described. A summary of human neonatal PK of supportive pharmacotherapy to improve neurodevelopmental outcomes is provided.

EXPERT OPINION

To support drug development for this population, knowledge from clinical observations (PK data, real-world data on physiology), preclinical (in vitro and in vivo (minipig)) data, and molecular and cellular biology insights can be integrated into a predictive physiologically-based PK (PBPK) framework, as illustrated by the I-PREDICT project (Innovative physiology-based pharmacokinetic model to predict drug exposure in neonates undergoing cooling therapy). Current knowledge, challenges, and expert opinion on the future directions of this research topic are provided.

Effects of Enzyme Induction and Polymorphism on the Pharmacokinetics of Isoniazid and Rifampin in Tuberculosis/HIV Patients

Tuberculosis is the most common cause of death in HIV-infected individuals. Rifampin and isoniazid are the backbones of the current first-line antitubercular therapy. The aim of the present study was to describe the time-dependent pharmacokinetics and pharmacogenetics of rifampin and isoniazid and to quantitatively evaluate the drug-drug interaction between rifampin and isoniazid in patients coinfected with HIV. Plasma concentrations of isoniazid, acetyl-isoniazid, isonicotinic acid, rifampin, and 25-desacetylrifampin from 40 HIV therapy-naive patients were determined by liquid chromatography-tandem mass spectrometry (LC-MS/MS) after the first dose and at steady state of antitubercular therapy. Patients were genotyped for determination of acetylator status and CYP2C19 phenotype. Nonlinear mixed-effects models were developed to describe the pharmacokinetic data. The model estimated an autoinduction of both rifampin bioavailability (0.5-fold) and clearance (2.3-fold). 25-Desacetylrifampin clearance was 2.1-fold higher at steady state than after the first dose. Additionally, ultrarapid CYP2C19 metabolizers had a 2-fold-higher rifampin clearance at steady state than intermediate or extensive metabolizers. An induction of isonicotinic acid formation from isoniazid dependent on total rifampin dose was estimated. Simulations indicated a 30% lower isoniazid exposure at steady state during administration of standard rifampin doses than isoniazid exposure in noninduced individuals. Rifampin exposure was correlated with CYP2C19 polymorphism, and rifampin administration may increase exposure to toxic metabolites by isoniazid in patients. Both findings may influence the risk of treatment failure, resistance development, and toxicity and require further investigation, especially with regard to ongoing high-dose rifampin trials.

Effect of autoinduction and food on the pharmacokinetics of furmonertinib and its active metabolite characterized by a population pharmacokinetic model

Furmonertinib (AST2818) is a novel third-generation irreversible EGFR TKI and recently has been approved in China for the treatment of non-small cell lung cancer (NSCLC) with EGFR-sensitizing and T790M resistance mutations. In the current study, we developed a semi-mechanistic population pharmacokinetic model to characterize the nonstationary pharmacokinetics (PK) of the furmonertinib and its active metabolite AST5902 simultaneously. The PK data of furmonertinib and AST5902 were obtained from 38 NSCLC patients and 16 healthy volunteers receiving 20-240 mg furmonertinib in three clinical trials. A nonlinear mixed-effects modeling approach was used to describe the PK data. The absorption process of furmonertinib was described by a transit compartment model. The disposition of both furmonertinib and AST5902 was described by a two-compartment model. An indirect response model accounted for the autoinduction of furmonertinib metabolism mediated by CYP3A4. The model-based simulation suggested that furmonertinib clearance was increased in one cycle of treatment (orally once daily for 21 days) compared to baseline, ranging from 1.1 to 1.8 fold corresponding to the dose range of 20-240 mg. The concentration of furmonertinib was decreased over time whereas that of AST5902 was increased. Interestingly, the concentration of the total active compounds (furmonertinib and AST5902) appeared to be stable. The food intake, serum alkaline phosphatase and body weight were identified as statistically significant covariates. The mechanism of food effect on PK was investigated, where the food intake might increase the bioavailability of furmonertinib via increasing the splanchnic blood flow. Overall, a population PK model was successfully developed to characterize the nonstationary PK of furmonertinib and AST5902 simultaneously. The concentrations of total active compounds were less affected by the autoinduction of furmonertinib metabolism.

A semimechanistic population pharmacokinetic and pharmacodynamic model incorporating autoinduction for the dose justification of TAS-114

TAS-114 is a dual deoxyuridine triphosphatase (dUTPase) and dihydropyrimidine dehydrogenase (DPD) inhibitor expected to widen the therapeutic index of capecitabine. Its maximum tolerated dose (MTD) was determined from a safety perspective in a combination study with capecitabine; however, its inhibitory effects on DPD activity were not assessed in the study. The dose justification to select its MTD as the recommended dose in terms of DPD inhibition has been required, but the autoinduction profile of TAS-114 made it difficult. To this end, an approach using a population pharmacokinetic (PPK)/pharmacodynamic (PD) model incorporating autoinduction was planned; however, the utility of this approach in the dose justification has not been reported. Thus, the aim of this study was to demonstrate the utility of a PPK/PD model incorporating autoinduction in the dose justification via a case study of TAS-114. Plasma concentrations of TAS-114 from 185 subjects and those of the endogenous DPD substrate uracil from 24 subjects were used. A two-compartment model with first-order absorption with lag time and an enzyme turnover model were selected for the pharmacokinetic (PK) model. Moreover, an indirect response model was selected for the PD model to capture the changes in plasma uracil concentrations. Model-based simulations provided the dose justification that DPD inhibition by TAS-114 reached a plateau level at the MTD, whereas exposures of TAS-114 increased dose dependently. Thus, the utility of a PPK/PD model incorporating autoinduction in the dose justification was demonstrated via this case study of TAS-114.

Gutsy science: In vitro systems of the human intestine to model oral drug disposition

The intestine has important gate-keeping functions that can profoundly affect the systemic blood exposure of orally administered drugs. Thus, characterizing a new molecular entity's (NME) disposition within the intestine is of utmost importance in drug development. While currently used in vitro systems, such as Ussing chamber, precision-cut intestinal slices, immortalized cell lines, and primary enterocytes provide substantial knowledge about drug absorption and the intestinal first-pass effect, they remain sub-optimal for quantitatively predicting this process and the oral bioavailability of many drugs. Use of novel in vitro systems such as intestinal organoids and intestinal microphysiological systems have provided substantial advances over the past decade, expanding our understanding of intestinal physiology, pathology, and development. However, application of these emerging in vitro systems in the pharmaceutical science is in its infancy. Preliminary work has demonstrated that these systems more accurately recapitulate the physiology and biochemistry of the intact intestine, as it relates to oral drug disposition, and thus they hold considerable promise as preclinical testing platforms of the future. Here we review currently used and emerging in vitro models of the human intestine employed in pharmaceutical science research. We also highlight aspects of these emerging tools that require further study.

Sedative-hypnotic, anxiolytic and possible side effects of Salvia limbata C. A. Mey. Extracts and the effects of phenological stage and altitude on the rosmarinic acid content

ETHNOPHARMACOLOGICAL RELEVANCE

Salvia limbata C. A. Mey. (Persian name: Maryam Goli-e-labeh dar) has been used for treating central nervous disorders such as insomnia, anxiety and depression in Persian traditional medicine. S. limbata is known for its pharmacological activities which could be at least in a part, upon the presence of rosmarinic acid (RA). However, the sedative-hypnotic effect, anxiolytic activity, possible side effects, and the mechanism of action of S. limbata extract has not yet been examined.

AIM OF THE STUDY

In the current study the sedative-hypnotic effect, anxiolytic activity, possible side effects, and the mechanism of action of S. limbata extracts were evaluated. Besides, the effects of altitude and phenological stage on the RA content of S. limbata were investigated.

MATERIALS AND METHODS

Sedative-hypnotic and anxiolytic effects were evaluated through the pentobarbital induced loss of righting reflex test and open field test, respectively. Flumazenil was used to reveal the mechanism of action. Possible side effects were investigated in the passive avoidance and grip strength tests. Besides, the effects of altitude and phenological stage (vegetative, flowering, and seed setting) on the RA content of S. limbata were evaluated using reversed-phase high-performance liquid chromatography (RP-HPLC).

RESULTS

Following behavioral tests, sedative-hypnotic and anxiolytic effects were observed. Since the observed effects were reversed by flumazenil and no side effect on the memory and muscle strength was reported, modulation of the α1-containing GABA-A receptors could be proposed as one of the involved mechanisms. According to the RP-HPLC analysis, harvesting S. limbata in the vegetative stage at the altitude of 2500 m led to the highest content of RA (8.67 ± 0.13 mg/g dry matter). Among different extract of the plant samples collected in the vegetative stage at the altitude of 2500 m, the hydroalcoholic extract showed the highest rosmarinic acid content.

CONCLUSION

The obtained results help to find the optimum situation to gain the highest content of RA as well as the pharmacological activity that could be economically important for the pharmaceutical industries.

Time-dependent enzyme inactivation: Numerical analyses of in vitro data and prediction of drug-drug interactions

Cytochrome P450 (CYP) enzyme kinetics often do not conform to Michaelis-Menten assumptions, and time-dependent inactivation (TDI) of CYPs displays complexities such as multiple substrate binding, partial inactivation, quasi-irreversible inactivation, and sequential metabolism. Additionally, in vitro experimental issues such as lipid partitioning, enzyme concentrations, and inactivator depletion can further complicate the parameterization of in vitro TDI. The traditional replot method used to analyze in vitro TDI datasets is unable to handle complexities in CYP kinetics, and numerical approaches using ordinary differential equations of the kinetic schemes offer several advantages. Improvement in the parameterization of CYP in vitro kinetics has the potential to improve prediction of clinical drug-drug interactions (DDIs). This manuscript discusses various complexities in TDI kinetics of CYPs, and numerical approaches to model these complexities. The extrapolation of CYP in vitro TDI parameters to predict in vivo DDIs with static and dynamic modeling is discussed, along with a discussion on current gaps in knowledge and future directions to improve the prediction of DDI with in vitro data for CYP catalyzed drug metabolism.

Improved Predictions of Drug-Drug Interactions Mediated by Time-Dependent Inhibition of CYP3A

Time-dependent inactivation (TDI) of cytochrome P450s (CYPs) is a leading cause of clinical drug-drug interactions (DDIs). Current methods tend to overpredict DDIs. In this study, a numerical approach was used to model complex CYP3A TDI in human-liver microsomes. The inhibitors evaluated included troleandomycin (TAO), erythromycin (ERY), verapamil (VER), and diltiazem (DTZ) along with the primary metabolites N-demethyl erythromycin (NDE), norverapamil (NV), and N-desmethyl diltiazem (NDD). The complexities incorporated into the models included multiple-binding kinetics, quasi-irreversible inactivation, sequential metabolism, inhibitor depletion, and membrane partitioning. The resulting inactivation parameters were incorporated into static in vitro-in vivo correlation (IVIVC) models to predict clinical DDIs. For 77 clinically observed DDIs, with a hepatic-CYP3A-synthesis-rate constant of 0.000 146 min-1, the average fold difference between the observed and predicted DDIs was 3.17 for the standard replot method and 1.45 for the numerical method. Similar results were obtained using a synthesis-rate constant of 0.000 32 min-1. These results suggest that numerical methods can successfully model complex in vitro TDI kinetics and that the resulting DDI predictions are more accurate than those obtained with the standard replot approach.

Derivation of CYP3A4 and CYP2B6 degradation rate constants in primary human hepatocytes: A siRNA-silencing-based approach

The first-order degradation rate constant (k_deg) of cytochrome P450 (CYP) enzymes is a known source of uncertainty in the prediction of time-dependent drug-drug interactions (DDIs) in physiologically-based pharmacokinetic (PBPK) modelling. This study aimed to measure CYP k_deg using siRNA to suppress CYP expression in primary human hepatocytes followed by incubation over a time-course and tracking of protein expression and activity to observe degradation. The magnitude of gene knockdown was determined by qPCR and activity was measured by probe substrate metabolite formation and CYP2B6-Glo™ assay. Protein disappearance was determined by Western blotting. During a time-course of 96 and 60 h of incubation, over 60% and 76% mRNA knockdown was observed for CYP3A4 and CYP2B6, respectively. The k_deg of CYP3A4 and CYP2B6 protein was 0.0138 h^-1 (±0.0023) and 0.0375 h^-1 (±0.025), respectively. The k_deg derived from probe substrate metabolism activity was 0.0171 h^-1 (±0.0025) for CYP3A4 and 0.0258 h^-1 (±0.0093) for CYP2B6. The CYP3A4 k_deg values derived from protein disappearance and metabolic activity were in relatively good agreement with each other and similar to published values. This novel approach can now be used for other less well-characterised CYPs.

Effects of Strong CYP3A Inhibition and Induction on the Pharmacokinetics of Ixazomib, an Oral Proteasome Inhibitor: Results of Drug-Drug Interaction Studies in Patients With Advanced Solid Tumors or Lymphoma and a Physiologically Based Pharmacokinetic Analysis

At clinically relevant ixazomib concentrations, in vitro studies demonstrated that no specific cytochrome P450 (CYP) enzyme predominantly contributes to ixazomib metabolism. However, at higher than clinical concentrations, ixazomib was metabolized by multiple CYP isoforms, with the estimated relative contribution being highest for CYP3A at 42%. This multiarm phase 1 study (Clinicaltrials.gov identifier: NCT01454076) investigated the effect of the strong CYP3A inhibitors ketoconazole and clarithromycin and the strong CYP3A inducer rifampin on the pharmacokinetics of ixazomib. Eighty-eight patients were enrolled across the 3 drug-drug interaction studies; the ixazomib toxicity profile was consistent with previous studies. Ketoconazole and clarithromycin had no clinically meaningful effects on the pharmacokinetics of ixazomib. The geometric least-squares mean area under the plasma concentration-time curve from 0 to 264 hours postdose ratio (90%CI) with vs without ketoconazole coadministration was 1.09 (0.91-1.31) and was 1.11 (0.86-1.43) with vs without clarithromycin coadministration. Reduced plasma exposures of ixazomib were observed following coadministration with rifampin. Ixazomib area under the plasma concentration-time curve from time 0 to the time of the last quantifiable concentration was reduced by 74% (geometric least-squares mean ratio of 0.26 [90%CI 0.18-0.37]), and maximum observed plasma concentration was reduced by 54% (geometric least-squares mean ratio of 0.46 [90%CI 0.29-0.73]) in the presence of rifampin. The clinical drug-drug interaction study results were reconciled well by a physiologically based pharmacokinetic model that incorporated a minor contribution of CYP3A to overall ixazomib clearance and quantitatively considered the strength of induction of CYP3A and intestinal P-glycoprotein by rifampin. On the basis of these study results, the ixazomib prescribing information recommends that patients should avoid concomitant administration of strong CYP3A inducers with ixazomib.

Accidental Overdose of Oral Cyclosporine in Haematopoietic Stem Cell Transplantation: A Case Report and Literature Review

A 26-year-old woman developed symptoms of acute toxicity during cyclosporine (CsA) therapy for graft-versus-host disease prophylaxis. The standard regimen included CsA in a dose of 1.5 mg/kg (120 mg) every 12 h, but, as a medication error, she received a high dose of 500 mg of oral CsA. After 2 h, she developed nausea and vomiting and, subsequently, flushing, chest tightness, tremor and vertigo. Laboratory and clinical examinations revealed high blood CsA concentrations (1000 ng/mL after 12 h) with a mild increase in blood pressure. Therefore, the patient was diagnosed with an acute CsA overdose. Before confirmation of the overdose by measurement of drug concentrations, the second dose was administered at its routine time because of uncertainty about the aetiology of the symptoms. The third dose was withheld, and the patient was monitored closely for clinical and laboratory presentations until the time when the abnormalities were relieved. CsA administration was then resumed with the correct prescription. The patient was discharged with successful engraftment and normal biochemical laboratory results after 1 month. Evaluation with the Naranjo assessment score indicated a probable relationship between the patient's symptoms and overdosage with the suspected drug. Currently, detailed presentations of acute CsA toxicity cases due to overdose are limited in the medical literature. Evaluation of the patient's medical and laboratory records, with cooperation of all responsible clinical staff, along with a review of the literature, were very helpful in discovery of the toxicity incident. Vigilance of health care providers with regard to medication errors and early detection of toxicity symptoms can decrease CsA-related morbidity and mortality in the future.

Pattern Recognition in Pharmacokinetic Data Analysis

Pattern recognition is a key element in pharmacokinetic data analyses when first selecting a model to be regressed to data. We call this process going from data to insight and it is an important aspect of exploratory data analysis (EDA). But there are very few formal ways or strategies that scientists typically use when the experiment has been done and data collected. This report deals with identifying the properties of a kinetic model by dissecting the pattern that concentration-time data reveal. Pattern recognition is a pivotal activity when modeling kinetic data, because a rigorous strategy is essential for dissecting the determinants behind concentration-time courses. First, we extend a commonly used relationship for calculation of the number of potential model parameters by simultaneously utilizing all concentration-time courses. Then, a set of points to consider are proposed that specifically addresses exploratory data analyses, number of phases in the concentration-time course, baseline behavior, time delays, peak shifts with increasing doses, flip-flop phenomena, saturation, and other potential nonlinearities that an experienced eye catches in the data. Finally, we set up a series of equations related to the patterns. In other words, we look at what causes the shapes that make up the concentration-time course and propose a strategy to construct a model. By practicing pattern recognition, one can significantly improve the quality and timeliness of data analysis and model building. A consequence of this is a better understanding of the complete concentration-time profile.

Severe decrease of cyclosporine levels in a heart transplant recipient receiving the direct thrombin inhibitor argatroban

This case report is about a suspected interaction between argatroban, a direct thrombin inhibitor, and cyclosporine, which occurred in a 60-year-old patient after a second heart transplantation. We explored 4 possible mechanisms of interaction, which are an analytical interference, an idiopathic hemodilution, an increase of renal and hepatic clearance, and a metabolic drug-drug interaction.

Population pharmacokinetics of artesunate and dihydroartemisinin during long-term oral administration of artesunate to patients with metastatic breast cancer

PURPOSE

The purpose of this study were firstly to characterize the population pharmacokinetics of artesunate (ARS) and its active metabolite dihydroartemisinin (DHA) in patients with metastatic breast cancer during long-term (>3 weeks) daily oral ARS administration and secondly to study the relationship between salivary and plasma concentrations of DHA.

METHODS

Drug concentration-time data from 23 patients, receiving oral ARS (100, 150, or 200 mg OD), was analyzed using nonlinear mixed effects modeling. A combined drug-metabolite population pharmacokinetic model was developed to describe the plasma pharmacokinetics of ARS and DHA in plasma. Saliva drug concentrations were incorporated as being directly proportional to plasma concentrations.

RESULTS

A first-order absorption model for ARS linked to a combined two-compartment disposition model for ARS and one-compartment disposition model for DHA provided the best fit to the data. No covariates were identified that could explain between-subject variability. A time-dependent increase in apparent elimination clearance of DHA was observed. Salivary DHA concentrations were proportionally correlated with total DHA plasma concentrations, with an estimated slope factor of 0.116.

CONCLUSIONS

Population pharmacokinetics of ARS and DHA in patients with breast cancer was well described by a combined drug-metabolite model without any covariates and with an increase in apparent elimination clearance of DHA over time. The estimated DHA saliva/plasma ratio was in good agreement with the reported DHA unbound fraction in human plasma. Saliva ARS concentrations correlated poorly with plasma concentrations. This suggests the use of saliva sampling for therapeutic drug monitoring of DHA. However, further studies are warranted to investigate the robustness of this approach.

Pharmacokinetic interaction between nevirapine and nortriptyline in rats: inhibition of nevirapine metabolism by nortriptyline

One of the most frequent comorbidities of HIV infection is depression, with a lifetime prevalence of 22 to 45%. Therefore, it was decided to study a potential pharmacokinetic interaction between the nonnucleoside reverse transcriptase inhibitor nevirapine (NVP) and the tricyclic antidepressant nortriptyline (NT). NVP and NT were administered to rats either orally, intraduodenally, or intravenously, and the changes in plasma levels and pharmacokinetic parameters were analyzed. Experiments with rat and human hepatic microsomes were carried out to evaluate the inhibitory effects of NT on NVP metabolism. NVP plasma concentrations were significantly higher when this drug was coadministered with NT. The maximum plasma concentrations of NVP were increased 2 to 5 times and the total plasma clearance was decreased 7-fold in the presence of NT. However, statistically significant differences in the pharmacokinetic parameters of NT in the absence and presence of NVP were not found. In vitro studies with rat and human hepatic microsomes confirmed the inhibition of NVP hepatic metabolism by NT in a concentration-dependent way, with the inhibition being more intense in the case of rat microsomes. In conclusion, a pharmacokinetic interaction between NVP and NT was detected. This interaction was a consequence of the inhibition of hepatic metabolism of NVP by NT. In vivo human studies are required to evaluate the effects of this interaction on the pharmacokinetics of NVP before it can be taken into account for patients receiving NVP.

Risk assessment of accidental nortriptyline poisoning: the importance of cytochrome P450 for nortriptyline elimination investigated using a population-based pharmacokinetic simulator

It is not possible to make a prospective clinical study that reveals the importance of the nortriptyline metabolising cytochrome P450 (CYP) isoforms (CYP1A2, CYP2C19, CYP2D6, and CYP3A4) in relation to attaining potential toxic nortriptyline concentrations with a possibly fatal outcome. Therefore to study this we have applied the population based pharmacokinetic simulator Simcyp. The objective was to estimate how important CYP2C19 and CYP2D6 phenotype status, hepatic activity of CYP3A4, body weight, CYP2D6 phenotype dose adjustment, and drug-drug interactions are with regard to accidental poisoning in a virtual population receiving a daily dose of 100mg nortriptyline. Accidental poisoning is here defined as intake of a normal dose which because of slow metabolism may lead to potentially toxic concentrations. The input parameters values for Simcyp were based on average literature in vitro and in vivo data. The Simcyp simulations of nortriptyline pharmacokinetics reflected reported clinical concentration-time profiles, therapeutic drug monitoring data, and the consequence of CYP2D6 poor metaboliser (PM) and ultrarapid metaboliser status. Of the investigated factors, the simulations indicate that having CYP2D6 PM status is a major risk factor for attaining high concentrations and thereby possibly becoming poisoned by nortriptyline. Of the CYP2D6 PM subjects 16% would attain plasma concentrations exceeding the toxic limit. Individuals with the combination of CYP2D6 PM status and 10% of the average liver CYP3A4 expression had a 90% risk of becoming poisoned. The results point towards the combination of low CYP3A4 activity and CYP2D6 PM status of major importance for attaining possibly toxic nortriptyline concentrations. In a forensic toxicological context, the results indicate that both the activity of CYP3A4, information on possible drug-drug interactions, and the genotype of CYP2D6 are needed in order to elucidate whether an individual might have been accidentally poisoned because of slow metabolism. In a clinical context, the simulations suggest that precise individual dose adjustment of nortriptyline requires information regarding the activity of both CYP3A4 and CYP2D6. This underlines the value of therapeutic drug monitoring for nortriptyline. Population based pharmacokinetic simulations are considered useful tools for risk assessment in clinical and forensic toxicology.

A mechanism-based population pharmacokinetic model for characterizing time-dependent pharmacokinetics of midostaurin and its metabolites in human subjects

BACKGROUND AND OBJECTIVE

Midostaurin, a novel potent inhibitor of protein kinase C enzyme and class III receptor tyrosine kinases, including Fms-like tyrosine kinase-3 (FLT3) and c-KIT, shows time-dependent pharmacokinetics in human subjects, presumably due to enzyme auto-induction. The purpose of this study was to develop a mechanism-based population pharmacokinetic model to describe the plasma concentration profiles of midostaurin and its metabolites and to characterize the time course of auto-induction.

SUBJECTS AND METHODS

Data from 37 diabetic patients who received oral doses of midostaurin (25 mg twice daily, 50 mg twice daily or 75 mg twice daily) for 28 days were analysed using nonlinear mixed-effects modelling. The structural model included a gut compartment for drug input and central and peripheral compartments for midostaurin, with drug output from the central compartment to either of two compartments for the midostaurin metabolites CGP62221 and CGP52421. Different enzyme induction sub-models were evaluated to account for the observed time-dependent decrease in midostaurin concentrations.

RESULTS

An enzyme turnover model, with CGP62221 formation (CL(1)) being a linear process but CGP52421 formation (CL(2)) being inducible, was found to be most appropriate. In the pre-induced state, CL(1) and CL(2) of midostaurin were determined to be 1.47 L/h and 0.501 L/h, respectively. At the end of 28 days of dosing, CL(2) was increased by 5.2-, 6.6- and 6.9-fold in the 25 mg, 50 mg and 75 mg groups, respectively, resulting in a 2.1- to 2.5-fold increase in total clearance of midostaurin. The final model estimated a mean maximum fold of induction (E(max)) of 8.61 and a concentration producing 50% of the E(max) (EC(50)) of 1700 ng/mL (approximately 2.9 micromol/L) for CGP52421-mediated enzyme induction.

CONCLUSIONS

The population pharmacokinetic model that was developed was able to describe the time-dependent pharmacokinetic profiles of midostaurin and its auto-induction mechanism. Thus it may be useful for designing an appropriate dosage regimen for midostaurin. The unique feature of this model included a precursor compartment that was able to capture the time delays of auto-induction. The use of such precursor extension in the model may be applicable to other drugs showing long time delays in enzyme auto-induction.

Pharmacodynamics of carbamazepine-mediated induction of CYP3A4, CYP1A2, and Pgp as assessed by probe substrates midazolam, caffeine, and digoxin

The aim of this study was to develop a model describing the carbamazepine autoinduction and the carbamazepine-mediated induction of CYP3A4, CYP1A2, and P-glycoprotein. Seven healthy volunteers were dosed with carbamazepine over 16 consecutive days. The CYP3A4, CYP1A2, and P-glycoprotein activities were assessed, using midazolam, caffeine, and digoxin as probe substrates, on 12 occasions, covering the preinduced state and the onset and termination of the induction process. The data were evaluated using a mechanistic pharmacokinetic approach in NONMEM. The induction processes were described using turnover models, with carbamazepine and carbamazepine-10,11-epoxide as the driving force of the induction. The half-lives of CYP3A4 and CYP1A2 were estimated to be 70 and 105 h, respectively. P-glycoprotein was not affected by the carbamazepine treatment. The possibility of modeling the pharmacodynamics of enzyme induction using a turnover model was illustrated, and the time course of the process was estimated with good precision.

Dealing with time-dependent pharmacokinetics during the early clinical development of a new leukotriene B4 synthesis inhibitor

PURPOSE

The aim of this study was to explore the possibility of achieving a practical dosing regimen for 2,4,6-triiodophenol (AM-24), a new leukotriene B4 (LTB4) synthesis inhibitor. First, a model capable of dealing with the nonlinearity in its pharmacokinetic profile was built, and then it was combined with a pharmacodynamic model previously established with data from earlier phase I trials.

METHODS

One week after the first 240-, 350-, or 500-mg oral dose of AM-24, six additional doses were given to 24 healthy volunteers once daily. A total of 33 blood samples were obtained from each individual. Different models, including enzyme turnover models, were fitted to the data by using the software NONMEM.

RESULTS

Drug absorption was modeled with a first-order process. Drug disposition was described with a one-compartment model, and elimination with an (auto)inhibited and a noninhibited clearance. AM-24 inhibited the enzyme production rate to a maximum of 98%. Relative bioavailability was independent of the decrease in the amount of enzyme. The estimate of the enzyme turnover half-life was 8.5 h.

CONCLUSIONS

Simulations have shown that steady-state conditions eliciting 90% of maximal LTB4 synthesis inhibition can be reached after 3 weeks during an oral treatment with AM-24 administered at the dosage of 500 mg once daily.

A mechanism-based integrated pharmacokinetic enzyme model describing the time course and magnitude of phenobarbital-mediated enzyme induction in the rat

PURPOSE

To characterize the magnitude, time course, and specificity of phenobarbital (PB)-mediated enzyme induction, and further, to develop an integrated pharmacokinetic (PK)-enzyme model describing the changes in the activities of CYP enzymes as well as in the PK of PB.

METHODS

PB plasma concentrations and in vitro activities of several CYP enzymes were measured in rats treated with PB between 0 and 14 days. A PB PK-enzyme induction model was developed using the program NONMEM: .

RESULTS

PB treatment both induces and reduces the activity of CYP enzymes by stimulating the enzymes' formation or elimination rates. Certain CYP enzymes affected the PB PK through autoinduction. The half-life of the induction process was estimated to be 2 days for CYP1A2, CYP3A1/2, and CYP2B1/2, and 3 days for androstenedione producing enzymes. The CYP2C11 activity was rapidly reduced by PB treatment. A lag time for the PB autoinduction was observed. This lag time is explained by the rate difference between induction and reduction in CYP activities.

CONCLUSION

To our knowledge, this is the first example of an induction model that simultaneously describes plasma PK and in vitro data. It does so by integrating the bidirectional interaction between drug and enzymes in a mechanistic manner.

The Metabolic Fate of Amitriptyline, Nortriptyline and Amitriptylinoxide in Man

Amitriptyline (AT), the most widely used tricyclic antidepressant, undergoes oxidative metabolism in the side chain with production of the secondary amine nortriptyline (NT), a primary amine, and the N-oxide amitriptylinoxide (AT-NO); in addition, direct conjugation leads to a quaternary ammonium-linked glucuronide. Hydroxylation of AT or NT at the ethylene bridge of the central seven-membered ring results in four isomeric alcohols and occurs with high stereo- and enantioselectivity, the (-)-(E)-10-hydroxy compounds usually being the major products. The disposition of the alcohols is also partially enantioselective, for instance with regard to glucuronidation and reversible oxidation to ketones. Introduction of a second hydroxy group results in isomeric glycols. Oxidative attack at an aromatic ring is a minor pathway leading to dihydrodiols and phenols. Numerous metabolites originate by combinations of reactions in the ring system and the side chain. AT-NO is by about one-third excreted in unchanged form or as 10-hydroxy derivative; the major part is reduced to AT and metabolized further. The review covers current knowledge on the enzymes participating in the individual pathways. Their quantitative importance is inferred from kinetic studies in volunteers and patients and from experiments in vitro. Clinical consequences of biochemical findings mainly derive from the impact of the polymorphic CYP2D6 mediating (-)-(E)-10-hydroxylation and from its potential inhibition by other psychoactive drugs.

Role of variability in explaining ethanol pharmacokinetics: research and forensic applications

Variability in the rate and extent of absorption, distribution and elimination of ethanol has important ramifications in clinical and legal medicine. The speed of absorption of ethanol from the gut depends on time of day, drinking pattern, dosage form, concentration of ethanol in the beverage, and particularly the fed or fasting state of the individual. During the absorption phase, a concentration gradient exists between the stomach, portal vein and the peripheral venous circulation. First-pass metabolism and bioavailability are difficult to assess because of dose-, time- and flow-dependent kinetics. Ethanol is transported by the bloodstream to all parts of the body. The rate of equilibration is governed by the ratio of blood flow to tissue mass. Arterial and venous concentrations differ as a function of time after drinking. Ethanol has low solubility in lipids and does not bind to plasma proteins, so volume of distribution is closely related to the amount of water in the body, contributing to sex- and age-related differences in disposition. The bulk of ethanol ingested (95-98%) is metabolised and the remainder is excreted in breath, urine and sweat. The rate-limiting step in oxidation is conversion of ethanol into acetaldehyde by cytosolic alcohol dehydrogenase (ADH), which has a low Michaelis-Menten constant (Km) of 0.05-0.1 g/L. Moreover, this enzyme displays polymorphism, which accounts for racial and ethnic variations in pharmacokinetics. When a moderate dose is ingested, zero-order elimination operates for a large part of the blood-concentration time course, since ADH quickly becomes saturated. Another ethanol-metabolising enzyme, cytochrome P450 2E1, has a higher Km (0.5-0.8 g/L) and is also inducible, so that the clearance of ethanol is increased in heavy drinkers. Study design influences variability in blood ethanol pharmacokinetics. Oral or intravenous administration, or fed or fasted state, might require different pharmacokinetic models. Recent work supports the need for multicompartment models to describe the disposition of ethanol instead of the traditional one-compartment model with zero-order elimination. Moreover, appropriate statistical analysis is needed to isolate between- and within-subject components of variation. Samples at low blood ethanol concentrations improve the estimation of parameters and reduce variability. Variability in ethanol pharmacokinetics stems from a combination of both genetic and environmental factors, and also from the nonlinear nature of ethanol disposition, experimental design, subject selection strategy and dose dependency. More work is needed to document variability in ethanol pharmacokinetics in real-world situations.

Decreased plasma levels of amitriptyline and its metabolites on comedication with an extract from St. John's wort ( Hypericum perforatum )

Extracts of St. John's wort ( Hypericum perforatum ) became increasingly popular as easily available remedies for mild to moderate depression. Comedication with hypericum extract was recently shown to drastically reduce plasma concentration of ciclosporin, digoxin, and indinavir. We investigated the possible interaction of hypericum extract LI160 with amitriptyline. Both antidepressants have a high probability of concomitant use. Twelve patients requiring amitriptyline treatment received a single dose of hypericum extract (900 mg) at day 1, continued by a 12-to 14-day treatment with retarded amitriptyline (75 mg twice daily). Then hypericum (900 mg/day) was added for another 14 to 16 days. Steady-state pharmacokinetics of amitriptyline were compared before and after multiple-dose treatment with hypericum extract. Furthermore, comparisons were made for single-dose kinetics of hypericum-extract ingredients hypericin, pseudohypericin, and hyperforin between the first day of concomitant treatment and LI160 alone. Multiple-dose comedication with LI160 led to a statistically significant decrease in the area under the plasma concentration-time curve within one dosing interval of amitriptyline by 22% ( p = 0.03) and nortriptyline by 41% ( p = 0.002), as well as of all hydroxylated metabolites, except for 10-E-hydroxynortriptyline. Plasma levels of amitriptyline and hydroxylated metabolites gradually decreased, whereas nortriptyline concentrations were already markedly decreased after 3 days of cotreatment with hypericum. Cumulative urinary amounts of amitriptyline and metabolites decreased to the same extent as plasma concentrations upon hypericum comedication. Induction of cytochrome P-450 enzymes or drug transporters (P-glycoprotein) by St. John's wort extract may explain this pharmacokinetic interaction. Physicians should be aware of this interaction when treating patients with amitriptyline.

Interindividual variability in inhibition and induction of cytochrome P450 enzymes

Drug interactions have always been a major concern in medicine for clinicians and patients. Inhibition and induction of cytochrome P450 (CYP) enzymes are probably the most common causes for documented drug interactions. Today, many pharmaceutical companies are predicting potential interactions of new drug candidates. Can in vivo drug interactions be predicted accurately from in vitro metabolic studies? Should the prediction be qualitative or quantitative? Although some scientists believe that quantitative prediction of drug interactions is possible, others are less optimistic and believe that quantitative prediction would be very difficult. There are many factors that contribute to our inability to quantitatively predict drug interactions. One of the major complicating factors is the large interindividual variability in response to enzyme inhibition and induction. This review examines the sources that are responsible for the interindividual variability in inhibition and induction of cytochrome P450 enzymes.

A mechanism-based pharmacokinetic model for the cytochrome P450 drug-drug interaction between cyclophosphamide and thioTEPA and the autoinduction of cyclophosphamide

Cyclophosphamide (CP) is widely used in high-dose chemotherapy regimens in combination with thioTEPA. CP is a prodrug and is activated by cytochrome P450 to 4-hydroxycyclophosphamide (HCP) which yields the final cytotoxic metabolite phosphoramide mustard (PM). The metabolism of CP into HCP exhibits autoinduction but is inhibited by thioTEPA. The aim of this study was to develop a population pharmacokinetic model for the bioactivation route of CP incorporating the phenomena of both autoinduction and the drug-drug interaction between CP and thioTEPA. Plasma samples were collected from 34 patients who received high-dose CP, thioTEPA and carboplatin in short infusions during 4 consecutive days. Elimination of CP was described by a noninducible route and an inducible route leading to HCP. The latter route was mediated by a hypothetical amount of enzyme. Autoinduction leads to a zero-order increase in amount of this enzyme during treatment. Inhibition by thioTEPA was modeled as a reversible, competitive, concentration-dependent inhibition. PM pharmacokinetics were described by first-order formation from HCP and first-order elimination. The final models for CP, HCP, and PM provided an adequate fit of the experimental data. The volume of distribution, noninducible and initial inducible clearances of CP were 31.0 L, 1.58 L/hr and 4.76 L/hr, respectively. The enzyme amount increased with a zero-order rate constant of 0.041 amount * hr-1. After each thioTEPA infusion, however, approximately 80% of the enzyme was inhibited. This inhibition was reversible with a half-life of 6.5 hr. The formation and elimination rate constants of PM were 1.58 and 0.338 hr-1, respectively. The developed model enabled the assessment of the complex pharmacokinetics of CP in combination with thio TEPA. This model provided an adequate description of enzyme induction and inhibition and can be used for treatment optimization in this combination.

Clinical pharmacokinetics and pharmacodynamics of ifosfamide and its metabolites

This review discusses several issues in the clinical pharmacology of the antitumour agent ifosfamide and its metabolites. Ifosfamide is effective in a large number of malignant diseases. Its use, however, can be accompanied by haematological toxicity, neurotoxicity and nephrotoxicity. Since its development in the middle of the 1960s, most of the extensive metabolism of ifosfamide has been elucidated. Identification of specific isoenzymes responsible for ifosfamide metabolism may lead to an improved efficacy/toxicity ratio by modulation of the metabolic pathways. Whether ifosfamide is specifically transported by erythrocytes and which activated ifosfamide metabolites play a key role in this transport is currently being debated. In most clinical pharmacokinetic studies, the phenomenon of autoinduction has been observed, but the mechanism is not completely understood. Assessment of the pharmacokinetics of ifosfamide and metabolites has long been impaired by the lack of reliable bioanalytical assays. The recent development of improved bioanalytical assays has changed this dramatically, allowing extensive pharmacokinetic assessment, identifying key issues such as population differences in pharmacokinetic parameters, differences in elimination dependent upon route and schedule of administration, implications of the chirality of the drug and interpatient pharmacokinetic variability. The mechanisms of action of cytotoxicity, neurotoxicity, urotoxicity and nephrotoxicity have been pivotal issues in the assessment of the pharmacodynamics of ifosfamide. Correlations between the new insights into ifosfamide metabolism, pharmacokinetics and pharmacodynamics will rationalise the further development of therapeutic drug monitoring and dose individualisation of ifosfamide treatment.

Evaluation of the autoinduction of ifosfamide metabolism by a population pharmacokinetic approach using NONMEM

AIMS

This study investigated the population pharmacokinetics of ifosfamide in 15 patients treated for soft tissue sarcoma with 9 or 12 g m-2 ifosfamide by means of a 72 h continuous i.v. infusion.

METHODS

A model was developed using nonlinear mixed effects modelling (NONMEM) to describe the nonlinear pharmacokinetics of ifosfamide by linking the ifosfamide plasma concentrations to the extent of the autoinduction.

RESULTS

The proposed model revealed the effect of autoinduction on the disposition of ifosfamide. The initial clearance, volume of distribution, rate constant for enzyme degradation, induction half-life of the enzyme and the ifosfamide concentration at 50% of the maximum inhibition of enzyme degradation were estimated at 2.94 +/- 0.27 l h-1, 43.5 +/- 2.9 l, 0.0546 +/- 0. 0078 h-1, 12.7 h and 30.7 +/- 4.8 microM, respectively. Interindividual variabilities of initial clearance, volume of distribution, rate constant for enzyme degradation were 24.5, 23.4 and 22.7%, respectively. Proportional and additive variability not explained by the model were 13.6% and 0.0763 microM, respectively.

CONCLUSIONS

The absence of a lag time for the autoinduction of ifosfamide metabolism could be the result of an immediate inhibition of the enzymatic degradation of CYP3A4 by ifosfamide. By application of the autoinduction model individual pharmacokinetic profiles of patients were described with adequate precision. This model may therefore be used in the future development of a model to individualize dose selection in patients.

A therapeutic use of the methadone fluvoxamine drug interaction

BACKGROUND

Methadone has been effectively used in the treatment of opiate dependence. Adequate dose and blood level have correlated with success in treatment. A number of factors including the regular use of alcohol, medications, and urinary pH can influence blood level and thereby effectiveness. Fluvoxamine has been shown to increase methadone blood levels.

METHODS

Single case report.

RESULTS

A patient unable to maintain an effective methadone blood level despite a dose of 200 mg per day was administered fluvoxamine with subsequent increase in her methadone blood level and reduction of opiate withdrawal symptoms.

CONCLUSIONS

In patients unable to maintain an effective methadone blood level throughout the dosing interval, fluvoxamine can help increase the methadone blood level and alleviate opiate withdrawal symptoms.

A mechanism-based pharmacokinetic-enzyme model for cyclophosphamide autoinduction in breast cancer patients

AIMS

This study investigated the pharmacokinetics of cyclophosphamide (CP) and its main metabolite 4-hydroxycyclophosphamide (4-OH-CP) in patients with breast cancer undergoing high dose chemotherapy prior to autologous stem cell transplantation. An enzyme turn-over model was also developed to study the time course of cyclophosphamide induction.

METHODS

Fourteen patients received a combination of CP (6 g m-2 ), thiotepum (500 mg m-2 ) and carboplatin (800 mg m-2 ) as a 96 h infusion. Plasma concentrations of CP and 4-OH-CP were determined with h.p.l.c. and a pharmacokinetic and enzyme turn-over model applied to data using NONMEM.

RESULTS

CP plasma concentrations were described by a two-compartment model with a noninducible and an inducible pathway, the latter forming 4-OH-CP. In the final enzyme model, CP affects the amount of enzymes by increasing the enzyme production rate. CP concentrations decreased during the infusion with no subsequent change in 4-OH-CP concentrations. CP inducible and noninducible clearance were estimated to 1.76 l h-1 (90% C.I. 0.92-2.58) and 1.14 l h-1 (0.31-1.85), respectively. The induction resulted in an approximately doubled CP clearance through the inducible pathway at the end of treatment. The model predicted the enzyme turn-over half-life to be 24 h.

CONCLUSIONS

The presented mechanism-based enzyme induction model where the pharmacokinetics of the inducer and the enzyme pool counterbalance each other successfully described CP autoinduction. It is reasonable to believe that CP affects its own elimination by increasing the enzyme production rate and thereby increasing the amount of enzyme by which CP is eliminated.

Population pharmacokinetics of methadone in opiate users: characterization of time-dependent changes

AIMS

Although methadone is widely used to treat opiate dependence, guidelines for its dosage are poorly defined. There is increasing evidence to suggest that a strategy based on plasma drug monitoring may be useful to detect non-compliance. Therefore, we have developed a population-based pharmacokinetic (POP-PK) model that characterises adaptive changes in methadone kinetics.

METHODS

Sparse plasma rac-methadone concentrations measured in 35 opiate-users were assessed using the P-Pharm software. The final structural model comprised a biexponential function with first-order input and allowance for time-dependent change in both clearance (CL) and initial volume of distribution (V ). Values of these parameters were allowed to increase or decrease exponentially to an asymptotic value.

RESULTS

Increase in individual values of CL and increase or decrease in individual values of V with time was observed in applying the model to the experimental data.

CONCLUSIONS

A time-dependent increase in the clearance of methadone is consistent with auto-induction of CYP3A4, the enzyme responsible for much of the metabolism of the drug. The changes in V with time might reflect both up- and down-regulation of alpha1-acid glycoprotein, the major plasma binding site for methadone. By accounting for adaptive kinetic changes, the POP-PK model provides an improved basis for forecasting plasma methadone concentrations to predict and adjust dosage of the drug and to monitor compliance in opiate-users on maintenance treatment.

Nortriptyline E-10-hydroxylation in vitro is mediated by human CYP2D6 (high affinity) and CYP3A4 (low affinity): implications for interactions with enzyme-inducing drugs

The human cytochrome P450 (CYP) isoforms mediating nortriptyline 10-hydroxylation have been identified using kinetic studies on heterologously expressed human CYPs and chemical inhibition studies on human liver microsomes. Nortriptyline was metabolized to E-10-hydroxynortriptyline by human lymphoblast-expressed CYPs 2D6 (Km 2.1 microM) and 3A4 (Km 37.4 microM) with high and low affinity, respectively, whereas CYPs 1A2, 2A6, 2B6, 2C9, 2C19, and 2E1 had no detectable activity. Human liver microsomal nortriptyline E-10-hydroxylation displayed biphasic kinetics. The high-affinity component (Km 1.3 +/- 0.4 microM, n = 11 livers) was selectively inhibited by the CYP 2D6 inhibitor quinidine, whereas the CYP3A4 inhibitor ketoconazole selectively inhibited the low-affinity component (K(m) 24.4 +/- 7 microM, n = 11 livers). Inhibition by ketoconazole increased with increasing substrate concentration, whereas the reverse was true for quinidine. The Vmax of the low-affinity component in human liver microsomes was significantly correlated (r2 = 0.84) with the relative activity factor for CYP3A4, a measure of the amount of catalytically active enzyme. A simulation of the relative contribution of CYPs 2D6 and 3A4 to net nortriptyline hydroxylation rate suggested that the relative contribution of CYP3A4 is only 20% even at the higher end of the therapeutic range. Induction of CYP3A4 will increase its importance and increase the net metabolic rate, whereas inhibition of CYP3A4 will be of little importance due to its minimal relative contribution under uninduced conditions. The identification of CYP3A4 as a low-affinity nortriptyline E-10-hydroxylase explains the ability of poor metabolizers of debrisoquin to hydroxylate nortriptyline, as well as the increased in vivo clearance via this pathway caused by CYP3A4-inducing drugs such as pentobarbital, carbamazepine, and rifampin.