Dissecting the CRISPR Cas1-Cas2 Protospacer Binding and Selection Mechanism by Using Molecular Dynamics Simulations

Cas1 and Cas2 are highly conserved proteins among the clustered regularly interspaced short palindromic repeat Cas (CRISPR-Cas) systems and play a crucial role in protospacer selection and integration. According to the double-forked CRISPR Cas1-Cas2 complex, we conducted extensive all-atom molecular dynamics simulations to investigate the protospacer DNA binding and recognition. Our findings revealed that single-point amino acid mutations in Cas1 or in Cas2 had little impact on the binding of the protospacer, both in the binding and precatalytic states. In contrast, multiple-point amino acid mutations, particularly G74A, P80L, and V89A mutations on Cas2 and Cas2' proteins (m-multiple1 system), significantly affected the protospacer binding and selection. Notably, mutations on Cas2 and Cas2' led to an increased number of hydrogen bonds (#HBs) between Cas2&Cas2' and the dsDNA in the m-multiple1 system compared with the wild-type system. And the strong electrostatic interactions between Cas1-Cas2 and the protospacer DNA (psDNA) in the m-multiple1 system again suggested the increase in the binding affinity of protospacer acquisition. Specifically, mutations in Cas2 and Cas2' can remotely make the protospacer adjacent motif complementary (PAMc) sequences better in recognition by the two active sites, while multiple mutations K211E, P202Q, P212L, R138L, V134A, A286T, P282H, and P294H on Cas1a/Cas1b&Cas1a'/Cas1b' (m-multiple2 system) decrease the #HBs and the electrostatic interactions and make the PAMc worse in recognition compared with the wild-type system.