Numbat: an interactive software tool for fitting Deltachi-tensors to molecular coordinates using pseudocontact shifts

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.

The evolution of paramagnetic NMR as a tool in structural biology

Paramagnetic NMR data contain extremely accurate long-range information on metalloprotein structures and, when used in the frame of integrative structural biology approaches, they allow for the retrieval of structural details to a resolution that is not achievable using other techniques. Paramagnetic data thus represent an extremely powerful tool to refine protein models in solution, especially when coupled to X-ray or cryoelectron microscopy data, to monitor the formation of complexes and determine the relative arrangements of their components, and to highlight the presence of conformational heterogeneity. More recently, theoretical and computational advancements in quantum chemical calculations of paramagnetic NMR observables are progressively opening new routes in structural biology, because they allow for the determination of the structure within the coordination sphere of the metal center, thus acting as a loupe on sites that are difficult to observe but very important for protein function.

Incorporation of residual chemical shift anisotropy into the treatment of 15N pseudocontact shifts for structural refinement

Paramagnetic NMR experiments, including the pseudocontact shift experiment, have seen increasing use due to recently developed probes and labeling strategies. The pseudocontact shift experiment can provide valuable intra- or inter-molecular distance and orientation information. However, the use of ¹H/¹³C or ¹H/¹⁵N PCS data in structure calculations is currently complicated by the contribution of residual chemical shift anisotropy to the ¹³C or ¹⁵N datasets. Here, we present a corrected PCS energy term for the software package Xplor-NIH with the appropriate residual chemical shift anisotropy correction and show its suitability for model refinements of ubiquitin labeled at residue 57 with a Tm-M8-SPy tag. For data taken at 800 MHz, the improvement with the corrected energy term is sufficient to make the quality of the fit for the ¹⁵N dataset comparable to that of the ¹H dataset, for which no correction is needed. The corrected energy term is expected to become more relevant with increased use of higher field instruments and as new paramagnetic probes with larger magnetic susceptibility tensors continue to be developed.

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.

DEER experiments reveal fundamental differences between calmodulin complexes with IQ and MARCKS peptides in solution

Calmodulin (CaM) is a calcium-binding protein that regulates the function of many proteins by indirectly conferring Ca²⁺ sensitivity, and it undergoes a large conformational change on partners' binding. We compared the solution binding mode of the target peptides MARCKS and IQ by double electron-electron resonance (DEER) distance measurements and paramagnetic NMR. We combined nitroxide and Gd(III) spin labels, including specific substitution of one of the Ca²⁺ ions in the CaM mutant N60D by a Gd(III) ion. The binding of MARCKS to holo-CaM resulted neither in a closed conformation nor in a unique relative orientation between the two CaM domains, in contrast with the crystal structure. Binding of IQ to holo-CaM did generate a closed conformation. Using elastic network modeling and 12 distance restraints obtained from multiple holo-CaM/IQ DEER data, we derived a model of the solution structure, which is in reasonable agreement with the crystal structure.

Structure determination of high-energy states in a dynamic protein ensemble

Macromolecular function frequently requires that proteins change conformation into high-energy states^1-4. However, methods for solving the structures of these functionally essential, lowly populated states are lacking. Here we develop a method for high-resolution structure determination of minorly populated states by coupling NMR spectroscopy-derived pseudocontact shifts⁵ (PCSs) with Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion⁶ (PCS-CPMG). Our approach additionally defines the corresponding kinetics and thermodynamics of high-energy excursions, thereby characterizing the entire free-energy landscape. Using a large set of simulated data for adenylate kinase (Adk), calmodulin and Src kinase, we find that high-energy PCSs accurately determine high-energy structures (with a root mean squared deviation of less than 3.5 angström). Applying our methodology to Adk during catalysis, we find that the high-energy excursion involves surprisingly small openings of the AMP and ATP lids. This previously unresolved high-energy structure solves a longstanding controversy about conformational interconversions that are rate-limiting for catalysis. Primed for either substrate binding or product release, the high-energy structure of Adk suggests a two-step mechanism combining conformational selection to this state, followed by an induced-fit step into a fully closed state for catalysis of the phosphoryl-transfer reaction. Unlike other methods for resolving high-energy states, such as cryo-electron microscopy and X-ray crystallography, our solution PCS-CPMG approach excels in cases involving domain rearrangements of smaller systems (less than 60 kDa) and populations as low as 0.5%, and enables the simultaneous determination of protein structure, kinetics and thermodynamics while proteins perform their function.

Organoarsenic probes to study proteins by NMR spectroscopy

Arsenical probes enable structural studies of proteins. We report the first organoarsenic probes for nuclear magnetic resonance (NMR) and electron paramagnetic resonance (EPR) spectroscopy to study proteins in solutions. These probes can be attached to irregular loop regions. A lanthanide-binding tag induces sizable pseudocontact shifts in protein NMR spectra of a magnitude never observed for small paramagnetic probes before.

Characterizing conformational ensembles of multi-domain proteins using anisotropic paramagnetic NMR restraints

It has been over two decades since paramagnetic NMR started to form part of the essential techniques for structural analysis of proteins under physiological conditions. Paramagnetic NMR has significantly expanded our understanding of the inherent flexibility of proteins, in particular, those that are formed by combinations of two or more domains. Here, we present a brief overview of techniques to characterize conformational ensembles of such multi-domain proteins using paramagnetic NMR restraints produced through anisotropic metals, with a focus on the basics of anisotropic paramagnetic effects, the general procedures of conformational ensemble reconstruction, and some representative reweighting approaches.

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.

Rigid, Highly Reactive and Stable DOTA-like Tags Containing a Thiol-Specific Phenylsulfonyl Pyridine Moiety for Protein Modification and NMR Analysis*

Site specific installation of a paramagnetic ion with magnetic anisotropy in a biomolecule generates valuable structural restraints, such as pseudocontact shifts (PCSs) and residual dipolar couplings (RDCs). These paramagnetic effects can be used to characterize the structures, interactions and dynamics of biological macromolecules and their complexes. Two single-armed DOTA-like tags, BrPSPy-DO3M(S)A-Ln and BrPSPy-6M-DO3M(S)A-Ln, each containing a thiol-specific reacting group, that is, a phenylsulfonyl pyridine moiety, are demonstrated as rigid, reactive and stable paramagnetic tags for protein modification by formation of a reducing resistant thioether bond between the protein and the tag. The two tags present high reactivity with the solvent exposed thiol group in aqueous solution at room temperature. The introduction of Br at the meta-position in pyridine enhances the reactivity of 4-phenylsulfonyl pyridine towards the solvent exposed thiol group in a protein, whereas the ortho-methyl group in pyridine increases the rigidity of the tag in the protein conjugates. The high performance of these two tags has been demonstrated in different cysteine mutants of ubiquitin and GB1. The high reactivity and rigidity of these two tags can be added in the toolbox of paramagnetic tags suitable for the high-resolution NMR measurements of biological macromolecules and their complexes.

Advances in NMR Spectroscopy of Weakly Aligned Biomolecular Systems

The measurement and application of residual dipolar couplings (RDCs) in solution NMR studies of biological macromolecules has become well established over the past quarter of a century. Numerous methods for generating the requisite anisotropic orientational molecular distribution have been demonstrated, each with its specific strengths and weaknesses. In parallel, an enormous number of pulse schemes have been introduced to measure the many different types of RDCs, ranging from the most widely measured backbone amide ¹⁵N-¹H RDCs, to ¹H-¹H RDCs and couplings between low-γ nuclei. Applications of RDCs range from structure validation and refinement to the determination of relative domain orientations, the measurement of backbone and domain motions, and de novo structure determination. Nevertheless, it appears that the power of the RDC methodology remains underutilized. This review aims to highlight the practical aspects of sample preparation and RDC measurement while describing some of the most straightforward applications that take advantage of the exceptionally precise information contained in such data. Some emphasis will be placed on more recent developments that enable the accurate measurement of RDCs in larger systems, which is key to the ongoing shift in focus of biological NMR spectroscopy from structure determination toward gaining improved understanding of how molecular flexibility drives protein function.

An automated iterative approach for protein structure refinement using pseudocontact shifts

NMR structure calculation using NOE-derived distance restraints requires a considerable number of assignments of both backbone and sidechains resonances, often difficult or impossible to get for large or complex proteins. Pseudocontact shifts (PCSs) also play a well-established role in NMR protein structure calculation, usually to augment existing structural, mostly NOE-derived, information. Existing refinement protocols using PCSs usually either require a sizeable number of sidechain assignments or are complemented by other experimental restraints. Here, we present an automated iterative procedure to perform backbone protein structure refinements requiring only a limited amount of backbone amide PCSs. Already known structural features from a starting homology model, in this case modules of repeat proteins, are framed into a scaffold that is subsequently refined by experimental PCSs. The method produces reliable indicators that can be monitored to judge about the performance. We applied it to a system in which sidechain assignments are hardly possible, designed Armadillo repeat proteins (dArmRPs), and we calculated the solution NMR structure of YM4A, a dArmRP containing four sequence-identical internal modules, obtaining high convergence to a single structure. We suggest that this approach is particularly useful when approximate folds are known from other techniques, such as X-ray crystallography, while avoiding inherent artefacts due to, for instance, crystal packing.

Conformational ensemble of a multidomain protein explored by Gd3+ electron paramagnetic resonance

Despite their importance in function, the conformational state of proteins and its changes are often poorly understood, mainly because of the lack of an efficient tool. MurD, a 47-kDa protein enzyme responsible for peptidoglycan biosynthesis, is one of those proteins whose conformational states and changes during their catalytic cycle are not well understood. Although it has been considered that MurD takes a single conformational state in solution as shown by a crystal structure, the solution nuclear magnetic resonance (NMR) study suggested the existence of multiple conformational state of apo MurD in solution. However, the conformational distribution has not been evaluated. In this work, we investigate the conformational states of MurD by the use of electron paramagnetic resonance (EPR), especially intergadolinium distance measurement using double electron-electron resonance (DEER) measurement. The gadolinium ions are fixed on specific positions on MurD via a rigid double-arm paramagnetic lanthanide tag that has been originally developed for paramagnetic NMR. The combined use of NMR and EPR enables accurate interpretation of the DEER distance information to the structural information of MurD. The DEER distance measurement for apo MurD shows a broad distance distribution, whereas the presence of the inhibitor narrows the distance distribution. The results suggest that MurD exists in a wide variety of conformational states in the absence of ligands, whereas binding of the inhibitor eliminates variation in conformational states. The multiple conformational states of MurD were previously implied by NMR experiments, but our DEER data provided structural characterization of the conformational variety of MurD.

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.

Dynamic Exchange of the Metal Chelating Moiety: A Key Factor in Determining the Rigidity of Protein-Tag Conjugates in Paramagnetic NMR

Site-specific labeling of proteins with a paramagnetic tag is an efficient way to provide atomic-resolution information about the dynamics, interactions, and structures of the proteins and protein-ligand complexes. The paramagnetic effects manifested in NMR spectroscopy generally contain paramagnetic relaxation enhancement, pseudocontact shifts (PCSs), and residual dipolar coupling (RDC), and these effects correlate closely with the flexibility of protein-tag conjugates. The rigidity of the paramagnetic tag is greatly important in decoding the structural details of macromolecular complexes, because paramagnetic averaging reduces the PCSs and RDCs. Here we show that the dynamic exchange of the metal chelating moiety is a key factor in determining the rigidity of the paramagnetic tag in the protein conjugates. Decreasing the conformational exchange rates in the metal chelating moiety greatly minimizes the paramagnetic averaging and thus increases PCSs and RDCs. This effect has been demonstrated in an open-chain tag, Py-l-Cys-DTPA, which generates large PCSs and RDCs that are comparable to those of the reported cyclic DOTA-like tags. The proposed route offers a unique way to design suitable paramagnetic tags for applications in biological systems.

NMR pseudocontact shifts in a symmetric protein homotrimer

NMR pseudocontact shifts are a valuable tool for structural and functional studies of proteins. Protein multimers mediate key functional roles in biology, but methods for their study by pseudocontact shifts are so far not available. Paramagnetic tags attached to identical subunits in multimeric proteins cause a combined pseudocontact shift that cannot be described by the standard single-point model. Here, we report pseudocontact shifts generated simultaneously by three paramagnetic Tm-M7PyThiazole-DOTA tags to the trimeric molecular chaperone Skp and provide an approach for the analysis of this and related symmetric systems. The pseudocontact shifts were described by a "three-point" model, in which positions and parameters of the three paramagnetic tags were fitted. A good correlation between experimental data and predicted values was found, validating the approach. The study establishes that pseudocontact shifts can readily be applied to multimeric proteins, offering new perspectives for studies of large protein complexes by paramagnetic NMR spectroscopy.

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.

Paramagpy: software for fitting magnetic susceptibility tensors using paramagnetic effects measured in NMR spectra

Paramagnetic metal ions with fast-relaxing electrons generate pseudocontact shifts (PCSs), residual dipolar couplings (RDCs), paramagnetic relaxation enhancements (PREs) and cross-correlated relaxation (CCR) in the nuclear magnetic resonance (NMR) spectra of the molecules they bind to. These effects offer long-range structural information in molecules equipped with binding sites for such metal ions. Here we present the new open-source software Paramagpy, which has been written in Python 3 with a graphic user interface. Paramagpy combines the functionalities of different currently available programs to support the fitting of magnetic susceptibility tensors using PCS, RDC, PRE and CCR data and molecular coordinates in Protein Data Bank (PDB) format, including a convenient graphical user interface. Paramagpy uses efficient fitting algorithms to avoid local minima and supports corrections to back-calculated PCS and PRE data arising from cross-correlation effects with chemical shift tensors. The source code is available from 10.5281/zenodo.3594568 .

A novel, rationally designed lanthanoid chelating tag delivers large paramagnetic structural restraints for biomolecular NMR

A novel, rationally designed lanthanoid chelating tag enables fast ligation to biomacromolecules and delivers long-range structural restraints by NMR.

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.

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.

Capturing a dynamically interacting inhibitor by paramagnetic NMR spectroscopy

Transient and fuzzy intermolecular interactions are fundamental to many biological processes. Despite their importance, they are notoriously challenging to characterize. Effects induced by paramagnetic ligands in the NMR spectra of interacting biomolecules provide an opportunity to amplify subtle manifestations of weak intermolecular interactions observed for diamagnetic ligands. Here, we present an approach to characterizing dynamic interactions between a partially flexible dimeric protein, HIV-1 protease, and a metallacarborane-based ligand, a system for which data obtained by standard NMR approaches do not enable detailed structural interpretation. We show that for the case where the experimental data are significantly averaged to values close to zero the standard fitting of pseudocontact shifts cannot provide reliable structural information. We based our approach on generating a large ensemble of full atomic models, for which the experimental data can be predicted, ensemble averaged and finally compared to the experiment. We demonstrate that a combination of paramagnetic NMR experiments, quantum chemical calculations, and molecular dynamics simulations offers a route towards structural characterization of dynamic protein-ligand complexes.

A Double-Armed, Hydrophilic Transition Metal Complex as a Paramagnetic NMR Probe

Synthetic metal complexes can be used as paramagnetic probes for the study of proteins and protein complexes. Herein, two transition metal NMR probes (TraNPs) are reported. TraNPs are attached through two arms to a protein to generate a pseudocontact shift (PCS) using cobalt(II), or paramagnetic relaxation enhancement (PRE) with manganese(II). The PCS analysis of TraNPs attached to three different proteins shows that the size of the anisotropic component of the magnetic susceptibility depends on the probe surroundings at the surface of the protein, contrary to what is observed for lanthanoid‐based probes. The observed PCS are relatively small, making cobalt‐based probes suitable for localized studies, such as of an active site. The obtained PREs are stronger than those obtained with nitroxide spin labels and the possibility to generate both PCS and PRE offers advantages. The properties of TraNPs in comparison with other cobalt‐based probes are discussed.

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.

How Do Nuclei Couple to the Magnetic Moment of a Paramagnetic Center? A New Theory at the Gauntlet of the Experiments

The recent derivation, based on pure quantum chemistry (QC) first-principles, of the pseudocontact shifts (PCSs) caused by a paramagnetic metal center on far away nuclei has cast doubts on the validity of the semiempirical (SE) theory, predicting PCSs to arise from the metal magnetic susceptibility anisotropy. The SE theory has been used and applied countless times, especially in the last 2 decades, to obtain structural information on proteins containing paramagnetic metal ions. We show here that the QC and SE predictions can be directly tested against experiments, provided a suitable macromolecular system is used. The SE approach yields a good prediction of the experimental PCSs while the QC one does not. It appears that the classic theory is able to grasp satisfactorily the underlying physics.

Protein Structure Determination in Living Cells

To date, in-cell NMR has elucidated various aspects of protein behaviour by associating structures in physiological conditions. Meanwhile, current studies of this method mostly have deduced protein states in cells exclusively based on 'indirect' structural information from peak patterns and chemical shift changes but not 'direct' data explicitly including interatomic distances and angles. To fully understand the functions and physical properties of proteins inside cells, it is indispensable to obtain explicit structural data or determine three-dimensional (3D) structures of proteins in cells. Whilst the short lifetime of cells in a sample tube, low sample concentrations, and massive background signals make it difficult to observe NMR signals from proteins inside cells, several methodological advances help to overcome the problems. Paramagnetic effects have an outstanding potential for in-cell structural analysis. The combination of a limited amount of experimental in-cell data with software for ab initio protein structure prediction opens an avenue to visualise 3D protein structures inside cells. Conventional nuclear Overhauser effect spectroscopy (NOESY)-based structure determination is advantageous to elucidate the conformations of side-chain atoms of proteins as well as global structures. In this article, we review current progress for the structure analysis of proteins in living systems and discuss the feasibility of its future works.

Localization of ligands within human carbonic anhydrase II using 19F pseudocontact shift analysis

Unraveling the native structure of protein-ligand complexes in solution enables rational drug design. We report here the use of ¹⁹F pseudocontact shift (PCS) NMR as a method to determine fluorine positions of high affinity ligands bound within the drug target human carbonic anhydrase II with high accuracy. Three different ligands were localized within the protein by analysis of the obtained PCS from simple one-dimensional ¹⁹F spectra with an accuracy of up to 0.8 Å. In order to validate the PCS, four to five independent magnetic susceptibility tensors induced by lanthanide chelating tags bound site-specifically to single cysteine mutants were refined. Least-squares minimization and a Monte-Carlo approach allowed the assessment of experimental errors on the intersection of the corresponding four to five PCS isosurfaces. By defining an angle score that reflects the relative isosurface orientation for different tensor combinations, it was established that the ligand can be localized accurately using only three tensors, if the isosurfaces are close to orthogonal. For two out of three ligands, the determined position closely matched the X-ray coordinates. Our results for the third ligand suggest, in accordance with previously reported ab initio calculations, a rotated position for the difluorophenyl substituent, enabling a favorable interaction with Phe-131. The lanthanide-fluorine distance varied between 22 and 38 Å and induced ¹⁹F PCS ranged from 0.078 to 0.409 ppm, averaging to 0.213 ppm. Accordingly, even longer metal-fluorine distances will lead to meaningful PCS, rendering the investigation of protein-ligand complexes significantly larger than 30 kDa feasible.

Structural Basis of Protein Kinase Cα Regulation by the C-Terminal Tail

Protein kinase C (PKC) isoenzymes are multi-modular proteins activated at the membrane surface to regulate signal transduction processes. When activated by second messengers, PKC undergoes a drastic conformational and spatial transition from the inactive cytosolic state to the activated membrane-bound state. The complete structure of either state of PKC remains elusive. We demonstrate, using NMR spectroscopy, that the isolated Ca²⁺-sensing membrane-binding C2 domain of the conventional PKCα interacts with a conserved hydrophobic motif of the kinase C-terminal region, and we report a structural model of the complex. Our data suggest that the C-terminal region plays a dual role in regulating the PKC activity: activating, through sensitization of PKC to intracellular Ca²⁺ oscillations; and auto-inhibitory, through its interaction with a conserved positively charged region of the C2 domain.

Conformationally locked lanthanide chelating tags for convenient pseudocontact shift protein nuclear magnetic resonance spectroscopy

Pseudocontact shifts (PCS) generated by lanthanide chelating tags yield valuable restraints for investigating protein structures, dynamics and interactions in solution. In this work, dysprosium-, thulium- and terbium-complexes of eight-fold methylated 1,4,7,10-tetraazacyclododecane-1,4,7,10-tetraacetic acid tags [DOTA-M8-(4R4S)-SSPy] are presented that induce large pseudocontact shifts up to 5.5 ppm and adopt exclusively the square antiprismatic conformation. This is in contrast to our earlier findings on complexes of the stereoisomeric DOTA-M8-(8S)-SSPy, where significant amounts of the twisted square antiprismatic conformer for the Dy tag were observed. The Dy-, Tm-, Tb- and Lu-complexes of DOTA-M8-(4R4S)-SSPy were conjugated to ubiquitin S57C and selectively ¹⁵N leucine labeled human carbonic anhydrase II S50C, resulting in only one set of signals. Furthermore, we investigated the conformation of the thulium- and dysprosium-complexes in vacuo and with implicit water solvent using density functional theory calculations. The calculated energy differences between the two different conformations (7.0-50.5 kJ/mol) and experimental evidence from the corresponding ytterbium- and yttrium-complexes clearly suggest a SAP [Λ(δδδδ)] geometry for the complexes presented in this study. The lanthanide chelating tag studied in this work offer insights into the solution structure of proteins by inducing strong pseudocontact shifts, show different tensor properties compared to its predecessor, enables a convenient assignment procedure, is accessed by a more economic synthesis than its predecessor and constitutes a highly promising starting point for further developments of lanthanide chelating tags.

Ligand Proton Pseudocontact Shifts Determined from Paramagnetic Relaxation Dispersion in the Limit of NMR Intermediate Exchange

Delineation of protein-ligand interaction modes is key for rational drug discovery. The availability of complex crystal structures is often limited by the aqueous solubility of the compounds, while lead-like compounds with micromolar affinities normally fall into the NMR intermediate exchange regime, in which severe line broadening to beyond the detection of interfacial resonances limits NMR applications. Here, we developed a new method to retrieve low-populated bound-state ¹H pseudocontact shifts (PCSs) using paramagnetic relaxation dispersion (RD). We evaluated using a ¹H PCS-RD approach in a BRM bromodomain lead-like inhibitor to filter molecular docking poses using multiple intermolecular structural restraints. Considering the universal presence of proton atoms in druglike compounds, our work will have wide application in structure-guided drug discovery even under an extreme condition of NMR intermediate exchange and low aqueous solubility of ligands.

Farseer-NMR: automatic treatment, analysis and plotting of large, multi-variable NMR data

We present Farseer-NMR ( https://git.io/vAueU ), a software package to treat, evaluate and combine NMR spectroscopic data from sets of protein-derived peaklists covering a range of experimental conditions. The combined advances in NMR and molecular biology enable the study of complex biomolecular systems such as flexible proteins or large multibody complexes, which display a strong and functionally relevant response to their environmental conditions, e.g. the presence of ligands, site-directed mutations, post translational modifications, molecular crowders or the chemical composition of the solution. These advances have created a growing need to analyse those systems' responses to multiple variables. The combined analysis of NMR peaklists from large and multivariable datasets has become a new bottleneck in the NMR analysis pipeline, whereby information-rich NMR-derived parameters have to be manually generated, which can be tedious, repetitive and prone to human error, or even unfeasible for very large datasets. There is a persistent gap in the development and distribution of software focused on peaklist treatment, analysis and representation, and specifically able to handle large multivariable datasets, which are becoming more commonplace. In this regard, Farseer-NMR aims to close this longstanding gap in the automated NMR user pipeline and, altogether, reduce the time burden of analysis of large sets of peaklists from days/weeks to seconds/minutes. We have implemented some of the most common, as well as new, routines for calculation of NMR parameters and several publication-quality plotting templates to improve NMR data representation. Farseer-NMR has been written entirely in Python and its modular code base enables facile extension.

Complete assignment of Ala, Ile, Leu, Met and Val methyl groups of human blood group A and B glycosyltransferases using lanthanide-induced pseudocontact shifts and methyl-methyl NOESY

Human blood group A and B glycosyltransferases (GTA, GTB) are highly homologous glycosyltransferases. A number of high-resolution crystal structures is available showing that these enzymes convert from an open conformation into a catalytically active closed conformation upon substrate binding. However, the mechanism of glycosyltransfer is still under debate, and the precise nature as well as the time scales of conformational transitions are unknown. NMR offers a variety of experiments to shine more light on these unresolved questions. Therefore, in a first step we have assigned all methyl resonance signals in MILVA labeled samples of GTA and GTB, still a challenging task for 70 kDa homodimeric proteins. Assignments were obtained from methyl-methyl NOESY experiments, and from measurements of lanthanide-induced pseudocontact shifts (PCS) using high resolution crystal structures as templates. PCSs and chemical shift perturbations, induced by substrate analogue binding, suggest that the fully closed state is not adopted in the presence of lanthanide ions.

Short two-armed lanthanide-binding tags for paramagnetic NMR spectroscopy based on chiral 1,4,7,10-tetrakis(2-hydroxypropyl)-1,4,7,10-tetraazacyclododecane scaffolds

A new pair of enantiomeric two-armed lanthanide-binding tags have been developed for paramagnetic NMR studies of proteins. The tags produce large and significantly different paramagnetic effects to one another when bound to the same tagging site. Additionally, they are less sensitive to sample pH than our previous two-armed tag designs.

Paramagnetic NMR as a new tool in structural biology

NMR (nuclear magnetic resonance) investigation through the exploitation of paramagnetic effects is passing from an approach limited to few specialists in the field to a generally applicable method that must be considered, especially for the characterization of systems hardly affordable with other techniques. This is mostly due to the fact that paramagnetic data are long range in nature, thus providing information for the structural and dynamic characterization of complex biomolecular architectures in their native environment. On the other hand, this information usually needs to be complemented by data from other sources. Integration of paramagnetic NMR with other techniques, and the development of protocols for a joint analysis of all available data, is fundamental for achieving a comprehensive characterization of complex biological systems. We describe here a few examples of the new possibilities offered by paramagnetic data used in integrated structural approaches.

Stable and rigid DTPA-like paramagnetic tags suitable for in vitro and in situ protein NMR analysis

Organic synthesis of a ligand with high binding affinities for paramagnetic lanthanide ions is an effective way of generating paramagnetic effects on proteins. These paramagnetic effects manifested in high-resolution NMR spectroscopy are valuable dynamic and structural restraints of proteins and protein-ligand complexes. A paramagnetic tag generally contains a metal chelating moiety and a reactive group for protein modification. Herein we report two new DTPA-like tags, 4PS-PyDTTA and 4PS-6M-PyDTTA that can be site-specifically attached to a protein with a stable thioether bond. Both protein-tag adducts form stable lanthanide complexes, of which the binding affinities and paramagnetic tensors are tunable with respect to the 6-methyl group in pyridine. Paramagnetic relaxation enhancement (PRE) effects of Gd(III) complex on protein-tag adducts were evaluated in comparison with pseudocontact shift (PCS), and the results indicated that both 4PS-PyDTTA and 4PS-6M-PyDTTA tags are rigid and present high-quality PREs that are crucially important in elucidation of the dynamics and interactions of proteins and protein-ligand complexes. We also show that these two tags are suitable for in-situ protein NMR analysis.

3D structure determination of a protein in living cells using paramagnetic NMR spectroscopy

Determining the three-dimensional structure of a protein in living cells remains particularly challenging. We demonstrated that the integration of site-specific tagging proteins and GPS-Rosetta calculations provides a fast and effective way of determining the structures of proteins in living cells, and in principle the interactions and dynamics of protein-ligand complexes.

Long-range paramagnetic NMR data can provide a closer look on metal coordination in metalloproteins

Paramagnetic NMR data can be profitably incorporated in structural refinement protocols of metalloproteins or metal-substituted proteins, mostly as distance or angle restraints. However, they could in principle provide much more information, because the magnetic susceptibility of a paramagnetic metal ion is largely determined by its coordination sphere. This information can in turn be used to evaluate changes occurring in the coordination sphere of the metal when ligands (e.g.: inhibitors) are bound to the protein. This gives an experimental handle on the molecular structure in the vicinity of the metal which falls in the so-called blind sphere. The magnetic susceptibility anisotropy tensors of cobalt(II) and nickel(II) ions bound to human carbonic anhydrase II in free and inhibited forms have been determined. The change of the magnetic susceptibility anisotropy is directly linked to the binding mode of different ligands in the active site of the enzyme. Indication about the metal coordination sphere in the presence of an inhibitor in pharmaceutically relevant proteins could be important in the design of selective drugs with a structure-based approach.

Methyl group assignment using pseudocontact shifts with PARAssign

A new version of the program PARAssign has been evaluated for assignment of NMR resonances of the 76 methyl groups in leucines, isoleucines and valines in a 25 kDa protein, using only the structure of the protein and pseudocontact shifts (PCS) generated with a lanthanoid tag at up to three attachment sites. The number of reliable assignments depends strongly on two factors. The principle axes of the magnetic susceptibility tensors of the paramagnetic centers should not be parallel so as to avoid correlated PCS. Second, the fraction of resonances in the spectrum of a paramagnetic sample that can be paired with the diamagnetic counterparts is critical for the assignment. With the data from two tag positions a reliable assignment could be obtained for 60% of the methyl groups and for many of the remaining resonances the number of possible assignments is limited to two or three. With a single tag, reliable assignments can be obtained for methyl groups with large PCS near the tag. It is concluded that assignment of methyl group resonances by paramagnetic tagging can be particularly useful in combination with some additional data, such as from mutagenesis or NOE-based experiments. Approaches to yield the best assignment results with PCS generating tags are discussed.

New Lanthanide Tag for the Generation of Pseudocontact Shifts in DNA by Site-Specific Ligation to a Phosphorothioate Group

Pseudocontact shifts (PCS) generated by paramagnetic lanthanides provide a rich source of long-range structural restraints that can readily be measured by nuclear magnetic resonance (NMR) spectroscopy. Many different lanthanide-binding tags have been designed for site-specific tagging of proteins, but established routes for tagging DNA with a single metal ion rely on difficult chemical synthesis. Here we present a simple and practical strategy for site-specific tagging of inexpensive phosphorothioate (PT) oligonucleotides. Commercially available PT oligonucleotides are diastereomers with S and R stereoconfiguration at the backbone PT site. The respective S_P and R_P diastereomers can readily be separated by HPLC. A new alkylating lanthanide-binding tag, C10, was synthesized that delivered quantitative tagging yields with both diastereomers. PCSs were observed following ligation with the complementary DNA strand to form double-stranded DNA duplexes. The PCSs were larger for the S_P than the R_P oligonucleotide and good correlation between back-calculated and experimental PCSs was observed. The C10 tag can also be attached to cysteine residues in proteins, where it generates a stable thioether bond. Ligated to the A28C mutant of ubiquitin, the tag produced excellent fits of magnetic susceptibility anisotropy (Δχ) tensors, with larger tensors than for the tagged PT oligonucleotides, indicating that the tag is not completely immobilized after ligation with a PT group.