Development and integration of block operations for data invariant automation of digital preprocessing and analysis of biological and biomedical Raman spectra

Enhancing real-time cell culture monitoring: Automated Raman model optimization with Taguchi method

Raman spectroscopy has found widespread usage in monitoring cell culture processes both in research and practical applications. However, commonly, preprocessing methods, spectral regions, and modeling parameters have been chosen based on experience or trial-and-error strategies. These choices can significantly impact the performance of the models. There is an urgent need for a simple, effective, and automated approach to determine a suitable procedure for constructing accurate models. This paper introduces the adoption of a design of experiment (DoE) method to optimize partial least squares models for measuring the concentration of different components in cell culture bioreactors. The experimental implementation utilized the orthogonal test table L25(56). Within this framework, five factors were identified as control variables for the DoE method: the window width of Savitzky-Golay smoothing, the baseline correction method, the order of preprocessing steps, spectral regions, and the number of latent variables. The evaluation method for the model was considered as a factor subject to noise. The optimal combination of levels was determined through the signal-to-noise ratio response table employing Taguchi analysis. The effectiveness of this approach was validated through two cases, involving different cultivation scales, different Raman spectrometers, and different analytical components. The results consistently demonstrated that the proposed approach closely approximated the global optimum, regardless of data set size, predictive components, or the brand of Raman spectrometer. The performance of models recommended by the DoE strategy consistently surpassed those built using raw data, underscoring the reliability of models generated through this approach. When compared to exhaustive all-combination experiments, the DoE approach significantly reduces calculation times, making it highly practical for the implementation of Raman spectroscopy in bioprocess monitoring.

Rapid Vector-Based Peak Fitting and Resolution Enhancement for Correlation Analyses of Raman Hyperspectra

Spectroscopic peak parameters are important since they provide information about the analyte under study. Besides obtaining these parameters, peak fitting also resolves overlapped peaks. Thus, the obtained parameters should permit the construction of a higher-resolution version of the original spectrum. However, peak fitting is not an easy task due to computational reasons and because the true nature of the analyte is often unknown. These difficulties are major impediments when large hyperspectral data sets need to be processed rapidly, such as for manufacturing process control. We have developed a novel and relatively fast two-part algorithm to perform peak fitting and resolution enhancement on such data sets. In the first part of the algorithm, estimates of the total number of bands and their parameters were obtained from a representative spectrum in the data set, using a combination of techniques. Starting with these parameter estimates, all the spectra were then iteratively and rapidly fitted with Gaussian bands, exploiting intrinsic features of the Gaussian distribution with vector operations. The best fits for each spectrum were retained. By reducing the obtained bandwidths and commensurately increasing their amplitudes, high-resolution spectra were constructed that greatly improved correlation-based analyses. We tested the performance of the algorithm on synthetic spectra to confirm that this method could recover the ground truth correlations between highly overlapped peaks. To assess effective peak resolution, the method was applied to low-resolution spectra of glucose and compared to results from high-resolution spectra. We then processed a larger spectral data set from mammalian cells, fixed with methanol or air drying, to demonstrate the resolution enhancement of the algorithm on complex spectra and the effects of resolution-enhanced spectra on two-dimensional correlation spectroscopy and principal component analyses. The results indicated that the algorithm would allow users to obtain high-resolution spectra relatively fast and permit the recovery of important aspects of the data's intrinsic correlation structure.

Blind source separation of coexisting background in Raman spectra

Due to the Raman signal coexists with other scatter spectra which leads to the low ratio of the wanted signal and high background, the appropriate method should be applied to enhance this ratio. The nature of raw spectra is a multi-source system, so its determinacy must be ensured by multi-input. Besides, the faithfulness of output should be provided. Then, the huge fall within the frequencies of Raman and background almost satisfies separating demand for independent component analysis (ICA), and this analysis can give help to the achievement of the two type signals classing and estimate the optimal number of source and match ICA output signals to Raman or background. Thus, based on ICA and the mixing-entropy criteria, the background and Raman adapting calibration kit (BRACK) method is proposed, which is a kind of multiple raw spectral inputs and multiple output (MIMO) method. This method firstly divides the raw data into two parts of Raman and background by ICA, identifies Raman signal by entropy criterion, then restores the part of Raman signal. BRACK method obtains several advantages, for instance, well-adapted, no need for any additional option or extra-intervention, high fidelity, and no unwanted external information. In principle, the correction of background and Raman signals can be expected to be completed by BRACK method.

Applications of Raman spectroscopy in the development of cell therapies: state of the art and future perspectives

This comprehensive review article discusses current and future perspectives of Raman spectroscopy-based analyses of cell therapy processes and products.

Developing Fully Automated Quality Control Methods for Preprocessing Raman Spectra of Biomedical and Biological Samples

Spectral preprocessing is frequently required to render Raman spectra useful for further processing and analyses. The various preprocessing steps, individually and sequentially, are increasingly being automated to cope with large volumes of data from, for example, hyperspectral imaging studies. Full automation of preprocessing is especially desirable when it produces consistent results and requires minimal user input. It is therefore essential to evaluate the "quality" of such preprocessed spectra. However, relatively few methods exist to evaluate preprocessing quality, and fully automated methods for doing so are virtually non-existent. Here we provide a brief overview of fully automated spectral preprocessing and fully automated quality assessment of preprocessed spectra. We follow this with the introduction of fully automated methods to establish figures-of-merit that encapsulate preprocessing quality. By way of illustration, these quantitative methods are applied to simulated and real Raman spectra. Quality factor and quality parameter figures-of-merit resulting from individual preprocessing step quality tests, as well as overall figures-of-merit, were found to be consistent with the quality of preprocessed spectra.

Multivariate Analysis of Mixed Lipid Aggregate Phase Transitions Monitored Using Raman Spectroscopy

The phase behavior of aqueous 1,2-dimyristoyl-sn-glycero-3-phosphorylcholine (DMPC)/1,2-dihexanoyl-sn-glycero-3-phosphocholine (DHPC) mixtures between 8.0 ℃ and 41.0 ℃ were monitored using Raman spectroscopy. Temperature-dependent Raman matrices were assembled from series of spectra and subjected to multivariate analysis. The consensus of pseudo-rank estimation results is that seven to eight components account for the temperature-dependent changes observed in the spectra. The spectra and temperature response profiles of the mixture components were resolved by applying a variant of the non-negative matrix factorization (NMF) algorithm described by Lee and Seung (1999). The rotational ambiguity of the data matrix was reduced by augmenting the original temperature-dependent spectral matrix with its cumulative counterpart, i.e., the matrix formed by successive integration of the spectra across the temperature index (columns). Successive rounds of constrained NMF were used to isolate component spectra from a significant fluorescence background. Five major components exhibiting varying degrees of gel and liquid crystalline lipid character were resolved. Hydrogen-bonded water networks exhibiting varying degrees of organization are associated with the lipid components. Spectral parameters were computed to compare the chain conformation, packing, and hydration indicated by the resolved spectra. Based on spectral features and relative amounts of the components observed, four components reflect long chain lipid response. The fifth component could reflect the response of the short chain lipid, DHPC, but there were no definitive spectral features confirming this assignment. A minor component of uncertain assignment that exhibits a striking response to the DMPC pre-transition and chain melting transition also was recovered. While none of the spectra resolved exhibit features unequivocally attributable to a specific aggregate morphology or step in the gelation process, the results are consistent with the evolution of mixed phase bicelles (nanodisks) and small amounts of worm-like DMPC/DHPC aggregates, and perhaps DHPC micelles, at low temperature to suspensions of branched and entangled worm-like aggregates above the DMPC gel phase transition and perforated multi-lamellar aggregates at high temperature.

Empirical Factors Affecting the Quality of Non-Negative Matrix Factorization of Mammalian Cell Raman Spectra

Mammalian cells contain various macromolecules that can be investigated non-invasively with Raman spectroscopy. The particular mixture of major macromolecules present in a cell being probed are reflected in the measured Raman spectra. Determining macromolecular identities and estimating their concentrations from these mixture Raman spectra can distinguish cell types and otherwise enable biological research. However, the application of canonical multivariate methods, such as principal component analysis (PCA), to perform spectral unmixing yields mathematical solutions that can be difficult to interpret. Non-negative matrix factorization (NNMF) improves the interpretability of unmixed macromolecular components, but can be difficult to apply because ambiguities produced by overlapping Raman bands permit multiple solutions. Furthermore, theoretically sound methods can be difficult to implement in practice. Here we examined the effects of a number of empirical approaches on the quality of NNMF results. These approaches were evaluated on simulated mammalian cell Raman hyperspectra and the results were used to develop an enhanced procedure for implementing NNMF. We demonstrated the utility of this procedure using a Raman hyperspectral data set measured from human islet cells to recover the spectra of insulin and glucagon. This was compared to the relatively inferior PCA of these data.