Improved cytochrome P450 3A4 molecular models accurately predict the Phe215 requirement for raloxifene dehydrogenation selectivity

Formation of potentially toxic metabolites of drugs in reactions catalyzed by human drug-metabolizing enzymes

Data are presented on the formation of potentially toxic metabolites of drugs that are substrates of human drug metabolizing enzymes. The tabular data lists the formation of potentially toxic/reactive products. The data were obtained from in vitro experiments and showed that the oxidative reactions predominate (with 96% of the total potential toxication reactions). Reductive reactions (e.g., reduction of nitro to amino group and reductive dehalogenation) participate to the extent of 4%. Of the enzymes, cytochrome P450 (P450, CYP) enzymes catalyzed 72% of the reactions, myeloperoxidase (MPO) 7%, flavin-containing monooxygenase (FMO) 3%, aldehyde oxidase (AOX) 4%, sulfotransferase (SULT) 5%, and a group of minor participating enzymes to the extent of 9%. Within the P450 Superfamily, P450 Subfamily 3A (P450 3A4 and 3A5) participates to the extent of 27% and the Subfamily 2C (P450 2C9 and P450 2C19) to the extent of 16%, together catalyzing 43% of the reactions, followed by P450 Subfamily 1A (P450 1A1 and P450 1A2) with 15%. The P450 2D6 enzyme participated in an extent of 8%, P450 2E1 in 10%, and P450 2B6 in 6% of the reactions. All other enzymes participate to the extent of 14%. The data show that, of the human enzymes analyzed, P450 enzymes were dominant in catalyzing potential toxication reactions of drugs and their metabolites, with the major role assigned to the P450 Subfamily 3A and significant participation of the P450 Subfamilies 2C and 1A, plus the 2D6, 2E1 and 2B6 enzymes contributing. Selected examples of drugs that are activated or proposed to form toxic species are discussed.

Unraveling the Abnormal Molecular Mechanism of Suicide Inhibition of Cytochrome P450 3A4

Suicide inhibition of the CYP3A4 enzyme by a drug inactivates the enzyme in the drug biotransformation process and often shows safety concerns about the drug. Despite extensive experimental studies, the abnormal molecular mechanism of a suicide inhibitor that forms a covalent bond with the residue far away from the catalytically active center of CYP3A4 inactivating the enzyme remains elusive. Here, the authors used molecular simulation approaches to study in detail how diquinone methide (DQR), the metabolite product of raloxifene, unbinds from CYP3A4 and inactivates the enzyme at the atomistic level. The results clearly indicate that in one of the intermediate states formed in its unbinding process, DQR covalently binds to Cys239, a residue far away from the catalytically active center of CYP3A4, and hinders the substrate from entering or leaving the enzyme. This work therefore provides an unprecedented way of clarifying the abnormal mechanism of suicide inhibition of the CYP3A4 enzyme.

Add-On Selective Estrogen Receptor Modulators for Methadone Maintenance Treatment

Methadone maintenance treatment (MMT) remains the cornerstone for the management of opiate abuse. However, MMT can be associated with complex factors, including complications during the tolerance phase, the inability of some patients to maintain treatment effects during the tapering or abstinence phases, and the development of methadone dependence. Previous studies have revealed a sex disparity in MMT efficacy, showing that women undergoing MMT experiencing an increase in psychological symptoms compared with men and suggesting a link between disparate responses and the effects of estrogen signaling on methadone metabolism. More specifically, estradiol levels are positively associated with MMT dosing, and the expression of a single-nucleotide polymorphism (SNP) associated with estrogen receptor (ER) regulation is also associated with MMT dosing. In addition to performing mechanistic dissections of estrogen signaling in the presence of methadone, past studies have also proposed the targeting of estrogen signaling during MMT. The present report provides an overview of the relevant literature regarding sex effects, including differences in sex hormones and their potential impacts on MMT regimens. Moreover, this article provides a pharmacological perspective on the targeting of estrogen signals through the use of selective ER modulators (SERMs) during MMT. Preliminary preclinical experiments were also performed to evaluate the potential effects of targeting estrogen signaling with tamoxifen on methadone metabolism.

Human Family 1-4 cytochrome P450 enzymes involved in the metabolic activation of xenobiotic and physiological chemicals: an update

This is an overview of the metabolic activation of drugs, natural products, physiological compounds, and general chemicals by the catalytic activity of cytochrome P450 enzymes belonging to Families 1-4. The data were collected from > 5152 references. The total number of data entries of reactions catalyzed by P450s Families 1-4 was 7696 of which 1121 (~ 15%) were defined as bioactivation reactions of different degrees. The data were divided into groups of General Chemicals, Drugs, Natural Products, and Physiological Compounds, presented in tabular form. The metabolism and bioactivation of selected examples of each group are discussed. In most of the cases, the metabolites are directly toxic chemicals reacting with cell macromolecules, but in some cases the metabolites formed are not direct toxicants but participate as substrates in succeeding metabolic reactions (e.g., conjugation reactions), the products of which are final toxicants. We identified a high level of activation for three groups of compounds (General Chemicals, Drugs, and Natural Products) yielding activated metabolites and the generally low participation of Physiological Compounds in bioactivation reactions. In the group of General Chemicals, P450 enzymes 1A1, 1A2, and 1B1 dominate in the formation of activated metabolites. Drugs are mostly activated by the enzyme P450 3A4, and Natural Products by P450s 1A2, 2E1, and 3A4. Physiological Compounds showed no clearly dominant enzyme, but the highest numbers of activations are attributed to P450 1A, 1B1, and 3A enzymes. The results thus show, perhaps not surprisingly, that Physiological Compounds are infrequent substrates in bioactivation reactions catalyzed by P450 enzyme Families 1-4, with the exception of estrogens and arachidonic acid. The results thus provide information on the enzymes that activate specific groups of chemicals to toxic metabolites.

Molecular Dynamics Simulation Framework to Probe the Binding Hypothesis of CYP3A4 Inhibitors

The Cytochrome P450 family of heme-containing proteins plays a major role in catalyzing phase I metabolic reactions, and the CYP3A4 subtype is responsible for the metabolism of many currently marketed drugs. Additionally, CYP3A4 has an inherent affinity for a broad spectrum of structurally diverse chemical entities, often leading to drug-drug interactions mediated by the inhibition or induction of the metabolic enzyme. The current study explores the binding of selected highly efficient CYP3A4 inhibitors by docking and molecular dynamics (MD) simulation protocols and their binding free energy calculated using the WaterSwap method. The results indicate the importance of binding pocket residues including Phe57, Arg105, Arg106, Ser119, Arg212, Phe213, Thr309, Ser312, Ala370, Arg372, Glu374, Gly481 and Leu483 for interaction with CYP3A4 inhibitors. The residue-wise decomposition of the binding free energy from the WaterSwap method revealed the importance of binding site residues Arg106 and Arg372 in the stabilization of all the selected CYP3A4-inhibitor complexes. The WaterSwap binding energies were further complemented with the MM(GB/PB)SA results and it was observed that the binding energies calculated by both methods do not differ significantly. Overall, our results could guide towards the use of multiple computational approaches to achieve a better understanding of CYP3A4 inhibition, subsequently leading to the design of highly specific and efficient new chemical entities with suitable ADMETox properties and reduced side effects.

In silico and in vitro analysis of genetic variants of the equine CYP3A94, CYP3A95 and CYP3A97 isoenzymes

Cytochrome P450 enzymes (CYPs) of the equine CYP3A subfamily are predominantly involved in drug metabolism. In this study, genetic variants of the equine CYP3A94, CYP3A95, and CYP3A97 were identified and characterized using in silico modeling and in vitro enzyme kinetics. The genomes of 81 horses were sequenced to obtain the genetic variants. Structural CYP modifications of the most frequent variants were analyzed in silico using the 3D-structures predicted by homology modeling. Enzyme kinetic analyses were performed using testosterone as substrate. Twenty genetic variants were found including five missense variants (CYP3A94:p.Asp217Asn, CYP3A95:p.Asp214His, CYP3A95:p.Ser392Thr, CYP3A97:p.Ile119Thr, CYP3A97:p.Met500Val) with a higher percentage of minor allele frequency (MAF) (range 0.2-0.4). A splice-site variant (c.798 + 1G > A) in CYP3A94, likely to generate a truncated protein, was found in 50% of the horses. CYP3A94:p.Asp217Asn and CYP3A95:p.Asp214His were localized on the CYP F-α-helix, an important region for the substrate interactions in the human CYP3A4. Testosterone 2β-hydroxylation was diminished in CYP3A94²¹⁷Asn and CYP3A95³⁹²Thr. Ketoconazole inhibited 2β-hydroxylation differently in the five variants with the most pronounced inhibition obtained for CYP3A95³⁹²Thr. In vitro and in silico analyses of genetic variants allow unraveling structural features in equine CYPs that correlate with changes in the CYP activity.

Four Major Channels Detected in the Cytochrome P450 3A4: A Step toward Understanding Its Multispecificity

We computed the network of channels of the 3A4 isoform of the cytochrome P450 (CYP) on the basis of 16 crystal structures extracted from the Protein Data Bank (PDB). The calculations were performed with version 2 of the CCCPP software that we developed for this research project. We identified the minimal cost paths (MCPs) output by CCCPP as probable ways to access to the buried active site. The algorithm of calculation of the MCPs is presented in this paper, with its original method of visualization of the channels. We found that these MCPs constitute four major channels in CYP3A4. Among the many channels proposed by Cojocaru et al. in 2007, we found that only four of them open in 3A4. We provide a refined description of these channels together with associated quantitative data.

Structure-function analysis of porcine cytochrome P450 3A29 in the hydroxylation of T-2 toxin as revealed by docking and mutagenesis studies

T-2 toxin, one of the type A trichothecenes, presents a potential hazard to human and animal health. Our previous work demonstrated that porcine cytochrome P450 3A29 (CYP3A29) played an important role in the hydroxylation of T-2 toxin. To identify amino acids involved in this metabolic process, T-2 toxin was docked into a homology model of CYP3A29 based on a crystal structure of CYP3A4 using AutoDock 4.0. Nine residues of CYP3A29, Arg105, Arg106, Phe108, Ser119, Lys212, Phe213, Phe215, Arg372 and Glu374, which were found within 5 Å around T-2 toxin were subjected to site-directed mutagenesis. In the oxidation of nifedipine, the CL_int value of R106A was increased by nearly two-folds compared with the wild-type CYP3A29, while the substrate affinities and CL_int values of S119A and K212A were significantly reduced. In the hydroxylation of T-2 toxin, the generation of 3′-OH-T-2 by R105A, S119A and K212A was significantly less than that by the wild-type, whereas R106A slightly increased the generation of 3′-OH-T-2. These results were further confirmed by isothermal titration calorimetry analysis, suggesting that these four residues are important in the hydroxylation of T-2 toxin and Arg105 may be a specific recognition site for the toxin. Our study suggests a possible structure-function relationship of CYP3A29 in the hydroxylation of T-2 toxin, providing with new insights into the mechanism of CYP3A enzymes in the biotransformation of T-2 toxin.

Structural and energetic analysis to provide insight residues of CYP2C9, 2C11 and 2E1 involved in valproic acid dehydrogenation selectivity

Docking and molecular dynamics (MD) simulation have been two computational techniques used to gain insight about the substrate orientation within protein active sites, allowing to identify potential residues involved in the binding and catalytic mechanisms. In this study, both methods were combined to predict the regioselectivity in the binding mode of valproic acid (VPA) on three cytochrome P-450 (CYP) isoforms CYP2C9, CYP2C11, and CYP2E1, which are involved in the biotransformation of VPA yielding reactive hepatotoxic intermediate 2-n-propyl-4-pentenoic acid (4nVPA). There are experimental data about hydrogen atom abstraction of the C4-position of VPA to yield 4nVPA, however, there are not structural evidence about the binding mode of VPA and 4nVPA on CYPs. Therefore, the complexes between these CYP isoforms and VPA or 4nVPA were studied to explore their differences in binding and energetic stabilization. Docking results showed that VPA and 4nVPA are coupled into CYPs binding site in a similar conformation, but it does not explain the VPA hydrogen atom abstraction. On the other hand, MD simulations showed a set of energetic states that reorient VPA at the first ns, then making it susceptible to a dehydrogenation reaction. For 4nVPA, multiple binding modes were observed in which the different states could favor either undergo other reaction mechanism or ligand expulsion from the binding site. Otherwise, the energetic and entropic contribution point out a similar behavior for the three CYP complexes, showing as expected a more energetically favorable binding free energy for the complexes between CYPs and VPA than with 4nVPA.

Computational Survey of FHIT, A Putative Human Tumor Suppressor, Truncates Structure

BACKGROUND

Fragile Histidine Triad protein (FHIT), as a known tumor suppressor protein, has been proposed to play crucial role in inhibiting p53 degradation by MDM2. Studies have confirmed FHIT interaction with p53 or MDM2, although functional interacting domains of FHIT with MDM2 and/or p53 are not completely defined. Thus, through determining the significant structural interacting domains of FHIT, information with regard to MDM2 and p53 would be provided. As there were no previous studies evaluating the interaction of optimized important parts of target molecules, docking study was employed.

METHODS

Truncated structures of FHIT were screened to reveal critical sections engaging in FHIT interaction. HEX program was used in order to study the interaction of target structures.

RESULTS

Given the total energy, FHIT structures (β5-7, α1) and (α1) of FHIT were showed to be better candidates in comparison with other structures in interaction with optimized MDM2 part. Furthermore, FHIT structures (β4-7, α1) and (β5-7, α1) were considered to be better than other structures in interaction with optimized p53 part. FHIT truncates which interact with MDM2 optimized part exhibited lower energy levels than FHIT truncates which interact with p53 optimized part.

CONCLUSION

Our results can be useful for designing new inhibitors of this protein complex interaction which would result in tumor repression.

Protein engineering of nirobenzene dioxygenase for enantioselective synthesis of chiral sulfoxides

Nitrobenzene dioxygenase (NBDO) from Comamonas sp. is shown here to perform enantioselective oxidation of aromatic sulfides. Several para-substituted alkyl aryl sulfides were examined and it was found that the activity of the enzyme is dependent on the size of the substrate. Saturation mutagenesis was performed on different residues in the active site in order to improve activity and selectivity. Mutagenesis at position 258 in the α-hydroxylase subunit of NBDO improved both activity and enantioselectivity. Substitutions in position 293 improved the activity on all substrates and had diverse influence on enantioselectivity. Mutagenesis in position 207 provided two interesting variants, V207I and V207A, with opposite enantioselectivities. Furthermore, combining two favorable mutations, N258A and F293H, provided an improved variant with both higher activity (5.20 ± 0.01, 2.12 ± 0.21, 2.64 ± 0.14 and 4.01 ± 0.34 nmol min(-1) mg protein(-1) on thioanisole, ptolyl, Cl-thioanisole and Br-thioanisole, respectively, which is 1.7, 4.6, 7.1 and 26.7-fold compared with wild type) and improved enantioselectivity (e.g. 67% enantiomeric excess for Cl-thioanisole vs. 5% for wild type). Molecular docking and active site volume calculations were used to correlate between the structure of the substrates and the function of the enzymes. The results from this work suggest that the location of pro-chiral sulfides in the active site is coordinated by hydrophobic interactions and by steric considerations, which in turn influences the activity and enantioselectivity of NBDO.

Understanding the mechanism of cytochrome P450 3A4: recent advances and remaining problems

Cytochromes P450 (CYPs) represent a diverse group of heme-thiolate proteins found in almost all organisms. CYPs share a common protein fold but differ in substrate selectivity and catalyze a wide variety of monooxygenation reactions via activation of molecular oxygen. Among 57 human P450s, the 3A4 isoform (CYP3A4) is the most abundant and the most important because it metabolizes the majority of administered drugs. A remarkable feature of CYP3A4 is its extreme promiscuity in substrate specificity and cooperative substrate binding, which often leads to undesirable drug-drug interactions and toxic side effects. Owing to its importance in drug development and therapy, CYP3A4 has been the most extensively studied mammalian P450. In this review we provide an overview on recent progress and remaining problems in the CYP3A4 research.

Conformational dynamics of CYP3A4 demonstrate the important role of Arg212 coupled with the opening of ingress, egress and solvent channels to dehydrogenation of 4-hydroxy-tamoxifen

BACKGROUND

Structure-based methods for P450 substrates are commonly used during drug development to identify sites of metabolism. However, docking studies using available X-ray structures for the major drug-metabolizing P450, CYP3A4, do not always identify binding modes supportive of the production of high-energy toxic metabolites. Minor pathways such as P450-catalyzed dehydrogenation have been experimentally shown to produce reactive products capable of forming biomolecular adducts which can lead to increased risk toxicities. 4-Hydroxy-tamoxifen (4OHT) is metabolized by CYP3A4 via competing hydroxylation and dehydrogenation reactions.

METHODS

Ab initio gas-phase electronic structural characterization of 4OHT was used to develop a docking scoring scheme. Conformational sampling of CYP3A4 with molecular dynamics simulations along multiple trajectories were used to generate representative structures for docking studies using recently published heme parameters. A key predicted binding mode was tested experimentally using site-directed mutagenesis of CYP3A4 and liquid chromatography-mass spectroscopy analysis.

RESULTS

Docking with MD-refined CYP3A4 structures incorporating hexa-coordinate heme parameters identifies a unique binding mode involving ARG212 and channel 4, unobserved in the starting PDB ID: 1TQN X-ray structure. The models supporting dehydrogenation are consistent with results from in vitro incubations.

GENERAL SIGNIFICANCE

Our models indicate that coupled structural contributions of the ingress, egress and solvent channels to the CYP3A4 active site geometries play key roles in the observed 4OHT binding modes. Thus adequate sampling of the conformational space of these drug-metabolizing promiscuous enzymes is important for substrates that may bind in malleable regions of the enzyme active-site.

Quantum mechanically derived AMBER-compatible heme parameters for various states of the cytochrome P450 catalytic cycle

Molecular mechanics (MM) methods are computationally affordable tools for screening chemical libraries of novel compounds for sites of P450 metabolism. One challenge for MM methods has been the absence of a consistent and transferable set of parameters for the heme within the P450 active site. Experimental data indicate that mammalian P450 enzymes vary greatly in the size, architecture, and plasticity of their active sites. Thus, obtaining X-ray-based geometries for the development of accurate MM parameters for the major classes of hepatic P450 remains a daunting task. Our previous work with preliminary gas-phase quantum mechanics (QM)-derived atomic partial charges greatly improved the accuracy of docking studies of raloxifene to CYP3A4. We have therefore developed and tested a consistent set of transferable MM parameters based on gas-phase QM calculations of two model systems of the heme-a truncated (T-HM) and a full (F-HM) for four states of the P450 catalytic cycle. Our results indicate that the use of the atomic partial charges from the F-HM further improves the accuracy of docked predictions for raloxifene to CYP3A4. Different patterns for substrate docking are also observed depending on the choice of heme model and state. Newly parameterized heme models are tested in implicit and explicitly solvated MD simulations in the absence and presence of enzyme structures, for CYP3A4, and appear to be stable on the nanosecond simulation timescale. The new force field for the various heme states may aid the community for simulations of P450 enzymes and other heme-containing enzymes.

Metabolism of 4-Aminopiperidine Drugs by Cytochrome P450s: Molecular and Quantum Mechanical Insights into Drug Design

4-Aminopiperidines are a variety of therapeutic agents that are extensively metabolized by cytochrome P450s with CYP3A4 as a major isoform catalyzing their N-dealkylation reaction. However, its catalytic mechanism has not been fully elucidated in a molecular interaction level. Here, we applied theoretical approaches including the molecular mechanics-based docking to study the binding patterns and quantum mechanics-based reactivity calculations. They were supported by the experimental human liver microsomal clearance and P450 isoform phenotyping data. Our results herein suggested that the molecular interactions between substrates and CYP3A4 active site residues are essential for the N-dealkylation of 4-aminopiperidines. We also found that the serine 119 residue of CYP3A4 may serve as a key hydrogen-bonding partner to interact with the 4-amino groups of the studied drugs. The reactivity of the side chain α-carbon hydrogens drives the direction of catalysis as well. As a result, structure-based drug design approaches look promising to guide drug discovery programs into the optimized drug metabolism space.

Substituted imidazole of 5-fluoro-2-[4-[(2-phenyl-1H-imidazol-5-yl)methyl]-1-piperazinyl]pyrimidine Inactivates cytochrome P450 2D6 by protein adduction

5-Fluoro-2-[4-[(2-phenyl-1H-imidazol-5-yl)methyl]-1-piperazinyl]pyrimidine (SCH 66712) is a potent mechanism-based inactivator of human cytochrome P450 2D6 that displays type I binding spectra with a K(s) of 0.39 ± 0.10 μM. The partition ratio is ~3, indicating potent inactivation that addition of exogenous nucleophiles does not prevent. Within 15 min of incubation with SCH 66712 and NADPH, ∼90% of CYP2D6 activity is lost with only ~20% loss in ability to bind CO and ~25% loss of native heme over the same time. The stoichiometry of binding to the protein was 1.2:1. SDS-polyacrylamide gel electrophoresis with Western blotting and autoradiography analyses of CYP2D6 after incubations with radiolabeled SCH 66712 further support the presence of a protein adduct. Metabolites of SCH 66712 detected by mass spectrometry indicate that the phenyl group on the imidazole ring of SCH 66712 is one site of oxidation by CYP2D6 and could lead to methylene quinone formation. Three other metabolites were also observed. For understanding the metabolic pathway that leads to CYP2D6 inactivation, metabolism studies with CYP2C9 and CYP2C19 were performed because neither of these enzymes is significantly inhibited by SCH 66712. The metabolites formed by CYP2C9 and CYP2C19 are the same as those seen with CYP2D6, although in different abundance. Modeling studies with CYP2D6 revealed potential roles of various active site residues in the oxidation of SCH 66712 and inactivation of CYP2D6 and showed that the phenyl group of SCH 66712 is positioned at 2.2 Å from the heme iron.