Predicting bioactive conformations and binding modes of macrocycles

Macrocyclic compounds experience increasing interest in drug discovery. It is often thought that these large and chemically complex molecules provide promising candidates to address difficult targets and interfere with protein–protein interactions. From a computational viewpoint, these molecules are difficult to treat. For example, flexible docking of macrocyclic compounds is hindered by the limited ability of current docking approaches to optimize conformations of extended ring systems for pose prediction. Herein, we report predictions of bioactive conformations of macrocycles using conformational search and binding modes using docking. Conformational ensembles generated using specialized search technique of about 70 % of the tested macrocycles contained accurate bioactive conformations. However, these conformations were difficult to identify on the basis of conformational energies. Moreover, docking calculations with limited ligand flexibility starting from individual low energy conformations rarely yielded highly accurate binding modes. In about 40 % of the test cases, binding modes were approximated with reasonable accuracy. However, when conformational ensembles were subjected to rigid body docking, an increase in meaningful binding mode predictions to more than 50 % of the test cases was observed. Electrostatic effects did not contribute to these predictions in a positive or negative manner. Rather, achieving shape complementarity at macrocycle-target interfaces was a decisive factor. In summary, a combined computational protocol using pre-computed conformational ensembles of macrocycles as a starting point for docking shows promise in modeling binding modes of macrocyclic compounds.     

Similarity metrics for ligands reflecting the similarity of the target proteins

In this study we evaluate how far the scope of similarity searching can be extended to identify not only ligands binding to the same target as the reference ligand(s) but also ligands of other homologous targets without initially known ligands. This "homology-based similarity searching" requires molecular representations reflecting the ability of a molecule to interact with target proteins. The Similog keys, which are introduced here as a new molecular representation, were designed to fulfill such requirements. They are based only on the molecular constitution and are counts of atom triplets. Each triplet is characterized by the graph distances and the types of its atoms. The atom-typing scheme classifies each atom by its function as H-bond donor or acceptor and by its electronegativity and bulkiness. In this study the Similog keys are investigated in retrospective in silico screening experiments and compared with other conformation independent molecular representations. Studied were molecules of the MDDR database for which the activity data was augmented by standardized target classification information from public protein classification databases. The MDDR molecule set was split randomly into two halves. The first half formed the candidate set. Ligands of four targets (dopamine D2 receptor, opioid delta-receptor, factor Xa serine protease, and progesterone receptor) were taken from the second half to form the respective reference sets. Different similarity calculation methods are used to rank the molecules of the candidate set by their similarity to each of the four reference sets. The accumulated counts of molecules binding to the reference target and groups of targets with decreasing homology to it were examined as a function of the similarity rank for each reference set and similarity method. In summary, similarity searching based on Unity 2D-fingerprints or Similog keys are found to be equally effective in the identification of molecules binding to the same target as the reference set. However, the application of the Similog keys is more effective in comparison with the other investigated methods in the identification of ligands binding to any target belonging to the same family as the reference target. We attribute this superiority to the fact that the Similog keys provide a generalization of the chemical elements and that the keys are counted instead of merely noting their presence or absence in a binary form. The second most effective molecular representation are the occurrence counts of the public ISIS key fragments, which like the Similog method, incorporates key counting as well as a generalization of the chemical elements. The results obtained suggest that ligands for a new target can be identified by the following three-step procedure: 1. Select at least one target with known ligands which is homologous to the new target. 2. Combine the known ligands of the selected target(s) to a reference set. 3. Search candidate ligands for the new targets by their similarity to the reference set using the Similog method. This clearly enlarges the scope of similarity searching from the classical application for a single target to the identification of candidate ligands for whole target families and is expected to be of key utility for further systematic chemogenomics exploration of previously well explored target families.     

An ontology for pharmaceutical ligands and its application for in silico screening and library design

Annotation efforts in biosciences have focused in past years mainly on the annotation of genomic sequences. Only very limited effort has been put into annotation schemes for pharmaceutical ligands. Here we propose annotation schemes for the ligands of four major target classes, enzymes, G protein-coupled receptors (GPCRs), nuclear receptors (NRs), and ligand-gated ion channels (LGICs), and outline their usage for in silico screening and combinatorial library design. The proposed schemes cover ligand functionality and hierarchical levels of target classification. The classification schemes are based on those established by the EC, GPCRDB, NuclearDB, and LGICDB. The ligands of the MDL Drug Data Report (MDDR) database serve as a reference data set of known pharmacologically active compounds. All ligands were annotated according to the schemes when attribution was possible based on the activity classification provided by the reference database. The purpose of the ligand-target classification schemes is to allow annotation-based searching of the ligand database. In addition, the biological sequence information of the target is directly linkable to the ligand, hereby allowing sequence similarity-based identification of ligands of next homologous receptors. Ligands of specified levels can easily be retrieved to serve as comprehensive reference sets for cheminformatics-based similarity searches and for design of target class focused compound libraries. Retrospective in silico screening experiments within the MDDR01.1 database, searching for structures binding to dopamine D2, all dopamine receptors and all amine-binding class A GPCRs using known dopamine D2 binding compounds as a reference set, have shown that such reference sets are in particular useful for the identification of ligands binding to receptors closely related to the reference system. The potential for ligand identification drops with increasing phylogenetic distance. The analysis of the focus of a tertiary amine based combinatorial library compared to known amine binding class A GPCRs, peptide binding class A GPCRs, and LGIC ligands constitutes a second application scenario which illustrates how the focus of a combinatorial library can be treated quantitatively. The provided annotation schemes, which bridge chem- and bioinformatics by linking ligands to sequences, are expected to be of key utility for further systematic chemogenomics exploration of previously well explored target families.     

A marriage made in torsional space: using GALAHAD models to drive pharmacophore multiplet searches

Pharmacophore multiplets are useful tools for 3D database searching, with the queries used ordinarily being derived from ensembles of random conformations of active ligands. It seems reasonable to expect that their usefulness can be augmented by instead using queries derived from single ligand conformations obtained from aligned ligands. Comparisons of pharmacophore multiplet searching using random conformations with multiplet searching using single conformations derived from GALAHAD (a genetic algorithm with linear assignment for hypermolecular alignment of datasets) models do indeed show that, while query hypotheses based on random conformations are quite effective, hypotheses based on aligned conformations do a better job of discriminating between active and inactive compounds. In particular, the hypothesis created from a neuraminidase inhibitor model was more similar to half of 18 known actives than all but 0.2% of the compounds in a structurally diverse subset of the World Drug Index. Similarly, a model developed from five angiotensin II antagonists yielded hypotheses that placed 65 known antagonists within the top 0.1–1% of decoy databases. The differences in discriminating power ranged from 2 to 20-fold, depending on the protein target and the type of pharmacophore multiplet used.     

Evaluation of complexes of DNA duplexes and novel benzoxazoles or benzimidazoles by electrospray ionization mass spectrometry

Electrospray ionization mass spectrometry is used to compare the metal ion binding and metal-mediated DNA binding of benzoxazole (1, 2, 3, 4) and benzimidazole (5) compounds and to elucidate the putative binding modes and stoichiometries. The observed metal versus non-metal-mediated DNA binding, as well as the specificity of DNA binding, is correlated with the biological activities of the analogs. The ESI-MS spectra for the antibacterial benzoxazole and benzimidazole analogs 4 and 5 demonstrated non-specific and non-metal-mediated binding to DNA, with the appearance of DNA complexes containing multiple ligands. The anticancer analog 2 demonstrates a clear preference for metal-mediated DNA interactions, with an apparent selectivity for Ni2+ -mediated binding over the more physiologically relevant Mg2+ or Zn2+ cations. Complexation between DNA and the biologically inactive analog 1 was not observed, either in the absence or presence of metal cations.     

Crystal molecular dynamics simulations to speed up MM/PB(GB)SA evaluation of binding free energies of di-mannose deoxy analogs with P51G-m4-Cyanovirin-N

Complexes of two Cyanovirin-N (CVN) mutants, m4-CVN and P51G-m4-CVN, with deoxy di-mannose analogs were employed as models to generate conformational ensembles using explicit water Molecular Dynamics (MD) simulations in solution and in crystal environment. The results were utilized for evaluation of binding free energies with the molecular mechanics Poisson-Boltzmann (or Generalized Born) surface area, MM/PB(GB)SA, methods. The calculations provided the ranking of deoxy di-mannose ligands affinity in agreement with available qualitative experimental evidences. This confirms the importance of the hydrogen-bond network between di-mannose 3'- and 4'-hydroxyl groups and the protein binding site B(M) as a basis of the CVN activity as an effective HIV fusion inhibitor. Comparison of binding free energies averaged over snapshots from the solution and crystal simulations showed high promises in the use of the crystal matrix for acceleration of the conformational ensemble generation, the most time consuming step in MM/PB(GB)SA approach. Correlation between energy values based on solution versus crystal ensembles is 0.95 for both MM/PBSA and MM/GBSA methods.     

In-silico guided discovery of novel CCR9 antagonists

Antagonism of CCR9 is a promising mechanism for treatment of inflammatory bowel disease, including ulcerative colitis and Crohn’s disease. There is limited experimental data on CCR9 and its ligands, complicating efforts to identify new small molecule antagonists. We present here results of a successful virtual screening and rational hit-to-lead campaign that led to the discovery and initial optimization of novel CCR9 antagonists. This work uses a novel data fusion strategy to integrate the output of multiple computational tools, such as 2D similarity search, shape similarity, pharmacophore searching, and molecular docking, as well as the identification and incorporation of privileged chemokine fragments. The application of various ranking strategies, which combined consensus and parallel selection methods to achieve a balance of enrichment and novelty, resulted in 198 virtual screening hits in total, with an overall hit rate of 18%. Several hits were developed into early leads through targeted synthesis and purchase of analogs.     

Balancing the influence of molecular complexity on fingerprint similarity searching

Differences in molecular complexity and size are known to bias the evaluation of fingerprint similarity. For example, complex molecules tend to produce fingerprints with higher bit density than simple ones, which often leads to artificially high similarity values in search calculations. We introduce here a variant of the Tversky coefficient that makes it possible to modulate or eliminate molecular complexity effects when evaluating fingerprint similarity. This has enabled us to study in detail the role of molecular complexity in similarity searching and the relationship between reference and active database compounds. Balancing complexity effects leads to constant distributions of similarity values for reference and database molecules, independent of how compound contributions are weighted. When searching for active compounds with varying complexity, hit rates can be optimized by modulating complexity effects, rather than eliminating them, and adjusting relative compound weights. For reference molecules and active database compounds having different complexity, preferred parameter settings are identified.     

New methods for ligand-based virtual screening: use of data fusion and machine learning to enhance the effectiveness of similarity searching

Similarity searching using a single bioactive reference structure is a well-established technique for accessing chemical structure databases. This paper describes two extensions of the basic approach. First, we discuss the use of group fusion to combine the results of similarity searches when multiple reference structures are available. We demonstrate that this technique is notably more effective than conventional similarity searching in scaffold-hopping searches for structurally diverse sets of active molecules; conversely, the technique will do little to improve the search performance if the actives are structurally homogeneous. Second, we make the assumption that the nearest neighbors resulting from a similarity search, using a single bioactive reference structure, are also active and use this assumption to implement approximate forms of group fusion, substructural analysis, and binary kernel discrimination. This approach, called turbo similarity searching, is notably more effective than conventional similarity searching.     

Icotinib,Almonertinib, and Olmutinib: A 2D Similarity/Docking-Based Study to Predict the Potential Binding Modes and Interactions into EGFR

In the current study, a 2D similarity/docking-based study was used to predict the potential binding modes of icotinib, almonertinib, and olmutinib into EGFR. The similarity search of icotinib, almonertinib, and olmutinib against a database of 154 EGFR ligands revealed the highest similarity scores with erlotinib (0.9333), osimertinib (0.9487), and WZ4003 (0.8421), respectively. In addition, the results of the docking study of the three drugs into EGFR revealed high binding free energies (ΔG_b = −6.32 to −8.42 kcal/mol) compared to the co-crystallized ligands (ΔG_b = −7.03 to −8.07 kcal/mol). Analysis of the top-scoring poses of the three drugs was done to identify their potential binding modes. The distances between Cys797 in EGFR and the Michael acceptor sites in almonertinib and olmutinib were determined. In conclusion, the results could provide insights into the potential binding characteristics of the three drugs into EGFR which could help in the design of new more potent analogs.     

QUASI: a novel method for simultaneous superposition of multiple flexible ligands and virtual screening using partial similarity

The structure of many receptors is unknown, and only information about diverse ligands binding to them is available. A new method is presented for the superposition of such ligands, derivation of putative receptor site models and utilization of the models for screening of compound databases. In order to generate a receptor model, the similarity of all ligands is optimized simultaneously taking into account conformational flexibility and also the possibility that the ligands can bind to different regions of the site and only partially overlap. Ligand similarity is defined with respect to a receptor site model serving as a common reference frame. The receptor model is dynamic and coevolves with the ligand alignment until an optimal self-consistent superposition is achieved. When ligand conformational flexibility is permitted, different superposition models are possible and consistent with the data. Clustering of the superposition solutions is used to obtain diverse models. When the models are used to screen a database of compounds, high enrichments are obtained, comparable to those obtained in docking studies.     

Simulation studies of the fidelity of biomolecular structure ensemble recreation

We examine the ability of Bayesian methods to recreate structural ensembles for partially folded molecules from averaged data. Specifically we test the ability of various algorithms to recreate different transition state ensembles for folding proteins using a multiple replica simulation algorithm using input from "gold standard" reference ensembles that were first generated with a Go-like Hamiltonian having nonpairwise additive terms. A set of low resolution data, which function as the "experimental" phi values, were first constructed from this reference ensemble. The resulting phi values were then treated as one would treat laboratory experimental data and were used as input in the replica reconstruction algorithm. The resulting ensembles of structures obtained by the replica algorithm were compared to the gold standard reference ensemble, from which those "data" were, in fact, obtained. It is found that for a unimodal transition state ensemble with a low barrier, the multiple replica algorithm does recreate the reference ensemble fairly successfully when no experimental error is assumed. The Kolmogorov-Smirnov test as well as principal component analysis show that the overlap of the recovered and reference ensembles is significantly enhanced when multiple replicas are used. Reduction of the multiple replica ensembles by clustering successfully yields subensembles with close similarity to the reference ensembles. On the other hand, for a high barrier transition state with two distinct transition state ensembles, the single replica algorithm only samples a few structures of one of the reference ensemble basins. This is due to the fact that the phi values are intrinsically ensemble averaged quantities. The replica algorithm with multiple copies does sample both reference ensemble basins. In contrast to the single replica case, the multiple replicas are constrained to reproduce the average phi values, but allow fluctuations in phi for each individual copy. These fluctuations facilitate a more faithful sampling of the reference ensemble basins. Finally, we test how robustly the reconstruction algorithm can function by introducing errors in phi comparable in magnitude to those suggested by some authors. In this circumstance we observe that the chances of ensemble recovery with the replica algorithm are poor using a single replica, but are improved when multiple copies are used. A multimodal transition state ensemble, however, turns out to be more sensitive to large errors in phi (if appropriately gauged) and attempts at successful recreation of the reference ensemble with simple replica algorithms can fall short.     

Similar,or dissimilar,that is the question. How different are methods for comparison of compounds similarity?

Comparison of compounds similarity is one of the main strategies of virtual screening protocols. Both similarity and dissimilarity concepts are of great importance during the search for new active compounds. Similarity is important due to the assumption that underlies the process of searching for new drug candidates: structurally similar compounds should induce similar biological response. On the other hand, we are also interested in dissimilarity, as we usually aim to find structurally novel ligands. In the study, we compared several approaches of evaluating compound similarity. Various representations and metrics were applied and we indicated the rate of variation of the results that can occur when shifting from one strategy to another. We compared both general similarity of datasets using different approaches, as well as examined the changes in the set of nearest neighbors when changing one compound representation into another, and the influence of representation/metric settings on the clustering outcome. We hope that the study will be of great help during the preparation of virtual screening experiments, stressing the need for careful selection of the way, the compound similarity is assessed. The differences in the results that can be obtained via the application of particular strategy can significantly influence the outcome of comparison studies; therefore, its settings should be carefully selected beforerunning the comparison.     

Molecular Recognition of DNA by Hoechst Benzimidazoles: Exploring Beyond the Pyrrole‐Imidazole‐Hydroxypyrrole Polyamide‐Pairing Code

A series of three‐ring analogs of the minor‐groove‐binding molecule Hoechst 33258 ( 1 ), consisting of benzimidazole (B), imidazopyridine (P), and hydroxybenzimidazole (H) monomers, have been synthesized in order to investigate both their sequence specificity and binding modes. MPE⋅Fe^II Footprinting has revealed the preference of both PBB and BBB ligands for 5′‐WGWWW‐3′ and 5′‐WCWWW‐3′ tracts, as well as A⋅T‐rich sequences. Affinity‐cleavage titrations show no evidence for a 2 : 1 binding mode of these Hoechst analogs. Importantly, all derivatives are oriented in one direction at each of their binding sites. The implications of these results for the design of minor‐groove‐binding small molecules is discussed.     

Variable selection and model validation of 2D and 3D molecular descriptors>

We have found that molecular shape and electrostatics, in conjunction with 2D structural fingerprints, are important variables in discriminating classes of active and inactive compounds. The subject of this paper is how to explore the selection of these variables and identify their relative importance in quantitative structure-activity relationships (QSAR) analysis. We show the use of these variables in a form of similarity searching with respect to a crystal structure of a known bound ligand. This analysis is then validated through k-fold cross-validation of enrichments via several common classifiers. Additionally, we show an effective methodology using the variables in hypothesis generation; namely, when the crystal structure of a bound ligand is not known.     

Homology model-guided 3D-QSAR studies of HIV-1 integrase inhibitors

In the present study, we report the exploration of binding modes of potent HIV-1 integrase (IN) inhibitors MK-0518 (raltegravir) and GS-9137 (elvitegravir) as well as chalcone and related amide IN inhibitors we recently synthesized and the development of 3D-QSAR models for integrase inhibition. Homology models of DNA-bound HIV-1 IN were constructed on the basis of the X-ray crystal structure of the foamy virus IN-DNA complex (PDB ID: 3L2T ) and used for docking. The binding modes of raltegravir and elvitegravir in our homology models are in accordance with those in the foamy virus structure revealing interactions important for inhibitor-IN binding. To gain further insights into the structural requirements for IN inhibition, three-dimensional quantitative structure activity relationship (3D-QSAR) studies were conducted using raltegravir, elvitegravir, and their analogs; our synthesized 3-keto salicylic acid IN inhibitor series; as well as other structurally related HIV-1 IN inhibitors. In the first part of the study with 103 compounds, atom-fit alignments, I and II, and docking-based alignment, III, were used to develop 3D-QSAR models 1, 2, and 3, respectively, each comprising comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) 3D-QSARs. This initial analysis indicated that the docking-based (structure-based) model 3 performed better than the atom-fit (ligand-based) models 1 and 2, in terms of statistical significance and robustness. Thus, the docking-based alignment was then subsequently used with an expanded data set of 296 compounds for building a more comprehensive 3D-QSAR, model 4. Model 4 afforded good q2 values of 0.70 and 0.75 for CoMFA and CoMSIA 3D-QSARs, respectively, and showed good predictive performance on an external validation test set of 59 compounds with predictive r2 values up to 0.71. The HIV IN-DNA homology model of biological relevance and the comprehensive 3D-QSAR models developed in the present study provide insights and new predictive tools for structure-based design and optimization of IN inhibitors.     

Exploration of 3D activity cliffs on the basis of compound binding modes and comparison of 2D and 3D cliffs

Activity cliffs are formed by pairs or groups of structurally similar compounds having large differences in potency and are focal points of structure-activity relationship (SAR) analysis. The choice of molecular representations is a critically important aspect of activity cliffs analysis. Thus far, activity cliffs have predominantly been defined on the basis of molecular graph or fingerprint representations. Herein we introduce 3D activity cliffs derived from comparisons of experimentally determined compound binding modes. The analysis of 3D activity cliffs is generally applicable to target proteins for which structures of multiple ligand complexes are available. For two popular targets, β-secretase 1 (BACE1) and factor Xa (FXa), public domain X-ray structures with bound inhibitors were collected. Crystallographic binding modes of inhibitors were systematically compared using a 3D similarity method taking conformational, positional, and atomic property differences into account. In addition, standard 2D similarity relationships were also determined. SAR information associated with individual compounds substantially changed when either bioactive conformations or 2D molecular graphs were used for similarity evaluation. 3D activity cliffs were identified for BACE1 and FXa inhibitor sets and systematically compared to 2D cliffs. It was found that less than 40% of 3D activity cliffs were conserved when 2D similarity was applied. The limited conservation of 3D and 2D cliffs provides further evidence for the strong molecule representation dependence of activity cliffs. Moreover, 3D cliffs represent a new class of activity cliffs that convey SAR information in ways that differ from graph-based similarity measures. In cases where sufficient structural information is available, the comparison of 3D and 2D cliffs is expected to aid in SAR analysis and mapping of critical binding determinants.     

Systematic identification and classification of three-dimensional activity cliffs

Activity cliffs were systematically extracted from public domain X-ray structures of targets for which complexes with multiple ligands were available, following the concept of three-dimensional (3D) cliffs. Binding modes of ligands with well-defined potency measurements were compared in a pairwise manner, and their 3D similarity was calculated using a previously reported property density function-based method taking conformational, positional, and chemical differences into account. Requiring the presence of at least 80% 3D similarity and a potency difference of at least 2 orders of magnitude as cliff criteria, a total of 216 well-defined 3D activity cliffs were detected in the Protein Data Bank (PDB). These 3D-cliffs involved a total of 269 ligands active against 38 different targets belonging to 17 protein families. For 255 of these compounds, binding modes were available at high crystallographic resolution. All 3D-cliffs were analyzed in detail and assigned to different categories on the basis of crystallographic interaction patterns. In many instances, differences in ligand-target interactions suggested plausible causes for origins of 3D-cliffs. In other cases, short-range interactions seen in X-ray structures were insufficient to deduce possible reasons for cliff formation. The 3D-cliffs described herein further advance the rationalization of activity cliffs at the level of ligand-target interactions and should also be useful for other applications such as the calibration of energy functions for structure-based design. The pool of identified activity cliffs is provided to enable subsequent structure-based analyses of cliffs.