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.

Levels and Characteristics of mRNAs in Spores of Firmicute Species.

Spores of firmicute species contain 100s of mRNAs, whose major function in Bacillus subtilis is to provide ribonucleotides for new RNA synthesis when spores germinate. To determine if this is a general phenomenon, RNA was isolated from spores of multiple firmicute species and relative mRNA levels determined by transcriptome sequencing (RNA-seq). Determination of RNA levels in single spores allowed calculation of RNA nucleotides/spore, and assuming mRNA is 3% of spore RNA indicated that only ∼6% of spore mRNAs were present at >1/spore. Bacillus subtilis, Bacillus atrophaeus, and Clostridioides difficile spores had 49, 42, and 51 mRNAs at >1/spore, and numbers of mRNAs at ≥1/spore were ∼10 to 50% higher in Geobacillus stearothermophilus and Bacillus thuringiensis Al Hakam spores and ∼4-fold higher in Bacillus megaterium spores. In all species, some to many abundant spore mRNAs (i) were transcribed by RNA polymerase with forespore-specific σ factors, (ii) encoded proteins that were homologs of those encoded by abundant B. subtilis spore mRNAs and are proteins in dormant spores, and (iii) were likely transcribed in the mother cell compartment of the sporulating cell. Analysis of the coverage of RNA-seq reads on mRNAs from all species suggested that abundant spore mRNAs were fragmented, as was confirmed by reverse transcriptase quantitative PCR (RT-qPCR) analysis of abundant B. subtilis and C. difficile spore mRNAs. These data add to evidence indicating that the function of at least the great majority of mRNAs in all firmicute spores is to be degraded to generate ribonucleotides for new RNA synthesis when spores germinate. IMPORTANCE Only ∼6% of mRNAs in spores of six firmicute species are at ≥1 molecule/spore, many abundant spore mRNAs encode proteins similar to B. subtilis spore proteins, and some abundant B. subtilis and C. difficile spore mRNAs were fragmented. Most of the abundant B. subtilis and other Bacillales spore mRNAs are transcribed under the control of the forespore-specific RNA polymerase σ factors, F or G, and these results may stimulate transcription analyses in developing spores of species other than B. subtilis. These findings, plus the absence of key nucleotide biosynthetic enzymes in spores, suggest that firmicute spores' abundant mRNAs are not translated when spores germinate but instead are degraded to generate ribonucleotides for new RNA synthesis by the germinated spore.

NUScon: a community-driven platform for quantitative evaluation of nonuniform sampling in NMR.

Although the concepts of nonuniform sampling (NUS​​​​​​​) and non-Fourier spectral reconstruction in multidimensional NMR began to emerge 4 decades ago , it is only relatively recently that NUS has become more commonplace. Advantages of NUS include the ability to tailor experiments to reduce data collection time and to improve spectral quality, whether through detection of closely spaced peaks (i.e., "resolution") or peaks of weak intensity (i.e., "sensitivity"). Wider adoption of these methods is the result of improvements in computational performance, a growing abundance and flexibility of software, support from NMR spectrometer vendors, and the increased data sampling demands imposed by higher magnetic fields. However, the identification of best practices still remains a significant and unmet challenge. Unlike the discrete Fourier transform, non-Fourier methods used to reconstruct spectra from NUS data are nonlinear, depend on the complexity and nature of the signals, and lack quantitative or formal theory describing their performance. Seemingly subtle algorithmic differences may lead to significant variabilities in spectral qualities and artifacts. A community-based critical assessment of NUS challenge problems has been initiated, called the "Nonuniform Sampling Contest" (NUScon), with the objective of determining best practices for processing and analyzing NUS experiments. We address this objective by constructing challenges from NMR experiments that we inject with synthetic signals, and we process these challenges using workflows submitted by the community. In the initial rounds of NUScon our aim is to establish objective criteria for evaluating the quality of spectral reconstructions. We present here a software package for performing the quantitative analyses, and we present the results from the first two rounds of NUScon. We discuss the challenges that remain and present a roadmap for continued community-driven development with the ultimate aim of providing best practices in this rapidly evolving field. The NUScon software package and all data from evaluating the challenge problems are hosted on the NMRbox platform.

Merging NMR Data and Computation Facilitates Data-Centered Research.

The Biological Magnetic Resonance Data Bank (BMRB) has served the NMR structural biology community for 40 years, and has been instrumental in the development of many widely-used tools. It fosters the reuse of data resources in structural biology by embodying the FAIR data principles (Findable, Accessible, Inter-operable, and Re-usable). NMRbox is less than a decade old, but complements BMRB by providing NMR software and high-performance computing resources, facilitating the reuse of software resources. BMRB and NMRbox both facilitate reproducible research. NMRbox also fosters the development and deployment of complex meta-software. Combining BMRB and NMRbox helps speed and simplify workflows that utilize BMRB, and enables facile federation of BMRB with other data repositories. Utilization of BMRB and NMRbox in tandem will enable additional advances, such as machine learning, that are poised to become increasingly powerful.

Analysis of 5'-NAD capping of mRNAs in dormant spores of Bacillus subtilis.

Spores of Gram-positive bacteria contain 10s-1000s of different mRNAs. However, Bacillus subtilis spores contain only ∼ 50 mRNAs at > 1 molecule/spore, almost all transcribed only in the developing spore and encoding spore proteins. However, some spore mRNAs could be stabilized to ensure they are intact in dormant spores, perhaps to direct synthesis of proteins essential for spores' conversion to a growing cell in germinated spore outgrowth. Recent work shows that some growing B. subtilis cell mRNAs contain a 5'-NAD cap. Since this cap may stabilize mRNA in vivo, its presence on spore mRNAs would suggest that maintaining some intact spore mRNAs is important, perhaps because they have a translational role in outgrowth. However, significant levels of only a few abundant spore mRNAs had a 5'-NAD cap, and these were not the most stable spore mRNAs and had likely been fragmented. Even higher levels of 5'-NAD-capping were found on a few low abundance spore mRNAs, but these mRNAs were present in only small percentages of spores, and had again been fragmented. The new data are thus consistent with spore mRNAs serving only as a reservoir of ribonucleotides in outgrowth.

DNA Damage Kills Bacterial Spores and Cells Exposed to 222-Nanometer UV Radiation.

This study examined the microbicidal activity of 222-nm UV radiation (UV222), which is potentially a safer alternative to the 254-nm UV radiation (UV254) that is often used for surface decontamination. Spores and/or growing and stationary-phase cells of Bacillus cereus, Bacillus subtilis, Bacillus thuringiensis, Staphylococcus aureus, and Clostridioides difficile and a herpesvirus were all killed or inactivated by UV222 and at lower fluences than with UV254B. subtilis spores and cells lacking the major DNA repair protein RecA were more sensitive to UV222, as were spores lacking their DNA-protective proteins, the α/ß-type small, acid-soluble spore proteins. The spore cores' large amount of Ca2+-dipicolinic acid (∼25% of the core dry weight) also protected B. subtilis and C. difficile spores against UV222, while spores' proteinaceous coat may have given some slight protection against UV222 Survivors among B. subtilis spores treated with UV222 acquired a large number of mutations, and this radiation generated known mutagenic photoproducts in spore and cell DNA, primarily cyclobutane-type pyrimidine dimers in growing cells and an α-thyminyl-thymine adduct termed the spore photoproduct (SP) in spores. Notably, the loss of a key SP repair protein markedly decreased spore UV222 resistance. UV222-treated B. subtilis spores germinated relatively normally, and the generation of colonies from these germinated spores was not salt sensitive. The latter two findings suggest that UV222 does not kill spores by general protein damage, and thus, the new results are consistent with the notion that DNA damage is responsible for the killing of spores and cells by UV222IMPORTANCE Spores of a variety of bacteria are resistant to common decontamination agents, and many of them are major causes of food spoilage and some serious human diseases, including anthrax caused by spores of Bacillus anthracis Consequently, there is an ongoing need for efficient methods for spore eradication, in particular methods that have minimal deleterious effects on people or the environment. UV radiation at 254 nm (UV254) is sporicidal and commonly used for surface decontamination but can cause deleterious effects in humans. Recent work, however, suggests that 222-nm UV (UV222) may be less harmful to people than UV254 yet may still kill bacteria and at lower fluences than UV254 The present work has identified the damage by UV222 that leads to the killing of growing cells and spores of some bacteria, many of which are human pathogens, and UV222 also inactivates a herpesvirus.

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

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.

Nonuniform sampling by quantiles.

A flexible strategy for choosing samples nonuniformly from a Nyquist grid using the concept of statistical quantiles is presented for broad classes of NMR experimentation. Quantile-directed scheduling is intuitive and flexible for any weighting function, promotes reproducibility and seed independence, and is generalizable to multiple dimensions. In brief, weighting functions are divided into regions of equal probability, which define the samples to be acquired. Quantile scheduling therefore achieves close adherence to a probability distribution function, thereby minimizing gaps for any given degree of subsampling of the Nyquist grid. A characteristic of quantile scheduling is that one-dimensional, weighted NUS schedules are deterministic, however higher dimensional schedules are similar within a user-specified jittering parameter. To develop unweighted sampling, we investigated the minimum jitter needed to disrupt subharmonic tracts, and show that this criterion can be met in many cases by jittering within 25-50% of the subharmonic gap. For nD-NUS, three supplemental components to choosing samples by quantiles are proposed in this work: (i) forcing the corner samples to ensure sampling to specified maximum values in indirect evolution times, (ii) providing an option to triangular backfill sampling schedules to promote dense/uniform tracts at the beginning of signal evolution periods, and (iii) providing an option to force the edges of nD-NUS schedules to be identical to the 1D quantiles. Quantile-directed scheduling meets the diverse needs of current NUS experimentation, but can also be used for future NUS implementations such as off-grid NUS and more. A computer program implementing these principles (a.k.a. QSched) in 1D- and 2D-NUS is available under the general public license.