Prediction of the pharmacokinetic parameters of reduced-dolasetron in man using in vitro-in vivo and interspecies allometric scaling

Discovery of SHR1653, a Highly Potent and Selective OTR Antagonist with Improved Blood-Brain Barrier Penetration

The oxytocin receptor (OTR) plays a major role in the control of male sexual responses. Antagonists of the OTR have been reported to inhibit ejaculation in animal models and serve as a potential treatment for premature ejaculation (PE). Herein, we describe a novel scaffold featuring an aryl substituted 3-azabicyclo [3.1.0] hexane structure. The lead compound, SHR1653, was shown to be a highly potent OTR antagonist, which exhibited excellent selectivity over V_1AR, V_1BR, and V₂R. This novel molecule was shown to have a favorable pharmacokinetic profile across species, as well as robust in vivo efficacy in a rat uterine contraction model. Interestingly, SHR1653 exhibited excellent blood-brain barrier penetration, which might be beneficial for the treatment of CNS-related PE.

Prediction of human pharmacokinetics—evaluation of methods for prediction of hepatic metabolic clearance

Methods for prediction of hepatic clearance (CLH) in man have been evaluated. A physiologically-based in-vitro to in-vivo (PB-IVIV) method with human unbound fraction in blood (fu,bl) and hepatocyte intrinsic clearance (CLint)-data has a good rationale and appears to give the best predictions (maximum ∼2-fold errors; < 25% errors for half of CL-predictions; appropriate ranking). Inclusion of an empirical scaling factor is, however, needed, and reasons include the use of cryopreserved hepatocytes with low activity, and inappropriate CLint- and fu,bl-estimation methods. Thus, an improvement of this methodology is possible and required. Neglect of fu,bl or incorporation of incubation binding does not seem appropriate. When microsome CLint-data are used with this approach, the CLH is underpredicted by 5- to 9-fold on average, and a 106-fold underprediction (attrition potential) has been observed. The poor performance could probably be related to permeation, binding and low metabolic activity. Inclusion of scaling factors and neglect of fu,bl for basic and neutral compounds improve microsome predictions. The performance is, however, still not satisfactory. Allometry incorrectly assumes that the determinants for CLH relate to body weight and overpredicts human liver blood flow rate. Consequently, allometric methods have poor predictability. Simple allometry has an average overprediction potential, > 2-fold errors for ∼1/3 of predictions, and 140-fold underprediction to 5800-fold overprediction (potential safety risk) range. In-silico methodologies are available, but these need further development. Acceptable prediction errors for compounds with low and high CLH should be ∼50 and ∼10%, respectively. In conclusion, it is recommended that PB-IVIV with human hepatocyte CLint and fu,bl is applied and improved, limits for acceptable errors are decreased, and that animal CLH-studies and allometry are avoided.

Physiologically based pharmacokinetics (PBPK)

Allometric scaling is widely used to predict human pharmacokinetic parameters from preclinical species, and many different approaches have been proposed over the years to improve its predictive performance. Nevertheless, prediction errors are commonly observed in the practical application of simple allometry, for example, in cases where the hepatic metabolic clearance is mainly determined by enzyme activities, which do not scale allometrically across species. Therefore, if good correlation was noted for some drugs, poor correlation was observed for others, highlighting the need for other conceptual approaches. Physiologically based pharmacokinetic (PBPK) models are now a well-established approach to conduct extrapolations across species and to generate simulations of pharmacokinetic profiles under various physiological conditions. While conventional pharmacokinetic models are defined by drug-related data themselves, PBPK models have richer information content and integrate information from various sources, including drug-dependent, physiological, and biological parameters as they vary in between species, subjects, or with age and disease state. Therefore, the biological and mechanistic bases of PBPK models allow the extrapolation of the kinetic behavior of drugs with regard to dose, route, and species. In addition, by providing a link between tissue concentrations and toxicological or pharmacological effects, PBPK modeling represents a framework for mechanistic pharmacokinetic-pharmacodynamic models.

A global examination of allometric scaling for predicting human drug clearance and the prediction of large vertical allometry**This work was presented at the American Association of Pharmaceutical Scientists Annual meeting, Salt Lake City, USA, Oct. 26, 2003

Allometrically scaled data sets (138 compounds) used for predicting human clearance were obtained from the literature. Our analyses of these data have led to four observations. (1) The current data do not provide strong evidence that systemic clearance (CL(s); n = 102) is more predictable than apparent oral clearance (CL(po); n = 24), but caution needs to be applied because of potential CL(po) prediction error caused by differences in bioavailability across species. (2) CL(s) of proteins (n = 10) can be more accurately predicted than that of non-protein chemicals (n = 102). (3) CL(s) is more predictable for compounds eliminated by renal or biliary excretion (n = 33) than by metabolism (n = 57). (4) CL(s) predictability for hepatically eliminated compounds followed the order: high CL (n = 11) > intermediate CL (n = 17) > low CL (n = 29). All examples of large vertical allometry (% error of prediction greater than 1000%) occurred only when predicting human CL(s) of drugs having very low CL(s). A qualitative analysis revealed the application of two potential rules for predicting the occurrence of large vertical allometry: (1) ratio of unbound fraction of drug in plasma (f(u)) between rats and humans greater than 5; (2) C logP greater than 2. Metabolic elimination could also serve as an additional indicator for expecting large vertical allometry.

A Novel Strategy for Physiologically Based Predictions of Human Pharmacokinetics

BACKGROUND

The major aim of this study was to develop a strategy for predicting human pharmacokinetics using physiologically based pharmacokinetic (PBPK) modelling. This was compared with allometry (of plasma concentration-time profiles using the Dedrick approach), in order to determine the best approaches and strategies for the prediction of human pharmacokinetics.

METHODS

PBPK and Dedrick predictions were made for 19 F. Hoffmann-La Roche compounds. A strategy for the prediction of human pharmacokinetics using PBPK modelling was proposed in this study. Predicted values (pharmacokinetic parameters, plasma concentrations) were compared with observed values obtained after intravenous and oral administration in order to assess the accuracy of the prediction methods.

RESULTS

By following the proposed strategy for PBPK, a prediction would have been made prospectively for approximately 70% of the compounds. The prediction accuracy for these compounds in terms of the percentage of compounds with an average-fold error of <2-fold was 83%, 50%, 75%, 67%, 92% and 100% for apparent oral clearance (CL/F), apparent volume of distribution during terminal phase after oral administration (V(z)/F), terminal elimination half-life (t(1/2)), peak plasma concentration (C(max)), area under the plasma concentration-time curve (AUC) and time to reach C(max) (t(max)), respectively. For the other 30% compounds, unacceptable prediction accuracy was obtained in animals; therefore, a prospective prediction of human pharmacokinetics would not have been made using PBPK. For these compounds, prediction accuracy was also poor using the Dedrick approach. In the majority of cases, PBPK gave more accurate predictions of pharmacokinetic parameters and plasma concentration-time profiles than the Dedrick approach.

CONCLUSIONS

Based on the dataset evaluated in this study, PBPK gave reasonable predictions of human pharmacokinetics using preclinical data and is the recommended approach in the majority of cases. In addition, PBPK modelling is a useful tool to gain insights into the properties of a compound. Thus, PBPK can guide experimental efforts to obtain the relevant information necessary to understand the compound's properties before entry into human, ultimately resulting in a higher level of prediction accuracy.

A novel model for prediction of human drug clearance by allometric scaling

Sixty-one sets of clearance (CL) values in animal species were allometrically scaled for predicting human clearance. Unbound fractions (f(u)) of drug in plasma in rats and humans were obtained from the literature. A model was developed to predict human CL: CL=33.35 ml/min x (a/Rf(u))(0.770), where Rf(u) is the f(u) ratio between rats and humans and a is the coefficient obtained from allometric scaling. The new model was compared with simple allometric scaling and the "rule of exponents" (ROE). Results indicated that the new model provided better predictability for human values of CL than did ROE. It is especially significant that for the first time the proposed model improves the prediction of CL for drugs illustrating large vertical allometry.

Allometric principles for interspecies extrapolation in toxicological risk assessment--empirical investigations

Four types of data (toxicokinetic data of pharmaceuticals from six species including humans, LD(50) values from eight animal species, long-term NOAEL values of pesticides from mice, rats, and dogs, and toxicity data on anti-neoplastic agents from six species including humans) were used for interspecies comparisons. Species differences with regard to kinetic parameters and toxicity were evaluated and the concordance with predictions by allometric scaling according to caloric demand (allometric exponent 0.75) or to body weight (allometric exponent 1) was checked. For LD(50) values, agreement was poor for both allometric concepts. Recently reported concordance of LD(50) species differences with body weight scaling could be traced back to biased data selection. The other three datasets are clearly in agreement with the allometric scaling according to caloric demand. Caloric demand scaling is thus proposed as a generic interspecies extrapolation method in the absence of substance-specific data. Moreover, the evaluated data make it possible to describe uncertainty associated with the process of interspecies extrapolation by allometric rules.

Physiologically based pharmacokinetic (PBPK) modeling of disposition of epiroprim in humans

The objective of this study was to use in synergy physiologically based and empirical approaches to estimate the drug-specific input parameters of PBPK models of disposition to simulate the plasma concentration-time profile of epiroprim in human. The estimated input parameters were the tissue:plasma partition coefficients (Pt:p) for distribution and the blood clearance (CL) for the in vivo conditions. Epiroprim represents a challenge for such methods, because it shows large interspecies differences in its pharmacokinetic properties. Two approaches were used to predict the human Pt:p values: the tissue composition model (TCM) and the "Arundel approach" based on the volume of distribution at steady state (Vdss) determined in vivo in the rat. CL in human was predicted by (1) conventional allometric scaling of in vivo animal clearances (CAS), (2) physiologically based direct scaling up of in vitro hepatocyte data (DSU), and (3) allometric scaling of animal intrinsic in vivo blood CL normalized by the ratios of animal:human intrinsic clearances determined in vitro with hepatocytes (NAS). The performance of prediction was assessed by comparing separately the above pharmacokinetic parameters (Vdss estimated from the Pt:p values and blood CL) with the corresponding in vivo data obtained from the plasma kinetic profiles. These input parameters were used in PBPK models, and the resulting plasma concentration-time profiles of epiroprim were compared with those observed in rat and human. Previously to the construction of the human PBPK model, a model for the rat was also developed to gain more confidence on the model structure and assumptions. Overall, using the TCM and the NAS for the parameterization of the distribution and clearance, respectively, the PBPK model gave the more accurate predictions of epiroprim's disposition in human. This study represents therefore an attractive approach, which may potentially help the clinical candidate selection.

Metabolism and disposition of novel des-fluoro quinolone garenoxacin in experimental animals and an interspecies scaling of pharmacokinetic parameters

Garenoxacin is a novel quinolone that does not have a fluorine substituent at the C-6 position in the quinoline ring. Garenoxacin or 14C-garenoxacin was intravenously or orally administered to rats, dogs, and monkeys. Metabolic profiles and pharmacokinetic parameters were investigated focusing on the species differences and the allometric scaling of pharmacokinetic parameters. Garenoxacin was well absorbed following oral administration then underwent phase II metabolism in all species tested. Major metabolites of garenoxacin were the sulfate of garenoxacin (M1) and glucuronide (M6). Oxidative metabolites were present in very minor concentrations in all species tested. Another minor route of metabolism was the formation of the carbamoyl glucuronide. Garenoxacin is characterized across species by the observation that it circulates systemically, is excreted renally as unchanged drug, and is metabolized to M1 and M6, which are excreted specifically into the bile. The total clearances (CL) were 12.1, 2.43, and 3.39 ml/min/kg for rats, dogs, and monkeys, respectively. The distribution volume values of garenoxacin (Vss) were 0.88, 1.29, and 0.96 l/kg for rats, dogs, and monkeys, respectively. In all animals tested, the extrarenal clearance was larger than the renal clearance, and neither of the clearances was limited by blood flow. Despite these conditions, garenoxacin showed a good correlation for CL and Vss for allometric interspecies scaling.

Allometric scaling of xenobiotic clearance: uncertainty versus universality

Statistical analysis and Monte Carlo simulation were used to characterize uncertainty in the allometric exponent (b) of xenobiotic clearance (CL). CL values for 115 xenobiotics were from published studies in which at least 3 species were used for the purpose of interspecies comparison of pharmacokinetics. The b value for each xenobiotic was calculated along with its confidence interval (CI). For 24 xenobiotics (21%), there was no correlation between log CL and log body weight. For the other 91 cases, the mean +/- standard deviation of the b values was 0.74 +/- 0.16; range: 0.29 to 1.2. Most (81%) of these individual b values did not differ from either 0.67 or 0.75 at P = 0.05. When CL values for the subset of 91 substances were normalized to a common body weight coefficient (a), the b value for the 460 adjusted CL values was 0.74; the 99% CI was 0.71 to 0.76, which excluded 0.67. Monte Carlo simulation indicated that the wide range of observed b values could have resulted from random variability in CL values determined in a limited number of species, even though the underlying b value was 0.75. From the normalized CL values, four xenobiotic subgroups were examined: those that were (i) protein, and those that were (ii) eliminated mainly by renal excretion, (iii) by metabolism, or (iv) by renal excretion and metabolism combined. All subgroups except (ii) showed a b value not different from 0.75. The b value for the renal excretion subgroup (21 xenobiotics, 105 CL values) was 0.65, which differed from 0.75 but not from 0.67.

Utility of hepatocytes to model species differences in the metabolism of loxtidine and to predict pharmacokinetic parameters in rat, dog and man

1. The metabolism of loxtidine (1-methyl-5-[3-[3-[(1-piperidinyl) methyl] phenoxy] propyl] amino-1H-1,2,4-triazole-3-methanol) was studied in freshly isolated rat, dog and human hepatocytes. Metabolism in vitro was comparable with previously available in vivo data in all three species with the marked species differences observed in vivo being reproduced in the hepatocyte model. 2. The major route for the metabolism of loxtidine by rat hepatocytes was N-dealkylation to form the propionic acid and hydroxymethyl triazole metabolites. A minor metabolic route was the oxidation of loxtidine to a carboxylic acid metabolite. The major route of metabolism for loxtidine in dog hepatocytes was glucuronidation with oxidation to the carboxylic acid metabolite being of minor importance. Incubation of loxtidine with human hepatocytes resulted in the drug remaining largely unchanged but with the carboxylic acid metabolite being produced in minor amounts. 3. In vitro studies were undertaken with rat, dog and human hepatocytes to determine the Michaelis-Menten parameters Vmax and Km for the sum of all the metabolic pathways. These kinetic parameters were used to calculate the intrinsic clearance of loxtidine. Using appropriate scaling factors, the predicted in vivo hepatic clearance was then calculated. The predicted intrinsic clearances were 51.4 +/- 12.4, 8.0 +/- 0.8 and 1.0 +/- 0.6 ml/min/kg for rat, dog and human hepatocytes respectively. These data were then used to calculate hepatic clearances of 24.5, 3.1 and 0.2 ml/min/kg for rat, dog and man respectively. 4. In vivo hepatic and intrinsic clearances for loxtidine were determined in rat, dog and human volunteers. The hepatic clearances of loxtidine were 26.6, 6.6 and 0.4 ml/min/kg in rat, dog and man respectively and intrinsic clearances were 58.5, 18.6 and 2.0 ml/min/kg in rat, dog and man respectively. 5. The present studies demonstrate that the hepatocyte model may be a valuable in vitro tool for predicting both qualitative and quantitative aspects of the metabolism of a drug in animals and man at an early stage of the drug development process.