Structure prediction and R115866 binding study of human CYP26A1: homology modelling,fold recognition,molecular docking and MD simulations

Two three-dimensional (3D) models of human cytochrome P450 26A1 (CYP26A1) were constructed using the programs Modeller and Sybyl-GeneFold, respectively. After refinement by molecular mechanics and molecular dynamics (MD) simulations, the two models were validated by structure analysis-validation online server. Subsequently, a flexible docking study was performed on the model constructed by GeneFold with the potent and specific inhibitor R115866 to examine the enzyme–inhibitor interactions. From the docking results, we can see R115866 interacts with amino acid residues at the active site by multiple hydrophobic interactions including the side chains of His111, Trp112, Ser115, Val116, Leu125, Ser126, Leu221, Phe222, Glu296, Phe299, Gly300, Glu303, Thr304, Pro371 and the cofactor heme. Trp112 and Thr304 form hydrogen bonds with R115866 and play important roles in stabilising the complex. This constructed CYP26A1 model may provide an opportunity to understand the action mode of the enzyme and could be useful in designing novel retinoic acid metabolism blocking agents (RAMBAs).     

Homology model of 1alpha,25-dihydroxyvitamin D3 24-hydroxylase cytochrome P450 24A1 (CYP24A1): active site architecture and ligand binding

Homology models of cytochrome P450 24A1 (CYP24A1) were constructed using three human P450 structures, CYP2C8, CYP2C9 and CYP3A4 as templates for the model building. Using molecular operating environment (MOE) software the lowest energy CYP24A1 model was then assessed for stereochemical quality and side chain environment. Further active site optimisation of the CYP24A1 model built using the CYP3A4 template was performed by molecular dynamics to generate a final CYP24A1 model. The natural substrate, 1,25-dihydroxyvitamin D(3) (calcitriol) and the CYP24 inhibitor (R)-N-(2-(1H-imidazol-1-yl)-2-phenylethyl)-4'-chlorobiphenyl-4-carboxamide ((R)-VID-400) were docked into the model allowing further validation of the active site architecture. Using the docking studies structurally and functionally important residues were identified with subsequent characterisation of secondary structure.     

A novel human cytochrome P450, CYP26C1, involved in metabolism of 9-cis and all-trans isomers of retinoic acid

Retinoids are potent regulators of cell proliferation, cell differentiation, and morphogenesis and are important therapeutic agents in oncology and dermatology. The gene regulatory activity of endogenous retinoids is effected primarily by retinoic acid isomers (all-trans and 9-cis) that are synthesized from retinaldehyde precursors in a broad range of tissues and act as ligands for nuclear retinoic acid receptors. The catabolism of all-trans-retinoic acid (atRA) is an important mechanism of controlling RA levels in cell and tissues. We have previously identified two cytochrome P450s, P450RAI-1 and P450RAI-2 (herein named CYP26A1 and CYP26B1), which were shown to be responsible for catabolism of atRA both in the embryo and the adult. In this report, we describe the identification, molecular cloning, and substrate characterization of a third member of the CYP26 family, named CYP26C1. Transiently transfected cells expressing CYP26C1 convert atRA to polar water-soluble metabolites similar to those generated by CYP26A1 and -B1. Competition studies with all-trans, 13-cis, and 9-cis isomers of retinoic acid demonstrated that atRA was the preferred substrate for CYP26C1. Although CYP26C1 shares extensive sequence similarity with CYP26A1 and CYP26B1, its catalytic activity appears distinct from those of other CYP26 family members. Specifically, CYP26C1 can also recognize and metabolize 9-cis-RA and is much less sensitive than the other CYP26 family members to the inhibitory effects of ketoconazole. CYP26C1 is not widely expressed in the adult but is inducible by RA in HPK1a, transformed human keratinocyte cell lines. This third CYP26 member may play a specific role in catabolizing both all-trans and 9-cis isomers of RA.     

Arabidopsis CYP72C1 is an atypical cytochrome P450 that inactivates brassinosteroids

Cytochrome P450 monooxygenases (P450s) are a diverse family of proteins that have specialized roles in secondary metabolism and in normal cell development. Two P450s in particular, CYP734A1 and CYP72C1, have been identified as brassinosteroid-inactivating enzymes important for steroid-mediated signal transduction in Arabidopsis thaliana. Genetic analyses have demonstrated that these P450s modulate growth throughout plant development. While members of the CYP734A subfamily inactivate brassinosteroids through C-26 hydroxylation, the biochemical activity of CYP72C1 is unknown. Because CYP734A1 and CYP72C1 in Arabidopsis diverge more than brassinosteroid inactivating P450s in other plants, this study examines the structure and biochemistry of each enzyme. Three-dimensional models were generated to examine the substrate binding site structures and determine how they might affect the function of each P450. These models have indicated that the active site of CYP72C1 does not contain several conserved amino acids typically needed for substrate hydroxylation. Heterologous expression of these P450s followed by substrate binding analyses have indicated that CYP734A1 binds active brassinosteroids, brassinolide and castasterone, as well as their upstream precursors whereas CYP72C1 binds precursors more effectively. Seedling growth assays have demonstrated that the genetic state of CYP734A1, but not CYP72C1, affected responsiveness to high levels of exogenous brassinolide supporting our observations that CYP72C1 acts on brassinolide precursors. Although there may be some overlap in their physiological function, the distinct biochemical functions of these proteins in Arabidopsis has significant potential to fine-tune the levels of different brassinosteroid hormones throughout plant growth and development.     

Distinct binding of cholesterol and 5beta-cholestane-3alpha,7alpha,12alpha-triol to cytochrome P450 27A1: evidence from modeling and site-directed mutagenesis studies

Cytochrome P450 27A1 (P450 27A1 or CYP27A1) is an important enzyme that participates in different pathways of cholesterol degradation as well as in the activation of vitamin D(3). Several approaches were utilized to investigate how two physiological substrates, cholesterol and 5beta-cholestane-3alpha,7alpha,12alpha-triol, interact with CYP27A1. The enzyme active site was first probed spectrally by assessing binding of the two substrates and five substrate analogues followed by computer modeling and site-directed mutagenesis. The computer models suggest that the spatial positions and orientations of cholesterol and 5beta-cholestane-3alpha,7alpha,12alpha-triol are different in the enzyme active site. As a result, some of the active site residues interact with both substrates, although they are situated differently relative to each steroid, and some residues bind only one substrate. Mutation of the overlapping substrate-contact residues (W100, H103, T110, M301C, V367, I481, and V482) affected CYP27A1 binding and enzyme activity in a substrate-dependent manner and allowed identification of several important side chains. T110 is proposed to interact with the 12alpha-hydroxyl of 5beta-cholestane-3alpha,7alpha,12alpha-triol, whereas V367 seems to be crucial for correct positioning of the cholesterol C26 methyl group and for regioselective hydroxylation of this substrate. Distinct binding of the CYP27A1 substrates may provide insight into why phenotypic manifestations of cerebrotendinous xanthomatosis, a disease associated with CYP27A1 deficiency, are so diverse.     

A role for the strained phenylalanine ring rotation induced by substrate binding to cytochrome CYP102A1

X-ray crystal structures of CYP102A1 (P450BM-3) have shown that PHE87 rotates upon substrate binding and is in contact with the heme cofactor. Analysis of substrate binding data combined with DFT calculations suggest that the ring is rotated into an unfavorable interaction with the heme and this could drive active site rearrangement.     

A second CYP26 P450 in humans and zebrafish: CYP26B1.

A new P450 gene has been found in humans. It has 44% sequence identity to CYP26A1 from human and mouse, which places it in a new subfamily, CYP26B. There is only one human EST from a retinal library (AA012833) that matches the coding region. No homologous ESTs are found in mouse. A zebrafish EST AI721901 shows 68% identity to the human protein. This zebrafish EST is only 41% identical to the zebrafish CYP26A1 protein sequence, so it represents the homolog of the human CYP26B1 sequence. It is not known if this gene product will act on all-trans-retinoic acid like the CYP26A1 protein or if it might hydroxylate the 9-cis- or 13-cis-retinoic acid isoforms not recognized by CYP26A1. The importance of the CYP26A1 P450 in mouse and zebrafish development flags the CYP26B1 gene as a potential developmental gene.     

Binding of a Cyano- and Fluoro-containing Drug Bicalutamide to Cytochrome P450 46A1: UNUSUAL FEATURES AND SPECTRAL RESPONSE*

Cytochrome P450 46A1 (CYP46A1) is the cholesterol 24-hydroxylase initiating the major pathways of cholesterol removal from the brain, and bicalutamide (BIC) is a drug of choice for the treatment of progressive androgen-dependent prostate cancer. We evaluated the interactions of BIC with CYP46A1 by x-ray crystallography and by conducting solution and mutagenesis studies. Because BIC is administered to patients as a racemic mixture of the S and R isomers, we studied all three, racemic BIC as well as the S and R isomers. Co-crystallization of CYP46A1 with racemic BIC led to structure determinations at 2.1 Å resolution with the drug complexes exhibiting novel properties. Both BIC isomers bind to the CYP46A1 active site in very similar single orientation and adopt an energetically unfavorable folded conformation. This folded BIC conformation is correlated with drug-induced localized shifts of amino acid side chains in CYP46A1 and unusual interactions in the CYP46A1-BIC complex. One of these interactions involves a water molecule that is coordinated to the P450 heme iron and also hydrogen-bonded to the BIC nitrile. Due to polarization of the water in this environment, the heme elicits previously unreported types of P450 spectral responses. We also observed that access to the P450 active site was affected by differential recognition of S versus R isomers at the CYP46A1 surface arising from BIC conformational polymorphism.     

Structural Analysis of the Streptomyces avermitilis CYP107W1-Oligomycin A Complex and Role of the Tryptophan 178 Residue

CYP107W1 from Streptomyces avermitilis is a cytochrome P450 enzyme involved in the biosynthesis of macrolide oligomycin A. A previous study reported that CYP107W1 regioselectively hydroxylated C12 of oligomycin C to produce oligomycin A, and the crystal structure of ligand free CYP107W1 was determined. Here, we analyzed the structural properties of the CYP107W1-oligomycin A complex and characterized the functional role of the Trp178 residue in CYP107W1. The crystal structure of the CYP107W1 complex with oligomycin A was determined at a resolution of 2.6 Å. Oligomycin A is bound in the substrate access channel on the upper side of the prosthetic heme mainly by hydrophobic interactions. In particular, the Trp178 residue in the active site intercalates into the large macrolide ring, thereby guiding the substrate into the correct binding orientation for a productive P450 reaction. A Trp178 to Gly mutation resulted in the distortion of binding titration spectra with oligomycin A, whereas binding spectra with azoles were not affected. The Gly178 mutant’s catalytic turnover number for the 12-hydroxylation reaction of oligomycin C was highly reduced. These results indicate that Trp178, located in the open pocket of the active site, may be a critical residue for the productive binding conformation of large macrolide substrates.     

Inhibition of human cytochrome P450 1A2 by flavones: A molecular modeling study

Cytochrome P450 1A2 metabolizes a number of important drugs, procarcinogens, and endogenous compounds. Several flavones, a class of phytochemicals consumed in the human diet, have been shown to differentially inhibit human P450 1A2-mediated methoxyresorufin demethylase. A molecular model of this P450 was constructed in order to elucidate the molecular basis of the P450-flavone interaction. Flavone and its 3,5,7-trihydroxy and 3,5,7-trimethoxy derivatives were docked into the active site to assess their mode of binding. The site is hydrophobic and includes several residues that hydrogen bond with substituents on the flavone nucleus. The binding interactions of these flavones in the modeled active side are consistent with their relative inhibitory potentials, namely 3,5,7-trihydroxylflavone > flavone >3,5,7-trimethoxylflavone, toward P450 1A2-mediated methoxyresorufin demethylation.     

Inhibition of human cytochrome P450 1A2 by flavones: A molecular modeling study

Cytochrome P450 1A2 metabolizes a number of important drugs, procarcinogens, and endogenous compounds. Several flavones, a class of phytochemicals consumed in the human diet, have been shown to differentially inhibit human P450 1A2-mediated methoxyresorufin demethylase. A molecular model of this P450 was constructed in order to elucidate the molecular basis of the P450-flavone interaction. Flavone and its 3,5,7-trihydroxy and 3,5,7-trimethoxy derivatives were docked into the active site to assess their mode of binding. The site is hydrophobic and includes several residues that hydrogen bond with substituents on the flavone nucleus. The binding interactions of these flavones in the modeled active side are consistent with their relative inhibitory potentials, namely 3,5,7-trihydroxylflavone > flavone >3,5,7-trimethoxylflavone, toward P450 1A2-mediated methoxyresorufin demethylation.     

Computer modeling of cytochrome P450 2E1 three-dimensional structure

A computer model of human cytochrome P450 2E1 (CYP2E1) three-dimensional structure and active site was constructed based on homology with crystallographic coordinates of CYP2C5 and CYP2C9. A high degree of secondary structure homology for human, mouse, rat and rabbit CYP2E1 was demonstrated. The location of heme and the supporting alpha-helices was established. CYP2E1, CYP2C5 and CYP2C9 active sites are distinguished by pocket size and their amino acid residues composition. Key amino acid residues forming the active site channel and substrate-binding cavity are presented. Active site surface area and volume for CYP2E1, CYP2C5 and CYP2C9 were calculated.     

Crystal structures of substrate-free and retinoic acid-bound cyanobacterial cytochrome P450 CYP120A1

The crystal structures of substrate-free and all-trans-retinoic acid-bound CYP120A1 from Synechocystis sp. PCC 6803 were determined at 2.4 and 2.1 A resolution, respectively, representing the first structural characterization of a cyanobacterial P450. Features of CYP120A1 not observed in other P450 structures include an aromatic ladder flanking the channel leading to the active site and a triple-glycine motif within SRS5. Using spectroscopic methods, CYP120A1 is shown to bind 13-cis-retinoic acid, 9-cis-retinoic acid, and retinal with high affinity and dissociation constants of less than 1 microM. Metabolism of retinoic acid by CYP120A1 suggests that CYP120A1 hydroxylates a variety of retinoid derivatives in vivo. On the basis of the retinoic acid-bound CYP120A1 crystal structure, we propose that either carbon 2 or the methyl groups (C16 or C17) of the beta-ionone ring are modified by CYP120A1.     

Exploring CYP1A1 as anticancer target: homology modeling and in silico inhibitor design

Microsomal cytochrome P450 family 1 enzymes has great importance in the bioactivation of mutagens. P450 catalyzed reactions involve a wide range of substrates, and this versatality is reflected in a structural diversity, evident in the active sites of available P450 structures. This structure offers a template to study further structure-function relationships of alternative substrates and other cytochrome P450 family 1 members. In this paper, we document a homology model of CYP P450 1A1 from Homo sapiens, developed on the basis of template crystal structure of human microsomal P450 1a2 in complex with inhibitor (PDB Id: 2HI4). Homology modeling is performed at the programs, both in the commercial and public realms. We tried to explore CYP1A1 as a potential target for anticancer chemotherapy. To gain an insight into the binding of ligands with enzyme, protein-ligand complex was developed by including information about the known ligand as spatial restraints and optimizing the mutual interactions between the ligand and the binding site. Active site characterization and the study for involvement of specific aminoacids in binding with ligand, facilitates structure based inhibitor design. This study should prove useful in the design and development of potential novel anticancer agents. Figure The refined structure of homology model of CYP1A1     

Structural basis for three-step sequential catalysis by the cholesterol side chain cleavage enzyme CYP11A1

Mitochondrial cytochrome P450 11A1 (CYP11A1 or P450 11A1) is the only known enzyme that cleaves the side chain of cholesterol, yielding pregnenolone, the precursor of all steroid hormones. Pregnenolone is formed via three sequential monooxygenation reactions that involve the progressive production of 22R-hydroxycholesterol (22HC) and 20α,22R-dihydroxycholesterol, followed by the cleavage of the C20-C22 bond. Herein, we present the 2.5-Å crystal structure of CYP11A1 in complex with the first reaction intermediate, 22HC. The active site cavity in CYP11A1 represents a long curved tube that extends from the protein surface to the heme group, the site of catalysis. 22HC occupies two-thirds of the cavity with the 22R-hydroxyl group nearest the heme, 2.56 Å from the iron. The space at the entrance to the active site is not taken up by 22HC but filled with ordered water molecules. The network formed by these water molecules allows the “soft” recognition of the 22HC 3β-hydroxyl. Such a mode of 22HC binding suggests shuttling of the sterol intermediates between the active site entrance and the heme group during the three-step reaction. Translational freedom of 22HC and torsional motion of its aliphatic tail are supported by solution studies. The CYP11A1–22HC co-complex also provides insight into the structural basis of the strict substrate specificity and high catalytic efficiency of the enzyme and highlights conserved structural motifs involved in redox partner interactions by mitochondrial P450s.     

2 D - QSAR studies on CYP26A1 inhibitory activity of 1-[benzofuran-2-yl-(4-alkyl/aryl-phenyl)-methyl]- 1 H-triazoles

The Quantitative Structure Activity Relationship (QSAR) study is performed over a set of 15, 4-alkyl/aryl-substituted 1- [benzofuran-2-yl-phenylmethyl]-1 H-triazoles derivatives. This study is based on the application of physicochemical parameters in QSAR. The parameters include (MR (molar refractivity), MW (molecular weight), Pc (parachor), St (surface tension), D (density), Ir (index of refraction) and log P (partition coefficient). The parameters describing physiochemical properties are used as independent variables and the biological activity (IC(50)) is considered as dependent variable in multiple regression analysis. Different models were generated with high co-efficient of determination (R(2)). The 2D-QSAR study identified compounds capable of inhibiting the metabolic breakdown of the retinoid (trans-retinoic acid (ATRA)) involved in the activation of specific nuclear Retinoic acid receptors (RARs). This study identifies R115866 as a potential inhibitor of the cytochrome P450 (CYP) mediated metabolism with increased RA levels for retinoid actions.     

Differences in flexibility of active sites of cytochromes P450 probed by resonance Raman and UV-Vis absorption spectroscopy.

Spectroscopic methods reveal differences in flexibility and stability of P450 forms. Among microsomal P450s, the most flexible active site has been found in the CYP3A4 enzyme as it is compressible and the heme vinyl side chains may adopt two different conformations. On the other hand, active site of this enzyme denatures quite easily upon hydrostatic pressure. The most rigid active site able to withstand the effect of high pressure has CYP1A2. The bacterial CYP102 (BM3) flavocytochrome has also a rather stable, but flexible active site. The differences between CYP3A4 and CYP1A2 active sites apparently reflect their ability to bind various substrates: whereas the CYP3A4 binds a vast variety of structures, the CYP1A2 preferentially binds planar, aromatic structures and its substrate specificity is relatively narrow.     

Retinoic acid and Cyp26b1 are critical regulators of osteogenesis in the axial skeleton

Retinoic acid (RA) plays important roles in diverse biological processes ranging from germ cell specification to limb patterning. RA ultimately exerts its effect in the nucleus, but how RA levels are being generated and maintained locally is less clear. Here, we have analyzed the zebrafish stocksteif mutant, which exhibits severe over-ossification of the entire vertebral column. stocksteif encodes cyp26b1, a cytochrome P450 member that metabolizes RA. The mutant is completely phenocopied by treating 4 dpf wild-type embryos with either RA or the pharmacological Cyp26 blocker R115866, thus identifying a previously unappreciated role for RA and cyp26b1 in osteogenesis of the vertebral column. Cyp26b1 is expressed within osteoblast cells, demonstrating that RA levels within these cells need to be tightly controlled. Furthermore, we have examined the effect of RA on osteoblasts in vivo. As numbers of osteoblasts do not change upon RA treatment, we suggest that RA causes increased activity of axial osteoblasts, ultimately resulting in defective skeletogenesis.     

Molecular modeling of cytochrome P450 1A1: enzyme-substrate interactions and substrate binding affinities

Human cytochrome P450 1A1, which is present in lungs, plays an important role in the metabolic activation of chemical carcinogens, and in particular, is thought to be linked to lung cancer. The mechanism of carcinogenesis is related to the enzyme's ability to oxidize highly toxic compounds, such as polycyclic aromatic hydrocarbons (PAHs), to their carcinogenic derivatives. In order to better understand P450 1A1 function, a homology model of this enzyme has been constructed. The model has been based on the structure of P450 2C5, the first mammalian P450 to be crystallized. The coordinates of the model have been calculated using a consensus strategy, and the resulting structure has been evaluated with the ProStat and Profiles-3D programs. P450 1A1 substrates, such as benzo[a]pyrene, ethoxyresorufin and methoxyresorufin, were then docked into the active site of the model, and key amino acid residues able to interact with the substrate, have been identified. The analysis of enzyme-substrate interactions indicated that hydrophobic interactions are mainly responsible for binding of these substrates in the active site. Moreover, the non-bond enzyme-substrate interaction energy for ethoxyresorufin was lower than that for methoxyresorufin, which is consistent with higher activity of 1A1 towards the former substrate. Key residue Val-382 may play an important role in these interactions. Additionally, we performed binding free energy calculations for the three substrates. The obtained values were similar to those observed experimentally, which suggests that this approach might be useful for prediction of binding constants.