Resolving batch chromatographic overlapping peaks of flavoring essence using stepwise key spectrum selection

Stepwise key spectrum selection (SKSS) was introduced to resolve batch overlapping peaks from gas chromatography-mass spectrometry (GC-MS) analysis of ten batch tobacco flavoring samples in different storage times. Resolution was implemented on a software platform that embedded the SKSS method. The data from GC-MS analysis of the samples were saved and prepared in ASCII files and then were inputted into the software platform for visual inspections. The data segment with overlapping peaks was precut for subsequent analysis. Spectral background in the data was removed using a linear fitting of the baseline. Four components in the overlapping peaks were automatically detected by the SKSS method. The resolution of the concentration profiles and spectra of the four components was conducted by setting only one parameter, the negative area ratio, as 0.01. The fixed size moving window evolving factor analysis and evolving factor analysis were applied to validate the resolved concentration profiles. The resolved mass spectra were validated by the searched standard through library search at the pure component regions revealed by the resolved concentration profiles. The results showed that the SKSS method could be a simple but powerful tool in resolving batch chromatographic overlapping peaks.     

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.     

Problem of mixtures with known compositions and IRONFLEA method for multivariate curve resolution

A special case of gray spectral data systems [(a) F.-T. Chau, Y.-Z. Liang, J. Gao, X.-G. Shao (Eds.), Chemometrics: From Basics to Wavelet Transform, Chemical Analysis Series, vol. 164, John Wiley & Sons, Inc., 2004; (b) Y.Z. Liang, O.M. Kvalheim, R. Manne, Chemom. Intell. Lab. Syst. 18 (1993) 235-250] is discussed here and the least-squares method for the multivariate curve resolution (MCR) named IRONFLEA is proposed. The system under consideration is the bilinear spectral data of the samples with known chemical compositions and unknown concentration matrix. If the spectra of samples (A_i) and (Q + A_i) (i = 1, …, n, n ≥ 2) are available, then the spectrum and the concentrations of Q could be found and the solution is unique. A practical chemical model for this problem could be mixtures, polymers, peptides, oligosaccharides, or supramolecular formations made of a limited number of monomeric components. In the cases of polymeric or oligomeric samples the spectral contributions and the concentrations of the particular monomeric units are extracted. The method is capable of extracting chemically meaningful spectra of components. The method is implemented in SAS IML code and tested for the deconvolution of spectra of polymers made of styrene derivatives with known monomeric compositions [(a) H. Fenniri, L. Ding, A.E. Ribbe, Y. Zyrianov, J. Am. Chem. Soc. 123 (2001) 8151-8152; (b) H. Fenniri, S. Chun, L. Ding, Y. Zyrianov, K. Hallenga, J. Am. Chem. Soc. 125 (2003) 10546-10560]. The method performs calculations fast enough to allow the incorporation of leave-one-out outlier removal procedure.     

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.     

Detection and use of selective ions for providing initial estimates of the concentration vectors in iterative target transformation factor analysis

The characteristic of discreteness in mass spectral direction distinguishes GC/MS from other spectral techniques. Based on this feature, we propose a new method to construct the initial concentration vectors for iterative target transformation factor analysis (ITTFA). For each chemical component, a search for a selective ion with good signal-to-noise ratio is first conducted using evolving factor analysis (EFA) and information entropy. The corresponding chromatogram of the selective ion is subsequently applied to construct an initial concentration vector for ITTFA. Special strategies are developed to cope with chromatographic patterns with embedded peaks and complex multicomponent structure. Results from three simulated and one real mixture and comparison with results from heuristic evolving latent projections (HELP) and a previously published method for definition of the initial profiles for ITTFA, indicate that selective ion chromatogram (SIC) ITTFA represents a fast, automatic and accurate method for resolution of GC/MS data.     

Soft modelling for the resolution of highly overlapped voltammetric peaks: application to some Pb-phytochelatin systems

A differential pulse polarographic (DPP) study of the Pb²⁺/Cys-Gly, Pb²⁺/γ-Glu-Cys, Pb²⁺/PC₂ and Pb²⁺/PC₃ systems is performed, being PC₂ and PC₃ the phytochelatins of general structure (γ-Glu-Cys)_n-Gly, with n = 2 and 3, respectively. Analysis of DPP data is assisted by multivariate curve resolution with alternating least squares (MCR-ALS) method in order to establish the complexes formation sequence and their final stoichiometries. DPP signals of these systems present, besides overlapping of peaks due to free metal ion and metal complexes, interference of mercury anodic signals. Despite these complications, MCR-ALS allows us to propose a model of complexation for each system, and some tentative structures for the complexes.     

Resolution of overlapping capillary electrophoresis peaks by using chemometric analysis: quantification of the components in compound reserpine tablets

The application of multivariate curve resolution with alternating least squares (MCR-ALS) methods to second-order data from capillary electrophoresis with diode array detector (CE-DAD) is reported. Initial qualitative solutions obtained by evolving factor analysis (EFA) and pure-variable detection method can be further optimized by a simultaneous analysis of multiple electrophoresis run data with ALS regression. While unknown samples are analyzed simultaneously against the corresponding standards in different composition ratios, the exact amounts of common components in different CE runs can be determined by the traditional calibration curve method, and quantification can thus be achieved. The above methods are applied to the determination of the components in compound reserpine tablets in overlapping peaks from CE. The quantification results are compared with those of the first derivative of the electropherogram method and artificial neural network (ANN) method.     

Rapid method for drawing overlapping resolution maps for thin-layer chromatography

Summary The plotting of an overlapping resolution map is a very long process. Derivation of an original resolution formula provides an equation that gives directly the solvent-front migration-distance required for the separation of various compounds. This equation is used to speed up 30-fold the drawing of overlapping maps.     

Second-order advantage from kinetic-spectroscopic data matrices in the presence of extreme spectral overlapping A multivariate curve resolution--alternating least-squares approach

Multivariate curve resolution coupled to alternating least-squares (MCR-ALS) has been employed to model kinetic-spectroscopic second-order data, with focus on the achievement of the important second-order advantage, under conditions of extreme spectral overlapping among sample components. A series of simulated examples shows that MCR-ALS can conveniently handle the studied analytical problem unlike other second-order multivariate calibration algorithms, provided matrix augmentation is implemented in the spectral mode instead of in the usual kinetic mode. The approach has also been applied to three experimental examples, which involve the determination of: (1) the antiparkinsonian carbidopa (analyte) in the presence of levodopa as a potential interferent, both reacting with cerium (IV) to produce the fluorescent species cerium (III) with different kinetics; (2) Fe(II) (analyte) in the presence of the interferent Zn(II), both catalyzing the oxidation of methyl orange with potassium bromate; and (3) tartrazine (analyte) in the presence of the interferent brilliant blue, both oxidized with potassium bromate, with the interferent leading to a product with an absorption spectrum very similar to tartrazine. The results indicate good analytical performance towards the analytes, despite the intense spectral overlapping and the presence of unexpected constituents in the test samples.     

Combined Kalman filtering and steepest descent minimization for quantitative analysis of unresolved misaligned chromatographic peaks. The resolution of alternariol and altenuisol mycotoxins in mixtures by HPLC

Summary The potential of a computational approach for the quantitative resolution of seriously overlapping chromatographic peaks when there is loss of collinearity between the pure component peaks and the mixture peak has been explored. The program makes iterative use of the Kalman filter algorithm for resolving the mixture peak with the component peaks aligned according to some values of the position parameters, and of a steepest descent minimization procedure to find the optimal alignment. This combined procedure has been applied to the quantitive resolution of the HPLC chromatograms of alternariol and altenuisol mycotoxins in synthetic mixtures and in real samples.     

Spline wavelet in the resolution of overlapping voltammetric peaks

Spline wavelet transform is used to resolve overlapping voltammetric peaks. A suitable resolving factor is chosen to multiply the filters of spline wavelet and make it a peak resoluter. Simulated overlapping voltammetric peaks are processed by the peak resoluter and satisfactory results are obtained. Base-line separation can be achieved, the relative errors of peak position are less than 3.0%, and the relative errors of peak area are less than 5.0%. The effect of different resolving factors and spline wavelet basis are discussed. To test the procedure, two systems, cadmium (Ⅱ)-indium (Ⅲ) and lead (Ⅱ)-thallium (Ⅰ), are used.     

Towards unsupervised analysis of second-order chromatographic data: automated selection of number of components in multivariate curve-resolution methods

A method to apply multivariate curve-resolution unattendedly is presented. The algorithm is suitable to perform deconvolution of two-way data (e.g. retrieving the individual elution profiles and spectra of co-eluting compounds from signals obtained from a chromatograph equipped with multiple-channel detection: LC-DAD or GC-MS). The method is especially adequate to achieve the advantages of deconvolution approaches when huge amounts of data are present and manual application of multivariate techniques is too time-consuming. The philosophy of the algorithm is to mimic the reactions of an expert user when applying the orthogonal projection approach--multivariate curve-resolution techniques. Basically, the method establishes a way to check the number of significant components in the data matrix. The performance of the method was superior to the Malinowski F-test. The algorithm was tested with HPLC-DAD signals.     

Optimization of supercritical fluid chromatographic separation of nitroaromatics by the overlapping resolution mapping scheme

The use of an interpretative optimization method to predict optimum conditions for the supercritical fluid chromatographic separation of nitroaromatics is demonstrated. The scheme utilizes the overlapping resolution mapping (ORM) procedure for multivariate prediction of the pressure gradient and temperature in order to separate eight nitroaromatics. Using the ORM scheme, only seven experiments were conducted before a global optimum was revealed. All the peaks were satisfactorily separated using experimental conditions derived from the scheme.     

Resolution method for overlapping peaks based on the fractional-order differential

Equations between the differential order and the maximum of the fractional-order differential for the specified peak signals are developed based on the variation of the maximum of the specified peak signals at different orders. Also, equations between the differential order and the zero-crossing of the fractional-order differential of the specified peak signals are proposed according to the variation of the zero-crossing of the specified peak signals at different orders. Characteristic paramters of the Gaus- sian peak, Lorentzian peak, and Tsallis peak can be estimated using estimator I and estimator II which are obtained by the equations above. As a result, a new method is presented to resolve the overlapped peaks signal. Firstly, a fractional-order differential of the specified peak signals is obtained with the fractional-order differentiation filter. Then, characteristic paramters of the specified peak signals can be extracted using estimator I and estimator II. Finally, the Tsallis peak is used as a model to assign the overlapping peak signals correctly. Experimental results show that the proposed method is efficient and effective for the simulated overlapping peaks and detected overlapping voltammetric peak signals.     

A new statistical method for the automated detection of peaks in UV-DAD chromatograms of a sample mixture

One of the major issues within the context of the fully automated development of chromatographic methods consists of the automated detection and identification of peaks coming from complex samples such as multi-component pharmaceutical formulations or stability studies of these formulations. The same problem can also occur with plant materials or biological matrices. This step is thus critical and time-consuming, especially when a Design of Experiments (DOE) approach is used to generate chromatograms. The use of DOE will often maximize the changes of the analytical conditions in order to explore an experimental domain. Unfortunately, this generally provides very different and “unpredictable” chromatograms which can be difficult to interpret, thus complicating peak detection and peak tracking (i.e. matching peaks among all the chromatograms). In this context, Independent Components Analysis (ICA), a new statistically based signal processing methods was investigated to solve this problem. The ICA principle assumes that the observed signal is the resultant of several phenomena (known as sources) and that all these sources are statistically independent. Under those assumptions, ICA is able to recover the sources which will have a high probability of representing the constitutive components of a chromatogram. In the present study, ICA was successfully applied for the first time to HPLC–UV-DAD chromatograms and it was shown that ICA allows differentiation of noise and artifact components from those of interest by applying clustering methods based on high-order statistics computed on these components. Furthermore, on the basis of the described numerical strategy, it was also possible to reconstruct a cleaned chromatogram with minimum influence of noise and baseline artifacts. This can present a significant advance towards the objective of providing helpful tools for the automated development of liquid chromatography (LC) methods. It seems that analytical investigations could be shortened when using this type of methodologies.     

Multivariate curve resolution of incomplete fused multiset data from chromatographic and spectrophotometric analyses for drug photostability studies

An advanced and powerful chemometric approach is proposed for the analysis of incomplete multiset data obtained by fusion of hyphenated liquid chromatographic DAD/MS data with UV spectrophotometric data from acid–base titration and kinetic degradation experiments. Column- and row-wise augmented data blocks were combined and simultaneously processed by means of a new version of the multivariate curve resolution-alternating least squares (MCR-ALS) technique, including the simultaneous analysis of incomplete multiset data from different instrumental techniques. The proposed procedure was applied to the detailed study of the kinetic photodegradation process of the amiloride (AML) drug. All chemical species involved in the degradation and equilibrium reactions were resolved and the pH dependent kinetic pathway described.     

Fast chromatographic method for the determination of dyes in beverages by using high performance liquid chromatography—Diode array detection data and second order algorithms

A fast chromatographic methodology is presented for the analysis of three synthetic dyes in non-alcoholic beverages: amaranth (E123), sunset yellow FCF (E110) and tartrazine (E102). Seven soft drinks (purchased from a local supermarket) were homogenized, filtered and injected into the chromatographic system. Second order data were obtained by a rapid LC separation and DAD detection. A comparative study of the performance of two second order algorithms (MCR-ALS and U-PLS/RBL) applied to model the data, is presented. Interestingly, the data present time shift between different chromatograms and cannot be conveniently corrected to determine the above-mentioned dyes in beverage samples. This fact originates the lack of trilinearity that cannot be conveniently pre-processed and can hardly be modelled by using U-PLS/RBL algorithm. On the contrary, MCR-ALS has shown to be an excellent tool for modelling this kind of data allowing to reach acceptable figures of merit. Recovery values ranged between 97% and 105% when analyzing artificial and real samples were indicative of the good performance of the method. In contrast with the complete separation, which consumes 10 mL of methanol and 3 mL of 0.08 mol L⁻¹ ammonium acetate, the proposed fast chromatography method requires only 0.46 mL of methanol and 1.54 mL of 0.08 mol L⁻¹ ammonium acetate. Consequently, analysis time could be reduced up to 14.2% of the necessary time to perform the complete separation allowing saving both solvents and time, which are related to a reduction of both the costs per analysis and environmental impact.     

Investigation of the equality constraint effect on the reduction of the rotational ambiguity in three-component system using a novel grid search method

The obtained results by soft modeling multivariate curve resolution methods often are not unique and are questionable because of rotational ambiguity. It means a range of feasible solutions equally fit experimental data and fulfill the constraints. Regarding to chemometric literature, a survey of useful constraints for the reduction of the rotational ambiguity is a big challenge for chemometrician. It is worth to study the effects of applying constraints on the reduction of rotational ambiguity, since it can help us to choose the useful constraints in order to impose in multivariate curve resolution methods for analyzing data sets. In this work, we have investigated the effect of equality constraint on decreasing of the rotational ambiguity. For calculation of all feasible solutions corresponding with known spectrum, a novel systematic grid search method based on Species-based Particle Swarm Optimization is proposed in a three-component system.     

Correction method for quantitative area determination of overlapping chromatographic peaks based on the exponentially modified Gaussian (EMG) model

A simple method is presented for peak area correction of overlapping peaks. This correction is necessary for the normal approach of dealing with overlapping peaks by a vertical line at the valley point. The relative area errors caused by this vertical line are calculated as the correction factors in three dimensions of peak separation, peak ratio, and peak tailing skew. The calculation is based on the exponentially modified Gaussian asymmetric peak model.     

Parametric signal fitting by gaussian peak adjustment: A new multivariate curve resolution method for non-bilinear voltammetric measurements

A new methodology based on the fitting of signals to parametric functions is proposed for the multivariate curve resolution (MCR) analysis of overlapping and peak-shaped voltammetric signals which progressively get broader or narrower and move along the potential axis, thus causing a dramatic loss of linearity. The method is based on the least squares fitting of gaussian functions at both sides of the peaks by using adjustable parameters for the peak height, position and symmetry. It consists of several home-made programs written in Matlab environment, which are freely available as supplementary material of the present work. The application to the systems Zn(II)–oxalate, and to the phytochelatin PC₅ in a wide pH range provides excellent results as compared to these of more conventional linear methods, which raises good expectations about future application to electrochemical and even non-electrochemical data.