Drugs targeting SARS-CoV-2 could have saved millions of lives during the COVID-19 pandemic, and it is now crucial to develop inhibitors of coronavirus replication in preparation for future outbreaks. We explored two virtual screening strategies to find inhibitors of the SARS-CoV-2 main protease in ultralarge chemical libraries. First, structure-based docking was used to screen a diverse library of 235 million virtual compounds against the active site. One hundred top-ranked compounds were tested in binding and enzymatic assays. Second, a fragment discovered by crystallographic screening was optimized guided by docking of millions of elaborated molecules and experimental testing of 93 compounds. Three inhibitors were identified in the first library screen, and five of the selected fragment elaborations showed inhibitory effects. Crystal structures of target–inhibitor complexes confirmed docking predictions and guided hit-to-lead optimization, resulting in a noncovalent main protease inhibitor with nanomolar affinity, a promising in vitro pharmacokinetic profile, and broad-spectrum antiviral effect in infected cells.

The interaction between HIV-1 protease and 58 structurally diverse transition-state analogue inhibitors has been analyzed by a surface plasmon resonance based biosensor. Association and dissociation rate constants and affinities were determined and displayed as k(on)-k(off)-K(D) maps. It was shown that different classes of inhibitors fall into distinct clusters in these maps. Significant changes in association and dissociation rates were found as a result of modifying the P1/P1' or P2/P2' side chains of a linear lead compound. Similarly, cyclic urea and cyclic sulfamide inhibitors displayed different kinetic features and the affinities of both classes of cyclic compounds were limited by fast dissociation rates. These results confirm that association and dissociation rates are important features of drug-target interactions and indicate that optimization of inhibitor efficacy may be guided by aiming for high association and low dissociation rates rather than high affinity alone. The present approach thus provides a new tool for structure-interaction kinetic analysis and drug discovery.

Ten C2-symmetric cyclic urea and sulfamide derivatives have been synthesized from L-mannonic gamma-lactone and D-mannitol. The results of experimental measurement of their inhibitory potencies against HIV-1 protease were compared to calculated free energies of binding derived from molecular dynamics (MD) simulations. The compounds were selected, firstly, to enable elucidation of the role of stereochemistry for binding affinity (1a-d) and, secondly, to allow evaluation of the effects of variation in the link to the P1 and P1' phenyl groups on affinity (1a and 2-5). Thirdly, compounds with hydrogen bond-accepting or-donating groups attached to the phenyl groups in the P2 and P2' side chains (6 and 7) were selected. Binding free energies were estimated by a linear response method, whose predictive power for estimating binding affinities from MD simulations was demonstrated.