3D-QSAR studies of farnesyltransferase inhibitors: a comparative molecular field analysis approach

3D-QSAR analysis has been performed on a series of previously synthesized benzonitrile derivatives, which were screened as farnesyltransferase inhibitors, using comparative molecular field analysis (CoMFA) with partial least-square fit to predict the steric and electrostatic molecular field interactions for the activity. The CoMFA study was carried out using a training set of 34 compounds. The predictive ability of the model developed was assessed using a test set of eight compounds (r(pred)(2) as high as 0.770). The analyzed 3D-QSAR CoMFA model has demonstrated a good fit, having r(2) value of 0.991 and cross-validated coefficient q(2) value as 0.619. The analysis of CoMFA contour maps provided insight into the possible modification of the molecules for better activity.     

3D-QSAR and molecular docking studies of 2-pyrimidinecarbonitrile derivatives as inhibitors against falcipain-3

The three dimensional-quantitative structure activity relationship (3D-QSAR) studies were performed on a series of falcipain-3 inhibitors using comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) techniques. A training set containing 42 molecules served to establish the QSAR models. The optimum CoMFA and CoMSIA models obtained for the training set were statistically significant with cross-validated correlation coefficients r(cv)(2) (q(2)) of 0.549 and 0.608, and conventional correlation coefficients (r(2)) of 0.976 and 0.932, respectively. An independent test set of 12 molecules validated the external predictive power of both models with predicted correlation coefficients (r(pred)(2)) for CoMFA and CoMSIA as 0.697 and 0.509, respectively. The docking of inhibitors into falcipain-3 active site using GOLD software revealed the vital interactions and binding conformation of the inhibitors. The CoMFA and CoMSIA field contour maps agree well with the structural characteristics of the binding pocket of falcipain-3 active site, which suggests that the information rendered by 3D-QSAR models and the docking interactions can provide guidelines for the development of improved falcipain-3 inhibitors.     

Comparative molecular field analysis and comparative molecular similarity indices analysis of human thymidine kinase 1 substrates

Thymidine kinase 1 (TK1) is a key target for antiviral and anticancer chemotherapy. Three-dimensional quantitative structure-activity relationship (3D-QSAR) using comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) techniques was applied to analyze the phosphorylation capacity of a series of 31 TK1 substrates. The optimal predictive CoMFA model with 26 molecules provided the following values: cross-validated r(2) (q(2))=0.651, non-cross-validated r(2)=0.980, standard error of estimate (s)=0.207, F=129.3. For the optimal CoMSIA model the following values were found: q(2)=0.619, r(2)=0.994, s=0.104, F=372.2. The CoMSIA model includes steric, electrostatic, and hydrogen bond donor fields. The predictive capacity of both models was successfully validated by calculating known phosphorylation rates of five TK1 substrates that were not included in the training set. Contour maps obtained from CoMFA and CoMSIA models correlated with the experimentally developed SAR.     

CoMFA and CoMSIA 3D-QSAR analysis of diaryloxy-methano-phenanthrene derivatives as anti-tubercular agents

Comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) based on three-dimensional quantitative structure-activity relationship (3D-QSAR) studies were conducted on a series (44 compounds) of diaryloxy-methano-phenanthrene derivatives as potent antitubercular agents. The best predictions were obtained with a CoMFA standard model (q (2)=0.625, r (2)=0.994) and with CoMSIA combined steric, electrostatic, and hydrophobic fields (q (2)=0.486, r (2)=0.986). Both models were validated by a test set of seven compounds and gave satisfactory predictive r (2) values of 0.999 and 0.745, respectively. CoMFA and CoMSIA contour maps were used to analyze the structural features of the ligands to account for the activity in terms of positively contributing physicochemical properties: steric, electrostatic, and hydrophobic fields. The information obtained from CoMFA and CoMSIA 3-D contour maps can be used for further design of phenanthrene-based analogs as anti-TB agents. The resulting contour maps, produced by the best CoMFA and CoMSIA models, were used to identify the structural features relevant to the biological activity in this series of analogs. Further analysis of these interaction-field contour maps also showed a high level of internal consistency. This study suggests that introduction of bulky and highly electronegative groups on the basic amino side chain along with decreasing steric bulk and electronegativity on the phenanthrene ring might be suitable for designing better antitubercular agents.     

Computational studies of COX-2 inhibitors: 3D-QSAR and docking

The 3D-QSAR (three-dimensional quantitative structure-activity relationships) studies for 88 selective COX-2 (cyclooxygenase-2) inhibitors belonging to three chemical classes (triaryl rings, diaryl cycloalkanopyrazoles, and diphenyl hydrazides) were conducted using comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA). Partial least squares analysis produced statistically significant models with q(2) values of 0.84 and 0.79 for CoMFA and CoMSIA, respectively. The binding energies calculated from flexible docking were correlated with inhibitory activities by the least-squares fit method. The three chemical classes of inhibitors showed reasonable internal predictability (r(2)=0.51, 0.49, and 0.54), but the sulfonyl-containing inhibitors demonstrated distinctively low binding energy compared to the others. The electrostatic interaction energy between the Arg513 of the COX-2 active site and sulfonyl group of the triaryl rings seemed to have the responsibility for difference in binding energy. Comparative binding energy (COMBINE) analyses gave q(2) values of 0.64, 0.63, and 0.50 for triaryl rings, diaryl cycloalkanopyrazoles, and diphenyl hydrazides, respectively. In this COMBINE model, some protein residues were highlighted as particularly important for inhibitory activity. The combination of ligand-based and structure-based models provided an improved understanding in the interaction between the three chemical classes and the COX-2.     

3D-QSAR of N-myristoyltransferase inhibiting antifungal agents by CoMFA and CoMSIA methods

A series of benzofuran antifungals was examined to determine the structural requirements of N-myristoyltransferase (Nmt) enzyme inhibition by three-dimensional quantitative structure-activity relationship (3D-QSAR) using comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) methods. Evaluation of 20 compounds (training set) served to establish the model, which was validated by evaluation of a set of 6 compounds (test set). The lowest energy conformer of the most active molecule obtained from systematic search was used as the template structure for the alignment. The best predictions were obtained with the CoMFA model from RMS fit, with r(2)(cv)=0.828, r(2)(conv)=0.989, r(2)(pred)=0.754 and with the CoMSIA model combining hydrophobic, hydrogen bond donor and hydrogen bond acceptor fields with r(2)(cv)=0.821, r(2)(conv)=0.978 and r(2)(pred)=0.747. The models obtained from the present study can be useful for the development of new Nmt inhibitors as potential antifungals. The docking studies were also carried out wherein the active and inactive molecules were docked into the active site of the recently reported Candida albicans Nmt (CaNmt) crystal structure to analyze enzyme-inhibitor interactions. The results obtained from the present 3D-QSAR and docking studies were found complimentary.     

3D-QSAR CoMFA studies on bis-coumarine analogues as urease inhibitors: a strategic design in anti-urease agents

A 3D-QSAR study has been performed on thirty (30) bis-coumarine derivatives to correlate their chemical structures with their observed urease inhibitory activity. Due to the absence of information on their active mechanism, comparative molecular field analysis (CoMFA) was used in the study. Two different properties: steric, electrostatic, assumed to cover the major contributions to ligand binding, were used to generate the 3D-QSAR model. Significant cross-validated correlation coefficients q(2) (0.558) and r(2) (0.992) for CoMFA were obtained, indicating the statistical significance of this class of compounds. The red electrostatic contour map highlighting those portion of compounds which may be interacting with nickel metal center in the active site of urease; while the blue contour map indicates positively charged groups in the ligands have improved biological activity and thus lower the IC(50)s. The steric contour map shows that bulkier substitutions at the 'R' position are detrimental to ligand receptor interaction. Actual urease inhibitory activities of this class and the predicted values were in good agreement with the experimental results. Moreover, from the contour maps, the key features vital to ligand binding have been identified, which are important for us to trace the important properties and gain insight into the potential mechanisms of intermolecular interactions between the ligand and receptor.     

3D QSAR studies on T-type calcium channel blockers using CoMFA and CoMSIA

Comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) were performed on a series of isoxazolyl compounds as a potent T-type calcium channel blockers. A set of 24 structurally similar compounds served to establish the model. Four different conformations of the most active compound were used as template structures for the alignment, three of which were obtained from Catalyst pharmacophore modeling and one by using SYBYL random search option. All CoMFA and CoMSIA models gave cross-validated r(2) (q(2)) value of more than 0.5 and conventional r(2) value of more than 0.85. The predictive ability of the models was validated by an external test set of 10 compounds, which gave satisfactory pred r(2) values ranging from 0.577 to 0.866 for all models. Best predictions were obtained with CoMFA std model of Conformer no: 3 alignment (q(2)=0.756, r(2)=0.963), giving predictive r(2) value of 0.866 for the test set. CoMFA and CoMSIA contour maps were used to analyze the structural features of the ligands accounting for the activity in terms of positively contributing physicochemical properties: steric, electrostatic, hydrophobic and hydrogen bonding fields.     

3-D-QSAR CoMFA and CoMSIA studies on tetrahydrofuroyl-L-phenylalanine derivatives as VLA-4 antagonists

Comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) were performed on tetrahydrofuroyl-L-phenylalanine derivatives as VLA-4 antagonists. The best CoMFA and CoMSIA models that were generated using atom based alignment from a training set of twenty five tetrahydrofuroyl-L-phenylalanine derivatives, are six-component models with good statistics; CoMFA: r(2)(cv)=0.366, r(2)=0.983, s=0.099, F=172.661 and PRESS=4.435; CoMSIA: r(2)(cv)=0.528, r(2)=0.995, s=0.054, F=577.87 and PRESS=3.563. Both of these 3-D-QSAR models were validated using a test set of eleven compounds, whose predicted pIC(50) values fall within one log unit of the actual pIC(50). The contour diagrams obtained for the various CoMFA and CoMSIA field contributions can be mapped back onto structural features to explain the activity trends of the molecules analysed. Based on the spatial arrangement of the various field contributions, novel molecules with improved activity can be designed.     

3D-QSAR,molecular docking,and molecular dynamic simulations for prediction of new Hsp90 inhibitors based on isoxazole scaffold

Heat shock protein 90(Hsp90), as a molecular chaperone, play a crucial role in folding and proper function of many proteins. Hsp90 inhibitors containing isoxazole scaffold are currently being used in the treatment of cancer as tumor suppressers. Here in the present studies, new compounds based on isoxazole scaffold were predicted using a combination of molecular modeling techniques including three-dimensional quantitative structure–activity relationship (3D-QSAR), molecular docking and molecular dynamic (MD) simulations. Comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) were also done. The steric and electrostatic contour map of CoMFA and CoMSIA were created. Hydrophobic, hydrogen bond donor and acceptor of CoMSIA model also were generated, and new compounds were predicted by CoMFA and CoMSIA contour maps. To investigate the binding modes of the predicted compounds in the active site of Hsp90, a molecular docking simulation was carried out. MD simulations were also conducted to evaluate the obtained results on the best predicted compound and the best reported Hsp90 inhibitors in the 3D-QSAR model. Findings indicate that the predicted ligands were stable in the active site of Hsp90.     

Exploration of a binding mode of indole amide analogues as potent histone deacetylase inhibitors and 3D-QSAR analyses

Docking simulations and three-dimensional quantitative structure-activity relationship (3D-QSAR) analyses were conducted on a series of indole amide analogues as potent histone deacetylase inhibitors. The studies include comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA). Selected ligands were docked into the active site of human HDAC1. Based on the docking results, a novel binding mode of indole amide analogues in the human HDAC1 catalytic core is presented, and enzyme/inhibitor interactions are discussed. The indole amide group is located in the open pocket, and anchored to the protein through a pair of hydrogen bonds with Asp99 O-atom and amide NH group on ligand. Based on the binding mode, predictive 3D-QSAR models were established, which had conventional r2 and cross-validated coefficient values (r(cv)2) up to 0.982 and 0.601 for CoMFA and 0.954 and 0.598 for CoMSIA, respectively. A comparison of the 3D-QSAR field contributions with the structural features of the binding site showed good correlation between the two analyses. The results of 3D-QSAR and docking studies validate each other and provided insight into the structural requirements for activity of this class of molecules as HDAC inhibitors. The CoMFA and CoMSIA PLS contour maps and MOLCAD-generated active site electrostatic, lipophilicity, and hydrogen-bonding potential surface maps, as well as the docking studies, provided good insights into inhibitor-HDAC interactions at the molecular level. Based on these results, novel molecules with improved activity can be designed.     

Two- and three-dimensional quantitative structure-activity relationships for a series of purine nucleoside phosphorylase inhibitors

Comparative molecular field analysis (CoMFA), comparative molecular similarity indices analysis, and hologram quantitative structure-activity relationship (HQSAR) studies were conducted on a series of 52 training set inhibitors of calf spleen purine nucleoside phosphorylase (PNP). Significant cross-validated correlation coefficients (CoMFA, q(2)=0.68; CoMSIA, q(2)=0.66; and HQSAR, q(2)=0.70) were obtained, indicating the potential of the models for untested compounds. The models were then used to predict the inhibitory potency of 16 test set compounds that were not included in the training set, and the predicted values were in good agreement with the experimental results. The final QSAR models along with the information gathered from 3D contour and 2D contribution maps should be useful for the design of novel inhibitors of PNP having improved potency.     

3D-QSAR, homology modeling, and molecular docking studies on spiropiperidines analogues as agonists of nociceptin/orphanin FQ receptor

The nociceptin/orphanin FQ receptor (NOP) has been implicated in a wide range of biological functions, including pain, anxiety, depression and drug abuse. Especially, its agonists have a great potential to be developed into anxiolytics. However, the crystal structure of NOP is still not available. In the present work, both structure-based and ligand-based modeling methods have been used to achieve a comprehensive understanding on 67N-substituted spiropiperidine analogues as NOP agonists. The comparative molecular-field analysis method was performed to formulate a reasonable 3D-QSAR model (cross-validated coefficient q(2)=0.819 and conventional r(2)=0.950), whose robustness and predictability were further verified by leave-eight-out, Y-randomization, and external test-set validations. The excellent performance of CoMFA to the affinity differences among these compounds was attributed to the contributions of electrostatic/hydrogen-bonding and steric/hydrophobic interactions, which was supported by the Surflex-Dock and CDOCKER molecular-docking simulations based on the 3D model of NOP built by the homology modeling method. The CoMFA contour maps and the molecular docking simulations were integrated to propose a binding mode for the spiropiperidine analogues at the binding site of NOP.     

CoMFA and CoMSIA 3D-QSAR analysis of DMDP derivatives as anti-cancer agents

Comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) based on three dimensional quantitative structure-activity relationship (3D-QSAR) studies were conducted on a series (78 compounds) of 2, 4-diamino-5-methyl-5-deazapteridine (DMDP) derivatives as potent anticancer agents. The best prediction were obtained with a CoMFA standard model (q(2) = 0.530, r(2) = 0.903) and with CoMSIA combined steric, electrostatic, hydrophobic and hydrogen bond donor fields (q(2) = 0.548, r(2) = 0.909). Both models were validated by a test set of ten compounds producing very good predictive r(2) values of 0.935 and 0.842, respectively. CoMFA and CoMSIA contour maps were then used to analyze the structural features of ligands to account for the activity in terms of positively contributing physiochemical properties such as steric, electrostatic, hydrophobic and hydrogen bond donor fields. The resulting contour maps produced by the best CoMFA and CoMSIA models were used to identify the structural features relevant to the biological activity in this series of analogs. This study suggests that the highly electropositive substituents with low steric tolerance are required at 5 position of the pteridine ring and bulky electronegatve substituents are required at the meta-position of the phenyl ring. The information obtained from CoMFA and CoMSIA 3-D contour maps can be used for the design of deazapteridine-based analogs as anticancer agents.     

Investigation on the binding mode of benzothiophene analogues as potent factor IXa (FIXa) inhibitors in thrombosis by CoMFA, docking and molecular dynamic studies

Recently, benzothiophenes attract much attention of interest due to its possible inhibitory activity targeting FIXa, a blood coagulation factor that is essential for the amplification or consolidation phase of blood coagulation. To explore this inhibitory mechanism, three-dimensional quantitative structure-activity relationship (3D-QSAR), molecular docking and molecular dynamics (MD) studies on a series of 84 benzothiophene analogues, for the first time, were performed. As a result, a highly predictive CoMFA model was developed with the q(2) = 0.52, r(2) = 0.97 and r(2)(pred) = 0.81, respectively. The CoMFA contour maps, the docking analysis, as well as the MD simulation results are all in a good agreement, proving the reliability and robustness of the model. These models and the information, we hoped, would be helpful in screening and development of novel drugs against thrombosis prior to synthesis.