Lanthanide-Induced Pseudocontact Shifts for Solution Structure Refinements of Macromolecules in Shells up to 40 Å from the Metal Ion

One touch is all it takes: the supramolecular interaction between ubiquitin and lanthanide complexes revisited by paramagnetic NMR and molecular dynamics

The supramolecular interaction between lanthanide complexes and proteins is at the heart of numerous chemical and biological studies. Some of these complexes have demonstrated remarkable interaction properties with proteins or peptides in solution and in the crystalline state. Here we have used the paramagnetism of lanthanide ions to characterize the affinity of two lanthanide complexes for ubiquitin. As the interaction process is dynamic, the acquired NMR data only reflect the time average of the different steps. We have used molecular dynamics (MD) simulations to get a deeper insight into the detailed interaction scenario at the microsecond scale. This NMR/MD approach enabled us to establish that the tris-dipicolinate complex interacts specifically with arginines and lysines, while the crystallophore explores the protein surface through weak interactions with carboxylates. These observations shed new light on the dynamic interaction properties of these complexes, which will ultimately enable us to propose a crystallization mechanism.

Electrostatic interactions - key determinants of the metal selectivity in La3+ and Ca2+ binding proteins

Nearly half of all known proteins contain metal co-factors. In the course of evolution two dozen metal cations (mostly monovalent and divalent species) have been selected to participate in processes of vital importance for living organisms. Trivalent metal cations have also been selected, although to a lesser extent as compared with their mono- and divalent counterparts. Notably, factors governing the metal selectivity in trivalent metal centers in proteins are less well understood than those in the respective divalent metal centers. Thus, the source of high La³⁺/Ca²⁺ selectivity in lanthanum-binding proteins, as compared with that of calcium-binding proteins (i.e., calmodulin), is still shrouded in mystery. The well-calibrated thermochemical calculations, performed here, reveal the dominating role of electrostatic interactions in shaping the metal selectivity in La³⁺-binding centers. The calculations also disclose other (second-order) determinants of metal selectivity in these systems, such as the rigidity and extent of solvent exposure of the binding site. All these factors are also implicated in shaping the metal selectivity in Ca²⁺-binding proteins.

Site-Selective Solvation-Induced Conformational Switching of Heteroleptic Heteronuclear Tb(III) and Y(III) Trisphthalocyaninates for the Control of Their Magnetic Anisotropy

In the present work, we report the synthesis of isomeric heteronuclear terbium(III) and yttrium(III) triple-decker phthalocyaninates [(BuO)₈Pc]M[(BuO)₈Pc]M*[(15C5)₄Pc] (M = Tb, M* = Y or M = Y, M* = Tb, [(BuO)₈Pc]^2--octa-n-butoxyphthalocyaninato-ligand, [(15C5)₄Pc]^2--tetra-15-crown-5-phthalocyaninato-ligand). We show that these complexes undergo solvation-induced switching: the conformers in which both metal centers are in square-antiprismatic environments are stabilized in toluene, whereas in dichloromethane, the metal centers M and M* are in distorted prismatic and antiprismatic environments, respectively. This conclusion follows from the detailed analysis of lanthanide-induced shifts in ¹H NMR spectra, which makes it possible to extract the axial component of the magnetic susceptibility tensor χaxTb and to show that this term is particularly sensitive to conformational switching when terbium(III) ion is placed in the switchable "M" site. This result provides a new tool for controlling the magnetic properties of lanthanide complexes with phthalocyanine ligands.

Recent Advances in NMR Protein Structure Prediction with ROSETTA

Nuclear magnetic resonance (NMR) spectroscopy is a powerful method for studying the structure and dynamics of proteins in their native state. For high-resolution NMR structure determination, the collection of a rich restraint dataset is necessary. This can be difficult to achieve for proteins with high molecular weight or a complex architecture. Computational modeling techniques can complement sparse NMR datasets (<1 restraint per residue) with additional structural information to elucidate protein structures in these difficult cases. The Rosetta software for protein structure modeling and design is used by structural biologists for structure determination tasks in which limited experimental data is available. This review gives an overview of the computational protocols available in the Rosetta framework for modeling protein structures from NMR data. We explain the computational algorithms used for the integration of different NMR data types in Rosetta. We also highlight new developments, including modeling tools for data from paramagnetic NMR and hydrogen-deuterium exchange, as well as chemical shifts in CS-Rosetta. Furthermore, strategies are discussed to complement and improve structure predictions made by the current state-of-the-art AlphaFold2 program using NMR-guided Rosetta modeling.

Interdomain Dynamics via Paramagnetic NMR on the Highly Flexible Complex Calmodulin/Munc13-1

Paramagnetic NMR constraints are very useful to study protein interdomain motion, but their interpretation is not always straightforward. On the example of the particularly flexible complex Calmodulin/Munc13-1, we present a new approach to characterize this motion with pseudocontact shifts and residual dipolar couplings. Using molecular mechanics, we sampled the conformational space of the complex and used a genetic algorithm to find ensembles that are in agreement with the data. We used the Bayesian information criterion to determine the ideal ensemble size. This way, we were able to make an accurate, unambiguous, reproducible model of the interdomain motion of Calmodulin/Munc13-1 without prior knowledge about the domain orientation from crystallography.

Distinct stereospecific effect of chiral tether between a tag and protein on the rigidity of paramagnetic tag

Flexibility between the paramagnetic tag and its protein conjugates is a common yet unresolved issue in the applications of paramagnetic NMR spectroscopy in biological systems. The flexibility greatly attenuates the magnetic anisotropy and compromises paramagnetic effects especially for pseudocontact shift and residual dipolar couplings. Great efforts have been made to improve the rigidity of paramagnetic tag in the protein conjugates, however, the effect of local environment vicinal to the protein ligation site on the paramagnetic effects remains poorly understood. In the present work, the stereospecific effect of chiral tether between the protein and a tag on the paramagnetic effects produced by the tag attached via a D- and L-type linker between the protein and paramagnetic metal chelating moiety was assessed. The remarkable chiral effect of the D- and L-type tether between the tag and the protein on the rigidity of paramagnetic tag is disclosed in a number of protein-tag-Ln complexes. The chiral tether formed between the D-type tag and L-type protein surface minimizes the effect of the local environment surrounding the ligation site on the averaging of paramagnetic tag, which is helpful to preserve the rigidity of a paramagnetic tag in the protein conjugates.

Electron correlation and vibrational effects in predictions of paramagnetic NMR shifts

Electronic structure calculations are fundamentally important for the interpretation of nuclear magnetic resonance (NMR) spectra from paramagnetic systems that include organometallic and inorganic compounds, catalysts, or metal-binding sites in proteins. Prediction of induced paramagnetic NMR shifts requires knowledge of electron paramagnetic resonance (EPR) parameters: the electronic g tensor, zero-field splitting D tensor, and hyperfine A tensor. The isotropic part of A, called the hyperfine coupling constant (HFCC), is one of the most troublesome properties for quantum chemistry calculations. Yet, even relatively small errors in calculations of HFCC tend to propagate into large errors in the predicted NMR shifts. The poor quality of A tensors that are currently calculated using density functional theory (DFT) constitutes a bottleneck in improving the reliability of interpretation of the NMR spectra from paramagnetic systems. In this work, electron correlation effects in calculations of HFCCs with a hierarchy of ab initio methods were assessed, and the applicability of different levels of DFT approximations and the coupled cluster singles and doubles (CCSD) method was tested. These assessments were performed for the set of selected test systems comprising an organic radical, and complexes with transition metal and rare-earth ions, for which experimental data are available. Severe deficiencies of DFT were revealed but the CCSD method was able to deliver good agreement with experimental data for all systems considered, however, at substantial computational costs. We proposed a more computationally tractable alternative, where the A was computed with the coupled cluster theory exploiting locality of electron correlation. This alternative is based on the domain-based local pair natural orbital coupled cluster singles and doubles (DLPNO-CCSD) method. In this way the robustness and reliability of the coupled cluster theory were incorporated into the modern formalism for the prediction of induced paramagnetic NMR shifts, and became applicable to systems of chemical interest. This approach was verified for the bis(cyclopentadienyl)vanadium(II) complex (Cp₂V; vanadocene), and the metal-binding site of the Zn²⁺ → Co²⁺ substituted superoxide dismutase (SOD) metalloprotein. Excellent agreement with experimental NMR shifts was achieved, which represented a substantial improvement over previous theoretical attempts. The effects of vibrational corrections to orbital shielding and hyperfine tensor were evaluated and discussed within the second-order vibrational perturbation theory (VPT2) framework.

Pseudocontact Shifts in Biomolecular NMR Spectroscopy

Paramagnetic centers in biomolecules, such as specific metal ions that are bound to a protein, affect the nuclei in their surrounding in various ways. One of these effects is the pseudocontact shift (PCS), which leads to strong chemical shift perturbations of nuclear spins, with a remarkably long range of 50 Å and beyond. The PCS in solution NMR is an effect originating from the anisotropic part of the dipole-dipole interaction between the magnetic momentum of unpaired electrons and nuclear spins. The PCS contains spatial information that can be exploited in multiple ways to characterize structure, function, and dynamics of biomacromolecules. It can be used to refine structures, magnify effects of dynamics, help resonance assignments, allows for an intermolecular positioning system, and gives structural information in sensitivity-limited situations where all other methods fail. Here, we review applications of the PCS in biomolecular solution NMR spectroscopy, starting from early works on natural metalloproteins, following the development of non-natural tags to chelate and attach lanthanoid ions to any biomolecular target to advanced applications on large biomolecular complexes and inside living cells. We thus hope to not only highlight past applications but also shed light on the tremendous potential the PCS has in structural biology.

Lanthanide complexes with polyaminopolycarboxylates as prospective NMR/MRI diagnostic probes: peculiarities of molecular structure, dynamics and paramagnetic properties

The paramagnetic lanthanide complexes with polyaminopolycarboxylate (PAPC) ligands attract considerable attention from the standpoint of potential applications thereof as relaxation agents used in medical magnetic resonance imaging (MRI) and in luminescent materials, as well as owing to promising use thereof as paramagnetic labels for studying the properties of biopolymers since they exhibit thermodynamic stability, good solubility in aqueous media and moderate toxicity. For the last decades, the NMR methods have been used to determine the physical and chemical properties of paramagnetic Ln compounds. The studies concerning paramagnetic NMR lanthanide-induced shifts (LISs) in dissolved Ln complexes, as well as the analysis of band shape as a function of temperature make it possible to obtain valuable information on the structure, intra- and intermolecular dynamics and paramagnetic properties thereof. This review is devoted solely to the following features: firstly, the processes of intramolecular dynamics of lanthanide complexes with polyamino-polycarboxylate ligands such as DOTA, EDTA and DTPA and their derivatives studied by NMR; secondly, the LISs of lanthanide complexes with EDTA, DOTA, DTPA and some of their derivatives depending on temperature and pH. Moreover, in this review, for the first time, the dependence of the activation energy of molecular dynamics in complexes with polydentate ligands on the atomic number of the lanthanide cation is analyzed and a monotonic change in energy is detected, which is due to the effect of lanthanide contraction. It should be noted that this phenomenon is quite general and may also appear in the future in many other series of lanthanide complexes with both other multidentate ligands and with bidentate and monodentate ligands. In the future, it is possible to predict the dependence of the properties of certain lanthanide complexes on the ionic radius of the lanthanide cation based on the approaches presented in the review. In this review, we have also presented the dynamic NMR as the main research method widely used to analyze the processes of molecular dynamics, and the structural studies based on the NMR relaxation spectroscopy and LIS analysis.

A Chiral Lanthanide Tag for Stable and Rigid Attachment to Single Cysteine Residues in Proteins for NMR, EPR and Time-Resolved Luminescence Studies

A lanthanide-binding tag site-specifically attached to a protein presents a tool to probe the protein by multiple spectroscopic techniques, including nuclear magnetic resonance, electron paramagnetic resonance and time-resolved luminescence spectroscopy. Here a new stable chiral LnIII tag, referred to as C12, is presented for spontaneous and quantitative reaction with a cysteine residue to generate a stable thioether bond. The synthetic protocol of the tag is relatively straightforward, and the tag is stable for storage and shipping. It displays greatly enhanced reactivity towards selenocysteine, opening a route towards selective tagging of selenocysteine in proteins containing cysteine residues. Loaded with TbIII or TmIII ions, the C12 tag readily generates pseudocontact shifts (PCS) in protein NMR spectra. It produces a relatively rigid tether between lanthanide and protein, which is beneficial for interpretation of the PCSs by single magnetic susceptibility anisotropy tensors, and it is suitable for measuring distance distributions in double electron-electron resonance experiments. Upon reaction with cysteine or other thiol compounds, the TbIII complex exhibits a 100-fold enhancement in luminescence quantum yield, affording a highly sensitive turn-on luminescence probe for time-resolved FRET assays and enzyme reaction monitoring.

Capturing the dynamic association between a tris-dipicolinate lanthanide complex and a decapeptide: a combined paramagnetic NMR and molecular dynamics exploration

In the realm of biomolecules, peptides can present a large diversity of structures. Our study sheds new light on the structural interplay between a tris-dipicolinate lanthanide probe and a decapeptide SASYKTLPRG. Although a rather trivial, electrostatically driven interaction was expected, the combination of paramagnetic NMR and molecular dynamics simulations reveals a highly dynamic association process and allows for providing extensive insights into the interaction sites and their occupancy. This study highlights the importance of a large conformational sampling to reconcile characteristic time in NMR with molecular dynamics simulations, where sampling in the microsecond range is needed. This study opens the door for a detailed mechanistic elucidation of the early steps of lanthanide complex-peptide or lanthanide complex-protein interaction or self-assembly processes.

CeTiO2N oxynitride perovskite: paramagnetic 14N MAS NMR without paramagnetic shifts

14N magic-angle spinning (MAS) nuclear magnetic resonance (NMR) spectra of diamagnetic LaTiO2N perovskite oxynitride and its paramagnetic counterpart CeTiO2N are presented. The latter, to the best of our knowledge, constitutes the first high-resolution 14N MAS NMR spectrum collected from a paramagnetic solid material. The unpaired 4f-electrons in CeTiO2N do not induce a paramagnetic 14N NMR shift. This is remarkable given the direct Ce−N contacts in the structure for which ab initio calculations predict substantial Ce→14N contact shift interaction. The same effect is revealed with 14N MAS NMR for SrWO2N (unpaired 5d-electrons).

Characterization of lanthanoid-binding proteins using NMR spectroscopy

The variety of magnetic properties exhibited by paramagnetic lanthanoids provides outstanding information in NMR-based structural biology and therefore can be a very useful tool for characterizing lanthanoid-binding proteins. Because of their dependence on the relative positions of the protein nuclei and of the lanthanoid ion, the paramagnetic restraints (PCS, PRDC and PRE) provide information on structure and dynamics of proteins. In this Chapter, we cover the use of lanthanoids in structural biology including protein sample preparation, NMR experiments and data interpretation.

Paramagnetic NMR Spectroscopy Is a Tool to Address Reactivity, Structure, and Protein–Protein Interactions of Metalloproteins: The Case of Iron–Sulfur Proteins

The study of cellular machineries responsible for the iron–sulfur (Fe–S) cluster biogenesis has led to the identification of a large number of proteins, whose importance for life is documented by an increasing number of diseases linked to them. The labile nature of Fe–S clusters and the transient protein–protein interactions, occurring during the various steps of the maturation process, make their structural characterization in solution particularly difficult. Paramagnetic nuclear magnetic resonance (NMR) has been used for decades to characterize chemical composition, magnetic coupling, and the electronic structure of Fe–S clusters in proteins; it represents, therefore, a powerful tool to study the protein–protein interaction networks of proteins involving into iron–sulfur cluster biogenesis. The optimization of the various NMR experiments with respect to the hyperfine interaction will be summarized here in the form of a protocol; recently developed experiments for measuring longitudinal and transverse nuclear relaxation rates in highly paramagnetic systems will be also reviewed. Finally, we will address the use of extrinsic paramagnetic centers covalently bound to diamagnetic proteins, which contributed over the last twenty years to promote the applications of paramagnetic NMR well beyond the structural biology of metalloproteins.

Lanthanide Tagging of Oligonucleotides to Nucleobase for Paramagnetic NMR

Although lanthanide tags, which have large anisotropic magnetic susceptibilities, have already been introduced to enrich NMR parameters by long‐range pseudoconact shifts (PCSs) and residual dipolar couplings (RDCs) of proteins, their application to nucleotides has so far been limited to one previous report, due to the high affinities of lanthanides for the phosphodiester backbone of nucleotides and difficult organic synthesis. Herein, we report successful attachment of a lanthanide tag to a chemically synthesized oligonucleotide via a disulfide bond. NMR experiments reveal PCSs of up to 1 ppm and H−H RDCs of up to 8 Hz at 950 MHz. Although weaker magnetic alignment was achieved than with proteins, the paramagnetic data could be fitted to the known structure of the DNA, taking the mobility of the tag into account. While further rigidification of the tag is desirable, this tag could also be used to measure heteronuclear RDCs of ¹³C,¹⁵N‐labeled chemically synthesized DNA and RNA.

Cysteine oxidation to the sulfinic acid induces oxoform-specific lanthanide binding and fluorescence in a designed peptide

Cysteine sulfinic acid (Cys-SO₂^-) is a protein post-translational modification that is formed reversibly under oxidative conditions. A short, encodable peptide was developed whose metal binding and terbium luminescence are dependent on cysteine (Cys) oxidation to the sulfinic acid. The protein design is based on the modification of a key metal-binding aspartate (Asp) in a canonical EF-Hand motif (DKDADGWISPAEAK) to Cys. In this design, Cys in the thiol oxidation state does not mimic the native Asp, and thus the peptide binds terbium(III) (Tb³⁺) poorly and exhibits weak terbium luminescence (fluorescence). In contrast, when Cys is oxidized to the Cys sulfinic acid oxoform, the Cys sulfinate effectively mimics Asp, resulting in a significant increase in terbium affinity and luminescence. Asp residues at positions 1, 3, and 5 of the EF-Hand motif were examined as potential sites for Cys oxidation-responsive metal binding. The peptide with Cys at residue 1 exhibited the highest Tb³⁺ affinity in both oxidation states. The peptide with Cys at residue 3 exhibited a 4.2-fold distinction in affinity between the oxidation states. Most significantly, the peptide with Cys at residue 5 had only modest Tb³⁺ affinity as the Cys thiol, but exhibited a 30-fold increase in Tb³⁺ affinity and an 18-fold increase in Tb³⁺ luminescence on Cys oxidation to the sulfinic acid. This peptide (Ac-DKDACGWISPAEAK-NH₂) exhibited selective Tb³⁺ binding via Cys-SO₂^- over the thiol, S-glutathionyl, S-nitrosyl, and sulfonic acid oxoforms, indicating substantially greater Lewis basicity of the sulfinate than the sulfonate. NMR spectroscopy and quantum homology modeling indicated that the designed peptide binds metal with an overall geometry similar to that of an EF-Hand motif, with the Cys sulfinate effectively replacing Asp as a metal-binding ligand. This peptide was applied to detect Cys oxidation to the sulfinic acid by fluorescence spectroscopy, suggesting its broader application in understanding Cys sulfinic acid biology.

Paramagnetic solid-state NMR of proteins

The paramagnetic properties of metal ions and stable radicals can affect NMR spectra, which can lead to changes in peak intensities, relaxation times and chemical shifts. The changes from paramagnetic effects provide intriguing opportunities for solid-state NMR studies of proteins. In this review, we summarized the trends and progress of paramagnetic solid-state NMR of proteins in the past decade, and showed that paramagnetic effects have great potential applications for sensitivity enhancement, structure determination and topological analysis for microcrystalline proteins, protein complexes, protein aggregates and membrane proteins.

Three-Dimensional Protein Structure Determination Using Pseudocontact Shifts of Backbone Amide Protons Generated by Double-Histidine Co2+-Binding Motifs at Multiple Sites

Pseudocontact shifts (PCSs) generated by paramagnetic metal ions contribute highly informative long-range structure restraints that can be measured in solution and are ideally suited to guide structure prediction algorithms in determining global protein folds. We recently demonstrated that PCSs, which are relatively small but of high quality, can be generated by a double-histidine (dHis) motif in an α-helix, which provides a well-defined binding site for a single Co²⁺ ion. Here we show that PCSs of backbone amide protons generated by dHis-Co²⁺ motifs positioned in four different α-helices of a protein deliver excellent restraints to determine the three-dimensional (3D) structure of a protein in a way akin to the global positioning system (GPS). We demonstrate the approach with GPS-Rosetta calculations of the 3D structure of the C-terminal domain of the chaperone ERp29 (ERp29-C). Despite the relatively small size of the PCSs generated by the dHis-Co²⁺ motifs, the structure calculations converged readily. Generating PCSs by the dHis-Co²⁺ motif thus presents an excellent alternative to the use of lanthanide tags.

Design of a Protein Motif Responsive to Tyrosine Nitration and an Encoded Turn-Off Sensor of Tyrosine Nitration

Tyrosine nitration is a protein post-translational modification that is predominantly non-enzymatic and is observed to be increased under conditions of nitrosative stress and in numerous disease states. A small protein motif (14-18 amino acids) responsive to tyrosine nitration has been developed. In this design, nitrotyrosine replaced the conserved Glu12 of an EF-hand metal-binding motif. Thus, the non-nitrated peptide bound terbium weakly. In contrast, tyrosine nitration resulted in a 45-fold increase in terbium affinity. Nuclear magnetic resonance spectroscopy indicated direct binding of nitrotyrosine to the metal and EF-hand-like metal contacts in this designed peptide. Nitrotyrosine is an efficient quencher of fluorescence. To develop a sensor of tyrosine nitration, the initial design was modified to incorporate Glu residues at EF-hand positions 9 and 16 as additional metal-binding residues, to increase the terbium affinity of the peptide with unmodified tyrosine. This peptide with a tyrosine at residue 12 bound terbium and effectively sensitized terbium luminescence. Tyrosine nitration resulted in a 180-fold increase in terbium affinity ( Kd = 1.6 μM) and quenching of terbium luminescence. This sequence was incorporated as an encoded protein tag and applied as a turn-off fluorescent protein sensor of tyrosine nitration. The sensor was responsive to nitration by peroxynitrite, with fluorescence quenched upon nitration. The greater terbium affinity upon tyrosine nitration resulted in a large dynamic range and sensitivity to substoichiometric nitration. An improved approach for the synthesis of peptides containing nitrotyrosine was also developed, via the in situ silyl protection of nitrotyrosine. This work represents the first designed, encodable protein motif that is responsive to tyrosine nitration.

Effects of chelator lipids, paramagnetic metal ions and trehalose on liposomes by solid-state NMR

The effects of various lipid bound paramagnetic metal ions on liposomes prepared in the presence of trehalose and chelator lipids are evaluated to observe site-specific signal changes on liposome samples with optimal resolution in solid-state NMR spectroscopy. We found that Mn²⁺, Gd³⁺ and Dy³⁺ have different influences on the lipid ¹³C sites depending on their penetration depths into the bilayer, which can be extracted as distance information. The trehalose-liposome mixture is efficiently packed into solid-state NMR rotors and provides optimal resolution at reasonable instrument temperatures (10-50 °C). The effectiveness and convenience of the trehalose preparation for studying a membrane protein in liposomes are demonstrated by a membrane sample with a model membrane peptide to show that trehalose is useful to prepare consistent and stable membrane protein liposome samples for solid-state NMR.

Structure and dynamics of Helicobacter pylori nickel-chaperone HypA: an integrated approach using NMR spectroscopy, functional assays and computational tools

Helicobacter pylori HypA (HpHypA) is a metallochaperone necessary for maturation of [Ni,Fe]-hydrogenase and urease, the enzymes required for colonization and survival of H. pylori in the gastric mucosa. HpHypA contains a structural Zn(II) site and a unique Ni(II) binding site at the N-terminus. X-ray absorption spectra suggested that the Zn(II) coordination depends on pH and on the presence of Ni(II). This study was performed to investigate the structural properties of HpHypA as a function of pH and Ni(II) binding, using NMR spectroscopy combined with DFT and molecular dynamics calculations. The solution structure of apo,Zn-HpHypA, containing Zn(II) but devoid of Ni(II), was determined using 2D, 3D and 4D NMR spectroscopy. The structure suggests that a Ni-binding and a Zn-binding domain, joined through a short linker, could undergo mutual reorientation. This flexibility has no physiological effect on acid viability or urease maturation in H. pylori. Atomistic molecular dynamics simulations suggest that Ni(II) binding is important for the conformational stability of the N-terminal helix. NMR chemical shift perturbation analysis indicates that no structural changes occur in the Zn-binding domain upon addition of Ni(II) in the pH 6.3-7.2 range. The structure of the Ni(II) binding site was probed using ¹H NMR spectroscopy experiments tailored to reveal hyperfine-shifted signals around the paramagnetic metal ion. On this basis, two possible models were derived using quantum-mechanical DFT calculations. The results provide a comprehensive picture of the Ni(II) mode to HpHypA, important to rationalize, at the molecular level, the functional interactions of this chaperone with its protein partners.

Fluorinated Paramagnetic Complexes: Sensitive and Responsive Probes for Magnetic Resonance Spectroscopy and Imaging

Fluorine magnetic resonance spectroscopy (MRS) and magnetic resonance imaging (MRI) of chemical and physiological processes is becoming more widespread. The strength of this technique comes from the negligible background signal in in vivo¹⁹F MRI and the large chemical shift window of ¹⁹F that enables it to image concomitantly more than one marker. These same advantages have also been successfully exploited in the design of responsive ¹⁹F probes. Part of the recent growth of this technique can be attributed to novel designs of ¹⁹F probes with improved imaging parameters due to the incorporation of paramagnetic metal ions. In this review, we provide a description of the theories and strategies that have been employed successfully to improve the sensitivity of ¹⁹F probes with paramagnetic metal ions. The Bloch-Wangsness-Redfield theory accurately predicts how molecular parameters such as internuclear distance, geometry, rotational correlation times, as well as the nature, oxidation state, and spin state of the metal ion affect the sensitivity of the fluorine-based probes. The principles governing the design of responsive ¹⁹F probes are subsequently described in a "how to" guide format. Examples of such probes and their advantages and disadvantages are highlighted through a synopsis of the literature.

Modulation of RNA primer formation by Mn(II)-substituted T7 DNA primase

Lagging strand DNA synthesis by DNA polymerase requires RNA primers produced by DNA primase. The N-terminal primase domain of the gene 4 protein of phage T7 comprises a zinc-binding domain that recognizes a specific DNA sequence and an RNA polymerase domain that catalyzes RNA polymerization. Based on its crystal structure, the RNA polymerase domain contains two Mg(II) ions. Mn(II) substitution leads to elevated RNA primer synthesis by T7 DNA primase. NMR analysis revealed that upon binding Mn(II), T7 DNA primase undergoes conformational changes near the metal cofactor binding site that are not observed when the enzyme binds Mg(II). A machine-learning algorithm called linear discriminant analysis (LDA) was trained by using the large collection of Mn(II) and Mg(II) binding sites available in the protein data bank (PDB). Application of the model to DNA primase revealed a preference in the enzyme's second metal binding site for Mn(II) over Mg(II), suggesting that T7 DNA primase activity modulation when bound to Mn(II) is based on structural changes in the enzyme.

NMR characterization of HtpG, the E. coli Hsp90, using sparse labeling with 13C-methyl alanine

A strategy for acquiring structural information from sparsely isotopically labeled large proteins is illustrated with an application to the E. coli heat-shock protein, HtpG (high temperature protein G), a 145 kDa dimer. It uses ¹³C-alanine methyl labeling in a perdeuterated background to take advantage of the sensitivity and resolution of Methyl-TROSY spectra, as well as the backbone-centered structural information from ¹H-¹³C residual dipolar couplings (RDCs) of alanine methyl groups. In all, 40 of the 47 expected crosspeaks were resolved and 36 gave RDC data. Assignments of crosspeaks were partially achieved by transferring assignments from those made on individual domains using triple resonance methods. However, these were incomplete and in many cases the transfer was ambiguous. A genetic algorithm search for consistency between predictions based on domain structures and measurements for chemical shifts and RDCs allowed 60% of the 40 resolved crosspeaks to be assigned with confidence. Chemical shift changes of these crosspeaks on adding an ATP analog to the apo-protein are shown to be consistent with structural changes expected on comparing previous crystal structures for apo- and complex- structures. RDCs collected on the assigned alanine methyl peaks are used to generate a new solution model for the apo-protein structure.

Insights into Carbohydrate Recognition by 3D Structure Determination of Protein–Carbohydrate Complexes Using NMR

This chapter provides an overview of protein–carbohydrate complex structures determined with NMR spectroscopy and deposited in the Protein Data Bank (PDB). These 14 structures include protein–carbohydrate interactions ranging from nanomolar to millimolar affinities. Two complexes are discussed in detail, one representing a tightly bound complex and one a weak but specific interaction. This review illustrates that NMR spectroscopy is a competitive method for three-dimensional structure determination of protein–carbohydrate complexes, especially in the case of weak interactions. The number of biological functions in which protein–carbohydrate interactions are involved is steadily growing. Essential functions of the immune system such as the distinction between self and non-self, or the resolution of inflammation, involve critical protein–carbohydrate recognition events. It is therefore expected that by providing atomic details, NMR spectroscopy can make a significant contribution in the near future to unexplored pathways of the immune system and of many other biological processes.

Recent excitements in protein NMR: Large proteins and biologically relevant dynamics

The advent of Transverse Relaxation Optimized SpectroscopY (TROSY) and perdeuteration allowed biomolecular NMR spectroscopists to overcome the size limitation barrier (approx. 20 kDa) in de novo structure determination of proteins. The utility of these techniques was immediately demonstrated on large proteins and protein complexes (e.g. GroELGroES, ClpP protease, Hsp90-p53, 20S proteasome, etc.). Further, recent methodological developments such as Residual Dipolar Couplings and Paramagnetic Relaxation Enhancement allowed accurate measurement of long-range structural restraints. Additionally, Carr-Purcell-Meiboom-Gill (CPMG), rotating frame relaxation experiments (R1(rho)) and saturation transfer experiments (CEST and DEST) created never-before accessibility to the (mu)s-ms timescale dynamic parameters that led to the deeper understanding of biological processes. Meanwhile, the excitement in the field continued with a series of developments in the fast data acquisition methods allowing rapid structural studies on less stable proteins. This review aims to discuss important developments in the field of biomolecular NMR spectroscopy in the recent past, i.e., in the post TROSY era. These developments not only gave access to the structural studies of large protein assemblies, but also revolutionized tools in the arsenal of today's biomolecular NMR and point to a bright future of biomolecular NMR spectroscopy.

Model-free extraction of spin label position distributions from pseudocontact shift data

Not a point, but a cloud: advanced PCS data analysis using 3D probability density reconstruction provides more information.

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³⁺ 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.

3D Computational Modeling of Proteins Using Sparse Paramagnetic NMR Data

Computational modeling of proteins using evolutionary or de novo approaches offers rapid structural characterization, but often suffers from low success rates in generating high quality models comparable to the accuracy of structures observed in X-ray crystallography or nuclear magnetic resonance (NMR) spectroscopy. A computational/experimental hybrid approach incorporating sparse experimental restraints in computational modeling algorithms drastically improves reliability and accuracy of 3D models. This chapter discusses the use of structural information obtained from various paramagnetic NMR measurements and demonstrates computational algorithms implementing pseudocontact shifts as restraints to determine the structure of proteins at atomic resolution.

Computing distance distributions from dipolar evolution data with overtones: RIDME spectroscopy with Gd(iii)-based spin labels

Extraction of distance distributions between high-spin paramagnetic centers from relaxation induced dipolar modulation enhancement (RIDME) data is affected by the presence of overtones of dipolar frequencies.

A protocol for the refinement of NMR structures using simultaneously pseudocontact shift restraints from multiple lanthanide ions

The binding of paramagnetic metal ions to proteins produces a number of different effects on the NMR spectra of the system. In particular, when the magnetic susceptibility of the metal ion is anisotropic, pseudocontact shifts (PCSs) arise and can be easily measured. They constitute very useful restraints for the solution structure determination of metal-binding proteins. In this context, there has been great interest in the use of lanthanide(III) ions to induce PCSs in diamagnetic proteins, e.g. through the replacement native calcium(II) ions. By preparing multiple samples in each of which a different ion of the lanthanide series is introduced, it is possible to obtain multiple independent PCS datasets that can be used synergistically to generate protein structure ensembles (typically called bundles). For typical NMR-based determination of protein structure, it is necessary to perform an energetic refinement of such initial bundles to obtain final structures whose geometric quality is suitable for deposition in the PDB. This can be conveniently done by using restrained molecular dynamics simulations (rMD) in explicit solvent. However, there are no available protocols for rMD using multiple PCS datasets as part of the restraints. In this work, we extended the PCS module of the AMBER MD package to handle multiple datasets and tuned a previously developed protocol for NMR structure refinement to achieve consistent convergence with PCS restraints. Test calculations with real experimental data show that this new implementation delivers the expected improvement of protein geometry, resulting in final structures that are of suitable quality for deposition. Furthermore, we observe that also initial structures generated only with traditional restraints can be successfully refined using traditional and PCS restraints simultaneously.

Pseudo-Contact NMR Shifts over the Paramagnetic Metalloprotein CoMMP-12 from First Principles

Long-range pseudo-contact NMR shifts (PCSs) provide important restraints for the structure refinement of proteins when a paramagnetic metal center is present, either naturally or introduced artificially. Here we show that ab initio quantum-chemical methods and a modern version of the Kurland-McGarvey approach for paramagnetic NMR (pNMR) shifts in the presence of zero-field splitting (ZFS) together provide accurate predictions of all PCSs in a metalloprotein (high-spin cobalt-substituted MMP-12 as a test case). Computations of 314 ¹³ C PCSs using g- and ZFS tensors based on multi-reference methods provide a reliable bridge between EPR-parameter- and susceptibility-based pNMR formalisms. Due to the high sensitivity of PCSs to even small structural differences, local structures based either on X-ray diffraction or on various DFT optimizations could be evaluated critically by comparing computed and experimental PCSs. Many DFT functionals provide insufficiently accurate structures. We also found the available 1RMZ PDB X-ray structure to exhibit deficiencies related to binding of a hydroxamate inhibitor. This has led to a newly refined PDB structure for MMP-12 (5LAB) that provides a more accurate coordination arrangement and PCSs.

Rapid identification of amino acid types in proteins using phase modulated 2D HN(CACB) and 2D HN(COCACB)

We present a simple approach to rapidly identify amino acid types in proteins from a 2D spectrum. The method is based on the fact that (13)C(β) chemical shifts of different amino acid types fall in distinct spectral regions. By evolving the (13)C chemical shifts in the conventional HNCACB or HN(CO)CACB type experiment for a single specified delay period, the phase of the cross peaks of different amino acid residues are modulated depending on their (13)C(β) shift values. Following this specified evolution period, the 2D HN projections of these experiments are acquired. The (13)C evolution period can be chosen such that all residues belonging to a given set of amino acid types have the same phase pattern (positive or negative) facilitating their identification. This approach does not require the preparation of any additional samples, involves the analysis of 2D [(15)N-(1)H] HSQC-type spectra obtained from the routinely used triple resonance experiments with minor modifications, and is applicable to deuterated proteins. The method will be useful for quick assignment of signals that shift during ligand binding or in combination with selective labeling/unlabeling approaches for identification of amino acid types to aid the sequential assignment process.

Paramagnetic effects on the NMR spectra of isotropic bicelles with headgroup modified chelator lipids and metal ions

We characterized the paramagnetic effects of nine metal ions on NMR signals of isotropic bicelles with headgroup-modified lipids. We found that Mn(2+), Gd(3+) and Dy(3+) show evidence for influencing NMR signals on the surface more than inside and on the disc edge, providing distance information in the bilayers.

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.

Nitrogen detected TROSY at high field yields high resolution and sensitivity for protein NMR

Detection of (15)N in multidimensional NMR experiments of proteins has sparsely been utilized because of the low gyromagnetic ratio (γ) of nitrogen and the presumed low sensitivity of such experiments. Here we show that selecting the TROSY components of proton-attached (15)N nuclei (TROSY (15)NH) yields high quality spectra in high field magnets (>600 MHz) by taking advantage of the slow (15)N transverse relaxation and compensating for the inherently low (15)N sensitivity. The (15)N TROSY transverse relaxation rates increase modestly with molecular weight but the TROSY gain in peak heights depends strongly on the magnetic field strength. Theoretical simulations predict that the narrowest line width for the TROSY (15)NH component can be obtained at 900 MHz, but sensitivity reaches its maximum around 1.2 GHz. Based on these considerations, a (15)N-detected 2D (1)H-(15)N TROSY-HSQC ((15)N-detected TROSY-HSQC) experiment was developed and high-quality 2D spectra were recorded at 800 MHz in 2 h for 1 mM maltose-binding protein at 278 K (τc ~ 40 ns). Unlike for (1)H detected TROSY, deuteration is not mandatory to benefit (15)N detected TROSY due to reduced dipolar broadening, which facilitates studies of proteins that cannot be deuterated, especially in cases where production requires eukaryotic expression systems. The option of recording (15)N TROSY of proteins expressed in H2O media also alleviates the problem of incomplete amide proton back exchange, which often hampers the detection of amide groups in the core of large molecular weight proteins that are expressed in D2O culture media and cannot be refolded for amide back exchange. These results illustrate the potential of (15)NH-detected TROSY experiments as a means to exploit the high resolution offered by high field magnets near and above 1 GHz.

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.

Paramagnetic NMR probes for characterization of the dynamic conformations and interactions of oligosaccharides

Paramagnetism-assisted nuclear magnetic resonance (NMR) techniques have recently been applied to a wide variety of biomolecular systems, using sophisticated immobilization methods to attach paramagnetic probes, such as spin labels and lanthanide-chelating groups, at specific sites of the target biomolecules. This is also true in the field of carbohydrate NMR spectroscopy. NMR analysis of oligosaccharides is often precluded by peak overlap resulting from the lack of variability of local chemical structures, by the insufficiency of conformational restraints from nuclear Overhauser effect (NOE) data due to low proton density, and moreover, by the inherently flexible nature of carbohydrate chains. Paramagnetic probes attached to the reducing ends of oligosaccharides cause paramagnetic relaxation enhancements (PREs) and/or pseudocontact shifts (PCSs) resolve the peak overlap problem. These spectral perturbations can be sources of long-range atomic distance information, which complements the local conformational information derived from J couplings and NOEs. Furthermore, paramagnetic NMR approaches, in conjunction with computational methods, have opened up possibilities for the description of dynamic conformational ensembles of oligosaccharides in solution. Several applications of paramagnetic NMR techniques are presented to demonstrate their utility for characterizing the conformational dynamics of oligosaccharides and for probing the carbohydrate-recognition modes of proteins. These techniques can be applied to the characterization of transient, non-stoichiometric interactions and will contribute to the visualization of dynamic biomolecular processes involving sugar chains.

Capturing conformational States in proteins using sparse paramagnetic NMR data

Capturing conformational changes in proteins or protein-protein complexes is a challenge for both experimentalists and computational biologists. Solution nuclear magnetic resonance (NMR) is unique in that it permits structural studies of proteins under greatly varying conditions, and thus allows us to monitor induced structural changes. Paramagnetic effects are increasingly used to study protein structures as they give ready access to rich structural information of orientation and long-range distance restraints from the NMR signals of backbone amides, and reliable methods have become available to tag proteins with paramagnetic metal ions site-specifically and at multiple sites. In this study, we show how sparse pseudocontact shift (PCS) data can be used to computationally model conformational states in a protein system, by first identifying core structural elements that are not affected by the environmental change, and then computationally completing the remaining structure based on experimental restraints from PCS. The approach is demonstrated on a 27 kDa two-domain NS2B-NS3 protease system of the dengue virus serotype 2, for which distinct closed and open conformational states have been observed in crystal structures. By changing the input PCS data, the observed conformational states in the dengue virus protease are reproduced without modifying the computational procedure. This data driven Rosetta protocol enables identification of conformational states of a protein system, which are otherwise difficult to obtain either experimentally or computationally.

Integral membrane protein structure determination using pseudocontact shifts

Obtaining enough experimental restraints can be a limiting factor in the NMR structure determination of larger proteins. This is particularly the case for large assemblies such as membrane proteins that have been solubilized in a membrane-mimicking environment. Whilst in such cases extensive deuteration strategies are regularly utilised with the aim to improve the spectral quality, these schemes often limit the number of NOEs obtainable, making complementary strategies highly beneficial for successful structure elucidation. Recently, lanthanide-induced pseudocontact shifts (PCSs) have been established as a structural tool for globular proteins. Here, we demonstrate that a PCS-based approach can be successfully applied for the structure determination of integral membrane proteins. Using the 7TM α-helical microbial receptor pSRII, we show that PCS-derived restraints from lanthanide binding tags attached to four different positions of the protein facilitate the backbone structure determination when combined with a limited set of NOEs. In contrast, the same set of NOEs fails to determine the correct 3D fold. The latter situation is frequently encountered in polytopical α-helical membrane proteins and a PCS approach is thus suitable even for this particularly challenging class of membrane proteins. The ease of measuring PCSs makes this an attractive route for structure determination of large membrane proteins in general.