Structure of mammalian cytochrome P450 2C5 complexed with diclofenac at 2.1 A resolution: evidence for an induced fit model of substrate binding

P450 variations bifurcate the early terpene indole alkaloid pathway in Catharanthus roseus and Camptotheca acuminata

Divergent terpene indole alkaloid (TIA) pathways in Catharanthus roseus and Camptotheca acuminata generate vinblastine and vincristine, and camptothecin, respectively. In contrast to Catharanthus which feeds secologanin (from methylated loganin) into its species-specific late pathway, Camptotheca feeds secologanic acid (from unmethylated loganic acid) into its late pathway. Having identified putative Camptotheca secologanic acid synthases (SLASs) and cytochrome P450 reductases (CPRs) in transcriptome databases, we have demonstrated that two P450s, CYP72A564 and CYP72A565, are capable of utilizing both loganic acid and loganin to generate secologanic acid and secologanin. We have extended the previous report of these activities by CYP72A565 and CYP72A610 (Yang et al., 2019) by demonstrating that both Arabidopsis CPRs (ATR1, ATR2) couple with these CYP72A proteins in yeast microsomal assays and that purified Camptotheca CPR1 couples with them in in vitro reconstitution assays. Kinetic analyses of purified full-length Camptotheca SLASs have indicated that both process loganic acid with nearly identical catalytic rates and efficiencies as measured by their k_cat and k_cat/K_M. In contrast, CYP72A564 processes loganin with two-fold greater efficiency than CYP72A565 correlating with the former's 3-fold greater affinity for loganin. The closely-related CYP72A730 does not bind or process either compound. Molecular modeling of these three proteins and comparisons with Catharanthus secologanin synthase (SLS) have identified key differences that likely determine their SLAS versus SLS selectivities. Our ability to reconstitute these SLAS/SLS activities provides valuable tools for further examinations of the residues involved in substrate recognition and determinations of their unusual mechanism of C-C bond scission.

CYP2J2 Molecular Recognition: A New Axis for Therapeutic Design

Cytochrome P450 (CYP) epoxygenases are a special subset of heme-containing CYP enzymes capable of performing the epoxidation of polyunsaturated fatty acids (PUFA) and the metabolism of xenobiotics. This dual functionality positions epoxygenases along a metabolic crossroad. Therefore, structure-function studies are critical for understanding their role in bioactive oxy-lipid synthesis, drug-PUFA interactions, and for designing therapeutics that directly target the epoxygenases. To better exploit CYP epoxygenases as therapeutic targets, there is a need for improved understanding of epoxygenase structure-function. Of the characterized epoxygenases, human CYP2J2 stands out as a potential target because of its role in cardiovascular physiology. In this review, the early research on the discovery and activity of epoxygenases is contextualized to more recent advances in CYP epoxygenase enzymology with respect to PUFA and drug metabolism. Additionally, this review employs CYP2J2 epoxygenase as a model system to highlight both the seminal works and recent advances in epoxygenase enzymology. Herein we cover CYP2J2's interactions with PUFAs and xenobiotics, its tissue-specific physiological roles in diseased states, and its structural features that enable epoxygenase function. Additionally, the enumeration of research on CYP2J2 identifies the future needs for the molecular characterization of CYP2J2 to enable a new axis of therapeutic design.

Inhibitory Effect of AB-PINACA, Indazole Carboxamide Synthetic Cannabinoid, on Human Major Drug-Metabolizing Enzymes and Transporters

Indazole carboxamide synthetic cannabinoid, AB-PINACA, has been placed into Schedule I of the Controlled Substances Act by the US Drug Enforcement Administration since 2015. Despite the possibility of AB-PINACA exposure in drug abusers, the interactions between AB-PINACA and drug-metabolizing enzymes and transporters that play crucial roles in the pharmacokinetics and efficacy of various substrate drugs have not been investigated. This study was performed to investigate the inhibitory effects of AB-PINACA on eight clinically important human major cytochrome P450s (CYPs) and six uridine 5′-diphospho-glucuronosyltransferases (UGT) in human liver microsomes and the activities of six solute carrier transporters and two efflux transporters in transporter-overexpressing cells. AB-PINACA reversibly inhibited the metabolic activities of CYP2C8 (K_i, 16.9 µM), CYP2C9 (K_i, 6.7 µM), and CYP2C19 (K_i, 16.1 µM) and the transport activity of OAT3 (K_i, 8.3 µM). It exhibited time-dependent inhibition on CYP3A4 (K_i, 17.6 µM; k_inact, 0.04047 min⁻¹). Other metabolizing enzymes and transporters such as CYP1A2, CYP2A6, CYP2B6, CYP2D6, UGT1A1, UGT1A3, UGT1A4, UGT1A6, UGT1A9, UGT2B7, OAT1, OATP1B1, OATP1B3, OCT1, OCT2, P-glycoprotein, and BCRP, exhibited only weak interactions with AB-PINACA. These data suggest that AB-PINACA can cause drug-drug interactions with CYP3A4 substrates but that the significance of drug interactions between AB-PINACA and CYP2C8, CYP2C9, CYP2C19, or OAT3 substrates should be interpreted carefully.

Single-Molecule Mechanistic Study of Enzyme Hysteresis

Hysteresis is an important feature of enzyme-catalyzed reactions, as it reflects the influence of enzyme regulation in the presence of ligands such as substrates or allosteric molecules. In typical kinetic studies of enzyme activity, hysteretic behavior is observed as a "lag" or "burst" in the time course of the catalyzed reaction. These lags and bursts are due to the relatively slow transition from one state to another state of the enzyme molecule, with different states having different kinetic properties. However, it is difficult to understand the underlying mechanism of hysteresis by observing bulk reactions because the different enzyme molecules in the population behave stochastically. In this work, we studied the hysteretic behavior of mutant β-glucuronidase (GUS) using a high-throughput single-molecule array platform and investigated the effect of thermal treatment on the hysteresis.

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.

Specificity and mechanism of carbohydrate demethylation by cytochrome P450 monooxygenases

Degradation of carbohydrates by bacteria represents a key step in energy metabolism that can be inhibited by methylated sugars. Removal of methyl groups, which is critical for further processing, poses a biocatalytic challenge because enzymes need to overcome a high energy barrier. Our structural and computational analysis revealed how a member of the cytochrome P450 family evolved to oxidize a carbohydrate ligand. Using structural biology, we ascertained the molecular determinants of substrate specificity and revealed a highly specialized active site complementary to the substrate chemistry. Invariance of the residues involved in substrate recognition across the subfamily suggests that they are critical for enzyme function and when mutated, the enzyme lost substrate recognition. The structure of a carbohydrate-active P450 adds mechanistic insight into monooxygenase action on a methylated monosaccharide and reveals the broad conservation of the active site machinery across the subfamily.

Conformational Change Induced by Putidaredoxin Binding to Ferrous CO-ligated Cytochrome P450cam Characterized by 2D IR Spectroscopy

The importance of conformational dynamics to protein function is now well-appreciated. An outstanding question is whether they are involved in the effector role played by putidaredoxin (Pdx) in its reduction of the O₂ complex of cytochrome P450cam (P450cam), an archetypical member of the cytochrome P450 superfamily. Recent studies have reported that binding of Pdx induces a conformational change from a closed to an open state of ferric P450cam, but a similar conformational change does not appear to occur for the ferrous, CO-ligated enzyme. To better understand the effector role of Pdx when binding the ferrous, CO-ligated P450cam, we applied 2D IR spectroscopy to compare the conformations and dynamics of the wild-type (wt) enzyme in the absence and presence of Pdx, as well as of L358P P450cam (L358P), which has served as a putative model for the Pdx complex. The CO vibrations of the Pdx complex and L358P report population of two conformational states in which the CO experiences distinct environments. The dynamics among the CO frequencies indicate that the energy landscape of substates within one conformation are reflective of the closed state of P450cam, and for the other conformation, differ from the free wt enzyme, but are equivalent between the Pdx complex and L358P. The two states co-populated by the Pdx complex are postulated to reflect a loosely bound encounter complex and a more tightly bound state, as is commonly observed for the dynamic complexes of redox partners. Significantly, this study shows that the binding of Pdx to ferrous, CO-ligated P450cam does perturb the conformational ensemble in a way that might underlie the effector role of Pdx.

Kauniolide synthase is a P450 with unusual hydroxylation and cyclization-elimination activity

Guaianolides are an important class of sesquiterpene lactones with unique biological and pharmaceutical properties. They have been postulated to be derived from germacranolides, but for years no progress has been made in the elucidation of their biosynthesis that requires an unknown cyclization mechanism. Here we demonstrate the isolation and characterization of a cytochrome P450 from feverfew (Tanacetum parthenium), kauniolide synthase. Kauniolide synthase catalyses the formation of the guaianolide kauniolide from the germacranolide substrate costunolide. Unlike most cytochrome P450s, kauniolide synthase combines stereoselective hydroxylation of costunolide at the C3 position, with water elimination, cyclization and regioselective deprotonation. This unique mechanism of action is supported by in silico modelling and docking experiments. The full kauniolide biosynthesis pathway is reconstructed in the heterologous hosts Nicotiana benthamiana and yeast, paving the way for biotechnological production of guaianolide-type sesquiterpene lactones.

Molecular dynamics investigations of regioselectivity of anionic/aromatic substrates by a family of enzymes: a case study of diclofenac binding in CYP2C isoforms

The CYP2C subfamily is of particular importance in the metabolism of drugs, food toxins, and procarcinogens. Like other P450 subfamilies, 2C enzymes share a high sequence identity, but significantly contribute in different ways to hepatic capacity to metabolize drugs. They often metabolize the same substrate to more than one product with different catalytic sites. Because it is challenging to characterize experimentally, much still remains unknown about the reason for why the substrate regioselectivity of these closely related subfamily members is different. Here, we have investigated the structural features of CYP2C8, CYP2C9, and CYP2C19 bound with their shared substrate diclofenac to elucidate the underlying molecular mechanism for the substrate regioselectivity of CYP2C subfamily enzymes. The obtained results demonstrate how a sequence divergence for the active site residues causes heterogeneous variations in the secondary structures and in major tunnel selections, and further affects the shape and chemical properties of the substrate-binding site. Structural analysis and free energy calculations showed that the most important determinants of regioselectivity among the CYP2C isoforms are the geometrical features of the active sites, as well as the hydrogen bonds and the hydrophobic interactions, mainly presenting as the various locations of Arg108 and substitutions of Phe205 for Ile205 in CYP2C8. The MM-GB/SA calculations combined with PMF results accord well with the experimental KM values, bridging the gap between the theory and the experimentally observed results of binding affinity differences. The present study provides important insights into the structure-function relationships of CYP2C subfamily enzymes, the knowledge of ligand binding characteristics and key residue contributions could guide future experimental and computational work on the synthesis of drugs with better pharmacokinetic properties so that CYP interactions could be avoided.

Homology Modeling and Molecular Dynamics Simulation Combined with X-ray Solution Scattering Defining Protein Structures of Thromboxane and Prostacyclin Synthases

A combination of molecular dynamics (MD) simulations and X-ray scattering (SAXS) has emerged as the approach of choice for studying protein structures and dynamics in solution. This approach has potential applications for membrane proteins that neither are soluble nor form crystals easily. We explore the water-coupled dynamic structures of thromboxane synthase (TXAS) and prostacyclin synthase (PGIS) from scanning HPLC-SAXS measurements combined with MD ensemble analyses. Both proteins are heme-containing enzymes in the cytochrome P450 family, known as prostaglandin H₂ (PGH₂) isomerase, with counter-functions in regulation of platelet aggregation. Currently, the X-ray crystallographic structures of PGIS are available, but those for TXAS are not. The use of homology modeling of the TXAS structure with ns-μs explicit water solvation MD simulations allows much more accurate estimation of the configuration space with loop motion and origin of the protein behaviors in solution. In contrast to the stability of the conserved PGIS structure in solution, the pronounced TXAS flexibility has been revealed to have unstructured loop regions in connection with the characteristic P450 structural elements. The MD-derived and experimental-solution SAXS results are in excellent agreement. The significant protein internal motions, whole-molecule structures, and potential problems with protein folding, crystallization, and functionality are examined.

Functioning of drug-metabolizing microsomal cytochrome P450s: In silico probing of proteins suggests that the distal heme 'active site' pocket plays a relatively 'passive role' in some enzyme-substrate interactions

PURPOSE

The currently held mechanistic understanding of microsomal cytochrome P450s (CYPs) seeks that diverse drug molecules bind within the deep-seated distal heme pocket and subsequently react at the heme centre. To explain a bevy of experimental observations and meta-analyses, we indulge a hypothesis that involves a "diffusible radical mediated" mechanism. This new hypothesis posits that many substrates could also bind at alternate loci on/within the enzyme and be reacted without the pertinent moiety accessing a bonding proximity to the purported catalytic Fe-O enzyme intermediate.

METHODS

Through blind and heme-distal pocket centered dockings of various substrates and non-substrates (drug molecules of diverse sizes, classes, topographies etc.) of microsomal CYPs, we explored the possibility of access of substrates via the distal channels, its binding energies, docking orientations, distance of reactive moieties (or molecule per se) to/from the heme centre, etc. We investigated specific cases like- (a) large drug molecules as substrates, (b) classical marker drug substrates, (c) class of drugs as substrates (Sartans, Statins etc.), (d) substrate preferences between related and unrelated CYPs, (e) man-made site-directed mutants' and naturally occurring mutants' reactivity and metabolic disposition, (f) drug-drug interactions, (g) overall affinities of drug substrate versus oxidized product, (h) meta-analysis of in silico versus experimental binding constants and reaction/residence times etc.

RESULTS

It was found that heme-centered dockings of the substrate/modulator drug molecules with the available CYP crystal structures gave poor docking geometries and distances from Fe-heme centre. In conjunction with several other arguments, the findings discount the relevance of erstwhile hypothesis in many CYP systems. Consequently, the newly proposed hypothesis is deemed a viable alternate, as it satisfies Occam's razor.

CONCLUSIONS

The new proposal affords expanded scope for explaining the mechanism, kinetics and overall phenomenology of CYP mediated drug metabolism. It is now understood that the heme-iron and the hydrophobic distal pocket of CYPs serve primarily to stabilize the reactive intermediate (diffusible radical) and the surface or crypts of the apoprotein bind to the xenobiotic substrate (and in some cases, the heme distal pocket could also serve the latter function). Thus, CYPs enhance reaction rates and selectivity/specificity via a hitherto unrecognized modality.

The Halogen Bond

The halogen bond occurs when there is evidence of a net attractive interaction between an electrophilic region associated with a halogen atom in a molecular entity and a nucleophilic region in another, or the same, molecular entity. In this fairly extensive review, after a brief history of the interaction, we will provide the reader with a snapshot of where the research on the halogen bond is now, and, perhaps, where it is going. The specific advantages brought up by a design based on the use of the halogen bond will be demonstrated in quite different fields spanning from material sciences to biomolecular recognition and drug design.

Discovery of a novel binding pocket for CYP 2C9 inhibitors: crystallography, pharmacophore modelling and inhibitor SAR

CYP 2C9 co-crystal structures of compound1(green) in binding mode 1 (CYP 2C9 in green) and compound2(yellow) in binding mode 2 (CYP 2C9 in yellow).

Computational insights into the photocyclization of diclofenac in solution: effects of halogen and hydrogen bonding

The effects of noncovalent interactions, namely halogen and hydrogen bonding, on the photochemical conversion of the photosensitizing drug diclofenac (DCF) in solution were investigated computationally.

The Schistosoma mansoni Cytochrome P450 (CYP3050A1) Is Essential for Worm Survival and Egg Development

Schistosomiasis affects millions of people in developing countries and is responsible for more than 200,000 deaths annually. Because of toxicity and limited spectrum of activity of alternatives, there is effectively only one drug, praziquantel, available for its treatment. Recent data suggest that drug resistance could soon be a problem. There is therefore the need to identify new drug targets and develop drugs for the treatment of schistosomiasis. Analysis of the Schistosoma mansoni genome sequence for proteins involved in detoxification processes found that it encodes a single cytochrome P450 (CYP450) gene. Here we report that the 1452 bp open reading frame has a characteristic heme-binding region in its catalytic domain with a conserved heme ligating cysteine, a hydrophobic leader sequence present as the membrane interacting region, and overall structural conservation. The highest sequence identity to human CYP450s is 22%. Double stranded RNA (dsRNA) silencing of S. mansoni (Sm)CYP450 in schistosomula results in worm death. Treating larval or adult worms with antifungal azole CYP450 inhibitors results in worm death at low micromolar concentrations. In addition, combinations of SmCYP450-specific dsRNA and miconazole show additive schistosomicidal effects supporting the hypothesis that SmCYP450 is the target of miconazole. Treatment of developing S. mansoni eggs with miconazole results in a dose dependent arrest in embryonic development. Our results indicate that SmCYP450 is essential for worm survival and egg development and validates it as a novel drug target. Preliminary structure-activity relationship suggests that the 1-(2,4-dichlorophenyl)-2-(1H-imidazol-1-yl)ethan-1-ol moiety of miconazole is necessary for activity and that miconazole activity and selectivity could be improved by rational drug design.

Over 600 million people in endemic countries are at risk of contracting schistosomiasis, which results in over 200,000 deaths each year and significant illness to most people that are infected. There are concerns that the drug widely used for the treatment of schistosomiasis, praziquantel, may be losing efficacy due to evolution of drug resistant worms. Since the disease mainly affects the poor in developing countries, pharmaceutical companies have little interest in developing new drugs and none are currently being tested. In this paper we focus on a novel parasite protein, cytochrome P450, which we propose to be a new drug target. Worms are unusual in having only one cytochrome P450 gene; humans have 57 cytochrome P450 genes. By using reverse genetic and chemical approaches we found that the schistosome cytochrome P450 is essential for worm survival and egg development and, therefore, is an essential and druggable target. Drugs that target fungal cytochrome P450s and are already in use for treating several human diseases were identified as potential hits for further development for schistosomiasis treatment.

Topical diclofenac does not affect the antiplatelet properties of aspirin as compared to the intermediate effects of oral diclofenac: A prospective, randomized, complete crossover study

Nonsteroidal anti-inflammatory drugs (NSAIDs) adversely interact with aspirin, diminishing its antiplatelet effect and potentially placing patients at an increased risk for recurrent thrombotic events. This crossover study aimed to determine whether the topical NSAID diclofenac epolamine 1.3% patch or oral diclofenac 50 mg interfered with the antiplatelet effects of aspirin 325 mg. Twelve healthy men and women aged 18-50 were included. Participants were randomized into 5 treatment arms: aspirin, diclofenac potassium 50 mg, diclofenac patch, diclofenac potassium plus ASA 325 mg, and diclofenac patch plus aspirin. Platelet responsiveness was determined using whole-blood impedance aggregation (WBA) to collagen 1 μg/mL and arachidonic acid (AA) 0.5 mM and was sampled every 2 hours. No significant difference in platelet function was observed following the diclofenac patch and aspirin vs aspirin alone. Oral diclofenac produced a mixed effect with significant reduction in platelet inhibition at hour 2 and hour 8 following aspirin administration. Topical diclofenac does not significantly interfere with the antiplatelet effects of aspirin and may be a safer alternative to the oral formulation.

Molecular basis for nonspecificity of nonsteroidal anti-inflammatory drugs (NSAIDs)

Inhibition of the production of inflammatory mediators by the action of nonsteroidal anti-inflammatory drugs (NSAIDs) is highly accredited to their recognition of cyclooxygenase enzymes. Along with inflammation relief, however, NSAIDs also cause adverse effects. Although NSAIDs strongly inhibit enzymes of the prostaglandin synthesis pathways, several other proteins also serve as fairly potent targets for these drugs. Based on their recognition pattern, these receptors are categorised as enzymes modifying NSAIDs, noncatalytic proteins binding to NSAIDs and enzymes with catalytic functions that are inhibited by NSAIDs. The extensive binding of NSAIDs is responsible for their limited in vivo efficacy as well as the large spectrum of their effects. The biochemical nature of drugs binding to multiple protein targets and its implications on physiology are discussed.

High warfarin sensitivity in carriers of CYP2C9*35 is determined by the impaired interaction with P450 oxidoreductase

Cytochrome P450 2C9 (CYP2C9) metabolizes many clinically important drugs including warfarin and diclofenac. We have recently reported a new allelic variant, CYP2C9*35, found in a warfarin hypersensitive patient with Arg125Leu and Arg144Cys mutations. Here, we have investigated the molecular basis for the functional consequences of these polymorphic changes. CYP2C9.1 and CYP2C9-Arg144Cys expressed in human embryonic kidney 293 cells effectively metabolized both S-warfarin and diclofenac in NADPH-dependent reactions, whereas CYP2C9-Arg125Leu or CYP2C9.35 were catalytically silent. However, when NADPH was replaced by a direct electron donor to CYPs, cumene hydroperoxide, hereby bypassing the CYP oxidoreductase (POR), all variant enzymes were active, indicating unproductive interactions between CYP2C9.35 and POR. In silico analysis revealed a decrease of the electrostatic potential of CYP2C9-Arg125Leu-POR interacting surface and the loss of stabilizing salt bridges between these proteins. In conclusion, our data strongly suggest that the Arg125Leu mutation in CYP2C9.35 prevents CYP2C9-POR interactions resulting in the absence of NADPH-dependent CYP2C9-catalyzed activity in vivo, thus influencing the warfarin sensitivity in the carriers of this allele.

Correlating structure and function of drug-metabolizing enzymes: progress and ongoing challenges

This report summarizes a symposium sponsored by the American Society for Pharmacology and Experimental Therapeutics at Experimental Biology held April 20-24 in Boston, MA. Presentations discussed the status of cytochrome P450 (P450) knowledge, emphasizing advances and challenges in relating structure with function and in applying this information to drug design. First, at least one structure of most major human drug-metabolizing P450 enzymes is known. However, the flexibility of these active sites can limit the predictive value of one structure for other ligands. A second limitation is our coarse-grain understanding of P450 interactions with membranes, other P450 enzymes, NADPH-cytochrome P450 reductase, and cytochrome b5. Recent work has examined differential P450 interactions with reductase in mixed P450 systems and P450:P450 complexes in reconstituted systems and cells, suggesting another level of functional control. In addition, protein nuclear magnetic resonance is a new approach to probe these protein/protein interactions, identifying interacting b5 and P450 surfaces, showing that b5 and reductase binding are mutually exclusive, and demonstrating ligand modulation of CYP17A1/b5 interactions. One desired outcome is the application of such information to control drug metabolism and/or design selective P450 inhibitors. A final presentation highlighted development of a CYP3A4 inhibitor that slows clearance of human immunodeficiency virus drugs otherwise rapidly metabolized by CYP3A4. Although understanding P450 structure/function relationships is an ongoing challenge, translational advances will benefit from continued integration of existing and new biophysical approaches.

[Metabolism of xenobiotics: beneficial and adverse effects]

The systems developed by living organisms for the metabolism of xenobiotics play a key role in the adaptation of living species to their chemical environment. Recent data about mammalian cytochrome P450 structures have led to a better understanding of the molecular basis for the adaptability of these enzymes to xenobiotics exhibiting highly variable structures. The action of these enzymes on xenobiotics leads to other beneficial effects such as the bioactivation of some drugs, but also to adverse effects with the formation of aggressive metabolites for the cell that are responsible for the appearance of many toxicities.

Structural diversity of eukaryotic membrane cytochrome p450s

X-ray crystal structures are available for 29 eukaryotic microsomal, chloroplast, or mitochondrial cytochrome P450s, including two non-monooxygenase P450s. These structures provide a basis for understanding structure-function relations that underlie their distinct catalytic activities. Moreover, structural plasticity has been characterized for individual P450s that aids in understanding substrate binding in P450s that mediate drug clearance.

Carboetomidate: an analog of etomidate that interacts weakly with 11β-hydroxylase

BACKGROUND

Carboetomidate is a pyrrole etomidate analog that is 3 orders of magnitude less potent an inhibitor of in vitro cortisol synthesis than etomidate (an imidazole) and does not inhibit in vivo steroid production. Although carboetomidate's reduced functional effect on steroid synthesis is thought to reflect lower binding affinity to 11β-hydroxylase, differential binding to this enzyme has never been experimentally demonstrated. In the current study, we tested the hypothesis that carboetomidate and etomidate bind with differential affinity to 11β-hydroxylase by comparing their abilities to inhibit photoaffinity labeling of purified enzyme by a photoactivatable etomidate analog and to modify the enzyme's absorption spectrum in a way that is indicative of ligand binding. In addition, we made a preliminary exploration of the manner in which etomidate and carboetomidate might differentially interact with this site using spectroscopic methods as well as molecular modeling techniques to better understand the structural basis for their selectivity.

METHODS

The ability of azi-etomidate to inhibit cortisol synthesis was tested by assessing its ability to inhibit cortisol synthesis by H295R cells. The binding affinities of etomidate and carboetomidate to 11β-hydroxylase were compared by assessing their abilities to (1) inhibit photoincorporation of the photolabile etomidate analog [(3)H]azi-etomidate into the enzyme and (2) modify the absorption spectrum of the enzyme's heme group. In silico docking studies of etomidate, carboetomidate, and azi-etomidate binding to 11β-hydroxylase were performed using the computer software GOLD.

RESULTS

Similar to etomidate, azi-etomidate potently inhibits in vitro cortisol synthesis. Etomidate inhibited [(3)H]azi-etomidate photolabeling of 11β-hydroxylase in a concentration-dependent manner. At a concentration of 40 µM, etomidate reduced photoincorporation of [(3)H]azi-etomidate by 96% ± 1% whereas carboetomidate had no experimentally detectable effect. On addition of etomidate to 11β-hydroxylase, a type 2 difference spectrum was produced indicative of etomidate complexation with the enzyme's heme iron; carboetomidate had no effect whereas azi-etomidate produced a reverse type 1 spectrum. Computer modeling studies predicted that etomidate, carboetomidate, and azi-etomidate can fit into the heme-containing pocket that forms 11β-hydroxylase's active site and pose with their carbonyl oxygens interacting with the heme iron and their phenyl rings stacking with phenylalanine-80. However, additional unique poses were identified for etomidate and azi-etomidate that likely account for their higher affinities.

CONCLUSIONS

Carboetomidate's reduced ability to suppress in vitro and in vivo steroid synthesis as compared with etomidate reflects its lower binding affinity to 11β-hydroxylase and may be attributed to carboetomidate's inability to form a coordination bond with the heme iron located at the enzyme's active site.

Structural characterization of human cytochrome P450 2C19: active site differences between P450s 2C8, 2C9, and 2C19

To identify the structural features underlying the distinct substrate and inhibitor profiles of P450 2C19 relative to the closely related human enzymes, P450s 2C8 and 2C9, the atomic structure (Protein Data Bank code 4GQS) of cytochrome P450 2C19 complexed with the inhibitor (2-methyl-1-benzofuran-3-yl)-(4-hydroxy-3,5-dimethylphenyl)methanone (Protein Data Bank chemical component 0XV) was determined to 2.87 Å resolution by x-ray crystallography. The conformation of the peptide backbone of P450 2C19 is most similar to that of P450 2C8, but the substrate-binding cavity of P450 2C8 is much larger than that of P450 2C19 due to differences in the amino acid residues that form the substrate-binding cavities of the two enzymes. In contrast, the substrate-binding cavity of P450 2C19 is much more similar in size to that of the structure of the P450 2C9 flurbiprofen complex than to that of a modified P450 2C9 or that of P450 2C8. The cavities of the P450 2C19 0XV complex and the P450 2C9 flurbiprofen complex differ, however, because the helix B-C loops of the two enzymes are dissimilar. These conformational differences reflect the effects of adjacent structural elements that interact with the B-C loops and that differ between the two enzymes. The availability of a structure for 2C19 will facilitate computational approaches for predictions of substrate and inhibitor binding to this enzyme.

Cytochrome P450 1A1-mediated anticancer drug discovery: in silico findings

INTRODUCTION

Target-specific drugs may offer fewer side/adverse effects in comparison with other anticancer agents and thus save normal healthy cells to a greater extent. The selective overexpression of cytochrome P450 1A1 (CYP1A1) in tumor cells induces the metabolism of benzothiazole and aminoflavone compounds to their reactive species, which are responsible for DNA adduct formation and cell death. This review encompasses the novelty of CYP1A1 as an anticancer drug target and explores the possible in silico strategies that would be applicable in the discovery and development of future antitumor compounds.

AREAS COVERED

This review highlights the various ligand-based and target-based in silico methodologies that were efficiently used in exploration of CYP1A1 as a novel antitumor target. These methodologies include electronic structure analysis, CoMFA studies, homology modeling, molecular docking, molecular dynamics analysis, pharmacophore mapping and quantitative structure activity relationship (QSAR) studies. It also focuses on the various approaches used in the development of the lysyl amide prodrug of 5F-203 (NSC710305) and dimethanesulfonate salt of 5-aminoflavone (NSC710464) as clinical candidates from their less potent analogues.

EXPERT OPINION

Selective overexpression of CYP1A1 in cancer cells offers tumor-specific drug design to ameliorate the current adverse effects associated with existing antitumor agents. Medicinal chemistry and in vitro driven approaches, in combination with knowledge-based drug design and by using the currently available tools of in silico methodologies, would certainly make it possible to design and develop novel anticancer compounds targeting CYP1A1.

Plasticity of CYP2B enzymes: structural and solution biophysical methods

In the past three years, major advances in understanding cytochrome P450 2B (CYP2B) structure-function relationships have been made through determination of multiple ligand-bound and one ligand-free X-ray crystal structure of CYP2B4 and one ligand-bound X-ray crystal structure of CYP2B6. These structures have provided insight into the features that provide the high degree of plasticity of the enzymes. A combination of a phenylalanine cluster that allows for concerted movement of helices F through G and a conserved set of electrostatic interactions involving Arg(262) facilitates movement of this region to accommodate binding of ligands of various sizes without perturbing most of the P450 fold. Integrating solution based techniques such as NMR or deuterium exchange mass spectrometry (DXMS) with computational methods including molecular docking has provided further insight into enzyme behavior upon ligand binding. In addition, extended molecular dynamics simulations have provided a link between an open and a closed conformation of ligand-free CYP2B4 found in crystal structures. Other studies revealed the utility of rational engineering in improving stability of P450s to facilitate structural studies. The solution and computational results combined with the X-ray crystal structures yield a comprehensive picture of how these enzymes adopt different conformations to bind various ligands.

Crystal structure of human cytochrome P450 2D6 with prinomastat bound

Human cytochrome P450 2D6 contributes to the metabolism of >15% of drugs used in clinical practice. This study determined the structure of P450 2D6 complexed with a substrate and potent inhibitor, prinomastat, to 2.85 Å resolution by x-ray crystallography. Prinomastat binding is well defined by electron density maps with its pyridyl nitrogen bound to the heme iron. The structure of ligand-bound P450 2D6 differs significantly from the ligand-free structure reported for the P450 2D6 Met-374 variant (Protein Data Bank code 2F9Q). Superposition of the structures reveals significant differences for β sheet 1, helices A, F, F', G", G, and H as well as the helix B-C loop. The structure of the ligand complex exhibits a closed active site cavity that conforms closely to the shape of prinomastat. The closure of the open cavity seen for the 2F9Q structure reflects a change in the direction and pitch of helix F and introduction of a turn at Gly-218, which is followed by a well defined helix F' that was not observed in the 2F9Q structure. These differences reflect considerable structural flexibility that is likely to contribute to the catalytic versatility of P450 2D6, and this new structure provides an alternative model for in silico studies of substrate interactions with P450 2D6.

Structural features of cytochromes P450 and ligands that affect drug metabolism as revealed by X-ray crystallography and NMR

Cytochromes P450 (P450s) play a major role in the clearance of drugs, toxins, and environmental pollutants. Additionally, metabolism by P450s can result in toxic or carcinogenic products. The metabolism of pharmaceuticals by P450s is a major concern during the design of new drug candidates. Determining the interactions between P450s and compounds of very diverse structures is complicated by the variability in P450-ligand interactions. Understanding the protein structural elements and the chemical attributes of ligands that dictate their orientation in the P450 active site will aid in the development of effective and safe therapeutic agents. The goal of this review is to describe P450-ligand interactions from two perspectives. The first is the various structural elements that microsomal P450s have at their disposal to assume the different conformations observed in X-ray crystal structures. The second is P450-ligand dynamics analyzed by NMR relaxation studies.

Identification of H-Ras-specific motif for the activation of invasive signaling program in human breast epithelial cells

Increased expression and/or activation of H-Ras are often associated with tumor aggressiveness in breast cancer. Previously, we showed that H-Ras, but not N-Ras, induces MCF10A human breast epithelial cell invasion and migration, whereas both H-Ras and N-Ras induce cell proliferation and phenotypic transformation. In an attempt to determine the sequence requirement directing the divergent phenotype induced by H-Ras and N-Ras with a focus on the induction of human breast cell invasion, we investigated the structural and functional relationships between H-Ras and N-Ras using domain-swap and site-directed mutagenesis approaches. Here, we report that the hypervariable region (HVR), consisting of amino acids 166 to 189 in H-Ras, determines the invasive/migratory signaling program as shown by the exchange of invasive phenotype by swapping HVR sequences between H-Ras and N-Ras. We also demonstrate that the H-Ras-specific additional palmitoylation site at Cys184 is not responsible for the signaling events that distinguish between H-Ras and N-Ras. Importantly, this work identifies the C-terminal HVR, especially the flexible linker domain with two consecutive proline residues Pro173 and Pro174, as a critical domain that contributes to activation of H-Ras and its invasive potential in human breast epithelial cells. The present study sheds light on the structural basis for the Ras isoform-specific invasive program of breast epithelial cells, providing information for the development of agents that specifically target invasion-related H-Ras pathways in human cancer.

Structure and function of cytochromes P450 2B: from mechanism-based inactivators to X-ray crystal structures and back

This article reviews work from the author dating back to 1978 and focuses on the structural basis of cytochrome P450 (P450) function using available contemporary techniques. Early studies used mechanism-based inactivators that bound to the protein moiety of hepatic P450s to try to localize the active site. Subsequent studies used cDNA cloning, heterologous expression, site-directed mutagenesis, and homology modeling based on multiple bacterial P450 X-ray crystal structures to predict the active sites of CYP2B enzymes with considerable accuracy. Breakthroughs in engineering and expression of mammalian P450s enabled us to determine our first X-ray crystal structure of ligand-free rabbit CYP2B4. To date, we have solved 11 CYP2B4 and three human CYP2B6 structures, which represent four significantly different conformations. The plasticity of CYP2B4 has been confirmed by deuterium exchange mass spectrometry and is substantiated by molecular dynamics simulations. In addition to major movement of secondary structure elements, more subtle reorientation of active site side chains, especially Phe206, Phe297, and Glu301, contributes to the ability of CYP2B enzymes to bind various ligands. Isothermal titration calorimetry has proven to be a useful tool for studying the thermodynamics of ligand binding to CYP2B4 and CYP2B6, and NMR has enabled study of ligand binding orientation in solution as an adjunct to X-ray crystallography. A major challenge remains to harness the power of the various approaches to facilitate prediction of CYP2B specificity and inhibition.

Structural characterization of the complex between alpha-naphthoflavone and human cytochrome P450 1B1

The atomic structure of human P450 1B1 was determined by x-ray crystallography to 2.7 Å resolution with α-naphthoflavone (ANF) bound in the active site cavity. Although the amino acid sequences of human P450s 1B1 and 1A2 have diverged significantly, both enzymes exhibit narrow active site cavities, which underlie similarities in their substrate profiles. Helix I residues adopt a relatively flat conformation in both enzymes, and a characteristic distortion of helix F places Phe(231) in 1B1 and Phe(226) in 1A2 in similar positions for π-π stacking with ANF. ANF binds in a distinctly different orientation in P450 1B1 from that observed for 1A2. This reflects, in part, divergent conformations of the helix B'-C loop that are stabilized by different hydrogen-bonding interactions in the two enzymes. Additionally, differences between the two enzymes for other amino acids that line the edges of the cavity contribute to distinct orientations of ANF in the two active sites. Thus, the narrow cavity is conserved in both P450 subfamily 1A and P450 subfamily 1B with sequence divergence around the edges of the cavity that modify substrate and inhibitor binding. The conservation of these P450 1B1 active site amino acid residues across vertebrate species suggests that these structural features are conserved.

Brief historical overview and recent progress on cytochromes P450: adaptation of aerobic organisms to their chemical environment and new mechanisms of prodrug bioactivation

The present brief overview of the history of the development of our knowledge on cytochromes P450 (P450s) illustrates the spectacular progress that have been made on P450 mechanisms and structures especially during these last 20 years. Recently published structures of mammalian P450-substrate complexes have shown the great diversity of size, shape, and binding modes that are offered by the conformationally flexible substrate binding sites of xenobiotic-metabolizing P450s. They have also shown that these binding sites can adapt themselves to the great structural diversity of xenobiotics, to facilitate their oxidation and elimination. Our present detailed knowledge of the mechanisms and chemistry of P450s allows us to understand, at the molecular level, the origin of the various consequences of P450-dependent metabolism of drugs in pharmacology and toxicology. This is here illustrated by recent data on the detailed mechanism of bioactivation of the anti-thrombotic prodrugs ticlopidine, clopidogrel, and prasugrel.

Conformational plasticity and structure/function relationships in cytochromes P450

The cytochrome P450s are a superfamily of enzymes that are found in all kingdoms of living organisms, and typically catalyze the oxidative addition of atomic oxygen to an unactivated C-C or C-H bond. Over 8000 nonredundant sequences of putative and confirmed P450 enzymes have been identified, but three-dimensional structures have been determined for only a small fraction of these. While all P450 enzymes for which structures have been determined share a common global fold, the flexibility and modularity of structure around the active site account for the ability of P450 enzymes to accommodate a vast number of structurally dissimilar substrates and support a wide range of selective oxidations. In this review, known P450 structures are compared, and some structural criteria for prediction of substrate selectivity and reaction type are suggested. The importance of dynamic processes such as redox-dependent and effector-induced conformational changes in determining catalytic competence and regio- and stereoselectivity is discussed, and noncrystallographic methods for characterizing P450 structures and dynamics, in particular, mass spectrometry and nuclear magnetic resonance spectroscopy are reviewed.

Studies on 16α-Hydroxylation of Steroid Molecules and Regioselective Binding Mode in Homology-Modeled Cytochrome P450-2C11

We investigated the 16α-hydroxylation of steroid molecules and regioselective binding mode in homology-modeled cytochrome P450-2C11 to correlate the biological study with the computational molecular modeling. It revealed that there was a positive relationship between the observed inhibitory potencies and the binding free energies. Docking of steroid molecules into this homology-modeled CYP2C11 indicated that 16α-hydroxylation is favored with steroidal molecules possessing the following components, (1) a bent A-B ring configuration (5β-reduced), (2) C-3 α-hydroxyl group, (3) C-17β-acetyl group, and (4) methyl group at both the C-18 and C-19. These respective steroid components requirements were defined as the inhibitory contribution factor. Overall studies of the male rat CYP2C11 metabolism revealed that the above-mentioned steroid components requirements were essential to induce an effective inhibition of [³H]progesterone 16α-hydroxylation. As far as docking of homology-modeled CYP2C11 against investigated steroids is concerned, they are docked at the active site superimposed with flurbiprofen. It was also found that the distance between heme iron and C16α-H was between 4 to 6 Å and that the related angle was in the range of 180 ± 45°.

Using a homology model of cytochrome P450 2D6 to predict substrate site of metabolism

CYP2D6 is an important enzyme that is involved in first pass metabolism and is responsible for metabolizing ~25% of currently marketed drugs. A homology model of CYP2D6 was built using X-ray structures of ligand-bound CYP2C5 complexes as templates. This homology model was used in docking studies to rationalize and predict the site of metabolism of known CYP2D6 substrates. While the homology model was generally found to be in good agreement with the recently solved apo (ligand-free) X-ray structure of CYP2D6, significant differences between the structures were observed in the B' and F-G helical region. These structural differences are similar to those observed between ligand-free and ligand-bound structures of other CYPs and suggest that these conformational changes result from induced-fit adaptations upon ligand binding. By docking to the homology model using Glide, it was possible to identify the correct site of metabolism for a set of 16 CYP2D6 substrates 85% of the time when the 5 top scoring poses were examined. On the other hand, docking to the apo CYP2D6 X-ray structure led to a loss in accuracy in predicting the sites of metabolism for many of the CYP2D6 substrates considered in this study. These results demonstrate the importance of describing substrate-induced conformational changes that occur upon binding. The best results were obtained using Glide SP with van der Waals scaling set to 0.8 for both the receptor and ligand atoms. A discussion of putative binding modes that explain the distribution of metabolic sites for substrates, as well as a relationship between the number of metabolic sites and substrate size, are also presented. In addition, analysis of these binding modes enabled us to rationalize the typical hydroxylation and O-demethylation reactions catalyzed by CYP2D6 as well as the less common N-dealkylation.

A single amino acid residue, Ala 105, confers 16alpha-hydroxylase activity to human cytochrome P450 17alpha-hydroxylase/17,20 lyase

In adrenal steroidogenesis, CYP17 catalyses the 17alpha-hydroxylation of pregnenolone and progesterone and the subsequent 17,20-lyase reaction, yielding adrenal androgens. The enzyme exhibits distinctly different selectivities towards these substrates in various species. CYP17 has also been shown to exhibit 16alpha-hydroxylase activity towards progesterone in some species, with only human and chimp CYP17 catalysing the biosynthesis of substantial amounts of 16-OHprogesterone. The 16alpha-hydroxylase activity was investigated by introducing an Ala105Leu substitution into human CYP17. The converse mutation, Leu105Ala was introduced into the baboon, goat and pig enzymes. Wt human CYP17 converted approximately 30% progesterone to 16-OHprogesterone while the Ala105Leu mutant converted negligible amounts to 16-OHprogesterone ( approximately 9%), comparable to wt CYP17 of the other three species when expressed in COS-1 cells. The ratio of 17-hydroxylated products to 16-OHprogesterone of human CYP17 was 2.7 and that of the mutant human construct 10.5. Similar ratios were observed for human and goat CYP17 with the corresponding Ala or Leu residues. Although the Leu105Ala mutation of both baboon and pig CYP17 exhibited the same trend regarding the ratios, the rate of progesterone conversion was reduced. Coexpression with cytochrome b(5) significantly decreased the ratio of 17-hydroxylated products to 16-OHprogesterone in the Leu105 constructs, while effects were negligible with Ala at this position. Homology models show that Ala105 faces towards the active pocket in the predicted B'-C domain of CYP17. The smaller residue allows more flexibility of movement in the active pocket than Leu, presenting both the C16 and C17 of progesterone to the iron-oxy complex.

Two-dimensional NMR and all-atom molecular dynamics of cytochrome P450 CYP119 reveal hidden conformational substates

Cytochrome P450 enzymes are versatile catalysts involved in a wide variety of biological processes from hormonal regulation and antibiotic synthesis to drug metabolism. A hallmark of their versatility is their promiscuous nature, allowing them to recognize a wide variety of chemically diverse substrates. However, the molecular details of this promiscuity have remained elusive. Here, we have utilized two-dimensional heteronuclear single quantum coherence NMR spectroscopy to examine a series of mutants site-specific labeled with the unnatural amino acid, [(13)C]p-methoxyphenylalanine, in conjunction with all-atom molecular dynamics simulations to examine substrate and inhibitor binding to CYP119, a P450 from Sulfolobus acidocaldarius. The results suggest that tight binding hydrophobic ligands tend to lock the enzyme into a single conformational substate, whereas weak binding low affinity ligands bind loosely in the active site, resulting in a distribution of localized conformers. Furthermore, the molecular dynamics simulations suggest that the ligand-free enzyme samples ligand-bound conformations of the enzyme and, therefore, that ligand binding may proceed largely through a process of conformational selection rather than induced fit.