Comparative modeling of UDP-N-acetylmuramoyl-glycyl-D-glutamate-2, 6-diaminopimelate ligase from Mycobacterium leprae and analysis of its binding features through molecular docking studies

Leprosy is an infectious disease caused by Mycobacterium leprae. The increasing drug and multi-drug resistance of M. leprae enforce the importance of finding new drug targets. Mycobacterium has unusually impermeable cell wall that contributes to considerable resistance to many drugs. Peptidoglycan is an important component of the cell wall of M. leprae. UDP-N-acetylmuramoyl-glycyl-D-glutamate-2, 6-diaminopimelate ligase (MurE) plays a crucial role in the peptidoglycan biosynthesis and hence it could be considered as a potential drug target for leprosy. Structure of this enzyme for M. leprae has not yet been elucidated. We modeled the three-dimensional structure of MurE from M. leprae using comparative modeling methods based on the X-ray crystal structure of MurE from E. coli and validated. The 3D-structure of M. leprae MurE enzyme was docked with its substrates meso-diaminopimelic acid (A2pm) and UDP-N-acetyl muramoyl-glycyl-D- glutamate (UMGG) and its product UDP-N-acetyl muramoyl-glycyl-D-glu-meso-A₂pm (UTP) and also with ATP. The docked complexes reveal the amino acids responsible for binding the substrates. Superposition of these complex structures suggests that carboxylic acid group of UMGG is positioned in proximity to γ-phosphate of the ATP to facilitate the formation of acylphosphate intermediate. The orientation of an amino group of A₂pm facilitates the nucleophilic attack to form the product. Overall, the proposed model together with its binding features gained from docking studies could help to design a truly selective ligand inhibitor specific to MurE for the treatment of leprosy.     

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.     

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.     

Characterization of binding site of closed-state KCNQ1 potassium channel by homology modeling, molecular docking, and pharmacophore identification

This investigation was performed to assess the importance of interaction in the binding of blockers to KCNQ1 potassium using molecular modeling. This work could be considered made up by three main steps: (1) the construction of closed-state structure of KCNQ1 through homology modeling; (2) the automated docking of three blockers: IKS-142, L-735821, and BMS-IKS, using DOCK program; (3) the generation and validation of pharmacophore for KCNQ1 ligands using Catalyst/HypoGen. The obtained results highlight the hydrophobic or aromatic residues involved in S6 transmembrane domain and the base of the pore helix of KCNQ1, confirming the mutagenesis data and pharmacophore model, and giving new suggestions for the rational design of novel KCNQ1 ligands.     

Discovery of FIXa inhibitors by combination of pharmacophore modeling,molecular docking,and 3D-QSAR modeling

Human Coagulation Factor IXa (FIXa), specifically inhibited at the initiation stage of the blood coagulation cascade, is an excellent target for developing selective and safe anticoagulants. To explore this inhibitory mechanism, 86 FIXa inhibitors were selected to generate pharmacophore models and subsequently SAR models. Both best pharmacophore model and ROC curve were built through the Receptor–Ligand Pharmacophore Generation module. CoMFA model based on molecular docking and PLS factor analysis methods were developed. Model propagations values are q²?=?0.709, r²?=?0.949, and r²_pred?=?0.905. The satisfactory q² value of 0.609, r² value of 0.962, and r²_pred value of 0.819 for CoMSIA indicated that the CoMFA and CoMSIA models are both available to predict the inhibitory activity on FIXa. On the basis of pharmacophore modeling, molecular docking, and 3D-QSAR modeling screening, six molecules are screened as potential FIXa inhibitors.     

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.     

Combined pharmacophore modeling, 3D-QSAR and docking studies to identify novel HDAC inhibitors using drug repurposing

Abstract

Histone deacetylases (HDACs), a critical family of epigenetic enzymes, has emerged as a promising target for antitumor drugs. Here, we describe our protocol of virtual screening in identification of novel potential HDAC inhibitors through pharmacophore modeling, 3D-QSAR, molecular docking and molecular dynamics (MD) simulation. Considering the limitation of current virtual screening works, drug repurposing strategy was applied to discover druggable HDAC inhibitor. The ligand-based pharmacophore and 3D-QSAR models were established, and their reliability was validated by different methods. Then, the DrugBank database was screened, followed by molecular docking. MD simulation (100?ns) was performed to further study the stability of ligand binding modes. Finally, results indicated the hit DB03889 with high in silico inhibitory potency was suitable for further experimental analysis.

Communicated by Ramaswamy H. Sarma     

Discovery of novel inhibitors of 11beta-hydroxysteroid dehydrogenase type 1 by docking and pharmacophore modeling

11beta-hydroxysteroid dehydrogenase type 1 (11beta-HSD1) is a potential target for treatment of diabetes and metabolic syndrome. Docking and pharmacophore modeling have been used to discover novel inhibitors of 11beta-HSD1. Several compounds, with large structural diversity and good potency against 11beta-HSD1, have been found and their potency was determined by the enzyme assay. New scaffolds of 11beta-HSD1 inhibitors are also reported.     

Rational design,synthesis, pharmacophore modeling,and docking studies for identification of novel potent DNA-PK inhibitors

Drugs of cancer based upon ionizing radiation or chemotherapeutic treatment may affect breaking of DNA double strand in cell. DNA-PK enzyme has emerged as an attractive target for drug discovery efforts toward DNA repair pathways. Hence, the search for potent and selective DNA-PK inhibitors has particularly considered state-of-the art and several series of inhibitors have been designed. In this article, a novel benchmark DNA-PK database of 43 compounds was built and described. Ligand-based approaches including pharmacophore and QSAR modeling were applied and novel models were introduced and analyzed for predicting activity test for DNA-PK drug candidates. Based upon the modeling results, we gave a report of synthesis of fifteen novel 2-((8-methyl-2-morpholino-4-oxo-4H-benzo[e][1,3]oxazin-7-yl)oxy)acetamide derivatives and in vitro evaluation for DNA-PK inhibitory and antiproliferative activities. These fifteen compounds overall are satisfied with Lipinski's rule of five. The biological testing of target compounds showed five promising active compounds 7c, 7d, 7f, 9e and 9f with micromolar DNA-PK activity range from 0.25 to 5 µM. In addition, SAR of the compounds activity was investigated and confirmed that the terminal aryl moiety was found to be quite crucial for DNA-PK activity. Moreover flexible docking simulation was done for the potent compounds into the putative binding site of the 3D homology model of DNA-PK enzyme and the probable interaction model between DNA-PK and the ligands was investigated and interpreted.     

Computational screening of dipeptidyl peptidase IV inhibitors from micoroalgal metabolites by pharmacophore modeling and molecular docking

Dipeptidyl peptidase IV (DPP‐IV) catalyzes conversion of GLP‐1 (glucagon like peptide 1) to inert structure, which results in insufficient secretion of insulin and increase in postprandial blood glucose level. The present study attempts to identify novel inhibitors from bioactive metabolites present in microalgae against DPP‐IV through virtual screening, molecular docking, and pharmacophore modeling for the active target. Possible binding modes of all 60 ligands against DPP‐IV receptor were constructed using MTiOpenScreen virtual screening server. Pharmacophore model was built based on identified 38 DPP‐IV test ligands by using the web‐based PharmaGist program which encompasses hydrogen‐bond acceptors, hydrophobic groups, spatial features, and aromatic rings. The pharmacophore model having highest scores was selected to screen active DPP‐IV ligands. Highest scoring model was used as a query in ZincPharmer screening. All identified ligands were filtered, based on the Lipinski's rule‐of‐five and were subjected to docking studies. In the process of docking analyses, we considered different bonding modes of one ligand with multiple active cavities of DPP‐IV with the help of AutoDock 4.0. The docking analyses indicate that the bioactive constituents, namely, β‐stigmasterol, barbamide, docosahexaenoic acid, arachidonic acid, and harman showed the best binding energies on DPP‐IV receptor and hydrogen bonding with ASP545, GLY741, TYR754, TYR666, ARG125, TYR547, SER630, and LYS554 residues. This study concludes that docosahexaenoic acid, arachidonic acid, β‐stigmasterol, barbamide, harman, ZINC58564986, ZINC56907325, ZINC69432950, ZINC69431828, ZINC73533041, ZINC84287073, ZINC69849395, and ZINC10508406 act as possible DPP‐IV inhibitors.     

The discovery of potential acetylcholinesterase inhibitors: A combination of pharmacophore modeling,virtual screening,and molecular docking studies

Background  

Alzheimer's disease (AD) is the most common cause of dementia characterized by progressive cognitive impairment in the elderly people. The most dramatic abnormalities are those of the cholinergic system. Acetylcholinesterase (AChE) plays a key role in the regulation of the cholinergic system, and hence, inhibition of AChE has emerged as one of the most promising strategies for the treatment of AD.     

Homology modeling of human CCR5 and analysis of its binding properties through molecular docking and molecular dynamics simulation

In this study, homology modeling, molecular docking and molecular dynamics simulation were performed to explore structural features and binding mechanism of some inhibitors of chemokine receptor type 5 (CCR5), and to construct a model for designing new CCR5 inhibitors for preventing HIV attachment to the host cell. A homology modeling procedure was employed to construct a 3D model of CCR5. For this procedure, the X-ray crystal structure of bovine rhodopsin (1F88A) at 2.80? resolution was used as template. After inserting the constructed model into a hydrated lipid bilayer, a 20ns molecular dynamics (MD) simulation was performed on the whole system. After reaching the equilibrium, twenty-four CCR5 inhibitors were docked in the active site of the obtained model. The binding models of the investigated antagonists indicate the mechanism of binding of the studied compounds to the CCR5 obviously. Moreover, 3D pictures of inhibitor-protein complex provided precious data regarding the binding orientation of each antagonist into the active site of this protein. One additional 20 ns MD simulation was performed on the initial structure of the CCR5-ligand 21 complex, resulted from the previous docking calculations, embedded in a hydrated POPE bilayer to explore the effects of the presence of lipid bilayer in the vicinity of CCR5-ligand complex. This article is part of a Special Issue entitled Protein translocation across or insertion into membranes.     

Homology modeling and docking studies of human Bcl-2L10 protein

Cancer, an unrestrained proliferation of cells, is one of the lead cause of death. Nearly 12.5 million people are diagnosed with cancer worldwide, 7.5 million people die of which 2.5 million cases are from India. Major cause for cancer is restriction of programmed cell death (apoptosis). Multiple signaling pathways regulate apoptosis. Bcl-2 (B - Cell Lymphomas-2) family proteins play a vital role as central regulators of apoptosis. Bcl-2L10, a novel anti-apoptotic protein, blocks apoptosis by mitochondrial dependent mechanism. The present study evaluates the 3D structure of Bcl-2L10 protein using homology modeling and aims to understand plausible functional and binding interactions between Bcl-2L10 with BH3 domain of BAX using protein - protein docking. The docking studies show binding of BH3 domain at Lys 110, Trp-111, Pro-115, Glu-119 and Asp-127 in the groove of BH 1, 2 and 3 domains of Bcl-2L10. Heterodimerization of anti-apoptotic Bcl-2 and BH3 domain of pro-apoptotic Bcl-2 proteins instigates apoptosis. Profound understanding of Bcl-2 pathway may prove useful in identification of future therapeutic targets for cancer.     

Fluorescence resonance energy transfer analysis of secretin docking to its receptor: mapping distances between residues distributed throughout the ligand pharmacophore and distinct receptor residues

Full structural characterization of G protein-coupled receptors has been limited to rhodopsin, with its uniquely stable structure and ability to be crystallized. For other members of this important superfamily, direct structural insights have been limited to NMR structures of soluble domains. Two members of the Class II family have recently had the structures of their isolated amino-terminal regions solved by NMR, yet it remains unclear how that domain is aligned with the heptahelical transmembrane bundle domain of those receptors. Indeed, three distinct orientations have been suggested for different members of this family. In the current work, we have utilized fluorescence resonance energy transfer to establish the distances between four residues distributed throughout fully biologically active, high affinity analogues of secretin and distinct residues in each of four extracellular regions of the intact secretin receptor. These 16 distance constraints were utilized along with nine photoaffinity labeling spatial approximation constraints to study the three proposed orientations of the peptide-binding amino terminus and helical bundle domains of this receptor. In the best model, the carboxyl terminus of secretin was found to bind in a groove above the beta-hairpin region of the receptor amino terminus, with its amino-terminal end adjacent to the third extracellular loop and top of transmembrane segment VI. This refined model of the intact receptor was also fully consistent with the spatial approximation of the Trp(48)-Asp(49)-Asn(50) endogenous agonist segment with the third extracellular loop region that it has been shown to photolabel. This provides strong evidence for the orientation of peptide-binding and signaling domains of a prototypic Class II G protein-coupled receptor.     

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.     

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.     

New insights about HERG blockade obtained from protein modeling, potential energy mapping, and docking studies

We created a homology model of the homo-tetrameric pore domain of HERG using the crystal structure of the bacterial potassium channel, KvAP, as a template. We docked a set of known blockers with well-characterized effects on channel function into the lumen of the pore between the selectivity filter and extracellular entrance using a novel docking and refinement procedure incorporating Glide and Prime. Key aromatic groups of the blockers are predicted to form multiple simultaneous ring stacking and hydrophobic interactions among the eight aromatic residues lining the pore. Furthermore, each blocker can achieve these interactions via multiple docking configurations. To further interpret the docking results, we mapped hydrophobic and hydrophilic potentials within the lumen of each refined docked complex. Hydrophilic iso-potential contours define a 'propeller-shaped' volume at the selectivity filter entrance. Hydrophobic contours define a hollow 'crown-shaped' volume located above the 'propeller', whose hydrophobic 'rim' extends along the pore axis between Tyr652 and Phe656. Blockers adopt conformations/binding orientations that closely mimic the shapes and properties of these contours. Blocker basic groups are localized in the hydrophilic 'propeller', forming electrostatic interactions with Ser624 rather than a generally accepted pi-cation interaction with Tyr652. Terfenadine, cisapride, sertindole, ibutilide, and clofilium adopt similar docked poses, in which their N-substituents bridge radially across the hollow interior of the 'crown' (analogous to the hub and spokes of a wheel), and project aromatic/hydrophobic portions into the hydrophobic 'rim'. MK-499 docks with its longitudinal axis parallel to the axis of the pore and 'crown', and its hydrophobic groups buried within the hydrophobic 'rim'.