Toll-Like Receptor 7 Agonists: Chemical Feature Based Pharmacophore Identification and Molecular Docking Studies

Chemical feature based pharmacophore models were generated for Toll-like receptors 7 (TLR7) agonists using HypoGen algorithm, which is implemented in the Discovery Studio software. Several methods tools used in validation of pharmacophore model were presented. The first hypothesis Hypo1 was considered to be the best pharmacophore model, which consists of four features: one hydrogen bond acceptor, one hydrogen bond donor, and two hydrophobic features. In addition, homology modeling and molecular docking studies were employed to probe the intermolecular interactions between TLR7 and its agonists. The results further confirmed the reliability of the pharmacophore model. The obtained pharmacophore model (Hypo1) was then employed as a query to screen the Traditional Chinese Medicine Database (TCMD) for other potential lead compounds. One hit was identified as a potent TLR7 agonist, which has antiviral activity against hepatitis virus in vitro. Therefore, our current work provides confidence for the utility of the selected chemical feature based pharmacophore model to design novel TLR7 agonists with desired biological activity.     

Pharmacophore modeling study based on known Spleen tyrosine kinase inhibitors together with virtual screening for identifying novel inhibitors

In this investigation, chemical features based 3D pharmacophore models were developed based on the known inhibitors of Spleen tyrosine kinase (Syk) with the aid of hiphop and hyporefine modules within catalyst. The best quantitative pharmacophore model, Hypo1, was used as a 3D structural query for retrieving potential inhibitors from chemical databases including Specs, NCI, MayBridge, and Chinese Nature Product Database (CNPD). The hit compounds were subsequently subjected to filtering by Lipinski’s rule of five and docking studies to refine the retrieved hits. Finally 30 compounds were selected from the top ranked hit compounds and conducted an in vitro kinase inhibitory assay. Six compounds showed a good inhibitory potency against Syk, which have been selected for further investigation.     

Pharmacophore modeling and virtual screening for designing potential PLK1 inhibitors

Pharmacophore models of Polo-like kinase-1 (PLK1) inhibitors have been established by using the HipHop and HypoGen algorithms implemented in the Catalyst software package. The best quantitative pharmacophore model, Hypo1, which has the highest correlation coefficient (0.9895), consists of one hydrogen bond acceptor, one hydrogen bond donor, one hydrophobic feature, and one hydrophobic aliphatic feature. Hypo1 was further validated by test set and cross validation method. Then Hypo1 was used as a 3D query to screen several databases including Specs, NCI, Maybridge, and Chinese Nature Product Database (CNPD). The hit compounds were subsequently subjected to filtering by Lipinski's rule of five and docking study to refine the retrieved hits and as a result to reduce the rate of false positive. Finally, a total of 20 compounds were selected and have been shifted to in vitro and in vivo studies. As far as we know, this is the first report on the pharmacophore modeling even the first publicly reported virtual screening study of PLK1 inhibitors.     

Pharmacophore modeling and hybrid virtual screening for the discovery of novel IκB kinase 2 (IKK2) inhibitors

IKK2 (IκB kinase 2) inhibitors have been identified as potential drug candidates in the treatment of various immune/inflammatory disorders as well as cancer. So far more than one hundred small molecule inhibitors against IKK2 have been reported publicly. In this investigation, pharmacophore modeling was carried out to clarify the essential structure-activity relationship for the known IKK2 inhibitors. One of the established pharmacophore hypotheses, namely Hypo8, which has the best prediction ability to an external test data set, was suggested as a template for virtual screening. Evaluation of the performances of Hypo8 and a hybrid method (Hypo81docking) in virtual screening indicated that the use of the hybrid virtual screening considerably increased the hit rate and enrichment factor. The hybrid method was therefore adopted for screening several commercially available chemical databases, including Specs, NCI, Maybridge and Chinese Nature Product Database (CNPD), for novel potent IKK2 inhibitors. The hit compounds were subsequently subjected to filtering by Lipinski's rule of five. Finally some of the final hit compounds were selected and suggested for further experimental investigations.     

Pharmacophore modeling and in silico screening for new KDR kinase inhibitors

In order to elucidate the essential structural features for KDR kinase inhibitors, three-dimensional pharmacophore hypotheses were built on the basis of a set of known KDR kinase inhibitors selected from the literature with CATALYST program. Several methods tools used in validation of pharmacophore hypothsis were presented, and the first hypothesis (Hypo1) was considered to be the best pharmacophore hypothesis. The model (Hypo1) was then employed as 3D search query to screen the Traditional Chinese Medicine Database (TCMD) for other potential lead compounds. One hit illustrated high binding affinity with KDR kinase measured by the surface plasmon resonance biosensor. Docking studies may help elucidate the mechanisms of KDR kinase receptor-ligand interactions.     

Ligand-based and structure-based approaches in identifying ideal pharmacophore against c-Jun N-terminal kinase-3

Structure and ligand based pharmacophore modeling and docking studies carried out using diversified set of c-Jun N-terminal kinase-3 (JNK3) inhibitors are presented in this paper. Ligand based pharmacophore model (LBPM) was developed for 106 inhibitors of JNK3 using a training set of 21 compounds to reveal structural and chemical features necessary for these molecules to inhibit JNK3. Hypo1 consisted of two hydrogen bond acceptors (HBA), one hydrogen bond donor (HBD), and a hydrophobic (HY) feature with a correlation coefficient (r²) of 0.950. This pharmacophore model was validated using test set containing 85 inhibitors and had a good r² of 0.846. All the molecules were docked using Glide software and interestingly, all the docked conformations showed hydrogen bond interactions with important hinge region amino acids (Gln155 and Met149) and these interactions were compared with Hypo1 features. The results of ligand based pharmacophore model (LBPM) and docking studies are validated each other. The structure based pharmacophore model (SBPM) studies have identified additional features, two hydrogen bond donors and one hydrogen bond acceptor. The combination of these methodologies is useful in designing ideal pharmacophore which provides a powerful tool for the discovery of novel and selective JNK3 inhibitors.     

Pharmacophore modeling and virtual screening for designing potential 5-Lipoxygenase inhibitors

Inhibitors of the 5-Lipoxygenase (5-LOX) pathway have a therapeutic potential in a variety of inflammatory disorders such as asthma. In this study, chemical feature based pharmacophore models of inhibitors of 5-LOX have been developed with the aid of HipHop and HypoGen modules within Catalyst program package. The best quantitative pharmacophore model, Hypo1, which has the highest correlation coefficient (0.97), consists of two hydrogen-bond acceptors, one hydrophobic feature and one ring aromatic feature. Hypo1 was further validated by test set and cross validation method. The application of the model shows great success in predicting the activities of 65 known 5-LOX inhibitors in our test set with a correlation coefficient of 0.85 with a cross validation of 95% confidence level, proving that the model is reliable in identifying structurally diverse compounds for inhibitory activity against 5-LOX. Furthermore, Hypo1 was used as a 3D query for screening Maybridge and NCI databases within catalyst and also drug like compounds obtained from Enamine Ltd, which follow Lipinski’s rule of five. The hit compounds were subsequently subjected to filtering by docking and visualization, to identify the potential lead molecules. Finally 5 potential lead compounds, identified in the above process, were evaluated for their inhibitory activities. These studies resulted in the identification of two compounds with potent inhibition of 5-LOX activity with IC₅₀ of 14 μM and 35 μM, respectively. These studies thus validate the pharmacophore model generated and suggest the usefulness of the model in screening of various small molecule libraries and identification of potential lead compounds for 5-LOX inhibition.     

3D QSAR pharmacophore-based virtual screening for the identification of potential inhibitors of tyrosinase

Abstract

Tyrosinase plays an important role in melanin biosynthesis and protects skin against ultraviolet radiations. Functional deficiency of tyrosinase results in serious dermatological diseases. Tyrosinase also participates in neuromelanin formation in the human brain, which leads to neurodegeneration resulting in Parkinson’s disease. In fruits and vegetables, tyrosinase plays a critical role in senescence, causing undesired browning that results in faster deterioration and shorter shelf lines. The only commercially available tyrosinase is mushroom tyrosinase and it shows the highest homology to the mammalian tyrosinase. Although kojic acid is currently used as a tyrosinase inhibitor, they have serious side effects such as dermatitis, carcinogenesis and hepatotoxicity. Therefore, in order to develop a more active and safer tyrosinase inhibitor, 3D QSAR pharmacophore models were generated based on experimentally known inhibitors. The pharmacophore model, Hypo1, was developed with a large cost difference, high correlation coefficient and low RMS deviation. Hypo1 showed a good spatial arrangement; consisting of five-point features including two hydrogen bond acceptor, one hydrogen bond donor and two hydrophobic features. Hypo1 was further validated by cost analysis, test set and Fisher’s randomisation method. Hypo1 was used as a 3D query for screening the in-house drug-like databases, and the hits were further selected by applying ADMET, Lipinski’s rule of five and fit value criteria. To identify binding conformations, the obtained hits were subjected to molecular docking. Finally, molecular dynamics simulations revealed the appropriate binding modes of hit compounds. To conclude, we propose the final three hit compounds with new structural scaffolds as a virtual candidate as tyrosinase inhibitors.

Communicated by Ramaswamy H. Sarma     

Pharmacophore Modeling and Virtual Screening for the Discovery of New type 4 cAMP Phosphodiesterase (PDE4) Inhibitors

Type 4 cAMP phosphodiesterase (PDE4) inhibitors show a broad spectrum of anti-inflammatory effects in almost all kinds of inflamed cells, by an increase in cAMP levels which is a pivotal second messenger responsible for various biological processes. These inhibitors are now considered as the potential drugs for treatment of chronic inflammatory diseases. However, some recently marketed inhibitors e.g., roflumilast, have shown adverse effects such as nausea and emesis, thus restricting its use. In order to identify novel PDE4 inhibitors with improved therapeutic indexes, a highly correlating (r = 0.963930) pharmacophore model (Hypo1) was established on the basis of known PDE4 inhibitors. Validated Hypo1 was used in database screening to identify chemical with required pharmacophoric features. These compounds are further screened by using the rule of five, ADMET and molecular docking. Finally, twelve hits which showed good results with respect to following properties such as estimated activity, calculated drug-like properties and scores were proposed as potential leads to inhibit the PDE4 activity. Therefore, this study will not only assist in the development of new potent hits for PDE4 inhibitors, but also give a better understanding of their interaction with PDE4. On a wider scope, this will be helpful for the rational design of novel potent enzyme inhibitors.     

Novel virtual lead identification in the discovery of hematopoietic cell kinase (HCK) inhibitors: application of 3D QSAR and molecular dynamics simulation

High level of hematopoietic cell kinase (Hck) is associated with drug resistance in chronic myeloid leukemia. Additionally, Hck activity has also been connected with the pathogenesis of HIV-1 and chronic obstructive pulmonary disease. In this study, three-dimensional (3D) QSAR pharmacophore models were generated for Hck based on experimentally known inhibitors. A best pharmacophore model, Hypo1, was developed with high correlation coefficient (0.975), Low RMS deviation (0.60) and large cost difference (49.31), containing three ring aromatic and one hydrophobic aliphatic feature. It was further validated by the test set (r?=?0.96) and Fisher’s randomization method (95%). Hypo 1 was used as a 3D query for screening the chemical databases, and the hits were further screened by applying Lipinski’s rule of five and ADMET properties. Selected hit compounds were subjected to molecular docking to identify binding conformations in the active site. Finally, the appropriate binding modes of final hit compounds were revealed by molecular dynamics (MD) simulations and free energy calculation studies. Hence, we propose the final three hit compounds as virtual candidates for Hck inhibitors.     

Design checkpoint kinase 2 inhibitors by pharmacophore modeling and virtual screening techniques

Damage to DNA is caused by ionizing radiation, genotoxic chemicals or collapsed replication forks. When DNA is damaged or cells fail to respond, a mutation that is associated with breast or ovarian cancer may occur. Mammalian cells control and stabilize the genome using a cell cycle checkpoint to prevent damage to DNA or to repair damaged DNA. Checkpoint kinase 2 (Chk2) is one of the important kinases, which strongly affects DNA-damage and plays an important role in the response to the breakage of DNA double-strands and related lesions. Therefore, this study concerns Chk2. Its purpose is to find potential inhibitors using the pharmacophore hypotheses (PhModels) and virtual screening techniques. PhModels can identify inhibitors with high biological activities and virtual screening techniques are used to screen the database of the National Cancer Institute (NCI) to retrieve compounds that exhibit all of the pharmacophoric features of potential inhibitors with high interaction energy. Ten PhModels were generated using the HypoGen best algorithm. The established PhModel, Hypo01, was evaluated by performing a cost function analysis of its correlation coefficient (r), root mean square deviation (RMSD), cost difference, and configuration cost, with the values 0.955, 1.28, 192.51, and 16.07, respectively. The result of Fischer’s cross-validation test for the Hypo01 model yielded a 95% confidence level, and the correlation coefficient of the testing set (r_test) had a best value of 0.81. The potential inhibitors were then chosen from the NCI database by Hypo01 model screening and molecular docking using the cdocker docking program. Finally, the selected compounds exhibited the identified pharmacophoric features and had a high interaction energy between the ligand and the receptor. Eighty-three potential inhibitors for Chk2 are retrieved for further study.     

Identification of potent virtual leads to design novel PLK1 inhibitors: pharmacophore modelling,virtual screening and molecular docking studies

The aim of this study was to identify novel scaffolds and utilise them in designing potent PLK1 inhibitors. Three-dimensional pharmacophore models on the basis of chemical features were developed for PLK1 on the basis of the known inhibitors. The best pharmacophore model, Hypo 1, which has the highest correlation (0.96), the highest cost difference (75.7494), the lowest total cost and RMSD (75.7494, 0.5458), contains two hydrophobics, one ring aromatic and one hydrogen donor. Hypo 1 was validated by the test set, decoy set and the Fischer's randomisation method. Then it was used for chemical database virtual screening. The hit compounds were filtered by Lipinski's rule of five and absorption, distribution, metabolism, elimination and toxicity properties. Finally, 24 compounds with good estimated activity values were used for docking studies. These results will be used to develop new inhibitors of PLK1 as leads.     

Pharmacophore modeling, molecular docking, and molecular dynamics simulation approaches for identifying new lead compounds for inhibiting aldose reductase 2

Aldose reductase 2 (ALR2), which catalyzes the reduction of glucose to sorbitol using NADP as a cofactor, has been implicated in the etiology of secondary complications of diabetes. A pharmacophore model, Hypo1, was built based on 26 compounds with known ALR2-inhibiting activity values. Hypo1 contains important chemical features required for an ALR2 inhibitor, and demonstrates good predictive ability by having a high correlation coefficient (0.95) as well as the highest cost difference (128.44) and the lowest RMS deviation (1.02) among the ten pharmacophore models examined. Hypo1 was further validated by Fisher's randomization method (95%), test set (r = 0.91), and the decoy set shows the goodness of fit (0.70). Furthermore, during virtual screening, Hypo1 was used as a 3D query to screen the NCI database, and the hit leads were sorted by applying Lipinski's rule of five and ADME properties. The best-fitting leads were subjected to docking to identify a suitable orientation at the ALR2 active site. The molecule that showed the strongest interactions with the critical amino acids was used in molecular dynamics simulations to calculate its binding affinity to the candidate molecules. Thus, Hypo1 describes the key structure-activity relationship along with the estimated activities of ALR2 inhibitors. The hit molecules were searched against PubChem to find similar molecules with new scaffolds. Finally, four molecules were found to satisfy all of the chemical features and the geometric constraints of Hypo1, as well as to show good dock scores, PLPs and PMFs. Thus, we believe that Hypo1 facilitates the selection of novel scaffolds for ALR2, allowing new classes of ALR2 inhibitors to be designed.     

Identification of small-molecule inhibitors of the steroidogenic acute regulatory protein (STARD1) by structure-based design

A homology model of the steroidogenic acute regulatory protein (STAR)-related lipid transfer (START) domain of STARD1 was built, and the cholesterol binding site was identified. Structure-based design studies were performed to identify small molecule inhibitors of the START domain. The lead compounds were selected based on cAMP-induced, but not 22R-hydroxycholesterol-supported, inhibition of steroid synthesis by 50% at 10 μM. The results obtained by molecular docking & dynamics show a good correlation between bioactivity, docking scores and calculated binding energies of ligand-protein complexes. The best active compounds will be optimized further and used to develop potential drugs to control excessive steroid formation.     

基于处方挖掘与分子动力学模拟筛选严重急性呼吸综合征冠状病毒2潜在抑制剂分子的研究

目的 从中药筛选具有潜在抑制严重急性呼吸综合征冠状病毒2 （SARS-CoV-2） 活性的成分，进一步从原子水平揭 示其抑制SARS-CoV-2 表面刺突蛋白（S 蛋白） 受体结合域（RBD） 与血管紧张素转化酶2 （ACE2） 结合的内在机制。 方法 检索新型冠状病毒（简称“新冠肺炎”） 治疗中药处方，构建“新冠肺炎中药候选活性成分数据库”。用具有ACE2 抑制活性的小分子化合物构建HipHop药效团模型，并对“新冠肺炎中药候选活性成分数据库”中活性成分筛选。采用分子 对接和分子动力学模拟方法研究候选活性成分与ACE2 的结合方式及其对SARS-CoV-2 S 蛋白与ACE2 识别的影响。 结果 本文通过中药处方挖掘和分子动力学模拟，从143 个新冠肺炎治疗中药处方中筛选出10 种可与SARS-CoV-2 S 蛋白/ 人源ACE2 识别位点结合的中药成分。其中，枇杷叶主要活性成分23-trans-p-coumaryhormentic acid 与ACE2 具有最高的亲和 力，且23-trans-p-coumaryhormentic acid 的结合可有效阻断SARS-CoV-2 S蛋白与宿主细胞ACE2 的结合。结论 本文通过虚 拟筛选发现了SARS-CoV-2 潜在抑制剂分子23-trans-p-coumaryhormentic acid，同时从原子水平预测了其抑制SARS-CoV-2 S 蛋白与ACE2 结合的内在机制，这将为SARS-CoV-2 特异性抗病毒药物的研发提供理论依据。     

Application of molecular modeling to analysis of inhibition of kinesin motor proteins of the BimC subfamily by monastrol and related compounds

Application of molecular modeling approaches has potential to contribute to rational drug design. These approaches may be especially useful when attempting to elucidate the structural features associated with novel drug targets. In this study, molecular docking and molecular dynamics were applied to studies of inhibition of the human motor protein denoted HsEg5 and other homologues in the BimC subfamily. These proteins are essential for mitosis, so compounds that inhibit their activity may have potential as anticancer therapeutics. The discovery of a small-molecule cell-permeable inhibitor, monastrol, has stimulated research in this area. Interestingly, monastrol is reported to inhibit the human and Xenopus forms of Eg5, but not those from Drosophila and Aspergillus. In this study, homology modeling was used to generate models of the Xenopus, Drosophila, and Aspergillus homologues, using the crystal structure of the human protein in complex with monastrol as a template. A series of known inhibitors was docked into each of the homologues, and the differences in binding energies were consistent with reported experimental data. Molecular dynamics revealed significant changes in the structure of the Aspergillus homologue that may contribute to its relative insensitivity to monastrol and related compounds.     

DFT study of the Lewis acid mediated synthesis of 3-acyltetramic acids

The recent emergence and re-emergence of alphaviruses, in particular the chikungunya virus (CHIKV), in numerous countries has invoked a worldwide threat to human health, while simultaneously generating an economic burden on affected countries. There are currently no vaccines or effective drugs available for the treatment of the CHIKV, and with few lead compounds reported, the vital medicinal chemistry is significantly more challenging. This study reports on the discovery of potential inhibitors for the nsP3 macro domain of CHIKV using molecular docking, virtual screening, and molecular dynamics simulations, as well as work done to evaluate and confirm the active site of nsP3. Virtual screening was carried out based on blind docking as well as focused docking, using the database of 1541 compounds from NCI Diversity Set II, to identify hit compounds for nsP3. The top hit compounds were further subjected to molecular dynamic simulations, yielding a greater understanding of the dynamic behavior of nsP3 and its complexes with various ligands, concurrently confirming the outcomes of docking, and establishing in silico lead compounds which target the CHIKV nsP3 enzyme.

Molecular insights of SAH enzyme catalysis and implication for inhibitor design

Biological transmethylation reaction is a key step in the duplication of virus life cycle, in which S-adenosylmethionine plays as the methyl donor. The product of this reactions, S-adenosylhomocysteine (AdoHcy) inhibits the transmethylation process. AdoHcy is hydrolysed to adenosine and L-homocysteine by the action of S-adenosylhomocysteine hydrolase (SAH). Thus the virus life cycle should be cut off once the action of SAH is inhibited. Our study was focussed on the discovery of potential inhibitor against SAH. We performed a similarity search in Traditional Chinese Medicine Database and retrieved 17 hits with high similarity. After that we virtually docked the 17 compounds as well as the natural substrates to the hydrolase using Autodock 3.0.1 software. Then we discussed about the mechanism of the inhibition reaction, followed by proposing the potential inhibitors by comparing best docked solutions and possible modification for the best inhibitors.     

Structure-based screening and molecular dynamics simulations offer novel natural compounds as potential inhibitors of Mycobacterium tuberculosis isocitrate lyase

Mycobacterium tuberculosis is the etiological agent of tuberculosis in humans and is responsible for more than two million deaths annually. M. tuberculosis isocitrate lyase (MtbICL) catalyzes the first step in the glyoxylate cycle, plays a pivotal role in the persistence of M. tuberculosis, which acts as a potential target for an anti-tubercular drug. To identify the potential anti-tuberculosis compound, we conducted a structure-based virtual screening of natural compounds from the ZINC database (n = 1,67,748) against the MtbICL structure. The ligands were docked against MtbICL in three sequential docking modes that resulted in 340 ligands having better docking score. These compounds were evaluated for Lipinski and ADMET prediction, and 27 compounds were found to fit well with re-docking studies. After refinement by molecular docking and drug-likeness analyses, three potential inhibitors (ZINC1306071, ZINC2111081, and ZINC2134917) were identified. These three ligands and the reference compounds were further subjected to molecular dynamics simulation and binding energy analyses to compare the dynamic structure of protein after ligand binding and the stability of the MtbICL and bound complexes. The binding free energy analyses were calculated to validate and capture the intermolecular interactions. The results suggested that the three compounds had a negative binding energy with ?96.462, ?143.549, and ?122.526 kJ mol^?1 for compounds with IDs ZINC1306071, ZINC2111081, and ZINC2134917, respectively. These lead compounds displayed substantial pharmacological and structural properties to be drug candidates. We concluded that ZINC2111081 has a great potential to inhibit MtbICL and would add to the drug discovery process against tuberculosis.     

Preventing the return of smallpox: molecular modeling studies on thymidylate kinase from Variola virus

Smallpox was one of the most devastating diseases in the human history and still represents a serious menace today due to its potential use by bioterrorists. Considering this threat and the non-existence of effective chemotherapy, we propose the enzyme thymidylate kinase from Variola virus (VarTMPK) as a potential target to the drug design against smallpox. We first built a homology model for VarTMPK and performed molecular docking studies on it in order to investigate the interactions with inhibitors of Vaccinia virus TMPK (VacTMPK). Subsequently, molecular dynamics (MD) simulations of these compounds inside VarTMPK and human TMPK (HssTMPK) were carried out in order to select the most promising and selective compounds as leads for the design of potential VarTMPK inhibitors. Results of the docking and MD simulations corroborated to each other, suggesting selectivity towards VarTMPK and, also, a good correlation with the experimental data.