QSPR approaches to elucidate the stability constants between β-cyclodextrin and some organic compounds: Docking based 3D conformer

Phenolic Acid-β-Cyclodextrin Complexation Study to Mask Bitterness in Wheat Bran: A Machine Learning-Based QSAR Study

The need to solvate and encapsulate hydro-sensitive molecules drives noticeable trends in the applications of cyclodextrins in the pharmaceutical industry, in foods, polymers, materials, and in agricultural science. Among them, β-cyclodextrin is one of the most used for the entrapment of phenolic acid compounds to mask the bitterness of wheat bran. In this regard, there is still a need for good data and especially for a robust predictive model that assesses the bitterness masking capabilities of β-cyclodextrin for various phenolic compounds. This study uses a dataset of 20 phenolic acids docked into the β-cyclodextrin cavity to generate three different binding constants. The data from the docking study were combined with topological, topographical, and quantum-chemical features from the ligands in a machine learning-based structure-activity relationship study. Three different models for each binding constant were computed using a combination of the genetic algorithm (GA) and multiple linear regression (MLR) approaches. The developed ML/QSAR models showed a very good performance, with high predictive ability and correlation coefficients of 0.969 and 0.984 for the training and test sets, respectively. The models revealed several factors responsible for binding with cyclodextrin, showing positive contributions toward the binding affinity values, including such features as the presence of six-membered rings in the molecule, branching, electronegativity values, and polar surface area.

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.

Development of QSRR model for hydroxamic acids using PCA-GA-BP algorithm incorporated with molecular interaction-based features

As a potent zinc chelator, hydroxamic acid has been applied in the design of inhibitors of zinc metalloenzyme, such as histone deacetylases (HDACs). A series of hydroxamic acids with HDAC inhibitory activities were subjected to the QSRR (Quantitative Structure-Retention Relationships) study. Experimental data in combination with calculated molecular descriptors were used for the development of the QSRR model. Specially, we employed PCA (principal component analysis) to accomplish dimension reduction of descriptors and utilized the principal components of compounds (16 training compounds, 4 validation compounds and 7 test compounds) to execute GA (genetic algorithm)-BP (error backpropagation) algorithm. We performed double cross-validation approach for obtaining a more convincing model. Moreover, we introduced molecular interaction-based features (molecular docking scores) as a new type of molecular descriptor to represent the interactions between analytes and the mobile phase. Our results indicated that the incorporation of molecular interaction-based features significantly improved the accuracy of the QSRR model, (R2 value is 0.842, RMSEP value is 0.440, and MAE value is 0.573). Our study not only developed QSRR model for the prediction of the retention time of hydroxamic acid in HPLC but also proved the feasibility of using molecular interaction-based features as molecular descriptors.

Predicting Value of Binding Constants of Organic Ligands to Beta-Cyclodextrin: Application of MARSplines and Descriptors Encoded in SMILES String

The quantitative structure–activity relationship (QSPR) model was formulated to quantify values of the binding constant (lnK) of a series of ligands to beta–cyclodextrin (β-CD). For this purpose, the multivariate adaptive regression splines (MARSplines) methodology was adopted with molecular descriptors derived from the simplified molecular input line entry specification (SMILES) strings. This approach allows discovery of regression equations consisting of new non-linear components (basis functions) being combinations of molecular descriptors. The model was subjected to the standard internal and external validation procedures, which indicated its high predictive power. The appearance of polarity-related descriptors, such as XlogP, confirms the hydrophobic nature of the cyclodextrin cavity. The model can be used for predicting the affinity of new ligands to β-CD. However, a non-standard application was also proposed for classification into Biopharmaceutical Classification System (BCS) drug types. It was found that a single parameter, which is the estimated value of lnK, is sufficient to distinguish highly permeable drugs (BCS class I and II) from low permeable ones (BCS class II and IV). In general, it was found that drugs of the former group exhibit higher affinity to β-CD then the latter group (class III and IV).

3D-QSAR for binding constants of β-cyclodextrin host-guest complexes by utilising spectrophores as molecular descriptors

The increasing use of cyclodextrins (CDs) as host-molecules for host-guest complexes and their remediation application in environmental science requires to establish easily accessible models of "quantitative structure-activity relationships" (QSARs) to predict the binding constants. A new open-source molecular descriptor, so called spectrophores, was utilised to build 3D-QSAR models which have R² of 0.95 and RMSE of 0.20.

Free-energy patterns in inclusion complexes: the relevance of non-included moieties in the stability constants

A MD/PMF-based procedure is designed for quantification of the interaction and respective components, guiding complex formation in water between β-CD and several naphthalene derivatives, highlighting the relevance of substituents.