DFT Modelling of Molecular Structure, Vibrational and UV-Vis Absorption Spectra of T-2 Toxin and 3-Deacetylcalonectrin

This paper discusses the applicability of optical and vibrational spectroscopies for the identification and characterization of the T-2 mycotoxin. Vibrational states and electronic structure of the T-2 toxin molecules are simulated using a density-functional quantum-mechanical approach. A numerical experiment aimed at comparing the predicted structural, vibrational and electronic properties of the T-2 toxin with analogous characteristics of the structurally similar 3-deacetylcalonectrin is performed, and the characteristic spectral features that can be used as fingerprints of the T-2 toxin are determined. It is shown that theoretical studies of the structure and spectroscopic features of trichothecene molecules facilitate the development of methods for the detection and characterization of the metabolites.


Introduction
The T-2 toxin and 3-deacetylcalonectrin are A-type trichothecene micotoxines [1][2][3]. These toxins are produced by various Fusarium fungi [1,2]. In particular, 3-deacetylcalonectrin has been extracted from F. graminearum and F. culmorum [4], while the T-2 toxin has been excreted by F. sporotrichioides, F. tricinctum and F. poae [5]. The T-2 is the most toxic of the trichothecenes and poses a serious danger to humans and animals by ingestion of contaminated grains and their processed products [5]. Unlike most biological toxins, T-2 mycotoxin can be absorbed through intact skin. The toxin-releasing moulds grow on a large variety of cereals. Certain toxins produced by the genus Fusarium fungi are characteristic of certain types of cereals, for example, the T-2 toxin often appears in oats [6]. Trichothecenes can cause severe skin disease or intestinal mucosa irritation and diarrhoea. The suppression of animal immune systems has also been noted as a chronic effect of T-2 toxin contamination.
The production of the T-2 mycotoxin occurs during plant infection by the fungus genus Fusarium, followed by deep toxin penetration into the grain [7]. The particular importance in T-2 mycotoxin identification is associated with extensive food production from grains that are potentially susceptible to Fusarium spp. infection, such as barley, wheat and oat [2].
Theoretical studies of the spectroscopic features of the trichothecene structure could facilitate characterization and methods of detection development. To date, the identification vibrational frequencies were multiplied by a scaling factor of 0.983 for both T-2 and 3deacetylcalonectrin. The IR peak half-width at half-height for the demonstrated spectra was 4 cm −1 . For convenient comparison with the experiment in this paper, the scaled frequencies are discussed.
The excitation energies and oscillator strength for 60 singlet excited states were simulated with the use of the TD-DFT approach [26]. The UV-Visible absorption spectra were simulated by assuming a Gaussian band shape with a peak width of 0.3 eV. The Polarizable Continuum Model (PCM) [27] was used for simulating the solvent effect in the methanol media.
All calculations were carried out for isolated molecules, whereas the computational results were compared with experimental data on the structure and IR absorption of crystalline materials. Such inconsistency is due to two circumstances. First, both of the studied compounds are rather complicated. There are 536 atoms in the unit cell of the T-2 toxin crystal. Hence, rigorous periodic DFT modelling is hardly realistic. Second, the intermolecular binding in the T-2 toxin crystal is mainly governed by weak van-der Waals forces [11]. Hence, periodic DFT modelling is unlikely to provide any fundamentally new information on the spectra of the internal vibrational modes. Our main goal is to determine characteristic lines in the IR absorption spectra. In relation to this, internal vibrational modes are mostly related to the characteristic lines in the middle-frequency spectral region. Previous molecular crystal studies have evidenced that the frequencies of the intramolecular vibrational modes weakly depend on the crystalline environment in the absence of strong inter-molecular interaction. This assumption is widely used in discussing the vibrational states of the molecular crystals [28][29][30][31].

Structural Peculiarities of Fused Rings in Free Molecules
The molecular structures of T-2 and 3-deacetylcalonectrin are schematically represented in Figure 1. One can see that the core structures of both molecules consist of three fused rings: cyclohexene (A-ring), tetrahydropyran (B-ring) and cyclopentane (C-ring). The biological activity of the toxins is related to characteristic structural features, which are the epoxy ring, the hydroxyl groups and the acetate group (R2) [24]. The T-2 toxin molecule is bigger-it contains two additional functional groups, R1 and R3. This leads to additional conformational diversity as compared to 3-deacetylcalonectrin [7]. The crystal structure of the T-2 toxin was investigated by the single crystal XRD method [11]. It was found that the lattice is orthorhombic and corresponds to the P2 1 2 1 2 1 space group, with eight molecules in the unit cell. The asymmetric unit contains two molecules (assigned as molecule I and molecule II). They differ in the dihedral angles of the 3-methylbutanoate group (R1). At the same time, it was noted that the molecular structure strongly depends on the environment and on the state of matter [7].
In view of the conformational diversity, we have decided to look for the lowest energy configuration for both molecules. Geometry optimization was carried out for isolated molecules. This corresponds to non-interacting molecule approximation. Such an approach allows us to compare the structural characteristics of the T-2 toxin and the closely related 3-deacetylcalonectrin molecules. Moreover, it is possible to discuss differences in the structure and vibrational spectra of these molecules.
The XRD data for the T-2 crystals indicate the existence of two different molecular configurations. We performed geometry optimization starting from both of them and obtained two stable configurations. One (conformer I) corresponds to the lowest energy structure and is shown in Figure 2b. The other (conformer II) differs in the orientation of the functional groups and is shown in Figure 2c.
green and blue contours highlight the R1, R2 and R3 groups.
The crystal structure of the T-2 toxin was investigated by the single crystal XRD method [11]. It was found that the lattice is orthorhombic and corresponds to the P212121 space group, with eight molecules in the unit cell. The asymmetric unit contains two molecules (assigned as molecule I and molecule II). They differ in the dihedral angles of the 3-methylbutanoate group (R1). At the same time, it was noted that the molecular structure strongly depends on the environment and on the state of matter [7].
In view of the conformational diversity, we have decided to look for the lowest energy configuration for both molecules. Geometry optimization was carried out for isolated molecules. This corresponds to non-interacting molecule approximation. Such an approach allows us to compare the structural characteristics of the T-2 toxin and the closely related 3-deacetylcalonectrin molecules. Moreover, it is possible to discuss differences in the structure and vibrational spectra of these molecules.
The XRD data for the T-2 crystals indicate the existence of two different molecular configurations. We performed geometry optimization starting from both of them and obtained two stable configurations. One (conformer I) corresponds to the lowest energy structure and is shown in Figure 2b. The other (conformer II) differs in the orientation of the functional groups and is shown in Figure 2c. One can see that the core structures (i.e., the structure of the fused rings) of all of the molecules are similar. The T-2 conformers are similar in the following: if a plane is conventionally drawn along the A-ring, the dipole moment of the C=O bond in 3methylbutanoate is directed downward, while in both acetate groups this dipole moment is directed upward in both optimized geometries. The most prominent difference between two conformers of T2 concerns the orientation of the 3-methylbutanoate group (R1). This is in line with experimental data [11] for two molecules in the asymmetric unit. The result allows us to suggest that the two conformers found in our calculations correspond to two nonequivalent molecules in the crystal structure of T-2 . The dihedral angles C23-C31-O46-C47, C22-C36-O39-C40 and C8-C10-O18-C19 for conformer I are equal to  The above-mentioned angles for conformer I are close to those experimentally observed in the XRD study [11]; for molecule I, the angles are equal to 83.4 • , 169.9 • and −138.8 • respectively. However, the experimental values for the second molecule are equal to 118.0 • , 169.0 • and −116.7 • . The divergence in the C8-C10-O18-C19 angle for conformer II is clearly manifested (see Table S1). The observed divergence can be explained by intermolecular interaction, which is not taken into account when performing calculations in the gas phase. Based on NMR studies, conformational flexibility in the 3-methylbutanoate and acetate was noted [7]. This structural dissimilarity moderately affects the dipole moment. The total dipole moment of conformer I is equal to 2.257 Debye, and it equals 2.3315 Debye for conformer II, which is approximately 3.3% higher. The difference in total energy is approximately 0.146 kcal/mol. This is relatively small comparing to the thermal energy at 300 K, which is equal to approximately 0.6 kcal/mol. Such a low energy difference makes the thermodynamic probabilities for the two conformers almost equal.
The core structures (i.e., the three fused rings) in the 3-deacetylcalonectrin and T-2 toxin molecules have several common structural features.
According to Table S1, the shortest bond in the cyclohexene ring (A-ring), corresponding as expected to the C1=C2 double bond, is approximately equal to 1.330 molecules are similar. The T-2 conformers are similar in the following: if a plane is conventionally drawn along the A-ring, the dipole moment of the C=O bond in 3-methylbutanoate is directed downward, while in both acetate groups this dipole moment is directed upward in both optimized geometries. The most prominent difference between two conformers of T2 concerns the orientation of the 3-methylbutanoate group (R1). This is in line with experimental data [11] for two molecules in the asymmetric unit. The result allows us to suggest that the two conformers found in our calculations correspond to two nonequivalent molecules in the crystal structure of T-2 . The dihedral angles C23-C31-O46-C47,  C22-C36-O39-C40 and C8-C10-O18-C19 for conformer I are equal to 87.3°, 172.7° and  −143.5°, and for conformer II are equal to 119.1°, 161.8° and −142.4°, respectively. The above-mentioned angles for conformer I are close to those experimentally observed in the XRD study [11]; for molecule I, the angles are equal to 83.4°, 169.9° and −138.8° respectively. However, the experimental values for the second molecule are equal to 118.0°, 169.0° and −116.7°. The divergence in the C8-C10-O18-C19 angle for conformer II is clearly manifested (see Table S1). The observed divergence can be explained by intermolecular interaction, which is not taken into account when performing calculations in the gas phase. Based on NMR studies, conformational flexibility in the 3-methylbutanoate and acetate was noted [7]. This structural dissimilarity moderately affects the dipole moment. The total dipole moment of conformer I is equal to 2.257 Debye, and it equals 2.3315 Debye for conformer II, which is approximately 3.3% higher. The difference in total energy is approximately 0.146 kcal/mol. This is relatively small comparing to the thermal energy at 300 K, which is equal to approximately 0.6 kcal/mol. Such a low energy difference makes the thermodynamic probabilities for the two conformers almost equal. The core structures (i.e., the three fused rings) in the 3-deacetylcalonectrin and T-2 toxin molecules have several common structural features.
According to Table S1, the shortest bond in the cyclohexene ring (A-ring), corresponding as expected to the C1=C2 double bond, is approximately equal to 1.330 Ǻ for both T-2 toxin and 3-deacetylcalonectrin. The longest bond among the single carbon bonds is predicted at the border of the tetrahydropyran type ring (B-ring); it equals 1.555 Ǻ for the T-2 toxin (C22-C26) and 3-deacetylcalonectrin (C23-C19). According to XRD, the longest C-C bonds are C22-C26 and C22-C23 for the T-2 toxin and C19-C20 and C19-C23 for 3-deacetylcalonectrin. The lengths of the remaining C-C bonds take intermediate values between 1.501 and 1.555 Ǻ. The corresponding planar angles in both molecules have similar values. The large values of the plane angles refer to the C-C=C and C=C-C parts of the A-ring, which is explained by the relatively higher electron density localized in the double carbon bond compared to the single bond. The smallest flat angle is approximately 108° and is formed by two carbon atoms adjoining the B-ring and carbon in the CH2 group (C23-C22-C26).
According to theory, the greatest relative differences between the T-2 toxin molecule and the 3-deacetylcalonectrin occur in the dihedral angles with the carbon atom to which the 3-methylbutanoate group is attached. In the T-2 molecule the electrostatic repulsion between the R1 and R2 groups leads to a shift of the C31 atom. This changes the dihedral angle in the A-ring from almost zero (С23-С1=С2-С28 angle in 3-deacetylcalonectrin) to approximately 1.6-1.7° (С26-С1-С2-С31 angle in T-2 toxin conformers). According to calculations, the largest difference in dihedral angles takes place for the C2-C31-C23-C22 angle in both conformers of the T-2 toxin (about −38.7°) and the corresponding C2-C28-C20-C19 angle in 3-deacetylcalonectrin (approximately −46.1°); see Table S1.
In the XRD experiment, the values of the C2-C31-C23-C22 angle were −41.2° and −36.7° for the two molecules in the asymmetric unit [11]. The calculated values of the angle for both T-2 toxin and 3-deacetylcalonectrin. The longest bond among the single carbon bonds is predicted at the border of the tetrahydropyran type ring (B-ring); it equals 1.555 conformer II (c).
One can see that the core structures (i.e., the structure of the fused rings) of all of the molecules are similar. The T-2 conformers are similar in the following: if a plane is conventionally drawn along the A-ring, the dipole moment of the C=O bond in 3-methylbutanoate is directed downward, while in both acetate groups this dipole moment is directed upward in both optimized geometries. The most prominent difference between two conformers of T2 concerns the orientation of the 3-methylbutanoate group (R1). This is in line with experimental data [11] for two molecules in the asymmetric unit. The result allows us to suggest that the two conformers found in our calculations correspond to two nonequivalent molecules in the crystal structure of T-2. The dihedral angles C23-C31-O46-C47, C22-C36-O39-C40 and C8-C10-O18-C19 for conformer I are equal to 87.3°, 172.7° and −143.5°, and for conformer II are equal to 119.1°, 161.8° and −142.4°, respectively. The above-mentioned angles for conformer I are close to those experimentally observed in the XRD study [11]; for molecule I, the angles are equal to 83.4°, 169.9° and −138.8° respectively. However, the experimental values for the second molecule are equal to 118.0°, 169.0° and −116.7°. The divergence in the C8-C10-O18-C19 angle for conformer II is clearly manifested (see Table S1). The observed divergence can be explained by intermolecular interaction, which is not taken into account when performing calculations in the gas phase. Based on NMR studies, conformational flexibility in the 3-methylbutanoate and acetate was noted [7]. This structural dissimilarity moderately affects the dipole moment. The total dipole moment of conformer I is equal to 2.257 Debye, and it equals 2.3315 Debye for conformer II, which is approximately 3.3% higher. The difference in total energy is approximately 0.146 kcal/mol. This is relatively small comparing to the thermal energy at 300 K, which is equal to approximately 0.6 kcal/mol. Such a low energy difference makes the thermodynamic probabilities for the two conformers almost equal.
The core structures (i.e., the three fused rings) in the 3-deacetylcalonectrin and T-2 toxin molecules have several common structural features.
According to Table S1, the shortest bond in the cyclohexene ring (A-ring), corresponding as expected to the C1=C2 double bond, is approximately equal to 1.330 Ǻ for both T-2 toxin and 3-deacetylcalonectrin. The longest bond among the single carbon bonds is predicted at the border of the tetrahydropyran type ring (B-ring); it equals 1.555 Ǻ for the T-2 toxin (C22-C26) and 3-deacetylcalonectrin (C23-C19). According to XRD, the longest C-C bonds are C22-C26 and C22-C23 for the T-2 toxin and C19-C20 and C19-C23 for 3-deacetylcalonectrin. The lengths of the remaining C-C bonds take intermediate values between 1.501 and 1.555 Ǻ. The corresponding planar angles in both molecules have similar values. The large values of the plane angles refer to the C-C=C and C=C-C parts of the A-ring, which is explained by the relatively higher electron density localized in the double carbon bond compared to the single bond. The smallest flat angle is approximately 108° and is formed by two carbon atoms adjoining the B-ring and carbon in the CH2 group (C23-C22-C26).
According to theory, the greatest relative differences between the T-2 toxin molecule and the 3-deacetylcalonectrin occur in the dihedral angles with the carbon atom to which the 3-methylbutanoate group is attached. In the T-2 molecule the electrostatic repulsion between the R1 and R2 groups leads to a shift of the C31 atom. This changes the dihedral angle in the A-ring from almost zero (С23-С1=С2-С28 angle in 3-deacetylcalonectrin) to approximately 1.6-1.7° (С26-С1-С2-С31 angle in T-2 toxin conformers). According to calculations, the largest difference in dihedral angles takes place for the C2-C31-C23-C22 angle in both conformers of the T-2 toxin (about −38.7°) and the corresponding C2-C28-C20-C19 angle in 3-deacetylcalonectrin (approximately −46.1°); see Table S1.
In the XRD experiment, the values of the C2-C31-C23-C22 angle were −41.2° and −36.7° for the two molecules in the asymmetric unit [11]. The calculated values of the angle for the T-2 toxin (C22-C26) and 3-deacetylcalonectrin (C23-C19). According to XRD, the longest C-C bonds are C22-C26 and C22-C23 for the T-2 toxin and C19-C20 and C19-C23 for 3-deacetylcalonectrin. The lengths of the remaining C-C bonds take intermediate values between 1.501 and 1.555 for T-2 molecule (blue, red and black numbers are representation of carbon, atoms, respectively) (a) and stereochemical pictures for conformer I (b) and t the core structures (i.e., the structure of the fused rings) of all of the r. The T-2 conformers are similar in the following: if a plane is conlong the A-ring, the dipole moment of the C=O bond in 3-methylbutanward, while in both acetate groups this dipole moment is directed mized geometries. The most prominent difference between two conns the orientation of the 3-methylbutanoate group (R1). This is in line ata [11] for two molecules in the asymmetric unit. The result allows e two conformers found in our calculations correspond to two nonein the crystal structure of T-2. The dihedral angles C23-C31-O46-C47, nd C8-C10-O18-C19 for conformer I are equal to 87.3°, 172.7° and former II are equal to 119.1°, 161.8° and −142.4°, respectively. The gles for conformer I are close to those experimentally observed in the molecule I, the angles are equal to 83.4°, 169.9° and −138.8° respece experimental values for the second molecule are equal to 118.0°, he divergence in the C8-C10-O18-C19 angle for conformer II is clearly le S1). The observed divergence can be explained by intermolecular s not taken into account when performing calculations in the gas R studies, conformational flexibility in the 3-methylbutanoate and ]. This structural dissimilarity moderately affects the dipole moment. ent of conformer I is equal to 2.257 Debye, and it equals 2.3315 Debye ich is approximately 3.3% higher. The difference in total energy is kcal/mol. This is relatively small comparing to the thermal energy at l to approximately 0.6 kcal/mol. Such a low energy difference makes probabilities for the two conformers almost equal. ures (i.e., the three fused rings) in the 3-deacetylcalonectrin and T-2 e several common structural features. able S1, the shortest bond in the cyclohexene ring (A-ring), corred to the C1=C2 double bond, is approximately equal to 1.330 Ǻ for deacetylcalonectrin. The longest bond among the single carbon bonds order of the tetrahydropyran type ring (B-ring); it equals 1.555 Ǻ for 26) and 3-deacetylcalonectrin (C23-C19). According to XRD, the long-22-C26 and C22-C23 for the T-2 toxin and C19-C20 and C19-C23 for n. The lengths of the remaining C-C bonds take intermediate values .555 Ǻ. The corresponding planar angles in both molecules have sime values of the plane angles refer to the C-C=C and C=C-C parts of explained by the relatively higher electron density localized in the compared to the single bond. The smallest flat angle is approximately y two carbon atoms adjoining the B-ring and carbon in the CH2 group eory, the greatest relative differences between the T-2 toxin molecule lonectrin occur in the dihedral angles with the carbon atom to which te group is attached. In the T-2 molecule the electrostatic repulsion R2 groups leads to a shift of the C31 atom. This changes the dihedral from almost zero (С23-С1=С2-С28 angle in 3-deacetylcalonectrin) to .7° (С26-С1-С2-С31 angle in T-2 toxin conformers). According to caldifference in dihedral angles takes place for the C2-C31-C23-C22 anrs of the T-2 toxin (about −38.7°) and the corresponding C2-C28-C20tylcalonectrin (approximately −46.1°); see Table S1. eriment, the values of the C2-C31-C23-C22 angle were −41.2° and lecules in the asymmetric unit [11]. The calculated values of the angle . The corresponding planar angles in both molecules have similar values. The large values of the plane angles refer to the C-C=C and C=C-C parts of the A-ring, which is explained by the relatively higher electron density localized in the double carbon bond compared to the single bond. The smallest flat angle is approximately 108 • and is formed by two carbon atoms adjoining the B-ring and carbon in the CH 2 group (C23-C22-C26).
According to theory, the greatest relative differences between the T-2 toxin molecule and the 3-deacetylcalonectrin occur in the dihedral angles with the carbon atom to which the 3-methylbutanoate group is attached. In the T-2 molecule the electrostatic repulsion between the R1 and R2 groups leads to a shift of the C31 atom. This changes the dihedral angle in the A-ring from almost zero (C23-C1=C2-C28 angle in 3-deacetylcalonectrin) to approximately 1.6-1.7 • (C26-C1-C2-C31 angle in T-2 toxin conformers). According to calculations, the largest difference in dihedral angles takes place for the C2-C31-C23-C22 angle in both conformers of the T-2 toxin (about −38.7 • ) and the corresponding C2-C28-C20-C19 angle in 3-deacetylcalonectrin (approximately −46.1 • ); see Table S1.
In the XRD experiment, the values of the C2-C31-C23-C22 angle were −41.2 • and −36.7 • for the two molecules in the asymmetric unit [11]. The calculated values of the angle for both T-2 conformers is approximately −38.7 o , which is closer to the averaged value over two experimental molecules.
The difference between the optimized bond lengths in the B-rings of the two toxin divergences is within a few thousandths of an angstrom. It is less than 0.62 • for the planar angles.
Let us consider dihedral angles involving the C19 atom (in the notations of Figure 2a) and the C22 atom (in the notations of Figure 3). These angles have similar values (within 0.3 • divergence) for the two T-2 toxin conformers. At the same time, the divergence increases up to 1-1.5 • when comparing the angles of the T-2 toxin with the corresponding angles of 3-deacetylcalonectrin. divergences is within a few thousandths of an angstrom. It is less than 0.62° for the planar angles.
Let us consider dihedral angles involving the C19 atom (in the notations of Figure  2a) and the C22 atom (in the notations of Figure 3). These angles have similar values (within 0.3° divergence) for the two T-2 toxin conformers. At the same time, the divergence increases up to 1-1.5° when comparing the angles of the T-2 toxin with the corresponding angles of 3-deacetylcalonectrin. Such a small difference could be due to the R2 group attached to these carbon atoms. The electrostatic interaction of this group with the R1 and R3 groups nearby is small. The shortest bonds correspond to C-O in the B-ring. The C-C bonds in both types of toxins lie in the range of 1.51 to 1.59 Ǻ.
A slight overestimation in the calculated bond lengths in comparison with the X-ray structural data can be noted for the B-ring of the T-2 toxin. The largest error is approximately 1.7% for the C21-C6 bond. The highest divergences (up to 2°) were noted for the planar angles. The highest divergence for dihedral angles takes place for the C21-C6-O5-C26 angle.
According to theory, the values of bond lengths and planar and dihedral angles in the cyclopentyl ring (C-ring) are close in both T-2 conformers. Difference within 0.66% were manifested between the C-C bonds of the T-2 conformers and 3-deacetylcalonectrin. The largest value takes place for the C3-C10 bonds. The smallest difference between dihedral angles in the C-ring of the T-2 conformers and 3-deacetylcalonectrin takes place for the angle with the C10 atom (approximately 0.05°). For the remaining dihedral angles, the difference varies from 2.2° to 4.2°; see Table S1. The divergence between the theoretical and experimental values of the dihedral angles in the C-ring lies within 1-1.5° for the T-2 toxin molecule. The largest variance (approximately 2%) in bond lengths is noted for the C8-C10 bond.
The calculated C-O bond in the epoxy ring (E-ring) is shorter than in the experiment. We suppose that this is due to the absence of the hydrogen bond through the O63 atom. Such a small difference could be due to the R2 group attached to these carbon atoms. The electrostatic interaction of this group with the R1 and R3 groups nearby is small. The shortest bonds correspond to C-O in the B-ring. The C-C bonds in both types of toxins lie in the range of 1.51 to 1.59 atoms, respectively) (a) and stereochemical pictures for conformer I (b) and t the core structures (i.e., the structure of the fused rings) of all of the ar. The T-2 conformers are similar in the following: if a plane is conlong the A-ring, the dipole moment of the C=O bond in 3-methylbutawnward, while in both acetate groups this dipole moment is directed imized geometries. The most prominent difference between two conrns the orientation of the 3-methylbutanoate group (R1). This is in line ata [11] for two molecules in the asymmetric unit. The result allows e two conformers found in our calculations correspond to two nonein the crystal structure of T-2. The dihedral angles C23-C31-O46-C47, nd C8-C10-O18-C19 for conformer I are equal to 87.3°, 172.7° and former II are equal to 119.1°, 161.8° and −142.4°, respectively. The gles for conformer I are close to those experimentally observed in the molecule I, the angles are equal to 83.4°, 169.9° and −138.8° respece experimental values for the second molecule are equal to 118.0°, he divergence in the C8-C10-O18-C19 angle for conformer II is clearly le S1). The observed divergence can be explained by intermolecular is not taken into account when performing calculations in the gas R studies, conformational flexibility in the 3-methylbutanoate and ]. This structural dissimilarity moderately affects the dipole moment. ent of conformer I is equal to 2.257 Debye, and it equals 2.3315 Debye hich is approximately 3.3% higher. The difference in total energy is kcal/mol. This is relatively small comparing to the thermal energy at al to approximately 0.6 kcal/mol. Such a low energy difference makes probabilities for the two conformers almost equal. ures (i.e., the three fused rings) in the 3-deacetylcalonectrin and T-2 e several common structural features. able S1, the shortest bond in the cyclohexene ring (A-ring), correed to the C1=C2 double bond, is approximately equal to 1.330 Ǻ for -deacetylcalonectrin. The longest bond among the single carbon bonds order of the tetrahydropyran type ring (B-ring); it equals 1.555 Ǻ for 26) and 3-deacetylcalonectrin (C23-C19). According to XRD, the long-22-C26 and C22-C23 for the T-2 toxin and C19-C20 and C19-C23 for n. The lengths of the remaining C-C bonds take intermediate values .555 Ǻ. The corresponding planar angles in both molecules have sime values of the plane angles refer to the C-C=C and C=C-C parts of s explained by the relatively higher electron density localized in the compared to the single bond. The smallest flat angle is approximately y two carbon atoms adjoining the B-ring and carbon in the CH2 group eory, the greatest relative differences between the T-2 toxin molecule lonectrin occur in the dihedral angles with the carbon atom to which ate group is attached. In the T-2 molecule the electrostatic repulsion R2 groups leads to a shift of the C31 atom. This changes the dihedral from almost zero (С23-С1=С2-С28 angle in 3-deacetylcalonectrin) to .7° (С26-С1-С2-С31 angle in T-2 toxin conformers). According to calt difference in dihedral angles takes place for the C2-C31-C23-C22 aners of the T-2 toxin (about −38.7°) and the corresponding C2-C28-C20etylcalonectrin (approximately −46.1°); see Table S1. periment, the values of the C2-C31-C23-C22 angle were −41.2° and olecules in the asymmetric unit [11]. The calculated values of the angle . A slight overestimation in the calculated bond lengths in comparison with the X-ray structural data can be noted for the B-ring of the T-2 toxin. The largest error is approximately 1.7% for the C21-C6 bond. The highest divergences (up to 2 • ) were noted for the planar angles. The highest divergence for dihedral angles takes place for the C21-C6-O5-C26 angle.
According to theory, the values of bond lengths and planar and dihedral angles in the cyclopentyl ring (C-ring) are close in both T-2 conformers. Difference within 0.66% were manifested between the C-C bonds of the T-2 conformers and 3-deacetylcalonectrin. The largest value takes place for the C3-C10 bonds. The smallest difference between dihedral angles in the C-ring of the T-2 conformers and 3-deacetylcalonectrin takes place for the angle with the C10 atom (approximately 0.05 • ). For the remaining dihedral angles, the difference varies from 2.2 • to 4.2 • ; see Table S1. The divergence between the theoretical and experimental values of the dihedral angles in the C-ring lies within 1-1.5 • for the T-2 toxin molecule. The largest variance (approximately 2%) in bond lengths is noted for the C8-C10 bond.
The calculated C-O bond in the epoxy ring (E-ring) is shorter than in the experiment. We suppose that this is due to the absence of the hydrogen bond through the O63 atom. A significantly short contact between the OH group of molecule 1 and the oxygen atoms of the epoxy group of molecule 2 was noted in the real crystal structure [7].

Arrangement of 3-Methylbutanoate and Acetate Groups
As stated above, the 3-methylbutanoate group (R1 group) has conformational flexibility. The calculated values of the dihedral angles near the A-ring are in good agreement with the experimental data. Additionally, the value of the C23-C31-C46-C47 angle fits the experiment for two conformers relatively well (the divergence lies within 4 • ). At the same time, the hydrogen bonds with the terminal CH 3 group, which affects the dihedral angles C46-C47-C49-C52, were not taken into account in the gas phase calculations. In particular, these dihedral angles are equal to approximately −154 • for both conformers, while the experimental data provide value of −77.7 • for molecule 1 and −87.2 • for molecule 2. At the same time, the difference between the experiment and theory for the dihedral angle C47-C49-C52-C53 lies within 8-10 • . Presumably, the effect of this divergence involves the frequencies of stretching and bending vibrations in CH 3 , slightly increasing them due to the absence of hydrogen interaction. At the same time, the characteristic vibrational frequencies in the torsion angles with carbon atoms have relatively low values (usually less than 300 cm −1 ) and low polarity, which is expected to have a low influence on IR active modes.
For the R2 functional group, the highest divergence between theory and experiment concerns the value of the C22-C36-O39-C40 angle in conformers II of the T-2 toxin. The calculated value of this angle is equal to 172.7 • for conformer I, whereas the experiment gives 169.9 • for molecule 1. In turn, the calculated value of this angle is equal to 161.8 • for conformer II, whereas the experiment gives 169.
In the case of 3-deacetylcalonectrin, a similar functional group has flatter dihedral angles C22-C36-O39-C40 and C36-O39-C40-C42, which are close to 180 o . The C3-C22-C36-O39 angle is less by 7-8 • , which means a reversal of this radical towards the A-ring rather than in the direction of the methyl group attached to the B-ring.
The calculated parameters of the acetate group attached to the C-ring are very similar for both T-2 toxin conformers. The main differences concern the CO bond length, linking the acetate group and the C-ring, as well as the dihedral angle C8-C10-O18-C19, which determines the angle of rotation relative to the ring plane. For the T-2 toxin conformer I, the dihedral angle C8-C10-O18-C19 is equal to −143.5 • , and the length of the C10-C18 bond is equal to 1.442 • . This is close to experimental XRD data from which the corresponding dihedral angle is −138.8 • and the bond length is 1.439 One can see that the core structures (i.e., the structure of the fused rings) of all of the molecules are similar. The T-2 conformers are similar in the following: if a plane is conventionally drawn along the A-ring, the dipole moment of the C=O bond in 3-methylbutanoate is directed downward, while in both acetate groups this dipole moment is directed upward in both optimized geometries. The most prominent difference between two conformers of T2 concerns the orientation of the 3-methylbutanoate group (R1). This is in line with experimental data [11] for two molecules in the asymmetric unit. The result allows us to suggest that the two conformers found in our calculations correspond to two nonequivalent molecules in the crystal structure of T-2. The dihedral angles C23-C31-O46-C47, C22-C36-O39-C40 and C8-C10-O18-C19 for conformer I are equal to 87.3°, 172.7° and −143.5°, and for conformer II are equal to 119.1°, 161.8° and −142.4°, respectively. The above-mentioned angles for conformer I are close to those experimentally observed in the XRD study [11]; for molecule I, the angles are equal to 83.4°, 169.9° and −138.8° respectively. However, the experimental values for the second molecule are equal to 118.0°, 169.0° and −116.7°. The divergence in the C8-C10-O18-C19 angle for conformer II is clearly manifested (see Table S1). The observed divergence can be explained by intermolecular interaction, which is not taken into account when performing calculations in the gas phase. Based on NMR studies, conformational flexibility in the 3-methylbutanoate and acetate was noted [7]. This structural dissimilarity moderately affects the dipole moment. The total dipole moment of conformer I is equal to 2.257 Debye, and it equals 2.3315 Debye for conformer II, which is approximately 3.3% higher. The difference in total energy is approximately 0.146 kcal/mol. This is relatively small comparing to the thermal energy at 300 K, which is equal to approximately 0.6 kcal/mol. Such a low energy difference makes the thermodynamic probabilities for the two conformers almost equal.
The core structures (i.e., the three fused rings) in the 3-deacetylcalonectrin and T-2 toxin molecules have several common structural features.
According to Table S1, the shortest bond in the cyclohexene ring (A-ring), corresponding as expected to the C1=C2 double bond, is approximately equal to 1.330 Ǻ for both T-2 toxin and 3-deacetylcalonectrin. The longest bond among the single carbon bonds is predicted at the border of the tetrahydropyran type ring (B-ring); it equals 1.555 Ǻ for the T-2 toxin (C22-C26) and 3-deacetylcalonectrin (C23-C19). According to XRD, the longest C-C bonds are C22-C26 and C22-C23 for the T-2 toxin and C19-C20 and C19-C23 for 3-deacetylcalonectrin. The lengths of the remaining C-C bonds take intermediate values between 1.501 and 1.555 Ǻ. The corresponding planar angles in both molecules have similar values. The large values of the plane angles refer to the C-C=C and C=C-C parts of the A-ring, which is explained by the relatively higher electron density localized in the double carbon bond compared to the single bond. The smallest flat angle is approximately 108° and is formed by two carbon atoms adjoining the B-ring and carbon in the CH2 group (C23-C22-C26).
According to theory, the greatest relative differences between the T-2 toxin molecule and the 3-deacetylcalonectrin occur in the dihedral angles with the carbon atom to which the 3-methylbutanoate group is attached. In the T-2 molecule the electrostatic repulsion between the R1 and R2 groups leads to a shift of the C31 atom. This changes the dihedral angle in the A-ring from almost zero (С23-С1=С2-С28 angle in 3-deacetylcalonectrin) to approximately 1.6-1.7° (С26-С1-С2-С31 angle in T-2 toxin conformers). According to calculations, the largest difference in dihedral angles takes place for the C2-C31-C23-C22 angle in both conformers of the T-2 toxin (about −38.7°) and the corresponding C2-C28-C20-C19 angle in 3-deacetylcalonectrin (approximately −46.1°); see Table S1.
In the XRD experiment, the values of the C2-C31-C23-C22 angle were −41.2° and −36.7° for the two molecules in the asymmetric unit [11]. The calculated values of the angle  Atomic labels for T-2 molecule (blue, red and black numbers are representation of carbon, oxygen and hydrogen atoms, respectively) (a) and stereochemical pictures for conformer I (b) and conformer II (c).
One can see that the core structures (i.e., the structure of the fused rings) of all of the molecules are similar. The T-2 conformers are similar in the following: if a plane is conventionally drawn along the A-ring, the dipole moment of the C=O bond in 3-methylbutanoate is directed downward, while in both acetate groups this dipole moment is directed upward in both optimized geometries. The most prominent difference between two conformers of T2 concerns the orientation of the 3-methylbutanoate group (R1). This is in line with experimental data [11] for two molecules in the asymmetric unit. The result allows us to suggest that the two conformers found in our calculations correspond to two nonequivalent molecules in the crystal structure of T-2. The dihedral angles C23-C31-O46-C47, C22-C36-O39-C40 and C8-C10-O18-C19 for conformer I are equal to 87.3°, 172.7° and −143.5°, and for conformer II are equal to 119.1°, 161.8° and −142.4°, respectively. The above-mentioned angles for conformer I are close to those experimentally observed in the XRD study [11]; for molecule I, the angles are equal to 83.4°, 169.9° and −138.8° respectively. However, the experimental values for the second molecule are equal to 118.0°, 169.0° and −116.7°. The divergence in the C8-C10-O18-C19 angle for conformer II is clearly manifested (see Table S1). The observed divergence can be explained by intermolecular interaction, which is not taken into account when performing calculations in the gas phase. Based on NMR studies, conformational flexibility in the 3-methylbutanoate and acetate was noted [7]. This structural dissimilarity moderately affects the dipole moment. The total dipole moment of conformer I is equal to 2.257 Debye, and it equals 2.3315 Debye for conformer II, which is approximately 3.3% higher. The difference in total energy is approximately 0.146 kcal/mol. This is relatively small comparing to the thermal energy at 300 K, which is equal to approximately 0.6 kcal/mol. Such a low energy difference makes the thermodynamic probabilities for the two conformers almost equal.
The core structures (i.e., the three fused rings) in the 3-deacetylcalonectrin and T-2 toxin molecules have several common structural features.
According to Table S1, the shortest bond in the cyclohexene ring (A-ring), corresponding as expected to the C1=C2 double bond, is approximately equal to 1.330 Ǻ for both T-2 toxin and 3-deacetylcalonectrin. The longest bond among the single carbon bonds is predicted at the border of the tetrahydropyran type ring (B-ring); it equals 1.555 Ǻ for the T-2 toxin (C22-C26) and 3-deacetylcalonectrin (C23-C19). According to XRD, the longest C-C bonds are C22-C26 and C22-C23 for the T-2 toxin and C19-C20 and C19-C23 for 3-deacetylcalonectrin. The lengths of the remaining C-C bonds take intermediate values between 1.501 and 1.555 Ǻ. The corresponding planar angles in both molecules have similar values. The large values of the plane angles refer to the C-C=C and C=C-C parts of the A-ring, which is explained by the relatively higher electron density localized in the double carbon bond compared to the single bond. The smallest flat angle is approximately 108° and is formed by two carbon atoms adjoining the B-ring and carbon in the CH2 group (C23-C22-C26).
According to theory, the greatest relative differences between the T-2 toxin molecule and the 3-deacetylcalonectrin occur in the dihedral angles with the carbon atom to which the 3-methylbutanoate group is attached. In the T-2 molecule the electrostatic repulsion between the R1 and R2 groups leads to a shift of the C31 atom. This changes the dihedral angle in the A-ring from almost zero (С23-С1=С2-С28 angle in 3-deacetylcalonectrin) to approximately 1.6-1.7° (С26-С1-С2-С31 angle in T-2 toxin conformers). According to calculations, the largest difference in dihedral angles takes place for the C2-C31-C23-C22 angle in both conformers of the T-2 toxin (about −38.7°) and the corresponding C2-C28-C20-C19 angle in 3-deacetylcalonectrin (approximately −46.1°); see Table S1.
In the XRD experiment, the values of the C2-C31-C23-C22 angle were −41.2° and −36.7° for the two molecules in the asymmetric unit [11]. The calculated values of the angle , which are close to conformer I. At the same time, the XRD data for molecule 2 provide a bond length of 1.46 Figure 2. Atomic labels for T-2 molecule (blue, red and black numbers are representation of carbon, oxygen and hydrogen atoms, respectively) (a) and stereochemical pictures for conformer I (b) and conformer II (c).
One can see that the core structures (i.e., the structure of the fused rings) of all of the molecules are similar. The T-2 conformers are similar in the following: if a plane is conventionally drawn along the A-ring, the dipole moment of the C=O bond in 3-methylbutanoate is directed downward, while in both acetate groups this dipole moment is directed upward in both optimized geometries. The most prominent difference between two conformers of T2 concerns the orientation of the 3-methylbutanoate group (R1). This is in line with experimental data [11] for two molecules in the asymmetric unit. The result allows us to suggest that the two conformers found in our calculations correspond to two nonequivalent molecules in the crystal structure of T-2. The dihedral angles C23-C31-O46-C47, C22-C36-O39-C40 and C8-C10-O18-C19 for conformer I are equal to 87.3°, 172.7° and −143.5°, and for conformer II are equal to 119.1°, 161.8° and −142.4°, respectively. The above-mentioned angles for conformer I are close to those experimentally observed in the XRD study [11]; for molecule I, the angles are equal to 83.4°, 169.9° and −138.8° respectively. However, the experimental values for the second molecule are equal to 118.0°, 169.0° and −116.7°. The divergence in the C8-C10-O18-C19 angle for conformer II is clearly manifested (see Table S1). The observed divergence can be explained by intermolecular interaction, which is not taken into account when performing calculations in the gas phase. Based on NMR studies, conformational flexibility in the 3-methylbutanoate and acetate was noted [7]. This structural dissimilarity moderately affects the dipole moment. The total dipole moment of conformer I is equal to 2.257 Debye, and it equals 2.3315 Debye for conformer II, which is approximately 3.3% higher. The difference in total energy is approximately 0.146 kcal/mol. This is relatively small comparing to the thermal energy at 300 K, which is equal to approximately 0.6 kcal/mol. Such a low energy difference makes the thermodynamic probabilities for the two conformers almost equal.
The core structures (i.e., the three fused rings) in the 3-deacetylcalonectrin and T-2 toxin molecules have several common structural features.
According to Table S1, the shortest bond in the cyclohexene ring (A-ring), corresponding as expected to the C1=C2 double bond, is approximately equal to 1.330 Ǻ for both T-2 toxin and 3-deacetylcalonectrin. The longest bond among the single carbon bonds is predicted at the border of the tetrahydropyran type ring (B-ring); it equals 1.555 Ǻ for the T-2 toxin (C22-C26) and 3-deacetylcalonectrin (C23-C19). According to XRD, the longest C-C bonds are C22-C26 and C22-C23 for the T-2 toxin and C19-C20 and C19-C23 for 3-deacetylcalonectrin. The lengths of the remaining C-C bonds take intermediate values between 1.501 and 1.555 Ǻ. The corresponding planar angles in both molecules have similar values. The large values of the plane angles refer to the C-C=C and C=C-C parts of the A-ring, which is explained by the relatively higher electron density localized in the double carbon bond compared to the single bond. The smallest flat angle is approximately 108° and is formed by two carbon atoms adjoining the B-ring and carbon in the CH2 group (C23-C22-C26).
According to theory, the greatest relative differences between the T-2 toxin molecule and the 3-deacetylcalonectrin occur in the dihedral angles with the carbon atom to which the 3-methylbutanoate group is attached. In the T-2 molecule the electrostatic repulsion between the R1 and R2 groups leads to a shift of the C31 atom. This changes the dihedral angle in the A-ring from almost zero (С23-С1=С2-С28 angle in 3-deacetylcalonectrin) to approximately 1.6-1.7° (С26-С1-С2-С31 angle in T-2 toxin conformers). According to calculations, the largest difference in dihedral angles takes place for the C2-C31-C23-C22 angle in both conformers of the T-2 toxin (about −38.7°) and the corresponding C2-C28-C20-C19 angle in 3-deacetylcalonectrin (approximately −46.1°); see Table S1.
In the XRD experiment, the values of the C2-C31-C23-C22 angle were −41.2° and −36.7° for the two molecules in the asymmetric unit [11]. The calculated values of the angle and a dihedral angle of −116.7 • , which differ markedly from the analogue dihedral angle in conformer II.
Based on the aforementioned reasons, conformer I of the T-2 toxin, which is structurally close to one of the molecules obtained in XRD [3] and slightly more stable than conformer II, was chosen to model the IR absorption spectrum.

Vibrational Properties
Currently, experimental information on the vibrational properties of the T-2 toxin is rather scarce. Inversely, the vibrational states of the structural analogue of 3-deacetylcalonectrin have been thoroughly studied [3]. Thus, we have decided to perform a comparative analysis of the simulated vibrational spectra of these compounds, taking into account the experimental data of [3].

3-Deacetylcalonectrin Molecule
The experimental IR absorption spectrum in the range 750-4000 cm −1 was reported in [3]. The spectral curve contains several absorption bands. Their frequency positions are listed in Table 1. Some of these bands are quite wide, which suggest that several vibrational modes contribute to these bands. This circumstance was taken into account when investigating the calculated and observed spectra. The simulated IR absorbance spectrum is shown in Figure 4a. Supposed assignment of the observed IR absorption bands to the calculated normal modes is presented in the Table 1. The highest-frequency peak of 3626 cm −1 in the experimental IR absorption spectra is attributed to the O-H bond stretching vibration. The calculated frequency of this mode is 3707 cm −1 . Such overestimation of the frequency value is typically for these modes due to the high anharmonicity of the X-H bond stretching vibrations.
The experimental IR absorption spectrum exhibits a wide feature with a maximum at 2962 cm −1 . Various C-H stretching vibrations are located in the range of 2800-3100 cm −1 . The 2962 cm −1 absorption maximum corresponds to the antisymmetric vibrations of the methyl groups. The peak has well-pronounced shoulders at 2920 cm −1 and approximately 3030 cm −1 . The former is related to the antisymmetric vibrations in the CH2 group. Interpretation of the latter requires further calculations. The additional shoulder near 2870 cm −1 corresponds to symmetric stretching vibrations in the CH2 and CH3 groups [32].
The corresponding vibrational modes in the calculated spectrum are located at 2949-3114 cm −1 . The highest frequency mode at 3167 cm −1 corresponds to the antisymmetric ν(CH2) mode in the epoxy ring. Its frequency is predicted to be higher than the stretching vibration in ν(=C-H). The shoulder at 3030 cm −1 observed in the IR spectrum from [3] was interpreted as absorption associated with antisymmetric vibration in the epoxy group CH2 and the stretching vibration of hydrogen in the =C-H unit.
According to theory, the frequencies of the symmetric stretching vibrations in the CH2 and CH3 groups, as well as CH stretching vibrations, are 2949-3044 cm −1 . For the vast majority of stretching hydrogen vibrational modes, the following relation (1) is predicted: νas(CH2) epoxy > ν(=C-H) > νas(CH3), νas(CH2) > νs(CH3), νs(CH2), ν(CH) (1)  The 3-deacetylcalonectrin molecule contains 46 atoms, resulting in 132 vibrational modes. In order to gain a better understanding of the spectral region of the most IR active modes, a set of the simple models with characteristic polar groups were tested. In the case of 3-deacetylcalonectrin, these models are ethyl acetate (see Figure 4), tetrahydropyran and cyclopentanol (see Figure S1). Vibrational modes with close frequencies form characteristic spectral features.
The highest-frequency peak of 3626 cm −1 in the experimental IR absorption spectra is attributed to the O-H bond stretching vibration. The calculated frequency of this mode is 3707 cm −1 . Such overestimation of the frequency value is typically for these modes due to the high anharmonicity of the X-H bond stretching vibrations.
The experimental IR absorption spectrum exhibits a wide feature with a maximum at 2962 cm −1 . Various C-H stretching vibrations are located in the range of 2800-3100 cm −1 . The 2962 cm −1 absorption maximum corresponds to the antisymmetric vibrations of the methyl groups. The peak has well-pronounced shoulders at 2920 cm −1 and approximately 3030 cm −1 . The former is related to the antisymmetric vibrations in the CH 2 group. Interpretation of the latter requires further calculations. The additional shoulder near 2870 cm −1 corresponds to symmetric stretching vibrations in the CH 2 and CH 3 groups [32].
The corresponding vibrational modes in the calculated spectrum are located at 2949-3114 cm −1 . The highest frequency mode at 3167 cm −1 corresponds to the antisymmetric ν(CH 2 ) mode in the epoxy ring. Its frequency is predicted to be higher than the stretching vibration in ν(=C-H). The shoulder at 3030 cm −1 observed in the IR spectrum from [3] was interpreted as absorption associated with antisymmetric vibration in the epoxy group CH 2 and the stretching vibration of hydrogen in the =C-H unit.
According to theory, the frequencies of the symmetric stretching vibrations in the CH 2 and CH 3 groups, as well as CH stretching vibrations, are 2949-3044 cm −1 . For the vast majority of stretching hydrogen vibrational modes, the following relation (1) is predicted: The experimental IR spectrum of [3] contains an intense absorption peak at 1767 cm −1 and a weak satellite peak at 1709 cm −1 , which were assigned to the v(C=O) and v(C=C), vibrations, respectively. The calculation predicts similar absorption spectrum with analogous peaks at 1772 cm −1 and 1703 cm −1 . Calculations show that the vibration in the ethyl acetate molecule localized within the C=O bond has close frequency of 1772 cm −1 (see Figure 5).

Materials 2022, 15, x FOR PEER REVIEW
The experimental IR spectrum of [3] contains an intense absorption peak at 1 and a weak satellite peak at 1709 cm −1 , which were assigned to the v(C=O) and vibrations, respectively. The calculation predicts similar absorption spectrum wi gous peaks at 1772 cm −1 and 1703 cm −1 . Calculations show that the vibration in acetate molecule localized within the C=O bond has close frequency of 1772 cm −1 ure 5). The broad band near 1468 cm −1 in the experimental spectrum corresponds well-resolved bands, with maxima at 1495 cm −1 , 1470 cm −1 and 1454 cm −1 in the c spectrum. These bands are related to the δ(CH2) modes and the asymmetric δ(CH The coupling of such vibrational modes is due to the closeness of their frequencie to situations in other organic compounds [32]. The broad band near 1468 cm −1 in the experimental spectrum corresponds to three well-resolved bands, with maxima at 1495 cm −1 , 1470 cm −1 and 1454 cm −1 in the calculated spectrum. These bands are related to the δ(CH 2 ) modes and the asymmetric δ(CH 3 ) modes. The coupling of such vibrational modes is due to the closeness of their frequencies, similar to situations in other organic compounds [32].
The absorption feature at 1387 cm −1 in the experimental spectrum is rather wide (approximately 100 cm −1 ). Ten modes contribute to the IR absorption within this frequency interval in the calculated spectrum. Their cumulative effect gives rise to a spectral band with two maxima located at 1371 cm -1 and 1430 cm -1 .
Several modes contribute to the IR absorption band at 1430 cm −1 . These modes have a mixed character, with the main contributions from δ(HCH) and symmetric δ s (CH 3 ) vibrations. Furthermore, a small addition of δ(CCH), δ(HCO) and δ(COH) was noticed. Notice that a value of 1430 cm −1 lies between the characteristic frequencies of the δ(CH 2 ) (1450 ± 10 cm −1 ) and δ s (CH 3 ) (1375 ± 10 cm −1 ) modes [32]. In the low frequency shoulder of the theoretical spectrum (1320-1357 cm −1 ), a significant contribution arises from the δ(CCH) vibrations. In these modes, the w(HCH) and τ(HCH) vibrations, having lower vibration frequencies, also begin to participate.
The most intense band in the experimental absorbance spectrum is located at 1227 cm −1 . The analogous band in the calculated IR spectrum corresponds to the vibrational mode, with a frequency of 1233 cm −1 . This is the stretching vibration of the single C41-O40 bond in the acetate group. Due to the high polarity of the C-O bond, this mode is very active in IR absorption. More precisely, this is an asymmetric bond stretching mode localized within the C-C-O bridge. The mode also includes some contribution from the δ(H 44 C 43 C 41 ) bending vibration, with the small addition of δ(H 45 C 43 H 46 ) and w(H 30 C 29 H 31 ) vibrations. The most IR-active mode of the ethyl acetate molecule has exactly the same atomic displacements and frequency. The difference is within 0.7%.
Near the intense IR peak, there is a weak band at 1177 cm −1 . It can be assigned to the superposition of three lines around 1160 cm −1 in the calculated IR spectrum. These modes are delocalized over the whole molecule. A characteristic feature of these modes is the contribution of the C-C bond stretching vibrations of the four-fold coordinated C atoms.
Another significant band in the experimental absorption spectrum of 3-deactylcalonectrin is located at 1054 cm −1 . The corresponding vibrational modes in the theoretical spectrum are located at 1013-1070 cm −1 . Analysis of the calculated eigenvectors shows that the participating modes have two characteristic contributions (see Table 1). The first one is the hydrogen atom vibrations of a mostly rocking type, with the smaller addition of twisting and wagging types in the CH 2 groups. Another contribution is due to non-hydrogen vibrations, especially localized in stretching single C-C and C-O bonds. The calculated vector of the dipole derivative is collinear with the C-O bonds.
It is predicted that the highest frequency vibrational mode in the region of 1013-1070 cm −1 is located at 1062 cm −1 . The significant contribution to this mode comes from vibrations of the oxygen atom in the B-ring. The atomic displacements in the COC group are asymmetric and close to those observed for corresponding modes localized in the tetrahydropyran ring (see Figure S2 in Supporting Information).
According to calculations, there are two intense peaks in the IR spectrum of 3deacetylcalonectrin in the region of 1013-1060 cm −1 . Their frequencies are 1040 and 1044 cm −1 . The 1040 cm −1 mode has a contribution of the v(C 29 O 40 ) mode, localized in the acetate group. Similar atomic displacements were observed for a mode located at 1047 cm −1 in the ethyl acetate molecule. The greater difference, as compared with another modes in Table S2, can be explained by the significant atomic displacements that occur closer to the boundary carbon. This carbon atom, in the case of the 3-deacetylcalonectrin molecule, is terminal in R2, and through it the fused rings are attached.
The 1044 cm −1 mode has a significant contribution of the v(O16-C8) bond stretching vibration with the O atom within the hydroxyl group. Calculations predict a similar vibrational mode with a frequency of approximately 1076 cm −1 in the cyclopentanol (see Figure S3a).
Slightly weaker IR-active modes with frequencies of 1024, 1032 and 1061 cm −1 are predicted for the IR spectrum of 3-deacetylcalonectrin. They involve different stretching vibrations localized in the C-O bonds of the B-ring. A vibrational mode similar to the 1032 cm −1 mode is observed in the cyclopentanol model. It has atomic displacements localized in the C-O(H) bond and a close frequency of 1054 cm −1 (see Figure S3b).
The experimental absorption band at 974 cm −1 can be assigned to the mode with a frequency of 987 cm −1 in the calculated spectrum. In this mode, the atomic displacements occur mainly within the A-ring, with a characteristic v(C28C20) bond stretching vibration.
It is also worth noting that, in the calculated spectrum, an additional three modes are predicted with frequencies of 920, 955 and 962 cm −1 . Their IR activity is comparable with that of the 987 cm −1 mode. The 962 cm −1 mode involves contributions of the v(C33O34), v(C18O34) and δ(C18O33C34) vibrations, localized in the epoxide ring. Both 920 and 955 cm −1 modes include contributions of the v(C-O) vibrations localized in the B-ring. The 920 cm −1 mode also contains significant contributions of the v(C41O40) vibration in the acetate group and the v(C29C19) and v(C34O33) vibrations in the epoxide group.
Modes in the region of 630-882 cm −1 include significant contributions of the v(C-C) and v(C-O) vibrations, delocalized over the fused rings accompanied by the rocking modes. These modes have low polarity and therefore low activity in the predicted IR spectrum.
Additional contributions from δ(CCC), δ(COC), δ(OCC) and δ(OCO) bending modes start to occur. These contributions become much more significant in the region of 630-300 cm −1 , in which the ρ(HCH) contributions are also present.
Predominantly torsion modes are present below 300 cm −1 . The most active ones in this region are those localized in the CCOH groups, i.e., of the tors(CCOH) type. The cyclopentanol model predicts that the τ(OH) torsional mode is located at 280 cm −1 . This is close to the analogous vibrational mode (256 cm −1 ) in the calculated spectra of 3-deacetylcalonectrin. The relatively lower frequency of the mode in 3-deacetylcalonectrin can be explained by mixing with other torsional modes. Furthermore, some contributions of τ(OH) torsional vibrations can be found for the 370 and 388 cm −1 modes in the calculated spectrum of 3-deacetylcalonectrin.

T-2 Toxin Molecule
As mentioned previously, the T-2 toxin molecule has several structural similarities with 3-deacetylcalonectrin and the model structures shown in Figure S1. The one common feature of the T-2 toxin is the presence of the 3-methylbutanoate group. In order to reveal a manifestation of structural similarities, it was decided to compare the predicted IR absorbance spectrum of the T-2 toxin with that predicted for methyl 3-methylbutanoate (see Figure S4). The predicted spectra for the T-2 toxin, 3-deacetylcalonectrin, ethyl acetate and methyl 3-methylbutanoate are compared in Figure 4. Table 2 summarizes the most IR-intensive vibrational modes. According to theory, the acetate group attached to the C-ring has a low influence on v(OH) vibration. The corresponding mode in the calculated scaled IR absorbance spectrum is located at 3710 cm −1 . As in the case of 3-deacetylcalonectrin, the anharmonicity of the X-H bond stretching vibrations may also take place. Moreover, in [7], the short contact between the OH group of one molecule and the oxygen atom of another was highlighted.
For the T-2 toxin, various carbon-hydrogen stretching vibrations are located in the 2944-3116 cm −1 range, close to the similar region in 3-deacetylcalonectrin. For both substances, the frequencies of these modes are also overestimated due to neglect of anharmonicity. Table 2. Assignment of selected vibrational modes of the T-2 toxin molecule for calculated IR absorbance spectra. In the table, the following notations have been used: v-stretching vibrations, δ-bending vibrations, w-wagging, tors-torsional, τ-twisting, ρ-rocking modes, asymasymmetric, sym-symmetric vibrations. In the case of water-like >CH 2 structures, v sym and v asym correspond to symmetric and asymmetric water-like stretching vibrations, with corresponding atoms in parentheses. The upper index denotes the atom number and lower index denotes the number of elements attached to the previous atom.   The stretching C-H vibrational mode in the epoxy group possesses the highest frequency in this region. Its value is higher than in the case of the stretching vibrations of hydrogen bonded to sp 2 carbon (v(=C-H)). As mentioned above, the relation (1) is valid for the majority of vibrational C-H stretching modes. However, the relation can be violated due to weak coupling of modes with the boundary frequency values. This takes place for the modes localized in the R1 and R3 groups of the T-2 toxin molecule. The predominantly asymmetric vibrations in CH 2 and CH 3 are located in the 3017-3170 cm −1 range.

Mode
In the 2944-3116 cm −1 range the main difference between the calculated IR spectra of the two types of toxin concerns the resulting band shape. In the case of the T-2 toxin, the maximal absorption is related to the vibrational modes localized in the 3-methylbutanoate group with antisymmetric stretching atomic displacements. These modes are located near 3029, 3033 and 3043 cm −1 . Moreover, in the region of symmetric stretching vibrations among the three vibrational modes with the greatest absorption, two of them related the atomic displacements in the 3-methylbutanoate group, and the other one occurs due to symmetric vibrations in the CH 3 group attached to the A-ring. The presence of similar highly active vibrational modes in the calculated IR absorbance spectrum of methyl 3methylbutanoate leads to similar contour and maxima in that region (see Figure 6). However, in the 2944-3116 сm −1 region, the contribution from other C-H bonds of the T-2 toxin should not be forgotten. In this region, the stretching vibrations v(CH) with carbon atoms situated at the A-and B-ring boundaries possess the lowest frequency at 2944 cm −1 .
The theory predicts three intense IR absorption peaks near 1760 cm −1 . These originate from the v(C=O) vibrations. Such vibration localized in the R2 group corresponds to the peak located at 1775 cm −1 and has the highest IR-intensity. A similar feature group is present in the spectrum of 3-deacetylcalonectrin close at 1772 cm −1 . The v(C=O) vibration localized in the R3 group has a lower frequency and is located near to 1762 cm −1 . The v(C=O) vibration with the lowest frequency (1751 cm −1 ) is localized in the 3-methylbutanoate group.
Analysis of the v(C=O) modes in the ethyl acetate and methyl 3-methylbutanoate model molecules provides frequencies of 1772 and 1764 cm −1 , respectively. This is in qualitative agreement with the situation in the T-2 toxin molecule in which the v(C=O) mode However, in the 2944-3116 cm −1 region, the contribution from other C-H bonds of the T-2 toxin should not be forgotten. In this region, the stretching vibrations v(CH) with carbon atoms situated at the A-and B-ring boundaries possess the lowest frequency at 2944 cm −1 .
The theory predicts three intense IR absorption peaks near 1760 cm −1 . These originate from the v(C=O) vibrations. Such vibration localized in the R2 group corresponds to the peak located at 1775 cm −1 and has the highest IR-intensity. A similar feature group is present in the spectrum of 3-deacetylcalonectrin close at 1772 cm −1 . The v(C=O) vibration localized in the R3 group has a lower frequency and is located near to 1762 cm −1 . The v(C=O) vibration with the lowest frequency (1751 cm −1 ) is localized in the 3-methylbutanoate group.
Analysis of the v(C=O) modes in the ethyl acetate and methyl 3-methylbutanoate model molecules provides frequencies of 1772 and 1764 cm −1 , respectively. This is in qualitative agreement with the situation in the T-2 toxin molecule in which the v(C=O) mode localized in the R1 group has a lower frequency compared to those localized in the R2 and R3 groups.
Considering the T-2 toxin, the theory also predicts a double carbon stretching vibrational mode located at 1710 cm −1 .
The modes with the main contribution of the δ(HCH) and δ as (CH 3 ) deformational vibrations are located in the 1436-1509 cm −1 range. They give rise to a complex spectral feature, with three well-pronounced maxima at 1453, 1471 and 1490 cm −1 . The modes which significantly contribute to IR-absorbance in this region are listed in Table 2.
The region of 1362-1418 cm −1 in the calculated IR absorption spectra exhibits two bands with maxima near 1372 and 1390 cm −1 (see Figure 6). The main contributions to the 1372 cm −1 band give two vibrational modes. The more IR-active one is located at 1372 cm −1 . It is attributed to δ s (CH 3 ) vibrations in the acetate groups. The second mode (less IR-active) is located near 1376 cm −1 . It is assigned for various deformation vibrations with hydrogen atoms of δ(CCH) type. These δ(CCH) vibrations are distributed within the T-2 toxin molecule; several are present in the fused rings and other ones are present in the 3-methylbutanoate group. A similar vibration located at 1373 cm −1 can be found in the spectrum of the methyl 3-methylbutanoate molecule (see Table S3).
The 1390 cm −1 band in the spectrum of the T-2 toxin consists of two main contributions. The first of these is vibrational mode, with a frequency of 1391 cm −1 . Its predominant atomic displacements are δ s (C 12 H 3 ) and δ s (C 27 H 3 ), and also δ(O 16 C 8 H 9 ) and δ(C 8 O 16 H 17 ). The second contribution is the vibrational mode located at 1387 cm −1 . This is predominantly associated with δ(C 8 O 16 H 17 ), δ(O 16 C 8 H 9 ) and δ(O 18 C 10 H 11 ) vibrations.
The theory predicts for the T-2 toxin, a characteristic IR-peak near 1300 cm −1 , which is absent in the spectrum of 3-deacetylcalonectrin. It is related to the mode localized predominantly in the 3-methylbutanoate group and includes the w(H 51 C 49 H 50 ), δ(C 53 C 52 H 55 ), v(C 47 C 49 ) and v(C 47 O 46 ) vibrations. Comparison with the spectrum of the methyl 3methylbutanoate model supports this hypothesis. In the calculated IR absorbance spectrum of the model molecule, there is a 1302 cm −1 mode which has a similar atomic pattern (see Table S3 and Figure 6).
The most intensive band in the theoretical IR absorbance spectrum of the T-2 toxin is located at 1237 cm −1 . It is related to three modes with frequencies of 1239, 1236 and 1230 cm −1 . Their characteristic feature is the prominent contribution of v(CO) in both acetate groups. In the case of the 1239 cm −1 vibrational mode, the predominant v(CO) occurs in the R2. A similar mode exists in the spectrum of 3-deacetylcalonectrin. Additionally, in this mode, τ(H 51 C 49 H 50 ) atomic displacements take place. The 1230 cm −1 mode also has prominent v(CO) as well as (δ(C 40 C 42 H 45 ) and δ(C 8 O 16 H 17 )), τ(H 50 C 49 H 51 ) and tors(H 7 C 6 C 8 H 9 ) contributions. In the case of 1236 cm −1 , the predominant contribution is from v(CO) in the acetate groups of R3 as well as τ(H 50 C 49 H 51 ), with a smaller addition of various bending deformational modes with hydrogen atoms. The calculated dipole derivative unit vector in such vibrational modes has a direction close to the CO bond in which predominant displacements occur. Hence, these modes were interpreted as predominantly v(CO) absorption in the R2 and R3 groups (not in fused rings and not in the 3-methylbutanoate group). This is supported by the presence of a highly IR active mode in the spectrum of ethyl acetate and at the same time the absence of IR active modes is in the same region in the spectrum of methyl 3-methylbutanoate (see Figure 6).
Characteristic spectral features for the T-2 toxin are the bands in the 1090-1200 cm −1 region. The common feature of intense IR absorbance modes in this region is the significant contribution of stretching vibrations in the 3-methylbutanoate group, which include stretching vibrations v(C 47 O 46 ). The calculated frequencies of these modes are 1109, 1115, 1163, 1166 and 1183 cm −1 . The 1183 cm −1 mode in the T-2 toxin has similar atomic displacements in the 3-methylbutanoate groups to the 1196 cm −1 mode in methyl 3-methylbutanoate.
Additionally, there are significant contributions for the 1109 cm −1 mode from stretching vibrations in the A-ring and τ(HCH) in the epoxy ring; for the 1163 cm −1 mode, there are w(HCH) in the epoxy ring and v(C 12 C 3 ), and for the 1166 cm −1 mode there are w(HCH) between excited and ground state were performed. As a result, information on absorption for the T-2 toxin in the energy range from 2.6053 eV (475.90 nm) to 7.6618 eV (161.82 nm) and for 3-deacetylcalonectrin in the energy range from 2.4984 eV (496.25 nm) to 8.7945 eV (140.98 nm) was acquired.
The first excited state for T-2 and 3-deacetylcalonectrin possesses the longest wavelength and corresponds to a linear combination of transitions between orbitals, with the maximal contribution from the transition between Highest Occupied Molecular Orbital (HOMO) and Lowest Unoccupied Molecular Orbital (LUMO) (see Figures 7 and 8). This contribution is approximately 87% for the T-2 toxin and 75% for 3-deacetylcalonectrin. Nevertheless, the oscillator strength of this transition for both mycotoxins is extremely small and equals 0.0009 and 0.0005 for T-2 and 3-deacetylcalonectrin, respectively. It should be noted that, for most of the calculated transitions, the oscillator strength does not exceed 0.071 and 0.047 for T-2 and 3-deacetylcalonectrin, respectively.
The most intense 3-deacetylcalonectrin calculated transitions is 6.7390 eV (183.98 nm); its oscillator strength is 0.6168, which is more than an order of magnitude higher than for other transitions. The most pronounced contribution is related to the transition between orbitals 82 (HOMO−1) and 85 (LUMO+1) ( Table 3). The highest change in the transition between these orbitals is associated with the electron density in the double carbon bond. The π bonding character is noted for the 82 orbital, and the π* antibonding character is inherent to the 85 orbital. Thus, the discussed transition can be attributed as a π-π* transition. The 3-deacetylcalonectrin absorption spectrum is shown in Figure 7 (bottom).
The most intense transition among those calculated for the T-2 toxin is 6.4710 eV (191.60 nm); its oscillator strength, equal to 0.4432, is more than six times higher than the maximum oscillator strength for other transitions. In this case, the maximum contribution was achieved due to transitions from orbitals 124 (HOMO−1), 123 (HOMO−2) and 121 (HOMO−4) to orbital 127 (LUMO+1) (see Table 3 and Figure 8).
A characteristic similarity of the 124, 123 and 121 orbitals is the bonding nature of the electron density in the region of the C=C bond. In the 124 orbital, which has the greatest contribution to the transition, the most delocalized bonding region extends along the atoms from C=C in the A-ring, through C-O-C in the B-ring and ends at the carbon atom in the C-ring bound to the OH group. The bonding pattern at the double carbon bond corresponds to π bonding. During the transition to the 127 orbital at C=C, a distribution of electron density is noted similar to π* antibonding.
The transition with the highest oscillator strength in the T-2 toxin is similar to that in 3-deacetylcalonectrine in terms of the change in the C=C bond electron density. The main differences for the T-2 toxin are linked to a greater degree of delocalization within the trichothecene rings, as well as the presence of certain electron densities at the 3-methylbutanoate group, which results in a longer wavelength shift in absorption.

Conclusions
The molecular structure and vibrational spectra of 3-deacetylcalonectrin and the T-2 toxin were studied by computer modelling based on DFT calculation. Geometry optimization of 3-deacetylcalonectrin and the T-2 toxin provided stable configurations that were in good agreement in similar structural parts with experimental XRD data. It is determined that core units in both molecules have close structural parameters (bond length, angles). This supports the idea of the rigidity of fused rings. Several conformer configurations differing in orientations of the attached functional groups were found for the T-2 molecule.
The simulated IR spectrum of 3-deacetylcalonectrin is in good agreement with the experimental data. This result made it possible to assign all significant spectral feature to particular vibrational modes localized in different structural units. It was noticed that the most intense IR bands are associated with vibrations localized in the C-O and C=O bonds of the methyl acetate and acetate radicals. The modes with vibrations in the C-O and C=O bonds in the 3-methylbutanoate group have lower IR activity. Based on the computational results, it was possible to predict several spectral regions (approximately 1160-1185 cm −1 and 1300 cm −1 ) in which the characteristic vibrational modes of this group are expected.
The modelling of optical absorption spectra of 3-deacetylcalonectrin and T-2 toxins predicts only one significant π-π* transition associated with the electron density change localized in the C=C bond. The transition wavelength calculated for the T-2 toxin is close to the experimental value reported in several papers.