- Open Access
Applying internal coordinate mechanics to model the interactions between 8R-lipoxygenase and its substrate
BMC Bioinformatics volume 11, Article number: S2 (2010)
Lipoxygenases (LOX) play pivotal roles in the biosynthesis of leukotrienes and other biologically active potent signalling compounds. Developing inhibitors for LOX is of high interest to researchers. Modelling the interactions between LOX and its substrate arachidonic acid is critical for developing LOX specific inhibitors. Currently, there are no LOX-substrate structures. Recently, the structure of a coral LOX, 8R-LOX, which is 41% sequence identical to the human 5-LOX was solved to 1.85Å resolution. This structure provides a foundation for modelling enzyme-substrate interactions.
In this research, we applied a computational method, Internal Coordinate Mechanics (ICM), to model the interactions between 8R-LOX and its substrate arachidonic acid. Docking arachidonic acid to 8R-LOX was performed. The most favoured docked ligand conformations were retained. We compared the results of our simulation with a proposed model and concluded that the binding pocket identified in this study agrees with the proposed model partially.
The results showed that the conformation of arachidonic acid docked into the ICM-identified docking site has less energy than that docked into the manually defined docking site for pseudo wild type 8R-LOX. The mutation at I805 resulted in no docking pocket found near Fe atom. The energy of the arachidonic acid conformation docked into the manually defined docking site is higher in mutant 8R-LOX than in wild type 8R-LOX. The arachidonic acid conformations are not productive conformations.
We concluded that, for the wild type 8R-LOX, the conformation of arachidonic acid docked into the ICM-identified docking site is more stable than that docked into the manually defined docking site. Mutation affects the structure of the putative active site pocket of 8R-LOX, and leads no docking pockets around the catalytic Fe atom. The docking simulation in a mutant 8R-LOX demonstrated that the structural change due to the mutation impacts the enzyme activity. Further research and analysis is required to obtain the 8R-LOX-substrate model.
Lipoxygenases (LOX) are non-heme iron dioxygenases that catalyze the stereo- and regio- specific formation of fatty acid hydroperoxides from polyunsaturated fatty acids, which are commonly found in plants and animals [1, 2]. Human lipoxygenases play pivotal roles in the biosynthesis of leukotrienes and other biologically active potent signalling compounds in the inflammatory response . Due to this, human lipoxygenases are targets for developing inhibitors to modulate the effects of the potent signalling compounds [4, 5]. Leukotrienes are derived from arachidonic acid, which is the substrate of lipoxygenase . An understanding of substrate arachidonic acid recognition is helpful in the design of enzyme-specific inhibitors . However, there are no animal LOX structures available to provide a model for how the substrate binds in the active site. The structure of a coral lipoxygenase, 8R-LOX, from Plexaura homomalla, with 41% identical to a human lipoxygenase, 5-LOX, was recently solved to 1.85 Å . There is no human isozyme with better than 40% sequence identity to human 5-LOX. The high resolution structure of 8R-LOX can provide a strong foundation for modelling enzyme-substrate interactions.
Many computational methods have been used to model the interactions between enzyme and substrate. Internal Coordinate Mechanics (ICM) is one of the computational methods used by researchers to predict enzyme-substrate interactions. ICM is a stochastic global optimization methodology with biased probability Monte Carlo procedure which can be used to effectively model substrate docking and predict structure [8–17]. It has been used to identify the active site of enzyme and acknowledged to be an accurate predictive tool of binding geometry today [13, 18].
In this research, we will use the high resolution structure of 8R-LOX as the model protein and apply ICM to model how the substrate arachidonic acid, interacts with 8R-LOX. This research will lead us to develop substrate-LOX models for the lipoxygenase superfamily and facilitate the development of specific anti-LOX inhibitors.
Retrieval of high resolution structure of 8R-LOX and pre-processing
The high resolution 1.85Å structure of pseudo-wild type 8R-LOX (psWT) and 1.9Å structure of a deletion mutant of psWT 8R-LOX (I805W:psWT) were obtained from RCSB Protein Data Bank with access code 3FG1 and 3FG3, respectively . ICM Pro version 3.6 was used for docking arachidonic acid into the protein.
First, the pdb files of psWT and I805W:psWT were converted into ICM objects, by optimizing hydrogens, optimizing side-chains and deleting water molecules. Since there are four identical molecules in psWT and I805W:psWT and we would need only one to perform the docking, we removed three of them to avoid the noises and inconveniences caused by the existence of the other three molecules. For the same purpose, the unrelated small molecules were also removed from the ICM objects of psWT and I805W:psWT (Figure 1).
Docking arachidonic acid into 8R-LOX
ICM Pro version 3.6 was used to perform docking the substrate, arachidonic acid into psWT and I805W:psWT. We first set psWT and I805W:psWT as the receptors and found the potential ligand binding pockets in psWT and I805W:psWT. Tolerance value we used was 3 to avoid missing any docking pockets (the default value is 5).
Of the potential ligand binding pockets found, the one which is closest to the catalytic Fe was selected to perform the docking. We also manually defined a docking site for psWT and I805W:psWT based on the binding site proposed by Neau . This manually defined docking site includes 13 residues (Table 1).
Reviewing and adjusting docking site was done. Then the receptor map was generated. This step was to construct energy map of the environment within the docking area box. Interactive docking was performed for the chosen docking pocket and manually defined docking pocket. The thoroughness which represents the length of simulation was set as 1 as suggested by ICM Pro. Properties of conformations were saved for viewing and analysis.
Docking of arachidonic acid into psWT 8R-LOX
With tolerance value 3, we found 26 potential ligand binding pockets (Figure 2). Of the 26 potential ligand binding pockets, the one which was closest to Fe atom and located in the target area was chosen to dock arachidonic acid into psWT 8R-LOX (Figure 2). All other binding pockets were not considered for doing docking because either they were too far from the catalytic Fe atom or the volume was too big to properly accommodate the small ligand, arachidonic acid. The volume of the selected binding pocket is 171.162 Å3, and the radius of it is 3.4Å. This pocket is composed of 15 residues plus Fe atom. The 15 residues are Q753, L754, H757, L758, H762, I796, D797, G800, L804, I810, V811, L815, T996, L1000 and I1066 (Table 1).
The docking simulation resulted in seventy-three conformations of arachidonic acid docked in the psWT 8R-LOX. The conformation of arachidonic acid with the lowest energy (-63.623kcal/mol) is shown in Figure 3. From the conformation, we can see arachidonic acid is docked into the binding pocket with carboxyl terminal pointing inside of psWT 8R-LOX. However, there are other conformations with the carboxyl terminal pointing to the outside of psWT 8R-LOX although the energy is high.
Interactive docking of arachidonic acid into the manually defined docking site (Table 1) resulted in seventy-six conformations of arachidonic acid docked into the psWT 8R-LOX. We found that the conformation of arachidonic acid with the lowest energy is not located around the Fe atom. Putting together information about location, energy, and feasibility, we identified one conformation as the most favoured conformation (Figure 4). This conformation is located near the Fe and with low energy (2.706kcal/mol). In this conformation, the carboxyl terminal of arachidonic acid points to outside of 8R-LOX, and the methyl C terminal is located inside of the psWT 8R-LOX.
Docking of arachidonic acid to mutant I805W:psWT 8R-LOX
We found twenty-five potential ligand binding pockets in I805W:psWT 8R-LOX (Figure 5). There are none of these pockets enclosing the catalytic Fe atom. Data shows that the closest pocket to Fe atom is as far as 4.5Å from the catalytic Fe atom, not containing Fe atom. Since no pocket sits around Fe atom, no ICM docking was performed for any of these 25 pockets.
However, we did run ICM docking using the manually defined docking site for I805W:psWT (Table 1). Interactive docking simulation for the manually defined docking site resulted in 73 conformations of arachidonic acid docked to I805W:psWT 8R-LOX. The most favoured arachidonic acid conformation is near the Fe atom, 3.1 Å, and has the lowest energy, 39.947kcal/mol (Figure 6).
Conformations of arachidonic acid docked into psWT 8R-LOX
Comparison of the conformations of arachidonic acid docked into the ICM-identified docking site and that docked into the manually defined docking site shows that the location of arachidonic acid for both conformations is overlapped. However, there is difference between the two conformations. In the conformation of arachidonic acid docked into the ICM-identified docking site, the carboxyl terminal of arachidonic acid points inside of 8R-LOX, while in the conformation of arachidonic acid docked into the manually defined docking site, the carboxyl terminal of arachidonic acid points outside of 8R-LOX (Figures 3 and 4).
The conformation of arachidonic acid in Figure 4 implies a possible access for arachidonic acid to Fe atom from ambience. However, if we compare the energy for the two conformations, the conformation of arachidonic acid in Figure 3 is much lower than that in Figure 4, indicating the conformation of arachidonic acid in Figure 3 is more stable and favoured. We realized that there are 6 residues in common for ICM-identified docking pocket and the manually defined docking pocket. They are L758, I796, G800, L804, L1000, and I1066. These common residues explain why the locations of the two conformations are overlapped. The common residues for the two docking pockets indicate that the binding site proposed by Neau et al  makes sense. Our simulation confirmed the proposed binding site partially. Further simulation has been designed to verify/modify the proposed binding site in 8R-LOX.
Conformations of arachidonic acid docked into psWT 8R-LOX and I805W:psWT 8R-LOX
We found that 25 potential docking pockets were identified in I805W:psWT (Figure 5), which is 1 less than that in psWT 8R-LOX (Figure 2). We can find that the pockets in Figure 2 and Figure 5 are basically the same except the area around Fe atom, indicating that the mutation at I805 changed the 3-D structure of the area around the catalytic Fe atom in 8R-LOX. The interesting thing is that there is none of the 25 pockets containing the catalytic Fe atom. This tells us that the mutation made at I805 caused a fundamental change in the binding site of 8R-LOX, leaving no docking pocket around Fe atom.
The docking results for the manually defined docking site show that the energy of the arachidonic acid conformation in I805W:psWT, 39.947kcal/mol (Figure 6) is much greater than that in the psWT 8R-LOX, 2.706kcal/mol (Figure 4). This tells us that the arachidonic acid conformation docked into the mutant 8R-LOX is much less stable than that docked into the wild-type 8R-LOX. The structure of the manually defined binding pocket is changed by the mutation, causing the binding ability of 8R-LOX for arachidonic acid weakened, the binding of arachidonic acid to 8R-LOX not stable, and therefore the activity of the enzyme reduced. This consists with the experimental data that mutations significantly impact enzyme activity .
Productivity of arachidonic acid conformations
We analyzed the arachidonic acid conformations for all the docking simulations. We found that no docking simulation produces productive conformation considering the stereo- and region- specificity of the product and the process of catalysis. The 10th carbon atom in all arachidonic acid conformations does not position against the catalytic Fe atom (Figures 3, 4, 6). In the catalysis, the first step of the reaction is hydrogen abstraction at the 10th carbon of a pentadiene moiety of arachidonic acid by the activated Fe3+ to produce a free radical intermediate, which requires the 10th carbon to face to Fe atom. All our arachidonic acid conformations do not show this. The reason could be that we removed all the water molecules before doing the docking. However, in the hydrogen abstraction at the 10th carbon of arachidonic acid, a water molecule is required [1, 6]. We will perform the docking simulation with the water molecule present and hope to obtain productive conformations.
We applied the computational method “Internal Coordinate Mechanics” in docking arachidonic acid to 8R-LOX. The results showed that, for the wild type 8R-LOX, the conformation of arachidonic acid docked into ICM-identified docking site is more stable than that docked into the manually defined docking site. Mutation at I805 affects the structure of the putative binding site of 8R-LOX, which leads no docking pockets around Fe atom. The docking simulation in I805W:psWT 8R-LOX explains and consists with the experimental data presented by Neau . Further research and analysis is required to obtain productive conformations, identify the binding site of 8R-LOX and model the interactions between arachidonic acid and 8R-LOX.
Internal coordinate mechanics
- psWT- Pseudo-wild type 8R-LOX:
A deletion mutant of psWT 8R-LOX.
Brash AR: Lipoxygenases: occurrence, functions, catalysis, and acquisition of substrate. J Biol Chem 1999, 274: 23679–23682. 10.1074/jbc.274.34.23679
Kuhn H, Thiele BJ: The diversity of the lipoxygenase family. Many sequence data but little information on biological significance. FEBS Lett 1999, 449: 7–11. 10.1016/S0014-5793(99)00396-8
Dahlen SE, Hansson G, Hedqvist P, Bjorck T, Granstrom E, Dahlen B: Allergen challenge of lung tissue from asthmatics elicits bronchial contraction that correlates with the release of leukotrienes C4, D4, and E4. Proc Natl Acad Sci USA 1983, 80: 1712–1716. 10.1073/pnas.80.6.1712
Werz O, Steinhilber D: Therapeutic options for 5lipoxygenase inhibitors. Pharmacol Ther 2006, 112: 701–718. 10.1016/j.pharmthera.2006.05.009
Dahlen SE: Treatment of asthma with antileukotrienes: first line or last resort therapy? Eur J Pharmacol 2006, 533: 40–56. 10.1016/j.ejphar.2005.12.070
Liavonchanka A, Feussner I: Lipoxygenases: Occurrence, functions and catalysis. Journal of Plant Physiology 2006, 163: 348–357. 10.1016/j.jplph.2005.11.006
Neau DB, Gilbert NC, Bartlett SG, Boeglin W, Brash AR, Newcomer ME: The 1.85 Å Structure of an 8R-Lipoxygenase Suggests a General Model for Lipoxygenase Product Specificity. Biochemistry 2009, 48: 7906–7915. 10.1021/bi900084m
Ruben A, Maxim T, Dmitry K: ICM -A new method for protein modeling and design: Applications to docking and structure prediction from the distorted native conformation. Journal of Computational Chemistry 1994, 15: 488–506. 10.1002/jcc.540150503
Abagyan R, Totrov M: Biased Probability Monte Carlo Conformational Searches and Electrostatic Calculations for Peptides and Proteins. Journal of Molecular Biology 1994, 235: 983–1002. 10.1006/jmbi.1994.1052
Totrov M, Abagyan R: Flexible protein-ligand docking by global energy optimization in internal coordinates. Proteins 1997, 29(Suppl 1):215–220. Publisher Full Text 10.1002/(SICI)1097-0134(1997)1+%3C;215::AID-PROT29%3E;3.0.CO;2-Q
Stigler RD, Hoffmann B, Abagyan R, Schneider-Mergener J: Soft docking an and a peptide to an anticholera toxin antibody using internal coordinate mechanics. Structure 1999, 7: 663–670. 10.1016/S0969-2126(99)80087-2
Zhou Y, Abagyan R: How and why phosphotyrosine-containing peptides bind to the SH2 and PTB domains. Fold Des 1998, 3: 513–522. 10.1016/S1359-0278(98)00067-4
Farag NA, Mohamed SR, Soliman GAH: Design, synthesis, and docking studies of novel benzopyrone derivatives as H1-antihistaminic agents. Bioorganic & Medicinal Chemistry 2008, 16: 9009–9017. 10.1016/j.bmc.2008.08.039
Khan MTH, Orhan I, Senol FS, Kartal M, Sener B, Dvorsk M, Smejkal K, Slapetov T: Cholinesterase inhibitory activities of some flavonoid derivatives and chosen xanthone and their molecular docking studies. Chemico-Biological Interactions 2009, 181: 383–389. 10.1016/j.cbi.2009.06.024
Eissa AA, Farag NAH, Soliman GAH: Synthesis, biological evaluation and docking studies of novel benzopyranone congeners for their expected activity as anti-inflammatory, analgesic and antipyretic agents. Bioorganic & Medicinal Chemistry 2009, 17: 5059–5070. 10.1016/j.bmc.2009.05.073
Gani OABSM, Adekoya OA, Giurato L, Spyrakis F, Cozzini P, Guccione S, Winberg JO, Sylte I: Theoretical Calculations of the Catalytic Triad in Short-Chain Alcohol Dehydrogenases/Reductases. Biophysical Journal 2008, 94: 1412–1427. 10.1529/biophysj.107.111096
Abagyan R, Totrov M, Kuznetsov D: Icm - a New Method for Protein Modeling and Design - Applications to Docking and Structure Prediction from the Distorted Native Conformation. J Comput Chem 1994, 15: 488–506. 10.1002/jcc.540150503
Prasad JC, Goldstone JV, Camacho CJ, Vajda S, Stegeman JJ: Ensemble modelling of substrate binding to cytochromes P450: analysis of catalytic differences between CYP1A orthologs. Biochemistry 2007, 46: 2640–2654. 10.1021/bi062320m
RCSB Protein Data Bank.[http://www.pdb.org]
Hughes MA, Brash AR: Investigation of the mechanism of biosynthesis of 8 -hydroxyeicosatetraenoic acid in mouse skin. Biochim Biophys Acta 1991, 1081: 347–354.
The project described was supported by Grant Number P20RR016456 from the National Center For Research Resources. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Center For Research Resources or the National Institutes of Health.
We thank Dr. Marcia E. Newcomer and Dr. Konstantin G. Kousoulas for their valuable comments, suggestions, and mentorship.
This article has been published as part of BMC Bioinformatics Volume 11 Supplement 6, 2010: Proceedings of the Seventh Annual MCBIOS Conference. Bioinformatics: Systems, Biology, Informatics and Computation. The full contents of the supplement are available online at http://www.biomedcentral.com/1471-2105/11?issue=S6.
The authors declare that they have no competing interests.
SB conceived of the study, and participated in its design and coordination and drafted the manuscript. TD performed the docking simulation, helped to draft the manuscript. EK participated in the design. All authors read and approved the final manuscript.
About this article
Cite this article
Bai, S., Du, T. & Khosravi, E. Applying internal coordinate mechanics to model the interactions between 8R-lipoxygenase and its substrate. BMC Bioinformatics 11 (Suppl 6), S2 (2010). https://doi.org/10.1186/1471-2105-11-S6-S2
- Arachidonic Acid
- Docking Simulation
- Docking Site
- High Resolution Structure
- Interactive Docking