Stereospecific assignments in proteins using exact NOEs

Advances in the exact nuclear Overhauser effect 2018-2022

The introduction of the exact nuclear Overhauser enhancement (eNOE) methodology to solution-state nuclear magnetic resonance (NMR) spectroscopy results in tighter distance restraints from NOEs than in convention analysis. These improved restraints allow for higher resolution in structure calculation and even the disentanglement of different conformations of macromolecules. While initial work primarily focused on technical development of the eNOE, structural studies aimed at the elucidation of spatial sampling in proteins and nucleic acids were published in parallel prior to 2018. The period of 2018-2022 saw a continued series of technical innovation, but also major applications addressing biological questions. Here, we review both aspects, covering topics from the implementation of non-uniform sampling of NOESY buildups, novel pulse sequences, adaption of the eNOE to solid-state NMR, advances in eNOE data analysis, and innovations in structural ensemble calculation, to applications to protein, RNA, and DNA structure elucidation.

Reducing the measurement time of exact NOEs by non-uniform sampling

We have previously reported on the measurement of exact NOEs (eNOEs), which yield a wealth of additional information in comparison to conventional NOEs. We have used these eNOEs in a variety of applications, including calculating high-resolution structures of proteins and RNA molecules. The collection of eNOEs is challenging, however, due to the need to measure a NOESY buildup series consisting of typically four NOESY spectra with varying mixing times in a single measurement session. While the 2D version can be completed in a few days, a fully sampled 3D-NOESY buildup series can take 10 days or more to acquire. This can be both expensive as well as problematic in the case of samples that are not stable over such a long period of time. One potential method to significantly decrease the required measurement time of eNOEs is to use non-uniform sampling (NUS) to decrease the number of points measured in the indirect dimensions. The effect of NUS on the extremely tight distance restraints extracted from eNOEs may be very pronounced. Therefore, we investigated the fidelity of eNOEs measured from three test cases at decreasing NUS densities: the 18.4 kDa protein human Pin1, the 4.1 kDa WW domain of Pin1 (both in 3D), and a 4.6 kDa 14mer RNA UUCG tetraloop (2D). Our results show that NUS imparted negligible error on the eNOE distances derived from good quality data down to 10% sampling for all three cases, but there is a noticeable decrease in the eNOE yield that is dependent upon the underlying sparsity, and thus complexity, of the sample. For Pin1, this transition occurred at roughly 40% while for the WW domain and the UUCG tetraloop it occurred at lower NUS densities of 20% and 10%, respectively. We rationalized these numbers through reconstruction simulations under various conditions. The extent of this loss depends upon the number of scans taken as well as the number of peaks to be reconstructed. Based on these findings, we have created guidelines for choosing an optimal NUS density depending on the number of peaks needed to be reconstructed in the densest region of a 2D or 3D NOESY spectrum.

Efficient Stereospecific Hβ2/3 NMR Assignment Strategy for Mid-Size Proteins

We present a strategy for stereospecific NMR assignment of H^β2 and H^β3 protons in mid-size proteins (~150 residues). For such proteins, resonance overlap in standard experiments is severe, thereby preventing unambiguous assignment of a large fraction of β-methylenes. To alleviate this limitation, assignment experiments may be run in high static fields, where higher decoupling power is required. Three-bond H^α–H^β J-couplings (³J_Hα–Hβ) are critical for stereospecific assignments of β-methylene protons, and for determining rotameric χ₁ states. Therefore, we modified a pulse sequence designed to measure accurate ³J_Hα–Hβ couplings such that probe heating was reduced, while the decoupling performance was improved. To further increase the resolution, we applied non-uniform sampling (NUS) schemes in the indirect ¹H and ¹³C dimensions. The approach was applied to two medium-sized proteins, odorant binding protein 22 (OBP22; 14.4 kDa) and Pin1 (18.2 kDa), at 900 MHz polarizing fields. The coupling values obtained from NUS and linear sampling were extremely well correlated. However, NUS decreased the overlap of H^β2/3 protons, thus supplying a higher yield of extracted ³J_Hα-Hβ coupling values when compared with linear sampling. A similar effect could be achieved with linear prediction applied to the linearly sampled data prior to the Fourier transformation. Finally, we used ³J_Hα–Hβ couplings from Pin1 in combination with either conventional or exact nuclear Overhauser enhancement (eNOE) restraints to determine the stereospecific assignments of β-methylene protons. The use of eNOEs further increased the fraction of unambiguously assigned resonances when compared with procedures using conventional NOEs.

eNORA2 Exact NOE Analysis Program

We have recently developed an NMR protocol to extract exact distances between nuclei in proteins from an exact interpretation of NOESY buildup intensities (eNOEs). This enabled us to calculate multistate structural ensembles that exhibit realistic spatial sampling and long-range correlations. Our initial studies were laborious and required a deep understanding of the underlying spin dynamics. Here, we present a MatLab package that integrates all data processing steps required to convert intensities of assigned peaks in NOESY series into upper and lower distance limits for structure calculation. Those steps include organization of the data in object format, extraction of autorelaxation and cross-relaxation rate constants by fitting of diagonal peak decays and cross peak buildups, validation of the data, correction for spin diffusion, graphical display of the results, and generation of distance limits in CYANA compatible format. The analysis may be carried out using a full relaxation matrix or a simplified "divide and conquer" approach that allows for partial deuteration of protons. As the program does not require expertise beyond that of standard resonance assignment/structure calculation, it is suitable for experts and nonexperts alike.

The Exact Nuclear Overhauser Enhancement: Recent Advances

Although often depicted as rigid structures, proteins are highly dynamic systems, whose motions are essential to their functions. Despite this, it is difficult to investigate protein dynamics due to the rapid timescale at which they sample their conformational space, leading most NMR-determined structures to represent only an averaged snapshot of the dynamic picture. While NMR relaxation measurements can help to determine local dynamics, it is difficult to detect translational or concerted motion, and only recently have significant advances been made to make it possible to acquire a more holistic representation of the dynamics and structural landscapes of proteins. Here, we briefly revisit our most recent progress in the theory and use of exact nuclear Overhauser enhancements (eNOEs) for the calculation of structural ensembles that describe their conformational space. New developments are primarily targeted at increasing the number and improving the quality of extracted eNOE distance restraints, such that the multi-state structure calculation can be applied to proteins of higher molecular weights. We then review the implications of the exact NOE to the protein dynamics and function of cyclophilin A and the WW domain of Pin1, and finally discuss our current research and future directions.

The Exact NOE as an Alternative in Ensemble Structure Determination

The structure-function paradigm is increasingly replaced by the structure-dynamics-function paradigm. All protein activity is steered by the interplay between enthalpy and entropy. Conformational dynamics serves as a proxy of conformational entropy. Therefore, it is essential to study not only the average conformation but also the spatial sampling of a protein on all timescales. To this purpose, we have established a protocol for determining multiple-state ensembles of proteins based on exact nuclear Overhauser effects (eNOEs). We have recently extended our previously reported eNOE data set for the protein GB3 by a very large set of backbone and side-chain residual dipolar couplings and three-bond J couplings. Here, we demonstrate that at least four structural states are required to represent the complete data set by dissecting the contributions to the CYANA target function, which quantifies restraint violations in structure calculation. We present a four-state ensemble of GB3, which largely preserves the characteristics obtained from eNOEs only. Due to the abundance of the input data, the ensemble and χ(1) angles in particular are well suited for cross-validation of the input data and comparison to x-ray structures. Principal component analysis is used to automatically identify and validate relevant states of the ensembles. Overall, our findings suggest that eNOEs are a valuable alternative to traditional NMR probes in spatial elucidation of proteins.

The Dynamic Basis for Signal Propagation in Human Pin1-WW

Allostery is the structural manifestation of information transduction in biomolecules. Its hallmark is conformational change induced by perturbations at a distal site. An increasing body of evidence demonstrates the presence of allostery in very flexible and even disordered proteins, encouraging a thermodynamic description of this phenomenon. Still, resolving such processes at atomic resolution is difficult. Here we establish a protocol to determine atomistic thermodynamic models of such systems using high-resolution solution state nuclear magnetic resonance data and extensive molecular simulations. Using this methodology, we study information transduction in the WW domain of a key cell-cycle regulator Pin1. Pin1 binds promiscuously to phospho-Ser/Thr-Pro motifs, however, disparate structural and dynamic responses have been reported upon binding different ligands. Our model consists of two topologically distinct states whose relative population may be specifically skewed by an incoming ligand. This model provides a canonical basis for the understanding of multi-functionality in Pin1.

Nuclear overhauser spectroscopy of chiral CHD methylene groups

Nuclear magnetic resonance spectroscopy (NMR) can provide a great deal of information about structure and dynamics of biomolecules. The quality of an NMR structure strongly depends on the number of experimental observables and on their accurate conversion into geometric restraints. When distance restraints are derived from nuclear Overhauser effect spectroscopy (NOESY), stereo-specific assignments of prochiral atoms can contribute significantly to the accuracy of NMR structures of proteins and nucleic acids. Here we introduce a series of NOESY-based pulse sequences that can assist in the assignment of chiral CHD methylene protons in random fractionally deuterated proteins. Partial deuteration suppresses spin-diffusion between the two protons of CH2 groups that normally impedes the distinction of cross-relaxation networks for these two protons in NOESY spectra. Three and four-dimensional spectra allow one to distinguish cross-relaxation pathways involving either of the two methylene protons so that one can obtain stereospecific assignments. In addition, the analysis provides a large number of stereospecific distance restraints. Non-uniform sampling was used to ensure optimal signal resolution in 4D spectra and reduce ambiguities of the assignments. Automatic assignment procedures were modified for efficient and accurate stereospecific assignments during automated structure calculations based on 3D spectra. The protocol was applied to calcium-loaded calbindin D9k. A large number of stereospecific assignments lead to a significant improvement of the accuracy of the structure.

The experimental accuracy of the uni-directional exact NOE

We have established protocols to calculate exact NOEs (eNOE) from NOE data. eNOEs lend unprecedented precision to the calculation of distance restraints used for structure calculation. Moreover, as eNOEs are averaged quantities over all conformations of a molecule, they may contain accessible information of the sampled conformational space. In practice, a prerequisite for an exact interpretation is the evaluation of both NOESY cross-peak buildups. For large molecular sizes, the fraction of NOEs which can only be obtained from one cross peak typically increases. Distance restraints derived from such NOEs must be used with a tolerance for errors associated with the broken symmetry of the individual magnetization transfer pathways. The correct choice of upper and lower limits is particularly important for multiple-state ensemble calculation, where too narrow tolerances may lead to incorrect spatial sampling. In order to dissect these pathways in heavy-atom resolved 3D NOESY experiments, we analyze 2D [(1)H, (1)H]-NOESY experiments, which are the fundamental building blocks of the former. In combination with an analysis of excitation and inversion profiles of pulses on heavy atoms and relaxation effects during HXQC elements, we derive a rule for the correct choice of upper and lower distance limits derived from such uni-directional NOEs. We show that normalization of the cross- to the diagonal-peak intensities of the spins of magnetization destination rather than origin leads to similar errors of the distance restraints. This opens up the prospect of extended collection of unidirectional eNOEs.

Compiled data set of exact NOE distance limits, residual dipolar couplings and scalar couplings for the protein GB3

We compiled an NMR data set consisting of exact nuclear Overhauser enhancement (eNOE) distance limits, residual dipolar couplings (RDCs) and scalar (J) couplings for GB3, which forms one of the largest and most diverse data set for structural characterization of a protein to date. All data have small experimental errors, which are carefully estimated. We use the data in the research article Vogeli et al., 2015, Complementarity and congruence between exact NOEs and traditional NMR probes for spatial decoding of protein dynamics, J. Struct. Biol., 191, 3, 306–317, doi:10.1016/j.jsb.2015.07.008 [1] for cross-validation in multiple-state structural ensemble calculation. We advocate this set to be an ideal test case for molecular dynamics simulations and structure calculations.

Combined automated NOE assignment and structure calculation with CYANA

The automated assignment of NOESY cross peaks has become a fundamental technique for NMR protein structure analysis. A widely used algorithm for this purpose is implemented in the program CYANA. It has been used for a large number of structure determinations of proteins in solution but was so far not described in full detail. In this paper we present a complete description of the CYANA implementation of automated NOESY assignment, which differs extensively from its predecessor CANDID by the use of a consistent probabilistic treatment, and we discuss its performance in the second round of the critical assessment of structure determination by NMR.

Towards a true protein movie: a perspective on the potential impact of the ensemble-based structure determination using exact NOEs

Confined by the Boltzmann distribution of the energies of the states, a multitude of structural states are inherent to biomolecules. For a detailed understanding of a protein's function, its entire structural landscape at atomic resolution and insight into the interconversion between all the structural states (i.e. dynamics) are required. Whereas dedicated trickery with NMR relaxation provides aspects of local dynamics, and 3D structure determination by NMR is well established, only recently have several attempts been made to formulate a more comprehensive description of the dynamics and the structural landscape of a protein. Here, a perspective is given on the use of exact NOEs (eNOEs) for the elucidation of structural ensembles of a protein describing the covered conformational space.