Angiogenesis is a physiological process involved in blood vessel formation related to several developmental process in embryo,placenta and wound healing. It is controlled by a tenuous balance between angiogenic and angiostatic factors.1 In cancer, this imbalance is getting activated in the endothelial cell leading to metastatic proliferation. When a small dormant tumor initiates angiogenesis, referred to as the angiogenic switch, it secretes factors that induce sprouting and chemotaxis of endothelial cells (ECs) towards the tumor mass.2 Within this hypoxic environment the transcription factor Hypoxia-Inducible-Factor-1- (HIF-1) is stabilized that activates the expression of multiple genes contributing to the angiogenic process. HIF-1 induced proteins include Vascular Endothelial Growth 2 Factor (VEGF) and Basic Fibroblast Growth Factor (bFGF), which promote vascular permeability and growth, respectively. Thus HIF expression is considered to be the key signalling factor in cancer angiogenesis. Other HIF- 1 induced gene products include matrix metal-loproteinases (MMPs) that breakdown the extracellular matrix to facilitate extracellular migration release growth factors such as integrins and tumor-associated macrophages (TAMs).3 Cruciferous vegetables were extensively studied for their anticancer properties due o the presence of glucosinolates and sulforaphane.4 Glu-cosinolates are hydrolyzed to isothiocyanates (ITCs). (ITCs), a bioactive component in cruciferous vegetables that are being investigated for their chemopre-ventive and chemotherapeutic effects. These chemicals are responsible for the pungent aroma and bitter flavour of cruciferous vegetables. Indole-3- carbinol (an indole) and sulforaphane (an isothiocyanate) are the most frequently examined phytochemicals for their anticancer effects.6 The purpose of the study is to design a pharmacophore model based on the interaction of different phytoconstituents of Brassicaceae with HIF and optimising it as a lead against angiogenesis. The model was developed using common pharmacophore and ligand mapping followed by its validation using Molecular Dynamics (MD) simulation.
MATERIALS AND METHODS
The study is designed for pharmacophore candidate from set of potent ligands with reported activities. Pharmacophore features were devised from set of ligands using PHARMAGIST an online web server.
It involves an indirect approach which requires a 3D training set that serves as a pivot point for alignment. The input set include the phyto-chemicals that are known to bind to HIF. The output consists of candidate pharmacophores that are generated by MFA (molecular field analy-sis).It employs the flexibility of input ligands in a deterministic manner within the alignment process.
Identification and screening of phytoconstituents and their studies
Anticancer and antiangiogenic compounds from different members of Brassicaceae were collected and reviewed. From the collection the compounds with reported anti angiogenic properties were separated and thoroughly studied using wide literature collection. From the list, constituents were screened for drug like properties using Lipinski rule of five. The resultant ligands were devised for their pharmacophore features using PHARMAGIST.
Pharmacophore model development and virtual screening
The compounds with potent antiangiogenic properties were used as base template for the generation of pharmacophore using COMMON PHARMACOPHORE GENERATION of Accelyrs Discovery Studio 2.0. (ADS). Different hypotheses were developed with diffierent parameters such as H -bond acceptor, donor and hydrophobic characteristics in common. The ligand pharmacophore and mapping protocol of ADS was used for mapping the best compounds with distinct pharmacophore features. Virtual screening of compounds was done with the developed pharmacophore in MINIMAY bridge database of ADS. The best phar-macophore models were chosen based on their structure likeness with the developed pharmacophore. All retrieved compounds were validated using drug likeness, ADMET descriptors and TOPKAT module of ADS.7,8
Binding site analysis and Molecular docking
The binding site was identified from receptor cavities using flood-filling algorithms as implemented in ADS.9 A grid resolution of 0.5 which is an indication of grid spacing and minimum number of grid points used for this analysis.10- The compounds with good drug likeness and ADMET properties were chosen and docked to the binding site of HIF using AUTODOCK of ADS 2.0. The resulting docked complexes were selected based on dock score and PMF (potential of mean force).
Molecular dynamics simulation
Molecular dynamics simulation is a computational method that calculates the time-dependent behaviour of a molecular system. Molecular dynamics simulations provide detailed information on the fluctuations and conformational changes of proteins and nucleic acids. These methods are used to investigate the structure, and thermodynamics of biological molecules and their complexes. The complex of drug like compound with target protein HIF was analyzed by producing simulation at 200 K for 1000 ns of the complex. The dynamics was then analyzed and the interacting residues were found.12–13
RESULTS AND DISCUSSION
Target protein preparation
The 3 D structure of HIF was retrieved from PDB and stabilised for docking. The structure was minimised with conjugate gradient to resolve steric clashes and unwanted interactions.
Identification and screening of phytoconstituents and their studies
From the literature search different ligands from Brassicaceae were selected.29 ligands were screened using Lipinski filter and the structures retrieved from PUBCHEM were used for the analysis. The list of ligands used for pharmacophore mapping were shown in Table 1.
Pharmacophore model generation
The compounds with good drug likeness were used as base structures for the generation of pharmacophore using PHARMAGIST. The phar-macophoric features were analysed for the 19 ligands using PHARMA-GIST. Pharmacophore model generation was done using COMMON FEATURE PHARMACOPHORE TOOL IN ADS.14 Four ligands such as Kaempferol, Glucoraphanin, Flavanol, Progoitrin were found to have best pairwise alignment with flavanol as the pivotal molecule. Ligands such as Progoitrin and Kaempferol were found to be best mapped with generated pharmacophore. Drug likeness of the above two compounds were further screened for their analogues using database. Compounds with best fit values were subjected to toxicity prediction and pharmacokinetics using ADS as shown in Figure 1(a) and (b).
HIF was found to have 10 binding sites and the docking was done using ligand fit module by Monte Carlo methods.15 The final docked poses were evaluated based on dock scores and hydrogen bonding interaction with the binding site residues as shown in Table 2 and Figure 2.15
|Site number||Molecule||PMF||Dock score|
Molecular Dynamics simulation
Molecular dynamics (MD) is a method of studying the motions in the conformational space of molecular systems. MD simulation was performed over period of 1000 ns with the temperature at 200 K for all selected top ten drug molecules with HIF. The model produced after simulation shows that the HIS115, SER 11 and TYR 98 of HIF interact with HIS 56, ALA 3 and TYR 5 of the 1st drug like compound. The 10 mol of compound showed the interaction of ALA 3 and LEU 2. The binding complex of 1st mole with HIF complex showed stable conformation throughout the simulation process. The drug like compound HTS 0115 (C15H21BrN2O3) found to have the highest dock score with the receptor as shown in Figure 3. It has the same pharmacophoric features as PROGOITRIN.
Progoitrin is a biochemical form of the glucosinolate (2 hydroxy 3 butenyl) derivative predominantly found in cabbage, broccoli, brussel sprouts and other members of Brassicaceae.16 It is reported to possess anticancer and antiangiogenic properties. It is hydrolysed by myrosinase to form cyclic thiocarbamate goitrin which is found to reduce the production of thyroid hormone. There exists an inverse relationship between the thyroid vasculature and the local availability of iodine mediated by trace element dependant regulation of VEGF (vascular endothelial growth factor). HIF is responsible for VEGF activation by perceiving decreased iodide concentration as pseudo hypoxic stimulus.17–18 Therefore, we conclude that molecule HTS 0115 (C15H21BrN2O3) has ANTI ANGIOGENIC property and further be optimised as an angiogenic inhibitor using in vitro and in vivo studies.
Natural product research recently focuses on optimising the phytocon-stituents as druggable leads in different signalling mechanism associated with cancer. Targeting intervening mechanism such as angiogenesis associated with cancer may repurpose existing leads for new therapeutic indications. In our study we devised pharmacophoe mapping from Bras-sicaceae and their molecular interaction to explore suitable HIF inhibitors. Among the drug like compounds reported from database, HTS 0115 (C15H21BrN2O3) with best dock score was validated using simulation studies. It shares the same pharmacophoric features of Progoitrin a biochemical form of glucosinolate. Hence it may be optimised in future as an effective natural lead in tumor angiogenesis.