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.     

CoMFA 3D-QSAR analysis of HIV-1 RT nonnucleoside inhibitors, TIBO derivatives based on docking conformation and alignment

HIV-1 RT is one of the key enzymes in the duplication of HIV-1. Inhibitors of HIV-1 RT are classified as nonnucleoside RT inhibitors (NNRTIs) and nucleoside analogues. NNRTIs bind in a region not associated with the active site of the enzyme. Within the NNRTI category, there is a set of inhibitors commonly referred to as TIBO inhibitors. Fifty TIBO inhibitors were used in the work to build 3-D QSAR models. The two known crystal structures of complexes are used to investigate and validate the docking protocol. The results show that the docking simulations reproduce the crystal complexes very well with RMSDs of approximately 1 A and approximately 0.6 A for 1REV and 1COU, respectively. The alignment of molecules and "active" conformation selection are the key to a successful 3D-QSAR model by CoMFA. The flexible docking (Autodock3) was used on determination of "active" conformation and molecular alignment, and CoMFA and CoMSIA were used to develop 3D-QSAR models of 50 TIBOs in the work. The 3D-QSAR models demonstrate a good ability to predict the activity of studied compounds (r2 = 0.972, 0.944, q2 = 0.704, 0.776). It is shown that the steric and electrostatic properties predicted by CoMFA contours can be related to the binding structure of the complex. The results demonstrate that the combination of ligand-based and receptor-based modeling is a powerful approach to build 3D-QSAR models.     

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.

     

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.     

3D-QSAR Study of Melittin and Amoebapore Analogues by CoMFA and CoMSIA Methods

Peptides are one of the indispensable substances in life. The use of computer aided drug design(CADD) methods to design peptides and peptiodmimetics can short the design cycle, save research funding, improve the level of whole research to a large extent and guide the discovery of new drugs. In this paper, Melittin and amoebapore three-dimensional quantitative structureactivity relationship(3D-QSAR) models were established by using comparative molecular field analysis(CoMFA) and comparative molecular similarity indices analysis(CoMSIA) method. The result shows that, the correlation coefficient(q~2) was 0.583 and non-cross-validation correlation coefficient(r~2) was 0.972 for the melittin CoMFA model. The q~2 and r~2 were 0.630 and 0.995 for the best CoMSIA model, 0.645 and 0.993 for the amoebapore CoMFA model, and 0.738 and 0.996 for the best CoMSIA model. The statistical parameters demonstrated that the CoMFA and CoMSIA models had both good predictive ability and high statistical stability, and can provide theoretical basis for designing new high activity polypeptide drugs.     

CoMFA, CoMSIA and eigenvalue analysis on dibenzodioxepinone and dibenzodioxocinone derivatives as cholesteryl ester transfer protein inhibitors

CoMFA, CoMSIA and eigenvalue analysis (EVA) were performed to study the structural features of 61 diverse dibenzodioxepinone and dibenzodioxocinone analogues to probe cholesteryl ester transfer protein (CETP) inhibitory activity. Three methods yielded statistically significant models upon assessment of cross-validation, bootstrapping, and progressive scrambling. This was further validated by an external set of 13 derivatives. Our results demonstrate that three models have a good interpolation as well as extrapolation. The hydrophobic features were confirmed to contribute significantly to inhibitor potencies, while a pre-oriented hydrogen bond provided by the hydroxyl group at the 3-position indicated a good correlation with previous SAR, and a hydrogen bond acceptor may play a crucial role in CETP inhibition. These derived models may help us to gain a deeper understanding of the binding interaction of these lactone-based compounds and aid in the design of new potent compounds against CETP.     

HIV-1逆转录酶抑制剂的三维定量构效关系

采用比较分子力场分析法(CoMFA)和比较分子相似性指数分析法(CoMSIA)对一系列非核苷类HIV-1逆转录酶抑制剂(苯磺酰基亚胺噻唑类化合物)进行了三维定量构效关系的研究,获得了高可靠性的CoMFA和CoMSIA模型,其交叉验证相关系数q2值分别为0.748和0.607.通过对CoMFA和CoMSIA模型三维等势图的分析,确定了该类化合物抗HIV-1活性的结构要求.研究结果表明,对苯磺酰基亚胺噻唑类化合物而言,在苯环的C-5位引入体积大和电负性强的基团能增加其抑制活性;苯环的C-2位的氢键给体基团对活性有利;噻唑环的R2取代基疏水性增大会降低生物活性.研究结果表明,可以指导新HIV-1逆转录酶抑制剂的设计和合成.     

Syntheses of new difluoromethylene benzoxazole and 1,2,4-oxadiazole derivatives, as potent non-nucleoside HIV-1 reverse transcriptase inhibitors

In an effort to prepare new fluorine-containing compounds, which are active against HIV, several chemical modifications of a series of benzoxazole and 1,2,4-oxadiazole-CF₂CHOHAr derivatives were carried out. The products (9-30) which all have one or two CF₂ groups were tested for activity against HIV-1; they were devoid of significant activity, one of them being cytotoxic.     

噁唑-苯甲酰胺类FtsZ受体抑制剂的3D-QSAR研究

FtsZ蛋白和细菌的繁殖过程有关,是药物设计理想的靶点。本文采用Sybyl分子模拟软件,利用比较分子场方法(CoMFA)及比较分子相似性指数的方法(CoMSIA)对已报道的33个噁唑-苯甲酰胺类FtsZ受体抑制剂进行了分析,建立了三维定量构效关系(3D-Qsar)模型。CoMFA模型的交互验证系数q2为0.619,线性回归系数r2为0.988;CoMSIA模型的q2为0.633,r2为0.936。模型具有较好的预测能力,为今后噁唑-苯甲酰胺类化合物的设计和改造提供了理论依据。     

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.     

Structural and electronic properties of new fullerene derivatives and their possible application as HIV-1 protease inhibitors

Density functional theory (DFT) calculations have been carried out at the hybrid Becke 3-Lee–Yang–Parr; B3LYP/3-21G** level of theory to study two series of hydroxy-chalca-acetic acid-(4-pyrrolidin-1-yl-phenyl) ester [C₆₀-C₂H₄N-(4-XCOCH₂OH)C₆H₄] and hydroxy-chalcoacetic acid-[2-(2-hydroxy-acetylchalcanyl)-4-pyrrolidin-1-yl-phenyl] ester[C₆₀-C₂H₄N-(3,4-XCOCH₂OH)C₆H₄]. The X atom is O, S or Se for the two series. The vibrational spectra, physical, chemical, thermodynamics and Quantitative Structure Activity Relationship (QSAR) properties of the studied molecules are calculated and discussed. We have evaluated these molecules as HIV-1 protease inhibitors based on the hydrogenation interaction between the hydroxymethylcarbonyl (HMC) groups and the two aspartic acid of the HIV-1 protease active site. Results show that some of the investigated fullerene-based derivatives can be considered promising as HIV-1 protease inhibitors.     

A 3D-QSAR Study of HIV-1 Integrase Inhibitors Using RASMS and Topomer CoMFA

Acquired Immunodeficiency Syndrome(AIDS) is a significant human health threat around the world. Therefore, the study of anti-human immunodeficiency virus(HIV) drug design has become an important task for today's society. In this paper, a three-dimensional quantitative structure-activity relationships study(3 D-QSAR) was conducted on 53 HIV-1 integrase inhibitors(IN) using random sampling analysis on molecular surface(RASMS) and Topomer comparative molecular field analysis(Topomer CoMFA). The multiple correlation coefficients of fitting, cross-validation, and external validation of two models were 0.926, 0.815 and 0.908 and 0.930, 0.726 and 0.855, respectively. The results indicated that two models obtained had both favorable estimation stability and good prediction capability. Topomer Search was used to search appropriate R groups from ZINC database, and 28 new compounds were designed thereby. The Topomer CoMFA model was subsequently used to predict the biological activity of these compounds, showing that 24 of the new compounds were more active than the template molecule. Ligands of the template molecule and new designed compounds were used for molecular docking to study the interaction of these compounds with the protein receptor. The results show that the ligands would form hydrogen-bonding interactions with the residues LEU58, THR83, GLN62, MET155, LYS119 and ALA154 of the protein receptor generally, thereby providing additional insights for the design of even more effective HIV/AIDS drugs.     

3D-QSAR and molecular modeling of HIV-1 integrase inhibitors

Three-dimensional quantitative structure-activity relationship (3D QSAR) methods were applied on a series of inhibitors of HIV-1 integrase with respect to their inhibition of 3-processing and 3-end joining steps in vitro.The training set consisted of 27 compounds belonging to the class of thiazolothiazepines. The predictive ability of each model was evaluated using test set I consisting of four thiazolothiazepines and test set II comprised of seven compounds belonging to an entirely different structural class of coumarins. Maximum Common Substructure (MCS) based method was used to align the molecules and this was compared with other known methods of alignment. Two methods of 3D QSAR: comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) were analyzed in terms of their predictive abilities. CoMSIA produced significantly better results for all correlations. The results indicate a strong correlation between the inhibitory activity of these compounds and the steric and electrostatic fields around them. CoMSIA models with considerable internal as well as external predictive ability were obtained. A poor correlation obtained with hydrophobic field indicates that the binding of thiazolothiazepines to HIV-1 integrase is mainly enthalpic in nature. Further the most active compound of the series was docked into the active site using the crystal structure of integrase. The binding site was formed by the amino acid residues 64-67, 116, 148, 151-152, 155-156, and 159. The comparison of coefficient contour maps with the steric and electrostatic properties of the receptor shows high level of compatibility.     

3D-QSAR study on heterocyclic topoisomerase II inhibitors using CoMSIA

Selective topoisomerase II (Topo II) inhibitors have interested to a great extent for the design of new antitumoral compounds in recent years. Comparative molecular similarity indices analysis (CoMSIA) was performed on a series of previously synthesized benzoxazole, benzimidazole, and oxazolo(4,5-b)pyridine derivatives as eukaryotic Topo II inhibitors. A training set of 16 heterocyclic compounds was used to establish the CoMSIA model. They were constructed and geometrically optimized using SYBYL v7.0. The predictive ability of the model was assessed using a test set of 7 compounds. The best model has demonstrated a good fit having r2 value of 0.968 and cross-validated coefficient q2 value as 0.562 including steric and hydrophobic fields. The hydrophobic interactions showed a dominant role for increasing Topo II inhibitor activity and hydrophilic substituent was found more important than hydrophobic one on the 5 or 6 position of benzazole moiety. The model obtained from the present study can be useful for the modification and/or evaluation of the development of new Topo II inhibitors as potential antitumor compounds.     

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.     

Homology model-guided 3D-QSAR studies of HIV-1 integrase inhibitors

In the present study, we report the exploration of binding modes of potent HIV-1 integrase (IN) inhibitors MK-0518 (raltegravir) and GS-9137 (elvitegravir) as well as chalcone and related amide IN inhibitors we recently synthesized and the development of 3D-QSAR models for integrase inhibition. Homology models of DNA-bound HIV-1 IN were constructed on the basis of the X-ray crystal structure of the foamy virus IN-DNA complex (PDB ID: 3L2T ) and used for docking. The binding modes of raltegravir and elvitegravir in our homology models are in accordance with those in the foamy virus structure revealing interactions important for inhibitor-IN binding. To gain further insights into the structural requirements for IN inhibition, three-dimensional quantitative structure activity relationship (3D-QSAR) studies were conducted using raltegravir, elvitegravir, and their analogs; our synthesized 3-keto salicylic acid IN inhibitor series; as well as other structurally related HIV-1 IN inhibitors. In the first part of the study with 103 compounds, atom-fit alignments, I and II, and docking-based alignment, III, were used to develop 3D-QSAR models 1, 2, and 3, respectively, each comprising comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) 3D-QSARs. This initial analysis indicated that the docking-based (structure-based) model 3 performed better than the atom-fit (ligand-based) models 1 and 2, in terms of statistical significance and robustness. Thus, the docking-based alignment was then subsequently used with an expanded data set of 296 compounds for building a more comprehensive 3D-QSAR, model 4. Model 4 afforded good q2 values of 0.70 and 0.75 for CoMFA and CoMSIA 3D-QSARs, respectively, and showed good predictive performance on an external validation test set of 59 compounds with predictive r2 values up to 0.71. The HIV IN-DNA homology model of biological relevance and the comprehensive 3D-QSAR models developed in the present study provide insights and new predictive tools for structure-based design and optimization of IN inhibitors.