Abstract

Allicin and its derivatives are physiologically active molecules with many potential health benefits. It is not clear if they have a direct effect on DNA or protein. In order to elucidate the allicin and its derivatives’ effect on DNA base pairs, we investigated various complexes of base pair and its bridge with an allicin and its derivatives. The investigated allicin derivatives were (E)-Ajoene, (Z)-Ajoene, Amz121, and Bmz73 radicals. The DFT calculation results revealed that the investigated molecules would favor binding to bridge of the base pairs instead of directly affecting the base pairs. The Bmz73 radical could break DNA by change bonding in it because the Bmz73 radical significantly affected the P-O bond of the bridge of base pair.

1. Introduction

Infections that develop into chronic conditions are a fast-growing problem; especially it is important in the developed world [1, 2]. In relation to emerging multiresistant bacteria, development of antimicrobials and new treatment strategies of infections should be expected to become a high-priority research area [2]. Indeed, garlic has been used for centuries because of its therapeutic and health-promoting properties [3, 4]. Allicin is almost exclusively responsible for the antimicrobial activity of freshly crushed garlic [4]. Allicin is a physiologically active molecule with many potential health benefits [4] but there are many unanswered questions about its functions [1, 513]. Allicin is active against a broader spectrum of microorganisms than most of the commonly used antibiotics. For example, allicin is active against both Gram positive (e.g., Bacillus and Streptococcus) and Gram negative bacteria (e.g., Salmonella typhimurium and Agrobacterium tumefaciens), whereas penicillin is practically not effective against the latter [14, 15]. Allicin is also active against human pathogens that are resistant against certain antibiotics [1, 16]. Beyond its well documented strong antibacterial properties allicin also shows toxic effects towards fungal cells and is able to inhibit spore germination and hyphal growth in vivo and in vitro [1, 17]. Allicin showed promising activity against many plant-pathogenic fungal species, including pests, which are of economic importance also [4].

When allicin is inside the cell, the antibiotic efficiency depends upon reaching and reacting with its targets and on the importance of those targets to the cell [4]. Sulfonic acids have been studied in many chemical reactions [1820]. However, the transformation rate of allicin depends on the applied techniques of transformation, temperature, and solvents nature [21]. There are additional potent antithrombotic agent from garlic which are (E,Z)-Ajoene [2224] which have neuroprotective effects also [23]. The Ajoene structures are formed from a chemical reaction involving two allicin molecules. On the other hand the allicin molecule is unstable structure which forms radicals over chemical reactions depending on the environment. Depending on homolytic cleavage or heterolytic cleavage of the S-S bonds in the allicin there are allicin derivative radicals after dissociation of the S-S bond [25]. The most stable radicals are labeled as Amz121 and Bmz73 [25] (Figure 2). Allicin radicals are good candidates for possible DNA damage also [25].

Does allicin have an effect on DNA, RNA, and/or protein synthesis (Figure 1) [4]? In 1988, Feldberg et al. investigated the effect of allicin on the synthesis of DNA, RNA, and proteins in Salmonella typhimurium [26]. Unfortunately, it remains unclear if the effect observed really showed a direct decrease in DNA, RNA, and protein synthesis (Figure 1) [4]. The possible allicin attacks on proteins are over disulfide exchange reaction what can lead to misfolding or break-off of the protein. Another possible attack is directly to DNA which can be two types directly affecting DNA or by structures after disulfide exchange reaction what can lead to break-off or blocking DNA from its usage.

According to Durlak et al. [9] CPMD simulations and DFT and MP2 methodology analysis there are about 20 conformers of allicin which have quite similar energies. In solvent any of allicin conformer can easy become from one to the other type conformer and the trans-allicin structure does not exist [9]. According to other DFT calculations together with experiment the alkali metal ions can strengthen the S1S1′ bond in the allicin, while protonation or the loss of an electron weakens the S1S1′ bond [27]. In the gas phase the dehydration reaction mechanism takes place easier for the allicin radical cation than for the protonated allicin [25]. According to other DFT analyses [28] it is more difficult in the intramolecular to transfer H atom during decomposition of allicin in ethanol.

The allicin derivatives which may be important for biological activity are allicin, Ajoene structures, or allicin radicals while possible effects are on protein or on DNA. This work was dedicated for computational analysis for the possible DNA damage models (Figure 2). In the present study, in order to elucidate the allicin and its derivatives effect on DNA base pairs, the stable structures and their electronic properties for the complexes of G-C, A-T base pairs, and base pairs bridge with an allicin, (E)-Ajoene, (Z)-Ajoene, and allicin radicals (Figure 3) by using the DFT calculations in vacuum were investigated. The results revealed that the allicin, Ajoene, and its allicin radical prefer to interact with the bridge of the base pairs. The calculation results showed that allicin and its derivative molecules would favor binding to bridge of the base pairs instead of directly affecting the C-G or A-T base pairs.

2. Calculations and Investigated Structures

2.1. Calculation Details

Despite several chemical reaction studies using HandHLYP/6-311++G(d,p) [29], (U)B3LYP/6-31+G (d,p) [25], or B3LYP/6-31+G (d,p) [27] levels of DFT theory there is lack of allicin or allicin radicals modeling. The allicin complexes are used for modeling of the complexes in [9, 30]. In [9] for DFT calculations, B2PLYP/Def2-TZVPPD are used for single point calculations while optimization are performed using MP2/aug-cc-pVTZ. The vibrational analyses with M062X and PBE1PBE functional are not good enough according to experimental data [9]. The agreement with experimental data for on the dehydration reaction of the allicin calculations is achieved in [28] where (U)B3LYP/6-31+g(d,p) level calculation is used. The same calculation level (U)B3LYP/6-31+G(d,p) is used for theoretical studies on dissociation of the S-S bond in the allicin in [25] where results agree with the experimental data. The intermolecular interaction energy between allicin and solvent (pentane and ethanol) is calculated by DFT at the B3LYP/6-31+G(d,p) level in [27]. The B3LYP functional with 6-31G(d) [31], 6-311G(d,p) [3234], 6-311+G(2d,p) [32], 6-311++G(2d,p) [32], and cc-pVTZ [33] basis sets is good methods for modeling several small molecules as supermolecules in analyzing the proton transfer while the best results are by using cc-pVTZ basis set with B3LYP functional [32, 33]. In these calculations all allicin type structures and DNA parts were computed using DFT theory with B3LYP/cc-pVTZ. The calculations were performed using the Gaussian 09 program package (Rev. D.01) [35].

According to [36] the optimization calculations should be performed by increasing basis set. The calculations were done in such way. Zero-point energy (ZPVE) corrections were not applied for the modeling of the complexes. BSSE-corrected energies were not taken for the analysis of final 45 ground-state complexes because the results should not change significantly as it is shown in [33]. All stationary points were modeled in the gas phase (vacuum). Additionally, the allicin molecules were characterized by frequency analysis where frequency calculations having all real frequencies identify minimum structures.

There were analyses of chemical bonds of calculated complexes in the paper. The bond analysis was performed by analyzing Wiberg bond index. Primary focus was on the Wiberg bond index values between complex structures by looking for largest values. Also the distances between S atom of the interacting structure and nearest atom of the investigated structures were measured. Mulliken and NBO charges and HUMO-LUMO gap calculations for all the structures and complexes were performed.

Binding energies in the paper were calculated according to the following formula: where is final binding energy, is energy of optimized interacting and investigated structures complex, is optimized interacting structure energy, and is optimized investigated structure energy.

2.2. Complexes of the Interacting Structures with the DNA Parts

For the modeling in this paper the primary idea was that the S-S must be active part for the possible interactions. Thus it was enough to have one stable conformer which had open form of the S-S bond. The interacting structures were allicin, (E)-Ajoene, (Z)-Ajoene, and allicin radicals labeled as Amz121 and Bmz73 similar to [24] (Figure 3). Frequency calculations showed the trans-allicin structure had five imaginary frequencies. The cis-allicin structures had not imaginary frequencies what meant they were in the minimum. Before modeling complexes all structures where optimized separately. The primary DNA structure was taken from PDB bank [37]. The DNA was divided into three major parts: G-C base pair, A-T base pair, and bridge of base pairs as it is shown in Figure 2. After optimization of DNA structural parts (G-C base pair, A-T base pair, and bridge of base pairs) the final structures were compared with crystallographic data [37]. The structures were quite similar and the difference was up to 5% for the bond lengths, angles, and dihedral angles. In order to find the most preferable interacting site of possible allicin derivatives around the DNA, artificial initial model structures were made by patterns as shown in Figure 2. The S atom or S atoms from S-S bond were put as close to the DNA parts as possible. There were three major starting model positions for optimizations where the interacting structure was from different side of the investigated structure. In Figure 2 allicin cases with G-C base pair (Figure 2(a)), A-T base pair (Figure 2(b)), and bridge of base pairs (Figure 2(c)) are shown after the final optimizations. In Figure 2 the X sign represents the initial position before the optimization. Positions 1, 2, and 3 (Figure 2) were the patterns of artificial made initial model structures for the investigated complexes. The model structural parts were taken from separately optimized structures. Similar complexes were modeled and calculated with all interacting structures. The ground-state differences between same type complexes with different positions are shown in Table 1. There was no particular one position for interacting structure where the all complexes would be in the lowest ground state (Table 1).

3. Results and Discussion

3.1. Allicin Interaction with the DNA Structural Parts

Firstly, allicin molecule was put at three positions around G-C, A-T, and the bridge of the base pairs as shown in Figure 2. The results are presented in Table 1. The nearest position between allicin S and G-C part is 3.57 Å. The nearest distance to the A-T pair is 4.01 Å. The nearest O of the base pairs bridge from the allicin is 3.74 Å. In all cases the Wiberg bond index between dimer parts was almost zero: the largest value was 0.015 for S and H atoms in the case with G-C base pair (Table 1). The S atoms had 0.48 and −0.13 Mulliken charges and 1.09 and −0.04 NBO charges in vacuum allicin. The charge values for the S atoms were similar in complexes and they varied in 0.5 range. The HOMO-LUMO gaps in all cases were about 4 eV. The largest variations there were in the bridge of the base pairs case: the HOMO-LUMO gap was 4.79 eV in position 2 and 3.66 eV in position 1. There was suggestion that the something would change if we put several allicin molecules. So two and three allicin complexes were used but there were no significant differences in results. There were several cases when one of the three types of allicin moved between A-T and G-C groups but this is unrealistic: in investigated model the A-T (or G-C) groups could freely move in the calculation space what is not true for the real DNA sequence. The binding energies between allicin and G-C base pairs varied form 14.87 kJ/mol to 31.75 kJ/mol. Similar results were between allicin and A-T base pairs where the binding energies varied from 21.48 kJ/mol to 28.32 kJ/mol. The binding energies were almost twice stronger for the base pair bridge then for G-C or A-T base pairs, which varied from 42.8 kJ/mol to 68.61 kJ/mol. Thus the results suggested that allicin molecules would favor binding to bridge of the base pairs and the calculated binding energies could be up to 70 kJ/mol. The calculation results showed that allicin molecule did not change significantly bonding between atoms of the DNA structural parts.

3.2. Ajoene Interaction with the DNA Structural Parts

With the same procedures as with allicin molecule, the calculations were performed for Ajoene molecules and those molecules are labeled as (E)-Ajoene and (Z)-Ajoene (Figure 3). The modeling was performed according to Figure 2 schema and results are presented in Table 1 where nearest positions between the S atom of Ajoene and O and N atoms of G-C, A-T, and the base pairs bridge are shown. In all cases the Wiberg bond index between dimer parts were almost zero: the largest value was 0.028 between H and O atoms in the case of (E)-Ajoene and bridge of the base pair complex (Table 1). The S atoms of (E)-Ajoene had zero and −0.06 Mulliken charges and 0.14 and 0.08 NBO charges. The charge values for the S atoms of (E)-Ajoene were similar in complexes and they varied in 0.05 range. The S atoms of (Z)-Ajoene had −0.05 and −0.07 Mulliken charges and 0.09 and 0.08 NBO charges. The charge values for the S atoms of (Z)-Ajoene were similar in complexes and they varied in 0.5 range. The HOMO-LUMO gaps for complexes varied from 3.37 eV to 4.75 eV. The smallest HOMO-LUMO 3.37 eV gap was for (E)-Ajoene with C-G base pair while the largest HOMO-LUMO 4.75 eV gap was for (E)-Ajoene with A-T base pair. The binding energies between Ajoene molecules and G-C or A-T base pairs varied from 8 kJ/mol to 38 kJ/mol where the largest bonding energy was between (Z)-Ajoene and G-C base pair which was 37.68 kJ/mol. The binding energies were almost twice stronger for the base pair bridge and then for the G-C or A-T base pairs which were from 46.57 kJ/mol to 85.82 kJ/mol. The results suggested that Ajoene molecules similar to allicin would favor binding to bridge of the base pairs and the binding was up to 85.82 kJ/mol. The calculation results showed that the Ajoene molecules did not change significantly bonding of the DNA structural parts.

3.3. Allicin Radicals Interaction with the DNA Structural Parts

Additional similar calculations were performed using allicin radicals: Amz121 and Bmz73 (Table 1). The Amz121 structure and G-C or A-T base pairs had binding energy from 251.51 kJ/mol to 314.62 kJ/mol which were very strong. The Wiberg bond indexes were from 0.476 to 0.638. There were two types of bonding: between S and O (3 positions in G-C base pair; all positions in A-T base pair cases); between S and N (1 and 2 positions in G-C base pair). The Amz121 structure with bridge of the base pairs had binding energy twice larger: from 683.4 kJ/mol to 686.12 kJ/mol. Also the Wiberg bond indexes were larger which were about 0.65 and they were between S and O atoms. The S atom of Amz121 had −0.06 Mulliken charge and 0.14 NBO charge. In all cases the S atom donated electron because the charges became positive: S atoms of Amz121 in complexes had 0.62–0.66 Mulliken charges and 1.34–1.47 NBO charges. The smallest HOMO-LUMO gaps from 2.96 eV to 3.76 eV were for complexes with Amz121 and G-C base pair structures while the largest HOMO-LUMO gaps from 5.9 eV to 6.24 eV were for complexes with Amz121 and bridge of the base pair structures.

The Bmz73 structure and G-C or A-T base pairs had binding energy from 375.12 kJ/mol to 433.7 kJ/mol. The Wiberg bond indexes were from 0.77 to 0.91. There were two types of bonding: between S and O (3 positions in G-C base pair; 1 position in A-T base pair cases); between S and N (1 and 2 positions in G-C base pair; 2 and 3 positions in A-T base pair). The Bmz73 structure with bridge of the base pairs had binding energy almost twice larger: from 674 kJ/mol to 802 kJ/mol. The Wiberg bond indexes were about 0.8 and they were between S and O atoms. The S atom of Bmz73 had 0.32 Mulliken charge and 0.56 NBO charge. The NBO charges on the S atom were quite similar for the A-T base pair complex and for the bridge of the base pair complex: S atoms of Bmz73 had 0.07–0.16 Mulliken charges and 0.42–0.52 NBO charges. In C-G base pair complex case the S atom the charges became less positive: S atoms of Bmz73 in C-G base pair complex had 0.03–0.17 Mulliken charges, 0.09–0.37 NBO charges. The smallest HOMO-LUMO gaps from 3.67 eV to 4.38 eV were for complexes with Bmz73 and G-C base pair structures while the largest HOMO-LUMO gaps from 5.28 eV to 5.43 eV were for complexes with Bmz73 and bridge of the base pair structures.

Bmz73 radical significantly changes P-O bond length of the bridge from 1.68 Å to 2.18 Å in 1 position complex. The P-O bond enthalpy was 389.84 kJ/mol without Bmz73 radical. The P-O bond enthalpy became 9 times smaller (42.52 kJ/mol) in complex of Bmz73 and bridge of the pair base. The P-O bond Wiberg bond index was 0.58 for the bridge of the base pair without Bmz73. The P-O bond Wiberg bond index became 0.15 in complex of the Bmz73 and the bridge of the base pair. In this case the O atom of the bridge of the base pair had new 0.79 Wiberg bond index with the S atom of the Bmz73 radical. Thus Bmz73 can break or significantly affect the P-O bonding in the bridge of DNA.

4. Conclusions

The DFT calculations with B3LYP/cc-pVTZ were performed between interacting structures and DNA parts (G-C and A-T base pairs and the base pairs bridge) where interacting structures were allicin, Ajoene structures, and allicin radicals. The allicin analyses suggested the allicin structure was in minimum when it has cis-configuration. In all case of the allicin, Ajoene molecules with the DNA parts had almost zero Wiberg bond indexes. Allicin radicals had Wiberg bond indexes from 0.5 to 0.9 with the DNA structural parts. In Bmz73 radical case it significantly affected the P-O bonding of the bridge of base pair. This means that the Bmz73 radical could break DNA by change bonding in it. All calculation results showed that allicin and its derivative molecules would favor binding to the bridge of the base pairs instead of directly affecting C-G or A-T base pairs.

Conflicts of Interest

The author declares that there are no conflicts of interest regarding the publication of this article.

Acknowledgments

The public access supercomputer from the High Performance Computing Center (HPCC) of the Lithuanian National Center of Physical and Technology Sciences (NCPTS) at Vilnius University was used. Mindaugas Macernis acknowledge the prize for research of young scientists by the Lithuanian Academy of Science. This work is funded by the Research Council of Lithuania (LMT Grant no. MIP-080/2015).