QSAR and Pharmacophore Studies of Thiazolidine-4-carboxylic Acid Derivatives as Novel Influenza Neuraminidase Inhibitors Using HQSAR, Topomer CoMFA and CoMSIA

In order to understand the chemical-biological interactions governing their activities toward neuraminidase (NA), QSAR models of 28 thiazolidine-4-carboxylic acid derivatives with inhibitory influenza A virus were developed. The obtained HQSAR (hologram quantitative structure activity relationship), Topomer CoMFA and CoMSIA (comparative molecular similarity indices analysis) models were robust and had good exterior predictive capabilities. Moreover, QSAR modeling results elucidated that hydrogen bonds highly contributed to the inhibitory activity, then electrostatic and hydrophobic factors. Squared multiple correlation coefficients (R2) of HQSAR, Topomer CoMFA and CoMSIA models were 0.994, 0.978 and 0.996, respectively. Squared cross-validated correlation coefficients (Q2) of HQSAR, Topomer CoMFA and CoMSIA models were in turn 0.951, 919 and 0.820. Furthermore, squared multiple correlation coefficients for the test set (R2test) of HQSAR, CoMFA and CoMSIA models were 0.879, 0.912 and 0.953, respectively. Squared cross-validated correlation coefficients for the test set (Q2ext) of HQSAR, Topomer CoMFA and CoMSIA models were 0.867, 0.884 and 0.899, correspondingly.     

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.     

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.     

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.     

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.     

Studies on molecular mechanism between SHP2 and pyridine derivatives by 3D-QSAR,molecular docking and MD simulations

BackgroundSrc homology 2 (SH2)-containing protein tyrosine phosphatase 2 (SHP2) as a major phosphatase would affect the development of tumors by regulating several cellular processes, and is a significant potential target for cancer treatment.MethodsIn the present work, a series of pyridine derivatives possessing a wide range of inhibitory activity was employed to investigate the structural requirements by developing three dimensional quantitative structure–activity relationship (3D-QSAR) models using comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) methods. The results show that CoMFA (R²_cv = 0.646, R²_pred = 0.5587) and CoMSIA (R²_cv = 0.777, R²_pred = 0.7131) have excellent stability and predictability. The relationship between the inhibitory activity and structure of the inhibitors was analyzed by the derived contour maps. Furthermore, the QSAR models were validated by molecular docking and molecular dynamics simulations, which were also applied to reveal the potential molecular mechanism of these inhibitors.FindingsIt was found that Arg110, Asn216, Thr218, Thr252 and Pro490 play a crucial role in stabilizing the inhibitors. Additionally, MM/PBSA calculations provided the binding free energy were also conducted to explain the discrepancy of binding activities. Overall, the outcomes of this work could provide useful information and theoretical guidance for the development of novel and potent SHP2 inhibitors.     

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.     

In silico study directed towards identification of the key structural features of GyrB inhibitors targeting MTB DNA gyrase: HQSAR,CoMSIA and molecular dynamics simulations

ABSTRACT

Mycobacterium tuberculosis DNA gyrase subunit B (GyrB) has been identified as a promising target for rational drug design against fluoroquinolone drug-resistant tuberculosis. In this study, we attempted to identify the key structural feature for highly potent GyrB inhibitors through 2D-QSAR using HQSAR, 3D-QSAR using CoMSIA and molecular dynamics (MD) simulations approaches on a series of thiazole urea core derivatives. The best HQSAR and CoMSIA models based on IC₅₀ and MIC displayed the structural basis required for good activity against both GyrB enzyme and mycobacterial cell. MD simulations and binding free energy analysis using MM-GBSA and waterswap calculations revealed that the urea core of inhibitors has the strongest interaction with Asp79 via hydrogen bond interactions. In addition, cation-pi interaction and hydrophobic interactions of the R₂ substituent with Arg82 and Arg141 help to enhance the binding affinity in the GyrB ATPase binding site. Thus, the present study provides crucial structural features and a structural concept for rational design of novel DNA gyrase inhibitors with improved biological activities against both enzyme and mycobacterial cell, and with good pharmacokinetic properties and drug safety profiles.     

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.

     

Docking-based CoMFA and CoMSIA study of azaindole carboxylic acid derivatives as promising HIV-1 integrase inhibitors

Three-dimensional quantitative structure–activity relationship (3D-QSAR) studies were performed based on a series of azaindole carboxylic acid derivatives that had previously been reported as promising HIV-1 integrase inhibitors. Docking studies to explore the binding mode were performed based on the highly active molecule 36. The best docked conformation of molecule 36 was used as template for alignment. The comparative molecular field analysis (CoMFA) model (including steric and electrostatic fields) yielded the cross validation q ² = 0.655, non-cross validation r ² = 0.989 and predictive r ² _pred = 0.979. The best comparative molecular similarity indices analysis (CoMSIA) model (including steric, electrostatic, hydrophobic and hydrogen-bond acceptor fields) yielded the cross validation q ² = 0.719, non-cross validation r ² = 0.992 and predictive r ² _pred = 0.953. A series of new azaindole carboxylic acid derivatives were designed and the HIV-1 integrase inhibitory activities of these designed compounds were predicted based on the CoMFA and CoMSIA models.     

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.     

Design of novel SHP2 inhibitors using Topomer CoMFA,HQSAR analysis,and molecular docking

Structural Chemistry - The normal expression of SHP2 protein is a key factor in the production and action of cancer cells. Highly active SHP2 inhibitors could inhibit the promotion effect of SHP2...     

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.     

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.     

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

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.     

Identification of novel Plasmodium falciparum PI4KB inhibitors as potential anti-malarial drugs: Homology modeling,molecular docking and molecular dynamics simulations

The current study was set to discover selective Plasmodium falciparum phosphatidylinositol-4-OH kinase type III beta (pfPI4KB) inhibitors as potential antimalarial agents using combined structure-based and ligand-based drug discovery approach. A comparative model of pfPI4KB was first constructed and validated using molecular docking techniques. Performance of Autodock4.2 and Vina4 software in predicting the inhibitor-PI4KB binding mode and energy was assessed based on two Test Sets: Test Set I contained five ligands with resolved crystal structures with PI4KB, while Test Set II considered eleven compounds with known IC50 value towards PI4KB. The outperformance of Autodock as compared to Vina was reported, giving a correlation coefficient (R²) value of 0.87 and 0.90 for Test Set I and Test Set II, respectively. Pharmacophore-based screening was then conducted to identify drug-like molecules from ZINC database with physicochemical similarity to two potent pfPI4KB inhibitors –namely cpa and cpb. For each query inhibitor, the best 1000 hits in terms of TanimotoCombo scores were selected and subjected to molecular docking and molecular dynamics (MD) calculations. Binding energy was then estimated using molecular mechanics–generalized Born surface area (MM-GBSA) approach over 50 ns MD simulations of the inhibitor-pfPI4KB complexes. According to the calculated MM-GBSA binding energies, ZINC78988474 and ZINC20564116 were identified as potent pfPI4KB inhibitors with binding energies better than those of cpa and cpb, with ΔG_binding ≥ −34.56 kcal/mol. The inhibitor-pfPI4KB interaction and stability were examined over 50 ns MD simulation; as well the selectivity of the identified inhibitors towards pfPI4KB over PI4KB was reported.