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.     

Pharmacophore based virtual screening,molecular docking and biological evaluation to identify novel PDE5 inhibitors with vasodilatory activity

Prompted by the role of PDE5 and its closely associated cAMP and cGMP in hypertension, we have attempted to discover novel PDE5 inhibitors through ligand based virtual screening. Rigorously validated model comprising of one HBA, one HY and one RA was used as a query to search the NCI database leading to retrieval of many compounds which were screened on the basis of estimated activity, fit value and Lipinski’s violation. Selected compounds were subjected to docking studies which resulted into visualization of potential interaction capabilities of NCI compounds in line to pharmacophoric features. Finally three compounds were subjected to in vitro evaluation using the isolated rat aortic model. The results showed that all three compounds are potent and novel PDE5 inhibitors with vasodilatory activity range from 10⁻² to 10⁻⁵ M.     

Pharmacophore modeling,multiple docking,and molecular dynamics studies on Wee1 kinase inhibitors

Wee1-like protein kinase (Wee1) is a tyrosine kinase that regulates the G₂ checkpoint and prevents entry into mitosis in response to DNA damage. Based on a series of signaling pathways initiated by Wee1, Wee1 has been recognized as a potential target for cancer therapy. To discover potent Wee1 inhibitors with novel scaffolds, ligand-based pharmacophore model has been built based on 101 known Wee1 inhibitors. Then the best pharmacophore model, AADRRR.340, with good partial least square (PLS) statistics (R²?=?0.9212, Q²?=?0.7457), was selected and validated. The validated model was used as a three-dimensional (3D) search query for databases virtual screening. The filtered molecules were further analyzed and refined by Lipinski’s rule of 5, multiple docking procedures (high throughput virtual screening (HTVS), standard precision (SP), genetic optimization for ligand docking (GOLD), extra precision (XP), and unique quantum polarized ligand docking (QPLD)); absorption, distribution, metabolism, excretion, and toxicity (ADMET) screening; and the Prime/molecular mechanics generalized born surface area (MM-GBSA) method binding free energy calculations. Eight leads were identified as potential Wee1 inhibitors, and a 50?ns molecular dynamics (MD) simulation was carried out for top four inhibitors to predict the stability of ligand–protein complex. Molecular mechanics Poisson–Boltzmann surface area (MM-PBSA) based on MD simulation and the energy contribution per residue to the binding energy were calculated. In the end, three hits with good stabilization and affinity to protein were identified.

Communicated by Ramaswamy H. Sarma     

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.     

Discovery of novel 5α-reductase type II inhibitors by pharmacophore modelling,virtual screening,molecular docking and molecular dynamics simulations

Benign prostatic hyperplasia (BPH) is caused by augmented levels of androgen dihydrotestosterone (DHT) which is involved in the growth of the prostate in humans. 5α-Reductase type II (5αR2) is an intracellular enzyme that catalyses the formation of DHT from testosterone; hence, the inhibition of 5αR2 has emerged as one of the most promising strategies for the treatment of BPH. In this study, a computational approach that integrates ligand-based pharmacophore modelling, virtual screening, molecular docking and molecular dynamics (MD) simulations was adopted to discover novel 5αR2 inhibitors with less side effects. After validating by Fischer's randomisation and Güner–Henry test, the best quantitative pharmacophore model (Hypo1), consisting of two hydrogen-bond acceptors and three hydrophobic features, was subsequently used as a three-dimensional-query in virtual screening to identify potential hits from Maybridge and National Cancer Institute databases. These hits were further filtered by ADMET (absorption, distribution, metabolism, elimination and toxicology) and molecular docking experiments, and their binding stabilities were validated by 10-ns MD simulations. Finally, only one hit was identified as a potential lead based on higher predicted inhibitory activity to 5αR2 compared with the most active inhibitor (finasteride). Our results further suggest that this potential lead could easily be synthesised and has structural novelty, making it a promising candidate for treating BPH.     

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 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 modeling,molecular docking and molecular dynamics studies on natural products database to discover novel skeleton as non-purine xanthine oxidase inhibitors

Gout is a common inflammatory arthritis caused by the deposition of urate crystals within joints. It is increasingly in prevalence during the past few decades as shown by the epidemiological survey results. Xanthine oxidase (XO) is a key enzyme to transfer hypoxanthine and xanthine to uric acid, whose overproduction leads to gout. Therefore, inhibiting the activity of xanthine oxidase is an important way to reduce the production of urate. In the study, in order to identify the potential natural products targeting XO, pharmacophore modeling was employed to filter databases. Here, two methods, pharmacophore based on ligand and pharmacophore based on receptor-ligand, were constructed by Discovery Studio. Then GOLD was used to refine the potential compounds with higher fitness scores. Finally, molecular docking and dynamics simulations were employed to analyze the interactions between compounds and protein. The best hypothesis was set as a 3D query to screen database, returning 785 and 297 compounds respectively. A merged set of the above 1082 molecules was subjected to molecular docking, which returned 144 hits with high-fitness scores. These molecules were clustered in four main kinds depending on different backbones. What is more, molecular docking showed that the representative compounds established key interactions with the amino acid residues in the protein, and the RMSD and RMSF of molecular dynamics results showed that these compounds can stabilize the protein. The information represented in the study confirmed previous reports. And it may assist to discover and design new backbones as potential XO inhibitors based on natural products.     

Pharmacophore and structure-based drug design,molecular dynamics and admet/tox studies to design novel potential pad4 inhibitors

We have used docking (GLIDE), pharmacophore modeling (Discovery Studio), long trajectory molecular dynamics (Discovery Studio) and ADMET/Tox (QikProp and DEREK) to investigate PAD4 in order to determine potential novel inhibitors and hits. We have carried out virtual screening in the ZINC natural compounds database. Pharmacokinetics and Toxicity of the best hits were assessed using databases implemented in softwares that create models based on chemical structures taking into account consideration about the toxicophoric groups. A wide variety of pharmaceutical relevant properties are determined in order to make decisions about molecular suitability. After screening and analysis, the 6 most promising PAD4 inhibitors are suggested, with strong interactions (pi-stacking, hydrogen bonds, hydrophobic contacts) and suitable pharmacotherapeutic profile as well.     

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.     

Discovery of promising FtsZ inhibitors by E-pharmacophore, 3D-QSAR,molecular docking study,and molecular dynamics simulation

Filamentous temperature-sensitive protein Z (FtsZ), playing a key role in bacterial cell division, is regarded as a promising target for the design of antimicrobial agent. This study is looking for potential high-efficiency FtsZ inhibitors. Ligand-based pharmacophore and E-pharmacophore, virtual screening and molecular docking were used to detect promising FtsZ inhibitors, and molecular dynamics simulation was used to study the stability of protein-ligand complexes in this paper. Sixty-three inhibitors from published literatures with pIC₅₀ ranging from 2.483 to 5.678 were collected to develop ligand-based pharmacophore model. 4DXD bound with 9PC was selected to develop the E-pharmacophore model. The pharmacophore models validated by test set method and decoy set were employed for virtual screening to exclude inactive compounds against ZINC database. After molecular docking, ADME analysis, IFD docking and MM-GBSA, 8 hits were identified as potent FtsZ inhibitors. A 50?ns molecular dynamics simulation was implemented on the compounds to assess the stability between potent inhibitors and FtsZ. The results indicated that the candidate compounds had a high docking score and were strongly combined with FtsZ by forming hydrogen bonding interactions with key amino acid residues, and van der Waals forces and hydrophobic interactions had significant contribution to the stability of the binding. Molecular dynamics simulation results showed that the protein-ligand compounds performed well in both the stability and flexibility of the simulation process.     

Pharmacophore identification of PAK4 inhibitors

p21-activated kinase 4 (PAK4) is a serine–threonine protein kinase which plays an important role in a wide variety of human diseases including cancer. The inhibition of this kinase is of great interest for the treatment of cancer. In the present study, we report three pharmacophore models of PAK4 based on a small set of PAK4 inhibitors, a crystal structure of PAK4 and a docked complex between PAK4 and a potent inhibitor. These results might provide useful and reliable tools in identifying structurally diverse compounds with desired biological activity.     

Pharmacophore modeling and molecular dynamics approach to identify putative DNA Gyrase B inhibitors for resistant tuberculosis

One of the major mechanisms followed by the therapeutic agents to target the causative organism of TB, mycobacterium tuberculosis (Mtb), involves disruption of the replication cycle of the pathogen DNA. The process involves two steps that occur simultaneously, ie, breakage and reunion of DNA at gyrase A (GyrA) domain and ATP hydrolysis at gyrase B (GyrB) domain. Current therapy for multi-drug resistant TB involves FDA approved, Fluoroquinolone-based antibiotics, which act by targeting the replication process at GyrA domain. However, resistance against fluoroquinolones due to mutations in the GyrA domain has limited the use of this therapy and shifted the focus of the research community on the GyrB domain. Thus, this study involves in silico designing of chemotherapeutic agents for resistant TB by targeting GyrB domain. In the current study, a pharmacophore model for GyrB domain was generated using reported inhibitors. It was utilized as a query search against three commercial databases to identify GyrB domain inhibitors. Additionally, a qualitative Hip-Hop pharmacophore model for GyrA was also developed on the basis of some marketed fluoroquinolone-based GyrA inhibitors, to remove non-selective gyrase inhibitors obtained in virtual screening. Further, molecular dynamic simulations were carried out to determine the stability of the obtained molecules in complex with both the domains. Finally, Molecular mechanics with generalized Born and surface area solvation score was calculated to determine the binding affinity of obtained molecule with both domains to determine the selectivity of the obtained molecules that resulted in seven putative specific inhibitors of GyrB domain.     

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.     

Development of hydroxylated naphthylchalcones as polyphenol oxidase inhibitors: Synthesis,biochemistry and molecular docking studies

Polyphenol oxidase (Tyrosinase) has received great attention, since it is the key enzyme in melanin biosynthesis. In this study, novel hydroxy naphthylchalcone compounds were synthesized, and their inhibitory effects on mushroom tyrosinase activity were evaluated. The structures of the compounds synthesized were confirmed by ¹H NMR, ¹³C NMR, FTIR and HRMS. Two of the compounds synthesized inhibited the diphenolase activity of tyrosinase in a dose dependent manner and exhibited much higher tyrosinase inhibitory activities (IC₅₀ values of 10.4 μM and 14.4 μM, respectively) than the positive control, kojic acid (IC₅₀: 27.5 μM). Kinetic analysis showed that their inhibition was reversible. Both the novel compounds displayed competitive inhibition with their K_i values of 3.8 μM and 4.5 μM, respectively. Docking results confirmed that the active inhibitors strongly interacted with the mushroom tyrosinase residues. This study suggests hydroxy naphthylchalcone compounds to serve as promising candidates for use as depigmentation agents.     

Identification of potential tumour-associated carbonic anhydrase isozyme IX inhibitors: atom-based 3D-QSAR modelling,pharmacophore-based virtual screening and molecular docking studies

Abstract

Tumour hypoxia results in dramatic changes in the gene expression, proliferation and survival of tumour cells. The tumour cells shift towards anaerobic glycolysis which results in change of pH in their microenvironment. In response to this stress, over expression of carbonic anhydrase IX (CA IX) genes is observed in many solid tumours. So, selective inhibition of CA IX can be a promising target for anti-cancer drugs. In this work in silico tools like atom-based 3D-QSAR modelling, pharmacophore-based virtual screening and molecular docking were used to identify potential CA IX inhibitors. Based on the training set used in the QSAR model, twenty pharmacophore models were generated. Out of these, HHHR_1, AHHR_1, DHHHR_1, AHHHR_1 model was used to screen a database of 1,50,000 compounds retrieved from ZINC 15 database. R² and Q² was 0.9864 and 0.8799, respectively, for the developed QSAR model. 163 compounds showed a phase screen score above 2.4 in which ZINC02260669 was the highest ranked (screen score, 2.852058) compound in all the four models. Built QSAR model was used to predict the activity of all these 163 compounds and ZINC72370966 showed the highest predicted activity with pKi value of 7.649. These compounds were docked against CA IX (human) protein (PDB ID 5FL6) and molecular docking results showed favourable binding interactions for the best ten identified hits. This work gives design insights and some potential scaffolds which can be developed as CA IX inhibitors.

Communicated by Ramaswamy H. Sarma     

Pharmacophore feature prediction and molecular docking approach to identify novel anti‐HCV protease inhibitors

Discovering a potential drug for HCV treatment is a challenging task in the field of drug research. This study initiates with computational screening and modeling of promising ligand molecules. The foremost modeling method involves the identification of novel compound and its molecular interaction based on pharmacophore features. A total of 197 HCV compounds for NS3/4A protein target were screened for our study. The pharmacophore models were generated using PHASE module implemented in Schrodinger suite. The pharmacophore features include one hydrogen bond acceptor, one hydrogen bond donor, and three hydrophobic sites. As a result, based on mentioned hypothesis the model ADHHH.159 corresponds to the CID 59533233. Furthermore, docking was performed using maestro for all the 197 compounds. Among these, the CID 59533313 and 59533233 possess the best binding energy of ?11.75 and ?10.40 kcal/mol, respectively. The interactions studies indicated that the CID complexed with the NS3/4A protein possess better binding affinity with the other compounds. Further the compounds were subjected to calculate the ADME properties. Therefore, it can be concluded that these two compounds could be a potential alternative drug for the development of HCV.     

Multicomplex-based pharmacophore modeling in conjunction with multi-target docking and molecular dynamics simulations for the identification of PfDHFR inhibitors

Abstract

Plasmodium falciparum dihydrofolate reductase enzyme (PfDHFR) is counted as one of the attractive and validated antimalarial drug targets. However, the point mutations in the active site of wild-type PfDHFR have developed resistance against the well-known antifolates. Therefore, there is a dire need for the development of inhibitors that can inhibit both wild-type and mutant-type DHFR enzyme. In the present contribution, we have constructed the common feature pharmacophore models from the available PfDHFR. A representative hypothesis was prioritized and then employed for the screening of natural product library to search for the molecules with complementary features responsible for the inhibition. The screened candidates were processed via drug-likeness filters and molecular docking studies. The docking was carried out on the wild-type PfDHFR (3QGT); double-mutant PfDHFR (3UM5 and 1J3J) and quadruple-mutant PfDHFR (1J3K) enzymes. A total of eight common hits were obtained from the docking calculations that could be the potential inhibitors for both wild and mutant type DHFR enzymes. Eventually, the stability of these candidates with the selected proteins was evaluated via molecular dynamics simulations. Except for SPECS14, all the prioritized candidates were found to be stable throughout the simulation run. Overall, the strategy employed in the present work resulted in the retrieval of seven candidates that may show inhibitory activity against PfDHFR and could be further exploited as a scaffold to develop novel antimalarials.

Communicated by Ramaswamy H. Sarma     

Biological evaluation and docking studies of natural isocoumarins as inhibitors for human kallikrein 5 and 7

Human kallikrein 5 and 7 (KLK5 and KLK7) are trypsin-like and chymotrypsin-like serine proteases, respectively, and promising targets for the treatment of skin desquamation, inflammation and cancer. In an effort to develop new inhibitors for these enzymes, we carried out enzymatic inhibition assays and docking studies with three isocoumarin compounds. Some promising inhibitors were uncovered, with vioxanthin and 8,8'-paepalantine being the most potent competitive inhibitors of KLK5 (K(i)=22.9 μM) and KLK7 (K(i)=12.2 μM), respectively. Our docking studies showed a good correlation with the experimental results, and revealed a distinct binding mode for the inhibitors at the binding sites of KLK5 and KLK7. In addition, the docking results suggested that the formation of hydrogen bonds at the oxyanion hole is essential for a good inhibitor.