Albert Marten
Schulte‡
a,
Georgios
Alachouzos‡
*a,
Wiktor
Szymanski
*ab and
Ben L.
Feringa
*a
aCentre for Systems Chemistry, Stratingh Institute for Chemistry, Faculty for Science and Engineering, University of Groningen, Nijenborgh 4, 9747 AG Groningen, The Netherlands. E-mail: alachouzos.georgios@gmail.com; w.c.szymanski@rug.nl; b.l.feringa@rug.nl
bDepartment of Radiology, Medical Imaging Center, University Medical Center Groningen, University of Groningen, Hanzeplein 1, 9713 GZ Groningen, The Netherlands
First published on 5th January 2024
Photocleavable protecting groups (PPGs) enable the precise spatiotemporal control over the release of a payload of interest, in particular a bioactive substance, through light irradiation. A crucial parameter that determines the practical applicability of PPGs is the efficiency of payload release, largely governed by the quantum yield of photolysis (QY). Understanding which parameters determine the QY will prove crucial for engineering improved PPGs and their effective future applications, especially in the emerging field of photopharmacology. The Contact Ion Pair (CIP) has been recognized as an important intermediate in the uncaging process, but the key influence of its fate on the quantum yield has not been explored yet, limiting our ability to design improved PPGs. Here, we demonstrate that the CIP escape mechanism of PPGs is crucial for determining their payload- and solvent-dependent photolysis QY, and illustrate that an intramolecular type of CIP escape is superior over diffusion-dependent CIP escape. Furthermore, we report a strong correlation of the photolysis QY of a range of coumarin PPGs with the DFT-calculated height of all three energy barriers involved in the photolysis reaction, despite the vastly different mechanisms of CIP escape that these PPGs exhibit. Using the insights obtained through our analysis, we were able to predict the photolysis QY of a newly designed PPG with particularly high accuracy. The level of understanding of the factors determining the QY of PPGs presented here will move the ever-expanding field of PPG applications forward and provides a blueprint for the development of PPGs with QYs that are independent of payload-topology and solvent polarity.
One cornerstone class of such photochemical tools are photolabile protecting groups (PPGs)6–13 that enable the use of light to control one of the simplest chemical processes: the breaking of a covalent sigma bond. Through enabling this primary photochemical process, PPGs have found applications as protecting groups of functional moieties in organic synthesis of small molecule building blocks and in total syntheses of natural products, and have shown to be highly versatile for the release of (bioactive) payloads in molecular biology, materials science, and in photopharmacology.14–23
One of the most widely used classes of PPGs is that of the heterolytic PPGs, containing among others coumarin-24–28 and BODIPY-PPGs.29–32 Their shared mechanism of payload release relies on excited state heterolysis of the bond between the payload and the α-carbon of the PPG (Fig. 1a, k1). In this process, a contact ion pair (CIP) intermediate is formed, in which the positive charge is located on the α-carbon of the PPG, and the negative charge on the payload. The initial heterolysis step is reversible, and – depending on its stability – the CIP will either recombine to reform the substrate (Fig. 1a, k−1), or react further through diffusion of the CIP and trapping by the solvent (such as water, Fig. 1a, k2) to release the payload and the PPG-alcohol.33–36
Since the applications of PPGs have exponentially expanded in the last decade,37–40 much of their development has focused on tuning their efficiency of payload release,41,42 which hinges on the mathematical product ε × Φ, where the molar absorptivity ε defines the efficiency of photon absorption at a given irradiation wavelength and a given PPG concentration, and where the quantum yield of photolysis (QY, Φ) defines the chance that a PPG in the excited state will release the payload (instead of returning to the ground state through other processes). The QY of a PPG is a crucial parameter, as it describes the ability of the PPG to carry out the very function it was designed for. A low QY significantly hinders the application of PPGs in fields where light-delivery is challenging, such as photopharmacology, where insufficient uncaging efficiency would require prolonged irradiation or use of higher intensity light, which is not feasible in a clinical setting.43 Furthermore, for the successful application of PPGs in fields that require near-instant release of a payload, such as time-resolved crystallography or materials chemistry, a high QY is crucial.44,45
A typical approach to increase the QY of heterolytic PPGs is to alter the nature of the leaving group, e.g. by employing a less Brønsted basic payload that efficiently stabilizes its incipient anionic charge in the CIP, increasing the ratio of the k1/k−1 rates and driving the uncaging process in the forward direction (Fig. 1b, CIP stabilization through anion stabilization),33,46 however it has to be noted that this way of QY engineering is unique to heterolytic PPGs, and does not necessarily work for PPGs that rely on different mechanisms of photolysis.47 When using this approach, often a correlation can be found between the pKa of the payload and the uncaging QY, with PPGs bearing payloads with lower pKa featuring higher QYs (Fig. 1c).33 However, this strategy involves direct chemical alterations of the payload, potentially affecting the desired application (for example, by diminishing the bioactivity of a biological molecule transiently “photocaged” by the PPG).
Recently, our group reported on a novel strategy to boost the QY of heterolytic PPGs by stabilizing the generated PPG cationic component of the CIP via hyperconjugation and/or positive charge delocalization (Fig. 1b, CIP stabilization through cation stabilization and Fig. 1d, compounds 3–6). Stabilization of the CIP was thus achieved by retarding the unproductive k−1 process, and allowing the CIP to undergo the productive chemical process k2 (Fig. 1b).48 In addition, we found a clear correlation between the experimentally recorded QY and the TD-DFT-calculated barriers of photolysis (k1) and CIP recombination (k−1), offering a strong computational predictive method for the design of efficient PPGs (Fig. 1e).
In our continuous effort towards understanding the full spectrum of chemical parameters that affect PPG QY, we summarize two puzzling observations from our own and other research groups' studies:
(1) Stabilizing the incipient anionic payload component of the CIP by employing a more weakly basic payload (i.e. with a stronger conjugate acid) does not always result in a more efficient uncaging process (Fig. 1c, stable anions).33,46 Instead, a complex interplay of payload conjugate acid pKaand overall payload topology (size, lipophilicity, etc.), and the polarity of the employed solvent49all appear to influence QY.
(2) Stabilizing the incipient cationic PPG component of the CIP as a tertiary cation (compound 6, i.e. methyl-versus allyl-substituted cations) resulted in a deviation from the correlation between experimental PPG QY and TD-DFT-calculated k−1/k1 energy barrier ratio, showing a higher than expected QY (Fig. 1e), and yielding a different product of PPG release: an alkene, which is a result of deprotonation instead of the usual solvent-trapped alcohol photoproduct.48 Similarly, coumarin- or BODIPY-PPGs engineered to give different photoproducts (e.g. via intramolecular capture of the cationic PPG component of the CIP) also appear follow this general trend of increased QY (Fig. 1c, alternative photoproducts).42,50,51
As a consequence of these puzzling observations, we hypothesized that the rate of CIP formation k1 and the rate of unproductive CIP recombination k−1, though critical parameters in engineering efficient QYs of PPGs,34,48 do not provide the full picture for overall PPG photolysis efficiency. Therefore, we set out to further investigate the productive process k2 that allows the payload to escape the CIP. While the crucial influence of the CIP intermediate itself on PPG photolysis has been recognized,33,34,37,48 its fate has so far not been explored as a key PPG design parameter.
Herein we outline the limitations of current PPG QY optimization methods and provide a robust framework for the future QY engineering methods for these valuable photochemical tools. We demonstrate the significant importance of the fate of the CIP in determining solvent- and payload-dependent photochemical properties, and showcase factors other than CIP stability that influence the QY of heterolytic diffusion-dependent PPGs. Furthermore, we provide more detailed insight in the photochemical behavior of the previously reported class of di-methyl-substituted coumarin PPGs 6, and demonstrate their unique photochemical properties with an overall uncaging mechanism independent of the molecular topology of the payload and the nature of employed solvent.
At the onset of this study, we hypothesized that the newly observed alkene photoproduct in the photochemical release of a model payload from 6 was due to an alternative CIP escape mechanism k2, that relied on deprotonation of the incipient coumarin cation from the payload anion (Fig. 2b). In the tertiary coumarin 6, as opposed to the secondary coumarins 1–5, this reaction would be promoted by three effects: firstly, the tertiary coumarin PPG has six β-protons available for deprotonation; secondly the elimination results in a double bond with a higher degree of substitution, making it more stable; thirdly, the increased steric bulk on the α-carbon resulting from the two methyl groups could potentially decrease the rate of the competing process, i.e. the water trapping of the cation.
We hypothesized that this new CIP escape mechanism would be an intramolecular -or ‘intra-ion pair’-process, circumventing the need for diffusion of both ions in the CIP. Being liberated from diffusion-dependency, we hypothesized that for tertiary coumarin, the k2 CIP escape energy barrier, and therefore the QY of 6, would not be affected by solvent polarity and payload size/hydrophobicity (Fig. 2b). Hypothetically, the QY of tertiary coumarin 6 would be significantly more independent from these factors than that of diffusion-dependent PPGs 1–5.
Throughout this paper, we measure the QYs of the PPGs using UV-Vis spectroscopy, following the consumption of the PPG substrate upon irradiation. Given that CIP recombination (k−1) happens on the nanosecond time-scale,34 it is undetectable by our spectrophotometer and CIP recombination events will appear as the substrate remaining intact. Only when the CIP successfully undergoes CIP escape (k2) the substrate will be consumed. Therefore, through following the efficiency of substrate conversion (i.e. the QY) the efficiency of CIP escape is also captured.
In a similar fashion, the photochemical QY of prenyl-coumarin 5 was also strongly influenced by the addition of other organic solvents (Fig. 3c, grey bars). The addition of ethanol already lowered the photochemical QY of 5 by 60%, while even less polar solvents further lowered its QY (Fig. 3c, grey bars). As expected, the same solvents had a minimal effect on the QY of tertiary coumarin 6 (Fig. 3c, red bars), hypothetically due to the solvent independent nature of its intramolecular CIP escape mechanism. Even the addition of comparatively apolar solvents such as dioxane showed a negligible effect on the QY of 6, and 6 could remarkably also be uncaged with identical efficiency in pure methyl tert-butyl ether (MTBE) (Fig. 3c, red bars). Unexpectedly, methanol and MeCN negatively affected the QY of 6 to a certain extent, despite their relatively high polarity.
Having established that the solvent polarity drastically influences the QY of diffusion-dependent PPG 5, we further compared the influence of the ionic strength on the uncaging efficiency of PPGs 5 and 6. We have chosen lithium perchlorate to control the ionic strength of the solution, since the perchlorate anion is non-nucleophilic, very weakly basic and also photochemically innocent. Especially the last aspect is of crucial importance, as illustrated by the difference in uncaging rates between the heterolytic PPG system Cy7-PPG derived from cationic heptamethine cyanine dyes,20 which we engineered to contain chloride as the counteranion to the cationic photocage, and its iodide salt counterpart reported by Stacko et al.53 This comparison reveals the non-innocent photochemical nature of the iodide counterion, which promotes oxygen-dependent processes, an effect well-known in photodynamic therapy.54,55 Interestingly, we have not observed any influence of the increase of ionic strength on the photouncaging efficiency of both PPGs 5 and 6, indicating that for both PPGs the process is independent on the presence of non-reactive ions (see ESI Section 3.5†).
Phenylacetic acid and Ibufenac were chosen as the larger payloads with increasing size and lipophilicity. Both these payloads are (phenyl)acetic acid analogues and are expected to have highly similar pKa values. This is crucial because payloads with a lower pKa value will stabilize the anionic component of the CIP intermediate and increase the QY through altering the k1 and k−1 energy barriers.33,48 Choosing payloads with (near) identical pKa allowed us to achieve isolated alteration of the k2 barrier, and therefore study its effect on the QY.
Steglich esterification of the carboxylic acids with prenyl-coumarin alcohol yielded prenyl-coumarins 7 and 9 bearing the two lipophilic payloads (Fig. 3d). For tertiary coumarin, the increased steric hindrance required longer reaction times at elevated temperature of the tertiary alcohol with the anhydrides of the respective payloads, to yield tertiary coumarins 8 and 10 (Fig. 3d). The quantum yields of photolysis of coumarins 5–10 were determined in a water/MeCN mixture. For diffusion-dependent prenyl-coumarin PPG (Fig. 3d, grey bars), the larger phenylacetic acid payload resulted in a significant decrease in its QY, presumably through increasing the diffusion k2 CIP escape barrier (Fig. 3d, 7). For prenyl-coumarin loaded with the Ibufenac payload (Fig. 3d, 9) the decrease in QY was even more pronounced (Fig. 3d, 5vs.9). The QY of tertiary coumarin uncaging was affected minimally by the increased size of its payload (Fig. 3d, red bars 6, 8, 10). No decrease in QY was found when moving from an acetic acid payload to the larger lipophilic phenylacetic acid payload (Fig. 3d, 6vs.8), while an Ibufenac payload minimally lowered the QY (Fig. 3d, 10). These results demonstrate the contrast between both CIP escape mechanisms, and illustrate that for tertiary coumarin uncaging does not depend on diffusion of both ions in the CIP, and therefore no relationship is found between the size or hydrophobicity of these ions and the photolysis QY. In contrast, for a PPG employing the conventional diffusion dependent CIP escape, these factors influence the QY significantly.
As the size of the payloads increases in the order acetic acid – phenylacetic acid – Ibufenac, so does their lipophilicity. Therefore, using these payloads alone, the combined influence of these two factors on the QY cannot be separated, and no conclusion can be made whether only one of them contributes to the drop in QY, or both. Therefore, we attempted to separate these factors through introducing polarity into the payloads, while keeping the size unaffected. We loaded prenyl-coumarin with aza-analogues of phenylacetic acid and Ibufenac (Fig. 3e, 11 and 12, respectively). While being very similar in size, the polarity of these payloads differs extensively. In agreement with our rationale regarding payload polarity and diffusion barrier height, for both payloads the introduction of a nitrogen atom increased the QY as compared to their carbon-analogues (Fig. 3e, blue bars compared to grey bars). These results demonstrate that irrespective of size, payload polarity influences the height of the diffusion barrier and the introduction of polarity has a favorable effect on the QY, likely through the reduction of this barrier. An interplay between both payload size and polarity determines the overall height of the k2 barrier and therefore both factors influence the final photochemical QY.
Overall, the results described in this section demonstrate the shortcomings of conventional diffusion-dependent heterolytic PPGs. The QYs of these PPGs can be strongly negatively affected by factors such as payload topology and polarity of the environment. The deprotonation mechanism – an intramolecular type of CIP escape – showcases itself as being superior, illustrated by the QY of tertiary coumarin that is minimally influenced by these factors. Furthermore, the tertiary coumarin PPGs 6, 8 and 10 showed excellent hydrolytic stability. No hydrolysis products were observed after incubation at 25 °C in water/DMSO mixtures (see ESI Section 3.10†). Lastly, uncaging chemical yields were determined for model compounds 5–8, that all showed high chemical yields ranging from 90–100% (see ESI Section 3.9†).
To expand our analysis that correlates photolysis QYs with the individual reaction barriers, regardless of the mechanism of CIP escape, we set out to include the mechanism of CIP escape k2 barrier in our analysis. Given that k2, like k1, is a productive step in the photolysis reaction, we included the value of the associated energy barrier into the denominator to reach a fraction that includes the height of all three barriers (eqn (1) (right)). This energy barrier ratio is related to the probability that a compound undergoing photoheterolysis travels its energy landscape productively. Therefore, we believed it should correlate to the photolysis QY as well. Eqn (1) shows the previously published fraction used for correlating the QY to the energy barriers (left)48 and a new energy barrier ratio that takes into account all steps in the photolysis process and respective energy barriers (right)
(1) |
Similar to our previous analysis,48 we sought out to compute all excited state energy barriers by TD-DFT calculations. For these calculations, we made the assumption that any productive photochemical processes would be occurring from the first singlet excited state (S1), which is consistent with what is known of coumarin-PPG photoheterolysis.33,48 TD-DFT has been previously used to compute photoheterolysis barriers on both the S1 and T1 excited state energy landscapes of coumarin- and BODIPY-PPGs.20,32,41,48,56 It should be noted that, while heterolysis (k1) occurs in S1,33 relaxation to S0 might happen at any time after CIP formation. Therefore, the extent to which CIP recombination (k−1) and CIP escape (k2) occur at S1 -if at all- is unclear. However, using TD-DFT we were unable to find a stable CIP in S0 and therefore unable to perform any computations on the k−1 and k2 energy barriers in S0. Therefore, our analyses are constrained to photochemical processes on the S1 potential energy landscape, and are to be seen as an approximation of the actual photolysis energy landscape. Nonetheless, the strong correlation of the k−1/k1 energy barrier ratio to the photolysis QY obtained in our previous analysis48 suggests that the computed heights of the energy barriers relative to one another among different PPGs are highly relevant.
To supplement our previous analysis that relied solely on k1 an k−1, the k2 barrier for the deprotonation that is proposed for the tertiary PPG 6 was also computed by TD-DFT in the first singlet excited state, and predicted to be 9.4 kcal mol−1 (see ESI Section 4†). For the other PPGs 1–5 that relied on CIP escape via diffusion of the chromophore cation and payload anion and water trapping, their k2 diffusion energy barriers could not be computed by TD-DFT. Therefore, we set out to numerically approximate their diffusion barrier using eqn (1). In this equation, we filled in the DFT calculated 9.4 kcal mol−1 for the k2 barrier height of 6. The value of the acetate diffusion barrier for all other coumarins was left as a variable. Through maximizing R2 numerically (Fig. 4), a k2 barrier height for diffusion of ∼14.6 kcal mol−1 was obtained.
Fig. 4 Exponential plot of the QY of photolysis versus the energy barrier ratio as reported in our previous analysis48 (left). Including in the energy barrier ratio allows for an estimation of the diffusion barrier through R2-optimization (right). |
Additionally, performing the same analysis without including the calculated k1 heterolysis barriers also gave an excellent correlation with the QY, resulting in an R2 of 0.9975. Using this method, a slightly lower k2 acetate diffusion barrier was found (see ESI, Section 5†).
The obtained correlation of the ratio of the photolysis energy barriers with the QY is strong (R2 > 0.99). Still, it has to be noted that the predicted height of the k2 diffusion barrier relies on the accuracy of all energy barriers computed by TD-DFT, which as discussed above, were all performed exclusively for the first singlet excited state. Therefore, the obtained diffusion barrier height of 14.6 kcal mol−1 should be seen as a qualitative estimation, rather than a quantitative one.
(2) |
Eqn (2) shows the formula describing the relationship between the QY of coumarin PPGs 1–6 and the energy barrier ratio. This formula connects the QY to the photolysis energy barriers alone, and does not account for any other photophysical processes that may occur after excitation. For example, eqn (2) does not account for excited state relaxation processes such as fluorescence and internal conversion, the rate of which do influence the QY of competing photolysis. If the rates of these processes would differ significantly between each coumarin variant 1–6, we would be unable to correlate the photolysis QYs solely to the computed photolysis energy barrier ratio. However, judging from the experimentally obtained strong correlation (Fig. 4 and previous work48), we conclude that it is likely not the case. This is supported by the fluorescence QYs we determined for the parent alcohols of 1–6, that lack the acetic acid payload and therefore the PPG functionality itself. These alcohols all displayed similar fluorescence QYs, ranging from 9–12% (see ESI Section 3.8†). Therefore, interestingly, modification of the coumarin alpha-carbon appears to not affect the rate of fluorescence -and likely other photophysical processes- to a significant extent. A possible explanation for this could be that the performed synthetic modifications do not result in alteration of the chromophore itself, but importantly stabilize the cation that is generated exclusively after photoheterolysis and not in any of the other competing processes.
Since our analysis solely focused on the photolysis barriers, we were interested in testing the reliability of the newly found correlation. We hypothesized that if the obtained formula (eqn (2)) adequately describes the system, it could allow us to predict the photolysis QY of an isotopologue PPG of 6. The compound we designed for this purpose was deuterated-coumarin 6-D6 (Fig. 5a). Being a tertiary coumarin, we expected it to undergo deprotonation as CIP escape mechanism. And by virtue of being an isotopologue of 6, 6-D6 would have identical (photo)-chemical and topological properties as 6. The primary isotope effect of deuteron rather than proton abstraction would result in a larger k2 barrier for 6-D6 than for tertiary coumarin 6, effectively resulting in a third k2 barrier (diffusion, H+-abstraction, D+-abstraction). DFT calculations confirmed this, yielding a predicted energy barrier of deuteron abstraction of 11.5 kcal mol−1, as compared to 9.4 kcal mol−1 for proton abstraction of 6. Eqn (3) shows the energy barrier ratio ‘x’ as calculated for compound 6-D6 based on the DFT-calculated height of the three energy barriers, and the predicted QY of 6-D6 through filling in ‘x’ in the previously established formula correlating the QY to the energy barrier ratio.
(3) |
Fig. 5 (a) Synthetic scheme for the formation of deuterated tertiary coumarin 6-D6. (b) Predicted QY of coumarin 6-D6 using eqn (2), and measured QY of 6-D6 (55 µM, water 2.7% MeCN, 25 °C, shown is the average and SD-value of a triplicate measurement). |
Feeding the DFT calculated energy barriers for k1, k−1 and k2 to eqn (2) allowed us to predict the photolysis QY of 6-D6 at 7.06% (eqn (3)), roughly 30% lower than the QY of its protonated analogue 6. To our delight, experimental determination of the QY yielded an almost identical number of 7.0 ± 0.4%.
The fact that eqn (2) is highly predictive indicates that our analysis comprises the most important factors that determine coumarin PPG QY. We expect equations analogous to eqn (2) to allow for the optimization of QY in the design of new PPGs, through revealing trends correlating PPG design to expected QY. These trends can be established using acetate as a model cargo, but we expect them to be transferable to other types of cargos with different diffusion barriers. Our general approach could be expanded to other heterolytic PPGs and aid in their efficient design, bearing in mind that the specific equation would have to be adapted due to those PPGs having different rates of competing photophysical processes.
Fig. 6 Plot of eqn (2), displaying the energy barrier ratio (x-axis) vs. the QY (y-axis). An increase of the k2 barrier height is predicted to lower the QY of efficient diffusion-dependent PPGs drastically (prenyl-coumarin with acetic acid vs. a large payload (red sphere)). The same relative increase of the k2 barrier height has a minimal effect on the QY of inefficient PPGs (e.g., primary coumarin). |
Therefore, when scientific fields employing PPGs as photochemical tools move towards the use of more efficient PPGs that still rely on diffusion as CIP escape, the size of the payload becomes a highly important factor and a larger payload will likely lower the QY significantly, negating the efforts of QY optimization. Since almost all relevant payloads are larger than acetic acid (note that most payloads used in the field of photopharmacology are larger bioactive compounds) this effect for diffusion-dependent PPGs will become highly relevant. It emphasizes the importance of engineering the CIP escape mechanism k2 and further highlights the advantage that PPGs with an intramolecular type of CIP escape – like 6, that circumvent these issues- have over conventional diffusion-dependents PPGs, most notably when moving towards higher QYs. For example, for the field of photopharmacology, the development of water-soluble PPGs with an intramolecular type of CIP escape and a high QY will be a major step forward.
Footnotes |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d3sc05725a |
‡ These authors contributed equally to this work. |
This journal is © The Royal Society of Chemistry 2024 |