QSAR/QSPR Modeling: Introduction

Temperature-Dependent Density and Viscosity Prediction for Hydrocarbons: Machine Learning and Molecular Dynamics Simulations

Machine learning-based predictive models allow rapid and reliable prediction of material properties and facilitate innovative materials design. Base oils used in the formulation of lubricant products are complex hydrocarbons of varying sizes and structure. This study developed Gaussian process regression-based models to accurately predict the temperature-dependent density and dynamic viscosity of 305 complex hydrocarbons. In our approach, strongly correlated/collinear predictors were trimmed, important predictors were selected by least absolute shrinkage and selection operator (LASSO) regularization and prior domain knowledge, hyperparameters were systematically optimized by Bayesian optimization, and the models were interpreted. The approach provided versatile and quantitative structure-property relationship (QSPR) models with relatively simple predictors for determining the dynamic viscosity and density of complex hydrocarbons at any temperature. In addition, we developed molecular dynamics simulation-based descriptors and evaluated the feasibility and versatility of dynamic descriptors from simulations for predicting the material properties. It was found that the models developed using a comparably smaller pool of dynamic descriptors performed similarly in predicting density and viscosity to models based on many more static descriptors. The best models were shown to predict density and dynamic viscosity with coefficient of determination (R2) values of 99.6% and 97.7%, respectively, for all data sets, including a test data set of 45 molecules. Finally, partial dependency plots (PDPs), individual conditional expectation (ICE) plots, local interpretable model-agnostic explanation (LIME) values, and trimmed model R2 values were used to identify the most important static and dynamic predictors of the density and viscosity.

MolPredictX: Online Biological Activity Predictions by Machine Learning Models

Here we report the development of MolPredictX, an innovate and freely accessible web interface for biological activity predictions of query molecules. MolPredictX utilizes in-house QSAR models to provide 27 qualitative predictions (active or inactive), and quantitative probabilities for bioactivity against parasitic (Trypanosoma and Leishmania), viral (Dengue, Sars-CoV and Hepatitis C), pathogenic yeast (Candida albicans), bacterial (Salmonella enterica and Escherichia coli), and Alzheimer disease enzymes. In this article, we introduce the methodology and usability of this webtool, highlighting its potential role in the development of new drugs against a variety of diseases. MolPredictX is undergoing continuous development and is freely available at https://www.molpredictx.ufpb.br/.

Methodology-Centered Review of Molecular Modeling, Simulation, and Prediction of SARS-CoV-2

Despite tremendous efforts in the past two years, our understanding of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), virus-host interactions, immune response, virulence, transmission, and evolution is still very limited. This limitation calls for further in-depth investigation. Computational studies have become an indispensable component in combating coronavirus disease 2019 (COVID-19) due to their low cost, their efficiency, and the fact that they are free from safety and ethical constraints. Additionally, the mechanism that governs the global evolution and transmission of SARS-CoV-2 cannot be revealed from individual experiments and was discovered by integrating genotyping of massive viral sequences, biophysical modeling of protein-protein interactions, deep mutational data, deep learning, and advanced mathematics. There exists a tsunami of literature on the molecular modeling, simulations, and predictions of SARS-CoV-2 and related developments of drugs, vaccines, antibodies, and diagnostics. To provide readers with a quick update about this literature, we present a comprehensive and systematic methodology-centered review. Aspects such as molecular biophysics, bioinformatics, cheminformatics, machine learning, and mathematics are discussed. This review will be beneficial to researchers who are looking for ways to contribute to SARS-CoV-2 studies and those who are interested in the status of the field.

QSAR study of antituberculosis activity of oxadiazole derivatives using DFT calculations

Mycobacterium tuberculosis (Mtb) is the causative agent of infectious diseases worldwide. Oxadiazole derivatives have many biological activities and can be a good alternative to antimicrobial drugs. In this study, the quantitative structure-activity relationship (QSAR) of fifty-one novel oxadiazoles derivatives has been studied using the density functional theory (DFT) and statistical methods. Becke's three-parameter hybrid method and the Lee-Yang-Parr B3LYP functional employing 6-31++G (d) basis set are used to calculated quantum chemical descriptors using Gaussian09 software. The other descriptors including Lipinski, physicochemistry, topological, etc. were calculated using Chembio3d software. Statistically, the best correlation between the independent variables and the PMIC as the dependent variable was a 6-variable equation for which the correlation coefficient were as follows R² = 0.86 and R = 0.93. Also, the values of MAE = 0.003 and Q²_CV = 0.9 confirm the acceptability of the obtained model. The obtained equation shows that NRB, energy gap (ΔE), Henry's law constant, O-C, and C-N bonds length, and the Free Gibbs energy have the highest correlation with the anti-Tb activity.

Ecological hazard assessment via species sensitivity distributions: The non-exchangeability issue

Probabilistic approaches to hazard assessment use species sensitivity distributions (SSDs) to characterize hazard for toxicants exposure for different species within a community. Many of the assumptions at the core of SSDs are unrealistic, among them the assumption that the tolerance levels of all species in a specific ecological community are a priori exchangeable for each new toxic substance. Here we propose the use of a particular test to detect situations where such an assumption is violated. Then, a new method based on non-nested random effects model is required to identify novel SSDs capable of taking into account species non-exchangeability. Credible intervals, representing SSD uncertainty, could be determined based on our procedure. This leads to new and reliable estimates of the environmental hazard. We present a Bayesian modeling approach to address model inference issues, using Markov chain Monte Carlo sampling.

In silico screening of FDA approved drugs on AXL kinase and validation for breast cancer cell line

AXL kinase has been over expressed in many tumors and its involvement in cell proliferation, migration, survival, and resistance makes the kinase as attractive therapeutic target for many cancers. In this study, we performed a virtual screening of the food and drug administration (FDA) approved drug molecule database against AXL kinase for repurposing studies of breast cancer. We have identified three non-cancer drugs with good binding energies were subjected to in vitro breast cancer MCF-7 cell lines. Three drug molecules showing the activity with good IC50 values toward the cancer cell line. We also carried out a 2 dimensional (2 D) quantitative structure activity relation (QSAR) studies on N-[4-(Quinolin-4-yloxy)phenyl]benzenesulfonamides derivatives to design potent inhibitors for AXL kinase. The final QSAR equation was robust with good predictivity and the statistical validation having R² and Q² values are 0.91 and 0.86, respectively. QSAR equation descriptors informs about the chemical properties of AXL inhibitors and helpful for designing novel inhibitors. Communicated by Ramaswamy H. Sarma.

Probing the opportunities for designing anthelmintic leads by sub-structural topology-based QSAR modelling

A quantitative structure-activity (QSAR) model has been developed for enriched tubulin inhibitors, which were retrieved from sequence similarity searches and applicability domain analysis. Using partial least square (PLS) method and leave-one-out (LOO) validation approach, the model was generated with the correlation statistics of [Formula: see text] and [Formula: see text] of 0.68 and 0.69, respectively. The present study indicates that topological descriptors, viz. BIC, CH_3_C, IC, JX and Kappa_2 correlate well with biological activity. ADME and toxicity (or ADME/T) assessment showed that out of 260 molecules, 255 molecules successfully passed the ADME/T assessment test, wherein the drug-likeness attributes were exhibited. These results showed that topological indices and the colchicine binding domain directly influence the aetiology of helminthic infections. Further, we anticipate that our model can be applied for guiding and designing potential anthelmintic inhibitors.

QSAR modeling in ecotoxicological risk assessment: application to the prediction of acute contact toxicity of pesticides on bees (Apis mellifera L.)

Despite their indisputable importance around the world, the pesticides can be dangerous for a range of species of ecological importance such as honeybees (Apis mellifera L.). Thus, a particular attention should be paid to their protection, not only for their ecological importance by contributing to the maintenance of wild plant diversity, but also for their economic value as honey producers and crop-pollinating agents. For all these reasons, the environmental protection requires the resort of risk assessment of pesticides. The goal of this work was therefore to develop a validated QSAR model to predict contact acute toxicity (LD₅₀) of 111 pesticides to bees because the QSAR models devoted to this species are very scarce. The analysis of the statistical parameters of this model and those published in the literature shows that our model is more efficient. The QSAR model was assessed according to the OECD principles for the validation of QSAR models. The calculated values for the internal and external validation statistic parameters (Q ² and [Formula: see text] are greater than 0.85. In addition to this validation, a mathematical equation derived from the ANN model was used to predict the LD₅₀ of 20 other pesticides. A good correlation between predicted and experimental values was found (R ² = 0.97 and RMSE = 0.14). As a result, this equation could be a means of predicting the toxicity of new pesticides.

Ezqsar: An R Package for Developing QSAR Models Directly From Structures

Background:

Quantitative Structure Activity Relationship (QSAR) is a difficult computational chemistry approach for beginner scientists and a time consuming one for even more experienced researchers.

Method and Materials:

Ezqsar which is introduced here addresses both the issues. It considers important steps to have a reliable QSAR model. Besides calculation of descriptors using CDK library, highly correlated descriptors are removed, a provided data set is divided to train and test sets, descriptors are selected by a statistical method, statistical parameter for the model are presented and applicability domain is investigated.

Results:

Finally, the model can be applied to predict the activities for an extra set of molecules for a purpose of either lead optimization or virtual screening. The performance is demonstrated by an example.

Conclusion:

The R package, ezqsar, is freely available viahttps://github.com/shamsaraj/ezqsar, and it runs on Linux and MS-Windows.