The clustered regularly interspaced short palindromic
repeats (CRISPR)
technology is an RNA-guided targeted genome-editing tool using Cas
family proteins. Two magnesium-dependent nuclease domains of the
Cas9 enzyme, termed HNH and RuvC, are responsible for cleaving the
target DNA (t-DNA) and nontarget DNA strands, respectively. The HNH
domain is believed to determine the DNA cleavage activity of both
endonuclease domains and is sensitive to complementary RNA-DNA base
pairing. However, the underlying molecular mechanisms of CRISPR-Cas9,
by which it rebukes or accepts mismatches, are poorly understood.
Thus, investigation of the structure and dynamics of the catalytic
state of Cas9 with either matched or mismatched t-DNA can provide
insights into improving its specificity by reducing off-target cleavages.
Here, we focus on a recently discovered catalytic-active form of the Streptococcus pyogenes Cas9 (SpCas9) and employ classical
molecular dynamics and coupled quantum mechanics/molecular mechanics
simulations to study two possible mechanisms of t-DNA cleavage reaction
catalyzed by the HNH domain. Moreover, by designing a mismatched t-DNA
structure called MM5 (C to G at the fifth position from the protospacer
adjacent motif region), the impact of single-guide RNA (sgRNA) and
t-DNA complementarity on the catalysis process was investigated. Based
on these simulations, our calculated binding affinities, minimum energy
paths, and analysis of catalytically important residues provide atomic-level
details of the differences between matched and mismatched cleavage
reactions. In addition, several residues exhibit significant differences
in their catalytic roles for the two studied systems, including K253,
K263, R820, K896, and K913.