Synthesis and Computational and X-ray Structure of 2, 3, 5-Triphenyl Tetrazolium, 5-Ethyl-5-phenylbarbituric Acid Salt

: The title compound triphenyl tetrazolium (TPT) of phenobarbital, 5-Ethyl-5-phenylbarbituric acid triphenyl tetrazolium salt (4) was prepared by the reaction of 5-Ethyl-5-phenyl-2,4,6( 1H, 3H, 5H )-pyrimidinetrione, monosodium salt (1) with triphenyl tetrazolium chloride (3) in deionized water at an ambient temperature through a cation exchange reaction. Colorless crystals of compound four suitable for an X-ray structural analysis were obtained by slow evaporation from acetonitrile. Compound four had crystallized in the monoclinic space group, P 21/c , with a = 15.3678 (9) Å, b = 12.2710 (7) Å, c = 21.8514 (13) Å, β = 109.867 (2) ◦ , V = 3875.5 (4) Å 3 , and Z = 4. A Through density functional theory (DFT) calculations, the probable molecular association structure in the phenobarbitone -triphenyl tetrazolium solution was studied. With the 6-311G-(d,p) basis set, the gas phase features of the phenobarbital-triphenyl tetrazolium clusters with a phenobarbitone dimer and water molecules, including an optimum structure and intermolecular hydrogen bonding, were investigated in detail. In addition, the positions and strengths of the intermolecular hydrogen bond interactions between the phenobarbitone and triphenyl tetrazolium molecules were analyzed using atoms in molecule (AIM) analysis, reduced density gradient (RDG) methods, the XRD method, and the non-covalent interaction (NCI) index method. In addition, the molecular electrostatic potential (MEP) surfaces were analyzed to determine the electrophilic and nucleophilic centers.


Introduction
Phenobarbital sodium (5-Ethyl-5-phenylbarbituric acid), or sodium salt, has the chemical structure C 12 H 11 N 2 NaO 3 , and its relative molecular mass 254. 22. Sodium phenobarbital is available as 30, 60, 65, and 130 mg/mL injections and as a sterile powder in 120 mg ampules [1]. Phenobarbital, also known as phenobarbitone, is a medication recommended by the world Health Organization in developing countries for the treatment of certain types of epilepsy [2]. Phenobarbital is widely used in the treatment of partial and generalized tonic-colonic seizures in all age groups [3] and is considered a first-line drug for treating seizures and status epilepticus in newborns.
On the other hand, 2, 3, 5-triphenyltetrazolium chloride is a heterocyclic compound with a five-member ring which contains four nitrogen atoms; one of these atoms bears a positive charge. The use of 2,3,5-tetraphenyl tetrazolium salt for the extraction and spectrophotometric and potentiometric determination of various elements and ions has been reported [4]. Recently, triphenyl tetrazolium salt was used as an ion-pair reagent for PVC membrane sensors for many target analytes [5][6][7].
In this study, we hope to report herein the synthesis and X-ray structure of the tetraphenyl tetrazolium salt of phenobarbital. Noncovalent interactions dominate the chemical interactions between a protein and a drug and a catalyst and its substrate, as well as within the self-assembly of nanomaterials [8,9] and even some chemical processes [10][11][12]. This class of interactions includes hydrogen bonding, dipole-dipole interactions, steric repulsion, and London dispersion throughout a broad range of binding energies [13]. The molecular structure is determined by covalent, noncovalent, and electrostatic interactions, the latter two of which drive a majority of biological reactions. Covalent bonds are easily identifiable in a three-dimensional molecule structure, while noncovalent interactions are obscured by the absence of bonds. Although there are numerous methods for viewing and analyzing covalent and electrostatic interactions, there is no corresponding approach for noncovalent interactions. This technique would be useful in many fields, including the study of self-assembled materials and the development of new medicines [14].
In this study, we provide a method for mapping and analyzing noncovalent interactions that requires only information on molecular geometry, complementing existing approaches for covalent and electrostatic interactions. On the basis of the aforementioned factors, one of the objectives of this work is to provide a fundamental understanding of the main interaction between phenobarbitone and triphenyl tetrazolium that may be used to predict the reactivity of and provide an understanding for this interaction. The optimal structures of phenobarbitone and triphenyl tetrazolium were determined using DFT calculations. Hydrogen bonding and intermolecular interactions were analyzed using the atoms in molecules (AIM), reduced density gradient (RDG), and non-covalent interaction (NCI) techniques.

General
The melting point (uncorrected.) was determined using a Gallenkamp melting point apparatus. X-ray crystallography was measured on a Bruker APEX-II D8 venture diffractometer equipped with graphite monochromatic Mo Kα radiation, with λ = 0.71073 Å at 100 (2) K. The IR spectra were recorded with a Perkin-Elmer FTIR spectrometer. Bruker 500 and 700 MHz and 125 and 176 MHz instruments were used to record 1H NMR and 13 C NMR in DMSOd6, respectively, using TMS as an internal standard (with chemical shifts in δ ppm). The mass spectrum was measured on an Agilent Triple Quadrupole 6410 QQQ LC/MS equipped with an ESI (electrospray ionization) source.

X-ray Crystallography General
Single crystals of compound four were obtained by slow evaporation from acetonitrile. A suitable crystal was selected for X-ray diffraction analysis. The data were collected by a Bruker APEX-II D8 Venture diffractometer equipped with graphite monochromatic Mo Kα radiation, with λ = 0.71073 Å at 100 (2) K. Cell refinement and data reduction were completed by a Bruker SAINT, and the program used to solve the structure and refine the structure was SHELXS-97 [15]. The final refinement of the collected data was performed by full-matrix least-squares techniques with anisotropic thermal data for the non-hydrogen atoms on F 2 . All the hydrogen atoms were placed in calculated positions and constrained to ride on their parent atoms. Absorption correction by a multi-scan method was performed using SADABS software. The crystal data and refinement data are shown in Table 1. 1.00 and −0.59

Computational Study
The Gaussian 09W [16] software package was used to complete all quantum-chemical computations. To optimize the geometry of the ion pairs and calculate their energies, the B3LYP functional [17,18] method with the GD3 version of Grimme's dispersion correction [19] was used in conjunction with the 6-311G-(d,p) [20] basis set. The geometries of the phenobarbitone, triphenyl tetrazolium, and phenobarbitone dimer, and those of the water and triphenyl tetrazolium complexes, were employed as starting points for constructing the basic structures of the complex of the formations ion pairs. Numerous configurations of the ion pairs were produced as a result of these interactions. One of the generated geometries was consistent with and was employed in the same-level vibrational frequency computations to describe all stationary locations as minima (no imaginary frequencies) and to evaluate their thermodynamic properties [21].
The proton affinities (PA) of the triphenyl tetrazolium and acid conjugated bases were calculated as the differences between the enthalpy values of the cations and acids and their corresponding phenobarbitone and anions, though with the sign reversed. The differences between the free energy of the ion pair and the sum of the free energies of the triphenyl tetrazolium, phenobarbitone, phenobarbitone ion and water molecules were used to compute the change in the Gibbs free energy G 298 associated with ion pair formation (at the standard condition). The interaction energy E int between the ions in the ion pair was calculated using the super molecule approach [19,22], that is, it was calculated as the difference in energy between the ion pair and the ions that comprised it. The Boys and Bernardi counterpoise methods [23] were used to the optimized structures of the ion pairs in order to determine the basis set superposition error (BSSE). When the BSSE was considered, the resultant interaction energy was slightly lower (by no more than 2%), but the relative order of the energies in the analyzed compounds' series remained the same.
To find the hydrogen bonding interaction between the ion pairs, the following criteria were used: (i) the geometric properties of the H···O fragment [24], (ii) the predicted energy of hydrogen bonds using Espinosa's equation [25,26], (iii) the topological parameters of the bond critical point at the H···O contacts [27,28], and (iv) the quantity of charges transferred from the electron donor to the electron acceptor [29]. Bader's quantum theory of atoms in molecules (QTAIM) technique was used to analyze the topological electron density [30,31]. This approach has also been applied in number of studies aimed at elucidating noncovalent interactions in task-specific ILs [21,32]. Within the QTAIM framework, the interactions between atoms are intimately related to the topological properties of the electron density ρ(r), specifically, the set and types of critical points at which its gradient is zero. The bond critical point BCP (3, −1) and the bond path going through it are of particular importance in our research since they are required for the chemical bond or, in the general case, for the stabilization of the interatomic interaction between the two bonded atoms. As a result, the higher the ρ(r) value at the BCP, the greater the concentration of electronic charge in the surface at this point and the stronger the considered contact. In comparison to covalent bonding (where ρ(r) is~10 −1 au), the value of ρ(r) for hydrogen bonding [27] and the van der Waals interactions are fairly tiny, with 102 au for hydrogen bonding [27] and 10 −3 au for the van der Waals interactions [33]. The (r) value of the interacting atoms at their BCP and the value of their Laplacian ∇2ρ(r), the total energy density H(r) [27], and the ratio of the absolute potential energy density to the kinetic energy density |V(r)|/G(r) [28] are the most frequently utilized parameters used to describe the nature and strength of bonding interactions.
GaussView 06 [37] was used to show the molecular structure and molecular electrostatic potential (MEP) surface. Mulliken atomic charge and local reactivity descriptors were derived using Mulliken population analysis by computing the single point energies of the N, (N − 1), and (N + 1) species of the molecule using the 6-311G (d, p) basis set. All DFT calculations were performed at the 5-AU molecule's ground state energy level, with no constraints on the potential energy surface.

X-ray Crystallography
The title compound (5-Ethyl-5-phenylbarbituric acid triphenyl tetrazolium salt (4)) was crystallized as the salt of one anion-cation in the presence of one neutral phenobarbital molecule and one water solvent molecule in an asymmetric unit ( Figure 1). Comparison with other crystal structures of 2,3,5-triphenyltetrazolium salts revealed that the internal dimensions of the [TPT]+ cation are largely insensitive to the local environment [39]. The five atoms in the tetrazolium ring (N5-N6-N7-N8-C25) were coplanar, and N5-N6 and N7-N8 were almost equivalent (1.314 and 1.317 Å) ( Table 2). The N6-N7 bond (1.328 Å) was longer than the former two bonds. This may be ascribed to the repulsion between the relevant phenyl groups. The torsion angles of the phenyl rings were 116.63, 117.72, and 171.35 • . In the phenobarbital molecules, the C-N and C=O bond lengths in the two barbituric rings did not vary significantly between the structures ( Table 2).
This crystal structure had two types of hydrogen bonds that stabilized the structure: the first type involved ten intermolecular H bonds (between separate molecules), and the second type involved one intramolecular H bond (between parts of the same molecule). These hydrogen bonds were detected in a three-dimensional framework structure, forming a chain extending along the a-axis (Table 3) (the Cambridge Crystallographic Data Center (CCDC) CCDC:1436012).
The cation TPT and the anion PBT were linked by two hydrogen bonds (Figure 2), where the carbonyl groups in the anion PBT interacted as donor groups with the hydrogen atoms C27 and C34 of the cation as acceptors. While the cation (TPT) and the neutral PBT were connected by two hydrogen bonds, the carbonyl group of the neutral PBT interacted as a donor, with the hydrogen atoms of C28 and C29 as acceptors, forming three hydrogen bonds. On the other hand, the anion PBT and the neutral PBT interacted through four hydrogen bonds. The first two hydrogen bonds existed between the carbonyl groups (O5 and O6) in the anion PBT as donors, with the hydrogen atoms (N2 and C20) in the neutral PBT as acceptors. The other hydrogen bonds were between the carbonyl groups (O1) and the amine (N4) in the neutral PBT, with the hydrogen atoms (N2 and C18) in the anion PBT as acceptors.     (6) 167 (6) N2-H1N2···O6 iii 0.91 (7) 1.94 (7) 2.827 (5) 163 (5) N2-H1N2···N4 iii 0.91 (7) 2.60 (6) 3.361 (5)  In the anion PBT molecule, a hydrogen bond was formed by the intramolecular interaction between N4 and O6. In addition, the water molecules interacted with the anion in O6 and N4 as acceptors. Finally, the bond in the anion molecule was an intermolecular bond between H3 and O5.

Molecular Geometry
The geometries depicted in Figure 3 were selected to best reflect the energetically favored ion pair combinations and their thermodynamic stability and to maximize the hydrogen bond interaction, as well as to agree with the experimental crystal geometries. Additionally, all phenobarbitones were optimized for both protonated and deprotonated structures. Figure 3 shows the optimized geometries of the molecules. All these geometries were the local minima on the potential energy surface. The mean absolute errors for the (21) bond distances and the (36) bond angles are listed in Table 4. The MAE values ranged from 0.0213 to 0.0003 for the complex bond distances. The complex bond angles had MAE values of below 2.455 • . B3LYP was utilized to obtain these MAE values.

Interaction Energies (IE)
To obtain the most stable configurations and binding energies of structures, the DFT approach at the B3LYP level using the 6-311G(d,p) basis set has been widely used. Phenobarbitone, water, triphenyl tetrazolium, and the phenobarbitone dimer, water, and triphenyl tetrazolium were investigated as a complex in this research. The complexes (molar ratio of 1:1) in the gas phase have complexation energy (raw), BSSE energy, complexation energy (corrected), and Gibbs free energy changes (∆G) of −104.59 kcal/mole, 0.0234, −89.88, and −11.814 kcal/mol, respectively. The phenobarbitone dimer, water, and triphenyl tetrazolium complex has the lowest complexation energy, indicating that it is the most stable phenobarbitone dimer, water, and triphenyl tetrazolium complex yet discovered. It also has a negative ∆G value, indicating that it is simple to generate spontaneously, and this result agreed with the experimental method.
To obtain the most stable configurations and binding energies of the structures, the DFT approach at the B3LYP level using the 6-311G(d,p) basis set was used. Before investigating the complexation energy of the different crystal systems, the optimization of the complex structures at different ratios of TPT-PBT -(PBT-1) was studied, resulting in three equilibrium structures without imaginary frequencies (frequencies of less than zero): TPT:PBT (1:1), TPT:  (Table 5). Consequently, these complexation energy results for the complex of the phenobarbitone dimer and water molecules are lower than those for the complexes with a phenobarbitone dimer without water molecules and a single phenobarbitone, indicating that it is more stable than the other complexes. Moreover, the complexation energies of these complexes have negative values, indicating that they can be easily formed spontaneously and that the results are consistent with the experimental method. The results of the complexation energy tests are consistent with the results of the XRD crystallography tests on the real crystal.

Non-Covalent Interaction (NCI) Index
As with the AIM method, the reduced density gradient (RDG) is another effective technique for accounting for non-covalent interactions. Non-covalent interactions between molecules can be visualized using RDG scatter plots and non-covalent interaction (NCI) plots [11,37]. In this method, RDG is plotted against electron density and multiplied by the sign of the second eigenvalue (sign(λ 2 )ρ) [11], and both the inter-and intramolecular weak interactions may be as shown in Figure 4. On the negative scale (blue color), the RDG dispersed spots represent the H bonding interactions, while the spikes (green color) and positive scale of (sign(λ 2 )ρ) reflect the van der Waals interactions and the steric repulsions, respectively.
As seen in Figure 4A, the H bond of the studied complex can be shown at (0.025 a.u.). Therefore, the RDG area, which is shown by black circles, corresponds to the complexes O84· · · H98, O96· · · H86, and O57· · · H9.
In addition, as seen in Figure 4B, the NCI plot shows the H bonds, van der Waals interaction, and steric effect of the studied complex. The blue circles represent non-covalent bonds (H bonds), which suggests three H bonds. It is worth noting that the ions in the complex comprise aromatic moieties that appear to be involved in the van der Waals interaction with the strong RDG spikes, which are shown with a green color between the ions in Figure 4B. The steric effect, which is responsible for repulsion between ions, is represented by the red circle in Figure 4B.

Quantum Theory of Atoms in Molecules (QTAIM)
The quantum theory of atoms in molecules, developed by the late Prof. R. F. W. Bader, has been extensively employed in several studies to determine the topological features of various forms of π-hole interactions [29,30]. AIM analysis was completed on all compounds using the DFT:B3LYP/6-311G(d,p) basis set to gain a better understanding of the noncovalent H· · · O bonds.
The classification of the bonds between molecules in the complex at the bond critical points (BCP) 121, 138, 156, 163, 174, 177, 187, 191, 197, 204, 207, 209, and 225 was determined in accordance with Koch and Popelier [33,34], based on the findings of the topological analysis of the electron density for the complex depicted in Table 6 and Figure 5. Because their Laplacian ∇ 2 ρ(r) values are positive and in the range of (0.0112-0.1025 a.u.), and their total energy density H(r) values are positive and in the range of (0.0006-0.002 a.u.), as well as, the values of the ratios of absolute potential energy density to kinetic energy density |V(r)|/G(r) in the interval (0.7356-0.957) are less than 1. Therefore, all bonds are closed-shell interactions (van der Waals interactions or weak H bonds). Moreover, the bonds at the BCP were separated into two parts based on the value of the ρ(r) and ∇ 2 ρ(r): the first part, at BCP 163, 207, and 225, was characterized as van der Waals interactions since the (ρ(r)) was less than (0.002 a.u.) and the ∇ 2 ρ(r) was also less than (0.024 a.u.), which is within the limits of van der Waals interactions. The second part, at the bond critical points (BCP) 121, 138,156,174,177,187,191,197,204, and 209, were classified as hydrogen bond interactions because the (ρ(r)) value and the Laplacian ∇ 2 ρ(r) at those bond critical points were in the ranges of (0.0056-0.0285 a.u.) and (0.0214-0.1025 a.u.), respectively, which is typical for H bond interactions. Table 6. AIM parameters of the chosen H bonds at the bond critical points (BCPs) for the interaction between the phenobarbitone and triphenyl tetrazolium complexes.

MESP Analysis
The molecular electrostatic potential (MEP) surfaces [20,31] are the best graphical representations of the electrostatic potential across the surface of a molecule for the purpose of identifying the electrophilic and nucleophilic centers. The MEP surfaces are colorcoded, with blue and red representing the most positive and negative regions, respectively, and with green representing the neutral region. As seen in Figure 6, the MEP maps of phenobarbitone, water, triphenyl tetrazolium, and the complex of phenobarbitone dimer, water, and triphenyl tetrazolium complex are depicted in the gas phase. The MEP plot of the acceptor (triphenyl tetrazolium) is characterized by a positive zone (blue) in the center (surface map value of 0.152 au), which is regarded as an electrophile. The negative area originates from the C=O (−0.04088 and −0.207 au,) groups of phenobarbitone, respectively. With regards to the phenobarbitone primary negative area (red), which is located on the O15, O18, and O19 atoms (−0.0406, −0.0337, and −0.0336 au, respectively), it can be called an n-donor (nucleophile). Following the creation of a complex between the donor and acceptor, the C=O values are increased and the N2 atom of triphenyl tetrazolium is decreased to obtain a lower value. These results imply an n-electron transfer from the donor's N2 to the phenobarbitone's C=O groups. The PCM model's surface map values likewise produced close results. As a result, the ESP map surfaces exhibited excellent agreement with the experimental results.

Reactivity Descriptors
Numerous reactivity descriptors, including ionization potential (Ip), electron affinity (A), chemical potential (µ), hardness (η), electrophilicity index (ω), and softness (σ), were estimated from the HOMO(N), HOMO (N + 1), and HOMO (N − 1) surfaces, providing insight into the reactivity of the chemical reactions. Equations are used to describe these characteristics. Table 7 summarizes the electrical interactions of triphenyl tetrazolium and phenobarbitone with the generated CT complex characteristics. The electrical properties of the phenobarbitone and triphenyl tetrazolium molecules are deduced from this table. When determining a molecule's HOMO−LUMO energies, a high EHOMO indicates a good electron donor, whereas a low ELUMO indicates a good electron acceptor. Because phenobarbitone has a lower ELUMO than triphenyl tetrazolium in gas and PCM analysis, it is regarded as an electron acceptor; yet, because triphenyl tetrazolium has a higher EHOMO than phenobarbitone, it is considered an electron donor. Additionally, the chemical potential is an index of potential that specifies the direction of electron flow between molecules. Electrons flow from a structure with the highest chemical potential to one with the lowest chemical potential. With regards to this interpretation, triphenyl tetrazolium

Reactivity Descriptors
Numerous reactivity descriptors, including ionization potential (Ip), electron affinity (A), chemical potential (µ), hardness (η), electrophilicity index (ω), and softness (σ), were estimated from the HOMO(N), HOMO (N + 1), and HOMO (N − 1) surfaces, providing insight into the reactivity of the chemical reactions. Equations are used to describe these characteristics. Table 7 summarizes the electrical interactions of triphenyl tetrazolium and phenobarbitone with the generated CT complex characteristics. The electrical properties of the phenobarbitone and triphenyl tetrazolium molecules are deduced from this table. When determining a molecule's HOMO−LUMO energies, a high E HOMO indicates a good electron donor, whereas a low E LUMO indicates a good electron acceptor. Because phenobarbitone has a lower E LUMO than triphenyl tetrazolium in gas and PCM analysis, it is regarded as an electron acceptor; yet, because triphenyl tetrazolium has a higher E HOMO than phenobarbitone, it is considered an electron donor. Additionally, the chemical potential is an index of potential that specifies the direction of electron flow between molecules. Electrons flow from a structure with the highest chemical potential to one with the lowest chemical potential. With regards to this interpretation, triphenyl tetrazolium possesses a stronger chemical potential than phenobarbitone. Additionally, the electrophilicity of phenobarbitone is greater than that of triphenyl tetrazolium, indicating that phenobarbitone is the better electrophile and should be considered an e-acceptor, while triphenyl tetrazolium is an e-donor. Additionally, the softness values and these results established that triphenyl tetrazolium is an electron donor in gas and PCM analyses, whereas phenobarbitone is an electron acceptor. Table 7. Calculated HOMO(N), HOMO (N + 1), and HOMO (N − 1) energy bands, as well as the chemical potential (µ), electronegativity (χ), global hardness (η), global softness (S), and global electrophilicity indexes (ω) for tetrahydrofuran (IE in eV) and its derivatives at the B3LYP/6-311G(d,p) level.  Figure 6 shows that HOMO is localized only on phenobarbitone, while LUMO is localized on the triphenyl tetrazolium moiety. The electron affinity (EA) and ionization potential (IP) can be calculated using Koopman's theorem [10]: EA = −ELUMO and IP = −EHOMO. Table 7 provides an overview of some of the parameters calculated using DFT. The EHOMO of the complex has a value of 7.37 eV, which is comparable to the EHOMO of the donor phenobarbitone, which has a value of 9.08 eV. On the other hand, the ELUMO of the complex has a value of 1.34 eV, which is closely related to the ELUMO of the acceptor triphenyl tetrazolium (−0.05 eV). The HOMOs and LUMOs follow the same trend in other donor acceptor systems [4]. The HOMO-LUMO plots and the energy gap show the charge transfer between phenobarbitone and triphenyl tetrazolium within the complex.

Conclusions
In conclusion, the title compound four was prepared efficiently by the reaction of phenobarbital sodium (1) with triphenyl tetrazolium chloride (3) in deionized water at an ambient temperature. The structure of compound four was established on the basis of its X-ray single crystal analysis. In the present study, a TPT, a PBT, and their complex were described, and their geometric structures and electrical characteristics were fully examined using experimental-and DFT-level theoretical calculations. This demonstrated that the estimated geometric characteristics, such as bond length and angle, were in excellent agreement with the XRD crystallography results. The HOMO and LUMO tests were used to determine the drug's energy gap, chemical activity, and charge transfer. MEP testing has also been shown to identify electrophilicity and nucleophilicity zones. NCI analysis identifies van der Walls interactions and weak interactions. All these analyses provided results which concluded that the complex, with a ratio of 1:2 (TPT: PBT), was more stable, and this suggestion agreed with the XRD results. In addition, the complex interacted through charge transfer interaction.