Predicting liquid chromatography-electrospray ionization/mass spectrometry signal from the structure of model compounds and experimental factors; case study of aripiprazole and its impurities.

A priori estimation of analyte response is crucial for the efficient development of liquid chromatography-electrospray ionization/mass spectrometry (LC-ESI/MS) methods, but remains a demanding task given the lack of knowledge about the factors affecting the experimental outcome. In this research, we address the challenge of discovering the interactive relationship between signal response and structural properties, method parameters and solvent-related descriptors throughout an approach featuring quantitative structure-property relationship (QSPR) and design of experiments (DoE). To systematically investigate the experimental domain within which QSPR prediction should be undertaken, we varied LC and instrumental factors according to the Box-Behnken DoE scheme. Seven compounds, including aripiprazole and its impurities, were subjected to 57 different experimental conditions, resulting in 399 LC-ESI/MS data endpoints. To obtain a more standard distribution of the measured response, the peak areas were log-transformed before modeling. QSPR predictions were made using features selected by Genetic Algorithm (GA) and providing Gradient Boosted Trees (GBT) with training data. Proposed model showed satisfactory performance on test data with a RMSEP of 1.57 % and a of 96.48 %. This is the first QSPR study in LC-ESI/MS that provided a holistic overview of the analyte's response behavior across the experimental and chemical space. Since intramolecular electronic effects and molecular size were given great importance, the GA-GBT model improved the understanding of signal response generation of model compounds. It also highlighted the need to fine-tune the parameters affecting desolvation and droplet charging efficiency.

Analytical Quality by Design: Achieving Robustness of an LC-CAD Method for the Analysis of Non-Volatile Fatty Acids.

An alternative to the time-consuming and error-prone pharmacopoeial gas chromatography method for the analysis of fatty acids (FAs) is urgently needed. The objective was therefore to propose a robust liquid chromatography method with charged aerosol detection for the analysis of polysorbate 80 (PS80) and magnesium stearate. FAs with different numbers of carbon atoms in the chain necessitated the use of a gradient method with a Hypersil Gold C18 column and acetonitrile as organic modifier. The risk-based Analytical Quality by Design approach was applied to define the Method Operable Design Region (MODR). Formic acid concentration, initial and final percentages of acetonitrile, gradient elution time, column temperature, and mobile phase flow rate were identified as critical method parameters (CMPs). The initial and final percentages of acetonitrile were fixed while the remaining CMPs were fine-tuned using response surface methodology. Critical method attributes included the baseline separation of adjacent peaks (α-linolenic and myristic acid, and oleic and petroselinic acid) and the retention factor of the last compound eluted, stearic acid. The MODR was calculated by Monte Carlo simulations with a probability equal or greater than 90%. Finally, the column temperature was set at 33 °C, the flow rate was 0.575 mL/min, and acetonitrile linearly increased from 70 to 80% (v/v) within 14.2 min.

The secret of reversed-phase/weak cation exchange retention mechanisms in mixed-mode liquid chromatography applied for small drug molecule analysis.

Resolving complex sample mixtures by liquid chromatography in a single run is challenging. The so-called mixed-mode liquid chromatography (MMLC) which combines several retention mechanisms within a single column, can provide resource-efficient separation of solutes of diverse nature. The Acclaim Mixed-Mode WCX-1 column, encompassing hydrophobic and weak cation exchange interactions, was employed for the analysis of small drug molecules. The stationary phase's interaction abilities were assessed by analysing molecules of different ionisation potentials. Mixed Quantitative Structure-Retention Relationship (QSRR) models were developed for revealing significant experimental parameters (EPs) and molecular features governing molecular retention. According to the plan of Face-Centred Central Composite Design, EPs (column temperature, acetonitrile content, pH and buffer concentration of aqueous mobile phase) variations were included in QSRR modelling. QSRRs were developed upon the whole data set (global model) and upon discrete parts, related to similarly ionized analytes (local models) by applying gradient boosted trees as a regression tool. Root mean squared errors of prediction for global and local QSRR models for cations, anions and neutrals were respectively 0.131; 0.105; 0.102 and 0.042 with the coefficient of determination 0.947; 0.872; 0.954 and 0.996, indicating satisfactory performances of all models, with slightly better accuracy of local ones. The research showed that influences of EPs were dependant on the molecule's ionisation potential. The molecular descriptors highlighted by models pointed out that electrostatic and hydrophobic interactions and hydrogen bonds participate in the retention process. The molecule's conformation significance was evaluated along with the topological relationship between the interaction centres, explicitly determined for each molecular species through local models. All models showed good molecular retention predictability thus showing potential for facilitating the method development.

Charged aerosol detector response modeling for fatty acids based on experimental settings and molecular features: a machine learning approach.

The charged aerosol detector (CAD) is the latest representative of aerosol-based detectors that generate a response independent of the analytes' chemical structure. This study was aimed at accurately predicting the CAD response of homologous fatty acids under varying experimental conditions. Fatty acids from C12 to C18 were used as model substances due to semivolatile characterics that caused non-uniform CAD behaviour. Considering both experimental conditions and molecular descriptors, a mixed quantitative structure-property relationship (QSPR) modeling was performed using Gradient Boosted Trees (GBT). The ensemble of 10 decisions trees (learning rate set at 0.55, the maximal depth set at 5, and the sample rate set at 1.0) was able to explain approximately 99% (Q2: 0.987, RMSE: 0.051) of the observed variance in CAD responses. Validation using an external test compound confirmed the high predictive ability of the model established (R2: 0.990, RMSEP: 0.050). With respect to the intrinsic attribute selection strategy, GBT used almost all independent variables during model building. Finally, it attributed the highest importance to the power function value, the flow rate of the mobile phase, evaporation temperature, the content of the organic solvent in the mobile phase and the molecular descriptors such as molecular weight (MW), Radial Distribution Function-080/weighted by mass (RDF080m) and average coefficient of the last eigenvector from distance/detour matrix (Ve2_D/Dt). The identification of the factors most relevant to the CAD responsiveness has contributed to a better understanding of the underlying mechanisms of signal generation. An increased CAD response that was obtained for acetone as organic modifier demonstrated its potential to replace the more expensive and environmentally harmful acetonitrile.

Comparative chemometric and quantitative structure-retention relationship analysis of anisotropic lipophilicity of 1-arylsuccinimide derivatives determined in high-performance thin-layer chromatography system with aprotic solvents.

Numerous structurally different amides and imides including succinimide derivatives exhibit diverse bioactive potential. The development of new compounds requires rationalization in the design in order to provide structural changes that guarantee favorable physico-chemical properties, pharmacological activity and safety. In the present research, a comprehensive study with comparison of the chromatographic lipophilicity and other physico-chemical properties of five groups of 1-arylsuccinimide derivatives was conducted. The chemometric analysis of their physico-chemical properties was carried out by using unsupervised (hierarchical cluster analysis and principal component analysis) and supervised pattern recognition methods (linear discriminant analysis), while the correlations between the in silico molecular features and chromatographic lipophilicity were examined applying linear and non-linear Quantitative Structure-Retention Relationship (QSRR) approaches. The main aim of the conducted research was to determine similarities and dissimilarities among the studied 1-arylsuccinimides, to point out the molecular features which have significant influence on their lipophilicity, as well as to establish high-quality QSRR models which can be used in prediction of chromatographic lipophilicity of structurally similar 1-arylsuccinimides. This study is a continuation of analysis and determination of the physico-chemical properties of 1-arylsuccinimides which could be important guidelines in further in vitro and eventually in vivo studies of their biological potential.

Performance comparison of nonlinear and linear regression algorithms coupled with different attribute selection methods for quantitative structure - retention relationships modelling in micellar liquid chromatography.

In micellar liquid chromatography (MLC), the addition of a surfactant to the mobile phase in excess is accompanied by an alteration of its solubilising capacity and a change in the stationary phase's properties. As an implication, the prediction of the analytes' retention in MLC mode becomes a challenging task. Mixed Quantitative Structure - Retention Relationships (QSRR) modelling represents a powerful tool for estimating the analytes' retention. This study compares 48 successfully developed mixed QSRR models with respect to their ability to predict retention of aripiprazole and its five impurities from molecular structures and factors that describe the Brij - acetonitrile system. The development of the models was based on an automatic combining of six attribute (feature) selection methods with eight predictive algorithms and the optimization of hyper-parameters. The feature selection methods included Principal Component Analysis (PCA), Non-negative Matrix Factorization (NMF), ReliefF, Multiple Linear Regression (MLR), Mutual Info and F-Regression. The series of investigated predictive algorithms comprised Linear Regressions (LR), Ridge Regression, Lasso Regression, Artificial Neural Networks (ANN), Support Vector Regression (SVR), Random Forest (RF), Gradient Boosted Trees (GBT) and K-Nearest neighbourhood (k-NN). A sufficient amount of data for building the model (78 cases in total) was provided by conducting 13 experiments for each of the 6 analytes and collecting the target responses afterwards. Different experimental settings were established by varying the values of the concentration of Brij L23, pH of the aqueous phase and acetonitrile content in the mobile phase according to the Box-Behnken design. In addition to the chromatographic parameters, the pool of independent variables was expanded by 27 molecular descriptors from all major groups (physicochemical, quantum chemical, topological and spatial structural descriptors). The best model was chosen by taking into consideration the Root Mean Square Error (RMSE) and cross-validation (CV) correlation coefficient (Q2) values. Interestingly, the comparative analysis indicated that a change in the set of input variables had a minor impact on the performance of the final models. On the other hand, different regression algorithms showed great diversity in the ability to learn patterns conserved in the data. In this regard, testing many regression algorithms is necessary in order to find the most suitable technique for model building. In the specific case, GBT-based models have demonstrated the best ability to predict the retention factor in the MLC mode. Steric factors and dipole-dipole interactions have proven to be relevant to the observed retention behaviour. This study, although being of a smaller scale, is a most promising starting point for comprehensive MLC retention prediction.

Quantitative structure-property relationship modeling of polar analytes lacking UV chromophores to charged aerosol detector response.

In this study, a quantitative structure-property relationship model was built in order to link molecular descriptors and chromatographic parameters as inputs towards CAD responsiveness. Aminoglycoside antibiotics, sugars, and acetylated amino sugars, which all lack a UV/vis chromophore, were selected as model substances due to their polar nature that represents a challenge in generating a CAD response. Acetone, PFPA, flow rate, data rate, filter constant, SM5_B(s), ATS7s, SpMin1_Bh(v), Mor09e, Mor22e, E1u, R7v+, and VP as the most influential inputs were correlated with the CAD response by virtue of ANN applying a backpropagation learning rule. External validation on previously unseen substances showed that the developed 13-6-3-1 ANN model could be used for CAD response prediction across the examined experimental domain reliably (R2 0.989 and RMSE 0.036). The obtained network was used to reveal CAD response correlations. The impact of organic modifier content and flow rate was in accordance with the theory of the detector's functioning. Additionally, the significance of SpMin1_Bh(v) aided in emphasizing the often neglected surface-dependent CAD character, while the importance of Mor22e as a molecular descriptor accentuated its dependency on the number of electronegative atoms taking part in charging the formed particles. The significance of PFPA demonstrated the possibility of using evaporative chaotropic reagents in CAD response improvement when dealing with highly polar substances that act as kosmotropes. The network was also used in identifying possible interactions between the most significant inputs. A joint effect of PFPA and acetone was shown, representing a good starting point for further investigation with different and, especially, eco-friendly organic solvents and chaotropic agents in the routine application of CAD.

Multicriteria Optimization Methodology in Stability-Indicating Method Development of Cilazapril and Hydrochlorothiazide.

Multicriteria optimization methodology was applied in development of UHPLC-UV-MS method for separation of cilazapril, hydrochlorothiazide and their degradation products. This method is also applicable for analysis of cilazapril, hydrochlorothiazide and their degradation products in combined tablet formulation. Prior to method optimization forced degradation studies were conducted. Cilazapril and hydrochlorothiazide were subjected to acidic (0.1, 0.5 and 1.0 M HCl), basic (0.1, 0.5 and 1.0 M NaOH), thermal (70°C), oxidative (3-30% H2O2) degradation and photodegradation (day light). Cilazapril appeared to be unstable toward acid and base and resulted in formation of cilazaprilat. Hydrochlorothiazide significantly degraded after acid, base and thermal hydrolysis and formed degradation product was 4-amino-6-chlorobenzene-1.3-disulfonamide. For both substances, after oxidative degradation unknown products have arisen. Initial percentage of acetonitrile in mobile phase, final percentage of acetonitrile in mobile phase, time of gradient elution and column temperature were defined as variables to be optimized toward two chromatographic responses by means of central composite design and Derringer's desirability function. The satisfactory chromatographic analysis was achieved on Kinetex C18 (2.6 µm, 50 × 2.1 mm) column with temperature set at 25°C. The final mobile phase consisted of acetonitrile and 20 mM ammonium formate buffer (pH adjusted to 8.5). The flow rate of the mobile phase was 400 µL min-1 and it was pumped in a gradient elution mode.