Antioxidant Properties of Lapachol and Its Derivatives and Their Ability to Chelate Iron (II) Cation: DFT and QTAIM Studies

Pajoudoro, Djafarou Ngouh; Lissouck, Daniel; Ateba Amana, Baruch; Mfomo, Joseph Zobo; Abdallah, A. E. B.; Toze, Alfred Aristide Flavien; Mama, Désiré Bikele

doi:https://doi.org/10.1155/2020/2103239

Bioinorganic Chemistry and Applications

On this page

Abstract Introduction Results and Discussion Conclusion Data Availability Conflicts of Interest Acknowledgments Supplementary Materials References Copyright Related Articles

Research Article | Open Access

Volume 2020 | Article ID 2103239 | https://doi.org/10.1155/2020/2103239

Antioxidant Properties of Lapachol and Its Derivatives and Their Ability to Chelate Iron (II) Cation: DFT and QTAIM Studies

Djafarou Ngouh Pajoudoro,¹Daniel Lissouck,²Baruch Ateba Amana,¹Joseph Zobo Mfomo,³A. E. B. Abdallah,²Alfred Aristide Flavien Toze,¹and Désiré Bikele Mama¹

Academic Editor: Concepción López

Received13 Dec 2018

Revised09 Jun 2019

Accepted21 Aug 2019

Published31 Mar 2020

Abstract

The elucidation of the complexation of lapachol and its derivatives to Fe²⁺ cation has been done using the density functional theory (DFT). This complexation has been limited to bidentate and tridentate to Fe²⁺ cation. Geometry optimizations have been implemented in gas and solution phase (water, acetonitrile, chlorobenzene, benzene, and toluene) for ligands at B3LYP/6-311++G (d,p) level of theory using B3LYP/6-31+G(d,p) optimized data as starting point. But, the geometrical optimizations in solution phase of the 22 complexes analyzed of lapachol and its derivatives to Fe²⁺ cation were restricted to acetonitrile and benzene. The complexation energy and the metal ion affinity (MIA) have also been calculated using the B3LYP method. The results obtained indicated a proportionality between the MIA values and the retained charge on Fe²⁺ cation for k²-(O₁,O₂) modes. But, an inverse proportionality has been yielded between these two parameters for k³-(O₂, C=C) tridentate modes. For k³-(O₃,C=C) tridentate mode coordination, the higher stability has been obtained. In this latter tridentate coordination in gas phase, the topological analysis of complexes exhibits the fact that the electron density is concentrated between the O₃ oxygen atom of the ligand attached to Fe²⁺ and this metal cation. Moreover, the hydrogen bond strength calculated for isolated ligands (situated between 23.92 and 30.15 kJ/mol) is in the range of normal HBs. Collectively, all the complexation processes have shown to be highly exothermic. Our results have also shown that the electron extraction from Fe²⁺...La_i complexes is more difficult compared to that from free ligands.

1. Introduction

Lapachol (2-hydroxy-3-(3-methylbut-2-enyl) naphthalene-1,4 dione is a natural quinone isolated from Bignoniaceae family found in Brazil [1] (more specifically from bark extracts of Tabebuia sarratifolia [2]). The interest of our attention resulted from its extensive range of biological and pharmaceutic activities [3–5]: antitumor, antibiotic, antimalaria, trypanocidal, and leishmanicidal activities. Its ability to interfere with the bioactivities of topoisomerases (set of enzymes critical for DNA replication in the cell) justifies the intense investigation of its derivatives [6]. Moreover, its higher activity against carcinosarcoma tumor (Walker 256) [7] has prompted Hartwell and Abbot to examine the efficiency of its derivatives. In the same vein, studies on the activity against carcinosarcoma (Walker 256) and that on the structural modifications of these derivatives have been reported experimentally and theoretically [8]. Silva et al. [9] suggested that this evident cytotoxic property is due to the aromaticity of naphthalenic ring of the molecular systems examined. The antioxidant ability of lapachol has been experimentally confirmed by Sadananda et al. through the isolation and characterization of lapachol from different parts of Tabebuia argentea [10]. The excellent chelating keto-enol function of lapachol and its derivatives has been proved [11, 12]. In addition, the enhancement of antiplasmodial activity has been observed by the ruthenium chelation of free lapachol by Barbosa et al. [13]. Other modifications of the activity of lapachol resulting from the metal chelation have been experimentally tested using the manganese (II) lapacholate polymer. In the same vein, Tabrizi et al. have synthesized and characterized 3^rd row metal (Cu, Co, and Ni) complexes of lapachol that have exhibited a promising antitumor activity [14]. These three complexes were previously characterized by cyclic voltammetry [15]. In addition, the antimony (V) and bismuth (V) complexes of lapachol have shown a significant cytotoxic activity [16]. Despite this fact, a deep analysis of the influence of the metal chelation on the biological activities of free lapachol is lacking. Such an examination (of the antioxidant activity for instance) may give plausible explanations of the effectiveness of the intake of lapachol tea for skin antiaging effect [17]. To the best of our knowledge, few computational studies [8] concerning the antioxidant activity of lapachol and its derivatives have been reported. Until now, there is no study on the impact of the metal chelation of these molecular systems on their antioxidant powers.

The characterization of the antioxidant activity for large number of antioxidants RX-H, [18–22] has been done by the calculation of the RX-H (X = donor atom such as O, N, ...) bond dissociation enthalpy (BDE). From these calculations, the prediction of the relationship structure-activity and the design of novel possible antioxidant have been carried out. The radical formed from the hydrogen atom transfer (HAT) can be obtained through two-step mechanisms: the single electron transfer-proton transfer (SET-PT) mechanism (expressed using ionization potential (IP) and proton dissociation enthalpy (PDE) of the RX-H^•+ radical cation) and the sequential proton loss-electron transfer (SPLET) mechanism that matches up with the proton affinity (PA) of the anion RX^- and the reaction enthalpy of abstraction of an electron from this anion denoted as ETE (electron transfer enthalpy). The ability of the potential antioxidant to chelate transition metal ions, precisely iron and copper cations (in order to prevent the damage of oxidative stress from the participation in free radicals generation [23]), is also considered as antioxidant mechanism. The metal chelation reduces the production of toxic hydroxyl radical (OH•) from Fenton reaction:

The current study reports the results of a computational examination of the antioxidant power of molecular system in different media (gas and solution phase) and the influence of the iron (II) chelation. A systematic study of reaction enthalpies related to SET mechanism for lapachol and its seven derivatives (Figure 1) denoted as La_n (n = 1–8) is performed. A direct verification of accuracy of calculated reaction enthalpies cannot be performed due to the nonavailability of their experimental homologues. The calculated physiochemical antioxidant descriptors (IP) for this molecular system are compared to that of the classical antioxidants (gallic acid, caffeic acid, ferulic acid, ascorbic acid (vitamin C), and trolox). The metal cation affinity for Fe²⁺ cation and the impact of this metal chelation on the antioxidant activities through the examination of the SET mechanism are also inserted. This last point must be helpful for the clarification of the properties of La_n−Fe²⁺ that may facilitate the selection of appropriate chelating and scavenging agent for the chelating therapy of Alzheimer [24].

2. Computational Details

The geometry optimizations of lapachol and its seven derivatives molecular systems were performed at the B3LYP/6-311++G(d,p) level using Gaussian 09 [25] program from B3LYP/6-31+G(d,p) data considered as starting point. Collectively, these two series of geometry optimizations were carried out with no symmetry constraints in the ground state. The frequency calculations have been done to select optimized structures without negative frequencies. Free ligands (RX-H) were optimized in gas phase and in five solvents (water, acetonitrile, chlorobenzene, benzene, and toluene). For these latter cases, the integral equation formalism of polarizable continuum model (IEF-PCM) [26, 27] was used. The ionization potential (IP) of these free ligands (enthalpy change of the reaction (2)) is considered as descriptor of the SET (single electron transfer) antioxidant mechanism:

Higher antioxidant activity of the molecular library is the molecule that exhibits lower (IP) values. This molecular descriptor was therefore calculated using the following formula

The gas phase enthalpy for electron (e⁻) is equal to 3.145 kJ·mol⁻¹ [28]. We have adopted the solution phase enthalpies proposed by Markovic et al. [29]. The ionization potential free energy (IPFE) which is the reaction free energies of reaction (2) is calculated according to the following formula:where the calculated gas free energy of an electron G(e⁻) is equal to −3.72 kJ/mol [30]. Electron solvation free energies were taken from the literature [29].

From the enthalpy change of the complexation (5) of lapachol and its seven derivatives by Fe²⁺, the La_i−Fe²⁺ binding energy (i = 1–8) denoted by E_int can be predicted:where , , and represent, respectively, total energy of the optimized complex (), isolated ligand (), and free iron (II) cation. The enthalpy (ΔH) and free enthalpy ((ΔG) change were also predicted. The ionization potential of these complexes (IP_C) as enthalpy change of the reaction (7) was calculated according to (8). This fact has favored the assessment of such a complexation on the antioxidant mechanism in various media through its influence on SET mechanism:where H(La_i...Feⁿ⁺) is the formation enthalpy of the complex La_i...Feⁿ⁺ (n = 2 and 3).

In general, the ionization values can be directly obtained as the negative of the eigenvalue of highest occupied molecular orbital (HOMO) (IP for La_i ligands and IP_C for La_i…Fe²⁺ complexes) as recommended by Koopman’s theorem [31]. But, this procedure has been adopted exceptionally for La_i…Fe²⁺ complexes (IP_C) due to the difficult estimation of the geometry changes provoked by the reflections or the relaxations of ligand fragment during the geometrical optimization progression.

The quantum mechanics atom in molecule (QMAIM) theory was consequently performed from the B3LYP/6-311++G(d,p) optimized structures to investigate these two parameters (bonding properties and bond critical points). More specifically, the clarification of the strength of hydrogen bond interactions (in isolated ligands) and metal-ligand interactions in complexes has been done using the Multiwfn program [32]. For bonding interactions, the indicators are electron densities and its Laplacian at bond critical points (CP). The total number of these CPs obtained is in accordance with the Poincare–Hopf rule [33]. For deeper examination, the density of the total energy of electrons (H) defined as the sum of the Lagrangian kinetic electron density (G) and the potential electron density (V) at bond critical points (BCPs) has been estimated:where G_BCP and V_BCP are, respectively, defined in (10) and (11):

The interatomic interaction energy (denoted as E_in) in isolated ligands and in complexes was predicted by Espinosa approach [34]:

3. Results and Discussion

3.1. General Considerations on Isolated Ligands

The optimized structure for lapachol (La₁) is characterized by planarity of the naphthalenic ring. This plane also contained the carbon atom through which the 3-methylbut-2-enyl substituent is connected to this ring. But the whole molecular system is not planar due to the free rotation about 113.7° around the C-C single bond bound to the C=C double bond. The substitution of the two methyl groups by the bromide atom on the lateral chain (La₃) provokes an augmentation of this rotation angle about 11°. A migration of the C=C double bond toward the naphthalenic ring of the lapachol leads to an optimized structure (La₂) characterized by the C=C double bond contained in the plane formed by this ring. Such rigidity induces the symmetrical position of the two methyl groups according to the preceding plan. The substitution of the unsaturated substituent by the 3-methylbutyl (La₄) or 3-hydroxy-3-methylbutyl substituent has yielded optimized structures with a lateral chain almost perpendicular to naphthalenic ring. On the whole, optimized structures obtained in gas phase are characterized by the formation of hydrogen bond (HB) between the hydrogen atom (H₁) of hydroxyl group (H₁-O₂) and the oxygen atom of the carbonyl group (C₁=O₁). This formation of hydrogen bonds is in agreement with the cutoff definition of X-H...A hydrogen bonds proposed by Steiner [35] (H...A > 3.0 Å and X-H...A angle > 110°: Table 1S).

3.2. Metal Chelation Mechanism

Our results reveal the fact that Fe²⁺ cation is bound to La_1-8 ligand through bidentate or tridentate connection. The bidentate attachment (denoted k²-O₁,O₂) is observed when Fe²⁺ cation is placed in vicinity of sp² oxygen (O₁) and sp³ oxygen (O₂) of the O₂-H₁ group (Figure 2). The optimized structures of Fe²⁺...La_i complexes obtained denoted as aLa_i (i = 1–8) exhibited a Fe²⁺ cation contained in the plane formed by the naphthalenic ring. This structure is characterized by Fe-O₁ bond distances longer than Fe-O₂ ones (average difference is 0.087 Å). Collectively, the optimization of k²-O₁, O₂ complexes leads to structures where the hydrogen atom (H₁) shifts to be 10.2–29.4° off-plane. We have also integrated two possibilities to bind Fe²⁺ cation for La₅ ligand: (i) to bind to Fe²⁺ cation to sp³ O₂ oxygen atom of the O₂-H₁ group and sp³ O oxygen atom of the O-H group of the R substituent of La₅ ligand (k²-O₃,O: ala₅₂) and (ii) to bind to Fe²⁺ cation to sp² O₃ oxygen atom of the naphthalenic system and sp³ O oxygen atom of the O-H group of the R substituent of the ligand concerned (k²-O₃,O: ala₅₃). In the former case, the breakage of O₁⋯H₁-O₂ hydrogen bond is not observed, but the H₁ proton is transferred to O₁ and then becomes closer to O₁ atom than to O₂ atom. The simultaneous binding of Fe²⁺ cation to sp² O₂ oxygen atom of the O₂-H₁ group and C=C π bond of the substituent R (k³-O₂,C=C) (Figure 2) has also been observed. Similar tridentate attachment (k³-O₃,C=C) of Fe²⁺ cation to sp² O₃ oxygen of the naphthalenic moiety and C=C π bond of the same La_1-8 ligand is possible. For sake of clarity, these complexes have been, respectively, symbolized by bLa_i and cLa_i (i = 1-8). From the input geometry of the bLa₂ complex, the optimized structure is characterized by Fe²⁺ ion bound to La₂ ligand through a bidentate attachment. In order to evaluate the importance of hydrogen bonds, we have built initial geometrical structures with hydrogen bond or not (in both cases). The introspection of optimized structures obtained reveals the fact that the variation of Fe²⁺-ligand bonds (0.010–0.014 Å) is minor. This fact showed that the insertion of the O₁⋯H₁-O₂ hydrogen bond in the kLa_i (k = b,c) does not give any additional information on geometrical parameters. On the whole, Fe-O_i (i = 1, 2) bond lengths are longer than Fe-C homologues.

For La₃ ligand, the lowest total energy corresponding to the cLa₃ complex indicates the fact that k³-(O₃,C=C) coordination mode yields the most stable complex. Similar observations have also been obtained for other ligands (La₂, La₆, and La₇) in gas phase. This higher stability obtained for such a coordination mode is attributed to simultaneous attachment of Fe²⁺ cation to C=C and sp² O₃ oxygen atom. The latter exhibits greater electron donor nature compared to sp³ O₂ oxygen atom. This fact is in accordance to Cu²⁺ binding to 2-methyl-1-(2,4,6-trihydroxy-3-(3-hydroxy-3,7-dimethyloct-6-enyl)phenyl)propan-1-one, called Hyperjovinol extract from Hypericum jovis [36]. The significant relative energy obtained for k²-O₁,O₂ confirmed the augmentation of the stability of complexes in multidentate chelation of Fe(II) previously enlightened by Mwadham et al. [37]. The higher stability of aLa₅₃ compared to that of other Fe⋯La₅ complexes (aLa₅ (219,782 kJ/mol higher in total energy than aLa₅₃) and aLa₅₂ (28,209 kJ/mol higher in total energy than aLa₅₃) also displayed the sensitive influence of the type of electron donor group on which the Fe²⁺ cation is connected. In order to have much improved understanding of the stability, the interaction energies, interaction enthalpies, and the Gibbs energies of complexation have been calculated at B3LYP/6-311++G(d,p) level. The obtained results in the two solvents adopted (acetonitrile and benzene) for complexation and that obtained for gas phase are shown in Table 1. Collectively, all the interaction energies are highly negative, demonstrating that the complexation process is highly exothermic. The sensitive drop of binding energies observed when passing from gas phase to solvent used is attributed to solvent effect that obstructs the interaction between Fe²⁺ cation and La_i ligands. This drop is more significant for acetonitrile. For a more measureable description of the ligands’ affinity toward Fe²⁺ cation, Table 1 displays that the selectivity on the basis of the complexation free energy is very versatile when passing to one coordination to another: an increasing ligand affinity is ordered as La₄ < La₅ < La₃ < La₂ < La₆ < La₇ (for k²-O₁,O₂ mode), La₇ < La₂ < La₁ < La₃ < La₆ (for k³-O₂, C=C mode), and La₂ < La₃ < La₁ < La₆ < La₇ (for k²-O₃, C=C mode) in gas phase. The versatility observed is due to diverse geometry changes of La_i ligands within the geometrical structures of Fe²⁺...La_i complexes during the optimization geometry. On the whole, the negative values of complexation free energy support the finding that this process is spontaneous in gas phase and two solvents used (acetonitrile and benzene).

Figure 3 indicates a possible correlation between the retained NBO charge on the Fe atom and the metal ion affinity (MIA) which is assumed to be the negative of the enthalpy change during the complexation process: in the case of k²-(O₁,O₂) mode, MIA values vary proportionately with the retained charge on Fe²⁺ cation. On the other hand, for k³-(O₂, C=C) mode, MIA values vary inversely with the retained charge on Fe²⁺ cation. In the both cases, an exception is made for aLa_i (i = 6 and 7) complexes. Figure 4 exhibits the fact that there is no correlation between MIA values and retained charge on Fe²⁺ cation for k³-(O₃, C=C) complexes. The values of the net charge of Fe²⁺ cation in the range + (1.432–0.513)e approve the fact that the Fe²⁺ cation experiences a reduction while the La_i ligands undergo an oxidation. The electron transfer from the latter to Fe²⁺ cation which augments with the attachment of Fe²⁺ cation to C=C confirms the oxidation role of the metal cation toward the molecular system examined. This result is in accordance with previous observations done on 1,2,3-triazol derivatives [38], curcumin [24], and Hyperjovinol [36]. From Table 1, one could notice that the diminution of the electron transfer is more sensitive in acetonitrile. This can be attributed to its polar property. On the whole, this is a suggestion that electron transfer from ligand to Fe²⁺ cation happens during the formation of coordination bond.

(a)

(b)

3.3. Topological Properties

The values of Laplacian of electron density () calculated (BCPs) for X-H...A bond critical points (BCPs) are a second derivative. The local electronic charge (depletion) is observed in the regions in which the values of Laplacian of electron density are positive (Table 1S). This reveals the fact that the interaction is characterized by the depletion of the density in the region of contact of hydrogen atom H₁ and oxygen atom O_1. The local excess in the kinetic energy is predominant in such an interaction. The average ratio greater than 1 in various media (gas (1.15), water (1.14), acetonitrile (1.18), benzene (1.19), and toluene (1.16)) confirms the noncovalent bonding interaction. The formation of hydrogen bonds is confirmed by positive values of small value of and positive values of for H_BCP (Table 1S). Our data exhibit a linear correlation between these two parameters and the bond distance (Figure 1S) in gas phase. The correlation coefficients (R²) were, respectively, equal to 0.98836 and 0.98096. This trend has been obtained in solution phase. In benzene for instance, these coefficients become, respectively, equal to 0.99951 and 0.99592. The inverse proportionality between hydrogen bond distances and these two topological parameters indicates an attenuation of orbital overlap that provokes an augmentation of the interatomic distances. As it is seen in Table 1S, HB strengths calculated in the 23.92–30.15 kJ/mol range are in agreement with that of normal HBs (8.36–41.8 kJ/mol) [39].

The nature of interactions between the donor ligand and Fe²⁺ cation in the Fe²⁺...La_i complexes has been examined using the topological parameters listed in Table 2: the electron density (ρ(r)), its Laplacian , kinetic energy density , potential energy density , total energy of electron (H(r), and ellipticity of the metal-ligand bond (ɛ)) at their bond critical points (BCP). The electron density of BCP for coordination bonds is found to be varied in the range: 0.069–0.111 Å^′, 0.119–0.330 Å^′ and 0.049–0.330 Å^′, respectively, for k²-(O₁,O₂), k²-(O_i,O) and k³-(O_i, C=C): i = 2, 3. A survey of results of Table 2 demonstrates that the values of Laplacian of the electron density are collectively positive with exception made for aLa₅₂ complex. This implies the preeminent interaction of the contraction of the electron density ρ(r) toward each nucleus. Therefore, the parallel gradient and the curvature of ρ(r) are great. The net forces of repulsion act on the nuclei. It is noteworthy to mention that ρ(r) values of Fe–O_i (i = 2, 3) bonds are higher than those Fe–C in the case of k³-(O₂, C=C) and k³-(O₃,C=C) coordination mode, supporting the fact that the bond distances of the latter are greatest.

The calculated average ratio is equal to 0.9670, 0.8108, and 0.8081 for k²-(O₁,O₂), k³-(O₂, C=C), and k³-(O₃,C=C) coordination mode, respectively. These ratios comprised between 0.5 and 1 jointed to low values of and collective positive values of indicate the partially covalent nature of the metal-ligand bonding in gas phase. Similar results have been computationally obtained by Nkungli and Ghogomu through an investigation of the binding Fe³⁺ cation (protoporphyrin IX) to 4-methoxyacetophenone thiosemicarbazone) [40]. Additive investigations based on the values at BCP (, , and : eigenvalue of the Hessian matrix [40]) have been done to have further qualitative information on type of metal-ligand bonding interaction. A bond critical point is represented by one positive and two negative eigenvalues. For k²-(O₂,O) coordination mode, the values that range from 0.552 to 0.621 confirm the noncovalent nature of the Fe²⁺-ligand interaction. Contrary to k²-(O₃,C=C) mode in which the for Fe²⁺-O₃ interaction (in almost complexes): 8,605 (cLa₁), 6,341 (cLa₃), 5,100 (cLa₆) and 6,641 (cLa₇).

Therefore, the concentration of electron density is located between the ligand atoms attached to Fe²⁺ and the latter. For all interaction analyzed, the negative values of obtained from our calculations exhibited the accumulation of electron density at the bond critical point. The highest value of ellipticity of the bond ɛ () is obtained for Fe-O₃ bond in aLa₅₃ complex ((ɛ = 1.893919).

3.4. Electron Transfer Mechanism

The calculated IP values obtained from equation (3) and Koopman’s theorem for isolated ligand examined in various media are listed in Table 2S. The introspection of Table 2S indicates that IP values obtained for direct procedure (Koopman’s theorem [31]) are lower in gas phase. The difference in IP ranges from 58.362 to 137.224 kJ/mol. Similar observations have been found out for nonpolar solvent (benzene and toluene) with lower differences in IP: an average difference in IP is 32.451 and 27.367 kJ/mol, respectively, for benzene and toluene. Nonetheless, higher average differences in IP (155.622 and 135,733 for water and acetonitrile respectively) have been yielded in the benefit of that obtained from Koopman’s theorem. On the whole, we can conclude that the IP gap between results obtained from one procedure to another varies with the dielectric constant of the solvent. This may be attributed to the variation of the enthalpy of the solvation of the electron. Arbitrarily, we concentrate the further analysis of the antioxidant power of isolated ligand according to equation (3). From Table 2S, it can be perceived that the shift of the C=C double bond toward the naphthalenic ring of La₁ leads to a 9.84 kJ·mol⁻¹ drop in IP (La₂). This can be attributed to the extension of the conjugation system. This extension associated to conjugation of the π-electron has previously shown their influence on IP values of polyphenolic deoxybenzoins [41]. Recently, the effect of the double bond position on antioxidant power of 2-phenyl-benzofuran has been experimentally enlightened by Li et al. [42]. On the contrary, the substitution of the two methyl groups of La₁ by bromide atoms provokes a 36.66 kJ/mol augmentation in IP (La₃). In the same vein, the saturation of the C=C double bond of La₁ results in a 46.70 kJ/mol increase in IP (La₄). This augmentation illustrates the important role of this C=C double in improving antioxidant capacity. This is in agreement with enhancement of biological effect by the presence of C=C bond demonstrated by Chen [43]. The difference in IP for La₂ and La₆ equal to 4.776 kJ/mol exhibits the effect of steric hindrance on antioxidant activity. The alteration of the antioxidant activity of some Schiff base ligands and their copper (II) complexes by this steric hindrance has been reviewed by Salga et al. [44]. Similarly, Sun et al. have experimentally demonstrated that the antioxidant activity of flavonoids isolated from Hebei balmy chrysanthemum is sensitive to the steric hindrance of glucoside [45]. From Table 2S, it should be observed that molecules with lower IP value than that of gallic acid (792.9 kJ/mol [46]) and caffeic acid (762.9 kJ/mol) calculated at the same level of theory are La₁, La₂, and La₆. Furthermore, the IPs for La_1-3 and La₆ are lower than that of ascorbic acid in gas phase and in water (804. 000 kJ/mol (gas phase) and 517.484 (water) kJ/mol at B3LYP/6-311++G(2d,2p)//B3LYP/6-31G(d,p) [41]. On the whole, this denotes the higher antioxidant power of molecules (La_1-3 and La₆) compared to that of ascorbic acid in these two media.

In solution phase, IP values are lower than the matching gas phase values. The average differences between IP in gas and that in other media are equal to 267.611 kJ/mol (water), 253.931 kJ/mol (acetonitrile), 90.930 kJ/mol (benzene), and 97.808 kJ/mol (toluene). This demonstrates the sensitivity of the cation radical forms of isolated ligand (La_i) to the polarity of different solvents used [21, 47]. This may be attributed to the lowest value of the enthalpies of solvation of the electron. The increasing order of IP in media is water < acetonitrile < toluene < benzene < gas. From Table 2S, the IPFE values calculated are indicating that SET mechanism is not spontaneous regardless of the environment considered. The computed IP_C values of Fe²⁺ chelates in three media (gas, water, and benzene) are reported according to Koopman’s theorem (as above indicated) in Table 3S. The comparison of this IPC values of complexes to those of free ligands (Tables 2S and 3S) displays the facts that the values of the former are greatly higher. The difference in IP is in the following range: 107.934–712.718 (gas phase), 36.810–632.220 (acetonitrile), and 280.561–407.583 kJ/mol (benzene).

This is an indication that the electron extraction from Fe²⁺...La_i complexes is more difficult than in that from free ligand. This close-fitting binding of electron in cationic Fe²⁺complex system is similar to the computational prediction done for Fe²⁺-caffeic acid phenethyl ester derivative complexes (in gas phase) [48] and Fe²⁺-curcumin complexes (in gas and DMSO solvent phase) [24]. In the latter case, other divalent metal cations (Mn²⁺ and Zn²⁺) have shown their ability to enhance the antioxidant activity of isolated curcumin.

4. Conclusion

In the present paper, we have presented the computational calculations at B3LYP/6-311++G(d,p) level of theory of the lapachol and seven of its derivatives. The integral equation formalism of polarizable continuum model (IEF-PCM) was used to examine for predictions in solution phase (water, acetonitrile, benzene, and toluene). We limited our examination on SET mechanism. It turns out that this mechanism is not spontaneous regardless of the environment considered. The increasing order of IP yielded in media is ordered as water < acetonitrile < toluene < benzene < gas. IPs of Fe²⁺ complexes are greater than those of free isolated ligands. This fact exhibits the difficult electron extraction from Fe²⁺...La_i complexes. In general, the electron transfer from ligand to Fe²⁺ cation is more pronounced in tridentate mode coordination. In isolated ligands, the hydrogen bonding interactions have shown to be noncovalent in all the media, whereas the metal-ligand bonding is intermediate type interaction in gas phase. On the whole, the complexation process has shown to be highly exothermic. The optimized structures of Fe²⁺...La_i complexes at B3LYP/6-311++G(d,p) are characterized by Fe-O_i (i = 1, 2) bond lengths that are longer than Fe-C homologues in all the media.

Data Availability

The data used to support the findings of this study are available.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The authors are grateful to Dr. Cyril Assongo Kenfack (Cepamoq, University of Douala, Cameroon) for helpful discussion. Scientific support and access to computing resources from the High Performance Computing Center of the University of Strasbourg (funded by Equipex: Programme d’Investissements d’Avenir) is also gratefully acknowledged.

Supplementary Materials

Table 1S: intermolecular hydrogen bond parameters; BCP parameters: electron density ( ) and Laplacian of the electron density ( ) and hydrogen bond strength (kJ/mol) for lapachol and its derivatives at B3LYP/6-311++G(d,p)//B3LYP/6-31+G(d,p) level in other media. Table 2S: values in kJ/mol of ionization potential (bold face numbers are calculated using Koopman theorem) and ionization potential free energy of free ligands in various media calculated at B3LYP/6-311++G(d,p). Table 3S: values in kJ/mol of ionization potential and ionization potential free energy of Fe2± ligand complexes in various media calculated at B3LYP/6- 311++G(d,p)//B3LYP/6-31+G(d,p). Figure 1S: linear dependence of H-bond distances at bond critical points in the gas phase with electron density and Laplacian of electron density. (Supplementary Materials)

References

S. G. C. Fonseca, R. M. C. Braga, and D. P. Santana, “Lapachol—chemistry, pharmacology and assay methods,” Revista Brasileira de Ciências Farmacêuticas, vol. 84, pp. 9–16, 2003.
View at: Google Scholar
M. Fernandes-de-Oliveira, “Contribuição ao conhecimento químico das espécies Tabebuia serratifolia Nichols e Tabebuia rosa bertol,” Universidade Federal do Ceará, Fortaleza, Ceará, Brasil, 2000, Ph.D. thesis.
View at: Google Scholar
L. H. Carvalho, E. M. M. Rocha, D. S. Raslan, A. B. Oliveira, and A. V. Krettli, “Natural and synthetic naphthoquinones against erythrocytic stages of Plasmodium falciparum,” Brazilian Journal of Medical and Biological Research, vol. 21, pp. 485–487, 1988.
View at: Google Scholar
J. S. Driscoll, G. F. Hazard, H. B. Wood, and A. Goldin, “Structure-antitumor activity relationships among quinone derivatives,” Cancer Chemotherapy, vol. 4, pp. 1–362, 1974.
View at: Google Scholar
J. N. Lopes, F. S. Cruz, R. Docampo et al., “In vitro and in vivo evaluation of the toxicity of 1,4-naphthoquinone and 1,2-naphthoquinone derivatives against Trypanosoma cruzi,” Annals of Tropical Medicine & Parasitology, vol. 72, no. 6, pp. 523–531, 1978.
View at: Publisher Site | Google Scholar
S. M. Wuerzberger, J. J. Pink, S. M. Planchon, K. L. Byers, W. G. Bornmann, and D. A. Boothman, “Induction of apoptosis in MCF-7:WS8 breast cancer cells by/3-Lapachone1,” Cancer Research, vol. 58, pp. 1876–1885, 1998.
View at: Google Scholar
S. Subramanian, M. M. C. Ferreira, and M. Trsic, “A structure-activity relationship study of lapachol and some derivatives of 1,4-naphthoquinones against carcinosarcoma walker 256,” Structural Chemistry, vol. 9, no. 1, pp. 47–57, 1998.
View at: Google Scholar
K. V. Rao, T. J. McBride, and J. J. Oleson, “Recognition and evaluation of lapachol as an antitumor agent,” Cancer Resources, vol. 28, pp. 1952–1954, 1968.
View at: Google Scholar
M. N. d. Silva, V. F. Ferreira, and M. C. B. V. d. Souza, “Um panorama atual da química e da farmacologia de naftoquinonas, com ênfase na beta-lapachona e derivados,” Química Nova, vol. 26, no. 3, pp. 407–416, 2003.
View at: Publisher Site | Google Scholar
T. S. Sadananda, R. Nirupama, K. Chaithra, M. Govindappa, C. P. Chandrappa, and B. Vinay, “Antimicrobial and antioxidant activities of endophytes from Tabebuia argentea and identification of anticancer agent (lapachol),” Journal of Medicinal Plants Research, vol. 5, no. 16, pp. 3643–3652, 2011.
View at: Google Scholar
K. H. Thompson, C. A. Barta, and C. Orvig, “Metal complexes of maltol and close analogues in medicinal inorganic chemistry,” Journal of the Chemical Society, vol. 35, p. 545, 2006.
View at: Google Scholar
T. A. Annan, C. Peppe, and D. G. Tuck, “The direct electrochemical synthesis of d10 metal ion derivatives of some anionic bidentate oxygen donors,” Canadian Journal of Chemistry, vol. 68, no. 3, pp. 423–430, 1990.
View at: Publisher Site | Google Scholar
M. I. F. Barbosa, R. S. Corrêa, K. M. de Oliveira et al., “Antiparasitic activities of novel ruthenium/lapachol complexes,” Journal of Inorganic Biochemistry, vol. 136, pp. 33–39, 2014.
View at: Publisher Site | Google Scholar
L. Tabrizi, F. Talaie, and H. Chiniforoshan, “Copper (II), cobalt (II) and nickel (II) complexes of lapachol: synthesis, DNA interaction, and cytotoxicity,” Journal of Biomolecular Structure and Dynamics, vol. 35, no. 15, 2017.
View at: Publisher Site | Google Scholar
R. Hernández-Molina, I. Kalinina, P. Esparza et al., “Complexes of Co(II), Ni(II) and Cu(II) with lapachol,” Polyhedron, vol. 26, no. 17, pp. 4860–4864, 2007.
View at: Publisher Site | Google Scholar
L. G. d. Oliveira, M. M. Silva, F. C. S. d. Paula et al., “Antimony(V) and bismuth(V) complexes of lapachol: synthesis, crystal structure and cytotoxic activity,” Molecules, vol. 16, no. 12, pp. 10314–10323, 2011.
View at: Publisher Site | Google Scholar
D. B. Mowrey, Muira–Puama (Liriosma Ovata) Herbal Tonic Therapies, McGraw-Hill Education, New York, NY, USA, 1996.
M. J. Zobo, M. D. Bikélé, D. Lissouck et al., “Thermodynamics-antioxidant activity relationships of some 4-benzylidenamino-4,5-dihydro-1h-1,2,4-triazol-5-one derivatives: theoretical evaluation,” International Journal of Food Properties, vol. 20, no. 9, pp. 1935–1948, 2017.
View at: Publisher Site | Google Scholar
J. S. Wright, E. R. Johnson, and G. A. DiLabio, “Predicting the activity of phenolic antioxidants: theoretical method, analysis of substituent effects, and application to major families of antioxidants,” Journal of the American Chemical Society, vol. 123, no. 6, pp. 1173–1183, 2001.
View at: Publisher Site | Google Scholar
A. J. Javan and M. J. Javan, “Electronic structure of some thymol derivatives correlated with the radical scavenging activity: theoretical study,” Food Chemistry, vol. 165, pp. 451–459, 2014.
View at: Publisher Site | Google Scholar
A. J. Javan, M. J. Javan, and Z. A. Tehrani, “Theoretical investigation on antioxidant activity of bromophenols from the marine red alga rhodomela confervoides: H-atom vs. electron transfer mechanism,” Journal of Agricultural and Food Chemistry, vol. 61, no. 7, pp. 1534–1541, 2013.
View at: Publisher Site | Google Scholar
J. Lengyel, J. Rimarčík, A. Vagánek, and E. Klein, “On the radical scavenging activity of isoflavones: thermodynamics of O-H bond cleavage,” Physical Chemistry Chemical Physics, vol. 15, no. 26, pp. 10895–10903, 2013.
View at: Publisher Site | Google Scholar
V. S. Jovanovic, S. Steenken, G. M. Simic, and Y. Hara, Flavonoids in Health and Disease, Marcel Dekker, New York, NY, USA, 1998.
C. P. V. Mary, S. Vijayakumar, and R. Shankar, “Metal chelating ability and antioxidant properties of curcumin-metal complexes—a DFT approach,” Journal of Molecular Graphics and Modelling, vol. 79, pp. 1–14, 2018.
View at: Publisher Site | Google Scholar
M. J. Frisch, G. W. Trucks, H. B. Schlegel et al., Gaussian 09, Revision A.02, Gaussian, Inc, Wallingford, UK, 2009.
E. Cances and B. Mennucci, “New applications of integral equations methods for solvation continuum models: ionic solutions and liquid crystals,” Journal of Mathematical Chemistry, vol. 23, no. 3-4, pp. 309–326, 1998.
View at: Google Scholar
E. Cancès, B. Mennucci, and J. Tomasi, “A new integral equation formalism for the polarizable continuum model: theoretical background and applications to Isotropic and anisotropic dielectrics,” The Journal of Chemical Physics, vol. 107, no. 8, pp. 3032–3041, 1997.
View at: Publisher Site | Google Scholar
L. A. Kelly and M. A. J. Rodgers, “Photoreduction of methyl viologen mediated by tris(bipyridyl)ruthenium(II) in inert colloidal suspensions,” Journal of Physical Chemistry, vol. 98, no. 25, pp. 6386–6391, 1994.
View at: Publisher Site | Google Scholar
Z. Markovic, J. Dorovic, J. M. Dimitric Markovic, R. Biocanin, and D. Amic, “Comparative density functional study of antioxidative activity of the hydroxybenzoic acids and their anions,” Turkish Journal of Chemistry, vol. 40, no. 3, pp. 499–509, 2016.
View at: Publisher Site | Google Scholar
M. D. Tissandier, K. A. Cowen, W. Y. Feng et al., “The proton’s absolute aqueous enthalpy and Gibbs free energy of solvation from cluster-ion solvation data,” The Journal of Physical Chemistry A, vol. 102, no. 40, pp. 7787–7794, 1998.
View at: Publisher Site | Google Scholar
M. Piris, J. M. Matxain, X. Lopez, and J. M. Ugalde, “The extended Koopmans’ theorem: vertical ionization potentials from natural orbital functional theory,” The Journal of Chemical Physics, vol. 136, no. 17, Article ID 174116, 2012.
View at: Publisher Site | Google Scholar
T. Lu and F. Chen, “Multiwfn: a multifunctional wavefunction analyzer,” Journal of Computational Chemistry, vol. 33, no. 5, pp. 580–592, 2012.
View at: Publisher Site | Google Scholar
S. Jenkins, Z. Liu, and S. R. Kirk, “A bond, ring and cage resolved Poincaré-Holf relationship for isomerization pathways,” Molecular Physics: An International Journal at the Interface between Chemistry and Physics, vol. 111, no. 20, pp. 3104–3116, 2013.
View at: Publisher Site | Google Scholar
E. Espinosa, E. Molins, and C. Lecomte, “Hydrogen bond strengths revealed by topological analyses of experimentally observed electron densities,” Chemical Physics Letters, vol. 285, no. 3-4, pp. 170–173, 1998.
View at: Publisher Site | Google Scholar
T. Steiner, “The hydrogen bond in the solid state,” Angewandte Chemie International Edition, vol. 41, no. 1, pp. 48–76, 2002.
View at: Publisher Site | Google Scholar
L. Mammino, “Investigation of the antioxidant properties of hyperjovinol a through its Cu(II) coordination ability,” Journal of Molecular Modeling, vol. 19, no. 5, pp. 2127–2142, 2013.
View at: Publisher Site | Google Scholar
M. K. Mwadham, T. T. Van, S. M. Kamogelo et al., “Conformational, electronic and antioxidant properties of lucidone, linderone and methyllinderone: DFT, QTAIM and NBO studies,” Molecular Physics, vol. 113, no. 7, pp. 683–697, 2015.
View at: Publisher Site | Google Scholar
H. J. Nomo, D. M. Bikélé, J. N. Ghogomu, and E. Younang, “A DFT study of structural and bonding properties of complexes obtained from first-row transition metal chelation by 3-alkyl-4-phenylacetylamino-4,5-dihydro-1H-1,2,4-triazol-5-one and its derivatives,” Bioinorganic Chemistry and Applications, vol. 2017, Article ID 5237865, 15 pages, 2017.
View at: Publisher Site | Google Scholar
T. Mihaylov, N. Trendafilova, I. Kostova, I. Georgieva, and G. Bauer, “DFT modeling and spectroscopic study of metal-ligand bonding in La(III) complex of coumarin-3-carboxylic acid,” Chemical Physics, vol. 327, no. 2-3, pp. 209–219, 2006.
View at: Publisher Site | Google Scholar
N. K. Nkungli and J. N. Ghogomu, “Theoretical analysis of the binding of iron(III) protoporphyrin IX to 4-methoxyacetophenone thiosemicarbazone via DFT-D3, MEP, QTAIM, NCI, ELF, and LOL studies,” Journal of Molecular Modeling, vol. 23, no. 7, pp. 1–20, 2017.
View at: Publisher Site | Google Scholar
X. Yunsheng, Z. Youguang, A. Lin, D. Yunyan, and L. Yi, “Density functional theory of the structure-antioxidant activity of polyphenolic deoxybenzoins,” Food Chemistry, vol. 151, pp. 198–206, 2014.
View at: Publisher Site | Google Scholar
X. Li and H. XieR. Zhan and D. Chen, “Effect of double bond position on 2-phenyl-benzofuran antioxidants: a comparative study of moracin C and iso-moracin C,” Molecules, vol. 23, no. 4, p. 754, 2018.
View at: Publisher Site | Google Scholar
L. Chen, “Dietary phenolic compound with the presence of C2=C3 double bond take the pre-emptive opportunities to enhance its biological effects,” Journal of Food Microbiology, vol. 2, no. 1, pp. 4–6, 2018.
View at: Google Scholar
M. Salga, I. Sada, and A. Mustapha, “Influence of steric hindrance on the antioxidant activity of some schiff base ligands and their copper(II) complexes,” Oriental Journal of Chemistry, vol. 30, no. 4, pp. 1529–1534, 2014.
View at: Publisher Site | Google Scholar
J. Sun, Y. Huang, G. Sun, X. Sun, M. Qin, and D. Zhao, “Study on in vitro antioxidant activity of flavonoids contained in Hebei balmy chrysanthemum and structure-activity relationship,” Zhongguo Zhong Yao Za Zhi, vol. 37, no. 13, pp. 1958–1962, 2012.
View at: Google Scholar
Y. Chen, H. Xiao, J. Zheng, and G. Liang, “Structure thermodynamics-antioxidant activity relationships of selected natural phenolic acids and derivatives: an experimental and theoretical evaluation,” PLoS One, vol. 10, no. 3, Article ID e0121276, 2015.
View at: Publisher Site | Google Scholar
A. D. T. Fouegue, D. Bikélé Mama, J. N. Ghogomu, Y. Elie, and M.-A. Etoh, “The substitution effect on reaction, enthalpies of antioxidant mechanisms of juglone and its derivatives in gas and solution phase: DFT study,” Journal of Chemistry, vol. 2018, Article ID 1958047, 10 pages, 2018.
View at: Publisher Site | Google Scholar
O. Holtomo, M. Nsangou, J. J. Fifen, and O. Motapon, “Antioxidative potency and UV-Vis spectra features of the compounds resulting from the chelation of Fe²⁺ by caffeic acid phenethyl ester and two of its derivatives,” Computational and Theoretical Chemistry, vol. 1067, pp. 135–147, 2015.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2020 Djafarou Ngouh Pajoudoro et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies