Identification of potential CRAC channel inhibitors: Pharmacophore mapping, 3D-QSAR modelling,and molecular docking approach

Upregulation of store-operated Ca²⁺ influx via ORAI1, an integral component of the CRAC channel, is responsible for abnormal cytokine release in active rheumatoid arthritis, and therefore ORAI1 has been proposed as an attractive molecular target. In this study, we attempted to predict the mechanical insights of ORAI1 inhibitors through pharmacophore modelling, 3D-QSAR, molecular docking and free energy analysis. Various hypotheses of pharmacophores were generated and from that, a pharmacophore hypothesis with two hydrogen bond acceptors, one hydrogen bond donor and two aromatic rings (AADRR) resulted in a statistically significant 3D-QSAR model (r² = 0.84 and q² = 0.74). We believe that the obtained statistical model is a reliable QSAR model for the diverse dataset of inhibitors against the IL-2 production assay. The visualization of contours in active and inactive compounds generated from the 3D-QSAR models and molecular docking studies revealed major interaction with GLN108, HIS113 and ASP114, and interestingly, these residues are located near the Ca²⁺ selectivity filter region. Free energy binding analysis revealed that Coulomb energy, van der Waals energy and non-polar solvation terms are more favourable for ligand binding. Thus, the present study provides the physical and chemical requirements for the development of novel ORAI1 inhibitors with improved biological activity.     

Studying the Binding Modes of Novel 2-Aminopyridine Derivatives as Effective and Selective c-Met Kinase Type 1 Inhibitors Using Molecular Modeling Approaches

The mesenchymal epithelial cell transforming factor c-Met, encoded by c-Met proto-oncogene and known as a high-affinity receptor for Hepatocyte Growth Factor (HGF), is one of the receptor tyrosine kinases (RTKs) members. The HGF/c-Met signaling pathway has close correlation with tumor growth, invasion and metastasis. Thus, c-Met kinase has emerged as a prominent therapeutic target for cancer drug discovery. Recently a series of novel 2-aminopyridine derivatives targeting c-Met kinase with high biological activity were reported. In this study, 3D quantitative structure-activity relationship (QSAR), molecular docking and molecular dynamics simulations (MD) were employed to research the binding modes of these inhibitors.The results show that both the atom-based and docking-based CoMFA (Q² = 0.596, R² = 0.950 in atom-based model and Q² = 0.563, R² = 0.985 in docking-based model) and CoMSIA (Q² = 0.646, R² = 0.931 in atom-based model and Q² = 0.568, R² = 0.983 in docking-based model) models own satisfactory performance with good reliabilities and powerful external predictabilities. Molecular docking study suggests that Tyr1230 and Arg1208 might be the key residues, and electrostatic and hydrogen bond interactions were shown to be vital to the activity, concordance with QSAR analysis. Then MD simulation was performed to further explore the binding mode of the most potent inhibitor. The obtained results provide important references for further rational design of c-Met Kinase type I inhibitors.     

Identification of potential HIV-1 integrase strand transfer inhibitors: In silico virtual screening and QM/MM docking studies

HIV-1 integrase (IN) is a retroviral enzyme that catalyses integration of the reverse-transcribed viral DNA into the host genome, which is necessary for efficient viral replication. In this study, we have performed an in silico virtual screening for the identification of potential HIV-1 IN strand transfer (ST) inhibitors. Pharmacophore modelling and atom-based 3D-QSAR studies were carried out for a series of compounds belonging to 3-Hydroxypyrimidine-2,4-diones. Based on the ligand-based pharmacophore model, we obtained a five-point pharmacophore with two hydrogen bond acceptors (A), one hydrogen bond donor (D), one hydrophobic group (H) and one aromatic ring (R) as pharmacophoric features. The pharmacophore hypothesis AADHR was used as a 3D query in a sequential virtual screening study to filter small molecule databases Maybridge, ChemBridge and Asinex. Hits matching with pharmacophore hypothesis AADHR were retrieved and passed progressively through Lipinski’s rule of five filtering, molecular docking and hierarchical clustering. The five compounds with best hits with novel and diverse chemotypes were subjected to QM/MM docking, which showed improved docking accuracy. We further performed molecular dynamics simulation and found three compounds that form stable interactions with key residues. These compounds could be used as a leads for further drug development and rational design of HIV-1 IN inhibitors.     

Pharmacophore modelling,atom-based 3D-QSAR generation and virtual screening of molecules projected for mPGES-1 inhibitory activity

COX-2 inhibitors exhibit anticancer effects in various cancer models but due to the adverse side effects associated with these inhibitors, targeting molecules downstream of COX-2 (such as mPGES-1) has been suggested. Even after calls for mPGES-1 inhibitor design, to date there are only a few published inhibitors targeting the enzyme and displaying anticancer activity. In the present study, we have deployed both ligand and structure-based drug design approaches to hunt novel drug-like candidates as mPGES-1 inhibitors. Fifty-four compounds with tested mPGES-1 inhibitory value were used to develop a model with four pharmacophoric features. 3D-QSAR studies were undertaken to check the robustness of the model. Statistical parameters such as r² = 0.9924, q² = 0.5761 and F test = 1139.7 indicated significant predictive ability of the proposed model. Our QSAR model exhibits sites where a hydrogen bond donor, hydrophobic group and the aromatic ring can be substituted so as to enhance the efficacy of the inhibitor. Furthermore, we used our validated pharmacophore model as a three-dimensional query to screen the FDA-approved Lopac database. Finally, five compounds were selected as potent mPGES-1 inhibitors on the basis of their docking energy and pharmacokinetic properties such as ADME and Lipinski rule of five.     

Imidazo[1,2-a]pyrazine inhibitors of phosphoinositide 3-kinase alpha (PI3Kα): 3D-QSAR analysis utilizing the Hybrid Monte Carlo algorithm to refine receptor-ligand complexes for molecular alignment

Phosphoinositide 3-kinase alpha (PI3Kα) is a lipid kinase involved in several cellular functions such as cell growth, proliferation, differentiation and survival, and its anomalous regulation leads to cancerous conditions. PI3Kα inhibition completely blocks the cancer signalling pathway, hence it can be explored as an important therapeutic target for cancer treatment. In the present study, docking analysis of 49 selective imidazo[1,2-a]pyrazine inhibitors of PI3Kα was carried out using the QM-Polarized ligand docking (QPLD) program of the Schrödinger software, followed by the refinement of receptor–ligand conformations using the Hybrid Monte Carlo algorithm in the Liaison program, and alignment of refined conformations of inhibitors was utilized for the development of an atom-based 3D-QSAR model in the PHASE program. Among the five generated models, the best model was selected corresponding to PLS factor 2, displaying the highest value of Q²_test (0.650). The selected model also displayed high values of r²_train (0.917), F-value (166.5) and Pearson-r (0.877) and a low value of SD (0.265). The contour plots generated for the selected 3D-QSAR model were correlated with the results of docking simulations. Finally, this combined information generated from 3D-QSAR and docking analysis was used to design new congeners.     

Three-dimensional quantitative structure-activity relationship (3D QSAR) and pharmacophore elucidation of tetrahydropyran derivatives as serotonin and norepinephrine transporter inhibitors

Three-dimensional quantitative structure-activity relationship (3D QSAR) using comparative molecular field analysis (CoMFA) was performed on a series of substituted tetrahydropyran (THP) derivatives possessing serotonin (SERT) and norepinephrine (NET) transporter inhibitory activities. The study aimed to rationalize the potency of these inhibitors for SERT and NET as well as the observed selectivity differences for NET over SERT. The dataset consisted of 29 molecules, of which 23 molecules were used as the training set for deriving CoMFA models for SERT and NET uptake inhibitory activities. Superimpositions were performed using atom-based fitting and 3-point pharmacophore-based alignment. Two charge calculation methods, Gasteiger-Hückel and semiempirical PM3, were tried. Both alignment methods were analyzed in terms of their predictive abilities and produced comparable results with high internal and external predictivities. The models obtained using the 3-point pharmacophore-based alignment outperformed the models with atom-based fitting in terms of relevant statistics and interpretability of the generated contour maps. Steric fields dominated electrostatic fields in terms of contribution. The selectivity analysis (NET over SERT), though yielded models with good internal predictivity, showed very poor external test set predictions. The analysis was repeated with 24 molecules after systematically excluding so-called outliers (5 out of 29) from the model derivation process. The resulting CoMFA model using the atom-based fitting exhibited good statistics and was able to explain most of the selectivity (NET over SERT)-discriminating factors. The presence of −OH substituent on the THP ring was found to be one of the most important factors governing the NET selectivity over SERT. Thus, a 4-point NET-selective pharmacophore, after introducing this newly found H-bond donor/acceptor feature in addition to the initial 3-point pharmacophore, was proposed.     

Molecular dynamics and pharmacophore modelling studies of different subtype (ALK and EGFR (T790M)) inhibitors in NSCLC

Extensively validated 3D pharmacophore models for ALK (anaplastic lymphoma kinase) and EGFR (T790M) (epithelial growth factor receptor with acquired secondary mutation) were developed. The pharmacophore model for ALK (r² = 0.96, q² = 0.692) suggested that two hydrogen bond acceptors and three hydrophobic groups arranged in 3-D space are essential for the binding affinity of ALK inhibitors. Similarly, the pharmacophore model for EGFR (T790M) (r² = 0.92, q² = 0.72) suggested that the presence of a hydrogen bond acceptor, two hydrogen bond donors and a hydrophobic group plays vital role in binding of an inhibitor of EGFR (T790M). These pharmacophore models allowed searches for novel ALK and EGFR (T790M) dual inhibitors from multiconformer 3D databases (Asinex, Chembridge and Maybridge). Finally, the eight best hits were selected for molecular dynamics simulation, to study the stability of their complexes with both proteins and final binding orientations of these molecules. After molecular dynamics simulations, one hit has been predicted to possess good binding affinity for both ALK and EGFR (T790M), which can be further investigated for its experimental in-vitro/in-vivo activities.     

Pharmacophore modelling and electronic feature analysis of hydroxamic acid derivatives,the HIV integrase inhibitors

Hydroxamic acid derivatives with metal ion binding properties were collected from the literature to generate a pharmacophore and 3D-QSAR model for HIV strand transfer inhibition. The derived pharmacophore model (AAAHRR) recognizes both metal ion binding site and hydrophobic group. The QSAR model generated using this hypothesis expressed statistical significance (r ² = 0.971 for the training set and q ² = 0.913 for the test set). The ability of this pharmacophore model to retrieve other metal ion binding inhibitors was examined by screening the ChemBank database (ligandinfo) incorporated with 10 known strand transfer inhibitors. The studied favourable and unfavourable contours of chemical features (H-bond donor, acceptor and hydrophobic sites) revealed the role of hydrophobic substitution at the fluorobenzene ring and cyclization of the metal ion binding hydroxamic acid in effective integrase inhibition. Analysis of the frontier orbitals, HOMO and LUMO revealed that the nucleophilic / electrophilic interactions depend on the significant overlapping observed at the azaindole and hydroxamic acid groups. In essence, the generated pharmacophore model is competent enough to disclose the essential site-specific interactions involved in the inhibition of HIV integrase, and hence can be used in virtual screening to identify novel scaffolds as leads with increased anti-viral potency.     

Modeling ligand-receptor interaction for some MHC class II HLA-DR4 peptide mimetic inhibitors using several molecular docking and 3D QSAR techniques

The ligand-receptor interaction between some peptidomimetic inhibitors and a class II MHC peptide presenting molecule, the HLA-DR4 receptor, was modeled using some three-dimensional (3D) quantitative structure-activity relationship (QSAR) methods such as the Comparative Molecular Field Analysis (CoMFA), Comparative Molecular Similarity Indices Analysis (CoMSIA), and a pharmacophore building method, the Catalyst program. The structures of these peptidomimetic inhibitors were generated theoretically, and the conformations used in the 3D QSAR studies were defined by docking them into the known structure of HLA-DR4 receptor through the GOLD, GLIDE Rigidly, GLIDE Flexible, and Xscore programs. Some of the parameters used in these docking programs were selected by docking an X-ray ligand into the receptor and comparing the root-means-square difference (RMSD) computed between the coordinates of the X-ray and docked structure. However, the goodness of a docking result for docking a series of peptidomimetic inhibitors into the HLA-DR4 receptor was judged by comparing the Spearman's rank correlation coefficient computed between each docking result and the activity data taken from the literature. The best CoMFA and CoMSIA models were constructed using the aligned structures of the best docking result. The CoMSIA was conducted in a stepwise manner to identify some important molecular features that were further employed in a pharmacophore building process by the Catalyst program. It was found that most inhibitors of the training set were accurately predicted by the best pharmacophore model, the Hypo1 hypothesis constructed. The deviation or conflict found between the actual and predicted activities of some inhibitors of both the training and the test sets were also investigated by mapping the Hypo1 hypothesis onto the corresponding structures of the inhibitors.     

Computational Approaches for the Design of Novel Anticancer Compounds Based on Pyrazolo[3,4-d]pyrimidine Derivatives as TRAP1 Inhibitor

In the present in-silico study, various computational techniques were applied to determine potent compounds against TRAP1 kinase. The pharmacophore hypothesis DHHRR_1 consists of important features required for activity. The 3D QSAR study showed a statistically significant model with R² = 0.96 and Q² = 0.57. Leave one out (LOO) cross-validation (R² CV = 0.58) was used to validate the QSAR model. The molecular docking study showed maximum XP docking scores (−11.265, −10.532, −10.422, −10.827, −10.753 kcal/mol) for potent pyrazole analogs (42, 46, 49, 56, 43), respectively, with significant interactions with amino acid residues (ASP 594, CYS 532, PHE 583, SER 536) against TRAP1 kinase receptors (PDB ID: 5Y3N). Furthermore, the docking results were validated using the 100 ns MD simulations performed for the selected five docked complexes. The selected inhibitors showed relatively higher binding affinities than the TRAP1 inhibitor molecules present in the literature. The ZINC database was used for a virtual screening study that screened ZINC05297837, ZINC05434822, and ZINC72286418, which showed similar binding interactions to those shown by potent ligands. Absorption, distribution, metabolism, and excretion (ADME) analysis showed noticeable results. The results of the study may be helpful for the further development of potent TRAP1 inhibitors     

Docking Study and Three‐Dimensional Quantitative Structure‐Activity Relationship (3D‐QSAR) Analyses and Novel Molecular Design of a Series of 4‐Aminoquinazolines as Inhibitors of Aurora B Kinase

The Aurora proteins are critical regulators of major mitotic events and attractive targets for anticancer therapy. 3D‐QSAR studies based on molecular docking were performed on a dataset of 40 4‐aminoquinazolines compounds. The CoMSIA model produced significantly better results than CoMFA model, with q²=0.652 and r²=0.991. The contours analysis provides useful information about the structural requirements for 4‐aminoquinazolines for inhibiting Aurora B. Scaffold hopping method was used to generate new structures based on the maximum common substructure of the training and test set compounds. The ADMET property, binding affinity and inhibitory activity of the new designed compounds were predicted, respectively. Finally 16 compounds were identified as the novel inhibitors for Aurora B kinase.     

Binding mode of inhibitors and Cryptosporidium parvum IMP dehydrogenase: A combined ligand- and receptor-based study

A combined ligand- and target-based approach was used to analyse the interaction models of Cryptosporidium parvum inosine 5’-monophosphate dehydrogenase (CpIMPDH) with selective inhibitors. First, a ligand-based pharmacophore model was generated from 20 NAD⁺ competitive CpIMPDH inhibitors with the HipHop module. The characteristic of the NAD⁺ binding site of CpIMPDH was then described, and the binding modes of the representative inhibitors were studied by molecular docking. The combination of the pharmacophore model and the docking results allowed us to evaluate the pharmacophore features and structural information of the NAD⁺ binding site of CpIMPDH. This research supports the proposal of an interaction model inside the NAD⁺ binding site of CpIMPDH, consisting of four key interaction points: two hydrophobic-aromatic groups, a hydrophobic-aliphatic group and a hydrogen bond donor. This study also provides guidance for the design of more potent CpIMPDH inhibitors for the treatment of Cryptosporidium infections.     

Strategies to design pyrazolyl urea derivatives for p38 kinase inhibition: a molecular modeling study

The p38 protein kinase is a serine–threonine mitogen activated protein kinase, which plays an important role in inflammation and arthritis. A combined study of 3D-QSAR and molecular docking has been undertaken to explore the structural insights of pyrazolyl urea p38 kinase inhibitors. The 3D-QSAR studies involved comparative molecular field analysis (CoMFA) and comparative molecular similarity indices (CoMSIA). The best CoMFA model was derived from the atom fit alignment with a cross-validated r ² (q ²) value of 0.516 and conventional r ² of 0.950, while the best CoMSIA model yielded a q ² of 0.455 and r ² of 0.979 (39 molecules in training set, 9 molecules in test set). The CoMFA and CoMSIA contour maps generated from these models provided inklings about the influence of interactive molecular fields in the space on the activity. GOLD, Sybyl (FlexX) and AutoDock docking protocols were exercised to explore the protein–inhibitor interactions. The integration of 3D-QSAR and molecular docking has proffered essential structural features of pyrazolyl urea inhibitors and also strategies to design new potent analogues with enhanced activity.     

Binding Free Energy (BFE) Calculations and Quantitative Structure–Activity Relationship (QSAR) Analysis of Schistosoma mansoni Histone Deacetylase 8 (smHDAC8) Inhibitors

Histone-modifying proteins have been identified as promising targets to treat several diseases including cancer and parasitic ailments. In silico methods have been incorporated within a variety of drug discovery programs to facilitate the identification and development of novel lead compounds. In this study, we explore the binding modes of a series of benzhydroxamates derivatives developed as histone deacetylase inhibitors of Schistosoma mansoni histone deacetylase (smHDAC) using molecular docking and binding free energy (BFE) calculations. The developed docking protocol was able to correctly reproduce the experimentally established binding modes of resolved smHDAC8–inhibitor complexes. However, as has been reported in former studies, the obtained docking scores weakly correlate with the experimentally determined activity of the studied inhibitors. Thus, the obtained docking poses were refined and rescored using the Amber software. From the computed protein–inhibitor BFE, different quantitative structure–activity relationship (QSAR) models could be developed and validated using several cross-validation techniques. Some of the generated QSAR models with good correlation could explain up to ~73% variance in activity within the studied training set molecules. The best performing models were subsequently tested on an external test set of newly designed and synthesized analogs. In vitro testing showed a good correlation between the predicted and experimentally observed IC₅₀ values. Thus, the generated models can be considered as interesting tools for the identification of novel smHDAC8 inhibitors.     

Discovery and design of new PI3K inhibitors through pharmacophore-based virtual screening,molecular docking,and binding free energy analysis

Phosphoinositide 3-kinase (PI3K)/protein kinase B (AKT)/the mammalian target of rapamycin (mTOR) signaling pathway is one of the sought after therapeutic target for treating human cancers. This pathway is often hyper activated in cancers. In the present study, pharmacophore-based virtual screening, molecular docking, and binding free energy calculations were performed on a series of quinoline derivatives which were reported to be effective against PI3Kα. A five-point pharmacophore hypothesis with one hydrogen bond acceptor (A), one hydrogen bond donor (D), one hydrophobic group (H), and two aromatic rings (R) was developed with acceptable R² and Q² values of 0.93 and 0.60 respectively. Eventually, common pharmacophore hypothesis-based screening was conducted against TOSLab, CPP, and ASINEX macrocylce databases, and potential hits were identified which were further subjected to rigorous docking process in order to screen out drug like molecules having crucial interactions with the target PI3Kα. Finally, binding free energy analysis was carried out for the top hits obtained from docking process. We also designed new 1, 3, 4-oxadiazole-based cyclic peptides by incorporating the structural features of the hits obtained from the above databases. Among the designed cyclic peptides, the cyclic peptide with tryptophan moiety showed good interactions and free binding energy values. On the whole, this study helped us in identifying new promising molecules as PI3Kα inhibitors which can be explored further to generate greater number of compounds with better pharmacokinetic properties.     

Structure-based 3D QSAR and design of novel acetylcholinesterase inhibitors

The paper describes the construction, validation and application of a structure-based 3D QSAR model of novel acetylcholinesterase (AChE) inhibitors. Initial use was made of four X-ray structures of AChE complexed with small, non-specific inhibitors to create a model of the binding of recently developed aminopyridazine derivatives. Combined automated and manual docking methods were applied to dock the co-crystallized inhibitors into the binding pocket. Validation of the modelling process was achieved by comparing the predicted enzyme-bound conformation with the known conformation in the X-ray structure. The successful prediction of the binding conformation of the known inhibitors gave confidence that we could use our model to evaluate the binding conformation of the aminopyridazine compounds. The alignment of 42 aminopyridazine compounds derived by the docking procedure was taken as the basis for a 3D QSAR analysis applying the GRID/GOLPE method. A model of high quality was obtained using the GRID water probe, as confirmed by the cross-validation method (q² _LOO=0.937, q² _{L50% O}=0.910). The validated model, together with the information obtained from the calculated AChE-inhibitor complexes, were considered for the design of novel compounds. Seven designed inhibitors which were synthesized and tested were shown to be highly active. After performing our modelling study the X-ray structure of AChE complexed with donepezil, an inhibitor structurally related to the developed aminopyirdazines, has been made available. The good agreement found between the predicted binding conformation of the aminopyridazines and the one observed for donepezil in the crystal structure further supports our developed model.     

Pharmacophore modelling,molecular docking and virtual screening for EGFR (HER 1) tyrosine kinase inhibitors

A pharmacophore model has been developed using diverse classes of epidermal growth factor receptor (EGFR) tyrosine kinase (TK) inhibitors useful in the treatment of human tumours. Among the top 10 generated hypotheses, the second hypothesis, with one hydrogen bond acceptor, one ring aromatic and three hydrophobic features, was found to be the best on the basis of Cat Scramble validation as well as test set prediction (r _training?=?0.89, r _test?=?0.82). The model also maps well to the external test set molecules as well as clinically active molecules and corroborates the docking studies. Finally, 10 hits were identified as potential leads after virtual screening of ZINC database for EGFR TK inhibition. The study may facilitate the designing and discovery of novel EGFR TK inhibitors.