Validation of a lanthanide tag for the analysis of protein dynamics by paramagnetic NMR spectroscopy

Paramagnetic NMR restraints for the characterization of protein structural rearrangements

Mobility is a common feature of biomacromolecules, often fundamental for their function. Thus, in many cases, biomacromolecules cannot be described by a single conformation, but rather by a conformational ensemble. NMR paramagnetic data demonstrated quite informative to monitor this conformational variability, especially when used in conjunction with data from different sources. Due to their long-range nature, paramagnetic data can, for instance, i) clearly demonstrate the occurrence of conformational rearrangements, ii) reveal the presence of minor conformational states, sampled only for a short time, iii) indicate the most representative conformations within the conformational ensemble sampled by the molecule, iv) provide an upper limit to the weight of each conformation.

Distinct mechanism of Tb3+ and Eu3+ binding to NCS1

Lanthanides have been frequently used as biomimetic compounds for NMR and fluorescence studies of Ca2+ binding proteins due to having similar physical properties and coordination geometry to Ca2+ ions. Here we report that a member of the neuronal calcium sensor family, neuronal calcium sensor 1, complexes with two lanthanide ions Tb3+ and Eu3+. The affinity for Tb3+ is nearly 50 times higher than that for Ca2+ (Kd,Tb3+ = 0.002 ± 0.0001 μM and Kd, Ca2+ = 91 nM) whereas Eu3+ binding is notably weaker, Kd,Eu3+ = 26 ± 1 μM. Interestingly, despite having identical charge and similar ionic radii, Tb3+ and Eu3+ ions exhibit a distinct binding stoichiometry for NCS1 with one Eu3+ and two Tb3+ ions bound per NCS1 monomer, as demonstrated in fluorescence titration and mass spectrometry studies. These results suggest that the lanthanides' affinity for the individual EF hands is fine-tuned by a small variation in the ion charge density as well as EF hand binding loop amino acid sequence. As observed previously for other lanthanide:protein complexes, the emission intensity of Ln3+ is enhanced upon complexation with the protein, likely due to the displacement of water molecules by oxygen atoms from the coordinating amino acid residues. The overall shape of the Tb3+NCS1 and Eu3+NCS1 monomer shows high levels of similarity compared to the Ca2+ bound protein based on their collision cross section. However, the distinct occupation of EF hands impacts NCS1 oligomerization and affinity for the D2R peptide that mimics the NCS1 binding site on the D2R receptor. Specifically, the Tb3+NCS1 complex populates the dimer and has comparable affinity for the D2R peptide, whereas Eu3+ bound NCS1 remains in the monomeric form with a negligible affinity for the D2R peptide.

Rigidified and Hydrophilic DOTA-like Lanthanoid Ligands: Design, Synthesis, and Dynamic Properties

Limiting the dynamics of paramagnetic tags is crucial for the accuracy of the structural information derived from paramagnetic nuclear magnetic resonance (NMR) experiments. A hydrophilic rigid 2,2',2″,2‴-(1,4,7,10-tetraazacyclododecane-1,4,7,10-tetrayl)tetraacetic acid (DOTA)-like lanthanoid complex was designed and synthesized following a strategy that allows the incorporation of two sets of two adjacent substituents. This resulted in a C₂ symmetric hydrophilic and rigid macrocyclic ring, featuring four chiral hydroxyl-methylene substituents. NMR spectroscopy was used to investigate the conformational dynamics of the novel macrocycle upon complexation with europium and compared to DOTA and its derivatives. The twisted square antiprismatic and square antiprismatic conformers coexist, but the former is favored, which is different from DOTA. Two-dimensional ¹H exchange spectroscopy shows that ring flipping of the cyclen-ring is suppressed due to the presence of the four chiral equatorial hydroxyl-methylene substituents at proximate positions. The reorientation of the pendant arms causes conformational exchange between two conformers. The reorientation of the coordination arms is slower when the ring flipping is suppressed. This indicates that these complexes are suitable scaffolds to develop rigid probes for paramagnetic NMR of proteins. Due to their hydrophilic nature, it is anticipated that they are less likely to cause protein precipitation than their more hydrophobic counterparts.

Paramagnetic Chemical Probes for Studying Biological Macromolecules

Paramagnetic chemical probes have been used in electron paramagnetic resonance (EPR) and nuclear magnetic resonance (NMR) spectroscopy for more than four decades. Recent years witnessed a great increase in the variety of probes for the study of biological macromolecules (proteins, nucleic acids, and oligosaccharides). This Review aims to provide a comprehensive overview of the existing paramagnetic chemical probes, including chemical synthetic approaches, functional properties, and selected applications. Recent developments have seen, in particular, a rapid expansion of the range of lanthanoid probes with anisotropic magnetic susceptibilities for the generation of structural restraints based on residual dipolar couplings and pseudocontact shifts in solution and solid state NMR spectroscopy, mostly for protein studies. Also many new isotropic paramagnetic probes, suitable for NMR measurements of paramagnetic relaxation enhancements, as well as EPR spectroscopic studies (in particular double resonance techniques) have been developed and employed to investigate biological macromolecules. Notwithstanding the large number of reported probes, only few have found broad application and further development of probes for dedicated applications is foreseen.

Phosphoserine for the generation of lanthanide-binding sites on proteins for paramagnetic nuclear magnetic resonance spectroscopy

Pseudocontact shifts (PCSs) generated by paramagnetic lanthanide ions provide valuable long-range structural information in nuclear magnetic resonance (NMR) spectroscopic analyses of biological macromolecules such as proteins, but labelling proteins site-specifically with a single lanthanide ion remains an ongoing challenge, especially for proteins that are not suitable for ligation with cysteine-reactive lanthanide complexes. We show that a specific lanthanide-binding site can be installed on proteins by incorporation of phosphoserine in conjunction with other negatively charged residues, such as aspartate, glutamate or a second phosphoserine residue. The close proximity of the binding sites to the protein backbone leads to good immobilization of the lanthanide ion, as evidenced by the excellent quality of fits between experimental PCSs and PCSs calculated with a single magnetic susceptibility anisotropy (Δ χ ) tensor. An improved two-plasmid system was designed to enhance the yields of proteins with genetically encoded phosphoserine, and good lanthanide ion affinities were obtained when the side chains of the phosphoserine and aspartate residues are not engaged in salt bridges, although the presence of too many negatively charged residues in close proximity can also lead to unfolding of the protein. In view of the quality of the Δ χ tensors that can be obtained from lanthanide-binding sites generated by site-specific incorporation of phosphoserine, this method presents an attractive tool for generating PCSs in stable proteins, particularly as it is independent of cysteine residues.

Paramagnetic NMR in drug discovery

The presence of an unpaired electron in paramagnetic molecules generates significant effects in NMR spectra, which can be exploited to provide restraints complementary to those used in standard structure-calculation protocols. NMR already occupies a central position in drug discovery for its use in fragment screening, structural biology and validation of ligand-target interactions. Paramagnetic restraints provide unique opportunities, for example, for more sensitive screening to identify weaker-binding fragments. A key application of paramagnetic NMR in drug discovery, however, is to provide new structural restraints in cases where crystallography proves intractable. This is particularly important at early stages in drug-discovery programs where crystal structures of weakly-binding fragments are difficult to obtain and crystallization artefacts are probable, but structural information about ligand poses is crucial to guide medicinal chemistry. Numerous applications show the value of paramagnetic restraints to filter computational docking poses and to generate interaction models. Paramagnetic relaxation enhancements (PREs) generate a distance-dependent effect, while pseudo-contact shift (PCS) restraints provide both distance and angular information. Here, we review strategies for introducing paramagnetic centers and discuss examples that illustrate the utility of paramagnetic restraints in drug discovery. Combined with standard approaches, such as chemical shift perturbation and NOE-derived distance information, paramagnetic NMR promises a valuable source of information for many challenging drug-discovery programs.

Design and applications of lanthanide chelating tags for pseudocontact shift NMR spectroscopy with biomacromolecules

In this review, lanthanide chelating tags and their applications to pseudocontact shift NMR spectroscopy as well as analysis of residual dipolar couplings are covered. A complete overview is presented of DOTA-derived and non-DOTA-derived lanthanide chelating tags, critical points in the design of lanthanide chelating tags as appropriate linker moieties, their stability under reductive conditions, e.g., for in-cell applications, the magnitude of the anisotropy transferred from the lanthanide chelating tag to the biomacromolecule under investigation and structural properties, as well as conformational bias of the lanthanide chelating tags are discussed. Furthermore, all DOTA-derived lanthanide chelating tags used for PCS NMR spectroscopy published to date are displayed in tabular form, including their anisotropy parameters, with all employed lanthanide ions, C_B-Ln distances and tagging reaction conditions, i.e., the stoichiometry of lanthanide chelating tags, pH, buffer composition, temperature and reaction time. Additionally, applications of lanthanide chelating tags for pseudocontact shifts and residual dipolar couplings that have been reported for proteins, protein-protein and protein-ligand complexes, carbohydrates, carbohydrate-protein complexes, nucleic acids and nucleic acid-protein complexes are presented and critically reviewed. The vast and impressive range of applications of lanthanide chelating tags to structural investigations of biomacromolecules in solution clearly illustrates the significance of this particular field of research. The extension of the repertoire of lanthanide chelating tags from proteins to nucleic acids holds great promise for the determination of valuable structural parameters and further developments in characterizing intermolecular interactions.

Recent advances in measuring the kinetics of biomolecules by NMR relaxation dispersion spectroscopy

Protein function can be modulated or dictated by the amplitude and timescale of biomolecular motion, therefore it is imperative to study protein dynamics. Nuclear Magnetic Resonance (NMR) spectroscopy is a powerful technique capable of studying timescales of motion that range from those faster than molecular reorientation on the picosecond timescale to those that occur in real-time. Across this entire regime, NMR observables can report on the amplitude of atomic motion, and the kinetics of atomic motion can be ascertained with a wide variety of experimental techniques from real-time to milliseconds and several nanoseconds to picoseconds. Still a four orders of magnitude window between several nanoseconds and tens of microseconds has remained elusive. Here, we highlight new relaxation dispersion NMR techniques that serve to cover this "hidden-time" window up to hundreds of nanoseconds that achieve atomic resolution while studying the molecule under physiological conditions.

How to tackle protein structural data from solution and solid state: An integrated approach

Long-range NMR restraints, such as diamagnetic residual dipolar couplings and paramagnetic data, can be used to determine 3D structures of macromolecules. They are also used to monitor, and potentially to improve, the accuracy of a macromolecular structure in solution by validating or "correcting" a crystal model. Since crystal structures suffer from crystal packing forces they may not be accurate models for the macromolecular structures in solution. However, the presence of real differences should be tested for by simultaneous refinement of the structure using both crystal and solution NMR data. To achieve this, the program REFMAC5 from CCP4 was modified to allow the simultaneous use of X-ray crystallographic and paramagnetic NMR data and/or diamagnetic residual dipolar couplings. Inconsistencies between crystal structures and solution NMR data, if any, may be due either to structural rearrangements occurring on passing from the solution to solid state, or to a greater degree of conformational heterogeneity in solution with respect to the crystal. In the case of multidomain proteins, paramagnetic restraints can provide the correct mutual orientations and positions of domains in solution, as well as information on the conformational variability experienced by the macromolecule.

Determination of pseudocontact shifts of low-populated excited states by NMR chemical exchange saturation transfer

Angular and distance restraints for low populated excited conformations are studied using PCS–CEST NMR spectroscopy.

Information content of long-range NMR data for the characterization of conformational heterogeneity

Long-range NMR data, namely residual dipolar couplings (RDCs) from external alignment and paramagnetic data, are becoming increasingly popular for the characterization of conformational heterogeneity of multidomain biomacromolecules and protein complexes. The question addressed here is how much information is contained in these averaged data. We have analyzed and compared the information content of conformationally averaged RDCs caused by steric alignment and of both RDCs and pseudocontact shifts caused by paramagnetic alignment, and found that, despite the substantial differences, they contain a similar amount of information. Furthermore, using several synthetic tests we find that both sets of data are equally good towards recovering the major state(s) in conformational distributions.

Visualizing transient dark states by NMR spectroscopy

Myriad biological processes proceed through states that defy characterization by conventional atomic-resolution structural biological methods. The invisibility of these 'dark' states can arise from their transient nature, low equilibrium population, large molecular weight, and/or heterogeneity. Although they are invisible, these dark states underlie a range of processes, acting as encounter complexes between proteins and as intermediates in protein folding and aggregation. New methods have made these states accessible to high-resolution analysis by nuclear magnetic resonance (NMR) spectroscopy, as long as the dark state is in dynamic equilibrium with an NMR-visible species. These methods - paramagnetic NMR, relaxation dispersion, saturation transfer, lifetime line broadening, and hydrogen exchange - allow the exploration of otherwise invisible states in exchange with a visible species over a range of timescales, each taking advantage of some unique property of the dark state to amplify its effect on a particular NMR observable. In this review, we introduce these methods and explore two specific techniques - paramagnetic relaxation enhancement and dark state exchange saturation transfer - in greater detail.

A minor conformation of a lanthanide tag on adenylate kinase characterized by paramagnetic relaxation dispersion NMR spectroscopy

NMR relaxation dispersion techniques provide a powerful method to study protein dynamics by characterizing lowly populated conformations that are in dynamic exchange with the major state. Paramagnetic NMR is a versatile tool for investigating the structures and dynamics of proteins. These two techniques were combined here to measure accurate and precise pseudocontact shifts of a lowly populated conformation. This method delivers valuable long-range structural restraints for higher energy conformations of macromolecules in solution. Another advantage of combining pseudocontact shifts with relaxation dispersion is the increase in the amplitude of dispersion profiles. Lowly populated states are often involved in functional processes, such as enzyme catalysis, signaling, and protein/protein interactions. The presented results also unveil a critical problem with the lanthanide tag used to generate paramagnetic relaxation dispersion effects in proteins, namely that the motions of the tag can interfere severely with the observation of protein dynamics. The two-point attached CLaNP-5 lanthanide tag was linked to adenylate kinase. From the paramagnetic relaxation dispersion only motion of the tag is observed. The data can be described accurately by a two-state model in which the protein-attached tag undergoes a 23° tilting motion on a timescale of milliseconds. The work demonstrates the large potential of paramagnetic relaxation dispersion and the challenge to improve current tags to minimize relaxation dispersion from tag movements.

FANTEN: a new web-based interface for the analysis of magnetic anisotropy-induced NMR data

Pseudocontact shifts (PCSs) and residual dipolar couplings (RDCs) arising from the presence of paramagnetic metal ions in proteins as well as RDCs due to partial orientation induced by external orienting media are nowadays routinely measured as a part of the NMR characterization of biologically relevant systems. PCSs and RDCs are becoming more and more popular as restraints (1) to determine and/or refine protein structures in solution, (2) to monitor the extent of conformational heterogeneity in systems composed of rigid domains which can reorient with respect to one another, and (3) to obtain structural information in protein-protein complexes. The use of both PCSs and RDCs proceeds through the determination of the anisotropy tensors which are at the origin of these NMR observables. A new user-friendly web tool, called FANTEN (Finding ANisotropy TENsors), has been developed for the determination of the anisotropy tensors related to PCSs and RDCs and has been made freely available through the WeNMR ( http://fanten-enmr.cerm.unifi.it:8080 ) gateway. The program has many new features not available in other existing programs, among which the possibility of a joint analysis of several sets of PCS and RDC data and the possibility to perform rigid body minimizations.

A bis(pyridine N-oxide) analogue of DOTA: relaxometric properties of the Gd(III) complex and efficient sensitization of visible and NIR-emitting lanthanide(III) cations including Pr(III) and Ho(III)

We report the synthesis of a cyclen-based ligand (4,10-bis[(1-oxidopyridin-2-yl)methyl]-1,4,7,10-tetraazacyclododecane-1,7-diacetic acid=L1) containing two acetate and two 2-methylpyridine N-oxide arms anchored on the nitrogen atoms of the cyclen platform, which has been designed for stable complexation of lanthanide(III) ions in aqueous solution. Relaxometric studies suggest that the thermodynamic stability and kinetic inertness of the Gd(III) complex may be sufficient for biological applications. A detailed structural study of the complexes by (1) H NMR spectroscopy and DFT calculations indicates that they adopt an anti-Δ(λλλλ) conformation in aqueous solution, that is, an anti-square antiprismatic (anti-SAP) isomeric form, as demonstrated by analysis of the (1) H NMR paramagnetic shifts induced by Yb(III) . The water-exchange rate of the Gd(III) complex is ${k{{298\hfill \atop {\rm ex}\hfill}}}$=6.7×10(6)  s(-1) , about a quarter of that for the mono-oxidopyridine analogue, but still about 50 % higher than the ${k{{298\hfill \atop {\rm ex}\hfill}}}$ of GdDOTA (DOTA=1,4,7,10-tetraazacyclododecane-1,4,7,10-tetraacetic acid). The 2-methylpyridine N-oxide chromophores can be used to sensitize a wide range of Ln(III) ions emitting in both the visible (Eu(III) and Tb(III) ) and NIR (Pr(III) , Nd(III) , Ho(III) , Yb(III) ) spectral regions. The emission quantum yield determined for the Yb(III) complex (${Q{{{\rm L}\hfill \atop {\rm Yb}\hfill}}}$=7.3(1)×10(-3) ) is among the highest ever reported for complexes of this metal ion in aqueous solution. The sensitization ability of the ligand, together with the spectroscopic and relaxometric properties of its complexes, constitute a useful step forward on the way to efficient dual probes for optical imaging (OI) and MRI.

Simultaneous use of solution NMR and X-ray data in REFMAC5 for joint refinement/detection of structural differences

The program REFMAC5 from CCP4 was modified to allow the simultaneous use of X-ray crystallographic data and paramagnetic NMR data (pseudocontact shifts and self-orientation residual dipolar couplings) and/or diamagnetic residual dipolar couplings. Incorporation of these long-range NMR restraints in REFMAC5 can reveal differences between solid-state and solution conformations of molecules or, in their absence, can be used together with X-ray crystallographic data for structural refinement. Since NMR and X-ray data are complementary, when a single structure is consistent with both sets of data and still maintains reasonably `ideal' geometries, the reliability of the derived atomic model is expected to increase. The program was tested on five different proteins: the catalytic domain of matrix metalloproteinase 1, GB3, ubiquitin, free calmodulin and calmodulin complexed with a peptide. In some cases the joint refinement produced a single model consistent with both sets of observations, while in other cases it indicated, outside the experimental uncertainty, the presence of different protein conformations in solution and in the solid state.

Transient protein encounters characterized by paramagnetic NMR

Multiple encounters, sampled by interacting proteins in search of the dominant, functionally active orientation, are visualized by paramagnetic NMR.

Interdomain dynamics explored by paramagnetic NMR

An ensemble-based approach is presented to explore the conformational space sampled by a multidomain protein showing moderate interdomain dynamics in terms of translational and rotational motions. The strategy was applied on a complex of calmodulin (CaM) with the IQ-recognition motif from the voltage-gated calcium channel Ca(v)1.2 (IQ), which adopts three different interdomain orientations in the crystal. The N60D mutant of calmodulin was used to collect pseudocontact shifts and paramagnetically induced residual dipolar couplings for six different lanthanide ions. Then, starting from the crystal structure, pools of conformations were generated by free MD. We found the three crystal conformations in solution, but four additional MD-derived conformations had to be included into the ensemble to fulfill all the paramagnetic data and cross-validate optimally against unused paramagnetic data. Alternative approaches led to similar ensembles. Our "ensemble" approach is a simple and efficient tool to probe and describe the interdomain dynamics and represents a general method that can be used to provide a proper ensemble description of multidomain proteins.

Small-molecule binding sites on proteins established by paramagnetic NMR spectroscopy

Determining the three-dimensional structure of a small molecule-protein complex with weak affinity can be a significant challenge. We present a paramagnetic NMR method to determine intermolecular structure restraints based on pseudocontact shifts (PCSs). Since the ligand must be in fast exchange between free and bound states and the fraction bound can be as low as a few percent, the method is ideal for ligands with high micromolar to millimolar dissociation constants. Paramagnetic tags are attached, one at a time, in a well-defined way via two arms at several sites on the protein surface. The ligand PCSs were measured from simple 1D (1)H spectra and used as docking restraints. An independent confirmation of the complex structure was carried out using intermolecular NOEs. The results show that structures derived from these two approaches are similar. The best results are obtained if the magnetic susceptibility tensors of the tags are known, but it is demonstrated that with two-armed probes, the magnetic susceptibility tensor can be predicted with sufficient accuracy to provide a low-resolution model of the ligand orientation and the location of the binding site in the absence of isotope-labeled protein. This approach can facilitate fragment-based drug discovery in obtaining structural information on the initial fragment hits.

Functional dynamics of proteins revealed by solution NMR

Solution NMR spectroscopy can analyze the dynamics of proteins on a wide range of timescales, from picoseconds to even days, in a site-specific manner, and thus its results are complementary to the detailed but largely static structural information obtained by X-ray crystallography. We review recent progresses in a variety of NMR techniques, including relaxation dispersion and paramagnetic relaxation enhancement (PRE), that permit the observation of the low-populated states, which had been 'invisible' with other techniques. In addition, we review how NMR spectroscopy can be used to elucidate functionally relevant protein dynamics.

MaxOcc: a web portal for maximum occurrence analysis

The MaxOcc web portal is presented for the characterization of the conformational heterogeneity of two-domain proteins, through the calculation of the Maximum Occurrence that each protein conformation can have in agreement with experimental data. Whatever the real ensemble of conformations sampled by a protein, the weight of any conformation cannot exceed the calculated corresponding Maximum Occurrence value. The present portal allows users to compute these values using any combination of restraints like pseudocontact shifts, paramagnetism-based residual dipolar couplings, paramagnetic relaxation enhancements and small angle X-ray scattering profiles, given the 3D structure of the two domains as input. MaxOcc is embedded within the NMR grid services of the WeNMR project and is available via the WeNMR gateway at http://py-enmr.cerm.unifi.it/access/index/maxocc . It can be used freely upon registration to the grid with a digital certificate.

Solid-state NMR crystallography through paramagnetic restraints

Pseudocontact shifts (PCSs) measured by solid-state NMR spectroscopy (SS-NMR) on microcrystalline powders of a paramagnetic metalloprotein permit NMR crystallography. Along with other restraints for SS-NMR experiments, the protein molecular structure as well as the correct crystal packing are obtained.

Maximum occurrence analysis of protein conformations for different distributions of paramagnetic metal ions within flexible two-domain proteins

Multidomain proteins are composed of rigid domains connected by (flexible) linkers. Therefore, the domains may experience a large degree of reciprocal reorientation. Pseudocontact shifts and residual dipolar couplings arising from one or more paramagnetic metals successively placed in a single metal binding site in the protein can be used as restraints to assess the degree of mobility of the different domains. They can be used to determine the maximum occurrence (MO) of each possible protein conformation, i.e. the maximum weight that such conformations can have independently of the real structural ensemble, in agreement with the provided restraints. In the case of two-domain proteins, the metal ions can be placed all in the same domain, or distributed between the two domains. It has been demonstrated that the quantity of independent information for the characterization of the system is larger when all metals are bound in the same domain. At the same time, it has been shown that there are practical advantages in placing the metals in different domains. Here, it is shown that distributing the metals between the domains provides a tool for defining a coefficient of compatibility among the restraints obtained from different metals, without a significant decrease of the capability of the MO values to discriminate among conformations with different weights.

Narrowing the conformational space sampled by two-domain proteins with paramagnetic probes in both domains

Calmodulin is a two-domain protein which in solution can adopt a variety of conformations upon reorientation of its domains. The maximum occurrence (MO) of a set of calmodulin conformations that are representative of the overall conformational space possibly sampled by the protein, has been calculated from the paramagnetism-based restraints. These restraints were measured after inclusion of a lanthanide binding tag in the C-terminal domain to supplement the data obtained by substitution of three paramagnetic lanthanide ions to the calcium ion in the second calcium binding loop of the N-terminal domain. The analysis shows that the availability of paramagnetic restraints arising from metal ions placed on both domains, reduces the MO of the conformations to different extents, thereby helping to identify those conformations that can be mostly sampled by the protein.