Fernanda I. Saldívar-González
Gabriel Navarrete-Vázquez
José L. Medina-Franco- 1Department of Pharmacy, DIFACQUIM Research Group, School of Chemistry, Universidad Nacional Autónoma de México, Mexico City, Mexico
- 2Facultad de Farmacia, Universidad Autónoma Del Estado de Morelos, Cuernavaca, Morelos, Mexico
Virtual small molecule libraries are valuable resources for identifying bioactive compounds in virtual screening campaigns and improving the quality of libraries in terms of physicochemical properties, complexity, and structural diversity. In this context, the computational-aided design of libraries focused against antidiabetic targets can provide novel alternatives for treating type II diabetes mellitus (T2DM). In this work, we integrated the information generated to date on compounds with antidiabetic activity, advances in computational methods, and knowledge of chemical transformations available in the literature to design multi-target compound libraries focused on T2DM. We evaluated the novelty and diversity of the newly generated library by comparing it with antidiabetic compounds approved for clinical use, natural products, and multi-target compounds tested in vivo in experimental antidiabetic models. The designed libraries are freely available and are a valuable starting point for drug design, chemical synthesis, and biological evaluation or further computational filtering. Also, the compendium of 280 transformation rules identified in a medicinal chemistry context is made available in the linear notation SMIRKS for use in other chemical library enumeration or hit optimization approaches.
1 Introduction
Type 2 diabetes mellitus (T2DM) is a metabolic disorder characterized by hyperglycemia caused by defects in insulin secretion and/or action due to a complex network of pathological conditions (Galicia-Garcia et al., 2020). Currently, T2DM is one of the diseases with the highest socio-health impact and prevalence worldwide. Although pharmacotherapeutic options include different mechanisms of action, they are limited by side effects and lack of blood glucose control in the diabetic population (Shah et al., 2021). Another problem is that patients with this disease are prone to polypharmacy, which increases the risk of adverse effects and makes it difficult for patients to adhere to their treatment and receive proper follow-up from healthcare professionals (Dobrică et al., 2019). For this reason, new biological targets have been explored in multi-target approaches (Makhoba et al., 2020). Similarly, virtual libraries focused on single therapeutic targets have been developed using various computational approaches and their application in multi-target approaches is emerging.
Recent advances in computational methods and the incorporation of synthetic expert knowledge have inspired research groups to develop de novo and “make-on-demand” chemical libraries (Walters, 2019). Several companies use the so-called “novel molecular matter” in early-phase drug discovery (Korn et al., 2023). Specifically in the design of antidiabetic compounds, Otava released the chemical structures of ten libraries focused on DM-related targets designed under ligand- and structure-based approaches or combinations of both (OTAVAchemicals, Ltd. - synthetic organic compounds for research and drug discovery, n.d.). ChemDiv developed methods for screening diverse and highly specialized focused compounds. Recently, ChemDiv released an “Annotated space library” with more than 18,000 chemical compounds covering 38 validated targets (including targets for T2DM) across 900 drugs launched in the last 10 years. Academic groups have also generated virtual libraries focused on T2DM. For example, Chen et al. used a generative method to design compounds targeted for GPR40 (Chen et al., 2021). To our knowledge, no in silico multi-target libraries have been disclosed for T2DM. However, there is published information on active compounds and pharmacophore models that can guide the design of multi-target compounds for T2DM (Artasensi et al., 2020; Lillich et al., 2020; Tassopoulou et al., 2022).
Figure 1 shows the structures of pharmacophores and multi-target compounds studied in vivo models for T2DM and metabolic syndrome (MetS). For example, dual peroxisome proliferator-activated receptor (PPAR) α/γ agonists can improve insulin sensitivity and reduce triglyceride and blood glucose levels without the PPARγ-related weight gain since the latter effect is balanced by PPARα agonistic activity (Balakumar et al., 2019). Attempts towards developing dual agonists for PPAR α/γ yielded promising molecules that have reached clinical trials. However, many of these compounds have failed due to significant side effects. Compounds such as 1 (MHY908) and 2 (LT175) continue to be explored, as they have been shown to have beneficial effects on blood glucose and insulin resistance in animal models of T2DM (Gilardi et al., 2014; Park et al., 2016).
FIGURE 1. Pharmacophores and chemical structures of multi-target compounds studied in vivo models for T2DM and MEtS. (A) PPARα/γ agonists, (B) PPARα, PPARγ, GPR40, AR and PTP1B (C) PPARγ/SUR agonists (D) PPARγ/FFA1 agonists.
The development of multi-target compounds for GPR40, PTP1B, AR, PPARα, and PPARγ may provide additional therapeutic benefits in preventing or delaying the development of diabetic complications (Domínguez-Mendoza et al., 2021). Navarrete’s group designed compounds 3 and 4 (Figure 1) that showed robust in vivo antihyperglycemic activity (Domínguez-Mendoza et al., 2021). Dual PPARγ/SUR agonists, such as compound 5 (Ibrahim et al., 2017a), can improve insulin sensitivity and stimulate insulin secretion simultaneously, making them an attractive therapeutic option for patients with T2DM who have insulin resistance and decreased insulin secretion. Another combination that has been explored is that of PPARy and FFA1, also known as GPR40. Dual PPARγ/FFA1 agonists (compounds 6 and 7, Darwish et al., 2018) can improve insulin sensitivity, stimulate insulin secretion, and lower blood glucose levels.
Table 1 shows multi-target compounds for T2DM based on pharmacophores and virtual screening approaches that have been reported and tested in vivo and can guide the design of new chemical libraries. The table includes information on the combination of targets, the number of lead compounds that hit all targets, and the implications or outcome in T2DM and MEtS. The table indicated that there are ninety-one multi-target compounds with in vivo antidiabetic activity. These have been tested in a total of twenty target combinations. The most studied biological target in multi-target approaches is PPARγ. The most successful target combination is PPARα/γ as several compounds directed to these two targets are in clinical studies such as ragaglitazar, imiglitazar, muraglitazar, tesaglitazar, naveglitazar, aleglitazar, saroglitazar, netoglitazone and lobeglitazone. In particular, PPARs have demonstrated clinical efficacy in metabolic diseases such as T2DM and those related to lipid metabolism, which has led to testing and optimization of different compounds such as fibrates and TZDs to simultaneously modulate multiple targets, with a synergistic effect on T2DM and MetS (Ammazzalorso et al., 2019). Other compounds under investigation or in any of the experimental phases by the FDA include elafibranor for PPARα/δ (Schattenberg et al., 2021), telmisartan and fimasartan for PPARγ/AT1 (Seo et al., 2022), lanifibranor and sodelglitazar for PPARα/δ/γ (Kamata et al., 2023) and licogliflozin and sotagliflozin for SGLT1-2, the latter already approved by the FDA in 2023 to reduce the risk of cardiovascular death and heart failure in adults with heart failure, T2DM, chronic kidney disease, and other cardiovascular risk factors (Packer, 2023).
TABLE 1. Multi-target ligands studied in vivo models for type II diabetes mellitus.
An attractive approach to exploring and expanding the chemical space around the first hit compounds of single and multi-target compounds is the computational generation of chemical libraries that can be used in virtual screening campaigns (Walters, 2019). Through chemical library enumeration, it is possible to find new bioactive compounds and generate therapeutic options for emerging and challenging molecular targets and complex diseases. It is also possible to control features such as library size, complexity, physicochemical properties, and structural diversity (Ruddigkeit et al., 2012). The goal is to help design high-quality analog series overcoming issues such as low potency, off-target activities, metabolic stability, or poor physicochemical properties for oral administration. One approach that can aid rigorous exploration of the chemistry around first hit compounds is using approaches based on transformation rules from empirical observation and systematic identification using chemoinformatics methods. A recent application of transformation rules was presented with DrugSpaceX, a database with more than 100 million compounds transformed from approved drug molecules (Yang et al., 2021). Although transformation rules are useful for generating in silico libraries, the list of rules currently available in the public domain is limited, in many cases due to the difficulty in collecting, curating, and annotating such information (Rarey et al., 2022). To address this problem, we compiled, organized, and made freely available an extensive list of transformation rules for generating compound libraries, as part of this work.
This study aimed to design a multi-target focused library for T2DM using a comprehensive set of chemical transformation rules herein collected, curated, and annotated. As shown in Figure 2, we evaluated the novelty and diversity of the focused library by comparing it with antidiabetic compounds approved for clinical use, natural products, and multi-target compounds reported with in vivo activity. To focus the library on attractive and synthetically feasible compounds, computational filters based on medicinal chemistry criteria were applied. Finally, virtual screening using molecular docking for PTP1B and AR was performed at Molecular Operating Environment (MOE) version 2022.02 (Chemical Computing Group (CCG), 2023) and the ADME-Tox properties of compounds with potential multi-target activity were calculated using ADMElab 2.0 (Xiong et al., 2021). Here, we selected PTP1B and AR considering the reference compounds used for enumeration and the currently available hypotheses and information from molecular dynamics models (Domínguez-Mendoza et al., 2021). In particular, this combination could be attractive to modulate insulin signaling, reducing insulin resistance and preventing or delaying diabetic complications such as nephropathies, neuropathies, and cardiomyopathies (Maccari and Ottanà, 2015). The significance of this work is twofold: making available a focused multi-targeted library for T2DM with full details of the methodology used to construct the compounds and making publicly available a general and extensive list of transformation rules to explore the chemical space of targets of therapeutic relevance.
FIGURE 2. Overview of the methodology implemented in this study to design a multi-target focused library for T2DM using a comprehensive set of chemical transformation rules. First, multi-target compounds evaluated for in vivo antidiabetic activity and transformation rules were retrieved from the literature. The latter were encoded into SMIRKS linear notation. Using MOE, the new chemical library was enumerated and compared to reference libraries of antidiabetic compounds to assess their chemical diversity. The compounds in the enumerated library were filtered using criteria as follows: QED > 0.67, SAScore ≤ 6, rsynth > 0.5 and zero structural alerts included in the RDKit molecule catalog filter node in KNIME. These compounds were further filtered in docking-based virtual screening with protein tyrosine phosphatase 1B (PTP1B) and aldose reductase (AR). Finally, the ADME-T properties of the best-scoring compounds in both targets were calculated.
2 Methods
2.1 Encoding chemical transformations in linear notation (SMIRKS format)
Transformation rules associated with hit optimization were retrieved from the literature. These transformation rules included modifications associated with molecular, physicochemical, pharmacological, ADME, safety, and toxicity parameters. Other rules are associated with structural diversity and bioisosteric changes. Most of the transformation rules found in the literature come from the addition, substitution or removal of a functional group. There are also rules that include cyclization and ring substitution by other rings or intramolecular hydrogen bonding groups. In total, 280 transformation rules were collected and converted into SMIRKS notation (Daylight theory: SMIRKS - A reaction transform language, 2023) using MarvinSketck (Chemaxon, 2023). The transformation rules were considered in a protonation state of 7.4. Table 2 summarizes examples of the rules identified and their SMIRKS. The Supplementary material includes the complete list of transformation rules and the literature reference.
TABLE 2. Examples of transformation rules retrieved from the literature.
2.2 Enumeration of a new multi-target focused library
Compounds 3 and 4 (Figure 1) were selected as reference compounds to enumerate the new library focused on multi-target compounds. These compounds have robust antihyperglycemic activity in vivo and molecular dynamics studies with PTP1B and AR provide relevant structure ligand-protein interaction information for structure-based optimization studies (Domínguez-Mendoza et al., 2021). The best predicted scoring conformations of these compounds docked with a crystallographic structure of the PTP1B protein obtained from the Protein Data Bank (Berman et al., 2000) (PDB ID: 4Y14 (Krishnan et al., 2015)) and that maintained protein-ligand interactions reported in literature were used as the basis for the enumeration of the new library. To compare the effect on the number of compounds and the molecular diversity generated, we used the 175 transformation rules integrated into the MedChem module of MOE version 2022.02 and added 273 transformation rules that we constructed as part of this study. Two iterations were used. Compounds that had a molecular weight (MW) < 630 and topological surface area (TPSA) between 40 and 150 Å were kept. The threshold values were established based on the minimum and maximum values of the multi-target compounds designed for T2DM.
2.3 Evaluation of the chemical diversity of the multi-target focused library
The compounds generated in Section 2.2 were compared in terms of physicochemical properties with antidiabetic compounds approved for clinical use retrieved from DrugBank (Wishart et al., 2018); antidiabetic compounds from medicinal plants obtained from DiaNat-DB (Madariaga-Mazón et al., 2021), and multi-target compounds for T2DM studied in vivo assays. Compounds in SMILES format can be found in the Supplementary material. Prior to analysis, each molecule was prepared using the open-source cheminformatics toolkits RDKit (Landrum, n.d.) and MolVS (MolVS: Molecule Validation and Standardization—MolVS 0.1.1 documentation, n.d.). Compounds were standardized and those containing multiple components were split, keeping the largest component. Compounds with valence errors or chemical elements other than H, B, C, N, O, F, Si, P, S, Cl, Se, Br, and I were removed. The remaining compounds were neutralized and reionized to subsequently generate a canonical tautomer. Finally, duplicated compounds were deleted and canonical simplified molecular-input line-entry system (SMILES) (ignoring stereochemistry information) were generated as molecular representation (Weininger, 1988). Table 3 summarizes the number and source of compounds used in the comparison.
TABLE 3. Reference compound datasets to assess the novelty and properties of the newly designed libraries.
2.3.1 Property rules and synthetic accessibility
To profile the generated molecules based on common property rules and synthetic accessibility, we calculated seven descriptors for each molecular entity with the RDKit library: molecular weight (MW), octanol–water partition coefficient (logP), number of hydrogen bond acceptors (HBA), number of hydrogen bond donors (HBD), topological polar surface area (TPSA), the number of rotatable bonds (RotB) and the distribution of quantitative estimate of drug-likeness (QED) (Bickerton et al., 2012).
Among the several computational scores available to estimate the synthetic accessibility of generated compounds, and based on previous comparisons of scores. (Skoraczyński et al., 2023), in this work, we calculated rsynth (a ligand-based score implemented in MOE, version 2022.02) and SASscore (Ertl and Schuffenhauer, 2009) as a structure-based approach.
2.3.2 Chemical multiverse: visual representation and analysis
To have a comprehensive analysis of the chemical space, we used three well-established visualization methods and different types of descriptors Noteworthy, analyzing the chemical space with distinct and complementary descriptors is crucial because each one provides a different and complementary perspective of the chemical space (aka, a “chemical multiverse” as discussed elsewhere (Medina-Franco et al., 2022). In this study, we used principal component analysis (PCA) and t-distributed stochastic neighbor embedding t-SNE based on six molecular properties of pharmacological interest, namely,; HBD, HBA, logP, MW, RB, TPSA. As a third method to characterize the chemical space, we used was the Tree MAP (TMAP) algorithm (Probst and Reymond, 2020) based on ECFP4 as a structural fingerprint. Additionally, the molecular shape of the compounds in the newly generated library was evaluated using the principal moments of inertia (PMI) graph, which was carried out by calculating the lowest energy conformation of each compound using MMFF94x as a force field in MOE, version 2022.22. Once the lowest energy conformer was calculated, values of normalized PMI ratios, npr1 (I1/I3) and npr2 (I1/I3), were determined in MOE. Then, npr1 and npr2 were plotted on a triangular graph with the vertices (0,1), (0.5,0.5), and (1,1) representing a perfect rod, disc, and sphere, respectively.
2.4 Filtering of compounds with relevance in pharmaceutical chemistry
Various filters (e.g., calculated descriptors) can be used to improve the selection of enumerated compounds, including the removal of compounds with undesirable pharmaceutical properties and molecules very difficult to synthesize. In this work, we chose to filter by SAScore, rsynth, QED, and the filters included in the RDKit molecule catalog filter node in KNIME. This node removes compounds with Pan-assay interference compounds (PAINS) (Baell and Holloway, 2010), unwanted functionality due to potential toxicity or unfavorable pharmacokinetic properties, and problematic functional groups. The compounds in the enumerated library were filtered using criteria as follows: QED > 0.67, SAScore ≤ 6, rsynth >0.5 and zero structural alerts. It should be noted that other cutoff values could be used for other applications.
2.5 Structure-based filtering
Although enumeration of compounds can be performed using the protein structure to obtain a score value in MOE, this can be computationally expensive considering the number of compounds that can potentially be generated. In addition, although the search algorithm can be selected, neither an algorithm to perform the rescoring nor a specific number of conformations to be obtained can be selected. For this reason, we followed a structure-based filtering to select the most promising compound subset as a multi-target library. In this study, compounds that showed a docking score better than the reference compounds for PTP1B (PDB ID: 4Y14 (Krishnan et al., 2015)), were further filtered considering docking with another therapeutic target relevant to T2DM: AR (PDB ID: 4XZH (Ruiz et al., 2015)). Of note, other relevant targets in T2DM can be considered during the design of the reference compounds. Here, we selected PTP1B and AR considering the current available hypothesis and model information based on molecular dynamics (Domínguez-Mendoza et al., 2021). Also, we want to point out that several other docking programs could be used (including post consensus scoring analysis). However, testing different docking programs and exploring multiple consensus analysis schemes is out of the scope of this study that is focused on proposing a general approach to design multi-target focused libraries.
The protein preparation of PDB ID: 4Y14 and 4XZH was made with default settings of the QuickPrep module of MOE v. 2022.02. This module carries out the following processes: calibration of the structure by protonation, addition of all the lacking hydrogen atoms, elimination of water molecules 4.5 Å farther from the protein, addition of missing amino acids residues, neutralization of the endpoints adjoining empty residues and energy minimization. We used AMBER14:EHT forcefield (ff14SB (Maier et al., 2015) as parameter for the energy minimization stage for the protein. The ligands were also prepared in MOE, we used MFF94x as forcefield. For docking, the receptor was considered rigid and the ligands flexible. We used the default settings for placement (method: triangle matcher, score function: London dG) and refinement (method: rigid receptor, score function: GBVI/WSA dG) (Vilar et al., 2008).
2.6 Prediction of ADME-Tox properties
The ADME-tox properties of the compounds that had successfully passed the filters described in Section 2.6 and Section 2.7 were evaluated using ADMETlab 2.0 platform (Xiong et al., 2021). This platform has been compared with other free online ADMET tools and has significant advantages (Dulsat et al., 2023). Based on these comparisons we used this tool because ADMETlab is a complete platform in terms of the large number of relevant parameters that can be predicted including elimination parameters such as clearance and half-life (t1/2). The latter two descriptors are particularly relevant in chronic diseases such as T2DM, where ideally compounds with a long half-life are sought to reduce the dose. Finally, ADMETlab 2.0 allows evaluating up to 500 compounds at the same time.
3 Results and discussion
3.1 Encoding chemical transformations in SMIRKS format
Based on the structures exploited in the multi-target pharmacophores for DMT2, 280 transformation rules were collected from the literature, of which 113 were bioisosteric modifications of the phenyl group (Subbaiah and Meanwell, 2021), 25 for the amide bond (Kumari et al., 2020), 36 for the carboxyl group (Bredael et al., 2022), 57 for the phosphate group, and 49 for other moieties including ester, alcohol, alkyl, aminophenyl, and nitro to name a few examples. Compared to the preloaded MOE transformations (version 2022.22), these included only 10 for the phenyl group, 7 for the amide bond and 4 for the carboxyl group. The remaining transformations in MOE mainly concern cyclization and substitution of rings by other rings or intramolecular hydrogen bonding groups. It is important to mention that some examples of the rule transformations present in Table 2 have been applied in the optimization of antidiabetic compounds. For example, Huang et al. reported the effect of bioisosteric replacement of a phenyl ring in the biphenyl moiety with cyclohexyl motif in a GPR40 agonist (Huang et al., 2019). Another bioisosteric replacement of the phenyl ring in an antidiabetic compound was reported by (Piemontese et al., 2015). In that work, the authors replaced a terminal phenyl ring in a dual PPAR α/γ agonist with the n-propyl oxime moiety. For its part, Dr. Navarrete-Vázquez’s group has proposed various series of thiazolidine-2,4-dione and barbituric acid derivatives with robust antihyperglycemic activity in vivo as acid bioisosteres (Hidalgo-Figueroa et al., 2013; Domínguez-Mendoza et al., 2021).
3.2 Enumeration of new compounds
Using the transformation rules preloaded in MOE and performing two iterations, 6,838 molecules were obtained from compound 3 and 1834 from compound 4. After adding the transformation rules collected from the literature and keeping two iterations, the number of compounds increased to 52,185 (from compound 3) and to 32,593 (from compound 4), after curation of the enumerated library. This large increase in the number of compounds was expected, as it followed an iterative process, and the transformations considered in this work are focused on moieties that contain the reference compounds.
3.3 Evaluation of the chemical diversity of the multi-target focused library
3.3.1 Property rules and synthetic accessibility
Figure 3 shows the distributions of each descriptor for the compounds of the generated library and the reference libraries using rain cloud plots (Allen et al., 2019). These plots allow visualization of the probability density and typical boxplot statistics such as median, mean, and confidence intervals. According to these plots, DiaNat-DB (329 antidiabetic compounds from medicinal plants) has a wider distribution in terms of properties of pharmaceutical importance. Since MW and TPSA values were used as filters to generate the new compounds, the distributions of these values for generated compounds are skewed and resemble the distribution of multi-target compounds and approved drugs. Comparing the plots of the other descriptors, we can see that the transformations used increased the range of properties such as LogP and HBD, HBA and chirality (Supplementary Table S1 in the Supplementary Material).
FIGURE 3. Rain cloud plots of the six physicochemical properties of pharmaceutical relevance: (A) molecular weight (MW), (B) topological polar surface area (TPSA), (C) partition coefficient octanol/water (log P), (D) number of rotatable bonds (RB), (E) hydrogen bond acceptors (HBA) and (F) hydrogen bond donors (HBD).
Figure 4 shows the distributions of the calculated descriptors to quantity synthetic accessibility and drug-likeness through QED. As can be seen in Figure 4, most of the generated compounds have a SAScore value of less than 6, indicating that they are, in principle, synthesizable (Ertl and Schuffenhauer, 2009). Regarding the QED value, which is a measure of drug-likeness based on the concept of desirability, a value greater than or equal to 0.67 represents attractive compounds, and the lower this value, the compounds are considered unattractive (QED 0.49∼0.67) and too complex (QED < 0.49). These reference values were established based on the calculated physicochemical properties of marketed oral drugs and published human data (Bickerton et al., 2012). When we compare the QED distribution of all antidiabetic compounds, we find that the compounds generated by transformation rules have a higher mean value (0.49) than the multi-target compounds obtained from the literature (0.46). The summary statistics of these plots can be found in the Supplementary material.
FIGURE 4. Distribution of (A) quantitative estimate of drug likeness (QED) and (B) SAscore of all the antidiabetic compounds contained in approved drugs (red) DiaNat database (green), multi-target compounds (orange) and multi-target compounds generated by medicinal chemistry transformation rules (purple). Vertical dashed lines represent the mean of the distributions.
3.3.2 Chemical multiverse of generated compounds
The concept of chemical multiverse (e.g., for the same compound datasets, different chemical spaces, each based on a different set of descriptors) was used to compare comprehensively the visual representation of the chemical space of the generated compounds with collections of reference compounds (Figure 5). The PCA of six physicochemical properties: MW, HBD, HBA, logP, TPSA and RB shows DiaNat-DB is the most diverse database in terms of physicochemical properties (Figure 5A). Using a non-linear method such as t-SNE and the same descriptors, we obtain a different visualization that allows us to visualize clusters or groups of data points and their relative proximities (5b). In terms of molecular fingerprints (ECFP4) (5c) and molecular shapes (5d), the library generated with transformation rules exhibited the largest structural and shape diversity, which is even larger than that of the DiaNat-DB database.
FIGURE 5. Visual representation of the chemical multiverse of antidiabetic compounds contained in approved drugs (red) DiaNat-DB (green), multi-target compounds (orange) and multi-target compounds generated by medicinal chemistry transformation rules (purple). (A) PCA of six physicochemical properties: MW, HBD, HBA, logP, TPSA and RB (B) t-SNE of six physicochemical properties (vide supra), (C) TMAP based on molecular ECFP4 fingerprint. (D) PMI space. Each corner on the triangular PMI plot indicates compounds with certain shape characteristics. The top left corner of the PMI represents compounds with rodlike shape, the top right corner represents compounds with spherical shape, and the bottom corner represents compounds with disc-like shape.
3.4 Filtering of compounds with relevance in pharmaceutical chemistry
To focus the library on attractive and synthetically feasible compounds, computational filters were applied as indicated in Section 2.6 of the Methods section. Table 4 summarizes the criteria used and the number of compounds remaining after applying the filters. As can be seen, the filter that considerably reduced the number of compounds was the QED value. The filtered compounds are characterized by not having chiral centers. If this is an important feature, you will need to consider it when applying the filters. Filtered compounds are listed in the Supplementary material.
TABLE 4. Filters applied and the number of compounds remaining for virtual screening.
3.5 Structure-based filtering
The compounds selected from the filtering described in Section 3.6 were subjected to a docking-based virtual screening with PTP1B and AR. The docking protocols are described in Section 2.6. The docking scores calculated with MOE for compounds 3 (PTP1B score: −7.04 kcal/mol and AR score = −8.84 kcal/mol) and 4 (PTP1B score: −7.91 kcal/mol and AR score = −8.98 kcal/mol) were used as cut-offs to select potential multi-target compounds. It was also checked whether the docking of these compounds reproduced the interactions previously reported in molecular dynamics studies (Domínguez-Mendoza et al., 2021). Supplementary Table S2 shows the score values and interaction plots of compounds 3 and 4 with PTP1B and AR. The first virtual screening with PTP1B yielded 1,655 compounds with a lower score value than the reference compounds: 1,543 from compound 3 and 112 from compound 4. The virtual screening hit compounds were docked with AR yielding 816 hit compounds: 792 were from compound 3 and 24 from compound 4. Figure 6 shows examples of designed compounds that have the potential to be used in multi-target approaches. The figure shows the parent structure (3 and 4), the transformation rules and the compounds designed with their corresponding calculated docking scores with PTP1B and AR. We also include the 3D docking models for the proposed compounds, and their overlap with the reference compounds. The docking results for PTP1B and AR of the 816 compounds can be found in the Supplementary material.
FIGURE 6. Examples of compounds selected from the multi-target virtual library, transformation rules used, and the calculated docking scores. 3D docking models for the proposed compounds (green) that overlap with the reference compounds (orange) for PTP1B and AR are included.
3.6 Prediction of ADME-Tox properties
Table 5 shows the average and standard deviation for different ADME-Tox properties of approved DMT2 drugs, multi-target compounds with reported in vivo activity and the 816 compounds that showed the potential to interact with PTP1B and AR. The values described to make an empirical decision are taken directly from the ADMElab 2.0 documentation (https://admetmesh.scbdd.com/explanatthat overlapion/index). As can be seen, the subset of compounds selected from the multi-target library generated in this work exhibit optimal adsorption and distribution properties. It is also noteworthy that in terms of metabolism, ADMETlab 2.0 predicts that the compounds generated in this work have a lower probability of being inhibitors of various CYPs compared to approved drugs and multi-target compounds that have a moderate to high probability of being inhibitors of CYP2C9, CYP2C19, and CYP3A4. Knowing this information is important because the inhibition of some CYPs is associated with the risk of hypoglycemia. Knowing this information is important because the inhibition of some CYPs is associated with the risk of hypoglycemia. For example, CYP2C9 metabolizes nateglinide, repaglinida, rosiglitazone, and most sulfonylureas, such as glibenclamide, glimepiride and glipizide (Holstein et al., 2012). Other examples are pioglitazone and repaginate metabolized with CYP3A4. ADMETlab 2.0 also predicts the probability of being substrates of various CYPs. This data and other properties calculated with ADMETlab 2.0 for each compound can be found in the Supplementary material. It is important to mention that the prediction of inhibitory activity in ADMETlab 2.0 is based on a dataset containing information on inhibitory activity of compounds obtained by high-throughput quantitative screening with an in vitro bioluminescence assay (Veith et al., 2009). However, in the description of these data, it is mentioned that the inhibitory activity in the assay may be due to compounds acting as substrates and that some weakly binding substrates may be classified as “inactive,” so that the predictions may need further confirmation.
In terms of excretion, clearance (CL) and half-life (T1/2) are important pharmacokinetic parameters that allow defining a drug´s dosing frequency. In the case of ADMETlab 2.0, the half-life is not measured in units of time. The output value is the probability of falling into category 1 (T1/2 ≤ 3). That is, the greater the probability of falling into category 1, the more likely the substance is to be classified as “poor” because its T1/2 would be lower (T1/2 ≤ 3). For antidiabetic drugs, the average CL is 4.4716 mL/min/kg (poor) and the T1/2 is 0.4080 (medium). The discrepancy in predictions could be due to the difference in models or datasets. In the case of the generated multi-target compounds, they may not be optimal for reducing the frequency of administration. Finally, although the probability of compounds being hERG blockers is reduced, the likelihood of drug-induced liver injury (DILI) would remain a challenge to optimize.
4 Conclusion
Designing multi-target compounds is an attractive approach to develop therapeutic treatments for complex diseases such as T2DM and MetS. Herein, we collected from the literature and analyzed ninety-one multi-target compounds for which in vivo antidiabetic activity has been reported, with a total of twenty target combinations. Following an enumeration based on transformation rules, we expanded the relevant chemical space of two of these multi-target hit compounds. More than 450 transformation rules were applied, of which 280 are made openly available to the scientific community. We concluded that the compounds generated with transformation rules have similar physicochemical properties to antidiabetic drugs and multi-target compounds reported in literature. Of the 84,778 generated compounds with valid structures, 85% are predicted to be synthetically feasible. The enumerated compounds are also attractive considering structural and shape diversity.
To focus on attractive and synthetically feasible compounds to perform virtual screening, various drug-likeness and quality filters were applied, yielding a multi-target virtual library with 2037 compounds. After a docking-based virtual screening with PTP1B and AR, 816 multi-target compounds were selected. Compounds in this library have favorable ADME properties, making the library an attractive source of promising candidates for further research and development.
In line with open and democratization of science, the newly designed multi-target focused library is freely available as a valuable source of starting points for chemical synthesis, biological evaluation, or further computational analysis such as virtual screening or reference libraries in library design.
Data availability statement
The original contributions presented in the study are publicly available. This data can be found here: https://figshare.com/projects/Design_of_a_multi-target_focused_library_for_antidiabetic_targets_using_a_comprehensive_set_of_chemical_transformation_rules/175194.
Author contributions
FS-G: Writing–original draft, Writing–review and editing, Conceptualization, Data curation, Formal Analysis, Investigation, Methodology, Software. GN-V: Conceptualization, Investigation, Writing–review and editing, Supervision. JM-F: Funding acquisition, Project administration, Resources, Supervision, Writing–original draft, Writing–review and editing, Formal Analysis, Conceptualization.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. We thank the support of DGAPA, UNAM, Programa de Apoyo a Proyectos de Investigación e Innovación Tecnológica (PAPIIT), Grant No. IV200121.
Acknowledgments
FS-G thanks Consejo Nacional de Humanidades, Ciencias y Tecnologías (CONAHCyT), Mexico for the Ph.D. scholarship 848061. The technical assistance and enthusiastic participation of Jenny Serrano López, Vanesa Garduño Ibarra, Leonardo Saravia Fuentes and Diana Prado Romero (School of Chemistry, UNAM), is greatly acknowledged. We also acknowledge the innovation space UNAM-HUAWEI the computational resources to use their supercomputer under project-7 “Desarrollo y aplicación de algoritmos de inteligencia artificial para el diseño de fármacos aplicables al tratamiento de diabetes mellitus y cáncer”.
Conflict of interest
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.
The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphar.2023.1276444/full#supplementary-material
Abbreviations
ADME, absorption, distribution, metabolism, and excretion; AR, aldose reductase; DM, diabetes mellitus; FFA1, Free fatty acid receptor 1; GPR40, G-protein-coupled receptor 40; MetS, metabolic syndrome; MOE, Molecular Operating Environment; MW, molecular weight; PCA, principal component analysis; PMI, principal moments of inertia; PPARs, peroxisome proliferator-activated receptors; PTP1B, protein tyrosine phosphatase 1B; QED, quantitative estimate of drug-likeness; RB, rotatable bonds; T2DM, Type II diabetes mellitus; TMAP, tree map; TPSA, topological polar surface area; t-SNE, t-distributed stochastic neighbor embedding; SMILES, Simplified Molecular Input Line Entry Specification, SMIRKS: line notation for specifying reaction transforms; SUR, Sulfonylurea receptor.
References
Allen, M., Poggiali, D., Whitaker, K., Marshall, T. R., and Kievit, R. A. (2019). Raincloud plots: a multi-platform tool for robust data visualization. Wellcome Open Res. 4, 63. doi:10.12688/wellcomeopenres.15191.1
Ammazzalorso, A., Maccallini, C., Amoia, P., and Amoroso, R. (2019). Multitarget PPARγ agonists as innovative modulators of the metabolic syndrome. Eur. J. Med. Chem. 173, 261–273. doi:10.1016/j.ejmech.2019.04.030
Artasensi, A., Pedretti, A., Vistoli, G., and Fumagalli, L. (2020). Type 2 diabetes mellitus: a Review of multi-target drugs. Molecules 25, 1987. doi:10.3390/molecules25081987
Baell, J. B., and Holloway, G. A. (2010). New substructure filters for removal of pan assay interference compounds (PAINS) from screening libraries and for their exclusion in bioassays. J. Med. Chem. 53, 2719–2740. doi:10.1021/jm901137j
Balakumar, P., Rose, M., Ganti, S. S., Krishan, P., and Singh, M. (2007). PPAR dual agonists: are they opening Pandora’s Box? Pharmacol. Res. 56, 91–98. doi:10.1016/j.phrs.2007.03.002
Balakumar, P., Mahadevan, N., and Sambathkumar, R. (2019). A contemporary overview of pparα/γ dual agonists for the management of diabetic dyslipidemia. Curr. Mol. Pharmacol. 12, 195–201. doi:10.2174/1874467212666190111165015
Benson, S. C., Pershadsingh, H. A., Ho, C. I., Chittiboyina, A., Desai, P., Pravenec, M., et al. (2004). Identification of telmisartan as a unique angiotensin II receptor antagonist with selective PPARgamma-modulating activity. Hypertension 43, 993–1002. doi:10.1161/01.HYP.0000123072.34629.57
Berman, H. M., Westbrook, J., Feng, Z., Gilliland, G., Bhat, T. N., Weissig, H., et al. (2000). The protein Data Bank. Nucleic Acids Res. 28, 235–242. doi:10.1093/nar/28.1.235
Bhattarai, B. R., Kafle, B., Hwang, J.-S., Khadka, D., Lee, S.-M., Kang, J.-S., et al. (2009). Thiazolidinedione derivatives as PTP1B inhibitors with antihyperglycemic and antiobesity effects. Bioorg. Med. Chem. Lett. 19, 6161–6165. doi:10.1016/j.bmcl.2009.09.020
Bhattarai, B. R., Kafle, B., Hwang, J.-S., Ham, S. W., Lee, K.-H., Park, H., et al. (2010). Novel thiazolidinedione derivatives with anti-obesity effects: dual action as PTP1B inhibitors and PPAR-γ activators. Bioorg. Med. Chem. Lett. 20, 6758–6763. doi:10.1016/j.bmcl.2010.08.130
Bickerton, G. R., Paolini, G. V., Besnard, J., Muresan, S., and Hopkins, A. L. (2012). Quantifying the chemical beauty of drugs. Nat. Chem. 4, 90–98. doi:10.1038/nchem.1243
Blöcher, R., Lamers, C., Wittmann, S. K., Merk, D., Hartmann, M., Weizel, L., et al. (2016). N-Benzylbenzamides: a novel merged scaffold for orally available dual soluble epoxide hydrolase/peroxisome proliferator-activated receptor γ modulators. J. Med. Chem. 59, 61–81. doi:10.1021/acs.jmedchem.5b01239
Boubia, B., Poupardin, O., Barth, M., Binet, J., Peralba, P., Mounier, L., et al. (2018). Design, synthesis, and evaluation of a novel series of indole sulfonamide peroxisome proliferator activated receptor (PPAR) α/γ/δ triple activators: discovery of lanifibranor, a new antifibrotic clinical candidate. J. Med. Chem. 61, 2246–2265. doi:10.1021/acs.jmedchem.7b01285
Bredael, K., Geurs, S., Clarisse, D., De Bosscher, K., and D’hooghe, M. (2022). Carboxylic acid bioisosteres in medicinal chemistry: synthesis and properties. J. Chem. Chem. Eng. 2022, 1–21. doi:10.1155/2022/2164558
Casimiro-Garcia, A., Filzen, G. F., Flynn, D., Bigge, C. F., Chen, J., Davis, J. A., et al. (2011). Discovery of a series of imidazo[4,5-b]pyridines with dual activity at angiotensin II type 1 receptor and peroxisome proliferator-activated receptor-γ. J. Med. Chem. 54, 4219–4233. doi:10.1021/jm200409s
Casimiro-Garcia, A., Heemstra, R. J., Bigge, C. F., Chen, J., Ciske, F. A., Davis, J. A., et al. (2013). Design, synthesis, and evaluation of imidazo[4,5-c]pyridin-4-one derivatives with dual activity at angiotensin II type 1 receptor and peroxisome proliferator-activated receptor-γ. Bioorg. Med. Chem. Lett. 23, 767–772. doi:10.1016/j.bmcl.2012.11.088
Chadha, N., and Silakari, O. (2017). Identification of low micromolar dual inhibitors for aldose reductase (ALR2) and poly (ADP-ribose) polymerase (PARP-1) using structure based design approach. Bioorg. Med. Chem. Lett. 27, 2324–2330. doi:10.1016/j.bmcl.2017.04.038
Chemaxon (2023). Chemaxon. Available at: https://www.chemaxon.com (Accessed May 22, 2023).
Chen, X., Chen, Z., Xu, D., Lyu, Y., Li, Y., Li, S., et al. (2021). De novo design of g protein-coupled receptor 40 peptide agonists for type 2 diabetes mellitus based on artificial intelligence and site-directed mutagenesis. Front. Bioeng. Biotechnol. 9, 694100. doi:10.3389/fbioe.2021.694100
Choung, W., Jung, H. J., Nam, E. H., Yang, D., Yoo, B., Choi, H., et al. (2018). Discovery of the bifunctional modulator of angiotensin II type 1 receptor (AT1R) and PPARγ derived from the AT1R antagonist, Fimasartan. Fimasartan. Bioorg. Med. Chem. Lett. 28, 3155–3160. doi:10.1016/j.bmcl.2018.08.036
Colín-Lozano, B., Estrada-Soto, S., Chávez-Silva, F., Gutiérrez-Hernández, A., Cerón-Romero, L., Giacoman-Martínez, A., et al. (2018). Design, synthesis and in combo antidiabetic bioevaluation of multitarget phenylpropanoic acids. Molecules 23, 340. doi:10.3390/molecules23020340
Darwish, K. M., Salama, I., Mostafa, S., Gomaa, M. S., Khafagy, E.-S., and Helal, M. A. (2018). Synthesis, biological evaluation, and molecular docking investigation of benzhydrol- and indole-based dual PPAR-γ/FFAR1 agonists. Bioorg. Med. Chem. Lett. 28, 1595–1602. doi:10.1016/j.bmcl.2018.03.051
Daylight theory: SMIRKS - A reaction transform language (2023). Daylight theory: SMIRKS - a reaction transform language. Available at: https://www.daylight.com/dayhtml/doc/theory/theory.smirks.html (Accessed July 1, 2023).
Dobrică, E.-C., Găman, M.-A., Cozma, M.-A., Bratu, O. G., Pantea Stoian, A., and Diaconu, C. C. (2019). Polypharmacy in type 2 diabetes mellitus: insights from an internal medicine department. Medicina 55, 436. doi:10.3390/medicina55080436
Domínguez-Mendoza, E. A., Galván-Ciprés, Y., Martínez-Miranda, J., Miranda-González, C., Colín-Lozano, B., Hernández-Núñez, E., et al. (2021). Design, synthesis, and in silico multitarget pharmacological simulations of acid bioisosteres with a validated in vivo antihyperglycemic effect. Molecules 26, 799. doi:10.3390/molecules26040799
Dubois, M. A. J., Croft, R. A., Ding, Y., Choi, C., Owen, D. R., Bull, J. A., et al. (2021). Investigating 3,3-diaryloxetanes as potential bioisosteres through matched molecular pair analysis. RSC Med. Chem. 12, 2045–2052. doi:10.1039/d1md00248a
Dulsat, J., López-Nieto, B., Estrada-Tejedor, R., and Borrell, J. I. (2023). Evaluation of free online ADMET tools for academic or small biotech environments. Molecules 28, 776. doi:10.3390/molecules28020776
Ertl, P., and Schuffenhauer, A. (2009). Estimation of synthetic accessibility score of drug-like molecules based on molecular complexity and fragment contributions. J. Cheminform. 1, 8. doi:10.1186/1758-2946-1-8
Estrada-Soto, S., Cerón-Romero, L., Navarrete-Vázquez, G., Rosales-Ortega, E., Gómez-Zamudio, J., Cruz, M., et al. (2022). PPARα/γ, adiponectin, and GLUT4 overexpression induced by moronic acid methyl ester influenced glucose and triglyceride levels of experimental diabetic mice. Can. J. Physiol. Pharmacol. 100, 295–305. doi:10.1139/cjpp-2021-0526
Fang, Y., Zhang, S., Wu, W., Liu, Y., Yang, J., Li, Y., et al. (2020). Design and synthesis of tetrahydropyridopyrimidine derivatives as dual GPR119 and DPP-4 modulators. Bioorg. Chem. 94, 103390. doi:10.1016/j.bioorg.2019.103390
Galicia-Garcia, U., Benito-Vicente, A., Jebari, S., Larrea-Sebal, A., Siddiqi, H., Uribe, K. B., et al. (2020). Pathophysiology of type 2 diabetes mellitus. Int. J. Mol. Sci. 21, 6275. doi:10.3390/ijms21176275
Gilardi, F., Giudici, M., Mitro, N., Maschi, O., Guerrini, U., Rando, G., et al. (2014). LT175 is a novel PPARα/γ ligand with potent insulin-sensitizing effects and reduced adipogenic properties. J. Biol. Chem. 289, 6908–6920. doi:10.1074/jbc.M113.506394
Gutierréz-Hernández, A., Galván-Ciprés, Y., Domínguez-Mendoza, E. A., Aguirre-Vidal, Y., Estrada-Soto, S., Almanza-Pérez, J. C., et al. (2019). Design, synthesis, antihyperglycemic studies, and docking simulations of benzimidazole-thiazolidinedione hybrids. J. Chem. Chem. Eng. 2019, 1–8. doi:10.1155/2019/1650145
Hanf, R., Millatt, L. J., Cariou, B., Noel, B., Rigou, G., Delataille, P., et al. (2014). The dual peroxisome proliferator-activated receptor alpha/delta agonist GFT505 exerts anti-diabetic effects in db/db mice without peroxisome proliferator-activated receptor gamma-associated adverse cardiac effects. Diab. Vasc. Dis. Res. 11, 440–447. doi:10.1177/1479164114548027
He, B. K., Ning, Z. Q., Li, Z. B., Shan, S., Pan, D. S., Ko, B. C. B., et al. (2012). In vitro and in vivo characterizations of chiglitazar, a newly identified PPAR pan-agonist. PPAR Res. 2012, 546548. doi:10.1155/2012/546548
He, Y.-L., Haynes, W., Meyers, C. D., Amer, A., Zhang, Y., Mahling, P., et al. (2019). The effects of licogliflozin, a dual SGLT1/2 inhibitor, on body weight in obese patients with or without diabetes. Diabetes Obes. Metab. 21, 1311–1321. doi:10.1111/dom.13654
Herrera-Rueda, M. Á., Tlahuext, H., Paoli, P., Giacoman-Martínez, A., Almanza-Pérez, J. C., Pérez-Sánchez, H., et al. (2018). Design, synthesis, in vitro, in vivo and in silico pharmacological characterization of antidiabetic N-Boc-l-tyrosine-based compounds. Biomed. Pharmacother. 108, 670–678. doi:10.1016/j.biopha.2018.09.074
Hidalgo-Figueroa, S., Ramírez-Espinosa, J. J., Estrada-Soto, S., Almanza-Pérez, J. C., Román-Ramos, R., Alarcón-Aguilar, F. J., et al. (2013). Discovery of thiazolidine-2,4-dione/biphenylcarbonitrile hybrid as dual PPAR α/γ modulator with antidiabetic effect: in vitro, in silico and in vivo approaches. Chem. Biol. Drug Des. 81, 474–483. doi:10.1111/cbdd.12102
Hidalgo-Figueroa, S., Rodríguez-Luévano, A., Almanza-Pérez, J. C., Giacoman-Martínez, A., Ortiz-Andrade, R., León-Rivera, I., et al. (2021). Synthesis, molecular docking, dynamic simulation and pharmacological characterization of potent multifunctional agent (dual GPR40-PPARγ agonist) for the treatment of experimental type 2 diabetes. Eur. J. Pharmacol. 907, 174244. doi:10.1016/j.ejphar.2021.174244
Holstein, A., Beil, W., and Kovacs, P. (2012). CYP2C metabolism of oral antidiabetic drugs--impact on pharmacokinetics, drug interactions and pharmacogenetic aspects. Expert Opin. Drug Metab. Toxicol. 8, 1549–1563. doi:10.1517/17425255.2012.722619
Hu, L., Zhou, Z., Deng, L., Ren, Q., Cai, Z., Wang, B., et al. (2020). HWL-088, a new and highly effective FFA1/PPARδ dual agonist, attenuates nonalcoholic steatohepatitis by regulating lipid metabolism, inflammation and fibrosis. J. Pharm. Pharmacol. 72, 1564–1573. doi:10.1111/jphp.13342
Huan, Y., Jiang, Q., Li, G., Bai, G., Zhou, T., Liu, S., et al. (2017). The dual DPP4 inhibitor and GPR119 agonist HBK001 regulates glycemic control and beta cell function ex and in vivo. Sci. Rep. 7, 4351. doi:10.1038/s41598-017-04633-5
Huang, H., Meegalla, S. K., Lanter, J. C., Winters, M. P., Zhao, S., Littrell, J., et al. (2019). Discovery of a GPR40 superagonist: the impact of aryl propionic acid α-fluorination. ACS Med. Chem. Lett. 10, 16–21. doi:10.1021/acsmedchemlett.8b00444
Hye Khan, M. A., Kolb, L., Skibba, M., Hartmann, M., Blöcher, R., Proschak, E., et al. (2018). A novel dual PPAR-γ agonist/sEH inhibitor treats diabetic complications in a rat model of type 2 diabetes. Diabetologia 61, 2235–2246. doi:10.1007/s00125-018-4685-0
Ibrahim, M. K., Eissa, I. H., Abdallah, A. E., Metwaly, A. M., Radwan, M. M., and ElSohly, M. A. (2017a). Design, synthesis, molecular modeling and anti-hyperglycemic evaluation of novel quinoxaline derivatives as potential PPARγ and SUR agonists. Bioorg. Med. Chem. 25, 1496–1513. doi:10.1016/j.bmc.2017.01.015
Ibrahim, M. K., Eissa, I. H., Alesawy, M. S., Metwaly, A. M., Radwan, M. M., and ElSohly, M. A. (2017b). Design, synthesis, molecular modeling and anti-hyperglycemic evaluation of quinazolin-4(3H)-one derivatives as potential PPARγ and SUR agonists. Bioorg. Med. Chem. 25, 4723–4744. doi:10.1016/j.bmc.2017.07.015
Kamata, S., Honda, A., and Ishii, I. (2023). Current clinical trial status and future prospects of PPAR-targeted drugs for treating nonalcoholic fatty liver disease. Biomolecules 13, 1264. doi:10.3390/biom13081264
Korn, M., Ehrt, C., Ruggiu, F., Gastreich, M., and Rarey, M. (2023). Navigating large chemical spaces in early-phase drug discovery. Curr. Opin. Struct. Biol. 80, 102578. doi:10.1016/j.sbi.2023.102578
Krishnan, N., Krishnan, K., Connors, C. R., Choy, M. S., Page, R., Peti, W., et al. (2015). PTP1B inhibition suggests a therapeutic strategy for Rett syndrome. J. Clin. Invest. 125, 3163–3177. doi:10.1172/JCI80323
Kumar, M., Kumar Singh, P., Choudhary, S., and Silakari, O. (2022). Hydantoin based dual inhibitors of ALR2 and PARP-1: design, synthesis, in-vitro and in-vivo evaluation. Bioorg. Chem. 129, 106108. doi:10.1016/j.bioorg.2022.106108
Kumari, S., Carmona, A. V., Tiwari, A. K., and Trippier, P. C. (2020). Amide bond bioisosteres: strategies, synthesis, and successes. J. Med. Chem. 63, 12290–12358. doi:10.1021/acs.jmedchem.0c00530
Kuo, G.-H., Gaul, M. D., Liang, Y., Xu, J. Z., Du, F., Hornby, P., et al. (2018). Synthesis and biological evaluation of benzocyclobutane-C-glycosides as potent and orally active SGLT1/SGLT2 dual inhibitors. Bioorg. Med. Chem. Lett. 28, 1182–1187. doi:10.1016/j.bmcl.2018.02.057
Lamotte, Y., Faucher, N., Sançon, J., Pineau, O., Sautet, S., Fouchet, M.-H., et al. (2014). Discovery of novel indazole derivatives as dual angiotensin II antagonists and partial PPARγ agonists. Bioorg. Med. Chem. Lett. 24, 1098–1103. doi:10.1016/j.bmcl.2014.01.004
Landrum, G. RDKit. Available at: http://rdkit.org/ (Accessed July 11, 2023).
Lapuerta, P., Zambrowicz, B., Strumph, P., and Sands, A. (2015). Development of sotagliflozin, a dual sodium-dependent glucose transporter 1/2 inhibitor. Diab. Vasc. Dis. Res. 12, 101–110. doi:10.1177/1479164114563304
Lei, L., Liu, Q., Liu, S., Huan, Y., Sun, S., Chen, Z., et al. (2015). Antidiabetic potential of a novel dual-target activator of glucokinase and peroxisome proliferator activated receptor-γ. Metabolism 64, 1250–1261. doi:10.1016/j.metabol.2015.06.014
Li, Y., Tian, K., Qin, A., Zhang, L., Huo, L., Lei, L., et al. (2014). Discovery of novel urea derivatives as dual-target hypoglycemic agents that activate glucokinase and PPARγ. Eur. J. Med. Chem. 76, 182–192. doi:10.1016/j.ejmech.2014.02.024
Li, Z., Zhou, Z., Deng, F., Li, Y., Zhang, D., and Zhang, L. (2018). Design, synthesis, and biological evaluation of novel pan agonists of FFA1, PPARγ and PPARδ. Eur. J. Med. Chem. 159, 267–276. doi:10.1016/j.ejmech.2018.09.071
Li, Z., Chen, Y., Zhou, Z., Deng, L., Xu, Y., Hu, L., et al. (2019a). Discovery of first-in-class thiazole-based dual FFA1/PPARδ agonists as potential anti-diabetic agents. Eur. J. Med. Chem. 164, 352–365. doi:10.1016/j.ejmech.2018.12.069
Li, Z., Hu, L., Wang, X., Zhou, Z., Deng, L., Xu, Y., et al. (2019b). Design, synthesis, and biological evaluation of novel dual FFA1 (GPR40)/PPARδ agonists as potential anti-diabetic agents. Bioorg. Chem. 92, 103254. doi:10.1016/j.bioorg.2019.103254
Li, Z., Zhou, Z., Hu, L., Deng, L., Ren, Q., and Zhang, L. (2020). ZLY032, the first-in-class dual FFA1/PPARδ agonist, improves glucolipid metabolism and alleviates hepatic fibrosis. Pharmacol. Res. 159, 105035. doi:10.1016/j.phrs.2020.105035
Li, Z., Ren, Q., Zhou, Z., Cai, Z., Wang, B., Han, J., et al. (2021). Discovery of the first-in-class dual PPARδ/γ partial agonist for the treatment of metabolic syndrome. Eur. J. Med. Chem. 225, 113807. doi:10.1016/j.ejmech.2021.113807
Lillich, F. F., Imig, J. D., and Proschak, E. (2020). Multi-target approaches in metabolic syndrome. Front. Pharmacol. 11, 554961. doi:10.3389/fphar.2020.554961
Liu, K., Zhao, X., Qi, X., Hou, D.-L., Li, H.-B., Gu, Y.-H., et al. (2021). Design, synthesis, and biological evaluation of a novel dual peroxisome proliferator-activated receptor alpha/delta agonist for the treatment of diabetic kidney disease through anti-inflammatory mechanisms. Eur. J. Med. Chem. 218, 113388. doi:10.1016/j.ejmech.2021.113388
Maccari, R., and Ottanà, R. (2015). Targeting aldose reductase for the treatment of diabetes complications and inflammatory diseases: new insights and future directions. J. Med. Chem. 58, 2047–2067. doi:10.1021/jm500907a
Madariaga-Mazón, A., Naveja, J. J., Medina-Franco, J. L., Noriega-Colima, K. O., and Martinez-Mayorga, K. (2021). DiaNat-DB: a molecular database of antidiabetic compounds from medicinal plants. RSC Adv. 11, 5172–5178. doi:10.1039/d0ra10453a
Mahindroo, N., Huang, C.-F., Peng, Y.-H., Wang, C.-C., Liao, C.-C., Lien, T.-W., et al. (2005). Novel indole-based peroxisome proliferator-activated receptor agonists: design, SAR, structural biology, and biological activities. J. Med. Chem. 48, 8194–8208. doi:10.1021/jm0506930
Maier, J. A., Martinez, C., Kasavajhala, K., Wickstrom, L., Hauser, K. E., and Simmerling, C. (2015). ff14SB: improving the accuracy of protein side chain and backbone parameters from ff99SB. J. Chem. Theory Comput. 11, 3696–3713. doi:10.1021/acs.jctc.5b00255
Makhoba, X. H., Viegas, C., Mosa, R. A., Viegas, F. P. D., and Pooe, O. J. (2020). Potential impact of the multi-target drug approach in the treatment of some complex diseases. Drug Des. devel. Ther. 14, 3235–3249. doi:10.2147/DDDT.S257494
Medina-Franco, J. L., Chávez-Hernández, A. L., López-López, E., and Saldívar-González, F. I. (2022). Chemical multiverse: an expanded view of chemical space. Mol. Inf. 41, e2200116. doi:10.1002/minf.202200116
MolVS Molecule validation and standardization — MolVS 0.1.1 documentation. Available at: https://molvs.readthedocs.io/en/latest/ (Accessed July 11, 2023).
Otake, K., Azukizawa, S., Takahashi, K., Fukui, M., Shibabayashi, M., Kamemoto, H., et al. (2011). 2-Acyl-tetrahydroisoquinoline-3-carboxylic acids: lead compounds with triple actions, peroxisome proliferator-activated receptor α/γ agonist and protein-tyrosine phosphatase 1B inhibitory activities. Chem. Pharm. Bull. 59, 876–879. doi:10.1248/cpb.59.876
Otake, K., Azukizawa, S., Fukui, M., Kunishiro, K., Kamemoto, H., Kanda, M., et al. (2012). Novel (S)-1,2,3,4-tetrahydroisoquinoline-3-carboxylic acids: peroxisome proliferator-activated receptor γ selective agonists with protein-tyrosine phosphatase 1B inhibition. Bioorg. Med. Chem. 20, 1060–1075. doi:10.1016/j.bmc.2011.11.035
Otake, K., Azukizawa, S., Takeda, S., Fukui, M., Kawahara, A., Kitao, T., et al. (2015). Novel 2,7-substituted (S)-1,2,3,4-tetrahydroisoquinoline-3-carboxylic acids: peroxisome proliferator-activated receptor γ partial agonists with protein-tyrosine phosphatase 1B inhibition. Chem. Pharm. Bull. 63, 998–1014. doi:10.1248/cpb.c15-00508
OTAVAchemicals Ltd - synthetic organic compounds for research and drug discovery. Available at: https://www.otavachemicals.com (Accessed February 23, 2023).
Packer, M. (2023). Dual SGLT1 and SGLT2 inhibitor sotagliflozin achieves FDA approval: landmark or landmine? Nat. Cardiovasc. Res. 2, 705–707. doi:10.1038/s44161-023-00306-x
Park, M. H., Kim, D. H., Kim, M. J., Lee, E. K., An, H. J., Jeong, J. W., et al. (2016). Effects of MHY908, a new synthetic PPARα/γ dual agonist, on inflammatory responses and Insulin resistance in aged rats. J. Gerontol. A Biol. Sci. Med. Sci. 71, 300–309. doi:10.1093/gerona/glv043
Piemontese, L., Fracchiolla, G., Carrieri, A., Parente, M., Laghezza, A., Carbonara, G., et al. (2015). Design, synthesis and biological evaluation of a class of bioisosteric oximes of the novel dual peroxisome proliferator-activated receptor α/γ ligand LT175. Eur. J. Med. Chem. 90, 583–594. doi:10.1016/j.ejmech.2014.11.044
Press, N. J., Taylor, R. J., Fullerton, J. D., Tranter, P., McCarthy, C., Keller, T. H., et al. (2012). Solubility-driven optimization of phosphodiesterase-4 inhibitors leading to a clinical candidate. J. Med. Chem. 55, 7472–7479. doi:10.1021/jm300459a
Press, N. J., Taylor, R. J., Fullerton, J., Tranter, P., McCarthy, C., Keller, T. H., et al. (2015). Discovery and optimization of 4-(8-(3-Fluorophenyl)-1,7-naphthyridin-6-yl)transcyclohexanecarboxylic acid, an improved PDE4 inhibitor for the treatment of chronic obstructive pulmonary disease (COPD). J. Med. Chem. 58, 6747–6752. doi:10.1021/acs.jmedchem.5b00902
Probst, D., and Reymond, J.-L. (2020). Visualization of very large high-dimensional data sets as minimum spanning trees. J. Cheminform. 12, 12. doi:10.1186/s13321-020-0416-x
Rarey, M., Nicklaus, M. C., and Warr, W. (2022). Special issue on reaction informatics and chemical space. J. Chem. Inf. Model. 62, 2009–2010. doi:10.1021/acs.jcim.2c00390
Ren, Q., Deng, L., Zhou, Z., Wang, X., Hu, L., Xie, R., et al. (2020). Design, synthesis, and biological evaluation of novel dual PPARα/δ agonists for the treatment of T2DM. Bioorg. Chem. 101, 103963. doi:10.1016/j.bioorg.2020.103963
Ruddigkeit, L., van Deursen, R., Blum, L. C., and Reymond, J.-L. (2012). Enumeration of 166 billion organic small molecules in the chemical universe database GDB-17. J. Chem. Inf. Model. 52, 2864–2875. doi:10.1021/ci300415d
Ruiz, F. X., Cousido-Siah, A., Porté, S., Domínguez, M., Crespo, I., Rechlin, C., et al. (2015). Structural determinants of the selectivity of 3-benzyluracil-1-acetic acids toward human enzymes aldose reductase and AKR1B10. ChemMedChem 10, 1989–2003. doi:10.1002/cmdc.201500393
Schattenberg, J. M., Pares, A., Kowdley, K. V., Heneghan, M. A., Caldwell, S., Pratt, D., et al. (2021). A randomized placebo-controlled trial of elafibranor in patients with primary biliary cholangitis and incomplete response to UDCA. J. Hepatol. 74, 1344–1354. doi:10.1016/j.jhep.2021.01.013
Seo, S. M., Ihm, S. H., Yi, J.-E., Jeong, S. H., and Kim, B.-S. (2022). Comparative efficacy and safety of fimasartan in patients with hypertension: a network meta-analysis of randomized controlled trials. J. Clin. Hypertens. 24, 971–983. doi:10.1111/jch.14536
Shah, N., Abdalla, M. A., Deshmukh, H., and Sathyapalan, T. (2021). Therapeutics for type-2 diabetes mellitus: a glance at the recent inclusions and novel agents under development for use in clinical practice. Ther. Adv. Endocrinol. Metab. 12, 20420188211042145. doi:10.1177/20420188211042145
Skoraczyński, G., Kitlas, M., Miasojedow, B., and Gambin, A. (2023). Critical assessment of synthetic accessibility scores in computer-assisted synthesis planning. J. Cheminform. 15, 6–9. doi:10.1186/s13321-023-00678-z
Sodano, T. M., Combee, L. A., and Stephenson, C. R. J. (2020). Recent advances and outlook for the isosteric replacement of anilines. ACS Med. Chem. Lett. 11, 1785–1788. doi:10.1021/acsmedchemlett.9b00687
Song, H.-P., Tian, K., Lei, L., Shen, Z.-F., Liu, S.-X., Zhang, L.-J., et al. (2011). Novel N-(pyrimidin-4-yl)thiazol-2-amine derivatives as dual-action hypoglycemic agents that activate GK and PPARγ. Yao Xue Xue Bao 1, 166–171. doi:10.1016/j.apsb.2011.07.002
Subbaiah, M. A. M., and Meanwell, N. A. (2021). Bioisosteres of the phenyl ring: recent strategic applications in lead optimization and drug design. J. Med. Chem. 64, 14046–14128. doi:10.1021/acs.jmedchem.1c01215
Tassopoulou, V.-P., Tzara, A., and Kourounakis, A. P. (2022). Design of improved antidiabetic drugs: a journey from single to multitarget Agents. ChemMedChem 17, e202200320. doi:10.1002/cmdc.202200320
Veith, H., Southall, N., Huang, R., James, T., Fayne, D., Artemenko, N., et al. (2009). Comprehensive characterization of cytochrome P450 isozyme selectivity across chemical libraries. Nat. Biotechnol. 27, 1050–1055. doi:10.1038/nbt.1581
Vilar, S., Cozza, G., and Moro, S. (2008). Medicinal chemistry and the molecular operating environment (MOE): application of QSAR and molecular docking to drug discovery. Curr. Top. Med. Chem. 8, 1555–1572. doi:10.2174/156802608786786624
Walters, W. P. (2019). Virtual chemical libraries. J. Med. Chem. 62, 1116–1124. doi:10.1021/acs.jmedchem.8b01048
Weininger, D. (1988). SMILES, a chemical language and information system. 1. Introduction to methodology and encoding rules. J. Chem. Inf. Comput. Sci. 28, 31–36. doi:10.1021/ci00057a005
Wishart, D. S., Feunang, Y. D., Guo, A. C., Lo, E. J., Marcu, A., Grant, J. R., et al. (2018). DrugBank 5.0: a major update to the DrugBank database for 2018. Nucleic Acids Res. 46, D1074–D1082. doi:10.1093/nar/gkx1037
Xiong, G., Wu, Z., Yi, J., Fu, L., Yang, Z., Hsieh, C., et al. (2021). ADMETlab 2.0: an integrated online platform for accurate and comprehensive predictions of ADMET properties. Nucleic Acids Res. 49, W5–W14. doi:10.1093/nar/gkab255
Xu, G., Gaul, M. D., Kuo, G.-H., Du, F., Xu, J. Z., Wallace, N., et al. (2018). Design, synthesis and biological evaluation of (2S,3R,4R,5S,6R)-5-fluoro-6-(hydroxymethyl)-2-aryltetrahydro-2H-pyran-3,4-diols as potent and orally active SGLT dual inhibitors. Bioorg. Med. Chem. Lett. 28, 3446–3453. doi:10.1016/j.bmcl.2018.09.025
Xu, G., Du, F., Kuo, G.-H., Xu, J. Z., Liang, Y., Demarest, K., et al. (2020). 5,5-Difluoro- and 5-fluoro-5-methyl-hexose-based c-glucosides as potent and orally bioavailable SGLT1 and SGLT2 dual inhibitors. Bioorg. Med. Chem. Lett. 30, 127387. doi:10.1016/j.bmcl.2020.127387
Yang, T., Li, Z., Chen, Y., Feng, D., Wang, G., Fu, Z., et al. (2021). DrugSpaceX: a large screenable and synthetically tractable database extending drug space. Nucleic Acids Res. 49, D1170–D1178. doi:10.1093/nar/gkaa920
Zhang, C., Crawford, J. J., Landry, M. L., Chen, H., Kenny, J. R., Khojasteh, S. C., et al. (2020). Strategies to mitigate the bioactivation of aryl amines. Chem. Res. Toxicol. 33, 1950–1959. doi:10.1021/acs.chemrestox.0c00138
Keywords: computer-aided drug design, diabetes, target-focused libraries, library enumeration, ligandbased design, multi-target design, transformation rules, reaction informatics
Citation: Saldívar-González FI, Navarrete-Vázquez G and Medina-Franco JL (2023) Design of a multi-target focused library for antidiabetic targets using a comprehensive set of chemical transformation rules. Front. Pharmacol. 14:1276444. doi: 10.3389/fphar.2023.1276444
Received: 12 August 2023; Accepted: 16 October 2023;
Published: 02 November 2023.
Edited by:
Concepción González-Bello, Universidade de Santiago de Compostela, SpainReviewed by:
Elena Cichero, University of Genoa, ItalyPranav Shah, National Center for Advancing Translational Sciences (NIH), United States
Copyright © 2023 Saldívar-González, Navarrete-Vázquez 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: Fernanda I. Saldívar-González, fer.saldivarg@gmail.com; José L. Medina-Franco, medinajl@unam.mx