Computational design, synthesis and evaluation of new sulphonamide derivatives targeting HIV-1 gp120.

Attachment of envelope glycoprotein gp120 to the host cell receptor CD4 is the first step during the human immunodeficiency virus-1 (HIV-1) entry into the host cells that makes it a promising target for drug design. To elucidate the crucial three dimensional (3D) structural features of reported HIV-1 gp120 CD4 binding inhibitors, 3D pharmacophores were generated and receptor based approach was employed to quantify these structural features. A four-partial least square factor model with good statistics and predictive ability was generated for the dataset of 100 molecules. To further ascertain the structural requirement for gp120-CD4 binding inhibition, molecular interaction studies of inhibitors with gp120 was carried out by performing molecular docking using Glide 5.6. Based on these studies, structural requirements were drawn and new molecules were designed accordingly to yield new sulphonamides derivatives. A water based green synthetic approach was adopted to obtain these compounds which were evaluated for their HIV-1 gp120 CD4 binding inhibition. The newly synthesized compounds exhibited remarkable activity (10-fold increase) when compared with the standard BMS 806. Further the stability of newly synthesized derivatives with HIV-1 gp120 was also investigated through molecular dynamics simulation studies. This provides a proof of concept for molecular modeling based design of new inhibitors for inhibition of HIV-1 gp120 CD4 interaction.

Structural insights of Staphylococcus aureus FtsZ inhibitors through molecular docking, 3D-QSAR and molecular dynamics simulations.

Filamentous temperature-sensitive protein Z (FtsZ) is a protein encoded by the FtsZ gene that assembles into a Z-ring at the future site of the septum of bacterial cell division. Structurally, FtsZ is a homolog of eukaryotic tubulin but has low sequence similarity; this makes it possible to obtain FtsZ inhibitors without affecting the eukaryotic cell division. Computational studies were performed on a series of substituted 3-arylalkoxybenzamide derivatives reported as inhibitors of FtsZ activity in Staphylococcus aureus. Quantitative structure-activity relationship models (QSAR) models generated showed good statistical reliability, which is evident from r2ncv and r2loo values. The predictive ability of these models was determined and an acceptable predictive correlation (r2Pred) values were obtained. Finally, we performed molecular dynamics simulations in order to examine the stability of protein-ligand interactions. This facilitated us to compare free binding energies of cocrystal ligand and newly designed molecule B1. The good concordance between the docking results and comparative molecular field analysis (CoMFA)/comparative molecular similarity indices analysis (CoMSIA) contour maps afforded obliging clues for the rational modification of molecules to design more potent FtsZ inhibitors.

Molecular docking, 3D-QSAR, molecular dynamics, synthesis and anticancer activity of tyrosine kinase 2 (TYK 2) inhibitors.

A therapeutic rationale is proposed by selectively targeting tyrosine kinase 2 (TYK 2) to obtain potent TYK 2 inhibitors by molecular modeling studies. In the present study, we have taken tyrosine kinase (TYK 2) inhibitors and carried out molecular docking, 3 D quantitative structure-activity relationship (3D-QSAR) analysis and molecular dynamics (MD). Based on the 3D-QSAR results thirteen new compounds (R-1 to R-13) were designed and synthesized in good yields. The synthesized molecules were evaluated for their in vitro anticancer activity against LnCap and A549 cell lines. The molecules R-1, R-3, R-5, R-7, and R-10 exhibited considerable anti cancer activity.

Molecular modeling-driven approach for identification of Janus kinase 1 inhibitors through 3D-QSAR, docking and molecular dynamics simulations.

Janus kinase 1 (JAK 1) belongs to the JAK family of intracellular nonreceptor tyrosine kinase. JAK-signal transducer and activator of transcription (JAK-STAT) pathway mediate signaling by cytokines, which control survival, proliferation and differentiation of a variety of cells. Three-dimensional quantitative structure activity relationship (3 D-QSAR), molecular docking and molecular dynamics (MD) methods was carried out on a dataset of Janus kinase 1(JAK 1) inhibitors. Ligands were constructed and docked into the active site of protein using GLIDE 5.6. Best docked poses were selected after analysis for further 3 D-QSAR analysis using comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) methodology. Employing 60 molecules in the training set, 3 D-QSAR models were generate that showed good statistical reliability, which is clearly observed in terms of r2ncv and q2loo values. The predictive ability of these models was determined using a test set of 25 molecules that gave acceptable predictive correlation (r2Pred) values. The key amino acid residues were identified by means of molecular docking, and the stability and rationality of the derived molecular conformations were also validated by MD simulation. The good consonance between the docking results and CoMFA/CoMSIA contour maps provides helpful clues about the reasonable modification of molecules in order to design more efficient JAK 1 inhibitors. The developed models are expected to provide some directives for further synthesis of highly effective JAK 1 inhibitors.

Computational analysis of ABL kinase mutations allows predicting drug sensitivity against selective kinase inhibitors.

The ABL kinase inhibitor imatinib has been used as front-line therapy for Philadelphia-positive chronic myeloid leukemia. However, a significant proportion of imatinib-treated patients relapse due to occurrence of mutations in the ABL kinase domain. Although inhibitor sensitivity for a set of mutations was reported, the role of less frequent ABL kinase mutations in drug sensitivity/resistance is not known. Moreover, recent reports indicate distinct resistance profiles for second-generation ABL inhibitors. We thus employed a computational approach to predict drug sensitivity of 234 point mutations that were reported in chronic myeloid leukemia patients. Initial validation analysis of our approach using a panel of previously studied frequent mutations indicated that the computational data generated in this study correlated well with the published experimental/clinical data. In addition, we present drug sensitivity profiles for remaining point mutations by computational docking analysis using imatinib as well as next generation ABL inhibitors nilotinib, dasatinib, bosutinib, axitinib, and ponatinib. Our results indicate distinct drug sensitivity profiles for ABL mutants toward kinase inhibitors. In addition, drug sensitivity profiles of a set of compound mutations in ABL kinase were also presented in this study. Thus, our large scale computational study provides comprehensive sensitivity/resistance profiles of ABL mutations toward specific kinase inhibitors.

An integrated molecular modeling approach for in silico design of new tetracyclic derivatives as ALK inhibitors.

Anaplastic lymphoma kinase (ALK), a promising therapeutic target for treatment of human cancers, is a receptor tyrosine kinase that instigates the activation of several signal transduction pathways. In the present study, in silico methods have been employed in order to explore the structural features and functionalities of a series of tetracyclic derivatives displaying potent inhibitory activity toward ALK. Initially docking was performed using GLIDE 5.6 to probe the bioactive conformation of all the compounds and to understand the binding modes of inhibitors. The docking results revealed that ligand interaction with Met 1199 plays a crucial role in binding of inhibitors to ALK. Further to establish a robust 3D-QSAR model using CoMFA and CoMSIA methods, the whole dataset was divided into three splits. Model obtained from Split 3 showed high accuracy ([Formula: see text] of 0.700 and 0.682, [Formula: see text] of 0.971 and 0.974, [Formula: see text] of 0.673 and 0.811, respectively for CoMFA and CoMSIA). The key structural requirements for enhancing the inhibitory activity were derived from CoMFA and CoMSIA contours in combination with site map analysis. Substituting small electronegative groups at Position 8 by replacing either morpholine or piperidine rings and maintaining hydrophobic character at Position 9 in tetracyclic derivatives can enhance the inhibitory potential. Finally, we performed molecular dynamics simulations in order to investigate the stability of protein ligand interactions and MM/GBSA calculations to compare binding free energies of co-crystal ligand and newly designed molecule N1. Based on the coherence of outcome of various molecular modeling studies, a set of 11 new molecules having potential predicted inhibitory activity were designed.

3D QSAR based design of novel oxindole derivative as 5HT7 inhibitors.

To understand the structural requirements of 5-hydroxytryptamine (5HT7) receptor inhibitors and to design new ligands against 5HT7 receptor with enhanced inhibitory potency, a three-dimensional quantitative structure-activity relationship study with comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) for a data set of 56 molecules consisting of oxindole, tetrahydronaphthalene, aryl ketone substituted arylpiperazinealkylamide derivatives was performed. Derived model showed good statistical reliability in terms of predicting 5HT7 inhibitory activity of the molecules, based on molecular property fields like steric, electrostatic, hydrophobic, hydrogen bond donor and hydrogen bond acceptor fields. This is evident from statistical parameters like conventional r2 and a cross validated (q2) values of 0.985, 0.743 for CoMFA and 0.970, 0.608 for CoMSIA, respectively. Predictive ability of the models to determine 5HT7 antagonistic activity is validated using a test set of 16 molecules that were not included in the training set. Predictive r2 obtained for the test set was 0.560 and 0.619 for CoMFA and CoMSIA, respectively. Steric, electrostatic fields majorly contributed toward activity which forms the basis for design of new molecules. Absorption, distribution, metabolism and elimination (ADME) calculation using QikProp 2.5 (Schrodinger 2010, Portland, OR) reveals that the molecules confer to Lipinski's rule of five in majority of the cases.

Multiple receptor conformation docking, dock pose clustering and 3D QSAR studies on human poly(ADP-ribose) polymerase-1 (PARP-1) inhibitors.

Poly(ADP-ribose) polymerase-1 (PARP-1) functions as a DNA damage sensor and signaling molecule. It plays a vital role in the repair of DNA strand breaks induced by radiation and chemotherapeutic drugs; inhibitors of this enzyme have the potential to improve cancer chemotherapy or radiotherapy. Three-dimensional quantitative structure activity relationship (3D QSAR) models were developed using comparative molecular field analysis, comparative molecular similarity indices analysis and docking studies. A set of 88 molecules were docked into the active site of six X-ray crystal structures of poly(ADP-ribose)polymerase-1 (PARP-1), by a procedure called multiple receptor conformation docking (MRCD), in order to improve the 3D QSAR models through the analysis of binding conformations. The docked poses were clustered to obtain the best receptor binding conformation. These dock poses from clustering were used for 3D QSAR analysis. Based on MRCD and QSAR information, some key features have been identified that explain the observed variance in the activity. Two receptor-based QSAR models were generated; these models showed good internal and external statistical reliability that is evident from the [Formula: see text], [Formula: see text] and [Formula: see text]. The identified key features enabled us to design new PARP-1 inhibitors.

Molecular docking guided structure based design of symmetrical N,N'-disubstituted urea/thiourea as HIV-1 gp120-CD4 binding inhibitors.

Induced fit molecular docking studies were performed on BMS-806 derivatives reported as small molecule inhibitors of HIV-1 gp120-CD4 binding. Comprehensive study of protein-ligand interactions guided in identification and design of novel symmetrical N,N'-disubstituted urea and thiourea as HIV-1 gp120-CD4 binding inhibitors. These molecules were synthesized in aqueous medium using microwave irradiation. Synthesized molecules were screened for their inhibitory ability by HIV-1 gp120-CD4 capture enzyme-linked immunosorbent assay (ELISA). Designed compounds were found to inhibit HIV-1 gp120-CD4 binding in micromolar (0.013-0.247 µM) concentrations.

Synthesis, Antifungal Ergosterol Inhibition, Antibiofilm Activities, and Molecular Docking on ß-Tubulin and Sterol 14-Alpha Demethylase along with DFT-Based Quantum Mechanical Calculation of Pyrazole Containing Fused Pyridine-Pyrimidine Derivatives.

Multidrug-resistant fungal infections have become much more common in recent years, especially in immune-compromised patients. Therefore, researchers and pharmaceutical professionals have focused on the development of novel antifungal agents that can tackle the problem of resistance. In continuation to this, a novel series of pyrazole-bearing pyrido[2,3-d]pyrimidine-2,4(1H,3H)-dione derivatives (4a-4o) have been developed. These compounds have been screened against Candida albicans, Aspergillus niger, and Aspergillus clavatus. The synthesized compounds were characterized by well-known spectroscopic techniques, i.e., IR, 1H NMR, 13C NMR, and mass spectrometry. In vitro antifungal results revealed that compound 4n showed activity against C. albicans having MIC value of 200 µg/mL. To know the plausible mode of action, the active derivatives were screened for anti-biofilm and ergosterol biosynthesis inhibition activities. The compounds 4h, 4j, 4k, and 4n showed greater ergosterol biosynthesis inhibition than the control DMSO. To comprehend how molecules interact with the receptor, studies of molecular docking of 4k and 4n have been performed on the homology-modeled protein of ß-tubulin. The molecular docking revealed that the active compounds 4h, 4j, 4k, 4l, and 4n interacting with the active site amino acid of sterol 14-alpha demethylase (PDB ID: 5v5z) indicate one of the possible modes of action of ergosterol inhibition activity. The synthesized compounds 4c, 4e, 4h, 4i, 4j, 4k, 4l, and 4n inhibited biofilm formation and possessed the potential for anti-biofilm activity. DFT-based quantum mechanical calculations were carried out to optimize, predict, and compare the vibration modes of the molecule 4a.

Molecular docking and 3D-QSAR studies on inhibitors of DNA damage signaling enzyme human PARP-1.

Poly (ADP-ribose) polymerase-1 (PARP-1) operates in a DNA damage signaling network. Molecular docking and three dimensional-quantitative structure activity relationship (3D-QSAR) studies were performed on human PARP-1 inhibitors. Docked conformation obtained for each molecule was used as such for 3D-QSAR analysis. Molecules were divided into a training set and a test set randomly in four different ways, partial least square analysis was performed to obtain QSAR models using the comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA). Derived models showed good statistical reliability that is evident from their r², q²(loo) and r²(pred) values. To obtain a consensus for predictive ability from all the models, average regression coefficient r²(avg) was calculated. CoMFA and CoMSIA models showed a value of 0.930 and 0.936, respectively. Information obtained from the best 3D-QSAR model was applied for optimization of lead molecule and design of novel potential inhibitors.

Identification of potential Aurora kinase-C protein inhibitors: an amalgamation of energy minimization, virtual screening, prime MMGBSA and AutoDock.

Genomic instability is a trademark of cancer evolution, there is still a wide necessity to discover safe and effective anti-cancer drugs. Aurora kinase-C protein is a member of the Aurora kinase family involved in regulating mitosis phases of the cell cycle. Aurora kinase-C protein has an aberrant expression at spindle assembly checkpoint leading to human cervical cancer and colorectal cancer. Knowledge of 3D structure of the protein is vital for identify possible inhibitors. In the present study, 3D structure prediction of Aurora kinase-C protein is taken up by using comparative homology modeling techniques and minimized by NAMD-VMD Software. The quality of modeled protein was validated using ProSA and Ramachandran plot. Probable binding site in protein was identified from SiteMap and molecular docking based virtual screening was carried out using Asinex database to identify lead molecules. The obtained lead molecules were further optimized using PrimeMM-GBSA and AutoDock. ADME properties of the leads were calculated to ascertain the drug-like properties at lead molecules.Communicated by Ramaswamy H. Sarma.

Rational design of methicillin resistance staphylococcus aureus inhibitors through 3D-QSAR, molecular docking and molecular dynamics simulations.

Staphylococcus aureus is a gram positive bacterium. It is the leading cause of skin and respiratory infections, osteomyelitis, Ritter's disease, endocarditis, and bacteraemia in the developed world. We employed combined studies of 3D QSAR, molecular docking which are validated by molecular dynamics simulations and in silico ADME prediction have been performed on Isothiazoloquinolones inhibitors against methicillin resistance Staphylococcus aureus. Three-dimensional quantitative structure-activity relationship (3D-QSAR) study was applied using comparative molecular field analysis (CoMFA) with Q2 of 0.578, R2 of 0.988, and comparative molecular similarity indices analysis (CoMSIA) with Q2 of 0.554, R2 of 0.975. The predictive ability of these model was determined using a test set of molecules that gave acceptable predictive correlation (r2 Pred) values 0.55 and 0.57 of CoMFA and CoMSIA respectively. Docking, simulations were employed to position the inhibitors into protein active site to find out the most probable binding mode and most reliable conformations. Developed models and Docking methods provide guidance to design molecules with enhanced activity.

Molecular dynamics and MM/GBSA-integrated protocol probing the correlation between biological activities and binding free energies of HIV-1 TAR RNA inhibitors.

The interaction of HIV-1 transactivator protein Tat with its cognate transactivation response (TAR) RNA has emerged as a promising target for developing antiviral compounds and treating HIV infection, since it is a crucial step for efficient transcription and replication. In the present study, molecular dynamics (MD) simulations and MM/GBSA calculations have been performed on a series of neamine derivatives in order to estimate appropriate MD simulation time for acceptable correlation between ΔGbind and experimental pIC50 values. Initially, all inhibitors were docked into the active site of HIV-1 TAR RNA. Later to explore various conformations and examine the docking results, MD simulations were carried out. Finally, binding free energies were calculated using MM/GBSA method and were correlated with experimental pIC50 values at different time scales (0-1 to 0-10 ns). From this study, it is clear that in case of neamine derivatives as simulation time increased the correlation between binding free energy and experimental pIC50 values increased correspondingly. Therefore, the binding energies which can be interpreted at longer simulation times can be used to predict the bioactivity of new neamine derivatives. Moreover, in this work, we have identified some plausible critical nucleotide interactions with neamine derivatives that are responsible for potent inhibitory activity. Furthermore, we also provide some insights into a new class of oxadiazole-based back bone cyclic peptides designed by incorporating the structural features of neamine derivatives. On the whole, this approach can provide a valuable guidance for designing new potent inhibitors and modify the existing compounds targeting HIV-1 TAR RNA.

Computational Analysis of Epidermal Growth Factor Receptor Mutations Predicts Differential Drug Sensitivity Profiles toward Kinase Inhibitors.

INTRODUCTION: A significant proportion of patients with lung cancer carry mutations in the EGFR kinase domain. The presence of a deletion mutation in exon 19 or L858R point mutation in the EGFR kinase domain has been shown to cause enhanced efficacy of inhibitor treatment in patients with NSCLC. Several less frequent (uncommon) mutations in the EGFR kinase domain with potential implications in treatment response have also been reported. The role of a limited number of uncommon mutations in drug sensitivity was experimentally verified. However, a huge number of these mutations remain uncharacterized for inhibitor sensitivity or resistance. METHODS: A large-scale computational analysis of clinically reported 298 point mutants of EGFR kinase domain has been performed, and drug sensitivity profiles for each mutant toward seven kinase inhibitors has been determined by molecular docking. In addition, the relative inhibitor binding affinity toward each drug as compared with that of adenosine triphosphate was calculated for each mutant. RESULTS: The inhibitor sensitivity profiles predicted in this study for a set of previously characterized mutants correlated well with the published clinical, experimental, and computational data. Both the single and compound mutations displayed differential inhibitor sensitivity toward first- and next-generation kinase inhibitors. CONCLUSIONS: The present study provides predicted drug sensitivity profiles for a large panel of uncommon EGFR mutations toward multiple inhibitors, which may help clinicians in deciding mutant-specific treatment strategies.

Identification of novel anti cancer agents by applying insilico methods for inhibition of TSPO protein.

Cancer is a genomic disease characterised as impaired cellular energy metabolism. Cancer cells derive most of their energy from oxidative phosphorylation unlike normal ones during cell progression TSPO protein present in external mitochondrial membrane, is involved in various cellular functions like Cell proliferation, mitochondrial respiration, synthesis of steroids and also participates in import of cholesterol into the inner mitochondrial membrane from outside of the membrane of mitochondria. The 3D model of TSPO protein is built using comparative homology modelling techniques and validated by proSA, Ramachandran plot and ERRAT in the present work. Active site prediction is carried out using SiteMap and literature, which allows the prediction of the important binding pockets for the identification of putative active site. New molecular entities as TSPO inhibitors were obtained from Virtual screening using MS Spectrum databank in Schrodinger suite and were prioritised based on Glide Score. Docking was performed using Autodock to identify molecules with different scaffolds and were prioritised based on binding energy and RMSD values. Qikprop is used to calculate pharmacokinetic properties of the screened molecules which are found to be in permissible range as possible novel inhibitors of TSPO protein to supress cell proliferation.

Molecular docking, 3D QSAR and dynamics simulation studies of imidazo-pyrrolopyridines as janus kinase 1 (JAK 1) inhibitors.

Janus kinase 1 (JAK 1) plays a critical role in initiating responses to cytokines by the JAK-signal transducer and activator of transcription (JAK-STAT). This controls survival, proliferation and differentiation of a variety of cells. Docking, 3D quantitative structure activity relationship (3D-QSAR) and molecular dynamics (MD) studies were performed on a series of Imidazo-pyrrolopyridine derivatives reported as JAK 1 inhibitors. QSAR model was generated using 30 molecules in the training set; developed model showed good statistical reliability, which is evident from r2ncv and r2loo values. The predictive ability of this model was determined using a test set of 13 molecules that gave acceptable predictive correlation (r2Pred) values. Finally, molecular dynamics simulation was performed to validate docking results and MM/GBSA calculations. This facilitated us to compare binding free energies of cocrystal ligand and newly designed molecule R1. The good concordance between the docking results and CoMFA/CoMSIA contour maps afforded obliging clues for the rational modification of molecules to design more potent JAK 1 inhibitors.

Multiple receptor conformation docking and dock pose clustering as tool for CoMFA and CoMSIA analysis - a case study on HIV-1 protease inhibitors.

Multiple receptors conformation docking (MRCD) and clustering of dock poses allows seamless incorporation of receptor binding conformation of the molecules on wide range of ligands with varied structural scaffold. The accuracy of the approach was tested on a set of 120 cyclic urea molecules having HIV-1 protease inhibitory activity using 12 high resolution X-ray crystal structures and one NMR resolved conformation of HIV-1 protease extracted from protein data bank. A cross validation was performed on 25 non-cyclic urea HIV-1 protease inhibitor having varied structures. The comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) models were generated using 60 molecules in the training set by applying leave one out cross validation method, r (loo) (2) values of 0.598 and 0.674 for CoMFA and CoMSIA respectively and non-cross validated regression coefficient r(2) values of 0.983 and 0.985 were obtained for CoMFA and CoMSIA respectively. The predictive ability of these models was determined using a test set of 60 cyclic urea molecules that gave predictive correlation (r (pred) (2) ) of 0.684 and 0.64 respectively for CoMFA and CoMSIA indicating good internal predictive ability. Based on this information 25 non-cyclic urea molecules were taken as a test set to check the external predictive ability of these models. This gave remarkable out come with r (pred) (2) of 0.61 and 0.53 for CoMFA and CoMSIA respectively. The results invariably show that this method is useful for performing 3D QSAR analysis on molecules having different structural motifs.

Molecular docking and 3D-QSAR studies on triazolinone and pyridazinone, non-nucleoside inhibitor of HIV-1 reverse transcriptase.

Nonnucleoside reverse transcriptase inhibitors (NNRTIs) are allosteric inhibitors of the HIV-1 reverse transcriptase. Recently a series of Triazolinone and Pyridazinone were reported as potent inhibitors of HIV-1 wild type reverse transcriptase. In the present study, docking and 3D quantitative structure activity relationship (3D QSAR) studies involving comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) were performed on 31 molecules. Ligands were built and minimized using Tripos force field and applying Gasteiger-Hückel charges. These ligands were docked into protein active site using GLIDE 4.0. The docked poses were analyzed; the best docked poses were selected and aligned. CoMFA and CoMSIA fields were calculated using SYBYL6.9. The molecules were divided into training set and test set, a PLS analysis was performed and QSAR models were generated. The model showed good statistical reliability which is evident from the r2 nv, q2 loo and r2 pred values. The CoMFA model provides the most significant correlation of steric and electrostatic fields with biological activities. The CoMSIA model provides a correlation of steric, electrostatic, acceptor and hydrophobic fields with biological activities. The information rendered by 3D QSAR model initiated us to optimize the lead and design new potential inhibitors.