A Paramagnetic NMR Spectroscopy Toolbox for the Characterisation of Paramagnetic/Spin-Crossover Coordination Complexes and Metal-Organic Cages

The large paramagnetic shifts and short relaxation times resulting from the presence of a paramagnetic centre complicate NMR data acquisition and interpretation in solution. As a result, NMR analysis of paramagnetic complexes is limited in comparison to diamagnetic compounds and often relies on theoretical models. We report a toolbox of 1D (1H, proton-coupled 13C, selective 1H-decoupling 13C, steady-state NOE) and 2D (COSY, NOESY, HMQC) paramagnetic NMR methods that enables unprecedented structural characterisation and in some cases, provides more structural information than would be observable for a diamagnetic analogue. We demonstrate the toolbox's broad versatility for fields from coordination chemistry and spin-crossover complexes to supramolecular chemistry through the characterisation of CoII and high-spin FeII mononuclear complexes as well as a Co4L6 cage.

Measuring transverse relaxation in highly paramagnetic systems

The enhancement of nuclear relaxation rates due to the interaction with a paramagnetic center (known as Paramagnetic Relaxation Enhancement) is a powerful source of structural and dynamics information, widely used in structural biology. However, many signals affected by the hyperfine interaction relax faster than the evolution periods of common NMR experiments and therefore they are broadened beyond detection. This gives rise to a so-called blind sphere around the paramagnetic center, which is a major limitation in the use of PREs. Reducing the blind sphere is extremely important in paramagnetic metalloproteins. The identification, characterization, and proper structural restraining of the first coordination sphere of the metal ion(s) and its immediate neighboring regions is key to understand their biological function. The novel HSQC scheme we propose here, that we termed R2-weighted, HSQC-AP, achieves this aim by detecting signals that escaped detection in a conventional HSQC experiment and provides fully reliable R2 values in the range of 1H R2 rates ca. 50-400 s-1. Independently on the type of paramagnetic center and on the size of the molecule, this experiment decreases the radius of the blind sphere and increases the number of detectable PREs. Here, we report the validation of this approach for the case of PioC, a small protein containing a high potential 4Fe-4S cluster in the reduced [Fe4S4]2+ form. The blind sphere was contracted to a minimal extent, enabling the measurement of R2 rates for the cluster coordinating residues.

Magnetic susceptibility and paramagnetism-based NMR

The magnetic interactions between the nuclear magnetic moment and the magnetic moment of unpaired electron(s) depend on the structure and dynamics of the molecules where the paramagnetic center is located and of their partners. The long-range nature of the magnetic interactions is thus a reporter of invaluable information for structural biology studies, when other techniques often do not provide enough data for the atomic-level characterization of the system. This precious information explains the flourishing of paramagnetism-assisted NMR studies in recent years. Many paramagnetic effects are related to the magnetic susceptibility of the paramagnetic metal. Although these effects have been known for more than half a century, different theoretical models and new approaches have been proposed in the last decade. In this review, we have summarized the consequences for NMR spectroscopy of magnetic interactions between nuclear and electron magnetic moments, and thus of the presence of a magnetic susceptibility due to metals, and we do so using a unified notation.

Using Paramagnetism to Slow Down Nuclear Relaxation in Protein NMR

Paramagnetic metal ions accelerate nuclear spin relaxation; this effect is widely used for distance measurement and called paramagnetic relaxation enhancement (PRE). Theoretical predictions established that, under special circumstances, it is also possible to achieve a reduction in nuclear relaxation rates (negative PRE). This situation would occur if the mechanism of nuclear relaxation in the diamagnetic state is counterbalanced by a paramagnetic relaxation mechanism caused by the metal ion. Here we report the first experimental evidence for such a cross-correlation effect. Using a uniformly ¹⁵N-labeled mutant of calbindin D_9k loaded with either Tm³⁺ or Tb³⁺, reduced R₁ and R₂ relaxation rates of backbone ¹⁵N spins were observed compared with the diamagnetic reference (the same protein loaded with Y³⁺). The effect arises from the compensation of the chemical shift anisotropy tensor by the anisotropic dipolar shielding generated by the unpaired electron spin.

Asymmetry in 13C-13C COSY spectra provides information on ligand geometry in paramagnetic proteins

The relative intensity of Calpha-C' cross-peaks in homonuclear 13C COSY spectra depends on the relaxation properties of Calpha and C' spins, which, in the proximity of a paramagnetic center, are related to the metal-to-carbon distance. Their quantitative analysis has lead, for the cerium-substituted dicalcium protein, calbindin D9k, to the straightforward identification of peaks arising from metal-coordinating groups. The monodentate or bidentate metal binding mode of carboxylates was identified directly via NMR.

NMR Spectroscopy of Paramagnetic Metalloproteins

This article deals with the solution structure determination of paramagnetic metalloproteins by NMR spectroscopy. These proteins were believed not to be suitable for NMR investigations for structure determination until a decade ago, but eventually novel experiments and software protocols were developed, with the aim of making the approach suitable for the goal and as user-friendly and safe as possible. In the article, we also give hints for the optimization of experiments with respect to each particular metal ion, with the aim of also providing a handy tool for nonspecialists. Finally, a section is dedicated to the significant progress made on 13C direct detection, which reduces the negative effects of paramagnetism and may constitute a new chapter in the whole field of NMR spectroscopy.

13C direct detected experiments: optimization for paramagnetic signals

To optimize 13C direct detected experiments for the observation of signals close to a paramagnetic center, we have assessed the sensitivity of different sequences based on CO-Cali coherence transfer. Features of CACO experiments were tested for Calbindin D9k, in which one of the two native Ca2+ ions is replaced by the paramagnetic Ce3+ ion. We have studied the comparison of single vs multiple quantum coherence transfer evolution as well as the influence of in-phase vs anti-phase detection of 13CO signals and finally the comparison of a coherence transfer step based on a CyO in plane with respect to a Cy ali in plane. The acquisition of the anti-phase component of the signal, accomplished by the removal of the last refocusing steps, allowed the identification of some signals unobserved with other pathways. The structural dependency of paramagnetism-induced nuclear relaxation is such that the identification of the most suitable coherence transfer pathway is not known "a priori" but it is driven by the relative proximity of Cali and CO to the paramagnetic center.

New Routes to the Detection of Relaxation Allowed Coherence Transfer in Paramagnetic Molecules

A new pulse sequence is proposed to measure cross correlation rates between 1H Curie spin relaxation and 1H-15N or 1H-13C dipole-dipole coupling (%@mt;sys@%Gamma%@sx@%H,HX%@be@%CS,DD%@sxx@%%@mx@% ) in paramagnetic systems. The new sequence has been conceived to obtain quantitative measurements of cross correlation rates in the close proximity to the paramagnetic center, preventing the loss of information due to fast transverse relaxation. The approach was tested on the dicalcium protein calbindin D9k in which Ca2+ has been substituted at site II with Ce3+. At variance with previously reported experiments, all peaks observed in HSQC experiments tailored to paramagnetic signals give quantitative estimates of %@mt;sys@%Gamma%@sx@%H%@ital@%,%@rsf@%HX%@be@%CS%@ital@%,%@rsf@%DD%@sxx@%%@mx@% . This is crucial to refine the immediate proximity of the metal ion.

Tailored HCCH-TOCSY experiment for resonance assignment in the proximity of a paramagnetic center

The presence of a paramagnetic center may disturb both coherent and incoherent communication between nuclear spins that are affected, to some extent, by the hyperfine interaction. This is a limiting factor to an extensive use of paramagnetic probes in NMR spectroscopy to enhance partial alignment and to exploit cross correlation effects and pseudocontact shifts. We propose here an HCCH-TOCSY experiment tailored to identify spin systems involving resonances that are partly or completely affected by hyperfine interaction. The efficiency of polarization transfer steps when fast relaxing nuclei are involved is discussed. The sequence is tested for the protein Calbindin D(9k), in which one of the two native Ca2+ ions is replaced by the paramagnetic Ce3+ ion as well as for the oxidized form of cytochrome b(562).

Cross-correlation effects involving curie spin relaxation in methyl groups

Cross-correlation effects arising in methyl protons due to the simultaneous presence of dipole-dipole, chemical shift anisotropy, and Curie spin relaxation mechanisms in paramagnetic systems are analyzed. We assess the potential of obtaining structural constraints from the cross-correlation of Curie spin relaxation with dipolar relaxation mechanisms among methyl proton spins. By theoretical analysis and numerical simulations we characterize the transfer functions describing the interconversion processes of different ranks of multispin order. The time dependence of these processes contains a new type of structural information, the orientation of the methyl C(3)-axis with respect to the electron center. Experimental confirmation is found for selected methyl groups in low spin Fe(3+) sperm whale myoglobin.

Cross correlation between the dipole-dipole interaction and the Curie spin relaxation: the effect of anisotropic magnetic susceptibility

Cross-correlated relaxation caused by the interference of nuclear dipole-dipole interaction and the Curie spin relaxation (DD-CSR cross relaxation) is generalized to treat the case of anisotropic magnetic susceptibility, including the important case where the latter originates from zero-field splitting. It is shown that the phenomenon of DD-CSR cross relaxation is absolutely general and to be expected under any electronic configuration. The results of the generalization are presented for a model system, and the consequences for paramagnetic metalloproteins are illustrated with an example of cerium(III)-substituted calbindin. The effects of the magnetic anisotropy are found to be substantial.

Locating the metal ion in calcium-binding proteins by using cerium(III) as a probe

The detection and assignment of NMR spectroscopic signals of carbon atoms from carbonyl and carboxylate groups in the loop hosting the Ce(III) ion was performed for the cerium-substituted calcium-binding protein calbindin D9k. This provided a tool to characterize in solution the first coordination sphere of the metal ion. Due to the well-documented possibility of replacing calcium with metal ions of the Ln(III) series, this approach turns out to be extremely efficient for characterizing in solution the coordination of calcium ions in proteins, independently of the availability of X-ray crystal structures. The present approach completes the structural characterization of lanthanide-substituted calcium-binding proteins, for which the role of long-range constraints arising from hyperfine interaction and self-orientation has already been assessed.

An 1H-13C-13C-edited 1H NMR experiment for making resonance assignments in the active site of heme proteins

In paramagnetic heme proteins, it is often problematic to make proton resonance assignments for heme substituents that do not have large isotropic shifts and consequently lie under the large envelope of polypeptide resonances. Furthermore, assignments that would normally be performed with the aid of HMBC experiments in diamagnetic molecules can prove difficult in the active site of paramagnetic heme proteins if T2(-1) > 2JCH. To circumvent this problem, a new method is presented to selectively detect 1H in 1Hn-13C-13C fragments biosynthetically introduced into the active site of heme proteins. The pulse sequence combines well-known building blocks such as INEPT to transfer 1H spin magnetization to bonded 13C nuclei, followed by INADEQUATE to generate 13C-13C double-quantum coherence that is selected with pulsed field gradients, and finally reverse-INEPT to transfer magnetization back to 1H nuclei for subsequent observation. The new 1Hn-13C-13C edited experiment takes advantage of the relatively large values of 1JCH and 1JCC, avoiding the long interpulse delays in HMBC that compromise the detectability of rapidly relaxing nuclei. The potential applicability of the pulse sequence is demonstrated by its contribution to the unambiguous assignment of the carbonyl carbons in the heme propionates of ferricytochrome b5.

Pseudocontact shifts as constraints for energy minimization and molecular dynamics calculations on solution structures of paramagnetic metalloproteins

The pseudocontact shifts of NMR signals, which arise from the magnetic susceptibility anisotropy of paramagnetic molecules, have been used as structural constraints under the form of a pseudopotential in the SANDER module of the AMBER 4.1 molecular dynamics software package. With this procedure, restrained energy minimization (REM) and restrained molecular dynamics (RMD) calculations can be performed on structural models by using pseudocontact shifts. The structure of the cyanide adduct of the Met80Ala mutant of the yeast iso-1-cytochrome c has been used for successfully testing the calculations. For this protein, a family of structures is available, which was obtained by using NOE and pseudocontact shifts as constraints in a distance geometry program. The structures obtained by REM and RMD calculations with the inclusion of pseudocontact shifts are analyzed.