Chemometric processing of second-order liquid chromatographic data with UV–vis and fluorescence detection. A comparison of multivariate curve resolution and parallel factor analysis 2

Second-order liquid chromatographic data with multivariate spectral (UV–vis or fluorescence) detection usually show changes in elution time profiles from sample to sample, causing a loss of trilinearity in the data. In order to analyze them with an appropriate model, the latter should permit a given component to have different time profiles in different samples. Two popular models in this regard are multivariate curve resolution-alternating least-squares (MCR-ALS) and parallel factor analysis 2 (PARAFAC2). The conditions to be fulfilled for successful application of the latter model are discussed on the basis of simple chromatographic concepts. An exhaustive analysis of the multivariate calibration models is carried out, employing both simulated and experimental chromatographic data sets. The latter involve the quantitation of benzimidazolic and carbamate pesticides in fruit and juice samples using liquid chromatography with diode array detection, and of polycyclic aromatic hydrocarbons in water samples, in both cases in the presence of potential interferents using liquid chromatography with fluorescence spectral detection, thereby achieving the second-order advantage. The overall results seem to favor MCR-ALS over PARAFAC2, especially in the presence of potential interferents.     

Quantification from highly drifted and overlapped chromatographic peaks using second-order calibration methods

For determining low levels of pesticides and phenolic compounds in river and wastewater samples by high performance liquid chromatography (HPLC), solid phase extraction (SPE) is commonly used before the chromatographic separation. This preconcentration step is not necessarily selective for the analytes of interest and it may retain other compounds of similar characteristics as well. In this case, we present, humic and fulvic acids caused a large baseline drift and overlapped the analytes to be quantified. The inaccurate determinations of the area of the peaks of these analytes made it difficult to quantify them with univariate calibration. Here we compare three second-order calibration algorithms (generalized rank annihilation method (GRAM), parallel factor analysis (PARAFAC) and multivariate curve resolution-alternating least squares (MCR-ALS)) which efficiently solve this problem. These methods use second-order data, i.e., a matrix of responses for each peak, which is easily obtained with a high performance liquid chromatography-diode array detector (HPLC-DAD). With these methods, the area does not need to be directly measured and predictions are more accurate. They also save time and resources because they can quantify analytes even if the peaks are not resolved. GRAM and PARAFAC require trilinear data. Biased and imprecise concentrations (relative standard deviation, %R.S.D. = 34) were obtained without correcting the time-shift. Hence, a time-shift correction algorithm to align the peaks was needed to obtain accurate predictions. MCR-ALS was the most robust to the time-shift. All three algorithms provided similar mean predictions, which were comparable to those obtained when sulfite was added to the samples. However, the predictions for the different replicates were more similar for the second-order algorithms (%R.S.D. = 3) than the ones obtained by univariate calibration after the sulfite addition (%R.S.D. = 13).     

Modeling four and three-way fast high-performance liquid chromatography with fluorescence detection data for quantitation of fluoroquinolones in water samples

This paper presents a study regarding the acquisition and analytical utilization of four and three-way data, acquired by following the excitation–emission fluorescence matrices at different elution times, in a fast liquid chromatographic HPLC procedure. This kind of data were implemented for first time for quantitative purposes, and applied to the determination of two fluoroquinolones in tap water samples, as a model to show the potentiality of the proposed strategy of four-way data generation. The data were modeled with three well-known algorithms: PARAFAC, U-PLS/RTL and MCR-ALS, the latter conveniently adapted to model third-order data. The second-order advantage was exploited when analyzing samples containing uncalibrated interferences. PARAFAC and MCR-ALS were the algorithms that better exploited the second-order advantage when no peak time shifts occurred among samples. On the other hand, when the quadrilinearity was lost due to the occurrence of temporal shifts, MCR-ALS furnished the better results. Relative error of prediction (REP%) obtained were 9.9% for ofloxacin and 14.0% for ciprofloxacin. In addition, a significant enhancement in the analytical figures of merit was observed when going from second- to third-order data (reduction of ca. 70% in LODs).     

Novel combination of non-aqueous capillary electrophoresis and multivariate curve resolution-alternating least squares to determine phenolic acids in virgin olive oil

This paper presents the development of a non-aqueous capillary electrophoresis method coupled to UV detection combined with multivariate curve resolution-alternating least-squares (MCR-ALS) to carry out the resolution and quantitation of a mixture of six phenolic acids in virgin olive oil samples. p-Coumaric, caffeic, ferulic, 3,4-dihydroxyphenylacetic, vanillic and 4-hydroxyphenilacetic acids have been the analytes under study. All of them present different absorption spectra and overlapped time profiles with the olive oil matrix interferences and between them. The modeling strategy involves the building of a single MCR-ALS model composed of matrices augmented in the temporal mode, namely spectra remain invariant while time profiles may change from sample to sample. So MCR-ALS was used to cope with the coeluting interferences, on accounting the second order advantage inherent to this algorithm which, in addition, is able to handle data sets deviating from trilinearity, like the data herein analyzed. The method was firstly applied to resolve standard mixtures of the analytes randomly prepared in 1-propanol and, secondly, in real virgin olive oil samples, getting recovery values near to 100% in all cases. The importance and novelty of this methodology relies on the combination of non-aqueous capillary electrophoresis second-order data and MCR-ALS algorithm which allows performing the resolution of these compounds simplifying the previous sample pretreatment stages.     

Multivariate selectivity as a metric for evaluating comprehensive two-dimensional gas chromatography-time-of-flight mass spectrometry subjected to chemometric peak deconvolution

Two-dimensional gas chromatography (GC x GC) coupled to time-of-flight mass spectrometry (TOFMS) [GC x GC-TOFMS)] is a highly selective technique well suited to analyzing complex mixtures. The data generated is information-rich, making it applicable to multivariate quantitative analysis and pattern recognition. One separation on a GC x GC-TOFMS provides retention times on two chromatographic columns and a complete mass spectrum for each component within the mixture. In this report, we demonstrate how GC x GC-TOFMS combined with trilinear chemometric techniques, specifically parallel factor analysis (PARAFAC) initiated by trilinear decomposition (TLD), results in a powerful analytical methodology for multivariate deconvolution. Using PARAFAC, partially resolved components in complex mixtures can be deconvoluted and identified without requiring a standard data set, signal shape assumptions or any fully selective mass signals. A set of four isomers (iso-butyl, sec-butyl, tert-butyl, and n-butyl benzenes) is used to investigate the practical limitations of PARAFAC for the deconvolution of isomers at varying degrees of chromatographic resolution and mass spectral selectivity. In this report, multivariate selectivity was tested as a metric for evaluating GC x GC-TOFMS data that is subjected to PARAFAC peak deconvolution. It was found that deconvolution results were best with multivariate selectivities over 0.18. Furthermore, the application of GC x GC-TOFMS followed by TLD/PARAFAC is demonstrated for a plant metabolite sample. A region of GC x GC-TOFMS data from a complex natural sample of a derivatized metabolic plant extract from Huilmo (Sisyrinchium striatum) was analyzed using TLD/PARAFAC, demonstrating the utility of this analytical technique on a natural sample containing overlapped analytes without selective ions or peak shape assumptions.     

Multivariate curve resolution-particle swarm optimization: A high-throughput approach to exploit pure information from multi-component hyphenated chromatographic signals

Multivariate curve resolution-particle swarm optimization (MCR-PSO) algorithm is proposed to exploit pure chromatographic and spectroscopic information from multi-component hyphenated chromatographic signals. This new MCR method is based on rotation of mathematically unique PCA solutions into the chemically meaningful MCR solutions. To obtain a proper rotation matrix, an objective function based on non-fulfillment of constraints is defined and is optimized using particle swarm optimization (PSO) algorithm. Initial values of rotation matrix are calculated using local rank analysis and heuristic evolving latent projection (HELP) method. The ability of MCR-PSO in resolving the chromatographic data is evaluated using simulated gas chromatography–mass spectrometry (GC–MS) and high-performance liquid chromatography–diode array detection (HPLC–DAD) data. To present a comprehensive study, different number of components and various levels of noise under proper constraints of non-negativity, unimodality and spectral normalization are considered. Calculation of the extent of rotational ambiguity in MCR solutions for different chromatographic systems using MCR-BANDS method showed that MCR-PSO solutions are always in the range of feasible solutions like true solutions. In addition, the performance of MCR-PSO is compared with other popular MCR methods of multivariate curve resolution-objective function minimization (MCR-FMIN) and multivariate curve resolution-alternating least squares (MCR-ALS). The results showed that MCR-PSO solutions are rather similar or better (in some cases) than other MCR methods in terms of statistical parameters. Finally MCR-PSO is successfully applied in the resolution of real GC–MS data. It should be pointed out that in addition to multivariate resolution of hyphenated chromatographic signals, MCR-PSO algorithm can be straightforwardly applied to other types of separation, spectroscopic and electrochemical data.     

Determination of pesticides and metabolites in wine by high performance liquid chromatography and second-order calibration methods

The models parallel factor analysis (PARAFAC) and the recently introduced bilinear least squares (BLLS) were applied to develop second-order calibration methods to high performance liquid chromatography with diode array detection (HPLC-DAD) data, where overlap of interferences with the compounds of interest was observed, making the determination and resolution of the analytes possible. In this work, the simultaneous determination of five pesticides and two metabolites in wine samples by HPLC-DAD was performed, using the second-order advantage. The results of two chromatographic methods were compared, involving either isocratic or gradient elution. An appropriate preprocessing method was necessary to correct the effects of time shifts, baseline variations and background. BLLS presented results that were of the same quality as PARAFAC in five cases, but in two other situations only PARAFAC enabled analyte quantitation. Relative errors of prediction lower than 10% for all compounds were obtained, indicating that the methodology employing HLPC-DAD and second-order calibration can handle complex analytical systems.     

Automated resolution of overlapping peaks in chromatographic data

Parallel factor analysis 2 (PARAFAC2) has been shown to be a powerful tool for resolution of complex overlapping peaks in chromatographic analyses. It is particularly useful because of its ability to handle shifts in the elution time mode and peak shape changes. Like all curve resolution techniques, PARAFAC2 will only find chemically meaningful parameters (elution time profiles and mass spectra) if the correct number of factors are determined. So far, the primary way to determine an appropriate number of factors, when using PARAFAC2, is to calculate models with different number of factors and then inspect the models manually. This approach is time consuming, and the result may be biased because of the manual assessment of the model quality, making PARAFAC2 inaccessible for analytical chemists in general. Here, we develop a method that can determine an appropriate number of factors in an automated way. The automation is based on a number of model diagnostics (quality criteria) collected from models with different numbers of factors. Combining these diagnostics, it is possible to assess what the appropriate number of components is. In this work, only gas chromatography–mass spectrometry data are considered. However, it will most likely be fairly straightforward to expand the work to also cover liquid chromatography data (with a multivariate detector). Automating the model quality evaluation of the PARAFAC2 model enables both the inexperienced and trained user to perform comprehensive and advanced analysis of chromatographic data with a minimum of manual work. © 2013 The Authors. Journal of Chemometrics Published by John Wiley & Sons Ltd.     

Photoinduced spectrofluorimetric determination of fluoroquinolones in human urine by using three- and two-way spectroscopic data and multivariate calibration

Multivariate methods comprise of a group of chemometric tools allowing the analysis of different analytical data, i.e., spectroscopic, chromatographic obtained from multichannel detector systems. Second-way data are widely used in analytical applications in combination with multivariate calibration methods, but three- and higher-way data are yet not as widely applied. In complex biological samples, the employment of the three-way data is of special interest, as they may be combined with methods that exploit the second-order advantage allowing calculating individual concentrations of the analytes of interest in the presence of unknown interferences in untreated samples. A very sensitive and selective method is proposed, by coupling photoinduced fluorescence and multivariate analysis of the three-way data excitation-emission fluorescence matrices (EEMs), of the photoproducts obtained from UV irradiation of three fluoroquinolones: enoxacin (ENO), norfloxacin (NOR) and ofloxacin (OFLO). The application of a previous photoirrradiation process allows the determination of mixtures of ENO, NOR and OFLO, in urine samples at biological levels without sample pretreatments. The resolution ability of N-way partial least squares (N-PLS), parallel factor analysis (PARAFAC) and self weighted alternating trilinear decomposition (SWATLD), is compared with partial least squares (PLS) and unfolded-PLS (U-PLS), in the analysis of ENO, NOR and OFLO in human urine samples.     

Validation of an analytical method to determine sulfamides in kidney by HPLC-DAD and PARAFAC2 with first-order derivative chromatograms

Six sulfamides were extracted from kidney and analysed by high-performance liquid chromatography with diode array detection (HPLC-DAD): sulfadiazine, sulfamethazine, sulfamethoxypyridazine, sulfamethoxazole, sulfadimethoxine and sulfaquinoxaline.Two main difficulties arose in identifying and quantifying the analytes. Firstly, the chromatographic peaks of the matrix interferences overlapped with those of the analytes. The uniqueness property of PARAFAC2 solved this problem. Secondly, the gradient elution caused a baseline drift. The first-order derivative of the chromatograms minimized its effect.The analytical method was validated. As the performance criteria detailed in the European Decision 2002/657/EC are based on specific signals, this paper generalizes those criteria for higher-order and non-specific signals. In this sense the proposed methodology is general and can be applied to any chromatographic method (HPLC or GC) with a detector that provide a multivariate signal (MS, DAD, EC, etc.).     

Simultaneous determination of tyrosine and dopamine in urine samples using excitation-emission matrix fluorescence coupled with second-order calibration

A "green" and quick analytical method for complex compounds was developed for simultaneous determination of tyrosine (Tyr) and dopamine (DA) in urine samples in this paper. The three-way responsive data recorded by excitation-emission matrix fluorescence (EEM) spectrometer was analyzed using second-order calibration methods based on both parallel factor analysis (PARAFAC) and selfweighted alternating trilinear decomposition (SWATLD) algorithms. The EEM spectra of the analytes were overlapped with the background in urine samples. However the second-order advantage of both PARAFAC and SWATLD methods was exploited, even in the presence of unknown interferences and the satisfactory results can be obtained. Furthermore, the linear ranges of Tyr and DA were determined to be 0.042-6.42 μg/mL and 0.18-4.43 μmg/mL, respectively, and the accuracies of both methods were validated by the analytical figures of merit (FOM).     

Multiway chemometric decomposition of EEM of fluorescence of CdTe quantum dots obtained as function of pH

The effect of the pH (from 3 to 10) on the excitation emission matrices (EEMs) of fluorescence of CdTe quantum dots (QDs), capped with mercaptopropionic acid (MPA), were analyzed by multiway decomposition methods of parallel factor analysis (PARAFAC), a variant of the parallel factor analysis method (PARAFAC2) and multivariate curve resolution alternating least squares (MCR-ALS). Three different sized CdTe QDs with emission maximum at 555 nm (QDa), 594 nm (QDb) and 628 nm (QDc) were selected for analysis. The three-way data structures composed of sets of EEMs obtained as function of the pH (EEMs, pH) do not have a trilinear structure. A marked deviation to the trilinearity is observed in the emission wavelength order—the emission spectra suffers wavelength shift as the pH is varied. The pH-induced variation of the fluorescence properties of QDs is described with only one-component PARAFAC2 or MCR-ALS models—other components are necessary to model scattering and/or other background signals in (EEMs, pH) data structures. Bigger sized QDs are more suitable tools for analytical methodologies because they show higher Stokes shifts (resulting in simpler models) and higher pH range sensitivity. The pH dependence of the maximum wavelength of the emission spectra is particularly suitable for the development of QDs/EEMs wavelength-encoded pH sensor bioimaging or biological label methodologies when coupled to multiway chemometric decomposition.     

Algorithm combination strategy to obtain the second‐order advantage: simultaneous determination of target analytes in plasma using three‐dimensional fluorescence spectroscopy

Although a number of algorithms have established to obtain the well‐known second‐order advantage that quantifies analytes of interest in the presence of interferents, each has associated problems. In this work, for the first time, the optimization procedure of trilinear decomposition has been divided into three subparts, and a novel strategy is developed for assembling the advantages of the alternating trilinear decomposition (ATLD) algorithm, the self‐weighted alternating trilinear decomposition (SWATLD) algorithm, and the parallel factor analysis (PARAFAC) algorithm. The performance of the proposed strategy was evaluated using a simulated data set, a published fluorescence data set together with a new fluorescence data set that simultaneously quantifies procaine and tetracaine in plasma. Results show that the novel method can accurately and effectively estimate the qualitative and quantitative information of analytes of interest. Besides, the resolved profiles are very stable with respect to the number of components as long as the employed number is chosen to be equal or larger than the underlying one. Additionally, the study confirms that better prediction can be obtained by the new strategy when compared with ATLD, SWATLD, and PARAFAC as well as the strategy that employs direct trilinear decomposition method as initial values for PARAFAC. Moreover, the strategy can be directly extended to third‐order or higher‐order data analysis. Copyright © 2012 John Wiley & Sons, Ltd.     

Ultrafast quantitation of six quinolones in water samples by second-order capillary electrophoresis data modeling with multivariate curve resolution–alternating least squares

This paper presents the development of a capillary electrophoresis method with diode array detector coupled to multivariate curve resolution–alternating least squares (MCR-ALS) to conduct the resolution and quantitation of a mixture of six quinolones in the presence of several unexpected components. Overlapping of time profiles between analytes and water matrix interferences were mathematically solved by data modeling with the well-known MCR-ALS algorithm. With the aim of overcoming the drawback originated by two compounds with similar spectra, a special strategy was implemented to model the complete electropherogram instead of dividing the data in the region as usually performed in previous works. The method was first applied to quantitate analytes in standard mixtures which were randomly prepared in ultrapure water. Then, tap water samples spiked with several interferences were analyzed. Recoveries between 76.7 and 125 % and limits of detection between 5 and 18 μg L^?1 were achieved.     

Matrix-free analysis of aflatoxins in pistachio nuts using parallel factor modeling of liquid chromatography diode-array detection data

In the present study a second-order calibration strategy for high performance liquid chromatography with diode-array detection (HPLC-DAD) has been developed using parallel factor analysis (PARAFAC) and has been applied for simultaneous determination of aflatoxins B₁, B₂, G₁ and G₂ in pistachio nuts in the presence of matrix interferences. Sample preparation was based on solvent extraction (SE) followed by solid phase extraction (SPE) on Bond Elut C18 cartridges. Since the sample preparation procedure was not selective to the analytes of interest, exploiting second-order advantage to obtain concentrations of individual analytes in the presence of uncalibrated interfering compounds seemed necessary. Appropriate pre-processing steps have been applied to correct background signals and the effect of retention time shifts. Transferred calibration data set obtained from standardization of solvent based calibration data has been used in prediction step. The results of PARAFAC on a set of spiked and naturally contaminated pistachio nuts indicated that the four aflatoxins could be successfully determined. The method was validated and multivariate analytical figures of merit were calculated. The advantages of the proposed method are using a low-cost SPE step relative to standard method of aflatoxin analysis (immune affinity column assay), a unique and simple isocratic elution program for all samples and a calibration transfer for saving both chemicals and time of analysis. This study show that coupling of SPE-HPLC-DAD with PARAFAC as a powerful second-order calibration method can be considered as an alternative method for resolution and quantification of aflatoxins in the presence of unknown interferences obtained through analysis of highly complex matrix of pistachio samples and cost per analysis can be reduced significantly.     

The application to wastewaters of chemometric approaches to handling problems of highly complex matrices

This overview covers current chemometric methodologies using second-order advantage to solve problems of analyzing highly complex matrices. Among the existing algorithms, it focuses on those most frequently used (e.g., the standard for second-order approaches to data analysis, PARAFAC (parallel factor analysis), and MCR-ALS (multivariate curve resolution alternating least squares), as well as the most recently implemented BLLS (bilinear least-squares), and U-PLS/RBL (unfolded partial least squares/residual bilinearization)). All of these are based on linear models. The overview also covers ANN/RBL (artificial neural networks followed by residual bilinearization), which achieves the second-order advantage in systems involving non-linear behavior. In addition, the overview deals with the drawbacks of these approaches, as well as other drawbacks that are inherent in the analytical techniques to question.     

A novel strategy on the spectrochromatographic analysis of a quaternary mixture by parallel factor analysis model

Poor chromatographic resolution is one of the main challenges in chromatographic analysis. Partially separated chromatographic peaks frequently occur, due to the nature of analytes and the demand for fast analysis using high flow rates and shorter columns. Modeling of chromatographic three-way data using suitable chemometric tools enables determining co-eluted peaks without using additional experimental efforts. In this paper, parallel factor analysis (PARAFAC) was applied to chromatographic data for the quantitative resolution of a quaternary mixture at the co-elution condition of acetaminophen, aspirin, ascorbic acid, and guaifenesin in a spectrochromatogram. The spectrochromatograms of the calibration set, validation set, and real samples were arranged as a three-way array. In the next step, the PARAFAC model was implemented to decompose the spectrochromatographic array into trilinear components, corresponding to spectral, chromatographic, and relative concentration profiles of the analytes. The chromatographic and spectral modes were used for the qualitative analysis of components, whereas the analytes in commercial tablets were quantified from their individual profiles in their concentration mode. This study indicated that the application of the PARAFAC model provided a novel strategy for determining overlapping peaks in a chromatogram to perform the analysis of multicomponent mixtures with reduced runtime and without additional efforts.     

Second-order multivariate calibration procedures applied to high-performance liquid chromatography coupled to fast-scanning fluorescence detection for the determination of fluoroquinolones

Different second-order multivariate calibration algorithms, namely parallel factor analysis (PARAFAC), N-dimensional partial least-squares (N-PLS) and multivariate curve resolution-alternating least-squares (MCR-ALS) have been compared for the analysis of four fluoroquinolones in aqueous solutions, including some human urine samples (additional four fluoroquinolones were simultaneously determined by univariate calibration). Data were measured in a short time with a chromatographic system operating in the isocratic mode. The detection system consisted of a fast-scanning spectrofluorimeter, which allows one to obtain second-order data matrices containing the fluorescence intensity as a function of retention time and emission wavelength. The developed approach enabled us to determine eight analytes, some of them with overlapped profiles, without the necessity of applying an elution gradient, and thus significantly reducing both the experimental time and complexity. The study was employed for the discussion of the scopes of the applied second-order chemometric tools. The quality of the proposed technique coupled to each of the evaluated algorithms was assessed on the basis of the figures of merit for the determination of fluoroquinolones in the analyzed water and urine samples. Univariate calibration of four analytes led to limits of detection in the range 20–40 ng mL⁻¹ and root mean square errors for the validation samples in the range 30–60 ng mL⁻¹ (corresponding to relative prediction errors of 3–8%). The ranges for second-order multivariate calibration (using PARAFAC and N-PLS) of the remaining four analytes were: limit of detection, 2–8 ng mL⁻¹, root mean square errors, 3–50 ng mL⁻¹ and relative prediction errors, 1–5%.     

Determination of phenolic compounds and authentication of PDO Lambrusco wines by HPLC-DAD and chemometric techniques

This work proposes a fast and simple method for detection and quantification of phenolic compounds in PDO Lambrusco wines using HPLC-DAD and chemometric techniques. Samples belonging to three different varieties of Lambrusco (Grasparossa, Salamino and Sorbara) were analyzed to provide a methodology appropriate for routine analysis. Given the high complexity of the sample and the coelution among chromatographic peaks, the use of chemometric techniques to extract the information of the individual eluting compounds was needed. Multivariate curve resolution-alternating least squares (MCR-ALS) allowed the resolution of the chromatographic peaks obtained and the use of this information for the quantification of the phenolic analytes in the presence of interferences. Use of multiset analysis and local rank/selectivity information was proven to be crucial for the correct resolution and quantification of compounds. The quantitative data provided by MCR-ALS about the phenolic targets and additional compounds present in the samples analyzed provided wine composition profiles, which were afterwards used to distinguish among wine varieties. Principal component analysis applied to the wine profiles allowed characterizing the wines according to their varieties.