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.     

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.     

3D-QSAR Studies on Barbituric Acid Derivatives as Urease Inhibitors and the Effect of Charges on the Quality of a Model

Urease enzyme (EC 3.5.1.5) has been determined as a virulence factor in pathogenic microorganisms that are accountable for the development of different diseases in humans and animals. In continuance of our earlier study on the helicobacter pylori urease inhibition by barbituric acid derivatives, 3D-QSAR (three dimensional quantitative structural activity relationship) advance studies were performed by Comparative Molecular Field Analysis (CoMFA) and Comparative Molecular Similarity Indices Analysis (CoMSIA) methods. Different partial charges were calculated to examine their consequences on the predictive ability of the developed models. The finest developed model for CoMFA and CoMSIA were achieved by using MMFF94 charges. The developed CoMFA model gives significant results with cross-validation (q²) value of 0.597 and correlation coefficients (r²) of 0.897. Moreover, five different fields i.e., steric, electrostatic, and hydrophobic, H-bond acceptor and H-bond donors were used to produce a CoMSIA model, with q² and r² of 0.602 and 0.98, respectively. The generated models were further validated by using an external test set. Both models display good predictive power with r²_pred ≥ 0.8. The analysis of obtained CoMFA and CoMSIA contour maps provided detailed insight for the promising modification of the barbituric acid derivatives with an enhanced biological activity.     

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.     

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.     

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.     

Docking-Based 3D-QSAR Studies for 1,3,4-oxadiazol-2-one Derivatives as FAAH Inhibitors

This work aimed to construct 3D-QSAR CoMFA and CoMSIA models for a series of 31 FAAH inhibitors, containing the 1,3,4-oxadiazol-2-one moiety. The obtained models were characterized by good statistical parameters: CoMFA Q² = 0.61, R² = 0.98; CoMSIA Q² = 0.64, R² = 0.93. The CoMFA model field contributions were 54.1% and 45.9% for steric and electrostatic fields, respectively. In the CoMSIA model, electrostatic, steric, hydrogen bond donor, and hydrogen acceptor properties were equal to 34.6%, 23.9%, 23.4%, and 18.0%, respectively. These models were validated by applying the leave-one-out technique, the seven-element test set (CoMFA r²_test-set = 0.91; CoMSIA r²_test-set = 0.91), a progressive scrambling test, and external validation criteria developed by Golbraikh and Tropsha (CoMFA r²₀ = 0.98, k = 0.95; CoMSIA r²₀ = 0.98, k = 0.89). As the statistical significance of the obtained model was confirmed, the results of the CoMFA and CoMSIA field calculation were mapped onto the enzyme binding site. It gave us the opportunity to discuss the structure–activity relationship based on the ligand–enzyme interactions. In particular, examination of the electrostatic properties of the established CoMFA model revealed fields that correspond to the regions where electropositive substituents are not desired, e.g., in the neighborhood of the 1,3,4-oxadiazol-2-one moiety. This highlights the importance of heterocycle, a highly electronegative moiety in this area of each ligand. Examination of hydrogen bond donor and acceptor properties contour maps revealed several spots where the implementation of another hydrogen-bond-donating moiety will positively impact molecules’ binding affinity, e.g., in the neighborhood of the 1,3,4-oxadiazol-2-one ring. On the other hand, there is a large isopleth that refers to the favorable H-bond properties close to the terminal phenoxy group of a ligand, which means that, generally speaking, H-bond acceptors are desired in this area.     

The Three Dimensional Quantitative Structure Activity Relationships (3D-QSAR) and Docking Studies of Curcumin Derivatives as Androgen Receptor Antagonists

Androgen receptor antagonists have been proved to be effective anti-prostate cancer agents. 3D-QSAR and Molecular docking methods were performed on curcumin derivatives as androgen receptor antagonists. The bioactive conformation was explored by docking the potent compound 29 into the binding site of AR. The constructed Comparative Molecular Field Analysis (CoMFA) and Comparative Similarity Indices Analysis (CoMSIA) models produced statistically significant results with the cross-validated correlation coefficients q(2) of 0.658 and 0.567, non-cross-validated correlation coefficients r(2) of 0.988 and 0.978, and predicted correction coefficients r(2) (pred) of 0.715 and 0.793, respectively. These results ensure the CoMFA and CoMSIA models as a tool to guide the design of novel potent AR antagonists. A set of 30 new analogs were proposed by utilizing the results revealed in the present study, and were predicted with potential activities in the developed models.     

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.     

Molecular Modeling Studies of 4,5-Dihydro-1H-pyrazolo[4,3-h] quinazoline Derivatives as Potent CDK2/Cyclin A Inhibitors Using 3D-QSAR and Docking

CDK2/cyclin A has appeared as an attractive drug targets over the years with diverse therapeutic potentials. A computational strategy based on comparative molecular fields analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) followed by molecular docking studies were performed on a series of 4,5-dihydro-1H-pyrazolo[4,3-h]quinazoline derivatives as potent CDK2/cyclin A inhibitors. The CoMFA and CoMSIA models, using 38 molecules in the training set, gave r(2) (cv) values of 0.747 and 0.518 and r(2) values of 0.970 and 0.934, respectively. 3D contour maps generated by the CoMFA and CoMSIA models were used to identify the key structural requirements responsible for the biological activity. Molecular docking was applied to explore the binding mode between the ligands and the receptor. The information obtained from molecular modeling studies may be helpful to design novel inhibitors of CDK2/cyclin A with desired activity.     

Investigation of the Structure Requirement for 5-HT(6) Binding Affinity of Arylsulfonyl Derivatives: A Computational Study

5-HT(6) receptor has been implicated in a series of diseases including anxiety, depression, schizophrenia and cognitive dysfunctions. 5-HT(6) ligands have been reported to play a significant role in the treatment for central nervous system (CNS) diseases. Presently, a large series of 223 5-HT(6) ligands were studied using a combinational method by 3D-QSAR, molecular docking and molecular dynamics calculations for further improvement of potency. The optimal 3D models exhibit satisfying statistical results with r(2) (ncv), q(2) values of 0.85 and 0.50 for CoMFA, 0.81 and 0.53 for CoMSIA, respectively. Their predictive powers were validated by external test set, showing r(2) (pred) of 0.71 and 0.76. The contour maps also provide a visual representation of contributions of steric, electrostatic, hydrophobic and hydrogen bond fields as well as the prospective binding models. In addition, the agreement between 3D-QSAR, molecular docking and molecular dynamics simulation proves the rationality of the developed models. These results, we hope, may be helpful in designing novel and potential 5-HT(6) ligands.     

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.     

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.     

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.     

靛玉红类衍生物的三维定量构效关系研究

采用比较分子场方法(CoMFA)对16个靛玉红类衍生物进行三维定量构效关系(3D-QSAR)研究,探讨靛玉红类衍生物的抗癌活性与其结构的关系。建立了3D-QSAR的CoMFA模型,得到了非交叉验证相关系数r2=0.987,标准偏差SD=2.141,统计方差比F=117.136,交叉验证相关系数q2=0.496。该模型合理、可信,具有一定的预测能力,表明改变取代基R1基团的体积才是提高靛玉红类衍生物抗癌活性的主要途径,为开发活性更好的抗癌药物提供一定的理论参考。     

3D-QSAR Studies of Polyazaheterocyclic Ligands Used in Lanthanide and Actinide Extraction Processes

Three-dimensional quantitative structure activity relationships (3D-QSAR) methods of comparative molecular field analysis (CoMFA), CoMFA region focusing (CoMFA-RF), and comparative molecular similarity indices analysis (CoMSIA) were successfully employed to predict and model of the separation factors of N-donor heterocyclic extractants in the separation of americium(III) from europium(III) ions from acidic aqueous solutions. An all-orientation search (AOS) strategy was used to acquire the best orientation and minimize the effect of the initial orientation of aligned compounds. The obtained 3D-QSAR models showed R² > 0.97 and an absolute error of > 0.1 in log unit of the separation factors. Validation, reliability, and robustness of 3D-QSAR models were evaluated by the prediction of external test sets, leave-one-out, leave-four-out, bootstrapping, and progressive scrambling approaches. The results confirmed that steric hindrance and electrostatic interactions besides hydrophobic effects of N-donor extractants play an important role in Am(III)/Eu(III) separation processes. The obtained information could be very useful in designing the most efficient ligands and finding new extractants for lanthanide and actinide separation.     

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.     

Modelling Atypical Small‐Molecule Mimics of an Important Stem Cell Cytokine,Thrombopoietin

We report the first comprehensive 3D QSAR study of a large, structurally diverse set of compounds that act as atypical thrombopoietin (TPO) mimics by interacting with the transmembrane domain of the TPO receptor, c‐MPL. These agonists of c‐MPL were superimposed according to a pharmacophore hypothesis, resulting in 3D QSAR models of high statistical significance. The pharmacophore‐based superimposition and comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) were used to derive the QSAR models relating structure to the published in vitro bioactivities of the TPO mimics. The CoMFA and CoMSIA models gave high correlation coefficients of the bioactivities with the derived fields, resulting in robust prediction of agonist activity of the superimposed compounds. The models have been interpreted in terms of the requirements for binding to the transmembrane domain of the TPO receptor.