Abstract
CRISPR-Cas9 is a powerful tool for target genome editing in living cells. Significant advances have been made to understand how this system cleaves target DNA. HNH is a nuclease domain, which shares structural similarity with the HNH endonuclease characterzied by a beta-beta-alpha-metal fold. Therefore, based on one- and two-metal-ion mechanisms, homology modeling and molecular dynamics (MD) simulation are suitable tools for building an atomic model of Cas9 in the DNA cleavage state. Here, by modeling and MD, we presented an atomic model of SpCas9-sgRNA-DNA complex with the cleavage state. This model shows that the HNH and RuvC conformations resemble their DNA cleavage state where the active-sites in the complex coordinate with DNA, Mg2+ ions, and water. Among them, residues D10, E762, H983, and D986 locate at the first shell of the RuvC active-site and interact with the ions directly, residues H982 or/and H985 are general (Lewis) bases, and the coordinated water is located at the positions for nucleophilic attack of the scissile phosphate. Meanwhile, this catalytic model led us to engineer a new SpCas9 variant (SpCas9-H982A + H983D) with reduced off-target effects. Thus, our study provided new mechanistic insights into the CRISPR-Cas9 system in the DNA cleavage state and offered useful guidance for engineering new CRISPR-Cas9 editing systems with improved specificity.