SAnDReS: A Computational Tool for Docking

Deep Learning in Drug Design: Protein-Ligand Binding Affinity Prediction

Computational drug design relies on the calculation of binding strength between two biological counterparts especially a chemical compound, i.e., a ligand, and a protein. Predicting the affinity of protein-ligand binding with reasonable accuracy is crucial for drug discovery, and enables the optimization of compounds to achieve better interaction with their target protein. In this paper, we propose a data-driven framework named DeepAtom to accurately predict the protein-ligand binding affinity. With 3D Convolutional Neural Network (3D-CNN) architecture, DeepAtom could automatically extract binding related atomic interaction patterns from the voxelized complex structure. Compared with the other CNN based approaches, our light-weight model design effectively improves the model representational capacity, even with the limited available training data. We carried out validation experiments on the PDBbind v.2016 benchmark and the independent Astex Diverse Set. We demonstrate that the less feature engineering dependent DeepAtom approach consistently outperforms the other baseline scoring methods. We also compile and propose a new benchmark dataset to further improve the model performances. With the new dataset as training input, DeepAtom achieves Pearson's R=0.83 and RMSE=1.23 pK units on the PDBbind v.2016 core set. The promising results demonstrate that DeepAtom models can be potentially adopted in computational drug development protocols such as molecular docking and virtual screening.

Computational Prediction of Binding Affinity for CDK2-ligand Complexes. A Protein Target for Cancer Drug Discovery

BACKGROUND

CDK2 participates in the control of eukaryotic cell-cycle progression. Due to the great interest in CDK2 for drug development and the relative easiness in crystallizing this enzyme, we have over 400 structural studies focused on this protein target. This structural data is the basis for the development of computational models to estimate CDK2-ligand binding affinity.

OBJECTIVE

This work focuses on the recent developments in the application of supervised machine learning modeling to develop scoring functions to predict the binding affinity of CDK2.

METHOD

We employed the structures available at the protein data bank and the ligand information accessed from the BindingDB, Binding MOAD, and PDBbind to evaluate the predictive performance of machine learning techniques combined with physical modeling used to calculate binding affinity. We compared this hybrid methodology with classical scoring functions available in docking programs.

RESULTS

Our comparative analysis of previously published models indicated that a model created using a combination of a mass-spring system and cross-validated Elastic Net to predict the binding affinity of CDK2-inhibitor complexes outperformed classical scoring functions available in AutoDock4 and AutoDock Vina.

CONCLUSION

All studies reviewed here suggest that targeted machine learning models are superior to classical scoring functions to calculate binding affinities. Specifically for CDK2, we see that the combination of physical modeling with supervised machine learning techniques exhibits improved predictive performance to calculate the protein-ligand binding affinity. These results find theoretical support in the application of the concept of scoring function space.

Machine Learning-Based Scoring Functions, Development and Applications with SAnDReS

BACKGROUND

Analysis of atomic coordinates of protein-ligand complexes can provide three-dimensional data to generate computational models to evaluate binding affinity and thermodynamic state functions. Application of machine learning techniques can create models to assess protein-ligand potential energy and binding affinity. These methods show superior predictive performance when compared with classical scoring functions available in docking programs.

OBJECTIVE

Our purpose here is to review the development and application of the program SAnDReS. We describe the creation of machine learning models to assess the binding affinity of protein-ligand complexes.

METHODS

SAnDReS implements machine learning methods available in the scikit-learn library. This program is available for download at https://github.com/azevedolab/sandres. SAnDReS uses crystallographic structures, binding and thermodynamic data to create targeted scoring functions.

RESULTS

Recent applications of the program SAnDReS to drug targets such as Coagulation factor Xa, cyclin-dependent kinases and HIV-1 protease were able to create targeted scoring functions to predict inhibition of these proteins. These targeted models outperform classical scoring functions.

CONCLUSION

Here, we reviewed the development of machine learning scoring functions to predict binding affinity through the application of the program SAnDReS. Our studies show the superior predictive performance of the SAnDReS-developed models when compared with classical scoring functions available in the programs such as AutoDock4, Molegro Virtual Docker and AutoDock Vina.

Docking Paradigm in Drug Design

Docking is in demand for the rational computer aided structure based drug design. A review of docking methods and programs is presented. Different types of docking programs are described. They include docking of non-covalent small ligands, protein-protein docking, supercomputer docking, quantum docking, the new generation of docking programs and the application of docking for covalent inhibitors discovery. Taking into account the threat of COVID-19, we present here a short review of docking applications to the discovery of inhibitors of SARS-CoV and SARS-CoV-2 target proteins, including our own result of the search for inhibitors of SARS-CoV-2 main protease using docking and quantum chemical post-processing. The conclusion is made that docking is extremely important in the fight against COVID-19 during the process of development of antivirus drugs having a direct action on SARS-CoV-2 target proteins.

Structural analyses and classification of novel isoniazid resistance coupled mutational landscapes in Mycobacterium tuberculosis: a combined molecular docking and MD simulation study

Drug resistance in Mycobacterium tuberculosis has become a major challenge to the current regime of treatment as well as to the containment of the disease globally. The molecular and genetic studies identified frequently occurring point mutations in the virulent protein such as KatG of M. tuberculosis resulted in the development of isoniazid tolerance in the pathogen. This study aims to analyze the structural basis of the disease mutations available in the literature as well as to predict novel alteration in the KatG which may cause similar deleterious effects. Around 15 experimentally derived mutations were included in this study and pathogenic mutational landscapes containing 60 site-specific alterations were predicted using the available in silico techniques. The effects of these mutations on the stability of the protein were studied and an exhaustive docking study was conducted for each classified perturbations, which identify the highest changes in the binding energies in p.Meth255Ile among experimental and p.Ala222Arg in computationally predicted mutations. Furthermore, the structural effects on these substitutions were analyzed using the principles of molecular dynamic simulations each for a 100 ns time scale, which validated the interaction studies. The outcome of this study may enable the identification of the novel drug resistance-associated point mutations which were not previously reported and may contribute significantly in a variety of experimental studies as well as facilitate the process of drug design and discovery.Communicated by Ramaswamy H. Sarma.

Inhibiting CDK6 Activity by Quercetin Is an Attractive Strategy for Cancer Therapy

Cyclin-dependent kinase 6 (CDK6) is a potential drug target that plays an important role in the progression of different types of cancers. We performed in silico and in vitro screening of different natural compounds and found that quercetin has a high binding affinity for the CDK6 and inhibits its activity with an IC50 = 5.89 μM. Molecular docking and a 200 ns whole atom simulation of the CDK6-quercetin complex provide insights into the binding mechanism and stability of the complex. Binding parameters ascertained by fluorescence and isothermal titration calorimetry studies revealed a binding constant in the range of 107 M-1 of quercetin to the CDK6. Thermodynamic parameters associated with the formation of the CDK6-quercetin complex suggested an electrostatic interaction-driven process. The cell-based protein expression studies in the breast (MCF-7) and lung (A549) cancer cells revealed that the treatment of quercetin decreases the expression of CDK6. Quercetin also decreases the viability and colony formation potential of selected cancer cells. Moreover, quercetin induces apoptosis, by decreasing the production of reactive oxygen species and CDK6 expression. Both in silico and in vitro studies highlight the significance of quercetin for the development of anticancer leads in terms of CDK6 inhibitors.