Pharmacophore modeling for identification of anti-IGF-1R drugs and in-vitro validation of fulvestrant as a potential inhibitor

PLOS ONE, May 2018

Insulin-like growth factor 1 receptor (IGF-1R) is an important therapeutic target for breast cancer treatment. The alteration in the IGF-1R associated signaling network due to various genetic and environmental factors leads the system towards metastasis. The pharmacophore modeling and logical approaches have been applied to analyze the behaviour of complex regulatory network involved in breast cancer. A total of 23 inhibitors were selected to generate ligand based pharmacophore using the tool, Molecular Operating Environment (MOE). The best model consisted of three pharmacophore features: aromatic hydrophobic (HyD/Aro), hydrophobic (HyD) and hydrogen bond acceptor (HBA). This model was validated against World drug bank (WDB) database screening to identify 189 hits with the required pharmacophore features and was further screened by using Lipinski positive compounds. Finally, the most effective drug, fulvestrant, was selected. Fulvestrant is a selective estrogen receptor down regulator (SERD). This inhibitor was further studied by using both in-silico and in-vitro approaches that showed the targeted effect of fulvestrant in ER+ MCF-7 cells. Results suggested that fulvestrant has selective cytotoxic effect and a dose dependent response on IRS-1, IGF-1R, PDZK1 and ER-α in MCF-7 cells. PDZK1 can be an important inhibitory target using fulvestrant because it directly regulates IGF-1R.

A PDF file should load here. If you do not see its contents the file may be temporarily unavailable at the journal website or you do not have a PDF plug-in installed and enabled in your browser.

Alternatively, you can download the file locally and open with any standalone PDF reader:

http://journals.plos.org/plosone/article/file?id=10.1371/journal.pone.0196312&type=printable

Pharmacophore modeling for identification of anti-IGF-1R drugs and in-vitro validation of fulvestrant as a potential inhibitor

May Pharmacophore modeling for identification of anti-IGF-1R drugs and in-vitro validation of fulvestrant as a potential inhibitor Samra Khalid 0 1 Rumeza Hanif 0 1 Ishrat Jabeen 1 Qaisar Mansoor 1 Muhammad Ismail 1 0 Atta-ur-Rahman School of Applied Biosciences, National University of Sciences and Technology , Islamabad , Pakistan , 2 Northern Institute for Cancer Research, Newcastle upon Tyne Hospitals NHS Foundation Trust, The Medical School, University of Newcastle upon Tyne , Newcastle upon Tyne , United Kingdom , 3 Research Center for Modeling & Simulation (RCMS), National University of Sciences and Technology , Islamabad , Pakistan , 4 Institute of Biomedical and Genetic Engineering (IBGE), KRL Hospital , Islamabad , Pakistan 1 Editor: Aamir Ahmad, University of South Alabama Mitchell Cancer Institute , UNITED STATES Insulin-like growth factor 1 receptor (IGF-1R) is an important therapeutic target for breast cancer treatment. The alteration in the IGF-1R associated signaling network due to various genetic and environmental factors leads the system towards metastasis. The pharmacophore modeling and logical approaches have been applied to analyze the behaviour of complex regulatory network involved in breast cancer. A total of 23 inhibitors were selected to generate ligand based pharmacophore using the tool, Molecular Operating Environment (MOE). The best model consisted of three pharmacophore features: aromatic hydrophobic (HyD/Aro), hydrophobic (HyD) and hydrogen bond acceptor (HBA). This model was validated against World drug bank (WDB) database screening to identify 189 hits with the required pharmacophore features and was further screened by using Lipinski positive compounds. Finally, the most effective drug, fulvestrant, was selected. Fulvestrant is a selective estrogen receptor down regulator (SERD). This inhibitor was further studied by using both in-silico and in-vitro approaches that showed the targeted effect of fulvestrant in ER+ MCF-7 cells. Results suggested that fulvestrant has selective cytotoxic effect and a dose dependent response on IRS-1, IGF-1R, PDZK1 and ER-α in MCF-7 cells. PDZK1 can be an important inhibitory target using fulvestrant because it directly regulates IGF-1R. - Data Availability Statement: All relevant data are within the paper and its Supporting Information files. Funding: No specific funding was secured to carry out this study. Competing interests: The authors have declared that no competing interests exist. Introduction Insulin-like growth factor type-1 receptor (IGF-1R), a trans-membrane tyrosine kinase, is involved in normal body growth and development [ 1 ]. It has two extracellular ligand binding domains, alpha (α) and beta (β) [ 2, 3 ]. IGF-1R is regulated by the binding of ligands, insulinlike growth factors such as IGF-1, to process cell proliferation and differentiation [4±6]. Previous in vivo and in vitro studies have linked higher levels of IGF-1R and its ligands with various types of cancer development and progression including breast cancer [7±10], prostate cancer [11], myeloma [ 12 ] and colon cancer [ 13, 14 ]. About 50% of the breast tumors have been reported with an over expression of IGF-1R [15]. Although several clinical trials inhibiting this receptor have been completed but unfortunately monoclonal antibodies and tyrosine kinase inhibitors targeting IGF-1R failed in phase III clinical trials for several reasons [16±18]. The activation of IGF-1R upon ligand binding induces phosphorylation of an adopter protein insulin receptor substrate-1 (IRS-1) which is also linked to various cancer subtypes [ 6, 19 ]. The signaling cascade of IGF-1R begins by the activation of several downstream mediators such as phosphoinositide3 kinase-serine/threonine protein kinases (PI3k-Akt), mitogen activated kinase-extracellular signal regulated kinase (MEK-ERK) and ataxia telangiectasia mutatedataxia telangiectasia Rad3 related (ATM-ATR) pathways [19±23]. Deregulation of these pathways induce over-expression of estrogen receptor-alpha (ER-α) which indirectly stimulates the activation of PDZ domain containing 1 (PDZK1) gene expression [24]. PDZK1 protein, also known as NHERF (Na+/H+ exchange regulatory factor), interacts with phospholipase C-β (PLC-β) and contributes to the regulation of G-protein coupled receptor (GPCR)-mediated signaling [ 25 ]. The increased expression of PDZK1 leads to the subsequent phosphorylation of ERK1/2 and calcium ions (Ca2+) signaling in response to somatostatin (SST) and IGF-1R [ 25, 26 ]. The direct molecular interaction between IGF-1R and PDZK1 enhances expression of ER-α associated with breast cancer metastasis [26]. The IGF-1R pathway facilitates loss of function mutations of multiple tumor suppressor and oncogenes including breast cancer susceptibility genes 1/2 (BRCA1/2), p53 and mouse double minute 2 homolog (Mdm2) which drastically influence resistance to apoptosis [ 20, 27 ]. This study focused on the identification of inhibitors against IGF-1R by using well-known in-silico approaches, i.e. pharmacophore modeling [28], virtual screening (VS) [ 29 ] and continuous hybrid Petri net (PN) [ 30 ]. Ligand based pharmacophore modeling is used to generate a set of chemical compounds with required pharmacophore features such as hydrophobic (HyD), aromatic (Aro), hydrogen bond acceptors (HBAs) or donors (HBDs), cations, and anions [31±33]. This ligand based modeling defines the supramolecular interactions of the above mentioned features with the desired molecular target to block its biological activity [ 32 ]. In order to identify the potential inhibitory drugs that can bind to the target, virtual screening (VS) is performed. VS is a computational drug discovery technique used to screen these chemical structures which are most likely to bind to one or more active ligands [ 33, 34 ]. This study was further enriched with continuous PN modeling [35], which allows us to analyze the delay parameters of the involved entities (proteins/genes). PN is a graph theoretical approach, which has been successfully implemented for the models and analysis of homeostatic/pathological response of IGF-1R associated network with breast cancer. The computational modeling provides a new insight to analyze the complex dynamical interactions among genes and proteins related to multifactorial diseases such as cancer. We have deployed a molecular drug screening approach which screened the drugs that bind to the active site of target molecules and inhibit their activity. The purpose of this study is to identify new IGF-1R inhibitors by using bioinformatics tools for breast cancer treatment. One of the inhibitors, fulvestrant, was validated by in-vitro experiments to understand the changes in expression levels of genes and proteins are involved in the breast cancer signaling pathway. Experimental methods Bioinformatics deals with computational approaches which have been used to identify the complex biological functions and drug designing. Such in-silico approaches are used to design the drug targets and solve various biomedical problems [ 36 ]. Furthermore, the activity of 2 / 29 potent inhibitor is further analyzed by in-vitro experiments to validate the effect of such inhibitors on genes/proteins involved in signaling pathway related to breast cancer. The methodology of both in-silico and in-vitro approaches of the current study are explained below. Dataset collection Pharmacophore modeling was first introduced by Paul Ehrlich to demonstrate the effectiveness of ligand-protein interaction by providing chemical features of the compounds [ 28 ]. The pharmacophore model generation demands for an accurate and precise input data [ 28 ]. In dataset collection, the following criteria was considered: (i) To characterize the activity of binding interaction of compounds with receptor; (ii) inhibitory potency (IC50) in a micromolar (μM) range; (iii) the most, least and inactive compounds must be included in the training set. According to the above mentioned criteria, a series of active chemical compounds were retrieved from the literature along with their inhibitory potencies (IC50) range from 0.04μM to 200μM [37±52]. The IC50 values of all the desired compounds were measured using cell viability and tyrosine kinase biological assays [37±52]. The 2-Dimensional (2D) structure of the compounds was drawn using ChemDraw program (version 8.0) [ 53 ] as shown in Fig 1. The collected dataset was classified on the basis of IC50 value, that is, most active (IC50 0.5μM), moderately active (IC50 20μM), less active (IC50 60μM) and inactive (IC50 > 60μM). Conformational dataset generation The conformations of the dataset were generated by using default settings in software MOE version 2007.09 [ 54 ]. MOE is comprehensive drug discovery software used ligand and structure-based pharmacophore modeling. However, in this study we used ligand-based pharmacophore modeling to reveal the chemical features important for the compounds activity against IGF-1R. The Scientific Vector Language (SVL) command-line is provided by MOE that integrates with database browsers which helps in the formation of packed conformation dataset [ 55 ]. The conformational search module in MOE was used to generate minimum energy conformations of each active and inactive compound. Ligand-based pharmacophore modeling A conformational training dataset of 23 chemical compounds was implemented in MOE to generate pharmacophore hypothesis. The methodology for ligand based pharmacophore modeling is explained in Fig 2. The dataset consisted of 21 active and 2 inactive compounds were employed to influence the quality of pharmacophore model depending on the two basic values are specificity and sensitivity. Active compounds have well known pharmacokinetic properties such as drug absorption, distribution, metabolism and excretion that directly bind to targeted receptors while inactive compounds had poor binding affinity [ 55 ]. In the present study, only active compounds were selected with common features such as HBDs, HBAs, Aro, HyD, number of non-polar atoms and rotatable bonds. The developed pharmacophore model was selected based on molecular sequence and low RMSD (see supporting Information S1 Table). Virtual screening (VS) The pharmacophore model developed using the most active compound picropodophyllin (0.04 μM) was further used for screening the large library of World Drug Bank (WDB) database [ 29 ]. After VS, 2534 compounds were first screened by Lipinski' positive compounds < 5 HBD groups, < 10 HBA groups and Lipinski drug-likeness [ 56 ]. Finally, drug-like compounds 3 / 29 Fig 1. Dataset: Chemical structures and activities of IGF-1R inhibitors in the training set. The dataset consists of chemical structures and activities of IGF-1R inhibitors of four most active (IC50 0.5μM), thirteen moderately active (IC50 20μM), four less active (IC50 60μM) and two were considered as inactive (IC50 > 60μM) compounds. with the features such as total number of rotatable bonds, number of rings and topological polar surface area [57±59] were further screened to examine final hit compounds. Petri net (PN) model generation The PN model was generated using the SNOOPY tool (version 2.0) [ 60 ] to model continuous dynamics based on trajectories involved in the signaling pathway related to breast cancer. In this study, we have implemented a graph theoretical hybrid PN approach which uses the kinetic logic formalism based on signaling proteins (represented as places) and interactions among proteins/genes (represented as transitions). These dynamics are specified as mass action and Michaelis-Menten equations [ 60 ]. The value of kinetic parameter is developed by the firing of transition with infinite number of tokens from source place on a target place. The firing rate is enabled to produce the expression dynamics with the same initial marking executed by multiple time simulations. The marking of places within tokens to describe the concentration of proteins are modeled to monitor the dynamics of parameters in a signaling devised by [ 61 ]. We analyzed the PN model with 10, 50 and 100 time units in simulation runs. 4 / 29 Fig 2. Workflow diagram presenting the ligand based pharmacophore modeling. (A) In the first part of pharmacophore modeling is the building of conformationaldatabase of 23 chemical compounds in MOE (i) to compute confirmations by confirmation import. (ii) Pack the conformational data-base in the MOE commandline. (iii) After packing merge the conformational data-base with original dataset of 23 compounds. (B) In the second part of building of pharmacophore modeling (i) choose an active compound such as picropodophyllin (ii) to build a pharmacophore query. (iii) Open the pharmacophore query editor to select the feature by (iv) changing the ligand annotation, number and radius of pharmacophore features under query rendering to build the pharmacophore model. (v) A click on search starts the query search to confirm the (vi) results of 21 active and 2 inactive compounds. The data was obtained by high throughput technologies (western blots, microarrays, immunohistochemistry) of several studies [6, 21, 62±66] used to validate the expression levels of individual entity involved in ER-α associated breast cancer progression. In-vitro cell culture MCF-7 cells were maintained in Dulbecco's Modified Eagle's medium (DMEM) (Invitrogen, Paisley, UK) supplemented with 10% (v/v) foetal calf serum (FCS). Cells were incubated at 37ÊC with 5% CO2 in a humidified incubator. Before treatment with fulvestrant the medium was changed and cells were cultured after 3±5 days. Cell viability assay 1×105 MCF-7 cells, obtained from the American Type Culture Collection (ATCC, Manassas, VA), were platted in a 96-well plate for 24 hours (hrs). The cells were washed with phosphate buffer saline (PBS) and treated with different concentrations (0.001±1 μM mL−1) of fulvestrant 5 / 29 Sequences (5'-3') GGGAATGGAGTGCTGTATG CACAGAAGCTTCGTTGAGAA CCCACAGTACAGCCTCACATT CACATGGTGAATGGTTTCCA CCACCAACCAGTGCACCATT GGTCTTTTCGTATCCCACCTTTC ACCTTCAACACCCCAGCCATGTACG CTGATCCACATCTGCTGGAAGGTGG (Faslodex, AstraZeneca, USA) for 48 hrs. Later, 250μg mL−1 MTT (Sigma-Aldrich, Biotechnology, USA) was added to each well and incubated for 4 hrs at 37ÊC. The cultured media was replaced with 100μL of Dimethyl sulfoxide (DMSO) (Sigma-Aldrich, Biotechnology, USA) after incubation. The absorbance was recorded at a wavelength of 570nm in a microplate reader (Techno Service, AMP Platos R496, Egypt). The %age viability of cells were calculated by using the following equation: Cell viability rate (%) = [1- (absorbance of test samples/absorbance of the control group)] x 100%. The 50% inhibitory concentration (IC50 values) was calculated from the plotted absorbance data for the dose±response curve. Reverse transcription and quantitative real-time PCR (q-RT-PCR) RNA was isolated from treated MCF-7 cells by Trizol reagent (Invitrogen, Waltham, USA). From 1 μg purified RNA, cDNA was synthesized using RevertAid™ Reverse Transcriptase (Thermo Scientific) and oligodT (Fermentas, Waltham, USA) in cDNA kit (Invitrogen, Los Angeles, USA) according to manufacturerʹs instructions. For real time PCR detection of IGF1R [ 67 ], PDZK1 [ 68 ] and ER-α [ 69 ] (primer sequences shown in Table 1), 1 μg of cDNA was amplified using a SYBR Green PCR kit (Invitrogen, Los Angeles, USA) with 0.3μM of forward and reverse primers. The PCR program was run on real time PCR (Sansure Biotech Inc, Changsha, China) and conditions were started with initial denaturation at 95ÊC for 10 min followed by 40 PCR cycles (95ÊC for 15 seconds, 60ÊC for 20 seconds and 72ÊC for 40 seconds) and final extension at 72ÊC for 5 minutes. In all q-RT-PCR experiments, β-actin [ 70 ] was amplified as a housekeeping gene to the relative mRNA of target genes IRS-1, IGF-1R, PDZK1 and ER-α to calculate ΔΔCt values. Western blotting MCF-7 cells were treated with fulvestrant (0.001±0.01 μM) for 48hrs. Cells were washed with PBS and collected to prepare cell lysate with radioimmunoprecipitate assay (RIPA) buffer: 50mM Tris-hydrochloride (HCl) pH 7.5, 150mM sodium chloride (NaCl), 1mM Ethylenediaminetetraacetic acid (EDTA), 1% Nonyl phenoxypolyethoxylethanol (NP-40) (v/v), 0.25% sodium deoxycholate (C24H40O4) (w/v), 1μgml-1 pepstatin (C34H63N5O9), 1μgml-1 aprotinin, 1μgml-1 leupeptin (C20H38N6O4), 2mM sodium orthovanadate (Na3VO4), 2mM sodium fluoride (NaF) and 2 mM PMSF (Phenylmethylsulphonyl fluoride) (Sigma-Aldrich, Dorset, United Kingdom) [ 71 ]. The protein was quantified by bicinchonic acid (BCA) assay as a standard (Thermo Scientific, Lough borough, UK) [ 72 ]. Proteins in the cell lysates and media were separated by 12% SDS-PAGE (sodium dodecyl sulfate-polyacrylamide gel electrophoresis). Proteins were transferred to 0.45μm nitrocellulose membrane (Schleicher & Schuell, Krackeler, USA) and probed with antibodies: esrogen receptor-alpha (ER-α) 1:2000 (sc-8005) (Santa 6 / 29 Cruz, Biotechnology, California); anti-PDZK1 1:3000 (#HPA006155); IRS-1 1:2000 (#3407); type I IGF receptor (#3027) (Cell Signaling Technologies, Hitchin, United Kingdom) and GAPDH (sc-25778) (Santa Cruz Biotechnology, Heidelberg, Germany). Membranes were blocked with 5% milk-TBST (25 mM Tris pH7.4, 150 mM NaCl, 0.1% Tween 20) and incubated with 1:5000 secondary antibodies (Cell Signaling Technologies, Hitchin, United Kingdom). Proteins were visualized by SuperSignal West Dura Substrate (Thermo Scientific, Rockford, USA) and exposed to X-ray film. Statistics The OD of X-ray film protein bands was measured by densitometric quantification with LabWorks 4.0 software, adjusted and normalized to GAPDH or total corresponding proteins (UVP, Inc, Cambridge, United Kingdom). Data was expressed as %age for the maximum amount of molecule detected in each experiment. Statistical differences between groups were performed by one-way analysis of variance (ANOVA) or un-paired `t' test using the GraphPad PRISM statistic software (version 7.0) (GraphPad Software, La Jolla, California USA). Statistically significant P-values (< 0.05) were considered. Results Pharmacophore model generation and evaluation Pharmacophore model was generated from the dataset of 23 compounds with inhibitory potency ranging between 0.04 to 200μM (Fig 3). These compounds have the potential to inhibit IGF-1R which interacts with various transcription factors including ER-α, ER-β, Vascular Endothelial Growth Factor Receptor-3 (VEGFR-3), human epidermal growth factor receptor 2 (HER2), IRS-1, Mitogen-activated protein kinase (MAPK) and PI3K/Akt [37±52]. Several studies have elucidated over-expression of IGF-1R in breast cancer progression [ 4, 6, 20, 21, 73, 74 ]. The statistical parameters, root mean square deviation (RMSD) (0.01±0.89) (shown in S1 Table as supporting information) of 21 active compounds were generated to build the pharmacophore model (shown in Fig 3) using the tool, molecular operating environment (MOE) version 2007.09 [54]. The parametric values included three biological features along with radius (R): F1; HyD/Aro (1.4), F2; HyD (0.8) and F3; HBA (0.9) (shown in Fig 3) represented by orange, green and purple coloured circles, respectively (shown in Fig 4). The dataset is divided into three categories: the most active (IC50 0.5μM), moderately active (IC50 20μM) and less active (IC50 60μM). The inactive compounds had poor pharmacokinetic properties such as drug absorption, distribution, metabolism, excretion and drug-drug interaction which were undetected to prove the effectiveness of our model. The selectivity of our model has been evaluated by the set of both active (truely identified as positive) and inactive (truely identified as negative) compounds. The active compounds with diverse scaffolds were matched to features of pharmacophore model while the inactive compounds were not retrieved by the model. The enrichment of pharmacophore model has been assessed by two basic values, selectivity and specificity [ 75 ]. The model with accuracy rate of 1, showed the best predictive ability with high selectivity and specificity which are defined by the retrieval of active and inactive compounds, respectively. The distances between pharmacophore features have been identified as: (F1-HyD/Aro-F2- HyD: 3.75Å), (F1-HyD/AroF3-HBA: 4.99Å) and (F2-HyD-F3-HBA: 6.79Å), shown in Table 2. The distances shown here represent how pharmacophore features are participating to interact with different chemical groups to identify potential ligands against the target IGF-1R. 7 / 29 Fig 3. Pharmacophore model and aligned ligands. The pharmacophore model of IGF-1R inhibitors are generated by MOE module with set of aligned active compounds. Pharmacophore model was generated with three chemical features such as aromatic hydrophobic (HyD/Aro); hydrophobic (HyD) and hydrogen bond acceptor (HBA) represented by orange, green and purple, respectively. Database screening Virtual screening (VS) was perform using the developed pharmacophore model against world drug bank (WDB) database [ 76 ] to identify final eight hit compounds as shown in Fig 5. Some of these compounds such as etonogestrol (DB00294), desogestrol (DB00304) and fulvestrant (DB00947) have anti-estrogenic activity while others belong to the class of organic compounds: (9aS)-4-bromo-9a-butyl-7-hydroxy-1,2,9,9a-tetrahydro-3H-fluoren-3-one (DB07757), 4-(4-hydroxyphenyl)-1-naphthaldehyde oxime (DB07150), 3-bromo-6-hydroxy2-(4-hydroxyphenyl)-1h-inden-1-one (DB07230), 4,4'-propane-2,2-diyldiphenol (DB06973) and 3-ethyl-2-(4-hydroxyphenyl)-2h-indazol-5-ol (DB07712). Etonogestrol and desogestrol are synthetically active metabolites that have high affinity for estrogen receptor (ER) and progesterone receptor (PR) [ 77, 78 ]. In the treatment of estrogen receptor positive (ER+) breast cancer, both of these drugs restrain fertility by suppressing the release of luteinizing hormone (LH), one of the contraceptive hormones imperative in ovulation [ 77, 78 ]. Fulvestrant has been used as a selective hormone receptor positive (HR+) breast cancer treatment and 8 / 29 Fig 4. Pharmacophore features. Pharmacophore features are defined as orange, green and purple colored represented by F1-aromatic hydrophobic (HyD/Aro), F2-hydrophobic (HyD) and F3-hydrogen bond acceptor (HBA). provides greater control over endocrine therapy resistance [ 79 ]. The organic compounds such as DB07757, DB07150, DB07230, DB06973 and DB07712 have transcriptional efficacy against ER-α and ER-β [80±82]. Due to unavailability of these compounds, we finally selected the most active drug fulvestrant in this study which was further evaluated by using both in-silico and in-vitro experiments. We have also found that fulvestrant is an anti-IGF-1R compound against breast cancer cells and can control tumerogenesis. This identification indicates that insilico based approaches are used to save our time and resources to identify the potential drugs against specific target to control complex diseases such as cancer. 9 / 29 Fig 5. Potential hits against IGF-1R. Out of Potential 8 hits, fulvestrant is selected for biological evaluation. Analysis of hybrid PN modeling of IGF-1R associated diseased regulatory network The hybrid PN models were constructed to observe the continuous dynamics of key proteins involved in IGF-1R associated un-treated (diseased) and treated network. Two hybrid PN models and their simulated graphs were observed to represent the mutated (given in Figs 6 and 7) and normal behaviour after treatment (given in Figs 8 and 9). Further the role of entities involved in IGF-1R signaling was evaluated. The analysis of PN model was performed to reveal the time dependent behaviour of each entity involved in diseased regulatory network (Fig 6). The signaling cascade begins by the binding of ligand IGF (token number of 5) with the receptor (IGF-1R) which leads to the phosphorylation of IRS-1. It further initiates the signaling of protein kinase PI3k which is involved in the activation of ER-α through phosphorylation of Akt. Previous studies have demonstrated that the hyper-activity of ER-α is enhanced during the pathogenesis by transcriptional activation of IGF-1 [ 20, 63, 74, 83 ]. Fig 6 illustrates that ER-α has positive feedback loop to IGF-1R which is switched on to inhibit the activity of tumor suppressor genes (TSGs) including p53 and Mdm2. Mdm2 is a key negative regulator of 10 / 29 Fig 6. Illustration of diseased IGF-1R associated hybrid PN model. In this PN model, circles (red, blue and white colored) represent continuous places which described the behaviour of entities (IGF, IGF-1R, IRS-1, Akt, PI3k, ER-α and TSGs), while squares (pink, black and white colored) represent the continuous transitions to illustrate the processes of phosphorylation, inhibition and activation. Activation (represented by directed arrows) starts signal from continuous place towards transition while inhibition (represented by inhibitory arcs) stops signal towards continuous transition. The ligand IGF-1 is given with an arbitrary token no. of five. The rate of mass action of continuous transitions of the involved entities is taken as 1. 11 / 29 Fig 7. Simulation of un-treated hybrid PN model. The simulation graph represents the relative expression level (x-axis) of IGF-1R associated entities with respect to time (y-axis). The mutated behaviour of IGF-1R associated breast cancer signaling is observed by the relative low expression levels of p53, BRCA1 and Mdm2 (dash, round dot and square dot) with the over-expression of IGF-1R (red). p53 and regulates the activity of BRCA1 which maintains the homeostatic function of system. In pathological conditions, loss of function mutations of TSGs can disrupt the function of hormonal and growth factor receptors in signaling network. The continuous parameters of the constructed IGF-1R associated signaling network were selected using tool SNOOPY (version 2.0) [ 60 ] by encoding the wet lab observed behaviours in PN modeling. The hybrid PN analysis obtained in this study resulted to identify the increased expression of IGF-1R, IRS-1 and ER-α drastically involved in increased risk of breast cancer metastasis. Our model has predicted that deregulation of IGF-1R associated signaling pathway should be controlled by the inhibition of multiple targets in order to treat breast cancer. The simulation results of IGF-1R associated diseased hybrid PN model is demonstrated in Fig 7. It shows the continous dynamics of each entity through a simulation graph to observe the relative expression levels with respect to time. It has been observed that increased expression of IGF-1R limits the homeostatic behaviour of TSGs by the inhibitory effect of ER-α [6, 20, 84±86]. The higher levels of IGF-1R and ER-α (represented by red and cyan sigmoid curves) lead to suppression of p53, BRCA1 and Mdm2 (represented by dash, round dot and 12 / 29 square dot curves) which results in the disruption of normal body functions. The IGF-1R is activated (increased level by 2 fold) through the phosphorylated activity of substrate IRS-1 (represented by pink sigmoid curve) to continue processes such as cell proliferation and DNA replication [ 6, 19 ]. The feedback regulation of ER-α has achieved by the transcriptional activation/autophosphorylation of kinases including Akt and PI3k (represented by orange and green sigmoid curves). The biological effect of E2-mediated PI3k/Akt activation is up-regulated through ER-α dependent mechanism which stimulates growth in breast cancer cells [87]. So, it is important to note that IGF-1R, IRS-1 and ER-α serve as important inhibitory targets in breast cancer treatment. Fig 8. Illustration of IGF-1R associated hybrid PN model demonstrates the inhibitory effect of inhibitor (fulvestrant). In this PN model, circles (red, blue and white colored) represent continuous places which describes the behaviour of entities (IGF, IGF-1R, IRS-1, Akt, PI3k, ER-α and TSGs) and the squares (pink, black and white colored) represent the continuous transitions to illustrate the processes of phosphorylation, inhibition and activation. Activation (represented by directed arrows) starts signal from continuous place towards transition and inhibition (represented by inhibitory arcs) stops signal towards continuous transition. The ligand IGF is given with an arbitrary token number of five. The rate of mass action of continuous transitions of the involved entities is taken as 1. 13 / 29 In-silico treatment with fulvestrant changed the relative expression levels of proteins involved in IGF-1R associated Hybrid PN modeling Previously, it has been shown that fulvestrant down-regulates the activity of ER-α and has excellent anti-proliferative efficacy in several cancers dealing with ER+ breast [ 88 ], ovarian [ 89 ], non-small cell lung [ 90 ] and gastric cancer cells [ 91 ]. The hybrid PN model was constructed to demonstrate the inhibitory effect of fulvestrant on key entities such as IGF-1R, IRS1 and ER-α involved in IGF-1R associated regulatory network (shown in Fig 8). Fulvestrant has also been to be an anti-IGF-1R drug which further controls the transcriptional activation and auto-phosphorylation of signaling kinases such as Akt (shown in Fig 8). In homeostatic conditions, ER-α and PI3k involved in direct or indirect (through E2) activation of TSGs including p53, BRCA1 and PTEN. PTEN act as a negative regulator of Akt that preferentially de-phosphorylate PI3k [ 92 ] and results in the regulation of Mdm2 which induces apoptosis in cancerous cells [ 93 ]. The regulatory effect of Mdm2 has also been achieved by the down-regulating the expression of ER-α which leads the system towards homeostasis. The treated behaviours of entities were simulated to observe the relative expression levels of IGF-1R associated hybrid PN model (Fig 8) and graph plotted with respect to time is shown in Fig 9). It shows positive network of interaction of PTEN (represented by black sigmoid curve) with p53 which is directly involved in the regulation of cell cycle [ 94 ]. The PTEN/p53 complex inhibits PI3k and Akt (represented by green and orange sigmoid curves) signaling pathway that promotes transcriptional activity of Mdm2 [ 94, 95 ]. The inhibitory actions of p53 (represented by dash sigmoid curve) and Mdm2 (represented by square dot sigmoid curve) towards ER-α (represented by cyan sigmoid curve) can also be activated by up-regulated expression of BRCA1 (represented by round dot sigmoid curve). Previous studies have confirmed that regulation of TSGs results in respective activation or deactivation of ER-α [ 63, 65, 86, 96, 97 ]. The increased expression of PTEN and TSGs has been observed in down-regulating the expression of IGF-1R (represented by red sigmoid curve), IRS-1 (represented by pink sigmoid curve) and ER-α with the help of fulvestrant. Comparison analysis of un-treated and treated hybrid PN models The comparison analysis was performed to show the simulation results of IGF-1R associated signaling genes/proteins (ER-α, IGF-1R, p53, BRCA1, Mdm2, Akt, PI3k and IRS-1) with both un-treated (represented by black sigmoid curves) and treated with fulvestrant (represented by green sigmoid curves) hybrid PN models with respect to time (Fig 10). Previously, it has been observed that IGF-1R and ER-α serve as an important inhibitory targets to control breast cancer pathogenesis [ 20, 98, 99 ]. Our results are in line with these previous experimental findings [ 6, 19, 23, 73, 85, 86, 100 ]. The differences in expression levels of proteins shown in Fig 10A± 10F based on our interpretation of the results performed by hybrid PN modeling. In the untreated model, proteins (ER-α, IGF-1R, Akt and IRS1) were shown to be over-expressed which leads towards breast cancer metastasis (Fig 10A, 10B, 10F and 10H). While the treated model, were shown the controlled expression of proteins which results in the regulation of TSGs (Fig 10C, 10D and 10E). These TSG are helped in maintaining homeostasis. Fulvestrant reduces the %age viability of breast cancer cells In the present study, cell viability was determined in MCF-7 (ATCC, Manassas, VA) cell line by MTT assay, treated with various concentrations (0.001±1μM) of fulvestrant (Faslodex, AstraZeneca, USA). Previous studies have described the mechanism of interaction of fulvestrant in MCF-7 cells to stop cell cycle proliferation [101±103] and reduced expression levels of estrogen regulated genes/proteins involved in breast cancer progression [ 79 ]. IC50 value of 14 / 29 Fig 9. Simulation of treated hybrid PN model. The simulation graph represents the relative expression level (x-axis) of IGF-1R associated entities with respect to time (y-axis). The effect of inhibitor (fulvestrant) is observed by the relative high expression levels of p53, BRCA1 and Mdm2 (dash, round dot and square dot) with the downregulation of IGF-1R (red). fulvestrant was calculated as 0.02 μM and its regression analysis showed that large R-square (R2) value of 0.88 for linear model was established under good fit. The p value was considered statistically significant (p < 0.0001) (Fig 11). It shows 50% significant reduction in the viability of MCF-7 cells in a dose dependent manner as compared to control. In contrast, fulvestrant at same concentration was only able to reduce cell viability to 110, 101, 95 and 92% respectively, in human corneal epithelial cells (HCEC) (data not shown). These results provide that fulvestrant has selective cytotoxic effect in MCF-7 cells without damaging normal cells. Treatment with fulvestrant decreased mRNA expression of target genes IGF-1R, PDZK1 and ER-α in MCF-7 cells MCF-7 cells were treated with drug fulvestrant and mRNA levels of IGF-1R, PDZK1 and ER-α gene was quantified by quantitative real-time (q-RT-PCR). A significant decrease (1.6, 1.55 and 1.4) was observed in the expressions of IGF-1R, PDZK1 and ER-α genes treated with 15 / 29 Fig 10. Comparison of simulated graphs of both un-treated and treated IGF-1R associated HPN models. The simulation graph represents the relative expression level (x-axis) of key entities with respect to time (y-axis). The black sigmoid curves represent the un-treated behaviours and green curves represent the treated behaviours of entities involved in IGF-1R associated signaling. Fig 10A±10H represents the relative change in dynamical behaviours of key proteins (ER-α, IGF-1R, p53, BRCA1, Mdm2, Akt, PI3k and IRS-1) before and after treated with anti-estrogen to be occurred. fulvestrant (shown in Fig 12). The role of ER-α in metastasis of breast cancer has been investigated and considered as important inhibitory target to control cell proliferation [ 20 ]. The metastatic response to PDZK1 and ER-α was reduced significantly in MCF-7 cells in which overexpression of IGF-1R had been decreased (Fig 12). Subsequent to the reduction in IGF-1R expression levels resulted in reduce expressions of ER-α and PDZK1, associated with decreased risk of developing breast cancer metastasis. Treatment with fulvestrant decreased protein expression of IRS-1, IGF-1R, PDZK1 and ER-α in MCF-7 cells Dose dependent treatment of ER+ MCF-7 cells with fulvestrant (0.001±0.01μM) for 48hrs shows the reduced expression levels of proteins such as IRS-1, IGF-1R, PDZK1 and ER-α (Fig 13A). Previously, fulvestrant has been reported to down-regulate the expression of PI3k in combination with the inhibitors such as BYL719, GDC-0941, GDC-0980 and BKM120 [ 104 ]. 16 / 29 Fig 11. Dose dependent effect of fulvestrant on cell viability. MTT assay of MCF-7 cells treated with various concentration of fulvestrant. Statistically significant reduced the % viability in MCF-7 cells (One-way ANOVA; p < 0.0001). Fig 13B represents the relative expression levels of un-treated proteins which were significantly higher: IRS-1 (p = 0.0002), IGF-1R (p = 0.001), PDZK1 (p = 0.0004) and ER-α (p = 0.0008) in breast cancer cells. This result is consistent with hybrid PN modeling (see section analysis of hybrid PN modeling of IGF-1R associated diseased regulatory network) showed in this study which observed the over-expression of IRS-1, IGF-1R and ER-α by down-regulating the expression of TSGs. PDZK1 is one of the important therapeutic target in breast cancer treatment that connects strong relationship with ER-α stimulated through direct interaction with IGF-1R [ 26 ]. Hence, it is important to note that the up-regulated expressions of all these proteins based on both in-silico and in-vitro analysis, serves as an important inhibitory targets for metastatic breast cancer treatment. A dose dependent response of fulvestrant was measured statistically significant to p<0.0001 on relative abundance of proteins compared with loading control (Fig 13C). It shows the up-regulated expressions of all proteins are relatively controlled by low amount of fulvestrant tend to reduce the cellular proliferation and metastasis. Our 17 / 29 Fig 12. IGF-1R, PDZK1 and ER-α mRNA levels in MCF-7 cell line treated with fulvestrant: MCF-7 cells were treated with 0.01μM fulvestrant for 48hrs prior to RNA preparation. β-actin was assessed as internal control. Data was analyzed by Graph Pad Prism and asterisks shows the statistical significance of genes IGF-1R, PDZK1 and ER-α (p<0.05) in treated cultured medium in MCF-7 cells. results suggest that fulvestrant is an effective drug that inhibits the pathogenic/carcinogenic effects of estrogen dependent IGF-1R, IRS-1, ER-α and PDZK1 signaling pathways to control breast cancer progression. Comparison of expression levels of both In-silico and In-vitro analysis of un-treated and treated behaviours of proteins in breast cancer signaling The comparison study of both in-silico (given in Figs 7, 9 and 10) and in-vitro analysis (given in Fig 13) of un-treated and treated behaviours of proteins IGF-1R, IRS-1 and ER-α showed a strong correlation in our results as presented in Table 3. The similar expression levels of protein IRS-1 was observed through in-silico simulations and in-vitro western analysis of both untreated and treated experiments. In breast cancer cells, transcriptional activity of IRS-1 can be increased by steroidal hormones such as estrogen and progesterone respectively in response to IGF-1 [ 105, 106 ]. Previously studies have showed that IRS-1 was highly expressed in MCF-7 cells and sensitized to specific chemotherapeutic agents [106±108]. We found that reducing 18 / 29 Fig 13. Western immunoblots shows dose dependent response of fulvestrant on relative abundance of proteins (IRS-1, IGF-1R, PDZK1 and ER-α) in MCF-7 cell lysates. (A) Lane 1 shows the control. Lane 2, 3 represents the protein expression profile of breast cancer cell line and following treatment with fulvestrant (0.001±0.01) for 48hrs. GAPDH was assessed as a loading control. Cells were incubated with antibodies and fixed against IRS-1, IGF-1R, PDZK1 and ER-α. Protein bands were measured and normalized to GAPDH or total individual protein expression. (B) Relative expressions of proteins compared with loading control. Asterisks shows the statistical significance of proteins cultured in serum containing medium in MCF-7 cells (unpaired t test; IRS-1 p < 0.0002, IGF-1R p < 0.0011, PDZK1 p < 0.0004 and ER-α p < 0.0008). (C) Relative abundance of dose dependent proteins was expressed as %age of maximum value measured as means ± SEM from three replicate for each blot experiment. Asterisks indicate the statistical significance of proteins cultured in serum containing treated medium in MCF-7 cells (One-way ANOVA; p < 0.01). Proteins IGF-1R IRS-1 ER-α Un-treated Treated In-silico simulations In-vitro western analysis In-silico simulations In-vitro western analysis +++ +++ ++ ++ +++ +++ ------------19 / 29 IRS-1 levels in breast cancer cells resulted in a reduced expression of ER-α. In un-treated cells, the over-expression of IGF-1R (represented by +++) considered as most active protein in insilico modeling as compared to western analysis (represented by ++) was observed. In contrast, ER-α shows the moderately active protein (represented by ++) compared to wet-lab experiments which means that there is a crosstalk between IGF-1R and ER-α, associated with an increased risk of breast cancer metastasis. These results confirm the previous studies demonstrating cross-regulation between IGF-1R and IRS-1 [ 6, 19 ], IGF-1R and ER-α [6], IRS-1 and ER-α [ 107, 109 ]. High expression of IGF-1R due to phosphorylated activity of IRS-1 can enhanced activity of ER-α which leads the tumor towards metastasis. Furthermore, the association of IGF-1R with ER-α in response to IRS-1 within the nucleus and modulates its transcriptional activity at estrogen responsive genes [ 20, 107, 109 ]. Our in-silico and in-vitro studies have shown relative decrease in expression levels of IGF-1R, IRS-1 and ER-α (represented byÐ) by the biological activity of targeted drug. Fulvestrant should be used to inhibit multiple targets (IGF-1R, IRS-1, PDZK1 and ER-α) involved in breast cancer progression. Discussion Worldwide, the most frequently diagnosed cancer among females is breast cancer [ 110 ]. It is a heterogeneous disease caused by a combination of environmental and genetic factors [ 111 ]. In ER+ breast cancer treatment, tamoxifen is used as a targeted adjuvant therapy for all the stages [ 112 ]. It is a selective ER modulator (SERM) that inhibits estrogen from binding the ligand but development of resistance towards tamoxifen in any stages of breast cancer is very common [ 112 ]. The mechanism of resistance to SERMs and other steroidal and non-steroidal classes of drugs such as aromatase inhibitors are increased by disruption in the growth factors (IGFs, EGF) signaling pathways [113±115]. Previously, in the clinical trials, various strategies (monoclonal antibodies against IGF-1R, monoclonal antibodies against IGF-1R ligands (IGF-1 and IGF-2), and IGF-1R tyrosine kinase inhibitors) have been used to control the expression levels of IGF-1R which gets over expressed in about 50% of the cases of breast cancer [ 15 ]. The aim of current study was to adopt a more comprehensive strategy by identifying new or/and known compounds (the hits mentioned in the paper) which have the potential to complete blockade IGF-1R and its signaling network (ligand and downstream signaling molecules). In the present study, we selected predictive features to build a proficient pharmacophore model (Fig 3). It was used to screen the WDB database to short list potent inhibitory drugs (DB00294, DB00304, DB00947, DB07757, DB07150, DB07230, DB06973 and DB07712) (Fig 5) against IGF-1R. Etonogestrol also known as implanon has been revealed to induce mild insulin resistance (IR) but effects are not clinically significant for healthy females [ 116 ]. One of the oral contraceptive methods is sub-dermal contraceptive implant (SCI) used in combination with etonogestrol, the synthetically active metabolite of desogestrol [ 117 ]. However, in contrast to etonogestrol and desogestrol, fulvestrant is selective estrogen receptor down regulator (SERD) that prevents dimerization of ER and proves to be effective in breast cancer treatment [ 79, 118 ]. This effect leads to abridged crosstalk signaling between ER and ER-independent growth factor signaling pathway, IGF, thus delaying resistance to endocrine treatment [ 118 ]. The results of MTT assay showed that fulvestrant significantly decreased cell viability (%) of MCF-7 cells in a dose dependent manner (Fig 11). It was observed that fulvestrant exhibits significantly selective cytotoxicity in breast cancer as compared to normal HCECs at the same concentration. Previous studies have reported that ER antagonist fulvestrant has the ability to enhance immune and chemotherapeutic-mediated cytotoxicity in lung carcinoma cells [ 119 ]. The intense collaboration of in-silico methods in drug repurposing (through ligand based pharmacophore modeling and VS) and hybrid PN modeling [ 120 ] was performed to find the 20 / 29 significance of fulvestrant on activity of proteins involved in IGF-1R signaling. The powerful synergy of PN model (Fig 6) was constructed based on various studies performed in wet lab using techniques such as: reverse phase protein analysis, DNA sequencing and copy number, immunohistochemistry, micro RNA, polymorphisms and western analysis [6, 21, 62±66, 83, 121]. Studies showed that increased expression of estrogen by ligand binding interaction (IGFs with IGF-1R) can regulate the downstream signaling mediators such as IRS-1, Akt and PI3k in breast cancer cells. Our simulation results suggested that the mutated behaviours of p53, BRCA1 and Mdm2 (represented by dash, round dot and square dot curves) were down-regulated (Fig 7). The significant up-regulated expression of IGF-1R is increased by the suppression of TSGs which lead tumor towards metastasis. It was concluded that high levels of IGF-1R, IRS-1 and ER-α can be controlled by inhibitor, fulvestrant. The therapeutic role of fulvestrant has demonstrated to be effective in the treatment of locally advanced ER or/and PR positive and HER2 negative breast cancer patients [ 79 ]. It gives a more prominent control over endocrine treatment consistence, decreasing oral ingestion and pharmacokinetic interactions with sustenance or other medications, which are vital perspectives to be considered in patients with breast cancer [ 79 ]. Due to preferable understanding of biological signaling pathways engrossed in the tumor growth and development, a few trials are assessing the role of new organic medications to blocking regulatory network. Some of the new drugs such as everolimus (mTOR inhibitor) [ 122 ], BYL719, GDC-0941, GDC-0980, and BKM120 (PI3k inhibitors) are being evaluated in combination with fulvestrant for the treatment of breast cancer metastasis [ 79 ]. Fig 8 describes a treatment algorithms of IGF-1R associated signaling network involved in breast cancer metastasis. The simulation results showed the effective behaviour of inhibitor against IGF-1R, IRS-1 and ER-α which leads the system towards homeostasis. The increased expression of PTEN controlling the autophosphorylation and over-expression of hormonal receptors (Fig 9) were observed by hybrid PN modeling. The differences between expression level of both un-treated and treated behaviours of IGF-1R associated entities were clearly seen in comparison study to analyse the inhibitory effect of fulvestrant (Fig 10). The levels of IGF-1R, IRS-1, ER-α and Akt are relatively decreased (represented by green sigmoid curves) by the up-regulated expressions of TSGs as compared to un-treated dynamical behaviours (represented by black sigmoid curves) with respect to time. It is based on our interpretation obtained in this study through in-silico modeling and simulations that fulvestrant act as an effective inhibitor against multiple breast cancer inhibitory targets to control the effect of metastasis. After ensemble PN model, which is constructed through Kyoto Encyclopedia of Genes and Genomes (KEGG) [ 123 ] or literature databases of interactions among genes and proteins [6, 20, 21, 62±66, 83, 113, 121], we selected most effective compound fulvestrant for subsequent wet experiments. The IC50 value of fulvestrant was 0.02μM as determined by MTT assay therefore we selected this compound for further investigations. Fulvestrant is currently used in long term estrogen deprivation (LTED) therapy of ER+ breast cancer cells. It was recently identified as a potential anti-cancer which inhibits PI3k-Akt, ERK/MAPK and Janus kinase (JAK)-signal transducer and activator of transcription (STAT) pathways involved in cell survival, proliferation, invasion and drug resistance [ 103 ]. Based upon our data, we constructed a resultant mechanism of IGF-1R associated breast cancer cellular pathway by which targeted growth factors to overcome endocrine therapy resistance (Fig 14). The q-RT-PCR and western blot analysis illustrated in Figs 12 and 13, acquire a more aggressive biological phenotype and how fulvestrant interplay in this signaling pathway. In summary, breast cancer cells exhibits activation of multiple growth promoting factors: IGF-1R, IRS-1, PDZK1 and ER-α. They also show an enhanced response to IGF, EGF, ER and antagonist, which in turn activates PI3k-Akt, ATM/ATR, MAPK/ERK1/2 and G-protein coupled receptor (GPCR) signaling pathways, 21 / 29 Fig 14. Pathways by which targeting growth factor and estrogen receptor to overcome endocrine therapy resistance in breast cancer. A mechanism begins by the cell signals (represented by stars) transduction through trans-membrane receptors (IGF-1R, EGFR and estrogen) and antagonist from extracellular matrix to nucleus. Estrogen up-regulates the transcription factors (TF) which regulates the activity of growth factor IGF-1. Binding of ligands IGF and EGF with targeted receptors (IGF1R and EGFR) involve in phosphorylation and trans-activation of several downstream mediator proteins lead the system towards metastasis. G-protein coupled receptor (GPCR) transmits the signal by the binding of isomers (Gα) which enhance the transcription of cyclin D1 by over-expressing PDZK1. The up-regulated expression of growth factors could lead to the activation of estrogen and progesterone receptors (ER, PR). Fulvestrant (Ful) inhibits the phospho-IRS-1, IGF-1R, ER-α and PDZK1, involved in breast cancer invasion, cell survival, metastasis, angiogenesis and cell proliferation. which cross-talk with elevated estrogens. Treatment of cells with fulvestrant significantly inhibits ER-α dependent IGF-1R, phospho-IRS-1 and PDZK1 pathways regulating cancer cell invasion, survival, metastasis, angiogenesis and proliferation. The signaling alterations of MCF-7 cells in turn feed into multiple TSGs including BRCA1, p53, Mdm2 and PTEN which are activated as a result of activated signaling. There is a strong correlation of PDZK1 with multiple signaling pathways including estrogen dependent IGF-1R and chemokine (C-X-C motif) receptor 4 (CXCR4) signaling which can be 22 / 29 targeted therapeutically to treat breast cancer [ 24, 26, 124 ]. It has been reported that the overexpression of PDZK1 was associated with resistance to paclitaxel-5-fluorouracil-etoposide at low concentration [ 24 ]. In this study, we reported for the first time that fulvestrant exhibited significant inhibitory activity against PDZK1. More importantly, our results demonstrate a novel function for PDZK1 and find a molecular crosstalk between biological regulatory growth factors that are involved in breast cancer signaling. Although, we tested fulvestrant in-silico and in-vitro in MCF-7 cells, same biological results were observed to analyse the deregulated expression of growth promoting factors (Table 3). It shows the inhibitory effect of fulvestrant on IGF-1R, IRS-1, ER-α and PDZK1 at low concentration. We believe, that fulvestrant may exhibit as an effective and safe anti-cancer drug against multiple targets specially IGF-1R but further experimental validation is needed using inhibition assay against this receptor's activity and its downstream signaling molecules which should be elucidated in the future studies. Supporting information S1 Table. Pharmacophore model was validated by the significance of statistical parameters with respect to molecular sequence and root mean square deviation (RMSD). It shows the RMSD of 21 set of active compounds with respect to IC50 values were generated to build the pharmacophore model. The low RMSD values have ability to predict the activity of the conformational dataset compounds. (DOCX) Acknowledgments We thankful to Dr. Felicity May from Northern Institute for Cancer Research, Newcastle upon Tyne Hospitals NHS Foundation Trust, The Medical School, University of Newcastle upon Tyne, Framlington Place, Newcastle upon Tyne, UK for providing me research lab facilities for this work. Data curation: Samra Khalid, Qaisar Mansoor. Formal analysis: Samra Khalid. Investigation: Samra Khalid. Methodology: Samra Khalid. Project administration: Samra Khalid. Resources: Samra Khalid, Ishrat Jabeen, Qaisar Mansoor, Muhammad Ismail. Software: Ishrat Jabeen. Supervision: Rumeza Hanif, Ishrat Jabeen, Muhammad Ismail. Validation: Samra Khalid, Qaisar Mansoor. Visualization: Samra Khalid, Qaisar Mansoor. Writing ± original draft: Samra Khalid. Writing ± review & editing: Samra Khalid. 23 / 29 24 / 29 22. Watters JJ, Chun T-Y, Kim Y-N, Bertics PJ, Gorski J. Estrogen modulation of prolactin gene expression requires an intact mitogen-activated protein kinase signal transduction pathway in cultured rat pituitary cells. Molecular Endocrinology. 2000; 14(11):1872±81. https://doi.org/10.1210/mend.14.11. 0551 PMID: 11075818 25 / 29 breast cancer cells. Breast cancer research and treatment. 2005; 94(1):37±46. https://doi.org/10. 1007/s10549-005-6939-z PMID: 16142439 26 / 29 27 / 29 Korean women. PloS one. 2014; 9(1):e84532. https://doi.org/10.1371/journal.pone.0084532 PMID: 24392142 28 / 29 1. LeRoith D , Werner H , Beitner-Johnson D , ROBERTS JR AT. Molecular and cellular aspects of the insulin-like growth factor I receptor . Endocrine reviews. 1995 ; 16 ( 2 ): 143 ± 63 . https://doi.org/10.1210/ edrv-16-2-143 PMID: 7540132 2. Bhaumick B , Bala R , Hollenberg M. Somatomedin receptor of human placenta: solubilization, photolabeling, partial purification, and comparison with insulin receptor . Proceedings of the National Academy of Sciences . 1981 ; 78 ( 7 ): 4279 ± 83 . 3. Chernausek SD , Jacobs S , Van Wyk JJ. Structural similarities between human receptors for somatomedin C and insulin: analysis by affinity labeling . Biochemistry . 1981 ; 20 ( 26 ): 7345 ± 50 . PMID: 6275879 4. Pandini G , Vigneri R , Costantino A , Frasca F , Ippolito A , Fujita-Yamaguchi Y , et al. Insulin and insulinlike growth factor-I (IGF-I) receptor overexpression in breast cancers leads to insulin/IGF-I hybrid receptor overexpression: evidence for a second mechanism of IGF-I signaling . Clinical Cancer Research . 1999 ; 5 ( 7 ): 1935 ± 44 . PMID: 10430101 5. Adams TE , Epa V , Garrett T , Ward C . Structure and function of the type 1 insulin-like growth factor receptor . Cellular and Molecular Life Sciences CMLS . 2000 ; 57 ( 7 ): 1050 ± 93 . https://doi.org/10.1007/ PL00000744 PMID: 10961344 6. Fagan DH , Yee D . Crosstalk between IGF1R and estrogen receptor signaling in breast cancer . Journal of mammary gland biology and neoplasia . 2008 ; 13 ( 4 ): 423 ±9. https://doi.org/10.1007/s10911-008 - 9098-0 PMID: 19003523 7. Cullen KJ , Yee D , Sly WS , Perdue J , Hampton B , Lippman ME , et al. Insulin-like growth factor receptor expression and function in human breast cancer . Cancer research . 1990 ; 50 ( 1 ): 48 ± 53 . PMID: 2152773 8. Gooch JL , Van Den Berg CL , Yee D . Insulin-like growth factor (IGF)-I rescues breast cancer cells from chemotherapy-induced cell death±proliferative and anti-apoptotic effects . Breast cancer research and treatment. 1999 ; 56 ( 1 ):1± 10 . PMID: 10517338 9. Lee A , Yee D . Insulin-like growth factors and breast cancer . Biomedicine & pharmacotherapy . 1995 ; 49 ( 9 ): 415 ± 21 . 10. Peyrat JP , Bonneterre J . Type 1 IGF receptor in human breast diseases . Breast cancer research and treatment. 1992 ; 22 ( 1 ): 59 ± 67 . PMID: 1421425 11. Nickerson T , Chang F , Lorimer D , Smeekens SP , Sawyers CL , Pollak M. In vivo progression of LAPC9 and LNCaP prostate cancer models to androgen independence is associated with increased expression of insulin-like growth factor I (IGF-I) and IGF-I receptor (IGF-IR) . Cancer research . 2001 ; 61 ( 16 ): 6276 ± 80 . PMID: 11507082 12. Ge N-L , Rudikoff S. Insulin-like growth factor I is a dual effector of multiple myeloma cell growth . Blood . 2000 ; 96 ( 8 ): 2856 ± 61 . PMID: 11023522 13. Hassan A , Macaulay V. The insulin-like growth factor system as a therapeutic target in colorectal cancer . Annals of Oncology . 2002 ; 13 ( 3 ): 349 ± 56 . PMID: 11996463 14. Wu Y , Yakar S , Zhao L , Hennighausen L , LeRoith D. Circulating insulin -like growth factor-I levels regulate colon cancer growth and metastasis . Cancer research . 2002 ; 62 ( 4 ): 1030 ± 5 . PMID: 11861378 15. Farabaugh SM , Boone DN , Lee AV . Role of IGF1R in breast cancer subtypes, stemness, and lineage differentiation . Frontiers in endocrinology. 2015 ; 6 : 59 . https://doi.org/10.3389/fendo. 2015 .00059 PMID: 25964777 16. Puzanov I , Lindsay CR , Goff L , Sosman J , Gilbert J , Berlin J , et al. A phase I study of continuous oral dosing of OSI-906, a dual inhibitor of insulin-like growth factor-1 and insulin receptors, in patients with advanced solid tumors . Clinical cancer research . 2014 . 17. Jones RL , Kim ES , Nava-Parada P , Alam S , Johnson FM , Stephens A , et al. Phase I study of intermittent oral dosing of the insulin-like growth factor-1 and insulin receptors inhibitor OSI-906 in patients with advanced solid tumors . Clinical cancer research . 2014:clincanres. 0265 . 2014 . 18. Guha M. Anticancer IGF1R classes take more knocks . Nature Publishing Group; 2013 . 19. Law JH , Habibi G , Hu K , Masoudi H , Wang MY , Stratford AL , et al. Phosphorylated insulin-like growth factor-i/insulin receptor is present in all breast cancer subtypes and is related to poor survival . Cancer research . 2008 ; 68 ( 24 ): 10238 ± 46 . https://doi.org/10.1158/ 0008 - 5472 .CAN- 08 -2755 PMID: 19074892 20. Khalid S , Hanif R , Tareen SH , Siddiqa A , Bibi Z , Ahmad J . Formal modeling and analysis of ER-α associated Biological Regulatory Network in breast cancer . PeerJ . 2016 ; 4:e2542 . https://doi.org/10.7717/ peerj.2542 PMID: 27781158 21. Pollak M. Insulin and insulin-like growth factor signalling in neoplasia . Nature Reviews Cancer . 2008 ; 8 ( 12 ): 915 ± 28 . https://doi.org/10.1038/nrc2536 PMID: 19029956 23. Vivanco I , Sawyers CL . The phosphatidylinositol 3-kinase±AKT pathway in human cancer . Nature Reviews Cancer . 2002 ; 2 ( 7 ): 489 ± 501 . https://doi.org/10.1038/nrc839 PMID: 12094235 24. Kim H , Abd Elmageed ZY , Davis C , El-Bahrawy AH , Naura AS , Ekaidi I , et al. Correlation between PDZK1, Cdc37, Akt and breast cancer malignancy: the role of PDZK1 in cell growth through Akt stabilization by increasing and interacting with Cdc37 . Mol Med . 2014 ; 20 : 270 ±9. https://doi.org/10.2119/ molmed. 2013 .00166 PMID: 24869908 25. Kim JK , Kwon O , Kim J , Kim E -K, Park HK , Lee JE , et al. PDZ domain-containing 1 (PDZK1) protein regulates phospholipase C-β3 ( PLC-β3 ) -specific activation of somatostatin by forming a ternary complex with PLC-β3 and somatostatin receptors . Journal of Biological Chemistry . 2012 ; 287 ( 25 ): 21012 ± 24 . https://doi.org/10.1074/jbc. M111.337865 PMID: 22528496 26. Kim H , Elmageed ZYA , Ju J , Naura AS , Abdel-Mageed AB , Varughese S , et al. PDZK1 is a novel factor in breast cancer that is indirectly regulated by estrogen through IGF-1R and promotes estrogenmediated growth . Molecular Medicine . 2013 ; 19 ( 1 ): 253 . 27. Shelton JG , Steelman LS , White ER , McCubrey JA . Synergy between PI3K/Akt and Raf/MEK/ERK pathways in IGF-1R mediated cell cycle progression and prevention of apoptosis in hematopoietic cells . Cell Cycle . 2004 ; 3 ( 3 ): 370 ± 7 . 28. Ehrlich P. UÈber die constitution des diphtheriegiftes . DMW-Deutsche Medizinische Wochenschrift . 1898 ; 24 ( 38 ): 597 ± 600 . 29. Bajorath J . Selected concepts and investigations in compound classification, molecular descriptor analysis, and virtual screening . Journal of chemical information and computer sciences . 2001 ; 41 ( 2 ): 233 ± 45 . PMID: 11277704 30. Alla H , David R . Continuous and hybrid Petri nets . Journal of Circuits, Systems, and Computers . 1998 ; 8 ( 01 ): 159 ± 88 . 32. Wolber G , Seidel T , Bendix F , Langer T . Molecule-pharmacophore superpositioning and pattern matching in computational drug design . Drug discovery today . 2008 ; 13 ( 1 ): 23 ± 9 . Wermuth C-G , Ganellin CR , Lindberg P , Mitscher LA . Glossary of terms used in medicinal chemistry (IUPAC Recommendations 1997 ). Annual reports in medicinal chemistry . 1998 ; 33 : 385 ± 95 . 33. Yang S-Y. Pharmacophore modeling and applications in drug discovery: challenges and recent advances . Drug discovery today . 2010 ; 15 ( 11 ): 444 ± 50 . 34. Gao Q , Yang L , Zhu Y. Pharmacophore based drug design approach as a practical process in drug discovery. Current computer-aided drug design . 2010 ; 6 ( 1 ): 37 ± 49 . PMID: 20370694 35. David R , editor Modeling of hybrid systems using continuous and hybrid Petri nets . Petri Nets and Performance Models , 1997 , Proceedings of the Seventh International Workshop on; 1997: IEEE. 36. Chenna R , Sugawara H , Koike T , Lopez R , Gibson TJ , Higgins DG , et al. Multiple sequence alignment with the Clustal series of programs . Nucleic acids research . 2003 ; 31 ( 13 ): 3497 ± 500 . PMID: 12824352 37. Girnita A , Girnita L , del Prete F , Bartolazzi A , Larsson O , Axelson M . Cyclolignans as inhibitors of the insulin-like growth factor-1 receptor and malignant cell growth . Cancer research . 2004 ; 64 ( 1 ): 236 ± 42 . PMID: 14729630 38. Bumke-Vogt C , Osterhoff MA , Borchert A , Guzman-Perez V , Sarem Z , Birkenfeld AL , et al. The flavones apigenin and luteolin induce FOXO1 translocation but inhibit gluconeogenic and lipogenic gene expression in human cells . PloS one . 2014 ; 9 ( 8 ):e104321. https://doi.org/10.1371/journal.pone. 0104321 PMID: 25136826 39. Li M , He Z , Ermakova S , Zheng D , Tang F , Cho Y-Y , et al. Direct inhibition of insulin-like growth factorI receptor kinase activity by (− ) − epigallocatechin - 3 -gallate regulates cell transformation . Cancer Epidemiology Biomarkers & Prevention . 2007 ; 16 ( 3 ): 598 ± 605 . 40. Marconett CN , Singhal AK , Sundar SN , Firestone GL . Indole-3-carbinol disrupts estrogen receptoralpha dependent expression of insulin-like growth factor-1 receptor and insulin receptor substrate-1 and proliferation of human breast cancer cells . Molecular and cellular endocrinology . 2012 ; 363 ( 1 ): 74 ± 84 . 41. Stahl P , Kissau L , Mazitschek R , Giannis A , Waldmann H. Natural Product Derived Receptor Tyrosine Kinase Inhibitors: Identification of IGF1R, Tie-2, and VEGFR-3 Inhibitors. Angewandte Chemie International Edition . 2002 ; 41 ( 7 ): 1174 ± 8 . PMID: 12491250 42. Youngren JF , Gable K , Penaranda C , Maddux BA , Zavodovskaya M , Lobo M , et al. Nordihydroguaiaretic acid (NDGA) inhibits the IGF-1 and c-erbB2/HER2/neu receptors and suppresses growth in 43. Chen J , Hou R , Zhang X , Ye Y , Wang Y , Tian J . Calycosin suppresses breast cancer cell growth via ERβ-dependent regulation of IGF-1R, p38 MAPK and PI3K/Akt pathways . PloS one . 2014 ; 9 ( 3 ): e91245. https://doi.org/10.1371/journal.pone. 0091245 PMID: 24618835 44. Gable KL , Maddux BA , Penaranda C , Zavodovskaya M , Campbell MJ , Lobo M , et al. Diarylureas are small-molecule inhibitors of insulin-like growth factor I receptor signaling and breast cancer cell growth . Molecular cancer therapeutics . 2006 ; 5 ( 4 ): 1079 ± 86 . https://doi.org/10.1158/ 1535 - 7163 .MCT- 05 -0397 PMID: 16648580 45. Anderson MO , Yu H , Penaranda C , Maddux BA , Goldfine ID , Youngren JF , et al. Parallel synthesis of diarylureas and their evaluation as inhibitors of insulin-like growth factor receptor . Journal of combinatorial chemistry . 2006 ; 8 ( 5 ): 784 ± 90 . https://doi.org/10.1021/cc050136z PMID: 16961415 46. Chen L-H , Fang J , Sun Z , Li H , Wu Y , Demark-Wahnefried W , et al. Enterolactone inhibits insulin-like growth factor-1 receptor signaling in human prostatic carcinoma PC-3 cells . The Journal of nutrition . 2009 ; 139 ( 4 ): 653 ±9. https://doi.org/10.3945/jn.108.101832 PMID: 19211828 47. Ryan CJ , Harzstark AH , Rosenberg J , Lin A , Claros C , Goldfine ID , et al. A pilot dose-escalation study of the effects of nordihydroguareacetic acid on hormone and prostate specific antigen levels in patients with relapsed prostate cancer . BJU international . 2008 ; 101 ( 4 ): 436 ±9. https://doi.org/10.1111/j. 1464 - 410X . 2007 . 07330 . x PMID : 18234062 48. Anastasius N , Boston S , Lacey M , Storing N , Whitehead SA . Evidence that low-dose, long-term genistein treatment inhibits oestradiol-stimulated growth in MCF-7 cells by down-regulation of the PI3- kinase/Akt signalling pathway . The Journal of steroid biochemistry and molecular biology . 2009 ; 116 ( 1 ): 50 ± 5 . 49. Kang J-H , Song K-H , Woo J-K, Park MH , Rhee MH , Choi C , et al. Ginsenoside Rp1 from Panax ginseng exhibits anti-cancer activity by down-regulation of the IGF-1R/Akt pathway in breast cancer cells . Plant Foods for Human Nutrition . 2011 ; 66 ( 3 ): 298 ± 305 . https://doi.org/10.1007/s11130-011 -0242-4 PMID: 21748437 50. Xia Y , Jin L , Zhang B , Xue H , Li Q , Xu Y. The potentiation of curcumin on insulin-like growth factor-1 action in MCF-7 human breast carcinoma cells . Life sciences . 2007 ; 80 ( 23 ): 2161 ±9. https://doi.org/10. 1016/j.lfs. 2007 . 04 .008 PMID: 17499312 51. El Amrani M , Lai D , Debbab A , Aly AH , Siems K , Seidel C , et al. Protein kinase and HDAC inhibitors from the endophytic fungus Epicoccum nigrum . Journal of natural products . 2013 ; 77 ( 1 ): 49 ± 56 . https://doi.org/10.1021/np4005745 PMID: 24328302 52. Tong L-j, Xie H , Peng T , Liu X-f , Xin X-l, Huang X , et al. Establishment of platform for screening insulin-like growth factor-1 receptor inhibitors and evaluation of novel inhibitors . Acta Pharmacologica Sinica . 2011 ; 32 ( 7 ): 930 ±8. https://doi.org/10.1038/aps. 2011 .23 PMID: 21643004 53. Mendelsohn LD . ChemDraw 8 ultra, windows and macintosh versions . Journal of chemical information and computer sciences . 2004 ; 44 ( 6 ): 2225 ± 6 . 54. Labute P. Molecular Operating Environment . Chemical Computing Group. Inc, Montreal. 2008 . 55. Chen I-J , Foloppe N. Conformational sampling of druglike molecules with MOE and catalyst: implications for pharmacophore modeling and virtual screening . Journal of chemical information and modeling . 2008 ; 48 ( 9 ): 1773 ± 91 . https://doi.org/10.1021/ci800130k PMID: 18763758 56. Lipinski CA , Lombardo F , Dominy BW , Feeney PJ . Experimental and computational approaches to estimate solubility and permeability in drug discovery and development settings . Advanced drug delivery reviews . 2012 ; 64 :4± 17 . 57. Ertl P , Rohde B , Selzer P . Fast calculation of molecular polar surface area as a sum of fragmentbased contributions and its application to the prediction of drug transport properties . Journal of medicinal chemistry . 2000 ; 43 ( 20 ): 3714 ± 7 . PMID: 11020286 58. Veber DF , Johnson SR , Cheng H -Y, Smith BR , Ward KW , Kopple KD . Molecular properties that influence the oral bioavailability of drug candidates . Journal of medicinal chemistry . 2002 ; 45 ( 12 ): 2615 ± 23 . PMID: 12036371 59. Keller TH , Pichota A , Yin Z . A practical view of `druggability' . Current opinion in chemical biology . 2006 ; 10 ( 4 ): 357 ± 61 . https://doi.org/10.1016/j.cbpa. 2006 . 06 .014 PMID: 16814592 60. Heiner M , Herajy M , Liu F , Rohr C , Schwarick M , editors. Snoopy± a unifying Petri net tool . International Conference on Application and Theory of Petri Nets and Concurrency ; 2012: Springer. 61. Ruths D , Muller M , Tseng J-T, Nakhleh L , Ram PT . The signaling petri net-based simulator: a nonparametric strategy for characterizing the dynamics of cell-specific signaling networks . PLoS Comput Biol . 2008 ; 4 ( 2 ):e1000005. https://doi.org/10.1371/journal.pcbi. 1000005 PMID: 18463702 62. Network CGA . Comprehensive molecular portraits of human breast tumours . Nature . 2012 ; 490 ( 7418 ): 61 ± 70 . https://doi.org/10.1038/nature11412 PMID: 23000897 63. Sotiriou C , Neo S-Y , McShane LM , Korn EL , Long PM , Jazaeri A , et al. Breast cancer classification and prognosis based on gene expression profiles from a population-based study . Proceedings of the National Academy of Sciences . 2003 ; 100 ( 18 ): 10393 ± 8 . 64. Kang H , Yi Y , Kim H , Hong Y , Seong Y , Bae I. BRCA1 negatively regulates IGF-1 expression through an estrogen-responsive element-like site . Cell death & disease . 2012 ; 3 ( 6 ): e336 . 65. Bailey ST , Shin H , Westerling T , Liu XS , Brown M. Estrogen receptor prevents p53-dependent apoptosis in breast cancer . Proceedings of the National Academy of Sciences . 2012 ; 109 ( 44 ): 18060 ± 5 . 66. Caldon CE . Estrogen signaling and the DNA damage response in hormone dependent breast cancers . Molecular mechanisms of cellular stress responses in cancer and their therapeutic implications . 2007 : 137 . 67. Koczorowska MM , KWAŚNIEWSKA A , GO ZÂDZICKA-JOÂ ZEFIAK A. Expression of IGF1, IGF1R, and SP1 factors in HPV infected cervical cancer cells . Bull Vet Inst Pulawy . 2011 ; 55 : 281 ± 6 . 68. Inoue J , Otsuki T , Hirasawa A , Imoto I , Matsuo Y , Shimizu S , et al. Overexpression of PDZK1 within the 1q12-q22 amplicon is likely to be associated with drug-resistance phenotype in multiple myeloma . The American journal of pathology . 2004 ; 165 ( 1 ): 71 ± 81 . https://doi.org/10.1016/S0002- 9440 ( 10 ) 63276 - 2 PMID: 15215163 69. BIECHE I , PARFAIT B , LAURENDEAU I , GIRAULT I , VIDAUD M , LIDEREAU R. Quantification of estrogen receptor α and β expression in sporadic breast cancer . Oncogene . 2001 ; 20 ( 56 ): 8109 ± 15 . https://doi.org/10.1038/sj.onc. 1204917 PMID: 11781824 70. Hanif R . Molecular mechanism of regulation of iron transport across placenta : UCL (University College London); 2012 . 71. Davison Z , de Blacquière GE , Westley BR , May FE . Insulin-like growth factor-dependent proliferation and survival of triple-negative breast cancer cells: implications for therapy . Neoplasia . 2011 ; 13 ( 6 ): 504 ± 15 . PMID: 21677874 72. Molloy CA , May FE , Westley BR . Insulin receptor substrate-1 expression is regulated by estrogen in the MCF-7 human breast cancer cell line . Journal of Biological Chemistry . 2000 ; 275 ( 17 ): 12565 ± 71 . PMID: 10777546 73. Christopoulos PF , Msaouel P , Koutsilieris M. The role of the insulin-like growth factor-1 system in breast cancer . Molecular cancer . 2015 ; 14 ( 1 ): 1 . 74. Yerushalmi R , Gelmon KA , Leung S , Gao D , Cheang M , Pollak M , et al. Insulin-like growth factor receptor (IGF-1R) in breast cancer subtypes . Breast cancer research and treatment. 2012 ; 132 ( 1 ): 131 ± 42 . https://doi.org/10.1007/s10549-011 -1529-8 PMID: 21574055 75. Triballeau N , Acher F , Brabet I , Pin J-P , Bertrand H-O . Virtual screening workflow development guided by the ªreceiver operating characteristicº curve approach. Application to high-throughput docking on metabotropic glutamate receptor subtype 4 . Journal of medicinal chemistry . 2005 ; 48 ( 7 ): 2534 ± 47 . https://doi.org/10.1021/jm049092j PMID: 15801843 76. Law V , Knox C , Djoumbou Y , Jewison T , Guo AC , Liu Y , et al. DrugBank 4 . 0: shedding new light on drug metabolism . Nucleic acids research . 2013 ; 42 ( D1 ): D1091 ± D7 . 77. Maddox DD , Rahman Z. Etonogestrel (Implanon), another treatment option for contraception . Pharmacy and Therapeutics . 2008 ; 33 ( 6 ): 337 . 78. Barkfeldt J , Virkkunen A , Dieben T. The effects of two progestogen-only pills containing either desogestrel (75 μg/day) or levonorgestrel (30 μg/day) on lipid metabolism . Contraception . 2001 ; 64 ( 5 ): 295 ± 9 . PMID: 11777489 79. Ciruelos E , Pascual T , Vozmediano MLA , Blanco M , Manso L , Parrilla L , et al. The therapeutic role of fulvestrant in the management of patients with hormone receptor-positive breast cancer . The Breast . 2014 ; 23 ( 3 ): 201 ±8. https://doi.org/10.1016/j.breast. 2014 . 01 .016 PMID: 24589524 80. Peters K , Peters E-M , Panitzsch T , Tochtermann W. Crystal structure of (3aR*, 6aS* , 9aS*, 9bR*)- decahydro-8 , 8 - dimethyl- 6 , 9 -dioxoazuleno [4, 5-c] furan - 1 ( 3H ) -one,(CH3) 2C12H12O4 . Zeitschrift fuÈr Kristallographie-New Crystal Structures . 1999 ; 214 ( 2 ): 269 ± 70 . 81. Shaefer CF . Lixisenatide: a new member of the Glucagon-Like peptide 1 receptor agonist class of incretin therapies . Clinical Diabetes . 2016 ; 34 ( 2 ): 81 ±5. https://doi.org/10.2337/diaclin.34.2.81 PMID: 27092017 82. Wasilewski A , Lewandowska U , Mosinska P , Watala C , Storr M , Fichna J , et al. Cannabinoid receptor type 1 and mu-opioid receptor polymorphisms are associated with cyclic vomiting syndrome . The American journal of gastroenterology . 2017 . 83. Kang H-S , Ahn SH , Mishra SK , Hong K-M , Lee ES , Shin KH , et al. Association of polymorphisms and haplotypes in the insulin-like growth factor 1 receptor (IGF1R) gene with the risk of breast cancer in 84. Alluri PG , Speers C , Chinnaiyan AM . Estrogen receptor mutations and their role in breast cancer progression . Breast Cancer Research . 2014 ; 16 ( 6 ): 1 . 85. Angeloni S , Martin M , Garcia-Morales P , Castro-Galache M , Ferragut J , Saceda M. Regulation of estrogen receptor-alpha expression by the tumor suppressor gene p53 in MCF-7 cells . Journal of endocrinology . 2004 ; 180 ( 3 ): 497 ± 504 . PMID: 15012604 86. Kim K , Burghardt R , Barhoumi R , Lee S-o , Liu X , Safe S. MDM2 regulates estrogen receptor α and estrogen responsiveness in breast cancer cells . Journal of molecular endocrinology . 2011 ; 46 ( 2 ): 67 ± 79 . https://doi.org/10.1677/JME-10-0110 PMID: 21169420 87. Lee Y-R , Park J , Yu H-N , Kim J-S , Youn HJ , Jung SH . Up-regulation of PI3K/Akt signaling by 17β- estradiol through activation of estrogen receptor-α, but not estrogen receptor-β, and stimulates cell growth in breast cancer cells . Biochemical and biophysical research communications . 2005 ; 336 ( 4 ): 1221 ±6. https://doi.org/10.1016/j.bbrc. 2005 . 08 .256 PMID: 16169518 88. Di Leo A , Jerusalem G , Petruzelka L , Torres R , Bondarenko IN , Khasanov R , et al. Results of the CONFIRM phase III trial comparing fulvestrant 250 mg with fulvestrant 500 mg in postmenopausal women with estrogen receptor±positive advanced breast cancer . Journal of Clinical Oncology . 2010 ; 28 ( 30 ): 4594 ± 600 . https://doi.org/10.1200/JCO. 2010 . 28 .8415 PMID: 20855825 89. Park S-H , Cheung LW , Wong AS , Leung PC . Estrogen regulates Snail and Slug in the down-regulation of E-cadherin and induces metastatic potential of ovarian cancer cells through estrogen receptor α . Molecular Endocrinology . 2008 ; 22 ( 9 ): 2085 ±98. https://doi.org/10.1210/me.2007-0512 PMID: 18550773 90. Stabile LP , Davis ALG , Gubish CT , Hopkins TM , Luketich JD , Christie N , et al. Human non-small cell lung tumors and cells derived from normal lung express both estrogen receptor α and β and show biological responses to estrogen . Cancer research . 2002 ; 62 ( 7 ): 2141 ± 50 . PMID: 11929836 91. Kameda C , Nakamura M , Tanaka H , Yamasaki A , Kubo M , Tanaka M , et al. Oestrogen receptor-α contributes to the regulation of the hedgehog signalling pathway in ERα-positive gastric cancer . British journal of cancer . 2010 ; 102 ( 4 ): 738 ± 47 . https://doi.org/10.1038/sj.bjc. 6605517 PMID: 20087349 92. Georgescu M-M. PTEN tumor suppressor network in PI3K-Akt pathway control . Genes & cancer . 2010 ; 1 ( 12 ): 1170 ± 7 . 93. Capoulade C , Mir LM , Carlier K , LeÂcluse Y , TeÂtaud C , Mishal Z , et al. Apoptosis of tumoral and nontumoral lymphoid cells is induced by both mdm2 and p53 antisense oligodeoxynucleotides . Blood . 2001 ; 97 ( 4 ): 1043 ± 9 . PMID: 11159535 94. Nakanishi A , Kitagishi Y , Ogura Y , Matsuda S. The tumor suppressor PTEN interacts with p53 in hereditary cancer (Review) . International journal of oncology . 2014 ; 44 ( 6 ): 1813 ± 9 . https://doi.org/10. 3892/ijo. 2014 .2377 PMID: 24718924 95. Mayo LD , Donner DB . The PTEN , Mdm2, p53 tumor suppressor±oncoprotein network . Trends in biochemical sciences . 2002 ; 27 ( 9 ): 462 ± 7 . PMID: 12217521 96. Berger CE , Qian Y , Liu G , Chen H , Chen X. p53 , a target of estrogen receptor (ER) α, modulates DNA damage-induced growth suppression in ER-positive breast cancer cells . Journal of Biological Chemistry . 2012 ; 287 ( 36 ): 30117 ± 27 . https://doi.org/10.1074/jbc. M112.367326 PMID: 22787161 97. Ma Y , Fan S , Hu C , Meng Q , Fuqua SA , Pestell RG , et al. BRCA1 regulates acetylation and ubiquitination of estrogen receptor-α . Molecular Endocrinology . 2010 ; 24 ( 1 ): 76 ± 90 . https://doi.org/10.1210/me. 2009-0218 PMID: 19887647 98. Renoir J-M , Marsaud V , Lazennec G . Estrogen receptor signaling as a target for novel breast cancer therapeutics . Biochemical pharmacology . 2013 ; 85 ( 4 ): 449 ± 65 . https://doi.org/10.1016/j.bcp. 2012 . 10 . 018 PMID: 23103568 99. Ariazi EA , Ariazi JL , Cordera F , Jordan VC . Estrogen receptors as therapeutic targets in breast cancer . Current topics in medicinal chemistry . 2006 ; 6 ( 3 ): 181 ± 202 . PMID: 16515478 100. Abramovitch S , Werner H. Functional and physical interactions between BRCA1 and p53 in transcriptional regulation of the IGF-IR gene . Hormone and metabolic research . 2003 ; 35 ( 11 /12): 758 ± 62 . 101. McGowan EM , Alling N , Jackson EA , Yagoub D , Haass NK , Allen JD , et al. Evaluation of cell cycle arrest in estrogen responsive MCF-7 breast cancer cells: pitfalls of the MTS assay . PloS one . 2011 ; 6 ( 6 ):e20623. https://doi.org/10.1371/journal.pone. 0020623 PMID: 21673993 102. Doisneau-Sixou S , Sergio C , Carroll J , Hui R , Musgrove E , Sutherland R . Estrogen and antiestrogen regulation of cell cycle progression in breast cancer cells . Endocrine-related cancer . 2003 ; 10 ( 2 ): 179 ± 86 . PMID: 12790780 103. Liu S , Meng X , Chen H , Liu W , Miller T , Murph M , et al. Targeting tyrosine-kinases and estrogen receptor abrogates resistance to endocrine therapy in breast cancer . Oncotarget . 2014 ; 5 ( 19 ): 9049 . https:// doi.org/10.18632/oncotarget.2022 PMID: 24979294 104. Juric D , Rodon J , Gonzalez-Angulo A , Burris H , Bendell J , Berlin J , et al. BYL719, a next generation PI3K alpha specific inhibitor: Preliminary safety . PK, and efficacy results from the first-in-human study . 2012 . 105. Dearth RK , Cui X , Kim H-J , Hadsell DL , Lee AV . Oncogenic transformation by the signaling adaptor proteins insulin receptor substrate (IRS)-1 and IRS-2 . Cell Cycle . 2007 ; 6 ( 6 ): 705 ± 13 . https://doi.org/ 10.4161/cc.6.6.4035 PMID: 17374994 106. Porter HA , Perry A , Kingsley C , Tran NL , Keegan AD . IRS1 is highly expressed in localized breast tumors and regulates the sensitivity of breast cancer cells to chemotherapy, while IRS2 is highly expressed in invasive breast tumors . Cancer letters . 2013 ; 338 ( 2 ): 239 ± 48 . https://doi.org/10.1016/j. canlet. 2013 . 03 .030 PMID: 23562473 107. Mardilovich K , Pankratz SL , Shaw LM . Expression and function of the insulin receptor substrate proteins in cancer . Cell Communication and Signaling . 2009 ; 7 ( 1 ): 1 . 108. Schnarr B , Strunz K , Ohsam J , Benner A , Wacker J , Mayer D . Down-regulation of insulin-like growth factor-I receptor and insulin receptor substrate-1 expression in advanced human breast cancer . International journal of cancer . 2000 ; 89 ( 6 ): 506 ± 13 . PMID: 11102895 109. Sisci D , Morelli C , Garofalo C , Romeo F , Morabito L , Casaburi F , et al. Expression of nuclear insulin receptor substrate 1 in breast cancer . Journal of clinical pathology . 2007 ; 60 ( 6 ): 633 ± 41 . https://doi. org/10.1136/jcp. 2006 .039107 PMID: 16882697 110. Siegel RL , Miller KD , Jemal A . Cancer statistics, 2016 . CA: a cancer journal for clinicians . 2016 ; 66 ( 1 ):7± 30 . 111. Vargo-Gogola T , Rosen JM . Modelling breast cancer: one size does not fit all . Nat Rev Cancer . 2007 ; 7 ( 9 ): 659 ± 72 . Epub 2007/08/28. nrc2193 [pii] https://doi.org/10.1038/nrc2193 PMID: 17721431 . 112. Macgregor JI , Jordan VC . Basic guide to the mechanisms of antiestrogen action . Pharmacological reviews . 1998 ; 50 ( 2 ): 151 ± 96 . PMID: 9647865 113. Zhao M , Ramaswamy B . Mechanisms and therapeutic advances in the management of endocrineresistant breast cancer . World J Clin Oncol . 2014 ; 5 ( 3 ): 248 ± 62 . https://doi.org/10.5306/wjco. v5.i3.248 PMID: 25114842 114. Knowlden JM , Hutcheson IR , Barrow D , Gee JMW , Nicholson RI . Insulin-like growth factor-I receptor signaling in tamoxifen-resistant breast cancer: a supporting role to the epidermal growth factor receptor . Endocrinology . 2005 ; 146 ( 11 ): 4609 ± 18 . https://doi.org/10.1210/en.2005-0247 PMID: 16037379 115. Van Asten K , Neven P , Lintermans A , Wildiers H , Paridaens R . Aromatase inhibitors in the breast cancer clinic: focus on exemestane . Endocrine-related cancer . 2014 ; 21 ( 1):R31±R49 . https://doi.org/10. 1530/ERC-13-0269 PMID: 24434719 116. Biswas A , Viegas OA , Bennink HJC , Korver T , Ratnam SS . Implanon® contraceptive implants: effects on carbohydrate metabolism . Contraception . 2001 ; 63 ( 3 ): 137 ± 41 . PMID: 11368985 117. Croxatto HB . Mechanisms that explain the contraceptive action of progestin implants for women . Contraception . 2002 ; 65 ( 1 ): 21 ± 7 . PMID: 11861052 118. Banhegyi R , Laczo I , Fulop F , Mellar E , Piko B . The Role of Fulvestrant in the Treatment of Metastatic Breast Cancer: A Case Report . Journal of Steroids & Hormonal Science . 2014 ; 2014 . 119. Hamilton DH , Griner LM , Keller JM , Hu X , Southall N , Marugan J , et al. Targeting estrogen receptor signaling with fulvestrant enhances immune and chemotherapy-mediated cytotoxicity of human lung cancer . Clinical Cancer Research . 2016:clincanres. 3059 . 2015 . 120. Brauer W , Reisig W , Rozenberg G. Petri Nets: Central Models and Their Properties: Advances in Petri Nets 1986 , Part I Proceedings of an Advanced Course Bad Honnef , 8 .±19. September 1986: Springer; 2006 . 121. Liao XH , Lu DL , Wang N , Liu LY , Wang Y , Li YQ , et al. Estrogen receptor α mediates proliferation of breast cancer MCF±7 cells via a p21/PCNA/E2F1-dependent pathway . FEBS Journal . 2014 ; 281 ( 3 ): 927 ± 42 . https://doi.org/10.1111/febs.12658 PMID: 24283290 122. Beeram M , Tan Q-T , Tekmal R , Russell D , Middleton A . Akt-induced endocrine therapy resistance is reversed by inhibition of mTOR signaling . Annals of Oncology . 2007 ; 18 ( 8 ): 1323 ±8. https://doi.org/10. 1093/annonc/mdm170 PMID: 17693645 123. Kanehisa M , Goto S. KEGG : kyoto encyclopedia of genes and genomes . Nucleic acids research . 2000 ; 28 ( 1 ): 27 ± 30 . PMID: 10592173 124. Khalid S , Hanif R . Association of rs1801157 single nucleotide polymorphism of CXCL12 gene in breast cancer in Pakistan and in-silico expression analysis of CXCL12±CXCR4 associated biological regulatory network . PeerJ . 2017 ; 5:e3822 . https://doi.org/10.7717/peerj.3822 PMID: 28929029


This is a preview of a remote PDF: http://journals.plos.org/plosone/article/file?id=10.1371/journal.pone.0196312&type=printable

Samra Khalid, Rumeza Hanif, Ishrat Jabeen, Qaisar Mansoor, Muhammad Ismail. Pharmacophore modeling for identification of anti-IGF-1R drugs and in-vitro validation of fulvestrant as a potential inhibitor, PLOS ONE, 2018, DOI: 10.1371/journal.pone.0196312