Simulation of interindividual differences in inactivation of reactive para-benzoquinone imine metabolites of diclofenac by glutathione S-transferases in human liver cytosol.

Diclofenac is a widely prescribed NSAID that causes severe idiosyncratic drug induced liver injury (IDILI) in a small part of the patient population. Formation of protein-reactive metabolites is considered to play a role in the development of diclofenac-induced IDILI. Therefore, a high hepatic activity of enzymes involved in bioactivation of diclofenac is expected to increase the risk for liver injury. However, the extent of covalent protein binding may also be determined by activity of protective enzymes, such as glutathione S-transferases (GSTs). This is supported by an association study in which a correlation was found between NSAID-induced IDILI and the combined null genotypes of GSTM1 and GSTT1. In the present study, the activity of 10 different recombinant human GSTs in inactivation of protein-reactive quinoneimine (QI) metabolites of diclofenac was tested. Both at low and high GSH concentrations, high activities of GSTA1-1, A2-2, A3-3, M1-1, M3-3 and P1-1 in the inactivation of these QIs were found. By using the expression levels of GSTs in livers of 22 donors, a 6-fold variation in GST-dependent inactivation of reactive diclofenac metabolites was predicted. Moreover, it was shown in vitro that GSTs can strongly increase the efficiency of GSH to protect against the alkylation of the model thiol N-acetylcysteine by reactive diclofenac metabolites. The results of this study demonstrate that variability of GST expression may significantly contribute to the inter-individual differences in susceptibility to diclofenac-induced liver injury. In addition, expression levels of GSTs in in vitro models for hepatotoxicity may be important factors determining sensitivity to diclofenac cytotoxicity.

Biosynthesis of a steroid metabolite by an engineered Rhodococcus erythropolis strain expressing a mutant cytochrome P450 BM3 enzyme.

In the present study, the use of Rhodococcus erythropolis mutant strain RG9 expressing the cytochrome P450 BM3 mutant M02 enzyme has been evaluated for whole-cell biotransformation of a 17-ketosteroid, norandrostenedione, as a model substrate. Purified P450 BM3 mutant M02 enzyme hydroxylated the steroid with >95 % regioselectivity to form 16-ß-OH norandrostenedione, as confirmed by NMR analysis. Whole cells of R. erythropolis RG9 expressing P450 BM3 M02 enzyme also converted norandrostenedione into the 16-ß-hydroxylated product, resulting in the formation of about 0.35 g/L. Whole cells of strain RG9 itself did not convert norandrostenedione, indicating that metabolite formation is P450 BM3 M02 enzyme mediated. This study shows that R. erythropolis is a novel and interesting host for the heterologous expression of highly selective and active P450 BM3 M02 enzyme variants able to perform steroid bioconversions.

Cytochrome P450-mediated bioactivation of mefenamic acid to quinoneimine intermediates and inactivation by human glutathione S-transferases.

Mefenamic acid (MFA) has been associated with rare but severe cases of hepatotoxicity, nephrotoxicity, gastrointestinal toxicity, and hypersensitivity reactions that are believed to result from the formation of reactive metabolites. Although formation of protein-reactive acylating metabolites by phase II metabolism has been well-studied and proposed to be the cause of these toxic side effects, the oxidative bioactivation of MFA has not yet been competely characterized. In the present study, the oxidative bioactivation of MFA was studied using human liver microsomes (HLM) and recombinant human P450 enzymes. In addition to the major metabolite 3'-OH-methyl-MFA, resulting from the benzylic hydroxylation by CYP2C9, 4'-hydroxy-MFA and 5-hydroxy-MFA were identified as metabolites resulting from oxidative metabolism of both aromatic rings of MFA. In the presence of GSH, three GSH conjugates were formed that appeared to result from GSH conjugation of the two quinoneimines formed by further oxidation of 4'-hydroxy-MFA and 5-hydroxy-MFA. The major GSH conjugate was identified as 4'-OH-5'-glutathionyl-MFA and was formed at the highest activity by CYP1A2 and to a lesser extent by CYP2C9 and CYP3A4. Two minor GSH conjugates resulted from secondary oxidation of 5-hydroxy-MFA and were formed at the highest activity by CYP1A2 and to a lesser extent by CYP3A4. Additionally, the ability of seven human glutathione S-transferases (hGSTs) to catalyze the GSH conjugation of the quinoneimines formed by P450s was also investigated. The highest increase of total GSH conjugation was observed with hGSTP1-1, followed by hepatic hGSTs hGSTA2-2 and hGSTM1-1. The results of this study show that, next to phase II metabolites, reactive quinoneimines formed by oxidative bioactivation might also contribute to the idiosyncratic toxicity of MFA.

Application of engineered cytochrome P450 mutants as biocatalysts for the synthesis of benzylic and aromatic metabolites of fenamic acid NSAIDs.

Cytochrome P450 BM3 mutants are promising biocatalysts for the production of drug metabolites. In the present study, the ability of cytochrome P450 BM3 mutants to produce oxidative metabolites of structurally related NSAIDs meclofenamic acid, mefenamic acid and tolfenamic acid was investigated. A library of engineered P450 BM3 mutants was screened with meclofenamic acid (1) to identify catalytically active and selective mutants. Three mono-hydroxylated metabolites were identified for 1. The hydroxylated products were confirmed by NMR analysis to be 3'-OH-methyl-meclofenamic acid (1a), 5-OH-meclofenamic acid (1b) and 4'-OH-meclofenamic acid (1c) which are human relevant metabolites. P450 BM3 variants containing V87I and V87F mutation showed high selectivity for benzylic and aromatic hydroxylation of 1 respectively. The applicability of these mutants to selectively hydroxylate structurally similar drugs such as mefenamic acid (2) and tolfenamic acid (3) was also investigated. The tested variants showed high total turnover numbers in the order of 4000-6000 and can be used as biocatalysts for preparative scale synthesis. Both 1 and 2 could undergo benzylic and aromatic hydroxylation by the P450 BM3 mutants, whereas 3 was hydroxylated only on aromatic rings. The P450 BM3 variant M11 V87F hydroxylated the aromatic ring at 4' position of all three drugs tested with high regioselectivity. Reference metabolites produced by P450 BM3 mutants allowed the characterisation of activity and regioselectivity of metabolism of all three NSAIDs by thirteen recombinant human P450s. In conclusion, engineered P450 BM3 mutants that are capable of benzylic or aromatic hydroxylation of fenamic acid containing NSAIDs, with high selectivity and turnover numbers have been identified. This shows their potential use as a greener alternative for the generation of drug metabolites.

Human NAD(P)H:quinone oxidoreductase 1 (NQO1)-mediated inactivation of reactive quinoneimine metabolites of diclofenac and mefenamic acid.

NAD(P)H: quinone oxidoreductase 1 (NQO1) is an enzyme capable of reducing a broad range of chemically reactive quinones and quinoneimines (QIs) and can be strongly upregulated by Nrf2/Keap1-mediated stress responses. Several commonly used drugs implicated in adverse drug reactions (ADRs) are known to form reactive QI metabolites upon bioactivation by P450, such as acetaminophen (APAP), diclofenac (DF), and mefenamic acid (MFA). In the present study, the reductive activity of human NQO1 toward the QI metabolites derived from APAP and hydroxy-metabolites of DF and MFA was studied, using purified bacterial P450 BM3 (CYP102A1) mutant M11 as a bioactivation system. The NQO1-catalyzed reduction of the QI metabolites was quantified relative to spontaneous glutathione (GSH) conjugation. Addition of NQO1 to the incubations strongly reduced the formation of all corresponding GSH conjugates, and this activity could be prevented by dicoumarol, a selective NQO1 inhibitor. The GSH conjugation was strongly increased by adding human GSTP1-1 in a wide range of GSH concentrations. Still, NQO1 could effectively compete with the GST catalyzed GSH conjugation by reducing the QIs. In conclusion, we identified the QI metabolites of the 4'- and 5-hydroxy-metabolites of DF and MFA as novel substrates for human NQO1. NQO1-mediated reduction proves to be an effective pathway to detoxify these QI metabolites in addition to GSH conjugation. Genetically determined deficiency of NQO1 therefore might be a risk factor for ADRs induced by reactive QI drug metabolites.

Effect of human glutathione S-transferase hGSTP1-1 polymorphism on the detoxification of reactive metabolites of clozapine, diclofenac and acetaminophen.

Recent association studies suggest that genetically determined deficiencies in GSTs might be a risk factor for idiosyncratic adverse drug reactions resulting from the formation of reactive drug metabolites. hGSTP1-1 is polymorphic in the human population with a number of single nucleotide polymorphisms that yield an amino acid change in the encoded protein. Three allelic variants of hGSTP1-1 containing an Ile105Val or Ala114Val substitution, or a combination of both, have been the most widely studied and showed different activity when compared to wild-type hGSTP1-1*A (Ile105/Ala114). In the present study, we studied the ability of these allelic variants to catalyze the GSH conjugation of reactive metabolites of acetaminophen, clozapine, and diclofenac formed by bioactivation in in vitro incubations by human liver microsomes and drug metabolizing P450 BM3 mutants. The results show that effects of the change of amino acid at residue 105 and 114 on conjugation reactions were substrate dependent. A single substitution at residue 105 affects the ability to catalyze GSH conjugation, while when both residue 105 and 114 were substituted the effect was additionally enhanced. Single mutation at position 114 did not show a significant effect. The different hGSTP1-1 mutants showed slightly altered regioselectivities in formation of individual GSH conjugates of clozapine which suggests that the binding orientation of the reactive nitrenium ion of clozapine is affected by the mutations. For diclofenac, a significant decrease in activity in GSH-conjugation of diclofenac 1',4'-quinone imine was observed for variants hGSTP1-1*B (Val105/Ala114) and hGSTP1-1*C (Val105/Val114). However, since the differences in total GSH conjugation activity catalyzed by these allelic variants were not higher than 30%, differences in inactivation of reactive intermediates by hGSTP1-1 are not likely to be a major factor in determining interindividual difference in susceptibility to adverse drug reactions induced by the drugs studied.

Free energy calculations give insight into the stereoselective hydroxylation of α-ionones by engineered cytochrome P450 BM3 mutants.

Previously, stereoselective hydroxylation of α-ionone by Cytochrome P450 BM3 mutants M01 A82W and M11 L437N was observed. While both mutants hydroxylate α-ionone in a regioselective manner at the C3 position, M01 A82W catalyzes formation of trans-3-OH-α-ionone products whereas M11 L437N exhibits opposite stereoselectivity, producing trans-(3S,6S)-OH-α-ionone and cis-(3S,6R)-OH-α-ionone. Here, we explore the stereoselective C3 hydroxylation of α-ionone by Cytochrome P450 BM3 mutants M01 A82W and M11 L437N using molecular dynamics-based free energy calculations to study the interaction between the enzyme and both the substrates and the products. The one-step perturbation approach is applied using an optimized reference state for substrates and products. While the free energy differences between the substrates free in solution amount to ~0 kJ mol(-1), the differences in mutant M01 A82W agree with the experimentally obtained dissociation constants K(d). Moreover, a correlation with experimentally observed trends in product formation is found in both mutants. The trans isomers show the most favorable relative binding free energy in the range of all four possible hydroxylated diastereomers for mutant M01 A82W, while the trans product from (6S)-α-ionone and the cis product from (6R)-α-ionone show highest affinity for mutant M11 L437N. Marcus theory is subsequently used to relate the thermodynamic stability to transition state energies and rates of formation.

A single active site mutation inverts stereoselectivity of 16-hydroxylation of testosterone catalyzed by engineered cytochrome P450 BM3.

Inversion of stereoselectivity: screening of a minimal mutant library revealed a cytochrome P450 BM3 variant M01 A82W S72I capable of producing 16 α-OH-testosterone. Remarkably, a single active site mutation S72I in M01 A82W inverted the stereoselectivity of hydroxylation from 16 ß to 16 α. Introduction of S72I mutation in another 16 ß-OH-selective variant M11 V87I, also resulted in similar inversion of stereoselectivity.

The role of protein plasticity in computational rationalization studies on regioselectivity in testosterone hydroxylation by cytochrome P450 BM3 mutants.

Recently, it was found that mutations in the binding cavity of drug-metabolizing Cytochrome P450 BM3 mutants can result in major changes in regioselectivity in testosterone (TES) hydroxylation. In the current work, we report the intrinsic reactivity of TES' C-H bonds and our attempts to rationalize experimentally observed changes in TES hydroxylation using a protein structure-based in silico approach, by setting up and employing a combined Molecular Dynamics (MD) and ligand docking approach to account for the flexibility and plasticity of BM3 mutants. Using this approach, about 100,000 TES binding poses were obtained per mutant. The predicted regioselectivity in TES hydroxylation by the mutants was found to be in disagreement with experiment. As revealed in a detailed structural analysis of the obtained docking poses, this disagreement is due to limitations in correctly scoring hydrogen-bonding and steric interactions with specific active-site residues, which could explain the experimentally observed trends in regioselectivity in TES hydroxylation.