QSAR study of α1β4 integrin inhibitors by GA-MLR and GA-SVM methods

QSAR analysis of VEGFR-2 inhibitors based on machine learning, Topomer CoMFA and molecule docking

VEGFR-2 kinase inhibitors are clinically approved drugs that can effectively target cancer angiogenesis. However, such inhibitors have adverse effects such as skin toxicity, gastrointestinal reactions and hepatic impairment. In this study, machine learning and Topomer CoMFA, which is an alignment-dependent, descriptor-based method, were employed to build structural activity relationship models of potentially new VEGFR-2 inhibitors. The prediction ac-curacy of the training and test sets of the 2D-SAR model were 82.4 and 80.1%, respectively, with KNN. Topomer CoMFA approach was then used for 3D-QSAR modeling of VEGFR-2 inhibitors. The coefficient of q2 for cross-validation of the model 1 was greater than 0.5, suggesting that a stable drug activity-prediction model was obtained. Molecular docking was further performed to simulate the interactions between the five most promising compounds and VEGFR-2 target protein and the Total Scores were all greater than 6, indicating that they had a strong hydrogen bond interactions were present. This study successfully used machine learning to obtain five potentially novel VEGFR-2 inhibitors to increase our arsenal of drugs to combat cancer.

A study of the interaction space of two lactate dehydrogenase isoforms (LDHA and LDHB) and some of their inhibitors using proteochemometrics modeling

Lactate dehydrogenase (LDH) is a tetramer enzyme that converts pyruvate to lactate reversibly. This enzyme becomes important because it is associated with diseases such as cancers, heart disease, liver problems, and most importantly, corona disease. As a system-based method, proteochemometrics does not require knowledge of the protein's three-dimensional structure, but rather depends on the amino acid sequence and protein descriptors. Here, we applied this methodology to model a set of LDHA and LDHB isoenzyme inhibitors. To implement the proteochemetrics method, the camb package in the R Studio Server programming environment was used. The activity of 312 compounds of LDHA and LDHB isoenzyme inhibitors from the valid Binding DB database was retrieved. The proteochemometrics method was applied to three machine learning algorithms gradient amplification model, random forest, and support vector machine as regression methods to find the best model. Through the combination of different models into an ensemble (greedy and stacking optimization), we explored the possibility of improving the performance of models. For the RF best ensemble model of inhibitors of LDHA and LDHB isoenzymes, and were 0.66 and 0.62, respectively. LDH inhibitory activation is influenced by Morgan fingerprints and topological structure descriptors.

Molecular Modeling Insights into the Structure and Behavior of Integrins: A Review

Integrins are heterodimeric glycoproteins crucial to the physiology and pathology of many biological functions. As adhesion molecules, they mediate immune cell trafficking, migration, and immunological synapse formation during inflammation and cancer. The recognition of the vital roles of integrins in various diseases revealed their therapeutic potential. Despite the great effort in the last thirty years, up to now, only seven integrin-based drugs have entered the market. Recent progress in deciphering integrin functions, signaling, and interactions with ligands, along with advancement in rational drug design strategies, provide an opportunity to exploit their therapeutic potential and discover novel agents. This review will discuss the molecular modeling methods used in determining integrins' dynamic properties and in providing information toward understanding their properties and function at the atomic level. Then, we will survey the relevant contributions and the current understanding of integrin structure, activation, the binding of essential ligands, and the role of molecular modeling methods in the rational design of antagonists. We will emphasize the role played by molecular modeling methods in progress in these areas and the designing of integrin antagonists.

Density function theory calculation to study the oxidation potential of electron-donating compounds; affirming the oxidation mechanism by NICS calculations

The manuscript describes a method for understanding the correlation of structural features and first oxidation potentials [Formula: see text] of electron-donating compounds (EDCs) with tetrathiafulvalene (TTF), dithiadiazafulvalenes (DTDAF), and tetraazafulvalene (TAF) frameworks. The density functional theory (DFT) procedure at B3LYP (6-31 + g(d)) was used for geometric optimization, given the large dimensions of the molecules studied, and their high structural similarity. First of all, the correlation between the oxidation potential and the highest occupied molecular orbital (HOMO) energy level as an effective quantum chemical descriptor was examined. Then, nucleus-independent chemical shifts (NICSs) calculation was applied to affirm the oxidation mechanism and interpret the effect of replacing the sulfur atoms by nitrogen, on the oxidation process. Finally, a more comprehensive investigation of structural features that affect the oxidation potential, topological, geometrical, constitutional, as well as, electrostatic, charged partial surface area, quantum-chemical, molecular orbital, and thermodynamic descriptors was calculated. A predictive model was developed based on the genetic algorithm multivariate linear regression (GA-MLR). There was an outstanding agreement between the theoretical and the experimental values obtained for the first oxidation potentials of the test set (Q²_Ext = 0.981).

Insight into the structural requirements of gelatinases (MMP-2 and MMP-9) inhibitors by multiple validated molecular modelling approaches: Part II

Inhibition of the matrix metalloproteinases (MMPs) is effective against metastasis of secondary tumours. Previous MMP inhibitors have failed in clinical trials due to their off-target toxicity in solid tumours. Thus, newer MMP inhibitors now have paramount importance. Here, different molecular modelling techniques were applied on a dataset of 110 gelatinase (MMP-2 and MMP-9) inhibitors. The objectives of the present study were to identify structural fingerprints for gelatinase inhibition and also to develop statistically validated QSAR models for the screening and prediction of different derivatives as MMP-2 (gelatinase A) and MMP-9 (gelatinase B) inhibitors. The Bayesian classification study provided the ROC values for the training set of 0.837 and 0.815 for MMP-2 and MMP-9, respectively. The linear model also produced the leave-one-out cross-validated Q² of 0.805 (eq. 1, MMP-2) and 0.724 (eq. 2, MMP-9), an r² of 0.845 (eq. 1, MMP-2) and 0.782 (eq. 2, MMP-9), an r²_Pred of 0.806 (eq. 1, MMP-2) and 0.732 (eq. 2, MMP-9). Similarly, non-linear learning models were also statistically significant and reliable. Overall, this study may help in the rational design of newer compounds with higher gelatinase inhibition to fight against both primary and secondary cancers in future.

The application of molecular topology for ulcerative colitis drug discovery

INTRODUCTION

Although the therapeutic arsenal against ulcerative colitis has greatly expanded (including the revolutionary advent of biologics), there remain patients who are refractory to current medications while the safety of the available therapeutics could also be improved. Molecular topology provides a theoretic framework for the discovery of new therapeutic agents in a very efficient manner, and its applications in the field of ulcerative colitis have slowly begun to flourish. Areas covered: After discussing the basics of molecular topology, the authors review QSAR models focusing on validated targets for the treatment of ulcerative colitis, entirely or partially based on topological descriptors. Expert opinion: The application of molecular topology to ulcerative colitis drug discovery is still very limited, and many of the existing reports seem to be strictly theoretic, with no experimental validation or practical applications. Interestingly, mechanism-independent models based on phenotypic responses have recently been reported. Such models are in agreement with the recent interest raised by network pharmacology as a potential solution for complex disorders. These and other similar studies applying molecular topology suggest that some therapeutic categories may present a 'topological pattern' that goes beyond a specific mechanism of action.

3D-QSAR and molecular docking study of LRRK2 kinase inhibitors by CoMFA and CoMSIA methods

Three-dimensional quantitative structure-activity relationship (3D-QSAR) modelling was conducted on a series of leucine-rich repeat kinase 2 (LRRK2) antagonists using CoMFA and CoMSIA methods. The data set, which consisted of 37 molecules, was divided into training and test subsets by using a hierarchical clustering method. Both CoMFA and CoMSIA models were derived using a training set on the basis of the common substructure-based alignment. The optimum PLS model built by CoMFA and CoMSIA provided satisfactory statistical results (q(2) = 0.589 and r(2) = 0.927 and q(2) = 0.473 and r(2) = 0.802, respectively). The external predictive ability of the models was evaluated by using seven compounds. Moreover, an external evaluation set with known experimental data was used to evaluate the external predictive ability of the porposed models. The statistical parameters indicated that CoMFA (after region focusing) has high predictive ability in comparison with standard CoMFA and CoMSIA models. Molecular docking was also performed on the most active compound to investigate the existence of interactions between the most active inhibitor and the LRRK2 receptor. Based on the obtained results and CoMFA contour maps, some features were introduced to provide useful insights for designing novel and potent LRRK2 inhibitors.