Discovery of novel inhibitors of vascular endothelial growth factor-A–Neuropilin-1 interaction by structure-based virtual screening

Neuropilin-1 (NRP-1), one of the most important co-receptors of vascular endothelial growth factor-A (VEGF-A), increases its angiogenic action in several chronic diseases including cancer by increasing the activity of associated tyrosine kinase receptors, VEGFR1 and VEGFR2. Binding of VEGF-A to NRP-1 plays a critical role in pathological angiogenesis and tumor progression. Today, targeting this interaction is a validated approach to fight against angiogenesis-dependent diseases. Only anti-NRP-1 antibodies, peptide and peptidomimetic drug-candidates or hits have been developed thus far. In order to identify potent orally active small organic molecules various experimental and in silico approaches can be used. Here we report, novel promising small drug-like molecules disrupting the binding of VEGF-A₁₆₅ to NRP-1. We carried out structure-based virtual screening experiments using the ChemBridge compound collection on the VEGF-A₁₆₅ binding pocket of NRP-1. After docking and two rounds of similarity search computations, we identified 4 compounds that inhibit the biotinylated VEGF-A₁₆₅ binding to recombinant NRP-1 with K_i of about 10 μM. These compounds contain a common chlorobenzyloxy alkyloxy halogenobenzyl amine scaffold that can serve as a base for further development of new NRP-1 inhibitors.     

VEGF-A165 induces human aortic smooth muscle cell migration by activating neuropilin-1-VEGFR1-PI3K axis

Vascular smooth muscle cells (SMCs), one of the major cell types of the vascular wall, play a critical role in the process of angiogenesis under both physiological and pathophysiological conditions, including the cancer microenvironment. Previous studies have shown that VEGF-A 165 augments vascular SMC migration via VEGFR2 (KDR/Flk1) pathways. In this study, we found that VEGF-A 165 (recombinant protein or breast tumor cell-secreted) is also capable of inducing migration of VEGFR2-negative human aortic smooth muscle cells (hAOSMCs), and this induction is mediated through a molecular cross-talk of neuropilin-1 (NRP-1), VEGFR1 (Flt-1), and phosphoinositide 3-kinase (PI3K)/Akt signaling kinase. We found that VEGF-A 165 induces hAOSMC migration parallel with the induction of NRP-1 and VEGFR1 expressions and their associations along with the activation of PI3K/Akt. Neutralization of VEGF action by its antibody or inhibition of VEGF-induced PI3K/Akt kinase activation by wortmannin, a PI3K/Akt specific inhibitor, results in inhibition of VEGF-induced hAOSMC migration. Moreover, RNAi-mediated elimination of the NRP-1 expression or blocking of the activity of VEGFR1 by its antibody in hAOSMCs impairs the VEGF-A 165-induced migration of these cells as well as activation of PI3K/Akt kinase. Collectively, these results establish, for the first time, a mechanistic link among VEGF-A 165, NRP-1, VEGFR1, and PI3K/Akt in the regulation of migration of human vascular smooth muscle cells that eventually could be involved in the angiogenic switch.     

Exogenous recombinant dimeric neuropilin-1 is sufficient to drive angiogenesis

Neuropilin-1 (NRP-1) is present on the cell surface of endothelial cells, or as a soluble truncated variant. Membrane NRP-1 is proposed to enhance angiogenesis by promoting the formation of a signaling complex between vascular endothelial growth factor-A(165) (VEGF-A(165)), VEGF receptor-2 (VEGFR-2) and heparan sulfate, whereas the soluble NRP-1 is thought to act as an antagonist of signaling complex formation. We have analyzed the angiogenic potential of a chimera comprising the entire extracellular NRP-1 region dimerized through an Fc IgG domain and a monomeric truncated NRP-1 variant. Both NRP-1 proteins stimulated tubular morphogenesis and cell migration in HDMECs and HUVECs. Fc rNRP-1 was able to induce VEGFR-2 phosphorylation and expression of the VEGFR-2 specific target, regulator of calcineurin-1 (RCAN1.4). siRNA mediated gene silencing of VEGFR-2 revealed that VEGFR-2 was required for Fc rNRP-1 mediated activation of the intracellular signaling proteins PLC-γ, AKT, and MAPK and tubular morphogenesis. The stimulatory activity was independent of VEGF-A(165). This was evidenced by depleting the cell culture of exogenous VEGF-A(165), and using instead for routine culture VEGF-A(121), which does not interact with NRP-1, and by the inability of VEGF-A sequestering antibodies to inhibit the angiogenic activity of the NRP proteins. Analysis of angiogenesis over a period of 6 days in an in vitro fibroblast/endothelial co-culture model revealed that Fc rNRP-1 could induce endothelial cell tubular morphogenesis. Thus, we conclude that soluble Fc rNRP-1 is a VEGF-A(165)-independent agonist of VEGFR-2 and stimulates angiogenesis in endothelial cells.     

Pharmacophore modeling,multiple docking,and molecular dynamics studies on Wee1 kinase inhibitors

Wee1-like protein kinase (Wee1) is a tyrosine kinase that regulates the G₂ checkpoint and prevents entry into mitosis in response to DNA damage. Based on a series of signaling pathways initiated by Wee1, Wee1 has been recognized as a potential target for cancer therapy. To discover potent Wee1 inhibitors with novel scaffolds, ligand-based pharmacophore model has been built based on 101 known Wee1 inhibitors. Then the best pharmacophore model, AADRRR.340, with good partial least square (PLS) statistics (R²?=?0.9212, Q²?=?0.7457), was selected and validated. The validated model was used as a three-dimensional (3D) search query for databases virtual screening. The filtered molecules were further analyzed and refined by Lipinski’s rule of 5, multiple docking procedures (high throughput virtual screening (HTVS), standard precision (SP), genetic optimization for ligand docking (GOLD), extra precision (XP), and unique quantum polarized ligand docking (QPLD)); absorption, distribution, metabolism, excretion, and toxicity (ADMET) screening; and the Prime/molecular mechanics generalized born surface area (MM-GBSA) method binding free energy calculations. Eight leads were identified as potential Wee1 inhibitors, and a 50?ns molecular dynamics (MD) simulation was carried out for top four inhibitors to predict the stability of ligand–protein complex. Molecular mechanics Poisson–Boltzmann surface area (MM-PBSA) based on MD simulation and the energy contribution per residue to the binding energy were calculated. In the end, three hits with good stabilization and affinity to protein were identified.

Communicated by Ramaswamy H. Sarma     

Carbohydrate-based peptidomimetics targeting neuropilin-1: Synthesis,molecular docking study and in vitro biological activities

Neuropilin-1 (NRP-1), a transmembrane glycoprotein acting as a co-receptor of VEGF-A, is expressed by cancer and angiogenic endothelial cells and is involved in the angiogenesis process. Taking advantage of functionalities and stereodiversities of sugar derivatives, the design and the synthesis of carbohydrate based peptidomimetics are here described. One of these compounds (56) demonstrated inhibition of VEGF-A₁₆₅ binding to NRP-1 (IC₅₀ = 39 μM) and specificity for NRP-1 over VEGF-R2. Biological evaluations were performed on human umbilical vein endothelial cells (HUVECs) through activation of downstream proteins (AKT and ERK phosphorylation), viability/proliferation assays and in vitro measurements of anti-angiogenic abilities.     

Identification of Pak1 inhibitors using water thermodynamic analysis

Abstract

p21-activated kinases (Paks) play an integral component in various cellular diverse processes. The full activation of Pak is dependent upon several serine residues present in the N-terminal region, a threonine present at the activation loop, and finally the phosphorylation of these residues ensure the complete activation of Pak1. The present study deals with the identification of novel potent candidates of Pak1 using computational methods as anti-cancer compounds. A diverse energy based pharmacophore (e-pharmacophore) was developed using four co-crystal inhibitors of Pak1 having pharmacophore features of 5 (DRDRR), 6 (DRHADR), and 7 (RRARDRP and DRRDADH) hypotheses. These models were used for rigorous screening against e-molecule database. The obtained hits were filtered using ADME/T and molecular docking to identify the high affinity binders. These hits were subjected to hierarchical clustering using dendritic fingerprint inorder to identify structurally diverse molecules. The diverse hits were scored against generated water maps to obtain WM/MM ΔG binding energy. Furthermore, molecular dynamics simulation and density functional theory calculations were performed on the final hits to understand the stability of the complexes. Five structurally diverse novel Pak1 inhibitors (4835785, 32198676, 32407813, 76038049, and 32945545) were obtained from virtual screening, water thermodynamics and WM/MM ΔG binding energy. All hits revealed similar mode of binding pattern with the hinge region residues replacing the unstable water molecules in the binding site. The obtained novel hits could be used as a platform to design potent drugs that could be experimentally tested against cancer patients having increased Pak1 expression.     

Discovery of promising FtsZ inhibitors by E-pharmacophore, 3D-QSAR,molecular docking study,and molecular dynamics simulation

Filamentous temperature-sensitive protein Z (FtsZ), playing a key role in bacterial cell division, is regarded as a promising target for the design of antimicrobial agent. This study is looking for potential high-efficiency FtsZ inhibitors. Ligand-based pharmacophore and E-pharmacophore, virtual screening and molecular docking were used to detect promising FtsZ inhibitors, and molecular dynamics simulation was used to study the stability of protein-ligand complexes in this paper. Sixty-three inhibitors from published literatures with pIC₅₀ ranging from 2.483 to 5.678 were collected to develop ligand-based pharmacophore model. 4DXD bound with 9PC was selected to develop the E-pharmacophore model. The pharmacophore models validated by test set method and decoy set were employed for virtual screening to exclude inactive compounds against ZINC database. After molecular docking, ADME analysis, IFD docking and MM-GBSA, 8 hits were identified as potent FtsZ inhibitors. A 50?ns molecular dynamics simulation was implemented on the compounds to assess the stability between potent inhibitors and FtsZ. The results indicated that the candidate compounds had a high docking score and were strongly combined with FtsZ by forming hydrogen bonding interactions with key amino acid residues, and van der Waals forces and hydrophobic interactions had significant contribution to the stability of the binding. Molecular dynamics simulation results showed that the protein-ligand compounds performed well in both the stability and flexibility of the simulation process.     

Inhibitor design against JNK1 through e-pharmacophore modeling docking and molecular dynamics simulations

c-Jun-NH2 terminal kinases (JNKs) come under a class of serine/threonine protein kinases and are encoded by three genes, namely JNK1, JNK2 and JNK3. Human JNK1 is a cytosolic kinase belonging to mitogen-activated protein kinase (MAPK) family, which plays a major role in intracrinal signal transduction cascade mechanism. Overexpressed human JNK1, a key kinase interacts with other kinases involved in the etiology of many cancers, such as skin cancer, liver cancer, breast cancer, brain tumors, leukemia, multiple myeloma and lymphoma. Thus, to unveil a novel human JNK1 antagonist, receptor-based pharmacophore modeling was performed with the available eighteen cocrystal structures of JNK1 in the protein data bank. Eighteen e-pharmacophores were generated from the 18 cocrystal structures. Four common e-pharmacophores were developed from the 18 e-pharmacophores, which were used as three-dimensional (3D) query for shape-based similarity screening against more than one million small molecules to generate a JNK1 ligand library. Rigid receptor docking (RRD) performed using GLIDE v6.3 for the 1683 compounds from in-house library and 18 cocrystal ligands with human JNK1 from lower stringency to higher stringency revealed 17 leads. Further to derive the best leads, dock complexes obtained from RRD were studied further with quantum-polarized ligand docking (QPLD), induced fit docking (IFD) and molecular mechanics/generalized Born surface area (MM-GBSA). Four leads have showed lesser binding free energy and better binding affinity towards JNK1 compared to 18 cocrystal ligands. Additionally, JNK1–lead1 complex interaction stability was reasserted using 50?ns MD simulations run and also compared with the best resolute cocrystal structure using Desmond v3.8. Thus, the results obtained from RRD, QPLD, IFD and MD simulations indicated that lead1 might be used as a potent antagonist toward human JNK1 in cancer therapeutics.     

Structure-activity relationship study of tetrapeptide inhibitors of the Vascular Endothelial Growth Factor A binding to Neuropilin-1

Neuropilin-1 is considered as one of the key receptors responsible for signaling pathways involved in pathological angiogenesis necessary for tumor progression, therefore targeting of VEGF₁₆₅ binding to NRP-1 could be a relevant strategy for antiangiogenic treatment. It was shown before that the VEGF₁₆₅/NRP-1 interaction can be inhibited by short tetrapeptides with K/RXXR sequence.Here, we present a structure–activity relationship study of the systematic optimization of amino acid residues in positions 1–3 in the above tetrapeptides. All the 13 synthesized analogs possessed C-terminal arginine that is a necessary element for interaction with NRP-1. The obtained results of the inhibitory activity and modeling by molecular dynamics indicate that simultaneous interactions of the basic amino acid residues in position 1 and 4 (Arg) with Neuropilin-1 are crucial and their cooperation strongly affects the inhibitory activity. In addition, the binding strength is modulated by the flexibility of the peptide backbone (in the central part of the peptide), and the nature of the side chain of the amino acids at the second or third position. A dramatic decrease in the activity to the receptor is observed in flexible derivatives that are missing proline residues. The results described in this paper should prove useful for future studies aimed at establishing the best pharmacophore for inhibitors of VEGF₁₆₅ binding to NRP-1.     

Exploration of potential RSK2 inhibitors by pharmacophore modelling,structure-based 3D-QSAR,molecular docking study and molecular dynamics simulation

Abstract

The p90 ribosomal s6 kinase 2 (RSK2) is a promising target because of its over expression and activation in human cancer cells and tissues. Over the last few years, significant efforts have been made in order to develop RSK2 inhibitors to treat myeloma, prostatic cancer, skin cancer and etc., but with limited success so far. In this paper, pharmacophore modelling, molecular docking study and molecular dynamics (MD) simulation have been performed to explore the novel inhibitors of RSK2. Pharmacophore models were developed by 95 molecules having pIC₅₀ ranging from 4.577 to 9.000. The pharmacophore model includes one hydrogen bond acceptor (A), one hydrogen bond donor (D), one hydrophobic feature (H) and one aromatic ring (R). It is the best pharmacophore hypothesis that has the highest correlation coefficient (R² = 0.91) and cross validation coefficient (Q² = 0.71) at 5 component PLS factor. It was evaluated using enrichment analysis and the best model was used for virtual screening. The constraints used in this study were docking score, ADME properties, binding free energy estimates and IFD Score to screen the database. Ultimately, 12 hits were identified as potent and novel RSK2 inhibitors. A 15 ns molecular dynamics (MD) simulation was further employed to validate the reliability of the docking results.     

Antiangiogenic and antitumor activities of peptide inhibiting the vascular endothelial growth factor binding to neuropilin-1

Neuropilin-1 (NRP-1), a non-tyrosine kinase receptor of vascular endothelial growth factor-165 (VEGF165), was found expressed on endothelial and some tumor cells. Since its overexpression is correlated with tumor angiogenesis and progression, the targeting of NRP-1 could be a potential anti-cancer strategy. To explore this hypothesis, we identified a peptide inhibiting the VEGF165 binding to NRP-1 and we tested whether it was able to inhibit tumor growth and angiogenesis. To prove the target of peptide action, we assessed its effects on binding of radiolabeled VEGF165 to recombinant receptors and to cultured cells expressing only VEGFR-2 (KDR) or NRP-1. Antiangiogenic activity of the peptide was tested in vitro in tubulogenesis assays and in vivo in nude mice xenotransplanted in fat-pad with breast cancer MDA-MB-231 cells. Tumor volumes, vascularity and proliferation indices were determined. The selected peptide, ATWLPPR, inhibited the VEGF165 binding to NRP-1 but not to tyrosine kinase receptors, VEGFR-1 (flt-1) and KDR; nor did it bind to heparin. It diminished the VEGF-induced human umbilical vein endothelial cell proliferation and tubular formation on Matrigel and in co-culture with fibroblasts. Administration of ATWLPPR to nude mice inhibited the growth of MDA-MB-231 xenografts, and reduced blood vessel density and endothelial cell area but did not alter the proliferation indices of the tumor. In conclusion, ATWLPPR, a previously identified KDR-interacting peptide, was shown to inhibit the VEGF165 interactions with NRP-1 but not with KDR and to decrease the tumor angiogenesis and growth, thus validating, in vivo, NRP-1 as a possible target for antiangiogenic and antitumor agents.     

Structural insights into the binding mode of flavonols with the active site of matrix metalloproteinase-9 through molecular docking and molecular dynamic simulations studies

Cartilage degradation in rheumatoid arthritis is mediated principally by the collagenases and gelatinases. Gelatinase B (also called matrix metalloproteinase 9 – MMP-9), is a valid target molecule which is known to participate in cartilage degradation as well as angiogenesis associated with the disease and inhibition of its activity shall prevent cartilage damage and angiogenesis. The focus of this study is to investigate the possibilities of MMP-9 inhibition by flavonol class of bioflavonoids by studying their crucial binding interactions at the active site of MMP 9 using molecular docking (Glide XP and QPLD) and further improvisation by post-docking MM-GBSA and molecular dynamic (MD) simulations. The results show that flavonols can convincingly bind to active site of MMP-9 as demonstrated by their stable interactions at the S1′ specificity pocket and favourable binding energies. Gossypin has emerged as a promising candidate with a docking score of ?14.618 kcal/mol, binding energy of ?79.97 kcal/mol and a stable MD pattern over 15 ns. In addition, interaction mechanisms with respect to catalytic site zinc are also discussed. Further, the drug-like characters of the ligands were also analysed using ADME analysis.     

Vascular endothelial growth factor (VEGF)-A165-induced prostacyclin synthesis requires the activation of VEGF receptor-1 and -2 heterodimer

We previously reported that vascular endothelial growth factor (VEGF)-A(165) inflammatory effect is mediated by acute platelet-activating factor synthesis from endothelial cells upon the activation of VEGF receptor-2 (VEGFR-2) and its coreceptor, neuropilin-1 (NRP-1). In addition, VEGF-A(165) promotes the release of other endothelial mediators including nitric oxide and prostacyclin (PGI(2)). However, it is unknown whether VEGF-A(165) is mediating PGI(2) synthesis through VEGF receptor-1 (VEGFR-1) and/or VEGF receptor-2 (VEGFR-2) activation and whether the coreceptor NRP-1 potentiates VEGF-A(165) activity. In this study, PGI(2) synthesis in bovine aortic endothelial cells (BAEC) was assessed by quantifying its stable metabolite (6-keto prostaglandin F(1alpha), 6-keto PGF(1alpha)) by enzyme-linked immunosorbent assay. Treatment of BAEC with VEGF analogs, VEGF-A(165) (VEGFR-1, VEGFR-2 and NRP-1 agonist) and VEGF-A(121) (VEGFR-1 and VEGFR-2 agonist) (up to 10(-9) m), increased PGI(2) synthesis by 70- and 40-fold within 15 min. Treatment with VEGFR-1 (placental growth factor and VEGF-B) or VEGFR-2 (VEGF-C) agonist did not increase PGI(2) synthesis. The combination of VEGFR-1 and VEGFR-2 agonists did not increase PGI(2) release. Pretreatment with a VEGFR-2 inhibitor abrogated PGI(2) release mediated by VEGF-A(165) and VEGF-A(121), and pretreatment of BAEC with antisense oligomers targeting VEGFR-1 or VEGFR-2 mRNA reduced PGI(2) synthesis mediated by VEGF-A(165) and VEGF-A(121) up to 79%. In summary, our data demonstrate that the activation of VEGFR-1 and VEGFR-2 heterodimer (VEGFR-1/R-2) is essential for PGI(2) synthesis mediated by VEGF-A(165) and VEGF-A(121), which cannot be reproduced by the parallel activation of VEGFR-1 and VEGFR-2 homodimers with corresponding agonists. In addition, the binding of VEGF-A(165) to NRP-1 potentiates its capacity to promote PGI(2) synthesis.     

Pharmacoinformatics analysis to identify inhibitors of Mtb-ASADH

Aspartate-semialdehyde dehydrogenase (ASADH; EC 1.2.1.11) is a key enzyme in the biosynthesis of essential amino acids in prokaryotes and fungi, inhibition of ASADH leads to the development of novel antitubercular agents. In the present work, a combined structure and ligand-based pharmacophore modeling, molecular docking, and molecular dynamics (MD) approaches were employed to identify potent inhibitors of mycobacterium tuberculosis (Mtb)-ASADH. The structure-based pharmacophore hypothesis consists of three hydrogen bond acceptor (HBA), two negatively ionizable, and one positively ionizable center, while ligand-based pharmacophore consists of additional one HBA and one hydrogen bond donor features. The validated pharmacophore models were used to screen the chemical databases (ZINC and NCI). The screened hits were subjected to ADME and toxicity filters, and subsequently to the molecular docking analysis. Best-docked 25 compounds carry the characteristics of highly electronegative functional groups (–COOH and –NO₂) on both sides and exhibited the H-bonding interactions with highly conserved residues Arg99, Arg249, and His256. For further validation of docking results, MD simulation studies were carried out on two representative compounds NSC51108 and ZINC04203124. Both the compounds remain bound to the key active residues of Mtb-ASADH during the MD simulations. These identified hits can be further used for lead optimization and in the design more potent inhibitors against Mtb-ASADH.     

Stepwise development of structure–activity relationship of diverse PARP-1 inhibitors through comparative and validated in silico modeling techniques and molecular dynamics simulation

Inhibitors of poly (ADP-ribose) polymerase-1 (PARP-1) enzyme are useful for the treatment of various diseases including cancer. Comparative in silico studies were performed on different ligand-based (2D-QSAR, Kernel-based partial least square (KPLS) analysis, Pharmacophore Search Engine (PHASE) pharmacophore mapping), and structure-based (molecular docking, MM-GBSA analyses, Gaussian-based 3D-QSAR analyses on docked poses) modeling techniques to explore the structure–activity relationship of a diverse set of PARP-1 inhibitors. Two-dimensional (2D)-QSAR highlighted the importance of charge topological index (JGI7), fractional polar surface area (Jurs_FPSA3), and connectivity index (CIC2) along with different molecular fragments. Favorable and unfavorable fingerprints were demonstrated in KPLS analysis, whereas important pharmacophore features (one acceptor, one donor, and two ring aromatic) along with favorable and unfavorable field effects were demonstrated in PHASE-based pharmacophore model. MM-GBSA analyses revealed significance of different polar, non-polar, and solvation energies. Docking-based alignment of ligands was used to perform Gaussian-based 3D-QSAR study that further demonstrated importance of different field effects. Overall, it was found that polar interactions (hydrogen bonding, bridged hydrogen bonding, and pi–cation) play major roles for higher activity. Steric groups increase the total contact surface area but it should have higher fractional polar surface area to adjust solvation energy. Structure-based pharmacophore mapping spotted the positive ionizable feature of ligands as the most important feature for discriminating highly active compounds from inactives. Molecular dynamics simulation, conducted on highly active ligands, described the dynamic behaviors of the protein complexes and supported the interpretations obtained from other modeling analyses. The current study may be useful for designing PARP-1 inhibitors.     

Combined pharmacophore modeling, 3D-QSAR and docking studies to identify novel HDAC inhibitors using drug repurposing

Abstract

Histone deacetylases (HDACs), a critical family of epigenetic enzymes, has emerged as a promising target for antitumor drugs. Here, we describe our protocol of virtual screening in identification of novel potential HDAC inhibitors through pharmacophore modeling, 3D-QSAR, molecular docking and molecular dynamics (MD) simulation. Considering the limitation of current virtual screening works, drug repurposing strategy was applied to discover druggable HDAC inhibitor. The ligand-based pharmacophore and 3D-QSAR models were established, and their reliability was validated by different methods. Then, the DrugBank database was screened, followed by molecular docking. MD simulation (100?ns) was performed to further study the stability of ligand binding modes. Finally, results indicated the hit DB03889 with high in silico inhibitory potency was suitable for further experimental analysis.

Communicated by Ramaswamy H. Sarma     

Discovery of novel 5α-reductase type II inhibitors by pharmacophore modelling,virtual screening,molecular docking and molecular dynamics simulations

Benign prostatic hyperplasia (BPH) is caused by augmented levels of androgen dihydrotestosterone (DHT) which is involved in the growth of the prostate in humans. 5α-Reductase type II (5αR2) is an intracellular enzyme that catalyses the formation of DHT from testosterone; hence, the inhibition of 5αR2 has emerged as one of the most promising strategies for the treatment of BPH. In this study, a computational approach that integrates ligand-based pharmacophore modelling, virtual screening, molecular docking and molecular dynamics (MD) simulations was adopted to discover novel 5αR2 inhibitors with less side effects. After validating by Fischer's randomisation and Güner–Henry test, the best quantitative pharmacophore model (Hypo1), consisting of two hydrogen-bond acceptors and three hydrophobic features, was subsequently used as a three-dimensional-query in virtual screening to identify potential hits from Maybridge and National Cancer Institute databases. These hits were further filtered by ADMET (absorption, distribution, metabolism, elimination and toxicology) and molecular docking experiments, and their binding stabilities were validated by 10-ns MD simulations. Finally, only one hit was identified as a potential lead based on higher predicted inhibitory activity to 5αR2 compared with the most active inhibitor (finasteride). Our results further suggest that this potential lead could easily be synthesised and has structural novelty, making it a promising candidate for treating BPH.     

Overexpression of VEGF189 in breast cancer cells induces apoptosis via NRP1 under stress conditions

The existence of multiple VEGF-A isoforms raised the possibility that they may have distinct functions in tumor growth. We have previously published that VEGF189 and VEGF165 contribute to breast cancer progression and angiogenesis, but VEGF165 induced the most rapid tumor uptake. Since VEGF165 has been described as a survival factor for breast tumor cells, we questioned here the effects of VEGF189 on the survival/apoptosis of MDA-MB-231 cells. We used clones that overexpress VEGF189 (V189) or VEGF165 (V165) isoforms and compared them to a control one (cV). Overexpression of VEGF189 resulted in increased cell apoptosis, as determined by Annexin-V apoptosis assay, under serum starvation and doxorubicin treatment, while VEGF 165 was confirmed to be a survival factor. Since MDA-MB-231 highly express NRP1 (a co-receptor for VEGF-A), we used short hairpin RNA (shRNA) to knock down NRP1 expression. V189shNRP1 clones were characterized by reduced apoptosis and higher necrosis, as compared with V189shCtl, under stress conditions. Unexpectedly, NRP1 knockdown had no effect on the survival or apoptosis of V165 cells. VEGF189 showed greater affinity toward NRP1 than VEGF165 using a BIAcore binding assay. Finally, since endogenously produced urokinase-type plasminogen (uPA) has been found to prevent apoptosis in breast cancers, we analyzed the level of uPA activity in our clones. An inhibition of uPA activity was observed in V189shNRP1 clones. Altogether, these results suggest a major role of NRP1 in apoptosis induced by VEGF189 in stress conditions and confirm VEGF165 as a survival factor.Key words: VEGF isoforms, survival, apoptosis, NRP-1, breast cancer cells     

Structure-function analysis of the antiangiogenic ATWLPPR peptide inhibiting VEGF(165) binding to neuropilin-1 and molecular dynamics simulations of the ATWLPPR/neuropilin-1 complex

Heptapeptide ATWLPPR (A7R), identified in our laboratory by screening a mutated phage library, was shown to bind specifically to neuropilin-1 (NRP-1) and then to selectively inhibit VEGF₁₆₅ binding to this receptor. In vivo, treatment with A7R resulted in decreasing breast cancer angiogenesis and growth. The present work is focused on structural characterization of A7R. Analogs of the peptide, obtained by substitution of each amino acid with alanine (alanine-scanning) or by amino acid deletion, have been systematically assayed to determine the relative importance of the side chains of each residue with respect to the inhibitory effect of A7R on VEGF₁₆₅ binding to NRP-1. We show here the importance of the C-terminal sequence LPPR and particularly the key role of C-terminal arginine. In solution, A7R displays significant secondary structure of the backbone adopting an extended conformation. However, the functional groups of arginine are very flexible in the absence of NRP-1 pointing to an induced fit upon binding to the receptor. A MD trajectory of the A7R/NRP-1 complex in explicit water, based on the recent tuftsin/NRP-1 crystal structure, has revealed the hydrogen-bonding network that contributes to A7R's binding activity.     

Ligand- and structure-based in silico studies to identify kinesin spindle protein (KSP) inhibitors as potential anticancer agents

Kinesin spindle protein (KSP) belongs to the kinesin superfamily of microtubule-based motor proteins. KSP is responsible for the establishment of the bipolar mitotic spindle which mediates cell division. Inhibition of KSP expedites the blockade of the normal cell cycle during mitosis through the generation of monoastral MT arrays that finally cause apoptotic cell death. As KSP is highly expressed in proliferating/cancer cells, it has gained considerable attention as a potential drug target for cancer chemotherapy. Therefore, this study envisaged to design novel KSP inhibitors by employing computational techniques/tools such as pharmacophore modelling, virtual database screening, molecular docking and molecular dynamics. Initially, the pharmacophore models were generated from the data-set of highly potent KSP inhibitors and the pharmacophore models were validated against in house test set ligands. The validated pharmacophore model was then taken for database screening (Maybridge and ChemBridge) to yield hits, which were further filtered for their drug-likeliness. The potential hits retrieved from virtual database screening were docked using CDOCKER to identify the ligand binding landscape. The top-ranked hits obtained from molecular docking were progressed to molecular dynamics (AMBER) simulations to deduce the ligand binding affinity. This study identified MB-41570 and CB-10358 as potential hits and evaluated these experimentally using in vitro KSP ATPase inhibition assays.