Discovery of Potential, Dual-Active Histamine H3 Receptor Ligands with Combined Antioxidant Properties

In an attempt to find new dual acting histamine H3 receptor (H3R) ligands, we designed a series of compounds, structurally based on previously described in our group, a highly active and selective human histamine H3 receptor (hH3R) ligand KSK63. As a result, 15 obtained compounds show moderate hH3R affinity, the best being the compound 17 (hH3R Ki = 518 nM). Docking to the histamine H3R homology model revealed two possible binding modes, with key interactions retained in both cases. In an attempt to find possible dual acting ligands, selected compounds were tested for antioxidant properties. Compound 16 (hH3R Ki = 592 nM) showed the strongest antioxidant properties at the concentration of 10−4 mol/L. It significantly reduced the amount of free radicals presenting 50–60% of ascorbic acid activity in the 2,2-diphenyl-1-picrylhydrazyl (DPPH) assay, as well as showed antioxidative properties in the ferric reducing antioxidant power (FRAP) assay. Despite the yet unknown antioxidation mechanism and moderate hH3R affinity, 16 (QD13) constitutes a starting point for the search of potential dual acting H3R ligands-promising tools for the treatment of neurological disorders associated with increased neuronal oxidative stress.


Introduction
Histamine H 3 receptor (H 3 R) is an outstanding member of the histamine receptor family. Since its discovery and gene cloning, it has served as a widely explored drug target due to its broad spectrum of neuromodulatory effects in the central nervous system (CNS) [1]. On the one hand, activation of presynaptic receptors on histaminergic neurons results in regulation of synthesis and release of the natural ligand histamine. On the other hand, activation of H 3 R as a heteroreceptor influences the whole range of other neurotransmitters, including acetylcholine, dopamine, serotonin, norepinephrine, γ-aminobutyric acid, glutamate, and substance P, to name a few [2][3][4]. Anatomically speaking, H 3 Rs are prominently expressed in humans in the basal ganglia, globus pallidus, hippocampus, and cortex-regions associated with cognition, sleep, wakefulness, and homeostatic regulation [5]. Consequently, H 3 R ligands, either antagonist or inverse agonists, as they also demonstrate constitutive activity, contribute to one of the largest groups of contemporary GPCR ligand research, with implications to treat a number of CNS-based diseases, including Alzheimer's (AD) and Parkinson's disease (PD), schizophrenia, ADHD, obesity, and narcolepsy [6,7]. However, while several big pharma companies were working on development of H 3 R inhibiting ligand and reached the Phase III studies in these indications, In light of the presented findings, we designed the series of ligands with modifications in the hydrophilic and distal regulatory regions, namely by introducing different substituents in the piperazine (R 1 ) and benzophenone (R 2 ) moieties, respectively, based on the general structure depicted in Figure 2. The high activity associated with the basic part of the molecule, prompts to look for new scaffolds, even if many attempts lead to  In light of the presented findings, we designed the series of ligands with modifications in the hydrophilic and distal regulatory regions, namely by introducing different substituents in the piperazine (R 1 ) and benzophenone (R 2 ) moieties, respectively, based on the general structure depicted in Figure 2. The high activity associated with the basic part of the molecule, prompts to look for new scaffolds, even if many attempts lead to inactive compounds. For this series, we have chosen five different (cyclo)alkyl carbonyl and 4-arylpiperazine substituents. Using direct (cyclo)alkyl carbonyl piperazine N4 nitrogen substituents, we wanted to test whether such reduction in piperazine N4 nitrogen negative charge would still be acceptable for the retention of H 3 R affinity [27]. On the other hand, in order to avoid the possible overlap of inductive and mesomeric effects of p-substituted methoxy group, we have chosen m-substituted derivatives, so as to take only the former one into account. On the eastern end, we tested whether the increase of the hydrophobicity, by introduction of halogen atoms to benzophenone moiety, is also tolerated.

Chemistry
The desired final compounds 4-18 were obtained through the synthetic route presented in Scheme 1, according to the previously reported methods [23,30]. Starting phenoxy-alkyl bromides were obtained in one-step O-alkylation of commercially available (substituted) p-hydroxy benzophenones with 1,3-dibromopropane, refluxed either in sodium propanolate (1) or acetone with powdered potassium carbonate, and a catalytic amount of potassium iodide (2,3). Obtained precursor bromides were then coupled with commercially available 4-substituted piperazinyl derivatives in the mixture of ethanol/water with powdered potassium carbonate and a catalytic amount of potassium iodide. Final products were obtained either as quick crystallizing free bases (5, 10-15), or were isolated as oxalic acid salts (4, 6-9, [16][17][18]. For detailed information, please refer to the Materials and Methods section. Hence, it is relevant to search for novel bioactive compounds able to combine in one molecule multiple action properties. The search for possible dual active H 3 R ligands, which might act as neuroprotectives, both directly (by expressing anti-oxidative properties) and indirectly (through the effect of H 3 R antagonism), seems necessary, as their applications might be beneficial in future therapies. Even though the obtained compounds lack the phenolic protons believed to be essential for the antioxidant activity, paying attention to, e.g., the structure of the known antioxidant melatonin [28], with respect to recent findings on the tetratargeting Contilisant [26,29], selected compounds were tested for possible antioxidative properties. Last, but not least, biological test results were complemented by the docking studies to the histamine H 3 receptor homology model.

Chemistry
The desired final compounds 4-18 were obtained through the synthetic route presented in Scheme 1, according to the previously reported methods [23,30]. Starting phenoxyalkyl bromides were obtained in one-step O-alkylation of commercially available (substituted) p-hydroxy benzophenones with 1,3-dibromopropane, refluxed either in sodium propanolate (1) or acetone with powdered potassium carbonate, and a catalytic amount of potassium iodide (2,3). Obtained precursor bromides were then coupled with commercially available 4-substituted piperazinyl derivatives in the mixture of ethanol/water with powdered potassium carbonate and a catalytic amount of potassium iodide. Final products were obtained either as quick crystallizing free bases (5, 10-15), or were isolated as oxalic acid salts (4, 6-9, 16-18). For detailed information, please refer to the Materials and Methods section. amount of potassium iodide (2,3). Obtained precursor bromides were then coupled with commercially available 4-substituted piperazinyl derivatives in the mixture of ethanol/water with powdered potassium carbonate and a catalytic amount of potassium iodide. Final products were obtained either as quick crystallizing free bases (5, 10-15), or were isolated as oxalic acid salts (4, 6-9, 16-18). For detailed information, please refer to the Materials and Methods section. Scheme 1. General synthetic pathway for compounds 4-18.

Biological Evaluation
All compounds, in either basic or oxalate forms, were then tested in H3R binding studies in vitro, as described previously [31,32]. In brief, compounds were tested at five to eleven appropriate concentrations in a [ 3 H]N α -methylhistamine (KD = 3.08 nM) radioligand depletion assay to determine the affinity for the human recombinant histamine H3R Scheme 1. General synthetic pathway for compounds 4-18.

Biological Evaluation
All compounds, in either basic or oxalate forms, were then tested in H 3 R binding studies in vitro, as described previously [31,32]. In brief, compounds were tested at five to eleven appropriate concentrations in a [ 3 H]N α -methylhistamine (K D = 3.08 nM) radioligand depletion assay to determine the affinity for the human recombinant histamine H 3 R stably expressed in HEK-293 cells. For poorly soluble compounds, 16 and 17, an addition to 0.025% HCl in the binding test buffer was used.
Human H 3 R binding characterization of the tested compounds are presented in Table 1. Most of the tested compounds show relatively low to no histamine H 3 receptor affinity, with the exception of p-nitro-(14, 15) and m-methoxyphenyl derivatives (16)(17)(18), which display affinities at moderate levels, with compound 17 being the most affine (hH 3 R K i = 518 nM) of the group. Regarding the piperazine substituents, one could observe that direct substitution of piperazine N4 nitrogen with simple (cyclo)alkyl carbonyl moieties has negative influence on hH 3 R affinity. Slightly higher values were demonstrated for p-acetylphenyl derivatives; however, still at an unacceptable level. It appears that such a group is not suitable for histamine H 3 receptor binding. The reasoning behind these observations might be the low pK a value of N4 nitrogen of the obtained compounds. With, e.g., pK a = −6.0 for compounds 4-6, it might be assumed that the majority of the compounds at physiological pH are in non-ionized form. This might result in a lack of crucial interactions with binding pocket amino acids, resulting in very low or no affinity at all. It appears that only m-methoxy- (16)(17)(18) and p-nitrophenyl derivatives (14, 15) express the pK a at a level that allows for abundance of ionized forms. Calculated pK a values can be found in Table S1 of the Supplementary Material. Antioxidant activity at the level above 50% of ascorbic acid activity at a concentration of 10 −4 mol/L might be considered noteworthy. In the DPPH assay, compound 16 (QD13) at a concentration of 10 −4 mol/L caused a decrease in absorbance by 12.37 ± 1.31% when compared to the sample containing the solvent + reaction mixture (blank, maximum concentration of the DPPH radical). Moreover, 16, at a concentration of 10 −4 , had 60.28% of the ascorbic acid activity (10 −4 mol/L, calculated from Table 2). In fact, of all the tested compounds, 16 expressed the strongest antioxidant properties. Antioxidant activity at the level above 50% of ascorbic acid activity at a concentration of 10 −4 mol/L might be considered noteworthy. In the DPPH assay, compound 16 (QD13) at a concentration of 10 −4 mol/L caused a decrease in absorbance by 12.37 ± 1.31% when compared to the sample containing the solvent + reaction mixture (blank, maximum concentration of the DPPH radical). Moreover, 16, at a concentration of 10 −4 , had 60.28% of the ascorbic acid activity (10 −4 mol/L, calculated from Table 2). In fact, of all the tested compounds, 16 expressed the strongest antioxidant properties. Antioxidant activity at the level above 50% of ascorbic acid activity at a concentration of 10 −4 mol/L might be considered noteworthy. In the DPPH assay, compound 16 (QD13) at a concentration of 10 −4 mol/L caused a decrease in absorbance by 12.37 ± 1.31% when compared to the sample containing the solvent + reaction mixture (blank, maximum concentration of the DPPH radical). Moreover, 16, at a concentration of 10 −4 , had 60.28% of the ascorbic acid activity (10 −4 mol/L, calculated from Table 2). In fact, of all the tested compounds, 16 expressed the strongest antioxidant properties. of 10 −4 mol/L might be considered noteworthy. In the DPPH assay, compound 16 (QD13) at a concentration of 10 −4 mol/L caused a decrease in absorbance by 12.37 ± 1.31% when compared to the sample containing the solvent + reaction mixture (blank, maximum concentration of the DPPH radical). Moreover, 16, at a concentration of 10 −4 , had 60.28% of the ascorbic acid activity (10 −4 mol/L, calculated from Table 2). In fact, of all the tested compounds, 16 expressed the strongest antioxidant properties. of 10 −4 mol/L might be considered noteworthy. In the DPPH assay, compound 16 (QD13) at a concentration of 10 −4 mol/L caused a decrease in absorbance by 12.37 ± 1.31% when compared to the sample containing the solvent + reaction mixture (blank, maximum concentration of the DPPH radical). Moreover, 16, at a concentration of 10 −4 , had 60.28% of the ascorbic acid activity (10 −4 mol/L, calculated from Table 2). In fact, of all the tested compounds, 16 expressed the strongest antioxidant properties. Eight of the other compounds showed a slightly smaller but pronounced decrease in absorbance (antioxidant activity >20% of ascorbic acid activity at 10 −4 mol/L); thus, scavenging DPPH free radicals ( Table 2, left panel).
On the other hand, in the FRAP assay, compounds 16-18 showed highest activity of all the tested compounds. The level of absorbance at a concentration 10 −4 mol/L was above 8% of ascorbic acid activity at a concentration of 10 −4 mol/L ( Eight of the other compounds showed a slightly smaller but pronounced decrease in absorbance (antioxidant activity >20% of ascorbic acid activity at 10 −4 mol/L); thus, scavenging DPPH free radicals ( Table 2, left panel).
On the other hand, in the FRAP assay, compounds 16-18 showed highest activity of all the tested compounds. The level of absorbance at a concentration 10 −4 mol/L was above 8% of ascorbic acid activity at a concentration of 10 −4 mol/L ( Eight of the other compounds showed a slightly smaller but pronounced decrease in absorbance (antioxidant activity >20% of ascorbic acid activity at 10 −4 mol/L); thus, scavenging DPPH free radicals ( Table 2, left panel).
On the other hand, in the FRAP assay, compounds 16-18 showed highest activity of all the tested compounds. The level of absorbance at a concentration 10 −4 mol/L was above 8% of ascorbic acid activity at a concentration of 10 −4 mol/L ( Eight of the other compounds showed a slightly smaller but pronounced decrease in absorbance (antioxidant activity >20% of ascorbic acid activity at 10 −4 mol/L); thus, scavenging DPPH free radicals ( Table 2, left panel).
On the other hand, in the FRAP assay, compounds 16-18 showed highest activity of all the tested compounds. The level of absorbance at a concentration 10 −4 mol/L was above 8% of ascorbic acid activity at a concentration of 10 −4 mol/L ( Table 2, right panel).
Considering both tests results, compound 16 showed the strongest antioxidant properties. At a concentration of 10 −4 mol/L, it significantly reduced the number of free radicals (DPPH test) and exhibited antioxidative properties in general (FRAP test) despite yet unknown mechanism. Furthermore, other compounds from the m-methoxyphenyl subgroup, 17 and 18, also showed antioxidant activity in FRAP assay. Among three various Cl >1000

QD3
Molecules 2021, 25 Eight of the other compounds showed a slightly smaller but pronounced decrease in absorbance (antioxidant activity >20% of ascorbic acid activity at 10 −4 mol/L); thus, scavenging DPPH free radicals ( Table 2, left panel).
On the other hand, in the FRAP assay, compounds 16-18 showed highest activity of all the tested compounds. The level of absorbance at a concentration 10 −4 mol/L was above 8% of ascorbic acid activity at a concentration of 10 −4 mol/L ( Eight of the other compounds showed a slightly smaller but pronounced decrease in absorbance (antioxidant activity >20% of ascorbic acid activity at 10 −4 mol/L); thus, scavenging DPPH free radicals ( Table 2, left panel).
On the other hand, in the FRAP assay, compounds 16-18 showed highest activity of all the tested compounds. The level of absorbance at a concentration 10 −4 mol/L was above 8% of ascorbic acid activity at a concentration of 10 −4 mol/L ( Eight of the other compounds showed a slightly smaller but pronounced decrease in absorbance (antioxidant activity >20% of ascorbic acid activity at 10 −4 mol/L); thus, scavenging DPPH free radicals ( Table 2, left panel).
On the other hand, in the FRAP assay, compounds 16-18 showed highest activity of all the tested compounds. The level of absorbance at a concentration 10 −4 mol/L was above 8% of ascorbic acid activity at a concentration of 10 −4 mol/L ( Eight of the other compounds showed a slightly smaller but pronounced decrease in absorbance (antioxidant activity >20% of ascorbic acid activity at 10 −4 mol/L); thus, scavenging DPPH free radicals ( Table 2, left panel).
On the other hand, in the FRAP assay, compounds 16-18 showed highest activity of all the tested compounds. The level of absorbance at a concentration 10 −4 mol/L was above 8% of ascorbic acid activity at a concentration of 10 −4 mol/L ( Eight of the other compounds showed a slightly smaller but pronounced decrease in absorbance (antioxidant activity >20% of ascorbic acid activity at 10 −4 mol/L); thus, scavenging DPPH free radicals ( Table 2, left panel).
On the other hand, in the FRAP assay, compounds 16-18 showed highest activity of all the tested compounds. The level of absorbance at a concentration 10 −4 mol/L was above 8% of ascorbic acid activity at a concentration of 10 −4 mol/L ( Eight of the other compounds showed a slightly smaller but pronounced decrease in absorbance (antioxidant activity >20% of ascorbic acid activity at 10 −4 mol/L); thus, scavenging DPPH free radicals ( Table 2, left panel).
On the other hand, in the FRAP assay, compounds 16-18 showed highest activity of all the tested compounds. The level of absorbance at a concentration 10 −4 mol/L was above 8% of ascorbic acid activity at a concentration of 10 −4 mol/L ( Eight of the other compounds showed a slightly smaller but pronounced decrease in absorbance (antioxidant activity >20% of ascorbic acid activity at 10 −4 mol/L); thus, scavenging DPPH free radicals ( Table 2, left panel).
On the other hand, in the FRAP assay, compounds 16-18 showed highest activity of all the tested compounds. The level of absorbance at a concentration 10 −4 mol/L was above 8% of ascorbic acid activity at a concentration of 10 −4 mol/L ( Table 2,   The overall affinity row for 4-piperazinyl substituents can be summarized as mmethoxyphenyl > p-nitrophenyl > p-acetylphenyl > cyclopropylmethanone > acyl (best to worst) for the herein presented series. Moreover, no clear structure-affinity relationship among the three benzophenone modified derivatives can be drawn.

Antioxidative Properties
In order to determine the possible antioxidant properties, selected compounds were tested in 2,2-diphenyl-1-picrylhydrazyl (DPPH) and ferric reducing antioxidant power (FRAP) assay with ascorbic acid as a reference compound (Table 2). For details, please refer to the Materials and Methods Section. Graphical representation of tests results can be found in Supplementary Material. Antioxidant activity at the level above 50% of ascorbic acid activity at a concentration of 10 −4 mol/L might be considered noteworthy. In the DPPH assay, compound 16 (QD13) at a concentration of 10 −4 mol/L caused a decrease in absorbance by 12.37 ± 1.31% when compared to the sample containing the solvent + reaction mixture (blank, maximum concentration of the DPPH radical). Moreover, 16, at a concentration of 10 −4 , had 60.28% of the ascorbic acid activity (10 −4 mol/L, calculated from Table 2). In fact, of all the tested compounds, 16 expressed the strongest antioxidant properties.
Eight of the other compounds showed a slightly smaller but pronounced decrease in absorbance (antioxidant activity > 20% of ascorbic acid activity at 10 −4 mol/L); thus, scavenging DPPH free radicals ( Table 2, left panel).
On the other hand, in the FRAP assay, compounds 16-18 showed highest activity of all the tested compounds. The level of absorbance at a concentration 10 −4 mol/L was above 8% of ascorbic acid activity at a concentration of 10 −4 mol/L ( Table 2, right panel).
Considering both tests results, compound 16 showed the strongest antioxidant properties. At a concentration of 10 −4 mol/L, it significantly reduced the number of free radicals (DPPH test) and exhibited antioxidative properties in general (FRAP test) despite yet unknown mechanism. Furthermore, other compounds from the m-methoxyphenyl subgroup, 17 and 18, also showed antioxidant activity in FRAP assay. Among three various phenyl (Ph) substituents, only methoxy derivatives appeared to possess antioxidant properties. This might be due to the fact that the p-NO 2 -Ph, p-COCH 3 -Ph are capable of attracting electrons, instead of giving them, while p-OCH 3 -Ph are proficient in giving electrons rather than capturing them, which in turn is an indication of their antioxidant activity. However, with such an electron withdrawing substituent as m-OCH 3 , antioxidant activity appears in the studied group of compounds.
Overall, these results prompt for future investigation of compounds 16-18 antioxidant properties in various disease entities. Despite its moderate affinity for histamine H 3 receptors, 16 (QD13) could be investigated as an adjuvant to known therapies in many diseases, to which contributes an oxidative stress and where modulation of the brain levels of histamine and other neurotransmitters provide improvement.

Docking Studies
Three out of four known and described histamine receptors still remain a Gordian knot by means of crystallography. To date, only histamine H 1 receptor crystal structure was resolved (PDB ID: 3RZE [33]). Therefore, we used a previously described histamine H 3 R homology model in our study, which was constructed on the template of crystal structure of M 2 muscarinic acetylcholine receptor (PDB ID: 3UON) [22]. Considering that, due to very low basicity, the abundance of compounds at physiological pH might appear in non-ionized forms, for docking studies, we prepared both protonated and non-protonated conformers.
All of the compounds were docked and, independently of the form, fit the H 3 R binding pocket. In case of non-protonated conformers, compounds were placed shallower in the binding pocket, and mostly stabilized through π-π stacking or additional H-bonds with Y374 (3.51) or R381 (6.58). However, no interactions with either of the key amino acids, E206 (5.46) and D114 (3.32), were found. For possible non-ionized forms, this might explain very low or no affinity for the histamine H 3 receptor. Calculated binding poses for non-protonated conformers can be found in the Supplementary Material (Table S2).
Nonetheless, two different binding modes were surprisingly observed for protonated conformers: "standard (ST)", with the eastern region pointing to the extracellular space, and "upside-down (UD)" (Figure 3). The positions of protonated basic moieties of the ligands were mostly similar, independently of the mode. In both cases, a crucial histamine H 3 R antagonist/inverse agonist interaction-salt bridge and/or hydrogen bond formation between protonated amine nitrogen and E206 (5.46) was, however, preserved due to the approximate location of the protonated nitrogen (superscripts denote Ballesteros-Weinstein numbering [34]). and W371 (6.48), in the case most active for this series of m-methoxyphenyl derivatives, was also present. These interactions might explain the overall higher H3R affinity of the ligands from this subgroup. On the other hand, the UD mode was observed for p-acetylphenyl and p-nitrophenyl derivatives, with the exception of 18, for which both modes were observed. Although, the ST mode expressed slightly lower glide energy (−58.78 vs −56.08 kcal/mol). In this mode, the benzophenone fragments reached deeper towards the narrow sub pocket localized among the transmembrane regions (TM) 1, TM3, and TM7. The latter is formed by lateral L117 (3.35) and W371 (6.48), and S121 (3.39) and A122 (3.40) from the bottom. Additional pose stabilization through cation-π interactions between protonated piperazine nitrogen and F193 (ECL2), as well as the hydrogen bond between the carbonyl/nitro group and R381 (6.58) was found. Along with the novel compounds, reference structure KSK63 was also docked. In this case the pyridyl moiety nitrogen is protonated instead of the piperazine one. This allows for a shift of the ligand towards TM3 and formation of ionic interaction with key aminer- In the ST mode, benzophenone fragments occupied the space fenced by aromatic features of F193, Y189 (ECL2), and Y394 (7.35) on the sides, and Y91 (2.61) and Y94 (2.64) on the top. In most cases, a hydrogen bond between the carbonyl group and R381(6.58) was found. The latter one was not observed for inactive compounds 4-9. Additional stabilization through halogen bond formation with Y91 (2.61) and π-π stacking with Y394 (7.35) and W371 (6.48), in the case most active for this series of m-methoxyphenyl derivatives, was also present. These interactions might explain the overall higher H 3 R affinity of the ligands from this subgroup.
On the other hand, the UD mode was observed for p-acetylphenyl and p-nitrophenyl derivatives, with the exception of 18, for which both modes were observed. Although, the ST mode expressed slightly lower glide energy (−58.78 vs −56.08 kcal/mol). In this mode, the benzophenone fragments reached deeper towards the narrow sub pocket localized among the transmembrane regions (TM) 1, TM3, and TM7. The latter is formed by lateral L117 (3.35) and W371 (6.48), and S121 (3.39) and A122 (3.40) from the bottom. Additional pose stabilization through cation-π interactions between protonated piperazine nitrogen and F193 (ECL2), as well as the hydrogen bond between the carbonyl/nitro group and R381 (6.58) was found.
Along with the novel compounds, reference structure KSK63 was also docked. In this case the pyridyl moiety nitrogen is protonated instead of the piperazine one. This allows for a shift of the ligand towards TM3 and formation of ionic interaction with key aminergic GPCR anchor, D114 (3.32) (Figure 3). Lack of such interactions observed for the novel, described compounds, might also be responsible for their much lower affinity towards histamine H 3 receptor.
The stability of the calculated poses for this groups' most affine compounds 15-17 and both orientations of 18 along with the reference compound KSK63 was further evaluated by means of molecular dynamics (MD) simulations. From each simulation, 7 poses were selected (starting pose, and after each 100ps up to 600ps). In the case of 15-17 ( Figure 4) and KSK63, complexes appeared stable through the whole 600ps simulation, retaining the key interactions, and the potential energy (U) of the atomic system at the level of 1000 kcal/mol. In all cases, starting poses (marked in grey, Figure 4) were very similar to its orientation at the end of simulation (marked yellow). However, a shift in the binding site with retained conformation appeared in the first 100ps of the simulation. On the other hand, the analysis of KSK63 behavior allowed for possible explanation of its' high affinity. The conformation remained quite stable during the simulation with consistent set of interactions. What was not observed for the remaining structures, the protonated pyridyl nitrogen retained the interactions with not only the key ionic anchor D114 (3.32) but also with E206 (5.46) , through the whole simulation. Detailed analysis of compounds behavior during MD simulations is provided by examining the changes in their interactions with H 3 R ( Figure 5). The results indicate a relatively consistent set of ligand-protein interactions occurring during the whole MD simulation. Most of the interactions occur within the TM3, TM5, TM6 and TM7 and ECL2, with consistent interactions with Y115 (3.33), Y189 (45.51), E206 (5.46), Y394 (7.35). In the case of 18, the standard orientation seems more stable during the simulation, with most of the starting key interactions retained in the last frame. were selected (starting pose, and after each 100ps up to 600ps). In the case of 15-17 ( Figure  4) and KSK63, complexes appeared stable through the whole 600ps simulation, retaining the key interactions, and the potential energy (U) of the atomic system at the level of ~1000 kcal/mol. In all cases, starting poses (marked in grey, Figure 4) were very similar to its orientation at the end of simulation (marked yellow). However, a shift in the binding site with retained conformation appeared in the first 100ps of the simulation. On the other hand, the analysis of KSK63 behavior allowed for possible explanation of its' high affinity. The conformation remained quite stable during the simulation with consistent set of interactions. What was not observed for the remaining structures, the protonated pyridyl nitrogen retained the interactions with not only the key ionic anchor D114 (3.32) but also with E206 (5.46), through the whole simulation. Detailed analysis of compounds behavior during MD simulations is provided by examining the changes in their interactions with H3R ( Figure 5). The results indicate a relatively consistent set of ligand-protein interactions occurring during the whole MD simulation. Most of the interactions occur within the TM3, TM5, TM6 and TM7 and ECL2, with consistent interactions with Y115 (3.33), Y189 (45.51), E206 (5.46), Y394 (7.35). In the case of 18, the standard orientation seems more stable during the simulation, with most of the starting key interactions retained in the last frame.

Chemistry
All reagents were purchased from commercial suppliers (Alfa Aesar, Sigma-Aldrich, CHESS GmbH, Chempur) and were used without further purification. Melting points

Chemistry
All reagents were purchased from commercial suppliers (Alfa Aesar, Sigma-Aldrich, CHESS GmbH, Chempur) and were used without further purification. Melting points (mp) were determined on a Büchi Melting Point M560 apparatus (BÜCHI Labortechnik AG, Flawil, Switzerland) and were uncorrected. Liquid chromatography/mass spectra (LC/MS) were obtained using Waters TQ Detector mass spectrometer (Milford, MA, USA). 1 H NMR spectra were recorded on a Varian Mercury 300 MHz PFG (Palo Alto, CA, USA) spectrometer in DMSO-d 6 . Chemical shifts were expressed in parts per million (ppm) using the solvent signal as an internal standard. Data have been reported in the following order: multiplicity (s, singlet; d, doublet; t, triplet; qi, quintet; m, multiplet; br, broad; Ac, acetyl; Prop, propionyl; t-b, tert-butyl; t-p, tert-pentyl; Ph, phenyl; Pip, piperazine), approximate coupling constants J expressed in Hertz (Hz), number of protons. 13 C NMR spectra were recorded on Varian-Mercury-VX 300 MHz PFG or Bruker 400 MHz spectrometer (Berlin, Germany) at 75 MHz in DMSO-d 6 . LC-MS separation was carried out on a system consisting of a Waters ACQUITY UPLC, coupled to a Waters TQD mass spectrometer. Retention times (t R ) were stated in minutes. For flash chromatography (FC) purification, silica gel 60 (0.063-0.20 mm; Merck) and the following solvent systems were used: Ipetroleum ether:EtOAc (97.5:1.25); II-CH 2 Cl 2 :MeOH (95:5).

Procedure 1
To a solution of proper hydroxy benzophenone (0.01 mol) in acetone (60 mL), anhydrous potassium carbonate (1.12 g, 0.008 mol) and catalytic amount of potassium iodide were added. The reaction mixture was stirred for 10 min at room temperature, after which 1,3-dibromopropane (12 g, 0.06 mol) was added. The reaction mixture was then heated at reflux for 24 h, after which the inorganic solids were filtered, and acetone was evaporated. Remaining oils were then purified by flash chromatography (system I). Procedure was performed according to the literature reference [30]. To a solution of freshly prepared sodium propanolate (25 mL propanol; 0.575 g, 0.025 mol Na), (4fluorophenyl) (4-hydroxyphenyl)methanone (5.4 g, 0.025 mol) was added and stirred at room temperature (RT) for 15 min; 1,3-dibromopropane (15.14 g, 0.075 mol) was then added dropwise over one hour. The reaction mixture was stirred at 60 • C for 3 h, and then refluxed for another 3 h. After cooling down to RT reaction mixture was filtrated and evaporated. To a resulting brown oil, MeOH was added, and purified by flash chromatography (system II); obtaining 6 g of yellowish, quickly crystallizing oil. Final products were obtained according to the procedure described in [22] and purified by flash chromatography (CH 2 Cl 2 /MeOH, 95:5). Resulting oils were either: transformed into oxalic acid salts, using 10% excess of oxalic acid solution in absolute ethanol in RT, and then precipitated by addition of ethyl ether or, in the case of self-crystalizing oils, remained free bases.
Scintillation data (c.p.m.) corrected for non-specific binding were analyzed using GraphPad Prism (V6.01, San Diego, CA, USA) software, using non-linear least squares/regression fit. K i values were calculated from IC 50 values according to the Cheng-Prusoff equation [35]. Statistical analyses were performed on pK i values from at least three experiments, each performed, at least, in duplicates. Mean affinity values (K i ) were transferred into nanomolar concentrations (with 95% confidence interval).

In Vitro Antioxidant Activity
Antioxidant properties of compounds were tested in vitro in two different assays: the 2,2-diphenyl-1-picrylhydrazyl (DPPH) assay [36,37] and the FeCl 3 reduction activity assay (FRAP, ferric reducing antioxidant power) [37]. The tested compounds and ascorbic acid were dissolved in DMSO at the 10 −2 mol/L concentration, and in ethanol at the concentration of 10 −4 mol/L respectively. Moreover, 30% and 50% ascorbic acid activity cut-off for DPPH and FRAP assays respectively were set ad hoc.

DPPH Assay
DPPH is a molecule containing a stable free radical. The reduction of DPPH (purple) in ethanol solution takes place in the presence of a hydrogen-donating antioxidant due to the formation of the non-radical form of DPPH-H (yellow). This transformation results in color change from purple to yellow, which can be measured spectrophotometrically. Thus, the decrease in absorbance is proportional to the decrease in concentration of DPPH (free radical) and discoloration of reaction mixture indicates the scavenging of free radicals by a tested antioxidant.
In order to determine antioxidant capacity, 20 µL of the compound solutions (dissolved in 96% ethanol) was mixed with 180 µL of 0.3 mM ethanolic DPPH solution. The change in the absorbance was detected at 517 nm after 30 min of incubation. Results were expressed as percentage decrease in absorbance of the tested sample, compared to the sample containing the solvent + reaction mixture (blank with maximum concentration of the DPPH radical). L-ascorbic acid was used as a reference compound. Calculations were made using MS Excel, from the proportion where the change of absorbance in the vitamin C (vit. C) group was treated as 100%, according to Equation (1): Calculation of % Vit. C activity: ∆A b -the change in the absorbance of the test sample and ∆A Vit. C is change in the absorbance of reference sample, or with vit. C.

FRAP Assay
A modified method of Benzie and Strain [38] was adopted for the FRAP assay. The stock solutions included 300 mM acetate buffer (3.1 g of C 2 H 3 NaO 2 × 3H 2 O and 16 mL of C 2 H 4 O 2 ), pH 3.6, 10 mM TPTZ (2,4,6-tripyridyl-s-triazine) solution in 40 mM HCl, and 20 mM FeCl 3 × 6H 2 O solution. The working solution was prepared by mixing 10 parts of acetate buffer, 1 part of TPTZ solution, and 1 part of FeCl 3 × 6H 2 O solution. A portion of 180 µL of the FRAP solution was mixed with 20 µL of the tested compound solution and incubated at room temperature for 10 min in the dark. Readings of the color product (ferrous tripyridyl triazine complex) were then taken at 593 nm against ethanol. Results are expressed as an increase in absorbance of the test sample compared to a sample containing the solvent. L-ascorbic acid was used as reference compound.

Molecular Modeling
For docking purposes, Schrödinger Maestro Suite (v. 11.5.011, Release 2019-01) was used [39]. Ligands were built in either their forms (protonated N1 piperazine nitrogen, structure charge +1) or non-ionized forms. Bioactive conformations were generated using ConfGen module [40,41] (force field: OPLS3, water environment; minimization method: Polak-Ribier Conjugate Gradient (PRCG) with maximum iterations of 2500 and a convergence threshold of 0.05; conformational search of 100 steps per rotational bond). For all compounds, the five lowest energy conformers were selected for docking studies. Possible binding pocket adaptation in the presence of certain ligands was examined using an induced fit refinement protocol [42] with 17. In order to validate the methods used, the ligand was then re-docked with high confidence. The site was then centered on the ligand. Docking to rigid the form of the receptor was performed using the Glide module [43,44] (extra precision, flexible ligand sampling, and a maximum of five poses per conformer). Ligands were rated according to their position in the binding pocket, interactions with the binding pocket amino acids, as well as the docking score value. Ligand-receptor visualizations were generated using Schrödinger Maestro.

Conclusions
Looking for a new, active structures, is not an easy task. It even becomes harder when the reference structure is of one of the highest affinity hH 3 R ligands obtained to date; thus, increasing the possibility of ending up with lower affinity compounds. However, obtaining such compounds teaches us a lesson on structural uniqueness.
In an attempt to find new possibly dual acting H 3 R ligands, in this study we obtained a series of compounds of moderate hH 3 R affinity, ranging around K i of 500 nM (17), on the one hand, and on the other hand, expressing antioxidative properties (e.g., 16, hH 3 R K i = 592 nM, showing 50-60% of vitamin C activity at a 10 −4 mol/L). To our knowledge, alongside a recently published Tetratarget ligand Contilisant and its predecessor ASS234 [27], obtained herein derivatives are some of the few histamine H 3 R ligands with described antioxidative properties. From a structural point of view, it appeared that simple (cyclo)alkyl substituents at the 4-position of piperazine negatively impact the ligand's activity at the desired target, when compared to KSK63 (Figure 1). Nonetheless, in the case of simple acyl moiety, a benzene ring separating the acyl from the basic core, through additional interactions, might force slightly higher affinity within our novel series (e.g., 10 vs. 4), yet still on the low level. Interestingly, docking into the histamine H 3 R homology model revealed two putative binding modes, with ionic key interactions retained in both cases. The stability of such obtained complexes was demonstrated by MD simulations. This also allowed for possible explanation of KSK63 high affinity. Last, but not least, in an attempt to find possible dual acting ligands, selected compounds were tested for antioxidant properties. Among others, compound 16 (hH 3 R K i = 592 nM) showed the strongest antioxidant properties at a 10 −4 mol/L concentration. Compound 16 not only significantly reduced the amount of free radicals showing 50-60% of AA activity, but had also exhibited antioxidative properties in general. Despite yet unknown mechanisms and moderate hH 3 R affinity, compound 16 (QD13) constitutes a novel, potential dual acting H 3 R ligand-promising starting point towards treatment of neurologic disorders associated with increased neuronal oxidative stress.
In conclusion, although the obtained series of compounds appeared to be weak histamine H 3 receptors ligands, they might be valuable pharmacological tools in the search for novel molecules, targeting oxidative stress-based diseases. Identification of strong histamine H 3 receptor ligands, with antioxidative properties, will drive our further, planned studies.
Supplementary Materials: The following are available online, The Supplementary Material contain Table S1: Calculated pKa values for piperazine N1 and N4 nitrogen, Table S2: Putative binding poses for non-protonated conformers, Table S3: Putative binding poses for protonated conformers, Table  S4: Molecular dynamics frames alignment for selected ligands and Figure S1: Antioxidant activity of the tested compounds in DPPH and FRAP assay.

Institutional Review Board Statement: Not applicable.
Informed Consent Statement: Not applicable.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.