Next Article in Journal
Experimental and Theoretical Insights into the Intermolecular Interactions in Saturated Systems of Dapsone in Conventional and Deep Eutectic Solvents
Previous Article in Journal
Recycling and Degradation of Polyamides
Previous Article in Special Issue
Structural Modeling of TRPA1 Ion Channel—Determination of the Binding Site for Antagonists
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Discovery of a Novel Chemo-Type for TAAR1 Agonism via Molecular Modeling

1
Department of Pharmacy, Section of Medicinal Chemistry, School of Medical and Pharmaceutical Sciences, University of Genoa, Viale Benedetto XV, 3, 16132 Genoa, Italy
2
Institute of Translational Biomedicine, St. Petersburg State University, 199034 St. Petersburg, Russia
3
St. Petersburg University Hospital, St. Petersburg State University, 199034 St. Petersburg, Russia
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Equally contributing senior Authors.
Molecules 2024, 29(8), 1739; https://doi.org/10.3390/molecules29081739
Submission received: 16 February 2024 / Revised: 27 March 2024 / Accepted: 9 April 2024 / Published: 11 April 2024

Abstract

:
The search for novel effective TAAR1 ligands continues to draw great attention due to the wide range of pharmacological applications related to TAAR1 targeting. Herein, molecular docking studies of known TAAR1 ligands, characterized by an oxazoline core, have been performed in order to identify novel promising chemo-types for the discovery of more active TAAR1 agonists. In particular, the oxazoline-based compound S18616 has been taken as a reference compound for the computational study, leading to the development of quite flat and conformationally locked ligands. The choice of a “Y-shape” conformation was suggested for the design of TAAR1 ligands, interacting with the protein cavity delimited by ASP103 and aromatic residues such as PHE186, PHE195, PHE268, and PHE267. The obtained results allowed us to preliminary in silico screen an in-house series of pyrimidinone-benzimidazoles (1a10a) as a novel scaffold to target TAAR1. Combined ligand-based (LBCM) and structure based (SBCM) computational methods suggested the biological evaluation of compounds 1a10a, leading to the identification of derivatives 1a3a (hTAAR1 EC50 = 526.3–657.4 nM) as promising novel TAAR1 agonists.

1. Introduction

Trace amines (TAs) are a group of endogenous chemical messengers closely related to the biogenic amine neurotransmitters, i.e., dopamine (DA), serotonin (5-HT), and norepinephrine (NE). Among them, β-phenylethylamine (β-PEA), p-tyramine, tryptamine, p-octopamine, and others are found in both invertebrate and vertebrate species [1,2,3]. The Trace Amine-Associated Receptor (TAAR) family belongs to the class A G-protein coupled receptors (GPCRs) and consists of six functional members in humans (TAAR1, 2, 5, 6, 8, and 9), with TAAR1 being the most investigated. TAARs are widely expressed at low levels in the brain and periphery [4]. In the central nervous system (CNS), it was originally believed that the TAAR2–TAAR9 subtypes were expressed primarily in the olfactory system, and indeed a role in olfaction has been identified for most of them [5,6,7,8]. However, recent data indicated they are also expressed in limbic brain areas receiving olfactory input and involved in the regulation of emotional behaviors and adult neurogenesis [9,10,11]. By contrast, TAAR1 has no known role in olfaction, and in rodents is highly expressed in the ventral tegmental area and dorsal raphe nuclei [12,13,14], brain regions involved in dopaminergic and serotonergic signaling, respectively. Recently, an association between mutations in the TAAR1 gene and schizophrenia was reported [15,16], suggesting the implication of TAAR1 in the pathways driving schizophrenia [17].
In this scenario, targeting TAAR1 could lead to novel approaches for the treatment of several disorders [18,19,20] including not only schizophrenia but also depression, attention deficit hyperactivity disorder, addiction [21], and metabolic diseases by means of agonist compounds, and Parkinson’s disease by antagonists [22,23,24,25,26,27].
During recent years, the development of TAAR1 ligands has arisen from the main chemical scaffold of the endogenous TAs, leading to the identification of several analogues often endowed with a more promising mouse TAAR1 (mTAAR1) affinity, rather than towards the human ortholog (hTAAR1) [3,28]. This preliminary information highlighted the species-specificity issues complicating virtual screening campaigns and the rational design processes aimed towards novel hTAAR1 effective ligands [29,30,31,32].
In this context, we recently explored via molecular modeling studies, followed by chemical synthesis and biological evaluation, different series of more selective m/hTAAR1 ligands, obtaining in silico the main features for selective compounds [33,34,35].
Hoffmann-La Roche was an early leader in TAAR1 ligand discovery, relying on high-throughput screening (HTS) and repositioning approaches. Indeed, the discovery of novel TAAR1 agonists exploited libraries of known dopaminergic, serotonergic, and adrenergic drugs [36,37,38,39,40] and a number of promising TAAR1 agonists has been discovered [41,42].
Among them, compound S18616 1 (Figure 1A) was reported in the literature as a potent alpha2-adrenoreceptor (α2-ADR) agonist, being then also evaluated as a TAAR1 agonist [43].
To pursue more selective TAAR1 or α2-ADR ligands, structural simplification of the main S18616 1 tricyclic ring has been afforded, leading to different series of more selective TAAR1 ligands [36,37] such as oxazoline derivatives (237) [38]. As shown in Figure 1B, these compounds can be mainly divided into benzyl-based ones (24), phenyl (hetero)alkyl-containing compounds (521), and phenyl-based derivatives (2237).
The most promising identified TAAR1 agonists exhibit: (i) a flexible main spacer tethering the oxazoline core and the terminal aromatic ring (see 5 as prototype of the alkyl-phenyl-based analogues; Figure 1C), (ii) a folded/branched conformation due to a heteroalkyl spacer in the molecule (see 11 as prototype of the heteroaryl-containing analogues; Figure 1C), (iii) rigid positioning due to a phenyl ring directly connected to the main oxazoline core (see 35 as prototype of the phenyl-based analogues; Figure 1C).
These data provide useful guidelines for the design of novel putative TAAR1 agonists. Molecular modeling approaches, such as ligand- and structure-based computational methods, represent widely exploited complemental tools to guide the discovery of novel bioactive molecules [44,45,46,47,48,49,50]. While a number of homology models of m/hTAAR1 have been already reported by us or described in the literature [51,52,53,54], the human TAAR1 model was recently determined thanks to the AlphaFold protein structure database (AF-Q96RJ0-F1) [55,56]. Notably, AlphaFold is used to predict the structures of almost all of the proteins in the human proteome giving high-confidence predicted structures towards new avenues of investigation from a structural perspective [57]. However, homology modeling strategies and deep learning AlphaFold could be complementary, running in different ways and giving different applicability based on the protein of interest [58].
This information prompted us to explore in silico via molecular docking studies the aforementioned hTAAR1 agonist prototype S18616 1 and the expanded series of oxazoline derivatives 137 (see the chemical structures in Table S1). The results were expected to provide more information on the putative binding mode of flexible or not TAAR1 ligands at the protein binding cavity and to pave the way for the following in silico evaluation of novel TAAR1 agonists. Indeed, the in-house small library of compounds 1a10a (Figure 1D), previously reported as human platelet aggregation inhibitors [59], has been herein studied as putative TAAR1 ligands (see the chemical structures in Table S2). The final aims of this study focused on the discovery of novel chemo-types endowed with TAAR1 agonism ability.
Compounds 1a10a are characterized by a flat and mainly aromatic tricyclic ring (Figure 1D) mimicking the main core of the reference compound S18616 1 (Figure 1A). In particular, the 1a10a pyrimidinone ring is expected to mimic the oxazoline ring of the TAAR1 agonists 137. In addition, most of the compounds 1a10a maintained a key basic moiety in R1, to simulate the compound S18616 1 primary amine group.
Herein, in silico investigation of 1a10a has been performed combining ligand-based computational methods (LBCMs) (FLAP2.2.1. software) [60,61] and structure-based ones (SBCMs), such as molecular docking calculations by means of the AlphaFold hTAAR1 model. Using the results of the aforementioned 137 docking mode as a filtering tool and combining the results of the two computational strategies led to the biological evaluation of ten compounds as putative hTAAR1 ligands.
The subsequent in vitro evaluation confirmed the hTAAR1 agonism ability featured by six out of the ten proposed compounds (1a6a), showing for the 1a derivative (hTAAR1 EC50 = 526.3 nM) the lowest EC50 value among the in silico evaluated analogues (hTAAR1 EC50 = 526.3–1.43 μM). Finally, the whole study has been completed by in silico prediction of descriptors related to absorption, distribution, metabolism, and excretion properties (ADME) of 1a6a, in order to prioritize the most promising newly discovered TAAR1 ligands. The obtained results in this study allowed us to validate the computational protocol towards the identification of a novel scaffold for the design of novel TAAR1 agonists. The biological data in tandem with the reported ADME prediction have supported the proposed scaffold as being worthy of further optimization.

2. Results

During the last years, a number of structure-based studies and ligand-based approaches have been attempted to rationalize the structure–activity relationship (SAR) of different series of TAAR1 agonists, such as imidazole- and imidazoline-based compounds [62], as well as a limited number of oxazolines endowed with hTAAR1 agonism ability [35]. Most of the previously investigated derivatives shared the key structural features required for a TAAR1-targeting activity, as a basic core forming a key salt-bridge with a conserved hTAAR1 aspartic acid (ASP103) and an aromatic moiety forming π–π stacking and van der Waals interactions with the surrounding residues.
Very recently, we pointed out the relevant role played by two aromatic cores endowed with H-bonding features within the TAAR1 agonist, to be folded as a “Y-shape” conformation [63], as it is in the case of other GPCR ligands [64,65].
Based on the hTAAR1 AlphaFold protein model, herein we preliminary explored the effectiveness of structural variations afforded to the main scaffold of the conformationally locked hTAAR1 agonist S18616 1, via molecular docking of the analogues 237 (see the chemical structures in Table S1). Details of the calculated scoring functions are listed in Table S3. Most of the oxazolines 137 featured S values spanning from −4.6935 Kcal/mol, such as 15 (hTAAR1 EC50 = 1540 nM), to −7.1321 Kcal/mol (21; hTAAR1 EC50 = 9 nM), with the reference compound S18616 being endowed with S = −6.8545 Kcal/mol.
As shown in Table 1, a perspective of compound 1 and of compounds 221 [38,43], as flexible oxazoline-based TAAR1 agonists, is reported. The predicted ΔG value of each best-scored protein–agonist complex has been listed, based on the performed molecular docking calculations.
Regarding the phenyl-based oxazoline compounds 2237 [38,43], the predicted ΔG value of each best-scored protein–agonist complex is reported in Table 2.
The obtained results are expected to (i) reveal more information about the protein–agonist interactions in presence of differently flexible ligands as those herein considered, and (ii) point out useful guidelines to identify novel promising chemo-types for the development of novel TAAR1 agonists. Indeed, compounds 1a10a have been evaluated in silico and described as follows, thanks to LBCM (FLAP software) [60,61] in tandem with molecular docking studies, as a structure-based tool. Based on the combined LB- and SBCM results, ten compounds have been assayed as putative TAAR1 agonists. Finally, the final biological evaluation of compounds 1a10a allowed us to confirm the reliability of the computational studies.

2.1. Molecular Docking Studies of Oxazoline-Based TAAR1 Agonists

According to our molecular docking calculations, the compound 1 (S18616, EC50 = 15 nM) docking pose revealed an H-bond between the amino group of the ligand and the side chain of residue ASP103 (Figure 2A). This is probably an interaction key for bioactivity, as was also hypothesized for other TAAR1-targeting molecules [66].
Several other non-polar interactions were found between the aromatic six-membered ring of the molecule and residues PHE186, PHE195, and PHE268 (staggered π-stacking interaction). In addition, the cyclohexyl group was engaged in Van der Waals contacts with VAL184, ILE 104, PHE267, and ILE290. A perspective of the ligand and receptor hydrophobic/polar properties is reported in Figure 2B and Figure 2C,D, respectively.
Herein, compound 1 (S18616) was taken as a reference compound for a three-dimensional SAR analysis involving the other analogues 237, based on compound 1’s high potency as TAAR1 agonist and its conformationally locked binding pose.
The TAAR1 ligands 237 were derived by structural simplification of the previous derivative 1 opening the tricyclic ring and obtaining oxazoline-based phenyl-, benzyl-, phenyl (amino)alkyl-containing derivatives.
The benzyl-based ones (24; hTAAR1 EC50 = 154–2900 nM) have been less explored due to their lower potency trend compared to 1, whereas a number of phenyl (hetero)alkyl-containing compounds (521; EC50 = 9–2260 nM) and phenyl-based compounds (2237; hTAAR1 EC50 = 11–10,000 nM) exhibited promising EC50 values.
Within the benzyl-based compound series, the opening of the central cycle of the reference 1 in tandem or not with the elimination of the chlorine atom on the phenyl ring was disadvantageous to the TAAR1 agonism ability, as featured by 34 (hTAAR1 EC50 = 330–2900 nM) and 2 (hTAAR1 EC50 = 154 nM), respectively. In addition, modifying the S enantiomer to the R one impaired the compound potency values, as confirmed by the two enantiomer compounds 3 (S) and 4 (R).
The introduction of a longer spacer tethering the main amino-oxazoline ring and the terminal phenyl group led to the effective phenylethyl- and propylethyl-containing analogues 5 (hTAAR1 EC50 = 18 nM), 6 (hTAAR1 EC50 = 27 nM), and 19 (hTAAR1 EC50 = 27 nM), respectively. This modification was favorable for activity, leading to a decrease in EC50 values of around one order of magnitude with respect to compounds 2 and 3, being quite comparable with the prototype 1 (hTAAR1 EC50 = 15 nM). Moreover, it should be noticed that compounds 5 and 6 differ from each other by a chlorine atom on the phenyl ring (compound 5), which seems to be further advantageous for activity.
As shown in Figure 3A, the best-scored docking pose of compound 6 (hTAAR1 EC50 = 27 nM) highlights a folded ligand conformation that allowed the compound to display the pivotal H-bond with the key residue ASP103, in tandem with additional π-π stacking with the surrounding aromatic residues. Indeed, a favorable interaction between the phenyl ring and PHE268 (T-shaped π-interaction) was highlighted, together with the insertion of the phenyl in a sub-pocket defined by SER198, SER108, SER 197, and PHE268. On the other hand, no hydrophobic contacts involving PHE186 were present.
The introduction of additional alkyl groups at the stereo-center or involving the ethyl-chain spacer was not advantageous, as shown by compounds 7 (EC50 = 330 nM) and 1517 (hTAAR1 EC50 = 730–2260 nM), except for 18 (hTAAR1 EC50 = 18 nM). The lower activity of oxazoline 7 compared to the reference compound 1 was probably due to an unfavorable hindrance between the introduced group and the surrounding residues (VAL184, ILE104, and the key residue ASP103), as depicted in Figure 3B. This event moved the compound quite far from the ASP103 residue, making the corresponding H-bond weaker than that featured by the analogue 6.
Conversely, compound 18 proved to better mimic the docking positioning observed for 1, moving the primary amino group near ASP103 (see Figure 3C).
The choice of the (S)-ethyl substituent branching the main ethyl spacer allowed simulation of the required hydrophobic features of the S18616 1 cyclohexyl group, while the terminal phenyl ring was projected towards the aromatic residue PHE268. Accordingly, the corresponding (R)-ethyl substituted analogue 17 (hTAAR1 EC50 = 2260 nM) was differently folded (Figure 3D), exhibiting higher EC50 values than 18 (hTAAR1 EC50 = 18 nM).
As regards the phenyl heteroalkyl-containing compounds (8–14, 2021; hTAAR1 EC50 = 9–580 nM), the introduction of H-bond acceptor moieties is preferred to that of H-bond donor ones, rendering the presence of tertiary amino groups and phenoxy moieties encouraged. Accordingly, compound 9 exhibiting a secondary amino group tethered to the main oxazoline core displayed lower EC50 values (hTAAR1 EC50 = 580 nM) than the congeners 10–11 (hTAAR1 EC50 = 27 nM), featuring a tertiary amino-containing substituent. The corresponding compound 11 docking pose revealed that the overall pattern of interactions with respect to the reference compound S18616 1 was quite well maintained (Figure 4A).
The H-bond with ASP103 is present and involved both the side chain and the oxygen backbone of the residue. It is possible to notice that the chloro-phenyl moiety exhibited a small displacement with respect to compound 1, but with the chlorine atom pointing in the same direction. The presence of the ethyl substituent on the linker amine partially reproduced the hindrance and hydrophobicity of the junction between the aromatic and central cycle of the reference compound S18616.
According to these data, it is possible to hypothesize that a certain hydrophobicity in this area may be of relevance in ligand affinity, leading the ligand to establish better contacts with the surrounding residues (e.g., PHE195, PHE186, VAL184, ILE104, PHE267), or at least to avoid contacts between polar and non-polar groups. Notably, increasing the hindrance of the aliphatic substituent to iPr dimension as shown by 13 (hTAAR1 EC50 = 140 nM) resulted in decreased activity, possibly because of steric clashes with the mentioned residues. The presence of the chlorine atom, although not fundamental, was advantageous for activity [compare 11 (hTAAR1 EC50 = 29 nM) and 12 (hTAAR1 EC50 = 59 nM)], as it can insert in a small volume delimited by residues SER107, SER108, and PHE268.
The presence of the phenoxy moiety combined with the overall three-atom spacer allowed compound 21 (hTAAR1 EC50 = 9 nM) to exhibit the pivotal folded conformation guiding the proper positioning of the ligand primary amino group near ASP103, as described for all the previous effective TAAR1 agonists (Figure 4B). In addition, stabilizing interactions were observed between the phenyl-ring of compound 21 and the surrounding aromatic residues (PHE195, PHE268, PHE267).
The elimination of the linker, leading to a direct bond between the two cyclic systems, resulted in acceptable activity as shown by a consistent number of phenyl-based compounds (2237; hTAAR1 EC50 = 11-10.000 nM), as reported for 2224, 26, 30, 3337 (EC50 = 21–67 nM). The docking pose of compound 23 (hTAAR1 EC50 = 23 nM), taken as representative of these effective TAAR1 agonists, revealed a similar binding mode with respect to the reference compound S18616. As shown in Figure 4C, the phenyl ring is oriented similarly, and the H-bond with ASP103 is maintained. Small differences are highlighted, such as a variation in the chlorine orientation, and a different conformation of the oxazoline ring, lacking contacts with residue TRP264.
Meta- (24; hTAAR1 EC50 = 21 nM) and para- (25, hTAAR1 EC50 = 143 nM) chloro substitution were explored in addition to the previously cited ortho- analogue 23. While the meta-substitution was well tolerated, the para-substitution led to a considerable increase in EC50. Then, para-bromo substitution (27; hTAAR1 EC50 = 150 nM) leads to a comparable activity with the previous compound 25. Notably, when the para-Br substitution is combined with the stereocenter inversion to R configuration, a strong loss in activity occurred (28; hTAAR1 EC50 = 10,000 nM). Interestingly, dichloro substitution at the meta- and para- positions (26; hTAAR1 EC50 = 31 nM) led to a good activity. Compound 32 (hTAAR1 EC50 = 490 nM), with a fluorine in meta- position, exhibited a strong decrease in activity compared to the meta-chloro-containing analogue 24 (hTAAR1 EC50 = 21 nM). However, if this modification is accompanied by the presence of a methyl substituent in ortho- position, the activity is comparable to reference compound 1 (37; hTAAR1 EC50 = 17 nM).
Accordingly, further to the halogen atom, the additional introduction of electron-donor substituents endowed with hydrophobic properties, such as aliphatic chains, led to improved TAAR1 agonists, as reported for 3436 (hTAAR1 EC50 = 11–26 nM). Based on the molecular docking pose of the very potent compound 34 (hTAAR1 EC50 = 11 nM), the oxazoline is shifted towards residue VAL184, establishing good non-polar interactions between this residue and ILE290 (Figure 4D). The chlorine atom is inserted in a sub-pocket delimited by residues THR196, SER198, and PHE195.
The presence of an ethyl- or cyclopropyl- group in place of the methyl in ortho also maintained an excellent activity, as reported by compounds 35 (hTAAR1 EC50 = 26 nM) and 36 (hTAAR1 EC50 = 12 nM). Bulkier substituents such as the diaryl moiety were not tolerated, as described for compound 31 (hTAAR1 EC50 = 2670 nM). Finally, the only insertion of the ortho-methyl substitution led to quite effective EC50 values [see 33 (hTAAR1 EC50 = 67 nM)].

2.2. In Silico Evaluation of Pyrimidinone-Benzimidazole-Based Compounds as TAAR1 Agonists

According to the aforementioned 137 molecular docking studies, H-bonding to ASP103 and featuring π–π stacking with aromatic residues such as PHE186, PHE195, TRP264, and PHE268 was advantageous for TAAR1 agonism. This pharmacophore feature was achieved by the most effective compounds of the 137 series, exhibiting a branched structure endowed with H-bonding and aromatic features. As a result, the TAAR1 agonist “Y-shape” conformation was encouraged. On this basis, we focused on a small library of pyrimidinone-benzimidazole-based compounds 1a10a, previously reported by us as human platelet aggregation inhibitors (see the chemical structures in Table S2) [59], to be evaluated in silico as putative TAAR1 ligands, prior to their biological evaluation. As shown in Figure 5, all compounds are characterized by a flat and mainly aromatic tricyclic ring mimicking the main core of the reference compound S18616 1.
While the pyrimidinone ring A featured comparable electrostatic properties and electron-rich atoms to the corresponding 1 oxazoline ring (ring A), the terminal six-membered ring of the benzimidazole portion guaranteed the proper aromatic moieties displayed by the prototype ring C. On the other hand, the central ring B of the 1a10a series was endowed with additional nitrogen atoms if compared to the S18616 1 ring B. The previously reported oxazoline series 237 pointed out the effectiveness of proper heteroatoms within the main spacer tethering the two terminal A and C rings, therefore supporting this structural variation.
As shown in Figure 5, most of the compounds maintained a key basic moiety in R1, in order to simulate the compound S18616 1 primary amine group, while only three compounds bearing a halogen atom (7a, 10a) or a hydroxyl group (9a) in R1 have been taken into account.
Finally, the steric hindrance displayed by 1a10a should be in good agreement with that of 1, even if the molecule volume and shape are different. Based on the above, we assessed in silico the putative TAAR1 agonism ability of the pyrimidinone-benzimidazole scaffold, using ligand- and structure-based virtual screening methods.
Initially, we proceeded with a LBCM approach, using compound S18616 1 as reference compound. The results of the LBCM were ordered according to the Glob-Prod parameter (see experimental section), being the corresponding score values reported in Table S4. As control compounds, the TAAR1 agonists 237 have been also evaluated via LBCM. As shown in Table S5, the benzyl-based ones (24; EC50 = 154–2900 nM) have been predicted as very promising (Glob-Prod = 0.5961–0.6494), due to the similar volume and dimension compared to S18616; most of the phenyl (hetero)alkyl-containing compounds (6, 7, 9, 1315, 1721; EC50 = 9–2260 nM) were shown as effective, featuring Glob-Prod values spanning from 0.4090 to 0.5047. The phenyl-based compounds (2237; EC50 = 11–10,000 nM) have been predicted as the less comparable oxazolines with respect to the template S18616, being in any case endowed with hTAAR1 agonism ability. In this case, the corresponding Glob-Prod values spanned from 0.3051 to 0.4076. By a perspective of the activity values featured by all the S18616 analogues 137, a Glob-Prod value of 0.3000 has been herein considered as acceptable for the ligand-based evaluation of further hTAAR1 agonists.
As regard the pyrimidinone-benzimidazoles 1a10a, three derivatives (1a3a) were endowed with Glob-Prod values higher than 0.3000, while the analogues 4a6a spanned from 0.2700 to 0.3000. Compounds 7a10a displayed the lowest values of Glob-Prod in this series of compounds (Glob-Prod = 0.1051–0.2692). Then, we proceed with a 3D-analysis of the results via visual comparison of the pyrimidinone-benzimidazole candidates’ molecular interaction fields (MIFs) with respect to those of the template S18616 1, following the score order. The following MIF probes have been taken into account: (i) DRY to evaluate the role played by hydrophobic features, (ii) N1 to assess favorable interaction with the H-bond donor (HBD) probe, (iii) O to assess favorable interaction with the H-bond acceptor (HBA) probe.
As regards the top-scored pyrimidinone benzimidazoles 1a3a, compound 1a (Glob-Prod = 0.3499) displayed the highest Glob-Prod value (see Figure 6A). Good compatibility with the template S18616 1 was observed for the three analyzed interactions (DRY, N1, and O MIFs, shown in yellow, blue, and red, respectively).
The O MIF (red) of 1a is well matching the template one in the area of the piperazine-oxazoline moieties, as the protonated piperazine ring of this compound mimics the amino-group of the amino-oxazoline of S18686. A small common area is also present for N1 MIF (blue MIF), as the oxazoline nitrogen displayed a similar interaction field as the piperazine tertiary amine of the candidate. DRY MIFs (yellow) largely overlap both on the front and back faces of the compounds, highlighting common interaction areas given by the aromatic/cyclic moieties and their hydrophobic substituents.
Then, compound 3a (Glob-Prod = 0.3167) exhibited overall good compatibility with the template, especially for the hydrophobic MIF (see Figure 6B). Again, the O MIF superposition only involved one of the two volumes shown by the template, as only one HBD of the candidate is coherently oriented with the S18686 amine group. N1 MIF overlaps are negligible. Compound 2a (Glob-Prod = 0.3154) exhibited similar behavior, with only one O MIF superimposed on the template one (see Figure 6C). In this case, the N1 MIF of the candidate was partially superimposed on the template one, as the piperazine tertiary amine (2a) mimics the oxazoline basic nitrogen (S18686). The DRY MIFs exhibit good superimposition on the two faces of the compounds.
The following set of best-scored compounds included the analogues 4a6a (Glob-Prod = 0.2700 to 0.3000), which were predicted as the most promising of the last ones 7a10a (Glob-Prod = 0.1051–0.2692).
According to the Glob-Prod ranking, the pyrimidinone-benzimidazole 6a (Glob-Prod = 0.2962) exhibited modest compatibility with the template, especially in terms of HBD and HBA interactions, while hydrophobic fields maintained a certain coherence with the template (see Figure 6D).
Compound 4a (Glob-Prod = 0.2841) displayed a small superposition between the O MIFs, in particular due to a similar orientation of the amino moiety of S18686, and the protonated pyrrolidine of the aforementioned 4a (see Figure 6E). An overlap between N1 MIF was also present: the involved moieties were the oxazoline basic nitrogen and the tertiary amine of the piperazine ring of 4a. However, such overlap derived from different areas of the molecules, and should not be regarded as a signal of ligand similarity. The overlap between DRY MIFs is strongly decreased, possibly due to the altered positioning of the tricyclic plane of 4a. The analogue 5a (Glob-Prod = 0.2838) displayed O MIFs superposition due to consistent positioning of one of the NH of the amino group of S18686 and one of the NH of the protonated piperazine of 5a (see Figure 6F). No N1 MIFs superpositions are highlighted. DRY MIFs exhibited relatively large areas of overlap.
Then, we found a lesser compatibility of the MIFs. Compound 9a (Glob-Prod = 0.2692) exhibited no superposition of O MIFs, but a certain overlap of N1 fields, due to a similar positioning of S18686 oxazoline nitrogen, and the pyridine one 9a (Figure S1A). DRY MIFs overlap quite widely in the front area of the molecules. Compound 10a (Glob-Prod = 0.2463) showed no significant MIF superpositions. In particular, the O MIF deriving from the protonated amine of the candidate partially overlaps with the N1 MIF generated by the oxazoline nitrogen of S18686 (Figure S1B).
Similarly, the two last-scored compounds, 7a (Glob-Prod = 0.1051) and 8a (Glob-Prod = 0.2435) (Figure S1C,D), showed very small or no superpositions of HBD/A or DRY MIFs with the template, respectively, which is the reason for the low positioning of these candidates in the score list. In Figure S2 an overall perspective of the main MIFs shared by the best-scored candidate 1a and the lowest analogue 7a, with respect to the template S18616, was reported.
While ligand-based studies allowed us to compare the main steric and electrostatic features of 1a10a with those of 1, taken as reference agonist, molecular docking studies of the in-house pyrimidinone-benzimidazole series were exploited to assess the ability of this chemo-type to simulate the docking mode of the 137 agonists.
Molecular docking studies of compounds 1a10a have been performed, relying on the AlphaFold protein structure. Details of the molecular docking calculations are shown in the experimental section. The obtained scoring functions for the pyrimidinone-benzimidazole-based compounds are listed in Table S6.
As shown in Table 3, a perspective of the chemical structure of 1a10a [59] and of the predicted ΔG value of each best-scored protein–ligand complex is listed. The corresponding hTAAR1 EC50 values have been added, based on the following biological assays.
The derived docking poses of only the most promising TAAR1 agonists based on the obtained scoring functions, with respect to the control compound S186161, are discussed. As a result, the putative most interesting derivatives were 1a3a, exhibiting S values < −3.0000 Kcal/mol, while 4a6a spanned from −3.1141 to −2.6083 Kcal/mol. Conversely, the analogues 7a10a were predictive as poorly active/inactive exhibiting S values > −1.4568 Kcal/mol. On this basis, the docking modes of only 1a6a are discussed as follows.
According to our calculations, the area of interaction of most of the aforementioned 1a6a is the same of the reference one S18616 1.
The choice of a piperazine group in R1 in tandem with a small hydrophobic chain in R is predicted as the most beneficial substitution involving the pyrimidinone-benzimidazole core, as shown by 1a and 2a, based on their obtained scoring functions.
In the compound 1a the tri-cyclic core established non-polar interactions with residues PHE186, THR194 (methyl), ILE104, PHE268, PHE267, PHE195, VAL184, and ILE290, with the N-methyl substituent inserted within PHE186, THR100, and VAL184 sidechains (Figure 7A). The piperazine NH established an H-bond with ASP103 carboxylic group, while the aliphatic part of the ring contacts ILE290, VAL184, THR100 (methyl), and HIS99.
For compound 2a, the orientation of the core is more similar to the 1a case, with the phenyl ring on the left side, and the cyclo-amide on the right. However, the result is that the compound is upside-down, as the N-aliphatic substituent is pointing downwards, while the carbonyl is pointing upwards (Figure 7B). The core of the compound established non-polar interactions with ILE104, THR194, VAL184, PHE268, PHE267, and PHE195 (offset π-stacking). London forces were present between the iPr moiety of the ligand and residues PHE267, ILE290, and TRP264. The piperazine ring again made one H-bond with ASP103, while the non-polar part of the ring interacted with residues VAL184, THR100, HIS99, ILE290.
The introduction of an aromatic pendant at the N (10) atom led to the derivative 3a and 6a, exhibiting the benzyl and the pyridinil-methyl-groups in R, respectively.
In case of compound 3a, the interaction pattern was coherent with the previous case, contacting PHE195, ILE104, PHE268, TRP264, with the benzyl inserted among residues PHE267 (staggered π-interaction), VAL184, THR271, PHE195 (T-shaped π-stacking), and PHE186. The protonated piperazine ring established one H-bond with ASP103, in addition to non-polar contacts with residues PHE267, ILE290, TYR294, VAL184 (see Figure S3A).
Regarding 6a, the tri-cyclic core orientation slightly differed from the previous compounds (1a and 3a). In this case, in fact, the phenyl group was positioned in the upper part of the binding pocket, within residues PHE186, THR194, PHE195 (staggered π-interaction), THR271, VAL184, and PHE267, while the rest of the core contacted residues PHE268, and PHE199 (see Figure S3B). Interestingly, the chlorine atom was inserted in a volume delimited by residues PHE199, SER108, and SER107, as in the case of some of the amino-oxazoline derivatives previously examined. When a pyridine ring was introduced in place of a piperazine, no H-bonds were possible with residue ASP103. However, this moiety established several other non-polar interactions with residues VAL184, ILE290, ILE104, TRP264, and PHE267 (staggered π-interaction). Interestingly, the pyridine ring was also involved in an intra-molecular T-shaped π-stacking with the phenyl ring of the tri-cyclic system.
Finally, the introduction of heteroalkyl-based chains in N (10) led to the pyrimidinone benzimidazoles 4a and 5a, bearing the pyrrolidine and morpholine substituents in R, respectively.
As regards compound 4a, the piperazine moiety established one H-bond with ASP103 (Figure 7C).
The tri-cyclic core established non-polar interactions with the surrounding residues (ILE290, PHE267, PHE268, PHE195, ILE104). The ethyl-N-cyclopentyl- moiety protrudes towards ASN164, and established non-polar contacts with residues PHE186, ILE159, VAL184, THR100.
The orientation of the 5a core allowed the morpholine ring to be H-bonded with TYR294 (through the oxygen atom) and ASP103 (by means of the protonated tertiary amine) (see Figure 7D). In addition, this moiety also established non-polar interactions with VAL184, PHE267, TRP264, ILE290. The tri-cyclic core interacted with the surrounding residues (ILE104, THR194, PHE186, PHE267, and PHE195 (offset π-interaction with the phenyl ring of the ligand). The piperazine ring was inserted in a hydrophobic sub-pocket, interacting with PHE268, TRP264, ILE111, and PHE260.

2.3. Biological Evaluation of Compounds 1a10a as Novel TAAR1 Agonists

According to the results obtained via LBCM in tandem with those suggested by the previously cited structure-based results, we proceeded with the biological evaluation of the ten explored pyrimidinone-benzimidazoles 1a10a. While compounds 1a6a proved to be the most promising in silico, based on the previously cited computational studies, in vitro screening of the derivatives 7a10a has been also performed to elucidate the SAR of the series.
Thus, all the novel compounds 1a10a were tested at 10 μM to evaluate their potential hTAAR1 agonism ability (Figure 8), by an in vitro screening method based on the bioluminescence resonance energy transfer technique (BRET) (see the experimental section) [63].
The activity of these compounds was determined by means of HEK-293 cells transfected with hTAAR1, or empty vector as control, and a cAMP BRET biosensor. The applied positive control compound for agonism was tyramine hydrochloride (1 μM). The EC50 value of the positive control derivative was 41.6 nM under the experimental conditions.
For the compounds that were confirmed as active, a dose-response experiment was performed using concentrations ranging from 10 nM to 100 μM to derive the related EC50 values. All the compounds displayed an Emax > 85% at m/hTAAR1, and thus they were regarded as full agonists.
Among the active compounds 1a6a, the piperazine-containing derivatives 1a5a (hTAAR1 EC50 = 525–1060 nM) turned out to be more potent than 6a (hTAAR1 EC50 = 1430 nM). In particular, the choice of the piperazine basic ring in the presence of a small hydrophobic chain in N (10) led to the most promising pyrimidinone-benzimidazoles 1a2a (hTAAR1 EC50 = 526–657 nM), endowed with comparable or higher potency values with respect to the 18% of the oxazoline derivatives 137 previously reported. These effective structural variations were in accordance with the reported molecular docking poses revealing beneficial Van der Waals contacts and π-π stacking with the hTAAR1 cavity. Removing the piperazine ring in R1 led to compounds 6a10a, most of them (7a10a) being inactive as TAAR1 agonists.
The introduction of a flexible aromatic pendant in N (10) was well-tolerated, especially in the case of a maintained piperazine ring in R1. Accordingly, compound 3a (hTAAR1 EC50 = 756 nM) displayed lower hTAAR1 EC50 values than 6a (hTAAR1 EC50 = 1430 nM). Finally, including electron-rich atoms within the R substituents could impact the overall positioning of the compound, leading to mildly potent hTAAR1 agonists (see 4a, 5a; hTAAR1 EC50 = 1010–1060 nM).
Based on the obtained biological assays, a perspective of the calculated ADMET properties for the most promising pyrimidinone-benzimidazole derivatives has been performed.
Nowadays, in silico prediction of absorption, distribution, metabolism, excretion, and toxicity (ADMET) properties is thought to be a valuable supporting tool in medicinal chemistry and these techniques are widely exploited in the literature [34,68,69].
Herein, ADMET properties explaining the drug-like profile of compounds 1a-6a and S186116 have been predicted, via the ACD/Lab Percepta platform [70] and SwissADME website [71]. Both of them are well-known tools applied in the prediction of pharmacokinetic properties [72,73,74,75,76,77,78].
Based on Veber’s [79] and Lipinski’s rules [80], the following descriptors have been calculated: (i) logarithmic ratio of the octanol–water partitioning coefficient (cLogP), (ii) molecular weight (MW) of compounds, (iii) H-bonding acceptor number (HBA), and H-bonding donor moieties (HBD), (iv) number of rotatable bonds (nRot_bond), (v) topological polar surface area (TPSA) (see Table 4).
All the herein explored TAAR1 agonists 1a6a and the reference S18616 1 fulfilled most of Lipinski’s rule and Veber’s rule.
Further ADME parameters have been taken into account such as: (i) human intestinal absorption (HIA), estimation of the plasmatic protein binding event (% PPB), volume of distribution (Vd), ligand affinity toward human serum albumin (LogKa HSA) and of putative oral bioavailability, as a percentage (F %) (see Table 4).
Based on the above, among the newly developed TAAR1 agonists 1a6a the most interesting, 1a (EC50 = 526 nM), featuring a small alkyl chain in R, showed quite comparable MW and TPSA values to those of the reference ligand S18161. On the other hand, it was endowed with modest HIA and F% values, if compared to the S18161.
The choice of H-bonding features in R (4a, 5a; EC50 = 1010–1060 nM) led to reduced LogP values and increased the TPSA ones. This kind of variation impaired the HIA levels as well as the predicted F% parameter. Excessive TPSA values are also predicted to compromise blood–brain barrier permeability, making 5a probably inactive in the CNS.
Conversely, the introduction of the benzyl substituent in R, as shown by 3a (EC50 = 756 nM), allowed an increase in the LogP value with respect to the alkyl-based analogues 1a2a, turning in ameliorated HIA and F% values. Indeed, 3a was endowed with quite comparable F values with respect to S18616. Accordingly, the pyridine containing- derivative 6a (EC50 = 1430 nM) also exhibits high HIA and modest F% values, being predicted as able to pass the blood–brain barrier.
In silico evaluation of toxicity properties in terms of cytochrome inhibition and of lethal dose via mouse oral administration have been explored as well as PAINS (Pan-Assay Interference Structures) analysis (see Table S7). The results pointed out no cytochrome inhibition events or PAINS ability. Finally, additional in silico evaluations of the most promising compounds 1a3a have been performed via SwissTarget [71,81,82,83,84], suggesting a prominent role as GPCRs or enzyme-targeting compounds for the three aforementioned pyrimidinone-benzimidazoles 1a3a (see Figure S4). This information is in agreement with the results herein reported regarding their TAAR1 agonism ability and also for the previously described mechanism of action as human platelet aggregation inhibitors via phosphodiesterase inhibition [59].

3. Discussion

The development of TAAR1 ligands represents an intriguing strategy for the treatment of several disorders including schizophrenia, depression, attention deficit hyperactivity disorder, Parkinson’s disease, addiction, and metabolic diseases [22]. As shown in the literature, repositioning approaches involving different series of GPCR-targeting compounds have been explored [23,24,25], such as dopaminergic, serotonergic, and adrenergic drugs, such as the agonist S18616.
Recently, structural variations of the main S18616 tricyclic ring were afforded, leading to different series of more selective TAAR1 ligands featuring the amino-oxazoline main core [36,37,38].
Compounds 137, previously disclosed in the literature and herein investigated in silico, displayed effective substitutions linked to the main five-membered ring, such as oxazoline-based phenyl-, benzyl-, phenyl (amino)alkyl-containing derivatives. Among them, a number of phenyl (hetero)alkyl-containing compounds (521; EC50 = 9–2260 nM) and phenyl-based (2237; EC50 = 11–10,000 nM) compounds showed promising hTAAR1 EC50 values.
Herein, we explored via molecular docking studies the effectiveness of these structural variations in terms of TAAR1 targeting ability, thanks to the protein AlphaFold model. The results pointed out the relevant role played by aromatic cores linked to the oxazoline ring, being also endowed with H-bonding features. The introduction of branched and/or folded groups was promising, being in accordance with the pivotal role played by the TAAR1 agonist “Y-shape” conformation [62,64,65].
In particular, halogen and aliphatic substitutions on the additional phenyl ring exhibited complex effects, and the combination of substituents seems relevant. For the compounds without a linker, p-substitution alone was not very well tolerated, but activity was restored when it was combined with o- or m-substitutions. In the case of fluorine, meta-substitution alone was not tolerated, but further o-methyl-substitution led to excellent activity (compound 37). For compounds bearing a linker, the information on phenyl substitution was not very abundant, being that unsubstituted rings exhibited good activity.
In terms of protein–ligand interactions, the key feature for successful compounds was the unperturbed establishment of the key H-bond between the amine-oxazoline moiety and ASP103, together with at least one hydrophobic moiety. In particular, aromatic rings can establish π-interactions, in addition to London forces, due to the presence of several aromatic residues in the TAAR1 pocket. The presence of an additional steric hindrance in the sub-pocket formed by SER107, SER108, and PHE268 was favorable.
On this basis, following molecular modeling, studies have been developed based on LBCM and SBCM to explore compounds 1a10a [59] as putative TAAR1 agonists (the chemical structures in Table S2). As regards the structure-based results, the flat and mainly aromatic tricyclic ring of all the pyrimidinone-benzimidazole derivatives 1a10a mimics the main core of the reference compound S18616, compounds 1a3a being predicted as the most promising as putative hTAAR1 ligands. In particular, the choice of a piperazine group in R1 in tandem with a small hydrophobic chain in R allowed the compounds to exhibit key contacts with the protein cavity, such as: (i) non-polar interactions with residues PHE186, THR194, ILE104, PHE268, PHE267, PHE195, VAL184, and ILE290, and (ii) an H-bond with the ASP103 carboxylic group.
The results obtained by the concomitant ligand-based studies were in agreement with the following biological assays. Indeed, the best-scored compounds 1a3a (hTAAR1 EC50 = 526.3–756.80 nM) were properly predicted as the most effective TAAR1 agonists within the pyrimidinone-benzimidazole series, exhibiting the highest Glob-Prod score values (Glob-Prod = 0.3154–0.3499). The calculated 1a3a MIFs were in good agreement with those featured by the template S18616, thanks to the corresponding basic moieties and aromatic features. While the analogues 4a6a partially mimic the template MIFs due to the amine groups, the lowest-scored congeners 7a10a poorly simulated the reference MIFs. Accordingly, the following biological evaluation revealed 4a6a (hTAAR1 = 1010–1430 nM) as modest TAAR1 agonists. A perspective of this information pointed out the pyrimidinone-benzimidazole series as a novel chemo-type for the development of novel TAAR1-targeting compounds, exhibiting a new core if compared to the limited number of the so far explored TAAR1 agonists.
Finally, in silico prediction of ADMET properties and PAINS (Pan Assay Interference structures) analysis gave further support to the most promising compounds 1a3a as drug-like derivatives. Among them, compound 3a was endowed with the highest predicted LD50 values and oral bioavailability parameters.

4. Materials and Methods

4.1. Computational Studies

All the studied compounds, including the oxazoline 137 and the pyrimidinone benzimidazoles 1a10a, were manually built by the MOE Builder program included in the MOE software [67] and assigned the prevalent protonation state at pH = 7.4 through the wash function. Ligand parametrization was carried out with the AM1 partial charges as calculation method. Energy minimization step has been performed by the Energy Minimize Program using MMFF94x forcefield of MOE and RMS (root mean square) gradient equal to 0.0001, being the root mean square gradient of the norm of the gradient times the square root of the number of (unfixed) atoms. This approach allowed us to build in silico a single low-energy conformation for each compound [67]. The tridimensional model of the hTAAR1 receptor has been retrieved via the AlphaFold protein model database [56], using the AF-Q96RJ0-F1 code. The model was processed by means of the “Structure Preparation” tool from MOE 2019.01 suite and refined via the “Protonate3D” tool as previously mentioned [61]. This allowed us to assign the most probable protonation state to each residue being the partial charges determined according to the AMBER10:EHT force field, included in the MOE software.
All the molecular docking simulations at the hTAAR1 AlphaFold protein model were performed by means of the DOCK tool included in the MOE software, via a template-based approach using the previously described S18616-TAAR1 complex [61].
This was developed applying the Alpha Triangle method and Affinity ΔG prediction as final scoring function to obtain from 50 to 10 best score conformers. Details of the template-based docking calculations as well as of the scoring functions are shown in our previous papers [85,86].
The LBCM studies were performed with FLAP2.2.1 Ligand-Based module [60,61], using as a template the minimized form of compound S18686 [43], a well-known TAAR1 agonist endowed with high potency and a rigid scaffold. For each of the ten candidates, 25 conformers were generated starting from their optimized geometry (ligand preparation step previously described). Four probes were selected for the analysis, namely the H probe (expressing shape similarity), DRY (hydrophobic probe), O (HBA probe, corresponding to a carbonyl-group), and N1 (HBD probe, corresponding to amide NH). Briefly, the ligand is inserted in a grid with nodes at a defined distance from each other (here 0.75Å). Each probe is positioned in correspondence with the grid nodes, and the interaction energies with the ligand are calculated. The output of this analysis can be visualized as 3D-objects called Molecular Interaction Fields (MIFs), identifying areas of favored interactions between the ligand and the considered probe. For each molecule, several sets of quadruplets representing the most favorable interactions are extracted by the MIFs, enabling an easily-comparable representation. All the quadruplets of a molecule are stored together, representing a pharmacophore fingerprint of the molecule itself. Such quadruplets are used to generate different superimpositions to the template, and the Probe scores and Distance scores are produced. The Probe scores represent the degree of overlap of the MIFs of the test molecule and the template for each probe individually. The Distance score indicates the overall difference of probe score between the ligand and template. The product of the individual Probe scores gives the Glob-Prod global score, ranging from 0 (bad) to 1 (good) [60,87]. In the present case, the Glob-Prod score was considered for compound evaluation, as different types of interactions are expected based on the shape and electrostatic/steric properties of the reference compound 1. Moreover, MIF information was used to visually compare the in-house set of ligands with the known active compound, the template (S18686) [43] In fact, the superposition of the same type of MIF between the target molecule and the template is regarded as sign of ligand similarity, enhancing the probability of giving a similar biological output [60].

4.2. Chemicals

Compounds 1a10a were synthesized as previously described [46]. The synthetic route to compounds 1a5a and 7a, 10a included the reaction of 2-aminobenzimidazole with proper alkyl halides in refluxing acetone in the presence of anhydrous K2CO3/ KOH to obtain the corresponding 1-substituted 2-aminobenzimidazoles. On the contrary, in the case of the pyridine-containing compounds 6a and 9a, 2-aminobenzimidazole was treated with an equimolar amount of EtONa in anhydrous EtOH to give the sodium salt derivative, prior to treatment with the proper alkyl halide.
The 1-substituted 2-aminobenzimidazoles were then condensed with diethyl malonate in the presence of EtONa to afford good yields of crystalline sodium salts of the corresponding 10-substituted 2-hydroxypyrimido [1,2-a]benzimidazol-4(10H)-ones. These compounds were subsequently treated with excess refluxing POCl3 affording the corresponding 2-chloroderivatives (such as 6a, 7a). The following reaction of 2-chloroderivatives with excess piperazine in refluxing ethanol led to the desired 10-substituted 2-(1-piperazinyl) pyrimido [1,2-a]benzimidazol-4(10H)-ones (such as 1a5a).

4.3. Biological Evaluation of Compounds 1a10a as TAAR1 Agonists

Bioluminescence Resonance Energy Transfer (BRET) Measurement

BRET screening has been described in detail elsewhere [63,88,89,90,91,92,93,94]. Briefly, human embryonic kidney 293 cells (HEK293T; ATCC CRL-11268) were transiently co-transfected with plasmids encoding hTAAR1 and a cAMP BRET biosensor using Lipofectamine® 2000 reagent (ThermoFisher, Waltham, MA, USA), and then plated in 96-well white opaque microplates (Corning, Corning, NY, USA) at a density of 50 × 104 cells per well. After 24-h incubation, culture medium was removed and 70 µL of Hank’s balanced salt solution containing calcium and magnesium was added to each well, followed by the addition of 10 µL 200 mM 3-isobutyl-1-methylxanthine solution (Sigma, St. Louis, MO, USA), a phosphodiesterase inhibitor, and 10 µL 50 µM coelenterazine-h solution (Promega, Madison, WI, USA), as luciferase substrate. All tested compounds were dissolved in DMSO to yield 10 mM stock solutions. After 10 min incubation, either 10 µL of vehicle or 10× of the concentrated solution of compound to be tested (10 µM final concentration) was added. After an additional 5-min incubation, luminescence readings were collected using a Mithras LB943 multimodal plate reader (Berthold Technologies, Wildbad, Germany). The BRET signal was determined by calculating the ratio of the light emitted at 505 to 555 nm to the light emitted at 465 to 505 nm. To confirm the specificity of positive responses, parallel screening on cells transfected with empty core vector was carried out. For active compounds, separate dose–response experiments were performed in order to calculate the EC50 values. Curves were fitted by applying non-linear regression models in GraphPad Prism 6 (GraphPad Software, San Diego, CA, USA). Data are representative of four independent experiments and are expressed as means (errors in EC50 are within 10%).
The prediction of all the reported ADMET parameters was developed by means of the Advanced Chemistry Development (ACD) Percepta platform [70]. The software prediction is performed based on the software implemented training libraries, which include experimentally determined pharmacokinetic and safety properties for different series of compounds. Prediction of PAINS (Pan Assay Interference structures) was derived via the SwissADME website [71].

5. Conclusions

Recently, a number of combined structure-based studies and ligand-based studies have been attempted to explore the selectivity and species-specificity profile of different series of TAAR1 ligands.
Herein, based on the hTAAR1 AlphaFold protein model, we explored the effectiveness of structural variations involving the main scaffold of the conformationally locked hTAAR1 agonist S18616. The analogues 137 have been investigated via molecular docking studies, prior to proceeding with the biological evaluation of a series of in-house pyrimidinone-benzimidazole derivatives (1a10a). Based on the combined LB- and SBCM results, 1a10a have been assayed as putative TAAR1 agonists. Final biological evaluation revealed the piperazine-containing derivatives 1a3a (hTAAR1 EC50 = 526.3–657.4 nM) as the most promising of the series and allowed us to confirm the reliability of the computational studies. Then, in silico prediction of the compounds ADMET properties suggested the future development of the pyrimidinone-benzimidazole scaffold in the search for novel effective hTAAR1 agonists.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/molecules29081739/s1: Figure S1: MIF Comparison of the template S18686 and of the candidates 9a, 10a, 8a, and 7a; Figure S2: Common MIFs between the template S18686 and the candidates 1a and 7a; Figure S3: Molecular docking pose of compound 3a (C atom; violet) (A) and 6a (C atom; violet) (B) within the hTAAR1 binding site; Figure S4: Prediction of putative off-target preferences featured by 1a3a; Table S1: Chemical structure and biological activity as hTAAR1 agonists of compounds 137; Table S2: Chemical structure of in-house compounds 1a10a; Table S3: Ten top-scored docking positioning of 137 at the hTAAR1; Table S4: The calculated scoring function Glob-Prod values based on the MIFs H, N1, DRY, O are reported (compounds 1a10a as candidates, S18616 being the template); Table S5: The calculated scoring function Glob-Prod values based on the MIFs H, N1, DRY, O are reported (compounds 237 as candidates, S18616 being the template); Table S6: Ten top scored docking positioning of 1a10a at the hTAAR1; Table S7: Calculated descriptors related to toxicity properties, cytochrome inhibition, and lethal dose via mouse oral administration have been explored for 1a6a and for the reference S18616.

Author Contributions

Conceptualization, E.C. and S.S.; methodology, E.C. and R.R.G.; data curation, A.C., F.M. and E.K.; software, E.C. and P.F; funding acquisition, E.C., S.S., P.F. and R.R.G.; investigation, G.G. and N.S.; writing—original draft preparation, E.C. and R.R.G.; writing—review and editing, S.S., F.M. and M.T.; visualization, N.S., G.G. and E.K.; supervision, E.C. and A.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by FRA2022 from University of Genoa. E.K. and R.R.G. were supported by the Russian Science Foundation (19-75-30008-P).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data are contained within the article.

Conflicts of Interest

The authors declare no conflicts of interest.

Correction Statement

This article has been republished with minor corrections to the image quality of the chemical structures in Tables 1–4. This change does not affect the scientific content of the article.

References

  1. Gainetdinov, R.R.; Hoener, M.C.; Berry, M.D. Trace Amines and Their Receptors. Pharmacol. Rev. 2018, 70, 549–620. [Google Scholar] [CrossRef]
  2. Dedic, N.; Dworak, H.; Zeni, C.; Rutigliano, G.; Howes, O.D. Therapeutic Potential of TAAR1 Agonists in Schizophrenia: Evidence from Preclinical Models and Clinical Studies. Int. J. Mol. Sci. 2021, 22, 13185. [Google Scholar] [CrossRef] [PubMed]
  3. Hart, M.E.; Suchland, K.L.; Miyakawa, M.; Bunzow, J.R.; Grandy, D.K.; Scanlan, T.S. Trace Amine-Associated Receptor Agonists: Synthesis and Evaluation of Thyronamines and Related Analogues. J. Med. Chem. 2006, 49, 1101–1112. [Google Scholar] [CrossRef] [PubMed]
  4. Rutigliano, G.; Accorroni, A.; Zucchi, R. The Case for TAAR1 as a Modulator of Central Nervous System Function. Front. Pharmacol. 2017, 8, 987. [Google Scholar] [CrossRef]
  5. Xu, Z.; Li, Q. TAAR Agonists. Cell. Mol. Neurobiol. 2020, 40, 257–272. [Google Scholar] [CrossRef] [PubMed]
  6. Halff, E.F.; Rutigliano, G.; Garcia-Hidalgo, A.; Howes, O.D. Trace Amine-Associated Receptor 1 (TAAR1) Agonism as a New Treatment Strategy for Schizophrenia and Related Disorders. Trends Neurosci. 2023, 46, 60–74. [Google Scholar] [CrossRef] [PubMed]
  7. Guo, L.; Dai, W.; Xu, Z.; Liang, Q.; Miller, E.T.; Li, S.; Gao, X.; Baldwin, M.W.; Chai, R.; Li, Q. Evolution of Brain-Expressed Biogenic Amine Receptors into Olfactory Trace Amine-Associated Receptors. Mol. Biol. Evol. 2022, 39, msac006. [Google Scholar] [CrossRef] [PubMed]
  8. Guo, L.; Cheng, J.; Lian, S.; Liu, Q.; Lu, Y.; Zheng, Y.; Zhu, K.; Zhang, M.; Kong, Y.; Zhang, C.; et al. Structural Basis of Amine Odorant Perception by a Mammal Olfactory Receptor. Nature 2023, 618, 193–200. [Google Scholar] [CrossRef]
  9. Espinoza, S.; Sukhanov, I.; Efimova, E.V.; Kozlova, A.; Antonova, K.A.; Illiano, P.; Leo, D.; Merkulyeva, N.; Kalinina, D.; Musienko, P.; et al. Trace Amine-Associated Receptor 5 Provides Olfactory Input Into Limbic Brain Areas and Modulates Emotional Behaviors and Serotonin Transmission. Front. Mol. Neurosci. 2020, 13, 18. [Google Scholar] [CrossRef]
  10. Efimova, E.V.; Kozlova, A.A.; Razenkova, V.; Katolikova, N.V.; Antonova, K.A.; Sotnikova, T.D.; Merkulyeva, N.S.; Veshchitskii, A.S.; Kalinina, D.S.; Korzhevskii, D.E.; et al. Increased Dopamine Transmission and Adult Neurogenesis in Trace Amine-Associated Receptor 5 (TAAR5) Knockout Mice. Neuropharmacology 2021, 182, 108373. [Google Scholar] [CrossRef] [PubMed]
  11. Efimova, E.V.; Kuvarzin, S.R.; Mor, M.S.; Katolikova, N.V.; Shemiakova, T.S.; Razenkova, V.; Ptukha, M.; Kozlova, A.A.; Murtazina, R.Z.; Smirnova, D.; et al. Trace Amine-Associated Receptor 2 Is Expressed in the Limbic Brain Areas and Is Involved in Dopamine Regulation and Adult Neurogenesis. Front. Behav. Neurosci. 2022, 16, 847410. [Google Scholar] [CrossRef] [PubMed]
  12. Borowsky, B.; Adham, N.; Jones, K.A.; Raddatz, R.; Artymyshyn, R.; Ogozalek, K.L.; Durkin, M.M.; Lakhlani, P.P.; Bonini, J.A.; Pathirana, S.; et al. Trace Amines: Identification of a Family of Mammalian G Protein-Coupled Receptors. Proc. Natl. Acad. Sci. USA 2001, 98, 8966–8971. [Google Scholar] [CrossRef] [PubMed]
  13. Bunzow, J.R.; Sonders, M.S.; Arttamangkul, S.; Harrison, L.M.; Zhang, G.; Quigley, D.I.; Darland, T.; Suchland, K.L.; Pasumamula, S.; Kennedy, J.L.; et al. Amphetamine, 3,4-Methylenedioxymethamphetamine, Lysergic Acid Diethylamide, and Metabolites of the Catecholamine Neurotransmitters Are Agonists of a Rat Trace Amine Receptor. Mol. Pharmacol. 2001, 60, 1181–1188. [Google Scholar] [CrossRef] [PubMed]
  14. Lindemann, L.; Meyer, C.A.; Jeanneau, K.; Bradaia, A.; Ozmen, L.; Bluethmann, H.; Bettler, B.; Wettstein, J.G.; Borroni, E.; Moreau, J.-L.; et al. Trace Amine-Associated Receptor 1 Modulates Dopaminergic Activity. J. Pharmacol. Exp. Ther. 2008, 324, 948–956. [Google Scholar] [CrossRef]
  15. John, J.; Kukshal, P.; Bhatia, T.; Chowdari, K.V.; Nimgaonkar, V.L.; Deshpande, S.N.; Thelma, B.K. Possible Role of Rare Variants in Trace Amine Associated Receptor 1 in Schizophrenia. Schizophr. Res. 2017, 189, 190–195. [Google Scholar] [CrossRef] [PubMed]
  16. Mühlhaus, J.; Dinter, J.; Jyrch, S.; Teumer, A.; Jacobi, S.F.; Homuth, G.; Kühnen, P.; Wiegand, S.; Grüters, A.; Völzke, H.; et al. Investigation of Naturally Occurring Single-Nucleotide Variants in Human TAAR1. Front. Pharmacol. 2017, 8, 807. [Google Scholar] [CrossRef] [PubMed]
  17. Imbriglio, T.; Alborghetti, M.; Bruno, V.; Battaglia, G.; Nicoletti, F.; Cannella, M. Up-Regulation of the Trace Amine Receptor, TAAR-1, in the Prefrontal Cortex of Individuals Affected by Schizophrenia. Schizophr. Bull. 2023, 50, 374–381. [Google Scholar] [CrossRef] [PubMed]
  18. De Gregorio, D.; Posa, L.; Ochoa-Sanchez, R.; McLaughlin, R.; Maione, S.; Comai, S.; Gobbi, G. The Hallucinogen D-Lysergic Diethylamide (LSD) Decreases Dopamine Firing Activity through 5-HT1A, D2 and TAAR1 Receptors. Pharmacol. Res. 2016, 113, 81–91. [Google Scholar] [CrossRef]
  19. Liberles, S.D. Trace Amine-Associated Receptors: Ligands, Neural Circuits, and Behaviors. Curr. Opin. Neurobiol. 2015, 34, 1–7. [Google Scholar] [CrossRef]
  20. Polini, B.; Ricardi, C.; Bertolini, A.; Carnicelli, V.; Rutigliano, G.; Saponaro, F.; Zucchi, R.; Chiellini, G. T1AM/TAAR1 System Reduces Inflammatory Response and β-Amyloid Toxicity in Human Microglial HMC3 Cell Line. Int. J. Mol. Sci. 2023, 24, 11569. [Google Scholar] [CrossRef]
  21. Frycz, B.A.; Nowicka, K.; Konopka, A.; Hoener, M.C.; Bulska, E.; Kaczmarek, L.; Stefaniuk, M. Activation of Trace Amine-associated Receptor 1 (TAAR1) Transiently Reduces Alcohol Drinking in Socially Housed Mice. Addict. Biol. 2023, 28, e13285. [Google Scholar] [CrossRef] [PubMed]
  22. Berry, M.D.; Gainetdinov, R.R.; Hoener, M.C.; Shahid, M. Pharmacology of Human Trace Amine-Associated Receptors: Therapeutic Opportunities and Challenges. Pharmacol. Ther. 2017, 180, 161–180. [Google Scholar] [CrossRef] [PubMed]
  23. Jia, L.; Li, S.; Dai, W.; Guo, L.; Xu, Z.; Scott, A.M.; Zhang, Z.; Ren, J.; Zhang, Q.; Dexheimer, T.S.; et al. Convergent Olfactory Trace Amine-Associated Receptors Detect Biogenic Polyamines with Distinct Motifs via a Conserved Binding Site. J. Biol. Chem. 2021, 297, 101268. [Google Scholar] [CrossRef] [PubMed]
  24. Cisneros, I.E.; Ghorpade, A.; Borgmann, K. Methamphetamine Activates Trace Amine Associated Receptor 1 to Regulate Astrocyte Excitatory Amino Acid Transporter-2 via Differential CREB Phosphorylation During HIV-Associated Neurocognitive Disorders. Front. Neurol. 2020, 11, 593146. [Google Scholar] [CrossRef]
  25. Freyberg, Z.; Saavedra, J.M. Trace Amines and Trace Amine-Associated Receptors: A New Frontier in Cell Signaling. Cell. Mol. Neurobiol. 2020, 40, 189–190. [Google Scholar] [CrossRef]
  26. Zhai, L.; Xiao, H.; Lin, C.; Wong, H.L.X.; Lam, Y.Y.; Gong, M.; Wu, G.; Ning, Z.; Huang, C.; Zhang, Y.; et al. Gut Microbiota-Derived Tryptamine and Phenethylamine Impair Insulin Sensitivity in Metabolic Syndrome and Irritable Bowel Syndrome. Nat. Commun. 2023, 14, 4986. [Google Scholar] [CrossRef]
  27. Zheng, Y.; Yasuda, M.; Yamao, M.; Gokan, T.; Sejima, Y.; Nishikawa, T.; Katayama, S. Fermented Soybean Foods (Natto) Ameliorate Age-Related Cognitive Decline by Hippocampal TAAR1-Mediated Activation of the CaMKII/CREB/BDNF Signaling Pathway in Senescence-Accelerated Mouse Prone 8 (SAMP8). Food Funct. 2023, 14, 10097–10106. [Google Scholar] [CrossRef] [PubMed]
  28. Liu, J.; Wu, R.; Li, J.-X. TAAR1 as an Emerging Target for the Treatment of Psychiatric Disorders. Pharmacol. Ther. 2024, 253, 108580. [Google Scholar] [CrossRef]
  29. Cichero, E.; Tonelli, M. Targeting Species-Specific Trace Amine-Associated Receptor 1 Ligands: To Date Perspective of the Rational Drug Design Process. Future Med. Chem. 2017, 9, 1507–1527. [Google Scholar] [CrossRef]
  30. Wainscott, D.B.; Little, S.P.; Yin, T.; Tu, Y.; Rocco, V.P.; He, J.X.; Nelson, D.L. Pharmacologic Characterization of the Cloned Human Trace Amine-Associated Receptor1 (TAAR1) and Evidence for Species Differences with the Rat TAAR1. J. Pharmacol. Exp. Ther. 2007, 320, 475–485. [Google Scholar] [CrossRef]
  31. Reese, E.A.; Bunzow, J.R.; Arttamangkul, S.; Sonders, M.S.; Grandy, D.K. Trace Amine-Associated Receptor 1 Displays Species-Dependent Stereoselectivity for Isomers of Methamphetamine, Amphetamine, and Para-Hydroxyamphetamine. J. Pharmacol. Exp. Ther. 2007, 321, 178–186. [Google Scholar] [CrossRef] [PubMed]
  32. Tan, E.S.; Miyakawa, M.; Bunzow, J.R.; Grandy, D.K.; Scanlan, T.S. Exploring the Structure−Activity Relationship of the Ethylamine Portion of 3-Iodothyronamine for Rat and Mouse Trace Amine-Associated Receptor 1. J. Med. Chem. 2007, 50, 2787–2798. [Google Scholar] [CrossRef] [PubMed]
  33. Guariento, S.; Tonelli, M.; Espinoza, S.; Gerasimov, A.S.; Gainetdinov, R.R.; Cichero, E. Rational Design, Chemical Synthesis and Biological Evaluation of Novel Biguanides Exploring Species-Specificity Responsiveness of TAAR1 Agonists. Eur. J. Med. Chem. 2018, 146, 171–184. [Google Scholar] [CrossRef]
  34. Tonelli, M.; Espinoza, S.; Gainetdinov, R.R.; Cichero, E. Novel Biguanide-Based Derivatives Scouted as TAAR1 Agonists: Synthesis, Biological Evaluation, ADME Prediction and Molecular Docking Studies. Eur. J. Med. Chem. 2017, 127, 781–792. [Google Scholar] [CrossRef]
  35. Francesconi, V.; Cichero, E.; Kanov, E.V.; Laurini, E.; Pricl, S.; Gainetdinov, R.R.; Tonelli, M. Novel 1-Amidino-4-Phenylpiperazines as Potent Agonists at Human TAAR1 Receptor: Rational Design, Synthesis, Biological Evaluation and Molecular Docking Studies. Pharmaceuticals 2020, 13, 391. [Google Scholar] [CrossRef] [PubMed]
  36. Galley, G.; Stalder, H.; Goergler, A.; Hoener, M.C.; Norcross, R.D. Optimisation of Imidazole Compounds as Selective TAAR1 Agonists: Discovery of RO5073012. Bioorg. Med. Chem. Lett. 2012, 22, 5244–5248. [Google Scholar] [CrossRef] [PubMed]
  37. Simmler, L.D.; Rickli, A.; Schramm, Y.; Hoener, M.C.; Liechti, M.E. Pharmacological Profiles of Aminoindanes, Piperazines, and Pipradrol Derivatives. Biochem. Pharmacol. 2014, 88, 237–244. [Google Scholar] [CrossRef]
  38. Galley, G.; Beurier, A.; Décoret, G.; Goergler, A.; Hutter, R.; Mohr, S.; Pähler, A.; Schmid, P.; Türck, D.; Unger, R.; et al. Discovery and Characterization of 2-Aminooxazolines as Highly Potent, Selective, and Orally Active TAAR1 Agonists. ACS Med. Chem. Lett. 2016, 7, 192–197. [Google Scholar] [CrossRef] [PubMed]
  39. Kotańska, M.; Marcinkowska, M.; Kuder, K.J.; Walczak, M.; Bednarski, M.; Siwek, A.; Kołaczkowski, M. Metabolic and Cardiovascular Benefits and Risks of 4-Hydroxy Guanabenz Hydrochloride: A2-Adrenoceptor and Trace Amine-Associated Receptor 1 Ligand. Pharmacol. Rep. 2023, 75, 1211–1229. [Google Scholar] [CrossRef] [PubMed]
  40. Gowtham, L.; Halder, N.; Angmo, D.; Singh, S.B.; Jayasundar, R.; Dada, T.; Velpandian, T. Untargeted Metabolomics in the Aqueous Humor Reveals the Involvement of TAAR Pathway in Glaucoma. Exp. Eye Res. 2023, 234, 109592. [Google Scholar] [CrossRef]
  41. Dedic, N.; Jones, P.G.; Hopkins, S.C.; Lew, R.; Shao, L.; Campbell, J.E.; Spear, K.L.; Large, T.H.; Campbell, U.C.; Hanania, T.; et al. SEP-363856, a Novel Psychotropic Agent with a Unique, Non-D2 Receptor Mechanism of Action. J. Pharmacol. Exp. Ther. 2019, 371, 1–14. [Google Scholar] [CrossRef]
  42. Achtyes, E.D.; Hopkins, S.C.; Dedic, N.; Dworak, H.; Zeni, C.; Koblan, K. Ulotaront: Review of Preliminary Evidence for the Efficacy and Safety of a TAAR1 Agonist in Schizophrenia. Eur. Arch. Psychiatry Clin. Neurosci. 2023, 273, 1543–1556. [Google Scholar] [CrossRef] [PubMed]
  43. Millan, M.J.; Dekeyne, A.; Newman-Tancredi, A.; Cussac, D.; Audinot, V.; Milligan, G.; Duqueyroix, D.; Girardon, S.; Mullot, J.; Boutin, J.A.; et al. S18616, a Highly Potent, Spiroimidazoline Agonist at Alpha(2)-Adrenoceptors: I. Receptor Profile, Antinociceptive and Hypothermic Actions in Comparison with Dexmedetomidine and Clonidine. J. Pharmacol. Exp. Ther. 2000, 295, 1192–1205. [Google Scholar] [PubMed]
  44. Pérez-Peña, H.; Abel, A.-C.; Shevelev, M.; Prota, A.E.; Pieraccini, S.; Horvath, D. Computational Approaches to the Rational Design of Tubulin-Targeting Agents. Biomolecules 2023, 13, 285. [Google Scholar] [CrossRef]
  45. Dorahy, G.; Chen, J.Z.; Balle, T. Computer-Aided Drug Design towards New Psychotropic and Neurological Drugs. Molecules 2023, 28, 1324. [Google Scholar] [CrossRef]
  46. Bon, C.; Chern, T.-R.; Cichero, E.; O’Brien, T.E.; Gustincich, S.; Gainetdinov, R.R.; Espinoza, S. Discovery of Novel Trace Amine-Associated Receptor 5 (TAAR5) Antagonists Using a Deep Convolutional Neural Network. Int. J. Mol. Sci. 2022, 23, 3127. [Google Scholar] [CrossRef]
  47. Galati, S.; Di Stefano, M.; Bertini, S.; Granchi, C.; Giordano, A.; Gado, F.; Macchia, M.; Tuccinardi, T.; Poli, G. Identification of New GSK3β Inhibitors through a Consensus Machine Learning-Based Virtual Screening. Int. J. Mol. Sci. 2023, 24, 17233. [Google Scholar] [CrossRef] [PubMed]
  48. Galati, S.; Sainas, S.; Giorgis, M.; Boschi, D.; Lolli, M.L.; Ortore, G.; Poli, G.; Tuccinardi, T. Identification of Human Dihydroorotate Dehydrogenase Inhibitor by a Pharmacophore-Based Virtual Screening Study. Molecules 2022, 27, 3660. [Google Scholar] [CrossRef]
  49. Poli, G.; Di Stefano, M.; Estevez, J.A.; Minutolo, F.; Granchi, C.; Giordano, A.; Parisi, S.; Mauceri, M.; Canzonieri, V.; Macchia, M.; et al. New PIN1 Inhibitors Identified through a Pharmacophore-Driven, Hierarchical Consensus Docking Strategy. J. Enzym. Inhib. Med. Chem. 2022, 37, 145–150. [Google Scholar] [CrossRef]
  50. Baselious, F.; Hilscher, S.; Robaa, D.; Barinka, C.; Schutkowski, M.; Sippl, W. Comparative Structure-Based Virtual Screening Utilizing Optimized AlphaFold Model Identifies Selective HDAC11 Inhibitor. Int. J. Mol. Sci. 2024, 25, 1358. [Google Scholar] [CrossRef]
  51. Liao, S.; Pino, M.J.; Deleon, C.; Lindner-Jackson, M.; Wu, C. Interaction Analyses of HTAAR1 and MTAAR1 with Antagonist EPPTB. Life Sci. 2022, 300, 120553. [Google Scholar] [CrossRef]
  52. Cichero, E.; Espinoza, S.; Franchini, S.; Guariento, S.; Brasili, L.; Gainetdinov, R.R.; Fossa, P. Further Insights Into the Pharmacology of the Human Trace Amine-Associated Receptors: Discovery of Novel Ligands for TAAR1 by a Virtual Screening Approach. Chem. Biol. Drug Des. 2014, 84, 712–720. [Google Scholar] [CrossRef] [PubMed]
  53. Cichero, E.; Espinoza, S.; Gainetdinov, R.R.; Brasili, L.; Fossa, P. Insights into the Structure and Pharmacology of the Human Trace Amine-Associated Receptor 1 (HTAAR1): Homology Modelling and Docking Studies. Chem. Biol. Drug Des. 2013, 81, 509–516. [Google Scholar] [CrossRef] [PubMed]
  54. Reese, E.A.; Norimatsu, Y.; Grandy, M.S.; Suchland, K.L.; Bunzow, J.R.; Grandy, D.K. Exploring the Determinants of Trace Amine-Associated Receptor 1′s Functional Selectivity for the Stereoisomers of Amphetamine and Methamphetamine. J. Med. Chem. 2014, 57, 378–390. [Google Scholar] [CrossRef]
  55. Jumper, J.; Evans, R.; Pritzel, A.; Green, T.; Figurnov, M.; Ronneberger, O.; Tunyasuvunakool, K.; Bates, R.; Žídek, A.; Potapenko, A.; et al. Highly Accurate Protein Structure Prediction with AlphaFold. Nature 2021, 596, 583–589. [Google Scholar] [CrossRef] [PubMed]
  56. Varadi, M.; Anyango, S.; Deshpande, M.; Nair, S.; Natassia, C.; Yordanova, G.; Yuan, D.; Stroe, O.; Wood, G.; Laydon, A.; et al. AlphaFold Protein Structure Database: Massively Expanding the Structural Coverage of Protein-Sequence Space with High-Accuracy Models. Nucleic Acids Res. 2022, 50, D439–D444. [Google Scholar] [CrossRef]
  57. Tunyasuvunakool, K.; Adler, J.; Wu, Z.; Green, T.; Zielinski, M.; Žídek, A.; Bridgland, A.; Cowie, A.; Meyer, C.; Laydon, A.; et al. Highly Accurate Protein Structure Prediction for the Human Proteome. Nature 2021, 596, 590–596. [Google Scholar] [CrossRef]
  58. Kosoglu, K.; Aydin, Z.; Tuncbag, N.; Gursoy, A.; Keskin, O. Structural Coverage of the Human Interactome. Brief. Bioinform. 2023, 25, bbad496. [Google Scholar] [CrossRef]
  59. Di Braccio, M.; Grossi, G.; Signorello, M.G.; Leoncini, G.; Cichero, E.; Fossa, P.; Alfei, S.; Damonte, G. Synthesis, in Vitro Antiplatelet Activity and Molecular Modelling Studies of 10-Substituted 2-(1-Piperazinyl)Pyrimido[1,2-a]Benzimidazol-4(10H)-Ones. Eur. J. Med. Chem. 2013, 62, 564–578. [Google Scholar] [CrossRef] [PubMed]
  60. Baroni, M.; Cruciani, G.; Sciabola, S.; Perruccio, F.; Mason, J.S. A Common Reference Framework for Analyzing/Comparing Proteins and Ligands. Fingerprints for Ligands And Proteins (FLAP): Theory and Application. J. Chem. Inf. Model. 2007, 47, 279–294. [Google Scholar] [CrossRef]
  61. Cross, S.; Baroni, M.; Goracci, L.; Cruciani, G. GRID-Based Three-Dimensional Pharmacophores I: FLAPpharm, a Novel Approach for Pharmacophore Elucidation. J. Chem. Inf. Model. 2012, 52, 2587–2598. [Google Scholar] [CrossRef] [PubMed]
  62. Cichero, E.; Francesconi, V.; Casini, B.; Casale, M.; Kanov, E.; Gerasimov, A.S.; Sukhanov, I.; Savchenko, A.; Espinoza, S.; Gainetdinov, R.R.; et al. Discovery of Guanfacine as a Novel TAAR1 Agonist: A Combination Strategy through Molecular Modeling Studies and Biological Assays. Pharmaceuticals 2023, 16, 1632. [Google Scholar] [CrossRef] [PubMed]
  63. Barak, L.S.; Salahpour, A.; Zhang, X.; Masri, B.; Sotnikova, T.D.; Ramsey, A.J.; Violin, J.D.; Lefkowitz, R.J.; Caron, M.G.; Gainetdinov, R.R. Pharmacological Characterization of Membrane-Expressed Human Trace Amine-Associated Receptor 1 (TAAR1) by a Bioluminescence Resonance Energy Transfer CAMP Biosensor. Mol. Pharmacol. 2008, 74, 585–594. [Google Scholar] [CrossRef] [PubMed]
  64. Moro, S.; Bacilieri, M.; Cacciari, B.; Bolcato, C.; Cusan, C.; Pastorin, G.; Klotz, K.-N.; Spalluto, G. The Application of a 3D-QSAR (AutoMEP/PLS) Approach as an Efficient Pharmacodynamic-Driven Filtering Method for Small-Sized Virtual Library: Application to a Lead Optimization of a Human A3 Adenosine Receptor Antagonist. Bioorg. Med. Chem. 2006, 14, 4923–4932. [Google Scholar] [CrossRef] [PubMed]
  65. Cichero, E.; Menozzi, G.; Spallarossa, A.; Mosti, L.; Fossa, P. Exploring the Binding Features of Rimonabant Analogues and Acyclic CB1 Antagonists: Docking Studies and QSAR Analysis. J. Mol. Model. 2008, 14, 1131–1145. [Google Scholar] [CrossRef] [PubMed]
  66. Nair, P.C.; Chalker, J.M.; McKinnon, R.A.; Langmead, C.J.; Gregory, K.J.; Bastiampillai, T. Trace Amine-Associated Receptor 1 (TAAR1): Molecular and Clinical Insights for the Treatment of Schizophrenia and Related Comorbidities. ACS Pharmacol. Transl. Sci. 2022, 5, 183–188. [Google Scholar] [CrossRef]
  67. Chemical Computing Group ULC Molecular Operating Environment (MOE2019.01) 2021. Available online: http://www.chemcomp.com/ (accessed on 15 February 2024).
  68. Dulsat, J.; López-Nieto, B.; Estrada-Tejedor, R.; Borrell, J.I. Evaluation of Free Online ADMET Tools for Academic or Small Biotech Environments. Molecules 2023, 28, 776. [Google Scholar] [CrossRef] [PubMed]
  69. Parodi, A.; Righetti, G.; Pesce, E.; Salis, A.; Tasso, B.; Urbinati, C.; Tomati, V.; Damonte, G.; Rusnati, M.; Pedemonte, N.; et al. Discovery of Novel VX-809 Hybrid Derivatives as F508del-CFTR Correctors by Molecular Modeling, Chemical Synthesis and Biological Assays. Eur. J. Med. Chem. 2020, 208, 112833. [Google Scholar] [CrossRef]
  70. Advanced Chemistry Development, Inc. ACD/Percepta Platform; Advanced Chemistry Development, Inc.: Toronto, Canada, 2015. [Google Scholar]
  71. Gfeller, D.; Michielin, O.; Zoete, V. Shaping the Interaction Landscape of Bioactive Molecules. Bioinformatics 2013, 29, 3073–3079. [Google Scholar] [CrossRef]
  72. Chai, Y.; Zhang, E.; Cai, Z.; Xu, D.; Zhu, C.; Sun, B. Isolation, Synthesis and Identification of Degraded Impurities in Letermovir. J. Pharm. Biomed. Anal. 2023, 236, 115691. [Google Scholar] [CrossRef]
  73. Vezenkov, L.T.; Tsekova, D.S.; Kostadinova, I.; Mihaylova, R.; Lozanov, V.; Vassilev, N.G.; Danchev, N.; Tsakovska, I.; Pajeva, I. Synthesis and Biological Study of 4-Aminopyridine-Peptide Derivatives Designed for the Treatment of Neurodegenerative Disorders. Curr. Alzheimer Res. 2023, 20, 120–129. [Google Scholar] [CrossRef] [PubMed]
  74. Maciejewska, K.; Czarnecka, K.; Kręcisz, P.; Niedziałek, D.; Wieczorek, G.; Skibiński, R.; Szymański, P. Novel Cyclopentaquinoline and Acridine Analogs as Multifunctional, Potent Drug Candidates in Alzheimer’s Disease. Int. J. Mol. Sci. 2022, 23, 5876. [Google Scholar] [CrossRef] [PubMed]
  75. Chen, F.; Yin, X.; Wang, Y.; Lv, Y.; Sheng, S.; Ouyang, S.; Zhong, Y. Pharmacokinetics, Tissue Distribution, and Druggability Prediction of the Natural Anticancer Active Compound Cytisine N-Isoflavones Combined with Computer Simulation. Biol. Pharm. Bull. 2020, 43, 976–984. [Google Scholar] [CrossRef] [PubMed]
  76. Huang, C.; Zhang, Y.; Xu, Y.; Wei, S.; Yang, T.; Wang, S.; Li, C.; Lin, H.; Li, X.; Zhao, S.; et al. Prepared Radix Polygoni Multiflori and Emodin Alleviate Lipid Droplet Accumulation in Nonalcoholic Fatty Liver Disease through MAPK Signaling Pathway Inhibition. Aging 2024, 16, 2362–2384. [Google Scholar] [CrossRef] [PubMed]
  77. Bepari, A.K.; Shatabda, S.; Reza, H.M. Virtual Screening of Flavonoids as Potential RIPK1 Inhibitors for Neurodegeneration Therapy. PeerJ 2024, 12, e16762. [Google Scholar] [CrossRef]
  78. Ahmad, V.; Khan, M.I.; Jamal, Q.M.S.; Alzahrani, F.A.; Albiheyri, R. Computational Molecular Docking and Simulation-Based Assessment of Anti-Inflammatory Properties of Nyctanthes Arbor-Tristis Linn Phytochemicals. Pharmaceuticals 2023, 17, 18. [Google Scholar] [CrossRef] [PubMed]
  79. Veber, D.F.; Johnson, S.R.; Cheng, H.-Y.; Smith, B.R.; Ward, K.W.; Kopple, K.D. Molecular Properties That Influence the Oral Bioavailability of Drug Candidates. J. Med. Chem. 2002, 45, 2615–2623. [Google Scholar] [CrossRef] [PubMed]
  80. Lipinski, C.A.; Lombardo, F.; Dominy, B.W.; Feeney, P.J. Experimental and Computational Approaches to Estimate Solubility and Permeability in Drug Discovery and Development Settings. Adv. Drug Deliv. Rev. 2012, 64, 4–17. [Google Scholar] [CrossRef]
  81. Hou, C.; Wen, X.; Yan, S.; Gu, X.; Jiang, Y.; Chen, F.; Liu, Y.; Zhu, Y.; Liu, X. Network-Based Pharmacology-Based Research on the Effect and Mechanism of the Hedyotis Diffusa–Scutellaria Barbata Pair in the Treatment of Hepatocellular Carcinoma. Sci. Rep. 2024, 14, 963. [Google Scholar] [CrossRef] [PubMed]
  82. Mao, J.; Tang, L.; Fang, L.; Tian, C.; Zhu, Z.; Li, Y. Systematic Pharmacology-Based Strategy to Explore the Mechanism of Semen Strychni for Treatment of Papillary Thyroid Carcinoma. Sci. Rep. 2023, 13, 18492. [Google Scholar] [CrossRef] [PubMed]
  83. Wu, X.; Zheng, X.; Tang, H.; Zhao, L.; He, C.; Zou, Y.; Song, X.; Li, L.; Yin, Z.; Ye, G. A Network Pharmacology Approach to Identify the Mechanisms and Molecular Targets of Curcumin against Alzheimer Disease. Medicine 2022, 101, e30194. [Google Scholar] [CrossRef]
  84. Yang, X.; Guan, Y.; Yan, B.; Xie, Y.; Zhou, M.; Wu, Y.; Yao, L.; Qiu, X.; Yan, F.; Chen, Y.; et al. Evidence-based Complementary and Alternative Medicine Bioinformatics Approach through Network Pharmacology and Molecular Docking to Determine the Molecular Mechanisms of Erjing Pill in Alzheimer’s Disease. Exp. Ther. Med. 2021, 22, 1252. [Google Scholar] [CrossRef] [PubMed]
  85. Asproni, B.; Catto, M.; Loriga, G.; Murineddu, G.; Corona, P.; Purgatorio, R.; Cichero, E.; Fossa, P.; Scarano, N.; Martínez, A.L.; et al. Novel Thienocycloalkylpyridazinones as Useful Scaffolds for Acetylcholinesterase Inhibition and Serotonin 5-HT6 Receptor Interaction. Bioorg. Med. Chem. 2023, 84, 117256. [Google Scholar] [CrossRef] [PubMed]
  86. Djeujo, F.M.; Ragazzi, E.; Urettini, M.; Sauro, B.; Cichero, E.; Tonelli, M.; Froldi, G. Magnolol and Luteolin Inhibition of α-Glucosidase Activity: Kinetics and Type of Interaction Detected by In Vitro and In Silico Studies. Pharmaceuticals 2022, 15, 205. [Google Scholar] [CrossRef] [PubMed]
  87. Poongavanam, V.; Kongsted, J. Virtual Screening Models for Prediction of HIV-1 RT Associated RNase H Inhibition. PLoS ONE 2013, 8, e73478. [Google Scholar] [CrossRef]
  88. Zheng, C.; Javitch, J.A.; Lambert, N.A.; Donthamsetti, P.; Gurevich, V.V. In-Cell Arrestin-Receptor Interaction Assays. Curr. Protoc. 2023, 3, e890. [Google Scholar] [CrossRef] [PubMed]
  89. Ti, R.; Pang, B.; Yu, L.; Gan, B.; Ma, W.; Warshel, A.; Ren, R.; Zhu, L. Fine-tuning activation specificity of G-protein-coupled receptors via automated path searching. Proc. Natl. Acad. Sci. USA 2024, 121, e2317893121. [Google Scholar] [CrossRef] [PubMed]
  90. Calabretta, M.M.; Michelini, E. Current advances in the use of bioluminescence assays for drug discovery: An update of the last ten years. Expert. Opin. Drug Discov. 2024, 19, 85. [Google Scholar] [CrossRef] [PubMed]
  91. Compan, V.; Pelegrín, P. Measuring IL-1β Processing by Bioluminescence Sensors: Using a Bioluminescence Resonance Energy Transfer Biosensor. Methods Mol. Biol. 2023, 2696, 47. [Google Scholar] [CrossRef]
  92. Janicot, R.; Park, J.C.; Garcia-Marcos, M. Detecting GPCR Signals with Optical Biosensors of Gα-GTP in Cell Lines and Primary Cell Cultures. Curr. Protoc. 2023, 3, e796. [Google Scholar] [CrossRef]
  93. Cho, K.F.; Javier, N.; Choi, K. BRET measurement on CCD camera-based microtiter plate readers. SLAS Discov. 2022, 27, 413. [Google Scholar] [CrossRef] [PubMed]
  94. Mujawar, A.; De, A. In Vivo Assessment of Protein-Protein Interactions Using BRET Assay. Methods Mol. Biol. 2022, 2525, 239. [Google Scholar] [CrossRef]
Figure 1. Structural variations performed at the S18616 1 (A) [43] compound as TAAR1 prototype towards further oxazoline-based agonists 237 (B) [38]. The promising compounds 5, 11, and 35 have been pointed out (C). The chemical structure of the herein proposed TAAR1 agonists 1a10a is shown (D). The corresponding hTAAR1 EC50 values are reported.
Figure 1. Structural variations performed at the S18616 1 (A) [43] compound as TAAR1 prototype towards further oxazoline-based agonists 237 (B) [38]. The promising compounds 5, 11, and 35 have been pointed out (C). The chemical structure of the herein proposed TAAR1 agonists 1a10a is shown (D). The corresponding hTAAR1 EC50 values are reported.
Molecules 29 01739 g001
Figure 2. Molecular docking pose of compound S18616 1 [43] (C atom; light green) within the hTAAR1 binding site (A). The most important residues are shown and labeled. H-bonds are highlighted as dashed yellow lines. The TAAR1 agonist chemical structure and the ligand hydrophobicity/electrostatic properties are also reported (hydrophobic and mild polar areas are shown in green and cyan, respectively; electron-rich and -poor groups are highlighted in red and blue). (B). The top view and side view of the receptor in the presence of S18616 1 are shown in (C,D), respectively. The H-bonding and hydrophobic areas of the protein are colored magenta and green. All the surface representations have been produced via MOE2019.01 software [67].
Figure 2. Molecular docking pose of compound S18616 1 [43] (C atom; light green) within the hTAAR1 binding site (A). The most important residues are shown and labeled. H-bonds are highlighted as dashed yellow lines. The TAAR1 agonist chemical structure and the ligand hydrophobicity/electrostatic properties are also reported (hydrophobic and mild polar areas are shown in green and cyan, respectively; electron-rich and -poor groups are highlighted in red and blue). (B). The top view and side view of the receptor in the presence of S18616 1 are shown in (C,D), respectively. The H-bonding and hydrophobic areas of the protein are colored magenta and green. All the surface representations have been produced via MOE2019.01 software [67].
Molecules 29 01739 g002
Figure 3. Molecular docking pose of compound 6 (C atom; violet) (A), 7 (C atom; violet) (B), 18 (C atom; violet) (C), and 17 (C atom; violet) (D) within the hTAAR1 binding site. The docking pose of the reference compound S18616 1 [43] is reported (C atom; light green). The most important residues are shown and labeled. H-bonds are highlighted as dashed yellow lines. The chemical structures of the TAAR1 agonists 6, 7 [38] and 17, 18 [38] are reported.
Figure 3. Molecular docking pose of compound 6 (C atom; violet) (A), 7 (C atom; violet) (B), 18 (C atom; violet) (C), and 17 (C atom; violet) (D) within the hTAAR1 binding site. The docking pose of the reference compound S18616 1 [43] is reported (C atom; light green). The most important residues are shown and labeled. H-bonds are highlighted as dashed yellow lines. The chemical structures of the TAAR1 agonists 6, 7 [38] and 17, 18 [38] are reported.
Molecules 29 01739 g003
Figure 4. Molecular docking pose of compound 11 (C atom; violet) (A), 21 (C atom; violet) (B), 23 (C atom; violet) (C), and 34 (C atom; violet) (D) within the hTAAR1 binding site. The docking pose of the reference compound S18616 1 [43] is reported (C atom; light green). The most important residues are shown and labeled. H-bonds are highlighted as dashed yellow lines. The chemical structures of the TAAR1 agonists 11, 21 [38] and 23, 34 [38] are reported.
Figure 4. Molecular docking pose of compound 11 (C atom; violet) (A), 21 (C atom; violet) (B), 23 (C atom; violet) (C), and 34 (C atom; violet) (D) within the hTAAR1 binding site. The docking pose of the reference compound S18616 1 [43] is reported (C atom; light green). The most important residues are shown and labeled. H-bonds are highlighted as dashed yellow lines. The chemical structures of the TAAR1 agonists 11, 21 [38] and 23, 34 [38] are reported.
Molecules 29 01739 g004
Figure 5. Schematic comparison of the herein investigated series of compounds 1a10a [59], with respect to the hTAAR1 prototype S18616 1 [43].
Figure 5. Schematic comparison of the herein investigated series of compounds 1a10a [59], with respect to the hTAAR1 prototype S18616 1 [43].
Molecules 29 01739 g005
Figure 6. MIF Comparison of the template S18686 1 [43] (C atom; green) MIFs (solid) and those (wireframes) of the candidates [59] (C atom; magenta). Blue: N1 MIF, red: O MIF, yellow: DRY MIF. The analysis is reported for compounds 1a (A), 3a (B), 2a (C), 6a (D), 4a (E), and 5a (F) based on the calculated Glob-Prod ranking (1a, Glob-Prod = 0.3499; 3a, Glob-Prod = 0.3167; 2a, Glob-Prod = 0.3154; 6a, Glob-Prod = 0.2962; 4a, Glob-Prod = 0.2841; 5a, Glob-Prod = 0.2838).
Figure 6. MIF Comparison of the template S18686 1 [43] (C atom; green) MIFs (solid) and those (wireframes) of the candidates [59] (C atom; magenta). Blue: N1 MIF, red: O MIF, yellow: DRY MIF. The analysis is reported for compounds 1a (A), 3a (B), 2a (C), 6a (D), 4a (E), and 5a (F) based on the calculated Glob-Prod ranking (1a, Glob-Prod = 0.3499; 3a, Glob-Prod = 0.3167; 2a, Glob-Prod = 0.3154; 6a, Glob-Prod = 0.2962; 4a, Glob-Prod = 0.2841; 5a, Glob-Prod = 0.2838).
Molecules 29 01739 g006
Figure 7. Molecular docking pose of compound 1a (C atom; violet) (A), 2a (C atom; violet) (B), 4a (C atom; violet) (C), and 5a (C atom; violet) (D) within the hTAAR1 binding site. The docking pose of the reference compound S18616 1 [43] is reported (C atom; light green). The most important residues are shown and labeled. H-bonds are highlighted as dashed yellow lines. The chemical structures of compounds 1a, 2a, 5a, and 6a [59] are reported.
Figure 7. Molecular docking pose of compound 1a (C atom; violet) (A), 2a (C atom; violet) (B), 4a (C atom; violet) (C), and 5a (C atom; violet) (D) within the hTAAR1 binding site. The docking pose of the reference compound S18616 1 [43] is reported (C atom; light green). The most important residues are shown and labeled. H-bonds are highlighted as dashed yellow lines. The chemical structures of compounds 1a, 2a, 5a, and 6a [59] are reported.
Molecules 29 01739 g007
Figure 8. Chemical structure and hTAAR1 EC50 values of the herein screened agonists 1a10a [59].
Figure 8. Chemical structure and hTAAR1 EC50 values of the herein screened agonists 1a10a [59].
Molecules 29 01739 g008
Table 1. Chemical structure and biological activity as hTAAR1 agonists (EC50 values) of compound 1 and of compounds (Comp.) 221, featuring a flexible linker between the oxazoline ring and the terminal phenyl group [38,43]. The predicted ΔG value of each best-scored protein–ligand complex has been reported, as calculated in terms of molecular docking final scoring function (S, as Kcal/mol).
Table 1. Chemical structure and biological activity as hTAAR1 agonists (EC50 values) of compound 1 and of compounds (Comp.) 221, featuring a flexible linker between the oxazoline ring and the terminal phenyl group [38,43]. The predicted ΔG value of each best-scored protein–ligand complex has been reported, as calculated in terms of molecular docking final scoring function (S, as Kcal/mol).
Comp.Chemical StructureEC50 (nM)S Comp.Chemical StructureEC50 (nM)S
1Molecules 29 01739 i00115−6.854512Molecules 29 01739 i00259−6.1163
2Molecules 29 01739 i003154−6.252513Molecules 29 01739 i004140−5.1443
3Molecules 29 01739 i005330−6.154914Molecules 29 01739 i006230−6.3802
4Molecules 29 01739 i0072900−6.797215Molecules 29 01739 i0081540−4.6935
5Molecules 29 01739 i00918−6.886416Molecules 29 01739 i010730−7.2651
6Molecules 29 01739 i01127−7.132117Molecules 29 01739 i0122260−5.9924
7Molecules 29 01739 i013330−5.615618Molecules 29 01739 i01418−5.8995
8Molecules 29 01739 i015270−6.196319Molecules 29 01739 i01627−5.8673
9Molecules 29 01739 i017580−4.903420Molecules 29 01739 i018360−5.1471
10Molecules 29 01739 i01927−6.593421Molecules 29 01739 i0209−7.1321
11Molecules 29 01739 i02129−6.5403
Table 2. Chemical structure and biological activity as hTAAR1 agonists (EC50 values) of compounds (Comp.) 2237, as phenyl-based derivatives [38,43]. The predicted ΔG value of each best-scored protein–ligand complex has been reported, as calculated in terms of molecular docking final scoring function (S, as Kcal/mol).
Table 2. Chemical structure and biological activity as hTAAR1 agonists (EC50 values) of compounds (Comp.) 2237, as phenyl-based derivatives [38,43]. The predicted ΔG value of each best-scored protein–ligand complex has been reported, as calculated in terms of molecular docking final scoring function (S, as Kcal/mol).
Comp.Chemical StructureEC50 (nM)SComp.Chemical StructureEC50 (nM)S
22Molecules 29 01739 i02267−6.296730Molecules 29 01739 i02341−5.7542
23Molecules 29 01739 i02423−6.637631Molecules 29 01739 i0252670−5.2125
24Molecules 29 01739 i02621−6.766032Molecules 29 01739 i027490−5.6820
25Molecules 29 01739 i028143−6.058733Molecules 29 01739 i02967−6.7540
26Molecules 29 01739 i03031−6.517234Molecules 29 01739 i03111−6.9372
27Molecules 29 01739 i032150−5.924135Molecules 29 01739 i03326−6.4499
28Molecules 29 01739 i034>10,000−5.087636Molecules 29 01739 i03512−6.0151
29Molecules 29 01739 i036165−5.568137Molecules 29 01739 i03717−6.6828
Table 3. Chemical structure and biological activity as hTAAR1 agonists (EC50 values) of the in-house compounds (Comp.) 1a10a [59]. The predicted ΔG value of each best-scored protein–ligand complex has been reported, as calculated in terms of molecular docking final scoring function (S, as Kcal/mol).
Table 3. Chemical structure and biological activity as hTAAR1 agonists (EC50 values) of the in-house compounds (Comp.) 1a10a [59]. The predicted ΔG value of each best-scored protein–ligand complex has been reported, as calculated in terms of molecular docking final scoring function (S, as Kcal/mol).
Comp.Chemical StructureEC50 (nM)SComp.Chemical StructureEC50 (nM)S
1aMolecules 29 01739 i038526−3.16956aMolecules 29 01739 i0391430−2.9174
2aMolecules 29 01739 i040657−3.12547aMolecules 29 01739 i041inactive3.9286
3aMolecules 29 01739 i042756−3.14088aMolecules 29 01739 i043inactive0.4922
4aMolecules 29 01739 i0441060−2.60839aMolecules 29 01739 i045inactive−1.4568
5aMolecules 29 01739 i0461010−3.114110aMolecules 29 01739 i047inactive0.4922
Table 4. Calculated parameters related to Lipinski’s rules, Veber’s rules, and absorption/distribution properties, referred to the compounds (Comp.) 1a6a [59] and to the reference S18616 [43], are reported. The corresponding chemical structure and hTAAR1 EC50 values (EC50) are also shown. Compounds predicted as unable to pass the blood–brain barrier are listed in italics. Reliability index values for a number of descriptors are shown as R.I. Values higher than 0.30 are ranked as reliable by the Advanced Chemistry Development (ACD) Percepta platform (https://www.acdlabs.com/products/percepta-platform/, accessed on 8 April 2024) [70].
Table 4. Calculated parameters related to Lipinski’s rules, Veber’s rules, and absorption/distribution properties, referred to the compounds (Comp.) 1a6a [59] and to the reference S18616 [43], are reported. The corresponding chemical structure and hTAAR1 EC50 values (EC50) are also shown. Compounds predicted as unable to pass the blood–brain barrier are listed in italics. Reliability index values for a number of descriptors are shown as R.I. Values higher than 0.30 are ranked as reliable by the Advanced Chemistry Development (ACD) Percepta platform (https://www.acdlabs.com/products/percepta-platform/, accessed on 8 April 2024) [70].
Comp.MW aN. H-Bond
Acceptor b
N. H-Bond
Donor c
N. Rotatable
Bonds d
CLogP e
(R.I. ≥ 0.45)
TPSA fLogBB gLogPS hHIA (%) iVd (L/kg) j%PPB k
(R.I. ≥ 0.40)
LogKa HAS l
(R.I. ≥ 0.30)
F% m (Oral)
1a297.366120.8751.18−0.58−3.1462.184.094.4827.1
2a311.386131.2151.18−0.51−2.9522.384.544.4930.7
3a359.426132.2151.18−0.53−2.5952.895.794.9580.8
4a366.467140.8254.42−0.463.3174.677.894.328.9
5a382.468140.5263.65−0.48−3.5293.275.844.2315.7
6a310.745022.5348.80−0.77−1.31001.798.255.2937.7
S18616236.703203.1147.610.56−1.41003.271.823.4999.3
Chemical structure
1a
Molecules 29 01739 i048
EC50 = 526 nM
2a
Molecules 29 01739 i049
EC50 = 657 nM
3a
Molecules 29 01739 i050
EC50 = 756 nM
4a
Molecules 29 01739 i051
EC50 = 1060 nM
5a
Molecules 29 01739 i052
EC50 = 1010 nM
6a
Molecules 29 01739 i053
EC50 = 1430 nM
S18616
Molecules 29 01739 i054
EC50 = 15 nM
a Molecular weight; b Number of H-bond acceptors; c Number of H-bond donors; d Number of rotatable bonds; e Logarithmic ratio of the octanol–water partitioning coefficient; f Topological polar surface area; g Extent of brain penetration based on ratio of total drug concentrations in tissue and plasma at steady-state conditions; h Rate of brain penetration. PS represents Permeability-Surface area product and is derived from the kinetic equation of capillary transport; i HIA represents the human intestinal absorption, expressed as percentage of the molecule able to pass through the intestinal membrane; j prediction of Volume of Distribution (Vd) of the compound in the body; k estimation of the plasmatic protein binding event (%PPB); l ligand affinity toward human serum albumin; m Percentage oral bioavailability.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Grossi, G.; Scarano, N.; Musumeci, F.; Tonelli, M.; Kanov, E.; Carbone, A.; Fossa, P.; Gainetdinov, R.R.; Cichero, E.; Schenone, S. Discovery of a Novel Chemo-Type for TAAR1 Agonism via Molecular Modeling. Molecules 2024, 29, 1739. https://doi.org/10.3390/molecules29081739

AMA Style

Grossi G, Scarano N, Musumeci F, Tonelli M, Kanov E, Carbone A, Fossa P, Gainetdinov RR, Cichero E, Schenone S. Discovery of a Novel Chemo-Type for TAAR1 Agonism via Molecular Modeling. Molecules. 2024; 29(8):1739. https://doi.org/10.3390/molecules29081739

Chicago/Turabian Style

Grossi, Giancarlo, Naomi Scarano, Francesca Musumeci, Michele Tonelli, Evgeny Kanov, Anna Carbone, Paola Fossa, Raul R. Gainetdinov, Elena Cichero, and Silvia Schenone. 2024. "Discovery of a Novel Chemo-Type for TAAR1 Agonism via Molecular Modeling" Molecules 29, no. 8: 1739. https://doi.org/10.3390/molecules29081739

APA Style

Grossi, G., Scarano, N., Musumeci, F., Tonelli, M., Kanov, E., Carbone, A., Fossa, P., Gainetdinov, R. R., Cichero, E., & Schenone, S. (2024). Discovery of a Novel Chemo-Type for TAAR1 Agonism via Molecular Modeling. Molecules, 29(8), 1739. https://doi.org/10.3390/molecules29081739

Article Metrics

Back to TopTop