an increasing number of scientific investigations involve the analysis of large sample sets , often assembled in a range of divergent compositions .
one of the best methods for atomic - level characterization of molecules and mixtures is solution - state nuclear magnetic resonance ( nmr ) spectroscopy ( maher et al .
the unparalleled capabilities of nmr to acquire useful data , however , require a non - trivial level of expertise with an nmr spectrometer and familiarity with its underlying principles . to set up even the simplest one - dimensional experiments
for example , the spectrometer probe must be  tuned and matched  for each new sample placed in the spectrometer to maximize the efficiency of radio - frequency ( rf ) signals sent and received from that sample .
in addition , the magnetic field passing through the sample needs to be made as homogeneous as possible in order to optimize spectral lineshapes .
this process , called shimming , is achieved by adjusting the electrical current in a multitude of  shim  coils directly adjacent to the sample .
also , the pulse program that dictates which nuclei are probed for each experiment contains radio frequency ( rf ) 
pulses  that need to be calibrated for optimal signal - to - noise ( s / n ) .
other parameters , such as the range of frequencies to sample ( spectral - width , or sw ) , are difficult to determine a priori , and must be manually deduced for subsequent data collections .
hence , experiment setups for a large set of samples can potentially consume a significant portion of an investigator s time and effort . in recent years
bruker biospin has introduced several hardware and software products that soften the requirements for user technical expertise and promote high - throughput nmr spectroscopy .
hardware accessories for automated probe tuning and matching ( atm ) and sample - tube changers ( samplejet ) provide users the convenience of manipulating the probe or sample , respectively , from the newest versions of topspin , bruker s software for nmr data acquisition and analysis ( soininen et al . 2009 ) .
topspin provides the interface for these hardware accessories , as well as automated procedures for sample shimming , pulse calibrations , and receiver gain optimization .
topspin includes a legacy software suite , icon - nmr , for high - throughput data acquisition that incorporates many of the software features described previously . in attempting to perform several studies at the national magnetic resonance facility at madison ( nmrfam ) involving large sample sets ( e.g. metabolomics , crude extracts , small - molecule libraries , protein  ligand screening ) , we encountered limitations to icon - nmr that hinder accessibility , preclude sample set heterogeneity , and limit the quality of acquired data . the specific limitations in icon - nmr include : a restrictive interface with access separate from topspin ; a complicated and inflexible menu system for sample information entry ; reduced performance of automated sample shimming ; difficulty in accommodating different solvents in the same sample set ; and an inability to adapt experiment parameters to each sample .
these limitations provided the impetus to develop a more straightforward , intuitive , high - throughput methodology for automated data acquisition across diverse sample sets .
our goals were to simplify the setup procedure for data acquisition , provide easily customizable functionality , and improve the quality of data acquired over what was previously obtained from icon - nmr .
automated , high - throughput nmr data collection for large sample sets first requires access to a spectrometer equipped with the automated hardware and software features described in the introduction . for our purposes
, we developed and tested this new methodology on four bruker biospin nmr spectrometers at nmrfam ; a 500  mhz avance iii with 5  mm triple resonance cryoprobe , a 600  mhz avance iii with 5  mm quadruple resonance cryoprobe , a 600  mhz avance iii with 1.7  mm triple resonance cryoprobe , and a 700  mhz avance iii with 5  mm quadruple resonance cryoprobe .
each spectrometer was equipped with samplejet and atm accessories , running topspin v. 3.0 under centos 5 . to develop our methodology we utilized the python programming language ( conway 1995 ) interpreter recently added in version 2.0 of topspin .
it should be noted that we began development of nmrbot using topspin v. 2.0 , hence nmrbot is backward compatible with this earlier version of topspin .
the interpreter currently accepts functions from the python v. 2.7.3 standard library and a number of modules designed by bruker to access specific spectrometer functions .
methodology development focused on three areas : ( 1 ) design of an intuitive and flexible user interface for entry of sample information and experimental parameters , ( 2 ) automated operation of the spectrometer and cursory spectral analysis , and ( 3 ) sample data tracking and error handling . to provide users with a more straightforward and flexible interface to setup automated data acquisition , we developed two approaches to input sample and experiment information .
execution of the python script starts the setup wizard that provides access to both approaches ( fig .  1 ) .
as such , the first window of the setup wizard asks users to select between manual or text file modes for input .
the manual input mode involves a series of input windows prompting the user for pertinent details ( fig .  1 , left ) .
a python function was developed for each window , designed to verify user input or identify input errors .
alternatively , the text file input mode circumvents the manual input windows and prompts the user to enter the name of a text file that contains sample and experiment details ( fig .  1 , right ) .
the details in this file must be enumerated in the self - defining text archive and retrieval ( star ) format ( hall 1991 ) , with sample specific information denoted separately from optional folder format and series parameters . a python function to decode this type of text file was designed to validate the inputs and identify potential errors before acquisition begins . for both input modes ,
requisite input details are the number of samples , names , solvents , rack position , and experiment parameter set names ( fig .  1 , white boxes in center ) .
for experiment details , our method relies on predefined experiment parameter sets , a convention employed by bruker to easily recall all spectrometer parameters for a specific experiment .
any number of available optional inputs , as described below , for customized optimization and operation of the spectrometer can be easily appended in the initialization script.fig . 
flowchart of the setup wizard ( left ) user interface showing start and end points ( diamonds ) , requisite inputs for manual path ( white boxes ) , optional manual inputs ( shaded boxes ) , and optional course for text file input ( right ) .
the expansion shows an example of an input text file in star format ( hall 1991 ) containing relevant sample and experiment parameters nmrbot sample and experiment parameter input methods .
flowchart of the setup wizard ( left ) user interface showing start and end points ( diamonds ) , requisite inputs for manual path ( white boxes ) , optional manual inputs ( shaded boxes ) , and optional course for text file input ( right ) .
the expansion shows an example of an input text file in star format ( hall 1991 ) containing relevant sample and experiment parameters the optional inputs are described here in the order they appear in the manual input mode of the setup wizard ( fig .  1 , left ) .
the option to include a water standard sample ( 90  % h2o , 10  % d2o ) is offered , and is the first to appear to differentiate from other samples in the series .
standard sample to establish optimal three - dimensional sample shims with topshim , the automated shimming routine in topspin .
after the requisite sample set information described above is collected , further optional inputs are presented .
the first two provide simple bookkeeping preferences for sample folder name format and sample condition notes .
subsequent input options allow the user to tailor certain acquisition parameters for specific samples or the entire set .
ideal for low concentration samples in the sample series , an option is offered to multiply the number of scans ( ns ) for every nmr experiment of that sample by a given factor .
the next optional input allows the user to define which one - dimensional experiments should be used to identify the range of observable sample peaks , as well as to define any subsequent experiments that should have their sw parameters adapted to this range along with a change of the spectrum center ( offset ) corresponding to the center of this range .
another option asks if the user wants common rf pulse parameters ( length and power ) to be automatically loaded from the  prosol  table , another bruker convention in topspin , rather than using those in the loaded parameter set .
if a sample in the series has a large solvent peak such as h2o , the user can opt to have the position of that peak automatically determined and entered as the offset .
the next option can be selected to automatically tune and match the spectrometer for each new sample in the series .
another optional input allows the user to define the desired topshim command used to optimize the shims for each sample .
a final optional input allows the user to toggle the use of the  rga  function in topspin that automatically optimizes the receiver gain parameter for each experiment .
if setup details are manually entered into the setup wizard , those details are compiled and output as a star formatted text file .
this file , or altered versions thereof , can be used as text file input to future executions of nmrbot .
once all information has been input into the setup wizard , either manually or by text file , a final window allows the user to review the parameters entered . upon confirmation
we have developed several python functions to automate the process of data acquisition and analysis .
many of these functions rely on core spectrometer interface functions included in the topspin python interpreter .
the most commonly utilized function in our methodology passes commands directly to the command line of topspin .
this allows our methodology to largely follow the standard series of commands for manual operation of the spectrometer . in this way samples
are inserted and shimmed , a deuterium lock is established , the probe is tuned and matched , sample data folders are created , and experiment parameters sets are loaded .
these steps can be circumvented or modified by any optional inputs submitted by the user , such as skipping the tune / match step , changing experiment parameters from prosol or ns options , or optimizing the receiver gain . if the user opts to find the solvent peak for a sample , a one - dimensional h nmr experiment is loaded separately and a spectrum is acquired , followed by automated peak - picking and peak analysis to determine the position of the large solvent peak .
this value is then passed internally to all experiments loaded for that sample to define the h dimension offset .
the function developed for peak list analysis is also used to determine optimum sw ( and corresponding offset in the center of this range ) for any spectral dimension , if the user included this option .
this information is passed along to any experiments flagged by the user for adapted sw .
as with the setup wizard input functions , the automated data acquisition functions are designed to validate each step of data acquisition ; if an error is encountered , it is logged and the acquisition proceeds to the next experiment in the list . during development and testing , we found it advantageous to audit the progress of sample setups and automated data acquisition .
this feature provides an accounting of each step of the method with real - time updates of the software s activity displayed in the terminal window associated with topspin .
also , these updates are appended to an audit text file output in each sample s data directory , and all updates of the sample series are likewise output to the user s experiment directory .
we made every effort to design the acquisition functions so that they would continue on to the next procedure upon error detection .
all the functions are combined in a single script file , named fam_tools.py , and placed in the directory < topspinhome>/exp / stan
the same directory contains a short script file , named nmrbot.py , which begins the initial process of information collection from the user .
every session is invoked by typing the name of the short script file in the command line of topspin . for direct comparison
, identical nmr experiments were acquired on several complex mixture samples using the two automated methods , icon - nmr and nmrbot .
the h 1d , c 1d and h - c 2d hsqc experiments used the same parameter sets and shimming routine .
an additional 2d hsqc was acquired with nmrbot , employing the adaptive spectral - width feature to automatically determine the optimal c dimension sw parameter using peaks observed in the c 1d experiment .
the nmrbot method was able to reproduce the basic behavior of icon - nmr for automated data acquisition on a series of samples .
for example , the time to acquire the data using the icon - nmr and nmrbot methods was comparable , however the setup of nmrbot required less time from the user and several novel nmrbot methods did add small amounts of time in certain circumstances , as described below .
the setup time saving for nmrbot as compared to icon - nmr varied between about 5  min ( for manual input ) and 20  min ( for text file input ) .
user feedback during development indicated that the two nmrbot input methods , especially the text file input method , were found to be more accessible and straightforward than the icon - nmr interface .
users attributed this to two main differences of the input interface : ( 1 ) the successive prompts of the nmrbot setup wizard present the user with all pertinent setup variables , reducing the chance of a missed input and distilling the functional complexity of the interface ; ( 2 ) eliminating the requirement for familiarity with drop - down menus and other non - obvious input methods , as are employed in icon - nmr , thus lowering the learning - curve to use nmrbot and reducing the cognitive complexity of the interface .
in fact , the nmrbot text file input mode obviates most of the setup wizard interface .
several important differences were noted in the performance of the two automated methodologies ( see table  1 ) .
first , icon - nmr displayed difficulty in automatically determining optimum solvent lock parameters for samples using different solvents than proceeding ones in the series .
presumably , icon - nmr uses the same software protocol as nmrbot , which is to call the  lock 
procedure as is done manually , so this difference in performance occurred for reasons unknown .
show a clear improvement in resolution for the nmrbot acquired spectrum , and hence allowed more complete analysis of this sample .
again , this difference in performance is inexplicable given that both methods rely on the topshim procedure in topspin .
these performance limitations for icon - nmr impacted the quality of acquired spectra , potentially wasting spectrometer time and requiring manual reacquisition .
unusable  when acquired with icon - nmr . all spectra collected thus far with the fully developed nmrbot method have met quality criteria for each user.table  1qualitative assessment of nmrbot features as compared to icon - nmrfeaturenmrboticon - nmrt ( min)user interface5 , 20sample shimmingprobe tuningsample handlingadapted spectral - widthx  scan multiplierx 
shimmingx<10text file audit trailx<1stars in the method columns indicate the presence of a feature and , if applicable , the number of stars indicates feature performance as determined by nmrbot user feedback .
the time difference ( t ) column indicates any nmrbot feature time difference as compared to icon - nmrtime savings from nmrbot manual or text file input methodswater standard sample shimming is lengthy , but can reduce time for later procedures ( see text)fig . 
2h nmr spectra of a complex mixture collected by two automated methods , ( bottom ) icon - nmr and ( top ) nmrbot , which use the same automated shimming method ( topshim ) .
the spectrum shimmed under the nmrbot protocol shows slightly better resolution qualitative assessment of nmrbot features as compared to icon - nmr stars in the method columns indicate the presence of a feature and , if applicable , the number of stars indicates feature performance as determined by nmrbot user feedback .
the time difference ( t ) column indicates any nmrbot feature time difference as compared to icon - nmr time savings from nmrbot manual or text file input methods water standard sample shimming is lengthy , but can reduce time for later procedures ( see text ) h nmr spectra of a complex mixture collected by two automated methods , ( bottom ) icon - nmr and ( top ) nmrbot , which use the same automated shimming method ( topshim ) .
the spectrum shimmed under the nmrbot protocol shows slightly better resolution several features developed for nmrbot provided other distinct advantages over icon - nmr . the option to initially perform three - dimensional shimming on a water standard sample
did add up to 10 additional minutes to the overall acquisition time of nmrbot , but with a slight reduction in shim times and improved resolution for subsequent samples in the series .
other optional inputs allowed specific parameters to be automatically determined and modified on - the - fly , increasing data quality and consistency .
for example , the input ns multiple improved the s / n for all data acquired on low concentration samples using the same parameter sets as other more concentrated samples in the series . for samples with large solvent peaks that needed to be suppressed by presaturation , the ability of nmrbot to automatically determine the offset allowed a general parameter set to be used and did not require the investigator s time to predetermine the offset .
the method to automatically determine the optimum offset , however , did add up to 60  s to the overall method per sample .
moreover , use of the adaptive spectral - width feature enhanced data resolution and promoted more effective decoupling in certain experiments .
this advantage is shown in fig .  3 , which compares spectra initiated from the same parameter set , one with default settings and the other with the adaptive spectral - width feature enabled for the c dimension .
the time required for automated peaklist analysis to determine adapted spectra - width parameters was negligible ( 12  s ) , and slightly increased total acquisition time when applied to the direct dimension of subsequent experiments.fig . 
the 2d hsqc spectrum on the left was collected with general acquisition parameters ( i.e. c - dimension sw parameter of 200  ppm ) .
the spectrum on the right was collected using the same acquisition parameters , except for a 50  % reduction in spectral - width and a corresponding change in c - offset ( dimension center ) as determined by automated analysis of a 1d c spectrum acquired previously in the experiment set .
the gains in c resolution and decoupling efficiency are highlighted by the trace along the left edge of each spectrum improvements in 2d h- c hsqc spectral quality due to adaptive spectral - width .
the 2d hsqc spectrum on the left was collected with general acquisition parameters ( i.e. c - dimension sw parameter of 200  ppm ) .
the spectrum on the right was collected using the same acquisition parameters , except for a 50  % reduction in spectral - width and a corresponding change in c - offset ( dimension center ) as determined by automated analysis of a 1d c spectrum acquired previously in the experiment set .
the gains in c resolution and decoupling efficiency are highlighted by the trace along the left edge of each spectrum the creation of audit trails for each sample and the entire sample series provides investigators with an accounting of all spectrometer activity and a means of validating data acquisition procedures .
this feature is also useful for discerning the time and nature of acquisition errors , helping nmrbot users and developers alike .
other convenience features , such as folder name format , condition notes , and automated spectrum title details , aid users in organizing and tracking acquired data .
our motivation for developing nmrbot as a custom - built application arose from perceived limitations in the available automated data acquisition software for bruker nmr spectrometers .
its initial development was meant to bypass these limitations , but we quickly determined that other improvements would assist the needs of investigators at nmrfam .
nmrbot provides an accessible , robust , time - saving setup interface for spectrometer users of all stripes .
optional features expand the functionality of automated acquisition and further save investigator time by automating the determination of several parameters that enhance data quality and consistency .
this methodology currently stands as an alternative to bruker s icon - nmr . in future releases of nmrbot
, we plan to include features such as series completion time calculations and a run - time interface for more user control during acquisition .
the python scripts for nmrbot are currently available from the nmrfam website ( www.nmrfam.wisc.edu/software/nmrbot/ ) , along with simple installation and usage instructions .
 to facilitate the high - throughput acquisition of nuclear magnetic resonance ( nmr ) experimental data on large sets of samples , we have developed a simple and straightforward automated methodology that capitalizes on recent advances in bruker biospin nmr spectrometer hardware and software . given the daunting challenge for non - nmr experts to collect quality spectra , our goal was to increase user accessibility , provide customized functionality , and 
 improve the consistency and reliability of resultant data . 
 this methodology , nmrbot , is encoded in a set of scripts written in the python programming language accessible within the bruker biospin topspin software . 
 nmrbot improves automated data acquisition and offers novel tools for use in optimizing experimental parameters on the fly . 
 this automated procedure has been successfully implemented for investigations in metabolomics , small - molecule library profiling , and protein  ligand titrations on four bruker biospin nmr spectrometers at the national magnetic resonance facility at madison . 
 the investigators reported benefits from ease of setup , improved spectral quality , convenient customizations , and overall time savings . 
Introduction
Methods
Results and discussion
Concluding remarks