A Competition between Hydrogen, Stacking, and Halogen Bonding in N-(4-((3-Methyl-1,4-dioxo-1,4-dihydronaphthalen-2-yl)selanyl)phenyl)acetamide: Structure, Hirshfeld Surface Analysis, 3D Energy Framework Approach, and DFT Calculation

N-(4-((3-Methyl-1,4-dioxo-1,4-dihydronaphthalen-2-yl)selanyl)phenyl)acetamide (5), C19H15NO3Se, was prepared in two steps from 4,4′-diselanediyldianiline (3) via reduction and subsequent nucleophilic reaction with 2-methyl-3-bromo-1,4-naphthalenedione, followed by acetylation with acetic anhydride. The cytotoxicity was estimated against 158N and 158JP oligodendrocytes and the redox profile was also evaluated using different in vitro assays. The technique of single-crystal X-ray diffraction is used to confirm the structure of compound 5. The enantiopure 5 crystallizes in space group P21 with Flack parameter 0.017 (8), exhibiting a chiral layered absolute structure. Molecular structural studies showed that the crystal structure is foremost stabilized by N-H···O and relatively weak C-H···O contacts between molecules, and additionally stabilized by weak C-H···π and Se···N interactions. Hirshfeld surface analysis is used to quantitatively investigate the noncovalent interactions that stabilize crystal packing. Framework energy diagrams were used to graphically represent the stabilizing interaction energies for crystal packing. The analysis of the energy framework shows that the interactions energies of and C-H···π and C-O···π are primarily dispersive and are the crystal’s main important forces. Density functional theory (DFT) calculations were used to determine the compound’s stability, chemical reactivity, and other parameters by determining the HOMO-LUMO energy differences. The determination of its optimized surface of the molecular electrostatic potential (MEP) was also carried out. This study was conducted to demonstrate both the electron-rich and electron-poor sites.

Additionally, we are not only interested in the synthesis and biological evaluation of Se-based quinones but also in the investigation of their crystal structure elucidation. Herein, we report the crystal structure of compound 5. The cytotoxicity was estimated against 158 N and 158 JP oligodendrocytes and the redox profile was also evaluated using different in vitro assays.
In crystal packing, the nature of intermolecular interactions was characterized through the Hirshfeld surface analysis description. In addition, DFT calculations were employed to optimize the structure of compound 5 in its isolated state. Moreover, Frontier molecular orbitals (FMO) and MEP mapping studies were achieved to identify global reactivity descriptors and the complementary interacting sites in compound 5.

Crystal Structure Measurement
X-ray data were estimated at 115 K from a single crystal of compound 5 on a Bruker diffractometer (APEX-II CCD) using a single wavelength (λ = 0.71073). X-ray MoK α . SAINT [59,60] and SADABS [61] were used to reveal the multi-scan absorption correction and Data reduction. Olex2 [62] was used to solve the structure together with other programs, such as the SHELXT. Additionally, refinement of the structure was implemented with the 2018/3 package of SHELXL [63]. The hydrogens were fixed in their computed places and a riding model was used for the refinement with hydrogen lengths of 0.88 Å (NH) and 0.95-0.98 Å (CH). Displacement of the isotropic parameters was placed to 1.2 (CH, NH) or 1.5 (CH3) times U eq of the parent atom. With riding coordinates, aromatic/amide H has been refined as N(H), C5(H5), C6(H6), C7(H7), C8(H8), C13(H13), C14(H14), C16(H16), and C17(H17). Idealized CH3 was refined into rotating group C11(H11A,H11B,H11C) and C19(H19A,H19B,H19C). The refinement of the twin data was made with a 0.983(8) scale factor, since this resolved compound has an absolute structure. The proof for the chirality is also given from the Flack parameter [64], which is 0.017 (8). PLATON was used to check the results of the X-ray analysis [65]. Mercury 4.0 software was utilized for the single-crystal graphical representation. The relevant crystallographic information for compound 5 is tabulated in Table 1 [66].

Theoretical Calculations
The Gaussian 09 program was used for all the calculations presented in this work [67]. First, the geometries of the compound considered were entirely optimized at the 6-311+G(d,p) level [68]. In the literature, the B3LYP functional is widely used and leads to dependable results regarding the organic compounds' ground state properties [69]. The calculations for normal modes of vibration were performed using optimized geometries, which led to real frequencies that indicated that these geometries were minimums on the potential energy surfaces. We also used the energy quantities to research and confirm the stability of these compounds. The results obtained were compared with the experimental data. Empirical corrections for dispersion were evaluated using DFT-D corrections [70], which necessarily deepens the interaction well and is an essential correction to the DFT energy of this compound.

Design and Synthesis of Compound 5
Selenium dioxide and malonitrile are being used as an efficient selenocyanating medium for several aromatic hydrocarbons, such as aniline [59]. In this case, the corresponding 4-selenocyanatoaniline (2) was obtained with a good yield (88%). The 4,4 -Diselanediyldianiline (3) was prepared with an 82% yield via the hydrolysis of 4-selenocyanatoaniline (2) using strong alkali (e.g., NaOH). The 4,4 -Diselanediyldianiline (3) was used for the synthesis of diverse organoselenium compounds by the reduction of the diselenide bond followed by a nucleophilic reaction with different halogenated compounds. Within this context, the reduction of diselenide compound 3 with NaBH 4 and the succeeding treatment with 2-methyl-3-bromo-l,4-naphthoquinone in a two-phase solvent (ethyl acetate and water 1:1) system and 5% mol Aliquat 336 as the phase transfer agent afforded compound 4 an 89% yield in Scheme 1. The reaction of compound 4 with acetic anhydride furnished compound 5 with a 75% yield (Scheme 1) [54,58].
The levels of ROS were monitored with the H2-DCFDA and DHE assays using flow cytometry. Briefly, 158 N cells were mixed with serial concentrations (0-50 µM) of compound 4 and compound 5 using Vitamin E as the control.
In the H2-DCFDA assay, compound 4 showed a pronounced pro-oxidant activity. On the other hand, compound 5 was able to decrease the H2-DCF intensity, thus diminishing the ROS level at 10 µM compared to vitamin E. In the DHE assay, compound 4 was able to promote the production of O 2 − in 158 N cells [54,58]. Moreover, compound 4 and compound 5 exhibited moderate GPx-like activity. Additionally, bleomycin-induced DNA damage, DPPH, and ABTS in vitro assays were employed to estimate compounds 4 and 5's radical scavenging activities. Compounds 4 and 5 manifested good pro-oxidant properties (up to 60% compared to vitamin C). Interestingly, these data were in good agreement with that obtained from the DHE and H2DCF assays [54,58].

Analysis of the Molecular Packing
Compound 5 crystallizes in the monoclinic system with the chiral space group P2 1 . Its asymmetric unit (Figure 2) is composed of one chiral molecule. The Flack parameter of the absolute configuration of compound 5 is found + 0.017 (8), indicating that the configuration is R. The FLACK parameter with a value of 0 indicates the correct structure and 1 indicates the inverted structure. The main body of compound 5's structure consists of two nearly coplanar fused rings, a benzene ring (A) and phenyl ring (B), and another phenyl ring (C) attached to the Se atom-see Figure 2. The two planes composed by (A, B) and C rings are twisted to each other by 65.35 • . The C1-Se-C12 angle is 97.36 (11) Å (Table 2), which is like the values reported in the literature for the three related compounds [71][72][73]. The Se-C1 and Se-C12 bond lengths are found to be 1.916(3) Å and 1.924(3) Å, respectively, and may be regarded as normal (Table 2) [74,75]. The PLATON analysis suggests that these molecules are bonded together by several types of noncovalent interactions, such as hydrogen bonds (C-H···O, N-H···O, C-H···π), C-O···π, and unusual short contacts such as Se···N contacts. The C(19)-H(19A)···O(1) and N-H···O(3) hydrogen bonds connect molecules into double layers that stack along the a and c axis (Figure 3a,b, Table 3

Hirshfeld Surface Analysis
The idea of Hirshfeld surface (HS) results come from a desire to interpret the space engaged by a molecule in a crystal to divide the crystal's electron density into molecular fragments of electron densities. The nature of intermolecular interactions can be quantified using an HS investigation. CrystalExplorer17.5 was used to perform the HS calculations [76]. Figure 5 shows the HS plotted over d norm , a property with values ranging from −0.4854 (red) to 1.2838 (blue) a.u. [77]. The intermolecular contacts are visualized using color-coding. Intermolecular contacts with distances less than, equal to, and greater than van der Waal radii are shown in red, white, and blue spots, respectively. The presence of the N-H···O interaction in crystal packing is indicated by the intense red regions around the NH and O atoms of the carbonyl groups of the cation, whereas the C-H···O interaction is indicated by the pale-red spot-on HS. Additionally, two-dimensional (2D) finger-plots are employed to give quantitative data regarding the type and nature of the intermolecular interactions that are candidates in the crystal packing [78]. 2D Fingerplots are computed for each interatomic contact and the overall interactions. The reciprocal contact of individually interatomic contact is also comprised in the computation of each interatomic contact. The 2D Fingerplots for the overall interactions is shown in Figure 6a. Interatomic contacts that contribute significantly to crystal packing are represented by spikes in this graph. H···H, C···H, and O···H are the interatomic contacts that contribute the most to crystal packing, with percentage contributions of 40% (Figure 5b), 24. 7% (Figure 6c), and 20.7% (Figure 6d), respectively. Se···H, O···C, Se···C, C···C, N···H, Se···N, C···N, and Se···O are the other interatomic contacts that play a smaller role in crystal packing, with percentage contributions of 4.8%, 4.3%, 1.9%, 1.5%, 0.9%, 0.8%, 0.2% and 0.2%, respectively (Figure 7).

Energy Framework
Energy frameworks could be simulated using Crystal Explorer17 software. The latter is a graphical representation of individual energy components as cylinders connecting the centroids of interacting molecular pairs, with E ele , E dis , and E tot color-coded in red, green, and blue, respectively, and the interaction energy's magnitude proportional to the radius of the corresponding cylinders. Therefore, the intermolecular interaction energies for compound 5 are calculated using the CE-HF . . . HF/3-21G energy model presented in Crystal-Explorer [76]. 4b-between the central molecule (grey molecule) and the x, y, z symmetry-related molecule (red)-are, by far, the strongest interaction among near neighbors with an interaction energy of −83.6 kJ mol −1 . Hydrogen bonds with C-H···O interactions energies (shamrock green) have an energy value of −35.2 kJ/mol. The total energy diagram (Figure 9c) showed a strong resemblance to the dispersion energy frameworks (Figures 9b), indicating that they play a significant role in the total forces in crystal packing.

Geometric Structures
The DFT structure optimization of compound 5 was performed using the B3LYP/6-311+ G(d,p) level of theory in the gaseous phase ( Figure 10). Calculations were performed to estimate the stability of the studied compound, the reactivity chemical by determining the HOMO-LUMO energy differences, the ionization potential (I), the affinity electronics (A), the electrophilicity index (ω), the chemical potential (µ), and the hardness (η) and the softness (S) of compound 5. The determination of their optimized potential surface molecular electrostatic (MEP), to characterize the influence of the different substitution groups, has also been carried out. Finally, a study was carried out to demonstrate the rich sites and the poor sites of electrons. The visualized HOMO and LUMO of compound 5 are shown in Figure 11.  The results are presented in different, numbered tables. In addition, it seems interesting to us to start with the results obtained from the geometric parameters molecules studied. The results obtained concerning this calculation are grouped in Table 4 according to the following numbering. Compare these theoretically calculated geometric parameters with experimental data from a crystallographic analysis [79], knowing that the relative deviation ∆ for a geometric parameter X is expressed as a percentage by the relation: X theo : the theoretical value of the quantity X. X exp : the experimental value of the quantity X. Analysis of the results reported in Table 4 shows that the average deviation of the distances and angles obtained by the DFT method is on the order of less than (3%) overall, this lets us say that the different theoretical results obtained are in very good agreement with those obtained experimentally by crystallographic analysis.

Surfaces with Molecular Electrostatic Potential (MEP)
Molecular electrostatic potentials (MEP) are amongst the descriptors to correlate physicochemical properties with molecular structure. Therefore, MEP can differentiate sites of electrophilic attack (electron-rich areas) from nucleophilic attack (electron-poor areas). MEP are useful tools for predicting biological processes, especially sites of electrophile or nucleophile attack by a reagent. In Figure 12, the different amplitudes of the MEP are shown by different colors as follows: red < orange < yellow < green < blue [80]. The region close to the N atoms of the pyridine and imidazole ring is characterized by red and yellow and exhibits high electron density (DE). It, therefore, is attributed to a negative potential. Thus, the electrophilic attack is more likely preferred at these sites.

Study of Overall Global Reactivities
The HOMO-LUMO gap is an essential factor in the chemistry of quantum; it grants the characterization of molecular stability. Furthermore, frontier molecular orbitals are essential in the prediction of molecular reactivity [81,82]. Molecules with a small HOMO-LUMO energy gap could be distorted to increase this gap. The results of the descriptors calculated by DFT (Table 5) are the HOMO and LUMO energies, energy gap (∆E), electronegativity (χ), the hardness (η), the chemical potential (µ), the global softness (S), the global electrophilicity index (ω), and the softness (σ), respectively, and given by the Equations (1)-(8) [83]: Vertical Electron Affinity (VEA) = E neutral -E anion The compound studied has a hard character because it carries a high value of chemical hardness. On the other hand, the compound carries the low value of the chemical potential followed by the high value of the electrophilicity index, hence the compound promotes its electrophilic behavior. From a biological point of view, selenoquinone compounds (e.g., compounds 4 and 5) demonstrated potential cytotoxicity. In theory, the Se and quinone centers complement each other, i.e., the quinone redox cycling process leads to the formation of superoxide radicals [52][53][54][55]. The latter, in turn, activates the Se center to further induce oxidation of redox-sensitive cellular compartments (e.g., enzymes and proteins) and, thus, causes cell death [52][53][54][55]. These results were in agreement with the DFT calculations in which the HOMO is mainly located on the nucleophilic centers (e.g., quinone and Se) of compound 5, however, it is still too early to reach conclusions at this point. An unambiguous QSAR requires a diverse and huge number of compounds to get a clear understanding of the relation between the DFT calculation results and potential cytotoxicity. Ultimately, this justifies the requirements of more studies, such as additional cytotoxicity experiments and computational calculations to identify the possible parameters and the underlying structural requirements leading to cytotoxicity.

Conclusions
The structure of the novel organoselenium compound 5 was confirmed through employing the single-crystal X-ray diffraction technique. The layers that form parallel to a crystallographic axis are produced by H-bonding types N-H···O and C-H···O. The presence of relatively weak C···H, Se···N, and C···O interactions contributes to further crystal packing stabilization. The H···H interatomic contact, which contributes the most to crystal packing with a percentage contribution of 40%, is revealed by HS analysis for further noncovalent interaction exploration. The crystal packing's three-dimensional topology was analyzed and visualized using energy framework calculations. The total interaction energy in the crystal packing is dispersive, according to the results.
Compound 5 exhibited potential cytotoxicity against the 158 N and 158 JP oligodendrocytes and decreased the ROS level at 10 µM compared to vitamin E. It also showed moderate GPx-like activity.
Our work is a theoretical contribution to the study of the structure; chemical reactivity of the studied compound, which plays an important role in biology; pharmacology; and synthesis. We carried out two parts of this research work. Excellent agreement between the calculated and experimental data was achieved. It is in this context that Frontier Molecular Orbital Theory, and, ultimately, Conceptual DFT emerged and developed.

Data Availability Statement:
The raw/processed data generated in this work are available upon request from the corresponding author.