Next Article in Journal
In Vitro Antioxidant, Anti-Inflammatory Activity and Bioaccessibility of Ethanolic Extracts from Mexican Moringa oleifera Leaf
Next Article in Special Issue
Polyphenol-Rich Aronia melanocarpa Fruit Beneficially Impact Cholesterol, Glucose, and Serum and Gut Metabolites: A Randomized Clinical Trial
Previous Article in Journal
Isolation, Identification and Antibacterial Characteristics of Lacticaseibacillus rhamnosus YT
Previous Article in Special Issue
Vasorelaxant Effect and Blood Pressure Reduction Potential of Pitaya Juice Concentrate (Stenocereus huastecorum) Associated with Calcium Channel Blockade
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

QM-CSA: A Novel Quantum Mechanics-Based Protocol for Evaluation of the Carcinogen-Scavenging Activity of Polyphenolic Compounds

1
Faculty of Chemistry and Chemical Engineering, University of Maribor, Smetanova 17, SI-2000 Maribor, Slovenia
2
Institute of Environmental Protection and Sensors, Beloruska Ulica 7, SI-2000 Maribor, Slovenia
3
Faculty of Mathematics, Natural Sciences and Information Technologies, University of Primorska, Glagoljaška 8, SI-6000 Koper, Slovenia
*
Authors to whom correspondence should be addressed.
Foods 2024, 13(17), 2708; https://doi.org/10.3390/foods13172708
Submission received: 25 July 2024 / Revised: 20 August 2024 / Accepted: 22 August 2024 / Published: 27 August 2024
(This article belongs to the Special Issue Polyphenols and Health Benefits: 2nd Edition)

Abstract

:
In this study, a novel quantum mechanics-based protocol for the evaluation of carcinogen-scavenging activity (QM-CSA) is developed. The QM-CSA protocol represents a universal and quantitative approach to evaluate and compare the activation-free energies for alkylation reactions between individual polyphenolic compounds and chemical carcinogens of the epoxy type at physiological conditions by applying two scales: the absolute scale allowing for the comparison with guanine and the relative scale allowing the comparison with glutathione as a reference compound. The devised quantum mechanical methodology was validated by comparing the activation-free energies calculated with 14 DFT functionals in conjunction with two implicit solvation models (SMD and CPCM) and the experimental activation-free energies for reactions between nine investigated chemical carcinogens and guanine. According to the obtained results, the best agreement with experimental data was achieved by applying DFT functionals M11-L and MN12-L in conjunction with the flexible 6-311++G(d,p) basis set and implicit solvation model SMD, and the obtained uncertainties were proven to be similar to the experimental ones. To demonstrate the applicability of the QM-CSA protocol, functionals M11-L, and MN12-L in conjunction with the SMD implicit solvation model were applied to calculate activation-free energies for the reactions of nine investigated chemical carcinogens of the epoxy type with three catechins, namely EGCG, EGC, and (+)-catechin. The order of CSA in this series of catechins in comparison to guanine and glutathione was determined as (+)-catechin > EGC > EGCG. The obtained results, for the first time, demonstrated the evaluation and comparison of CSA in a series of selected catechins with respect to glutathione and guanine. Moreover, the presented results provide valuable insights into the reaction mechanisms and configurations of the corresponding transition states. The novel QM-CSA protocol is also expected to expand the kinetic data for alkylation reactions between various polyphenolic compounds and chemical carcinogens of the epoxy type, which is currently lacking in the scientific literature.

Graphical Abstract

1. Introduction

Cancer represents the second deadliest disease of the modern era after cardiovascular disorders. Carcinogenesis is usually triggered by DNA damage, which is, in most cases, caused by genotoxic chemical carcinogens that enter our body from the environment or are ingested with food. The term chemical carcinogens refers to compounds that contribute to cancer initiation by causing damage or mutations in DNA. If DNA damage results in oncogene activation or the deactivation of tumor suppressor genes, normal cells become neoplastic, which leads to cancer development [1]. The current study is focused on nine exogenous chemical carcinogens, which are implicated in the etiology of various cancer types. Aflatoxin B1 represents the most potent natural chemical carcinogen known to man, synthesized by molds of the Aspergillus genus, which grow on cereals, nuts, and oils in storage units with high temperatures and humidity. Aflatoxin B1 is, after ingestion, metabolized by cytochrome P450 3A enzymes to the highly reactive ultimate chemical carcinogen aflatoxin B1-exo-8,9-epoxide (AFB1), which can damage DNA and initiate cancer development [2]. Ethylene oxide (EO) is the ultimate chemical carcinogen of ethylene and is used for the sterilization of medical instruments in third-world countries [3]. Chloroethylene oxide (CEO) represents an electrophilic metabolite of vinyl chloride, which is used in the farming, plastic, and electrical industries [4]. Vinyl carbamate epoxide (VCE) is the ultimate chemical carcinogen of urethane, which is formed as a by-product of fermentation [5]. The ultimate chemical carcinogen glycidamide (GA) represents a metabolite of acrylamide, which develops while frying and baking various foods [6]. Styrene oxide (SO) is metabolized from styrene, which is often applied in the plastic industry [7]. Propylene oxide (PO) is the ultimate chemical carcinogen of propylene used as a medium for sterilization as well as for the production of polyurethane foams, solvents, and coolants [8]. 2-Cyanoethylene oxide (CyEO) is formed by the metabolic oxidation of acrylonitrile, which is applied to the manufacturing of acrylic fibers, nitrile rubber, and adhesives [9]. Finally, beta-propiolactone (BPL) represents a chemical carcinogen used in vaccines to prevent viral infections [10]. Biotransformations of nine procarcinogens into investigated genotoxic ultimate chemical carcinogens are depicted in Figure 1. Except for beta-propiolactone, the investigated chemical carcinogens are metabolized in the human body by cytochrome P450 enzymes into genotoxic ultimate chemical carcinogens of the epoxy type, which react with DNA primarily at the N7 atom of the most nucleophilic base guanine, leading to DNA damage and cancer initiation [2,4].
Glutathione represents a cellular tripeptide and an important hydrophilic biomolecule that participates in cellular detoxification processes. Glutathione regulates oxidative stress and protects cells from lipid peroxides, reactive oxygen and nitrogen species, and electrophilic xenobiotics [11]. Glutathione acts as the strongest natural non-enzymatic scavenger of chemical carcinogens in human cells by reducing their reactivity and subsequent toxicity. The conjugation of chemical carcinogens with glutathione facilitates their greater polarity and easier elimination from the body. According to recent studies, glutathione also represents the key regulator of cancer cell differentiation, proliferation, apoptosis, and immune function [12].
Polyphenols form a large group of more than 10,000 secondary plant metabolites with one or more aromatic rings containing one or more hydroxyl groups. Polyphenols are abundantly present in a majority of fruits, vegetables, and herbs and exert several beneficial health effects [13,14]. According to their chemical structure, polyphenols can be divided into five classes, namely flavonoids, phenolic acids, stilbenes, lignans, and tannins, among which flavonoids and polyphenolic acids are the most abundant [15]. Catechins (flavan-3-ols) represent one of the most important subclasses of flavonoids, gaining recognition as chemopreventive agents due to their various beneficial health effects as well as their safety [16]. Green tea (Camellia sinensis) represents one of the most consumed beverages, which is a rich source of epistructured gallocatechins, namely epigallocatechin-3-gallate (EGCG) and epigallocatechin (EGC), with proven health benefits [17,18,19]. On the other hand, cocoa represents a rich source of nonepistructured catechins, especially (+)-catechin, exerting various beneficial biological effects [20]. Several in vitro and in vivo studies reported that green tea and cocoa consumption exhibits antioxidative and anti-inflammatory effects, lowers the risk of hypertension, prevents cardiovascular disease and type II diabetes, aa well as improves brain function, although catechins are mainly investigated for their cancer chemopreventive and chemotherapeutic effects [21,22,23,24,25].
EGCG and EGC represent important catechins obtained from green tea with various known beneficial health effects, including antioxidant [26,27], anti-inflammatory [28], antimicrobial [29], and anticarcinogenic [30] effects. EGCG and EGC exhibited inhibitory activities on cyclooxygenase 2 (COX-2), leading to suppressed cell proliferation, tumor invasion, and anti-inflammatory effects [31,32]. The strong anticarcinogenic character of the green tea catechins EGCG and ECG was demonstrated through the reduction of carcinogen-induced reactive oxygen species (ROS), antiproliferative, and pro-apoptotic effects in MCF7 breast cancer cells [33]. The most prevalent green tea catechin, EGCG, also exhibited antiangiogenic and antimetastatic effects as well as suppressed proliferation and invasion of cancer cells through down-regulation of matrix metalloproteinase (MMP) enzymes MMP-2 [34] and MMP-9 [35]. EGCG also inhibited the activation of hypoxia-inducible factor-1 (HIF-1α), nuclear factor kappa B (NF-κB), and vascular endothelial growth factor (VEGF) expression, resulting in the suppression of angiogenesis and breast cancer progression [36]. Moreover, EGCG also demonstrated inhibition of B-cell lymphoma 2 (Bcl-2) and NF-κB in hepatocellular cancer cells LM6 [37], and deactivated activator protein-1 (AP-1) expression, resulting in apoptosis of AGS gastric cancer cells [38].
Anticarcinogenic effects of important catechin from cocoa (+)-catechin have also been extensively investigated [39,40,41,42]. (+)-Catechin exhibited anticancer effects against A549 lung cancer cells through the dose-dependent induction of cyclin kinase inhibitor p21 and suppression of cyclin E1, leading to cell cycle arrest [42]. It was also reported that (+)-catechin up-regulated p53 in colorectal cancer cells HT-29, promoting cell cycle arrest and inhibition of cell proliferation. Moreover, (+)-catechin demonstrated cytotoxicity in colorectal cancer cells HT-29 and activated their apoptosis through up-regulation of caspases 3 and 8 [41]. (+)-Catechin can also induce the inhibition of angiogenesis by regulating nitric oxide (NO) and tumor necrosis factor alpha (TNF-α) production [39]. Furthermore, (+)-catechin reduced VEGF mRNA levels in B16F-10 melanoma cells. The anticancer effects of (+)-catechin against human colon adenocarcinoma cells HCT 15 and HCT 116 as well as against human larynx carcinoma HepG-2 cells were also investigated [40]. The obtained results showed that (+)-catechin treatment resulted in decreased cell proliferation and increased the presence of apoptotic bodies in the investigated cancer cells.
Although the anticarcinogenic effects of catechins EGCG, EGC, and (+)-catechin (Figure 2) were reported in several experimental studies, their potential as direct scavengers of investigated chemical carcinogens remains unexplored.
Nowadays, it is possible to perform quantum mechanical calculations and obtain reliable kinetic data for small- or medium-sized systems up to 100 atoms at a reasonable computational cost. Quantum mechanics-based methodologies currently play an important role in investigations of reaction mechanisms and represent promising alternatives to experimental approaches, particularly in the case of difficult reaction conditions and expensive or (geno)toxic chemicals. In our previous studies, we reported that the polyphenols 6-gingerol and xanthohumol and its derivatives isoxhanthohumol, 8-prenylnaringenin, and 6-prenylnaringenin represent potential scavengers of ultimate chemical carcinogens of the epoxy type [43,44]. The obtained results at the Hertree–Fock/6-311++G(d,p) level of theory, in combination with two implicit solvation models, confirmed the proposed SN2 reaction mechanism of the most abundant monoanionic forms of the studied polyphenols at physiological conditions. Taking advantage of quantum mechanical methodologies to systematically investigate the carcinogen-scavenging activity (CSA) of polyphenolic compounds, therefore, represents a safe, fast, and inexpensive alternative to experimental assays.
It must be noted that plant extracts contain a complex composition of various polyphenolic compounds, which results in challenging chemical analyses and exhaustive evaluations of compounds or combinations of compounds that are responsible for the specific biological activity [45]. By using quantum mechanics-based computational approaches, the potential of individual polyphenolic compounds, such as EGCG, EGC, and (+)-catechin, to scavenge chemical carcinogens can be quantitatively evaluated and compared. The main purpose of this comparison is to identify polyphenolic compounds with the highest CSA and hopefully guide and accelerate the design of efficient strategies to prevent cancer initiation. However, a universal and quantitative quantum mechanical methodology to evaluate and compare the carcinogen-scavenging activity of individual polyphenolic compounds by applying an appropriate combination of Density Functional Theory (DFT) functionals and implicit solvation models is still not established. The applicability of theoretical quantum mechanics-based approaches to study chemical reactions is mainly influenced by the selected method, flexible basis set, and implicit solvation model, which determine the accuracy of the obtained results. Moreover, the influence of the pH on the protonation state of polyphenolic compounds should also be taken into account. pKa values of hydroxyl groups of polyphenolic compounds determine the portion of neutral, mono-anionic, and dianionic forms of polyphenolic compounds at a specific pH value [46].
In this study, a novel quantum mechanics-based protocol for the evaluation of carcinogen-scavenging activity (QM-CSA) by applying the advanced Density Functional Theory (DFT) functionals M11-L and MN12-L in combination with the flexible 6-311++G(d,p) basis set and implicit solvation model density (SMD) is presented. The developed QM-CSA protocol allows for the evaluation of the molecular mechanisms and kinetics of the alkylation reactions involved in the CSA of polyphenolic compounds by applying two scales: the absolute scale yielding the comparison with guanine and the relative scale yielding the comparison with glutathione as a reference compound. To demonstrate the applicability of the devised QM-CSA protocol, the CSA of the important epicatechins EGCG and EGC, as well as that of nonepistructured (+)-catechin, was evaluated. The selected catechins have not been previously studied in the role of scavengers of chemical carcinogens; therefore, their mechanistic and kinetic data are provided for the first time.

2. Computational Methods

The aim of this study was to develop a standardized quantum mechanics-based methodology for the evaluation of carcinogen-scavenging activity (QM-CSA). The QM-CSA protocol is intended to establish a universal strategy to evaluate and compare the primary CSA of natural polyphenols by applying an appropriate combination of DFT functional, flexible basis set, and implicit solvation model. The QM-CSA protocol is based on kinetic evaluations and considers the biological conditions (e.g., pH or solvent) that may influence the CSA.
The QM-CSA protocol aims to provide:
  • The identification of the main molecular mechanism of alkylation reactions between nine chemical carcinogens and polyphenolic compounds at physiological conditions (pH 7.4 and an aqueous environment).
  • Two scales for the quantification of CSA: the absolute scale based on the comparison of the calculated activation-free energies for alkylation reactions of polyphenols with guanine and the relative scale using glutathione as a reference compound.
  • The establishment of the order of CSA in a series of polyphenolic compounds at physiological conditions.
  • The identification of the most effective scavengers of chemical carcinogens of the epoxy type and the structural features responsible for their activity.
The QM-CSA protocol has been devised to systematically address the important aspects, which are presented in the following subsections.

2.1. Investigation of Molecular Mechanisms of Alkylation Reactions of Guanine, Glutathione, and Polyphenolic Compounds with Chemical Carcinogens at Physiological Conditions

All nine studied chemical carcinogens of the epoxy type represent electrophilic species that primarily react with the most nucleophilic DNA base guanine at the atom N7, resulting in the covalent chemical carcinogen–guanine adduct formation. The rate-limiting step for the reactions of chemical carcinogens with nucleophilic DNA bases is the epoxy ring opening [6]. Using the ab initio Hertree–Fock method in combination with the flexible 6-311++G(d,p) basis set and self-consistent reaction field implicit solvation models, our research group successfully reproduced experimentally determined activation-free energies of alkylation reactions between the most reactive DNA base guanine and nine genotoxic chemical carcinogens, namely aflatoxin B1 exo-8,9-epoxide [2], 2-cyanoethylene oxide [9], glycidamide [6], chloroethylene oxide [4], ethylene oxide [3], styrene oxide [7], propylene oxide [8], beta-propiolactone [10], and vinyl carbamate epoxide [5], which corresponded to the proposed SN2 reaction mechanism. Subsequent protonation is a fast process because of the proton-rich microenvironment surrounding the negatively charged DNA. The depurination of formed adducts then leads to mutations and chromosomal aberrations [6]. The proposed SN2 substitution mechanisms for the reactions of nine investigated chemical carcinogens with (a) guanine and (b) glutathione are presented in Scheme A1, Scheme A2, Scheme A3, Scheme A4, Scheme A5, Scheme A6, Scheme A7, Scheme A8 and Scheme A9 in Appendix A, which were prepared in the ChemDraw program. For the rate-limiting step of alkylation reactions between polyphenolic compounds and the nine investigated chemical carcinogens at physiological conditions (pH 7.4 and an aqueous solution), the SN2 reaction mechanism is proposed as well.

2.2. The Determination of an Appropriate Combination of Density Functional, Flexible Basis Set, and Implicit Solvation Model

In order to determine the best combination of quantum mechanical method, flexible basis set, and implicit solvation model to calculate activation-free energies of alkylation reactions between the nine chemical carcinogens and guanine, 14 advanced DFT functionals in combination with the 6-311++G(d,p) basis set and two implicit solvation models, namely the Conductor-like Polarizable Continuum Model (CPCM) and the Solvation Model Density (SMD), were applied. The selected set of DFT functionals included 10 frequently applied meta-Generalized Gradient Approximations (GGA) exchange–correlation functionals, namely M05-2X [47], M06-2X [48], M06-HF [49], M08-HX [50], M11 [51], M11-L [52], MN12-L [53], MN12-SX [54], MN15 [55], and MN15-L [56], 2 meta-Non-separable Generalized Gradient Approximations (NGA) exchange–correlation functionals N12 [57] and N12-SX [54], as well as 2 hybrid-GGA exchange–correlation functionals B3LYP-D3 [58] and WB97XD [59]. All quantum mechanical calculations of activation-free energies were performed with the program package Gaussian 16 on the supercomputer system HPC Vega located at the Institute of Information Science (IZUM). The computationally obtained activation-free energies for reactions between guanine and investigated chemical carcinogens were then compared to their corresponding experimental counterparts. The latter were determined from the experimental reaction rate constants k [2,3,4,5,6,7,8,9,10], which can be obtained by applying the Transition State Theory (TST) of Eyring (Equation (1)) [60]:
k = k B T h e ( Δ G k B T )
where k B and h are the Boltzmann and Planck constants, respectively; T is the temperature (K); and Δ G is the activation-free energy, representing the free energy difference between the transition states and reactants (kcal/mol). The main advantage of using TST is that it requires only the free energy of reactants and transition state structures for the calculation of the corresponding Δ G , which allows for its application to various chemical reactions. TST is applicable for medium- to large-size systems, such as chemical carcinogens of the epoxy type reacting with guanine. To evaluate the accuracy of the devised QM-CSA methodology, the calculated activation-free energies for the reactions of the investigated chemical carcinogens with guanine were compared to their experimental counterparts.
Based on the analysis of the obtained results (Section 3), the combination of meta-GGA functional M11-L [52] or MN12-L [53] with flexible basis set 6–311++G(d,p) and SMD implicit solvation model provided the best agreement with experimental data. Meta-GGA exchange-correlation functional M11-L is based on a dual-range local exchange to ensure accuracy for single-configurational as well as multiconfigurational molecules, solid-state lattice constants, and barrier heights [52]. A local exchange-correlation meta-GGA functional MN12-L was designed to study transition metal chemistry, atomization energies, ionization potentials, noncovalent interactions, solid-state lattice constants, and barrier heights [53]. Both M11-L and MN12-L functionals have also been recommended for calculations of activation-free energies by their developers [52,53] and successfully applied by independent authors to predict reaction barrier heights [61,62,63].
In terms of the choice of a flexible basis set, it is crucial to include polarization and diffuse functions for nonhydrogen atoms. The application of diffuse functions is important to investigate polyphenolic compounds in anionic forms, which play a key role in CSA. The Gaussian basis set 6-311++G(d,p) represents the flexible basis set of choice in terms of more satisfactory results and lower computational cost when compared with the augmented correlation consistent polarized valence double zeta (aug-cc-pVDZ) basis set [64].
As the presence of the solvent significantly affects the reactivity of molecules, it is important to include solvation effects in quantum mechanical calculations of alkylation reactions between chemical carcinogens of the epoxy type and polyphenolic compounds. Different solvation models and their advantages were comprehensively addressed in review articles [65,66,67]. In this study, the most widely applied implicit solvation models, namely SMD and CPCM, were included in quantum mechanical calculations of activation-free energies to mimic the aqueous environment surrounding chemical carcinogens and guanine in human cells. The combination of DFT functionals M11-L or MN12-L with flexible 6-311++G(d,p) basis sets and the SMD implicit solvation model proved to be the most successful at reproducing the experimental activation-free energies. The implicit solvation model SMD employs the quantum mechanical charge density of a solute that interacts with a polarizable continuum of the solvent. SMD represents a universal implicit solvation model that can be applied to any charged or uncharged solute in various polar and organic solvents [67]. To conclude, DFT functionals M11-L and MN12-L in combination with the flexible 6-311++G(d,p) basis set and SMD solvation model are recommended to study alkylation reactions with chemical carcinogens of the epoxy type.

2.3. Determination of Anticarcinogenic Potential of Various Polyphenolic Compounds as Scavengers of Nine Investigated Chemical Carcinogens of the Epoxy Type Using QM-CSA Methodology

For comparison and evaluation of the carcinogen-scavenging activity of polyphenolic compounds, it is recommended to apply the same quantum mechanical methodology, namely DFT functional M11-L or MN12-L, in conjunction with the flexible 6-311++G(d,p) basis set and SMD implicit solvation model. It should also be noted that the reliability of the obtained results depends on the chosen functional, and that only functionals M11-L and MN12-L in combination with the SMD solvation model, proved to accurately reproduce the experimental activation-free energies for alkylation reactions of the nine studied chemical carcinogens with guanine. Moreover, the studied alkylation reactions should take place in an aqueous solution at a physiological pH in order to rightfully compare the calculated activation-free energies. By applying the same computational methodology, a comprehensive comparison of the CSA of different polyphenolic compounds can be performed by investigating from one to all nine chemical carcinogens, namely aflatoxin B1-exo-8,9-epoxide, ethylene oxide, chloroethylene oxide, vinyl carbamate epoxide, glycidamide, styrene oxide, propylene oxide, 2-cyanoethylene oxide, and beta-propiolactone.
In this study, the CSA of three investigated catechins, EGCG, EGC, and (+)-catechin, as well as glutathione under physiological conditions (pH = 7.4, aqueous solution), was determined by applying DFT functionals M11-L and MN12-L in combination with the flexible 6-311++G(d,p) basis set and SMD solvation model. EGCG, EGC, and (+)-catechin were studied in their nucleophilic (monoanionic) forms at physiological conditions. The pKa values of the most acidic protons of EGCG, EGC, and (+)-catechin at the pH of 7.4 were calculated with MarvinSketch software (20.19, ChemAxon, Budapest, Hungary) [68], and their values were in agreement with the experimentally determined ones [46]. Geometry optimization and vibrational analysis of reactant structures were performed to obtain the corresponding minima on the potential energy surface. To obtain the initial geometries of corresponding transition states, a relaxed potential surface scan was performed [35]. The identified initial geometries were subsequently subjected to geometry optimization and vibrational analysis [37]. Quantum mechanical calculations of activation-free energies were performed with Gaussian 16 on the supercomputer HPC Vega located at IZUM. Only the first step of the reaction between the individual chemical carcinogen and the investigated catechins, as well as guanine and glutathione, was examined, as subsequent protonation is believed to be a fast process due to a proton-rich microenvironment surrounding the product of the initial reaction step [43].

2.4. Evaluation of the CSA of Polyphenolic Compounds by Applying Two Scales

A comparison of different polyphenolic compounds in the role of scavengers of nine chemical carcinogens can be established by applying two scales.
The first scale allows for an absolute comparison between activation-free energies Δ G for alkylation reactions of the investigated polyphenol and the most nucleophilic DNA base, guanine, with the same chemical carcinogens. For this purpose, the calculated Δ G value of the reaction between polyphenol and the chemical carcinogen is compared with the Δ G value of the corresponding competing reaction between guanine and the same chemical carcinogen. If the Δ G value for the reaction of the chemical carcinogen with polyphenol is lower than the Δ G value of the competing reaction of the same chemical carcinogen with guanine, the investigated polyphenol represents an efficient scavenger of this chemical carcinogen (Figure 3). The lower Δ G value directly corresponds to a higher reaction rate constant and indicates a higher CSA of the polyphenolic compound [43]. The absolute scale is, therefore, defined by Equation (2):
Δ Δ G = Δ G p o l y p h e n o l Δ G g u a n i n e
In the second scale, glutathione, the strongest natural scavenger of chemical carcinogens in human cells, is applied as a reference compound [69,70]. Alkylation reactions between glutathione and the same nine chemical carcinogens were investigated under identical conditions (pH value of 7.4, aqueous solution) as the reactions with polyphenols and guanine by applying the DFT functionals M11-L and MN12-L in combination with the flexible 6-311++G(d,p) basis set and SMD solvation model. The greater or lower activation-free energies of the investigated polyphenols compared to glutathione determine their relative CSA. The formula to calculate the factor r g l u t a t h i o n e is provided in Equation (3):
r g l u t a t h i o n e = Δ G p o l y p h e n o l Δ G g l u t a t h i o n e
where Δ G p o l y p h e n o l represents the activation-free energy for the reaction of polyphenol with the chemical carcinogen, and Δ G g l u t a t h i o n e represents the activation-free energy for the competing reaction of glutathione with the same chemical carcinogen. Therefore, the relative CSA of polyphenolic compounds can be expressed as r-times lower or higher than that of glutathione ( r g l u t a t h i o n e = 1). In this way, the most effective polyphenolic scavengers of chemical carcinogens can be identified based on the comparison with guanine and glutathione, and a ranked series of activity can also be established. This approach is also expected to maximize the cancelation of errors inherent to all calculated activation-free energies.

3. Results and Discussion

3.1. Validation of the Developed QM-CSA Protocol

Due to the absence of experimental data on the activation-free energies of alkylation reactions between the nine studied chemical carcinogens and polyphenolic compounds, a test set of nine experimental activation-free energies for the competing reactions between the nine chemical carcinogens and guanine was used to validate the devised QM-CSA protocol. In Figure 4, the differences between the experimental activation-free energies and the ones calculated with 14 most frequently applied DFT functionals, namely M05-2X, M06-2X, M06-HF, M08-HX, M11, M11-L, MN12-L, MN12-SX, MN15, MN15-L, N12, N12-SX, B3LYP-D3, and wB97XD, in combination with the flexible 6-311++G(d,p) basis set and two implicit solvation models, namely SMD and CPCM, are collected. The dark green color designates the best agreement with the corresponding experimental activation-free energy for the alkylation reaction between individual chemical carcinogens and guanine, while the light green color designates a poor agreement with the corresponding experimental value.
As can be observed from Figure 4, a large majority of studied functionals accurately reproduced the experimental activation-free energy for the reactions between guanine and BPL with a minimum Δ(ΔG) difference obtained by the combination N12-SX/SMD (−0.04 kcal/mol) and a maximum Δ(ΔG) difference in the case of MN15-L/CPCM (3.96 kcal/mol). On the other hand, a large majority of the applied functionals very poorly reproduced the experimental activation-free energy in the case of AFB1, with a minimum Δ(ΔG) difference obtained with MN12-L/SMD (1.78 kcal/mol) and a maximum Δ(ΔG) difference in the case of MN15/SMD (6.98 kcal/mol). Interestingly, the most frequently applied hybrid meta-GGA functional M06-2X in combination with solvation models SMD or CPCM poorly reproduced the experimental activation-free energies for reactions between eight investigated chemical carcinogens and guanine, with Δ(ΔG) differences ranging from 2.80 kcal/mol in the case of BPL (M06-2X/SMD) to 6.88 kcal/mol in the case of EO (M06-2X/CPCM), while VCE posed the only exception with Δ(ΔG) below 1 kcal/mol for combination with either solvation model. From Figure 4, it is also evident that the activation-free energies calculated using M11-L and MN12-L DFT functionals in conjunction with the 6-311++G(d,p) basis set and SMD solvation model provide the best overall agreement with the experimental data for all nine investigated chemical carcinogens. The obtained Δ(ΔG) values range from 0.31 kcal/mol (BPL) to 1.89 kcal/mol (AFB1) in the case of M11-L/SMD and from 0.05 kcal/mol (PO) to −1.86 kcal/mol (CEO) in the case of MN12-L/SMD.
In Table 1, the calculated activation-free energies with M11-L/SMD and MN12-L/SMD methods, the lowest vibrational frequencies and distances between the reactive centers corresponding to reactant and transition states, as well as absolute and relative errors, together with the corresponding experimental activation-free energies for reactions of guanine with nine investigated chemical carcinogens, are collected. The calculated activation-free energies with the corresponding absolute and relative errors as well as vibrational frequencies of reactants and transition states for all tested DFT functionals are provided in Supplementary Material Tables S1–S9.
From Table 1, it can be observed that the lowest obtained real frequencies of optimized reactant structures of guanine and nine investigated chemical carcinogens were between 8.28 cm−1 and 45.33 cm−1. The obtained single imaginary frequencies of transition states ranging from 207.05 i cm−1 to 654.07 i cm−1 corresponded to the cleavage of the covalent bond inside the epoxy ring of chemical carcinogen and to the formation of a new covalent bond between the most nucleophilic atom N7 of guanine and the achiral epoxy carbon of the nine chemical carcinogens (Table 1). The SN2 substitution mechanisms for the formation of the nine chemical carcinogen–guanine adducts were proposed (Appendix A). The distances between the reactive centers of the reactants dR lied between 2.79 Å and 4.24 Å and were approximately 1 Å to 2 Å shorter than the corresponding distances in transition states ranging from 1.94 Å to 2.19 Å. The shorter distances correspond to stronger intermolecular interactions in the transition states than in the reactant states. Similar distances between reactive centers in transition state structures dTS obtained with the M11-L/SMD and MN12-L/SMD methods also indicate similar geometries of optimized transition state structures. The predicted activation-free energies by applying the M11-L and MN12-L DFT functionals in combination with the flexible 6-311++G(d,p) basis set and SMD solvation model were in the range of 16.99 kcal/mol to 27.16 kcal/mol and from 16.88 kcal/mol to 26.49 kcal/mol, respectively. Among all chemical carcinogens, the Δ G s m d obtained with the M11-L/SMD and MN12-L/SMD methods was the lowest for the reaction of AFB1 with guanine (16.99 kcal/mol and 16.88 kcal/mol, respectively), indicating that the reaction of guanine with AFB1 is the fastest. The Δ G s m d obtained with M11-L/SMD and MN12-L/SMD is the highest for the reaction of guanine with styrene oxide (27.16 kcal/mol and 25.50 kcal/mol, respectively), indicating the slowest reaction of guanine with this chemical carcinogen. From Table 1, it is also evident that the DFT functionals M11-L/SMD and MN12-L/SMD predicted very similar values of activation-free energies for reactions of all nine investigated chemical carcinogens with guanine (difference ~1 kcal/mol).
The largest discrepancy between the calculated and experimental activation-free energy was found for the reaction of AFB1 exo-8.9-epoxide and of vinyl carbamate epoxide with guanine in the case of M11-L/SMD, as well as for the reaction of chloroethylene oxide and AFB1 exo-8.9-epoxide in the case of MN12-L/SMD. The corresponding absolute errors | Δ ( Δ G ) | were 1.89 kcal/mol and 1.75 kcal/mol for M11-L/SMD, respectively, as well as 1.86 kcal/mol and 1.78 kcal/mol for MN12-L/SMD, respectively. In addition, the corresponding relative errors obtained with M11-L/SMD for guanine reactions with AFB1 exo-8.9-epoxide and vinyl carbamate were 0.12 and 0.08, respectively, as well as 0.09 and 0.12, respectively, for chloroethylene oxide and AFB1 exo-8.9-epoxide reactions with guanine obtained using MN12-L/SMD. However, the mean absolute error for the test set of nine investigated chemical carcinogens reacting with guanine was found to be 1.08 kcal/mol for M11-L/SMD and 1.11 kcal/mol for MN12-L/SMD, which corresponds to an acceptable computational accuracy of 1 kcal/mol. The mean absolute and relative errors for all tested functionals are provided in Supplementary Material Table S10. As can be observed from Figure 4 and Table S10, the DFT functionals N12SX/SMD, M06HF/SMD, and WB97XD/SMD may also provide reliable kinetic data for some chemical carcinogens (their mean absolute errors were 2.28, 2.39, and 2.41 kcal/mol, respectively); however, their application is not recommended due to the very poor agreement with experimental values in the case of AFB1 exo-8,9-epoxide (Δ(ΔG) values were 6.38, 6.51, and 5.22 kcal/mol, respectively).
The correlation between the calculated activation-free energies using (a) M11-L/SMD and (b) MN12-L/SMD and the experimental activation-free energies for the reactions between the nine investigated chemical carcinogens and guanine is presented in Figure 5. The obtained R2 values for M11-L/SMD and MN12-L/SMD were 0.88 and 0.87, respectively, indicating a high correlation between the calculated and experimental activation-free energies for reactions between guanine and the nine investigated chemical carcinogens.
The obtained low mean absolute and relative errors as well as high R2 values demonstrated the reliability of M11-L and MN12-L DFT functionals in combination with the flexible 6-311++G(d,p) basis set and SMD implicit solvation model to evaluate the reactivity of the nine investigated chemical carcinogens with guanine. The QM-CSA methodology is valid for all chemical carcinogens of the epoxy type and can reliably estimate the CSA. Moreover, the obtained activation-free energies for reactions between nine investigated chemical carcinogens and guanine will allow for the evaluation of the CSA of various natural polyphenolic compounds. The optimized transition state structures obtained with the M11-L/SMD method, corresponding to the studied alkylation reactions between the nine investigated chemical carcinogens and guanine, are presented in Figure 6.
Graphical representations of optimized transition state geometries revealed the formation of a new covalent bond between the N7 atom of guanine and the achiral epoxy carbon atom of the nine investigated chemical carcinogens. The simultaneous cleavage of a covalent bond between the achiral carbon and epoxy oxygen of the chemical carcinogen also confirms the correct transition state structure allocation. The obtained single imaginary frequencies with normal modes following this reaction path for all transition state structures (Table 1) correspond to the proposed SN2 reaction mechanism, which is also in agreement with the previously published scientific literature [2,3,4,5,6,7,8,9,10].

3.2. Application of the QM-CSA Protocol: Mechanistic Insights into Alkylation Reactions of Nine Investigated Chemical Carcinogens with EGCG, EGC, and (+)-Catechin

In this subsection, the application of the QM-CSA protocol to the alkylation reactions between the nine investigated chemical carcinogens and three catechins, namely EGCG, EGC, and (+)-catechin, is presented. The previously discussed results revealed that the most accurate Δ G s m d for the alkylation reactions of nine investigated chemical carcinogens with guanine can be obtained using the DFT functional M11-L in combination with the flexible 6-311++G(d,p) basis set and implicit solvation model SMD. The mean absolute difference between the experimental and calculated activation-free energies with the M11-L/SMD method for the reactions between the nine chemical carcinogens and guanine was approximately 1 kcal/mol (1.08 kcal/mol), which corresponds to an acceptable computational error. In addition, the MN12-L functional in combination with the flexible 6-311++G(d,p) basis set and SMD solvation model reproduced the experimental values for the reactions between the same nine chemical carcinogens and guanine with an acceptable mean computational error (1.11 kcal/mol) as well. Therefore, the DFT functionals M11-L/SMD and MN12-L/SMD were also applied to calculate the corresponding Δ G s m d for the reactions of nine investigated chemical carcinogens with EGCG, EGC, and (+)-catechin, as well as to calculate Δ G s m d for the reactions of the same nine chemical carcinogens with a reference compound glutathione. Δ G s m d values calculated with the M11-L/SMD and MN12-L/SMD methods are also expected to be in very good agreement with the experimental Δ G e x p values for the reactions of the same nine chemical carcinogens with EGCG, EGC, (+)-catechin, and glutathione. The calculated Δ G s m d values for the reactions of EGCG, EGC, and (+)-catechin with the nine chemical carcinogens were then compared to the calculated Δ G s m d values for the reactions of the same chemical carcinogens with guanine and glutathione. A lower Δ G s m d indicates a faster reaction of the polyphenolic compound with the chemical carcinogen. We hypothesized that the Δ G s m d of the reaction between the nine chemical carcinogens and EGCG, EGC, or (+)-catechin is lower than the Δ G s m d of the competing reaction between the same chemical carcinogens and guanine (Figure 3 in the Section 2).
In Table 2, the calculated Δ G s m d with the M11-L/SMD method for alkylation reactions between the nine investigated chemical carcinogens and EGCG, EGC, (+)-catechin as well as glutathione, the lowest vibrational frequencies of reactant states, imaginary frequencies of transition states, the corresponding distances between the reactive centers in reactants and transition states, and relative CSA values r g l u t a t h i o n e are collected. The calculated Δ G g u a n i n e values for the reactions between the nine chemical carcinogens and guanine are also provided. The corresponding values obtained with the MN12-L functional in combination with the flexible 6-311++G(d,p) basis set and SMD implicit solvation model are reported in Table S11 of the Supplementary Materials.
As can be observed from Table 2, the correctly optimized structures of the reactant state must exhibit only real vibrational frequencies. The lowest obtained real vibrational frequencies of reactants EGCG, EGC, (+)-catechin, and glutathione with nine investigated chemical carcinogens were between 3.85 cm−1 and 27.95 cm−1. The correctly optimized transition state structure is characterized by exactly one imaginary vibrational frequency, which in normal mode must represent the reaction coordinate, corresponding to the cleavage of the covalent bond inside the epoxy ring of the investigated chemical carcinogen and the formation of a new covalent bond between the nonchiral carbon atom of the chemical carcinogen and the most nucleophilic oxygen atom of EGCG, EGC, and (+)-catechin or the sulfur atom on glutathione. The obtained single imaginary frequencies ranging from 143.49 i cm−1 to 712.15 i cm−1 indeed corresponded to the formation of the new covalent bond between the most nucleophilic oxygen of EGCG, EGC, and (+)-catechin or the sulfur on glutathione and the achiral epoxy carbon of all nine investigated chemical carcinogens (Table 2).
Much weaker intermolecular interactions in the reactant than in the transition states are reflected in larger values of dR compared to dTS, indicating a much shallower potential hypersurface in the case of reactants. The distances between the reactive centers of the reactants lie between 3.06 Å and 4.82 Å and are approximately 1 Å to 2 Å shorter than in the corresponding transition states, with distances ranging from 1.99 Å to 3.06 Å. From Table 2, it can also be observed that the distances between reactive centers in transition states were between 0.5 Å and 0.7 Å shorter for the reactions of all investigated chemical carcinogens with the three studied catechins than in the case of glutathione due to a larger van der Waals radius of the nucleophilic sulfur.
The predicted activation-free energies using the M11-L/SMD method were in the range of 15.70 kcal/mol to 25.71 kcal/mol, 11.64 kcal/mol to 22.12 kcal/mol, and 8.02 kcal/mol to 20.32 kcal/mol for the reactions of chemical carcinogens with EGCG, EGC, and (+)-catechin, respectively. Among all investigated chemical carcinogens, the Δ G s m d obtained with the M11-L/SMD method is the lowest for the reaction of chloroethylene oxide with (+)-catechin (8.02 kcal/mol). Therefore, the reaction of chloroethylene oxide with (+)-catechin is the fastest. The Δ G s m d obtained with the M11-L/SMD is the highest for the reaction of EGCG with propylene oxide (25.71 kcal/mol). Therefore, this reaction is the slowest.
Moreover, quantum mechanical calculations with the M11-L/SMD method predicted significantly lower activation-free energies (by 2 to 10 kcal/mol) for reactions of all investigated chemical carcinogens with (+)-catechin and EGC than for the competing reactions with the most reactive DNA base, guanine. EGCG demonstrated lower activation-free energies (by >1 kcal/mol) for reactions with six investigated chemical carcinogens than guanine and similar activation-free energies for reactions with propylene oxide, ethylene oxide, and AFB1 exo-8.9-epoxide (difference <1 kcal/mol).
The calculated ratios r g l u t a t h i o n e were higher than one for reactions of all nine investigated chemical carcinogens with EGCG, six investigated chemical carcinogens with EGC, and five investigated chemical carcinogens with (+)-catechin, indicating a higher reactivity of glutathione. In the case of reactions of EGC with glycidamide, vinyl carbamate epoxide, and beta-propiolactone as well as of (+)-catechin with glycidamide, beta-propiolactone, chloroethylene oxide, and 2-cyanoethylene oxide, the calculated ratios r g l u t a t h i o n e were lower than one, indicating a higher relative scavenging potential of EGC and (+)-catechin towards these chemical carcinogens than glutathione.
From the results collected in Table 2, the following general order of reactivity with the nine investigated chemical carcinogens can be established: (+)-catechin > EGC > EGCG. Nonepistructured (+)-catechin, therefore, represents the most efficient scavenger of investigated chemical carcinogens, although epistructured EGC and EGCG also represent potential carcinogen scavengers. The observed order of reactivity may in part be a consequence of the least sterically hindered structure of nonepistructured (+)-catechin, which allowed for the closest contact with investigated chemical carcinogens, whereas the reactions of chemical carcinogens with larger epistructured catechins EGC and EGCG were more sterically hindered. Graphical representations of transition state structures for the reactions of the most efficient polyphenolic scavenger (+)-catechin and the reference compound glutathione with investigated chemical carcinogens are depicted in Figure 7.
From Figure 7, the formation of a new covalent bond between the most nucleophilic phenolic oxygen atom of (+)-catechin or the sulfur atom of glutathione and the achiral epoxy carbon of nine investigated chemical carcinogens can be observed. The correct transition state structure allocation was also confirmed by the simultaneous cleavage of the covalent bond connecting this achiral carbon with the epoxy oxygen of the chemical carcinogen. The exclusively obtained single imaginary vibrational frequencies (Table 2) corresponding to this reaction path additionally confirm the validity of the proposed SN2 reaction mechanism for alkylation reactions of polyphenolic compounds with investigated chemical carcinogens [43,44].
Based on the obtained results, we can safely assume that the reactions between a large majority of investigated chemical carcinogens and (+)-catechin, EGC, as well as EGCG, would occur faster than the competing reactions between guanine and the same chemical carcinogens, meaning that the studied catechins could successfully prevent chemical carcinogen-induced DNA adduct formation. The obtained results also revealed a novel possible mechanism of the chemoprotective effects of EGCG, EGC, and (+)-catechin against nine investigated chemical carcinogens, potentially contributing to the prevention of chemical carcinogen–DNA adduct formation.

4. Conclusions

In this study, a novel quantum mechanics-based QM-CSA protocol to evaluate and compare the potential of polyphenolic compounds as scavengers of chemical carcinogens of the epoxy type under physiological conditions is developed. The QM-CSA protocol is based on kinetic calculations and represents a unique and quantitative approach to predicting the carcinogen-scavenging activity (CSA) of polyphenolic compounds by applying two scales: the absolute scale based on the comparison of activation-free energies of a series of polyphenols with guanine and a relative scale based on the comparison with glutathione as a reference compound. The reliability of QM-CSA was demonstrated by comparison of the calculated activation-free energies with the corresponding experimental activation-free energies for alkylation reactions between nine investigated chemical carcinogens and the most reactive DNA base, guanine.
The importance of the selection of an appropriate combination of DFT functionals, namely M11-L or MN12-L, flexible basis set 6-311++G(d,p), and implicit solvation model SMD to achieve the best agreement with experimental data under physiological conditions was demonstrated. The functionals M11-L and MN12-L, in conjunction with the SMD solvation model, were identified as the most reliable DFT methods to evaluate the reactivity of nine investigated chemical carcinogens of the epoxy type with guanine among the 14 tested functionals. The uncertainties of the proposed QM-CSA methodology have been proven to be similar to experimental ones with mean absolute errors of 1.08 kcal/mol and 1.11 kcal/mol for M11-L/SMD and MN12-L/SMD, respectively.
Moreover, the QM-CSA protocol was also applied to evaluate the CSA of a series of three catechins, namely EGCG, EGC, and (+)-catechin, and to compare their potential as scavengers of nine investigated chemical carcinogens at physiological conditions with guanine and glutathione. The following order of reactivity in the series of three catechins was determined: (+)-catechin > EGC > EGCG, indicating that (+)-catechin represents the most efficient scavenger of the investigated chemical carcinogens in comparison with guanine and glutathione. The proposed SN2 reaction mechanism was also confirmed.
To conclude, the QM-CSA protocol allows for the estimation of activation-free energies, the identification of polyphenolic compounds with the highest carcinogen scavenging potential, the prediction of order/trend of scavenging activity, and insights into the molecular mechanisms of the studied alkylation reactions. The applications of the QM-CSA protocol are intended to increase the kinetic data on alkylation reactions between polyphenolic compounds and chemical carcinogens of the epoxy type, which are currently lacking in the scientific literature, and to facilitate the design of efficient novel pharmacological strategies for cancer prevention.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/foods13172708/s1, Tables S1–S9. The activation-free energies for the reaction between nine investigated chemical carcinogens and guanine calculated with 14 DFT functionals in conjunction with the 6-311++G(d,p) flexible basis set and the SMD or CPCM solvation model, together with the corresponding absolute and relative errors, as well as the frequencies corresponding to reactant and transition state structures; Table S10. The average absolute and relative errors for the reactions between nine investigated chemical carcinogens and guanine calculated with 14 tested DFT functionals in conjunction with the SMD and CPCM solvation models; Table S11. Computational results obtained with the MN12-L functional in conjunction with the 6-311++G(d,p) flexible basis set and SMD implicit solvation model for the reactions of nine investigated chemical carcinogens with EGCG, EGC, (+)-catechin, and glutathione.

Author Contributions

Conceptualization, V.F., J.T., and U.B.; methodology, V.F., J.T., and U.B.; validation, V.F.; formal analysis, V.F.; investigation, V.F., and J.T.; data curation, V.F. and J.T., resources, U.B.; Writing—Original draft preparation, V.F.; Writing—Review and editing, V.F., J.T., and U.B.; visualization, V.F.; supervision, U.B.; funding acquisition, U.B. All authors have read and agreed to the published version of the manuscript.

Funding

Financial support was provided through the Slovenian Research and Innovation Agency (ARIS) program, and project grants Z4-2654, J1-2471, P2-0046, L2-3175, J4-4633, J1-4398, L2-4430, J3-4498, J7-4638, J1-4414, J3-4497, P2-0438, and I0-E015 are gratefully acknowledged.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The original contributions presented in the study are included in the article and Supplementary Material, further inquiries can be directed to the corresponding authors.

Acknowledgments

The authors gratefully acknowledge the HPC RIVR consortium for computing resources HPC Vega.

Conflicts of Interest

The authors declare no conflicts of interest.

Appendix A

The proposed SN2 mechanisms for the alkylation reactions of nine investigated chemical carcinogens with (a) guanine and (b) glutathione are depicted in Scheme A1, Scheme A2, Scheme A3, Scheme A4, Scheme A5, Scheme A6, Scheme A7, Scheme A8 and Scheme A9.
Scheme A1. The proposed SN2 mechanism for the reaction of styrene oxide with (a) guanine and (b) glutathione.
Scheme A1. The proposed SN2 mechanism for the reaction of styrene oxide with (a) guanine and (b) glutathione.
Foods 13 02708 sch0a1
Scheme A2. The proposed SN2 mechanism for the reaction of propylene oxide with (a) guanine and (b) glutathione.
Scheme A2. The proposed SN2 mechanism for the reaction of propylene oxide with (a) guanine and (b) glutathione.
Foods 13 02708 sch0a2
Scheme A3. The proposed SN2 mechanism for the reaction of ethylene oxide with (a) guanine and (b) glutathione.
Scheme A3. The proposed SN2 mechanism for the reaction of ethylene oxide with (a) guanine and (b) glutathione.
Foods 13 02708 sch0a3
Scheme A4. The proposed SN2 mechanism for the reaction of glycidamide with (a) guanine and (b) glutathione. * represents the chiral center.
Scheme A4. The proposed SN2 mechanism for the reaction of glycidamide with (a) guanine and (b) glutathione. * represents the chiral center.
Foods 13 02708 sch0a4
Scheme A5. The proposed SN2 mechanism for the reaction of vinyl carbamate epoxide with (a) guanine and (b) glutathione. * represents the chiral center.
Scheme A5. The proposed SN2 mechanism for the reaction of vinyl carbamate epoxide with (a) guanine and (b) glutathione. * represents the chiral center.
Foods 13 02708 sch0a5
Scheme A6. The proposed SN2 mechanism for the reaction of beta-propiolactone with (a) guanine and (b) glutathione.
Scheme A6. The proposed SN2 mechanism for the reaction of beta-propiolactone with (a) guanine and (b) glutathione.
Foods 13 02708 sch0a6
Scheme A7. The proposed SN2 mechanism for the reaction of chloroethylene oxide with (a) guanine and (b) glutathione. * represents the chiral center.
Scheme A7. The proposed SN2 mechanism for the reaction of chloroethylene oxide with (a) guanine and (b) glutathione. * represents the chiral center.
Foods 13 02708 sch0a7
Scheme A8. The proposed SN2 mechanism for the reaction of 2-cyanoethylene oxide with (a) guanine and (b) glutathione. * represents the chiral center.
Scheme A8. The proposed SN2 mechanism for the reaction of 2-cyanoethylene oxide with (a) guanine and (b) glutathione. * represents the chiral center.
Foods 13 02708 sch0a8
Scheme A9. The proposed SN2 mechanism for the reaction of AFB1 exo-8,9-epoxide with (a) guanine and (b) glutathione.
Scheme A9. The proposed SN2 mechanism for the reaction of AFB1 exo-8,9-epoxide with (a) guanine and (b) glutathione.
Foods 13 02708 sch0a9

References

  1. Tomasetti, C.; Li, L.; Vogelstein, B. Stem cell divisions, somatic mutations, cancer etiology, and cancer prevention. Science 2017, 355, 1330–1334. [Google Scholar] [CrossRef] [PubMed]
  2. Bren, U.; Guengerich, F.P.; Mavri, J. Guanine alkylation by the potent carcinogen aflatoxin B1: Quantum chemical calculations. Chem. Res. Toxicol. 2007, 20, 1134–1140. [Google Scholar] [CrossRef] [PubMed]
  3. Kranjc, A.; Mavri, J. Guanine alkylation by ethylene oxide: Calculation of chemical reactivity. J. Phys. Chem. A 2006, 110, 5740–5744. [Google Scholar] [CrossRef] [PubMed]
  4. Bren, U.; Zupan, M.; Guengerich, F.P.; Mavri, J. Chemical reactivity as a tool to study carcinogenicity: Reaction between chloroethylene oxide and guanine. J. Org. Chem. 2006, 71, 4078–4084. [Google Scholar] [CrossRef]
  5. Lajovic, A.; Nagy, L.D.; Guengerich, F.P.; Bren, U. Carcinogenesis of urethane: Simulation versus experiment. Chem. Res. Toxicol. 2015, 28, 691. [Google Scholar] [CrossRef] [PubMed]
  6. Galesa, K.; Bren, U.; Kranjc, A.; Mavri, J. Carcinogenicity of acrylamide: A computational study. J. Agric. Food Chem. 2008, 56, 8720–8727. [Google Scholar] [CrossRef]
  7. Kržan, M.; Mavri, J. Carcinogenicity of styrene oxide: Calculation of chemical reactivity. Croat. Chem. Acta 2009, 82, 317–322. [Google Scholar]
  8. Mavri, J. Can the chemical reactivity of an ultimate carcinogen be related to its carcinogenicity? An application to propylene oxide. Toxicol. In Vitro 2013, 27, 479–485. [Google Scholar] [CrossRef] [PubMed]
  9. Gladović, M.; Španinger, E.; Bren, U. Nucleic bases alkylation with acrylonitrile and cyanoethylene oxide: A computational study. Chem. Res. Toxicol. 2017. [Google Scholar] [CrossRef] [PubMed]
  10. Španinger, E.; Bren, U. Carcinogenesis of β-propiolactone: A computational study. Chem. Res. Toxicol. 2020, 33, 769–781. [Google Scholar] [CrossRef]
  11. Pompella, A.; Visvikis, A.; Paolicchi, A.; De Tata, V.; Casini, A.F. The changing faces of glutathione, a cellular protagonist. Biochem. Pharmacol. 2003, 66, 1499–1503. [Google Scholar] [CrossRef] [PubMed]
  12. Kennedy, L.; Sandhu, J.K.; Harper, M.-E.; Cuperlovic-Culf, M. Role of glutathione in cancer: From mechanisms to therapies. Biomolecules 2020, 10, 1429. [Google Scholar] [CrossRef] [PubMed]
  13. González-Vallinas, M.; González-Castejón, M.; Rodríguez-Casado, A.; de Molina, A.R. Dietary phytochemicals in cancer prevention and therapy: A complementary approach with promising perspectives. Nutr. Rev. 2013, 71, 585–599. [Google Scholar] [CrossRef]
  14. Li, A.-N.; Li, S.; Zhang, Y.-J.; Xu, X.-R.; Chen, Y.-M.; Li, H.-B. Resources and biological activities of natural polyphenols. Nutrients 2014, 6, 6020–6047. [Google Scholar] [CrossRef] [PubMed]
  15. Brglez Mojzer, E.; Knez Hrnčič, M.; Škerget, M.; Knez, Ž.; Bren, U. Polyphenols: Extraction methods, antioxidative action, bioavailability and anticarcinogenic effects. Molecules 2016, 21, 901. [Google Scholar] [CrossRef]
  16. Cerezo-Guisado, M.I.; Zur, R.; Lorenzo, M.J.; Risco, A.; Martín-Serrano, M.A.; Alvarez-Barrientos, A.; Cuenda, A.; Centeno, F. Implication of Akt, ERK1/2 and alternative p38MAPK signalling pathways in human colon cancer cell apoptosis induced by green tea EGCG. Food Chem. Toxicol. 2015, 84, 125–132. [Google Scholar] [CrossRef] [PubMed]
  17. Cheng, Z.; Zhang, Z.; Han, Y.; Wang, J.; Wang, Y.; Chen, X.; Shao, Y.; Cheng, Y.; Zhou, W.; Lu, X. A review on anti-cancer effect of green tea catechins. J. Funct. Foods 2020, 74, 104172. [Google Scholar] [CrossRef]
  18. Chen, X.; Du, Y.; Wu, L.; Xie, J.; Chen, X.; Hu, B.; Wu, Z.; Yao, Q.; Li, Q. Effects of tea-polysaccharide conjugates and metal ions on precipitate formation by epigallocatechin gallate and caffeine, the key components of green tea infusion. J. Agric. Food Chem. 2019, 67, 3744–3751. [Google Scholar] [CrossRef] [PubMed]
  19. Haratifar, S.; Meckling, K.; Corredig, M. Antiproliferative activity of tea catechins associated with casein micelles, using HT29 colon cancer cells. J. Dairy Sci. 2014, 97, 672–678. [Google Scholar] [CrossRef]
  20. Sorrenti, V.; Ali, S.; Mancin, L.; Davinelli, S.; Paoli, A.; Scapagnini, G. Cocoa polyphenols and gut microbiota interplay: Bioavailability, prebiotic effect, and impact on human health. Nutrients 2020, 12, 1908. [Google Scholar] [CrossRef]
  21. Cao, J.; Han, J.; Xiao, H.; Qiao, J.; Han, M. Effect of tea polyphenol compounds on anticancer drugs in terms of anti-tumor activity, toxicology, and pharmacokinetics. Nutrients 2016, 8, 762. [Google Scholar] [CrossRef] [PubMed]
  22. Rahmani, A.H.; Allemailem, K.S.; Aly, S.M.; Khan, M.A. Implications of green tea and its constituents in the prevention of cancer via the modulation of cell signalling pathway. BioMed Res. Int. 2015, 2015, 925640. [Google Scholar] [CrossRef]
  23. Shirakami, Y.; Shimizu, M. Possible mechanisms of green tea and its constituents against cancer. Molecules 2018, 23, 2284. [Google Scholar] [CrossRef] [PubMed]
  24. Suganuma, M.; Takahashi, A.; Watanabe, T.; Iida, K.; Matsuzaki, T.; Yoshikawa, H.Y.; Fujiki, H. Biophysical approach to mechanisms of cancer prevention and treatment with green tea catechins. Molecules 2016, 21, 1566. [Google Scholar] [CrossRef] [PubMed]
  25. Andújar, I.; Recio, M.C.; Giner, R.M.; Ríos, J. Cocoa polyphenols and their potential benefits for human health. Oxid. Med. Cell. Longev. 2012, 2012, 906252. [Google Scholar] [CrossRef]
  26. Singh, B.N.; Shankar, S.; Srivastava, R.K. Green tea catechin, epigallocatechin-3-gallate (EGCG): Mechanisms, perspectives and clinical applications. Biochem. Pharmacol. 2011, 82, 1807–1821. [Google Scholar] [CrossRef] [PubMed]
  27. Yang, C.S.; Lambert, J.D.; Sang, S. Antioxidative and anti-carcinogenic activities of tea polyphenols. Arch. Toxicol. 2009, 83, 11–21. [Google Scholar] [CrossRef]
  28. Ohishi, T.; Goto, S.; Monira, P.; Isemura, M.; Nakamura, Y. Anti-inflammatory action of green tea. Antiinflamm. Antiallergy Agents Med. Chem. 2016, 15, 74–90. [Google Scholar] [CrossRef]
  29. Reygaert, W.C. The antimicrobial possibilities of green tea. Front. Microbiol. 2014, 5, 434. [Google Scholar] [CrossRef]
  30. Ann Beltz, L.; Kay Bayer, D.; Lynn Moss, A.; Mitchell Simet, I. Mechanisms of cancer prevention by green and black tea polyphenols. Curr. Med. Chem. Anticancer Agents 2006, 6, 389–406. [Google Scholar] [CrossRef] [PubMed]
  31. Kundu, J.K.; Na, H.-K.; Chun, K.-S.; Kim, Y.-K.; Lee, S.J.; Lee, S.S.; Lee, O.-S.; Sim, Y.-C.; Surh, Y.-J. Inhibition of phorbol ester–induced COX-2 expression by epigallocatechin gallate in mouse skin and cultured human mammary epithelial cells. J. Nutr. 2003, 133, 3805S–3810S. [Google Scholar] [CrossRef]
  32. Fechtner, S.; Ahmed, S. In vitro Comparison of Anti-Inflammatory Activity of Green Tea Catechins. FASEB J. 2016, 30, 916.4. [Google Scholar] [CrossRef]
  33. Rathore, K.; Choudhary, S.; Odoi, A.; Wang, H.-C.R. Green tea catechin intervention of reactive oxygen species-mediated ERK pathway activation and chronically induced breast cell carcinogenesis. Carcinogenesis 2012, 33, 174–183. [Google Scholar] [CrossRef]
  34. Ho, H.C.; Huang, C.C.; Lu, Y.T.; Yeh, C.M.; Ho, Y.T.; Yang, S.F.; Hsin, C.H.; Lin, C.W. Epigallocatechin-3-gallate inhibits migration of human nasopharyngeal carcinoma cells by repressing MMP-2 expression. J. Cell. Physiol. 2019, 234, 20915–20924. [Google Scholar] [CrossRef] [PubMed]
  35. Sen, T.; Chatterjee, A. Epigallocatechin-3-gallate (EGCG) downregulates EGF-induced MMP-9 in breast cancer cells: Involvement of integrin receptor α5β1 in the process. Eur. J. Nutr. 2011, 50, 465–478. [Google Scholar] [CrossRef] [PubMed]
  36. Gu, J.-W.; Makey, K.L.; Tucker, K.B.; Chinchar, E.; Mao, X.; Pei, I.; Thomas, E.Y.; Miele, L. EGCG, a major green tea catechin suppresses breast tumor angiogenesis and growth via inhibiting the activation of HIF-1α and NFκB, and VEGF expression. Vascular Cell 2013, 5, 9. [Google Scholar] [CrossRef] [PubMed]
  37. Zhang, Y.; Duan, W.; Owusu, L.; Wu, D.; Xin, Y. Epigallocatechin-3-gallate induces the apoptosis of hepatocellular carcinoma LM6 cells but not non-cancerous liver cells. Int. J. Mol. Med. 2015, 35, 117–124. [Google Scholar] [CrossRef] [PubMed]
  38. Kim, H.S.; Kim, M.H.; Jeong, M.; Hwang, Y.S.; Lim, S.H.; Shin, B.A.; Ahn, B.W.; Jung, Y.D. EGCG blocks tumor promoter-induced MMP-9 expression via suppression of MAPK and AP-1 activation in human gastric AGS cells. Anticancer Res. 2004, 24, 747–754. [Google Scholar]
  39. Guruvayoorappan, C.; Kuttan, G. (+)-Catechin inhibits tumour angiogenesis and regulates the production of nitric oxide and TNF-α in LPS-stimulated macrophages. Innate Immun. 2008, 14, 160–174. [Google Scholar] [CrossRef]
  40. Manikandan, R.; Beulaja, M.; Arulvasu, C.; Sellamuthu, S.; Dinesh, D.; Prabhu, D.; Babu, G.; Vaseeharan, B.; Prabhu, N. Synergistic anticancer activity of curcumin and catechin: An in vitro study using human cancer cell lines. Microsc. Res. Tech. 2012, 75, 112–116. [Google Scholar] [CrossRef]
  41. Carla Cadoná, F.; Kolinski Machado, A.; Farina Azzolin, V.; Barbisan, F.; Bortoluzzi Dornelles, E.; Glanzner, W.; Bayard Gonçalves, P.D.; Elias Assmann, C.; Esteves Ribeiro, E.; Beatrice Mânica da Cruz, I. Guaraná a caffeine-rich food increases oxaliplatin sensitivity of colorectal HT-29 cells by apoptosis pathway modulation. Curr. Med. Chem. Anticancer Agents 2016, 16, 1055–1065. [Google Scholar] [CrossRef] [PubMed]
  42. Sun, H.; Yin, M.; Hao, D.; Shen, Y. Anti-cancer activity of catechin against A549 lung carcinoma cells by induction of cyclin kinase inhibitor P21 and suppression of cyclin E1 and P–AKT. Appl. Sci. 2020, 10, 2065. [Google Scholar] [CrossRef]
  43. Furlan, V.; Bren, U. Protective Effects of [6]-Gingerol Against Chemical Carcinogens: Mechanistic Insights. Int. J. Mol. Sci. 2020, 21, 695. [Google Scholar] [CrossRef] [PubMed]
  44. Štern, A.; Furlan, V.; Novak, M.; Štampar, M.; Kolenc, Z.; Kores, K.; Filipič, M.; Bren, U.; Žegura, B. Chemoprotective Effects of Xanthohumol against the Carcinogenic Mycotoxin Aflatoxin B1. Foods 2021, 10, 1331. [Google Scholar] [CrossRef] [PubMed]
  45. Furlan, V.; Bren, U. Helichrysum italicum: From Extraction, Distillation, and Encapsulation Techniques to Beneficial Health Effects. Foods 2023, 12, 802. [Google Scholar] [CrossRef] [PubMed]
  46. Muzolf, M.; Szymusiak, H.; Gliszczyńska-Świgło, A.; Rietjens, I.M.; Tyrakowska, B.E. pH-dependent radical scavenging capacity of green tea catechins. J. Agric. Food Chem. 2008, 56, 816–823. [Google Scholar] [CrossRef]
  47. Zhao, Y.; Schultz, N.E.; Truhlar, D.G. Design of density functionals by combining the method of constraint satisfaction with parametrization for thermochemistry, thermochemical kinetics, and noncovalent interactions. J. Chem. Theory Comput. 2006, 2, 364–382. [Google Scholar] [CrossRef] [PubMed]
  48. Zhao, Y.; Truhlar, D.G. The M06 suite of density functionals for main group thermochemistry, thermochemical kinetics, noncovalent interactions, excited states, and transition elements: Two new functionals and systematic testing of four M06-class functionals and 12 other functionals. Theor. Chem. Acc. 2008, 120, 215–241. [Google Scholar]
  49. Zhao, Y.; Truhlar, D.G. Density functional for spectroscopy: No long-range self-interaction error, good performance for Rydberg and charge-transfer states, and better performance on average than B3LYP for ground states. J. Phys. Chem. A 2006, 110, 13126–13130. [Google Scholar] [CrossRef]
  50. Zhao, Y.; Truhlar, D.G. Exploring the limit of accuracy of the global hybrid meta density functional for main-group thermochemistry, kinetics, and noncovalent interactions. J. Chem. Theory Comput. 2008, 4, 1849–1868. [Google Scholar] [CrossRef]
  51. Peverati, R.; Truhlar, D.G. Improving the accuracy of hybrid meta-GGA density functionals by range separation. J. Phys. Chem. Lett. 2011, 2, 2810–2817. [Google Scholar] [CrossRef]
  52. Peverati, R.; Truhlar, D.G. M11-L: A local density functional that provides improved accuracy for electronic structure calculations in chemistry and physics. J. Phys. Chem. Lett. 2012, 3, 117–124. [Google Scholar] [CrossRef]
  53. Peverati, R.; Truhlar, D.G. An improved and broadly accurate local approximation to the exchange–correlation density functional: The MN12-L functional for electronic structure calculations in chemistry and physics. Phys. Chem. Chem. Phys. 2012, 14, 13171–13174. [Google Scholar] [CrossRef] [PubMed]
  54. Peverati, R.; Truhlar, D.G. Screened-exchange density functionals with broad accuracy for chemistry and solid-state physics. Phys. Chem. Chem. Phys. 2012, 14, 16187–16191. [Google Scholar] [CrossRef] [PubMed]
  55. Haoyu, S.Y.; He, X.; Li, S.L.; Truhlar, D.G. MN15: A Kohn–Sham global-hybrid exchange–correlation density functional with broad accuracy for multi-reference and single-reference systems and noncovalent interactions. Chem. Sci. 2016, 7, 5032–5051. [Google Scholar]
  56. Yu, H.S.; He, X.; Truhlar, D.G. MN15-L: A new local exchange-correlation functional for Kohn–Sham density functional theory with broad accuracy for atoms, molecules, and solids. J. Chem. Theory Comput. 2016, 12, 1280–1293. [Google Scholar] [CrossRef]
  57. Peverati, R.; Truhlar, D.G. Exchange–correlation functional with good accuracy for both structural and energetic properties while depending only on the density and its gradient. J. Chem. Theory Comput. 2012, 8, 2310–2319. [Google Scholar] [CrossRef] [PubMed]
  58. Schneebeli, S.T.; Bochevarov, A.D.; Friesner, R.A. Parameterization of a B3LYP specific correction for noncovalent interactions and basis set superposition error on a gigantic data set of CCSD(T) quality noncovalent interaction energies. J. Chem. Theory Comput. 2011, 7, 658–668. [Google Scholar] [CrossRef]
  59. Chai, J.-D.; Head-Gordon, M. Long-range corrected hybrid density functionals with damped atom–atom dispersion corrections. Phys. Chem. Chem. Phys. 2008, 10, 6615–6620. [Google Scholar] [CrossRef]
  60. Eyring, H. The activated complex in chemical reactions. J. Chem. Phys. 1935, 3, 107–115. [Google Scholar] [CrossRef]
  61. Peverati, R.; Truhlar, D.G. Quest for a universal density functional: The accuracy of density functionals across a broad spectrum of databases in chemistry and physics. Philos. Trans. A Math. Phys. Eng. Sci. 2014, 372, 20120476. [Google Scholar] [CrossRef] [PubMed]
  62. Jana, K.; Ghosh, S.; Wakchaure, P.D.; Bandyopadhyay, T.; Ganguly, B. Probing the role of imidazopyridine and imidazophosphorine scaffolds to design novel proton pump inhibitor for H+, K+-ATPase: A DFT study. ACS Omega 2019, 4, 1311–1321. [Google Scholar] [CrossRef]
  63. Mardirossian, N.; Head-Gordon, M. How accurate are the Minnesota density functionals for noncovalent interactions, isomerization energies, thermochemistry, and barrier heights involving molecules composed of main-group elements? J. Chem. Theory Comput. 2016, 12, 4303–4325. [Google Scholar] [CrossRef] [PubMed]
  64. Wiberg, K.B. Basis set effects on calculated geometries: 6-311++G** vs. aug-cc-pVDZ. J. Comput. Chem. 2004, 25, 1342–1346. [Google Scholar] [CrossRef]
  65. Tomasi, J.; Persico, M. Molecular interactions in solution: An overview of methods based on continuous distributions of the solvent. Chem. Rev. 1994, 94, 2027–2094. [Google Scholar] [CrossRef]
  66. Cramer, C.J.; Truhlar, D.G. Implicit solvation models: Equilibria, structure, spectra, and dynamics. Chem. Rev. 1999, 99, 2161–2200. [Google Scholar] [CrossRef]
  67. Marenich, A.V.; Cramer, C.J.; Truhlar, D.G. Universal solvation model based on solute electron density and on a continuum model of the solvent defined by the bulk dielectric constant and atomic surface tensions. J. Phys. Chem. B 2009, 113, 6378–6396. [Google Scholar] [CrossRef]
  68. MarvinSketch; Version 20.19; ChemAxon: Budapest, Hungary, 2021; Available online: https://www.chemaxon.com/products/marvin/marvinsketch/ (accessed on 20 August 2024).
  69. Hayes, J.D.; McLellan, L.I. Glutathione and glutathione-dependent enzymes represent a co-ordinately regulated defence against oxidative stress. Free Radic. Res. 1999, 31, 273–300. [Google Scholar] [CrossRef] [PubMed]
  70. Dringen, R.; Brandmann, M.; Hohnholt, M.C.; Blumrich, E.-M. Glutathione-dependent detoxification processes in astrocytes. Neurochem. Res. 2015, 40, 2570–2582. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Biotransformations of nine procarcinogens into investigated genotoxic ultimate chemical carcinogens, among which eight require metabolic activation by cytochrome P450 enzymes: (a) aflatoxin B1 exo-8,9-epoxide, (b) ethylene oxide, (c) chloroethylene oxide, (d) vinyl carbamate epoxide, (e) glycidamide, (f) styrene oxide, (g) propylene oxide, and (h) 2-cyanoethylene oxide, while (i) beta-propiolactone needs no metabolic activation by cytochrome P450 enzymes and acts as a genotoxic chemical carcinogen by itself.
Figure 1. Biotransformations of nine procarcinogens into investigated genotoxic ultimate chemical carcinogens, among which eight require metabolic activation by cytochrome P450 enzymes: (a) aflatoxin B1 exo-8,9-epoxide, (b) ethylene oxide, (c) chloroethylene oxide, (d) vinyl carbamate epoxide, (e) glycidamide, (f) styrene oxide, (g) propylene oxide, and (h) 2-cyanoethylene oxide, while (i) beta-propiolactone needs no metabolic activation by cytochrome P450 enzymes and acts as a genotoxic chemical carcinogen by itself.
Foods 13 02708 g001
Figure 2. Structural formulas of (a) epigallocatechin-3-gallate (EGCG), (b) epigallocatechin (EGC), and (c) (+)-catechin.
Figure 2. Structural formulas of (a) epigallocatechin-3-gallate (EGCG), (b) epigallocatechin (EGC), and (c) (+)-catechin.
Foods 13 02708 g002
Figure 3. Activation-free energy, Δ G , for the reaction between the ultimate chemical carcinogen (AFB1 exo-8,9-epoxide) and its polyphenolic scavenger ((+)-catechin) (green arrow) in comparison with the activation-free energy, Δ G , for a competing reaction between the same chemical carcinogen and guanine (red arrow).
Figure 3. Activation-free energy, Δ G , for the reaction between the ultimate chemical carcinogen (AFB1 exo-8,9-epoxide) and its polyphenolic scavenger ((+)-catechin) (green arrow) in comparison with the activation-free energy, Δ G , for a competing reaction between the same chemical carcinogen and guanine (red arrow).
Foods 13 02708 g003
Figure 4. The heat map of differences Δ Δ G = ( Δ G c a l Δ G e x p ) between the activation-free energies calculated with 14 DFT functionals in combination with the flexible 6-311++G(d,p) basis set and implicit solvation models CPCM or SMD ( Δ G c a l ), and the corresponding experimental activation-free energies ( Δ G e x p ) for reactions between the nine investigated chemical carcinogens and guanine. The dark green color designates a good agreement with the corresponding experimental activation-free energy for the alkylation reaction between individual chemical carcinogens and guanine (ideally Δ Δ G = ( Δ G c a l Δ G e x p ) = 0 ), while the light green color designates a poor agreement with the corresponding experimental value ( Δ Δ G >> 0 or Δ Δ G << 0).
Figure 4. The heat map of differences Δ Δ G = ( Δ G c a l Δ G e x p ) between the activation-free energies calculated with 14 DFT functionals in combination with the flexible 6-311++G(d,p) basis set and implicit solvation models CPCM or SMD ( Δ G c a l ), and the corresponding experimental activation-free energies ( Δ G e x p ) for reactions between the nine investigated chemical carcinogens and guanine. The dark green color designates a good agreement with the corresponding experimental activation-free energy for the alkylation reaction between individual chemical carcinogens and guanine (ideally Δ Δ G = ( Δ G c a l Δ G e x p ) = 0 ), while the light green color designates a poor agreement with the corresponding experimental value ( Δ Δ G >> 0 or Δ Δ G << 0).
Foods 13 02708 g004
Figure 5. The correlation between the calculated activation-free energies using (a) M11-L/SMD and (b) MN12-L/SMD methods, respectively, and the experimental activation-free energies for reactions between nine investigated chemical carcinogens and guanine.
Figure 5. The correlation between the calculated activation-free energies using (a) M11-L/SMD and (b) MN12-L/SMD methods, respectively, and the experimental activation-free energies for reactions between nine investigated chemical carcinogens and guanine.
Foods 13 02708 g005
Figure 6. Transition state structures of the nucleophilic attack of the N7 atom of guanine on the achiral carbon in the epoxy ring of the nine investigated chemical carcinogens as predicted by the DFT functional M11-L in combination with the flexible 6-311++G(d,p) basis set and SMD implicit solvation model. The carbon atoms are depicted in gray; the oxygen atoms are red, the nitrogen atoms are blue; the chlorine atom is green; and the hydrogen atoms are white.
Figure 6. Transition state structures of the nucleophilic attack of the N7 atom of guanine on the achiral carbon in the epoxy ring of the nine investigated chemical carcinogens as predicted by the DFT functional M11-L in combination with the flexible 6-311++G(d,p) basis set and SMD implicit solvation model. The carbon atoms are depicted in gray; the oxygen atoms are red, the nitrogen atoms are blue; the chlorine atom is green; and the hydrogen atoms are white.
Foods 13 02708 g006
Figure 7. Transition state structures of alkylation reactions of the most efficient polyphenolic scavenger (+)-catechin and the reference compound glutathione with investigated chemical carcinogens obtained using the DFT functional M11-L in combination with the flexible 6-311++G(d,p) basis set and SMD implicit solvation model. The carbon atoms are gray; the oxygen atoms are red; the nitrogen atoms are blue; the sulfur atom is yellow; the chlorine atom is green; and the hydrogen atoms are white.
Figure 7. Transition state structures of alkylation reactions of the most efficient polyphenolic scavenger (+)-catechin and the reference compound glutathione with investigated chemical carcinogens obtained using the DFT functional M11-L in combination with the flexible 6-311++G(d,p) basis set and SMD implicit solvation model. The carbon atoms are gray; the oxygen atoms are red; the nitrogen atoms are blue; the sulfur atom is yellow; the chlorine atom is green; and the hydrogen atoms are white.
Foods 13 02708 g007
Table 1. Activation-free energies, lowest vibrational frequencies, and distances between the reactive centers of reactant and transition states, as well as absolute and relative errors together with corresponding experimental activation-free energies for the reactions between nine studied chemical carcinogens and guanine at the M11-L/SMD and MN12-L/SMD levels of theory.
Table 1. Activation-free energies, lowest vibrational frequencies, and distances between the reactive centers of reactant and transition states, as well as absolute and relative errors together with corresponding experimental activation-free energies for the reactions between nine studied chemical carcinogens and guanine at the M11-L/SMD and MN12-L/SMD levels of theory.
Functional/Solvation Model Δ G S M D
(kcal/mol) a
ωTS
(i cm−1) b
ωR
(cm−1) c
dTS
(Å) d
dR
(Å) e
| Δ ( Δ G ) |
(kcal/mol) f
Δ Δ G Δ G e x p
(/) g
Δ G e x p
(kcal/mol) h
Styrene Oxide
M11-L/SMD27.16640.0817.631.943.340.660.0226.5
MN12-L/SMD25.50577.3522.261.943.131.000.0426.5
Propylene Oxide
M11-L/SMD26.34629.2329.992.063.820.940.0425.4
MN12-L/SMD25.45569.0430.072.053.170.050.00225.4
Ethylene Oxide
M11-L/SMD24.08613.1711.772.072.790.620.0224.7
MN12-L/SMD23.29580.768.282.063.241.410.0524.7
Glycidamide
M11-L/SMD23.76654.0724.062.083.700.960.0422.8
MN12-L/SMD23.35598.0324.222.063.260.550.0222.8
Vinyl Carbamate Epoxide
M11-L/SMD20.65543.7529.152.153.491.750.0822.4
MN12-L/SMD20.95515.4425.842.133.361.450.0622.4
Beta-propiolactone
M11-L/SMD21.11583.9516.852.133.570.310.0120.8
MN12-L/SMD21.34562.6917.172.133.390.540.0320.8
Chloroethylene Oxide
M11-L/SMD18.35539.2623.652.193.591.150.0619.5
MN12-L/SMD17.64515.2936.832.163.191.860.0919.5
2-Cyanoethylene Oxide
M11-L/SMD20.65 646.5120.002.103.321.450.0819.2
MN12-L/SMD20.57616.9245.332.083.501.370.0719.2
AFB1 Exo-8.9-Epoxide
M11-L/SMD16.99225.8721.532.124.211.890.1215.1
MN12-L/SMD16.88207.0526.032.124.141.780.1215.1
a Activation-free energy is calculated with the M11-L or MN12-L functionals in combination with the flexible 6-311++G(d,p) basis set and SMD solvation model. b The exact imaginary vibrational frequency of the transition state structure. c The lowest real vibrational frequency of the reactant state structure. d The distance in the transition state structure between the most nucleophilic nitrogen atom N7 on guanine and the achiral electrophilic epoxy carbon of the chemical carcinogen. e The distance in the reactant structure between the most nucleophilic atom N7 of guanine and the chiral electrophilic epoxy carbon of chemical carcinogen. f The absolute error between the calculated and experimental activation-free energies for reactions between investigated chemical carcinogens and guanine. g The relative error between the calculated and experimental activation-free energies for reactions between investigated chemical carcinogens and guanine. h Experimental activation-free energies for reactions of investigated chemical carcinogens with guanine.
Table 2. Computational results obtained with the M11-L functional in combination with the flexible 6-311++G(d,p) basis set and SMD implicit solvation model for the alkylation reactions of the nine investigated chemical carcinogens with EGCG, EGC, (+)-catechin, and glutathione.
Table 2. Computational results obtained with the M11-L functional in combination with the flexible 6-311++G(d,p) basis set and SMD implicit solvation model for the alkylation reactions of the nine investigated chemical carcinogens with EGCG, EGC, (+)-catechin, and glutathione.
Functional/Solvation Model
M11-L/SMD
Δ G S M D
(kcal/mol) a
ωTS
(i cm−1) b
ωR
(cm−1) c
dTS
(Å) d
dR
(Å) e
r
(/) f
Δ G g u a n i n e
(kcal/mol) g
Styrene Oxide
EGCG25.65585.1610.481.993.351.5327.16
EGC19.80577.896.192.043.671.1827.16
(+)-Catechin20.32591.1211.052.033.141.2127.16
Glutathione16.81677.0915.612.604.141.0027.16
Propylene Oxide
EGCG25.71566.0314.462.003.251.5226.34
EGC20.73581.6125.132.033.411.2326.34
(+)-Catechin18.85555.3616.312.023.351.1226.34
Glutathione16.87663.7522.002.583.471.0026.34
Ethylene Oxide
EGCG24.27566.5614.552.013.731.7024.08
EGC22.12586.5927.952.053.851.5524.08
(+)-Catechin17.37574.6621.412.033.361.2224.08
Glutathione14.25675.2327.412.603.691.0024.08
Glycidamide
EGCG22.07552.2816.342.024.481.0923.76
EGC19.35562.818.972.064.640.9523.76
(+)-Catechin19.46583.3411.722.053.480.9623.76
Glutathione20.30694.059.162.574.131.0023.76
Vinyl Carbamate Epoxide
EGCG17.64564.0410.752.103.581.4720.65
EGC11.64591.4211.132.133.540.9720.65
(+)-Catechin13.61590.7721.982.123.701.1420.65
Glutathione11.97655.0514.842.654.821.0020.65
Beta-propiolactone
EGCG16.97581.6912.282.093.061.2021.11
EGC12.85615.4117.642.103.230.9121.11
(+)-Catechin10.65630.2516.372.083.240.7521.11
Glutathione14.19653.9412.492.623.651.0021.11
Chloroethylene Oxide
EGCG15.70580.3315.142.143.821.6818.35
EGC12.86607.7417.612.163.521.3718.35
(+)-Catechin8.02591.699.262.153.470.8618.35
Glutathione9.36633.7923.682.743.691.0018.35
2-Cyanoethylene Oxide
EGCG20.01612.516.302.043.301.4320.65
EGC17.97634.525.882.073.371.2920.65
(+)-Catechin13.95644.963.852.053.450.9920.65
Glutathione13.99712.159.042.633.681.0020.65
AFB1 Exo-8.9-Epoxide
EGCG17.15156.0514.122.403.811.3316.99
EGC13.15204.1410.802.683.401.0216.99
(+)-Catechin14.37143.4910.842.393.321.1116.99
Glutathione12.92245.548.883.063.771.0016.99
a Activation-free energies obtained with the DFT functional M11-L in combination with the flexible 6-311++G(d,p) basis set and SMD implicit solvation model. b The exact imaginary vibrational frequency of the transition state structure. c The lowest vibrational frequency of the reactant state structure. d The distance in the transition state structure between the most nucleophilic phenolic oxygen of EGCG, EGC, (+)-catechin, or the sulfur atom of glutathione, and the achiral electrophilic epoxy carbon of the chemical carcinogen. e The distance in the reactant structure between the most nucleophilic phenolic oxygen of EGCG, EGC, (+)-catechin, or the sulfur atom of glutathione, and the achiral electrophilic epoxy carbon of the chemical carcinogen. f The ratio of activation-free energies between the investigated polyphenols, namely EGCG, EGC, (+)-catechin, and glutathione, allowing for the determination of the relative CSA. g Activation-free energy for the reaction with guanine was calculated using the DFT functional M11-L in combination with the flexible 6-311++G(d,p) basis set and SMD implicit solvation model.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Furlan, V.; Tošović, J.; Bren, U. QM-CSA: A Novel Quantum Mechanics-Based Protocol for Evaluation of the Carcinogen-Scavenging Activity of Polyphenolic Compounds. Foods 2024, 13, 2708. https://doi.org/10.3390/foods13172708

AMA Style

Furlan V, Tošović J, Bren U. QM-CSA: A Novel Quantum Mechanics-Based Protocol for Evaluation of the Carcinogen-Scavenging Activity of Polyphenolic Compounds. Foods. 2024; 13(17):2708. https://doi.org/10.3390/foods13172708

Chicago/Turabian Style

Furlan, Veronika, Jelena Tošović, and Urban Bren. 2024. "QM-CSA: A Novel Quantum Mechanics-Based Protocol for Evaluation of the Carcinogen-Scavenging Activity of Polyphenolic Compounds" Foods 13, no. 17: 2708. https://doi.org/10.3390/foods13172708

APA Style

Furlan, V., Tošović, J., & Bren, U. (2024). QM-CSA: A Novel Quantum Mechanics-Based Protocol for Evaluation of the Carcinogen-Scavenging Activity of Polyphenolic Compounds. Foods, 13(17), 2708. https://doi.org/10.3390/foods13172708

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop