DUckCov: a Dynamic Undocking-Based Virtual Screening Protocol for Covalent Binders

Modeling receptor flexibility in the structure-based design of KRASG12C inhibitors

KRAS has long been referred to as an ‘undruggable’ target due to its high affinity for its cognate ligands (GDP and GTP) and its lack of readily exploited allosteric binding pockets. Recent progress in the development of covalent inhibitors of KRAS^G12C has revealed that occupancy of an allosteric binding site located between the α3-helix and switch-II loop of KRAS^G12C—sometimes referred to as the ‘switch-II pocket’—holds great potential in the design of direct inhibitors of KRAS^G12C. In studying diverse switch-II pocket binders during the development of sotorasib (AMG 510), the first FDA-approved inhibitor of KRAS^G12C, we found the dramatic conformational flexibility of the switch-II pocket posing significant challenges toward the structure-based design of inhibitors. Here, we present our computational approaches for dealing with receptor flexibility in the prediction of ligand binding pose and binding affinity. For binding pose prediction, we modified the covalent docking program CovDock to allow for protein conformational mobility. This new docking approach, termed as FlexCovDock, improves success rates from 55 to 89% for binding pose prediction on a dataset of 10 cross-docking cases and has been prospectively validated across diverse ligand chemotypes. For binding affinity prediction, we found standard free energy perturbation (FEP) methods could not adequately handle the significant conformational change of the switch-II loop. We developed a new computational strategy to accelerate conformational transitions through the use of targeted protein mutations. Using this methodology, the mean unsigned error (MUE) of binding affinity prediction were reduced from 1.44 to 0.89 kcal/mol on a set of 14 compounds. These approaches were of significant use in facilitating the structure-based design of KRAS^G12C inhibitors and are anticipated to be of further use in the design of covalent (and noncovalent) inhibitors of other conformationally labile protein targets.

CovPDB: a high-resolution coverage of the covalent protein-ligand interactome

In recent years, the drug discovery paradigm has shifted toward compounds that covalently modify disease-associated target proteins, because they tend to possess high potency, selectivity, and duration of action. The rational design of novel targeted covalent inhibitors (TCIs) typically starts from resolved macromolecular structures of target proteins in their apo or holo forms. However, the existing TCI databases contain only a paucity of covalent protein–ligand (cP–L) complexes. Herein, we report CovPDB, the first database solely dedicated to high-resolution cocrystal structures of biologically relevant cP–L complexes, curated from the Protein Data Bank. For these curated complexes, the chemical structures and warheads of pre-reactive electrophilic ligands as well as the covalent bonding mechanisms to their target proteins were expertly manually annotated. Totally, CovPDB contains 733 proteins and 1,501 ligands, relating to 2,294 cP–L complexes, 93 reactive warheads, 14 targetable residues, and 21 covalent mechanisms. Users are provided with an intuitive and interactive web interface that allows multiple search and browsing options to explore the covalent interactome at a molecular level in order to develop novel TCIs. CovPDB is freely accessible at http://www.pharmbioinf.uni-freiburg.de/covpdb/ and its contents are available for download as flat files of various formats.

Structure-based inhibitor design of mutant RAS proteins-a paradigm shift

As a member of small GTPase family, KRAS protein is a key physiological modulator of various cellular activities including proliferation. However, mutations of KRAS present in numerous cancer types, most frequently in pancreatic (> 60%), colorectal (> 40%), and lung cancers, drive oncogenic processes through overactivation of proliferation. The G12C mutation of KRAS protein is especially abundant in the case of these types of malignancies. Despite its key importance in human disease, KRAS was assumed to be non-druggable for a long time since the protein seemingly lacks potential drug-binding pockets except the nucleotide-binding site, which is difficult to be targeted due to the high affinity of KRAS for both GDP and GTP. Recently, a new approach broke the ice and provided evidence that upon covalent targeting of the G12C mutant KRAS, a highly dynamic pocket was revealed. This novel targeting is especially important since it serves with an inherent solution for drug selectivity. Based on these results, various structure-based drug design projects have been launched to develop selective KRAS mutant inhibitors. In addition to the covalent modification strategy mostly applicable for G12C mutation, different innovative solutions have been suggested for the other frequently occurring oncogenic G12 mutants. Here we summarize the latest advances of this field, provide perspectives for novel approaches, and highlight the special properties of KRAS, which might issue some new challenges.

An automatic pipeline for the design of irreversible derivatives identifies a potent SARS-CoV-2 Mpro inhibitor

Designing covalent inhibitors is increasingly important, although it remains challenging. Here, we present covalentizer, a computational pipeline for identifying irreversible inhibitors based on structures of targets with non-covalent binders. Through covalent docking of tailored focused libraries, we identify candidates that can bind covalently to a nearby cysteine while preserving the interactions of the original molecule. We found ∼11,000 cysteines proximal to a ligand across 8,386 complexes in the PDB. Of these, the protocol identified 1,553 structures with covalent predictions. In a prospective evaluation, five out of nine predicted covalent kinase inhibitors showed half-maximal inhibitory concentration (IC₅₀) values between 155 nM and 4.5 μM. Application against an existing SARS-CoV M^pro reversible inhibitor led to an acrylamide inhibitor series with low micromolar IC₅₀ values against SARS-CoV-2 M^pro. The docking was validated by 12 co-crystal structures. Together these examples hint at the vast number of covalent inhibitors accessible through our protocol.

Discovery of fragments inducing conformational effects in dynamic proteins using a second-harmonic generation biosensor

Biophysical screening of compound libraries for the identification of ligands that interact with a protein is efficient, but does typically not reveal if (or how) ligands may interfere with its functional properties. For this a biochemical/functional assay is required. But for proteins whose function is dependent on a conformational change, such assays are typically complex or have low throughput. Here we have explored a high-throughput second-harmonic generation (SHG) biosensor to detect fragments that induce conformational changes upon binding to a protein in real time and identify dynamic regions. Multiwell plate format SHG assays were developed for wild-type and six engineered single-cysteine mutants of acetyl choline binding protein (AChBP), a homologue to ligand gated ion channels (LGICs). They were conjugated with second harmonic-active labels via amine or maleimide coupling. To validate the assay, it was confirmed that the conformational changes induced in AChBP by nicotinic acetyl choline receptor (nAChR) agonists and antagonists were qualitatively different. A 1056 fragment library was subsequently screened against all variants and conformational modulators of AChBP were successfully identified, with hit rates from 9-22%, depending on the AChBP variant. A subset of four hits was selected for orthogonal validation and structural analysis. A time-resolved grating-coupled interferometry-based biosensor assay confirmed the interaction to be a reversible 1-step 1 : 1 interaction, and provided estimates of affinities and interaction kinetic rate constants (K D = 0.28-63 μM, k a = 0.1-6 μM-1 s-1, k d = 1 s-1). X-ray crystallography of two of the fragments confirmed their binding at a previously described conformationally dynamic site, corresponding to the regulatory site of LGICs. These results reveal that SHG has the sensitivity to identify fragments that induce conformational changes in a protein. A selection of fragment hits with a response profile different to known LGIC regulators was characterized and confirmed to bind to dynamic regions of the protein.

Fragment-based covalent ligand discovery

Targeted covalent inhibitors have regained widespread attention in drug discovery and have emerged as powerful tools for basic biomedical research. Fueled by considerable improvements in mass spectrometry sensitivity and sample processing, chemoproteomic strategies have revealed thousands of proteins that can be covalently modified by reactive small molecules. Fragment-based drug discovery, which has traditionally been used in a target-centric fashion, is now being deployed on a proteome-wide scale thereby expanding its utility to both the discovery of novel covalent ligands and their cognate protein targets. This powerful approach is allowing ‘high-throughput’ serendipitous discovery of cryptic pockets leading to the identification of pharmacological modulators of proteins previously viewed as “undruggable”. The reactive fragment toolkit has been enabled by recent advances in the development of new chemistries that target residues other than cysteine including lysine and tyrosine. Here, we review the emerging area of covalent fragment-based ligand discovery, which integrates the benefits of covalent targeting and fragment-based medicinal chemistry. We discuss how the two strategies synergize to facilitate the efficient discovery of new pharmacological modulators of established and new therapeutic target proteins.

Covalent fragment-based ligand discovery greatly facilitates the discovery of useful fragments for drug discovery and helps unveil chemical-tractable biological targets in native biological systems.

The latest automated docking technologies for novel drug discovery

INTRODUCTION

Molecular docking has been consolidated as one of the most important methods in the molecular modeling field. It has been recognized as a prominent tool in the study of protein-ligand complexes, to describe intermolecular interactions, to accurately predict poses of multiple ligands, to discover novel promising bioactive compounds. Molecular docking methods have evolved in terms of their accuracy and reliability; but there are pending issues to solve for improving the connection between the docking results and the experimental evidence.

AREAS COVERED

In this article, the author reviews very recent innovative molecular docking applications with special emphasis on reverse docking, treatment of protein flexibility, the use of experimental data to guide the selection of docking poses, the application of Quantum mechanics(QM) in docking, and covalent docking.

EXPERT OPINION

There are several issues being worked on in recent years that will lead to important breakthroughs in molecular docking methods in the near future These developments are related to more efficient exploration of large datasets and receptor conformations, advances in electronic description, and the use of structural information for guiding the selection of results.

Selective and Effective: Current Progress in Computational Structure-Based Drug Discovery of Targeted Covalent Inhibitors

Targeted covalent inhibitors are currently showing great promise for systems that are normally difficult to target with small molecule therapies. This renewed interest has spurred the refinement of existing computational methods as well as the designof new ones, expanding the toolbox for discovery and optimization of selectiveand effective covalent inhibitors. Commonly applied approaches are covalentdocking methods that predict the conformation of the covalent complex with known residues. More recently, a new predictive method, reactive docking, was developed, building on the growing corpus of data generated by large proteomics experiments. This method was successfully used in several 'inverse drug discovery' programs that use high-throughput techniques to isolate effective compounds based on screening of entire compound libraries based on desired phenotypes.

Combined strategies in structure-based virtual screening

The identification and optimization of lead compounds are inalienable components in drug design and discovery pipelines. As a powerful computational approach for the identification of hits with novel structural scaffolds, structure-based virtual screening (SBVS) has exhibited a remarkably increasing influence in the early stages of drug discovery. During the past decade, a variety of techniques and algorithms have been proposed and tested with different purposes in the scope of SBVS. Although SBVS has been a common and proven technology, it still shows some challenges and problems that are needed to be addressed, where the negative influence regardless of protein flexibility and the inaccurate prediction of binding affinity are the two major challenges. Here, focusing on these difficulties, we summarize a series of combined strategies or workflows developed by our group and others. Furthermore, several representative successful applications from recent publications are also discussed to demonstrate the effectiveness of the combined SBVS strategies in drug discovery campaigns.

Accelerating Rare Dissociative Processes in Biomolecules Using Selectively Scaled MD Simulations

Molecular dynamics (MD) simulations can be a powerful tool for modeling complex dissociative processes such as ligand unbinding. However, many biologically relevant dissociative processes occur on timescales that far exceed the timescales of typical MD simulations. Here, we implement and apply an enhanced sampling method in which specific energy terms in the potential energy function are selectively "scaled" to accelerate dissociative events during simulations. Using ligand unbinding as an example of a complex dissociative process, we selectively scaled up ligand-water interactions in an attempt to increase the rate of ligand unbinding. After applying our selectively scaled MD (ssMD) approach to several cyclin-dependent kinase-inhibitor complexes, we discovered that we could accelerate ligand unbinding, thereby allowing, in some cases, unbinding events to occur within as little as 2 ns. Moreover, we found that we could make realistic estimates of the initial unbinding times (τ_unbind^sim) as well as the accompanying change in free energy (ΔG^sim) of the inhibitors from our ssMD simulation data. To accomplish this, we employed a previously described Kramers'-based rate extrapolation method and a newly described free energy extrapolation method. Because our ssMD approach is general, it should find utility as an easy-to-deploy, enhanced sampling method for modeling other dissociative processes.

DUckCov: a Dynamic Undocking-Based Virtual Screening Protocol for Covalent Binders

Thanks to recent guidelines, the design of safe and effective covalent drugs has gained significant interest. Other than targeting non‐conserved nucleophilic residues, optimizing the noncovalent binding framework is important to improve potency and selectivity of covalent binders toward the desired target. Significant efforts have been made in extending the computational toolkits to include a covalent mechanism of protein targeting, like in the development of covalent docking methods for binding mode prediction. To highlight the value of the noncovalent complex in the covalent binding process, here we describe a new protocol using tethered and constrained docking in combination with Dynamic Undocking (DUck) as a tool to privilege strong protein binders for the identification of novel covalent inhibitors. At the end of the protocol, dedicated covalent docking methods were used to rank and select the virtual hits based on the predicted binding mode. By validating the method on JAK3 and KRas, we demonstrate how this fast iterative protocol can be applied to explore a wide chemical space and identify potent targeted covalent inhibitors.