Anna
Nalepa
a,
Marco
Malferrari
b,
Wolfgang
Lubitz
a,
Giovanni
Venturoli
bc,
Klaus
Möbius
ad and
Anton
Savitsky
*a
aMax-Planck-Institut für Chemische Energiekonversion, Stiftstrasse 34-36, D-45470 Mülheim an der Ruhr, Germany. E-mail: anton.savitsky@cec.mpg.de; Fax: +49-208-3063951; Tel: +49-208-3063555
bLaboratorio di Biochimica e Biofisica, Dipartimento di Farmacia e Biotecnologie, FaBiT, Università di Bologna, via Irnerio 42, I-40126 Bologna, Italy
cConsorzio Nazionale Interuniversitario per le Scienze Fisiche della Materia (CNISM), c/o Dipartimento di Fisica, Università di Bologna, via Irnerio 46, I-40126 Bologna, Italy
dDepartment of Physics, Free University Berlin, Arnimallee 14, D-14195 Berlin, Germany
First published on 25th September 2017
Using isotope labeled water (D2O and H217O) and pulsed W-band (94 GHz) high-field multiresonance EPR spectroscopies, such as ELDOR-detected NMR and ENDOR, the biologically important question of detection and quantification of local water in proteins is addressed. A bacterial reaction center (bRC) from Rhodobacter sphaeroides R26 embedded into a trehalose glass matrix is used as a model system. The bRC hosts the two native radical cofactor ions (primary electron donor) and
(primary electron acceptor) as well as an artificial nitroxide spin label site-specifically attached to the surface of the H-protein domain. The three paramagnetic reporter groups have distinctly different local environments. They serve as local probes to detect water molecules via magnetic interactions (electron–nuclear hyperfine and quadrupole) with either deuterons or 17O nuclei. bRCs were equilibrated in an atmosphere of different relative humidities allowing us to control precisely the hydration levels of the protein. We show that by using oxygen-17 labeled water quantitative conclusions can be made in contrast to using D2O which suffers from proton–deuterium exchange processes in the protein. From the experiments we also conclude that dry trehalose operates as an anhydrobiotic protein stabilizer in line with the “anchorage hypothesis” of bio-protection. It predicts selective changes in the first solvation shell of the protein upon trehalose–matrix dehydration with subsequent changes in the hydrogen-bonding network. Changes in hydrogen-bonding patterns usually have an impact on the global function of a biological system.
In biological systems water is usually divided into three distinctly different classes: (i) internal water molecules that are hydrogen-bonded to specific amino-acid residues in a water pocket or are mobile along inner-protein channels, they are often of key importance for protein function;5–7 (ii) surface water molecules in the hydration shell of the protein at the solute–solvent interface, and (iii) bulk water molecules randomly distributed in the protein matrix. Water molecules in the protein hydration layer have restricted dynamics with respect to water molecules in the bulk as shown using different experimental methods which probe different dynamic ranges of molecular motion, i.e. Nuclear Magnetic Resonance (NMR),8–11 Dielectric Spectroscopy in the microwave region,12 Terahertz Absorption Spectroscopy,13,14 Infrared Spectroscopy (IR),15,16 and Neutron Spectroscopy.1 Additionally, powerful extensions of molecular dynamics simulation techniques have been recently employed to unravel solvation dynamics of large molecules.13,17–19 At physiological temperatures solvated proteins fluctuate between a multitude of conformational sub-states in a rugged energy landscape that is hierarchically organized in energy tiers.20,21 For most proteins cooled below the glass-transition temperature (typically around 200 K) their biological function is blocked due to restricted conformational motion. The “freezing out” of conformational dynamics by lowering the temperature is thus a common strategy for studying function-dynamics relationships in proteins.
However, freezing proteins in the presence of a cryoprotectant (used to minimize freezing damage by ice crystals) is problematic because it aggravates disentanglement of the influences of solvent and temperature on the protein dynamics. An elegant alternative approach is to embed the protein into amorphous matrices formed by disaccharides like trehalose.22,23 This allows preservation of the native protein fold during extensive protein dehydration, even at temperatures well above room temperature. In nature, the extraordinary bio-protective capabilities of disaccharide glasses are exploited by specific organisms, which are able to survive extreme conditions of temperature and dehydration by entering a state of reversibly arrested metabolic activity, called anhydrobiosis or cryptobiosis.24 Extensive spectroscopic work, exploiting neutron scattering,25,26 Raman,27 optical laser-flash,28 FTIR,29 and EPR (Electron Paramagnetic Resonance),30 on different proteins incorporated into trehalose glasses revealed a tight protein–matrix dynamic coupling at a low water content, implying that the protein dynamics is controlled by that of the water–trehalose matrix coating the protein surface.25–31 At room temperature, the stepwise dehydration of the trehalose matrix results in increasingly inhibited dynamics of the embedded protein. This was observed in both small globular proteins like myoglobin22,25 and large membrane proteins like bacterial photosynthetic reaction centers23,32 and photosystem I33 of oxygenic photosynthesis.
Trehalose is the most efficient sugar for bio-protection against extreme dehydration and osmotic stress. The exact mechanistic details of the trehalose efficiency are not clear yet (for a recent review see ref. 34), but likely involve an interplay between several factors that include its extraordinarily high glass-transition temperature (385 K),35 its polymorphism adopting several crystalline and amorphous states,36 its distinct propensity for hydrogen-bonding and the pronounced rigidity of its dehydrated glass matrix.37
In this work, we apply pulsed multi-resonance EPR techniques, such as the double-resonance extensions ENDOR (electron–nuclear double resonance) and EDNMR (ELDOR-detected NMR), at a high microwave (mw) frequency (W-band, 94 GHz) and a correspondingly high magnetic Zeeman field (3.4 T). Both techniques are extremely powerful for measuring small electron–nuclear hyperfine couplings even in large low-symmetry biomolecules. They allow us to quantify the accessibility of local water at various sites in bacterial photosynthetic reaction centers (bRCs) from Rhodobacter (Rb.) sphaeroides R26 that were embedded in trehalose glasses at distinct hydration levels. The advantage of EPR spectroscopy in comparison with other techniques such as THz spectroscopy and IR spectroscopy is that site-specific information on local water of the solvation layer can be obtained with a minimum of assumptions in the spectral analysis. Employing isotope labeled water provides the necessary contrast from internal non-exchangeable water molecules.
A bRC is an integral membrane protein that catalyzes the initial photochemical processes to convert light energy into chemical free energy.38–40 The three protein subunits, L, M, and H, host several cofactors41 sequentially involved in light-driven electron transfer. The two bacteriochlorophyll a (BChl a) molecules near the periplasmic side of the membrane form the “special pair” P865. After photo-excitation to its singlet state, the subsequent electron transfer proceeds via the protein branch A38 to the QA acceptor (see Fig. 1). In the present study, the secondary electron transfer is blocked by using the inhibitor stigmatellin.42
![]() | ||
Fig. 1 (a) X-ray structural model of the bRC from Rb. sphaeroides R2643 with the protein subunits embedding the electron-transfer cofactors P865 (the “special pair” donor), two BChls, two BPhes (bacteriopheophytin a), two ubiquinone-10 acceptors (QA and QB) and a non-heme iron Fe2+. The paramagnetic Fe2+ cofactor can be replaced by diamagnetic Zn2+. QA and QB form the primary and secondary electron acceptor, respectively. Here, the quinone QB is replaced by stigmatellin to block electron-transfer from ![]() |
Using EPR spectroscopy one can measure the electron Zeeman and electron–nuclear hyperfine interactions and thereby probe the unpaired electron spin density distribution in the transient ionic radical states, and
, that are created during light-induced charge separation. These paramagnetic states have been thoroughly characterized previously using both cw (continuous wave) and pulsed EPR, for review see ref. 45 and 46. P865 and QA are located close to the periplasmic and cytoplasmic sides of the membrane, respectively, and are thus well suited to be used to study local water accessibility in the protein near its surface. In the bRC preparation of this work, with no cytochrome present and electron transfer from
to QB blocked, the charge-separated state
recombines to P865QA within 100 ms by direct electron tunneling.38 Hence, continuous illumination generates a steady-state EPR signal of the
state that rapidly decays once the light is switched off, see Fig. S1 in the ESI.† The W-band EPR spectra of both
and
reveal rhombic g-tensors47,48 resulting in characteristic line shapes, see Fig. 2a.
![]() | ||
Fig. 2 (a) Top: Simulated rigid limit W-band EPR absorption spectra of ![]() ![]() |
The native bRC contains five cysteine residues which are buried within the protein domains except for cysteine 156 which is moderately solvent exposed in subunit H. This residue can be site-specifically spin labeled with an external paramagnetic probe molecule. We used MTSSL (1-oxyl-2,2,5,5-tetramethylpyrroline-3-methyl-methanethiosulphonate) as the spin label.49,50 The nitroxide spin label51 provides a third paramagnetic probe within the bRC in addition to the “natural” radical ions and
. Due to the protein 3D structure the “artificial” nitroxide spin label experienced a water accessibility very different from that of the native cofactors, P865 and QA, see Fig. 1b.
Even under high-field (W-band) conditions, the close similarity of the g-tensor components results in strong overlap of the EPR spectra of ,
and nitroxide spin labels embedded in their protein matrix environment. For the present study, this dictates the necessity of two variants of bRC protein preparations to be investigated independently: (i) the nitroxide-labeled bRC, (ii) the bRC with the Fe2+ cofactor replaced by diamagnetic Zn2+ to enable resolved EPR detection of both native cofactors,
and
, without stranding because of excessive line broadening by paramagnetic Fe2+.53
The effect of bRC coatings on a trehalose glass matrix has been thoroughly studied previously using laser-flash optical and EPR spectroscopies.23,28,30,32,54 In bRCs, the coating with trehalose affects the electron-transfer step between the primary and secondary quinone acceptors32 which is a conformationally gated process.55 The trehalose matrix also influences the lifetime of the charge-separated radical-pair state, . Upon progressive dehydration of the trehalose matrix, the kinetics of
charge recombination becomes faster and exhibits widely distributed rate constants.23,28,54,56 Thus, it mimics at room temperature the bRC recombination kinetics observed at cryogenic temperatures in the water–glycerol system when frozen in the dark.57,58
A reasonable explanation for this matrix effect is the stronger dynamic coupling between sugar and protein via hydrogen-bonding networks in a dehydrated trehalose matrix. It blocks the fluctuations between conformational substates of the bRC,28 as put forward in the “anchorage hypothesis”.31 Two relative hydration levels, r, were chosen in this work that result in a different extent of slaving the protein dynamics to the embedding sugar matrix: (i) r = 11% for which the protein dynamics is arrested on the time scale of seconds,54 and (ii) r = 74% for which the dynamics is only mildly retarded as compared to solution.59
Using EDNMR we aim to probe directly if the dehydration of the trehalose matrix results in a changed hydration level of the inner protein core surrounding the native cofactors, and
, and of the protein surface at the labeled site of cysteine 156. To trace the exchange of water between the sugar matrix and the protein, the use of isotope labeled water is essential, with regard to both deuterium and oxygen-17.
Our recent studies59,60 showed that the molecular organization of the trehalose matrix is such that in the dehydrated state it forms a homogenous amorphous matrix. Therefore, in the present work the samples were always prepared by dehydration first before the process of rehydration was studied.
The present paper is organized as follows: before reporting on the ultimate subject of this study, the water exchange in bRCs embedded in the trehalose matrix of different water (D2O) contents, we have to characterize the molecular organization of the trehalose matrix without proteins at different hydration levels. Only with this information at hand meaningful conclusions can be drawn on prevailing trehalose–protein interaction mechanisms as mediated by surrounding water molecules. Hence, (i) we start with a careful EPR characterization of the trehalose matrix prepared at relevant relative D2O humidity levels in the absence of bRCs. (ii) Subsequently, the influence of bRC trehalose coatings at different regions of the protein complex on the water accessibility of the distinctly different binding sites of the three spin probes, the surface-localized nitroxide radical as well as the inner-protein localized native cofactor radical ions and
, is compared for the relative D2O humidity levels, r = 11% and r = 74%. To this end, the rehydration timescale of the trehalose film is determined at these three binding sites by comparing the deuterium resonance signal intensity in the respective W-band EDNMR recordings. We end the present paper with the H217O water-accessibility results obtained by W-band EDNMR and Davies ENDOR on nitroxide spin labeled bRC/trehalose samples in which 17O is used as a nuclear spin probe instead of 2H. The respective advantages of both approaches are highlighted.
Both Fe2+- and Zn2+-containing bRC preparations were concentrated to 335 μM and supplemented with stigmatellin, a potent inhibitor of the -to-QB electron transfer.42 By removal of the iron in Zn-bRCs and replacement of QB by stigmatellin fast spin relaxation of the
radical anion by high-spin Fe2+ (S = 2) is avoided as well as unwanted formation of the semiquinone biradical
in a subsequent electron-transfer step occurring in wild-type bRCs. Hence, the charge-separated primary donor–acceptor radical pair,
, is exclusively formed upon illumination and its transient lifetime and spin dynamics are not shortened by spin interactions with nearby Fe2+ and consecutive secondary electron transfer.
The bRC/trehalose glassy samples for EPR measurements were prepared as described previously.30 An amount of 25 μL of the concentrated (335 μM) bRC solution was mixed with 25 μL of 1.67 M trehalose dissolved in 10 mM Tris/HCl buffer containing 0.025% LDAO. This results in a sugar/bRC molar ratio of 5 × 103. The bRC/trehalose solution was layered on an optical glass and dried in a desiccator under N2 flow to form a film at room temperature. The desired glassy sample was obtained after 4 h of gas flux. Subsequently, the samples were transferred into a sealed box where further dehydration of the sample occurred via controlled exposition to an atmosphere of low relative humidity over a saturated aqueous salt solution at room temperature using H2O with natural isotope abundance. This isopiestic method of controlled sample dehydration has been employed successfully both in EPR30 and IR65–67 studies on bRC/trehalose glasses and bRC films dehydrated in the absence of sugar. For measurements on samples equilibrated to the desired relative humidity with either D2O or H217O vapor, the saturated solutions of an adequate salt68 were prepared using the respective isotope-labeled water at room temperature. The slight change in vapor pressure due to the use of D2O instead of H2O is negligible for our purpose.69 The dehydrated bRC/trehalose films containing about 0.5 mM of protein were then crumbled into small flakes and inserted into the W-band sample quartz capillaries (i.d. 0.6 mm). In the present study two specific hydration levels of bRC/trehalose glasses were obtained when using LiCl (r = 11%) or NaCl (r = 74%) saturated solutions.68 Throughout the text, the type of water vapor used in the equilibration procedure is indicated in brackets following the value of r, e.g., 11% (LiCl/D2O) corresponds to equilibration at r = 11% using a saturated LiCl solution in D2O.
The R1 nitroxide/trehalose glasses were prepared in an analogous way, except for using a higher nitroxide concentration (≈1 mM) and omitting the detergent in the starting trehalose solution. Two independent preparations were used for two EPR experiments in order to prove the consistency of the obtained data.
The EDNMR experiments employed two-pulse Hahn-echo detection at mw frequency ν1 that is matched to the cavity resonance frequency. The required high-turning angle (HTA) mw pulse was applied at frequency ν2. The EDNMR spectra were acquired by continuously sweeping ν2 at a fixed magnetic field, usually in steps of 70 kHz. The ν2 frequency-sweep range was adjusted individually, depending on the resonance signal of interest. The width of the central “blind spot”, which occurs when ν2 approaches ν1, was used to measure the mw field amplitude, ω1, at frequency ν1 of the detection pulse sequence, as described in ref. 70. The HTA pulse length and mw field amplitude were optimized for the respective nucleus of interest. For details of the mw pulse scheme settings, see the ESI.† Similar to ENDOR, in the EDNMR spectra the hyperfine or/and quadrupole coupled proton or deuterium nuclei yield signals centered at the nuclear Zeeman frequency νn(H) or νn(D), split by the hyperfine or quadrupole coupling constant.
EDNMR spectra were processed as described previously.70 In short, the spectra were inverted and normalized with respect to the central blind-spot depth. The blind spot was subsequently fitted by using either a Lorentzian line shape or a cubic spline and subtracted (see Fig. S2 in the ESI†). The resulting EDNMR signals have intensities between 0 and 1. Since the information content for weakly coupled D and 17O signals is equivalent for left-hand (ν2 < ν1) and right-hand (ν2 > ν1) sides of the EDNMR spectra,71 only signals with Δν = (ν2 − ν1) > 0 are shown.
ENDOR spectra were recorded using the Mims ENDOR pulse sequence, tp–τ–tp–trf–T–tp–τ–echo, for deuterium or the Davies ENDOR pulse sequence, tinv–trf–T–tp–τ–2tp–τ–echo, for 17O nuclei. A τ value of 320 ns was chosen for Mims ENDOR. The length of the selective π pulse, tinv, in Davies ENDOR was set to 300 ns. The radio frequency pulse lengths, trf, of 35 μs and 23 μs were used for 2H and 17O ENDOR, respectively.
The intensity of the D line in EDNMR recordings continuously increases during the first 40 hours of equilibration until reaching a plateau after 3 days with a D-line intensity of 0.07 ± 0.01, see Fig. 3c and Table S1 in the ESI.† The deuterium signal originates mainly from D2O, but presumably also from exchangeable hydroxyl groups within the trehalose molecule73,74 which would add intensity to the deuterium line in the EDNMR spectrum. Therefore, we included, as a reference, the nitroxide/trehalose glass prepared under identical conditions but using D2O as the solvent. Such a preparation procedure allows for fast exchange of hydroxyl groups in the trehalose molecules. In this sample the intensity of the deuterium signal in the EDNMR spectrum is two times larger, see Fig. 3a. This suggests that in the dehydrated sugar glass the exchange of hydroxyl groups in trehalose molecules is slower compared to liquid solution. Thus, the exchanged hydroxyl groups do not contribute significantly to the deuterium signal intensity in the EDNMR spectrum on the equilibration timescale employed. The observed signal can be ascribed to D2O molecules, proving that the exchange of residual water retained in the trehalose takes place at a constant low relative humidity level. The exchange process occurs in spite of the high rigidity of the trehalose glass at r = 11% (which corresponds to only 0.52 ± 0.03 water per trehalose molecule59). The equilibration process is, however, slower as compared to rehydration from r = 11% to r = 74%, see Fig. S4 in the ESI.† In the latter case the new hydration level is reached within 24 hours as compared to 48 hours at r = 11%.
From quantum-chemical calculations and EPR observations it is known that nitroxide radicals form H-bonds with nearby hydrogen donor groups.70,75–77 Hence, the deuterium line in the EDNMR spectrum is assigned to either the deuterated water molecules directly coordinated to nitroxide radicals and/or to D2O molecules beyond the first solvation shell. The coordinated water molecules yield an EDNMR signal centered at the nuclear Zeeman frequency νn(D) split by the hyperfine, A, and quadrupole coupling constants (nuclear spin of deuterium I = 1). The distant water molecules would exhibit small, unresolved hyperfine couplings that merge to a Gaussian line centered at νn(D), the so called “matrix line”. This means that in the weak-coupling regime (A < 2·|νn(D)|), the spectral contributions from both coordinated and distant water molecules are centered at νn(D) overlapping partially. However, in all our EDNMR spectra the signal at νn(D) could satisfactorily be fitted to a Gaussian line shape with a full width at half maximum (FWHM) of ca. 1.1 MHz (at a 1.2 μs echo integration window) – with no hint of additional structure or splitting.
To obtain additional information whether both coordinated and distant water molecules contribute, W-band Mims ENDOR spectra were recorded. This technique has proven to be useful for the detection of narrow spectral line shapes.78Fig. 3d shows Mims ENDOR spectra recorded after different equilibration times for the R1/trehalose glass at r = 11% (LiCl/D2O), analogous to the above EDNMR recordings. The ENDOR spectra were taken at the magnetic field position in the gxx–gyy EPR spectral range where the largest hyperfine coupling is expected in the case of planar H-bonded nitroxide complexes. After 6 hours of equilibration, both the distant and H-bonded deuterium signals are seen in the ENDOR spectrum (bottom trace in Fig. 3d). The hyperfine coupling at this field position is 1.8 ± 0.3 MHz, but the lines are relatively broad and the quadrupole splitting from deuterium is not resolved. The rigidity of the dry trehalose glass is likely to prevent structural relaxation to the energetically most favorable H-bond geometry, resulting in broader ENDOR lines. In the course of equilibration, the signal intensity of the matrix line increases as compared to hyperfine-coupled deuterium, especially between 6 and 27 hours of equilibration, but the overall spectral changes are negligible. The fact that no change in the ENDOR line shape occurs between 2, 3 or 7 days suggests that equilibration is complete. This is in line with the results from our EDNMR experiments.
For the trehalose glass prepared directly in D2O the ENDOR line shape does not change. The same H-bond donor group, i.e., the hydroxyl group of water and trehalose is likely to have a similar hyperfine coupling constant. The only difference is a slightly broader matrix line which is explained by the larger number of deuterium nuclei in the vicinity of the spin probe. These originate from exchanged hydroxyl protons in the trehalose molecule.
Overall, the lack of changes in the coupled deuterium signal strongly suggests that the structural organization of dry amorphous trehalose (r = 11%) does not change at the molecular level in the course of rehydration equilibration. This is in contrast to rehydration of the trehalose glass at r = 74% (NaCl/D2O), see Fig. S3 and S5 in the ESI.† The structural reorganization and partial polycrystallization of trehalose films at r = 74% result in nitroxide molecules preferentially located in water-rich domains.59 The increased plasticity of the rehydrated trehalose sugar allows for relaxation to an optimum geometry of nitroxide H-bonds with the surrounding water molecules. As a result, the hyperfine coupled deuterium lines in the ENDOR spectra are narrower and the quadrupole splitting is partially resolved, see Fig. S5 in the ESI.† Additionally, 14N EDNMR experiments reveal an increased fraction of nitroxide hydrogen bonded to two water molecules upon prolonged rehydration, see Fig. S3 and Table S2 in the ESI.†
![]() | ||
Fig. 4 (a) W-band EDNMR spectra of SL-bRC/trehalose glasses equilibrated at r = 11% (LiCl/D2O) and r = 74% (NaCl/D2O) (red and black traces, respectively). The orange and grey traces show the corresponding EDNMR spectra of the samples subsequently re-equilibrated at r = 11% (LiCl/H2O). The spectra were recorded at the gzz, MI = 0 spectral field position at 60 K indicated by the green arrow in Fig. 3b. Note that the spectra are shown in a frequency range different from that of Fig. 3a. (b) W-band Mims ENDOR spectra of SL-bRC/trehalose glasses equilibrated at r = 11% (LiCl/D2O) and r = 74% (NaCl/D2O). The spectra were recorded at the magnetic field position indicated in Fig. 3b with the blue arrow. |
Equilibrated to | r = 11% (D2O) from r = 11% (H2O) | r = 74% (D2O) from r = 11% (H2O) | r = 11% (H2O) from r = 11% (D2O) | r = 11% (H2O) from r = 74% (D2O) |
D-line EDNMR intensity | 0.097 ± 0.003 | 0.21 ± 0.08 | 0.037 ± 0.01 | 0.107 ± 0.006 |
The protein environment introduces additional complexity since the H/D exchange in amino acid residues on the protein surface is likely to add intensity to the deuterium line. This contribution is probably the reason for the only-a-factor-of-two difference in deuterium line intensity in the EDNMR spectra of SL-bRC/trehalose films at r = 11% and 74%. At both hydration levels, the water present in the embedding sugar clearly exchanges with the first solvation shell of the protein. To support this notion, SL-bRC/trehalose glasses at r = 74% (NaCl/D2O) and 11% (LiCl/D2O) were subsequently re-equilibrated at r = 11% (LiCl/H2O). The intensity of the D line significantly decreased in both cases (see the gray and orange traces in Fig. 4a and Table 1). The reason for the slower equilibration kinetics (not shown) and the presence of a residual D signal might originate in the kinetic isotope effect, but will not be further addressed in this work.
The first spectral feature of interest is the broad line centered at the 14N Larmor frequency (10.4 MHz). It stems from the peptide nitrogen NH of Ala M260 and/or the imidazole nitrogen Nδ(1)–H of His M219 in close proximity to the two carbonyl oxygens of QA.80,81 An additional signal is observed that is composed of a pair of lines centered at 2·νn(14N), i.e., 20.8 MHz at the present magnetic field, and split by twice the hyperfine coupling. It originates from double-quantum nitrogen transitions and, due to their lower probability, it has a significantly lower intensity.71 The next signal of interest centered at the deuterium Larmor frequency (21.9 MHz) is evident in the EDNMR spectra of Zn-bRCs in frozen D2O solution or Zn-bRC/trehalose rehydrated at r = 74% (NaCl/D2O). The intensity of the deuterium line in Zn-bRCs in the rehydrated trehalose glass and in frozen D2O solution is approximately the same. Clearly, the diffusion of water within the protein embedded in the rehydrated trehalose glass is not hampered in comparison to that of the protein in solution. Interestingly, in the case of the dry protein/trehalose glass (r = 11%, LiCl/D2O), the deuterium line is also observed, although its intensity is only ≈1/10 of that in frozen solution. This shows that the retarded conformational flexibility of the protein, resulting from trehalose coating, does not prevent the exchange of water molecules between the protein interior and the surrounding trehalose matrix. Moreover, the water exchange between the trehalose matrix and protein takes place on the same time scale as the water exchange between trehalose and the vapor atmosphere.
The large EDNMR linewidth of 2.4 ± 0.1 MHz (FWHM) does not allow us to draw conclusions on the origin of the D signal. Fig. 5c shows the comparison between the ENDOR spectra of Zn-bRCs equilibrated at r = 74% and in frozen D2O solution. Narrow ENDOR lines are observed in both cases.
The narrow ENDOR lines originate from hydrogen-bonded deuterium in well-defined geometry of the binding site. The ENDOR spectrum of can be analyzed in more detail than that of the nitroxide radical in a rigid matrix. At the chosen EPR field position of
the powder-like deuterium ENDOR spectrum is expected to consist of sharp features corresponding to parallel and perpendicular components of an axially symmetric hyperfine tensor. Two pairs of ENDOR lines are therefore seen, each of which is additionally split into a doublet due to the nuclear quadrupole interaction with deuterium. The hyperfine and quadrupole coupling constants for
in the Zn-bRC/trehalose glass equilibrated at r = 74% (NaCl/D2O) are the same as in the frozen D2O solution, i.e., they originate from the same type of coupled deuterium. Additionally, field-frequency swept W-band Mims ENDOR spectra were recorded across the
EPR spectrum of this sample (see Fig. S6 in the ESI†). Indeed, a well-defined dipolar ENDOR pattern is seen, suggesting a well-defined geometry of the hydrogen bond between the quinone and a H-bond donor of the binding site. From earlier ENDOR experiments of the Feher group80,81 it is known that
is stabilized in the protein pocket by two H-bonds between the quinone oxygens and protein residues Ala (M260) and His (M219). Notably, previous EPR and FTIR studies have shown that hydrogen bonds to the oxygens of quinone carbonyls are not equivalent.80,82–84 This allowed selectively deuterating one of the H-bonds to QA and extracting hyperfine and quadrupole coupling values for each of the deuterons.80,81 And from our earlier FTIR and W-band EPR work we know that embedding the bRC into the trehalose glass has no influence on the geometry of H-bonds stabilizing QA in the protein matrix.30
The comparison of the hyperfine/quadrupole values of this work with those reported earlier80 reveals that the observed ENDOR and EDNMR lines are most likely due to deuterons on O1 (Ala), see Table S3 in the ESI.† At first sight, it is surprising that no matrix water is observed in the ENDOR spectra of both frozen-solution and hydrated bRCs at r = 74% (NaCl/D2O). This finding is, however, in accordance with the high-resolution X-ray structure of bRCs that reveals the closest water cluster, consisting of 10 water molecules, is located only in distances between 0.5 and 1.2 nm from QA.44
Next, we report the results for as a spin probe. Fig. 5d shows the W-band EDNMR spectra of Zn-bRCs in frozen D2O solution, recorded at the maximum of the
EPR intensity and close to the gzz spectral position. In this spectral region there is no overlap of the
and
spectra. The broad line at 10.4 MHz, seen in both spectra, originates from 14N nuclei of the bacteriochlorophyll special-pair dimer85 and, possibly, from more distant amino acids located in its vicinity. The line pair centered at 2·νn(14N) and split by 6.4 MHz, i.e., twice the hyperfine coupling constant, stems from 14N double-quantum transitions.71,85 The high-frequency line of this pair partially overlaps with the deuterium line at 21.8 MHz, as can be seen in the EDNMR spectrum recorded at the maximum of the
EPR intensity. However, the EDNMR spectrum recorded at the shoulder of the
EPR spectrum shows only double-quantum 14N lines but no deuterium lines. This finding can be rationalized on the basis of the high-resolution X-ray structure of bRCs which reveals several water molecules with long occupancy times in the close vicinity of P865.44 These water molecules were ascribed to play a role in the primary electron transfer steps.86 The lack of a distinct deuterium signal in the EDNMR spectrum can be explained by the aromaticity of the spin probe
: there is an extended π-electron delocalization of unpaired spin density over the radical cation state
.87,88 This results in only a small π-electron spin density at any given nuclear position within the bacteriochlorophyll macrocyclic framework and, thus, in unresolved hyperfine coupling constants in frozen-solution samples (in contrast to liquid–solution
samples).87
![]() | ||
Fig. 6 (a) W-band EDNMR spectra of SL-bRC/trehalose glass samples equilibrated for 70 hours at r = 74% (NaCl/H217O) black trace; r = 11% (LiCl/H217O) magenta trace; and r = 11% (LiCl/H2O) red trace. The spectra were recorded at the same nitroxide spectral position (gzz, MI = 0) indicated in Fig. 3b by the green arrow. EDNMR spectra are shown for Δν < 0 to obtain a better separation between 17O and 14N signal contributions (see Fig. S2b in ESI†). (b) The 17O EDNMR spectra of r = 74% (NaCl/H217O) (black trace) and r = 11% (LiCl/H217O) (magenta trace) samples. The gray dashed lines show the best fit simulation of experimental recordings to a Gaussian line. (c) Microwave field amplitude dependence of the 17O EDNMR spectra of SL-bRC/trehalose glasses equilibrated for three days at r = 11% (LiCl/H217O). The black spectra were recorded using a HTA pulse length of tHTA = 20 μs. The blue trace shows the EDNMR spectrum acquired with tHTA = 7 μs. The green line shows the Davies ENDOR spectrum recorded at the same spectral position. |
The 17O line intensity at r = 11% is 6 ± 2 times smaller than that at r = 74% (intensities of 0.017 ± 0.005 and 0.11 ± 0.03, respectively). The difference in the water content is more pronounced than in the case of D2O equilibration where H/D exchange contributes to the D line intensity. At the 11% hydration level, the EDNMR signal is likely to originate from distant 17O nuclei, as judged from the Gaussian line shape. But because of the low signal-to-noise ratio of the spectrum, the presence of coupled 17O nuclei cannot be excluded. At the 74% hydration level, the 17O line in the EDNMR spectrum is not described any more by a Gaussian line shape as is the case of the deuterium line. Rather, resolved lines from matrix and hyperfine-coupled 17O nuclei can be discerned (see the upper trace in Fig. 6b).
It has been previously shown that by measuring the power dependence of EDNMR lines the different chemical environments of the nucleus of interest can be revealed.89 A similar approach was employed here to achieve a better distinction between hyperfine/quadrupole coupled and distant 17O nuclei. The coupled nuclei exhibit significantly larger hyperfine/quadrupole interactions with the unpaired electron spin of the nitroxide and, therefore, have intrinsically a larger transition probability for an EDNMR response, i.e., can be driven by a lower amplitude of the HTA mw pulse.89 The distant nuclei, on the other hand, have a much smaller hyperfine interaction with the unpaired electron spin and, thus, a high-amplitude HTA pulse is required to drive their transitions. Fig. 6c shows EDNMR spectra recorded for SL-RC/trehalose at r = 74% (H217O) with variable mw field amplitudes. At the highest mw field amplitude used, the intensity of the 17O EDNMR line is about 0.1, and both distant and coupled 17O nuclei are detected. When decreasing the mw field amplitude the overall signal intensity decreases, but the matrix 17O line is affected to a larger extent. To suppress the matrix line almost completely, a two times lower mw field amplitude and shorter HTA pulse length were used, see the blue trace in Fig. 6c.
Subsequently, W-band Davies ENDOR measurements were performed to confirm the 17O hyperfine coupling constant as extracted from the EDNMR spectrum, see green trace in Fig. 6c. In Davies ENDOR exclusively coupled water is probed since the signal from distant 17O nuclei is suppressed because of the blind spot at νn(17O).78 The presence of the 17O line in the EDNMR spectrum of the dehydrated protein/trehalose glass (r = 11%) unambiguously shows that, despite extensive dehydration and the high rigidity of the trehalose, the water retained in the trehalose sugar and the water on the protein surface can exchange.
The bRCs were embedded into a trehalose glass and equilibrated in an atmosphere of either low (r = 11%) or high (r = 74%) relative humidity. The resulting sugar matrix has very different properties: at low humidity and the corresponding low hydration level the protein–matrix system is rigid with only 0.5 water molecules per sugar molecule; at higher humidity and hydration level the sugar matrix is less rigid and contains more water molecules as potential partners for hydrogen bonding to the protein and to the sugar matrix. Different hydration levels, therefore, affect the protein dynamics differently: at r = 11%, the dynamics is blocked,54,59,60 at least on the timescale of charge recombination, but is not affected at r = 74%. Such an effect of the dry trehalose matrix, to accomplish inhibition of internal protein dynamics and, thereby, leading to bio-protection of proteins against thermal denaturation, reflects a strong “slaving” of the protein dynamics to that of the embedding matrix.20,21 The detailed molecular mechanism of bio-protection is still unclear and leaves open questions, for instance: why does the hydration of different embedding disaccharide sugars like trehalose or sucrose have a strongly different effect on the solvation shell of the protein and on the dynamics of the water molecules inside the protein?59
Several existing hypotheses (to some extent complementary and not mutually exclusive) try to explain the remarkable trehalose efficacy in bio-protection by predicting different changes in the first solvation shell of the protein upon dehydration. The “water replacement hypothesis”90 states that upon drying the hydrogen-bond interaction of polar groups at the protein interface to water is replaced by hydrogen bonds to the trehalose sugar by which the protein structure is strongly stabilized. The “preferential hydration hypothesis”,91 on the other hand, states that water molecules become entrapped at the interface between the glassy sugar matrix and the protein surface. The “high viscosity hypothesis”92 assumes that the inhibition of dynamic processes responsible for protein denaturation is mainly due to the high viscosity of the dried disaccharide matrix impeding diffusion. In an attempt to generalize the above models by incorporating some of their features, the “anchorage hypothesis”31 has been proposed stating that hindering of protein dynamics is caused by a tight anchoring of the protein surface to the rigid matrix through the formation of an extended H-bond network. According to this view, the disaccharide matrix locks the protein surface either directly (as assumed in the water replacement hypothesis) or through bridging residual water molecules at the protein–matrix interface, thus preventing large-scale protein internal motions. To differentiate experimentally between these hypotheses, a sensitive spectroscopic method capable of direct water detection in a site-specific manner was needed.
In the present work, pulsed high-field EPR techniques such as W-band ELDOR-detected NMR and ENDOR were used to unravel hydration–dehydration processes of bRC protein preparations embedded in trehalose matrices with controlled 2H and 17O isotope-labeled water contents. Specifically, the EDNMR method proved to be exceptionally sensitive in detecting hyperfine-coupled deuterium or 17O nuclei originating from D2O or H217O water, respectively, which were introduced into the sugar matrix in a controlled fashion by means of the isopiestic technique.
In trehalose-embedded SL-bRCs both EDNMR and ENDOR spectra show that water molecules, detected by a nitroxide spin label attached to the bRC surface, are retained in the first and second solvation shell, even under extensively drying conditions (r = 11%) which block the internal protein dynamics. Additionally, surface water molecules probed in SL-bRCs by the nitroxide radical, as well as inner water molecules probed in Zn-bRCs in the vicinity of the cofactor ion, exchange with water that is in equilibrium with the trehalose matrix. These observations strongly argue against the “water replacement hypothesis” sensu stricto, which assumes that the sugar binds directly to the polar groups of the protein surface, serving as a “water substitute” upon removal of the hydration shell.
As reported in the Results section, for SL-bRCs the ratio of 17O EDNMR line intensities for the hydrated state (r = 74%) to the dehydrated state (r = 11%) is 6 ± 2 reflecting the number of distant water molecules in the vicinity of the nitroxide spin label. This ratio is close to the H2O/trehalose molar ratios of both hydration states of 4.7 ± 0.2 as reported for the pure trehalose matrix.59 Thus, the scaling of the 17O signal intensity follows the changes in the overall water content in the trehalose matrix. This fact is in quantitative conflict with the “water replacement hypothesis”. At the same time, the above observation weakens the applicability of the “preferential hydration hypothesis”, because it does not give evidence for a preferential hydration of the protein surface in comparison with the average hydration of the whole sugar matrix. Thus, our findings appear to be consistent with the “anchorage hypothesis” assumption that the trehalose matrix immobilizes the protein by locking its surface through direct H-bonds between exposed protein residues and sugar molecules, as well as through indirect H-bonds mediated by bridging water molecules of the residual hydration shell. Clearly, a quantification of direct protein–trehalose versus water-mediated protein–matrix contacts would require the investigation of a large number of site-specifically nitroxide-labeled sites distributed over the protein surface, which was beyond the scope of the present work.
The exchange of protons to deuterons that are H-bonded to oxygens in the quinone even at a low hydration level is a surprising finding. Clearly, while progressive dehydration effectively blocks large-scale protein motions, the residual protein mobility still allows for exchange of stronger-bound water that is retained within the protein – and for eventual H/D exchange at the QA binding site.
While the spectral sensitivity of EDNMR relative to ENDOR is higher by up to several orders of magnitude, thus making it the method of choice, the limited resolution of deuterium EDNMR using D2O poses a disadvantage for local water sensing in a quantitative manner. However, deuterium EDNMR is well-suited for qualitative studies, like those presented here to determine the sugar rehydration kinetics for which no differentiation between coupled and distant water is necessary. When quantitative information on the different types of water in proteins is required, 17O EDNMR using H217O can be used instead, taking advantage of the high sensitivity of the method and the possibility to resolve coupled from distant 17O contributions. In this respect, W-band EDNMR in conjunction with D and /or 17O isotope labeling is indeed a very strong methodology for local water sensing.
As to the significance of such studies we want to emphasize that hydration water plays a crucial role in protein dynamics and structural relaxation on all time scales. Analysis suggests that changes in the amount of hydration water affect not only the protein's energy landscape but also significantly influence structural fluctuations between conformational substates – thereby controlling the biological function. To understand the functional difference between surface-bound water and bulk water is a key issue for understanding and controlling biological macromolecules, including proteins, DNA and RNA, right up to controlling anhydrobiotic food preservation and storage in dehydrated trehalose matrices at room temperature. Water is just far more than merely a solvent. Changes in water-mediated hydrogen-bonding patterns usually have a crucial impact on the global function of a biological system.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/c7cp03942e |
This journal is © the Owner Societies 2017 |