Integration of fingerprint-based similarity searching and kernel-based partial least squares analysis to predict inhibitory activity against CSK, HER2, JAK1, JAK2, and JAK3.

Fingerprint-based similarity searching is an important strategy for virtual screening in drug discovery. In the present study, we carried out a systematic virtual screening study, followed by the establishment of kernel-based partial least square (KPLS) analysis prediction models for five tyrosine kinase drug targets, C-terminal SRC kinase (CSK), human epidermal growth factor 2 (HER2), and Janus kinases 1, 2, and 3 (JAK1, JAK2, and JAK3), using a dataset of 3688 compounds. These kinases are important drug discovery targets, particularly as HER2 has been validated for the treatment of metastatic breast cancer, JAK inhibitors have been validated for the clinical management of arthritis and autoimmune diseases, and CSK has been found to play an important role in bone remodeling in arthritis. We conducted similarity screenings with the most active molecule for each target in the dataset as a query using eight (8) types of two-dimensional (2D) molecular fingerprints, comprising seven Hashed fingerprints, Linear, Dendritic, Radial, Pairwise, Triplet, Torsion, and MOLSPRINT2D, and one Structural keys fingerprint, MACCS. The top ranked 1% of compounds from each target's similarity screening results was used to set up kernel-based partial least square (KPLS) prediction models, with q2 values up to 0.8. The best KPLS model for each target was selected based on its predictive ability and boot strapping results and used for prediction. This integrated study approach combining similarity screening with KPLS analysis has a high potential to enhance the accuracy and efficiency of virtual screening and thus improve the drug discovery process.

SS-4 is a highly selective small molecule inhibitor of STAT3 tyrosine phosphorylation that potently inhibits GBM tumorigenesis in vitro and in vivo.

Glioblastoma (GBM) is a highly aggressive cancer with a dismal prognosis. Constitutively active STAT3 has a causal role in GBM progression and is associated with poor patient survival. We rationally designed a novel small molecule, SS-4, by computational modeling to specifically interact with STAT3. SS-4 strongly and selectively inhibited STAT3 tyrosine (Y)-705 phosphorylation in MT330 and LN229 GBM cells and inhibited their proliferation and induced apoptosis with an IC50 of ∼100 nM. The antiproliferative and apoptotic actions of SS-4 were Y-705 phosphorylation dependent, as evidenced by its lack of effects on STAT3 knockout (STAT3KO) cells or STAT3KO cells that overexpressed a phospho-Y705 deficient (STAT3Y705F) mutant, and the recovery of effects when wild-type STAT3 or a phospho-serine (S)727 deficient mutant was expressed in STAT3KO cells. SS-4 increased the expression of STAT3 repressed genes, while decreasing the expression of STAT3 promoted genes. Importantly, SS-4 markedly reduced the growth of GBM intracranial tumor xenografts. These data together identify SS-4 as a potent STAT3 inhibitor that selectively blocks Y705-phosphorylation, induces apoptosis, and inhibits growth of human GBM models in vitro and in vivo.

Targeting MDM2 for Neuroblastoma Therapy: In Vitro and In Vivo Anticancer Activity and Mechanism of Action.

BACKGROUND: Neuroblastoma is an aggressive pediatric solid tumor with an overall survival rate of <50% for patients with high-risk disease. The majority (>98%) of pathologically-diagnosed neuroblastomas have wild-type p53 with intact functional activity. However, the mouse double minute 2 (MDM2) homolog, an E3 ubiquitin ligase, is overexpressed in neuroblastoma and leads to inhibition of p53. MDM2 also exerts p53-independent oncogenic functions. Thus, MDM2 seems to be an attractive target for the reactivation of p53 and attenuation of oncogenic activity in neuroblastoma. METHODS: In this study, we evaluated the anticancer activities and underlying mechanisms of action of SP141, a first-in-class MDM2 inhibitor, in neuroblastoma cell lines with different p53 backgrounds. The findings were confirmed in mouse xenograft models of neuroblastoma. RESULTS: We demonstrate that SP141 reduces neuroblastoma cell viability, induces apoptosis, arrests cells at the G2/M phase, and prevents cell migration, independent of p53. In addition, in neuroblastoma xenograft models, SP141 inhibited MDM2 expression and suppressed tumor growth without any host toxicity at the effective dose. CONCLUSIONS: MDM2 inhibition by SP141 results in the inhibition of neuroblastoma growth and metastasis, regardless of the p53 status of the cells and tumors. These findings provide proof-of-concept that SP141 represents a novel treatment option for both p53 wild-type and p53 null neuroblastoma.

Targeted Brain Tumor Therapy by Inhibiting the MDM2 Oncogene: In Vitro and In Vivo Antitumor Activity and Mechanism of Action.

There is a desperate need for novel and efficacious chemotherapeutic strategies for human brain cancers. There are abundant molecular alterations along the p53 and MDM2 pathways in human glioma, which play critical roles in drug resistance. The present study was designed to evaluate the in vitro and in vivo antitumor activity of a novel brain-penetrating small molecule MDM2 degrader, termed SP-141. In a panel of nine human glioblastoma and medulloblastoma cell lines, SP-141, as a single agent, potently killed the brain tumor-derived cell lines with IC50 values ranging from 35.8 to 688.8 nM. Treatment with SP-141 resulted in diminished MDM2 and increased p53 and p21cip1 levels, G2/M cell cycle arrest, and marked apoptosis. In intracranial xenograft models of U87MG glioblastoma (wt p53) and DAOY medulloblastoma (mutant p53) expressing luciferase, treatment with SP-141 caused a significant 4- to 9-fold decrease in tumor growth in the absence of discernible toxicity. Further, combination treatment with a low dose of SP-141 (IC20) and temozolomide, a standard anti-glioma drug, led to synergistic cell killing (1.3- to 31-fold) in glioma cell lines, suggesting a novel means for overcoming temozolomide resistance. Considering that SP-141 can be taken up by the brain without the need for any special delivery, our results suggest that SP-141 should be further explored for the treatment of tumors of the central nervous system, regardless of the p53 status of the tumor.

Inhibiting ß-Catenin by ß-Carboline-Type MDM2 Inhibitor for Pancreatic Cancer Therapy.

The ß-catenin and MDM2 oncoproteins are overexpressed and constitutively activated in human pancreatic cancer and contribute to its initiation, progression, and metastasis. The Wnt/ß-catenin signaling pathway strongly interacts with the MDM2-p53 signaling pathway, accelerating the tumorigenesis and its development. Therefore, therapies inhibiting both ß-catenin and MDM2 are suggested to be ideal treatments for patients with advanced pancreatic cancer. We have recently identified a novel class of ß-carboline compounds as the specific and potent MDM2 inhibitors, including a lead compound SP141. In the present study, we utilized SP141 as an exemplary ß-carboline compound to characterize ß-catenin as a molecular target of the ß-carboline compounds and to demonstrate an important role of ß-catenin in the anticancer activity of ß-carboline. We found that the silencing of either ß-catenin or MDM2 largely reduced the anticancer activity of SP141 while the double silencing of both genes almost completely blocked SP141's activity. SP141 directly bound to ß-catenin and inhibited its expression and activity in pancreatic cancer cells in vitro and in vivo. The inhibitory effects of SP141 on ß-catenin were mediated by the ubiquitin-proteasome system in an MDM2-independent manner. In conclusion, these results suggest that SP141 exerts its anticancer activity by dually inhibiting ß-catenin and MDM2. We envision that ß-carboline derivatives can be developed as promising dual inhibitors of ß-catenin and MDM2 for the treatment of advanced pancreatic cancer.

QSAR Studies of New Pyrido[3,4-b]indole Derivatives as Inhibitors of Colon and Pancreatic Cancer Cell Proliferation.

We have discovered a new class of pyrido[b]bindole derivatives that show potent and broad spectrum anticancer activity with IC50 values down to submicromolar levels. Structure-activity relationship data acquired with the compounds as antiproliferative agents against several cancer cell lines, i.e. human HCT116 colon cancer cell line, and HPAC, Mia-PaCa2 and Panc-1 pancreatic cancer cell lines, were subjected to two different QSAR modeling methods. A kernel-based partial least squares (KPLS) regression analysis with chemical 2D fingerprint descriptors, and a PHASE pharmacophore alignment with 3D-QSAR study. The KPLS method afforded successful predictive QSAR models for antiproliferative activity of the HCT116 colon cell line and on two of the pancreatic cancer cell lines HPAC and Mia-PaCa2, with the following statistics: R 2s of 0.99, 0.99 and 0.98, for training set coefficients of determination, and external test set predictive r 2s of 0.70, 0.58 and 0.70, respectively. The best 2D fingerprint descriptor for both the HCT116 and HPAC data out of the eight finger prints utilized was the atom triplet fingerprint; whereas the one that worked best for the Mia-PaCa2 data was the linear fingerprint descriptor. The PHASE pharmacophore based 3D-QSAR study afforded a four-point pharmacophore model comprising one hydrogen bond donor (D) and three ring (R) elements, which yielded a successful 3D-QSAR model only with the HCT116 cell line data with training set R 2 of 0.683, and an external test set predictive r 2 of 0.562. With the PHASE 3D-QSAR, the influence of electronic effects and hydrophobicity were visualized, and were in agreement with the observed SAR of substitutions, while the KPLS method the relative extent of contribution of each atom in a compound to the activity. These models will foster the lead optimization process for this potent series of anticancer pyrido [3,4-b]indole compounds.

Synthesis, Biological Evaluation and Modeling Studies of New Pyrido[3,4-b]indole Derivatives as Broad-Spectrum Potent Anticancer Agents.

OBJECTIVE: There is an urgent need drugs against particularly difficult to treat solid tumors such as pancreatic, triple negative breast, lung, colon, metastatic prostate cancers and melanoma. Thus, the objective of this study was to synthesize compounds based computational modeling that indicated the pyrido[3,4-b]indole class bind to MDM2, a new cancer target for which there are still no drug on the market. METHODS: Compounds were synthesized by established methods and tested for antiproliferative activity against a broad range of human cancer cell lines, comprising HCT116 colon, HPAC, MIA PaCa-2 and Panc-1 pancreatic, MCF-7 and MDA-MB-468 breast, A375 and WM164 melanoma, A549 lung, and LNCaP, DU145 and PC3 prostate cancer lines. Computational docking was also undertaken. RESULTS: The novel pyrido[3,4-b]indoles synthesized exhibited a clear SAR with regards to antiproliferative activity, with potent broad-spectrum anticancer activity with IC50s down to 80, 130, 130 and 200 nM for breast, colon, melanoma and pancreatic cancer cells, respectively. 1-Naphthyl at C1 combined with methoxy at C6 provided the best antiproliferative activity. Thus, compound 11 (1-naphthyl-6-methoxy-9H-pyrido[3,4-b]indole) showed the highest potency. A mechanistic feature of the compounds as a group is a strongly selective G2/M cell cycle phase arrest. Docking at on MDM2 suggested a hydrogen bond interaction between the 6-methoxy Tyr106, hydrophobic interaction with Val93, pi-pi stacking interactions with Tyr100 and His96 and hydrophobic interactions with Leu54 and Ile99. An N9-methyl group disrupted binding interactions, such as H-bond interactions involving the N9 hydrogen. CONCLUSION: We have identified a novel series of pyrido[3,4-b]indoles with potent broad spectrum anticancer activity towards the most aggressive and difficult to treat cancers including metastatic pancreatic cancer, non-small cell lung cancer, triple negative breast cancers, and BRAFV600E mutant melanoma, as well as metastatic colon and prostate cancers. There was also evidence of selectivity towards cancer cells relative to normal cells. These compounds will serve as new leads from which novel therapeutics and molecular tools can be developed for a wide variety of cancers.

A quantitative LC-MS/MS method for determination of SP-141, a novel pyrido[b]indole anticancer agent, and its application to a mouse PK study.

In the present study, a specific and sensitive liquid chromatography-triple quadrupole mass spectrometry method was developed and validated for the determination of SP-141, a novel pyrido[b]indole anticancer agent. After a liquid-liquid extraction with n-hexane-dichloromethane-2-propanol (20:10:1, v/v/v) mixture, the analyte was separated on a Kinetex C18 column (50×2.1mm, 2.6µm) with mobile phases comprising of water (0.1% formic acid, v/v) and acetonitrile (0.1% formic acid, v/v) at a flow rate of 0.4mL/min. The test compound (SP-141) and the internal standard (SP-157) were analyzed in the multiple reaction-monitoring mode using the mass transitions m/z 325.1 → 282.0. The method was linear in the concentration range of 0.648-162ng/mL with coefficients of determination (R(2)) of 0.999 in mouse plasma. The lower limit of quantification was 0.648ng/mL. The intra- and inter-day assay precisions (coefficient of variation, %CV) were less than 4.2% and accuracies (relative error, %RE) ranged from -6.1% to 2.1%. The extraction recoveries were between 97.1 and 103.1% and the relative matrix effect was minimal. In addition, SP-141 was found to be stable in the plasma after three freeze-thaw cycles, at 37°C and 4°C for 24h, and at -80°C for 4 weeks. It was also stable in the stock solution at room temperature for 24h and after preparation in the autosampler for 36h. The validated method was successfully applied to an initial pharmacokinetic study of SP-141 in CD-1 mice following intraperitoneal and intravenous administrations.

Synthesis and structure-activity relationship study of substituted caffeate esters as antinociceptive agents modulating the TREK-1 channel.

The TWIK-related K(+) channel, TREK-1, has recently emerged as an attractive therapeutic target for the development of a novel class of analgesic drugs. It has been reported that TREK-1 -/- mice were more sensitive than wild-type mice to painful stimuli, suggesting that activation of TREK-1 could result in pain inhibition. Here we report the synthesis of a series of substituted caffeate esters (12a-u) based on the hit compound CDC 2 (cinnamyl 3,4-dihydroxyl-α-cyanocinnamate). These analogs were evaluated for their ability to modulate TREK-1 channel by electrophysiology and for their in vivo antinociceptive activity (acetic acid induced-writhing assay) leading to the identification a series of novel molecules able to activate TREK-1 and displaying potent analgesic activity in vivo.

Identification of PDZ ligands by docking-based virtual screening for the development of novel analgesic agents.

Disrupting the interaction between the PDZ protein, PSD-95, and its target ligands (such as the glutamate NMDA receptor or the serotonin 5-HT2A receptor) was found to reduce hyperalgesia in various models of neuropathic pain. Here, we set out to identify lead molecules which would interact with PSD-95, and hence, would potentially display analgesic activity. We describe the virtual screening of the Asinex and Cambridge databases which together contain almost one million molecules. Using three successive docking filters and visual inspection, we identified three structural classes of molecules and synthesized a potential lead compound from each class. The binding of the molecules with the PDZ domains of PSD-95 was assessed by (1)H-(15)N HSQC NMR experiments. The analgesic activity of the best ligand, quinoline 2, was evaluated in vivo in a model of neuropathic pain and showed promising results.

Inhibition of gastric tumor growth by a novel Hsp90 inhibitor.

Heat shock protein 90 (Hsp90) is a molecular chaperone engaging in multiple cellular signaling by stabilizing oncoproteins (e.g. Akt and c-Raf) in tumor cells. Whereas Hsp90 inhibitors such as 17-AAG exert promising antitumor effects in clinical trials, current efforts focus on developing agents targeting Hsp90 with improved efficacy and lower toxicity. Using a fluorescence polarization assay, we screened over a hundred of synthetic small molecules and identified a resorcinol derivative LD053 that bound the Hsp90 ATP-binding pocket. The binding of LD053 to Hsp90 dissociated the co-chaperone protein cdc37 from Hsp90, resulting in destabilization of Akt and c-Raf and subsequent inhibition of PI3K/Akt and c-Raf/Mek/Erk signaling in BGC823 gastric cancer cells. As a consequence, LD053 decreased cancer cell viability and induced apoptosis evidenced by increased subG0/G1 cell population and increased cleavage of caspase 3 and PARP. Interestingly, normal human cells appeared insensitive to LD053 treatments. Consistent with its in vitro anticancer activities, LD053 significantly inhibited growth of BGC823 xenografts in nude mice without apparent body weight loss. These results thus demonstrate that LD053 is a novel Hsp90 inhibitor and has potential to be used to treat gastric cancer.

Synthesis and evaluation of bidentate ligands designed to interact with PDZ domains.

We designed bidentate ligands to target PDZ domains through two binding sites: site S0, delimited by the GLGF loop, and site S1, a zone situated around loop ß(B)/ß(C). A molecular docking study allowed us to design a generic S0 binder, to which was attached a variable size linker, itself linked to an amino acid aimed to interact with the S1 site of PDZ domains. A series of 15 novel bidentate ligands was prepared in 6-11 steps in good overall yield (24-43%). Some of these ligands showed an inhibitory activity against serotonin 5-HT2A receptor/PSD-95 interaction. This was assessed by pull-down assay using a synthetic decapeptide corresponding to the C-terminal residues of the receptor as a bait.

NMR evaluation of interactions between substituted-indole and PDZ1 domain of PSD-95.

We synthesized small organic molecules designed as PDZ ligands. These indole-based compounds were evaluated for their interaction with the PDZ1 domain of the post-synaptic density 95 (PSD-95) protein. Three molecules were found to interact with the targeted PDZ protein by NMR. One of them showed chemical shift perturbations closely related to the natural ligands.

Novel quinoline and naphthalene derivatives as potent antimycobacterial agents.

We have designed and synthesized both the quinoline and naphthalene based molecules influenced by the unique structural make-up of mefloquine and TMC207, respectively. These compounds were evaluated for their anti-mycobacterial activity against drug sensitive Mycobacterium tuberculosis H37Rv in vitro at single-dose concentration (6.25 microg/mL). The compounds 22, 23, 26 and 27 inhibited the growth of M. tuberculosis H37Rv 99%, 90%, 98% and 91% respectively. Minimum inhibitory concentration of compounds 22, 23, 26 and 27 was found to be 6.25 microg/mL. Our molecular modeling and docking studies of designed compounds showed hydrogen bonding with Glu-61, Tyr-64 and Asn-190 amino acid residues at the putative binding site of ATP synthase, these interactions were coherent as shown by Mefloquine and TMC207, where hydrogen bonding was found with Tyr-64 and Glu-61 respectively. SAR analysis indicates importance of hydroxyl group and nature of substituents on piperazinyl-phenyl ring was critical in dictating the biological activity of newly synthesized compounds.

Pharmacophore generation and atom-based 3D-QSAR of novel 2-(4-methylsulfonylphenyl)pyrimidines as COX-2 inhibitors.

Cyclooxygenase-2 (COX-2) inhibitors are widely used for the treatment of pain and inflammatory disorders such as rheumatoid arthritis and osteoarthritis. A series of novel 2-(4-methylsulfonylphenyl)pyrimidine derivatives has been reported as COX-2 inhibitors. In order to understand the structural requirement of these COX-2 inhibitors, a ligand-based pharmacophore and atom-based 3D-QSAR model have been developed. A five-point pharmacophore with four hydrogen bond acceptors (A) and one hydrogen bond donor (D) was obtained. The pharmacophore hypothesis yielded a 3D-QSAR model with good partial least-square (PLS) statistics results. The training set correlation is characterized by PLS factors (r (2) = 0.642, SD = 0.65, F = 82.7, P = 7.617 e - 12). The test set correlation is characterized by PLS factors (Q (2) (ext) = 0.841, RMSE = 0.24,Pearson-R = 0.91). A docking study revealed the binding orientations of these inhibitors at active site amino acid residues (Arg513, Val523, Phe518, Ser530, Tyr355, His90) of COX-2 enzyme. The results of ligand-based pharmacophore hypothesis and atom-based 3D-QSAR give detailed structural insights as well as highlights important binding features of novel 2-(4-methylsulfonylphenyl)pyrimidine derivatives as COX-2 inhibitors which can provide guidance for the rational design of novel potent COX-2 inhibitors.

3D-QSAR of histone deacetylase inhibitors: hydroxamate analogues.

The histone deacetylase enzyme has increasingly become an attractive target for developing novel anticancer drugs. Hydroxamates are a new class of anticancer agents reported to act by selective inhibition of the histone deacetylase (HDAC) enzyme. Comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) were employed to study three-dimensional quantitative structure-activity relationships (3D-QSARs). QSAR models were derived from a training set of 40 molecules. An external test set consisting of 17 molecules was used to validate the CoMFA and CoMSIA models. All molecules were superimposed on the template structure by atom-based, multifit and the SYBYL QSAR rigid body field fit alignments. The statistical quality of the QSAR models was assessed using the parameters r(2)(conv), r(2)(cv) and r(2)(pred). In addition to steric and electronic fields, ClogP was also taken as descriptor to account for lipophilicity. The resulting models exhibited a good conventional r(2)(conv) and cross-validated r(2)(cv) values up to 0.910 and 0.502 for CoMFA and 0.987 and 0.534 for CoMSIA. Robust cross-validation by 2 groups was performed 25 times to eliminate chance correlation. The CoMFA models exhibited good external predictivity as compared to that of CoMSIA models. These 3D-QSAR models are very useful for design of novel HDAC inhibitors.

3D-QSAR of histone deacetylase inhibitors as anticancer agents by genetic function approximation.

Histone deacetylases (HDACs) play a critical role in gene transcription and are implicated in cancer therapy and other diseases. Inhibitors of HDACs induce cell differentiation and suppress cell proliferation in the tumor cells. Although many such inhibitors have been designed and synthesized, but selective inhibitors for HDAC isoforms are lacking. Various hydroxamic acid analogues have been reported as HDAC inhibitors. Here, we report a three-dimensional quantitative structure-activity relationship (3D-QSAR) study performed using genetic function approximation (GFA) for this class of molecules. QSAR models were generated using a training set of 39 molecules and the predictive ability of final model was assessed using a test set of 17 molecules. The internal consistency of the final QSAR model was 0.712 and showed good external predictivity of 0.585. The results of the present QSAR study indicated that molecular shape analysis (MSA). thermodynamic and structural descriptors are important for inhibition of HDACs.