Framework for and evaluation of bursts in random sampling of multidimensional NMR experiments

nus-tool: A unified program for generating and analyzing sample schedules for nonuniformly sampled NMR experiments

Increases in digital resolution achieved by high-field NMR require increases in spectral width. Additionally, the ability to resolve two overlapping peaks requires a sufficiently long acquisition time. These constraints combine, so that achieving high resolution spectra on high-field magnets requires long experiment times when employing uniform sampling and Fourier Transform processing. These limitations may be addressed by using nonuniform sampling (NUS), but the complexity of the parameter space across the variety of available NUS schemes greatly hinders the establishment of optimal approaches and best practices. We address these challenges with nus-tool, which is a software package for generating and analyzing NUS schedules. The nus-tool software internally implements random sampling and exponentially biased sampling. Through pre-configured plug-ins, it also provides access to quantile sampling and Poisson gap sampling. The software computes the relative sensitivity, mean evolution time, point spread function, and peak-to-sidelobe ratio; all of which can be determined for a candidate sample schedule prior to running an experiment to verify expected sensitivity, resolution, and artifact suppression. The nus-tool package is freely available on the NMRbox platform through an interactive GUI and via the command line, which is especially useful for scripted workflows that investigate the effectiveness of various NUS schemes.

Revisiting aliasing noise to build more robust sparsity in nonuniform sampling 2D-NMR

A continuing priority is to better understand and resolve the barriers to using nonuniform sampling (NUS) in challenging small molecule 2D NMR with subsampling of the Nyquist grid (a.k.a. coverage) below 50%. Possible causes for artifacts, often termed sampling noise, in 1D-NUS of 2D-NMR are revisited here, where weak aliasing artifacts are a growing concern as NUS becomes sparser. As NUS schedules become sparser, repeat sequences are shown to occur in the dense sampling regions early in the sampling schedule, causing aliasing artifacts in resulting spectra. An intuitive screening approach that detects patterns in sampling schedules based on a convolutional filter was implemented. Sampling schedules that have low proportions of repeat sequences show significantly reduced artifacts. Another route to remediate early repeat sequences is a short period of uniform sampling at the beginning of the schedule, which also leads to a significant suppression of unwanted sampling noise. Combining the repeat sequence filter with a survey of HSQC and LR-HSQMBC experiments, it is shown that very short initial uniform regions of about 2%-4% of the sampling space can ameliorate repeat sequences in sparser NUS and lead to robust spectral reconstructions by iterative soft thresholding (IST), even when the point spread function is unchanged. Using the principles developed here, a suite of 'one-click' schedules was developed for broader use.

Clustered sparsity and Poisson-gap sampling

Non-uniform sampling (NUS) is a popular way of reducing the amount of time taken by multidimensional NMR experiments. Among the various non-uniform sampling schemes that exist, the Poisson-gap (PG) schedules are particularly popular, especially when combined with compressed-sensing (CS) reconstruction of missing data points. However, the use of PG is based mainly on practical experience and has not, as yet, been explained in terms of CS theory. Moreover, an apparent contradiction exists between the reported effectiveness of PG and CS theory, which states that a "flat" pseudo-random generator is the best way to generate sampling schedules in order to reconstruct sparse spectra. In this paper we explain how, and in what situations, PG reveals its superior features in NMR spectroscopy. We support our theoretical considerations with simulations and analyses of experimental data from the Biological Magnetic Resonance Bank (BMRB). Our analyses reveal a previously unnoticed feature of many NMR spectra that explains the success of "blue-noise" schedules, such as PG. We call this feature "clustered sparsity". This refers to the fact that the peaks in NMR spectra are not just sparse but often form clusters in the indirect dimension, and PG is particularly suited to deal with such situations. Additionally, we discuss why denser sampling in the initial and final parts of the clustered signal may be useful.

Non-uniform sampling in pulse dipolar spectroscopy by EPR: the redistribution of noise and the optimization of data acquisition

Pulse dipolar spectroscopy (PDS) in Electron Paramagnetic Resonance (EPR) is the method of choice for determining the distance distribution function for mono-, bi- or multi- spin-labeled macromolecules and nanostructures. PDS acquisition schemes conventionally use uniform sampling of the dipolar trace, but non-uniform sampling (NUS) schemes can decrease the total measurement time or increase the accuracy of the resulting distance distributions. NUS requires optimization of the data acquisition scheme, as well as changes in data processing algorithms to accommodate the non-uniformly sampled data. We investigate in silico the applicability of the NUS approach in PDS, considering its effect on random, truncation and sampling noise in the experimental data. Each type of noise in the time-domain data propagates differently and non-uniformly into the distance spectrum as errors in the distance distribution. NUS schemes seem to be a valid approach for increasing sensitivity and/or throughput in PDS by decreasing and redistributing noise in the distance spectrum so that it has less impact on the distance spectrum.

Developing nonuniform sampling strategies to improve sensitivity and resolution in 1,1-ADEQUATE experiments

Nonuniform sampling (NUS) strategies are developed for acquiring highly resolved 1,1-ADEQUATE spectra, in both conventional and homodecoupled (HD) variants with improved sensitivity. Specifically, the quantile-directed and Poisson gap methods were critically compared for distributing the samples nonuniformly, and the quantile schedules were further optimized for weighting. Both maximum entropy and iterative soft thresholding spectral estimation algorithms were evaluated. All NUS approaches were robust when the degree of data reduction is moderate, on the order of a 50% reduction of sampling points. Further sampling reduction by NUS is facilitated by using weighted schedules designed by the quantile method, which also suppresses sampling noise well. Seed independence and the ability to specify the sample weighting in quantile scheduling are important in optimizing NUS for 1,1-ADEQUATE data acquisition. Using NUS yields an improvement in sensitivity, while also making longer evolution times accessible that would be difficult or impractical to attain by uniform sampling. Theoretical predictions for the sensitivity enhancements in these experiments are in the range of 5-20%; NUS is shown to disambiguate weak signals, reveal some ⁿ J_CC correlations obscured by noise, and improve signal strength relative to uniform sampling in the same experimental time. This work presents sample schedule development for applying NUS to challenging experiments. The schedules developed here are made available for general use and should facilitate the broader utilization of ADEQUATE experiments (including 1,1-, 1,n-, and HD- variants) for challenging structure elucidation problems.

Enhancing Compression Level for More Efficient Compressed Sensing and Other Lessons from NMR Spectroscopy

Modern nuclear magnetic resonance spectroscopy (NMR) is based on two- and higher-dimensional experiments that allow the solving of molecular structures, i.e., determine the relative positions of single atoms very precisely. However, rich chemical information comes at the price of long data acquisition times (up to several days). This problem can be alleviated by compressed sensing (CS)—a method that revolutionized many fields of technology. It is known that CS performs the most efficiently when measured objects feature a high level of compressibility, which in the case of NMR signal means that its frequency domain representation (spectrum) has a low number of significant points. However, many NMR spectroscopists are not aware of the fact that various well-known signal acquisition procedures enhance compressibility and thus should be used prior to CS reconstruction. In this study, we discuss such procedures and show to what extent they are complementary to CS approaches. We believe that the survey will be useful not only for NMR spectroscopists but also to inspire the broader signal processing community.

The influence of the probability density function on spectral quality in nonuniformly sampled multidimensional NMR

The goal of nonuniform sampling (NUS) is to select a subset of free induction decays (FIDs) from the conventional, uniform grid in a manner that sufficiently samples short evolution times needed for improved sensitivity and long evolution times needed for enhanced resolution. In addition to specifying the number of FIDs to be collected from a uniform grid, NUS schemes also specify the distribution of the selected FIDs, which directly impacts sampling-induced artifacts. Sampling schemes typically address these heuristic guidelines by utilizing a probability density function (PDF) to bias the distribution of sampled evolution times. Given this common approach, schemes differentiate themselves by how the evolution times are distributed within the envelope of the PDF. Here, we employ maximum entropy reconstruction and utilize in situ receiver operating characteristic (IROC) to conduct a critical comparison of the sensitivity and resolution that can be achieved by three types of biased sampling schemes: exponential (PDF is exponentially decaying), Poisson-gap (PDF derived from a sine function), and quantile-directed (PDF defined by simple polynomial decay). This methodology reveals practical insights and trends regarding how the sampling schemes and bias can provide the highest sensitivity and resolution for two nonuniformly sampled dimensions in a three-dimensional biomolecular NMR experiment. The IROC analysis circumvents the limitations of common metrics when used with nonlinear spectral estimation (a characteristic of all methods used with NUS) by quantifying the spectral quality via synthetic signals that are added to the empirical dataset. Recovery of these synthetic signals provides a proxy for the quality of the empirical portion of the spectrum. The central finding is that differences in spectral quality are primarily driven by the strength of bias in the PDF. In addition, a sampling coverage threshold is observed that appears to be connected to the dependence of each NUS method on its random seed. The differences between sampling schemes and biases are most relevant below 20% coverage where seed-dependence is high, whereas at higher coverages, the performance metrics for all of the sampling schemes begin to converge and approach a seed-independent regime. The results presented here show that aggressive sampling at low coverage can produce high-quality spectra by employing a sampling scheme that adheres to a decaying PDF with a bias to a broad range of short evolution times and includes relatively few FIDs at long evolution times.

Accelerated acquisition in pure-shift spectra based on prior knowledge from 1H NMR

Pure shift NMR with maximum performance by non-uniform sampling with prior knowledge.