Comparative binding energy COMBINE analysis for understanding the binding determinants of type II dehydroquinase inhibitors

Identification of new putative inhibitors of Mycobacterium tuberculosis 3-dehydroshikimate dehydratase from a combination of ligand- and structure-based and deep learning in silico approaches

The development of new drugs against Mycobacterium tuberculosis is an essential strategy for fighting drug resistance. Although 3-dehydroquinate dehydratase (MtDHQ) is known to be a highly relevant target for M. tuberculosis, current research shows new putative inhibitors of MtDHQ selected by a large-scale ensemble-docking strategy combining ligand- and target-based chemoinformatic methods to deep learning. Initial chemical library was reduced from 216 million to approximately 460 thousand after pharmacophore, toxicity and molecular weight filters. Final library was subjected to an ensemble-docking protocol in GOLD which selected the top 300 molecules (GHITS). GHITS displayed different structures and characteristics when compared to known inhibitors (KINH). GHITS were further screened by post-docking analysis in AMMOS2 and deep learning virtual screening in DeepPurpose. DeepPurpose predicted that a number of GHITS had comparable or better affinity for the target than KINH. The best molecule was selected by consensus ranking using GOLD, AMMOS2 and DeepPurpose scores. Molecular dynamics revealed that the top hit displayed consistent and stable binding to MtDHQ, making strong interactions with active-site loop residues. Results forward new putative inhibitors of MtDHQ and reinforce the potential application of artificial intelligence methods for drug design. This work represents the first step in the validation of these molecules as inhibitors of MtDHQ.

A Protocol to Use Comparative Binding Energy Analysis to Estimate Drug-Target Residence Time

Comparative Binding Energy (COMBINE) analysis is an approach for deriving a target-specific scoring function to compute binding free energy, drug-binding kinetics, or a related property by exploiting the information contained in the three-dimensional structures of receptor-ligand complexes. Here, we describe the process of setting up and running COMBINE analysis to derive a Quantitative Structure-Kinetics Relationship (QSKR) for the dissociation rate constants (k_off) of inhibitors of a drug target. The derived QSKR model can be used to estimate residence times (τ, τ=1/k_off) for similar inhibitors binding to the same target, and it can also help to identify key receptor-ligand interactions that distinguish inhibitors with short and long residence times. Herein, we demonstrate the protocol for the application of COMBINE analysis on a dataset of 70 inhibitors of heat shock protein 90 (HSP90) belonging to 11 different chemical classes. The procedure is generally applicable to any drug target with known structural information on its complexes with inhibitors.

Active learning strategies with COMBINE analysis: new tricks for an old dog

The COMBINE method was designed to study congeneric series of compounds including structural information of ligand-protein complexes. Although very successful, the method has not received the same level of attention than other alternatives to study Quantitative Structure Active Relationships (QSAR) mainly because lack of ways to measure the uncertainty of the predictions and the need for large datasets. Active learning, a semi-supervised learning approach that makes use of uncertainty to enhance models' performance while reducing the size of the training sets, has been used in this work to address both problems. We propose two estimators of uncertainty: the pool of regressors and the distance to the training set. The performance of the methods has been evaluated by testing the resulting active learning workflows in 3 diverse datasets: HIV-1 protease inhibitors, Taxol-derivatives and BRD4 inhibitors. The proposed strategies were successful in 80% of the cases for the taxol-derivatives and BRD4 inhibitors, while outperformed random selection in the case of the HIV-1 protease inhibitors time-split. Our results suggest that AL-COMBINE might be an effective way of producing consistently superior QSAR models with a limited number of samples.

QM/MM simulations identify the determinants of catalytic activity differences between type II dehydroquinase enzymes

Type II dehydroquinase enzymes (DHQ2), recognized targets for antibiotic drug discovery, show significantly different activities dependent on the species: DHQ2 from Mycobacterium tuberculosis (MtDHQ2) and Helicobacter pylori (HpDHQ2) show a 50-fold difference in catalytic efficiency. Revealing the determinants of this activity difference is important for our understanding of biological catalysis and further offers the potential to contribute to tailoring specificity in drug design. Molecular dynamics simulations using a quantum mechanics/molecular mechanics potential, with correlated ab initio single point corrections, identify and quantify the subtle determinants of the experimentally observed difference in efficiency. The rate-determining step involves the formation of an enolate intermediate: more efficient stabilization of the enolate and transition state of the key step in MtDHQ2, mainly by the essential residues Tyr24 and Arg19, makes it more efficient than HpDHQ2. Further, a water molecule, which is absent in MtDHQ2 but involved in generation of the catalytic Tyr22 tyrosinate in HpDHQ2, was found to destabilize both the transition state and the enolate intermediate. The quantification of the contribution of key residues and water molecules in the rate-determining step of the mechanism also leads to improved understanding of higher potencies and specificity of known inhibitors, which should aid ongoing inhibitor design.

Structural basis for the development of SARS 3CL protease inhibitors from a peptide mimic to an aza-decaline scaffold

Design of inhibitors against severe acute respiratory syndrome (SARS) chymotrypsin-like protease (3CL(pro) ) is a potentially important approach to fight against SARS. We have developed several synthetic inhibitors by structure-based drug design. In this report, we reveal two crystal structures of SARS 3CL(pro) complexed with two new inhibitors based on our previous work. These structures combined with six crystal structures complexed with a series of related ligands reported by us are collectively analyzed. To these eight complexes, the structural basis for inhibitor binding was analyzed by the COMBINE method, which is a chemometrical analysis optimized for the protein-ligand complex. The analysis revealed that the first two latent variables gave a cumulative contribution ratio of r(2)  = 0.971. Interestingly, scores using the second latent variables for each complex were strongly correlated with root mean square deviations (RMSDs) of side-chain heavy atoms of Met(49) from those of the intact crystal structure of SARS-3CL(pro) (r = 0.77) enlarging the S2 pocket. The substantial contribution of this side chain (∼10%) for the explanation of pIC50 s was dependent on stereochemistry and the chemical structure of the ligand adapted to the S2 pocket of the protease. Thus, starting from a substrate mimic inhibitor, a design for a central scaffold for a low molecular weight inhibitor was evaluated to develop a further potent inhibitor. © 2015 Wiley Periodicals, Inc. Biopolymers (Pept Sci) 106: 391-403, 2016.