Beta-amino acid substitutions and structure-based CoMFA modeling of hepatitis C virus NS3 protease inhibitors

In an effort to develop a new type of HCV NS3 peptidomimetic inhibitor, a series of tripeptide inhibitors incorporating a mix of alpha- and beta-amino acids has been synthesized. To understand the structural implications of beta-amino acid substitution, the P(1), P(2), and P(3) positions of a potent tripeptide scaffold were scanned and combined with carboxylic acid and acyl sulfonamide C-terminal groups. Inhibition was evaluated and revealed that the structural changes resulted in a loss in potency compared with the alpha-peptide analogues. However, several compounds exhibited muM potency. Inhibition data were compared with modeled ligand-protein binding poses to understand how changes in ligand structure affected inhibition potency. The P(3) position seemed to be the least sensitive position for beta-amino acid substitution. Moreover, the importance of a proper oxyanion hole interaction for good potency was suggested by both inhibition data and molecular modeling. To gain further insight into the structural requirements for potent inhibitors, a three-dimensional quantitative structure-activity relationship (3D-QSAR) model has been constructed using comparative molecular field analysis (CoMFA). The most predictive CoMFA model has q(2)=0.48 and r(pred)(2)=0.68.     

Structure-activity relationships and CoMFA of N-3 substituted phenoxypropyl piperidine benzimidazol-2-one analogues as NOP receptor agonists with analgesic properties

The N-3 position of a series of 3-phenoxypropyl piperidine benzimidazol-2-one analogues was optimised using the predictive power of a CoMFA model. The model was used to prioritise compounds for synthesis culminating in the triazole (+)-24. (+)-24 was found to be a high affinity, potent NOP agonist and demonstrated both antinociceptive and antiallodynic effects when administered iv to rodents.     

Synthesis and structure-activity relationships of N-substituted spiropiperidines as nociceptin receptor ligands

A series of N-substituted analogs based upon the spiropiperidine core of 1 was synthesized and exhibited high binding affinity to the nociceptin (NOP) receptor. The selectivities against other known opioid receptors were determined.     

3D-QSAR analysis of 2,4,5- and 2,3,4,5-substituted imidazoles as potent and nontoxic modulators of P-glycoprotein mediated MDR.

3D-Quantitative structure-activity relationships of 2,4,5- and 2,3,4,5-substituted imidazoles as a novel class of potent and nontoxic modulators of Pgp mediated MDR were investigated using CoMFA and COMSIA approaches. The best CoMFA model obtained from 46 imidazole analogues is a two-component model with the following statistics. R(2)(cv)=0.643, RMSE(cv)=0.360 for the cross-validation, and R(2)=0.767, RMSE=0.290 for the fitted. The best COMSIA model obtained is also a two-component model with the following statistics. R(2)(cv)=0.619, RMSE(cv)=0.372 for the cross-validation, and R(2)=0.765, RMSE=0.292 for the fitted.     

Molecular docking and 3-D-QSAR studies on the possible antimalarial mechanism of artemisinin analogues

Artemisinin (Qinghaosu) is a natural constituent found in Artemisia annua L, which is an effective drug against chloroquine-resistant Plasmodium falciparum strains and cerebral malaria. The antimalarial activities of artemisinin and its analogues appear to be mediated by the interactions of the drugs with hemin. In order to understand the antimalarial mechanism and the relationship between the physicochemical properties and the antimalarial activities of artemisinin analogues, we performed molecular docking simulations to probe the interactions of these analogues with hemin, and then performed three-dimensional quantitative structure-activity relationship (3-D-QSAR) studies on the basis of the docking models employing comparative molecular force fields analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA). Molecular docking simulations generated probable 'bioactive' conformations of artemisinin analogues and provided a new insight into the antimalarial mechanism. The subsequent partial least squares (PLS) analysis indicates that the calculate binding energies correlate well with the experimental activity values. The CoMFA and CoMSIA models based on the bioactive conformations proved to have good predictive ability and in turn match well with the docking result, which further testified the reliability of the docking model. Combining these results, that is molecular docking and 3-D-QSAR, together, the binding model and activity of new synthesized artemisinin derivatives were well explained.     

Essential structural profile of a dual functional inhibitor against cyclooxygenase-2 (COX-2) and 5-lipoxygenase (5-LOX): molecular docking and 3D-QSAR analyses on DHDMBF analogues

It is recently proposed that compounds with equal capabilities of inhibiting COX and 5-LOX, both are key enzymes involved in the arachidonic acid (AA) cascade, are expected to be safer non-steroidal anti-inflammatory drugs (NSAIDs). To dig out helpful information in designing dual functional inhibitors against the two enzymes, homology modeling, molecular dynamics (MD) simulations, automated docking, and 3D-QSAR analyses were performed in this study on 21 COX-2/5-LOX dual inhibitors, namely, 7-tert-butyl-2,3-dihydro-3,3-dimethylbenzofuran (DHDMBF) analogues. A 3D-model of 5-LOX was built based on the high-resolution X-ray structure of rabbit reticulocyte 15-lipoxygenase. Molecular docking was then applied to locate the binding orientations and conformations of DHDMBF analogues with COX-2 and 5-LOX, respectively, leading to highly predictive CoMFA models constructed on the basis of the binding conformations with q2 values of 0.782 and 0.634 for COX-2 and 5-LOX, respectively. In addition, CoMFA field distributions were found in good agreement with the structural characteristics of the corresponding binding sites. Both the docking simulations and QSAR analyses suggest that new potent dual inhibitors should share a structural feature with a moderately bulky group at R2 position and a rather negatively charged group around the position of the carbonyl group of DHDMBFs. Therefore, the final 3D-QSAR models and the information of the inhibitor-enzyme interaction should be useful in developing new NSAIDs as anti-inflammation drugs with favorable safety profile.     

3-D-QSAR analysis of N-(3-acyloxy-2-benzylpropyl)-N'-dihydroxytetrahydrobenzazepine and tetrahydroisoquinoline and N-(3-acyloxy-2-benzylpropyl)-N'-(4-hydroxy-3-methoxybenzyl) thioureas analogues as potent vanilloid receptor ligands

3-D-Quantitative structure--activity relationships of N-(3-acyloxy-2-benzylpropyl)-N'-dihydroxytetrahydro-benzazepine and tetrahydroisoquinoline and N-(3-acyloxy-2-benzylpropyl)-N'-(4-hydroxy-3-methoxybenzyl) thiourea analogues as potent vanilloid receptor ligands were investigated using the CoMFA and the COMSIA methods. The best CoMFA model obtained in this study from 29 substituted thiourea analogues is a two-component model with the following statistics. R(2)((cv))=0.407 and RMSE((cv))=0.532 for the cross-validation, and R(2)=0.705 and RMSE=0.375 for the fitted. The best COMSIA model obtained from the same 29 compounds is a two-component model with the following statistics: R(2)((cv))=0.336 and RMSE((cv))=0.563 for the cross-validation, and R(2)=0.693 and RMSE=0.382 for the fitted.     

Synthesis,biological activity and receptor-based 3-D QSAR study of 3'-N-substituted-3'-N-debenzoylpaclitaxel analogues

3'-N-Substituted-3'-N-debenzoylpaclitaxel analogues were synthesized and investigated for their 3-D QSAR by using comparative molecular field analysis (CoMFA). The CoMFA model obtained from receptor(microtubule)-paclitaxel binding structure displays an excellent predictive power to forecast the biological activity of new 3'-N-substituted-3'-N-debenzoylpaclitaxel analogues as well as the ability to explain the activity of the known paclitaxel analogues. The cross-validated r(2)(cv) values of the selected models are 0.835 and 0.616 for A549 and SK-OV-3, respectively, and the non-cross-validated r(2)(ncv) values of them are 0.992 and 0.974.     

Synthesis and evaluation of N-3 substituted phenoxypropyl piperidine benzimidazol-2-one analogues as NOP receptor agonists with analgesic and sedative properties

A series of 3-phenoxypropyl piperidine benzimidazol-2-one analogues have been discovered as novel NOP receptor agonists. Structure-activity relationships have been explored via N-3 substitution of the benzimidazol-2-one with a range of functionality. The N-methyl acetamide derivative (+)-7f was found to be a high-affinity, potent NOP agonist with greater than 100-fold selectivity over the MOP receptor. Furthermore (+)-7f was shown to be both antinociceptive and sedative when administered iv to rodents.     

Structural basis for androgen receptor agonists and antagonists: interaction of SPEED 98-listed chemicals and related compounds with the androgen receptor based on an in vitro reporter gene assay and 3D-QSAR

The androgen receptor (AR) activity of listed chemicals, so called SPEED 98, by the Ministry of the Environment, Japan, and structurally related chemicals was characterized using MDA-kb2 human breast cancer cells stably expressing an androgen-responsive luciferase reporter gene, MMTV-luc. Since our results suggested that chemicals with diverse chemical structures were capable of disrupting the endocrine systems mediated by AR, a comparative molecular field analysis (CoMFA) model was developed to analyze the structural requirements necessary to disrupt AR function. A significant CoMFA model with r(2)=0.825 and q(2)=0.332 was developed for AR antagonist activity of 35 pure antagonists excluding procymidone. On the other hand, a good CoMFA model with r(2)=0.983 and q(2)=0.555 was obtained for antagonist activity of 13 chemicals with both agonist and antagonist activities. The steric and electrostatic properties were sufficient to describe the structural requirements for AR antagonist activity. In addition, the structural difference of AR agonists and antagonists was explained based on CoMFA results and the AR-LBD crystal structure. As several ERalpha agonists such as diethylstilbestrol (DES) acted as AR antagonists, the surface area of the AR ligand-binding domain (LBD) was compared with that of the ERalpha-LBD based on their reported crystal structures to analyze how those ligands interact with LBDs. The surface area of AR-LBD was shown to be smaller than that of ERalpha-LBD and therefore compounds with both estrogenic and antiandrogenic activities can fit well into the ERalpha-LBD but may protrude from the AR-LBD. It is likely that this subtle difference of the surface areas of the LBDs determines whether an ERalpha agonist acts as an AR antagonist or an agonist.     

High-affinity carbamate analogues of morphinan at opioid receptors

A series of carbamate analogues were synthesized from levorphanol (1a), cyclorphan (2a) or butorphan (3a) and evaluated in vitro for their binding affinity at mu, delta, and kappa opioid receptors. Functional activities of these compounds were measured in the [(35)S]GTPgammaS binding assay. Phenyl carbamate derivatives 2d and 3d showed the highest binding affinity for kappa receptor (K(i)=0.046 and 0.051 nM) and for mu receptor (K(i)=0.11 and 0.12 nM). Compound 1c showed the highest mu selectivity. The preliminary assay for agonist and antagonist properties of these ligands in stimulating [(35)S]GTPgammaS binding mediated by the kappa opioid receptor illustrated that all of these ligands were kappa agonists. At the mu receptor, compounds 1b, 1c, 2b, and 3b were agonists, while compounds 2c-e and 3c-e were mu agonists/antagonists.     

Application of comparative molecular field analysis to dopamine D2 partial agonists

Comparative molecular field analysis (CoMFA) has been applied to novel D2 partial agonists. Due to the predictability of the CoMFA model across different series of D2 partial agonists, we believe these compounds are binding to the D2 agonist receptor in the conformation and alignment described herein.     

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.     

3D QSAR models for alpha2a-adrenoceptor agonists

Three-dimensional structure-activity relationship studies of alpha2a-adrenoceptor agonists were carried out by Distance Comparison (DISCOthech) and Comparative Molecular Field Analysis (CoMFA) methods to define the pharmacophore and a quantitative model, respectively, of this class of compounds. The statistical validation of the CoMFA model indicates its high predictive performance for binding affinities of new alpha2a-adrenoceptor agonists.     

3D-QSAR design of novel antiepileptic sulfamides

A three-dimensional quantitative structure-activity relationship method, the comparative molecular field analysis (CoMFA), was applied to design new anticonvulsant symmetric sulfamides. The training set (27 structures) was comprised by traditional and new-generation anticonvulsant (AC) ligands that exhibit a potent activity in MES test. Physicochemical determinants of binding, such as steric and electrostatic properties, were mapped onto the molecular structures of the set, in order to interpret graphically the CoMFA results in terms of field contribution maps. The 3D-QSAR models demonstrate a good ability to predict the activity of the designed compounds (r(2)=0.967, q(2)=0.756).     

Design, synthesis and CoMFA studies of N1-amino acid substituted 2,4,5-triphenyl imidazoline derivatives as p53-MDM2 binding inhibitors

A series of novel N1-amino-acid substituted 2,4,5-triphenyl imidazoline derivatives was designed and synthesized based on our previous studies. All synthesized target compounds were screened for their p53-MDM2 binding inhibitory activities and anti-proliferative activities against five cancer cell lines. Among them, twelve compounds displayed improved binding inhibitory activities and most compounds showed higher cell growth inhibition activities with IC(50) values in the low micromolar range. Compound 6c exhibited marked p53-MDM2 binding inhibitory activity (IC(50)=0.59 μM) which was eightfold more potent than that of Nutlin-1 (IC(50)=4.78 μM). CoMFA analysis was performed based on obtained biological data and resulted in a statistically significant CoMFA model with high predict abilities (q(2)=0.645, r(2)=0.979).     

Synthesis and structure–activity relationships of N-substituted spiropiperidines as nociceptin receptor ligands: Part 2

A series of N-8 substituted analogs based upon the spiropiperidine core of the original lead compound 1 was synthesized. This lead has been elaborated to compounds to give compounds 2 and 3 (R = H) that exhibited high NOP binding affinity as well as selectivity against other known opioid receptors. These two series have been further functionalized at the amido nitrogen. The synthesis and structure–activity relationship (SAR) of these and related compounds are discussed.     

Comparative molecular field analysis of colchicine inhibition and tubulin polymerization for combretastatins binding to the colchicine binding site on beta-tubulin

A molecular modeling study using Comparative Molecular Field Analysis (CoMFA) was undertaken to develop a predictive model for combretastatin binding to the colchicine binding site of tubulin. Furthermore, we examined the potential contribution of lipophilicity (log P) and molecular dipole moment and were unable to correlate these properties to the observed biological data. In this study we first confirmed that tubulin polymerization inhibition (IC50) correlated (R2 = 0.92) with [3H]colchicine displacement. Although these data correlated quite well, we developed two independent models for each set of data to quantify structural features that may contribute to each biological property independently. To develop our predictive model we first examined a series of molecular alignments for the training set and ultimately found that overlaying the respective trimethoxyphenyl rings (A ring) of the analogues generated the best correlated model. The CoMFA yielded a cross-validated R2 = 0.41 (optimum number of components equal to 5) for the tubulin polymerization model and an R2 = 0.38 (optimum number of components equal to 5) for [3H]colchicine inhibition. Final non-cross-validation generated models for tubulin polymerization (R2 of 0.93) and colchicine inhibition (R2 of 0.91). These models were validated by predicting both biological properties for compounds not used in the training set. These models accurately predicted the IC50 for tubulin polymerization with an R2 of 0.88 (n = 6) and those of [3H]colchicine displacement with an R2 of 0.80 (n = 7). This study represents the first predictive model for the colchicine binding site over a wide range of combretastatin analogues.