The Masamune-Bergman reaction[1],[2] is an example of a highly unusual class of chemical mechanism[3] involving the presumed formation of the biradical species shown as Int1 below by cyclisation of a cycloenediyne reactant. Such a species is so reactive that it will be quickly trapped, as for example by dihydrobenzene to form the final product. This cycloenediyne is not just an obscure chemical curiosity, the motif is incorporated into the natural product Calicheamicin, which is a potent antitumor antibiotic discovered in the 1980s. This drug owes its activity to the cyclisation TS1 shown below, which for n=2 occurs at the low temperature of 310K. The resulting biradical Int1 is a potent hydrogen abstractor, the species acting this way for hydrogen atoms associated with deoxyribose of DNA, ultimately leading to strand scission. Although I have explored many a mechanism on this blog using computational methods, I have never included any biradical examples. Here I explore the computational aspects of this reaction, and also include a pathway proceeding vis TS2- Int2 – TS3 in which hydrogen abstraction precedes cyclisation, in order to see how competitive such an alternative might be as a function of the ring size (n in scheme below).
The computational procedure was ωB97XD/Def2-TZVPP and the FAIR data is collected at DOI: 10.14469/hpc/14546 [4]. A spin unrestricted procedure is adopted using an approximation to allow for biradicaloid species, namely an initial first guess at the wavefunction using the keyword guess(mixed) which mixes what would be the HOMO and the LUMO of the molecule in a closed shell sense to allow a combination which includes an open shell singlet with one electron in the HOMO and one electron in the LUMO (a biradical). Part of the purpose of this approach is to try to find out if it gives reasonable results for such a mechanism. I will introduce the spin expectation operator <S2> to help identify biradicals. For closed shell singlets it has the value 0.0, for a pure biradical it has the value 1.0. Thus for species Int1, the values are typically ~0.995 and for the preceding TS1 ~ 0.3 to 0.57. IRC (Intrinsic reaction coordinate) calculations for TS1 show a smooth transition from values of <S2> = 0.0 (Reactant) through to 1.0 (Int1).
The results are shown below for three values of n, revealing that as the ring size increases (ending with an acyclic system Et2) the free energy barrier increases significantly, as indeed is reported[1],[2]. The alternative pathway proceeding via TS2 is always higher in free energy and varies much less with ring size. This route can therefore be firmly excluded from contention.
Table. Free energies for two mechanistic routes | ||||||
---|---|---|---|---|---|---|
System | Reactant | TS1 | TS2 | Int2 | TS3 | Int1 |
n=1 | -580.843968 0.0 | -580.806441 23.6 | -580.771042 45.8 | -580.797829 29.0 | -580.795875 30.2 | -580.838524 3.5 |
n=2 | -620.142895 0.0 | -620.090298 33.0† | -620.068740 46.5 | -620.094239 30.5 | -620.088955 33.8 | -620.131731 7.0‡ |
n=3 | -659.434065 0.0 | -659.370635 39.8 | -659.356146 48.9 | -659.384469 31.1 | -659.375211 36.9 | -659.413969 12.6 |
Et2 | -621.348992 0.0 | -621.278904 44.0 | -621.265041 52.7 | -621.292104 35.7 | -621.280607 42.9 | -621.319521 18.5 |
†<S2> =0.27. ‡Final Product (n=2) = -620.327868 (-116.1 kcal/mol)
This computational modelling largely agrees with the observations made for this reaction, with just one inconsistency. For n=2, the reaction is reported as taking place at 37°C, for which a typical free energy barrier would be in the region of ~24±2 kcal/mol,[5] around 9 kcal/mol lower than the computed value at this level of theory. This could originate from either a deficiency in the computational model, possibly in the handling of the open shell biradicaloid character by use of a simple spin unrestricted model,[6] or incursion of some lower energy process into the mechanism (free radical involvement?). I will continue probing this issue to see if its origins can be identified.
In the next part of this blog, I will investigate the mechanism as applied to Calicheamicin to see how the more complex bicycloenediyne nature of this natural product affects it.
References
- N. Darby, C.U. Kim, J.A. Salaün, K.W. Shelton, S. Takada, and S. Masamune, "Concerning the 1,5-didehydro[10]annulene system", J. Chem. Soc. D, vol. 0, pp. 1516-1517, 1971. http://dx.doi.org/10.1039/C29710001516
- R.R. Jones, and R.G. Bergman, "p-Benzyne. Generation as an intermediate in a thermal isomerization reaction and trapping evidence for the 1,4-benzenediyl structure", Journal of the American Chemical Society, vol. 94, pp. 660-661, 1972. http://dx.doi.org/10.1021/ja00757a071
- R.K. Mohamed, P.W. Peterson, and I.V. Alabugin, "Concerted Reactions That Produce Diradicals and Zwitterions: Electronic, Steric, Conformational, and Kinetic Control of Cycloaromatization Processes", Chemical Reviews, vol. 113, pp. 7089-7129, 2013. http://dx.doi.org/10.1021/cr4000682
- K.C. Nicolaou, G. Zuccarello, C. Riemer, V.A. Estevez, and W.M. Dai, "Design, synthesis, and study of simple monocyclic conjugated enediynes. The 10-membered ring enediyne moiety of the enediyne anticancer antibiotics", Journal of the American Chemical Society, vol. 114, pp. 7360-7371, 1992. http://dx.doi.org/10.1021/ja00045a005
- E.M. Greer, C.V. Cosgriff, and C. Doubleday, "Computational Evidence for Heavy-Atom Tunneling in the Bergman Cyclization of a 10-Membered-Ring Enediyne", Journal of the American Chemical Society, vol. 135, pp. 10194-10197, 2013. http://dx.doi.org/10.1021/ja402445a