Quantitative structure –retention relationship modeling of selected antipsychotics and their impurities in green liquid chromatography using cyclodextrin mobile phases

A Review of Machine Learning and QSAR/QSPR Predictions for Complexes of Organic Molecules with Cyclodextrins

Cyclodextrins are macrocyclic rings composed of glucose residues. Due to their remarkable structural properties, they can form host-guest inclusion complexes, which is why they are frequently used in the pharmaceutical, cosmetic, and food industries, as well as in environmental and analytical chemistry. This review presents the reports from 2011 to 2023 on the quantitative structure-activity/property relationship (QSAR/QSPR) approach, which is primarily employed to predict the thermodynamic stability of inclusion complexes. This article extensively discusses the significant developments related to the size of available experimental data, the available sets of descriptors, and the machine learning (ML) algorithms used, such as support vector machines, random forests, artificial neural networks, and gradient boosting. As QSAR/QPR analysis only requires molecular structures of guests and experimental values of stability constants, this approach may be particularly useful for predicting these values for complexes with randomly substituted cyclodextrins, as well as for estimating their dependence on pH. This work proposes solutions on how to effectively use this knowledge, which is especially important for researchers who will deal with this topic in the future. This review also presents other applications of ML in relation to CD complexes, including the prediction of physicochemical properties of CD complexes, the development of analytical methods based on complexation with CDs, and the optimisation of experimental conditions for the preparation of the complexes.

Antiemetic effects of sclareol, possibly through 5-HT3 and D2 receptor interaction pathways: In-vivo and in-silico studies

BACKGROUND

Emesis is a complex physiological phenomenon that serves as a defense against numerous toxins, stressful situations, adverse medication responses, chemotherapy, and movement. Nevertheless, preventing emesis during chemotherapy or other situations is a significant issue for researchers. Hence, the majority view contends that successfully combining therapy is the best course of action. In-vivo analysis offers a more comprehensive grasp of how compounds behave within a complex biological environment, whereas in-silico evaluation refers to the use of computational models to forecast biological interactions.

OBJECTIVES

The objectives of the present study were to evaluate the effects of Sclareol (SCL) on copper sulphate-induced emetic chicks and to investigate the combined effects of these compounds using a conventional co-treatment approach and in-silico study.

METHODS

SCL (5, 10, and 15 mg/kg) administered orally with or without pre-treatment with anti-emetic drugs (Ondansetron (ODN): 24 mg/kg, Domperidone (DOM): 80 mg/kg, Hyoscine butylbromide (HYS): 100 mg/kg, and Promethazine hydrochloride (PRO): 100 mg/kg) to illustrate the effects and the potential involvement with 5HT3, D2, M3/AChM, H1, or NK1 receptors by SCL. Furthermore, an in-silico analysis was conducted to forecast the role of these receptors in the emetic process.

RESULTS

The results suggest that SCL exerted a dose-dependent anti-emetic effect on the chicks. Pretreatment with SCL-10 significantly minimized the number of retches and lengthened the emesis tendency of the experimental animals. SCL-10 significantly increased the anti-emetic effects of ODN and DOM. However, compared to the ODN-treated group, (SCL-10 + ODN) group considerably (p < 0.0001) extended the latency duration (109.40 ± 1.03 s) and significantly (p < 0.01) decreased the number of retches (20.00 ± 0.70), indicating an anti-emetic effect on the test animals. In in-silico analysis, SCL exhibited promising binding affinities with suggesting receptors.

CONCLUSION

SCL-10 exerted an inhibitory-like effect on emetic chicks, probably through the interaction of the 5HT3 and D2 receptors. Further studies are highly appreciated to validate this study and determine the precise mechanism(s) behind the anti-emetic effects of SCL. We expect that SCL-10 may be utilized as an antiemetic treatment in a single dosage form or that it may function as a synergist with other traditional medicines.

Development and optimization of stability-indicating method of ethinylestradiol, levonorgestrel, and their main impurities using quality by design approach

The association of Ethinylestradiol 0.03 mg and Levonorgestrel 0.15 mg is a hormonal contraceptive that combines estrogen and progestogen. According to a bibliographic survey, these combined drugs present at least 18 known degradation products, which are required to control the potential impurities harmful to human health. The high number of impurities and the low concentrations of the active pharmaceutical ingredients (APIs) and their respective degradation products increase the complexity of the stability-indicating method development for this medicine. Thus, this work aimed to develop and optimize the stability-indicating method using the quality by design (QbD) approach and in-silico tools for application in samples of oral contraceptives sold in Brazil. The analysis samples were initially subjected to a forced degradation study through 7 days of exposure under acid and alkali hydrolysis, oxidative condition, and oxidation by metal ions. In addition to the chemical exposure, the sample was subjected to physical stress through 10 days of exposure under dry heat, moisture, and photolytic degradation. These exposure samples were analyzed in the development and optimization of chromatographic conditions. As a result, the developed method was able to separate 20 known substances, including the two APIs and their respective 18 degradation products, as well as unknown degradation products obtained by the forced degradation study. Finally, this stability-indicating method was successfully applied for comparative analysis of contraceptive drugs marketed in Brazil, newly purchased and subjected to accelerated stability condition at 40 °C and 75% RH over the 6-month period.

Transfer of gas chromatographic retention data among poly(siloxane) columns by quantitative structure-retention relationships based on molecular descriptors of both solutes and stationary phases

In the present study, computational molecular descriptors of 90 saturated esters and seven poly(siloxane) stationary phases with different polarity (SE-30, OV-7, DC-710, OV-25, XE-60, OV-225 and Silar-5CP) were combined into quantitative structure-retention relationship (QSRR) models aimed at predicting the Kováts retention indices (RIs) of the solutes. The molecular descriptors (174) of the stationary phases included in the models were computed using Dragon software from poly(siloxane) oligomers made of 20 siloxane units reflecting the nominal composition of the stationary phase, whereas 439 molecular descriptors were adopted to represent the esters. Different QSRR models were generated by means of Partial Least Squares (PLS) regression to assess the accuracy of this approach in predicting the RIs of unexplored solutes both in known and external stationary phases. After calibration of each PLS model, the descriptors were selected/discarded according to their relevance, evaluated by Covariance Selection (CovSel), and the PLS models were re-built, which resulted in a noticeable improvement of their predictive ability. Firstly, all the available data were equally divided into a training and a test set; the model built on the calibration set was used to predict the RIs of the validation observations. Successively, seven diverse PLS models were created following a "leave-one-column-out" fashion procedure, each one finalized to the estimation of the RIs of the 90 esters associated with a single stationary phase, whereas the calibration model was calculated on the remaining data. All the estimated models provided successful results on the external stationary phase, and predictive performance further increased after variable selection based on CovSel analysis. The final models provided a Root Mean Square Error in Cross Validation (RMSECV) in the range 12-20, a Root Mean Square Error in Prediction (RMSEP) in the range 11-26, and Mean Absolute Percentage Errors in Prediction (MAMEPs) in the range 0.7-1.5, revealing accurate cross-column prediction. Eventually, to test the robustness of the proposed approach, the 90 solutes were equally partitioned into a calibration and a test set and two further QSSR strategies were applied. The first PLS model was calibrated on all the seven stationary phases and the RIs of the 45 external solutes in the same seven columns were simultaneously predicted. The last QSRR approach followed a "leave-one-column-out" scheme and RI of 45 test solutes on an external stationary phase was predicted by a PLS model calibrated with the data of the 45 remaining solutes and the six left stationary phases. After selection of the significant molecular descriptors, PLS regression provided RMSECV values in the range 6-19, RMSEPs in the range 10-14, and MAPEPs in the range 0.9-2.4, revealing the suitability of the approach to deduce the RI of unknown solutes in uncharted stationary phases.

Deep Q-learning for the selection of optimal isocratic scouting runs in liquid chromatography

An important challenge in chromatography is the development of adequate separation methods. Accurate retention models can significantly simplify and expedite the development of adequate separation methods for complex mixtures. The purpose of this study was to introduce reinforcement learning to chromatographic method development, by training a double deep Q-learning algorithm to select optimal isocratic scouting runs to generate accurate retention models. These scouting runs were fit to the Neue-Kuss retention model, which was then used to predict retention factors both under isocratic and gradient conditions. The quality of these predictions was compared to experimental data points, by computing a mean relative percentage error (MRPE) between the predicted and actual retention factors. By providing the reinforcement learning algorithm with a reward whenever the scouting runs led to accurate retention models and a penalty when the analysis time of a selected scouting run was too high (> 1h); it was hypothesized that the reinforcement learning algorithm should by time learn to select good scouting runs for compounds displaying a variety of characteristics. The reinforcement learning algorithm developed in this work was first trained on simulated data, and then evaluated on experimental data for 57 small molecules - each run at 10 different fractions of organic modifier (0.05 to 0.90) and four different linear gradients. The results showed that the MRPE of these retention models (3.77% for isocratic runs and 1.93% for gradient runs), mostly obtained via 3 isocratic scouting runs for each compound, were comparable in performance to retention models obtained by fitting the Neue-Kuss model to all (10) available isocratic datapoints (3.26% for isocratic runs and 4.97% for gradient runs) and retention models obtained via a "chromatographer's selection" of three scouting runs (3.86% for isocratic runs and 6.66% for gradient runs). It was therefore concluded that the reinforcement learning algorithm learned to select optimal scouting runs for retention modeling, by selecting 3 (out of 10) isocratic scouting runs per compound, that were informative enough to successfully capture the retention behavior of each compound.

Experimental design in HPLC separation of pharmaceuticals

Design of Experiments (DoE) is an indispensable tool in contemporary drug analysis as it simultaneously balances a number of chromatographic parameters to ensure optimal separation in High Pressure Liquid Chromatography (HPLC). This manuscript briefly outlines the theoretical background of the DOE and provides step-by-step instruction for its implementation in HPLC pharmaceutical practice. It particularly discusses the classification of various design types and their possibilities to rationalize the different stages of HPLC method development workflow, such as the selection of the most influential factors, factors optimization and assessment of the method robustness. Additionally, the application of the DOE-based Analytical Quality by Design (AQbD) concept in the LC method development has been summarized. Recent achievements in the use of DOE in the development of stability-indicating LC and hyphenated LC-MS methods have also been briefly reported. Performing of Quantitative structure retention relationship (QSRR) study enhanced with DOE-based data collection was recomended as a future perspective in description of retention in HPLC system.

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 (Q²) 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.

Characterization of bonded stationary phase performance as a function of qualitative and quantitative chromatographic factors in chaotropic chromatography with risperidone and its impurities as model substances

Numerous stationary phases have been developed with the aim to provide desired performances during chromatographic analysis of the basic solutes in their protonated form. In this work, the procedure for the characterization of bonded stationary phase performance, when both qualitative and quantitative chromatographic factors were varied in chaotropic chromatography, was proposed. Risperidone and its three impurities were selected as model substances, while acetonitrile content in the mobile phase (20-30%), the pH of the aqueous phase (3.00-5.00), the content of chaotropic agents in the aqueous phase (10-100 mM), type of chaotropic agent (NaClO₄, CF₃COONa), and stationary phase type (Zorbax Eclipse XDB, Zorbax Extend) were studied as chromatographic factors. The proposed procedure implies the combination of D-optimal experimental design, indirect modeling, and polynomial-modified Gaussian model, while grid point search method was selected for the final choice of the experimental conditions which lead to the best possible stationary phase performance for basic solutes. Good agreement between experimentally obtained chromatogram and simulated chromatogram for chosen experimental conditions (25% acetonitrile, 75 mM of NaClO₄, pH 4.00 on Zorbax Eclipse XDB column) confirmed the applicability of the proposed procedure. The additional point was selected for the verification of proposed procedure ability to distinguish changes in solutes' elution order. Simulated chromatogram for 21.5% acetonitrile, 85 mM of NaClO₄, pH 5.00 on Zorbax Eclipse XDB column was in line with experimental data. Furthermore, the values of left and right peak half-widths obtained from indirect modeling were used in order to evaluate performances of differently modified stationary phases applying a half-width plots approach. The results from half-width plot approach as well as from the proposed procedure indicate higher efficiency and better separation performance of the stationary phase extra densely bonded and double end-capped with trimethylsilyl group than the stationary phase with the combination of end-capping and bidentate silane bonding for chromatographic analysis of basic solutes in RP-HPLC systems with chaotropic agents. Graphical abstract ᅟ.