Nature-Inspired O-Benzyl Oxime-Based Derivatives as New Dual-Acting Agents Targeting Aldose Reductase and Oxidative Stress

Aldose reductase (ALR2) is the enzyme in charge of developing cellular toxicity caused by diabetic hyperglycemia, which in turn leads to the generation of reactive oxygen species triggering oxidative stress. Therefore, inhibiting ALR2 while pursuing a concomitant anti-oxidant activity through dual-acting agents is now recognized as the gold standard treatment for preventing or at least delaying the progression of diabetic complications. Herein we describe a novel series of (E)-benzaldehyde O-benzyl oximes 6a–e, 7a–e, 8a–e, and 9–11 as ALR2 inhibitors endowed with anti-oxidant properties. Inspired by the natural products, the synthesized derivatives are characterized by a different polyhydroxy substitution pattern on their benzaldehyde fragment, which proved crucial for both the enzyme inhibitory activity and the anti-oxidant capacity. Derivatives (E)-2,3,4-trihydroxybenzaldehyde O-(3-methoxybenzyl) oxime (7b) and (E)-2,3,4-trihydroxybenzaldehyde O-(4-methoxybenzyl) oxime (8b) turned out to be the most effective dual-acting products, proving to combine the best ALR2 inhibitory properties with significant anti-oxidant efficacy.


Introduction
Aldose reductase (alditol:NADP + oxidoreductase, EC 1.1.1.21, AKR1B1, ALR2) is commonly acknowledged as the enzyme in charge of the development of long-term diabetic complications. Therefore, its inhibition represents the key tool to prevent or at least delay the progression and severity of the pathological changes that diabetics face over time [1].
This consideration has prompted medicinal chemists to design many structurally different compounds as potent ALR2 inhibitors (ARIs), with a goal to develop a gold standard treatment for diabetic diseases affecting the heart, nerve, eye, and kidney. However, a true and widespread ARI therapy has not yet become a reality despite expectations. Epalrestat (1, Figure 1) is the only drug successfully marketed in Japan, India, and China to treat a true and widespread ARI therapy has not yet become a reality despite expectations. Epalrestat (1, Figure 1) is the only drug successfully marketed in Japan, India, and China to treat diabetic patients with neuropathy. Most of the compounds under investigation displayed a limited clinical efficacy, often due to unfavorable pharmacokinetic properties or even unexpected toxicities. Accordingly, the obtainment of effective and safe ARIs is still an unmet research goal deserving to be continued. A thorough survey of the most recent literature in the ARI field highlights an increasingly interest in natural compounds, mainly represented by polyhydroxylated flavonoids as quercetin (2, Figure 1), chalcones as butein (3, Figure 1), and coumarins as scopoletin (4, Figure 1), obtained from different plant species [2]. Besides fully addressing the current policy of the World Health Organization (WHO), which promotes the integration of traditional herbal medicines into the health care practices of countries [3], the use of these derivatives is justified by their multifunctional profile, which combines a promising ARI activity with a potent anti-oxidant efficacy. This latter aspect, in particular, is emerging as highly desirable in the setting up of novel ARIs. In fact, diabetes is a pathology characterized by both an increased formation of free radicals and a decreased antioxidant capacity, and the resulting oxidative stress plays a crucial role in the onset and progression of cell damage. Moreover, it may also trigger post-translational modification of key proteins, including ALR2, whose oxidation of the CYS299 residue leads to a marked reduction in sensitivity to ARIs, thus compromising the therapeutic effectiveness of these compounds [4,5].
Despite their biological versatility, polyhydroxylated natural derivatives generally struggle to become therapeutic agents, as most of them show low selectivity for the desired target and poor pharmacokinetic and pharmacodynamic properties. Therefore, instead of being claimed as active ingredients per se, these compounds should be better utilized as inspiring scaffolds by medicinal chemists, who have the tools to improve their functional profile through suitable structural modifications [6].
Pursuing our interest in the ARI field [7][8][9][10][11][12] and prompted by our successful experience of turning the flavonoid core into the bioisosteric pyrido [1,2-a]pyrimidin-4-one nucleus, by obtaining potent and selective anti-oxidant ARIs such as 5 ( Figure 1) [13], we developed a novel series of polyhydroxylated (E)-benzaldehyde O-benzyl oximes 6-11 ( Figure 1) as prototypes of dual-acting ARIs inspired by the stilbene derivative resveratrol A thorough survey of the most recent literature in the ARI field highlights an increasingly interest in natural compounds, mainly represented by polyhydroxylated flavonoids as quercetin (2, Figure 1), chalcones as butein (3, Figure 1), and coumarins as scopoletin (4, Figure 1), obtained from different plant species [2]. Besides fully addressing the current policy of the World Health Organization (WHO), which promotes the integration of traditional herbal medicines into the health care practices of countries [3], the use of these derivatives is justified by their multifunctional profile, which combines a promising ARI activity with a potent anti-oxidant efficacy. This latter aspect, in particular, is emerging as highly desirable in the setting up of novel ARIs. In fact, diabetes is a pathology characterized by both an increased formation of free radicals and a decreased antioxidant capacity, and the resulting oxidative stress plays a crucial role in the onset and progression of cell damage. Moreover, it may also trigger post-translational modification of key proteins, including ALR2, whose oxidation of the CYS299 residue leads to a marked reduction in sensitivity to ARIs, thus compromising the therapeutic effectiveness of these compounds [4,5].
Despite their biological versatility, polyhydroxylated natural derivatives generally struggle to become therapeutic agents, as most of them show low selectivity for the desired target and poor pharmacokinetic and pharmacodynamic properties. Therefore, instead of being claimed as active ingredients per se, these compounds should be better utilized as inspiring scaffolds by medicinal chemists, who have the tools to improve their functional profile through suitable structural modifications [6].
Pursuing our interest in the ARI field [7][8][9][10][11][12] and prompted by our successful experience of turning the flavonoid core into the bioisosteric pyrido [1,2-a]pyrimidin-4-one nucleus, by obtaining potent and selective anti-oxidant ARIs such as 5 ( Figure 1) [13], we developed a novel series of polyhydroxylated (E)-benzaldehyde O-benzyl oximes 6-11 ( Figure 1) as prototypes of dual-acting ARIs inspired by the stilbene derivative resveratrol (12, Figure 1), a natural compound overlooked by the ARI community. Different polyhydroxylated substitution patterns were explored on the A phenyl ring of the compounds ( Figure 1), being these groups decisive for the anti-oxidant activity and, possibly, for the interaction with the catalytic binding site of the enzyme, as demonstrated for selected natural derivatives [13]. In addition, methoxy residues and chlorine atoms were inserted on the B phenyl ring (Figure 1), as they are reported to increase the lipophilicity of ARIs improving their penetration through membranes and tissues, including the blood-brain barrier [14].
The novel compounds were investigated in vitro for their capacity to inhibit the human recombinant ALR2 and their anti-oxidant properties. To this end, both the TBARS and the DPPH assays were exploited to account for the ability of the compounds to protect from lipid peroxidation, which is significantly increased in diabetes even in wellcontrolled patients, and to scavenge free radicals, generated in a significant amount by the hyperglycemic conditions [15,16]. Moreover, docking simulations and molecular dynamic studies of the most active derivatives into the human ALR2 binding site were also carried out to rationalize the structure-activity relationships (SARs) observed.

Chemistry
Analytical grade reagents and solvents were purchased from Sigma-Aldrich (St. Louis, MO, USA), and were used as supplied. 1 H-NMR spectra were determined with a Bruker Ultrashield TM 400 MHz (Fällander, Switzerland). The following abbreviations are being used in this paper: singlet (s), doublet (d), triplet (t), doublet of doublets (dd), doublet of triplet (dt), and multiplet (m). Melting points (m.p.) were measured with a Leica Galen III microscope (Leica/Cambridge Instruments) and were uncorrected. Reactions were monitored by thin-layer chromatography (TLC) on silica gel plates containing a fluorescent indicator (Merck Silica Gel 60 F254); spots were detected under UV light (254 nm). Na 2 SO 4 was always used as the drying agent. Evaporation was carried out under reduced pressure (rotating evaporator). IR spectra of the synthesized compounds were recorded using a Cary 600 Series FTIR Spectrometer, (Agilent Technologies, Mulgrave, Victoria, Australia). The purity of the target inhibitors was determined by HPLC analysis, using a Shimadzu LC-20AD liquid chromatograph (PDA, 250-500 nm) and a Luna C18 column (250 mm × 4.6 mm, 5 µm, Phenomenex), with a gradient of 100% methanol and a flow rate of 0.2 mL/min. All the compounds showed percent purity values ≥95%.
General Procedure for the Synthesis of (E)-Polyhydroxy Substituted Benzaldehyde O-(benzyl) Oximes 6a-e, 7a-e, 8a-e, and 9-11 To a solution of the commercially available polyhydroxybenzaldehyde 17-21 (0.264 mmol) in MeOH (4 mL), an aqueous solution (1 mL) of the opportune hydroxylamine hydrochloride 16a-e (0.264 mmol) [17][18][19] was added. The reaction mixture was stirred for 1 h at room temperature until the disappearance of the starting material (TLC analysis). The mixture was evaporated to dryness and the resulting crude was added with water and extracted with EtOAc. The organic phase was washed (3x) with water, dried (Na 2 SO 4 ), filtered, and evaporated under reduced pressure to obtain a crude solid.
Expression of aldose reductase.
The expression strain BL21 (DE3) pLysS was obtained after a heat-shock transformation of competent cells with the expression plasmid. Transformed cells were used to inoculate pre-cultures containing 15 mL of LB media supplemented with ampicillin overnight at 37 • C. The following morning, a Fernbach flask containing 1l of ZYP-5052 automedium (25 mM Na 2 HPO 4 , 25 mM KH 2 PO 4 , 50 mM NH 4 Cl, 5 mM Na 2 SO 4 , pH 6.7, 5% glycerol, 0.5 g glucose, 2 g α-lactose, 10 g N-Z-amine AS, 5 g yeast extract, and 1 mM MgSO 4 ) and 100 µg/mL ampicillin was inoculated with the overnight pre-culture of BL21 (DE3) pLysS at 0.05 O.D. at 600 nm. ZYP-5052 medium is a buffered complex medium containing glucose, lactose, and glycerol formulated to induce protein expression after glucose depletion [21]. Expression was performed using a standardized two-step process. In the first part of fermentation, cells were grown at 37 • C to quickly reach the glucose depletion phase just before the induction. After that step at 0.9 O.D. at 600 nm (4 h), the temperature was lowered to 20 • C for 18 h to favor soluble protein expression. After 18 h, cells were pelleted by centrifugation and re-suspended in 80 mL of lysis buffer containing 100 mM Tris HCl pH 8, 150 mM NaCl, and 5% glycerol. 2 µL of benzonase and 0.1 M MgCl 2 were added to the suspension. After 30 min at 4 • C with gentle stirring, the cells were lysed by two passes at 1.5 kbar on EMULSIFLEX. After centrifugation for 45 min at 18.000 rpm at 4 • C, the soluble extract contained in the supernatant was separated from the insoluble fraction contained in the pellet.
Purification and TEV cleavage of DSBC-aldose reductase.
The soluble extract was filtered through 0.22 µM and loaded on a 5 mL HisTrap™FF column (GE Healthcare) with a 1 mL/min flow rate. Buffers used for purification were a binding buffer containing 100 mM tris-HCl pH 8, 150 mM NaCl, 5 mM Imidazole, 5% glycerol, and elution buffer containing 100 mM Tris-HCl pH 8, 150 mM NaCl, 500 mM imidazole, and 5% glycerol. The 6His-tagged fusion protein was eluted with a linear gradient (0 to 100% B in 30 min at a flow rate of 1 mL/min. The fractions containing the 6His-tagged fusion protein were pooled and dialyzed for 3 h against a buffer containing 100 mM Tris HCl pH 8, 150 mM NaCl using a Spectra/Por ® Dialysis Membrane (MWCO: 3500). The protein of interest was then cleaved with 10% (w/w) TEV protease overnight at 4 • C.
Purification of aldose reductase.
The TEV cleavage was loaded on a 5 mL HisTrap TM FF column (GE Healthcare) with a 1 mL/min flow rate. The protein DSBC and TEV protease with 6His-tag were bound on the column. The ALR2 was recovered in the unrestrained fraction and dialyzed for 3 h against a buffer containing 100 mM Tris HCl pH 8, 150 mM NaCl using a Spectra/Por ® Dialysis Membrane (MWCO: 3500). The last step was a purification of ALR2 on size-exclusion chromatography (SEC) Sephacryl ® S-100 HR GE Healthcare in buffer 100 mM Tris HCl pH 8, 150 mM NaCl (Figures S1 and S2, Supplementary Materials). ALR2 was concentrated to 1.47 mg/mL (41 µM).
The human recombinant ALR2 enzyme activity was determined spectrophotometrically by monitoring the change in absorbance at 340 nm, due to the oxidation of the NADPH cofactor catalyzed by the enzyme. According to a previously reported protocol, the enzymatic assay was performed in a 96-multiwell plate (250 µL) [7]. The reaction mixture, containing sodium phosphate buffer pH = 6.2 (0.1 M), NADPH (0.15 mM), ALR2 (0.050 µM), and D,L-glyceraldehyde (10 mM), was monitored for 5 min at 30 • C using a Victor3 Perkin Elmer Multilabel Plate Reader. The variation of pyridine coenzyme concentration versus time was analyzed through the WorkOut program. All the new synthesized compounds were solubilized in DMSO (10 mM) and screened at a final concentration of 100 µM. Assays were made in triplicate, and Tolrestat was used as the reference drug. Compounds 6b, 7b, and 8b, showing the best inhibitory activity at 100 µM, were further tested at additional concentrations within 50 µM and 10 µM.

Anti-Oxidant Activity Thiobarbituric Acid Reactive Substances (TBARS) Assay
A TBARS assay was carried out according to a previously reported protocol [13], using freshly isolated rat brain. Stock solutions of 10 mM test compounds were obtained in DMSO, then diluted to test samples with phosphate buffer pH 7.4. In a falcon tube, 100 µL of FeCl 3 (20 µM) and 100 µL of ascorbic acid (100 µM) were added to 100 µL of rat brain homogenate, in the absence or presence of 100 µL of test compounds. The resulting solution was diluted to 1 mL with phosphate buffer pH 7.4 and incubated at 37 • C for 30 min under slight stirring. Then, 0.5 mL of thiobarbituric acid (1% w/v in 0.05 N NaOH) and 0.5 mL of 25% v/v HCl were added, and the mixture was boiled for 10 min. After being cooled in an ice-cold water bath, the mixture was added with 3 mL of n-butanol and centrifuged at 2000× g for 10 min. The absorbance of the organic phase was recorded at 532 nm using a SPECTROstarNano (200-1000 nm) UV/Vis spectrophotometer. The synthesized compounds were tested in triplicate at 100 µM, 10 µM, and 1µM concentrations. Resveratrol and its dimethoxyderivative pterostilbene (Sigma-Aldrich) were used as the reference compound. TBARS was expressed as nmoles of malondialdehyde (MDA)/10 mg of rat brain tissue (wet weight), using a reference curve obtained with 1,1,3,3-tetramethoxypropane as the standard [22].

DPPH Assay
The DPPH assay was carried out in a 96-microwell plate (200 µL), following a previously described protocol [23]. Briefly, 195 µL of freshly prepared DPPH solution (60 µM in MeOH) was added with 5 µL of either test compounds or MeOH as the control. The plate was incubated in the dark for 1 h, at r.t. The absorbance was measured at 517 nm using a SPECTROstarNano (200-1000 nm) UV/Vis spectrophotometer. The synthesized compounds were tested in triplicate at 100 µM and 10 µM. Resveratrol (Sigma-Aldrich) was used as the reference standard. The percentage of the antioxidant activity (AA) of the compounds was determined according to the following Equation (1): where Abs DPPH is the absorbance of DPPH solution, Abs Blank1 is the absorbance of MeOH, Abs sample is the absorbance of DPPH solution containing the test compound, and Abs Blank2 is the absorbance of test compound solution without DPPH.

Molecular Docking Calculations
Docking studies were performed on the prospective inhibitors, including 6b, 7b, and 8b ( Table 1). The ligand compounds were drawn in MarvinSketch application-21.17 [28,29] and converted to mol2 format. Each compound was manually assessed in the Molegro Molecular View-7.0.0 (MMV) [30] to ascertain correct bond angles and hybridization state. Ligands were minimized using the steepest descent method and GAFF force field in Avogadro-1.2.0n [31,32]. The compounds were individually docked into the active site of the 2FZB, 4JIR, 4IGS, and 4JII aldose reductase protein. However, this study reported the grid box parameters for the systems involving 4JII because the conformation gave better binding energies with the inhibitors of interest. These were specified using AutodockTool-1.5.6 software [33,34] and defined the center as X = −6.304, Y = 6.647, and Z = 19.394 with dimensions X = 58, Y = 46, and Y = 66. The grid box was centered around TRP112 native conformation as it corresponds to the residual trio (TYR49, HIS111, and TRP112) of the anion pocket and specificity binding site of human ALR2 protein [24]. The structural and visual analyses of the 4JII ALR2-6b, ALR2-7b, and ALR2-8b were carried out using BIOVIA Discover studio Visualizer Software-4.0 [35,36]. Table 1. Substitution pattern of (E)-benzaldehyde O-benzyl oximes 6a-e, 7a-e, 8a-e, and 9-11. ter binding energies with the inhibitors of interest. These were specified using Autodock-Tool-1.5.6 software [33,34] and defined the center as X = −6.304, Y = 6.647, and Z = 19.394 with dimensions X = 58, Y = 46, and Y = 66. The grid box was centered around TRP112 native conformation as it corresponds to the residual trio (TYR49, HIS111, and TRP112) of the anion pocket and specificity binding site of human ALR2 protein [24]. The structural and visual analyses of the 4JII ALR2-6b, ALR2-7b, and ALR2-8b were carried out using BIOVIA Discover studio Visualizer Software-4.0 [35,36]. Table 1. Substitution pattern of (E)-benzaldehyde O-benzyl oximes 6a-e, 7a-e, 8a-e, and 9-11.

Retrieval and Preparation of Aldose Reductase
The X-ray crystallized ALR2 monomeric structure 4JII [25] was obtained from the RCSB Protein Data Bank. Compounds 6b (working code FB7), 7b (working code FB8), and 8b (working code FB6) were drawn in the MarvinSketch tool, and their geometries were optimized with the General Amber Force Field (GAFF) in Avogadro tool-1.2. On [31]. Both

. Molecular Dynamics Studies Retrieval and Preparation of Aldose Reductase
The X-ray crystallized ALR2 monomeric structure 4JII [25] was obtained from the RCSB Protein Data Bank. Compounds 6b (working code FB7), 7b (working code FB8), and 8b (working code FB6) were drawn in the MarvinSketch tool, and their geometries were optimized with the General Amber Force Field (GAFF) in Avogadro tool-1.2. On [31]. Both protein and inhibitors were prepared using the Chimeratool-1.14 [26] and the AutoDockTool-1.5.6 [37]. Moreover, the systems were opened and saved on the Molegro Molecular Viewer-2.5 [38] for the molecular dynamics simulation study. The preparation of the unbound and ligand-bound complexes involves removing nonstandard molecules, including water, NAP+, and sodium [27]. Figure S3 (Supplementary Materials) provides the structures of the ALR2 conformation showing some active site residues (A) and the chemical structures of the test inhibitors.

Molecular Dynamics (MD) Simulations
The MD simulation of unbound and bound ALR2 systems were carried out with Amber code using the graphic processor unit (GPU)-Compute Unified Device Architecture (CUDA) accelerated version of the AMBER18 package at the Lengau CHPC server [39] with integrated PMEND engine and LEAP modules [40]. The FF18SB force field variant was applied to describe the inhibitor-ALR2 enzyme complexes. The antechamber was used to add the atomic partial charges for the ligands using the restrained electrostatic potential (RESP) and the General Amber Force Field (GAFF) procedures [41]. The Leap module of Amber 18 enabled the addition of hydrogen atoms, sodium (Na + ), or chloride (Cl − ) counter ions to the systems for neutralization [42]. The complexes, implicitly, were suspended within an orthorhombic box of TIP3P water molecules to contain all the atoms within 8 Å of any box edges.
The systems were initially minimized for 2500 steps with 500 kcal/mol Å 2 restraint potential, and a whole minimization step of 5000 steps was further run without restraint using the conjugate algorithm. The heating step was carried out in a canonical ensemble condition (NVT) by a gradual heating simulation (0 to 300 K) executed for 5 ps while several atoms and volume were fixed. Moreover, the systems-containing solutes were imposed with a potential harmonic restraint of 10 kcal/mol Å 2 and a collision frequency of 1.0 ps −1 . An equilibration estimation of 1 ns was carried out while keeping the operating temperature (300 K) constant. Moreover, the number of atoms and pressure were constant, depicting an isobaric-isothermal ensemble (NPT). Using the Berendsen barostat, the complex pressure was kept at 1 bar [43,44]. The SHAKE algorithm was employed to construct the hydrogen atom bonds [43]. Coordinates and trajectories were printed and analyzed every 1 ps using the PTRAJ module available in AMBER18GPU [42,45]. The stabilities and flexibilities of the free protein and ligand-bound complexes were investigated statistically by estimating the averages RMSD [39] and RMSF [34], respectively. Furthermore, principal component analysis (PCA) [46] was calculated to unravel the extent of the enzyme's atomic displacement. In addition, the dynamic cross-correlation matrix data were computed to gain insights into the strength and similarities of the linear relationship between the interresidues of the bound and unbound systems. The detailed methods for obtaining RMSD, RMSF, PCA, and SASA [39,47] have been frequently reported in the literature, including our previous research works.

Binding Free Energy Computations
Over the years, molecular dynamics have been widely applied to study the mechanism of several biological processes by computing differences in the binding free energy [48]. The binding affinities of the bound and unbound systems were obtained by computing the binding free energies using the molecular mechanics/GB surface area method (MMGBSA) [49]. The detailed protocols of the methodology can be found in our previous works [39]. The free energies were calculated based on the average of over 40,000 snapshots extracted from the simulation trajectories. The estimated free binding energy, ∆G, for each molecular system, including the complex, inhibitor, and protein, can be given as: ∆G bind = E gas + G sol − TS E gas = E int + E vdw + E ele (4) ∆G bind is the summation of the gas phase, while E gas denotes the gas-phase energy. The solvation energy (G sol ) < the entropy (T∆S). E int term = the internal energy, E vdw term = the van der Waals energies, and E ele = the electrostatic interaction energies. The E gas was obtained from the AMBER18 FF14SB force fields module. The energy of the polar and the non-polar states was used to account for the solvation free energy. The product of the solvent accessible surface area (SASA) and a water probe radius of 1.4 Å was obtained to compute the non-polar solvation energy (G SA ). The S term is the total entropy of the system, which is often calculated using the normal mode method incorporated in the AMBER18 module [50].
tem, which is often calculated using the normal mode method incorporated in the A BER18 module [50].

Scheme 1. Synthesis of the (E)-benzaldehyde O-benzyl
The O-(arylmethyl)hydroxylamine hydrochlorides 16a-e were synthetized as pr ously described by a Mitsunobu reaction between the suitably substituted benzyl brom 13a-e and the N-hydroxyphthalimide 14, followed by the removal of the protec phthalimido group from the key intermediates 15a-e, carried out with ammonia solu 7N in MeOH [17,51,52]. Treatment of 16a-e with the opportune aldehydes 17-21, in m anolic solution and at room temperature, afforded the desired (E)-benzaldehyde O-ben oxime derivatives 6a-e, 7a-e, 8a-e, and 9-11 as the only geometric isomers. The cor configuration of the products was attributed by their 1 HNMR spectra on the basis of chemical shift value of the diagnostic iminic proton, ranging from 8.46 to 8.03 ppm reported by the literature for the E form [52]. The IR spectra of compounds 6a-e, 7a-e, Scheme 1. Synthesis of the (E)-benzaldehyde O-benzyl oximes 6a-e, 7a-e, 8a-e, 9-11.
The O-(arylmethyl)hydroxylamine hydrochlorides 16a-e were synthetized as previously described by a Mitsunobu reaction between the suitably substituted benzyl bromide 13a-e and the N-hydroxyphthalimide 14, followed by the removal of the protecting phthalimido group from the key intermediates 15a-e, carried out with ammonia solution 7N in MeOH [17,51,52]. Treatment of 16a-e with the opportune aldehydes 17-21, in methanolic solution and at room temperature, afforded the desired (E)-benzaldehyde O-benzyl oxime derivatives 6a-e, 7a-e, 8a-e, and 9-11 as the only geometric isomers. The correct configuration of the products was attributed by their 1 HNMR spectra on the basis of the chemical shift value of the diagnostic iminic proton, ranging from 8.46 to 8.03 ppm as reported by the literature for the E form [52]. The IR spectra of compounds 6a-e, 7a-e, 8a-e, and 9-11 were recorded as well. The obtained spectral signatures are listed in Table S1 (Supplementary Materials). All the synthesized derivatives were investigated in vitro for their dual activity. As for the inhibitory efficacy against the target ALR2, compounds were tested at 100 µM using a human recombinant enzyme. Results obtained, summarized in Table 2, highlight the prominent role of the polyhydroxy decoration of the A ring over the substituents on the B ring in determining the activity of the compounds. In fact, the 3,5-dihydroxy substitution pattern was detrimental, regardless of the residue on the B phenyl ring. See for example derivatives 6a, 7a, 8a, and 11, which were devoid of any activity. The 2,4-, 2,5-, and 3,4dihydroxy functionalization conferred some inhibitory efficacy, which was shown to be of the same order regardless of hydroxy groups' position on the A ring and the substituent of B ring. On the contrary, the 2,3,4-trihydroxy functionalization proved to be the most promising. Indeed, compound 6b, carrying a 2-methoxybenzyl moiety, almost halved ALR2 activity, and an even better efficacy was shown by 7b and 8b, bearing 3-methoxy and 4-methoxy residues on the B phenyl ring, respectively. They both stood out as the best performing compounds in the whole series, thus also emphasizing the importance of the substituent on the B ring to achieve the best interaction with the catalytic binding site of the enzyme. Regarding the anti-oxidant properties, compounds were tested in both TBARS and DPPH assays.
The TBARS assay highlights the ability of the compounds to block the hydroxy radicaldependent lipid peroxidation induced by the oxidant system Fe(III)/ascorbic acid in rat brain homogenate, and results in malondialdehyde (MDA) as the ultimate metabolite. Tested at 100 µM, all the synthesized compounds showed a relevant anti-oxidant efficacy ( Figure 2). Moreover, in this case, the reciprocal arrangement of the hydroxy substituents on the A phenyl ring of the derivatives proved to be a crucial determinant, remarkably affecting their functional potency. In particular, the 2,3,4-trihydroxy, 3,4-dihydroxy, and 2,5-dihydroxy combinations turned out to perform the best. Indeed, for 6b-d, 7b-d, and 8b-d, the obtained amounts of MDA (nmol per mg) in brain tissue were comparable or even lower than those attained with both resveratrol and its methoxylated derivatives, pterostilbene, used as the reference compounds ( Figure 2). Significantly, 6b-d, 7b,c, and 8b-d showed an excellent activity when tested at 10 µM (Figure 2), with 7b still active even at 1 µM ( Figure S4, Supplementary Materials). On the contrary, the anti-oxidant potency of the natural compounds came down on a dose-dependent basis. In light of the relevant results from the TBARS assay, compounds 6b-d, 7b-d, and 8b-d were also investigated for their free radical scavenging capacities by the DPPH assay. The selected derivatives showed remarkable efficacy (Figure 3). Notably, when tested at 100 μM concentration, they proved to be more potent than the reference resveratrol. As for derivatives 6c, 7c, and 8c, the same efficacy was maintained even at 10 μM test concentration, demonstrating that, among the explored polyhydroxy combinations, the 3,4-  6a-e, 7a-e, 8a-e, and 9-11, resveratrol, and pterostilbene on the production of thiobarbituric acid reactive substances (TBARS) in rat brain homogenate.
In light of the relevant results from the TBARS assay, compounds 6b-d, 7b-d, and 8b-d were also investigated for their free radical scavenging capacities by the DPPH assay. The selected derivatives showed remarkable efficacy ( Figure 3). Notably, when tested at 100 µM concentration, they proved to be more potent than the reference resveratrol. As for derivatives 6c, 7c, and 8c, the same efficacy was maintained even at 10 µM test concentration, demonstrating that, among the explored polyhydroxy combinations, the 3,4-dihydroxy catechol arrangement on the A ring ensures the achievement of the best free-radical quenching properties. O-benzyl oximes 6a-e, 7a-e, 8a-e, and 9-11, resveratrol, and pterostilbene on the production of thiobarbituric acid reactive substances (TBARS) in rat brain homogenate.

Figure 2. Effects of (E)-benzaldehyde
In light of the relevant results from the TBARS assay, compounds 6b-d, 7b-d, and 8b-d were also investigated for their free radical scavenging capacities by the DPPH assay. The selected derivatives showed remarkable efficacy ( Figure 3). Notably, when tested at 100 μM concentration, they proved to be more potent than the reference resveratrol. As for derivatives 6c, 7c, and 8c, the same efficacy was maintained even at 10 μM test concentration, demonstrating that, among the explored polyhydroxy combinations, the 3,4dihydroxy catechol arrangement on the A ring ensures the achievement of the best freeradical quenching properties. The promising ALR2 inhibitory activity of 6b, 7b, and 8b was rationalized through molecular docking studies, thereby representing an increasingly valuable tool in drug design [53]. The catalytic binding site of ALR2 is characterized by a high degree of plasticity [54][55][56]. From a functional point of view, such a feature provides a broad substrate The promising ALR2 inhibitory activity of 6b, 7b, and 8b was rationalized through molecular docking studies, thereby representing an increasingly valuable tool in drug design [53]. The catalytic binding site of ALR2 is characterized by a high degree of plasticity [54][55][56]. From a functional point of view, such a feature provides a broad substrate specificity to the enzyme. However, in structural terms, it implies the existence of different crystal arrangements obtained in the presence of different compounds, ALR2 being able to modify the conformation of its binding site by an induced-fit adaptability to the ligands. Therefore, to accomplish a reliable docking investigation, different conformations should be considered.
As reported in Table 3, docking calculations showed that 4JII conformation allowed to achieve the highest binding energy, indicating that ligands are more embedded in this active pocket than the other ALR2 reported structures. Moreover, the highest binding affinities were obtained for the ALR2-7b complex across all the investigated enzymeligand complexes. The resulting energies, including −7.61 kcal/mol, −8.1 kcal/mol, −8.4 kcal/mol, and −9.4 kcal/mol, demonstrate the highest affinity of the compound compared to the parents 6b and 8b, regardless of conformation used. Starting from the evidence, we focused on the 4JII ALR2 structure to investigate the binding forces and interaction types of the protein-inhibitor systems, to gain further insights into the binding and inhibitory mechanisms of the prospective inhibitors 6b, 7b, and 8b. Figure S5 (Supplementary Materials) highlights the ALR2 pocket residues interacting with the pharmacophoric portions of inhibitor 6b, including TRP21, VAL48, TRP112, VAL301, LEU302, and PHE321. The observed interactions were highly hydrophobic, indicating the extent to which the ligand was buried within the enzyme active site.
Similar results were also achieved by docking compound 8b into the same enzyme pocket ( Figure S6, Supplementary Materials). Besides the already mentioned VAL48, TRP112, VAL301, and LEU302, ASP300 and PHE123 were recruited to obtain the complex. Compound 8b was kept within the ALR2 binding site by highly hydrophobic interactions. Still, additional hydrogen bonds could also be formed between the 3-hydroxy group on the A ring of the compound and both VAL301 and LEU302, which provides the reason for the better inhibitory efficacy of 8b compared to 6b.
As for derivative 7b ( Figure S7, Supplementary Materials), the ALR2 residues contacting the inhibitor included TRP21, VAL48, TYR49, HIS111, TRP112, CYS299, VAL301, and LEU302. Similar to the parents 6b and 8b, the observed interactions were mainly hydrophobic. However, in addition to those previously described, classical hydrogen bond interactions involving VAL301 and a sulphur bond engaging CYS299 were also observed, which allowed 7b to be placed into the enzyme binding site. Moreover, differently from 6b and 8b, compound 7b can also interact with the NAD + co-factor, due to the 3-methoxy group on the B phenyl ring. On the whole, the occurring interactions maximize the binding of 7b within the target enzyme, rationalizing the observed inhibitory efficacy and justifying the highest binding energy resulting from the ALR2-7b complex.
The obtained ALR2-inhibitor complexes were compared with the unbound ALR2 system by evaluating the motion and dynamic of their structural behaviors using the all-atom molecular dynamics (MD) simulation for 200 ns. The per-residues energy decomposition (PRED) analysis identified the interacting residues and their contributing energies to the binding affinity [57]. The contributions of the ALR2 residues were unraveled on the Discovery Studio Viewer, which individually allowed insights into the interacting network of residue-inhibitor complexes ALR2-6b, ALR2-7b, and ALR2-8b (Figures 4A-6A). The PRED of ALR2 residues revealed some interactions with the ligands at the active site during the 200 ns MD simulations. Figures 4B-6B show that the inhibitors adopted favorable orientations, whereby the functional groups such as phenyl (-C 6 H 5 ), methoxy (-OCH 3 ), and hydroxyl (-OH) interacted with the residues by various covalent hydrogen bonds (C-H, pi-donor H-bond, H-H) and hydrophobic interactions, including pi-pi stacked, pi-pi T-shaped, alkyl, and pi-alkyl bonds. In particular, in the case of 8b, the methoxy group of the ligand formed strong alkyl and pi-alkyl interactions with the backbones of LYS22 and VAL48 and, to expatiate, the pi-donor H-bond with VAL48. The CYS299 oxygen interacted with the CH 2 -O group of the ligand, to form a C-H bond, while TRP80 and TRP112 formed pi-pi T-shaped and pi-pi stacked, with phenyl of trihydroxybenzaldehyde and methoxybenzyl groups of 8b, respectively. bonds (C-H, pi-donor H-bond, H-H) and hydrophobic interactions, including pi stacked, pi-pi T-shaped, alkyl, and pi-alkyl bonds. In particular, in the case of 8b, methoxy group of the ligand formed strong alkyl and pi-alkyl interactions with the ba bones of LYS22 and VAL48 and, to expatiate, the pi-donor H-bond with VAL48. T CYS299 oxygen interacted with the CH2-O group of the ligand, to form a C-H bond, wh TRP80 and TRP112 formed pi-pi T-shaped and pi-pi stacked, with phenyl of tri droxybenzaldehyde and methoxybenzyl groups of 8b, respectively.  The notable differences in the interactions across the three inhibitors are the type and the quality of hydrogen bonds formed within each complex. For instance, ALR2-8b and ALR2-7b showed more H-bonds and hydrophobic interactions than ALR2-6b, due to the higher residue numbers involved in the interaction. We also observed that 7b formed three strong non-polar covalent H-bonds with VAL301, LEU302, and SER304, which may have contributed to the high binding energy of the ALR2-7b complex ( Figure 4B). Although the inference may not be accurate for all protein-ligand interactions, the binding free energy calculated for each complex provided in Table 1 justified this assertion. The resonating phenyl rings of the trihydroxybenzaldehyde and methoxybenzyl moieties formed pi-alkyl bonds, pi-pi stacked, and T-shaped bonds. Other strong interactions were the different polar and non-polar hydrogen bonds formed, especially with compound 7b.  Furthermore, we obtained residue-based contributions from the decomposition of total binding energy (BPE) using the MM/GB(PB) method. The van der Waals (vdW) and electrostatic (Elec.) interactions were computed. The two forces showed that the residues and energy contributed more to the total binding energy. The stability of a protein-ligand system depends on intermolecular forces between the interacting residues and the ligands and, subsequently, the binding free energy (delta G) [58]. Figures 4C, 5C and 6C showed the individual hydrophobic residue pocketing 8b, 6b, and 7b, respectively, revealing their energy contributions. The 8b interacted with TRP21, LYS22, VAL48, TRP80, TRP112, and CYS299. Similarly, 6b and 7b interacted with TRP21, TRP112, and CYS299. Other residues that interacted with 6b included LYS22, TRP80, ASN300, VAL301, and LEU302, while 7b interacted with more residues, including PHE116, SER304, and PRO311.
The residues found in TRP21 and TRP112 and the energy contributions to the binding of the ligands were significant, as seen in Figures 4C, 5C and 6C. In contrasting these MD results with the molecular docking results, the findings showed that 8b-and 7bbound ALR2 complexes had higher average energies than the 6b-bound complexes; the contributed residues from TRP21 and TRP112 had higher binding energies to than 6b. The notable differences in the interactions across the three inhibitors are the type and the quality of hydrogen bonds formed within each complex. For instance, ALR2-8b and ALR2-7b showed more H-bonds and hydrophobic interactions than ALR2-6b, due to the higher residue numbers involved in the interaction. We also observed that 7b formed three strong non-polar covalent H-bonds with VAL301, LEU302, and SER304, which may have contributed to the high binding energy of the ALR2-7b complex ( Figure 4B). Although the inference may not be accurate for all protein-ligand interactions, the binding free energy calculated for each complex provided in Table 1 justified this assertion. The resonating phenyl rings of the trihydroxybenzaldehyde and methoxybenzyl moieties formed pi-alkyl bonds, pi-pi stacked, and T-shaped bonds. Other strong interactions were the different polar and non-polar hydrogen bonds formed, especially with compound 7b.
Furthermore, we obtained residue-based contributions from the decomposition of total binding energy (BPE) using the MM/GB(PB) method. The van der Waals (vdW) and electrostatic (Elec.) interactions were computed. The two forces showed that the residues and energy contributed more to the total binding energy. The stability of a protein-ligand system depends on intermolecular forces between the interacting residues and the ligands and, subsequently, the binding free energy (delta G) [58]. Figures 4C-6C showed the individual hydrophobic residue pocketing 8b, 6b, and 7b, respectively, revealing their energy contributions. The 8b interacted with TRP21, LYS22, VAL48, TRP80, TRP112, and CYS299. Similarly, 6b and 7b interacted with TRP21, TRP112, and CYS299. Other residues that interacted with 6b included LYS22, TRP80, ASN300, VAL301, and LEU302, while 7b interacted with more residues, including PHE116, SER304, and PRO311. The computational MM/GB(PB)SA approach [59] was used to obtain the total binding free energy (BFE) of the complexes (ALR2-6b, ALR2-7b, and ALR2-8b). The BFE analysis estimates distinct energy contributions within the protein binding pocket and orientation with the best intermolecular interactions. The approach provides insights into the binding affinities of the compounds. Table 4 showed the total binding free energy (kcal/mol) values calculated for the ALR2-6b (−25.40), ALR2-7b (−34.51), and ALR2-8b (−31.48) complexes. The binding energies for 6b, 7b, and 8b calculated from the molecular docking comparatively disagreed with the experimental data, unlike the binding energies calculated from the molecular dynamics studies. Simply put, the inhibition of the ALR2 enzyme by 8b (67% inhibition) and 7b (63% inhibition) are more potent and promising than 6b (44% inhibition). The molecular docking results revealed that binding energies of 8b, 6b, and 7b increase in ascending order, respectively. The MD studies showed that 8b and 7b have higher binding energies. Likewise, the PRED analysis and the ALR2-ligand interaction networks expatiated these findings, as seen in the immediate section. The molecular dynamics simulations method is more suitable to achieve a stable protein-ligand complex than molecular docking, while combining the two methods is more reliable [34,60]. Hence, we inferred that the binding affinity of the ligands, especially 8b, was not stable during molecular docking studies since the binding energy of the ALR2-8b was greater than 6b in the MD findings. Therefore, the results showed that the systems' electrostatic energies (∆Eelec) were lower than the van der Waal (∆E vdW ) forces across the three complexes.
The molecular dynamics (MD) simulation method was carried out to determine the atomic and molecular movements over 200 ns. Root mean square deviation [61] (RMSD) was used to estimate the differences in the C-α backbones of the Apo (initial) structural conformation and the bound ALR2 complexes. The observed deviations quantified the conformational stability of the protein during the MD simulations. Recall, more minor deviations indicate a more stable protein or protein-complex structure. The C-α backbone RMSD was determined for 200 ns simulation time to evaluate the stability of the unbound and bound ALR2 systems. Figure 7A shows that the four systems equilibrated immediately after simulation. Still, ALR2-Apo and ALR2-8b visibly fluctuated at about 60 ns with a small RMS deviation increment and remained stable again from 90 ns. In contrast, ALR2-6b and ALR2-7b appeared to maintain overall stability almost throughout the 200 ns. However, the average RMSD values for ALR2-Apo (1.41 Å), ALR2-8b (1.3 Å), ALR2-6b (1.25 Å), and ALR2-7b (1.21 Å) showed deviations in descending order, although this was minimal.
A C-α atom of proteins having the highest radius of gyration for the system size range suggests less tight packing proteins and vice versa. Figure 8A showed the RoG values of the inhibitor-bound ALR2 systems over 200 ns simulations decreased insignificantly at a pressure of 1 bar compared to the unbound protein. The average C-α atoms RoG of systems Apo-ALR2, ALR2-8b, ALR2-6b, and ALR2-7b were 19.36, 19.33, 19.33, 19.33, respectively. There was no difference in the RoG of complexes. Therefore, the Apo system seems less compact than the bound systems. Furthermore, the interactions and mechanisms of inhibition for the ligands were investigated by computing the principal component analysis (PCA). The analysis helps identify the most strongly correlated variables with each component or the numbers with the large magnitude farthest from zero in either direction [46]. The level at which a Furthermore, the interactions and mechanisms of inhibition for the ligands were investigated by computing the principal component analysis (PCA). The analysis helps identify the most strongly correlated variables with each component or the numbers with the large magnitude farthest from zero in either direction [46]. The level at which a correlation is essential is determined, for instance, a correlation above a given value may be deemed necessary. PCA reveals the displacement of the ALR2 C-α atoms during MD simulations performed using the CPPTRAJ in AMBER18 GPU by computing the first two components, PC1 and PC2. The MD simulations' convergence in the space was obtained from the averaged structure covariance matrix using Cartesian coordinates. PC1 and PC2 were plotted along the vertical and horizontal axes, respectively, which indicate covariance matrix after eigenvector elimination. Each point between the single-directional motions was uniquely orientated due to the overlapping of the typical structural conformations. Figure 8B illustrates the scatter plot generated by the unbound protein (black), ALR2-8b (red), ALR2-6b (green), and ALR2-7b (blue). The plot showed significant differences among the systems, indicating variations in the overall motion between the unbound system and the complexes. The critical protein motions along specific directions were represented by eigenvectors, given by the normal modes.
Further investigation into the inhibition actions of the 6b, 7b, and 8b against ALR2 revealed observation of significant visual structural changes in the loops proximal to and covering the pocket of the enzyme. The study of loops provides insights into various functions of proteins, including binding properties, dynamics, and shape [62]. Loop structures consist of 218-228 and 300-318 residues. We can refer to these residue regions as loops A and B. The trajectory snapshots were collected in nanoseconds at simulation time intervals (1, 5, 50, 100, 150, and 200). Figure 9 reveals the differential changes in the loop structures of the Apo, ALR2-6b, ALR2-7b, and ALR2-8b systems.
Biomolecules 2022, 12, x FOR PEER REVIEW 22 of 28 correlation is essential is determined, for instance, a correlation above a given value may be deemed necessary. PCA reveals the displacement of the ALR2 C-α atoms during MD simulations performed using the CPPTRAJ in AMBER18 GPU by computing the first two components, PC1 and PC2. The MD simulations' convergence in the space was obtained from the averaged structure covariance matrix using Cartesian coordinates. PC1 and PC2 were plotted along the vertical and horizontal axes, respectively, which indicate covariance matrix after eigenvector elimination. Each point between the single-directional motions was uniquely orientated due to the overlapping of the typical structural conformations. Figure 8B illustrates the scatter plot generated by the unbound protein (black), ALR2-8b (red), ALR2-6b (green), and ALR2-7b (blue). The plot showed significant differences among the systems, indicating variations in the overall motion between the unbound system and the complexes. The critical protein motions along specific directions were represented by eigenvectors, given by the normal modes. Further investigation into the inhibition actions of the 6b, 7b, and 8b against ALR2 revealed observation of significant visual structural changes in the loops proximal to and covering the pocket of the enzyme. The study of loops provides insights into various functions of proteins, including binding properties, dynamics, and shape [62]. Loop structures consist of 218-228 and 300-318 residues. We can refer to these residue regions as loops A and B. The trajectory snapshots were collected in nanoseconds at simulation time intervals (1, 5, 50, 100, 150, and 200). Figure 9 reveals the differential changes in the loop structures of the Apo, ALR2-6b, ALR2-7b, and ALR2-8b systems.  Interestingly, loop A snapshots showed that loop A moved to cover the ALR2 active site, especially in the ALR2-8b complex as the simulation progressed.
On the other hand, loop B displayed twisting behavior throughout the 200 ns simulations. Such structural changes are critical for the functional activity of proteins, including aldose reductase. For instance, the opening and closing of loop A could impact the ligandbinding terrain of the ALR2 protein to enhance the binding affinity of the inhibitors.
Therefore, we further investigated the dynamics of these two loops around the ALR2 active site over 200 ns MD simulations by computing the RMSF for the loop structures to ascertain the observed visual changes. Figure 10 shows the RMSF plots of loops A and B. Of course, the loops RMSF showed fluctuations. The fluctuations displayed by the Apo and ALR2-8b are significant, while ALR2-6b and ALR2-7b changes were almost the same. Indeed, these dynamic loop analyses further confirm the ligands' inhibition of ALR2 activity. Interestingly, loop A snapshots showed that loop A moved to cover the ALR2 active site, especially in the ALR2-8b complex as the simulation progressed.
On the other hand, loop B displayed twisting behavior throughout the 200 ns simulations. Such structural changes are critical for the functional activity of proteins, including aldose reductase. For instance, the opening and closing of loop A could impact the ligand-binding terrain of the ALR2 protein to enhance the binding affinity of the inhibitors.
Therefore, we further investigated the dynamics of these two loops around the ALR2 active site over 200 ns MD simulations by computing the RMSF for the loop structures to ascertain the observed visual changes. Figure 10 shows the RMSF plots of loops A and B. Of course, the loops RMSF showed fluctuations. The fluctuations displayed by the Apo and ALR2-8b are significant, while ALR2-6b and ALR2-7b changes were almost the same. Indeed, these dynamic loop analyses further confirm the ligands' inhibition of ALR2 activity. The molecular dynamic (MD) simulations technique and solvent accessible surface area (SASA) have been explored to study proteins' and solvents' conformational movements [63]. The comparison of the MD simulations at low and high pressure gives rise to microscopic pressure-induced denaturation, findings that corroborate experimental data. We performed 200 ns simulations to observe the unfolding conformational events in the Apo and inhibitor-bound systems by investigating the SASA. The SASA analysis of ALR2 systems' dynamics and motions further corroborates the above findings. The average SASA (Å) values for the Apo, ALR2-8b, ALR2-6b, and ALR2-7b were 13,552.9, 13,054.3, The molecular dynamic (MD) simulations technique and solvent accessible surface area (SASA) have been explored to study proteins' and solvents' conformational movements [63]. The comparison of the MD simulations at low and high pressure gives rise to microscopic pressure-induced denaturation, findings that corroborate experimental data. We performed 200 ns simulations to observe the unfolding conformational events in the Apo and inhibitor-bound systems by investigating the SASA. The SASA analysis of ALR2 systems' dynamics and motions further corroborates the above findings. The average SASA (Å) values for the Apo, ALR2-8b, ALR2-6b, and ALR2-7b were 13,552.9, 13,054.3, 12,277.8, and 13,037.6, respectively. The SASA difference of the unbound and bound with the highest SASA value was about 500 Å. These results ( Figure 11) revealed significantly reduced solvent values in the bound systems, with ALR2-6b having the lowest SASA. The inhibitors induced reorientation of the proximal residues to the active site and residues that make up the site. The more the time evolution of simulation and the more burying of the inhibitors, the more likely the solvent reduction was at the protein-ligand interface. Moreover, the SASA trajectories of aldose reductase for polar residues include charged groups, non-polar species, and whole protein. Hence, we could infer that the hydrophobic residues in ALR2 surfaces are more exposed than the polar or charged residues, accounting for the decrease in the solvent accessibility of the bound ALR2 systems.
Biomolecules 2022, 12, x FOR PEER REVIEW 24 of 28 12,277.8, and 13,037.6, respectively. The SASA difference of the unbound and bound with the highest SASA value was about 500 Å. These results ( Figure 11) revealed significantly reduced solvent values in the bound systems, with ALR2-6b having the lowest SASA. The inhibitors induced reorientation of the proximal residues to the active site and residues that make up the site. The more the time evolution of simulation and the more burying of the inhibitors, the more likely the solvent reduction was at the protein-ligand interface. Moreover, the SASA trajectories of aldose reductase for polar residues include charged groups, non-polar species, and whole protein. Hence, we could infer that the hydrophobic residues in ALR2 surfaces are more exposed than the polar or charged residues, accounting for the decrease in the solvent accessibility of the bound ALR2 systems. Figure 11. Solvent accessible surface area (SASA) plots of unbound protein (black, working code AR-Apo), and ALR2-8b (red, working code AR-FB6), ALR2-6b (green, working code AR-FB7), and ALR2-7b (blue, working code AR-FB8) systems.

Conclusions
This work described a novel series of (E)-benzaldehyde O-(benzyl)-oximes as ALR2 inhibitors possessing ancillary anti-oxidant properties. Inspired by natural compounds, they are characterized by an (E)-benzaldehyde O-benzyl oxime core structure, which has already proven to be a key fragment for the obtainment of effective enzymatic inhibitors [64,65]. A comprehensive polyhydroxy substitution pattern on the benzaldehyde residue has been investigated, in a combination with suitable substituents on the benzyl moiety. As for the ALR2 inhibitory efficacy, the 2,3,4-trihydroxy substituents turned out to be the most promising, once combined with a 3-or 4-methoxy group on the benzyl ring. In fact, compound (E)-2,3,4-trihydroxybenzaldehyde O-(3-methoxybenzyl) oxime (7b) and (E)-2,3,4-trihydroxybenzaldehyde O-(4-methoxybenzyl) oxime (8b) displayed favorable ALR2 inhibitory activity, corroborated by both molecular docking and molecular dynamics studies. The ligand-bound ALR2 systems showed insignificant average RMSD and RMSF, implying relatively stable simulations and thoroughly equilibrated bound systems. The PRED analysis revealed the detailed residue energy contributions for each Figure 11. Solvent accessible surface area (SASA) plots of unbound protein (black, working code AR-Apo), and ALR2-8b (red, working code AR-FB6), ALR2-6b (green, working code AR-FB7), and ALR2-7b (blue, working code AR-FB8) systems.

Conclusions
This work described a novel series of (E)-benzaldehyde O-(benzyl)-oximes as ALR2 inhibitors possessing ancillary anti-oxidant properties. Inspired by natural compounds, they are characterized by an (E)-benzaldehyde O-benzyl oxime core structure, which has already proven to be a key fragment for the obtainment of effective enzymatic inhibitors [64,65]. A comprehensive polyhydroxy substitution pattern on the benzaldehyde residue has been investigated, in a combination with suitable substituents on the benzyl moiety. As for the ALR2 inhibitory efficacy, the 2,3,4-trihydroxy substituents turned out to be the most promising, once combined with a 3-or 4-methoxy group on the benzyl ring. In fact, compound (E)-2,3,4-trihydroxybenzaldehyde O-(3-methoxybenzyl) oxime (7b) and (E)-2,3,4-trihydroxybenzaldehyde O-(4-methoxybenzyl) oxime (8b) displayed favorable ALR2 inhibitory activity, corroborated by both molecular docking and molecular dynamics studies. The ligand-bound ALR2 systems showed insignificant average RMSD and RMSF, implying relatively stable simulations and thoroughly equilibrated bound systems. The PRED analysis revealed the detailed residue energy contributions for each complex. As for the anti-oxidant efficacy, both the 2,3,4-trihydroxy and the 3,4-dihydroxy substitution patterns conferred a remarkable activity, which proved to be well above the one shown by natural compounds such as resveratrol and pterostilbene, used as the reference compounds. Results obtained confirm the relevance of hydroxy groups to give anti-oxidant properties to heterocyclic compounds, as already demonstrated by several authors for similar scaffolds [66,67]. Although far from being conclusive, results obtained with the novel series demonstrated the usefulness of the proposed scaffold for the obtainment of active compounds, suggesting the (E)-benzaldehyde O-(benzyl)-oximes as prototypes of novel anti-oxidant ARIs deserving to be further investigated.