Synthesis, DFT Calculations, In Silico Studies, and Antimicrobial Evaluation of Benzimidazole-Thiadiazole Derivatives

In this study, a series of new benzimidazole-thiadiazole hybrids were synthesized, and the synthesized compounds were screened for their antimicrobial activities against eight species of pathogenic bacteria and three fungal species. Azithromycin, voriconazole, and fluconazole were used as reference drugs in the mtt assay. Among them, compounds 5f and 5h showed potent antifungal activity against C. albicans with a MIC of 3.90 μg/mL. Further, the results of the antimicrobial assay for compounds 5a, 5b, 5f, and 5h proved to be potent against E. faecalis (ATCC 2942) on the basis of an acceptable MIC value of 3.90 μg/mL. The cytotoxic effects of compounds that are effective as a result of their antimicrobial activity on healthy mouse fibroblast cells (L929) were evaluated. According to HOMO–LUMO analysis, compound 5h (with the lower ΔE = 3.417 eV) is chemically more reactive than the other molecules, which is compatible with the highest antibacterial and antifungal activity results. A molecular docking study was performed to understand their binding modes within the sterol 14-α demethylase active site and to interpret their promising fungal inhibitory activities. Molecular dynamics (MD) simulations of the most potent compounds 5f and 5h were found to be quite stable in the active site of the 14-α demethylase (5TZ1) protein.


INTRODUCTION
Antimicrobial resistance (AMR) is a significant public health concern of the twenty-first century.It threatens the effective prevention and treatment of an increasing variety of microbial diseases that are resistant to the traditional drugs used to treat them. 1 Microbial resistance is one of the most burning problems in clinical practice, and one of the main objectives of current biomedical research is to discover new, powerful drugs that can combat multiresistant bacteria.−4 Benzimidazole is known as a lucky structure in pharmaceutical chemistry and has a variety of biological functions.Benzimidazole has a benzene ring system in which the benzene ring is attached to a five-member imidazole ring having nitrogen atoms at positions 1 and 3; thus it is known as a heterocyclic aromatic compound. 5,6Additionally, benzimidazole is a tremendous scaffold of therapeutic importance with promising pharmacological properties.The safety and efficacy profiles of benzimidazole medications in clinical use are well established.Because of their therapeutic effects, benzimidazole derivatives have attracted a lot of interest in the medical field, providing excellent results as anticancer, 7 antiviral, 8 antimicrobial, anti-inflammatory, antihistamine, 9 antihypertensive, anti-tubercular, analgesics, antiulcer, and anthelmintics. 10The benzimidazole ring is present in many significant medications that are therapeutically employed in the research field.Albendazole (anthelminthic), Bendamustine (anticancer), Omeprazole (antiulser), and Astemizole (antihistamine) are examples of drugs with a benzimidazole structure. 11,12hiadiazoles are five-membered heterocyclic compounds with two nitrogen atoms and one sulfur atom.They are a type of azole molecule. 13During the last decades, 1,3,4-thiadiazole derivatives have drawn much attention due to their biological and pharmaceutical activities and have been investigated increasingly due to their numerous therapeutic and industrial applications, which is due to the presence of =N−C−Smoiety. 14A variety of 1,3,4-thiadiazole is in use, like Acetazolamide (diuretic), Cefazolin, Cefazedone (antibiotics), Megazol (antiprotozoal), Timolol maleate (NSAIDs), Methazolamide (carbonic anhydrase inhibitor), and Sulphamethizole (antibacterial). 15n this study, a new series of benzimidazole-thiadiazole derivatives were synthesized and characterized by 1 H NMR, 13 C NMR, and HRMS.Synthesized compounds were screened for their antimicrobial activities against eight species of pathogenic bacteria [Escherichia coli (ATCC 25 922), Serratia marcescens (ATCC 8100), Klebsiella pneumoniae (ATCC 13 883), Pseudomonas aeruginosa (ATCC 27 853), Enterococcus faecalis (ATCC 2942), Bacillus subtilis, Staphylococcus.aureus (ATCC 29 213), S. epidermidis (ATCC 12 228)] and three fungal species Candida albicans (ATCC 24 433), C. krusei (ATCC 6258), C. parapsilosis (ATCC 22 019)].The cytotoxic effects of the final compounds that are effective as a result of antimicrobial activity on healthy mouse fibroblast cells (L929) were evaluated.With the use of Candida's 14α-demethylase (CYP51), molecular docking investigations were also carried out.Using 100 ns molecular dynamics (MD) simulations, we investigated the stability of compounds containing CYP51 was investigated.It is crucial to know the precise structure of the ligand with the lowest minimum energy when conducting a molecular docking investigation.Density functional theory (DFT) was used to model the eight newly synthesized chemicals.

RESULTS AND DISCUSSION
2.1.Chemistry.The target molecules were synthesized in five processes, as shown in Figure 1.To obtain the sodium metabisulfite addition product of the aldehyde, the methyl 4formylbenzoate compound aldehyde was first treated with sodium metabisulfite in ethanol.In the second step, methyl 4-(1H-benzo[d]imidazole-2-yl)benzoate (2) was produced as a consequence of the condensation reaction between the sodium metabisulfite product and benzen-1,2-diamine under reflux.Compound 2 was treated with hydrazine hydrate in ethanol in the following step to produce compound (3).Ethanol was refluxed with the hydrazide derivative compound and the corresponding isothiocyanate derivatives, and the precipitated product was filtered out.The thiosemicarbazide molecule was cyclized in the presence of strong sulfuric acid to produce the thiadiazole derivatives (5a−i) in the last step.1 and 2 provide a summary of the antimicrobial activity test results and the minimum inhibitory concentrations (MICs) of the compounds (5a−i).
Among the series, compounds 5a, 5b, 5f, and 5h were found to be the most active molecules with an MIC value of 3.90 μg/ mL, and they are specific toward the Gram-positive bacterial, E. faecalis (ATCC 2942).Compounds 5a exhibited noteworthy activity with an MIC value of 7.81 μg/mL against S. aureus (ATCC 29213) and S. epidermidis.Compound 5d When the antifungal activity test results of the compounds were examined, it was found that the compounds were generally more effective against C. albicans.According to the antifungal activity test, compounds 5f and 5h exhibited better potency with an MIC value 3.90 μg/mL against C. albicans comparable to the reference drug fluconazole (MIC 7.81 μg/ mL) and voriconazole (MIC 3.90 μg/mL).Compounds 5b, 5d, 5e, and 5g were observed to have significant antifungal activity on C. albicans.The MIC values for these compounds were determined to be 7.81 μg/mL.
Structure−activity relationship (SAR), after closely reviewing the data for antibacterial activity, the following conclusions were reached: • Activity data revealed that compared to Gram-negative bacteria, the majority of the synthesized benzimidazolethiadiazole derivatives shown superior potency against Gram-positive bacteria.• In most of the cases, ethyl substituted showed better antibacterial inhibitory activity than 2-chloroethyl substituted against E. faecalis (ATCC 2942).• N-isopropyl substituted showed better antibacterial inhibitory activity than N-propyl and N-isobutyl substituted.
• It was determined that the synthesized compounds were more sensitive to C. albicans.2.3.Cytotoxicity Assay.In order to evaluate the cytotoxic effects of compounds (5a, 5b, 5d, 5e, 5f, 5h) on healthy cells, the compounds that are effective as a result of antimicrobial activity were selected and their cytotoxic effects on healthy mouse fibroblast cells (L929) were evaluated.The cell line appearances for all compounds are shown in Figure 2. The IC 50 values of the compounds are given in Table 3.
As a result of cytotoxicity studies conducted on healthy cell cultures, the IC 50 values of all compounds except compound ″5h″ were found to be higher than 100 μM.Considering the MIC values of the compounds in antimicrobial activity, it is seen that they are much lower than the IC 50 values at which the compounds have a toxic effect on healthy cells.From this information, it is concluded that the compounds show antimicrobial activity and are not toxic at the MIC values.

In Silico Studies. 2.4.1. Quantum Chemical
Calculations.The density functional theory (DFT) calculations were performed by using the Gaussian 09 program 16 with the B3LYP exchange correlation functional with the 6-311G(d,p) basis set.GaussView 5.0 program 17 was used to generate the input geometries and visualize the results.The optimized geometries of all structures correspond to true minima as no imaginary frequencies are observed in the vibration frequency investigation.The molecular electrostatic potential (MEP) and HOMO−LUMO analyses were performed at the B3LYP/6-311G(d,p) level in order to   investigate the electronic properties of current molecules.MEP and HOMO−LUMO diagrams are shown in Figure 3.
The molecular geometry parameters regarding the optimization results of all molecules are in good agreement with the experimental studies possessing similar main building blocks in the literature. 18The main molecule in the building blocks of compound 5h consists of a 1, 3, 4-thiadiazole ring with an amine group attached to a methoxy-phenyl ring and a phenyl-benzimidazole group at the second and fifth locations of this ring, respectively.N−C single and N =C double bonds of benzimidazole ring are 1.379 and 1.316 Å, respectively, which are in line with the experimental values of 1.381 (3) Å and 1.324 (3) Å in literature. 18The rotation of the phenylbenzimidazole group with respect to the 1, 3, 4-thiadiazole ring is measured as −179.6°,while the rotation of amine group attached to methoxy-phenyl ring to this ring is 180.0°.These values indicate that the molecule is planar.All molecules showed the same trend.
According to the HOMO−LUMO analysis, 5h is chemically more reactive than the other molecules (Table 4).The chemical reactivity order is 5h> 5d > 5i> 5e > 5g> 5a > 5f> 5b > 5c.5c, having the higher energy gap value (ΔE = 3.756 eV), is more stable and less reactive than the others, which is in line with the lowest antibacterial and antifungal activity results.
As known from the literature, a lower value of ionization potential (IP) indicates that it has a better property of electron donor. 19According to Table 4, it is seen that the ionization potential of 5h is the lowest than the others, which refers to a better property of electron donor.A good description of donor−acceptor properties sheds light on hydrogen bonding interactions.This context provides a preliminary assessment for ligand-protein interactions.The electrophilic indexes (ω) of all molecules belong to the strong electrophiles group as their values are bigger than 1.50 eV. 20UMOs are distributed throughout the whole molecule throughout the conjugated π system, while HOMOs are localized in electronegative nitrogen atoms and aromatic rings, as shown in Figure 3.In the 5h molecule, since the methoxy group resonated with the phenyl ring, the electron density of the ring increased, and HOMO is localized in this part, while LUMO distribution is located in aromatic rings except the methoxy phenyl ring.
According to the MEP diagram (Figure 3), in all molecules, negative regions possessing the high electron density are seen around the N atoms in the benzimidazole and thiadiazole ring, which are responsible for electrophilic attacks.In addition to these negative regions, 5h has also a negative zone around the oxygen of methoxy.Positive regions with the low electron density of all molecules are formed around both the N−H group in the benzimidazole ring and N−H group bonded to the thiadiazole ring, which are responsible for nucleophilic attacks.
2.4.2.Molecular Docking Studies.Molecular docking studies of synthesized thiadiazole derivatives (5a−i) were performed using Glide (Grid-based Ligand Docking with Energetics) to understand their binding modes within the sterol 14-α demethylase active site and to interpret their promising fungal inhibitory activities.This target was chosen since it has been noted that mostly antifungals inhibit the enzyme lanosterol 14-α − demethylase, which is dependent on cytochrome P450 and depletes Ergosterol, a crucial component of the fungal cell membrane.Table 5 21,22 It was observed that all the synthesized compounds nicely docked into the active site of C. albicans sterol 14-α demethylase with a good docking score ranging from −9.15 to −8.209 kcal/mol. 23,24According to the docking result analysis, the most active molecules 5f (MIC= 3.90 μg/mL) and 5h (MIC= 3.90 μg/mL) were effective against C. albicans and had significant binding affinity toward the 14-α demethylase protein with docking scores of −8.25 and −9.077 kcal/mol, respectively.The standard drugs Voriconazole, Fluconazole, and the cocrystallized inhibitor VT1161 had docking scores of −7.004, −6.633, and −8.528 kcal/mol, respectively.In previous studies and X-ray crystallographic structures, the binding pocket of VT1161 was identified with residues, such as Tyr64, Tyr118, Leu121, Thr122, Phe126, Ile131, Tyr132, Phe228, Pro230, Phe233, Gly303, Ile304, Gly307, Gly308, Thr311, Leu376, His377, Ser378, Phe380, Tyr505, Ser507, and Met508.Our docking results of synthesized thiadiazole derivatives with the 14-α demethylase protein also showed a similar docking profile.The 2D and 3D visual representations of the representative compound 5f and 5h interactions depicted in Figure 4 were generated using Maestro's ligand interaction tool.Both of these compounds lodge in the active site with a similar approach to cocrystallized ligands and show hydrophobic and hydrogen bonding with different crucial residues.The binding interaction of compound 5f shows that it forms hydrogen bonds with Tyr132 and ionic interactions with Hem 601 through the benzimidazole scaffold.At a distance of 3.63 Å, Tyr118 tacking with the central hydrophobic phenyl ring of compound 5f (Figure 4A).In the case of compound 5h, three hydrophobic interactions are visible with the benzimidazole (His377) and thiadiazole (Tyr1188) scaffolds and one hydrogen bond with Tyr132 in the active site of the 14-α demethylase protein (Figure 4B).The tetrazole-based antifungal drug candidate VT1161 demonstrates significant interactions with key residues within the 14-α demethylase protein, notably, Tyr118, Tyr132, and His377.Specifically, representative compounds 5f and 5h establish contact with these critical residues, underscoring their importance in the antifungal mechanism.The comparable docking scores of these compounds further suggest robust binding affinity, reinforcing the hypothesis that they exert antifungal effects by effectively suppressing the activity of the 14-α demethylase protein.This interaction profile highlights the potential of compounds 5f and 5h as promising candidates for further exploration in the development of antifungal therapeutics targeting C. albicans.

Result Obtained from Molecular Dynamics
Simulation.MD simulation for compounds 5f and 5h in complex with 14-α demethylase (5TZ1) protein selected for MD simulations along a time scale of 100 ns to validate the results of molecular docking and assess their conformational stability.The RMSD of α carbon atoms in all systems is analyzed to understand their stability.The RMSD of the 5f-5TZ1 complex reaches to ∼1.73 Å from 0 to 50 ns, and after that, the system maintains the average RMSD value of 1.77 Å until the end of the simulation (Figure 5A).While assessing the RMSD of the 5h-5TZ1 complex, a steady decrease in MSD is observed after 15 ns, and after 51 ns until 100 ns, a slight fluctuation is observed, indicating that the trajectories generated during these times are stable.The average RMSD for 5f is 1.71 Å, whereas the RMSD for 5h is detected around 1.60 Å, indicating higher structural stability in both complexes.The average ligand RMSD values for 5f and the 5h are 2.54 and 2.76 Å, respectively.For the 5f-5TZ1 complex, the ligand RMSD shows an initial lower value until 20 ns, after which small fluctuations occur.The RMSD remains consistent overall, indicating minor adaptations or adjustments in the ligand conformation during the simulation.On the other hand, in the case of the 5h-5TZ1 complex, the ligand RMSD exhibits no major fluctuations throughout the simulation, except for the initial time fluctuations attributed to equilibration.This suggests that the conformation of the 5h ligand remains stable and undergoes fewer structural changes compared to those of the 5f ligand during the entire simulation period (Figure 5B).
RMSF values are also calculated to determine the dynamic behavior of both protein residues.Maximum fluctuations of 3.20 and 3.04 Å are detected in residues Leu239 and Asp269, respectively.The visual analysis of MD simulation trajectories suggests that all the drugs engaged in significant binding interactions with the hotspot residues, namely, Gly307, Ile131, Leu121, Leu376, Leu87, Met508, Phe126, Phe233, Phe380, Pro230, Ser378, Tyr118, and Tyr132 of the 14-α demethylase protein.It was found that all of these interacting residues had  RMSF values of 0.411 and 1.628 Å, indicating the stability of compounds 5f and 5h during simulation (Figure 5C).Less fluctuation is indicative of a stronger protein structural stability.The structural compactness of both complexes is analyzed by using the radius of gyration (RGyr).The radius of gyration is the estimated distance that a body's mass is concentrated around an axis.It is frequently used to learn more about the conformational and folding behaviors of macromolecules. 23,24It is an evaluation of a macromolecule's overall compactness and 3-D structure under diverse settings.The wide range of RGyr values indicated that the simulation may have caused the protein to unfold.RGyr is found in the range of 21.004−21.232and 20.820−21.116Å for 5f−5TZ1 and 5h−5TZ1 complexes, respectively, which indicates that the 14α demethylase protein is almost stable during the whole run (Figure 5D).In a protein−ligand complex, the strength of a ligand molecule's binding to the target protein is determined by the presence of hydrogen bonds.
Compounds 5f (orange) and 5h (green) show significant intermolecular hydrogen bonding with the 14-α demethylase protein in the MD simulation at 100 ns, with a range from 0 to 3. The quantity and distribution of hydrogen bonds are displayed in Figure 5D.The maximum and mean hydrogen bond numbers of the 3a: CDK5A1 complex, 5f−5TZ1 complex, and 5h−5TZ1 complex were 3, 3, and 1.1 and 0.7, respectively (Figure 5E).Hydrogen bond analysis not only confirms the docking results but also shows that no conformational changes occurred in binding positions during the MD simulation.

CONCLUSION
In summary, a class of novel benzimidazole-1,3,4-thiadiazole hybrids was designed and synthesized and evaluated for antimicrobial activity.Activity data disclosed that most of the synthesized benzimidazole-thiadiazole derivatives showed better potency against Gram-positive bacteria as compared to Gram-negative bacteria.When the antifungal activity test results of the compounds were examined, it was found that the compounds were generally more effective against C. albicans.Compounds 5a, 5b, 5f, and 5h were found to be the most active molecules with MIC value 3.90 μg/mL against E. faecalis (ATCC 2942).According to the antifungal activity test, compounds 5f and 5h exhibited better potency with an MIC value of 3.90 μg/mL against C. albicans.According to the HOMO−LUMO analysis, compound 5h (with the lower ΔE = 3.417 eV) is chemically more reactive than the other molecules, which is compatible with the highest antibacterial and antifungal activity results.The effects of the compounds on the L929 mouse fibroblast (normal) cell line were studied to determine the site of cytotoxicity.Furthermore, molecular docking studies are used to predict how designed or synthesized compounds interact with the target protein/ enzyme.Molecular dynamics (MD) simulation of the most potent compounds 5f and 5h were found to be quite stable in the active site of the 14-α demethylase (5TZ1) protein.

MATERIAL AND METHOD
The complete chemicals used in the synthetic process were acquired from Merck Chemicals (Merck KGaA, Darmstadt, Germany) or Sigma-Aldrich Chemicals (Sigma-Aldrich Corp., St. Louis, MO, USA).The compounds' uncorrected melting points were ascertained using an MP90 digital melting point instrument (Mettler Toledo, OH, USA).A Bruker digital FT-NMR spectrometer (Bruker Bioscience, Billerica, MA, USA) operating at 300 and 75 MHz registered the 1 H-and 13 C NMR spectra of the compounds produced in DMSO-d 6 , respectively.In the NMR spectra, splitting patterns were denoted as follows: s for singlet, d for doublet, t for triplet, and m for multiplet.The reported unit of coupling constants (J) was Hertz.The Shimadzu LC/MSMS system (Shimadzu, Tokyo, Japan) was used to determine the M + 1 peaks.Using Silica Gel 60 F254 TLC for thin-layer chromatography (TLC), all reactions were observed.
4.1.Chemistry.4.1.1.Synthesis of Sodium Metabisulfite Salt of Benzaldehyde (1) Derivative.Ethanol was used to dissolve methyl 4-formyl benzoate (5g, 0.03 mol).Drop by drop, sodium metabisulfite (6.84 g, 0.036 mol) in ethanol was added to the benzaldehyde solution.The reaction mixture was agitated for 1 h at room temperature following the completion of the dripping process.The precipitated product was removed by using a filter.(2).After dissolving benzen-1,2-diamine (0.022 mol) in DMF, sodium metabisulfite salt of benzaldehyde derivative (7.09 g, 0.026 mol) was added.Pouring the reaction contents into freezing water at the end of the reaction allowed the result to precipitate.After being filtered out, the precipitated product crystallized from the ethanol.3).Compound 2 (0.018 mol) was added to the same vial, along with an excess of hydrazine hydrate (5 mL) and 15 mL of ethanol.For 12 h, the mixture was refluxed.Following the completion of the reaction, the mixture was placed in freezing water, and the result was filtered.standard methods CLSI M07-A9 (2012) and NCCLS M27-A2 (2002), as described in the previous study. 25,26.3.Cytotoxicity Assay.The effect of the compounds between 5a−i on the viability of the L929 cell line was analyzed by MTT assay.The MTT method was performed as previously described. 27.4.In Silico Studies.4.4.1.Quantum Chemical Calculations.Using the gradient-corrected correlation functional of Lee, Yang, and Parr (LYP) and Becke's threeparameter exchange functional (B3), electronic characterization of the produced compounds was performed in the gas phase.The optimization was validated using frequency calculations, and the structures' minimum energy was found when imaginary frequencies were absent.28 The same technique was used to perform MESP analysis using these optimized structure HOMO−LUMO energies (electronic characteristics, such as the ionization potential, electronegativity, electrophilic index, nucleophilic index, and chemical potential produced from these energies).Quantum computations were carried out using Jaguar software, and the molecular orbitals were examined using the maestro interface.

Molecular Docking and Molecular Dynamic (MD) Simulation
Studies.The Schrodinger Glide was used to conduct molecular docking studies for the newly thiadiazole derivatives (5a−i) on the crystal structure of sterol 14-α demethylase (CYP51) from C. albicans cocrystallized with the tetrazole-based antifungal drug candidate VT1161(PDB ID: 5TZ1) receptor to evaluate their insilico inhibitory effects.The protein preparation wizard was used to protonate the protein under physiological pH and to apply the OPLS3e force field to these two sensors. 29,30As previously mentioned, each protein was rectified and 3D hydrogenated, and energy was minimized.Chem.Draw was used to create the 2D structures of the pyrimidine derivatives (5a−i).Then, LigPrep was used to generate 3D structures, add charges, minimize energy, and compile all structures into a single molecular database file. 31inally, docking investigations were carried out utilizing the Glide docking wizard and Standard Procedure (SP) as the docking procedure, with docking results visualized using the Maestro Graphical User Interface (GUI).To determine the atomic-level binding stability of the top-ranked compounds and gain insights into their molecular interactions, molecular dynamics (MD) simulations were conducted using the Desmond module of Schrodinger.The 5f−5TZ1 complex and 5h−5TZ1 complexes were solvated under orthorhombic periodic boundary conditions, maintaining a 10 Å buffer region between protein atoms and box edges with the explicit SPC water model.In the system builder, sodium and chloride ions were introduced to neutralize charges, and a 0.15 M NaCl salt concentration was added to mimic human physiological conditions. 32The built system was then minimized using the fixed parameters of the OPLS3e force field to eliminate electronic clashes and appropriately align the protein structure within the simulation boundaries. 33Long-range electrostatic interactions were evaluated using the smooth particle mesh Ewald approach with a tolerance of 1e−09, while short-range van der Waals and Coulomb interactions were computed with a cutoff radius of 9.0 Å.After importing the minimized build system (.cms file) into the molecular dynamics module, a 100 ns simulation was conducted under an'isothermal−isobaric ensemble' (NPT) at a temperature of 300 K and a pressure of 1 bar.The'Nose-Hoover chain thermostat' and'Martyna-Tobias-Klein barostat' techniques were employed at 100 and 200 ps intervals for isothermal−isobaric conditions, respectively. 34Simulation snapshots were retrieved at 100 ps intervals, and the resulting trajectories were analyzed.All computational modeling is performed on a workstation, featuring Ubuntu 22.04.2LTS 64-bit configuration, Intel Xeon W-2245 @ 3.90 GHz, 8 cores, CUDA 12, and NVIDIA RTX A4000 graphics processing unit.

Figure 3 .
Figure 3. Molecular Electrostatic Potential (MEP) and HOMO− LUMO diagrams of compounds 5a−i at the B3LYP/6-311G(d,p) level.Atom colors: carbon in gray, nitrogen in blue, chlorine in green, oxygen in red, sulfur in yellow, and hydrogen in white.The surfaces plotted by the 0.0004 electrons/b3 contour of the electronic density.For 5a molecule: Color ranges: in Au: blue, more positive than 0.0630; green, between 0.0630 and 0; yellow, between 0 and − 0.0630; red, more negative than − 0.0630.

Figure 5 .
Figure 5. MD simulation trajectory analysis of ligand−protein complexes.A. Time-dependent RMSD plot; B. ligand RMSD is C. individual amino acids RMSF plot; D. time-dependent radius of gyration plot; E. time-dependent hydrogen bond analysis.

Table 4 .
HOMO−LUMO Energies (eV) and Calculated Global Reactivity Parameters of the Best Stable States of Compounds 5a−i at the B3LYP/6-311G(d,p) Level in the Gas Phase a