DockingApp RF: A State-of-the-Art Novel Scoring Function for Molecular Docking in a User-Friendly Interface to AutoDock Vina

Machine Learning for Sequence and Structure-Based Protein-Ligand Interaction Prediction

Developing new drugs is too expensive and time -consuming. Accurately predicting the interaction between drugs and targets will likely change how the drug is discovered. Machine learning-based protein-ligand interaction prediction has demonstrated significant potential. In this paper, computational methods, focusing on sequence and structure to study protein-ligand interactions, are examined. Therefore, this paper starts by presenting an overview of the data sets applied in this area, as well as the various approaches applied for representing proteins and ligands. Then, sequence-based and structure-based classification criteria are subsequently utilized to categorize and summarize both the classical machine learning models and deep learning models employed in protein-ligand interaction studies. Moreover, the evaluation methods and interpretability of these models are proposed. Furthermore, delving into the diverse applications of protein-ligand interaction models in drug research is presented. Lastly, the current challenges and future directions in this field are addressed.

Molecular Dynamics as a Tool for Virtual Ligand Screening

Rational drug design is essential for new drugs to emerge, especially when the structure of a target protein or nucleic acid is known. To that purpose, high-throughput virtual ligand screening campaigns aim at discovering computationally new binding molecules or fragments to modulate particular biomolecular interactions or biological activities, related to a disease process. The structure-based virtual ligand screening process primarily relies on docking methods which allow predicting the binding of a molecule to a biological target structure with a correct conformation and the best possible affinity. The docking method itself is not sufficient as it suffers from several and crucial limitations (lack of full protein flexibility information, no solvation and ion effects, poor scoring functions, and unreliable molecular affinity estimation).At the interface of computer techniques and drug discovery, molecular dynamics (MD) allows introducing protein flexibility before or after a docking protocol, refining the structure of protein-drug complexes in the presence of water, ions, and even in membrane-like environments, describing more precisely the temporal evolution of the biological complex and ranking these complexes with more accurate binding energy calculations. In this chapter, we describe the up-to-date MD, which plays the role of supporting tools in the virtual ligand screening (VS) process.Without a doubt, using docking in combination with MD is an attractive approach in structure-based drug discovery protocols nowadays. It has proved its efficiency through many examples in the literature and is a powerful method to significantly reduce the amount of required wet experimentations (Tarcsay et al, J Chem Inf Model 53:2990-2999, 2013; Barakat et al, PLoS One 7:e51329, 2012; De Vivo et al, J Med Chem 59:4035-4061, 2016; Durrant, McCammon, BMC Biol 9:71-79, 2011; Galeazzi, Curr Comput Aided Drug Des 5:225-240, 2009; Hospital et al, Adv Appl Bioinforma Chem 8:37-47, 2015; Jiang et al, Molecules 20:12769-12786, 2015; Kundu et al, J Mol Graph Model 61:160-174, 2015; Mirza et al, J Mol Graph Model 66:99-107, 2016; Moroy et al, Future Med Chem 7:2317-2331, 2015; Naresh et al, J Mol Graph Model 61:272-280, 2015; Nichols et al, J Chem Inf Model 51:1439-1446, 2011; Nichols et al, Methods Mol Biol 819:93-103, 2012; Okimoto et al, PLoS Comput Biol 5:e1000528, 2009; Rodriguez-Bussey et al, Biopolymers 105:35-42, 2016; Sliwoski et al, Pharmacol Rev 66:334-395, 2014).

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.

Dockey: a modern integrated tool for large-scale molecular docking and virtual screening

Molecular docking is a structure-based and computer-aided drug design approach that plays a pivotal role in drug discovery and pharmaceutical research. AutoDock is the most widely used molecular docking tool for study of protein-ligand interactions and virtual screening. Although many tools have been developed to streamline and automate the AutoDock docking pipeline, some of them still use outdated graphical user interfaces and have not been updated for a long time. Meanwhile, some of them lack cross-platform compatibility and evaluation metrics for screening lead compound candidates. To overcome these limitations, we have developed Dockey, a flexible and intuitive graphical interface tool with seamless integration of several useful tools, which implements a complete docking pipeline covering molecular sanitization, molecular preparation, paralleled docking execution, interaction detection and conformation visualization. Specifically, Dockey can detect the non-covalent interactions between small molecules and proteins and perform cross-docking between multiple receptors and ligands. It has the capacity to automatically dock thousands of ligands to multiple receptors and analyze the corresponding docking results in parallel. All the generated data will be kept in a project file that can be shared between any systems and computers with the pre-installation of Dockey. We anticipate that these unique characteristics will make it attractive for researchers to conduct large-scale molecular docking without complicated operations, particularly for beginners. Dockey is implemented in Python and freely available at https://github.com/lmdu/dockey.

Identification of novel human nicotinamide N-methyltransferase inhibitors: a structure-based pharmacophore modeling and molecular dynamics approach

Human nicotinamide N-methyltransferase (hNNMT) is a cytosolic enzyme associated in the phase-II metabolism, belonging to the S-adenosyl-L-methionine (SAM)-dependent methyltransferases family. Overexpression of hNNMT was observed in diseases such as metabolic disorders and different types of cancers, which suggest NNMT as a prospective therapeutic target. In this study we propose a structure-based pharmacophore model to understand the structural features responsible for the pharmacological activity. The generated model was validated using the ROC curve (AUC), goodness of hit score (GH), specificity, sensitivity and enrichment factor (EF). The pharmacophore was employed to retrieve active molecules from the ZINC database, followed by virtual-screening and molecular docking. Six molecules with the best pharmfit score, binding energy and ADMET properties were identified in this study. A 150 ns molecular dynamics simulation was performed on the selected molecules complexed with hNNMT protein to validate the results. The molecules ZINC35464499, ZINC13311192, ZINC31159282, ZINC14650833, ZINC14819515 and ZINC00303881 were identified, which could be act as the potential hNNMT inhibitors and can also be used as direct hits for developing novel hNNMT antagonists.Communicated by Ramaswamy H. Sarma.

Systematic Improvement of the Performance of Machine Learning Scoring Functions by Incorporating Features of Protein-Bound Water Molecules

Water molecules at the ligand-protein interfaces play crucial roles in the binding of the ligands, but the behavior of protein-bound water is largely ignored in many currently used machine learning (ML)-based scoring functions (SFs). In an attempt to improve the prediction performance of existing ML-based SFs, we estimated the water distribution with a HydraMap (HM) method and then incorporated the features extracted from protein-bound waters obtained in this way into three ML-based SFs: RF-Score, ECIF, and PLEC. It was found that a combination of HM-based features can consistently improve the performance of all three SFs, including their scoring, ranking, and docking power. HydraMap-based features show consistently good performance with both crystal structures and docked structures, demonstrating their robustness for SFs. Overall, HM-based features, which are a statistical representation of hydration sites at protein-ligand interfaces, are expected to improve the prediction performance for diverse SFs.

Antagonistic Effect and In Vitro Activity of Dauricine on Glucagon Receptor

Abnormal increases in glucagon (GCG) are the primary cause of type II diabetes mellitus. When GCG interacts with a glucagon receptor (GCGR), GCG can increase the blood glucose level. In this paper, a compound that could interfere with the binding of GCG and GCGR to inhibit the increase of blood glucose was investigated. First, molecular docking was used to conduct preliminary screening of compounds whose active components could combine with GCGR by AutoDock Vina. The binding of the receptor-ligand complex was analyzed by PyMOL. Results showed that dauricine could tightly bind to the receptor pocket. Second, the plasmid pcDNA3.1(+)-GCGR containing the target gene was transfected into HEK293 cells for expression, which was the cell model established to screen GCGR antagonist. Dauricine, the lead compound of glucagon receptor antagonist (GRA), was screened using the GRA screening model in vitro. Finally, using [Des-His1, Glu9]-Glucagon amide as the positive control, flow cytometry was used to express the antagonistic effect of the compound. Consequently, dauricine can antagonize the GCGR.

Artificial intelligence in virtual screening: models versus experiments

A typical drug discovery project involves identifying active compounds with significant binding potential for selected disease-specific targets. Experimental high-throughput screening (HTS) is a traditional approach to drug discovery, but is expensive and time-consuming when dealing with huge chemical libraries with billions of compounds. The search space can be narrowed down with the use of reliable computational screening approaches. In this review, we focus on various machine-learning (ML) and deep-learning (DL)-based scoring functions developed for solving classification and ranking problems in drug discovery. We highlight studies in which ML and DL models were successfully deployed to identify lead compounds for which the experimental validations are available from bioassay studies.

Tianma Formula Alleviates Dementia via ACER2-Mediated Sphingolipid Signaling Pathway Involving Aβ

Objective

To reveal the molecular mechanism of the antagonistic effect of traditional Chinese medicine Tianma formula (TF) on dementia including vascular dementia (VaD) and Alzheimer's disease (AD) and to provide a scientific basis for the study of traditional Chinese medicine for prevention and treatment of dementia.

Method

The TF was derived from the concerted application of traditional Chinese medicine. We detected the pharmacological effect of TF in VaD rats. The molecular mechanism of TF was examined by APP/PS1 mice in vivo, Caenorhabditis elegans (C. elegans) in vitro, ELISA, pathological assay via HE staining, and transcriptome. Based on RNA-seq analysis in VaD rats, the differentially expressed genes (DEGs) were identified and then verified by quantitative PCR (qPCR) and ELISA. The molecular mechanisms of TF on dementia were further confirmed by network pharmacology and molecular docking finally.

Results

The Morris water maze showed that TF could improve the cognitive memory function of the VaD rats. The ELISA and histological analysis suggested that TF could protect the hippocampus via reducing tau and IL-6 levels and increasing SYN expression. Meanwhile, it could protect the neurological function by alleviating Aβ deposition in APP/PS1 mice and C. elegans. In the RNA-seq analysis, 3 sphingolipid metabolism pathway-related genes, ADORA3, FCER1G, and ACER2, and another 5 nerve-related genes in 45 key DEGs were identified, so it indicated that the protection mechanism of TF was mainly associated with the sphingolipid metabolism pathway. In the qPCR assay, compared with the model group, the mRNA expressions of the 8 genes mentioned above were upregulated, and these results were consistent with RNA-seq. The protein and mRNA levels of ACER2 were also upregulated. Also, the results of network pharmacology analysis and molecular docking were consistent with those of RNA-seq analysis.

Conclusion

TF alleviates dementia by reducing Aβ deposition via the ACER2-mediated sphingolipid signaling pathway.

TSSF-hERG: A machine-learning-based hERG potassium channel-specific scoring function for chemical cardiotoxicity prediction

The human ether-à-go-go-related gene (hERG) encodes the Kv11.1 voltage-gated potassium ion (K⁺) channel that conducts the rapidly activating delayed rectifier current (I_Kr) in cardiomyocytes to regulate the repolarization process. Some drugs, as blockers of hERG potassium channels, cannot be marketed due to prolonged QT intervals, as well known as cardiotoxicity. Predetermining the binding affinity values between drugs and hERG through in silico methods can greatly reduce the time and cost required for experimental verification. In this study, we collected 9,215 compounds with AutoDock Vina's docking structures as training set, and collected compounds from four references as test sets. A series of models for predicting the binding affinities of hERG blockers were built based on five machine learning algorithms and combinations of interaction features and ligand features. The model built by support vector regression (SVR) using the combination of all features achieved the best performance on both tenfold cross-validation and external verification, which was selected and named as TSSF-hERG (target-specific scoring function for hERG). TSSF-hERG is more accurate than the classic scoring function of AutoDock Vina and the machine-learning-based generic scoring function RF-Score, with a Pearson's correlation coefficient (Rp) of 0.765, a Spearman's rank correlation coefficient (Rs) of 0.757, a root-mean-square error (RMSE) of 0.585 in a tenfold cross-validation study. All results demonstrated that TSSF-hERG would be useful for improving the power of binding affinity prediction between hERG and compounds, which can be further used for prediction or virtual screening of the hERG-related cardiotoxicity of drug candidates.