Journal of Biophysical Chemistry
Vol.08 No.01(2017), Article ID:76364,11 pages

Molecular Docking of the Inhibitory Activities of Triterpenoids of Lonchocarpus cyanescens against Ulcer

Isaiah A. Adejoro1, Sodiq O. Waheed1,2, Omolara O. Adeboye1,3*, Frederick U. Akhigbe1

1Department of Chemistry, University of Ibadan, Ibadan, Nigeria

2Department of Chemistry and Biochemistry, Faculty of Science, University of Porto, Porto, Portugal

3Department of Chemistry, Emmanuel Alayande College of Education, Oyo, Nigeria

Copyright © 2017 by authors and Scientific Research Publishing Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).

Received: December 15, 2016; Accepted: February 25, 2017; Published: February 28, 2017


Ulcer is one of the life threatening diseases. It is an open sore on an external or internal surface of the body caused by a break in the skin or mucous membrane which fails to heal. In this work, specific ligands that are suitable for ulcer have been studied computationally. Docking of the triterpenoids of Lonchocarpus cyanescens with target proteins of PDB codes 1AFC, 1AXM and 2AXM were performed using AutoDock Vina and Pymol for docking and post-docking analysis, respectively. In this study, the triterpenoid ligands with binding affinity/inhibitory constants −7.2/5.21, −6.5/16.99 and −6.2/28.20 for OH, and −6.7/12.12, −6.3/23.82 and −6.1/33.40 for OCH3 were compared with the standard ligands. Our study indicates that the results corresponding to triterpenoid ligands are close to standard ligands.


Ulcer, Molecular Docking, AutoDock Vina, Pymol, Binding Affinity

1. Introduction

Ulcers are holes that form in the lining of the upper part of the gastrointestinal tract (GIT). The different types of ulcer are peptic, gastric, and duodenal ulcers. A peptic ulcer is the umbrella term for all types of gastrointestinal ulcers. Technically, an ulcer is at least 0.5 centimeters wide, but can be much larger [1] . Duodenal ulcers mainly occur in people between 20 and 50 years old, and gastric ulcers typically appear in people over age 40 [2] [3] . Duodenal ulcers are about four times more common than gastric ulcers (Nematollahi et al., 2011). Ulcers can be caused by varieties of factors which include: Helicobacter pylori infection, excessive gastric acid secretion, the use of Non Steroidal Anti-Inflammatory Drugs (NSAIDs), heavy use of alcohol, smoking and decreased mucosal defence against gastric acid [4] [5] . Lonchocarpus cyanescens is a species of shrub from Fabaceae family. Various studies have been carried out on the bioactivity, phythotherapeutic and anti-psychotic properties [6] - [12] of Lonchocarpus cy- anescens. Experimental studies on Lonchocarpus cyanescens triterpenoids showed that they can be used to treat ulcer [9] . Bioassay investigation carried out on the leaf extract of Lonchocarpus cyanescens triterpenoids for anti oxidant activities revealed that the anti oxidant properties correlated with its phenolic and flavonoid contents [11] . Molecular docking studies of the triterpenoids (Figure 1) using three (3) different malaria targets with PDB codes 3QS1, 1LS5 and 1SME were investigated using AutoDock vina. The docking studies showed that the ligands docked well with the targets [12] . Study on the molecular mechanism of anti ulcer action of Lonchocarpus cyanescens triterpenoids has not been reported. In this study, Lonchocarpus cyanescens triterpenoids are used as the ligands in molecular docking to validate its inhibitory activity of against ulcer.

2. Computational Procedures

The initial structures of proteins which are responsible for ulcer (1AFC, I AXM and 2AXM), were obtained in the form of pdb files from the Protein Data Bank [13] (, and then to obtain the desired chains, multiple ligands and non-protein parts were deleted from the pdb files using Discovery Studio 4.1 visualizer. After that OpenBabel GUI version 2.3.2a and Spartan 14 version 1.18 were used to convert the pdb file format and optimize the geometry of the ligands, respectively. AutoDock Tools 1.5.6 and AutoDock Vina version 1.1.2

Figure 1. Structures of the ligands.

were used for molecular docking process [14] ( and to analyse the output of docking process, EduPymol version was used. The exhaustiveness which determines how comprehensive the software search for the best binding mode was set to the default value of 8 which accompanied the ligands from the list of protein residues reported in the literature to characterise their binding sites for all the docking runs. Autodock tools (ADT) was used to add polar hydrogen to the prepared receptor. The three-dimensional affinity and electrostatic grid boxes were generated to cover the entire active site of the receptor. The number of grid points in x, y, z, was set to 40, 40, 40 with the spacing 1.000 Å. The grid centre (x, y and z) was varied from one receptor to another (see, Table 1). After the successful completion of the docking runs, different conformations of the ligands known as binding modes were obtained with their respective binding affinity and the stable one which happens to be the one with the lowest binding affinity was picked as the pose and was employed in the post-docking analysis using Edupymol by opening the receptors in the form of pdbqt file as well as the different binding modes of the ligands obtained after running AutoDock vina which were also in pdbqt file format to check the best binding mode that fitted well with the binding site cavity. The protein-ligand hydrogen interactions were also obtained in terms of atomic distance. The inhibition constants Ki were calculated using Equation (1)


where ΔG is the binding affinity in kcal/mol, R is gas constant, 1.987 cal/mol/K and T is absolute temperature, assumed to be room temperature, 298.15 K [15] .

ADME and Drug Likeness Analysis

Molecular properties such as membrane permeability and bioavailability of leading compounds are always associated with some basic molecular descriptor such as logP (partition coefficient), molecular weight (mw), or number of hydrogen bond acceptors and donors in a molecule [16] . These molecular properties were used in formulating “rule of five” [17] . Lipinski’s rule states that molecule with good membrane permeability has MW ≤ 500, hydrogen bond donors ≤ 5 and hydrogen bond acceptor ≤ 10. Therefore, Lipinski’s Rule of Five was used to test the bioavailability characteristics such as Adsorption, Distribution, Metabolism and Elimination (ADME) of the ligands used in this study. In the present study, these molecular properties and druglikeness analysis were done using of Spartan 14 software.

Table 1. Grid centres (X, Y and Z) used for the docking runs.

3. Results and Discussion

Docking Score Results

From Table 2, the results of the studies revealed that the (the ligands used in this study) synthesized ligands compared favourably well with the standard ligands as shown in their binding affinities. The results also showed that OH derivative ligand has a better inhibitory effect (Figure 2(a), Figure 3(a) and Figure 4(a) for re-docking and Figure 2(b), Figure 2(c), Figure 3(b), Figure 3(c), Figure 4(b) and Figure 4(c) for OH and OCH3 respectively) than the OCH3 derivative which implies that OH derivative ligand will be more potent than OCH3 derivative in the treatment of ulcer. The OH and OCH3 derivatives triterpenoids gave binding affinities values (−7.2 Kcal/mol for 1AFC, −6.5 Kcal/mol for 1AXM and −6.2 Kcal/mol for 2AXM receptors) and (−6.7 Kcal/mol for 1AFC, −6.3 Kcal/ mol for 1AXM and −6.6 Kcal/mol for 2AXM receptors) respectively against the standard ligands value (−6.3 Kcal/mol for 1AFC, −6.2 Kcal/mol for 1AXM and −5.9 Kcal/mol for 2AXM receptors). Moreover, the smaller inhibition constant values (Ki) indicated that the ligands docked well to the receptor and that a relatively low concentration of the ligands is adequate to maximally occupy a ligand binding sites and trigger a physiological response. The binding affinity and inhibition constant for the OCH3 are higher than for OH. This is as a result of the inductive (electron releasing) effect of the methyl group. However, both derivatives have the capacity to inhibit ulcer, and as reported, this may be due to the high antioxidant properties of 14.83% at 0.50 mg/ml with scavenging effect of 63.89% and a strong correlation between antioxidant activity, phenol and flavonoid contents with correlation coefficient of 0.9906 and 0.9926 [18] . The findings from preliminary spectroscopic investigations that L. cyanescens extract contains essentially phyto-phenolic compounds was said to account for the use of the plants in the treatment of arthritic conditions and ulcers [6] [19] . Anti- ulcer & antioxidant activities of Hedranthera barteri [20] with possible involvement of H+, K+ ATPase inhibitory activity was assessed in cold-restraint (CRU), aspirin (ASP), alcohol (AL), pyloric ligation (PL) induced gastric ulcer models in rats and histamine-induced duodenal ulcer (HST) in guinea pigs. The effect of DMHBR (100 mg/kg) on gastric juice for free and total acidity, peptic activity and mucin secretion, using the pylorus ligated model, were evaluated. The H+, K+ and ATPase activity was assayed in gastric microsomes, spectrophotometri

Table 2. Docking score results for ulcer receptors.

Ei = Binding affinity (kcal/mol). Ki = Inhibition constant (µM).


Figure 2. (a) Re-docking Analysis of 1AFC Receptor (Ribbon and surface Representation). (b) Docking Analysis of 1AFC Receptor vs. OH derivative ligand (Ribbon and surface Representation). (c) Docking Analysis of 1AFC Receptor vs. OCH3 derivative ligand (Ribbon and Representation).

cally. The in vitro anti-oxidant assays were explored through DPPH, nitric oxide, hydroxyl radical, superoxide anion scavenging assays. It was reported that DMHBR reduced the incidence of ulcers in CRU (63.3%), PL (58.5%), ASP (52.7%), HST (75.0%) and AL (53.87%). Also, reductions were observed in the free acidity (49.4%), total acidity (45.8%) and peptic activity (32.9%) with increase in the mucin secretion by 81.6%. DMHBR (60 - 100 μg/ml) inhibited the H+, K+ and ATPase activity with IC50 of 89.64 μg/ml compared with omeprazole


Figure 3. (a) Re-docking Analysis of 1AXM Receptor (Ribbon and surface Representation). (b) Docking Analysis of 1AXM Receptor vs. OH derivative ligand (Ribbon and surface Representation). (c) Docking Analysis of 1AXM Receptor vs. OCH3 derivative ligand (Ribbon and surface Representation).

(10 - 50 μg/ml ) with IC50 of 32.26 μg/ml. DMHBR showed antioxidant activity with IC50 values of DPPH (397.69 μg/ml), nitric oxide (475.88 μg/ml), hydroxyl radical (244.22 μg/ml) and superoxide anion radical (285.20 μg/ml). It was observed that DMHBR showed anti-ulcer activity against experimentally-induced


Figure 4. (a) Re-docking Analysis of 2AXM Receptor (Ribbon and surface Representation). (b) Docking Analysis of 2AXM Receptor vs. OH derivative ligand (Ribbon and surface Representation). (c) Docking Analysis of 2AXM Receptor vs. OCH3 derivative ligand (Ribbon and Representations).

peptic ulcer models and exhibited both cytoprotective and anti-secretory property. It exhibited a proton pump inhibition activity and its anti-ulcer properties was said to be partly ascribed to its antioxidant activities [20] . The results obtained from the molecular docking of the inhibitory activity of L. cyanescens active ligands follow the same pattern with the results obtained from molecular docking of L. cyanescens ligands as malaria inhibitors [12] .

The molecular descriptor (see Table 2) such as molecular weight, polar surface area volume, hydrogen bond donor, hydrogen bond acceptor, polarizability, lipophilicity and lipophobicity are vital to describe the bioactivity of a molecule. The Lipinski properties such as molecular weight (358.610 for OH and 372.637 for OCH3 ligand derivative), partition coefficient (3.20 for OH and 3.42 for OCH3 ligand derivative) which indicates the lipophilicity of the ligand (LogP), an estimate of a compound’s overall lipophilicity, it is a factor that influences the behaviour of a compound in a range of biological membranes, hepatetic clearance, selectivity and non-specific toxicity [21] . For good absorption, log P of a compound must be ≤5 therefore, these compounds are effective in terms of lipophilicity. Based on the Lipinski’s rule, the drug must have molecular weight value ≤ 500, Hydrogen Bond Donor ≤ 5, Hydrogen Bond Acceptor ≤ 10 and partition coefficient (LogP) value ≤ 5. From the results in Table 2, the ligands passed the Lipinski’s rule of five. Also, the Polar Surface Area (PSA) which is an indicator of the ligand hydrophilicity plays an important role in shaping the protein-ligand interaction by affecting the non-bonded contribution to the binding energy. Molecules with PSA greater than 140 Å2 are usually believe to be poor at penetrating cell membranes. For molecules to penetrate the blood-brain barrier, PSA should be less than 60 Å2 [16] [22] which is also fulfilled. The number of Hydrogen Bond Donor (1 for OH and 0 for OCH3 ligand derivative) and Hydrogen Bond Acceptor (1 for OH and 1 for OCH3 ligand derivative) obtained are also in good agreement with the Lipinski value for HBD and HBA (Table 2). The HOMO and LUMO are important properties that offer realistic and qualitative information about the excitation properties of molecules [23] [24] . The calculated HOMO are −6.53 eV for OH, and −6.18 eV for OCH3, and the calculated LUMO are −4.76 eV for OH, and −4.34 eV for OCH3 (Table 3). Therefore, the calculated band gaps are 1.77 eV for OH and 1.84 eV for OCH3. Band gap plays important role in the protein-ligand interaction, the lower the band gap, the easier the excitation and the better the ability of a molecule to donate electrons to the surrounding. The calculated ovality which measure the degree of deviation [24] of the cross section of the core from prefect circularity are 1.42 and 1.46 for OH and OCH3 derivatives, respectively. The calculated dipole moments are 1.18 and 0.92 for OH and OCH3, respectively. The dipole moment for OH derivative shows that the non-bonded interactions, the dipole-dipole interactions are essential in ligand-receptor interaction which has been shown to contribute about 3 to 5 kJ/mol [25] to the ligand-receptor energy.

Furthermore, Asn-18, Gln-127, Lys-118, Arg-122, Lys-113, Asn-114 residues are involved in the binding of the ligands to the 1AFC receptor. Similarly, Arg-122, Lys-118, Gln-127, Ala-129, Asn-114, Lys-128 residues are also involved in the binding of the ligands to 1AXM receptor while Arg-122, Lys-118, Leu-111, Lys-113, Gly-126, Arg-119 and Gln-127 residues are involved in the binding of the ligands to the 2AXM receptor binding cavity. The Hydrogen bond distance between the ligands and the proteins are shown in Table 4.

Table 3. Properties of the ligands used.

Table 4. Interaction between ligands and the receptors.

4. Conclusion

This study shows that Lonchocarpus cyanescens triterpenoids as the ligands in molecular docking to validate its inhibitory activity of against ulcer are effective and can be used to design drugs to treat ulcer because the results of the triterpenoid ligands are in good agreement with the standard ligands.

Cite this paper

Adejoro, I.A., Waheed, S.O., Adeboye, O.O. and Akhigbe, F.U. (2017) Molecular Docking of the Inhibitory Activities of Triterpenoids of Lonchocarpus cyanescens against Ulcer. Journal of Biophysical Chemistry, 8, 1-11.


  1. 1. Nematollahi, A., Noushin, A. and Wiart, C. (2011) Design and Modelling Studies on Liriodenine Derivatives as Novel Topoisomerase II Inhibitors. International Journal of ChemTech Research, 3, 1622-1627.

  2. 2. Oliveira, M., Cenzi, G., Nunes, R.R., Andrighetti, C.R., Valadao, D.M., Reis, C., Simoes, C.M., Nunes, R.J., Junior, M.C., Taranto, A.G., Sanchez, B.A., Viana, G.H. and Varotti, F.P. (2013) Antimicorbial Activity of 4-Metoxychalcone: Docking Studies as Falcipain/Plasmepsin Inhibitors, ADMET and Lipophilic Efficiency Analysis to Identify a Putative Oral Lead Candidate. Molecules, 18, 15276-15287.

  3. 3. Coutinho, J., Cortopassi, W., Oliveira, A.A., Franca, T.C. and Krettli, A.U. (2011) Antimalarial Activity of Potential Inhibitors of Plasmodium Falciparum Lactate Dehydrogenase Enzyme Selected for Docking Studies. PLoS ONE, 6, 1-7.

  4. 4. Das, B., Mondal, U., Sharma, S., Ghosh, P. and Bethra, A.K. (2012) Molecular Docking of Triazine Analogues. Journal of Chemical and Pharmaceutical Research, 4, 1595-1600.

  5. 5. Kaladhar, D., Satyanarayana, K., Chaitanya, A. and Hussain, S. (2010) Clinical Analysis Drug Designing and QSAR Studies on Rheumatoid Arthritis. International Journal of Pharma and Bio Sciences, 1, 132-138.

  6. 6. Iwu, M.M. and Anyanwu, B.N. (1982) Phytotherapeutic Profile of Nigerian Herbs. I: Anti Inflammatory and Anti-Arthritic Agents. Journal of Ethnopharmacology, 6, 263-274.

  7. 7. Shuaibu, M.N., Wuyep, P.A., Yanagi, T., Hirayama, K., Tanaka, T., et al. (2008) The Use of Microfluorometric Method for Activity-Guided Isolation of Antiplasmodial Compound from Plant Extracts. Parasitology Research, 102, 1119-1127.

  8. 8. Sonibare, M.A., Umukoro, S. and Shonibare, E.T. (2012) Antipsychotic Property of Aqueous and Ethanolic Extracts of Lonchocarpus cyanescens (Schumach and Thonn.) Benth. (Fabaceae) in Rodents. Journal of Natural Medicines, 66, 127-132.

  9. 9. Moronkola, D.O. and Oladosu, I.A. (2013) Chemical Composition of Lonchocarpus cyanescens Benth., (Fabaceae)—Case Study of Its Volatile Oils, and Two Triterpenoids. American Journal of Plant Sciences, 4, 1653-1659.

  10. 10. Arowona, I.T., Sonibare, M.A. and Umukoro, S. (2014) Antipsychotic Property of Solvent-Partitioned Fractions of Lonchocarpus cyanescens Leaf Extract in Mice. Journal of Basic and Clinical Physiology and Pharmacology, 25, 235-240.

  11. 11. Samuel, B., Adigun, O. and Adaramoye, O. (2014) Bioassay-Guided Investigation of Lonchocarpus cyanescens Benth Leaves Extracts for Antioxidant Activities. Academic Journals, 13, 2240-2247.

  12. 12. Adejoro, I.A., Waheed, S.O. and Adeboye, O.O. (2016) Molecular Docking Studies of Lonchocarpus cyanescens Triterpenoids as Inhibitors for Malaria. Journal of Physical Chemistry & Biophysics, 6, 213.

  13. 13. Protein Data Bank (2015)

  14. 14. Trott, O. and Olson, A.J. (2010) AutoDock Vina: Improving the Speed and Accuracy of Docking with a New Scoring Function, Efficient Optimization and Multithreading. Journal of Computational Chemistry, 31, 455-461.

  15. 15. Atkins, P.W. (1986) Physical Chemistry. Oxford University Press, Oxford, 263.

  16. 16. Ertl, P., Rohde, B. and Selzer, P. (2000) Fast Calculation of Molecular Polar Surface Surface Area as a Sum of Fragment-Based Contributions and Its Application to the Prediction of Drug Transport Properties. Journal of Medicinal Chemistry, 43, 3714-3717.

  17. 17. Lipinski, C.A., Lombardo, F., Dominy, B.W. and Feeney, P.J. (2001) Experimental and Computational Approaches to Estimate Solubility and Permeability in Drug Discovery and Development Settings. Advanced Drug Delivery Reviews, 46, 3-26.

  18. 18. Babatunde, S., Adigun, O. and Adaramoye, O. (2014) Bioassay-Guided Investigation of Lonchocarpus cyanescens Benth Leaves for Antioxidant Activities. African Journal of Biothechnology, 13, 2240-2247.

  19. 19. Onasanwo, S.A., Singh, N., Olaleye, S.B., Mishra, V. and Palit, G. (2010) Anti-Ulcer & Antioxidant Activities of Hedranthera barteri {(Hook F.) Pichon} with Possible Involvement of H+, K+ ATPase Inhibitory Activity. Indian Journal of Medical Research, 132, 442-449.

  20. 20. Hughes, J.D., Blagg, J., Price, D.A., Bailey, S., Decrescenzo, G.A., Devraj, R.V., Ellsworth, E., Obian, Y.M., Gibbs, M.E., Gilles, R.W., Greene, N., Huang, E., Krieger-Burke, T., Loesel, J., Wager, T., Whiteley, L. and Zhang, Y. (2008) Physiochemical Drug Properties Associated with In Vivo Toxicological Outcomes. Bioorganic & Medicinal Chemistry Letters, 18, 4872.

  21. 21. Singh, Y.P. and Singh, R.A. (2008) Theoretical Studies of Different Tautomers of Anti Cancer Drug: Dichloroacetate. Pakistan Journal of Pharmaceutical Science, 21, 390-395.

  22. 22. Bouachrine, M., Hamidi, M., Bouzzine, S.M. and Taoufik, H. (2009) Theoretical Study on the Structure and Electronic Properties of New Materials Based on Thiophene and Oxadiazole. Journal of Chemical Research, 101, 29-37.

  23. 23. Semire, B., Oyebamiji, A. and Ahmad, M. (2012) Theoretical Study on Structure and Electronic Properties of 2,5-Bis[4-N,N-Diethylaminostyryl] Thiophene and Its Furan and Pyrrole Derivatives Using Density Functional Theory (Dft). Pakistan Journal of Chemistry, 2, 166-173.

  24. 24. Leach, A.R. (2001) Molecular Modelling: Principle and Applications. Prentice Hall, Upper Saddle River.

  25. 25. David, F.L.L. and Howard, B.B. (2002) Molecular Binding Interactions, Their Estimation and Rationalization in QSARs in Terms of Theoretically Derived Parameters. The Scientific World Journal, 2, 1776-1802.