Improved CoMFA modeling by optimization of settings

The possibility of improving the predictive ability of comparative molecular field analysis (CoMFA) by settings optimization has been evaluated to show that CoMFA predictive ability can be improved. Ten different CoMFA settings are evaluated, producing a total of 6120 models. This method has been applied to nine different data sets, including the widely used benchmark steroid data set, as well as eight other data sets proposed as QSAR benchmarking data sets by Sutherland et al. (J. Med. Chem. 2004, 47, 5541-5554). All data sets have been studied using training and test sets to allow for both internal (q(2)) and external (r(2)(pred)) predictive ability assessment. CoMFA settings optimization was successful in developing models with improved q(2) and r(2)(pred) as compared to default CoMFA modeling. Optimized CoMFA is compared with comparative molecular similarity indices analysis (CoMSIA) and holographic quantitative structure-activity relationship (HQSAR) models and found to consistently produce models with improved or equivalent q(2) and r(2)(pred). The ability of settings optimization to improve model predictive ability has been validated using both internal and external predictions, and the risk of chance correlation has been evaluated using response variable randomization tests.     

Three-dimensional quantitative structure-activity relationships of steroid aromatase inhibitors

Summary Inhibition of aromatase, a cytochrome P450 that converts androgens to estrogens, is relevant in the therapeutic control of breast cancer. We investigate this inhibition using a three-dimensional quantitative structure-activity relationship (3D QSAR) method known as Comparative Molecular Field Analysis, CoMFA [Cramer III, R.D. et al., J. Am. Chem. Soc., 110 (1988) 5959]. We analyzed the data for 50 steroid inhibitors [Numazawa, M. et al., J. Med. Chem., 37 (1994) 2198, and references cited therein] assayed against androstenedione on human placental microsomes. An initial CoMFA resulted in a three-component model for log(1/K_i), with an explained variance r² of 0.885, and a cross-validated q² of 0.673. Chemometric studies were performed using GOLPE [Baroni, M. et al., Quant. Struct.-Act. Relatsh., 12 (1993) 9]. The CoMFA/GOLPE model is discussed in terms of robustness, predictivity, explanatory power and simplicity. After randomized exclusion of 25 or 10 compounds (repeated 25 times), the q² for one component was 0.62 and 0.61, respectively, while r² was 0.674. We demonstrate that the predictive r² based on the mean activity (Y_m) of the training set is misleading, while the test set Y_m-based predictive r² index gives a more accurate estimate of external predictivity. Using CoMFA, the observed differences in aromatase inhibition among C6-substituted steroids are rationalized at the atomic level. The CoMFA fields are consistent with known, potent inhibitors of aromatase, not included in the model. When positioned in the same alignment, these compounds have distinct features that overlap with the steric and electrostatic fields obtained in the CoMFA model. The presence of two hydrophobic binding pockets near the aromatase active site is discussed: a steric bulk tolerant one, common for C4, C6-alpha and C7-alpha substitutents, and a smaller one at the C6-beta region.     

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 and molecular docking studies of benzimidazole derivatives as hepatitis C virus NS5B polymerase inhibitors

The urgent need for novel HCV antiviral agents has provided an impetus for understanding the structural requisites of NS5B polymerase inhibitors at the molecular level. Toward this objective, comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) of 67 HCV NS5B polymerase inhibitors were performed using two methods. First, ligand-based 3D QSAR studies were performed based on the lowest energy conformations employing the atom fit alignment method. Second, receptor-based 3D QSAR models were derived from the predicted binding conformations obtained by docking all NS5B inhibitors at the allosteric binding site of NS5B (PDB ID: 2dxs). Results generated from the ligand-based model were found superior (r2cv values of 0.630 for CoMFA and 0.668 for CoMSIA) to those obtained by the receptor-based model (r2cv values of 0.536 and 0.561 for CoMFA and CoMSIA, respectively). The predictive ability of the models was validated using a structurally diversified test set of 22 compounds that had not been included in a preliminary training set of 45 compounds. The predictive r2 values for the ligand-based CoMFA and CoMSIA models were 0.734 and 0.800, respectively, while the corresponding predictive r2 values for the receptor-based CoMFA and CoMSIA models were 0.538 and 0.639, respectively. The greater potency of the tryptophan derivatives over that of the tyrosine derivatives was interpreted based on CoMFA steric and electrostatic contour maps. The CoMSIA results revealed that for a NS5B inhibitor to have appreciable inhibitory activity it requires hydrogen bond donor and acceptor groups at the 5-position of the indole ring and an R substituent at the chiral carbon, respectively. Interpretation of the CoMFA and CoMSIA contour maps in context of the topology of the allosteric binding site of NS5B provided insight into NS5B-inhibitor interactions. Taken together, the present 3D QSAR models were found to accurately predict the HCV NS5B polymerase inhibitory activity of structurally diverse test set compounds and to yield reliable clues for further optimization of the benzimidazole derivatives in the data set.     

Exploration of a binding mode of benzothiazol-2-yl acetonitrile pyrimidine core based derivatives as potent c-Jun N-terminal kinase-3 inhibitors and 3D-QSAR analyses

C-Jun N-terminal kinase (JNK) is a therapeutic target for inhibitors which may provide clinical benefit in the pathogenesis of rheumatoid arthritis (RA) as well as in various apoptosis-related disorders. The benzothiazol-2-yl acetonitrile derivatives, recently reported by Pascale et al. (J. Med. Chem. 2005, 48, 4596-4607), are the first generation JNK inhibitors of this class. To understand inhibitory mechanisms and elucidate pharmacophoric properties of these derivatives molecular docking and 3D-QSAR studies were performed on a set of 44 compounds. Ligand Fit module of Cerius2 (4.9) was employed to locate the binding orientations of all the compounds within the JNK-3 ATP binding site. A good correlation (r2=0.810) between the calculated binding free energies (-PMF score) and the experimental inhibitory activities suggests that the identified binding conformations of these potential inhibitors are reliable. Based on the binding conformations, robust and highly predictive 3D-QSAR models were developed with conventional r2 0.886 and 0.802, full cross-validation r2 0.980 and 0.788, and predictive r2 0.965 and 0.968 for MFA and MSA, respectively. The interaction mode was demonstrated taking into consideration inhibitor conformation, hydrogen bonding, and electrostatic interaction. The 3D-QSAR model built in this study will provide clear guidelines for a novel inhibitor design based on the benzothiazole derivatives against JNK-3 for the treatment of inflammatory disorders.     

Mapping the binding site of a large set of quinazoline type EGF-R inhibitors using molecular field analyses and molecular docking studies

In the current work, three-dimensional QSAR studies for one large set of quinazoline type epidermal growth factor receptor (EGF-R) inhibitors were conducted using two types of molecular field analysis techniques: comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA). These compounds belonging to six different structural classes were randomly divided into a training set of 122 compounds and a test set of 13 compounds. The statistical results showed that the 3D-QSAR models derived from CoMFA were superior to those generated from CoMSIA. The most optimal CoMFA model after region focusing bears significant cross-validated r(2)(cv) of 0.60 and conventional r(2) of 0.92. The predictive power of the best CoMFA model was further validated by the accurate estimation to these compounds in the external test set, and the mean agreement of experimental and predicted log(IC(50)) values of the inhibitors is 0.6 log unit. Separate CoMFA models were conducted to evaluate the influence of different partial charges (Gasteiger-Marsili, Gasteiger-Hückel, MMFF94, ESP-AM1, and MPA-AM1) on the statistical quality of the models. The resulting CoMFA field map provides information on the geometry of the binding site cavity and the relative weights of various properties in different site pockets for each of the substrates considered. Moreover, in the current work, we applied MD simulations combined with MM/PBSA (Molecular mechanics/Possion-Boltzmann Surface Area) to determine the correct binding mode of the best inhibitor for which no ligand-protein crystal structure was present. To proceed, we define the following procedure: three hundred picosecond molecular dynamics simulations were first performed for the four binding modes suggested by DOCK 4.0 and manual docking, and then MM/PBSA was carried out for the collected snapshots. The most favorable binding mode identified by MM/PBSA has a binding free energy about 10 kcal/mol more favorable than the second best one. The most favorable binding mode identified by MM/PBSA can give satisfactory explanation of the SAR data of the studied molecules and is in good agreement with the contour maps of CoMFA. The most favorable binding mode suggests that with the quinazoline-based inhibitor, the N3 atom is hydrogen-bonded to a water molecule which, in turn, interacts with Thr 766, not Thr 830 as proposed by Wissner et al. (J. Med. Chem. 2000, 43, 3244). The predicted complex structure of quinazoline type inhibitor with EGF-R as well as the pharmacophore mapping from CoMFA can interpret the structure activities of the inhibitors well and afford us important information for structure-based drug design.     

3D-QSAR Analysis of a Series of Dihydroquinolizinone Derivatives as a Hepatitis B Virus Expression Inhibitor

In this study, we explored a three-dimensional quantitative structure-activity relationship(3D-QSAR) model of 63 HBV viral gene expression inhibitors containing dihydroquinolizinones. Two high predictive QSAR models have been built, including comparative molecular field analysis(CoMFA) and comparative molecular similarity indices analysis(CoMSIA). The internal validation parameter(CoMFA, q~2 = 0.701, r~2 = 0.999; CoMSIA, q~2 = 0.721, r~2 = 0.998) and external validation parameter(CoMFA, r~2_(pred = 0.999); CoMSIA, r~2_(pred = 0.999)) indicated that the models have good predictive abilities and significant statistical reliability. We designed several molecules with potentially higher predicted activity on the basis of the result of the models. This work might provide useful information to design novel HBV viral gene expression inhibitors.     

A comparative study of quantitative structure-activity relationship methods based on gallic acid derivatives

By using hologram quantitative structure-activity relationship (HQSAR) and comparative molecular field analysis (CoMFA) methods, the relationships between the structures of 49 gallic acid derivatives and their analgesic activity have been investigated to yield statistically reliable models with considerable predictive power. The best HQSAR model was generated using atoms, bond and connectivity as fragment distinction parameters and fragment size 5-7 from a hologram length of 307 with 3 components. High conventional r2 (r2 = 0.825) and cross-validation r2 (r2(cv) = 0.726) values were obtained. CoMFA analyses varying lattice size and location, grid spacing, probe charges and using, Tripos standard and Indicator force field were performed. The best model was developed with 4 components using sp3-hybridized carbon atom with +1.0 charge as probe, grid spacing (2 A), lattice offset (1.0, 3.0, -2.5). The CoMFA model showed a conventional correlation coefficient r2 of 0.889 and across-validation r2(cv) equals to 0.633. The robustness and predictive ability of the HQSAR and CoMFA models have been validated by means of an external test set. The results indicate that both models possess high statistical quality in the prediction of analgesic potency of novel gallic acid analogs.     

3D QSAR studies on protein tyrosine phosphatase 1B inhibitors: comparison of the quality and predictivity among 3D QSAR models obtained from different conformer-based alignments

A set of 65 flexible peptidomimetic competitive inhibitors (52 in the training set and 13 in the test set) of protein tyrosine phosphatase 1B (PTP1B) has been used to compare the quality and predictive power of 3D quantitative structure-activity relationship (QSAR) comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) models for the three most commonly used conformer-based alignments, namely, cocrystallized conformer-based alignment (CCBA), docked conformer-based alignment (DCBA), and global minima energy conformer-based alignment (GMCBA). These three conformers of 5-[(2S)-2-({(2S)-2-[(tert-butoxycarbonyl)amino]-3-phenylpropanoyl}amino)3-oxo-3-pentylamino)propyl]-2-(carboxymethoxy)benzoic acid (compound number 66) were obtained from the X-ray structure of its cocrystallized complex with PTP1B (PDB ID: 1JF7), its docking studies, and its global minima by simulated annealing. Among the 3D QSAR models developed using the above three alignments, the CCBA provided the optimal predictive CoMFA model for the training set with cross-validated r2 (q2)=0.708, non-cross-validated r2=0.902, standard error of estimate (s)=0.165, and F=202.553 and the optimal CoMSIA model with q2=0.440, r2=0.799, s=0.192, and F=117.782. These models also showed the best test set prediction for the 13 compounds with predictive r2 values of 0.706 and 0.683, respectively. Though the QSAR models derived using the other two alignments also produced statistically acceptable models in the order DCBA>GMCBA in terms of the values of q2, r2, and predictive r2, they were inferior to the corresponding models derived using CCBA. Thus, the order of preference for the alignment selection for 3D QSAR model development may be CCBA>DCBA>GMCBA, and the information obtained from the CoMFA and CoMSIA contour maps may be useful in designing specific PTP1B inhibitors.     

GA strategy for variable selection in QSAR studies: application of GA-based region selection to a 3D-QSAR study of acetylcholinesterase inhibitors

Comparative molecular field analysis (CoMFA) with partial least squares (PLS) is one of the most frequently used tools in three-dimensional quantitative structure-activity relationships (3D-QSAR) studies. Although many successful CoMFA applications have proved the value of this approach, there are some problems in its proper application. Especially, the inability of PLS to handle the low signal-to-noise ratio (sample-to-variable ratio) has attracted much attention from QSAR researchers as an exciting research target, and several variable selection methods have been proposed. More recently, we have developed a novel variable selection method for CoMFA modeling (GARGS: genetic algorithm-based region selection), and its utility has been demonstrated in the previous paper (Kimura, T., et al. J. Chem. Inf. Comput. Sci. 1998, 38, 276-282). The purpose of this study is to evaluate whether GARGS can pinpoint known molecular interactions in 3D space. We have used a published set of acetylcholinesterase (AChE) inhibitors as a test example. By applying GARGS to a data set of AChE inhibitors, several improved models with high internal prediction and low number of field variables were obtained. External validation was performed to select a final model among them. The coefficient contour maps of the final GARGS model were compared with the properties of the active site in AChE and the consistency between them was evaluated.     

Evaluation of designed ligands by a multiple screening method: Application to glycogen phosphorylase inhibitors constructed with a variety of approaches

Glycogen phosphorylase (GP) is an important enzyme that regulates blood glucose level and a key therapeutic target for the treatment of type II diabetes. In this study, a number of potential GP inhibitors are designed with a variety of computational approaches. They include the applications of MCSS, LUDI and CoMFA to identify additional fragments that can be attached to existing lead molecules; the use of 2D and 3D similarity-based QSAR models (HQSAR and SMGNN) and of the LUDI program to identify novel molecules that may bind to the glucose binding site. The designed ligands are evaluated by a multiple screening method, which is a combination of commercial and in-house ligand-receptor binding affinity prediction programs used in a previous study (So and Karplus, J. Comp.-Aid. Mol. Des., 13 (1999), 243–258). Each method is used at an appropriate point in the screening, as determined by both the accuracy of the calculations and the computational cost. A comparison of the strengths and weaknesses of the ligand design approaches is made.     

Molecular Modeling and Design of Arylthioindole Derivatives as Tubulin Inhibitors

Three-dimensional quantitative structure activity relationship （3D-QSAR） and docking studies of a series of arylthioindole derivatives as tubulin inhibitors against human breast cancer cell line MCF-7 have been carried out. An optimal 3D-QSAR model from the comparative molecular field analysis （CoMFA） for training set with significant statistical quality （R2=0.898） and predictive ability （q2=0.654） was established. The same model was further applied to predict pIC50 values of the compounds in test set, and the resulting predictive correlation coefficient R2（pred） reaches 0.816, further showing that this CoMFA model has high predictive ability. Moreover, the appropriate binding orientations and conformations of these compounds interacting with tubulin are located by docking study, and it is very interesting to find the consistency between the CoMFA field distribution and the 3D topology structure of active site of tubulin. Based on CoMFA along with docking results, some important factors improving the activities of these compounds were discussed in detail and were summarized as follows： the substituents R3-R5 （on the phenyl ring） with higher electronegativity, the substituent R6 with higher eleetropositivity and bigger bulk, the substituent R7 with smaller bulk, and so on. In addition, five new compounds with higher activities have been designed. Such results can offer useful theoretical references for experimental works.     

Structural investigations of anthranilimide derivatives by CoMFA and CoMSIA 3D-QSAR studies reveal novel insight into their structures toward glycogen phosphorylase inhibition

In the present work, three-dimensional quantitative structure–activity relationship (3-D QSAR) studies on a set of 70 anthranilimide compounds has been performed using docking-based as well as substructure-based molecular alignments. This resulted in the selection of more statistically relevant substructure-based alignment for further studies. Further, molecular models with good predictive power were derived using CoMFA (r ^2?=?0.997; Q ^2?=?0.578) and CoMSIA (r ^2?=?0.976; Q ^2?=?0.506), for predicting the biological activity of new compounds. The so-developed contour plots identified several key features of the compounds explaining wide activity ranges. Based on the information derived from the CoMFA contour maps, novel leads were proposed which showed better predicted activity with respect to the already reported systems. Thus, the present study not only offers a highly significant predictive QSAR model for anthranilimide derivatives as glycogen phosphorylase (GP) inhibitors which can eventually assist and complement the rational drug-design attempts, but also proposes a highly predictive pharmacophore model as a guide for further development of selective and more potent GP inhibitors as anti-diabetic agents.     

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.     

磺酰脲类乙酰羟基酸合成酶抑制剂Topomer CoMFA研究

磺酰脲类除草剂是一类高选择性、广谱、低毒的化合物,在世界范围内得到了广泛的应用。本文采用易位体-比较分子力场法（Topomer CoMFA）对75个磺酰脲类化合物与植物源野生型拟南芥AHAS酶的离体相互作用进行了三维定量构效关系研究,快速准确地构建了Topomer CoMFA模型,该模型具有较强的预测能力（交叉验证相关系数q2为0.890,非交叉验证相关系数r2为0.967）。此模型对测试集的10个化合物的pKi值进行预测,其预测值与实际值一致。     

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.     

A ligand-based comparative molecular field analysis (CoMFA) and homology model based molecular docking studies on 3′, 4′-dihydroxyflavones as rat 5-lipoxygenase inhibitors: Design of new inhibitors

In this study, ligand based comparative molecular field analysis (CoMFA) with five principal components was performed on class of 3′, 4′-dihydroxyflavone derivatives for potent rat 5-LOX inhibitors. The percentage contributions in building of CoMFA model were 91.36% for steric field and 8.6% for electrostatic field. R² values for training and test sets were found to be 0.9320 and 0.8259, respectively. In case of LOO, LTO and LMO cross validation test, q² values were 0.6587, 0.6479 and 0.5547, respectively. These results indicate that the model has high statistical reliability and good predictive power. The extracted contour maps were used to identify the important regions where the modification was necessary to design a new molecule with improved activity. The study has developed a homology model for rat 5-LOX and recognized the key residues at the binding site. Docking of most active molecule to the binding site of 5-LOX confirmed the stability and rationality of CoMFA model. Based on molecular docking results and CoMFA contour plots, new inhibitors with higher activity with respect to the most active compound in data set were designed.