Quantum Computational, Spectroscopic (FT-IR, FT-Raman, NMR, and UV–Vis) Hirshfeld Surface and Molecular Docking-Dynamics Studies on 5-Hydroxymethyluracil (Monomer and Trimer)

For many decades, uracil has been an antineoplastic agent used in combination with tegafur to treat various human cancers, including breast, prostate, and liver cancer. Therefore, it is necessary to explore the molecular features of uracil and its derivatives. Herein, the molecule’s 5-hydroxymethyluracil has been thoroughly characterized by NMR, UV–Vis, and FT-IR spectroscopy by means of experimental and theoretical analysis. Density functional theory (DFT) using the B3LYP method at 6-311++G(d,p) was computed to achieve the optimized geometric parameters of the molecule in the ground state. For further investigation and computation of the NLO, NBO, NHO analysis, and FMO, the improved geometrical parameters were utilized. The potential energy distribution was used to allocate the vibrational frequencies using the VEDA 4 program. The NBO study determined the relationship between the donor and acceptor. The molecule’s charge distribution and reactive regions were highlighted using the MEP and Fukui functions. Maps of the hole and electron density distribution in the excited state were generated using the TD-DFT method and PCM solvent model in order to reveal electronic characteristics. The energies and diagrams for the lowest unoccupied molecular orbital (LUMO) and the highest occupied molecular orbital (HOMO) were also provided. The HOMO–LUMO band gap estimated the charge transport within the molecule. When examining the intermolecular interactions in 5-HMU, Hirshfeld surface analysis was used, and fingerprint plots were also produced. The molecular docking investigation involved docking 5-HMU with six different protein receptors. Molecular dynamic simulation has given a better idea of the binding of the ligand with protein.


Introduction
One of the important derivatives of methyl uracil is 5-hydroxymethyluracil (5-HMU). The heterocyclic category includes uracil and its derivatives. Since some of them are the primary building blocks of DNA and RNA such as pyrimidine, cytosine, and uracil, and their derivatives are of prime value. Researchers in the fields of spectroscopy, medicine,

Introduction
One of the important derivatives of methyl uracil is 5-hydroxymethyluracil (5-HMU). The heterocyclic category includes uracil and its derivatives. Since some of them are the primary building blocks of DNA and RNA such as pyrimidine, cytosine, and uracil, and their derivatives are of prime value. Researchers in the fields of spectroscopy, medicine, and industry have all paid close attention to pyridine-derived biomolecules because of their biological activities [1][2][3][4][5][6].
Reactive oxygen species are created by cells, which damage biomolecules and are crucial to the development of cancers and aging. The creation of 8-oxo-7,8-dihydroguanine, 5-HMU, 5-hydroxyuracil (5-HU), and other major forms occurs as a result of oxidative DNA damage when nucleotides are hydroxylated. As indicators of oxidative DNA damage, hydroxylated DNA adducts have been studied [7]. Although the biological effects of 5-HMU in humans are not completely known, certain glycosylase activities that remove 5-HMU from DNA have been found [8]. During the T-transition mutations, the conversion of 5-methylcytosine to 5-HMU is frequently detected in cancers. High level of 5-HMU formation in tissue DNA has been linked to human disorders and, hence, is regarded as a biomarker for breast cancer and oxidative DNA damage [9,10].
The pharmacological properties of uracil and its derivatives (Scheme 1) have served as the motivation for recent spectroscopic studies [11]. The formation of 5-HMU takes place by the attack of carbonyl oxygen on the DNA's methyl group. Hydroxy radicals may bring on different modifications to DNA bases. A big aspect of the development of matters is to explore the hydrogen bonding interaction between the 5-HMU and base pairs of DNA [11]. The excited state processes (internal conversion and intersystem crossover) in the fields of intramolecular charge transfer have significantly contributed to analyzing their kinetics of hydrogen bonding resulting from vibrations of the hydrogen donor and acceptor groups [12][13][14]. Scheme 1. Uracil and its derivatives.
As far as researchers know, there has not been any literature-reported computational vibrational research on the 5-HMU's structure. In this study, DFT/ B3LYP theoretical methods with the "6-311++G(d,p) basis set" were employed to gain a comprehensive description of the molecular structure and vibrational pattern of 5-HMU along with experimental analysis (FT-IR, NMR and UV-Vis spectroscopy). Finally, the findings were compared with the theoretical outcomes. The chemicals ampyra, amine dimedone, sulfanilic acid, and sulfathiazole were the focus of similar experiments that we also conducted [15][16][17][18][19].

Optimized Geometry
The optimized geometry of 5-HMU was acquired using the B3LYP methods with the 6-311++G (d, p) basis set, and their comparisons have been made with the experimental analysis (Table 1A and B). The optimized molecular geometry of 5-HMU is shown in Figure 1A and the crystal structure analysis with the diamond program for a trimer ( Figure 1B) and two dimers ( Figure 1C and D), which was in close resemblance with the structural parameters of the refcode AYEBIP [20,21]. Cagrı et al. have calculated the bond Scheme 1. Uracil and its derivatives.
As far as researchers know, there has not been any literature-reported computational vibrational research on the 5-HMU's structure. In this study, DFT/ B3LYP theoretical methods with the "6-311++G(d,p) basis set" were employed to gain a comprehensive description of the molecular structure and vibrational pattern of 5-HMU along with experimental analysis (FT-IR, NMR and UV-Vis spectroscopy). Finally, the findings were compared with the theoretical outcomes. The chemicals ampyra, amine dimedone, sulfanilic acid, and sulfathiazole were the focus of similar experiments that we also conducted [15][16][17][18][19].

Optimized Geometry
The optimized geometry of 5-HMU was acquired using the B3LYP methods with the 6-311++G (d, p) basis set, and their comparisons have been made with the experimental analysis (Table 1A,B). The optimized molecular geometry of 5-HMU is shown in Figure 1A and the crystal structure analysis with the diamond program for a trimer ( Figure 1B) and two dimers ( Figure 1C,D), which was in close resemblance with the structural parameters of the refcode AYEBIP [20,21] 51, and between 1.08 and 1.1 Å, respectively, which were close to the experimental values of uracil [22]. Each 5-HMU molecule is linked to each other via a hydrogen bond to form a trimer unit. In this case, the 2-N-H...O bonds hold within the 5-HMU trimer. The carbonyl and hydroxyl groups drive forming a trimer between monomer complexes in their close proximity. Figure 1 depicts the intermolecular hydrogen bond N-H...O formation based on the optimized structure (bond angles and distances) with the least amount of energy (Table 1B). Because of the inter-molecular interactions, the energy (E) of the trimer was found to be −1588.58 a.u. with N-H . . . O bond lengths of 1.94 and 1.92 Å (Supplementary Materials, Figure S1). The bond angle of 174.17 • and 174.54 • in the trimer unit is almost linear. Two independent dimers exhibited by 5-HMU are shown in Figure 1C,D, and the parameters were summarized in Supplementary Materials; Tables S1 and S2. proximity. Figure 1 depicts the intermolecular hydrogen bond N-H...O formation based on the optimized structure (bond angles and distances) with the least amount of energy (Table 1B). Because of the inter-molecular interactions, the energy (E) of the trimer was found to be −1588.58 a.u. with N-H…O bond lengths of 1.94 and 1.92 Å (Supplementary Materials, Figure S1). The bond angle of 174.17º and 174.54º in the trimer unit is almost linear. Two independent dimers exhibited by 5-HMU are shown in Figures 1C and D, and the parameters were summarized in Supplementary materials; Table S1 and S2.

Vibrational Spectral Analysis
IR spectroscopy may be used to identify chemical groups, and it provides microscopic data on their conformations. IR spectroscopy is a vital method to locate the actual binding regions. Both observed and calculated FT-IR and FT-Raman spectra of 5-HMU were compared in Figure 2A,B. 5-HMU is composed of 14 atoms, resulting in 36 fundamental vibrations modes with C1 point group symmetry, all of which are active in FT-IR and FT-Raman. The estimated frequencies (scaled and unscaled) and PED assignments for 5-HMU on the gaseous phase were assigned by the VEDA software (Table 2).
Molecules 2023, 28, x FOR PEER REVIEW 5 of 22 By using the B3LYP method and 6-311++G(d,p) method, the bands at 3495 and 3455 cm −1 correspond to N-H stretching frequencies, consistent with the experimental value at 3497 cm −1 [23]. These pure modes provide 50% and 50% PED, respectively.
The stretching vibration in the FT-Raman spectrum observed at 3068, 2980, and 2917 cm −1 for C-H was noticed to be in the range 3083, 2954, and 2919 cm −1 in the FT-IR by calculated methods, and their PED assignments contribute to 99%, 94%, and 95% ( Table  2).

C-N, H-C-C-C, and C-C Vibrations
The C-N stretching frequency occurs in a composite, which permits the vibrational of various bands [24] In benzamide, the designated C-N stretching band at 1368 cm -1 [25]. The stretching vibrations of C-N with a maximum PED assignment of 34% are allocated to the theoretically scaled wavenumbers in FT-IR estimated at 1444, 1438, 1369, 1354, 1166, 1142, 937, 746, and 525 cm −1 . The FT-Raman band of C-N vibrations was displayed at 1432, 1365, 1175, 1148, 1090, and 514 cm −1 .
The H-C-C-C torsion vibrations were observed in this investigation, with FT-IR spectra maxima at 1626, 1369, 1290, and 1142 cm −1 .
The bands in the range of 1400 and 1650 cm -1 have been ascribed to vibrations of the C-C heterocyclic ring [26]. In the range of 1300-1000 cm -1 , stretching vibrations (C-C) are predictable [27]. The C-C stretching vibrations have theoretical frequencies of 1626, 1290, 1166, and 1091 cm −1 . PED contributes 65%, 15%, 23%, and 11% to various modalities. In the present study, the FT-Raman spectrum shows the stretching vibration of C-C was observed at 1269 and 791 cm −1 . It demonstrates that theoretical values and experimental data are in good accord.

N-H and C-H Vibrations
By using the B3LYP method and 6-311++G(d,p) method, the bands at 3495 and 3455 cm −1 correspond to N-H stretching frequencies, consistent with the experimental value at 3497 cm −1 [23]. These pure modes provide 50% and 50% PED, respectively.
The stretching vibration in the FT-Raman spectrum observed at 3068, 2980, and 2917 cm −1 for C-H was noticed to be in the range 3083, 2954, and 2919 cm −1 in the FT-IR by calculated methods, and their PED assignments contribute to 99%, 94%, and 95% (Table 2).

C-N, H-C-C-C, and C-C Vibrations
The C-N stretching frequency occurs in a composite, which permits the vibrational of various bands [24] In benzamide, the designated C-N stretching band at 1368 cm −1 [25]. The stretching vibrations of C-N with a maximum PED assignment of 34% are allocated to the theoretically scaled wavenumbers in FT-IR estimated at 1444, 1438, 1369, 1354, 1166, 1142, 937, 746, and 525 cm −1 . The FT-Raman band of C-N vibrations was displayed at 1432, 1365, 1175, 1148, 1090, and 514 cm −1 .
The H-C-C-C torsion vibrations were observed in this investigation, with FT-IR spectra maxima at 1626, 1369, 1290, and 1142 cm −1 .
The bands in the range of 1400 and 1650 cm −1 have been ascribed to vibrations of the C-C heterocyclic ring [26]. In the range of 1300-1000 cm −1 , stretching vibrations (C-C) are predictable [27]. The C-C stretching vibrations have theoretical frequencies of 1626, 1290, 1166, and 1091 cm −1 . PED contributes 65%, 15%, 23%, and 11% to various modalities. In the present study, the FT-Raman spectrum shows the stretching vibration of C-C was observed at 1269 and 791 cm −1 . It demonstrates that theoretical values and experimental data are in good accord.

NMR Analysis
NMR is a reliable method for obtaining detailed information regarding the structure and chemical environment of organic compounds [28]. Hence, the overall geometry optimization of 5-HMU was analyzed employing the B3LYP method and 6-311++G(d,p) basis set, and with the help of GIAO methodology using the aug-cc-pV5Z def2/J basis set, their chemical shifts value was calculated. The input files have been created with the atomistica online tool [29]. Table 3 shows the actual and estimated chemical shifts for 1 H and 13 C NMR. In DMSO solvent, the experimental NMR of 5-HMU were recorded and displayed in Figure 3A,B. In 1 H NMR, the peaks for hydrogen attached to nitrogen (H11 and H13) were observed at 10.66 and 7.20 ppm and theoretically appeared at 7.06 ppm and 6.75 ppm, respectively. In 13 C NMR, the carbonyl carbons (C1 and C5) were observed in the area of 164. 31

Hirshfeld Surface (HS) Analysis
HS analysis was established to divide the crystal's electron density into covalently electron density pieces. The HS is plotted over d norm , a property that has values between −0.6245 (red) and 0.9768 (blue) a.u. [30][31][32][33][34] along with d i , d e , shape index, curvedness, and fragment patch, respectively ( Figure 4). The values of d i ranged from 0.7259 to 2.3342, of d e from 0.7281 to 2.2269, of shape index from −0.9922 to 0.9975, of curvedness from −4.043 to 0.306, and of fragment patch from 0 to 13 (Supplementary Material, Table S3). Intermolecular interactions could be seen using color labeling. Spots in red, white, and blue represent intermolecular interactions that occur at distances that are less than, equal to, and greater than the van der Waal radii, respectively. set, their chemical shifts value was calculated. The input files have been created with the atomistica online tool [29]. Table 3 shows the actual and estimated chemical shifts for 1 H and 13 C NMR. In DMSO solvent, the experimental NMR of 5-HMU were recorded and displayed in Figure 3 (A and B). In 1 H NMR, the peaks for hydrogen attached to nitrogen (H11 and H13) were observed at 10.66 and 7.20 ppm and theoretically appeared at 7.06 ppm and 6.75 ppm, respectively. In 13 C NMR, the carbonyl carbons (C1 and C5) were observed in the area of 164. 31

Hirshfeld Surface (HS) Analysis
HS analysis was established to divide the crystal's electron density into covalently electron density pieces. The HS is plotted over dnorm, a property that has values between −0.6245 (red) and 0.9768 (blue) a.u. [30][31][32][33][34] along with di, de, shape index, curvedness, and fragment patch, respectively ( Figure 4). The values of di ranged from 0.7259 to 2.3342, of de from 0.7281 to 2.2269, of shape index from −0.9922 to 0.9975, of curvedness from −4.043 to 0.306, and of fragment patch from 0 to 13 (Supplementary Material, Table  S3). Intermolecular interactions could be seen using color labeling. Spots in red, white, and blue represent intermolecular interactions that occur at distances that are less than, equal to, and greater than the van der Waal radii, respectively. The HS and its inside closest nucleus are separated by a distance of di, whereas the closest nucleus outside is separated by a distance of de. The type of weak interactions that the molecules in the crystal packing experience are also quantified using 2-D fingerprint designs [35][36][37]. When calculating individual interatomic interactions, the mu- The HS and its inside closest nucleus are separated by a distance of d i , whereas the closest nucleus outside is separated by a distance of d e . The type of weak interactions that the molecules in the crystal packing experience are also quantified using 2-D fingerprint designs [35][36][37]. When calculating individual interatomic interactions, the mutual contact for each interatomic contact is also taken into account. The spikes of 2D FP plot show the interatomic interactions that significantly influence crystal packing (Supplementary Material, Figure S2

Molecular Electrostatic Potential Surface (MEPs)
The MEP's crucial points can be utilized to identify both the structural and electronic characteristics of molecules. It also indicates the shape and size and its positive, negative, and neutral potentials of the molecular system with the help of a color gradient. It is a very useful factor for its electrophilic and nucleophilic properties as well as hydrogen bonding interactions [38]. It also holds importance to verify the proof of the drug's reaction as an inhibitor and its reaction toward the electrophilic and nucleophilic attacks. This can also be used to predict the relationship between the physicochemical characteristics and molecular structure of the drug [39] such as (a) molecular recognition, drug-receptor interactions; (b) reactive properties of biomolecules; and (c) chemical carcinogenesis, etc. [40,41].
The negative and positive potential sites are preferred for an electrophilic and nucleophilic attack, respectively, corresponding to red and blue colors, respectively. The zero potential sites are encoded in a green color.
The MEP of 5-HMU and its trimer was calculated using DFT and B3LY/6-311++G(d,p) level of theory and ranged 6.918 × 10 −2 , -6.256 to 6.918 × 10 −2 , 6.256 × 10 −2 a.u, respectively. The 3D plots of monomer and trimer show that the negatively charged potential. Which is a site of the electrophilic attack, is seen near the oxygen atom while the positive potential (blue color) is found near the carbon and hydrogen atoms ( Figure 5A,B). It can be inferred that the reactive sites of the 5-HMU are C=O and nitrogen atoms. The MEP's crucial points can be utilized to identify both the structural and electronic characteristics of molecules. It also indicates the shape and size and its positive, negative, and neutral potentials of the molecular system with the help of a color gradient. It is a very useful factor for its electrophilic and nucleophilic properties as well as hydrogen bonding interactions [38]. It also holds importance to verify the proof of the drug's reaction as an inhibitor and its reaction toward the electrophilic and nucleophilic attacks. This can also be used to predict the relationship between the physicochemical characteristics and molecular structure of the drug [39] such as (a) molecular recognition, drug-receptor interactions; (b) reactive properties of biomolecules; and (c) chemical carcinogenesis, etc [40,41].
The negative and positive potential sites are preferred for an electrophilic and nucleophilic attack, respectively, corresponding to red and blue colors, respectively. The zero potential sites are encoded in a green color.
The MEP of 5-HMU and its trimer was calculated using DFT and B3LY/6-311++G(d,p) level of theory and ranged 6.918 × 10 −2 , -6.256 to 6.918 × 10 −2 , 6.256 × 10 −2 a.u, respectively. The 3D plots of monomer and trimer show that the negatively charged potential. Which is a site of the electrophilic attack, is seen near the oxygen atom while the positive potential (blue color) is found near the carbon and hydrogen atoms ( Figure 5A and B). It can be inferred that the reactive sites of the 5-HMU are C=O and nitrogen atoms.

Electron Localization Function
If an electron pair is localized, it means that there is a high possibility of paired electron system with the opposite spin in a certain area of space, leading to the low probability of interaction with outer electrons of this area. A surplus of kinetic energy density generated as a result of Pauli repulsion serves as the foundation for ELF measurement. The baseline of the electron localization function is formed by the electron pairs caused by covalent bonds in molecular space.
ELF can also reveal some important information regarding reactivity, molecular bonding, and chemical structure [42]. Tiny, localized electrons comprise atomic shells, chemical bonds, and lone pair electrons, which governed the strong Pauli repulsion zone. The graph represents high values in a red color (maximum Pauli repulsion) and lower values in a blue color (no Pauli repulsion). The ELF values can be represented

Electron Localization Function
If an electron pair is localized, it means that there is a high possibility of paired electron system with the opposite spin in a certain area of space, leading to the low probability of interaction with outer electrons of this area. A surplus of kinetic energy density generated as a result of Pauli repulsion serves as the foundation for ELF measurement. The baseline of the electron localization function is formed by the electron pairs caused by covalent bonds in molecular space.
ELF can also reveal some important information regarding reactivity, molecular bonding, and chemical structure [42]. Tiny, localized electrons comprise atomic shells, chemical bonds, and lone pair electrons, which governed the strong Pauli repulsion zone. The graph represents high values in a red color (maximum Pauli repulsion) and lower values in a blue color (no Pauli repulsion). The ELF values can be represented graphically in two and three dimensions as a shaded surface map with projection exposing the hydrogen bonding ( Figure 6). As we can see, the region around the N-H and C-H bonds has a high ELF value (red) and, hence, a greater chance of finding a spin pair electron while the region around the nucleus of the atoms in the rings shows blue rings.

Non-linear Optical (NLO) Analysis
DFT methods employing B3LYP/6-311++G(d,p) basis set have been utilized exte sively to study the NLO characteristics of organic materials [43,44]. For this purpose, t dipole moment (μD), linear polarizability (α), and first-order hyperpolarizability (β were computed; electron correlation, which is dependent on the basis set used, can al the hyperpolarizability value [45][46][47]. Table 4 lists the computed dipole moment as w as the linear polarizability and first-order hyperpolarizability. The values α and β from the Gaussian output files are given in atomic units (a.u.), which were th translated into electrostatic units (esu) using the formulas 1 a.u. = 0.184210 −24 esu α and 1 a.u. = 8.639610 −33 esu for β.
The calculated dipole moment for 5-HMU was found to be 3.5115 Debye, and t highest dipole moment is seen in the μz component. The linear polarizability a first-order hyperpolarizabilities were -0.914 ×10 −23 and 0.3545×10 −30 esu, respective compared to the values of urea which are 0.9771×10 −23 and 0.927 ×10 −30 esu, respectively The NBO technique offers knowledge about intra-and inter-molecular interactio between the bonds, as well as a suitable framework for dealing interactions in each fill and virtual orbital area, as well as charge transfer and conjugative interactions in mol ular systems [48]. The second-order Fock matrix was employed in the NBO analysis

Non-Linear Optical (NLO) Analysis
DFT methods employing B3LYP/6-311++G(d,p) basis set have been utilized extensively to study the NLO characteristics of organic materials [43,44]. For this purpose, the dipole moment (µD), linear polarizability (α), and first-order hyperpolarizability (β tot ) were computed; electron correlation, which is dependent on the basis set used, can alter the hyperpolarizability value [45][46][47]. Table 4 lists the computed dipole moment as well as the linear polarizability and first-order hyperpolarizability. The values of α and β from the Gaussian output files are given in atomic units (a.u.), which were then translated into electrostatic units (e.s.u) using the formulas 1 a.u. = 0.184210 −24 e.s.u for α and 1 a.u. = 8.639610 −33 e.s.u for β. The calculated dipole moment for 5-HMU was found to be 3.5115 Debye, and the highest dipole moment is seen in the µ z component. The linear polarizability and first-order hyperpolarizabilities were −0.914 × 10 −23 and 0.3545 × 10 −30 e.s.u, respectively, compared to the values of urea which are 0.9771 × 10 −23 and 0.927 × 10 −30 e.s.u, respectively.

Natural Bond Orbital (NBO) and Natural Hybrid Orbital (NHO) Analysis
The NBO technique offers knowledge about intra-and inter-molecular interactions between the bonds, as well as a suitable framework for dealing interactions in each filled and virtual orbital area, as well as charge transfer and conjugative interactions in molecular systems [48]. The second-order Fock matrix was employed in the NBO analysis to evaluate the donor-acceptor interactions [49]. An empty non-Lewis orbital replaces the localized NBO of the perfect Lewis structure as a result of the interaction. The stabilization energy E (2) associated with the delocalization i → j is calculated for each donor (i) and acceptor (j) as: 2 , ε σ , ε σ *, and n σ are the Fock matrix element between donor (i) and acceptor (j), the energies of bonding NBO, anti-bonding NBO, and population in the s (donor) orbital, respectively [50]. The more intense the engagement between electron donors and electron acceptors, or the greater the propensity for electron donors to donate to electron acceptors, the stronger the E (2) value, and the more conjugation throughout the system. A stabilizing donor-acceptor interaction takes place when electron density delocalizes between occupied Lewis type (bond or lone pair) NBO orbitals and formally unoccupied (antibonding and Rydberg) non-Lewis NBO orbitals. NBO analysis was analyzed on a Gaussian 09W [51] package at DFT-B3LYP/6-311++G(d,p) theory level (Supplementary Material, Table S4). The weak departure from precisely localized natural Lewis structure that signifies the primary noncovalent effect is designated by the delocalization of σ-σ* electrons. In 5-HMU, the highest stabilization is shown by the delocalization of π (C3-C4) to π*(C5-O7) of 22.5 kJ/mol. The electron contributing from the lone pair(1) of the nitrogen atom (N2) and N6 to the (C1-O8), which offers a good stabilizing energy of 53.48 and 58 kJ/mol, respectively, are another important interaction connected to the resonance in the molecule. Similarly, the electron contributing from lone pair(1) N6 to π*(C5-O7) gives a moderate stabilization energy of 48.98 kJ/mol.
The NBO, which is composed of its natural atomic hybrids (h A and h B ), displays sp hybridization, their polarization coefficients (C A and C B ), and atomic orbital (percentage s-character, p-character) (Supplementary Material, Table S5). The σ(C1-N2) is constituted by the sp 2.16 of carbon (C1) and sp 1.9 of nitrogen (N2). The C1 here has 31.61% s character and 68.28% p character, and N2 has 34.43% s character and 65.25% p character. As a result, the NBO C1-N2 is a bond between a C1 sp 2.16 hybrid and a N2 sp 1.9 hybrid. The larger polarization coefficient (C B = 0.7915) of the N2 atom compared to the C1 atom (C A = 0.6112) reflects its stronger electronegativity. The p-component of the hybrid described by the polar (θ) and azimuthal (Φ) angles is utilized to specify its orientation. The hybrid's directionality is determined by comparing the central line direction between the two nuclei to the maximum angular amplitude. The bond bends, which is then stated as the angle of divergence between these two directions. Because of the optimization of geometry, the orientation of geometry changes (Supplementary Material, Tables S6 and S7). The carbon NHOs of NBO σ C1-N2 (NBO 1) and σ C1-N6 (NBO 2) are bent away from the line of centers of C-N by 5.3 • and 5 • , respectively, while the nitrogen NHOs are bent away by 1.2 • and 1.4 • , respectively.

Population Analysis
Analysis of atomic charge is essential because it determines the atoms NMR chemical shifts and impacts the system's dipole moment, polarizability of molecules, electronic structure, and molecular reactivity. Using the B3LYP/6-311++G(d,p) method and basis set together with the Mulliken population analysis (MPA) method, the charges on the atoms of the 5-HMU were determined (Supplementary Material, Table S8 and Figure S3). There are five carbon atoms, three oxygen, and two nitrogen atoms as depicted in Figure S3. Almost all carbon atoms are positively charged except C9. Out of the carbon atoms, C1 attached with the hydroxy group has the highest positive charge due to the electron-withdrawing nature of oxygen. C9 was the most negatively charged of all the atoms. The oxygens are almost similar in charge, same with the case of nitrogen atoms. It is usual that all the hydrogen atoms have the same charges with a little difference due to the presence of different groups attached to them.
The Fukui function is a useful tool based on Mulliken population analysis. It gives quantitative information such as the electron densities and determines the most electrophilic and nucleophilic sites on a molecule. The electron density decreases upon addition of electron to some spots providing the negative value of the Fukui function. On the contrary, positive values mean an electron is removed and electron density is increased.
As per Parr and Yang [52], the Fukui function's value can be used to determine the reactivity. As a result, the condensed Fukui function's value increases and so does the atomic center's level of reactivity. When the molecule is attacked by an electrophile or when it loses electrons, the reactivity is denoted by the symbol f-, whereas when it is attacked by a nucleophile or when it gains electrons, the reactivity is denoted by f+. The value of Fukui functions along with the local softness of each atom obtained using Mulliken population analysis has been summarized in Supplementary Material, Table S8. It can be seen that C9 and O8 have the highest values for f + and f − , respectively. These sites are likely to undergo nucleophilic and electrophilic attacks, respectively. The reactivity order for electrophilic attack is O8 > C4 > O7 > N2 > O10 > N6 > C5 while the order for nucleophilic attack is C9 > C1 > N2 > C3 > O10 > N6, and the order for radical attack is C9 > C1 > C3. The electrophilic attack is shown to have higher reactivity than the other two.

Frontier Molecular Orbital (FMO) and UV-Vis Spectra
The UV-Vis spectrum of 5-HMU was measured experimentally in DMSO solvent. Using the TD-DFT/B3LYP/6-311++G(d,p) basis set and the PCM model polarizable continuum model, the UV-Vis spectrum was estimated in the DMSO and methanol solvents. The excitation energy, oscillatory strength, and bandgap are summarized in Supplementary Materials; Table S9. The calculated absorption maximum in gas, DMSO, and methanol were 261.67, 252.24, and 251.81 nm, respectively, while the experimental value in DMSO was observed at 275 nm. The calculated and experimental UV-Vis spectra of the titled molecule were given in Figure 7. It is very much clear from the spectra that the molecule shows a (red and hyperchromic) shift in DMSO and MeOH because of secondary interaction with the solvent.
There are various methods for calculating the excitation energy. The first is the difference between the lowest unoccupied molecular orbital (LUMO) and highest occupied molecular orbital (HOMO) of a neutral system (LUMO). These orbitals regulate the manner in which medications interact with other molecules, such as interactions between pharmaceuticals and their receptors. The HOMO electrons' propensity to transfer to the LUMO can be accurately determined qualitatively using the frontier molecular orbitals (FMO). The HOMO, LUMO, and FMO gaps, which are crucial aspects of quantum chemistry [53], serve as illustrations of the chemical reactivity and kinetic stability of molecules. While the LUMO orbital mostly functions as an electron acceptor, the HOMO orbital primarily serves as an electron donor. At B3LYP/6-311++G(d,p), the HOMO, LUMO, and orbital energy gaps for the 5-HMU were determined. The highest molecular orbital value that is occupied is −7.12 eV, and the lowest molecular orbital value that is empty is −1.68 eV (Supplementary Materials; Table S10). HOMO-LUMO, HOMO-1-LUMO+1, and HOMO-2-LUMO+2 have energy gaps of 5.4338 eV, 7.0066 eV, and 7.9288 eV, respectively (Figure 8). It can be seen from the plot of HOMO that the positive and negative regions are spread equally all over the molecule same as in the case of LUMO. The charge transfer interaction between the molecules is explained by the energy gap between HOMO and LUMO. The smaller the energy gap (5.44 eV), the easier the charge transfer between the two orbitals due to the electron acceptor's strong electron-accepting capacity. The bioactivity of the molecule is due to the high charge transfer contact. There are some important tools such as ionization potential (I.E), electron affinity (E.A), electronegativity (χ), chemical hardness (η), and chemical potential to check the order of stability of the molecular system (Supplementary Materials; Table S10). HOMO and LUMO energies are used to compute chemical hardness and chemical potential (µ). A negative chemical potential is defined as electronegativity. Chemical hardness could be a good predictor of chemical stability. Soft molecules are those with a little low energy gap, while hard molecules are those with a big energy gap [54,55].  Table  S10). HOMO and LUMO energies are used to compute chemical hardness and chemical potential (μ). A negative chemical potential is defined as electronegativity. Chemical hardness could be a good predictor of chemical stability. Soft molecules are those with a little low energy gap, while hard molecules are those with a big energy gap [54,55].  When the chemical hardness of a substance is raised, it is said to have higher stability. A hard molecule is said to have a wide energy gap when compared to soft molecules [56]. The calculated chemical hardness of 5-HMU is 2.71. The electrophilicity index of can be used to describe a compound's biological potential, and the calculated electrophilicity index of 5-HMU is 3.55. Furthermore, it has very low softness value of 0.369. It provides information on the toxicity of various contaminants based on their reactivity and site selectivity [57].  Table  S10). HOMO and LUMO energies are used to compute chemical hardness and chemical potential (μ). A negative chemical potential is defined as electronegativity. Chemical hardness could be a good predictor of chemical stability. Soft molecules are those with a little low energy gap, while hard molecules are those with a big energy gap [54,55].  When the chemical hardness of a substance is raised, it is said to have higher stability. A hard molecule is said to have a wide energy gap when compared to soft molecules [56]. The calculated chemical hardness of 5-HMU is 2.71. The electrophilicity index of can be used to describe a compound's biological potential, and the calculated electrophilicity index of 5-HMU is 3.55. Furthermore, it has very low softness value of 0.369. It provides information on the toxicity of various contaminants based on their reactivity and site selectivity [57]. When the chemical hardness of a substance is raised, it is said to have higher stability. A hard molecule is said to have a wide energy gap when compared to soft molecules [56]. The calculated chemical hardness of 5-HMU is 2.71. The electrophilicity index of can be used to describe a compound's biological potential, and the calculated electrophilicity index of 5-HMU is 3.55. Furthermore, it has very low softness value of 0.369. It provides information on the toxicity of various contaminants based on their reactivity and site selectivity [57].

Electron Excitation (Electron and Hole Density Distribution) Analysis
Electron and hole density distribution provides the concept of multi-molecular orbital excitation that delivered obvious features of the excited state for the receptors and the corresponding receptor anion. Once an electron is delocalized from A to B as a part of single electron excitation processes (where A and B are real space functions), the excitation of an electron from an occupied to virtual molecular orbital (MO) is taken as a model for excited state. The ground state of MO can be compared to the various regions of photoexcited states shown by the density distribution maps for electrons (ρ ele (r)) and holes (ρ hole (r)). According to T. Lu and F. Chen, the electron and hole density distribution (EDD and HDD) can be explained by MO wavefunction (Φ) and the configuration coefficient (w), which resembles to the transition of an electron from occupied MO(i) to a virtual MO (1) as shown in equations below [58]: i→l i→l j =i→l The electron and hole density distribution maps for 5-HMU were calculated using the aforementioned method. The B3LYP/6-311++G(d,p) level calculations were performed using PCM in DMSO and MeOH and the Multiwfn 3.3.4 software [58].
The strong absorption band maxima at 275 nm (4.51 eV) corresponding to HOMO-LUMO were predicted by calculations for the electronic structure of 5-HMU at the TD-DFT-B3LYP/6-311++G(d,p) level with PCM in DMSO. This transition is shown in Figure 9, whereas the absorption at max = 251.81 nm (4.92 eV) for MeOH. The estimated EDD map in Figure 9A depicts a significantly thicker surface on the ring and C=O groups, while the HDD map depicts a denser isosurface on the whole molecule in DMSO. Figure 9B shows EDD and HDD maps created in MeOH; the EDD map exhibits a thicker surface on the entire ring and C=O groups, and the HDD map displays a denser surface on the entire molecule (Supplementary Materials; Table S11).
Molecules 2023, 28, x FOR PEER REVIEW 14 of 22 Electron and hole density distribution provides the concept of multi-molecular orbital excitation that delivered obvious features of the excited state for the receptors and the corresponding receptor anion. Once an electron is delocalized from A to B as a part of single electron excitation processes (where A and B are real space functions), the excitation of an electron from an occupied to virtual molecular orbital (MO) is taken as a model for excited state. The ground state of MO can be compared to the various regions of photoexcited states shown by the density distribution maps for electrons (ρ ele (r)) and holes (ρ hole (r)). According to T. Lu and F. Chen, the electron and hole density distribution (EDD and HDD) can be explained by MO wavefunction (Φ) and the configuration coefficient (w), which resembles to the transition of an electron from occupied MO(i) to a virtual MO (1) as shown in equations below [58]: The electron and hole density distribution maps for 5-HMU were calculated using the aforementioned method. The B3LYP/6-311++G(d,p) level calculations were performed using PCM in DMSO and MeOH and the Multiwfn 3.3.4 software [58].
The strong absorption band maxima at 275 nm (4.51 eV) corresponding to HO-MO-LUMO were predicted by calculations for the electronic structure of 5-HMU at the TD-DFT-B3LYP/6-311++G(d,p) level with PCM in DMSO. This transition is shown in Figure 9, whereas the absorption at max = 251.81 nm (4.92 eV) for MeOH. The estimated EDD map in Figure 9(A) depicts a significantly thicker surface on the ring and C=O groups, while the HDD map depicts a denser isosurface on the whole molecule in DMSO. Figure 9(B) shows EDD and HDD maps created in MeOH; the EDD map exhibits a thicker surface on the entire ring and C=O groups, and the HDD map displays a denser surface on the entire molecule (Supplementary Materials; Table S11).

Thermodynamic Properties
The B3LYP/6-311++G(d,p) technique and basis set by ORCA software were used to calculate the thermodynamic characteristics of 5-HMU. The values of the enthalpy (H 0 ), entropy, and Gibbs free energy were computed between 100 K and 700 K. The intensity of molecular oscillations also rises with temperature because of higher translational and rotational energy, which is supported by uniform distribution [59,60]. In light of the fact that thermal stability is crucial both at low and high temperatures, physicochemical effects are observed in the industry to preserve food or medicine storage. As a result, temperature is important in chemical thermodynamics and pharmacodynamics. From the Gibbs free energy, biological thermodynamic applications such as protein stability, membrane, DNA, and enzyme kinetics can be calculated.
The Gibbs free energy (G 0 ) shows an inverse relation with temperature, and hence it decreases with an increase in temperature (Supplementary Materials; Figure S4 and Table S12). The correlation between the temperature and the thermodynamic function was determined using quadratic and linear equations. The correlation (R 2 ) between the two for entropy, enthalpy, and Gibbs free energy are 1, 0.9999, and 0.9747, respectively. The relationships are plotted in Supplementary Material; Figure  The thermodynamic properties also influence how non-covalent interactions such as electrostatic interactions, multiple hydrogen bonds, van der Waals interactions, and hydrophobic effects bind to one another [61]. Utilizing the relationship between thermodynamic functions and the second rule of thermodynamics, these features can be used in the thermochemical field to determine the flow of chemical reactions and identify other thermodynamic energy [62].

Molecular Docking
The ligand-protein docking method is used to predict the ligand's binding sites in proteins using a three-dimensional network. Predicting chemical interactions between the lead compound and biologically relevant targets (proteins, enzymes, etc.) has become critical in the science of drug discovery. Furthermore, by analyzing the docking of the lead molecule with a variety of protein targets, one can gain insight into the underlying molecular mechanisms of selectivity [63].
The two fundamental steps of the docking procedure are the prediction of the ligand structure as well as its orientation within these sites and the evaluation of the binding affinity. The Swiss ADME online source predicted the target protein that needs to dock with the 5-HMU and downloaded the file in RCSB protein database (pdb) format. The relevant target protein ID (docked with a similar ligand as hydantoin) was picked from the Swiss ADME-target prediction site and retrieved from the protein data bank with the aid of a smile string of 5-HMU ligand, which was added to the site of the Swiss target prediction (PDB). The active sites for docking are ligands that have already docked inside proteins. The proteins 4P6T 4P6S, 6QXD, 6EI4, 5OAE, and 5M67 all belonging to the oxidoreductase class except 5M67 and were docked with 5-HMU showing the bond distance of hydrogen bond interaction between 5-HMU and the target ( Figure 10A-F). In this docking, all the metal centers of the proteins are avoided in PDB's. The binding energies have been calculated using AutoDock vina [64,65]. The ligand's characteristics, including its binding energy and inhibition constant, are collected in Table 5. The lowest binding energy (∆G) and intermolecular energy were found to be −5.0 kcal/mol between 5-HMU and 6QXD in the protein-ligand interaction. The highest binding energy was found for 5M67. The calculated bond distance for the same object is 2.270. The solid lines in the figure depict the hydrogen bonds between the ligand and the targeted protein. Inhibition constants were calculated by K i = exp(∆G/RT). The two fundamental steps of the docking procedure are the prediction of the ligand structure as well as its orientation within these sites and the evaluation of the binding affinity. The Swiss ADME online source predicted the target protein that needs to dock with the 5-HMU and downloaded the file in RCSB protein database (pdb) format. The relevant target protein ID (docked with a similar ligand as hydantoin) was picked from the Swiss ADME-target prediction site and retrieved from the protein data bank with the aid of a smile string of 5-HMU ligand, which was added to the site of the Swiss target prediction (PDB). The active sites for docking are ligands that have already docked inside proteins. The proteins 4P6T 4P6S, 6QXD, 6EI4, 5OAE, and 5M67 all belonging to the oxidoreductase class except 5M67 and were docked with 5-HMU showing the bond distance of hydrogen bond interaction between 5-HMU and the target ( Figure 10A-F). In this docking, all the metal centers of the proteins are avoided in PDB's. The binding energies have been calculated using AutoDock vina [64,65]. The ligand's characteristics, including its binding energy and inhibition constant, are collected in Table 5. The lowest binding energy (∆G) and intermolecular energy were found to be -5.0 kcal/mol between 5-HMU and 6QXD in the protein-ligand interaction. The highest binding energy was found for 5M67. The calculated bond distance for the same object is 2.270. The solid lines in the figure depict the hydrogen bonds between the ligand and the targeted protein. Inhibition constants were calculated by Ki =exp(∆G/RT).

Drug-Likeness
To decide the structural appearances of ligands in order to achieve effective and efficient results, a set of rules known as "drug-likeness" is employed following the Lipinski's rule of five [66,67]. We looked at the polar surface area (TPSA), molar refractivity, and the quantity of hydrogen bond donors and acceptors (HBD and HBA, respectively) (Supplementary Materials; Table S13), which shows that these metrics have values within acceptable ranges for the medication candidate under evaluation. The HBD and HBA should be less than 5. The title molecule satisfies these requirements. The MLogP, one of the most critical characteristics (showing the hydrophobic/lipophilic nature of the molecule), is equal to −1.22 which is far below Lipinski's rule of five threshold value of 5. The TPSA value for a compound is 85.95 Å 2 , which is much below the maximum cut-off value of 140 Å 2 (Supplementary Materials; Figure S6). The molar refractivity was calculated to be 33.81, which falls within the range of pharmaceutical candidates [68].

Molecular Dynamic (MD) Simulation
The efficacy of ligand binding on protein structure was analyzed by using MD modeling. The GROMACS 2022 software was utilized to study MD simulation [69]. The 'CHARMM36 force field' was used to generate the protein topology. Each system was solvated using the TIP3P water model before being neutralized with the proper concentrations of Na + and Cl − ions [70]. The docked protein-ligand combination of the molecule with the PDB ID: 6QXD was used to evaluate the stability of the 5-HMU. To confirm the stability of the complex structures, the simulation was run for 20 nanoseconds at room temperature. The system was designed to use as little energy as possible, employing the steepest descent technique. Position restraint together with NVT and NPT ensembles were used after the energy minimization process. Following temperature stabilization, an isothermal isobaric ensemble (NPT) was performed. The trajectory indicates that the docked ligand was stable throughout the simulation in the protein's binding site where the protein-ligand is relatively stable and exhibits minor conformational changes. The data were examined using the graphical tool origin once all of the MD simulation phases had been completed. The built-in functions of GROMACS 2022 were used to find the different vital parameters, such as root mean square deviation (RMSD), number of hydrogen bonds occurring between protein and the ligand, and interaction energy. The stability of the complex is represented by the RMSD stabilization ( Figure 11A). After alignment, further RMSD values were determined for the protein with an inbuilt option in the VMD software; therefore, the protein showed a smaller variation in the RMSD values. Thirteen stable hydrogen bonds are formed between the protein and the ligand, and these bonds define how strong the complex is ( Figure 11B). However, in this instance, the total interaction energy is useful. According to the common formula for adding two values, this value is equal to −137.34 ± 11.8 kJ mol −1 after propagating the error. According to every outcome of the MD simulation, 5-HMU can bind to the active sites of proteins to create a stable protein complex. The resultant values of MD simulations are more precise and correct in comparison to docking results. We may predict that there is more fluctuation in the RMSD values of ligand arisen due to movement all 13 hydrogen bonds of ligand during the simulation, but not at once. Hydrogen bonds occupancy with their donor-acceptor sides showed in Supplementary Materials; Table S14.

Materials and Instrumentations
The spectral study was conducted using the 5-hydroxymethyl uracil (Sigma-Aldrich Delhi, India). The infrared spectrum was measured with FT-IR PerkinElmer (Serial No. 105627) using a KBr disc method (FT-RAMAN at IIT Kanpur, Kanpur, India). NMR spectra were achieved with the help of a Bruker Advance 400/Avll HD-300 (FT-NMR) apparatus. The UV-Vis analysis was carried out on a UV-1280 Spectrophotometer (Shimadzu, Mumbai, India).

Computational Analysis
The supporting information has provided computational details (See the Supplementary Materials).

Materials and Instrumentations
The spectral study was conducted using the 5-hydroxymethyl uracil (Sigma-Aldrich Delhi, India). The infrared spectrum was measured with FT-IR PerkinElmer (Serial No. 105627) using a KBr disc method (FT-RAMAN at IIT Kanpur, Kanpur, India). NMR spectra were achieved with the help of a Bruker Advance 400/Avll HD-300 (FT-NMR) apparatus. The UV-Vis analysis was carried out on a UV-1280 Spectrophotometer (Shimadzu, Mumbai, India).

Computational Analysis
The supporting information has provided computational details (See the Supplementary Materials).

Conclusions
The title chemical 5-HMU was examined through the spectroscopic analysis using FTIR, FT-Raman, NMR ( 1 H and 13 C), and UV-Visible in accordance with the theoretical DFT method. The optimized crystal structure of 5-HMU and its trimer that was produced with the least amount of energy were compared to the experimental crystal structure. The development of an intermolecular H-bond with adjacent monomers was found to stabilize the trimer, and its ideal molecular structure was established using the 6-311++G(d,p) basis set and the DFT-B3LYP approach. The complete vibrational assignment was investigated based on potential energy distribution (PED). The nucleophilic and electrophilic sites, as well as the H-bonding interactions of the 5-HMU, were predicted via MEP mapping. The positive potential sites were found on hydrogen and carbon atoms, whereas the negative potential sites were found on oxygen atoms. The electron localization function (ELF) was employed to represent the degree of electron localization. The resonance interaction energy of 5-HMU was 48.98 kJ/mol, as measured by the maximum stabilization of electron leading LP(1) N6-π* C5-O7. The deviation angle, coupled with the polar (θ) and azimuthal (φ) angles of the vector from the nucleus, were utilized to determine the orientation of a hybrid. The energy gap between highest occupied and lowest unoccupied was determined to be 5.4338 eV. Three metrics, such as dipole moment, polarizability, and hyperpolarizability, were used to validate a material nonlinear optical characteristic. High chemical hardness (2.72) indicated the compound stability, while low chemical softness (0.37) was designated as non-toxicity, and the electrophilicity index (3.56) indicated 5-HMU biological activity. EDD and HDD maps for the first excited-state were created with DMSO and MeOH to indicate single MO pair excitation. In the temperature range of 100 K to 700 K, the ORCA program was utilized to determine various thermodynamic features. It depicts how the thermodynamic properties change as the temperature increases. The Fukui function describes the intramolecular interaction between atoms in a molecule. The 5-HMU receptor has the lowest binding energy of −5.0 kcal/mol, indicating that further research into the therapeutic potential of the mentioned molecule is recommended. The molecular dynamic simulation gave better idea of binding and showed 11 stable hydrogen bonds.

Conflicts of Interest:
There are no conflict to declare.