Pharmacophore generation,atom-based 3D-QSAR,molecular docking and molecular dynamics simulation studies on benzamide analogues as FtsZ inhibitors

FtsZ is an appealing target for the design of antimicrobial agent that can be used to defeat the multidrug-resistant bacterial pathogens. Pharmacophore modelling, molecular docking and molecular dynamics (MD) simulation studies were performed on a series of three-substituted benzamide derivatives. In the present study a five-featured pharmacophore model with one hydrogen bond acceptors, one hydrogen bond donors, one hydrophobic and two aromatic rings was developed using 97 molecules having MIC values ranging from .07 to 957 μM. A statistically significant 3D-QSAR model was obtained using this pharmacophore hypothesis with a good correlation coefficient (R² = .8319), cross validated coefficient (Q² = .6213) and a high Fisher ratio (F = 103.9) with three component PLS factor. A good correlation between experimental and predicted activity of the training (R² = .83) and test set (R² = .67) molecules were displayed by ADHRR.1682 model. The generated model was further validated by enrichment studies using the decoy test and MAE-based criteria to measure the efficiency of the model. The docking studies of all selected inhibitors in the active site of FtsZ protein showed crucial hydrogen bond interactions with Val 207, Asn 263, Leu 209, Gly 205 and Asn-299 residues. The binding free energies of these inhibitors were calculated by the molecular mechanics/generalized born surface area VSGB 2.0 method. Finally, a 15 ns MD simulation was done to confirm the stability of the 4DXD–ligand complex. On a wider scope, the prospect of present work provides insight in designing molecules with better selective FtsZ inhibitory potential.     

A combination of pharmacophore modeling,atom-based 3D-QSAR,molecular docking and molecular dynamics simulation studies on PDE4 enzyme inhibitors

Phosphodiesterases 4 enzyme is an attractive target for the design of anti-inflammatory and bronchodilator agents. In the present study, pharmacophore and atom-based 3D-QSAR studies were carried out for pyrazolopyridine and quinoline derivatives using Schrödinger suite 2014-3. A four-point pharmacophore model was developed using 74 molecules having pIC₅₀ ranging from 10.1 to 4.5. The best four feature model consists of one hydrogen bond acceptor, two aromatic rings, and one hydrophobic group. The pharmacophore hypothesis yielded a statistically significant 3D-QSAR model, with a high correlation coefficient (R^2?=?.9949), cross validation coefficient (Q^2?=?.7291), and Pearson-r (.9107) at six component partial least square factor. The external validation indicated that our QSAR model possessed high predictive power with R² value of .88. The generated model was further validated by enrichment studies using the decoy test. Molecular docking, free energy calculation, and molecular dynamics (MD) simulation studies have been performed to explore the putative binding modes of these ligands. A 10-ns MD simulation confirmed the docking results of both stability of the 1XMU–ligand complex and the presumed active conformation. Outcomes of the present study provide insight in designing novel molecules with better PDE4 inhibitory activity.     

Pharmacophore models generation by catalyst and phase consensus-based virtual screening protocol against PI3Kα inhibitors

Phosphoinositide 3-kinases (PI3Ks) family has emerged as promising targets for novel therapeutic agents against neoplastic diseases. Pharmacophore and 3D-quantitative structure–activity relationship modelling were applied to study the structure–activity relationship of PI3K inhibitors. The best HypoGen pharmacophore hypothesis Hypo1 with a correlation coefficient of 0.961 consists of one hydrogen-bond acceptor, one hydrogen-bond donor and two hydrophobic features, whereas the best phase hypothesis AADRRR.378 with favourable statistics (q² = 0.7368, r² = 0.9863) has two hydrogen-bond acceptors, one hydrogen-bond donor and three ring aromatic features. Multiple methods, such as Fischer validation, molecular docking and mapping of test set molecules, were carried out to validate these pharmacophore models. Furthermore, a comparative molecular similarity indices analysis candidate hypothesis model was generated as a supplement of pharmacophore hypothesis. Detailed protein–ligand binding information obtained by Glide was utilised in compound optimisation and virtual screening. A molecular database of 133 known inhibitors and 6179 decoys was built for a screening test to quantitatively analyse various hypotheses and scoring parameters. Finally, we designed a workflow integrating HypoGen pharmacophore searching, phase pharmacophore searching and molecular docking for screening the database. With an improved criterion of enrichment factor (EF = 17.43) and ROC curve (AUC = 0.946), this workflow would provide us an original method for novel PI3K inhibitors.     

Ligand-based pharmacophore modeling and docking studies on vitamin D receptor inhibitors

In recent years, pharmacophore modeling and molecular docking approaches have been extensively used to characterize the structural requirements and explore the conformational space of a ligand in the binding pocket of the selected target protein. Herein, we report a pharmacophore modeling and molecular docking of 45 compounds comprising of the indole scaffold as vitamin D receptor (VDR) inhibitors. Based on the selected best hypothesis (DRRRR.61), an atom-based three-dimensional quantitative structure-activity relationships model was developed to rationalize the structural requirement of biological activity modulating components. The developed model predicted the binding affinity for the training set and test set with R²_(training) = 0.8869 and R²_(test) = 0.8139, respectively. Furthermore, molecular docking and dynamics simulation were performed to understand the underpinning of binding interaction and stability of selected VDR inhibitors in the binding pocket. In conclusion, the results presented here, in the form of functional and structural data, agreed well with the proposed pharmacophores and provide further insights into the development of novel VDR inhibitors with better activity.     

A profound computational study to prioritize the natural compound inhibitors against the P. falciparum orotidine-5-monophosphate decarboxylase enzyme

Abstract

In the current contribution, a multicomplex-based pharmacophore modeling approach was employed on the structural proteome of Plasmodium falciparum orotidine-5-monophosphate decarboxylase enzyme (PfOMPDC). Among the constructed pharmacophore models, the representative hypotheses were selected as the primary filter to screen the molecules with the complementary features responsible for showing inhibition. Thereafter, auxiliary evaluations were performed on the screened candidates via drug-likeness and molecular docking studies. Subsequently, the stability of the docked protein-ligand complexes was scrutinized by employing molecular dynamics simulations and molecular mechanics-Poisson Boltzmann surface area based free binding energy calculations. The stability the docked candidates was compared with the highly active crystallized inhibitor (3S9Y-FNU) to seek more potential candidates. All the docked molecules displayed stable dynamic behavior and high binding free energy in comparison to 3S9Y-FNU. The employed workflow resulted in the retrieval of five drug-like candidates with diverse scaffolds that may show inhibitory activity against PfOMPDC and could be further used as the novel scaffold to develop novel antimalarials.

Communicated by Ramaswamy H. Sarma     

Identification of InhA inhibitors: A combination of virtual screening,molecular dynamics simulations and quantum chemical studies

In the present work, multiple pharmacophore-based virtual screening of the SPECS natural product database was carried out to identify novel inhibitors of the validated biological target, InhA. The pharmacophore models were built from the five different groups of the co-crystallized ligands present within the active site. The generated models with the same features from each group were pooled and subjected to the test set validation, receiver–operator characteristic analysis and Güner–Henry studies. A set of five hypotheses with sensitivity > 0.5, specificity > 0.5, area under curve (AUC) > 0.7, and goodness of hit score > 0.7 were retrieved and exploited for the virtual screening. The common hits (87 molecules) obtained from these hypotheses were processed via drug-likeness filters. The filtered molecules (27 molecules) were compared for the binding modes and the top scored molecules (12 molecules) along with the reference (triclosan (TCL), docking score = ?11.65 kcal/mol) were rescored and reprioritized via molecular mechanics-generalized Born surface area approach. Eventually, the stability of reprioritized (10 molecules) docked complexes was scrutinized via molecular dynamics simulations. Moreover, the quantum chemical studies of the dynamically stable compounds (9 molecules) were performed to understand structural features essential for the activity. Overall, the protocol resulted in the recognition of nine lead compounds that can be targeted against InhA.     

Exploration and validation of diphosphate‐based Plasmodium LytB inhibitors using computational approaches

LytB or IspH is an indispensable enzyme and a suitable drug target of Plasmodium falciparum that participate in isoprenoid biosynthesis of nonmevalonate pathway (MEP). Recently, we have investigated the structural dynamics of Plasmodium LytB and proposed some novel diphosphate‐based inhibitors using molecular modeling and docking studies. Here, we have tried to characterize those previously screened molecules by quantitative structure activity relationships and pharmacophore‐based analyses, as well as validated the dynamics of their interactions with LytB protein. Five total compounds having PubChem CID 516 , 125696 , 165275 , 448012 , and 9921431 were predicted with significant inhibitory activity by quantitative structure activity relationships and pharmacophore models. Again, the molecular dynamics simulation results showed that these five compounds are able to form stable complexes with the receptor through many direct and water mediated interactions. The binding free energies calculated by Poisson‐Boltzmann surface area method resulted within the range between ?99.77 and ? 43.74 kcal/mol, which favoured their profound inhibitory affinity. Residues of LytB like His41, His74, Ser222, Ser223, and Asn224 in LytB were the main protagonists in contributing the majority of interaction energies to the ligands. Finally, the ADMET, toxicity, and drug‐likeness scores also affirm these compounds to be considered for further development of new antimalarial inhibitor in the future.     

Structural insights of JAK2 inhibitors: pharmacophore modeling and ligand-based 3D-QSAR studies of pyrido-indole derivatives

Abstract

In this study we have performed pharmacophore modeling and built a 3D QSAR model for pyrido-indole derivatives as Janus Kinase 2 inhibitors. An efficient pharmacophore has been identified from a data set of 51 molecules and the identified pharmacophore hypothesis consisted of one hydrogen bond acceptor, two hydrogen bond donors and three aromatic rings, i.e. ADDRRR. A powerful 3D-QSAR model has also been constructed by employing Partial Least Square regression analysis with a regression coefficient of 0.97 (R²) and Q² of 0.95, and Pearson-R of 0.98.     

Identification of novel urease inhibitors: pharmacophore modeling,virtual screening and molecular docking studies

Abstract

Pharmacophore modeling and atom-based three-dimensional quantitative structure–activity relationship (3D-QSAR) have been developed on N-acylglycino- and hippurohydroxamic acid derivatives, which are known potential inhibitors of urease. This is followed by virtual screening and ADMET (absorption, distribution, metabolism, excretion and toxicity) studies on a large library of known drugs in order to get lead molecules as Helicobacter pylori urease inhibitors. A suitable three-featured pharmacophore model comprising one H-bond acceptor and two H-bond donor features (ADD.10) has been found to be the best QSAR model. An external library of compounds (～3000 molecules), pre-filtered using Lipinski’s rule of five, has been further screened using the pharmacophore model ADD.10. By analyzing the fitness of the hits with respect to the pharmacophore model and their binding interaction inside the urease active site, four molecules have been predicted to be extremely good urease inhibitors. Two of these have significant potential and should be taken up for further drug-designing process.     

Pharmacophore mapping of a series of pyrrolopyrimidines,indolopyrimidines and their congeners as multidrug-resistance-associated protein (MRP1) modulators

Pharmacophore mapping studies were undertaken for a series of molecules belonging to pyrrolopyrimidines, indolopyrimidines and their congeners as multidrug resistance-associated protein (MRP1) modulators. A five-point pharmacophore with two hydrogen bond acceptors (A), one lipophilic/hydrophobic group (H), one positive ionic feature (P) and one aromatic ring (R) as pharmacophoric features was developed. The pharmacophore hypothesis yielded a statistically significant 3D-QSAR model, with a correlation coefficient of r ² = 0.799 for training set molecules. The model generated showed excellent predictive power, with a correlation coefficient Q ² = 0.679 for an external test set of 20 molecules. The pharmacophore was further validated using four structurally diverse compounds with MRP1 modulatory activity. These compounds mapped well onto four of the five features of the pharmacophore. The pharmacophore proposed here was then utilised for the successful retrieval of active molecules with diverse chemotypes from database search. The geometry and features of pharmacophore are expected to be useful for the design of selective MRP1 inhibitors. Figure Alignment of multidrug resistance-associated protein (MRP1) inhibitors with the developed pharmacophore. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

In silico studies on potential MCF-7 inhibitors: a combination of pharmacophore and 3D-QSAR modeling,virtual screening,molecular docking,and pharmacokinetic analysis

Gallic acid and its derivatives exhibit a diverse range of biological applications, including anti-cancer activity. In this work, a data-set of forty-six molecules containing the galloyl moiety, and known to show anticarcinogenic activity against the MCF-7 human cancer cell line, have been chosen for pharmacophore modeling and 3D-Quantitative Structure Activity Relationship (3D-QSAR) studies. A tree-based partitioning algorithm has been used to find common pharmacophore hypotheses. The QSAR model was generated for three, four, and five featured hypotheses with increasing PLS factors and analyzed. Results for five featured hypotheses with three acceptors and two aromatic rings were the best out of all the possible combinations. On analyzing the results, the most robust (R^2?=?.8990) hypothesis with a good predictive power (Q^2?=?.7049) was found to be AAARR.35. A good external validation (R² = .6109) was also obtained. In order to design new MCF-7 inhibitors, the QSAR model was further utilized in pharmacophore-based virtual screening of a large database. The predicted IC₅₀ values of the identified potential MCF-7 inhibitors were found to lie in the micromolar range. Molecular docking into the colchicine domain of tubulin was performed in order to examine one of the probable mechanisms. This revealed various interactions between the ligand and the active site protein residues. The present study is expected to provide an effective guide for methodical development of potent MCF-7 inhibitors.     

Insight into the structural requirements of thiophene-3-carbonitriles-based MurF inhibitors by 3D-QSAR,molecular docking and molecular dynamics study

The discovery of clinically relevant inhibitors against MurF enzyme has proven to be a challenging task. In order to get further insight into the structural features required for the MurF inhibitory activity, we performed pharmacophore and atom-based three-dimensional quantitative structure–activity relationship studies for novel thiophene-3-carbonitriles based MurF inhibitors. The five-feature pharmacophore model was generated using 48 inhibitors having IC₅₀ values ranging from 0.18 to 663?μm. The best-fitted model showed a higher coefficient of determination (R²?=?0.978), cross-validation coefficient (Q²?=?0.8835) and Pearson coefficient (0.9406) at four component partial least-squares factor. The model was validated with external data set and enrichment study. The effectiveness of the docking protocol was validated by docking the co-crystallized ligand into the catalytic pocket of MurF enzyme. Further, binding free energy calculated by the molecular mechanics generalized Born surface area approach showed that van der Waals and non-polar solvation energy terms are the main contributors to ligand binding in the active site of MurF enzyme. A 10-ns molecular dynamic simulation was performed to confirm the stability of the 3ZM6-ligand complex. Four new molecules are also designed as potent MurF inhibitors. These results provide insights regarding the development of novel MurF inhibitors with better binding affinity.     

Molecular docking and pharmacophore model studies of Rho kinase inhibitors

Molecular docking and pharmacophore model approaches were used to characterise the binding features of four different series of Rho kinase (ROCK) inhibitors. Docking simulation of 20 inhibitors with ROCK was performed. The binding conformations and binding affinities of these inhibitors were obtained using AutoDock 4.0 software. The predicted binding affinities correlate well with the activities of these inhibitors (R ² = 0.904). 3D pharmacophore models were generated for ROCK based on highly active inhibitors implemented in Catalyst 4.11 program. The best pharmacophore model consists of one hydrogen bond acceptor feature and two hydrophobic features, and they all seemed to be essential for inhibitors in terms of their binding activities. It is anticipated that the findings reported in this paper may provide very useful information for designing new ROCK inhibitors.     

Pharmacophore mapping,molecular docking and QSAR studies of structurally diverse compounds as CYP2B6 inhibitors

Pharmacophore mapping, molecular docking and quantitative structure–activity relationship (QSAR) studies were carried out for a structurally diverse set of 48 compounds as CYP2B6 inhibitors. The generated best pharmacophore hypotheses from the three methods of conformer generation (FAST, BEST and conformer algorithm based on energy screening and recursive buildup) indicate the importance of two features, namely, hydrogen bond acceptor [electron-rich centre] and ring aromaticity. The distance between the two centres of the important features for ideal inhibitors varied from 5.82 to 6.03 Å. The chemometric tools used for the QSAR analysis were genetic function approximation (GFA) and genetic partial least squares. The developed QSAR models indicate the importance of an electron-rich centre, size of molecule, impact of branching and ring system and distribution of charges in the molecular surface. The docking study confirms the importance of an electron-rich centre for binding with the iron atom of the cytochrome enzyme. A GFA model with spline option was found to be the best model based on internal validation as well as the r ² _{m (overall)} criterion (Q ² = 0.772, r ² _{m (overall)} = 0.774). According to the external prediction statistics (R ² _pred = 0.876), another GFA-derived model with spline option outperforms the remaining models.     

Design of antimalarial transmission blocking agents: Pharmacophore mapping of ligands active against stage-V mature gametocytes of Plasmodium falciparum

Abstract

The discovery of transmission-blocking (T-B) agents is crucial for preventing and complete removal of malaria infection. However, most of the existing antimalarials are only active against the asexual stages of Plasmodium parasite, but ineffective against the sexual stage (gametocytes). In this background, we have developed pharmacophore models against the stage-V mature gametocytes of P. falciparum parasites. The pharmacophore model (Hypo-1) showed five pharmacophoric features namely, one hydrogen bond donor (HBD), one hydrophobic aliphatic (HYAl), one ring aromatic (RA), and two hydrophobic aromatic (HYAr) essential for the anti-gametocytic activity. The amino, methyl, fused phenyl ring of the quinazoline heterocycle, two phenyl rings of biphenyl moiety (HBD, HYAl, HYAr1, HYAr2 and RA) are the crucial features responsible for the non-specific anti-gametocytic activity (PfG). Subsequently, the model (Hypo-2) developed against the stage-V female gametocytes (PffG) showed the contribution of three pharmacophoric features namely, two hydrogen bond acceptor (HYA) and one RA required for the anti-gametocytic activity. The sulfhydryl, imine and pyridyl groups are observed to be essential for anti-gametocytic activity against female gametocytes. Both the models (PfG and PfGG) showed the classification accuracies of 78.26 and 71.64% for training set compounds and 60.80 and 60.18% for the test set compounds, respectively, for classification of compounds into higher and lower active classes. Also, both the models were found to retain the higher active compounds (IC₅₀ <100?nM) in top 1% of total compounds (actives and decoys) as observed after screening the decoy set compounds.

Communicated by Ramaswamy H Sarma     

Identification of some novel pyrazolo[1,5-a]pyrimidine derivatives as InhA inhibitors through pharmacophore-based virtual screening and molecular docking

The InhA inhibitors play key role in mycolic acid synthesis by preventing the fatty acid biosynthesis pathway. In this present article, Pharmacophore modelling and molecular docking study followed by in silico virtual screening could be considered as effective strategy to identify newer enoyl-ACP reductase inhibitors. Pyrrolidine carboxamide derivatives were opted to generate pharmacophore models using HypoGen algorithm in Discovery studio 2.1. Further it was employed to screen Zinc and Minimaybridge databases to identify and design newer potent hit molecules. The retrieved newer hits were further evaluated for their drug likeliness and docked against enoyl acyl carrier protein reductase. Here, novel pyrazolo[1,5-a]pyrimidine analogues were designed and synthesized with good yields. Structural elucidation of synthesized final molecules was perform through IR, MASS, ¹H-NMR, ¹³C-NMR spectroscopy and further tested for its in vitro anti-tubercular activity against H37Rv strain using Microplate Alamar blue assay (MABA) method. Most of the synthesized compounds displayed strong anti-tubercular activities. Further, these potent compounds were gauged for MDR-TB, XDR-TB and cytotoxic study.     

Molecular docking and CoMFA studies of thiazoloquin(az)olin(on)es as CD38 inhibitors: determination of inhibitory mechanism,pharmacophore interactions,and design of new inhibitors

In this research, molecular docking and 3D-QSAR studies were carried out on a series of 79 thiazoloquin(az)olin(on)es as CD38 inhibitors. Based on docking results, four interactions including hydrogen bonding with main chain of GLU-226 (H-M-GLU-226), Van der Waals interactions with side chain of TRP-125 (V-S-TRP-125), TRP-189 (V-S-TRP-189), and THR-221 (V-S-THR-221) were considered as pharmacological interactions. Active conformation of each ligand was extracted from docking studies and was used for carrying out 3D-QSAR modeling. Comparative molecular field analysis (CoMFA) was performed on CD38 inhibitory activities of these compounds on human and mouse. We developed CoMFA models with five components as optimum models for both data-sets. For human data-set, a model with high predictive power was developed. R², RMSE, and F-test values for training set of this model were .94, .24, and 179.58, respectively, and R² and RMSE for its test set were .92 and .32, respectively. The q² and RMSE values for leave-one-out cross validation test on training set were .78 and .46, respectively, that demonstrate created model is robust. Based on extracted steric and electrostatic contour maps for this model, three inhibitors with pIC₅₀ larger than 8.85 were designed.     

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.