Unbinding Kinetics of Muscarinic M3 Receptor Antagonists Explained by Metadynamics Simulations

The residence time (RT), the time for which a drug remains bound to its biological target, is a critical parameter for drug design. The prediction of this key kinetic property has been proven to be challenging and computationally demanding in the framework of atomistic simulations. In the present work, we setup and applied two distinct metadynamics protocols to estimate the RTs of muscarinic M3 receptor antagonists. In the first method, derived from the conformational flooding approach, the kinetics of unbinding is retrieved from a physics-based parameter known as the acceleration factor α (i.e., the running average over time of the potential deposited in the bound state). Such an approach is expected to recover the absolute RT value for a compound of interest. In the second method, known as the t_META-D approach, a qualitative estimation of the RT is given by the time of simulation required to drive the ligand from the binding site to the solvent bulk. This approach has been developed to reproduce the change of experimental RTs for compounds targeting the same target. Our analysis shows that both computational protocols are able to rank compounds in agreement with their experimental RTs. Quantitative structure-kinetics relationship (SKR) models can be identified and employed to predict the impact of a chemical modification on the experimental RT once a calibration study has been performed.

Computationally driven discovery of SARS-CoV-2 Mpro inhibitors: from design to experimental validation

We report a fast-track computationally driven discovery of new SARS-CoV-2 main protease (M^pro) inhibitors whose potency ranges from mM for the initial non-covalent ligands to sub-μM for the final covalent compound (IC₅₀ = 830 ± 50 nM). The project extensively relied on high-resolution all-atom molecular dynamics simulations and absolute binding free energy calculations performed using the polarizable AMOEBA force field. The study is complemented by extensive adaptive sampling simulations that are used to rationalize the different ligand binding poses through the explicit reconstruction of the ligand–protein conformation space. Machine learning predictions are also performed to predict selected compound properties. While simulations extensively use high performance computing to strongly reduce the time-to-solution, they were systematically coupled to nuclear magnetic resonance experiments to drive synthesis and for in vitro characterization of compounds. Such a study highlights the power of in silico strategies that rely on structure-based approaches for drug design and allows the protein conformational multiplicity problem to be addressed. The proposed fluorinated tetrahydroquinolines open routes for further optimization of M^pro inhibitors towards low nM affinities.

The dominant binding mode of the QUB-00006-Int-07 main protease inhibitor during absolute binding free energy simulations.

SkeleDock: A Web Application for Scaffold Docking in PlayMolecule

SkeleDock is a scaffold docking algorithm which uses the structure of a protein-ligand complex as a template to model the binding mode of a chemically similar system. This algorithm was evaluated in the D3R Grand Challenge 4 pose prediction challenge, where it achieved competitive performance. Furthermore, we show that if crystallized fragments of the target ligand are available then SkeleDock can outperform rDock docking software at predicting the binding mode. This Application Note also addresses the capacity of this algorithm to model macrocycles and deal with scaffold hopping. SkeleDock can be accessed at https://playmolecule.org/SkeleDock/.

Can We Still Trust Docking Results? An Extension of the Applicability of DockBench on PDBbind Database

The number of entries in the Protein Data Bank (PDB) has doubled in the last decade, and it has increased tenfold in the last twenty years. The availability of an ever-growing number of structures is having a huge impact on the Structure-Based Drug Discovery (SBDD), allowing investigation of new targets and giving the possibility to have multiple structures of the same macromolecule in a complex with different ligands. Such a large resource often implies the choice of the most suitable complex for molecular docking calculation, and this task is complicated by the plethora of possible posing and scoring function algorithms available, which may influence the quality of the outcomes. Here, we report a large benchmark performed on the PDBbind database containing more than four thousand entries and seventeen popular docking protocols. We found that, even in protein families wherein docking protocols generally showed acceptable results, certain ligand-protein complexes are poorly reproduced in the self-docking procedure. Such a trend in certain protein families is more pronounced, and this underlines the importance in identification of a suitable protein–ligand conformation coupled to a well-performing docking protocol.

PathwayMap: Molecular Pathway Association with Self-Normalizing Neural Networks

Drug discovery suffers from high attrition because compounds initially deemed as promising can later show ineffectiveness or toxicity resulting from a poor understanding of their activity profile. In this work, we describe a deep self-normalizing neural network model for the prediction of molecular pathway association and evaluate its performance, showing an AUC ranging from 0.69 to 0.91 on a set of compounds extracted from ChEMBL and from 0.81 to 0.83 on an external data set provided by Novartis. We finally discuss the applicability of the proposed model in the domain of lead discovery. A usable application is available via PlayMolecule.org .

Synthesis and preliminary structure-activity relationship study of 2-aryl-2H-pyrazolo[4,3-c]quinolin-3-ones as potential checkpoint kinase 1 (Chk1) inhibitors

The serine-threonine checkpoint kinase 1 (Chk1) plays a critical role in the cell cycle arrest in response to DNA damage. In the last decade, Chk1 inhibitors have emerged as a novel therapeutic strategy to potentiate the anti-tumour efficacy of cytotoxic chemotherapeutic agents. In the search for new Chk1 inhibitors, a congeneric series of 2-aryl-2 H-pyrazolo[4,3-c]quinolin-3-one (PQ) was evaluated by in-vitro and in-silico approaches for the first time. A total of 30 PQ structures were synthesised in good to excellent yields using conventional or microwave heating, highlighting that 14 of them are new chemical entities. Noteworthy, in this preliminary study two compounds 4e₂ and 4h₂ have shown a modest but significant reduction in the basal activity of the Chk1 kinase. Starting from these preliminary results, we have designed the second generation of analogous in this class and further studies are in progress in our laboratories.

DockBench as docking selector tool: the lesson learned from D3R Grand Challenge 2015

Structure-based drug design (SBDD) has matured within the last two decades as a valuable tool for the optimization of low molecular weight lead compounds to highly potent drugs. The key step in SBDD requires knowledge of the three-dimensional structure of the target-ligand complex, which is usually determined by X-ray crystallography. In the absence of structural information for the complex, SBDD relies on the generation of plausible molecular docking models. However, molecular docking protocols suffer from inaccuracies in the description of the interaction energies between the ligand and the target molecule, and often fail in the prediction of the correct binding mode. In this context, the appropriate selection of the most accurate docking protocol is absolutely relevant for the final molecular docking result, even if addressing this point is absolutely not a trivial task. D3R Grand Challenge 2015 has represented a precious opportunity to test the performance of DockBench, an integrate informatics platform to automatically compare RMDS-based molecular docking performances of different docking/scoring methods. The overall performance resulted in the blind prediction are encouraging in particular for the pose prediction task, in which several complex were predicted with a sufficient accuracy for medicinal chemistry purposes.

New Trends in Inspecting GPCR-ligand Recognition Process: the Contribution of the Molecular Modeling Section (MMS) at the University of Padova

In this review, we present a survey of the recent advances carried out by our research groups in the field of ligand-GPCRs recognition process simulations recently implemented at the Molecular Modeling Section (MMS) of the University of Padova. We briefly describe a platform of tools we have tuned to aid the identification of novel GPCRs binders and the better understanding of their binding mechanisms, based on two extensively used computational techniques such as molecular docking and MD simulations. The developed methodologies encompass: (i) the selection of suitable protocols for docking studies, (ii) the exploration of the dynamical evolution of ligand-protein interaction networks, (iii) the detailed investigation of the role of water molecules upon ligand binding, and (iv) a glance at the way the ligand might go through prior reaching the binding site.

DockBench: An Integrated Informatic Platform Bridging the Gap between the Robust Validation of Docking Protocols and Virtual Screening Simulations

Virtual screening (VS) is a computational methodology that streamlines the drug discovery process by reducing costs and required resources through the in silico identification of potential drug candidates. Structure-based VS (SBVS) exploits knowledge about the three-dimensional (3D) structure of protein targets and uses the docking methodology as search engine for novel hits. The success of a SBVS campaign strongly depends upon the accuracy of the docking protocol used to select the candidates from large chemical libraries. The identification of suitable protocols is therefore a crucial step in the setup of SBVS experiments. Carrying out extensive benchmark studies, however, is usually a tangled task that requires users' proficiency in handling different file formats and philosophies at the basis of the plethora of existing software packages. We present here DockBench 1.0, a platform available free of charge that eases the pipeline by automating the entire procedure, from docking benchmark to VS setups. In its current implementation, DockBench 1.0 handles seven docking software packages and offers the possibility to test up to seventeen different protocols. The main features of our platform are presented here and the results of the benchmark study of human Checkpoint kinase 1 (hChk1) are discussed as validation test.

ALK kinase domain mutations in primary anaplastic large cell lymphoma: consequences on NPM-ALK activity and sensitivity to tyrosine kinase inhibitors

ALK inhibitor crizotinib has shown potent antitumor activity in children with refractory Anaplastic Large Cell Lymphoma (ALCL) and the opportunity to include ALK inhibitors in first-line therapies is oncoming. However, recent studies suggest that crizotinib-resistance mutations may emerge in ALCL patients. In the present study, we analyzed ALK kinase domain mutational status of 36 paediatric ALCL patients at diagnosis to identify point mutations and gene aberrations that could impact on NPM-ALK gene expression, activity and sensitivity to small-molecule inhibitors. Amplicon ultra-deep sequencing of ALK kinase domain detected 2 single point mutations, R335Q and R291Q, in 2 cases, 2 common deletions of exon 23 and 25 in all the patients, and 7 splicing-related INDELs in a variable number of them. The functional impact of missense mutations and INDELs was evaluated. Point mutations were shown to affect protein kinase activity, signalling output and drug sensitivity. INDELs, instead, generated kinase-dead variants with dominant negative effect on NPM-ALK kinase, in virtue of their capacity of forming non-functional heterocomplexes. Consistently, when co-expressed, INDELs increased crizotinib inhibitory activity on NPM-ALK signal processing, as demonstrated by the significant reduction of STAT3 phosphorylation. Functional changes in ALK kinase activity induced by both point mutations and structural rearrangements were resolved by molecular modelling and dynamic simulation analysis, providing novel insights into ALK kinase domain folding and regulation. Therefore, these data suggest that NPM-ALK pre-therapeutic mutations may be found at low frequency in ALCL patients. These mutations occur randomly within the ALK kinase domain and affect protein activity, while preserving responsiveness to crizotinib.

Exploring the recognition pathway at the human A2A adenosine receptor of the endogenous agonist adenosine using supervised molecular dynamics simulations

The putative adenosine-hA_2A AR recognition pathway is suggested by a series of Supervised Molecular Dynamics (SuMD) simulations.

Alternative quality assessment strategy to compare performances of GPCR-ligand docking protocols: the human adenosine A(2A) receptor as a case study

The progress made in the field of G protein-coupled receptors (GPCRs) structural determination has increased the adoption of docking-driven approaches for the identification or optimization of novel potent and selective ligands. In this work, we compared the performances of the 16 different docking/scoring combinations using the recently released crystal structures of the human A2A AR (hA2A AR) in complex with both agonists and antagonists. The proposed evaluation strategy encompasses the use of three complementary "quality descriptors": (a) the number of conformations generated by a docking algorithm having a RMSD value lower than the crystal structure resolution (R); (b) a novel consensus-based function defined as "protocol score"; and (c) the interaction energy maps (IEMs) analysis, based on the identification of key ligand-receptor interactions observed in the crystal structures.

Implementing the "Best Template Searching" tool into Adenosiland platform

Background

Adenosine receptors (ARs) belong to the G protein-coupled receptors (GCPRs) family. The recent release of X-ray structures of the human A_2A AR (h A_2A AR ) in complex with agonists and antagonists has increased the application of structure-based drug design approaches to this class of receptors. Among them, homology modeling represents the method of choice to gather structural information on the other receptor subtypes, namely A₁, A_2B, and A₃ ARs. With the aim of helping users in the selection of either a template to build its own models or ARs homology models publicly available on our platform, we implemented our web-resource dedicated to ARs, Adenosiland, with the “Best Template Searching” facility. This tool is freely accessible at the following web address: http://mms.dsfarm.unipd.it/Adenosiland/ligand.php.

Findings

The template suggestions and homology models provided by the “Best Template Searching” tool are guided by the similarity of a query structure (putative or known ARs ligand) with all ligands co-crystallized with hA_2A AR subtype. The tool computes several similarity indexes and sort the outcoming results according to the index selected by the user.

Conclusions

We have implemented our web-resource dedicated to ARs Adenosiland with the “Best Template Searching” facility, a tool to guide template and models selection for hARs modelling. The underlying idea of our new facility, that is the selection of a template (or models built upon a template) whose co-crystallized ligand shares the highest similarity with the query structure, can be easily extended to other GPCRs.