Tetraenone A: A New β-Ionone Derivative from Tetraena aegyptia

In this study, the chemical investigation of Tetraena aegyptia (Zygophyllaceae) led to the identification of a new megastigmene derivative, tetraenone A ((2S, 5R, 6R, 7E)-2-hydroxy-5,6-dihydro-β-ionone) (1), along with (3S, 5R, 6S, 7E)-3-hydroxy-5,6-epoxy-5,6-dihydro-β-ionone- (2), 3,4-dihydroxy-cinnamyl alcohol-4-glucoside (3), 3β,19α-dihydroxy-ursan-28-oic acid (4), quinovic acid (5), p-coumaric acid (6), and ferulic acid (7), for the first time. The chemical structures of 1–7 were confirmed by analysis of their 1D and 2D NMR and HRESIMS spectra and by their comparison with the relevant literature. The absolute configurations of 1 and 2 were assigned based on NOESY interactions and ECD spectra. Conformational analysis showed that 1 existed exclusively in one of the two theoretically possible chair conformers with a predominant s-trans configuration for the 3-oxobut-1-en-1-yl group with the ring, while the half-chair conformer had a pseudo-axial hydroxy group that was predominant over the other half-chair conformation. Boat conformations were not among the most stable conformations, and the s-trans isomerism was in favor of s-cis configuration. In silico investigation revealed that 1 and 2 had more favorable binding interactions with Mpro rather than with TMPRSS2. Accordingly, molecular dynamic simulations were performed on the complexes of compounds 1 and 2 with Mpro to explore the stability of their interaction with the target protein structure. Compounds 1 and 2 might offer a possible starting point for developing covalent inhibitors of Mpro of SARS-CoV-2.


Introduction
The genus Tetraena belongs to the caltrop family (Zygophyllaceae), which is distributed in the deserts and salt marshes of Egypt and is represented by nine species, including Tetraena aegyptia [1,2].It has been used in traditional medicine for a long time to treat various ailments, such as gout and rheumatism [1].Many studies have documented several biological activities for Tetraena aegyptia, such as anti-inflammatory, antidiabetic, and antitumor activities [3][4][5][6][7][8].Chemical studies on Tetraena aegyptia have revealed that quinovic acid and flavonoid glycosides are the main constituents in addition to a sulfonyl epoxy-lignan, which has been recently reported [6][7][8].
In the present chemical study on Tetraena aegyptia, seven compounds, including two dihydro-β-ionones (1 and 2), two triterpenes (4 and 5), and three phenolic compounds (3 and 6-7), were separated (Figure 1) for the first time from the plant.The absolute configurations of 1 and 2 were determined using circular dichroism (CD).CD is a powerful technique used for assigning the stereochemistry of organic compounds.The CD value corresponds to the difference in the absorption of left and right circularly polarized light by a chiral molecule and is expressed in ∆ε (ε l − ε r , differential molar extinction coefficient) [9,10].In this paper, the absolute configurations were obtained by comparing the experimental CD spectra with the theoretical CD spectra, and the conformational analysis was also studied.In 2022, a US patent (US2022/0175868) claimed the usefulness of the extract of the Middle Eastern herb black calla lily (Arum palaestinum) as a treatment for COVID-19, which is a viral disease caused by SARS-CoV-2 infection.This patent showed that the extract contains several phytocompounds, including β-ionone epoxide, compound 2.However, no further information was provided on the possible molecular target.The accumulated research on COVID-19 has demonstrated that two prominent proteins might serve as potential targets to develop possible treatments for COVID-19: the SARS-CoV-2 main protease, which is a viral protein involved in the replication of the virus, and the transmembrane serine protease 2 (TMPRSS2), which is a human protein responsible for viral entry into the host cells and, hence, viral infection.In light of this information, an in silico study of the isolated β-ionone-derived compounds 1 and 2 was conducted against these two molecular targets.M pro , a cysteine protease of SARS-CoV-2, is a key target for developing anti-SARS-CoV-2 therapeutics for the treatment of COVID-19.Similar to other cysteine proteases, the active site incorporates catalytic dyad residues, Cys145/His41, which provide a basis for the development of inhibitors.As of September 2021, the pharmaceutical giant Pfizer has initiated phase II/III clinical trials of PF-07321332, an orally administered peptidomimetic SARS-CoV-2 M pro inhibitor, which was later approved as Nirmatrelvir.

Conformational Analysis of Compounds 1 and 2
Conformational analyses of the forms of compounds 1 and 2 were conducted by adopting the systematic conformational search method using the Merk molecular force field (MMFF).The results are shown in the following sections.

Conformational Analysis of Compound 1
Compound 1 has a 1,1-dimethylcyclohexane scaffold bearing a 3-oxobut-1-en-1-yl group at the 6-position in a trans relationship with the hydroxy group at the 2-position and the methyl group at the 5-position.Conformational analysis showed two major conformers that constitute more than 70% of the calculated conformers' population (0.705 cumulative Boltzmann weight) and, together with the other three conformers, constitute more than 95% of the calculated conformers' population (0.959 cumulative Boltzmann weight).Theoretically, the six-membered cyclohexane ring can show two minima, namely chair and twist-boat conformations, and two maxima, namely boat and half-chair conformations.However, the conformational analysis showed that the five most energetically favorable conformers for compound 1 exclusively adopt chair conformations (relative energy difference within 7.17 kJ/mol).In addition, the hydroxy group at the 2-position, the 3-oxobut-1-en-1-yl group at the 6-position, and the methyl group at the 5-position of all the five conformers are exclusively in axial, equatorial, and equatorial positions, respectively.This indicates the absence of ring flipping to accommodate the other theoretically possible chair conformation in which the hydroxy group, the 3-oxobut-1-en-1-yl, and the methyl group at 5-position might accommodate equatorial, axial, and axial positions, respectively.Accordingly, the different conformers obtained arose from the different s-cis and s-trans conformations arising from rotation around sigma-bonds in the 3-oxobut-1-en-1-yl group as well as the position of the hydroxyl group's hydrogen atom.Interestingly, the s-trans relationship between the olefinic double bond and the cyclohexane ring was predominant in the most energetically favorable four conformers, constituting 93.6% of the conformers' population.Meanwhile, the s-cis relationship between the olefinic double bond and the cyclohexane ring existed in only the fifth most energetically favorable conformer, which constitutes only 2.3% of the conformers' populations.In contrast to the strong preference of the s-trans conformer considering the olefinic double bond and the cyclohexane ring, the s-cis/trans isomerization of the olefinic double bond and the carbonyl double bond was less impactful and was encountered with almost equal probabilities in the identified most energetically favorable conformers.On the other side, the position of the hydrogen atom of the hydroxyl group was clearly found in the most energetically favorable five conformers.This was indicated by three-quarters of the conformers' population having an almost 180 • dihedral angle (i.e., antiperiplanar), involving H-O-C 2 -C 1 atoms (conformers 1, 2, and 5), while the dihedral angle was almost 80 • (i.e., synclinal) in the remaining one-quarter of the conformers' population (Table 1).Compound 2 has a cyclohexene oxide (AKA 7-oxabicyclo [4.1.0]heptane)scaffold bearing geminal dimethyl substituents at the 1-position as well as a 3-oxobut-1-en-1-yl group at the 6-position in a cis relationship with the hydroxy group at the 3-position and the methyl group at the 5-position.Accordingly, the oxide group would be trans-configured to the 3-oxobut-1-en-1-yl, hydroxy, and 5-methyl groups.In contrast to compound 1's cyclohexane scaffold, which showed a limited number of conformers, the condensation of the small oxirane ring had pronounced effects on the conformational properties of the six-membered ring, resulting in a significantly larger number of conformers for  Compound 2 has a cyclohexene oxide (AKA 7-oxabicyclo [4.1.0]heptane)scaffold bearing geminal dimethyl substituents at the 1-position as well as a 3-oxobut-1-en-1-yl group at the 6-position in a cis relationship with the hydroxy group at the 3-position and the methyl group at the 5-position.Accordingly, the oxide group would be trans-configured to the 3-oxobut-1-en-1-yl, hydroxy, and 5-methyl groups.In contrast to compound 1's cyclohexane scaffold, which showed a limited number of conformers, the condensation of the small oxirane ring had pronounced effects on the conformational properties of the six-membered ring, resulting in a significantly larger number of conformers for  Compound 2 has a cyclohexene oxide (AKA 7-oxabicyclo [4.1.0]heptane)scaffold bearing geminal dimethyl substituents at the 1-position as well as a 3-oxobut-1-en-1-yl group at the 6-position in a cis relationship with the hydroxy group at the 3-position and the methyl group at the 5-position.Accordingly, the oxide group would be trans-configured to the 3-oxobut-1-en-1-yl, hydroxy, and 5-methyl groups.In contrast to compound 1's cyclohexane scaffold, which showed a limited number of conformers, the condensation of the small oxirane ring had pronounced effects on the conformational properties of the six-membered ring, resulting in a significantly larger number of conformers for  Compound 2 has a cyclohexene oxide (AKA 7-oxabicyclo [4.1.0]heptane)scaffold bearing geminal dimethyl substituents at the 1-position as well as a 3-oxobut-1-en-1-yl group at the 6-position in a cis relationship with the hydroxy group at the 3-position and the methyl group at the 5-position.Accordingly, the oxide group would be trans-configured to the 3-oxobut-1-en-1-yl, hydroxy, and 5-methyl groups.In contrast to compound 1's cyclohexane scaffold, which showed a limited number of conformers, the condensation of the small oxirane ring had pronounced effects on the conformational properties of the six-membered ring, resulting in a significantly larger number of conformers for  Compound 2 has a cyclohexene oxide (AKA 7-oxabicyclo [4.1.0]heptane)scaffold bearing geminal dimethyl substituents at the 1-position as well as a 3-oxobut-1-en-1-yl group at the 6-position in a cis relationship with the hydroxy group at the 3-position and the methyl group at the 5-position.Accordingly, the oxide group would be trans-configured to the 3-oxobut-1-en-1-yl, hydroxy, and 5-methyl groups.In contrast to compound 1's cyclohexane scaffold, which showed a limited number of conformers, the condensation of the small oxirane ring had pronounced effects on the conformational properties of the six-membered ring, resulting in a significantly larger number of conformers for  Compound 2 has a cyclohexene oxide (AKA 7-oxabicyclo [4.1.0]heptane)scaffold bearing geminal dimethyl substituents at the 1-position as well as a 3-oxobut-1-en-1-yl group at the 6-position in a cis relationship with the hydroxy group at the 3-position and the methyl group at the 5-position.Accordingly, the oxide group would be trans-configured to the 3-oxobut-1-en-1-yl, hydroxy, and 5-methyl groups.In contrast to compound 1's cyclohexane scaffold, which showed a limited number of conformers, the condensation of the small oxirane ring had pronounced effects on the conformational properties of the six-membered ring, resulting in a significantly larger number of conformers for  Compound 2 has a cyclohexene oxide (AKA 7-oxabicyclo [4.1.0]heptane)scaffold bearing geminal dimethyl substituents at the 1-position as well as a 3-oxobut-1-en-1-yl group at the 6-position in a cis relationship with the hydroxy group at the 3-position and the methyl group at the 5-position.Accordingly, the oxide group would be trans-configured to the 3-oxobut-1-en-1-yl, hydroxy, and 5-methyl groups.In contrast to compound 1's cyclohexane scaffold, which showed a limited number of conformers, the condensation of the small oxirane ring had pronounced effects on the conformational properties of the six-membered ring, resulting in a significantly larger number of conformers for  Compound 2 has a cyclohexene oxide (AKA 7-oxabicyclo [4.1.0]heptane)scaffold bearing geminal dimethyl substituents at the 1-position as well as a 3-oxobut-1-en-1-yl group at the 6-position in a cis relationship with the hydroxy group at the 3-position and the methyl group at the 5-position.Accordingly, the oxide group would be trans-configured to the 3-oxobut-1-en-1-yl, hydroxy, and 5-methyl groups.In contrast to compound 1's cyclohexane scaffold, which showed a limited number of conformers, the condensation of the small oxirane ring had pronounced effects on the conformational properties of the six-membered ring, resulting in a significantly larger number of conformers for  Compound 2 has a cyclohexene oxide (AKA 7-oxabicyclo [4.1.0]heptane)scaffold bearing geminal dimethyl substituents at the 1-position as well as a 3-oxobut-1-en-1-yl group at the 6-position in a cis relationship with the hydroxy group at the 3-position and the methyl group at the 5-position.Accordingly, the oxide group would be trans-configured to the 3-oxobut-1-en-1-yl, hydroxy, and 5-methyl groups.In contrast to compound 1's cyclohexane scaffold, which showed a limited number of conformers, the condensation of the small oxirane ring had pronounced effects on the conformational properties of the six-membered ring, resulting in a significantly larger number of conformers for  Compound 2 has a cyclohexene oxide (AKA 7-oxabicyclo [4.1.0]heptane)scaffold bearing geminal dimethyl substituents at the 1-position as well as a 3-oxobut-1-en-1-yl group at the 6-position in a cis relationship with the hydroxy group at the 3-position and the methyl group at the 5-position.Accordingly, the oxide group would be trans-configured to the 3-oxobut-1-en-1-yl, hydroxy, and 5-methyl groups.In contrast to compound 1's cyclohexane scaffold, which showed a limited number of conformers, the condensation of the small oxirane ring had pronounced effects on the conformational properties of the six-membered ring, resulting in a significantly larger number of conformers for  Compound 2 has a cyclohexene oxide (AKA 7-oxabicyclo [4.1.0]heptane)scaffold bearing geminal dimethyl substituents at the 1-position as well as a 3-oxobut-1-en-1-yl group at the 6-position in a cis relationship with the hydroxy group at the 3-position and the methyl group at the 5-position.Accordingly, the oxide group would be trans-configured to the 3-oxobut-1-en-1-yl, hydroxy, and 5-methyl groups.In contrast to compound 1's cyclohexane scaffold, which showed a limited number of conformers, the condensation of the small oxirane ring had pronounced effects on the conformational properties of the six-membered ring, resulting in a significantly larger number of conformers for  Compound 2 has a cyclohexene oxide (AKA 7-oxabicyclo [4.1.0]heptane)scaffold bearing geminal dimethyl substituents at the 1-position as well as a 3-oxobut-1-en-1-yl group at the 6-position in a cis relationship with the hydroxy group at the 3-position and the methyl group at the 5-position.Accordingly, the oxide group would be trans-configured to the 3-oxobut-1-en-1-yl, hydroxy, and 5-methyl groups.In contrast to compound 1's cyclohexane scaffold, which showed a limited number of conformers, the condensation of the small oxirane ring had pronounced effects on the conformational properties of the six-membered ring, resulting in a significantly larger number of conformers for

Conformational Analysis of Compound 2
Compound 2 has a cyclohexene oxide (AKA 7-oxabicyclo [4.1.0]heptane)scaffold bearing geminal dimethyl substituents at the 1-position as well as a 3-oxobut-1-en-1-yl group at the 6-position in a cis relationship with the hydroxy group at the 3-position and the methyl group at the 5-position.Accordingly, the oxide group would be trans-configured to the 3oxobut-1-en-1-yl, hydroxy, and 5-methyl groups.In contrast to compound 1's cyclohexane scaffold, which showed a limited number of conformers, the condensation of the small oxirane ring had pronounced effects on the conformational properties of the six-membered ring, resulting in a significantly larger number of conformers for compound 2. Referring to the six-membered cyclohexane ring, two interconverting half-chair conformations are theoretically possible in addition to endo/exo-boat conformations [25].Despite the previous studies suggesting the cyclohexene oxide to have boat conformations in equilibria with the most stable half-chair conformations [25,26], only half-chair conformations could be detected in the retrieved twelve most energetically favorable conformers (relative energy difference within 7.50 kJ/mol).Together, these twelve conformers constitute 96.3% of the conformers' population (cumulative Boltzmann weights of 0.963).However, 52.0% of the conformers' population exists in the form of only two conformers (conformers 1 and 2).Notably, all of the retrieved twelve conformers accommodate the 3-oxobut-1-en-1-yl, and the 5-methyl groups in pseudo-equatorial positions.Meanwhile, the hydroxyl group accommodated a pseudo-axial in nine conformers constituting 81.5% of the conformers' population (combined Boltzmann weights of 0.815 for conformers 1-3 and 6-11) and pseudo-equatorial in only three conformers constituting 14.8% of the conformers' population (combined Boltzmann weights of 0.148 for conformers 4, 5, and 12) [27,28].
This indicates that ring flipping equilibrium is more in favor of the half-chair conformation having the hydroxy group as pseudo-axial.In comparison with compound 1, the position of the hydrogen of the hydroxyl group was less influential, showing 60% of the retrieved conformers (  4 -atoms, respectively.While the s-cis/trans relationship (arising from rotation around sigma-bonds in the 3-oxobut-1-en-1-yl group) between the olefinic double bond and the ring was much more influential than the s-cis/trans relationship between the olefinic and the carbonyl double bonds, both have a potential impact in compound 2. Thus, at least 83.2% of the conformers' population (conformers 1, 2, 3, 4, 5, 6, 9, and 12) have the olefinic double bond and the oxirane ring in s-cis relationship compared with at least 13.1% of the conformers' population (conformers 7, 8, 10, and 11) having them in s-trans relationship.Meanwhile, at least 75.9% of the conformers' population (conformers 1, 2, 4, 5, 7, 9, 11, and 12) have the olefinic and carbonyl double bonds in s-cis relationship compared with at least 20.4% of the conformers' population (conformers 3, 6, 8, and 10) having them in s-trans relationship.Together, these results indicate the preference of compound 2 to accommodate the half-chair conformation in which the hydroxyl group accommodates the pseudo-axial position and has s-cis relationships for the olefinic double bond/oxirane ring and the olefinic/carbonyl double bonds (Table 2).2).   4 -atoms, respectively.While the scis/trans relationship (arising from rotation around sigma-bonds in the 3-oxobut-1-en-1-yl group) between the olefinic double bond and the ring was much more influential than the s-cis/trans relationship between the olefinic and the carbonyl double bonds, both have a potential impact in compound 2. Thus, at least 83.2% of the conformers' population (conformers 1, 2, 3, 4, 5, 6, 9, and 12) have the olefinic double bond and the oxirane ring in scis relationship compared with at least 13.1% of the conformers' population (conformers 7, 8, 10, and 11) having them in s-trans relationship.Meanwhile, at least 75.9% of the conformers' population (conformers 1, 2, 4, 5, 7, 9, 11, and 12) have the olefinic and carbonyl double bonds in s-cis relationship compared with at least 20.4% of the conformers' population (conformers 3, 6, 8, and 10) having them in s-trans relationship.Together, these results indicate the preference of compound 2 to accommodate the half-chair conformation in which the hydroxyl group accommodates the pseudo-axial position and has s-cis relationships for the olefinic double bond/oxirane ring and the olefinic/carbonyl double bonds (Table 2).compound 2. Referring to the six-membered cyclohexane ring, two interconverting halfchair conformations are theoretically possible in addition to endo/exo-boat conformations [25].Despite the previous studies suggesting the cyclohexene oxide to have boat conformations in equilibria with the most stable half-chair conformations [25,26], only half-chair conformations could be detected in the retrieved twelve most energetically favorable conformers (relative energy difference within 7.50 kJ/mol).Together, these twelve conformers constitute 96.3% of the conformers' population (cumulative Boltzmann weights of 0.963).However, 52.0% of the conformers' population exists in the form of only two conformers (conformers 1 and 2).Notably, all of the retrieved twelve conformers accommodate the 3oxobut-1-en-1-yl, and the 5-methyl groups in pseudo-equatorial positions.Meanwhile, the hydroxyl group accommodated a pseudo-axial in nine conformers constituting 81.5% of the conformers' population (combined Boltzmann weights of 0.815 for conformers 1-3 and 6-11) and pseudo-equatorial in only three conformers constituting 14.8% of the conformers' population (combined Boltzmann weights of 0.148 for conformers 4, 5, and 12) [27,28].This indicates that ring flipping equilibrium is more in favor of the half-chair conformation having the hydroxy group as pseudo-axial.In comparison with compound 1, the position of the hydrogen of the hydroxyl group was less influential, showing 60% of the retrieved conformers (1, 3, 4, 7, and 8) accommodating antiperiplanar/synclinal dihedral angles (almost 180/60°) for H-O-C 3 -C 4 and H-O-C 3 -C 2 -atoms, respectively, and 39% of the retrieved conformers (2, 5, 6, 10, and 11) have the antiperiplanar/synclinal dihedral angles (almost 180/60°) for H-O-C 3 -C 2 and H-O-C 3 -C 4 -atoms, respectively.While the scis/trans relationship (arising from rotation around sigma-bonds in the 3-oxobut-1-en-1-yl group) between the olefinic double bond and the ring was much more influential than the s-cis/trans relationship between the olefinic and the carbonyl double bonds, both have a potential impact in compound 2. Thus, at least 83.2% of the conformers' population (conformers 1, 2, 3, 4, 5, 6, 9, and 12) have the olefinic double bond and the oxirane ring in scis relationship compared with at least 13.1% of the conformers' population (conformers 7, 8, 10, and 11) having them in s-trans relationship.Meanwhile, at least 75.9% of the conformers' population (conformers 1, 2, 4, 5, 7, 9, 11, and 12) have the olefinic and carbonyl double bonds in s-cis relationship compared with at least 20.4% of the conformers' population (conformers 3, 6, 8, and 10) having them in s-trans relationship.Together, these results indicate the preference of compound 2 to accommodate the half-chair conformation in which the hydroxyl group accommodates the pseudo-axial position and has s-cis relationships for the olefinic double bond/oxirane ring and the olefinic/carbonyl double bonds (Table 2).compound 2. Referring to the six-membered cyclohexane ring, two interconverting halfchair conformations are theoretically possible in addition to endo/exo-boat conformations [25].Despite the previous studies suggesting the cyclohexene oxide to have boat conformations in equilibria with the most stable half-chair conformations [25,26], only half-chair conformations could be detected in the retrieved twelve most energetically favorable conformers (relative energy difference within 7.50 kJ/mol).Together, these twelve conformers constitute 96.3% of the conformers' population (cumulative Boltzmann weights of 0.963).However, 52.0% of the conformers' population exists in the form of only two conformers (conformers 1 and 2).Notably, all of the retrieved twelve conformers accommodate the 3oxobut-1-en-1-yl, and the 5-methyl groups in pseudo-equatorial positions.Meanwhile, the hydroxyl group accommodated a pseudo-axial in nine conformers constituting 81.5% of the conformers' population (combined Boltzmann weights of 0.815 for conformers 1-3 and 6-11) and pseudo-equatorial in only three conformers constituting 14.8% of the conformers' population (combined Boltzmann weights of 0.148 for conformers 4, 5, and 12) [27,28].This indicates that ring flipping equilibrium is more in favor of the half-chair conformation having the hydroxy group as pseudo-axial.In comparison with compound 1, the position of the hydrogen of the hydroxyl group was less influential, showing 60% of the retrieved conformers (1, 3, 4, 7, and 8) accommodating antiperiplanar/synclinal dihedral angles (almost 180/60°) for H-O-C 3 -C 4 and H-O-C 3 -C 2 -atoms, respectively, and 39% of the retrieved conformers (2, 5, 6, 10, and 11) have the antiperiplanar/synclinal dihedral angles (almost 180/60°) for H-O-C 3 -C 2 and H-O-C 3 -C 4 -atoms, respectively.While the scis/trans relationship (arising from rotation around sigma-bonds in the 3-oxobut-1-en-1-yl group) between the olefinic double bond and the ring was much more influential than the s-cis/trans relationship between the olefinic and the carbonyl double bonds, both have a potential impact in compound 2. Thus, at least 83.2% of the conformers' population (conformers 1, 2, 3, 4, 5, 6, 9, and 12) have the olefinic double bond and the oxirane ring in scis relationship compared with at least 13.1% of the conformers' population (conformers 7, 8, 10, and 11) having them in s-trans relationship.Meanwhile, at least 75.9% of the conformers' population (conformers 1, 2, 4, 5, 7, 9, 11, and 12) have the olefinic and carbonyl double bonds in s-cis relationship compared with at least 20.4% of the conformers' population (conformers 3, 6, 8, and 10) having them in s-trans relationship.Together, these results indicate the preference of compound 2 to accommodate the half-chair conformation in which the hydroxyl group accommodates the pseudo-axial position and has s-cis relationships for the olefinic double bond/oxirane ring and the olefinic/carbonyl double bonds (Table 2).double bonds in s-cis relationship compared with at least 20.4% of the conformers' population (conformers 3, 6, 8, and 10) having them in s-trans relationship.Together, these results indicate the preference of compound 2 to accommodate the half-chair conformation in which the hydroxyl group accommodates the pseudo-axial position and has s-cis relationships for the olefinic double bond/oxirane ring and the olefinic/carbonyl double bonds (Table 2).double bonds in s-cis relationship compared with at least 20.4% of the conformers' population (conformers 3, 6, 8, and 10) having them in s-trans relationship.Together, these results indicate the preference of compound 2 to accommodate the half-chair conformation in which the hydroxyl group accommodates the pseudo-axial position and has s-cis relationships for the olefinic double bond/oxirane ring and the olefinic/carbonyl double bonds (Table 2).double bonds in s-cis relationship compared with at least 20.4% of the conformers' population (conformers 3, 6, 8, and 10) having them in s-trans relationship.Together, these results indicate the preference of compound 2 to accommodate the half-chair conformation in which the hydroxyl group accommodates the pseudo-axial position and has s-cis relationships for the olefinic double bond/oxirane ring and the olefinic/carbonyl double bonds (Table 2).Electronic circular dichroism (ECD) spectra of (R)-and (S)-enantiomers of compounds 1 and 2 were calculated to assign their absolute configurations.As ab initio calculations are computationally demanding, a proper compromise between computational cost and accuracy is needed.Consequently, the time-dependent density functional theory (TDDFT) method was implemented for its satisfactory accuracy/computational load.The implementation of a proper combination of functional/basis is pivotal for TDDFT calculations.Benchmarking studies have indicated that the calculation of ECD spectra using range-separated hybrid functions such as CAM-B3LYP is more accurate than only hybrid Electronic circular dichroism (ECD) spectra of (R)-and (S)-enantiomers of compounds 1 and 2 were calculated to assign their absolute configurations.As ab initio calculations are computationally demanding, a proper compromise between computational cost and accuracy is needed.Consequently, the time-dependent density functional theory (TDDFT) method was implemented for its satisfactory accuracy/computational load.The implementation of a proper combination of functional/basis is pivotal for TDDFT calculations.Benchmarking studies have indicated that the calculation of ECD spectra using range-separated hybrid functions such as CAM-B3LYP is more accurate than only hybrid Electronic circular dichroism (ECD) spectra of (R)-and (S)-enantiomers of compounds 1 and 2 were calculated to assign their absolute configurations.As ab initio calculations are computationally demanding, a proper compromise between computational cost and accuracy is needed.Consequently, the time-dependent density functional theory (TDDFT) method was implemented for its satisfactory accuracy/computational load.The implementation of a proper combination of functional/basis is pivotal for TDDFT calculations.Benchmarking studies have indicated that the calculation of ECD spectra using range-separated hybrid functions such as CAM-B3LYP is more accurate than only hybrid Electronic circular dichroism (ECD) spectra of (R)-and (S)-enantiomers of compounds 1 and 2 were calculated to assign their absolute configurations.As ab initio calculations are computationally demanding, a proper compromise between computational cost and accuracy is needed.Consequently, the time-dependent density functional theory (TDDFT) method was implemented for its satisfactory accuracy/computational load.The implementation of a proper combination of functional/basis is pivotal for TDDFT calculations.Benchmarking studies have indicated that the calculation of ECD spectra using range-separated hybrid functions such as CAM-B3LYP is more accurate than only hybrid Electronic circular dichroism (ECD) spectra of (R)-and (S)-enantiomers of compounds 1 and 2 were calculated to assign their absolute configurations.As ab initio calculations are computationally demanding, a proper compromise between computational cost and accuracy is needed.Consequently, the time-dependent density functional theory (TDDFT) method was implemented for its satisfactory accuracy/computational load.The implementation of a proper combination of functional/basis is pivotal for TDDFT calculations.Benchmarking studies have indicated that the calculation of ECD spectra using range-separated hybrid functions such as CAM-B3LYP is more accurate than only hybrid Electronic circular dichroism (ECD) spectra of (R)-and (S)-enantiomers of compounds 1 and 2 were calculated to assign their absolute configurations.As ab initio calculations are computationally demanding, a proper compromise between computational cost and accuracy is needed.Consequently, the time-dependent density functional theory (TDDFT) method was implemented for its satisfactory accuracy/computational load.The implementation of a proper combination of functional/basis is pivotal for TDDFT calculations.Benchmarking studies have indicated that the calculation of ECD spectra using range-separated hybrid functions such as CAM-B3LYP is more accurate than only hybrid 2.2.3.Assignment of the Absolute Stereochemical Configuration of Compounds 1 and 2 (ECD Spectra) Electronic circular dichroism (ECD) spectra of (R)-and (S)-enantiomers of compounds 1 and 2 were calculated to assign their absolute configurations.As ab initio calculations are computationally demanding, a proper compromise between computational cost and accuracy is needed.Consequently, the time-dependent density functional theory (TDDFT) method was implemented for its satisfactory accuracy/computational load.The implementation of a proper combination of functional/basis is pivotal for TDDFT calculations.Benchmarking studies have indicated that the calculation of ECD spectra using range-separated hybrid functions such as CAM-B3LYP is more accurate than only hybrid functions such as B3LYP.Benchmarking studies also indicated more reliability for the implementation of Ahlrichs basis sets such as SVP for TDDFT calculations rather than split-valence Pople basis sets such as 6-31+G.Therefore, CAM-B3LYP/SVP as a functional/basis set combination was employed in calculating ECD spectra using a conductor-like polarizable continuum model (CPCM) as a solvent model.
The experimental CD spectrum of compound 1 (Figure 3) showed a positive Cotton effect at 219.5 nm and a negative Cotton effect at 208.13 nm.This experimental CD spectrum matched the calculated ECD spectrum for the 2S, 5R, 6R stereoisomer; thus, compound 1 was established as (2S, 5R, 6R, 7E)-2-hydroxy-5,6-dihydro-β-ionone.Meanwhile, the experimental CD spectrum of compound 2 showed a negative Cotton effect at 221.13 nm, which matched the calculated ECD spectrum for the 3R, 5S, 6R stereoisomer; thus, the absolute configuration of compound 2 was concluded as (3S, 5R, 6S, 7E)-3-hydroxy-5,6-epoxy-5,6dihydro-β-ionone (also named as (3S, 5R, 6S, 7E)-5,6-epoxy-3-hydroxy-7-megastigmene-9one).As a validated drug target, developing the small molecule inhibitors of M pro rather than the peptidomimetic inhibitor might offer multiple advantages in the fight against COVID-19.As is clear from the crystal structure of M pro (PDB: 6lu7), the nucleophilic key residue Cys145 can establish a covalent interaction with inhibitors compromising an electrophilic moiety.Considering that the structures of compounds 1 and 2 incorporate an α,β-unsaturated carbonyl moiety, which is a well-known electrophilic might undergo nucleophilic addition, it was interesting to explore whether these compounds can bind to M pro and establish a covalent interaction with Cys145.Therefore, in silico covalent docking was addressed to predict the abilities of compounds 1 and 2 to bind and inhibit M of SARS-CoV-2.As shown in Figure 4, compound 1 could successfully dock into the active site of M pro of SARS-CoV-2 with a calculated favorable docking score of −4.89957 kcal/mol.As a validated drug target, developing the small molecule inhibitors of M pro rather than the peptidomimetic inhibitor might offer multiple advantages in the fight against COVID-19.As is clear from the crystal structure of M pro (PDB: 6lu7), the nucleophilic key residue Cys145 can establish a covalent interaction with inhibitors compromising an electrophilic moiety.Considering that the structures of compounds 1 and 2 incorporate an α,β-unsaturated carbonyl moiety, which is a well-known electrophilic moiety that might undergo nucleophilic addition, it was interesting to explore whether these compounds can bind to M pro and establish a covalent interaction with Cys145.Therefore, in silico covalent docking was addressed to predict the abilities of compounds 1 and 2 to bind and inhibit M pro of SARS-CoV-2.As shown in Figure 4, compound 1 could successfully dock into the active site of M pro of SARS-CoV-2 with a calculated favorable docking score of −4.89957 kcal/mol.α,β-unsaturated carbonyl moiety, which is a well-known electrophilic moiety that might undergo nucleophilic addition, it was interesting to explore whether these compounds can bind to M pro and establish a covalent interaction with Cys145.Therefore, in silico covalent docking was addressed to predict the abilities of compounds 1 and 2 to bind and inhibit M pro of SARS-CoV-2.As shown in Figure 4, compound 1 could successfully dock into the active site of M pro of SARS-CoV-2 with a calculated favorable docking score of −4.89957 kcal/mol.In addition to establishing a covalent interaction with the catalytic cys145, it established a network of favorable interactions involving two hydrogen-donor bonding interactions between the compound 1 oxygen atom of the carbonyl group and the backbone NH of both of Cys145 and Gly143 and one hydrogen bond acceptor between the compound 1 OH group and the side chain carbonyl group of Asn142, as well as two hydrophobic π-alkyl interactions between both vicinal methyl groups of compound 1 and the imidazole ring of His163.Compound 2, as shown in Figure 5, also docked successfully into the active site of M pro of SARS-CoV-2 with an almost similar docking score to compound 1 (docking score of −4.85742 kcal/mol for compound 2 versus docking score of −4.89957 kcal/mol for compound 1).As illustrated in Figure 2, the oxygen atom of the carbonyl group of 2 still maintains the two hydrogen donor bonding interactions with the backbone NH of both Cys145 and Gly143.In addition to establishing a covalent interaction with the catalytic cys145, it established a network of favorable interactions involving two hydrogen-donor bonding interactions between the compound 1 oxygen atom of the carbonyl group and the backbone NH of both of Cys145 and Gly143 and one hydrogen bond acceptor between the compound 1 OH group and the side chain carbonyl group of Asn142, as well as two hydrophobic π-alkyl interactions between both vicinal methyl groups of compound 1 and the imidazole ring of His163.Compound 2, as shown in Figure 5, also docked successfully into the active site of M pro of SARS-CoV-2 with an almost similar docking score to compound 1 (docking score of −4.85742 kcal/mol for compound 2 versus docking score of −4.89957 kcal/mol for compound 1).As illustrated in Figure 2, the oxygen atom of the carbonyl group of 2 still maintains the two hydrogen donor bonding interactions with the backbone NH of both Cys145 and Gly143.However, only one methyl group of the two vicinal methyl groups of 2 is in a hydrophobic π-alkyl interaction with the imidazole ring of His41, while the methyl group vicinal to the epoxy moiety is in a hydrophobic π-alkyl interaction with imidazole ring of His163.Meanwhile, the OH group of 2 is involved in a hydrogen-bond donor interaction with the backbone NH of Glu166.The epoxy group of 2 did not contribute any significant binding interactions, and hence, it might be insignificant for activity.However, only one methyl group of the two vicinal methyl groups of 2 is in a hydrophobic π-alkyl interaction with the imidazole ring of His41, while the methyl group vicinal to the epoxy moiety is in a hydrophobic π-alkyl interaction with imidazole ring of His163.Meanwhile, the OH group of 2 is involved in a hydrogen-bond donor interaction with the backbone NH of Glu166.The epoxy group of 2 did not contribute any significant binding interactions, and hence, it might be insignificant for activity.

In Silico Evaluation of Compounds 1 and 2 against Transmembrane Serine Protease 2 (TMPRSS2)
TMPRSS2 is a member of the serine protease family that is highly expressed in nasal, bronchial, and gastrointestinal epithelial cells.It was found that TMPRSS2 is a key factor involved in viral entry and spread inside the human body.SARS-CoV-2 utilizes this enzyme to cleave its spike protein (S protein) and prime the virus for cell entry.Inhibitors of TMPRSS2 are expected to be promising treatments for SARS-CoV-2 infections.The substrate site of TMPRSS2 encompasses a catalytic triad that involves the amino acid residues Ser441, Asp435, and His296.As found in the crystal structure of TMPRSS2 (PDB ID: 7meq), covalent inhibitors undergo nucleophilic attack by Ser441 to establish a covalent bond, resulting in inactivation and inhibition of TMPRSS2.As mentioned above, the structures of compounds 1 and 2 incorporate an α,β-unsaturated carbonyl moiety, which is a well-known electrophilic moiety that might undergo nucleophilic addition.Accordingly, a covalent docking was carried out to explore the ability of compounds 1 and 2 to dock into the active site of TMPRSS2 and establish stable bonding interactions with it.Both compounds 1 and 2 docked successfully with docking scores of −4.1051 and −3.8134 kcal/mol, respectively.As illustrated in Figure 6, compound 1 was able to successfully form a covalent bonding interaction with Ser441.
respectively.As illustrated in Figure 6, compound 1 was able to successfully form a covalent bonding interaction with Ser441.In addition, it formed a favorable hydrogen-bonding acceptor interaction between the oxygen atom of the carbonyl moiety and the backbone NH of the Gly439 residue.Also, it formed a carbon-hydrogen interaction between the OH group and Cys437.Regarding compound 2, successful covalent docking was also established into the active site of the enzyme as shown in Figure 7, with a maintained covalent bond with Ser441 via its α,βunsaturated carbonyl warhead moiety.In addition, it formed a favorable hydrogen-bonding acceptor interaction between the oxygen atom of the carbonyl moiety and the backbone NH of the Gly439 residue.Also, it formed a carbon-hydrogen interaction between the OH group and Cys437.Regarding compound 2, successful covalent docking was also established into the active site of the enzyme as shown in Figure 7, with a maintained covalent bond with Ser441 via its α,βunsaturated carbonyl warhead moiety.
Additionally, it formed other favorable interactions with other amino acid residues inside the active site, where a hydrogen-bond donor interaction was formed between the oxygen atom of the carbonyl moiety, as in compound 1, but here with Ser441 of the catalytic triad of the active site.Moreover, the epoxy group of 2 contributed successfully to a hydrogen-bond donor interaction with Gly439.This interaction suggested that the epoxy functionality was of interesting value, regarding the interaction of compound 2, in contrast to compound 1.Furthermore, two favorable hydrophobic π-alkyl interactions were formed between His296 of the catalytic triad of the enzyme with each cyclohexane ring and one methyl group of the two vicinal methyl groups of 2.
In addition, it formed a favorable hydrogen-bonding acceptor interaction between the oxygen atom of the carbonyl moiety and the backbone NH of the Gly439 residue.Also, it formed a carbon-hydrogen interaction between the OH group and Cys437.Regarding compound 2, successful covalent docking was also established into the active site of the enzyme as shown in Figure 7, with a maintained covalent bond with Ser441 via its α,βunsaturated carbonyl warhead moiety.Additionally, it formed other favorable interactions with other amino acid residues inside the active site, where a hydrogen-bond donor interaction was formed between the oxygen atom of the carbonyl moiety, as in compound 1, but here with Ser441 of the catalytic triad of the active site.Moreover, the epoxy group of 2 contributed successfully to a hydrogen-bond donor interaction with Gly439.This interaction suggested that the epoxy functionality was of interesting value, regarding the interaction of compound 2, in contrast to compound 1.Furthermore, two favorable hydrophobic π-alkyl interactions were formed between His296 of the catalytic triad of the enzyme with each cyclohexane ring and one methyl group of the two vicinal methyl groups of 2.

Molecular Dynamic Simulations
A molecular dynamic simulation was performed to evaluate the stability of targetligand complexes for compounds 1 and 2 with their target protein.By examination of the docking scores of compounds 1 and 2, it was found that their docking scores with M pro were better than those in the case of docking with the TMPRSS2 enzyme.Accordingly, the best docking pose for each compound with M pro was subjected to molecular dynamic simulation over a period of 10 ns.The whole ligand-protein complex was utilized to evaluate the stability of the complex formed during the docking experiment.Regarding compound 1, the simulation revealed that there was an initial increase in the root mean square deviation (RMSD) value for the protein backbone until it reached about 4.5-5 Å in a period of about 2 ns (Figure 8).The RMSD value regarding the ligand of the same complex showed initial fluctuation until about 8 ns, and then it was slightly stabilized for the rest of the simulation time.Similarly, compound 2 was found to have almost the same pattern (Figure 8).The average number of hydrogen bonds formed between the ligand and target protein was calculated and plotted against time frames.Again, compounds 1 and 2 were found to form about the same average number of hydrogen bonds (Figure 9).Collectively, molecular dynamic simulation studies have revealed that both investigated compounds reach stability at about the same time, and they are able to form nearly the same average number of hydrogen bonds with their target, M pro of SARS-CoV2.

Molecular Dynamic Simulations
A molecular dynamic simulation was performed to evaluate the stability of targetligand complexes for compounds 1 and 2 with their target protein.By examination of the docking scores of compounds 1 and 2, it was found that their docking scores with M pro were better than those in the case of docking with the TMPRSS2 enzyme.Accordingly, the best docking pose for each compound with M pro was subjected to molecular dynamic simulation over a period of 10 ns.The whole ligand-protein complex was utilized to evaluate the stability of the complex formed during the docking experiment.Regarding compound 1, the simulation revealed that there was an initial increase in the root mean square deviation (RMSD) value for the protein backbone until it reached about 4.5-5 Å in a period of about 2 ns (Figure 8).The RMSD value regarding the ligand of the same complex showed initial fluctuation until about 8 ns, and then it was slightly stabilized for the rest of the simulation time.Similarly, compound 2 was found to have almost the same pattern (Figure 8).The average number of hydrogen bonds formed between the ligand and target protein was calculated and plotted against time frames.Again, compounds 1 and 2 were found to form about the same average number of hydrogen bonds (Figure 9).Collectively, molecular dynamic simulation studies have revealed that both investigated compounds reach stability at about the same time, and they are able to form nearly the same average number of hydrogen bonds with their target, M pro of SARS-CoV2.

Plant Material
Aerial parts of T. aegyptia were collected from the coastal desert in Egypt, and a voucher specimen (ZA-35-PD) was kept in the Department of Pharmacognosy, Faculty of Pharmacy, Mansoura University.

compound 2 .
Referring to the six-membered cyclohexane ring, two interconverting halfchair conformations are theoretically possible in addition to endo/exo-boat conformations[25].Despite the previous studies suggesting the cyclohexene oxide to have boat conformations in equilibria with the most stable half-chair conformations[25,26], only half-chair conformations could be detected in the retrieved twelve most energetically favorable conformers (relative energy difference within 7.50 kJ/mol).Together, these twelve conformers constitute 96.3% of the conformers' population (cumulative Boltzmann weights of 0.963).However, 52.0% of the conformers' population exists in the form of only two conformers (conformers 1 and 2).Notably, all of the retrieved twelve conformers accommodate the 3oxobut-1-en-1-yl, and the 5-methyl groups in pseudo-equatorial positions.Meanwhile, the hydroxyl group accommodated a pseudo-axial in nine conformers constituting 81.5% of the conformers' population (combined Boltzmann weights of 0.815 for conformers 1-3 and 6-11) and pseudo-equatorial in only three conformers constituting 14.8% of the conformers' population (combined Boltzmann weights of 0.148 for conformers 4, 5, and 12)[27,28].This indicates that ring flipping equilibrium is more in favor of the half-chair conformation having the hydroxy group as pseudo-axial.In comparison with compound 1, the position of the hydrogen of the hydroxyl group was less influential, showing 60% of the retrieved conformers (1, 3, 4, 7, and 8) accommodating antiperiplanar/synclinal dihedral angles (almost 180/60°) for H-O-C 3 -C 4 and H-O-C 3 -C 2 -atoms, respectively, and 39% of the retrieved conformers (2, 5, 6, 10, and 11) have the antiperiplanar/synclinal dihedral angles (almost 180/60°) for H-O-C 3 -C 2 and H-O-C 3 -C 4 -atoms, respectively.While the scis/trans relationship (arising from rotation around sigma-bonds in the 3-oxobut-1-en-1-yl group) between the olefinic double bond and the ring was much more influential than the s-cis/trans relationship between the olefinic and the carbonyl double bonds, both have a potential impact in compound 2. Thus, at least 83.2% of the conformers' population (conformers 1, 2, 3, 4, 5, 6, 9, and 12) have the olefinic double bond and the oxirane ring in scis relationship compared with at least 13.1% of the conformers' population (conformers 7, 8, 10, and 11) having them in s-trans relationship.Meanwhile, at least 75.9% of the conformers' population (conformers 1, 2, 4, 5, 7, 9, 11, and 12) have the olefinic and carbonyl double bonds in s-cis relationship compared with at least 20.4% of the conformers' population (conformers 3, 6, 8, and 10) having them in s-trans relationship.Together, these results indicate the preference of compound 2 to accommodate the half-chair conformation in which the hydroxyl group accommodates the pseudo-axial position and has s-cis relationships for the olefinic double bond/oxirane ring and the olefinic/carbonyl double bonds (Table

compound 2 .
Referring to the six-membered cyclohexane ring, two interconverting halfchair conformations are theoretically possible in addition to endo/exo-boat conformations[25].Despite the previous studies suggesting the cyclohexene oxide to have boat conformations in equilibria with the most stable half-chair conformations[25,26], only half-chair conformations could be detected in the retrieved twelve most energetically favorable conformers (relative energy difference within 7.50 kJ/mol).Together, these twelve conformers constitute 96.3% of the conformers' population (cumulative Boltzmann weights of 0.963).However, 52.0% of the conformers' population exists in the form of only two conformers (conformers 1 and 2).Notably, all of the retrieved twelve conformers accommodate the 3oxobut-1-en-1-yl, and the 5-methyl groups in pseudo-equatorial positions.Meanwhile, the hydroxyl group accommodated a pseudo-axial in nine conformers constituting 81.5% of the conformers' population (combined Boltzmann weights of 0.815 for conformers 1-3 and 6-11) and pseudo-equatorial in only three conformers constituting 14.8% of the conformers' population (combined Boltzmann weights of 0.148 for conformers 4, 5, and 12)[27,28].This indicates that ring flipping equilibrium is more in favor of the half-chair conformation having the hydroxy group as pseudo-axial.In comparison with compound 1, the position of the hydrogen of the hydroxyl group was less influential, showing 60% of the retrieved conformers(1, 3, 4, 7, and 8)  accommodating antiperiplanar/synclinal dihedral angles (almost 180/60°) for H-O-C 3 -C 4 and H-O-C 3 -C 2 -atoms, respectively, and 39% of the retrieved conformers (2, 5, 6, 10, and 11) have the antiperiplanar/synclinal dihedral angles (almost 180/60°) for H-O-C 3 -C 2 and H-O-C 3 -C

3 .
Assignment of the Absolute Stereochemical Configuration of Compounds 1 and 2 (ECD Spectra) Metabolites 2023, 13, x FOR PEER REVIEW 10 of 17

Figure 3 .
Figure 3. CD and ECD spectra of compounds 1 and 2.

Figure 3 .
Figure 3. CD and ECD spectra of compounds 1 and 2.

Figure 4 .
Figure 4. Compound 1 covalently docked into the active site of M pro of SARS-CoV-2.Figure 4. Compound 1 covalently docked into the active site of M pro of SARS-CoV-2.

Figure 4 .
Figure 4. Compound 1 covalently docked into the active site of M pro of SARS-CoV-2.Figure 4. Compound 1 covalently docked into the active site of M pro of SARS-CoV-2.

Figure 5 .
Figure 5. Compound 2 covalently docked into the active site of M pro of SARS-CoV-2.

Figure 5 .
Figure 5. Compound 2 covalently docked into the active site of M pro of SARS-CoV-2.

Figure 6 .
Figure 6.Compound 1 covalently docked into the active site of TMPRSS2.

Figure 6 .
Figure 6.Compound 1 covalently docked into the active site of TMPRSS2.

Figure 7 .
Figure 7. Compound 2 covalently docked into the active site of TMPRSS2.

Figure 7 .
Figure 7. Compound 2 covalently docked into the active site of TMPRSS2.

Figure 8 .
Figure 8.The ligand-protein RMSD plot of compounds 1 and 2 complexed with M pro of SARS-CoV2: (A) compound 1; (B) compound 2. Left and right Y-axes represent the RMSD values of protein and ligand in molecular distance unit Angstrom (Å).X-axis demonstrates the time in picoseconds (ps).

Figure 8 .
Figure 8.The ligand-protein RMSD plot of compounds 1 and 2 complexed with M pro of SARS-CoV2: (A) compound 1; (B) compound 2. Left and right Y-axes represent the RMSD values of protein and ligand in molecular distance unit Angstrom (Å).X-axis demonstrates the time in picoseconds (ps).

Figure 8 .
Figure 8.The ligand-protein RMSD plot of compounds 1 and 2 complexed with M pro of SARS-CoV2: (A) compound 1; (B) compound 2. Left and right Y-axes represent the RMSD values of protein and ligand in molecular distance unit Angstrom (Å).X-axis demonstrates the time in picoseconds (ps).

Figure 9 .
Figure 9.The timeline plot of hydrogen bonding interactions between compounds 1 and 2 and M pro of SARS-CoV2: (A) compound 1; (B) compound 2. Y-axis represents the number of hydrogen bonds, and X-axis demonstrates the time in picoseconds (ps).

Table 1 .
Calculated conformers of compound 1 that constitute more than 95% of the conformations' population.

Table 1 .
Calculated conformers of compound 1 that constitute more than 95% of the conformations' population.

Table 1 .
Calculated conformers of compound 1 that constitute more than 95% of the conformations' population.

Table 1 .
Calculated conformers of compound 1 that constitute more than 95% of the conformations' population.

Table 1 .
Calculated conformers of compound 1 that constitute more than 95% of the conformations' population.

Table 1 .
Calculated conformers of compound 1 that constitute more than 95% of the conformations' population.

Table 1 .
Calculated conformers of compound 1 that constitute more than 95% of the conformations' population.

Table 1 .
Calculated conformers of compound 1 that constitute more than 95% of the conformations' population.

Table 1 .
Calculated conformers of compound 1 that constitute more than 95% of the conformations' population.

Table 1 .
Calculated conformers of compound 1 that constitute more than 95% of the conformations' population.

Table 1 .
Calculated conformers of compound 1 that constitute more than 95% of the conformations' population.

Table 1 .
Calculated conformers of compound 1 that constitute more than 95% of the conformations' population.

Table 1 .
Calculated conformers of compound 1 that constitute more than 95% of the conformations' population.

Table 2 .
Calculated conformers of compound 2 that constitute more than 95% of the conformations' population.

Table 2 .
Calculated conformers of compound 2 that constitute more than 95% of the conformations' population.

Table 2 .
Calculated conformers of compound 2 that constitute more than 95% of the conformations' population.

Table 2 .
Calculated conformers of compound 2 that constitute more than 95% of the conformations' population.

Table 2 .
Calculated conformers of compound 2 that constitute more than 95% of the conformations' population.

Table 2 .
Calculated conformers of compound 2 that constitute more than 95% of the conformations' population.
double bonds in s-cis relationship compared with at least 20.4% of the conformers' population (conformers 3, 6, 8, and 10) having them in s-trans relationship.Together, these results indicate the preference of compound 2 to accommodate the half-chair conformation in which the hydroxyl group accommodates the pseudo-axial position and has s-cis relationships for the olefinic double bond/oxirane ring and the olefinic/carbonyl double bonds (Table2).

Table 2 .
Calculated conformers of compound 2 that constitute more than 95% of the conformations' population.

Table 2 .
Calculated conformers of compound 2 that constitute more than 95% of the conformations' population.

Table 2 .
Calculated conformers of compound 2 that constitute more than 95% of the conformations' population.