Pharmacokinetics of a 503B outsourcing facility-produced theophylline in dogs.

Theophylline is an important drug for treatment of canine chronic bronchitis and bradyarrhythmias, but new products require validation since pharmacokinetics in dogs can vary by formulation. A new, 503B outsourcing facility-produced theophylline product (OFT) is available for veterinary use. Outsourcing facilities have many advantages over traditional compounding sources including current good manufacturing practice compliance. The purpose of this study was to establish the pharmacokinetics of OFT in dogs. Eight healthy dogs received 11 mg/kg intravenous aminophylline and 10 mg/kg oral OFT followed by serial blood sampling in a two-way, randomized, crossover design with 7-day washout. Plasma theophylline concentrations were quantified by liquid chromatography-mass spectrometry. Bioavailability, maximum concentration, time to maximum concentration, half-life and area under the curve were: 97 ± 10%, 7.13 ± 0.71 µg/mL, 10.50 ± 2.07 h, 9.20 ± 2.87 h, and 141 ± 37.6 µg*h/mL, respectively. Steady-state predictions supported twice daily dosing of the OFT, but specific dosage recommendations are hindered by lack of a canine-specific therapeutic range for plasma theophylline concentration. These findings suggest that the OFT is well absorbed and can likely be dosed twice daily in dogs, but future pharmacodynamic and clinical studies are needed to establish a definitive therapeutic range for theophylline in this species.

Theophylline dosing and pharmacokinetics for renal protection in neonates with hypoxic-ischemic encephalopathy undergoing therapeutic hypothermia.

BACKGROUND: Theophylline, a non-selective adenosine receptor antagonist, improves renal perfusion in the setting of hypoxia-ischemia and may offer therapeutic benefit in neonates with hypoxic-ischemic encephalopathy (HIE) undergoing hypothermia. We evaluated the pharmacokinetics and dose-exposure relationships of theophylline in this population to guide dosing strategies. METHODS: A population pharmacokinetic analysis was performed in 22 neonates with HIE undergoing hypothermia who were part of a prospective study or retrospective chart review. Aminophylline (intravenous salt form of theophylline) was given per institutional standard of care for low urine output and/or rising serum creatinine (5 mg/kg intravenous (i.v.) load then 1.8 mg/kg i.v. q6h). The ability of different dosing regimens to achieve target concentrations (4-10 mg/L) associated with clinical response was examined. RESULTS: Birth weight was a significant predictor of theophylline clearance and volume of distribution (p < 0.05). The median half-life was 39.5 h (range 27.2-50.4). An aminophylline loading dose of 7 mg/kg followed by 1.6 mg/kg q12h was predicted to achieve target concentrations in 84% of simulated neonates. CONCLUSIONS: In neonates with HIE undergoing hypothermia, theophylline clearance was low with a 50% longer half-life compared to full-term normothermic neonates without HIE. Dosing strategies need to consider the unique pharmacokinetic needs of this population. IMPACT: Theophylline is a potential renal-protective therapy in neonates with HIE undergoing therapeutic hypothermia; however, the pharmacokinetics and dose needs in this population are not known. Theophylline clearance was low in neonates with HIE undergoing therapeutic hypothermia with a 50% longer half-life compared to full-term normothermic neonates without HIE. As theophylline is advanced in clinical development, dosing strategies will need to consider the unique pharmacokinetic needs of neonates with HIE undergoing therapeutic hypothermia.

Well-tempered MCMC simulations for population pharmacokinetic models.

A full Bayesian statistical treatment of complex pharmacokinetic or pharmacodynamic models, in particular in a population context, gives access to powerful inference, including on model structure. Markov Chain Monte Carlo (MCMC) samplers are typically used to estimate the joint posterior parameter distribution of interest. Among MCMC samplers, the simulated tempering algorithm (TMCMC) has a number of advantages: it can sample from sharp multi-modal posteriors; it provides insight into identifiability issues useful for model simplification; it can be used to compute accurate Bayes factors for model choice; the simulated Markov chains mix quickly and have assured convergence in certain conditions. The main challenge when implementing this approach is to find an adequate scale of auxiliary inverse temperatures (perks) and associated scaling constants. We solved that problem by adaptive stochastic optimization and describe our implementation of TMCMC sampling in the GNU MCSim software. Once a grid of perks is obtained, it is easy to perform posterior-tempered MCMC sampling or likelihood-tempered MCMC (thermodynamic integration, which bridges the joint prior and the posterior parameter distributions, with assured convergence of a single sampling chain). We compare TMCMC to other samplers and demonstrate its efficient sampling of multi-modal posteriors and calculation of Bayes factors in two stylized case-studies and two realistic population pharmacokinetic inference problems, one of them involving a large PBPK model.

Assessment of CYP-Mediated Drug Interactions for Evocalcet, a New Calcimimetic Agent, Based on In Vitro Investigations and a Cocktail Study in Humans.

Evocalcet is a novel calcimimetic agent for the treatment of secondary hyperparathyroidism (SHPT). This study evaluated the effects of evocalcet on inhibition and induction of cytochrome P450 (CYP) isozymes. Although drug interactions arising from reversible inhibition of CYP isozymes by evocalcet were considered unlikely based on the results of in vitro studies and static model analyses, the potential for evocalcet to cause time-dependent inhibition of CYP3A or induction of several CYP isozymes could not be ruled out. Therefore, a clinical drug-drug interaction (DDI) study to evaluate the effects of evocalcet on the pharmacokinetics (PKs) of probe substrates for CYP isozymes (CYP1A2, CYP2B6, CYP2C8, CYP2C9, and CYP3A) was conducted in healthy male volunteers using a novel cocktail combination. Evocalcet did not significantly affect the PKs of the probe substrates, confirming that CYP-mediated interactions were unlikely.

Pharmacometrics models with hidden Markovian dynamics.

The aim of this paper is to provide an overview of pharmacometric models that involve some latent process with Markovian dynamics. Such models include hidden Markov models which may be useful for describing the dynamics of a disease state that jumps from one state to another at discrete times. On the contrary, diffusion models are continuous-time and continuous-state Markov models that are relevant for modelling non observed phenomena that fluctuate continuously and randomly over time. We show that an extension of these models to mixed effects models is straightforward in a population context. We then show how the forward-backward algorithm used for inference in hidden Markov models and the extended Kalman filter used for inference in diffusion models can be combined with standard inference algorithms in mixed effects models for estimating the parameters of the model. The use of these models is illustrated with two applications: a hidden Markov model for describing the epileptic activity of a large number of patients and a stochastic differential equation based model for describing the pharmacokinetics of theophyllin.

Control of methylxanthines in the competition horse: pharmacokinetic/pharmacodynamic studies on caffeine, theobromine and theophylline for the assessment of irrelevant concentrations.

Methylxanthines positives in competition samples have challenged doping control laboratories and racing jurisdictions since methylxanthines are naturally occurring prohibited substances and often constituents of feed. For theobromine, an international threshold (renamed in International Residue Limit, IRL) of 2 µg/mL in urine has been established. On the basis of the data presented herein, a threshold or rather an IRL for theobromine in plasma of 0.3 µg/mL was proposed and was thereupon approved by the International Federation of Horseracing Authorities (IFHA). Official recommendations for reporting caffeine and theophylline are still lacking. The aim of the study was to investigate IRLs for theobromine in blood and for caffeine and theophylline in blood and urine. Therefore, a set of six administrations were carried out including both single i.v. and single oral administrations of caffeine, theobromine and theophylline. Plasma and urine concentrations were determined using a validated liquid chromatography-tandem mass spectrometry (LC-MS/MS). Applying the Toutain model approach an effective plasma concentration (EPC) of caffeine was estimated at 3.05 µg/mL, irrelevant concentrations in blood (IPC) and urine (IUC) approached 6 and 12 ng/mL, respectively. EPC of theobromine was calculated with 3.80 µg/mL, and irrelevant concentrations of theobromine were determined at 8 ng/mL in plasma and at 142 ng/mL in urine. Toutain modelling of the theophylline data produced an EPC, IPC, and IUC of 3.20 µg/mL, 6 ng/mL, and 75 ng/mL, respectively. The obtained irrelevant concentrations were used to postulate IRLs for theobromine in plasma and for caffeine and theophylline in plasma and urine. Copyright © 2016 John Wiley & Sons, Ltd.

Prediction of inter-individual variability on the pharmacokinetics of CYP1A2 substrates in non-smoking healthy volunteers.

The activity of CYP1A2, a major drug-metabolizing enzyme, is known to be affected by various environmental factors. Our study aimed to predict inter-individual variability of AUC/Dose of CYP1A2 substrates in non-smoking healthy volunteers using the Monte Carlo simulation. Inter-individual variability in hepatic intrinsic clearance of CYP1A2 substrates (CLint,h,1A2) was estimated using dispersion model based on the inter-individual variability (N = 96) of the AUC of caffeine, a major CYP1A2 substrate. The estimated coefficient of variation (CV) of CLint,h,1A2 was 55%, similar to previously reported CLint,h,2D6 (60%) but larger than CLint,h,3A4 (33%). Then, this estimated CV was validated by predicting the CVs of AUC/Dose of tizanidine and phenacetin, which are mainly metabolized by CYP1A2 and have negligible renal clearance. As a result, reported CVs were successfully predicted within 2.5-97.5 percentile range of predicted values. Moreover, CVs for AUC/Dose of the CYP1A2 substrates theophylline and lidocaine, which are affected by other CYPs and renal clearance, were also successfully predicted. The inter-individual variability of AUC/Dose of CYP1A2 substrates was successfully predicted using 55% CV for CLint,h,1A2, and the results, along with those reported by our group for other CYPs, support the prediction of inter-individual variability of pharmacokinetics in the clinical setting.

Bayesian Population Physiologically-Based Pharmacokinetic (PBPK) Approach for a Physiologically Realistic Characterization of Interindividual Variability in Clinically Relevant Populations.

Interindividual variability in anatomical and physiological properties results in significant differences in drug pharmacokinetics. The consideration of such pharmacokinetic variability supports optimal drug efficacy and safety for each single individual, e.g. by identification of individual-specific dosings. One clear objective in clinical drug development is therefore a thorough characterization of the physiological sources of interindividual variability. In this work, we present a Bayesian population physiologically-based pharmacokinetic (PBPK) approach for the mechanistically and physiologically realistic identification of interindividual variability. The consideration of a generic and highly detailed mechanistic PBPK model structure enables the integration of large amounts of prior physiological knowledge, which is then updated with new experimental data in a Bayesian framework. A covariate model integrates known relationships of physiological parameters to age, gender and body height. We further provide a framework for estimation of the a posteriori parameter dependency structure at the population level. The approach is demonstrated considering a cohort of healthy individuals and theophylline as an application example. The variability and co-variability of physiological parameters are specified within the population; respectively. Significant correlations are identified between population parameters and are applied for individual- and population-specific visual predictive checks of the pharmacokinetic behavior, which leads to improved results compared to present population approaches. In the future, the integration of a generic PBPK model into an hierarchical approach allows for extrapolations to other populations or drugs, while the Bayesian paradigm allows for an iterative application of the approach and thereby a continuous updating of physiological knowledge with new data. This will facilitate decision making e.g. from preclinical to clinical development or extrapolation of PK behavior from healthy to clinically significant populations.

Assessment of the hepatic CYP reductase null mouse model and its potential application in drug discovery.

CYP Oxidoreductase (Por) is the essential electron donor for all CYP enzymes and is responsible for the activation of CYP. The Taconic Hepatic CYP Reductase Null (HRN) mouse model possesses a targeted mutation that results in liver-specific deletion of the Por gene thereby resulting in a disruption of CYP metabolism in the liver. The objectives of these studies were to further characterize the HRN mouse using probe drugs metabolized by CYP. In addition, tumor exposure in xenograft tumor bearing HRN immune-compromised (nude) mice was also determined. In HRN mice following intravenous (iv) administration of midazolam, clearance (CL) was reduced by ∼ 80% compared to wild-type mice (WT). After oral administration, the AUC of midazolam was increased by ∼ 20-fold in HRN mice compared to WT mice; this greater effect suggests that hepatic first pass plays a role in the oral CL of midazolam. A 50% and an 80% decrease in CL were also observed in HRN mice following iv administration of docetaxel and theophylline, respectively, compared to WT mice. In addition, a 2- to 3-fold increase in tumor concentrations of G4222, a tool compound, were observed in tumor bearing HRN nude mice compared to tumor bearing nude WT mice. The observations from these experiments demonstrate that, for compounds that are extensively metabolized by hepatic CYP, the HRN mouse model could potentially be valuable for evaluating in vivo efficacy of tool compounds in drug discovery where high hepatic CL and low exposure may prevent in vivo evaluation of a new chemical entity.

Application of a systems approach to the bottom-up assessment of pharmacokinetics in obese patients: expected variations in clearance.

BACKGROUND AND OBJECTIVES: The maintenance dose of a drug is dependent on drug clearance, and thus any biochemical and physiological changes in obesity that affect parameters such as cardiac output, renal function, expression of drug-metabolizing enzymes and protein binding may result in altered clearance compared with that observed in normal-weight subjects (corrected or uncorrected for body weight). Because of the increasing worldwide incidence of obesity, there is a need for more information regarding the optimal dosing of drug therapy to be made available to prescribers. This is usually provided via clinical studies in obese people; however, such studies are not available for all drugs that might be used in obese subjects. Incorporation of the relevant physiological and biochemical changes into predictive bottom-up pharmacokinetic models in order to optimize dosage regimens may offer a logical way forward for the cases where no clinical data exist. The aims of the current report are to apply such a 'systems approach' to identify the likelihood of observing variations in the clearance of drugs in obesity and morbid obesity for a set of compounds for which clinical data, as well as the necessary in vitro information, are available, and to provide a framework for assessing other drugs in the future. METHODS: The population-specific changes in demographic, physiological and biochemical parameters that are known to be relevant to obese and morbidly obese subjects were collated and incorporated into two separate population libraries. These libraries, together with mechanistic in vitro-in vivo extrapolations (IVIVE) within the Simcyp Population-based Simulator™, were used to predict the clearance of oral alprazolam, oral caffeine, oral chlorzoxazone, oral ciclosporin, intravenous and oral midazolam, intravenous phenytoin, oral theophylline and oral triazolam. The design of the simulated studies was matched as closely as possible with that of the clinical studies. Outcome was measured by the predicted ratio of the clearance of the drug in obese and lean subjects ± its 90% confidence interval, compared with observed values. The overall statistical measures of the performance of the model to detect differences in compound clearance between obese and lean populations were investigated by measuring sensitivity, specificity, positive predictive value (PPV) and negative predictive value (NPV). A power calculation was carried out to investigate the impact of the sample size on the overall outcome of clinical studies. RESULTS: The model was successful in predicting clearance in obese subjects, with the degree to which simulations could mimic the outcome of in vivo studies being greater than 60% for six of the eight drugs. A clear difference in the clearance of chlorzoxazone was correctly picked up via simulation. The overall statistical measures of the performance of the Simcyp Simulator were 100% sensitivity, 66% specificity, 60% PPV and 100% NPV. Studies designed on the basis of the ratio of the absolute values required substantial numbers of participants in order to detect a significant difference, except for phenytoin and chlorzoxazone, where the ratios of the weight-normalized clearances generally showed statistically significant differences with a smaller number of subjects. CONCLUSION: Extension of a mechanistic predictive pharmacokinetic model to accommodate physiological and biochemical changes associated with obesity and morbid obesity allowed prediction of changes in drug clearance on the basis of in vitro data, with reasonable accuracy across a range of compounds that are metabolized by different enzymes. Prediction of the effects of obesity on drug clearance, normalized by various body size scalars, is of potential value in the design of clinical studies during drug development and in the introduction of dosage adjustments that are likely to be needed in clinical practice.

Caffeine and paraxanthine HPLC assay for CYP1A2 phenotype assessment using saliva and plasma.

Caffeine has been extensively used as a probe to measure CYP1A2 activity in humans with caffeine clearance or the paraxanthine (major metabolite of caffeine) to caffeine concentration ratio being regarded as the preferred metric. A simple reverse-phased C(18) HPLC assay using ethyl acetate liquid-liquid extraction was developed to quantitate caffeine and paraxanthine concentrations in saliva and plasma. The mobile phase consisted of acetonitrile-acetic acid-H(2)O (100:1:899) and analytes were quantitated with UV detection at 280 nm. The extraction recovery for paraxanthine and caffeine was approximately 70% in both saliva and plasma. The assay was linear over the concentration ranges 0.05-2.50 and 0.05-5.00 µg/mL, for paraxanthine and caffeine, respectively, in saliva. In plasma the assay was linear over the ranges 0.025-2.50 and 0.025-5.00 µg/mL for paraxanthine and caffeine, respectively. Intra- and inter-assay precision and accuracy were less than 15%. Detection limits were 0.015 µg/mL for paraxanthine and caffeine in saliva, while it was 0.005 µg/mL for paraxanthine and caffeine in plasma. Utility was established in samples collected from two healthy volunteers who abstained from caffeine for 24 h and received a single 100 mg oral dose of caffeine. The assay developed is a robust, simple and precise technique to measure caffeine and paraxanthine in saliva and plasma of healthy volunteers after a single oral dose of caffeine.

A new definition of pharmaceutical quality: assembly of a risk simulation platform to investigate the impact of manufacturing/product variability on clinical performance.

The absence of a unanimous, industry-specific definition of quality is, to a certain degree, impeding the progress of ongoing efforts to "modernize" the pharmaceutical industry. This work was predicated on requests by Dr. Woodcock (FDA) to re-define pharmaceutical quality in terms of risk by linking production characteristics to clinical attributes. A risk simulation platform that integrates population statistics, drug delivery system characteristics, dosing guidelines, patient compliance estimates, production metrics, and pharmacokinetic, pharmacodynamic, and in vitro-in vivo correlation models to investigate the impact of manufacturing variability on clinical performance of a model extended-release theophylline solid oral dosage system was developed. Manufacturing was characterized by inter- and intra-batch content uniformity and dissolution variability metrics, while clinical performance was described by a probabilistic pharmacodynamic model that expressed the probability of inefficacy and toxicity as a function of plasma concentrations. Least-squares regression revealed that both patient compliance variables, percent of doses taken and dosing time variability, significantly impacted efficacy and toxicity. Additionally, intra-batch content uniformity variability elicited a significant change in risk scores for the two adverse events and, therefore, was identified as a critical quality attribute. The proposed methodology demonstrates that pharmaceutical quality can be recast to explicitly reflect clinical performance.

Substitution of theophylline slow-release formulations according to the rebate contracts in the German statutory health insurance.

On the basis of the rebate contracts between individual statutory health insurance funds and pharmaceutical enterprises, the generic substitution of prescribed medications is economically attractive and is advocated for statutory health insurees in Germany. In addition to the drugs whose substitution can be considered to be uncritical, rebate contracts also include controversial substances such as the bronchodilator theophylline (CAS 58-55-9), which has a narrow therapeutic range and should only be substituted under certain conditions. The objective of this article was to check the safety of the substitution of theophylline by means of a comparative evaluation of bioequivalence studies carried out on theophylline slow-release preparations. A systematic literature search was carried out in the MEDLINE database. The search terms used were combinations of the following key words: theophylline, generics, bioequivalence, substitution, brand and non-brand. In addition, a manual search was performed in the reference lists of the relevant articles. Only articles that were published between January 1, 1988 and August 30, 2008 were to be included. Five studies conformed to the inclusion and exclusion criteria. Two of the studies came to the conclusion that the preparations analysed were bioequivalent. In the remaining three studies there was no bioequivalence found between the preparations and the reference product. Because of the heterogeneity of study outcomes no metanalysis could be performed. On the basis of the studies analysed the conclusion can be drawn that a theophylline slow-release preparation should only be substituted under close monitoring by a physician because of the many factors which can adversely affect serum levels, such as the narrow therapeutic range of the active ingredient, the patient's metabolisation rate or the different galenics of the preparations. Nevertheless, the question remains as to whether the costs saved by the rebate contracts would not be significantly outweighed.

Extension of the SAEM algorithm for nonlinear mixed models with 2 levels of random effects.

This article focuses on parameter estimation of multilevel nonlinear mixed-effects models (MNLMEMs). These models are used to analyze data presenting multiple hierarchical levels of grouping (cluster data, clinical trials with several observation periods, ...). The variability of the individual parameters of the regression function is thus decomposed as a between-subject variability and higher levels of variability (e.g. within-subject variability). We propose maximum likelihood estimates of parameters of those MNLMEMs with 2 levels of random effects, using an extension of the stochastic approximation version of expectation-maximization (SAEM)-Monte Carlo Markov chain algorithm. The extended SAEM algorithm is split into an explicit direct expectation-maximization (EM) algorithm and a stochastic EM part. Compared to the original algorithm, additional sufficient statistics have to be approximated by relying on the conditional distribution of the second level of random effects. This estimation method is evaluated on pharmacokinetic crossover simulated trials, mimicking theophylline concentration data. Results obtained on those data sets with either the SAEM algorithm or the first-order conditional estimates (FOCE) algorithm (implemented in the nlme function of R software) are compared: biases and root mean square errors of almost all the SAEM estimates are smaller than the FOCE ones. Finally, we apply the extended SAEM algorithm to analyze the pharmacokinetic interaction of tenofovir on atazanavir, a novel protease inhibitor, from the Agence Nationale de Recherche sur le Sida 107-Puzzle 2 study. A significant decrease of the area under the curve of atazanavir is found in patients receiving both treatments.

Application of a sample pooling method for the accelerated assessment of the rate of uptake of drugs by the brain in rats.

The purpose of this study was to examine the feasibility of applying a sample pooling method to the accelerated estimation of the uptake clearance of drugs to the brain in rats. Brain uptake clearances (CL(uptake)) were estimated for five model compounds using the sample pooling method and an integration plot analysis. CL(uptake) was also evaluated for caffeine and theophylline by brain microdialysis. The parameters and throughput of the pooling method were compared with those of typically used standard methods. The correlation for CL(uptake) was statistically significant (P<0.005) between the integration plot and the current method; the throughput of evaluation was 15-fold higher for the sample pooling method. A comparison of CL(uptake) values indicated that the three methods showed comparable results for caffeine while the CL(uptake) of theophylline using the proposed method was significantly different from those of the other methods. A kinetic analysis indicated that a compound with a slower CL(uptake) and longer half-life (e.g., theophylline) is more prone to error and that the lower limit of the CL(uptake) of 0.17 mL min(-1) (g brain)(-1) may be set so as to have an error less than 20% of the estimation. These results suggest that the sample pooling method is applicable for use in the accelerated estimation of the uptake clearance of compounds in the brain for which the value is greater than 0.17 mL min(-1) (g brain)(-1).

Assessment of CYP1A2 activity in clinical practice: why, how, and when?

The cytochrome P450 enzyme CYP1A2 mediates the rate-limiting step in the metabolism of many drugs including theophylline, clozapine, and tacrine as well as in the bioactivation of procarcinogens. CYP1A2 activity shows both pronounced intra- and interindividual variability, which is, among other factors, related to smoking causing enzyme induction, to drug intake and to dietary factors which may result in induction or inhibition. In contrast to these exogenous factors, genetic influences on enzyme activity seem to be less pronounced. Therefore, phenotyping of CYP1A2, i.e. the determination of the actual activity of the enzyme in vivo, represents a useful approach both for scientific and clinical applications. CYP1A2 is almost exclusively expressed in the liver. Since liver tissue cannot be obtained for direct phenotyping, a probe drug which is metabolized by CYP1A2 has to be given. Proposed probe drugs include caffeine, theophylline, and melatonin. Caffeine is most often used because of the predominating role of CYP1A2 in its overall metabolism and the excellent tolerability. Various urinary, plasma, saliva, and breath based CYP1A2 caffeine metrics have been applied. While caffeine clearance is considered as the gold standard, the salivary or plasma ratio of paraxanthine to caffeine in a sample taken approximately 6 hr after a defined dose of caffeine is a more convenient, less expensive but also fully validated CYP1A2 phenotyping metric. CYP1A2 phenotyping is applied frequently in epidemiologic and drug-drug interaction studies, but its clinical use and usefulness remains to be established.

Non-linear mixed-effects pharmacokinetic/pharmacodynamic modelling in NLME using differential equations.

The standard software for non-linear mixed-effect analysis of pharmacokinetic/pharmacodynamic (PK/PD) data is NONMEM while the non-linear mixed-effects package NLME is an alternative as long as the models are fairly simple. We present the nlmeODE package which combines the ordinary differential equation (ODE) solver package odesolve and the non-linear mixed effects package NLME thereby enabling the analysis of complicated systems of ODEs by non-linear mixed-effects modelling. The pharmacokinetics of the anti-asthmatic drug theophylline is used to illustrate the applicability of the nlmeODE package for population PK/PD analysis using the available data analysis tools in R for model inspection and validation. The nlmeODE package is numerically stable and provides accurate parameter estimates which are consistent with NONMEM estimates.

Uncertainty factors for chemical risk assessment: interspecies differences in the in vivo pharmacokinetics and metabolism of human CYP1A2 substrates.

The 100-fold default uncertainty factor is used to convert a no-observed-adverse-effect level (NOAEL) from a animal toxicity study, to a "safe" value for human intake. The composite uncertainty factor (100) has to allow for interspecies (10-fold) and interindividual (10-fold) differences in toxicokinetics and toxicodynamics. The aim of the current study was to assess the validity of the interspecies default for toxicokinetics (4.0) for each of the test species (dog, rabbit, rat and mouse), using published data for compounds eliminated by CYP1A2 in humans (caffeine, theobromine, theophylline and paraxanthine). An analysis of the published literature showed that the absorption, bioavailability and route of excretion were generally similar between humans and the test species, for each probe substrate. However, interspecies differences in the route of metabolism, and the enzymes involved in this process, were identified. The magnitude of difference in the internal dose, between species, showed that values for the mouse (10.6) and rat (5.4) exceed the 4.0-fold default, whereas the rabbit (2.6) and dog (1.6) were below this value. This work supports the need to replace the generic default factors by a compound-related value derived from specific, relevant, quantitative data; this would result in more relevant and reliable non-cancer risk assessments.

Uncertainty factors for chemical risk assessment. human variability in the pharmacokinetics of CYP1A2 probe substrates.

A 100-fold uncertainty factor is used to derive acceptable daily intakes for compounds causing thresholded toxicity. The 10-fold factor for human variability can be further subdivided into two factors of 10(0.5) (3.16) to allow for toxicokinetics and toxicodynamics. The validity of the human kinetic subfactor has been analysed in relation to CYP1A2 metabolism using published in vivo pharmacokinetic parameters selected to reflect chronic exposure (metabolic and total clearances and area under the plasma concentration-time curve: CLm, CL and AUC) and acute exposure (the peak plasma concentration, C(max)). The variability in CYP1A2 activity in healthy adults, based on data after oral and intravenous dosage (CLm, CL and AUC), ranged from 34 to 42%. The variability in C(max) was 21%. The default kinetic factor of 3.16 would cover at least 99% of the healthy adult population, assuming that the data were log-normally distributed, but would give lower protection for some subgroups (pregnant women at term, healthy elderly, patients with liver disease), and was inadequate for neonates. This analysis of in vivo kinetic data for CYP1A2 substrates illustrates the importance of quantifying human variability in specific metabolic pathways, and of identifying potentially susceptible subgroups of the human population, in order to determine the scientific validity of uncertainty factors.

Drug interactions avoided-a useful indicator of good prescribing practice.

AIMS: To develop an index of quality prescribing in general practice by investigating the incidence of potential drug interactions when medicines were coprescribed within the State supported General Medical Services (GMS) in Ireland. METHODS: We determined an odds ratio (OR), as a measure of the relative risk of being exposed to a potential interaction, comparing the use of the H2-receptor antagonist, cimetidine, with that of the noninteracting agents ranitidine, famotidine and nizatidine in users and nonusers of warfarin, phenytoin and theophylline. RESULTS AND CONCLUSIONS: In 86 510 prescriptions for the H2-receptor antagonists potentially interacting drugs were dispensed to 8188 (9%) patients in the Eastern Health Board Region of the GMS. We found that prescribers were significantly less likely to use cimetidine (OR = 0.20,95% CI 0. 17-0.21, P < 0.001) in those patients who were coprescribed warfarin, suggesting good prescribing practice within the GMS. Similarly there was preferential use of the noninteracting H2-receptor antagonists in patients receiving phenytoin or theophylline and the extent of this selective prescribing was in keeping with the rank order of severity of interaction with these drugs. This novel pharmacological index may be a sensitive marker of good prescribing practice.