|Home | About | Journals | Submit | Contact Us | Français|
WaterLOGSY is a popular ligand-observed NMR technique to screen for protein-ligand interactions, yet when applied to measure dissociation constants (KD) through ligand titration, the results were found to be strongly dependent on sample conditions. Herein, we show that accurate KDs can be obtained by waterLOGSY with optimised experimental setup.
Modern drug discovery is often based on finding novel small molecule inhibitors that can bind disease-related proteins and inhibit their biological activities and interactions1,2,3. Understanding how small molecule inhibitors bind their target proteins is therefore one of the major focuses of biochemical and medicinal research, and knowing the strength of this interaction and is an important step in any drug discovery programme as it enables ranking of inhibitors based on their binding affinities. A common way to quantify the strength of protein-ligand interactions is by the dissociation constant (or binding constant, KD).
Nuclear magnetic resonance (NMR) spectroscopy is an established method for the screening of protein ligands and for the quantification of protein-ligand binding4,5,6,7,8,9,10. Broadly, there are two NMR-based protocols to study these binding interactions. The first involves the observation of proteins that are isotopically labelled, either in the backbone amides (e.g. 15N) or with unnatural amino acids (e.g. fluorinated analogues)7,11,12. With these, binding constants can be measured in a ligand titration experiment by following peak intensity changes or chemical shift perturbations of the labelled amino acid resonances4,5,6. A second class of NMR techniques utilises the observation of ligand resonances8,9,10. The ligands of interest are usually small molecules with molecular weight <1kDa enabling the methods to exploit the large differences in size between the protein and ligand. These ligand-observe experiments may be further subdivided into two categories. The first involves observing changes in NMR parameters of the ligand resonances directly, such as relaxation times, linewidths and chemical shifts. The second category involves the observation of ligand binding indirectly through the nuclear Overhauser effect (NOE). The water-ligand observed via gradient spectroscopy (waterLOGSY) experiment belongs to this category.
WaterLOGSY is a method initially designed for the screening of protein ligands from a mixture of potential binders13,14. In the waterLOGSY experiment, compounds that bind to the target protein often give positive resonances (i.e. have the same sign as protein resonances), and compounds that do not interact with the protein usually give negative resonances (Supplementary Fig. S1). The experiment relies on fast dissociation of bound ligands to carry into the free state magnetisation originating from irradiated solvent water that has been transmitted via the protein complex, and so typically works best for weak binding systems, such as those with dissociation constants in the high μM to mM region15.
WaterLOGSY may, in principle, be applied to measure KDs through ligand titration experiments14. The observed waterLOGSY signal of a binding ligand is an average between the NOE that arises for the free population and that which arises for the bound population, which have opposing signs. Corrections to the observed waterLOGSY signals are needed to be made so that the binding isotherm only reflects the bound ligand population. This can be achieved by conducting a control titration that contains only the ligand (i.e. in the absence of the protein; Supplementary Fig. S2). Dalvit et al. demonstrated the use of waterLOGSY to obtain KD with human serum albumin (HSA) and L-tryptophan. By titrating L-tryptophan (50μM to 600μM) into a fixed concentration (10μM) of HSA, a KD of 290μM was obtained14. The measured KD was in agreement with reported binding constants that were obtained by NMR and also other biophysical methods16,17. The use of waterLOGSY to measure KDs has since been applied by several research groups18,19,20.
There is, however, a caveat about the accuracy of such titration experiments as it was discovered that the observed KD values (KDobs) yielded by waterLOGSY may be influenced by experimental conditions21. Fielding et al. demonstrated that KDobs values became greater with increasing protein concentrations by conducting a series of L-tryptophan titration experiments to varying concentrations of bovine serum albumin (BSA) ranging from 15 to 75μM. In this, a KDobs of 230μM was obtained with 15μM BSA, but a ten-fold greater value of 3.0mM was obtained with 75μM BSA. The affinity between a protein and a ligand ought to be independent of the protein concentration, thus the KDobs values obtained by waterLOGSY at higher protein concentrations were overestimated.
A similar influence of experimental conditions on the KDobs values has been reported for saturation transfer difference (STD) NMR22, another NOE-based ligand observe NMR technique. Conceptually STD NMR is similar to waterLOGSY, with magnetisation transferred via NOE from the saturated protein to the bound ligands23. It was observed that KDobss obtained from STD NMR experiments vary significantly with saturation time (i.e. time allowed for the magnetisation to transfer), and accurate KD values could only be determined when the experiments were conducted at a saturation times that were close to zero22. The deviation of KDobs to the ‘true’ KD value was attributed to the rebinding of partially saturated ligands to the protein during the protein irradiation period.
Herein, we report on systematic investigations of the use of waterLOGSY to determine protein-ligand binding constants. Our data show that KDobs values obtained by waterLOGSY titration experiments were indeed strongly dependent on solution conditions, leading to an overestimation of KD values. We further show that accurate binding constants can be determined if the likelihood of ligand rebinding is minimised through use of appropriate experimental conditions.
We hypothesised that the deviation from the ‘true’ KD of KDobs obtained from a waterLOGSY titration experiment originates from a similar mechanism to that reported for STD NMR, namely ligand rebinding during the mixing period22. In this time there exists the probability that a ligand that is already carrying saturation from a binding event may re-enter the binding site before its perturbed magnetisation has fully relaxed back to equilibrium22. Consequently, the observed waterLOGSY signals become reduced in intensity and the resultant binding titration profile does not accurately reflect the 1:1 binding model that was originally assumed14, yielding overestimated dissociation constants. The probability of such fast rebinding processes is increased at higher protein concentrations22, suggesting the greatest perturbation from the ‘true’ KD value may be anticipated for systems with high protein concentrations.
In order to validate this proposal, we repeated the observations made by Fielding et al., choosing the same model system of serum albumin and L-tryptophan as used by Dalvit et al. and Fielding et al.14,21. The reported KD values for L-tryptophan and HSA range from ~100μM to ~250μM16,17,21. As serum albumins are known to have multiple low affinity L-tryptophan binding sites24,25, the standard model that assumes binding to a single site was not used because this does not take into account non-specific binding4,5,6,26. Taira and Terada proposed an alternative model by assuming one high affinity binding site and an unlimited number of low affinity non-specific sites, which they tested with several serum albumin-ligand systems27. We found that this model gave better non-linear curve fitting to our titration data with both HSA and BSA than the standard 1:1 binding model (Supplementary Fig. S3), so was applied throughout this study. The H-δ signal of the bicyclic ring of L-tryptophan was integrated for this study because it is a singlet and gives the highest intensity (Supplementary Fig. S4).
We first measured the KDs of L-tryptophan to HSA at a long mixing time (the period in which magnetisation transfer onto the ligand occurs) of 2.0seconds. Three different HSA concentrations were used and in agreement with previous reports21, we observed that the KDobs increases as the protein concentration increases, ranging from 510μM to 6.5mM (Table 1, Fig. 1 and Supplementary Figs S5–S7). Similar phenomena were observed when STD NMR was applied for direct KD measurements, in which the variation in KDobs was attributed to the rebinding of partially-saturated ligands to the protein during the saturation time period22.
We conducted similar experiments with BSA, observing that at 20μM, the KDobs for L-tryptophan was ~120μM, but this increased to ~1.5mM when the titration was conducted with 75μM BSA (Supplementary Figs S8 and S9). These results confirmed previous observations that the use of waterLOGSY to measure protein-ligand KDs may lead to their overestimation21. Deviation from the true KD value was again most severe at high protein concentrations, which is consistent with our proposal that this inaccuracy originates from rebinding to the protein of partially saturated ligands22.
We also sought to investigate the influence on measured KDs of waterLOGSY mixing times. It was envisaged that at short mixing times there will be less opportunity for partially saturated ligands to rebind, and therefore the KDobs would approach the true KD value. Six different mixing times were used ranging from 0.15seconds to 2.0seconds with a HSA concentration fixed at 100μM. Our results showed that KDobs deviated strongly from the true KD value even with modest mixing times (e.g. <1second). At very short mixing time (0.15seconds), a KDobs of 760μM was obtained, but a linear increase up to 6.5mM was observed as the mixing times became longer (Table 1, Fig. 1 and Supplementary Figs S7, S10–S14). These results show that under conditions that are typically employed in waterLOGSY experiments (mixing time of ~1second), an overestimation of KD may result, so limiting the applicability of waterLOGSY for quantitative ligand binding analyses.
High protein concentrations may encourage the unwanted rebinding of previously saturated ligands, and in accordance with this, we found that the deviation of KDobs with mixing time was less apparent when the protein concentration was reduced. When L-tryptophan titrations were conducted with 50μM HSA and a short 0.5s mixing time, a KDobs of 480μM was obtained (approximately twice the reported KD value), in contrast to 1.5mM observed with 100μM protein at this mixing time. The increase in KDobs as a function of mixing time was also less significant at reduced HSA concentrations (Table 1 and Fig. 1); it was not possible to measure meaningful KDobs at mixing time less than 0.5seconds due to poor signal-to-noise (Supplementary Figs S6, S15–S17).
Most significantly, there was relatively little variation in the KDobs values with mixing times between 0.5 and 1.0second when the L-tryptophan titrations were conducted with 25μM HSA, ranging from 190μM up to 280μM (Table 1, Fig. 1 and Supplementary Figs S5, S18–S20). These results suggest that the derivations of KDobs from the ‘true’ KD value as a function of mixing times will reduce when a lower protein concentration is used. They also suggest that it may be possible to obtain an accurate KD value even with a long mixing time (e.g. 1.0second) if the experiment was conducted with a sufficiently low protein concentration.
To explore this notion further, we selected caffeine as an alternative, weaker binding ligand with a reported KD value for binding to HSA varying between 0.9mM and 1.5mM21,28,29,30. Similar to L-tryptophan, there also exists multiple weak non-specific caffeine binding sites on HSA28,29,30. Titration experiments were again conducted with varying mixing times (0.15seconds to 2.0seconds) and varying protein concentrations (25μM, 50μM and 100μM). When the caffeine titration was conducted with 100μM HSA, even at a very short mixing time (0.15seconds), an erroneous KDobs of 3.2mM was determined, and we again observed a progressive increase in the value with longer mixing times (Table 2, Fig. 2 and Supplementary Figs S21–S26). Similar to L-tryptophan, these deviations of KDobs from the true KD value with different mixing times decreased when the titrations were conducted at lower protein concentrations, as observed in similar work using STD NMR for direct KD measurements22. In fact, when the experiments were performed with 25μM HSA, no variation of KDobs was observed as a function of mixing time and all experiments yielded a KDobs of ~1.5mM, consistent with reported data (Fig. 2 and Supplementary Figs S27–S34). This again indicates that it is possible to obtain an accurate dissociation constant by waterLOGSY titration provided a sufficiently low protein concentration is used, and suggests a means of determining whether the protein concentration used was appropriately ‘low’. Conducting titration experiments using (at least) two mixing times (e.g. 0.5 and 1.0second) should yield similar KD values in the optimised situation, which reflect the true dissociation constant. It is also worth noting that with L-tryptophan at 25μM HSA concentration, meaningful KDobs values could be obtained with a mixing time of up to 1.0second, whilst for caffeine at 25μM HSA concentration, it was possible to obtain meaningful KDobs values even at a very long mixing time (e.g. 2.0seconds). The relationship between KDobs and mixing time is likely complex due to multiple competing factors including (but not limited to) the extent of ligand rebinding and ligand dissociation rates. It is therefore advisable not to conduct quantitative measurements at excessively long mixing times (>1.0sec) to avoid any derivation from the true KD value.
Finally, in order to demonstrate that this protocol is also applicable to 1:1 protein-ligand system in the absence of multiple non-specific binding sites, the binding affinity of a boronic acid with α-chymotrypsin was measured. α-Chymotrypsin is a serine protease known to bind boronic acids as a reversible covalent complex through the formation of a tetrahedral adduct via its nucleophilic serine residue (Supplementary Fig. S35)31. 3-Fluorophenylboronic acid (Supplementary Fig. S35) was chosen as the model ligand system because its binding constant to α-chymotrypsin could be determined independently using 19F NMR spectroscopy.
Thus, by titrating 3-fluorophenylboronic acid to 100μM α-chymotrypsin and monitoring changes in the boronic acid’s 19F chemical shift, a KD of 630μM was determined (Supplementary Fig. S36). WaterLOGSY was then conducted with varying mixing times (0.5 to 2.0seconds) and varying protein concentrations (10μM, 25μM, 50μM and 100μM; Supplementary Fig. S37). Similar to our observations with the serum albumins, deviations in KDobs were observed at high α-chymotrypsin concentrations. For example, millimolar KDobs values were obtained when the titrations were conducted with 50μM and 100μM α-chymotrypsin concentrations (Supplementary Figs S38–S49). However, when the titrations were conducted at low α-chymotrypsin concentration (10μM), the KDobs values obtained at both short and long mixing times were similar (~630, ~780 and ~700μM at 0.5, 1.0 and 2.0s mixing time respectively), and in agreement with the KD value obtained by 19F NMR (630μM), again suggesting that an accurate KD value may be obtained by waterLOGSY when a ‘low’ protein concentration is used.
Taken together, these results confirmed previous observations that measured KD values are influenced directly by the protein concentrations employed in ligand titrations21. We hypothesise that this is likely due to the rebinding of ligands that carry the negative NOE obtained from a previous binding event, as it was proposed for STD NMR measurements, another NOE-based method for the measurement of protein-ligand interactions22. These studies show that the deviation of KDobs from the true KD value is particularly severe at high protein concentrations and long mixing times, in agreement with this proposal. They also showed that it is possible to eliminate this detrimental influence by using a ‘low’ protein concentration, which may be determined by conducting titration experiments with (at least) two different mixing times. Any increase in measured KD with mixing time would suggest the KD values determined may be higher than the true value, at least placing an upper limit on the dissociation constant. There already exist many biophysical techniques, including NMR, able to measure binding affinities, although, most are designed to measure relatively strong protein-ligand binding interactions. If KD determination is necessary for a weak protein-ligand binding system (e.g. during a fragment screen), the choice of biophysical techniques for such a measurement is limited. Whilst the use of waterLOGSY to measure KD may not be required on a regular basis, we believe it fills an important gap in the arsenal of techniques available to cover the whole range of ligand binding affinities and the protocols suggested here should aid in its application, or at least help avoid its inappropriate use.
Unless otherwise stated, all chemicals were from Sigma-Aldrich. Tris-d11 was from Cortecnet. D2O and DMSO-d6 were from Cambridge Isotope Laboratories. HSA and BSA (both fatty acid free and globulin free, ≥99%) and α-chymotrypsin (from bovine pancreas, Type II, ≥40 units/mg protein) were from Sigma-Aldrich.
HSA waterLOGSY experiments were conducted at a 1H frequency of 600MHz using a Bruker Avance spectrometer equipped with a BBI probe. All experiments were conducted at 298K. 5mm diameter NMR tubes with a sample volume of 500μL were used in all experiments. Solutions were buffered using 50mM sodium phosphate (pH 7.5) dissolved in 90% H2O and 10% D2O.
BSA waterLOGSY experiments were conducted at a 1H frequency of 700MHz using a Bruker Avance III spectrometer equipped with a TCI cryoprobe. All experiments were conducted at 293K. 5mm diameter NMR tubes with a sample volume of 500μL were used in all experiments. Solutions were buffered using 50mM Tris-d11 (pH 7.5) and 0.02% NaN3 dissolved in 90% H2O and 10% D2O.
α-Chymotrypsin waterLOGSY experiments were conducted at a 1H frequency of 700MHz using a Bruker Avance III spectrometer equipped with a TCI cryoprobe. All experiments were conducted at 298K. 5mm diameter NMR tubes with a sample volume of 500μL were used in all experiments. Solutions were buffered using 100mM MES buffer (pH 6.5) dissolved in 90% H2O and 10% D2O.
α-Chymotrypsin 19F experiments were conducted at a 19F frequency of 470MHz using a multinuclear BBFO probe. All experiments were conducted at 298K. 5mm diameter NMR tubes with a sample volume of 500μL were used in all experiments. Solutions were buffered using 100mM MES buffer (pH 6.5) dissolved in 90% H2O and 10% D2O. Trifluoroacetone was added to the sample for chemical shift reference at −87ppm vs CFCl3.
WaterLOGSY experiments were conducted using the pulse sequence as described by Dalvit et al.13,14. The pulse tip-angle calibration using the single-pulse nutation method (Bruker pulsecal routine) was undertaken for each sample33. Typical experimental parameters were as follows: relaxation delay 15seconds (BBI probe) or 5seconds (cryoprobe), number of transients 64–256 (BBI probe) or 16 (cryoprobe). Solvent excitation was achieved using a 7.5ms 180 degree selective Gaus1_180r.1000 pulse and water suppression was achieved by the excitation sculpting method using a 2ms 180 degree selective Sinc1.1000 pulse at the H2O frequency, supplemented with a 4ms 90 degree selective Sinc1.1000 water flip-back pulse.
I indicates changes in corrected waterLOGSY intensity from the titrations. Imax indicates the maximum intensity change. [LT] is the titrated ligand concentration and Ns is the non-specific term (where applicable). This uses the approximation [LT] ~ [L] (the free ligand concentration) which is valid when the ligand is used in large excess over the protein. Curve fitting process was conducted using SigmaPlot 12.5 (Systat Software, USA).
How to cite this article: Huang, R. et al. Protein-ligand binding affinity determination by the waterLOGSY method: An optimised approach considering ligand rebinding. Sci. Rep. 7, 43727; doi: 10.1038/srep43727 (2017).
Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
R.H. and I.K.H.L. acknowledge the financial support from Neurological Foundation of New Zealand, Health Research Council of New Zealand, The University of Auckland and Biochemical Society, and instrument support from Dr. M. Schmitz. AB acknowledges Le Conseil Regional d’Auvergne, La Fondation d’Auvergne, and The Erasmus Programme for financial support.
The authors declare no competing financial interests.
Author Contributions R.H. conducted experiments with HSA, I.K.H.L. conducted experiments with BSA and A.B. conducted experiments with α-chymotrypsin. T.D.W.C. and I.K.H.L. designed the studies. T.D.W.C. and I.K.H.L. wrote the manuscript with the help of R.H. All authors reviewed the manuscript.