Scaling Drug Clearance from Adults to the Young Children for Drugs Undergoing Hepatic Metabolism: A Simulation Study to Search for the Simplest Scaling Method

Previous research showed that scaling drug clearance from adults to children based on body weight alone is not accurate for all hepatically cleared drugs in very young children. This study systematically assesses the accuracy of scaling methods that, in addition to body weight, also take age-based variables into account for drugs undergoing hepatic metabolism in children younger than five years, namely scaling with (1) a body weight-based function using an age-dependent exponent (ADE) and (2) a body weight-based function with fixed exponent of 0.75 (AS0.75) combined with isoenzyme maturation functions (MF_PBPK) similar to those implemented in physiologically based pharmacokinetic (PBPK) models (AS0.75 + MF_PBPK). A PBPK-based simulation workflow was used, including hypothetical drugs with a wide range of properties and metabolized by different isoenzymes. Adult clearance values were scaled to seven typical children between one day and four years. Prediction errors of ± 50% were considered reasonably accurate. Isoenzyme maturation was found to be an important driver of changes in hepatic metabolic clearance in children younger than five years, which prevents the systematic accuracy of ADE scaling. AS0.75 + MF_PBPK, when accounting for maturation of isoenzymes and microsomal protein per gram of liver (MPPGL), can reasonably accurately scale hepatic metabolic clearance for all low and intermediate extraction ratio drugs except for drugs binding to alpha-1-acid glycoprotein in neonates. As differences in the impact of changes in system-specific parameters on drugs with different properties yield differences in clearance ontogeny, it is unlikely that for the remaining drugs, scaling methods that do not take drug properties into account will be systematically accurate.

Application of the MechPeff model to predict passive effective intestinal permeability in the different regions of the rodent small intestine and colon

A major component of physiologically based pharmacokinetic (PBPK) models is the prediction of the rate and extent of absorption of orally dosed drugs for which knowledge of effective passive intestinal permeability (P_eff ) is essential. Single-pass intestinal perfusion (SPIP) studies are used to establish effective permeability in vivo but are difficult to perform in rodents, while mechanistic models to predict drug P_eff in rat and mouse have not been published. This work evaluates the predictive performance of the 'MechPeff' model to predict P_eff in the rodent intestine based upon knowledge of regional gut physiology and drug-specific physicochemical parameters. The 'MechPeff' model, built-in to the Simcyp Rat and Mouse Simulators, predicts transcellular, paracellular and mucus layer permeabilities and combines these to give the overall P_eff . The jejunal and/or ileal P_eff was predicted for 12 (4) acidic, 13 (12) basic, 10 (8) neutral and 2 (0) ampholytic drugs in the rat (mouse), spanning a wide range of MW and logP_o:w , and compared with experimental P_eff obtained using SPIP. A key input is the intrinsic transcellular permeability (P_trans,0 ) which can be derived from modelling of appropriate in vitro permeability experiments or predicted from physicochemical properties. The P_eff predictions were reasonably good when experimentally derived P_trans,0 was used; from 42 P_eff,rat values, 24 (57%) were within 3-fold, and of 19 P_eff,mouse values, 12 (63%) were within 3-fold, of observed P_eff . Considering the lack of alternative models to predict P_eff in preclinical species, and the minimal drug-specific inputs required, this model provides a valuable tool within drug discovery and development programmes. Copyright © 2017 John Wiley & Sons, Ltd.

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.

Systematic and quantitative assessment of the effect of chronic kidney disease on CYP2D6 and CYP3A4/5

Recent reviews suggest that chronic kidney disease (CKD) can affect the pharmacokinetics of nonrenally eliminated drugs, but the impact of CKD on individual elimination pathways has not been systematically evaluated. In this study we developed a comprehensive dataset of the effect of CKD on the pharmacokinetics of CYP2D6‐ and CYP3A4/5‐metabolized drugs. Drugs for evaluation were selected based on clinical drug–drug interaction (CYP3A4/5 and CYP2D6) and pharmacogenetic (CYP2D6) studies. Information from dedicated CKD studies was available for 13 and 18 of the CYP2D6 and CYP3A4/5 model drugs, respectively. Analysis of these data suggested that CYP2D6‐mediated clearance is generally decreased in parallel with the severity of CKD. There was no apparent relationship between the severity of CKD and CYP3A4/5‐mediated clearance. The observed elimination‐route dependency in CKD effects between CYP2D6 and CYP3A4/5 may inform the need to conduct clinical CKD studies with nonrenally eliminated drugs for optimal use of drugs in patients with CKD.

Semiphysiologically based pharmacokinetic model for midazolam and CYP3A mediated metabolite 1-OH-midazolam in morbidly obese and weight loss surgery patients

This study aimed to describe the pharmacokinetics of midazolam and its cytochrome P450 3A (CYP3A) mediated metabolite 1‐OH‐midazolam in morbidly obese patients receiving oral and i.v. midazolam before (n = 20) and one year after weight loss surgery (n = 18), thereby providing insight into the influence of weight loss surgery on CYP3A activity in the gut wall and liver. In a semiphysiologically based pharmacokinetic (semi‐PBPK) model in which different blood flow scenarios were evaluated, intrinsic hepatic clearance of midazolam (CL_int,H) was 2 (95% CI 1.40–1.64) times higher compared to morbidly obese patients before surgery (P < 0.01). Midazolam gut wall clearance (CL_int,G) was slightly lower in patients after surgery (P > 0.05), with low values for both groups. The results of the semi‐PBPK model suggest that, in patients after weight loss surgery, CYP3A hepatic metabolizing capacity seems to recover compared to morbidly obese patients, whereas CYP3A mediated CL_int,G was low for both populations and showed large interindividual variability.

Physiologically Based Pharmacokinetics Is Impacting Drug Development and Regulatory Decision Making

It is no coincidence that the reports of two meetings, one organized by the US Food and Drug Administration (FDA), in March 2014, and the other by the UK Medicines and Healthcare Products Regulatory (MHRA), in collaboration with ABPI (the Association of British Pharmaceutical Industry), in June 2014, have been published in tandem in CPT-PSP.12 Both reports deal with the same topic, namely, the impact of physiologically based pharmacokinetics (PBPK) in clinical drug development and the best practices for such applications. This reflects the transition of PBPK from academic curiosity to industrial norm, manifested by the regulatory agencies encouraging its use and receiving an increasing number of submissions containing PBPK models. The goal of both meetings was to help determine the need and facilitate the development of regulatory guidances on this subject within the conceptual framework of model informed drug development and regulatory decision-making. A further reflection of this intent is the publication by the European Medicines Agency of a Concept Paper on PBPK.3 One is reminded of a similar train of events surrounding the introduction of population PK/PD and nonlinear mixed effects modeling in the early-late 1990s, again with encouragement and receptivity of regulatory agencies leading to FDA guidance on the topic.4 Indeed, the intention of PBPK modeling and simulation is to complement other approaches, such as compartmental modeling, or, in some cases, replace them with a more mechanistic approach. PBPK models represent an important class of models that characterize absorption, distribution, metabolism, excretion (ADME) processes and their underlying biological and physiological drivers. An increased understanding of these drivers and their unique interactions with drug substance and formulation factors provides critical insights into how drugs will behave in healthy volunteers and patients with disease.

Application of an LC-MS/MS method for the simultaneous quantification of human intestinal transporter proteins absolute abundance using a QconCAT technique

Transporter proteins expressed in the gastrointestinal tract play a major role in the oral absorption of some drugs, and their involvement may lead to drug-drug interaction (DDI) susceptibility when given in combination with drugs known to inhibit gut wall transporters. Anticipating such liabilities and predicting the magnitude of the impact of transporter proteins on oral drug absorption and DDIs requires quantification of their expression in human intestine, and linking these to data obtained through in vitro experiments. A quantitative targeted absolute proteomic method employing liquid chromatography coupled to tandem mass spectrometry (LC-MS/MS) together with a quantitative concatenation (QconCAT) strategy to provide proteotypic peptide standards has been applied to quantify ATP1A1 (sodium/potassium-ATPase; Na/K-ATPase), CDH17 (human peptide transporter 1; HPT1), ABCB1 (P-glycoprotein; P-gp), ABCG2 (breast cancer resistance protein; BCRP), ABCC2 (multidrug resistance-associated protein 2; MRP2) and SLC51A (Organic Solute Transporter subunit alpha; OST-α), in human distal jejunum (n=3) and distal ileum (n=1) enterocyte membranes. Previously developed selected reaction monitoring (SRM) schedules were optimised to enable quantification of the proteotypic peptides for each transporter. After harvesting enterocytes by calcium chelation elution and generating a total membrane fraction, the proteins were subjected to proteolytic digestion. To account for losses of peptides during the digestion procedure, a gravimetric method is also presented. The linearity of quantifying the QconCAT from an internal standard (correlation coefficient, R(2)=0.998) and quantification of all target peptides in a pooled intestinal quality control sample (R(2)≥ 0.980) was established. The assay was also assessed for within and between-day precision, demonstrating a <15% coefficient of variation for all peptides across 3 separate analytical runs, over 2 days. The methods were applied to obtain the absolute abundances for all targeted proteins. In all samples, Na/K-ATPase, HPT1, P-gp and BCRP were detected above the lower limit of quantitation (i.e., >0.2 fmol/μg membrane protein). MRP2 abundance could be quantified in distal jejunum but not in the distal ileum sample. OST-α was not detected in 2 out of 3 jejunum samples. This study highlights the utility of a QconCAT strategy to quantify absolute transporter abundances in human intestinal tissues.

A mechanistic framework for in vitro-in vivo extrapolation of liver membrane transporters: prediction of drug-drug interaction between rosuvastatin and cyclosporine

Background and Objectives

The interplay between liver metabolising enzymes and transporters is a complex process involving system-related parameters such as liver blood perfusion as well as drug attributes including protein and lipid binding, ionisation, relative magnitude of passive and active permeation. Metabolism- and/or transporter-mediated drug–drug interactions (mDDIs and tDDIs) add to the complexity of this interplay. Thus, gaining meaningful insight into the impact of each element on the disposition of a drug and accurately predicting drug–drug interactions becomes very challenging. To address this, an in vitro–in vivo extrapolation (IVIVE)-linked mechanistic physiologically based pharmacokinetic (PBPK) framework for modelling liver transporters and their interplay with liver metabolising enzymes has been developed and implemented within the Simcyp Simulator^®.

Methods

In this article an IVIVE technique for liver transporters is described and a full-body PBPK model is developed. Passive and active (saturable) transport at both liver sinusoidal and canalicular membranes are accounted for and the impact of binding and ionisation processes is considered. The model also accommodates tDDIs involving inhibition of multiple transporters. Integrating prior in vitro information on the metabolism and transporter kinetics of rosuvastatin (organic-anion transporting polypeptides OATP1B1, OAT1B3 and OATP2B1, sodium-dependent taurocholate co-transporting polypeptide [NTCP] and breast cancer resistance protein [BCRP]) with one clinical dataset, the PBPK model was used to simulate the drug disposition of rosuvastatin for 11 reported studies that had not been used for development of the rosuvastatin model.

Results

The simulated area under the plasma concentration–time curve (AUC), maximum concentration (C_max) and the time to reach C_max (t_max) values of rosuvastatin over the dose range of 10–80 mg, were within 2-fold of the observed data. Subsequently, the validated model was used to investigate the impact of coadministration of cyclosporine (ciclosporin), an inhibitor of OATPs, BCRP and NTCP, on the exposure of rosuvastatin in healthy volunteers.

Conclusion

The results show the utility of the model to integrate a wide range of in vitro and in vivo data and simulate the outcome of clinical studies, with implications for their design.

Electronic supplementary material

The online version of this article (doi:10.1007/s40262-013-0097-y) contains supplementary material, which is available to authorized users.

Interaction Between Domperidone and Ketoconazole: Toward Prediction of Consequent QTc Prolongation Using Purely In Vitro Information

We aimed to investigate the application of combined mechanistic pharmacokinetic (PK) and pharmacodynamic (PD) modeling and simulation in predicting the domperidone (DOM) triggered pseudo-electrocardiogram modification in the presence of a CYP3A inhibitor, ketoconazole (KETO), using in vitro-in vivo extrapolation. In vitro metabolic and inhibitory data were incorporated into physiologically based pharmacokinetic (PBPK) models within Simcyp to simulate time course of plasma DOM and KETO concentrations when administered alone or in combination with KETO (DOM+KETO). Simulated DOM concentrations in plasma were used to predict changes in gender-specific QTcF (Fridericia correction) intervals within the Cardiac Safety Simulator platform taking into consideration DOM, KETO, and DOM+KETO triggered inhibition of multiple ionic currents in population. Combination of in vitro-in vivo extrapolation, PBPK, and systems pharmacology of electric currents in the heart was able to predict the direction and magnitude of PK and PD changes under coadministration of the two drugs although some disparities were detected.

Evaluation of an In Silico PBPK Post-Bariatric Surgery Model through Simulating Oral Drug Bioavailability of Atorvastatin and Cyclosporine

An increasing prevalence of morbid obesity has led to dramatic increases in the number of bariatric surgeries performed. Altered gastrointestinal physiology following surgery can be associated with modified oral drug bioavailability (Foral). In the absence of clinical data, an indication of changes to Foral via systems pharmacology models would be of value in adjusting dose levels after surgery. A previously developed virtual "post-bariatric surgery" population was evaluated through mimicking clinical investigations on cyclosporine and atorvastatin after bariatric surgery. Cyclosporine simulations displayed a reduced fraction absorbed through gut wall (fa) and Foral after surgery, consistent with reported observations. Simulated atorvastatin Foral postsurgery was broadly reflective of observed data with indications of counteracting interplay between reduced fa and an increased fraction escaping gut wall metabolism (FG). Inability to fully recover observed atorvastatin exposure after biliopancreatic diversion with duodenal switch highlights the current gap regarding the knowledge of associated biological changes.CPT: Pharmacometrics & Systems Pharmacology (2013) 2, e47; doi:10.1038/psp.2013.23; advance online publication 12 June 2013.

Are there differences in the catalytic activity per unit enzyme of recombinantly expressed and human liver microsomal cytochrome P450 2C9? A systematic investigation into inter-system extrapolation factors

The 'relative activity factor' (RAF) compares the activity per unit of microsomal protein in recombinantly expressed cytochrome P450 enzymes (rhCYP) and human liver without separating the potential sources of variation (i.e. abundance of enzyme per mg of protein or variation of activity per unit enzyme). The dimensionless 'inter-system extrapolation factor' (ISEF) dissects differences in activity from those in CYP abundance. Detailed protocols for the determination of this scalar, which is used in population in vitro-in vivo extrapolation (IVIVE), are currently lacking. The present study determined an ISEF for CYP2C9 and, for the first time, systematically evaluated the effects of probe substrate, cytochrome b5 and methods for assessing the intrinsic clearance (CL(int) ). Values of ISEF for S-warfarin, tolbutamide and diclofenac were 0.75 ± 0.18, 0.57 ± 0.07 and 0.37 ± 0.07, respectively, using CL(int) values derived from the kinetic values V(max) and K(m) of metabolite formation in rhCYP2C9 + reductase + b5 BD Supersomes™. The ISEF values obtained using rhCYP2C9 + reductase BD Supersomes™ were more variable, with values of 7.16 ± 1.25, 0.89 ± 0.52 and 0.50 ± 0.05 for S-warfarin, tolbutamide and diclofenac, respectively. Although the ISEF values obtained from rhCYP2C9 + reductase + b5 for the three probe substrates were statistically different (p < 0.001), the use of the mean value of 0.54 resulted in predicted oral clearance values for all three substrates within 1.4 fold of the observed literature values. For consistency in the relative activity across substrates, use of a b5 expressing recombinant system, with the intrinsic clearance calculated from full kinetic data is recommended for generation of the CYP2C9 ISEF. Furthermore, as ISEFs have been found to be sensitive to differences in accessory proteins, rhCYP system specific ISEFs are recommended.

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.

Prediction ofin vivodrug clearance fromin vitrodata. 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.

Mechanism-based inactivation of human cytochrome P450 enzymes: strategies for diagnosis and drug-drug interaction risk assessment

Among drugs that cause pharmacokinetic drug-drug interactions, mechanism-based inactivators of cytochrome P450 represent several of those agents that cause interactions of the greatest magnitude. In vitro inactivation kinetic data can be used to predict the potential for new drugs to cause drug interactions in the clinic. However, several factors exist, each with its own uncertainty, that must be taken into account in order to predict the magnitude of interactions reliably. These include aspects of in vitro experimental design, an understanding of relevant in vivo concentrations of the inactivator, and the extent to which the inactivated enzyme is involved in the clearance of the affected drug. Additionally, the rate of enzyme degradation in vivo is also an important factor that needs to be considered in the prediction of the drug interaction magnitudes. To address mechanism-based inactivation for new drugs, various in vitro experimental approaches have been employed. The selection of approaches for in vitro kinetic characterization of inactivation as well as in vitro-in vivo extrapolation should be guided by the purpose of the exercise and the stage of drug discovery and development, with an increase in the level of sophistication throughout the research and development process.

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.

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 <or=2. IVIVE predictions were accurate in 14 of 15 cases (mean-fold error range: 1.02-4.00). All five AS methods were accurate in 13, 11, 10, 10 and 14 cases, respectively. However, in some cases the error of AS exceeded fivefold. On the basis of the current results, IVIVE is more reliable than AS in predicting human clearance values for drugs mainly metabolized by CYP450 enzymes. This suggests that the place of AS methods in pre-clinical drug development warrants further scrutiny.

Influence of beta2-adrenoceptor gene polymorphisms on beta2-adrenoceptor-mediated responses in human lung mast cells

BACKGROUND AND PURPOSE

Previous studies have shown that beta(2)-adrenoceptor-mediated responses in human lung mast cells are highly variable. The aims of the present study were to establish whether polymorphisms of the beta (2)-adrenoceptor gene (ADRB2) influence this variability in (a) beta(2)-adrenoceptor-mediated inhibition and (b) desensitization of beta(2)-adrenoceptor-mediated responses in human lung mast cells.

EXPERIMENTAL APPROACH

Mast cells were isolated from human lung tissue. The inhibitory effects of the beta-adrenoceptor agonist, isoprenaline (10(-10)-10(-5) M), on IgE-mediated histamine release from mast cells were determined (n=92). Moreover, the inhibitory effects of isoprenaline were evaluated following a desensitizing treatment involving long-term (24 h) incubation of mast cells with isoprenaline (10(-6) M) (n=65). A potential influence of polymorphisms on these functional responses was determined by genotyping 11 positions, in the promoter and coding regions, of ADRB2 previously reported as polymorphic.

KEY RESULTS

There was no influence of any of the polymorphic positions of ADRB2 on the potency of isoprenaline to inhibit histamine release from mast cells with the exception of position 491C>T (Thr164Ile). There was no influence of any of the polymorphic positions of ADRB2 on the extent of desensitization of the isoprenaline-mediated response following a desensitizing treatment except for position 46G>A (Gly16Arg). Analyses at the haplotype level indicated that there was no influence of haplotype on beta (2)-adrenoceptor-mediated responses in mast cells.

CONCLUSIONS AND IMPLICATIONS

These data indicate that certain polymorphisms in ADRB2 influence beta(2)-adrenoceptor-mediated responses in human lung mast cells.

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.

Determination by liquid chromatography-mass spectrometry of clomiphene isomers in the plasma of patients undergoing treatment for the induction of ovulation

A rapid, sensitive and selective LC-MS method is described for the simultaneous determination of zuclomiphene and enclomiphene in plasma from patients undergoing treatment with clomiphene citrate for the induction of ovulation. Samples spiked with N-didesmethyltamoxifen, the internal standard, were extracted into methyl tertiary butyl ether. The compounds were separated on a Luna C(18) analytical column, and a mobile phase of methanol-water (70:30 v/v) containing 0.05% trifluoroacetic acid at a flow rate of 1ml/min. The limits of determination were 35pg/ml and 7pg/ml for zu- and enclomiphene, respectively. Within-day coefficients of variation ranged from 2.1% to 7.2%.

A six-year evaluation of methadone prescribing practices at a substance misuse treatment centre in the UK

OBJECTIVES

To document changes in prescribing practice at a specialized substance misuse service in the UK occurring since the introduction of the 1999 UK National Guidelines on the management of drug misuse, and to explore a possible link between the length of time spent in methadone maintenance therapy (MMT) and the dosage prescribed.

METHODS

A retrospective analysis of a computerized prescription database between 1996 and 2002 obtained from Sheffield Care Trust Substance Misuse Service was performed. The relationship between various measures of dosage and the length of time spent in MMT was investigated.

RESULTS

In accordance with the 1999 UK National Guidelines, the proportion of injectable methadone prescribed decreased from 22% to 16%. This was offset by an increase in the prescribing of methadone elixir from 74% to 79%. The 'maximum dose' of methadone prescribed correlated significantly with patient retention, explaining 14% of the variation in time spent in MMT.

CONCLUSIONS

Our findings indicate that publication of the UK National Guidelines had a measurable effect on prescribing practice at the Service. We found that a higher methadone dose is associated with increased patient retention in MMT. However, as only a maximum of 14% of the variation in the length of stay is related to methadone dose, the importance of other aspects of treatment such as counselling and rehabilitation programmes, should be considered for the successful treatment of opioid abusers.

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.

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.

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.

Timing and type of glucocorticoid replacement in adult congenital adrenal hyperplasia

Congenital adrenal hyperplasia (CAH) is one of the commonest inherited diseases. Treatment during childhood is directed towards obtaining normal growth. There is an extensive literature on management of CAH during childhood but little published on how patients should be treated as adults. CAH results in increased adrenocorticotropic hormone levels driving the adrenal to produce cortisol and thereby excessive cortisol precursors which are predominantly androgenic. In normal individuals there is a marked circadian rhythm in cortisol release with the lowest levels occurring shortly after midnight and rising between 02.00 and 03.00 h to peak at around 06.00-08.00 h after waking. The treatment of choice for children with CAH is hydrocortisone three times daily. There is no consensus on treating adults. A recent survey by the UK Society for Endocrinology showed that hydrocortisone was the most popular steroid, followed by dexamethasone and then prednisolone. Sixty percent of respondents used a reverse circadian pattern of glucocorticoid treatment and only 16% used a dose based on weight or surface area. In a recent study we demonstrated that the most important variable determining hydrocortisone bioavailability is weight and we advocated a weight-related dosing regimen monitored using a nomogram for serum cortisol. .

Preliminary evidence of the cardiovascular effects of polysubstance misuse in nightclubs

Adverse reactions to polysubstance misuse are common in nightclubs, yet little is known of the physiological effects of polysubstance misuse in this environment. This study examined the heart rate, blood pressure and oral temperature of 50 participants recruited in a nightclub on four separate nights. In addition, the increase in environmental temperature was recorded throughout each night. There were no differences in oral temperature between polysubstance misusers (i.e. those who used Ecstasy, amphetamine, cocaine, alcohol and cannabis) and alcohol/cannabis misusers. On the other hand, there were significant differences in both heart rate and blood pressure between the two groups. These data suggest that polysubstance misusers may be at risk from cardio- and cerebrovascular toxicity. Further field/on-site work is clearly needed to investigate the effects of polysubstance misuse.

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.

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.

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%.

Contribution of midazolam and its 1-hydroxy metabolite to preoperative sedation in children: a pharmacokinetic-pharmacodynamic analysis

BACKGROUND

Oral midazolam is widely used for preoperative sedation in children. We have studied the pharmacokinetics (PK) of both midazolam and its active 1-hydroxy metabolite and their contributions to sedative effect in 45 children attending for day surgery.

METHODS

Blood samples (two per individual) were collected at the beginning and end of the surgical procedure. Plasma midazolam and 1-hydroxymidazolam (1-OHMDZ) were measured by HPLC. Sedation score (score: 1 = awake, 2 = drowsy/asleep) was recorded at the same time as the first blood sample. The population-PK software P-Pharm was used to analyse the data. Age, weight, sex, concomitant drugs, and the metabolic ratio, 1-OHMDZ/midazolam were investigated as co-variates of the PK of midazolam and 1-OHMDZ. The pharmacokinetic-pharmacodynamic (PK-PD) modelling of the score in relation to plasma midazolam and 1-OHMDZ was performed using logistic regression analysis.

RESULTS

A median dose of 0.5 mg kg-1 was given to the children, median age 5 yr (range from 9 months to 12 yr) and weight 21 kg (range 8-75 kg). Average concentrations of midazolam 150 ng ml-1 and 1-OHMDZ 90 ng ml-1 were observed in the first plasma samples. These concentrations resulted in an odds ratio of 4 in favour of score 2 vs 1. The best PK-PD model included both midazolam and 1-OHMDZ as active moieties and predicted correct scores in 86% of cases.

CONCLUSION

Studies of midazolam should evaluate the contribution of 1-OHMDZ to the overall PD effect. The metabolite 1-OHMDZ has approximately half the activity of the parent drug and can compensate for at least part of the decreased effect due to increased midazolam metabolism.

A new rapidly absorbed paracetamol tablet containing sodium bicarbonate. II. Dissolution studies and in vitro/in vivo correlation

The objective of this study was to compare the in vitro dissolution profile of a new rapidly absorbed paracetamol tablet containing sodium bicarbonate (PS) with that of a conventional paracetamol tablet (P), and to relate these by deconvolution and mapping to in vivo release. The dissolution methods used include the standard procedure described in the USP monograph for paracetamol tablets, employing buffer at pH 5.8 or 0.05 M HCl at stirrer speeds between 10 and 50 rpm. The mapping process was developed and implemented in Microsoft Excel worksheets that iteratively calculated the optimal values of scale and shape factors which linked in vivo time to in vitro time. The in vitro-in vivo correlation (IVIVC) was carried out simultaneously for both formulations to produce common mapping factors. The USP method, using buffer at pH 5.8, demonstrated no difference between the two products. However, using an acidic medium the rate of dissolution of P but not of PS decreased with decreasing stirrer speed. A significant correlation (r = 0.773; p < .00001) was established between in vivo release and in vitro dissolution using the profiles obtained with 0.05 M HCl and a stirrer speed of 30 rpm. The scale factor for optimal simultaneous IVIVC in the fasting state was 2.54 and the shape factor was 0.16; corresponding values for mapping in the fed state were 3.37 and 0.13 (implying a larger in vitro-in vivo time difference but reduced shape difference in the fed state). The current IVIVC explains, in part, the observed in vivo variability of the two products. The approach to mapping may also be extended to different batches of these products, to predict the impact of any changes of in vitro dissolution on in vivo release and plasma drug concentration-time profiles.

A new rapidly absorbed paracetamol tablet containing sodium bicarbonate. I. A four-way crossover study to compare the concentration-time profile of paracetamol from the new paracetamol/sodium bicarbonate tablet and a conventional paracetamol tablet in fed and fasted volunteers

The primary objective of this four-way crossover study was to compare the concentration-time profile of paracetamol from a new rapidly absorbed paracetamol tablet containing sodium bicarbonate (PS) with a conventional paracetamol tablet (P), in a panel of 28 fed and fasted healthy volunteers. The results demonstrated that paracetamol was absorbed more rapidly from tablets containing sodium bicarbonate compared to conventional tablets, as indicated by a shorter tmax in both the fed and fasted state and a higher Cmax in the fasted state. The two formulations were bioequivalent with respect to area under curve (AUC). Food did not affect the extent of absorption from either formulation, as indicated by AUC, however, food did reduce the rate of absorption from both formulations, as indicated by a longer tmax and a lower Cmax. Metabolic activation of paracetamol to its oxidation metabolites, as assessed by combined partial clearances to subsequent secondary metabolites cysteine and mercapturic acid conjugates, indicated that the two formulations were bioequivalent in this respect.

Analysis of solvent central nervous system toxicity and ethanol interactions using a human population physiologically based kinetic and dynamic model

The effect of acute ethanol-mediated inhibition of m-xylene metabolism on central nervous system (CNS) depression in the human worker population was investigated using physiologically based pharmacokinetic (PBPK) models and probabilistic random (Monte Carlo) sampling. PBPK models of inhaled m-xylene and orally ingested ethanol were developed and combined by a competitive enzyme (CYP2E1) inhibition model. Human interindividual variability was modeled by combining estimated statistical distributions of model parameters with the deterministic PBPK models and multiple random or Monte Carlo simulations. A simple threshold pharmacodynamic model was obtained by simulating m-xylene kinetics in human studies where CNS effects were observed and assigning the peak venous blood m-xylene concentration (C(V,max)) as the dose surrogate of toxicity. Probabilistic estimates of an individual experiencing CNS disturbances given exposure to the current UK occupational exposure standard (100 ppm time-weighted average over 8 h), with and without ethanol ingestion, were obtained. The probability of experiencing CNS effects given this scenario increases markedly and nonlinearly with ethanol dose. As CYP2E1-mediated metabolism of other occupationally relevant organic compounds may be inhibited by ethanol, simulation studies of this type should have an increasingly significant role in the chemical toxicity risk assessment.

The placebo response to citric acid-induced cough: pharmacodynamics and gender differences

Characteristics of the response to placebo in a citric acid-induced cough challenge were investigated as part of a randomized, double-blind crossover trial to assess the antitussive effect of dextromethorphan. Baseline cough responses were established on two occasions in 22 healthy subjects. They received 60 ml placebo antitussive syrup and cough frequency following five inhalations of 10% citric acid over 5 min was measured at regular intervals up to 12 h. Response-time models of varying complexity were used to describe the placebo cough suppression data. The cough response to placebo was also compared to that of the untreated state. The placebo cough response was best characterized by a non-linear increase in cough suppression up to a maximum reduction of 1.6 coughs from baseline at 4-4.5 h, followed by a non-linear return to baseline. The cough response in the untreated state was not different from that of placebo (P=0.99). Females coughed more frequently than males (median number of coughs=10.5 vs. 9.0, respectively P<0.001; Mann-Whitney U test), and adaptation to the cough stimulus was significantly more rapid in females (P<0.025). Accordingly, in trials that use citric acid-induced cough, gender should be considered in study design, particularly in relation to the timing of measurements.