Exploiting high-energy hydration sites for the discovery of potent peptide aldehyde inhibitors of the SARS-CoV-2 main protease with cellular antiviral activity.

Small-molecule antivirals that prevent the replication of the SARS-CoV-2 virus by blocking the enzymatic activity of its main protease (Mpro) are and will be a tenet of pandemic preparedness. However, the peptidic nature of such compounds often precludes the design of compounds within favorable physical property ranges, limiting cellular activity. Here we describe the discovery of peptide aldehyde Mpro inhibitors with potent enzymatic and cellular antiviral activity. This structure-activity relationship (SAR) exploration was guided by the use of calculated hydration site thermodynamic maps (WaterMap) to drive potency via displacement of waters from high-energy sites. Thousands of diverse compounds were designed to target these high-energy hydration sites and then prioritized for synthesis by physics- and structure-based Free-Energy Perturbation (FEP+) simulations, which accurately predicted biochemical potencies. This approach ultimately led to the rapid discovery of lead compounds with unique SAR that exhibited potent enzymatic and cellular activity with excellent pan-coronavirus coverage.

Discovery of Cycloalkyl[c]thiophenes as Novel Scaffolds for Hypoxia-Inducible Factor-2α Inhibitors.

Hypoxia-inducible factors (HIFs) are heterodimeric transcription factors induced in diverse pathophysiological settings. Inhibition of HIF-2α has become a strategy for cancer treatment since the discovery that small molecules, upon binding into a small cavity of the HIF-2α PAS B domain, can alter its conformation and disturb the activity of the HIF dimer complex. Herein, the design, synthesis, and systematic SAR exploration of cycloalkyl[c]thiophenes as novel HIF-2α inhibitors are described, providing the first chemotype featuring an alkoxy-aryl scaffold. X-ray data confirmed the ability of these inhibitors to induce perturbation of key amino acids by appropriately presenting key pharmacophoric elements in the hydrophobic cavity. Selected compounds showed inhibition of VEGF-A secretion in cancer cells and prevention of Arg1 expression and activity in IL4-stimulated macrophages. Moreover, in vivo target gene modulation was demonstrated with compound 35r. Thus, the disclosed HIF-2α inhibitors represent valuable tools for investigating selective HIF-2α inhibition and its effect on tumor biology.

An ABSINTH-Based Protocol for Predicting Binding Affinities between Proteins and Small Molecules.

The core task in computational drug discovery is to accurately predict binding free energies in receptor-ligand systems for large libraries of putative binders. Here, the ABSINTH implicit solvent model and force field are extended to describe small, organic molecules and their interactions with proteins. We show that an automatic pipeline based on partitioning arbitrary molecules into substructures corresponding to model compounds with known free energies of solvation can be combined with the CHARMM general force field into a method that is successful at the two important challenges a scoring function faces in virtual screening work flows: it ranks known binders with correlation values rivaling that of comparable state-of-the-art methods and it enriches true binders in a set of decoys. Our protocol introduces innovative modifications to common virtual screening workflows, notably the use of explicit ions as competitors and the integration over multiple protein and ligand species differing in their protonation states. We demonstrate the value of modifications to both the protocol and ABSINTH itself. We conclude by discussing the limitations of high-throughput implicit methods such as the one proposed here.

Large-Scale Assessment of Binding Free Energy Calculations in Active Drug Discovery Projects.

Accurate ranking of compounds with regards to their binding affinity to a protein using computational methods is of great interest to pharmaceutical research. Physics-based free energy calculations are regarded as the most rigorous way to estimate binding affinity. In recent years, many retrospective studies carried out both in academia and industry have demonstrated its potential. Here, we present the results of large-scale prospective application of the FEP+ method in active drug discovery projects in an industry setting at Merck KGaA, Darmstadt, Germany. We compare these prospective data to results obtained on a new diverse, public benchmark of eight pharmaceutically relevant targets. Our results offer insights into the challenges faced when using free energy calculations in real-life drug discovery projects and identify limitations that could be tackled by future method development. The new public data set we provide to the community can support further method development and comparative benchmarking of free energy calculations.

Site-Specific Conjugation of Native Antibodies Using Engineered Microbial Transglutaminases.

Site-specific bioconjugation technologies are frequently employed to generate homogeneous antibody-drug conjugates (ADCs) and are generally considered superior to stochastic approaches like lysine coupling. However, most of the technologies developed so far require undesired manipulation of the antibody sequence or its glycan structures. Herein, we report the successful engineering of microbial transglutaminase enabling efficient, site-specific conjugation of drug-linker constructs to position HC-Q295 of native, fully glycosylated IgG-type antibodies. ADCs generated via this approach demonstrate excellent stability in vitro as well as strong efficacy in vitro and in vivo. As it employs different drug-linker structures and several native antibodies, our study additionally proves the broad applicability of this approach.

Discovery of a Locally and Orally Active CXCL12 Neutraligand (LIT-927) with Anti-inflammatory Effect in a Murine Model of Allergic Airway Hypereosinophilia.

We previously reported Chalcone-4 (1) that binds the chemokine CXCL12, not its cognate receptors CXCR4 or CXCR7, and neutralizes its biological activity. However, this neutraligand suffers from limitations such as poor chemical stability, solubility, and oral activity. Herein, we report on the discovery of pyrimidinone 57 (LIT-927), a novel neutraligand of CXCL12 which displays a higher solubility than 1 and is no longer a Michael acceptor. While both 1 and 57 reduce eosinophil recruitment in a murine model of allergic airway hypereosinophilia, 57 is the only one to display inhibitory activity following oral administration. Thereby, we here describe 57 as the first orally active CXCL12 neutraligand with anti-inflammatory properties. Combined with a high binding selectivity for CXCL12 over other chemokines, 57 represents a powerful pharmacological tool to investigate CXCL12 physiology in vivo and to explore the activity of chemokine neutralization in inflammatory and related diseases.

Concurrent cooperativity and substrate inhibition in the epoxidation of carbamazepine by cytochrome P450 3A4 active site mutants inspired by molecular dynamics simulations.

Cytochrome P450 3A4 (CYP3A4) is the major human P450 responsible for the metabolism of carbamazepine (CBZ). To explore the mechanisms of interactions of CYP3A4 with this anticonvulsive drug, we carried out multiple molecular dynamics (MD) simulations, starting with the complex of CYP3A4 manually docked with CBZ. On the basis of these simulations, we engineered CYP3A4 mutants I369F, I369L, A370V, and A370L, in which the productive binding orientation was expected to be stabilized, thus leading to increased turnover of CBZ to the 10,11-epoxide product. In addition, we generated CYP3A4 mutant S119A as a control construct with putative destabilization of the productive binding pose. Evaluation of the kinetics profiles of CBZ epoxidation demonstrate that CYP3A4-containing bacterial membranes (bactosomes) as well as purified CYP3A4 (wild-type and mutants I369L/F) exhibit substrate inhibition in reconstituted systems. In contrast, mutants S119A and A370V/L exhibit S-shaped profiles that are indicative of homotropic cooperativity. MD simulations with two to four CBZ molecules provide evidence that the substrate-binding pocket of CYP3A4 can accommodate more than one molecule of CBZ. Analysis of the kinetics profiles of CBZ metabolism with a model that combines the formalism of the Hill equation with an allowance for substrate inhibition demonstrates that the mechanism of interactions of CBZ with CYP3A4 involves multiple substrate-binding events (most likely three). Despite the retention of the multisite binding mechanism in the mutants, functional manifestations reveal an exquisite sensitivity to even minor structural changes in the binding pocket that are introduced by conservative substitutions such as I369F, I369L, and A370V.

Experimental and computational study of BODIPY dye-labeled cavitand dynamics.

Understanding the distance distribution and dynamics between moieties attached to the walls of a resorcin[4]arene cavitand, which is switchable between an expanded kite and a contracted vase form, might enable the use of this molecular system for the study of fundamental distance-dependent interactions. Toward this goal, a combined experimental and molecular dynamics (MD) simulation study on donor/acceptor borondipyrromethene (BODIPY) dye-labeled cavitands present in the vase and kite forms was performed. Direct comparison between anisotropy decays calculated from MD simulations with experimental fluorescence anisotropy data showed excellent agreement, indicating that the simulations provide an accurate representation of the dynamics of the system. Distance distributions between the BODIPY dyes were established by comparing time-resolved Förster resonance energy transfer experiments and MD simulations. Fluorescence intensity decay curves emulated on the basis of the MD trajectories showed good agreement with the experimental data, suggesting that the simulations present an accurate picture of the distance distributions and dynamics in this molecular system and provide an important tool for understanding the behavior of extended molecular systems and designing future applications.

Structure-guided fragment-based in silico drug design of dengue protease inhibitors.

An in silico fragment-based drug design approach was devised and applied towards the identification of small molecule inhibitors of the dengue virus (DENV) NS2B-NS3 protease. Currently, no DENV protease co-crystal structure with bound inhibitor and fully formed substrate binding site is available. Therefore a homology model of DENV NS2B-NS3 protease was generated employing a multiple template spatial restraints method and used for structure-based design. A library of molecular fragments was derived from the ZINC screening database with help of the retrosynthetic combinatorial analysis procedure (RECAP). 150,000 molecular fragments were docked to the DENV protease homology model and the docking poses were rescored using a target-specific scoring function. High scoring fragments were assembled to small molecule candidates by an implicit linking cascade. The cascade included substructure searching and structural filters focusing on interactions with the S1 and S2 pockets of the protease. The chemical space adjacent to the promising candidates was further explored by neighborhood searching. A total of 23 compounds were tested experimentally and two compounds were discovered to inhibit dengue protease (IC(50) = 7.7 µM and 37.9 µM, respectively) and the related West Nile virus protease (IC(50) = 6.3 µM and 39.0 µM, respectively). This study demonstrates the successful application of a structure-guided fragment-based in silico drug design approach for dengue protease inhibitors providing straightforward hit generation using a combination of homology modeling, fragment docking, chemical similarity and structural filters.