Quantum Chemistry Calculations for Metabolomics.

A primary goal of metabolomics studies is to fully characterize the small-molecule composition of complex biological and environmental samples. However, despite advances in analytical technologies over the past two decades, the majority of small molecules in complex samples are not readily identifiable due to the immense structural and chemical diversity present within the metabolome. Current gold-standard identification methods rely on reference libraries built using authentic chemical materials ("standards"), which are not available for most molecules. Computational quantum chemistry methods, which can be used to calculate chemical properties that are then measured by analytical platforms, offer an alternative route for building reference libraries, i.e., in silico libraries for "standards-free" identification. In this review, we cover the major roadblocks currently facing metabolomics and discuss applications where quantum chemistry calculations offer a solution. Several successful examples for nuclear magnetic resonance spectroscopy, ion mobility spectrometry, infrared spectroscopy, and mass spectrometry methods are reviewed. Finally, we consider current best practices, sources of error, and provide an outlook for quantum chemistry calculations in metabolomics studies. We expect this review will inspire researchers in the field of small-molecule identification to accelerate adoption of in silico methods for generation of reference libraries and to add quantum chemistry calculations as another tool at their disposal to characterize complex samples.

DEIMoS: An Open-Source Tool for Processing High-Dimensional Mass Spectrometry Data.

We present DEIMoS: Data Extraction for Integrated Multidimensional Spectrometry, a Python application programming interface (API) and command-line tool for high-dimensional mass spectrometry data analysis workflows that offers ease of development and access to efficient algorithmic implementations. Functionality includes feature detection, feature alignment, collision cross section (CCS) calibration, isotope detection, and MS/MS spectral deconvolution, with the output comprising detected features aligned across study samples and characterized by mass, CCS, tandem mass spectra, and isotopic signature. Notably, DEIMoS operates on N-dimensional data, largely agnostic to acquisition instrumentation; algorithm implementations simultaneously utilize all dimensions to (i) offer greater separation between features, thus improving detection sensitivity, (ii) increase alignment/feature matching confidence among data sets, and (iii) mitigate convolution artifacts in tandem mass spectra. We demonstrate DEIMoS with LC-IMS-MS/MS metabolomics data to illustrate the advantages of a multidimensional approach in each data processing step.

Exploring the Impacts of Conformer Selection Methods on Ion Mobility Collision Cross Section Predictions.

The prediction of structure dependent molecular properties, such as collision cross sections as measured using ion mobility spectrometry, are crucially dependent on the selection of the correct population of molecular conformers. Here, we report an in-depth evaluation of multiple conformation selection techniques, including simple averaging, Boltzmann weighting, lowest energy selection, low energy threshold reductions, and similarity reduction. Generating 50 000 conformers each for 18 molecules, we used the In Silico Chemical Library Engine (ISiCLE) to calculate the collision cross sections for the entire data set. First, we employed Monte Carlo simulations to understand the variability between conformer structures as generated using simulated annealing. Then we employed Monte Carlo simulations to the aforementioned conformer selection techniques applied on the simulated molecular property: the ion mobility collision cross section. Based on our analyses, we found Boltzmann weighting to be a good trade-off between precision and theoretical accuracy. Combining multiple techniques revealed that energy thresholds and root-mean-squared deviation-based similarity reductions can save considerable computational expense while maintaining property prediction accuracy. Molecular dynamic conformer generation tools like AMBER can continue to generate new lowest energy conformers even after tens of thousands of generations, decreasing precision between runs. This reduced precision can be ameliorated and theoretical accuracy increased by running density functional theory geometry optimization on carefully selected conformers.

A Practical Guide to Metabolomics Software Development.

A growing number of software tools have been developed for metabolomics data processing and analysis. Many new tools are contributed by metabolomics practitioners who have limited prior experience with software development, and the tools are subsequently implemented by users with expertise that ranges from basic point-and-click data analysis to advanced coding. This Perspective is intended to introduce metabolomics software users and developers to important considerations that determine the overall impact of a publicly available tool within the scientific community. The recommendations reflect the collective experience of an NIH-sponsored Metabolomics Consortium working group that was formed with the goal of researching guidelines and best practices for metabolomics tool development. The recommendations are aimed at metabolomics researchers with little formal background in programming and are organized into three stages: (i) preparation, (ii) tool development, and (iii) distribution and maintenance.

Ligand- and Structure-Based Analysis of Deep Learning-Generated Potential α2a Adrenoceptor Agonists.

The α2a adrenoceptor is a medically relevant subtype of the G protein-coupled receptor family. Unfortunately, high-throughput techniques aimed at producing novel drug leads for this receptor have been largely unsuccessful because of the complex pharmacology of adrenergic receptors. As such, cutting-edge in silico ligand- and structure-based assessment and de novo deep learning methods are well positioned to provide new insights into protein-ligand interactions and potential active compounds. In this work, we (i) collect a dataset of α2a adrenoceptor agonists and provide it as a resource for the drug design community; (ii) use the dataset as a basis to generate candidate-active structures via deep learning; and (iii) apply computational ligand- and structure-based analysis techniques to gain new insights into α2a adrenoceptor agonists and assess the quality of the computer-generated compounds. We further describe how such assessment techniques can be applied to putative chemical probes with a case study involving proposed medetomidine-based probes.

Mass Spectrometry Adduct Calculator.

We describe the Mass Spectrometry Adduct Calculator (MSAC), an automated Python tool to calculate the adduct ion masses of a parent molecule. Here, adduct refers to a version of a parent molecule [M] that is charged due to addition or loss of atoms and electrons resulting in a charged ion, for example, [M + H]+. MSAC includes a database of 147 potential adducts and adduct/neutral loss combinations and their mass-to-charge ratios (m/z) as extracted from the NIST/EPA/NIH Mass Spectral Library (NIST17), Global Natural Products Social Molecular Networking Public Spectral Libraries (GNPS), and MassBank of North America (MoNA). The calculator relies on user-selected subsets of the combined database to calculate expected m/z for adducts of molecules supplied as formulas. This tool is intended to help researchers create identification libraries to collect evidence for the presence of molecules in mass spectrometry data. While the included adduct database focuses on adducts typically detected during liquid chromatography-mass spectrometry analyses, users may supply their own lists of adducts and charge states for calculating expected m/z. We also analyzed statistics on adducts from spectra contained in the three selected mass spectral libraries. MSAC is freely available at https://github.com/pnnl/MSAC.

Application and assessment of deep learning for the generation of potential NMDA receptor antagonists.

Uncompetitive antagonists of the N-methyl d-aspartate receptor (NMDAR) have demonstrated therapeutic benefit in the treatment of neurological diseases such as Parkinson's and Alzheimer's, but some also cause dissociative effects that have led to the synthesis of illicit drugs. The ability to generate NMDAR antagonists in silico is therefore desirable for both new medication development and preempting and identifying new designer drugs. Recently, generative deep learning models have been applied to de novo drug design as a means to expand the amount of chemical space that can be explored for potential drug-like compounds. In this study, we assess the application of a generative model to the NMDAR to achieve two primary objectives: (i) the creation and release of a comprehensive library of experimentally validated NMDAR phencyclidine (PCP) site antagonists to assist the drug discovery community and (ii) an analysis of both the advantages conferred by applying such generative artificial intelligence models to drug design and the current limitations of the approach. We apply, and provide source code for, a variety of ligand- and structure-based assessment techniques used in standard drug discovery analyses to the deep learning-generated compounds. We present twelve candidate antagonists that are not available in existing chemical databases to provide an example of what this type of workflow can achieve, though synthesis and experimental validation of these compounds are still required.

Chespa: Streamlining Expansive Chemical Space Evaluation of Molecular Sets.

Thousands of chemical properties can be calculated for small molecules, which can be used to place the molecules within the context of a broader "chemical space." These definitions vary based on compounds of interest and the goals for the given chemical space definition. Here, we introduce a customizable Python module, chespa, built to easily assess different chemical space definitions through clustering of compounds in these spaces and visualizing trends of these clusters. To demonstrate this, chespa currently streamlines prediction of various molecular descriptors (predicted chemical properties, molecular substructures, AI-based chemical space, and chemical class ontology) in order to test six different chemical space definitions. Furthermore, we investigated how these varying definitions trend with mass spectrometry (MS)-based observability, that is, the ability of a molecule to be observed with MS (e.g., as a function of the molecule ionizability), using an example data set from the U.S. EPA's nontargeted analysis collaborative trial, where blinded samples had been analyzed previously, providing 1398 data points. Improved understanding of observability would offer many advantages in small-molecule identification, such as (i) a priori selection of experimental conditions based on suspected sample composition, (ii) the ability to reduce the number of candidate structures during compound identification by removing those less likely to ionize, and, in turn, (iii) a reduced false discovery rate and increased confidence in identifications. Factors controlling observability are not fully understood, making prediction of this property nontrivial and a prime candidate for chemical space analysis. Chespa is available at github.com/pnnl/chespa.

ISiCLE: A Quantum Chemistry Pipeline for Establishing in Silico Collision Cross Section Libraries.

High-throughput, comprehensive, and confident identifications of metabolites and other chemicals in biological and environmental samples will revolutionize our understanding of the role these chemically diverse molecules play in biological systems. Despite recent technological advances, metabolomics studies still result in the detection of a disproportionate number of features that cannot be confidently assigned to a chemical structure. This inadequacy is driven by the single most significant limitation in metabolomics, the reliance on reference libraries constructed by analysis of authentic reference materials with limited commercial availability. To this end, we have developed the in silico chemical library engine (ISiCLE), a high-performance computing-friendly cheminformatics workflow for generating libraries of chemical properties. In the instantiation described here, we predict probable three-dimensional molecular conformers (i.e., conformational isomers) using chemical identifiers as input, from which collision cross sections (CCS) are derived. The approach employs first-principles simulation, distinguished by the use of molecular dynamics, quantum chemistry, and ion mobility calculations, to generate structures and chemical property libraries, all without training data. Importantly, optimization of ISiCLE included a refactoring of the popular MOBCAL code for trajectory-based mobility calculations, improving its computational efficiency by over 2 orders of magnitude. Calculated CCS values were validated against 1983 experimentally measured CCS values and compared to previously reported CCS calculation approaches. Average calculated CCS error for the validation set is 3.2% using standard parameters, outperforming other density functional theory (DFT)-based methods and machine learning methods (e.g., MetCCS). An online database is introduced for sharing both calculated and experimental CCS values ( metabolomics.pnnl.gov ), initially including a CCS library with over 1 million entries. Finally, three successful applications of molecule characterization using calculated CCS are described, including providing evidence for the presence of an environmental degradation product, the separation of molecular isomers, and an initial characterization of complex blinded mixtures of exposure chemicals. This work represents a method to address the limitations of small molecule identification and offers an alternative to generating chemical identification libraries experimentally by analyzing authentic reference materials. All code is available at github.com/pnnl .

SLIM Ultrahigh Resolution Ion Mobility Spectrometry Separations of Isotopologues and Isotopomers Reveal Mobility Shifts due to Mass Distribution Changes.

We report on separations of ion isotopologues and isotopomers using ultrahigh-resolution traveling wave-based Structures for Lossless Ion Manipulations with serpentine ultralong path and extended routing ion mobility spectrometry coupled to mass spectrometry (SLIM SUPER IMS-MS). Mobility separations of ions from the naturally occurring ion isotopic envelopes (e.g., [M], [M+1], [M+2], ... ions) showed the first and second isotopic peaks (i.e., [M+1] and [M+2]) for various tetraalkylammonium ions could be resolved from their respective monoisotopic ion peak ([M]) after SLIM SUPER IMS with resolving powers of ∼400-600. Similar separations were obtained for other compounds (e.g., tetrapeptide ions). Greater separation was obtained using argon versus helium drift gas, as expected from the greater reduced mass contribution to ion mobility described by the Mason-Schamp relationship. To more directly explore the role of isotopic substitutions, we studied a mixture of specific isotopically substituted (15N, 13C, and 2H) protonated arginine isotopologues. While the separations in nitrogen were primarily due to their reduced mass differences, similar to the naturally occurring isotopologues, their separations in helium, where higher resolving powers could also be achieved, revealed distinct additional relative mobility shifts. These shifts appeared correlated, after correction for the reduced mass contribution, with changes in the ion center of mass due to the different locations of heavy atom substitutions. The origin of these apparent mass distribution-induced mobility shifts was then further explored using a mixture of Iodoacetyl Tandem Mass Tag (iodoTMT) isotopomers (i.e., each having the same exact mass, but with different isotopic substitution sites). Again, the observed mobility shifts appeared correlated with changes in the ion center of mass leading to multiple monoisotopic mobilities being observed for some isotopomers (up to a ∼0.04% difference in mobility). These mobility shifts thus appear to reflect details of the ion structure, derived from the changes due to ion rotation impacting collision frequency or momentum transfer, and highlight the potential for new approaches for ion structural characterization.

New mass spectrometry technologies contributing towards comprehensive and high throughput omics analyses of single cells.

Mass-spectrometry based omics technologies - namely proteomics, metabolomics and lipidomics - have enabled the molecular level systems biology investigation of organisms in unprecedented detail. There has been increasing interest for gaining a thorough, functional understanding of the biological consequences associated with cellular heterogeneity in a wide variety of research areas such as developmental biology, precision medicine, cancer research and microbiome science. Recent advances in mass spectrometry (MS) instrumentation and sample handling strategies are quickly making comprehensive omics analyses of single cells feasible, but key breakthroughs are still required to push through remaining bottlenecks. In this review, we discuss the challenges faced by single cell MS-based omics analyses and highlight recent technological advances that collectively can contribute to comprehensive and high throughput omics analyses in single cells. We provide a vision of the potential of integrating pioneering technologies such as Structures for Lossless Ion Manipulations (SLIM) for improved sensitivity and resolution, novel peptide identification tactics and standards free metabolomics approaches for future applications in single cell analysis.

Evaluation of In Silico Multifeature Libraries for Providing Evidence for the Presence of Small Molecules in Synthetic Blinded Samples.

The current gold standard for unambiguous molecular identification in metabolomics analysis is comparing two or more orthogonal properties from the analysis of authentic reference materials (standards) to experimental data acquired in the same laboratory with the same analytical methods. This represents a significant limitation for comprehensive chemical identification of small molecules in complex samples. The process is time consuming and costly, and the majority of molecules are not yet represented by standards. Thus, there is a need to assemble evidence for the presence of small molecules in complex samples through the use of libraries containing calculated chemical properties. To address this need, we developed a Multi-Attribute Matching Engine (MAME) and a library derived in part from our in silico chemical library engine (ISiCLE). Here, we describe an initial evaluation of these methods in a blinded analysis of synthetic chemical mixtures as part of the U.S. Environmental Protection Agency's (EPA) Non-Targeted Analysis Collaborative Trial (ENTACT, Phase 1). For molecules in all mixtures, the initial blinded false negative rate (FNR), false discovery rate (FDR), and accuracy were 57%, 77%, and 91%, respectively. For high evidence scores, the FDR was 35%. After unblinding of the sample compositions, we optimized the scoring parameters to better exploit the available evidence and increased the accuracy for molecules suspected as present. The final FNR, FDR, and accuracy were 67%, 53%, and 96%, respectively. For high evidence scores, the FDR was 10%. This study demonstrates that multiattribute matching methods in conjunction with in silico libraries may one day enable reduced reliance on experimentally derived libraries for building evidence for the presence of molecules in complex samples.

Isolation of Tryptanthrin and Reassessment of Evidence for Its Isobaric Isostere Wrightiadione in Plants of the Wrightia Genus.

A series of Wrightia hanleyi extracts was screened for activity against Mycobacterium tuberculosis H37Rv. One active fraction contained a compound that initially appeared to be either the isoflavonoid wrightiadione or the alkaloid tryptanthrin, both of which have been previously reported in other Wrightia species. Characterization by NMR and MS, as well as evaluation of the literature describing these compounds, led to the conclusion that wrightiadione (1) was misidentified in the first report of its isolation from W. tomentosa in 1992 and again in 2015 when reported in W. pubescens and W. religiosa. Instead, the molecule described in these reports and in the present work is almost certainly the isobaric (same nominal mass) and isosteric (same number of atoms, valency, and shape) tryptanthrin (2), a well-known quinazolinone alkaloid found in a variety of plants including Wrightia species. Tryptanthrin (2) is also accessible synthetically via several routes and has been thoroughly characterized. Wrightiadione (1) has been synthesized and characterized and may have useful biological activity; however, this compound can no longer be said to be known to exist in Nature. To our knowledge, this misidentification of wrightiadione (1) has heretofore been unrecognized.

Improving network inference algorithms using resampling methods.

BACKGROUND: Relatively small changes to gene expression data dramatically affect co-expression networks inferred from that data which, in turn, can significantly alter the subsequent biological interpretation. This error propagation is an underappreciated problem that, while hinted at in the literature, has not yet been thoroughly explored. Resampling methods (e.g. bootstrap aggregation, random subspace method) are hypothesized to alleviate variability in network inference methods by minimizing outlier effects and distilling persistent associations in the data. But the efficacy of the approach assumes the generalization from statistical theory holds true in biological network inference applications. RESULTS: We evaluated the effect of bootstrap aggregation on inferred networks using commonly applied network inference methods in terms of stability, or resilience to perturbations in the underlying expression data, a metric for accuracy, and functional enrichment of edge interactions. CONCLUSION: Bootstrap aggregation results in improved stability and, depending on the size of the input dataset, a marginal improvement to accuracy assessed by each method's ability to link genes in the same functional pathway.

Salmonella-Mediated Inflammation Eliminates Competitors for Fructose-Asparagine in the Gut.

Salmonella enterica elicits intestinal inflammation to gain access to nutrients. One of these nutrients is fructose-asparagine (F-Asn). The availability of F-Asn to Salmonella during infection is dependent upon Salmonella pathogenicity islands 1 and 2, which in turn are required to provoke inflammation. Here, we determined that F-Asn is present in mouse chow at approximately 400 pmol/mg (dry weight). F-Asn is also present in the intestinal tract of germfree mice at 2,700 pmol/mg (dry weight) and in the intestinal tract of conventional mice at 9 to 28 pmol/mg. These findings suggest that the mouse intestinal microbiota consumes F-Asn. We utilized heavy-labeled precursors of F-Asn to monitor its formation in the intestine, in the presence or absence of inflammation, and none was observed. Finally, we determined that some members of the class Clostridia encode F-Asn utilization pathways and that they are eliminated from highly inflamed Salmonella-infected mice. Collectively, our studies identify the source of F-Asn as the diet and that Salmonella-mediated inflammation is required to eliminate competitors and allow the pathogen nearly exclusive access to this nutrient.

PIXiE: an algorithm for automated ion mobility arrival time extraction and collision cross section calculation using global data association.

MOTIVATION: Drift tube ion mobility spectrometry coupled with mass spectrometry (DTIMS-MS) is increasingly implemented in high throughput omics workflows, and new informatics approaches are necessary for processing the associated data. To automatically extract arrival times for molecules measured by DTIMS at multiple electric fields and compute their associated collisional cross sections (CCS), we created the PNNL Ion Mobility Cross Section Extractor (PIXiE). The primary application presented for this algorithm is the extraction of data that can then be used to create a reference library of experimental CCS values for use in high throughput omics analyses. RESULTS: We demonstrate the utility of this approach by automatically extracting arrival times and calculating the associated CCSs for a set of endogenous metabolites and xenobiotics. The PIXiE-generated CCS values were within error of those calculated using commercially available instrument vendor software. AVAILABILITY AND IMPLEMENTATION: PIXiE is an open-source tool, freely available on Github. The documentation, source code of the software, and a GUI can be found at https://github.com/PNNL-Comp-Mass-Spec/PIXiE and the source code of the backend workflow library used by PIXiE can be found at https://github.com/PNNL-Comp-Mass-Spec/IMS-Informed-Library . CONTACT: erin.baker@pnnl.gov or thomas.metz@pnnl.gov. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Structural and metabolic responses of Staphylococcus aureus biofilms to hyperosmotic and antibiotic stress.

Biofilms alter their metabolism in response to environmental stress. This study explores the effect of a hyperosmotic agent-antibiotic treatment on the metabolism of Staphylococcus aureus biofilms through the use of nuclear magnetic resonance (NMR) techniques. To determine the metabolic activity of S. aureus, we quantified the concentrations of metabolites in spent medium using high-resolution NMR spectroscopy. Biofilm porosity, thickness, biovolume, and relative diffusion coefficient depth profiles were obtained using NMR microimaging. Dissolved oxygen concentration was measured to determine the availability of oxygen within the biofilm. Under vancomycin-only treatment, the biofilm communities switched to fermentation under anaerobic condition, as evidenced by high concentrations of formate (7.4 ± 2.7 mM), acetate (13.1 ± 0.9 mM), and lactate (3.0 ± 0.8 mM), and there was no detectable dissolved oxygen in the biofilm. In addition, we observed the highest consumption of pyruvate (0.19 mM remaining from an initial 40 mM concentration), the sole carbon source, under the vancomycin-only treatment. On the other hand, relative effective diffusion coefficients increased from 0.73 ± 0.08 to 0.88 ± 0.08 under vancomycin-only treatment but decreased from 0.71 ± 0.04 to 0.60 ± 0.07 under maltodextrin-only and from 0.73 ± 0.06 to 0.56 ± 0.08 under combined treatments. There was an increase in biovolume, from 2.5 ± 1 mm3 to 7 ± 1 mm3 , under the vancomycin-only treatment, while the maltodextrin-only and combined treatments showed no significant change in biovolume over time. This indicated that physical biofilm growth was halted during maltodextrin-only and combined treatments.

Enhancing glycan isomer separations with metal ions and positive and negative polarity ion mobility spectrometry-mass spectrometry analyses.

Glycomics has become an increasingly important field of research since glycans play critical roles in biology processes ranging from molecular recognition and signaling to cellular communication. Glycans often conjugate with other biomolecules, such as proteins and lipids, and alter their properties and functions, so glycan characterization is essential for understanding the effects they have on cellular systems. However, the analysis of glycans is extremely difficult due to their complexity and structural diversity (i.e., the number and identity of monomer units, and configuration of their glycosidic linkages and connectivities). In this work, we coupled ion mobility spectrometry with mass spectrometry (IMS-MS) to characterize glycan standards and biologically important isomers of synthetic αGal-containing O-glycans including glycotopes of the protozoan parasite Trypanosoma cruzi, which is the causative agent of Chagas disease. IMS-MS results showed significant differences for the glycan structural isomers when analyzed in positive and negative polarity and complexed with different metal cations. These results suggest that specific metal ions or ion polarities could be used to target and baseline separate glycan isomers of interest with IMS-MS. Graphical abstract Glycan isomers, such as fructose and glucose, show distinct separations in positive and negative ion mode.

Toward high-resolution NMR spectroscopy of microscopic liquid samples.

A longstanding limitation of high-resolution NMR spectroscopy is the requirement for samples to have macroscopic dimensions. Commercial probes, for example, are designed for volumes of at least 5 µL, in spite of decades of work directed toward the goal of miniaturization. Progress in miniaturizing inductive detectors has been limited by a perceived need to meet two technical requirements: (1) minimal separation between the sample and the detector, which is essential for sensitivity, and (2) near-perfect magnetic-field homogeneity at the sample, which is typically needed for spectral resolution. The first of these requirements is real, but the second can be relaxed, as we demonstrate here. By using pulse sequences that yield high-resolution spectra in an inhomogeneous field, we eliminate the need for near-perfect field homogeneity and the accompanying requirement for susceptibility matching of microfabricated detector components. With this requirement removed, typical imperfections in microfabricated components can be tolerated, and detector dimensions can be matched to those of the sample, even for samples of volume ≪5 µL. Pulse sequences that are robust to field inhomogeneity thus enable small-volume detection with optimal sensitivity. We illustrate the potential of this approach to miniaturization by presenting spectra acquired with a flat-wire detector that can easily be scaled to subnanoliter volumes. In particular, we report high-resolution NMR spectroscopy of an alanine sample of volume 500 pL.

High-resolution microstrip NMR detectors for subnanoliter samples.

We present the numerical optimization and experimental characterization of two microstrip-based nuclear magnetic resonance (NMR) detectors. The first detector, introduced in our previous work, was a flat wire detector with a strip resting on a substrate, and the second detector was created by adding a ground plane on top of the strip conductor, separated by a sample-carrying capillary and a thin layer of insulator. The dimensional parameters of the detectors were optimized using numerical simulations with regards to radio frequency (RF) sensitivity and homogeneity, with particular attention given to the effect of the ground plane. The influence of copper surface finish and substrate surface on the spectral resolution was investigated, and a resolution of 0.8-1.5 Hz was obtained on 1 nL deionized water depending on sample positioning. For 0.13 nmol sucrose (0.2 M in 0.63 nL H2O) encapsulated between two Fluorinert plugs, high RF homogeneity (A810°/A90° = 70-80%) and high sensitivity (expressed in the limit of detection nLODm = 0.73-1.21 nmol s1/2) were achieved, allowing for high-performance 2D NMR spectroscopy of subnanoliter samples.