Abstract
A comprehensive understanding of the atomic level mechanism governing the binding nonstructural protein 1 of SARS-CoV-1 (SARS-CoV-1 NSP1) and SARS-CoV-2 (SARS-CoV-2 NSP1) to Pol α-primase is important to advance the development of small molecule inhibitors for the treatment COVID-19. In this study, we use both all-atom steered molecular dynamics (all-atom SMD) and coarse-grained umbrella sampling (coarse-grained US) simulations to assess the binding affinity of SARS-CoV-1 NSP1 and SARS-CoV-2 NSP1 to Pol α-primase. Our all-atom SMD and coarse-grained US simulations consistently indicate that SARS-CoV-2 NSP1 exhibits stronger affinity for Pol α-primase compared to SARS-CoV-1 NSP1, implying that SARS-CoV-2 poses a greater risk than SARS-CoV-1 in impeding DNA replication for DNA synthesis. Through an energetic decomposition analysis of the interaction energies within these complexes, we identify electrostatic interactions as the primary contributors to the observed difference in binding affinity. We found that hydrogen bonds between Asp33 and Arg616 in SARS1 NSP1-Pol α-primase, and Asp33 with Arg616 and Lys655 in SARS2 NSP1-Pol α-primase, are critical for the interaction of both SARS-CoV-1 NSP1 and SARS-CoV-2 NSP1 with Pol α-primase. Asp33 in SARS-CoV-2 NSP1 shows increased solubility and stability compared to SARS-CoV-1 NSP1, enhancing its association with Pol α-primase. This finding lays the groundwork for innovative strategies aimed at inhibiting the interaction between these entities, offering promising avenues for therapeutic intervention against COVID-19. We also estimated the binding free energy of DNA to Pol α-primase, SARS1 NSP1-Pol α-primase, and SARS2 NSP1-Pol α-primase using the MM-PBSA method. The results show the order: Pol α-primase-DNA < SARS1 NSP1-Pol α-primase-DNA < SARS2 NSP1-Pol α-primase-DNA, indicating that both SARS-CoV-1 NSP1 and SARS-CoV-2 NSP1 reduce DNA binding to Pol α-primase, suggesting impaired DNA synthesis.