Inhibition of the hexamerization of SARS-CoV-2 endoribonuclease and modeling of RNA structures bound to the hexamer

Non-structural protein 15 (Nsp15) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) forms a homo hexamer and functions as an endoribonuclease. Here, we propose that Nsp15 activity may be inhibited by preventing its hexamerization through drug binding. We first explored the stable conformation of the Nsp15 monomer as the global free energy minimum conformation in the free energy landscape using a combination of parallel cascade selection molecular dynamics (PaCS-MD) and the Markov state model (MSM), and found that the Nsp15 monomer forms a more open conformation with larger druggable pockets on the surface. Targeting the pockets with high druggability scores, we conducted ligand docking and identified compounds that tightly bind to the Nsp15 monomer. The top poses with Nsp15 were subjected to binding free energy calculations by dissociation PaCS-MD and MSM (dPaCS-MD/MSM), indicating the stability of the complexes. One of the identified pockets, which is distinctively bound by inosine analogues, may be an alternative binding site to stabilize viral RNA binding and/or an alternative catalytic site. We constructed a stable RNA structure model bound to both UTP and alternative binding sites, providing a reasonable proposed model of the Nsp15/RNA complex.