Roles of phenylalanine at position 120 and glutamic acid at position 222 in the oxidation of chiral substrates by cytochrome P450 2D6

Tryptophan-75 Is a Low-Energy Channel-Gating Residue that Facilitates Substrate Egress/Access in Cytochrome P450 2D6

CYP2D6 is a major drug metabolizing enzyme with a buried active site. Channels leading to the active site from various enzyme surfaces are believed to facilitate ligand egress and access to the active site. The present study used molecular dynamics (MD) and in vitro studies with CYP2D6*1 and a Trp75-to-Ala mutant to examine channel gating in CYP2D6 by Trp75. MD simulations measured energy landscapes of Trp75 conformations and simulated substrate passage within channel 2b using bufuralol as a model substrate. Trp75 alternated between multiple stable states that supported substrate transport along channel 2b with low-energy barriers between states (∼ -1 kcal/mol). Trp75 conformations were stabilized primarily by hydrogen bonding between Trp75 and Glu222, Asn226, Ala225, or Gln72. Energy barriers were low between Trp75 conformations, allowing Trp75 to easily move between various conformations over time and to function in both binding to and moving substrates in the 2b channel of CYP2D6. Michaelis-Menten kinetic studies completed with purified enzyme in a reconstituted system showed overall reduced enzyme efficiency for metabolism of bufuralol and dextromethorphan by the Trp75Ala mutant compared with CYP2D6*1. In stopped-flow measurements, k off for dextromethorphan was decreased in the absence of Trp75. Our results support a role for Trp75 in substrate shuttling to the active site of CYP2D6. SIGNIFICANCE STATEMENT: Using combined molecular dynamics and in vitro assays, this study shows for the first time a role for Trp75 as a channel entrance gating residue in the mechanism of substrate binding/unbinding in CYP2D6. Energy landscapes derived from molecular dynamics were used to quantitate the strength of gating, and kinetics assays showed the impact on enzyme efficiency and k off of a Trp75Ala mutation.

Deciphering Reaction Determinants of Altered-Activity CYP2D6 Variants by Well-Tempered Metadynamics Simulation and QM/MM Calculations

The xenobiotic metabolizing enzyme CYP2D6 is the P450 cytochrome family member with the highest rate of polymorphism. This causes changes in the enzyme activity and specificity, which can ultimately lead to adverse reactions during drug treatment. To avoid or lower CYP-related toxicity risks, prediction of the most likely positions within a molecule where a metabolic reaction might occur is paramount. In order to obtain accurate predictions, it is crucial to understand all phenomena within the active site of the enzyme that contribute to an efficient substrate recognition and the subsequent catalytic reaction together with their relative weight within the overall thermodynamic context. This study aims to define the weight of the driving forces upon the C-H bond activation within CYP2D6 wild-type and a clinically relevant allelic variant with increased activity (CYP2D6*53) featuring two amino acid mutations in close vicinity of the heme. First, we investigated the steric and electrostatic complementarity of the substrate bufuralol using well-tempered metadynamics simulations with the aim to obtain the free energy profiles for each site of metabolism (SoM) within the different active sites. Second, the stereoelectronic complementarity was determined for each SoM within the two different active-site environments. Relying on the well-tempered metadynamics simulation energy profiles of each SoM, we identified the binding mode that was closest to the preferred transition-state geometry for efficient C-H bond activation. The binding modes were then used as starting structures for the quantum mechanics/molecular mechanics calculations performed to quantify the corresponding activation barriers. Our results show the relevance of the steric component in orienting the SoM in an energetically accessible position toward the heme. However, the corresponding intrinsic reactivity and electronic complementarity within the active site must be accurately evaluated in order to obtain a meaningful reaction prediction, from which the predominant SoM can be determined. The F120I mutation lowered the activation barrier for the major site and one of the minor SoMs. However, it had an impact neither on the CYP2D6 enantioselectivity preference of the oxidation reaction nor on the stereoselectivity from the substrate point of view.

New insights into the structural characteristics and functional relevance of the human cytochrome P450 2D6 enzyme

To date, the crystal structures of at least 12 human CYPs (1A2, 2A6, 2A13, 2C8, 2C9, 2D6, 2E1, 2R1, 3A4, 7A1, 8A1, and 46A1) have been determined. CYP2D6 accounts for only a small percentage of all hepatic CYPs (< 2%), but it metabolizes approximately 25% of clinically used drugs with significant polymorphisms. CYP2D6 also metabolizes procarcinogens and neurotoxins, such as 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine, 1,2,3,4-tetrahydroquinoline, and indolealkylamines. Moreover, the enzyme utilizes hydroxytryptamines and neurosteroids as endogenous substrates. Typical CYP2D6 substrates are usually lipophilic bases with an aromatic ring and a nitrogen atom, which can be protonated at physiological pH. Substrate binding is generally followed by oxidation (5-7 A) from the proposed nitrogen-Asp301 interaction. A number of homology models have been constructed to explore the structural features of CYP2D6, while antibody studies also provide useful structural information. Site-directed mutagenesis studies have demonstrated that Glu216, Asp301, Phe120, Phe481, and Phe483 play important roles in determining the binding of ligands to CYP2D6. The structure of human CYP2D6 has been recently determined and shows the characteristic CYP fold observed for other members of the CYP superfamily. The lengths and orientations of the individual secondary structural elements in the CYP2D6 structure are similar to those seen in other human CYP2 members, such as CYP2C9 and 2C8. The 2D6 structure has a well-defined active-site cavity located above the heme group with a volume of approximately 540 A(3), which is larger than equivalent cavities in CYP2A6 (260 A(3)), 1A2 (375 A(3)), and 2E1 (190 A(3)), but smaller than those in CYP3A4 (1385 A(3)) and 2C8 (1438 A(3)). Further studies are required to delineate the molecular mechanisms involved in CYP2D6 ligand interactions and their implications for drug development and clinical practice.

The mechanism causing the difference in kinetic properties between rat CYP2D4 and human CYP2D6 in the oxidation of dextromethorphan and bufuralol

The capacity to oxidize bufuralol (BF) and dextromethorphan (DEX) was compared kinetically between human CYP2D6 and four rat CYP2D (CYP2D1, -2D2, -2D3 and -2D4) isoenzymes in a yeast cell expression system. In BF 1''-hydroxylation and DEX O-demethylation, only CYP2D4 showed hook-shaped Eadie-Hofstee plots, the other four CYP2D enzymes exhibiting linear plots. In DEX N-demethylation, rat CYP2D2 did not show any detectable activity under the conditions used, whereas the other four enzymes yielded linear Eadie-Hofstee plots. To elucidate the mechanisms causing the nonlinear kinetics, four CYP2D4 mutants, CYP2D4-F109I, -V123F, -L216F and -A486F, were prepared. CYP2D4-V123F, -L216F and -A486F yielded linear or linear-like Eadie-Hofstee plots for BF 1''-hydroxylation, whereas only CYP2D4-A486F exhibited linear plots for DEX O-demethylation. The substitution of Phe-109 by isoleucine did not have any effect on the oxidative capacity of CYP2D4 for either BF or DEX. These results suggest that the introduction of phenylalanine in the active-site cavity of CYP2D4 simplifies complicated interactions between the substrates and the amino acid residues, but the mechanisms causing the simplification differ between BF and DEX.

Analysis of CYP2D6 substrate interactions by computational methods

Cytochrome P450 CYP2D6 is involved in the oxidation of well over 150 drugs and, in general, those which contain a basic nitrogen atom in the molecule. To clarify how the residues of CYP2D6 are utilized for orientating a wide range of its specific substrates and distinguishing them from a variety of other organic compounds, docking studies by AutoDock and molecular dynamics (MD) simulations were conducted. Specific ligands were docked to both the homology model and crystal structures optimally to estimate the site of reaction on the ligand molecule and the binding energy for the complex, which were generally in good agreement with the experimental data. MD simulation for the CYP2D6-propranolol complex was then carried out to reveal the amino acid residues interacting with the substrate at the active site. Phe-120, Glu-216, Asp-301, and Phe-483 are identified as the substrate-binding residues in agreement with previously reported site-directed mutagenesis data and the crystal structure reported recently (PDB code: 2F9Q). As well as these residues, our theoretical prediction suggests that Phe-219 and Glu-222 are also important residues for mediating oxidation of substrates, especially propranolol.

Stereoselective hexobarbital 3'-hydroxylation by CYP2C19 expressed in yeast cells and the roles of amino acid residues at positions 300 and 476

We examined the enzymatic function of recombinant CYP2C19 in enantiomeric hexobarbital (HB) 3'-hydroxylation, and searched the roles of amino acid residues, such as Phe-100, Phe-114, Asp-293, Glu-300, and Phe-476 of CYP2C19 in the stereoselective HB 3'-hydroxylation, using a yeast cell expression system and site-directed mutagenesis method. CYP2C19 wild-type exerted substrate enantioselectivity of (R)-HB>>(S)-HB and metabolite diastereoselectivity of 3'(R)<3'(S) in 3'-hydroxylation of HB enantiomers. The substitution of Asp-293 by alanine failed to yield an observable peak at 450 nm in its reduced carbon monoxide-difference spectrum. CYP2C19-E300A and CYP2C19-E300V with alanine and valine, respectively, in place of Glu-300 exerted total HB 3'-hydroxylation activities of 45 and 108%, respectively, that of the wild-type. Interestingly, these two mutants showed substrate enantioselectivity of (R)-HB<(S)-HB, which is opposite to that of the wild-type, while metabolite diasteroselectivity remained unchanged. The replacement of Phe-476 by alanine increased total HB 3'-hydroxylation activity to approximately 3-fold that of the wild-type. Particularly, 3'(S)-OH-(S)-HB-forming activity elevated to 7-fold that of the wild-type, resulting in the reversal of the substrate enantioselectivity. In contrast, the substitution of phenylalanine at positions 100 and 114 by alanine did not produce a remarkable change in the total activity or the substrate enantioselectivity. These results indicate that Glu-300 and Phe-476 are important in stereoselective oxidation of HB enantiomers by CYP2C19.

Generation of a homology model for the human cytochrome P450, CYP24A1, and the testing of putative substrate binding residues by site-directed mutagenesis and enzyme activity studies

A systematic analysis of conserved H-bonding patterns and tertiary structural motifs from 13 crystal structures was used to create a homology model for the human multicatalytic cytochrome P450, CYP24A1, involved in catabolism of 1alpha,25-dihydroxyvitamin D3. The substrate was docked in the active site and used to identify potential substrate contact residues in the B' helix, B'/C loop, F-helix and the beta-5 hairpin. Seven CYP24A1 mutants were created and studied by mammalian cell transfection and CYP24A1 activity assay. Mutants showed reduced metabolic rates and altered metabolite patterns compared to wild-type. We conclude that: Ile-131 positions substrate via A-ring and cis-triene contacts; Trp-134 and Gly-499 are determinants of substrate access; Leu-148 contacts the substrate side-chain; Met-246 is important in mediating regioselectivity. Our findings validate the new model of CYP24A1, which can now be used to predict structural modifications for rational vitamin D drug design.

Cloning of a cDNA encoding a novel marmoset CYP2C enzyme, expression in yeast cells and characterization of its enzymatic functions

We cloned a cDNA encoding a novel CYP2C enzyme, called P450 M-2C, from a marmoset liver. The deduced amino acid sequence showed high identities to those of human CYP2C8 (87%), CYP2C9 (78%) and CYP2C19 (77%). The P450 M-2C enzyme expressed in yeast cells catalyzed p-methylhydroxylation of only tolbutamide among four substrates tested, paclitaxel as a CYP2C8 substrate, diclofenac and tolbutamide as CYP2C9 substrates and S-mephenytoin as a CYP2C19 substrate. p-Methylhydroxylation of tolbutamide by marmoset liver microsomes showed monophasic kinetics, and the apparent K(m) value (1.2 mM) for the substrate was similar to that of the recombinant P450 M-2C (1.8 mM). Although all of the recombinant human CYP2C8, CYP2C9 and CYP2C19 expressed in yeast cells catalyzed tolbutamide p-methylhydroxylation, the kinetic profile of CYP2C8 was most similar to that of P450 M-2C. Tolbutamide oxidation by the marmoset liver microsomes and the recombinant P450 M-2C was inhibited most effectively by quercetin, a CYP2C8 inhibitor, followed by omeprazole, a CYP2C19 inhibitor, whereas sulfaphenazole, a CYP2C9 inhibitor, was less potent under the conditions used. These results indicate that P450 M-2C is the major tolbutamide p-methylhydroxylase in the marmoset liver.