Journal Pre-proof Measuring conformational equilibria in allosteric proteins with time-resolved tmFRET William N. Zagotta, Eric G.B. Evans, Pierce Eggan, Maxx H. Tessmer, Kyle D. Shaffer, E. James Petersson, Stefan Stoll, Sharona E. Gordon PII: S0006-3495(24)00091-2
DOI: https://doi.org/10.1016/j.bpj.2024.01.033 Reference: BPJ 12946 To appear in: Biophysical Journal Received Date: 9 October 2023 Accepted Date: 29 January 2024
Please cite this article as: Zagotta WN, Evans EGB, Eggan P, Tessmer MH, Shaffer KD, Petersson EJ, Stoll S, Gordon SE, Measuring conformational equilibria in allosteric proteins with time-resolved tmFRET, Biophysical Journal (2024), doi: https://doi.org/10.1016/j.bpj.2024.01.033. This is a PDF file of an article that has undergone enhancements after acceptance, such as the addition of a cover page and metadata, and formatting for readability, but it is not yet the definitive version of record. This version will undergo additional copyediting, typesetting and review before it is published in its final form, but we are providing this version to give early visibility of the article. Please note that, during the production process, errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain. © 2024 Biophysical Society.
Measuring conformational equilibria in allosteric proteins with time-resolved tmFRET William N. Zagottaa, Eric G. B. Evansa,b, Pierce Eggana, Maxx H. Tessmerb, Kyle D. Shafferc, E. James Peterssonc, Stefan Stollb, Sharona E. Gordona
a Department of Physiology and Biophysics, University of Washington, Seattle, Washington 98195 b Department of Chemistry, University of Washington, Seattle, Washington 98195 c of Department of Chemistry, University of Pennsylvania, 231 South 34th Street, Philadelphia, PA 19104
re lP ur and na
Sharona E. Gordon, Ph.D. Dept. of Physiology and Biophysics Box 357290 University of Washington Seattle, WA 98195-7290 Tel: (206) 616-4861 Email: seg@uw.edu -p ro Address correspondence to: Jo
William N. Zagotta, Ph.D. Dept. of Physiology and Biophysics Box 357290 University of Washington Seattle, WA 98195-7290 Tel: (206) 685-3878 Email: zagotta@uw.edu Running Title: Distributions with time-resolved tmFRET
ABSTRACT lP re -p ro of
Proteins are the workhorses of biology, orchestrating a myriad of cellular functions through intricate conformational changes. Protein allostery, the phenomenon where binding of ligands or environmental changes induce conformational rearrangements in the protein, is fundamental to these processes. We have previously shown that transition metal Förster resonance energy transfer (tmFRET) can be used to interrogate the conformational rearrangements associated with protein allostery and have recently introduced novel FRET acceptors utilizing metal-bipyridyl derivatives to measure long (>20 Å) intramolecular distances in proteins. Here, we combine our tmFRET system with fluorescence lifetime measurements to measure the distances, conformational heterogeneity, and energetics of maltose binding protein (MBP), a model allosteric protein. Time-resolved tmFRET captures near-instantaneous snapshots of distance distributions, offering insights into protein dynamics. We show that time-resolved tmFRET can accurately determine distance distributions and conformational heterogeneity of proteins. Our results demonstrate the sensitivity of time-resolved tmFRET in detecting subtle conformational or energetic changes in protein conformations, which are crucial for understanding allostery. In addition, we extend the use of metal-bipyridyl compounds, showing Cu(phen)2+ can serve as a spin label for pulse dipolar electron paramagnetic resonance (EPR) spectroscopy, a method which also reveals distance distributions and conformational heterogeneity. The EPR studies both establish Cu(phen)2+ as a useful spin label for pulse dipolar EPR and validate our time-resolved tmFRET measurements. Our approach offers a versatile tool for deciphering conformational landscapes and understanding the regulatory mechanisms governing biological processes.
Responsible for the regulation of virtually all biological processes, allosteric proteins are fundamental to all life. To understand the mechanisms for their vital functions, we must determine the structure and energetics of these proteins under physiological conditions. In this work, we have developed a fluorescence method that allows for the simultaneous measurement of protein structure and energetics under physiological conditions. These studies pave the way for future advances in physiology and medicine.
INTRODUCTION Protein allostery plays a pivotal role in the regulation of virtually all biological processes. In response to the binding of specific molecules or changes in environmental conditions, allosteric proteins undergo distinct changes in structure that regulate the protein’s activity or interaction with other molecules. This dynamic behavior allows proteins to function as molecular switches, orchestrating a wide range of biological functions such as enzymatic catalysis, signal transduction, gene regulation, and cellular motility. The mechanism of allostery involves a choreography of the protein’s structure and energetics. A stereotypical allosteric protein might have two conformations, a resting state and an active state (Figure 1A). Whereas the active state might be energetically unfavorable (positive ΔG) in the absence of ligand, it becomes more favorable (negative ΔG) in the presence of ligand. In this way, the conformation of the protein, and therefore its activity, is coupled to the binding of a ligand.
One general way to determine this ΔG is to measure conformational equilibria based on distance distributions between a pair of probes, attached to the protein of interest, which undergo a change in separation between resting and active states (Figure 1A). The hypothetical distance distributions in Figure 1A show that, in the absence of ligand, or “apo” condition, the resting conformation (longer distance) dominates, and, in saturating ligand or “holo” condition, the active conformation (shorter distance) dominates. The relative occupancy of the two conformations can be quantified by the relative area of the probability distribution for each conformational state. The ratio of the areas reveals the equilibrium constant and therefore the ΔG for the transition between the apo and holo states. These distance distributions reveal two types of heterogeneity that must be considered for any allosteric protein: 1) heterogeneity in functional state, as both resting and active conformations are present in both the absence and presence of ligand, and 2) for any given state, heterogeneity in the protein backbone and probe rotameric ensembles produce a distribution of distances between probes.
While Förster resonance energy transfer (FRET) has been used as a “molecular ruler” to measure distances in proteins, standard steady-state FRET experiments provide just a single number, the apparent FRET efficiency, from which one can calculate only a single weighted-average distance (1,2). Time-resolved FRET experiments, however, generate richer data, from which the distribution of distances can be recovered (2-11). Time-resolved FRET utilizes fluorescence lifetimes, the latencies between absorption and emission of photons from a fluorophore. In the simplest cases, fluorescence lifetimes are single exponentially distributed with a time constant of a few nanoseconds. The presence of a FRET acceptor accelerates the fluorescence decay of the donor fluorophore in a manner that is highly dependent on the distance between the donor and acceptor (Figure 1B). If, for example, there are two states with different distances between the donor and acceptor, the decay will be double exponential with the fraction of each component representing the prevalence of that state. Importantly the interconversion among states is generally slower than the nanosecond fluorescence lifetime, so time-resolved FRET captures a near-instantaneous snapshot of the distance distribution without the averaging of distances, in contrast to steady-state or single-molecule fluorescence methods. The utility of time-resolved FRET for deciphering protein dynamics has been limited by several experimental factors: 1) Labeling proteins site-specifically with donor and acceptor fluorophores can be a challenge. 2) The large size of most visible-light fluorophores and the length of the linkers that connect the fluorophores to the protein make it difficult to discriminate backbone dynamics and energetics from those of the fluorophore/linker. 3) Most fluorophores demonstrate multi-exponential lifetimes, which makes analysis of FRET data more complicated. 4) The fluorescence lifetimes of most dyes in common use are in the few-nanosecond range, giving a limited dynamic range for measuring time-resolved FRET. We recently developed a novel system for time-resolved FRET that overcomes these limitations by combining a noncanonical amino acid fluorophore donor and a transition metal ion acceptor (12). Here we employ this novel system with metal-bipyridyl acceptors, developed in a companion paper in this issue (13), to measure longer distance distributions in a model protein, maltose binding protein (MBP). We show that time-resolved FRET can quantify both the heterogeneity of a given conformational state and the energetics that govern the distribution of a protein among conformational states, collectively referred to as protein dynamics.
MATERIALS AND METHODS Expression and purification of MBP For tmFRET experiments, the expression and purification of MBP was done as described previously (13). Briefly, MBP TAG constructs with a C-terminal twin-strep tag were cotransfected with a plasmid containing the AcdA9 aminoacyl tRNA synthetase and its cognate tRNA (14) in BL-21(DE3) cells. Cultures were induced in the presence of 0.6 mM Acd in the media, and MBP was purified on a Streptactin column (IBA Life Sciences, Göttingen, Germany) column.
For RIDME experiments, dual cysteine constructs of MBP-295C-211C and MBP-322C-278C with Nterminal 6×His tags were expressed from a pETM11 vector in E. coli C43(DE3) and subsequently purified by Co2+ affinity chromatography as previously described (15). The 6xHis tag was removed by incubation with a 1:50 (TEV:MBP) weight ratio of TEV protease (4 h at room temperature, then 12 h at 4 °C) in K+Tris buffer (130 mM KCl, 30 mM Tris, pH 7.4) containing 0.5 mM EDTA and 1 mM TCEP. The reaction was desalted into K+-Tris buffer (pH 7.4) plus 5 mM imidazole and 50 μM TCEP and further purified by reverse IMAC over TALON resin. Flow through containing cleaved MBP was supplemented with 1 mM TCEP and 5 mM EDTA, concentrated (30 kDa MWCO), and stored at 4 °C.
PhenM and [Ru(bpy)2phenM]2+ were prepared as stocks in DMSO and used within minutes of final dilution into aqueous solution. 2 M hydroxylamine hydrochloride in water was prepared for use in Fe2+ experiments and used for only one day. FeCl2 was prepared as a 100 mM stock with 1 M hydroxylamine hydrochloride in water and made fresh for each experiment day. CuCl2 was prepared as a 100 mM stock in water.
ur [Ru(bpy)2phenM]2+ labeling Jo
Time-resolved fluorescence measurements required a higher protein concentration (i.e., higher concentration of donor) and therefore a higher concentration of [Ru(bpy)2phenM]2+ was required for labeling. To reduce background absorption due to [Ru(bpy)2phenM]2+ in solution, we first labelled with [Ru(bpy)2phenM]2+ and then column purified the protein. Specifically, 100 mM [Ru(bpy)2phenM]2+ stock in DMSO was diluted to 1 mM in the concentrated protein in K+-Tris buffer (pH 7.4). After 10 minutes, the solution was passed over a Bio-Rad Micro Bio-Spin 6 column that had been equilibrated with K+-Tris buffer (pH 7.4) to remove unreacted label. K+-Tris buffer (pH 7.4) was also used to elute the protein.
[Fe(phenM)3]2+ labeling For [Fe(phenM)3]2+ experiments, a solution of 920 µM FeCl2 in 9.2 mM hydroxylamine hydrochloride with 2.3 mM phenM was prepared in water. After recording the donor-only fluorescence lifetime, this [Fe(phenM)3]2+ solution was added to the protein drop to a final concentration of 76.8 µM Fe2+, 768 µM hydroxylamine hydrochloride, and 192 µM phenM. For these brief experiments, additional hydroxylamine hydrochloride was not required to prevent oxidation of Fe2+ to Fe3+.
[Fe(phenM)]2+ labeling For [Fe(phenM)]2+ experiments, a 20 mM phenM stock in K+-Tris buffer (130 mM KCl, 30 mM Tris, pH 8.3) was added to MBP protein to achieve a final concentration of 2 mM phenM. After 10 minutes, the 4
solution was passed over a Bio-Rad Micro Bio-Spin 6 column that had been equilibrated with K+-Tris buffer (pH 8.3) to remove unreacted phenM label. K+-Tris buffer (pH 8.3) was also used to elute the protein. After measuring the lifetime of the purified protein in the absence Fe2+, we added Fe2+ (with a 10-fold excess of hydroxylamine hydrochloride) to a final concentration of 800 μM.
Measurement of fluorescence lifetime using FLIM -p ro of
The theory underlying our FRET measurements with fluorescence lifetimes is well described elsewhere (2). Briefly, FRET decreases the fluorescence lifetime of a donor fluorophore by providing an additional path by which an excited state electron can lose its energy. When using a pulse excitation source and measuring fluorescence in the time domain, the decrease in lifetime is readily apparent as a faster decay in fluorescence intensity after excitation (Figure 1B). When using a frequency (𝜔) modulated excitation source, the lifetimes of donor in the absence and presence of acceptor are determined from the phase delays (the phase shift between the excitation and emission, 𝜑𝜔 ) and modulation ratios (the fractional decrease in the amplitude of the emission, 𝑚𝜔 ) at each frequency (c.f. Figure 2A). With our frequency domain instrument, the frequency dependence of both 𝜑𝜔 and 𝑚𝜔 are required to resolve complex lifetimes. A similar analysis can be performed using a time domain instrument (2).
Frequency domain fluorescence lifetime data were collected using a Q2 laser scanner and A320 FastFLIM system (ISS, Inc., Champaign, IL, USA) mounted on a Nikon TE2000U microscope (Melville, NY, USA) and VistaVision software (ISS, Inc.). Acd or Atto 425 (the standard for calibration of the fluorescence lifetime) were excited using a 375 nm pulse diode laser (ISS, Inc.), driven by FastFLIM at the repetition rate of 10 MHz, with a 387 nm long-pass dichroic mirror, and emission was collected using a 451/106 nm band-pass emission filter and Hamamatsu model H7422P PMT detector. Affinity purified protein was used after about a 1:10 dilution in K+-Tris buffer. For each experiment, 11 μl of fluorescent sample was pipetted onto an ethanol-cleaned #1.5 glass coverslip mounted directly above the 10x 0.5 NA objective. Other reagents (maltose, [Ru(bpy)2phenM]2+, [Fe(phenM)3]2+, Fe2+, or EDTA) were pipetted directly into the sample drop and mixed at the final concentrations indicated in the text. For each condition, 256x256 confocal images were collected with a pinhole of 200 µm and a pixel dwell time of 1 ms. The pixels were averaged together for analysis, except as described for the phasor plot. The experimental phase delays (𝜑𝜔 ) and the modulation ratios (𝑚𝜔 ) of the fluorescence signal in response to an oscillatory stimulus with frequency 𝜔 were obtained using VistaVision software from the sine and cosine Fourier transform of the phase histogram H(p), subject to the instrument response function (IRF) calibrated with 2 µM Atto 425 in water with a lifetime of 3.6 ns (2,16,17). These data were fit with the Gaussian model described in the Supporting Material to obtain distance distributions.
In silico labeling and distance distribution simulations Computational modeling of Acd and metal-phenM labels, as well as distance distribution predictions, were performed using chiLife (18) with the accessible-volume sampling method (19,20). Acd, and cysteine conjugates of [Ru(bpy)2phenM]2+, [Fe(phenM)3]2+, [Fe(phenM)]2+, and [Cu(phenM)]2+ were added as custom labels in chiLife. Briefly, starting label structures were constructed in PyMOL and energy minimized with the GFN force field (GFN-FF) in xTB (21). Custom labels were superimposed onto labeling sites of the target pdb structure, and mobile dihedral angles were uniformly sampled. Rotamers with internal clashes (< 2 Å) were discarded. External clashes were evaluated using a modified, repulsive-only Lennard-Jones potential and used to weight rotamers as previously described (19). The 5
lowest weighted rotamers cumulatively accounting for a fraction of 0.005 of the total rotamer weights were discarded. Sampling was terminated once 10,000 samples had been attempted, generating between 400 and 2,500 rotamers, depending on the specific label and protein site. To calculate a simulated distance distribution between two rotamer ensembles, a weighted histogram was made for pairwise distances between the spin or fluorescent centers of each pair of rotamers from the two ensembles. For Acd, the center coordinates were defined by the mean position of all atoms in the central acridone ring. For the metal-phenM labels, the center coordinates were on the transition metal ion. Histograms were then convolved with gaussian distributions with a 1 Å standard deviation, and the resulting distributions were normalized.
[Cu(phenM)]2+ spin-labeling and EPR sample preparation Jo EPR measurements ur na lP re -p ro of
Purified MBP-295C-211C and MBP-322C-278C (~ 50 μM) were desalted (G-25) into K+-Tris buffer (pH 7.4) and immediately reacted with 0.5 mM phenM, freshly prepared from DMSO stock as a 5 mM solution in K+-Tris buffer with 1 mM EDTA. The reaction was nutated at 4 °C for 1 h, desalted (G-25) into K+-Tris buffer (pH 7.4), and concentrated (5 kDa MWCO). phenM-labeled MBP solutions were then incubated with 1 mM CuSO4 for 10 minutes at room temperature and loaded into 10 kDa MWCO microdialysis units (Thermo) and dialyzed against K+-Tris buffer prepared in deuterium oxide (D2O). Dialysis was carried out for ~ 18 h, replacing dialysis buffer with fresh deuterated K+-Tris buffer twice. RIDME samples were prepared with ~ 10 μM labeled MBP supplemented with 30 % (v/v) d8-glycerol. Holo MBP samples were additionally supplemented with 5 mM maltose from a stock solution in D2O. Samples were loaded into 1.5 mm OD/1.1 mm ID quartz tubes (Sutter) with flame-sealed bottoms and flash frozen in liquid nitrogen (LN2). Samples were stored at −80 °C until measurement. Samples for CW EPR were prepared similarly, but without use of deuterated buffers and with 25 % (v/v) glycerol. CW EPR samples were loaded into 4 mm OD quartz EPR tubes (Wilmad), frozen in LN2, and measured on the same day.
Continuous-wave EPR spectra were recorded at 112 K on a Bruker EMX spectrometer operating at Xband frequency (∼9.3 GHz) with a Bruker ER 4102SHQE resonator. Spectra were recorded with 100 kHz field modulation with a sweep rate of 3.6 G/s and a modulation amplitude of 5 G. Spectra were background subtracted and baseline corrected in LabVIEW™. For visual comparison between samples (Figure S4), spectra were normalized by the double integral of the respected field-modulated spectrum. Magnetic parameters g and A were determined by least-squares fitting spectra using EasySpin 6.0 (22), assuming axial g and A tensors and including anisotropic line broadenings as additional fit parameters (Figure S4). Pulse EPR experiments were performed at Q-band frequency (∼34 GHz) using a Bruker EleXsys E580 spectrometer with an overcoupled Bruker EN 5107D2 resonator. Pulses were amplified with a 300 W TWT amplifier (Applied Systems Engineering) and sample temperatures of 20 K or 10 K were maintained using a variable-temperature cryogen-free system (Bruker/ColdEdge). RIDME was performed using the established 5-pulse sequence (𝜋/2) ― 𝜏1 ― (𝜋) ― 𝜏1 + t ― (𝜋/2) ― TR ― (𝜋/2) ― (𝜏2 – t) ― (𝜋) ― 𝜏2 ― (echo) (23). 𝜋 /2 and 𝜋 pulses were 12 and 24 ns, respectively, and were applied at frequency and magnetic field values corresponding to the maximum of the Cu2+ echo-detected field swept spectrum. To avoid dynamic decoupling artifacts in the RIDME time-trace, 𝜏1 was chosen to be longer than 𝜏2, with 6
values of 4 μs and 3.5 μs for 𝜏1 and 𝜏2, respectively (24). Echo modulations due to solvent deuterium were suppressed by averaging over 𝜏1 and 𝜏2 with 16 ns increments over 8 steps. The relaxation interval TR was 195 μs and was selected to be ~ 0.75 that of the spin-lattice relaxation time (T1e) of the Cu2+ spin label at 20 K determined by inversion recovery experiments. Echo crossings were removed with a 32step phase cycle (25); however, a small echo crossing artifact at t ≈ 0 could not be removed by phase cycling. This artifact, along with residual ESEEM contributions, were removed by recording a second RIDME time-trace with identical pulse lengths and delays at 10 K, where TR ≈ 0.04T1e. Division of the 20 K data sets by the 10 K data sets gave the artifact-free dipolar evolution time traces used in all distance distribution analyses (Figure S5).
Divided and phase-corrected RIDME time-traces were analyzed with DeerLab version 1.1 (26) using Tikhonov regularization and compactness regularization (27). Residual intermolecular background was simultaneously modeled with the RIDME foreground using a homogenous 3-dimensional background model. The dipolar kernel was modified by replacing the default free electron g-value with an effective g-value of 2.1203 for each Cu2+ ion, which was determined from fits to the dual labeled MBP[Cu(phenM)]2+ CW EPR spectra. Uncertainty estimations were determined by the asymptotic method in DeerLab and 95% confidence intervals were plotted as error bands on the RIDME probability distributions. All plots were generated in KaleidaGraph version 5.0 (Synergy Software) and visualized with KaleidaGraph and Inkscape version 1.2.
Time-resolved FRET can be measured in either the time domain, usually with time correlated singlephoton counting (TCSPC), or in the frequency domain, with both approaches yielding equivalent information (2). Here, we measured fluorescence lifetimes using a frequency-domain lifetime instrument (see Materials and Methods). For frequency-domain measurements, the lifetime data are visualized with a Weber plot that shows the phase shift of the response (phase delay) and the decrease in the amplitude of the response (modulation ratio) as a function of the modulation frequency of the excitation light (Figure 2A). These data can then be fit with models for the lifetimes that assume that they have single-exponential, multiexponential, or nonexponential decays. To investigate the utility of time-resolved FRET to measure distance distributions, we incorporated the fluorescent noncanonical amino acid acridon-2-ylalanine (Acd) into MBP using amber codon suppression in bacteria as previously described (13,14,28). MBP is a clamshell-shaped protein that undergoes a significant closure of the clamshell upon binding maltose. For these experiments, we used two donor sites for specific incorporation of Acd, amino acid 295 at the outer lip of the clamshell and 322 on the back side of the clamshell. These donor fluorophore sites were then paired with single cysteine mutations for incorporation of transition metal chelates as FRET acceptors. Using MBP allowed us to test if time-resolved tmFRET could measure distances, distance distributions, and state energetics over a range of distances in a protein with a well-characterized structure and conformational rearrangement. Acd incorporated at both donor sites (MBP-322Acd and MBP-295Acd) exhibited long, single-exponential fluorescence lifetimes similar to free Acd (Figure 2A and 5A, gray symbols). For Acd at 295 in wildtype MBP, the lifetime was nonexponential due to quenching by proximal Y307; therefore, all of our MBP295Acd constructs also contain a Y307S mutation (referred to here as MBP-295Acd) (12,29). For both sites, the lifetimes were slightly longer in the presence of maltose (MBP-295Acd: apo, 14.7 ± 0.02 ns 7
(n=13); holo 15.3 ± 0.1 ns (n=6); MBP-322Acd: apo, 15.4 ± 0.01 ns (n=20); holo 15.7 ± 0.01 ns (n=9)). This is likely due to a small change in the environment of the incorporated Acd in the presence of maltose and was factored into our subsequent analysis.
Time-resolved FRET with [Ru(bpy)2phenM]2+ and [Fe(phenM)3]2+produce accurate average distances and narrow distance distributions ro of
In the previous paper, we have shown that the cysteine-reactive metal chelate [Ru(2,2′-bpy)2(1,10phenanthroline-5-maleimide)]2+ ([Ru(bpy)2phenM]2+) can act as a long-distance tmFRET acceptor for Acd (13). [Ru(bpy)2phenM]2+ exhibits a substantial absorption in the visible range that overlaps with the emission spectrum of Acd, giving an R0, the distance producing 50% FRET efficiency, of 43.5 Å. Consistent with its R0, labeling of MBP-322Acd-278C with [Ru(bpy)2phenM]2+ produced a substantial maltose-dependent decrease in steady-state Acd fluorescence indicating the presence of FRET between Acd and [Ru(bpy)2phenM]2+ (13). These steady-state FRET measurements, however, do not reveal the conformational heterogeneity in the sample.
To determine if time-resolved FRET could be used to measure distance distributions, we measured the fluorescence lifetimes of MBP-322Acd-278C labelled with [Ru(bpy)2phenM]2+. As shown in Figure 2A, labeling with [Ru(bpy)2phenM]2+ caused a substantial decrease in the average lifetime (manifesting as a shift in the phase delay and modulation ratio curves to higher frequencies in the Weber plot). Furthermore, subsequent addition of maltose increased the average lifetime (shifting the phase delay and modulation ratio curves to lower frequencies) reflecting an increase in the average distance between Acd at 322 and [Ru(bpy)2phenM]2+ at 278C, as predicted from structural modeling (see below). No change in lifetime was observed for MBP-322Acd without the cysteine mutation (data not shown; no native cysteines are present in MBP). These data suggest that time-resolved tmFRET could be used to reveal the intramolecular distance distributions in our samples.
To quantify the distance distributions from the fluorescence lifetime data, we fit the data with a model that predicts the lifetimes for a distribution of distances. The model assumes the following: 1) The fluorescence lifetimes of the donor-only protein (i.e., in the absence of acceptor) is single-exponentially distributed with a time constant 𝜏𝐷 (though 𝜏𝐷 can be different in the resting and active states). 2) The decrease in lifetime in the presence of the attached acceptor is due to a FRET mechanism with a known R0. 3) Donor and acceptor dipoles are randomly oriented relative to each other (κ2 = 2/3), a reasonable assumption when one member of the FRET pair is a metal ion (30). 4) There is only a single acceptor for each donor. 5) The distance distribution for each state can be approximated by a Gaussian distribution with a distinct mean distance and standard deviation. And 6) the distances do not change appreciably on the time scale of the fluorescence lifetime. Most of these assumptions can be experimentally verified in our sample. This Gaussian model was globally fit to the phase delay and modulation ratio data across multiple conditions in the same experiment (apo, holo, and intermediate concentrations of ligand). The values of 10 to 12 free parameters for each experiment were determined using 𝜒 2 minimization. Because of how the distributions were parameterized, the fits generally had fewer free parameters than fitting with a sum of exponentials. This approach was pioneered in the 1970s, primarily by Steinberg and coworkers (31,32). Global fits of the Gaussian model to the fluorescence lifetime data of MBP-322Acd-278C labelled with [Ru(bpy)2phenM]2+ in the apo and holo states are shown in Figure 2A. The values of the mean distances 8
and standard deviations in the apo and holo states for 14 different experiments are shown in a spaghetti plot in Figure2B and a scatter plot in Figure 2C and were very reproducible. To compare our data to the predictions of structural modeling, rotameric ensembles of Acd and [Ru(bpy)2phenM]2+ were modelled onto crystal structures of apo and holo MBP (33,34) using the accessible-volume approach in chiLife (18), and used to predict distance distributions for MBP-322Acd-278C labeled with [Ru(bpy)2phenM]2+ (Figure 2B). In both the absence and presence of maltose, the distance distributions were remarkably similar to those predicted by chiLife. The average experimental distances and maltose-dependent delta distance were within 1 Å of those predicted by chiLife, and the widths of the distributions were also similar (though the experimentally determined width was consistently larger for the holo state). This similarity suggests that specific interactions between the probes and the protein (not accounted for by chiLife) do not play a significant role at these sites. Both the experimentally determined and predicted distance distributions between Acd and [Ru(bpy)2phenM]2+ were surprisingly narrow, indicating that [Ru(bpy)2phenM]2+ did not add appreciably to the measured heterogeneity under these conditions.
To determine the ability of tmFRET between Acd and [Ru(bpy)2phenM]2+ to resolve probability distributions that are a mixture of resting and active states, we performed lifetime experiments with subsaturating maltose concentrations applied to MBP-322Acd-278C labelled with [Ru(bpy)2phenM]2+, one just below the KD and another just above the KD. The intermediate concentrations produced intermediate curves on the Weber plot (Figure 3A). To analyze the lifetimes in a model-independent way, we graphed the data on a phasor plot, a plot of the out of-phase versus the in-phase components of the fluorescence in frequency domain experiments (16). This plot revealed that the data for these intermediate concentrations fall on a line between the zero and saturating maltose concentrations, indicating that these intermediate maltose concentrations are a mixture of the same resting and active states produced by zero and saturating maltose (Figure 3B). We therefore performed global fits of the Gaussian model to the data on the Weber plot at all four maltose concentrations (0, 5 μM, 9.2 μM, and 3 mM), constraining the mean distance and standard deviation for each state to be the same for all conditions and allowing the fraction in the active state to vary (Fig 3C). The maltose dependence of the fraction of the active state nicely conformed to a binding isotherm, yielding an affinity of 6.3 μM, consistent with that previously measured using steady-state fluorescence (13) (Figure 3D). This is remarkable given that the difference in the average distance between the resting and active states is only about 6 Å. We also performed similar time-resolved FRET experiments on MBP-322Acd-278C with [Fe(phen maleimide)3]2+ ([Fe(phenM)3]2+). Like Ru2+, Fe2+ forms complexes with three bipyridyls or phenanthrolines that are highly absorbent in the visible range and can be used as FRET acceptors with visible fluorophores (13). The R0 of Acd with [Fe(phenM)3]2 is 41.8 Å. The addition of [Fe(phenM)3]2+ caused a dramatic decrease in the average fluorescence lifetime of MBP-322Acd-278C as expected for FRET between Acd and [Fe(phenM)3]2+ (Figure 4A). Furthermore, the average lifetime systematically increased upon addition of increasing concentrations of maltose, consistent with the maltosedependent increase in distance predicted for MBP-322Acd-278C. No change in lifetime was observed for MBP-322Acd without the cysteine mutation (data not shown). Global fits of the Gaussian model to the data on the Weber plot at five different maltose concentrations (0, 5 μM, 9.2 μM, 12.9 μM, and 3 mM) yielded mean distances and standard deviations similar to those for [Ru(bpy)2phenM]2+, even though the R0 was somewhat smaller (Figure 4D). The phasor plot (Figure 4B) and average distributions (Figure 4C) display a gradual shift in the equilibrium between the resting and active state with increasing maltose 9
concentrations. Finally, the maltose dependence of the fraction of the active state nicely conformed to a binding isotherm, yielding an affinity of 9.2 μM for maltose (Figure 4E). These experiments demonstrate that both [Ru(bpy)2phenM]2+ and [Fe(phenM)3]2+ make good acceptors when using time-resolved tmFRET to determine distance distributions among conformational states and the free energy difference between conformational states.
Time-resolved FRET with [Fe(phenM)]2+ ro of
For MBP-295Acd-211C, the predicted donor-acceptor distances are shorter (on the order of 30 Å) and, instead of maltose increasing the donor-acceptor distance like for MBP-322Acd-278C, maltose decreases the distance. At these distances, the FRET efficiency with [Ru(bpy)2phenM]2+ and [Fe(phenM)3]2+ is nearly one, making these acceptors inadequate for detecting distance changes in MBP-295Acd-211C (13). However, we have found that Fe2+ bound to a single phenanthroline maleimide ([Fe(phenM)]2+) exhibits a much lower absorption and therefore a much lower R0 (24.4 Å), making it ideal for the donoracceptor distances in MBP-295Acd-211C (13).
The addition of Fe2+ caused a dramatic decrease in the average fluorescence lifetime of MBP-295Acd211C-phenM as expected for FRET between Acd and [Fe(phenM)]2+ (Figure 5A). The decrease in lifetime was even greater in the presence of saturating concentrations of maltose, as expected from the maltose-dependent decrease in distance between the donor and acceptor in MBP-295Acd-211C. For [Fe(phenM)]2+, but not for [Ru(bpy)2phenM]2+ or [Fe(phenM)3]2+, the FRET was fully reversible with EDTA, as seen by the open symbols in the Weber plot surrounding the original donor-only points (Figure 5A). Fitting the Gaussian model to the lifetime data revealed that MBP-295Acd-211C-phenM had a somewhat broader distance distribution than MBP-322Acd-278C, but with a maltose-dependent decrease in average distance (Figure 5B and C). Once again, the distance distributions were similar to those predicted by chiLife, albeit with a somewhat shorter distance than predicted for the apo state. These data indicate that [Fe(phenM)]2+ can be used as an effective FRET acceptor for lifetime measurements of distance distributions with mid-range donor-acceptor distances.
Pulse dipolar EPR with [Cu(phenM)]2+ To further explore the heterogeneity contributed by the phenM side chain, we turned to pulse dipolar EPR spectroscopy. Pulse dipolar EPR methods such as double electron-electron resonance (DEER) and relaxation-induced dipolar modulation enhancement (RIDME) measure distance distributions between unpaired electrons introduced into proteins via site-directed spin labeling (35). These data are typically analyzed using non-parametric models and therefore place no underlying assumptions on the shape of the distance distributions, only that they are smooth. For these experiments, we introduced cysteine substitutions at both donor and acceptor sites in our previous MBP constructs, generating MBP-295C211C and MBP-322C-278C. Nitroxide spin labels introduced at these site pairs have previously been employed to detect maltose-dependent conformational changes in MBP using DEER (15). To examine distance distributions obtained using the phenM label, MBP-295C-211C and MBP-322C-278C were labelled with phenM and subsequently labelled with Cu2+ (Figure 6A). Cu2+ is a d9 transition metal ion containing one unpaired electron (S = 1/2) and can therefore be used as a spin label for EPR experiments. Indeed, both MBP-295C-211C and MBP-322C-278C dual labeled with [Cu(phenM)]2+ displayed continuous-wave EPR spectra consistent with each Cu2+ coordinated by a single phenanthroline, indicating specific labeling of the phenM side chains with Cu2+ (Figure S4). 10
To determine distance distributions between [Cu(phenM)]2+ labels on MBP directly, we performed RIDME experiments in the absence and presence of saturating maltose (Figure 6B, S5). For MBP-295C211C labeled with [Cu(phenM)]2+, RIDME in the absence of maltose reveals a broad distance distribution centered at 35.1 Å, which narrowed and shifted to a most probable distance of 28.6 Å in saturating maltose (Figure 6B). These distances are in excellent agreement with those predicted by chiLife (36.1 Å and 28.4 Å for apo and holo, respectively) (Figure 6B dashed line). RIDME data on MBP-322C-278C labeled with [Cu(phenM)]2+, positioned on the backside of the clamshell, reveal the expected increase in Cu2+‒Cu2+ distance upon addition of maltose, with distance distributions centered at 43.5 Å and 49.5 Å for apo and holo conditions, respectively (Figure 6C). Again, these distances are within 1 Å of the chiLife predictions (43.1 Å and 49.8 Å for apo and holo, respectively). These results establish [Cu-phenM]2+ as a useful and commercially available spin label for determining conformational distributions in proteins using pulse dipolar EPR spectroscopy. Moreover, the similarity of both the RIDME distributions (Figure 6B,C) and the tmFRET distributions (Figure 2B and 5B) to their respective chiLife predictions, and to each other, suggests the distributions at room temperature are well captured by the rapid freezing of the sample in the RIDME experiments. Overall, these experiments support the accuracy of RIDME, tmFRET, and chiLife for estimating the rotameric ensembles and resulting distance distributions involving [metal(phenM)]2+ labels.
This paper applies time-resolved tmFRET to study protein allostery and conformational dynamics. tmFRET utilizes a fluorescent noncanonical amino acid as the donor and our new metal-bipyridyl derivatives as the acceptor to overcome limitations of traditional FRET methods (13). We applied this method to MBP and demonstrated it can accurately determine distances, conformational heterogeneity, and energetics. The results highlight the utility of time-resolved tmFRET in characterizing protein dynamics and conformational changes, offering valuable insights into the mechanisms of allosteric regulation.
In addition to its utility in tmFRET, we showed that the cysteine-reactive bipyridyl derivative phenM can be used with Cu2+ as a spin label in pulse dipolar EPR spectroscopy. Distance distributions from pulsed dipolar EPR are commonly determined using non-parametric models which have the advantage of requiring no underlying assumptions about the shape of the distance distribution. Using the same sites on MBP to which we introduced the donor and acceptor for tmFRET, we show that the distance distributions produced by RIDME were similar to the predictions of chiLife and also to the distance distributions measured with time-resolved tmFRET. These results both establish [Cu(phenM)]2+ as a spin label for pulse dipolar EPR spectroscopy and validate the distance distributions measured with timeresolved tmFRET. The determination of distance distributions directly from the lifetime data is an ill-posed problem and, therefore, we needed to parameterize the distance distributions for our time-resolved tmFRET experiments. For our model, the heterogeneous distances between the donor and acceptor for each state were described by a Gaussian distribution with a distinct mean distance and standard deviation, although other distance distributions such as Lorentzians have also been used (36-38). Given these assumptions, values for the parameters were well determined from an analysis of the 𝜒 2 surfaces (Figure S2A-D) and were fairly consistent across multiple experiments (Figures 2C, 4D, 5C). In addition, the distances, conformational heterogeneity, and energetics determined from the Gaussian model were 11
consistent with molecular modeling (chiLife, Figure 2B, 5B) and EPR spectroscopy (RIDME, Figure 6), and subsaturating maltose concentrations (Figure 3 and 4). lP