Crystal structure of the human prostacyclin synthase

In silico design and in vitro assessment of anti-Helicobacter pylori compounds as potential small-molecule arginase inhibitors

Related to a variety of gastrointestinal disorders ranging from gastric ulcer to gastric adenocarcinoma, the infection caused by the gram-negative bacteria Helicobacter pylori (H. pylori) poses as a great threat to human health; hence, the search for new treatments is a global priority. The H. pylori arginase (HPA) protein has been widely studied as one of the main virulence factors of this bacterium, being involved in the prevention of nitric oxide-mediated bacterial cell death, which is a central component of innate immunity. Given the growing need for the development of new drugs capable of combating the infection by H. pylori, the present work describes the search for new HPA inhibitors, using virtual screening techniques based on molecular docking followed by the evaluation of the proposed modes of interaction at the HPA active site. In vitro studies of minimum inhibitory concentration (MIC) and minimum bactericidal concentration (MBC), followed by cytotoxicity activity in gastric adenocarcinoma and non-cancer cells, were performed. The results highlighted compounds 6, 11, and 13 as potential inhibitors of HPA; within these compounds, the results indicated 13 presented an improved activity toward H. pylori killing, with MIC and MBC both at 64 µg/mL. Moreover, compound 13 also presented a selectivity index of 8.3, thus being more selective for gastric adenocarcinoma cells compared to the commercial drug cisplatin. Overall, the present work demonstrates the search strategy based on in silico and in vitro techniques is able to support the rational design of new anti-H. pylori drugs.

Human cytochrome P450 enzymes 5-51 as targets of drugs and natural and environmental compounds: mechanisms, induction, and inhibition - toxic effects and benefits

Cytochrome P450 (P450, CYP) enzymes have long been of interest due to their roles in the metabolism of drugs, pesticides, pro-carcinogens, and other xenobiotic chemicals. They have also been of interest due to their very critical roles in the biosynthesis and metabolism of steroids, vitamins, and certain eicosanoids. This review covers the 22 (of the total of 57) human P450s in Families 5-51 and their substrate selectivity. Furthermore, included is information and references regarding inducibility, inhibition, and (in some cases) stimulation by chemicals. We update and discuss important aspects of each of these 22 P450s and questions that remain open.

A combined computational strategy of sequence and structural analysis predicts the existence of a functional eicosanoid pathway in Drosophila melanogaster

This study reports on a putative eicosanoid biosynthesis pathway in Drosophila melanogaster and challenges the currently held view that mechanistic routes to synthesize eicosanoid or eicosanoid-like biolipids do not exist in insects, since to date, putative fly homologs of most mammalian enzymes have not been identified. Here we use systematic and comprehensive bioinformatics approaches to identify most of the mammalian eicosanoid synthesis enzymes. Sensitive sequence analysis techniques identified candidate Drosophila enzymes that share low global sequence identities with their human counterparts. Twenty Drosophila candidates were selected based upon (a) sequence identity with human enzymes of the cyclooxygenase and lipoxygenase branches, (b) similar domain architecture and structural conservation of the catalytic domain, and (c) presence of potentially equivalent functional residues. Evaluation of full-length structural models for these 20 top-scoring Drosophila candidates revealed a surprising degree of conservation in their overall folds and potential analogs for functional residues in all 20 enzymes. Although we were unable to identify any suitable candidate for lipoxygenase enzymes, we report structural homology models of three fly cyclooxygenases. Our findings predict that the D. melanogaster genome likely codes for one or more pathways for eicosanoid or eicosanoid-like biolipid synthesis. Our study suggests that classical and/or novel eicosanoids mediators must regulate biological functions in insects–predictions that can be tested with the power of Drosophila genetics. Such experimental analysis of eicosanoid biology in a simple model organism will have high relevance to human development and health.

Molecular characterization of Penicillium oxalicum 6R,8R-linoleate diol synthase with new regiospecificity

Diol synthase-derived metabolites are involved in the sexual and asexual life cycles of fungi. A putative diol synthase from Penicillium oxalicum was found to convert palmitoleic acid (16:1n-7), oleic acid (18:1n-9), linoleic acid (18:2n-6), and α-linolenic acid (18:3n-3) to 6S,8R-dihydroxy-9(Z)-hexadecenoic acid, 6R,8R-dihydroxy-9(Z)-octadecenoic acid, 6R,8R-dihydroxy-9,12(Z,Z)-octadecadienoic acid, and 6S,8R-dihydroxy-9,12,15(Z,Z,Z)-octadecatrienoic acid, respectively, which were identified by liquid chromatography-tandem mass spectrometry (LC-MS/MS) and nuclear magnetic resonance (NMR) spectroscopy analyses. The specific activity and catalytic efficiency of P. oxalicum 6,8-diol synthase were the highest for 18:2n-6, indicating that the enzyme is a 6R,8R-linoleate diol synthase (6R,8R-LDS) with new regiospecificity. This is the first report of a 6R,8R-LDS. LDS is a fusion protein consisting of a dioxygenase domain at the N-terminus and a cytochrome P450/hydroperoxide isomerase (P450/HPI) domain at the C-terminus. The putative active-site residues in the C-terminal domain of P. oxalicum 6R,8R-LDS were proposed based on a substrate-docking homology model. The results of the site-directed mutagenesis within C-terminal P450 domain suggested that Asn⁸⁸⁶, Arg⁷⁰⁷, and Arg⁹³⁴, are catalytic importance and belong to the catalytic groove. Phe⁷⁹⁴ and Gln⁸⁸⁹ were found to be involved in the regiospecific rearrangement of hydroperoxide, while the F794E and Q889A variants of P. oxalicum 6,8-LDS acted as 7,8- and 8,11-LDSs, respectively. All these mutations critically affected the HPI activity of P. oxalicum 6R,8R-LDS.

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.

Prostaglandin terminal synthases as novel therapeutic targets

Non-steroidal anti-inflammatory drugs (NSAIDs) exert their anti-inflammatory and anti-tumor effects by reducing prostaglandin (PG) production via the inhibition of cyclooxygenase (COX). However, the gastrointestinal, renal and cardiovascular side effects associated with the pharmacological inhibition of the COX enzymes have focused renewed attention onto other potential targets for NSAIDs. PGH₂, a COX metabolite, is converted to each PG species by species-specific PG terminal synthases. Because of their potential for more selective modulation of PG production, PG terminal synthases are now being investigated as a novel target for NSAIDs. In this review, I summarize the current understanding of PG terminal synthases, with a focus on microsomal PGE synthase-1 (mPGES-1) and PGI synthase (PGIS). mPGES-1 and PGIS cooperatively exacerbate inflammatory reactions but have opposing effects on carcinogenesis. mPGES-1 and PGIS are expected to be attractive alternatives to COX as therapeutic targets for several diseases, including inflammatory diseases and cancer.

Prostacyclin-producing human mesenchymal cells target H19 lncRNA to augment endogenous progenitor function in hindlimb ischaemia

Promoting the paracrine effects of human mesenchymal stem cell (hMSC) therapy may contribute to improvements in patient outcomes. Here we develop an innovative strategy to enhance the paracrine effects of hMSCs. In a mouse hindlimb ischaemia model, we examine the effects of hMSCs in which a novel triple-catalytic enzyme is introduced to stably produce prostacyclin (PGI₂-hMSCs). We show that PGI₂-hMSCs facilitate perfusion recovery and enhance running capability as compared with control hMSCs or iloprost (a stable PGI₂ analogue). Transplanted PGI₂-hMSCs do not incorporate long term into host tissue, but rather they mediate host regeneration and muscle mass gain in a paracrine manner. Mechanistically, this involves long noncoding RNA H19 in promoting PGI₂-hMSC-associated survival and proliferation of host progenitor cells under hypoxic conditions. Together, our data reveal the novel ability of PGI₂-hMSCs to stimulate host regenerative processes and improve physical function by regulating long noncoding RNA in resident progenitor cells.

Human mesenchymal stem cells (hMSC) expressing paracrine factors may enhance therapeutic benefits when transplanted. Here, the authors show that hMSCs stably expressing prostacyclin enhance host regeneration and muscle mass gain in a mouse hindlimb ischaemia model, mediated by the long noncoding RNA H19.

Comparison of intrinsic dynamics of cytochrome p450 proteins using normal mode analysis

Cytochrome P450 enzymes are hemeproteins that catalyze the monooxygenation of a wide-range of structurally diverse substrates of endogenous and exogenous origin. These heme monooxygenases receive electrons from NADH/NADPH via electron transfer proteins. The cytochrome P450 enzymes, which constitute a diverse superfamily of more than 8,700 proteins, share a common tertiary fold but < 25% sequence identity. Based on their electron transfer protein partner, cytochrome P450 proteins are classified into six broad classes. Traditional methods of pro are based on the canonical paradigm that attributes proteins' function to their three-dimensional structure, which is determined by their primary structure that is the amino acid sequence. It is increasingly recognized that protein dynamics play an important role in molecular recognition and catalytic activity. As the mobility of a protein is an intrinsic property that is encrypted in its primary structure, we examined if different classes of cytochrome P450 enzymes display any unique patterns of intrinsic mobility. Normal mode analysis was performed to characterize the intrinsic dynamics of five classes of cytochrome P450 proteins. The present study revealed that cytochrome P450 enzymes share a strong dynamic similarity (root mean squared inner product > 55% and Bhattacharyya coefficient > 80%), despite the low sequence identity (< 25%) and sequence similarity (< 50%) across the cytochrome P450 superfamily. Noticeable differences in Cα atom fluctuations of structural elements responsible for substrate binding were noticed. These differences in residue fluctuations might be crucial for substrate selectivity in these enzymes.

Functional role of the conserved i-helix residue I346 in CYP5A1-Nanodiscs

Thromboxane synthase (CYP5A1) is a non-classical cytochrome P450 (CYP) expressed in human platelets that mediates vascular homeostasis by producing thromboxane A2 (TXA2) through the isomerization of prostaglandin H2 (PGH2). A homology alignment of CYP5A1 with human CYPs indicates that a highly conserved I-helix threonine residue is occupied by an isoleucine at position 346 in CYP5A1. We find that reverse-engineering CYP5A1 to contain either threonine or serine in this position dramatically increases TXA2 formation. Interestingly, the levels of malondialdehyde (MDA), a homolytic fragmentation product of PGH2 formed via a pathway independent of TXA2 formation, remain constant. Furthermore, spectral analysis using two PGH2 substrate analogs supports the observed activity changes in the hydroxyl-containing mutants. The more constrained active site of the I346T mutant displays altered PGH2 substrate analog binding properties. Together these studies provide new mechanistic insights into CYP5A1 mediated isomerization of PGH2 with respect to a critical active site residue.

Analysis of genetic polymorphism and biochemical characterization of a functionally decreased variant in prostacyclin synthase gene (CYP8A1) in humans

Prostacyclin synthase (CYP8A1) is an enzyme responsible for the biosynthesis of prostacyclin (PGI2) which inhibits platelet activation and exhibits anti-inflammatory effect. The objectives of this study were to identify CYP8A1 genetic variants and characterize functional consequences of CYP8A1 variants. In total, 27 variants including four previously unidentified single-nucleotide polymorphisms (SNPs) were identified by direct DNA sequencing in Koreans (n=48). Among them, CYP8A1 A447T and E314Stop were newly assigned as CYP8A1(∗)5 and CYP8A1(∗)6 by the Human Cytochrome P450 Allele Nomenclature Committee, respectively. CYP8A1(∗)5 was found in the heme binding area in three individuals as a heterozygous mutation. To investigate the functional change of CYP8A1(∗)5, CYP8A1(∗)5 and wild-type CYP8A1 protein were overexpressed in an Escherichia coli expression system and purified. Metabolism of PGH2 by the CYP8A1(∗)5 protein exhibited significantly decreased activity, resulting in a 45% decrease in Vmax and a 1.8-fold decrease in intrinsic clearance compared to the wild-type. Based on the predicted crystal structure of CYP8A1(∗)5 using the Molecular Operating Environment platform, the distance from CYP8A1 Cys441 to the heme was altered with a significantly changed binding free energy for the mutant protein. Further studies would be needed to determine the effect of CYP8A1(∗)5 on PGI2 levels in humans.

Structural characterization of human cholesterol 7α-hydroxylase

Hepatic conversion to bile acids is a major elimination route for cholesterol in mammals. CYP7A1 catalyzes the first and rate-limiting step in classic bile acid biosynthesis, converting cholesterol to 7α-hydroxycholesterol. To identify the structural determinants that govern the stereospecific hydroxylation of cholesterol, we solved the crystal structure of CYP7A1 in the ligand-free state. The structure-based mutation T104L in the B' helix, corresponding to the nonpolar residue of CYP7B1, was used to obtain crystals of complexes with cholest-4-en-3-one and with cholesterol oxidation product 7-ketocholesterol (7KCh). The structures reveal a motif of residues that promote cholest-4-en-3-one binding parallel to the heme, thus positioning the C7 atom for hydroxylation. Additional regions of the binding cavity (most distant from the access channel) are involved to accommodate the elongated conformation of the aliphatic side chain. Structural complex with 7KCh shows an active site rigidity and provides an explanation for its inhibitory effect. Based on our previously published data, we proposed a model of cholesterol abstraction from the membrane by CYP7A1 for metabolism. CYP7A1 structural data provide a molecular basis for understanding of the diversity of 7α-hydroxylases, on the one hand, and cholesterol-metabolizing enzymes adapted for their specific activity, on the other hand.

7,8- and 5,8-Linoleate diol synthases support the heterolytic scission of oxygen-oxygen bonds by different amide residues

Linoleate diol synthases (LDS) are fungal dioxygenase-cytochrome P450 fusion enzymes. They oxidize 18:2n-6 sequentially to 8R-hydroperoxylinoleic acid (8R-HPODE) and 7S,8S- or 5S,8R-dihydroxylinoleic acids (DiHODE) by intramolecular oxygen transfer. The P450 domains contain a conserved sequence, Ala-Asn-Gln-Xaa-Gln, presumably located in the I-helices. The Asn938Leu replacement of 7,8-LDS of Gaeumannomyces graminis virtually abolished and the Asn938Asp and Asn938Gln replacements reduced the hydroperoxide isomerase activity. Gln941Leu and Gln941Glu substitutions had little effects. Replacements of the homologous Asn(887) and Gln(890) residues of 5,8-LDS of Aspergillus fumigatus yielded the opposite results. Asn887Leu and Asn887Gln of 5,8-LDS retained 5,8-DiHODE as the main metabolite with an increased formation of 6,8- and 8,11-DiHODE, whereas Gln890Leu almost abolished the 5,8-LDS activity. Replacement of Gln(890) with Glu also retained 5,8-DiHODE as the main product, but shifted oxygenation from C-5 to C-7 and C-11 and to formation of epoxyalcohols by homolytic scission of 8R-HPODE. P450 hydroxylases usually contain an "acid-alcohol" pair in the I-helices for the heterolytic scission of O2 and formation of compound I (Por(+) Fe(IV)=O) and water. The function of the acid-alcohol pair appears to be replaced by two different amide residues, Asn(938) of 7,8-LDS and Gln(890) of 5,8-LDS, for heterolysis of 8R-HPODE to generate compound I.

Conformational changes of enzymes and DNA in molecular dynamics: influenced by pH, temperature, and ligand

Protein conformation, which has been a research hotspot for human diseases, is an important factor of protein properties. Recently, a series of approaches have been utilized to investigate the conformational changes under different conditions. Some of them have gained promising achievements, but it is still deficient in the detail researches at the atomic level. In this chapter, a series of computational examples of protein conformational changes under different pH environment, temperature, and ligand binding are described. We further show some useful methods, such as constant pH molecular dynamics simulations, molecular docking, and molecular mechanics Poisson-Boltzmann surface area/generalized Born surface area calculations. In comparison with the experimental results, the methods mentioned above are reasonable to detect and predict the interaction between residue and residue, residue and DNA, and residue and ligand. Additionally, some crucial interactions that cause protein conformational changes are discovered and discussed in this chapter. In summary, our work can give penetrating information to understand the pH-, temperature-, and ligand-induced conformational change mechanisms.

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.

Expression of fusion proteins of Aspergillus terreus reveals a novel allene oxide synthase

Aspergilli oxidize C18 unsaturated fatty acids by dioxygenase-cytochrome P450 fusion proteins to signal molecules involved in reproduction and host-pathogen interactions. Aspergillus terreus expresses linoleate 9R-dioxygenase (9R-DOX) and allene oxide synthase (AOS) activities in membrane fractions. The genome contains five genes (ATEG), which may code for a 9R-DOX-AOS fusion protein. The genes were cloned and expressed, but none of them oxidized 18:2n-6 to 9R-hydroperoxy-10(E),12(Z)-octadecadienoic acid (9R-HPODE). ATEG_02036 transformed 9R-HPODE to an unstable allene oxide, 9(R),10-epoxy-10,12(Z)-octadecadienoic acid. A substitution in the P450 domain (C1073S) abolished AOS activity. The N964V and N964D mutants both showed markedly reduced AOS activity, suggesting that Asn(964) may facilitate homolytic cleavage of the dioxygen bond of 9R-HPODE with formation of compound II in analogy with plant AOS (CYP74) and prostacyclin synthase (CYP8A1). ATEG_03992 was identified as 5,8-linoleate diol synthase (5,8-LDS). Replacement of Asn(878) in 5,8-LDS with leucine (N878L) mainly shifted ferryl oxygen insertion from C-5 toward C-6, but replacements of Gln(881) markedly affected catalysis. The Q881L mutant virtually abolished the diol synthase activity. Replacement of Gln(881) with Asn, Glu, Asp, or Lys residues augmented the homolytic cleavage of 8R-HPODE with formation of 10-hydroxy-8(9)-epoxy-12(Z)-octadecenoic acid (erythro/threo, 1-4:1) and/or shifted ferryl oxygen insertion from C-5 toward C-11. We conclude that homolysis and heterolysis of the dioxygen bond with formation of compound II in AOS and compound I in 5,8-LDS are influenced by Asn and Gln residues, respectively, of the I-helices. AOS of A. terreus appears to have evolved independently of CYP74 but with an analogous reaction mechanism.

Biochemical pathogenesis of aspirin exacerbated respiratory disease (AERD)

Aspirin exacerbated respiratory disease (AERD) is a distinct clinical entity characterized by eosinophilic rhinosinusitis, asthma and often nasal polyposis. Exposure to aspirin or other nonsteroid anti-inflammatory drugs (NSAIDs) exacerbates bronchospasms with asthma and rhinitis. Disease progression suggests a skewing towards TH2 type cellular response along with moderate to severe eosinophil and mast cell infiltration. Alterations in upper and lower airway cellular milieu with abnormalities in eicosanoid metabolism and altered eicosanoid receptor expression are the key features underlying AERD pathogenesis. Dysregulation of arachidonic acid (AA) metabolism, notably reduced prostaglandin E2 (PGE2) synthesis compared to their aspirin tolerant counterpart and relatively increased PGD2 production, a TH2/eosinophil chemoattractant are reported in AERD. Underproduced PGE2 is metabolized by overexpression of 15 prostaglandin dehydrogenase (15-PGDH) to inactive products further reducing PGE2 at real time. This relives the inhibitory effect of PGE2 on 5-lipoxygenase (5-LOX) resulting in overproduction of cysteinyl leukotrienes (CysLTs). Diminished formation of CysLT antagonists called lipoxins (LXs) also augments CysLTs responsiveness. Occasional intake of NSAIDs favors even more 5-LOX product formation, further narrowing the bronchoconstrictive bottle neck, resulting in acute asthmatic exacerbations along with increased mucus production. This review focuses on abnormalities in biochemical and molecular mechanisms in eicosanoid biosynthesis, eicosanoid receptor dysregulation and associated polymorphisms with special reference to arachidonic acid metabolism in AERD.

CYP2W1, CYP4F11 and CYP8A1 polymorphisms and interaction of CYP2W1 genotypes with risk factors in Mexican women with breast cancer

Breast cancer (BCa) is the leading type of cancer in Mexican women. Genetic factors, such as single nucleotide polymorphisms (SNP) of P450 system, have been reported in BCa. In this report, and for the first time in the literature, we analyzed the rs3735684 (7021 G>A), rs11553651 (15016 G>T) and rs56195291 (60020 C>G) polymorphisms in the CYP2W1, 4F11 and 8A1 genes in patients with BCa and in healthy Mexican women to identify a potential association between these polymorphisms and BCa risk. Patients and controls were used for polymorphism analysis using an allelic discrimination assay with TaqMan probes and confirmed by DNA sequencing. Links with clinic-pathological characteristics were also analyzed. Statistical analysis was performed using the standard χ2 or Fisher exact test statistic. No significant differences were observed in the distributions of CYP2W1 (OR 8.6, 95%CI 0.43-172.5 P>0.05; OR 2.0, 95%CI 0.76-5.4, P>0.05) and CYP4F11 (OR 0.3, 95%CI 0.01-8.4 P>0.05) genotypes between the patients and controls. Only the CYP8A1 CC genotype was detected in patients with BCa and the controls. All polymorphism frequencies were in Hardy-Weinberg Equilibrium (HWE) in the controls (P>0.05). We found a significant association between BCa risk and smoking, use of oral contraceptives or hormonal replacement therapy (HRT), obesity, hyperglycemia, chronic diseases, family history of cancer and menopausal status in the population studied (P<0.05). Tobacco, oral contraceptive or HRT, chronic diseases and obesity or overweight were strongly associated with almost eight, thirty-five, nine and five-fold increased risk for BCa. Tobaco, obesity and hyperglycemia significantly increased the risk of BCa in the patients carrying variant genotypes of CYP2W1 (P<0.05). These results indicate that the CYP2W1 rs3735684, CYP4F11 rs11553651 and CYP8A1 rs56195291 SNPs are not a key risk factor for BCa in Mexican women. This study did not detect an association between the CYP2W1, 4F11 and 8A1 genes polymorphisms and BCa risk in a Mexican population. However, some clinico-pathological risk factors interact with CYP2W1 genotypes and modifies susceptibility to BCa.

Prediction and analysis of the modular structure of cytochrome P450 monooxygenases

Background

Cytochrome P450 monooxygenases (CYPs) form a vast and diverse family of highly variable sequences. They catalyze a wide variety of oxidative reactions and are therefore of great relevance in drug development and biotechnological applications. Despite their differences in sequence and substrate specificity, the structures of CYPs are highly similar. Although being in research focus for years, factors mediating selectivity and activity remain vague.

Description

This systematic comparison of CYPs based on the Cytochrome P450 Engineering Database (CYPED) involved sequence and structure analysis of more than 8000 sequences. 31 structures have been applied to generate a reliable structure-based HMM profile in order to predict structurally conserved regions. Therefore, it was possible to automatically transfer these modules on CYP sequences without any secondary structure information, to analyze substrate interacting residues and to compare interaction sites with redox partners.

Conclusions

Functionally relevant structural sites of CYPs were predicted. Regions involved in substrate binding were analyzed in all sequences among the CYPED. For all CYPs that require a reductase, two reductase interaction sites were identified and classified according to their length. The newly gained insights promise an improvement of engineered enzyme properties for potential biotechnological application. The annotated sequences are accessible on the current version of the CYPED. The prediction tool can be applied to any CYP sequence via the web interface at http://www.cyped.uni-stuttgart.de/cgi-bin/strpred/dosecpred.pl.

Protein-peptide affinity determination using an h/d exchange dilution strategy: application to antigen-antibody interactions

A new methodology using hydrogen/deuterium amide exchange (HDX) to determine the binding affinity of protein-peptide interactions is reported. The method, based on our previously established approach, protein ligand interaction by mass spectrometry, titration, and H/D exchange (PLIMSTEX) [J. Am. Chem. Soc.2003, 125, 5252-5253], makes use of a dilution strategy (dPLIMSTEX) for HDX, using the mass of the peptide ligand as readout. We employed dPLIMSTEX to study the interaction of calcium-saturated calmodulin with the opioid peptide β-endorphin as a model system; the affinity results are in good agreement with those from traditional PLIMSTEX and with literature values obtained by using other methods. We show that the dPLIMSTEX method is feasible to quantify an antigen-antibody interaction involving a 3-nitrotyrosine modified peptide in complex with a monoclonal anti-nitrotyrosine antibody. A dissociation constant in the low nanomolar range was determined, and a binding stoichiometry of antibody/peptide of 1:2 was confirmed. In addition, we determined that the epitope in the binding interface contains a minimum of five amino acids. The dPLIMSTEX approach is a sensitive and powerful tool for the quantitative determination of peptide affinities with antibodies, complementary to conventional immuno-analytical techniques.

Assessment of ligand binding residue predictions in CASP8

Here we detail the assessment process for the binding site prediction category of the eighth Critical Assessment of Protein Structure Prediction experiment (CASP8). Predictions were only evaluated for those targets that bound biologically relevant ligands and were assessed using the Matthews Correlation Coefficient. The results of the analysis clearly demonstrate that three predictors from two groups (Lee and Sternberg) stand out from the rest. A further two groups perform well over subsets of metal binding or nonmetal ligand binding targets. The best methods were able to make consistently reliable predictions based on model structures, though it was noticeable that the two targets that were not well predicted were also the hardest targets. The number of predictors that submitted new methods in this category was highly encouraging and suggests that current technology is at the level that experimental biochemists and structural biologists could benefit from what is clearly a growing field.

Comparison of PGH2 binding site in prostaglandin synthases

BACKGROUND

Prostaglandin H2 (PGH2) is a common precursor for the synthesis of five different Prostanoids via specific Prostanoid Synthases. The binding of this substrate with these Synthases is not properly understood. Moreover, currently no crystal structure of complexes bound with PGH2 has been reported. Hence, understanding the interactions of PGH2 and characterizing its binding sites in these synthases is crucial for developing novel therapeutics based on these proteins as targets.

RESULTS

Shape and physico-chemical properties of the PGH2 binding sites of the four prostanoid synthases were analyzed and compared in order to understand the molecular basis of the specificity. This study provides models with predicted pockets for the binding of PGH2 with PGD, PGE, PGF and PGI Synthases. The results closely match with available experimental data. The comparison showed seven physico-chemical features that are common to the four PGH2 binding sites. However this common pattern is not statistically unique and is not specific enough to distinguish between proteins that can or cannot bind PGH2. A large scale search in ASTRAL data bank, a non redundant Protein Data Bank, for a similar pattern showed the uniqueness of each of the PGH2 binding site in these Synthases.

CONCLUSION

The binding pockets in PGDS, PGES, PGFS and PGIS are unique and do not share significant commonality which can be characterized as a PGH2 binding site. Local comparison of these protein structures highlights a case of convergent evolution in analogous functional sites.

On the mechanism of microsomal prostaglandin E synthase type-2--a theoretical study of endoperoxide reaction with MeS(-)

The reaction pathways of deprotonation versus nucleophilic substitution involving mPGES-2 enzyme catalysis were investigated by ab initio molecular orbital theory calculations for the reaction of methylthiolate with the endoperoxide core of PGH(2) and by the combined quantum mechanical molecular mechanical methods. The calculations showed that deprotonation mechanism is energetically more favorable than the nucleophilic substitution pathway.

Thematic Review Series: Proteomics. An integrated omics analysis of eicosanoid biology

Eicosanoids have been implicated in a vast number of devastating inflammatory conditions, including arthritis, atherosclerosis, pain, and cancer. Currently, over a hundred different eicosanoids have been identified, with many having potent bioactive signaling capacity. These lipid metabolites are synthesized de novo by at least 50 unique enzymes, many of which have been cloned and characterized. Due to the extensive characterization of eicosanoid biosynthetic pathways, this field provides a unique framework for integrating genomics, proteomics, and metabolomics toward the investigation of disease pathology. To facilitate a concerted systems biology approach, this review outlines the proteins implicated in eicosanoid biosynthesis and signaling in human, mouse, and rat. Applications of the extensive genomic and lipidomic research to date illustrate the questions in eicosanoid signaling that could be uniquely addressed by a thorough analysis of the entire eicosanoid proteome.

Identification of PpoA from Aspergillus nidulans as a fusion protein of a fatty acid heme dioxygenase/peroxidase and a cytochrome P450

The homothallic ascomycete Aspergillus nidulans serves as model organism for filamentous fungi because of its ability to propagate with both asexual and sexual life cycles, and fatty acid-derived substances regulate the balance between both cycles. These so-called psi (precocious sexual inducer) factors are produced by psi factor-producing oxygenases (Ppo enzymes). Bioinformatic analysis predicted the presence of two different heme domains in Ppo proteins: in the N-terminal region, a fatty acid heme dioxygenase/peroxidase domain is predicted, whereas in the C-terminal region, a P450 heme thiolate domain is predicted. To analyze the reaction catalyzed by Ppo enzymes, PpoA was expressed in Escherichia coli as an active enzyme. The protein was purified by 62-fold and identified as a homotetrameric ferric heme protein that metabolizes mono- as well as polyunsaturated C(16) and C(18) fatty acids at pH approximately 7.25. The presence of thiolate-ligated heme was confirmed on the basis of sequence alignments and the appearance of a characteristic 450 nm CO-binding spectrum. Studies on its reaction mechanism revealed that PpoA uses different heme domains to catalyze two separate reactions. Within the heme peroxidase domain, linoleic acid is oxidized to (8R)-hydroperoxyoctadecadienoic acid by abstracting a H-atom from C-8 of the fatty acid, yielding a carbon-centered radical that reacts with molecular dioxygen. In the second reaction step, 8-hydroperoxyoctadecadienoic acid is isomerized within the P450 heme thiolate domain to 5,8-dihydroxyoctadecadienoic acid. We identify PpoA as a bifunctional P450 fusion protein that uses a previously unknown reaction mechanism for forming psi factors.

Characterization of the peroxidase mechanism upon reaction of prostacyclin synthase with peracetic acid. Identification of a tyrosyl radical intermediate

Prostacyclin synthase (PGIS) is a membrane-bound class III cytochrome P450 that catalyzes an isomerization of prostaglandin H(2), an endoperoxide, to prostacyclin. We report here the characterization of the PGIS intermediates in reactions with other peroxides, peracetic acid (PA), and iodosylbenzene. Rapid-scan stopped-flow experiments revealed an intermediate with an absorption spectrum similar to that of compound ES (Cpd ES), which is an oxo-ferryl (Fe(IV)O) plus a protein-derived radical. Cpd ES, formed upon reaction with PA, has an X-band (9 GHz) EPR signal of g = 2.0047 and a half-saturation power, P(1/2), of 0.73 mW. High-field (130 GHz) EPR reveals the presence of two species of tyrosyl radicals in Cpd ES with their g-tensor components (g(x), g(y), g(z)) of 2.00970, 2.00433, 2.00211 and 2.00700, 2.00433, 2.00211 at a 1:2 ratio, indicating that one is involved in hydrogen bonding and the other is not. The line width of the g = 2 signal becomes narrower, while its P(1/2) value becomes smaller as the reaction proceeds, indicating migration of the unpaired electron to an alternative site. The rate of electron migration ( approximately 0.2 s(-1)) is similar to that of heme bleaching, suggesting the migration is associated with the enzymatic inactivation. Moreover, a g = 6 signal that is presumably a high-spin ferric species emerges after the appearance of the amino acid radical and subsequently decays at a rate comparable to that of enzymatic inactivation. This loss of the g = 6 species thus likely indicates another pathway leading to enzymatic inactivation. The inactivation, however, was prevented by the exogenous reductant guaiacol. The studies of PGIS with PA described herein provide a mechanistic model of a peroxidase reaction catalyzed by the class III cytochromes P450.

Density functional studies on thromboxane biosynthesis: mechanism and role of the heme-thiolate system

Reaction mechanisms for the isomerization of prostaglandin H(2) to thromboxane A(2), and degradation to 12-L-hydroxy-5,8,10-heptadecatrienoic acid (HHT) and malondialdehyde (MDA), catalyzed by thromboxane synthase, were investigated using the unrestricted Becke-three-parameter plus Lee-Yang-Parr (UB3LYP) density functional level theory. In addition to the reaction pathway through Fe(IV)-porphyrin intermediates, a new reaction pathway through Fe(III)-porphyrin pi-cation radical intermediates was found. Both reactions proceed with the homolytic cleavage of endoperoxide O-O to give an alkoxy radical. This intermediate converts into an allyl radical intermediate by a C-C homolytic cleavage, followed by the formation of thromboxane A(2) having a 6-membered ring through a one electron transfer, or the degradation into HHT and MDA. The proposed mechanism shows that an iron(III)-containing system having electron acceptor ability is essential for the 6-membered ring formation leading to thromboxane A(2). Our results suggest that the step of the endoperoxide O-O homolytic bond cleavage has the highest activation energy following the binding of prostaglandin H(2) to thromboxane synthase.

Modes of heme binding and substrate access for cytochrome P450 CYP74A revealed by crystal structures of allene oxide synthase

Cytochrome P450s exist ubiquitously in all organisms and are involved in many biological processes. Allene oxide synthase (AOS) is a P450 enzyme that plays a key role in the biosynthesis of oxylipin jasmonates, which are involved in signal and defense reactions in higher plants. The crystal structures of guayule (Parthenium argentatum) AOS (CYP74A2) and its complex with the substrate analog 13(S)-hydroxyoctadeca-9Z,11E-dienoic acid have been determined. The structures exhibit a classic P450 fold but possess a heme-binding mode with an unusually long heme binding loop and a unique I-helix. The structures also reveal two channels through which substrate and product may access and leave the active site. The entrances are defined by a loop between beta3-2 and beta3-3. Asn-276 in the substrate binding site may interact with the substrate's hydroperoxy group and play an important role in catalysis, and Lys-282 at the entrance may control substrate access and binding. These studies provide both structural insights into AOS and related P450s and a structural basis to understand the distinct reaction mechanism.

Crystallization and preliminary X-ray analysis of allene oxide synthase, cytochrome P450 CYP74A2, from Parthenium argentatum

Oxylipins are oxygenated derivatives of fatty acids and pivotal signaling molecules in plants and animals. Allene oxide synthase (AOS) is a key cytochrome P450 CYP74 enzyme involved in the biosynthesis of plant oxylipin jasmonates to convert 13(S)-hydroperoxide to allene oxide. Guayule (Parthenium argentatum) AOS, CYP74A2, was expressed in Escherichia coli. Protein was purified using affinity chromatography and size exclusion chromatography, and then crystallized. Two different crystal forms were obtained from 0.2 M (NH(4))H(2)PO(4), 50% MPD, 0.1 M Tris, pH 8.5 at 277 K using the hanging-drop vapor-diffusion method. Preliminary X-ray analysis was carried out, and the crystals were found to belong to the tetragonal space group I422 with cell parameters a = b = 126.5, c = 163.9 A, and the monoclinic space group C2 with cell parameters a = 336.5, b = 184.2, c = 159.0 A, beta = 118.6 degrees . Diffraction data were collected to 2.4 A resolution from a tetragonal form of crystal using a home X-ray source.

Spectroscopic characterization of the oxyferrous complex of prostacyclin synthase in solution and in trapped sol-gel matrix

Prostacyclin synthase (PGIS) is a member of the cytochrome P450 family in which the oxyferrous complexes are generally labile in the absence of substrate. At 4 degrees C, the on-rate constants and off-rate constants of oxygen binding to PGIS in solution are 5.9 x 10(5) m(-1).s(-1) and 29 s(-1), respectively. The oxyferrous complex decays to a ferric form at a rate of 12 s(-1). We report, for the first time, a stable oxyferrous complex of PGIS in a transparent sol-gel monolith. The encapsulated ferric PGIS retained the same spectroscopic features as in solution. The binding capabilities of the encapsulated PGIS were demonstrated by spectral changes upon the addition of O-based, N-based and C-based ligands. The peroxidase activity of PGIS in sol-gel was three orders of magnitude slower than that in solution owing to the restricted diffusion of the substrate in sol-gel. The oxyferrous complex in sol-gel was observable for 24 h at room temperature and displayed a much red-shifted Soret peak. Stabilization of the ferrous-carbon monoxide complex in sol-gel was observed as an enrichment of the 450-nm species over the 420-nm species. This result suggests that the sol-gel method may be applied to other P450s to generate a stable intermediate in the di-oxygen activation.

Structures of prostacyclin synthase and its complexes with substrate analog and inhibitor reveal a ligand-specific heme conformation change

Prostacyclin synthase (PGIS) is a cytochrome P450 (P450) enzyme that catalyzes production of prostacyclin from prostaglandin H(2). PGIS is unusual in that it catalyzes an isomerization rather than a monooxygenation, which is typical of P450 enzymes. To understand the structural basis for prostacyclin biosynthesis in greater detail, we have determined the crystal structures of ligand-free, inhibitor (minoxidil)-bound and substrate analog U51605-bound PGIS. These structures demonstrate a stereo-specific substrate binding and suggest features of the enzyme that facilitate isomerization. Unlike most microsomal P450s, where large substrate-induced conformational changes take place at the distal side of the heme, conformational changes in PGIS are observed at the proximal side and in the heme itself. The conserved and extensive heme propionate-protein interactions seen in all other P450s, which are largely absent in the ligand-free PGIS, are recovered upon U51605 binding accompanied by water exclusion from the active site. In contrast, when minoxidil binds, the propionate-protein interactions are not recovered and water molecules are largely retained. These findings suggest that PGIS represents a divergent evolution of the P450 family, in which a heme barrier has evolved to ensure strict binding specificity for prostaglandin H(2), leading to a radical-mediated isomerization with high product fidelity. The U51605-bound structure also provides a view of the substrate entrance and product exit channels.

Unusual regioselectivity and active site topology of human cytochrome P450 2J2

The oxidation of six derivatives of terfenadone by recombinant human CYP2J2 (CYP = cytochrome P450) was studied by high-performance liquid chromatography coupled to mass spectrometry (MS) using tandem MS techniques and by 1H NMR spectroscopy. CYP2J2 exhibited a surprising regioselectivity in favor of the hydroxylation of the substrate terminal chain at the weakly reactive homobenzylic position. In contrast, hydroxylation of the same substrates by CYP3A4 mainly occurred on the most chemically reactive sites of the substrates (N-oxidation and benzylic hydroxylation). A 3D homology model of CYP2J2 was constructed using recently published structures of CYP2A6, CYP2B4, CYP2C8, CYP2C9, and CYP2D6 as templates. In contrast with other CYP2 structures, it revealed an active site cavity with a severely restricted access of substrates to the heme through a narrow hydrophobic channel. Dynamic docking of terfenadone derivatives in the CYP2J2 active site allowed one to interpret the unexpected regioselectivity of the hydroxylation of these substrates by CYP2J2, which is mainly based on this restricted access to the iron. The structural features that have been found to be important for recognition of substrates or inhibitors by CYP2J2 were also interpreted on the basis of CYP2J2-substrate interactions in this model.

Conformational dynamics in the F/G segment of CYP51 from Mycobacterium tuberculosis monitored by FRET

A cysteine was introduced into the FG-loop (P187C) of CYP51 from Mycobacterium tuberculosis (MT) for selective labeling with BODIPY and fluorescence energy transfer (FRET) analysis. Förster radius for the BODIPY-heme pair was calculated assuming that the distance between the heme and Cys187 in solution corresponds to that in the crystal structure of ligand free MTCYP51. Interaction of MTCYP51 with azole inhibitors ketoconazole and fluconazole or the substrate analog estriol did not influence the fluorescence, but titration with the substrate lanosterol quenched BODIPY emission, the effect being proportional to the portion of substrate bound MTCYP51. The detected changes correspond to approximately 10A decrease in the calculated distance between BODIPY-Cys187 and the heme. The results confirm (1) functional importance of conformational motions in the MTCYP51 F/G segment and (2) applicability of FRET to monitor them in solution.

Reaction mechanisms of 15-hydroperoxyeicosatetraenoic acid catalyzed by human prostacyclin and thromboxane synthases

Prostacyclin synthase (PGIS) and thromboxane synthase (TXAS) are atypical cytochrome P450s. They do not require NADPH or dioxygen for isomerization of prostaglandin H(2) (PGH(2)) to produce prostacyclin (PGI(2)) and thromboxane A(2) (TXA(2)). PGI(2) and TXA(2) have opposing actions on platelet aggregation and blood vessel tone. In this report, we use a lipid hydroperoxide, 15-hydroperoxyeicosatetraenoic acid (15-HPETE), to explore the active site characteristics of PGIS and TXAS. The two enzymes transformed 15-HPETE not only into 13-hydroxy-14,15-epoxy-5,8,11-eicosatrienoic acid (13-OH-14,15-EET), like many microsomal P450s, but also to 15-ketoeicosatetraenoic acid (15-KETE) and 15-hydroxyeicosatetraenoic acid (15-HETE). 13-OH-14,15-EET and 15-KETE result from homolytic cleavage of the O-O bond, whereas 15-HETE results from heterolytic cleavage, a common peroxidase pathway. About 80% of 15-HPETE was homolytically cleaved by PGIS and 60% was homolytically cleaved by TXAS. The V(max) of homolytic cleavage is 3.5-fold faster than heterolytic cleavage for PGIS-catalyzed reactions (1100 min(-1)vs. 320 min(-1)) and 1.4-fold faster for TXAS (170 min(-1)vs. 120 min(-1)). Similar K(M) values for homolytic and heterolytic cleavages were found for PGIS ( approximately 60 microM 15-HPETE) and TXAS ( approximately 80 microM 15-HPETE), making PGIS a more efficient catalyst for the 15-HPETE reaction.