Optimization of resolution and sensitivity of 4D NOESY using multi-dimensional decomposition

Leveraging the HMBC to Facilitate Metabolite Identification

The accuracy and ease of metabolite assignments from a complex mixture are expected to be facilitated by employing a multispectral approach. The two-dimensional (2D) 1H-13C heteronuclear single quantum coherence (HSQC) and 2D 1H-1H-total correlation spectroscopy (TOCSY) are the experiments commonly used for metabolite assignments. The 2D 1H-13C HSQC-TOCSY and 2D 1H-13C heteronuclear multiple-bond correlation (HMBC) are routinely used by natural products chemists but have seen minimal usage in metabolomics despite the unique information, the nearly complete 1H-1H and 1H-13C and spin systems provided by these experiments that may improve the accuracy and reliability of metabolite assignments. The use of a 13C-labeled feedstock such as glucose is a routine practice in metabolomics to improve sensitivity and to emphasize the detection of specific metabolites but causes severe artifacts and an increase in spectral complexity in the HMBC experiment. To address this issue, the standard HMBC pulse sequence was modified to include carbon decoupling. Nonuniform sampling was also employed for rapid data collection. A dataset of reference 2D 1H-13C HMBC spectra was collected for 94 common metabolites. 13C-13C spin connectivity was then obtained by generating a covariance pseudo-spectrum from the carbon-decoupled HMBC and the 1H-13C HSQC-TOCSY spectra. The resulting 13C-13C pseudo-spectrum provides a connectivity map of the entire carbon backbone that uniquely describes each metabolite and would enable automated metabolite identification.

Enhanced Nuclear Magnetic Resonance Spectroscopy with Isotropic Mixing as a Pseudodimension

Chemical analysis based on liquid-state nuclear magnetic resonance spectroscopy exploits numerous observables, mainly chemical shifts, relaxation rates, and internuclear coupling constants. Regarding the latter, the efficiencies of internuclear coherence transfers may be encoded in spectral peak intensities. The dependencies of these intensities on the experimental parameter that influences the transfer, for example, mixing time, are an important source of structural information. Yet, they are costly to measure and difficult to analyze. Here, we show that peak intensity build-up curves in two-dimensional total correlation spectroscopy (2D TOCSY) experiments may be quickly measured by employing nonuniform sampling and that their analysis can be effective if supported by quantum mechanical calculations. Thus, such curves can be used to form a new, third pseudodimension of the TOCSY spectrum. Similarly to the other two frequency dimensions, this one also resolves ambiguities and provides characteristic information. We show how the approach supports the analysis of a fragment of protein Tau Repeat-4 domain. Yet, its potential applications are far broader, including the analysis of complex mixtures or other polymers.

Codon Harmonization of a Kir3.1-KirBac1.3 Chimera for Structural Study Optimization

The expression of functional, folded, and isotopically enriched membrane proteins is an enduring bottleneck for nuclear magnetic resonance (NMR) studies. Indeed, historically, protein yield optimization has been insufficient to allow NMR analysis of many complex Eukaryotic membrane proteins. However, recent work has found that manipulation of plasmid codons improves the odds of successful NMR-friendly protein production. In the last decade, numerous studies showed that matching codon usage patterns in recombinant gene sequences to those in the native sequence is positively correlated with increased protein yield. This phenomenon, dubbed codon harmonization, may be a powerful tool in optimizing recombinant expression of difficult-to-produce membrane proteins for structural studies. Here, we apply this technique to an inward rectifier K⁺ Channel (Kir) 3.1-KirBac1.3 chimera. Kir3.1 falls within the G protein-coupled inward rectifier K⁺ (GIRK) channel family, thus NMR studies may inform on the nuances of GIRK gating action in the presence and absence of its G Protein, lipid, and small molecule ligands. In our hands, harmonized plasmids increase protein yield nearly two-fold compared to the traditional ‘fully codon optimized’ construct. We then employ a fluorescence-based functional assay and solid-state NMR correlation spectroscopy to show the final protein product is folded and functional.

Comparison of backbone dynamics of the p50 dimerization domain of NFκB in the homodimeric transcription factor NFκB1 and in its heterodimeric complex with RelA (p65)

The nuclear factor of kappa light polypeptide gene enhancer in B‐cells (NFκB) transcription factors play a critical role in human immune response. The family includes homodimers and heterodimers of five component proteins, which mediate different transcriptional responses and bind preferentially to different DNA sequences. Crystal structures of DNA complexes show that the dimers of the Rel‐homology regions are structurally very similar. Differing DNA sequence preference together with structural similarity suggests that the dimers may differ in their dynamics. In this study, we present the first near‐complete ¹⁵N, ¹³C_α/β, and H_N backbone resonance assignments of two dimers of the dimerization domain (DD) of the NFκB1 (p50) protein (residues 241–351): the homodimer of two p50 domains and a heterodimer of the p50 DD with the p65 DD. As expected, the two dimers behave very similarly, with chemical shift differences between them largely concentrated in the dimer interface and attributable to specific differences in the amino acid sequences of p50 and p65. A comparison of the picosecond‐nanosecond dynamics of the homo‐ and heterodimers also shows that the environment of p50 is similar, with an overall slightly reduced correlation time for the homodimer compared to the heterodimer, consistent with its slightly smaller molecular weight. These results demonstrate that NMR spectroscopy can be used to explore subtle changes in structure and dynamics that have the potential to give insights into differences in specificity that can be exploited in the design of new therapeutic agents.

The Human Fecal Microbiota Metabolizes Foodborne Heterocyclic Aromatic Amines by Reuterin Conjugation and Further Transformations

SCOPE

Heterocyclic aromatic amines (HAAs) are process-induced food contaminants with high mutagenic and/or carcinogenic potential. Although the human gut microbiota is known to affect the metabolism of dietary constituents, its impact on HAA metabolism and toxicity has been little studied. Here, the glycerol-dependent metabolism of seven foodborne HAAs (AαC, Trp-P-1, harman, norharman, PhIP, MeIQx, and MeIQ) by the human fecal microbiota is investigated.

METHODS AND RESULTS

As analyzed by HPLC-DAD/FLD, the extent of conversion is strongly dependent on glycerol supplementation and HAA structure. AαC (60-100%) and the 2-aminoimidazoazarenes (up to 58%) are especially prone to microbial conversion. Based on high-resolution MS and/or NMR spectroscopy data, 70 fecal metabolites are identified in total, mainly formed by chemical reactions with one or two molecules of microbially derived reuterin. Moreover, it has been demonstrated that the human fecal microbiota can further transform reuterin adducts by reduction and/or hydroxylation reactions. Upon isolation, some reuterin-induced HAA metabolites appear to be partially unstable, complicating structural identification.

CONCLUSION

The formation of microbial metabolites needs to be incorporated into risk assessment considerations for HAAs in human health. In this study, several HAA metabolites, mainly reuterin-dependent, are identified in vitro, providing the basis for future human studies investigating microbial HAA metabolism.

Synthetic Biology-Based Solution NMR Studies on Membrane Proteins in Lipid Environments

Although membrane proteins are in the focus of biochemical research for many decades the general knowledge of this important class is far behind soluble proteins. Despite several recent technical developments, the most challenging feature still is the generation of high-quality samples in environments suitable for the selected application. Reconstitution of membrane proteins into lipid bilayers will generate the most native-like environment and is therefore commonly desired. However, it poses tremendous problems to solution-state NMR analysis due to the dramatic increase in particle size resulting in high rotational correlation times. Nevertheless, a few promising strategies for the solution NMR analysis of membrane inserted proteins are emerging and will be discussed in this chapter. We focus on the generation of membrane protein samples in nanodisc membranes by cell-free systems and will describe the characteristic advantages of that platform in providing tailored protein expression and folding environments. We indicate frequent problems that have to be overcome in cell-free synthesis, nanodisc preparation, and customization for samples dedicated for solution-state NMR. Detailed instructions for sample preparation are given, and solution NMR approaches suitable for membrane proteins in bilayers are compiled. We further discuss the current strategies applied for signal detection from such difficult samples and describe the type of information that can be extracted from the various experiments. In summary, a comprehensive guideline for the analysis of membrane proteins in native-like membrane environments by solution-state NMR techniques will be provided.

SOFAST-HMQC-an efficient tool for metabolomics

Nuclear magnetic resonance (NMR)-based metabolomics relies mostly on 1D NMR; however, the technique is limited by overlap of the signals from the metabolites. In order to circumvent this problem, 2D ¹H-¹³C correlation spectroscopy techniques are often used. However owing to poorer natural abundance and gyromagnetic ratio of ¹³C, the acquisition time for 2D ¹H-¹³C heteronuclear single quantum coherence spectroscopy (HSQC) is long. This makes it almost impossible to be used in high throughput study. We have reported the application of selective optimized flip angle short transient (SOFAST) technique coupled to heteronuclear multiple quantum correlation (HMQC) along with nonlinear sampling (NUS) in urine and serum samples. This technique takes sevenfold less experimental time than the conventional ¹H-¹³C HSQC experiment with retention of almost all molecular information. Hence, this can be used for high throughput study. Graphical abstract SOFAST-HMQC is a two-dimensional NMR technique that significantly decreases experimental time without loss of information. This technique is applied in complex biofluid samples that are used for high throughput metabolomics studies and shows promise of better information recovery than conventional two-dimensional NMR technique in shorter time.

Acceleration of protein backbone NMR assignment by combinatorial labeling: Application to a small molecule binding study

Selective labeling with stable isotopes has long been recognized as a valuable tool in protein NMR to alleviate signal overlap and sensitivity limitations. In this study, combinatorial 15 N-, 13 Cα -, and 13 C'-selective labeling has been used during the backbone assignment of human cyclophilin D to explore binding of an inhibitor molecule. Using a cell-free expression system, a scheme that involves 15 N, 1-13 C, 2-13 C, fully 15 N/13 C, and unlabeled amino acids was optimized to gain a maximum of assignment information from three samples. This scheme was combined with time-shared triple-resonance NMR experiments, which allows a fast and efficient backbone assignment by giving the unambiguous assignment of unique amino acid pairs in the protein, the identity of ambiguous pairs and information about all 19 non-proline amino acid types. It is therefore well suited for binding studies where de novo assignments of amide 1 H and 15 N resonances need to be obtained, even in cases where sensitivity is the limiting factor.

Automated assignment of NMR chemical shifts based on a known structure and 4D spectra

Apart from their central role during 3D structure determination of proteins the backbone chemical shift assignment is the basis for a number of applications, like chemical shift perturbation mapping and studies on the dynamics of proteins. This assignment is not a trivial task even if a 3D protein structure is known and needs almost as much effort as the assignment for structure prediction if performed manually. We present here a new algorithm based solely on 4D [(1)H,(15)N]-HSQC-NOESY-[(1)H,(15)N]-HSQC spectra which is able to assign a large percentage of chemical shifts (73-82 %) unambiguously, demonstrated with proteins up to a size of 250 residues. For the remaining residues, a small number of possible assignments is filtered out. This is done by comparing distances in the 3D structure to restraints obtained from the peak volumes in the 4D spectrum. Using dead-end elimination, assignments are removed in which at least one of the restraints is violated. Including additional information from chemical shift predictions, a complete unambiguous assignment was obtained for Ubiquitin and 95 % of the residues were correctly assigned in the 251 residue-long N-terminal domain of enzyme I. The program including source code is available at https://github.com/thomasexner/4Dassign .

Solution structure of the calmodulin-like C-terminal domain of Entamoeba α-actinin2

Cell motility is dependent on a dynamic meshwork of actin filaments that is remodelled continuously. A large number of associated proteins that are severs, cross-links, or caps the filament ends have been identified and the actin cross-linker α-actinin has been implied in several important cellular processes. In Entamoeba histolytica, the etiological agent of human amoebiasis, α-actinin is believed to be required for infection. To better understand the role of α-actinin in the infectious process we have determined the solution structure of the C-terminal calmodulin-like domain using NMR. The final structure ensemble of the apo form shows two lobes, that both resemble other pairs of calcium-binding EF-hand motifs, connected with a mobile linker.

Applications of high dimensionality experiments to biomolecular NMR

High dimensionality NMR experiments facilitate resonance assignment and precise determination of spectral parameters such as coupling constants. Sparse non-uniform sampling enables acquisition of experiments of high dimensionality with high resolution in acceptable time. In this review we present and compare some significant applications of NMR experiments of dimensionality higher than three in the field of biomolecular studies in solution.

HN-NCA heteronuclear TOCSY-NH experiment for (1)H(N) and (15)N sequential correlations in ((13)C, (15)N) labelled intrinsically disordered proteins

A simple triple resonance NMR experiment that leads to the correlation of the backbone amide resonances of each amino acid residue 'i' with that of residues 'i-1' and 'i+1' in ((13)C, (15)N) labelled intrinsically disordered proteins (IDPs) is presented. The experimental scheme, {HN-NCA heteronuclear TOCSY-NH}, exploits the favourable relaxation properties of IDPs and the presence of (1) J CαN and (2) J CαN couplings to transfer the (15)N x magnetisation from amino acid residue 'i' to adjacent residues via the application of a band-selective (15)N-(13)C(α) heteronuclear cross-polarisation sequence of ~100 ms duration. Employing non-uniform sampling in the indirect dimensions, the efficacy of the approach has been demonstrated by the acquisition of 3D HNN chemical shift correlation spectra of α-synuclein. The experimental performance of the RF pulse sequence has been compared with that of the conventional INEPT-based HN(CA)NH pulse scheme. As the availability of data from both the HCCNH and HNN experiments will make it possible to use the information extracted from one experiment to simplify the analysis of the data of the other and lead to a robust approach for unambiguous backbone and side-chain resonance assignments, a time-saving strategy for the simultaneous collection of HCCNH and HNN data is also described.

Non-Uniform Sampling and J-UNIO Automation for Efficient Protein NMR Structure Determination

High-resolution structure determination of small proteins in solution is one of the big assets of NMR spectroscopy in structural biology. Improvements in the efficiency of NMR structure determination by advances in NMR experiments and automation of data handling therefore attracts continued interest. Here, non-uniform sampling (NUS) of 3D heteronuclear-resolved [(1)H,(1)H]-NOESY data yielded two- to three-fold savings of instrument time for structure determinations of soluble proteins. With the 152-residue protein NP_372339.1 from Staphylococcus aureus and the 71-residue protein NP_346341.1 from Streptococcus pneumonia we show that high-quality structures can be obtained with NUS NMR data, which are equally well amenable to robust automated analysis as the corresponding uniformly sampled data.

Development of a method for reconstruction of crowded NMR spectra from undersampled time-domain data

NMR is a unique methodology for obtaining information about the conformational dynamics of proteins in heterogeneous biomolecular systems. In various NMR methods, such as transferred cross-saturation, relaxation dispersion, and paramagnetic relaxation enhancement experiments, fast determination of the signal intensity ratios in the NMR spectra with high accuracy is required for analyses of targets with low yields and stabilities. However, conventional methods for the reconstruction of spectra from undersampled time-domain data, such as linear prediction, spectroscopy with integration of frequency and time domain, and analysis of Fourier, and compressed sensing were not effective for the accurate determination of the signal intensity ratios of the crowded two-dimensional spectra of proteins. Here, we developed an NMR spectra reconstruction method, "conservation of experimental data in analysis of Fourier" (Co-ANAFOR), to reconstruct the crowded spectra from the undersampled time-domain data. The number of sampling points required for the transferred cross-saturation experiments between membrane proteins, photosystem I and cytochrome b 6 f, and their ligand, plastocyanin, with Co-ANAFOR was half of that needed for linear prediction, and the peak height reduction ratios of the spectra reconstructed from truncated time-domain data by Co-ANAFOR were more accurate than those reconstructed from non-uniformly sampled data by compressed sensing.

NMR Methods for the Study of Instrinsically Disordered Proteins Structure, Dynamics, and Interactions: General Overview and Practical Guidelines

Thanks to recent improvements in NMR instrumentation, pulse sequence design, and sample preparation, a panoply of new NMR tools has become available for atomic resolution characterization of intrinsically disordered proteins (IDPs) that are optimized for the particular chemical and spectroscopic properties of these molecules. A wide range of NMR observables can now be measured on increasingly complex IDPs that report on their structural and dynamic properties in isolation, as part of a larger complex, or even inside an entire living cell. Herein we present basic NMR concepts, as well as optimised tools available for the study of IDPs in solution. In particular, the following sections are discussed hereafter: a short introduction to NMR spectroscopy and instrumentation (Sect. 3.1), the effect of order and disorder on NMR observables (Sect. 3.2), particular challenges and bottlenecks for NMR studies of IDPs (Sect. 3.3), 2D HN and CON NMR experiments: the fingerprint of an IDP (Sect. 3.4), tools for overcoming major bottlenecks of IDP NMR studies (Sect. 3.5), 13C detected experiments (Sect. 3.6), from 2D to 3D: from simple snapshots to site-resolved characterization of IDPs (Sect. 3.7), sequential NMR assignment: 3D experiments (Sect. 3.8), high-dimensional NMR experiments (nD, with n>3) (Sect. 3.9) and conclusions and perspectives (Sect. 3.10).

Practical aspects of NMR signal assignment in larger and challenging proteins

NMR has matured into a technique routinely employed for studying proteins in near physiological conditions. However, applications to larger proteins are impeded by the complexity of the various correlation maps necessary to assign NMR signals. This article reviews the data analysis techniques traditionally employed for resonance assignment and describes alternative protocols necessary for overcoming challenges in large protein spectra. In particular, simultaneous analysis of multiple spectra may help overcome ambiguities or may reveal correlations in an indirect manner. Similarly, visualization of orthogonal planes in a multidimensional spectrum can provide alternative assignment procedures. We describe examples of such strategies for assignment of backbone, methyl, and nOe resonances. We describe experimental aspects of data acquisition for the related experiments and provide guidelines for preliminary studies. Focus is placed on large folded monomeric proteins and examples are provided for 37, 48, 53, and 81 kDa proteins.

Time-resolved multidimensional NMR with non-uniform sampling

Time-resolved experiments demand high resolution both in spectral dimensions and in time of the studied kinetic process. The latter requirement traditionally prohibits applications of the multidimensional experiments, which, although capable of providing invaluable information about structure and dynamics and almost unlimited spectral resolution, require too lengthy data collection. Our work shows that the problem has a solution in using modern methods of NMR data collection and signal processing. A continuous fast pulsing three-dimensional experiment is acquired using non-uniform sampling during full time of the studied reaction. High sensitivity and time-resolution of a few minutes is achieved by simultaneous processing of the full data set with the multi-dimensional decomposition. The method is verified and illustrated in realistic simulations and by measuring deuterium exchange rates of amide protons in ubiquitin. We applied the method for characterizing kinetics of in vitro phosphorylation of two tyrosine residues in an intrinsically disordered cytosolic domain of the B cell receptor protein CD79b. Signals of many residues including tyrosines in both phosphorylated and unmodified forms of CD79b are found in a heavily crowded region of 2D ¹H-¹⁵N correlation spectrum and the significantly enhanced spectral resolution provided by the 3D time-resolved approach was essential for the quantitative site-specific analysis.

In situ study of reaction kinetics using compressed sensing NMR

CS-NMR improves the temporal resolution of conventional multi-dimensional NMR for species identification and study of reaction kinetics.

Sampling scheme and compressed sensing applied to solid-state NMR spectroscopy

We describe the incorporation of non-uniform sampling (NUS) compressed sensing (CS) into oriented sample (OS) solid-state NMR for stationary aligned samples and magic angle spinning (MAS) Solid-state NMR for unoriented 'powder' samples. Both simulated and experimental results indicate that 25-33% of a full linearly sampled data set is required to reconstruct two- and three-dimensional solid-state NMR spectra with high fidelity. A modest increase in signal-to-noise ratio accompanies the reconstruction.

NMR methods for structural studies of large monomeric and multimeric proteins

NMR structural studies of large monomeric and multimeric proteins face distinct challenges. In large monomeric proteins, the common occurrence of frequency degeneracies between residues impedes unambiguous assignment of NMR signals. To overcome this barrier, nonuniform sampling (NUS) is used to measure spectra with optimal resolution within reasonable time, new correlation maps resolve previous impasses in assignment strategies, and novel selective methyl labeling schemes provide additional structural probes without cluttering NMR spectra. These advances push the limits of NMR studies of large monomeric proteins. Large multimeric and multidomain proteins are studied by NMR when individual components can also be studied by NMR and have known structures. The structural properties of large assemblies are obtained by identifying binding surfaces, by orienting domains, and employing limited distance constraints. Segmental labeling and the combination of NMR with other methods have helped popularize NMR studies of such systems.

4D Non-uniformly sampled C,C-NOESY experiment for sequential assignment of 13C, 15N-labeled RNAs

A 4D (13)C(aromatic),(13)C(ribose)-edited NOESY experiment is introduced to improve sequential assignment of non-coding RNA, often hampered by a limited dispersion of (1)H and (13)C chemical shifts. The (13)C-labeling of RNA is fully utilized by inclusion of two (13)C evolution periods. These dimensions provide enhanced dispersion of resonances in the 4D spectrum. High spectral resolution is obtained using random non-uniform sampling in three indirect dimensions. The autocorrelation peaks are efficiently suppressed using band-selective pulses. Since the dynamic range of observed resonances is significantly decreased, the reconstruction of the 4D spectrum is greatly simplified. The experiment can replace two conventionally sampled 3D NOESY spectra (either ribose-(13)C- or aromatic-(13)C-separated), and remove most ambiguities encountered during sequential walks. The assignment strategy based on a homonuclear and 4D C,C-edited NOESY experiments is proposed and verified on a 34-nt RNA showing typical structure elements.

Extrapolation and phase correction of non-uniformly broadened signals

The initial part of FID-signals cannot always be acquired experimentally. This is particularly true for signals characterised by strong inhomogeneous broadening, such as those in porous materials, e.g. cements, soils and rocks, those measured by portable NMR-apparatus, or EPR-signals. Here we report on a numerical method we designed to extrapolate those initial missing parts, i.e. to retrieve their amplitude and phase. Should the entire signal be available from an experiment, the algorithm can still be used as an automatic phase-corrector and a low-pass filter. The method is based on the use of cardinal series, applies to any oversampled signals and requires no prior knowledge of the system under study. We show that the method can also be used to restore entire one-dimensional MRI-data sets from those in which less than half of the k-space was sampled, thus not only potentially allowing to speed up data acquisition - when extended to two or three dimensions, but also to circumvent phase-distortions usually encountered when exploring the k-space near its origin.

Effects of NMR spectral resolution on protein structure calculation

Adequate digital resolution and signal sensitivity are two critical factors for protein structure determinations by solution NMR spectroscopy. The prime objective for obtaining high digital resolution is to resolve peak overlap, especially in NOESY spectra with thousands of signals where the signal analysis needs to be performed on a large scale. Achieving maximum digital resolution is usually limited by the practically available measurement time. We developed a method utilizing non-uniform sampling for balancing digital resolution and signal sensitivity, and performed a large-scale analysis of the effect of the digital resolution on the accuracy of the resulting protein structures. Structure calculations were performed as a function of digital resolution for about 400 proteins with molecular sizes ranging between 5 and 33 kDa. The structural accuracy was assessed by atomic coordinate RMSD values from the reference structures of the proteins. In addition, we monitored also the number of assigned NOESY cross peaks, the average signal sensitivity, and the chemical shift spectral overlap. We show that high resolution is equally important for proteins of every molecular size. The chemical shift spectral overlap depends strongly on the corresponding spectral digital resolution. Thus, knowing the extent of overlap can be a predictor of the resulting structural accuracy. Our results show that for every molecular size a minimal digital resolution, corresponding to the natural linewidth, needs to be achieved for obtaining the highest accuracy possible for the given protein size using state-of-the-art automated NOESY assignment and structure calculation methods.

Selective diagonal-free (13)C, (13)C-edited aliphatic-aromatic NOESY experiment with non-uniform sampling

A band-selective aromatic-aliphatic C,C-edited four-dimensional NOESY experiment is proposed here. Its key advantage is the absence of auto-correlation signals which makes it very attractive for joint use with non-uniform sampling. It is demonstrated here that the sensitivity of the experiment is not significantly affected by utilization of selective pulses (for either aromatic-13C or aliphatic-13C spins). The method was applied to the sample of E32Q mutant of human S100A1 protein, a homodimer of total molecular mass ~20 kDa. High-resolution 4D spectra were obtained from ~1.5 % of sampling points required conventionally. It is shown that superior resolution facilitates unambiguous assignment of observed aliphatic-aromatic cross-peaks. Additionally, the addition of aliphatic-13C dimension enables to resolve peaks with degenerated aliphatic (1)H chemical shifts. All observed cross-peaks were validated against previously determined 3D structure of E32Q mutant of S100A1 protein (PDB 2LHL). The increased reliability of structural constraints obtained from the proposed high-resolution 4D 13C(ali),13C(aro)-edited NOESY can be exploited in the automated protocols of structure determination of proteins.

Speeding up the measurement of one-bond scalar (1J) and residual dipolar couplings (1D) by using non-uniform sampling (NUS)

The accurate and precise measurement of one-bond scalar and residual dipolar coupling (RDC) constants is of prime importance to be able to use RDCs for structure determination. If coupling constants are to be extracted from the indirect dimension of HSQC spectra a significant saving of measurement time can be achieved by non-uniform sampling (NUS). Coupling constants can either be obtained with the same precision as in traditionally acquired spectra in a fraction of the measurement time or the precision can be significantly improved if the same amount of measurement time as for traditionally acquired spectra is invested. The application of NUS for the measurement of (1)J (scalar coupling constants) and (1)T (total couplings constants) from different kinds of ω(1)-coupled spectra (including also J-scaled ones) is examined in detail and the possible gains in time or resolution are discussed. When using the newly proposed compressed sensing (CS) algorithm for processing, the quality of the spectra is comparable to the traditionally sampled ones.

Suppression of sampling artefacts in high-resolution four-dimensional NMR spectra using signal separation algorithm

The development of non-uniform sampling (NUS) strategies permits to obtain high-dimensional spectra with increased resolution in significantly reduced experimental time. We extended a previously proposed signal separation algorithm (SSA) to process sparse four-dimensional NMR data. It is employed for two experiments carried out for a partially unstructured 114-residue construct of chicken Engrailed 2 protein, namely 4D HCCH-TOCSY and 4D C,N-edited NOESY. The SSA allowed us to obtain high-quality spectra using only as little as 0.16% of the available samples, with low sampling artefacts approaching the thermal noise level in most spectral regions. It is demonstrated that NUS 4D HCCH-TOCSY is dominated by sampling noise and requires efficient artefact suppression. On the other hand, 4D C,N-edited NOESY is a particularly attractive experiment for NUS, as the absence of diagonal peaks renders the problem of artefacts less critical. We also present a transverse-relaxation optimized sequence for HMQC that is especially designed for longer evolution periods in the indirectly detected proton dimension in high-dimensional pulse sequences. In conjunction with novel sampling strategies and efficient processing methods, this improvement enabled us to obtain unique structural information about aliphatic-amide contacts.

Deterministic schedules for robust and reproducible non-uniform sampling in multidimensional NMR

We show that a simple, general, and easily reproducible method for generating non-uniform sampling (NUS) schedules preserves the benefits of random sampling, including inherently reduced sampling artifacts, while removing the pitfalls associated with choosing an arbitrary seed. Sampling schedules are generated from a discrete cumulative distribution function (CDF) that closely fits the continuous CDF of the desired probability density function. We compare random and deterministic sampling using a Gaussian probability density function applied to 2D HSQC spectra. Data are processed using the previously published method of Spectroscopy by Integration of Frequency and Time domain data (SIFT). NUS spectra from deterministic sampling schedules were found to be at least as good as those from random schedules at the SIFT critical sampling density, and significantly better at half that sampling density. The method can be applied to any probability density function and generalized to greater than two dimensions.

Cardinal series to restore NMR-signals dominated by strong inhomogeneous broadening

We have devised two numerical methods of restoring incomplete band-limited NMR-signals to integrity by either interpolating or extrapolating them. Both methods are based on use of the finite cardinal series, whose filtering properties were discussed previously, to model signals. They require no prior knowledge about the system under study, but only that the available parts of the signal were oversampled enough. The methods were tested on two types of computer-simulated signal. It proved superior to the linear prediction methods and Lagrange interpolation when applied to signals measured in highly inhomogeneous magnetic fields. The extrapolation method was then applied to restore experimentally-measured refocused FID-signals of a porous medium. The missing parts of the signal of up to several times the size of its Nyquist period could be recovered by either method.

Exclusively NOESY-based automated NMR assignment and structure determination of proteins

A fully automated method is presented for determining NMR solution structures of proteins using exclusively NOESY spectra as input, obviating the need to measure any spectra only for obtaining resonance assignments but devoid of structural information. Applied to two small proteins, the approach yielded structures that coincided closely with conventionally determined structures.

Contemporary methods in structure determination of membrane proteins by solution NMR

Integral membrane proteins are vital to life, being responsible for information and material exchange between a cell and its environment. Although high-resolution structural information is needed to understand how these functions are achieved, membrane proteins remain an under-represented subset of the protein structure databank. Solution NMR is increasingly demonstrating its ability to help address this knowledge shortfall, with the development of a diverse array of techniques to counter the challenges presented by membrane proteins. Here we document the advances that are helping to define solution NMR as an effective tool for membrane protein structure determination. Developments introduced over the last decade in the production of isotope-labeled samples, reconstitution of these samples into the growing selection of NMR-compatible membrane-mimetic systems, and the approaches used for the acquisition and application of structural restraints from these complexes are reviewed.

A novel strategy for NMR resonance assignment and protein structure determination

The quality of protein structures determined by nuclear magnetic resonance (NMR) spectroscopy is contingent on the number and quality of experimentally-derived resonance assignments, distance and angular restraints. Two key features of protein NMR data have posed challenges for the routine and automated structure determination of small to medium sized proteins; (1) spectral resolution - especially of crowded nuclear Overhauser effect spectroscopy (NOESY) spectra, and (2) the reliance on a continuous network of weak scalar couplings as part of most common assignment protocols. In order to facilitate NMR structure determination, we developed a semi-automated strategy that utilizes non-uniform sampling (NUS) and multidimensional decomposition (MDD) for optimal data collection and processing of selected, high resolution multidimensional NMR experiments, combined it with an ABACUS protocol for sequential and side chain resonance assignments, and streamlined this procedure to execute structure and refinement calculations in CYANA and CNS, respectively. Two graphical user interfaces (GUIs) were developed to facilitate efficient analysis and compilation of the data and to guide automated structure determination. This integrated method was implemented and refined on over 30 high quality structures of proteins ranging from 5.5 to 16.5 kDa in size.

Automated projection spectroscopy and its applications

This chapter presents the NMR technique APSY (automated projection spectroscopy) and its applications for sequence-specific resonance assignments of proteins. The result of an APSY experiment is a list of chemical shift correlations for an N-dimensional NMR spectrum (N≥3). This list is obtained in a fully automated way by the dedicated algorithm GAPRO (geometric analysis of projections) from a geometric analysis of experimentally recorded, low-dimensional projections. Because the positions of corresponding peaks in multiple projections are correlated, thermal noise and other uncorrelated artifacts are efficiently suppressed. We describe the theoretical background of the APSY method and discuss technical aspects that guide its optimal use. Further, applications of APSY-NMR spectroscopy for fully automated sequence-specific backbone and side chain assignments of proteins are described. We discuss the choice of suitable experiments for this purpose and show several examples. APSY is of particular interest for the assignment of soluble unfolded proteins, which is a time-consuming task by conventional means. With this class of proteins, APSY-NMR experiments with up to seven dimensions have been recorded. Sequence-specific assignments of protein side chains in turn are obtained from a 5D TOCSY-APSY-NMR experiment.

NleG Type 3 effectors from enterohaemorrhagic Escherichia coli are U-Box E3 ubiquitin ligases

NleG homologues constitute the largest family of Type 3 effectors delivered by pathogenic E. coli, with fourteen members in the enterohaemorrhagic (EHEC) O157:H7 strain alone. Identified recently as part of the non-LEE-encoded (Nle) effector set, this family remained uncharacterised and shared no sequence homology to other proteins including those of known function. The C-terminal domain of NleG2-3 (residues 90 to 191) is the most conserved region in NleG proteins and was solved by NMR. Structural analysis of this structure revealed the presence of a RING finger/U-box motif. Functional assays demonstrated that NleG2-3 as well as NleG5-1, NleG6-2 and NleG9' family members exhibited a strong autoubiquitination activity in vitro; a characteristic usually expressed by eukaryotic ubiquitin E3 ligases. When screened for activity against a panel of 30 human E2 enzymes, the NleG2-3 and NleG5-1 homologues showed an identical profile with only UBE2E2, UBE2E3 and UBE2D2 enzymes supporting NleG activity. Fluorescence polarization analysis yielded a binding affinity constant of 56+/-2 microM for the UBE2D2/NleG5-1 interaction, a value comparable with previous studies on E2/E3 affinities. The UBE2D2 interaction interface on NleG2-3 defined by NMR chemical shift perturbation and mutagenesis was shown to be generally similar to that characterised for human RING finger ubiquitin ligases. The alanine substitutions of UBE2D2 residues Arg5 and Lys63, critical for activation of eukaryotic E3 ligases, also significantly decreased both NleG binding and autoubiquitination activity. These results demonstrate that bacteria-encoded NleG effectors are E3 ubiquitin ligases analogous to RING finger and U-box enzymes in eukaryotes.

Iterative algorithm of discrete Fourier transform for processing randomly sampled NMR data sets

Spectra obtained by application of multidimensional Fourier Transformation (MFT) to sparsely sampled nD NMR signals are usually corrupted due to missing data. In the present paper this phenomenon is investigated on simulations and experiments. An effective iterative algorithm for artifact suppression for sparse on-grid NMR data sets is discussed in detail. It includes automated peak recognition based on statistical methods. The results enable one to study NMR spectra of high dynamic range of peak intensities preserving benefits of random sampling, namely the superior resolution in indirectly measured dimensions. Experimental examples include 3D (15)N- and (13)C-edited NOESY-HSQC spectra of human ubiquitin.

Speeding up (13)C direct detection biomolecular NMR spectroscopy

After the exploitation of (1)H polarization as a starting source for (13)C direct detection experiments, pulse sequences are designed which exploit the accelerated (1)H longitudinal relaxation to expedite (13)C direct detection experiments. We show here that 2D experiments based on (13)C direct detection on a 0.5 mM water sample of ubiquitin can be recorded in a few minutes and 3D experiments in a few hours. We also show that fast methods like nonuniform sampling can be easily implemented. As overall experimental time has always been a counter indication for the use of (13)C direct detection experiments, this research opens new avenues for the application of (13)C NMR to biological molecules.

Coupled decomposition of four-dimensional NOESY spectra

Four-dimensional (4D) NOESY spectra provide unambiguous distance information at a resolution that cannot be achieved in fewer dimensions and thus increase the quality of biomolecular structure determination substantially. Since the degree of chemical shift degeneracy increases with protein size, the use of 4D NOESY spectra is particularly important for large proteins. The potential high resolution in 4D spectra cannot be achieved in a reasonable time with conventional acquisition routines that sample the Nyquist grid uniformly. It can, however, be obtained with nonuniform sampling of the data grid, but optimal processing of such data has not yet been established. Here we describe a processing method for a pair of sparsely sampled 4D NOESY spectra, a methyl-methyl and an amide-methyl NOESY, recorded on a perdeuterated protein with protonated isoleucine, leucine, and valine methyl groups. The coupled multidimensional decomposition (Co-MDD) of these two spectra together with a 2D template spectrum results in a substantial increase in sensitivity, evidenced by 50-100% additional cross peaks, when compared to alternative processing schemes. At the same time, Co-MDD allows the use of low sparse levels of 10-15% of the full data grid for NOESY spectra. For the 283-residue integral human membrane protein VDAC-1, which has a rotational correlation time of about 70 ns in detergent micelles, the two 4D Co-MDD NOESYs yielded a total of 366 NOEs, resulting in 139 unambiguous upper limit distance constraints for the structure calculation.

Automated NMR structure determination of stereo-array isotope labeled ubiquitin from minimal sets of spectra using the SAIL-FLYA system

Stereo-array isotope labeling (SAIL) has been combined with the fully automated NMR structure determination algorithm FLYA to determine the three-dimensional structure of the protein ubiquitin from different sets of input NMR spectra. SAIL provides a complete stereo- and regio-specific pattern of stable isotopes that results in sharper resonance lines and reduced signal overlap, without information loss. Here we show that as a result of the superior quality of the SAIL NMR spectra, reliable, fully automated analyses of the NMR spectra and structure calculations are possible using fewer input spectra than with conventional uniformly 13C/15N-labeled proteins. FLYA calculations with SAIL ubiquitin, using a single three-dimensional "through-bond" spectrum (and 2D HSQC spectra) in addition to the 13C-edited and 15N-edited NOESY spectra for conformational restraints, yielded structures with an accuracy of 0.83-1.15 A for the backbone RMSD to the conventionally determined solution structure of SAIL ubiquitin. NMR structures can thus be determined almost exclusively from the NOESY spectra that yield the conformational restraints, without the need to record many spectra only for determining intermediate, auxiliary data of the chemical shift assignments. The FLYA calculations for this report resulted in 252 ubiquitin structure bundles, obtained with different input data but identical structure calculation and refinement methods. These structures cover the entire range from highly accurate structures to seriously, but not trivially, wrong structures, and thus constitute a valuable database for the substantiation of structure validation methods.

Spectroscopy by integration of frequency and time domain information for fast acquisition of high-resolution dark spectra

A simple and effective method, SIFT (spectroscopy by integration of frequency and time domain information), is introduced for processing nonuniformly sampled multidimensional NMR data. Applying the computationally efficient Gerchberg-Papoulis (G-P) algorithm, used previously in picture processing and medical imaging, SIFT supplements data at nonuniform points in the time domain with the information carried by known "dark" points (i.e., empty regions) in the frequency domain. We demonstrate that this rapid integration not only removes the severe pseudonoise characteristic of the Fourier transforms of nonuniformly sampled data, but also provides a robust procedure for using frequency information to replace time measurements. The latter can be used to avoid unnecessary sampling in sampling-limited experiments, and the former can be used to take advantage of the ability of nonuniformly sampled data to minimize trade-offs between the signal-to-noise ratio and the resolution in sensitivity-limited experiments. Processing 2D and 3D data sets takes about 0.1 and 2 min, respectively, on a personal computer. With these several attractive features, SIFT offers a novel, model-independent, flexible, and user-friendly tool for efficient and accurate processing of multidimensional NMR data.

APSY-NMR with proteins: practical aspects and backbone assignment

Automated projection spectroscopy (APSY) is an NMR technique for the recording of discrete sets of projection spectra from higher-dimensional NMR experiments, with automatic identification of the multidimensional chemical shift correlations by the dedicated algorithm GAPRO. This paper presents technical details for optimizing the set-up and the analysis of APSY-NMR experiments with proteins. Since experience so far indicates that the sensitivity for signal detection may become the principal limiting factor for applications with larger proteins or more dilute samples, we performed an APSY-NMR experiment at the limit of sensitivity, and then investigated the effects of varying selected experimental parameters. To obtain the desired reference data, a 4D APSY-HNCOCA experiment with a 12-kDa protein was recorded in 13 min. Based on the analysis of this data set and on general considerations, expressions for the sensitivity of APSY-NMR experiments have been generated to guide the selection of the projection angles, the calculation of the sweep widths, and the choice of other acquisition and processing parameters. In addition, a new peak picking routine and a new validation tool for the final result of the GAPRO spectral analysis are introduced. In continuation of previous reports on the use of APSY-NMR for sequence-specific resonance assignment of proteins, we present the results of a systematic search for suitable combinations of a minimal number of four- and five-dimensional APSY-NMR experiments that can provide the input for algorithms that generate automated protein backbone assignments.

Maximum Entropy Spectral Reconstruction of Non-Uniformly Sampled Data

The time required to complete a multidimensional NMR experiment is directly proportional to the number of evolution times sampled in the indirect dimensions. A consequence when utilizing conventional methods of data acquisition and spectrum analysis is that resolution in the indirect dimensions is frequently sample-limited. The problem becomes more acute at very high magnetic fields, where increased chemical shift dispersion requires shorter time increments to avoid aliasing. It has long been recognized that a way to avoid this limitation is to utilize methods of spectrum analysis that do not require data to be sampled at uniform intervals, permitting the collection of data at long evolution times requisite for high resolution without requiring collection of data at all intervening multiples of the sampling interval. Several promising methods have evolved that are seemingly quite different, yet can be shown to yield similar results when applied to similar sampling strategies, emphasizing the importance of the choice of samples, regardless of the technique used to compute the spectrum. Maximum entropy (MaxEnt) reconstruction is a very general method for spectrum analysis of non-uniformly sampled data (NUS), and because it can be used with essentially arbitrary sampling strategies and makes no assumptions about the nature of the signal, it provides a convenient basis for exploring the influence of the choice of samples on spectral quality. In this article we use this versatility of MaxEnt reconstruction to compare different approaches to NUS in multidimensional NMR and suggest strategies for improving spectral quality by careful choice of sample times.

Automatic assignment of protein backbone resonances by direct spectrum inspection in targeted acquisition of NMR data

The necessity to acquire large multidimensional datasets, a basis for assignment of NMR resonances, results in long data acquisition times during which substantial degradation of a protein sample might occur. Here we propose a method applicable for such a protein for automatic assignment of backbone resonances by direct inspection of multidimensional NMR spectra. In order to establish an optimal balance between completeness of resonance assignment and losses of cross-peaks due to dynamic processes/degradation of protein, assignment of backbone resonances is set as a stirring criterion for dynamically controlled targeted nonlinear NMR data acquisition. The result is demonstrated with the 12 kDa (13)C,(15) N-labeled apo-form of heme chaperone protein CcmE, where hydrolytic cleavage of 29 C-terminal amino acids is detected. For this protein, 90 and 98% of manually assignable resonances are automatically assigned within 10 and 40 h of nonlinear sampling of five 3D NMR spectra, respectively, instead of 600 h needed to complete the full time domain grid. In addition, resonances stemming from degradation products are identified. This study indicates that automatic resonance assignment might serve as a guiding criterion for optimal run-time allocation of NMR resources in applications to proteins prone to degradation.

Hyperdimensional NMR spectroscopy with nonlinear sampling

An approach is described for joint interleaved recording, real-time processing, and analysis of NMR data sets. The method employs multidimensional decomposition to find common information in a set of conventional triple-resonance spectra recorded in the nonlinear sampling mode, and builds a model of hyperdimensional (HD) spectrum. While preserving sensitivity per unit of measurement time and allowing for maximal spectral resolution, the approach reduces data collection time on average by 2 orders of magnitude compared to the conventional method. The 7-10 dimensional HD spectrum, which is represented as a set of deconvoluted 1D vectors, is easy to handle and amenable for automated analysis. The method is exemplified by automated assignment for two protein systems of low and high spectral complexity: ubiquitin (globular, 8 kDa) and zetacyt (naturally disordered, 13 kDa). The collection and backbone assignment of the data sets are achieved in real time after approximately 1 and 10 h, respectively. The approach removes the most critical time bottlenecks in data acquisition and analysis. Thus, it can significantly increase the value of NMR spectroscopy in structural biology, for example, in high-throughput structural genomics applications.

Lineshapes and artifacts in Multidimensional Fourier Transform of arbitrary sampled NMR data sets

The comprehensive description of Multidimensional Fourier Transform applied to arbitrary sampled NMR data is presented. Lineshapes and signal-to-artifact ratio are discussed in detail with regard to time domain sampling scheme and applied data weighting. It is demonstrated that transformation method with simple summation instead of numerical integration is most useful for significantly undersampled experiments. Additionally, the optimized random sampling schedule which enables significant improvement of obtained spectra is proposed. The new procedure of cleaning spectra is presented, it is based on predictability of artifacts pattern when sampling scheme and amplitude of intense signals are known. The results enable observation of high dynamic range spectra as for example heteronuclear edited NOESY. We show the application of new approach to the 3D (15)N-edited NOESY-HSQC spectrum acquired for (13)C, (15)N labeled ubiquitin sample with random time domain sampling.