2.1 Biological data We used 53 noscapine derivatives (TABLE 1)(commonly called as noscapinoids) that were synthesized in our laboratory10,11,12,13 for QSAR model building. Anti-cancer activities of these compounds were measured against CEM, a multi-drug resistant human T-cell lymphoblastoid cell line under identical experimental conditions (to minimize bias)16 . Cells were grown in RPMI 1640 medium (Mediatech, http:// cellgro.com) supplemented with 10% fetal calf serum, 1% antibiotics (penicillin/streptomycin), 2 mM 1-glutamine, at 37? C in a humidified atmosphere with 5% CO2. Cell proliferation assay was performed in 96- well plates as previously described12. Briefly, cells were seeded in 100 ?L growth medium at a density of 5 x 103 cells per well in 96-well plates and allowed to establish for 24 h. Serially diluted concentrations (0.01?500 ?M) of noscapinoids were then added and cells were incubated for 72 h. We measured the inhibition of cell proliferation in a colorimeter by 3-(4,5-dimethylthiazol- 2-yl)-5-(3-carboxymethoxyphenyl)- 2-(4-sulfophenyl)-2 Htetrazolium inner salt (MTS) assay using the CellTiter96 AQueous One Solution Reagent (Promega, Madison, WI). Cells were incubated with MTS for 3 h and absorbance was measured at a wavelength of 490 nm using a microplate reader. The percentage of cell survival as a function of drug concentration was then plotted to determine the IC50 values (the drug concentration needed to prevent cell proliferation by 50%). 2.2 Molecular modeling The builder feature in Maestro (Schr?dinger package, version 8.5)17 was used to build the molecular models of noscapine and its analogs. An appropriate bond was assigned to each structure order using ligprep (Schr?dinger package, version 2.4)18, and a unique low-energy ring conformation with correct chirality was generated. All structures were subjected to molecular mechanics energy minimizations using Macromodel (Schr?dinger pack-age, version 9.8)19 with default settings. OPLS 2005 force field was used to assign partial atomic charges to the molecular structures.We performed complete geometric optimiza-tion of these structures using Jaguar (Schr?dinger package, ver-sion 7.7)20 to ensure that the geometry of the structure was fairly reasonable. We have used hybrid density functional theory with Becke?s three-parameter exchange potential and the Lee-Yang-Parr correlation functional (B3LYP)21,22 and basis set 3-21G*23,24,25for geometrical optimization. 2.3 Compound alignment and 3D QSAR model building. All the optimized structures along with their respected IC50 values were imported into Forge (version 10, Cresset) 26 for building field based 3D QSAR model.Out of 53 molecules, we selected 43molecules as a training set to build the model while the remaining 10 compounds served as a test set to evaluate the model.Amino noscapine was considered as template molecule for the calculation of field points. The field points were generated using XED (eXtended Electron Distribution) force field. Four differentmolecular fields such as positive and negative electrostatic, ?shape? (van der Waals), and ?hydrophobic? (a density function correlated with steric bulk and hydrophobicity) were calculated. A maximum of 200 conformations were generated for each molecule with longdistance electrostatics and attractive vdW forces turned off, forbetter conformation populations. Duplicate conformers were filtered using RMS threshold value of 0.5.All conformers found were minimized using the XED (eXtended Electron Distribution) force field with a gradient cut-off value of 0.1. The training set molecules were aligned with the template molecule using maximum common substructure conformations generated for each molecule. A special conformation hunt was performed for all the training set molecules, where the common substructure with the reference molecule was held in the same conformation and groups that were not a part of the common substructure were again conformation hunted. The resulting conformations were then scored against the reference molecule using a scoring function which uses 50% field similarity and 50% Dice volume similarity. The top scoring conformations were presented as the top scoring alignments. All the alignments were visually inspected to ensure the best possible model. The collective field points of the training set were used to derive a gauge invariant set of sampling points, which reduced the number of descriptors that needed to be considered, while ensuring that all regions around the molecule which might contribute to activity are adequately sampled. With a distance of 1? between the sample points, sample values were calculated. A sphere exclusion algorithm was run to filter this set down, and then the field value for each molecule was calculated to generate a data matrix, on which partial least squares (PLS) regression was applied to generate a PLS model . 2.4 Model validation The model was subjected to validation using test set data. The predictive capability of the developed QSAR model was validated based on several statistical tests. Prediction error sum of squares (PRESS) is a standard index to measure the accuracy of a modeling method based on the cross validation technique. The cross validation regression coefficient (R2 LOO) was calculated based on the PRESS and SSY (Sum of squares of deviation of the experimental values from their mean) using W here, Yexp, Ypred and are the experimental, predicted and mean values of experimental activity, respectively of the training set compounds. Also, the accuracy of the prediction of Y S Santosi & PK Nair Int. J. Fundamental Applied Sci. Vol. 1, No. 4 (2012) 81-86 Table I: Chemical structures of noscapine and its congeners used in the present study, along with their observed inhibitory activity of cell proliferation (CEM cell line) Sl No. Compound Structure IC50 (M) Sl No. Compound Structure IC50 (M) 1 N O O H3CO CH3 O OCH3 OCH3 O H H 16.59 x 10-6 2 N O O H3CO CH3 O OCH3 OCH3 H H 28.3 x 10- 6 3 N O O H3CO CH3 O OCH3 OCH3 O H H Cl 1.2 x 10-6 4 N O O H3CO CH3 O OCH3 OCH3 H H Cl 45.2 x 10- 6 5 N O O H3CO CH3 O OCH3 OCH3 O H H Br 1.9 x 10-6 6 N O O H3CO CH3 O OCH3 OCH3 H H Br 2.8 x 10-6 7 N O O H3CO CH3 O OCH3 OCH3 O H H F 2.3 x 10-6 8 N O O H3CO CH3 O OCH3 OCH3 H H F 15.5 x 10- 6 9 N CH3 N3 O O OCH H 3 O O OCH3 OCH3 H 2.6 x 10-6 10 N O O H3CO CH3 O OCH3 OCH3 H H N3 2.6 x 10-6 11 N CH3 NH2 O O OCH H 3 O O OCH3 OCH3 H 3.0 x 10-6 12 N O O H3CO CH3 O OCH3 OCH3 H H NH2 3.0 x 10-6 13 N CH3 NO2 O O OCH H 3 O O OCH3 OCH3 H 10.0 x 10-6 14 N O O H3CO CH3 O OCH3 OCH3 H H NO2 10.0 x 10- 6 15 N O O H3CO CH3 O OCH3 OCH3 O H H I 38.9 x 10-6 16 N O O H3CO CH3 O OCH3 OCH3 H H I 30.5 x 10- 6 Sl No. Compound Structure IC50(M) Sl No. Compound Structure IC50 (M) 17 N O O OCH H 3 N O O Cl O OCH3 OCH3 H 48.0 x 10-6 18 N O O H OCH3 N O Cl O O OCH3 OCH3 H 45.2 x 10-6 19 N O O OCH H 3 N O Cl Cl O O OCH3 OCH3 H 44.2 x 10-6 20 N O O H OCH3 N O F F O O OCH3 OCH3 H 40.6 x 10-6 21 N O O OCH H 3 N S O O OCH3 OCH3 H 41.9 x 10-6 22 N O O H OCH3 N S Cl O O OCH3 OCH3 H 45.1 x 10-6 23 N O O OCH H 3 N S OCH3 O O OCH3 OCH3 H 46.5 x 10-6 24 N O O H OCH3 N S CH3 O O OCH3 OCH3 H 44.4 x 10-6 25 N O O OCH H 3 N S O O OCH3 OCH3 H 44.1 x 10-6 26 N CH3 O O OCH H 3 O O OH OCH3 H 51.8 x 10-6 27 N O O OCH H 3 CH3 O O OH OH H 39.5 x 10-6 28 N O O H HO CH3 O O OCH3 OCH3 H 42.6 x 10-6 29 N O O H HO O O OH OH H CH3 41.0 x 10-6 30 N HO HO OCH H 3 CH3 O O OCH3 OCH3 H 42.3 x 10-6 31 N CH3 HO HO H HO O O OH OH H 35.5 x 10-6 32 N O O H OCH3 CH3 OCH3 OCH3 H HO 53.6 x 10-6 33 N CH3 O O OCH H 3 O OCH3 OCH3 H HO HO 56.0 x 10-6 34 N CH3 O O H OCH3 OCH3 OCH3 H HO HO 52.8 x 10-6 35 N O O H OCH3 O O OCH3 OCH3 H 37.4 x 10-6 36 N O O H OCH3 O CH3 O O O OCH3 OCH3 H 43.3 x 10-6 2 e x p 2 1 2 e x p 1 ( ) 1 1 ( ) n p r e d i L O O n i Y Y P R E S S R S S Y Y Y = = - = - = - - ? ? 84 the QSAR equation was validated by R2 and R2 adj. The determination coefficient in prediction (R2 test) was calculated using the following equation27: (2) where and Ytest are the predicted values based on the QSAR equation (model response) and experimental activity of the test set compounds. is the mean activity value of the training set compounds.
This body of work represents the first report of a highly predictive Field based 3D QSAR model of noscapine and its derivatives (noscapinoids). Modification of the properties in accordance with the results of this work by incorporating some structural features into the scaffold structure of noscapine will certainly achieve the final goal of improved anti-cancer noscapinoids.
The IC50 value of noscapinoids determined against human lymphoblastoid cells (CEM) was used as anti-cancer activity in evaluating the structure?activity relationships of noscapinoids quantitatively. Given a normal distribution (bellshapedcurve) of activity values of a set of compounds, the rule of thumb in developing fairly accurate QSAR models suggests that the difference between the highest and lowest biological activity of the compounds should be three to four orders of magnitude28 . However, with the array of noscapinoids examined in this study, many of the compounds (22 out of 53) have weaker activity. Given this asymmetric distribution of activity, the actual range of raw experimental IC50data in our case is, in fact is quite impressive (ranging from 1.2 to 56 ?M, Table I) and can be reasonably used to assess the structure?activity relationship among noscapinoids. Noscapine derivatives 3, 5?7, 9?14and 37-53 have significantly better activities (IC50< 16.0 ?M) than the other compounds. On the contrary, derivatives 17 ?36 (aryl substituted N-carbamoyl/N-thiocarbamoylnoscapineanalogues) generally showed very weak or no activity. We have divided the data set into training and test set compising of randomly selected 43 and 10 molecules respectively. Amino noscapine (molecule 11) is considered as template.The positive and negative electrostatic, ?shape? (van der Waals), and ?hydrophobic? (a density function correlated with steric bulk and hydrophobicity) fields were calculated for the template molecule (Figure IIa-d). The training set molecules were aligned with the template molecule using maximum common substructure conformations generated for each molecule. After the conformation hunt the top scoring conformations are presented as the top scoring alignments.Molecular field points are calculated for the given set of molecules, the positive, negative, surface and hydrophobic molecular field points of the training set molecules were showed a reasonably good alignment with the reference field points (Figure IIIa-f). Further themolecular electrostatic Potential, hydrophobic and vander walls surface map generated for the aligned training set are in synchronous with the template. A sample of the calculated field point positions of the alignment is sampled out using a sphere exclusion algorithm. These are used to calculate the descriptors to be used for model generation using PLS. The model was validated against the test set and determination coefficient (R2 test for the test set was calculated to be 0.713 indicating high predictive ability of the model. The quality of the prediction models for the training set and test set com-pounds is shown in figure IVa-b. The alignment of the calculated molecular field points of the test set with the template molecule confirmed the reliability of the models shown in figure Va-f. The QSAR model developed in this study is statistically best fitted (Rtrain 2 = 0.884)for the prediction of antitumor activities (IC50) of training and test sets of molecules, as reported in Table I. The quality of the prediction models for the training set and test set com-pounds as shown in Figure 4a-bis appreciable. The R2 and R2 LOO values (0.884 and 0.875) of the model corroborate the criteria for a highlypredictive QSAR model29,30,31and determination coefficient (R2 test for the test set was calculated to be 0.713 indicating the reliability of the model. The standard error of estimate for the model was 0.159, which is an indicator of the robustness of the fit and suggests that the predicted IC50 based on QSAR model is reliable. Therefore, the QSAR model is highly predictive.
Phthalideisoquinoline[(S)-3-((R)-4-methoxy-6- methyl-5,6,7,8-tetrahydro [1,3]dioxolo [4,5-g]isoquinolin-5-yl) -6,7-dimethoxyiso-benzofuran-1(3H)-one], also commonly known as noscapine (Figure I) is an opium alkaloid, isolated from Papaversomniferum. It is a non-addictive derivative of opium. This agent is orally available and primarily used for its antitussive (cough-suppressing) effects1, 2, 3.Recently noscapine (Figure I) and its derivatives are shown to binds with microtubules (MT) stoichiometrically and possesses anti-cancer activity4,5.However, unlike other clinically used MT-binding drugs such as vincas (MT depolymerizing) and taxanes (MT over-polymerizing and bundling), noscapine simply suppresses the dynamics of MT assembly without changing the dimer/ polymer ratio. As a result of very slow MT dynamics, noscapine blocks cell cycle progression at mitosis (prometaphase), leading eventually to apoptotic cell death in many cancer cell types4,5,6,7,8. In addition, noscapine seems to be a poor substrate for the drug pumps because it remains effective against multidrug-resistant cancer cells9 .Although noscapine is cytotoxic in a variety of different cancer cell lines in the public library of the U.S. National Cancer Institute (60- cell screen), the IC50 values remain in the high micromolar ranges (21.1 to 100 ?M). Opportunities must now be explored to acquire better and more effective deriva-tives. Indeed, some initial efforts have been quite encouraging in which some more effective derivative of the lead compound noscapine have been designed , synthesized and tested 10,11,12,13 . With recent advances in combinatorial synthesis technology several thousand compounds can be rapidly synthesized, thus expanding the relevant chemical space. Hence, there is a need for predictive models like QSAR (Quantitative Structure Activity Relationship) that will prioritize compounds for screening, aid rational synthesis and facilitate lead identification. A QSAR model is useful in relating biological activity to physico-chemical and structural descriptors of compounds. Specifically 3D-QSAR model seems more appropriate with better predictive accuracy for the newly designed compounds. A 3D-QSAR model is built by aligning three dimensional conformers of active compounds and subsequently used to score a candidate compound on the basis to be of a fitting function that evaluates the alignment of three dimensional chemical features to the model14 . To investigate the structural determinants of noscapinoids responsible for anti-cancer activity, attempts were made to develop a 3D QSAR model utilizing ?field point? descriptors using Forge V10 software (Cresset group).ForgeV10(Cresset group) is a new integrated approach to producing 3D QSAR models, utilizing ?field point? descriptors15.It uses molecular field sampling technique which has advantages over grid-based sampling methods used by many 3D QSAR model building programs. No parameterization of grid size or spacing needed. Moreover, the method is gauge invariant and reproducible and edge effects are reduced as the sample positions are placed at the most important regions of space. We believe that our ligand -based approaches in identifying ?field point? descriptors to develop 3D QSAR model will provide important insights that will lead to the design of novel noscapinoids as anti-cancer agents.
Empey DW, Laitinen LA, Young GA, Bye CE & Hughes, DT, Comparison of the antitussive effects of codeine phosphate 20 mg, dextromethorphan 30 mg and noscapine 30 mg using citric acid-induced cough in normal subjects. Eur j clin pharmacol, 16 Chopra RN, & Mukherjee BI, Dikshit BB Narcotine: its pharmacological action and therapeutic uses. Indian J Med Res, 18 (1930) 35?49. 3. Karlsson M O, Dahlstr?m B, Eckern?s SA, Johansson M & Alm AT, Pharmacokinetics of oral noscapine. Eur j clin pharmacol 39 (1990) 275-279. 4. Ye K, Ke Y, Keshava N, Shanks J, Kapp JA, Tekmal RR, Petros J, & Joshi HC, Opium alkaloid noscapine is an antitumor agent that arrests metaphase and induce apoptosis in dividing cells. Proc Natl Acad Sci USA, 95 (1998) 1601?1606. 5. Pradeep K N, Chatterji BP, Surya N V, Aneja R, Chandra R, Kanteveri S, & Joshi HC, .Rational design, synthesis and biological evaluations of amino-noscapine: A high affinity tubulin -binding noscapinoid. J Comp Aided Drug Design, 25 (2011) 443 -454 6. Ye K, Zhou J, Landen JW, Bradbury EM, & Joshi HC Sustained activation of p34(cdc2) is required for noscapineinduced apoptosis. J BiolChem 276 (2001) 46697?46700. 7. Zhou J, Gupta K, Yao J, Ye K, Panda D, Giannakakou P, & Joshi HC, Paclitaxel resistant human ovarian cancer cells undergo cJun Nh2-terminal kinase-mediated apoptosis in response to noscapine. J BiolChem 277 (2002) 39777?39785. 8. Zhou J, Panda D, Landen JW, Wilson L, & Joshi HC Minor alteration of microtubule dynamics causes loss of tension across kinetochore pairs and activates the spindle checkpoint. J Biol Chem, 277 (2002) 17200?17208. 9. Aneja R, Liu M, Yates C, Gao J, Dong X, Zhou B, Vangapandu, S N, Zhou, J, & Joshi, H C, Multidrug Resistance-Associated Protein-Over Expressing Teniposide-Resistant Human Lymphomas undergo Apoptosis by a Tubulin-Binding Agent. Cancer Res, 68 (2008) 1495?1503. 10. Aneja R, Vangapandu SN, Lopus M, Viswesarappa V G, Dhiman N, Verma A, Chandra R, Panda D, & Joshi HC, Synthesis of Microtubule-Interfering Halogenated Noscapine Analogs Perturb Mitosis in Cancer Cells followed by Cell Death. Biochem Pharmacol, 72 (2006) 415?426. 11. Aneja R, Vangapandu SN, & Joshi HC, Synthesis and Biological Evaluation of a Cyclic Ether Fluorinated Noscapine Analog. Bioorg Med Chem, 14(2006) 8352?8358. 12. Aneja R, Vangapandu SN, Lopus M, Chandra R, Panda D, & Joshi HC, Development of a Novel Nitro-Derivative of Noscapine for the Potential Treatment of Drug-Resistant Ovarian Cancer and T-Cell Lymphoma. Mol. Pharmacol, 69 (2006) 1801?1809. 13. Aggarwal S, Ghosh NN, Aneja R, Joshi HC, & Chandra R, A Convenient Synthesis of Aryl-Substituted N-Carbamoyl/NThiocarbamoylNarcotine and Related Compounds. Helvetica Chimica Acta, 85 (2002) 2458?2462. 14. Venkatarajan S M, Patel N, Bachmeier C, Mullan M, & Paris D, A 3D-QSAR model based screen for dihydropyridine-like compound library to identify inhibitors of amyloid beta (A?) production. Bioinformation , 5 (2010) 122-127. 15. Cheeseright T, Mackey M, Rose S, Vinter A, Molecular Field Extrema as Descriptors of Biological Activity: Definition and Validation. J Chem Inf Model, 46 (2006) 665-676. 16. Beck WT, & Cirtain MC, Continued expression of vinca alkaloid resistance by CCRF-CEM cells after treatment with tunicamycin or pronase. Cancer Res 42 (1982)184?189. 17. Schr?dinger, LLC. Maestro, version 8.5, Schr?dinger, LLC: New York, NY. 18. LigPrep, version 2.4, Schr?dinger, LLC, New York, NY, 2010. 19. Macro Model, version 9.8, Schr?dinger, LLC, New York, NY, 2010. 20. Schr?dinger, LLC. Jaguar, version 7.7, Schr?dinger, LLC: New York, NY. 21. Lee, C, Yang W, & Parr RG, Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density. Phys Rev B, 37 (1988) 785-789. 22. Becke A D, A new mixing of Hartree-Fock and local densityfunctional theories. J Chem Phys, 98 (1993) 1372-1377. 23. Binkley J S, Pople JA, & Hehre W J, Self-consistent molecular orbital methods. 21. Small split-valence basis sets for first-row elements. J Am Chem Soc, 102 (1980) 939-947. 24. Gordon M S, Binkley J S, Pople J A, Pietro W J, & Hehre W J, Self-consistent molecular-orbital methods. 22. Small split-valence basis sets for second-row elements. J Am Chem Soc, 104 (1982) 2797-2803. 25. Pietro W J, Francl M M, Hehre W J, Defrees D J, Pople J A, & Binkley J S, Self-consistent molecular orbital methods. Supplemented small split-valence basis sets for second-row elements. J Am Chem Soc, 104 (1982) 5039-5048. 26. Cresset Biomolecular discovery Ltd Forge V10 27. Roy P P, & Roy K, On some aspects of variable selection for partial least squares regression models. QSAR Comb Sci, 27 (2008) 302?313. 28. Checchi PM, Nettles JH, Zhou J, Snyder JP, & Joshi HC Microtubule-interacting drugs for cancer treatment. Trends Pharmacol Sci, 24 (2003)361?365 29. Naik PK, Alam A, Malhotra A, & Rizvi O, Molecular Modeling and Structure-Activity Relationship of the Podophyllotoxin and Its Congeners. J Biomol Screen, 15 (2010) 528?540. 30. Santoshi S, Naik PK & Joshi HC, Rational design of novel antimicrotubule agent (9-azido-noscapine) from quantitative structure activity relationship (QSAR) evaluation of noscapinoids. J Biomol Screen, 9 (2011) 1047-1058. 31. Naik PK, Sindhura S T, & Singh H, Quantitative StructureActivity Relationship (QSAR) of the Insecticides: The Development of Predictive In Vivo Insecticide Activity Models. SAR QSAR Env Res, 20 (2009) 551?556