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

Gene Editing in Dimorphic Fungi Using CRISPR/Cas9

Dimorphic fungi in the genera Blastomyces, Histoplasma, Coccidioides, and Paracoccidioides are important human pathogens that affect human health in many countries throughout the world. Understanding the biology of these fungi is important for the development of effective treatments and vaccines. Gene editing is a critically important tool for research into these organisms. In recent years, gene targeting approaches employing RNA-guided DNA nucleases, such as clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated nuclease 9 (Cas9), have exploded in popularity. Here, we provide a detailed description of the steps involved in applying CRISPR/Cas9 technology to dimorphic fungi, with Blastomyces dermatitidis in particular as our model fungal pathogen. We discuss the design and construction of single guide RNA and Cas9-expressing targeting vectors (including multiplexed vectors) as well as introduction of these plasmids into Blastomyces using Agrobacterium-mediated transformation. Finally, we cover the outcomes that may be expected in terms of gene-editing efficiency and types of gene alterations produced. © 2020 Wiley Periodicals LLC. Basic Protocol 1: Construction of CRISPR/Cas9 targeting vectors Support Protocol 1: Choosing protospacers in the target gene Basic Protocol 2: Agrobacterium-mediated transformation of Blastomyces Support Protocol 2: Preparation of electrocompetent Agrobacterium Support Protocol 3: Preparation and recovery of Blastomyces frozen stocks.

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.

Cas9 Cuts and Consequences; Detecting, Predicting, and Mitigating CRISPR/Cas9 On- and Off-Target Damage: Techniques for Detecting, Predicting, and Mitigating the On- and off-target Effects of Cas9 Editing

Large deletions and genomic re-arrangements are increasingly recognized as common products of double-strand break repair at Clustered Regularly Interspaced, Short Palindromic Repeats - CRISPR associated protein 9 (CRISPR/Cas9) on-target sites. Together with well-known off-target editing products from Cas9 target misrecognition, these are important limitations, that need to be addressed. Rigorous assessment of Cas9-editing is necessary to ensure validity of observed phenotypes in Cas9-edited cell-lines and model organisms. Here the mechanisms of Cas9 specificity, and strategies to assess and mitigate unwanted effects of Cas9 editing are reviewed; covering guide-RNA design, RNA modifications, Cas9 modifications, control of Cas9 activity; computational prediction for off-targets, and experimental methods for detecting Cas9 cleavage. Although recognition of the prevalence of on- and off-target effects of Cas9 editing has increased in recent years, broader uptake across the gene editing community will be important in determining the specificity of Cas9 across diverse applications and organisms.

Defining the mutation signatures of DNA polymerase θ in cancer genomes

DNA polymerase theta (POLQ)-mediated end joining (TMEJ) is a distinct pathway for mediating DNA double-strand break (DSB) repair. TMEJ is required for the viability of BRCA-mutated cancer cells. It is crucial to identify tumors that rely on POLQ activity for DSB repair, because such tumors are defective in other DSB repair pathways and have predicted sensitivity to POLQ inhibition and to cancer therapies that produce DSBs. We define here the POLQ-associated mutation signatures in human cancers, characterized by short insertions and deletions in a specific range of microhomologies. By analyzing 82 COSMIC (Catalogue of Somatic Mutations in Cancer) signatures, we found that BRCA-mutated cancers with a higher level of POLQ expression have a greatly enhanced representation of the small insertion and deletion signature 6, as well as single base substitution signature 3. Using human cancer cells with disruptions of POLQ, we further show that TMEJ dominates end joining of two separated DSBs (distal EJ). Templated insertions with microhomology are enriched in POLQ-dependent distal EJ. The use of this signature analysis will aid in identifying tumors relying on POLQ activity.

Review of Artificial Intelligence Applications and Algorithms for Brain Organoid Research

The human brain organoid is a miniature three-dimensional tissue culture that can simulate the structure and function of the brain in an in vitro culture environment. Although we consider that human brain organoids could be used to understand brain development and diseases, experimental models of human brain organoids are so highly variable that we apply artificial intelligence (AI) techniques to investigate the development mechanism of the human brain. Therefore, this study briefly reviewed commonly used AI applications for human brain organoid-magnetic resonance imaging, electroencephalography, and gene editing techniques, as well as related AI algorithms. Finally, we discussed the limitations, challenges, and future study direction of AI-based technology for human brain organoids.

Determinants of Base Editing Outcomes from Target Library Analysis and Machine Learning

Although base editors are widely used to install targeted point mutations, the factors that determine base editing outcomes are not well understood. We characterized sequence-activity relationships of 11 cytosine and adenine base editors (CBEs and ABEs) on 38,538 genomically integrated targets in mammalian cells and used the resulting outcomes to train BE-Hive, a machine learning model that accurately predicts base editing genotypic outcomes (R ≈ 0.9) and efficiency (R ≈ 0.7). We corrected 3,388 disease-associated SNVs with ≥90% precision, including 675 alleles with bystander nucleotides that BE-Hive correctly predicted would not be edited. We discovered determinants of previously unpredictable C-to-G, or C-to-A editing and used these discoveries to correct coding sequences of 174 pathogenic transversion SNVs with ≥90% precision. Finally, we used insights from BE-Hive to engineer novel CBE variants that modulate editing outcomes. These discoveries illuminate base editing, enable editing at previously intractable targets, and provide new base editors with improved editing capabilities.

SpCas9- and LbCas12a-Mediated DNA Editing Produce Different Gene Knockout Outcomes in Zebrafish Embryos

CRISPR/Cas (Clustered Regularly Interspaced Short Palindromic Repeats/CRISPR associated protein) genome editing is a powerful technology widely used in current genetic research. In the most simple and straightforward way it can be applied for a gene knockout resulting from repair errors, induced by dsDNA cleavage by Cas nuclease. For decades, zebrafish (Danio rerio) has been known as a convenient model object of developmental biology. Both commonly used nucleases SpCas9 (Streptococcus pyogenes Cas9) and LbCas12a (Lachnospiraceae bacterium Cas12a) are extensively used in this model. Among them, LbCas12a is featured with higher specificity and efficiency of homology-directed editing in human cells and mouse. But the editing outcomes for these two nucleases in zebrafish are still not compared quantitatively. Therefore, to reveal possible advantages of one nuclease in comparison to the other in the context of gene knockout generation, we compare here the outcomes of repair of the DNA breaks introduced by these two commonly used nucleases in zebrafish embryos. To address this question, we microinjected the ribonucleoprotein complexes of the both nucleases with the corresponding guide RNAs in zebrafish zygotes and sequenced the target gene regions after three days of development. We found that LbCas12a editing resulted in longer deletions and more rare inserts, in comparison to those generated by SpCas9, while the editing efficiencies (percentage of mutated copies of the target gene to all gene copies in the embryo) of both nucleases were the same. On the other hand, overlapping of protospacers resulted in similarities in repair outcome, although they were cut by two different nucleases. Thus, our results indicate that the repair outcome depends both on the nuclease mode of action and on protospacer sequence.

Precise, predictable multi-nucleotide deletions in rice and wheat using APOBEC-Cas9

Short insertions and deletions can be produced in plant genomes using CRISPR-Cas editors, but reliable production of larger deletions in specific target sites has proven difficult to achieve. We report the development of a series of APOBEC-Cas9 fusion-induced deletion systems (AFIDs) that combine Cas9 with human APOBEC3A (A3A), uracil DNA-glucosidase and apurinic or apyrimidinic site lyase. In rice and wheat, AFID-3 generated deletions from 5'-deaminated C bases to the Cas9-cleavage site. Approximately one-third of deletions produced using AFID-3 in rice and wheat protoplasts (30.2%) and regenerated plants (34.8%) were predictable. We show that eAFID-3, in which the A3A in AFID-3 is replaced with truncated APOBEC3B (A3Bctd), produced more uniform deletions from the preferred TC motif to the double-strand break. AFIDs could be applied to study regulatory regions and protein domains to improve crop plants.

Highly Parallel Profiling of Cas9 Variant Specificity

Determining the off-target cleavage profile of programmable nucleases is an important consideration for any genome editing experiment, and a number of Cas9 variants have been reported that improve specificity. We describe here tagmentation-based tag integration site sequencing (TTISS), an efficient, scalable method for analyzing double-strand breaks (DSBs) that we apply in parallel to eight Cas9 variants across 59 targets. Additionally, we generated thousands of other Cas9 variants and screened for variants with enhanced specificity and activity, identifying LZ3 Cas9, a high specificity variant with a unique +1 insertion profile. This comprehensive comparison reveals a general trade-off between Cas9 activity and specificity and provides information about the frequency of generation of +1 insertions, which has implications for correcting frameshift mutations.

Protein-Mutation-Induced Conformational Changes of the DNA and Nuclease Domain in CRISPR/Cas9 Systems by Molecular Dynamics Simulations

Class 2 CRISPR (clustered regularly interspaced short palindromic repeats) systems offer a unique protocol for genome editing in eukaryotic cells. The nuclease activity of Cas9 has been harnessed to perform precise genome editing by creating double-strand breaks. However, the nuclease activity of Cas9 can be triggered when there is imperfect complementarity between the RNA guide sequence and an off-target genomic site, which is a major limitation of the CRISPR technique for practical applications. Hence, understanding the binding mechanisms in CRISPR/Cas9 for predicting ways to increase cleavage specificity is a timely research target. One way to understand and tune the binding strength is to study wild-type and mutant Cas9, in complex with a guide RNA and a target DNA. We have performed classical all-atom MD simulations over a cumulative time scale of 13.5 μs of CRISPR/Cas9 ternary complexes with the wild-type Cas9 from Streptococcus pyogenes and three of its mutants: K855A, H982A, and the combination K855A+H982A, selected from the outcome of experimental work. Our results reveal significant structural impact of the mutations, with implications for specificity. We find that the "unwound" part of the nontarget DNA strand exhibits enhanced flexibility in complexes with Cas9 mutants and tries to move away from the HNH/RuvC interface, where it is otherwise stabilized by electrostatic couplings in the wild-type complex. Our findings refine an electrostatic model by which cleavage specificity can be optimized through protein mutations.

Enhanced CRISPR-Cas9 correction of Duchenne muscular dystrophy in mice by a self-complementary AAV delivery system

Duchenne muscular dystrophy (DMD) is a lethal neuromuscular disease caused by mutations in the dystrophin gene (DMD). Previously, we applied CRISPR-Cas9-mediated "single-cut" genome editing to correct diverse genetic mutations in animal models of DMD. However, high doses of adeno-associated virus (AAV) are required for efficient in vivo genome editing, posing challenges for clinical application. In this study, we packaged Cas9 nuclease in single-stranded AAV (ssAAV) and CRISPR single guide RNAs in self-complementary AAV (scAAV) and delivered this dual AAV system into a mouse model of DMD. The dose of scAAV required for efficient genome editing were at least 20-fold lower than with ssAAV. Mice receiving systemic treatment showed restoration of dystrophin expression and improved muscle contractility. These findings show that the efficiency of CRISPR-Cas9-mediated genome editing can be substantially improved by using the scAAV system. This represents an important advancement toward therapeutic translation of genome editing for DMD.

Recent advances in the CRISPR genome editing tool set

Genome editing took a dramatic turn with the development of the clustered regularly interspaced short palindromic repeats (CRISPR)-CRISPR-associated proteins (Cas) system. The CRISPR-Cas system is functionally divided into classes 1 and 2 according to the composition of the effector genes. Class 2 consists of a single effector nuclease, and routine practice of genome editing has been achieved by the development of the Class 2 CRISPR-Cas system, which includes the type II, V, and VI CRISPR-Cas systems. Types II and V can be used for DNA editing, while type VI is employed for RNA editing. CRISPR techniques induce both qualitative and quantitative alterations in gene expression via the double-stranded breakage (DSB) repair pathway, base editing, transposase-dependent DNA integration, and gene regulation using the CRISPR-dCas or type VI CRISPR system. Despite significant technical improvements, technical challenges should be further addressed, including insufficient indel and HDR efficiency, off-target activity, the large size of Cas, PAM restrictions, and immune responses. If sophisticatedly refined, CRISPR technology will harness the process of DNA rewriting, which has potential applications in therapeutics, diagnostics, and biotechnology.

Structural and functional insights into the bona fide catalytic state of Streptococcus pyogenes Cas9 HNH nuclease domain

The CRISPR-associated endonuclease Cas9 from Streptococcus pyogenes (SpyCas9), along with a programmable single-guide RNA (sgRNA), has been exploited as a significant genome-editing tool. Despite the recent advances in determining the SpyCas9 structures and DNA cleavage mechanism, the cleavage-competent conformation of the catalytic HNH nuclease domain of SpyCas9 remains largely elusive and debatable. By integrating computational and experimental approaches, we unveiled and validated the activated Cas9-sgRNA-DNA ternary complex in which the HNH domain is neatly poised for cleaving the target DNA strand. In this catalysis model, the HNH employs the catalytic triad of D839-H840-N863 for cleavage catalysis, rather than previously implicated D839-H840-D861, D837-D839-H840, or D839-H840-D861-N863. Our study contributes critical information to defining the catalytic conformation of the HNH domain and advances the knowledge about the conformational activation underlying Cas9-mediated DNA cleavage.

Deciphering Off-Target Effects in CRISPR-Cas9 through Accelerated Molecular Dynamics

CRISPR-Cas9 is the state-of-the-art technology for editing and manipulating nucleic acids. However, the occurrence of off-target mutations can limit its applicability. Here, all-atom enhanced molecular dynamics (MD) simulations-using Gaussian accelerated MD (GaMD)-are used to decipher the mechanism of off-target binding at the molecular level. GaMD reveals that base pair mismatches in the target DNA at distal sites with respect to the protospacer adjacent motif (PAM) can induce an extended opening of the RNA:DNA heteroduplex, which leads to newly formed interactions between the unwound DNA and the L2 loop of the catalytic HNH domain. These conserved interactions constitute a "lock" effectively decreasing the conformational freedom of the HNH domain and hampering its activation for cleavage. Remarkably, depending on their positions at PAM distal sites, DNA mismatches responsible for off-target cleavages are unable to "lock" the HNH domain, thereby leading to the unselective cleavage of DNA sequences. In consistency with the available experimental data, the ability to "lock" the catalytic HNH domain in an inactive "conformational checkpoint" is shown to be a key determinant in the onset of off-target effects. This mechanistic rationale contributes in clarifying a long lasting open issue in the CRISPR-Cas9 function and poses the foundation for designing novel and more specific Cas9 variants, which could be obtained by magnifying the "locking" interactions between HNH and the target DNA in the presence of any incorrect off-target sequence, thus preventing undesired cleavages.

Multiplexed Cas9 targeting reveals genomic location effects and gRNA-based staggered breaks influencing mutation efficiency

Understanding the impact of guide RNA (gRNA) and genomic locus on CRISPR-Cas9 activity is crucial to design effective gene editing assays. However, it is challenging to profile Cas9 activity in the endogenous cellular environment. Here we leverage our TRIP technology to integrate ~ 1k barcoded reporter genes in the genomes of mouse embryonic stem cells. We target the integrated reporters (IRs) using RNA-guided Cas9 and characterize induced mutations by sequencing. We report that gRNA-sequence and IR locus explain most variation in mutation efficiency. Predominant insertions of a gRNA-specific nucleotide are consistent with template-dependent repair of staggered DNA ends with 1-bp 5' overhangs. We confirm that such staggered ends are induced by Cas9 in mouse pre-B cells. To explain observed insertions, we propose a model generating primarily blunt and occasionally staggered DNA ends. Mutation patterns indicate that gRNA-sequence controls the fraction of staggered ends, which could be used to optimize Cas9-based insertion efficiency.

Single and multiple gene knockouts by CRISPR-Cas9 in maize

The analysis of 93 mutant alleles in 18 genes demonstrated that CRISPR-Cas9 is a robust tool for targeted mutagenesis in maize, permitting efficient generation of single and multiple knockouts. CRISPR-Cas9 technology is a simple and efficient tool for targeted mutagenesis of the genome. It has been implemented in many plant species, including crops such as maize. Here we report single- and multiple-gene mutagenesis via stably transformed maize plants. Two different CRISPR-Cas9 vectors were used allowing the expression of multiple guide RNAs and different strategies to knockout either independent or paralogous genes. A total of 12 plasmids, representing 28 different single guide RNAs (sgRNAs), were generated to target 20 genes. For 18 of these genes, at least one mutant allele was obtained, while two genes were recalcitrant to sequence editing. 19% (16/83) of mutant plants showed biallelic mutations. Small insertions or deletions of less than ten nucleotides were most frequently observed, regardless of whether the gene was targeted by one or more sgRNAs. Deletions of defined regions located between the target sites of two guide RNAs were also reported although the exact deletion size was variable. Double and triple mutants were created in a single step, which is especially valuable for functional analysis of genes with strong genetic linkage. Off-target effects were theoretically limited due to rigorous sgRNA design and random experimental checks at three potential off-target sites did not reveal any editing. Sanger chromatograms allowed to unambiguously class the primary transformants; the majority (85%) were fully edited plants transmitting systematically all detected mutations to the next generation, generally following Mendelian segregation.

Structure and Dynamics of the CRISPR-Cas9 Catalytic Complex

CRISPR-Cas9 is a bacterial immune system with exciting applications for genome editing. In spite of extensive experimental characterization, the active site chemistry of the RuvC domain-which performs DNA cleavages-has remained elusive. Its knowledge is key for structure-based engineering aimed at improving DNA cleavages. Here, we deliver an in-depth characterization by using quantum-classical (QM/MM) molecular dynamics (MD) simulations and a Gaussian accelerated MD method, coupled with bioinformatics analysis. We disclose a two-metal aided architecture in the RuvC active site, which is poised to operate DNA cleavages, in analogy with other DNA/RNA processing enzymes. The conformational dynamics of the RuvC domain further reveals that an "arginine finger" stably contacts the scissile phosphate, with the function of stabilizing the active complex. Remarkably, the formation of a catalytically competent state of the RuvC domain is only observed upon the conformational activation of the other nuclease domain of CRISPR-Cas9-i.e., the HNH domain-such allowing concerted cleavages of double stranded DNA. This structure is in agreement with the available experimental data and remarkably differs from previous models based on classical mechanics, demonstrating also that only quantum mechanical simulations can accurately describe the metal-aided active site in CRISPR-Cas9. This fully catalytic structure-in which both the HNH and RuvC domains are prone to perform DNA cleavages-constitutes a stepping-stone for understanding DNA cleavage and specificity. It calls for novel experimental verifications and offers the structural foundations for engineering efforts aimed at improving the genome editing capability of CRISPR-Cas9.

Strategies for Efficient Genome Editing Using CRISPR-Cas9

The targetable DNA endonuclease CRISPR-Cas9 has transformed analysis of biological processes by enabling robust genome editing in model and nonmodel organisms. Although rules directing Cas9 to its target DNA via a guide RNA are straightforward, wide variation occurs in editing efficiency and repair outcomes for both imprecise error-prone repair and precise templated repair. We found that imprecise and precise DNA repair from double-strand breaks (DSBs) is asymmetric, favoring repair in one direction. Using this knowledge, we designed RNA guides and repair templates that increased the frequency of imprecise insertions and deletions and greatly enhanced precise insertion of point mutations in Caenorhabditis elegans We also devised strategies to insert long (10 kb) exogenous sequences and incorporate multiple nucleotide substitutions at a considerable distance from DSBs. We expanded the repertoire of co-conversion markers appropriate for diverse nematode species. These selectable markers enable rapid identification of Cas9-edited animals also likely to carry edits in desired targets. Lastly, we explored the timing, location, frequency, sex dependence, and categories of DSB repair events by developing loci with allele-specific Cas9 targets that can be contributed during mating from either male or hermaphrodite germ cells. We found a striking difference in editing efficiency between maternally and paternally contributed genomes. Furthermore, imprecise repair and precise repair from exogenous repair templates occur with high frequency before and after fertilization. Our strategies enhance Cas9-targeting efficiency, lend insight into the timing and mechanisms of DSB repair, and establish guidelines for achieving predictable precise and imprecise repair outcomes with high frequency.

Assessing the Performance of the Nonbonded Mg2+ Models in a Two-Metal-Dependent Ribonuclease

Magnesium ions (Mg²⁺), abundant in living cells, are essential for biomolecular structure, dynamics, and function. The biological importance of Mg²⁺ has motivated continuous development and improvement of various Mg²⁺ models for molecular dynamics (MD) simulations during the last decades. There are four types of nonbonded Mg²⁺ models: the point charge models based on a 12-6 or 12-6-4 type Lennard-Jones (LJ) potential, and the multisite models based on a 12-6 or 12-6-4 LJ potential. Here, we systematically assessed the performance of these four types of nonbonded Mg²⁺ models (21 models in total) in terms of maintaining a challenging intermediate state configuration captured in the structure of a prototypical two-metal-ion RNase H complex with an RNA/DNA hybrid. Our data demonstrate that the 12-6-4 multisite models, which account for charge-induced dipole interactions, perform the best in reproducing all the unique coordination modes in this intermediate state and maintaining the correct carboxylate denticity. Our benchmark work provides a useful guideline for MD simulations and structural refinement of Mg²⁺-containing biomolecular systems.

Predicting the mutations generated by repair of Cas9-induced double-strand breaks

The DNA mutation produced by cellular repair of a CRISPR-Cas9-generated double-strand break determines its phenotypic effect. It is known that the mutational outcomes are not random, but depend on DNA sequence at the targeted location. Here we systematically study the influence of flanking DNA sequence on repair outcome by measuring the edits generated by >40,000 guide RNAs (gRNAs) in synthetic constructs. We performed the experiments in a range of genetic backgrounds and using alternative CRISPR-Cas9 reagents. In total, we gathered data for >109 mutational outcomes. The majority of reproducible mutations are insertions of a single base, short deletions or longer microhomology-mediated deletions. Each gRNA has an individual cell-line-dependent bias toward particular outcomes. We uncover sequence determinants of the mutations produced and use these to derive a predictor of Cas9 editing outcomes. Improved understanding of sequence repair will allow better design of gene editing experiments.

Harnessing accurate non-homologous end joining for efficient precise deletion in CRISPR/Cas9-mediated genome editing

Background

Many applications of CRISPR/Cas9-mediated genome editing require Cas9-induced non-homologous end joining (NHEJ), which was thought to be error prone. However, with directly ligatable ends, Cas9-induced DNA double strand breaks may be repaired preferentially by accurate NHEJ.

Results

In the repair of two adjacent double strand breaks induced by paired Cas9-gRNAs at 71 genome sites, accurate NHEJ accounts for about 50% of NHEJ events. This paired Cas9-gRNA approach underestimates the level of accurate NHEJ due to frequent + 1 templated insertions, which can be avoided by the predefined Watson/Crick orientation of protospacer adjacent motifs (PAMs). The paired Cas9-gRNA strategy also provides a flexible, reporter-less approach for analyzing both accurate and mutagenic NHEJ in cells and in vivo, and it has been validated in cells deficient for XRCC4 and in mouse liver. Due to high frequencies of precise deletions of defined “3n”-, “3n + 1”-, or “3n + 2”-bp length, accurate NHEJ is used to improve the efficiency and homogeneity of gene knockouts and targeted in-frame deletions. Compared to “3n + 1”-bp, “3n + 2”-bp can overcome + 1 templated insertions to increase the frequency of out-of-frame mutations. By applying paired Cas9-gRNAs to edit MDC1 and key 53BP1 domains, we are able to generate predicted, precise deletions for functional analysis. Lastly, a Plk3 inhibitor promotes NHEJ with bias towards accurate NHEJ, providing a chemical approach to improve genome editing requiring precise deletions.

Conclusions

NHEJ is inherently accurate in repair of Cas9-induced DNA double strand breaks and can be harnessed to improve CRISPR/Cas9 genome editing requiring precise deletion of a defined length.

Electronic supplementary material

The online version of this article (10.1186/s13059-018-1518-x) contains supplementary material, which is available to authorized users.

Molecular dynamics provides insight into how N251A and N251Y mutations in the active site of Bacillus licheniformis RN-01 levansucrase disrupt production of long-chain levan

Produced by levansucrase, levan and levan oligosaccharides (GF_n) have potential applications in food and pharmaceutical industries such as prebiotics, anti-tumor and anti-inflammatory agents. Previous study reported that Bacillus licheniformis RN-01 levansucrase could produce levan oligosaccharides and long-chain levan. However, its N251A and N251Y mutants could effectively produce short-chain oligosaccharides upto GF_3, but they could not produce long-chain levan. We hypothesized that these mutations probably reduced GF₃ binding affinity in levansucrase active site that contains fructosyl-Asp93 intermediate and caused GF₃ to be in an unfavorable orientation for transfructosylation; therefore, levansucrase could not effectively extend GF₃ by one fructosyl residue to produce GF₄ and subsequently long-chain levan. However, these mutations probably did not significantly reduce binding affinity or drastically change orientation of GF₂; therefore, levansucrase could still extend GF₂ to produce GF₃. Using this hypothesis, we employed molecular dynamics to investigate effects of these mutations on GF₂/GF₃ binding in levansucrase active site. Our results reasonably support this hypothesis as N251A and N251Y mutations did not significantly reduce GF₂ binding affinity, as calculated by MM-GBSA technique and hydrogen bond occupations, or drastically change orientation of GF₂ in levansucrase active site, as measured by distance between atoms necessary for transfructosylation. However, these mutations drastically decreased GF₃ binding affinity and caused GF₃ to be in an unfavorable orientation for transfructosylation. Furthermore, the free energy decomposition and hydrogen bond occupation results suggest the importance of Arg255 in GF₂/GF₃ binding in levansucrase active site. This study provides important and novel insight into the effects of N251A and N251Y mutations on GF₂/GF₃ binding in levansucrase active site and how they may disrupt production of long-chain levan. This knowledge could be beneficial in designing levansucrase to efficiently produce levan oligosaccharides with desired length.

Decoding non-random mutational signatures at Cas9 targeted sites

The mutation patterns at Cas9 targeted sites contain unique information regarding the nuclease activity and repair mechanisms in mammalian cells. However, analytical framework for extracting such information are lacking. Here, we present a novel computational platform called Rational InDel Meta-Analysis (RIMA) that enables an in-depth comprehensive analysis of Cas9-induced genetic alterations, especially InDels mutations. RIMA can be used to quantitate the contribution of classical microhomology-mediated end joining (c-MMEJ) pathway in the formation of mutations at Cas9 target sites. We used RIMA to compare mutational signatures at 15 independent Cas9 target sites in human A549 wildtype and A549-POLQ knockout cells to elucidate the role of DNA polymerase θ in c-MMEJ. Moreover, the single nucleotide insertions at the Cas9 target sites represent duplications of preceding nucleotides, suggesting that the flexibility of the Cas9 nuclease domains results in both blunt- and staggered-end cuts. Thymine at the fourth nucleotide before protospacer adjacent motif (PAM) results in a two-fold higher occurrence of single nucleotide InDels compared to guanine at the same position. This study provides a novel approach for the characterization of the Cas9 nucleases with improved accuracy in predicting genome editing outcomes and a potential strategy for homology-independent targeted genomic integration.

Developmental barcoding of whole mouse via homing CRISPR

In vivo barcoding using nuclease-induced mutations is a powerful approach for recording biological information, including developmental lineages; however, its application in mammalian systems has been limited. We present in vivo barcoding in the mouse with multiple homing guide RNAs that each generate hundreds of mutant alleles and combine to produce an exponential diversity of barcodes. Activation upon conception and continued mutagenesis through gestation resulted in developmentally barcoded mice wherein information is recorded in lineage-specific mutations. We used these recordings for reliable post hoc reconstruction of the earliest lineages and investigation of axis development in the brain. Our results provide an enabling and versatile platform for in vivo barcoding and lineage tracing in a mammalian model system.

Insight into the substrate specificity change caused by the Y227H mutation of α-glucosidase III from the European honeybee (Apis mellifera) through molecular dynamics simulations

Honey from the European honeybee, Apis mellifera, is produced by α-glucosidases (HBGases) and is widely used in food, pharmaceutical, and cosmetic industries. Categorized by their substrate specificities, HBGases have three isoforms: HBGase I, II and III. Previous experimental investigations showed that wild-type HBGase III from Apis mellifera (WT) preferred sucrose to maltose as a substrate, while the Y227H mutant (MT) preferred maltose to sucrose. This mutant can potentially be used for malt hydrolysis because it can efficiently hydrolyze maltose. In this work, to elucidate important factors contributing to substrate specificity of this enzyme and gain insight into how the Y227H mutation causes substrate specificity change, WT and MT homology models were constructed, and sucrose/maltose was docked into active sites of the WT and MT. AMBER14 was employed to perform three independent molecular dynamics runs for these four complexes. Based on the relative binding free energies calculated by the MM-GBSA method, sucrose is better than maltose for WT binding, while maltose is better than sucrose for MT binding. These rankings support the experimentally observed substrate specificity that WT preferred sucrose to maltose as a substrate, while MT preferred maltose to sucrose, suggesting the importance of binding affinity for substrate specificity. We also found that the Y227H mutation caused changes in the proximities between the atoms necessary for sucrose/maltose hydrolysis that may affect enzyme efficiency in the hydrolysis of sucrose/maltose. Moreover, the per-residue binding free energy decomposition results show that Y227/H227 may be a key residue for preference binding of sucrose/maltose in the WT/MT active site. Our study provides important and novel insight into the binding of sucrose/maltose in the active site of Apis mellifera HBGase III and into how the Y227H mutation leads to the substrate specificity change at the molecular level. This knowledge could be beneficial in the design of this enzyme for increased production of desired products.

Identification of a unique Ca2+-binding site in rat acid-sensing ion channel 3

Acid-sensing ion channels (ASICs) evolved to sense changes in extracellular acidity with the divalent cation calcium (Ca²⁺) as an allosteric modulator and channel blocker. The channel-blocking activity is most apparent in ASIC3, as removing Ca²⁺ results in channel opening, with the site's location remaining unresolved. Here we show that a ring of rat ASIC3 (rASIC3) glutamates (Glu435), located above the channel gate, modulates proton sensitivity and contributes to the formation of the elusive Ca²⁺ block site. Mutation of this residue to glycine, the equivalent residue in chicken ASIC1, diminished the rASIC3 Ca²⁺ block effect. Atomistic molecular dynamic simulations corroborate the involvement of this acidic residue in forming a high-affinity Ca²⁺ site atop the channel pore. Furthermore, the reported observations provide clarity for past controversies regarding ASIC channel gating. Our findings enhance understanding of ASIC gating mechanisms and provide structural and energetic insights into this unique calcium-binding site.

RNA Structural Dynamics As Captured by Molecular Simulations: A Comprehensive Overview

With both catalytic and genetic functions, ribonucleic acid (RNA) is perhaps the most pluripotent chemical species in molecular biology, and its functions are intimately linked to its structure and dynamics. Computer simulations, and in particular atomistic molecular dynamics (MD), allow structural dynamics of biomolecular systems to be investigated with unprecedented temporal and spatial resolution. We here provide a comprehensive overview of the fast-developing field of MD simulations of RNA molecules. We begin with an in-depth, evaluatory coverage of the most fundamental methodological challenges that set the basis for the future development of the field, in particular, the current developments and inherent physical limitations of the atomistic force fields and the recent advances in a broad spectrum of enhanced sampling methods. We also survey the closely related field of coarse-grained modeling of RNA systems. After dealing with the methodological aspects, we provide an exhaustive overview of the available RNA simulation literature, ranging from studies of the smallest RNA oligonucleotides to investigations of the entire ribosome. Our review encompasses tetranucleotides, tetraloops, a number of small RNA motifs, A-helix RNA, kissing-loop complexes, the TAR RNA element, the decoding center and other important regions of the ribosome, as well as assorted others systems. Extended sections are devoted to RNA-ion interactions, ribozymes, riboswitches, and protein/RNA complexes. Our overview is written for as broad of an audience as possible, aiming to provide a much-needed interdisciplinary bridge between computation and experiment, together with a perspective on the future of the field.

Real-time observation of flexible domain movements in CRISPR-Cas9

The CRISPR‐associated protein Cas9 is widely used for genome editing because it cleaves target DNA through the assistance of a single‐guide RNA (sgRNA). Structural studies have revealed the multi‐domain architecture of Cas9 and suggested sequential domain movements of Cas9 upon binding to the sgRNA and the target DNA. These studies also hinted at the flexibility between domains; however, it remains unclear whether these flexible movements occur in solution. Here, we directly observed dynamic fluctuations of multiple Cas9 domains, using single‐molecule FRET. We found that the flexible domain movements allow Cas9 to adopt transient conformations beyond those captured in the crystal structures. Importantly, the HNH nuclease domain only accessed the DNA cleavage position during such flexible movements, suggesting the importance of this flexibility in the DNA cleavage process. Our FRET data also revealed the conformational flexibility of apo‐Cas9, which may play a role in the assembly with the sgRNA. Collectively, our results highlight the potential role of domain fluctuations in driving Cas9‐catalyzed DNA cleavage.

Mechanisms of improved specificity of engineered Cas9s revealed by single-molecule FRET analysis

Cas9 (from Streptococcus pyogenes) in complex with a guide RNA targets complementary DNA for cleavage. Here, we developed a single-molecule FRET analysis to study the mechanisms of specificity enhancement of two engineered Cas9s (eCas9 and Cas9-HF1). A DNA-unwinding assay showed that mismatches affect cleavage reactions through rebalancing the unwinding-rewinding equilibrium. Increasing PAM-distal mismatches facilitates rewinding, and the associated cleavage impairment shows that cleavage proceeds from the unwound state. Engineered Cas9s depopulate the unwound state more readily with mismatches. The intrinsic cleavage rate is much lower for engineered Cas9s, preventing cleavage from transiently unwound off-targets. Engineered Cas9s require approximately one additional base pair match for stable binding, freeing them from sites that would otherwise sequester them. Therefore, engineered Cas9s achieve their improved specificity by inhibiting stable DNA binding to partially matching sequences, making DNA unwinding more sensitive to mismatches and slowing down the intrinsic cleavage reaction.

Harnessing "A Billion Years of Experimentation": The Ongoing Exploration and Exploitation of CRISPR-Cas Immune Systems

The famed physicist-turned-biologist, Max Delbrück, once remarked that, for physicists, "the field of bacterial viruses is a fine playground for serious children who ask ambitious questions." Early discoveries in that playground helped establish molecular genetics, and half a century later, biologists delving into the same field have ushered in the era of precision genome engineering. The focus has of course shifted-from bacterial viruses and their mechanisms of infection to the bacterial hosts and their mechanisms of immunity-but it is the very same evolutionary arms race that continues to awe and inspire researchers worldwide. In this review, we explore the remarkable diversity of CRISPR-Cas adaptive immune systems, describe the molecular components that mediate nucleic acid targeting, and outline the use of these RNA-guided machines for biotechnology applications. CRISPR-Cas research has yielded far more than just Cas9-based genome-editing tools, and the wide-reaching, innovative impacts of this fascinating biological playground are sure to be felt for years to come.

CRISPR/Cas9 cleavages in budding yeast reveal templated insertions and strand-specific insertion/deletion profiles

Harnessing CRISPR-Cas9 technology provides an unprecedented ability to modify genomic loci via DNA double-strand break (DSB) induction and repair. We analyzed nonhomologous end-joining (NHEJ) repair induced by Cas9 in budding yeast and found that the orientation of binding of Cas9 and its guide RNA (gRNA) profoundly influences the pattern of insertion/deletions (indels) at the site of cleavage. A common indel created by Cas9 is a 1-bp (+1) insertion that appears to result from Cas9 creating a 1-nt 5' overhang that is filled in by a DNA polymerase and ligated. The origin of +1 insertions was investigated by using two gRNAs with PAM sequences located on opposite DNA strands but designed to cleave the same sequence. These templated +1 insertions are dependent on the X-family DNA polymerase, Pol4. Deleting Pol4 also eliminated +2 and +3 insertions, which are biased toward homonucleotide insertions. Using inverted PAM sequences, we also found significant differences in overall NHEJ efficiency and repair profiles, suggesting that the binding of the Cas9:gRNA complex influences subsequent NHEJ processing. As with events induced by the site-specific HO endonuclease, CRISPR-Cas9-mediated NHEJ repair depends on the Ku heterodimer and DNA ligase 4. Cas9 events are highly dependent on the Mre11-Rad50-Xrs2 complex, independent of Mre11's nuclease activity. Inspection of the outcomes of a large number of Cas9 cleavage events in mammalian cells reveals a similar templated origin of +1 insertions in human cells, but also a significant frequency of similarly templated +2 insertions.

Bidirectional Degradation of DNA Cleavage Products Catalyzed by CRISPR/Cas9

Since the initial characterization of Streptococcus pyogenes CRISPR/Cas9 as a powerful gene-editing tool, it has been widely accepted that Cas9 generates blunt-ended DNA products by concerted cleavage of the target (tDNA) and non-target (ntDNA) strands three nucleotides away from the protospacer adjacent motif (PAM) by HNH and RuvC nuclease active sites, respectively. Following initial DNA cleavage, RuvC catalyzes 3'→5' degradation of the ntDNA resulting in DNA products of various lengths. Here, we found that Cas9 selects multiple sites for initial ntDNA cleavage and preferentially generates staggered-ended DNA products containing single-nucleotide 5'-overhangs. We also quantitatively evaluated 3'→5' post-cleavage trimming (PCT) activity of RuvC to find that ntDNA degradation continues up to the -10 position on the PAM distal DNA product and is kinetically significant when compared to extremely slow DNA product release. We also discovered a previously unidentified 5'→3' PCT activity of RuvC which can shorten the PAM proximal ntDNA product by precisely one nucleotide with a comparable rate as the 3'→5' PCT activity. Taken together, our results demonstrate that RuvC-catalyzed PCT ultimately generates DNA fragments with heterogeneous ends following initial DNA cleavage including a PAM proximal fragment with a blunt end and a PAM distal fragment with a staggered-end, 3'-recessed on the ntDNA strand. These kinetic and biochemical findings underline the importance of temporal control of Cas9 during gene-editing experiments and help explain the patterns of nucleotide insertions at sites of Cas9-catalyzed gene modification in vivo.

Dynamics changes of CRISPR-Cas9 systems induced by high fidelity mutations

CRISPR-Cas9, a powerful genome editing tool, has widely been applied in biological fields.

Structure and Dynamics of Cas9 HNH Domain Catalytic State

The bacterial CRISPR-Cas9 immune system has been harnessed as a powerful and versatile genome-editing tool and holds immense promise for future therapeutic applications. Despite recent advances in understanding Cas9 structures and its functional mechanism, little is known about the catalytic state of the Cas9 HNH nuclease domain, and identifying how the divalent metal ions affect the HNH domain conformational transition remains elusive. A deeper understanding of Cas9 activation and its cleavage mechanism can enable further optimization of Cas9-based genome-editing specificity and efficiency. Using two distinct molecular dynamics simulation techniques, we have obtained a cross-validated catalytically active state of Cas9 HNH domain primed for cutting the target DNA strand. Moreover, herein we demonstrate the essential roles of the catalytic Mg²⁺ for the active state formation and stability. Importantly, we suggest that the derived catalytic conformation of the HNH domain can be exploited for rational engineering of Cas9 variants with enhanced specificity.

CRISPR-Cas9 conformational activation as elucidated from enhanced molecular simulations

CRISPR-Cas9 has become a facile genome editing technology, yet the structural and mechanistic features underlying its function are unclear. Here, we perform extensive molecular simulations in an enhanced sampling regime, using a Gaussian-accelerated molecular dynamics (GaMD) methodology, which probes displacements over hundreds of microseconds to milliseconds, to reveal the conformational dynamics of the endonuclease Cas9 during its activation toward catalysis. We disclose the conformational transition of Cas9 from its apo form to the RNA-bound form, suggesting a mechanism for RNA recruitment in which the domain relocations cause the formation of a positively charged cavity for nucleic acid binding. GaMD also reveals the conformation of a catalytically competent Cas9, which is prone for catalysis and whose experimental characterization is still limited. We show that, upon DNA binding, the conformational dynamics of the HNH domain triggers the formation of the active state, explaining how the HNH domain exerts a conformational control domain over DNA cleavage [Sternberg SH et al. (2015) Nature, 527, 110-113]. These results provide atomic-level information on the molecular mechanism of CRISPR-Cas9 that will inspire future experimental investigations aimed at fully clarifying the biophysics of this unique genome editing machinery and at developing new tools for nucleic acid manipulation based on CRISPR-Cas9.