Energy decomposition and waterswapping analysis to investigate the SNP associated DPD mediated 5-FU resistance.

5-fluorouracil is an essential component of systemic chemotherapy for colon, breast, head, and neck cancer patients. However, tumoral overexpression of the dihydropyrimidine dehydrogenase has rendered 5-FU clinically ineffective by inactivating it to 5'-6'-dihydro fluorouracil. The responses to 5-FU in terms of efficacy and toxicity greatly differ depending upon the population group, because of variability in the DPD activity levels. In the current study, key active site amino acids involved in the 5-FU inactivation were investigated by modelling the 3D structure of human DPD in a complex with 5-FU. The identified amino acids were analyzed for their possible missense mutations available in dbSNP database. Out of 12 missense SNPs, four were validated either by sequencing in the 1000 Genomes project or frequency/genotype data. The recorded validated missense SNPs were further considered to analyze the effect of their respective alterations on 5-FU binding. Overall findings suggested that population bearing the Glu611Val DPD mutation (rs762523739) is highly vulnerable to 5-FU resistance. From the docking, electrostatic complementarity, dynamics, and energy decomposition analyses it was found that the above mutation showed superior scores than the wild DPD -5FU complex. Therefore, prescribing prodrug NUC-3373 or DPD inhibitors (Gimeracil/3-Cyano-2,6-Dihydroxypyridines) as adjuvant therapy may overcome the 5-FU resistance.

In silico guided design of non-covalent inhibitors of DprE1: synthesis and biological evaluation.

DprE1 is a potential target of resistant tuberculosis (TB), especially multidrug-resistant (MDR) and extensively drug-resistant (XDR) TB. 2-benzoxazolinone is a closely related bioisostere of some scaffolds such as benzoxazoles, benzimidazole, benzothiazolinone, and benzothiazoles that have been previously explored against DprE1. Thus, a ligand-based quantitative pharmacophore model (AHRR.8) of DprE1 was developed and this pharmacophore model was utilized in activity profiling of some 2-benzoxazolinones from an in-house database using virtual screening. Obtained hits were subject to molecular docking, molecular dynamics (MD), and MM/GBSA calculations, which resulted in benzoyl-substituted derivatives of 2-benzoxazolinone showing strong interactions with the key amino acid residues in the active site of DprE1. Based on in silico results, the top five hits were duly synthesized and evaluated against the XDR-TB strain. This study is an initial effort to explore 2-benzoxazolinones against XDR-TB, which can be submitted further to lead optimization for refining the results.

In silico guided designing of 4-(1H-benzo[d]imidazol-2-yl)phenol-based mutual-prodrugs of NSAIDs: synthesis and biological evaluation.

The free COOH group of conventional NSAIDs is a structural feature for non-selective cyclooxygenase (COX) inhibition and the molecular cause of their gastrointestinal (GI) toxicity. In this context, an in house database of synthesizable ester prodrugs of some well-known NSAIDs was developed by combining their -COOH group with -OH of a newly identified antioxidant 4-(1H-benzo[d]imidazol-2-yl)phenol (BZ). The antioxidant potential of BZ was unveiled through in silico PASS prediction and in vitro/in vivo evaluation. The in house database of NSAIDs-BZ prodrugs was first subjected to screening with our previously reported pharmacophore models of hCES1 (AAHRR.430) and hCES2 (AHHR.21) for determining hydrolytic susceptibility. Biotransformation behaviour of screened prodrugs was then assessed by using QM/MM and sterimol parameterization, followed by ADMET calculations to predict the drug likeness. On the basis of in silico results, five prodrugs were duly synthesized and the best three were subject to the in vivo evaluation for their anti-inflammatory, analgesic, antioxidant activities, and ulcerogenic index. Among these prodrugs, BN2 and BN5 displayed better anti-inflammatory and analgesics potential in comparison to their parent drugs. All the prodrugs were found to be gastro sparing in the rat model and significantly improved the levels of oxidative stress biomarkers in both blood plasma as well as gastric homogenate.

Structure based designing of benzimidazole/benzoxazole derivatives as anti-leishmanial agents.

Folates are essential biomolecules required to carry out many crucial processes in leishmania parasite. Dihydrofolate reductase-thymidylate synthase (DHFR-TS) and pteridine reductase 1 (PTR1) involved in folate biosynthesis in leishmania have been established as suitable targets for development of chemotherapy against leishmaniasis. In the present study, various computational tools such as homology modelling, pharmacophore modelling, docking, molecular dynamics and molecular mechanics have been employed to design dual DHFR-TS and PTR1 inhibitors. Two designed molecules, i.e. 2-(4-((4-nitrobenzyl)oxy)phenyl)-1H-benzo[d]imidazole and 2-(4-((2,4-dichlorobenzyl)oxy)phenyl)-1H-benzo[d]oxazolemolecules were synthesized. MTT (3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyl tetrazolium bromide) assay was performed to evaluate in vitro activity of molecules against promastigote form of Leishmania donovani using Miltefosine as standard. 2-(4-((4-nitrobenzyl)oxy)phenyl)-1H-benzo[d]imidazole and 2-(4-((2,4-dichlorobenzyl)oxy)phenyl)-1H-benzo[d]oxazolemolecules were found to be moderately active with showed IC50 = 68 ± 2.8 µM and 57 ± 4.2 µM, respectively.

Molecular dynamics and pharmacophore modelling studies of different subtype (ALK and EGFR (T790M)) inhibitors in NSCLC.

Extensively validated 3D pharmacophore models for ALK (anaplastic lymphoma kinase) and EGFR (T790M) (epithelial growth factor receptor with acquired secondary mutation) were developed. The pharmacophore model for ALK (r2 = 0.96, q2 = 0.692) suggested that two hydrogen bond acceptors and three hydrophobic groups arranged in 3-D space are essential for the binding affinity of ALK inhibitors. Similarly, the pharmacophore model for EGFR (T790M) (r2 = 0.92, q2 = 0.72) suggested that the presence of a hydrogen bond acceptor, two hydrogen bond donors and a hydrophobic group plays vital role in binding of an inhibitor of EGFR (T790M). These pharmacophore models allowed searches for novel ALK and EGFR (T790M) dual inhibitors from multiconformer 3D databases (Asinex, Chembridge and Maybridge). Finally, the eight best hits were selected for molecular dynamics simulation, to study the stability of their complexes with both proteins and final binding orientations of these molecules. After molecular dynamics simulations, one hit has been predicted to possess good binding affinity for both ALK and EGFR (T790M), which can be further investigated for its experimental in-vitro/in-vivo activities.

Identification of new dual spleen tyrosine kinase (Syk) and phosphoionositide-3-kinase δ (PI3Kδ) inhibitors using ligand and structure-based integrated ideal pharmacophore models.

Owing to the complex pathophysiology of autoimmune disorders, it is very challenging to develop successful treatment strategies. Single-target agents are not desired therapeutics for such multi-factorial disorders. Considering the current need for the treatment of complex autoimmune disorders, dual inhibitors of Syk and PI3Kδ have been designed using ligand and structure-based molecular modelling strategies. In the present work, structure and ligand-based pharmacophore modelling was implemented for a varied set of Syk and PI3Kδ inhibitors. Ligand-based pharmacophore models (LBPMs) were developed for two kinases: ADPR.14 (r(2)train = 0.809) for Syk, comprising one hydrogen bond acceptor, one hydrogen bond donor, one positive ionisable and one ring aromatic feature, and for PI3Kδ: AAARR.45 (r(2)train = 0.942) consisting of three hydrogen bond acceptor and two ring aromatic features. The generated e-pharmacophore models revealed an additional ring aromatic and hydrophobic feature important for Syk and PI3Kδ inhibition, respectively. Subsequently, LBPMs were modified resulting in APDRR.14 hypothesis for Syk inhibitors and AAAHRR.45 hypothesis for PI3Kδ inhibitors employed for virtual screening. Thus, the combination of ligand and structure-based pharmacophore modelling helped in developing ideal pharmacophore models that may be an efficient tool for the designing of novel dual inhibitors of Syk and PI3Kδ.

Pharmacophore and docking-based virtual screening approach for the design of new dual inhibitors of Janus kinase 1 and Janus kinase 2.

Janus kinase 1 and 2, non-receptor protein tyrosine kinases, are implicated in various cancerous diseases. Involvement of these two enzymes in the pathways that stimulate cell proliferation in cancerous conditions makes them potential therapeutic targets for designing new dual-targeted agents for the treatment of cancer. In the present study, two separate pharmacophore models were developed and the best models for JAK1 (AAADH.25) and JAK2 (ADRR.92) were selected on the basis of their external predictive ability. Both models were subjected to a systematic virtual screening (VS) protocol using a PHASE database of 1.5 million molecules. The hits retrieved in VS were investigated for ADME properties to avoid selection of molecules with a poor pharmacokinetic profile. The molecules considered to be within the range of acceptable limits of ADME properties were further employed for docking simulations with JAK1 and JAK2 proteins to explore the final hits that possess structural features of both pharmacophore models as well as display essential interactions with both of them. Thus, the new molecules obtained in this way should show inhibitory activity against JAK1 and JAK2 and may serve as novel therapeutic agents for the treatment of cancerous disease conditions.

Imidazo[1,2-a]pyrazine inhibitors of phosphoinositide 3-kinase alpha (PI3Kα): 3D-QSAR analysis utilizing the Hybrid Monte Carlo algorithm to refine receptor-ligand complexes for molecular alignment.

Phosphoinositide 3-kinase alpha (PI3Kα) is a lipid kinase involved in several cellular functions such as cell growth, proliferation, differentiation and survival, and its anomalous regulation leads to cancerous conditions. PI3Kα inhibition completely blocks the cancer signalling pathway, hence it can be explored as an important therapeutic target for cancer treatment. In the present study, docking analysis of 49 selective imidazo[1,2-a]pyrazine inhibitors of PI3Kα was carried out using the QM-Polarized ligand docking (QPLD) program of the Schrödinger software, followed by the refinement of receptor-ligand conformations using the Hybrid Monte Carlo algorithm in the Liaison program, and alignment of refined conformations of inhibitors was utilized for the development of an atom-based 3D-QSAR model in the PHASE program. Among the five generated models, the best model was selected corresponding to PLS factor 2, displaying the highest value of Q(2)test (0.650). The selected model also displayed high values of r(2)train (0.917), F-value (166.5) and Pearson-r (0.877) and a low value of SD (0.265). The contour plots generated for the selected 3D-QSAR model were correlated with the results of docking simulations. Finally, this combined information generated from 3D-QSAR and docking analysis was used to design new congeners.

Glutamine: fructose-6-phosphate amidotransferase (GFAT): homology modelling and designing of new inhibitors using pharmacophore and docking based hierarchical virtual screening protocol.

Glutamine: fructose-6-phosphate amidotransferase (GFAT), also termed GFPT1 and GFAT1, catalyzes the first committed step of the hexosamine biosynthesis pathway in mammals and consequently plays an important role in type 2 diabetes. In the present study, a combination of pharmacophore modelling, homology modelling, and molecular docking analysis was performed to design new glutamine competitive inhibitors of human GFAT, and to investigate important interaction details of inhibitor molecules. A pharmacophore model of GFAT inhibitors was developed, subsequently validated, and utilized for the screening by the PHASE database to identify new molecules. Afterwards, homology modelling was performed to construct the glutamine-binding site of the GFAT protein. The modelled active site was utilized to dock the studied molecules to investigate important receptor-ligand interactions and to scrutinize database-screened molecules on the basis of essential interactions. This systematic in silico protocol helped us to identify new molecules that would be explored for the treatment of type 2 diabetes and its complications.

A three-dimensional pharmacophore modelling of ITK inhibitors and virtual screening for novel inhibitors.

Interleukin-2-inducible T-cell kinase (ITK) is a key member of the Tec family of non-receptor tyrosine kinases, and has been found to be a novel target for a number of inflammatory and autoimmune diseases. A three-dimensional pharmacophore model has been generated for protein ITK from its known inhibitors. The best HypoGen model consisted of four pharmacophore features: one hydrogen bond acceptor, one hydrogen bond donor and two hydrophobic rings. This model showed a correlation coefficient of 0.947, a root mean square deviation of 0.914 and a configuration cost of 16.866. The model was validated using test set prediction and Fischer's test. A test set containing 204 compounds showed an r(2) of 0.745 between estimated activity and activity measured experimentally. Fisher's test gave a confidence level of 95%. The best pharmacophore model (Hypo1) was then employed for virtual screening (3D database searching), including Lipinsiki's filter, to obtain a pool of more drug-like molecules. The molecular pool thus retrieved was subjected to docking analysis with a study protein to remove any molecules showing false positive activity for ITK.