Potentiality of proton nuclear magnetic resonance and multivariate calibration methods for the determination of dermatan sulfate contamination in heparin samples

A 1H NMR method for the quantification of dermatan sulfate impurities in heparin industrial samples is proposed. The method is based on the analysis of 1H NMR spectral data by multivariate calibration. The 1H NMR spectra of heparin and dermatan sulfate standards showed characteristic profiles. Thus, differences in the methyl peaks of acetamido groups of heparin and dermatan sulfate were greatly advantageous for the analysis. Other hydrogens of the sugar ring were also relevant in this study. Thus, the determination of dermatan sulfate by multivariate calibration depended on all these differences. Partial least squares regression (PLS) was chosen as the calibration method. In addition, a data standardization procedure was developed in order that 1H NMR spectra registered with different instruments operating under different measurement conditions were comparable. The quantification of dermatan sulfate in the samples was satisfactory, with an overall prediction error of 6%.     

High-sensitivity visualisation of contaminants in heparin samples by spectral filtering of 1H NMR spectra

A novel application of two-dimensional correlation analysis has been employed to filter (1)H NMR heparin spectra distinguishing acceptable natural variation and the presence of foreign species. Analysis of contaminated heparin samples, compared to a dataset of accepted heparin samples using two-dimensional correlation spectroscopic analysis of their 1-dimensional (1)H NMR spectra, allowed the spectral features of contaminants to be recovered with high sensitivity, without having to resort to more complicated NMR experiments. Contaminants, which exhibited features distinct from those of heparin and those with features normally hidden within the spectral mass of heparin could be distinguished readily. A heparin sample which had been pre-mixed with a known contaminant, oversulfated chondroitin sulfate (OSCS), was tested against the heparin reference library. It was possible to recover the (1)H NMR spectrum of the OSCS component through difference 2D-COS power spectrum analysis of as little as 0.25% (w/w) with ease, and of 2% (w/w) for more challenging contaminants, whose NMR signals fell under those of heparin. The approach shows great promise for the quality control of heparin and provides the basis for greatly improved regulatory control for the analysis of heparin, as well as other intrinsically heterogeneous and varied products.     

Identification of heparin samples that contain impurities or contaminants by chemometric pattern recognition analysis of proton NMR spectral data

Chemometric analysis of a set of one-dimensional (1D) (1)H nuclear magnetic resonance (NMR) spectral data for heparin sodium active pharmaceutical ingredient (API) samples was employed to distinguish USP-grade heparin samples from those containing oversulfated chondroitin sulfate (OSCS) contaminant and/or unacceptable levels of dermatan sulfate (DS) impurity. Three chemometric pattern recognition approaches were implemented: classification and regression tree (CART), artificial neural network (ANN), and support vector machine (SVM). Heparin sodium samples from various manufacturers were analyzed in 2008 and 2009 by 1D (1)H NMR, strong anion-exchange high-performance liquid chromatography, and percent galactosamine in total hexosamine tests. Based on these data, the samples were divided into three groups: Heparin, DS ≤ 1.0% and OSCS = 0%; DS, DS > 1.0% and OSCS = 0%; and OSCS, OSCS > 0% with any content of DS. Three data sets corresponding to different chemical shift regions (1.95-2.20, 3.10-5.70, and 1.95-5.70 ppm) were evaluated. While all three chemometric approaches were able to effectively model the data in the 1.95-2.20 ppm region, SVM was found to substantially outperform CART and ANN for data in the 3.10-5.70 ppm region in terms of classification success rate. A 100% prediction rate was frequently achieved for discrimination between heparin and OSCS samples. The majority of classification errors between heparin and DS involved cases where the DS content was close to the 1.0% DS borderline between the two classes. When these borderline samples were removed, nearly perfect classification results were attained. Satisfactory results were achieved when the resulting models were challenged by test samples containing blends of heparin APIs spiked with non-, partially, or fully oversulfated chondroitin sulfate A, heparan sulfate, or DS at the 1.0%, 5.0%, and 10.0% (w/w) levels. This study demonstrated that the combination of 1D (1)H NMR spectroscopy with multivariate chemometric methods is a nonsubjective, statistics-based approach for heparin quality control and purity assessment that, once standardized, minimizes the need for expert analysts.     

Improving prediction selectivity for on-line near-infrared monitoring of components in etchant solution by spectral range optimization

The components (H₃PO₄, HNO₃, CH₃COOH and water) in an etchant solution have been accurately measured in an on-line manner using near-infrared (NIR) spectroscopy by directly illuminating NIR radiation through a Teflon line. In particular, the spectral features according to the change of H₃PO₄ or HNO₃ concentrations were not mainly from NIR absorption themselves, but from the perturbation (or displacement) of water bands; therefore, the resulting spectral variations were quite similar to each other. Consequently partial least squares (PLS) prediction selectivity among the components should be the most critical issue for continuous on-line compositional monitoring by NIR spectroscopy. To improve selectivity of the calibration model, we have optimized the calibration models by finding selective spectral ranges with the use of moving window PLS. Using the optimized PLS models for each component, the resulting prediction accuracies were substantially improved. Furthermore, on-line prediction selectivity was evaluated by spiking individual pure components step by step and examining the resulting prediction trends. When optimized PLS models were used, each concentration was selectively and sensitively varied at each spike; meanwhile, when whole or non-optimized ranges were used for PLS, the prediction selectivity was greatly degraded. This study verifies that the selection of an optimal spectral range for PLS is the most important factor to make Teflon-based NIR measurements successful for on-line and real-time monitoring of etching solutions.     

Partial least squares modeling of combined infrared, H NMR and C NMR spectra to predict long residue properties of crude oils

Research has been carried out to determine the potential of partial least squares (PLS) modeling of mid-infrared (IR) spectra of crude oils combined with the corresponding ¹H and ¹³C nuclear magnetic resonance (NMR) data, to predict the long residue (LR) properties of these substances. The study elaborates further on a recently developed and patented method to predict this type of information from only IR spectra. In the present study, PLS modeling was carried out for 7 different LR properties, i.e., yield long-on-crude (YLC), density (D_LR), viscosity (V_LR), sulfur content (S), pour point (PP), asphaltenes (Asph) and carbon residue (CR). Research was based on the spectra of 48 crude oil samples of which 28 were used to build the PLS models and the remaining 20 for validation. For each property, PLS modeling was carried out on single type IR, ¹³C NMR and ¹H NMR spectra and on 3 sets of merged spectra, i.e., IR + ¹H NMR, IR + ¹³C NMR and IR + ¹H NMR + ¹³C NMR. The merged spectra were created by considering the NMR data as a scaled extension of the IR spectral region. In addition, PLS modeling of coupled spectra was performed after a Principal Component Analysis (PCA) of the IR, ¹³C NMR and ¹H NMR calibration sets. For these models, the 10 most relevant PCA scores of each set were concatenated and scaled prior to PLS modeling. The validation results of the individual IR models, expressed as root-mean-square-error-of-prediction (RMSEP) values, turned out to be slightly better than those obtained for the models using single input ¹³C NMR or ¹H NMR data. For the models based on IR spectra combined with NMR data, a significant improvement of the RMSEP values was not observed neither for the models based on merged spectra nor for those based on the PCA scores. It implies, that the commonly accepted complementary character of NMR and IR is, at least for the crude oil and bitumen samples under study, not reflected in the results of PLS modeling. Regarding these results, the absence of sample preparation and the straightforward way of data acquisition, IR spectroscopy is preferred over NMR for the prediction of LR properties of crude oils at site.     

Highly routinely reproducible alignment of 1H NMR spectral peaks of metabolites in huge sets of urines

A method to obtain high reproducibility of (1)H NMR chemical shift of peaks of biofluid metabolites, by simple acidification with HCl is evaluated. Biofluid (1)H NMR analysis is indeed spoiled by a strong chemical shift dependence of metabolite peaks on parameters such as ionic strength, concentration of some earth alkali cations and, mostly, on pH of samples. The resulting chemical shift variations, as large as 0.1 ppm, generate misalignments of homogeneous peaks, artifacts and misinterpretations. Reproducible alignment is essential in (1)H NMR based metabonomics, where peak misalignments prevent even very wide bins (i.e., 0.04 ppm, as elsewhere proposed) from being used to integrate spectral data for multivariate statistical analysis. Here is demonstrated that routine acidification with HCl to 1.2≤pH≤2.0 ensures highly reproducible peak alignment of urine (1)H NMR spectra. In this respect, simple inspection of citrate peaks in the urine can be used to measure pH, as it will be extensively discussed, in that at such low pH they show no dependency on other urine components as reported at higher pH. Under these conditions, in as many as 493 urine samples, in which concentrations of Ca(2+), Mg(2+), K(+), Na(+), Cl(-), phosphate, and creatinine and ionic strength measured by means of well standardized conventional procedures, showed very wide ranges, peaks align within a SD always lower than 0.002 ppm, thus allowing the use of integration bins at least five times narrower than 0.04 ppm.     

Comparison of Artificial Neural Networks with Partial Least Squares Regression for Simultaneous Determinations by ICP-AES

Simultaneous determination of several elements （U, Ta, Mn, Zr and W） with inductively coupled plasma atomic emission spectrometry （ICP-AES） in the presence of spectral interference was performed using chemometrics methods. True comparison between artificial neural network （ANN） and partial least squares regression （PLS） for simultaneous determination in different degrees of overlap was investigated. The emission spectra were recorded at uranium analytical line （263.553 nm） with a 0.06 nm spectral window by ICP-AES. Principal component analysis was applied to data and scores on 5 dominant principal components were subjected to ANN. A 5-5-5 （input, hidden and output neurons） network was used with linear transfer function after both hidden and output layers. The PI,S model was trained with five latent variables and 20 samples in calibration set. The relative errors of predictions （REP） in test set were 3.75% and 3.56% for ANN and PLS respectively.     

Partial least square analysis of lysozyme near-infrared spectra

Proteins possess strong absorption features in the combination range (5000-4000 cm⁻¹) of the near infrared (NIR) spectrum. These features can be used for quantitative analysis. Partial least squares (PLS) regression was used to analyze NIR spectra of lysozyme with the leave-one-out, full cross-validation method. A strategy for spectral range optimization with cross-validation PLS calibration was presented. A five-factor PLS model based on the spectral range between 4720 and 4540 cm⁻¹ provided the best calibration model for lysozyme in aqueous solutions. For 47 samples ranging from 0.01 to 10 mg/mL, the root mean square error of prediction was 0.076 mg/mL. This result was compared with values reported in the literature for protein measurements by NIR absorption spectroscopy in human serum and animal cell culture supernatants.     

Partial least squares modeling of combined infrared, 1H NMR and 13C NMR spectra to predict long residue properties of crude oils

Research has been carried out to determine the potential of partial least squares (PLS) modeling of mid-infrared (IR) spectra of crude oils combined with the corresponding ¹H and ¹³C nuclear magnetic resonance (NMR) data, to predict the long residue (LR) properties of these substances. The study elaborates further on a recently developed and patented method to predict this type of information from only IR spectra. In the present study, PLS modeling was carried out for 7 different LR properties, i.e., yield long-on-crude (YLC), density (D_LR), viscosity (V_LR), sulfur content (S), pour point (PP), asphaltenes (Asph) and carbon residue (CR). Research was based on the spectra of 48 crude oil samples of which 28 were used to build the PLS models and the remaining 20 for validation. For each property, PLS modeling was carried out on single type IR, ¹³C NMR and ¹H NMR spectra and on 3 sets of merged spectra, i.e., IR + ¹H NMR, IR + ¹³C NMR and IR + ¹H NMR + ¹³C NMR. The merged spectra were created by considering the NMR data as a scaled extension of the IR spectral region. In addition, PLS modeling of coupled spectra was performed after a Principal Component Analysis (PCA) of the IR, ¹³C NMR and ¹H NMR calibration sets. For these models, the 10 most relevant PCA scores of each set were concatenated and scaled prior to PLS modeling. The validation results of the individual IR models, expressed as root-mean-square-error-of-prediction (RMSEP) values, turned out to be slightly better than those obtained for the models using single input ¹³C NMR or ¹H NMR data. For the models based on IR spectra combined with NMR data, a significant improvement of the RMSEP values was not observed neither for the models based on merged spectra nor for those based on the PCA scores. It implies, that the commonly accepted complementary character of NMR and IR is, at least for the crude oil and bitumen samples under study, not reflected in the results of PLS modeling. Regarding these results, the absence of sample preparation and the straightforward way of data acquisition, IR spectroscopy is preferred over NMR for the prediction of LR properties of crude oils at site.     

Feasibility of monitoring acetic acid process using near-infrared spectroscopy

Near-infrared (NIR) spectroscopy has been utilized to demonstrate its feasibility for the measurement of major components in the acetic acid process. In order to simulate the acetic acid process, synthetic mixtures were prepared from five different components: acetic acid, methyl acetate, methyl iodide, water, and potassium iodide. Partial least squares (PLS) regression was utilized to differentiate the spectral characteristics as well as to quantify each component for the mixtures. The spectral features of acetic acid, methyl acetate, methyl iodide, and water are noticeably different with each other over the entire NIR region. The quantity of iodide ion, which does not absorb NIR radiation, was determined using the wavelength shift and intensity change of water absorption band caused by the change of iodide ion concentration. The PLS calibration results of the five components show good correlation with reference data. They also demonstrate the technical feasibility of NIR spectroscopy for monitoring important components in the acetic acid process.     

Comparison of near-infrared and Raman spectroscopy for on-line monitoring of etchant solutions directly through a Teflon tube

Both near-infrared (NIR) and Raman spectroscopy have been studied for the quantitative measurement of components (H(3)PO(4), HNO(3), and CH(3)COOH) in an etchant solution and the corresponding prediction robustness has been evaluated. Both measurements were accomplished by illuminating radiation directly through a Teflon tube. Raman spectral features of each component were much clearer and more selective than those observed in the NIR spectrum. Especially, NIR spectral variation pertinent to H(3)PO(4) and HNO(3) were mostly based on the displacement and perturbation of water bands rather than due solely to NIR absorption. Therefore, the resulting spectral variations were not highly specific. When the validation set contained minor spectral variations resulting from a slight instrumental change, NIR prediction performance for all three components degraded substantially by showing obvious prediction bias. However, the accuracies of Raman predictions were maintained. Since partial least squares (PLS) models for each component were built using NIR spectra of poor specificity with broadly overlapping features, even minor spectral differences introduced by instrumental variations sensitively influenced the prediction performance of the NIR models. Overall, the selectivity (specificity) of a targeting spectroscopic method should be considered critically to secure prediction robustness for monitoring components in an etchant solution.     

近红外光谱结合膜富集技术测定饮料中微量邻苯二甲酸二异辛酯的含量

提出了一种基于在线膜富集的近红外漫反射光谱技术,对饮料中的微量塑化剂邻苯二甲酸二异辛酯(DEHP)进行快速检测。采用聚醚砜膜对饮料中的DEHP进行富集,将富集DEHP的膜直接进行近红外漫反射检测。参考DEHP的透射近红外光谱,对波数进行选择,以4 420～4 060、4 700～4 540、6 040～5 600cm-1作为建模的波数区间。通过比较原始光谱、多元散射校正、一阶求导、二阶求导及其组合,考察了光谱预处理方法对模型的影响,用去一交互验证法建立了偏最小二乘(PLS)模型,并用所建立的校正模型对校正集样品进行了预测。结果表明,在选定的波数区间,当用一阶求导对校正集光谱进行预处理时,所建立的模型对校正集的预测效果最佳,在隐变量数为7时,对校正集所有样品的校正均方根误差(RMSEC)为0.188 7mg/L。用此模型对预测集样品进行预测时,DEHP的质量浓度在0.5～5.0 mg/L范围内,预测均方根误差(RMSEP)为0.232 4 mg/L,平均相对预测误差为6.29%。     

Determination of ethanol in fuel ethanol and beverages by Fourier transform (FT)-near infrared and FT-Raman spectrometries

Fourier transform-near infrared (FT-NIR) and FT-Raman spectrometries have been used to design partial least squares (PLS) calibration models for the determination of the ethanol content of ethanol fuel and alcoholic beverages. In the FT-NIR measurements the spectra were obtained using air as reference, and the spectral region for PLS modeling were selected based on the spectral distribution of the relative standard deviation in concentration. In the FT-Raman measurements hexachloro-1,3-butadiene (HCBD) has been used as an external standard. In the PLS/FT-NIR modeling for ethanol fuel analysis 50 ethanol fuel standards (84.9-100% (w/w)) were used (25 in the calibration, 25 in the validation). In the PLS/FT-Raman modeling 25 standards were used (13 in the calibration, 12 in the validation). The PLS/FT-NIR and FT-Raman models for beverage analysis made use of 24 standards (0-100% (v/v)). Twelve of them contained sugars (1-5% (w/w)), one-half was used in the calibration and the other half in the validation. Different spectral pre-processing were used in the PLS modeling, depending on the type of sample investigated. In the ethanol fuel analysis the FT-NIR pre-processing was a 17 points smoothed first derivative and for beverages no spectral pre-processing was used. The FT-Raman spectra were pre-processed by vector normalization in the ethanol fuel analysis and by a second derivative (17 points smoothing) in the beverage analysis. The PLS models were used in the analysis of real ethanol fuel and beverage samples. A t-test has shown that the FT-NIR model has an accuracy equivalent to that of the reference method (ASTM D4052) in the analysis of ethanol fuel, while in the analysis of beverages, the FT-Raman model presents an accuracy equivalent to the reference method. The limits of detection for NIR and Raman calibration models were 0.05 and 0.2% (w/w), respectively. It has also been shown that both techniques, present better results than gas chromatography (GC) in evaluating the ethanol content of beverages.     

Sustained prediction ability of net analyte preprocessing methods using reduced calibration sets. Theoretical and experimental study involving the spectrophotometric analysis of multicomponent mixtures

A newly developed multivariate method involving net analyte preprocessing (NAP) was tested using central composite calibration designs of progressively decreasing size regarding the multivariate simultaneous spectrophotometric determination of three active components (phenylephrine, diphenhydramine and naphazoline) and one excipient (methylparaben) in nasal solutions. Its performance was evaluated and compared with that of partial least-squares (PLS-1). Minimisation of the calibration predicted error sum of squares (PRESS) as a function of a moving spectral window helped to select appropriate working spectral ranges for both methods. The comparison of NAP and PLS results was carried out using two tests: (1) the elliptical joint confidence region for the slope and intercept of a predicted versus actual concentrations plot for a large validation set of samples and (2) the D-optimality criterion concerning the information content of the calibration data matrix. Extensive simulations and experimental validation showed that, unlike PLS, the NAP method is able to furnish highly satisfactory results when the calibration set is reduced from a full four-component central composite to a fractional central composite, as expected from the modelling requirements of net analyte based methods.     

Quantification of paracetamol in intact tablets using near-infrared transmittance spectroscopy.

Production batch samples of paracetamol tablets and specially prepared out-of-specification batches covering the range 90-110% of the stated amount (500 mg) were analysed by the BP official UV assay and by NIR transmittance spectroscopy. NIR measurements were made on 20 intact tablets from each batch, scanned five times each (10 min measurement time per batch) over the spectral range 6000-11,520 cm-1. An average spectrum was calculated for each batch. Partial least squares (PLS) regression models were set up using a calibration set (20 batches) between the NIR response and the reference tablet paracetamol content (UV). Various pre-treatments of the spectra were examined; the smallest relative standard error of prediction (0.73%) was obtained using the first derivative of the absorbance over the full spectrum. Only two principal components were required for the PLS model to give a good relationship between the spectral information and paracetamol content. Applying this model to the validation set (15 batches) gave a mean bias of -0.08% and a mean accuracy of 0.59% with relative standard deviations of 0.75 and 0.44%, respectively. The proposed method is non-destructive and therefore lends itself to on-line/at-line production control purposes. The method is easy to use and does not require a knowledge of the mass of the tablets.     

Spectrophotometric determination of ternary mixtures of thiamin, riboflavin and pyridoxal in pharmaceutical and human plasma by least-squares support vector machines.

Ternary mixtures of thiamin, riboflavin and pyridoxal have been simultaneously determined in synthetic and real samples by applications of spectrophotometric and least-squares support vector machines. The calibration graphs were linear in the ranges of 1.0 - 20.0, 1.0 - 10.0 and 1.0 - 20.0 microg ml(-1) with detection limits of 0.6, 0.5 and 0.7 microg ml(-1) for thiamin, riboflavin and pyridoxal, respectively. The experimental calibration matrix was designed with 21 mixtures of these chemicals. The concentrations were varied between calibration graph concentrations of vitamins. The simultaneous determination of these vitamin mixtures by using spectrophotometric methods is a difficult problem, due to spectral interferences. The partial least squares (PLS) modeling and least-squares support vector machines were used for the multivariate calibration of the spectrophotometric data. An excellent model was built using LS-SVM, with low prediction errors and superior performance in relation to PLS. The root mean square errors of prediction (RMSEP) for thiamin, riboflavin and pyridoxal with PLS and LS-SVM were 0.6926, 0.3755, 0.4322 and 0.0421, 0.0318, 0.0457, respectively. The proposed method was satisfactorily applied to the rapid simultaneous determination of thiamin, riboflavin and pyridoxal in commercial pharmaceutical preparations and human plasma samples.     

Improving the robustness of a partial least squares (PLS) model based on pure component selectivity analysis and range optimization: case study for the analysis of an etching solution containing hydrogen peroxide

Pure component selectivity analysis (PCSA) was successfully utilized to enhance the robustness of a partial least squares (PLS) model by examining the selectivity of a given component to other components. The samples used in this study were composed of NH₄OH, H₂O₂ and H₂O, a popular etchant solution in the electronic industry. Corresponding near-infrared (NIR) spectra (9000-7500 cm⁻¹) were used to build PLS models. The selective determination of H₂O₂ without influences from NH₄OH and H₂O was a key issue since its molecular structure is similar to that of H₂O and NH₄OH also has a hydroxyl functional group. The best spectral ranges for the determination of NH₄OH and H₂O₂ were found with the use of moving window PLS (MW-PLS) and corresponding selectivity was examined by pure component selectivity analysis. The PLS calibration for NH₄OH was free from interferences from the other components due to the presence of its unique NH absorption bands. Since the spectral variation from H₂O₂ was broadly overlapping and much less distinct than that from NH₄OH, the selectivity and prediction performance for the H₂O₂ calibration were sensitively varied depending on the spectral ranges and number of factors used. PCSA, based on the comparison between regression vectors from PLS and the net analyte signal (NAS), was an effective method to prevent over-fitting of the H₂O₂ calibration. A robust H₂O₂ calibration model with minimal interferences from other components was developed. PCSA should be included as a standard method in PLS calibrations where prediction error only is the usual measure of performance.     

Simultaneous determination of theophylline and caffeine by proton magnetic resonance spectroscopy using partial least squares regression techniques.

A 1H-NMR procedure based on an analysis of its data by a multivariate calibration method was conducted for the simultaneous determination of theophylline and caffeine in synthetic and real samples. Partial least squares regression (PLS) was chosen as the calibration method. The methyl signals of theophilline at 3.36 and 3.54 ppm that overlapped with those of caffeine were significant characteristics which were employed in this study for their analyses. The proposed method was successfully applied to recovery studies of theophylline and caffeine from real tablet samples.     

The prediction of total anthocyanin concentration in red-grape homogenates using visible-near-infrared spectroscopy and artificial neural networks

This study compares the performance of partial least squares (PLS) regression analysis and artificial neural networks (ANN) for the prediction of total anthocyanin concentration in red-grape homogenates from their visible-near-infrared (Vis-NIR) spectra. The PLS prediction of anthocyanin concentrations for new-season samples from Vis-NIR spectra was characterised by regression non-linearity and prediction bias. In practice, this usually requires the inclusion of some samples from the new vintage to improve the prediction. The use of WinISI LOCAL partly alleviated these problems but still resulted in increased error at high and low extremes of the anthocyanin concentration range. Artificial neural networks regression was investigated as an alternative method to PLS, due to the inherent advantages of ANN for modelling non-linear systems. The method proposed here combines the advantages of the data reduction capabilities of PLS regression with the non-linear modelling capabilities of ANN. With the use of PLS scores as inputs for ANN regression, the model was shown to be quicker and easier to train than using raw full-spectrum data. The ANN calibration for prediction of new vintage grape data, using PLS scores as inputs, was more linear and accurate than global and LOCAL PLS models and appears to reduce the need for refreshing the calibration with new-season samples. ANN with PLS scores required fewer inputs and was less prone to overfitting than using PCA scores. A variation of the ANN method, using carefully selected spectral frequencies as inputs, resulted in prediction accuracy comparable to those using PLS scores but, as for PCA inputs, was also prone to overfitting with redundant wavelengths.     

Metabolic trajectory characterisation of xenobiotic-induced hepatotoxic lesions using statistical batch processing of NMR data

Multivariate statistical batch processing (BP) analysis of 1H NMR urine spectra was employed to establish time-dependent metabolic variations in animals treated with the model hepatotoxin, alpha-naphthylisothiocyanate (ANIT). ANIT (100 mg kg(-1)) was administered orally to rats (n = 5) and urine samples were collected from dosed and matching control rats at time-points up to 168 h post-dose. Urine samples were measured via 1H NMR spectroscopy and partial least squares (PLS) based batch processing analysis was used to interpret the spectral data, treating each rat as an individual batch comprising a series of timed urine samples. A model defining the mean urine profile over the 7 day study period was established, together with model confidence limits (+/-3 standard deviation), for the control group. Samples obtained from ANIT treated animals were evaluated using the control model. Time-dependent deviations from the control model were evident in all ANIT treated animals consisting of glycosuria, bile aciduria, an initial decrease in taurine levels followed by taurinuria and a reduction of tricarboxylic acid cycle intermediate excretion. BP provided an efficient means of visualising the biochemical response to ANIT in terms of both inter-animal variation and net variation in metabolite excretion profiles. BP also allowed multivariate statistical limits for normality to be established and provided a template for defining the sequence of time-dependent metabolic consequences of toxicity in NMR based metabonomic studies.