- 1Department of Chemistry, University of Louisiana at Lafayette, Lafayette, LA, United States
- 2Department of Chemical Engineering, University of Louisiana at Lafayette, Lafayette, LA, United States
- 3School of Chemistry, University of Bristol, Bristol, United Kingdom
We present a contemporary mechanistic description of the light-driven conversion of cyclopropenone containing enediyne (CPE) precusors to ring-opened species amenable to further Bergman cyclization and formation of stable biradical species that have been proposed for use in light-induced cancer treatment. The transformation is rationalized in terms of (purely singlet state) Norrish type-I chemistry, wherein photoinduced opening of one C–C bond in the cyclopropenone ring facilitates non-adiabatic coupling to high levels of the ground state, subsequent loss of CO and Bergman cyclization of the enediyne intermediate to the cytotoxic target biradical species. Limited investigations of substituent effects on the ensuing photochemistry serve to vindicate the experimental choices of Popik and coworkers (J. Org. Chem., 2005, 70, 1297–1305). Specifically, replacing the phenyl moiety in the chosen model CPE by a 1,4-benzoquinone unit leads to a stronger, red-shifted parent absorption, and increases the exoergicity of the parent → biradical conversion.
Introduction
Drug discovery is one of the fastest growing and most potent fields in modern day science (Elion, 1993; Takenaka, 2001; Paul et al., 2010; Rask-Andersen et al., 2011; Swinney and Anthony, 2011; Warren, 2011; Lee et al., 2012). The use and development of state-of-the-art experimental methods, coupled with molecular modeling and artificial intelligence, are revolutionizing contemporary drug discovery (Chan et al., 2019; Henstock, 2019; Mak and Pichika, 2019). These methods have also been vital for developing thorough understandings of the physiological activities of a diverse range of drugs.
Photoinitiated drug activation, wherein UV/visible light is used to drive the action of a drug, already finds use in the context of psoralens and the treatment of a range of skin conditions (Stern et al., 1980; Menter et al., 2009). The more widespread use of light-sensitive drugs, however, is still in its infancy (Farrer et al., 2010; Vernooij et al., 2018; Shi et al., 2019, 2020; Imberti et al., 2020). Such photoinitiation therapies are based on the premise that light-irradiation of a chromophore leads to formation of an excited state, with a total energy well above that of the ground state. In most cases, the total energy of the excited state species exceeds many of the activation barriers associated with ground state reactions. Electronic excitation thus provides a means of activating a reaction that is otherwise unfavorable in the ground state configuration. The excited states formed by electronic excitation can decay in many ways, both radiative (e.g., fluorescence) and non-radiative (examples of which include internal conversion, reaction (e.g., isomerization) and dissociation). Such non-radiative decay processes are usually controlled (either in part or exclusively) by conical intersections (CIs) between potential energy surfaces, which arise when two or more electronic states—e.g., an excited state and the ground state—become degenerate along a given reaction path (Domcke et al., 2004, 2011; Credo Chung et al., 2007; Xie et al., 2018). The region of degeneracy is then viewed as a CI when orthogonal nuclear motions are considered. Within the chemical physics community, CIs are now widely recognized as mediators of ultrafast internal conversion between electronic states. However, the concept, involvement and importance of CIs has yet to be fully appreciated within the wider scientific community when discussing and utilizing molecular photochemistry. This manuscript is intended, in part, to help bridge the gap between synthetic photochemistry and mechanistic photophysics—by highlighting the crucial roles of CIs in a system wherein, hitherto, such processes have been generally unrecognized and thus ignored.
Photodissociations are a class of non-radiative excited state decay processes available to molecules with labile leaving groups. The initial photoexcitation in such cases may directly populate a so-called 1nσ* or 1πσ* state, i.e., a state formed by promoting an electron from a non-bonding (n) or bonding π valence orbital to an antibonding σ* orbital (Cronin et al., 2004; Ashfold et al., 2006, 2010; Devine et al., 2006; Nix et al., 2007; Roberts and Stavros, 2014; Marchetti et al., 2016). In many other cases, the primary photoexcitation may involve a strongly absorbing π*←π transition; light-driven bond fission in such cases only occurs after non-adiabatic coupling between the ππ* state and an (n/π)σ* continuum. In either scenario, the σ* character lowers the overall bond order of the bond around which the orbital is localized and encourages dissociation (or predissociation) along the relevant bond-stretch coordinate.
For completeness, we note that light irradiation already finds many other roles in medicine—most notably in photodynamic therapy (PDT). PDT relies on the sensitized formation of the highly reactive excited singlet (1Δg) state of molecular oxygen (henceforth 1O2), which can then oxidize biomolecules in the environment of cancerous cells and thereby restrict their biochemistry and eventual cell division. The sensitized formation of 1O2 requires the introduction and photoexcitation of a strongly absorbing chromophore with a high propensity for intersystem crossing (ISC). Resonant energy transfer between the photo-produced triplet state chromophore molecules and O2 elevates the latter to its reactive singlet state. Though effective, PDT has recognized shortcomings: e.g., (i) the eventual formation of 1O2 is indirect (relying on the excitation of, and energy transfer from, the chromophore species) and (ii) cancerous cells are usually hypoxic and successful PDT thus relies on the presence of dissolved O2 in the chromophore solution at the point that it is administered—thus limiting the control possible when performing PDT.
Light-activated chemotherapy drugs offer a direct means of targeted cancer therapy. Cyclopropenone-containing enediyne precursors are one such class of photoinitiated drugs that have been identified as potential candidates for photoinduced cancer therapy (Poloukhtine and Popik, 2005a,b; Poloukhtine et al., 2008; Pandithavidana et al., 2009). Upon photoexcitation, these molecules photodissociate, eliminating CO and forming an enediyne which then undergoes a Bergman cyclization to form a cycloalkene biradical product (Jones and Bergman, 1972; Luxon et al., 2018) that is capable of oxidizing DNA. The overall mechanism for one such precursor−2,3-benzobicyclo[8.1.0]undec-1(10)-en-4-yn-11-one—henceforth variously abbreviated as CPE and as A—is illustrated in Scheme 1.
Scheme 1. Molecular structures associated with the photoinduced dissociation of the model cyclopropenone-containing enediyne (A) and the subsequent Bergman rearrangement of (B) to the biradical (C).
The present manuscript reports the use of contemporary multi-reference electronic structure methods to explore mechanistic details of the reaction paths that lead to the observed (Poloukhtine and Popik, 2005a) photoinduced CO elimination from CPE and the ensuing Bergman rearrangement of B to C and how these processes might be affected by selected changes to the phenyl ring.
Computational Methodology
All calculations were conducted for the isolated (i.e., gas phase) molecule. The ground state minimum energy geometry of CPE (A) was optimized using the Coulomb-Attenuated Model Becke-3rd parameter-Lee-Yang-Parr (CAM-B3LYP) functional (Yanai et al., 2004) of Density Functional Theory (DFT), coupled to the 6-31G(d) Pople basis set (Hehre et al., 1970). Since the reaction of interest (Scheme 1) starts with photoinduced elimination of CO, two-dimensional potential energy (PE) profiles of the ground and first excited electronic states of this CPE along the two C–CO stretch coordinates of the cyclopropenone moiety were first investigated using Time-Dependent Density Functional Theory (TD-DFT) (Van Caillie and Amos, 2000; Furche and Ahlrichs, 2002; Scalmani et al., 2006), with the above functional and basis set. As shown in Supplementary Figure 1 in the Electronic Supplementary Information (ESI), these test calculations returned a maximum on the excited state PE surface when both C–CO bond distances were extended in tandem, but minima when either C–CO bond was extended while the other was held fixed. This hinted that the photoinduced CO elimination involves an initial ring-opening prior to release of the CO moiety.
Guided by these computations, the lowest triplet spin configuration of the ring-opened species formed by breaking the C–C bond nearer the phenyl ring was optimized at the CAM-B3LYP/6-31G(d) level of theory and used as a proxy for the singlet ring-opened structure and energy (since calculations for the ground singlet state inevitably resulted in reformation of the ring-closed CPE molecule).
PE profiles connecting the ground state CPE structure to the optimized ring-opened biradical structure and then onwards to the bicyclic structure B were then computed using the complete active space second-order perturbation theory (CASPT2) method, (Roos et al., 1982; Andersson et al., 1990, 1992; Park and Shiozaki, 2017) based on a state-averaged complete active space self-consistent field (SA-CASSCF) reference wavefunction (comprising five singlet and four triplet states) and coupled to the cc-pVDZ basis set (Dunning, 1989). The PE path was constructed using successive linear interpolations in internal coordinates (LIICs) and each LIIC was checked to ensure that the intermediate geometries represented a sensible interpolation between the optimized ground state CPE, the ring-opened species and then B. All of these CASSCF and CASPT2 computations employed an active space involving 10 electrons distributed in 10 orbitals—the five highest occupied (1–5) and the five lowest virtual (6–10) orbitals shown in Supplementary Figure 2 of the ESI, where 5 and 6 are, respectively, the highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) in the ground state. A further more limited set of PE profiles mapping the Bergman rearrangement of B to the biradical C were then computed using CASPT2/cc-pVDZ methods, based on a SA-CASSCF reference wavefunction comprising just two singlet and two triplet states (since our primary interest was in determining the magnitude of any energy barrier to B → C conversion on the ground state PE surface) and an active space comprising 10 electrons in 10 orbitals (shown in Supplementary Figure 3 of the ESI where, again, orbitals 5 and 6 are, respectively the HOMO and LUMO of the ground state molecule). In all cases, an imaginary level shift of 0.3 Eh was used to aid convergence and mitigate against the involvement of intruder states.
Vertical excitation energies (VEEs) and transition dipole moments were extracted from the associated CASSCF/CASPT2 calculations at the ground state equilibrium geometries of A (and B) to derive oscillator strengths (f) for transitions to the first few excited states (just the respective S1 states in the case of B) using Equation (1).
where Ei and E0 are, respectively, the energies of the excited state of interest and the ground state and the μ0i are the associated transition dipole moments along α (α = x, y, and z). The different x and y coordinates used in describing the transitions in A and B are shown by the red and green arrows in Supplementary Figures 2, 3; the z coordinate in both cases is out of the page.
Additional CASPT2 computations using the same basis set and imaginary level shift together with a (10,8) active space were undertaken for three analogs of CPE wherein the phenyl ring had been substituted symmetrically with methoxy (OMe) or cyano (CN) groups in the 1- and 4-positions, or replaced with a 1,4-benzoquinone ring, in order to explore substituent effects on the VEEs and transition strengths. Supplementary Figure 4 in the ESI shows the highest occupied and lowest virtual orbitals in the ground state of each system. These variants will henceforth be referred to as 1,4-MeO-CPE, 1,4-CN-CPE and 1,4-O=CPE. The latter was selected as (2,3-dimethyl substituted) 1,4-O=CPE has been synthesized (from the dimethylated analog of 1,4-MeO-CPE), shown to form the corresponding enediyne by photoinduced CO elimination and touted as a potential cytotoxin (Poloukhtine and Popik, 2005b). The other two substituted CPEs were selected to compare and contrast the effect of strongly π-accepting (i.e., CN) and π-donating (i.e., MeO) substitutents. The ground state energies of the A, B, and C analogs of 1,4-MeO-CPE, 1,4-CN-CPE and 1,4-O=CPE, relative to bare CPE, were optimized using the high-level CBS-QB3 method (Saracino et al., 2003; Cysewski, 2007; Rayne and Forest, 2010). All DFT and CBS-QB3 calculations were performed using Gaussian 16 (Frisch et al., 2016), while all CASSCF/CASPT2 computations were undertaken in Molpro 2018 (Werner et al., 2018).
Results and Discussion
Ground State Structure and Vertical Excitation Energies
Figures 1A–C display plan and side-on views of the ground state optimized geometry of CPE. The phenyl, cyclopropenone and alkyne moieties lie in a common plane, while the C4H8 alkyl link adopts a non-planar configuration. Given the cyclic nature of CPE, the geometry around the cyclopropenone unit can be viewed as a strained form of bare cyclopropenone—with an intra-large-ring C–C=C bond angle of ~175° (cf. an H–C=C bond angle of ~125° in cyclopropenone)–and the alkyne moiety adopts a strained form of butyne.
Figure 1. Plan (A) and side views (B,C) of the optimized molecular structure of CPE in its ground electronic state.
Table 1 lists the vertical excitation energies (VEEs) to the first few excited singlet and triplet states of CPE and the dominant orbital promotions involved in the first two singlet excitations. The calculated VEE of the more intense S2-S0 transition agrees well with that reported for CPE in Poloukhtine and Popik (2005a) and Poloukhtine and Popik (2006a), but the calculated S1-S0 transition energy is below that associated with the observed weaker, longer wavelength absorption maximum.
Table 1. Vertical excitation energies (VEEs) and oscillator strengths (f) of transitions to the first few singlet and triplet excited states of CPE, calculated at the CASPT2/AVDZ level of theory.
Reference to Table 1 and the orbitals displayed in Supplementary Figure 2 shows that the first excited singlet (S1) state of CPE has ππ* character and is formed by electron promotion between orbitals that include π-electron density distributed over the phenyl ring, cyclopropenone and alkyne moieties. The S2 state has both nπ* and ππ* character and is formed by promotion from a ‘mixed’ bonding orbital that involves a non-bonding pσ component (localized on the O atom) and a π-bonding component that is mainly localized on the phenyl ring but also conjugates across the C=C bond of the cyclopropenone unit. The participating n/π and π* orbitals show significant spatial overlap, reflected in a calculated S2-S0 oscillator strength some 2.5× larger than that of the S1-S0 transition (Table 1). This suggests a significant contribution from the phenyl-ring centered π-bonding region and, for brevity, we will henceforth refer to S2 as a ππ* excited state also.
For completeness, Table 1 also lists the calculated VEEs to the first three triplet states. Excitations to these states from the singlet ground state are spin-forbidden and thus calculated to have zero oscillator strength. We also note the near degeneracy of the T1 and T2 states in the vertical region, and that the T2 and S1 states share a common dominant electronic configuration at this geometry.
Reaction Path Following Photoexcitation
Figure 2 displays the PE profiles along the reaction path depicted in Scheme 1. Panel (a) shows the profiles returned from a LIIC along the dimensionless coordinate Qa connecting the minimum energy geometries of ground state CPE and the (triplet) intermediate formed via ring-opening of the cyclopropenone unit. As noted previously (Poloukhtine and Popik, 2006b), the alternative ring-opened intermediate that would be formed by cleaving the other C–C bond in the cyclopropenone unit (i.e., the bond further from the phenyl ring) lies at significantly higher energy (~0.35 eV in the present calculations) and is not considered further. As such, Qa is predominantly the C–C bond stretch indicated by the black arrow in Figure 2A. Along this path, the PE of the S0 state shows a steady increase—as would be expected for the rupture of a C–C σ-bond. The PE of the S1 state, in contrast, shows a small barrier and then decreases along this ring-opening coordinate. At this point, it is important to recall that LIIC pathways must, by construction, return upper estimates of any barriers that would be derived by locating and optimizing the true transition states, and we henceforth assume minimal barrier to distortion along Qa on the S1 PE surface. The topography of the S1 potential can be understood by recognizing that the electronic configuration of the S1 state (predominantly ππ* in the vertical region) gains increasing nσ* character upon C–C bond extension. The decline in the S1 state PE, together with the progressive increase in the S0 state PE, leads to an inevitable crossing at Qa = 0.9. This region of conical intersection is likely to encourage non-adiabatic coupling between the S1 and S0 states, i.e., to facilitate internal conversion to high vibrational levels of the S0 state of CPE. These topological details serve to reinforce the earlier conclusion that the light-induced decarbonylation in cyclopropenones is a stepwise process initiated on the S1 PE surface (Poloukhtine and Popik, 2006b).
Figure 2. CASPT2 potential energy profiles for the ground and first few excited states of CPE along the sequence of dimensionless LIICs (panels (a) through (c) as discussed in the text) involved in the photoinduced conversion of A → B (+ CO) products, with the molecular geometries at various critical points depicted. The color coding of the states shown in the inset in panel (A) defines the state ordering at Qa = 0. Note, the energetic ordering of the states labeled T1 and T2 switches once Qa > 0, and thereafter the adiabatic PE profiles of the T1 and T2 states are described by, respectively, the yellow and pale blue curves—here and in later figures.
This photoinduced ring opening can be viewed as another example of a singlet-mediated Norrish type I α-C–C bond cleavage in a ketone. Norrish type I reactions have traditionally been assumed to start with efficient S1 → T1 ISC (Norrish and Bamford, 1936, 1937; Diau et al., 2001, 2002), and indeed, reference to Figure 2A suggests that S1 → T1 ISC in CPE could constitute another route to the ring-opened species. However, the recent literature includes a growing number of predictions and/or demonstrations of purely singlet-state enabled fission of the C–C bond in the α-position (relative to the carbonyl group) (Maeda et al., 2010; Nádasdi et al., 2010; Marchetti et al., 2019), notably in strained cyclic ketones like cyclobutanone (Xia et al., 2015; Shemesh et al., 2016; Kao et al., 2020). In all such cases investigated to date, the C–C=O moiety remaining after photochemical cleavage of the other α-C–C bond is predicted to be near linear at the CI (Marchetti et al., 2019)—as is found in the present case also (see Figure 2). The photoinduced C–C bond fission in CPE is another such example, not least because—as we now show—Scheme 1 requires further chemistry that can only occur on the S0 PE surface. As Supplementary Figure 5 in the ESI shows, the Mulliken charge distributions for the optimized ring-opened structure of the isolated molecule in its T1 state and for the S0 state molecule at that same geometry are similar, and far from limiting biradical or zwitterionic in character. Clearly, however, the degree of charge separation is likely to be different in solution and to be solvent-dependent.
The highly vibrationally-excited S0 molecules formed by non-adiabatic coupling from the S1 state may relax to reform the starting CPE molecule or follow a rival path on the S0 PE surface. Given the experimental observations (Poloukhtine and Popik, 2005a, 2006b), one possible rival path must involve CO loss from the ring-opened biradical to form intermediate B in Scheme 1. To model this process, PE profiles were first computed simply by stepping the C–CO bond stretch coordinate (RC−CO), but this revealed an obvious discontinuity in the S0 and T1 PE profiles at RC−CO ~ 1.8 Å. This was traced to a localized linear to bent change in the geometry of the C–C=O unit, that we assume to be a signature of the onset of a π-stacking interaction between the emerging CO product and the alkyne group. Thus, the overall C–CO bond fission following opening of the cyclopropenone ring was modeled using two successive LIICs—Qb, to cover the region from the minimum energy geometry of the triplet ring-opened structure to the relaxed bent geometry at RC−CO ~1.8 Å, then Qc, from the bent geometry at RC−CO ~1.8 Å to long range.
The S0 and excited state PE profiles associated with Qb in Figure 2B show several features of note. The PE profile of the T1 state (which is the most stable configuration of the ring-opened species) increases steadily. Unsurprisingly, given that the S1 and T1 states at these geometries are well-approximated as spin-flip variations of the same electronic configuration, the S1 PE profile shows a similar increase. The PE of the S0 state, in contrast, declines steadily along Qb. This finding can be understood by inspecting the HOMO of the ring-opened species, which shows antibonding π* character around the C–CO bond. Occupancy of this orbital reduces the C–CO bond strength and, as with the initial ring-opening step, this HOMO gains progressive σ* character upon C–CO bond extension—leading to eventual bond fission. Figure 2C shows the corresponding PE profiles along Qc. Of particular relevance to the current narrative, the PE profile of the S0 state shows a continued steady decrease en route to forming the alkyne intermediate (structure B in Scheme 1).
Scheme 1 shows this alkyne intermediate undergoing a Bergman cyclization to form a ring-closed bicyclic biradical (structure C). Figure 3 shows the calculated (by CASPT2) energy profile for this final step on the S0 PE surface, along a LIIC (Qd) linking structures B and C. This shows the di-alkyne intermediate B and the tricyclic biradical species C having essentially the same minimum energies, and separated by an intervening energy barrier with ΔE ~0.8 eV. Reference to Figure 2 shows that the minimum energy of A is also very similar to that of the B (+CO) limit. Noting the earlier comment that any energy barrier returned by a LIIC calculation will necessarily be an upper estimate of the true transition state energy and that the calculated barrier energy in Figure 3 is well below that of the CI between the S1 and S0 PE surfaces of CPE (Figure 2A) suggests that both B and C could be thermodynamically viable products following photoexcitation of CPE, within which population may partition—consistent with the experimental observation of B formation [and the ensuing Bergman cyclisation to C at quite modest temperatures (84°C)] (Poloukhtine and Popik, 2005a).
Figure 3. CASPT2 potential energy profiles for the ground and first three excited states along the dimensionless LIIC Qd describing the B → C conversion.
Traditional organic photosyntheses such as that used to drive the transformation A → C employ prolonged illumination, so it is also prudent to consider possible photoinduced chemistry of the product B (Poloukhtine and Popik, 2005a,b, 2006a,b; Poloukhtine et al., 2008). Table 2 lists the calculated VEEs to the S1, T1, and T2 states of B (all of which are best viewed as having ππ* character) and the S1-S0 oscillator strength and dominant contributing promotions based on the orbitals displayed in Supplementary Figure 3. Figure 3 shows the calculated PE profiles of these three excited states along Qd. As Table 2 shows, the S1-S0 transition of B is predicted to have a similar oscillator strength to that of CPE but to have an appreciably larger VEE—which matches reasonably with the absorption maximum of the product formed by 300 nm irradiation of CPE reported in Poloukhtine and Popik (2005a).
Table 2. Vertical excitation energies (VEEs) and oscillator strengths (f) of transitions to the lowest singlet and triplet excited states of B, calculated at the CASPT2/AVDZ level of theory.
The PE profiles shown in Figure 3 hint that secondary photoexcitation of B could feasibly aid C formation. Specifically, we note the near degeneracy of the S1 and T2 states at small Qd and the possibility of strong non-adiabatic coupling around Qd ~ 0.4, where the S1 potential samples a region of configuration space that supports a conical intersection between the T2 and (not shown) T3 potentials. Any S1 population undergoing (spin-orbit induced) transfer to the T2 potential in this region could thereafter follow an energetically “down-hill” path via the T2/T1 conical intersection at Qd ~ 0.55 and relax into the potential well of the low-lying triplet form of C. Higher level theory could inform on the possible importance of spin-orbit enabled transfer to the T1 potential, but electron paramagnetic resonance (EPR) studies as a function of illumination time might be the best way of exploring the relative importances of one vs. sequential two photon production of C by near UV irradiation under synthetically relevant solution phase conditions.
Substituent Effects on the Electronic Absorption and Potential Reactivity of CPEs
Table 3 compares the VEEs, the x, y, and z components of the respective transition dipole moments and the overall oscillator strengths for forming the S1 and S2 states of bare CPE and the 1,4-MeO-CPE, 1,4-O=CPE and 1,4-CN-CPE analogs returned by the CASPT2 calculations. The orbital promotions that make the greatest contributions to the S1-S0 and S2-S0 transitions in each of the substituted molecules are included in Supplementary Figure 4 in the ESI.
Table 3. Vertical excitation energies (VEEs), x, y, and z components of the respective transition dipole moments [TDMs (in atomic units)] and oscillator strengths (f) for the S1-S0 and S2-S0 transitions of (a) bare CPE, (b) 1,4-methoxy substituted CPE, (c) the 1,4-benzoquinone derivative, and (d) 1,4-cyano substituted CPE, calculated at the CASPT2/AVDZ level of theory.
Even a cursory inspection of Table 3 suffices to show that such symmetric changes at the 1- and 4-positions on the phenyl ring have a substantial effect on the electronic absorption of CPE. All are predicted to enhance the S1-S0 oscillator strength (cf. bare CPE), most strikingly in the cases of 1,4-MeO-CPE and the 1,4-benzoquinone derivative. Each involves net loss of electron density from the cyclopropenone carbonyl group (reflected in the negative TDMy values in Table 3); the S1 molecules have some zwitterionic character—in accord with conclusions reached in earlier theoretical studies of alkyl- or aryl-substituted cyclopropenones (Poloukhtine and Popik, 2006b). Substituting the phenyl ring with π-accepting (i.e., CN) or π-donating (i.e., MeO) groups increases the VEE of the S1-S0 transition (cf. bare CPE), but the enhanced conjugation introduced by replacing the phenyl ring by a 1,4-benzoquinone ring causes an obvious red-shift in the S1-S0 absorption.
The effects on the S2-S0 transitions are more substituent specific, with the calculated TDMs suggesting accumulation of electron density on the substituted phenyl ring in the case of 1,4-MeO-CPE and on the cyclopropenone sub-unit in the case of 1,4-O=CPE. In all cases, substitution is predicted to reduce the VEE of the S2-S0 transition (cf. bare CPE)—in accord with the available experimental data for 1,4-O=CPE and CPE (Poloukhtine and Popik, 2005a,b)
Table 4 shows the effect of substituents on the relative stabilities of the bicyclic enediyne and the Bergman cyclized product (henceforth labeled B′ and C′) relative to the corresponding B and C products from photoexcitation of bare CPE. CO elimination (i.e., A′ → B′) is favored, energetically, by either CN substitution or by replacing the phenyl ring with a 1,4-benzoquinone unit, whereas 1,4-MeO substitution leads to a B′ product that is marginally less stable than for bare CPE. The former stabilizations can be plausibly understood by considering the extended π-systems in the enediyne product. Thermodynamically, 1,4-O=CPE is in a class of its own with regard to the final B′ → C′ electrocyclization reaction, consistent with its adoption as a precursor of choice by (Poloukhtine and Popik, 2006a).
Table 4. CBS-QB3 calculated enthalpy changes associated with the respective A′ → B′ and B′ → C′ conversions in (b) 1,4-methoxy substituted CPE, (c) the 1,4-benzoquinone derivative, and (d) 1,4-cyano substituted CPE, each referenced to the corresponding A → B and B → C conversions in bare CPE (a).
Conclusions
This study offers a contemporary mechanistic description of the light-driven conversion of cyclopropenone containing enediyne precusors like CPE to ring-opened species amenable to further Bergman cyclization, resulting in stable biradical species that have been proposed for future use in light-induced cancer treatment (Poloukhtine and Popik, 2005a,b, Poloukhtine and Popik, 2006a; Poloukhtine et al., 2008). The initial photoinduced cleavage of the C–C bond in the cyclopropenone ring closer to the phenyl group is shown to follow Norrish type-I chemistry. In the present case, this can be accommodated by purely singlet-state chemistry, i.e., efficient non-adiabatic coupling from the photo-prepared S1 state to high levels of the ground (S0) PE surface, upon which subsequent nuclear rearrangements lead to CO loss and formation of B, followed by the electrocyclization of B to yield the cytotoxic target biradical species C.
Viewed from the perspective of the photoexcited A* precursor in Scheme 1, the reactions through to ground state B (and C) products are exoergic processes, the efficiencies of which are likely to be sensitive to the relative rates of reaction and internal energy loss to the solvent. Deliberate (or otherwise) secondary photoexcitation of the enediyne intermediate B is suggested as a way of boosting the yield of the target diradical species C. The present limited investigation of the effects of chemically modifying the phenyl moiety in A serves to vindicate the experimental choices made by Popik and coworkers (Poloukhtine and Popik, 2006a). Specifically, replacing the phenyl moiety in CPE by a 1,4-benzoquinone unit is shown to both red-shift and boost the absorption of the A′ precursor and to increase the exoergicity of the required A′ → C′ conversion—complementing current efforts to use visible-light-responsive photocatalysts to trigger alkyne generation from (non-visible light absorbing) cyclopropenones (Mishiro et al., 2019).
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Author Contributions
TK and MA designed the project. SL, BM and TK undertook the underlying data. All authors contributed to drafting the manuscript.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The handling editor declared a past collaboration with the authors TK and MA.
Acknowledgments
TK thanks the University of Louisiana at Lafayette for start-up funds.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fchem.2020.596590/full#supplementary-material
References
Andersson, K., Malmqvist, P., and Roos, B. O. (1992). Second-order perturbation theory with a complete active space self-consistent field reference function. J. Chem. Phys. 96, 1218–1226. doi: 10.1063/1.462209
Andersson, K., Malmqvist, P. A., Roos, B. O., and Sadlej, A. (1990). Second-order perturbation theory with a CASSCF reference function. J. Phys. Chem. 94, 5483–5488. doi: 10.1021/j100377a012
Ashfold, M. N. R., Cronin, B., Devine, A. L., and Dixon, R. (2006). The role of πσ* excited states in the photodissociation of heteroaromatic molecules. Science 312, 1637–1640. doi: 10.1126/science.1125436
Ashfold, M. N. R., King, G. A., Murdock, D., and Nix, M. G. (2010). πσ* excited states in molecular photochemistry. Phys. Chem. Chem. Phys. 12, 1218–1238. doi: 10.1039/B921706A
Chan, H. C. S., Shan, H., Dahoun, T., Vogel, H., and Yuan, S. (2019). Advancing drug discovery via artificial intelligence. Trends Pharmacol. Sci. 40, 592–604. doi: 10.1016/j.tips.2019.06.004
Credo Chung, W., Lan, Z., Ohtsuki, Y., Shimakura, N., Domcke, W., and Fujimura, Y. (2007). Conical intersections involving the dissociative 1πσ* state in 9H-adenine: a quantum chemical ab initio study. Phys. Chem. Chem. Phys. 9, 2075–2084. doi: 10.1039/B618745E
Cronin, B., Nix, M. G. D., Qadiri, R. H., and Ashfold, M. N. (2004). High resolution photofragment translational spectroscopy studies of the near ultraviolet photolysis of pyrrole. Phys. Chem. Chem. Phys. 6, 5031–5041. doi: 10.1039/b411589a
Cysewski, P. (2007). Accurate gas phase acidities of carboxylic acids estimated by scaling the vibrational contribution of ab initio gibbs free energies. J. Mol. Model. 13, 801–808. doi: 10.1007/s00894-007-0200-0
Devine, A. L., Cronin, B., and Nix, M. G. D. (2006). High resolution photofragment translational spectroscopy studies of the near ultraviolet photolysis of imidazole. J. Chem. Phys. 125:184302. doi: 10.1063/1.2364504
Diau, E. W. G., Kötting, C., Sølling, T. I., and Zewail, A. H. (2002). Femtochemistry of norrish type-I reactions: III. highly excited ketones—theoretical. ChemPhysChem 3, 57–78. doi: 10.1002/1439-7641(20020118)3:1<57::AID-CPHC57>3.0.CO;2-F
Diau, E. W. G., Kötting, C., and Zewail, A. H. (2001). Femtochemistry of norrish type-I reactions: experimental, I. and theoretical studies of acetone and related ketones on the S1 surface. ChemPhysChem 2, 273–293. doi: 10.1002/1439-7641(20010518)2:5<273::AID-CPHC273>3.0.CO;2-H
Domcke, W., Yarkony, D. R., and Koeppel, H. (2004). Conical Intersections: Electronic Structure, Dynamics & Spectroscopy. doi: 10.1142/5406
Domcke, W., Yarkony, D. R., and Koeppel, H. (2011). Conical Intersections: Theory, Computation and Experiment. Singapore: World Scientific. doi: 10.1142/7803
Dunning, T. H. (1989). Gaussian basis sets for use in correlated molecular calculations. I. The atoms boron through neon and hydrogen. J. Chem. Phys. 90, 1007–1023. doi: 10.1063/1.456153
Elion, G. B. (1993). The quest for a cure. Annu. Rev. Pharmacol. Toxicol. 33, 1–25. doi: 10.1146/annurev.pa.33.040193.000245
Farrer, N. J., Woods, J. A., Salassa, L., and Zhao, Y. (2010). A potent trans-diimine platinum anticancer complex photoactivated by visible light. Angew. Chemie Int. Ed. 49, 8905–8908. doi: 10.1002/anie.201003399
Frisch, M. J., Trucks, G. W., and Schlegel, H. B. (2016). Gaussian 16, Revision A.03. Wallingford, CT: Gaussian Inc.
Furche, F., and Ahlrichs, R. (2002). Adiabatic time-dependent density functional methods for excited state properties. J. Chem. Phys. 117, 7433–7447. doi: 10.1063/1.1508368
Hehre, W. J., Ditchfield, R., Stewart, R. F., and Pople, J. (1970). Self-consistent molecular orbital methods. IV. use of gaussian expansions of slater-type orbitals. Extension to second-row molecules. J. Chem. Phys. 52, 2769–2773. doi: 10.1063/1.1673374
Henstock, P. V. (2019). Artificial intelligence for pharma: time for internal investment. Trends Pharmacol. Sci. 40, 543–546. doi: 10.1016/j.tips.2019.05.003
Imberti, C., Zhang, P., Huang, H., and Sadler, P. J. (2020). New designs for phototherapeutic transition metal complexes. Angew. Chemie Int. Ed. 59, 61–73. doi: 10.1002/anie.201905171
Jones, R. R., and Bergman, R. G. (1972). p-Benzyne. Generation as an intermediate in a thermal isomerization reaction and trapping evidence for the 1,4-benzenediyl structure. J. Am. Chem. Soc. 94, 660–661. doi: 10.1021/ja00757a071
Kao, M.-H., Venkatraman, R. K., and Ashfold, M. N. R. (2020). Effects of ring-strain on the ultrafast photochemistry of cyclic ketones. Chem. Sci. 11, 1991–2000. doi: 10.1039/C9SC05208A
Lee, J. A., Uhlik, M. T., and Moxham, C. M. (2012). Modern phenotypic drug discovery is a viable, neoclassic pharma strategy. J. Med. Chem. 55, 4527–4538. doi: 10.1021/jm201649s
Luxon, A. R., Orms, N., Kanters, R., and Krylov, A. I. (2018). An ab initio exploration of the bergman cyclization. J. Phys. Chem. A 122, 420–430. doi: 10.1021/acs.jpca.7b10576s
Maeda, S., Ohno, K., and Morokuma, K. (2010). A theoretical study on the photodissociation of acetone: insight into the slow intersystem crossing and exploration of nonadiabatic pathways to the ground state. J. Phys. Chem. Lett. 1, 1841–1845. doi: 10.1021/jz100551y
Mak, K.-K., and Pichika, M. R. (2019). Artificial intelligence in drug development: present status and future prospects. Drug Discov. Today 24, 773–780. doi: 10.1016/j.drudis.2018.11.014
Marchetti, B., Karsili, T. N. V., and Ashfold, M. N. R. (2016). A ‘bottom up’, ab initio computational approach to understanding fundamental photophysical processes in nitrogen containing heterocycles, DNA bases and base pairs. Phys. Chem. Chem. Phys. 18, 20007–20027. doi: 10.1039/C6CP00165C
Marchetti, B. V., Karsili, T. N., and Ashfold, M. N. R. (2019). Exploring norrish type I and type II reactions: an ab initio mechanistic study highlighting singlet-state mediated chemistry. Phys. Chem. Chem. Phys. 21, 14418–14428. doi: 10.1039/C8CP07292B
Menter, A., Korman, N. J., Elmets, C. A., and Feldman, S. (2009). Guidelines of care for the management of psoriasis and psoriatic arthritis: section 3. Guidelines of care for the management and treatment of psoriasis with topical therapies. J. Am. Acad. Dermatol. 60, 643–659. doi: 10.1016/j.jaad.2008.12.032
Mishiro, K., Kimura, T., Furuyama, T., and Kunishima, M. (2019). Phototriggered active alkyne generation from cyclopropenones with visible light-responsive photocatalysts. Org. Lett. 21, 4101–4105. doi: 10.1021/acs.orglett.9b01280
Nádasdi, R., Zügner, G. L., Farkas, M., Dóbé, S., and Maeda, S. (2010). Photochemistry of methyl ethyl ketone: quantum yields and S1/S0-diradical mechanism of photodissociation. ChemPhysChem 11, 3883–3895. doi: 10.1002/cphc.201000522
Nix, M. G. D., Devine, A. L., Cronin, B., and Ashfold, M. N. (2007). Ultraviolet photolysis of adenine: dissociation via the 1πσ* state. J. Chem. Phys. 126:124312. doi: 10.1063/1.2712843
Norrish, R. G. W., and Bamford, C. H. (1936). Photodecomposition of aldehydes and ketones. Nature 138:1016a0. doi: 10.1038/1381016a0
Norrish, R. G. W., and Bamford, C. H. (1937). Photo-decomposition of aldehydes and ketones. Nature 140, 195–196. doi: 10.1038/140195b0
Pandithavidana, D. R., Poloukhtine, A., and Popik, V. V. (2009). Photochemical generation and reversible cycloaromatization of a nine-membered ring cyclic enediyne. J. Am. Chem. Soc. 131, 351–356. doi: 10.1021/ja8077076
Park, J. W., and Shiozaki, T. (2017). Analytical derivative coupling for multistate CASPT2 theory. J. Chem. Theory Comput. 13, 2561–2570. doi: 10.1021/acs.jctc.7b00018
Paul, S. M., Mytelka, D. S., and Dunwiddie, C. T. (2010). How to improve R&D productivity: the pharmaceutical industry's grand challenge. Nat. Rev. Drug Discov. 9, 203–214. doi: 10.1038/nrd3078
Poloukhtine, A., Karpov, G., and Popik, V. V. (2008). Towards photoswitchable enediyne antibiotics. Curr. Top. Med. Chem. 8, 460–469. doi: 10.2174/156802608783955700
Poloukhtine, A., and Popik, V. V. (2005a). Application of photochemical decarbonylation of cyclopropenones for the in situ generation of reactive enediynes. Construction of a cyclopropenone-containing enediyne precursor by using a cyclopropenone acetal building block. J. Org. Chem. 70, 1297–1305. doi: 10.1021/jo048065y
Poloukhtine, A., and Popik, V. V. (2006a). Mechanism of the cyclopropenone decarbonylation reaction. A density functional theory and transient spectroscopy study. J. Phys. Chem. A 110, 1749–1757. doi: 10.1021/jp0563641
Poloukhtine, A., and Popik, V. V. (2006b). Two-photon photochemical generation of reactive enediyne. J. Org. Chem. 71, 7417–7421. doi: 10.1021/jo061285m
Poloukhtine, A., and Popik, V. V. (2005b). Photoswitchable enediynes: use of cyclopropenone as photocleavable masking group for the enediyne triple bond. Chem. Commun. 617–619. doi: 10.1039/b414951c
Rask-Andersen, M., Almén, M. S., and Schiöth, H. B. (2011). Trends in the exploitation of novel drug targets. Nat. Rev. Drug Discov. 10, 579–590. doi: 10.1038/nrd3478
Rayne, S., and Forest, K. (2010). Theoretical studies on the pKa values of perfluoroalkyl carboxylic acids. J. Mol. Struct. Theochem. 949, 60–69. doi: 10.1016/j.theochem.2010.03.003
Roberts, G. M., and Stavros, V. G. (2014). The role of πσ* states in the photochemistry of heteroaromatic biomolecules and their subunits: insights from gas-phase femtosecond spectroscopy. Chem. Sci. 5, 1698–1722. doi: 10.1039/c3sc53175a
Roos, B. O., Linse, P., and Siegbahn, P. E. M. (1982). A simple method for the evaluation of the second-order-perturbation energy from external double-excitations with a CASSCF reference wavefunction. Chem. Phys. 66, 197–207. doi: 10.1016/0301-0104(82)88019-1
Saracino, G. A. A., Improta, R., and Barone, V. (2003). Absolute pKa determination for carboxylic acids using density functional theory and the polarizable continuum model. Chem. Phys. Lett. 373, 411–415. doi: 10.1016/S0009-2614(03)00607-9
Scalmani, G., Frisch, M. J., Mennucci, B., Tomasi, J., and Cammi, R. (2006). Geometries and properties of excited states in the gas phase and in solution: theory and application of a time-dependent density functional theory polarizable continuum model. J. Chem. Phys. 124:94107. doi: 10.1063/1.2173258
Shemesh, D., Nizkorodov, S. A., and Gerber, R. B. (2016). Photochemical reactions of cyclohexanone: mechanisms and dynamics. J. Phys. Chem. A 120, 7112–7120. doi: 10.1021/acs.jpca.6b06184
Shi, H., Imberti, C., Huang, H., Hands-Portman, I., and Sadler, P. J. (2020). Biotinylated photoactive Pt(IV) anticancer complexes. Chem. Commun. 56, 2320–2323. doi: 10.1039/C9CC07845B
Shi, H., Imberti, C., and Sadler, P. J. (2019). Diazido platinum(IV) complexes for photoactivated anticancer chemotherapy. Inorg. Chem. Front. 6, 1623–1638. doi: 10.1039/C9QI00288J
Stern, R. S., Kleinerman, R. A., and Parrish, J. A. (1980). Phototoxic reactions to photoactive drugs in patients treated with PUVA. Arch. Dermatol. 116, 1269–1271. doi: 10.1001/archderm.116.11.1269
Swinney, D. C., and Anthony, J. (2011). How were new medicines discovered? Nat. Rev. Drug Discov. 10, 507–519. doi: 10.1038/nrd3480
Takenaka, T. (2001). Classical vs reverse pharmacology in drug discovery. BJU Int. 88, 7–10. doi: 10.1111/j.1464-410X.2001.00112.x
Van Caillie, C., and Amos, R. D. (2000). Geometric derivatives of density functional theory excitation energies using gradient-corrected functionals. Chem. Phys. Lett. 317, 159–164. doi: 10.1016/S0009-2614(99)01346-9
Vernooij, R. R., Joshi, T., Horbury, M. D., and Graham, B. (2018). Spectroscopic studies on photoinduced reactions of the anticancer prodrug, trans,trans,trans-[Pt(N3)2(OH)2(py)2]. Chem. A Eur. J. 24, 5790–5803. doi: 10.1002/chem.201705349
Warren, J. (2011). Drug discovery: lessons from evolution. Br. J. Clin. Pharmacol. 71, 497–503. doi: 10.1111/j.1365-2125.2010.03854.x
Werner, H.-J., Knowles, P. J., Knizia, G., Manby, F., R Schütz, M., Celani, P., et al. (2018). MOLPRO, Version 1, a Package of ab initio Programs. Available online at: https://www.molpro.net
Xia, S.-H., Liu, X.-Y., Fang, Q., and Cui, G. (2015). Excited-state ring-opening mechanism of cyclic ketones: a MS-CASPT2//CASSCF study. J. Phys. Chem. A 119, 3569–3576. doi: 10.1021/acs.jpca.5b00302
Xie, C., Malbon, C. L., Yarkony, D. R., and Xie, D. (2018). Signatures of a conical intersection in adiabatic dissociation on the ground electronic state. J. Am. Chem. Soc. 140, 1986–1989. doi: 10.1021/jacs.7b11489
Keywords: photodissociation, photostability, photophysics, physical organic chemistry, photochemistry, photoactive drugs
Citation: Léger SJ, Marchetti B, Ashfold MNR and Karsili TNV (2020) The Role of Norrish Type-I Chemistry in Photoactive Drugs: An ab initio Study of a Cyclopropenone-Enediyne Drug Precursor. Front. Chem. 8:596590. doi: 10.3389/fchem.2020.596590
Received: 19 August 2020; Accepted: 25 November 2020;
Published: 22 December 2020.
Edited by:
Michael Staniforth, University of Warwick, United KingdomReviewed by:
Samer Gozem, Georgia State University, United StatesArkaprabha Konar, Kent State University, United States
Copyright © 2020 Léger, Marchetti, Ashfold and Karsili. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Tolga N. V. Karsili, dG9sZ2Eua2Fyc2lsaSYjeDAwMDQwO2xvdWlzaWFuYS5lZHU=