Combined 3D-QSAR Modeling and Molecular Docking Studies on Pyrrole-Indolin-2-ones as Aurora A Kinase Inhibitors

Aurora kinases have emerged as attractive targets for the design of anticancer drugs. 3D-QSAR (comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA)) and Surflex-docking studies were performed on a series of pyrrole-indoline-2-ones as Aurora A inhibitors. The CoMFA and CoMSIA models using 25 inhibitors in the training set gave r(2) (cv) values of 0.726 and 0.566, and r(2) values of 0.972 and 0.984, respectively. The adapted alignment method with the suitable parameters resulted in reliable models. The contour maps produced by the CoMFA and CoMSIA models were employed to rationalize the key structural requirements responsible for the activity. Surflex-docking studies revealed that the sulfo group, secondary amine group on indolin-2-one, and carbonyl of 6,7-dihydro-1H-indol-4(5H)-one groups were significant for binding to the receptor, and some essential features were also identified. Based on the 3D-QSAR and docking results, a set of new molecules with high predicted activities were designed.     

3D-QSAR and Molecular Docking Studies on Fused Pyrazoles as p38α Mitogen-Activated Protein Kinase Inhibitors

The p38α mitogen-activated protein kinase (MAPK) has become an attractive target for the treatment of many diseases such as rheumatoid arthritis, inflammatory bowel disease and Crohn's disease. In this paper, 3D-QSAR and molecular docking studies were performed on 59 p38α MAPK inhibitors. Comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) were applied to determine the structural requirements for potency in inhibiting p38α MAPK. The resulting model of CoMFA and CoMSIA exhibited good r(2) (cv) values of 0.725 and 0.609, and r(2) values of 0.961 and 0.905, respectively. Molecular docking was used to explore the binding mode between the inhibitors and p38α MAPK. We have accordingly designed a series of novel p38α MAPK inhibitors by utilizing the structure-activity relationship (SAR) results revealed in the present study, which were predicted with excellent potencies in the developed models. The results provided a useful guide to design new compounds for p38α MAPK inhibitors.     

Combined pharmacophore modeling, docking, and 3D-QSAR studies of PLK1 inhibitors

Polo-like kinase 1, an important enzyme with diverse biological actions in cell mitosis, is a promising target for developing novel anticancer drugs. A combined molecular docking, structure-based pharmacophore modeling and three-dimensional quantitative structure-activity relationship (3D-QSAR) study was performed on a set of 4,5-dihydro-1H-pyrazolo[4,3-h]quinazoline derivatives as PLK1 inhibitors. The common substructure, molecular docking and pharmacophore-based alignment were used to develop different 3D-QSAR models. The comparative molecular field analysis (CoMFA) and comparative molecule similarity indices analysis (CoMSIA) models gave statistically significant results. These models showed good q(2) and r(2) (pred) values and revealed a good response to test set validation. All of the structural insights obtained from the 3D-QSAR contour maps are consistent with the available crystal structure of PLK1. The contour maps obtained from the 3D-QSAR models in combination with the structure based pharmacophore model help to better interpret the structure-activity relationship. These satisfactory results may aid the design of novel PLK1 inhibitors. This is the first report on 3D-QSAR study of PLK1 inhibitors.     

Comparative Molecular Field Analysis (CoMFA) and Comparative Molecular Similarity Indices Analysis (CoMSIA) Studies on α(1A)-Adrenergic Receptor Antagonists Based on Pharmacophore Molecular Alignment

The α(1A)-adrenergic receptor (α(1A)-AR) antagonist is useful in treating benign prostatic hyperplasia, lower urinary tract symptoms, and cardiac arrhythmia. Three-dimensional quantitative structure-activity relationship (3D-QSAR) studies were performed on a set of α(1A)-AR antagonists of N-aryl and N-nitrogen class. Statistically significant models constructed from comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) were established based on a training set of 32 ligands using pharmacophore-based molecular alignment. The leave-oneout cross-validation correlation coefficients were q(2) (CoMFA) = 0.840 and q(2) (CoMSIA) = 0.840. The high correlation between the cross-validated/predicted and experimental activities of a test set of 12 ligands revealed that the CoMFA and CoMSIA models were robust (r(2) (pred) (/CoMFA) = 0.694; r(2) (pred) (/CoMSIA) = 0.671). The generated models suggested that electrostatic, hydrophobic, and hydrogen bonding interactions play important roles between ligands and receptors in the active site. Our study serves as a guide for further experimental investigations on the synthesis of new compounds. Structural modifications based on the present 3D-QSAR results may lead to the discovery of other α(1A)-AR antagonists.     

Inhibitory mode of N-phenyl-4-pyrazolo[1,5-b] pyridazin-3-ylpyrimidin-2-amine series derivatives against GSK-3: molecular docking and 3D-QSAR analyses

Glycogen synthase kinase 3 (GSK-3) inhibition is an important research topic because of its wide range of associated health implications. The interaction mode of a series of N-phenyl-4-pyrazolo[1,5-b]pyridazin-3-ylpyrimidin-2-amine compounds with human GSK-3 has been studied using molecular docking and 3D-QSAR approaches. In the 3D-QSAR studies, the molecular alignment and conformation determination are so important that they affect the success of a model. Flexible docking (AutoDock3.0.5) was used for the determination of 'active' conformation and molecular alignment. Comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) were used to develop 3D-QSAR models of 80 N-phenyl-4-pyrazolo[1,5-b]pyridazin-3-ylpyrimidin-2-amine compounds. The r(2) values were 0.870 and 0.861 for CoMFA and CoMSIA models, respectively. The predictive ability of these models was validated by 10 compounds of the test set. Mapping these models back to the topology of the active site of GSK-3 led to a better understanding of the vital N-phenyl-4-pyrazolo[1,5-b]pyridazin-3-ylpyrimidin-2-amines-GSK-3 interactions. The results demonstrate that combination of ligand-based and receptor-based modeling is a powerful approach to build 3D-QSAR models. The interaction mode from this study may be helpful for the design of a novel inhibitor and guide the selection of candidate sites for further experimental studies on site-directed mutagenesis.     

Pyrrolo[3,2-d]pyrimidine Derivatives as Type II Kinase Insert Domain Receptor (KDR) Inhibitors: CoMFA and CoMSIA Studies

Kinase insert domain receptor (KDR) inhibitors have been proved to be very effective anticancer agents. Molecular docking, 3D-QSAR methods, CoMFA and CoMSIA were performed on pyrrolo[3,2-d]pyrimidine derivatives as non-ATP competitive KDR inhibitors (type II). The bioactive conformation was explored by docking one potent compound 20 into the active site of KDR in its DFG-out inactive conformation. The constructed CoMFA and CoMSIA models produced statistically significant results with the cross-validated correlation coefficients q(2) of 0.542 and 0.552, non-cross-validated correlation coefficients r(2) of 0.912 and 0.955, and predicted correction coefficients r(2) (pred) of 0.913 and 0.897, respectively. These results ensure the CoMFA and CoMSIA models as a tool to guide the design of a series of new potent KDR inhibitors.     

Combined 3D-QSAR and Docking Modelling Study on Indolocarbazole Series Compounds as Tie-2 Inhibitors

Tie-2, a kind of endothelial cell tyrosine kinase receptor, is required for embryonic blood vessel development and tumor angiogenesis. Several compounds that showed potent activity toward this attractive anticancer drug target in the assay have been reported. In order to investigate the structure-activity correlation of indolocarbazole series compounds and modify them to improve their selectivity and activity, 3D-QSAR models were built using CoMFA and CoMSIA methods and molecular docking was used to check the results. Based on the common sketch align, two good QSAR models with high predictabilities (CoMFA model: q(2) = 0.823, r(2) = 0.979; CoMSIA model: q(2) = 0.804, r(2) = 0.967) were obtained and the contour maps obtained from both models were applied to identify the influence on the biological activity. Molecular docking was then used to confirm the results. Combined with the molecular docking results, the detail binding mode between the ligands and Tie-2 was elucidated, which enabled us to interpret the structure-activity relationship. These satisf actory results not only offered help to comprehend the action mechanism of indolocarbazole series compounds, but also provide new information for the design of new potent inhibitors.     

Structural Determination of Three Different Series of Compounds as Hsp90 Inhibitors Using 3D-QSAR Modeling, Molecular Docking and Molecular Dynamics Methods

Hsp90 is involved in correcting, folding, maturation and activation of a diverse array of client proteins; it has also been implicated in the treatment of cancer in recent years. In this work, comparative molecular field analysis (CoMFA), comparative molecular similarity indices analysis (CoMSIA), molecular docking and molecular dynamics were performed on three different series of Hsp90 inhibitors to build 3D-QSAR models, which were based on the ligand-based or receptor-based methods. The optimum 3D-QSAR models exhibited reasonable statistical characteristics with averaging internal q(2) > 0.60 and external r(2) (pred) > 0.66 for Benzamide tetrahydro-4H-carbazol-4-one analogs (BT), AT13387 derivatives (AT) and Dihydroxylphenyl amides (DA). The results revealed that steric effects contributed the most to the BT model, whereas H-bonding was more important to AT, and electrostatic, hydrophobic, H-bond donor almost contributed equally to the DA model. The docking analysis showed that Asp93, Tyr139 and Thr184 in Hsp90 are important for the three series of inhibitors. Molecular dynamics simulation (MD) further indicated that the conformation derived from docking is basically consistent with the average structure extracted from MD simulation. These results not only lead to a better understanding of interactions between these inhibitors and Hsp90 receptor but also provide useful information for the design of new inhibitors with a specific activity.     

Binding Studies and Lead Generation of Pteridin-7(8H)-one Derivatives Targeting FLT3

Ligand modification by substituting chemical groups within the binding pocket is a popular strategy for kinase drug development. In this study, a series of pteridin-7(8H)-one derivatives targeting wild-type FMS-like tyrosine kinase-3 (FLT3) and its D835Y mutant (FL3_D835Y) were studied using a combination of molecular modeling techniques, such as docking, molecular dynamics (MD), binding energy calculation, and three-dimensional quantitative structure-activity relationship (3D-QSAR) studies. We determined the protein–ligand binding affinity by employing molecular mechanics Poisson–Boltzmann/generalized Born surface area (MM-PB/GBSA), fast pulling ligand (FPL) simulation, linear interaction energy (LIE), umbrella sampling (US), and free energy perturbation (FEP) scoring functions. The structure–activity relationship (SAR) study was conducted using comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA), and the results were emphasized as a SAR scheme. In both the CoMFA and CoMSIA models, satisfactory correlation statistics were obtained between the observed and predicted inhibitory activity. The MD and SAR models were co-utilized to design several new compounds, and their inhibitory activities were anticipated using the CoMSIA model. The designed compounds with higher predicted pIC₅₀ values than the most active compound were carried out for binding free energy evaluation to wild-type and mutant receptors using MM-PB/GBSA, LIE, and FEP methods.     

A Combined Pharmacophore Modeling, 3D QSAR and Virtual Screening Studies on Imidazopyridines as B-Raf Inhibitors

B-Raf kinase is an important target in treatment of cancers. In order to design and find potent B-Raf inhibitors (BRIs), 3D pharmacophore models were created using the Genetic Algorithm with Linear Assignment of Hypermolecular Alignment of Database (GALAHAD). The best pharmacophore model obtained which was used in effective alignment of the data set contains two acceptor atoms, three donor atoms and three hydrophobes. In succession, comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) were performed on 39 imidazopyridine BRIs to build three dimensional quantitative structure-activity relationship (3D QSAR) models based on both pharmacophore and docking alignments. The CoMSIA model based on the pharmacophore alignment shows the best result (q² = 0.621, r²_pred = 0.885). This 3D QSAR approach provides significant insights that are useful for designing potent BRIs. In addition, the obtained best pharmacophore model was used for virtual screening against the NCI2000 database. The hit compounds were further filtered with molecular docking, and their biological activities were predicted using the CoMSIA model, and three potential BRIs with new skeletons were obtained.     

Studies on [5,6]-Fused Bicyclic Scaffolds Derivatives as Potent Dual B-RafV600E/KDR Inhibitors Using Docking and 3D-QSAR Approaches

Research and development of multi-target inhibitors has attracted increasing attention as anticancer therapeutics. B-Raf^V600E synergistically works with vascular endothelial growth factor receptor 2 (KDR) to promote the occurrence and progression of cancers, and the development of dual-target drugs simultaneously against these two kinds of kinase may offer a better treatment advantage. In this paper, docking and three-dimensional quantitative structure activity relationship (3D-QSAR) studies were performed on a series of dual B-Raf/KDR inhibitors with a novel hinge-binding group, [5,6]-fused bicyclic scaffold. Docking studies revealed optimal binding conformations of these compounds interacting with both B-Raf and KDR. Based on these conformations, comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) 3D-QSAR models were constructed, and the best CoMFA (q² = 0.542, r² = 0.989 for B-Raf; q² = 0.768, r² = 0.991 for KDR) and CoMSIA models (q² = 0.519, r² = 0.992 for B-Raf; q² = 0.849, r² = 0.993 for KDR) were generated. Further external validations confirmed their predictability, yielding satisfactory correlation coefficients (r²_pred = 0.764 (CoMFA), r²_pred = 0.841 (CoMSIA) for B-Raf, r²_pred = 0.912 (CoMFA), r²_pred = 0.846 (CoMSIA) for KDR, respectively). Through graphical analysis and comparison on docking results and 3D-QSAR contour maps, key amino acids that affect the ligand-receptor interactions were identified and structural features influencing the activities were discussed. New potent derivatives were designed, and subjected to preliminary pharmacological evaluation. The study may offer useful references for the modification and development of novel dual B-Raf/KDR inhibitors.     

Pharmacophore and Molecular Docking Guided 3D-QSAR Study of Bacterial Enoyl-ACP Reductase (FabI) Inhibitors

Enoyl acyl carrier protein (ACP) reductase (FabI) is a potential target for the development of antibacterial agents. Three-dimensional quantitative structure-activity relationships (3D-QSAR) for substituted formamides series of FabI inhibitors were investigated using comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) techniques. Pharmacophore and molecular docking methods were used for construction of the molecular alignments. A training set of 36 compounds was performed to create the 3D-QSAR models and their external predictivity was proven using a test set of 11 compounds. Graphical interpretation of the results revealed important structural features of the formamides related to the active site of FabI. The results may be exploited for further optimization of the design of new potent FabI inhibitors.     

Three Dimensional Quantitative Structure-Activity Relationship of 5H-Pyrido[4,3-b]indol-4-carboxamide JAK2 Inhibitors

Janus kinase 2 (JAK2) is an intracellular nonreceptor tyrosine kinase that belongs to the JAK family of kinases, which play an important role in survival, proliferation, and differentiation of a variety of cells. JAK2 inhibitors are potential drugs for the treatment of myeloproliferative neoplasms. The three dimensional quantitative structure-activity relationships have been studied on a series of JAK2 inhibitors by comparative molecular field analysis (CoMFA), and comparative molecular similarity indices analysis (CoMSIA). The CoMFA model had a cross-validated coefficient q² of 0.633, and the relation non-cross-validated coefficient r² of 0.976. The F value is 225.030. The contributions of steric and electrostatic fields to the activity are 55.2% and 44.8%, respectively. For the CoMSIA study, the q², r², and F values of the model are 0.614, 0.929, and 88.771, respectively. The contributions of steric, electrostatic, hydrophobic, hydrogen bond donor, and hydrogen bond donor fields to the activity are 27.3%, 23.9%, 16.4%, 21.7%, and 10.7%, respectively. The CoMFA and CoMSIA models showed strong predictive ability, and the 3D contour plots give the basis on the structure modification of JAK2 inhibitors.     

3D-QSAR and Molecular Docking Studies on Derivatives of MK-0457, GSK1070916 and SNS-314 as Inhibitors against Aurora B Kinase

Development of anticancer drugs targeting Aurora B, an important member of the serine/threonine kinases family, has been extensively focused on in recent years. In this work, by applying an integrated computational method, including comparative molecular field analysis (CoMFA), comparative molecular similarity indices analysis (CoMSIA), homology modeling and molecular docking, we investigated the structural determinants of Aurora B inhibitors based on three different series of derivatives of 108 molecules. The resultant optimum 3D-QSAR models exhibited (q(2) = 0.605, r(2) (pred) = 0.826), (q(2) = 0.52, r(2) (pred) = 0.798) and (q(2) = 0.582, r(2) (pred) = 0.971) for MK-0457, GSK1070916 and SNS-314 classes, respectively, and the 3D contour maps generated from these models were analyzed individually. The contour map analysis for the MK-0457 model revealed the relative importance of steric and electrostatic effects for Aurora B inhibition, whereas, the electronegative groups with hydrogen bond donating capacity showed a great impact on the inhibitory activity for the derivatives of GSK1070916. Additionally, the predictive model of the SNS-314 class revealed the great importance of hydrophobic favorable contour, since hydrophobic favorable substituents added to this region bind to a deep and narrow hydrophobic pocket composed of residues that are hydrophobic in nature and thus enhanced the inhibitory activity. Moreover, based on the docking study, a further comparison of the binding modes was accomplished to identify a set of critical residues that play a key role in stabilizing the drug-target interactions. Overall, the high level of consistency between the 3D contour maps and the topographical features of binding sites led to our identification of several key structural requirements for more potency inhibitors. Taken together, the results will serve as a basis for future drug development of inhibitors against Aurora B kinase for various tumors.     

Investigation of Antigen-Antibody Interactions of Sulfonamides with a Monoclonal Antibody in a Fluorescence Polarization Immunoassay Using 3D-QSAR Models

A three-dimensional quantitative structure-activity relationship (3D-QSAR) model of sulfonamide analogs binding a monoclonal antibody (MAb(SMR)) produced against sulfamerazine was carried out by Distance Comparison (DISCOtech), comparative molecular field analysis (CoMFA), and comparative molecular similarity indices analysis (CoMSIA). The affinities of the MAb(SMR), expressed as Log(10)IC(50), for 17 sulfonamide analogs were determined by competitive fluorescence polarization immunoassay (FPIA). The results demonstrated that the proposed pharmacophore model containing two hydrogen-bond acceptors, two hydrogen-bond donors and two hydrophobic centers characterized the structural features of the sulfonamides necessary for MAb(SMR) binding. Removal of two outliers from the initial set of 17 sulfonamide analogs improved the predictability of the models. The 3D-QSAR models of 15 sulfonamides based on CoMFA and CoMSIA resulted in q(2) (cv) values of 0.600 and 0.523, and r(2) values of 0.995 and 0.994, respectively, which indicates that both methods have significant predictive capability. Connolly surface analysis, which mainly focused on steric force fields, was performed to complement the results from CoMFA and CoMSIA. This novel study combining FPIA with pharmacophore modeling demonstrates that multidisciplinary research is useful for investigating antigen-antibody interactions and also may provide information required for the design of new haptens.     

Structural Determinants of CX-4945 Derivatives as Protein Kinase CK2 Inhibitors: A Computational Study

Protein kinase CK2, also known as casein kinase-2, is involved in a broad range of physiological events including cell growth, proliferation and suppression of apoptosis which are related to human cancers. A series of compounds were identified as CK2 inhibitors and their inhibitory activities varied depending on their structures. In order to explore the structure-activity correlation of CX-4945 derivatives as inhibitors of CK2, in the present study, a set of ligand- and receptor-based 3D-QSAR models were developed employing Comparative Molecular Field Analysis (CoMFA) and Comparative Molecular Similarity Index Analysis (CoMSIA). The optimum CoMFA (R(cv) (2) = 0.618, R(pred) (2) = 0.892) and CoMSIA (R(cv) (2) = 0.681, R(pred) (2) = 0.843) models exhibited reasonable statistical characteristics for CX-4945 derivatives. The results indicated that electrostatic effects contributed the most to both CoMFA and CoMSIA models. The combination of docking analysis and molecular dynamics (MD) simulation showed that Leu45, Lys68, Glu81, Val116, Asp175 and Trp176 of CK2 which formed several direct or water-bridged H-bonds with CX-4945 are crucial for CX-4945 derivatives recognition to CK2. These results can offer useful theoretical references for designing more potent CK2 inhibitors.     

3D QSAR Studies,Pharmacophore Modeling and Virtual Screening on a Series of Steroidal Aromatase Inhibitors

Aromatase inhibitors are the most important targets in treatment of estrogen-dependent cancers. In order to search for potent steroidal aromatase inhibitors (SAIs) with lower side effects and overcome cellular resistance, comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) were performed on a series of SAIs to build 3D QSAR models. The reliable and predictive CoMFA and CoMSIA models were obtained with statistical results (CoMFA: q² = 0.636, r²_ncv = 0.988, r²_pred = 0.658; CoMSIA: q² = 0.843, r²_ncv = 0.989, r²_pred = 0.601). This 3D QSAR approach provides significant insights that can be used to develop novel and potent SAIs. In addition, Genetic algorithm with linear assignment of hypermolecular alignment of database (GALAHAD) was used to derive 3D pharmacophore models. The selected pharmacophore model contains two acceptor atoms and four hydrophobic centers, which was used as a 3D query for virtual screening against NCI2000 database. Six hit compounds were obtained and their biological activities were further predicted by the CoMFA and CoMSIA models, which are expected to design potent and novel SAIs.     

Probing Structural Features and Binding Mode of 3-Arylpyrimidin-2,4-diones within Housefly γ-Aminobutyric Acid (GABA) Receptor

In order to obtain structural features of 3-arylpyrimidin-2,4-diones emerged as promising inhibitors of insect γ-aminobutyric acid (GABA) receptor, a set of ligand-/receptor-based 3D-QSAR models for 60 derivatives are generated using Comparative Molecular Field Analysis (CoMFA) and Comparative Molecular Similarity Index Analysis (CoMSIA). The statistically optimal CoMSIA model is produced with highest q(2) of 0.62, r(2) (ncv) of 0.97, and r(2) (pred) of 0.95. A minor/bulky electronegative hydrophilic polar substituent at the 1-/6-postion of the uracil ring, and bulky substituents at the 3'-, 4'- and 5'-positions of the benzene ring are beneficial for the enhanced potency of the inhibitors as revealed by the obtained 3D-contour maps. Furthermore, homology modeling, molecular dynamics (MD) simulation and molecular docking are also carried out to gain a better understanding of the probable binding modes of these inhibitors, and the results show that residues Ala-183(C), Thr-187(B), Thr-187(D) and Thr-187(E) in the second transmembrane domains of GABA receptor are responsible for the H-bonding interactions with the inhibitor. The good correlation between docking observations and 3D-QSAR analyses further proves the model reasonability in probing the structural features and the binding mode of 3-arylpyrimidin-2,4-dione derivatives within the housefly GABA receptor.     

聚氧乙烯类和硫酸盐类表面活性剂HLB值的研究

亲水亲油平衡值(Hydrophile/Lipophile Balance,HLB)是表征表面活性剂性质的重要参数。采用比较分子力场分析(comparative molecular field analysis,CoMFA)方法和比较分子相似性指数分析(comparative similarity indices analysis,CoMSIA)方法,研究了30个烷基酚聚氧乙烯醚系列非离子表面活性剂的三维定量结构与pHLB值的关系。建立的CoMFA模型交叉验证系数为0.705,非交叉验证系数为0.805:CoMSIA模型交叉验证系数为0.685,非交叉验证系数为0.813,都有较好的预侧能力。此外,运用多元线性回归方法建立一种预报硫酸盐类表面活性剂HLB的QSPR模型,模型的复相关系数为0.711。所有选取的测试集和训练集的预报结果较为吻合,说明硫酸盐类表面活性剂HLB值QSPR模型具有良好的预报能力和普适性。     

Structural requirements of N-substituted spiropiperidine analogues as agonists of nociceptin/orphanin FQ receptor

The nociceptin/orphanin FQ (NOP) receptor is involved in a wide range of biological functions, including pain, anxiety, depression and drug abuse. Especially, its agonists have great potential to be developed into anxiolytics. In this work, both the ligand- and receptor-based three-dimensional quantitative structure-activity relationship (3D-QSAR) studies were carried out using comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) techniques on 103 N-substituted spiropiperidine analogues as NOP agonists. The resultant optimal ligand-based CoMSIA model exhibited Q(2) of 0.501, R(2) (ncv) of 0.912 and its predictive ability was validated by using an independent test set of 26 compounds which gave R(2) (pred) value of 0.818. In addition, docking analysis and molecular dynamics simulation (MD) were also applied to elucidate the probable binding modes of these agonists. Interpretation of the 3D contour maps, in the context of the topology of the active site of NOP, provided insight into the NOP-agonist interactions. The information obtained from this work can be used to accurately predict the binding affinity of related agonists and also facilitate the future rational design of novel agonists with improved activity.