Molecular description of pyrimidine-based inhibitors with activity against FAK combining 3D-QSAR analysis,molecular docking and molecular dynamics

Focal adhesion kinase (FAK) is a promising target for developing more effective anticancer drugs. To better understand the structure-activity relationships and mechanism of actions of FAK inhibitors, a molecular modeling study using 3D-QSAR, molecular docking, molecular dynamics simulations, and binding free energy analysis were conducted. Two types of satisfactory 3D-QSAR models were generated, comprising the CoMFA model (R²_cv = 0.528, R²_pred = 0.7557) and CoMSIA model (R²_cv = 0.757, R²_pred = 0.8362), for predicting the inhibitory activities of novel inhibitors. The derived contour maps indicate structural characteristics for substituents on the template. Molecular docking, molecular dynamic simulations and binding free energy calculations further reveal that the binding of inhibitors to FAK is mainly contributed from hydrophobic, electrostatic and hydrogen bonding interactions. In addition, some key residues (Arg14, Glu88, Cys90, Arg138, Asn139, Leu141, and Leu155) responsible for ligand-receptor binding are highlighted. All structural information obtained from 3D-QSAR models and molecular dynamics is consist with the available experimental activities. All the results will facilitate the optimization of this series of FAK inhibitors with higher inhibitory activities.     

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,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.     

3D-QSAR study,docking molecular and simulation dynamic on series of benzimidazole derivatives as anti-cancer agents

Benzimidazole is an important heterocyclic organic compound which has a structural analogy to nucleotides found in human body and hence is an important pharmacophore in medicinal chemistry. The anti-cancer activities for a diverse set of benzimidazole as anti-cancer agents against breast cancer cell line (MCF7) assay have been subjected to 3D-QSAR (3-Dimensional Quantitative Structural-Activity Relationship) studies. Both CoMFA and CoMSIA models exhibit significant results in terms of statistical parameters as determination coefficients R² > 0.9 and Leave One Out cross-validation determination coefficients Q²> 6. The predictive quality of both 3D QSAR models have been assessed by external validation and Y-randomization test. Five new compounds have been designed and predicted by in silico ADMET method. In the second part, we have used the docking molecular and simulation dynamics (MD) to investigate the bonding interactions and stability of the designed compounds into the Pin1. Then, we have compared them to Trastuzumab and Tamoxifen as a standard inhibitors drug of breast cancer. The designed compounds form stable hydrogen and hydrophobic bonding interactions with the residues Lys63, Gln131, Ser154, Arg 68 and Arg69 of Pin1 receptor during 100 ns as a time of the simulation. The obtained results showed that the new benzimidazole are useful as a template for future design of more potent inhibitors against breast cancer cell lines (MCF7).     

Molecular modelling studies on cinnoline-based BTK inhibitors using docking and structure-based 3D-QSAR

ABSTRACT

BTK inhibitors have been proved as an effective target for B-cell malignancies. Ibrutinib is the most advanced irreversible BTK inhibitor for treating mantle cell lymphoma/chronic lymphocytic leukaemia but with existing drug resistance and adverse effects. To design novel effective and safety reversible BTK inhibitors, 115 newly cinnoline analogues were selected to perform molecular docking and 3D-QSAR study because of the main scaffold similarity to Ibrutinib. Both established CoMFA and CoMSIA models obtained high predictive and satisfactory value. CoMFA/CoMSIA contour maps demonstrated that bulky substitutions are preferred at R₁ and R₃ positions, and introducing hydrophilic and negative electrostatic substitutions at R₁ positions is important for improving BTK inhibitory activities. These results will be useful to provide clues for rationally designing novel and high potency BTK inhibitors.     

马来酰胺类糖原合成酶激酶-3β抑制剂的分子对接和三维定量构效关系

通过分子对接和三维定量构效关系(3D-QSAR)两种方法来确定两类马来酰胺类的糖原合成酶激酶-3β(GSK-3β)抑制剂的结合方式. 首先, 用分子对接确定抑制剂与GSK-3β结合模式及其相互作用; 然后用比较分子力场分析法(CoMFA)与比较分子相似性指数分析法(CoMSIA)对48个化合物做三维定量构效关系的分析. 两种方法得出的交互验证回归系数分别为0.669(CoMFA)和0.683(CoMSIA), 证明该模型具有很好的统计相关性, 同时也说明该模型具有较高的预测能力.根据该模型提供的信息, 设计出9个预测活性较好的分子.     

Molecular docking,dynamics simulations and 3D-QSAR modeling of arylpiperazine derivatives of 3,5-dioxo-(2H,4H)-1,2,4-triazine as 5-HT1AR agonists

Serotonin receptor, 5-HT_1AR, agonists and partial agonists have established drug candidates for psychiatric and neurologic disorders. Recently, we reported the synthesis and evaluation of arylpiperazine derivatives of 3,5-dioxo-(2H,4H)-1,2,4-triazine as 5-HT_1AR ligands. Herein, we generated a homology model of the receptor and docked the ligands against it, predicted the stability of the receptor model and complexes by molecular dynamics and generated a 3D-QSAR model for the arylpiperazine derivatives of 3,5-dioxo-(2H,4H)-1,2,4-triazine. The model suggests the hydrophobic part that arises from the aromatic region and the electron withdrawing parts play a vital role in the agonist activity of the lead molecules.     

3D-QSAR,molecular docking and ADMET studies of thioquinazolinone derivatives against breast cancer

Breast cancer is a deadly disease and the second largest cause of mortality on a worldwide platform. Despite the availability of several cancer treatments, life expectancies stay relatively poor. Consequently, the medicinal chemistry community prioritizes the quick discovery of novel anticancer drugs. In recent years, computational approaches have been widely used to accelerate the drug development process. In light of this, in the current work, we performed three-dimensional quantitative structure-activity relationships (3D-QSAR) and molecular docking analyses on thioquinazolinone derivatives with aromatase enzyme (PDB: 3S7S). External validation was used to validate the prediction capabilities of the generated model. The best CoMSIA (comparative molecular similarity indices analysis) model exhibited the significant values of $Q^2$, $R^{2}and$ $R_{pred}^2$. These findings suggested that the electrostatic, hydrophobic and hydrogen bond donor and acceptor fields have a significant effect on inhibition of breast cancer. Thus, a number of innovative potent aromatase inhibitors were designed and their biological activities were predicted based on the best model. Furthermore, molecular docking studies were carried out for the designed compounds against breast cancer. Additionally, ADMET proprieties were used to evaluate drug-likeness of these novel drug candidates. The most active compounds found by these computational studies could be helpful for synthesis and testing as prospective future anti-cancer treatments.     

The novel 4-hydroxyphenylpyruvate dioxygenase inhibitors in vivo and in silico approach: 3D-QSAR analysis,molecular docking,bioassay and molecular dynamics

4-Hydroxyphenylpyruvate dioxygenase (HPPD) is not only an important target enzyme for the treatment of type I tyrosinemia, but also a new target for design bleaching herbicides, and it plays key role in the biosynthesis of tocopherol and plastoquinone. Thirty-six known active pyridine derivatives were collected, and comparative molecular field analysis (CoMFA) and comparative molecular similarity index analysis (CoMSIA) models based on common skeleton were constructed to obtain novel HPPD herbicides with higher activity. Two new HPPD inhibitors were rationally designed and synthesized according to the CoMFA and CoMSIA models and verified by enzyme activity, biological assays, and molecular docking. The promising compound W1 ((E)-5-(3-(4-bromophenyl)acryloyl)-6-hydroxy-2,3-dihydropyridin-4(1H)-one) showed better AtHPPD inhibitory activity, and the bioassay results revealed that some weeds showed bleaching symptoms. The good binding stability of W1 and protein was confirmed by molecular dynamics simulation in 100 ns. These results would be highly useful in the progress of new HPPD inhibitors discovery.     

Pharmacophore modeling,atom based 3D-QSAR,molecular docking and molecular dynamics studies on Escherichia coli ParE inhibitors

ATP dependent ParE enzyme is as an attractive target for the development of antibacterial agents. Atom based 3D-QSAR model AADHR.187 was developed based on the thirty eight Escherichia coli ParE inhibitors. The generated model showed statistically significant coefficient of determinations for the training (R² = 0.985) and test (R² = 0.86) sets. The cross-validated correlation coefficient (q2) was 0.976. The utility of the generated model was validated by the enrichment study. The model was also validated with structurally diverse external test set of ten compounds. Contour plot analysis of the generated model unveiled the chemical features necessary for the E. coli ParE enzyme inhibition. Extra-precision docking result revealed that hydrogen bonding and ionic interactions play a major role in ParE protein-ligand binding. Binding free energy was computed for the data set inhibitors to validate the binding affinity. A 30-ns molecular dynamics simulation showed high stability and effective binding of inhibitor 34 within the active site of ParE enzyme. Using the best fitted model AADHR.187, pharmacophore-based high-throughput virtual screening was performed to identify virtual hits. Based on the above studies three new molecules are proposed as E. coli ParE inhibitors with high binding affinity and favourable ADME properties.     

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 studies on indole and 7-azoindole derivatives as ROCK-2 inhibitors: An integrative computational approach

Rho Kinases (ROCK) has been found to regulate a wide range of fundamental cell functions such as contraction, motility, proliferation, and apoptosis. Recent experiments have defined new functions of ROCKs in cells, including centrosome positioning and cell-size regulation, which might contribute to various physiological and pathological states. In this study, we have performed pharmacophore modeling and 3D QSAR studies on a series of 36 indoles and 7-azoindoles derivatives as ROCK2 inhibitors to elucidate the structural variations with their inhibitory activities. Ligand based CoMFA and CoMSIA models were generated based on three different alignment methods such as systematic search, simulated annealing and pharmacophore. A total of 15 CoMFA models and 27 CoMSIA were generated using different alignments. One model from each alignment is selected based on the statistical values. Contour maps of the selected models were compared, analysed and reported. The 3D QSAR study revealed that electro positive group linked to the methoxy-benzene ring position of the structure will enhance the biological activity and bulkier substitutions are preferred in the methyl dihydroindole region. Also, it is found that the hydrogen bond donor substituted at the R₁ position enhances the inhibitory activity. In future, this study would give proper guidelines to further enhance the activity of novel inhibitors for ROCK2.     

Glycogen synthase kinase-3 inhibition by 3-anilino-4-phenylmaleimides: insights from 3D-QSAR and docking

Glycogen synthase kinase-3, a serine/threonine kinase, has been implicated in a wide variety of pathological conditions such as diabetes, Alzheimer’s disease, stroke, bipolar disorder, malaria and cancer. Herein we report 3D-QSAR analyses using CoMFA and CoMSIA and molecular docking studies on 3-anilino-4-phenylmaleimides as GSK-3α inhibitors, in order to better understand the mechanism of action and structure-activity relationship of these compounds. Comparison of the active site residues of GSK-3α and GSK-3β isoforms shows that all the key amino acids involved in polar interactions with the maleimides for the β isoform are the same in the α isoform, except that Asp133 in the β isoform is replaced by Glu196 in the α isoform. We prepared a homology model for GSK-3α, and showed that the change from Asp to Glu should not affect maleimide binding significantly. Docking studies revealed the binding poses of three subclasses of these ligands, namely anilino, N-methylanilino and indoline derivatives, within the active site of the β isoform, and helped to explain the difference in their inhibitory activity.     

雌激素β受体喹啉类配体的分子对接及3D-QSAR研究

对33个喹啉衍生物的雌激素β受体活性进行了分子对接以及比较分子力场分析(CoMFA)和比较分子相似性指数分析(CoMSIA). 对接结果显示氢键和疏水作用是配体与受体结合的主要因素,同时结果亦显示对接结合能与观测值pIC50具有极显著的线性相关性. 根据对接后各优势构象将33个样本进行叠合并进行CoMFA与CoMSIA研究,均得到了较优的结果,其中以选用立体场、静电场和疏水场建立的CoMSIA模型结果最优,其主成分数,r2,q2（LOO）和r2pred分别为2, 0.894, 0.708和0.802. 构效关系模型分析显示基团的空间位阻、电性及疏水作用是影响活性的主要因素     

In silico design novel (5-imidazol-2-yl-4-phenylpyrimidin-2-yl)[2-(2-pyridylamino)ethyl]amine derivatives as inhibitors for glycogen synthase kinase 3 based on 3D-QSAR,molecular docking and molecular dynamics simulation

Glycogen Synthase Kinase 3 (GSK-3) is a member of cellular kinase with various functions, such as glucose regulation, cellular differentiation, neuronal function and cell apoptosis. It has been proved as an important therapeutic target in type 2 diabetes mellitus and Alzheimer's disease. To better understand their structure–activity relationships and mechanism of action, an integrated computational study, including three dimensional quantitative structure-activity relationship (3D-QSAR), molecular docking, and molecular dynamics (MD), was performed on 79 (5-Imidazol-2-yl-4-phenylpyrimidin-2-yl)[2-(2-pyridylamino)ethyl]amine GSK-3 inhibitors. In this paper, we constructed 3D-QSAR using comparative molecular field analysis (CoMFA) and comparative molecular similarity index analysis (CoMSIA) method. The results showed that the CoMFA model (q ² = 0.743, r² = 0.980) and the CoMSIA model (q² = 0.813, r² = 0.976) had stable and reliable predictive ability. The electrostatic and H-bond donor fields play important roles in the models. The contour maps of the model visually showed the relationship between the activity of compounds and their three-dimensional structure. Molecular docking was used to identify the key amino acid residues at the active site of GSK-3 and explore its binding mode with ligands. Based on 3D-QSAR models, contour maps and the binding feature between GSK-3 and inhibitor, we designed 10 novel compounds with good potential activity and ADME/T profile. Molecular dynamics simulation results validated that Ile62, Val70 and Lys85 located in the active site play a key role for GSK-3 complexed with inhibitors. These results might provide important information for designing GSK-3 inhibitors with high activity.     

Molecular Modeling Studies of Vascular Endothelial Growth Factor Receptor Tyrosine Kinase Inhibitors Combining Molecular Docking and 3D-QSAR Methods

The vascular endothelial growth factor (VEGF) and its receptor tyrosine kinases VEGFR-2 or kinase insertdomain receptor (KDR) have emerged as attractive targets for the design of novel anticancer agents. In the present work, molecular docking method combined with three dimensional quantitative structure-activity relationships (comparative molecular field analysis (CoMFA) and comparative molecular similarity indice analysis (CoMSIA)) to analyze the possible interactions between KDR and those derivatives which acted as selective inhibitors. The CoMFA and CoMSIA models gave a cross-validated coefficient Q2 of 0.713 and 0.549, non-cross-validated R2 values of 0.974 and 0.878, and predicted R2 values of 0.966 and 0.823, respectively. The 3D contour maps generated by the CoMFA and CoMSIA models were used to identify the key structural requirements responsible for the biological activity. The information obtained from 3D-QSAR and docking studies were very helpful to design novel selective inhibitors of KDR with desired activity and good chemical property.