Probing Protein Dynamics Using Multifield Variable Temperature NMR Relaxation and Molecular Dynamics Simulation

Probing the Dynamics of Yersinia Adhesin A (YadA) in Outer Membranes Hints at Requirements for β-Barrel Membrane Insertion

The vast majority of cells are protected and functionalized by a dense surface layer of glycans, proteoglycans, and glycolipids. This surface represents an underexplored space in structural biology that is exceedingly challenging to recreate in vitro. Here, we investigate β-barrel protein dynamics within an asymmetric outer membrane environment, with the trimeric autotransporter Yersinia adhesin A (YadA) as an example. Magic-angle spinning NMR relaxation data and a model-free approach reveal increased mobility in the second half of strand β2 after the conserved G72, which is responsible for membrane insertion and autotransport, and in the subsequent loop toward β3. In contrast, the protomer-protomer interaction sites (β1i-β4i-1) are rigid. Intriguingly, the mobility in the β-strand section following G72 is substantially elevated in the outer membrane and less so in the detergent environment of microcrystals. A possible source is revealed by molecular dynamics simulations that show the formation of a salt bridge involving E79 and R76 in competition with a dynamic interplay of calcium binding by E79 and the phosphate groups of the lipids. An estimation of overall barrel motion in the outer membrane and detergent-containing crystals yields values of around 41 ns for both. The global motion of YadA in the outer membrane has a stronger rotational component orthogonal to the symmetry axis of the trimeric porin than in the detergent-containing crystal. In summary, our investigation shows that the mobility in the second half of β2 and the loop to β3 required for membrane insertion and autotransport is maintained in the final folded form of YadA.

Nuclear spin diffusion under fast magic-angle spinning in solid-state NMR

Solid-state nuclear spin diffusion is the coherent and reversible process through which spin order is transferred via dipolar couplings. With the recent increases in magic-angle spinning (MAS) frequencies and magnetic fields becoming routinely applied in solid-state nuclear magnetic resonance, understanding how the increased 1H resolution obtained affects spin diffusion is necessary for interpretation of several common experiments. To investigate the coherent contributions to spin diffusion with fast MAS, we have developed a low-order correlation in Liouville space model based on the work of Dumez et al. (J. Chem. Phys. 33, 224501, 2010). Specifically, we introduce a new method for basis set selection, which accounts for the resonance-offset dependence at fast MAS. Furthermore, we consider the necessity of including chemical shift, both isotropic and anisotropic, in the modeling of spin diffusion. Using this model, we explore how different experimental factors change the nature of spin diffusion. Then, we show case studies to exemplify the issues that arise in using spin diffusion techniques at fast spinning. We show that the efficiency of polarization transfer via spin diffusion occurring within a deuterated and 100% back-exchanged protein sample at 60 kHz MAS is almost entirely dependent on resonance offset. We additionally identify temperature-dependent magnetization transfer in beta-aspartyl L-alanine, which could be explained by the influence of an incoherent relaxation-based nuclear Overhauser effect.

Slice & Dice: nested spin-lattice relaxation measurements

Spin-lattice relaxation rate (R₁) measurements are commonly used to characterize protein dynamics. However, the time needed to collect the data can be quite long due to long relaxation times of the low-gamma nuclei, especially in the solid state. We present a method to collect backbone heavy atom relaxation data by nesting the collection of datasets in the solid state. This method results in a factor of 2 to 2.5 times faster data acquisition for backbone R₁ relaxation data for the ¹³C and ¹⁵N sites of proteins.

1H-Detected Biomolecular NMR under Fast Magic-Angle Spinning

Since the first pioneering studies on small deuterated peptides dating more than 20 years ago, 1H detection has evolved into the most efficient approach for investigation of biomolecular structure, dynamics, and interactions by solid-state NMR. The development of faster and faster magic-angle spinning (MAS) rates (up to 150 kHz today) at ultrahigh magnetic fields has triggered a real revolution in the field. This new spinning regime reduces the 1H-1H dipolar couplings, so that a direct detection of 1H signals, for long impossible without proton dilution, has become possible at high resolution. The switch from the traditional MAS NMR approaches with 13C and 15N detection to 1H boosts the signal by more than an order of magnitude, accelerating the site-specific analysis and opening the way to more complex immobilized biological systems of higher molecular weight and available in limited amounts. This paper reviews the concepts underlying this recent leap forward in sensitivity and resolution, presents a detailed description of the experimental aspects of acquisition of multidimensional correlation spectra with fast MAS, and summarizes the most successful strategies for the assignment of the resonances and for the elucidation of protein structure and conformational dynamics. It finally outlines the many examples where 1H-detected MAS NMR has contributed to the detailed characterization of a variety of crystalline and noncrystalline biomolecular targets involved in biological processes ranging from catalysis through drug binding, viral infectivity, amyloid fibril formation, to transport across lipid membranes.

Dynamic Nuclear Polarization for Sensitivity Enhancement in Biomolecular Solid-State NMR

Solid-state NMR with magic-angle spinning (MAS) is an important method in structural biology. While NMR can provide invaluable information about local geometry on an atomic scale even for large biomolecular assemblies lacking long-range order, it is often limited by low sensitivity due to small nuclear spin polarization in thermal equilibrium. Dynamic nuclear polarization (DNP) has evolved during the last decades to become a powerful method capable of increasing this sensitivity by two to three orders of magnitude, thereby reducing the valuable experimental time from weeks or months to just hours or days; in many cases, this allows experiments that would be otherwise completely unfeasible. In this review, we give an overview of the developments that have opened the field for DNP-enhanced biomolecular solid-state NMR including state-of-the-art applications at fast MAS and high magnetic field. We present DNP mechanisms, polarizing agents, and sample constitution methods suitable for biomolecules. A wide field of biomolecular NMR applications is covered including membrane proteins, amyloid fibrils, large biomolecular assemblies, and biomaterials. Finally, we present perspectives and recent developments that may shape the field of biomolecular DNP in the future.

Dipolar Order Parameters in Large Systems With Fast Spinning

Order parameters are a useful tool for quantifying amplitudes of molecular motions. Here we measure dipolar order parameters by recoupling heteronuclear dipole-dipole couplings under fast spinning. We apply symmetry based recoupling methods to samples spinning under magic angle at 60 kHz by employing a variable flip angle compound inversion pulse. We validate the methods by measuring site-specific ¹⁵N-¹H order parameters of a microcrystalline protein over a small temperature range and the same protein in a large, precipitated complex with antibody. The measurements of the order parameters in the complex are consistent with the observed protein undergoing overall motion within the assembly.

Accelerating 15N and 13C R1 and R1ρ relaxation measurements by multiple pathway solid-state NMR experiments

Magic angle spinning (MAS) Solid-state NMR is a powerful technique to probe dynamics of biological systems at atomic resolution. R₁ and R_1ρ relaxation measurements can provide detailed insight on amplitudes and time scales of motions, especially when information from several different site-specific types of probes is combined. However, such experiments are time-consuming to perform. Shortening the time necessary to record relaxation data for different nuclei will greatly enhance practicality of such approaches. Here, we present staggered acquisition experiments to acquire multiple relaxation experiments from a single excitation to reduce the overall experimental time. Our strategy enables one to collect ¹⁵N and ¹³C relaxation data in a single experiment in a fraction of the time necessary for two separate experiments, with the same signal to noise ratio.

Observation of Sub-Microsecond Protein Methyl-Side Chain Dynamics by Nanoparticle-Assisted NMR Spin Relaxation

Amino-acid side-chain properties in proteins are key determinants of protein function. NMR spin relaxation of side chains is an important source of information about local protein dynamics and flexibility. However, traditional solution NMR relaxation methods are most sensitive to sub-nanosecond dynamics lacking information on slower ns-μs time-scale motions. Nanoparticle-assisted NMR spin relaxation (NASR) of methyl-side chains is introduced here as a window into these ns-μs dynamics. NASR utilizes the transient and nonspecific interactions between folded proteins and slowly tumbling spherical nanoparticles (NPs), whereby the increase of the relaxation rates reflects motions on time scales from ps all the way to the overall tumbling correlation time of the NPs ranging from hundreds of ns to μs. The observed motional amplitude of each methyl group can then be expressed by a model-free NASR S² order parameter. The method is demonstrated for ²H-relaxation of CH₂D methyl moieties and cross-correlated relaxation of CH₃ groups for proteins Im7 and ubiquitin in the presence of anionic silica-nanoparticles. Both types of relaxation experiments, dominated by either quadrupolar or dipolar interactions, yield highly consistent results. Im7 shows additional dynamics on the intermediate time scales taking place in a functionally important loop, whereas ubiquitin visits the majority of its conformational substates on the sub-ns time scale. These experimental observations are in good agreement with 4-10 μs all-atom molecular dynamics trajectories. NASR probes side-chain dynamics on a much wider range of motional time scales than previously possible, thereby providing new insights into the interplay between protein structure, dynamics, and molecular interactions that govern protein function.

Dynamics of uniformly labelled solid proteins between 100 and 300 K: A 2D 2H-13C MAS NMR approach

We describe a ²H based MAS nuclear magnetic resonance (NMR) method to obtain site-specific molecular dynamics of biomolecules. The method utilizes the use of deuterium nucleus as a spin label that is proven to be very useful in dynamics studies of solid biological and functional materials. The aim is to understand overall characteristics of protein backbone and side-chain motions for CD₃, CD₂ and CD groups, in terms of timescale, type and activation energy of the underlying processes. Variable temperature two-dimensional (2D) ²H-¹³C correlation MAS NMR spectra were recorded for the uniformly ²H,¹³C,¹⁵N labelled Alanine and microcrystalline SH3 at a broad temperature range, from 320 K down to 100 K. First, the deuterium quadrupolar-coupling constant from specific D-C sites is obtained with the 2D experiment by utilizing carbon chemical shifts. Second, the static quadrupolar patterns are obtained at 100 K. Third, variable temperature approach enabled the observation of quadrupolar pattern over different motional regimes; slow, intermediate and fast. And finally, the apparent activation energies for C-D sites are determined and compared, by evaluating the temperature induced signal intensities. This information led to the determination of the dynamic processes for different D-C sites at a broad range of temperature and motional timescales. This is a first representation of 2D ²H-¹³C MAS NMR approach applied to fully isotope labelled deuterated protein covering 220 K temperature range.

Protein structural dynamics by Magic-Angle Spinning NMR

Magic-Angle Spinning (MAS) Nuclear Magnetic Resonance (NMR) is a fast-developing technique, capable of complementing solution NMR, X-ray crystallography, and electron microscopy for the biophysical characterization of microcrystalline, poorly crystalline or disordered protein samples, such as enzymes, biomolecular assemblies, membrane-embedded systems or fibrils. Beyond structures, MAS NMR is an ideal tool for the investigation of dynamics, since it is unique in its ability to distinguish static and dynamic disorder, and to characterize not only amplitudes but also timescales of motion. Building on seminal work on model proteins, the technique is now ripe for widespread application in structural biology. This review briefly summarizes the recent evolutions in biomolecular MAS NMR and accounts for the growing number of systems where this spectroscopy has provided a description of conformational dynamics over the very last few years.

Atomic Resolution Map of Hierarchical Self-Assembly for an Amyloidogenic Protein Probed through Thermal 15N-R2 Correlation Matrices

Soluble oligomers formed by amyloidogenic intrinsically disordered proteins are some of the most cytotoxic species linked to neurodegeneration. Due to the transient and heterogeneous nature of such oligomeric intermediates, the underlying self-association events often remain elusive. NMR relaxation measurements sensitive to zero-frequency spectral densities (J(0)), such as the ¹⁵N - R₂ rates, are ideally suited to map sites of self-association at atomic resolution without the need of exogenous labels. Such experiments exploit the dynamic exchange between NMR visible monomers and slowly tumbling oligomers. However,¹⁵N - R₂ rates are also sensitive to intrinsic monomer dynamics, and it is often difficult to discern these contributions from those arising from exchange with oligomers. Another challenge pertains to defining a hierarchy of self-association. Here, using the archetypical amyloidogenic protein alpha synuclein (αS), we show that the temperature-dependence of ¹⁵N - R₂ effectively identifies self-association sites with reduced bias from internal dynamics. The key signature of the residues involved in self-association is a nonlinear temperature-dependence of ¹⁵N - R₂ with a positive ΔR₂/ΔT slope. These two hallmarks are systematically probed through a thermal R₂ correlation matrix, from which the network of residues involved in self-association as well as the hierarchy of αS self-association sites is extracted through agglomerative clustering. We find that aggregation is initiated by residues within the NAC region that is solvent inaccessible in αS fibrils and eventually extends to the N-terminal segment harboring familial PD mutations. These hierarchical self-association maps help dissect the essential drivers of oligomerization and reveal how amyloid inhibitors affect oligomer formation.

Similarities and Differences among Protein Dynamics Studied by Variable Temperature Nuclear Magnetic Resonance Relaxation

Understanding and describing the dynamics of proteins is one of the major challenges in biology. Here, we use multifield variable-temperature NMR longitudinal relaxation (R₁) measurements to determine the hierarchical activation energies of motions of four different proteins: two small globular proteins (GB1 and the SH3 domain of α-spectrin), an intrinsically disordered protein (the C-terminus of the nucleoprotein of the Sendai virus, Sendai Ntail), and an outer membrane protein (OmpG). The activation energies map the motions occurring in the side chains, in the backbone, and in the hydration shells of the proteins. We were able to identify similarities and differences in the average motions of the proteins. We find that the NMR relaxation properties of the four proteins do share similar features. The data characterizing average backbone motions are found to be very similar, the same for methyl group rotations, and similar activation energies are measured. The main observed difference occurs for the intrinsically disordered Sendai Ntail, where we observe much lower energy of activation for motions of protons associated with the protein-solvent interface as compared to the others. We also observe variability between the proteins regarding side chain ¹⁵N relaxation of lysine residues, with a higher activation energy observed in OmpG. This hints at strong interactions with negatively charged lipids in the bilayer and provides a possible mechanistic clue for the "positive-inside" rule for helical membrane proteins. Overall, these observations refine the understanding of the similarities and differences between hierarchical dynamics in proteins.

Structural Dynamics by NMR in the Solid State: The Unified MOMD Perspective Applied to Organic Frameworks with Interlocked Molecules

The microscopic-order-macroscopic-disorder (MOMD) approach for NMR lineshape analysis has been applied to the University of Windsor Dynamic Materials (UWDM) of types 1, 2, α-3, β-3, and 5, which are metal-organic frameworks (MOFs) comprising mobile mechanically interlocked molecules (MIMs). The mobile MIM components are selectively deuterated crown ether macrocycles - 24C6, 22C6, and B24C6. Their motion is described in MOMD by an effective/collective dynamic mode characterized by a diffusion tensor, R, a restricting/ordering potential, u, expanded in the Wigner rotation matrix elements, D_0, K^L, and features of local geometry. Experimental ²H lineshapes are available over 220 K (on average) and in some cases 320 K. They are reproduced with axial R, u given by the terms D_0,0² and D_0,|2|², and established local geometry. For UWDM of types 1, β-3, and 5, where the macrocycle resides in a relatively loose space, u is in the 1-3 kT, R_∥ in the (1.0-2.5) × 10⁶ s^-1, and R_⊥ in the (0.4-2.5) × 10⁴ s^-1 range; the deuterium atom is bonded to a carbon atom with tetrahedral coordination character. For UWDM of types 2 and α-3, where the macrocycle resides in a much tighter space, a substantial change in the symmetry of u and the coordination character of the ²H-bonded carbon are detected at higher temperatures. The activation energies for R_∥ and R_⊥ are characteristic of each system. The MOMD model is general; effective/collective dynamic modes are treated. The characteristics of motion, ordering, and geometry are physically well-defined; they differ from case to case in extent and symmetry but not in essence. Physical clarity and consistency provide new insights. A previous interpretation of the same experimental data used models consisting of collections of independent simple motions. These models are specific to each case and temperature. Within their scope, generating consistent physical pictures and comparing cases are difficult; possible collective modes are neglected.

77Se NMR Probes the Protein Environment of Selenomethionine

Sulfur is critical for the correct structure and proper function of proteins. Yet, lacking a sensitive enough isotope, nuclear magnetic resonance (NMR) experiments are unable to deliver for sulfur in proteins the usual wealth of chemical, dynamic, and structural information. This limitation can be circumvented by substituting sulfur with selenium, which has similar physicochemical properties and minimal impact on protein structures but possesses an NMR compatible isotope (⁷⁷Se). Here we exploit the sensitivity of ⁷⁷Se NMR to the nucleus' chemical milieu and use selenomethionine as a probe for its proteinaceous environment. However, such selenium NMR spectra of proteins currently resist a reliable interpretation because systematic connections between variations of system variables and changes in ⁷⁷Se NMR parameters are still lacking. To start narrowing this knowledge gap, we report here on a biological ⁷⁷Se magnetic resonance data bank based on a systematically designed library of GB1 variants in which a single selenomethionine was introduced at different locations within the protein. We recorded the resulting isotropic ⁷⁷Se chemical shifts and relaxation times for six GB1 variants by solution-state ⁷⁷Se NMR. For four of the GB1 variants we were also able to determine the chemical shift anisotropy tensor of SeM by solid-state ⁷⁷Se NMR. To enable interpretation of the NMR data, the structures of five of the GB1 variants were solved by X-ray crystallography to a resolution of 1.2 Å, allowing us to unambiguously determine the conformation of the selenomethionine. Finally, we combine our solution- and solid-state NMR data with the structural information to arrive at general insights regarding the execution and interpretation of ⁷⁷Se NMR experiments that exploit selenomethionine to probe proteins.

Aromatic Ring Dynamics, Thermal Activation, and Transient Conformations of a 468 kDa Enzyme by Specific 1H-13C Labeling and Fast Magic-Angle Spinning NMR

Aromatic residues are located at structurally important sites of many proteins. Probing their interactions and dynamics can provide important functional insight but is challenging in large proteins. Here, we introduce approaches to characterize the dynamics of phenylalanine residues using ¹H-detected fast magic-angle spinning (MAS) NMR combined with a tailored isotope-labeling scheme. Our approach yields isolated two-spin systems that are ideally suited for artifact-free dynamics measurements, and allows probing motions effectively without molecular weight limitations. The application to the TET2 enzyme assembly of ∼0.5 MDa size, the currently largest protein assigned by MAS NMR, provides insights into motions occurring on a wide range of time scales (picoseconds to milliseconds). We quantitatively probe ring-flip motions and show the temperature dependence by MAS NMR measurements down to 100 K. Interestingly, favorable line widths are observed down to 100 K, with potential implications for DNP NMR. Furthermore, we report the first ¹³C R_1ρ MAS NMR relaxation-dispersion measurements and detect structural excursions occurring on a microsecond time scale in the entry pore to the catalytic chamber and at a trimer interface that was proposed as the exit pore. We show that the labeling scheme with deuteration at ca. 50 kHz MAS provides superior resolution compared to 100 kHz MAS experiments with protonated, uniformly ¹³C-labeled samples.