All relevant data are within the paper.

Introduction {#sec001}
============

Specialized proteins known as aquaporins are key to transport of water across biological membranes \[[@pone.0177273.ref001]\]. Interesting historical perspectives on the discovery of aquaporins can be found in \[[@pone.0177273.ref002]\] and \[[@pone.0177273.ref003]\]. Nuclear magnetic resonance (NMR) has played an important role in quantification and understanding of water exchange across cell membranes. A system that has been extensively studied by NMR are red blood cells \[[@pone.0177273.ref004], [@pone.0177273.ref005]\]. To measure water exchange in cell suspensions, NMR experiments require a contrast between intracellular and extracellular water compartments. Most commonly, this contrast is provided by differences in relaxation rates \[[@pone.0177273.ref006]\] or apparent diffusion coefficients (ADC) \[[@pone.0177273.ref007]\].

In relaxation based NMR approaches, the *T*~2~ \[[@pone.0177273.ref006]\] or *T*~1~ \[[@pone.0177273.ref008], [@pone.0177273.ref009]\] relaxation times of the extracellular water are selectively altered by a relaxation agent that cannot penetrate the cell membrane. The signal decay depends on the concentration of relaxation agent and on the exchange rate constant. The exchange can be estimated based on the two-compartment exchange model \[[@pone.0177273.ref008]\]. This approach is particularly suited to quantify relatively fast exchange, occurring on a 10 ms time-scale, like in the case of red blood cells \[[@pone.0177273.ref004], [@pone.0177273.ref005]\]. However, the invasive nature of this approach has been highlighted in a recent comparative study \[[@pone.0177273.ref010]\].

Diffusion NMR \[[@pone.0177273.ref011], [@pone.0177273.ref012]\] offers a non invasive alternative to measure exchange \[[@pone.0177273.ref013]\], relying on a naturally occurring difference in water ADC between the intra- and extracellular compartments. As measured by a pulsed gradient spin echo (PGSE) sequence, the intracellular water ADC, *D*~i~, is lower than that of the extracellular water, *D*~e~, due to restricted diffusion caused by the cell membranes \[[@pone.0177273.ref014]\]. The effects of exchange can be analysed using the Kärger model \[[@pone.0177273.ref013], [@pone.0177273.ref015]\]. Although the Kärger model \[[@pone.0177273.ref015]\] considers Gaussian diffusion in two exchanging compartments, it is used in several applications \[[@pone.0177273.ref007], [@pone.0177273.ref016]--[@pone.0177273.ref024]\]. The Kärger model can be modified to account for restricted diffusion \[[@pone.0177273.ref025]\] and for effects of different relaxation rate constants in the different compartments \[[@pone.0177273.ref026]--[@pone.0177273.ref028]\]. The rate of exchange is proportional to the membrane permeability, *P*. The average residence time for water molecules inside the cells, known as the intracellular lifetime, *τ*~i~, is inversely proportional to permeability, so that 1/*τ*~i~ = *PA*/*V*, where *A* is the membrane surface area and *V* is the intracellular volume \[[@pone.0177273.ref029], [@pone.0177273.ref030]\].

In contrast to the relaxation based approach, the diffusion exchange NMR is better suited for systems with relatively slow exchange, occurring on a time-scale of several 100 ms, which is inaccessible by the relaxation-based approach. In such conditions, a particularly advantageous method is the filter-exchange spectroscopy (FEXSY) introduced by Åslund *et al*. \[[@pone.0177273.ref019]\]. FEXSY is inspired by the diffusion-diffusion correlation experiment \[[@pone.0177273.ref017]\] using double diffusion encoding (DDE) \[[@pone.0177273.ref031]\] separated by a mixing time, *t*~m~. This allows to disentangle the effects of time-dependent diffusion and exchange. Compared to the full diffusion-diffusion correlation, FEXSY uses a much faster sparse sampling protocol \[[@pone.0177273.ref032]\], where the first diffusion encoding is applied as a low-pass diffusion filter that selectively attenuates the signal from compartments with high ADC. Furthermore, FEXSY analysis does not rely on inverse Laplace transform and thus avoids its inherent pitfalls \[[@pone.0177273.ref032]\].

FEXSY is suited to quantify intracellular lifetimes of approximately 200 to 2000 ms, such as for example in baker's yeast cells \[[@pone.0177273.ref019]--[@pone.0177273.ref021]\] or in human breast cells \[[@pone.0177273.ref020], [@pone.0177273.ref022]\]. Due to it's non-invasive nature, the FEXSY protocol can be adapted for constraints of human clinical and preclinical imaging systems \[[@pone.0177273.ref020], [@pone.0177273.ref033]\]. The imaging implementation of FEXSY, called filter exchange imaging (FEXI), can be applied *in vivo* to map apparent exchange rates in e.g. human brain \[[@pone.0177273.ref021], [@pone.0177273.ref024]\] or in human embryonic kidney cells \[[@pone.0177273.ref023]\].

The methodology proposed in this article was motivated by experiments designed to study water transport abilities of different types of human aquaporins. The results of this study will be presented elsewhere (Palmgren M, Hernebring M, Eriksson S, Elbing K, Geijer C, Lasič S, Dahl P, Hansen JS, Topgaard D, Lindkvist-Petersson K, unpublished work, 2017). Several model systems have been used to investigate the properties of water transport across the membranes of cells. In particular, *Xenopus laevis* oocytes and liposomes are systems commonly used to study the water transport activity of aquaporins \[[@pone.0177273.ref034], [@pone.0177273.ref035]\]. In such studies, water transport is commonly measured by detecting the swelling/shrinkage of the oocytes/liposomes upon an applied change in osmolality. In these experiments, water transport is quantified indirectly in terms of the osmotic water permeability coefficient. In contrast, the NMR methodology detects diffusional permeability coefficients. In our study, we use yeast cells (*Pichia pastoris*) as a model system. These cells can be modified using a methanol induced promoter to heterologously express human aquaporins \[[@pone.0177273.ref036], [@pone.0177273.ref037]\]. The advantage of using yeast cells compared to for instance red blood cells is the possibility to comparably easy remove all endogenous genes and insert genes of interest. In addition, yeast is a suitable system to over-express heterologous proteins at sufficient levels \[[@pone.0177273.ref037]\]. The intracellular lifetimes of these cells are expected to be at a 100 ms time-scale and extending beyond 1 s for cells without aquaporins, which is 10 to 100 times longer than in red blood cells. In such conditions, the FEXSY method offers an ideal opportunity to non-invasively study the transport mechanisms in equilibrium-exchange conditions.

In a multi-compartment system, the diffusion and relaxation effects can be different for each compartment. They can be quantitatively resolved only by encoding for both effects in a single two-dimensional *D* − *T*~2~ correlation experiment \[[@pone.0177273.ref032], [@pone.0177273.ref038]--[@pone.0177273.ref043]\]. In systems previously studied with FEXSY, the *T*~2~ relaxation was considered identical in all compartments \[[@pone.0177273.ref019]\]. However, different *T*~2~ relaxation for the intra- and extracellular water can cause bias when estimating the intracellular lifetime. In our preliminary experiments on genetically modified yeast cells, we observed large differences in the intra- and extracellular *T*~2~ values. For a quantitative comparative study, addressing the bias due to different relaxation times is crucial.

Here we introduce an augmented FEXSY protocol, which can be applied in comparative studies of aquaporins and cell cultures in general. To allow measurements of samples with short *T*~2~, we implemented the FEXSY sequence with a double pulsed gradient stimulated echo (double PGSTE) \[[@pone.0177273.ref044]\]. Most importantly, we present a method that can account for the confounding effects of different values of *T*~2~ and yield a bias-free estimate of the intracellular lifetime. The key to our novel method is a combination of the FEXSY and the PGSE based *D* − *T*~2~ correlation measurements. For our proof-of-concept, we used genetically modified *Pichia pastoris* yeast cells with deletion of the intrinsic aquaporins (*aqy1*Δ*agp1*Δ). For comparison, also cells with deletion of intrinsic and expression of human aquaporin (*aqy1*Δ*agp1*Δ + hAQP1) where used.

Theory {#sec002}
======

A schematic of the cell suspension system is shown in [Fig 1](#pone.0177273.g001){ref-type="fig"}. The system parameters for the intra- and extracellular compartments affecting the signal attenuation in the PGSE and FEXSY sequences are listed in [Fig 1C](#pone.0177273.g001){ref-type="fig"}. Cell membranes with and without aquaporins represent the *aqy1*Δ*agp1*Δ + hAQP1 and *aqy1*Δ*agp1*Δ yeast cells used in this article (see [Fig 1A and 1B](#pone.0177273.g001){ref-type="fig"}). The trajectories illustrate the different intracellular lifetimes for the two types of cells.

![Illustration of the intracellular and extracellular compartments separated by cell membranes in a cell suspension.\
(A) *aqy1*Δ*agp1*Δ + hAQP1 are yeast cells with aquaporins and (B) *aqy1*Δ*agp1*Δ are yeast cells without aquaporins. Diffusion trajectory of one water molecule is shown inside the cells. (C) The different parameters in the intracellular/extracellular compartments are: apparent diffusion coefficients, *D*~i~/*D*~e~, transverse relaxation times, *T*~2i~/*T*~2e~, lifetimes of water molecules, *τ*~i~/*τ*~e~.](pone.0177273.g001){#pone.0177273.g001}

Our analysis uses an approximate and phenomenological approach to quantify the effects of diffusion, relaxation and exchange in a two-compartment system. We apply a similar mathematical formalism as in the stochastic treatment of relaxation in multiple phase systems \[[@pone.0177273.ref045]\] or in the analysis of multi-site chemical exchange based on the matrix formulation of the Bloch equations \[[@pone.0177273.ref046]\]. The time evolution of a spin-system, e.g. undergoing relaxation, is rigorously described in terms of the density operator subjected to the Liouville---von Neuman equation using evolution superoperators \[[@pone.0177273.ref047]\], often also called propagators \[[@pone.0177273.ref048]\]. In our analysis, we use operators in a more general mathematical sense, which here take the form of matrix exponentials. The signal evolution can thus be described as a composite of sequential operators, which include effects of diffusion, relaxation and exchange.

During an NMR pulse sequence, the magnetization experiences several coherence pathways, which are selected by phase cycling of the radio frequency pulses \[[@pone.0177273.ref011]\]. Rather than following the longitudinal and traverse components of different magnetization pathways, we consider the resulting signal magnitudes from intracellular and extracellular compartments, denoted with subscripts i and e, respectively. The acquired signal is given by contributions from the two compartments as *S* = **1** ⋅ **S**, where **S** = (*S*~i~, *S*~e~)^T^ and **1** = (1, 1). For convenience of notation, we use the transpose operator, denoted by the superscript *T*, to change row vectors into column vectors. With **S**~0~ = (*S*~0i~, *S*~0e~)^T^ we denote signal contributions that are not affected by diffusion, *T*~2~ relaxation or exchange. Assuming equal *T*~1~ for both contributions, the effects of *T*~1~ can be factored into **S**~0~. We further assume that the gradient pulses are short relative to the relaxation and exchange time-scales. This is a good approximation in our experimental conditions considering gradient pulses of 1.2 ms and it is used in the well-known and commonly applied Kärger model \[[@pone.0177273.ref013], [@pone.0177273.ref015]\].

The theoretical section is divided in three parts: (i) introducing operators to calculate the signal evolution due to diffusion, relaxation and exchange; (ii) applying the operators to predict signal attenuation in PGSE ([Fig 2](#pone.0177273.g002){ref-type="fig"}) and FEXSY ([Fig 3](#pone.0177273.g003){ref-type="fig"}) sequences, which is used in data fitting; (iii) obtaining analytic expressions to estimate the bias of the intracellular lifetime as a function of the difference in relaxation rate constants.

![Schematic representation of the PGSE experiment.\
(A) The PGSE pulse sequence with gradient pulse duration *δ*, interval between the onset of the gradient pulses Δ, effective diffusion encoding time *t*~D~, gradient strength *g*, and echo time *t*~E~. (B) The different patterned blocks mark the time intervals when the signal is subjected to the evolution operators for exchange, *T*~2~-relaxation and diffusion. The coloured blocks show the time intervals during which the operators defined in Eqs [9](#pone.0177273.e010){ref-type="disp-formula"} and [10](#pone.0177273.e011){ref-type="disp-formula"} apply.](pone.0177273.g002){#pone.0177273.g002}

![Schematic representation of the FEXSY experiment.\
(A) The FEXSY sequence with gradient pulse duration *δ*, interval between the onset of the gradient pulses Δ, effective diffusion encoding time *t*~D~, gradient strength of the filter block *g*~f~, gradient strength of the diffusion encoding block *g*, and the mixing time *t*~m~. Magnetization is affected by longitudinal relaxation during times *t* and *t*~1~ and by transverse relaxation during echo time, *t*~E~. (B) The different patterned blocks mark the time intervals when the signal is subjected to the evolution operators for exchange, *T*~2~-relaxation, and diffusion. The coloured blocks show the time intervals during which the operators defined in Eqs [12](#pone.0177273.e013){ref-type="disp-formula"}--[15](#pone.0177273.e016){ref-type="disp-formula"} apply.](pone.0177273.g003){#pone.0177273.g003}

Operators for diffusion, relaxation and exchange {#sec003}
------------------------------------------------

### Diffusion {#sec004}

To encode diffusion information, both the PGSE ([Fig 2](#pone.0177273.g002){ref-type="fig"}) and the FEXSY ([Fig 3](#pone.0177273.g003){ref-type="fig"}) sequences employ pulsed magnetic field gradients during the diffusion encoding blocks. The signal attenuation is controlled by the pulse duration, *δ*, separation in time, Δ, and amplitude, *g*, of the gradient pulses (see Figs [2](#pone.0177273.g002){ref-type="fig"} and [3](#pone.0177273.g003){ref-type="fig"}). For low *g*, the signal attenuation is proportional to the diffusion weighting *b* = *q*^2^ *t*~D~, where *q* = *γgδ*, *γ* the gyromagnetic ratio and *t*~D~ = (Δ − *δ*/3) the effective diffusion time. Considering only diffusion, the signal evolution is given by $$\begin{array}{r}
{\mathbf{S} = e^{- q^{2}\mathbf{D}t_{D}}\mathbf{S}_{0},} \\
\end{array}$$ where $$\begin{array}{r}
{\mathbf{D} = \begin{pmatrix}
D_{i} & 0 \\
0 & D_{e} \\
\end{pmatrix}.} \\
\end{array}$$ Note that the restricted diffusion effects are accounted for by the ADCs, *D*~i~ and *D*~e~, which depend on *δ* and *t*~D~ \[[@pone.0177273.ref049]\].

### Relaxation {#sec005}

When the intracellular transverse relaxation times, *T*~2i~ and *T*~2e~ are different, the signal evolution during the transverse encoding echo time, *t*~E~, is given by $$\begin{array}{r}
{\mathbf{S} = e^{- \mathbf{R}t_{E}}\mathbf{S}_{0},} \\
\end{array}$$ where $$\begin{array}{r}
{\mathbf{R} = \begin{pmatrix}
R_{2i} & 0 \\
0 & R_{2e} \\
\end{pmatrix}.} \\
\end{array}$$ Here *R*~2i,e~ = 1/*T*~2i,e~ are the relaxation rate constants for the respective compartments. The difference in relaxation rate constants is given by Δ*R*~2~ = *R*~2i~ − *R*~2e~.

### Exchange {#sec006}

We describe the effects of exchange in a two-compartment system according to the first order kinetics \[[@pone.0177273.ref045], [@pone.0177273.ref050]\]. Due to exchange effects alone, occurring during a time interval *t*, the signal evolution is given by $$\begin{array}{r}
{\mathbf{S} = e^{- \mathbf{K}t}\mathbf{S}_{0},} \\
\end{array}$$ where $$\begin{array}{r}
{\mathbf{K} = \begin{bmatrix}
k_{\text{ie}} & {- k_{\text{ei}}} \\
{- k_{\text{ie}}} & k_{\text{ei}} \\
\end{bmatrix}.} \\
\end{array}$$ Here *k*~ie~ and *k*~ei~ are the intracellular to extracellular and the reverse exchange rate constants, respectively. The corresponding lifetimes are given by *τ*~i~ = 1/*k*~ie~ and *τ*~e~ = 1/*k*~ei~. The exchange matrix **K** fulfills the detailed equilibrium condition \[[@pone.0177273.ref051]\] $$\begin{array}{r}
{{\mathbf{K}\mathbf{S}}_{0} = 0.} \\
\end{array}$$

The intra- and extracellular signal fractions are given by *f*~i~ = *S*~i~/(*S*~e~ + *S*~i~) and *f*~e~ = 1 − *f*~i~. Considering the equilibrium condition in [Eq 7](#pone.0177273.e007){ref-type="disp-formula"} and the total exchange rate constant, *k* = *k*~ie~ + *k*~ei~, the exchange matrix can be expressed as $$\begin{array}{r}
{\mathbf{K} = k\begin{bmatrix}
f_{e}^{\text{eq}} & {f_{e}^{\text{eq}} - 1} \\
{- f_{e}^{\text{eq}}} & {1 - f_{e}^{\text{eq}}} \\
\end{bmatrix},} \\
\end{array}$$ where $f_{e}^{\text{eq}}$ is the extracellular signal fraction at equilibrium.

Signal evolution during PGSE and FEXSY {#sec007}
--------------------------------------

For the PGSE sequence, the signal evolution is given by the composite of operators *O*~**P1**~ and *O*~**P2**~ acting during three intervals: before(*O*~**P1**~), during(*O*~**P2**~) and after(*O*~**P1**~) the diffusion encoding (see [Fig 2](#pone.0177273.g002){ref-type="fig"}). We assume gradient pulses to be short, which is reflected in the instantaneous changes in the coloured blocks shown at the bottom of [Fig 2](#pone.0177273.g002){ref-type="fig"} representing the three time intervals mentioned above. The *O*~**P1**~ operator includes the effects of relaxation and exchange, $$\begin{array}{r}
{\mathit{O}_{\mathbf{P}1} = e^{- \frac{1}{2}{(\mathbf{R} + \mathbf{K})}{(t_{E} - t_{D})}}.} \\
\end{array}$$ The operator *O*~**P2**~ includes the effects of relaxation, exchange and diffusion, $$\begin{array}{r}
{\mathit{O}_{\mathbf{P}2} = e^{- {(\mathbf{R} + \mathbf{K} + q^{2}\mathbf{D})}t_{D}}.} \\
\end{array}$$ The signal evolution in PGSE is given by $$\begin{array}{r}
{\mathbf{S}_{\text{PGSE}} = \mathit{O}_{\mathbf{P}1}\mathit{O}_{\mathbf{P}2}\mathit{O}_{\mathbf{P}1}\mathbf{S}_{0}.} \\
\end{array}$$ The PGSE sequence can be used to estimate Δ*R*~2~ based on the diffusion contrast between the two compartments. The diffusion and relaxation effects can be disentangled in a *D* − *T*~2~ correlation experiment \[[@pone.0177273.ref032], [@pone.0177273.ref043], [@pone.0177273.ref052], [@pone.0177273.ref053]\], where the diffusion weighting is kept constant, by using constant *δ* and Δ, while varying the relaxation weighting through *t*~E~.

For the FEXSY sequence, the signal evolution is given by the composite of operators *O*~**F1**~, *O*~**F2**~, *O*~**F3**~ and *O*~**F4**~. The time intervals during which these operators act and the effects that they include are depicted in [Fig 3](#pone.0177273.g003){ref-type="fig"}. Also in this case the assumption of short gradient pulses is reflected in the coloured blocks shown at the bottom of [Fig 3](#pone.0177273.g003){ref-type="fig"}. The operator *O*~**F1**~ includes the effects of relaxation and exchange, $$\begin{array}{r}
{\mathit{O}_{\mathbf{F}1} = e^{- \frac{1}{2}{(\mathbf{R} + \mathbf{K})}{(t_{E}/2 + t_{1} - t_{D})}},} \\
\end{array}$$ the operator *O*~**F2**~ includes the effects of relaxation, exchange and diffusion $$\begin{array}{r}
{\mathit{O}_{\mathbf{F}2} = e^{- \frac{1}{2}{(\mathbf{R} + \mathbf{K} + q^{2}\mathbf{D})}{(t_{D} - t_{1})}},} \\
\end{array}$$ the operator *O*~**F3**~ includes the effects of exchange and diffusion $$\begin{array}{r}
{\mathit{O}_{\mathbf{F}3} = e^{- {(\mathbf{K} + q^{2}\mathbf{D})}t_{1}},} \\
\end{array}$$ and the operator *O*~**F4**~ includes only the effects of exchange $$\begin{array}{r}
{\mathit{O}_{\mathbf{F}4} = e^{- \mathbf{K}t}.} \\
\end{array}$$ Note that FEXSY is specifically sensitive to exchange due to the isolated exchange effects in operator *O*~**F4**~. The signal evolution in FEXSY is given by $$\begin{array}{r}
{\mathbf{S}_{\text{FEXSY}} = \mathit{O}_{\mathbf{F}1}\mathit{O}_{\mathbf{F}2}\mathit{O}_{\mathbf{F}3}\mathit{O}_{\mathbf{F}2}\mathit{O}_{\mathbf{F}1}\mathit{O}_{\mathbf{F}4}\mathit{O}_{\mathbf{F}1}\mathit{O}_{\mathbf{F}2}\mathit{O}_{\mathbf{F}3}\mathit{O}_{\mathbf{F}2}\mathit{O}_{\mathbf{F}1}\mathbf{S}_{0}.} \\
\end{array}$$

Bias of the intracellular lifetime in the presence of relaxation rate differences {#sec008}
---------------------------------------------------------------------------------

Let's examine the bias of the exchange rate constant and the intracellular lifetime estimated by the FEXSY sequence in the presence of the difference in relaxation rate constants Δ*R*~2~. To obtain tractable analytical expressions in this step of analysis, we ignore the effects of exchange during the diffusion encoding blocks, similar as in the standard FEXSY analysis by Åslund et al. \[[@pone.0177273.ref019]\]. Hence, the patterned block representing exchange in [Fig 3B](#pone.0177273.g003){ref-type="fig"} does not apply in this analytical derivation. The diffusion effects can be accounted for implicitly by considering the extracellular signal fraction at equilibrium, $f_{e}^{\text{eq}}$, and the fraction perturbed by the diffusion filter block, *f*~e0~ at *t*~m~ = 0. Note that diffusion encoding is necessary to experimentally determine the signal fractions.

The signal contributions are given by $$\begin{array}{r}
{\mathbf{S}_{r} = S_{0}e^{- \mathbf{R}t_{E}}e^{- \mathbf{K}t_{m}}e^{- \mathbf{R}t_{E}}\left( 1 - f_{e0},f_{e0} \right)^{T},} \\
\end{array}$$ where *S*~0~ is the normalization factor. The relaxation-affected signal fraction is given by [Eq 17](#pone.0177273.e019){ref-type="disp-formula"} as $$\begin{array}{r}
{f_{e,r} = \mathbf{S}_{r}\left( 2 \right)/\left( \mathbf{S}_{r}\left( 1 \right) + \mathbf{S}_{r}\left( 2 \right) \right).} \\
\end{array}$$

Eqs [8](#pone.0177273.e008){ref-type="disp-formula"}, [17](#pone.0177273.e019){ref-type="disp-formula"} and [18](#pone.0177273.e020){ref-type="disp-formula"} yield the expression for the relaxation-affected extracellular fraction given by $$\begin{array}{r}
\begin{array}{cl}
 & {f_{e,r}\left( f_{e}^{\text{eq}},f_{e0},kt_{m},r \right) =} \\
 & \left\lbrack \left( f_{e}^{\text{eq}} - f_{e}^{\text{eq}}e^{- kt_{m}} \right)\left( f_{e0} - 1 \right) - f_{e0}e^{r}\left( f_{e}^{\text{eq}} + e^{- kt_{m}} - f_{e}^{\text{eq}}e^{- kt_{m}} \right) \right\rbrack \\
 & \left\lbrack \left( f_{e}^{\text{eq}} - f_{e}^{\text{eq}}e^{- kt_{m}} \right)\left( f_{e0} - 1 \right) - f_{e0}e^{r}\left( f_{e}^{\text{eq}} + e^{- kt_{m}} - f_{e}^{\text{eq}}e^{- kt_{m}} \right) \right. \\
 & {+ e^{r}\left( f_{e0} - 1 \right)\left( f_{e}^{\text{eq}}e^{- kt_{m}} - f_{e}^{\text{eq}} + 1 \right) + f_{e0}e^{- kt_{m}}\left( e^{kt_{m}} - 1 \right)\left( f_{e}^{\text{eq}} - 1 \right)\rbrack^{- 1},} \\
\end{array} \\
\end{array}$$ where *r* = Δ*R*~2~*t*~E~/2. The function in [Eq 19](#pone.0177273.e021){ref-type="disp-formula"} is shown graphically in [Fig 4](#pone.0177273.g004){ref-type="fig"}. For Δ*R*~2~ = 0, the asymptotic value of *f*~e~ at long *t*~m~ corresponds to the experimentally accessible value of *f*~e~ at short *t*~m~ and without the diffusion filtering \[[@pone.0177273.ref019]\]. This correspondence of *f*~e~ at short and long *t*~m~ no longer exists for \|Δ*R*~2~\| \> 0 (compare solid and dashed lines in [Fig 4](#pone.0177273.g004){ref-type="fig"}).

![Relaxation-affected extracellular signal fraction as a function of *kt*~m~ for different values of Δ*R*~2~*t*~E~.\
Solid lines represent the case when diffusion filter gradients are applied in the first PGSTE block, while dashed lines are for the case without the filter gradients (see [Fig 3](#pone.0177273.g003){ref-type="fig"}). The results for different values of Δ*R*~2~*t*~E~ are labeled on the right ranging from −5 to +5 and color coded from blue to red for negative values and from red to green for positive values. In all cases, *f*~e0~ = 0.05 and $f_{e}^{\text{eq}} = 0.5$.](pone.0177273.g004){#pone.0177273.g004}

In contrast to the case with Δ*R*~2~ = 0, when *f*~e~(*t*~m~) has an exponential dependence on *t*~m~ (see Eq 7 in \[[@pone.0177273.ref019]\]), the function *f*~e,*r*~(*t*~m~) does not have such a simple form. However, we can estimate the relaxation-affected exchange rate constant, *k*~r~, from the time dependence of *f*~e,*r*~ at short *t*~m~, $$\begin{array}{r}
{k_{r} = - \lim\limits_{t_{m}\rightarrow 0}\frac{\partial}{\partial t_{m}}\ln\left\lbrack \frac{f_{e,r}^{\text{eq}} - f_{e,r}\left( t_{m} \right)}{f_{e,r}^{\text{eq}} - f_{e0,r}} \right\rbrack.} \\
\end{array}$$

The [expression 20](#pone.0177273.e022){ref-type="disp-formula"} resembles the definition of the apparent exchange rate constant given by the time dependence of the ADC at short *t*~m~ (see Eq 11 in \[[@pone.0177273.ref020]\]). The $f_{e,r}^{\text{eq}}$ in [Eq 20](#pone.0177273.e022){ref-type="disp-formula"}, determined as the asymptotic value of the relaxation-affected extracellular fraction at long *t*~m~, can be expressed with the relaxation non-affected value $f_{e}^{\text{eq}}$ as $$\begin{array}{r}
{f_{e\infty,r}^{\text{eq}} = \lim\limits_{t_{m}\rightarrow\infty}f_{e,r}\left( t_{m} \right) = \frac{f_{e}^{\text{eq}}}{f_{e}^{\text{eq}} + \left( 1 - f_{e}^{\text{eq}} \right)e^{- r}}.} \\
\end{array}$$ When, on the other hand, $f_{e,r}^{\text{eq}}$ is determined at *t*~m~ = 0 and no diffusion filtering is applied in the first PGSTE block, as in our experiments, we have $$\begin{array}{r}
{f_{e0,r}^{\text{eq}} = \lim\limits_{t_{m}\rightarrow 0}f_{e,r}\left( t_{m} \right) = \frac{f_{e}^{\text{eq}}}{f_{e}^{\text{eq}} + \left( 1 - f_{e}^{\text{eq}} \right)e^{- 2r}}.} \\
\end{array}$$ A similar expression is obtained for the relaxation-affected filtered fraction $$\begin{array}{r}
{f_{e0,r} = \frac{f_{e0}}{f_{e0} + \left( 1 - f_{e0} \right)e^{- 2r}}.} \\
\end{array}$$

It is interesting that the effect of relaxation difference, expressed by the factor exp(−2*r*) in Eqs [22](#pone.0177273.e027){ref-type="disp-formula"} and [23](#pone.0177273.e028){ref-type="disp-formula"}, is reduced to the factor exp(−*r*) in [Eq 21](#pone.0177273.e025){ref-type="disp-formula"}. This is a consequence of exchange between two compartments with different *R*~2~ occurring during long *t*~m~, which cancels the effect of Δ*R*~2~ during the first PGSTE block. Compare the dashed/solid lines at short/long *t*~m~ in [Fig 4](#pone.0177273.g004){ref-type="fig"}. Note that even though the effects of Δ*R*~2~ are different in the filtered and non-filtered experiments, the difference vanishes at long *kt*~m~.

To investigate the bias of *k*~*r*~ caused by Δ*R*~2~, we consider two cases. In the first case, we evaluate [Eq 20](#pone.0177273.e022){ref-type="disp-formula"} using the non-filtered signal fraction at long *t*~m~, $f_{e\infty,r}^{\text{eq}}$ from [Eq 21](#pone.0177273.e025){ref-type="disp-formula"} and the filtered fraction *f*~e0,*r*~ from [Eq 23](#pone.0177273.e028){ref-type="disp-formula"} yielding $$\begin{array}{r}
{k_{r}^{\infty} = k\frac{\left( f_{e0}e^{r} - f_{e0} + 1 \right)\left( f_{e}^{\text{eq}}e^{r} - f_{e}^{\text{eq}} + 1 \right)}{\left( f_{e0}e^{2r} - f_{e0} + 1 \right)}.} \\
\end{array}$$

In the second case, when the experimentally accessible $f_{e0,r}^{\text{eq}}$ from [Eq 22](#pone.0177273.e027){ref-type="disp-formula"} is considered, we get $$\begin{array}{r}
{k_{r}^{0} = \frac{- ke^{r}\left( f_{e0}e^{r} - f_{e0} + 1 \right)\left( f_{e}^{\text{eq}}e^{r} - f_{e}^{\text{eq}} + 1 \right)\left( f_{e}^{\text{eq}} - f_{e0}f_{e}^{\text{eq}} - f_{e0}e^{r} + f_{e0}f_{e}^{\text{eq}}e^{r} \right)}{\left( f_{e0} - f_{e}^{\text{eq}} \right)\left( f_{e0}e^{r} - f_{e0} + 1 \right)}.} \\
\end{array}$$

From Eqs [24](#pone.0177273.e031){ref-type="disp-formula"} and [25](#pone.0177273.e033){ref-type="disp-formula"} we can calculate how the relaxation-affected intracellular lifetime, *τ*~i,*r*~ is related to the ground truth value *τ*~i~, i.e. the unbiased value that we wish to experimentally quantify. In [Fig 5](#pone.0177273.g005){ref-type="fig"}, the ratio of *τ*~i,*r*~/*τ*~i~ is shown as a function of Δ*R*~2~*t*~E~ for two different values of *f*~e0~ (rows) and three different values of $f_{e}^{\text{eq}}$ (columns). The red and blue curves correspond to Eqs [24](#pone.0177273.e031){ref-type="disp-formula"} and [25](#pone.0177273.e033){ref-type="disp-formula"}, respectively.

![Relative bias of the intracellular lifetime as a function of Δ*R*~2~*t*~E~.\
The ratio of *τ*~i,*r*~/*τ*~i~ as a function of Δ*R*~2~*t*~E~, is shown for six different sets of *f*~e0~ and $f_{e}^{\text{eq}}$ values. The solid blue lines are calculated from [Eq 25](#pone.0177273.e033){ref-type="disp-formula"}, while the dashed red lines are calculated from [Eq 24](#pone.0177273.e031){ref-type="disp-formula"}. The framed plot shows the bias from low *f*~e0~ and 50% extracellular fraction, which approximately corresponds to our experimental conditions. The vertical dotted line is at Δ*R*~2~*t*~E~ = 0.6, which corresponds to echo time *t*~E~ = 17 ms and the difference in relaxation rate constants Δ*R*~2~ = 35 s^−1^. Here the relative bias (solid blue line) is about −20%.](pone.0177273.g005){#pone.0177273.g005}

It is clear from [Fig 5](#pone.0177273.g005){ref-type="fig"} that the relative bias *τ*~i,*r*~/*τ*~i~ as a function of Δ*R*~2~*t*~E~ depends on *f*~e0~ and $f_{e}^{\text{eq}}$. This dependence is particularly pronounced when the equilibrium fraction is estimated at short mixing times (solid blue curves). The gray dotted line in the framed plot corresponds to our experimental value of Δ*R*~2~*t*~E~ = 0.6. When using the standard FEXSY analysis, ignoring the differences in relaxation rate constants and estimating the equilibrium fraction at short mixing times and for low *f*~e0~, the intracellular lifetime would in this case be underestimated by approximately 20%. However, for larger *f*~e0~ (see [Fig 5E](#pone.0177273.g005){ref-type="fig"}) the bias would only be −8%.

Samples {#sec09}
-------

Two different yeast cell strains *aqy1*Δ*agp1*Δ + hAQP1 and *aqy1*Δ*agp1*Δ were prepared by cultivating and inducing the yeast cells in methanol and harvesting after several hours. *Pichia pastoris* strains used in this study derive from the X33 strain (Invitrogen). The *aqy1*Δ*agp1*Δ + hAQP1 strain expresses human aquaporins while the *aqy1*Δ*agp1*Δ strain does not (see [Fig 1](#pone.0177273.g001){ref-type="fig"}).

The human aquaporins AQP1 with a His6 purification tag were cloned into pPICZB vector using restriction enzyme EcoRI and XbaI and transformed into *Pichia pastoris*. To identify clones that express the protein of interest, AQP transformants were spotted on YPD-agar plates supplemented with increasing zeocin concentrations (0, 0.5, 1.0 and 2.0 mg/ml). Cells were pre-grown in BMGY for 24 hours and protein production was induced with methanol to a final concentration of 0.25% (v/v) for 1 hour. The cells were harvested and washed with 20 mM HEPES pH 7.5, re-suspended in a 1:1 (v/v) solution with buffer, and transferred to 5 mm disposable NMR tubes. The tubes were centrifuged at 1000 *g* for 10 minutes to achieve a pellet with high concentration of cells (approximately 40% intracellular water). The excess fluid was removed from the tubes. The samples were kept on ice until measurement. Measurements were performed at 0 and 20°C.

NMR experiment {#sec010}
--------------

The NMR experiments were performed on a Bruker 200 MHz Avance II system equipped with a DIF-25 probe capable of z-gradients up to 9.6 Tm^−1^. The temperature was controlled by a Bruker variable temperature unit with thermostatic air flow (±0.1°C). The samples were allowed to equilibrate inside the magnet 15 minutes before the experiments started.

The *D* − *T*~2~ correlation experiment ([Fig 2](#pone.0177273.g002){ref-type="fig"}) and the FEXSY experiment ([Fig 3](#pone.0177273.g003){ref-type="fig"}) were used for all samples. For the PGSE experiment, 10 different *t*~E~ values were used, logarithmically spaced from 5.5 to 82 ms. For each *t*~E~, the gradients were incremented logarithmically in 10 steps from 0.096 to 9.6 T m^−1^, yielding *b* values from 6.1 ⋅ 10^6^ to 2.2 ⋅ 10^10^ s m^−2^. The diffusion encoding timing parameters *δ* and Δ were kept constant during the experiment at 1.2 ms and 2.7 ms, respectively. The short Δ used here was necessary to achieve the shortest echo time, *t*~E~ = 5.5 ms, which allowed estimating short *T*~2~ values in the intracellular compartment.

For the FEXSY experiment, a double stimulated echo was used to minimize signal attenuation due to short *T*~2~. The sequence was performed with a 16-step phase cycle. The 8 steps of the cycle were used to eliminate artifacts that arise from unwanted coherence pathways as demonstrated by Khrapitchev and Callaghan \[[@pone.0177273.ref054]\]. The cycle was repeated 2 times with alternating direction of the magnetization z-storage by alternating the phase of the fourth 90° pulse. The signals from these two cycles were subtracted from each other to eliminate the contribution from the relaxed magnetization. The diffusion encoding time parameters were equal in both blocks and kept constant, *δ* = 1.2 ms and Δ = 8.4 ms. The strength of the filter gradient, *g*~f~, was set to 2.88 T m^−1^, resulting in *b*~f~ = 4.4 ⋅ 10^9^ s m^−2^. In the second diffusion encoding block the gradients were incremented logarithmically in 16 steps between 0.23 and 7.7 T m^−1^ yielding *b* values from 4.4 ⋅ 10^7^ to 4.9 ⋅ 10^10^ s m^−2^. The mixing times were 14, 54, 94 and 214 ms. For the non-filtered acquisition, used to estimate $f_{e0}^{\text{eq}}$, a filter gradient of 1% of *g*~f~ was used at *t*~m~ = 14 ms. The low filter gradient was used in addition to the 16 phase cycles to avoid contributions from unwanted magnetization pathways. An inversion recovery experiment with 16 recovery times between 10 and 2000 ms was performed to check that *T*~1~ was single-exponential.

Data analysis {#sec011}
-------------

All data analysis was performed with Matlab (The MathWorks Inc., Natick, MA). Signal intensities were obtained by integrating over the water peak in the NMR spectra. A constrained global fit of Eqs [11](#pone.0177273.e012){ref-type="disp-formula"} and [16](#pone.0177273.e017){ref-type="disp-formula"} was made to the signal intensities from PGSE and FEXSY using a *lsqcurvefit* Matlab routine. Since the acquisition parameters were identical in both experiments, no signal weighting of the two datasets was necessary in the optimization routine. The common fit parameters, constrained to be identical for both datasets, were *τ*~i~, *τ*~e~, *T*~2i~ and *T*~2e~. Through the relationship shown in [Eq 8](#pone.0177273.e008){ref-type="disp-formula"} the value of $f_{e}^{\text{eq}}$ can be retrieved by $f_{e}^{\text{eq}} = \tau_{e}/\left( \tau_{i} + \tau_{e} \right)$. To account for different diffusion times in the two experiments $D_{i}^{\text{PGSE}}$, $D_{e}^{\text{PGSE}}$, $D_{i}^{\text{FEXSY}}$ and $D_{e}^{\text{FEXSY}}$ were used as free fit parameters. A single *S*~0~ value was used in PGSE, while *S*~0~(*t*~m~) was used in FEXSY to account for *T*~1~ relaxation. *S*~0~ here refers to the total initial signal, *S*~0~ = *S*~0i~ + *S*~0e~.

The errors in the fit parameters were estimated by least-square analysis of 1500 datasets obtained by bootstrap re-sampling \[[@pone.0177273.ref055]\] of the two experimental datasets, together consisting of 150 data points. To eliminate outlier fits, only the 85% of the fits with lowest *χ*^2^ were retained. From the retained fits, the means and standard deviations of the fit parameters were calculated. All the errors presented in this article refer to the standard deviation.

Results and Discussion {#sec012}
======================

Due to the difference in relaxation rate constants, Δ*R*~2~, the standard FEXSY protocol yields biased estimates of *τ*~i~, as shown by the dashed red line in [Fig 5](#pone.0177273.g005){ref-type="fig"}. To avoid the bias, we introduced a complementary measurement of *D* − *T*~2~ correlation by a PGSE sequence with varying echo time. This measurement is specifically sensitive to *T*~2~ in the intracellular and extracellular compartments and can thus be used in combination with FEXSY to yield accurate quantification of *τ*~i~. The analysis is done by a global constrained fitting of Eqs [11](#pone.0177273.e012){ref-type="disp-formula"} and [16](#pone.0177273.e017){ref-type="disp-formula"} to the entire dataset consisting of PGSE and FEXSY data. The model accounts for both exchange and relaxation, *τ*~i,e~ and *T*~2i,e~, which are constrained to be identical for both experiments. The data and the fitting results for three experiments (*aqy1*Δ*agp1*Δ @ 0°C, *aqy1*Δ*agp1*Δ @ 20°C and *aqy1*Δ*agp1*Δ + hAQP1 @ 0°C) are shown in three columns of [Fig 6](#pone.0177273.g006){ref-type="fig"}. The upper row refers to the PGSE experiment and the lower row to the FEXSY experiment.

![Signal vs. *b* from PGSE and FEXSY.\
Top row shows results of the PGSE experiment for increasing echo times, *t*~E~, indicated by the arrow. Bottom row shows results of the FEXSY experiment, where the arrow indicates the increasing mixing times, *t*~m~. Columns show results for different yeast strains at different temperatures. Columns 1 and 2: *aqy1*Δ*agp1*Δ measured at T = 20°C and at T = 0°C. Column 3: *aqy1*Δ*agp1*Δ + hAQP1 measured at T = 0°C. The result of a constrained global fit of Eqs [11](#pone.0177273.e012){ref-type="disp-formula"} and [16](#pone.0177273.e017){ref-type="disp-formula"} to the entire dataset (upper and lower plots) is shown by the solid black lines. For the PGSE experiment, all the signal intensities are normalized by the maximum *S*~0~ value at shortest *t*~E~, while for the FEXSY experiment, different normalization constants, *S*~0~(*b*~f~, *t*~m~), are used for different *b*~f~ and *t*~m~ values.](pone.0177273.g006){#pone.0177273.g006}

To maximize the contrast between *D*~i~ and *D*~e~ in FEXSY, longer *t*~D~ is favourable. On the other hand, a short *t*~D~ is favoured in the PGSE sequence to minimize the effect of exchange during the diffusion encoding and allow for the shortest possible *t*~E~ to be used. In our experiments, we empirically optimized the *t*~D~ in PGSE and FEXSY sequences separately. However, the use of different *t*~D~ values might be disadvantageous for the global fit analysis, since this yields different ADCs for the two experiments, not allowing constraining *D*~i~ and *D*~e~ to be identical for both datasets. Alternatively, a PGSTE sequence with fixed *t*~D~ and variable *t*~E~ could be used instead of the PGSE sequence to ensure identical *D*~i~ and *D*~e~ values.

To aid visualization in [Fig 6](#pone.0177273.g006){ref-type="fig"}, all the signal intensities for the PGSE experiment were normalized by the maximum *S*~0~ value at the shortest *t*~E~, while for the FEXSY experiment, different normalization constants, *S*~0~(*b*~f~, *t*~m~), were used for different *b*~f~ and *t*~m~ values. The effect of *T*~1~ relaxation was included in the fit parameters *S*~0~ for PGSE and *S*~0~(*b*~f~, *t*~m~) for FEXSY. Although a single value of *T*~1~ was not assumed in our analysis, complementary inversion recovery experiments showed monoexponential signal relaxation with *T*~1~ ≈ 400 ms for all studied samples (data not shown).

The bi-modal attenuations observed in the first row of [Fig 6](#pone.0177273.g006){ref-type="fig"} are due to the diffusion contrast between the intra- and extracellular compartments. For the PGSE experiment, signal intensities decrease with increasing *t*~E~ due to *T*~2~ relaxation, with a more pronounced decrease for the intracellular (high *b* values) compartment than the extracellular compartment (low *b* values). This clearly demonstrates that the intracellular compartment has a shorter *T*~2~ than the extracellular compartment. Varying the diffusion encoding, *b*, allows identifying the relaxation weighted contributions from the two compartments, while varying *t*~E~ at constant diffusion weighting allows quantifying *T*~2~ values for each of the two compartments.

For the FEXSY experiment, different attenuation curves correspond to different mixing times. The data with black markers was measured without the diffusion filter and has the highest signal-to-noise ratio. The equilibrium extracellular signal fraction, $f_{e}^{\text{eq}}$, estimated from these data, as in the standard protocol \[[@pone.0177273.ref019]\], is biased due to the difference in relaxation rate constants. In fact, these data yield the relaxation-affected value $f_{e0,r}^{\text{eq}}$, given by [Eq 22](#pone.0177273.e027){ref-type="disp-formula"}. After diffusion filtering at short *t*~m~, the extracellular signal is significantly more reduced than the intracellular signal and consequently the attenuation curve becomes almost monoexponential. As *t*~m~ increases, the extracellular signal fraction increases through exchange. However, for *aqy1*Δ*agp1*Δ measured at 0°C the exchange is so slow that there is hardly any recovery of the extracellular signal at longer *t*~m~. At 0°C, *aqy1*Δ*agp1*Δ + hAQP1 exhibit an exchange rate constant that is comparable to that in the *aqy1*Δ*agp1*Δ at 20°C. For *aqy1*Δ*agp1*Δ + hAQP1 at 20°C the exchange was too fast to be measured with FEXSY (data not shown).

The reduction of the exchange rate constant with temperature, observed in our experiments, is consistent with the results on baker's yeast presented in Ref. \[[@pone.0177273.ref019]\]. There, the Arrhenius analysis of cell membrane permeability indicates that a significant change in the activation energy occurs at low temperatures. This could be caused by conformational changes in aquaporins or, more likely, due to a phase transition in the membrane lipids, consistent with the permeability decrease of lipid bilayer \[[@pone.0177273.ref056]\].

In our experiments, the temperatures where chosen so that the same set of *t*~m~ in FEXSY allowed a precise quantification of exchange for both *aqy1*Δ*agp1*Δ + hAQP1 and *aqy1*Δ*agp1*Δ yeast strains. The quantitative results for *τ*~i~, *τ*~e~, *T*~2i~ and *T*~2e~ are summarized in [Table 1](#pone.0177273.t001){ref-type="table"}. In *aqy1*Δ*agp1*Δ, which have no aquaporins, the exchange occurs by diffusion through the lipid bilayer or possibly through other transport proteins that might be present in the membrane \[[@pone.0177273.ref057]\]. In this case, the exchange rate constant was reliably quantified at 20°C (see left column of [Fig 6](#pone.0177273.g006){ref-type="fig"}), but it was evidently very slow at 0°C (see middle column of [Fig 6](#pone.0177273.g006){ref-type="fig"}). Although the exchange rate constant could not be quantified in this case, we can conclude that the *τ*~i~ and *τ*~e~ were longer than 1 second. The aquaporin called AQP1 is known to be a very good water transporter \[[@pone.0177273.ref058]\], which is confirmed by our results shown in the right column of [Fig 6](#pone.0177273.g006){ref-type="fig"}. For the cells with aquaporin measured at 20°C, the exchange rate constant was out of range accessible by our method, which is approximately 0.5--20 *s*^−1^, corresponding to *τ*~i~ values approximately in the range of 50--2000 ms. This result clearly confirms a very pronounced difference in the exchange rate constants between *aqy1*Δ*agp1*Δ + hAQP1 and *aqy1*Δ*agp1*Δ. Further interpretation of results in terms of permeability per aquaporin and the discussion about aquaporin properties will be presented elsewhere.

10.1371/journal.pone.0177273.t001

###### Estimated lifetimes and *T*~2~ values for the intra- and extracellular water in different yeast samples.

The lifetimes, *τ*~i~ and *τ*~e~ and relaxation times, *T*~2i~ and *T*~2e~ resulting from a global fit of Eqs [11](#pone.0177273.e012){ref-type="disp-formula"} and [16](#pone.0177273.e017){ref-type="disp-formula"} to the PGSE and FEXSY datasets shown in [Fig 6](#pone.0177273.g006){ref-type="fig"}. Uncertainties correspond to one standard deviation determined by the bootstrapping error analysis.

![](pone.0177273.t001){#pone.0177273.t001g}

  -----------------------------------------------------------------------------------
                        *aqy1*Δ*agp1*Δ\   *aqy1*Δ*agp1*Δ\   *aqy1*Δ*agp1*Δ + hAQP1\
                          *T* = 20°C        *T* = 0°C         *T* = 0°C
  --------------------- ----------------- ----------------- -------------------------
  *τ*~i~ (ms)           260 ± 30          \> 1000           100 ± 20

  *τ*~e~ (ms)           380 ± 60          \> 1000           140 ± 30

  $f_{e}^{\text{eq}}$   0.59 ± 0.04       0.6 ± 0.2         0.58 ± 0.07

  *T*~2i~ (ms)          18.1 ± 0.7        16.1 ± 0.3        16.7 ± 0.4

  *T*~2e~ (ms)          35 ± 1            43.8 ± 0.6        52 ± 3
  -----------------------------------------------------------------------------------

The intra- and extracellular *T*~2~ values shown in [Table 1](#pone.0177273.t001){ref-type="table"} correspond to Δ*R*~2~ values of 22--35 s^−1^. The vertical dotted line in the framed plot of [Fig 5](#pone.0177273.g005){ref-type="fig"} indicates the Δ*R*~2~*t*~E~ = 0.6, which corresponds to Δ*R*~2~ = 35 s^−1^ in our experiments with *t*~E~ = 17 ms. At this Δ*R*~2~, we estimate that the intracellular lifetime would be underestimated by about 20%, if determined from the standard FEXSY protocol ignoring differences in relaxation rate constants. From measurements on 24 cell suspension with different strains of yeast cells (data will be published separately), we see that Δ*R*~2~ values are not specific to a certain yeast strain. In addition, Δ*R*~2~ can vary for different cultivations of the same strain. We speculate that *T*~2~ values are very sensitive to sample preparation. It is clear from our study that an unbiased quantification of exchange requires to properly account for Δ*R*~2~. This is particularly important when comparing different cell cultures, e.g. when transport properties of different aquaporins are of interest.

Our simplified analysis includes two key assumptions: (1) exchange between compartments with restricted and hindered diffusion occurs on a time scale that is slow compared to the diffusion encoding times and (2) relaxation times are slow compared to the gradient pulse duration. The accuracy and applicable range of our analysis should in future studies be thoroughly investigated, e.g. by Monte Carlo simulations, accounting for the effects of restricted diffusion, relaxation and exchange.

Conclusions {#sec013}
===========

Diffusion NMR offers a non-invasive means to quantify molecular exchange in equilibrium-exchange conditions. In several biological applications of interest, where exchange takes place on a relatively slow time-scale of several 100 ms, the FEXSY method is particularly advantageous. An example of such an application involves studying the transport properties of different types of human aquaporins, which can be expressed in genetically modified yeast cells.

In our experiments on genetically modified yeast cells, different *T*~2~ values in the intra- and extracellular compartments were observed. In this case, using the standard FEXSY protocol would yield intracellular lifetimes that are underestimated by approximately 20%. To avoid this significant bias, we combine FEXSY, a method that is specifically sensitive to exchange, with a *D* − *T*~2~ correlation experiment, which is specifically sensitive to transverse relaxation and can detect subtle differences in *T*~2~ values between compartments with different apparent diffusivities. The effects of diffusion, exchange and relaxation can be accounted for in both experiments by using the signal evolution operators for the intra- and extracellular contributions. We analyse the data by global fitting, constraining the lifetimes, *τ*~i~ and *τ*~e~, as well as the relaxation times, *T*~2i~ and *T*~2e~, to be identical for both datasets. We demonstrate this novel method in our proof-of-concept experiments. The results of an extended study involving several types of aquaporins will be presented elsewhere.

When using FEXSY, accounting for differences between intra- and extracellular relaxation times is necessary whenever an unbiased quantification of exchange is required. This is particularly crucial in studies comparing different cell cultures, where it cannot be safely assumed that the values of *T*~2~ or the differences between *T*~2i~ and *T*~2e~ are the same across different cultivations.

The work was financially supported by the Swedish Foundation for Strategic Research (SSF) (AM13-0090), the Cancer foundation (2010/1171 and 2014/575), the Swedish Research Council (VR) (2011-2891, 2014-3910 and the Linnaeus grant), the Organizing Molecular Matter (OMM) center of excellence and CR Development AB (Lund, Sweden). The funder provided support in the form of salaries for authors \[SL\], but did not have any additional role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript. The specific roles of these authors are articulated in the 'author contributions' section.

*aqy1*Δ*agp1*Δ

:   Genetically modified *Pichia pastoris* yeast cells with deletion of the intrinsic aquaporins

*aqy1*Δ*agp1*Δ + hAQP1

:   Genetically modified *Pichia pastoris* yeast cells with deletion of the intrinsic aquaporins and expression of human aquaporin

NMR

:   Nuclear magnetic resonance

PGSE

:   Pulsed gradient spin-echo

PGSTE

:   Pulsed gradient stimulated echo

ADC

:   Apparent diffusion coefficient

DDE

:   Double diffusion encoding

FEXSY

:   Filter exchange spectroscopy

FEXI

:   Filter exchange imaging

[^1]: **Competing Interests:**DT is co-owner of CR Development AB, holding patents related to the described method: Method And System For Magnetic Resonance Imaging, And Use Thereof (WO2008/147326, WO2010/134870). SL is employed by CR Development AB. This does not alter our adherence to PLOS ONE policies on sharing data and materials.

[^2]: **Conceptualization:** SE KE OS KL DT SL.**Data curation:** SE KE KL DT.**Formal analysis:** SE KE SL.**Funding acquisition:** OS KL DT.**Investigation:** SE KE.**Methodology:** SE DT SL.**Project administration:** SE KE KL DT SL.**Resources:** KE KL DT.**Software:** SE SL DT.**Supervision:** KL DT SL.**Validation:** SE KE OS KL DT SL.**Visualization:** SE SL.**Writing -- original draft:** SE SL.**Writing -- review & editing:** SE KE OS KL DT SL.
