Next Article in Journal
New Advances in the Exploration of Esterases with PET and Fluorescent Probes
Next Article in Special Issue
On the Jahn–Teller Effect in Silver Complexes of Dimethyl Amino Phenyl Substituted Phthalocyanine
Previous Article in Journal
Functionalization of ZnAl-Layered Double Hydroxide with Ensulizole and Its Application as a UV-Protective Agent in a Transparent Polymer Coating
Previous Article in Special Issue
The Synthesis and Base-Induced Breakdown of Triaryl 1,4-Oxathiins—An Experimental and DFT Study
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On the Potential Role of the (Pseudo-) Jahn–Teller Effect in the Membrane Transport Processes: Enniatin B and Beauvericin

Institute of Physical Chemistry and Chemical Physics, Slovak University of Technology in Bratislava (STU), Radlinskeho 9, SK-81237 Bratislava, Slovakia
*
Author to whom correspondence should be addressed.
Molecules 2023, 28(17), 6264; https://doi.org/10.3390/molecules28176264
Submission received: 9 August 2023 / Revised: 20 August 2023 / Accepted: 24 August 2023 / Published: 27 August 2023
(This article belongs to the Special Issue Multiconfigurational and DFT Methods Applied to Chemical Systems)

Abstract

:
The molecular structure of mycotoxins enniatin B and beauvericin, which are used as ionophores, was studied using density functional theory in various symmetry groups and singly charged states. We have shown that the charge addition or removal causes significant structural changes. Unlike the neutral C3 molecules, the stability of the charged C1 structures was explained by the Jahn–Teller or Pseudo-Jahn–Teller effect. This finding agrees with the available experimental X-ray structures of their metal complexes where electron density transfer from the metal can be expected. Hence, the membrane permeability of metal sandwich-structure complexes possessing antimicrobial activities is modulated by the conformational changes.

1. Introduction

Symmetry, a common feature found in nature and chemistry, seems to bypass compounds produced by living systems. In general, high-symmetric organic molecules are toxic and highly carcinogenic. However, there are few types of biomolecules in plants and fungi that possess symmetry, such as the plant pigment curcumine (diferuloylmethane) with antioxidant, antibacterial, anti-inflammatory, and antitumor activity [1,2]. Among the small number of symmetric molecules produced by fungi, those of the Fusarium genus attract extensive attention due to their frequent contamination of food commodities, e.g., cereal grains [3,4]. They produce a number of secondary metabolites, which are questioned due to their potential health hazards [5,6]. Among them, enniatin B (ENB) and beauvericin (BEA) belong to the most studied mycotoxins due to their prevalent occurrence in cereals [7,8,9]. Although toxic, they have shown numerous biological activities, including antibiotic, antimicrobial, insecticidal, anthelmintic, anticancer and apoptosis-inducing effects [10,11,12,13,14,15,16,17,18,19,20].
Their molecular structure consists of three D-hydroxyisovaleric acids in alternation with three N-methyl branched-chain L-amino acids. ENB contains three N-methyl valine residues, whereas in BEA, there are three N-methyl phenylalanine moieties (see Figure 1). The interior of their structures contains strategically spaced oxygens with free electron pairs capable of forming ion–dipole interactions with mono- and divalent cations. The lipophilic exterior of their molecules enables them to incorporate into biological membranes and thus create cation-selective pores. As a result, disturbances of the action potential, metabolic state, and cell homeostasis that lead to cell death may occur [21,22]. Due to the size of small depsipeptide ring sizes, their six carbonyl oxygen atoms must be alternately arranged above and below the mean molecular plane, while the peptide and ester groups are in trans-positions and approximately planar [23]. Further studies revealed that the biological action of these compounds is related to their ability to induce transmembrane cation transport, i.e., act as ionophores [24,25]. Related experiments showed [26] that the membrane effect is associated with their ability to form stable complexes with the cations of alkaline and alkaline-earth metals. As experimentally shown by Ovchinnikov et al. [26], the complexes formed undergo conformational changes depending on the solvent used and on the central atom metal. The flexible depsipeptide chain of enniatins can be adjusted to the size of the metal ion and causes the low selectivity of these complexones. The 2:1 and 3:2 complexes are formed where the sandwiched cations on the symmetry axis are coordinated prevailingly by the carbonyl oxygen atoms. It implies better screening of the metal cation from the anion and solvent in comparison with the 1:1 complexes and higher solubility in organic solvents.
The schematic structure of ENB and BEA depicted in Figure 1 reveals that the highest possible symmetry point group of their molecules is C3. This fact opens the question of whether the Jahn–Teller (JT) or Pseudo-Jahn–Teller (PJT) effect can particularly modulate the energy levels of these symmetric compounds, and which consequences can be awaited from this effect. For the symmetric neutral enol and keto forms of curcumine, we have recently demonstrated the potential significance of the PJT effect [27], which may rationalize the photoprotection action and activity of naturally occurring symmetric dyes in plants.
Although many published works deal with the X-ray structure of naturally occurring ENB and BEA, which are collected in the Cambridge Structural Database (CSD), analysis of their symmetry space group is not available. Therefore, the partial aims of this study are: (1) to analyze the symmetry of ENB and BEA structures available in CSD; (2) to perform optimization of the gas-phase structures in various symmetric groups and charged states using Density functional theory; and (3) to explain the stable structures of lower symmetry using (P)JT treatment. Finally, the possible (P)JT effect for the investigated molecules will be estimated.

2. Theoretical Background

The JT theorem [28] states that any non-linear configuration of atomic nuclei in a degenerate electron state is unstable. Therefore, at least one stable nuclear configuration of lower symmetry must be obtained during a symmetry decrease where the electron degeneracy is removed. The adiabatic potential surfaces (APSs) of such systems can be described by an analytic function based on perturbation theory. Its minimization produces the atomic coordinates of the corresponding stable systems [29]. For large systems, this treatment is too complicated. This problem can be simplified using group theory, and the symmetry of the stable structures can be predicted.
In the method of step-by-step symmetry descent [30,31] the removal of some symmetry elements from the parent ‘unstable’ high-symmetry group (i.e., symmetry decrease) with a multidimensional irreducible representation, which corresponds to a molecule in a degenerate electron state, causes consecutive splitting of this representation into its ‘stable’ subgroups until a nondegenerate electron state (i.e., one-dimensional representation) is obtained. In our case, molecules of the C3 symmetry group in a double-degenerate E electron state undergo a symmetry descent to the C1 symmetry (the C3 symmetry axis is removed).
Alternatively, the epikernel-principle method [32] is based on the JT active distortion coordinate Q of Λ representation for a degenerate electron state Ψ of Γ representation within the parent symmetry group. Λ is the non-totally symmetric part of the symmetrized direct product [Γ ⊗ Γ], which corresponds to a non-vanishing value of < Ψ   H ^ Q Ψ > integrals where H ^ denotes Hamiltonian. For the molecules of the C3 symmetry group in a double-degenerate E electron state:
[E ⊗ E] = A ⊕ E
Here, we obtain the JT active coordinate of the E representation. According to the epikernel principle, the extrema of a JT energy surface correspond to the kernel K(G, Λ) or epikernel E(G, Λ) subgroups of the parent group G. Kernels contain symmetry operations that leave the Λ representation invariant, whereas epikernels leave invariant-only some components of the degenerate Λ representation. In our case,
K(C3, E) = E(C3, E) = C1
is in agreement with the above-mentioned method of step-by-step symmetry descent. The energy difference between the high-symmetry unstable and low-symmetry stable structures of the same compound is denoted as the Jahn–Teller stabilization energy EJT.
A similar instability of a high-symmetry structure, known as the PJT effect [33], may be observed in the case of sufficiently strong vibronic coupling between two pseudodegenerate electronic states (usually ground and excited). Using perturbation theory, the energy surface for the simplest case of two interacting electronic states Ψ1 and Ψ2 of different space symmetries can be described by the formula
E Q = 1 2 K Q 2 ± [ 2 4 + F 2 Q 2 ] 1 / 2
where E is the energy of the electronic state, Q is the distortion coordinate, Δ is the difference in energy between both electronic states in the undistorted geometry, K is the primary force constant (without vibronic coupling) and F is the vibronic coupling constant. The curvature of the lower state is negative if the energy difference is
< 2 F 2 K
Therefore, an instability in the Q direction can be concluded. Otherwise, the stable structure corresponds to Q = 0, that is, the high-symmetry structure is preserved despite the diminished curvature of the lower energy state. Higher Δ values imply a weaker PJT interaction. Vibronic interactions between the states of various spin multiplicities are forbidden because their spin states are orthogonal. In general, PJT interactions are possible in any system. Only if their consequences are observed (in spectra, structure, reactivity, etc.) can the corresponding vibronic interaction be denoted as the PJT effect.
The symmetry of stable PJT structures can be predicted using the epikernel principle method [34,35]. For non-vanishing values of < Ψ 1   H ^ Q Ψ 2 > integrals, the representation Λ of the JT active coordinate Q is corresponding to the non-totally symmetric part of the direct product of Γ1 and Γ2 representations of electronic states Ψ1 and Ψ2, respectively. The symmetry groups of stable structures correspond to the kernel K(G, Λ) or the epikernel E(G, Λ) subgroups of the parent group G (see above). For the molecules of the C3 symmetry group in a non-degenerate A ground electron state, excited E states can produce non symmetric JT active coordinates because of the direct product
A ⊗ E = E
According to Equation (2), we again obtain stable structures of the C1 symmetry group.

3. Results

In the Cambridge Structural Database (CSD) [36], we have found five X-ray structures that contain a neutral ENB molecule (Table 1). Moreover, three additional X-ray structures of enniatin B complexes with alkaline metals were acquired as well. We checked the molecular symmetry using the distance between neighboring carbonyl oxygens OC=O (which might be bonded to metal cations), as well as bridging oxygens Obridge and nitrogens N and confirmed the results by the (O–O–O)C=O, (O–O–O)bridge and N–N–N angles between the atoms mentioned (Figure 2). Despite too high R-factors (and missing standard deviations) of some structures, we confirmed the C3 symmetry axis in neutral ENB molecules in the BICMEF, DESYIJ, and EROPIM structures, as well as in its potassium complex in the IHECUT structure. Fewer results were obtained for beauvericin. Only one X-ray structure containing its neutral molecule and a single beauvericin complex with barium was found in CSD (Table 1). Using the treatment mentioned above, there is no C3 symmetry axis in these structures (Table 2). This indicates the symmetry descent, which may be ascribed to vibronic interactions, i.e., some sort of (P)JT effect.
In the first step of our model study, we tried to optimize all systems within the C3 symmetry group. However, we have not found any C3 structures of cationic ENB and charged BEA species because they adopt a double-degenerate ground electron state that is not accessible within DFT treatment. Therefore, the JT effect reduces their molecular structures to C1 symmetry, as illustrated in Table 3. Neutral ENB and BEA molecules are the most stable, whereas their cations are the least stable. The calculated energies of their anions are close to those of the neutral ones. This reflects their affinity to form metal complexes (see Table 2). We have found anionic ENB structures of C3 symmetry, which are unstable (imaginary vibration of e representation), and the stable ones of C1 symmetry. Their existence might be explained by the PJT effect (see below).
The data collected in Tables S1 and S2 in Supplementary Information reveal that the electronic excited states of the systems under study are too high. The lowest energy vertical electron energy transition (S0→S1) is 5.368 eV (231 nm) for ENB and 5.363 eV (231 nm) for BEA. These energy values correspond to the UV region in agreement with experimental observations [26,37]. Vibronic interactions with double-degenerate excited states of much lower energies (Table S1 in Supplementary Information) produce an imaginary vibration (i.e., PJT active coordinate) of e symmetry that leads to the stable structure of the C1 symmetry group with EJT of 0.61 eV. However, the correlation of excited states of anionic ENB between the C3 and C1 symmetries is problematic. The electronic excitation energies in the stable C1 structure are much higher, and the corresponding oscillator strengths are very low. If we assume that small PJT perturbation does not cause significant changes in oscillator strengths, the 12E electron state in the C3 structure, with very high oscillator strength, has no counterparts among low excited states in the C1 structure. This implies that either this assumption is incorrect (i.e., the PJT effect does not correspond to a small perturbation), or that the 12E electron state does not undergo the above vibronic interaction.
We checked the molecular symmetry of all ENB and BEA species by interatomic distances and angles between heteroatoms (Table 4), as well as by their natural charges and spin populations (Table 5). The mutual distances between neighboring carbonyl O atoms in neutral molecules are shorter than those in cationic species and longer than those in anionic species. The reverse trend holds for the distances between the remaining heteroatoms. This implies a higher C=O perpendicularity to the main molecular plane in anionic species. A large benzyl substituent in BEA repels C=O groups more toward the center of the depsipeptide ring than the isopropyls in ENB.
Bridging O atoms in our systems are less negative than the carbonyl ones (Table 5). As expected, the negative charges of heteroatoms are higher in the anionic species, whereas in the cationic ones, they are lower than in the neutral molecules. The N atoms in the cationic C1 species have the highest spin populations. The spin density distribution in charged species is highly non-symmetric (Figure 3 and Figure 4). In the beauvericin anion, the spin density is located mainly at aromatic rings, whereas in remaining systems it is prevailingly located in the area of a single N atom and neighboring carbonyl. This asymmetric distribution should increase the reactivity of charged species.

4. Method

The geometries of the neutral (charge q = 0), cationic (q = +1) and anionic (q = −1) molecules of ENB and BEA (Figure 2) in the spin states of the ground singlet (q = 0) or doublet (q = +1, −1) were optimized within the C3 and C1 symmetry point groups using the M06-2X hybrid functional [38] (a very recent functional with dispersion corrections), combined with standard cc-pVDZ basis sets for all atoms taken from the Gaussian library (we were not able to use larger basis sets because of technical restrictions). Stability of the obtained geometries was checked on imaginary vibrations by vibrational analysis. The Gibbs free energies were computed for room temperature (298.15 K, see Supplementary Materials for details). Natural charges and natural spin populations on atoms were evaluated using Natural Bond Orbitals (NBO) analysis [39,40], which are more reliable than the analogous results obtained by Mulliken population analysis. Time-dependent DFT (TD-DFT) treatment [41,42] (more frequently used than single-excitation Configuration Interaction methods) for up to 30 vertical electronic states was used for excited state calculations. Gaussian16 (Revision B.01) software [43] was used for all quantum-chemical calculations because it brings more new methods, property predictions and performance enhancements than similar quantum-chemical software. The MOLDRAW (https://www.moldraw.software.informer.com, accessed on 9 September 2019) [44] and Molekel (Version 5.4.0.8) [45] software were used for geometry manipulation and visualization purposes.

5. Conclusions

We have shown that the charged species of compounds under study exhibit JT symmetry descent, unlike neutral molecules. It agrees with the experimental X-ray structures of their metal complexes (Table 2), where electron density transfer from the metal can be expected. Analysis of the symmetry changes reveals the existence of three stable C1 geometries as a consequence of the JT or PJT effect in their C3 parent structures. This fluxional behavior might explain the experimentally observed equilibria between various forms of metal complexes [26]. Our study indicates that the membrane permeability of the sandwich structure complexes mentioned above is modulated by conformational changes between various forms, which can be explained by the (P)JT effect. Further theoretical study of alkaline and alkaline-earth metal complexes is desirable to verify this hypothesis. Our approach can be used for any symmetric configuration of atomic nuclei undergoing the (P)JT effect.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/molecules28176264/s1, Tables S1 and S2: Lowest excited electronic states, corresponding excitation energies (Eexc) and oscillator strengths (f) in atomic units for the optimized structures of enniatin B and beauvericin, respectively, in various charge states (q), Gibbs free energy in Gaussian. Ref. [46] is cited in the Supplementary Materials.

Author Contributions

Investigation, D.Š., V.L. and M.B.; writing—original draft preparation, D.Š., V.L. and M.B.; writing—review and editing, V.L. and M.B. All authors have read and agreed to the published version of the manuscript.

Funding

The work has been supported by the Slovak Research and Development Agency (nos. APVV-19-0024 and APVV-20-0213), by the Scientific Grant Agency of the Ministry of Education, Science, Research, and Sport of the Slovak Republic VEGA (nos. 1/0461/21 and 1/0139/20), and by the Operational Program Integrated Infrastructure for the project: “ Support of research activities of Excellence laboratories STU in Bratislava”, Project no. 313021BXZ1, co-financed by the European Regional Development Fund.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data is contained within the article or Supplementary Materials.

Acknowledgments

The authors thank the HPC center at the Slovak University of Technology in Bratislava, which is a part of the Slovak Infrastructure of High Performance Computing (SIVVP project ITMS 26230120002 funded by the European Region Development Funds) for the computational time and resources made available.

Conflicts of Interest

The authors declare no conflict of interest.

Sample Availability

Not applicable.

References

  1. Sun, Y.M.; Zhang, H.Y.; Chen, D.Z.; Liu, C.B. Theoretical Elucidation on the Antioxidant Mechanism of Curcumin:  A DFT Study. Org. Lett. 2002, 4, 2909–2911. [Google Scholar] [PubMed]
  2. Aggarwal, B.B.; Sung, B. Pharmacological basis for the role of curcumin in chronic diseases: An age-old spice with modern targets. Trends Pharmacol. 2009, 30, 85–94. [Google Scholar]
  3. Logrieco, A.; Moretti, A.; Mule, G.; Paciolla, C.; Ritieni, A. Advances on the toxicity of the cereal contaminant Fusarium esadepsipeptides. Cereal Res. Commun. 2008, 36, 303–313. [Google Scholar]
  4. Ezekiel, C.N.; Ayeni, K.I.; Akinyemi, M.O.; Sulyok, M.; Oyedele, O.A.; Babalola, D.A.; Ogara, I.M.; Krska, R. Dietary Risk Assessment and Consumer Awareness of Mycotoxins among Household Consumers of Cereals, Nuts and Legumes in North-Central Nigeria. Toxins 2021, 13, 635. [Google Scholar] [PubMed]
  5. Jestoi, M. Emerging fusarium-mycotoxins fusaproliferin, beauvericin, enniatins, and moniliformin: A review. Crit. Rev. Food Sci. Nutr. 2008, 48, 21–49. [Google Scholar]
  6. Firáková, S.; Proksa, B.; Sturdíkova, M. Biosynthesis and biological activity of enniatins. Pharmazie 2007, 62, 563–568. [Google Scholar] [PubMed]
  7. Ivanova, L.; Egge-Jacobsen, W.M.; Solhaug, A.; Thoen, E.; Faeste, C.K. Lysosomes as a possible target of enniatin B-induced toxicity in Caco-2 cells. Chem. Res. Toxicol. 2012, 25, 1662–1674. [Google Scholar] [PubMed]
  8. Yuan, Y.; Meng, G.; Li, Y.; Wu, C. Study on In Vitro Metabolism and In Vivo Pharmacokinetics of Beauvericin. Toxins 2022, 14, 477. [Google Scholar]
  9. Logrieco, A.; Moretti, A.; Castella, G.; Kostecki, M.; Golinski, P.; Ritieni, A.; Chelkowski, J. Beauvericin production by Fusarium species. Appl. Environ. Microb. 1998, 64, 3084–3088. [Google Scholar]
  10. Jestoi, M.; Rokka, M.; Yli-Mattila, T.; Parikka, P.; Rizzo, A.; Peltonen, K. Presence and concentrations of the Fusarium-related mycotoxins beauvericin, enniatins and moniliformin in finnish grain species. Food Addit. Contamin. 2004, 21, 794–802. [Google Scholar]
  11. Meca, G.; Sospedra, I.; Valero, M.A.; Mañes, J.; Font, G.; Ruiz, M.J. Antibacterial activity of the enniatin B, produced by Fusarium tricinctum in liquid culture, and cytotoxic effects on Caco-2 cells. Toxicol. Mech. Methods 2011, 21, 503–512. [Google Scholar]
  12. Lin, H.I.; Lee, Y.J.; Chen, B.F.; Tsai, M.C.; Lu, J.L.; Chou, C.J.; Jow, G.M. Involvement of Bcl-2 family, cytochrome c and caspase 3 in induction of apoptosis by beauvericin in human non-small cell lung cancer cells. Cancer Lett. 2005, 239, 248–259. [Google Scholar]
  13. Jow, G.M.; Chou, C.J.; Chen, B.F.; Tsai, J.H. Beauvericin induces cytotoxic effects in human acute lymphoblastic leukemia cells through cytochrome c release, caspase 3 activation: The causative role of calcium. Cancer Lett. 2004, 216, 165–173. [Google Scholar] [PubMed]
  14. Dornetshuber-Fleiss, R.; Heilos, D.; Mohr, T.; Richter, L.; Süssmuth, R.D.; Zlesak, M.; Novicky, A.; Heffeter, P.; Lemmens-Gruber, R.; Berger, W. The naturally born fusariotoxin enniatin B and sorafenib exert synergistic activity against cervical cancer in vitro and in vivo. Biochem. Pharmacol. 2015, 93, 318–331. [Google Scholar]
  15. Dornetshuber, R.; Heffeter, P.; Kamyar, M.R.; Peterbauer, T.; Berger, W.; Lemmens-Gruber, R. Enniatin exerts p53-dependent cytostatic and p53-independent cytotoxic activities against human cancer cells. Chem. Res. Toxicol. 2007, 20, 465–473. [Google Scholar] [PubMed]
  16. Wu, Q.; Patocka, J.; Nepovimova, E.; Kuca, K. A review on the synthesis and bioactivity aspects of beauvericin, a Fusarium mycotoxin. Front. Pharmacol. 2018, 9, 1338. [Google Scholar]
  17. Sood, S.; Sandhu, S.S.; Mukherjee, T.K. Pharmacological and Therapeutic Potential of Beauvericin: A Short Review. J. Proteom. Bioinform. 2017, 10, 18–23. [Google Scholar]
  18. Castlebury, L.A.; Sutherland, J.B.; Tanner, L.A.; Henderson, A.L.; Cerniglia, C.E. Use of a bioassay to evaluate the toxicity of beauvericin to bacteria. World J. Microb. Biot. 1999, 15, 119–121. [Google Scholar] [CrossRef]
  19. Zhang, L.; Yan, K.; Zhang, Y.; Huang, R.; Bian, J.; Zheng, C.; Sun, H.; Chen, Z.; Sun, N.; An, R.; et al. High-throughput synergy screening identifies microbial metabolites as combination agents for the treatment of fungal infections. Proc. Natl. Acad. Sci. USA 2007, 104, 4606–4611. [Google Scholar]
  20. Fukuda, T.; Arai, M.; Yamaguchi, Y.; Masuma, R.; Tomoda, H.; Omura, S. New beauvericins, potentiators of antifungal miconazole activity, Produced by Beauveria sp. FKI-1366. I. Taxonomy, fermentation, isolation and biological properties. J. Antibiot. 2004, 57, 110–116. [Google Scholar]
  21. Kouri, K.; Lemmens, M.; Lemmens-Gruber, R. Beauvericininduced channels in ventricular myocytes and liposomes. Biochim. Biophys. Acta 2003, 1609, 203–210. [Google Scholar] [PubMed]
  22. Kamyar, M.; Rawnduzi, P.; Studenik, C.R.; Kouri, K.; Lemmens-Gruber, R. Investigation of the electrophysiological properties of enniatins. Arch. Biochem. Biophys. 2004, 429, 215–223. [Google Scholar] [PubMed]
  23. Dobler, M.; Dunitz, I.D.; Krajewski, J. Structure of the K+ complex with enniatin B, a macrocyclic antibiotic with K+ transport properties. J. Mol. Biol. 1969, 42, 603. [Google Scholar] [PubMed]
  24. Mueller, P.; Rudin, D.O. Development of K+-Na+ discrimination in experimental bimolecular lipid membranes by macrocyclic antibiotics. Biochem. Biophys. Res. Commun. 1967, 26, 398404. [Google Scholar]
  25. Pressman, B.C. Biological applications of ionophores. Annu. Rev. Biochem. 1976, 45, 501–530. [Google Scholar] [CrossRef] [PubMed]
  26. Ovchinnikov, Y.A.; Ivanov, V.T.; Evstratov, A.V.; Mikhaleva, I.I.; Bystrov, V.F.; Portnova, S.L.; Balashova, T.A.; Meshcheryakova, E.N.; Tulchinsky, V.M. The Enniatin Ionophoress. Conformation and ion binding properties. J. Pept. Protein Res. 1974, 6, 465–498. [Google Scholar]
  27. Štellerová, D.; Lukeš, V.; Breza, M. How does pseudo-Jahn-Teller effect induce the photoprotective potential of curcumin? Molecules 2023, 28, 2946. [Google Scholar]
  28. Jahn, H.A.; Teller, E. Stability of polyatomic molecules in degenerate electronic states. I. Orbital degeneracy. Proc. R. Soc. Lond. A 1937, 161, 220–235. [Google Scholar]
  29. Opik, U.; Pryce, M.H.L. Studies of the Jahn-Teller effect I. A survey of the static problem. Proc. R. Soc. Lond. A 1957, 238, 425–447. [Google Scholar]
  30. Pelikán, P.; Breza, M. Classification of the possible symmetries of the Jahn—Teller systems. Chem. Pap. 1984, 39, 255–270. [Google Scholar]
  31. Breza, M. Group-Theoretical Analysis of Jahn-Teller Systems. In The Jahn-Teller Effect. Fundamentals and Implications for Physics and Chemistry; Köppel, H., Yarkony, D.R., Barentzen, H., Eds.; Springer: Berlin/Heidelberg, Germany, 2009; pp. 51–76. ISBN 978-3642034312. [Google Scholar]
  32. Ceulemans, A.; Beyens, D.; Vanquickenborne, L.G. Symmetry aspects of Jahn-Teller activity: Structure and reactivity. J. Am. Chem. Soc. 1984, 106, 5824–5837. [Google Scholar] [CrossRef]
  33. Bersuker, I.B. Pseudo-Jahn–Teller Effect—A Two-State Paradigm in Formation, Deformation, and Transformation of Molecular Systems and Solids. Chem. Rev. 2013, 113, 1351–1390. [Google Scholar] [CrossRef] [PubMed]
  34. Ceulemans, A.; Vanquickenborne, L.G. The Epikernel Principle. Struct. Bond. 1989, 71, 125–159. [Google Scholar]
  35. Breza, M. Group-Theoretical Treatment of Pseudo-Jahn-Teller Systems. In Vibronic Interactions and the Jahn-Teller Effect: Theory and Application. (Prog. Theor. Chem. Phys. B 23); Atanasov, M., Daul, C., Tregenna-Piggott, P.L.W., Eds.; Springer: Dordrecht, The Netherlands; Berlin/Heidelberg, Germany; London, UK; New York, NY, USA, 2012; pp. 59–82. ISBN 1567-7354. [Google Scholar]
  36. Groom, C.R.; Bruno, I.J.; Lightfoot, M.P.; Ward, S.C. The Cambridge Structural Database. Acta Cryst. 2016, B72, 171–179. [Google Scholar] [CrossRef] [PubMed]
  37. Xu, L.; Wang, J.; Zhao, J.; Li, P.; Shan, T.; Wang, J.; Li, X.; Zhou, L. Beauvericin from the Endophytic Fungus, Fusarium redolens, Isolated from Dioscorea zingiberensis and Its Antibacterial Activity. Nat. Prod. Commun. 2010, 5, 811–814. [Google Scholar] [CrossRef] [PubMed]
  38. 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 function. Theor. Chem. Acc. 2008, 120, 215–241. [Google Scholar] [CrossRef]
  39. Carpenter, J.E.; Weinhold, F. Analysis of the geometry of the hydroxymethyl radical by the different hybrids for different spins natural bond orbital procedure. J. Mol. Struct. 1988, 46, 41–62. [Google Scholar] [CrossRef]
  40. Lee, C.; Yang, W.; Parr, R.G. Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density. Phys. Rev. B 1988, 37, 785–789. [Google Scholar] [CrossRef]
  41. Bauernschmitt, R.; Ahlrichs, R. Treatment of electronic excitations within the adiabatic approximation of time dependent density functional theory. Chem. Phys. Lett. 1996, 256, 454–464. [Google Scholar] [CrossRef]
  42. Scalmani, G.; Frisch, M.J.; Mennucci, B.; Tomasi, J.; Cammi, R.; Barone, V. Geometries and properties of excited states in the gas phase and in solution: Theory and application of a time-dependent density functional theory polarizable continuum model. J. Chem. Phys. 2006, 124, 94107. [Google Scholar] [CrossRef]
  43. Frisch, G.W.M.J.; Schlegel, B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Petersson, G.A.; Nakatsuji, H.; Li, X.; et al. Gaussian 16, Revision B.01; Gaussian, Inc.: Wallingford, CT, USA, 2016. [Google Scholar]
  44. Ugliengo, P. MOLDRAW: A Program to Display and Manipulate Molecular and Crystal Structures, University Torino, Torino. 2012. Available online: https://www.moldraw.software.informer.com (accessed on 9 September 2019).
  45. Varetto, U. Molekel, Ver. 5.4.0.8; Swiss National Supercomputing Centre: Lugano, Switzerland, 2009; Available online: https://www.molekel.software.informer.com (accessed on 8 July 2017).
  46. Ochterski, J.W. Thermochemistry in Gaussian; Gaussian, Inc. 2020. Available online: https://gaussian.com/thermo/ (accessed on 18 August 2023).
Figure 1. Schematic structure and atom notation of enniatin B (a) and beauvericin (b).
Figure 1. Schematic structure and atom notation of enniatin B (a) and beauvericin (b).
Molecules 28 06264 g001
Figure 2. M06-2X/cc-pVDZ optimized gas phase structure in C3 symmetry group of ENB (a) and BEA (b) (O—red, N—blue, C—black, H—white).
Figure 2. M06-2X/cc-pVDZ optimized gas phase structure in C3 symmetry group of ENB (a) and BEA (b) (O—red, N—blue, C—black, H—white).
Molecules 28 06264 g002
Figure 3. M06-2X/cc-pVDZ calculated spin density distribution over the cation of (a) enniatin B and (b) beauvericin in C1 symmetry (O—red, N—blue, C—grey, H—white, positive, and negative spin density at the 0.0005 a.u. iso-surface is indicated as dark blue and green, respectively).
Figure 3. M06-2X/cc-pVDZ calculated spin density distribution over the cation of (a) enniatin B and (b) beauvericin in C1 symmetry (O—red, N—blue, C—grey, H—white, positive, and negative spin density at the 0.0005 a.u. iso-surface is indicated as dark blue and green, respectively).
Molecules 28 06264 g003
Figure 4. M06-2X/cc-pVDZ calculated spin density distribution over the anion of (a) enniatin B and (b) beauvericin in C1 symmetry (O—red, N—blue, C—grey, H—white, positive, and negative spin density at the 0.0005 a.u. iso-surface is indicated as dark blue and green, respectively).
Figure 4. M06-2X/cc-pVDZ calculated spin density distribution over the anion of (a) enniatin B and (b) beauvericin in C1 symmetry (O—red, N—blue, C—grey, H—white, positive, and negative spin density at the 0.0005 a.u. iso-surface is indicated as dark blue and green, respectively).
Molecules 28 06264 g004
Table 1. X-ray structures of the molecules under study found in CSD [36].
Table 1. X-ray structures of the molecules under study found in CSD [36].
CSD CodeChemical FormulaChemical NameSpace GroupR-Factor
Neutral ENB molecules
BICMEFC33H57N3O9·1.5H2OCyclo-tris(l-methylvalyl-d-2-hydroxyisovaleryl) sesquihydrateR30.06
CIKJAHC33H57N3O9·1.67H2OEnniatin B hydrateR30.057
DESYIJC33H57N3O9Enniatin B unknown solvateR30.113
EROPIM(C11H19NO3)3·0.2238(CH3O)3· 0.328677O3,6,9,12,15,18-hexaisopropyl-4,10,16-trimethyl-1,7,13-trioxa-4,10,16-triazacyclooctadecane-2,5,8,11,14,17-hexone methanol solvate hydrateR30.0397
ZASQOZC33H57N3O9Enniatin BP11210.0714
ENB complexes
IHECUT(C132H229K4N12O36)n(CNS)4ncatena-(tetrakis((μ2-Enniatin B)-potassium thiocyanate))P30.0643
MVHIRB10(C34H57N4O9RbS)ncatena-(l,d,l,l,d,l-Enniatin B)-rubidium isothiocyanateP432120.085
PEKFEQ(C66H114KN6O18)I·H2Obis(Enniatin B)-potassium iodide monohydrateP610.0766
Neutral BEA molecule
BEVERCC45H57N3O9·H2OCyclo-tri(l-N-methylphenylalanyl-d-alpha-hydroxyisovaleryl) monohydrateP210.096
BEA complex
BEAVBA(C45H57N3O9)(C6H2N3O7)2Ba·2C7H8Beauvericin barium picrate toluene solvateP212120.15
Table 2. Molecular point groups Γ, neighboring carbonyl O–O distances d(O–O)C=O, bridging O–O distances d(O–O)bridge, N–N distances d(N–N), neighboring carbonyl O–O–O angles Θ(O–O–O)C=O, bridging O–O–O angles Θ(O–O–O)bridge, and N–N–N angles Θ(N–N–N) in X-ray structures of ENB and BEA. All distances are in Å and all angles are in degrees.
Table 2. Molecular point groups Γ, neighboring carbonyl O–O distances d(O–O)C=O, bridging O–O distances d(O–O)bridge, N–N distances d(N–N), neighboring carbonyl O–O–O angles Θ(O–O–O)C=O, bridging O–O–O angles Θ(O–O–O)bridge, and N–N–N angles Θ(N–N–N) in X-ray structures of ENB and BEA. All distances are in Å and all angles are in degrees.
SystemΓd(O–O)C=Od(O–O)bridged(N–N)Θ(O–O–O)C=OΘ(O–O–O)bridgeΘ(N–N–N)
Neutral ENB molecules
BICMEFC34.827(3×)5.636(3×)5.497(3×)74.6(3×)60.0(3×)60.0(3×)
4.725(3×) 71.9(3×)
CIKJAHC14.86(1), 4.70(1)5.577(9)5.570(7)74.7(2), 82.6(2)60.0(1) (3×)60.0(1) (3×)
4.31(1), 3.66(1)5.58(1)5.570(9)74.7(2), 81.6(2)
4.86(1), 4.70(1)5.577(9)5.57(1)74.7(2), 73.1(2)
DESYIJC34.01(3) (3×)5.87(2) (3×)5.71(3) (3×)67.5(5) (3×)60.0(3) (3×)60.0(3) (3×)
3.89(3) (3×)66.0(5) (3×)
EROPIMC34.888(2) (3×)5.556(3) (3×)5.444(3) (3×)75.07(4) (3×)60.00(3) (3×)60.00(4) (3×)
4.776(3) (3×)74.57(4) (3×)
ZASQOZ (a)C14.746(9), 3.41(1)5.829(9)5.989(8)82.8(2), 74.2(2)56.38(9)52.8(1)
3.69(1), 4.26(1)5.90(1)5.510(9)57.4(2), 48.2(1)61.2(1)68.4(1)
4.907(9),4.969(9)5.540(7)5.13(1)72.7(2), 83.5(2)62.4(1)58.8(1)
4.92(1), 4.805(9)5.631(8)5.455(9)84.8(2), 74.0(2)59.8(1)51.8(1)
4.10(1), 3.72(1)5.91(1)5.97(1)56.5(2), 73.5(3)57.7(1)58.8(1)
3.64(1), 4.77(1)5.762(8)5.01(1)86.1(3), 68.1(2)62.5(1)69.4(1)
ENB complexes
IHECUT (b)C33.40(1) (3×)5.87(2) (3×)6.07(2) (3×)68.6(3) (3×)60.0(2) (12×)60.0(2)(12×)
3.37(2) (3×)64.7(3) (3×)
3.42(1) (3×)6.04(1) (3×)5.93(2) (3×)67.4(3) (3×)
3.52(2) (3×)67.3(3) (3×)
3.24(1) (3×)6.27(2) (3×)6.02(2) (3×)67.9(4) (3×)
3.59(2) (3×)64.6(4) (3×)
3.45(1) (3×)5.80(2) (3×)5.86(2) (3×)66.7(3) (3×)
3.46(1) (3×)63.2(3) (3×)
MVHIRB10C13.26, 3.525.855.8454.3, 91.456.865.5
3.70, 3.285.805.4690.9, 66.865.758.3
3.63, 4.246.325.3263.3, 123.657.556.1
PEKFEQ (a)C13.64(5), 3.64(4)5.93(4)5.93(4)68.5(7), 65.9(7)60.8(4)61.1(5)
3.55(3), 3.80(3)6.01(3)5.82(3)68.0(7), 62.9(6)58.9(4)59.6(5)
3.67(4), 3.63(4)5.94(3)6.05(5)70.5(7), 62.2(7)60.3(4)59.3(5)
3.62(4), 3.55(4)5.90(2)5.99(3)77.3(9), 61.6(9)60.8(4)58.8(4)
3.69(3), 3.48(3)6.02(4)6.20(5)77.7(9), 64.6(7)59.5(4)62.3(5)
3.76(4), 3.62(3)5.99(3)5.98(3)68.5(7), 56.0(6)59.7(4)58.9(4)
Neutral BEA molecule
BEVERCC13.21, 4.036.245.7052.5, 78.958.064.2
3.55, 4.066.305.9462.4, 70.760.556.0
3.95, 3.496.086.1968.5, 73.961.659.8
BEA complex
BEAVBAC13.33, 3.466.335.9267.6, 69.461.259.6
3.36, 3.826.205.9063.9, 74.259.160.0
3.11, 3.496.235.8869.2, 66.459.760.4
Remarks: (a) two independent ENB molecules; (b) four independent ENB molecules.
Table 3. Electronic DFT energies (EDFT), Gibbs free energies at room temperature (G298), Jahn–Teller stabilization energies (EJT), and symmetries of imaginary vibrations for optimized geometries of enniatin B and beauvericin in various charge (q) and spin states.
Table 3. Electronic DFT energies (EDFT), Gibbs free energies at room temperature (G298), Jahn–Teller stabilization energies (EJT), and symmetries of imaginary vibrations for optimized geometries of enniatin B and beauvericin in various charge (q) and spin states.
CompoundqSymmetry GroupGround Electron StateEDFT [Hartree]G298 [Hartree]EJT [eV]Imaginary Vibration
Symmetry
ENB0C31A−2132.49018−2131.682300.000-
+1C12A−2132.19790−2131.39270unknown-
−1C32A−2132.45687−2131.65855-e
C12A−2132.47930−2131.674060.610-
BEA0C31A−2589.67257−2588.792820.000-
+1C12A−2589.38017−2588.50454unknown-
−1C12A−2589.62864−2588.76092unknown-
Table 4. Molecular point groups, Γ, neighboring carbonyl O–O distances, d(O–O)C=O, bridging O–O distances, d(O–O)bridge, N–N distances, d(N–N), neighboring carbonyl O–O–O angles, Θ(O–O–O)C=O, bridging O–O–O angles, Θ(O–O–O)bridge, and N–N–N angles, Θ(N–N–N) of DFT optimized structures of enniatin B and beauvericin in various charge states q. All distances are in Å, all angles are in degrees.
Table 4. Molecular point groups, Γ, neighboring carbonyl O–O distances, d(O–O)C=O, bridging O–O distances, d(O–O)bridge, N–N distances, d(N–N), neighboring carbonyl O–O–O angles, Θ(O–O–O)C=O, bridging O–O–O angles, Θ(O–O–O)bridge, and N–N–N angles, Θ(N–N–N) of DFT optimized structures of enniatin B and beauvericin in various charge states q. All distances are in Å, all angles are in degrees.
CompoundqΓd(O–O)C=Od(O–O)bridged(N–N)Θ(O–O–O)C=OΘ(O–O–O)bridgeΘ(N–N–N)
ENB0C34.062 (3×)6.037 (3×)5.721 (3×)67.3 (3×)60.0 (3×)60.0 (3×)
3.845 (3×)65.7 (3×)
+1C13.380, 3.2056.0246.11468.1, 65.963.457.7
3.234, 3.0096.3935.77580.4, 69.357.758.8
3.742, 3.1226.5335.84656.7, 64.461.163.5
−1C35.026 (3×)5.389 (3×)5.224 (3×)78.0 (3×)60.0 (3×)60.0 (3×)
4.943 (3×)79.0 (3×)
C15.040, 5.0275.2615.18476.1, 80.466.759.3
4.983, 5.0654.9035.01788.8, 78.854.562.6
5.140, 5.0434.6654.95482.7, 84.858.858.1
BEA0C33.194 (3×)6.357 (3×)5.920 (3×)73.8 (3×)60.0 (3×)60.0 (3×)
3.593 (3×)60.9 (3×)
+1C13.122, 3.3806.3926.11469.3, 80.463.457.7
3.205, 3.2346.5335.77565.9, 68.161.163.5
3.009, 3.7426.0245.84664.4, 56.755.558.8
−1C13.211, 3.6416.3805.87762.6, 73.060.059.9
3.192, 3.6076.4095.87861.8, 73.160.459.9
3.184, 3.6086.3605.89961.4, 71.059.660.2
Table 5. Molecular point groups, Γ, natural charges of carbonyl oxygens, Q (O)C=O, bridging oxygens, Q(O)bridge, and nitrogens, Q(N), natural spin populations of carbonyl oxygens, ρ(O)C=O, bridging oxygens, ρ(O)bridge, and nitrogens, ρ(N) of DFT optimized structures of enniatin B and beauvericin in various charged states q.
Table 5. Molecular point groups, Γ, natural charges of carbonyl oxygens, Q (O)C=O, bridging oxygens, Q(O)bridge, and nitrogens, Q(N), natural spin populations of carbonyl oxygens, ρ(O)C=O, bridging oxygens, ρ(O)bridge, and nitrogens, ρ(N) of DFT optimized structures of enniatin B and beauvericin in various charged states q.
CompoundqΓQ(O)C=OQ(O)bridgeQ(N)ρ(O)C=Oρ(O)bridgeρ(N)
ENB0C3−0.639 (3×)−0.585 (3×)−0.554 (3×)---
−0.632 (3×)
+1C1−0.580, −0.454−0.554−0.2320.078, 0.1160.0140.657
−0.647, −0.623−0.576−0.5660.005, 0.0010.0000.000
−0.639, −0.617−0.583−0.5510.000, 0.0010.0000.001
−1C3−0.704 (3×)−0.594 (3×)−0.544 (3×)0.068 (3×)0.006 (3×)0.003 (3×)
−0.681 (3×)0.029 (3×)
C1−0.634, −0.684−0.591−0.5490.000 (4×)
0.265, 0.003
0.0000.001
−0.824, −0.688−0.669−0.5300.0310.036
−0.642, −0.651−0.579−0.5460.0000.000
BEA0C3−0.607 (3×)−0.598 (3×)−0.560 (3×)---
−0.638 (3×)
+1C1−0.560, −0.665−0.580−0.2130.040, 0.0030.0190.665
0.000 (2×)
−0.612, −0.641−0.593−0.5320.000 (2×)0.000
−0.628, −0.442−0.584−0.559−0.002, 0.174−0.001
−1C1−0.619, −0.632−0.600 (2×)−0.5590.007, 0.0030.000 (2×)
0.001
0.001
0.000 (2×)
−0.616, −0.631−0.601−0.5600.004, 0.001
−0.624, −0.634 −0.6010.012, 0.005
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

Štellerová, D.; Lukeš, V.; Breza, M. On the Potential Role of the (Pseudo-) Jahn–Teller Effect in the Membrane Transport Processes: Enniatin B and Beauvericin. Molecules 2023, 28, 6264. https://doi.org/10.3390/molecules28176264

AMA Style

Štellerová D, Lukeš V, Breza M. On the Potential Role of the (Pseudo-) Jahn–Teller Effect in the Membrane Transport Processes: Enniatin B and Beauvericin. Molecules. 2023; 28(17):6264. https://doi.org/10.3390/molecules28176264

Chicago/Turabian Style

Štellerová, Dagmar, Vladimír Lukeš, and Martin Breza. 2023. "On the Potential Role of the (Pseudo-) Jahn–Teller Effect in the Membrane Transport Processes: Enniatin B and Beauvericin" Molecules 28, no. 17: 6264. https://doi.org/10.3390/molecules28176264

APA Style

Štellerová, D., Lukeš, V., & Breza, M. (2023). On the Potential Role of the (Pseudo-) Jahn–Teller Effect in the Membrane Transport Processes: Enniatin B and Beauvericin. Molecules, 28(17), 6264. https://doi.org/10.3390/molecules28176264

Article Metrics

Back to TopTop