The Open Microbiology Journal




ISSN: 1874-2858 ― Volume 13, 2019

Drug-likeness of Phytic Acid and Its Analogues



Amitha Joy1, 2, S. Balaji*, 3
1 Department of Biotechnology, Sahrdaya College of Engineering and Technology, Kodakara-680684, Thrissur, India;
2 R&D Centre, Bharathiar University, Coimbatore, Tamilnadu, 641046, India;
3 Department of Biotechnology, Manipal Institute of Technology, Manipal 576104, Karnataka, India

Abstract

Inositol hexakisphosphate is known to be the phosphorous reserve in plants particularly in the seeds. Though it has been known for its antinutrient properties for many years, recent research shed light to reveal it as a novel anticancer agent. Hence the present study investigates the drug-likeness of phytic acid and its analogues through bioinformatics methods. Two potential cancer drug targets such as mitogen activated kinase and inositol 1,4,5-triphosphate receptor are included in the study. Out of 50 selected analogues of phytic acid, 42 structures interact well with the chosen drug targets. The best interacting structures are 1-diphosinositol pentakisphosphate and 2,3,4,5,6-pentaphosphonooxycyclohexyl dihydrogen phosphate. For both of these structures, the negative binding energy obtained was -49.5 KJ/mol; this affirms the stability of the complex. ADME properties are also predicted to assess the drug-like properties of the compounds. The structure activity relationship model is generated for 12 compounds with experimental IC50 values.

Keywords: Analogues, binding energy, bioactivity, docking, druglikeness, phytic acid.


Article Information


Identifiers and Pagination:

Year: 2015
Volume: 9
First Page: 141
Last Page: 149
Publisher Id: TOMICROJ-9-141
DOI: 10.2174/1874285801509010141

Article History:

Received Date: 20/11/2014
Revision Received Date: 18/4/2015
Acceptance Date: 21/4/2015
Electronic publication date: 3/11/2015
Collection year: 2015

Article Metrics:

CrossRef Citations:
0

Total Statistics:

Full-Text HTML Views: 1703
Abstract HTML Views: 986
PDF Downloads: 257
Total Views/Downloads: 2946

Unique Statistics:

Full-Text HTML Views: 762
Abstract HTML Views: 631
PDF Downloads: 169
Total Views/Downloads: 1562
Geographical View

© Joy and Balaji; Licensee Bentham Open.

open-access license: This is an open access article licensed under the terms of the (https://creativecommons.org/licenses/by/4.0/legalcode), which permits unrestricted, noncommercial use, distribution and reproduction in any medium, provided the work is properly cited.


* Address correspondence to this author at the Department of Biotechnology, Manipal Institute of Technology, Manipal 576104, Karnataka, India; E-mail: biobalagi@gmail.com, s.balaji@manipal.edu




1. INTRODUCTION

Inositol hexakisphosphate (InsP6), also known as phytic acid is the foremost inositol phosphate discovered [1Vucenik I, Shamsuddin AM. Cancer inhibition by inositol hexaphosphate (IP6) and inositol: from laboratory to clinic J Nutr 2003; 133: 3778-84.]. It is the most abundant organic phosphate in the world [2Irvine FR, Schell MJ. Back in the water: the return of the inositol phosphates Reviews 2001; 2: 327-38.]. In plants, InsP6 acts as a phosphate reserve, especially in the seeds [3Mohammad A, Nawab A, Mariya K. Quantitative detection of inositol hexakisphosphate (InsP6) in crop plants using polyacry-lamide gel electrophoresis (PAGE) American J Plant Sci 2013; 4: 1-6.
[http://dx.doi.org/10.4236/ajps.2013.412A3001]
]. In nature phytic acid (C6H18O24P6) can exist as free acid, phytate or phytin, depending on the physiological pH and metal salts [4George TT, Yizhou Z, Jean L, Cheng S. Gong. Absorption of Heavy metal Ions by Immobilized Phytic acid Appl Biochem Biotechnol 1997; •••: 63-5.]. The percentage of phytate in many of the plants has been quantified [5Andargachew G, Admasu T, Girma A, et al. Nutrient composition and effect of processing on antinutritional factors and mineral bioavailability of cultivated amochi in Ethiopia J Sci Develop 2011; 1: 21-9.] and their antinutrient properties are well studied [6Bohn L, Meyer AS, Rasmussen SK. Phytate: impact on environment and human nutrition. A challenge for molecular breeding J Zhejiang Univ Sci 2008; 9: 165-91.
[http://dx.doi.org/10.1631/jzus.B0710640]
]. Recent research investigations have demonstrated its activity against cancer cell lines such as human colon cancer cell line (HT-29), erythroleukemic cell line (K-562), human liver cancer cell line (HepG2), human prostate carcinoma cell line (DU 145) and rhabdomyosarcoma cell line (RD) [7Fox CH, Eberl M. Phytic acid (IP6), novel broad spectrum anti-neoplastic agent: a systematic review Complement Ther Med 2002; 10: 229-34.
[http://dx.doi.org/10.1016/S0965-2299(02)00092-4]
, 8Norhaizan ME, Norashareena MS, Abdah MA. Antioxidant and cytotoxicity effect of ricebranphytic acid as an anticancer agent on ovarian, breast and liver cancer cell lines Malays J Nutr 2011; 17: 367-75.]. The pessimistic view of phytic acid in human diet is optimistically assigned for human health due to its anti-cancer activities in colon, prostate, metastatic and breast cancers. It is also used to chelate and clear renal stones [9Dost K, Tokul O. Determination of phytic acid in wheat and wheat products by reverse phase high performance liquid chromato-graphy Anal Chim Acta 2006; 558: 22-7.
[http://dx.doi.org/10.1016/j.aca.2005.11.035]
]. Besides in vitro studies, the in vivo studies using animal models support anticancer effects of phytic acid and hence phytic acid is referred to as ‘natural cancer fighter’ [10Norhaizan ME, Hairuszah I, Abdah MA, et al. Rice bran phytic acid (IP6) induces growth inhibition, cell cycle arrest and apoptosis on human colorectal adenocarcinoma cells J Med Plants Res 2010; 4(21): 2283-9.]. However the drug-likeness of phytic acid is not well documented in literature, some more research investigations and clinical trials need to be done to prove its drug-likeness. Computational methods play a vital role in the screening of lead molecules, designing novel drug molecules, identifying potent protein targets and also to understand the mechanism of drug action in greater detail. In particular docking and QSAR studies are effective for understanding the structural features of lead molecules. In silico approaches can assistin this regard to study the interactions of phytic acid with various drug targets. At present there is no published research work on the structure-activity relationship (SAR) of phytic acid and its analogues. Elucidating SAR can help in screening out novel inhibitors that can bind with the druggable targets. Hence the present study makes an attempt to bring forth the drug-like properties of phytic acid and its analogues through computational approach.

2. METHODOLOGY

2.1. Ligand Preparation

The 2D structures of phytic acid and its analogues were retrieved from the Pubchem repository of NCBI [11Bolton E, Wang Y, Thiessen PA, et al. PubChem: Integrated Platform of Small Molecules and Biological Activities In: Annual Reports in Computational Chemistry. Washington: American Chemical Society 2008; p. 4.
[http://dx.doi.org/10.1016/S1574-1400(08)00012-1]
]. Ligprep module [12Release S. 2013-1: LigPrep, version 31. New York, NY: Schrödinger, LLC 2014.] was used to prepare ligands. It employed series of steps that performed 3D conversions, corrections and optimization of the following structures. The Pubchem accession numbers are given below

CID477, CID890, CID107758, CID125004, CID127297, CID178749, CID439456, CID443266, CID2524165, CID4200706, CID4487899, CID10251645, CID10747577, CID14375662, CID16752671, CID16752673, CID18365880, CID21099914, CID44274820, CID44332437, CID45479488, CID46173206, CID46173281, CID46173316, CID46173429, CID46173525, CID46905360, CID52949527, CID53379838, CID53380009, CID53380097, CID53380098, CID53380099, CID53380100, CID53380197, CID53380199, CID53380200, CID53380300, CID53380301, CID53380721, CID53380834, CID53380835, CID53380836, CID53380837, CID53462026, CID53477671, CID54002314, CID54207372, CID57503246 and CID57773931.

2.2. Identification and Preparation of Protein Targets

Identification of a protein target is a crucial step in interaction studies. In the present work one such target was selected using PASS online server [13Lagunin A, Stepanchikova A, Filimonov D, et al. Pass: prediction of activity spectra for biologically active substances Bioinformatics 2000; 16: 747-8.
[http://dx.doi.org/10.1093/bioinformatics/16.8.747]
]. The server predicts several pharmacological effects and biochemical mecha-nisms on the basis of structural formula of a substance. The prediction results are analyzed based on the activity values. The compound is very likely to reveal the activity in experimentation if the activity value > 0.7. If the compound is between 0.5 and 0.7 it is expected to exhibit the activity in experiment but the likelihood is less. The substance is improbable to display the activity in experiment if the activity value is < 0.5. Based on the PASS results Inositol 1,4,5-triphosphate receptor was selected as a protein target. Literature survey suggested one more potent drug target for colon cancer mitogen activated protein kinase [14Arslan MA, Kutuk O, Basaga H. Protein kinases as drug targets in cancer Curr Cancer Drug Targets 2006; 7: 623-34.
[http://dx.doi.org/10.2174/156800906778742479]
, 15Wawszczyk J, Orchel A, Kapral M, et al. Phytic acid down-regulates IL-8 secretion from colonic Epithelial cells by influencing mitogen-activated protein Kinase signaling pathway Acta Pol Pharm 2012; 69: 1276-82.]. Their corresponding PDB ids 1N4K and 1PMQ were chosen for the study. The crystal structures of the two receptors were prepared using Proteinprep wizard of Schrödinger Release 2013-1 [16Release S. 2013-1: Schrödinger Suite 2013 Protein Preparation Wizard; Epik version 24, Schrödinger, LLC, New York, NY, 2013; Impact version 59, Schrödinger, LLC, New York, NY, 2013; Prime version 32. New York, NY: Schrödinger, LLC 2013.]. The PDB file included a co-crystallized ligand and does not include explicit hydrogens. The structures were prepared by adding hydrogens appropriately to get refined structures.

2.3. Grid Generation and Docking

Grid generation was done using the Glide module of Schrödinger Release 2013-1 [17Suite S. 2013 QM-Polarized Ligand Docking protocol; Glide version 59, Schrödinger, LLC, New York, NY, 2013; Jaguar version 80, Schrödinger, LLC, New York, NY, 2013; QSite version 59. New York, NY: Schrödinger, LLC 2013.]. The ligand poses produced are subjected to move across a chain of ranked screens that estimate the ligand’s interaction with the receptor. In grid generation, pose refers to a whole measurement of the ligand location and orientation relative to the receptor core conformation as well as rotamer set conformations. Receptor atom van der Waals radii can be scaled and either partial charges from the force field or from the input structure can be chosen. The maximum size of the enclosing grid box was set as 50Å. The process of docking a ligand into a binding site tries to generate low energy conformations. The molecules are brought together by assigning planes that contain the sites of interaction and then moving the planes while calculating binding energies of interaction [18Cohen NC. Molecular modelling in drug design. USA: Elsevier 2006.]. Here docking was done using Glide module. Glide uses a series of hierarchical filters to search for possible locations of the ligand in the active-site region of the receptor. The shape and properties of the receptor are represented on a grid by different sets of fields that provide progressively more accurate scoring of the ligand poses [19Friesner RA, Murphyv RB, Repasky MP, et al. Glide: A new approach for rapid accurate docking and scoring. 1. Method and assessment of docking accuracy J Med Chem 2002; 47: 1739-49.
[http://dx.doi.org/10.1021/jm0306430]
].

2.4. Prime MM-GBSA

The ligand binding energies were calculated using molecular mechanics with generalized Born and surface area solvation (Prime MM-GBSA) [20Suite S. 2013-1 Induced Fit Docking protocol; Glide version 59, Schrödinger, LLC, New York, NY, 2013; Prime version 32. New York, NY: Schrödinger, LLC 2013.]. Combinations of the docking method with other methods, such as MD simulation and free energy binding calculation enables rational drug design [21Godschalk F, Genheden S, Söderhjelm P. Comparison of MM/ GBSA calculations based on explicit and implicit solvent simulations Phys Chem Chem Phys 2013; 15: 7731-9.
[http://dx.doi.org/10.1039/c3cp00116d]
]. Solvation plays an important role in molecular recognition that was considered in the MM-GBSA scoring [22Kroeme RT. Structure-based drug design: docking and scoring Curr Protein Pept Sci 2007; 8: 312-28.
[http://dx.doi.org/10.2174/138920307781369382]
]. Prime MM-GBSA predicts various energies of the complex such as coulomb energy, covalent binding energy, van der Waals energy, surface area, solvation energy, prime

energy, mmgbsa energy, mmgbsa energy of the free ligand, mmgbsa energy of the uncomplexed receptor, mmgbsa energy of the ligand in the complex, mmgbsa ligand strain energy and mmgbsa free energy of binding excluding ligand strain.

2.5. Computational Prediction of ADME Properties

Absorption, distribution, metabolism, and excretion (ADME) properties of phytic acid and its derivatives were predicted by QuikProp [23Small-Molecule DD. 2013-1: QikProp, version 36. New York, NY: Schrödinger, LLC 2013.]. QuikProp also predicts physically significant descriptors and pharmaceutically relevant properties for organic structures, either individually or in batches. The QuikProp predictions are for orally delivered drugs and for non-active transport. It also checks the Lipinski’s rule of five [24Lipinski CA. Drug-like properties and the causes of poor solubility and poor permeability J Pharmacol Toxicol Methods 2008; 44: 235-49.
[http://dx.doi.org/10.1016/S1056-8719(00)00107-6]
].

2.6. QSAR Studies

Quantitative structure activity relationship of phytic acid and its derivatives were predicted by using GUSAR software. It helps to create QSAR models on the basis of the appropriate training sets [25Karthick J, Praveen Kumar PK. In-silico analysis of targeted drug delivery to hepatic cells using lipidnano-particles to treat liver diseases Asian J Pharm Tech 2013; 3: 93-7.]. It used a self-consistent regression (a regularized least-squares method) for building QSAR models [26Zakharov AV, Lagunin AA, Filimonov DA, et al. Quantitative prediction of antitarget interaction profiles for chemical compounds Chem Res Toxicol 2012; 25: 2378-85.
[http://dx.doi.org/10.1021/tx300247r]
]. Out of the 50 structures, 12 were selected for the studies for which the experimental IC50 values were available in Pubchem database. Based on the training set, IC50 values were predicted using quantitative neighbourhood of atoms (QNA model) and multilevel neighbourhood of atoms (MNA model) [27Masand VH, Mahajana DT, Patila KN, et al. General Unrestricted Structure Activity Relationships based evaluation of quinoxaline derivatives as potential influenza NS1A protein inhibitors Der Pharma Chemica 2011; 3: 517-25.]. On the basis of QNA descriptors the involvement of every atom corresponding to its experimental activity is displayed. The predicted value correlates with the activity of the whole molecule is indicated by green, blue implies that the predicted value is below the activity of the whole molecule and red indicates that the predicted value is greater than that of the whole molecule.

3. RESULTS AND DISCUSSION

The structure of phytic acid was very likely to be active against 46 targets, Table 1. It is instinctive as well as evident that phytic acid is an inositol 1,4,5-triphosphate receptor 1 antagonist. This is in concord with Sakakura et al., they suggested inositol 1,4,5-triphosphate receptor as a molecular target for gastric cancers [28Sakakura C, Hagiwara A, Fukuda K, et al. Possible involvement of inositol 1,4,5-trisphosphate receptor type 3 (IP3R3) in the peritoneal dissemination of gastric cancers Anticancer Res 2003; 23: 3691-7.]. The crystal structure of inositol 1,4,5-triphosphate receptor in complex with IP3 is available in PDB (1N4K). Inositol 1,4,5-triphosphate receptor is involved in calcium signalling process. Calcium signaling is a vibrant signalling pathway which regulates numerous cellular processes, including fertilization, cell growth, transformation, secretion, smooth muscle contraction, sensory perception and neuronal signaling [29Berridge MJ. Inositol trisphosphate and calcium signaling Nature 1993; 361: 315-25.
[http://dx.doi.org/10.1038/361315a0]
]. Binding in the IP3 receptor therefore affects the calcium signalling pathway. MAPKs (Mitogen activated kinases) are involved in directing cellular responses to a diverse array of stimuli. They regulate cell functions including proliferation, gene expression, differentiation, mitosis, cell survival, and apoptosis [30Pearson G, Robinson F, Beers Gibson T, et al. Mitogen-activated protein (MAP) kinase pathways: regulation and physiological functions Endocr Rev 2001; 22(2): 153-83.]. There are several cancer drugs which target on MAPKs [31Hilger RA, Scheulen ME, Strumberg D. The Ras-Raf-MEK-ERK pathway in the treatment of cancer Onkologie 2002; 25(6): 511-8.
[http://dx.doi.org/10.1159/000068621]
, 32Sebolt-Leopold JS. Advances in the development of cancer therapeutics directed against the RAS-mitogen-activated protein kinase pathway Clin Cancer Res 2008; 14(12): 3651-6.
[http://dx.doi.org/10.1158/1078-0432.CCR-08-0333]
].

Table 1

Predicted biological activity spectrum of phytic acid.




Table 2

Docked ligands in the order of XP score (1N4K).




Table 3

Docked ligands in the order of XP score (1PMQ).




Table 4

MM-GBSA prediction.




Table 5

ADME properties prediction.




Table 6

Cytotoxicity prediction using QNA, MNA and Combinatorial model.




3.1. Docking Results

Docking was carried out using phytic acid and its analogues using the targets, 1N4K and 1PMQ. The results show that all the selected compounds effectively bind with the selected targets. 1-diphosphoinositol pentakisphosphate (46173525) gave the best docking score of -14.29 for 1PMQ and 2,3,4,5,6-pentaphosphonooxycyclohexyl dihydrogen phosphate (178749) gave the best docking score of -14.02 for 1N4K. Table 2 and Table 3 show the top ten docked poses. Out of fifty compounds, for 1PMQ forty two gave docking scores less than -10 and for 1N4K forty scored less than -10 which shows both the protein targets have good binding affinity towards the selected compounds. All the ligands which scored less than -10 have hydrogen bond donor atom count ranging from 8 to 14, hydrogen bond acceptor atom count ranging from 12 to 30 and heavy atom count ranging from 24 to 44 according to Pubchem repository [11Bolton E, Wang Y, Thiessen PA, et al. PubChem: Integrated Platform of Small Molecules and Biological Activities In: Annual Reports in Computational Chemistry. Washington: American Chemical Society 2008; p. 4.
[http://dx.doi.org/10.1016/S1574-1400(08)00012-1]
]. These features would have influenced the binding affinity of the ligands. All the poor binders have defined and undefined steriocentred atom count ranging from 4 to 6. This sterioisomerism also would have been responsible for their poor binding affinity. Stereochemistry in drug action is of great importance in medical practice [33McConathy J, Owens MJ. Stereochemistry in drug action Prim Care Companion J Clin Psychiatry 2003; 5: 703.]. Chiral drugs have two structurally similar forms that can act in a different way in biological systems owing to their different shapes in 3-dimensional space. In a chiral environment, one enantiomer may exhibit different chemical and pharmacological behavior than the other enantiomer and the difference in 3-dimensional structure of the inactive enantiomer prevents them from having a biological effect at this binding site [34Wallace AC, Laskowski RA, Thornton JM. LIGPLOT: a program to generate schematic diagrams of protein-ligand interactions Protein Eng 1992; 8: 127-34.
[http://dx.doi.org/10.1093/protein/8.2.127]
]. This explains the impact of stereoisomerism on binding affinity.

The amino acids involved in the binding site of 1N4K include tyrosine 567, arginine 511, arginine 265, lysine 508, arginine 269, glycine 268, lysine 569, arginine 568 and threonine 267. The amino acids involved in the binding site of 1PMQ include lysine 93, tyrosine 223, lysine 191, valine 225, glutamine 255, alanine 74, aspartate 189, valine 224, arginine 107, threonine 226, serine 193, serine 72, glycine 71, valine 78 and aspargine 194. The amino acids in the binding pocket of both the proteins have some residues in common, like tyrosine, arginine, threonine, lysine and glycine. Glide scoring function is given by the equation XP GlideScore = Ecoul + EvdW + Ebind + Epenalty [35Friesner RA, Murphyv RB, Repasky MP, et al. Extra precision glide:docking and scoring incorporating a model of hydrophobic enclosurefor protein-ligand complexes J Med Chem 2006; 49: 6177-96.
[http://dx.doi.org/10.1021/jm051256o]
, 36Yadava U, Singh M, Roychoudhury M. Pyrazolo[3,4-d]pyrimidines as inhibitor of anti-coagulation and inflammation activities of phospholipase A2: insight from molecular docking studies J Biol Phys 2013; 7: 10867-967.] where XP refers to the extra precision mode, its scoring function is given by the sum of the coulombic interaction energy, van der Waals energy, binding energy and penalty which includes factors that hinders binding. A pictorical representation of the interaction between phytic acid and 1PMQ receptor Fig. (1).

Fig. (1)

Interaction of phytic acid (sticks with a mesh surface) with the binding site of the receptor 1PMQ (solid surface), rendered using Maestro (a product of Schrödinger, Inc. [16]).



3.2. MMGBSA Prediction

The result of Prime MMGBSA distinguishes strong and weak binders (Table 4). The free energy of binding was least for (1r,2R,3S,4ss,5R,6S)-2,3,4,5,6-pentakis (phosphonooxy) cyclohexyl tetrahydrogen triphosphate (53477671). This ligand has 8 phosphate groups available for binding when compared with others. It also pocess maximum numbers of hydrogen bond donor atom count, acceptor atom count, heavy atom count and complexity values of 14, 30, 44 and 1340 respectively [11Bolton E, Wang Y, Thiessen PA, et al. PubChem: Integrated Platform of Small Molecules and Biological Activities In: Annual Reports in Computational Chemistry. Washington: American Chemical Society 2008; p. 4.
[http://dx.doi.org/10.1016/S1574-1400(08)00012-1]
] when compared with the other ligands. The coulombic energy of the compound (-38.8kcalmol-1), van-der-waal’s energy (-40.1 kcalmol-1), co-valent energy (-21.6 kcalmol-1), coulombic binding energy (-11092.4 kcalmol-1) and solvation binding energy (-2168.6 kcalmol-1) were also significantly low. The presence of aromatic group and hydroxyl group bearing tyrosine, phosphate anion binder arginine and highly reactive amino group bearing lysine in the binding pocket of both the targets also account for the binding efficacy and binding energy minima [37 The Chemistry of amino acids The Biology Project, Department of Biochemistry and Molecular Physics, The University of Arizona 2004.]. The compound 2,3,4,5,6-pentaphospho-nooxycyclohexyl dihydrogen phosphate (CID178749) ranked in the top docked poses in both the protein targets have a satisfactory binding energy minima of -12.8 kcalmol-1. This compound is otherwise known as technetium-99m. Technetium-99m is a radioisotope most widely used in medicine employed in over half of all nuclear medicine procedures [38Technetium-99m Fact Sheet 2002; 2: 320-83.].

Thirty one of the compounds gave minimum energy values. The presence of hydroxyl and phosphate functional groups in these compounds enhanced their binding affinity. The negative binding free energies showed favourable ligand protein complexes.

3.3. ADME Screening

According to the Qikprop analysis 44 physically significant descriptors and pharmaceutically relevant properties of the ligands like molecular weight, total solvent accessible surface area (SASA), hydrophobic and hydrophilic component and volume were obtained. The Table 5 shows some of the predicted properties. Natural products are frequently cited as exclusion to Lipinski's rules because environment has learned to sustain low hydrophobicity and intermolecular H-bond donating ability when it needs to make biologically active compounds with high molecular weight and large numbers of rotatable bonds [39Ganesan A. The impact of natural products upon modern drug discovery Curr Opin Chem Biol 2008; 12: 306-17.
[http://dx.doi.org/10.1016/j.cbpa.2008.03.016]
].

3.4. Quantitative Structure Activity Relationship Study Using GUSAR

QSARs are mathematical models that attempt to relate the structure-derived features of a compound to its biological or physicochemical activity. QSAR works on the assumption that structurally similar compounds have similar activities [40Dehmer M, Varuza K, Bonchev D. Statistical Modelling of Molecular Descriptors in QSAR/QSPR In: Yee LC, Wei YC, Eds. Current Modeling Methods Used in QSAR/QSPR Wiley-Black well. US 2012; Vol. 2.
[http://dx.doi.org/10.1002/9783527645121]
]. In the present work, in the consensus QSAR model given by GUSAR. In the present study majority of the selected compounds correlate with experimental and predicted activity. Three models were generated by GUSAR- QNA model, MNA and Combinatorial model as given in Table 6. For QNA model and Combinatorial model the predicted IC50 values were very less deviated from the experimental IC50. QNA model gave r2 (goodness of fit) of 0.943 and Combinatorial model gave the best result with r2 value 0.960. For MNA model the r2 value is 0.660.

The statistical characteristics of the model with high value of r2, q2, F and low value of SD and the equation indicates that the model is statistically significant, could be used for analogue screening. The regression equations for the MNA model, QNA model and combinatorial model are 1.33x-7.298, 0.899x + 0.419 and 0.836x – 0.303 respectively indicates a linear relationship between experimental IC50 value and predicted IC50 value. Therefore phytic acid and the selected analogues for GUSAR studies have similar activities.

CONCLUSION

The present work identified mitogen activated kinase and inositol 1,4,5-triphosphate receptor as the protein targets for phytic acid. The interaction (docking scores) are -14.29 and -14.02, respectively. The protein ligand complexes had a minimum binding energy of -49.5kJ/mol. The parent compound, phytic acid and its similar compounds resemble similar mode of binding due to chemical similarity. The amino acid pattern involved in the binding interaction is almost same for all the compounds. The analyzed compounds are predicted to be drug-like. The structure activity relationship model correlated the atom level contributions to the activity. This would be made beneficial for designing novel drug compounds. Furthermore animal trials and clinical studies and can reveal more insights into the drug-likeness of phytic acid.

CONFLICT OF INTEREST

The authors declare that there are no conflicts of interest.

ACKNOWLEDGEMENTS

The authors show their sincere gratitude to the reviewers of this manuscript. They also extend their gratitude to the software provider Mr. Raghu Rangaswamy, Vice President, Schrödinger and Vinod Devaraji Junior Applications Scientist for the evaluation version and academic license.

REFERENCES

[1] Vucenik I, Shamsuddin AM. Cancer inhibition by inositol hexaphosphate (IP6) and inositol: from laboratory to clinic J Nutr 2003; 133: 3778-84.
[2] Irvine FR, Schell MJ. Back in the water: the return of the inositol phosphates Reviews 2001; 2: 327-38.
[3] Mohammad A, Nawab A, Mariya K. Quantitative detection of inositol hexakisphosphate (InsP6) in crop plants using polyacry-lamide gel electrophoresis (PAGE) American J Plant Sci 2013; 4: 1-6.
[http://dx.doi.org/10.4236/ajps.2013.412A3001]
[4] George TT, Yizhou Z, Jean L, Cheng S. Gong. Absorption of Heavy metal Ions by Immobilized Phytic acid Appl Biochem Biotechnol 1997; •••: 63-5.
[5] Andargachew G, Admasu T, Girma A, et al. Nutrient composition and effect of processing on antinutritional factors and mineral bioavailability of cultivated amochi in Ethiopia J Sci Develop 2011; 1: 21-9.
[6] Bohn L, Meyer AS, Rasmussen SK. Phytate: impact on environment and human nutrition. A challenge for molecular breeding J Zhejiang Univ Sci 2008; 9: 165-91.
[http://dx.doi.org/10.1631/jzus.B0710640]
[7] Fox CH, Eberl M. Phytic acid (IP6), novel broad spectrum anti-neoplastic agent: a systematic review Complement Ther Med 2002; 10: 229-34.
[http://dx.doi.org/10.1016/S0965-2299(02)00092-4]
[8] Norhaizan ME, Norashareena MS, Abdah MA. Antioxidant and cytotoxicity effect of ricebranphytic acid as an anticancer agent on ovarian, breast and liver cancer cell lines Malays J Nutr 2011; 17: 367-75.
[9] Dost K, Tokul O. Determination of phytic acid in wheat and wheat products by reverse phase high performance liquid chromato-graphy Anal Chim Acta 2006; 558: 22-7.
[http://dx.doi.org/10.1016/j.aca.2005.11.035]
[10] Norhaizan ME, Hairuszah I, Abdah MA, et al. Rice bran phytic acid (IP6) induces growth inhibition, cell cycle arrest and apoptosis on human colorectal adenocarcinoma cells J Med Plants Res 2010; 4(21): 2283-9.
[11] Bolton E, Wang Y, Thiessen PA, et al. PubChem: Integrated Platform of Small Molecules and Biological Activities In: Annual Reports in Computational Chemistry. Washington: American Chemical Society 2008; p. 4.
[http://dx.doi.org/10.1016/S1574-1400(08)00012-1]
[12] Release S. 2013-1: LigPrep, version 31. New York, NY: Schrödinger, LLC 2014.
[13] Lagunin A, Stepanchikova A, Filimonov D, et al. Pass: prediction of activity spectra for biologically active substances Bioinformatics 2000; 16: 747-8.
[http://dx.doi.org/10.1093/bioinformatics/16.8.747]
[14] Arslan MA, Kutuk O, Basaga H. Protein kinases as drug targets in cancer Curr Cancer Drug Targets 2006; 7: 623-34.
[http://dx.doi.org/10.2174/156800906778742479]
[15] Wawszczyk J, Orchel A, Kapral M, et al. Phytic acid down-regulates IL-8 secretion from colonic Epithelial cells by influencing mitogen-activated protein Kinase signaling pathway Acta Pol Pharm 2012; 69: 1276-82.
[16] Release S. 2013-1: Schrödinger Suite 2013 Protein Preparation Wizard; Epik version 24, Schrödinger, LLC, New York, NY, 2013; Impact version 59, Schrödinger, LLC, New York, NY, 2013; Prime version 32. New York, NY: Schrödinger, LLC 2013.
[17] Suite S. 2013 QM-Polarized Ligand Docking protocol; Glide version 59, Schrödinger, LLC, New York, NY, 2013; Jaguar version 80, Schrödinger, LLC, New York, NY, 2013; QSite version 59. New York, NY: Schrödinger, LLC 2013.
[18] Cohen NC. Molecular modelling in drug design. USA: Elsevier 2006.
[19] Friesner RA, Murphyv RB, Repasky MP, et al. Glide: A new approach for rapid accurate docking and scoring. 1. Method and assessment of docking accuracy J Med Chem 2002; 47: 1739-49.
[http://dx.doi.org/10.1021/jm0306430]
[20] Suite S. 2013-1 Induced Fit Docking protocol; Glide version 59, Schrödinger, LLC, New York, NY, 2013; Prime version 32. New York, NY: Schrödinger, LLC 2013.
[21] Godschalk F, Genheden S, Söderhjelm P. Comparison of MM/ GBSA calculations based on explicit and implicit solvent simulations Phys Chem Chem Phys 2013; 15: 7731-9.
[http://dx.doi.org/10.1039/c3cp00116d]
[22] Kroeme RT. Structure-based drug design: docking and scoring Curr Protein Pept Sci 2007; 8: 312-28.
[http://dx.doi.org/10.2174/138920307781369382]
[23] Small-Molecule DD. 2013-1: QikProp, version 36. New York, NY: Schrödinger, LLC 2013.
[24] Lipinski CA. Drug-like properties and the causes of poor solubility and poor permeability J Pharmacol Toxicol Methods 2008; 44: 235-49.
[http://dx.doi.org/10.1016/S1056-8719(00)00107-6]
[25] Karthick J, Praveen Kumar PK. In-silico analysis of targeted drug delivery to hepatic cells using lipidnano-particles to treat liver diseases Asian J Pharm Tech 2013; 3: 93-7.
[26] Zakharov AV, Lagunin AA, Filimonov DA, et al. Quantitative prediction of antitarget interaction profiles for chemical compounds Chem Res Toxicol 2012; 25: 2378-85.
[http://dx.doi.org/10.1021/tx300247r]
[27] Masand VH, Mahajana DT, Patila KN, et al. General Unrestricted Structure Activity Relationships based evaluation of quinoxaline derivatives as potential influenza NS1A protein inhibitors Der Pharma Chemica 2011; 3: 517-25.
[28] Sakakura C, Hagiwara A, Fukuda K, et al. Possible involvement of inositol 1,4,5-trisphosphate receptor type 3 (IP3R3) in the peritoneal dissemination of gastric cancers Anticancer Res 2003; 23: 3691-7.
[29] Berridge MJ. Inositol trisphosphate and calcium signaling Nature 1993; 361: 315-25.
[http://dx.doi.org/10.1038/361315a0]
[30] Pearson G, Robinson F, Beers Gibson T, et al. Mitogen-activated protein (MAP) kinase pathways: regulation and physiological functions Endocr Rev 2001; 22(2): 153-83.
[31] Hilger RA, Scheulen ME, Strumberg D. The Ras-Raf-MEK-ERK pathway in the treatment of cancer Onkologie 2002; 25(6): 511-8.
[http://dx.doi.org/10.1159/000068621]
[32] Sebolt-Leopold JS. Advances in the development of cancer therapeutics directed against the RAS-mitogen-activated protein kinase pathway Clin Cancer Res 2008; 14(12): 3651-6.
[http://dx.doi.org/10.1158/1078-0432.CCR-08-0333]
[33] McConathy J, Owens MJ. Stereochemistry in drug action Prim Care Companion J Clin Psychiatry 2003; 5: 703.
[34] Wallace AC, Laskowski RA, Thornton JM. LIGPLOT: a program to generate schematic diagrams of protein-ligand interactions Protein Eng 1992; 8: 127-34.
[http://dx.doi.org/10.1093/protein/8.2.127]
[35] Friesner RA, Murphyv RB, Repasky MP, et al. Extra precision glide:docking and scoring incorporating a model of hydrophobic enclosurefor protein-ligand complexes J Med Chem 2006; 49: 6177-96.
[http://dx.doi.org/10.1021/jm051256o]
[36] Yadava U, Singh M, Roychoudhury M. Pyrazolo[3,4-d]pyrimidines as inhibitor of anti-coagulation and inflammation activities of phospholipase A2: insight from molecular docking studies J Biol Phys 2013; 7: 10867-967.
[37] The Chemistry of amino acids The Biology Project, Department of Biochemistry and Molecular Physics, The University of Arizona 2004.
[38] Technetium-99m Fact Sheet 2002; 2: 320-83.
[39] Ganesan A. The impact of natural products upon modern drug discovery Curr Opin Chem Biol 2008; 12: 306-17.
[http://dx.doi.org/10.1016/j.cbpa.2008.03.016]
[40] Dehmer M, Varuza K, Bonchev D. Statistical Modelling of Molecular Descriptors in QSAR/QSPR In: Yee LC, Wei YC, Eds. Current Modeling Methods Used in QSAR/QSPR Wiley-Black well. US 2012; Vol. 2.
[http://dx.doi.org/10.1002/9783527645121]

Endorsements



"Open access will revolutionize 21st century knowledge work and accelerate the diffusion of ideas and evidence that support just in time learning and the evolution of thinking in a number of disciplines."


Daniel Pesut
(Indiana University School of Nursing, USA)

"It is important that students and researchers from all over the world can have easy access to relevant, high-standard and timely scientific information. This is exactly what Open Access Journals provide and this is the reason why I support this endeavor."


Jacques Descotes
(Centre Antipoison-Centre de Pharmacovigilance, France)

"Publishing research articles is the key for future scientific progress. Open Access publishing is therefore of utmost importance for wider dissemination of information, and will help serving the best interest of the scientific community."


Patrice Talaga
(UCB S.A., Belgium)

"Open access journals are a novel concept in the medical literature. They offer accessible information to a wide variety of individuals, including physicians, medical students, clinical investigators, and the general public. They are an outstanding source of medical and scientific information."


Jeffrey M. Weinberg
(St. Luke's-Roosevelt Hospital Center, USA)

"Open access journals are extremely useful for graduate students, investigators and all other interested persons to read important scientific articles and subscribe scientific journals. Indeed, the research articles span a wide range of area and of high quality. This is specially a must for researchers belonging to institutions with limited library facility and funding to subscribe scientific journals."


Debomoy K. Lahiri
(Indiana University School of Medicine, USA)

"Open access journals represent a major break-through in publishing. They provide easy access to the latest research on a wide variety of issues. Relevant and timely articles are made available in a fraction of the time taken by more conventional publishers. Articles are of uniformly high quality and written by the world's leading authorities."


Robert Looney
(Naval Postgraduate School, USA)

"Open access journals have transformed the way scientific data is published and disseminated: particularly, whilst ensuring a high quality standard and transparency in the editorial process, they have increased the access to the scientific literature by those researchers that have limited library support or that are working on small budgets."


Richard Reithinger
(Westat, USA)

"Not only do open access journals greatly improve the access to high quality information for scientists in the developing world, it also provides extra exposure for our papers."


J. Ferwerda
(University of Oxford, UK)

"Open Access 'Chemistry' Journals allow the dissemination of knowledge at your finger tips without paying for the scientific content."


Sean L. Kitson
(Almac Sciences, Northern Ireland)

"In principle, all scientific journals should have open access, as should be science itself. Open access journals are very helpful for students, researchers and the general public including people from institutions which do not have library or cannot afford to subscribe scientific journals. The articles are high standard and cover a wide area."


Hubert Wolterbeek
(Delft University of Technology, The Netherlands)

"The widest possible diffusion of information is critical for the advancement of science. In this perspective, open access journals are instrumental in fostering researches and achievements."


Alessandro Laviano
(Sapienza - University of Rome, Italy)

"Open access journals are very useful for all scientists as they can have quick information in the different fields of science."


Philippe Hernigou
(Paris University, France)

"There are many scientists who can not afford the rather expensive subscriptions to scientific journals. Open access journals offer a good alternative for free access to good quality scientific information."


Fidel Toldrá
(Instituto de Agroquimica y Tecnologia de Alimentos, Spain)

"Open access journals have become a fundamental tool for students, researchers, patients and the general public. Many people from institutions which do not have library or cannot afford to subscribe scientific journals benefit of them on a daily basis. The articles are among the best and cover most scientific areas."


M. Bendandi
(University Clinic of Navarre, Spain)

"These journals provide researchers with a platform for rapid, open access scientific communication. The articles are of high quality and broad scope."


Peter Chiba
(University of Vienna, Austria)

"Open access journals are probably one of the most important contributions to promote and diffuse science worldwide."


Jaime Sampaio
(University of Trás-os-Montes e Alto Douro, Portugal)

"Open access journals make up a new and rather revolutionary way to scientific publication. This option opens several quite interesting possibilities to disseminate openly and freely new knowledge and even to facilitate interpersonal communication among scientists."


Eduardo A. Castro
(INIFTA, Argentina)

"Open access journals are freely available online throughout the world, for you to read, download, copy, distribute, and use. The articles published in the open access journals are high quality and cover a wide range of fields."


Kenji Hashimoto
(Chiba University, Japan)

"Open Access journals offer an innovative and efficient way of publication for academics and professionals in a wide range of disciplines. The papers published are of high quality after rigorous peer review and they are Indexed in: major international databases. I read Open Access journals to keep abreast of the recent development in my field of study."


Daniel Shek
(Chinese University of Hong Kong, Hong Kong)

"It is a modern trend for publishers to establish open access journals. Researchers, faculty members, and students will be greatly benefited by the new journals of Bentham Science Publishers Ltd. in this category."


Jih Ru Hwu
(National Central University, Taiwan)


Browse Contents



Webmaster Contact: info@benthamopen.net
Copyright © 2019 Bentham Open