Next Article in Journal
Polyhydric Stigmastane-Type Steroids Derivative from Vernonia amygdalina and Their Anti-Neuroinflammatory Activity
Next Article in Special Issue
Structure-Guide Design and Optimization of Potential Druglikeness Inhibitors for TGFβRI with the Pyrrolopyrimidine Scaffold
Previous Article in Journal
Stability Study, Quantification Method and Pharmacokinetics Investigation of a Coumarin–Monoterpene Conjugate Possessing Antiviral Properties against Respiratory Syncytial Virus
Previous Article in Special Issue
Chelation of Zinc with Biogenic Amino Acids: Description of Properties Using Balaban Index, Assessment of Biological Activity on Spirostomum Ambiguum Cellular Biosensor, Influence on Biofilms and Direct Antibacterial Action
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Deconstructing Markush: Improving the R&D Efficiency Using Library Selection in Early Drug Discovery

by
Leticia Manen-Freixa
,
José I. Borrell
,
Jordi Teixidó
and
Roger Estrada-Tejedor
*
Molecular Design Laboratory, Grup de Química Farmacèutica, IQS School of Engineering, Universitat Ramon Llull, Via Augusta 390, E-08017 Barcelona, Spain
*
Author to whom correspondence should be addressed.
Pharmaceuticals 2022, 15(9), 1159; https://doi.org/10.3390/ph15091159
Submission received: 12 July 2022 / Revised: 2 September 2022 / Accepted: 14 September 2022 / Published: 18 September 2022
(This article belongs to the Special Issue Structure and Ligand Based Drug Design)

Abstract

:
Most of the product patents claim a large number of compounds based on a Markush structure. However, the identification and optimization of new principal active ingredients is frequently driven by a simple Free Wilson approach, leading to a highly focused study only involving the chemical space nearby a hit compound. This fact raises the question: do the tested compounds described in patents really reflect the full molecular diversity described in the Markush structure? In this study, we contrast the performance of rational selection to conventional approaches in seven real-case patents, assessing their ability to describe the patent’s chemical space. Results demonstrate that the integration of computer-aided library selection methods in the early stages of the drug discovery process would boost the identification of new potential hits across the chemical space.

1. Introduction

The pharmaceutical industry has always been challenged in improving the research and development (R&D) efficiency [1,2]. Following the fail early strategy, having a wide range of molecular structures in the initial stages is crucial to increase the chances of finding a compound with biological activity against the target under study. In this sense, the exploration of molecular diversity during the hit or lead discovery phases plays a pivotal role. The description and exploration of the drug-like chemical space (the so-called drugspace [3]) has always been of great concern given the overwhelming number of molecules that can be obtained by fragment combination (which was barely estimated to be 1060 small molecules) [4]. In light of this scenario, finding new molecules with optimal drug-like properties becomes as difficult as finding a needle in a haystack, at least without the support of computer-aided techniques.
Many cheminformatic tools have been developed for better mapping of the chemical space by defining each molecule in a chemical library as a point in the molecular descriptor’s space (or using principal component analysis): distance-based approaches [5,6,7,8], cluster-based selections [9], cell-based selections [10,11] or optimization-based selections [12,13]. These techniques allow the rational navigation through the chemical space (a process known as chemography [3,14,15,16,17]), contributing to hit discovery [18,19,20,21] and the identification of unexplored regions that might hide biologically active compounds. Previously reported results in our group [22,23] proved that rational selection methods did improve the hit identification step when studying a library containing 125,396 analogs of HEPT, an inhibitor of HIV-1 reverse transcriptase. The goodness of 11 available diversity selection methods was assessed and the final selection of 25 compounds that covered 90% of the chemical space showed a broad activity range (with EC50 between 0.05 μM and >90 μM), even better than the reference compound (HEPT, EC50 = 3.3 μM).
Markush structures stated in the claims part of a product patent delimit the chemical space covered by a patent. With the appearance of a combinatorial library, they describe the analogs of a certain drug by defining a molecular skeleton that contains one or more variable substructures listed and represented by particular or generic fragments [24,25]. It is worth noting that the description of a Markush structure does not imply the proven synthesis of all the compounds derived from the combinatorial library. In fact, just a small representation of it is commonly reported, and most of the compounds found in the literature are highly similar to the original hit, evincing the application of a Free Wilson [26,27,28,29] approach during the discovery process (which assumes an additive and independent contribution of specific substituent groups on the biological activity of a molecule). This methodology, widely used in drug discovery, leads to the exploration of a tiny space surrounding the initial bioactive molecule, leaving a great part of the chemical space unexplored, which may hide potential lead candidates that might surpass the activity or undesired side effects of the original hit. In this sense, the integration of rational selection algorithms in the early drug discovery pipeline could expand the exploration of the chemical space, allowing the selection of more distant representative compounds [10,30,31,32,33,34,35] or guiding repurposing campaigns [36]. This leads one to ask, what would have happened if the authors of a given patent had applied rational selection algorithms to determine which compounds should be synthesized and biologically tested? How much better could they have explored the chemical space? To answer these questions, we have applied diversity selection methods to combinatorial libraries of different sizes, derived from the Markush structure described on the patent of seven new chemical entities approved by the FDA. We have compared the representativeness between the compounds claimed in each patent with reported activity found in the literature and subsets of compounds rationally selected using cluster analysis and partition-based (or partitioning) methods.

2. Results and Discussion

2.1. Markush Combinatorial Library

To demonstrate that rational selection can improve the efficiency to describe a chemical space, we have studied seven drug patents approved by the FDA (six of them in the last 12 years) referring to different therapeutic targets, with different sizes (in terms of the number of analogs included in the Markush structure) and including on-patent (ONP) and off-patent (OFP) drugs. OFP drugs have been included in the study to discuss whether the exploration of the chemical space is limited due to the situation of their property rights. Accordingly, Leflunomide (the oldest approved drug considered) has been added to the study to verify that this exploratory deficiency is still present in an older patent (see Table 1).
Thus, we have tested the goodness of rational selection over the traditional approach using seven different datasets obtained from the Markush structure stated in the patents of Dacomitinib, Abemaciclib, Tafenoquine, Ertugliflozin, Rufinamide, Azilsartan Medoxomil and Leflunomide drugs [37,38,39,40,41,42,43,44], retrieved from the Orange Book [45]. Corresponding Markush combinatorial libraries (MCL) have been developed from the Markush structure stated in the patent of each drug (Figure 1). The complete list of substituents for each dataset can be found in the Supporting Information (Tables S1–S7).
As expected, some of the Markush structures use vague indeterminations tending to protect vast chemical spaces. As an example, in the present cases of study, the enumeration of Dacomitinib analogs would lead to a large library of more than 25 million compounds, hampering the computational manageability of the dataset. Therefore, in this unique case, the MCL database was reduced to the combination of its explicitly defined structures in the patent (see Section 3).

2.2. Bibliographical Database and Bibliographic Combinatorial Library

Bibliographical databases (BD) refer to all the compounds derived from the Markush structure that have been reported as synthesized and, in some cases, biologically tested. To describe all the explored space known until the date for each patent, we have included not only the individual molecules explicitly declared on the patent’s claims but also all ulterior derivatives found in PubChem [46] (even if their declared current application differs from the one described in the original patent). It is worth mentioning that in this study, we have only used publicly available compounds. As an example, the set of 58 Tafenoquine analogs found in PubChem includes the seven molecules claimed in patents [38,42] and 51 molecules reported in the literature (belonging to the same Markush structure) for different applications (27 in the context of malaria [47,48,49,50,51,52], 5 for parasitic diseases [53], 13 as inhibitors of monoamine oxidase A [54] and 6 with undefined biological activity).
The number of studied compounds in BD sets differs excessively from the size of the MCL library. On the one hand, the low number of derivatives would evince the existence of synthetical limitations that hamper the obtention of more compounds. On the other hand, it would confirm the application of a highly focused exploratory methodology. For this reason, we have studied a second combinatorial database (named Bibliographic Combinatorial Library, BCL) for each patent by combining only the substituents present in BD, aiming to represent the real combinatorial space synthetically accessible until the date.

2.3. Clustering Methods for Chemical Space Exploration

Noting that the performance of clustering methods depends on dataset distribution [55,56], eight clustering and two partitioning methods were applied to assess the chemical space described in each patent of study.
Hierarchical agglomerative clustering (HRC) with a bottom-up approach with six linkage methods (single, complete, median, average, centroid and Ward), K-Means (KMN) and K-Medoids (KMED), binning and optimum variance binning (OV binning) were firstly used to divide the chemical spaces into k-optimal clusters in order to see the grade of homogeneity in the population distribution per cluster. The k-optimal value was firstly calculated through the average silhouette criterion [57,58,59], using KMN as the standard clustering methodology for all cases. However, no optimal cluster was obtained for all cases as the average silhouette score measured showed a concave curve trend with large maximum values (see Figure S4). This may be caused by the high correlation of the data due to their combinatorial nature and the large size of the datasets, which commonly suffer from the curse of dimensionality and a lack of efficiency. Given the first inconclusive attempts, new standard values commonly used as a rule of thumb in clustering analysis were considered ( N 2 , N and 10 · N ) to find the most representative partitioning size to be applied in further comparative discussions.
The case of Tafenoquine’s chemical space (Table 2) serves as an explanatory example to discuss the population distribution using each different size. The aim of our approach was to find the most balanced distribution of data, avoiding as many singletons and hyper-populated clusters as possible, unveiled by high standard deviation values. Results evidenced that the N clusters seemed to be a good compromise. This size first shows more balanced population frequencies per cluster than N 2 , whose results show large values in their standard deviations. Moreover, compared to 10 · N , N clustering presented a much lesser number of singletons.
After identifying N as the best clustering size for our study, the population distribution was assessed for all the clustering methods. HRC single, median and centroid clustering methods were discarded as they tend to represent a unique or a few hyper-populated clusters and a high number of singletons, losing the homogeneous representativeness of the chemical space (Figure S5 in Supporting Information). This distribution is the common result of many hierarchical agglomerative clustering steps with a bottom-up approach. Therefore, HRC average, HRC complete and Ward have been considered as hierarchical agglomerative clustering methodologies, KMN and KMED as non-hierarchical relocation clustering algorithms and OV binning as a representation of cell-based partitioning methods.

2.4. Bibliographical Representativeness in Its Chemical Space

To assess the degree of representativeness of BD in the chemical space claimed in a patent (MCL), space and population coverage (SC and PC, respectively) values have been calculated when dividing the MCL space into N B D partitions. Results show that a selection of an equal-sized set using random sampling better represents the chemical space than BD compounds (Table 3 and Table 4).
For example, the 58 molecules of Tafenoquine found in the literature are placed in 10 HRC Ward clusters, leading to 17.2% SC (10 out 58 clusters), and they include 5339 analogs, resulting in 18.8% PC (4789 out 25,472 compounds). In these conditions, even a random selection is able to achieve 61.8% SC and 65.3% PC. This trend is observed regardless of the clustering or partition-based method used and in all the seven cases. Only the application of the HRC Complete cluster algorithm on Azilsartan Medoxomil leads to better results than averaged random selections—although, in this case, results are clearly determined by the low number of compounds to select (NBD = 4). These results evince that the chemical space claimed in a drug patent is poorly described.
Needless to say, the use of rational selection would always cover 100% of the current partitioned chemical space by selecting one molecule of each cluster. Then, following the abovementioned example, a rational selection of 58 molecules would represent the 100% SC and PC of Tafenoquine’s chemical space. However, this does not imply the synthetical feasibility of the chosen compounds.

2.5. Comparing the Chemical Space Described by MCL, BD, and BCL

Quantifying the reduction of the MCL space that implies the use of BCL might be of key importance since its compounds are the ones genuinely expected to be synthetically feasible. This confirmation relies on the fact that BCL compounds include only the fragments coming from real studied candidates, constituting a more legitimate representation of the drug patent. For this purpose, the coverage of the MCL, BD and BCL libraries was compared by assessing their distribution along with each principal component and their density plots when projected on the space of the first three principal components (Figure 2). As expected, BD and BCL progressively improve the description of MCL, although neither of them, except for Dacomitinib and Leflunomide’s analogs (later discussed), significantly cover the chemical space derived from the Markush structure.
To quantify this observation, cluster analysis was performed, setting the number of clusters to a standard size of N . Firstly, according to the results, aligned with the previous study, random selection once again was shown to better represent the overall space rather than compounds derived from the current R&D methodology (BD), even with a higher number of clusters. For the sake of clarity, only results for Tafenoquine and Dacomitinib, using HRC average clustering and OV binning, are shown (Table 5 and Table 6). The results of all the libraries can be found in the Supporting Information (Tables S8–S14).
Although having changed the number of partitions, results agree with the ones described in the previous section: both SC and PC percentages obtained using a random selection are better than those obtained with the bibliographic database, either BD or BCL. Again, bibliographical databases only show better results in the HRC average method applied to Azilsartan Medoxomil, due to its size (Table S13). Overall, contrasting the obtained results with N and NBD partitioned space, one may estimate an optimal and synthetically manageable number of samples that could be chosen to be synthesized.
Although the BCL library includes a higher number of compounds, in many cases, it is not able to represent the MCL space; hence, this result would serve as evidence that the synthetically accessible chemical space for each dataset is still poorly known as the rational R&D methodology has followed a mainly focused trend around some original hits. In fact, the use of a rational selection of BD compounds would increase the efficiency of both traditional and cherry-picking methodologies in terms of coverage. As an example, when applying a rational selection of 58 Tafenoquine analogs, an optimal value of 36.3% of SC could be achieved (58 out of 160 clusters).
In contrast to other libraries, Dacomitinib’s BD showed high SC and PC values and a very centered distribution. This is explained by the reduction of the substituents considered in library enumeration (only those in claim 5, see Materials and Methods). This affected not only the MCL size (obtaining 16,530 compounds) but also the BCL (798 compounds). Consequently, the BCL comprises a significant number of analogs of the used dataset of 16,530 compounds. Thus, much lower SC and PC values would be expected when considering the full combinatorial database with more than 25 million analogs.
Additionally, Leflunomide is the unique example in which the number of molecules found BD (NBD) was greater than N , so the size of the corresponding BCL has better representativeness in its chemical space, rather than a random choice of N .

2.6. Towards a More Efficient Methodology

In light of the results of our study, enough evidence has been exposed to prove the lack of chemical space exploration in the traditional R&D methodology. This is the result of a procedure that relies on the hit-to-lead optimization step, which commonly aims to find an optimal compound around the original hit, typically involving a Free Wilson approach during the process (Figure 3). The synthesis and biological evaluation of these analogs is performed, leading to a lead or directly a drug candidate and, accordingly, a Markush structure is settled, defending that its involved analogs may present the same biological behavior as the original hit. In fact, even the fragment combination of the reported structures (BCL) does not represent properly the drug’s chemical space derived from the Markush structure (MCL).
This fact can be observed in the aforementioned example, considering the division of Tafenoquine’s chemical space in 58 (NBD) HRC Ward clusters. BD compounds are located in focused regions of the chemical space, while the rational selection is spread throughout it (Figure 4), leading to a better representation of the whole space.
Hence, we defend an alternative or complementary approach that departs from the combinatorial library obtained from a theoretical Markush structure or from the fragment combination of and original scaffold explored in previous studies (which would ensure the synthetical feasibility of its analogs). Secondly, a computational study, involving space clustering or partitioning, is suggested to rationally choose a handleable number of compounds to synthesize and test that may unveil a better lead in unexplored regions. Thus, this methodology would better consider the chemical diversity of the original Markush set, being the rational selected compounds a significant representation for the issue of study or further repurposing approaches.

3. Materials and Methods

3.1. Enumeration of Combinatorial Libraries

Markush combinatorial libraries (MCL) were fully enumerated for all patents except for Dacomitinib, which was reduced to a computationally handleable size by only combining the fragments present in the 51 compounds described in claim 5 (including all the positional substitutions in the aryl ring, reducing the more than 25 million original compounds to 19 × 15 × 58 = 16,530). On the contrary, in the antimalarial Tafenoquine case, the library of analogs was extended to the combination of the structures found in several expired patents that included this drug. The BCL databases were further prepared following the protocol described.
MarvinSketch was used for drawing, displaying and enumerating the Markush analogs of the chemical libraries (Marvin 20.21, 2020, ChemAxon Ltd.) [60].
The PubChem search has been performed programmatically by using the PUG-REST [61] Application Programming Interface (API). The whole procedure was integrated into a Python script, which checked the presence of a given compound (entered as a SMILES) on PubChem.

3.2. Describing the Chemical Space

All databases were desalted, protonated at pH 7, and their partial charges were calculated using the MMFF94x forcefield in MOE2020.09 software [62]. A total of 206 1D and 2D molecular descriptors were calculated using MOE2020.09 (the complete list of molecular descriptors is available in the Supporting Information, Table S15) to describe the chemical space mathematically. For all combinatorial libraries derived from each Markush structure (MCL libraries), the space dimensionality was reduced using Principal Component Analysis (PCA), keeping 95% of the original variance. The PCA was performed through Python scripting using Scikit-Learn [63].

3.3. Clustering and Partitioning Methodologies

The selection of diverse subsets was performed by comparing cluster analysis and partition-based methods on the PCA chemical space. Hierarchical agglomerative clustering with a bottom-up approach (HRC) was performed using the Scikit-Learn toolkit [63], considering Euclidean metrics and assessing different linkage criteria (single, complete, median, average, centroid and Ward). The K-means (KMN) and K-medioids (KMED) non-hierarchical clustering relocation algorithms were performed using the Scikit-Learn [63] and pyclust toolkits, respectively (pyclust PyPI. available online: https://pypi.python.org/pypi/pyclust/0.2.0, accessed on 10 July 2021). To assess the goodness of partition-based methods, binning and optimum variance binning (OV binning) were implemented in Python, adapting the procedure reported by Pascual et al. [22] (see Supporting Information, Figures S1–S3).

3.4. Space and Population Coverage

The molecular representativeness and space description of selected diverse subsets were assessed using two standard parameters: the space coverage and the population coverage [64,65].
On the one hand, space coverage (SC) represents the percentage of selected (occupied) clusters or bins (koc) by a given number of selected compounds over the total number of partitions, K (Equation (1)).
k o c K · 100
On the other hand, given a database of N compounds, population coverage (PC) measures the SC weighted by the occupancy in each cluster or bin, by dividing the population of the occupied clusters (noc) among N (Equation (2)).
n o c N · 100
The procedure to calculate SC and PC is exemplified in Figure 5.

3.5. Number of Clusters

Different numbers of clusters have been studied in the present work to assess the chemical space in terms of population and space coverage. Actually, the number of partitions (whether they are clusters or bins) ultimately determines the number of compounds to pick via rational selection, since they will correspond to the representative compound from each partition.
Preliminary studies were carried out to assess the behavior of a k-optimal, N 2 , N and 10 · N space fragmentation, being N the total number of compounds in each dataset. Finally, N was taken as a reference in clustering to optimally represent large datasets. Moreover, BD and BCL sizes (NBD and NBCL, correspondingly) were also used to discuss the representativeness of the analogs known until the date, in contrast with a random selection of the same size. Random selections were calculated as the mean value of the coverages represented by NBD and NBCL random samples (with 5000 repetitions) as a contrasting result with data found in the bibliography or other rational selections.

4. Conclusions

The hit-to-lead process in drug discovery has been traditionally based on the application of the Free Wilson approach, according to which, after hit identification, the structure of the drug candidate is progressively modified, attempting to improve its biological activity. Hence, the resulting procedure allows for exploring the surrounding chemical space of the initial hit compound, but there could be regions that remain unexplored, compromising the R&D efficiency.
We have assessed how well the chemical space claimed in a patent is actually explored, using seven patents as examples. For all cases, the space explored in the literature (BD) for each combinatorial library is very small, approximately 20% on average when clustering the chemical space in N clusters. Moreover, results show that even a random selection (by cherry-picking) may lead to better coverage than the molecules reported in the literature. These results are in agreement with results previously reported by our research group regarding the study of the chemical space described by HEPT analogs [22,23].
It has also been evidenced that, in most cases (5 out of 7), even the synthetically accessible combinatorial library (BCL), resulting from the fragmental combination of the molecules described in the literature, is not entirely representative of the chemical space (with lower values than a random study of N molecules).
Neither the real space explored (BD) nor the fragment combination of the studied molecules (BCL) significantly represents the space defined by the combinatorial libraries derived from the Markush structure, especially when they are compared with the coverage obtained by a statistically random sampling of N molecules. Thus, there is a large part of the chemical space claimed on patents that remains unexplored, and it can hide potential leads that may surpass the activity of the original hit or reduce undesired side effects. Rational selection algorithms could assist the traditional methodology to optimize the selection of representative compounds of a given chemical space. This could be applied to explore many pharmacokinetic profiles, such as toxicity, biological activity or solubility, among others.
Results reinforce the proposal to integrate the rational selection in the R&D process in early drug discovery or combine it in a mixed methodology involving a local optimization around the original hit. Nevertheless, it should be noted that many large libraries derived from a Markush structure may present candidates with problematic or even unfeasible synthesis and, hence, proper data curation is mandatory before proceeding to a definitive rational selection.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ph15091159/s1, Figure S1: Description of the Binning and OV Binning partitioning methods implemented in this study; Figure S2: Application of binning and OV binning on Tafenoquine’s dataset to divide the chemical space into 8 bins, Figure S3: OV Binning partitioning representation of Tafenoquine’s case study dataset for 2, 4, 8, 16, 32 and 64 bins; Figure S4: Silhouette Score Elbow plot for KMN clustering of Tafenoquine database; Figure S5: Boxplots of population distribution for the 6 libraries of analogues and √N clusters; Table S1: Dacomitinib fragment library (16,530 analogues); Table S2: Abemaciclib fragment library (45,696 analogues); Table S3: Tafenoquine fragment library (25,472 analogues); Table S4: Ertugliflozin fragment library (14,194 analogues); Table S5: Rufinamide fragment library (8959 analogues); Table S6: Azilsartan Medoxomil fragment library (1110 analogues); Table S7: Leflunomide fragment library (5641 analogues); Table S8: SC and PC for Dacomitinib analogues ( k = N ); Table S9: SC and PC for Abemaciclib analogues ( k = N ); Table S10: SC and PC for Tafenoquine analogues ( k = N ); Table S11: SC and PC for Ertugliflozin analogues ( k = N ); Table S12: SC and PC for Rufinamide analogues ( k = N ); Table S13: SC and PC for Azilsartan Medoxomil analogues ( k = N ); Table S14: SC and PC for Leflunomide analogues ( k = N ); Table S15: List of 206 1D and 2D descriptors calculated using MOE2020.09.

Author Contributions

Conceptualization, R.E.-T. and J.T.; methodology, R.E.-T. and J.T.; investigation, L.M.-F.; data curation, L.M.-F.; writing—original draft, L.M.-F. and R.E.-T.; writing—review and editing, L.M.-F., R.E.-T., J.T. and J.I.B.; supervision, R.E.-T. and J.T.; project administration, R.E.-T. and J.T.; funding acquisition, J.I.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article and supplementary material.

Acknowledgments

L. Manen-Freixa acknowledges IQS School of Engineering for the predoctoral grant received.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ruffolo, R.R. Why Has R&D Productivity Declined in the Pharmaceutical Industry? Expert Opin. Drug Discov. 2006, 1, 99–102. [Google Scholar] [CrossRef] [PubMed]
  2. Scannell, J.W.; Blanckley, A.; Boldon, H.; Warrington, B. Diagnosing the Decline in Pharmaceutical R&D Efficiency. Nat. Rev. Drug Discov. 2012, 11, 191–200. [Google Scholar] [CrossRef] [PubMed]
  3. Oprea, T.I.; Gottfries, J. Chemography: The Art of Navigating in Chemical Space. J. Comb. Chem. 2001, 3, 157–166. [Google Scholar] [CrossRef] [PubMed]
  4. Bohacek, R.S.; Mcmartin, C.; Guida, W.C. The Art and Practice of Structure-Based Drug Design: A Molecular Modeling Perspective. Med. Res Rev. 1996, 16, 3–50. [Google Scholar] [CrossRef]
  5. Miller, M.A. Chemical Database Techniques in Drug Discovery. Nat. Rev. Drug Discov. 2002, 1, 220–227. [Google Scholar] [CrossRef]
  6. Willett, P.; Barnard, J.M.; Downs, G.M. Chemical Similarity Searching. J. Chem. Inf. Comput. Sci. 1998, 38, 983–996. [Google Scholar] [CrossRef]
  7. Vogt, M. How Do We Optimize Chemical Space Navigation? Expert Opin. Drug Discov. 2020, 15, 523–525. [Google Scholar] [CrossRef]
  8. Snarey, M.; Terrett, N.K.; Willett, P.; Wilton, D.J. Comparison of Algorithms for Dissimilarity-Based Compound Selection. J. Mol. Graph. Model. 1997, 15, 372–385. [Google Scholar] [CrossRef]
  9. Dunbar, J.B. Cluster-Based Selection. Perspect. Drug Discov. Des. 1997, 7, 51–63. [Google Scholar] [CrossRef]
  10. Rabal, O.; Pascual, R.; Borrell, J.I.; Teixidó, J. Cell-Integral-Diversity Criterion: A Proposal for Minimizing Cluster Artifact in Cell-Based Selections. J. Chem. Inf. Model. 2007, 47, 1886–1896. [Google Scholar] [CrossRef]
  11. Mason, J.S.; Picket, S.D. Partition-Based Selection. Perspect. Drug Discov. Des. 1997, 7, 85–114. [Google Scholar] [CrossRef]
  12. Gorse, A.-D. Diversity in Medicinal Chemistry Space. Curr. Top. Med. Chem. 2006, 6, 3–18. [Google Scholar] [CrossRef]
  13. Gilad, Y.; Nadassy, K.; Senderowitz, H. A Reliable Computational Workflow for the Selection of Optimal Screening Libraries. J. Cheminform. 2015, 7, 61. [Google Scholar] [CrossRef] [PubMed]
  14. Lipinski, C.; Hopkins, A. Navigating Chemical Space for Biology and Medicine. Nature 2004, 432, 855–861. [Google Scholar] [CrossRef] [PubMed]
  15. Oprea, T.I. Chemical Space Navigation in Lead Discovery. Curr. Opin. Chem. Biol. 2002, 6, 384–389. [Google Scholar] [CrossRef]
  16. Reymond, J.; Van Deursen, R.; Blum, L.C.; Ruddigkeit, L. Chemical Space as a Source for New Drugs. MedChemComm 2010, 1, 30–38. [Google Scholar] [CrossRef]
  17. Zabolotna, Y.; Lin, A.; Horvath, D.; Marcou, G.; Volochnyuk, D.M.; Varnek, A. Chemography: Searching for Hidden Treasures. J. Chem. Inf. Model. 2021, 61, 179–188. [Google Scholar] [CrossRef]
  18. Lin, A.; Horvath, D.; Afonina, V.; Marcou, G.; Reymond, J.L.; Varnek, A. Mapping of the Available Chemical Space versus the Chemical Universe of Lead-Like Compounds. ChemMedChem 2018, 13, 540–554. [Google Scholar] [CrossRef]
  19. Fink, T.; Raymond, J.L. Virtual Exploration of the Chemical Universe up to 11 Atoms of C, N, O, F: Assembly of 26.4 Million Structures (110.9 Million Stereoisomers) and Analysis for New Ring Systems, Stereochemistry, Physicochemical Properties, Compound Classes, and Drug Discove. J. Chem. Inf. Model. 2007, 47, 342–353. [Google Scholar] [CrossRef]
  20. Ertl, P. Cheminformatics Analysis of Organic Substituents: Identification of the Most Common Substituents, Calculation of Substituent Properties, and Automatic Identification of Drug-like Bioisosteric Groups. J. Chem. Inf. Comput. Sci. 2003, 43, 374–380. [Google Scholar] [CrossRef]
  21. Warr, W.A.; Nicklaus, M.C.; Nicolaou, C.A.; Rarey, M. Exploration of Ultralarge Compound Collections for Drug Discovery. J. Chem. Inf. Model. 2022, 62(9), 2021–2034. [Google Scholar] [CrossRef] [PubMed]
  22. Pascual, R.; Mateu, M.; Gasteiger, J.; Borrell, J.I.; Teixidó, J. Design and Analysis of a Combinatorial Library of HEPT Analogues: Comparison of Selection Methodologies and Inspection of the Actually Covered Chemical Space. J. Chem. Inf. Comput. Sci. 2003, 43, 199–207. [Google Scholar] [CrossRef]
  23. Puig-De-La-Bellacasa, R.; Giménez, L.; Pettersson, S.; Pascual, R.; Gonzalo, E.; Esté, J.A.; Clotet, B.; Borrell, J.I.; Teixidó, J. Diverse Combinatorial Design, Synthesis and in Vitro Evaluation of New HEPT Analogues as Potential Non-Nucleoside HIV-1 Reverse Transcription Inhibitors. Eur. J. Med. Chem. 2012, 54, 159–174. [Google Scholar] [CrossRef] [PubMed]
  24. Barnard, J.M. A Comparison of Different Approaches to Markush Structure Handling. J. Chem. Inf. Comput. Sci. 1991, 31, 64–68. [Google Scholar] [CrossRef]
  25. Simmons, E.S. The Grammar of Markush Structure Searching: Vocabulary vs.. Syntax. World Pat. Inf. 1991, 13, 173. [Google Scholar] [CrossRef]
  26. Free, S.M.; Wilson, J.W. A Mathematical Contribution to Structure-Activity Studies. J. Med. Chem. 1964, 7, 395–399. [Google Scholar] [CrossRef] [PubMed]
  27. Kubinyi, H. Free Wilson Analysis. Theory, Applications and Its Relationship to Hansch Analysis. Quant. Struct. Relatsh. 1988, 7, 121–133. [Google Scholar] [CrossRef]
  28. Chen, H.; Carlsson, L.; Eriksson, M.; Varkonyi, P.; Norinder, U.; Nilsson, I. Beyond the Scope of Free-Wilson Analysis: Building Interpretable QSAR Models with Machine Learning Algorithms. J. Chem. Inf. Model. 2013, 53, 1324–1336. [Google Scholar] [CrossRef] [PubMed]
  29. Craig, P.N. Comparison of the Hansch and Free-Wilson Approaches to Structure-Activity Correlation. Biol. Correl. Hansch Approach 1974, 114, 115–129. [Google Scholar]
  30. Pascual, R.; Borrell, J.I.; Teixidó, J. Analysis of Selection Methodologies for Combinatorial Library Design. Mol. Divers. 2003, 6, 121–133. [Google Scholar] [CrossRef]
  31. Willett, P. Dissimilarity-Based Algorithms for Selecting Structurally Diverse Sets of Compounds. J. Comput. Biol. 1999, 6, 447–457. [Google Scholar] [CrossRef] [PubMed]
  32. Chen, G.; Zheng, S.; Luo, X.; Shen, J.; Zhu, W.; Liu, H.; Gui, C.; Zhang, J.; Zheng, M.; Chum, M.P.; et al. Focused Combinatorial Library Design Based on Structural Diversity, Druglikeness and Binding Affinity Score. J. Comb. Chem. 2005, 7, 398–406. [Google Scholar] [CrossRef] [PubMed]
  33. Meinl, T.; Ostermann, C.; Berthold, M.R. Maximum-Score Diversity Selection for Early Drug Discovery. J. Chem. Inf. Model. 2011, 51, 237–247. [Google Scholar] [CrossRef] [PubMed]
  34. Holliday, J.D.; Ranade, S.S.; Willett, P.A. Fast Algorithm for Selecting Sets of Dissimilar Molecules from Large Chemical Databases. Quant. Struct.−Act. Relat. 1995, 14, 501. [Google Scholar] [CrossRef]
  35. Willett, P. Similarity-Based Virtual Screening Using 2D Fingerprints. Drug Discov. Today 2006, 11, 1046–1053. [Google Scholar] [CrossRef]
  36. Groza, V.; Udrescu, M.; Bozdog, A.; Udrescu, L. Drug Repurposing Using Modularity Clustering in Drug-Drug Similarity Networks Based on Drug–Gene Interactions. Pharmaceutics 2021, 13, 2117. [Google Scholar] [CrossRef]
  37. Fakhoury, S.A.; Lee, H.T.; Reed, J.E.; Schlosser, K.M.; Sexton, K.E.; Tecle, H.; Winters, R.T. 4-phenylamino-quinazolin-6-yl-amides. U.S. Patent No. 7,772.243 B2, 2010. [Google Scholar]
  38. Mcchesney, J.; Nanayyakkra, D.N.; Bartlerr, M.; Ager, A.L. 8-aminoquinolines. U.S. Patent No. 2001/0007031 A1, 5 July 2001. [Google Scholar]
  39. Mascitti, V.; Collman, B.M. Dioxa-bicyclo[3.2.1]octane-2,3,4-triol dervatives. U.S. Patent No. 2010/0056618 A1, 4 March 2010. [Google Scholar]
  40. Coates, D.A.; Lawrence, M.; Knobeloch, J.M.; De Dios Magna, A.; De Prado Gonzalez, A.; Del Prado Catalina, M.F.; Garcia Paredes, M.C.; Martin de la Nava, M.; Martin Ortega Finger, M.D.; Martinez Perez, J.A.; et al. Protein Kinase Inhibitors. U.S. Patent No. 2010/0160340 A1, 24 June 2010. [Google Scholar]
  41. Heubach, G. 5-methyl-isoxazole-4-carboxylic Acid Anilides. U.S. Patent No. 4,087,535, 2 May 1978. [Google Scholar]
  42. Blumbergs, P.; LaMontagne, M.P. 4-methyl-5-(unsubstituted and substituted phenoxy)-2,6-dimethoxy-8- (aminoalkylamino) quinolines. U.S. Patent 4,617,394, 14 October 1986. [Google Scholar]
  43. Kuroita, T.; Sakamoto, H.; Ojima, M. Benzimidazole derivative and use as a receptor antagonist. U.S. Patent No. 7,572,920 B2, 11 August 2009. [Google Scholar]
  44. Meier, R. Fluorinated Benzyl Triazole Compounds. Patent No. EP 0199262, 15 April 1986. [Google Scholar]
  45. Orange Book: Approved Drug Products with Therapeutic Equivalence Evaluations. Available online: https://www.accessdata.fda.gov/scripts/cder/ob/index.cfm (accessed on 20 July 2021).
  46. Kim, S.; Chen, J.; Cheng, T.; Gindulyte, A.; He, J.; He, S.; Li, Q.; Shoemaker, B.A.; Thiessen, P.A.; Yu, B.; et al. PubChem in 2021: New Data Content and Improved Web Interfaces. Nucleic Acids Res. 2019, 49, D1388–D1395. [Google Scholar] [CrossRef]
  47. Dutta, A.K.; Avery, B.A.; Wyandt, C.M. Development and Validation of a Stability-Indicating Reversed-Phase High Performance Liquid Chromatography Method for NPC 1161C, a Novel 8-Aminoquinoline Anti-Malarial Drug. J. Chromatogr. A 2006, 1110, 35–45. [Google Scholar] [CrossRef]
  48. Strube, R.E.; LaMontagne, M.P. 4-methyl-5-(Unsubstituted and Substituted Phenoxy)-6-methoxy-8-(aminoalkylamino)quinolines. U.S. Patent No. 4,431,807, 14 October 1984. [Google Scholar]
  49. Nodiff, E.A.; Tanabe, K.; Chen, E.H.; Saggiomo, A.J. Modifications of Primaquine as Antimalarials. 3. 5-Phenoxy Derivatives of Primaquine. J. Med. Chem. 1982, 4816, 1097–1101. [Google Scholar] [CrossRef]
  50. Azad, C.S.; Saxena, M.; Siddiqui, A.J.; Bhardwaj, J.; Puri, S.K.; Dutta, G.P.; Anand, N.; Saxena, A.K. Synthesis of Primaquine Glyco-Conjugates as Potential Tissue Schizontocidal Antimalarial Agents. Chem. Biol. Drug Des. 2017, 90, 254–261. [Google Scholar] [CrossRef]
  51. LaMontagne, M.P.; Blumbergs, P.; Strube, R.E. Antimalarials. 14. 5-(Aryloxy)-4-Methylprimaquine Analogues. A Highly Effective Series of Blood and Tissue Schizonticidal Agents. J. Med. Chem. 1982, 25, 1094–1097. [Google Scholar] [CrossRef] [PubMed]
  52. Ashley, E.A.; Phyo, A.P. Drugs in Development for Malaria. Drugs 2018, 78, 861–879. [Google Scholar] [CrossRef] [PubMed]
  53. McChesney, J.; Nanayyakkra, D.N.; Bartlerr, M.; Ager, A.L. 8-aminoquinolines. Patent No. WO 97/36590, 28 March 1997. [Google Scholar]
  54. Chaurasiya, N.D.; Ganesan, S.; Nanayakkara, N.P.D.; Dias, L.R.S.; Walker, L.A.; Tekwani, B.L. Inhibition of Human Monoamine Oxidase A and B by 5-Phenoxy 8-Aminoquinoline Analogs. Bioorganic Med. Chem. Lett. 2012, 22, 1701–1704. [Google Scholar] [CrossRef] [PubMed]
  55. Brown, N. Chemoinformatics—An Introduction for Computer Scientists. ACM Comput. Surv. 2009, 41, 1–38. [Google Scholar] [CrossRef]
  56. Halkidi, M.; Vazirgiannis, M.; Balislakis, V. Quality Scheme Assessment in the Clustering Process. Lect. Notes Comput. Sci 2000, 1910, 265–276. [Google Scholar] [CrossRef]
  57. Rousseeuw, P.J. Silhouettes: A Graphical Aid to the Interpretation and Validation of Cluster Analysis. J. Comput. Appl. Math. 1987, 20, 53–65. [Google Scholar] [CrossRef]
  58. Mur, A.; Dormido, R.; Duro, N.; Dormido-Canto, S.; Vega, J. Determination of the Optimal Number of Clusters Using a Spectral Clustering Optimization. Expert Syst. Appl. 2016, 65, 304–314. [Google Scholar] [CrossRef]
  59. Hruschka, E.R.; Covões, T.F. Feature Selection for Cluster Analysis: An Approach Based on the Simplified Silhouette Criterion. In Proceedings of the International Conference on Computational Intelligence for Modelling, Control and Automation and CIMCA 2005 International Conference on Intelligent Agents, Web Technologies and Internet Commerce, Vienna, Austria, 28–30 November 2005; Volume 1, pp. 32–37. [Google Scholar] [CrossRef]
  60. ChemAxon—Software Solutions and Devices for Chemistry & Biology. Available online: https://chemaxon.com/ (accessed on 20 July 2021).
  61. Kim, S.; Thiessen, P.A.; Cheng, T.; Yu, B.; Bolton, E.E. An Update on PUG-REST: RESTful Interface for Programmatic Access to PubChem. Nucleic Acids Res. 2018, 46, W563–W570. [Google Scholar] [CrossRef]
  62. Molecular Operating Environment (MOE); Chemical Computing Group ULC: Montreal, QC, Canada, 2020.
  63. Pedregosa, F.; Varoquaux, G.; Gramfort, A.; Michel, V.; Thirion, B.; Grisel, O.; Blondel, M.; Prettenhofer, P.; Weiss, R.; Dubourg, V.; et al. Scikit-Learn: Machine Learning in Python. J. Mach. Learn. Res. 2011, 12, 2825–2830. [Google Scholar]
  64. Jamois, E.A.; Hassan, M.; Waldman, M. Evaluation of Reagent-Based and Product-Based Strategies in the Design of Combinatorial Library Subsets. J. Chem. Inf. Comput. Sci. 2000, 40, 63–70. [Google Scholar] [CrossRef]
  65. Waldman, M.; Li, H.; Hassan, M. Novel Algorithms for the Optimization of Molecular Diversity of Combinatorial Libraries. J. Mol. Graph. Model. 2000, 18, 412–426. [Google Scholar] [CrossRef]
Figure 1. Markush structure reported for the seven drugs under study.
Figure 1. Markush structure reported for the seven drugs under study.
Pharmaceuticals 15 01159 g001
Figure 2. For each subfigure, in the first row, there is the graphical representation of the scatter plot contour of MCL (pale yellow), BD (blue) and BCL (green) libraries when projected on PC1-PC2 (left), PC2-PC3 (middle) and PC3-PC1 (right) planes. In the second row, we show the density histogram for the considered libraries along PC1 (left), PC2 (middle) and PC3 (right). This scheme is repeated for all patents: (A) Dacomitinib, (B) Abemaciclib, (C) Tafenoquine, (D) Ertugliflozin, (E) Rufinamide, (F) Azilsartan Medoxomil and (G) Leflunomide.
Figure 2. For each subfigure, in the first row, there is the graphical representation of the scatter plot contour of MCL (pale yellow), BD (blue) and BCL (green) libraries when projected on PC1-PC2 (left), PC2-PC3 (middle) and PC3-PC1 (right) planes. In the second row, we show the density histogram for the considered libraries along PC1 (left), PC2 (middle) and PC3 (right). This scheme is repeated for all patents: (A) Dacomitinib, (B) Abemaciclib, (C) Tafenoquine, (D) Ertugliflozin, (E) Rufinamide, (F) Azilsartan Medoxomil and (G) Leflunomide.
Pharmaceuticals 15 01159 g002
Figure 3. Hit-to-lead process comparative workflow between the traditional methodology (A) and the approach suggested by the authors after the reported study (B).
Figure 3. Hit-to-lead process comparative workflow between the traditional methodology (A) and the approach suggested by the authors after the reported study (B).
Pharmaceuticals 15 01159 g003
Figure 4. t-SNE representation of Tafenoquine’s chemical space, depicting the MCL compounds (grey), BD subset (orange), compounds selected by the rational selection of NBD analogs (green) and the original hit (black triangle).
Figure 4. t-SNE representation of Tafenoquine’s chemical space, depicting the MCL compounds (grey), BD subset (orange), compounds selected by the rational selection of NBD analogs (green) and the original hit (black triangle).
Pharmaceuticals 15 01159 g004
Figure 5. Considering the selection of 9 compounds among a database consisting of 50 molecules (N = 50, represented as dots). (A) After applying a clustering or partition-based method, the chemical space is divided into 16 partitions (clusters or occupied bins, K = 16). (B) If the 9 selected compounds (depicted as stars) are distributed in 7 clusters (koc = 7, in orange), this corresponds to 43.8% SC. (C) Moreover, in these colored clusters are involved 28 samples of the database (noc). Hence, the 7 selected molecules are representative of 28 out of the 50 compounds included in the database, which represents 56.0% PC.
Figure 5. Considering the selection of 9 compounds among a database consisting of 50 molecules (N = 50, represented as dots). (A) After applying a clustering or partition-based method, the chemical space is divided into 16 partitions (clusters or occupied bins, K = 16). (B) If the 9 selected compounds (depicted as stars) are distributed in 7 clusters (koc = 7, in orange), this corresponds to 43.8% SC. (C) Moreover, in these colored clusters are involved 28 samples of the database (noc). Hence, the 7 selected molecules are representative of 28 out of the 50 compounds included in the database, which represents 56.0% PC.
Pharmaceuticals 15 01159 g005
Table 1. Information related to drugs under study. The year of approval (Year), the name of the patent applicant (Applic), the name of the disease (Disease), the patent status (Status), the number of analogs included in MCL (N), in Bibliographical Data (NBD) and in Bibliographical Combinatorial Library (NBCL) and the root square of N.
Table 1. Information related to drugs under study. The year of approval (Year), the name of the patent applicant (Applic), the name of the disease (Disease), the patent status (Status), the number of analogs included in MCL (N), in Bibliographical Data (NBD) and in Bibliographical Combinatorial Library (NBCL) and the root square of N.
DrugYearApplicDiseaseStatus N N B D N B C L N
Dacomitinib2018Pfizer Inc.Metastatic Non-Small-Cell Lung CancerONP16,53060798129
Abemaciclib2017Eli LillyBreast CancerONP45,69641736214
Tafenoquine2018GSKMalariaONP25,47258600160
Ertugliflozin2017MerckDiabetesONP14,1942156120
Rufinamide2008ESAILennox-Gastaut SyndromeOFP89592214495
Azilsartan Medoxomil2011TakedaHypertensionONP11104934
Leflunomide1998SanofiRheumatoid ArthritisOFP5641114284476
Table 2. Population distribution per cluster (k) analysis in Tafenoquine’s analogs’ chemical space, varying the clustering standard sizes. The mean value ( x ¯ ) of the population frequencies per cluster and its subsequent standard deviation (σ), along with the ratio of singletons to number of clusters (%s), are considered in discussion.
Table 2. Population distribution per cluster (k) analysis in Tafenoquine’s analogs’ chemical space, varying the clustering standard sizes. The mean value ( x ¯ ) of the population frequencies per cluster and its subsequent standard deviation (σ), along with the ratio of singletons to number of clusters (%s), are considered in discussion.
k = N 2 = 80 k = N = 160 k = 10 · N = 1600
x ¯ ± σ % s x ¯ ± σ % s x ¯ ± σ % s
HRC single318 ± 258228.7159 ± 181626.216 ± 6318.1
HRC complete318 ± 1840.0159 ± 930.016 ± 110.0
HRC median318 ± 8470.0159 ± 3330.016 ± 235.7
HRC average322 ± 3480.0159 ± 1510.016 ± 140.7
HRC centroid322 ± 22947.6159 ± 14366.916 ± 257.4
HRC Ward318 ± 1130.0159 ± 540.016 ± 60.0
KMN318 ± 720.0159 ± 380.016 ± 50.0
KMED318 ± 1250.0159 ± 640.016 ± 70.2
Binning411 ± 2690.0209 ± 1650.020 ± 2510.8
OV binning509 ± 5550.0173 ± 2010.720 ± 259.8
Table 3. Comparison of SC and PC values obtained by random selection and BD compounds when dividing the MCL of each patent in NBD clusters for Dacomitinib, Abemacilib and Tafenoquine databases.
Table 3. Comparison of SC and PC values obtained by random selection and BD compounds when dividing the MCL of each patent in NBD clusters for Dacomitinib, Abemacilib and Tafenoquine databases.
DacomitinibAbemaciclibTafenoquine
BDRandomBDRandomBDRandom
HRC averageSC26.741.122.035.117.246.8
PC74.280.166.285.912.878.4
HRC
Complete
SC50.054.034.152.025.959.3
PC69.072.844.974.019.367.6
HRC WardSC41.753.534.158.617.261.8
PC71.573.243.968.818.865.3
KMNSC48.357.634.160.322.462.5
PC67.369.439.167.421.064.7
KMEDSC55.057.634.159.124.161.4
PC68.068.844.468.321.865.9
OV binningSC25.042.928.151.824.049.2
PC59.178.346.284.035.482.9
Table 4. Comparison of SC and PC values obtained by random selection and BD compounds when dividing the MCL of each patent in NBD clusters for Ertugliflozin, Rufinamide, Azilsartan Medoxomil and Leflunomide databases.
Table 4. Comparison of SC and PC values obtained by random selection and BD compounds when dividing the MCL of each patent in NBD clusters for Ertugliflozin, Rufinamide, Azilsartan Medoxomil and Leflunomide databases.
ErtugliflozinRufinamideAzilsartan
Medoxomil
Leflunomide
BDRandom BDRandomBDRandomBDRandom
HRC averageSC9.519.836.430.650.026.427.250.1
PC83.888.46.486.21.398.711.676.4
HRC
Complete
SC23.838.722.743.650.026.326.355.3
PC59.478.84.079.71.398.717.971.5
HRC WardSC19.057.713.670.850.041.015.860.6
PC26.969.76.267.182.787.015.766.8
KMNSC23.858.718.258.950.052.115.860.8
PC23.569.37.969.459.984.813.166.7
KMEDSC19.057.727.361.650.068.318.459.8
PC29.870.827.466.646.868.820.867.6
OV binningSC13.365.628.650.325.056.626.453.0
PC18.282.29.690.11.780.232.583.0
Table 5. SC and PC results for a number of k = N clusters for Tafenoquine analogs.
Table 5. SC and PC results for a number of k = N clusters for Tafenoquine analogs.
HRC AverageHRC CompleteOV Binning
TafenoquineSelection SizeSCPCSCPCSCPC
BD588.86.411.99.713.617.0
Random BD5827.045.429.137.127.153.6
BCL60038.833.233.828.336.049.3
Random BCL60080.495.690.696.174.796.1
Random ( N )16050.974.857.769.148.279.8
Table 6. SC and PC results for a number of k = N clusters for Dacomitinib analogs.
Table 6. SC and PC results for a number of k = N clusters for Dacomitinib analogs.
HRC AverageHRC CompleteOV Binning
DacomitinibSelection sizeSCPCSCPCSCPC
BD6020.950.327.142.116.852.0
Random BD6031.254.834.446.430.270.6
BCL79879.192.186.890.964.392.3
Random BCL79886.998.495.998.881.998.1
Random ( N )12949.474.756.269.945.183.5
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Manen-Freixa, L.; Borrell, J.I.; Teixidó, J.; Estrada-Tejedor, R. Deconstructing Markush: Improving the R&D Efficiency Using Library Selection in Early Drug Discovery. Pharmaceuticals 2022, 15, 1159. https://doi.org/10.3390/ph15091159

AMA Style

Manen-Freixa L, Borrell JI, Teixidó J, Estrada-Tejedor R. Deconstructing Markush: Improving the R&D Efficiency Using Library Selection in Early Drug Discovery. Pharmaceuticals. 2022; 15(9):1159. https://doi.org/10.3390/ph15091159

Chicago/Turabian Style

Manen-Freixa, Leticia, José I. Borrell, Jordi Teixidó, and Roger Estrada-Tejedor. 2022. "Deconstructing Markush: Improving the R&D Efficiency Using Library Selection in Early Drug Discovery" Pharmaceuticals 15, no. 9: 1159. https://doi.org/10.3390/ph15091159

APA Style

Manen-Freixa, L., Borrell, J. I., Teixidó, J., & Estrada-Tejedor, R. (2022). Deconstructing Markush: Improving the R&D Efficiency Using Library Selection in Early Drug Discovery. Pharmaceuticals, 15(9), 1159. https://doi.org/10.3390/ph15091159

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop