Human multidrug resistance protein 4 (MRP4) is a cellular efflux transporter for paracetamol glutathione and cysteine conjugates.

Paracetamol (acetaminophen, APAP) overdose is a leading cause of acute drug-induced liver failure. APAP hepatotoxicity is mediated by the reactive metabolite N-acetyl-p-benzoquinone imine (NAPQI). NAPQI is inactivated by conjugation with glutathione (GSH) to APAP-GSH, which is further converted into its cysteine derivative APAP-CYS. Before necrosis of hepatocytes occurs, APAP-CYS is measurable in plasma of the affected patient and it has been proposed as an early biomarker of acetaminophen toxicity. APAP-GSH and APAP-CYS can be extruded by hepatocytes, but the transporters involved are unknown. In this study we examined whether ATP-binding cassette (ABC) transporters play a role in the cellular efflux of APAP, APAP-GSH, and APAP-CYS. The ABC transport proteins P-gp/ABCB1, BSEP/ABCB11, BCRP/ABCG2, and MRP/ABCC1-5 were overexpressed in HEK293 cells and membrane vesicles were produced. Whereas P-gp, BSEP, MRP3, MRP5, and BCRP did not transport any of the compounds, uptake of APAP-GSH was found for MRP1, MRP2 and MRP4. APAP-CYS appeared to be a substrate of MRP4 and none of the ABC proteins transported APAP. The results suggest that the NAPQI metabolite APAP-CYS can be excreted into plasma by MRP4, where it could be a useful biomarker for APAP exposure and toxicity. Characterization of the cellular efflux of APAP-CYS is important for its development as a biomarker, because plasma concentrations might be influenced by drug-transporter interactions and upregulation of MRP4.

A Comparative Linear Interaction Energy and MM/PBSA Study on SIRT1-Ligand Binding Free Energy Calculation.

Binding free energy (ΔGbind) computation can play an important role in prioritizing compounds to be evaluated experimentally on their affinity for target proteins, yet fast and accurate ΔGbind calculation remains an elusive task. In this study, we compare the performance of two popular end-point methods, i.e., linear interaction energy (LIE) and molecular mechanics/Poisson-Boltzmann surface area (MM/PBSA), with respect to their ability to correlate calculated binding affinities of 27 thieno[3,2-d]pyrimidine-6-carboxamide-derived sirtuin 1 (SIRT1) inhibitors with experimental data. Compared with the standard single-trajectory setup of MM/PBSA, our study elucidates that LIE allows to obtain direct ("absolute") values for SIRT1 binding free energies with lower compute requirements, while the accuracy in calculating relative values for ΔGbind is comparable (Pearson's r = 0.72 and 0.64 for LIE and MM/PBSA, respectively). We also investigate the potential of combining multiple docking poses in iterative LIE models and find that Boltzmann-like weighting of outcomes of simulations starting from different poses can retrieve appropriate binding orientations. In addition, we find that in this particular case study the LIE and MM/PBSA models can be optimized by neglecting the contributions from electrostatic and polar interactions to the ΔGbind calculations.

Reduction and Scavenging of Chemically Reactive Drug Metabolites by NAD(P)H:Quinone Oxidoreductase 1 and NRH:Quinone Oxidoreductase 2 and Variability in Hepatic Concentrations.

Detoxicating enzymes NAD(P)H:quinone oxidoreductase 1 (NQO1) and NRH:quinone oxidoreductase 2 (NQO2) catalyze the two-electron reduction of quinone-like compounds. The protective role of the polymorphic NQO1 and NQO2 enzymes is especially of interest in the liver as the major site of drug bioactivation to chemically reactive drug metabolites. In the current study, we quantified the concentrations of NQO1 and NQO2 in 20 human liver donors and NQO1 and NQO2 activities with quinone-like drug metabolites. Hepatic NQO1 concentrations ranged from 8 to 213 nM. Using recombinant NQO1, we showed that low nM concentrations of NQO1 are sufficient to reduce synthetic amodiaquine and carbamazepine quinone-like metabolites in vitro. Hepatic NQO2 concentrations ranged from 2 to 31 µM. NQO2 catalyzed the reduction of quinone-like metabolites derived from acetaminophen, clozapine, 4'-hydroxydiclofenac, mefenamic acid, amodiaquine, and carbamazepine. The reduction of the clozapine nitrenium ion supports association studies showing that NQO2 is a genetic risk factor for clozapine-induced agranulocytosis. The 5-hydroxydiclofenac quinone imine, which was previously shown to be reduced by NQO1, was not reduced by NQO2. Tacrine was identified as a potent NQO2 inhibitor and was applied to further confirm the catalytic activity of NQO2 in these assays. While the in vivo relevance of NQO2-catalyzed reduction of quinone-like metabolites remains to be established by identification of the physiologically relevant co-substrates, our results suggest an additional protective role of the NQO2 protein by non-enzymatic scavenging of quinone-like metabolites. Hepatic NQO1 activity in detoxication of quinone-like metabolites becomes especially important when other detoxication pathways are exhausted and NQO1 levels are induced.

Acetaminophen reduces the protein levels of high affinity amino acid permeases and causes tryptophan depletion.

In yeast, toxicity of acetaminophen (APAP), a frequently used analgesic and antipyretic drug, depends on ubiquitin-controlled processes. Previously, we showed a remarkable overlap in toxicity profiles between APAP and tyrosine, and a similarity with drugs like rapamycin and quinine, which induce degradation of the amino acid permease Tat2. Therefore, we investigated in yeast whether APAP reduced the expression levels of amino acid permeases. The protein levels of Tat2, Tat1, Mup1 and Hip1 were reduced, while the expression of the general permease Gap1 was increased, consistent with a nutrient starvation response. Overexpression of Tat1 and Tat2, but not Mup1, Hip1 and Gap1 conferred resistance to APAP. A tryptophan auxotrophic strain trp1Δ was more sensitive to APAP than wild-type and addition of tryptophan completely restored the growth restriction of trp1∆ upon APAP exposure, while tyrosine had an additive effect on APAP toxicity. Furthermore, intracellular aromatic amino acid concentrations were reduced upon APAP exposure. This effect was less prominent in ubiquitin-deficient yeast strains that were APAP resistant and showed a reduced degradation of high affinity amino acid permeases. APAP-induced changes in intracellular amino acid concentrations were also detected in hepatoma HepG2 cells indicating significance for humans.

Binding free energy predictions of farnesoid X receptor (FXR) agonists using a linear interaction energy (LIE) approach with reliability estimation: application to the D3R Grand Challenge 2.

Computational protein binding affinity prediction can play an important role in drug research but performing efficient and accurate binding free energy calculations is still challenging. In the context of phase 2 of the Drug Design Data Resource (D3R) Grand Challenge 2 we used our automated eTOX ALLIES approach to apply the (iterative) linear interaction energy (LIE) method and we evaluated its performance in predicting binding affinities for farnesoid X receptor (FXR) agonists. Efficiency was obtained by our pre-calibrated LIE models and molecular dynamics (MD) simulations at the nanosecond scale, while predictive accuracy was obtained for a small subset of compounds. Using our recently introduced reliability estimation metrics, we could classify predictions with higher confidence by featuring an applicability domain (AD) analysis in combination with protein-ligand interaction profiling. The outcomes of and agreement between our AD and interaction-profile analyses to distinguish and rationalize the performance of our predictions highlighted the relevance of sufficiently exploring protein-ligand interactions during training and it demonstrated the possibility to quantitatively and efficiently evaluate if this is achieved by using simulation data only.

Linking cytochrome P450 enzymes from Mycobacterium tuberculosis to their cognate ferredoxin partners.

Mycobacterium tuberculosis (Mtb) codes for 20 cytochrome P450 enzymes (CYPs), considered potential drug-targets due to their essential roles in bacterial viability and host infection. Catalytic activity of mycobacterial CYPs is dependent on electron transfer from a NAD (P)H-ferredoxin-reductase (FNR) and a ferredoxin (Fd). Two FNRs (FdrA and FprA) and five ferredoxins (Fdx, FdxA, FdxC, FdxD, and Rv1786) have been found in the Mtb genome. However, as of yet, the cognate redox partnerships have not been fully established. This is confounded by the fact that heterologous redox partners are routinely used to reconstitute Mtb CYP metabolism. To this end, this study aimed to biochemically characterize and identify cognate redox partnerships for Mtb CYPs. Interestingly, all combinations of FNRs and ferredoxins were active in the reduction of oxidized cytochrome c, but steady-state kinetic assays revealed FdxD as the most efficient redox partner for FdrA, whereas Fdx coupled preferably with FprA. CYP121A1, CYP124A1, CYP125A1, and CYP142A1 metabolism with the cognate redox partners was reconstituted in vitro showing an unanticipated selectivity in the requirement for electron transfer partnership, which did not necessarily correlate with proximity in the genome. This is the first description of microbial P450 metabolism in which multiple ferredoxins are functionally linked to multiple CYPs.

Characterization of human cytochrome P450 mediated bioactivation of amodiaquine and its major metabolite N-desethylamodiaquine.

AIMS: Oxidative bioactivation of amodiaquine (AQ) by cytochrome P450s to a reactive quinoneimine is considered as an important mechanism underlying its idiosyncratic hepatotoxicity. However, because internal exposure to its major metabolite N-desethylamodiaquine (DEAQ) is up to 240-fold higher than AQ, bioactivation of DEAQ might significantly contribute to covalent binding. The aim of the present study was to compare the kinetics of bioactivation of AQ and DEAQ by human liver microsomes (HLM) and to characterize the CYPs involved in bioactivation of AQ and DEAQ. METHODS: Glutathione was used to trap reactive metabolites formed in incubations of AQ and DEAQ with HLM and recombinant human cytochrome P450s (hCYPs). Kinetics of bioactivation of AQ and DEAQ in HLM and involvement of hCYPs were characterized by measuring corresponding glutathione conjugates (AQ-SG and DEAQ-SG) using a high-performance liquid chromatography method. RESULTS: Bioactivation of AQ and DEAQ in HLM both exhibited Michaelis-Menten kinetics. For AQ bioactivation, enzyme kinetical parameters were Km , 11.5 ± 2.0 µmol l-1 , Vmax , 59.2 ± 3.2 pmol min-1  mg-1 and CLint , 5.15 µl min-1  mg-1 . For DEAQ, parameters for bioactivation were Km , 6.1 ± 1.3 µmol l-1 , Vmax , 5.5 ± 0.4 pmol min-1  mg-1 and CLint 0.90 µl min-1  mg-1 . Recombinant hCYPs and inhibition studies with HLM showed involvement of CYP3A4, CYP2C8, CYP2C9 and CYP2D6 in bioactivation. CONCLUSIONS: The major metabolite DEAQ is likely to be quantitatively more important than AQ with respect to hepatic exposure to reactive metabolites in vivo. High expression of CYP3A4, CYP2C8, CYP2C9, and CYP2D6 may be risk factors for hepatotoxicity caused by AQ-therapy.

Comprehensive and Automated Linear Interaction Energy Based Binding-Affinity Prediction for Multifarious Cytochrome P450 Aromatase Inhibitors.

Cytochrome P450 aromatase (CYP19A1) plays a key role in the development of estrogen dependent breast cancer, and aromatase inhibitors have been at the front line of treatment for the past three decades. The development of potent, selective and safer inhibitors is ongoing with in silico screening methods playing a more prominent role in the search for promising lead compounds in bioactivity-relevant chemical space. Here we present a set of comprehensive binding affinity prediction models for CYP19A1 using our automated Linear Interaction Energy (LIE) based workflow on a set of 132 putative and structurally diverse aromatase inhibitors obtained from a typical industrial screening study. We extended the workflow with machine learning methods to automatically cluster training and test compounds in order to maximize the number of explained compounds in one or more predictive LIE models. The method uses protein-ligand interaction profiles obtained from Molecular Dynamics (MD) trajectories to help model search and define the applicability domain of the resolved models. Our method was successful in accounting for 86% of the data set in 3 robust models that show high correlation between calculated and observed values for ligand-binding free energies (RMSE < 2.5 kJ mol-1), with good cross-validation statistics.

Insights into regioselective metabolism of mefenamic acid by cytochrome P450 BM3 mutants through crystallography, docking, molecular dynamics, and free energy calculations.

Cytochrome P450 BM3 (CYP102A1) mutant M11 is able to metabolize a wide range of drugs and drug-like compounds. Among these, M11 was recently found to be able to catalyze formation of human metabolites of mefenamic acid and other nonsteroidal anti-inflammatory drugs (NSAIDs). Interestingly, single active-site mutations such as V87I were reported to invert regioselectivity in NSAID hydroxylation. In this work, we combine crystallography and molecular simulation to study the effect of single mutations on binding and regioselective metabolism of mefenamic acid by M11 mutants. The heme domain of the protein mutant M11 was expressed, purified, and crystallized, and its X-ray structure was used as template for modeling. A multistep approach was used that combines molecular docking, molecular dynamics (MD) simulation, and binding free-energy calculations to address protein flexibility. In this way, preferred binding modes that are consistent with oxidation at the experimentally observed sites of metabolism (SOMs) were identified. Whereas docking could not be used to retrospectively predict experimental trends in regioselectivity, we were able to rank binding modes in line with the preferred SOMs of mefenamic acid by M11 and its mutants by including protein flexibility and dynamics in free-energy computation. In addition, we could obtain structural insights into the change in regioselectivity of mefenamic acid hydroxylation due to single active-site mutations. Our findings confirm that use of MD and binding free-energy calculation is useful for studying biocatalysis in those cases in which enzyme binding is a critical event in determining the selective metabolism of a substrate.

Different Reactive Metabolites of Nevirapine Require Distinct Glutathione S-Transferase Isoforms for Bioinactivation.

Nevirapine (NVP) is a non-nucleoside reverse transcriptase-inhibitor, which is associated with severe idiosyncratic skin rash and hepatotoxicity. These adverse drug reactions are believed to be mediated by the formation of epoxides and/or quinone methide formed by oxidative metabolism by P450s and 12-sulfoxyl-NVP formed by sequential 12-hydroxylation and O-sulfonation. Although different GSH-conjugates and corresponding mercapturic acids have been demonstrated previously in vitro and in vivo, the role of the glutathione S-transferases in the inactivation of the different reactive metabolites has not been studied so far. In the present study the activity of 10 recombinant human glutathione S-transferases (GSTs) in the detoxification of the different reactive metabolites of NVP was studied. The results show that GSTP1-1 is a highly active catalyst of GSH-conjugation of the oxidative metabolites of NVP, even at high GSH-concentration. Experiments with trideuterated NVP suggest involvement of a reactive epoxide rather than quinone methide in the formation of the GSH-conjugate formed after oxidative bioactivation. GSH-conjugation of 12-sulfoxyl-NVP forming NVP-12-GSH was only catalyzed by GSTM1-1, GSTA1-1, and GSTA3-3. Although the exact expression levels of these enzymes in the skin is unknown, the relatively low activity of this catalysis makes it unlikely that GSTs can provide significant protection against this metabolite. However, since NVP-12-GSH is specifically formed via the 12-sulfoxyl-NVP, its corresponding urinary mercapturic acid can be considered as a biomarker for recent internal exposure to this protein-reactive sulfate. However, it has to be taken into account that 12-sulfoxyl-NVP is not completely trapped by GSH and that rates of bioinactivation will differ between patients due to variability in expression of GSTM1, GSTA1, and GSTA3.

Application of a cocktail approach to screen cytochrome P450 BM3 libraries for metabolic activity and diversity.

In the present study, the validity of using a cocktail screening method in combination with a chemometrical data mining approach to evaluate metabolic activity and diversity of drug-metabolizing bacterial Cytochrome P450 (CYP) BM3 mutants was investigated. In addition, the concept of utilizing an in-house-developed library of CYP BM3 mutants as a unique biocatalytic synthetic tool to support medicinal chemistry was evaluated. Metabolic efficiency of the mutant library towards a selection of CYP model substrates, being amitriptyline (AMI), buspirone (BUS), coumarine (COU), dextromethorphan (DEX), diclofenac (DIC) and norethisterone (NET), was investigated. First, metabolic activity of a selection of CYP BM3 mutants was screened against AMI and BUS. Subsequently, for a single CYP BM3 mutant, the effect of co-administration of multiple drugs on the metabolic activity and diversity towards AMI and BUS was investigated. Finally, a cocktail of AMI, BUS, COU, DEX, DIC and NET was screened against the whole in-house CYP BM3 library. Different validated quantitative and qualitative (U)HPLC-MS/MS-based analytical methods were applied to screen for substrate depletion and targeted product formation, followed by a more in-depth screen for metabolic diversity. A chemometrical approach was used to mine all data to search for unique metabolic properties of the mutants and allow classification of the mutants. The latter would open the possibility of obtaining a more in-depth mechanistic understanding of the metabolites. The presented method is the first MS-based method to screen CYP BM3 mutant libraries for diversity in combination with a chemometrical approach to interpret results and visualize differences between the tested mutants.

Evidence-based selection of training compounds for use in the mechanism-based integrated prediction of drug-induced liver injury in man.

The current test systems employed by pharmaceutical industry are poorly predictive for drug-induced liver injury (DILI). The 'MIP-DILI' project addresses this situation by the development of innovative preclinical test systems which are both mechanism-based and of physiological, pharmacological and pathological relevance to DILI in humans. An iterative, tiered approach with respect to test compounds, test systems, bioanalysis and systems analysis is adopted to evaluate existing models and develop new models that can provide validated test systems with respect to the prediction of specific forms of DILI and further elucidation of mechanisms. An essential component of this effort is the choice of compound training set that will be used to inform refinement and/or development of new model systems that allow prediction based on knowledge of mechanisms, in a tiered fashion. In this review, we focus on the selection of MIP-DILI training compounds for mechanism-based evaluation of non-clinical prediction of DILI. The selected compounds address both hepatocellular and cholestatic DILI patterns in man, covering a broad range of pharmacologies and chemistries, and taking into account available data on potential DILI mechanisms (e.g. mitochondrial injury, reactive metabolites, biliary transport inhibition, and immune responses). Known mechanisms by which these compounds are believed to cause liver injury have been described, where many if not all drugs in this review appear to exhibit multiple toxicological mechanisms. Thus, the training compounds selection offered a valuable tool to profile DILI mechanisms and to interrogate existing and novel in vitro systems for the prediction of human DILI.

Characterization of human cytochrome P450s involved in the bioactivation of tri-ortho-cresyl phosphate (ToCP).

Tri-ortho-cresyl phosphate (ToCP) is a multipurpose organophosphorus compound that is neurotoxic and suspected to be involved in aerotoxic syndrome in humans. It has been reported that not ToCP itself but a metabolite of ToCP, namely, 2-(ortho-cresyl)-4H-1,2,3-benzodioxaphosphoran-2-one (CBDP), may be responsible for this effect as it can irreversibly bind to human butyrylcholinesterase (BuChE) and human acetylcholinesterase (AChE). The bioactivation of ToCP into CBDP involves Cytochrome P450s (P450s). However, the individual human P450s responsible for this bioactivation have not been identified yet. In the present study, we aimed to investigate the metabolism of ToCP by different P450s and to determine the inhibitory effect of the in vitro generated ToCP-metabolites on human BuChE and AChE. Human liver microsomes, rat liver microsomes, and recombinant human P450s were used for that purpose. The recombinant P450s 2B6, 2C18, 2D6, 3A4 and 3A5 showed highest activity of ToCP-bioactivation to BuChE-inhibitory metabolites. Inhibition experiments using pooled human liver microsomes indicated that P450 3A4 and 3A5 were mainly involved in human hepatic bioactivation of ToCP. In addition, these experiments indicated a minor role for P450 1A2. Formation of CBDP by in-house expressed recombinant human P450s 1A2 and 3A4 was proven by both LC-MS and GC-MS analysis. When ToCP was incubated with P450 1A2 and 3A4 in the presence of human BuChE, CBDP-BuChE-adducts were detected by LC-MS/MS which were not present in the corresponding control incubations. These results confirmed the role of human P450s 1A2 and 3A4 in ToCP metabolism and demonstrated that CBDP is the metabolite responsible for the BuChE inactivation. Interindividual differences at the level of P450 1A2 and 3A4 might play an important role in the susceptibility of humans in developing neurotoxic effects, such as aerotoxic syndrome, after exposure to ToCP.

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.

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.

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.

Towards automated binding affinity prediction using an iterative linear interaction energy approach.

Binding affinity prediction of potential drugs to target and off-target proteins is an essential asset in drug development. These predictions require the calculation of binding free energies. In such calculations, it is a major challenge to properly account for both the dynamic nature of the protein and the possible variety of ligand-binding orientations, while keeping computational costs tractable. Recently, an iterative Linear Interaction Energy (LIE) approach was introduced, in which results from multiple simulations of a protein-ligand complex are combined into a single binding free energy using a Boltzmann weighting-based scheme. This method was shown to reach experimental accuracy for flexible proteins while retaining the computational efficiency of the general LIE approach. Here, we show that the iterative LIE approach can be used to predict binding affinities in an automated way. A workflow was designed using preselected protein conformations, automated ligand docking and clustering, and a (semi-)automated molecular dynamics simulation setup. We show that using this workflow, binding affinities of aryloxypropanolamines to the malleable Cytochrome P450 2D6 enzyme can be predicted without a priori knowledge of dominant protein-ligand conformations. In addition, we provide an outlook for an approach to assess the quality of the LIE predictions, based on simulation outcomes only.

Characterization of human cytochrome P450s involved in the bioactivation of clozapine.

Clozapine is known to cause hepatotoxicity in a small percentage of patients. Oxidative bioactivation to reactive intermediates by hepatic cytochrome P450s (P450s) has be proposed as a possible mechanism. However, in contrast to their role in formation of N-desmethylclozapine and clozapine N-oxide, the involvement of individual P450s in the bioactivation to reactive intermediates is much less well characterized. The results of the present study show that 7 of 14 recombinant human P450s were able to bioactivate clozapine to a glutathione-reactive nitrenium ion. CYP3A4 and CYP2D6 showed the highest specific activity. Enzyme kinetical characterization of these P450s showed comparable intrinsic clearance of bioactivation, implicating that CYP3A4 would be more important because of its higher hepatic expression, compared with CYP2D6. Inhibition experiments using pooled human liver microsomes confirmed the major role of CYP3A4 in hepatic bioactivation of clozapine. By studying bioactivation of clozapine in human liver microsomes from 100 different individuals, an 8-fold variability in bioactivation activity was observed. In two individuals bioactivation activity exceeded N-demethylation and N-oxidation activity. Quinidine did not show significant inhibition of bioactivation in any of these liver fractions, suggesting that CYP2D6 polymorphism is not an important factor in determining susceptibility to hepatotoxicity of clozapine. Therefore, interindividual differences and drug-drug interactions at the level of CYP3A4 might be factors determining exposure of hepatic tissue to reactive clozapine metabolites.

Effect of human glutathione S-transferases on glutathione-dependent inactivation of cytochrome P450-dependent reactive intermediates of diclofenac.

Idiosyncratic adverse drug reactions due to the anti-inflammatory drug diclofenac have been proposed to be caused by the generation of reactive acyl glucuronides and oxidative metabolites. For the oxidative metabolism of diclofenac by cytochromes P450 at least five different reactive intermediates have been proposed previously based on structural identification of their corresponding GSH-conjugates. In the present study, the ability of four human glutathione S-transferases (hGSTs) to catalyze the GSH-conjugation of the different reactive intermediates formed by P450s was investigated. Addition of pooled human liver cytosol and recombinant hGSTA1-1, hGSTM1-1, and hGSTP1-1 to incubations of diclofenac with human liver microsomes or purified CYP102A1M11 L437N as a model system significantly increased total GSH-conjugation. The strongest increase of total GSH-conjugation was observed by adding hGSTP1-1, whereas hGSTM1-1 and hGSTA1-1 showed lower activity. Addition of hGSTT1-1 only showed a minor effect. When considering the effects of hGSTs on GSH-conjugation of the different quinoneimines of diclofenac, it was found that hGSTP1-1 showed the highest activity in GSH-conjugation of the quinoneimine derived from 5-hydroxydiclofenac (5-OH-DF). hGSTM1-1 showed the highest activity in inactivation of the quinoneimine derived from 4'-hydroxydiclofenac (4'-OH-DF). Separate incubations with 5-OH-DF and 4'-OH-DF as substrates confirmed these results. hGSTs also catalyzed GSH-conjugation of the o-iminemethide formed by oxidative decarboxylation of diclofenac as well as the substitution of one of the chlorine atoms of DF by GSH. hGSTP1-1 showed the highest activity for the formation of these minor GSH-conjugates. These results suggest that hGSTs may play an important role in the inactivation of DF quinoneimines and its minor reactive intermediates especially in stress conditions when tissue levels of GSH are decreased.

Urinary lipid and protein oxidation products upon halothane, isoflurane, or sevoflurane anesthesia in humans: potential biomarkers for a subclinical nephrotoxicity.

OBJECTIVE: To investigate whether lipid and protein oxidation products are elevated and correlated with routine clinical markers of hepatic and renal function in patients anesthetized with halothane, isoflurane, or sevoflurane. METHODS: Urine and blood samples were collected from patient groups. Excretion of aldehydes, acetone, and o,o'-dityrosine was measured before and after anesthesia. Blood samples were analysed for clinical markers. RESULTS: Urinary concentrations of aldehydes, acetone, o,o'-dityrosine and glucose were significantly increased after anesthesia in halothane and sevoflurane groups earlier than clinical markers. Significant correlations were found in sevoflurane group. CONCLUSION: Lipid and protein oxidation contributes to subclinical sevoflurane nephrotoxicity. Oxidation products may serve as early biomarkers.