Optimal Piecewise Linear Regression Algorithm for QSAR Modelling

A machine learning-based QSAR model reveals important molecular features for understanding the potential inhibition mechanism of ionic liquids to acetylcholinesterase

The broad application of ionic liquids (ILs) has been hindered by uncertainties surrounding their ecotoxicity. In this work, a Quantitative Structure-Activity Relationship (QSAR) model was devised to predict the inhibition of ILs towards the activity of AChE, employing both Random Forest (RF) and eXtreme Gradient Boosting (XGBoost) machine learning approaches. Fourteen kings of essential molecular feature descriptors were screened from an initial roster of 244 descriptors through the application of a feature importance index and they showed a significant impact on the activity of AChE activity. The two models based solely on the 14 most critical molecular descriptors could maintain model's robustness and reliability. The correlation analysis between these 14 descriptors and the inhibition of AChE activity revealed the potential impact of the molecular characteristics on ILs toxicity. The results underscored the main influence of cations in ILs on the inhibitory activity towards the AChE enzyme. Specifically, cations exhibiting hydrophobicity properties were found to exert more potent inhibitory effects on the AChE enzyme. In addition, some other properties of the cations, such as the degree of branching, atomic weight and partial charge also modulated their inhibition potential. This study enhances the comprehension of the structure-activity relationship between ILs and AChE inhibition, providing a reference for designing safer and greener ILs.

Optimisation-based modelling for explainable lead discovery in malaria

BACKGROUND

The search for new antimalarial treatments is urgent due to growing resistance to existing therapies. The Open Source Malaria (OSM) project offers a promising starting point, having extensively screened various compounds for their effectiveness. Further analysis of the chemical space surrounding these compounds could provide the means for innovative drugs.

METHODS

We report an optimisation-based method for quantitative structure-activity relationship (QSAR) modelling that provides explainable modelling of ligand activity through a mathematical programming formulation. The methodology is based on piecewise regression principles and offers optimal detection of breakpoint features, efficient allocation of samples into distinct sub-groups based on breakpoint feature values, and insightful regression coefficients. Analysis of OSM antimalarial compounds yields interpretable results through rules generated by the model that reflect the contribution of individual fingerprint fragments in ligand activity prediction. Using knowledge of fragment prioritisation and screening of commercially available compound libraries, potential lead compounds for antimalarials are identified and evaluated experimentally via a Plasmodium falciparum asexual growth inhibition assay (PfGIA) and a human cell cytotoxicity assay.

CONCLUSIONS

Three compounds are identified as potential leads for antimalarials using the methodology described above. This work illustrates how explainable predictive models based on mathematical optimisation can pave the way towards more efficient fragment-based lead discovery as applied in malaria.

Improving predictions and understanding of primary and ultimate biodegradation rates with machine learning models

This study aimed to develop machine learning based quantitative structure biodegradability relationship (QSBR) models for predicting primary and ultimate biodegradation rates of organic chemicals, which are essential parameters for environmental risk assessment. For this purpose, experimental primary and ultimate biodegradation rates of high consistency were compiled for 173 organic compounds. A significant number of descriptors were calculated with a collection of quantum/computational chemistry software and tools to achieve comprehensive representation and interpretability. Following a pre-screening process, multiple QSBR models were developed for both primary and ultimate endpoints using three algorithms: extreme gradient boosting (XGBoost), support vector machine (SVM), and multiple linear regression (MLR). Furthermore, a unified QSBR model was constructed using the knowledge transfer technique and XGBoost. Results demonstrated that all QSBR models developed in this study had good performance. Particularly, SVM models exhibited high level of goodness of fit (coefficient of determination on the training set of 0.973 for primary and 0.980 for ultimate), robustness (leave-one-out cross-validated coefficient of 0.953 for primary and 0.967 for ultimate), and external predictive ability (external explained variance of 0.947 for primary and 0.958 for ultimate). The knowledge transfer technique enhanced model performance by learning from properties of two biodegradation endpoints. Williams plots were used to visualize the application domains of the models. Through SHapley Additive exPlanations (SHAP) analysis, this study identified key features affecting biodegradation rates. Notably, MDEO-12, APC2D1_C_O, and other features contributed to primary biodegradation, while AATS0v, AATS2v, and others inhibited it. For ultimate biodegradation, features like No. of Rotatable Bonds, APC2D1_C_O, and minHBa were contributors, while C1SP3, Halogen Ratio, GGI4, and others hindered the process. Also, the study quantified the contributions of each feature in predictions for individual chemicals. This research provides valuable tools for predicting both primary and ultimate biodegradation rates while offering insights into the mechanisms.

Deimos: A novel automated methodology for optimal grouping. Application to nanoinformatics case studies

In this study we present deimos, a computational methodology for optimal grouping, applied on the read-across prediction of engineered nanomaterials' (ENMs) toxicity-related properties. The method is based on the formulation and the solution of a mixed-integer optimization program (MILP) problem that automatically and simultaneously performs feature selection, defines the grouping boundaries according to the response variable and develops linear regression models in each group. For each group/region, the characteristic centroid is defined in order to allocate untested ENMs to the groups. The deimos MILP problem is integrated in a broader optimization workflow that selects the best performing methodology between the standard multiple linear regression (MLR), the least absolute shrinkage and selection operator (LASSO) models and the proposed deimos multiple-region model. The performance of the suggested methodology is demonstrated through the application to benchmark ENMs datasets and comparison with other predictive modelling approaches. However, the proposed method can be applied to property prediction of other than ENM chemical entities and it is not limited to ENMs toxicity prediction.

The Estimation of Bent Line Expectile Regression Model Based on a Smoothing Technique

A bent line expectile regression model can describe the effect of a covariate on the response variable with two different straight lines, which intersect at an unknown change-point. Due to the existence of the change-point, the objective function of the model is not differentiable with respect to the change-point, so it cannot be solved by the method of the traditional linear expectile regression model. For this model, a new estimation method is proposed by a smoothing technique, that is, using Gaussian kernel function to approximate the indicator function in the objective function. It can not only estimate the regression coefficients and change-point location simultaneously, but also have better estimation effect, which compensates for the insufficiency of the previous estimation methods. Under the given regularity conditions, the theoretical proofs of the consistency and asymptotic normality of the proposed estimators are derived. There are two parts of numerical simulations in this paper. Simulation 1 discusses various error distributions at different expectile levels under different conditions, the results show that the mean values of the biases of the estimation method in this paper, and other indicators, are very small, which indicates the robust property of the new method. Simulation 2 considers the symmetric and asymmetric bent lien expectile regression models, the results show that the estimated values of the estimation method in this paper are similar to the true values, which indicates the estimation effect and large sample performance of the proposed method are excellent. In the application research, the method in this paper is applied to the Arctic annual average temperature data and the Nile annual average flow data. The research shows that the standard errors of the estimation method in this paper are very similar to 0, indicating that the parameter estimation accuracy of the new method is very high, and the location of the change-point can be accurately estimated, which further confirms that the new method is effective and feasible.

In Silico Antiprotozoal Evaluation of 1,4-Naphthoquinone Derivatives against Chagas and Leishmaniasis Diseases Using QSAR, Molecular Docking, and ADME Approaches

Chagas and leishmaniasis are two neglected diseases considered as public health problems worldwide, for which there is no effective, low-cost, and low-toxicity treatment for the host. Naphthoquinones are ligands with redox properties involved in oxidative biological processes with a wide variety of activities, including antiparasitic. In this work, in silico methods of quantitative structure–activity relationship (QSAR), molecular docking, and calculation of ADME (absorption, distribution, metabolism, and excretion) properties were used to evaluate naphthoquinone derivatives with unknown antiprotozoal activity. QSAR models were developed for predicting antiparasitic activity against Trypanosoma cruzi, Leishmania amazonensis, and Leishmania infatum, as well as the QSAR model for toxicity activity. Most of the evaluated ligands presented high antiparasitic activity. According to the docking results, the family of triazole derivatives presented the best affinity with the different macromolecular targets. The ADME results showed that most of the evaluated compounds present adequate conditions to be administered orally. Naphthoquinone derivatives show good biological activity results, depending on the substituents attached to the quinone ring, and perhaps the potential to be converted into drugs or starting molecules.

Novel drug-target interactions via link prediction and network embedding

BACKGROUND

As many interactions between the chemical and genomic space remain undiscovered, computational methods able to identify potential drug-target interactions (DTIs) are employed to accelerate drug discovery and reduce the required cost. Predicting new DTIs can leverage drug repurposing by identifying new targets for approved drugs. However, developing an accurate computational framework that can efficiently incorporate chemical and genomic spaces remains extremely demanding. A key issue is that most DTI predictions suffer from the lack of experimentally validated negative interactions or limited availability of target 3D structures.

RESULTS

We report DT2Vec, a pipeline for DTI prediction based on graph embedding and gradient boosted tree classification. It maps drug-drug and protein-protein similarity networks to low-dimensional features and the DTI prediction is formulated as binary classification based on a strategy of concatenating the drug and target embedding vectors as input features. DT2Vec was compared with three top-performing graph similarity-based algorithms on a standard benchmark dataset and achieved competitive results. In order to explore credible novel DTIs, the model was applied to data from the ChEMBL repository that contain experimentally validated positive and negative interactions which yield a strong predictive model. Then, the developed model was applied to all possible unknown DTIs to predict new interactions. The applicability of DT2Vec as an effective method for drug repurposing is discussed through case studies and evaluation of some novel DTI predictions is undertaken using molecular docking.

CONCLUSIONS

The proposed method was able to integrate and map chemical and genomic space into low-dimensional dense vectors and showed promising results in predicting novel DTIs.

Application of multivariate linear regression models for selection of deep eutectic solvent for extraction of apigenin and luteolin from Chrysanthemum indicum L

INTRODUCTION

Among a variety of compounds presented in chrysanthemum, apigenin and luteolin are the two main components that play a major role in numerous biological activities of this herb.

OBJECTIVES

We aimed to obtain linear models showing the dependence of the yield of extraction of apigenin and luteolin on the composition of deep eutectic solvent and investigate the extraction of these two ingredients from Chrysanthemum indicum L.

METHODS

Two models showing the dependence of luteolin and apigenin concentrations on the composition of the solvent were established using a multilinear regression algorithm and were applied to screen 119 different solvents. After that, the extraction process was optimized using response surface methodology and an artificial neural network. Apigenin and luteolin were recovered from the extract by the combination of distillation and addition of water.

RESULTS

The screening results on 119 solvents revealed that choline chloride-acetic acid (1:4) was the most suitable deep eutectic solvent. It was showed that both response surface methodology and the artificial neural network could accurately determine the optimal conditions of extraction of apigenin and luteolin from C. indicum L., including time of extraction (65 minutes), temperature of extraction (90°C) and water content (20%). By the combination of distillation and addition of water, apigenin and luteolin could be effectively recovered from the deep eutectic solvent extract with a recovery rate of over 80%.

CONCLUSIONS

Deep eutectic solvent could be used as an effective green alternative to the conventional solvents for the extraction of bioactive compounds from plants.

Automatic Decision-Making Style Recognition Method Using Kinect Technology

In recent years, somatosensory interaction technology, represented by Microsoft's Kinect hardware platform, has been widely used in various fields, such as entertainment, education, and medicine. Kinect technology can easily capture and record behavioral data, which provides new opportunities for behavioral and psychological correlation analysis research. In this paper, an automatic decision-style recognition method is proposed. Experiments involving 240 subjects were conducted to obtain face data and individual decision-making style score. The face data was obtained using the Kinect camera, and the decision-style score were obtained via a questionnaire. To realize automatic recognition of an individual decision-making style, machine learning was employed to establish the mapping relationship between the face data and a scaled evaluation of the decision-making style score. This study adopts a variety of classical machine learning algorithms, including Linear regression, Support vector machine regression, Ridge regression, and Bayesian ridge regression. The experimental results show that the linear regression model returns the best results. The correlation coefficient between the linear regression model evaluation results and the scale evaluation results was 0.6, which represents a medium and higher correlation. The results verify the feasibility of automatic decision-making style recognition method based on facial analysis.

Deep Learning Approach for Discovery of In Silico Drugs for Combating COVID-19

Early diagnosis of pandemic diseases such as COVID-19 can prove beneficial in dealing with difficult situations and helping radiologists and other experts manage staffing more effectively. The application of deep learning techniques for genetics, microscopy, and drug discovery has created a global impact. It can enhance and speed up the process of medical research and development of vaccines, which is required for pandemics such as COVID-19. However, current drugs such as remdesivir and clinical trials of other chemical compounds have not shown many impressive results. Therefore, it can take more time to provide effective treatment or drugs. In this paper, a deep learning approach based on logistic regression, SVM, Random Forest, and QSAR modeling is suggested. QSAR modeling is done to find the drug targets with protein interaction along with the calculation of binding affinities. Then deep learning models were used for training the molecular descriptor dataset for the robust discovery of drugs and feature extraction for combating COVID-19. Results have shown more significant binding affinities (greater than −18) for many molecules that can be used to block the multiplication of SARS-CoV-2, responsible for COVID-19.

Automated Grouping of Nanomaterials and Read-Across Prediction of Their Adverse Effects Based on Mathematical Optimization

In this study, a computational workflow is presented for grouping engineered nanomaterials (ENMs) and for predicting their toxicity-related end points. A mixed integer-linear optimization program (MILP) problem is formulated, which automatically filters out the noisy variables, defines the grouping boundaries, and develops specific to each group predictive models. The method is extended to the multidimensional space, by considering the ENM characterization categories (e.g., biological, physicochemical, biokinetics, image etc.) as different dimensions. The performance of the proposed method is illustrated through the application to benchmark data sets and comparison with alternative predictive modeling approaches. The trained models using the above data sets were made publicly available through a user-friendly web service.

Prediction of K562 Cells Functional Inhibitors Based on Machine Learning Approaches

BACKGROUND

β thalassemia is a common monogenic genetic disease that is very harmful to human health. The disease arises is due to the deletion of or defects in β-globin, which reduces synthesis of the β-globin chain, resulting in a relatively excess number of α-chains. The formation of inclusion bodies deposited on the cell membrane causes a decrease in the ability of red blood cells to deform and a group of hereditary haemolytic diseases caused by massive destruction in the spleen.

METHODS

In this work, machine learning algorithms were employed to build a prediction model for inhibitors against K562 based on 117 inhibitors and 190 non-inhibitors.

RESULTS

The overall accuracy (ACC) of a 10-fold cross-validation test and an independent set test using Adaboost were 83.1% and 78.0%, respectively, surpassing Bayes Net, Random Forest, Random Tree, C4.5, SVM, KNN and Bagging.

CONCLUSION

This study indicated that Adaboost could be applied to build a learning model in the prediction of inhibitors against K526 cells.

Network-based piecewise linear regression for QSAR modelling

Quantitative Structure-Activity Relationship (QSAR) models are critical in various areas of drug discovery, for example in lead optimisation and virtual screening. Recently, the need for models that are not only predictive but also interpretable has been highlighted. In this paper, a new methodology is proposed to build interpretable QSAR models by combining elements of network analysis and piecewise linear regression. The algorithm presented, modSAR, splits data using a two-step procedure. First, compounds associated with a common target are represented as a network in terms of their structural similarity, revealing modules of similar chemical properties. Second, each module is subdivided into subsets (regions), each of which is modelled by an independent linear equation. Comparative analysis of QSAR models across five data sets of protein inhibitors obtained from ChEMBL is reported and it is shown that modSAR offers similar predictive accuracy to popular algorithms, such as Random Forest and Support Vector Machine. Moreover, we show that models built by modSAR are interpretatable, capable of evaluating the applicability domain of the compounds and serve well tasks such as virtual screening and the development of new drug leads.