3-D QSAR analysis of inhibition of murine soluble epoxide hydrolase (MsEH) by benzoylureas, arylureas, and their analogues

Two hundred and seventy-one compounds including benzoylureas, arylureas and related compounds were assayed using recombinant murine soluble epoxide hydrolase (MsEH) produced from a baculovirus expression system. Among all the insect growth regulators assayed, 18 benzoylphenylurea congeners showed weak activity against MsEH. Newly synthesized cyclohexylphenylurea, 1-benzyl-3-phenylurea, and 1,3-dibenzylurea analogues were rather potent. The introduction of a methyl group at the para-position of the phenyl ring of cyclohexylphenylurea enhanced the activity 6-fold, though similar substituent effects were not seen for any of the benzoylphenylureas. The activities of these compounds, including several previously reported compounds, such as dicyclohexylurea, diphenylurea, and their related analogues (Morisseau et al., Proc. Natl. Acad. Sci., 1999, 96, 8849), were quantitatively analyzed using comparative molecular field analysis (CoMFA), a three-dimensional quantitative structure-activity relationship (3-D QSAR) method. Both steric and electrostatic factors contributing to variations in the activity were visualized using CoMFA. CoMFA results showed that one side of the cyclohexylurea moiety having a trans-amide conformation (A-ring moiety) is surrounded by large sterically unfavorable fields, while the other side of A-ring moiety and the other cyclohexyl group (B-ring moiety) is encompassed by sterically favored fields. Electrostatically negative fields were scattered around the entire molecule, and a positive field surrounds the carbon of the carbonyl group. Hydrophobic fields were visualized using Kellogg's hydropathic interaction (HINT) in conjunction with CoMFA. Hydrophobically favorable fields appeared beside the 4- and 4'-carbon atoms of the cyclohexyl groups, and hydrophobically unfavorable fields surrounded the urea bridge. The addition of the molecular hydrophobicity, log P [corrected], to CoMFA did not improve the correlation significantly. The ligand-binding interactions shown by X-ray crystallographic data were rationalized using the results of the CoMFA and HINT analyses, and the essential physicochemical parameters for the design of new MsEH inhibitors were disclosed.     

Molecular docking and QSAR analyses of aromatic heterocycle thiosemicarbazone analogues for finding novel tyrosinase inhibitors

A collection of 36 thiosemicarbazone analogues possessed a broad span of tyrosinase inhibitory activities was designed and obtained. Robust and reliable CoMFA and CoMSIA models were gained to predict the structure–activity relationship and the new modifier direction. Inhibitory activities of the compounds were found to greatly depend upon molecular shape, size, and charge. The sterically bulky group at the C-4 position of the thiophene ring contributed a high capacity for biological activity. Some bulky substituents at the C1-position and C12-position, and electron-negative groups at the C3-position, helped to improve the activity of these analogues. The molecular docking results provided visual evidence for QSAR analysis and detailed information about binding mode, affinity, and the principal mechanism between the ligands and tyrosinase. Based on these, a prospective structure modification and optimization of the most potent compound, T32, was suggested for further research.     

High-throughput screening of ecdysone agonists using a reporter gene assay followed by 3-D QSAR analysis of the molting hormonal activity

In this study, 172 diacylhydrazine analogs were examined for their ability to activate an ecdysone (molting hormone)-dependent reporter gene in a silkworm (Bombyx mori) cell-based high-throughput screening assay. The measured EC(50) values (concentration required to cause an effect in 50% of the cells) were used to construct a 3-D QSAR model that describes the ecdysone agonist activities of the diacylhydrazine analogs. Of these compounds, 14 exhibited no activity and were excluded from the 3-D QSAR analysis. The resulting equation described approximately 74% of the activity for 158 compounds. The final equation consisted of 42% electrostatic and 58% steric effects (r(2) = 0.74 and q(2) = 0.45). Comparative molecular field analysis (CoMFA) was used to visualize the steric and electrostatic potential fields that were favorable and unfavorable for biological activity. Of particular interest was the observation that the hydrophobic parameter (logP) was not necessary for describing the observed activities, although previous studies have cited the importance of hydrophobic parameters in both classical and 3-D QSAR analyses of these compounds. Modeling studies of the B. mori ecdysone receptor supported the observed physicochemical parameters required for activity reported by the CoMFA models. Comparison of the present analysis with those performed using other lepidopteran assay systems evidenced a high degree of correlation (r(2) = 0.81 for a Sf-9 cell-based assay and r(2) = 0.89 for a Chilo suppressalis integument-based assay), indicating that it is valid to compare the results generated with the B. mori cell-based system to those generated with previous lepidopteran assays. This novel assay system is amendable to a high-throughput screening format and should greatly increase our ability to discover novel agonists of molting hormone (ecdysone) activity.     

3D-QSAR analysis on benzazole derivatives as eukaryotic topoisomerase II inhibitors by using comparative molecular field analysis method

Selective topoisomerase II inhibitors have created a great deal of interest in recent years for the design of new antitumoral compounds. 3D-QSAR analysis has been performed on a series of previously synthesized benzoxazole, benzimidazole, and oxazolo(4,5-b)pyridine derivatives, which are screened as eukaryotic topoisomerase II inhibitors, using comparative molecular field analysis (CoMFA) with partial least squares fit to predict the steric and electrostatic molecular field interactions for the activity. The CoMFA study was carried out using a training set of 16 compounds. The predictive ability of the model was assessed using a test set of 7 compounds. The analyzed 3D-QSAR CoMFA model has demonstrated a good fit, having r(2) value of 0.997 and cross-validated coefficient q(2) value as 0.435 for the model. The obtained model reveals that the electronegatively charged substituents such as NO(2) or COOCH(3) group on position R and/or R(1) at the heterocyclic ring system and positively charged atom and/or atom groups located between the benzazole moiety and 2-substituted phenyl ring as a bridge element improve the activity. On the other hand, a bulky substituent, such as methoxy group, attached to the ortho position of 2-phenyl-5-nitro-benzoxazole (1) enhances the activity similar to compound 13, which is both a meta and para substituent of the phenyl group attached to the 2-position of benzimidazole ring system, fit into the favored steric region to improve the activity.     

3D-QSAR study of adamantyl N-benzylbenzamides as melanogenesis inhibitors

Three-dimensional quantitative structure–activity relationship (3D-QSAR) modeling, comparative molecular field analysis (CoMFA), and comparative molecular similarity indices analysis (CoMSIA) of polyhydroxylated N-benzylbenzamide derivatives containing an adamantyl moiety were performed to understand the mechanism of action and structure–activity relationship of these compounds. Contour map analysis indicated that steric contributions of the adamantyl moiety and electrostatic contributions of the hydroxyl group at the 3-position are important in the activity. Activities of the training set and test sets predicted by CoMFA fit well with actual activities, demonstrating that CoMFA, along with the best calculated q² value, has the best predictive ability.     

酚酸类化感自毒物质对枸杞种子萌发的抑制作用

酚酸类化感自毒物质是导致农作物连作障碍的重要诱因之一。枸杞作为多年生木本植物,连作障碍突出,成因复杂。为揭示酚酸类物质对枸杞的化感自毒作用,研究了22种常见酚酸对枸杞种子萌发的抑制作用,并运用比较分子场分析方法(CoMFA)进行构效关系研究。结果发现,酚酸化合物对枸杞种子萌发作用的半抑制浓度IC_(50)范围为39.94—115.97 mg/L;CoMFA结果表明,运用立体场和静电场能较好解释该类化合物对枸杞种子萌发的化感自毒特征,苯环1/2/3/4位具有大体积、1位具有强负电性取代基的酚酸化合物,对枸杞种子萌发具有较强抑制作用。研究结果可为客观评价酚酸类化合物对枸杞的化感自毒作用提供基础数据支撑。     

Synthesis, 3D-QSAR, and docking studies of 1-phenyl-1H-1,2,3-triazoles as selective antagonists for beta3 over alpha1beta2gamma2 GABA receptors

A series of 16 1-phenyl-1H-1,2,3-triazoles with substituents at both the 4- and 5-positions of the triazole ring were synthesized, and a total of 49 compounds, including previously reported 4- or 5-monosubstituted analogues, were examined for their ability to inhibit the specific binding of [(3)H]4'-ethynyl-4-n-propylbicycloorthobenzoate (EBOB), a non-competitive antagonist, to human homo-oligomeric beta3 and hetero-oligomeric alpha1beta2gamma2 gamma-aminobutyric acid (GABA) receptors. Among all tested compounds, the 4-n-propyl-5-chloromethyl analogue of 1-(2,6-dichloro-4-trifluoromethylphenyl)-1H-1,2,3-triazole showed the highest level of affinity for both beta3 and alpha1beta2gamma2 receptors, with K(i) values of 659pM and 266nM, respectively. Most of the tested compounds showed selectivity for beta3 over alpha1beta2gamma2 receptors. Among all 1-phenyl-1H-1,2,3-triazoles, the 4-n-propyl-5-ethyl analogue exhibited the highest (>1133-fold) selectivity, followed by the 4-n-propyl-5-methyl analogue of 1-(2,6-dibromo-4-trifluoromethylphenyl)-1H-1,2,3-triazole with a >671-fold selectivity. The 2,6-dichloro plus 4-trifluoromethyl substitution pattern on the benzene ring was found to be important for the high affinity for both beta3 and alpha1beta2gamma2 receptors. Comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) provided similar contour maps, revealing that an electronegative substituent at the 4-position of the benzene ring, a compact, hydrophobic substituent at the 4-position of the triazole ring, and a small, electronegative substituent at the 5-position of the triazole ring play significant roles for the high potency in beta3 receptors. Molecular docking studies suggested that the putative binding sites for 1-phenyl-1H-1,2,3-triazole antagonists are located in the channel-lining 2'-6' region of the second transmembrane segment of beta3 and alpha1beta2gamma2 receptors. A difference in the hydrophobic environment at the 2' position might underlie the selectivity of 1-phenyl-1H-1,2,3-triazoles for beta3 over alpha1beta2gamma2 receptors. The compounds that had high affinity for beta3 receptors with homology to insect GABA receptors showed insecticidal activity against houseflies with LD(50) values in the pmol/fly range. The information obtained in the present study should prove helpful for the discovery of selective insect control chemicals.     

Synthesis and insecticidal activity of novel N-oxydihydropyrrole derivatives with a substituted spirocyclohexyl group

A series of 5-spirocyclohexyl-3-(2,6-dimethylphenyl)-1,5-dihydro-2H-pyrrol-2-one derivatives (3) with various substituents on the spirocyclohexyl ring was synthesized and evaluated for its insecticidal activity against the aphid, Myzus persicae. Substituents at the 1- and 4-positions of the dihydropyrrole ring were also varied to optimize the activity. An investigation of the structure-activity relationship revealed that methoxy, alkoxyalkoxy, ethylenedioxy and methoxyimino groups were favorable as substituents at the 4-position of the spirocyclohexyl ring. The activity was optimized by the respective substitution of a methoxy or methoxymethoxy moiety and cyclopropylcarbonyloxy group at the 1- and 4-positions of the dihydropyrrole ring.     

Design, synthesis and biological evaluation of novel 2-methylpyrimidine-4-ylamine derivatives as inhibitors of Escherichia coli pyruvate dehydrogenase complex E1

As potential inhibitors of Escherichia coli pyruvate dehydrogenase complex E1 (PDHc E1), a series of novel 2-methylpyrimidine-4-ylamine derivatives were designed based on the structure of the active site of PDHc E1 and synthesized using 'click chemistry'. Their inhibitory activity in vitro against PDHc E1 and fungicidal activity were examined. Some of these compounds such as 3g, 3l, 3n, 3o, and 5b demonstrated to be effective inhibitors of PDHc E1 from E. coli and exhibited antifungal activity. SAR analysis indicated that both, the inhibitory potency against E. coli PDHc E1 and the antifungal activity of title compounds, could be increased greatly by optimizing substituent groups in the compounds. The structures of substituent group in 5-position on the 1,2,3-triazole and 4-position on the benzene ring in title compounds were found to play a pivotal role in both above-mentioned biological activities. Amongst all the compounds, compound 5b with iodine in the 5-position of 1,2,3-triazole and with nitryl group in the 4-position of benzene ring acted as the best inhibitor against PDHc E1 from E. coli. It was also found to be the most effective compound with higher antifungal activity against Rhizoctonia solani and Botrytis cinerea at the dosage of 100 μg mL(-1). Therefore, in this study, compound 5b was used as a lead compound for further optimization.     

Structure-activity relationship of natural and synthetic coumarins inhibiting the multidrug transporter P-glycoprotein

A set of 32 natural and synthetic coumarins were tested in order to evaluate their activity on human leukemic cells (K562/R7) overexpressing P-glycoprotein (P-gp). Their ability to reduce the P-gp-mediated drug efflux of daunorubicin out of cells was evaluated at 10 microM. Four natural compounds, previously isolated from Calophyllum dispar (Clusiaceae) and substituted by a common alpha-(hydroxyisopropyl)dihydrofuran moiety, exhibited a significant inhibitory effect on P-gp when compared to the positive control cyclosporin A. A 3D-quantitative structure-activity relationship (3D-QSAR) analysis of the coumarins was performed using the biological results obtained by comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) of P-gp. Results showed a favorable electrostatic and steric volume, like the alpha-(hydroxyisopropyl)dihydrofuran moiety, beside C(5)-C(6) or C(7)-C(8) positions. In addition, the analysis revealed an important hydrophobic, neutral charge group, like phenyl, in position C(4) on the coumarinic ring.     

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

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

Structure-activity relationships of a series of tariquidar analogs as multidrug resistance modulators

Tariquidar (XR9576) analogs, modulators of cancer multidrug resistance (MDR), were subjected to QSAR and 3D-QSAR analyses. The structural features contributing to anti-MDR activity were identified by the Free-Wilson analysis and pharmacophore search using Hoechst 33342 as a template. 3D-QSAR CoMFA and CoMSIA models were derived and tested. The best models yielded an external predictivity of 0.66-0.75 squared correlation coefficient and outlined HB-acceptor, steric, and hydrophobic fields as the most important 3D properties. On the basis of the QSAR and 3D-QSAR analyses it was suggested that the strong inhibitory potency of the compounds studied is related to the presence of a bulky aromatic ring system with a 3rd positioned heteroatom toward the anthranilamide nucleus in the opposite end of the tetrahydroquinoline group. The results can help in directing the rational design of new generations of potent P-glycoprotein MDR modulators.     

CoMFA and HQSAR studies on 6,7-dimethoxy-4-pyrrolidylquinazoline derivatives as phosphodiesterase10A inhibitors

Phosphodiesterase10A (PDE10A) is an important enzyme with diverse biological actions in intracellular signaling systems, making it an emerging target for diseases such as schizophrenia, Huntington's disease, and diabetes mellitus. The objective of the current 3D QSAR study is to uncover some of the structural parameters which govern PDE10A inhibitory activity over PDE3A/B. Thus, comparative molecular field analysis (CoMFA) and hologram quantitative structure-activity relationship (HQSAR) studies were carried out on recently reported 6,7-dimethoxy-4-pyrrolidylquinazoline derivatives as PDE10A inhibitors. The best CoMFA model using atom-fit alignment approach with the bound conformation of compound 21 as the template yielded the steric parameter as a major contributor (nearly 70%) to the observed variations in biological activity. The best CoMFA model produced statistically significant results, with the cross-validated (r(cv)(2)) and conventional correlation (r(ncv)(2)) coefficients being 0.557 and 0.991, respectively, for the 21 training set compounds. Validation of the model by external set of six compounds yielded a high (0.919) predictive value. The CoMFA models of PDE10A and PDE3A/B activity were compared in order to address the selectivity issue of these inhibitors. The best HQSAR model for PDE10A was obtained with an r(cv)(2) of 0.704 and r(ncv)(2) of 0.902 using atoms, bonds, connections, chirality, donor, and acceptor as fragment distinction and default fragment size of 4-7 with three components for the 21 compounds. The HQSAR model predicted the external test-set of compounds well since a good agreement between the experimental and predicted values was verified. Taken together, the present QSAR models were found to accurately predict the PDE10A inhibitory activity of the test-set compounds and to yield reliable clues for further optimization of the quinazoline derivatives in the dataset.     

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.     

Comparative quantitative structure-activity study of radical scavengers

Classic and three-dimensional (3-D) QSAR analyses of 13 radical scavengers (1-13) were performed to derive two classic, two Apex-3-D and one comparative field analysis (CoMFA) models. Two classical models with predictive cross-validated r2 (Q2) over 0.96 indicated that the activity was attributed to the electronic COH and ELUMO, steric molar refractivity (MR) and lipophilic log P. Three-dimensional quantitative structure-activity relationship (3-D-QSAR) studies were performed by 3-D pharmacophore generation (Apex-3-D) and CoMFA techniques. For Apex-3-D studies, two best models with high Q2 (0.94 and 0.97) were yielded. Structural properties contributing to the activity were not only lipophilic but also the optimum steric property and geometry of side-chain composition. For CoMFA studies, the sp3 C(+1) probe provided the best Q2 of 0.79 with steric and electrostatic contributions of 42.3 and 57.7%, respectively. The activity of four new compounds (14-17) not included in the derivation were predicted with these models. Although the derived models were from limited data, the statistic relation was predictive. The linear correlations between the experimental IC50 values and the predicted values from classical and Apex-3-D models were found to be high and significant. The predicted activity of 17 from CoMFA was much lower than the experimental value; this deviation occurred according to the missing of hydrophobic field in standard CoMFA study. In vitro and ex vivo antilipid peroxidation in mouse brain and ESR studies of 14-17 were investigated for the radical-scavenging ability. The difference between the in vitro results, antilipid peroxidation and electron spin resonance (ESR) and ex vivo results in coumarin series was found. Thus, other properties for good bioavailability besides log P should also be taken into consideration.     

The new generation dihydropyridine type calcium blockers, bearing 4-phenyl oxypropanolamine, display alpha-/beta-adrenoceptor antagonist and long-acting antihypertensive activities

A new series of dihydropyridine derivatives, bearing oxypropanolamine moiety on phenyl ring at the 4-position of the dihydropyridine base, were prepared. Oxypropanolamine was synthesized by replacing the phenolic OH of vanillin or other compounds, having a phenyl aldehyde group, with epichlorohydrin, followed by cleavaging the obtained epoxide compounds with tert-butylamine, n-butylamine or 2-methoxy-1-oxyethylamino benzene (guaiacoxyethylamine), respectively. Obtained various oxypropanolamine compounds, still remaining a phenyl aldehyde moiety, were then performed by Hantzsch condensation reaction with methylacetoacetate or ethylacetoacetate, respectively, to give our new series of dihydropyridine linked with the 4-phenyl ring. These compounds were evaluated for inotropic, chronotropic, and aorta contractility that associated with calcium channel and adrenoceptor antagonist activities. Dihydropyridine derivatives that with oxypropanolamine side chain on their 4-phenyl ring associated alpha-/beta-adrenoceptor blocking activities created a new family of calcium entry and the third generation beta-adrenoceptor blockers. Optimizing this research to obtain more potent alpha-/beta-adrenoceptor blocking and long-acting antihypertensive oxypropanolamine on the 4-phenyl ring of dihydropyridine series compounds was thus accomplished and classified as third generation dihydropyridine type calcium channel blockers, in comparison with previous short-acting type nifedipine and long-acting type amlodipine. We concluded that compounds 1a, 1b and 1g showed not only markedly high calcium-antagonistic activity but also the highest antihypertensive effect; compounds 1b, 1c, 1f, 1g, 1i and 1j induced sustained antihypertensive effects are major and attributed to their calcium entry and alpha-adrenoceptor blocking activities in the blood vessel due to their introduction of 2-methoxy, 1-oxyethylamino benzene moiety in the side chain on the 4-phenyl ring of dihydropyridine. Bradycardiac effects of all the compounds 1a-1j resulted from calcium entry and beta-adrenoceptor blocking, which attenuate the sympathetic activation-associated reflex tachycardia in the heart. We selected compound 1b as candidate compound for further pharmacological and pre-clinical evaluation studies.     

Novel diphenylalkyl piperazine derivatives with high affinities for the dopamine transporter

The novel diphenyl piperazine derivatives containing the phenyl substituted aminopropanol moiety, including 1-[4,4-bis(4-fluorophenyl)butyl]-4-[2-hydroxy-3-(phenylamino)propyl]piperazine 1, which were modified at the connective between the diphenyl and piperazine moieties, have been found to be potent dopamine uptake inhibitors. To study the further structure-activity relationship (SAR) of these compounds, a new series was synthesized, with modifications at the 2-hydroxy-3-phenylaminopropyl moiety of 1. The series was evaluated for dopamine transporter (DAT) binding affinity with [3H]GBR12935 in rat striatal membranes. Most of the compounds showed moderate to high DAT binding affinities and some were approximately equivalent in activity to compound 1 or GBR12909 as a dopamine uptake inhibitor, with IC(50) values of nanomolar range. The SAR suggested that on exhibiting a potent interaction with the DAT, there is probably a steric limitation around the benzene ring of the phenylamino moiety of 1, allowing only small-sized substituents with the exception of basic moieties at the 4-position. In addition, the SAR at the 3-amino-2-propanol moiety of 1 suggested that either the nitrogen atom with an electron donating substituent or the unsubstituted nitrogen atom and also the hydroxy group are desirable for elicitation of a potent DAT binding affinity.