Potent and selective covalent inhibition of the papain-like protease from SARS-CoV-2

Direct-acting antivirals are needed to combat coronavirus disease 2019 (COVID-19), which is caused by severe acute respiratory syndrome-coronavirus-2 (SARS-CoV-2). The papain-like protease (PLpro) domain of Nsp3 from SARS-CoV-2 is essential for viral replication. In addition, PLpro dysregulates the host immune response by cleaving ubiquitin and interferon-stimulated gene 15 protein from host proteins. As a result, PLpro is a promising target for inhibition by small-molecule therapeutics. Here we design a series of covalent inhibitors by introducing a peptidomimetic linker and reactive electrophile onto analogs of the noncovalent PLpro inhibitor GRL0617. The most potent compound inhibits PLpro with kinact/KI = 9,600 M-1 s-1, achieves sub-μM EC50 values against three SARS-CoV-2 variants in mammalian cell lines, and does not inhibit a panel of human deubiquitinases (DUBs) at >30 μM concentrations of inhibitor. An X-ray co-crystal structure of the compound bound to PLpro validates our design strategy and establishes the molecular basis for covalent inhibition and selectivity against structurally similar human DUBs. These findings present an opportunity for further development of covalent PLpro inhibitors.

AI-Accelerated Design of Targeted Covalent Inhibitors for SARS-CoV-2

Direct-acting antivirals for the treatment of the COVID-19 pandemic caused by the SARS-CoV-2 virus are needed to complement vaccination efforts. Given the ongoing emergence of new variants, automated experimentation, and active learning based fast workflows for antiviral lead discovery remain critical to our ability to address the pandemic's evolution in a timely manner. While several such pipelines have been introduced to discover candidates with noncovalent interactions with the main protease (Mpro), here we developed a closed-loop artificial intelligence pipeline to design electrophilic warhead-based covalent candidates. This work introduces a deep learning-assisted automated computational workflow to introduce linkers and an electrophilic "warhead" to design covalent candidates and incorporates cutting-edge experimental techniques for validation. Using this process, promising candidates in the library were screened, and several potential hits were identified and tested experimentally using native mass spectrometry and fluorescence resonance energy transfer (FRET)-based screening assays. We identified four chloroacetamide-based covalent inhibitors of Mpro with micromolar affinities (KI of 5.27 μM) using our pipeline. Experimentally resolved binding modes for each compound were determined using room-temperature X-ray crystallography, which is consistent with the predicted poses. The induced conformational changes based on molecular dynamics simulations further suggest that the dynamics may be an important factor to further improve selectivity, thereby effectively lowering KI and reducing toxicity. These results demonstrate the utility of our modular and data-driven approach for potent and selective covalent inhibitor discovery and provide a platform to apply it to other emerging targets.

Hepatitis C virus NS3/4A inhibitors and other drug-like compounds as covalent binders of SARS-CoV-2 main protease

Severe acute respiratory syndrome-coronavirus 2 (SARS-CoV-2), which causes coronavirus disease 2019 (COVID-19), threatens global public health. The world needs rapid development of new antivirals and vaccines to control the current pandemic and to control the spread of the variants. Among the proteins synthesized by the SARS-CoV-2 genome, main protease (M^pro also known as 3CL^pro) is a primary drug target, due to its essential role in maturation of the viral polyproteins. In this study, we provide crystallographic evidence, along with some binding assay data, that three clinically approved anti hepatitis C virus drugs and two other drug-like compounds covalently bind to the M^pro Cys145 catalytic residue in the active site. Also, molecular docking studies can provide additional insight for the design of new antiviral inhibitors for SARS-CoV-2 using these drugs as lead compounds. One might consider derivatives of these lead compounds with higher affinity to the M^pro as potential COVID-19 therapeutics for further testing and possibly clinical trials.

High-Throughput Virtual Screening and Validation of a SARS-CoV-2 Main Protease Noncovalent Inhibitor

Despite the recent availability of vaccines against the acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the search for inhibitory therapeutic agents has assumed importance especially in the context of emerging new viral variants. In this paper, we describe the discovery of a novel noncovalent small-molecule inhibitor, MCULE-5948770040, that binds to and inhibits the SARS-Cov-2 main protease (Mpro) by employing a scalable high-throughput virtual screening (HTVS) framework and a targeted compound library of over 6.5 million molecules that could be readily ordered and purchased. Our HTVS framework leverages the U.S. supercomputing infrastructure achieving nearly 91% resource utilization and nearly 126 million docking calculations per hour. Downstream biochemical assays validate this Mpro inhibitor with an inhibition constant (Ki) of 2.9 μM (95% CI 2.2, 4.0). Furthermore, using room-temperature X-ray crystallography, we show that MCULE-5948770040 binds to a cleft in the primary binding site of Mpro forming stable hydrogen bond and hydrophobic interactions. We then used multiple μs-time scale molecular dynamics (MD) simulations and machine learning (ML) techniques to elucidate how the bound ligand alters the conformational states accessed by Mpro, involving motions both proximal and distal to the binding site. Together, our results demonstrate how MCULE-5948770040 inhibits Mpro and offers a springboard for further therapeutic design.

Structural, Electronic, and Electrostatic Determinants for Inhibitor Binding to Subsites S1 and S2 in SARS-CoV-2 Main Protease

Creating small-molecule antivirals specific for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) proteins is crucial to battle coronavirus disease 2019 (COVID-19). SARS-CoV-2 main protease (Mpro) is an established drug target for the design of protease inhibitors. We performed a structure-activity relationship (SAR) study of noncovalent compounds that bind in the enzyme's substrate-binding subsites S1 and S2, revealing structural, electronic, and electrostatic determinants of these sites. The study was guided by the X-ray/neutron structure of Mpro complexed with Mcule-5948770040 (compound 1), in which protonation states were directly visualized. Virtual reality-assisted structure analysis and small-molecule building were employed to generate analogues of 1. In vitro enzyme inhibition assays and room-temperature X-ray structures demonstrated the effect of chemical modifications on Mpro inhibition, showing that (1) maintaining correct geometry of an inhibitor's P1 group is essential to preserve the hydrogen bond with the protonated His163; (2) a positively charged linker is preferred; and (3) subsite S2 prefers nonbulky modestly electronegative groups.

Selective class IIa histone deacetylase inhibition via a nonchelating zinc-binding group

In contrast to studies on class I histone deacetylase (HDAC) inhibitors, the elucidation of the molecular mechanisms and therapeutic potential of class IIa HDACs (HDAC4, HDAC5, HDAC7 and HDAC9) is impaired by the lack of potent and selective chemical probes. Here we report the discovery of inhibitors that fill this void with an unprecedented metal-binding group, trifluoromethyloxadiazole (TFMO), which circumvents the selectivity and pharmacologic liabilities of hydroxamates. We confirm direct metal binding of the TFMO through crystallographic approaches and use chemoproteomics to demonstrate the superior selectivity of the TFMO series relative to a hydroxamate-substituted analog. We further apply these tool compounds to reveal gene regulation dependent on the catalytic active site of class IIa HDACs. The discovery of these inhibitors challenges the design process for targeting metalloenzymes through a chelating metal-binding group and suggests therapeutic potential for class IIa HDAC enzyme blockers distinct in mechanism and application compared to current HDAC inhibitors.

Blind docking of pharmaceutically relevant compounds using RosettaLigand

It is difficult to properly validate algorithms that dock a small molecule ligand into its protein receptor using data from the public domain: the predictions are not blind because the correct binding mode is already known, and public test cases may not be representative of compounds of interest such as drug leads. Here, we use private data from a real drug discovery program to carry out a blind evaluation of the RosettaLigand docking methodology and find that its performance is on average comparable with that of the best commercially available current small molecule docking programs. The strength of RosettaLigand is the use of the Rosetta sampling methodology to simultaneously optimize protein sidechain, protein backbone and ligand degrees of freedom; the extensive benchmark test described here identifies shortcomings in other aspects of the protocol and suggests clear routes to improving the method.

Crystal structure of the kinase domain of serum and glucocorticoid-regulated kinase 1 in complex with AMP PNP

Serum and glucocorticoid-regulated kinase 1 (SGK1) is a serine/threonine protein kinase of the AGC family which participates in the control of epithelial ion transport and is implicated in proliferation and apoptosis. We report here the 1.9 A crystal structure of the catalytic domain of inactive human SGK1 in complex with AMP-PNP. SGK1 exists as a dimer formed by two intermolecular disulfide bonds between Cys258 in the activation loop and Cys193. Although most of the SGK1 structure closely resembles the common protein kinase fold, the structure around the active site is unique when compared to most protein kinases. The alphaC helix is not present in this inactive form of SGK1 crystal structure; instead, the segment corresponding to the C helix forms a beta-strand that is stabilized by the N-terminal segment of the activation loop through a short antiparallel beta-sheet. Since the differences from other kinases occur around the ATP binding site, this structure can provide valuable insight into the design of selective and highly potent ATP-competitive inhibitors of SGK1 kinase.

Discovery of Aminofurazan-azabenzimidazoles as Inhibitors of Rho-Kinase with High Kinase Selectivity and Antihypertensive Activity

The discovery, proposed binding mode, and optimization of a novel class of Rho-kinase inhibitors are presented. Appropriate substitution on the 6-position of the azabenzimidazole core provided subnanomolar enzyme potency in vitro while dramatically improving selectivity over a panel of other kinases. Pharmacokinetic data was obtained for the most potent and selective examples and one (6n) has been shown to lower blood pressure in a rat model of hypertension.

A critical assessment of docking programs and scoring functions

Docking is a computational technique that samples conformations of small molecules in protein binding sites; scoring functions are used to assess which of these conformations best complements the protein binding site. An evaluation of 10 docking programs and 37 scoring functions was conducted against eight proteins of seven protein types for three tasks: binding mode prediction, virtual screening for lead identification, and rank-ordering by affinity for lead optimization. All of the docking programs were able to generate ligand conformations similar to crystallographically determined protein/ligand complex structures for at least one of the targets. However, scoring functions were less successful at distinguishing the crystallographic conformation from the set of docked poses. Docking programs identified active compounds from a pharmaceutically relevant pool of decoy compounds; however, no single program performed well for all of the targets. For prediction of compound affinity, none of the docking programs or scoring functions made a useful prediction of ligand binding affinity.

Structure activity relationships of 5-, 6-, and 7-methyl-substituted azepan-3-one cathepsin K inhibitors

The syntheses, in vitro characterizations, and rat and monkey in vivo pharmacokinetic profiles of a series of 5-, 6-, and 7-methyl-substituted azepanone-based cathepsin K inhibitors are described. Depending on the particular regiochemical substitution and stereochemical configuration, methyl-substituted azepanones were identified that had widely varied cathepsin K inhibitory potency as well as pharmacokinetic properties compared to the 4S-parent azepanone analogue, 1 (human cathepsin K, K(i,app) = 0.16 nM, rat oral bioavailability = 42%, rat in vivo clearance = 49.2 mL/min/kg). Of particular note, the 4S-7-cis-methylazepanone analogue, 10, had a K(i,app) = 0.041 nM vs human cathepsin K and 89% oral bioavailability and an in vivo clearance rate of 19.5 mL/min/kg in the rat. Hypotheses that rationalize some of the observed characteristics of these closely related analogues have been made using X-ray crystallography and conformational analysis. These examples demonstrate the potential for modulation of pharmacological properties of cathepsin inhibitors by substituting the azepanone core. The high potency for inhibition of cathepsin K coupled with the favorable rat and monkey pharmacokinetic characteristics of compound 10, also known as SB-462795 or relacatib, has made it the subject of considerable in vivo evaluation for safety and efficacy as an inhibitor of excessive bone resorption in rat, monkey, and human studies, which will be reported elsewhere.

Azepanone-Based Inhibitors of Human Cathepsin L

The extension of a previously reported cathepsin K azepanone-based inhibitor template to the design and synthesis of potent and selective inhibitors of the homologous cysteine protease cathepsin L is detailed. Structure-activity studies examining the effect of inhibitor selectivity as a function of the P3 and P2 binding elements of the potent cathepsin K inhibitor 1 revealed that incorporation of either a P3 quinoline-8-carboxamide or a naphthylene-1-carboxamide led to increased selectivity for cathepsin L over cathepsin K. Substitution of the P2 leucine of 1 with either a phenylalanine or a beta-naphthylalanine also resulted in an increased selectivity for cathepsin L over cathepsin K. Molecular modeling studies with the inhibitors docked within the active sites of both cathepsins L and K have rationalized the observed selectivities. Optimization of cathepsin L binding by the combination of the P3 naphthylene-1-carboxamide with the P2 beta-naphthylalanine provided 15, which is a potent, selective, and competitive inhibitor of human cathepsin L with a K(i) = 0.43 nM.

Indole naphthyridinones as inhibitors of bacterial enoyl-ACP reductases FabI and FabK

Bacterial enoyl-ACP reductase (FabI) is responsible for catalyzing the final step of bacterial fatty acid biosynthesis and is an attractive target for the development of novel antibacterial agents. Previously we reported the development of FabI inhibitor 4 with narrow spectrum antimicrobial activity and in vivo efficacy against Staphylococcus aureus via intraperitoneal (ip) administration. Through iterative medicinal chemistry aided by X-ray crystal structure analysis, a new series of inhibitors has been developed with greatly increased potency against FabI-containing organisms. Several of these new inhibitors have potent antibacterial activity against multidrug resistant strains of S. aureus, and compound 30 demonstrates exceptional oral (po) in vivo efficacy in a S. aureus infection model in rats. While optimizing FabI inhibitory activity, compounds 29 and 30 were identified as having low micromolar FabK inhibitory activity, thereby increasing the antimicrobial spectrum of these compounds to include the FabK-containing pathogens Streptococcus pneumoniae and Enterococcus faecalis. The results described herein support the hypothesis that bacterial enoyl-ACP reductases are valid targets for antibacterial agents.

First X-ray cocrystal structure of a bacterial FabH condensing enzyme and a small molecule inhibitor achieved using rational design and homology modeling

The first cocrystal structure of a bacterial FabH condensing enzyme and a small molecule inhibitor is reported. The inhibitor was obtained by rational modification of a high throughput screening lead with the aid of a S. pneumoniae FabH homology model. This homology model was used to design analogues that would have both high affinity for the enzyme and appropriate aqueous solubility to facilitate cocrystallization studies.

Defining and combating the mechanisms of triclosan resistance in clinical isolates of Staphylococcus aureus

The MICs of triclosan for 31 clinical isolates of Staphylococcus aureus were 0.016 micro g/ml (24 strains), 1 to 2 micro g/ml (6 strains), and 0.25 micro g/ml (1 strain). All the strains for which triclosan MICs were elevated (>0.016 micro g/ml) showed three- to fivefold increases in their levels of enoyl-acyl carrier protein (ACP) reductase (FabI) production. Furthermore, strains for which triclosan MICs were 1 to 2 micro g/ml overexpressed FabI with an F204C alteration. Binding studies with radiolabeled NAD(+) demonstrated that this change prevents the formation of the stable triclosan-NAD(+)-FabI complex, and both this alteration and its overexpression contributed to achieving MICs of 1 to 2 micro g/ml for these strains. Three novel, potent inhibitors of FabI (50% inhibitory concentrations, < or =64 nM) demonstrated up to 1,000-fold better activity than triclosan against the strains for which triclosan MICs were elevated. None of the compounds tested from this series formed a stable complex with NAD(+)-FabI. Consequently, although the overexpression of wild-type FabI gave rise to an increase in the MICs, as expected, overexpression of FabI with an F204C alteration did not cause an additional increase in resistance. Therefore, this work identifies the mechanisms of triclosan resistance in S. aureus, and we present three compounds from a novel chemical series of FabI inhibitors which have excellent activities against both triclosan-resistant and -sensitive clinical isolates of S. aureus.

Discovery of a novel and potent class of FabI-directed antibacterial agents

Bacterial enoyl-acyl carrier protein (ACP) reductase (FabI) catalyzes the final step in each elongation cycle of bacterial fatty acid biosynthesis and is an attractive target for the development of new antibacterial agents. High-throughput screening of the Staphylococcus aureus FabI enzyme identified a novel, weak inhibitor with no detectable antibacterial activity against S. aureus. Iterative medicinal chemistry and X-ray crystal structure-based design led to the identification of compound 4 [(E)-N-methyl-N-(2-methyl-1H-indol-3-ylmethyl)-3-(7-oxo-5,6,7,8-tetrahydro-1,8-naphthyridin-3-yl)acrylamide], which is 350-fold more potent than the original lead compound obtained by high-throughput screening in the FabI inhibition assay. Compound 4 has exquisite antistaphylococci activity, achieving MICs at which 90% of isolates are inhibited more than 500 times lower than those of nine currently available antibiotics against a panel of multidrug-resistant strains of S. aureus and Staphylococcus epidermidis. Furthermore, compound 4 exhibits excellent in vivo efficacy in an S. aureus infection model in rats. Biochemical and genetic approaches have confirmed that the mode of antibacterial action of compound 4 and related compounds is via inhibition of FabI. Compound 4 also exhibits weak FabK inhibitory activity, which may explain its antibacterial activity against Streptococcus pneumoniae and Enterococcus faecalis, which depend on FabK and both FabK and FabI, respectively, for their enoyl-ACP reductase function. These results show that compound 4 is representative of a new, totally synthetic series of antibacterial agents that has the potential to provide novel alternatives for the treatment of S. aureus infections that are resistant to our present armory of antibiotics.

Discovery of aminopyridine-based inhibitors of bacterial enoyl-ACP reductase (FabI)

Bacterial enoyl-ACP reductase (FabI) catalyzes the final step in each cycle of bacterial fatty acid biosynthesis and is an attractive target for the development of new antibacterial agents. Our efforts to identify potent, selective FabI inhibitors began with screening of the GlaxoSmithKline proprietary compound collection, which identified several small-molecule inhibitors of Staphylococcus aureus FabI. Through a combination of iterative medicinal chemistry and X-ray crystal structure based design, one of these leads was developed into the novel aminopyridine derivative 9, a low micromolar inhibitor of FabI from S. aureus (IC(50) = 2.4 microM) and Haemophilus influenzae (IC(50) = 4.2 microM). Compound 9 has good in vitro antibacterial activity against several organisms, including S. aureus (MIC = 0.5 microg/mL), and is effective in vivo in a S. aureus groin abscess infection model in rats. Through FabI overexpressor and macromolecular synthesis studies, the mode of action of 9 has been confirmed to be inhibition of fatty acid biosynthesis via inhibition of FabI. Taken together, these results support FabI as a valid antibacterial target and demonstrate the potential of small-molecule FabI inhibitors for the treatment of bacterial infections.

Structure-based combinatorial library design: discovery of non-peptidic inhibitors of caspases 3 and 8

Structure-based design of a combinatorial array was carried out in order to identify non-peptidic thiomethylketone inhibitors of caspases 3 and 8. Five compounds from the designed array were active against caspase 3, and two were active against caspase 8. A 2.5-A resolution co-crystal structure of caspase 3 and a thiomethylketone array member is reported. The structure-based design strategy has proved useful for identifying caspase inhibitors.

Potent and selective nonpeptide inhibitors of caspases 3 and 7 inhibit apoptosis and maintain cell functionality

Caspases have been strongly implicated to play an essential role in apoptosis. A critical question regarding the role(s) of these proteases is whether selective inhibition of an effector caspase(s) will prevent cell death. We have identified potent and selective non-peptide inhibitors of the effector caspases 3 and 7. The inhibition of apoptosis and maintenance of cell functionality with a caspase 3/7-selective inhibitor is demonstrated for the first time, and suggests that targeting these two caspases alone is sufficient for blocking apoptosis. Furthermore, an x-ray co-crystal structure of the complex between recombinant human caspase 3 and an isatin sulfonamide inhibitor has been solved to 2.8-A resolution. In contrast to previously reported peptide-based caspase inhibitors, the isatin sulfonamides derive their selectivity for caspases 3 and 7 by interacting primarily with the S(2) subsite, and do not bind in the caspase primary aspartic acid binding pocket (S(1)). These inhibitors blocked apoptosis in murine bone marrow neutrophils and human chondrocytes. Furthermore, in camptothecin-induced chondrocyte apoptosis, cell functionality as measured by type II collagen promoter activity is maintained, an activity considered essential for cartilage homeostasis. These data suggest that inhibiting chondrocyte cell death with a caspase 3/7-selective inhibitor may provide a novel therapeutic approach for the prevention and treatment of osteoarthritis, or other disease states characterized by excessive apoptosis.

Nucleic acid base-pairing and N-methylacetamide self-association in chloroform: affinity and conformation

A recently developed computational method, 'mining minima', is used to examine the hydrogen-bonding interactions of nucleic acid base-pairs and of the N-methylacetamide homodimer in chloroform. The mining minima algorithm aggressively samples molecular conformations, identifies the most important local minima, and computes their contributions to the overall free energy of the system. Here, the CHARMM 98 parameter set is used for the potential energy and the generalized Born/surface area solvent model is used to account for the influence of the solvent. Good agreement with experiment is obtained for the non-covalent binding affinities of a series of complexes. The computational approach used here is applicable to a range of molecular systems.

Thermodynamic linkage between the binding of protons and inhibitors to HIV-1 protease

The aspartyl dyad of free HIV-1 protease has apparent pK(a)s of approximately 3 and approximately 6, but recent NMR studies indicate that the aspartyl dyad is fixed in the doubly protonated form over a wide pH range when cyclic urea inhibitors are bound, and in the monoprotonated form when the inhibitor KNI-272 is bound. We present computations and measurements related to these changes in protonation and to the thermodynamic linkage between protonation and inhibition. The Poisson-Boltzmann model of electrostatics is used to compute the apparent pK(a)s of the aspartyl dyad in the free enzyme and in complexes with four different inhibitors. The calculations are done with two parameter sets. One assigns epsilon = 4 to the solute interior and uses a detailed model of ionization; the other uses epsilon = 20 for the solute interior and a simplified representation of ionization. For the free enzyme, both parameter sets agree well with previously measured apparent pK(a)s of approximately 3 and approximately 6. However, the calculations with an internal dielectric constant of 4 reproduce the large pKa shifts upon binding of inhibitors, but the calculations with an internal dielectric constant of 20 do not. This observation has implications for the accurate calculation of pK(a)s in complex protein environments. Because binding of a cyclic urea inhibitor shifts the pK(a)s of the aspartyl dyad, changing the pH is expected to change its apparent binding affinity. However, we find experimentally that the affinity is independent of pH from 5.5 to 7.0. Possible explanations for this discrepancy are discussed.

A new class of models for computing receptor-ligand binding affinities

Models for predicting the binding affinities of molecules in solution are either very detailed, making them computationally intensive and hard to test, or very simple, and thus less informative than one might wish. A new class of models that focus on the predominant states of the binding molecules promise to capture the essential physics of binding at modest computational cost.