Comparison of PGH2 binding site in prostaglandin synthases
© Paragi-Vedanthi and Doble; licensee BioMed Central Ltd. 2010
Published: 18 January 2010
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.
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.
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
The Cyclooxygenase (COX) pathway is an important part of the arachidonic acid (AA) metabolism, generating five primary prostanoids. The biosynthesis of these prostanoids involves a sequence of three-steps namely 1) Release of arachidonic acid from phospholipids by secretory, cytoplasmic or from both types of phospholipase A2 (sPLA2 and cPLA2), 2) Oxygenation of AA by COX enzymes to form prostaglandin endoperoxide H2 (PGH2), and 3) the subsequent conversion of PGH2 to Prostaglandin D2 (PGD2), Prostaglandin E2 (PGE2), Prostaglandin F2α (PGF2α), prostacyclin (PGI2), and Thromboxane T2 (TXA2) via seven specific synthases [1, 2].
PGD Synthase is responsible for the production of PGD2 as an allergy or inflammation mediator in mast and Th2 cells . There are 3 isoforms of PGE Synthase (PGES) namely microsomal PGE Synthase -1 (mPGES-1), microsomal PGE Synthase -2 (mPGES-2) and cytoplasmic PGE Synthase (cPGES) responsible for the production of PGE2, which is an ultimate mediator of pain and inflammation. PGE2 also plays a critical role in regulating renal function and facilitating reproduction . Prostaglandin F2α produced from PGF Synthase (PGFS) is a hormone-like substance participating in a wide range of body functions including the contraction and relaxation of smooth muscle, the dilation and constriction of blood vessels, control of blood pressure, and modulation of inflammation. PGF2α is used for the induction of abortion, for evacuation of the uterus after a missed abortion . PGI2 produced by PGIS, and TXA2 produced by TXAS, are critical for the maintenance of homeostasis in the vascular tissue [6, 7]. Since these five synthases, are involved in various important biological processes, they are potential drug targets and drugs are already in the market for the inhibition of PGDS, PGFS, PGIS and TXAS. mPGES-1 is being sought after as a novel target to relieve pain and inflammation after the withdrawal of popular COX-2 inhibitors from the market .
Understanding the interactions of PGH2 with these synthases and characterizing their binding sites is crucial for developing novel drugs and also to check for cross reactivity. PGH2 is an unstable compound and there are no structures of synthases available in the Protein Data Bank (PDB)  with it. In this paper, the PGH2 binding sites in these proteins were predicted using the PatchDock algorithm . The predicted binding sites were then compared using MultiBind , a multiple binding site alignment tool to look for common pattern which might help us to characterize a PGH2 binding site.
The crystal structure of four of the proteins namely PGDS, PGES, PGFS and PGIS are available in the PBD, while such a structure is not available for mPGES-1, cPGES and TXAS. For the purpose of docking studies the following structures were used: (i) PGDS - The structure of human hematopoietic prostaglandin D synthase complexed with HQL-79 (PDB: 2cvd, ). (ii) PGES - The structure of Microsomal prostaglandin E synthase type-2 (PDB: 1z9h, ). (iii)PGFS - The structure of prostaglandin F synathase containing bimatoprost (PDB: 2f38, ) and (iv) PGIS - The structure of human prostacyclin synthase in complex with inhibitor minoxidil (PDB: 3b6h, ).
The dockings of PGH2 with these synthases were performed to predict the putative binding site in the proteins. Docking models are obtained using the PatchDock algorithm . This software takes two molecules as input and computes the three-dimensional transformation of one of them with respect to the other with the goal of maximizing the surface shape complementarities and at the same time minimizing the number of steric clashes. Given two molecules, PatchDock first divides their surfaces into patches according to their surface shape, such as concave, convex, or flat. Then, it applies the Geometric Hashing algorithm to match the concave patches of one with the convex patches of the other protein and flat patches with flat patches and generates a set of candidate transformations. A set of scoring functions based on the shape complementarities and the atomic desolvation energy of the transformed complex is evaluated. Finally, redundant solutions are discarded by the application of a RMSD (root-mean-square deviation) clustering. This program is tested and shown to successfully predict protein interactions for many examples [16–19].
Binding site alignment
The alignment of these predicted binding sites for PGH2 in the four synthases was performed using the MultiBind algorithm developed by Shulmana et. al .
This algorithm performs multiple alignments of the binding sites and recognizes the structurally conserved physicochemical and geometrical patterns that may be responsible for the binding. The physicochemical properties considered by the software are hydrophobic, aliphatic (ALI) and aromatic interactions (PII), hydrogen bond donors (DON), hydrogen bond acceptors (ACC), and mixed donor/acceptors (DAC). The algorithm finds a set of transformations which will superimpose the binding sites in a manner that will maximize the physicochemical score of the matched properties. This alignment between protein binding sites is performed even in the absence of overall sequence, fold, or binding partner similarity, and also it does not consider the location of the binding partners. The scoring function and the algorithm of MultiBind are described in detail elsewhere .
Evaluation of common binding patterns
The frequency of random occurrence of structural patterns as recognized by MultiBind is searched with proteins in the ASTRAL dataset (V 1.73) [20, 21]. This dataset consists of all known protein structures that have less than 40% sequence identity and hence it represents a non-redundant group. This dataset contains only the structures of PGDS and mPGES-2 that are known to bind to PGH2. The other two structures namely PGFS and PGIS were added to make the dataset of 7649 structures in which four of them were known to bind PGH2. The structural pattern was chosen from the pseudocenters of the first molecule in the input order. Each pattern that was recognized on the surface of some other protein was scored using the physicochemical scoring function of MultiBind as explained before. The frequency of occurrence of a pattern was calculated as the ratio between the numbers of times it was observed with a score higher than a reference score, relative to the total number of searched proteins. The reference score is defined to be the score of the outlier, i.e., lowest score of the most different binding site that participated in the pattern construction with MultiBind.(Example: If A, B, C and D are the protein compared, and the binding sites of A and D differ the most with the least MultiBind Score then that is taken as the reference score for comparisons of ABCD, ABD, ACD and AD with the ASTRAL dataset). The obtained ratio represents the estimation of the chances for a random occurrence of the recognized pattern. Using the score of the outlier as a reference score provides the highest possible ratio and the worst case estimation of the most frequent pattern . The ratio of the number of similar patterns observed relative to the size of the searched dataset provides an estimation of the probability of observing such a pattern by chance, on a randomly selected protein. The lower the frequency of occurrence, rare is the pattern.
Docking models of PGH2 with four synthases
PatchDock successfully detected the surface pockets of PGDS, PGES, PGFS and PGIS and they are in agreement with previously published data [22–30]on the putative binding site of PGH2. In all these analysis no apriori information was used as an input for the docking algorithms, i.e., the surface pockets on the receptor molecules were detected automatically.
Binding site alignment and analysis
After obtaining the docking models and information about the PGH2 binding pocket, the goal was to compare these four predicted binding pockets and to determine the common features which facilitate the binding of the same substrate, PGH2.
The four binding sites compared here have different overall sequences and conserved residue patterns and are also structurally not related. Due to the above-mentioned differences, the proteins in this study cannot be aligned by standard alignment methods, that assume similarity of either sequence or backbone patterns. Thus, to compare between the predicted binding sites of the modelled complexes, we used MultiBind which performs a multiple structure alignment between protein binding sites, in the absence of overall sequence, fold, or binding partner similarity and recognizes the common spatial arrangements of physicochemical properties shared between the compared binding sites. The ligand from the crystal structures were used to extract the binding site pocket with the MULTIBIND algorithm rather with the PGH2 docked models as the algorithm do not accept docked models as input.
Multiple alignment of all the four binding sites
Details of the common pattern calculated between the four PGH2 binding sites using MultiBind. PII (aromatic) interactions, hydrogen bond acceptor (ACC), or mixed donor-acceptor (DON, DAC) and Aliphatic interactions(ALI).
Site 1: PGES/1z9hB
Site 2: PGDS/2cvdD
Site 4: PGIS/3b6hA
Multiple alignments of three binding sites
The alignment between all the four binding sites resulted in a common pattern with seven physiochemical properties that is not unique enough for identification of the binding of PGH2. Therefore binding sites were aligned, eliminating one synthase at a time from the MultiBind during run time.
Comparison between three predicted binding sites of PGH2 at a time revealed a pattern of between 6-8 common physicochemical properties. The frequency of occurrence of this pattern in the ASTRAL dataset ranged between 17-38% showing that the frequency of occurrence of this pattern in this dataset is also too high to be statistically significant and once again it can be concluded that it is not specific enough to distinguish between proteins that can or cannot bind PGH2.
Pairwise alignments of PGH2 binding sites using MultiBind
No. of detected features
Based on the MultiBind similarity score, PGDS-PGFS is ranked second. PGFS is a dual acting enzyme leading to the formation of both PGF2α from PGH2 and 9α,11β-PGF2 (PGF2αβ) from PGD2. It can bind to both PGH2 and PGD2 . The former is converted to PGD2 by PGDS; the later remains bound to the enzyme before being released. Having two common binding partners possibly explains their high similarity. It is surprising to note that both PGDS and PGES belonging to the same family with similar catalytic mechanism has the least number (11) of common features and are among the most different from each other when compared with other synthases. This is followed by PGIS-PGDS having 12 common features. mPGES-2 on the other hand is found to be more similar to PGIS which is a heme bound enzyme, than to PGDS. This can be explained based on a more recently solved structure in which PGES (PDBID: 2pbj) is found to contain glutathione (GSH) and heme bound to it and it is involved in degradation reactions similar to that of cytochrome P450. It degrades PGH2 into 12(S)-hydroxy-5(Z), 8(E),10(E)-heptadecatrienoic acid and malondialdehyde rather than converting it to PGE2 .
But evolution of similar enzymatic function on different structural frameworks is not an entirely uncommon event. A classical example is that of serine proteases: the Ser-His-Asp triad is present in an almost identical three-dimensional constellation in the distinctly different structural frames of trypsin and subtilisin (and their relatives) [32, 33]. Another example is each of the three families of sugar kinases which appear to have a distinct three-dimensional fold, and conserved sequence patterns are strikingly different for the three families. Yet each catalyzes chemically equivalent reactions on similar or identical substrates. The enzymatic function of sugar phosphorylation appears to have evolved independently on the three distinct structural frameworks, by convergent evolution. Another aspect which needs consideration is that these proteins can rearrange and undergo conformational changes to accommodate the substrate. In practice, both the side-chains and the protein backbone can undergo conformational changes upon substrate binding. Even the PGH2 molecule with 14 torsional degrees of freedom of rotation is highly flexible and can fit into different active sites differently. Also both the docking and the alignment algorithms used in the current study considers rigid conformations and do not address the possibility of protein and substrate flexibility.
In summary, the Prostanoid Synthases present a remarkable diversity of specificities for the binding of PGH2. The discovery of this striking molecular dissimilarity, associated to a functional substrate similarity, may help in suggesting new experiments aimed at a deeper understanding of the cross reactivity of Prostanoid synthases known to be involved in many important biological processes and human diseases.
A computational approach was employed to understand the interaction of PGH2 with the prostaglandin synthases. Docking models were consistent with the available experimental data for the interaction of PGH2 with the synthases. The spatial and physicochemical properties of the suggested binding sites were compared. A patterns common to all the four synthases was detected but it was not specific enough and was not likely to represent the features essential for the binding of PGH2. The alignment results suggest that the PGH2 binding sites are different on different proteins and they also have no close similarity with any other binding site found in the proteins of the ASTRAL dataset.
Padma Priya P V. thanks the Department of Science and Technology, Government of India for financial support.
This article has been published as part of BMC Bioinformatics Volume 11 Supplement 1, 2010: Selected articles from the Eighth Asia-Pacific Bioinformatics Conference (APBC 2010). The full contents of the supplement are available online at http://www.biomedcentral.com/1471-2105/11?issue=S1.
- Helliwell RJA, Adams LF, Mitchell MD: Prostaglandin synthases: recent developments and a novel hypothesis. Prostaglandins, Leukotrienes and Essential Fatty Acids 2004, 70: 101–113. 10.1016/j.plefa.2003.04.002View ArticleGoogle Scholar
- William SL: Prostanoid biosynthesis and mechanisms of action. Am J Physiol Renal Physiol 1992, 263: 181–191.Google Scholar
- Murakami M, Nakatani Y, Tanioka T, Kudo I: Prostaglandin E synthase. Prostaglandins Other Lipid Media 2002, 68: 383–99. 10.1016/S0090-6980(02)00043-6View ArticleGoogle Scholar
- Watanabe K: Prostaglandin F synthase. Prostaglandins & other Lipid Mediators 2002, 68: 401–407. 10.1016/S0090-6980(02)00044-8View ArticleGoogle Scholar
- Kanaoka Y, Urade Y: Hematopoietic prostaglandin D synthase. Leukotrienes and Essential Fatty Acids 2003, 69: 63–167.View ArticleGoogle Scholar
- Vane J, Corin RE: Prostacyclin: A Vascular Mediator. Eur J Vasc Endovasc Surg 2003, 26: 571–578. 10.1016/S1078-5884(03)00385-XView ArticlePubMedGoogle Scholar
- Wang LH, Kulmacz RJ: Thromboxane synthase: structure and function of protein and gene. Prostaglandins & other Lipid Mediators 2002, 68: 409–422. 10.1016/S0090-6980(02)00045-XView ArticleGoogle Scholar
- Scholich K, Geisslinger G: Is mPGES-1 a promising target for pain therapy? Trends in Pharmacological Sciences 2006, 27: 399–401. 10.1016/j.tips.2006.06.001View ArticlePubMedGoogle Scholar
- Berman HM, Westbrook J, Feng G, Gilliland G, Bhat TN, Weissig H, Shindyalov IN, Bourne PE: The Protein Data Bank. Nucleic Acids Research 2000, 28: 235–242. 10.1093/nar/28.1.235PubMed CentralView ArticlePubMedGoogle Scholar
- Schneidman-Duhovny D, Inbar Y, Nussinov R, Wolfson HJ: PatchDock and SymmDock: servers for rigid and symmetric docking. Nucleic Acids Res 2005, 33: W363–367. 10.1093/nar/gki481PubMed CentralView ArticlePubMedGoogle Scholar
- Shatsky M, Shulman-peleg A, Nussinov R, Wolfson HJ: The Multiple Common Point Set Problem and Its Application to Molecule Binding Pattern Detection. Journal of Computational Biology 2006, 13: 407–428. 10.1089/cmb.2006.13.407View ArticlePubMedGoogle Scholar
- Aritake K, Kado Y, Inoue T, Miyano M, Urade Y: Structural and Functional Characterization of HQL-79, an Orally Selective Inhibitor of Human Hematopoietic Prostaglandin D Synthase. The Journal of Biological Chemistry 2006, 281: 5277–15286. 10.1074/jbc.M506431200View ArticleGoogle Scholar
- Yamada T, Komoto J, Watanabe K, Ohmiya Y, Takusagawa F: Crystal Structure and Possible Catalytic Mechanism of Microsomal Prostaglandin E Synthase Type 2(mPGES-2). J Mol Biol 2005, 348: 1163–1176. 10.1016/j.jmb.2005.03.035View ArticlePubMedGoogle Scholar
- Komoto J, Yamada T, Watanabe K, Woodward DF, Takusagawa F: Prostaglandin F2R Formation from Prostaglandin H2 by Prostaglandin F Synthase (PGFS): Crystal Structure of PGFS Containing Bimatoprost. Biochemistry 2006, 45: 1987–1996. 10.1021/bi051861tView ArticlePubMedGoogle Scholar
- Li YC, Chiang CW, Yeh HC, Hsu PY, Whitby FG, Wang LH, Chan NL: Structures of Prostacyclin Synthase and Its Complexes with Substrate Analog and Inhibitor Reveal a Ligand-specific Heme Conformation Change. J Biol Chem 2008, 283: 2917–2926. 10.1074/jbc.M707470200PubMed CentralView ArticlePubMedGoogle Scholar
- Benyamini H, Shulman-Peleg A, Wolfson HJ, Belgorodsky B, Fadeev L, Gozin M: Interaction of C60-Fullerene and Carboxyfullerene with Proteins: Docking and Binding Site Alignment. Bioconjugate Chem 2006, 17: 378–386. 10.1021/bc050299gView ArticleGoogle Scholar
- Inbar Y, Schneidman-Duhovny D, Halperin I, Oron A, Nussinov R, Wolfson HJ: Approaching the CAPRI challenge with an efficient geometry-based docking. Proteins 2005, 60: 217–223. 10.1002/prot.20561View ArticlePubMedGoogle Scholar
- Gidalevitz T, Biswas C, Ding H, Schneidman-Duhovny D, Wolfson HJ, Stevens F, Radford S, Argon Y: Identification of the N-terminal peptide binding site of glucoseregulated protein 94. J Biol Chem 2004, 279: 16543–16552. 10.1074/jbc.M313060200View ArticlePubMedGoogle Scholar
- Schneidman-Duhovny D, Inbar Y, Polak V, Shatsky M, Halperin I, Benyamini H, Barzilai A, Dror O, Haspel N, Nussinov R, Wolfson HJ: Taking geometry to its edge: fast unbound rigid (and hinge-bent) docking. Proteins 2003, 52: 107–112. 10.1002/prot.10397View ArticlePubMedGoogle Scholar
- Brenner SE, Koehl P, Levitt M: The ASTRAL compendium for protein structure and sequence analysis. Nucleic Acids Res 2000, 28: 254–256. 10.1093/nar/28.1.254PubMed CentralView ArticlePubMedGoogle Scholar
- Chandonia JM, Hon G, Walker NS, Lo Conte L, Koehl P, Levitt M, Brenner SE: The ASTRAL Compendium in 2004. Nucleic Acids Res 2004, 32: D189–192. 10.1093/nar/gkh034PubMed CentralView ArticlePubMedGoogle Scholar
- Inoue T, Okano Y, Kado Y, Aritake K, Irikura D, Uodome N, Okazaki N, Kinugasa S, Shishitani H, Matsumura H, et al.: First Determination of the Inhibitor Complex Structure of Human Hematopoietic Prostaglandin D Synthase. J Biochem 2004, 135: 279–283. 10.1093/jb/mvh033View ArticlePubMedGoogle Scholar
- Inoue T, Irikura D, Okazaki N, Kinugasa S, Matsumura H, Uodome N, Yamamoto M, Kumasaka T, Miyano M, Kai Y, et al.: Mechanism of metal activation of human hematopoietic prostaglandin D synthase. Nature structural biology 2003, 10: 291–296. 10.1038/nsb907View ArticlePubMedGoogle Scholar
- Kanaoka Y, Ago H, Inagaki E, Nanayama T, Miyano M, Kikuno R, Fujii Y, Eguchi N, Toh H, Urade Y, et al.: Cloning and Crystal Structure of Hematopoietic Prostaglandin D Synthase. Cell 1997, 90: 1085–1095. 10.1016/S0092-8674(00)80374-8View ArticlePubMedGoogle Scholar
- Huang X, Yan W, Gao D, Tong M, Tai HH, Zhan CG: Structural and functional characterization of human microsomal prostaglandin E synthase-1 by computational modeling and site-directed mutagenesis. Bioorganic & Medicinal Chemistry 2006, 14: 3553–3562. 10.1016/j.bmc.2006.01.010View ArticleGoogle Scholar
- Watanabe K, Ohkubo H, Niwa H, Tanikawa N, Koda N, Ito S, Ohmiya Y: Essential 110Cys in active site of membrane-associated prostaglandin E synthase-2. Biochemical and Biophysical Research Communications 2003, 306: 577–581. 10.1016/S0006-291X(03)01025-8View ArticlePubMedGoogle Scholar
- Suzuki T, Fujii Y, Miyano M, Chen LY, Takahashi T, Watanabe K: cDNA Cloning, Expression, and Mutagenesis Study of Liver-type Prostaglandin F Synthase. The Journal of Biological Chemistry 1999, 274: 1241–248.Google Scholar
- Komoto J, Yamada T, Watanabe K, Takusagawa F: Crystal Structure of Human Prostaglandin F Synthase (AKR1C3). Biochemistry 2004, 43: 2188–2198. 10.1021/bi036046xView ArticlePubMedGoogle Scholar
- Shyue SK, Ruan KH, Wang LH, Wu KK: Prostacyclin Synthase Active Sites, Identification by moleculer modelling-guided site-directed mutagenesis. The Journal of Biological Chemistry 1997, 272: 3657–3662. 10.1074/jbc.272.6.3657View ArticlePubMedGoogle Scholar
- Chiang CW, Yeh HC, Wang LH, Chan NL: Crystal Structure of the Human Prostacyclin Synthase. J Mol Biol 2006, 364: 266–274. 10.1016/j.jmb.2006.09.039PubMed CentralView ArticlePubMedGoogle Scholar
- Yamada T, Takusagawa F: PGH2 degradation pathway catalyzed by GSH-heme complex bound microsomal prostaglandin E2 synthase type 2: the first example of a dual-function enzyme. Biochemistry 2007, 46: 8414–8424. 10.1021/bi700605mView ArticlePubMedGoogle Scholar
- Wright CS, Alden RA, Kraut J: Structure of subtilisin BPN at 2.5 A resolution. Nature 1969, 221: 235–242. 10.1038/221235a0View ArticlePubMedGoogle Scholar
- Brok P, Sander C, Valencia A: Convergent evolution of similar enzymatic function on different protein folds: The hexokinase, ribokinase, and galactokinase families of sugar kinases. Protein Science 1993, 2: 31–40.View ArticleGoogle Scholar
- Park JY, Pillinger MH, Abramson SB: Prostaglandin E2 synthesis and secretion:The role of PGE2 synthases. Clinical Immunology 2006, 119: 229–240. 10.1016/j.clim.2006.01.016View ArticlePubMedGoogle Scholar
- Ruan KH, Li P, Kulmacz RJ, Wu KK: Characterization of the Structure and Membrane Interaction of NH2-terminal Domain of Thromboxane A Synthase. JBC 1994, 269: 2093–2094.Google Scholar
- Delano WL: The PyMOL Molecular Graphics System. DeLano Scientific LLC, San Carlos, CA; 2002.Google Scholar
- Kanehisa M, Goto S: KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Res 2000, 28: 27–30. 10.1093/nar/28.1.27PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.