QSAR models of cytochrome P450 enzyme 1A2 inhibitors using CoMFA,CoMSIA and HQSAR

Quantitative structure–activity relationship (QSAR) studies were conducted on an in-house database of cytochrome P450 enzyme 1A2 inhibitors using the comparative molecular field analysis (CoMFA), comparative molecular similarity analysis (CoMSIA) and hologram QSAR (HQSAR) approaches. The database consisted of 36 active molecules featuring varied core structures. The model based on the naphthalene substructure alignment incorporating 19 molecules yielded the best model with a CoMFA cross validation value q² of 0.667 and a Pearson correlation coefficient r² of 0.976; a CoMSIA q² value of 0.616 and r² value of 0.985; and a HQSAR q² value of 0.652 and r² value of 0.917. A second model incorporating 34 molecules aligned using the benzene substructure yielded an acceptable CoMFA model with q² value of 0.5 and r² value of 0.991. Depending on the core structure of the molecule under consideration, new CYP1A2 inhibitors will be designed based on the results from these models.     

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 conformations,QSAR, and molecular design of Alkene‐3‐quinolinecarbonitriles as Src inhibitors

A theoretical study on binding orientations and quantitative structure–activity relationship (QSAR) of a novel series of alkene‐3‐quinolinecarbonitriles acting as Src inhibitors has been carried out by using the docking study and three‐dimensional QSAR (3D‐QSAR) analyses. The appropriate binding orientations and conformations of these compounds interacting with Src kinase were revealed by the docking studies, and the established 3D‐QSAR models show significant statistical quality and satisfactory predictive ability, with high R² values and q² values: comparative molecular field analysis (CoMFA) model (q² = 0.748, R² = 0.972), comparative molecular similarity indices analysis (CoMSIA) model (q² = 0.731, R² = 0.987). The systemic external validation indicated that both CoMFA and CoMSIA models possessed high predictive powers with $ R{^2}_{\!\!\!\rm pred} $ values of 0.818 and 0.892, $ {r^2}_{\!\!\!\rm m} $ values of 0.879 and 0.886, $ {r^2}_{\!\!\!\rm m(LOO)} $ values of 0.874 and 0.874, $ r^2_{\rm m(overall)} $ values of 0.879 and 0.885, respectively. Several key structural features of the compounds responsible for inhibitory activity were discussed in detail. Based on these structural factors, eight new compounds with quite higher predicted Src‐inhibitory activities have been designed and presented. We hope these theoretical results can offer some valuable references for the pharmaceutical molecular design as well as the action mechanism analysis. © 2012 Wiley Periodicals, Inc.     

Insight into the structural requirements of benzimidazole derivatives as interleukin‐2 inducible t‐cell kinase inhibitors by computational explorations

In the present work, a set of ligand‐ and receptor‐based 3D‐QSAR models were developed to explore the structure–activity relationship of 109 benzimidazole‐based interleukin‐2‐inducible T‐cell kinase (ITK) inhibitors. In order to reveal the requisite 3D structural features impacting the biological activities, a variety of in silico modeling approaches including the comparative molecular field analysis (CoMFA), comparative molecular similarity indices analysis (CoMSIA), docking, and molecular dynamics were applied. The results showed that the ligand‐based CoMFA model (Q² = 0.552, R²_ncv = 0.908, R²_pred = 0.787, SEE = 0.252, SEP = 0.558) and CoMSIA model (Q² = 0.579, R²_ncv = 0.914, R²_pred = 0.893, SEE = 0.240, SEP = 0.538) were superior to other models with greater predictive power. In addition, a combined analysis between the 3D contour maps and docking results showed that: (1) Compounds with bulky or hydrophobic substituents near ring D and electropositive or hydrogen acceptor groups around rings C and D could increase the activity. (2) The key amino acids impacting the receptor–ligand interactions in the binding pocket are Met438, Asp500, Lys391, and Glu439. The results obtained from this work may provide helpful guidelines in design of novel benzimidazole analogs as inhibitors of ITK. © 2013 Wiley Periodicals, Inc.     

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.     

Development of docking-based 3D QSAR models for the design of substituted quinoline derivatives as human dihydroorotate dehydrogenase (hDHODH) inhibitors

Abstract

This study has investigated docking-based 3D quantitative structure–activity relationships (QSARs) for a range of quinoline carboxylic acid derivatives by comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA). A docking study has shown that most of the compounds formed H-bonds with Arg136 and Gln47, which have already been shown to be essential for the binding of ligands at the active site of the hydroorotate dehydrogenase adenovirus (hDHODH). Bioactive conformations of all the molecules obtained from the docking study were used for the 3D QSAR study. The best CoMFA and CoMSIA models were obtained for the training set and were found to be statistically significant, with cross-validated coefficients (q² ) of 0.672 and 0.613, r² _cv of 0.635 and 0.598 and coefficients of determination (r² ) of 0.963 and 0.896, respectively. Both models were validated by a test set of 15 compounds, giving satisfactory predicted correlation coefficients (r² _pred) of 0.824 and 0.793 for the CoMFA and CoMSIA models, respectively. From the docking-based 3D QSAR study we designed 34 novel quinoline-based compounds and performed structure-based virtual screening. Finally, in silico pharmacokinetics and toxicities were predicted for 24 of the best docked molecules. The study provides valuable information for the understanding of interactions between hDHODH and the novel compounds.     

3D CoMFA,CoMSIA, topomer CoMFA and HQSAR studies on aromatic acid esters for carbonic anhydrase inhibitory activity

Molecular modelling studies [comparative molecular field analysis (CoMFA), comparative molecular similarity indices analysis (CoMSIA), topomer CoMFA and hologram quantitative structure–activity relationship (HQSAR)] have been performed on the series of 28 molecules belonging to the series of aromatic acid ester derivatives for their carbonic anhydrase inhibitory activity. The model exhibited good correlation coefficient (r²) and cross‐validated correlation coefficient (q²) for CoMFA, CoMSIA and HQSAR methods. On the basis of the findings from all these studies, a structure–activity relationship was established. Copyright © 2013 John Wiley & Sons, Ltd.     

Design,synthesis, pharmacological evaluation and in silico ADMET prediction of novel substituted benzimidazole derivatives as angiotensin II–AT1 receptor antagonists based on predictive 3D QSAR models

In this study we designed novel substituted benzimidazole derivatives and predicted their absorption, distribution, metabolism, excretion and toxicity (ADMET) properties, based on a predictive 3D QSAR study on 132 substituted benzimidazoles as AngII–AT₁ receptor antagonists. The two best predicted compounds were synthesized and evaluated for AngII–AT₁ receptor antagonism. Three different alignment tools for comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) were used. The best 3D QSAR models were obtained using the rigid body (Distill) alignment method. CoMFA and CoMSIA models were found to be statistically significant with leave-one-out correlation coefficients (q²) of 0.630 and 0.623, respectively, cross-validated coefficients (r²_cv) of 0.651 and 0.630, respectively, and conventional coefficients of determination (r²) of 0.848 and 0.843, respectively. 3D QSAR models were validated using a test set of 24 compounds, giving satisfactory predicted results (r²_pred) of 0.727 and 0.689 for the CoMFA and CoMSIA models, respectively. We have identified some key features in substituted benzimidazole derivatives, such as lipophilicity and H-bonding at the 2- and 5-positions of the benzimidazole nucleus, respectively, for AT₁ receptor antagonistic activity. We designed 20 novel substituted benzimidazole derivatives and predicted their activity. In silico ADMET properties were also predicted for these designed molecules. Finally, the compounds with best predicted activity were synthesized and evaluated for in vitro angiotensin II–AT₁ receptor antagonism.     

3D-QSAR (CoMFA,CoMSIA, HQSAR and topomer CoMFA), MD simulations and molecular docking studies on purinylpyridine derivatives as B-Raf inhibitors for the treatment of melanoma cancer

As per the World Health Organization (WHO), cancer is the second most leading cause of death after cardiovascular diseases in worldwide with around 9.88 million total new cases and 1.08 million were observed due to skin cancer in 2018. Amongst two types of skin cancer, progression of melanoma cancer is increasing day by day due to the environmental changes than non-melanoma cancer. Most of B-Raf mutation, specifically B-RafV600E, is responsible for the progression of the melanoma cancer. Here, various 3D-QSAR techniques like comparative molecular field analysis (CoMFA), comparative molecular similarity indices analysis (CoMSIA), molecular hologram QSAR (HQSAR) and topomer CoMFA were used to design novel B-Raf inhibitors by using 28 synthetic B-Raf inhibitors. Except for topomer CoMFA model, remaining models were generated by three different alignment methods in which distil-based alignment method was found best and gave prominent statistical values. After performing N-fold statistical validation, in CoMFA, q², r² and r²_pred values were found to be 0.638, 0.969 and 0.848, respectively. Similarly, q², r² and r²_pred values were found to be 0.796, 0.978 and 0.891 in CoMSIA (SHD) and 0.761, 0.973 and 0.852 in CoMSIA (SH) by N-fold statistical validation. In HQSAR analysis, statistical values were found for q² as 0.984, r² as 0.999 and r²_pred as 0.634 with 97 as best hologram length (BHL). The results of topomer CoMFA showed the q² value of 0.663 and the r² value of 0.967. Important features of purinylpyridine were identified by contour map analysis of all 3D-QSAR techniques, which could be useful to design the novel molecules as B-Raf inhibitors for the treatment of melanoma cancer.

     

The 3D Structure of the Binding Pocket of the Human Oxytocin Receptor for Benzoxazine Antagonists,Determined by Molecular Docking,Scoring Functions and 3D-QSAR Methods

Summary Molecular docking and 3D-QSAR studies were performed to determine the binding mode for a series of benzoxazine oxytocin antagonists taken from the literature. Structural hypotheses were generated by docking the most active molecule to the rigid receptor by means of AutoDock 3.05. The cluster analysis yielded seven possible binding conformations. These structures were refined by using constrained simulated annealing, and the further ligands were aligned in the refined receptor by molecular docking. A good correlation was found between the estimated ΔG_bind and the pK_i values for complex F. The Connolly-surface analysis, CoMFA and CoMSIA models q_CoMFA² = 0.653, q_CoMSA² = 0.630 and r_pred,CoMFA² = 0.852 , r_pred,CoMSIA² = 0.815) confirmed the scoring function results. The structural features of the receptor–ligand complex and the CoMFA and CoMSIA fields are in closely connected. These results suggest that receptor–ligand complex F is the most likely binding hypothesis for the studied benzoxazine analogs.     

A 3D‐QSAR Study of Celebrex‐Based Pdk1 Inhibitors Using Comfa Method

A 3D‐QSAR study of celebrex‐based compounds of PDK1 inhibitors using comparative molecular field analysis (CoMFA) was carried out. The structures of the compounds were obtained using quantum chemistry calculation. CoMFA calculations for a number of grouped subsets of compounds gave q² values of correlation in the range from 0 to 0.8. The low q² values should be mainly due to the narrow span of biological activity. Calculations for several subsets of 11–13 compounds gave high q² values, with 0.5–0.8. Factors affecting the results of the calculations are discussed. Calculated results with high q² values suggest that further chemical modifications of the compounds could lead to enhanced activity and could be an aid in the design of celebrex‐based cancer drugs.     

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.     

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.     

3D‐QSAR studies and molecular design on a novel series of pyrimidine benzimidazoles as Lck inhibitors

The development of selective lymphocyte‐specific kinase (Lck) inhibitors has attracted much attention for the research of the treatment of T‐cell mediated autoimmune and inflammatory diseases. In the present work, three‐dimensional quantitative structure–activity relationship (3D‐QSAR) analyses are performed on a novel series of 4‐amino‐6‐benzimidazole‐pyrimidines acting as Lck inhibitors. The established 3D‐QSAR models show significant statistical quality and satisfactory predictive ability, with high q² and R² values: the comparative molecular field analysis (CoMFA) model (q² = 0.802, R² = 0.991), and the comparative molecular similarity indexes analysis (CoMSIA) model (q² = 0.731, R² = 0.982). The systemic external validation indicates that both CoMFA and CoMSIA models are quite robust and possess high predictive abilities with values of 0.881 and 0.877, values of 0.897 and 0.847, values of 0.897 and 0.850, and values of 0.897 and 0.854, respectively. Several key structural features accounting for the inhibitory activities of these compounds are discussed. Based on established models and design considerations, six new compounds with significantly improved activities are theoretically designed, which still await experimental confirmation and evaluation. These theoretical results may provide a useful reference for understanding the action mechanism and designing novel potential Lck inhibitors. © 2014 Wiley Periodicals, Inc.     

3-D QSAR for intrinsic activity of 5-HT1A receptor ligands by the method of comparative molecular field analysis

The affinity of a ligand for a receptor is usually expressed in terms of the dissociation constant (K_i) of the drug-receptor complex, conveniently measured by the inhibition of radioligand binding. However, a ligand can be an antagonist, a partial agonist, or a full agonist, a property largely independent of its receptor affinity. This property can be quantitated as intrinsic activity (1A), which can range from 0 for a full antagonist to 1 for a full agonist. Although quantitative structure–activity relationship (QSAR) methods have been applied to the prediction of receptor affinity with considerable success, the prediction of IA, even qualitatively, has rarely been attempted. Because most traditional QSAR methods are limited to congeneric series, and there are often major structural differences between agonists and antagonists, this lack of success in predicting IA is understandable. To overcome this limitation, we used the method of comparative molecular field analysis (CoMFA), which, unlike traditional Hansch analysis, permits the inclusion of structurally dissimilar compounds in a single QSAR model. A structurally diverse set of 5-hydroxytryptamine_1A (5-HT_1A) receptor ligands, with literature IA data (determined by the inhibition of 5-HT sensitive forskolin-stimulated adenylate cyclase), was used to develop a 3-D QSAR model correlating intrinsic activity with molecular structure properties of 5HT_1A receptor ligands. This CoMFA model had a crossvalidated r² of 0.481, five components and final conventional r² of 0.943. The receptor model suggests that agonist and antagonist ligands can share parts of a common binding site on the receptor, with a primary agonist binding region that is also occupied by antagonists and a secondary binding site accommodating the excess bulk present in the sidechains of many antagonists and partial agonists. The CoMFA steric field graph clearly shows that agonists tend to be “flatter” (more coplanar) than antagonists, consistent with the difference between the 5-HT_1A agonist and antagonist pharmacophores proposed by Hibert and coworkers. The CoMFA electrostatic field graph suggests that, in the region surrounding the essential protonated aliphatic amino group, the positive molecular electrostatic potential may be weaker in antagonists as compared to agonists. Together, the steric and electrostatic maps suggest that in the secondary binding site region increased hydrophobic binding may enhance antagonist activity. These results demonstrate that CoMFA is capable of generating a statistically crossvalidated 3-D QSAR model that can successfully distinguish between agonist and antagonist 5-HT_1A ligands. To the best of our knowledge, this is the first time this or any other QSAR method has been successfully applied to the correlation of structure with IA rather than potency or affinity. The analysis has suggested various structural features associated with agonist and antagonist behaviors of 5-HT_1A ligands and thus should assist in the future design of drugs that act via 5-HT_1A receptors. © 1993 John Wiley & Sons, Inc.     

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.