John P.
Otto‡
a,
Lili
Wang‡
a,
Igor
Pochorovski
b,
Samuel M.
Blau
c,
Alán
Aspuru-Guzik
cd,
Zhenan
Bao
*b,
Gregory S.
Engel
*a and
Melanie
Chiu§
*b
aDepartment of Chemistry, University of Chicago, Chicago, IL 60637, USA. E-mail: gsengel@uchicago.edu
bDepartment of Chemical Engineering, Stanford University, Stanford, CA 94305, USA. E-mail: zbao@stanford.edu; melanie.chiu@stonybrook.edu
cDepartment of Chemistry and Chemical Biology, Harvard University, Cambridge, MA 02138, USA
dSenior Fellow, Canadian Institute for Advanced Research, Toronto, Ontario M5G 1Z8, Canada
First published on 15th March 2018
Förster Resonance Energy Transfer (FRET) is the incoherent transfer of an electronic excitation from a donor fluorophore to a nearby acceptor. FRET has been applied as a probe of local chromophore environments and distances on the nanoscale by extrapolating transfer efficiencies from standard experimental parameters, such as fluorescence intensities or lifetimes. Competition from nonradiative relaxation processes is often assumed to be constant in these extrapolations, but in actuality, this competition depends on the donor and acceptor environments and can, therefore, be affected by conformational changes. To study the effects of nonradiative relaxation on FRET dynamics, we perform two-dimensional electronic spectroscopy (2DES) on a pair of azaboraindacene (BODIPY) dyes, attached to opposite arms of a resorcin[4]arene cavitand. Temperature-induced switching between two equilibrium conformations, vase at 294 K to kite at 193 K, increases the donor–acceptor distance from 0.5 nm to 3 nm, affecting both FRET efficiency and nonradiative relaxation. By disentangling different dynamics based on lifetimes extracted from a series of 2D spectra, we independently observe nonradiative relaxation, FRET, and residual fluorescence from the donor in both vase to kite conformations. We observe changes in both FRET rate and nonradiative relaxation when the molecule switches from vase to kite, and measure a significantly greater difference in transfer efficiency between conformations than would be determined by standard lifetime-based measurements. These observations show that changes in competing nonradiative processes must be taken into account when highly accurate measurements of FRET efficiency are desired.
In Förster's model, the coupling between the donor and acceptor electronic states is treated as a second-order perturbation in the electronic coupling. The interaction is governed by the transition dipole moments of the excited electronic states. The coupling is sufficiently weak (and the transfer rate therefore sufficiently slow) that the donor relaxes to thermal vibrational equilibrium in its electronic excited state before transfer can occur. Although no absorption or emission occurs in the transfer itself (FRET is nonradiative), the same transition dipoles that are responsible for donor emission and acceptor absorption also drive the transfer, so resonance can be determined by absorption and emission spectra, and transfer competes with donor fluorescence and nonradiative relaxation.
Experimentally, FRET efficiency can be determined by photon counts of donor or acceptor fluorescence or by changes in the donor's decay lifetime.17,29 The latter approach is often preferred due to relative immunity to cross-talk and concentration artifacts,7,16,22 and is based on kinetic competition of decay processes. In this scenario, efficiency can be expressed as
(1) |
1/τD = kD = kD,f + kD,nr | (2a) |
1/τDA = kDA = kD,f + kD,nr + ktr | (2b) |
To investigate the interactions between nonradiative relaxation and FRET, we sought a pair of chromophores that undergo both excited-state processes. Moreover, we sought a means of exerting control over both the distance separating the chromophore pair and the vibrational bath that defines the local environment of each chromophore. Quinoxaline-bridged resorcin[4]arene cavitands emerged as ideal molecular scaffolds for fulfilling the latter requirement, as these conformational switches can be toggled between a pursed vase state and a ringent, kite state using various external stimuli, including temperature.31,32 Due to solvation-related entropic effects, the cavitand exists in the vase conformation at room temperature and above; conformational switching commences at ca. 223 K, and at 193 K, only the kite conformation is observed.33 To adapt this scaffold for our studies, we sought to append each of a pair of chromophores to opposing walls of the cavitand.34–36 We selected a pair of arylethynyl-functionalized azaboraindacene (BODIPY) chromophores for this purpose because they exhibit the requisite photophysical properties,37–39 and because both a donor and an acceptor chromophore are conveniently accessible from a single, divergent synthetic route.40
Given the ability of 2DES measurements to elucidate ultrafast transfer and relaxation dynamics in both biological41–43 and synthetic systems,44–46 we sought to study the BODIPY-functionalized cavitand using this technique. By separating the signal based on time dynamics, we disentangle three relaxation pathways (nonradiative relaxation, energy transfer, and residual fluorescence) and isolate their lineshapes from the donor in both conformations (kite at 193 K and vase at 294 K). Based on our signal isolation analysis, we observe that the contribution of nonradiative relaxation to the overall donor decay changes with conformation. We show that, due to this change in nonradiative relaxation, fluorescence lifetime measurements underreport the transfer efficiency difference between kite and vase conformations, highlighting the need to test assumptions about nonradiative processes when measuring FRET in complex environments.
Temperature-dependent conformational behavior of cavitand 1 was characterized by variable-temperature 1H NMR,33 absorption, and fluorescence spectroscopies.34 As shown in Fig. 1, the 1H NMR spectrum of cavitand 1 at 323 K exhibits two sets of sharp multiplets between 5.7 and 5.8 ppm, corresponding to the resorcin[4]arene methine protons (★) and indicative of the vase conformation. Cooling the sample to 223 K resulted in a diagnostic upfield shift of these methine proton signals to ca. 3.8 ppm due to increased shielding from the proximity of the aromatic wall flaps in the kite conformation. Consistent with previous reports, the vase-to-kite transition was also observed upon addition of acid (see Fig. S27 in the ESI†).34
The temperature-dependent conformational switching of cavitand 1 was corroborated by fluorescence and absorption spectroscopies, as depicted in Fig. 2. The vase conformation of cavitand 1 predominates at 294 K (Fig. 2a), evidenced by a blue-shift of the absorption maximum of cavitand 1 relative to that of the dye mixture; interactions between the donor and acceptor dyes in the vase conformation of cavitand 1 result in shifts of their ground- and excited-state energies. At 193 K, the absorption spectrum of cavitand 1 resembles that of the dye mixture, which is consistent with the weak interaction between the donor and acceptor dyes in the kite conformation (Fig. 2b). Fluorescence spectra of cavitand 1 at 294 K and 193 K (excitation at 525 nm for both conformations) show a clear increase in donor emission at 610 nm with a concurrent decrease in acceptor emission at 670 nm as the temperature is lowered. These fluorescence changes indicate a decrease in FRET efficiency and an increase in the distance between the donor and acceptor dye at lower temperature, which is consistent with conversion of vase to kite with decreasing temperature.
Computational studies provided further structural insights into the two conformations of cavitand 1. The geometries of the vase and kite conformations were optimized using density functional theory calculations (B3LYP/3-21G(d)) in the Gaussian 09 software package,49 and are depicted in Fig. 2c and d, respectively. In both conformations, terminal C6H13 groups on the cavitand base have been replaced with CH3. In the vase conformation, the donor and acceptor are sterically constrained to an orientation that promotes dipole–dipole coupling. The donor–acceptor distance in the vase conformation is ca. 0.5 nm, which is too close for independent solvation of the two chromophores, suggesting that they interact vibrationally, as well as electronically. In contrast, the donor and acceptor are separated by ca. 3 nm in the kite conformation, providing the chromophores with relative orientational freedom. Although the donor and acceptor are separated by less than R0, they interact weakly enough that the approximations made by FRET theory accurately describe their interaction. Having successfully synthesized cavitand 1 and confirmed that its conformation and transfer dynamics could be externally controlled by temperature, we next sought to probe the excited-state dynamics of cavitand 1, focusing on the interactions between nonradiative relaxation and FRET processes.
2D spectra are frequency–frequency correlation maps, describing the optical response of a system (recorded on ωt axis) at a series of waiting times, T, after initial excitation at ωτ. Signals along the diagonal indicate features detected at the same energy as their excitation, while signals above or below the diagonal indicate coupling or population transfer between excited states. Additionally, the lineshapes of 2D signals provide information on the mechanisms that lead to broad spectral signals.50 Inhomogeneous broadening (extension along the diagonal) gives insight into static variation in chemical environments, which perturb the excitation energies of individual chromophores.51 In contrast, homogeneous broadening (which determines antidiagonal linewidth) occurs when frequency fluctuations due to bath interactions occur quickly compared with the experimental timescale, such that an individual chromophore experiences all configurations of its local environment during a given measurement and thus represents the entire ensemble.52 In the homogeneous limit, one can estimate the dephasing time of a transition (the time required to lose phase correlation in the ensemble of sample molecules) by its homogeneous linewidth.53
At T = 20 fs, 2D spectra of both vase (Fig. 3a) and kite (Fig. 3d) show an inhomogeneously broadened peak along the diagonal at roughly 17100 cm−1. This signal corresponds to the initial absorption of the donor; at 20 fs the system has not yet relaxed past this starting point. Vibrational relaxation (aka Stokes shift) from this initial absorption occurs within the first 100 fs, leading to the below-diagonal cross peak in all subsequent spectra (see ESI† for more detail).
At T = 260 ps (Fig. 3b and e), two dominant features are present in both 2D spectra—a diagonal peak at 17100 cm−1 (white star) and a crosspeak at 17100 cm−1 (585 nm) in ωτ and 16400 (610 nm) cm−1 in ωt (black star). The diagonal feature at T = 260 ps is due to ground state bleach on the donor, because Stokes shift from the initial excitation has completed by this time. The below-diagonal crosspeak corresponds primarily to stimulated emission from the Stokes-shifted excited state of the donor. Along the anti-diagonal direction, all features are very broad, indicating the homogeneity of excited donor molecules. The width of these homogeneously broadened features is caused by rapid dephasing on the scale of tens of fs (Fig. S30†). A red shift in early time dynamics of the complex occurs on the same timescale, (see Fig. S31†) and the processes driving the red shift may also be the cause of this rapid dephasing. Further discussion of dephasing lifetimes and ultrafast dynamics can be found in the ESI.†
The waiting time decay dynamics were fit to the following equation
S(T) = A1e−T/τ1 + A2e−T/τ2 + A3 | (3) |
The left (grey) column in Fig. 4 shows a fast (∼25 ps) decay for cavitand 1 in the vase (294 K, top) and kite (193 K, bottom) conformations. These lifetimes are consistent with reports of a decay in the 15–30 ps range for analogous meso-aryl substituted boron dipyrrin structures.37,38 We assign this fast decay to a nonradiative, vibrationally driven relaxation arising from a combination of flexing in the dipyrrin backbone and rotation of the aryl group, as suggested by Li et al.38 and Kee et al.37 The multi-exponential nature of the decay means that this process (corresponding to the fastest exponential decay) does not compete kinetically with the processes observed in the other two exponential decays. This observation may be explained by a multi-well potential energy surface for the excited donor, as illustrated in Fig. 4 and observed by Li et. al.38
While the lifetimes of this nonradiative relaxation pathway are very similar for both conformations, the corresponding lineshapes in the 2D lifetime map differ notably. In the kite spectrum (Fig. 4, lower left plot, grey), the signal is primarily due to one homogeneously broadened feature at 17000 cm−1. The antidiagonal full width at half maximum (FWHM) of this feature corresponds to a homogeneous process with a dephasing lifetime of 48 fs, coincident with the timescale of relaxation within the electronic excited state as seen via Stokes shift in fs-resolved 2DES (Fig. S31†). By comparison, the vase spectrum (Fig. 3, upper left plot) shows an inhomogeneously broadened diagonal feature at 17100 cm−1, as well as a crosspeak feature at an excitation frequency of 17000 cm−1 and detection frequency of 16400 cm−1. Inhomogeneous broadening indicates that subsections of the sample population have differing energetic environments, such as those arising from conformational variation, that persist on the timescale of the experiment.
The middle (orange) column in Fig. 4 represents the strength of signal that decays with a lifetime of several hundred ps. This decay lifetime is not observed in fluorescence lifetime measurements of the donor solution. We assign this signal to energy transfer between donor and acceptor (FRET). The energy of the diagonal feature matches with GSB from the donor, while the energy of the below-diagonal feature matches SE from the Stokes-shifted excited donor, as well as the absorption of the acceptor. Both features decay due to FRET transfer to the acceptor. Because the initial vibrational relaxation in the donor excited state is complete within 100 fs (Fig. S31†), all signal at the diagonal observed on this timescale is due to GSB of the donor, meaning that only the loss of donor excitation is observed. The below-diagonal features may contain signal from the acceptor and have longer lifetimes than features on the diagonal (Table TS4 in ESI†). Hence, the diagonal feature is used for donor lifetime measurements. This FRET transfer pathway contributes to the majority of the signal, particularly in the vase conformation where the chromophores are in close proximity to one another.
The right (purple) column in Fig. 4 shows the portion of the 2D signal that does not decay on the timescale of 800 ps. We assign this signal to residual donor fluorescence based on the measured fluorescence lifetime of the donor (2–3 ns, ESI Table TS2†). The double peak lineshape of this signal can be described by a combination of GSB at the diagonal (corresponding to the original absorption) and SE at the crosspeak (corresponding to Stokes-shifted emission). Simulated spectra generated from the absorption and fluorescence of the donor alone (Fig. S32†) qualitatively match the lineshape of this long-lived signal for both conformations, further supporting our assignment. Previous work on FRET between orthogonal dyes has highlighted that low frequency vibrations (analogous here to the difference in orientational freedom between vase and kite) can alter transfer efficiency dramatically.59
While we identify the second of the three lifetimes as corresponding to energy transfer, it should be noted that the measured lifetime is distinct from the time actually required for transfer. FRET relies upon the same dipole coupling between ground and excited states that leads to donor fluorescence, meaning that fluorescence always occurs along with FRET and that the branching ratio between them is determined by their relative rates. As such, the measured decay in excited state population will have a rate k2D such that
k2D = ktr + kd,0 | (4) |
In addition to determining transfer lifetimes, we can calculate the branching ratio of FRET in this system using eqn (1), where τDA is the lifetime we assign to energy transfer in DAS and τD is the fluorescence decay lifetime of the donor in solution, measured via time-correlated single photon counting (TCSPC). This method results in a calculated transfer efficiency of 92.1 ± 1.3% transfer for the vase conformation and 78.0 ± 2.9% for the kite conformation. This approach only accounts for kinetic competition, however, and therefore only determines the efficiency of the subpopulation measured in DAS. If we use the peak diagonal intensities from DAS to approximate branching between the three subpopulations, we obtain overall FRET efficiencies of 68.7 ± 1.1% for the vase conformation and 44.3 ± 1.8% for the kite conformation, a significant reduction in overall values and a dramatically greater difference between conformations than the difference measured viaeqn (1).
The results above highlight how conformation-dependent variations in nonradiative relaxation processes affect the accuracy of lifetime-based FRET measurements. These effects should be considered in situations where the chromophore structure(s) or their relative orientations result in either: (1) breakdown of the ideal dipole approximation (IDA); or (2) appreciable competition in the excited state population between FRET and nonradiative relaxation. The former is relevant to chromophore structures that exhibit delocalized transition densities60,61 or to systems in which chromophore orientations are restricted such that isotropic dipole distributions that might mitigate IDA breakdown cannot be achieved.62,63 The latter situation might arise in donor chromophore structures that exhibit either modest quantum yields or conformational freedom that exacerbates nonradiative relaxation processes.37 In addition to considering the effect of variations in nonradiative relaxation on FRET measurements, these variations might also be exploited productively as probes for changing viscosity in cell environments via motion-induced change in emission measurements (MICE).64
UV/vis absorption spectra were recorded using a Cary 5000 dual-beam UV/VIS/NIR spectrophotometer. Steady state fluorescence was measured via a Horiba Fluorolog 3, and fluorescence quantum yields were recorded using the Horiba Quanta-Phi attachment. Time-domain lifetimes were measured on a ChronosBH lifetime fluorometer (ISS, Inc.) using Time-Correlated Single Photon Counting (TCSPC) methods. The fluorometer contained Becker-Hickl SPC-130 detection electronics and an HPM-100-40 Hybrid PMT detector. Tunable picosecond pulsed excitation at 560 nm was provided by a Fianium SC400-2 supercontinuum laser source with integrated pulse picker and AOTF. Emission wavelengths were selected with bandpass filters (Chroma HQ615/85 nm). The Instrument Response Function (IRF) was measured to be approximately 120 ps FWHM in a 1% scattering solution of Ludox LS colloidal silica. Multi-component exponential decay lifetimes were fit via a forward convolution method in the Vinci control and analysis software.
Two-dimensional electronic spectra were acquired using a homebuilt instrument that has been described previously.45 A Ti:sapphire oscillator and regenerative amplifier (Micra and Legend Elite, Coherent) generate a 5 kHz, 100 fs pulse train that is directed into a noncollinear optical parametric amplifier (NOPA, TOPAS). The output of the NOPA is then compressed using a prism pair followed by a spatial light modulator array (FemtoJock, Biophotonics). The resulting pulse (12 fs, 530–650 nm FWHM) is then split into four via a beamsplitter and a transmissive grating. The four beams are directed to the sample with a boxcar geometry. The time delays are controlled by paired glass wedges (coherence time, τ) and a translational stage (waiting time, T). Beams 1, 2 and 3 interact with the sample to generate the third-order signal. Beam 4 is attenuated by three orders of magnitude and serves as local oscillator for heterodyne detection. Signal, along with the local oscillator, is directed to a spectrometer and the resulting interferogram is recorded on a CCD array (Andor). All 2D data in the main text were collected in the same 24 h period in the same cryostat, varying only the temperature between the two measurements.
Footnotes |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/c8sc00818c |
‡ These authors contributed equally to this work. |
§ Current address: Department of Chemistry, Stony Brook University, Stony Brook, NY 11794, USA. |
This journal is © The Royal Society of Chemistry 2018 |