Development and applications of a physiologically-based model of paediatric oral drug absorption.

There is increasing interest in paediatric drug absorption and the development of biopharmaceutics tools to facilitate the development of oral formulations for neonates, infants and children. We describe the development and application of a physiologically-based model of paediatric drug absorption applicable from full term birth onwards. Paediatric age-specific parameters were included for salivary flow, gastric pH, gastric emptying (and associated food effects) and duodenal bile salt concentrations and the associated algorithms were integrated into a dissolution, absorption and metabolism model as part of a PBPK platform. For other parameters, there was either evidence for no age-related changes or a lack of data, so that adult values were applied. An initial assessment of the model was carried out by simulating the oral absorption of theophylline, paracetamol and ketoconazole over a range of paediatric ages. The absorption of the first two drugs, both BCS class 1 compounds, was predicted to be slower in early neonates compared to older age groups (median tmax values of 3 vs 2h, respectively), but with invariant fraction absorbed (fa). This is in agreement with clinical observations. The tmax of ketoconazole, a BCS class 2 compound, was predicted to be about 1h in both neonates and adults, but the fa value was higher in the former (0.87 vs 0.69). There is clearly a need to expand the components of the model as new information on the ontogeny of GI tract parameters becomes available, and to assess it against more in vivo data with evidence of specific age-related changes in oral drug absorption.

Metformin and cimetidine: Physiologically based pharmacokinetic modelling to investigate transporter mediated drug-drug interactions.

Metformin is used as a probe for OCT2 mediated transport when investigating possible DDIs with new chemical entities. The aim of the current study was to investigate the ability of physiologically-based pharmacokinetic (PBPK) models to simulate the effects of OCT and MATE inhibition by cimetidine on metformin kinetics. PBPK models were developed, incorporating mechanistic kidney and liver sub-models for metformin (OCT and MATE substrate) and a mechanistic kidney sub-model for cimetidine. The models were used to simulate inhibition of the MATE1, MATE2-K, OCT1 and OCT2 mediated transport of metformin by cimetidine. Assuming competitive inhibition and using cimetidine Ki values determined in vitro, the predicted metformin AUC ratio was 1.0 compared to an observed value of 1.46. The observed AUC ratio could only be recovered with this model when the cimetidine Ki for OCT2 was decreased 1000-fold or the Ki's for both OCT1 and OCT2 were decreased 500-fold. An alternative description of metformin renal transport by OCT1 and OCT2, incorporating electrochemical modulation of the rate of metformin uptake together with 8-18-fold decreases in cimetidine Ki's for OCTs and MATEs, allowed recovery of the extent of the observed effect of cimetidine on metformin AUC. While the final PBPK model has limitations, it demonstrates the benefit of allowing for the complexities of passive permeability combined with active cellular uptake modulated by an electrochemical gradient and active efflux.

Impact of the pharmaceutical sciences on health care: a reflection over the past 50 years.

During the last century, particularly the latter half, spectacular progress has been made in improving the health and longevity of people. The reasons are many, but the development of medicines has played a critical role. This report documents and reflects on the impressive contribution that those working in the pharmaceutical sciences have made to healthcare over the past 50 years. It is divided into six sections (drug discovery; absorption, distribution, metabolism, and excretion; pharmacokinetics and pharmacodynamics; drug formulation; drug regulation; and drug utilization), each describing key contributions that have been made in the progression of medicines, from conception to use. A common thread throughout is the application of translational science to the improvement of drug discovery, development, and therapeutic application. Each section has been coordinated by a leading scientist who was asked, after consulting widely with many colleagues across the globe, to identify "The five most influential ideas/concepts/developments introduced by 'pharmaceutical scientists' (in their field) over the past 50 years?" Although one cannot predict where the important breakthroughs will come in the future to meet the unmet medical needs, the evidence presented in this report should leave no doubt that those engaged in the pharmaceutical sciences will continue to make their contributions heavily felt.

Prediction of time-dependent CYP3A4 drug-drug interactions by physiologically based pharmacokinetic modelling: impact of inactivation parameters and enzyme turnover.

Predicting the magnitude of time-dependent metabolic drug-drug (mDDIs) interactions involving cytochrome P-450 3A4 (CYP3A4) from in vitro data requires accurate knowledge of the inactivation parameters of the inhibitor (K(I), k(inact)) and of the turnover of the enzyme (k(deg)) in both the gut and the liver. We have predicted the magnitude of mDDIs observed in 29 in vivo studies involving six CYP3A4 probe substrates and five mechanism based inhibitors of CYP3A4 of variable potency (azithromycin, clarithromycin, diltiazem, erythromycin and verapamil). Inactivation parameters determined anew in a single laboratory under standardised conditions together with data from substrate and inhibitor files within the Simcyp Simulator (Version 9.3) were used to determine a value of the hepatic k(deg) (0.0193 or 0.0077h(-1)) most appropriate for the prediction of mDDIs involving time-dependent inhibition of CYP3A4. The higher value resulted in decreased bias (geometric mean fold error - 1.05 versus 1.30) and increased precision (root mean squared error - 1.29 versus 2.30) of predictions of mean ratios of AUC in the absence and presence of inhibitor. Depending on the k(deg) value used (0.0193 versus 0.0077h(-1)), predicted mean ratios of AUC were within 2-fold of the observed values for all (100%) and 27 (93%) of the 29 studies, respectively and within 1.5-fold for 24 (83%) and 17 (59%) of the 29 studies, respectively. Comprehensive PBPK models were applied for accurate assessment of the potential for mDDIs involving time-dependent inhibition of CYP3A4 using a hepatic k(deg) value of 0.0193h(-1) in conjunction with inactivation parameters determined by the conventional experimental approach.

Towards a quantitative framework for the prediction of DDIs arising from cytochrome P450 induction.

Although CYP induction is not generally considered to be as clinically relevant as CYP inhibition, there are important examples where induction has caused both therapeutic failure, due to insufficient exposure to parent drug, and toxicity, mediated by increased formation of reactive metabolites. Furthermore, while there has been considerable progress in the extrapolation of in vitro data to predict the in vivo consequences of enzyme inhibition, less attention has been given to the quantitative impact of enzyme induction as a mechanism of drug-drug interaction (DDI) and as a component of compound selection and early drug development. We discuss current approaches in the context of a mechanistic framework for the prediction of the extent and time-course of enzyme induction in vivo based on in vitro experimentation. Factors influencing the extent of DDI due to CYP induction are summarised, and areas deficient in information that would allow more accurate prediction within target populations are highlighted.

Modified-release hydrocortisone for circadian therapy: a proof-of-principle study in dexamethasone-suppressed normal volunteers.

BACKGROUND: All existing long-term glucocorticoid replacement therapy is suboptimal as the normal nocturnal rise and waking morning peak of serum cortisol is not reproduced. AIM: To test whether it is possible to reproduce the normal overnight rise and morning peak in serum cortisol using an oral delayed and sustained release preparation of hydrocortisone (Cortisol(ds)). SUBJECTS AND METHODS: Six healthy normal male volunteers attended on two occasions, in a single-dose, open-label, nonrandomized study. Endogenous cortisol secretion was suppressed by administration of dexamethasone. Cortisol(ds) (formulation A or B) was administered at 2200 h on day 1. Blood samples for measurement of cortisol were taken from 2200 h every 30 min until 0700 h, then hourly until 2200 h on day 2. Fifteen body mass index (BMI)-matched control subjects had serum cortisol levels measured at 20-min intervals for 24 h. Serum cortisol profiles and pharmacokinetics after Cortisol(ds) were compared with those in controls. RESULTS: Formulations A and B were associated with delayed drug release (by 2 h and 4 h, respectively), with median peak cortisol concentrations at 4.5 h (0245 h) and 10 h (0800 h), respectively, thereby reproducing the normal early morning rise in serum cortisol. Total cortisol exposure was not different from controls. CONCLUSIONS: For the first time we have shown that it is possible to mimic the normal circadian rhythm of circulating cortisol with an oral modified-release formulation of hydrocortisone, providing the basis for development of physiological circadian replacement therapy in patients with adrenal insufficiency.

Disparity in holoprotein/apoprotein ratios of different standards used for immunoquantification of hepatic cytochrome P450 enzymes.

An analysis of reported hepatic abundances of CYP3A4 and 3A5 indicated that values determined by immunoquantification using commercially available, unpurified recombinant enzymes as standards are significantly lower than those determined using purified enzymes or human liver microsomes characterized with lysosomal peptides (CYP3A4: mean 45 versus 121 pmol/mg protein, p < 0.01; CYP3A5: mean 28 versus 83 pmol/mg protein, p < 0.05). When immunoquantifying cytochromes P450 (P450s), it is assumed that the holoprotein (holo)/apoprotein ratio is the same in the samples and the standard. Estimates of holo/apoprotein ratios from data reported for a range of P450s purified from human liver and non-commercial recombinant systems indicated less than complete and variable heme coupling dependent on enzyme and system.

Prediction of in vivo drug clearance from in vitro data. II: potential inter-ethnic differences.

Potential differences in drug clearance between Japanese and Caucasians were investigated by integrating data on demography, liver size, the abundance of the major cytochromes P450 and in vitro metabolic parameters. Eleven drugs (alprazolam, caffeine, chlorzoxazone, cyclosporine, midazolam, omeprazole, sildenafil, tolbutamide, triazolam, S-warfarin and zolpidem) fulfilled the entry criteria of the study (i.e. the necessary in vitro metabolism data were available and clearance values had been reported both in Caucasians and Japanese). Values of relevant biological variables were obtained from the literature, and clearance predictions were made using the Simcyp Population-Based ADME Simulator. The ratios of observed oral clearance (CLp.o.) values in Caucasians compared with Japanese ranged from 0.6 to 2.8 (integrating data from 82 sources). The CLp.o. values for alprazolam, caffeine and zolpidem were not statistically different between Caucasian and Japanese (p>0.05), whereas those for chorzoxazone, cyclosporine, omeprazole, tolbutamide and triazolam were higher in Caucasians (p<0.05), and those for midazolam, sildenafil and S-warfarin were higher in Japanese (p<0.05). CLp.o. values, predicted from in vitro data, were within 3-fold of observed in vivo values for seven of the 11 drugs in Japanese. Values for the predicted ratios ranged from 1.6 to 4.9. The predicted ratios were not significantly different from observed ratios for cyclosporine, omeprazole, tolbutamide and triazolam. Only partial success in predicting ethnic differences in clearance indicates the need for larger and more reliable databases on relevant variables. With such information, in silico predictions might be used with more confidence to decrease the need for repeating pharmacokinetic studies in different ethnic groups.

Circadian hydrocortisone infusions in patients with adrenal insufficiency and congenital adrenal hyperplasia.

OBJECTIVE: Conventional hydrocortisone therapy in adrenal insufficiency cannot provide physiological replacement. We have explored the potential of circadian delivery of hydrocortisone as proof of concept for such therapy delivered in modified-release tablet formulation. METHODS: We investigated whether the circadian intravenous infusion of hydrocortisone could improve control of ACTH and androgen levels. Two healthy subjects, two patients with Addison's disease and two patients with congenital adrenal hyperplasia (CAH) were studied. RESULTS: In patients on thrice daily oral hydrocortisone, peak serum cortisol levels were higher than in normal subjects and overnight levels were very low. Patients had very high plasma ACTH levels before their morning dose of hydrocortisone, both at the beginning and at the end of their conventional oral therapy: mean +/- SEM 311.8 +/- 123.2 and 311.2 +/- 85.4 ng/l, respectively. In the patients with CAH, serum 17-hydroxyprogesterone levels were also elevated: 550 and 642 nmol/l at the beginning and 550 and 777 nmol/l at the end of conventional treatment, respectively. The overall 24-h mean cortisol levels were similar for conventional oral hydrocortisone and the circadian infusion. At 0700 h, ACTH levels were much higher on conventional treatment than after circadian infusion: mean +/- SEM 311.2 +/- 85.4 vs. 70.5 +/- 45.0 ng/l, respectively (P < 0.05). The same pattern was observed in 17-hydroxyprogesterone levels, which were 550 and 777 nmol/l after conventional treatment and 3 and 64 nmol/l after circadian infusion. CONCLUSIONS: In patients with poor biochemical control of Addison's disease and CAH, a 24-h circadian infusion of hydrocortisone can decrease morning ACTH and 17-hydroxyprogesterone levels to near normal.

Prediction of metabolic drug clearance in humans: in vitro-in vivo extrapolation vs allometric scaling.

Previously in vitro-in vivo extrapolation (IVIVE) with the Simcyp Clearance and Interaction Simulator has been used to predict the clearance of 15 clinically used drugs in humans. The criteria for the selection of the drugs were that they are used as probes for the activity of specific cytochromes P450 (CYPs) or have a single CYP isoform as the major or sole contributor to their metabolism and that they do not exhibit non-linear kinetics in vivo. Where data were available for the clearance of the drugs in at least three animal species, the predictions from IVIVE have now been compared with those based on allometric scaling (AS). Adequate data were available for estimating oral clearance (CLp.o.) in 9 cases (alprazolam, sildenafil, caffeine, clozapine, cyclosporine, dextromethorphan, midazolam, omeprazole and tolbutamide) and intravenous clearance in 6 cases (CLi.v.) (cyclosporine, diclofenac, midazolam, omeprazole, theophylline and tolterodine). AS predictions were based on five different methods: (1) simple allometry (clearance versus body weight); (2) correction for maximum life-span potential (CL x MLP); (3) correction for brain weight (CL x BrW); (4) the use of body surface area; and (5) the rule of exponents. A prediction accuracy was indicated by mean-fold error and the Pearson product moment correlation coefficient. Predictions were considered successful if the mean-fold error was

Prediction of in vivo drug clearance from in vitro data. I: impact of inter-individual variability.

The Simcyp Population-Based ADME Simulator was used to predict median drug clearances and their associated variance from in vitro data. Fifteen drugs satisfied the entry criteria for the study and the relevant information (in vitro metabolism data and in vivo human clearance values) were collated from the literature. Predicted values of median clearances fell within 2-fold of observed values for 73% of the drugs (oral route) and 78% of the drugs (intravenous route) when microsomal binding was disregarded, and for 93% (oral) and 100% (intravenous) when it was considered. Irrespective of whether microsomal binding was considered, the predicted fold variability fell within 2-fold of the observed variability for 80% (oral) and 67% (intravenous) of the drugs.

A critical evaluation of the experimental design of studies of mechanism based enzyme inhibition, with implications for in vitro-in vivo extrapolation.

The published literature on mechanism based inhibition (MBI) of CYPs was evaluated with respect to experimental design, methodology and data analysis. Significant variation was apparent in the dilution factor, ratio of preincubation to incubation times and probe substrate concentrations used, and there were some anomalies in the estimation of associated kinetic parameters (k(inact), K(I), r). The impact of the application of inaccurate values of k(inact) and K(I) when extrapolating to the extent of inhibition in vivo is likely to be greatest for those compounds of intermediate inhibitory potency, but this also depends on the fraction of the net clearance of substrate subject to MBI and the pre-systemic and systemic exposure to the inhibitor. For potent inhibitors, the experimental procedure is unlikely to have a material influence on the maximum inhibition. Nevertheless, the bias in the values of the kinetic parameters may influence the time for recovery of enzyme activity following re-synthesis of the enzyme. Careful attention to the design of in vitro experiments to obtain accurate kinetic parameters is necessary for a reliable prediction of different aspects of the in vivo consequences of MBI. The review calls for experimental studies to quantify the impact of study design in studies of MBI, with a view to better harmonisation of protocols.

Mechanism-based inactivation of CYP2D6 by methylenedioxymethamphetamine.

The potency of methylenedioxymethamphetamine (MDMA) as a mechanism-based inhibitor of CYP2D6 has been defined using microsomes prepared from yeast expressing the enzyme and from three human livers. The inhibitory effect was increased by preincubation through formation of a metabolic intermediate complex. Inactivation parameters (kinact and KI), defined with respect to the O-demethylation of dextromethorphan, were 0.29 +/- 0.03 (S.E.) min(-1) and 12.9 +/- 3.6 (S.E.) microM for yeast-expressed CYP2D6, and 0.26 +/- 0.02 min(-1) and 14.4 +/- 2.5 microM, 0.15 +/- 0.01 min(-1) and 8.8 +/- 2.6 microM, and 0.12 +/- 0.05 min(-1) and 45.3 +/- 32.1 microM for the liver microsomal preparations. The rate of inactivation of CYP2D6 by MDMA decreased when quinidine, a competitive inhibitor of CYP2D6, was added to the primary incubation mixture. However, inactivation was unaffected by the addition of glutathione. The results indicate that MDMA is a potent mechanism-based inhibitor of CYP2D6, with implications for understanding its in vivo disposition and drug interaction potential.

Predicting drug clearance from recombinantly expressed CYPs: intersystem extrapolation factors.

1. Recombinantly expressed human cytochromes P450 (rhCYPs) have been underused for the prediction of human drug clearance (CL). 2. Differences in intrinsic activity (per unit CYP) between rhCYP and human liver enzymes complicate the issue and these discrepancies have not been investigated systematically. We define intersystem extrapolation factors (ISEFs) that allow the use of rhCYP data for the in vitro-in vivo extrapolation of human drug CL and the variance that is associated with interindividual variation of CYP abundance due to genetic and environmental effects. 3. A large database (n = 451) of metabolic stability data has been compiled and used to derive ISEFs for the most commonly used expression systems and CYP enzymes. 4. Statistical models were constructed for the ISEFs to determine major covariates in order to optimize experimental design to increase prediction accuracy. 5. Suggestions have been made for the conduct of future studies using rhCYP to predict human drug clearance.

Inter-individual variability in levels of human microsomal protein and hepatocellularity per gram of liver.

AIMS: To determine levels of microsomal protein (MPPGL) and hepatocellularity (HPGL) per gram of human liver and their interindividual variability. METHODS: Triplicate liver samples were used to determine values of MPPGL (n = 20) and HPGL (n = 7) after accounting for the fractional loss of microsomal protein or hepatocytes during processing. Repeated measurements from each liver sample allowed the estimation of true interindividual variability in MPPGL and HPGL using ANOVA. RESULTS: The value of MPPGL ranged from 26 to 54 mg g(-1) (mean(geo)= 33 mg g(-1)). The value of HPGL ranged from 65 to 185 x 10(6) cells g(-1) (mean(geo)= 10(7) x 10(6) cells g(-1)). CONCLUSIONS: There is significant interindividual variability in MPPGL, which has implications for the accurate extrapolation of in vitro data on drug metabolism to predict in vivo metabolic clearance.

A discordance between cytochrome P450 2D6 genotype and phenotype in patients undergoing methadone maintenance treatment.

AIMS: To assess CYP2D6 activity and genotype in a group of patients undergoing methadone maintenance treatment (MMT). METHODS: Blood samples from 34 MMT patients were genotyped by a polymerase chain reaction-based method, and results were compared with CYP2D6 phenotype (n = 28), as measured by the molar metabolic ratio (MR) of dextromethorphan (DEX)/dextrorphan (DOR) in plasma. RESULTS: Whereas 9% of patients (3/34) were poor metabolizers (PM) by genotype, 57% (16/28) were PM by phenotype (P < 0.005). Eight patients, who were genotypically extensive metabolizers (EM), were assigned as PM by their phenotype. The number of CYP2D6*4 alleles and sex were significant determinants of CYP2D6 activity in MMT patients, whereas other covariates (methadone dose, age, weight) did not contribute to variation in CYP2D6 activity. CONCLUSIONS: There was a discordance between genotype and in vivo CYP2D6 activity in MMT patients. This finding is consistent with inhibition of CYP2D6 activity by methadone and may have implications for the safety and efficacy of other CYP2D6 substrates taken by MMT patients.

Physiologically based modelling of inhibition of metabolism and assessment of the relative potency of drug and metabolite: dextromethorphan vs. dextrorphan using quinidine inhibition.

AIMS: To define the relative antitussive effect of dextromethorphan (DEX) and its primary metabolite dextrorphan (DOR) after administration of DEX. METHODS: Data were analysed from a double-blind, randomized cross-over study in which 22 subjects received the following oral treatments: (i) placebo; (ii) 30 mg DEX hydro-bromide; (iii) 60 mg DEX hydro-bromide; and (iv) 30 mg DEX hydro-bromide preceded at 1 h by quinidine HCl (50 mg). Cough was elicited using citric acid challenge. Pharmacokinetic data from all non-placebo arms of the study were fitted simultaneously. The parameters were then used as covariates in a link PK-PD model of cough suppression using data from all treatment arms. RESULTS: The best-fit PK model assumed two- and one-compartment PK models for DEX and DOR, respectively, and competitive inhibition of DEX metabolism by quinidine. The intrinsic clearance of DEX estimated from the model ranged from 59 to 1536 l x h(-1), which overlapped with that extrapolated from in vitro data (12-261 l x h(-1)) and showed similar variation (26- vs. 21-fold, respectively). The inhibitory effect of quinidine ([I]/Ki) was 19 (95% confidence interval of mean: 18-20) with an estimated average Ki of 0.017 microM. Although DEX and DOR were both active, the potency of the antitussive effect of DOR was 38% that of DEX. A sustained antitussive effect was related to slow removal of DEX/DOR from the effect site (ke0 = 0.07 h(-1)). CONCLUSIONS: Physiologically based PK modelling with perturbation of metabolism using an inhibitor allowed evaluation of the antitussive potency of DOR without the need for separate administration of DOR.

Determination of midazolam and 1'-hydroxymidazolam by liquid chromatography-mass spectrometry in plasma of patients undergoing methadone maintenance treatment.

A rapid, sensitive and selective LC-MS method was developed for the simultaneous determination of midazolam (MDZ) and 1'-hydroxymidazolam (1'-OHMDZ) in plasma taken from 54 patients undergoing methadone maintenance therapy, most of whom were multidrug users. Samples spiked with prazepam, the internal standard, and were extracted into diethyl ether. Compounds were separated on a Phenomenex Luna C(18) column and a mobile phase of acetonitrile-ammonium acetate buffer (10 mM, pH 4.7) (52:48, v/v) at a flow-rate of 1 ml/min. The limit of detection was 0.65 and 0.68 (ng/ml) for MDZ and 1'-OHMDZ, respectively. Within-day relative standard deviations were less than 8%.