A comparative molecular field analysis and molecular modelling studies on pyridylimidazole type of angiotensin II antagonists

A large number of compounds known as “AII (Angiotensin II) antagonists” have been developed for the treatment of various heart diseases such as hypertension, congestive heart failure, and chronic renal failure. Most of the currently known AII receptor antagonists share a similar chemical structure, consisting of nitrogen atoms, a lipophilic alkyl side chain and an acidic group. As a new series, we have designed and synthesized various pyridylimidazole derivatives. In this report we would like to discuss the structure–activity relationship of these series of compounds using the comparative molecular field analysis (CoMFA) methods. We could come up with a good CoMFA model (cross-validated and conventional r² values equal to 0.702 and 0.991, respectively) and the validity of the model was confirmed by synthesizing and measuring their biological activity of additional 6 compounds suggested by the model. This result provides additional information on the structural requirement for structurally diverse group of AII receptor antagonists.     

Combined 3D QSAR and molecular docking studies to reveal novel cannabinoid ligands with optimum binding activity

The combination of NMR spectroscopy and molecular modeling studies provided the putative bioactive conformation for the analgesic cannabinoid (CB) ligand (−)-2-(6a,7,10,10a-tetrahydro-6,6,9-trimethylhydroxy-6H-dibenzo[b,d]pyranyl)-2-hexyl 1,3-dithiolane which served as a template in reported three-dimensional quantitative structure–activity relationship (3D QSAR) studies [Durdagi et al., J. Med. Chem. 2007, 50, 2875]. The reported 3D models of the CB1 receptor allowed us to construct a new 3D QSAR model based on theoretical calculations and molecular docking studies. Statistical comparison of the constructed two 3D QSAR studies showed the improvement of the new model. In addition, the new model can explain more effectively the experimental data and thus it can serve more efficiently in the rational drug design of pharmacologically optimized CB analogues.     

Three-dimensional structure-activity relationship modeling of cocaine binding to two monoclonal antibodies by comparative molecular field analysis

Successful immunotherapy of cocaine addiction and overdoses requires cocaine-binding antibodies with specific properties, such as high affinity and selectivity for cocaine. We have determined the affinities of two cocaine-binding murine monoclonal antibodies (mAb: clones 3P1A6 and MM0240PA) for cocaine and its metabolites by [3H]-radioligand binding assays. mAb 3P1A6 (K(d) = 0.22 nM) displayed a 50-fold higher affinity for cocaine than mAb MM0240PA (K(d) = 11 nM) and also had a greater specificity for cocaine. For the systematic exploration of both antibodies' binding specificities, we used a set of approximately 35 cocaine analogues as structural probes by determining their relative binding affinities (RBAs) using an enzyme-linked immunosorbent competition assay. Three-dimensional quantitative structure-activity relationship (3D-QSAR) models on the basis of comparative molecular field analysis (CoMFA) techniques correlated the binding data with structural features of the ligands. The analysis indicated that despite the mAbs' differing specificities for cocaine, the relative contributions of the steric (approximately 80%) and electrostatic (approximately 20%) field interactions to ligand-binding were similar. Generated three-dimensional CoMFA contour plots then located the specific regions about cocaine where the ligand/receptor interactions occurred. While the overall binding patterns of the two mAbs had many features in common, distinct differences were observed about the phenyl ring and the methylester group of cocaine. Furthermore, using previously published data, a 3D-QSAR model was developed for cocaine binding to the dopamine reuptake transporter (DAT) that was compared to the mAb models. Although the relative steric and electrostatic field contributions were similar to those of the mAbs, the DAT cocaine-binding site showed a preference for negatively charged ligands. Besides establishing molecular level insight into the interactions that govern cocaine binding specificity by biopolymers, the three-dimensional images obtained reflect the properties of the mAbs binding pockets and provide the initial information needed for the possible design of novel antibodies with properties optimized for immunotherapy.     

Application of quantitative structure‐activity relationship analysis on an antibody and alternariol‐like compounds interaction study

The antigen‐antibody interaction determines the sensitivity and specificity of competitive immunoassay for hapten detection. In this paper, the specificity of a monoclonal antibody against alternariol‐like compounds was evaluated through indirect competitive ELISA. The results showed that the antibody had cross‐reactivity with 33 compounds with the binding affinity (expressed by IC₅₀) ranging from 9.4 ng/mL to 12.0 μg/mL. All the 33 compounds contained a common moiety and similar substituents. To understand how this common moiety and substituents affected the recognition ability of the antibody, a three‐dimensional quantitative structure‐activity relationship (3D‐QSAR) between the antibody and the 33 alternariol‐like compounds was constructed using comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) methods. The q² values of the CoMFA and CoMSIA models were 0.785 and 0.782, respectively, and the r² values were 0.911 and 0.988, respectively, indicating that the models had good predictive ability. The results of 3D‐QSAR showed that the most important factor affecting antibody recognition was the hydrogen bond mainly formed by the hydroxyl group of alternariol, followed by the hydrophobic force mainly formed by the methyl group. This study provides a reference for the design of new hapten and the mechanisms for antibody recognition.     

A comparative molecular field analysis of the biotransformation of sulfides by Rhodococcus erythropolis

A comparative molecular field analysis (CoMFA) was used to model the efficacy with which the Rhodococcus erythropolis mono-oxygenase, DszC, catalyzes the enantioselective sulfoxidation of a broad range of substrates. Experimentally determined values of both the yield and enantiomeric excess for this reaction were employed to create these CoMFA models. A highly predictive CoMFA model was constructed for the prediction of enantiomeric excess of the sulfoxide product. The predictive ability of the model was demonstrated by both cross-validation of the training set (q² = 0.74) and for an external test set of substrates. The enantiomeric excesses of the members of the test set, which also included two amino acid sulfides that were structurally distinct from the membership of the training set, were predicted well by the CoMFA model. Product yield was not modelled well by any CoMFA model. Different models comparing the likely bioactive conformations of the substrates suggest that most compounds assume an ‘extended’ conformation upon binding. Contour diagrams illustrating significant substrate–enzyme interactions suggest that the model, which predicts the enantiomeric excess, is consistent with previous conclusions regarding the effect of various substrate substitutions on the enantiopurity of the product of the biotransformation.     

Physicochemical explanation of peptide binding to HLA-A*0201 major histocompatibility complex: a three-dimensional quantitative structure-activity relationship study

A three-dimensional quantitative structure-activity relationship method for the prediction of peptide binding affinities to the MHC class I molecule HLA-A*0201 was developed by applying the CoMSIA technique on a set of 266 peptides. To increase the self consistency of the initial CoMSIA model, the poorly predicted peptides were excluded from the training set in a stepwise manner and then included in the study as a test set. The final model, based on 236 peptides and considering the steric, electrostatic, hydrophobic, hydrogen bond donor, and hydrogen bond acceptor fields, had q2 = 0.683 and r2 = 0.891. The stability of this model was proven by cross-validations in two and five groups and by a bootstrap analysis of the non-cross-validated model. The residuals between the experimental pIC50 (-logIC50) values and those calculated by "leave-one-out" cross-validation were analyzed. According to the best model, 63.2% of the peptides were predicted with /residuals/ < or = 0.5 log unit; 29.3% with 1.0 < or = /residuals/ < 0.5; and 7.5% with /residuals/ > 1.0 log unit. The mean /residual/ value was 0.489. The coefficient contour maps identify the physicochemical property requirements at each position in the peptide molecule and suggest amino acid sequences for high-affinity binding to the HLA-A*0201 molecule.     

Synthesis of novel hetero ring fused pyridine derivatives; Their anticancer activity,CoMFA and CoMSIA studies

A series of novel furo[2,3-b]pyridine-2-carboxamide 4a–h/pyrido[3′,2′:4,5]furo[3,2-d] pyrimidin-4(3H)-one derivatives 5a–p were prepared from pyridin 2(1H) one 1 via selective O-alkylation with α-bromoethylester followed by cyclization, then reaction with different aliphatic primary amines to obtain 4 and further reaction with triethyl orthoacetate/triethyl orthoformate. Also prepared novel furo[2,3-b]pyridine-2-carbohydrazide Schiff’s bases 7a–h and pyrido [3′,2′:4,5]furo[3,2-d]pyrimidin-4(3H)-one derivatives 8a–h starting from furo[2,3-b]pyridine carboxylate derivatives 3 by reaction with hydrazine hydrate to form 6 and reaction with diverse substituted aldehydes and cyclization. Products 4a–h, 5a–p, 7a–h and 8a–h were screened against four human cancer cell lines (HeLa, COLO205, Hep G2 and MCF 7) and one normal cell line (HEK 293). Compounds 4e, 4f, 4g, 5h, 7c, 7d, 7e and 7f showed significant anticancer activity against all the cell lines at micro molar concentration and found to be non-toxic to normal cell line. Studies for HeLa, COLO205 and MCF-7 using CoMFA and CoMSIA. Models from 3D-QSAR provided a strong basis for future rational design of more active and selective HeLa, COLO205 and MCF-7 cell line inhibitors.     

3D QSAR studies on peroxisome proliferator-activated receptor γ agonists using CoMFA and CoMSIA

The peroxisome proliferator-activated receptors (PPARs) have increasingly become attractive targets for developing novel anti-type 2 diabetic drugs. We employed comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) to study three-dimensional quantitative structure–activity relationship (3D QSAR) based on existing agonists of PPAR (including five thiazolidinediones and 74 tyrosine-based compounds). Predictive 3D QSAR models with conventional r² and cross-validated coefficient (q²) values up to 0.974 and 0.642 for CoMFA and 0.979 and 0.686 for COMSIA were established using the SYBYL package. These models were validated by a test set containing 18 compounds. The CoMFA and CoMSIA field distributions are in general agreement with the structural characteristics of the binding pockets of PPAR, which demonstrates that the 3D QSAR models built here are very useful in predicting activities of novel compounds for activating PPAR.      

Synthesis and comparative molecular field analysis (CoMFA) of antitumor 3-arylisoquinoline derivatives

In this study a series of 3-arylisoquinoline derivatives were synthesized and cytotoxicity against human melanoma tumor cell evaluated, and a three dimensional quantitative structure—activity relationship was investigated using the comparative molecular field analysis (CoMFA). The results suggested that the electrostatic, steric and hydrophobic factors of 3-arylisoquinolines were strongly correlated with the antitumor activity. Considerable predictive ability (cross-validated r² as high as 0.721) was obtained through CoMFA.     

Exhaustive CoMFA and CoMSIA analyses around different chemical entities: a ligand-based study exploring the affinity and selectivity profiles of 5-HT1A ligands

The 5-hydroxytryptamine (5-HT_1A) receptors represent an attractive target in drug discovery. In particular, 5-HT_1A agonists and partial agonists are deeply investigated for their potential role in the treatment of anxiety, depression, ischaemic brain disorder and more recently, of pain. On the other hand, 5-HT_1A antagonists have been revealed promising compounds in cognition disorders and, lately, in cancer. Thus, the discovery of 5HT_1A ligands is nowadays an appealing research activity in medicinal chemistry. In this work, Comparative Molecular Fields Analysis (CoMFA) and Comparative Molecular Similarity Index Analysis (CoMSIA) were applied on an in-house library of 5-HT_1A ligands bearing different chemical scaffolds in order to elucidate their affinity and selectivity for the target_. Following this procedure, a number of structural modifications have been drawn for the development of much more effective 5-HT_1AR ligands.

     

3D-QSAR studies of Dipeptidyl peptidase IV inhibitors using a docking based alignment

Dipeptidyl peptidase IV (DPP-IV) deactivates the incretin hormones GLP-1 and GIP by cleaving the penultimate proline or alanine from the N-terminal (P1-position) of the peptide. Inhibition of this enzyme will prevent the degradation of the incretin hormones and maintain glucose homeostasis; this makes it an attractive target for the development of drugs for diabetes. This paper reports 3D-QSAR analysis of several DPP-IV inhibitors, which were aligned by the receptor-based technique. The conformation of the molecules in the active site was obtained through docking methods. The QSAR models were generated on two training sets composed of 74 and 25 molecules which included phenylalanine, thiazolidine, and fluorinated pyrrolidine analogs. The 3D-QSAR models are robust with statistically significant r², q², and values. The CoMFA and CoMSIA models were used to design some new inhibitors with several fold higher binding affinity. Figure The CoMFA contours around molecule D1T155 (a) steric contours - favored (green); disfavored (yellow) (b) electrostatic contours - electropositive (blue); electronegative (red)     

Structure-activity relationship analysis of 3-phenylpyrazole derivatives as androgen receptor antagonists

Abstract

Prostate cancer (PCa), an epithelial malignancy that occurs in the prostate, reminds the second leading cause of cancer-related incidence in men worldwidely. Androgen receptor antagonists are the main therapeutic strategy of PCa, which can block the binding of androgen to androgen receptors. However, the long-term treatment of marketed anti-androgens in patients can inevitably cause drug resistance problem. The research of searching for new drugs with novel skeleton is always on the way. Recently, a series of 3-phenylpyrazole derivatives were reported to antagonize the function of AR, but their efficiencies are not good enough and need to be improved. In this work, comparative molecular field analysis and comparative molecular similarity indices analysis methods were employed to study the structure activity relationships of these derivatives. Two different methods were used to obtain the optimal molecular conformation alignments, one is based on atomic alignment and the other is based on molecular docking. The final result shows that both these two strategies can obtain satisfactory results and the atomic alignment performs a little better than docking. The models illustrate the key structural features highly related with the androgenic bioactivity and provide valuable suggestions for the design of new androgen receptor antagonists in future.

Communicated by Ramaswamy H. Sarma     

Models of Monoamine Oxidase A and B Active Sites Obtained by using 3D QSAR with CoMFA Analysis

Abstract

The monoamine oxidase catalyses the oxidative deamination of neuroactive amines. This enzyme exists in two forms A and B, which differ by substrates preference and inhibitors specificity. Investigation of the structures of these enzymes and design new selective inhibitors are of greatly interesting since MAO A inhibitors are used in therapeutic practice as antidepressants and MAO B inhibitors – in the treatment Parkinson's diseases. The three dimension structures of monoamine oxidases are still unknown. Therefore, one of the most perspective approach to define significant features of structure active site is method based on analysis of structure-activity relationship (3D QSAR) with comparison of molecular fields analysis (CoMFA) allowing to get the spatial distribution of important properties affecting the activity.

In present study we investigate the structures of active sites MAO A and B using 16 pyrazinocarbazole derivatives in variant conformation. Majority of pyrazinocarbazole derivatives have a rigit conformation, but three of those is sufficiently flexible. The latters can be in two conformation types: long molecules (substitution accommodate along axis of main structure) and short molecules (substitution accommodate at acute angle about of main structure). Several 3D QSAR and CoMFA models of MAO A and B active sites were design for data sets containing various types of flexible molecules conformation. All obtained models are statistical reliable and have sufficient predictive power for tested compound tetrindole. The best MAO A model that include two flexible molecules in long conformations was obtained, and the longest one of those in short conformation. In contrast, for MAO B model containing all flexible molecules in the short conformations is more preferred.

On the basis of obtained data the schematic models of MAO A and B active sites structures are proposed. According to these models MAO A active site have the narrow long cavity that accommodate long molecules, while MAO B active site is broader and shorter.     

Structure-activity relationship study of human liver microsomes-catalyzed hydrolysis rate of ester prodrugs of MENT by comparative molecular field analysis (CoMFA)

A series of MENT esters (3-71) was designed, prepared and tested to study the structure-activity relationship (SAR) of the hydrolysis rate with human liver microsomes of these prodrugs. Compounds were obtained covering a wide range of metabolic stability. The results are useful for the proper selection of prodrugs for different pharmaceutical formulations to deliver the potent and prostate-sparing androgen MENT. The MENT esters can especially be administered for male hormone replacement therapy and male contraception. Comparative molecular field analysis (CoMFA) was applied to a dataset of 28 esters, for which ED50 values could be obtained. The CoMFA model where the electrostatic and H-bond molecular fields were combined turned out to be most predictive. Despite the limited size of the dataset, CoMFA can help to rationalize the SAR of the ester hydrolysis rate of ester prodrugs of MENT.     

Investigation of the diastereomerism of dihydrobenzoxathiin SERMs for ER alpha by molecular modeling

Molecular dynamics simulations were performed to investigate the distinct uterine activity of ten dihydrobenzoxathiin diastereomers against human estrogen receptor (ER) α. These diastereomers share similar binding mode to ER α ligand binding domain (LBD). Dihydrobenzoxathiin diastereomers with full antagonistic activity form more stable hydrogen bonds with Glu353 and His524 of ER α LBD than corresponding diastereomers. The molecular mechanics based generalized born surface area (MM-GBSA) analysis revealed that van der Waals interactions are predominant to the binding of dihydrobenzoxathiin diastereomers to ER α LBD. The per-residue free energy decomposition analysis revealed that the uterine activity difference is contributed mainly by electrostatic interactions. Our study provides mechanistic insights into the difference of uterine activity for dihydrobenzoxathiin diastereomers.     

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.     

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.     

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.     

Benzimidazole-based derivatives as privileged scaffold developed for the treatment of the RSV infection: a computational study exploring the potency and cytotoxicity profiles

Respiratory syncytial virus (RSV) has been identified as a main cause of hospitalisation in infants and children. To date, the current therapeutic arsenal is limited to ribavirin and palivizumab with variable efficacy. In this work, starting from a number of in-house series of previously described anti-RSV agents based on the benzimidazole scaffold, with the aim at gaining a better understanding of the related chemical features involved in potency and safety profiles, we applied a computational study including two focussed comparative molecular fields analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA). The results allowed us to derive useful suggestions for the design of derivatives and also to set up statistical models predicting the potency and selectivity index (SI?=?CC₅₀/EC₅₀) of any new analogue prior to synthesis. Accordingly, here, we discuss preliminary results obtained through the applied exhaustive QSAR analyses, leading to design and synthesise more effective anti-RSV agents.