The need to solvate and encapsulate hydro-sensitive molecules drives noticeable trends in the applications of cyclodextrins in the pharmaceutical industry, in foods, polymers, materials, and in agricultural science. Among them, β-cyclodextrin is one of the most used for the entrapment of phenolic acid compounds to mask the bitterness of wheat bran. In this regard, there is still a need for good data and especially for a robust predictive model that assesses the bitterness masking capabilities of β-cyclodextrin for various phenolic compounds. This study uses a dataset of 20 phenolic acids docked into the β-cyclodextrin cavity to generate three different binding constants. The data from the docking study were combined with topological, topographical, and quantum-chemical features from the ligands in a machine learning-based structure-activity relationship study. Three different models for each binding constant were computed using a combination of the genetic algorithm (GA) and multiple linear regression (MLR) approaches. The developed ML/QSAR models showed a very good performance, with high predictive ability and correlation coefficients of 0.969 and 0.984 for the training and test sets, respectively. The models revealed several factors responsible for binding with cyclodextrin, showing positive contributions toward the binding affinity values, including such features as the presence of six-membered rings in the molecule, branching, electronegativity values, and polar surface area.

It is well-known that the hepatotoxicity of drugs can significantly influence their clinical use. Despite their effective therapeutic efficacy, many drugs are severely limited in clinical applications due to significant hepatotoxicity. In response, researchers have created several machine learning-based hepatotoxicity prediction models for use in drug discovery and development. Researchers aim to predict the potential hepatotoxicity of drugs to enhance their utility. However, current hepatotoxicity prediction models often suffer from being unverified, and they fail to capture the detailed toxicological structures of predicted hepatotoxic compounds. Using the 56 chemical constituents of Gardenia jasminoides as examples, we validated the trained hepatotoxicity prediction model through literature reviews, principal component analysis (PCA), and structural comparison methods. Ultimately, we successfully developed a model with strong predictive performance and conducted visual validation. Interestingly, we discovered that the predicted hepatotoxic chemical constituents of Gardenia possess both toxic and therapeutic effects, which are likely dose-dependent. This discovery greatly contributes to our understanding of the dual nature of drug-induced hepatotoxicity.

Dipeptidyl peptidase-4 (DPP-4) inhibitors belong to a prominent group of pharmaceutical agents that are used in the governance of type 2 diabetes mellitus (T2DM). They exert their antidiabetic effects by inhibiting the incretin hormones like glucagon-like peptide-1 and glucose-dependent insulinotropic polypeptide which, play a pivotal role in the regulation of blood glucose homoeostasis in our body. DPP-4 inhibitors have emerged as an important class of oral antidiabetic drugs for the treatment of T2DM. Surprisingly, only a few 2D-QSAR studies have been reported on DPP-4 inhibitors. Here, fragment-based QSAR (Laplacian-modified Bayesian modelling and Recursive partitioning (RP) approaches have been utilized on a dataset of 108 DPP-4 inhibitors to achieve a deeper understanding of the association among their molecular structures. The Bayesian analysis demonstrated satisfactory ROC values for the training as well as the test sets. Meanwhile, the RP analysis resulted in decision tree 3 with 2 leaves (Tree 3: 2 leaves). This present study is an effort to get an insight into the pivotal fragments modulating DPP-4 inhibition.

This study aimed to use a computational approach that combined the classification-based QSAR model, molecular docking, ADME studies, and molecular dynamics (MD) to identify potential inhibitors of Fyn kinase. First, a robust classification model was developed from a dataset of 1,078 compounds with known Fyn kinase inhibitory activity, using the XGBoost algorithm. After that, molecular docking was performed between potential compounds identified from the QSAR model and Fyn kinase to assess their binding strengths and key interactions, followed by MD simulations. ADME studies were additionally conducted to preliminarily evaluate the pharmacokinetics and drug-like characteristics of these compounds. The results showed that our obtained model exhibited good predictive performance with an accuracy of 0.95 on the test set, affirming its reliability in identifying potent Fyn kinase inhibitors. Through the application of this model in conjunction with molecular docking and ADME studies, nine compounds were identified as potential Fyn kinase inhibitors, including 208 (ZINC70708110), 728 (ZINC8792432), 734 (ZINC8792187), 736 (ZINC8792350), 738 (ZINC8792286), 739 (ZINC8792309), 817 (ZINC33901069), 852 (ZINC20759145), and 1227 (ZINC100006936). MD simulations further demonstrated that the four most promising compounds, 728, 734, 736, and 852 exhibited stable binding with Fyn kinase during the simulation process. Additionally, a web-based platform ( https://fynkinase.streamlit.app/ ) has been developed to streamline the screening process. This platform enables users to predict the activity of their substances of interest on Fyn kinase from their SMILES, using our classification-based QSAR model and molecular docking.

This study aims to assess the potential bioactivity of newly designed benzodiazepine-1,2,3-triazole derivatives using in-silico methodologies, with a primary focus on elucidating their inhibitory interactions with the butyrylcholinesterase (BuChE) enzyme, which is implicated in Alzheimer\'s disease. We employed multiple linear regression (MLR) methods to conduct a quantitative structure-activity relationship (QSAR) analysis on a collection of 31 benzodiazepine-1,2,3-triazole derivatives, with the goal of investigating, assessing, and predicting their activities, as well as designing novel compounds. This approach yielded highly accurate results, with coefficients of determination (R²) of 0.77 and 0.81 for the training and test datasets, respectively. Additionally, the optimized compounds were subjected to an Absorption, Distribution, Metabolism, Excretion, and Toxicity (ADMET) analysis, demonstrating their potential as non-hepatotoxic agents with enhanced absorption and blood-brain barrier permeability. To further validate these findings, the most favorable docking conformations were analyzed using molecular dynamics (MD) simulations with GROMACS software, predicting the stability of the formed complexes. These simulations underscored the critical role of hydrogen bonds in stabilizing the compounds at the BuChE receptor binding site. The results hold great promise for the development of innovative benzodiazepine-1,2,3-triazole derivatives as effective BuChE inhibitors, potentially leading to therapeutic interventions for Alzheimer\'s disease.

Histone deacetylases constitute a group of enzymes that participate in several biological processes. Notably, inhibiting HDAC8 has become a therapeutic strategy for various diseases. The current inhibitors for HDAC8 lack selectivity and target multiple HDACs. Consequently, there is a growing recognition of the need for selective HDAC8 inhibitors to enhance the effectiveness of therapeutic interventions. In our current study, we have utilized a multi-faceted approach, including Quantitative Structure-Activity Relationship (QSAR) combined with Quantitative Read-Across Structure-Activity Relationship (q-RASAR) modeling, pharmacophore mapping, molecular docking, and molecular dynamics (MD) simulations. The developed q-RASAR model has a high statistical significance and predictive ability (Q2F1:0.778, Q2F2:0.775). The contributions of important descriptors are discussed in detail to gain insight into the crucial structural features in HDAC8 inhibition. The best pharmacophore hypothesis exhibits a high regression coefficient (0.969) and a low root mean square deviation (0.944), highlighting the importance of correctly orienting hydrogen bond acceptor (HBA), ring aromatic (RA), and zinc-binding group (ZBG) features in designing potent HDAC8 inhibitors. To confirm the results of q-RASAR and pharmacophore mapping, molecular docking analysis of the five potent compounds (44, 54, 82, 102, and 118) was performed to gain further insights into these structural features crucial for interaction with the HDAC8 enzyme. Lastly, MD simulation studies of the most active compound (54, mapped correctly with the pharmacophore hypothesis) and the least active compound (34, mapped poorly with the pharmacophore hypothesis) were carried out to validate the observations of the studies above. This study not only refines our understanding of essential structural features for HDAC8 inhibition but also provides a robust framework for the rational design of novel selective HDAC8 inhibitors which may offer insights to medicinal chemists and researchers engaged in the development of HDAC8-targeted therapeutics.

Cancer is a serious threat to human life and social development and the use of scientific methods for cancer prevention and control is necessary. In this study, HQSAR, CoMFA, CoMSIA and TopomerCoMFA methods are used to establish models of 65 imidazo[4,5-b]pyridine derivatives to explore the quantitative structure-activity relationship between their anticancer activities and molecular conformations. The results show that the cross-validation coefficients q2 of HQSAR, CoMFA, CoMSIA and TopomerCoMFA are 0.892, 0.866, 0.877 and 0.905, respectively. The non-cross-validation coefficients r2 are 0.948, 0.983, 0.995 and 0.971, respectively. The externally validated complex correlation coefficients r2pred of external validation are 0.814, 0.829, 0.758 and 0.855, respectively. The PLS analysis verifies that the QSAR models have the highest prediction ability and stability. Based on these statistics, virtual screening based on R group is performed using the ZINC database by the Topomer search technology. Finally, 10 new compounds with higher activity are designed with the screened new fragments. In order to explore the binding modes and targets between ligands and protein receptors, these newly designed compounds are conjugated with macromolecular protein (PDB ID: 1MQ4) by molecular docking technology. Furthermore, to study the nature of the newly designed compound in dynamic states and the stability of the protein-ligand complex, molecular dynamics simulation is carried out for N3, N4, N5 and N7 docked with 1MQ4 protease structure for 50 ns. A free energy landscape is computed to search for the most stable conformation. These results prove the efficient and stability of the newly designed compounds. Finally, ADMET is used to predict the pharmacology and toxicity of the 10 designed drug molecules.

This study aims to identify potential focal adhesion kinase (FAK) inhibitors through an integrated computational approach, combining mol2vec descriptor-based QSAR, molecular docking, ADMET study, and molecular dynamics simulation. A dataset of 437 compounds with known FAK inhibitory activities was used to develop QSAR models using machine learning algorithms combined with mol2vec descriptors. Subsequently, the most promising compounds were subjected to molecular docking against FAK to evaluate their binding affinities and key interactions. ADMET study and molecular dynamics simulation were also employed to investigate the pharmacokinetic, drug-like properties, and the stability of the protein-ligand complexes. The results showed that the mol2vec descriptor-based QSAR model established by support vector regression demonstrated good predictive performance (R2 = 0.813, RMSE = 0.453, MAE = 0.263 in case of training set, and R2 = 0.729, RMSE = 0.635, MAE = 0.477 in case of test set), indicating their reliability in identifying potent FAK inhibitors. Using this QSAR model and molecular docking, compound 21 (ZINC000004523722) was identified as the most potential compound, with predicted logIC50 value and binding energy of 2.59 and - 9.3 kcal/mol, respectively. The results of molecular dynamics simulation and ADMET study also further suggested its potential as a promising drug candidate. However, because our research was merely theoretical, additional in vitro and in vivo studies are required for the verification of these results.

The fluoroquinolones ciprofloxacin, danofloxacin, enoxacin, levofloxacin and lomefloxacin, occur in water bodies worldwide and therefore pose a threat to the aquatic environment. Advanced purification procedures, such as electrochemical oxidation, may act as a remedy since they contribute to eliminating contaminants and prevent micropollutants from entering open water bodies. By electrochemical treatment in a micro-flow reactor equipped with a boron-doped diamond (BDD) electrode, the fluoroquinolones were efficiently degraded. A total of 15 new products were identified using high-performance high-resolution chromatography coupled with high-resolution multifragmentation mass spectrometry. The ecotoxicity of the emerging transformation products was estimated through in silico quantitative structure activity relationship analysis. Almost all transformation products were predicted less ecotoxic than the initial compounds. The fluoroquinolone degradation followed three major mechanisms depending on the voltage during the electrochemical oxidation. At approximately 1 V, the reactions started with the elimination of molecular hydrogen from the piperazine moiety. At approx. 1.25 V, methyl and methylene groups were eliminated. At 1.5 V, hydroxyl radicals, generated at the BDD electrode, led to substitution at the piperazine ring. This novel finding of the three reactions depending on voltage contributes to the mechanistic understanding of electrochemical oxidation as potential remedy against fluoroquinolones in the aquatic environment.

To explore novel terpenoid repellents, 22 candidate terpenoid derivatives were synthesized and tested for their electroantennogram (EAG) responses and repellent activities against Aedes albopictus. The results from the EAG experiments revealed that 5-(2-hydroxypropan-2-yl)-2-methylcyclohex-2-en-1-yl formate (compound 1) induced distinct EAG responses in female Aedes albopictus. At concentrations of 0.1, 1, 10, 100, and 1000 mg/L, the EAG response values for compound 1 were 179.59, 183.99, 190.38, 193.80, and 196.66 mV, demonstrating comparable or superior effectiveness to DEET. Repellent activity analysis indicated significant repellent activity for compound 1, closest to the positive control DEET. The in silico assessment of the ADMET profile of compound 1 indicates that it successfully passed the ADMET evaluation. Molecular docking studies exhibited favourable binding of compound 1 to the active site of the odorant binding protein (OBP) of Aedes albopictus, involving hydrophobic forces and hydrogen bond interactions with residues in the OBP pocket. The QSAR model highlighted the influential role of hydrogen-bonding receptors, positively charged surface area of weighted atoms, polarity parameters of molecules, and maximum nuclear-nuclear repulsion force of carbon-carbon bonds on the relative EAG response values of the tested compounds. This study holds substantial significance for the advancement of new terpenoid repellents.