Functional analysis of phenylalanine residues in the active site of cytochrome P450 2C9

In silico characterization of the psilocybin biosynthesis pathway

Nearly all mushrooms of the Psilocybe genus contain the natural product psilocybin, which is a psychoactive alkaloid derived from l-tryptophan. Considering their use in ancient times, as well as their psychedelic properties, these mushrooms have re-emerged with psychotherapeutic potential for treating depression, which has triggered increased pharmaceutical interest. However, the psilocybin biosynthesis pathway was only recently defined and, as such, little exists in the way of structural data. Accordingly, the aim of this study was to structurally characterize this pathway by generating homology models for the four Psilocybe cubensis enzymes involved in psilocybin biosynthesis (PsiD, a decarboxylase; PsiH, a monooxygenase; PsiK, a phosphotransferase; PsiM, a methyltransferase). Following initial model generation and alignment with the identified structural templates, repeated refinement of the models was carried out using secondary structure prediction, geometry evaluation, energy minimization, and molecular dynamics simulations in water. The final models were then evaluated using molecular docking interactions with their substrates, i.e., psilocybin precursors (l-tryptophan, tryptamine, 4-hydroxytryptamine, and norbaeocystin/baeocystin), all of which generated feasible binding modes for the expected biotransformation. Further plausibility of the psilocybin → aeruginascin, 4-hydroxytryptamine → norpsilocin, and tryptamine → N,N-dimethyltryptamine conversions, all mediated by the generated model for PsiM, suggests valid routes of formation for these key secondary metabolites. The structural characterization of these enzymes and their binding modes which emerged from this study can lead to a better understanding of psilocybin synthesis, thereby paving the way for the development of novel substrates and selective inhibitors, as well as improved biotechnological manipulation and production of psilocybin in vitro.

Heterologous Expression and Functional Characterization of Novel CYP2C9 Variants Identified in the Alaska Native People

CYP2C9 is a major form of human liver cytochrome P450 that is responsible for the oxidative metabolism of several widely used low-therapeutic index drugs, including (S)-warfarin and phenytoin. In a cohort of Alaska Native people, ultrarare or novel CYP2C9 protein variants, M1L (rs114071557), N218I (rs780801862), and P279T (rs182132442, CYP2C9*29), are expressed with higher frequencies than the well characterized CYP2C9*2 and CYP2C9*3 alleles. We report here on their relative expression in lentivirus-infected HepG2 cells and the functional characterization of purified reconstituted enzyme variants expressed in Escherichia coli toward (S)-warfarin, phenytoin, flurbiprofen, and (S)-naproxen. In the infected HepG2 cells, robust mRNA and protein expression were obtained for wild-type, N218I, and P279T variants, but as expected, the M1L variant protein was not translated in this liver-derived cell line. His-tagged wild-type protein and the N218I and P279T variants, but not M1L, expressed well in E. coli and were highly purified after affinity chromatography. Upon reconstitution with cytochrome P450 oxidoreductase and cytochrome b5, the N218I and P279T protein variants metabolized (S)-warfarin, phenytoin, flurbiprofen, and (S)-naproxen to the expected monohydroxylated or O-demethylated metabolites. Steady-state kinetic analyses revealed that the relative catalytic efficiency ratios of (S)-warfarin metabolism by the P279T and N218I variants were 87% and 24%, respectively, of wild-type CYP2C9 protein. A similar rank ordering was observed for metabolism of phenytoin, flurbiprofen, and (S)-naproxen. We conclude that carriers of the variant N218I and, especially, the M1L alleles would be at risk of exacerbated therapeutic effects from drugs that rely on CYP2C9 for their metabolic clearance. SIGNIFICANCE STATEMENT: Novel gene variants of CYP2C9-M1L, and N218I, along with P279T (CYP2C9*29)-are expressed in Alaska Native people at relatively high frequencies. In vitro characterization of their functional effects revealed that each variant confers reduced catalytic efficiency toward several substrates, including the low-therapeutic index drugs (S)-warfarin and phenytoin. These data provide the first functional information for new, common CYP2C9 variants in this understudied population. The data may help guide dose adjustments in allele carriers, thus mitigating potential healthcare disparities.

The relationships between cytochromes P450 and H2O2: Production, reaction, and inhibition

In this review we address the relationship between cytochromes P450 (P450) and H2O2. This association can affect biology in three distinct ways. First, P450s produce H2O2 as a byproduct either during catalysis or when no substrate is present. This reaction, known as uncoupling, releases reactive oxygen species that may have implications in disease. Second, H2O2 is used as an oxygen-donating co-substrate in peroxygenase and peroxidase reactions catalyzed by P450s. This activity has proven to be important mainly in reactions involving prokaryotic P450s, and investigators have harnessed this reaction with the aim of adaptation for industrial use. Third, H2O2-dependent inhibition of human P450s has been studied in our laboratory, demonstrating heme destruction and also the inactivating oxidation of the heme-thiolate ligand to a sulfenic acid (-SOH). This reversible oxidative modification of P450s may have implications in the prevention of uncoupling and may give new insights into the oxidative regulation of these enzymes. Research has elucidated many of the chemical mechanisms involved in the relationship between P450 and H2O2, but the application to biology is difficult to evaluate. Further studies are needed reveal both the harmful and protective natures of reactive oxygen species in an organismal context.

Structural Basis of Single-Nucleotide Polymorphisms in Cytochrome P450 2C9

Single-nucleotide polymorphisms in drug-metabolizing cytochrome P450 (CYP) enzymes are important contributors to interindividual differences in drug metabolism leading to adverse drug reactions. Despite their extensive characterization and importance in pharmacogenetics of clinical drugs, the structural basis of CYP polymorphisms has remained scant. Here we report the crystal structures of human CYP2C9 and its polymorphic variants, *3 (I359L) and *30 (A477T), with an antihypertensive drug losartan. The structures show distinct interaction and occupation of losartan in the active site, the access channel, and the peripheral binding site. The I359L substitution located far from the active site remarkably altered the residue side chains near the active site and the access channel, whereas the T477 substitution illustrated hydrogen-bonding interaction with the reoriented side chain of Q214. The results yield structural insights into the reduced catalytic activity of the CYP2C9 variants and have important implications for understanding genetic polymorphisms in CYP-mediated drug metabolism.

Electrochemically driven drug metabolism via a CYP1A2-UGT1A10 bienzyme confined in a graphene nano-cage

A graphene nano-cage with regulatable space for the assembly of a cytochrome P450 1A2-UDP-glucuronosyltransferase 1A10 bienzyme complex has been constructed via a click reaction, and successfully used to study drug sequential metabolism using an electrochemically-driven method.

Mechanistic studies of the cationic binding pocket of CYP2C9 in vitro and in silico: metabolism of nonionizable analogs of tienilic acid

Tienilic acid (TA) is selectively oxidized at the C-5 position of the thiophene ring by the human liver enzyme cytochrome P450 2C9 (CYP2C9). This oxidation is mediated by the proximal positioning of the thiophene over the heme iron, which is proposed to be coordinated by an interaction of the TA carboxylic acid to a cationic binding pocket in the enzyme active site. In this study, we investigated how chemical modification of TA influences the bioactivation by CYP2C9. For this investigation, nine analogs of TA were chosen with substitutions on either side of the molecule. We tested three parameters, including CYP2C9 inhibition, metabolic profiling, and in silico docking. Of the 10 compounds tested, only two (TA and a noncarboxyl analog) resulted in competitive and time-dependent inhibition of CYP2C9. Metabolic profiling revealed a trend in which substitution of the carboxylate with nonionizable functional groups resulted in metabolic switching from oxidation of the aromatic ring to dealkylation reactions at the opposite side of the structure. The in silico modeling predicted an opposite binding orientation to that of TA for many analogs, including the 3-thenoyl regio-isomer analog, which contradicts previous models. Together these data show that disrupting interactions with the cationic binding pocket of CYP2C9 will impact the sites of metabolism and inhibition of the enzyme.

Preferred binding orientations of phenacetin in CYP1A1 and CYP1A2 are associated with isoform-selective metabolism

Human cytochromes P450 1A1 and 1A2 play important roles in drug metabolism and chemical carcinogenesis. Although these two enzymes share high sequence identity, they display different substrate specificities and inhibitor susceptibilities. In the present studies, we investigated the structural basis for these differences with phenacetin as a probe using a number of complementary approaches, such as enzyme kinetics, stoichiometric assays, NMR, and molecular modeling. Kinetic and stoichiometric analyses revealed that substrate specificity (k(cat)/K(m)) of CYP1A2 was approximately 18-fold greater than that of CYP1A1, as expected. Moreover, despite higher H₂O₂ production, the coupling efficiency of reducing equivalents to acetaminophen formation in CYP1A2 was tighter than that in CYP1A1. CYP1A1, in contrast to CYP1A2, displayed much higher uncoupling, producing more water. The subsequent NMR longitudinal (T₁) relaxation studies with the substrate phenacetin and its product acetaminophen showed that both compounds displayed similar binding orientations within the active site of CYP1A1 and CYP1A2. However, the distance between the OCH₂ protons of the ethoxy group (site of phenacetin O-deethylation) and the heme iron was 1.5 Å shorter in CYP1A2 than in CYP1A1. The NMR findings are thus consistent with our kinetic and stoichiometric results, providing a likely molecular basis for more efficient metabolism of phenacetin by CYP1A2.

Structure and dynamics of the membrane-bound cytochrome P450 2C9

The microsomal, membrane-bound, human cytochrome P450 (CYP) 2C9 is a liver-specific monooxygenase essential for drug metabolism. CYPs require electron transfer from the membrane-bound CYP reductase (CPR) for catalysis. The structural details and functional relevance of the CYP-membrane interaction are not understood. From multiple coarse grained molecular simulations started with arbitrary configurations of protein-membrane complexes, we found two predominant orientations of CYP2C9 in the membrane, both consistent with experiments and conserved in atomic-resolution simulations. The dynamics of membrane-bound and soluble CYP2C9 revealed correlations between opening and closing of different tunnels from the enzyme's buried active site. The membrane facilitated the opening of a tunnel leading into it by stabilizing the open state of an internal aromatic gate. Other tunnels opened selectively in the simulations of product-bound CYP2C9. We propose that the membrane promotes binding of liposoluble substrates by stabilizing protein conformations with an open access tunnel and provide evidence for selective substrate access and product release routes in mammalian CYPs. The models derived here are suitable for extension to incorporate other CYPs for oligomerization studies or the CYP reductase for studies of the electron transfer mechanism, whereas the modeling procedure is generally applicable to study proteins anchored in the bilayer by a single transmembrane helix.

We describe the first atomic-detail models and simulations of a full-length, membrane-bound mammalian cytochrome P450. To date, all the structural studies of microsomal, drug-metabolizing cytochrome P450s have been performed using engineered, solubilized forms of the enzymes and it is not yet understood how the membrane influences their structure, dynamics, and ability to bind substrates. We focused on CYP2C9, the second most abundant cytochrome P450 in the human liver which oxidizes 20% of all marketed drugs. Here, we have derived models of CYP2C9-membrane complexes with a modeling procedure based on molecular dynamics simulations started with arbitrary configurations of the protein in the membrane and performed using both coarse grained and atomic-detail molecular representations. This procedure is expected to be generally applicable to proteins that are anchored in the membrane with a single transmembrane helix. The models and simulations provide evidence for selective substrate access and product release routes in membrane-bound CYPs. This observation may contribute to new strategies to manipulate the activity of CYPs and other enzymes with buried active sites. We anticipate that this study will bring about a paradigm shift towards studying microsomal CYPs as dynamic structures in their natural, lipid environment rather than in artificially solubilized forms.

Intramolecular heme ligation of the cytochrome P450 2C9 R108H mutant demonstrates pronounced conformational flexibility of the B-C loop region: implications for substrate binding

A previous study [Dickmann, L., et al. (2004) Mol. Pharmacol. 65, 842-850] revealed some unusual properties of the R108H mutant of cytochrome P450 2C9 (CYP2C9), including elevated thermostability relative to that of CYP2C9, as well as a UV-visible absorbance spectrum that was indicative of nitrogenous ligation to the heme iron. In our study, size-exclusion chromatography and UV-visible absorbance spectroscopy of CYP2C9 R108H monomers demonstrated that nitrogen ligation is indeed intramolecular. Pulsed electron paramagnetic resonance of CYP2C9 R108H monomers showed that a histidine is most likely bound to the heme as previously hypothesized. An energy-minimized model of the R108H mutant maintained a CYP fold, despite substantial movement of several loop regions of the mutant, and, therefore, represents an extreme example of a closed conformation of the enzyme. Molecular dynamics (MD) simulations of CYP2C9 were performed to study the range of energetically accessible CYP2C9 conformations. These in silico studies showed that the B-C loop region of CYP2C9 moves away from the heme to a position resembling the putative open conformation described for rabbit CYP2B4. A model involving the movement of the B-C loop region and R108 between the open and closed conformations of CYP2C9 is presented, which helps to explain the enzyme's ability to regio- and stereospecifically metabolize some ligands while allosterically activating others.

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.

CYP2C9 amino acid residues influencing phenytoin turnover and metabolite regio- and stereochemistry

Phenytoin has been an effective anticonvulsant agent for over 60 years, although its clinical use is complicated by nonlinear pharmacokinetics, a narrow therapeutic index, and metabolically based drug-drug interactions. Although it is well established that CYP2C9 is the major cytochrome P450 enzyme controlling metabolic elimination of phenytoin through its oxidative conversion to (S)-5-(4-hydroxyphenyl)-5-phenylhydantoin (p-HPPH), nothing is known about the amino acid binding determinants within the CYP2C9 active site that promote metabolism and maintain the tight stereocontrol of hydroxy metabolite formation. This knowledge gap was addressed here through the construction of nine active site mutants at amino acid positions Phe100, Arg108, Phe114, Leu208, and Phe476 and in vitro analysis of the steady-state kinetics and stereochemistry of p-HPPH formation. The F100L and F114W mutants exhibited 4- to 5-fold increases in catalytic efficiency, whereas the F100W, F114L, F476L, and F476W mutants lost >90% of their phenytoin hydroxylation capacity. This pattern of effects differs substantially from that found previously for (S)-warfarin and (S)-flurbiprofen metabolism, suggesting that these three ligands bind within discrete locations in the CYP2C9 active site. Only the F114L, F476L, and L208V mutants altered phenytoin's orientation during catalytic turnover. The L208V mutant also uniquely demonstrated enhanced 6-hydroxylation of (S)-warfarin. These latter data provide the first experimental evidence for a role of the F-G loop region in dictating the catalytic orientation of substrates within the CYP2C9 active site.