A down-to-earth analyst view of rotational ambiguity in second-order calibration with multivariate curve resolution - a tutorial

Enhancing gas chromatography-mass spectrometry resolution and pure analyte discovery using intra-chromatogram elution profile matching

Chemometric decomposition methods like multivariate curve resolution-alternating least squares (MCR-ALS) are often employed in gas chromatography-mass spectrometry (GC-MS) to improve analyte identification and quantitation. However, these methods can perform poorly for analytes with a low chromatographic resolution (Rs) and a high degree of spectral contamination from noise and background interferences. Thus, we propose a novel computational algorithm, termed mzCompare, to improve analyte identification and quantitation when coupled to MCR-ALS. The mzCompare method utilizes an underlying requirement that the retention time and peak shape between mass channels (m/z) of the same analyte should be similar. By discovering the selective m/z for a given analyte in a chromatogram, a pure elution profile can be generated and used as an equality constraint in MCR-ALS. The performance of the mzCompare methodology is demonstrated with both experimental and simulated chromatograms. Experimentally, unresolved analytes with a Rs as low as 0.05 could be confidently identified with mzCompare assisted MCR-ALS. Furthermore, application of the mzCompare algorithm to a complex aerospace fuel resulted in the discovery of 335 analytes, a 44 % increase compared to conventional peak detection methods. GC-MS simulations of target-interferent analyte pairs demonstrated that the performance of MCR-ALS deteriorated below a Rs of ∼0.25. However, mzCompare assisted MCR-ALS showed excellent identification and acceptable quantitative accuracy at a Rs of ∼0.02. These results show that the mzCompare algorithm can help analysts overcome modeling ambiguities resulting from the chemometric multiplex disadvantage.

Application of smart chemometric models for spectra resolution and determination of challenging multi-action quaternary mixture: statistical comparison with greenness assessment

A multivariate spectrophotometric method is a potential approach that enables discrimination of spectra of components in complex matrices (e.g., pharmaceutical formulation) serving as a substitution method for chromatography. Four green smart multivariate spectrophotometric models were proposed and validated, including principal component regression (PCR), partial least-squares (PLS), multivariate curve resolution-alternating least squares (MCR-ALS), and artificial neural networks (ANN). The developed chemometric models were compared to resolve highly overlapping spectra of Paracetamol (PARA), Chlorpheniramine maleate (CPM), Caffeine (CAF), and Ascorbic acid (ASC). The four multivariate calibration models were assessed via recoveries percent, and root mean square error of prediction. Hence, the proposed models were efficiently applied with no need for any preliminary separation step. The models were utilized to analyze the studied components in their combined pharmaceutical formulation (Grippostad® C capsules). Analytical GREEnness Metric Approach (AGREE) and eco-scale tools were applied to assess the greenness of the established models and found to be 0.77 and 85, respectively. Moreover, the proposed models have been compared to official ones showing no considerable variations in accuracy and precision. Therefore, these models can be highly advantageous for conducting standard pharmaceutical analysis of the substances researched within product testing laboratories.

A new constraint to model background signals when processing chromatographic-spectral second-order data with multivariate curve resolution

BACKGROUND

the chemometric processing of second-order chromatographic-spectral data is usually carried out with the aid of multivariate curve resolution-alternating least-squares (MCR-ALS). When baseline contributions occur in the data, the background profile retrieved with MCR-ALS may show abnormal lumps or negative dips at the position of the remaining component peaks.

RESULTS

The phenomenon is shown to be due to remaining rotational ambiguity in the obtained profiles, as confirmed by the estimation of the boundaries of the range of feasible bilinear profiles. To avoid the abnormal features in the retrieved profile, a new background interpolation constraint is proposed and described in detail. Both simulated and experimental data are employed to support the need of the new MCR-ALS constraint. In the latter case, the estimated analyte concentrations agreed with those previously reported.

SIGNIFICANCE

The developed procedure helps to reduce the extent of rotational ambiguity in the solution and to better interpret the results on physicochemical grounds.

A tutorial on multi-way data processing of excitation-emission fluorescence matrices acquired from semiconductor quantum dots sensing platforms

This tutorial demonstrates how to exploit the second-order advantage on excitation-emission fluorescence matrices (EEFMs) acquired from sensing platforms based on analyte-triggered semiconductor quantum dots (QDs) fluorescence modulation (quenching/enhancing). The advantage in processing such second-order EEFMs data from complex samples, seeking successful quantification, is comprehensively addressed. It is worth emphasizing that, aiming to exploit the second-order advantage, the selection of the most appropriate advanced chemometric model should rely on the matching between the data structure and the physicochemical chemometric model assumption. In this sense, the achievement of second-order advantage after EEFMs' processing is extensively addressed throughout this tutorial taking into consideration three different analytical situations, each involving a specific data structure: i) parallel factor analysis (PARAFAC), which is applied in a real dataset stacked in a three-way data array containing a trilinear data structure acquired from QDs-based detection with non-selective species; ii) multivariate curve resolution - alternating least-squares (MCR-ALS), which is also employed in a real dataset arranged in an augmented data matrix containing non-trilinear data structure acquired from QDs-based detection with a single breaking mode caused by background signals; iii) unfolded partial least-squares with residual bilinearization (U-PLS/RBL), which is applied in a dataset containing non-trilinear data acquired from a classical fluorescence system with two breaking modes caused by inner filter effect (IFE) in both instrumental modes (excitation and emission). The latter challenging data structure can be acquired via fluorescence quenching from IFE-based sensing platforms and chemometrically handled in two main steps. First, a set of calibration EEFMs data is converted into an unfolded data matrix during the unfolding process, followed by applying U-PLS model. Second, a post-calibration procedure using RBL analysis is applied to a test sample of EEFM maintained in its matrix form, in order to handle potential interferents. In the last section, the state-of-the-art of second-order EEFMs data acquired from semiconductor QDs-based sensing platforms and coupled to multi-way fluorescence data processing to accomplish a successful quantification, even with substantial interfering species, is critically reviewed.

Eucommia ulmoides Oliver-Tribulus terrestris L. Drug Pair Regulates Ferroptosis by Mediating the Neurovascular-Related Ligand-Receptor Interaction Pathway- A Potential Drug Pair for Treatment Hypertension and Prevention Ischemic Stroke

Background

In this study, we used the network pharmacology approach to explore the potential disease targets of the Eucommia ulmoides Oliver (EUO)-Tribulus terrestris L. (TT) drug pair in the treatment of hypertension-associated neurovascular lesions and IS via the ferroptosis pathway.

Methods

We used the Traditional Chinese Medicine Systems Pharmacology Database and Analysis Platform to search for the key active compounds and targets of the drug pair. Based on the GeneCards database, the relevant targets for the drug pair were obtained. Then, we performed the molecular docking of the screened core active ingredients and proteins using the DAVID database and the R AutoDock Vina software. Based on the GSE22255 dataset, these screened target proteins were used to build random forest (RF) and support vector machine (SVM) models. Finally, a new IS nomogram prediction model was constructed and evaluated.

Results

There were 36 active compounds in the EUO-TT drug pair. CHRM1, NR3C1, ADRB2, and OPRD1 proteins of the neuroactive ligand-receptor interaction pathway interacted with the proteins related to the ferroptosis pathway. Molecular docking experiments identified 12 active ingredients of the drug pair that may tightly bind to those target proteins. We constructed a visual IS nomogram prediction model using four genes (CHRM1, NR3C1, ADRB2, and OPRD1). The calibration curve, DCA, and clinical impact curves all indicated that the nomogram model is clinically applicable and diagnostically capable. CHRM1, NR3C1, ADRB2, and OPRD1, the target genes of the four effective components of the EUO-TT drug pair, were considered as risk markers for IS.

Conclusions

The active ingredients of EUO-TT drug pair may act on proteins associated with the neuroactive ligand-receptor interaction pathway to regulate ferroptosis in vascular neurons cells, ultimately affecting the onset and progression of hypertension.

Exercise Modifies the Transcriptional Regulatory Features of Monocytes in Alzheimer's Patients: A Multi-Omics Integration Analysis Based on Single Cell Technology

Monocytes have been reported to be important mediators of the protective effect of exercise against the development of Alzheimer's disease (AD). This study aims explored the mechanism by which monocytes achieve this. Using single cell transcriptome analysis, results showed that CD14 + and CD16 + monocytes interacted with other cells in the circulating blood. TNF, CCR1, APP, and AREG, the key ligand-receptor-related genes, were found to be differentially expressed between exercise-treated and AD patients. The SCENIC analysis was performed to identify individual clusters of the key transcription factors (TFs). Nine clusters (M1-M9) were obtained from the co-expression network. Among the identified TFs, MAFB, HES4, and FOSL1 were found to be differentially expressed in AD. Moreover, the M4 cluster to which MAFB, HES4, and FOSL1 belonged was defined as the signature cluster for AD phenotype. Differential analysis by bulkRNA-seq revealed that the expression of TNF, CCR1, and APP were all upregulated after exercise (p < 0.05). And ATF3, MAFB, HES4, and KLF4 that were identified in M4 clusters may be the TFs that regulate TNF, CCR1, and APP in exercise prescription. After that, APP, CCR1, TNF, ATF3, KLF4, HES4, and MAFB formed a regulatory network in the ERADMT gene set, and all of them were mechanistically linked. The ERADMT gene set has been found to be a potential risk marker for the development of AD and can be used as an indicator of compliance to exercise therapy in AD patients. Using single-cell integration analysis, a network of exercise-regulating TFs in monocytes was constructed for AD disease. The constructed network reveals the mechanism by which exercise regulated monocytes to confer therapeutic benefits against AD and its complications. However, this study, as a bioinformatic research, requires further experimental validation.

Chromatographic Applications in the Multi-Way Calibration Field

In this review, recent advances and applications using multi-way calibration protocols based on the processing of multi-dimensional chromatographic data are discussed. We first describe the various modes in which multi-way chromatographic data sets can be generated, including some important characteristics that should be taken into account for the selection of an adequate data processing model. We then discuss the different manners in which the collected instrumental data can be arranged, and the most usually applied models and algorithms for the decomposition of the data arrays. The latter activity leads to the estimation of surrogate variables (scores), useful for analyte quantitation in the presence of uncalibrated interferences, achieving the second-order advantage. Recent experimental reports based on multi-way liquid and gas chromatographic data are then reviewed. Finally, analytical figures of merit that should always accompany quantitative calibration reports are described.

Achieving the analytical second-order advantage with non-bilinear second-order data

Multi-way calibration based on second-order data constitutes a revolutionary milestone for analytical applications. However, most classical chemometric models assume that these data fulfil the property of low rank bilinearity, which cannot be accomplished by all instrumental methods. Indeed, various techniques are able to generate non-bilinear data, which are all potentially useful for the development of novel second-order calibration methodologies. However, the achievement of the second-order advantage in these cases may be severely limited, since methods for comprehensive modelling of non-bilinear second-order data remain only partially explored. In this research, the analytical performance of three well-known second-order models, namely non-bilinear rank annihilation (NBRA), unfolded partial least-squares with residual bilinearization (U-PLS-RBL) and multivariate curve resolution - alternating least-squares (MCR-ALS) is systematically assessed through sets of simulated and experimental non-bilinear second-order data, involving one analyte and one interferent. Although it is not possible to establish a single strategy to model any type of non-bilinear second-order data with the studied methods, each approach may lead to successful predictions under certain circumstances. It is shown that the prediction capacity is severely affected by data properties such as the level of instrumental noise, the rank of the response matrices and the signal selectivity pattern of the analyte.