J. Jesús Naveja
José L. Medina-Franco- 1PECEM, School of Medicine, Universidad Nacional Autónoma de México, Mexico City, Mexico
- 2Department of Pharmacy, School of Chemistry, Universidad Nacional Autónoma de México, Mexico City, Mexico
Herein we introduce the constellation plots as a general approach that merges different and complementary molecular representations to enhance the information contained in a visual representation and analysis of chemical space. The method is based on a combination of a sub-structure based representation and classification of compounds with a “classical” coordinate-based representation of chemical space. A distinctive outcome of the method is that organizing the compounds in analog series leads to the formation of groups of molecules, aka “constellations” in chemical space. The novel approach is general and can be used to rapidly identify, for instance, insightful and “bright” Structure-Activity Relationships (StARs) in chemical space that are easy to interpret. This kind of analysis is expected to be especially useful for lead identification in large datasets of unannotated molecules, such as those obtained through high-throughput screening. We demonstrate the application of the method using two datasets of focused inhibitors designed against DNMTs and AKT1.
Introduction
The concept of chemical space is broadly used in drug discovery because of its multiple potential applications; for instance, in library design, compound or dataset classification, compound selection, exploration of structure-activity relationships (SAR), and navigation though structure-property relationships (SPR) in general. However, a precise unique definition of chemical space is not simple. An even more challenging task is the visual representation of this subjective concept.
Chemical space is usually defined as the set of all possible organic compounds (Lipinski and Hopkins, 2004). It is widely recognized that the virtual chemical space is more than astronomically large, as not even all atoms in the universe would suffice to synthesize a single molecule from each of all the 1063 possible organic compounds of a size up to 30 atoms (Clayden et al., 2012). Nevertheless, massive efforts have been undertaken to enumerate billions of hypothetical organic compounds, thus allowing large virtual screening campaigns to take place (Reymond, 2015; Lyu et al., 2019).
Along with the increasing size of the mapped chemical space, the interest of applying cartographic methods to visualize the space has expanded (Oprea and Gottfries, 2001). As a result, numerous visualization and conceptualization approaches into chemical space have emerged (Larsson et al., 2007; Osolodkin et al., 2015; Naveja and Medina-Franco, 2017). A cornerstone and key aspect of all proposed methods is the molecular representation and parameters used to define the space where the compounds will reside. Chemical space visualizations have to reduce the dimensionality of the problem of comparing molecular structures, which can be done through algorithms such as principal components analysis and t-distributed stochastic neighbor embedding (see Osolodkin et al., 2015).
In most chemical space approaches, it is desirable that chemical analogs are closer to each other than unrelated and dissimilar molecules since this allows machine learning methods to identify clusters of structurally-related molecules (Medina-Franco et al., 2008; Naveja and Medina-Franco, 2015; Naveja et al., 2016, 2018a). In addition, clustering analog series would allow, at least in principle, to map SAR/SPR into that space. However, due to the vast amplitude of the chemical space and the inevitable loss of information with an initially large space projected into lower dimensions, it is expected that non-analog compounds will end up in the same cluster. Also, when many points in the chemical space are considered at once, visualizations become harder to interpret. To address this issue, approaches such as virtual reality have emerged (Probst and Reymond, 2018).
In parallel to such chemical space approaches based on coordinates, scaffold analysis is a more consistent and chemically-intuitive approach for exploring and identifying collections of analogs (Hu et al., 2011). Ever since the pioneering work by Bemis and Murcko (1996), computational identification of chemical scaffolds has been refined. In this line, Stumpfe et al. (2016) recently introduced the analog series-based scaffold (ASBS), a revolutionary scaffold concept that is more flexible and chemically sound than its predecessors. In fact, the ASBS has proven to yield more biologically meaningful structure-activity/property relationships (SA/PR) than other scaffold definitions (Dimova et al., 2016; Kunimoto et al., 2017; Bajorath, 2018; Dimova and Bajorath, 2018).
Although the chemical space of single analog series can be effectively explored and used, for instance, to guide lead optimization programmes (Vogt et al., 2018), methods for analyzing the relationship among scaffolds of different analog series remain to be explored. Of note, a difficulty in this regard emerges as analog-series based scaffolds tend not to be as consistent as Bemis-Murcko scaffolds, since they result from the retrospective analysis of analog series (Bajorath, 2018). Accordingly, a core framework inspired in the design of the ASBS avoids the shortcoming of inconsistency by allowing molecules to be annotated with more than one putative core (Naveja et al., Submitted). Hence, large libraries containing analogs can be condensed into fewer cores. In this way, SA/PR can be preferentially analyzed for the most explored regions of the chemical space: analog series.
Herein, we present a general methodology for applying the putative core framework to produce more concise and meaningful representations of the chemical space. To our understanding, this is the first method designed for charting multiple analog series into a coordinate-based chemical space, thus combining in a single plot two general and useful approaches of molecular representation and mapping. Of note, since within this framework cores may share analogs (i.e., analog series are allowed to share compounds), such cores can be connected, thus resembling constellations in the chemical space. Therefore, we termed the resulting graphics “constellation plots.” As it will be discussed, activity data (or any property of interest) can be mapped into the constellation plot allowing to explore SA/PRs in the space and quickly identify interesting regions in the space. The rest of this methodological paper is organized as follows: first, the concept scheme is presented and the formalism explained through a toy example; thereafter, two case studies using exemplary datasets are presented; finally, we discuss the conclusions and perspectives of this novel approach for combining the scaffold and the chemical space concepts.
Methods
Datasets Used in the Examples
For illustrating the application of constellation plots in two different context of analysis, we used two benchmark datasets that have been previously explored with other analysis approaches. One set was a group of 827 AKT1 inhibitors extracted and curated from ChEMBL (Gaulton et al., 2017; Naveja et al., 2018b). The second dataset was a collection of 286 compounds tested as inhibitors of DNMT (DNA methyltransferases). This second data set was integrated from multiple sources of information as described in Naveja and Medina-Franco (2018). Since this dataset integrates qualitative (such as those containing crystallographic data) and quantitative databases (such as those containing experimental determination of inhibition curves), for this dataset, we use a categorical classification of activity in “active” or “inactive.” The files of the two datasets are included as Supplementary Information.
Chemical Space and Analog Series
As mentioned above, constellation plots fuse two ligand-based concepts: chemical space and core analysis. Standard chemical space analysis is carried out by computing descriptors for a collection of molecules (e.g., physicochemical properties and/or structural features) and then applying dimensionality reduction approaches (Rosén et al., 2009; Osolodkin et al., 2015; González-Medina et al., 2016; Prieto-Martínez et al., 2016; Naveja and Medina-Franco, 2017; Borrel et al., 2018). As a result, every data point represents a single molecule (see Figure 1). This can render many visualizations hard to read and analyze by the naked eye. Furthermore, the numerous descriptors used are combined, such that every axis in the visualization turns out to have a quite abstract meaning. Herein, for the purpose of charting chemical space, t-distributed stochastic neighbor embedding (t-SNE) is used. This methodology reduces the number of data points in the center of the map as compared to other approaches and has been used successfully in chemical space charting (Maaten and Hinton, 2008; Lewis et al., 2015). However, other coordinate-based representations of chemical space can be used in this general approach.
Figure 1. A hypothetical example of a typical chemical space representation based on coordinates. The axes represent the chemical space and have abstract meanings regarding a combination of descriptors. In this case, t-SNE using Morgan fingerprints was applied. Every dot represents a single molecule. In activity landscape modeling, color is used to indicate a property (potency in a particular biological endpoint).
In contrast to chemical space, standard scaffold and analog series analysis aims toward a clear and consistent picture of the relationships among compounds. For instance, a scaffold is a substructure shared by all compounds annotated with it. A state-of-the-art approach for defining analog series-based scaffolds was proposed by Stumpfe et al. (2016). They have reasoned that for a scaffold to be relevant in medicinal chemistry, it should not only be a substructure of a molecule, but it also has to comply with three additional criteria: (i) be a major component of the whole molecule, (ii) be derived from the molecule through retrosynthetic rules, and (iii) summarize an analog series in a particular dataset. A number of computational approaches for obtaining ASBS have been proposed (Dimova et al., 2016; Stumpfe et al., 2016; Bajorath, 2018; Naveja et al., 2019). Within these approaches, an analog series is defined as a subnetwork connected by matched molecular pairs (MMPs) (Griffen et al., 2011).
Chemical space analysis of individual analog series has been carried out to measure progression in lead optimization and saturation of analog series (Kunimoto et al., 2018; Vogt et al., 2018; Yonchev et al., 2018). Nevertheless, the fact that assumption (iii) makes analog series inconsistent in as much as the scaffold definition is dependent on the dataset used (Bajorath, 2018) is a limitation for the exploration of chemical space of multiple analog series at once. In a recent study (Naveja et al., Submitted), we discussed that by removing assumption (iii) two effects take place: first, every molecule is allowed to be annotated to more than a single core (equivalent to the term “scaffold”); and second, complete consistency is achieved as no core annotations are ever omitted for any molecule (see Figure 2). It is within this general core framework that we propose using constellation plots.
Figure 2. Two examples of putative cores computed for two molecules. Note that in this approach the same chemical structure can be its own core (structures at the bottom). After RECAP fragmentation, hydrogens are added to the core structure to avoid invalid valence (marked in red).
Summarizing Analog Series Information in a Dataset Within the General Core Framework
Since the general core framework can assign multiple cores to single molecules, a useful step prior to mapping cores in the chemical space would be summarizing analog series in the smallest number of cores possible. As illustrated in Figure 3, in some instances it is possible to summarize a whole analog series in a single core structure, while in other cases this cannot be done without loss of information. Hence, for avoiding such situations, we did not discard cores unless only one compound mapped to it. Furthermore, if two or more cores mapped to exactly the same compounds, then only the largest core was kept and the others were disregarded from the analysis.
Figure 3. Examples of two analog series with multiple compounds and cores. (A) Analog series that can be summarized in a single core; (B) Analog series formed by multiple cores. In case (B) a single core is not enough for summarizing all information in an analog series.
Constellation Plots
After processing a collection of compounds under the general core framework, information is obtained in multiple regards, namely: (a) the chemical structure of every core; (b) the sets of molecules mapping to each core; (c) the molecules annotated to multiple cores; and (d) the analog series to which each compound and core are annotated. We propose a visualization methodology summarizing these four dimensions in a single graphic: the constellation plot that is schematically illustrated in Figure 4.
Figure 4. Schematic representation of a general form of a constellation plot. Every circle in the plot represents a core; the axes comprise the coordinates of the chemical structure of the core projected into a 2D representation of the chemical space as computed by any of the standard approaches (e.g., generated using continuous properties or molecular fingerprints and applying t-SNE or principal components analysis); the size of the circles indicates the “n” number of compounds annotated to a given core; connected circles are cores sharing compounds; the labels indicate the analog series every point belongs to; the color scale represents the average of a property/activity of the compounds mapping to the core.
Essentially, in a constellation plot, the chemical structure of representative cores in a database (for example, those annotated with a predefined minimum number of compounds) is used to find descriptors and map them into a chemical space as if they were single molecules. The size of the circles is used to represent the relative number of compounds annotated to each core. Cores sharing compounds are connected by lines forming “constellations” in the chemical space. Every circle is labeled with an identifier for the analog series to which each core belongs. Additionally, a color scale can be used to represent an average of a given property or activity of the compounds annotated with each core, thereby turning constellation plots useful for activity landscape modeling (Waddell and Medina-Franco, 2012). Of note, the activity can be, for instance, measured for a single molecular target. However, the property could also be a representative measure of the selectivity or promiscuity profile of all the compounds sharing a core across multiple biological endpoints (see section Conclusions and Perspectives).
Figure 4, as opposed to Figure 1, is able to summarize a larger number of compounds than points depicted and contains information about actual analogs. For instance, analog series I, J, and L form separate clusters, but the cluster top right has multiple chemotypes of distinct analog series. This could not be inferred from clustering algorithms applied to the chemical space information only.
Implementation
All scripts required for producing the data herein reported use free Python code and are made freely available in Supplementary Information. RDkit was used for computing fingerprints and manipulation of chemical structures (http://www.rdkit.org). Scikit-learn was used for computing t-SNE (Pedregosa et al., 2011).
Results and Discussion
The construction of constellations plots and exemplary applications are illustrated with two case studies of general interest in drug discovery. As mentioned in the section Methods, the first example consists of a dataset of 827 AKT1 inhibitors obtained from ChEMBL (Gaulton et al., 2017) and cheminformatically described in Naveja et al. (2018b). The second example employs a data set of 286 DNA methyltransferase (DNMT) inhibitors obtained from the integration of several databases as described in Naveja and Medina-Franco (2018).
Case Study 1: AKT1 Inhibitors
Analogs in this library could be summarized in 144 cores as discussed in the section Methods. The cores were organized in 79 analog series and contained 440 compounds (about half of the initial dataset). Figure 5 is the constellation plot for these data, where it becomes apparent that chemical space and chemical substructure information play simultaneous roles in describing the SAR. For instance, although some inactive cores are close to active cores in chemical space, they are not usually contained in the same analog series. Therefore, these could be categorized as “scaffold cliffs” rather than simple activity cliffs conceptualized as two small molecules with similar structures and very different activities (Maggiora, 2006). In this case, collections of molecules, rather than single molecules, are being compared.
Figure 5. General constellation plot for a dataset of AKT1 inhibitors. It is possible to navigate this map, as observed in Figure 6, where the constellation framed within dashed lines is further explored.
Figure 6 is a zoomed-in picture into a single “bright” (or predominantly active) constellation comprising five analog series and 55 compounds. As it is readily observed, analog series close in the chemical space have only slight dissimilarities within their scaffolds; in this case, they all share a naphthyridine or naphthyridinone scaffold. Constellation plots allow for a more precise visual SAR analysis and generation of hypotheses. For instance, the core associated to analog series 62 has only a different position for the nitrogens in the rings as well as where substitutions occur. Structural studies could then be conducted to elucidate which are the most relevant features for this kind of scaffolds to be active against AKT1. In this regard, a recent publication co-crystallizing 1,6-naphthyridinone derivatives similar to those in analog series 20 has shown that this scaffold is relevant in forming a π-π stacking interaction with the side chain of Trp80 of the PH-domain (Uhlenbrock et al., 2019). Nonetheless, variation of the position of nitrogen atoms in the scaffold were not considered in the cited study. Indeed, previous SAR studies of these analogs have found the position of the nitrogen atoms in these scaffolds to be critical for the activity against AKT (Zhao et al., 2005;Bilodeau et al., 2008).
Figure 6. “Zoom-in” into the constellation plot for AKT1 inhibitors selected from Figure 5. The core in analog series 62 is not as active as the nearby cores in analog series 2 and 12. Few structural differences can be noted for the compounds in this constellation.
Case Study 2: DNMT Inhibitors
Analogs in this library could be summarized in 23 cores following the procedure discussed in the section Methods. The cores were organized in 13 analog series and contained 46 compounds (about 16% of the initial dataset). Compounds in this library have annotated activity with DNMT1, DNMT3A, and/or DNMT3B. Figure 7 shows three constellations plots, where chemical space is the same and colors change to represent the activities against each DNMT. As elaborated on the section Methods, each circle in the plot represents a core in which coordinates in the 2D graph is given by similarity measurements computed from Morgan fingerprints using t-SNE for dimensionality reduction. Labels indicate the analog series to which cores belong. The color represents the percentage of active compounds sharing that core using a continuous color scale from red (less active cores) to yellow (more active cores). For this example of use of the constellation plots, the definition of “active” was determined from integrating qualitative and quantitative data sources as described in Naveja and Medina-Franco (2018). Circles in gray indicate cores with no reported activity for that particular DNMT. The size of the circle indicates the number of compounds sharing the core. Connected circles are cores sharing compounds. Figure 7 also shows the chemical structures of representative cores.
Figure 7. Constellation plots for a dataset of DNA methyltransferase (DNMT) inhibitors tested with DNMT1, DNMT3A, and DNMT3B.
The constellation plots for DNMT inhibitors in Figure 7 allow for rapidly getting several interesting insights of the SAR. For instance, cores at the top left part of the plot from analog series “A” are a bright constellation against DNMT1, i.e., a region in chemical space with active analogs. However, these analogs have not been tested against the other DNMT isoenzymes, which would help determine whether these inhibitors are selective.
Of note, there is a “dark” (or predominantly inactive) constellation in the chemical space of DNMT1 formed by six cores from analog series “D.” This dark constellation, however, is more active overall against DNMT3A and appears to be active against DNMT3B. Furthermore, not all cores in this constellation have been tested against DNMT3A and DNMT3B, where they have greater chances of being active.
The plot also reveals a constellation of nucleoside analogs from series “B” at the bottom-right region of the plot that is, overall, selective toward DNMT3B vs. DNMT1. This series has not been tested against DNMT3A yet. Moreover, most of the cores have been tested in DNMT1 only, thus hindering discussions on selectivity. In this regard, analysis of constellation plots is visually helpful in guiding multitarget drug discovery campaigns and in finding opportunities for selectivity.
Conclusions and Perspectives
We introduced a novel approach for combining chemical space and analog series methodologies into a single descriptive analysis that can be summarized in a constellation plot. Adding the analog series concept into the chemical space facilitates discussions of regions in the space, as every point summarizes a collection of analogs. A so-called “constellation in chemical space” can be conceptualized as those regions in chemical space formed by core scaffolds with similar structure (as defined by a coordinate-based projection). Mapping activity on the plot readily uncovers active and inactive zones, e.g., bright or dark regions, in chemical space. Of note, constellation plots would be useful for exploring virtually any chemical property, such as biological activity (as demonstrated with two case studies), but also physicochemical properties, complexity or selectivity statistics. In this regard, constellation plots are a flexible approach with multiple potential applications in academia and industry, aiding in the quest of finding potential leads from large collections of screening data (e.g., such as that produced by high-throughput screening campaigns). One of the next steps of this work is the application of the constellations plots to navigate through cell selectivity data of a comprehensive screening dataset. Results will be disclosed in due course.
Data Availability
The datasets generated for this study can be found in Supplementary Information.
Author Contributions
JN participated in the conceptualization, data gathering, data analysis, and drafted the first version of the manuscript. JM-F participated in the conceptualization, data analysis, and revision of the manuscript.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
JN was thankful to CONACyT for the granted scholarship number 622969. We also thank Consejo Nacional de Ciencia y Tecnología (CONACyT) for grant 282785.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fchem.2019.00510/full#supplementary-material
References
Bajorath, J. (2018). Improving the utility of molecular scaffolds for medicinal and computational chemistry. Future Med. Chem. 10, 1645–1648. doi: 10.4155/fmc-2018-0106
Bemis, G. W., and Murcko, M. A. (1996). The properties of known drugs. 1. Molecular frameworks. J. Med. Chem. 39, 2887–2893. doi: 10.1021/jm9602928
Bilodeau, M. T., Balitza, A. E., Hoffman, J. M., Manley, P. J., Barnett, S. F., Defeo-Jones, D., et al. (2008). Allosteric inhibitors of Akt1 and Akt2: a naphthyridinone with efficacy in an A2780 tumor xenograft model. Bioorg. Med. Chem. Lett. 18, 3178–3182. doi: 10.1016/j.bmcl.2008.04.074
Borrel, A., Kleinstreuer, N. C., and Fourches, D. (2018). Exploring drug space with ChemMaps.com Bioinformatics 34, 3773–3775. doi: 10.1093/bioinformatics/bty412
Clayden, J., Greeves, N., and Warren, S. (2012). Organic Chemistry. Oxford, UK: Oxford University Press. Available online at: http://global.oup.com/ukhe/product/organic-chemistry-9780199270293?cc=mx&lang=en&
Dimova, D., and Bajorath, J. (2018). Collection of analog series-based scaffolds from public compound sources. Future Sci. OA 4:FSO287. doi: 10.4155/fsoa-2017-0135
Dimova, D., Stumpfe, D., Hu, Y., and Bajorath, J. (2016). Analog series-based scaffolds: computational design and exploration of a new type of molecular scaffolds for medicinal chemistry. Future Sci. OA 2:FSO149. doi: 10.4155/fsoa-2016-0058
Gaulton, A., Hersey, A., Nowotka, M., Bento, A. P., Chambers, J., Mendez, D., et al. (2017). The ChEMBL database in 2017. Nucleic Acids Res. 45, D945–D954. doi: 10.1093/nar/gkw1074
González-Medina, M., Prieto-Martínez, F. D., Naveja, J. J., Méndez-Lucio, O., El-Elimat, T., Pearce, C. J., et al. (2016). Chemoinformatic expedition of the chemical space of fungal products. Future Med. Chem. 8, 1399–1412. doi: 10.4155/fmc-2016-0079
Griffen, E., Leach, A. G., Robb, G. R., and Warner, D. J. (2011). Matched molecular pairs as a medicinal chemistry tool. J. Med. Chem. 54, 7739–7750. doi: 10.1021/jm200452d
Hu, Y., Stumpfe, D., and Bajorath, J. (2011). Lessons learned from molecular scaffold analysis. J. Chem. Inf. Model. 51, 1742–1753. doi: 10.1021/ci200179y
Kunimoto, R., Dimova, D., and Bajorath, J. (2017). Application of a new scaffold concept for computational target deconvolution of chemical cancer cell line screens. ACS Omega 2, 1463–1468. doi: 10.1021/acsomega.7b00215
Kunimoto, R., Miyao, T., and Bajorath, J. (2018). Computational method for estimating progression saturation of analog series. RSC Adv. 8, 5484–5492. doi: 10.1039/C7RA13748F
Larsson, J., Gottfries, J., Muresan, S., and Backlund, A. (2007). ChemGPS-NP: tuned for navigation in biologically relevant chemical space. J. Nat. Prod. 70, 789–794. doi: 10.1021/np070002y
Lewis, R., Guha, R., Korcsmaros, T., and Bender, A. (2015). Synergy maps: exploring compound combinations using network-based visualization. J. Cheminform. 7, 36. doi: 10.1186/s13321-015-0090-6
Lipinski, C., and Hopkins, A. (2004). Navigating chemical space for biology and medicine. Nature 432, 855–861. doi: 10.1038/nature03193
Lyu, J., Wang, S., Balius, T. E., Singh, I., Levit, A., Moroz, Y. S., et al. (2019). Ultra-large library docking for discovering new chemotypes. Nature 566, 224–229. doi: 10.1038/s41586-019-0917-9
Maaten, L., and Hinton, G. (2008). Visualizing data using t-SNE. J. Mach. Learn. Res. 9, 2579–2605. Available online at: http://www.jmlr.org/papers/volume9/vandermaaten08a/vandermaaten08a.pdf
Maggiora, G. M. (2006). On outliers and activity cliffs–why QSAR often disappoints. J. Chem. Inf. Model. 46, 1535. doi: 10.1021/ci060117s
Medina-Franco, J., Martinez-Mayorga, K., Giulianotti, M., Houghten, R., and Pinilla, C. (2008). Visualization of the chemical space in drug discovery. CAD 4, 322–333. doi: 10.2174/157340908786786010
Naveja, J. J., Cortés-Benítez, F., Bratoeff, E., and Medina-Franco, J. L. (2016). Activity landscape analysis of novel 5α-reductase inhibitors. Mol. Divers. 20, 771–780. doi: 10.1007/s11030-016-9659-x
Naveja, J. J., and Medina-Franco, J. L. (2015). Activity landscape sweeping: insights into the mechanism of inhibition and optimization of DNMT1 inhibitors. RSC Adv. 5, 63882–63895. doi: 10.1039/C5RA12339A
Naveja, J. J., and Medina-Franco, J. L. (2017). ChemMaps: Towards an approach for visualizing the chemical space based on adaptive satellite compounds. F1000Res. 6:Chem Inf Sci-1134. doi: 10.12688/f1000research.12095.2
Naveja, J. J., and Medina-Franco, J. L. (2018). Insights from pharmacological similarity of epigenetic targets in epipolypharmacology. Drug Discov. Today 23, 141–150. doi: 10.1016/j.drudis.2017.10.006
Naveja, J. J., Norinder, U., Mucs, D., López-López, E., and Medina-Franco, J. L. (2018a). Chemical space, diversity and activity landscape analysis of estrogen receptor binders. RSC Adv. 8, 38229–38237. doi: 10.1039/C8RA07604A
Naveja, J. J., Oviedo-Osornio, C. I., Trujillo-Minero, N. N., and Medina-Franco, J. L. (2018b). Chemoinformatics: a perspective from an academic setting in Latin America. Mol. Divers. 22, 247–258. doi: 10.1007/s11030-017-9802-3
Naveja, J. J., Vogt, M., Stumpfe, D., Medina-Franco, J. L., and Bajorath, J. (2019). Systematic extraction of analogue series from large compound collections using a new computational compound–core relationship method. ACS Omega 4, 1027–1032. doi: 10.1021/acsomega.8b03390
Oprea, T. I., and Gottfries, J. (2001). Chemography: the art of navigating in chemical space. J. Comb. Chem. 3, 157–166. doi: 10.1021/cc0000388
Osolodkin, D. I., Radchenko, E. V., Orlov, A. A., Voronkov, A. E., Palyulin, V. A., and Zefirov, N. S. (2015). Progress in visual representations of chemical space. Expert Opin. Drug Discov. 10, 959–973. doi: 10.1517/17460441.2015.1060216
Pedregosa, F., Varoquaux, G., Gramfort, A., Michel, V., Thirion, B., Grisel, O., et al. (2011). Scikit-learn: machine learning in Python. J. Mach. Learn. Res. 12, 2825–2830. Available online at: http://www.jmlr.org/papers/volume12/pedregosa11a/pedregosa11a.pdf
Prieto-Martínez, F. D., Gortari, E. F., Méndez-Lucio, O., and Medina-Franco, J. L. (2016). A chemical space odyssey of inhibitors of histone deacetylases and bromodomains. RSC Adv. 6, 56225–56239. doi: 10.1039/C6RA07224K
Probst, D., and Reymond, J.-L. (2018). Exploring drugbank in virtual reality chemical space. J. Chem. Inf. Model. 58, 1731–1735. doi: 10.1021/acs.jcim.8b00402
Reymond, J.-L. (2015). The chemical space project. Acc. Chem. Res. 48, 722–730. doi: 10.1021/ar500432k
Rosén, J., Lövgren, A., Kogej, T., Muresan, S., Gottfries, J., and Backlund, A. (2009). ChemGPS-NP(Web): chemical space navigation online. J. Comput. Aided Mol. Des. 23, 253–259. doi: 10.1007/s10822-008-9255-y
Stumpfe, D., Dimova, D., and Bajorath, J. (2016). Computational method for the systematic identification of analog series and key compounds representing series and their biological activity profiles. J. Med. Chem. 59, 7667–7676. doi: 10.1021/acs.jmedchem.6b00906
Uhlenbrock, N., Smith, S., Weisner, J., Landel, I., Lindemann, M., Le, T. A., et al. (2019). Structural and chemical insights into the covalent-allosteric inhibition of the protein kinase Akt. Chem. Sci. 10, 3573–3585. doi: 10.1039/c8sc05212c
Vogt, M., Yonchev, D., and Bajorath, J. (2018). Computational method to evaluate progress in lead optimization. J. Med. Chem. 61, 10895–10900. doi: 10.1021/acs.jmedchem.8b01626
Waddell, J., and Medina-Franco, J. L. (2012). Bioactivity landscape modeling: chemoinformatic characterization of structure-activity relationships of compounds tested across multiple targets. Bioorg. Med. Chem. 20, 5443–5452. doi: 10.1016/j.bmc.2011.11.051
Yonchev, D., Vogt, M., Stumpfe, D., Kunimoto, R., Miyao, T., and Bajorath, J. (2018). Computational assessment of chemical saturation of analogue series under varying conditions. ACS Omega 3, 15799–15808. doi: 10.1021/acsomega.8b02087
Keywords: analog series, data visualization, descriptor, scaffold, structure-property relationships
Citation: Naveja JJ and Medina-Franco JL (2019) Finding Constellations in Chemical Space Through Core Analysis. Front. Chem. 7:510. doi: 10.3389/fchem.2019.00510
Received: 14 May 2019; Accepted: 03 July 2019;
Published: 16 July 2019.
Edited by:
Simona Rapposelli, University of Pisa, ItalyCopyright © 2019 Naveja and Medina-Franco. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: J. Jesús Naveja, naveja@comunidad.unam.mx; José L. Medina-Franco, medinajl@unam.mx