High confidence predictions of drug-drug interactions: predicting affinities for cytochrome P450 2C9 with multiple computational methods

Four different models are used to predict whether a compound will bind to 2C9 with a K(i) value of less than 10 microM. A training set of 276 compounds and a diverse validation set of 50 compounds were used to build and assess each model. The modeling methods are chosen to exploit the differences in how training sets are used to develop the predictive models. Two of the four methods develop partitioning trees based on global descriptions of structure using nine descriptors. A third method uses the same descriptors to develop local descriptions that relate activity to structures with similar descriptor characteristics. The fourth method uses a graph-theoretic approach to predict activity based on molecular structure. When all of these methods agree, the predictive accuracy is 94%. An external validation set of 11 compounds gives a predictive accuracy of 91% when all methods agree.

Use of Density Functional Calculations To Predict the Regioselectivity of Drugs and Molecules Metabolized by Aldehyde Oxidase

Aldehyde oxidase is a molybdenum hydroxylase that catalyzes the oxidation of aldehydes and nitrogen-containing heterocycles. The enzyme plays a dual role in the metabolism of physiologically important endogenous compounds and the biotransformation of xenobiotics. Using density functional theory methods, geometry optimization of tetrahedral intermediates of drugs and druglike compounds was examined to predict the likely metabolites of aldehyde oxidase. The calculations suggest that the lowest energy tetrahedral intermediate resulting from the initial substrate corresponds to the observed metabolite >or=90% of the time. Additional calculations were performed on a series of heterocyclic compounds where the products resulting from metabolism by xanthine oxidase and aldehyde oxidase differ in many instances. Again, the lowest energy tetrahedral intermediate corresponded to the observed product of aldehyde oxidase metabolism >or=90% for the compounds examined, while the observed products of xanthine oxidase were not well predicted.

Violations of pesticide use and worker safety regulations in North Carolina

In North Carolina, responsibility for providing training and enforcing various regulations related to pesticide use and agricultural worker safety rests with three state agencies. This article summarizes an 11-year history of enforcement procedures concerning agricultural pesticide use, the Worker Protection Standard, the Hazard Communication Standard, the Migrant Housing Act of North Carolina, and field sanitation standards. The difficulty of linking specific types of violations with worker safety is discussed.

Formation of the Active Species of Cytochrome P450 by Using Iodosylbenzene: A Case for Spin-Selective Reactivity

The generation of the active species for the enzyme cytochrome P450 by using the highly versatile oxygen surrogate iodosylbenzene (PhIO) often produces different results compared with the native route, in which the active species is generated through O(2) uptake and reduction by NADPH. One of these differences that is addressed here is the deuterium kinetic isotope effect (KIE) jump observed during N-dealkylation of N,N-dimethylaniline (DMA) by P450, when the reaction conditions change from the native to the PhIO route. The paper presents a theoretical analysis targeted to elucidate the mechanism of the reaction of PhIO with heme, to form the high-valent iron-oxo species Compound I (Cpd I), and define the origins of the KIE jump in the reaction of Cpd I with DMA. It is concluded that the likely origin of the KIE jump is the spin-selective chemistry of the enzyme cytochrome P450 under different preparation procedures. In the native route, the reaction proceeds via the doublet spin state of Cpd I and leads to a low KIE value. PhIO, however, diverts the reaction to the quartet spin state of Cpd I, which leads to the observed high KIE values. The KIE jump is reproduced here experimentally for the dealkylation of N,N-dimethyl-4-(methylthio)aniline, by using intra-molecular KIE measurements that avoid kinetic complexities. The effect of PhIO is compared with N,N-dimethylaniline-N-oxide (DMAO), which acts both as the oxygen donor and the substrate and leads to the same KIE values as the native route.

CYP2C9 protein interactions with cytochrome b(5): effects on the coupling of catalysis

The hemoprotein cytochrome b(5) (cyt b5) has been demonstrated to affect the kinetics of drug oxidation by the microsomal cytochromes P450 (P450s). However, the mechanisms through which cyt b5 exerts these effects are variable and P450 isoform-dependent. Whereas the effects of cyt b5 on the major drug-metabolizing enzymes CYP2D6, CYP2E1, and CYP3A4 are well studied, fewer studies conducted over limited ranges of cyt b5 concentrations have been performed on CYP2C9. In the present study with CYP2C9, cyt b5 exerted complex actions upon P450 oxidative reactions by affecting the rate of metabolite formation, the consumption of NADPH by cytochrome P450 reductase, and uncoupling of the reaction cycle to hydrogen peroxide and water. Cytochrome b(5) devoid of the heme moiety (apo-b5) exhibited effects similar to those of native cyt b5. All rates were highly dependent on the cyt b5 to CYP2C9 enzyme ratio, suggesting that the amount of cyt b5 present in an in vitro incubation is an important factor that can have an impact on the reliability of extrapolating in vitro generated data to predict the in vivo condition. The major effects of cyt b5 are hypothesized to result from a cyt b5-induced conformational change in CYP2C9 that results in an increased collision frequency between the iron-oxygen species (Cpd I) and the substrate, and a decrease in the oxidase activity. Together, these findings suggest that cyt b5 can alter multiple steps in the P450 catalytic cycle via complex interactions with P450 and P450 reductase.

Line-walking method for predicting the inhibition of P450 drug metabolism

A new method, called line-walking recursive partitioning (LWRP), for partitioning diverse structures on the basis of chemical properties that uses only nine descriptors of the shape, polarizability, and charge of the molecule is described. We use a training set of over 600 compounds and a validation set of 100 compounds for the cytochrome P450 enzymes 2C9, 2D6, and 3A4. The LWRP algorithm itself incorporates elements from support vector machines (SVMs) and recursive partitioning, while circumventing the need for the linear or quadratic programming methods required in SVMs. We compare LWRP with a many-descriptor SVM model, using the same dataset as that described in the literature.(1) The line-walking method, using nine descriptors, predicted the validation set with about 84-90% accuracy, a success rate comparable to that of the SVM method. Furthermore, line-walking was able to find errors in the assignment of inhibitor values within the validation set for the 2C9 inhibitors. When these errors are corrected, the model predicts with an even higher level of accuracy. Although this method has been applied to P450 enzymes, it should be of general use in partitioning molecules on the basis of function.

KINETIC ISOTOPE EFFECTS IMPLICATE A SINGLE OXIDANT FOR CYTOCHROME P450-MEDIATEDO-DEALKYLATION,N-OXYGENATION, AND AROMATIC HYDROXYLATION OF 6-METHOXYQUINOLINE

One major point of controversy in the area of cytochrome P450 (P450)-mediated oxidation reactions is the nature of the active-oxygen species. A number of hypotheses have been advanced which implicate a second oxidant besides the iron-oxo species designated as compound I (Cpd 1). This oxygen is thought to be either an iron-hydroperoxy species (Cpd 0) or a second spin-state of Cpd 1. Very little information is available on what fraction of P450 oxidations is mediated by the two different oxidants. Herein, we report results on three cytochrome P450-mediated reactions: O-dealkylation, N-oxygenation, and aromatic hydroxylation, which occur by three distinct chemical mechanisms. We have used kinetic isotope effects to test for branching from O-demethylation to N-oxygenation and aromatic hydroxylation, using 6-methoxyquinoline and 2H3-6-methoxyquinoline as substrates for P4501A2. Identical large inverse isotope effects on Vmax/Km are obtained for the formation of both the N-oxide and the phenol. This indicates that all three reactions occur through the same enzyme-substrate complex and, thus, through a single iron-oxygen species. The nature of the iron-oxygen species is less certain but is more likely to be iron-oxo Cpd 1, given the energetics of these reactions.

Clinical and toxicological relevance of CYP2C9: drug-drug interactions and pharmacogenetics

CYP2C9 is a major cytochrome P450 enzyme that is involved in the metabolic clearance of a wide variety of therapeutic agents, including nonsteroidal antiinflammatories, oral anticoagulants, and oral hypoglycemics. Disruption of CYP2C9 activity by metabolic inhibition or pharmacogenetic variability underlies many of the adverse drug reactions that are associated with the enzyme. CYP2C9 is also the first human P450 to be crystallized, and the structural basis for its substrate and inhibitor selectivity is becoming increasingly clear. New, ultrapotent inhibitors of CYP2C9 have been synthesised that aid in the development of quantitative structure-activity relationship (QSAR) models to facilitate drug redesign, and extensive resequencing of the gene and studies of its regulation will undoubtedly help us understand interindividual variability in drug response and toxicity controlled by this enzyme.

Charge and Substituent Effects on Affinity and Metabolism of Benzbromarone-Based CYP2C19 Inhibitors

Human cytochrome P450 (CYP) 2C19 is one of the most important CYP2C family members responsible for metabolizing commonly prescribed drugs. This research describes synthetic modifications to benzbromarone (Bzbr) to create the most potent CYP2C19 inhibitor ever reported. The most important features enabling analogues of Bzbr to bind to CYP2C19 with high affinity are low acidity (high pK(a) or nonionizability) and hydrophobic substituents adjacent to the phenol moiety. Though CYP2C19 was known to prefer neutral substrates, the extent was perhaps not realized until the anionic, parent compound Bzbr (K(i) = 3.7 microM) was compared to a less acidic dimethyl analogue (K(i) = 0.033 microM). However, differences in affinity for anionic and neutral Bzbr analogues did not appear to affect the regiospecificity of their metabolism by CYP's 2C19 and 2C9. In addition, some Bzbr analogues were metabolized both on the phenol and benzofuran rings. By using a substrate with a methyl ether in place of the Bzbr phenol, it was shown that some Bzbr analogues must be able to freely reposition after binding and oxidize the more energetically favorable position. Normally, O-demethylation of this methyl ether is favored over benzofuran hydroxylation based on ion current from LC/MS. Deuterium substitution of the methyl ether results in an inverse isotope effect on benzofuran hydroxylation (i.e. increased oxidation of this less favorable site). Likewise, Bzbr-based CoMFA models of CYP2C19 demonstrated no clear preference for any one ligand alignment. This suggests results from this modeling method must be interpreted carefully for each CYP isoform. In summary, Bzbr analogues have demonstrated they can be adapted to other CYP2C enzymes in order to probe isoform-specific properties.

Pathogenicity of Leptographium Species Associated with Loblolly Pine Decline

Freshly lifted seedlings and 21-year-old trees of loblolly pine were wound-inoculated with Leptographium species recovered from the soil and/or roots of trees with loblolly decline symptoms in central Alabama. Seedlings inoculated with L. procerum in the greenhouse produced significantly fewer root initials and a smaller root mass than control seedlings. Vertical lesions produced in seedlings by L. serpens and L. terebrantis were significantly longer than in controls. Lesions produced in mature trees by L. serpens and L. lundbergii were significantly longer than in controls. Of the fungi tested, L. serpens, L. terebrantis, and L. lundbergii were the most aggressive and may pose the greatest threat to loblolly pines.

Metabolism and mutagenicity of source water contaminants 1,3-dichloropropane and 2,2-dichloropropane

Cytochrome P450-dependent oxidation and glutathione (GSH)-dependent conjugation are the primary routes of metabolism of haloalkanes. Using rat liver microsomes and cytosol, we investigated the metabolism of two halopropanes found on the U.S. Environmental Protection Agency Contaminant Candidate List, 1,3-dichloropropane (1,3-DCP) and 2,2-dichloropropane (2,2-DCP). An automated headspace technique using gas chromatography was developed to determine rates of metabolism. Additional dihaloalkanes (1,2-dichloroethane, 1,2-dichloropropane, 1,4-dichlorobutane, 1,2-dibromoethane, 1,2-dibromopropane, 1,4-dibromobutane) were evaluated to assess structure-activity relationships. In general, brominated dihaloalkanes were eliminated from rat cytosol faster than chlorinated dihaloalkanes, reflecting the expected halide order of reactivity (Br > Cl). Furthermore, the rate of GSH conjugation was proportional to alpha,omega-haloalkane chain length. The clearance of 1,3-DCP via the GSH conjugation pathway (1.6 x 10(-4) l/h/mg cytosol protein) was minor relative to the P450 pathway (2.8 x 10(-2) l/h/mg microsomal protein). In contrast, we did not observe metabolism of 2,2-DCP via the GSH-dependent conjugation pathway and observed only a minor level of clearance via the P450 pathway (7 x 10(-4) l/h/mg microsomal protein). Neither compound was mutagenic in various strains of Salmonella, including those containing GSTT1-1, indicating that GSTT1-1 does not metabolize 1,3-DCP or 2,2-DCP to mutagens. Analysis of the reaction products of 1,3-DCP and GSH in cytosol by liquid chromatography/mass spectrometry revealed significant production of S-(3-chloropropyl) glutathione conjugate, indicating that the conjugate half-mustard does not rearrange to form a sulfonium ion, as typically occurs with vicinal dihaloalkanes.

Active-site characteristics of CYP2C19 and CYP2C9 probed with hydantoin and barbiturate inhibitors

Three series of N-3 alkyl substituted phenytoin, nirvanol, and barbiturate derivatives were synthesized and their inhibitor potencies were tested against recombinant CYP2C19 and CYP2C9 to probe the interaction of these ligands with the active sites of these enzymes. All compounds were found to be competitive inhibitors of both enzymes, although the degree of inhibitory potency was generally much greater towards CYP2C19. Inhibitor stereochemistry did not markedly influence K(i) towards CYP2C9, and log P adequately predicted inhibitor potency for this enzyme. In contrast, stereochemistry was an important factor in determining inhibitor potency towards CYP2C19. (S)-(+)-N-3-Benzylnirvanol and (R)-(-)-N-3-benzylphenobarbital emerged as the most potent and selective CYP2C19 inhibitors, with K(i) values of < 250nM--at least two orders of magnitude greater inhibitor potency than towards CYP2C9. Both inhibitors were metabolized preferentially at their C-5 phenyl substituents, indicating that CYP2C19 prefers to orient the N-3 substituents away from the active oxygen species. These features were incorporated into expanded CoMFA models for CYP2C9, and a new, validated CoMFA model for CYP2C19.

Kinetic Isotope Effects Implicate the Iron−Oxene as the Sole Oxidant in P450-Catalyzed N-Dealkylation

Multiple oxidants have been implicated as playing a role in cytochrome P450-mediated oxidations. Herein, we report results on N-dealkylation, one of the most facile reactions mediated by P450 enzymes. We have employed the N-oxides of a series of para-substituted 13C2H2-labeled N,N-dimethylanilines to function as both substrates and surrogate oxygen atom donors for P450cam and P4502E1. Kinetic isotope effect profiles obtained using the N-oxide system were found to closely match the profiles produced using the complete NAD(P)H/NAD(P)-P450 reductase/O2 system. The results are consistent with oxidation occurring solely through an iron-oxene species.

Quantitative Binding Models for CYP2C9 Based on Benzbromarone Analogues†

The cytochrome P450 (CYP) isoforms involved in xenobiotic metabolism are enzymes whose substrate selectivity remains difficult to predict due to wide specificity and dynamic protein-substrate interactions. To uncover the determinants of specificity for cytochrome CYP2C9, a novel library of benzbromarone (bzbr) inhibitors was used to reevaluate its pharmacophore. CoMSIA was used with the bzbr ligands to generate both quantitative binding models and three-dimensional contour plots that pinpoint predicted interactions that are important for binding to 2C9. Since this class of compounds is more potent than any other toward 2C9, the small molecule properties deemed most ideal by the software were used to address protein-ligand interactions using new mutagenesis and structural data. Nine new bzbr analogues provide evidence that specific electrostatic and hydrophobic interactions contribute the most to 2C9's specificity. Three of the new analogues are better isosteres of bzbr that contain bulky groups adjacent to the phenol and have increased pK(a) values. These ligands test the hypothesis that anionic substrates bind with higher affinity to 2C9. Since they have higher affinity than the previous nonacidic analogues, the importance of bulky groups on the phenol ring appears to have been underestimated. CoMSIA models predict that these bulky groups are favorable for their hydrophobicity, while a negative charge is favored at the ketone oxygen rather than the phenol oxygen. The overlap of this ketone with electronegative groups of other 2C9 substrates suggests they act as key positive charge acceptors.

Differential Roles of Arg97, Asp293, and Arg108 in Enzyme Stability and Substrate Specificity of CYP2C9

CYP2C9 metabolizes a wide range of drugs, many of which are negatively charged at physiological pH. Therefore, it has been thought that complementarily charged amino acid(s) are critically involved in substrate binding. Previous studies have implicated arginine residues at positions 97, 105, and 108 and aspartate at position 293 in the normal catalytic function of the enzyme. To elucidate the role of these amino acids in the substrate specificity of CYP2C9, a series of mutants were constructed and analyzed for functional activity, thermal stability, and ligand binding. Charge-modifying mutations at positions 97, 105, and 293 decreased catalytic activity toward diclofenac, (S)-warfarin, and pyrene in a substrate-independent manner with Arg105 the least, and Arg97 the most, sensitive amino acids in this regard. Decreases in functional activity paralleled thermal instability of the mutants, suggesting that loss of function reflects more generalized structural changes rather than the absence of a specifically charged amino acid at these three positions. The R108H mutant was inactive toward all three substrates because of unexpected nitrogen ligation to the heme. Conversely, the R108F mutant exhibited substrate-dependent catalytic behavior, with almost complete loss of activity toward (S)-warfarin and diclofenac, but preservation of pyrene metabolism. In addition, the R108F mutation abrogated the Type I difference spectra induced by flurbiprofen and benzbromarone, obligate anions at physiological pH. These data identify critical roles for Arg97 and Asp293 in the structural stability of the enzyme and demonstrate a selective role for Arg108 in the binding and metabolism of negatively charged substrates of CYP2C9.

Crystallographic Studies on the Complex Behavior of Nicotine Binding to P450cam (CYP101)†

Crystallographic and spectroscopic studies have been undertaken to characterize the binding behavior of the non-native substrate nicotine in the active site of the monooxygenase hemoprotein cytochrome P450cam. Despite the existence of a theoretical model that is consistent with the observed distribution of monooxygenation products, the crystal structure of the complex indicates that the primary binding mode of nicotine is unproductive. The structure is confirmed by spectral data that indicate direct coordination of substrate pyridine nitrogen with the heme iron. This would be the proper structure for evaluating binding affinity and inhibition. Reduction of the heme from Fe(III) to Fe(II) and introduction of carbon monoxide into crystals of the nicotine-P450cam complex, to simulate molecular oxygen binding, produces reorientation of the nicotine. This orientation is the appropriate one for predicting regioselectivity and the kinetic features of substrate oxidation. While it is not clear that such complicated behavior will be exhibited for other enzyme-substrate interactions, it is clear that a single crystal structure for a given substrate-enzyme interaction may not provide a good description of the binding mode responsible for product formation.

A method for determining two substrates binding in the same active site of cytochrome P450BM3: an explanation of high energy omega product formation

A number of enzymes from the cytochrome P450 family show atypical (non-Michaelis-Menten) kinetic behavior resulting from substrate activation, inhibition, partial inhibition, biphasic saturation, or autoactivation. Herein, we provide a technique that can identify multiple substrate occupancy in the same active site of a P450 as a result of an altered kinetic profile. Using an isotope effect on product ratios confirms that the enzyme-substrate (ES) complex responsible for omega hydroxylation of palmitic acid (palmitate) is in rapid equilibrium with the ES complex that leads to omega-1 hydroxylation of palmitate. Co-incubation of a second substrate, lauric acid (laurate), results in a change in the ratio of omega to omega-1 hydroxylated palmitate. Furthermore, an isotope effect on palmitate is observed when deuterated laurate is co-incubated with non-deuterated palmitate. These results are only consistent with both substrates being in the same active site simultaneously. This mode of binding explains how the F87A mutant of P450BM3 is able to produce the omega alcohol, a product that arises from the high-energy primary radical.

A new class of CYP2C9 inhibitors: probing 2C9 specificity with high-affinity benzbromarone derivatives

Noncovalent forces, other than hydrophobic interactions, are important determinants of substrate bias exhibited by some cytochromes P450. The CYP2C9 pharmacophore is proposed to include either an anionic group or hydrogen bond donor in addition to its hydrophobic groups. By constructing analogs of benzbromarone, evidence supporting the existence of a 2C9 anion-binding site was revealed. A nonsubstituted phenol analog was determined to have a pKa of 8.4 and a Ki of 414 nM whereas those with dihalogenated benzoyl phenols had pKa values between 4.2 to 5.2 and Ki values as low as 1 nM. The nonhalogenated, nonionizable analog is the poorest binder at 796 nM. The Ki range covers around three orders of magnitude with even the weakest binder being a more potent inhibitor than 2C9 substrate phenytoin. Thus, benzbromarone derivatives represent a class of molecules with the potential to reveal more structural details of the 2C9 active site.

Evidence for two different active oxygen species in cytochrome P450 BM3 mediated sulfoxidation and N-dealkylation reactions

Herein, we report the results from two experiments that are consistent with sulfoxidation and N-dealkylation involving two different enzyme substrate complexes and thus two different active oxygen species that do not interchange. The first experiment involves the use of a mutant that may increase the amount of the hydroperoxy-iron species (FeIIIO2H).1 This mutant increases the amount of sulfoxidation relative to the amount of N-dealkylation by 4-fold. In a second experiment, deuterium substitution on the N-methyl groups of substrate does not result in an increase in sulfoxidation. This later result is consistent with N-dealkylation and sulfoxidation being mediated by two different active oxygen species. While the data indicate two active oxygen species, they do not distinguish between the different possibilities for the active oxygen species.