Spontaneous Embedding of DNA Mismatches Within the RNA:DNA Hybrid of CRISPR-Cas9

Secondary Conformational Checkpoint in CRISPR-Cas9

A specific checkpoint between target DNA binding and cleavage primarily governs the precision of Cas9 gene editing. Although various CRISPR-Cas9 variants have been developed to improve DNA cleavage accuracy, we still lack a comprehensive understanding of how they work at the molecular level. Herein, we have focused on studying the late-stage conformational transitions of Cas9 and an evolved Cas9 mutant (evoCas9) that start from the precleavage state. Our submilliseconds of dynamic simulations reveal that the presence of base mismatches leads the HNH nuclease domain of Cas9 to alter its principal functional modes of motion, thereby impairing its conformational activation. This observation suggests the existence of a secondary conformational checkpoint that fine-tunes the final DNA cleavage activation. Remarkably, evoCas9 is prone to deviating from the normal activation pathway with base mismatches. This is characterized by a noticeable shift in the positioning of the HNH domain and a significantly perturbed allosteric communication network within the enzyme. Therefore, the mutations evolved in evoCas9 also reinforce the secondary checkpoint in addition to the previously identified primary checkpoint, collectively ensuring this variant's high gene-editing accuracy. This mechanism should also apply to other Cas9-guide RNA variants with enhanced fidelity.

Effect of Cas9 Protein on the Seed-Target Base Pair of the sgRNA/DNA Hybrid Duplex

CRISPR (clustered regularly interspaced short palindromic repeats)/CRISPR-associated protein (Cas9) has been widely used for gene editing. Not all guide RNAs can cleave the DNA efficiently remains a major challenge to CRISPR/Cas9-mediated genome engineering. Therefore, understanding how the Cas9 complex successfully and efficiently identifies specific functional targets through base-pairing has great implications for such applications. The 10-nt seed sequence at the 3' end of the guide RNA is critical to target recognition and cleavage. Here, through stretching molecular dynamics simulation, we studied the thermodynamics and kinetics of the binding-dissociation process of the seed base and the target DNA base with the Cas9 protein. The results showed that in the presence of Cas9 protein, the enthalpy change and entropy change in binding-dissociation of the seed base with the target are smaller than those without the Cas9 protein. The reduction of entropy penalty upon association with the protein resulted from the pre-organization of the seed base in an A-form helix, and the reduction of enthalpy change was due to the electrostatic attraction of the positively charged channel with the negative target DNA. The binding barrier coming from the entropy loss and the dissociation barrier resulting from the destruction of the base pair in the presence of Cas9 protein were lower than those without protein, which indicates that the seed region is crucial for efficiently searching the correct target by accelerating the binding rate and dissociating fast from the wrong target.

Hybrid Multitask Learning Reveals Sequence Features Driving Specificity in the CRISPR/Cas9 System

CRISPR/Cas9 technology is capable of precisely editing genomes and is at the heart of various scientific and medical advances in recent times. The advances in biomedical research are hindered because of the inadvertent burden on the genome when genome editors are employed—the off-target effects. Although experimental screens to detect off-targets have allowed understanding the activity of Cas9, that knowledge remains incomplete as the rules do not extrapolate well to new target sequences. Off-target prediction tools developed recently have increasingly relied on machine learning and deep learning techniques to reliably understand the complete threat of likely off-targets because the rules that drive Cas9 activity are not fully understood. In this study, we present a count-based as well as deep-learning-based approach to derive sequence features that are important in deciding on Cas9 activity at a sequence. There are two major challenges in off-target determination—the identification of a likely site of Cas9 activity and the prediction of the extent of Cas9 activity at that site. The hybrid multitask CNN–biLSTM model developed, named CRISP–RCNN, simultaneously predicts off-targets and the extent of activity on off-targets. Employing methods of integrated gradients and weighting kernels for feature importance approximation, analysis of nucleotide and position preference, and mismatch tolerance have been performed.

Genetic engineering and genome editing in plants, animals and humans: Facts and myths

Human history is inextricably linked to the introduction of desirable heritable traits in plants and animals. Selective breeding (SB) predates our historical period and has been practiced since the advent of agriculture and farming more than ten thousand years ago. Since the 1970s, methods of direct plant and animal genome manipulation are constantly being developed. These are collectively described as "genetic engineering" (GE). Plant GE aims to improve nutritional value, insect resistance and weed control. Animal GE has focused on livestock improvement and disease control. GE applications also involve medical improvements intended to treat human disease. The scientific consensus built around marketed products of GE organisms (GEOs) is usually well established, noting significant benefits and low risks. GEOs are exhaustively scrutinized in the EU and many non-EU countries for their effects on human health and the environment, but scrutiny should be equally applied to all previously untested organisms derived directly from nature or through selective breeding. In fact, there is no evidence to suggest that natural or selectively bred plants and animals are in principle safer to humans than GEOs. Natural and selectively bred strains evolve over time via genetic mutations that can be as risky to humans and the environment as the mutations found in GEOs. Thus, previously untested plant and animal strains aimed for marketing should be proven useful or harmful to humans only upon comparative testing, regardless of their origin. Highlighting the scientific consensus declaring significant benefits and rather manageable risks provided by equitably accessed GEOs, can mitigate negative predispositions by policy makers and the public. Accordingly, we provide an overview of the underlying technologies and the scientific consensus to help resolve popular myths about the safety and usefulness of GEOs.

The Electronic Structure of Genome Editors from the First Principles

Genome editing based on the CRISPR-Cas9 system has paved new avenues for medicine, pharmaceutics, biotechnology, and beyond. This article reports the role of first-principles (ab-initio) molecular dynamics (MD) in the CRISPR-Cas9 revolution, achieving a profound understanding of the enzymatic function and offering valuable insights for enzyme engineering. We introduce the methodologies and explain the use of ab-initio MD simulations to characterize the two-metal dependent mechanism of DNA cleavage in the RuvC domain of the Cas9 enzyme, and how a second catalytic domain, HNH, cleaves the target DNA with the aid of a single metal ion. A detailed description of how ab-initio MD is combined with free-energy methods - i.e., thermodynamic integration and metadynamics - to break and form chemical bonds is given, explaining the use of these methods to determine the chemical landscape and establish the catalytic mechanism in CRISPR-Cas9. The critical role of classical methods is also discussed, explaining theory and application of constant pH MD simulations, used to accurately predict the catalytic residues' protonation states. Overall, first-principles methods are shown to unravel the electronic structure of the Cas9 enzyme, providing valuable insights that can serve for the design of genome editing tools with improved catalytic efficiency or controllable activity.

Insights into the Mechanism of CRISPR/Cas9-Based Genome Editing from Molecular Dynamics Simulations

The CRISPR/Cas9 system is a popular genome-editing tool with immense therapeutic potential. It is a simple two-component system (Cas9 protein and RNA) that recognizes the DNA sequence on the basis of RNA:DNA complementarity, and the Cas9 protein catalyzes the double-stranded break in the DNA. In the past decade, near-atomic resolution structures at various stages of the CRISPR/Cas9 DNA editing pathway have been reported along with numerous experimental and computational studies. Such studies have boosted knowledge of the genome-editing mechanism. Despite such advancements, the application of CRISPR/Cas9 in therapeutics is still limited, primarily due to off-target effects. Several studies aim at engineering high-fidelity Cas9 to minimize the off-target effects. Molecular Dynamics (MD) simulations have been an excellent complement to the experimental studies for investigating the mechanism of CRISPR/Cas9 editing in terms of structure, thermodynamics, and kinetics. MD-based studies have uncovered several important molecular aspects of Cas9, such as nucleotide binding, catalytic mechanism, and off-target effects. In this Review, the contribution of MD simulation to understand the CRISPR/Cas9 mechanism has been discussed, preceded by an overview of the history, mechanism, and structural aspects of the CRISPR/Cas9 system. These studies are important for the rational design of highly specific Cas9 and will also be extremely promising for achieving more accurate genome editing in the future.

Natural Nucleoside Modifications in Guide RNAs Can Modulate the Activity of the CRISPR-Cas9 System In Vitro

At the present time, the clustered regularly interspaced short palindromic repeats (CRISPR)-CRISPR-associated protein 9 (Cas9) system has been widely adopted as an efficient genomic editing tool. However, there are some actual problems such as the off-target effects, cytotoxicity, and immunogenicity. The incorporation of modifications into guide RNAs permits enhancing both the efficiency and the specificity of the CRISPR-Cas9 system. In this study, we demonstrate that the inclusion of N⁶-methyladenosine, 5-methylcytidine, and pseudouridine in trans-activating RNA (tracrRNA) or in single guide RNA (sgRNA) enables efficient gene editing in vitro. We found that the complexes of modified guide RNAs with Cas9 protein promoted cleavage of the target short/long duplexes and plasmid substrates. In addition, the modified monomers in guide RNAs allow increasing the specificity of CRISPR-Cas9 system in vitro and promote diminishing both the immunostimulating and the cytotoxic effects of sgRNAs.

An overview of structural approaches to study therapeutic RNAs

RNAs provide considerable opportunities as therapeutic agent to expand the plethora of classical therapeutic targets, from extracellular and surface proteins to intracellular nucleic acids and its regulators, in a wide range of diseases. RNA versatility can be exploited to recognize cell types, perform cell therapy, and develop new vaccine classes. Therapeutic RNAs (aptamers, antisense nucleotides, siRNA, miRNA, mRNA and CRISPR-Cas9) can modulate or induce protein expression, inhibit molecular interactions, achieve genome editing as well as exon-skipping. A common RNA thread, which makes it very promising for therapeutic applications, is its structure, flexibility, and binding specificity. Moreover, RNA displays peculiar structural plasticity compared to proteins as well as to DNA. Here we summarize the recent advances and applications of therapeutic RNAs, and the experimental and computational methods to analyze their structure, by biophysical techniques (liquid-state NMR, scattering, reactivity, and computational simulations), with a focus on dynamic and flexibility aspects and to binding analysis. This will provide insights on the currently available RNA therapeutic applications and on the best techniques to evaluate its dynamics and reactivity.

Structural basis for Cas9 off-target activity

The target DNA specificity of the CRISPR-associated genome editor nuclease Cas9 is determined by complementarity to a 20-nucleotide segment in its guide RNA. However, Cas9 can bind and cleave partially complementary off-target sequences, which raises safety concerns for its use in clinical applications. Here, we report crystallographic structures of Cas9 bound to bona fide off-target substrates, revealing that off-target binding is enabled by a range of noncanonical base-pairing interactions within the guide:off-target heteroduplex. Off-target substrates containing single-nucleotide deletions relative to the guide RNA are accommodated by base skipping or multiple noncanonical base pairs rather than RNA bulge formation. Finally, PAM-distal mismatches result in duplex unpairing and induce a conformational change in the Cas9 REC lobe that perturbs its conformational activation. Together, these insights provide a structural rationale for the off-target activity of Cas9 and contribute to the improved rational design of guide RNAs and off-target prediction algorithms.

Deciphering the QR Code of the CRISPR-Cas9 System: Synergy between Gln768 (Q) and Arg976 (R)

Markov state models (MSMs) and machine learning (ML) algorithms can extrapolate the long-time-scale behavior of large biomolecules from molecular dynamics (MD) trajectories. In this study, an MD-MSM-ML scheme has been applied to probe the large endonuclease (Cas9) in the bacterial adaptive immunity CRISPR-Cas9 system. CRISPR has become a programmable and state-of-the-art powerful genome editing tool that has already revolutionized life sciences. CRISPR-Cas9 is programmed to process specific DNA sequences in the genome. However, human/biomedical applications are compromised by off-target DNA damage. Characterization of Cas9 at the structural and biophysical levels is a prerequisite for the development of efficient and high-fidelity Cas9 variants. The Cas9 wild type and two variants (R63A-R66A-R70A, R69A-R71A-R74A-R78A) are studied herein. The configurational space of Cas9 is provided with a focus on the conformations of the side chains of two residues (Gln768 and Arg976). A model for the synergy between those two residues is proposed. The results are discussed within the context of experimental literature. The results and methodology can be exploited for the study of large biomolecules in general and for the engineering of more efficient and safer Cas9 variants for applications.

A Machine Learning Approach to Identify the Importance of Novel Features for CRISPR/Cas9 Activity Prediction

The reprogrammable CRISPR/Cas9 genome editing tool's growing popularity is hindered by unwanted off-target effects. Efforts have been directed toward designing efficient guide RNAs as well as identifying potential off-target threats, yet factors that determine efficiency and off-target activity remain obscure. Based on sequence features, previous machine learning models performed poorly on new datasets, thus there is a need for the incorporation of novel features. The binding energy estimation of the gRNA-DNA hybrid as well as the Cas9-gRNA-DNA hybrid allowed generating better performing machine learning models for the prediction of Cas9 activity. The analysis of feature contribution towards the model output on a limited dataset indicated that energy features played a determining role along with the sequence features. The binding energy features proved essential for the prediction of on-target activity and off-target sites. The plateau, in the performance on unseen datasets, of current machine learning models could be overcome by incorporating novel features, such as binding energy, among others. The models are provided on GitHub (GitHub Inc., San Francisco, CA, USA).

Implementing accelerated dynamics to unravel the effects of high-fidelity Cas9 mutants on target DNA and guide RNA hybrid stability

The clustered regularly interspersed short palindromic repeats (CRISPR) and its associated nuclease (Cas9) offers a unique and easily reprogrammable system for editing eukaryotic genomes. Cas9 is guided to the target by an RNA strand, and precise edits are created by introducing double-stranded breaks. However, nuclease activity of Cas9 is also triggered at other sites other than the target sit, which is a major limitation for various applications. Cas9 variants have been designed to improve the efficacy of the tool by introducing certain mutations. However, the on-target activity of such Cas9 variants is often seen as compromised. Hence, understanding the sub-molecular differences in the variants is essential to elucidate the factors that contribute to efficiency. The study reveals distortions in the PAM-distal regions of the nucleic hybrids as well as changes in the interactions between the Cas9 variants and RNA-DNA hybrid, contributing to the explanation for differences in on-target activity.Communicated by Ramaswamy H. Sarma.

Engineering and Design of Programmable Genome Editors

Programmable genome editors are enzymes that can be targeted to a specific location in the genome for making site-specific alterations or deletions. The engineering, design, and development of sequence-specific editors has resulted in a dramatic increase in the precision of editing for nucleotide sequences. These editors can target specific locations in a genome, in vivo. The genome editors are being deployed for the development of genetically modified organisms for agriculture and industry, and for gene therapy of inherited human genetic disorders, cancer, and immunotherapy. Experimental and computational studies of structure, binding, activity, dynamics, and folding, reviewed here, have provided valuable insights that have the potential for increasing the functional efficiency of these gene/genome editors. Biochemical and biophysical studies of the specificities of natural and engineered genome editors reveal that increased binding affinity can be detrimental because of the increase of off-target effects and that the engineering and design of genome editors with higher specificity may require modulation and control of the conformational dynamics.

Dynamics and mechanisms of CRISPR-Cas9 through the lens of computational methods

The clustered regularly interspaced short palindromic repeat (CRISPR) genome-editing revolution established the beginning of a new era in life sciences. Here, we review the role of state-of-the-art computations in the CRISPR-Cas9 revolution, from the early refinement of cryo-EM data to enhanced simulations of large-scale conformational transitions. Molecular simulations reported a mechanism for RNA binding and the formation of a catalytically competent Cas9 enzyme, in agreement with subsequent structural studies. Inspired by single-molecule experiments, molecular dynamics offered a rationale for the onset of off-target effects, while graph theory unveiled the allosteric regulation. Finally, the use of a mixed quantum-classical approach established the catalytic mechanism of DNA cleavage. Overall, molecular simulations have been instrumental in understanding the dynamics and mechanism of CRISPR-Cas9, contributing to understanding function, catalysis, allostery, and specificity.

Searching for Low Probability Opening Events in a DNA Sliding Clamp

The β subunit of E. coli DNA polymererase III is a DNA sliding clamp associated with increasing the processivity of DNA synthesis. In its free form, it is a circular homodimer structure that can accomodate double-stranded DNA in a nonspecific manner. An open state of the clamp must be accessible before loading the DNA. The opening mechanism is still a matter of debate, as is the effect of bound DNA on opening/closing kinetics. We use a combination of atomistic, coarse-grained, and enhanced sampling strategies in both explicit and implicit solvents to identify opening events in the sliding clamp. Such simulations of large nucleic acid and their complexes are becoming available and are being driven by improvements in force fields and the creation of faster computers. Different models support alternative opening mechanisms, either through an in-plane or out-of-plane opening event. We further note some of the current limitations, despite advances, in modeling these highly charged systems with implicit solvent.

Twisting and swiveling domain motions in Cas9 to recognize target DNA duplexes, make double-strand breaks, and release cleaved duplexes

The CRISPR-associated protein 9 (Cas9) has been engineered as a precise gene editing tool to make double-strand breaks. CRISPR-associated protein 9 binds the folded guide RNA (gRNA) that serves as a binding scaffold to guide it to the target DNA duplex via a RecA-like strand-displacement mechanism but without ATP binding or hydrolysis. The target search begins with the protospacer adjacent motif or PAM-interacting domain, recognizing it at the major groove of the duplex and melting its downstream duplex where an RNA-DNA heteroduplex is formed at nanomolar affinity. The rate-limiting step is the formation of an R-loop structure where the HNH domain inserts between the target heteroduplex and the displaced non-target DNA strand. Once the R-loop structure is formed, the non-target strand is rapidly cleaved by RuvC and ejected from the active site. This event is immediately followed by cleavage of the target DNA strand by the HNH domain and product release. Within CRISPR-associated protein 9, the HNH domain is inserted into the RuvC domain near the RuvC active site via two linker loops that provide allosteric communication between the two active sites. Due to the high flexibility of these loops and active sites, biophysical techniques have been instrumental in characterizing the dynamics and mechanism of the CRISPR-associated protein 9 nucleases, aiding structural studies in the visualization of the complete active sites and relevant linker structures. Here, we review biochemical, structural, and biophysical studies on the underlying mechanism with emphasis on how CRISPR-associated protein 9 selects the target DNA duplex and rejects non-target sequences.

A Review: Computational Approaches to Design sgRNA of CRISPR-Cas9

:

Clustered regularly interspaced short palindromic repeats along with CRISPR-associated protein mechanisms preserve the memory of previous experiences with DNA invaders, in particular spacers that are embedded in CRISPR arrays between coordinate repeats. There has been a fast progression in the comprehension of this immune system and its implementations; however, there are numerous points of view that anticipate explanations to make the field an energetic research zone. The efficiency of CRISPR-Cas depends upon well-considered single guide RNA; for this purpose, many bioinformatics methods and tools are created to support the design of greatly active and precise single guide RNA. Insilico single guide RNA architecture is a crucial point for effective gene editing by means of the CRISPR technique. Persistent attempts have been made to improve in-silico single guide RNA formulation having great on-target effectiveness and decreased off-target effects. This review offers a summary of the CRISPR computational tools to help different researchers pick a specific tool for their work according to pros and cons, along with new thoughts to make new computational tools to overcome all existing limitations.

Enhanced specificity mutations perturb allosteric signaling in CRISPR-Cas9

CRISPR-Cas9 (clustered regularly interspaced short palindromic repeat and associated Cas9 protein) is a molecular tool with transformative genome editing capabilities. At the molecular level, an intricate allosteric signaling is critical for DNA cleavage, but its role in the specificity enhancement of the Cas9 endonuclease is poorly understood. Here, multi-microsecond molecular dynamics is combined with solution NMR and graph theory-derived models to probe the allosteric role of key specificity-enhancing mutations. We show that mutations responsible for increasing the specificity of Cas9 alter the allosteric structure of the catalytic HNH domain, impacting the signal transmission from the DNA recognition region to the catalytic sites for cleavage. Specifically, the K855A mutation strongly disrupts the allosteric connectivity of the HNH domain, exerting the highest perturbation on the signaling transfer, while K810A and K848A result in more moderate effects on the allosteric communication. This differential perturbation of the allosteric signal correlates to the order of specificity enhancement (K855A > K848A ~ K810A) observed in biochemical studies, with the mutation achieving the highest specificity most strongly perturbing the signaling transfer. These findings suggest that alterations of the allosteric communication from DNA recognition to cleavage are critical to increasing the specificity of Cas9 and that allosteric hotspots can be targeted through mutational studies for improving the system's function.

Gaussian accelerated molecular dynamics (GaMD): principles and applications

Gaussian accelerated molecular dynamics (GaMD) is a robust computational method for simultaneous unconstrained enhanced sampling and free energy calculations of biomolecules. It works by adding a harmonic boost potential to smooth biomolecular potential energy surface and reduce energy barriers. GaMD greatly accelerates biomolecular simulations by orders of magnitude. Without the need to set predefined reaction coordinates or collective variables, GaMD provides unconstrained enhanced sampling and is advantageous for simulating complex biological processes. The GaMD boost potential exhibits a Gaussian distribution, thereby allowing for energetic reweighting via cumulant expansion to the second order (i.e., "Gaussian approximation"). This leads to accurate reconstruction of free energy landscapes of biomolecules. Hybrid schemes with other enhanced sampling methods, such as the replica exchange GaMD (rex-GaMD) and replica exchange umbrella sampling GaMD (GaREUS), have also been introduced, further improving sampling and free energy calculations. Recently, new "selective GaMD" algorithms including the ligand GaMD (LiGaMD) and peptide GaMD (Pep-GaMD) enabled microsecond simulations to capture repetitive dissociation and binding of small-molecule ligands and highly flexible peptides. The simulations then allowed highly efficient quantitative characterization of the ligand/peptide binding thermodynamics and kinetics. Taken together, GaMD and its innovative variants are applicable to simulate a wide variety of biomolecular dynamics, including protein folding, conformational changes and allostery, ligand binding, peptide binding, protein-protein/nucleic acid/carbohydrate interactions, and carbohydrate/nucleic acid interactions. In this review, we present principles of the GaMD algorithms and recent applications in biomolecular simulations and drug design.

Establishing the allosteric mechanism in CRISPR-Cas9

Allostery is a fundamental property of proteins, which regulates biochemical information transfer between spatially distant sites. Here, we report on the critical role of molecular dynamics (MD) simulations in discovering the mechanism of allosteric communication within CRISPR-Cas9, a leading genome editing machinery with enormous promises for medicine and biotechnology. MD revealed how allostery intervenes during at least three steps of the CRISPR-Cas9 function: affecting DNA recognition, mediating the cleavage and interfering with the off-target activity. An allosteric communication that activates concerted DNA cleavages was found to led through the L1/L2 loops, which connect the HNH and RuvC catalytic domains. The identification of these "allosteric transducers" inspired the development of novel variants of the Cas9 protein with improved specificity, opening a new avenue for controlling the CRISPR-Cas9 activity. Discussed studies also highlight the critical role of the recognition lobe in the conformational activation of the catalytic HNH domain. Specifically, the REC3 region was found to modulate the dynamics of HNH by sensing the formation of the RNA:DNA hybrid. The role of REC3 was revealed to be particularly relevant in the presence of DNA mismatches. Indeed, interference of REC3 with the RNA:DNA hybrid containing mismatched pairs at specific positions resulted in locking HNH in an inactive "conformational checkpoint" conformation, thereby hampering off-target cleavages. Overall, MD simulations established the fundamental mechanisms underlying the allosterism of CRISPR-Cas9, aiding engineering strategies to develop new CRISPR-Cas9 variants for improved genome editing.

Molecular Dynamics to Predict Cryo-EM: Capturing Transitions and Short-Lived Conformational States of Biomolecules

Single-particle cryogenic electron microscopy (cryo-EM) has revolutionized the field of the structural biology, providing an access to the atomic resolution structures of large biomolecular complexes in their near-native environment. Today's cryo-EM maps can frequently reach the atomic-level resolution, while often containing a range of resolutions, with conformationally variable regions obtained at 6 Å or worse. Low resolution density maps obtained for protein flexible domains, as well as the ensemble of coexisting conformational states arising from cryo-EM, poses new challenges and opportunities for Molecular Dynamics (MD) simulations. With the ability to describe the biomolecular dynamics at the atomic level, MD can extend the capabilities of cryo-EM, capturing the conformational variability and predicting biologically relevant short-lived conformational states. Here, we report about the state-of-the-art MD procedures that are currently used to refine, reconstruct and interpret cryo-EM maps. We show the capability of MD to predict short-lived conformational states, finding remarkable confirmation by cryo-EM structures subsequently solved. This has been the case of the CRISPR-Cas9 genome editing machinery, whose catalytically active structure has been predicted through both long-time scale MD and enhanced sampling techniques 2 years earlier than cryo-EM. In summary, this contribution remarks the ability of MD to complement cryo-EM, describing conformational landscapes and relating structural transitions to function, ultimately discerning relevant short-lived conformational states and providing mechanistic knowledge of biological function.

Which Factors Affect the Occurrence of Off-Target Effects Caused by the Use of CRISPR/Cas: A Systematic Review in Plants

CRISPR/Cas enables a targeted modification of DNA sequences. Despite their ease and efficient use, one limitation is the potential occurrence of associated off-target effects. This systematic review aims to answer the following research question: Which factors affect the occurrence of off-target effects caused by the use of CRISPR/Cas in plants? Literature published until March 2019 was considered for this review. Articles were screened for relevance based on pre-defined inclusion criteria. Relevant studies were subject to critical appraisal. All studies included in the systematic review were synthesized in a narrative report, but studies rated as high and medium/high validity were reported separately from studies rated as low and medium/low or unclear validity. In addition, we ran a binary logistic regression analysis to verify five factors that may affect the occurrence of off-target effects: (1) Number of mismatches (2) Position of mismatches (3) GC-content of the targeting sequence (4) Altered nuclease variants (5) Delivery methods. In total, 180 relevant articles were included in this review containing 468 studies therein. Seventy nine percentage of these studies were rated as having high or medium/high validity. Within these studies, 6,416 potential off-target sequences were assessed for the occurrence of off-target effects. Results clearly indicate that an increased number of mismatches between the on-target and potential off-target sequence steeply decreases the likelihood of off-target effects. The observed rate of off-target effects decreased from 59% when there is one mismatch between the on-target and off-target sequences toward 0% when four or more mismatches exist. In addition, mismatch/es located within the first eight nucleotides proximal to the PAM significantly decreased the occurrence of off-target effects. There is no evidence that the GC-content significantly affects off-target effects. The database regarding the impact of the nuclease variant and the delivery method is very poor as the majority of studies applied the standard nuclease SpCas9 and the CRISPR/Cas system was stably delivered in the genome. Hence, a general significant impact of these two factors on the occurrence of off-target effects cannot be proved. This identified evidence gap needs to be filled by systematic studies exploring these individual factors in sufficient numbers.

Catalytic Mechanism of Non-Target DNA Cleavage in CRISPR-Cas9 Revealed by Ab Initio Molecular Dynamics

CRISPR-Cas9 is a cutting-edge genome editing technology, which uses the endonuclease Cas9 to introduce mutations at desired sites of the genome. This revolutionary tool is promising to treat a myriad of human genetic diseases. Nevertheless, the molecular basis of DNA cleavage, which is a fundamental step for genome editing, has not been established. Here, quantum-classical molecular dynamics (MD) and free energy methods are used to disclose the two-metal-dependent mechanism of phosphodiester bond cleavage in CRISPR-Cas9. Ab initio MD reveals a conformational rearrangement of the Mg2+-bound RuvC active site, which entails the relocation of H983 to act as a general base. Then, the DNA cleavage proceeds through a concerted associative pathway fundamentally assisted by the joint dynamics of the two Mg2+ ions. This clarifies previous controversial experimental evidence, which could not fully establish the catalytic role of the conserved H983 and the metal cluster conformation. The comparison with other two-metal-dependent enzymes supports the identified mechanism and suggests a common catalytic strategy for genome editing and recombination. Overall, the non-target DNA cleavage catalysis described here resolves a fundamental open question in the CRISPR-Cas9 biology and provides valuable insights for improving the catalytic efficiency and the metal-dependent function of the Cas9 enzyme, which are at the basis of the development of genome editing tools.