Design, Semisynthesis, and Estrogenic Activity of Lignan Derivatives from Natural Dibenzylbutyrolactones.

Based on molecular docking studies on the ERα, a series of lignan derivatives (3-16) were designed and semisynthesized from the natural dibenzylbutyrolactones bursehernin (1) and matairesinol dimethyl ether (2). To examine their estrogenic and antiestrogenic potencies, the effects of these compounds on estrogen receptor element (ERE)-driven reporter gene expression and viability in human ER+ breast cancer cells were evaluated. Lignan compounds induced ERE-driven reporter gene expression with very low potency as compared with the pure agonist E2. However, coincubation of 5 µM of lignan derivatives 1, 3, 4, 7, 8, 9, 11, 13, and 14 with increasing concentrations of E2 (from 0.01 pM to 1 nM) reduced both the potency and efficacy of pure agonists. The binding to the rhERα-LBD was validated by TR-FRET competitive binding assay and lignans bound to the rhERα with IC50 values from 0.16 µM (compound 14) to 6 µM (compound 4). Induced fit docking (IFD) and molecular dynamics (MD) simulations for compound 14 were carried out to further investigate the binding mode interactions. Finally, the in silico ADME predictions indicated that the most potent lignan derivatives exhibited good drug-likeness.

Insight into glycogen synthase kinase-3ß inhibitory activity of phyto-constituents from Melissa officinalis: in silico studies.

Over activity of Glycogen synthase kinase-3ß (GSK-3ß), a serine/threonine-protein kinase has been implicated in a number of diseases including stroke, type II diabetes and Alzheimer disease (AD). This study aimed to find novel inhibitors of GSK-3ß from phyto-constituents of Melissa officinalis with the aid of computational analysis. Molecular docking, induced-fit docking (IFD), calculation of binding free energy via the MM-GBSA approach and Lipinski's rule of five (RO5) were employed to filter the compounds and determine their druggability. Most importantly, the compounds pIC50 were predicted by machine learning-based model generated by AutoQSAR algorithm. The generated model was validated to affirm its predictive model. The best model obtained was Model kpls_desc_38 (R2 = 0.8467 and Q2 = 0.8069), and this external validated model was utilized to predict the bioactivities of the lead compounds. While a number of characterized compounds from Melissa officinalis showed better docking score, binding free energy alongside adherence to RO5 than co-cystallized ligand, only three compounds (salvianolic acid C, ellagic acid and naringenin) showed more satisfactory pIC50. The results obtained in this study can be useful to design potent inhibitors of GSK-3ß.

Structure-Based Site of Metabolism (SOM) Prediction of Ligand for CYP3A4 Enzyme: Comparison of Glide XP and Induced Fit Docking (IFD).

Metabolism is one of the prime reasons where most of drugs fail to accomplish their clinical trials. The enzyme CYP3A4, which belongs to the superfamily of cytochrome P450 enzymes (CYP), helps in the metabolism of a large number of drugs in the body. The enzyme CYP3A4 catalyzes oxidative chemical processes and shows a very broad range of ligand specificity. The understanding of the compound's structure where oxidation would take place is crucial for the successful modification of molecules to avoid unwanted metabolism and to increase its bioavailability. For this reason, it is required to know the site of metabolism (SOM) of the compounds, where compounds undergo enzymatic oxidation. It can be identified by predicting the accessibility of the substrate's atom toward oxygenated Fe atom of heme in a CYP protein. The CYP3A4 enzyme is highly flexible and can take significantly different conformations depending on the ligand with which it is being bound. To predict the accessibility of substrate atoms to the heme iron, conventional protein-rigid docking methods failed due to the high flexibility of the CYP3A4 protein. Herein, we demonstrated and compared the ability of the Glide extra precision (XP) and Induced Fit docking (IFD) tool of Schrodinger software suite to reproduce the binding mode of co-crystallized ligands into six X-ray crystallographic structures. We extend our studies toward the prediction of SOM for compounds whose experimental SOM is reported but the ligand-enzyme complex crystal structure is not available in the Protein Data Bank (PDB). The quality and accuracy of Glide XP and IFD was determined by calculating RMSD of docked ligands over the corresponding co-crystallized bound ligand and by measuring the distance between the SOM of the ligand and Fe atom of heme. It was observed that IFD reproduces the exact binding mode of available co-crystallized structures and correctly predicted the SOM of experimentally reported compounds. Our approach using IFD with multiple conformer structures of CYP3A4 will be one of the effective methods for SOM prediction.

Identification of novel NAD(P)H dehydrogenase [quinone] 1 antagonist using computational approaches.

NAD(P)H: quinone oxidoreductase 1 (NQO1) inhibitors are proved as promising therapeutic agents against cancer. This study is to determine potent NAD(P)H-dependent NQO1 inhibitors with new scaffold. Pharmacophore-based three-dimensional (3D) QSAR model has been built based on 45 NQO1 inhibitors reported in the literature. The structure-function correlation coefficient graph represents the relationship between phase activity and phase predicted activity for training and test sets. A QSAR model statistics shows the excellent correlation of the generated model. Pharmacophore hypothesis (AARR) yielded a statistically significant 3D QSASR model with a correlation coefficient of r2 = 0.99 as well as an excellent predictive power. From the analysis of pharmacophore-based virtual screening using by SPEC database, 4093 hits were obtained and were further filtered using virtual screening filters (HTVS, SP, XP) through structure based molecular docking. Based on glide energy and docking score, seven lead compounds show better binding affinity compared to the co-crystal inhibitor. The results of induced fit docking and prime/MM-GBSA suggest that leads AN-153/J117103 and AT-138/KB09997 binding with the catalytic site. Further, to understanding the stability of identified lead compounds MD simulations were done. The lead AN-153/J117103 showed the strong binding stable of the protein-ligand complex. Also the computed drug likeness reveals potential of this compound to treat cancer. AbbreviationsNQO1NAD(P)H-quinine oxidoreductase 1CPHcommon pharmacophore hypothesisPLSpartial least squireHBDhydrogen bond donorSDstandard deviationXPextra precisionIFDinduced fit dockingMM-GBSAmolecular mechanics generalized born surface areaMDSmolecular dynamics simulationRMSDroot mean square deviationRMSFroot mean square fluctuationRMSEroot mean square errorADMEabsorption distribution metabolism excretionsCommunicated by Ramaswamy H. Sarma.

Pyridine derivatives as anticancer lead compounds with Fatty Acid Synthase as the target: An in silico-guided in vitro study.

For the past few decades, structure-based drug discovery (SBDD) has become an inevitable technique in the drug development process for screening hit compounds against therapeutic targets. Here, we have successfully used the SBDD approach viz. virtual high-throughput screening to identify potential inhibitors against the Ketoacyl synthase (KS) domain of Fatty acid synthase (FASN). Overexpression of FASN, and subsequent enhancement of de novo lipogenesis is a key survival strategy of cancer cells. Hence, targeting lipid metabolism using FASN inhibitors has been considered as a promising method to induce metabolic stress, thereby posing a survival disadvantage to cancer cells. In the present study, we have successfully identified eight FASN inhibitors from Asinex Elite database by implementing in silico tools. Five of the hit compounds share a common ring structure, which enables characteristic binding interactions with FASN-KS. Among them, in vitro validation showed that SFA 22637550 possesses significant FASN inhibitory activity and antiproliferative effect in human cancer cells of various origins. The maximum sensitivity was exhibited towards HepG2 hepatocellular carcinoma cells (IC50 = 28 µM). The mode of cell death was found to be apoptosis with a significant increase in SubG0 population without affecting any other phases of the cell cycle. The current study puts forward an excellent core structure for the development of potent FASN inhibitors for successfully targeting cancer cell metabolism, thereby causing selective cell death.

Synthesis, anticholinesterase activity and molecular modeling study of novel carbamate-substituted thymol/carvacrol derivatives.

New thymol and carvacrol derivatives with the carbamate moiety were synthesized and their inhibitory effects on acetylcholinesterase (AChE) and butyrylcholinesterase (BuChE) were evaluated. 5-isopropyl-2-methylphenyl(3-fluorophenyl)carbamate (29) was found to be the most potent AChE inhibitor with IC50 values of 2.22µM, and 5-isopropyl-2-methylphenyl (4-fluorophenyl)carbamate (30) exhibited the strongest inhibition against BuChE with IC50 value of 0.02µM. Additionally, the result of H4IIE hepatoma cell toxicity assay for compounds 18, 20, 29, 30 and 35 showed negligible cell death at 0.07-10µM. Moreover in order to better understand the inhibitory profiles of these molecules, molecular modeling studies were applied. Binding poses of studied compounds at the binding pockets of AChE and BuChE targets were determined. Predicted binding energies of these compounds as well as structural and dynamical profiles of molecules at the target sites were estimated using induced fit docking (IFD) algorithms and post-processing molecular dynamics (MD) simulations methods (i.e., Molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) approaches).

Synthesis and evaluation of naphthyl bearing 1,2,3-triazole analogs as antiplasmodial agents, cytotoxicity and docking studies.

Novel series of naphthyl bearing 1,2,3-triazoles (4a-t) were synthesized and evaluated for their in vitro antiplasmodial activity against pyrimethamine (Pyr)-sensitive and resistant strains of Plasmodium falciparum. The synthesized compounds were assessed for their cytotoxicity employing human embryonic kidney cell line (HEK-293), and none of them was found to be toxic. Among them 4j, 4k, 4l, 4m, 4n, 4t exhibited significant antiplasmodial activity in both strains, of which compounds 4m, 4n and 4t (∼3.0-fold) displayed superior activity to Pyr against resistant strain. Pyr and selected compounds (4n, 4p and 4t) that repressed parasite development also inhibited PfDHFR activity of the soluble parasite extract, suggesting that anti-parasitic activity of these compounds is a result of inhibition of the parasite DHFR. In silico studies suggest that activity of these compounds might be enhanced due to π-π stacking.

Targeting prostate cancer with compounds possessing dual activity as androgen receptor antagonists and HDAC6 inhibitors.

While enzalutamide and abiraterone are approved for treatment of metastatic castration-resistant prostate cancer (mCRPC), approximately 20-40% of patients have no response to these agents. It has been stipulated that the lack of response and the development of secondary resistance to these drugs may be due to the presence of AR splice variants. HDAC6 has a role in regulating the androgen receptor (AR) by modulating heat shock protein 90 (Hsp90) acetylation, which controls the nuclear localization and activation of the AR in androgen-dependent and independent scenarios. With dual-acting AR-HDAC6 inhibitors it should be possible to target patients who don't respond to enzalutamide. Herein, we describe the design, synthesis and biological evaluation of dual-acting compounds which target AR and are also specific towards HDAC6. Our efforts led to compound 10 which was found to have potent dual activity (HDAC6 IC50=0.0356µM and AR binding IC50=<0.03µM). Compound 10 was further evaluated for antagonist and other cell-based activities, in vitro stability and pharmacokinetics.

Design of selective TACE inhibitors using molecular docking studies: Synthesis and preliminary evaluation of anti-inflammatory and TACE inhibitory activity.

Tumor necrosis factor-α (TNF-α) converting enzyme (TACE) has been considered one of the principal therapeutic targets for the treatment of TNF-dependent pathologies. Several TACE inhibitors have been reported, but none of them has been successfully passed to phase II clinical trials. In the present work, we attempted to design highly selective new non-hydroxamate sulfonamide TACE inhibitors. The docking study was performed on one of the crystal structures of TACE, selected based on its resolution and R value, to tackle the flexibility issue of the active site. The results allowed us to distinguish the analogues with a higher binding affinity toward the active site of TACE and to identify the substituent of analogues needed for binding with the surrounding site of the enzyme. Finally the analogues were docked on crystal structures of six different matrix metalloproteinases (MMPs) for a selectivity study of TACE over MMPs. Some of these analogues were synthesized and subjected to preliminary testing for in vivo anti-inflammatory activity and TACE inhibitory activity.