A combined drug discovery strategy based on machine learning and molecular docking

In silico screening and in vivo experimental validation of 15-PGDH inhibitors from traditional Chinese medicine promoting liver regeneration

The enzyme 15-hydroxyprostaglandin dehydrogenase (15-PGDH), which acts as a negative regulator of prostaglandin E2 (PGE2) levels and activity, represents a promising pharmacological target for promoting liver regeneration. In this study, we collected data on 15-PGDH homologous family proteins, their inhibitors, and traditional Chinese medicine (TCM) compounds. Leveraging machine learning and molecular docking techniques, we constructed a prediction model for virtual screening of 15-PGDH inhibitors from TCM compound library and successfully screened genistein as a potential 15-PGDH inhibitor. Through further validation, it was discovered that genistein considerably enhances liver regeneration by inhibiting 15-PGDH, resulting in a significant increase in the PGE2 level. Genistein's effectiveness suggests its potential as a novel therapeutic agent for liver diseases, highlighting this study's contribution to expanding the clinical applications of TCM.

SchNetPack Hyperparameter Optimization for a More Reliable Top Docking Scores Prediction

Options to improve the extrapolation power of the neural network designed using the SchNetPack package with respect to top docking scores prediction are presented. It is shown that hyperparameter tuning of the atomistic model representation (in the schnetpack.representation) improves the prediction of the top scoring compounds, which have characteristically a low incidence in randomized data sets for training of machine learning models. The prediction robustness is evaluated according to the mean square error (MSE) and the entropy of the average loss landscape decrease. Admittedly, the improvement of the top scoring compounds' prediction accuracy comes with the penalty of worsening the overall prediction power. It is revealed that the most impactful hyperparameter is the cutoff (5 Å is reported as the optimal choice). Other parameters (e.g., number of radial basis functions, number of interaction layers of the neural network, feature vector size or its batch size) are found to not affect the prediction robustness of the top scoring compounds in any comparable way relative to the cutoff. The MSE of the best docking score prediction (below -13 kcal/mol) improves from ca. 3.5 to 0.9 kcal/mol, while the prediction of less potent compounds (-13 to -11 kcal/mol) shows a lesser improvement, i.e., a decrease of MSE from 1.6 to 1.3 kcal/mol. Additionally, oversampling and undersampling of the training set with respect to the top scoring compounds' abundance is presented. The results indicate that the cutoff choice performs better than over- or undersampling of the training set, with undersampling performing better than oversampling.

Protein-Protein Docking Approach to GPCR Oligomerization

G-protein-coupled receptors (GPCRs), the largest family of human membrane proteins, play a crucial role in cellular control and are the target of approximately one-third of all drugs on the market. Targeting these complexes with selectivity or formulating small molecules capable of modulating receptor-receptor interactions could potentially offer novel avenues for drug discovery, fostering the development of more refined and safer pharmacotherapies. Due to the lack of experimentally derived X-ray crystallography spectra of GPCR oligomers, there is growing evidence supporting the development of new in silico approaches for predicting GPCR self-assembling structures. The significance of GPCR oligomerization, the challenges in modeling these structures, and the potential of protein-protein docking algorithms to address these challenges are discussed. The study also underscores the use of various software solutions for modeling GPCR oligomeric structures and presents practical cases where these techniques have been successfully applied.

Use of tree-based machine learning methods to screen affinitive peptides based on docking data

Screening peptides with good affinity is an important step in peptide-drug discovery. Recent advancement in computer and data science have made machine learning a useful tool in accurately affinitive-peptide screening. In current study, four different tree-based algorithms, including Classification and regression trees (CART), C5.0 decision tree (C50), Bagged CART (BAG) and Random Forest (RF), were employed to explore the relationship between experimental peptide affinities and virtual docking data, and the performance of each model was also compared in parallel. All four algorithms showed better performances on dataset pre-scaled, -centered and -PCA than other pre-processed dataset. After model re-built and hyperparameter optimization, the optimal C50 model (C50O) showed the best performances in terms of Accuracy, Kappa, Sensitivity, Specificity, F1, MCC and AUC when validated on test data and an unknown PEDV datasets evaluation (Accuracy=80.4 %). BAG and RFO (the optimal RF), as two best models during training process, did not performed as expecting during in testing and unknown dataset validations. Furthermore, the high correlation of the predictions of RFO and BAG to C50O implied the high stability and robustness of their prediction. Whereas although the good performance on unknown dataset, the poor performance in test data validation and correlation analysis indicated CARTO could not be used for future data prediction. To accurately evaluate the peptide affinity, the current study firstly gave a tree-model competition on affinitive peptide prediction by using virtual docking data, which would expand the application of machine learning algorithms in studying PepPIs and benefit the development of peptide therapeutics.

Identification of Potential p38γ Inhibitors via In Silico Screening, In Vitro Bioassay and Molecular Dynamics Simulation Studies

Protein kinase p38γ is an attractive target against cancer because it plays a pivotal role in cancer cell proliferation by phosphorylating the retinoblastoma tumour suppressor protein. Therefore, inhibition of p38γ with active small molecules represents an attractive alternative for developing anti-cancer drugs. In this work, we present a rigorous and systematic virtual screening framework to identify potential p38γ inhibitors against cancer. We combined the use of machine learning-based quantitative structure activity relationship modelling with conventional computer-aided drug discovery techniques, namely molecular docking and ligand-based methods, to identify potential p38γ inhibitors. The hit compounds were filtered using negative design techniques and then assessed for their binding stability with p38γ through molecular dynamics simulations. To this end, we identified a promising compound that inhibits p38γ activity at nanomolar concentrations and hepatocellular carcinoma cell growth in vitro in the low micromolar range. This hit compound could serve as a potential scaffold for further development of a potent p38γ inhibitor against cancer.

Computer especially AI-assisted drug virtual screening and design in traditional Chinese medicine

BACKGROUND

Traditional Chinese medicine (TCM), as a significant part of the global pharmaceutical science, the abundant molecular compounds it contains is a valuable potential source of designing and screening new drugs. However, due to the un-estimated quantity of the natural molecular compounds and diversity of the related problems drug discovery such as precise screening of molecular compounds or the evaluation of efficacy, physicochemical properties and pharmacokinetics, it is arduous for researchers to design or screen applicable compounds through old methods. With the rapid development of computer technology recently, especially artificial intelligence (AI), its innovation in the field of virtual screening contributes to an increasing efficiency and accuracy in the process of discovering new drugs.

PURPOSE

This study systematically reviewed the application of computational approaches and artificial intelligence in drug virtual filtering and devising of TCM and presented the potential perspective of computer-aided TCM development.

STUDY DESIGN

We made a systematic review following the Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) guidelines. Then screening the most typical articles for our research.

METHODS

The systematic review was performed by following the PRISMA guidelines. The databases PubMed, EMBASE, Web of Science, CNKI were used to search for publications that focused on computer-aided drug virtual screening and design in TCM.

RESULT

Totally, 42 corresponding articles were included in literature reviewing. Aforementioned studies were of great significance to the treatment and cost control of many challenging diseases such as COVID-19, diabetes, Alzheimer's Disease (AD), etc. Computational approaches and AI were widely used in virtual screening in the process of TCM advancing, which include structure-based virtual screening (SBVS) and ligand-based virtual screening (LBVS). Besides, computational technologies were also extensively applied in absorption, distribution, metabolism, excretion and toxicity (ADMET) prediction of candidate drugs and new drug design in crucial course of drug discovery.

CONCLUSIONS

The applications of computer and AI play an important role in the drug virtual screening and design in the field of TCM, with huge application prospects.

Disease-Ligand Identification Based on Flexible Neural Tree

In order to screen the disease-related compounds of a traditional Chinese medicine prescription in network pharmacology research accurately, a new virtual screening method based on flexible neural tree (FNT) model, hybrid evolutionary method and negative sample selection algorithm is proposed. A novel hybrid evolutionary algorithm based on the Grammar-guided genetic programming and salp swarm algorithm is proposed to infer the optimal FNT. According to hypertension, diabetes, and Corona Virus Disease 2019, disease-related compounds are collected from the up-to-date literatures. The unrelated compounds are chosen by negative sample selection algorithm. ECFP6, MACCS, Macrocycle, and RDKit are utilized to numerically characterize the chemical structure of each compound collected, respectively. The experiment results show that our proposed method performs better than classical classifiers [Support Vector Machine (SVM), random forest (RF), AdaBoost, decision tree (DT), Gradient Boosting Decision Tree (GBDT), KNN, logic regression (LR), and Naive Bayes (NB)], up-to-date classifier (gcForest), and deep learning method (forgeNet) in terms of AUC, ROC, TPR, FPR, Precision, Specificity, and F1. MACCS method is suitable for the maximum number of classifiers. All methods perform poorly with ECFP6 molecular descriptor.

Classification models and SAR analysis on thromboxane A2 synthase inhibitors by machine learning methods

Thromboxane A2 synthase (TXS) is a promising drug target for cardiovascular diseases and cancer. In this work, we conducted a structure-activity relationship (SAR) study on 526 TXS inhibitors for bioactivity prediction. Three types of descriptors (MACCS fingerprints, ECFP4 fingerprints, and MOE descriptors) were utilized to characterize inhibitors, 24 classification models were developed by support vector machine (SVM), random forest (RF), extreme gradient boosting (XGBoost), and deep neural networks (DNN). Then we reduced the number of fingerprints according to the contribution of descriptors to the models, and constructed 16 extra models on simplified fingerprints. In general, Model_4D built by DNN algorithm and 67 bits MACCS fingerprints performs best. The prediction accuracy of the model on the test set is 0.969, and Matthews correlation coefficient (MCC) is 0.936. The distance between compound and model (dSTD-PRO) was used to characterize the application domain of the model. In the test set of Model_4D, dSTD-PRO of 91.5% compounds is lower than the corresponding training set threshold (threshold0.90 = 0.1055), and the accuracy of these compounds is 0.983. In addition, the important descriptors were summarized and further analyzed. It showed that aromatic nitrogenous heterocyclic groups were beneficial to improve the bioactivity of TXS inhibitors.

A review of machine learning approaches for drug synergy prediction in cancer

Combinational pharmacotherapy with the synergistic/additive effect is a powerful treatment strategy for complex diseases such as malignancies. Identifying synergistic combinations with various compounds and structures requires testing a large number of compound combinations. However, in practice, examining different compounds by in vivo and in vitro approaches is costly, infeasible and challenging. In the last decades, significant success has been achieved by expanding computational methods in different pharmacological and bioinformatics domains. As promising tools, computational approaches such as machine learning algorithms (MLAs) are used for prioritizing combinational pharmacotherapies. This review aims to provide the models developed to predict synergistic drug combinations in cancer by MLAs with various information, including gene expression, protein-protein interactions, metabolite interactions, pathways and pharmaceutical information such as chemical structure, molecular descriptor and drug-target interactions.

Machine learning prediction of 3CLpro SARS-CoV-2 docking scores

Molecular docking results of two training sets containing 866 and 8,696 compounds were used to train three different machine learning (ML) approaches. Neural network approaches according to Keras and TensorFlow libraries and the gradient boosted decision trees approach of XGBoost were used with DScribe’s Smooth Overlap of Atomic Positions molecular descriptors. In addition, neural networks using the SchNetPack library and descriptors were used. The ML performance was tested on three different sets, including compounds for future organic synthesis. The final evaluation of the ML predicted docking scores was based on the ZINC in vivo set, from which 1,200 compounds were randomly selected with respect to their size. The results obtained showed a consistent ML prediction capability of docking scores, and even though compounds with more than 60 atoms were found slightly overestimated they remain valid for a subsequent evaluation of their drug repurposing suitability.

Predicting New Anti-Norovirus Inhibitor With the Help of Machine Learning Algorithms and Molecular Dynamics Simulation-Based Model

Hepatitis C virus (HCV) inhibitors are essential in the treatment of human norovirus (HuNoV). This study aimed to map out HCV NS5B RNA-dependent RNA polymerase inhibitors that could potentially be responsible for the inhibitory activity of HuNoV RdRp. It is necessary to develop robust machine learning and in silico methods to predict HuNoV RdRp compounds. In this study, Naïve Bayesian and random forest models were built to categorize norovirus RdRp inhibitors from the non-inhibitors using their molecular descriptors and PubChem fingerprints. The best model observed had accuracy, specificity, and sensitivity values of 98.40%, 97.62%, and 97.62%, respectively. Meanwhile, an external test set was used to validate model performance before applicability to the screened HCV compounds database. As a result, 775 compounds were predicted as NoV RdRp inhibitors. The pharmacokinetics calculations were used to filter out the inhibitors that lack drug-likeness properties. Molecular docking and molecular dynamics simulation investigated the inhibitors' binding modes and residues critical for the HuNoV RdRp receptor. The most active compound, CHEMBL167790, closely binds to the binding pocket of the RdRp enzyme and depicted stable binding with RMSD 0.8-3.2 Å, and the RMSF profile peak was between 1.0-4.0 Å, and the conformational fluctuations were at 450-460 residues. Moreover, the dynamic residue cross-correlation plot also showed the pairwise correlation between the binding residues 300-510 of the HuNoV RdRp receptor and CHEMBL167790. The principal component analysis depicted the enhanced movement of protein atoms. Moreover, additional residues such as Glu510 and Asn505 interacted with CHEMBL167790 via water bridge and established H-bond interactions after the simulation. http://zinc15.docking.org/substances/ZINC000013589565.

Applications of machine-learning methods for the discovery of NDM-1 inhibitors

The emergence of New Delhi metal beta-lactamase (NDM-1)-producing bacteria and their worldwide spread pose great challenges for the treatment of drug-resistant bacterial infections. These bacteria can hydrolyze most β-lactam antibacterials. Unfortunately, there are no clinically useful NDM-1 inhibitors. In the current work, we manually collected NDM-1 inhibitors reported in the past decade and established the first NDM-1 inhibitor database. Four machine-learning models were constructed using the structural and property characteristics of the collected compounds as input training set to discover potential NDM-1 inhibitors. In order to distinguish between high active inhibitors and putative positive drugs, a three-classification strategy was introduced in our study. In detail, the commonly used positive and negative divisions are converted into strongly active, weakly active, and inactive. The accuracy of the best prediction model designed based on this strategy reached 90.5%, compared with 69.14% achieved by the traditional docking-based virtual screening method. Consequently, the best model was used to virtually screen a natural product library. The safety of the selected compounds was analyzed by the ADMET prediction model based on machine learning. Seven novel NDM-1 inhibitors were identified, which will provide valuable clues for the discovery of NDM-1 inhibitors.

The Differences Between Individuals Engaging in Nonsuicidal Self-Injury and Suicide Attempt Are Complex (vs. Complicated or Simple)

BACKGROUND

Why do some people engage in nonsuicidal self-injury (NSSI) while others attempt suicide? One way to advance knowledge about this question is to shed light on the differences between people who engage in NSSI and people who attempt suicide. These groups could differ in three broad ways. First, these two groups may differ in a simple way, such that one or a small set of factors is both necessary and sufficient to accurately distinguish the two groups. Second, they might differ in a complicated way, meaning that a specific set of a large number of factors is both necessary and sufficient to accurately classify them. Third, they might differ in a complex way, with no necessary factor combinations and potentially no sufficient factor combinations. In this scenario, at the group level, complicated algorithms would either be insufficient (i.e., no complicated algorithm produces good accuracy) or unnecessary (i.e., many complicated algorithms produce good accuracy) to distinguish between groups. This study directly tested these three possibilities in a sample of people with a history of NSSI and/or suicide attempt.

METHOD

A total of 954 participants who have either engaged in NSSI and/or suicide attempt in their lifetime were recruited from online forums. Participants completed a series of measures on factors commonly associated with NSSI and suicide attempt. To test for simple differences, univariate logistic regressions were conducted. One theoretically informed multiple logistic regression model with suicidal desire, capability for suicide, and their interaction term was considered as well. To examine complicated and complex differences, multiple logistic regression and machine learning analyses were conducted.

RESULTS

No simple algorithm (i.e., single factor or small set of factors) accurately distinguished between groups. Complicated algorithms constructed with cross-validation methods produced fair accuracy; complicated algorithms constructed with bootstrap optimism methods produced good accuracy, but multiple different algorithms with this method produced similar results.

CONCLUSIONS

Findings were consistent with complex differences between people who engage in NSSI and suicide attempts. Specific complicated algorithms were either insufficient (cross-validation results) or unnecessary (bootstrap optimism results) to distinguish between these groups with high accuracy.

Interaction between piperine and genes associated with sciatica and its mechanism based on molecular docking technology and network pharmacology

Abstract

Piperine is the main active component of Piper longum L., which is also the main component of anti-sciatica Mongolian medicine Naru Sanwei pill. It has many pharmacological activities such as anti-inflammatory and immune regulation. This paper aims to preliminarily explore the potential mechanism of piperine in the treatment of sciatica through network pharmacology and molecular docking. TCMSP, ETCM database and literature mining were used to collect the active compounds of Piper longum L. Swiss TargetPrediction and SuperPred server were used to find the targets of compounds. At the same time, CTD database was used to collect the targets of sciatica. Then the above targets were compared and analyzed to select the targets of anti-sciatica in Piper longum L. The Go (gene ontology) annotation and KEGG pathway of the targets were enriched and analyzed by Metascape database platform. The molecular docking between the effective components and the targets was verified by Autodock. After that, the sciatica model of rats was established and treated with piperine. The expression level of inflammatory factors and proteins in the serum and tissues of rat sciatic nerve were detected by ELISA and Western blot. HE staining and immunohistochemistry were carried out on the sciatica tissues of rats. The results showed that Piper longum L. can regulate the development of sciatica and affect the expressions of PPARG and NF-kB1 through its active ingredient piperine, and there is endogenous interaction between PPARG and NF-kB1.

Graphic abstract

Identification of Potential Type II Diabetes in a Large-Scale Chinese Population Using a Systematic Machine Learning Framework

BACKGROUND

An estimated 425 million people globally have diabetes, accounting for 12% of the world's health expenditures, and the number continues to grow, placing a huge burden on the healthcare system, especially in those remote, underserved areas.

METHODS

A total of 584,168 adult subjects who have participated in the national physical examination were enrolled in this study. The risk factors for type II diabetes mellitus (T2DM) were identified by p values and odds ratio, using logistic regression (LR) based on variables of physical measurement and a questionnaire. Combined with the risk factors selected by LR, we used a decision tree, a random forest, AdaBoost with a decision tree (AdaBoost), and an extreme gradient boosting decision tree (XGBoost) to identify individuals with T2DM, compared the performance of the four machine learning classifiers, and used the best-performing classifier to output the degree of variables' importance scores of T2DM.

RESULTS

The results indicated that XGBoost had the best performance (accuracy = 0.906, precision = 0.910, recall = 0.902, F-1 = 0.906, and AUC = 0.968). The degree of variables' importance scores in XGBoost showed that BMI was the most significant feature, followed by age, waist circumference, systolic pressure, ethnicity, smoking amount, fatty liver, hypertension, physical activity, drinking status, dietary ratio (meat to vegetables), drink amount, smoking status, and diet habit (oil loving).

CONCLUSIONS

We proposed a classifier based on LR-XGBoost which used fourteen variables of patients which are easily obtained and noninvasive as predictor variables to identify potential incidents of T2DM. The classifier can accurately screen the risk of diabetes in the early phrase, and the degree of variables' importance scores gives a clue to prevent diabetes occurrence.

A Machine Learning-Based Prediction Platform for P-Glycoprotein Modulators and Its Validation by Molecular Docking

P-glycoprotein (P-gp) is an important determinant of multidrug resistance (MDR) because its overexpression is associated with increased efflux of various established chemotherapy drugs in many clinically resistant and refractory tumors. This leads to insufficient therapeutic targeting of tumor populations, representing a major drawback of cancer chemotherapy. Therefore, P-gp is a target for pharmacological inhibitors to overcome MDR. In the present study, we utilized machine learning strategies to establish a model for P-gp modulators to predict whether a given compound would behave as substrate or inhibitor of P-gp. Random forest feature selection algorithm-based leave-one-out random sampling was used. Testing the model with an external validation set revealed high performance scores. A P-gp modulator list of compounds from the ChEMBL database was used to test the performance, and predictions from both substrate and inhibitor classes were selected for the last step of validation with molecular docking. Predicted substrates revealed similar docking poses than that of doxorubicin, and predicted inhibitors revealed similar docking poses than that of the known P-gp inhibitor elacridar, implying the validity of the predictions. We conclude that the machine-learning approach introduced in this investigation may serve as a tool for the rapid detection of P-gp substrates and inhibitors in large chemical libraries.

Computational Health Engineering Applied to Model Infectious Diseases and Antimicrobial Resistance Spread

Infectious diseases are the primary cause of mortality worldwide. The dangers of infectious disease are compounded with antimicrobial resistance, which remains the greatest concern for human health. Although novel approaches are under investigation, the World Health Organization predicts that by 2050, septicaemia caused by antimicrobial resistant bacteria could result in 10 million deaths per year. One of the main challenges in medical microbiology is to develop novel experimental approaches, which enable a better understanding of bacterial infections and antimicrobial resistance. After the introduction of whole genome sequencing, there was a great improvement in bacterial detection and identification, which also enabled the characterization of virulence factors and antimicrobial resistance genes. Today, the use of in silico experiments jointly with computational and machine learning offer an in depth understanding of systems biology, allowing us to use this knowledge for the prevention, prediction, and control of infectious disease. Herein, the aim of this review is to discuss the latest advances in human health engineering and their applicability in the control of infectious diseases. An in-depth knowledge of host–pathogen–protein interactions, combined with a better understanding of a host’s immune response and bacterial fitness, are key determinants for halting infectious diseases and antimicrobial resistance dissemination.