Skip to main content
  • Research article
  • Open access
  • Published:

Cheminformatics studies to analyze the therapeutic potential of phytochemicals from Rhazya stricta

Abstract

Rhazya stricta is a unique medicinal plant source for many indole alkaloids, non-alkaloids, flavonoids, triterpenes and other unknown molecules with tremendous potential for therapeutic applications against many diseases. In the present article, we generated computational data on predictive properties and activity across two key therapeutic areas of cancer and obesity, and corresponding cheminformatics studies were carried out to examine druggable properties of these alkaloids. Computed physiochemical properties of the 78 indole alkaloids from R. stricta plant using industry-standard scientific molecular modeling software and their predictive anti-cancer activities from reliable web-source technologies indicate their plausible therapeutic applications. Their predictive ADME properties are further indicative of their drug-like-ness. We believe that the top-ranked molecules with anti-cancer activity are clearly amenable to chemical modifications for creating potent, safe and efficacious compounds with the feasibility of generating new chemical entities after pre-clinical and clinical studies.

Background

Rhazya stricta Decsne (Apocynaceae family), a traditional herbal medicinal plant from Western and South Asia, has been shown to have multiple pharmacological effects due to the presence of over 100 alkaloids [1,2,3]. The chemical constituents of this plant (R. stricta) may possess biological activities of antifungal, antimicrobial, antioxidant, CNS, hypertension, metabolic, and inflammatory disorders. Rhazimine, an alkaloid isolated from R. stricta leaves, was shown to affect arachidonic acid metabolism in human blood [4]. This alkaloid was shown to be a dual and selective inhibitor of platelet activating factor (PAF)-induced platelet aggregation and arachidonic acid metabolism. Other effects of the lyophilized extract of R. stricta include an antispasmodic effect in rat muscles [5]. In another study, antioxidant effects were observed at higher doses, and it reduced the hepatic and renal concentrations of glutathione (GSH) and increased the ascorbic acid levels, whereas the degree of lipid peroxidation was reduced [6]. A recent study has shown that the basic alkaloid fraction from R. stricta significantly induces one of the chemopreventive enzyme-Nqo 1, through an Nrf 2-dependent mechanism, thereby establishing its role as an anti-tumor agent [7]. In another pharmacological study, the biochemical parameters including blood lipid profile concentrations, liver enzyme activities and kidney functions were analyzed in rats [8]. It was also found that aqueous extract of R. stricta and indole alkaloids caused a significant increase in serum adiponectin levels and resulted in significant improvements in insulin resistance [9]. In another follow up study, we observed indole-alkaloids of R. stricta improved not only the lipid profile and liver function but also led to improvements in the insulin levels in rats, most likely via modulating insulin resistance [10]. Indole-alkaloids of R. stricta had been reported to have anticancer properties [11]. Other studies by our departmental colleagues showed that alkaloid extract of R. stricta leaves inhibited proliferation, colony formation and anchorage-independent growth in various cancer cell lines such as colon cancer, breast cancer and lung cancer [12,13,14].

Understanding the chemical structure, physiochemical, and chemical-informatic properties of these natural product compounds will give clues for further modifications required in their structures responsible for their biological activities. Even though, there have been about 100 chemical entities of indole-based alkaloid constituents of R. stricta which have been reported but their chemical structures are yet to be clustered and identified, and moreover the pharmacological application of any one of these constituents towards human health is yet to be identified. Understanding qualitative correlation of structures to their chemical druggability, IP potential, and their applicability towards a therapeutic area would be worth exploring prior to pre-clinical studies. Availability of this plant (R. stricta), thus its phytochemical constituents largely in Arabian and South Asian region makes it worth studying through computational, synthetic, and biological view point. Indole based alkaloids such as vinblastine and vincristine are well known for their anti-cancer properties. From systematically generated informatics data analysis, one would be able to evaluate the physiochemical properties of the potential therapeutic compounds. These promising molecules which have “desirable pharmacophores” may provide obvious extension to a better targeted therapeutic benefit. Conventional drugs obey set of rules such as Lipinski’s Rule-of-Five (RO5) [15], wherein all orally administered molecules need to have certain physiochemical properties. Calculation of these cheminformatic properties has thus become essential for all projects of new drug discovery which go through oral route of administration. Along with RO5, the new molecules also have to adhere to certain parameters which yield favorable ADMET outcome of an oral drug. We further evaluated these molecules for therapeutic activity, including anticancer, anti-obesity, anti-inflammatory, and anti-bacterial properties. Although these predictions are indicative only, the value of predictions in various target classes and therapeutic areas would be very useful for future experimental studies. Moreover, their metabolic fate with key enzymes such as P450’s is also predicted for probable drug–drug and drug-target (P450) interactions (reviewed in [16, 17]).

Methods

For prediction of various therapeutic potential of these molecules, commercially and publicly available technologies as below were utilized.

  1. a.

    PharmaExpert (http://www.pharmaexpert.ru)—PASS [18]

  2. b.

    Superpred (http://prediction.charite.de)—Predictive Targets [19]

  3. c.

    SwissTargetPrediction (http://www.swisstargetprediction.ch)—Predictive Target [20]

  4. d.

    CDRUG (http://bsb.kiz.ac.cn/CDRUG)—Anti-cancer activity [21]

Schrodinger [22], a scientific software that predicts drug-like properties and liabilities (viz. HERG and CNS), and ACD/Labs [23] for physiochemical and cheminformatics studies were utilized. Details of the molecules, names, structures were obtained from the literature, commercial sources, and knowledge-based web sources. Tables 1 and 2 gives the details of these molecules together with their 2D SMILES notation, respectively.

Table 1 Chemical structures and names of Rhazya stricta compounds
Table 2 SMILES codes for Rhazya stricta compounds

Results and discussion

Physiochemical and cheminformatic studies

ACD/Laboratories informatics modules generated physiochemical and cheminformatics data of R. stricta indole and non-indole alkaloids. For all the selected 78 molecules in this study, it was observed that less than 20% of the molecules are having molecular weights >450, while most molecules range around 300–350, indicating their viability for additional medicinal chemistry amenable nature. Most of these molecules are also moderately to highly soluble—mainly due to the high value of pKa (leading to solubility at neutral pH). Additionally, many of these indole/non-indole molecules are also less lipophilic (~75% of them have logP ~3 to 4). Alkaloids that violate Lipinski’s Rule-of-5 are either due to molecular weight or logP, are tetrahydrosecamine; presecamine; beta-sitosterol; ursolic acid; stigmasterol; oleanolic acid; secamine; bis-strictidine; 3,14-dehydrorhazigine; 16-hydroxyrhazisidine; rhazisidine; rhazigine; dihydrosecamine; dihydropresecamine; tetrahydropresecamine; decarbomethoxy-15,17-tetrahydrosecodine;16s,16′-decarboxytetrahydro-secamine. Figures 1 and 2 give the plots of molecular weight and LogP (lipophilicity) of individual compounds, accordingly. Since most of the molecules have a basic nitrogen and sometimes, may be more than one, leading to a larger pKa at physiological pH—thus leading most molecules are highly to moderately soluble at physiological pH. Very few compounds and non-indole alkaloids have no basic nitrogen leading to highly insoluble compounds in water at physiological pH. As the acidity goes up (leading towards pH 1), most compounds become largely soluble. A qualitative and quantitative (computational) estimate of solubility of these compounds are given in Tables 3 and 4, respectively.

Fig. 1
figure 1

Variation of Molecular weight of compounds of Rhazya stricta

Fig. 2
figure 2

Variation of LogP of compounds of Rhazya stricta

Table 3 Qualitative assessment of Rhazya stricta compounds with respect to Lipinski’s Rule-of-5 and solubility
Table 4 Predicted solubility and pKa (acid and base) of various Rhazya stricta compounds

QUIKPROP calculations

Predicted Quikprop properties for potential cardiac liabilities such as HERG, and CNS liabilities (Blood–Brain-Barrier) and drug-like nature of these molecules indicate that many of these molecules are well within the boundaries of accepted hit-, and lead-like nature. QuikProp calculations were performed using Schrodinger’s Maestro for various alkaloids of R. stricta. These predictions not only give Rule-of-5 data, but also predict the cardiotoxicity predictions (HERG) and CNS penetration potential (logBBB) properties. More importantly, it also gives the prediction regarding cell-permeability (Caco2). All these models are well validated in literature, and most of them perform well within the reproducible results for training datasets. Results indicate that many of the molecules have decent permeation through Caco2 cell lines (>300), while the polar surface area (PSA) is not too high (>120) for oral absorption. For HERG toxicity prediction, below −5 (i.e. −6, −7 etc.) is not considered to be safe. Hence, those molecules whose logHERG values are well below -5 (such as geissoschizine, presecamine, tetrahydrosecamine) may exhibit cardioliability. The human intestinal absorption is also predicted, and it appears for most molecules, these values are larger. Any %HIA prediction >90% is expected to be well absorbed, and their polar surface area (PSA) is also a direct correlation to it. Those molecules whose molecular weights are >500 exhibit rule-of-5 violation and this violation goes beyond 1 to a maximum of 3. Those molecules appear structurally much larger and like dimers. Table 5 gives computed Quikprop computed values of various alkaloids of R. stricta. Table 6 also indicates various other physiochemical parameters including surface tension, parachor etc. of R. stricta indole and non-indole analogs.

Table 5 Quikprop calculation (for physiochemical properties) of Rhazya stricta compounds
Table 6 Surface related and ring-related properties of Rhazya stricta compounds

Predicted therapeutic area applications

PASS—prediction of activity spectra for substances

This web-based predictive server from Way2Drug, has variety of annotators of substances for their probability of active or inactive towards few targets. Out of all services and products of them, we utilized PASS method of predictions. More than 100 activities are predicted with their probability of activities and in-activities. Some of them include kinase inhibitors, GPCR antagonists, and some specific targets like adrenergic receptors, and their kinase inhibitors. We considered the probability of active (Pa) >0.3 (i.e. >30%), and should be greater than probability of inactive (Pi). Given these conditions, we observed many alkaloids have indicated Pa >0.8 in certain conditions (such as, anthrine has predicted Pa at 90% towards β-adrenergic receptor kinase inhibitor, 5-HTA release stimulant). Majority of them also is predicted to be substrate to CYP3A4 and CYP2D6 indicating their metabolic instability (Pa ~ 0.5, 0.4, respectively). Several such predictions for all 78 alkaloids has been computed—leaving predictions to be validated, experimentally. Similarly, dihydrocorynantheol and corynantheol were also predicted to be 5-HT release stimulants, and have been projected to be chemosensitizers. Eburnamenine is predicted to be a Nootropic agent at 90% Pa, while eburnamine is predicted to be a CNS (anti-depressant and mood disorder management agent at >96% Pa). Strictosidine is predicted to be an antiprotozoal at 86% Pa, β-sitosterol is anti-hypercholesterolemic agent with Pa ~98%, rhazidigenine (rhazidine) is an antidyskinetic at 60% Pa, secamine is a H1F1A expression inhibitor at 83% Pa (but a non-pharmaceutically acceptable molecule due to high MW and many RO5 violations). A similar observations is also made for 16-hydrorhazisidine (72% Pa for H1F1A expression inhibitor). Strictamine is predicted to be gluconate 2-dehydrogenase acceptor with 70% Pa, and 1,2-dehydroaspidospermine (which is a small molecule) has been predicted to be analeptic with 77% Pa. Dihydrosecamine is predicted to be a H1F1A expression inhibitor with 77% Pa, and rhazidigenine-N-oxide is predicted to be a cognition disorder agent with 64% Pa. Decarbomethoxy-15,20,16,17-tetrahydrosecodine is a small molecule with ~70% Pa for antidyskinetic and antineuronic agent, 1,2-dehydrospidospermidine-N-oxide is predicted to be 87% as analeptic.

Anticancer activity through CDRUG

This set of predictions using the structures and SMILES codes of the alkaloids, annotates the anti-cancer activity by predicting “Mean logGI50”. Most molecules that have Mean LogGI50 values lower than −5 are considered to have anti-cancer activity. It is interesting to know that all the molecules of R. stricta alkaloids (indole/non-indole) have predicted mean logGI50 values ranging between −4.95 and −6.50—indicating they all may have anti-cancer activities. There are about 10 compounds that have predicted logGI50 values less than −6, which indicate strong anti-cancer activity. Table 7 shows the predicted mean LogGI50 values of all the compounds considered in the present study.

Table 7 Predicted mean LogGI50 of Rhazya stricta compounds whose values lower than −6.0 are highlighted in italics may exhibit anti-cancer activity

SuperPred—predicted target interactions

From this server studies on R. stricta alkaloids, we observed that many of these molecules may interact with CYP2D6 or CYP3A4 as substrates. The indication of these results mean that their target may be unknown, but they do modify the drug metabolism, and affect drug–drug interactions.

SwissTarget prediction

While predictions from this web-server may suggest each molecule have certain target activity, they almost correlate well with the PASS server prediction—which gives additional probability of prediction for each molecule to be active or inactive against the target of interest.

Overall from the calculated cheminformatics studies and web-server predictions, we understand that few molecules like anthrine, condylocarpine, dihydrocorynantheol etc. have predicted GIC50 values in sub µM concentrations, while they also have predicted drug–drug activity towards CYP3A4, and CYP2D6 enzymes. Most molecules turnout to be modulators of membrane receptor ligands while some have predicted cholinesterase, CNS (5HT2x), adenosine (A2A/A2B) activity. Moreover, all molecules have predicted activity towards certain targets (Pa > 30%).

Conclusions

Table 8 indicates the top 10-best naturally occurring indole alkaloids of R. stricta that were predicted to be having decent anti-cancer activity and other good physiochemical properties together with cheminformatics properties—these molecules are antirhine, 3-epi-antirhine, condylocarpine, eburnamine, eburnamonine, taberonine, ursolic acid, stigmasterol, olenaolic acid, (+)-vincadiformine, (−)-vincadiformine, (−)-16R,21R-omethyleburmanine, 15-hydroxy-vincadifformine, and dihydroburnamenine.

Table 8 Key details of top molecules with predicted targets for anti-cancer and anti-obesity, probable rule-of-5, predicted LogGI50 with predicted H-, and p values

References

  1. Marwat SK, Fazal-ur-Rehman, Usman K, Shah SS, Anwar N, Ullah I (2012) A review of phytochemistry, bioactivities and ethno medicinal uses of Rhazya stricta Decsne (Apocynaceae). Afr J Microbiol Res 6(8):1629–1641

    CAS  Google Scholar 

  2. Ahmad Y, Fatima K, Le Quesne PPW, Atta-ur-Rahman (1983) Further alkaloidal constituents of the leaves of Rhazya stricta. Phytochemistry 22:1017–1019

    Article  CAS  Google Scholar 

  3. Atta-ur-Rahman, Zaman K, Habib-ur-Rehman, Malik S (1986) Studies on alkaloids of Rhazya stricta. J Nat Prod 49:1138–1139

    Article  CAS  Google Scholar 

  4. Saeed SA, Simjee RU, Mahmood F, Sultana N, Atta-ur-Rahman (1993) Rhazimine from Rhazya stricta: a dual inhibitor of arachidonic acid metabolism and platelet activating factor-induced platelet aggregation. Planta Med 59(6):566–568

    Article  CAS  Google Scholar 

  5. Tanira MO, Ali BH, Bashir AK, Wasfi IA, Chandranath I (1996) Evaluation of the relaxant activity of some United Arab Emirates plants on intestinal smooth muscles. J Pharm Pharmacol 47:457–459

    Google Scholar 

  6. Ali BH, Bashir AK, Tanira MO (2000) The anti-oxidant action of Rhazya stricta in rats. Phytothera Res 14:469–471

    Article  CAS  Google Scholar 

  7. El Gendy MAM, Ali BH, Michail K, Siraki AG, El Kadi AOS (2012) Induction of quinine oxidoreductase 1 enzyme by R. stricta through Nrf2-dependent mechanisms. J Ethnopharmacol 144:416–424

    Article  Google Scholar 

  8. Baeshen NA, Lari SA, Aldoghaither HA, Elkady AI (2010) Biochemical evaluation of the effect of Rhazya stricta aqueous leaves extract in the liver and kidney functions in rats. Nat Sci 8:136–142

    Google Scholar 

  9. Baeshen N, Lari S, Al Doghaither HA, Ramadan HAI (2010) Effect of Rhazya stricta extract on rat adiponectin gene and insulin resistance. J American Sci 6:1237–1245

    Google Scholar 

  10. Baeshina NA, Yaghmoor SS, Ashmaouia HM, Kumosani TA, Saini KS (2014) The indole-alkaloid fraction of Rhazya stricta improves key biochemical parameters associated with metabolic syndrome in rats. Bothalia J 44:358–371

    Google Scholar 

  11. Mukhopadhayay S, Handy GA, Funayama S, Cordell GA (1981) Anticancer indole alkaloids of Rhazya stricta. J Nat Prod 44:696–700

    Article  Google Scholar 

  12. Baeshen NA, Elkady AI, AbuZinadah OA, Mutwakil MH (2012) Potential anticancer activity of the medicinal herb, Rhazya stricta, against human breast cancer. Af J Biotechnol 11(37):8960–8972

    Article  Google Scholar 

  13. Elkady AI (2013) Crude alkaloid extract of Rhazya stricta inhibits cell growth and sensitizes human lung cancer cells to cisplatin through induction of apoptosis. Gen Mol Biol 36(1):12–21

    Article  Google Scholar 

  14. Elkady AI, Hussein RAH, Abu-Zinadah OA (2014) Differential control of growth, apoptotic activity and gene expression in human colon cancer cells by extracts derived from medicinal herbs, Rhazya stricta and Zingiber officinale and their combination. World J Gastroenterol 20(41):15275–15288

    Article  Google Scholar 

  15. Lipinski CA, Lombardo F, Dominy BW, Feeney PJ (2001) Experimental and computational approaches to estimate solubility and permeability in drug discovery and development settings. Adv Drug Deliv Rev 46(1–3):3–26

    Article  CAS  Google Scholar 

  16. Purnapatre K, Khattar SK, Saini KS (2008) Cytochrome P450s in the development of target-based anticancer drugs. Cancer Lett 259:1–15

    Article  CAS  Google Scholar 

  17. Singh D, Kashyap A, Pandey RV, Saini KS (2011) Novel advances in cytochrome P450 research. Drug Discovery Today 16:793–799

    Article  CAS  Google Scholar 

  18. Filimonov DA, Lagunin AA, Gloriozova TA, Rudik AV, Druzhilovskii DS, Pogodin PV, Poroikov VV (2014) Prediction of the biological activity spectra of organic compounds using the PASS Online web resource. Chem Heterocycl Compd 50(3):444–457

    Article  CAS  Google Scholar 

  19. Nickel J, Gohlke BO, Erehman J, Banerjee P, Rong WW, Goede A, Dunkel M, Preissner R (2014) SuperPred: update on drug classification and target prediction. Nucleic Acids Res 42:26–31

    Article  Google Scholar 

  20. Gfeller D, Grosdidier A, Wirth M, Daina A, Michielin O, Zoete V (2014) Swiss Target prediction: a web server for target prediction of bioactive small molecules. Nucleic Acids Res 42(1–5):W32–W38. doi:10.1093/nar/gku293

    Article  CAS  Google Scholar 

  21. Li GH, Huang JF (2012) CDRUG: a web server for predicting anticancer activity of chemical compounds. Bioinformatics 28(24):3334–3335

    Article  CAS  Google Scholar 

  22. Schrödinger Release 2015-4: Maestro, version 10.4, Schrödinger, LLC, New York, 2015

  23. http://www.acdlabs.com/products/enterprise

Download references

Authors’ contribution

AYO, SV, RSB were involved in generation of computational data on predictive properties of various Rhazya stricta’s alkaloids; NHH and AMSO participated in data acquisition. SV, JSMS and KSS were involved in overall research planning & supervision, data analysis and manuscript writing. All authors read and approved the final manuscript.

Acknowledgements

The authors gratefully acknowledge the financial support from KAU Vice President for Educational Affairs Prof. Dr. Abdulrahman O Alyoubi, The Deanship of Scientific Research (DSR) at King Abdulaziz University (KAU), Jeddah, Saudi Arabia, represented by the Unit of Strategic Technologies Research through the Project number (D 008/431) for the Project entitled: “Identification and Isolation of Salt and Heat-Tolerance Genes of R. stricta and Detection of Metabolites and their Therapeutic Effects via Cheminformatics”.

Competing interests

The authors declare that they have no competing interests.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Kulvinder Singh Saini.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Obaid, A.Y., Voleti, S., Bora, R.S. et al. Cheminformatics studies to analyze the therapeutic potential of phytochemicals from Rhazya stricta . Chemistry Central Journal 11, 11 (2017). https://doi.org/10.1186/s13065-017-0240-1

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13065-017-0240-1

Keywords