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

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 (1a–10a) as a novel scaffold to target TAAR1. Combined ligand-based (LBCM) and structure based (SBCM) computational methods suggested the biological evaluation of compounds 1a–10a, leading to the identification of derivatives 1a–3a (hTAAR1 EC50 = 526.3–657.4 nM) as promising novel TAAR1 agonists.


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].
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].
Molecules 2024, 29, x FOR PEER REVIEW 2 of 26 mutations in the TAAR1 gene and schizophrenia was reported [15,16], suggesting the implication of TAAR1 in the pathways driving schizophrenia [17].
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].
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 alkylphenyl-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 Structural variations performed at the S18616 1 (A) [43] compound as TAAR1 prototype towards further oxazoline-based agonists 2-37 (B) [38].The promising compounds 5, 11, and 35 have been pointed out (C).The chemical structure of the herein proposed TAAR1 agonists 1a-10a is shown (D).The corresponding hTAAR1 EC 50 values are reported.
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 1-37 (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 1a-10a (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 1a-10a 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 1a-10a pyrimidinone ring is expected to mimic the oxazoline ring of the TAAR1 agonists 1-37.In addition, most of the compounds 1a-10a maintained a key basic moiety in R 1 , to simulate the compound S18616 1 primary amine group.
Herein, in silico investigation of 1a-10a 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 1-37 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 (1a-6a), showing for the 1a derivative (hTAAR1 EC 50 = 526.3nM) the lowest EC 50 value among the in silico evaluated analogues (hTAAR1 EC 50 = 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 1a-6a, 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.

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 2-37 (see the chemical structures in Table S1).Details of the calculated scoring functions are listed in Table S3.Most of the oxazolines 1-37 featured S values spanning from −4.6935 Kcal/mol, such as 15 (hTAAR1 EC 50 = 1540 nM), to −7.1321 Kcal/mol (21; hTAAR1 EC 50 = 9 nM), with the reference compound S18616 being endowed with S = −6.8545Kcal/mol.
As shown in Table 1, a perspective of compound 1 and of compounds 2-21 [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.[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).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 2-37 (see the chemical structures in Table S1).Details of the calculated scoring functions are listed in Table S3.Most of the oxazolines 1-37 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.8545Kcal/mol.
As shown in Table 1, a perspective of compound 1 and of compounds 2-21 [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.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 2-37 (see the chemical structures in Table S1).Details of the calculated scoring functions are listed in Table S3.Most of the oxazolines 1-37 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.8545Kcal/mol.
As shown in Table 1, a perspective of compound 1 and of compounds 2-21 [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.
Table 1.Chemical structure and biological activity as hTAAR1 agonists (EC50 values) of compound 1 and of compounds (Comp.)2-21, 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 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 2-37 (see the chemical structures in Table S1).Details of the calculated scoring functions are listed in Table S3.Most of the oxazolines 1-37 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.8545Kcal/mol.
As shown in Table 1, a perspective of compound 1 and of compounds 2-21 [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.
Table 1.Chemical structure and biological activity as hTAAR1 agonists (EC50 values) of compound 1 and of compounds (Comp.)2-21, 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 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 2-37 (see the chemical structures in Table S1).Details of the calculated scoring functions are listed in Table S3.Most of the oxazolines 1-37 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.8545Kcal/mol.
As shown in Table 1, a perspective of compound 1 and of compounds 2-21 [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.
Table 1.Chemical structure and biological activity as hTAAR1 agonists (EC50 values) of compound 1 and of compounds (Comp.)2-21, 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 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 2-37 (see the chemical structures in Table S1).Details of the calculated scoring functions are listed in Table S3.Most of the oxazolines 1-37 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.8545Kcal/mol.
As shown in Table 1, a perspective of compound 1 and of compounds 2-21 [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.
Table 1.Chemical structure and biological activity as hTAAR1 agonists (EC50 values) of compound 1 and of compounds (Comp.)2-21, 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 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 2-37 (see the chemical structures in Table S1).Details of the calculated scoring functions are listed in Table S3.Most of the oxazolines 1-37 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.8545Kcal/mol.
As shown in Table 1, a perspective of compound 1 and of compounds 2-21 [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.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 2-37 (see the chemical structures in Table S1).Details of the calculated scoring functions are listed in Table S3.Most of the oxazolines 1-37 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.8545Kcal/mol.
As shown in Table 1, a perspective of compound 1 and of compounds 2-21 [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.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 2-37 (see the chemical structures in Table S1).Details of the calculated scoring functions are listed in Table S3.Most of the oxazolines 1-37 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.8545Kcal/mol.
As shown in Table 1, a perspective of compound 1 and of compounds 2-21 [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.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 2-37 (see the chemical structures in Table S1).Details of the calculated scoring functions are listed in Table S3.Most of the oxazolines 1-37 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.8545Kcal/mol.
As shown in Table 1, a perspective of compound 1 and of compounds 2-21 [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.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 2-37 (see the chemical structures in Table S1).Details of the calculated scoring functions are listed in Table S3.Most of the oxazolines 1-37 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.8545Kcal/mol.
As shown in Table 1, a perspective of compound 1 and of compounds 2-21 [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.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 2-37 (see the chemical structures in Table S1).Details of the calculated scoring functions are listed in Table S3.Most of the oxazolines 1-37 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.8545Kcal/mol.
As shown in Table 1, a perspective of compound 1 and of compounds 2-21 [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.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 2-37 (see the chemical structures in Table S1).Details of the calculated scoring functions are listed in Table S3.Most of the oxazolines 1-37 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.8545Kcal/mol.
As shown in Table 1, a perspective of compound 1 and of compounds 2-21 [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.
Table 1.Chemical structure and biological activity as hTAAR1 agonists (EC50 values) of compound 1 and of compounds (Comp.)2-21, 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 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 2-37 (see the chemical structures in Table S1).Details of the calculated scoring functions are listed in Table S3.Most of the oxazolines 1-37 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.8545Kcal/mol.
As shown in Table 1, a perspective of compound 1 and of compounds 2-21 [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.
Table 1.Chemical structure and biological activity as hTAAR1 agonists (EC50 values) of compound 1 and of compounds (Comp.)2-21, 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 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 2-37 (see the chemical structures in Table S1).Details of the calculated scoring functions are listed in Table S3.Most of the oxazolines 1-37 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.8545Kcal/mol.
As shown in Table 1, a perspective of compound 1 and of compounds 2-21 [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.
Table 1.Chemical structure and biological activity as hTAAR1 agonists (EC50 values) of compound 1 and of compounds (Comp.)2-21, 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 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 2-37 (see the chemical structures in Table S1).Details of the calculated scoring functions are listed in Table S3.Most of the oxazolines 1-37 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.8545Kcal/mol.
As shown in Table 1, a perspective of compound 1 and of compounds 2-21 [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.
Table 1.Chemical structure and biological activity as hTAAR1 agonists (EC50 values) of compound 1 and of compounds (Comp.)2-21, 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 2.Chemical structure and biological activity as hTAAR1 agonists (EC50 values) of compounds (Comp.)22-37, 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.)22-37, 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.)22-37, 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.)22-37, 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.)22-37, 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 Regarding the phenyl-based oxazoline compounds 22-37 [38,43], the predicted ∆G value of each best-scored protein-agonist complex is reported in Table 2. Regarding the phenyl-based oxazoline compounds 22-37 [38,43], the predicted ΔG value of each best-scored protein-agonist complex is reported in Table 2.
Table 2.Chemical structure and biological activity as hTAAR1 agonists (EC50 values) of compounds (Comp.)22-37, 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.)22-37, 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.)22-37, 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.)22-37, 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.)22-37, 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.)22-37, 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.)22-37, 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.)22-37, 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 Regarding the phenyl-based oxazoline compounds 22-37 [38,43], the predicted ΔG value of each best-scored protein-agonist complex is reported in Table 2. Regarding the phenyl-based oxazoline compounds 22-37 [38,43], the predicted ΔG value of each best-scored protein-agonist complex is reported in Table 2. Regarding the phenyl-based oxazoline compounds 22-37 [38,43], the predicted ΔG value of each best-scored protein-agonist complex is reported in Table 2. Regarding the phenyl-based oxazoline compounds 22-37 [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 proteinagonist 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 1a-10a 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 1a-10a allowed us to confirm the reliability of the computational studies.

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].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 1a-10a 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 1a-10a allowed us to confirm the reliability of the computational studies.

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].The obtained results are expected to (i) reveal more information about the proteinagonist 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 1a-10a 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 1a-10a allowed us to confirm the reliability of the computational studies.

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].The obtained results are expected to (i) reveal more information about the proteinagonist 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 1a-10a 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 1a-10a allowed us to confirm the reliability of the computational studies.

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].The obtained results are expected to (i) reveal more information about the proteinagonist 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 1a-10a 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 1a-10a allowed us to confirm the reliability of the computational studies.

Molecular Docking Studies of Oxazoline-Based TAAR1 Agonists
According to our molecular docking calculations, the compound 1 (S18616, EC 50 = 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 Figures 2B and 2C,D, respectively.
Herein, compound 1 (S18616) was taken as a reference compound for a three-dimensional SAR analysis involving the other analogues 2-37, based on compound 1's high potency as TAAR1 agonist and its conformationally locked binding pose.
The TAAR1 ligands 2-37 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.
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 3-4 (hTAAR1 EC 50 = 330-2900 nM) and 2 (hTAAR1 EC 50 = 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).

Molecular Docking Studies of Oxazoline-Based TAAR1 Agonists
According to our molecular docking calculations, the compound 1 (S18616, E nM) docking pose revealed an H-bond between the amino group of the ligand side chain of residue ASP103 (Figure 2A).This is probably an interaction key for b ity, as was also hypothesized for other TAAR1-targeting molecules [66].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 EC 50 = 18 nM), 6 (hTAAR1 EC 50 = 27 nM), and 19 (hTAAR1 EC 50 = 27 nM), respectively.This modification was favorable for activity, leading to a decrease in EC 50 values of around one order of magnitude with respect to compounds 2 and 3, being quite comparable with the prototype 1 (hTAAR1 EC 50 = 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 EC 50 = 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 pivotal H-bond with the key residue ASP103, in tandem with additional π-π s with the surrounding aromatic residues.Indeed, a favorable interaction between t nyl ring and PHE268 (T-shaped π-interaction) was highlighted, together with the in of the phenyl in a sub-pocket defined by SER198, SER108, SER 197, and PHE268.other hand, no hydrophobic contacts involving PHE186 were present.[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.
The introduction of additional alkyl groups at the stereo-center or involving the ethylchain spacer was not advantageous, as shown by compounds 7 (EC 50 = 330 nM) and 15-17 (hTAAR1 EC 50 = 730-2260 nM), except for 18 (hTAAR1 EC 50 = 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 EC 50 = 2260 nM) was differently folded (Figure 3D), exhibiting higher EC 50 values than 18 (hTAAR1 EC 50 = 18 nM).
As regards the phenyl heteroalkyl-containing compounds (8-14, 20-21; hTAAR1 EC 50 = 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 EC 50 values (hTAAR1 EC 50 = 580 nM) than the congeners 10-11 (hTAAR1 EC 50 = 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).[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.
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 EC 50 = 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 EC 50 = 29 nM) and 12 (hTAAR1 EC 50 = 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 EC 50 = 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 (22-37; hTAAR1 EC 50 = 11-10.000nM), as reported for 22-24, 26, 30, 33-37 (EC 50 = 21-67 nM).The docking pose of compound 23 (hTAAR1 EC 50 = 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 EC 50 = 21 nM) and para-(25, hTAAR1 EC 50 = 143 nM) chloro substitution were explored in addition to the previously cited orthoanalogue 23.While the meta-substitution was well tolerated, the para-substitution led to a considerable increase in EC 50 .Then, para-bromo substitution (27; hTAAR1 EC 50 = 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 EC 50 = 10,000 nM).Interestingly, dichloro substitution at the metaand parapositions (26; hTAAR1 EC 50 = 31 nM) led to a good activity.Compound 32 (hTAAR1 EC 50 = 490 nM), with a fluorine in metaposition, exhibited a strong decrease in activity compared to the meta-chloro-containing analogue 24 (hTAAR1 EC 50 = 21 nM).However, if this modification is accompanied by the presence of a methyl substituent in orthoposition, the activity is comparable to reference compound 1 (37; hTAAR1 EC 50 = 17 nM).
Accordingly, further to the halogen atom, the additional introduction of electrondonor substituents endowed with hydrophobic properties, such as aliphatic chains, led to improved TAAR1 agonists, as reported for 34-36 (hTAAR1 EC 50 = 11-26 nM).Based on the molecular docking pose of the very potent compound 34 (hTAAR1 EC 50 = 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 EC 50 = 26 nM) and 36 (hTAAR1 EC 50 = 12 nM).Bulkier substituents such as the diaryl moiety were not tolerated, as described for compound 31 (hTAAR1 EC 50 = 2670 nM).Finally, the only insertion of the ortho-methyl substitution led to quite effective EC 50 values [see 33 (hTAAR1 EC 50 = 67 nM)].

In Silico Evaluation of Pyrimidinone-Benzimidazole-Based Compounds as TAAR1 Agonists
According to the aforementioned 1-37 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 1-37 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 1a-10a, 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.
Molecules 2024, 29, x FOR PEER REVIEW 10 of 26 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)].

In Silico Evaluation of Pyrimidinone-Benzimidazole-Based Compounds as TAAR1 Agonists
According to the aforementioned 1-37 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 1-37 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 1a-10a, 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 1a-10a series was endowed with additional nitrogen atoms if compared to the S18616 1 ring B. The previously reported oxazoline series 2-37 pointed out the effectiveness of proper heteroatoms within Figure 5. Schematic comparison of the herein investigated series of compounds 1a-10a [59], with respect to the hTAAR1 prototype S18616 1 [43].
While the pyrimidinone ring A featured comparable electrostatic properties and electron-rich atoms to the corresponding 1 oxazoline ring (ring A), the terminal sixmembered 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 1a-10a series was endowed with additional nitrogen atoms if compared to the S18616 1 ring B. The previously reported oxazoline series 2-37 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 R 1 , 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 R 1 have been taken into account.
Finally, the steric hindrance displayed by 1a-10a 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.
As regard the pyrimidinone-benzimidazoles 1a-10a, three derivatives (1a-3a) were endowed with Glob-Prod values higher than 0.3000, while the analogues 4a-6a spanned from 0.2700 to 0.3000.Compounds 7a-10a displayed the lowest values of Glob-Prod in this series of compounds (Glob-Prod = 0.1051-0.2692).Then, we proceed with a 3Danalysis 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.
The O MIF (red) of 1a is well matching the template one in the area of the piperazineoxazoline moieties, as the protonated piperazine ring of this compound mimics the aminogroup 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.
with the H-bond donor (HBD) probe, (iii) O to assess favorable interaction with the Hbond acceptor (HBA) probe.
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.
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 1a-10a 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 1-37 agonists.
Molecular docking studies of compounds 1a-10a 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-benzimidazolebased compounds are listed in Table S6.
As shown in Table 3, a perspective of the chemical structure of 1a-10a [59] and of the predicted ∆G value of each best-scored protein-ligand complex is listed.The corresponding hTAAR1 EC 50 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 1a-3a, exhibiting S values < −3.0000 Kcal/mol, while 4a-6a spanned from −3.1141 to −2.6083 Kcal/mol.Conversely, the analogues 7a-10a were predictive as poorly active/inactive exhibiting S values > −1.4568 Kcal/mol.On this basis, the docking modes of only 1a-6a are discussed as follows.
According to our calculations, the area of interaction of most of the aforementioned 1a-6a is the same of the reference one S18616 1.
The choice of a piperazine group in R 1 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.Table 3.Chemical structure and biological activity as hTAAR1 agonists (EC 50 values) of the in-house compounds (Comp.)1a-10a [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).in the experimental section.The obtained scoring functions for the pyrimidinone-benzimidazole-based compounds are listed in Table S6.

Comp
As shown in Table 3, a perspective of the chemical structure of 1a-10a [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.Table 3.Chemical structure and biological activity as hTAAR1 agonists (EC50 values) of the in-house compounds (Comp.)1a-10a [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).S6.

Comp
As shown in Table 3, a perspective of the chemical structure of 1a-10a [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.Table 3.Chemical structure and biological activity as hTAAR1 agonists (EC50 values) of the in-house compounds (Comp.)1a-10a [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).S6.

Comp
As shown in Table 3, a perspective of the chemical structure of 1a-10a [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.Table 3.Chemical structure and biological activity as hTAAR1 agonists (EC50 values) of the in-house compounds (Comp.)1a-10a [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).S6.

Comp
As shown in Table 3, a perspective of the chemical structure of 1a-10a [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.Table 3.Chemical structure and biological activity as hTAAR1 agonists (EC50 values) of the in-house compounds (Comp.)1a-10a [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).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 1a-3a, exhibiting S values < −3.0000 Kcal/mol, while 4a-6a spanned from −3.1141 to −2.6083 Kcal/mol.Conversely, the analogues 7a-10a were predictive as poorly active/inactive exhibiting S values > −1.4568 Kcal/mol.On this basis, the docking modes of only 1a-6a are discussed as follows.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 1a-3a, exhibiting S values < −3.0000 Kcal/mol, while 4a-6a spanned from −3.1141 to −2.6083 Kcal/mol.Conversely, the analogues 7a-10a were predictive as poorly active/inactive exhibiting S values > −1.4568 Kcal/mol.On this basis, the docking modes of only 1a-6a are discussed as follows.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 1a-3a, exhibiting S values < −3.0000 Kcal/mol, while 4a-6a spanned from −3.1141 to −2.6083 Kcal/mol.Conversely, the analogues 7a-10a were predictive as poorly active/inactive exhibiting S values > −1.4568 Kcal/mol.On this basis, the docking modes of only 1a-6a are discussed as follows.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 1a-3a, exhibiting S values < −3.0000 Kcal/mol, while 4a-6a spanned from −3.1141 to −2.6083 Kcal/mol.Conversely, the analogues 7a-10a were predictive as poorly active/inactive exhibiting S values > −1.4568 Kcal/mol.On this basis, the docking modes of only 1a-6a are discussed as follows.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 1a-3a, exhibiting S values < −3.0000 Kcal/mol, while 4a-6a spanned from −3.1141 to −2.6083 Kcal/mol.Conversely, the analogues 7a-10a were predictive as poorly active/inactive exhibiting S values > −1.4568 Kcal/mol.On this basis, the docking modes of only 1a-6a are discussed as follows.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 1a-3a, exhibiting S values < −3.0000 Kcal/mol, while 4a-6a spanned from −3.1141 to −2.6083 Kcal/mol.Conversely, the analogues 7a-10a were predictive as poorly active/inactive exhibiting S values > −1.4568 Kcal/mol.On this basis, the docking modes of only 1a-6a are discussed as follows.
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.
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.
According to our calculations, the area of interaction of most of the aforementioned 1a-6a is the same of the reference one S18616 1.
The choice of a piperazine group in R 1 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.
As regards compound 4a, the piperazine moiety established one H-bond with ASP103 (Figure 7C).
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.

Biological Evaluation of Compounds 1a-10a 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 1a-10a.While compounds 1a-6a proved to be the most promising in silico, based on the previously cited computational studies, in vitro screening of the derivatives 7a-10a has been also performed to elucidate the SAR of the series.Thus, all the novel compounds 1a-10a 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].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 1a-6a, the piperazine-containing derivatives 1a-5a (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 1a-2a (hTAAR1 EC50 = 526-657 nM), endowed with comparable or higher potency values with respect to the 18% of the oxazoline derivatives 1-37 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 R 1 led to compounds 6a-10a, most of them (7a-10a) 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 R 1 .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].
Table 4. Calculated parameters related to Lipinski's rules, Veber's rules, and absorption/distribution properties, referred to the compounds (Comp.)1a-6a [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 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 EC 50 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 EC 50 values.All the compounds displayed an Emax > 85% at m/hTAAR1, and thus they were regarded as full agonists.
Among the active compounds 1a-6a, the piperazine-containing derivatives 1a-5a (hTAAR1 EC 50 = 525-1060 nM) turned out to be more potent than 6a (hTAAR1 EC 50 = 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 1a-2a (hTAAR1 EC 50 = 526-657 nM), endowed with comparable or higher potency values with respect to the 18% of the oxazoline derivatives 1-37 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 R 1 led to compounds 6a-10a, most of them (7a-10a) 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 R 1 .Accordingly, compound 3a (hTAAR1 EC 50 = 756 nM) displayed lower hTAAR1 EC 50 values than 6a (hTAAR1 EC 50 = 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 EC 50 = 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].
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 1a-3a 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 1a-3a (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].

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].
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 oor 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 1a-10a [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 1a-10a mimics the main core of the reference compound S18616, compounds 1a-3a being predicted as the most promising as putative hTAAR1 ligands.In particular, the choice of a piperazine group in R 1 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 1a-3a (hTAAR1 EC 50 = 526.3-756.80nM) 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 1a-3a MIFs were in good agreement with those featured by the template S18616, thanks to the corresponding basic moieties and aromatic features.While the analogues 4a-6a partially mimic the template MIFs due to the amine groups, the lowest-scored congeners 7a-10a poorly simulated the reference MIFs.Accordingly, the following biological evaluation revealed 4a-6a (hTAAR1 = 1010-1430 nM) as modest TAAR1 agonists.A perspective of this information pointed out the pyrimidinonebenzimidazole 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 1a-3a as druglike derivatives.Among them, compound 3a was endowed with the highest predicted LD 50 values and oral bioavailability parameters.

Computational Studies
All the studied compounds, including the oxazoline 1-37 and the pyrimidinone benzimidazoles 1a-10a, 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].

Chemicals
Compounds 1a-10a were synthesized as previously described [46].The synthetic route to compounds 1a-5a and 7a, 10a included the reaction of 2-aminobenzimidazole with proper alkyl halides in refluxing acetone in the presence of anhydrous K 2 CO 3 / 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.

Biological Evaluation of Compounds 1a-10a 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 EC 50 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 EC 50 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].

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 1-37 have been investigated via molecular docking studies, prior to proceeding with the biological evaluation of a series of in-house pyrimidinonebenzimidazole derivatives (1a-10a).Based on the combined LB-and SBCM results, 1a-10a have been assayed as putative TAAR1 agonists.Final biological evaluation revealed the piperazine-containing derivatives 1a-3a (hTAAR1 EC 50 = 526.3-657.4nM) 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.

Figure 2 .
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

Figure 2 .
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.01software [67].

Figure 2 .
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

Figure 2 .
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

Figure 2 .
Figure 2. Molecular docking pose of compound S18616 1 [43] (C atom; light green) wi hTAAR1 binding site (A).The most important residues are shown and labeled.H-bonds a lighted as dashed yellow lines.The TAAR1 agonist chemical structure and the ligand hydr city/electrostatic properties are also reported (hydrophobic and mild polar areas are shown and cyan, respectively; electron-rich and -poor groups are highlighted in red and blue).(B).view and side view of the receptor in the presence of S18616 1 are shown in (C,D), respectiv H-bonding and hydrophobic areas of the protein are colored magenta and green.All the representations have been produced via MOE2019.01software [67].

Figure 2 .
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.01software [67].
couraged.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 (Figure4A).
Molecules 2024, 29, x FOR PEER REVIEW 13 of 26 in the experimental section.The obtained scoring functions for the pyrimidinone-benzimidazole-based compounds are listed in Table

2a
Molecules 2024, 29, x FOR PEER REVIEW 13 of 26 in the experimental section.The obtained scoring functions for the pyrimidinone-benzimidazole-based compounds are listed in Table Molecules 2024, 29, x FOR PEER REVIEW 13 of 26 in the experimental section.The obtained scoring functions for the pyrimidinone-benzimidazole-based compounds are listed in Table

Table 1 .
Chemical structure and biological activity as hTAAR1 agonists (EC 50 values) of compound 1 and of compounds (Comp.)2-21, featuring a flexible linker between the oxazoline ring and the terminal phenyl group