Docking-based 3D-QSAR study of pyridyl aminothiazole derivatives as checkpoint kinase 1 inhibitors

Checkpoint kinase 1 (Chk1) is a promising target for the design of novel anticancer agents. In the present work, molecular docking simulations and three-dimensional quantitative structure–activity relationship (3D-QSAR) studies were performed on pyridyl aminothiazole derivatives as Chk1 inhibitors. AutoDock was used to determine the probable binding conformations of all the compounds inside the active site of Chk1. Comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) models were developed based on the docking conformations and alignments. The CoMFA model produced statistically significant results with a cross-validated correlation coefficient (q²) of 0.608 and a coefficient of determination (r²) of 0.972. The reliable CoMSIA model with q² of 0.662 and r² of 0.970 was obtained from the combination of steric, electrostatic and hydrogen bond acceptor fields. The predictive power of the models were assessed using an external test set of 14 compounds and showed reasonable external predictabilities (r²_pred) of 0.668 and 0.641 for CoMFA and CoMSIA models, respectively. The models were further evaluated by leave-ten-out cross-validation, bootstrapping and progressive scrambling analyses. The study provides valuable information about the key structural elements that are required in the rational design of potential drug candidates of this class of Chk1 inhibitors.     

3D-QSAR and molecular docking study of LRRK2 kinase inhibitors by CoMFA and CoMSIA methods

Three-dimensional quantitative structure–activity relationship (3D-QSAR) modelling was conducted on a series of leucine-rich repeat kinase 2 (LRRK2) antagonists using CoMFA and CoMSIA methods. The data set, which consisted of 37 molecules, was divided into training and test subsets by using a hierarchical clustering method. Both CoMFA and CoMSIA models were derived using a training set on the basis of the common substructure-based alignment. The optimum PLS model built by CoMFA and CoMSIA provided satisfactory statistical results (q² = 0.589 and r² = 0.927 and q² = 0.473 and r² = 0.802, respectively). The external predictive ability of the models was evaluated by using seven compounds. Moreover, an external evaluation set with known experimental data was used to evaluate the external predictive ability of the porposed models. The statistical parameters indicated that CoMFA (after region focusing) has high predictive ability in comparison with standard CoMFA and CoMSIA models. Molecular docking was also performed on the most active compound to investigate the existence of interactions between the most active inhibitor and the LRRK2 receptor. Based on the obtained results and CoMFA contour maps, some features were introduced to provide useful insights for designing novel and potent LRRK2 inhibitors.     

CoMFA and CoMSIA 3D-quantitative structure-activity relationship model on benzodiazepine derivatives, inhibitors of phosphodiesterase IV

Recently, we reported structurally novel PDE4 inhibitors based on 1,4-benzodiazepine derivatives. The main interest in developing bezodiazepine-based PDE4 inhibitors is in their lack of adverse effects of emesis with respect to rolipram-like compounds. A large effort has thus been made toward the structural optimization of this series. In the absence of structural information on the inhibitor binding mode into the PDE4 active site, 2D-QSAR (H-QSAR) and two 3D-QSAR (CoMFA and CoMSIA) methods were applied to improve our understanding of the molecular mechanism controlling the PDE4 affinity of the benzodiazepine derivatives. As expected, the CoMSIA 3D contour maps have provided more information on the benzodiazepine interaction mode with the PDE4 active site whereas CoMFA has built the best tool for activity prediction. The 2D pharmacophoric model derived from CoMSIA fields is consistent with the crystal structure of the PDE4 active site reported recently. The combination of the 2D and 3D-QSAR models was used not only to predict new compounds from the structural optimization process, but also to screen a large library of bezodiazepine derivatives.     

Structure-based CoMFA and CoMSIA study of indolinone inhibitors of PDK1

Phosphoinositide-dependent protein kinase-1 (PDK1) is a Ser/Thr kinase which phosphorylates and activates members of the AGC kinase group known to control processes such as tumor cell growth, protection from apoptosis, and tumor angiogenesis. In this paper, CoMFA and CoMSIA studies were carried out on a training set of 56 conformationally rigid indolinone inhibitors of PDK1. Predictive 3D QSAR models, established using atom fit alignment rule based on crystallographic-bound conformation, had cross-validated (r _cv²) values of 0.738 and 0.816 and non-cross-validated (r _ncv²) values of 0.912 and 0.949 for CoMFA and CoMSIA models, respectively. The predictive ability of the CoMFA and CoMSIA models was determined using a test set of 14 compounds, which gave predictive correlation coefficients (r _pred²) of 0.865 and 0.837, respectively. Structure-based interpretation of the CoMFA and CoMSIA field properties provided further insights for the rational design of new PDK1 inhibitors. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Three-dimensional Quantitative Structure-activity Relationship Models of HIV-1 Integrase Inhibitors of DKAs

As one of the three viral encoded enzymes of HIV-1 infection, HIV-1 integrase has become an attractive drug target for the treatment. Diketoacid compounds (DKAs) are one kind of potent and selective inhibitors of HIV-1 IN. In the present work, two three-dimensional QSAR techniques (CoMFA and CoMSIA) were employed to correlate the molecular structure with the activity of inhibiting the strand transfer for 147 DKAs. The all-oritation search (AOS) and all-placement search (APS) were used to optimize the CoMFA model. The diketo and keto-enol tautomers of DKAs were also used to establish the CoMFA models. The results indicated that the enol was the dominant conformation in the HIV-1 IN and DKAs complexes. It can provide a new method and reference to identify the bioactive conformation of drugs by using QSAR analysis. The best CoMSIA model, with five fields combined, implied that the hydrophobic field is very important as well as the steric and electrostatic fields. All models indicated favorable internal validation. A comparative analysis with the three models demonstrated that the CoMFA model seems to be more predictive. The contour maps could afford steric, electrostatic, hydrophobic and H-bond information about the interaction of ligand-receptor complex visually. The models would give some useful guidelines for designing novel and potent HIV-1 integrase inhibitors.     

Computational Analysis of CRTh2 receptor antagonist: A Ligand-based CoMFA and CoMSIA approach

CRTh2 receptor is an important mediator of inflammatory effects and has attracted much attention as a therapeutic target for the treatment of conditions such as asthma, COPD, allergic rhinitis and atopic dermatitis. In pursuit of better CRTh2 receptor antagonist agents, 3D-QSAR studies were performed on a series of 2-(2-(benzylthio)-1H-benzo[d]imidazol-1-yl) acetic acids. There is no crystal structure information available on this protein; hence in this work, ligand-based comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) were performed by atom by atom matching alignment using systematic search and simulated annealing methods. The 3D-QSAR models were generated with 10 different combinations of test and training set molecules, since the robustness and predictive ability of the model is very important. We have generated 20 models for CoMFA and 100 models for CoMSIA based on two different alignments. Each model was validated with statistical cut off values such as q² > 0.4, r² > 0.5 and r²_pred > 0.5. Based on better q² and r²_pred values, the best predictions were obtained for the CoMFA (model 5 q² = 0.488, r²_pred = 0.732), and CoMSIA (model 45 q² = 0.525, r²_pred = 0.883) from systematic search conformation alignment. The high correlation between the cross-validated/predicted and experimental activities of a test set revealed that the CoMFA and CoMSIA models were robust. Statistical parameters from the generated QSAR models indicated the data is well fitted and have high predictive ability. The generated models suggest that steric, electrostatic, hydrophobic, H-bond donor and acceptor parameters are important for activity. Our study serves as a guide for further experimental investigations on the synthesis of new CRTh2 antagonist.     

A combined CoMFA and CoMSIA 3D-QSAR study of benzamide type antibacterial inhibitors of the FtsZ protein in drug-resistant Staphylococcus aureus

A major problem today is bacterial resistance to antibiotics and the small number of new therapeutic agents approved in recent years. The development of new antibiotics capable of acting on new targets is urgently required. The filamenting temperature-sensitive Z (FtsZ) bacterial protein is a key biomolecule for bacterial division and survival. This makes FtsZ an attractive new pharmacological target for the development of antibacterial agents. There have been several attempts to develop ligands able to inhibit FtsZ. Despite the large number of synthesized compounds that inhibit the FtsZ protein, there are no quantitative structure–activity relationships (QSAR) that allow for the rational design and synthesis of promising new molecules. We present the first 3D-QSAR study of a large and diverse set of molecules that are able to inhibit the FtsZ bacterial protein. We summarize a set of chemical changes that can be made in the steric, electrostatic, hydrophobic and donor/acceptor hydrogen-bonding properties of the pharmacophore, to generate new bioactive molecules against FtsZ. These results provide a rational guide for the design and synthesis of promising new antibacterial agents, supported by the strong statistical parameters obtained from CoMFA (r²_pred = 0.974) and CoMSIA (r²_pred = 0.980) analyses.     

Comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) of thiazolone derivatives as hepatitis C virus NS5B polymerase allosteric inhibitors

Three-dimensional quantitative structure-activity relationship (3D-QSAR) models for a series of thiazolone derivatives as novel inhibitors bound to the allosteric site of hepatitis C virus (HCV) NS5B polymerase were developed based on CoMFA and CoMSIA analyses. Two different conformations of the template molecule and the combinations of different CoMSIA field/fields were considered to build predictive CoMFA and CoMSIA models. The CoMFA and CoMSIA models with best predictive ability were obtained by the use of the template conformation from X-ray crystal structures. The best CoMFA and CoMSIA models gave q (2) values of 0.621 and 0.685, and r (2) values of 0.950 and 0.940, respectively for the 51 compounds in the training set. The predictive ability of the two models was also validated by using a test set of 16 compounds which gave r (pred) (2) values of 0.685 and 0.822, respectively. The information obtained from the CoMFA and CoMSIA 3D contour maps enables the interpretation of their structure-activity relationship and was also used to the design of several new inhibitors with improved activity.     

QSAR analyses on avian influenza virus neuraminidase inhibitors using CoMFA, CoMSIA, and HQSAR

The recent wide spreading of the H5N1 avian influenza virus (AIV) in Asia, Europe and Africa and its ability to cause fatal infections in human has raised serious concerns about a pending global flu pandemic. Neuraminidase (NA) inhibitors are currently the only option for treatment or prophylaxis in humans infected with this strain. However, drugs currently on the market often meet with rapidly emerging resistant mutants and only have limited application as inadequate supply of synthetic material. To dig out helpful information for designing potent inhibitors with novel structures against the NA, we used automated docking, CoMFA, CoMSIA, and HQSAR methods to investigate the quantitative structure-activity relationship for 126 NA inhibitors (NIs) with great structural diversities and wide range of bioactivities against influenza A virus. Based on the binding conformations discovered via molecular docking into the crystal structure of NA, CoMFA and CoMSIA models were successfully built with the cross-validated q (2) of 0.813 and 0.771, respectively. HQSAR was also carried out as a complementary study in that HQSAR technique does not require 3D information of these compounds and could provide a detailed molecular fragment contribution to the inhibitory activity. These models also show clearly how steric, electrostatic, hydrophobicity, and individual fragments affect the potency of NA inhibitors. In addition, CoMFA and CoMSIA field distributions are found to be in well agreement with the structural characteristics of the corresponding binding sites. Therefore, the final 3D-QSAR models and the information of the inhibitor-enzyme interaction should be useful in developing novel potent NA inhibitors.     

CoMFA and CoMSIA studies on C-aryl glucoside SGLT2 inhibitors as potential anti-diabetic agents

SGLT2 has become a target of therapeutic interest in diabetes research. CoMFA and CoMSIA studies were performed on C-aryl glucoside SGLT2 inhibitors (180 analogues) as potential anti-diabetic agents. Three different alignment strategies were used for the compounds. The best CoMFA and CoMSIA models were obtained by means of Distill rigid body alignment of training and test sets, and found statistically significant with cross-validated coefficients (q ²) of 0.602 and 0.618, respectively, and conventional coefficients (r ²) of 0.905 and 0.902, respectively. Both models were validated by a test set of 36 compounds giving satisfactory predicted correlation coefficients (r ² _pred) of 0.622 and 0.584 for CoMFA and CoMSIA models, respectively. A comparison was made with earlier 3D QSAR study on SGLT2 inhibitors, which shows that our 3D QSAR models are better than earlier models to predict good inhibitory activity. CoMFA and CoMSIA models generated in this work can provide useful information to design new compounds and helped in prediction of activity prior to synthesis.     

Pharmacophore and Docking-based 3D-QSAR Studies on HIV-1 Integrase Inhibitors

Integrase（IN） plays an essential role in the process of HIV-1 replication.IN inhibitors of diketo acid derivatives（DKAs） were analysed by the Comparative Molecular Field Analysis（CoMFA） and Comparative Molecular Similarity Induces Analysis（CoMSIA） methods.A set of 42 compounds were randomly selected as the training set（35） and test set（7）.Firstly,a good pharmacophore（goodness of hit=0.787） was obtained and used to align ligands.Then,predictive models were constructed with the CoMFA and CoMSIA methods based on the pharmacophore alignment.As a result,the CoMS1A method yielded the best model with an r2 of 0.955 and a q2 of 0.665,which can predict the activities of the tested DKAs very well（r2=0.559）.Finally,DKAs were docked into IN,and the predicit modes were superimposed on the contour maps obtained from the best CoMSIA model.The superimposed maps gave a visualized and meaningful insight into the inhibitory behaviors,providing significantly useful information for the rational drug design of anti-IN agents.     

Comparative QSAR studies using HQSAR,CoMFA, and CoMSIA methods on cyclic sulfone hydroxyethylamines as BACE1 inhibitors

The inhibition of β-secretase (BACE1) is currently the main pharmacological strategy available for Alzheimer’s disease (AD). 2D QSAR and 3D QSAR analysis on some cyclic sulfone hydroxyethylamines inhibitors against β-secretase (IC₅₀: 0.002–2.75 μM) were carried out using hologram QSAR (HQSAR), comparative molecular field analysis (CoMFA), and comparative molecular similarity indices analysis (CoMSIA) methods. The best model based on the training set was generated with a HQSAR q² value of 0.693 and r² value of 0.981; a CoMFA q² value of 0.534 and r² value of 0.913; and a CoMSIA q² value of 0.512 and r² value of 0.973. In order to gain further understand of the vital interactions between cyclic sulfone hydroxyethylamines and the protease, the analysis was performed by combining the CoMFA and CoMSIA field distributions with the active sites of the BACE1. The final QSAR models could be helpful in the design and development of novel active BACE1 inhibitors.     

CoMFA,CoMSIA, Topomer CoMFA,HQSAR, molecular docking and molecular dynamics simulations study of triazine morpholino derivatives as mTOR inhibitors for the treatment of breast cancer

mTOR has become a promising target for many types of cancer like breast, lung and renal cell carcinoma. CoMFA, CoMSIA, Topomer CoMFA and HQSAR were performed on the series of 39 triazine morpholino derivatives. CoMFA analysis showed q² value of 0.735, r²_cv value of 0.722 and r²_pred value of 0.769. CoMSIA analysis (SEHD) showed q² value of 0.761, r²_cv value of 0.775 and r²_pred value of 0.651. Topomer CoMFA analysis showed q² value of 0.693, r² (conventional correlation coefficient) value of 0.940 and r²_pred value of 0.720. HQSAR analysis showed q²,r²and r²_pred values of 0.694, 0.920 and 0.750, respectively. HQSAR analysis with the combination of atomic number (A), bond type (B) and atomic connections showed q² and r² values of 0.655 and 0.891, respectively. Contour maps from all studies provided significant insights. Molecular docking studies with molecular dynamics simulations were carried out on the highly potent compound 36. Furthermore, four acridine derivatives were designed and docking results of these designed compounds showed the same interactions as that of the standard PI-103 which proved the efficiency of 3D-QSAR and MD/MS study. In future, this study might be useful prior to synthesis for the designing of novel mTOR inhibitors.     

De novo design based identification of potential HIV-1 integrase inhibitors: A pharmacoinformatics study

In the present study, pharmacoinformatics paradigms include receptor-based de novo design, virtual screening through molecular docking and molecular dynamics (MD) simulation are implemented to identify novel and promising HIV-1 integrase inhibitors. The de novodrug/ligand/molecule design is a powerful and effective approach to design a large number of novel and structurally diverse compounds with the required pharmacological profiles. A crystal structure of HIV-1 integrase bound with standard inhibitor BI-224436 is used and a set of 80,000 compounds through the de novo approach in LigBuilder is designed. Initially, a number of criteria including molecular docking, in-silico toxicity and pharmacokinetics profile assessments are implied to reduce the chemical space. Finally, four de novo designed molecules are proposed as potential HIV-1 integrase inhibitors based on comparative analyses. Notably, strong binding interactions have been identified between a few newly identified catalytic amino acid residues and proposed HIV-1 integrase inhibitors. For evaluation of the dynamic stability of the protein-ligand complexes, a number of parameters are explored from the 100 ns MD simulation study. The MD simulation study suggested that proposed molecules efficiently retained their molecular interaction and structural integrity inside the HIV-1 integrase. The binding free energy is calculated through the Molecular Mechanics Poisson-Boltzmann Surface Area (MM-PBSA) approach for all complexes and it also explains their thermodynamic stability. Hence, proposed molecules through de novo design might be critical to inhibiting the HIV-1 integrase.     

DATA类逆转录酶抑制剂的三维定量构效关系

采用对接方法得到HIV-1抑制剂DATA(二芳基三嗪类)分子的活性构象, 进一步用比较分子场分析(CoMFA)和比较分子相似性分析(CoMSIA)法对DATA类逆转录酶抑制剂(RTIs)的三维定量构效关系(3D-QSAR)进行了研究, 建立3D-QSAR模型, 以指导进一步结构修饰. 用此模型预测了5个DATA类似物, 预测偏差较小, 表明了所建立的模型具有较强的预测能力.