Multiple-state ensemble structure determination from eNOE spectroscopy

NMR tools to detect protein allostery

Allostery is a fundamental mechanism of cellular homeostasis by intra-protein communication between distinct functional sites. It is an internal process of proteins to steer interactions not only with each other but also with other biomolecules such as ligands, lipids, and nucleic acids. In addition, allosteric regulation is particularly important in enzymatic activities. A major challenge in structural and molecular biology today is unraveling allosteric sites in proteins, to elucidate the detailed mechanism of allostery and the development of allosteric drugs. Here we summarize the recently developed tools and approaches which enable the elucidation of regulatory hotspots and correlated motion in biomolecules, focusing primarily on solution-state nuclear magnetic resonance spectroscopy (NMR). These tools open an avenue towards a rational understanding of the mechanism of allostery and provide essential information for the design of allosteric drugs.

Differential Structural Features of Two Mutant ADAR1p150 Zα Domains Associated with Aicardi-Goutières Syndrome

The Zα domain of ADARp150 is critical for proper Z-RNA substrate binding and is a key factor in the type-I interferon response pathway. Two point-mutations in this domain (N173S and P193A), which cause neurodegenerative disorders, are linked to decreased A-to-I editing in disease models. To understand this phenomenon at the molecular level, we biophysically and structurally characterized these two mutated domains, revealing that they bind Z-RNA with a decreased affinity. Less efficient binding to Z-RNA can be explained by structural changes in beta-wing, part of the Z-RNA-protein interface, and alteration of conformational dynamics of the proteins.

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.

An improved, time-efficient approach to extract accurate distance restraints for NMR2 structure calculation

Exact nuclear Overhauser enhancement (eNOE) yields highly accurate, ensemble averaged 1 H-1 H distance restraints with an accuracy of up to 0.1 Å for the multi-state structure determination of proteins as well as for nuclear magnetic resonance molecular replacement (N MR2 ) to determine the structure of the protein-ligand interaction site in a time-efficient manner. However, in the latter application, the acquired eNOEs lack the obtainable precision of 0.1 Å because of the asymmetrical nature of the filtered nuclear Overhauser enhancement spectroscopy (NOESY) experiment used in N MR2 . This error is further propagated to the eNOE equations used to fit and extract the distance restraints. In this work, a new analysis method is proposed to obtain inter-molecular distance restraints from the filtered NOESY spectrum more accurately and intuitively by dividing the NOE cross peak by the corresponding diagonal peak of the ligand. The method termed diagonal-normalised eNOEs was tested on the data acquired by on the complex of PIN1 and a small, weak-binding phenylimidazole fragment. N MR2 calculations performed using the distances derived from diagonal-normalised eNOEs yielded the right orientation of the fragment in the binding pocket and produced a structure that more closely resembles the benchmark X-ray structure (2XP6) with an average heavy-atom root-mean-square deviation (RMSD) of 1.681 Å with respect to it, when compared to the one produced with traditional N MR2 with an average heavy atom RMSD of 3.628 Å. This is attributed to the higher precision of the evaluated distance restraints.

On the use of residual dipolar couplings in multi-state structure calculation of two-domain proteins

Residual dipolar couplings (RDCs) are powerful nuclear magnetic resonance (NMR) probes for the structure calculation of biomacromolecules. Typically, an alignment tensor that defines the orientation of the entire molecule relative to the magnetic field is determined either before refinement of individual bond vectors or simultaneously with this refinement. For single-domain proteins this approach works well since all bond vectors can be described within the same coordinate frame, which is given by the alignment tensor. However, novel approaches are sought after for systems where no universal alignment tensor can be used. Here, we present an approach that can be applied to two-domain proteins that enables the calculation of multiple states within each domain as well as with respect to the relative positions of the two domains.

Motional clustering in supra-τc conformational exchange influences NOE cross-relaxation rate

Biomolecular spin relaxation processes, such as the NOE, are commonly modeled by rotational τ_c-tumbling combined with fast motions on the sub-τ_c timescale. Motions on the supra-τ_c timescale, in contrast, are considered to be completely decorrelated to the molecular tumbling and therefore invisible. Here, we show how supra-τ_c dynamics can nonetheless influence the NOE build-up between methyl groups. This effect arises because supra-τ_c motions can cluster the fast-motion ensembles into discrete states, affecting distance averaging as well as the fast-motion order parameter and hence the cross-relaxation rate. We present a computational approach to estimate methyl-methyl cross-relaxation rates from extensive (>100×τ_c) all-atom molecular dynamics (MD) trajectories on the example of the 723-residue protein Malate Synthase G. The approach uses Markov state models (MSMs) to resolve transitions between metastable states and thus to discriminate between sub-τ_c and supra-τ_c conformational exchange. We find that supra-τ_c exchange typically increases NOESY cross-peak intensities. The methods described in this work extend the theory of modeling sub-μs dynamics in spin relaxation and thus contribute to a quantitative estimation of NOE cross-relaxation rates from MD simulations, eventually leading to increased precision in structural and functional studies of large proteins.

A B-factor for NOEs?

Nuclear Overhauser effects (NOEs) are influenced by motion. Here, we derive exact, analytical results for a model of isotropic, harmonic fluctuations of atom positions that corresponds to the one underlying crystallographic B-factors. The model includes steric repulsion and yields closed-form expressions for the expected value of general invertible functions of the distance between two atoms, with the special case r^-6 for NOEs. We discuss the implications for the definition of an NOE-based B-factor in solution NMR.

Reconstruction of Coupled Intra- and Interdomain Protein Motion from Nuclear and Electron Magnetic Resonance

Proteins composed of multiple domains allow for structural heterogeneity and interdomain dynamics that may be vital for function. Intradomain structures and dynamics can influence interdomain conformations and vice versa. However, no established structure determination method is currently available that can probe the coupling of these motions. The protein Pin1 contains separate regulatory and catalytic domains that sample "extended" and "compact" states, and ligand binding changes this equilibrium. Ligand binding and interdomain distance have been shown to impact the activity of Pin1, suggesting interdomain allostery. In order to characterize the conformational equilibrium of Pin1, we describe a novel method to model the coupling between intra- and interdomain dynamics at atomic resolution using multistate ensembles. The method uses time-averaged nuclear magnetic resonance (NMR) restraints and double electron-electron resonance (DEER) data that resolve distance distributions. While the intradomain calculation is primarily driven by exact nuclear Overhauser enhancements (eNOEs), J couplings, and residual dipolar couplings (RDCs), the relative domain distribution is driven by paramagnetic relaxation enhancement (PREs), RDCs, interdomain NOEs, and DEER. Our data support a 70:30 population of the compact and extended states in apo Pin1. A multistate ensemble describes these conformations simultaneously, with distinct conformational differences located in the interdomain interface stabilizing the compact or extended states. We also describe correlated conformations between the catalytic site and interdomain interface that may explain allostery driven by interdomain contact.

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.

Protein Motional Details Revealed by Complementary Structural Biology Techniques

Proteins depend on defined molecular plasticity for their functionality. How to comprehensively capture dynamics correctly is of ubiquitous biological importance. Approaches commonly used to probe protein dynamics include model-free elucidation of site-specific motion by NMR relaxation, molecular dynamics (MD)-based approaches, and capturing the substates within a dynamic ensemble by recent eNOE-based multiple-structure approaches. Even though MD is sometimes combined with ensemble-averaged NMR restraints, these approaches have largely been developed and used individually. Owing to the different underlying concepts and practical requirements, it has remained unclear how they compare, and how they cross-validate and complement each other. Here, we extract and compare the differential information contents of MD simulations, NMR relaxation measurements, and eNOE-based multi-state structures for the SH3 domain of chicken α-spectrin. The data show that a validated, consistent, and detailed picture is feasible both for timescales and actual conformational states sampled in the dynamic ensemble. This includes the biologically important side-chain plasticity, for which experimentally cross-validated assessment is a significant challenge.

Aedes aegypti Odorant Binding Protein 22 selectively binds fatty acids through a conformational change in its C-terminal tail

Aedes aegypti is the primary vector for transmission of Dengue, Zika and chikungunya viruses. Previously it was shown that Dengue virus infection of the mosquito led to an in increased expression of the odorant binding protein 22 (AeOBP22) within the mosquito salivary gland and that siRNA mediated knockdown of AeOBP22 led to reduced mosquito feeding behaviors. Insect OBPs are implicated in the perception, storage and transport of chemosensory signaling molecules including air-borne odorants and pheromones. AeOBP22 is unusual as it is additionally expressed in multiple tissues, including the antenna, the male reproductive glands and is transferred to females during reproduction, indicating multiple roles in the mosquito life cycle. However, it is unclear what role it plays in these tissues and what ligands it interacts with. Here we present solution and X-ray crystallographic studies that indicate a potential role of AeOBP22 binding to fatty acids, and that the specificity for longer chain fatty acids is regulated by a conformational change in the C-terminal tail that leads to creation of an enlarged binding cavity that enhances binding affinity. This study sheds light onto the native ligands for AeOBP22 and provides insight into its potential functions in different tissues.

Complete NMR chemical shift assignments of odorant binding protein 22 from the yellow fever mosquito, Aedes aegypti, bound to arachidonic acid

Aedes aegypti mosquitoes are the vector for transmission of Dengue, Zika and chikungunya viruses. These mosquitos feed exclusively on human hosts for a blood meal. Previous studies have established that Dengue virus infection of the mosquito results in increased expression of the odorant binding proteins 22 and 10 within the mosquito salivary gland and silencing of these genes dramatically reduces blood-feeding behaviors. Odorant binding proteins are implicated in modulating the chemosensory perception of external stimuli that regulate behaviors such as host location, feeding and reproduction. However, the role that AeOBP22 plays in the salivary gland is unclear. Here, as a first step to a more complete understanding of the function of AeOBP22, we present the complete backbone and side chain chemical shift assignments of the protein in the complex it forms with arachidonic acid. These assignments reveal that the protein consists of seven α-helices, and that the arachidonic acid is bound tightly to the protein. Comparison with the chemical shift assignments of the apo-form of the protein reveals that binding of the fatty acid is accompanied by a large conformational change in the C-terminal helix, which appears disordered in the absence of lipid. This NMR data provides the basis for determining the structure of AeOBP22 and understanding the nature of the conformational changes that occur upon ligand binding. This information will provide a path to discover novel compounds that can interfere with AeOBP22 function and impact blood feeding by this mosquito.

Extending the Applicability of Exact Nuclear Overhauser Enhancements to Large Proteins and RNA

Distance-dependent nuclear Overhauser enhancements (NOEs) are one of the most popular and important experimental restraints for calculating NMR structures. Despite this, they are mostly employed as semiquantitative upper distance bounds, and this discards the wealth of information that is encoded in the cross-relaxation rate constant. Information that is lost includes exact distances between protons and dynamics that occur on the sub-millisecond timescale. Our recently introduced exact measurement of the NOE (eNOE) requires little additional experimental effort relative to other NMR observables. So far, we have used eNOEs to calculate multistate ensembles of proteins up to approximately 150 residues. Here, we briefly revisit eNOE methodology and present two new directions for the use of eNOEs: applications to large proteins and RNA.

NOE-Derived Methyl Distances from a 360 kDa Proteasome Complex

Nuclear magnetic resonance spectroscopy is the prime tool to probe structure and dynamics of biomolecules at atomic resolution. A serious challenge for that method is the size limit imposed on molecules to be studied. Standard studies are typically restricted to ca. 30-40 kDa. More recent developments lead to spin relaxation measurements in methyl groups in single proteins or protein complexes as large as a mega-Dalton, which directly allow the extraction of angular information or experiments with paramagnetic samples. However, these probes are all of indirect nature in contrast to the most intuitive and easy-to-interpret structural/dynamics restraint, the internuclear distance, which can be measured by nuclear Overhauser enhancement (NOE). Herein, we demonstrate time-averaged distance measurements on the 360 kDa half proteasome from Thermoplasma acidophilium. The approach is based on exact quantification of the NOE (eNOE). Our findings open up an avenue for such measurements on very large molecules. These restraints will help in a detailed determination of conformational changes upon perturbation such as ligand binding.

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 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.

Extending the eNOE data set of large proteins by evaluation of NOEs with unresolved diagonals

The representation of a protein's spatial sampling at atomic resolution is fundamental for understanding its function. NMR has been established as the best-suited technique toward this goal for small proteins. However, the accessible information content rapidly deteriorates with increasing protein size. We have recently demonstrated that for small proteins distance restraints with an accuracy smaller than 0.1 Å can be obtained by replacing traditional semi-quantitative Nuclear Overhauser Effects (NOEs) with exact NOEs (eNOE). The high quality of the data allowed us to calculate structural ensembles of the small model protein GB3 consisting of multiple rather than a single state. The analysis has been limited to small proteins because NOEs of spins with unresolved diagonal peaks cannot be used. Here we propose a simple approach to translate such NOEs into correct upper distance restraints, which opens access to larger biomolecules. We demonstrate that for 16 kDa cyclophilin A the collection of such restraints extends the original 1254 eNOEs to 3471.

Solution NMR studies of recombinant Aβ(1-42): from the presence of a micellar entity to residual β-sheet structure in the soluble species

Amyloid-β (Aβ) peptide is the major component found in senile plaques of Alzheimer's disease patients. The 42-residue fragment Aβ(1-42) is proposed to be one of the most pathogenic species therein. Here, the soluble Aβ(1-42) species were analyzed by various liquid-state NMR methods. Transient formation of a micelle species was observed at the onset of the aggregation kinetics. This micelle is dissolved after approximately one day. Subsequent loss of this species and the formation of protofibrils are proposed to be the route of fibril formation. Consequently, the observed micelle species is suggested to be on an off-pathway mechanism. Furthermore, characterization of the NMR-observable soluble species shows that it is a random-coil-like entity with low propensities for four β-strands. These β-strands correlate with the β-strand segments observed in Aβ fibrils. This finding indicates that the 3D structure of the fibrils might already be predisposed in the soluble species.

The nuclear Overhauser effect from a quantitative perspective

The nuclear Overhauser enhancement or effect (NOE) is the most important measure in liquid-state NMR with macromolecules. Thus, the NOE is the subject of numerous reviews and books. Here, the NOE is revisited in light of our recently introduced measurements of exact nuclear Overhauser enhancements (eNOEs), which enabled the determination of multiple-state 3D protein structures. This review encompasses all relevant facets from the theoretical considerations to the use of eNOEs in multiple-state structure calculation. Important aspects include a detailed presentation of the relaxation theory relevant for the nuclear Overhauser effect, the estimation of the correction for spin diffusion, the experimental determination of the eNOEs, the conversion of eNOE rates into distances and validation of their quality, the distance-restraint classification and the protocols for calculation of structures and ensembles.

Stereospecific assignments in proteins using exact NOEs

Recently developed methods to measure distances in proteins with high accuracy by "exact" nuclear Overhauser effects (eNOEs) make it possible to determine stereospecific assignments, which are particularly important to fully exploit the accuracy of the eNOE distance measurements. Stereospecific assignments are determined by comparing the eNOE-derived distances to protein structure bundles calculated without stereospecific assignments, or an independently determined crystal structure. The absolute and relative CYANA target function difference upon swapping the stereospecific assignment of a diastereotopic group yields the respective stereospecific assignment. We applied the method to the eNOE data set that has recently been obtained for the third immunoglobulin-binding domain of protein G (GB3). The 884 eNOEs provide relevant data for 47 of the total of 75 diastereotopic groups. Stereospecific assignments could be established for 45 diastereotopic groups (96 %) using the X-ray structure, or for 27 diastereotopic groups (57 %) using structures calculated with the eNOE data set without stereospecific assignments, all of which are in agreement with those determined previously. The latter case is relevant for structure determinations based on eNOEs. The accuracy of the eNOE distance measurements is crucial for making stereospecific assignments because applying the same method to the traditional NOE data set for GB3 with imprecise upper distance bounds yields only 13 correct stereospecific assignments using the X-ray structure or 2 correct stereospecific assignments using NMR structures calculated without stereospecific assignments.