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

Probing Biomolecular Interactions with Paramagnetic Nuclear Magnetic Resonance Spectroscopy

Recent advances in computational methods like AlphaFold have transformed structural biology, enabling accurate modeling of protein complexes and driving applications in drug discovery and protein engineering. However, predicting the structure of systems involving weak, transient, or dynamic interactions, or of complexes with disordered regions, remains challenging. Nuclear Magnetic Resonance (NMR) spectroscopy offers atomic-level insights into biomolecular complexes, even in weakly interacting and dynamic systems. Paramagnetic NMR, in particular, provides long-range structural restraints, easily exceeding distances over 25 Å, making it ideal for studying large protein complexes. Advances in chemical tools for introducing paramagnetic tags into proteins, combined with progress in electron paramagnetic resonance (EPR) spectroscopy, have enhanced the method's utility. This perspective article discusses paramagnetic NMR approaches for analyzing biomolecular complexes in solution and in the solid state, emphasizing quantities like pseudocontact shifts, residual dipolar couplings, and paramagnetic relaxation enhancements. Additionally, dynamic nuclear polarization offers a promising method to amplify NMR signals of large complexes, even in complex environments. The integration of AlphaFold protein structure prediction with paramagnetic NMR holds great potential for advancing our understanding of biomolecular interactions.

Solution Structure of a Lanthanide-binding DNA Aptamer Determined Using High Quality pseudocontact shift restraints

In this contribution we report the high-resolution NMR structure of a recently identified lanthanide-binding aptamer (LnA). We demonstrate that the rigid lanthanide binding by LnA allows for the measurement of anisotropic paramagnetic NMR restraints which to date remain largely inaccessible for nucleic acids. One type of such restraints - pseudocontact shifts (PCS) induced by four different paramagnetic lanthanides - was extensively used throughout the current structure determination study and the measured PCS turned out to be exceptionally well reproduced by the final aptamer structure. This finding opens the perspective for a broader application of paramagnetic effects in NMR studies of nucleic acids through the transplantation of the binding site found in LnA into other DNA/RNA systems.

Methyl probes in proteins for determining ligand binding mode in weak protein-ligand complexes

Structures of protein–ligand complexes provide critical information for drug design. Most protein–ligand complex structures are determined using X-ray crystallography, but where crystallography is not able to generate a structure for a complex, NMR is often the best alternative. However, the available tools to enable rapid and robust structure determination of protein–ligand complexes by NMR are currently limited. This leads to situations where projects are either discontinued or pursued without structural data, rendering the task more difficult. We previously reported the NMR Molecular Replacement (NMR²) approach that allows the structure of a protein–ligand complex to be determined without requiring the cumbersome task of protein resonance assignment. Herein, we describe the NMR² approach to determine the binding pose of a small molecule in a weak protein–ligand complex by collecting sparse protein methyl-to-ligand NOEs from a selectively labeled protein sample and an unlabeled ligand. In the selective labeling scheme all methyl containing residues of the protein are protonated in an otherwise deuterated background. This allows measurement of intermolecular NOEs with greater sensitivity using standard NOESY pulse sequences instead of isotope-filtered NMR experiments. This labelling approach is well suited to the NMR² approach and extends its utility to include larger protein–ligand complexes.

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.

Localising nuclear spins by pseudocontact shifts from a single tagging site

Ligating a protein at a specific site with a tag molecule containing a paramagnetic metal ion provides a versatile way of generating pseudocontact shifts (PCSs) in nuclear magnetic resonance (NMR) spectra. PCSs can be observed for nuclear spins far from the tagging site, and PCSs generated from multiple tagging sites have been shown to enable highly accurate structure determinations at specific sites of interest, even when using flexible tags, provided the fitted effective magnetic susceptibility anisotropy (Δ χ ) tensors accurately back-calculate the experimental PCSs measured in the immediate vicinity of the site of interest. The present work investigates the situation where only the local structure of a protein region or bound ligand is to be determined rather than the structure of the entire molecular system. In this case, the need for gathering structural information from tags deployed at multiple sites may be queried. Our study presents a computational simulation of the structural information available from samples produced with single tags attached at up to six different sites, up to six different tags attached to a single site, and in-between scenarios. The results indicate that the number of tags is more important than the number of tagging sites. This has important practical implications, as it is much easier to identify a single site that is suitable for tagging than multiple ones. In an initial experimental demonstration with the ubiquitin mutant S57C, PCSs generated with four different tags at a single site are shown to accurately pinpoint the location of amide protons in different segments of the protein.

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.

Localising individual atoms of tryptophan side chains in the metallo-β-lactamase IMP-1 by pseudocontact shifts from paramagnetic lanthanoid tags at multiple sites

The metallo-β -lactamase IMP-1 features a flexible loop near the active site that assumes different conformations in single crystal structures, which may assist in substrate binding and enzymatic activity. To probe the position of this loop, we labelled the tryptophan residues of IMP-1 with 7-13 C-indole and the protein with lanthanoid tags at three different sites. The magnetic susceptibility anisotropy (Δ χ ) tensors were determined by measuring pseudocontact shifts (PCSs) of backbone amide protons. The Δ χ tensors were subsequently used to identify the atomic coordinates of the tryptophan side chains in the protein. The PCSs were sufficient to determine the location of Trp28, which is in the active site loop targeted by our experiments, with high accuracy. Its average atomic coordinates showed barely significant changes in response to the inhibitor captopril. It was found that localisation spaces could be defined with better accuracy by including only the PCSs of a single paramagnetic lanthanoid ion for each tag and tagging site. The effect was attributed to the shallow angle with which PCS isosurfaces tend to intersect if generated by tags and tagging sites that are identical except for the paramagnetic lanthanoid ion.

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.

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.

The molecular basis of OH-PCB estrogen receptor activation

Polychlorinated bisphenols (PCBs) continue to contaminate food chains globally where they concentrate in tissues and disrupt the endocrine systems of species throughout the ecosphere. Hydroxylated PCBs (OH-PCBs) are major PCB metabolites and high-affinity inhibitors of human estrogen sulfotransferase (SULT1E1), which sulfonates estrogens and thus prevents them from binding to and activating their receptors. OH-PCB inhibition of SULT1E1 is believed to contribute significantly to PCB-based endocrine disruption. Here, for the first time, the molecular basis of OH-PCB inhibition of SULT1E1 is revealed in a structure of SULT1E1 in complex with OH-PCB1 (4ʹ-OH-2,6-dichlorobiphenol) and its substrates, estradiol (E2), and PAP (3’-phosphoadenosine-5-phosphosulfate). OH-PCB1 prevents catalysis by intercalating between E2 and catalytic residues and establishes a new E2-binding site whose E2 affinity and positioning are greater than and competitive with those of the reactive-binding pocket. Such complexes have not been observed previously and offer a novel template for the design of high-affinity inhibitors. Mutating residues in direct contact with OH-PCB weaken its affinity without compromising the enzyme’s catalytic parameters. These OH-PCB resistant mutants were used in stable transfectant studies to demonstrate that OH-PCBs regulate estrogen receptors in cultured human cell lines by binding the OH-PCB binding pocket of SULT1E1.

Using NMR Spectroscopy in the Fragment-Based Drug Discovery of Small-Molecule Anticancer Targeted Therapies

Against the challenge of providing personalized cancer care, the development of targeted therapies stands as a promising approach. The discovery of these agents can benefit from fragment-based drug discovery (FBDD) methods that help guide ligand design and provide key structural information on the targets of interest. In particular, nuclear magnetic resonance spectroscopy is a promising biophysical tool in fragment discovery due to its detection capabilities and versatility. This review provides an overview of FBDD, describes the basis of NMR-based fragment screening, summarizes some exciting technical advances reported over the past decades, and closes with a discussion of selected case studies where this technique has been used as part of drug discovery campaigns to produce lead compounds towards the design of anti-cancer targeted therapies.

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.

CcpNmr AnalysisScreen, a new software programme with dedicated automated analysis tools for fragment-based drug discovery by NMR

Fragment-based drug discovery or FBDD is one of the main methods used by industry and academia for identifying drug-like candidates in early stages of drug discovery. NMR has a significant impact at any stage of the drug discovery process, from primary identification of small molecules to the elucidation of binding modes for guiding optimisations. The essence of NMR as an analytical tool, however, requires the processing and analysis of relatively large amounts of single data items, e.g. spectra, which can be daunting when managed manually. One bottleneck in FBDD by NMR is a lack of adequate and well-integrated resources for NMR data analysis that are freely available to the community. Thus, scientists typically resort to manually inspecting large datasets and relying predominantly on subjective interpretations. In this manuscript, we present CcpNmr AnalysisScreen, a software package that provides computational tools for automated analysis of FBDD data by NMR. We outline how the quality of collected spectra can be evaluated quickly, and how robust workflows can be optimised for reliable and rapid hit identification. With an intuitive graphical user interface and powerful algorithms, AnalysisScreen enables easy analysis of the large datasets needed in the early process of drug discovery by NMR.

Protein-ligand structure determination with the NMR molecular replacement tool, NMR2

We recently reported on a new method called NMR Molecular Replacement that efficiently derives the structure of a protein-ligand complex at the interaction site. The method was successfully applied to high and low affinity complexes covering ligands from peptides to small molecules. The algorithm used in the NMR Molecular Replacement program has until now not been described in detail. Here, we present a complete description of the NMR Molecular Replacement implementation as well as several new features that further reduce the time required for structure elucidation.

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.

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.

Gadolinium Complexes as Contrast Agent for Cellular NMR Spectroscopy

Aqua Gd³⁺ and Gd-DOTA (gadolinium-1,4,7,10-tetraazacyclododecane-1,4,7,10-tetraacete) complexes were studied as a contrast agent in cellular NMR (nuclear magnetic resonance) spectroscopy for distinguishing between intracellular and extracellular spaces. The contrast agents for this purpose should provide strong paramagnetic relaxation enhancement and localize in the extracellular space without disturbing biological functions. Cell membrane permeability to Gd complexes was evaluated from the concentrations of gadolinium complexes in the inside and outside of E. coli cells measured by the ¹H-NMR relaxation. The site-specific binding of the complexes to E. coli cells was also analyzed by high-resolution solid-state ¹³C-NMR. The aqua Gd³⁺ complex did not enhance T₁ relaxation in proportion to the amount of added Gd³⁺. This Gd³⁺ concentration dependence and the ¹³C-NMR indicated that its strong cytotoxicity should be due to the binding of the paramagnetic ions to cellular components especially at the lipid membranes. In contrast, Gd-DOTA stayed in the solution states and enhanced relaxation in proportion to the added amount. This agent exhibited strong T₁ contrast between the intra- and extracellular spaces by a factor of ten at high concentrations under which the cells were viable over a long experimental time of days. These properties make Gd-DOTA suitable for selectively contrasting the living cellular space in NMR spectroscopy primarily owing to its weak interaction with cellular components.

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.

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.

Integrative Protein Modeling in RosettaNMR from Sparse Paramagnetic Restraints

Computational methods to predict protein structure from nuclear magnetic resonance (NMR) restraints that only require assignment of backbone signals, hold great potential to study larger proteins. Ideally, computational methods designed to work with sparse data need to add atomic detail that is missing in the experimental restraints. We introduce a comprehensive framework into the Rosetta suite that uses NMR restraints derived from paramagnetic labeling. Specifically, RosettaNMR incorporates pseudocontact shifts, residual dipolar couplings, and paramagnetic relaxation enhancements. It continues to use backbone chemical shifts and nuclear Overhauser effect distance restraints. We assess RosettaNMR for protein structure prediction by folding 28 monomeric proteins and 8 homo-oligomeric proteins. Furthermore, the general applicability of RosettaNMR is demonstrated on two protein-protein and three protein-ligand docking examples. Paramagnetic restraints generated more accurate models for 85% of the benchmark proteins and, when combined with chemical shifts, sampled high-accuracy models (≤2Å) in 50% of the cases.

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.

Paramagnetic NMR in solution and the solid state

The field of paramagnetic NMR has expanded considerably in recent years. This review addresses both the theoretical description of paramagnetic NMR, and the way in which it is currently practised. We provide a review of the theory of the NMR parameters of systems in both solution and the solid state. Here we unify the different languages used by the NMR, EPR, quantum chemistry/DFT, and magnetism communities to provide a comprehensive and coherent theoretical description. We cover the theory of the paramagnetic shift and shift anisotropy in solution both in the traditional formalism in terms of the magnetic susceptibility tensor, and using a more modern formalism employing the relevant EPR parameters, such as are used in first-principles calculations. In addition we examine the theory first in the simple non-relativistic picture, and then in the presence of spin-orbit coupling. These ideas are then extended to a description of the paramagnetic shift in periodic solids, where it is necessary to include the bulk magnetic properties, such as magnetic ordering at low temperatures. The description of the paramagnetic shift is completed by describing the current understanding of such shifts due to lanthanide and actinide ions. We then examine the paramagnetic relaxation enhancement, using a simple model employing a phenomenological picture of the electronic relaxation, and again using a more complex state-of-the-art theory which incorporates electronic relaxation explicitly. An additional important consideration in the solid state is the impact of bulk magnetic susceptibility effects on the form of the spectrum, where we include some ideas from the field of classical electrodynamics. We then continue by describing in detail the solution and solid-state NMR methods that have been deployed in the study of paramagnetic systems in chemistry, biology, and the materials sciences. Finally we describe a number of case studies in paramagnetic NMR that have been specifically chosen to highlight how the theory in part one, and the methods in part two, can be used in practice. The systems chosen include small organometallic complexes in solution, solid battery electrode materials, metalloproteins in both solution and the solid state, systems containing lanthanide ions, and multi-component materials used in pharmaceutical controlled-release formulations that have been doped with paramagnetic species to measure the component domain sizes.

Real-Time Insights into Biological Events: In-Cell Processes and Protein-Ligand Interactions

FlowNMR has the aim of continuously monitoring processes that occur in conditions that are not compatible with being carried out within a closed tube. However, it is sample intensive and not suitable for samples, such as proteins or living cells, that are often available in limited volumes and possibly low concentrations. We here propose a dialysis-based modification of a commercial flowNMR setup that allows for recycling the medium while confining the sample (proteins and cells) within the active volume of the tube. This approach is demonstrated in the specific cases of in-cell NMR and protein-based ligand studies.

Selective Methyl Labeling of Proteins: Enabling Structural and Mechanistic Studies As Well As Drug Discovery Applications by Solution-State NMR

Escherichia coli expression protocols for selective labeling of methyl groups in proteins have been essential in expanding the size range of targets that can be studied by biomolecular NMR. Based on the initial work achieving selective labeling of isoleucine, leucine, and valine residues, additional methods were developed over the past years which enabled the individual and/or simultaneous combinatorial labeling of all methyl containing amino acids. Together with the introduction of new methyl-optimized NMR experiments, this now allows the detailed characterization of protein-ligand interactions as well as mechanistic and dynamic processes of protein-protein complexes up to 1MDa in size. In this chapter, we provide a general introduction to selective labeling of proteins using E. coli-based expression systems, describe the considerations taken into account prior to the selective labeling of a protein, and include the protocols used to produce such proteins. An overview of applications using selectively labeled proteins with an emphasis on examples relevant to the drug discovery process is then presented.

Methyl group reorientation under ligand binding probed by pseudocontact shifts

Liquid-state NMR spectroscopy is a powerful technique to elucidate binding properties of ligands on proteins. Ligands binding in hydrophobic pockets are often in close proximity to methyl groups and binding can lead to subtle displacements of methyl containing side chains to accommodate the ligand. To establish whether pseudocontact shifts can be used to characterize ligand binding and the effects on methyl groups, the N-terminal domain of HSP90 was tagged with caged lanthanoid NMR probe 5 at three positions and titrated with a ligand. Binding was monitored using the resonances of leucine and valine methyl groups. The pseudocontact shifts (PCS) caused by ytterbium result in enhanced dispersion of the methyl spectrum, allowing more resonances to be observed. The effects of tag attachment on the spectrum and ligand binding are small. Significant changes in PCS were observed upon ligand binding, indicating displacements of several methyl groups. By determining the cross-section of PCS iso-surfaces generated by two or three paramagnetic centers, the new position of a methyl group can be estimated, showing displacements in the range of 1-3 Å for methyl groups in the binding site. The information about such subtle but significant changes may be used to improve docking studies and can find application in fragment-based drug discovery.

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.

Investigating Protein-Ligand Interactions by Solution Nuclear Magnetic Resonance Spectroscopy

Protein-ligand interactions are of fundamental importance in almost all processes in living organisms. The ligands comprise small molecules, drugs or biological macromolecules and their interaction strength varies over several orders of magnitude. Solution NMR spectroscopy offers a large repertoire of techniques to study such complexes. Here, we give an overview of the different NMR approaches available. The information they provide ranges from the simple information about the presence of binding or epitope mapping to the complete 3 D structure of the complex. NMR spectroscopy is particularly useful for the study of weak interactions and for the screening of binding ligands with atomic resolution.

Fast NMR Screening of Macromolecular Complexes by a Paramagnetic Spin Relaxation Filter

The paramagnetic spin relaxation filter is described for the rapid NMR screening of intermolecular interactions between ligands and macromolecular anionic receptors with large transverse relaxation enhancements (R 2p). The addition of micromolar concentrations of Gd3+ to the mixture produces the immediate broadening/suppression of the NMR signals of interacting species while leaving unaffected those of noncompetitive binders (one-dimensional and two-dimensional experiments). The method is highly sensitive, unveiling interactions that are too weak to generate changes in chemical shifts or relaxation times. It is operationally very simple and hence, it is amenable to ready implementation by nonspecialists. Examples of application such as detecting the formation of interpolymer complexes, cyclodextrin host-guest interactions, and the screening of DNA ligands are included that demonstrate the reliability and broad applicability of the method.

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.

The NSAID allosteric site of human cytosolic sulfotransferases

Non-steroidal anti-inflammatory drugs (NSAIDs) are among the most commonly prescribed drugs worldwide-more than 111 million prescriptions were written in the United States in 2014. NSAIDs allosterically inhibit cytosolic sulfotransferases (SULTs) with high specificity and therapeutically relevant affinities. This study focuses on the interactions of SULT1A1 and mefenamic acid (MEF)-a potent, highly specific NSAID inhibitor of 1A1. Here, the first structure of an NSAID allosteric site-the MEF-binding site of SULT1A1-is determined using spin-label triangulation NMR. The structure is confirmed by site-directed mutagenesis and provides a molecular framework for understanding NSAID binding and isoform specificity. The mechanism of NSAID inhibition is explored using molecular dynamics and equilibrium and pre-steady-state ligand-binding studies. MEF inhibits SULT1A1 turnover through an indirect (helix-mediated) stabilization of the closed form of the active-site cap of the enzyme, which traps the nucleotide and slows its release. Using the NSAID-binding site structure of SULT1A1 as a comparative model, it appears that 11 of the 13 human SULT isoforms harbor an NSAID-binding site. We hypothesize that these sites evolved to enable SULT isoforms to respond to metabolites that lie within their metabolic domains. Finally, the NSAID-binding site structure offers a template for developing isozyme-specific allosteric inhibitors that can be used to regulate specific areas of sulfuryl-transfer metabolism.

NMR in structure-based drug design

NMR spectroscopy is a powerful technique that can provide valuable structural information for drug discovery endeavors. Here, we discuss the strengths (and limitations) of NMR applications to structure-based drug discovery, highlighting the different levels of resolution and throughput obtainable. Additionally, the emerging field of paramagnetic NMR in drug discovery and recent developments in approaches to speed up and automate protein-observed NMR data collection and analysis are discussed.

NMR studies of ligand binding

NMR spectroscopy is an established tool in drug discovery, but its strength is commonly regarded to be largely confined to the early stages of hit discovery and fragment based drug design, where NMR offers unique capabilities of characterizing the binding modes of ligand molecules that bind sufficiently weakly to be in rapid exchange between bound and free state. Here we, first, provide a meta-review of recent reviews on NMR studies of ligand binding and, second, review recent progress towards NMR characterization of the ligand binding mode in stable protein-ligand complexes, with particular emphasis on the global positioning system (GPS) approach enabled by paramagnetic lanthanide tags.

Perspectives on paramagnetic NMR from a life sciences infrastructure

The effects arising in NMR spectroscopy because of the presence of unpaired electrons, collectively referred to as "paramagnetic NMR" have attracted increasing attention over the last decades. From the standpoint of the structural and mechanistic biology, paramagnetic NMR provides long range restraints that can be used to assess the accuracy of crystal structures in solution and to improve them by simultaneous refinements through NMR and X-ray data. These restraints also provide information on structure rearrangements and conformational variability in biomolecular systems. Theoretical improvements in quantum chemistry calculations can nowadays allow for accurate calculations of the paramagnetic data from a molecular structural model, thus providing a tool to refine the metal coordination environment by matching the paramagnetic effects observed far away from the metal. Furthermore, the availability of an improved technology (higher fields and faster magic angle spinning) has promoted paramagnetic NMR applications in the fast-growing area of biomolecular solid-state NMR. Major improvements in dynamic nuclear polarization have been recently achieved, especially through the exploitation of the Overhauser effect occurring through the contact-driven relaxation mechanism: the very large enhancement of the ¹³C signal observed in a variety of liquid organic compounds at high fields is expected to open up new perspectives for applications of solution NMR.

Hot-spot identification on a broad class of proteins and RNA suggest unifying principles of molecular recognition

Chemically diverse fragments tend to collectively bind at localized sites on proteins, which is a cornerstone of fragment-based techniques. A central question is how general are these strategies for predicting a wide variety of molecular interactions such as small molecule-protein, protein-protein and protein-nucleic acid for both experimental and computational methods. To address this issue, we recently proposed three governing principles, (1) accurate prediction of fragment-macromolecule binding free energy, (2) accurate prediction of water-macromolecule binding free energy, and (3) locating sites on a macromolecule that have high affinity for a diversity of fragments and low affinity for water. To test the generality of these concepts we used the computational technique of Simulated Annealing of Chemical Potential to design one small fragment to break the RecA-RecA protein-protein interaction and three fragments that inhibit peptide-deformylase via water-mediated multi-body interactions. Experiments confirm the predictions that 6-hydroxydopamine potently inhibits RecA and that PDF inhibition quantitatively tracks the water-mediated binding predictions. Additionally, the principles correctly predict the essential bound waters in HIV Protease, the surprisingly extensive binding site of elastase, the pinpoint location of electron transfer in dihydrofolate reductase, the HIV TAT-TAR protein-RNA interactions, and the MDM2-MDM4 differential binding to p53. The experimental confirmations of highly non-obvious predictions combined with the precise characterization of a broad range of known phenomena lend strong support to the generality of fragment-based methods for characterizing molecular recognition.

Tetrahydrobiopterin regulates monoamine neurotransmitter sulfonation

Monoamine neurotransmitters are among the hundreds of signaling small molecules whose target interactions are switched "on" and "off" via transfer of the sulfuryl-moiety (-SO3) from PAPS (3'-phosphoadenosine 5'-phosphosulfate) to the hydroxyls and amines of their scaffolds. These transfer reactions are catalyzed by a small family of broad-specificity enzymes-the human cytosolic sulfotransferases (SULTs). The first structure of a SULT allosteric-binding site (that of SULT1A1) has recently come to light. The site is conserved among SULT1 family members and is promiscuous-it binds catechins, a naturally occurring family of flavanols. Here, the catechin-binding site of SULT1A3, which sulfonates monoamine neurotransmitters, is modeled on that of 1A1 and used to screen in silico for endogenous metabolite 1A3 allosteres. Screening predicted a single high-affinity allostere, tetrahydrobiopterin (THB), an essential cofactor in monoamine neurotransmitter biosynthesis. THB is shown to bind and inhibit SULT1A3 with high affinity, 23 (±2) nM, and to bind weakly, if at all, to the four other major SULTs found in brain and liver. The structure of the THB-bound binding site is determined and confirms that THB binds the catechin site. A structural comparison of SULT1A3 with SULT1A1 (its immediate evolutionary progenitor) reveals how SULT1A3 acquired high affinity for THB and that the majority of residue changes needed to transform 1A1 into 1A3 are clustered at the allosteric and active sites. Finally, sequence records reveal that the coevolution of these sites played an essential role in the evolution of simian neurotransmitter metabolism.

Impact of spin label rigidity on extent and accuracy of distance information from PRE data

Paramagnetic relaxation enhancement (PRE) is a versatile tool for NMR spectroscopic structural and kinetic studies in biological macromolecules. Here, we compare the quality of PRE data derived from two spin labels with markedly different dynamic properties for large RNAs using the I-A riboswitch aptamer domain (78 nt) from Mesoplamsa florum as model system. We designed two I-A aptamer constructs that were spin-labeled by noncovalent hybridization of short spin-labeled oligomer fragments. As an example of a flexible spin label, ^UreidoU-TEMPO was incorporated into the 3' terminal end of helix P1 while, the recently developed rigid spin-label Çm was incorporated in the 5' terminal end of helix P1. We determined PRE rates obtained from aromatic ¹³C bound proton intensities and compared these rates to PREs derived from imino proton intensities in this sizeable RNA (~78 nt). PRE restraints derived from both imino and aromatic protons yielded similar data quality, and hence can both be reliably used for PRE determination. For NMR, the data quality derived from the rigid spin label Çm is slightly better than the data quality for the flexible ^UreidoTEMPO as judged by comparison of the structural agreement with the I-A aptamer crystal structure (3SKI).

Structure restraints from heteronuclear pseudocontact shifts generated by lanthanide tags at two different sites

Pseudocontact shifts (PCS) encode long-range information on 3D structures of protein backbones and side-chains. The level of structural detail that can be obtained increases with the number of different sites tagged with a paramagnetic metal ion to generate PCSs. Here we show that PCSs from two different sites can suffice to determine the structure of polypeptide chains and their location and orientation relative to the magnetic susceptibility tensor χ, provided that PCSs are available for ¹H as well as heteronuclear spins. In addition, PCSs from two different sites are shown to provide detailed structural information on the conformation of methyl group-bearing amino-acid side-chains. A previously published ensemble structure of ubiquitin is shown to explain the magnetic susceptibility and alignment tensors slightly better than structures that try to explain the experimental data by a single conformation, illustrating the potential of PCSs as a tool to investigate small conformational changes.

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.

The structure of the catechin-binding site of human sulfotransferase 1A1

We are just beginning to understand the allosteric regulation of the human cytosolic sulfotransferase (SULTs) family-13 disease-relevant enzymes that regulate the activities of hundreds, if not thousands, of signaling small molecules. SULT1A1, the predominant isoform in adult liver, harbors two noninteracting allosteric sites, each of which binds a different molecular family: the catechins (naturally occurring flavonols) and nonsteroidal antiinflammatory drugs (NSAIDs). Here, we present the structure of an SULT allosteric binding site-the catechin-binding site of SULT1A1 bound to epigallocatechin gallate (EGCG). The allosteric pocket resides in a dynamic region of the protein that enables EGCG to control opening and closure of the enzyme's active-site cap. Furthermore, the structure offers a molecular explanation for the isozyme specificity of EGCG, which is corroborated experimentally. The binding-site structure was obtained without X-ray crystallography or multidimensional NMR. Instead, a SULT1A1 apoprotein structure was used to guide positioning of a small number of spin-labeled single-Cys mutants that coat the entire enzyme surface with a paramagnetic field of sufficient strength to determine its contribution to the bound ligand's transverse (T₂) relaxation from its 1D solution spectrum. EGCG protons were mapped to the protein surface by triangulation using the T₂ values to calculate their distances to a trio of spin-labeled Cys mutants. The final structure was obtained using distance-constrained molecular dynamics docking. This approach, which is readily extensible to other systems, is applicable over a wide range of ligand affinities, requires little protein, avoids the need for isotopically labeled protein, and has no protein molecular weight limitations.