In the present study, a new complex compound of zinc with LA (ZnLA) was synthesized, the structure of diaqua [5-(1,2-dithiolan-3-yl) pentanoato] acetate zinc [Zn(C8H13O2S2)(CH3COO)(H2O)2)] was determined by XRD technique and the spectroscopic properties were investigated by FT-IR. In addition, the thermal behavior of the complex was investigated. Theoretical calculations have been carried out by using Density Functional Theory (DFT) method in the ground state. The frontier orbitals (HOMO-LUMO) and Mulliken population method were obtained computationally. The def2-TZVPP was chosen as basis set for the theoretical calculations. Also docking investigation calculation of ZnLA calculated with using CB-Dock 2 server and Auto Dock Vina program. Comparison of the experimental and theoretical results shows that there is a good agreement between them.
2.1. IR-spectroscopic Analysis Results
Solid-state IR-Fure spectra of LA and ZnLA powders were obtained by using an IR Affinity-1S (Shimadzu, Japan) spectrometer (
Figure 1). The specific vibrational frequencies of the complex compound were determined by comparison with those of the ligand (
Table 1).
It is notable that new bands appeared at 3266-3000, 1357 and 464 cm
−1 and in the spectrum of ZnLA, which were not observed in LA (
Figure 1,
Table 1). Absorption band frequency between 3266-3000 cm
−1 is due to the presence of O-H stretching mode of the hydroxyl group [
33]. This is explained by the occurrence of the hydroxyl group of the water molecule in the complex compound. Concurrently, the absorption band at 464 cm
−1 corresponds to the
stretching vibration. In general, the Zn-O bond absorption band corresponds to the area of 400-600 cm
−1. In some studies, the Zn-O bond absorption band was reported to be at 426 [
34] and 470 cm
−1 [
35], which is similar results with our work. Also, the absorption band in the 1357 cm
−1 region observed in the IR spectrum of ZnLA corresponds to the valence vibration of the
CH3 group in the acetic acid residue. The bands appearing at 2925 and 2852 cm
−1 are characteristic to CH
2 stretching mode [
34,
35]. Also, in the spectrum of ZnLA, absorption bands characteristic of associated acids disappeared, and intense absorption bands characteristic of
and
appeared in the 1535 and 1448 cm
−1 regions, respectively. The expression ∆
=
-
is located in the interval 87 (1535-1448) cm
−1 in ZnLA. From this, we can conclude that metal ions in ZnLA are coordinated through oxygen atoms in the ligand (LA).
2.3. DTA-TGA analysis results
The lipoic acid crystallized melts at 64°C, according to the first endotherm in differential thermal analysis (DTA) analysis (
Figure 3A). The thermal gravimetric analysis (TGA) curve shows that the material is stable up to 150°C, experiencing weight loss of 64% at ∼160–320°C. This loss is associated with the second endotherm in the DTA curve. These correspond to the total decomposition of LA [
36]. The differential thermal analysis of ZnLA revealed the decomposition of water molecules bound by hydrogen bonds at 115.84°C, and the endothermic effect observed at 158.37°C revealed the decomposition of water coordinated with the central atom (
Figure 3B). Decomposition of the LA residue from the compound was observed at 442.96°C. This temperature is close to the melting point of zinc and is used in its identification. The endothermic effect observed at this moment was characterized by anion decomposition (mass loss of 45.124%) and exothermic peaks that arose due to the formation of CO
2 and SO
2 as a result of combustion and oxidation at a temperature of 458°C. At this moment in 160-300°C decomposition of acetate and from 300 to 840°C it was decomposition of lipoic acid residue. On the derivatogram, the decrease in mass continued up to a temperature of 839.75°C. Starting from a temperature of 840°C, no change in mass was observed. As a final product, ZnO and part of the carbon that remained unoxidized turned into coal (27.815%) and remained a solid residue [
33].
Based on the results of the physicochemical analysis, an approximate structural formula of ZnLA was proposed (
Figure 4A,B). In the given picture one can see that the coordination number of zinc is 6, and this complex compound consists of molecule LA, molecule of acetic acid residues, also 2 molecules of water, which are associated with zinc by a coordination bond.
2.5. Frontier Molecular Orbital (FMO) Analysis
It is known, both molecular orbitals of highest occupied (HOMO) and lowest unoccupied (LUMO), which are called frontier MO’s are significant quantum chemical factors to take part in the chemical reactivity and kinetic stability of a molecule. Usually, HOMO is known as donor, while LUMO is known acceptor [
39,
40]. On the basis of this MO’s can be determined significantly helpful quantum chemical parameters, such as chemical hardness (η) and softness (σ), energy bandgap (ΔE), electronegativity (𝜒) and electrophilicity index (ω). According to Koopmans' theorem, while HOMO energy is related to ionization potential, LUMO energy is related to electron affinity [
41]. Another significant parameter is energy band gap value (difference between HOMO and LUMO), which is also crucial for the stability of the molecules [
42]. The patterns of the HOMO and LUMO of the LA and ZnLA has been computed by a method of B3LYP [def2-TZVPP], which are displayed in
Figure 5A–C, and their related energy values are given in
Table 4. The positive and negative phases are symbolized with red and blue, respectively.
According to diagrams of FMO are presented in
Figure 5, which shows the electronic distribution in HOMO-2, HOMO-1 and HOMO with LUMO, LUMO+1 and LUMO+2 for ZnLA. The FMO phenomenon suggests that the intramolecular charge transfer (ICT) character of the HOMO to LUMO transition exist. The delocalization of HOMO on all the atoms of a molecule except carbonyl and OMe groups while concentration of LUMO at all the atoms except OMe and carbonyl group for ZnLA can be easily observed from
Figure 5. It is well known that HOMO orbitals represent bonding character, while LUMO reflects anti-bonding character. In addition to that, HOMO-LUMO gap (∆E) along with calculated energy levels associated to ZnLA are summarized in
Table 4.
As
Table 4 indicates that the ∆E value in respect of ZnLA is 4.637 eV. In general, compounds hold higher energy gap reveal lower chemical reactivity while compounds possess lower energy gap reveal higher chemical reactivity. EHOMO, ELUMO and energy gaps (∆E) associated to them are effectively employed to predict global reactivity descriptors [
43,
44,
45] which assist in a clear description of the chemical reactivity, internal charge transfer, and stability of ZnLA. The equations employed for calculation of global reactivity descriptors are given in
Table 5.
An inverse relation is found between energy gaps, softness as well as reactivity of species. Moreover, it has direct relation to the hardness and stability of the molecule. Consequently, the molecules with larger energy gaps are considered harder molecules which resist any change in electronic configuration and clearly links low reactivity with more kinetic stability. On the other side, the molecules with low energy gaps are soft molecules having tunable character, high reactivity and lower kinetic stability [
46]. The
ZnLA showed higher ionization potential values as compared to that of electron affinity values as presented in
Table 5. The observed values of chemical hardness for
ZnLA are higher than their softness. The greater chemical hardness, thus, indicates less reactivity or more stability of
ZnLA. The stability and reactivity of compounds also connected with chemical potential of the compound. There is a direct relation between the chemical potential of the compound and stability whereas it is inversely related to reactivity. The result with respect to stability and reactivity of
ZnLA on the basis of chemical potential is almost similar to that of the result obtained from chemical hardness data. However,
ZnLA quite stable overall as reflected by its chemical hardness.
2.6. Molecular Docking Analysis (Studies)
Molecular docking has been proved as an efficient tool to study and predict the binding mode, the binding affinity of a ligand with the proteins, which is an essential part of drug discovery. This method for drug design is a highly reliable program, cost-effective and time-saving [
47,
48]. In this work, to explore the biological activity of ZnLA and LA, molecular docking simulations have been performed using CB-Dock 2 server [
49] and AutoDock Vina software [
50], and the docked structures were visualized by using Biovia Discover Studio Visualizer software [
51]. The 3D geometries of ZnLA and LA have been built using Avogadro program package [
52]. The target protein structures of 2GPA and 1ZSX were taken from the PDB database [
53] as a
pdb file.
Docking analysis have been investigated by of CB-Dock 2 server and AutoDock Vina results.
In the first case, AutoDock Vina was used for docking investigation. At the beginning protein structures were cleaned from additional molecules and polar hydrogen atoms were added. Proteins were minimized applying Kollman’s all partial atomic charges. Once minimized, the protein is loaded in MGLTools [
54] creating a PDBQT file that contains a protein structure with hydrogens in all polar residues, and it is then used by AutoDock Vina software to obtain the affinity binding values for a ZnLA and LA. The docking site for ZnLA and LA structures on protein targets was defined by establishing a cube with the dimensions 40 × 40 × 40 Å, covering the binding site predicted for CB-Dock 2 with a grid spacing of 0.375 Å centered on the center of mass of the ZnLA and LA. Ten runs with AutoDock Vina were performed in all cases per each ZnLA and LA structures, and for each run the best pose was saved. The average affinity for best poses was taken as the final affinity value for ZnLA and LA.
In the second case, the structures of the cleaned proteins from additional molecules using Biovia DS visualizer program have been used for molecular docking studies by help of the CB-Dock2 server. CB-Dock automatically identifies the binding sites. By default, the number of docking cavities is set to 5. During analysis of docking results, the main attention was paid to the active center, in which tested compounds are localized. The molecular docking results have been visualized by Biovia DS visualizer software.
Molecular docking studies conducted by the authors revealed that α-lipoic acid binds well to 2GPA and 1ZSX and acts as a weak inhibitor [
55]. Proteins 2GPA and 1ZSX are proteins responsible for blood glucose levels. By inhibiting them, it is possible to reduce the amount of glucose in the blood. In order to predict the mechanism of action and initial biological activity of the substance synthesized by us, studies were conducted with 2GPA and 1ZSX proteins.
Glycogen phosphorylase - 2GPA, which is the main enzyme regulating glycogen metabolism, catalyzes the process of phosphorolysis, in which glycogen is broken down into glucose-1-phosphate. In muscle, glycogen is converted to glucose-1-phosphate, undergoes glycolysis in the liver to produce metabolic energy and is converted to glucose. By inhibiting the 2GPA enzyme, it is possible to prevent the phosphorylase process, which, in turn, leads to blocking the process of glucose formation from glycogen [
56].
The most suitable posture of ZnLa-2GPA, obtained through various interactions at the binding site, was found to have the highest binding energy of -7.95
0.35 kcal/mol (
Table 6) which is 1.35 times greater than the binding energy of LA to 2GPA. This indicates a favorable binding affinity. It should be noted that Arg49, Thr94, Gly135, Leu136, Asn187, Lys574, Asn284, Arg292, Phe285, His341, His377, Tyr573, Lys574, Glu672, Ser674, Gly675 and Gly677 were involved in the interaction with the ZnLA molecule via hydrogen bonds. Arg569 and Ala673 are formed hydrophobic bonds with ZnLa, Tyr185 formed pi-sulfur bonds with ZnLA, Tyr648 pi-sulfur, Glu126, Lys574, Ala383, and Glu672 were involved in the interaction through electrostatic attraction (
Figure 6A,B).
The amino acids Ser, Tyr and Phe in 2GPA play an important role in phosphorylation processes. In our study, 14 amino acids in 2GPA were connected with ZnLA by different bonds (
Table 6). These include Ser and Tyr. Thus, ZnLA may induce antidiabetic activity by inhibiting its allosteric moiety and binding more tightly to 2GPA. In silico studies have been carried out comparing ZnLA with LA as a ligand, which also inhibits 2GPA, but its binding energy is 5.9 kcal/mol. It can be seen that the complex of zinc ZnLA is connected relatively stronger (BE=--8.3 kcal/mol) with 2GPA than LA. It should be noted that this compound may also be a potential zinc transporter. Of course, this must be confirmed in in vitro and in vivo experiments.
The most suitable posture of ZnLa-1ZSX, obtained through various interactions at the binding site, was found to have the highest binding energy of -7.7
0.4 kcal/mol (
Table 6) which is 1.16 times greater than the binding energy of LA with 1ZSX. It is noteworthy that Trp43, Tyr76, Lys104, Arg175, Ser230 and Arg250 are involved in the interaction with the ZnLA molecule via hydrogen bonds. An electrostatic interaction of Asp71, Lys104, Arg175, Trp43, and Trp229 with the ZnLa molecule was observed. In particular, Trp43 interacts with zinc in the pi-cation and acetate residues in the pi-anion interaction. It can be seen that Trp43 and Trp43, Trp229, Trp258 with sulfur, carbon, and alkyl groups in the molecule create pi-sulfur, pi-sigma, and pi-alkyl hydrophobic effects (
Figure 7A,B).
1ZSX - Kv potassium channel beta subunit (KCNAB2) regulates cellular processes such as hormone secretion and repolarization of excitable cells. Prolongation of action potentials by blocking delayed rectifying potassium channels in pancreatic beta cells is known to increase intracellular free calcium and stimulate glucose-dependent insulin release [
57]. ZnLA may be a potential antidiabetic drug by blocking 1ZSX.