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).     

Comparative molecular field analysis of non-steroidal aromatase inhibitors related to fadrozole

Summary A series of non-steroidal inhibitors of aromatase, structurally related to fadrozole (2), was investigated with the aim of developing a 3D QSAR model using the Comparative Molecular Field Analysis (CoMFA) technique. The alignment of the molecules was performed following two approaches (atom-by-atom and field fit), both starting from an initial hypothesis of superimposition of fadrozole to a steroidal inhibitor (3). From a number of CoMFA models built with different characteristics, one was recognized as the most statistically relevant; this one is discussed in detail. The features of the 3D QSAR model are consistent with those of other 3D and QSAR models of aromatase and its inhibitors.     

光系统Ⅱ抑制剂: 顺式氰基丙烯酸酯化合物的 三维构效关系的电子结构研究

以距离比较法所获得的顺式氰基丙烯酸酯化合物的活性构象为模板,对39个该类化合物采用比较分子场方法进行了三维构效关系的研究。结果表明,所获得的药效团模型具有很好的预测能力。同时采用量子化学的方法对活性构象模板分子电子结构作了讨论。     

PharmID: pharmacophore identification using Gibbs sampling

The binding of a small molecule to a protein is inherently a 3D matching problem. As crystal structures are not available for most drug targets, there is a need to be able to infer from bioassay data the key binding features of small molecules and their disposition in space, the pharmacophore. Fingerprints of 3D features and a modification of Gibbs sampling to align a set of known flexible ligands, where all compounds are active, are used to discern possible pharmacophores. A clique detection method is used to map the features back onto the binding conformations. The complete algorithm is described in detail, and it is shown that the method can find common superimposition for several test data sets. The method reproduces answers very close to the crystal structure and literature pharmacophores in the examples presented. The basic algorithm is relatively fast and can easily deal with up to 100 compounds and tens of thousands of conformations. The algorithm is also able to handle multiple binding mode problems, which means it can superimpose molecules within the same data set according to two different sets of binding features. We demonstrate the successful use of this algorithm for multiple binding modes for a set of D2 and D4 ligands.     

A self-organizing algorithm for molecular alignment and pharmacophore development

We present a method for simultaneous three-dimensional (3D) structure generation and pharmacophore-based alignment using a self-organizing algorithm called Stochastic Proximity Embedding (SPE). Current flexible molecular alignment methods either start from a single low-energy structure for each molecule and tweak bonds or torsion angles, or choose from multiple conformations of each molecule. Methods that generate structures and align them iteratively (e.g., genetic algorithms) are often slow. In earlier work, we used SPE to generate good-quality 3D conformations by iteratively adjusting pairwise distances between atoms based on a set of geometric rules, and showed that it samples conformational space better and runs faster than earlier programs. In this work, we run SPE on the entire ensemble of molecules to be aligned. Additional information about which atoms or groups of atoms in each molecule correspond to points in the pharmacophore can come from an automatically generated hypothesis or be specified manually. We add distance terms to SPE to bring pharmacophore points from different molecules closer in space, and also to line up normal/direction vectors associated with these points. We also permit pharmacophore points to be constrained to lie near external coordinates from a binding site. The aligned 3D molecular structures are nearly correct if the pharmacophore hypothesis is chemically feasible; postprocessing by minimization of suitable distance and energy functions further improves the structures and weeds out infeasible hypotheses. The method can be used to test 3D pharmacophores for a diverse set of active ligands, starting from only a hypothesis about corresponding atoms or groups.     

QSAR modeling of in vitro inhibition of cytochrome P450 3A4

We report the QSAR modeling of cytochrome P450 3A4 (CYP3A4) enzyme inhibition using four large data sets of in vitro data. These data sets consist of marketed drugs and drug-like compounds all tested in four assays measuring the inhibition of the metabolism of four different substrates by the CYP3A4 enzyme. The four probe substrates are benzyloxycoumarin, testosterone, benzyloxyresorufin, and midazolam. We first show that using state-of-the-art QSAR modeling approaches applied to only one of these four data sets does not lead to predictive models that would be useful for in silico filtering of chemical libraries. We then present the development and the testing of a multiple pharmacophore hypothesis (MPH) that is formulated as a conceptual extension of the traditional QSAR approach to modeling the promiscuous binding of a large variety of drugs to CYP3A4. In the simplest form, the MPH approach takes advantage of the multiple substrate data sets and identifies the binding of test compounds as either proximal or distal relative to that of a given substrate. Application of the approach to the in silico filtering of test compounds for potential inhibitors of CYP3A4 is also presented. In addition to an improvement in the QSAR modeling for the inhibition of CYP3A4, the results from this modeling approach provide structural insights into the drug-enzyme interactions. The existence of multiple inhibition data sets in the BioPrint database motivates the original development of the concept of a multiple pharmacophore hypothesis and provides a unique opportunity for formulating alternative strategies of QSAR modeling of the inhibition of the in vitro metabolism of CYP3A4.     

3D-QSAR studies on indole and 7-azoindole derivatives as ROCK-2 inhibitors: An integrative computational approach

Rho Kinases (ROCK) has been found to regulate a wide range of fundamental cell functions such as contraction, motility, proliferation, and apoptosis. Recent experiments have defined new functions of ROCKs in cells, including centrosome positioning and cell-size regulation, which might contribute to various physiological and pathological states. In this study, we have performed pharmacophore modeling and 3D QSAR studies on a series of 36 indoles and 7-azoindoles derivatives as ROCK2 inhibitors to elucidate the structural variations with their inhibitory activities. Ligand based CoMFA and CoMSIA models were generated based on three different alignment methods such as systematic search, simulated annealing and pharmacophore. A total of 15 CoMFA models and 27 CoMSIA were generated using different alignments. One model from each alignment is selected based on the statistical values. Contour maps of the selected models were compared, analysed and reported. The 3D QSAR study revealed that electro positive group linked to the methoxy-benzene ring position of the structure will enhance the biological activity and bulkier substitutions are preferred in the methyl dihydroindole region. Also, it is found that the hydrogen bond donor substituted at the R₁ position enhances the inhibitory activity. In future, this study would give proper guidelines to further enhance the activity of novel inhibitors for ROCK2.     

一种新的De Novo生物活性分子设计方法

由药效团进行虚拟活性结构生成与3D-QSAR模型相结合,筛选出有前途的结构多样性的化合物,并从中寻找活性先导化合物,是一种新的分子设计方法。采用这种方法对抗小麦赤霉病类含氟农药进行了研究,共生成了53个虚拟活性结构,通过3D-QSAR模型筛选出其中10个活性较高的结构,在活性最高的化合物基础上进行了结构修饰,得到了活性更高且毒性较低的理想化合物。研究结果表明这种方法能突破原模型化合物结构模式的局限,可以找到结构新颖的活性先导化合物,是一种非常有前途的分子设计方法,而且具有较高的筛选效率。     

Pharmacophore mapping and electronic feature analysis for a series of nitroaromatic compounds with antitubercular activity

A five point pharmacophore was generated using PHASE for a series of nitroaromatic compounds and their congeners as antitubercular agents. The generated pharmacophore yielded significant 3D‐QSAR model with r² of 0.890 for a training set of 92 molecules. The model also showed excellent predictive power with correlation coefficient Q² of 0.857 for a test set of 31 compounds. The pharmacophore indicated that presence of a nitro group, a piperazine moiety, one aromatic ring feature and two acceptor features are necessary for potent antitubercular activity. The pharmacophore was supported by electronic property analysis using density functional theory (DFT) at B3LYP/3–21*G level. Molecular electrostatic profile of the compounds was consistent with the generated pharmacophore model, particularly appearance of localized negative potential regions near both the oxygen atoms of nitro group extending laterally to the isoxazole ring system/amide bond in the most active compounds. Calculated data further revealed that all active compounds have smaller LUMO energies located over the nitro group, furan ring, and isoxazole ring/amide bond attached to it. Higher negative values of LUMO energies concentrated over the nitro group are indicative of the electron acceptor capacity of the compounds, suggesting that these compounds are prodrugs and must be activated by TB‐nitroreductase. The results obtained from this study should aid in efficient design and development of nitroaromatic compounds as antitubercular agents. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

Pharmacophore generation,atom-based 3D-QSAR and molecular dynamics simulation analyses of pyridine-3-carboxamide-6-yl-urea analogues as potential gyrase B inhibitors

DNA gyrase subunit B (GyrB) is an attractive drug target for the development of antibacterial agents with therapeutic potential. In the present study, computational studies based on pharmacophore modelling, atom-based QSAR, molecular docking, free binding energy calculation and dynamics simulation were performed on a series of pyridine-3-carboxamide-6-yl-urea derivatives. A pharmacophore model using 49 molecules revealed structural and chemical features necessary for these molecules to inhibit GyrB. The best fitted model AADDR.13 was generated with a coefficient of determination (r²) of 0.918. This model was validated using test set molecules and had a good r² of 0.78. 3D contour maps generated by the 3D atom-based QSAR revealed the key structural features responsible for the GyrB inhibitory activity. Extra precision molecular docking showed hydrogen bond interactions with key amino acid residues of ATP-binding pocket, important for inhibitor binding. Further, binding free energy was calculated by the MM-GBSA rescoring approach to validate the binding affinity. A 10 ns MD simulation of inhibitor #47 showed the stability of the predicted binding conformations. We identified 10 virtual hits by in silico high-throughput screening. A few new molecules were also designed as potent GyrB inhibitors. The information obtained from these methodologies may be helpful to design novel inhibitors of GyrB.     

Triple-layered QSAR studies on substituted 1,2,4-trioxanes as potential antimalarial agents: superiority of the quantitative pharmacophore-based alignment over common substructure-based alignment

This study reports the utilization of three approaches – pharmacophore, CoMFA/CoMSIA and HQSAR studies – to identify the essential structural requirements in 3D chemical space for the modulation of the antimalarial activity of substituted 1,2,4-trioxanes. The superiority of quantitative pharmacophore-based alignment (QuantitativePBA) over global minima energy conformer-based alignment (GMCBA) has been reported in CoMFA and CoMSIA studies. The developed models showed good statistical significance in internal validation (q ², group cross-validation and bootstrapping) and performed very well in predicting the antimalarial activity of test set compounds. Structural features in terms of their steric, electrostatic and hydrophobic interactions in 3D space have been found to be important for the antimalarial activity of substituted 1,2,4-trioxanes. Further, the HQSAR studies based on the same training and test set acted as an additional tool to find the sub-structural fingerprints of substituted 1,2,4-trioxanes for their antimalarial activity. Together, these studies may facilitate the design and discovery of new substituted 1,2,4-trioxanes with potent antimalarial activity.