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

Functional significance of CYP2B6 gene rare allelic variants identified in Japanese individuals

Cytochrome P450 2B6 (CYP2B6) catalyzes the metabolism of many drugs, including efavirenz and propofol. Genetic polymorphisms in CYP2B6 alter its enzymatic activity and substantially affect its pharmacokinetics. High-frequency variants, such as CYP2B6*6, are associated with the risk of developing side effects due to reduced CYP2B6 activity. However, the impact of rare alterations on enzyme function remains unknown, and some of these variants may significantly decrease the CYP2B6 activity. Therefore, in this study, we evaluated in vitro the functional alterations in 29 missense variants of the CYP2B6 gene identified in 8,380 Japanese individuals. Wild-type CYP2B6 and 29 rare CYP2B6 variants were transiently expressed in mammalian cells. The expression levels of variant CYP2B6 proteins in the microsomal fractions extracted from 293FT cells were assessed using western blotting and reduced-carbon monoxide difference spectroscopy, and a specific peak at 450 nm was detected in the wild-type and 19 variants. Furthermore, kinetic parameters were determined by assaying the reactions with efavirenz and propofol and quantifying the metabolite concentrations. We found that 12 variants had significantly lower or abolished enzymatic activity with both the substrates. In silico three-dimensional docking and molecular-dynamics simulations suggested that these functional changes were due to conformational changes in essential regions, such as the heme-binding site and ligand channels involved in transporting substrates to the active site. These findings have implications for predicting the plasma concentrations of CYP2B6 substrates and controlling their side effects.

Massively parallel characterization of CYP2C9 variant enzyme activity and abundance

CYP2C9 encodes a cytochrome P450 enzyme responsible for metabolizing up to 15% of small molecule drugs, and CYP2C9 variants can alter the safety and efficacy of these therapeutics. In particular, the anti-coagulant warfarin is prescribed to over 15 million people annually and polymorphisms in CYP2C9 can affect individual drug response and lead to an increased risk of hemorrhage. We developed click-seq, a pooled yeast-based activity assay, to test thousands of variants. Using click-seq, we measured the activity of 6,142 missense variants in yeast. We also measured the steady-state cellular abundance of 6,370 missense variants in a human cell line by using variant abundance by massively parallel sequencing (VAMP-seq). These data revealed that almost two-thirds of CYP2C9 variants showed decreased activity and that protein abundance accounted for half of the variation in CYP2C9 function. We also measured activity scores for 319 previously unannotated human variants, many of which may have clinical relevance.

Multiple drug binding modes in Mycobacterium tuberculosis CYP51B1

The Mycobacterium tuberculosis (Mtb) genome encodes 20 different cytochrome P450 enzymes (CYPs), many of which serve essential biosynthetic roles. CYP51B1, the Mtb version of eukaryotic sterol demethylase, remains a potential therapeutic target. The binding of three drug fragments containing nitrogen heterocycles to CYP51B1 is studied here by continuous wave (CW) and pulsed electron paramagnetic resonance (EPR) techniques to determine how each drug fragment binds to the heme active-site. All three drug fragments form a mixture of complexes, some of which retain the axial water ligand from the resting state. Hyperfine sublevel correlation spectroscopy (HYSCORE) and electron-nuclear double resonance spectroscopy (ENDOR) observe protons of the axial water and on the drug fragments that reveal drug binding modes. Binding in CYP51B1 is complicated by the presence of multiple binding modes that coexist in the same solution. These results aid our understanding of CYP-inhibitor interactions and will help guide future inhibitor design.

Induced fit for cytochrome P450 3A4 based on molecular dynamics

The present study aims at numerically describing to what extent substrate - enzyme complexes in solution may change over time as a natural process of conformational changes for a liganded enzyme in comparison to those movements which occur independently from substrate interaction, i.e. without a ligand. To this end, we selected structurally known pairs of liganded / unliganded CYP450 3A4 enzymes with different geometries hinting at induced fit events. We carried out molecular dynamics simulations (MD) comparing the trajectories in a "cross-over" protocol: (i) we added the ligand to the unliganded crystal form which should adopt geometries similar to the known geometry of the liganded crystal structure during MD, and - conversely - (ii) we removed the bound ligand form the known liganded complex to test if a geometry similar to the known unliganded (apo-) form can be adopted during MD. To compare continues changes we measured root means square deviations and frequencies. Results for case (i) hint at larger conformational changes required for accepting the substrate during its approach to final position - in contrast to case (ii) when mobility is fairly reduced by ligand binding (strain energy). In conclusion, a larger conformational sampling prior to ligand binding and the freezing-in (rigidity) of conformations for bound ligands can be interpreted as two conditions linked to induced-fit.

Insights into molecular mechanisms of drug metabolism dysfunction of human CYP2C9*30

Cytochrome P450 2C9 (CYP2C9) metabolizes about 15% of clinically administrated drugs. The allelic variant CYP2C9*30 (A477T) is associated to diminished response to the antihypertensive effects of the prodrug losartan and affected metabolism of other drugs. Here, we investigated molecular mechanisms involved in the functional consequences of this amino-acid substitution. Molecular dynamics (MD) simulations performed for the active species of the enzyme (heme in the Compound I state), in the apo or substrate-bound state, and binding energy analyses gave insights into altered protein structure and dynamics involved in the defective drug metabolism of human CYP2C9.30. Our data revealed an increased rigidity of the key Substrate Recognition Sites SRS1 and SRS5 and shifting of the β turn 4 of SRS6 toward the helix F in CYP2C9.30. Channel and binding substrate dynamics analyses showed altered substrate channel access and active site accommodation. These conformational and dynamic changes are believed to be involved in the governing mechanism of the reduced catalytic activity. An ensemble of representative conformations of the WT and A477T mutant properly accommodating drug substrates were identified, those structures can be used for prediction of new CYP2C9 and CYP2C9.30 substrates and drug-drug interactions.

CW EPR parameters reveal cytochrome P450 ligand binding modes

Cytochrome P450 (CYP) monoxygenses utilize heme cofactors to catalyze oxidation reactions. They play a critical role in metabolism of many classes of drugs, are an attractive target for drug development, and mediate several prominent drug interactions. Many substrates and inhibitors alter the spin state of the ferric heme by displacing the heme's axial water ligand in the resting enzyme to yield a five-coordinate iron complex, or they replace the axial water to yield a nitrogen-ligated six-coordinate iron complex, which are traditionally assigned by UV-vis spectroscopy. However, crystal structures and recent pulsed electron paramagnetic resonance (EPR) studies find a few cases where molecules hydrogen bond to the axial water. The water-bridged drug-H₂O-heme has UV-vis spectra similar to nitrogen-ligated, six-coordinate complexes, but are closer to "reverse type I" complexes described in older liteature. Here, pulsed and continuous wave (CW) EPR demonstrate that water-bridged complexes are remarkably common among a range of nitrogenous drugs or drug fragments that bind to CYP3A4 or CYP2C9. Principal component analysis reveals a distinct clustering of CW EPR spectral parameters for water-bridged complexes. CW EPR reveals heterogeneous mixtures of ligated states, including multiple directly-coordinated complexes and water-bridged complexes. These results suggest that water-bridged complexes are under-represented in CYP structural databases and can have energies similar to other ligation modes. The data indicates that water-bridged binding modes can be identified and distinguished from directly-coordinated binding by CW EPR.

Pharmacogenomics of CYP2C9: Functional and Clinical Considerations

CYP2C9 is the most abundant CYP2C subfamily enzyme in human liver and the most important contributor from this subfamily to drug metabolism. Polymorphisms resulting in decreased enzyme activity are common in the CYP2C9 gene and this, combined with narrow therapeutic indices for several key drug substrates, results in some important issues relating to drug safety and efficacy. CYP2C9 substrate selectivity is detailed and, based on crystal structures for the enzyme, we describe how CYP2C9 catalyzes these reactions. Factors relevant to clinical response to CYP2C9 substrates including inhibition, induction and genetic polymorphism are discussed in detail. In particular, we consider the issue of ethnic variation in pattern and frequency of genetic polymorphisms and clinical implications. Warfarin is the most well studied CYP2C9 substrate; recent work on use of dosing algorithms that include CYP2C9 genotype to improve patient safety during initiation of warfarin dosing are reviewed and prospects for their clinical implementation considered. Finally, we discuss a novel approach to cataloging the functional capabilities of rare 'variants of uncertain significance', which are increasingly detected as more exome and genome sequencing of diverse populations is conducted.

Longitudinal characterisation of function and structure of Bietti crystalline dystrophy: report on a novel homozygous mutation in CYP4V2

BACKGROUND

Bietti crystalline dystrophy (BCD) is a rare inherited disorder characterised by fine crystalline deposits in the corneal limbus and retinal posterior pole. In 2004, mutations in the CYP4V2 gene were identified as the cause of BCD. Here, we describe the report of a homozygous point mutation in a patient with BCD and provide detailed characterisation of functional and structural changes over 20 years.

METHODS

At regular intervals, the patient underwent repeat ophthalmic evaluations. DNA was extracted from buccal swabs, amplified by standard PCR and analysed for homology to the CYP4V2 sequence. Homology modelling was conducted using Iterative Threading ASSEmbly Refinement and molecular dynamics simulations using GROningen MAchine for Chemical Simulations.

RESULTS

The proband, a 47-year-old woman of German ancestry was diagnosed with crystalline retinopathy at age 25. Over the next 20 years, visual acuity and function gradually declined with progression of retinal pigment epithelium and choroidal atrophy. When first tested at 39 years of age, the multifocal electroretinogram (ERG) was markedly abnormal, more so for the right eye, whereas the full-field ERG was more symmetrical and lagged other measures of visual function. Gene sequencing showed a single C>T point mutation in exon 9 encoding a R400C amino acid change. Computational modelling suggests the mutation impairs function due to loss of a hydrogen bonding interaction with the propionate side chains of the haeme prosthetic group.

CONCLUSION

This is the first report of a homozygous R400C mutation in CYP4V2 with protein modelling showing high likelihood of enzyme dysfunction. The comprehensive long-term clinical follow-up provides insight into disease progression and highlights possible anti-inflammatory modulation of disease severity.

The CYP79A1 catalyzed conversion of tyrosine to (E)-p-hydroxyphenylacetaldoxime unravelled using an improved method for homology modeling

The vast diversity and membrane-bound nature of plant P450s makes it challenging to study the structural characteristics of this class of enzymes especially with respect to accurate intermolecular enzyme-substrate interactions. To address this problem we here apply a modified hybrid structure strategy for homology modeling of plant P450s. This allows for structural elucidation based on conserved motifs in the protein sequence and secondary structure predictions. We modeled the well-studied Sorghum bicolor cytochrome P450 CYP79A1 catalyzing the first step in the biosynthesis of the cyanogenic glucoside dhurrin. Docking experiments identified key regions of the active site involved in binding of the substrate and facilitating catalysis. Arginine 152 and threonine 534 were identified as key residues interacting with the substrate. The model was validated experimentally using site-directed mutagenesis. The new CYP79A1 model provides detailed insights into the mechanism of the initial steps in cyanogenic glycoside biosynthesis. The approach could guide functional characterization of other membrane-bound P450s and provide structural guidelines for elucidation of key structure-function relationships of other plant P450s.

Physical Studies of P450-P450 Interactions: Predicting Quaternary Structures of P450 Complexes in Membranes from Their X-ray Crystal Structures

Cytochrome P450 enzymes, which catalyze oxygenation reactions of both exogenous and endogenous chemicals, are membrane bound proteins that require interaction with their redox partners in order to function. Those responsible for drug and foreign compound metabolism are localized primarily in the endoplasmic reticulum of liver, lung, intestine, and other tissues. More recently, the potential for P450 enzymes to exist as supramolecular complexes has been shown by the demonstration of both homomeric and heteromeric complexes. The P450 units in these complexes are heterogeneous with respect to their distribution and function, and the interaction of different P450s can influence P450-specific metabolism. The goal of this review is to examine the evidence supporting the existence of physical complexes among P450 enzymes. Additionally, the review examines the crystal lattices of different P450 enzymes derived from X-ray diffraction data to make assumptions regarding possible quaternary structures in membranes and in turn, to predict how the quaternary structures could influence metabolism and explain the functional effects of specific P450-P450 interactions.

The tetrahydrobiopterin radical interacting with high- and low-spin heme in neuronal nitric oxide synthase - A new indicator of the extent of NOS coupling

Reaction intermediates trapped during the single-turnover reaction of the neuronal ferrous nitric oxide synthase oxygenase domain (Fe(II)nNOSOX) show four EPR spectra of free radicals. Fully-coupled nNOSOX with cofactor (tetrahydrobiopterin, BH4) and substrate (l-arginine) forms the typical BH4 cation radical with an EPR spectrum ~4.0mT wide and hyperfine tensors similar to reports for a biopterin cation radical in inducible NOSOX (iNOSOX). With excess thiol, nNOSox lacking BH4 and l-arg is known to produce superoxide. In contrast, we find that nNOSOX with BH4 but no l-arg forms two radicals with rather different, fast (~250μs at 5K) and slower (~500μs at 20K), electron spin relaxation rates and a combined ~7.0mT wide EPR spectrum. Rapid freeze-quench CW- and pulsed-EPR measurements are used to identify these radicals and their origin. These two species are the same radical with identical nuclear hyperfine couplings, but with spin-spin couplings to high-spin (4.0mT component) or low-spin (7.0mT component) Fe(III) heme. Uncoupled reactions of nNOS leave the enzyme in states that can be chemically reduced to sustain unregulated production of NO and reactive oxygen species in ischemia-reperfusion injury. The broad EPR signal is a convenient indicator of uncoupled nNOS reactions producing low-spin Fe(III) heme.

Immobilized Cytochrome P450 for Monitoring of P450-P450 Interactions and Metabolism

Cytochrome P450 (P450) protein-protein interactions have been shown to alter their catalytic activity. Furthermore, these interactions are isoform specific and can elicit activation, inhibition, or no effect on enzymatic activity. Studies show that these effects are also dependent on the protein partner cytochrome P450 reductase (CPR) and the order of protein addition to purified reconstituted enzyme systems. In this study, we use controlled immobilization of P450s to a gold surface to gain a better understanding of P450-P450 interactions between three key drug-metabolizing isoforms (CYP2C9, CYP3A4, and CYP2D6). Molecular modeling was used to assess the favorability of homomeric/heteromeric P450 complex formation. P450 complex formation in vitro was analyzed in real time utilizing surface plasmon resonance. Finally, the effects of P450 complex formation were investigated utilizing our immobilized platform and reconstituted enzyme systems. Molecular modeling shows favorable binding of CYP2C9-CPR, CYP2C9-CYP2D6, CYP2C9-CYP2C9, and CYP2C9-CYP3A4, in rank order.KDvalues obtained via surface plasmon resonance show strong binding, in the nanomolar range, for the above pairs, with CYP2C9-CYP2D6 yielding the lowestKD, followed by CYP2C9-CYP2C9, CYP2C9-CPR, and CYP2C9-CYP3A4. Metabolic incubations show that immobilized CYP2C9 metabolism was activated by homomeric complex formation. CYP2C9 metabolism was not affected by the presence of CYP3A4 with saturating CPR concentrations. CYP2C9 metabolism was activated by CYP2D6 at saturating CPR concentrations in solution but was inhibited when CYP2C9 was immobilized. The order of addition of proteins (CYP2C9, CYP2D6, CYP3A4, and CPR) influenced the magnitude of inhibition for CYP3A4 and CYP2D6. These results indicate isoform-specific P450 interactions and effects on P450-mediated metabolism.

Drug modulation of water-heme interactions in low-spin P450 complexes of CYP2C9d and CYP125A1

Azoles and pyridines are commonly incorporated into small molecule inhibitor scaffolds that target cytochromes P450 (CYPs) as a strategy to increase drug binding affinity, impart isoform-dependent selectivity, and improve metabolic stability. Optical absorbance spectra of the CYP-inhibitor complex are widely used to infer whether these inhibitors are ligated directly to the heme iron as catalytically inert, low-spin (type II) complexes. Here, we show that the low-spin complex between a drug-metabolizing CYP2C9 variant and 4-(3-phenylpropyl)-1H-1,2,3-triazole (PPT) retains an axial water ligand despite exhibiting elements of "classic" type II optical behavior. Hydrogens of the axial water ligand are observed by pulsed electron paramagnetic resonance (EPR) spectroscopy for both inhibitor-free and inhibitor-bound species and show that inhibitor binding does not displace the axial water. A (15)N label incorporated into PPT is 0.444 nm from the heme iron, showing that PPT is also in the active site. The reverse type I inhibitor, LP10, of CYP125A1 from Mycobacterium tuberculosis, known from X-ray crystal structures to form a low-spin water-bridged complex, is found by EPR and by visible and near-infrared magnetic circular dichroism spectroscopy to retain the axial water ligand in the complex in solution.

Soluble and membrane-bound Drosophila melanogaster CYP6G1 expressed in Escherichia coli: purification, activity, and binding properties toward multiple pesticides

Cytochrome P450 CYP6G1 has been implicated in the resistance of Drosophila melanogaster to numerous pesticides. While in vivo and in vitro studies have provided insight to the diverse functions of this enzyme, direct studies on the isolated CYP6G1 enzyme have not been possible due to the need for a source of recombinant enzyme. In the current study, the Cyp6g1 gene was isolated from D. melanogaster and re-engineered for heterologous expression in Escherichia coli. Approximately 460 nmol L⁻¹ of P450 holoenzyme were obtained in 500 mL cultures. The recombinant enzyme was located predominantly within the bacterial cytosol. A two-step purification protocol using Ni-chelate affinity chromatography followed by removal of detergent on a hydroxyapatite column produced essentially homogenous enzyme from both soluble and membrane fractions. Recombinant CYP6G1 exhibited p-nitroanisole O-dealkylation activity but was not active against eleven other typical P450 marker substrates. Substrate-induced binding spectra and IC₅₀ values for inhibition of p-nitroanisole O-dealkylation were obtained for a wide selection of pesticides, namely DDT, imidacloprid, chlorfenvinphos, malathion, endosulfan, dieldrin, dicyclanil, lufenuron and carbaryl, supporting previous in vivo and in vitro studies on Drosophila that have suggested that the enzyme is involved in multi-pesticide resistance in insects.

1,2,3-Triazole-heme interactions in cytochrome P450: functionally competent triazole-water-heme complexes

In comparison to imidazole (IMZ) and 1,2,4-triazole (1,2,4-TRZ), the isosteric 1,2,3-triazole (1,2,3-TRZ) is unrepresented among cytochrome P450 (CYP) inhibitors. This is surprising because 1,2,3-TRZs are easily obtained via "click" chemistry. To understand this underrepresentation of 1,2,3-TRZs among CYP inhibitors, thermodynamic and density functional theory computational studies were performed with unsubstituted IMZ, 1,2,4-TRZ, and 1,2,3-TRZ. The results indicate that the lower affinity of 1,2,3-TRZ for the heme iron includes a large unfavorable entropy term likely originating in solvent-1,2,3-TRZ interactions; the difference is not solely due to differences in the enthalpy of heme-ligand interactions. In addition, the 1,2,3-TRZ fragment was incorporated into a well-established CYP3A4 substrate and mechanism-based inactivator, 17-α-ethynylestradiol (17EE), via click chemistry. This derivative, 17-click, yielded optical spectra consistent with low-spin ferric heme iron (type II) in contrast to 17EE, which yields a high-spin complex (type I). Furthermore, the rate of CYP3A4-mediated metabolism of 17-click was comparable to that of 17EE, with a different regioselectivity. Surprisingly, continuous-wave electron paramagnetic resonance (EPR) and HYSCORE EPR spectroscopy indicate that 17-click does not displace water from the sixth axial ligand position of CYP3A4 as expected for a type II ligand. We propose a binding model in which 17-click pendant 1,2,3-TRZ hydrogen bonds with the sixth axial water ligand. The results demonstrate the potential for 1,2,3-TRZ to form metabolically labile water-bridged low-spin heme complexes, consistent with recent evidence that nitrogenous type II ligands of CYPs can be efficiently metabolized. The specific case of [CYP3A4·17-click] highlights the risk of interpreting CYP-ligand complex structure on the basis of optical spectra.

Computational prediction of metabolism: sites, products, SAR, P450 enzyme dynamics, and mechanisms

Metabolism of xenobiotics remains a central challenge for the discovery and development of drugs, cosmetics, nutritional supplements, and agrochemicals. Metabolic transformations are frequently related to the incidence of toxic effects that may result from the emergence of reactive species, the systemic accumulation of metabolites, or by induction of metabolic pathways. Experimental investigation of the metabolism of small organic molecules is particularly resource demanding; hence, computational methods are of considerable interest to complement experimental approaches. This review provides a broad overview of structure- and ligand-based computational methods for the prediction of xenobiotic metabolism. Current computational approaches to address xenobiotic metabolism are discussed from three major perspectives: (i) prediction of sites of metabolism (SOMs), (ii) elucidation of potential metabolites and their chemical structures, and (iii) prediction of direct and indirect effects of xenobiotics on metabolizing enzymes, where the focus is on the cytochrome P450 (CYP) superfamily of enzymes, the cardinal xenobiotics metabolizing enzymes. For each of these domains, a variety of approaches and their applications are systematically reviewed, including expert systems, data mining approaches, quantitative structure–activity relationships (QSARs), and machine learning-based methods, pharmacophore-based algorithms, shape-focused techniques, molecular interaction fields (MIFs), reactivity-focused techniques, protein–ligand docking, molecular dynamics (MD) simulations, and combinations of methods. Predictive metabolism is a developing area, and there is still enormous potential for improvement. However, it is clear that the combination of rapidly increasing amounts of available ligand- and structure-related experimental data (in particular, quantitative data) with novel and diverse simulation and modeling approaches is accelerating the development of effective tools for prediction of in vivo metabolism, which is reflected by the diverse and comprehensive data sources and methods for metabolism prediction reviewed here. This review attempts to survey the range and scope of computational methods applied to metabolism prediction and also to compare and contrast their applicability and performance.

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.