A Combinational Strategy of Model Disturbance and Outlier Comparison to Define Applicability Domain in Quantitative Structural Activity Relationship

Unveiling G-protein coupled receptor kinase-5 inhibitors for chronic degenerative diseases: Multilayered prioritization employing explainable machine learning-driven multi-class QSAR, ligand-based pharmacophore and free energy-inspired molecular simulation

GRK5 holds a pivotal role in cellular signaling pathways, with its overexpression in cardiomyocytes, neuronal cells, and tumor cells strongly associated with various chronic degenerative diseases, which highlights the urgent need for potential inhibitors. In this study, multiclass classification-based QSAR models were developed using diverse machine learning algorithms. These models were built from curated compounds with experimentally derived GRK5 inhibitory activity. Additionally, a pharmacophore model was constructed using active compounds from the dataset. Among the models, the SVM-based approach proved most effective and was initially used to screen DrugBank compounds within the applicability domain. Compounds showing significant GRK5 inhibitory potential underwent evaluation for key pharmacophoric features. Prospective compounds were subjected to molecular docking to assess binding affinity towards GRK5's key active site amino acid residues. Stability at the binding site was analyzed through 200 ns molecular dynamics simulations. MM-GBSA analysis quantified individual free energy components contributing to the total binding energy with respect to binding site residues. Metadynamics analysis, including PCA, FEL, and PDF, provided crucial insights into conformational changes of both apo and holo forms of GRK5 at defined energy states. The study identifies DB02844 (S-Adenosyl-1,8-Diamino-3-Thiooctane) and DB13155 (Esculin) as promising GRK5 inhibitors, warranting further in vitro and in vivo validation studies.

Machine learning and traditional QSAR modeling methods: a case study of known PXR activators

Pregnane X receptor (PXR), extensively expressed in human tissues related to digestion and metabolism, is responsible for recognizing and detoxifying diverse xenobiotics encountered by humans. To comprehend the promiscuous nature of PXR and its ability to bind a variety of ligands, computational approaches, viz., quantitative structure-activity relationship (QSAR) models, aid in the rapid dereplication of potential toxicological agents and mitigate the number of animals used to establish a meaningful regulatory decision. Recent advancements in machine learning techniques accommodating larger datasets are expected to aid in developing effective predictive models for complex mixtures (viz., dietary supplements) before undertaking in-depth experiments. Five hundred structurally diverse PXR ligands were used to develop traditional two-dimensional (2D) QSAR, machine-learning-based 2D-QSAR, field-based three-dimensional (3D) QSAR, and machine-learning-based 3D-QSAR models to establish the utility of predictive machine learning methods. Additionally, the applicability domain of the agonists was established to ensure the generation of robust QSAR models. A prediction set of dietary PXR agonists was used to externally-validate generated QSAR models. QSAR data analysis revealed that machine-learning 3D-QSAR techniques were more accurate in predicting the activity of external terpenes with an external validation squared correlation coefficient (R2) of 0.70 versus an R2 of 0.52 in machine-learning 2D-QSAR. Additionally, a visual summary of the binding pocket of PXR was assembled from the field 3D-QSAR models. By developing multiple QSAR models in this study, a robust groundwork for assessing PXR agonism from various chemical backbones has been established in anticipation of the identification of potential causative agents in complex mixtures.

Exploring proteasome inhibition using atomic weighted vector indices and machine learning approaches

Ubiquitin-proteasome system (UPS) is a highly regulated mechanism of intracellular protein degradation and turnover. The UPS is involved in different biological activities, such as the regulation of gene transcription and cell cycle. Several researchers have applied cheminformatics and artificial intelligence methods to study the inhibition of proteasomes, including the prediction of UPP inhibitors. Following this idea, we applied a new tool for obtaining molecular descriptors (MDs) for modeling proteasome Inhibition in terms of EC₅₀ (µmol/L), in which a set of new MDs called atomic weighted vectors (AWV) and several prediction algorithms were used in cheminformatics studies. In the manuscript, a set of descriptors based on AWV are presented as datasets for training different machine learning techniques, such as linear regression, multiple linear regression (MLR), random forest (RF), K-nearest neighbors (IBK), multi-layer perceptron, best-first search, and genetic algorithm. The results suggest that these atomic descriptors allow adequate modeling of proteasome inhibitors despite artificial intelligence techniques, as a variant to build efficient models for the prediction of inhibitory activity.

Discovery of novel anti-cyanobacterial allelochemicals by multi-conformational QSAR approach

Microcystis aeruginosa causes cyanobacterial harmful algal blooms (cHABs) in various freshwater environments. Due to global climate change, the cHABs have even spread to estuaries and coasts. Plant-derived flavones have been reported as allelochemicals that efficiently inhibit the growth of M. aeruginosa. Quantitative structure-activity relationship (QSAR) was applied to investigate the factors affecting the M. aeruginosa inhibitory activity of flavones, and to discover novel allelochemicals against M. aeruginosa. We constructed 2D and 3D-QSAR models based on the half maximum inhibitory concentration (IC₅₀) of 22 flavones against M. aeruginosa, using molecular descriptors from multiple stable conformations. Both models showed satisfactory performances (2D-QSAR: r²=0.899, q²=0.596, r_test²=0.801; 3D-QSAR: r²=0.810, q²=0.516, r_test²=0.897). The 2D-QSAR model indicates that the anti-cyanobacterial activity is positively correlated with minimum and maximum surface electrostatic potential, and negatively correlated with polarity index and polar surface area. Through the 3D-QSAR approach, electronegative hydroxyl groups in 5- and 4^'-position were favorable for the anti-cyanobacterial activity. In addition, we selected six untested flavones that fit the "activity-favorable" pattern of the visualized 3D-QSAR model. Five of the external flavones exhibited significant cyanobacterial inhibitory ability at their predicted IC₅₀ by the 3D-QSAR model. In particular, diosmetin achieved an inhibition rate of 70.50±4.74%, which was much higher than expected. The flavones screened by the 3D-QSAR model are novel cyanobacterial inhibitors and should be fully exploited to mitigate cHABs.

Unsupervised machine learning, QSAR modelling and web tool development for streamlining the lead identification process of antimalarial flavonoids

Identification of lead compounds with the traditional laboratory approach is expensive and time-consuming. Nowadays, in silico techniques have emerged as a promising approach for lead identification. In this study, we aim to develop robust and predictive 2D-QSAR models to identify lead flavonoids by predicting the IC₅₀ against Plasmodium falciparum. We applied machine learning algorithms (Principal component analysis followed by K-means clustering) and Pearson correlation analysis to select 9 molecular descriptors (MDs) for model building. We selected and validated the three best QSAR models after execution of multiple linear regression (MLR) 100 times with different combinations of MDs. The developed models have fulfilled the five principles for QSAR models as specified by the Organization for Economic Co-operation and Development. The outcome of the study is a reliable and sustainable in silico method of IC₅₀ (Mean ± SD) prediction that will positively impact the antimalarial drug development process by reducing the money and time required to identify potential antimalarial lead compounds from the class of flavonoids. We also developed a web tool (JazQSAR, https://etflin.com/news/4) to offer an easily accessible platform for the developed QSAR models.

Nano-QSAR modeling for predicting the cytotoxicity of metallic and metal oxide nanoparticles: A review

Given the rapid development of nanotechnology, it is crucial to understand the effects of nanoparticles on living organisms. However, it is laborious to perform toxicological tests on a case-by-case basis. Quantitative structure-activity relationship (QSAR) is an effective computational technique because it saves time, costs, and animal sacrifice. Therefore, this review presents general procedures for the construction and application of nano-QSAR models of metal-based and metal-oxide nanoparticles (MBNPs and MONPs). We also provide an overview of available databases and common algorithms. The molecular descriptors and their roles in the toxicological interpretation of MBNPs and MONPs are systematically reviewed and the future of nano-QSAR is discussed. Finally, we address the growing demand for novel nano-specific descriptors, new computational strategies to address the data shortage, in situ data for regulatory concerns, a better understanding of the physicochemical properties of NPs with bioactivity, and, most importantly, the design of nano-QSAR for real-life environmental predictions rather than laboratory simulations.

An alignment-independent three-dimensional quantitative structure-activity relationship study on ron receptor tyrosine kinase inhibitors

Recepteur d'Origine Nantais known as RON is a member of the receptor tyrosine kinase (RTK) superfamily which has recently gained increasing attention as cancer target for therapeutic intervention. The aim of this work was to perform an alignment-independent three-dimensional quantitative structure-activity relationship (3D QSAR) study for a series of RON inhibitors. A 3D QSAR model based on GRid-INdependent Descriptors (GRIND) methodology was generated using a set of 19 compounds with RON inhibitory activities. The generated 3D QSAR model revealed the main structural features important in the potency of RON inhibitors. The results obtained from the presented study can be used in lead optimization projects for designing of novel compounds where inhibition of RON is needed.

Planning new Trypanosoma cruzi CYP51 inhibitors using QSAR studies

Chagas disease kills over 10,000 people per year, and approximately 8 million people are infected by Trypanosoma cruzi. The reference drug for treatment of the disease, benznidazole, is the same since the 70s. In recent years, many CYP51 inhibitors were tested against this parasite's target. One of them, posaconazole, was even tested in clinical trials that unfortunately were not successful. Nevertheless, there are still many evidences that CYP51 is a great potential target to treat T. cruzi infection. The research for new effective molecules that can cure the chronic phase of the disease is essential. 2D and 3D-quantitative structure activity relationship (QSAR) studies were conducted in this work to create three QSAR models using the chemical structures of 197 published compounds that already went through either in vivo or in vitro tests. After the analysis of the models, new analogues not yet synthesized were suggested here and had their biological activity and synthetic availability assessed.

Molecular modeling of a series of dehydroquinate dehydratase type II inhibitors of Mycobacterium tuberculosis and design of new binders

Tuberculosis, caused by Mycobacterium tuberculosis (M. tuberculosis), is still responsible for a large number of fatal cases, especially in developing countries with alarming rates of incidence and prevalence worldwide. Mycobacterium tuberculosis has a remarkable ability to develop new resistance mechanisms to the conventional antimicrobials treatment. Because of this, there is an urgent need for novel bioactive compounds for its treatment. The dehydroquinate dehydratase II (DHQase II) is considered a key enzyme of shikimate pathway, and it can be used as a promising target for the design of new bioactive compounds with antibacterial action. The aim of this work was the construction of QSAR models to aid the design of new potential DHQase II inhibitors. For that purpose, various molecular modeling approaches, such as activity cliff, QSAR models and computer-aided ligand design were utilized. A predictive in silico 4D-QSAR model was built using a database comprising 86 inhibitors of DHQase II, and the model was used to predict the activity of the designed ligands. The obtained model proved to predict well the DHQase II inhibition for an external validation dataset ([Formula: see text] = 0.72). Also, the Activity Cliff analysis shed light on important structural features applied to the ligand design.

HIVprotI: an integrated web based platform for prediction and design of HIV proteins inhibitors

A number of anti-retroviral drugs are being used for treating Human Immunodeficiency Virus (HIV) infection. Due to emergence of drug resistant strains, there is a constant quest to discover more effective anti-HIV compounds. In this endeavor, computational tools have proven useful in accelerating drug discovery. Although methods were published to design a class of compounds against a specific HIV protein, but an integrated web server for the same is lacking. Therefore, we have developed support vector machine based regression models using experimentally validated data from ChEMBL repository. Quantitative structure activity relationship based features were selected for predicting inhibition activity of a compound against HIV proteins namely protease (PR), reverse transcriptase (RT) and integrase (IN). The models presented a maximum Pearson correlation coefficient of 0.78, 0.76, 0.74 and 0.76, 0.68, 0.72 during tenfold cross-validation on IC₅₀ and percent inhibition datasets of PR, RT, IN respectively. These models performed equally well on the independent datasets. Chemical space mapping, applicability domain analyses and other statistical tests further support robustness of the predictive models. Currently, we have identified a number of chemical descriptors that are imperative in predicting the compound inhibition potential. HIVprotI platform (http://bioinfo.imtech.res.in/manojk/hivproti) would be useful in virtual screening of inhibitors as well as designing of new molecules against the important HIV proteins for therapeutics development.

Hyperspectral imaging for mapping of total nitrogen spatial distribution in pepper plant

Visible/near-infrared (Vis/NIR) hyperspectral imaging was employed to determine the spatial distribution of total nitrogen in pepper plant. Hyperspectral images of samples (leaves, stems, and roots of pepper plants) were acquired and their total nitrogen contents (TNCs) were measured using Dumas combustion method. Mean spectra of all samples were extracted from regions of interest (ROIs) in hyperspectral images. Random frog (RF) algorithm was implemented to select important wavelengths which carried effective information for predicting the TNCs in leaf, stem, root, and whole-plant (leaf-stem-root), respectively. Based on full spectra and the selected important wavelengths, the quantitative relationships between spectral data and the corresponding TNCs in organs (leaf, stem, and root) and whole-plant (leaf-stem-root) were separately developed using partial least-squares regression (PLSR). As a result, the PLSR model built by the important wavelengths for predicting TNCs in whole-plant (leaf-stem-root) offered a promising result of correlation coefficient (R) for prediction (RP = 0.876) and root mean square error (RMSE) for prediction (RMSEP = 0.426%). Finally, the TNC of each pixel within ROI of the sample was estimated to generate the spatial distribution map of TNC in pepper plant. The achievements of the research indicated that hyperspectral imaging is promising and presents a powerful potential to determine nitrogen contents spatial distribution in pepper plant.