- Research article
- Open Access
Cheminformatics studies to analyze the therapeutic potential of phytochemicals from Rhazya stricta
Chemistry Central Journal volume 11, Article number: 11 (2017)
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.
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 . 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 . 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 . 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 . In another pharmacological study, the biochemical parameters including blood lipid profile concentrations, liver enzyme activities and kidney functions were analyzed in rats . 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 . 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 . Indole-alkaloids of R. stricta had been reported to have anticancer properties . 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) , 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]).
For prediction of various therapeutic potential of these molecules, commercially and publicly available technologies as below were utilized.
Schrodinger , a scientific software that predicts drug-like properties and liabilities (viz. HERG and CNS), and ACD/Labs  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.
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.
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.
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.
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.
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%).
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.
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
Ahmad Y, Fatima K, Le Quesne PPW, Atta-ur-Rahman (1983) Further alkaloidal constituents of the leaves of Rhazya stricta. Phytochemistry 22:1017–1019
Atta-ur-Rahman, Zaman K, Habib-ur-Rehman, Malik S (1986) Studies on alkaloids of Rhazya stricta. J Nat Prod 49:1138–1139
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
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
Ali BH, Bashir AK, Tanira MO (2000) The anti-oxidant action of Rhazya stricta in rats. Phytothera Res 14:469–471
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
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
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
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
Mukhopadhayay S, Handy GA, Funayama S, Cordell GA (1981) Anticancer indole alkaloids of Rhazya stricta. J Nat Prod 44:696–700
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
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
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
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
Purnapatre K, Khattar SK, Saini KS (2008) Cytochrome P450s in the development of target-based anticancer drugs. Cancer Lett 259:1–15
Singh D, Kashyap A, Pandey RV, Saini KS (2011) Novel advances in cytochrome P450 research. Drug Discovery Today 16:793–799
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
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
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
Li GH, Huang JF (2012) CDRUG: a web server for predicting anticancer activity of chemical compounds. Bioinformatics 28(24):3334–3335
Schrödinger Release 2015-4: Maestro, version 10.4, Schrödinger, LLC, New York, 2015
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.
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”.
The authors declare that they have no competing interests.
About this article
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) doi:10.1186/s13065-017-0240-1
- Rhazya stricta
- Physiochemical properties
- Anticancer molecules
- Anti-obesity molecules