Exploring 2D-QSAR for prediction of beta-secretase 1 (BACE1) inhibitory activity against Alzheimer's disease

Network Models of BACE-1 Inhibitors: Exploring Structural and Biochemical Relationships

This study investigates the clustering patterns of human β-secretase 1 (BACE-1) inhibitors using complex network methodologies based on various distance functions, including Euclidean, Tanimoto, Hamming, and Levenshtein distances. Molecular descriptor vectors such as molecular mass, Merck Molecular Force Field (MMFF) energy, Crippen partition coefficient (ClogP), Crippen molar refractivity (MR), eccentricity, Kappa indices, Synthetic Accessibility Score, Topological Polar Surface Area (TPSA), and 2D/3D autocorrelation entropies are employed to capture the diverse properties of these inhibitors. The Euclidean distance network demonstrates the most reliable clustering results, with strong agreement metrics and minimal information loss, indicating its robustness in capturing essential structural and physicochemical properties. Tanimoto and Hamming distance networks yield valuable clustering outcomes, albeit with moderate performance, while the Levenshtein distance network shows significant discrepancies. The analysis of eigenvector centrality across different networks identifies key inhibitors acting as hubs, which are likely critical in biochemical pathways. Community detection results highlight distinct clustering patterns, with well-defined communities providing insights into the functional and structural groupings of BACE-1 inhibitors. The study also conducts non-parametric tests, revealing significant differences in molecular descriptors, validating the clustering methodology. Despite its limitations, including reliance on specific descriptors and computational complexity, this study offers a comprehensive framework for understanding molecular interactions and guiding therapeutic interventions. Future research could integrate additional descriptors, advanced machine learning techniques, and dynamic network analysis to enhance clustering accuracy and applicability.

Multi-method computational evaluation of the inhibitors against leucine-rich repeat kinase 2 G2019S mutant for Parkinson's disease

Leucine-rich repeat kinase 2 G2019S mutant (LRRK2 G2019S) is a potential target for Parkinson's disease therapy. In this work, the computational evaluation of the LRRK2 G2019S inhibitors was conducted via a combined approach which contains a preliminary screening of a large database of compounds via similarity and pharmacophore, a secondary selection via structure-based affinity prediction and molecular docking, and a rescoring treatment for the final selection. MD simulations and MM/GBSA calculations were performed to check the agreement between different prediction methods for these inhibitors. 331 experimental ligands were collected, and 170 were used to build the structure-activity relationship. Eight representative ligand structural models were employed in similarity searching and pharmacophore screening over 14 million compounds. The process for selecting proper molecular descriptors provides a successful sample which can be used as a general strategy in QSAR modelling. The rescoring used in this work presents an alternative useful treatment for ranking and selection.

Quantitative structure activity relationship (QSAR) modeling study of some novel thiazolidine 4-one derivatives as potent anti-tubercular agents

This study aims to develop a QSAR model for Antitubercular activity. The quantitative structure-activity relationship (QSAR) approach predicted the thiazolidine-4-ones derivative's Antitubercular activity. For the QSAR study, 53 molecules with Antitubercular activity on H37Rv were collected from the literature. Compound structures were drawn by ACD/Labs ChemSketch. The energy minimization of the 2D structure was done using the MM2 force field in Chem3D pro. PaDEL Descriptor software was used to construct the molecular descriptors. QSARINS software was used in this work to develop the 2D QSAR model. A series of thiazolidine 4-one with MIC data were taken from the literature to develop the QSAR model. These compounds were split into a training set (43 compounds) and a test set (10 compounds). The PaDEL software calculated 2300 descriptors for this series of thiazolidine 4-one derivatives. The best predictive Model 4, which has R2 of 0.9092, R2adj of 0.8950 and LOF parameter of 0.0289 identify a preferred fit. The QSAR study resulted in a stable, predictive, and robust model representing the original dataset. In the QSAR equation, the molecular descriptor of MLFER_S, GATSe2, Shal, and EstateVSA 6 positively correlated with Antitubercular activity. While the SpMAD_Dzs 6 is negatively correlated with Antitubercular activity. The high polarizability, Electronegativities, Surface area contributions and number of Halogen atoms in the thiazolidine 4-one derivatives will increase the Antitubercular activity.

Quick and efficient quantitative predictions of androgen receptor binding affinity for screening Endocrine Disruptor Chemicals using 2D-QSAR and Chemical Read-Across

Endocrine Disruptor Chemicals are synthetic or natural molecules in the environment that promote adverse modifications of endogenous hormone regulation in humans and/or in animals. In the present research, we have applied two-dimensional quantitative structure-activity relationship (2D-QSAR) modeling to analyze the structural features of these chemicals responsible for binding to the androgen receptors (logRBA) in rats. We have collected the receptor binding data from the EDKB database (https://www.fda.gov/science-research/endocrine-disruptor-knowledge-base/accessing-edkb-database) and then employed the DTC-QSAR tool, available from https://dtclab.webs.com/software-tools, for dataset division, feature selection, and model development. The final partial least squares model was evaluated using various stringent validation criteria. From the model, we interpreted that hydrophobicity, steroidal nucleus, bulkiness and a hydrogen bond donor at an appropriate position contribute to the receptor binding affinity, while presence of electron rich features like aromaticity and polar groups decrease the receptor binding affinity. Additionally we have also performed chemical Read-Across predictions using Read-Across-v3.1 available from https://sites.google.com/jadavpuruniversity.in/dtc-lab-software/home, and the results for the external validation metrics were found to be better than the QSAR-derived predictions. The best quality of external predictions emerged from the q-RASAR approach which combines both read-across and QSAR. To explore the essential features responsible for the receptor binding, pharmacophore mapping, molecular docking along with molecular dynamics simulation were also performed, and the results are in accordance with the QSAR/q-RASAR findings.

Identification of sulfonamide-based butyrylcholinesterase inhibitors using machine learning

Aim: This study reports the designing of BChE inhibitors through machine learning (ML), followed by in silico and in vitro evaluations. Methodology: ML technique was used to predict the virtual hit, and its derivatives were synthesized and characterized. The compounds were evaluated by using various in vitro tests and in silico methods. Results: The gradient boosting classifier predicted N-phenyl-4-(phenylsulfonamido) benzamide as an active BChE inhibitor. The derivatives of the inhibitor, i.e., compounds 34, 37 and 54 were potent BChE inhibitors and displayed blood-brain barrier permeability with no significant AChE inhibition. Conclusion: The ML prediction was effective, and the synthesized compounds showed the BChE inhibitory activity, which was also supported by the in silico studies.

Identification of potential antivirals against 3CLpro enzyme for the treatment of SARS-CoV-2: A multi-step virtual screening study

The severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) outbreak is posing a serious public health threat worldwide in the form of COVD-19. Herein, we have performed two-dimensional quantitative structure-activity relationship (2D-QSAR) and three-dimensional pharmacophore modelling analysis employing inhibitors of 3-chymotrypsin-like protease (3CLpro), the leading protease that is crucial for the replication of SARS-CoV-2. The investigation aims to identify the important structural features responsible for the enzyme inhibition and the search for novel 3CLpro enzyme inhibitors as effective therapeutics for treating SARS-CoV-2. Furthermore, we carried out molecular docking studies using the most and least active compounds in the dataset, aiming to validate the contributions of various features as appeared in the QSAR models. Later, the stringently validated 2D-QSAR model was used to estimate the 3CLpro inhibitory activity of compounds from five chemical databases. Compounds with the significant predicted activity were then subjected to pharmacophore-based virtual screening to screen the top-rated compounds, which were then further subjected to molecular docking analysis, absorption, distribution, metabolism, excretion - toxicity (ADMET) profiling, and molecular dynamics (MD) simulation. The multi-step virtual screening analyses suggested that compounds CASAntiV-865453-58-3, CASAntiV-865453-40-3, and CASAntiV-2043031-84-9 could be used as effective therapeutic agents for the treatment of SARS-CoV-2.

Classification of beta-site amyloid precursor protein cleaving enzyme 1 inhibitors by using machine learning methods

The beta-site amyloid precursor protein cleaving enzyme 1 (BACE1) is a transmembrane aspartyl-protease, that cleaves amyloid precursor protein (APP) at the β-site. The sequential proteolytic cleavage of APP, first by β-secretase and then by γ-secretase complex, leads to the production and release of amyloid-β peptide, a pathological hallmark of Alzheimer's disease (AD). BACE1 inhibitors are reported to possess considerable potential in decreasing the level of amyloid-β in brain and preventing the progression of AD. A classification study has been conducted on 3536 diverse BACE1 inhibitors, obtained from Binding DB database, by extracting two types of descriptors, that is molecular property (Mordred) and fingerprints (Pubchem, MACCS and KRFP). Furthermore, based on the descriptors, various machine learning algorithms such as Naïve Bayesian (NB), nearest known neighbours (kNN), support vector machine (SVM), random forest (RF) and gradient-boosted algorithms (XGB) were applied to develop classification models. The performance of models was evaluated by using accuracy, precision, recall and F1 score of test set. The best NB, kNN, SVM, RF and XGB classifiers had F1 score of 0.74, 0.85, 0.86, 0.87 and 0.87, respectively. The diverse 3536 BACE1 inhibitors were clustered into 11 subsets, and the structural features of each subset were evaluated. The important fragments present in active and inactive compounds were also identified. The model developed in the study would serve as a valuable tool for the designing of BACE1 inhibitors, and also in virtual screening of molecules to identify these.

Deciphering the Interactions of Bioactive Compounds in Selected Traditional Medicinal Plants against Alzheimer's Diseases via Pharmacophore Modeling, Auto-QSAR, and Molecular Docking Approaches

Neurodegenerative diseases, for example Alzheimer’s, are perceived as driven by hereditary, cellular, and multifaceted biochemical actions. Numerous plant products, for example flavonoids, are documented in studies for having the ability to pass the blood-brain barrier and moderate the development of such illnesses. Computer-aided drug design (CADD) has achieved importance in the drug discovery world; innovative developments in the aspects of structure identification and characterization, bio-computational science, and molecular biology have added to the preparation of new medications towards these ailments. In this study we evaluated nine flavonoid compounds identified from three medicinal plants, namely T. diversifolia, B. sapida, and I. gabonensis for their inhibitory role on acetylcholinesterase (AChE), butyrylcholinesterase (BChE) and monoamine oxidase (MAO) activity, using pharmacophore modeling, auto-QSAR prediction, and molecular studies, in comparison with standard drugs. The results indicated that the pharmacophore models produced from structures of AChE, BChE and MAO could identify the active compounds, with a recuperation rate of the actives found near 100% in the complete ranked decoy database. Moreso, the robustness of the virtual screening method was accessed by well-established methods including enrichment factor (EF), receiver operating characteristic curve (ROC), Boltzmann-enhanced discrimination of receiver operating characteristic (BEDROC), and area under accumulation curve (AUAC). Most notably, the compounds’ pIC₅₀ values were predicted by a machine learning-based model generated by the AutoQSAR algorithm. The generated model was validated to affirm its predictive model. The best models achieved for AChE, BChE and MAO were models kpls_radial_17 (R² = 0.86 and Q² = 0.73), pls_38 (R² = 0.77 and Q² = 0.72), kpls_desc_44 (R² = 0.81 and Q² = 0.81) and these externally validated models were utilized to predict the bioactivities of the lead compounds. The binding affinity results of the ligands against the three selected targets revealed that luteolin displayed the highest affinity score of −9.60 kcal/mol, closely followed by apigenin and ellagic acid with docking scores of −9.60 and −9.53 kcal/mol, respectively. The least binding affinity was attained by gallic acid (−6.30 kcal/mol). The docking scores of our standards were −10.40 and −7.93 kcal/mol for donepezil and galanthamine, respectively. The toxicity prediction revealed that none of the flavonoids presented toxicity and they all had good absorption parameters for the analyzed targets. Hence, these compounds can be considered as likely leads for drug improvement against the same.

Amalgamation of in-silico, in-vitro and in-vivo approach to establish glabridin as a potential CYP2E1 inhibitor

CYP2E1 is directly or indirectly involved in the metabolism of ethanol and endogenous fatty acids but it plays a major role in the bio-activation of toxic substances that produce reactive metabolites leading to hepatotoxicity. Therefore, identification of CYP2E1 inhibitor from bioflavonoids class having useful pharmacological properties has dual benefit regarding avoidance of severe food-drug/nutraceutical-drug interaction and scope to develop a phytotherapeutics through an intended pharmacokinetic interaction.In the present study, we aimed to identify CYP2E1 inhibitor from experimental bioflavonoids which are unexplored for CYP2E1 inhibition till date using in-silico, in-vitro and in-vivo approaches.Results of in-vitro CYP2E1 inhibitory studies using CYP2E1-mediated chlorzoxazone 6-hydroxylation in human liver microsomes showed that glabridin have the highest potential than fisetin, epicatechin, nobiletin, and chrysin to inhibit CYP2E1 enzyme. Mechanistic investigations indicate that glabridin is a competitive CYP2E1 inhibitor. Molecular docking study results demonstrate that glabridin strongly interacted with the active site of human CYP2E1 enzyme. Pharmacokinetics of a CYP2E1 substrate in mice model indicates a significant alteration of chlorzoxazone and 6-hydroxychlorzoxazone plasma levels in the presence of glabridin. Further studies are needed to confirm the results at clinical level.Overall, glabridin is found to be a potential CYP2E1 inhibitor.

Design of Curcumin and Flavonoid Derivatives with Acetylcholinesterase and Beta-Secretase Inhibitory Activities Using in Silico Approaches

Acetylcholinesterase (AChE) and beta-secretase (BACE-1) are the two crucial enzymes involved in the pathology of Alzheimer's disease. The former is responsible for many defects in cholinergic signaling pathway and the latter is the primary enzyme in the biosynthesis of beta-amyloid as the main component of the amyloid plaques. These both abnormalities are found in the brains of Alzheimer's patients. In this study, in silico models were developed, including 3D-pharmacophore, 2D-QSAR (two-dimensional quantitative structure-activity relationship), and molecular docking, to screen virtually a database of compounds for AChE and BACE-1 inhibitory activities. A combinatorial library containing more than 3 million structures of curcumin and flavonoid derivatives was generated and screened for drug-likeness and enzymatic inhibitory bioactivities against AChE and BACE-1 through the validated in silico models. A total of 47 substances (two curcumins and 45 flavonoids), with remarkable predicted pIC₅₀ values against AChE and BACE-1 ranging from 4.24-5.11 (AChE) and 4.52-10.27 (BACE-1), were designed. The in vitro assays on AChE and BACE-1 were performed and confirmed the in silico results. The study indicated that, by using in silico methods, a series of curcumin and flavonoid structures were generated with promising predicted bioactivities. This would be a helpful foundation for the experimental investigations in the future. Designed compounds which were the most feasible for chemical synthesis could be potential candidates for further research and lead optimization.

Exploring QSAR models for assessment of acute fish toxicity of environmental transformation products of pesticides (ETPPs)

Environmental transformation products of pesticides (ETPPs) have a great deal of ecological impact owing to their ability to cause toxicity to the aquatic organisms, which can then be translated to the humans. The limited experimental data on biochemical and toxic effects of ETPPs, the high test costs together with regulatory limitations and the international push to reduce animal testing encourage greater dependence on predictive in silico techniques like quantitative structure-activity relationship (QSAR) models. The aim of the present work was to explore the key structural features, which regulate the toxicity towards fishes, for 85 ETPPs using a partial least squares (PLS) regression based chemometric model developed according to Organisation for Economic Co-operation and Development (OECD) guidelines. The model was extensively validated using both internal and external validation metrics, and the results so obtained justify the reliability and usefulness of the developed model (Q² = 0.648, R²_pred or Q²_F1 = 0.734 and Q²_F2 = 0.733). From the developed model, we can conclude that lipophilicity, polarity, presence of branching and the functional form of O-atom in the transformed structures of pesticides are the important features that are to be considered during ecotoxicity assessment of ETPPs. The information obtained from the descriptors of the developed model could be utilized in the future for assessing ETPPs with the benefit of providing an early warning of their potentially detrimental effect on fishes for regulatory purposes.