Design and Synthesis of Arylpiperazine Serotonergic/Dopaminergic Ligands with Neuroprotective Properties

Long-chain arylpiperazine scaffold is a versatile template to design central nervous system (CNS) drugs that target serotonin and dopamine receptors. Here we describe the synthesis and biological evaluation of ten new arylpiperazine derivatives designed to obtain an affinity profile at serotonin 5-HT1A, 5-HT2A, 5-HT7 receptor, and dopamine D2 receptor of prospective drugs to treat the core symptoms of autism spectrum disorder (ASD) or psychosis. Besides the structural features required for affinity at the target receptors, the new compounds incorporated structural fragments with antioxidant properties to counteract oxidative stress connected with ASD and psychosis. All the new compounds showed CNS MultiParameter Optimization score predictive of desirable ADMET properties and cross the blood–brain barrier. We identified compound 12a that combines an affinity profile compatible with antipsychotic activity (5-HT1A Ki = 41.5 nM, 5-HT2A Ki = 315 nM, 5-HT7 Ki = 42.5 nM, D2 Ki = 300 nM), and compound 9b that has an affinity profile consistent with studies in the context of ASD (5-HT1A Ki = 23.9 nM, 5-HT2A Ki = 39.4 nM, 5-HT7 Ki = 45.0 nM). Both compounds also had antioxidant properties. All compounds showed low in vitro metabolic stability, the only exception being compound 9b, which might be suitable for studies in vivo.


Introduction
The long-chain arylpiperazine structure is a versatile model that has allowed the identification of several drugs over the years, such as antipsychotics (aripiprazole, ziprasidone, and lurasidone) [1], anxiolytics (buspirone, tandospirone) [2], and antiparkinsonian drugs (piribedil) [3] (Figure 1). This variety of actions is due to the possibility of modulating the pharmacological profile of the arylpiperazine derivatives through structural variations. Even if more than 30 years have passed since the introduction of buspirone in the market, the long-chain arylpiperazine model is still valid, so that, in the last five years, new drugs with this structure have been introduced on the market [4].
The above-mentioned drugs owe their action to the modulation at the central nervous system (CNS) level of one or more serotonin and/or dopamine receptors that are relevant to the specific pathology. For example, aripiprazole owes its antipsychotic action to the blocking of dopaminergic D 2 and serotonergic 5-HT 2A receptors. The anxiolytic action of buspirone is due to a partial agonist activity on the 5-HT 1A receptors, while the antiparkinsonian effect of piribedil is due to dopaminergic D 2 receptors activation.
Recently, we reported on a set of arylpiperazine derivatives with an activity profile on a subgroup of serotonin receptors (5-HT 7 and 5-HT 1A agonist and 5-HT 2A antagonist) specifically designed to obtain prospective drugs for the treatment of the core symptoms of autism spectrum disorder [5]. Over the last few years, several studies have highlighted that besides an adequate modulation of the neuroreceptors relevant in the pathology, a neuroprotective action due to a reduction of oxidative stress might have a beneficial effect on these pathologies. In fact, schizophrenia, depression, and autism spectrum disorder (ASD) have a significant neuroinflammatory component linked to oxidative stress [6]. Therefore, the advantage of intervening in targeting a set of receptors and oxidative stress is evident. Of note, recent studies have shown that the activation of the 5-HT 1A [7] and 5-HT 7 receptors [8] produces neuroprotective effects. Thus, for the present study, we have designed a set of new arylpiperazine derivatives with the structural characteristics to target a group of neuroreceptors relevant to treating schizophrenia or ASD and producing neuroprotective effects. reported on a set of arylpiperazine derivatives with an activity profile on a subgroup o serotonin receptors (5-HT7 and 5-HT1A agonist and 5-HT2A antagonist) specifically de signed to obtain prospective drugs for the treatment of the core symptoms of autism spec trum disorder [5]. Over the last few years, several studies have highlighted that beside an adequate modulation of the neuroreceptors relevant in the pathology, a neuroprotec tive action due to a reduction of oxidative stress might have a beneficial effect on thes pathologies. In fact, schizophrenia, depression, and autism spectrum disorder (ASD) hav a significant neuroinflammatory component linked to oxidative stress [6]. Therefore, th advantage of intervening in targeting a set of receptors and oxidative stress is evident. O note, recent studies have shown that the activation of the 5-HT1A [7] and 5-HT7 receptor [8] produces neuroprotective effects. Thus, for the present study, we have designed a se of new arylpiperazine derivatives with the structural characteristics to target a group o neuroreceptors relevant to treating schizophrenia or ASD and producing neuroprotectiv effects.

Study Design
The starting point for our study was AG-45 (compound 20b in Ref. [5], Figure 1 which has an affinity for 5-HT1A, 5-HT2A, and 5-HT7 receptors. We specifically searche this affinity profile to obtain prospective drugs for the treatment of the core symptoms o ASD. Thus, we introduced structural fragments known for their antioxidant propertie such as the pyridyl [9], 1,2,4-triazine [10], 1,4-benzoxazine [11], and benzopyrone [12], an left unchanged the features responsible for the affinity profile (i.e., the biphenyl-like struc ture linked to the piperazine ring, and a trimethylene or tetramethylene spacer betwee the piperazine ring and the terminal fragment). The design of the target compounds graphically illustrated in Figure 2. Interestingly, AG-45 showed a measurable affinity fo the D2 receptor, and this prompted us to evaluate structural modifications to increase D receptor affinity in the prospect of obtaining antipsychotic compounds with antioxidan properties. Of note, while the affinity for D2 and 5-HT2A receptors is a crucial feature o atypical antipsychotics, targeting 5-HT1A and 5-HT7 receptors have been variously corre lated to beneficial effects in the context of psychosis [13,14].
The discovery of CNS drugs also faces the challenge of controlling physicochemica properties. CNS Multi-Parameter Optimisation (MPO) scoring tool is an approach tha can facilitate the design of CNS drugs [15]. This tool gives a desirability score (0.05-1) fo six physicochemical properties: molecular weight (MW); calculated partition coefficien (cLogP); calculated distribution coefficient at pH = 7.4 (cLogD7.4); acid dissociation con stant (pKa); topological polar surface area (TPSA); number of hydrogen bond donor

Study Design
The starting point for our study was AG-45 (compound 20b in Ref. [5], Figure 1), which has an affinity for 5-HT 1A , 5-HT 2A , and 5-HT 7 receptors. We specifically searched this affinity profile to obtain prospective drugs for the treatment of the core symptoms of ASD. Thus, we introduced structural fragments known for their antioxidant properties, such as the pyridyl [9], 1,2,4-triazine [10], 1,4-benzoxazine [11], and benzopyrone [12], and left unchanged the features responsible for the affinity profile (i.e., the biphenyl-like structure linked to the piperazine ring, and a trimethylene or tetramethylene spacer between the piperazine ring and the terminal fragment). The design of the target compounds is graphically illustrated in Figure 2. Interestingly, AG-45 showed a measurable affinity for the D 2 receptor, and this prompted us to evaluate structural modifications to increase D 2 receptor affinity in the prospect of obtaining antipsychotic compounds with antioxidant properties. Of note, while the affinity for D 2 and 5-HT 2A receptors is a crucial feature of atypical antipsychotics, targeting 5-HT 1A and 5-HT 7 receptors have been variously correlated to beneficial effects in the context of psychosis [13,14].
The discovery of CNS drugs also faces the challenge of controlling physicochemical properties. CNS Multi-Parameter Optimisation (MPO) scoring tool is an approach that can facilitate the design of CNS drugs [15]. This tool gives a desirability score (0.05-1) for six physicochemical properties: molecular weight (MW); calculated partition coefficient (cLogP); calculated distribution coefficient at pH = 7.4 (cLogD 7.4 ); acid dissociation constant (pK a ); topological polar surface area (TPSA); number of hydrogen bond donors (HBD). The sum of these scores yields the CNS MPO score on a 0.3-6 scale. A CNS MPO desirability score higher than 4 is predictive of desirable ADMET properties and cross the blood-brain barrier [15]. As shown in Table 1, all the newly designed compounds display favorable characteristics for a prospective CNS drug. (HBD). The sum of these scores yields the CNS MPO score on a 0.3-6 scale. A CNS MPO desirability score higher than 4 is predictive of desirable ADMET properties and cross the blood-brain barrier [15]. As shown in Table 1, all the newly designed compounds display favorable characteristics for a prospective CNS drug.

Chemistry
The synthesis of the target compounds 8a,b-12a,b was accomplished through the reaction of the appropriate chloroalkyl derivatives 2a,b-4a,b with arylpiperazines 5-7 (Scheme 1). Both chloroalkyl and arylpiperazine intermediates were synthesized with literature methods (as detailed in the Materials and Methods) except alkyl chloride 2b which was prepared by alkylating 6-hydroxy-2-methyl-2H-benzo[b] [1,4]oxazin-3(4H)-one (1) [16] with 1-bromo-4-chlorobutane via the Williamson's synthesis under basic conditions.  (HBD). The sum of these scores yields the CNS MPO score on a 0.3-6 scale. A CNS MPO desirability score higher than 4 is predictive of desirable ADMET properties and cross the blood-brain barrier [15]. As shown in Table 1, all the newly designed compounds display favorable characteristics for a prospective CNS drug.

Chemistry
The synthesis of the target compounds 8a,b-12a,b was accomplished through the reaction of the appropriate chloroalkyl derivatives 2a,b-4a,b with arylpiperazines 5-7 (Scheme 1). Both chloroalkyl and arylpiperazine intermediates were synthesized with literature methods (as detailed in the Materials and Methods) except alkyl chloride 2b which was prepared by alkylating 6-hydroxy-2-methyl-2H-benzo[b] [1,4] (HBD). The sum of these scores yields the CNS MPO score on a 0.3-6 scale. A CNS MPO desirability score higher than 4 is predictive of desirable ADMET properties and cross the blood-brain barrier [15]. As shown in Table 1, all the newly designed compounds display favorable characteristics for a prospective CNS drug.

Chemistry
The synthesis of the target compounds 8a,b-12a,b was accomplished through the reaction of the appropriate chloroalkyl derivatives 2a,b-4a,b with arylpiperazines 5-7 (Scheme 1). Both chloroalkyl and arylpiperazine intermediates were synthesized with literature methods (as detailed in the Materials and Methods) except alkyl chloride 2b which was prepared by alkylating 6-hydroxy-2-methyl-2H-benzo[b] [1,4] (HBD). The sum of these scores yields the CNS MPO score on a 0.3-6 scale. A CNS MPO desirability score higher than 4 is predictive of desirable ADMET properties and cross the blood-brain barrier [15]. As shown in Table 1, all the newly designed compounds display favorable characteristics for a prospective CNS drug.

Chemistry
The synthesis of the target compounds 8a,b-12a,b was accomplished through the reaction of the appropriate chloroalkyl derivatives 2a,b-4a,b with arylpiperazines 5-7 (Scheme 1). Both chloroalkyl and arylpiperazine intermediates were synthesized with literature methods (as detailed in the Materials and Methods) except alkyl chloride 2b which was prepared by alkylating 6-hydroxy-2-methyl-2H-benzo[b] [1,4] (HBD). The sum of these scores yields the CNS MPO score on a 0.3-6 scale. A CNS MPO desirability score higher than 4 is predictive of desirable ADMET properties and cross the blood-brain barrier [15]. As shown in Table 1, all the newly designed compounds display favorable characteristics for a prospective CNS drug.

Chemistry
The synthesis of the target compounds 8a,b-12a,b was accomplished through the reaction of the appropriate chloroalkyl derivatives 2a,b-4a,b with arylpiperazines 5-7 (Scheme 1). Both chloroalkyl and arylpiperazine intermediates were synthesized with literature methods (as detailed in the Materials and Methods) except alkyl chloride 2b which was prepared by alkylating 6-hydroxy-2-methyl-2H-benzo[b] [1,4] (HBD). The sum of these scores yields the CNS MPO score on a 0.3-6 scale. A CNS MPO desirability score higher than 4 is predictive of desirable ADMET properties and cross the blood-brain barrier [15]. As shown in Table 1, all the newly designed compounds display favorable characteristics for a prospective CNS drug.

Chemistry
The synthesis of the target compounds 8a,b-12a,b was accomplished through the reaction of the appropriate chloroalkyl derivatives 2a,b-4a,b with arylpiperazines 5-7 (Scheme 1). Both chloroalkyl and arylpiperazine intermediates were synthesized with literature methods (as detailed in the Materials and Methods) except alkyl chloride 2b which was prepared by alkylating 6-hydroxy

Chemistry
The synthesis of the target compounds 8a,b-12a,b was accomplished through the reaction of the appropriate chloroalkyl derivatives 2a,b-4a,b with arylpiperazines 5-7 (Scheme 1). Both chloroalkyl and arylpiperazine intermediates were synthesized with literature methods (as detailed in the Materials and Methods) except alkyl chloride 2b which was prepared by alkylating 6-hydroxy-2-methyl-2H-benzo[b] [1,4]

Radioligand Binding Experiments
The final compounds were tested for their affinity at serotonin 5-HT1A, 5-HT2A, and 5-HT7 receptors and at dopamine D1, D2, D3, and D5 receptors. As a first step, competition radioligand binding assays were carried out for all the compounds at the single concentration of 10 µ M, on membrane preparations from stable cell lines heterologously expressing the human cloned 5-HT1A, 5-HT2A, or 5-HT7 receptors. Affinity was determined for those compounds showing % of displacement of the specific binding of the radioligand over 65% at 10 µ M concentration, by means of competition binding curves of six concentration datapoints (typically from 0.1 nM/1 nM to 10 µ M/100 µ M). The affinity of the compounds, expressed as equilibrium dissociation constant (Ki) value, is shown in Table 2. For the affinity determination at dopamine receptor subtypes, membrane fractions of human dopamine D2s, D3, D1, and D5 receptor-expressing cells were used. They were co-incubated with [ 3 H]spiperone (for D2s and D3 receptor) or [ 3 H]SCH23390 (for D1 and D5 receptor) and appropriate concentrations of test compound for 120 min. Bound radioligand was separated from the free ligand using filtration. Radioactivity was measured by liquid scintillation counting. The resulting Ki values with 95% confidence intervals are listed in Table 2.

Radioligand Binding Experiments
The final compounds were tested for their affinity at serotonin 5-HT 1A , 5-HT 2A , and 5-HT 7 receptors and at dopamine D 1 , D 2 , D 3 , and D 5 receptors. As a first step, competition radioligand binding assays were carried out for all the compounds at the single concentration of 10 µM, on membrane preparations from stable cell lines heterologously expressing the human cloned 5-HT 1A , 5-HT 2A , or 5-HT 7 receptors. Affinity was determined for those compounds showing % of displacement of the specific binding of the radioligand over 65% at 10 µM concentration, by means of competition binding curves of six concentration datapoints (typically from 0.1 nM/1 nM to 10 µM/100 µM). The affinity of the compounds, expressed as equilibrium dissociation constant (K i ) value, is shown in Table 2. For the affinity determination at dopamine receptor subtypes, membrane fractions of human dopamine D 2s , D 3 , D 1 , and D 5 receptor-expressing cells were used. They were co-incubated with [ 3 H]spiperone (for D 2s and D 3 receptor) or [ 3 H]SCH23390 (for D 1 and D 5 receptor) and appropriate concentrations of test compound for 120 min. Bound radioligand was separated from the free ligand using filtration. Radioactivity was measured by liquid scintillation counting. The resulting K i values with 95% confidence intervals are listed in Table 2. 47.9 ± 3.5 1 Average data from n = 3 independent experiments performed in duplicate. 2 Full displacement of specific binding not achieved at the maximum concentration assayed (100 µM), so K i value might not be accurately estimated; maximum displacement achieved was 84%).

Docking Studies
Docking studies were performed to rationalize the different binding profiles of the compounds. Representative compounds were docked at the serotonin 5-HT 1A and 5-HT 2A , and dopamine D 2 and D1 using the available crystal structures of the receptors (Figures 3 and 4).

Docking Studies
Docking studies were performed to rationalize the different binding profiles of the compounds. Representative compounds were docked at the serotonin 5-HT1A and 5-HT2A, and dopamine D2 and D1 using the available crystal structures of the receptors (Figures 3  and 4).

In Vitro Metabolic Stability
To evaluate the susceptibility to first-pass oxidative metabolism, the leading cause of metabolic degradation in vivo [17] the target compounds were incubated with rat liver microsomes. The turnover of the parent compound was assessed as the percentage of the parent compound recovered after 30 min of incubation with microsomes in the presence of an NADPH-regenerating system [18]. The percentages of the recovered unchanged parent compound are listed in Table 1.

Neuroprotection against H2O2 in SHSY-5Y Cell Line
Initially, the cytotoxicity of the target compounds was assessed by treating SH-SY5Y cells for 48 h with various concentrations of compounds (0.1 to 100 μM). Compounds 11a and 11b were not cytotoxic at the tested doses, while the remaining were cytotoxic, exhibiting EC50 values between 20-50 μM ( Table 2). Cells treatment with 1 μM or 5 μM of each compound for 24 h did not result in an effect on cell viability (<4% and <8% decrease at a concentration of 1 μM and 5 μM, respectively).
Treatment of SH-SY5Y cells with 400 μM H2O2 for 24 h led to a reduction of cell viability of 72-76% ( Figure 3). The protective effect of each compound was determined by pretreating the cells with the test compound for 3 h at the non-toxic doses of 1 μM and 5 μM, and subsequent treatment with 400 μM H2O2 for an additional 24 h. Pre-treatment of SH-SY5Y cells with compounds 8a,b, 9b, 10a, 11a, and 12a resulted in a 14-20% recovery of cell viability with respect to treatment with the H2O2 alone ( Figure 5).

In Vitro Metabolic Stability
To evaluate the susceptibility to first-pass oxidative metabolism, the leading cause of metabolic degradation in vivo [17] the target compounds were incubated with rat liver microsomes. The turnover of the parent compound was assessed as the percentage of the parent compound recovered after 30 min of incubation with microsomes in the presence of an NADPH-regenerating system [18]. The percentages of the recovered unchanged parent compound are listed in Table 1.

Neuroprotection against H 2 O 2 in SHSY-5Y Cell Line
Initially, the cytotoxicity of the target compounds was assessed by treating SH-SY5Y cells for 48 h with various concentrations of compounds (0.1 to 100 µM). Compounds 11a and 11b were not cytotoxic at the tested doses, while the remaining were cytotoxic, exhibiting EC 50 values between 20-50 µM ( Table 2). Cells treatment with 1 µM or 5 µM of each compound for 24 h did not result in an effect on cell viability (<4% and <8% decrease at a concentration of 1 µM and 5 µM, respectively).
Treatment of SH-SY5Y cells with 400 µM H 2 O 2 for 24 h led to a reduction of cell viability of 72-76% ( Figure 3). The protective effect of each compound was determined by pretreating the cells with the test compound for 3 h at the non-toxic doses of 1 µM and 5 µM, and subsequent treatment with 400 µM H 2 O 2 for an additional 24 h. Pre-treatment of SH-SY5Y cells with compounds 8a,b, 9b, 10a, 11a, and 12a resulted in a 14-20% recovery of cell viability with respect to treatment with the H 2 O 2 alone ( Figure 5).

Figure 5.
Cell viability of SH-SY5Y cells pretreated with 1 μM or 5 μM of the test compound for 3 h, and subsequently exposed to 400 μM H2O2 for 24 h. The data are shown as the mean ± S.E.M. of three independent experiments (** p < 0.001 * p < 0.05 vs. 400 μM H2O2-treated group).

Discussion
The binding affinities of the target compounds are listed in Table 2. The analysis of the binding data at the screened receptors provides information regarding the role of the three structural elements that compose the structure, i.e., biphenyl-like system, spacer, and terminal fragment.
Considering the affinities for the 5-HT1A receptor, we can note that both biphenyllike systems and terminal fragments selected for the study are compatible with high affinity at the receptor. In fact, four compounds out of 10 show Ki values lower than 17.6 nM (compounds 8a,b, 11b, 12b). The differences in affinity are related to the length of the spacer between the piperazine and the terminal fragment. Examining the compounds characterized by the oxotriazinyl terminal fragment (compounds 9a,b-11a,b), the derivatives with a tetramethylene spacer have higher affinity than the trimethylene counterparts. The influence of the spacer length on affinity is related to the nature of the terminal fragment. In fact, considering the compounds featuring the bipyridyl system linked to the piperazine, the Ki ratios of the trimethylene and tetramethylene homologs are variable (11a/11b = 52; 12a/12b = 3.7; 8a/8b = 0.7).
As for the 5-HT2A receptor, the compounds featuring the bipyridyl system linked to the piperazine display moderate to low affinity. Within this group of compounds, the terminal fragment plays a role on affinity. In fact, the compounds with the oxotriazinyl fragment have a low affinity, whereas a bicyclic terminal fragment imparts comparatively higher affinity. The replacement of the bipyridyl with a biphenyl system has a favorable effect on the affinity. In fact, for the compounds with a tetramethylene spacer, 9b has a slightly higher affinity than 10b which shows much higher affinity than 11b. A similar trend can be observed for the trimethylene spacer derivatives 9a, 10a, and 11a.
The affinity values of the target compounds at 5-HT7 receptor are distributed in a narrower range as compared to the other receptors under study. In fact, the highest affinity values were shown by 8b (Ki = 9.38 nM) and the lowest by 11b (Ki = 79.4 nM). Consequently, the differences in affinity between homologs (trimethylene or tetramethylene spacer derivatives) or analogs (i.e., pair of compounds featuring the same biphenyl-like system and a different terminal fragment) are always lower than 2-fold.

Discussion
The binding affinities of the target compounds are listed in Table 2. The analysis of the binding data at the screened receptors provides information regarding the role of the three structural elements that compose the structure, i.e., biphenyl-like system, spacer, and terminal fragment.
Considering the affinities for the 5-HT 1A receptor, we can note that both biphenyl-like systems and terminal fragments selected for the study are compatible with high affinity at the receptor. In fact, four compounds out of 10 show K i values lower than 17.6 nM (compounds 8a,b, 11b, 12b). The differences in affinity are related to the length of the spacer between the piperazine and the terminal fragment. Examining the compounds characterized by the oxotriazinyl terminal fragment (compounds 9a,b-11a,b), the derivatives with a tetramethylene spacer have higher affinity than the trimethylene counterparts. The influence of the spacer length on affinity is related to the nature of the terminal fragment. In fact, considering the compounds featuring the bipyridyl system linked to the piperazine, the K i ratios of the trimethylene and tetramethylene homologs are variable (11a/11b = 52; 12a/12b = 3.7; 8a/8b = 0.7).
As for the 5-HT 2A receptor, the compounds featuring the bipyridyl system linked to the piperazine display moderate to low affinity. Within this group of compounds, the terminal fragment plays a role on affinity. In fact, the compounds with the oxotriazinyl fragment have a low affinity, whereas a bicyclic terminal fragment imparts comparatively higher affinity. The replacement of the bipyridyl with a biphenyl system has a favorable effect on the affinity. In fact, for the compounds with a tetramethylene spacer, 9b has a slightly higher affinity than 10b which shows much higher affinity than 11b. A similar trend can be observed for the trimethylene spacer derivatives 9a, 10a, and 11a.
The affinity values of the target compounds at 5-HT 7 receptor are distributed in a narrower range as compared to the other receptors under study. In fact, the highest affinity values were shown by 8b (K i = 9.38 nM) and the lowest by 11b (K i = 79.4 nM). Consequently, the differences in affinity between homologs (trimethylene or tetramethylene spacer derivatives) or analogs (i.e., pair of compounds featuring the same biphenyl-like system and a different terminal fragment) are always lower than 2-fold.
As for the D 2 -like receptors (D 2 and D 3 ), it is immediately evident that the oxotriazinyl terminal fragment leads to compounds devoid of affinity (compounds 9a,b-11a,b). By contrast, the bicyclic terminal fragments lead from moderate to high affinity (compounds 8a,b, 12a,b) to both D 2 and D 3 receptors, being the affinity for the former receptor slightly higher. For these compounds, the tetramethylene derivatives have slightly higher affinity than the trimethylene analogs.
Regarding the D 1 -like receptors (D 1 and D 5 ), the biphenyl-like system influences the affinity to a greater extent than the terminal fragment. In fact, the bipyridyl fragment leads to low affinity for D 1 and D 5 receptors, whereas the biphenyl fragment to moderate affinity. Considering compounds 9a,b and 10a,b, the trimethylene or tetramethylene homologs display a difference in affinity lower than 2-fold.
Collectively, the binding data of the target compounds revealed that the three structural elements that compose the structure have different effects on the affinity for the receptors under study. For 5-HT 1A receptor, variations of length of the spacer result in the highest variations in affinity. For 5-HT 2A receptors, the nature of the biphenyl-like system has high influence on the affinity. For 5-HT 7 receptor, the biphenyl-like systems linked to the piperazine ring confirmed their favorable role on affinity, with no substantial effects of both the spacer length and the terminal fragments. For D 2 -like receptor, the terminal fragment has a pronounced impact on affinity, whereas for D 1 -like receptor biphenyl-like system has a greater effect.
To identify molecular contacts responsible for the observed affinity profile, compounds 9b and 12b were docked at serotonin 5-HT 1A and 5-HT 2A and dopamine D 2 and D 1 receptors. We found that compounds orientations were analogous for all the studied receptors, except for the dopamine D 1 receptor, with the piperazine moiety forming a strong charged-assisted hydrogen bond with the highly conserved aspartate residue on the third transmembrane helix and the aryl substituent linked to the piperazine ring deeply buried in the binding site. By comparing the docking poses at serotonin 5-HT 2A receptor of compounds 9b (5-HT 2A K i = 39.4 nM) and 12b (5-HT 2A K i = 1492 nM), it is evident that the dimension of the terminal fragment impacts the interaction of the molecule with the receptor. In fact, the oxotriazinyl terminal fragment of compound 9b formed an additional contact with Leu229 of ECL2. In contrast, the bicyclic terminal fragment of compound 12b is too big to allow a proper accommodation in the binding pocket and leads to a steric clash. Similarly, compounds 9b and 12b showed very different binding poses at the dopamine D 2 receptor, which agree with the very different K i values (compound 9b K i = >10,000 nM, compound 12b K i = 17.0 nM). In fact, while the arylpiperazine moiety of compound 12b is located deeper into the binding site and the terminal fragment is orientated towards the extracellular milieu forming an additional contact with Asn396, compound 9b showed an opposite orientation with the arylpiperazine pointing towards the extracellular milieu that is not in agreement with the observed structure-activity relationships.
Binding affinity data and docking poses indicated that, when the structural requirements for high affinity for a given set of receptors overlap, compounds with different profiles can be identified. In fact, compounds 9b and 10b have a balanced profile concerning the affinity for 5-HT 1A , 5-HT 2A , and 5-HT 7 receptors, which is the desired affinity profile to obtain prospective drugs for the treatment of the core symptoms of ASD. In addition, the selectivity over D 2 receptors of both 9b and 10b is indeed favorable in terms of potential extrapyramidal side effects. On the other hand, compounds 8b and 12a have an affinity profile compatible with antipsychotic activity, as they have an affinity for the D 2 receptor, being 12a preferable for the balanced affinity for D 2 and 5-HT 2A receptors, accompanied by a comparable affinity for 5-HT 1A and 5-HT 7 receptors.
Next, we studied the potential of the target compounds to protect against apoptosis induced by H 2 O 2 using SH-SY5Y neuroblastoma cells as an experimental model.
We first evaluated the cytotoxicity of the target compounds alone and found that compounds 11a and 11b were not cytotoxic in the dose range of 0.1 to 100 µM, while the remaining compounds showed moderate cytotoxicity, having EC 50 values between 20-50 µM ( Table 2). Of note, SH-SY5Y cells treatment with each compound at 1 µM or 5 µM concentration for 24 h had a marginal effect on cell viability, with a decrease lower than 4% or 8% at 1 µM or 5 µM, respectively.
Next, the protective effect of each compound was determined by pretreating the cells with the test compound for 3 h at the non-toxic doses of 1 µM and 5 µM, and subsequent treatment with 400 µM H 2 O 2 for an additional 24 h. Treatment of SH-SY5Y cells with 400 µM H 2 O 2 for 24 h caused a reduction of cell viability of 72-76% (Figure 3). Pretreatment of SH -SY5Y cells with compounds 8a,b, 9b, and 10a-12a resulted in a statistically significant recovery of cell viability with respect to treatment with the 400 µM H 2 O 2 alone at both concentrations studied (Figure 3). Thus, compounds 8b, 9b, and 12a stand out of the set because they combine a favorable affinity profile with antioxidant properties.
The target compounds were screened for their in vitro metabolic stability to evaluate their liability to metabolic degradation by first-pass oxidative metabolism-the main cause of metabolic degradation in vivo [17]. The metabolic stability was assessed as the percentage of the parent compound recovered after 30 min of incubation with rat microsomes in the presence of an NADPH-regenerating system (Table 1). Compound 9b was the most stable compound as it showed a percentage of recovery (45%) that is predictive of low clearance in vivo, based on our previous data on a broad set of analog arylpiperazine derivatives [5,18]. Among the remaining compounds, 10b and 12a showed borderline stability value, while the others were massively metabolized. By comparing the cLogP values and metabolic stability of each compound, it emerges that the compound's oxidative liability is not related to lipophilicity, as shown by pairs of compounds 9a,b and 10a,b (Table 1).

Chemicals
Chemicals were purchased from Sigma-Aldrich, Alfa Aesar, TCI Chemicals. Unless otherwise stated, all chemicals were used without further purification. Thin layer chromatography (TLC) was performed using plates from Merck (silica gel 60 F254). Column chromatography was performed with 1:30 Merck silica gel 60 Å (63-200 µm) as the stationary phase. Flash chromatographic separations were performed on a Biotage SP1 purification system using flash cartridges pre-packed with KP-Sil 32−63 µm, 60 Å silica. 1 H NMR spectra were recorded on a Varian Mercury-VX spectrometer (300 MHz) or on a 500-vnmrs500 Agilent spectrometer (500 MHz). All chemical shift values are reported in ppm (δ). Recording of mass spectra was done on an HP6890-5973 MSD gas chromatograph/mass spectrometer; only significant m/z peaks, with their percentage of relative intensity in parentheses, are reported. HRMS-ESI analyses were performed on a Bruker Datonics MicrOTOF-Q II mass spectrometer, mass range 50-800 m/z, electrospray ion source in positive or negative ion mode. All spectra were in accordance with the assigned structures. Elemental analyses (C,H,N) of the target compounds were performed on a Eurovector Euro EA 3000 analyzer. Analyses indicated by the symbols of the elements were within ± 0.4% of the theoretical values. The purity of the target compounds listed in Table 1 was assessed by RP-HPLC and combustion analysis. All compounds showed ≥95% purity. RP-HPLC analysis was performed on an Agilent 1260 Infinity Binary LC System equipped with a diode array detector using a Phenomenex Gemini C-18 column (250 mm × 4.6 mm, 5 µm particle size). All target compounds (Table 1) were eluted with CH 3 CN/ammonium formate 50 mM pH 5, 8:2 (v/v) at a flow rate of 1 mL/min.

5-HT 2A Receptor
The affinity of the compounds for serotonin 5-HT 2A receptor was evaluated in membrane preparations from CHO-K1 cells stably expressing the human cloned receptor, following previously described procedures [23] with minor modifications. Competition binding experiments were performed using 1 nM [ 3 H]-ketanserin (50.3 Ci/mmol, 1 mCi/mL, PerkinElmer NET791250UC) as radioligand. Nonspecific binding was assessed in the presence of 1 µM methysergide. Risperidone was included in the assays as reference compound. Assays were carried out in duplicate in 96-well assay plates, in assay buffer (50 mM Tris-HCl, pH = 7.4). Assay mixtures were incubated at 37 • C for 30 min, followed by filtration through GF/B glass filter plates, and wash with ice-cold wash buffer (50 mM Tris-HCl, pH = 6.6). Affinity (equilibrium dissociation constant, K i ) of the compounds was calculated from competition binding curves as indicated above. K i value of the reference ligand risperidone was 0.20 ± 0.04 nM.

5-HT 7 Receptor
The affinity of the compounds for serotonin 5-HT 7 receptor was evaluated in membrane preparations from HEK293 cells stably expressing the human cloned 5-HT 7A receptor, following previously described procedures [24] with minor modifications. Competition binding experiments were performed using 2 nM [ 3 H]-SB269970 (34.5 Ci/mmol, 0.25 mCi/mL, Perkin Elmer NET1198U250UC) as radioligand. Nonspecific binding was assessed in the presence of 25 µM clozapine. Methiothepine was included in the assays as reference compound. Assays were carried out in duplicate in 96-well assay plates, in assay buffer (50 mM Tris-HCl, 4 mM MgCl 2 , 1 mM ascorbic acid, 0.1 mM pargyline, pH = 7.4). Assay mixtures were incubated at 37 • C for 60 min, followed by filtration through GF/C glass filter plates, and washing with ice-cold assay buffer. Affinity (equilibrium dissociation constant, K i ) of the compounds was calculated from competition binding curves as indicated above. K i value of the reference ligand methiothepin was 1.33 ± 0.32 nM.

Dopamine Receptors
Competition binding experiments at dopamine D 2s and D 3 receptor were performed to determine the affinity of the compounds. Cell culture conditions are described in Frank et al. [25]. Membrane preparations of transfected CHO-K1 cells stably expressing the human D 2s or D 3 receptor were performed according to Bautista-Aguilera et al. [26]. The competition binding experiments were conducted as previously reported [26,27]. Briefly, membrane fractions (D 2s R: 25 µg/200 µL; D 3 R: 20 µg/200 µL) were incubated with [ 3 H]spiperone (0.2 nM) and test compound for 120 min at room temperature. Seven appropriate concentrations of compound between 100 µM and 0.01 nM were used. Nonspecific binding was determined in the presence of 10 µM haloperidol. Haloperidol was also used as reference compound. Separation of the bound ligand from the free ligand was conducted by filtration through GF/B filters using deionized water.
The affinity for dopamine D 1 and D 5 receptors was evaluated in membrane fractions of HEK-293 cells stably expressing the human dopamine D 1 or D 5 receptor. HEK-293-D 1 cells were cultured in DMEM/F12 (15 mM HEPES, 1.2 g/L NaHCO 3 ) with 10% (v/v) FBS and 1% (v/v) L-glutamine. HEK-293-D 5 cells were cultured in DMEM/F12 (with L-glutamine, 1.2 g/L NaHCO 3 ) with 20% (v/v) FBS. Preparation of membrane fractions and radioligand competition assays were performed as previously described [26]. Briefly, membrane fractions were co-incubated with [ 3 H]-SCH23390 (0.3 nM) and different concentrations of test compound for 120 min at room temperature. Fluphenazine (100 µM) was used to determine non-specific binding and as reference compound.
Data of at least three independent experiments in triplicate were analyzed using GraphPad Prism 8.1 (San Diego, CA, USA). Binding curves were fitted to a non-linear regression model (one-site competition). K i values were calculated from IC 50 values using the Cheng-Prusoff equation. The statistical calculations were performed on -log(K i ). Results are presented in mean K i (nM) with 95% confidence intervals (nM).

Docking Studies
The structures of compounds 8b, 9b, 12a, and 12b were loaded into MOE 2019.0102 (Chemical Computing Group, Montreal, Canada). Subsequently, the protonation state was adjusted to dominant at pH = 7.0 and the compounds were subjected to energy minimization. X-ray structures of 5-HT 1A (PDB code: 7E2Z [28]), 5-HT 2A (PDB code: 6WHA [29]), D 1 (PDB code: 7JOZ [30]), and D 2 (PDB code: 6CM4 [31]) were loaded into MOE 2019.0102 and subjected to "QuickPrep" routine, which includes correction of missing loops, adjustment of protonation state and energy minimization. For docking, the initial placement was performed by setting the central interaction towards the acidic aspartate as essential. Initial scoring of 30 poses was performed using London dG scoring function, while refinement was performed with a rigid receptor, and the five best poses were selected based on GBVI/WSA dG scoring function. The best pose was subjected to energy minimization and inspected manually.

Stability Assays in Rat Liver Microsomes
Test compounds were pre-incubated at 37 • C with rat liver microsomes (Tebu-Bio, Milan, Italy) (1.0 mg/mL microsomal protein) at 10 µM final concentration in 100 mM potassium phosphate buffer (pH 7.4) for 10 min. Metabolic reactions were initiated by the addition of the NADPH regenerating system (containing 10 mM NADP, 50 mM glucose-6phosphate, and 10 unit/mL glucose-6-phosphate dehydrogenase, final glucose-6-phosphate dehydrogenase concentration, 1 unit/mL). After 30 min incubation, the reaction mixture was quenched by adding an equal volume of cold acetonitrile containing the internal standard. A test compound incubated with microsomes without NADPH regenerating system was included. Quenched samples were centrifuged at 4500 rpm for 15 min and the supernatants were injected for quantification analysis. Samples (100 µL) were analyzed by using an Agilent 1260 Infinity Binary LC System equipped with a diode array detector (Open Lab software was used to analyze the chromatographic data) and a Phenomenex Gemini C-18 column (250 mm × 4.6 mm, 5 µm particle size). The samples were eluted using CH 3 CN/20 mM ammonium formate pH 5.5 (7:3, v/v) as eluent (volumetric flow rate =1 mL/min). Concentrations were quantified by measuring the area under the peak.
The percentage of the parent compound remaining after a 30-min incubation was calculated according to the equation: % of parent compound remaining after 30 min = C parent /C control × 100 (1) where C parent is ligand concentration after incubation with microsome fraction and NADPH regenerating system and C control is ligand concentration after incubation with microsome fraction only.

Cell Viability
Determination of cell growth was performed using the MTT assay at 48 h. [32] On day 1, 25,000 cells/well were seeded into 96-well plates in a volume of 100 µL. On day 2, the various drug concentrations (1 µM-100 µM) were added. In all the experiments, the various drug-solvents (EtOH, DMSO) were added in each control to evaluate possible solvent cytotoxicity. After the established incubation time with drugs (48 h), MTT (0.5 mg/mL) was added to each well, and after 3-4 h incubation at 37 • C, the supernatant was removed. The formazan crystals were solubilized using 100 µL of DMSO/EtOH (1:1, v/v) and the absorbance values at 570 and 630 nm were determined on the microplate reader Victor 3 (PerkinElmer Life Sciences, Waltham, MA, USA).

Evaluation of Cell Viability
Determination of neuroprotection against H 2 O 2 was performed using the MTT assay at 24 h. [33] On day 1, 25,000 cells/well were seeded into 96-well plates in a volume of 100 µL. On day 2, the cells were pretreated (3 h) with two drug concentrations (1 µM and 5 µM) before being incubated with 400 µM H 2 O 2 for 24 h. In all the experiments, the various drug-solvents (EtOH, DMSO) were added in each control to evaluate a possible solvent cytotoxicity. After 24 h incubation, MTT (0.5 mg/mL) was added to each well, and after 3-4 h incubation at 37 • C, the supernatant was removed. The formazan crystals were solubilized using 100 µL of DMSO/EtOH (1:1, v/v) and the absorbance values at 570 and 630 nm were determined on the microplate reader Victor 3 from PerkinElmer Life Sciences (Waltham, MA, USA).

Statistical Analysis
Data were analyzed by applying the one-way repeated measures analysis of variance, and Bonferroni's multiple comparison test followed as a post hoc test. Results are reported as mean SD of at least two to three independent experiments, performed in triplicate. Statistical significance was accepted at p < 0.05.

Conclusions
In the present study, we have described the design, synthesis, and biological evaluation of a set of long-chain arylpiperazine derivatives. The structural modification led to the identification of new compounds displaying an array of affinity for serotonin 5-HT 1A , 5-HT 2A , 5-HT 7 receptors, and dopamine D 2 receptor and, in some cases, antioxidant properties. Binding affinity data evidenced that: (i) the nature of the terminal fragment had an impact mostly on the affinity at dopamine D 1 -and D 2 -like receptors; (ii) the length of the linker influenced the affinity at serotonin 5-HT 1A receptor; (iii) the nature of the biphenyl-like system linked to the piperazine ring influenced the affinity at serotonin 5-HT 2A and 5-HT 7 receptors. The most interesting compounds were: 12a that combines an affinity profile compatible with antipsychotic activity (affinity for dopamine D 2 (K i = 300 nM) and 5-HT 2A (K i = 315 nM) receptors, accompanied by an affinity for 5-HT 1A (K i = 41.5 nM) and 5-HT 7 (K i = 42.5 nM) receptors) and antioxidant properties; 9b that has an affinity profile compatible for studies in the context of ASD (affinity for serotonin 5-HT 1A (K i = 23.9 nM), 5-HT 2A , (K i = 39.4 nM), 5-HT 7 (K i = 45.0 nM) receptors and selectivity over dopamine D 2 receptors) and antioxidant properties. Even though the new compounds showed CNS MultiParameter Optimization score predictive of desirable ADMET properties and cross the blood-brain barrier, they fail to achieve in vitro metabolic stability suitable for studies in vivo. The only exception is compound 9b that therefore deserves further characterization.