key: cord-0760047-cvjf172l
authors: Ke, Yuqing; Ghalandari, Behafarid; Huang, Shiyi; Li, Sijie; Huang, Chengjie; Zhi, Xiao; Cui, Daxiang; Ding, Xianting
title: 2′-O-Methyl modified guide RNA promotes the single nucleotide polymorphism (SNP) discrimination ability of CRISPR–Cas12a systems
date: 2022-02-01
journal: Chemical science
DOI: 10.1039/d1sc06832f
sha: 878d96689d674293577fa679a6ea04c62a7187b1
doc_id: 760047
cord_uid: cvjf172l

The CRISPR–Cas12a system has been widely applied to genome editing and molecular diagnostics. However, off-target cleavages and false-positive results remain as major concerns in Cas12a practical applications. Herein, we propose a strategy by utilizing the 2′-O-methyl (2′-OMe) modified guide RNA (gRNA) to promote the Cas12a's specificity. Gibbs free energy analysis demonstrates that the 2′-OMe modifications at the 3′-end of gRNA effectively suppress the Cas12a's overall non-specific affinity while maintaining high on-target affinity. For general application illustrations, HBV genotyping and SARS-CoV-2 D614G mutant biosensing platforms are developed to validate the enhanced Cas12a's specificity. Our results indicate that the 2′-OMe modified gRNAs could discriminate single-base mutations with at least two-fold enhanced specificity compared to unmodified gRNAs. Furthermore, we investigate the enhancing mechanisms of the 2′-OMe modified Cas12a systems by molecular docking simulations and the results suggest that the 2′-OMe modifications at the 3′-end of gRNA reduce the Cas12a's binding activity to off-target DNA. This work offers a versatile and universal gRNA design strategy for highly specific Cas12a system development.

The clustered regularly interspaced short palindromic repeats (CRISPR) and CRISPR-associated (Cas) proteins are unique for bacteria and archaea, constituting an adaptive immune system against foreign mobile genetic elements. 1,2 CRISPR-Cas systems are classied into Class 1 (use multiple Cas proteins) and Class 2 systems (use a single multidomain Cas protein), which are subdivided into six types (types I to VI) based on the complexity and feature proteins. 3 As a member of Class 2 systems, the type II-A CRISPR-Cas9 has been most extensively studied and developed as a genome editing and therapeutic tool. 4 Cas9 possesses two nuclease sites -His-Asn-His (HNH) and RuvClike domains to achieve the precise double-stranded DNA (dsDNA) cleavage at a specic site mediated by a dual-CRISPR RNA (crRNA) and trans-activating crRNA (transcrRNA) guide. 5, 6 The dual -transcrRNA:crRNA guide of Cas9 can also be engineered as a single RNA chimera fused by a short loop. 7 Another recently characterized Class 2 CRISPR-Cas system, the type V-A Cas12a (also known as Cpf1), is a single RNA (41-44 nt) guided endonuclease that can also recognize, unwind and induce a double-strand break (DSB) in the target DNA with high spec-icity. 8 Distinct from Cas9, Cas12a possesses a single nuclease pocket in the RuvC domain and has been shown to exhibit indiscriminate single-stranded DNase activity aer binding with target DNA. 9, 10 Besides, Cas12a shows a thymine (T)-rich (typically 5 0 -TTTV-3 0 , V ¼ A, C, or G) protospacer adjacent motif (PAM) preference, while Cas9 recognizes guanine (G)-rich sequences. [11] [12] [13] Based on the precise target DNA cleavage and indiscriminate ssDNA degradation activities, CRISPR-Cas12a systems are widely exploited for molecular diagnostics [14] [15] [16] [17] and gene editing. 12, 18 However, a recent research study using nuclease digestion and deep sequencing (NucleaSeq) to compare Streptococcus pyogenes (Sp)Cas9 and Acidaminococcus species (As)Cas12a cleavage activities revealed that Cas12a has similar cleavage specicities to Cas9 and tolerates similar mismatches in vitro, 19, 20 though Cas12a is reported to show less off-target effects in vivo. 21, 22 Cas12a proteins are more sensitive to mismatches at PAM-proximal positions 1-8 within the protospacer (target recognition sequence) but tolerable to PAMdistal mismatches, especially for single nucleotide polymorphisms (SNPs). 19, 23 Therefore, it's necessary to suppress the off-target cleavages of Cas12a in case of false positives or safety issues.

To increase the efficacy and specicity of CRISPR systems, developing high-delity Cas proteins 24, 25 and optimizing gRNA design [26] [27] [28] are two main strategies. The molecular structures and activation mechanisms of Cas12a have been claried, allowing the development of high-performance proteins. 23, 29, 30 Recently, an enhanced Acidaminococcus sp. Cas12a (enAs-Cas12a) variant has been engineered to expand the targeting range and increase the genome editing activity of Cas12a, and a high-delity version of enAsCas12a has also been developed to improve the specicity. 24, 31 Considering the safety issues and easy modications, (g)RNAs are more accessible to be engineered than proteins. Chemical modications 32-34 and hairpin structures 35 are applied to promote the stability of gRNAs and improve the editing efficacy of Cas proteins. For example, chemical modications (such as 2 0 -O-methyl, bridge nucleic acids, locked nucleic acids and 2 0 -O-methyl-3 0 -phosphonoacetate, etc.) incorporated at specic sites of Cas9-gRNAs can dramatically reduce off-target cleavages while maintaining high on-target performance and enhancing the Cas9 genome editing efficiency in human cells. 33, 36, 37 Besides, Cas12a-gRNAs with chemically modied terminal nucleotides are found to be exonuclease resistant while retaining gene-editing activity. 38 Studies show that chemical modications of RNA are effective to probe the RNA-protein relationships in ribonucleoprotein (RNP) enzymes 39 and increase the stability of RNA against degradation by RNase. 40 Moreover, the 2 0 -position modications of the sugar ring in RNA are benecial to enhance the stability and binding activity of RNA. 41 2 0 -O-Methyl (2 0 -OMe, also known as Nm) modications play an important role in regulating the gene expression at the post-transcriptional level in vivo [42] [43] [44] and can increase the stability of secondary structures of paired RNA. 45, 46 Due to the commercial availability and low cost, 2 0 -OMe modications have emerged as an attractive strategy to modulate RNA ensembles. 45 However, little is known about the effects of 2 0 -OMe modications for Cas12a-gRNA to improve the Cas12a's SNP discrimination ability in vitro, and the enhancing mechanisms of utilizing modied gRNAs for Cas12a systems remain to be elucidated yet. Herein, we propose a 2 0 -O-methyl modied gRNA (mgRNA) strategy to suppress the binding affinity between Lachnospiraceae bacterium ND2006 (Lba)Cas12a and off-target DNA, thus improving the overall on-target specicity of the CRISPR-Cas12a systems (mg-CRISPR). The enhanced system shows great potential for establishing highly specic molecular diagnostic tools such as pathogen genotyping platforms.

2 0 -OMe RNA modications at the 3 0 -end of gRNA suppressed Cas12a's affinity to off-target DNA while maintaining high affinity to on-target DNA For proof-of-concept, we rst selected HBV genotypes B and C as our targets. HBV remains a worldwide health concern and HBV genotyping is signicant for clinical therapy, whereas the B-type and C-type are major distributions in Oceania and Asia. 47 As shown in Fig. 1A , a 20 nt sequence of the B-type (on-target) containing two single nucleotide polymorphisms (SNPs) differing from the C-type (off-target) was determined as the target site. We designed an unmodied gRNA (ugRNA) and three various mgRNAs with 2 0 -OMe modications at different positions (5 0 -end for mgRNA-5 0 , the middle region for mgRNAm, and 3 0 -end for mgRNA-3 0 , respectively).

To evaluate the effects of 2 0 -OMe modications on the activity of Cas12a, we performed the Gibbs free energy analysis (Fig. 1B) . Fluorescence quenching methods have been widely exploited to study protein dynamics and interactions. 48, 49 Cas12a proteins have natural uorophores (tryptophan, tyrosine and phenylalanine residues), allowing the measurement of conformational changes through intrinsic uorescence spectroscopy. The Gibbs free energy alteration during Cas12a to Cas12a/gRNA complex transition is termed "DG 0 1 " (Fig. 1C , S1 and S2 †). For target DNA recognition, the Gibbs free energy alteration from Cas12a/gRNA to Cas12a/gRNA/on-target DNA is denoted as "DG 0 2-[on] " (Fig. S3 Fig. 1D and Table S1, † the DG value of mgRNA-3 0 (2.32 kJ mol À1 ) was higher than those of ugRNA (1.13 kJ mol À1 ), mgRNA-5 0 (1.24 kJ mol À1 ) and mgRNA-m (2.01 kJ mol À1 ), indicating that the 3 0 -end modications of mgRNA-3 0 suppressed the Cas12a's overall affinity to off-target DNA. As shown in Fig. 1D , the SG values of mgRNA-5 0 (À80.62 kJ mol À1 ) and mgRNA-3 0 (À80.09 kJ mol À1 ) were similar to that of ugRNA (À79.61 kJ mol À1 ). These observations demonstrated that mgRNA-5 0 and mgRNA-3 0 maintained high-level on-target affinity as ugRNA. According to the standard Gibbs free energy analysis, 2 0 -OMe modications at the 3 0 -end of mgRNA-3 0 suppressed the Cas12a's overall off-target affinity without sacricing the ontarget affinity.

We hypothesized that the Cas12a's affinity to on-target and off-target DNA may be associated with the sensitivity and specicity of the Cas12a system, so we established an HBV genotyping platform (Fig. 1E ). For signal readout, we used a single-stranded DNA linked quenched uorophore (5 0 -HEX-A (12) -BHQ1-3 0 ) as the reporter. The reactions were processed at 37 C for 1.5 h with uorescence intensity measured every 2 minutes. To evaluate the specicity of the Cas12a system, we calculated the discrimination factor (DF) using eqn (1):

where RFU [on] means the uorescence intensity for on-target DNA and RFU [off] for off-target DNA. As shown in Fig. 1F and Table 1 , for target DNA at high concentrations (10 nM and 1 nM), the DF values of mgRNA-m (6.54 AE 0.31 and 8.11 AE 0.18, respectively) and mgRNA-3 0 (5.77 AE 0.64 and 8.64 AE 0.81, respectively) were higher than those of ugRNA (2.93 AE 0.54 and 2.87 AE 0.18, respectively) and mgRNA-5 0 (3.2 AE 0.5 and 3.89 AE 0.77, respectively). For target DNA at a lower concentration (0.1 nM), the DF value of mgRNA-3 0 (8.52 AE 0.4) was higher than that of ugRNA (5.62 AE 0.61), mgRNA-5 0 (6.53 AE 0.11) and mgRNA-m (3.89 AE 0.28). These observations demonstrated that 2 0 -OMe modications at the 3 0 -end of gRNA increased the Cas12a's discrimination efficacy by at least two-fold compared to unmodied gRNA. Notably, for high target concentrations (10 nM and 1 nM), the DF values of ugRNA were reducing along with the nuclease assay time, which was caused by the increasing off-target signal and rapidly saturated on-target signal. For a lower target concentration (0.1 nM), the DF values of all gRNAs increased with the assay time due to the weak off-target signal and slowly increasing on-target signal (Fig. S7 †) . mgRNA-m showed high DF values for high concentration targets but decreased DF values for a low concentration target, which may be caused by the low on-target sensitivity of mgRNA-m. In summary, mgRNA-3 0 possessed the highest nal DF values at 90 minutes compared to ugRNA and other mgRNAs. The specicity result of mgRNA-3 0 was consistent with our Gibbs free energy analysis that mgRNA-3 0 suppressed the Cas12a's overall affinity to off-target DNA compared to other gRNAs.

Moreover, we investigated the sensitivity of the Cas12a system. Various concentrations of on-target DNA (0 to 10 nM) were used as input DNA. According to the sensitivity results shown in Fig. 1G , mgRNA-5 0 and mgRNA-3 0 possessed a low LOD of 10 pM like ugRNA, while mgRNA-m showed higher LOD (100 pM) than ugRNA, which could explain the decreased DF values of mgRNA-m at a lower target DNA concentration (Fig. 1F ). These observations were consistent with our Gibbs free energy analysis that mgRNA-3 0 maintained the Cas12a's high-level on-target affinity as ugRNA, while mgRNA-m decreased the overall on-target affinity compared to ugRNA. Our results illustrated that the DG values were positively correlated with the Cas12a's overall specicity while the SG values were positively correlated with the Cas12a's sensitivity.

2 0 -OMe modications at the 3 0 -end of gRNA increased the single nucleotide polymorphism (SNP) discrimination efficacy of the Cas12a system

To verify the feasibility of 2 0 -OMe modication at the 3 0 -end of gRNA to improve the SNP discrimination ability of the Cas12a system, we selected another two target sequences in HBV B-type (on-target) and C-type (off-target) genomes containing only one SNP. As shown in Fig. 2A , target-11th contains just one SNP at position 11, while target-19th has single SNP at position 19. Accordingly, we designed another two unmodied gRNAs (ugRNA-11th for target-11th and ugRNA-19th for target-19th) and related modied gRNAs (mgRNA-11th for target-11th and mgRNA-19th for target-19th). Similar to the targets described in Fig. 1A , we modied gRNAs at different positions: the 5 0 -end for mgRNAs-5 0 , middle region for mgRNAs-m, and 3 0 -end for mgRNAs-3 0 , respectively. Three adjacent RNAs were modied for each mgRNA (Fig. 2B) .

At rst, we investigated the specicity of the Cas12a system based on unmodied and modied gRNAs. The specicity results are provided in Fig. 2C and Table 1 . For target-11th DNA at 10 nM, mgRNA-11th-3 0 (DF value of 2.46 AE 0.09) increased the discriminating ability of Cas12a by two fold compared to ugRNA-11th (1.13 AE 0.03). For target-11th DNA at 1 nM, mgRNA-11th-3 0 (4.09 AE 0.18) improved the Cas12a's specicity by about three-fold compared to ugRNA-11th (1.41 AE 0.05). For target-11th DNA at 0.1 nM, mgRNA-11th-3 0 (4.85 AE 0.66) showed a two-fold higher specicity than ugRNA-11th (2.57 AE 0.11). These results illustrated that the 2 0 -OMe modications at the 3 0end of gRNA enhanced the Cas12a's specicity even for single SNP discrimination. Besides, for target-11th DNA at various concentrations (10 nM, 1 nM and 0.1 nM), mgRNA-11th-5 0 (1.64 AE 0.09, 2.07 AE 0.09 and 2.85 AE 0.08, respectively) and mgRNA-11th-m (1.48 AE 0.02, 2.22 AE 0.26 and 3.66 AE 0.16, respectively) showed a slightly improved specicity of about 1.5 fold compared to ugRNA-11th. These observations demonstrated that 2 0 -OMe modications at the 5 0 -end and middle region of gRNA caused minor effects on the Cas12a's activity for target discrimination with only one SNP at the middle position of the spacer (position 11 in a total of 20 nucleotides).

Moreover, we evaluated the Cas12a's activity for target-19th using ugRNA-19th and mgRNAs-19th. According to the results shown in Fig. 2C (bottom) and Table 1 , mgRNA-19th-3 0 (6.33 AE 0.3) showed about two-fold specicity compared to ugRNA-19th . These results indicated that 2 0 -OMe modications at the middle region could promote the Cas12a's specicity for targets containing one SNP proximal to the 3 0 -end of gRNA. Next, we tested the sensitivity of the Cas12a system for target-11th and target-19th. As shown in Fig. 2D , mgRNA-11th-3 0 showed a lower LOD of 1 pM than ugRNA-11th, mgRNA-11th-5 0 and mgRNA-11th-m (50 pM, 10 pM and 10 pM, respectively). For target-19th, mgRNA-19th-3 0 showed a similar LOD of 50 pM to ugRNA-19th. These results demonstrated that 2 0 -OMe modications at the 3 0 -end of mgRNA wouldn't sacrice the sensitivity of the Cas12a system even for various target sites containing different SNP locations (at position 11 or 19).

Furthermore, to explore the performance of the Cas12a system based on ugRNAs and mgRNAs for detection in complex samples, we mixed HBV on-target DNA with bacterial gDNA and tested the nal uorescence intensity. The concentrations of HBV target DNA were 1 nM, while the amounts of bacterial gDNA were 80 ng for each test. As shown in Fig. 2E , all mgRNAs and ugRNAs for the three target sites possessed high stability and specicity even in complex samples. Then, to validate the versatility of the 2 0 -OMe modications to improve the Cas12a's specicity for different applications, we developed a SARS-CoV-2 D614G mutant detection platform (Fig. 3A) . The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has caused a worldwide pandemic and the mutant strain with one amino acid change at position 614 (D614G) has been reported to be dominant. 50, 51 As shown in Fig. 3B , we selected the 20 nt sequence containing the D614G mutant from the S1 subunit of the SARS-CoV-2 genome as the target site. Accordingly, we designed one unmodied gRNA complementary to the SARS-CoV-2 D614G mutant (ugRNA-D614G). Three modied gRNAs (mgRNAs) with 2 0 -OMe modications at various positions were also designed (5 0 -end modications for mgRNA-D614G-5 0 , the middle region for mgRNA-D614G-m and 5 0 -end for mgRNA-D614G-3 0 , respectively). At rst, we tested the specicity of the Cas12a based SARS-CoV-2 biosensing platform. The results are provided in Fig. 3C . For target DNA at 1 nM and 0.1 nM, the nal DF values of mgRNA-D614G-3 0 (1.43 AE 0.08 and 1.8 AE 0.1 for 1 nM and 0.1 nM, respectively) were about twofold higher than those of ugRNA-D614G (0.94 AE 0.04 and 1.04 AE 0.08, respectively), illustrating that 2 0 -OMe modications improved the Cas12a's specicity for D614G mutant detection even at low target concentrations. Moreover, to investigate the sensitivity of the Cas12a based D614G mutant detection system, we used D614G mutant DNA (concentrations varying from 0 to 10 nM) as input DNA. As shown in Fig. 3D , the LOD of mgRNA-D614G-3 0 (10 pM) was lower than that of ugRNA-D614G, mgRNA-D614G-5 0 and mgRNA-D614G-m (50 pM), indicating that 2 0 -OMe modications at the 3 0 -end of gRNA enhanced the Cas12a's sensitivity for D614G mutant detection. Furthermore, to explore the performance of the Cas12a system in a complex sample, we mixed the SARS-CoV-2 D614G mutant DNA (200 pM per reaction) with HBV B-type DNA (600 pM per reaction) and analysed the uorescence intensity readout. According to the results shown in Fig. 3E , all unmodied and modied gRNAs-D614G possessed high specicity even in the complex sample, demonstrating the stability of the Cas12a based D614G mutant biosensing platform.

To investigate the mechanisms of the increased Cas12a's specicity based on the 2 0 -OMe modied gRNA, molecular docking simulations were further performed to illustrate the interactions between Cas12a, gRNA and off-target DNA. Herein, we selected unmodied gRNA (ugRNA), mgRNA-3 0 (with 2 0 -OMe modications at the 3 0 -end) and the target site as shown in Fig. 1 for demonstrations (Fig. 4A ). Protein Data Bank (PDB) ID-5XUS was selected for Cas12a. The structure of the Cas12a protein is provided in Fig. 4B . The docking result of Cas12a binding with ugRNA is shown in Fig. 4C , while that of mgRNA-3 0 is provided in Fig. 4D . In the binary system of Cas12a interacting with gRNAs, the docking analysis showed that gRNAs are located in the groove of the REC lobe (consisting of REC1 and REC2 domains) and the NUC lobe (WED and PI domains) of the Cas12a protein ( Fig. 4C and D) . The 3 0 anks of ugRNA and mgRNA-3 0 are oriented to the REC1 domain surface, close to the channel between the REC2, RuvC1 and Nuc domains. The 5 0 anks of gRNAs are oriented to the groove of the REC2, WED and PI domains, close to the BH domain. Previous studies 11, 23 suggest that gRNA distribution on the Cas12a's surface is signicant for the subsequent catalysis processes, since the REC, RuvC, WED, PI and BH domains act as a clamp to trap gRNA. The theoretical binding energies of Cas12a interacting with ugRNA and mgRNA-3 0 are À356.45 kJ mol À1 and À361.99 kJ mol À1 , respectively, demonstrating that mgRNA-3 0 possesses higher affinity to the Cas12a protein than ugRNA. To better elucidate the difference between ugRNA and mgRNA-3 0 interacting with Cas12a, we performed the docking calculations of ugRNA and mgRNA-3 0 in a single gure (Fig. 4E) . The green strand refers to ugRNA while the blue strand represents mgRNA-3 0 . As shown in Fig. 4E (bottom) , the 3 0 ank of mgRNA-3 0 is oriented towards penetration Cas12a protein, while the 3 0 ank of ugRNA gets more freedom on the REC1 surface. In other words, there is a greater tendency for mgRNA-3 0 to form the Cas12a/mgRNA-3 0 complex compared to ugRNA due to the high stability of the Cas12a/mgRNA-3 0 structure. We believe the difference between the Cas12a/ugRNA and Cas12a/mgRNA-3 0 structures are caused by the 2 0 -OMe modications at the 3 0 -end of mgRNA-3 0 .

Moreover, we performed the docking calculations between the Cas12a/gRNA complex and off-target DNA in the ternary system. As shown in Fig. 4F and G, the gRNAs recognize the 5 0 ank of target DNA through the 3 0 anks of gRNAs. The calculated binding energies of Cas12a/ugRNA/off-target DNA and Cas12a/mgRNA-3 0 /off-target DNA complexes are À388.45 kJ mol À1 and À383.09 kJ mol À1 , respectively, suggesting that the Cas12a/mgRNA-3 0 complex suppressed the affinity to off-target DNA compared to Cas12a/ugRNA. Besides, the relative distance between the off-target DNA and the surface of REC1 in the Cas12a/mgRNA-3 0 complex (21 angstroms) increases by 5 angstroms compared to that of the Cas12a/ ugRNA complex (16 angstroms) . In other words, the activity between Cas12a/mgRNA-3 0 /off-target DNA is weaker than that of Cas12a/ugRNA/off-target DNA due to the structural differences between ugRNA and mgRNA-3 0 . 

In this study, we proposed an enhanced CRISPR-Cas12a system by utilizing 2 0 -O-methyl modied gRNA, termed mg-CRISPR, to promote the specicity of the Cas12a protein.

Three main improvements were achieved. Firstly, modied gRNAs with 2 0 -OMe modications at various positions of the gRNAs were developed to increase the specicity of the Cas12a system. The enhanced Cas12a system suppressed the off-target binding activity while maintaining high on-target affinity compared to unmodied gRNAs. Secondly, HBV genotyping and SARS-CoV-2 D614G mutant biosensing platforms were established to validate the high specicity and versatility of the Cas12a system. 2 0 -OMe modied gRNAs increased the Cas12a's specicity by up to three fold with high sensitivity (10 pM) compared to unmodied gRNAs. Besides, the enhanced Cas12a system was capable of single SNP discrimination. Current technologies for SNP analysis mainly include enzymatic methods, 52 uorescence-based assays, 53 PCR-based detections 54,55 and electrochemical methods. 56 Existing SNP detection methods mostly suffer from complex procedures, expensive instruments or low specicity. Moreover, CRISPRbased methods for SNP analysis are capable of sensitive, specic, programmable and rapid detection. 16 By combining the ordinary CRISPR system with the 2 0 -OMe modied gRNA studied in this research, the discrimination factors have been increased by three fold for single SNP detection. Thirdly, docking simulations were performed to investigate the mechanisms of 2 0 -OMe modied gRNA suppressing the non-specic activity of the Cas12a protein. Our results demonstrate that 2 0 -OMe modications of gRNA are versatile strategies to improve the Cas12a's specicity for diverse applications. 

PCR primers and the uorophore quencher-labelled ssDNA reporter were ordered from Sangon Biotechnology Co., Ltd (Shanghai, China). Primers and ssDNA uorescence reporters used in this study are listed in Table S2 . † PCR reactions were performed following instructions of OneTaq Quick-Load 2X Master Mix with Standard Buffer from New England Biolabs (Ipswich MA, United States). All unmodied and 2 0 -OMe-modied gRNAs were synthesized by Shanghai DNA Bioscience Co., Ltd (https://dnabio.cn/, Shanghai, China). The sequences of gRNAs are provided in Table S3 . † HBV plasmids containing the whole genome of genotype B or C were provided by Prof. Yumei Wen (Fudan University, Shanghai, China). SARS-CoV-2 plasmids containing the wild-type and D614G mutant spike gene (catalog numbers: GS-200519_A001 and GS-200519_A012, respectively) were ordered from GENEWIZ Biotechnology Co., Ltd (Suzhou, China). The detailed information of HBV and SARS-CoV-2 plasmids used in this study are provided in Table S4 . † The PCR products were puried using a MinElute PCR purication kit from Qiagen (Hilden, Germany). The LbaCas12a protein and NEBuffer 2.1 were purchased from New England Biolabs (Ipswich MA, United States). RNase-free water was ordered from Thermo Fisher (Waltham, United States).

PCR experiments were carried out on a LightCycler 96® (Roche). The real-time uorescence intensity of Cas12a collateral cleavage assay was measured by the uorescence channel of a LightCycler 96® (Roche, Switzerland). Electrophoresis was conducted using a Tanon EPS-100 electrophoresis system. Fluorescence of the Cas12a protein was obtained by using an Agilent Technologies Cary Eclipse uorescence spectrometer.

Firstly, HBV and SARS-CoV-2 plasmids were amplied by PCR using a OneTaq Quick-Load 2Â Master Mix. PCR primers for HBV and SARS-CoV-2 targets were designed by using NCBI Primer -Blast. Secondly, the amplied products were evaluated by agarose gel electrophoresis and then puried by using a MinElute PCR purication kit (Qiagen, German). Finally, the concentration of puried target DNA was evaluated by a Nanodrop (ALLSHENG, China).

The SNP sites for HBV B-type DNA and HBV C-type DNA were searched using the NCBI Nucleotide -Blast (BLAST: Basic local alignment search tool (nih.gov)), then the sequences containing both PAM (5 0 -TTTN) and SNP sites were determined to be the candidate protospacers for HBV gRNA. Finally, three HBV target sequences were selected in this study. For the HBV target used in Fig. 1 , B-type (on-target) DNA contains two SNPs differing from C-type (off-target) DNA. HBV target-11th contains one SNP at position 11, and HBV target-19th has one SNP at position 19. HBV gRNAs are designed according to the selected target sequences, the gRNA spacers are complementary with the HBV B-type DNA. For SARS-CoV-2 gRNA design, the target sequence containing both the PAM (5 0 -TTTN) and the D614G mutant site is selected as the protospacer. For both HBV and SARS-CoV-2 gRNAs, chemically modied gRNAs were designed by inserting three adjacent 2 0 -OMe modications at the 5 0 -end, middle region and 3 0 -end of gRNAs. The detailed information of all gRNAs is provided in Table S3 . †

Firstly, the LbaCas12a protein (2.5 mM) was incubated with gRNA (concentrations varying from 0 to 1.25 mM) in a cuvette. Then, the intrinsic uorescence spectra of the Cas12a/gRNA complex were measured using a uorescence spectrometer (Cary Eclipse, Agilent, Australia). Aerwards, the Cas12a/gRNA complex was incubated with HBV DNA (0-200 nM) and emission spectra of the complex were recorded. The wavelength of excitation was set at 285 nm and the scan range of emission wavelength was 295-500 nm. All measurements were conducted at 25 C and 37 C.

Fluorescence quenching methods were used to analyse the binding energetics of Cas12a-gRNA and Cas12a-gRNA-target DNA complexes. The quenched uorescence is described using the Stern-Volmer equation:

where F 0 and F are the intrinsic uorescence intensities of the Cas12a protein in the absence and presence of gRNA or target DNA, respectively. K SV is the Stern-Volmer quenching constant, [Q] is the concentration of gRNA or target DNA, k q is the biomolecular quenching rate constant, and s 0 is the average lifetime of the excited uorophore in the absence of a quencher (average s 0 ¼ 10 À8 s). 57 The exact parameters of interaction were calculated by using the modied Stern-Volmer equation:

in which f a is the mole fraction of accessible uorescence. According to the quenching mechanism, the double logarithmic equation was used to obtain thermodynamic parameters for the interactions of the Cas13a protein with crRNAs and HBV RNA:

in which K a and n mean the binding constant and the number of binding sites, respectively. 58 Assuming unchanged entropy (DS 0 ) and enthalpy (DH 0 ) at different temperatures, the type of interactions between the Cas12a protein and gRNAs or target DNA can be calculated using the Van 't Hoff equation:

where R and T are the gas constant and absolute temperature, respectively. And changes in standard Gibbs free energy (DG 0 ) can be calculated using the following equation:

Cas12a collateral cleavage assay

The LbaCas12a based uorescence assay was conducted with a 50 nM LbaCas12a protein, 22.5 nM gRNA, a 100 nM uorophore quencher-labelled ssDNA reporter and various amounts of target DNA in 1 Â NEB2.1 buffer (containing 5 mM NaCl, 1 mM Tris-HCl, 1 mM MgCl 2 , and 10 mg mL À1 BSA, pH 7.9). The total volume was 20 mL for a single reaction. Assays were proceeded for 1.5 h at 37 C. Meanwhile, the uorescence intensity of each assay was collected every minute using the uorescence channel of a qPCR instrument (Roche, Light-Cycler® 96, Switzerland).

Molecular docking study of the Cas12a/gRNA and Cas12a/ gRNA-target DNA complex

The HDOCK webserver (http://hdock.phys.hust.edu.cn/) was used to perform the docking calculations. 59, 60 The 3D structure of the LbaCas12a protein was acquired from Protein Data Bank ID 5XUS. The RNAComposer server (http://rnacomposer.cs.put.poznan.pl/) was utilized to generate the 3D structures of gRNAs. 61 The 3D structure of HBV target DNA was obtained using the HyperChem v8.0.6 program. The gRNAs and target DNA geometries were optimized using the CHARMM36 force eld implemented in the NAMD 2.10 package. 62, 63 The gRNAs and target DNA structures were immersed within 60Å and 80Å cubic boxes of TIP3P water molecules and neutralized by 25 and 35 Na + ions, respectively. Classical molecular dynamics simulations were performed considering calculation parameters based on preparation instructions. 64 Next, docking calculations were conducted to search the complex formation between the Cas12a protein and guide RNA (binary system), followed for the Cas12a/gRNA complex and target DNA (ternary system). The Visual Molecular Dynamics (VMD) package was used to investigate the docking results. 65 The optimal position of each molecule was analyzed according to the energy score and the condition of interaction was selected based on the minimum energy HDOCK score. We obtained the binding energy using the PM7 semiempirical method implemented in the MOPAC 2016 package. 66, 67 Statistical analysis Statistical analysis was performed by one-way ANOVA-multiple comparisons using GraphPad Prism 8.0. All assays were repeated three times. The results were presented as mean AE S.D. 

There are no conicts to declare.

All data associated with the research are available from the Royal Society of Chemistry online or from the author.

Principles of uorescence spectroscopy

In vitro study on the alterations of brain tubulin structure and assembly affected by magnetite nanoparticles

MOPAC2016