Cas9-catalyzed DNA Cleavage Generates Staggered Ends: Evidence from Molecular Dynamics Simulations

The CRISPR-associated endonuclease Cas9 from Streptococcus pyogenes (spCas9) along with a single guide RNA (sgRNA) has emerged as a versatile toolbox for genome editing. Despite recent advances in the mechanism studies on spCas9-sgRNA-mediated double-stranded DNA (dsDNA) recognition and cleavage, it is still unclear how the catalytic Mg2+ ions induce the conformation changes toward the catalytic active state. It also remains controversial whether Cas9 generates blunt-ended or staggered-ended breaks with overhangs in the DNA. To investigate these issues, here we performed the first all-atom molecular dynamics simulations of the spCas9-sgRNA-dsDNA system with and without Mg2+ bound. The simulation results showed that binding of two Mg2+ ions at the RuvC domain active site could lead to structurally and energetically favorable coordination ready for the non-target DNA strand cleavage. Importantly, we demonstrated with our simulations that Cas9-catalyzed DNA cleavage produces 1-bp staggered ends rather than generally assumed blunt ends.

[1]  R. Barrangou,et al.  CRISPR Provides Acquired Resistance Against Viruses in Prokaryotes , 2007, Science.

[2]  Wenning Wang,et al.  Insights into the Inhibitory Mechanism of D13-9001 to the Multidrug Transporter AcrB through Molecular Dynamics Simulations. , 2016, The journal of physical chemistry. B.

[3]  Jennifer A. Doudna,et al.  Conformational control of DNA target cleavage by CRISPR–Cas9 , 2015, Nature.

[4]  J. Andrew McCammon,et al.  Striking Plasticity of CRISPR-Cas9 and Key Role of Non-target DNA, as Revealed by Molecular Simulations , 2016, ACS central science.

[5]  R. Barrangou,et al.  Cas9–crRNA ribonucleoprotein complex mediates specific DNA cleavage for adaptive immunity in bacteria , 2012, Proceedings of the National Academy of Sciences.

[6]  Jennifer A. Doudna,et al.  Structures of Cas9 Endonucleases Reveal RNA-Mediated Conformational Activation , 2014, Science.

[7]  Ross C. Walker,et al.  An overview of the Amber biomolecular simulation package , 2013 .

[8]  T. Darden,et al.  Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems , 1993 .

[9]  David R. Liu,et al.  Small Molecule-Triggered Cas9 Protein with Improved Genome-Editing Specificity , 2015, Nature chemical biology.

[10]  Philippe Horvath,et al.  The CRISPR/Cas bacterial immune system cleaves bacteriophage and plasmid DNA , 2010, Nature.

[11]  Pengfei Li,et al.  Rational Design of Particle Mesh Ewald Compatible Lennard-Jones Parameters for +2 Metal Cations in Explicit Solvent. , 2013, Journal of chemical theory and computation.

[12]  Marcello Maresca,et al.  Obligate Ligation-Gated Recombination (ObLiGaRe): Custom-designed nuclease-mediated targeted integration through nonhomologous end joining , 2013, Genome research.

[13]  James E. DiCarlo,et al.  RNA-Guided Human Genome Engineering via Cas9 , 2013, Science.

[14]  L. Marraffini,et al.  CRISPR interference: RNA-directed adaptive immunity in bacteria and archaea , 2010, Nature Reviews Genetics.

[15]  L Wang,et al.  Molecular dynamics and free-energy calculations applied to affinity maturation in antibody 48G7. , 1999, Proceedings of the National Academy of Sciences of the United States of America.

[16]  Wei Yang,et al.  Catalytic mechanism of RNA backbone cleavage by ribonuclease H from quantum mechanics/molecular mechanics simulations. , 2011, Journal of the American Chemical Society.

[17]  Wei Yang,et al.  Stepwise analyses of metal ions in RNase H catalysis from substrate destabilization to product release , 2006, The EMBO journal.

[18]  P. Kollman,et al.  Settle: An analytical version of the SHAKE and RATTLE algorithm for rigid water models , 1992 .

[19]  Le Cong,et al.  Multiplex Genome Engineering Using CRISPR/Cas Systems , 2013, Science.

[20]  Wei Yang,et al.  Crystal Structures of RNase H Bound to an RNA/DNA Hybrid: Substrate Specificity and Metal-Dependent Catalysis , 2005, Cell.

[21]  C. Simmerling,et al.  ff14SB: Improving the Accuracy of Protein Side Chain and Backbone Parameters from ff99SB. , 2015, Journal of chemical theory and computation.

[22]  Feng Zhang,et al.  Crystal Structure of Cas9 in Complex with Guide RNA and Target DNA , 2014, Cell.

[23]  P A Kollman,et al.  Molecular dynamics and continuum solvent studies of the stability of polyG-polyC and polyA-polyT DNA duplexes in solution. , 1998, Journal of biomolecular structure & dynamics.

[24]  J. Joung,et al.  Defining and improving the genome-wide specificities of CRISPR–Cas9 nucleases , 2016, Nature Reviews Genetics.

[25]  Lei S. Qi,et al.  Small molecules enhance CRISPR genome editing in pluripotent stem cells. , 2015, Cell stem cell.

[26]  A. Regev,et al.  Cpf1 Is a Single RNA-Guided Endonuclease of a Class 2 CRISPR-Cas System , 2015, Cell.

[27]  Wei Yang,et al.  An equivalent metal ion in one- and two-metal-ion catalysis , 2008, Nature Structural &Molecular Biology.

[28]  Gilbert Chu,et al.  Processing of DNA for nonhomologous end‐joining by cell‐free extract , 2005, The EMBO journal.

[29]  C. Desmaze,et al.  Impact of the KU80 pathway on NHEJ-induced genome rearrangements in mammalian cells. , 2004, Molecular cell.

[30]  Jennifer A. Doudna,et al.  Structures of a CRISPR-Cas9 R-loop complex primed for DNA cleavage , 2016, Science.

[31]  Holger Gohlke,et al.  MMPBSA.py: An Efficient Program for End-State Free Energy Calculations. , 2012, Journal of chemical theory and computation.

[32]  J. Doudna,et al.  RNA-guided genetic silencing systems in bacteria and archaea , 2012, Nature.

[33]  K Schulten,et al.  VMD: visual molecular dynamics. , 1996, Journal of molecular graphics.

[34]  Wei Yang,et al.  Nucleases: diversity of structure, function and mechanism , 2010, Quarterly Reviews of Biophysics.

[35]  Osamu Nureki,et al.  Structural Basis for the Altered PAM Specificities of Engineered CRISPR-Cas9. , 2016, Molecular cell.

[36]  Jae Young Lee,et al.  Making and breaking nucleic acids: two-Mg2+-ion catalysis and substrate specificity. , 2006, Molecular cell.

[37]  P. Kollman,et al.  Calculating structures and free energies of complex molecules: combining molecular mechanics and continuum models. , 2000, Accounts of chemical research.

[38]  J. Doudna,et al.  A Programmable Dual-RNA–Guided DNA Endonuclease in Adaptive Bacterial Immunity , 2012, Science.

[39]  K. Merz,et al.  Systematic Parameterization of Monovalent Ions Employing the Nonbonded Model. , 2015, Journal of chemical theory and computation.

[40]  M. Jinek,et al.  Structural basis of PAM-dependent target DNA recognition by the Cas9 endonuclease , 2014, Nature.

[41]  Jennifer A. Doudna,et al.  A Cas9–guide RNA complex preorganized for target DNA recognition , 2015, Science.

[42]  Youyong Li,et al.  Assessing the performance of MM/PBSA and MM/GBSA methods. 4. Accuracies of MM/PBSA and MM/GBSA methodologies evaluated by various simulation protocols using PDBbind data set. , 2014, Physical chemistry chemical physics : PCCP.

[43]  M. Jinek,et al.  Structural Plasticity of PAM Recognition by Engineered Variants of the RNA-Guided Endonuclease Cas9. , 2016, Molecular cell.

[44]  Tingjun Hou,et al.  Assessing the Performance of the MM/PBSA and MM/GBSA Methods. 1. The Accuracy of Binding Free Energy Calculations Based on Molecular Dynamics Simulations , 2011, J. Chem. Inf. Model..