
==== Front
J Org Chem
J Org Chem
jo
joceah
The Journal of Organic Chemistry
0022-3263
1520-6904
American Chemical Society

39248097
10.1021/acs.joc.4c00859
Article
Theoretical Study on the Biosynthesis of the Mandapamates: Mechanistic Insights Using Density Functional Theory
Wang Di †
Pattenden Gerald ‡
https://orcid.org/0000-0002-8165-9477
Fow Kam Loon †
https://orcid.org/0000-0003-3046-137X
Stocks Michael J. §
https://orcid.org/0000-0002-2726-0983
Hirst Jonathan D. *‡
https://orcid.org/0000-0002-4301-2496
Tang Bencan *†
† Department of Chemical and Environmental Engineering, Nottingham Ningbo China Beacons of Excellence Research and Innovation Institute, Key Laboratory for Carbonaceous Waste Processing and Process Intensification Research of Zhejiang Province, the University of Nottingham Ningbo China, 199 Taikang East Road, Ningbo 315100, P. R. China
‡ School of Chemistry, University of Nottingham, University Park, Nottingham NG7 2RD, U.K.
§ Nottingham University Biodiscovery Institute, School of Pharmacy, University of Nottingham, University Park, Nottingham NG7 2RD, U.K.
* Email: jonathan.hirst@nottingham.ac.uk.
* Email: bencan.tang@nottingham.edu.cn.
09 09 2024
20 09 2024
89 18 1294612956
10 04 2024
23 08 2024
28 07 2024
© 2024 The Authors. Published by American Chemical Society
2024
The Authors
https://creativecommons.org/licenses/by/4.0/ Permits the broadest form of re-use including for commercial purposes, provided that author attribution and integrity are maintained (https://creativecommons.org/licenses/by/4.0/).

Density functional theory (B3LYP-D3(BJ) and ωB97XD) calculations have been used to assess the stereochemical outcomes of the proposed transannular [4 + 2] cycloaddition pathway for the biosynthesis of mandapamate and isomandapamate from macrocyclic intermediates. Calculations reveal that the topological shift between macrocyclic conformers is vital in controlling the stereoselectivity of the downstream steps toward the isomeric mandapamates. A stepwise 4 + 2 type process is energetically favored over a concerted [4 + 2] pathway at room temperature, and is consistent with the stereochemistries found in the natural products.

Natural Science Foundation of Ningbo Municipality 10.13039/100007834 2022J171 Department of Science, Innovation and Technology NA NA Ningbo Municipal Key Laboratory on Clean Energy Conversion Technologies NA 2014A22010 National Key Research and Development Program of China 10.13039/501100012166 2018YFE0101700 Science and Technology Department of Zhejiang Province 10.13039/501100008990 NA Ningbo Municipal Bureau of Science and Technology 10.13039/501100007928 2020Z092 National Natural Science Foundation of China 10.13039/501100001809 22171153 National Natural Science Foundation of China 10.13039/501100001809 21502101 Royal Academy of Engineering 10.13039/501100000287 NA document-id-old-9jo4c00859
document-id-new-14jo4c00859
ccc-price
==== Body
pmcIntroduction

Mandapamate 1 and its stereoisomer isomandapamate 2 are novel polycyclic cembranoid diterpenes found in soft corals.1 Other polycyclic cembranoids found in soft corals include plumarellide 3, intricarene 4 and bielschowskysin 5, and the closely related C19 norditerpenes ineleganolide 6 and sinulochmodin C (7).2 Each of these polycyclic structures is purported to be derived in nature from macrocyclic furanocembranoid structures, e.g., 8(3) and 9,4 via sequences involving oxidation/hydration followed by transannular cyclizations (Figure 1).5 Thus, in the case of the mandapamates 1 and 2, they are thought to be biosynthesised from the hypothetical precursor 10 via 11/12, leading to the key exo-enol ether intermediate 15 which then undergoes transannular [4 + 2] cycloaddition (Scheme 1a). Stable exo-enol ether structures similar to 15 were first reported in corals by Fenical et al.6 Others were described later, e.g., 16 and 17,7 and exo enol ether structures are also thought to be implicated in biosynthetic pathways which lead to plumarellide 3 and to bielschowskysin 5.5

Figure 1 Representative polycyclic cembranoid diterpenes, C19 norcembranoids, and macrocyclic furanocembranes.

Scheme 1 (a) Biosynthetic Proposal for the Formation of Mandapamate 1 and Isomandapamate 2

(b and c) Reported Synthetic Attempts towards Mandapamate 1 and Isomandapamate 2 Ring System.

Mandapamate 1 and isomandapamate 2 share a common angular [5,6,7]-carbocyclic ring system with (3R, 6S)- and (3S, 6R)-stereochemistry respectively at their C3 and C6 centers (Figure 1). This stereogenicity originates from the relative orientation of the oxygen atoms on the bridged cycloheptenes in their structures. We postulated that during the proposed biosynthesis of the mandapamates there is likely to be restraints imposed on the conformations of their macrocyclic enol ether/furanoxonium ion intermediates, viz. 13–15 (Scheme 1a), which would predetermine the stereochemical outcomes of the subsequent transannular cyclizations, thereby leading to the different stereochemistries present in the two natural products. The origins of these proposed conformational preferences are not known, which has motivated us to establish some theoretical support for their existence, and thereby importance, using density functional theory (DFT).

The complexities of the ring systems and the oxidation patterns present in the naturally occurring polycycles 1–7, combined with the interesting biological properties they display, have made them challenging targets for the synthetic chemist. Pertinent to our study, there has been a recent novel synthesis of the core carbocyclic ring system in bielschowskysin 5 which was based on an intramolecular [2 + 2] cyclization involving an exo enol ether intermediate akin to that shown in compound 15.8,9 Other synthetic studies implicating enol ether intermediates have been carried out toward plumarellide 3,10 and total syntheses of both ineleganolide 6 and sinulochmodin C (7) have recently been accomplished.11 Earlier Pattenden et al.12 and Trauner et al.13 independently achieved a biomimetic synthesis of (+)-intricarene 4 from bipinnatin J (8) implicating an intramolecular dipolar [5 + 2] cycloaddition reaction as the key step, and in 2011, Li and Pattenden presented biomimetic syntheses of ineleganolide 6 and sinulochmodin C (7) from the macrocyclic 5-episinuleptolide 18 (believed to be derived from deoxypukalide 9 in vivo) via sequences of transannular Michael reactions.14 Some studies have been made toward the synthesis of the mandapamates 1 and 2 exploring an intramolecular [4 + 2] cyclization strategy to produce their [5,6,7] tricyclic ring systems, viz. 19 to 20 and 21 to 22 (Scheme 1b, c),15,16 but with limited success.

The challenges presented by the synthesis of the naturally occurring polycycles 1–7, implicating their likely biosynthesis precursors, have also prompted theoretical chemists to examine features of the mechanisms of these likely transformations. Quantum chemical calculations are now widely used to uncover new mechanisms17 as well as seek theoretical support for biosynthetic proposals.18 Related computational studies have recently been reviewed.18a Tang and Paton, for example, recently used DFT to probe the pathway to the cyclobutane rings in the furanocembranoid providencin19 and in bielschowskysin 5. Thus, based on quantum chemical calculations of the thermal and photochemical pathways toward the formation of the cyclobutane ring in bielschowskysin 5, they concluded that while a photo [2 + 2] cycloaddition pathway is computed to be highly efficient, a thermal ring-closure process in a single step in water is also energetically feasible at room temperature.20 More recently, alongside their synthetic work toward bielschowskysin 5, Scesa et al. also used DFT calculations to analyze the conformation of the likely furanocembranoid precursors, e.g., deoxypukalide 9 to its formation.8b Earlier, Lygo et al. used DFT calculations to study the steps leading to the polycyclic ring system found in plumarellide 3 and compared a stepwise cyclization pathway with a [4 + 2] cycloaddition of a furanoxonium ion intermediate.21 Hence, transannular cycloaddition reactions implicating furanoxonium ions as key intermediates have been verified both experimentally5b,5e,10,16 and computationally with DFT.22

In the proposed biosynthesis of the mandapamates 1 and 2 from the furanocembranoid precursor 10 (Scheme 1a) the key exo enol ether intermediate 15, can exhibit both E- and Z-geometry, each of which can undergo either a concerted or stepwise transannular 4 + 2 cycloaddition process leading to the isomeric mandapamates. Known furancembranes akin to the hypothetical precursor 10 have been isolated and/or synthesized, including bipinnatin J (8),3,12,13Z-/E-deoxypukalide 9,4 leptolide,23 pukalide aldehyde,23 molestin E,24 in addition to two model precursors designed by Mehta15 and Pattenden.16 The existence of these furancembranes gives us confidence that our proposed precursors could be derived from naturally existing furancembranes, if they themselves are not natural products.

To gain insight into this biosynthesis proposal we have carried out detailed quantum chemical calculations of the key steps and the likely stereochemical and conformational requirements of any macrocyclic furanoxonium ion intermediate implicated in competitive transannular 4 + 2 cycloaddition processes. The key point for this calculation is to understand the feasibility of the proposed transannular reactions leading to mandapamate and isomandapamate, despite the possibility that some side reactions might occur experimentally. For example, transannular cycloaddition between the furan ring and the exo C(15)=C is a possibility. However, our calculations revealed that the associated transition state (TS) has a very high activation energy barrier of over 59.0 kcal mol–1, which rules out this pathway (Figure S7). Similarly, selective dihydroxylation or epoxidation of the C7–C8 double bond with the existence of C15–C16 double bond could be problematic. Therefore, a good number of synthetic attempts/trails would be required to find the right reaction conditions to achieve these transformations.

Computational Methods

All the calculations reported were carried out in Gaussian 16.25 Unless stated otherwise, B3LYP-D3(BJ)26/def2-TZVP27//B3LYP/6-31G(d)28 and ωB97XD29/def2-TZVP//B3LYP/6-31G(d) relative free energies at 298.15 K, obtained from optimized species including an implicit SMD description of water,30 are quoted throughout, the latter in parentheses. The vibrational frequencies were computed at the same level of theory as for the geometry optimizations to confirm whether each optimized structure is an energy minimum (possessing zero imaginary frequencies) or a TS (possessing one single imaginary frequency), and to obtain the zero-point vibrational energy and thermal corrections under 298.15 K and 1 atm pressure. All transition states were confirmed to connect reactants and products by intrinsic reaction coordinate calculations.31 All Gibbs free energies discussed in the manuscript were corrected using the quasi-rigid rotor-harmonic oscillator approach proposed by Grimme32 and implemented in the GoodVibes code33 (Figures S5, S6 and Table S5). The 3D graphics of molecules were generated using CYLview.34 Electrostatic potential surface (EPS) analysis based on the noncovalent interactions was performed using Multiwfn (v3.8_dev)35 and visualized using VMD.36 The level of the theory that we have adopted was based on functional/basis set benchmark studies carried out in this work (Tables S1, S2 and S3). The B3LYP-D3(BJ) and ωB97XD functionals provided similar accuracy and qualitative conclusions. We selected B3LYP/6-31G(d) level for the geometry optimization in order to achieve the best balance between accuracy and computational cost. The 0.0 kcal mol–1 reference includes the single point energies of 11b, H2O, H3O+ and CH3OH. The structures of H2O, H3O+ and CH3OH were optimized and computed at the same level of 11b. Conformational analyses of 11 and 12 were carried out via random searching in the GMMX 3.1 module of Gaussian 16 using the MMFF9437 force field. The computational details are given in the Supporting Information.

Results and Discussion

Conformational Analysis of the Macrocyclic Furanocembranes 11/12 to 13

Referring to the biosynthesis proposal (Scheme 1a) we first analyzed the expected topological shift in the macrocyclic ring of the intermediates 11 to 13 resulting from conformational changes in the orientation of their furan rings. Using relaxed coordinate scans, the potential energy surfaces for these conformational changes were obtained at the B3LYP/6-31G(d) level through the incremental rotation of the O–C6–C7–C8 dihedral angles (10° increments each frame, see Supporting Information). As shown in Figure 2a, the energy minimum for conformer 11a is almost equi-energetic with conformer 11b (+1.1 kcal mol–1). These two conformers differ in the relative location of the furan-ring plane. The conformer 11a is parallel while the conformer 11b is orthogonal to the macrocycle. The activation Gibbs free energy toward the TS TS11a,11b is 3.4 kcal mol–1 with a reverse barrier of 2.3 kcal mol–1 (11b → TS11a,11b), so that the interconversion between conformers 11a and 11b occurs readily. Following a pivotal exothermic process (furan dearomatisation), the conformers 11a and 11b are converted into the more stable furanoxonium ion intermediates Int-13a and Int-13b respectively. Interestingly, the configurations (E- vs Z-) of the Δ6,7 double bonds in 13 are divergent at this stage and they will be retained in the downstream steps, thereby controlling the C3 and C6 stereocenters in the final mandapamates. The pathways for the interconversion of the epoxy compounds (cis-12a⇄cis-12b, trans-12a⇄trans-12b) presented similar transition states TScis-12a,cis-12b and TStrans-12a,trans-12b (Figure 2b). Compounds cis-12 and trans-12 differ in the orientations of their epoxy moieties. For the cis epoxide, the activation Gibbs free energy toward the TS TScis-12a,cis-12b is 6.6 kcal mol–1 with a reverse barrier of 3.3 kcal mol–1 (cis-12b → TScis-12a,cis-12b). In the corresponding trans epoxide, the activation Gibbs free energy toward the TS TStrans-12a,trans-12b is 4.6 kcal mol–1 with a reverse barrier of 4.0 kcal mol–1 (trans-12b → TStrans-12a,trans-12b). These two transition states (TScis-12a,cis-12b and TStrans-12a,trans-12b) also connect to the corresponding parallel conformers (cis-12a and trans-12a) and to the orthogonal conformers (cis-12b and trans-12b) which are in equilibrium and lead to the formation of the furanoxonium ion intermediates Int-13a and Int-13b with E- and Z-configured (Δ6,7)-double bonds, respectively. regioselectivity. To investigate further the origin of the conformational preference in the transition states TS11a,11b, TScis-12a,cis-12b and TStrans-12a,trans-12b, the EPS analysis was performed. As displayed in Figure 2c, d, in the transition states TS11a,11b and TScis-12a,cis-12b, the moiety containing a hydroxy or epoxy group at C7 shows a negative surface potential (blue, rO···O 2.70 Å), while the oxygen atom on furan ring also shows a negative surface potential (blue, rO···O 2.87 Å), indicating that there is electrostatic repulsion between them. However, in the TS TStrans-12a,trans-12b, the oxygen atom on the furan ring shows a neutral surface potential (white, rO···O 3.70 Å), indicating that there is less electrostatic interaction. This also explains why the activation energy barrier of TScis-12a,cis-12b is higher than that of TStrans-12a,trans-12b. Therefore, the EPS analysis suggests that the electrostatic interaction would be the driving force for the rotation about the furan ring.

Figure 2 (a, b) Calculated conformational energy profiles and transformations toward Int-13a and Int-13b respectively. Free energies calculated with SMD(H2O)-B3LYP-D3(BJ)/def2-TZVP//B3LYP/6-31G(d) and SMD(H2O)-ωB97XD/def2-TZVP//B3LYP/6-31G(d) (quoted in parentheses) respectively. All energies are relative to that of 11b. (c) Color scheme of key transition states: H, white; C, gray; O, red. Bond distances in Å. (d) EPS analysis of the transition states TS11a,11b, TScis-12a,cis-12b and TStrans-12a,trans-12b.

Pathway for the Formation of Mandapamate 1

As already discussed, the proposed biosynthesis of the mandapamates 1 and 2 (Scheme 1a) proceeds via the central stable furanoxonium ion intermediate Int-13b derived from 11b and 12b. The conformations of 11b and 12b are orthogonal and hence the intermediate Int-13b will have a Z-configured (Δ6,7)-double bond. This Z-configured (Δ6,7)-Int-13b is a direct optimization result of the pivotal furan dearomatisation process from 11b and 12b. Because extensive searches could not locate transition states connecting both sides of 11b/12b and int-13, we consider that this furan dearomatisation process is exothermic and is likely to be a barrierless process on the PES. The furanoxonium ion intermediate Int-13b is next quenched by methanol leading to the intermediate Int-14b (Figure 3). This process proceeds via TS1b with an activation Gibbs free energy of 9.2 kcal mol–1. The resulting protonated intermediate Int-23b then undergoes deprotonation to afford the stable exo enol ether-cyclic ketal 15b. Lastly, 15b undergoes a transannular [4 + 2] cycloaddition via TS2b leading to mandapamate 1 with an activation Gibbs free energy of 16.0 kcal mol–1. However, the overall energy barrier of this pathway (Int-13b → TS2b, 36.3 kcal mol–1) seems to be overly high for a biosynthesis route. We therefore investigated the alternative, stepwise, intramolecular cyclization pathway from 13b to mandapamate 1. Indeed, our calculations demonstrated that the Z-configured (Δ6,7) intermediate Int-13b could undergo transannular C–C bond formation via TS3b (Figure 3) with an activation Gibbs free energy of 6.6 kcal mol–1 to generate the thermodynamically stable cyclopentane-substituted intermediate Int-24b. A second transannular C–C bond formation through TS4b then leads to the intermediate Int-25b containing the tricyclic ring system found in mandapamate 1 with an activation Gibbs free energy of 25.5 kcal mol–1. Methanolysis of the intermediate Int-25b, to Int-26b, followed by deprotonation finally produces mandapamate 1. Attempts to locate the transition states connecting both sides of int-25 and int-26 (MeOH nucleophilic attack) were unsuccessful, indicating this is a barrierless process. Compared with the aforementioned concerted [4 + 2] cycloaddition pathway, the overall energy barrier in this stepwise route (Int-13b → TS4b, 24.9 kcal mol–1, without enzymatic intervention) is more reasonable. Our calculations revealed that B3LYP-D3(BJ) and ωB97XD functionals predicted the same rate-determining step, with M06-2X predicted a much lower Gibbs free activation energy for TS4b, see Table S3 in the Supporting Information, ΔG⧧(TS4b) is ca. 25 kcal mol–1 predicted by B3LYP-D3(BJ) and ωB97XD, but 21.4 kcal mol–1 predicted by M06-2X. One might anticipate that with enzymatic intervention, the energy barriers would typically decrease by at least a few kcal mol–1. The TS TS3b is favored over the TS1b by 12.6 kcal mol–1 and the TS TS4b is favored over the TS2b by 11.4 kcal mol–1 respectively.

Figure 3 Calculated energy profile for the formation of mandapamate 1. Geometry optimized with B3LYP/6-31G(d) in gas phase. Free energies calculated with SMD(H2O)-B3LYP-D3(BJ)/def2-TZVP//B3LYP/6-31G(d) and SMD(H2O)-ωB97XD/def2-TZVP//B3LYP/6-31G(d) (quoted in parentheses) respectively. Bond distances in Å.

Pathway for the Formation of Isomandapamate 2

As illustrated in Figure 4, isomandapamate 2 can also be obtained through a competing concerted [4 + 2] cycloaddition and stepwise 4 + 2 type pathway. Thus, following the acid-catalyzed conversions of the compounds 11a and 12a (Scheme 1a) the stable furanoxonium ion intermediate Int-13a, containing an E-configured (Δ6,7)-double bond, is produced. This intermediate can then either proceed via TS1a (ΔG⧧ of −4.7 kcal mol–1) and TS2a (ΔG⧧ of 13.4 kcal mol–1) or via TS3a (ΔG⧧ of −16.9 kcal mol–1) and TS4a (ΔG⧧ of 1.7 kcal mol–1) leading to the tricyclic ring system found in isomandapamate 2 (Figure 4). However, the overall energy barrier for the concerted [4 + 2] cycloaddition route (Int-13a → TS2a, 35.2 kcal mol–1) remains higher than the stepwise 4 + 2 type route (Int-13a → TS4a, 23.5 kcal mol–1), which further indicates that the stepwise transannular C–C bond formation sequence to isomandapamate 2, similar to that found in mandapamate 1, is more feasible at room temperature.

Figure 4 Calculated energy profile for the formation of isomandapamate 2. Geometry optimized with B3LYP/6-31G(d) in gas phase. Free energies calculated with SMD(H2O)-B3LYP-D3(BJ)/def2-TZVP//B3LYP/6-31G(d) and SMD(H2O)-ωB97XD/def2-TZVP//B3LYP/6-31G(d) (quoted in parentheses) respectively. Bond distances in Å.

Relative Stereochemistries of the Mandapamates 1 and 2

Our calculations reveal that in the proposed biosynthesis of the mandapamates, the orientations of the oxygen atoms in the furan ring of the furanoxonium ion intermediate Int-13 will be maintained once the latter is formed. Accordingly, the Z-configured (Δ6,7) intermediate Int-13b points toward (3R, 6S)-mandapamate 1 and, likewise, the E-configured (Δ6,7) intermediate Int-13a points toward (3S, 6R)-isomandapamate 2. Moreover, according to previous conformational analyses, the configuration of the furanoxonium ion intermediate Int-13 is controlled by the conformational changes shown in the dihydroxy and epoxy compounds 11 and 12 respectively. We can conclude therefore that the initial conformational changes in the macrocyclic furanocembrane precursors dictate the stereoselectivity in the formation of mandapamate 1 and isomandapamate 2.

Concerted [4 + 2] Cycloaddition vs Stepwise 4 + 2 Type Pathway

The distortion/interaction model originally conceived for bimolecular reactive systems, can also be applied to intramolecular reactions with a fragmentation scheme.38 As displayed in Figure 5a, we computed the strain-related energies of distorting various fragments in the intramolecular transition structures TS2, TS3 and TS4. The fragments (bold moieties) are obtained after being separated and filled valences with hydrogen atoms where the covalent bonds have been broken. The distortion/interaction analysis revealed that the distortion energies of TS3 (TS3a: ΔE⧧dist = 29.9 kcal mol–1, TS3b: ΔE⧧dist = 34.6 kcal mol–1) representing the first C7–C11 bond formation in the stepwise pathway are lower than that of TS2 representing concerted pathway (i.e., the concerted formation of the bonds C6–C14 and C7–C11). Compared with TS2, the bond distance of the second C6–C14 bond formation in TS4 are shorter by 0.85 and 0.63 Å respectively, indicating the stronger interaction in the transition states TS4 (TS4a: ΔE⧧int = −35.6 kcal mol–1, TS4b: ΔE⧧int = −24.4 kcal mol–1) of the stepwise pathway. Through constrained coordinate scans, we obtained the PES describing the mechanism of formation of the transannular C–C bonds C6–C14 and C7–C11 at the B3LYP/6-31G(d) level. The most striking feature of the [4 + 2] cycloaddition pathway toward mandapamate 1 is that there is only one single saddle point connecting the intermediate 15b and 1 (Figure 3), and clearly this corresponds to the concerted TS TS2b shown in Figure 5b (left). By contrast, the stepwise 4 + 2 type pathway to mandapamate 1 presents two saddle points connecting three species with much lower energy barriers in Figure 5b (right). Thus, the first transannular C–C bond formation occurs between C7 and C11 via TS TS3b (Figure 5b, right), and then the second transannular C–C bond forms immediately via TS TS4b, leading to formation of the key intermediate Int-25b (Figure 3) containing the tricyclic ring system in the natural product. Likewise, the above is also applicable to the pathway toward isomandapamate 2, due to the presence of the same carbocyclic skeleton, albeit with different orientations of the oxygen-bridged moiety.

Figure 5 (a) The distortion/interaction analysis of optimized transition states toward mandapamate 1 and isomandapamate 2 at the SMD(H2O)-B3LYP-D3(BJ)/def2-TZVP//B3LYP/6-31G(d) level. Energies in kcal mol–1. Bond distances in Å. Color scheme of key transition states: H, white; C, gray; O, red. (b) Contour plot of the PES of [4 + 2] concerted and stepwise pathways. Relative energies with respect to the lowest structure energy.

Conclusion

We have explored the concerted and the stepwise pathways leading to the formation of the tricyclic ring systems in the mandapamates 1 and 2 using quantum chemical calculations (Scheme 2). Our calculations reveal that the topological shift between conformers in the macrocyclic precursors to 1 and 2 is significant in controlling which one leads to the (3R, 6S)-mandapamate 1 and which to the (3S, 6R)- isomandapamate 2 via Z-configured (Δ6,7) Int-13b and E-configured (Δ6,7) Int-13a respectively. Our study therefore highlights the significance and importance of analyzing the conformations of macrocyclic molecules in order to appreciate their inherent ring-strained reactivity.39 The energy barriers between the conformers Int-13a and Int-13b (ΔΔG11a⇄11b = 1.1 kcal mol–1, ΔΔGcis-12a⇄cis-12b = 3.3 kcal mol–1, ΔΔGtrans-12a⇄trans-12b = 0.6 kcal mol–1) would allow facile conformational changes in conditions in vivo. Comparing the proposed cyclization routes to the ring systems in the mandapamates, calculations demonstrate that a stepwise sequence is energetically favored over an alternative concerted [4 + 2] cycloaddition pathway at room temperature. Our calculations also demonstrate that the stepwise pathway explains the formation of the different stereochemistries found in the naturally occurring isomeric mandapamates 1 and 2. Our study has provided some new mechanistic insight. However, the actual biosynthetic pathway may differ from that suggested by our calculations. For example, enzymes can alter reaction pathways beyond simply lowering a TS barrier. We can only rely on synthetic strategies in the lab to explore the chemistry of synthesizing macrocyclic precursors as well as biomimetically synthesizing mandapamates 1 and 2. Synthetic work toward the synthesis of macrocyclic furanocembranoid 11/12 for the biomimetic syntheses of mandapamates 1 and 2 is currently being carried out in our laboratory. Progress will be reported in due course.

Scheme 2 Summary of the Computed Biosyntheses of Mandapamate 1 and Isomandapamate 2

Data Availability Statement

The data underlying this study are available in the published article and its Supporting Information

Supporting Information Available

The Supporting Information is available free of charge at https://pubs.acs.org/doi/10.1021/acs.joc.4c00859.Details of DFT calculations (PDF)

Supplementary Material

jo4c00859_si_001.pdf

Author Contributions

The manuscript was written through contributions of all authors. All authors have given approval to the final version of the manuscript.

The authors declare no competing financial interest.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (Grant no. 22171153 & 21502101), Ningbo Science and Technology Bureau under CM2025 programme (Grant no. 2020Z092), Ningbo Natural Science Foundation Programme (Grant no. 2022J171), Ministry of Science and Technology of the People’s Republic of China under funding scheme National Key R&D Program of Intergovernmental Kay Projects (Grant no. 2018YFE0101700), and Ningbo Municipal Key Laboratory on Clean Energy Conversion Technologies (Grant no. 2014A22010) as well as the Zhejiang Provincial Key Laboratory for Carbonaceous Waste Processing and Process Intensification Research funded by the Zhejiang Provincial Department of Science and Technology (Grant no. 2020E10018). Professor Hirst acknowledges support from the Department of Science, Innovation and Technology (DSIT) and the Royal Academy of Engineering under the Chairs in Emerging Technologies scheme. The authors also acknowledge The University of Nottingham, England and the University of Nottingham Ningbo China for the use of their High-Performance Computing facilities.
==== Refs
References

For isolations, see:

a Venkateswarlu Y. ; Biabani M. A. F. ; Reddy M. V. R. ; Rao T. P. ; Kunwar A. C. ; Faulkner D. J. Mandapamate, a diterpenoid from the soft coral Sinularia dissecta. Tetrahedron Lett. 1994, 35 , 2249–2252.
b Anjaneyulu A. S. R. ; Sagar K. S. ; Venugopal M. J. R. V. Terpenoid and steroid constituents of the Indian ocean soft coral Sinularia maxima. Tetrahedron 1995, 51 , 10997–11010. 10.1016/0040-4020(95)00655-R.
For isolations, see:

a Stonik V. A. ; Kapustina I. I. ; Kalinovsky A. I. ; Dmitrenok P. S. ; Grebnev B. B. New diterpenoids from the far-eastern gorgonian coral Plumarella sp. Tetrahedron Lett. 2002, 43 , 315–317. 10.1016/S0040-4039(01)02114-1.
b Marrero J. ; Rodríguez A. D. ; Barnes C. L. Intricarene, an unprecedented trispiropentacyclic diterpene from the caribbean sea plume Pseudopterogorgia kallos. Org. Lett. 2005, 7 , 1877–1880. 10.1021/ol0505961.15844929
c Marrero J. ; Rodríguez A. D. ; Baran P. ; Raptis R. G. ; Sánchez J. A. ; Ortega-Barria E. ; Capson T. L. Bielschowskysin, a gorgonian-derived biologically active diterpene with an unprecedented carbon skeleton. Org. Lett. 2004, 6 , 1661–1664. 10.1021/ol049495d.15128261
d Duh C.-Y. ; Wang S.-K. ; Chia M.-C. ; Chiang M. Y. A novel cytotoxic norditerpenoid from the formosan soft coral Sinularia inelegans. Tetrahedron Lett. 1999, 40 , 6033–6035. 10.1016/S0040-4039(99)01194-6.
e Tseng Y.-J. ; Ahmed A. F. ; Dai C.-F. ; Chiang M. Y. ; Sheu J.-H. Sinulochmodins A–C, three novel terpenoids from the soft coral Sinularia lochmodes. Org. Lett. 2005, 7 , 3813–3816. 10.1021/ol051513j.16092882
Rodríguez A. D. ; Shi J.-G. The first cembrane–pseudopterane cycloisomerization. J. Org. Chem. 1998, 63 , 420–421. 10.1021/jo971884g.11672022
Dorta E. ; Díaz-Marrero A. R. ; Brito I. ; Cueto M. ; D’Croz L. ; Darias J. The oxidation profile at C-18 of furanocembranolides may provide a taxonomical marker for several genera of octocorals. Tetrahedron 2007, 63 , 9057–9062. 10.1016/j.tet.2007.06.080.
For reviews of the polycyclic cembranoid and norcembranoid diterpenes, see:

a Roethle P. A. ; Trauner D. The chemistry of marine furanocembranoids, pseudopteranes, gersolanes, and related natural products. Nat. Prod. Rep. 2008, 25 , 298–317. 10.1039/b705660p.18389139
b Li Y. ; Pattenden G. Perspectives on the structural and biosynthetic interrelationships between oxygenated furanocembranoids and their polycyclic congeners found in corals. Nat. Prod. Rep. 2011, 28 , 1269–1310. 10.1039/c1np00023c.21637894
c Li Y. ; Pattenden G. Novel macrocyclic and polycyclic norcembranoid diterpenes from Sinularia species of soft coral: structural relationships and biosynthetic speculations. Nat. Prod. Rep. 2011, 28 , 429–440. 10.1039/c0np00029a.21212904
d Craig R. A. ; Stoltz B. M. Polycyclic furanobutenolide-derived cembranoid and norcembranoid natural products: biosynthetic connections and synthetic efforts. Chem. Rev. 2017, 117 , 7878–7909. 10.1021/acs.chemrev.7b00083.28520418
e Palframan M. J. ; Pattenden G. Biosynthetic interrelationships within polycyclic cembranoids isolated from corals: conjecture, biomimetic synthesis and reality. Eur. J. Org Chem. 2020, 2020 , 2330–2349. 10.1002/ejoc.201901438.
Bandurraga M. M. Natural Product Studies of Selected East Pacific Gorgonians. Ph.D. Thesis, University of California, San Diego (USA), 1981.
a Venkateswarlu Y. ; Sridevi K. V. ; Rao M. R. New furanocembranoid diterpenes from the soft coral Sinularia maxima. J. Nat. Prod. 1999, 62 , 756–758. 10.1021/np9804076.10346962
b Epifanio R. d. A. ; Maia L. F. ; Fenical W. Chemical defenses of the endemic brazilian gorgonian Lophogorgia violacea pallas (octocorallia, gorgonacea). J. Braz. Chem. Soc. 2000, 11 , 584–591. 10.1590/S0103-50532000000600006.
c Sánchez M. C. ; Ortega M. J. ; Zubía E. ; Carballo J. L. Cembrane diterpenes from the gorgonian Lophogorgia peruana. J. Nat. Prod. 2006, 69 , 1749–1755. 10.1021/np060388x.17190454
d Kamel H. N. ; Ferreira D. ; Garcia-Fernandez L. F. ; Slattery M. Cytotoxic diterpenoids from the hybrid soft coral Sinularia maxima × Sinularia polydactyla. J. Nat. Prod. 2007, 70 , 1223–1227. 10.1021/np070074p.17676898
e Grote D. ; Dahse H.-M. ; Seifert K. Furanocembranoids from the soft corals Sinularia asterolobata and Litophyton arboretum. Chem. Biodiversity 2008, 5 , 2449–2456. 10.1002/cbdv.200890210.
f Thao N. P. ; Nam N. H. ; Cuong N. X. ; Quang T. H. ; Tung P. T. ; Tai B. H. ; Luyen B. T. T. ; Chae D. ; Kim S. ; Koh Y.-S. ; Kiem P. V. ; Minh C. V. ; Kim Y. H. Diterpenoids from the soft coral Sinularia maxima and their inhibitory effects on lipopolysaccharide-stimulated production of pro-inflammatory cytokines in bone marrow-derived dendritic cells. Chem. Pharm. Bull. 2012, 60 , 1581–1589. 10.1248/cpb.c12-00756.
g Chitturi B. R. ; Tatipamula V. B. ; Dokuburra C. B. ; Mangamuri U. K. ; Tuniki V. R. ; Kalivendi S. V. ; Bunce R. A. ; Yenamandra V. Pambanolides A–C from the south Indian soft coral Sinularia inelegans. Tetrahedron 2016, 72 , 1933–1940. 10.1016/j.tet.2016.02.056.
a Scesa P. ; Wangpaichitr M. ; Savaraj N. ; West L. ; Roche S. P. A kinetic dearomatization strategy for an expedient biomimetic route to the bielschowskysin skeleton. Angew. Chem., Int. Ed. 2018, 57 , 1316–1321. 10.1002/anie.201711780.
b Scesa P. D. ; West L. M. ; Roche S. P. Role of macrocyclic conformational steering in a kinetic route toward bielschowskysin. J. Am. Chem. Soc. 2021, 143 , 7566–7577. 10.1021/jacs.1c03336.33945689
For another construction of the carbocyclic ring system in 5 see

a Nicolaou K. C. ; Adsool V. A. ; Hale C. R. H. An expedient synthesis of a functionalized core structure of bielschowskysin. Angew. Chem., Int. Ed. 2011, 50 , 5149–5152. 10.1002/anie.201101360.
a Palframan M. J. ; Pattenden G. Elaboration of the carbocyclic ring systems in plumarellide and rameswaralide using a coordinated intramolecular cycloaddition approach, based on a common biosynthesis model. Tetrahedron Lett. 2013, 54 , 324–328. 10.1016/j.tetlet.2012.11.043.
b Palframan M. J. ; Pattenden G. The versatility of furfuryl alcohols and furanoxonium ions in synthesis. Chem. Commun. 2014, 50 , 7223–7242. 10.1039/c4cc01196a.
a Tuccinardi J. P. ; Wood J. L. Total syntheses of (+)-ineleganolide and (−)-sinulochmodin C. J. Am. Chem. Soc. 2022, 144 , 20539–20547. 10.1021/jacs.2c09826.36283051
b Gross B. M. ; Han S.-J. ; Virgil S. C. ; Stoltz B. M. A convergent total synthesis of (+)-ineleganolide. J. Am. Chem. Soc. 2023, 145 , 7763–7767. 10.1021/jacs.3c02142.36989438
c Zhang Y.-P. ; Du S. ; Ma Y. ; Zhan W. ; Chen W. ; Yang X. ; Zhang H. Structure-unit-based total synthesis of (−)-sinulochmodin C. Angew. Chem., Int. Ed. 2024, 63 , e202315481 10.1002/anie.202315481.
d For a recent synthesis of rameswaralide see Truax N. J. ; Ayinde S. ; Liu J. O. ; Romo D. Total synthesis of rameswaralide utilizing a pharmacophore-directed retrosynthetic strategy. J. Am. Chem. Soc. 2022, 144 , 18575–18585. 10.1021/jacs.2c08245.36166374
Tang B. ; Bray C. D. ; Pattenden G. A biomimetic total synthesis of (+)-intricarene. Tetrahedron Lett. 2006, 47 , 6401–6404. 10.1016/j.tetlet.2006.06.150.
Roethle P. A. ; Hernandez P. T. ; Trauner D. Exploring biosynthetic relationships among furanocembranoids: synthesis of (−)-bipinnatin J, (+)-intricarene, (+)-rubifolide, and (+)-isoepilophodione B. Org. Lett. 2006, 8 , 5901–5904. 10.1021/ol062581o.17134301
Li Y. ; Pattenden G. Biomimetic syntheses of ineleganolide and sinulochmodin C from 5-episinuleptolide via sequences of transannular Michael reactions. Tetrahedron 2011, 67 , 10045–10052. 10.1016/j.tet.2011.09.040.
Vasamsetty L. ; Khan F. A. ; Mehta G. A model approach towards the polycyclic framework present in cembranoid natural products dissectolide A, plumarellide and mandapamate. Tetrahedron Lett. 2014, 55 , 7068–7071. 10.1016/j.tetlet.2014.10.141.
a Li Y. ; Pattenden G. Exploration of a proposed biomimetic synthetic route to plumarellide. Development of a facile transannular Diels–Alder reaction from a macrocyclic enedione leading to a new 5,6,7-tricyclic ring system. Tetrahedron Lett. 2011, 52 , 2088–2092. 10.1016/j.tetlet.2010.10.154.
b Li Y. ; Palframan M. J. ; Pattenden G. ; Winne J. M. A strategy towards the synthesis of plumarellide based on biosynthesis speculation, featuring a transannular 4 + 2 type cyclisation from a cembranoid furanoxonium ion intermediate. Tetrahedron 2014, 70 , 7229–7240. 10.1016/j.tet.2014.06.090.
a Peng Q. ; Paton R. S. Catalytic control in cyclizations: from computational mechanistic understanding to selectivity prediction. Acc. Chem. Res. 2016, 49 , 1042–1051. 10.1021/acs.accounts.6b00084.27137131
b Cheng G.-J. ; Zhang X. ; Chung L. W. ; Xu L. ; Wu Y.-D. Computational organic chemistry: bridging theory and experiment in establishing the mechanisms of chemical reactions. J. Am. Chem. Soc. 2015, 137 , 1706–1725. 10.1021/ja5112749.25568962
c Nguyen Q. N. N. ; Tantillo D. J. The many roles of quantum chemical predictions in synthetic organic chemistry. Chem.—Asian J. 2014, 9 , 674–680. 10.1002/asia.201301452.24376231
d Peng Q. ; Duarte F. ; Paton R. S. Computing organic stereoselectivity – from concepts to quantitative calculations and predictions. Chem. Soc. Rev. 2016, 45 , 6093–6107. 10.1039/C6CS00573J.27722685
e Sato H. ; Saito K. ; Yamazaki M. Acceleration of mechanistic investigation of plant secondary metabolism based on computational chemistry. Front. Plant Sci. 2019, 10 , 802–817. 10.3389/fpls.2019.00802.31293608
f Houk K. N. ; Liu F. Holy grails for computational organic chemistry and biochemistry. Acc. Chem. Res. 2017, 50 , 539–543. 10.1021/acs.accounts.6b00532.28945400
a Wang D. ; Zhou T. ; Ren J. ; Hirst J. D. ; Gao Z. ; Tang B. Applications of quantum chemistry in biomimetic syntheses of polycyclic furanocembrane derivatives. Synlett 2024, 35 , 565–575. 10.1055/a-2239-6577.
b Hong Y. J. ; Tantillo D. J. A potential energy surface bifurcation in terpene biosynthesis. Nat. Chem. 2009, 1 , 384–389. 10.1038/nchem.287.21378892
c Hong Y. J. ; Tantillo D. J. Biosynthetic consequences of multiple sequential post-transition-state bifurcations. Nat. Chem. 2014, 6 , 104–111. 10.1038/nchem.1843.24451585
d Hornsby C. E. ; Paton R. S. It’s all downhill from here. Nat. Chem. 2014, 6 , 88–89. 10.1038/nchem.1852.24451579
e Tantillo D. J. Interrogating chemical mechanisms in natural products biosynthesis using quantum chemical calculations. Wiley Interdiscip. Rev.: Comput. Mol. Sci. 2020, 10 , e1453 10.1002/wcms.1453.
f Hayes C. J. ; Palframan M. J. ; Pattenden G. The impact of macrocycle conformation on the taxadiene-forming carbocation cascade: insight gained from sobralene, a recently discovered verticillene isomer. J. Org. Chem. 2020, 85 , 4507–4514. 10.1021/acs.joc.0c00369.32101002
Tang B. ; Paton R. S. Biosynthesis of providencin: understanding photochemical cyclobutane formation with density functional theory. Org. Lett. 2019, 21 , 1243–1247. 10.1021/acs.orglett.8b03838.30767534
Tang B. ; Simion R. ; Paton R. S. Thermal and photochemical mechanisms for cyclobutane formation in bielschowskysin biosynthesis. Synlett 2015, 26 , 501–507. 10.1055/s-0034-1379893.
Lygo B. ; Palframan M. J. ; Pattenden G. Investigation of transannular cycloaddition reactions involving furanoxonium ions using DFT calculations. Implications for the origin of plumarellide and rameswaralide and related polycyclic metabolites isolated from corals. Org. Biomol. Chem. 2014, 12 , 7270–7278. 10.1039/C4OB00734D.25105757
See also, computational studies of the biomimetic formation of intricarene 4

a Wang S. C. ; Tantillo D. J. Theoretical studies on synthetic and biosynthetic oxidopyrylium–alkene cycloadditions: pericyclic pathways to intricarene. J. Org. Chem. 2008, 73 , 1516–1523. 10.1021/jo7023762.18205383
Gutiérrez M. ; Capson T. L. ; Guzmán H. M. ; González J. ; Ortega-Barría E. ; Quiñoá E. ; Riguera R. Leptolide, a new furanocembranolide diterpene from Leptogorgiaalba. J. Nat. Prod. 2005, 68 , 614–616. 10.1021/np049745z.15844963
a Chu M.-J. ; Tang X.-L. ; Han X. ; Li T. ; Luo X.-C. ; Jiang M.-M. ; Van Ofwegen L. ; Luo L.-Z. ; Zhang G. ; Li P.-L. ; Li G.-Q. Metabolites from the Paracel Islands soft coral Sinularia cf. molesta. Mar. Drugs 2018, 16 , 517–531. 10.3390/md16120517.30572615
b Hoff O. ; Kratena N. ; Aynetdinova D. ; Christensen K. E. ; Donohoe T. J. A vicinal diol approach for the total synthesis of molestin E, ent-sinulacembranolide A and ent-sinumaximol A. Chem.—Eur. J. 2022, 28 , e202202464 10.1002/chem.202202464.35946550
Frisch M. J. ; Trucks G. W. ; Schlegel H. B. ; Scuseria G. E. ; Robb M. A. ; Cheeseman J. R. ; Scalmani G. ; Barone V. ; Petersson G. A. ; Nakatsuji H. ; Li X. ; Caricato M. ; Marenich A. V. ; Bloino J. ; Janesko B. G. ; Gomperts R. ; Mennucci B. ; Hratchian H. P. ; Ortiz J. V. ; Izmaylov A. F. ; Sonnenberg J. L. ; Williams-Young D. ; Ding F. ; Lipparini F. ; Egidi F. ; Goings J. ; Peng B. ; Petrone A. ; Henderson T. ; Ranasinghe D. ; Zakrzewski V. G. ; Gao J. ; Rega N. ; Zheng G. ; Liang W. ; Hada M. ; Ehara M. ; Toyota K. ; Fukuda R. ; Hasegawa J. ; Ishida M. ; Nakajima T. ; Honda Y. ; Kitao O. ; Nakai H. ; Vreven T. ; Throssell K. ; Montgomery J. A. Jr. ; Peralta J. E. ; Ogliaro F. ; Bearpark M. J. ; Heyd J. J. ; Brothers E. N. ; Kudin K. N. ; Staroverov V. N. ; Keith T. A. ; Kobayashi R. ; Normand J. ; Raghavachari K. ; Rendell A. P. ; Burant J. C. ; Iyengar S. S. ; Tomasi J. ; Cossi M. ; Millam J. M. ; Klene M. ; Adamo C. ; Cammi R. ; Ochterski J. W. ; Martin R. L. ; Morokuma K. ; Farkas O. ; Foresman J. B. ; Fox D. J. Gaussian 16, Revision A.03; Gaussian, Inc.: Wallingford CT, 2016
a Becke A. D. Density-functional thermochemistry. III. The role of exact exchange. J. Chem. Phys. 1993, 98 , 5648–5652. 10.1063/1.464913.
b Lee C. ; Yang W. ; Parr R. G. Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density. Phys. Rev. B 1988, 37 , 785–789. 10.1103/PhysRevB.37.785.
c Grimme S. ; Antony J. ; Ehrlich S. ; Krieg H. A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-D) for the 94 elements H-Pu. J. Chem. Phys. 2010, 132 , 154104–154122. 10.1063/1.3382344.20423165
d Grimme S. ; Ehrlich S. ; Goerigk L. Effect of the damping function in dispersion corrected density functional theory. J. Comput. Chem. 2011, 32 , 1456–1465. 10.1002/jcc.21759.21370243
Weigend F. ; Ahlrichs R. Balanced basis sets of split valence, triple zeta valence and quadruple zeta valence quality for H to Rn: Design and assessment of accuracy. Phys. Chem. Chem. Phys. 2005, 7 , 3297–3305. 10.1039/b508541a.16240044
a Ditchfield R. ; Hehre W. J. ; Pople J. A. Self-consistent molecular-orbital methods. IX. An extended Gaussian-type basis for molecular-orbital studies of organic molecules. J. Chem. Phys. 1971, 54 , 724–728. 10.1063/1.1674902.
b Hehre W. J. ; Ditchfield R. ; Pople J. A. Self-consistent molecular orbital methods. XII. Further extensions of Gaussian-type basis sets for use in molecular orbital studies of organic molecules. J. Chem. Phys. 1972, 56 , 2257–2261. 10.1063/1.1677527.
c Hariharan P. C. ; Pople J. A. The influence of polarization functions on molecular orbital hydrogenation energies. Theoret. Chim. Acta 1973, 28 , 213–222. 10.1007/BF00533485.
Chai J.-D. ; Head-Gordon M. Long-range corrected hybrid density functionals with damped atom-atom dispersion corrections. Phys. Chem. Chem. Phys. 2008, 10 , 6615–6620. 10.1039/b810189b.18989472
a Tomasi J. ; Mennucci B. ; Cammi R. Quantum mechanical continuum solvation models. Chem. Rev. 2005, 105 , 2999–3094. 10.1021/cr9904009.16092826
b Marenich A. V. ; Cramer C. J. ; Truhlar D. G. Universal solvation model based on solute electron density and on a continuum model of the solvent defined by the bulk dielectric constant and atomic surface tensions. J. Phys. Chem. B 2009, 113 , 6378–6396. 10.1021/jp810292n.19366259
a Fukui K. The path of chemical reactions - the IRC approach. Acc. Chem. Res. 1981, 14 , 363–368. 10.1021/ar00072a001.
b Hratchian H. P. ; Schlegel H. B. Theory and Applications of Computational Chemistry: The First 40 Years. In Dykstra C. E. , Frenking G. , Kim K. S. , Scuseria G. , Eds.; Elsevier: Amsterdam, 2005; pp 195–249.
Grimme S. Supramolecular binding thermodynamics by dispersion-corrected density functional theory. Chem.—Eur. J. 2012, 18 , 9955–9964. 10.1002/chem.201200497.22782805
Luchini G. ; Alegre-Requena J. V. ; Funes-Ardoiz I. ; Paton R. S. GoodVibes: Automated thermochemistry for heterogeneous computational chemistry data. F1000Research 2020, 9 , 291 10.12688/f1000research.22758.1.
Legault C. Y. CYLview20; Université de Sherbrooke, 2020. http://www.cylview.org (accessed on September 7, 2022).
Lu T. ; Chen F. Multiwfn A multifunctional wavefunction analyzer. J. Comput. Chem. 2012, 33 , 580–592. 10.1002/jcc.22885.22162017
Humphrey W. ; Dalke A. ; Schulten K. VMD - Visual molecular dynamics. J. Mol. Graphics 1996, 14 , 33–38. 10.1016/0263-7855(96)00018-5.
Halgren T. A. Merck molecular force field. I. Basis, form, scope, parameterization, and performance of MMFF94. J. Comput. Chem. 1996, 17 , 490–519. 10.1002/(sici)1096-987x(199604)17:5/6<490::aid-jcc1>3.0.co;2-p.
a Bickelhaupt F. M. ; Houk K. N. Analyzing reaction rates with the distortion/interaction-activation strain model. Angew. Chem., Int. Ed. 2017, 56 , 10070–10086. 10.1002/anie.201701486.
b Krenske E. H. ; Houk K. N. ; Holmes A. B. ; Thompson J. Entropy versus tether strain effects on rates of intramolecular 1,3-dipolar cycloadditions of N-alkenylnitrones. Tetrahedron Lett. 2011, 52 , 2181–2184. 10.1016/j.tetlet.2010.11.121.
Chen R. ; Shen Y. ; Yang S. ; Zhang Y. Conformational design principles in total synthesis. Angew. Chem., Int. Ed. 2020, 59 , 14198–14210. 10.1002/anie.202003735.
