Pre-clinical pharmacokinetics of UK-453,061, a novel non-nucleoside reverse transcriptase inhibitor (NNRTI), and use of in silico physiologically based prediction tools to predict the oral pharmacokinetics of UK-453,061 in man.

1. UK-453,061 is a novel second-generation non-nucleoside reverse transcriptase inhibitor (NNRTI). Following intravenous bolus administration of UK-453,061 in male rat and infusion administration in dog, UK-453,061 had the following mean pharmacokinetic properties: elimination T(1/2) of 1.6 and 2.4 h, CL(p) of 26 and 10 ml min(-1) kg(-1) and V(ss) of 1.6 and 2 l kg(-1), respectively. 2. The half-lives of UK-453,061 disappearance in recombinant human CYPs 2C8, 2C9, 2A6, 2E1, 1A2, 2C19, 2D6 and 3A4 were 71, 100, 56, 101, 61, 34, 60 and 8 min, respectively. The disappearance half-life of UK-453,061 in human liver microsomes in the presence of UDPGA was 90 min. 3. Human clearance values were predicted using single-species scaling from in vivo data and from in vitro data using SimCYP. The human distribution of UK-453,061 was estimated using an in silico physiologically based pharmacokinetics (PBPK) methodology and absorption was predicted from measured physicochemical, permeability, and solubility data using GastroPlus and SimCYP. The C(max) was predicted to be 68, 185, 149% of the actual mean value using rat, dog and in vitro predicted values of human clearance at 30 mg and 53, 150, 29% of actual at 500 mg. The area under the curve (AUC) was predicted to be 73, 285 and 142% of the actual mean value using rat, dog and in vitro predicted values of human clearance at 30 mg and 52, 212 and 35% of actual at 500 mg. 4. This study demonstrates the utility of using in silico PBPK approaches to make predictions of human pharmacokinetics before dosing for the first time in humans.

Quantitative structure-activity relationships within a homologous series of 7-alkoxyresorufins exhibiting activity towards CYP1A and CYP2B enzymes: molecular modelling studies on key members of the resorufin series with CYP2C5-derived models of human CYP1A1, CYP1A2, CYP2B6 and CYP3A4.

1. The results of quantitative structure-activity relationships (QSARs) within a homologous series of 7-n-alkoxyresorufins are reported. They are consistent with homology modelling of the relevant P450s involved in their metabolism. 2. QSARs were generated for activities involving CYP1A and CYP2B enzymes with structural descriptors relating to compound planarity and other shape parameters, together with certain features of the n-alkoxyresorufin electronic structure, especially electron densities and superdelocalizabilities. 3. A quadratic relationship between compound lipophilicity and binding to CYP2B enzymes was apparent, and which indicated maximal interaction for 7-pentoxyresorufin. Such indications help to explain enzyme selectivity in terms of optimal alkyl chain length for fitting within the relevant active site region. 4. Calculation of the binding affinities for methoxy-, ethoxy-, pentoxy- and benzyloxy-resorufins towards either CYP1A2 or CYP2B6 enzymes, depending on the 7-alkoxyresorufin agree favourably with experimental values obtained from K(m) determinations.

Homology modelling of CYP3A4 from the CYP2C5 crystallographic template: analysis of typical CYP3A4 substrate interactions.

1. The results of homology modelling of cytochrome P4503A4 (CYP3A4), which is a human enzyme of major importance for the Phase 1 metabolism of drug substrates, from the CYP2C5 crystal structure is reported. 2. The overall homology between the two protein sequences was generally good (46%) with 24% of amino acid residues being identical and a 22% similarity between matched pairs in the CYP3A4 and CYP2C5 aligned sequences, thus indicating that CYP2C5 represents a viable template for modelling CYP3A4 by homology. 3. The CYP3A4 model appears to show consistency with the reported findings from the extensive site-directed mutagenesis studies already published. 4. Typical CYP3A4 substrates, such as midazolam, testosterone, nifedipine and verapamil, are shown to fit the putative active site of the enzyme structure in a manner consistent with their known positions of metabolism.

Homology modelling of the nuclear receptors: human oestrogen receptorbeta (hERbeta), the human pregnane-X-receptor (PXR), the Ah receptor (AhR) and the constitutive androstane receptor (CAR) ligand binding domains from the human oestrogen receptor alpha (hERalpha) crystal structure, and the human peroxisome proliferator activated receptor alpha (PPARalpha) ligand binding domain from the human PPARgamma crystal structure.

We have generated by homology the three-dimensional structures of the ligand binding domain (LBD) of several interrelated human steroid hormone receptors (SHRs). These are the oestrogen receptor beta (hERbeta), the pregnane-X-receptor (PXR), the Ah receptor (AhR) and the constitutive androstane receptor (CAR). They were produced by homology modelling from the human oestrogen receptor alpha (hERalpha) crystallographic coordinates [Nature 389 (1997) 753] as a template together with the amino acid sequences for hERbeta [FEBS Lett. 392 (1996) 49], PXR [J. Clin. Invest. 102 (1998) 1016], AhR [Proc. Natl. Acad. Sci. U.S.A. 89 (1992) 815] and CAR [Nature 395 (1998) 612; Mol. Cell. Biol. 14 (1994) 1544], respectively. The selective endogenous ligand, in each case, was docked interactively within the putative ligand binding site using the position of oestradiol in hERalpha as a guide, and the total energy was calculated. In each receptor model a number of different ligands known to fit closely within the ligand binding site were interactively docked and binding interactions noted. Specific binding interactions included combinations of hydrogen bonding and hydrophobic contacts with key amino acid sidechains, which varied depending on the nature of the ligand and receptor concerned. We also produced the human peroxisome proliferator activated receptor alpha (PPARalpha) by homology modelling using the human PPARgamma (hPPARgamma) LBD crystallographic coordinates summarised in [Toxicol. In Vitro 12 (1998) 619] as a template together with the amino acid sequence for hPPARalpha [Toxicol. In Vitro 12 (1998) 619; Nature 395 (1998) 137]. The models will provide a useful tool in unravelling the complexity in the physiologic response to xenobiotics by examining the ligand binding interactions and differences between the steroid hormone receptors activation or inactivation by their ligands.

Homology modelling of human CYP1A2 based on the CYP2C5 crystallographic template structure.

1. The results of homology modelling of human cytochrome P4501A2 (CYP1A2) based on the CYP2C5 crystal structure are reported. It exhibits improved sequence homology relative to that of CYP102. 2. It was demonstrated that many selective substrates for CYP1A2 could fit within the putative active site of the enzyme, and in orientations which agree with documented evidence for CYP1A2-mediated metabolism. 3. Furthermore, a number of amino acid residues lining the haem pocket have been shown, via site-directed mutagenesis, to have an influence on substrate metabolism, and these experimental findings from the literature are consistent with the modelled interactions for selective substrates. 4. The binding affinities of several CYP1A2 substrates have also been calculated from the CYP1A2 active site interactions and they agree closely with experimental values.

Homology modelling of CYP2A6 based on the CYP2C5 crystallographic template: enzyme-substrate interactions and QSARs for binding affinity and inhibition.

The results of homology modelling of the human P450 enzyme CYP2A6, based on the CYP2C5 crystallographic template structure are reported. A substantial number of selective substrates of the CYP2A6 enzyme fit the putative active site in a manner that is consistent with their known metabolites. Moreover, the evidence from site-directed mutagenesis experiments is in accordance with the current model, particularly in relation to complementary amino acid contacts within the haem environment. The binding of substrates is rationalized in terms of QSAR analyses and from a consideration of the contributory factors affecting the binding affinity. The latter approach appears to represent a highly correlated (R=0.99) method for estimating the relative strength of enzyme-substrate binding within CYP2A6-selective compounds, albeit within a fairly limited dataset of substrates.

Homology modelling of human CYP2E1 based on the CYP2C5 crystal structure: investigation of enzyme-substrate and enzyme-inhibitor interactions.

The construction of a homology model of human cytochrome P450 2E1 (CYP2E1) is reported, based on the CYP2C5 crystallographic template. A relatively high degree of primary sequence homology (identity=59%), as expected for proteins of the same CYP family, ensured a straightforward generation of the 3-dimensional model due to relatively few deletions and insertions of amino acid residues with respect to the CYP2C5 crystal structure. Probing the CYP2E1 model with typical substrates of the enzyme showed a good agreement with experimental information in the form of positions of metabolism for substrates, and with site-directed mutagenesis data on certain residues. Furthermore, quantitative relationships between substrate binding affinity and various structural parameters associated with the substrate molecules facilitated the formulation of a procedure for estimating relative binding energy and, consequently, K(m) or K(D) values towards the CYP2E1 enzyme. This method has been based on a consideration of the active site interactions between substrates and key amino acid residues lining the haem pocket, together with compound lipophilicity data from partition coefficients.

Diversity of Campylobacter isolates from retail poultry carcasses and from humans as demonstrated by pulsed-field gel electrophoresis.

Campylobacter spp. are a major contaminant of poultry. Eating undercooked chicken and handling raw poultry have been identified as risk factors for campylobacteriosis in humans. Previous studies have found Campylobacter spp. on 90% of poultry carcasses. In the present study, pulsed-field gel electrophoresis (PFGE) was used to assess the genetic diversity of strains on retail poultry carcasses. PFGE patterns of isolates from campylobacteriosis cases were compared to those from the poultry isolates. Over a 1-year study period (March 2000 through February 2001), whole fresh young chickens (n = 72) were obtained from three retail outlets in an urban community in the south-central United States. Campylobacter spp. were isolated from 82% of these carcasses. Strains (n = 70) were defined on the basis of their PFGE pattern. Sixty-seven percent of the carcasses from which Campylobacter spp. were isolated were contaminated with more than one PFGE-distinguishable strain. During the 1-year study period, most of the PFGE patterns (59%) were limited to isolates obtained from a single carcass. Forty-one percent of the PFGE-distinguishable strains were recovered from more than one carcass. Ninety-seven percent of the carcasses contaminated with the same strain were purchased at the same time from the same store. To examine the degree of genetic stability, four strains were followed in vitro over an estimated 1,000 doublings. The PFGE pattern of one of these isolates underwent minor changes during in vitro growth. The data indicate extensive variability in the PFGE patterns of Campylobacter spp. isolated from humans and from poultry carcasses. In spite of difficulties caused by such diversity and the fact that some carcasses are contaminated with more than one strain, the pattern variation provides a useful method for linking a particular strain to its source.

Quantitative structure--activity relationships for inducers of cytochromes P450 and nuclear receptor ligands involved in P450 regulation within the CYP1, CYP2, CYP3 and CYP4 families.

The results of quantitative structure-activity relationships (QSARs) are reported for several series of cytochrome P450 inducers, including those which also act as ligands for the various nuclear receptors involved in regulation of the relevant P450 genes, namely, CYP1, CYP2, CYP3 and CYP4. In several examples presented, the QSARs are consistent with homology modelling studies of the nuclear receptor ligand-binding domains (LBDs) based on available crystal structures of the oestrogen and peroxisome proliferator-activated receptors' LBDs. Good correlations (R=0.91-0.99) are found between various structural parameters and biological activity (either in the form of P450 induction or ligand-binding affinity) for the Ah receptor (AhR), human estrogen receptor alpha (hER alpha), human glucocorticoid receptor (hGR) and the rat peroxisome proliferator-activated receptor alpha (rPPAR alpha).

Molecular modelling of the peroxisome proliferator-activated receptor alpha (PPAR alpha) from human, rat and mouse, based on homology with the human PPAR gamma crystal structure.

The generation of homology models of human, rat and mouse peroxisome proliferator-activated receptor alpha (PPAR alpha) are reported, based on the recently published crystal structure of the human PPAR gamma ligand-binding domain (LBD) with bound ligand, rosiglitazone. It is found that a template of peroxisome proliferating fibrate drugs and related compounds can fit within the putative ligand-binding site of rat PPAR alpha, via contacts with amino acid residues which are consistent with their biological potency for peroxisome proliferation, site-directed mutagenesis experiments and with quantitative structure-activity relationship (QSAR) analysis studies. The experimental binding affinity of leukotriene B(4) (LTB(4)) for the mouse PPAR alpha agrees closely with the calculated value based on the modelled interactions, whereas selective PPAR alpha ligands such as clofibric acid are able to fit the human PPAR alpha binding site in agreement with reported site-directed mutagenesis information.

Factors influencing rates and clearance in P450-mediated reactions: QSARs for substrates of the xenobiotic-metabolizing hepatic microsomal P450s.

Various contributory factors associated with the kinetics of cytochrome P450-mediated catalytic activity and the metabolic clearance of drug substrates are discussed and evaluated, based on literature data and physicochemical parameters. Quantitative relationships between molecular structure and biological activity for several series of P450 substrates are presented which point to certain commonalities in P450-catalyzed reactions. In particular, it appears that frontier orbital energies are especially important for the estimation of reaction rates and clearance for many P450 substrates, although occasionally these have to be combined with other descriptors, such as compound lipophilicity (in the form of logP or logD(74)).

An interaction between the cytochrome P450 probe substrates chlorzoxazone (CYP2E1) and midazolam (CYP3A).

AIMS: The use of multiple probe substrates to evaluate the activity of drug metabolizing enzymes requires that there are no inter-substrate interactions. As part of a series of studies to develop a clinically useful collection of probe substrates that could be given alone or in any combination, we observed an interaction between midazolam (MDZ) and another component of the six-drug cocktail. Published data indicated that the interacting component was likely to be chlorzoxazone. This was investigated as part of a second study. The data relating to the interaction from both studies are reported here. METHODS: Both studies were performed in 16 healthy subjects. All treatments were given orally after an overnight fast. In study 1, which was performed to a four-period, open, crossover design, subjects received on separate occasions MDZ 5 mg, diclofenac 25 mg, a four drug cocktail (caffeine 100 mg, mephenytoin 100 mg, debrisoquine 10 mg and chlorzoxazone 250 mg) and a six drug cocktail (caffeine 100 mg, mephenytoin 100 mg, debrisoquine 10 mg, chlorzoxazone 250 mg, diclofenac 25 mg and MDZ 5 mg). In study 2, which was performed to a two-period, open, crossover design, subjects received a five drug cocktail (as the six drug cocktail in the first study, but without chlorzoxazone and with diclofenac dose increased to 50 mg) and a six drug cocktail (as five drug cocktail, with chlorzoxazone 250 mg). In both studies, blood samples were taken for measurement of plasma MDZ and 1-hydroxy MDZ (1-OH MDZ) concentrations. In study 1, blood samples were taken up to 12 h post-dose while in study 2 a single sample was taken 2 h after dosing. In study 1, the potential interaction between MDZ and the other components of the six drug cocktail was assessed by comparing AUClast ratios (1-OH MDZ/MDZ) between the two treatments. Additionally, a single sampling timepoint of 2 h post-dose for determination of concentration, rather than AUC, ratios was established. The 2 h plasma concentration ratios from studies 1 and 2 were combined and a pooled analysis performed to compare ratios within each study (to determine the change in ratio when MDZ was dosed with and without chlorzoxazone) and between studies (to determine the consistency of the ratios when MDZ was given either as part of the two six drug cocktails or when given alone and as part of the five drug cocktail). RESULTS: In study 1, both the AUClast ratio and the 2 h post-dose plasma concentration ratio were reduced when MDZ was given as part of the six drug cocktail in comparison with those for MDZ alone. This was the result of an increase in MDZ, rather than decrease in 1-OH MDZ, concentrations and was considered to result from a reduction in first pass metabolism of MDZ. The geometric mean AUClast values (with 95% CI) for MDZ were 95.6 (79.0, 115.7) and 160.4 (133.6, 192.6) microg l(-1) h when given alone and as part of the six drug cocktail, respectively. The corresponding values for 1-OH MDZ were 789.6 (697.6, 893.6) and 791.4 (701.7, 892.6) microg l(-1) h. The ratio of adjusted geometric mean AUClast ratios for the two treatments was 1.82 (90% CI 1.48, 2.23, P < 0.001). The pooled plasma 1-OH MDZ/MDZ ratio data from both studies showed that the differences in MDZ metabolism observed in study 1 were replicated in study 2. The adjusted geometric mean 1-OH MDZ/MDZ ratios when MDZ was given alone and as part of the six drug cocktail were 7.79 and 4.59, respectively, for study 1 (ratio 1.70, 95% CI 1.36, 2.11, P < 0.001) and 7.64 and 4.60 for study 2 (ratio 1.66, 95% CI 1.34, 2.06, P < 0.001). These data indicate that when given orally chlorzoxazone interacts with MDZ, increasing plasma MDZ concentrations. In contrast, there was no difference between the plasma 1-OH MDZ/MDZ ratios when MDZ was given alone and as part of the five drug cocktail indicating that there were no interactions between MDZ and any of the other components of that cocktail. CONCLUSIONS: Chlorzoxazone appears to significantly influence the pharmacokinetics of oral MDZ, probably through inhibition of first pass metabolism by CYP3A in the GI tract. Data from these studies and literature evidence showing a further interaction between chlorzoxazone and CYP1A2 substrates and questions concerning the specificity of chlorzoxazone as a probe substrate for CYP2E1, indicate that the use of chlorzoxazone in multisubstrate probe cocktails should be avoided.

Quantitative structure-activity relationships (QSARs) within series of inhibitors for mammalian cytochromes P450 (CYPs).

The results of quantitative structure-activity relationship (QSAR) studies on series of P450 inhibitors are reported. Cytochrome P450 families CYP1, CYP2 and CYP51 have been investigated for QSAR analysis, including those of CYP2 subfamilies: CYP2A, CYP2B, CYP2C, CYP2D and CYP2E. The accumulated evidence indicates different structural descriptors being involved, depending on the P450 enzyme concerned, although compound lipophilicity in the form of either logP or logD(7.4) appears to represent a common factor in some cases. This is thought to represent desolvation of the P450 active site, although quadratic expressions in lipophilicity tend to suggest that membrane transport is important, especially for CYP2B and CYP2E isoforms. In general, there is close agreement (R = 0.95-0.99) between experimental pKi values and those calculated via QSAR analysis.

Quantitative structure-activity relationships (QSARs) within substrates of human cytochromes P450 involved in drug metabolism.

The results of quantitative structure-activity relationship (QSAR) analyses are reported for structurally diverse series of chemicals which act as substrates or inhibitors for human hepatic microsomal cytochromes P450 (CYP). In particular, this study focuses on the major catalysts of drug metabolism in man, namely CYP1A2, CYP2B6, CYP2C9, CYP2C19, CYP2D6 and CYP3A4. It is found that good correlations (with correlation coefficients ranging from R = 0.94 to 0.99) with P450 binding affinity (Km and K(D)) or competitive inhibition (Ki) values are obtained in each case, especially when consideration of hydrogen bonding parameters are included in the QSAR analysis, together with the number of pi-pi stacking interactions.

Molecular modelling of human CYP2E1 by homology with the CYP102 haemoprotein domain: investigation of the interactions of substrates and inhibitors within the putative active site of the human CYP2E1 isoform.

1. The construction of a three-dimensional model of human CYP2E1 is reported. It is based on homology with the haemoprotein domain of the unusual bacterial P450, CYP102, which is of known crystal structure. 2. Interactive docking of a number of human CYP2E1 substrates is consistent with their known positions of CYP2E1-mediated metabolism, where specific interactions with key active site amino acid side-chains appear to rationalize the binding and orientation of substrate molecules. 3. Amino acid residues within the putative active site of human CYP2E1, including those associated with the binding of substrates and inhibitors, are shown to correspond with those identified by site-directed mutagenesis experiments conducted on CYP2 family isoforms, and they are known to affect substrate metabolism regioselectivity. 4. Consequently, it was found that the CYP2E1 active site exhibits complementarity with the structural characteristics of known substrates and inhibitors of this enzyme, including their relatively low molecular weights and disposition of hydrogen bond-forming groups.

Molecular modelling of CYP1 family enzymes CYP1A1, CYP1A2, CYP1A6 and CYP1B1 based on sequence homology with CYP102.

Molecular modelling of a number of CYP1 family enzymes from rat, plaice and human is described based on amino acid sequence homology with the haemoprotein domain of CYP102, a unique bacterial P450 of known structure. The interaction of various substrates and inhibitors within the putative active sites of rat CYP1A1, human CYP1A2, a fish CYP1 enzyme CYP1A6 (from plaice) and human CYP1B1, is shown to be consistent with P450-mediated oxidation in each example or, in the case of inhibitors, mechanism of inhibition. It is reported that relatively small changes between the enzymes' active site regions assist in the rationalization of CYP1 enzyme preferences for particular substrate types, and a template of superimposed CYP1A2 substrates is shown to fit the putative active site of the human CYP1A2 enzyme.

Molecular modelling of the human cytochrome P450 isoform CYP2A6 and investigations of CYP2A substrate selectivity.

(1) The generation of a homology model of CYP2A6, the major catalyst of human hepatic coumarin 7-hydroxylase activity, involves the use of the recently published substrate-bound CYP102 crystal structure as a template. (2) A substantial number of structurally diverse CYP2A6 substrates are found to dock satisfactorily within the putative active site of the enzyme, leading to the formulation of a structural template (or pharmacophore) for CYP2A6 specificity/selectivity. (3) The CYP2A6 model is consistent with available evidence from site-directed mutagenesis studies carried out on CYP2A subfamily isoforms, and enables some explanation of species differences in CYP2A-mediated metabolism of certain substrates. (4) Quantitative structure-activity relationship (QSAR) analysis of CYP2A5 (the mouse orthologue) mutants yields statistically significant correlations between various properties of amino acid residues and coumarin 7-hydroxylase activity.

A comparison of basal and induced hepatic microsomal cytochrome P450 monooxygenase activities in the cynomolgus monkey (Macaca fascicularis) and man.

1. The specific activities of hepatic microsomal cortisol 6beta-hydroxylase, coumarin 7-hydroxylase, S-mephenytoin 4'-hydroxylase and phenoxazone hydroxylase and the O-dealkylations of seven homologous alkoxyresorufins were < 3-fold different between the untreated (UT) cynomolgus monkey and man. 2. Heptoxy- and octoxyresorufin O-dealkylase, S-mephenytoin N-demethylase and dextromethorphan O-demethylase specific activities were > 6-fold higher, whereas tolbutamide hydroxylase was almost 5-fold lower in the UT monkey than in man. 3. Phenobarbitone induced (2-6-fold) coumarin 7-hydroxylase, cortisol 6beta-hydroxylase, S-mephenytoin N-demethylase, phenoxazone hydroxylase and benzyloxyresorufin O-dealkylase activities, but not the O-dealkylations of pentoxyresorufin or other alkoxyresorufins, in monkey. 4. Rifampicin induced (2-3-fold) cortisol 6beta-hydroxylase, S-mephenytoin 4'-hydroxylase, S-mephenytoin N-demethylase and tolbutamide hydroxylase activities, the O-dealkylations of methoxy-, ethoxy- and propoxyresorufin and CYP2C- and CYP3A-immunorelated proteins in monkey. 5. Dextromethorphan O-demethylase was significantly reduced by both phenobarbitone and rifampicin treatment in monkey. 6. Beta-naphthoflavone induced (8-39-fold) the O-dealkylations of several alkoxyresorufins, the greatest effect being on propoxyresorufin, but had no effect on the other activities measured in monkey. 7. Constitutive hepatic microsomal CYP2D6-immunorelated proteins were expressed at apparently much higher levels in monkey than in man.

Molecular modelling of CYP2B6, the human CYP2B isoform, by homology with the substrate-bound CYP102 crystal structure: evaluation of CYP2B6 substrate characteristics, the cytochrome b5 binding site and comparisons with CYP2B1 and CYP2B4.

1. Molecular modelling studies of CYP2B isoforms from rat (CYP2B1), rabbit (CYP2B4) and man (CYP2B6) are reported, with particular emphasis on substrate interactions with the human CYP2B isoform, CYP2B6. 2. The findings represent an advance on our previous study that focused primarily on the rat CYP2B isoform, CYP2B1, and involved homology modelling with substrate-free CYP102. 3. The current work utilizes the recently published substrate-bound CYP102 crystal structure as a template for construction of the CYP2B subfamily isoforms and shows, in particular, that known CYP2B6 substrate specificity and regioselectivity can be rationalized by putative active site interactions.

Advances in drug metabolism screening.

Developments in automation, analytical technologies and molecular biology are being exploited by drug metabolism scientists in order to provide enhanced in vitro systems for the study of the metabolic disposition of potential drug candidates. Routine investigation of factors such as metabolic stability and induction and inhibition of drug metabolizing enzymes is now preferred in the early stages of drug discovery. This, in turn, should provide a greater understanding of the underlying principles governing these processes and allow a greater role for drug metabolism in the design of new drug molecules.