CoMFA, CoMSIA, and molecular hologram QSAR studies of novel neuronal nAChRs ligands-open ring analogues of 3-pyridyl ether

3-Pyridyl ethers are excellent nAChRs ligands, which show high subtype selectivity and binding affinity to alpha4beta2 nAChR. Although the quantitative structure-activity relationship (QSAR) of nAChRs ligands has been widely investigated using various classes of compounds, the open ring analogues of 3-pyridyl ethers have been less involved in these studies due to the greater flexibility of this kind of molecule. In this study, two three-dimensional QSAR techniques and one two-dimensional QSAR technique were used to correlate the molecular structure with the biological activity of 64 analogues of 3-pyridyl ethers. Three different QSAR models were established. Their performances in the QSAR studies of open ring analogues of 3-pyridyl ethers were evaluated by the statistical values in the corresponding models. All models exhibited satisfactory predictive power. Of these models, the HQSAR behaved optimally in terms of the statistical values with q2=0.845, r2=0.897. Finally, graphic interpretation of three different models provided coincident information about the interaction of the ligand-receptor complex and supplied useful guidelines for the synthesis of novel, potent ligands.     

Binary quantitative structure-activity relationship (QSAR) analysis of estrogen receptor ligands

The use of high throughput screening (HTS) to identify lead compounds has greatly challenged conventional quantitative structure-activity relationship (QSAR) techniques that typically correlate structural variations in similar compounds with continuous changes in biological activity. A new QSAR-like methodology that can correlate less quantitative assay data (i.e., "active" versus "inactive"), as initially generated by HTS, has been introduced. In the present study, we have, for the first time, applied this approach to a drug discovery problem; that is, the study of the estrogen receptor ligands. The binding affinities of 463 estrogen analogues were transformed into a binary data format, and a predictive binary QSAR model was derived using 410 estrogen analogues as a training set. The model was applied to predict the activity of 53 estrogen analogues not included in the training set. An overall accuracy of 94% was obtained.     

Using physics-based pose predictions and free energy perturbation calculations to predict binding poses and relative binding affinities for FXR ligands in the D3R Grand Challenge 2

Computer-aided drug design has become an integral part of drug discovery and development in the pharmaceutical and biotechnology industry, and is nowadays extensively used in the lead identification and lead optimization phases. The drug design data resource (D3R) organizes challenges against blinded experimental data to prospectively test computational methodologies as an opportunity for improved methods and algorithms to emerge. We participated in Grand Challenge 2 to predict the crystallographic poses of 36 Farnesoid X Receptor (FXR)-bound ligands and the relative binding affinities for two designated subsets of 18 and 15 FXR-bound ligands. Here, we present our methodology for pose and affinity predictions and its evaluation after the release of the experimental data. For predicting the crystallographic poses, we used docking and physics-based pose prediction methods guided by the binding poses of native ligands. For FXR ligands with known chemotypes in the PDB, we accurately predicted their binding modes, while for those with unknown chemotypes the predictions were more challenging. Our group ranked #1st (based on the median RMSD) out of 46 groups, which submitted complete entries for the binding pose prediction challenge. For the relative binding affinity prediction challenge, we performed free energy perturbation (FEP) calculations coupled with molecular dynamics (MD) simulations. FEP/MD calculations displayed a high success rate in identifying compounds with better or worse binding affinity than the reference (parent) compound. Our studies suggest that when ligands with chemical precedent are available in the literature, binding pose predictions using docking and physics-based methods are reliable; however, predictions are challenging for ligands with completely unknown chemotypes. We also show that FEP/MD calculations hold predictive value and can nowadays be used in a high throughput mode in a lead optimization project provided that crystal structures of sufficiently high quality are available.     

Possible allosteric interactions of monoindazole-substituted P2 cyclic urea analogues with wild-type and mutant HIV-1 protease

Our ongoing efforts to understand the difference in the binding pattern of HIV-1 protease inhibitor (HIVPI) with the wild-type and mutant HIV-1 protease (HIVPR) and to provide mechanistic insight are continued further. We report here the results of a recent quantitative structure-activity relationship (QSAR) study on monoindazole-substituted P2 analogues of cyclic urea HIVPIs. The QSAR models revealed an inverted parabolic relationship between biological activity and calculated molar refractivity (CMR). That is, biological activity first decreases with increase in CMR and at a certain minimum point (inversion point) it suddenly changes and increases with further increase in CMR. CMR is a measure of volume-dependent-polarizability and is an indication of the polar interactions between ligand and receptor. The results seem to be best rationalized by larger molecules inducing a change in a receptor unit that allows for a new mode of interaction. Similar QSAR models were also observed for the biological activity of these molecules tested against a panel of mutant viruses including mutant strains with single amino acid substitution (I84V), double amino acid substitutions (I84V/V82F), and multiple amino acid changes corresponding to mutations observed in clinical isolates of patients treated with Ritonavir((R)). Interestingly the inversion points for these mutant strains were found larger than for wild-type. The subtle but significant difference in the inversion point indicates change in the shape and size of the binding pocket. Earlier QSAR studies have shown that the correlation of biological activity with an inverted parabola is an indicative of the 'allosteric interaction' of the ligands with the receptor. This report presents a detail analysis of these observations.     

Combining 4D pharmacophore generation and multidimensional QSAR: modeling ligand binding to the bradykinin B2 receptor

We recently reported the development of two receptor-modeling concepts (software Quasar and Raptor) based on multidimensional quantitative structure-activity relationships (QSAR) and allowing for the explicit simulation of induced fit. As the identification of the bioactive configuration of ligand molecules in such studies is all but unambiguous, each compound may be represented by an ensemble of different conformations, orientations, stereoisomers, and protonation states, leading to a 4D data set. In this account, we present a novel technology (software Symposar) allowed to automatically generate a 4D pharmacophore as input for multidimensional QSAR. Symposar aligns ligands utilizing fuzzylike 2D-subfeature mapping and, subsequently, a Monte Carlo search on a 3D similarity grid. The two-step concept (4D pharmacophore generation and quantification of ligand binding by multidimensional QSAR) was applied to 186 compounds binding to the bradykinin B2 receptor. The prediction of their binding affinity by means of the Quasar and Raptor technologies allowed for consensus scoring and generated topologically and quantitatively consistent receptor models. These converged at a cross-validated r2 of 0.752 and 0.815 and yielded a predictive r2 of 0.784 and 0.853 for a test set (for Quasar and Raptor, respectively).     

Development of biologically active compounds by combining 3D QSAR and structure-based design methods

One of the major challenges in computational approaches to drug design is the accurate prediction of the binding affinity of novel biomolecules. In the present study an automated procedure which combines docking and 3D-QSAR methods was applied to several drug targets. The developed receptor-based 3D-QSAR methodology was tested on several sets of ligands for which the three-dimensional structure of the target protein has been solved – namely estrogen receptor, acetylcholine esterase and protein-tyrosine-phosphatase 1B. The molecular alignments of the studied ligands were determined using the docking program AutoDock and were compared with the X-ray structures of the corresponding protein-ligand complexes. The automatically generated protein-based ligand alignment obtained was subsequently taken as basis for a comparative field analysis applying the GRID/GOLPE approach. Using GRID interaction fields and applying variable selection procedures, highly predictive models were obtained. It is expected that concepts from receptor-based 3D QSAR will be valuable tools for the analysis of high-throughput screening as well as virtual screening data     

Identification of Novel Cannabinoid CB2 Receptor Agonists from Botanical Compounds and Preliminary Evaluation of Their Anti-Osteoporotic Effects

As cannabinoid CB2 receptors (CB2R) possess various pharmacological effects—including anti-epilepsy, analgesia, anti-inflammation, anti-fibrosis, and regulation of bone metabolism—without the psychoactive side effects induced by cannabinoid CB1R activation, they have become the focus of research and development of new target drugs in recent years. The present study was intended to (1) establish a double luciferase screening system for a CB2R modulator; (2) validate the agonistic activities of the screened compounds on CB2R by determining cAMP accumulation using HEK293 cells that are stably expressing CB2R; (3) predict the binding affinity between ligands and CB2 receptors and characterize the binding modes using molecular docking; (4) analyze the CB2 receptors–ligand complex stability, conformational behavior, and interaction using molecular dynamics; and (5) evaluate the regulatory effects of the screened compounds on bone metabolism in osteoblasts and osteoclasts. The results demonstrated that the screening system had good stability and was able to screen cannabinoid CB2R modulators from botanical compounds. Altogether, nine CB2R agonists were identified by screening from 69 botanical compounds, and these CB2R agonists exhibited remarkable inhibitory effects on cAMP accumulation and good affinity to CB2R, as evidenced by the molecular docking and molecular dynamics. Five of the nine CB2R agonists could stimulate osteoblastic bone formation and inhibit osteoclastic bone resorption. All these findings may provide useful clues for the development of novel anti-osteoporotic drugs and help elucidate the mechanism underlying the biological activities of CB2R agonists identified from the botanical materials.     

Fuzzy tricentric pharmacophore fingerprints. 2. Application of topological fuzzy pharmacophore triplets in quantitative structure-activity relationships

Topological fuzzy pharmacophore triplets (2D-FPT), using the number of interposed bonds to measure separation between the atoms representing pharmacophore types, were employed to establish and validate quantitative structure-activity relationships (QSAR). Thirteen data sets for which state-of-the-art QSAR models were reported in literature were revisited in order to benchmark 2D-FPT biological activity-explaining propensities. Linear and nonlinear QSAR models were constructed for each compound series (following the original author's splitting into training/validation subsets) with three different 2D-FPT versions, using the genetic algorithm-driven Stochastic QSAR sampler (SQS) to pick relevant triplets and fit their coefficients. 2D-FPT QSARs are computationally cheap, interpretable, and perform well in benchmarking. In a majority of cases (10/13), default 2D-FPT models validated better than or as well as the best among those reported, including 3D overlay-dependent approaches. Most of the analogues series, either unaffected by protonation equilibria or unambiguously adopting expected protonation states, were equally well described by rule- or pKa-based pharmacophore flagging. Thermolysin inhibitors represent a notable exception: pKa-based flagging boosts model quality, although--surprisingly--not due to proteolytic equilibrium effects. The optimal degree of 2D-FPT fuzziness is compound set dependent. This work further confirmed the higher robustness of nonlinear over linear SQS models. In spite of the wealth of studied sets, benchmarking is nevertheless flawed by low intraset diversity: a whole series of thereby caused artifacts were evidenced, implicitly raising questions about the way QSAR studies are conducted nowadays. An in-depth investigation of thrombin inhibition models revealed that some of the selected triplets make sense (one of these stands for a topological pharmacophore covering the P1 and P2 binding pockets). Nevertheless, equations were either unable to predict the activity of the structurally different ligands or tended to indiscriminately predict any compound outside the training family to be active. 2D-FPT QSARs do however not depend on any common scaffold required for molecule superimposition and may in principle be trained on hand of diverse sets, which is a must in order to obtain widely applicable models. Adding (assumed) inactives of various families for training enabled discovery of models that specifically recognize the structurally different actives.