Skip to main content
Download PDF

Volume 11 Supplement 6

Proceedings of the Seventh Annual MCBIOS Conference. Bioinformatics: Systems, Biology, Informatics and Computation

Applying internal coordinate mechanics to model the interactions between 8R-lipoxygenase and its substrate

BMC Bioinformatics volume 11, Article number: S2 (2010) Cite this article

Abstract

Background

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.

Methods

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.

Results

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.

Conclusions

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.

Background

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 [3]. 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 [6]. An understanding of substrate arachidonic acid recognition is helpful in the design of enzyme-specific inhibitors [7]. 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 Å [7]. 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.

Methods

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 [19]. 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).

Figure 1
figure 1

The 1.85Å resolution crystal structure of pseudo-wild type 8R-LOX and the 1.9Å resolution crystal structure of the mutant I805W:psWT 8R-LOX used in docking simulation. The structure of I805W:psWT superimposes on the structure of psWT 8R-LOX. Their structure is nearly identical except that in the area of the mutation. W805 is marked out in the figure. The catalytic Fe atom is in blue.

Full size image

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 [7]. This manually defined docking site includes 13 residues (Table 1).

Table 1 Residues of the ICM-identified docking site for psWT and residues of the manually defined docking site for psWT and I805W:psWT 8R-LOX
Full size table

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.

Results

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).

Figure 2
figure 2

Distribution of the 26 potential ligand binding pockets of psWT 8R-LOX. The pockets are marked in different colours. The pocket for simulating the substrate docking is marked in pink. The catalytic Fe atom is in blue.

Full size image

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.

Figure 3
figure 3

The conformation of arachidonic acid docked into the ICM-identified docking site of psWT 8R-LOX with lowest energy. Upper: An overall view of the conformation of arachidonic acid. Lower: A close up view of the conformation of acrachidonic acid. The arachidonic acid is in yellow, the carboxyl terminal of it is in red. The catalytic Fe atom is in blue. The energy for this conformation is -63.623kcal/mol.

Full size image

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.

Figure 4
figure 4

The conformation of arachidonic acid docked into the manually defined docking site of psWT 8R-LOX. Upper: An overall view of the conformation of arachidonic acid. Lower: A close up view of the conformation of acrachidonic acid. The arachidonic acid is in yellow, the carboxyl terminal of it is in red. The catalytic Fe atom is in blue. The energy for this conformation is 2.706kcal/mol.

Full size image

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.

Figure 5
figure 5

The 25 potential ligand binding pockets of I805WT:psWT 8R-LOX. Pockets are patches in various colours. No pocket is located around the catalytic Fe atom. The catalytic Fe atom is in blue.

Full size image

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).

Figure 6
figure 6

The conformation of arachidonic acid docked into the manually defined docking site of I805WT:psWT 8R-LOX. Upper: An overall view of the conformation of arachidonic acid. Lower: A close up view of the conformation of acrachidonic acid.The arachidonic acid is in yellow, the carboxyl terminal of it is in red. The catalytic Fe atom is in blue. The energy for this conformation is 39.947kcal/mol.

Full size image

Discussion

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 [7] 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 [7].

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+[14] 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.

Conclusions

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 [7]. 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.

Abbreviations

LOX:

Lipoxygenase

ICM:

Internal coordinate mechanics

psWT- Pseudo-wild type 8R-LOX:

I805W

psWT:

A deletion mutant of psWT 8R-LOX.

References

  1. Brash AR: Lipoxygenases: occurrence, functions, catalysis, and acquisition of substrate. J Biol Chem 1999, 274: 23679–23682. 10.1074/jbc.274.34.23679

    Article  CAS  PubMed  Google Scholar 

  2. 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

    Article  CAS  PubMed  Google Scholar 

  3. 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

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  4. Werz O, Steinhilber D: Therapeutic options for 5lipoxygenase inhibitors. Pharmacol Ther 2006, 112: 701–718. 10.1016/j.pharmthera.2006.05.009

    Article  CAS  PubMed  Google Scholar 

  5. 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

    Article  CAS  PubMed  Google Scholar 

  6. Liavonchanka A, Feussner I: Lipoxygenases: Occurrence, functions and catalysis. Journal of Plant Physiology 2006, 163: 348–357. 10.1016/j.jplph.2005.11.006

    Article  CAS  PubMed  Google Scholar 

  7. 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

    Article  CAS  PubMed  Google Scholar 

  8. 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

    Article  Google Scholar 

  9. 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

    Article  CAS  PubMed  Google Scholar 

  10. 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

    Article  Google Scholar 

  11. 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

    Article  CAS  PubMed  Google Scholar 

  12. 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

    Article  CAS  PubMed  Google Scholar 

  13. 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

    Article  CAS  Google Scholar 

  14. 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

    Article  CAS  PubMed  Google Scholar 

  15. 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

    Article  CAS  Google Scholar 

  16. 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

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  17. 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

    Article  CAS  Google Scholar 

  18. 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

    Article  CAS  PubMed  Google Scholar 

  19. RCSB Protein Data Bank.[http://www.pdb.org]

  20. Hughes MA, Brash AR: Investigation of the mechanism of biosynthesis of 8 -hydroxyeicosatetraenoic acid in mouse skin. Biochim Biophys Acta 1991, 1081: 347–354.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

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.

Author information

Authors and Affiliations

  1. Department of Computer Science, Southern University and A&M College, Baton Rouge, LA, 70813, USA

    Shuju Bai, Tianchuan Du & Ebrahim Khosravi

Authors
  1. Shuju Bai
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Tianchuan Du
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Ebrahim Khosravi
    View author publications

    You can also search for this author in PubMed Google Scholar

Corresponding author

Correspondence to Shuju Bai.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

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.

Rights and permissions

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.

Reprints and permissions

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

Download citation

  • Published:

  • DOI: https://doi.org/10.1186/1471-2105-11-S6-S2

Keywords

BMC Bioinformatics

ISSN: 1471-2105

Contact us