Similarity analyses of chromatographic fingerprints as tools for identification and quality control of green tea.

Similarity assessment of complex chromatographic profiles of herbal medicinal products is important as a potential tool for their identification. Mathematical similarity parameters have the advantage to be more reliable than visual similarity evaluations of often subtle differences between the fingerprint profiles. In this paper, different similarity analysis (SA) parameters are applied on green-tea chromatographic fingerprint profiles in order to test their ability to identify (dis)similar tea samples. These parameters are either based on correlation or distance measurements. They are visualised in colour maps and evaluation plots. Correlation (r) and congruence (c) coefficients are shown to provide the same information about the similarity of samples. The standardised Euclidean distance (ds) reveals less information than the Euclidean distance (de), while Mahalanobis distances (dm) are unsuitable for the similarity assessment of chromatographic fingerprints. The adapted similarity score (ss*) combines the advantages of r (or c) and de. Similarity analysis based on correlation is useful if concentration differences between samples are not important, whereas SA based on distances also detects concentration differences well. The evaluation plots including statistical confidence limits for the plotted parameter are found suitable for the evaluation of new suspected samples during quality assurance. The ss* colour maps and evaluation plots are found to be the best tools (in comparison to the other studied parameters) for the distinction between deviating and genuine fingerprints. For all studied data sets it is confirmed that adequate data pre-treatment, such as aligning the chromatograms, prior to the similarity assessment, is essential. Furthermore, green-tea samples chromatographed on two dissimilar High-Performance Liquid Chromatography (HPLC) columns provided the same similarity assessment. Combining these complementary fingerprints did not improve the similarity analysis of the studied data set.

Recognizing paracetamol formulations with the same synthesis pathway based on their trace-enriched chromatographic impurity profiles.

The development of a new drug substance is an expensive and time-consuming process. Therefore, the developers want to maximize the profit from the drug by patenting the concerned molecule as well as its synthesis pathway. In a later stage a faster or cheaper manufacturing process can be developed and patented. The aim of this study is to recognize paracetamol-containing drug formulations in relation to their synthesis pathways, in order to demonstrate the possibility to reveal fraudulently synthesized paracetamol. Since different synthesis pathways require different starting materials, solvents, reagents and catalysts and since they can produce different intermediates, it is expected that drug products originating from a different synthesis pathway will exhibit a different impurity profile. Therefore, in this study several paracetamol samples, synthesized in four different ways, are analyzed using trace-enrichment high-performance liquid chromatography (HPLC). The resulting chromatographic data were chemometrically treated in order to reveal clustering tendencies in the hope of distinguishing the different pathways. When performing principal component analysis (PCA) only 3 vaguely outlined clusters appeared. Projection pursuit (PP) was able to reveal 4 clusters and the samples with known synthesis pathway, except one, were classified in the different clusters. When hierarchical clustering and auto-associative multivariate regression trees (AAMRT) were applied, the samples of the four synthesis pathways could also be distinguished. AAMRT has an added value, since it can indicate the variables (peaks and thus also the impurities) that are responsible for the differences between the samples synthesized differently.

Isotopic ratios to detect infringements of patents or proprietary processes of pharmaceuticals: two case studies.

Because of the increasing problem of drug counterfeiting and the potential danger related as well as the economic losses involved, the pharmaceutical industry and the regulatory instances are interested in the development of anti-counterfeiting and patent protection methodologies. In this paper, the evaluation of measured isotopic ratios by means of explorative chemometric techniques was performed to distinguish groups in two data sets containing samples of acetyl salicylic acid and ibuprofen, respectively. The samples in the data sets originated from different countries and manufacturers. For both compounds a clear distinction of groups of samples could be obtained. These groups could be explained based on the origin of the samples, both geographically as well as based on the manufacturer. Hypotheses were formulated concerning the synthetic pathways of the molecules and they were linked to the groups obtained with the chemometric tools.

Exploration of linear multivariate calibration techniques to predict the total antioxidant capacity of green tea from chromatographic fingerprints.

Nowadays fingerprinting is a generally applied technique for the identification and quality assessment of herbal products. In this study it was aimed to predict a quantitative property, the antioxidant capacity of green tea, from chromatographic fingerprints. Different linear multivariate calibration techniques, commonly applied on spectral data, were explored and compared. When the chromatograms were appropriately pretreated, all tested techniques were able to predict the total antioxidant capacity with a precision comparable to that of the reference method (Trolox equivalent antioxidant capacity assay). Stepwise multiple linear regression (MLR) however is less recommended because of inadequate variable selection. Principal components regression (PCR) also seems less preferable, because large variations not correlated with the total antioxidant capacity were also included in the model. This problem does not occur with partial least squares (PLS) models. Of all tested PLS methods, orthogonal projections to latent structures (O-PLS) was preferred because of its simplicity, reproducibility, good interpretability of the compounds' contribution to the antioxidant capacity and its good predictive and describing abilities.

Increasing conclusiveness of metabonomic studies by chem-informatic preprocessing of capillary electrophoretic data on urinary nucleoside profiles.

Nowadays, bioinformatics offers advanced tools and procedures of data mining aimed at finding consistent patterns or systematic relationships between variables. Numerous metabolites concentrations can readily be determined in a given biological system by high-throughput analytical methods. However, such row analytical data comprise noninformative components due to many disturbances normally occurring in analysis of biological samples. To eliminate those unwanted original analytical data components advanced chemometric data preprocessing methods might be of help. Here, such methods are applied to electrophoretic nucleoside profiles in urine samples of cancer patients and healthy volunteers. The electrophoretic nucleoside profiles were obtained under following conditions: 100 mM borate, 72.5 mM phosphate, 160 mM SDS, pH 6.7; 25 kV voltage, 30 degrees C temperature; untreated fused silica capillary 70 cm effective length, 50 microm I.D. Different most advanced preprocessing tools were applied for baseline correction, denoising and alignment of electrophoretic data. That approach was compared to standard procedure of electrophoretic peak integration. The best results of preprocessing were obtained after application of the so-called correlation optimized warping (COW) to align the data. The principal component analysis (PCA) of preprocessed data provides a clearly better consistency of the nucleoside electrophoretic profiles with health status of subjects than PCA of peak areas of original data (without preprocessing).

A comparison of three algorithms for chromatograms alignment.

In this paper the performance of three alignment algorithms, correlation optimized warping, parametric time warping and semi-parametric time warping, is compared on real chromatograms. Among these, parametric time warping is the simplest and fastest; generally less than 1s is required to align two chromatograms. It does not require the optimization of input parameters and allows the alignment of peak shifts in only one direction, or non-complex peak shifts in both directions. With correlation optimized warping and semi-parametric time warping complex peak shifts in both directions can be corrected but at the expense of the optimization of two input parameters. Semi-parametric time warping requires the selection of the proper number of B-splines in the warping function and, if necessary, the optimization of the penalty parameter. Often the default values can be used to obtain aligned signals. The optimization of the input parameters for correlation optimized warping (section length, slack) is not easy and time-consuming. Moreover, dependent on the input parameters, the computation time of the correlation optimized warping algorithm can be twice as long as for semi-parametric time warping for which computation times up to 23 s are required. However, the performance of both algorithms is equally good considering the improvement of the precision of the peak retention times and correlation coefficients between the chromatograms, after alignment. For the data aligned in this study, the average retention time precision and the lowest correlation before warping were 14 and 0.17, and were improved to three and 0.83, and six and 0.87 after warping, with correlation optimized warping and semi-parametric time warping, respectively.

Chemometric treatment of vanillin fingerprint chromatograms. Effect of different signal alignments on principal component analysis plots.

This study describes the chemometric treatment of vanillin fingerprint chromatograms to distinguish vanillin from different sources. Prior to principal component analysis, which is used to discriminate vanillin from different origins, the fingerprints are aligned. Three alignment algorithms are tested, correlation optimized warping (COW), target peak alignment (TPA) and semi-parametric time warping (STW). The performance of the three algorithms is evaluated and the effect of the different alignments on the PCA score plots is investigated. The alignment obtained with STW differs somewhat from that with COW and TPA. However, equivalent score plots were obtained regarding the different vanillin groups.

Automatic program for peak detection and deconvolution of multi-overlapped chromatographic signals part I: peak detection.

A series of two papers describing a procedure for automated peak deconvolution is presented. The goal is to develop a package of routines that can be used by non-experienced users. Part I (this paper) concerns peak detection, whereas Part II is dedicated to the deconvolution itself. In this first part, the most interesting features of the peak detection algorithms, which precede the deconvolution step, are outlined. High-order derivatives provide valuable information to assess the number of underlying compounds under a given peak cluster. A smoothing technique was found essential to compute properly the derivatives, since the noise is amplified when differences are calculated. The Savitsky-Golay smoother was applied in combination with the Durbin-Watson criterion to automate the window size selection. This strategy removed the noise without loosing valuable information. In some cases, it was found preferable to split the chromatogram in different elution regions, and apply the Durbin-Watson test and the Savitsky-Golay smoother to each region, separately. The derivatives allowed obtaining estimates of both peak parameters and the corresponding ranges for each eluting compound to be used in the deconvolution. An algorithm oriented to compare peaks from different chromatograms is also presented to perform deconvolution, using information from several related chromatograms.

Automatic program for peak detection and deconvolution of multi-overlapped chromatographic signals part II: peak model and deconvolution algorithms.

Several interlinked algorithms for peak deconvolution by non-linear regression are presented. These procedures, together with the peak detection methods outlined in Part I, have allowed the implementation of an automatic method able to process multi-overlapped signals, requiring little user interaction. A criterion based on the evaluation of the multivariate selectivity of the chromatographic signal is used to auto-select the most efficient deconvolution procedure for each chromatographic situation. In this way, non-optimal local solutions are avoided in cases of high overlap, and short computation times are obtained in situations of high resolution. A new algorithm, fitting both the original signal and the second derivatives is proved to avoid local optima in intermediate coelution situations. This allows achieving the global optimum without the need of background knowledge by the user. A previously reported peak model, a Gaussian with a polynomial standard deviation whose complexity can be modulated to enhance the fitting quality, was applied. However, the original formulation was modified to account baseline outside the peak region. Also, the optimal model complexity was auto-selected via error propagation theory. The method is able to process simultaneously several related chromatograms. The software was tested with both simulated and experimental chromatograms obtained with monolithic silica columns.

Prediction of total green tea antioxidant capacity from chromatograms by multivariate modeling.

In this paper, a fast strategy for determining the total antioxidant capacity of Chinese green tea extracts is developed. This strategy includes the use of experimental techniques, such as fast high-performance liquid chromatography (HPLC) on monolithic columns and a spectrophotometric approach to determine the total antioxidant capacity of the extracts. To extract the chemically relevant information from the obtained data, chemometrical approaches are used. Among them there are correlation optimized warping (COW) to align the chromatograms, robust principal component analysis (robust PCA) to detect outliers, and partial least squares (PLS) and uninformative variable elimination partial least squares (UVE-PLS) to construct a reliable multivariate regression model to predict the total antioxidant capacity from the fast chromatograms.

Development of a ginkgo biloba fingerprint chromatogram with UV and evaporative light scattering detection and optimization of the evaporative light scattering detector operating conditions.

A fingerprint chromatogram of a standardized Ginkgo biloba extract is developed on a monolithic silica column using a ternary gradient containing water, iso-propanol and tetrahydrofuran. For the detection, UV and evaporative light scattering (ELS) detectors are used, the latter allowing detection of the poor UV absorbing compounds as ginkgolides (A-C and J) and bilobalide in the extract. The complementary information between the UV and ELS fingerprint is evaluated. The ELS detector used in this study can operate in an impactor 'on' or 'off' mode. For each mode, the operating conditions such as the nebulizing gas flow rate, the drift tube temperature and the gain are optimized by use of three-level screening designs to obtain the best signal-to-noise (S/N) ratio in the final ELS fingerprint chromatogram. In both impactor modes, very similar S/N ratios are obtained for the nominal levels of the design. However, optimization of the operating conditions resulted, for both impactor modes, in a significant increase in S/N ratios compared to the initial evaluated conditions, obtained from the detector software.

Fast separations on monolithic silica columns: method transfer, robustness and column ageing for some case studies.

Six separation methods, developed on conventional silica high performance liquid chromatography (HPLC) columns were transferred to monolithic silica columns of 5 and 10 cm length. The transferred methods include the separation of an alkylbenzene mixture, the separations of drugs from their impurities (nimesulide, tetracycline, phenoxymethylpenicillin and erythromycin) and the separation of a green tea extract. The transfer of the first three methods was successful while for the latter three it was not. Increasing the flow rate up to 9 ml/min (where possible) inversely decreased the analysis time of the successfully transferred methods to 48 s (alkylbenzene mixture) 1.8 min (nimesulide mixture) and 3 min (tetracycline mixture) while still reasonable well separated peaks were obtained. The robustness and repeatability of the transferred and accelerated separations was found to be acceptable. Despite the use of flow rates up to 9 ml/min and frequent mobile phase changes with pH values varying from 3.5 to 7, the column performance was found to be rather constant and the column ageing to be minimal.