Hairpin‐Spacer crRNA‐Enhanced CRISPR/Cas13a System Promotes the Specificity of Single Nucleotide Polymorphism (SNP) Identification

Abstract The Cas13a system has great potential in RNA interference and molecular diagnostic fields. However, lacking guidelines for crRNA design hinders practical applications of the Cas13a system in RNA editing and single nucleotide polymorphism identification. This study posits that crRNAs with hairpin spacers improve the specificity of CRISPR/Cas13a system (termed hs‐CRISPR). Gibbs free energy analysis suggests that the hairpin‐spacer crRNAs (hs‐crRNAs) suppress Cas13a's affinity to off‐target RNA. A hepatitis B virus DNA genotyping platform is established to further validate the high‐specificity of hs‐CRISPR/Cas13a system. Compared to ordinary crRNA, hs‐crRNAs increase the specificity by threefold without sacrificing the sensitivity of the CRISPR/Cas13a system. Furthermore, the mechanism of the Cas13a/hs‐crRNA/target RNA composition is elucidated with theoretical simulations. This work builds on the fundamental understanding of Cas13a activation and offers significant improvements for the rational design of crRNA for the CRISPR/Cas13a system.


Introduction
Clustered regularly interspaced short palindromic repeats (CRISPR) and CRISPR-associated (Cas) protein systems offer various possibilities in molecular diagnostics [1][2][3] and RNA interference [4][5][6] beyond genome editing [7,8] due to the programmability of nuclease activity. The activity of CRISPR/Cas13a, an RNA-guided RNase, is determined by higher eukaryotes and prokaryotes nucleotide-binding (HEPN) domains. [9,10] Cas13a can combine with CRISPR RNA (crRNA) to form the Cas13a/crRNA duplex, which is competent at recognizing target RNA with the guidance of crRNA-spacer. Once assembled with target RNA, the Cas13a/crRNA/target RNA ternary complex will activate the HEPN domains to create a general RNase active site in the Cas13a protein. [11,12] These properties impart CRISPR/Cas13a with advantages of rapid, accurate, highsensitivity, and programmability for molecular diagnostics. [13][14][15] However, CRISPR/Cas systems are still challenged by high offtarget rates that affect the accuracy and stability of practical applications. [5,16,17] Optimizing Cas proteins [18,19] and/or crRNA design [20][21][22] are prominent strategies that increase efficacy. The clarification of activation mechanisms via probing structures helps the development of high-performance Cas13a proteins. [9,23,24] The most common method to optimize crRNA is to introduce synthetic mismatch into the crRNA-spacer, which promotes the identification ability of single nucleotide polymorphisms (SNPs) toward Cas13a proteins. [15,20,21,25] However, Cas13a protein optimization remains challenging due to its sophisticated nature. Besides, the location of synthetic mismatch in crRNA is uncertain and its enhancing mechanism is unclear. In recent years, studies show that hairpin structures enhance the functions of RNA molecules in extensive applications, such as strand-displacement amplification, [26] RNA interference, [27] aptamer biosensing, [28][29][30] and RNA origami. [31,32] Furthermore, RNA secondary structures increase the specificity of Cas9 and Cas12a systems in genetic engineering. [33] However, little is known regarding the value of hairpin-spacer structures in the CRISPR/Cas13a system while applied to molecular diagnostics and RNA interference.
Herein, we propose a hairpin-spacer crRNA (hs-crRNA) strategy that improves the specificity of the CRISPR/Cas13a system (hs-CRISPR) by suppressing the binding affinity between Cas13a and off-target RNA ( Figure 1A). First, we investigated the enhancing mechanism of hs-CRISPR and analyzed the variation of Gibbs free energy between the interaction of Cas13a with hairpinspacer crRNA and ordinary crRNA. Next, we evaluated the practical efficacy of hs-CRISPR at genotyping hepatitis B virus (HBV) DNA ( Figure 1B). Last, we examine the structural information of the Cas13a/hs-crRNA duplex and Cas13a/hs-crRNA/target-RNA complex via simulation.

Assay Scheme of hs-CRISPR
The principle of hs-CRISPR is described in Figure 1. The main basis of this system is the target RNA recognition of Cas13a/crRNA duplex and the resulting collateral trans-cleavage ( Figure 1A). At first, Cas13a is in the ground state (GS) and shows no nuclease activity. Then it can bind with crRNA and step into the transient state (TS). Cas13a/crRNA duplex in TS can be guided by crRNA-spacer to recognize specific target RNA and shows RNase activity, being in the active state (AS). While Cas13a/crRNA interacts with off-target RNA, Cas13a will be in undesired transient state (uTS) with no nuclease activity. To investigate the effects of hairpin crRNA-spacer structures on Cas13a activity, we designed crRNAs with or without engineered hairpin-spacer. Hairpin structures were created by extending the 3′-flank of crRNA with several bases complementary to spacer sequence. Basic hairpins are composed of stems and loops. By tuning the location and length of stems, we obtained different hairpin-spacer crRNAs (hs-crRNAs). Furthermore, we measured the Gibbs free energy changes from GS to TS (termed "ΔG 0 1 ") and TS to AS ("ΔG 0 2-[on] ") or TS to uTS (termed "ΔG 0 2-[off] ") using fluorescence spectrometer. The difference in ΔG 0 2-[on] and ΔG 0 2-[off] was determined as "ΔG". For hs-CRISPR, crRNA with low free energy was used and therefore high ΔG was obtained, making off-target RNA unable to activate the nuclease activity of Cas13a.
To verify the advanced discrimination ability of hs-CRISPR on SNPs, we further compared hs-crRNAs with ordinary crRNA (o-crRNA) for HBV genotyping ( Figure 1B). Considering HBV genotypes B and C are major subpopulations in Oceania and Asia, we selected HBV B-type and C-type as our genotyping targets for proof-of-concept in this study. HBV DNA are amplified by recombinase polymerase amplification (RPA) reaction. Then amplified products are transcribed into related RNAs using in vitro T7 transcription system, which function as targets of Cas13a/crRNA duplex. Notably, owing to the nucleotide polymorphism caused by different genotypes, HBV C-type cannot activate the RNase activity of LwaCas13a binding with HBV B-type crRNA (B crRNA).

Hs-crRNAs Show Lower Affinity to Cas13a than Ordinary crRNA
To explore the effects of hairpin crRNA-spacers on Cas13a activity, we compared the capabilities of crRNAs with or without hair-pin spacer for HBV DNA genotyping. As shown in Figure 2A, a 28-nt sequence of HBV containing two SNPs was chosen as the target of o-crRNA. By extending the 3′-end of o-crRNA to form a 10-bp stem structure at the distal region to crRNA-repeat, we obtained hs-crRNA1. By changing the stem location and length of hairpin structure in crRNA-spacer, we further designed four different hs-crRNAs, named hs-crRNA2 to hs-crRNA5, respectively. The stem structures of hs-crRNA2 located at the proximal region to crRNA-repeat while that of hs-crRNA3 showed an 8nt distance from crRNA-repeat. The stem lengths of hs-crRNA4 and hs-crRNA5 were 5-bp and 15-bp, respectively. The sequences of all crRNAs are supplied in Table S1 of the Supporting Information. Furthermore, we predicted the minimum free energy structures of crRNA-spacers using RNAstructure online tools and RNAfold webserver ( Figure S1, Supporting Information). We discovered that hs-crRNAs had lower predicted minimum free energies (−68.20, −74.48, −43.5136, −23.012, and −111.7128 kJ mol −1 for hs-crRNA1 to hs-crRNA5, respectively) than o-crRNA (−3.68 kJ mol −1 ).
To clarify the mechanism of Cas13a binding with crRNA, we investigated the conformational changes of Cas13a by fluorescence spectroscopy which has been widely applied in protein dynamics, assembly and interaction researches. [34][35][36] The Cas13a protein has natural fluorophores (tryptophan and tyrosine residues), allowing fluorescence analysis of conformational changes. First, we incubated Cas13a (GS) with different concentrations of crRNA ranging from 0 × 10 −6 to 2.5 × 10 −6 m at 25 and 37°C for 10 min, respectively ( Figure 2B). Then we measured the fluorescence intensity of Cas13a/crRNA duplex (TS) and analyzed the fluorescence quenching data using the Stern-Volmer equation (Equation 1) where F 0 and F are the fluorescence intensity of Cas13a protein in the absence and presence of RNA, respectively. K SV is the Stern-Volmer quenching constant, [Q] is the concentration of RNA, k q is biomolecular quenching rate constant, and 0 is the average lifetime of excited fluorophore in the absence of quencher (average 0 = 10 −8 s). [37] A modified Stern-Volmer equation was used to calculate the exact parameters of interaction (Equation 2) in which f a is the mole fraction of accessible fluorescence and K SV is the Stern-Volmer quenching constant of accessible fluorophores. This equation can be applied to determine the percentage of quencher accessible fluorophores. [38] The fluorescence spectra of Cas13a incubating with crRNAs (0 × 10 −6 -2.5 × 10 −6 ) and the Stern-Volmer plots are displayed in Figures S2 and S3 of the Supporting Information, respectively. To investigate the types of interactions between Cas13a and crRNA, we further calculated the entropy (ΔS 0 ) and enthalpy (ΔH 0 ) of Cas13a/crRNA duplexes using Equation (3) log . Cas13a in TS can recognize specific target RNA guided by crRNA and shows nuclease activity, entering the active state (AS). crRNA consists of direct repeat (dot in green) and spacer (dot in purple) structures. We designed three crRNAs with or without hairpin spacer. Basic hairpins are composed of stems and loops. The design parameters include location and length of stems. The Gibbs free energy changes from GS to TS (termed "ΔG 0 1 ") and TS to AS (termed "ΔG 0 2 ") represent the binding affinity of Cas13a to crRNA and target RNA, respectively. ΔG 0 1 for the hs-CRISPR/Cas13a system is labeled as ΔG 0 1 ′. The difference in binding affinity of Cas13a to on-target and off-target RNA was labeled as "ΔG". ΔG for the hs-CRISPR/Cas13a system is labeled as ΔG′. B) Workflow of HBV genotyping platform based on hs-CRISPR. HBV DNA is amplified by RPA reaction, then transcribed into RNA, which can be recognized by Cas13a/crRNA duplex. Notably, due to the nucleotide polymorphisms caused by genotypes, collateral cleavage will not occur for off-target RNA.
in which K a and n mean the binding constant and the number of binding sites, respectively. [39] Assuming unchanged entropy (ΔS 0 ) and enthalpy (ΔH 0 ) at different temperatures, the type of interactions between Cas13a protein and crRNAs or HBV RNA can be investigated using Van 't Hoff equation (Equation 4) where R and T are gas constant and absolute temperature, respectively. Changes of Standard Gibbs free energy ΔG 0 1 were calculated using Equation (5) According to the thermodynamic results provided in Figure 2C and Table S2 (Supporting Information), the ΔG 0 1 value of o-crRNA (−39.31 kJ mol −1 ) was lower than hs-crRNAs (−37.46, −37.64, −38.64, −38.83, and −39.12 kJ mol −1 for hs-crRNA1 to hs-crRNA5, respectively) at 37°C, demonstrating that ordinary crRNA possessed higher binding affinity toward Cas13a protein than hairpin-spacer crRNAs. We believe that this phenomenon is caused by the RNA structure differences, since the molecular structure is one of the most critical factors to determine the affinity. It should be noted that both the structures of RNA molecules and the amounts of Cas13a's space occupied by crRNA would affect the affinity of the protein. To further investigate the effects of Cas13a's lower binding affinity on Cas13a/crRNA duplex formation, we incubated the Cas13a with o-crRNA and hs-crRNA1 for different time (0-30 min) and measured the fluorescence spectra of the duplex. As shown in Figure S4 of the Supporting Information, longer incubation time would not affect the fluorescence spectra of the duplex, indicating that both o-crRNA and hs-crRNA1 interacted with the Cas13a and formed the Cas13a/crRNA duplex spontaneously. Notably, the shapes of fluorescence spectra of the Cas13a/o-crRNA and Cas13a/hs-crRNA1 duplexes were different, suggesting that certain areas of Cas13a were occupied by o-crRNA and hs-crRNA, which may also influence the Cas13a's binding affinity to RNA.

Hs-crRNAs Suppress the Cas13a's Overall Affinity to Off-Target RNA
To investigate the interactions between Cas13a/crRNA duplex and target RNA, we further incubated Cas13a/crRNA duplex with various concentrations of on-target or off-target RNA (0 × 10 −6 -1 × 10 −6 m) ( Figure 2B). Similar to analyzing the interaction between Cas13a and crRNA, we recorded the fluorescence intensity of Cas13a/crRNA/target RNA complex (AS) and analyzed  Figure S5 of the Supporting Information. Standard Gibbs free energy change ΔG 0 2 represents the binding affinity of Cas13a/crRNA duplex to on-target or off-target RNA. ΔG 0 2-[on] means the ΔG 0 2 value of Cas13a/crRNA interacting with on-target RNA, whereas ΔG 0 2-[off] represents off-target RNA. As shown in Figure 2C and for hs-crRNA1 to hs-crRNA5, respectively) were lower than o-crRNA (−39.06 kJ mol −1 ), suggesting that hairpin-spacer crRNAs increased the binding affinity of Cas13a/hs-crRNA duplex to ontarget RNA. Accordingly, we attributed the improved binding affinity to optimal conformation of Cas13a induced by hairpinspacer structures of hs-crRNAs.
Based on the Gibbs free energy analysis, we summarized the thermodynamic results of ΔG 0 1 and ΔG 0 2 in Table 1. We observed that the sum of ΔG 0 1 and ΔG 0 2-[on] values for all crRNAs (<−77.48 kJ mol −1 ) were lower than that of ΔG 0 1 and ΔG 0 values (> −76.95 kJ mol −1 ), indicating that the Cas13a protein possessed an overall high-level affinity to on-target RNA to be in the active state (AS). Notably, the difference of Cas13a's binding affinity to on-target and off-target RNA was determined as "ΔG" ). As shown in Table 1, hs-crRNAs showed an increased ΔG values (2.38, 2.31, 1.6, and 1.54 kJ mol −1 for hs-crRNA1, hs-crRNA2, hs-crRNA3, and hs-crRNA5, respectively) compared to o-crRNA (1.42 kJ mol −1 ). This observation indicated that hs-crRNAs suppressed the Cas13a's overall affinity to off-target RNA, making the Cas13a/hs-crRNA/off-target RNA complex unable to be in the AS due to the insufficient Gibbs

Hs-crRNAs Increase the Specificity of CRISPR/Cas13a System
To investigate the relationship between ΔG and the specificity of Cas13a system, we applied ordinary crRNA and hairpin-spacer crRNAs for HBV genotyping. The workflow of HBV genotyping is shown in Figure 3A. All reagents, including RPA mix, T7 transcription mix, and Cas13a nuclease assay reagents, were encapsulated in a single pellet. The whole detection process takes less than 2 h at 37°C. At first, we tested the specificity of the hs-CRISPR/Cas13a system as well as ordinary CRISPR/Cas13a system for HBV genotyping by utilizing hs-crRNAs and o-crRNA, respectively. As shown in Table 2 and Figure 3B, the on-target signals for Cas13a/hs-crRNA1, Cas13a/hs-crRNA2 to Cas13a/hs-crRNA5 (1297. 28 where F.I. [on] and F.I. [off] mean the end-point fluorescence intensity of on-target and off-target signals, respectively. As shown in According to our experimental results of all hs-crRNAs to on-target and off-target HBV, we propose a hairpin structure design strategy to enhance the specificity of the hs-CRISPR/Cas13a system. By comparing the results of hs-crRNA1 and hs-crRNA2 with hs-crRNA3, we concluded that the SNP-induced mismatch should be contained in the stem structure of hs-crRNA to efficiently enhance the Cas13a's specificity ( Figure 3D,E; Figure S6, Supporting Information). According to the results of hs-crRNA1 (10-bp stem), hs-crRNA4 (5-bp stem), and hs-crRNA5 (15-bp stem), the stem length of hairpin structure should be set at a proper range which was also determined by the location of the stem and the target sequence. For our experiments, the optimal stem lengths were 8-bp to 10-bp. Since hs-crRNA1 and hs-crRNA2 possessed the best performance in HBV genotyping, we further explored the sensitivity of hs-crRNA1 and hs-crRNA2 as well as o-crRNA. The limits of detection (LOD) of o-crRNA, hs-crRNA1, and hs-crRNA2 were all 10 0 × 10 −18 m (100 copies mL −1 ) ( Figure 3C-E). These results manifested that hairpin-spacer crRNAs improved the specificity of Cas13a without sacrificing the sensitivity. Moreover, we analyzed the standard curves of the sensitivity of o-crRNA as well as hs-crRNA1 and hs-crRNA2 using the dose-response regression ( Figure S7, Supporting Information). To investigate the performance of the Cas13a/hs-crRNA duplex in complex sample, we mixed HBV target DNA with five kinds of food-sourced bacteria gDNA and tested the specificity of the hs-CRISPR/Cas13a system. The amounts of bacteria gDNA were all 10 ng. According to the results provided in Figure S8 of the Supporting Information, both the Cas13a/hs-crRNA1 and Cas13a/hs-crRNA2 showed high-specificity in the complex sample, demonstrating the stability of the hs-CRISPR/Cas13a system.
To validate the feasibility of the hs-CRISPR/Cas13a system, we selected another 28-nt sequence containing four random SNPs for HBV B-type and C-type as the target of the o-crRNA-a and hs-crRNA-a. As shown in Figure 3F, hs-crRNA-a was obtained by extending the 3′-end of o-crRNA-a to form a 10-bp stem at  the repeat-proximal region. We examined the specificity of both o-crRNA-a and hs-crRNA-a for HBV genotyping by single-pellet detection and the results are provided in Figure 3G and Figure  S9 (Supporting Information). The on-target signals of o-crRNA-a and hs-crRNA-a were 1561 ± 154.21 RFU, and 1555.86 ± 72.07 RFU, respectively, whereas the off-target signals were 804.67 ± 139.75 and 395.6 ± 59.08 RFU, respectively. The DF values of o-crRNA-a and hs-crRNA-a were calculated and displayed in Table 2. Notably, the hs-crRNA-a increased the discrimination factor of the Cas13a protein by two-fold compared to o-crRNA-a, demonstrating the versatility of hairpin-spacer crRNAs to enhance the specificity of the CRISPR/Cas13a system.

Molecular Docking Simulations of Cas13a and hs-crRNA2
Since hairpin-spacer crRNAs can increase the identification ability of Cas13a system to SNPs, we further performed the molecular docking simulations to illustrate the binding conformation of Cas13a protein with hs-crRNA2 and target RNA. Protein Data Bank ID 5W1H was selected for Cas13a. The structure of hs-crRNA2 is provided in Figure 4A. In the Cas13a/hs-crRNA binary system, the hs-crRNA is located in the groove of the REC lobe (NTD) and the NUC lobe (Helical-2, Helical-3, and HEPN-2) of Cas13a protein, as shown in Figure 4B. As shown in Table S2 of the Supporting Information, the interaction between Cas13a and hs-crRNA2 was mainly electrostatics due to ΔH 0 < 0 and ΔS 0 > 0. The 5′ handle of hs-crRNA localizes on the NTD surface ( Figure 4C), close to channel formed between the NTD and the Helical-1 domains. According to previous reports, [9,23,24] cr-RNA overspreading on the surface of the NTD domain is helpful for subsequent catalysis processes. The 5′ flank is located in the groove of the Helical-1 and the HEPN-2 domains ( Figure 4D). The 3′ flank is oriented to the interface groove of the Helical-2 and the Helical-3 domains, closing to the Helical-3 domain. The 3′ handle of hs-crRNA localizes on the Helical-3 surface (Figure 4E). These domains trapped the hs-crRNA as clamps, and Cas13a protein is accessible for nuclease activity via conformational changes.

Molecular Docking Study of Cas13a/hs-crRNA2 Duplex and Target RNA
In the ternary system, the docking calculations indicated that the recognition of target RNA was guided by the crRNA spacer for Cas13a/hs-crRNA2 duplex (Figure 5A,B). During RNA recognition of the Cas13a/hs-crRNA2 duplex, the primary force was hydrophobic and the secondary forces are hydrogen bond, electrostatic, and base stacking. The spacer of hs-crRNA interacted with target RNA by 37-42, 47-58, and 62-64 NTs, counting from the 5′ end ( Figure 5C). In addition, the hydrogen bond and base stacking were the secondary contributing forces. Besides, target RNA bond to Cas13a via interaction with the Helical-3 domain (Figure 5D). The interaction between target RNA and the Helical-3 domain was governed by electrostatic force. The Helical-3 domain was close to the HEPN-2 domain, and target RNA recognition facilitated the formation of HEPN domain that decided the nuclease activity of Cas13a protein. Our calculations were consistent with previous reports studying the mechanisms of Cas13a activity [9,23,24] . These results suggest that hairpin-spacer crRNAs can be applied to Cas13a system as an effective strategy to promote the specificity.

Conclusion
In this study, we presented a strategy termed hs-CRISPR to increase the specificity of CRISPR/Cas13a system. Three main enhancements were developed and achieved. First, hairpin-spacer crRNAs with engineered secondary structures were proposed to suppress the nonspecific activity of CRISPR/Cas13a system. CRISPR/Cas technology is reported to show high off-target effects which impedes its applications and there lacks a rational guideline for high-specificity Cas13a-crRNA design. Herein, we developed crRNAs with hairpin-spacer and evaluated the binding energetics of Cas13a protein and RNA by fluorescence quenching mechanism. Our enhanced system showed a high on-target affinity of Cas13a while suppressed the off-target activity compared to ordinary crRNA. Second, we established an HBV genotyping platform to evaluate the feasibility of hs-CRISPR to promote the specificity of CRISPR/Cas13a system. Our platform possessed a high sensitivity of 1 × 10 −18 m (100 copies mL −1 ) and hairpinspacer crRNAs increased the discrimination ability of Cas13a system by at least three-fold than ordinary crRNA. The whole assay takes about 1 h and can be encapsulated in a single pellet which is time-saving and easy to operate. Third, we performed theoretical simulations to illustrate the mechanisms of Cas13a nuclease activity by interacting with hairpin-spacer crRNA and target RNA. Our findings suggest that hairpin-spacer crRNAs are compatible with traditional Cas13a system.
In summary, this enhanced hs-CRISPR system was verified to increase the specificity of CRISPR/Cas13a system both theoretically and experimentally in this study. By comparison of ordinary crRNAs and hairpin-spacer crRNAs, the off-target nuclease activity of Cas13a could be suppressed by introducing an extra secondary structure to decrease the binding affinity to off-target RNA. Furthermore, genotyping platforms based on the enhanced hs-CRISPR system possessed several advantages such as fast, easy-operation and programmability. Collectively, our enhanced hs-CRISPR system will facilitate future point-of-care molecular diagnostics and RNA interference based on CRSIPR/Cas13a technology.

Experimental Section
Materials and Instruments: All DNA templates and primers were synthesized from Sangon Biotechnology Co., Ltd. (Shanghai, China). The fluorophore quencher-labeled RNA reporter was purchased from Takara Biotechnology Co., Ltd. (Dalian, China). For nucleic acid targets, plasmids containing HBV genome of genotype B and C were provided by Prof. Y.
M. Wen (Fudan University, Shanghai, China). All sequences and plasmids used in this study were provided in Table S1 of the Supporting Information. The LwaCas13a protein was expressed and obtained from GenScript Co., Ltd. (Jiangsu, China). RPA was performed following instruction of TwistAmp Basic kit from TwistDx Inc. (Cambridge, UK). The dsDNAs were purified using MinElute gel extraction kit from Qiagen (Shanghai, China). crRNAs were purified using RNAXP clean beads from Beckman Coulter (Shanghai, China). The HiScribe T7 Quick High Yield RNA Synthesis kit, murine RNase inhibitor, T7 RNA polymerase, and NTP mix buffer were purchased from New England Biolabs (Beijing, China). MEGAclear Transcription Clean-up kit and RNase-free water were ordered from Thermo Fisher (Shanghai, China).
The real-time fluorescence intensity of nuclease assay was determined by a microplate reader (BioTek Synergy H1). Electrophoresis was carried out using Tanon EPS-100 electrophoresis system. Fluorescence of Cas13a protein was measured using Agilent Technologies Cary Eclipse Fluorescence Spectrometer.
Preparation of crRNAs: The designed crRNA with T7 promoter sequence was used for DNA template and annealed with a short T7 primer (5′-TAATACGACTCACTATAGGG-3′) in 1× annealing buffer containing 100 × 10 −3 m Tris-HCl (pH 7.5), 10 × 10 −3 m EDTA, 1 m NaCl. In vitro transcription of crRNA was performed using the HiScribe T7 Quick High Yield RNA Synthesis kit overnight at 37°C. Transcription reactions contained NTP buffer mix (10 × 10 −3 m each NTP final), T7 RNA polymerase mix, template DNA (2 µg), murine RNase inhibitor (1 U µL −1 final). After the synthesis of crRNA, adequate RNase-free water, and DNase I (4 units final) were added into the system to digest the DNA template by incubating at 37°C for 15 min. The transcribed crRNAs were then purified with twofold volume of RNAXP clean beads and subsequently stored at −20°C. The obtained crRNAs were measured by Nanodrop and analyzed by denaturing urea polyacrylamide gel electrophoresis (urea/PAGE).
HBV RNA Preparation and Purification: First, HBV plasmids were amplified by RPA using TwistAmp Basic kit. RPA primers for HBV targets were designed using NCBI Primer -Blast. Notably, the HBV RPA forward primers should contain a short T7 sequence at its 5′-end for later transcription. Second, the amplified products were evaluated by agarose gel electrophoresis and then purified by MinElute PCR purification kit (Qiagen, China). Third, the purified HBV DNA amplicons contained extended T7 sequences were transcribed into RNA using the HiScribe T7 Quick High Yield RNA Synthesis kit overnight at 37°C. Finally, the transcribed HBV RNA were purified by the MEGAclear Transcription Clean-up kit and then evaluated by urea/PAGE.
Fluorescence Measurement: LwaCas13a proteins have natural fluorophores (tyrosine and tryptophan residues) and allow the study of conformational changes. Intrinsic fluorescence spectra of Cas13a protein (5 × 10 −6 m) incubated with various concentrations of crRNAs (0 × 10 −6 to 2.5 × 10 −6 m) were measured using Agilent Technologies Cary Eclipse Fluorescence Spectrometer in cuvette. Afterward, the Cas13a/crRNA duplex was incubated with transcribed HBV RNA (concentrations ranging from 0 × 10 −6 to 1 × 10 −6 m) and emission spectra were recorded. The excitation and emission wavelengths were set at 285 and 295-500 nm, respectively. Experiments were performed at 25 and 37°C for all RNAs.
Molecular Docking Study of Cas13a/crRNA and Cas13a/crRNA/Target RNA Complex: The molecular docking simulation was performed to predict the consensus pose corresponds to the interaction interface at the primary state in binary and ternary systems of Cas13a/crRNA and Cas13a/crRNA-target RNA, respectively. The docking calculations were processed by the HDOCK server (http://hdock.phys.hust.edu.cn/). [40,41] The 3D structure of Cas13a protein was obtained from the Protein Data Bank ID 5W1H. The 3D structures of crRNA and target RNA were generated by the RNAComposer server (http://rnacomposer.cs.put.poznan. pl/). [42] Rigid docking calculations were carried out to investigate the binding activity between Cas13a and crRNA, followed for the Cas13a/crRNA duplex and target RNA. Besides, possible interactions were evaluated at the primary state in binary and ternary systems. The docking results were analyzed using Visual Molecular Dynamics. [43] The positions of molecules were investigated based on the energy score. The best pose of interaction was selected based on the minimum energy HDOCK score.
Statistics Analysis: Statistical analysis was performed by two-tailed Student's t-test using Graphpad Prism 8.0 and Origin 2021. Each experiment was repeated three times. The results were presented as mean ± S.D.

Supporting Information
Supporting Information is available from the Wiley Online Library or from the author.