Available online on 15.04.2022 at http://jddtonline.info
Journal of Drug Delivery and Therapeutics
Open Access to Pharmaceutical and Medical Research
Copyright © 2011-2022 The Author(s): This is an open-access article distributed under the terms of the CC BY-NC 4.0 which permits unrestricted use, distribution, and reproduction in any medium for non-commercial use provided the original author and source are credited
Open Access Full Text Article Research Article
In-silico Absorption, Distribution, Metabolism, Elimination and Toxicity profile of 9,12,15-Octadecatrienoic acid (ODA) from Moringa oleifera
Parvathi K.1, Kandeepan C.2, Sabitha M.2, Senthilkumar N.3, Ramya S.4, Boopathi NM.5, Ramanathan L.6, Jayakumararaj R.6*
1 Department of Zoology, Erode Arts and Science College (Autonomous), Erode – 638009, TN, India.
2 PG and Research Department of Zoology, Arulmigu Palaniandavar College of Arts and Culture, Palani – 624601, Dindigul District, TamilNadu, India
3 Institute of Forest Genetics and Tree Breeding (IFGTB), Indian Council of Forestry Research and Education (ICFRE), Coimbatore – 641002, TamilNadu, India
4 PG Department of Zoology, Yadava College (Men), Thiruppalai - 625014, Madurai, TamilNadu, India
5 Centre for Plant Molecular Biology and Biotechnology, Department of Plant Biotechnology, Tamil Nadu Agricultural University, Coimbatore - 641 003, TN, India
6 Department of Botany, Government Arts College, Melur – 625106, Madurai District, TN, India
|
Article Info: ___________________________________________ Article History: Received 24 March 2022 Reviewed 03 April 2022 Accepted 09 April 2022 Published 15 April 2022 ___________________________________________ Cite this article as: Parvathi K., Kandeepan C, Sabitha M, Senthilkumar N, Ramya S, Boopathi NM, Ramanathan L, Jayakumararaj R, In-silico Absorption, Distribution, Metabolism, Elimination and Toxicity profile of 9,12,15-Octadecatrienoic acid (ODA) from Moringa oleifera, Journal of Drug Delivery and Therapeutics. 2022; 12(2-s):142-150 DOI: http://dx.doi.org/10.22270/jddt.v12i2-s.5289 ________________________________________ *Address for Correspondence: Jayakumararaj R., Department of Botany, Government Arts College, Melur – 625106, Madurai District, TN, India |
Abstract ____________________________________________________________________________________________________________ 9,12,15-Octadecatrienoic acid (ODA) a carboxylic acid composed of 18 carbon atoms and three cis double bonds. ODA is a plant derived essential fatty acid indispensable to the human system. ODA refers to many different structural and conformational isomers, that differ in the position of the double bonds along the backbone and on whether they are in cis ('Z') or trans ('E') conformation. It has been well established that ODA can only be outsourced from food and then converted into eicosa-pentaenoic acid (EPA) and docosa-hexaenoic acid (DHA) in the human system. However, this metabolic process is highly limited and the rate of conversion is influenced by several factors such as dose, gender, and disease. Studies suggest that ODA is associated with reduced risk of fatal ischemic heart disease. Further, higher intake may reduce the risk of sudden death among prevalent myocardial infarction in patients consistent with induced antiarrhythmic effect. ODA significantly reduces blood clots. Traditional usage of ODA is attributed to its cardiovascular-protective, anti-cancer, neuro-protective, anti-osteoporotic, anti-inflammatory, and anti-oxidative properties. Recent pharmacological indicate that ODA has anti-metabolic syndrome, anticancer, anti-inflammatory, anti-oxidant, anti-obesity, neuro-protective, and more specifically involved in the regulation of gut-micro-floral functionalities. Studies, both experimental and clinical trials indicate that ODA has anti-metabolic syndrome effects. In short, ODA is potentially used to treat many diseases, but in-depth ADMET studies are required to firmly re-establish its clinical efficacy and market potential. Keywords: ADMET; Moringa oleifera; Secondary Metabolites; Natural Products (NPs); Bioactive Substances; Octadecatrienoic acid (ODA); Eicosa-Pentaenoic Acid (EPA); Docosa-Hexaenoic Acid (DHA) |
INTRODUCTION
Octadecatrienoic acid (ODA), chemically (18:3n-3 or 3n-6) is a carboxylic acid with 18 carbons and three cis double bonds. ODA is present in either cis ('Z') or trans ('E') conformation in the plants. ODA is an essential fatty acid, meaning that human system cannot synthesize ODA and it must be supplied via the diet.1 ODA in general is in liquid state at room temperature. Table 1 shows physicochemical properties of ODA. It is commonly found in seedoil, beans, walnuts, and leafy vegetables.2,3 ODA is a precursor of two important long chain omega-3 fatty acids, EPA (eicosa pentaenoic acid, 20:3n-5) and DHA (docosa hexaenoic acid, 22:3n-6), both of them have vital role in the reduced prostate tumor growth,4 development of brain,5 cardio-vascular-health,6 metabolic-inflammatory-response,7 obesity,8 diabetics9 mental disorders10 etc. Long-chain, n-3 poly-unsaturated fatty acids (PUFAs) of EPA and DHA can be directly supplied by the intake of fish, although very few fish can provide an abundance of these PUFAs.8,9,10
Plant-derived n-3 PUFA, ODA, is an important dietary source for EPA and DHA, because of the biochemical conversion of ODA to EPA and EPA to DHA11. Although the conversion efficiency is relatively limited in humans (<8% of ODA to EPA and<4% of ODA to DHA), ODA produced by the plants is the only dietary source for the availability of n-3 PUFAs in human system. Metabolism and functional effects of plant-derived omega-3 fatty acids in humans has been studied in detail.11 Eating foods or taking supplements rich in ODA may not sometime reduce lipid levels in people with hyperlipidemia12. This boils down to the fact that there is a growing interest in using ODA for a number of purposes related to health and disease, but there isn't enough reliable information to whether it might be practically helpful.
There are 13 species in the genus Moringa, of all M. oleifera is the best known and most widely distributed species due to its multifarious uses13. M. oleifera is native to India, however, cultivated in other regions of the world14. It is a deciduous tree with brittle stem, whitish-gray corky bark with branches; leaves pale green, bipinnate/ tri-pinnate with opposite, ovate leaflets15. M. oleifera recognized as "The Miracle Tree" due to its versatile nutraceutical uses16, all parts of the plant including leaves, roots, pods, seeds, and flowers have been explored for their nutraceutical and pharmaceutical properties15. M. oleifera has been traditionally used in folk remedies across various indigenous systems of medicine. Pharmacological studies indicate that extracts obtained from the plant have antioxidants17, anti-diabetic18, anti-bacterial, anti-fungal19, and anti-carcinogenic20 properties. Interestingly, no adverse effects have been reported yet. Though, significant variation in composition of different species exists21 versatile nature of phytochemicals remains the key aspect of nutrient content. Due to overwhelming nutritive and medicinal use of the plant, it is indicated that Moringa can be widely exploited for its nutritionally important phytoconstituents in the development of functional foods, nutraceutical products and therapeutic agents as like similar medicinal plants22-44. Aim of this study is to bioprospect ODA from the leaves of MO for its molecular and biological properties to develop ODA as a novel lead in drug discovery.
MATERIALS AND METHODS
The vNN method rests on the premise that compounds with similar structures have similar activities. It is therefore reasonable to weight the contributions of neighbours so that closer neighbours contribute more to the predicted value. The vNN method calculates the similarity distance between molecules in terms of their structure, and uses a distance threshold to define a domain of applicability. This applicability domain ensures that the predictions generated are reliable. vNN models can be built within minutes and require no re-training when new assay information becomes available - an important feature when keeping quantitative structure–activity relationship (QSAR) models up-to-date to maintain their performance levels. The performance characteristics of vNN-based models are comparable, and often superior to, those of other more elaborate model constructs.45-48 One of the most widely used measures of the similarity distance between two small molecules is the Tanimoto distance, d, which is defined as:
where η(P∩Q) is the number of features common to molecules p and q, and η(P) and η(Q) are the total numbers of features for molecules p and q, respectively. The predicted biological activity y is then given by a weighted average across structurally similar neighbours:
where di denotes the Tanimoto distance between a query molecule for which a prediction is made and a molecule i of the training set; d0 is a Tanimoto-distance threshold, beyond which two molecules are no longer considered to be sufficiently similar to be included in the average; yi is the experimentally measured activity of molecule i; v denotes the total number of molecules in the training set that satisfies the condition di≤d0; and h is a smoothing factor, which dampens the distance penalty. The values of h and d0 are determined from cross-validation studies. To identify structurally similar compounds, Accelrys extended-connectivity fingerprints with a diameter of four chemical bonds (ECFP4) was used in the present study,49 as it has previously been reported to show good overall performance.48,50,51
Model Validation
A 10-fold cross-validation (CV) procedure was used to validate new models and to determine the values of the smoothing factor h and Tanimoto distance d0. In this procedure, the data was randomly divided into 10 sets, and used 9 to develop the model and the 10th to validate it, this process was repeated 10 times, leaving each set of molecules out once. In building new models, the reported averages of the 10-fold CV as the performance measures.
Performance Measures
The following metrics were used to assess model performance. (1) sensitivity measures a model’s ability to correctly detect true positives, (2) specificity measures a model’s ability to detect true negatives, (3) accuracy measures a model’s ability to make correct predictions and (4) kappa compares the probability of correct predictions to the probability of correct predictions by chance (its value ranges from +1 (perfect agreement between model prediction and experiment) to –1 (complete disagreement), with 0 indicating no agreement beyond that expected by chance).
where TP, TN, FP, and FN denote the numbers of true positives, true negatives, false positives, and false negatives, respectively. Kappa is a metric for assessing the quality of binary classifiers. Pr (e) is an estimate of the probability of a correct prediction by chance. It is calculated as:
The coverage is the proportion of test molecules with at least one nearest neighbour that meets the similarity criterion. The coverage is a measure of how many test compounds are within the applicability domain of a prediction model.
RESULTS AND DISCUSSION
3D structure, molecular and biological properties of 9,12,15-Octadecatrienoic acid in M. oleifera, its physicochemical properties, lipophilicity properties, water solubility properties, pharmacokinetic properties, druglikeness, medicinal chemistry properties, ADMET properties of ODA from M. oleifera is provided in Table 1-8 respectively. Performance measures of vNN models in 10-fold cross validation using a restricted/ unrestricted applicability domain for 9,12,15-Octadecatrienoic acid from M. oleifera (Table 9).
ADMET Predictions
The implemented Absorption, Distribution, Metabolism, Excretion and Toxicity (ADMET) prediction models, including their performance measures, have been reported previously.45 The 15 models cover a diverse set of ADMET endpoints. Some of the models have already been published, including those for Maximum Recommended Therapeutic Dose (MRTD),46 chemical mutagenicity,47 human liver microsomal (HLM),48 Pgp inhibitor/substrates.49
Liver Toxicity - DILI: Drug-induced liver injury (DILI) has been one of the most commonly cited reasons for drug withdrawals from the market. This application predicts whether a compound could cause DILI. The dataset of 1,431 compounds was obtained from four sources used by Xu et al.52 This dataset contains both pharmaceuticals and non-pharmaceuticals; A compound was classified as causing DILI if it was associated with a high risk of DILI and not if there was no such risk.
Liver Toxicity - Cytotoxicity (HepG2): Cytotoxicity is the degree to which a chemical causes damage to cells. A cytotoxicity prediction model was developed, using in vitro data on toxicity against HepG2 cells for 6,000 structurally diverse compounds, collected from ChEMBL. In developing the model, compounds with an IC50 ≤ 10 μM were considered in the in vitro assay as cytotoxic.
Metabolism - HLM: The human liver microsomal (HLM) stability assay is commonly used to identify and exclude compounds that are too rapidly metabolized. For a drug to achieve effective therapeutic concentrations in the body, it cannot be metabolized too rapidly by the liver. Compounds with a half-life of 30 minutes or longer in an HLM assay are considered as stable; otherwise they are considered unstable. HLM data was retrieved from the ChEMBL database, manually curated the data, and classified compounds as stable or unstable based on the reported half-life (T1/2 > 30 min was considered stable, and T1/2 < 30 min unstable. The final dataset contained 3,654 compounds. Of these, 2,313 leads were classified as stable and 1,341 as unstable.48
Metabolism - Cytochrome P450 enzyme (CYP) inhibition: CYPs constitute a superfamily of proteins that play an important role in the metabolism and detoxification of xenobiotics. In-vitro data was derived from five main drug-metabolizing CYPs—1A2, 3A4, 2D6, 2C9, and 2C19—to develop CYP inhibition models. CYP inhibitors data was retrieved from PubChem and classified a compound with an IC50 ≤ 10 μM for an enzyme as an inhibitor of the enzyme. Predictions for the following enzymes: CYP1A2, CYP3A4, CYP2D6, CYP2C9, and CYP2C19 have been provided.
Membrane Transporters - BBB: The blood-brain barrier (BBB) is a highly selective barrier that separates the circulating blood from the central nervous system. A vNN-based BBB model was developed, using 352 compounds whose BBB permeability values (logBB) were obtained from the literature respectively.50,51 The compounds with logBB values of less than –0.3 were classified and greater than +0.3 as BBB non-permeable and permeable.
Membrane Transporters - Pgp Substrates and Inhibitors: P-glycoprotein (Pgp) is an essential cell membrane protein that extracts many foreign substances from the cell. Cancer cells often overexpress Pgp, which increases the efflux of chemotherapeutic agents from the cell and prevents treatment by reducing the effective intracellular concentrations of such agents—a phenomenon known as multidrug resistance. For this reason, identifying compounds that can either be transported out of the cell by Pgp (substrates) or impair Pgp function (inhibitors) is of great interest. A model was developed to predict both Pgp substrates and Pgp inhibitors.29 The Pgp substrate dataset was collected by Hou and co-workers.53-55 This dataset consists of measurements of 422 substrates and 400 non-substrates. To generate a large Pgp inhibitor dataset, two datasets were combined,56,57 and removed duplicates to form a combined dataset consisting of a training set of 1,319 inhibitors and 937 non-inhibitors.
hERG (Cardiotoxicity): The human ether-à-go-go-related gene (hERG) codes for a potassium ion channel involved in the normal cardiac repolarization activity of the heart. Drug-induced blockade of hERG function can cause long QT syndrome, which may result in arrhythmia and death. As much as 282 known hERG blockers were retrieved from the literature and classified compounds with an IC50 cut-off value of 10 μM or less as blockers.54 A set of 404 compounds were collected with IC50 values greater than 10 μM from ChEMBL and classified them as non-blockers.
MMP (Mitochondrial Toxicity): Given the fundamental role of mitochondria in cellular energetics and oxidative stress, mitochondrial dysfunction has been implicated in cancer, diabetes, neurodegenerative disorders, and cardiovascular diseases. Largest dataset of chemical-induced changes in mitochondrial membrane potential (MMP) was used based on the assumption that a compound that causes mitochondrial dysfunction is also likely to reduce the MMP. A vNN-based MMP prediction model was developed using 6,261 compounds collected from a previous study that screened a library of 10,000 compounds (~8,300 unique chemicals) at 15 concentrations, each in triplicate, to measure changes in the MMP in HepG2 cells.55 The study found that 913 compounds decreased the MMP, whereas 5,395 compounds had no effect.
Mutagenicity (Ames test): Mutagens are chemicals that cause abnormal genetic mutations leading to cancer. A common way to assess a chemical’s mutagenicity is the Ames test. A prediction model was developed using a literature dataset of 6,512 compounds, of which 3,503 were Ames-positive.
MRTD: The Maximum Recommended Therapeutic Dose (MRTD) is an estimated upper daily dose that is safe. A prediction model was built based on a dataset of MRTD values publically disclosed by the FDA, mostly of single-day oral doses for an average adult with a body weight of 60 kg, for 1,220 compounds (most of which are small organic drugs). Organometallics were excluded, high-molecular weight polymers (>5,000 Da), nonorganic chemicals, mixtures of chemicals, and very small molecules (<100 Da). An external test set of 160 compounds that were collected by the FDA for validation was used in the present study. The total dataset of the model contained 1185 compounds.46
CONCLUSION
In the present study ODA from M. oleifera was ADMET predicted for functional properties. It has been well established that in the human system, ODA is converted to EPA/ DHA. EPA/ DHA endowed with cardioprotective potentials lowers blood cholesterol level and reduces the risk of heart disease. With limited data, it is not obvious to conclude that ODA of MO is safe as a dietary ingredient as evidence on risks associated with ODA remains inadequate as of now. In-silico ADMET prediction data presented in the paper is hopefully expected to assist the process of drug discovery by rapid design, evaluation, and prioritization of ODA owing to its remarkable biomedical applications.
REFERENCES
1. Burdge GC. Metabolism of α-linolenic acid in humans. Prostaglandins, leukotrienes and essential fatty acids. 2006 Sep 1; 75(3):161-8. https://doi.org/10.1016/j.plefa.2006.05.013
2. Ayerza R. Oil content and fatty acid composition of chia (Salvia hispanica L.) from five northwestern locations in Argentina. Journal of the American Oil Chemists' Society. 1995 Sep; 72(9):1079-81. https://doi.org/10.1007/BF02660727
3. Raper NR, Cronin FJ, Exler J. Omega-3 fatty acid content of the US food supply. Journal of the American College of Nutrition. 1992 Jun 1; 11(3):304-8. https://doi.org/10.1080/07315724.1992.10718231
4. Bilodeau JF, Gevariya N, Larose J, Robitaille K, Roy J, Oger C, Galano JM, Bergeron A, Durand T, Fradet Y, Julien P. Long chain omega-3 fatty acids and their oxidized metabolites are associated with reduced prostate tumor growth. Prostaglandins, Leukotrienes and Essential Fatty Acids. 2021 Jan 1; 164:102215. https://doi.org/10.1016/j.plefa.2020.102215
5. Layé S, Nadjar A, Joffre C, Bazinet RP. Anti-inflammatory effects of omega-3 fatty acids in the brain: physiological mechanisms and relevance to pharmacology. Pharmacological reviews. 2018 Jan 1; 70(1):12-38. https://doi.org/10.1124/pr.117.014092
6. Kleber ME, Delgado GE, Lorkowski S, März W, von Schacky C. Omega-3 fatty acids and mortality in patients referred for coronary angiography. The Ludwigshafen Risk and Cardiovascular Health Study. Atherosclerosis. 2016 Sep 1; 252:175-81. https://doi.org/10.1016/j.atherosclerosis.2016.06.049
7. Bazan NG. Docosanoids and elovanoids from omega-3 fatty acids are pro-homeostatic modulators of inflammatory responses, cell damage and neuroprotection. Molecular aspects of medicine. 2018 Dec 1; 64:18-33. https://doi.org/10.1016/j.mam.2018.09.003
8. Buckley JD, Howe PR. Long-chain omega-3 polyunsaturated fatty acids may be beneficial for reducing obesity-a review. Nutrients. 2010 Dec; 2(12):1212-30. https://doi.org/10.3390/nu2121212
9. Kaushik M, Mozaffarian D, Spiegelman D, Manson JE, Willett WC, Hu FB. Long-chain omega-3 fatty acids, fish intake, and the risk of type 2 diabetes mellitus. The American journal of clinical nutrition. 2009 Sep 1; 90(3):613-20. https://doi.org/10.3945/ajcn.2008.27424
10. Sanchez-Villegas A, Henríquez P, Figueiras A, Ortuño F, Lahortiga F, Martínez-González MA. Long chain omega-3 fatty acids intake, fish consumption and mental disorders in the SUN cohort study. European journal of nutrition. 2007 Sep; 46(6):337-46. https://doi.org/10.1007/s00394-007-0671-x
11. Baker EJ, Miles EA, Burdge GC, Yaqoob P, Calder PC. Metabolism and functional effects of plant-derived omega-3 fatty acids in humans. Progress in lipid research. 2016 Oct 1; 64:30-56. https://doi.org/10.1016/j.plipres.2016.07.002
12. Micallef MA, Garg ML. Beyond blood lipids: phytosterols, statins and omega-3 polyunsaturated fatty acid therapy for hyperlipidemia. The Journal of nutritional biochemistry. 2009 Dec 1; 20(12):927-39. https://doi.org/10.1016/j.jnutbio.2009.06.009
13. Padayachee B, Baijnath H. An overview of the medicinal importance of Moringaceae. Journal of Medicinal Plants Research. 2012; 6(48):5831-9.
14. Kou X, Li B, Olayanju JB, Drake JM, Chen N. Nutraceutical or pharmacological potential of Moringa oleifera Lam. Nutrients. 2018; 10(3):343. https://doi.org/10.3390/nu10030343
15. Saini RK, Sivanesan I, Keum YS. Phytochemicals of Moringa oleifera: a review of their nutritional, therapeutic and industrial significance. 3 Biotech. 2016; 6(2):1-4. https://doi.org/10.1007/s13205-016-0526-3
16. Rathnayake ARMHA, Navaratne SB, Uthpala TGG Moringa oleifera plant and the nutritional and medicinal properties of Moringa oleifera leaves Trends Pros. Process. Hort. Crops (2019), pp. 251-268
17. Luqman S, Srivastava S, Kumar R, Maurya AK, Chanda D. Experimental assessment of Moringa oleifera leaf and fruit for its antistress, antioxidant, and scavenging potential using in vitro and in vivo assays. Evidence-Based Complementary and Alternative Medicine. 2012; 2012. https://doi.org/10.1155/2012/519084
18. Gopalakrishnan L, Doriya K, Kumar DS. Moringa oleifera: A review on nutritive importance and its medicinal application. Food science and human wellness. 2016; 5(2):49-56. https://doi.org/10.1016/j.fshw.2016.04.001
19. Mehta J, Shukla A, Bukhariya V, Charde R. The magic remedy of Moringa oleifera: an overview. International Journal of Biomedical and Advance Research. 2011; 2(6):215-27. https://doi.org/10.7439/ijbar.v2i6.35
20. Jung IL. Soluble extract from Moringa oleifera leaves with a new anticancer activity. PloS one. 2014; 9(4):e95492. https://doi.org/10.1371/journal.pone.0095492
21. Asiedu-Gyekye IJ, Frimpong-Manso SA, Awortwe C, Antwi DA, Nyarko AK. Micro-and macro-elemental composition and safety evaluation of the nutraceutical Moringa oleifera leaves. Journal of Toxicology. 2014; 2014. https://doi.org/10.1155/2014/786979
22. Manavalan N, Boopathi NM, Raveendran M. Medicinal and Therapeutic Properties of Moringa. InThe Moringa Genome 2021 (pp. 31-39). Springer, Cham. https://doi.org/10.1007/978-3-030-80956-0_4
23. Boopathi NM, Raveendran M. Moringa and Its Importance. In The Moringa Genome 2021 (pp. 1-9). Springer, Cham. https://doi.org/10.1007/978-3-030-80956-0_1
24. Loganathan T, Barathinivas A, Soorya C, Balamurugan S, Nagajothi TG, Ramya S, Jayakumararaj R. Physicochemical, Druggable, ADMET Pharmacoinformatics and Therapeutic Potentials of Azadirachtin-a Prenol Lipid (Triterpenoid) from Seed Oil Extracts of Azadirachta indica A. Juss. Journal of Drug Delivery and Therapeutics. 2021; 11(5):33-46. https://doi.org/10.22270/jddt.v11i5.4981
25. Kalaimathi RV, Jeevalatha A, Basha AN, Kandeepan C, Ramya S, Loganathan T, Jayakumararaj R. In-silico Absorption, Distribution, Metabolism, Elimination and Toxicity profile of Isopulegol from Rosmarinus officinalis. Journal of Drug Delivery and Therapeutics. 2022; 12(1):102-8. https://doi.org/10.22270/jddt.v12i1.5188
26. Jeevalatha A, Kalaimathi RV, Basha AN, Kandeepan C, Ramya S, Loganathan T, Jayakumararaj R. Profile of bioactive compounds in Rosmarinus officinalis. Journal of Drug Delivery and Therapeutics. 2022; 12(1):114-22. https://doi.org/10.22270/jddt.v12i1.5189
27. Ramya S, Murugan M, Krishnaveni K, Sabitha M, Kandeepan C, Jayakumararaj R. In-silico ADMET profile of Ellagic Acid from Syzygium cumini: A Natural Biaryl Polyphenol with Therapeutic Potential to Overcome Diabetic Associated Vascular Complications. Journal of Drug Delivery and Therapeutics. 2022; 12(1):91-101. https://doi.org/10.22270/jddt.v12i1.5179
28. Soorya C, Balamurugan S, Ramya S, Neethirajan K, Kandeepan C, Jayakumararaj R. Physicochemical, ADMET and Druggable properties of Myricetin: A Key Flavonoid in Syzygium cumini that regulates metabolic inflammations. Journal of Drug Delivery and Therapeutics. 2021; 11(4):66-73 https://doi.org/10.22270/jddt.v11i4.4890
29. Sabitha M, Krishnaveni K, Murugan M, Basha AN, Pallan GA, Kandeepan C, Ramya S, Jayakumararaj R. In-silico ADMET predicated Pharmacoinformatics of Quercetin-3-Galactoside, polyphenolic compound from Azadirachta indica, a sacred tree from Hill Temple in Alagarkovil Reserve Forest, Eastern Ghats, INDIA. Journal of Drug Delivery and Therapeutics. 2021; 11(5-S):77-84. https://doi.org/10.22270/jddt.v11i5-S.5026
30. Boopathi NM, Abubakar BY. Botanical Descriptions of Moringa spp. In The Moringa Genome 2021 (pp. 11-20). Springer, Cham. https://doi.org/10.1007/978-3-030-80956-0_2
31. Loganathan T, Barathinivas A, Soorya C, Balamurugan S, Nagajothi TG, Jayakumararaj R. GCMS Profile of Bioactive Secondary Metabolites with Therapeutic Potential in the Ethanolic Leaf Extracts of Azadirachta indica: A Sacred Traditional Medicinal Plant of INDIA. Journal of Drug Delivery and Therapeutics. 2021; 11(4-S):119-26. https://doi.org/10.22270/jddt.v11i4-S.4967
32. Krishnaveni K, Sabitha M, Murugan M, Kandeepan C, Ramya S, Loganathan T, Jayakumararaj R. vNN model cross validation towards Accuracy, Sensitivity, Specificity and kappa performance measures of β-caryophyllene using a restricted-unrestricted applicability domain on Artificial Intelligence & Machine Learning approach based in-silico prediction. Journal of Drug Delivery and Therapeutics. 2022; 12(1-S):123-31. https://doi.org/10.22270/jddt.v12i1-S.5222
33. Ramya S, Jepachanderamohan PJ, Kalayanasundaram M, Jayakumararaj R. In vitro antibacterial prospective of crude leaf extracts of Melia azedarach Linn. against selected bacterial strains. Ethnobotanical Leaflets. 2009; 2009(1):32.
34. Kandeepan C, Kalaimathi RV, Jeevalatha A, Basha AN, Ramya S, Jayakumararaj R. In-silico ADMET Pharmacoinformatics of Geraniol (3, 7-dimethylocta-trans-2, 6-dien-1-ol)-acyclic monoterpene alcohol drug from Leaf Essential Oil of Cymbopogon martinii from Sirumalai Hills (Eastern Ghats), INDIA. Journal of Drug Delivery and Therapeutics. 2021; 11(4-S):109-18. https://doi.org/10.22270/jddt.v11i4-S.4965
35. Ramya S, Neethirajan K, Jayakumararaj R. Profile of bioactive compounds in Syzygium cumini-a review. Journal of Pharmacy research. 2012; 5(8):4548-53.
36. Ramya S, Krishnasamy G, Jayakumararaj R, Periathambi N, Devaraj A. Bioprospecting Solanum nigrum Linn.(Solanaceae) as a potential source of Anti-Microbial agents against selected Bacterial strains. Asian Journal of Biomedical and Pharmaceutical Sciences. 2012; 2(12):65.
37. Sundari A, Jayakumararaj R. Medicinal plants used to cure cuts and wounds in Athur region of Thoothukudi district in Tamil Nadu, India. Journal of Drug Delivery and Therapeutics. 2020; 10(6-s):26-30. https://doi.org/10.22270/jddt.v10i6-s.4429
38. Sundari A, Jayakumararaj R. Herbal remedies used to treat skin disorders in Arasankulam region of Thoothukudi District in Tamil Nadu, India. Journal of Drug Delivery and Therapeutics. 2020; 10(5):33-8. https://doi.org/10.22270/jddt.v10i5.4277
39. Rajasekaran C, Meignanam E, Vijayakumar V, Kalaivani T, Ramya S, Premkumar N, Siva R, Jayakumararaj R Investigations on antibacterial activity of leaf extracts of Azadirachta indica A. Juss (Meliaceae): a traditional medicinal plant of India. Ethnobotanical leaflets. 2008; 2008(1):161.
40. Ramya S. Jayakumararaj R In Vitro Evaluation of Antibacterial Activity Using Crude Extracts of Catharanthus roseus L.(G.) Don. Ethnobotanical Leaflets. 2008; 2008(1):140.
41. Shanmugam S, Sundari A, Muneeswaran S, Vasanth C, Jayakumararaj R, Rajendran K. Ethnobotanical Indices on medicinal plants used to treat poisonous bites in Thiruppuvanam region of Sivagangai district in Tamil Nadu, India. Journal of Drug Delivery and Therapeutics. 2020; 10(6-s):31-6. https://doi.org/10.22270/jddt.v10i6-s.4432
42. Ramya S, Alaguchamy N, Maruthappan VM, Sivaperumal R, Sivalingam M, Krishnan A, Govindaraji V, Kannan K, Jayakumararaj R. Wound healing ethnomedicinal plants popular among the Malayali tribes in Vattal Hills, Dharmapuri, TN, India. Ethnobotanical Leaflets. 2009; 2009(10):6.
43. Kadhirvel K, Ramya S, Sudha TS, Ravi AV, Rajasekaran C, Selvi RV, Jayakumararaj R. Ethnomedicinal survey on plants used by tribals in Chitteri Hills. Environ We Int J Sci Tech. 2010; 5:35-46.
44. Kandeepan C, Sabitha M, Parvathi K, Senthilkumar N, Ramya S, Boopathi NM, Jayakumararaj R. Phytochemical Screening, GCMS Profile, and In-silico properties of Bioactive Compounds in Methanolic Leaf Extracts of Moringa oleifera. Journal of Drug Delivery and Therapeutics. 2022 Mar 15; 12(2):87-99. https://doi.org/10.22270/jddt.v12i2.5250
45. Schyman, P., R. Liu, V. Desai, and A. Wallqvist. vNN web server for ADMET predictions. Frontiers in Pharmacology. 2017 December 4; 8:889. https://doi.org/10.3389/fphar.2017.00889
46. Liu, R., G. Tawa, and A. Wallqvist. Locally weighted learning methods for predicting dose-dependent toxicity with application to the human maximum recommended daily dose. Chemical Research in Toxicology. 2012; 25(10):2216-2226. https://doi.org/10.1021/tx300279f
47. Liu, R., and A. Wallqvist. Merging applicability domains for in silico assessment of chemical mutagenicity. Journal of Chemical Information and Modeling. 2014; 54(3):793-800. https://doi.org/10.1021/ci500016v
48. Liu, R., P. Schyman, and A. Wallqvist. Critically assessing the predictive power of QSAR models for human liver microsomal stability. Journal of Chemical Information and Modeling. 2015; 55(8):1566-1575. https://doi.org/10.1021/acs.jcim.5b00255
49. Schyman, P., R. Liu, and A. Wallqvist. Using the variable-nearest neighbor method to identify P-glycoprotein substrates and inhibitors. ACS Omega. 2016; 1(5):923-929. https://doi.org/10.1021/acsomega.6b00247
50. Muehlbacher, M., G. Spitzer, K. Liedl, J. Kornhuber. Qualitative prediction of blood-brain barrier permeability on a large and refined dataset. Journal of Computer-Aided Molecular Design. 2011; 25:1095. https://doi.org/10.1007/s10822-011-9478-1
51. Naef R. A generally applicable computer algorithm based on the group additivity method for the calculation of seven molecular descriptors: heat of combustion, logPO/W, logS, refractivity, polarizability, toxicity and logBB of organic compounds; scope and limits of applicability. Molecules. 2015; 20(10):18279-18351 https://doi.org/10.3390/molecules201018279
52. Xu, Y., Z. Dai, F. Chen, S. Gao, J. Pei, and L. Lai. Deep learning for drug-induced liver injury. 2015, 55(10):2085-2093. https://doi.org/10.1021/acs.jcim.5b00238
53. Schyman, P., R. Liu, and A. Wallqvist. General purpose 2D and 3D similarity approach to identify hERG blockers. Journal of Chemical Information and Modeling. 2016; 56(1):213-222 https://doi.org/10.1021/acs.jcim.5b00616
54. Attene-Ramos, M., R. Huang, S. Michael, K. Witt, A. Richard, R. Tice, A. Simeonov, C. Austin, M. Xia. Profiling of the Tox21 chemical collection for mitochondrial function to identify compounds that acutely decrease mitochondrial membrane potential. 2015; 123(1):49 https://doi.org/10.1289/ehp.1408642
55. Li, D., L. Chen, Y. Li, S. Tian, H. Sun, and T. Hou. ADMET evaluation in drug discovery. Development of in Silico prediction models for P-Glycoprotein substrates. 2014; 11(3):716-726. https://doi.org/10.1021/mp400450m
56. Broccatelli, F., E. Carosati, A. Neri, M. Frosini, L. Goracci, T. Oprea, and G. Cruciani. A novel approach for predicting P-Glycoprotein (ABCB1) inhibition using molecular interaction fields. 2011; 54(6):1740-1751. https://doi.org/10.1021/jm101421d
57. Chen, L., Y. Li, Q. Zhao, H. Peng, and T. Hou. ADME evaluation in drug discovery. 10. Predictions of P-Glycoprotein inhibitors using recursive partitioning and naive bayesian classification techniques. 2011; 8(3):889-900. https://doi.org/10.1021/mp100465q
Table 1: Structure (3D), molecular and biological properties of 9,12,15-Octadecatrienoic acid in M. oleifera
|
originalSMILES CCC=CCC=CCC=CCCCCCCCC(=O)O
|
Molecular Properties |
Calculated Values |
|
5.84 |
||
|
37.30 |
||
|
Natoms |
20 |
|
|
MW |
278.44 |
|
|
nON |
2 |
|
|
nOHNH |
1 |
|
|
Nviolations |
1 |
|
|
Nrotb |
13 |
|
|
306.47 |
||
|
Biological Properties |
Bioactivity Scores |
|
|
GPCR ligand |
0.33 |
|
|
Ion channel modulator |
0.23 |
|
|
Kinase inhibitor |
0.19 |
|
|
Nuclear receptor ligand |
0.35 |
|
|
Protease inhibitor |
0.13 |
|
|
Enzyme inhibitor |
0.42 |
Table 2 Physicochemical Properties of Octadecatrienoic Acid
|
PROPERTY |
VALUE |
|
Formula |
C18H30O2 |
|
Molecular weight |
278.43 g/mol |
|
Num. heavy atoms |
20 |
|
Num. arom. heavy atoms |
0 |
|
Fraction Csp3 |
0.61 |
|
Num. rotatable bonds |
13 |
|
Num. H-bond acceptors |
2 |
|
Num. H-bond donors |
1 |
|
Molar Refractivity |
88.99 |
|
TPSA |
37.30 Ų |
Table 3 Lipophilicity Properties of Octadecatrienoic Acid
|
PROPERTY |
VALUE |
|
Log Po/w (iLOGP) |
3.36 |
|
Log Po/w (XLOGP3) |
6.46 |
|
Log Po/w (WLOGP) |
5.66 |
|
Log Po/w (MLOGP) |
4.38 |
|
Log Po/w (SILICOS-IT) |
5.59 |
|
Consensus Log Po/w |
5.09 |
Table 4 Water Solubility Properties of Octadecatrienoic Acid
|
PROPERTY |
VALUE |
|
Log S (ESOL) |
-4.78 |
|
Solubility |
4.64e-03 mg/ml ; 1.67e-05 mol/l |
|
Class |
Moderately soluble |
|
Log S (Ali) |
-7.04 |
|
Solubility |
2.55e-05 mg/ml ; 9.16e-08 mol/l |
|
Class |
Poorly soluble |
|
Log S (SILICOS-IT) |
-3.96 |
|
Solubility |
3.08e-02 mg/ml ; 1.11e-04 mol/l |
|
Class |
Soluble |
Table 5 Pharmacokinetic Properties of Octadecatrienoic Acid
|
PROPERTY |
VALUE |
|
GI absorption |
High |
|
BBB permeant |
Yes |
|
P-gp substrate |
No |
|
CYP1A2 inhibitor |
Yes |
|
CYP2C19 inhibitor |
No |
|
CYP2C9 inhibitor |
Yes |
|
CYP2D6 inhibitor |
No |
|
CYP3A4 inhibitor |
No |
|
Log Kp (skin permeation) |
-3.41 cm/s |
Table 6 Druglikeness Properties of Octadecatrienoic Acid
|
PROPERTY |
VALUE |
|
Lipinski |
Yes; 1 violation: MLOGP>4.15 |
|
Ghose |
No; 1 violation: WLOGP>5.6 |
|
Veber |
No; 1 violation: Rotors>10 |
|
Egan |
Yes |
|
Muegge |
No; 1 violation: XLOGP3>5 |
|
Bioavailability Score |
0.85 |
Table 7 Medicinal Chemistry Properties of Octadecatrienoic Acid
|
PROPERTY |
VALUE |
|
PAINS |
0 alert |
|
Brenk |
1 alert: isolated_alkene |
|
Leadlikeness |
No; 2 violations: Rotors>7, XLOGP3>3.5 |
|
Synthetic accessibility |
3.03 |
Table 8 ADMET properties of ODA from Moringa oleifera
|
Property |
Model Name |
Predicted Value |
Unit |
|
Absorption |
Water solubility |
-5.787 |
Numeric (log mol/L) |
|
Absorption |
Caco2 permeability |
1.577 |
N (log Papp in 10-6 cm/s) |
|
Absorption |
Intestinal absorption (human) |
92.836 |
Numeric (% Absorbed) |
|
Absorption |
Skin Permeability |
-2.722 |
Numeric (log Kp) |
|
Absorption |
P-glycoprotein substrate |
No |
Categorical (Yes/No) |
|
Absorption |
P-glycoprotein I inhibitor |
No |
Categorical (Yes/No) |
|
Absorption |
P-glycoprotein II inhibitor |
No |
Categorical (Yes/No) |
|
Distribution |
VDss (human) |
-0.617 |
Numeric (log L/kg) |
|
Distribution |
Fraction unbound (human) |
0.056 |
Numeric (Fu) |
|
Distribution |
BBB permeability |
-0.115 |
Numeric (log BB) |
|
Distribution |
CNS permeability |
-1.547 |
Numeric (log PS) |
|
Metabolism |
CYP2D6 substrate |
No |
Categorical (Yes/No) |
|
Metabolism |
CYP3A4 substrate |
Yes |
Categorical (Yes/No) |
|
Metabolism |
CYP1A2 inhibitior |
Yes |
Categorical (Yes/No) |
|
Metabolism |
CYP2C19 inhibitior |
No |
Categorical (Yes/No) |
|
Metabolism |
CYP2C9 inhibitior |
No |
Categorical (Yes/No) |
|
Metabolism |
CYP2D6 inhibitior |
No |
Categorical (Yes/No) |
|
Metabolism |
CYP3A4 inhibitior |
Yes |
Categorical (Yes/No) |
|
Excretion |
Total Clearance |
1.991 |
Numeric (log ml/min/kg) |
|
Excretion |
Renal OCT2 substrate |
No |
Categorical (Yes/No) |
|
Toxicity |
AMES toxicity |
No |
Categorical (Yes/No) |
|
Toxicity |
Max. tolerated dose (human) |
-0.84 |
Numeric (log mg/kg/day) |
|
Toxicity |
hERG I inhibitor |
No |
Categorical (Yes/No) |
|
Toxicity |
hERG II inhibitor |
No |
Categorical (Yes/No) |
|
Toxicity |
Oral Rat Acute Toxicity (LD50) |
1.441 |
Numeric (mol/kg) |
|
Toxicity |
Oral Rat Chronic Toxicity (LOAEL) |
3.115 |
N (log mg/kg_bw/day) |
|
Toxicity |
Hepatotoxicity |
Yes |
Categorical (Yes/No) |
|
Toxicity |
Skin Sensitisation |
Yes |
Categorical (Yes/No) |
|
Toxicity |
T.Pyriformis toxicity |
0.722 |
Numeric (log ug/L) |
|
Toxicity |
Minnow toxicity |
-1.183 |
Numeric (log mM) |
Table 9 Performance measures of vNN models in 10-fold cross validation using a restricted/ unrestricted applicability domain for 9,12,15-Octadecatrienoic acid from M. oleifera
|
Model |
Dataa |
d0b |
hc |
Accuracy |
Sensitivity |
Specificity |
kappa |
Rd |
Coverage |
|
DILI |
1427 |
0.60 |
0.50 |
0.71 |
0.70 |
0.73 |
0.42 |
0.66 |
|
|
1.00 |
0.20 |
0.67 |
0.62 |
0.72 |
0.34 |
1.00 |
|||
|
Cytotox (hep2g) |
6097 |
0.40 |
0.20 |
0.84 |
0.88 |
0.76 |
0.64 |
0.89 |
|
|
1.00 |
0.20 |
0.84 |
0.73 |
0.89 |
0.62 |
1.00 |
|||
|
HLM |
3219 |
0.40 |
0.20 |
0.81 |
0.72 |
0.87 |
0.59 |
0.91 |
|
|
1.00 |
0.20 |
0.81 |
0.70 |
0.87 |
0.57 |
1.00 |
|||
|
CYP1A2 |
7558 |
0.50 |
0.20 |
0.90 |
0.70 |
0.95 |
0.66 |
0.75 |
|
|
1.00 |
0.20 |
0.89 |
0.61 |
0.95 |
0.60 |
1.00 |
|||
|
CYP2C9 |
8072 |
0.50 |
0.20 |
0.91 |
0.55 |
0.96 |
0.54 |
0.76 |
|
|
1.00 |
0.20 |
0.90 |
0.44 |
0.96 |
0.46 |
1.00 |
|||
|
CYP2C19 |
8155 |
0.55 |
0.20 |
0.87 |
0.64 |
0.93 |
0.58 |
0.76 |
|
|
1.00 |
0.20 |
0.86 |
0.52 |
0.94 |
0.50 |
1.00 |
|||
|
CYP2D6 |
7805 |
0.50 |
0.20 |
0.89 |
0.61 |
0.94 |
0.57 |
0.75 |
|
|
1.00 |
0.20 |
0.88 |
0.52 |
0.95 |
0.51 |
1.00 |
|||
|
CYP3A4 |
10373 |
0.50 |
0.20 |
0.88 |
0.76 |
0.92 |
0.68 |
0.78 |
|
|
1.00 |
0.20 |
0.88 |
0.69 |
0.93 |
0.64 |
1.00 |
|||
|
BBB |
353 |
0.60 |
0.20 |
0.90 |
0.94 |
0.86 |
0.80 |
0.61 |
|
|
1.00 |
0.10 |
0.82 |
0.88 |
0.75 |
0.64 |
1.00 |
|||
|
Pgp Substrate |
822 |
0.60 |
0.20 |
0.79 |
0.80 |
0.79 |
0.58 |
0.66 |
|
|
1.00 |
0.20 |
0.73 |
0.73 |
0.74 |
0.47 |
1.00 |
|||
|
Pgp Inhibitor |
2304 |
0.50 |
0.20 |
0.85 |
0.91 |
0.73 |
0.66 |
0.76 |
|
|
1.00 |
0.10 |
0.81 |
0.86 |
0.74 |
0.61 |
1.00 |
|||
|
hERG |
685 |
0.70 |
0.70 |
0.84 |
0.84 |
0.83 |
0.68 |
0.80 |
|
|
1.00 |
0.20 |
0.82 |
0.82 |
0.83 |
0.64 |
1.00 |
|||
|
MMP |
6261 |
0.50 |
0.40 |
0.89 |
0.64 |
0.94 |
0.61 |
0.69 |
|
|
1.00 |
0.20 |
0.87 |
0.52 |
0.94 |
0.50 |
1.00 |
|||
|
AMES |
6512 |
0.50 |
0.40 |
0.82 |
0.86 |
0.75 |
0.62 |
0.79 |
|
|
1.00 |
0.20 |
0.79 |
0.82 |
0.75 |
0.57 |
1.00 |
|||
|
MRTDe |
1184 |
0.60 |
0.20 |
0.79 |
0.69 |
||||
|
1.00 |
0.20 |
0.74 |
1.00 |
aNumber of compounds in the dataset; bTanimoto-distance threshold value; cSmoothing factor; dPearson’s correlation coefficient ; eRegression model.