Mesodynamics in the SARS nucleocapsid measured by NMR field cycling

Comprehensive analysis of relaxation decays from high-resolution relaxometry

Relaxometry consists in measuring relaxation rates over orders of magnitude of magnetic fields to probe motions of complex systems. High-resolution relaxometry (HRR) experiments can be performed on conventional high-field NMR magnets equipped with a sample shuttle. During the experiment, the sample shuttle transfers the sample between the high-field magnetic center and a chosen position in the stray field for relaxation during a variable delay, thus using the stray field as a variable field. As the relaxation delay occurs outside of the probe, HRR experiments cannot rely on the control of cross-relaxation pathways, which is standard in high-field relaxation pulse sequences. Thus, decay rates are not pure relaxation rates, which may impair a reliable description of the dynamics. Previously, we took into account cross-relaxation effects in the analysis of high-resolution relaxometry data by applying a correction factor to relaxometry decay rates in order to estimate relaxation rates. These correction factors were obtained from the iterative simulation of the relaxation decay while the sample lies outside of the probe and a preceding analysis of relaxation rates which relies on the approximation of a priori multi-exponential decays by mono-exponential functions. However, an analysis protocol matching directly experimental and simulated relaxometry decays should be more self consistent and more generally applicable as it can accommodate deviations from mono-exponential decays. Here, we introduce Matching INtensities for the Optimization of Timescales and Amplitudes of motions Under Relaxometry (MINOTAUR), a framework for the analysis of high-resolution relaxometry that takes as input the intensity decays at all fields. This approach uses the full relaxation matrix to calculate intensity decays, allowing complex relaxation pathways to be taken into account. Therefore, it eliminates the need for a correction of decay rates and for fitting multi-exponential decays with mono-exponential functions. The MINOTAUR software is designed as a flexible framework where relaxation matrices and spectral density functions corresponding to various models of motions can be defined on a case-by-case basis. The agreement with our previous analyses of protein side-chain dynamics from carbon-13 relaxation is excellent, while providing a more robust analysis tool. We expect MINOTAUR to become the tool of choice for the analysis of high-resolution relaxometry.

The preference signature of the SARS-CoV-2 Nucleocapsid NTD for its 5'-genomic RNA elements

The nucleocapsid protein (N) of SARS-CoV-2 plays a pivotal role during the viral life cycle. It is involved in RNA transcription and accounts for packaging of the large genome into virus particles. N manages the enigmatic balance of bulk RNA-coating versus precise RNA-binding to designated cis-regulatory elements. Numerous studies report the involvement of its disordered segments in non-selective RNA-recognition, but how N organizes the inevitable recognition of specific motifs remains unanswered. We here use NMR spectroscopy to systematically analyze the interactions of N's N-terminal RNA-binding domain (NTD) with individual cis RNA elements clustering in the SARS-CoV-2 regulatory 5'-genomic end. Supported by broad solution-based biophysical data, we unravel the NTD RNA-binding preferences in the natural genome context. We show that the domain's flexible regions read the intrinsic signature of preferred RNA elements for selective and stable complex formation within the large pool of available motifs.

How does it really move? Recent progress in the investigation of protein nanosecond dynamics by NMR and simulation

Nuclear magnetic resonance (NMR) spin relaxation experiments currently probe molecular motions on timescales from picoseconds to nanoseconds. The detailed interpretation of these motions in atomic detail benefits from complementarity with the results from molecular dynamics (MD) simulations. In this mini-review, we describe the recent developments in experimental techniques to study the backbone dynamics from ¹⁵N relaxation and side-chain dynamics from ¹³C relaxation, discuss the different analysis approaches from model-free to dynamics detectors, and highlight the many ways that NMR relaxation experiments and MD simulations can be used together to improve the interpretation and gain insights into protein dynamics.

Convergent views on disordered protein dynamics from NMR and computational approaches

Intrinsically disordered proteins (IDPs) or intrinsically disordered regions (IDRs) is a class of biologically important proteins exhibiting specific biophysical characteristics. They lack a hydrophobic core, and their conformational behavior is strongly influenced by electrostatic interactions. IDPs and IDRs are highly dynamic, and a characterization of the motions of IDPs and IDRs is essential for their physically correct description. NMR together with molecular dynamics simulations are the methods best suited to such a task because they provide information about dynamics of proteins with atomistic resolution. Here, we present a study of motions of a disordered C-terminal domain of the delta subunit of RNA polymerase from Bacillus subtilis. Positively and negatively charged residues in the studied domain form transient electrostatic contacts critical for the biological function. Our study is focused on investigation of ps-ns dynamics of backbone of the delta subunit based on analysis of amide 15N NMR relaxation data and molecular dynamics simulations. In order to extend an informational content of NMR data to lower frequencies, which are more sensitive to slower motions, we combined standard (high-field) NMR relaxation experiments with high-resolution relaxometry. Altogether, we collected data reporting the relaxation at 12 different magnetic fields, resulting in an unprecedented data set. Our results document that the analysis of such data provides a consistent description of dynamics and confirms the validity of so far used protocols of the analysis of dynamics of IDPs also for a partially folded protein. In addition, the potential to access detailed description of motions at the timescale of tens of ns with the help of relaxometry data is discussed. Interestingly, in our case, it appears to be mostly relevant for a region involved in the formation of temporary contacts within the disordered region, which was previously proven to be biologically important.

The Role of Disordered Regions in Orchestrating the Properties of Multidomain Proteins: The SARS-CoV-2 Nucleocapsid Protein and Its Interaction with Enoxaparin

Novel and efficient strategies need to be developed to interfere with the SARS-CoV-2 virus. One of the most promising pharmaceutical targets is the nucleocapsid protein (N), responsible for genomic RNA packaging. N is composed of two folded domains and three intrinsically disordered regions (IDRs). The globular RNA binding domain (NTD) and the tethered IDRs are rich in positively charged residues. The study of the interaction of N with polyanions can thus help to elucidate one of the key driving forces responsible for its function, i.e., electrostatics. Heparin, one of the most negatively charged natural polyanions, has been used to contrast serious cases of COVID-19 infection, and we decided to study its interaction with N at the molecular level. We focused on the NTR construct, which comprises the NTD and two flanking IDRs, and on the NTD construct in isolation. We characterized this interaction using different nuclear magnetic resonance approaches and isothermal titration calorimetry. With these tools, we were able to identify an extended surface of NTD involved in the interaction. Moreover, we assessed the importance of the IDRs in increasing the affinity for heparin, highlighting how different tracts of these flexible regions modulate the interaction.

NMR Reveals Specific Tracts within the Intrinsically Disordered Regions of the SARS-CoV-2 Nucleocapsid Protein Involved in RNA Encountering

The SARS-CoV-2 nucleocapsid (N) protein is crucial for the highly organized packaging and transcription of the genomic RNA. Studying atomic details of the role of its intrinsically disordered regions (IDRs) in RNA recognition is challenging due to the absence of structure and to the repetitive nature of their primary sequence. IDRs are known to act in concert with the folded domains of N and here we use NMR spectroscopy to identify the priming events of N interacting with a regulatory SARS-CoV-2 RNA element. ¹³C-detected NMR experiments, acquired simultaneously to ¹H detected ones, provide information on the two IDRs flanking the N-terminal RNA binding domain (NTD) within the N-terminal region of the protein (NTR, 1–248). We identify specific tracts of the IDRs that most rapidly sense and engage with RNA, and thus provide an atom-resolved picture of the interplay between the folded and disordered regions of N during RNA interaction.

Phospholipids in Motion: High-Resolution 31P NMR Field Cycling Studies

Diverse phospholipid motions are key to membrane function but can be quite difficult to untangle and quantify. High-resolution field cycling ³¹P NMR spin-lattice relaxometry, where the sample is excited at high field, shuttled in the magnet bore for low-field relaxation, then shuttled back to high field for readout of the residual magnetization, provides data on phospholipid dynamics and structure. This information is encoded in the field dependence of the ³¹P spin-lattice relaxation rate (R₁). In the field range from 11.74 down to 0.003 T, three dipolar nuclear magnetic relaxation dispersions (NMRDs) and one due to ³¹P chemical shift anisotropy contribute to R₁ of phospholipids. Extraction of correlation times and maximum relaxation amplitudes for these NMRDs provides (1) lateral diffusion constants for different phospholipids in the same bilayer, (2) estimates of how additives alter the motion of the phospholipid about its long axis, and (3) an average ³¹P-¹H angle with respect to the bilayer normal, which reveals that polar headgroup motion is not restricted on a microsecond timescale. Relative motions within a phospholipid are also provided by comparing ³¹P NMRD profiles for specifically deuterated molecules as well as ¹³C and ¹H field dependence profiles to that of ³¹P. Although this work has dealt exclusively with phospholipids in small unilamellar vesicles, these same NMRDs can be measured for phospholipids in micelles and nanodisks, making this technique useful for monitoring lipid behavior in a variety of structures and assessing how additives alter specific lipid motions.

The Inherent Dynamics and Interaction Sites of the SARS-CoV-2 Nucleocapsid N-Terminal Region

The nucleocapsid protein is one of four structural proteins encoded by SARS-CoV-2 and plays a central role in packaging viral RNA and manipulating the host cell machinery, yet its dynamic behavior and promiscuity in nucleotide binding has made standard structural methods to address its atomic-resolution details difficult. To begin addressing the SARS-CoV-2 nucleocapsid protein interactions with both RNA and the host cell along with its dynamic behavior, we have specifically focused on the folded N-terminal domain (NTD) and its flanking regions using nuclear magnetic resonance solution studies. Studies performed here reveal a large repertoire of interactions, which includes a temperature-dependent self-association mediated by the disordered flanking regions that also serve as binding sites for host cell cyclophilin-A while nucleotide binding is largely mediated by the central NTD core. NMR studies that include relaxation experiments have revealed the complicated dynamic nature of this viral protein. Specifically, while much of the N-terminal core domain exhibits micro-millisecond motions, a central β-hairpin shows elevated inherent flexibility on the pico-nanosecond timescale and the serine/arginine-rich region of residues 176-209 undergoes multiple exchange phenomena. Collectively, these studies have begun to reveal the complexities of the nucleocapsid protein dynamics and its preferred interaction sites with its biological targets.

How wide is the window opened by high-resolution relaxometry on the internal dynamics of proteins in solution?

The dynamics of molecules in solution is usually quantified by the determination of timescale-specific amplitudes of motions. High-resolution nuclear magnetic resonance (NMR) relaxometry experiments-where the sample is transferred to low fields for longitudinal (T₁) relaxation, and back to high field for detection with residue-specific resolution-seeks to increase the ability to distinguish the contributions from motion on timescales slower than a few nanoseconds. However, tumbling of a molecule in solution masks some of these motions. Therefore, we investigate to what extent relaxometry improves timescale resolution, using the "detector" analysis of dynamics. Here, we demonstrate improvements in the characterization of internal dynamics of methyl-bearing side chains by carbon-13 relaxometry in the small protein ubiquitin. We show that relaxometry data leads to better information about nanosecond motions as compared to high-field relaxation data only. Our calculations show that gains from relaxometry are greater with increasing correlation time of rotational diffusion.

Theoretical and computational framework for the analysis of the relaxation properties of arbitrary spin systems. Application to high-resolution relaxometry

A wide variety of nuclear magnetic resonance experiments rely on the prediction and analysis of relaxation processes. Recently, innovative approaches have been introduced where the sample travels through a broad range of magnetic fields in the course of the experiment, such as dissolution dynamic nuclear polarization or high-resolution relaxometry. Understanding the relaxation properties of nuclear spin systems over orders of magnitude of magnetic fields is essential to rationalize the results of these experiments. For example, during a high-resolution relaxometry experiment, the absence of control of nuclear spin relaxation pathways during the sample transfers and relaxation delays leads to systematic deviations of polarization decays from an ideal mono-exponential decay with the pure longitudinal relaxation rate. These deviations have to be taken into account to describe quantitatively the dynamics of the system. Here, we present computational tools to (1) calculate analytical expressions of relaxation rates for a broad variety of spin systems and (2) use these analytical expressions to correct the deviations arising in high-resolution relaxometry experiments. These tools lead to a better understanding of nuclear spin relaxation, which is required to improve the sensitivity of many pulse sequences, and to better characterize motions in macromolecules.

Boosting the resolution of low-field [Formula: see text] relaxation experiments on intrinsically disordered proteins with triple-resonance NMR

Improving our understanding of nanosecond motions in disordered proteins requires the enhanced sampling of the spectral density function obtained from relaxation at low magnetic fields. High-resolution relaxometry and two-field NMR measurements of relaxation have, so far, only been based on the recording of one- or two-dimensional spectra, which provide insufficient resolution for challenging disordered proteins. Here, we introduce a 3D-HNCO-based two-field NMR experiment for measurements of protein backbone [Formula: see text] amide longitudinal relaxation rates. The experiment provides accurate longitudinal relaxation rates at low field (0.33 T in our case) preserving the resolution and sensitivity typical for high-field NMR spectroscopy. Radiofrequency pulses applied on six different radiofrequency channels are used to manipulate the spin system at both fields. The experiment was demonstrated on the C-terminal domain of [Formula: see text] subunit of RNA polymerase from Bacillus subtilis, a protein with highly repetitive amino-acid sequence and very low dispersion of backbone chemical shifts.

Protein Dynamics from Accurate Low-Field Site-Specific Longitudinal and Transverse Nuclear Spin Relaxation

Nuclear magnetic relaxation provides invaluable quantitative site-specific information on the dynamics of complex systems. Determining dynamics on nanosecond time scales requires relaxation measurements at low magnetic fields incompatible with high-resolution NMR. Here, we use a two-field NMR spectrometer to measure carbon-13 transverse and longitudinal relaxation rates at a field as low as 0.33 T (proton Larmor frequency 14 MHz) in specifically labeled side chains of the protein ubiquitin. The use of radiofrequency pulses enhances the accuracy of measurements as compared to high-resolution relaxometry approaches, where the sample is moved in the stray field of the superconducting magnet. Importantly, we demonstrate that accurate measurements at a single low magnetic field provide enough information to characterize complex motions on low nanosecond time scales, which opens a new window for the determination of site-specific nanosecond motions in complex systems such as proteins.

High sensitivity high-resolution full range relaxometry using a fast mechanical sample shuttling device and a cryo-probe

Field-dependent NMR studies of bio-molecular systems using a sample shuttling hardware operating on a high-field NMR apparatus have provided valuable structural and dynamic information. We have recently published a design of a compact sample transportation device, called "field-cycler", which was installed in a commercial spectrometer and which provided highly precise positioning and stability during high speed shuttling. In this communication, we demonstrate the first use of a sample shuttling device on a commercial high field standard bore NMR spectrometer, equipped with a commercial triple resonance cryogenically cooled NMR probe. The performance and robustness of the hardware operating in 1D and 2D field cycling experiments, as well as the impact of the sample shuttling time on the signal intensity are discussed.

Substrate and Cofactor Dynamics on Guanosine Monophosphate Reductase Probed by High Resolution Field Cycling 31P NMR Relaxometry

Guanosine-5'-monophosphate reductase (GMPR) catalyzes the reduction of GMP to IMP and ammonia with concomitant oxidation of NADPH. Here we investigated the structure and dynamics of enzyme-bound substrates and cofactors by measuring 31P relaxation rates over a large magnetic field range using high resolution field cycling NMR relaxometry. Surprisingly, these experiments reveal differences in the low field relaxation profiles for the monophosphate of GMP compared with IMP in their respective NADP+ complexes. These complexes undergo partial reactions that mimic different steps in the overall catalytic cycle. The relaxation profiles indicate that the substrate monophosphates have distinct interactions in E·IMP·NADP+ and E·GMP·NADP+ complexes. These findings were not anticipated by x-ray crystal structures, which show identical interactions for the monophosphates of GMP and IMP in several inert complexes. In addition, the motion of the cofactor is enhanced in the E·GMP·NADP+ complex. Last, the motions of the substrate and cofactor are coordinately regulated; the cofactor has faster local motions than GMP in the deamination complex but is more constrained than IMP in that complex, leading to hydride transfer. These results show that field cycling can be used to investigate the dynamics of protein-bound ligands and provide new insights into how portions of the substrate remote from the site of chemical transformation promote catalysis.

Structural characterization of the N-terminal part of the MERS-CoV nucleocapsid by X-ray diffraction and small-angle X-ray scattering

The N protein of coronaviruses is a multifunctional protein that is organized into several domains. The N‐terminal part is composed of an intrinsically disordered region (IDR) followed by a structured domain called the N‐terminal domain (NTD). In this study, the structure determination of the N‐terminal region of the MERS‐CoV N protein via X‐ray diffraction measurements is reported at a resolution of 2.4 Å. Since the first 30 amino acids were not resolved by X‐ray diffraction, the structural study was completed by a SAXS experiment to propose a structural model including the IDR. This model presents the N‐terminal region of the MERS‐CoV as a monomer that displays structural features in common with other coronavirus NTDs.

Protein dynamics from nuclear magnetic relaxation

Protein dynamics are explored by a variety of methods designed to measure nuclear magnetic relaxation rates.

High-resolution two-field nuclear magnetic resonance spectroscopy

Two-field NMR provides correlations of nuclear spins at the most favourable magnetic fields in a single experiment.

The SARS coronavirus nucleocapsid protein--forms and functions

•

Coronavirus N proteins share the same modular organization.

•

Structures of SARS-CoV N protein provide insight into nucleocapsid formation.

•

N protein binds to nucleic acid at multiple sites in a coupled-allostery manner.

•

A RNP packaging model highlighting the importance of disorder and modularity is proposed.

The nucleocapsid phosphoprotein of the severe acute respiratory syndrome coronavirus (SARS-CoV N protein) packages the viral genome into a helical ribonucleocapsid (RNP) and plays a fundamental role during viral self-assembly. It is a protein with multifarious activities. In this article we will review our current understanding of the N protein structure and its interaction with nucleic acid. Highlights of the progresses include uncovering the modular organization, determining the structures of the structural domains, realizing the roles of protein disorder in protein–protein and protein–nucleic acid interactions, and visualizing the ribonucleoprotein (RNP) structure inside the virions. It was also demonstrated that N-protein binds to nucleic acid at multiple sites with a coupled-allostery manner. We propose a SARS-CoV RNP model that conforms to existing data and bears resemblance to the existing RNP structures of RNA viruses. The model highlights the critical role of modular organization and intrinsic disorder of the N protein in the formation and functions of the dynamic RNP capsid in RNA viruses. This paper forms part of a symposium in Antiviral Research on “From SARS to MERS: 10 years of research on highly pathogenic human coronaviruses.”

Cytotoxic amphiphiles and phosphoinositides bind to two discrete sites on the Akt1 PH domain

The mechanism of binding of two promising anticancer agents (the cytotoxic alkylphospholipids perifosine and miltefosine) to the Akt PH domain is investigated by high-resolution field-cycling (31)P nuclear magnetic resonance (NMR) spectroscopy using a spin-labeled recombinant PH domain. These results strongly indicate that there are two discrete amphiphile binding sites on the domain: (i) the cationic site that binds phosphoinositides and some alkylphospholipids and (ii) a second site that is occupied by only the alkylphospholipids. The identification of this second site for amphiphiles on the Akt1 PH domain provides a new target for drug development as well as insights into the regulation of the activity of the intact Akt1 protein. The field-cycling NMR methodology could be used to define discrete phospholipid or amphiphile binding sites on a wide variety of peripheral membrane proteins.

Nanosecond time scale motions in proteins revealed by high-resolution NMR relaxometry

Understanding the molecular determinants underlying protein function requires the characterization of both structure and dynamics at atomic resolution. Nuclear relaxation rates allow a precise characterization of protein dynamics at the Larmor frequencies of spins. This usually limits the sampling of motions to a narrow range of frequencies corresponding to high magnetic fields. At lower fields one cannot achieve sufficient sensitivity and resolution in NMR. Here, we use a fast shuttle device where the polarization builds up and the signals are detected at high field, while longitudinal relaxation takes place at low fields 0.5 < B₀ < 14.1 T. The sample is propelled over a distance up to 50 cm by a blowgun-like system in about 50 ms. The analysis of nitrogen-15 relaxation in the protein ubiquitin over such a wide range of magnetic fields offers unprecedented insights into molecular dynamics. Some key regions of the protein feature structural fluctuations on nanosecond time scales, which have so far been overlooked in high-field relaxation studies. Nanosecond motions in proteins may have been underestimated by traditional high-field approaches, and slower supra-τ_c motions that have no effect on relaxation may have been overestimated. High-resolution relaxometry thus opens the way to a quantitative characterization of nanosecond motions in proteins.

High-resolution NMR field-cycling device for full-range relaxation and structural studies of biopolymers on a shared commercial instrument

Improvements are described in a shuttling field-cycling device (Redfield in Magn Reson Chem 41:753-768, 2003), designed to allow widespread access to this useful technique by configuring it as a removable module to a commercial 500 MHz NMR instrument. The main improvements described here, leading to greater versatility, high reliability and simple construction, include: shuttling provided by a linear motor driven by an integrated-control servomotor; provision of automated bucking magnets to allow fast two-stage cycling to nearly zero field; and overall control by a microprocessor. A brief review of history and publications that have used the system is followed by a discussion of topics related to such a device including discussion of some future applications. A description of new aspects of the shuttling device follows. The minimum round trip time to 1T and above is less than 0.25 s and to 0.002 T is 0.36 s. Commercial probes are used and sensitivity is that of the host spectrometer reduced only by relaxation during travel. A key element is development of a linkage that prevents vibration of the linear motor from reaching the probe.

Functional transcriptional regulatory sequence (TRS) RNA binding and helix destabilizing determinants of murine hepatitis virus (MHV) nucleocapsid (N) protein

Coronavirus (CoV) nucleocapsid (N) protein contains two structurally independent RNA binding domains. These are denoted N-terminal domain (NTD) and C-terminal domain and are joined by a charged linker region rich in serine and arginine residues (SR linker). In mouse hepatitis virus (MHV), the NTD binds the transcriptional regulatory sequence (TRS) RNA, a conserved hexanucleotide sequence required for subgenomic RNA synthesis. The NTD is also capable of disrupting a short RNA duplex. We show here that three residues on the β3 (Arg-125 and Tyr-127) and β5 (Tyr-190) strands play key roles in TRS RNA binding and helix destabilization with Ala substitutions of these residues lethal to the virus. NMR studies of the MHV NTD·TRS complex revealed that this region defines a major RNA binding interface in MHV with site-directed spin labeling studies consistent with a model in which the adenosine-rich 3'-region of TRS is anchored by Arg-125, Tyr-127, and Tyr-190 in a way that is critical for efficient subgenomic RNA synthesis in MHV. Characterization of CoV N NTDs from infectious bronchitis virus and from severe acute respiratory syndrome CoV revealed that, although detailed NTD-TRS determinants are distinct from those of MHV NTD, rapid helix destabilization activity of CoV N NTDs is most strongly correlated with CoV function and virus viability.

A compact high-speed mechanical sample shuttle for field-dependent high-resolution solution NMR

Analysis of NMR relaxation data has provided significant insight on molecular dynamic, leading to a more comprehensive understanding of macromolecular functions. However, traditional methodology allows relaxation measurements performed only at a few fixed high fields, thus severely restricting their potential for extracting more complete dynamic information. Here we report the design and performance of a compact high-speed servo-mechanical shuttle assembly adapted to a commercial 600 MHz high-field superconducting magnet. The assembly is capable of shuttling the sample in a regular NMR tube from the center of the magnet to the top (fringe field ∼0.01 T) in 100 ms with no loss of sensitivity other than that due to intrinsic relaxation. The shuttle device can be installed by a single experienced user in 30 min. Excellent 2D-(15)N-HSQC spectra of (u-(13)C, (15)N)-ubiquitin with relaxation at low fields (3.77 T) and detection at 14.1T were obtained to illustrate its utility in R(1) measurements of macromolecules at low fields. Field-dependent (13)C-R(1) data of (3,3,3-d)-alanine at various field strengths were determined and analyzed to assess CSA and (1)H-(13)C dipolar contributions to the carboxyl (13)C-R(1).

Dynamics of biomolecules from picoseconds to seconds at atomic resolution

Although biomolecular dynamics has been investigated using NMR for at least 40 years, only in the past 20 years have internal motions been characterized at atomic resolution throughout proteins and nucleic acids. This development was made possible by multidimensional heteronuclear NMR approaches that provide near complete sequential signal assignments of uniformly labeled biomolecules. Recent methodological advances have enabled characterization of internal dynamics on timescales ranging from picoseconds to seconds, both in solution and in the solid state. The size, complexity and functional significance of biomolecules investigated by NMR continue to grow, as do the insights that have been obtained about function. In this article I review a number of recent advances that have made such studies possible, and provide a few examples of where NMR either by itself or in combination with other approaches has paved the way to a better understanding of the complex relationship between dynamics and biomolecular function. Finally, I discuss prospects for further advances in this field.

Electron-nuclear interactions as probes of domain motion in proteins

Long range interactions between nuclear spins and paramagnetic ions can serve as a sensitive monitor of internal motion of various parts of proteins, including functional loops and separate domains. In the case of interdomain motion, the interactions between the ion and NMR-observable nuclei are modulated in direction and magnitude mainly by a combination of overall and interdomain motions. The effects on observable parameters such as paramagnetic relaxation enhancement (PRE) and pseudocontact shift (PCS) can, in principle, be used to characterize motion. These parameters are frequently used for the purpose of structural refinements. However, their use to probe actual domain motions is less common and is lacking a proper theoretical treatment from a motional perspective. In this work, a suitable spin Hamiltonian is incorporated in a two body diffusion model to produce the time correlation function for the nuclear spin-paramagnetic ion interactions. Simulated observables for nuclei in different positions with respect to the paramagnetic ion are produced. Based on these simulations, it demonstrated that both the PRE and the PCS can be very sensitive probes of domain motion. Results for different nuclei within the protein sense different aspects of the motions. Some are more sensitive to the amplitude of the internal motion, others are more sensitive to overall diffusion rates, allowing separation of these contributions. Experimentally, the interaction strength can also be tuned by substitution of different paramagnetic ions or by varying magnetic field strength (in the case of lanthanides) to allow the use of more detailed diffusion models without reducing the reliability of data fitting.