Local intersection volume: a new 3D descriptor applied to develop a 3D-QSAR pharmacophore model for benzodiazepine receptor ligands

QSPR model for Caco-2 cell permeability prediction using a combination of HQPSO and dual-RBF neural network

The aim of this study is to establish a promising QSPR model for the Caco-2 permeability prediction.

Ligand-based 3D-QSAR studies of diaryl acyl-sulfonamide analogues as human umbilical vein endothelial cells inhibitors stimulated by VEGF

Diaryl acylsulfonamide derivatives were reported as Human Umbilical Vein Endothelial cell (HUVEC) inhibitors, stimulated by Vascular Endothelial Growth Factor (VEGF). VEGF has angiogenic property to cause colorectal cancer. A ligand-based 3D-QSAR technique was carried out on diaryl acylsulfonamide derivatives by using Comparative Molecular Field Analysis (CoMFA) studies to find relations between biological activities of inhibitors and their structures. In absence of binding mechanism for the ligand with VEGF receptor, current study hopes to shed some light on the inhibition mechanism of the ligands with HUVEC. 3D-QSAR technique was applied to a set of fifty ligands in order to facilitate the design of more potent inhibitors. However, the maximum cross-validated correlation coefficient value was found to be 0.417. The value is relatively low when compared to the usual acceptable cross-validated correlation coefficient, but no further improvements were observed by applying different available options. Therefore, the final model was used for further analysis. Additionally, the resulted CoMFA model was validated by an external set of 10 compounds yielding surprisingly, a satisfactory correlation coefficient value (r(2) (pred) ) 0.80. Moreover, the analysis of the individual generated 3D contours helped in understanding the possible structural modifications of molecules to improve the inhibitory potency.

Benzopyrazine derivatives: A novel class of growth factor receptor bound protein 7 antagonists

Growth factor receptor bound protein 7 (Grb7) is an adapter protein that functions as a downstream effector of growth factor mediated signal transduction. Over-expression of Grb7 has been implicated in a variety of cancers such as breast, blood, pancreatic, esophageal, and gastric carcinomas. Inhibition of Grb7 has been shown to reduce the migratory and proliferative potential of these cancers, making it an attractive therapeutic target. Starting with a known peptide antagonist, the present work reports the application of a succession of computational ligand design tools comprising a ligand shape based similarity search, molecular docking and a 2D-similarity search to identify small molecular antagonists of the Grb7-SH2 domain from the NCI chemical database. Binding to the Grb7-SH2 domain was then experimentally tested using melting point shift assays and isothermal titration calorimetry. Overall, a total of 11 benzopyrazine based small molecular antagonists were identified with affinity for the Grb7-SH2 domain. Representative compounds tested using ITC were revealed to possess moderate binding affinity in the low micromolar range. Finally, the lead compound (NSC642056) was found to reduce the growth of a Grb7-expressing breast cancer cell line with an IC(50) of 86μM. It is expected that the identified antagonists will be useful additions to further explore the function of Grb7 and for the development of inhibitors with therapeutic potential.

Docking of 1,4-benzodiazepines in the alpha1/gamma2 GABA(A) receptor modulator site

Positive allosteric modulation of the GABA(A) receptor (GABA(A)R) via the benzodiazepine recognition site is the mechanism whereby diverse chemical classes of therapeutic agents act to reduce anxiety, induce and maintain sleep, reduce seizures, and induce conscious sedation. The binding of such therapeutic agents to this allosteric modulatory site increases the affinity of GABA for the agonist recognition site. A major unanswered question, however, relates to how positive allosteric modulators dock in the 1,4-benzodiazepine (BZD) recognition site. In the present study, the X-ray structure of an acetylcholine binding protein from the snail Lymnea stagnalis and the results from site-directed affinity-labeling studies were used as the basis for modeling of the BZD binding pocket at the alpha(1)/gamma(2) subunit interface. A tethered BZD was introduced into the binding pocket, and molecular simulations were carried out to yield a set of candidate orientations of the BZD ligand in the binding pocket. Candidate orientations were refined based on known structure-activity and stereospecificity characteristics of BZDs and the impact of the alpha(1)H101R mutation. Results favor a model in which the BZD molecule is oriented such that the C5-phenyl substituent extends approximately parallel to the plane of the membrane rather than parallel to the ion channel. Application of this computational modeling strategy, which integrates site-directed affinity labeling with structure-activity knowledge to create a molecular model of the docking of active ligands in the binding pocket, may provide a basis for the design of more selective GABA(A)R modulators with enhanced therapeutic potential.

The importance of discerning shape in molecular pharmacology

Shape is a fundamentally important molecular feature that often determines the fate of a compound in terms of molecular interactions with preferred and non-preferred biological targets. Complementarity of binding in small-molecule-protein, peptide-receptor, antigen-antibody and protein-protein interactions is the key to life and survival and also to targeting molecules with bioactivity. We review the application of shape in various biological systems such as substrate recognition, ligand specificity or selectivity and antibody recognition in the context of computational methods such as docking, quantitative structure-activity relationships, classification models and similarity-search algorithms. These in silico pharmacology methods have recently demonstrated the importance and applicability of determining molecular shape in drug discovery, virtual screening and predictive toxicology. The results from recently published studies show that shape and shape-based descriptors are at least as useful as other traditional molecular descriptors.

Pharmacophore mapping of arylbenzothiophene derivatives for MCF cell inhibition using classical and 3D space modeling approaches

Considering the worth of developing non-steroidal estrogen analogs, the present study explores the pharmacophore features of arylbenzothiophene derivatives for inhibitory activity to MCF-7 cells using classical QSAR and 3D space modeling approaches. The analysis shows that presence of phenolic hydroxyl group and ketonic linkage in the basic side chain of 2-arylbenzothiophene core of raloxifene derivatives are crucial. Additionally piperidine ring connected through ether linkage is favorable for inhibition of breast cancer cell line. These features for inhibitory activity are also highlighted through 3D space modeling approach that explored importance of critical inter features distance among HB-acceptor lipid, hydrophobic and HB-donor features in the arylbenzothiophene scaffold for activity.

Evaluation of anticonvulsant and analgesic effects of benzyl- and benzhydryl ureides

The anticonvulsant effects of benzyl- and benzhydryl ureides in mice models of seizures (maximal electroshock seizure test, pentylenetetrazol test, picrotoxin-induced seizure test) and the influence on spontaneous locomotor activity has been assessed. Furthermore, the analgesic effect of ureide derivatives was studied in the hot-plate test in mice. Selected compounds were investigated for their in vitro interaction with adenosine receptors as well as the benzodiazepine binding site of GABA(A) receptors. This study demonstrated the strong anticonvulsant activity of several ureides in electrically or chemically induced seizure models, and structure-activity relationships were discussed. 1-Benzyl-3-butyrylurea (9) was found to be equipotent to ethosuximide in the pentylenetetrazol test with regard to the number of attacks as well as the time of the onset of seizures. The ureide 9 also revealed the highest protective activity against seizures in the other models, maximal electroshock seizure and picrotoxin test. Moreover, 1-benzyl-3-butyrylurea was not neurotoxic at doses up to 200 mg/kg. Benzylureides 8-10 showed affinity to the adenosine A1 receptors at low micromolar concentrations. However, the apparent anticonvulsant activity in different seizure models does not appear to result from direct activation of adenosine A1 receptors or GABA(A) receptors, respectively. In the hot-plate test, the majority of investigated compounds exhibited analgesic activity. Again, compound 9 was superior to the other substances investigated, suggesting a potential therapeutic value of that ureide derivative.

Refinement and Use of the Approximate Similarity in QSAR Models for Benzodiazepine Receptor Ligands

Several considerations for refining the approximate similarity measurements have been introduced in this paper: the use of topological invariants for the calculation of similarity indexes and the development of new similarity correction processes. The quality of the new similarity measurements obtained with the proposed methods has permitted the development of fast, cheap, and simple quantitative structure-activity relationship models for the prediction of biological activities of nonbenzodiazepine gamma-aminobutyric acid(A)/benzodiazepine receptor ligands (58 compounds). Internal and external validations were carried out for the approximate similarity matrices computed using different approaches. Satisfactory results which compare reasonably well with a 3D approach were obtained: Q2= 0.65 and standard error in cross validation SECV= 0.83 for the training stage; r = 0.79 and error in external prediction = 0.82 for the test step. In addition, the method proposed was compared with other topological approaches based on constitutional similarity and on fingerprints. Satisfactory results were obtained.

Free-energy force-field three-dimensional quantitative structure–activity relationship analysis of a set of p38-mitogen activated protein kinase inhibitors

The p38-mitogen-activated protein kinases (p38-MAPKs) belong to a family of serine-threonine kinases activated by pro-inflammatory or stressful stimuli that are known to be involved in several diseases. Their biological importance, related to the release of inflammatory pro-cytokines such as tumor necrosis factor-alpha (TNF-alpha) and interleukin-1 (IL-1), has generated many studies aiming at the development of selective inhibitors for the treatment of inflammatory diseases. In this work, we developed receptor-based three dimensional (3D) quantitative structure-activity relationship (QSAR) models for a series of 33 pyridinyl imidazole compounds [Liverton et al. (1999) 42:2180], using a methodology named free-energy force-field (FEFF) [Tokarski and Hopfinger (1997) 37:792], in which scaled intra- and intermolecular energy terms of the Assisted Model Building Energy Refinement (AMBER) force field combined with a hydration-shell solvation model are the independent variables used in the QSAR studies. Multiple temperature molecular-dynamics simulations (MDS) of ligand-protein complexes and genetic-function approximation (GFA) were employed using partial least squares (PLS) as the fitting functions to develop FEFF-3D-QSAR models for the binding process. The best model obtained in the FEFF-3D-QSAR receptor-dependent (RD) method shows the importance of the van der Waals energy change upon binding and the electrostatic energy in the interaction of ligands with the receptor. The QSAR equations described here show good predictability and may be regarded as representatives of the binding process of ligands to p38-MAPK. Additionally, we have compared the top FEFF-3D-QSAR model with receptor independent (RI) 4D-QSAR models developed in a recent study [Romeiro et al. (2005) 19:385].

Design, synthesis, and pharmacological evaluation of new neuroactive pyrazolo[3,4-b]pyrrolo[3,4-d]pyridine derivatives with in vivo hypnotic and analgesic profile

The present study describes the synthesis and pharmacological profiles of four novel pyrazolo[3,4-b]pyrrolo[3,4-d]pyridine derivatives 2-5, which were structurally designed by using the sedative and analgesic drug zolpidem 1 as lead compound. The heterotricyclic system present in the target compounds 2-5 was constructed in good yields, exploiting a regioselective hetero Diels-Alder reaction of the key azabutadiene derivative 7 and functionalized N-phenylmaleimides 9-12. Additionally, we identified that 1-methyl-7-(4-nitrophenyl)-3-phenyl-3,6,7,8-tetrahydropyrazolo[3,4-b]pyrrolo[3,4-d]pyridine-6,8-dione derivative (LASSBio-873, 5) presented not only the most potent ability to promote sedation, which was similar to that induced by the standard benzodiazepine drug midazolam, but also potent central antinociceptive effect.

New optimized piperamide analogues with potent in vivo hypotensive properties

We describe herein the structural optimization of new piperamide analogues, designed from two natural prototypes, piperine 1 and piperdardine 2, obtained from Piper tuberculatum Jacq. (Piperaceae). Molecular modeling studies using semiempirical AM1 method were made in order to establish rational modifications to optimize them by molecular simplification. The targeted compounds (10) and (11) were respectively obtained using benzaldehyde (12) and para-anisaldehyde (13) as starting materials. 1H NMR spectra showed that the target compounds were diastereoselectively obtained as the (E)-isomer, the same geometry of the natural prototypes. These new synthetic amides presented significant hypotensive effects in cardiovascular essays using in vivo methodologies. Compound 11 (N-[5-(4'-methoxyphenyl)-2(E)-pentenoyl]thiomorpholine) showed a potency 10,000 times greater than its prototype 5, evidencing an optimization of the molecular architecture for this class of hypotensive drug candidates.

A novel 3D-QSAR comparative molecular field analysis (CoMFA) model of imidazole and quinazolinone functionalized p38 MAP kinase inhibitors

In this study we describe a new comparative molecular field analysis (CoMFA) model of dihydroquinazolinone and tetrasubstituted imidazole compounds with p38 MAPK inhibitory activity. A series of 51 (a training set of 40 and a test set of 11) dihydroquinazolinone [Bioorg. Med. Chem. Lett. 2003, 13, 277.] and tetrasubstituted imidazole [J. Med. Chem. 1999, 42, 2180.] derivatives known as p38 mitogen-activated protein kinase (p38 MAPK) selective inhibitors was studied by quantitative structure-activity relationship (3D-QSAR) analysis using comparative molecular field analysis. The 3D-QSAR models were generated and evaluated by a scheme that combines a genetic algorithm (GA) optimization with partial least squares (PLS) regression and by crossvalidation using the leave-one-out technique. The model was able to efficiently predict the activities of the compounds of the test set, suggesting that it can be used for the planning of new p38 MAPK inhibitor candidates useful to treat chronic inflammatory states.

LIV-3D-QSAR models for PGI2 receptor ligands using multiple conformations

A new 3D descriptor, the local intersection volume (LIV), was developed by our group and applied to the construction of 3D-QSAR models for ligands of the PGI(2) receptor (IP). The target compounds are a set of 42 aromatic heterocyclic derivatives [Meanwell et al., J. Med. Chem. 36 (1993), 3884], which show agonist activities in the IP receptor and are inhibitors of platelet aggregation. The LIV-3D-QSAR models were obtained through the analysis of 30% of the generated conformations for each compound, using a combined Genetic Algorithm (GA) and Partial Least Square (PLS) approach [Rogers and Hopfinger, J. Inf. Comput. Sci. 34 (1994) 854]. Statistically, Model 3 is the best as well as the most comprehensive in a mechanistic sense. Furthermore, it can be applied to design new IP ligands.