Open Access Article

This Open Access Article is licensed under a

Creative Commons Attribution 3.0 Unported Licence

Elizaveta A.
Suturina
^{a},
Daniel
Häussinger
^{b},
Kaspar
Zimmermann
^{b},
Luca
Garbuio
^{c},
Maxim
Yulikov
^{c},
Gunnar
Jeschke
^{c} and
Ilya
Kuprov
*^{a}
^{a}School of Chemistry, University of Southampton, Highfield Campus, Southampton, SO17 1BJ, UK. E-mail: i.kuprov@soton.ac.uk
^{b}Department of Chemistry, University of Basel, St. Johanns Ring 19, CH-4056 Basel, Switzerland
^{c}Department of Chemistry and Applied Biosciences, Swiss Federal Institute of Technology in Zurich, Vladimir Prelog Weg 1-5/10, CH-8093 Zürich, Switzerland

Received
19th August 2016
, Accepted 18th January 2017

First published on 20th January 2017

A significant problem with paramagnetic tags attached to proteins and nucleic acids is their conformational mobility. Each tag is statistically distributed within a volume between 5 and 10 Angstroms across; structural biology conclusions from NMR and EPR work are necessarily diluted by this uncertainty. The problem is solved in electron spin resonance, but remains open in the other major branch of paramagnetic resonance – pseudocontact shift (PCS) NMR spectroscopy, where structural biologists have so far been reluctantly using the point paramagnetic centre approximation. Here we describe a new method for extracting probability densities of lanthanide tags from PCS data. The method relies on Tikhonov-regularised 3D reconstruction and opens a new window into biomolecular structure and dynamics because it explores a very different range of conditions from those accessible to double electron resonance work on paramagnetic tags: a room-temperature solution rather than a glass at cryogenic temperatures. The method is illustrated using four different Tm^{3+} DOTA-M8 tagged mutants of human carbonic anhydrase II; the results are in good agreement with rotamer library and DEER data. The wealth of high-quality pseudocontact shift data accumulated by the biological magnetic resonance community over the last 30 years, and so far only processed using point models, could now become a major source of useful information on conformational distributions of paramagnetic tags in biomolecules.

The subject has a long-standing problem – lanthanide-containing protein tags have significant conformational mobility.^{14} Even DOTA-M8,^{15} which uses a sterically overcrowded – and therefore rigid – metal cage,^{16} still has a flexible linker. The conformational mobility of lanthanide tags is visible in the distance distributions measured by double electron resonance,^{17} and in molecular dynamics simulations.^{18} In this situation the commonly used point paramagnetic centre approximation^{3,19} for PCS is not expected to be valid,^{14,20,21} but quantum chemical calculations^{6,7} are prohibitively expensive.

The problem is solved by the recently discovered partial differential equation that treats the probability density of the metal ion and the resulting pseudocontact shift as scalar fields in three dimensions.^{20} We demonstrate here that it may be used to recover the spin label position distribution from the experimental PCS data. This creates a new window into protein structure and dynamics.

(1) |

(2) |

(3) |

The task of recovering ρ(r) from point measurements of σ(r) at the nuclei may be formulated as finding the paramagnetic centre probability density that minimises the following functional:

(4) |

State-of-the art numerical optimisation algorithms^{23} require first and second variations of the error functional with respect to the probability density. The first variations are:

(5) |

(6) |

The trust region reflective Newton–Raphson minimiser, as implemented in the Optimisation Toolbox^{24} supplied with Matlab, was used to obtain the optimum paramagnetic centre probability density on a finite grid, subject to the non-negativity constraint. Numerical implementation details are discussed in our recent paper^{21} and the associated Matlab source code is available in the paramagnetic NMR module supplied with versions 1.8 and later of the Spinach library.^{25}

Pseudocontact shifts were obtained by comparing ^{1}H–^{15}N HSQC spectra of diamagnetic Lu^{3+} and paramagnetic Tm^{3+} DOTA-M8 tagged mutants. PCS assignment was performed in two stages. At the first stage, at least eight of the 26 leucine peaks in the selectively ^{15}N-Leu labelled mutants were identified manually and supplied to NUMBAT,^{28} along with ^{15}N positions from the X-ray structure (PDB:3KS3) of hCA-II.^{29} This enabled the identification of all leucine ^{15}N atoms and the extraction of the approximate point model parameters (metal position and the anisotropic part of the magnetic susceptibility tensor). At the second stage, these parameters were used in another round of NUMBAT calculations on the uniformly ^{15}N labelled mutants to assist in the identification of the rest of the shifted signals, yielding a total of 364 (S217C), 366 (S50C, S220C) and 397 (S166C) unambiguous ^{1}H and ^{15}N PCS assignments. The relevant data is included into the ESI.†

Spin label distance distributions were measured at 10 K using the four-pulse DEER sequence.^{33} All pulses were 12 ns long; the frequency offset between the pump pulse and the detection pulse was 300 MHz. The duration of all DEER traces was at least 3.0 μs – long enough to sample and subtract the intermolecular background. Distance distributions were extracted using the DeerAnalysis package.^{34} The optimal values of the Tikhonov regularization parameter were found using the L-curve method.^{22}

(7) |

Fig. 1 Evolution of the probability density of the Tm^{3+} ion attached to C220 of the S220C mutant of human carbonic anhydrase II with a DOTA-M8 tag during the error functional optimisation process. The initial guess is a uniform distribution within the volume that is at least 2.0 Å from all atoms of the protein itself and at most 12 Å from any of its atoms, corresponding to the region of the space realistically accessible by the Tm^{3+} ion in a tag attached anywhere on the protein surface. As the optimisation proceeds, the probability density gradually becomes zero in the locations that are not consistent with the experimental PCS data. At the end of the optimisation, the probability density is localised, subject to the standard accuracy conditions associated with Tikhonov regularisation,^{22} in the region of space actually accessible to the Tm^{3+} ion. |

Once the approximate location of the paramagnetic centre becomes clear, the refinement of its distribution on a finer grid can proceed with a much reduced variational volume that only involves the region of space immediately adjacent to the approximate location of the tag. An example of such a volume is given in Fig. 2 (red cube in the right panel); a 20 × 20 × 20 Å cube is in practice sufficient.

Fig. 2 Diagnostic information and the outcome of a typical paramagnetic centre probability density reconstruction run. After the initial localisation stage (Fig. 1), the region of space in which the probability density is allowed to vary is chosen (right panel, red cube). Multiple reconstruction runs with different values of the regularisation parameter are performed to obtain the L-curve (left panel, cut-in). The optimum regularisation parameter is extracted as the maximum curvature point on the L-curve (middle panel). The final reconstruction is performed to obtain the probability density (right panel, red cloud) and the fitting plot (left panel). Blue circles in the left panel correspond to the point model fit and the red dots to the probability density fit. |

The optimisation is then performed repeatedly for different values of the regularisation parameter λ in eqn (4). The resulting L-curve is shown as an inset in the left panel of Fig. 2; the corresponding curvature plot is in the middle panel. The optimum value of the regularisation parameter (indicated with a red circle) is calculated and the optimisation is performed again with that value. This yields the paramagnetic centre probability density (red cloud in the right panel) on a fine grid, as well as the plot of the back-calculated pseudocontact shifts against the experimental ones (Fig. 2, left panel). A 256 × 256 × 256 point grid is in practice sufficient; the calculation takes a few hours on a Tesla K40 card.

Once the draft probability density is obtained, a different least squares optimisation is run, this time with respect to the five independent elements of the effective magnetic susceptibility tensor. The tensor is updated and the probability density reconstruction procedure described above is repeated. The whole procedure is performed multiple times until self-consistency is achieved between χ and ρ(r).

The procedure described above relies on two significant assumptions. Firstly, the protein structure is treated as rigid and only the paramagnetic centre is assumed to be delocalised. This is an approximation – in a real protein structure the pseudocontact shifts are also averaged over the distributions in the nuclear positions. From NMR data, in well-defined structures these have position distributions within about 0.4 Å for backbone atoms and 1.0 Å for all heavy atoms.^{42} It is therefore to be expected that the paramagnetic centre distribution obtained from the PCS data would be broader than the real one by approximately that amount. Secondly, eqn (2)–(6) rely on the magnetic susceptibility tensor being the same at each point in the tag distribution. This is not necessarily true because the orientation of the tag can vary. This matter has recently been studied in detail by Shishmarev and Otting;^{14} their conclusion was that a single effective χ tensor can describe the PCS field reasonably well, even in the presence of significant tag mobility. A recent experimental study by Abdelkader et al. has also concluded that using an effective magnetic susceptibility tensor to mask its orientational distribution is a good approximation.^{43} The algebraic structure of eqn (2) suggests that local variations in χ can be compensated by local variations in the probability density – the practical consequences of the constant effective magnetic susceptibility tensor assumption are therefore minor ripples in the probability density. A technical analysis of the accuracy of this approximation is given in the ESI;† the conclusion is that the resulting uncertainty is multiplicative – it would never generate probability density where there was none; it can only scale the true density by a factor related to the norms of the susceptibility tensors involved.

The results of the paramagnetic centre probability density reconstructions for S50C, S166C, S217C and S220C mutants of hCA-II with a Tm^{3+} containing DOTA-M8 tag attached to the corresponding cysteines are presented in Fig. 3 and 4. As could be expected, the paramagnetic centre locations predicted by the point model fits (dark grey three-dimensional crosses in Fig. 3) are located close to the centroids of the probability density distributions computed by regularization (coloured translucent bubbles). The distributions also overlap significantly with the Tm^{3+} ion positions predicted by the rotamer library (swarms of coloured spheres), providing an independent experimental confirmation of the validity of the rotamer library approach.^{40,41}

An important secondary contribution to the chemical shift in paramagnetic systems arises from the residual anisotropic chemical shifts (RACS) that are caused by the weak alignment of the magnetic susceptibility tensor by the applied magnetic field. This effect was recently studied in detail by Otting et al., who estimated the RACS correction magnitude for backbone ^{15}N nuclei to be about 0.1 ppm for a dysprosium ion rigidly coordinated inside a protein structure.^{44} In the context of this work, the RACS correction was assumed to be negligible – a thulium ion (weaker PCS than dysprosium) at the end of a flexible linker would generate a much smaller chemical shift correction than 0.1 ppm, which is itself smaller than the scatter observed in Fig. 2. In situations when RACS are suspected to be significant, we recommend running the reconstruction using only proton PCS data because the effect is negligible for protons.

Fig. 4 presents the comparison between the distance distributions between lanthanide ions in tags at different sites obtained using three physically different methods: DEER,^{34} PCS (this work) and rotamer libraries.^{40,41} For three out of four tagging site pairs, the agreement of the PCS data with the DEER data is very good and significantly better than the agreement of the rotamer library prediction with the DEER data. In the remaining case, both the rotamer library and the PCS prediction deviate from the DEER data by the same amount. This indicates that PCS-based reconstruction of the spatial distribution of the paramagnetic centre performs better than rotamer libraries, although testing on a broader range of proteins would be necessary to make that conclusion in a definitive way. One of the possible explanations for the difference between PCS and DEER reconstructions for the 50–166 dataset is the presence of structural changes caused by the double mutation – individual S50C and S166C mutations (used for PCS) might not have influenced the overall protein geometry in a detectable way, but in the double mutant (used for DEER) the changes could be significant. This conjecture is supported by the fact that a related 166–220 double mutant is completely unstable and could not be expressed in a non-degraded form. The observed difference should not therefore be held against either method; it yields a useful structural insight.

The question of tag probability density reconstruction is particularly pertinent to the many ongoing efforts to characterise domain mobility in proteins.^{45–48} Pseudocontact shift is a convenient parameter for those studies because the timescale of its emergence (i.e. the unpaired electron magnetisation equilibration time in a lanthanide) is in the picoseconds, and the time scale of its observation (i.e. the reciprocal frequency difference between the signals in the paramagnetic NMR spectrum) is in the milliseconds. The former is much faster than protein domain mobility, and the latter is much slower, meaning that the probability density is well defined and simply reflects structural heterogeneity – the corresponding theory is not troubled by the local dynamics effects that make nuclear spin relaxation theory so complicated.

The following would be a reasonable usage scenario for the method described above. In a multi-domain protein or nucleic acid, one of the domains should be tagged with a lanthanide. Pseudocontact shifts measured in the same domain should be used to run a probability density reconstruction. The resulting cloud would be a measure of how rigidly the tag is immobilised relative to its home domain. At the second stage, pseudocontact shifts measured in the other domains should be used to reconstruct the volume that is available to the tag; that volume is an indication of the volume explored by its home domain relative to other domains. The width of the tag distribution in its home domain would then be a measure of the uncertainty in the resulting conformational mobility conclusions. Rigidly immobilised and highly predictable tags^{49} are therefore likely to be beneficial.

On the detailed map of protein mobility analysis methods recently published by Ravera et al.,^{50} the PCS technique described in this paper belongs to the L-curve class, with a significant difference that the penalty functional is not molecular energy (of which there is no notion in the probability density formalism), but the more traditional Laplacian norm.^{22} There exists a possibility of introducing a contrast functional similar to the maximum entropy one,^{20} but we would not recommend using it because it is hard to justify on physical grounds, and also because the distance distribution widths are already in good agreement with other methods (Fig. 4).

Because the extracted distributions have the physical meaning of probability densities, multiple independent datasets (for example, from different metals or different structures in a bundle) may be combined by multiplication. We did not explore this matter further, but it bears notice that the possibility exists.

The probability density reconstruction technique described above is also important because it provides an independent experimental validation for the DEER method – so far, the distributions of the tag at each labelling site could only be modelled, and no experimental technique was available to check the results, except for DEER itself. The good agreement on both the centres and the widths of the distance distributions shown in Fig. 4 is a strong endorsement of the two-electron dipolar spectroscopy results.

- G. N. L. Mar, W. D. Horrocks and R. H. Holm, NMR of Paramagnetic Molecules: Principles and Applications, Elsevier Science, 1973 Search PubMed .
- I. Bertini, C. Luchinat and G. Parigi, Solution NMR of Paramagnetic Molecules: Applications to metallobiomolecules and models, Elsevier Science, 2001 Search PubMed .
- H. M. McConnell, J. Chem. Phys., 1957, 27, 226–229 CrossRef CAS .
- A. D. Buckingham and P. J. Stiles, Mol. Phys., 1972, 24, 99–108 CrossRef CAS .
- R. M. Golding and L. C. Stubbs, J. Magn. Reson., 1980, 40, 115–133 CAS .
- F. Gendron, K. Sharkas and J. Autschbach, J. Phys. Chem. Lett., 2015, 6, 2183–2188 CrossRef CAS PubMed .
- J. Vaara, S. A. Rouf and J. Mares, J. Chem. Theory Comput., 2015, 11, 4840–4849 CrossRef CAS PubMed .
- L. Banci, I. Bertini, G. G. Savellini, A. Romagnoli, P. Turano, M. A. Cremonini, C. Luchinat and H. B. Gray, Proteins, 1997, 29, 68–76 CrossRef CAS .
- L. Banci, I. Bertini, L. K. Bren, A. M. Cremonini, B. H. Gray, C. Luchinat and P. Turano, JBIC, J. Biol. Inorg. Chem., 1996, 1, 117–126 CrossRef CAS .
- L. Banci, I. Bertini, M. A. Cremonini, G. Gori-Savellini, C. Luchinat, K. Wüthrich and P. Güntert, J. Biomol. NMR, 1998, 12, 553–557 CrossRef CAS PubMed .
- G. Otting, Annu. Rev. Biophys., 2010, 39, 387–405 CrossRef CAS PubMed .
- F. Rodriguez-Castañeda, P. Haberz, A. Leonov and C. Griesinger, Magn. Reson. Chem., 2006, 44, S10–S16 CrossRef PubMed .
- X.-C. Su and G. Otting, J. Biomol. NMR, 2010, 46, 101–112 CrossRef CAS PubMed .
- D. Shishmarev and G. Otting, J. Biomol. NMR, 2013, 56, 203–216 CrossRef CAS PubMed .
- D. Häussinger, J.-R. Huang and S. Grzesiek, J. Am. Chem. Soc., 2009, 131, 14761–14767 CrossRef PubMed .
- A. C. L. Opina, M. Strickland, Y.-S. Lee, N. Tjandra, R. A. Byrd, R. E. Swenson and O. Vasalatiy, Dalton Trans., 2016, 4673–4687 RSC .
- H. Yagi, D. Banerjee, B. Graham, T. Huber, D. Goldfarb and G. Otting, J. Am. Chem. Soc., 2011, 133, 10418–10421 CrossRef CAS PubMed .
- A. Potapov, H. Yagi, T. Huber, S. Jergic, N. E. Dixon, G. Otting and D. Goldfarb, J. Am. Chem. Soc., 2010, 132, 9040–9048 CrossRef CAS PubMed .
- I. Bertini, C. Luchinat and G. Parigi, Prog. Nucl. Magn. Reson. Spectrosc., 2002, 40, 249–273 CrossRef CAS .
- G. T. P. Charnock and I. Kuprov, Phys. Chem. Chem. Phys., 2014, 16, 20184–20189 RSC .
- E. Suturina and I. Kuprov, Phys. Chem. Chem. Phys., 2016, 18, 26412–26422 RSC .
- P. C. Hansen, SIAM Rev., 1992, 34, 561–580 CrossRef .
- J. Nocedal and S. Wright, Numerical optimization, Springer Science & Business Media, 2006 Search PubMed .
- MATLAB Optimisation Toolbox, Release 2016a, The MathWorks, Inc., Natick, Massachusetts, United States Search PubMed .
- H. Hogben, M. Krzystyniak, G. Charnock, P. Hore and I. Kuprov, J. Magn. Reson., 2011, 208, 179–194 CrossRef CAS PubMed .
- S. K. Nair, T. L. Calderone, D. W. Christianson and C. Fierke, J. Biol. Chem., 1991, 266, 17320–17325 CAS .
- L. Zheng, U. Baumann and J.-L. Reymond, Nucleic Acids Res., 2004, 32, e115 CrossRef PubMed .
- C. Schmitz, M. J. Stanton-Cook, X.-C. Su, G. Otting and T. Huber, J. Biomol. NMR, 2008, 41, 179–189 CrossRef CAS PubMed .
- B. S. Avvaru, C. U. Kim, K. H. Sippel, S. M. Gruner, M. Agbandje-McKenna, D. N. Silverman and R. McKenna, Biochemistry, 2009, 49, 249–251 CrossRef PubMed .
- I. Gromov, J. Shane, J. Forrer, R. Rakhmatoullin, Y. Rozentzwaig and A. Schweiger, J. Magn. Reson., 2001, 149, 196–203 CrossRef CAS PubMed .
- R. Tschaggelar, B. Kasumaj, M. G. Santangelo, J. Forrer, P. Leger, H. Dube, F. Diederich, J. Harmer, R. Schuhmann and I. García-Rubio, J. Magn. Reson., 2009, 200, 81–87 CrossRef CAS PubMed .
- Y. Polyhach, E. Bordignon, R. Tschaggelar, S. Gandra, A. Godt and G. Jeschke, Phys. Chem. Chem. Phys., 2012, 14, 10762–10773 RSC .
- M. Pannier, S. Veit, A. Godt, G. Jeschke and H. W. Spiess, J. Magn. Reson., 2011, 213, 316–325 CrossRef CAS PubMed .
- G. Jeschke, V. Chechik, P. Ionita, A. Godt, H. Zimmermann, J. Banham, C. Timmel, D. Hilger and H. Jung, Appl. Magn. Reson., 2006, 30, 473–498 CrossRef CAS .
- R. R. Sokal, Univ. Kans. Sci. Bull., 1958, 38, 1409–1438 Search PubMed .
- MATLAB Statistics and Machine Learning Toolbox, Release 2016a, The MathWorks, Inc., Natick, Massachusetts, United States Search PubMed .
- F. Neese, Wiley Interdiscip. Rev.: Comput. Mol. Sci., 2012, 2, 73–78 CrossRef CAS .
- A. K. Rappé, C. J. Casewit, K. Colwell, W. Goddard Iii and W. Skiff, J. Am. Chem. Soc., 1992, 114, 10024–10035 CrossRef .
- L. Boltzmann, Sitzungsber. Akad. Wiss. Wien, Math.-Naturwiss. Kl., Abt. 1, 1872, 66, 275–370 Search PubMed .
- B. Joseph, A. Sikora, E. Bordignon, G. Jeschke, D. S. Cafiso and T. F. Prisner, Angew. Chem., Int. Ed., 2015, 54, 6196–6199 CrossRef CAS PubMed .
- Y. Polyhach, E. Bordignon and G. Jeschke, Phys. Chem. Chem. Phys., 2011, 13, 2356–2366 RSC .
- A. Bhattacharya, R. Tejero and G. T. Montelione, Proteins: Struct., Funct., Bioinf., 2007, 66, 778–795 CrossRef CAS PubMed .
- E. H. Abdelkader, X. Yao, A. Feintuch, L. A. Adams, L. Aurelio, B. Graham, D. Goldfarb and G. Otting, J. Biomol. NMR, 2016, 64, 39–51 CrossRef CAS PubMed .
- M. John, A. Y. Park, G. Pintacuda, N. E. Dixon and G. Otting, J. Am. Chem. Soc., 2005, 127, 17190–17191 CrossRef CAS PubMed .
- J. R. Huth, C. A. Bewley, G. M. Clore, A. M. Gronenborn, B. M. Jackson and A. G. Hinnebusch, Protein Sci., 1997, 6, 2359–2364 CrossRef CAS PubMed .
- M. W. F. Fischer, J. A. Losonczi, J. L. Weaver and J. H. Prestegard, Biochemistry, 1999, 38, 9013–9022 CrossRef CAS PubMed .
- R. W. Montalvao, A. Cavalli, X. Salvatella, T. L. Blundell and M. Vendruscolo, J. Am. Chem. Soc., 2008, 130, 15990–15996 CrossRef CAS PubMed .
- M. Sevvana, V. Vijayan, M. Zweckstetter, S. Reinelt, D. R. Madden, R. Herbst-Irmer, G. M. Sheldrick, M. Bott, C. Griesinger and S. Becker, J. Mol. Biol., 2008, 377, 512–523 CrossRef CAS PubMed .
- P. H. J. Keizers, A. Saragliadis, Y. Hiruma, M. Overhand and M. Ubbink, J. Am. Chem. Soc., 2008, 130, 14802–14812 CrossRef CAS PubMed .
- E. Ravera, L. Sgheri, G. Parigi and C. Luchinat, Phys. Chem. Chem. Phys., 2016, 18, 5686–5701 RSC .
- M. Allegrozzi, I. Bertini, M. B. L. Janik, Y. M. Lee, G. H. Lin and C. Luchinat, J. Am. Chem. Soc., 2000, 122, 4154–4161 CrossRef CAS .
- I. Bertini, C. Luchinat, G. Parigi and R. Pierattelli, ChemBioChem, 2005, 6, 1536–1549 CrossRef CAS PubMed .

## Footnote |

† Electronic supplementary information (ESI) available. See DOI: 10.1039/c6sc03736d |

This journal is © The Royal Society of Chemistry 2017 |