Impact of Heterotropic Allosteric Modulation on the Time-Dependent Inhibition of Cytochrome P450 3A4.

The current study was designed to investigate the influence of allosteric effectors on the metabolism of the prototypical cytochrome P450 (CYP) 3A4 substrate midazolam (MDZ), and on the determination in vitro time-dependent inhibition (TDI) of CYP3A4 using human liver microsomes (HLM). As the concentration of midazolam increased to 250 µM in HLMs, homotropic cooperativity resulted in a decrease in the 1'-hydroxymidazolam to 4-hydroxymidazolam ratio to a maximum of 1.1. The presence of varying concentrations of testosterone, progesterone (PGS), or carbamazepine (CBZ) in HLMs with MDZ could recapitulate the effect of homotropic cooperativity such that the formation rates of the 1'hydroxymidazolam and 4-hydroxymidazolam were equal even at low concentrations of MDZ. The presence of PGS (10 or 100 µM) and CBZ (100 or 1000 µM) in in vitro TDI determination of four known CYP3A4 time-dependent inactivators (clarithromycin, troleandomycin, mibefradil, raloxifene) simultaneously decreased potency and inactivation rate constant, resulting in fold changes in inactivation efficiency on average of 1.6-fold and 13-fold for the low and high concentrations of allosteric modulator tested, respectively. The formation of a metabolic-intermediate complex (MIC) for clarithromycin and troleandomycin decreased in the presence of the allosteric modulators in a concentration-dependent manner, reaching a new steady state formation that could not be overcome with increased incubation time. Maximum reduction of the MIC formed by clarithromycin was up to ∼91%, while troleandomycin MIC decreased up to ∼31%. These findings suggest that the absence of endogenous allosteric modulators may contribute to the poor translation of HLM-based drug-drug interaction predictions. SIGNIFICANCE STATEMENT: The reported overprediction of in vitro human liver microsome time-dependent inhibition of CYP3A4 and observed drug interactions in vivo remains an issue in drug development. We provide characterization of allosteric modulators on the CYP3A4 metabolism of the prototypical substrate midazolam, demonstrating the ability of the modulators to recapitulate the homotropic cooperativity of midazolam. Furthermore, we demonstrate that allosteric heterotropic cooperativity of CYP3A4 can impact the time-dependent inhibition kinetics of known mechanisms-based inhibitors, providing a potential mechanism to explain the overprediction.

Prediction of the Renal Organic Anion Transporter 1 (OAT1)- Mediated Drug Interactions for LY404039, the Active Metabolite of Pomaglumetad Methionil.

PURPOSE: The objective of this work was to demonstrate that clinical OAT1-mediated DDIs can be predicted using physiologically based pharmacokinetic (PBPK) modeling. METHODS: LY404039 is a metabotropic glutamate receptor 2/3 agonist and the active moiety of the prodrug pomaglumetad methionil (LY2140023). After oral administration, pomaglumetad methionil is rapidly taken up by enterocytes via PEPT1 and once absorbed, converted to LY404039 via membrane dehydropeptidase 1 (DPEP1). LY404039 is renally excreted by both glomerular filtration and active secretion and in vitro studies showed that the active secretion of LY404039 was mediated by the organic anion transporter 1 (OAT1). Both clinical and in vitro data were used to build a PBPK model to predict OAT1-mediated DDIs. RESULTS: In vitro inhibitory potencies (IC50) of the known OAT inhibitors, probenecid and ibuprofen, were determined to be 4.00 and 2.63 µM, respectively. Subsequently, clinical drug-drug interaction (DDI) study showed probenecid reduced the renal clearance of LY404039 by 30 to 40%. The PBPK bottom-up model, predicted a renal clearance that was approximately 20% lower than the observed one. The middle-out model, using an OAT1 relative activity factor (RAF) of 3, accurately reproduced the renal clearance of LY404039 and pharmacokinetic (PK) changes of LY404039 in the presence of probenecid. CONCLUSIONS: OAT1- mediated DDIs can be predicted using in vitro measured IC50 and PBPK modeling. The effect of ibuprofen was predicted to be minimal (AUC ratio of 1.15) and not clinically relevant.

Abemaciclib Does Not Have a Clinically Meaningful Effect on Pharmacokinetics of CYP1A2, CYP2C9, CYP2D6, and CYP3A4 Substrates in Patients with Cancer.

Abemaciclib is an orally administered, potent inhibitor of cyclin-dependent kinases 4 and 6 and is metabolized extensively by CYP3A4. The effects of abemaciclib on several CYPs were qualified in vitro and subsequently evaluated in a clinical study. In vitro, human hepatocytes were treated with vehicle, abemaciclib, or abemaciclib metabolites [N-desethylabemaciclib (M2) or hydroxyabemaciclib (M20)]. mRNA levels for eight CYPs were measured using reverse-transcription quantitative polymerase chain reaction, and, additionally, catalytic activities for three CYPs were determined. In the clinical study, adult patients with cancer received a drug cocktail containing CYP substrates [midazolam (3A), warfarin (2C9), dextromethorphan (2D6), and caffeine (1A2)] either alone or in combination with abemaciclib. Plasma pharmacokinetics (PK) samples were analyzed for all substrates, caffeine metabolite paraxanthine, and abemaciclib; polymorphisms of CYP2C9, CYP2D6, CYP3A4, and CYP3A5 were evaluated. In vitro, downregulation of CYP mRNA, including 1A2, 2B6, 2C8, 2C9, 2D6, and 3A, by abemaciclib and/or M2 and M20 was observed at clinically relevant concentrations. In humans, abemaciclib did not affect the PK of CYP2D6 or CYP2C9 substrates. Minor statistically significant but clinically irrelevant changes were observed for midazolam [area under the concentration versus time curve from zero to infinity (AUC0-inf) (13% lower), Cmax (15% lower)], caffeine [AUC0-inf (56% higher)], and paraxanthine: caffeine [area under the concentration versus time curve from 0 to 24 hours ratio (was approximately 30% lower)]. However, given the magnitude of the effect, these changes are not considered clinically relevant. In conclusion, the downregulation of CYP mRNA mediated by abemaciclib in vitro did not translate into clinically meaningful drug-drug interactions in patients with cancer. SIGNIFICANCE STATEMENT: Despite observations that abemaciclib alters the mRNA of various CYP isoforms in vitro, a clinical study using a drug cocktail approach found no clinically meaningful drug-drug interactions between abemaciclib and a range of CYP substrates [midazolam (CYP3A4), S-warfarin (CYP2C9), dextromethorphan (CYP2D6), and caffeine (CYP1A2)]. This lack of translation suggests greater understanding of mechanisms of CYP downregulation is needed to accurately predict clinical drug-drug interaction risk from in vitro data.

Hepatic Organic Anion Transporting Polypeptide-Mediated Clearance in the Beagle Dog: Assessing In Vitro-In Vivo Relationships and Applying Cross-Species Empirical Scaling Factors to Improve Prediction of Human Clearance.

In the present study, the beagle dog was evaluated as a preclinical model to investigate organic anion transporting polypeptide (OATP)-mediated hepatic clearance. In vitro studies were performed with nine OATP substrates in three lots of plated male dog hepatocytes ± OATP inhibitor cocktail to determine total uptake clearance (CLuptake) and total and unbound cell-to-medium concentration ratio (Kpuu). In vivo intrinsic hepatic clearances (CLint,H) were determined following intravenous drug administration (0.1 mg/kg) in male beagle dogs. The in vitro parameters were compared with those previously reported in plated human, monkey, and rat hepatocytes; the ability of cross-species scaling factors to improve prediction of human in vivo clearance was assessed. CLuptake in dog hepatocytes ranged from 9.4 to 135 µl/min/106 cells for fexofenadine and telmisartan, respectively. Active process contributed >75% to CLuptake for 5/9 drugs. Rosuvastatin and valsartan showed Kpuu > 10, whereas cerivastatin, pitavastatin, repaglinide, and telmisartan had Kpuu < 5. The extent of hepatocellular binding in dog was consistent with other preclinical species and humans. The bias (2.73-fold) obtained from comparison of predicted versus in vivo dog CLint,H was applied as an average empirical scaling factor (ESFav) for in vitro-in vivo extrapolation of human CLint,H The ESFav based on dog reduced underprediction of human CLint,H for the same data set (geometric mean fold error = 2.1), highlighting its utility as a preclinical model to investigate OATP-mediated uptake. The ESFav from all preclinical species resulted in comparable improvement of human clearance prediction, in contrast to drug-specific empirical scalars, rationalized by species differences in expression and/or relative contribution of particular transporters to drug hepatic uptake.

The Effect of Promiscuous Aggregation on in Vitro Drug Metabolism Assays.

PURPOSE: Many bioactive molecules show a type of solution phase behavior, termed promiscuous aggregation, whereby at micromolar concentrations, colloidal drug-rich aggregates are formed in aqueous solution. These aggregates are known to be a major cause of false positives and false negatives in select enzymatic high-throughput screening assays. The goal of this study was to investigate the impact of drug-rich aggregates on in vitro drug screening metabolism assays. METHODS: Cilnidipine was selected as an aggregate former and its impact on drug metabolism was evaluated against rCYP2D6, rCYP1A2, rCYP2C9 and human liver microsomes. RESULTS: The cilnidipine aggregates were shown to non-specifically inhibit multiple cytochrome P450 enzymes with an IC50 comparable with the IC50 of potent model inhibitors. CONCLUSIONS: This newly demonstrated mode of "promiscuous inhibition" is of great importance as it can lead to false positives during drug metabolism evaluations and thus it needs to be considered in the future to better predict in vivo drug-drug interactions.

Transient Alpha and Beta Synchrony Underlies Preparatory Recruitment of Directional Motor Networks.

Modulations in motor cortical beta and alpha activity have been implicated in the preparation, execution, and termination of voluntary movements. The functional role of motor cortex beta activity is yet to be defined, though two opposing theories prevail. The idling cortex theory suggests that large-scale motor networks, in the absence of input, revert to an intrinsic oscillatory state. The alternative theory proposes that beta activity promotes postural tone at the expense of voluntary movement. These theories are primarily based on observations of event-related desynchronization associated with movement onset. Here, we explore the changes in alpha and beta oscillatory activity associated with the specific behavioral patterns during an established directional uncertainty paradigm. We demonstrate that, consistent with current proposals, alpha and beta desynchronization reflects a process of disengagement from existing networks to enable the creation of functional assemblies. We demonstrate that, following desynchronization, a novel signature of transient alpha synchrony underlies the recruitment of functional assemblies required for directional control. Although alpha and beta desynchronization are dependent upon the number of cues presented, they are not predictive of movement preparation. However, the transient alpha synchrony occurs only when participants have sufficient information to prepare for movement and shows a direct relationship with behavioral performance measures.

The Impact of the Hepatocyte-to-Plasma pH Gradient on the Prediction of Hepatic Clearance and Drug-Drug Interactions for CYP2C9 and CYP3A4 Substrates.

Surrogate assays for drug metabolism and inhibition are traditionally performed in buffer systems at pH 7.4, despite evidence that hepatocyte intracellular pH is 7.0. This pH gradient can result in a pKa-dependent change in intracellular/extracellular concentrations for ionizable drugs that could affect predictions of clearance and P450 inhibition. The effect of microsomal incubation pH on in vitro enzyme kinetic parameters for CYP2C9 (diclofenac, (S)-warfarin) and CYP3A4 (midazolam, dextromethorphan, testosterone) substrates, enzyme specific reversible inhibitors (amiodarone, desethylamiodarone, clozapine, nicardipine, fluconazole, fluvoxamine, itraconazole) and a mechanism-based inhibitor (amiodarone) was investigated. Intrinsic clearance through CYP2C9 significantly increased (25% and 50% for diclofenac and (S)-warfarin respectively) at intracellular pH 7.0 compared with traditional pH 7.4. The CYP3A4 substrate dextromethorphan intrinsic clearance was decreased by 320% at pH 7.0, while midazolam and testosterone remained unchanged. Reversible inhibition of CYP2C9 was less potent at pH 7.0 compared with 7.4, while CYP3A4 inhibition potency was variably affected. Maximum enzyme inactivation rate of amiodarone toward CYP2C9 and CYP3A4 decreased at pH 7.0, while the irreversible inhibition constant remained unchanged for CYP2C9, but decreased for CYP3A4 at pH 7.0. Predictions of clearance and drug-drug interactions made through physiologically based pharmacokinetic models were improved with the inclusion of predicted intracellular concentrations based at pH 7.0 and in vitro parameters determined at pH 7.0. No general conclusion on the impact of pH could be made and therefore a recommendation to change buffer pH to 7.0 cannot be made at this time. It is recommended that the appropriate hepatocyte intracellular pH 7.0 be used for in vitro determinations when in vivo predictions are made.

The Impact of the Hepatocyte-to-Plasma pH Gradient on the Prediction of Hepatic Clearance and Drug-Drug Interactions for CYP2D6 Substrates.

The proton gradient from the intracellular space to plasma creates an unbound drug gradient for weak acids and bases that could modulate apparent drug clearance and drug-drug interactions. Cytochrome P450 intrinsic clearance and inhibitor potency are routinely determined in vitro at the plasma pH of 7.4 rather than the intrahepatocyte pH of 7.0. We determined the impact of pH on in vitro enzyme kinetic parameters and inhibition potency for substrates (bufuralol, dextromethorphan), reversible inhibitors (quinidine, amiodarone, desethylamiodarone, clozapine), and mechanism-based inhibitors (paroxetine, desethylamiodarone) of the major drug metabolizing-enzyme CYP2D6. The lower intracellular pH 7.0 compared with pH 7.4 resulted in a 60 and 50% decrease in intrinsic clearance for the substrates bufuralol and dextromethorphan, respectively. Reversible inhibition constants for three of the four inhibitors tested were unaffected by pH, whereas for the inhibitor quinidine, a 2-fold increase in the inhibition constant was observed at pH 7.0. For time-dependent inhibitors desethylamiodarone and paroxetine, changes in time-dependent inhibition parameters were different for each inhibitor. These results were incorporated into physiologically based pharmacokinetic models indicating that the changes in in vitro parameters determined at pH 7.0 offset the effect of increased unbound intracellular concentrations on apparent clearance and extent of drug-drug interactions. However, this offset between concentration and enzyme activity cannot be generalized for all substrates, inhibitors, and enzymes, as the effect of a lower pH in vitro varied significantly; therefore, it would be prudent to determine in vitro enzyme parameters at the hepatocyte-appropriate pH 7.0.

Prediction of in vivo clearance and associated variability of CYP2C19 substrates by genotypes in populations utilizing a pharmacogenetics-based mechanistic model.

It is important to examine the cytochrome P450 2C19 (CYP2C19) genetic contribution to drug disposition and responses of CYP2C19 substrates during drug development. Design of such clinical trials requires projection of genotype-dependent in vivo clearance and associated variabilities of the investigational drug, which is not generally available during early stages of drug development, but is essential for CYP2C19 substrates with multiple clearance pathways. This study evaluated the utility of pharmacogenetics-based mechanistic modeling in predicting such parameters. Hepatic CYP2C19 activity and variability within genotypes were derived from in vitro S-mephenytoin metabolic activity in genotyped human liver microsomes (N = 128). These data were then used in mechanistic models to predict genotype-dependent disposition of CYP2C19 substrates (i.e., S-mephenytoin, citalopram, pantoprazole, and voriconazole) by incorporating in vivo clearance or pharmacokinetics of wild-type subjects and parameters of other clearance pathways. Relative to the wild-type, the CYP2C19 abundance (coefficient of variation percentage) in CYP2C19*17/*17, *1/*17, *1/*1, *17/null, *1/null, and null/null microsomes was estimated as 1.85 (117%), 1.79 (155%), 1.00 (138%), 0.83 (80%), 0.38 (130%), and 0 (0%), respectively. The subsequent modeling and simulations predicted, within 2-fold of the observed, the means and variabilities of urinary S/R-mephenytoin ratio (36 of 37 genetic groups), the oral clearance of citalopram (9 of 9 genetic groups) and pantoprazole (6 of 6 genetic groups), and voriconazole oral clearance (4 of 4 genetic groups). Thus, relative CYP2C19 genotype-dependent hepatic activity and variability were quantified in vitro and used in a mechanistic model to predict pharmacokinetic variability, thus allowing the design of pharmacogenetics and drug-drug interaction trials for CYP2C19 substrates.

Prediction of renal transporter mediated drug-drug interactions for pemetrexed using physiologically based pharmacokinetic modeling.

Pemetrexed, an anionic anticancer drug with a narrow therapeutic index, is eliminated mainly by active renal tubular secretion. The in vitro to in vivo extrapolation approach used in this work was developed to predict possible drug-drug interactions (DDIs) that may occur after coadministration of pemetrexed and nonsteroidal anti-inflammatory drugs (NSAIDs), and it included in vitro assays, risk assessment models, and physiologically based pharmacokinetic (PBPK) models. The pemetrexed transport and its inhibition parameters by several NSAIDs were quantified using HEK-PEAK cells expressing organic anion transporter (OAT) 3 or OAT4. The NSAIDs were ranked according to their DDI index, calculated as the ratio of their maximum unbound concentration in plasma over the concentration inhibiting 50% (IC50) of active pemetrexed transport. A PBPK model for ibuprofen, the NSAID with the highest DDI index, was built incorporating active renal secretion in Simcyp Simulator. The bottom-up model for pemetrexed underpredicted the clearance by 2-fold. The model we built using a scaling factor of 5.3 for the maximal uptake rate (Vmax) of OAT3, which estimated using plasma concentration profiles from patients given a 10-minute infusion of 500 mg/m(2) of pemetrexed supplemented with folic acid and vitamin B12, recovered the clinical data adequately. The observed/predicted increases in Cmax and the area under the plasma-concentration time curve (AUC0-inf) of pemetrexed when ibuprofen was coadministered were 1.1 and 1.0, respectively. The coadministration of all other NSAIDs was predicted to have no significant impact on the AUC0-inf based on their DDI indexes. The PBPK model reasonably reproduced pemetrexed concentration time profiles in cancer patients and its interaction with ibuprofen.

Oscillatory beta activity mediates neuroplastic effects of motor cortex stimulation in humans.

Continuous theta burst stimulation (cTBS) is a repetitive transcranial magnetic stimulation protocol that can inhibit human motor cortex (M1) excitability and impair movement for ≤ 1 h. While offering valuable insights into brain function and potential therapeutic benefits, these neuroplastic effects are highly variable between individuals. The source of this variability, and the electrophysiological mechanisms underlying the inhibitory after-effects, are largely unknown. In this regard, oscillatory activity at beta frequency (15-35 Hz) is of particular interest as it is elevated in motor disorders such as Parkinson's disease and modulated during the generation of movements. Here, we used a source-level magnetoencephalography approach to investigate the hypothesis that the presence of neuroplastic effects following cTBS is associated with concurrent changes in oscillatory M1 beta activity. M1 cortices were localized with a synthetic aperture magnetometry beamforming analysis of visually cued index finger movements. Virtual electrode analysis was used to reconstruct the spontaneous and movement-related oscillatory activity in bilateral M1 cortices, before and from 10 to 45 min after cTBS. We demonstrate that 40 s of cTBS applied over left M1 reduced corticospinal excitability in the right index finger of 8/16 participants. In these responder participants only, cTBS increased the power of the spontaneous beta oscillations in stimulated M1 and delayed reaction times in the contralateral index finger. No further changes were observed in the latency or power of movement-related beta oscillations. These data provide insights into the electrophysiological mechanisms underlying cTBS-mediated impairment of motor function and demonstrate the association between spontaneous oscillatory beta activity in M1 and the inhibition of motor function.

Pharmacokinetic and pharmacodynamic alterations in the Roux-en-Y gastric bypass recipients.

OBJECTIVE: We conducted a pharmacokinetic (PK) study and a pharmacodynamic (PD) study to assess whether Roux-en-Y gastric bypass (RYGB) surgery is associated with significant changes to PK and PD of oral medications. BACKGROUND: The effect of RYGB on oral drug disposition is not well understood. METHODS: An oral cocktail of probe drugs for major drug-metabolizing enzymes (caffeine, tolbutamide, omeprazole, dextromethorphan, and oral and intravenous midazolam) was administered to 18 RYGB recipients and 18 controls. Timed blood and urine samples were obtained for PK analyses. Forty mg of oral furosemide was administered to 13 RYGB recipients and 14 controls, and urine and blood samples were collected for assessing furosemidePK, and urine volume and urine sodium excretion for PD analyses. RESULTS: Compared with controls, the RYGB group had significantly lower time to maximum plasma concentration (tmax) for caffeine (0.58 ± 0.5 vs 2.1 ± 2.2 hours, P < 0.0001), tolbutamide (1.4 ± 1.8 vs 2.1 ± 2.2 hours, P = 0.0001), omeprazole (1.1 ± 1.1 vs 4.4 ± 1.3 hours, P < 0.0001), and oral midazolam (0.5 ± 0.2 vs 0.7 ± 0.4 hours, P < 0.01). However, maximum plasma concentration, half-life, area under the curve, and oral bioavailability were not different. Compared with controls, the RYGB group had brisk natriuresis, with significantly lower tmax for urine sodium (1.3 ± 0.5 vs 3.1 ± 2.3 hours, P < 0.02) and correspondingly lower tmax for furosemide (1.8 ± 0.3 vs 4.2 ± 1.2 hours, P = 0.006). However, 6-hour urine sodium and 6-hour urine volume were not different between the two groups. CONCLUSIONS: RYGB recipients have significantly shorter tmax for the studied orally administered medications, but otherwise no other significant changes in PK were reported.

Optimization of drug-drug interaction study design: comparison of minimal physiologically based pharmacokinetic models on prediction of CYP3A inhibition by ketoconazole.

Ketoconazole is a potent CYP3A inhibitor used to assess the contribution of CYP3A to drug clearance and quantify the increase in drug exposure due to a strong inhibitor. Physiologically based pharmacokinetic (PBPK) models have been used to evaluate treatment regimens resulting in maximal CYP3A inhibition by ketoconazole but have reached different conclusions. We compare two PBPK models of the ketoconazole-midazolam interaction, model 1 (Chien et al., 2006) and model 2 implemented in Simcyp (version 11), to predict 16 published treatment regimens. With use of model 2, 41% of the study point estimates of area under the curve (AUC) ratio and 71% of the 90% confidence intervals were predicted within 1.5-fold of the observed, but these increased to 82 and 100%, respectively, with model 1. For midazolam, model 2 predicted a maximal midazolam AUC ratio of 8 and a hepatic fraction metabolized by CYP3A (f(m)) of 0.97, whereas model 1 predicted 17 and 0.90, respectively, which are more consistent with observed data. On the basis of model 1, ketoconazole (400 mg QD) for at least 3 days and substrate administration within 2 hours is required for maximal CYP3A inhibition. Ketoconazole treatment regimens that use 200 mg BID underestimate the systemic fraction metabolized by CYP3A (0.86 versus 0.90) for midazolam. The systematic underprediction also applies to CYP3A substrates with high bioavailability and long half-lives. The superior predictive performance of model 1 reflects the need for accumulation of ketoconazole at enzyme site and protracted inhibition. Model 2 is not recommended for inferring optimal study design and estimation of fraction metabolized by CYP3A.

CYP2B6 pharmacogenetics-based in vitro-in vivo extrapolation of efavirenz clearance by physiologically based pharmacokinetic modeling.

Efavirenz is mainly cleared by CYP2B6. The CYP2B6*6 allele is associated with lower efavirenz clearance. Efavirenz clearance was predictable using in vitro data for carriers of the CYP2B6*1/*1 genotype, but the prediction in carriers of the CYP2B6*6 allele was poor. To test the hypothesis that incorporation of mechanism of reduced efavirenz metabolism by the CYP2B6*6 allele can predict the genetic effect on efavirenz pharmacokinetics, in vitro-in vivo extrapolation of efavirenz clearance was performed by physiologically based pharmacokinetic modeling (Simcyp Simulator; Simcyp Ltd., Sheffield, UK) using data obtained from expressed CYP2B6.1 and CYP2B6.6 as well as human liver microsomes (HLMs) with CYP2B6*1/*1, *1/*6, and *6/*6 genotypes. Simulated pharmacokinetics of a single 600-mg oral dose of efavirenz for individuals with each genotype was compared with data observed in healthy subjects genotyped for the CYP2B6*6 allele (n = 20). Efavirenz clearance for carriers of the CYP2B6*1/*1 genotype was predicted reasonably well using HLM data, but the clearance in carriers of the CYP2B6*6 allele was underpredicted using both expressed and HLM systems. Improved prediction of efavirenz clearance was obtained from expressed CYP2B6 after recalculating intersystem extrapolation factors for CYP2B6.1 and CYP2B6.6 based on in vitro intrinsic clearance of bupropion 4-hydroxylation. These findings suggest that genetic effect on both CYP2B6 protein expression and catalytic efficiency needs to be taken into account for the prediction of pharmacokinetics in individuals carrying the CYP2B6*6/*6 genotype. Expressed CYP2B6 proteins may be a reliable in vitro system to predict effect of the CYP2B6*6 allele on the metabolism of CYP2B6 substrates.

The relative contributions of CYP3A4 and CYP3A5 to the metabolism of vinorelbine.

Vinorelbine is a semisynthetic vinca alkaloid used in the treatment of advanced breast and non-small cell lung cancers. Vincristine, a related vinca alkaloid, is 9-fold more efficiently metabolized by CYP3A5 than by CYP3A4 in vitro. This study quantified the relative contribution of CYP3A4 and CYP3A5 to the metabolism of vinorelbine in vitro using cDNA-expressed human cytochrome P450s (P450s) and human liver microsomes (HLMs). CYP3A4 and CYP3A5 were identified as the P450s capable of oxidizing vinorelbine using a panel of human enzymes and selective P450 inhibitors in HLMs. For CYP3A4 coexpressed with cytochrome b5 (CYP3A4+b5) and CYP3A5+b5, the Michaelis-Menten constants for vinorelbine were 2.6 and 3.6 µM, respectively, but the Vmax of 1.4 pmol/min/pmol was common to both enzymes. In HLMs, the intrinsic clearance of vinorelbine metabolism was highly correlated with CYP3A4 activity, and there was no significant difference in intrinsic clearance between CYP3A5 high and low expressers. When radiolabeled vinorelbine substrate was used, there were clear qualitative differences in metabolite formation fingerprints between CYP3A4+b5 and CYP3A5+b5 as determined by NMR and mass spectrometry analysis. One major metabolite (M2), a didehydro-vinorelbine, was present in both recombinant and microsomal systems but was more abundant in CYP3A4+b5 incubations. We conclude that despite the equivalent efficiency of recombinant CYP3A4 and CYP3A5 in vinorelbine metabolism the polymorphic expression of CYP3A5, as shown by the kinetics with HLMs, may have a minimal effect on systemic clearance of vinorelbine.

Pharmacogenomics of gemcitabine metabolism: functional analysis of genetic variants in cytidine deaminase and deoxycytidine kinase.

Gemcitabine (dFdC, 2',2'-difluorodeoxycytidine) is metabolized by cytidine deaminase (CDA) and deoxycytidine kinase (DCK), but the contribution of genetic variation in these enzymes to the variability in systemic exposure and response observed in cancer patients is unclear. Wild-type enzymes and variants of CDA (Lys27Gln and Ala70Thr) and DCK (Ile24Val, Ala119Gly, and Pro122Ser) were expressed in and purified from Escherichia coli, and enzyme kinetic parameters were estimated for cytarabine (Ara-C), dFdC, and its metabolite 2',2'-difluorodeoxyuridine (dFdU) as substrates. All three CDA proteins showed similar K(m) and V(max) for Ara-C and dFdC deamination, except for CDA70Thr, which had a 2.5-fold lower K(m) and 6-fold lower V(max) for Ara-C deamination. All four DCK proteins yielded comparable metabolic activity for Ara-C and dFdC monophosphorylation, except for DCK24Val, which demonstrated an approximately 2-fold increase (P < 0.05) in the intrinsic clearance of dFdC monophosphorylation due to a 40% decrease in K(m) (P < 0.05). DCK did not significantly contribute to dFdU monophosphorylation. In conclusion, the Lys27Gln substitution does not significantly modulate CDA activity toward dFdC, and therefore would not contribute to interindividual variability in response to gemcitabine. The higher in vitro catalytic efficiency of DCK24Val toward dFdC monophosphorylation may be relevant to dFdC clinical response. The substrate-dependent alterations in activities of CDA70Thr and DCK24Val in vitro were observed for the first time, and demonstrate that the in vivo consequences of these genetic variations should not be extrapolated from one substrate of these enzymes to another.

Literature based drug interaction prediction with clinical assessment using electronic medical records: novel myopathy associated drug interactions.

Drug-drug interactions (DDIs) are a common cause of adverse drug events. In this paper, we combined a literature discovery approach with analysis of a large electronic medical record database method to predict and evaluate novel DDIs. We predicted an initial set of 13197 potential DDIs based on substrates and inhibitors of cytochrome P450 (CYP) metabolism enzymes identified from published in vitro pharmacology experiments. Using a clinical repository of over 800,000 patients, we narrowed this theoretical set of DDIs to 3670 drug pairs actually taken by patients. Finally, we sought to identify novel combinations that synergistically increased the risk of myopathy. Five pairs were identified with their p-values less than 1E-06: loratadine and simvastatin (relative risk or RR = 1.69); loratadine and alprazolam (RR = 1.86); loratadine and duloxetine (RR = 1.94); loratadine and ropinirole (RR = 3.21); and promethazine and tegaserod (RR = 3.00). When taken together, each drug pair showed a significantly increased risk of myopathy when compared to the expected additive myopathy risk from taking either of the drugs alone. Based on additional literature data on in vitro drug metabolism and inhibition potency, loratadine and simvastatin and tegaserod and promethazine were predicted to have a strong DDI through the CYP3A4 and CYP2D6 enzymes, respectively. This new translational biomedical informatics approach supports not only detection of new clinically significant DDI signals, but also evaluation of their potential molecular mechanisms.

Quantitative structure-activity relationship models of clinical pharmacokinetics: clearance and volume of distribution.

Reliable prediction of two fundamental human pharmacokinetic (PK) parameters, systemic clearance (CL) and apparent volume of distribution (Vd), determine the size and frequency of drug dosing and are at the heart of drug discovery and development. Traditionally, estimated CL and Vd are derived from preclinical in vitro and in vivo absorption, distribution, metabolism, and excretion (ADME) measurements. In this paper, we report quantitative structure-activity relationship (QSAR) models for prediction of systemic CL and steady-state Vd (Vdss) from intravenous (iv) dosing in humans. These QSAR models avoid uncertainty associated with preclinical-to-clinical extrapolation and require two-dimensional structure drawing as the sole input. The clean, uniform training sets for these models were derived from the compilation published by Obach et al. (Drug Metab. Disp. 2008, 36, 1385-1405). Models for CL and Vdss were developed using both a support vector regression (SVR) method and a multiple linear regression (MLR) method. The SVR models employ a minimum of 2048-bit fingerprints developed in-house as structure quantifiers. The MLR models, on the other hand, are based on information-rich electro-topological states of two-atom fragments as descriptors and afford reverse QSAR (RQSAR) analysis to help model-guided, in silico modulation of structures for desired CL and Vdss. The capability of the models to predict iv CL and Vdss with acceptable accuracy was established by randomly splitting data into training and test sets. On average, for both CL and Vdss, 75% of test compounds were predicted within 2.5-fold of the value observed and 90% of test compounds were within 5.0-fold of the value observed. The performance of the final models developed from 525 compounds for CL and 569 compounds for Vdss was evaluated on an external set of 56 compounds. The predictions were either better or comparable to those predicted by other in silico models reported in the literature. To demonstrate the practical application of the RQSAR approach, the structure of vildagliptin, a high-CL and a high-Vdss compound, is modified based on the atomic contributions to its predicted CL and Vdss to propose compounds with lower CL and lower Vdss.

Rate of onset of inhibition of gut-wall and hepatic CYP3A by clarithromycin.

AIMS: To determine the extent and time-course of hepatic and intestinal cytochrome P450 3A (CYP3A) inactivation due to the mechanism-based inhibitor clarithromycin. METHODS: Intestinal and hepatic CYP3A inhibition was examined in 12 healthy volunteers following the administration of single and multiple doses of oral clarithromycin (500 mg). Intestinal biopsies were obtained under intravenous midazolam sedation at baseline and after the first dose, on days 2-4, and on days 6-8 of the clarithromycin treatment. The formation of 1'-hydroxymidazolam in biopsy tissue and the serum 1'-hydroxymidazolam:midazolam ratio were indicators of intestinal and hepatic CYP3A activity, respectively. RESULTS: Intestinal CYP3A activity decreased by 64 % (p = 0.0029) following the first dose of clarithromycin, but hepatic CYP3A activity did not significantly decrease. Repeated dosing of clarithromycin caused a significant decrease in hepatic CYP3A activity (p = 0.005), while intestinal activity showed little further decline. The CYP3A5 or CYP3A4*1B genotype were unable to account for inter-individual variability in CYP3A activity. CONCLUSIONS: Following the administration of clarithromycin, the onset of hepatic CYP3A inactivation is delayed compared to that of intestinal CYP3A. The time-course of drug-drug interactions due to clarithromycin will vary with the relative contribution of intestinal and hepatic CYP3A to the clearance and bioavailability of a victim substrate.

Sustained neuronal activity generated by glial plasticity.

Astrocytes release gliotransmitters, notably glutamate, that can affect neuronal and synaptic activity. In particular, astrocytic glutamate release results in the generation of NMDA receptor (NMDA-R)-mediated slow inward currents (SICs) in neurons. However, factors underlying the emergence of SICs and their physiological roles are essentially unknown. Here we show that, in acute slices of rat somatosensory thalamus, stimulation of lemniscal or cortical afferents results in a sustained increase of SICs in thalamocortical (TC) neurons that outlasts the duration of the stimulus by 1 h. This long-term enhancement of astrocytic glutamate release is induced by group I metabotropic glutamate receptors and is dependent on astrocytic intracellular calcium. Neuronal SICs are mediated by extrasynaptic NR2B subunit-containing NMDA-Rs and are capable of eliciting bursts. These are distinct from T-type Ca(2+) channel-dependent bursts of action potentials and are synchronized in neighboring TC neurons. These findings describe a previously unrecognized form of excitatory, nonsynaptic plasticity in the CNS that feeds forward to generate local neuronal firing long after stimulus termination.