The Optimization and Characterization of an RNA-Cleaving Fluorogenic DNAzyme Probe for MDA-MB-231 Cell Detection

Breast cancer is one of the most frequently diagnosed cancers in females worldwide and lacks specific biomarkers for early detection. In a previous study, we obtained a selective RNA-cleaving Fluorogenic DNAzyme (RFD) probe against MDA-MB-231 cells, typical breast cancer cells, through the systematic evolution of ligands by exponential process (SELEX). To improve the performance of this probe for actual application, we carried out a series of optimization experiments on the pH value of a reaction buffer, the type and concentration of cofactor ions, and sequence minimization. The length of the active domain of the probe reduced to 25 nt from 40 nt after optimization, which was synthesized more easily and economically. The detection limit of the optimized assay system was 2000 MDA-MB-231 cells in 30 min, which is more sensitive than the previous one (almost 5000 cells). The DNAzyme probe was also capable of distinguishing MDA-MB-231 cell specifically from 3 normal cells and 10 other tumor cells. This probe with high sensitivity, selectivity, and economic efficiency enhances the feasibility for further clinical application in breast cancer diagnosis. Herein, we developed an optimization system to produce a general strategy to establish an easy-to-use DNAzyme-based assay for other targets.


Introduction
Breast cancer, one of the most frequently diagnosed cancers, is the second leading cause of death by cancer in women [1,2]. The American Cancer Society (http://www.cancer.org) estimates that about 252,710 new cases of invasive breast cancer will be diagnosed in women, and about 40,610 women will die from breast cancer in the United States in 2017. In addition, the World Health Organization (WHO) reports an estimated 17 million deaths from breast cancer by 2030 [3]. Therefore, breast cancer remains a tremendous threat to the human health system. Clinical studies have demonstrated that the survival rate of patients is closely related to the time of disease diagnosis [4].
Early detection and timely therapy of breast cancer has long been considered to be a vital strategy in the control of the disease's development, which can lead to timely care of patients and prevent potential tumor metastasis [5,6]. Considering the importance of early detection, many conventional techniques, such as touch detection, imaging methods, biopsy, mammography, and serum tumor chemicals were of analytical reagent grade and used without further purification. The water used in this study was double deionized (ddH 2 O) and autoclaved with a Milli-Q Synthesis from a Millipore system. All DNA oligonucleotides were synthesized by Takara Biotechnology Co. Ltd. (Dalian, China). All DNA oligonucleotides were purified by 10% denatured PAGE before use except that the FAM-labeled DNA was purified by HPLC. The selection buffer (2×) contained 400 mM NaCl, 100 mM HEPES (pH 7.5), and 10 mM MgCl 2 . Lysis buffer contained 10 mM HEPES (pH 7.9), 10 mM KCl, 1 mM EDTA, 0.1% NP40, and 1× complete Protease Inhibitor Cocktail. All solutions were freshly prepared before use.

Sample Preparation
All kinds of cells were washed twice with 1× PBS and scraped into cell lysis buffer. The cell resuspensions were sonicated for 4 s and suspended for 9 s, this operation was repeated three times. The suspensions were incubated on ice for 30 min, while being inverted every 5 min. Finally, the cell lysate was centrifuged for 15 min (20,000 g, 4 • C). The supernatant was collected and stored in −80 • C for further use.

Sequences Ligation Process
Standard and modified DNA oligonucleotides were synthesized by Takara Biotechnology Co. Ltd. (Dalian, China). The ligation reaction was performed as follows and the DNA sequences used in the ligation process are summarized in Table 1. Substrate (100 pmol) was first phosphorylated with a small amount of ATP (100 nmol) and PNK (0.2 U/µL, final concentration, the total reaction volume is 100 µL) for 30 min at 37 • C. Then the mixture was heated to 90 • C for 5 min in order to inactivate the PNK. Next, this mixture was directly added into a solution containing 120 pmol AL (the sequence obtained through SELEX and not ligated with substrate) and 120 pmol DNA template. The whole system was heated at 90 • C for 50 s to completely inactivate the PNK and despiralize the DNA sequence and cooled to room temperature. The ligation was initiated by the addition of T4 DNA ligase (0.2 U/µL, final concentration and final volume was 200 µL) and continued at room temperature for 90 min. The ligated DNA was subsequently purified by 10% denatured polyacrylamide gel electrophoresis (dPAGE), followed by ethanol precipitation (repeated three times). The purpose of the operation was to separate the ligated DNA from unligated DNA. It should be noted that the substrate contained a fluorescein-labeled thymine deoxyribonucleotide (Fluorescein-dT, F), a adenine ribonucleotide (R) as the cleavage site and a DABCYL-labeled thymine deoxyribonucleotide (DABCYL-dT, Q).

Results and Discussion
According to a previous study [29], the AAI2-5 probe was modified with fluorescein and DABCYL to show green fluorescence excitation at 485 nm and to induce fluorescence quenching, respectively. Under natural conditions, fluorescein-dT became closer to DABCYL-dT, leading to quenched fluorescence. When the target was introduced, the catalytic ability of the probe was activated, followed by the cleaving of the substrate, increasing the distance between DABCYL-dT and fluorescein-dT and increasing the fluorescence signal (Scheme 1). It was found that the cleavage of the substrates could generate a nearly six-fold increase in fluorescence in the PAGE-based experiment [27,29,34,[38][39][40][41]. We calculated the percentage of cleavage using the formula [ thymine deoxyribonucleotide (Fluorescein-dT, F), a adenine ribonucleotide (R) as the cleavage site and a DABCYL-labeled thymine deoxyribonucleotide (DABCYL-dT, Q).

Results and Discussion
According to a previous study [29], the AAI2-5 probe was modified with fluorescein and DABCYL to show green fluorescence excitation at 485 nm and to induce fluorescence quenching, respectively. Under natural conditions, fluorescein-dT became closer to DABCYL-dT, leading to quenched fluorescence. When the target was introduced, the catalytic ability of the probe was activated, followed by the cleaving of the substrate, increasing the distance between DABCYL-dT and fluorescein-dT and increasing the fluorescence signal (Scheme 1). It was found that the cleavage of the substrates could generate a nearly six-fold increase in fluorescence in the PAGE-based experiment [27,29,34,[38][39][40][41]. We calculated the percentage of cleavage using the formula [

Optimization of the Reaction Buffer Conditions
The reaction condition was mainly affected by the pH value of the reaction buffer as well as the type and concentration of cofactor ions. We optimized these conditions one by one. For all cleavage reactions conducted in this part, the reaction system contains 1× selection buffer, 0.5 pmol/μL RFD probe, 0.2 mM EGTA, and 30 μg/mL cell lysate. The mixture was incubated for 60 min at room temperature. The cleavage product was separated by 10% dPAGE. As shown in Figure 1a, the catalytic activity only occurred in the pH range from 7 to 10, which is an alkaline condition. Because more OH − are present in alkaline conditions, the inherent transesterification reaction may occur more easily between a phosphodiester linkage and the nearby 2-hydroxyl group [28]. Therefore, the catalytic ability of the DNAzyme was stronger under the selected alkaline conditions [42]. However, we did not test a reaction in pH greater than 10 considering that the probe was cleaved in higher pH conditions without the target [28]. According to the data, the maximum catalytic rate occurred at pH 8, which was similar to the bio-system and was chosen as the best pH condition in subsequent experiments.
In the RFD catalytic reaction, bivalent ions play an important and necessary role as metal cofactors [43][44][45][46][47][48][49]. In this work, we studied eight other common bivalent ions (Ba 2+ , Ca 2+ , Cu 2+ , Fe 2+ , Ni 2+ , Zn 2+ , Mn 2+ , Hg 2+ ) as replacements of Mg 2+ (10 mM) in the selection buffer. Generally speaking, Mg 2+ (or another divalent metal cation) as a cofactor can accelerate the cleavage of a particular RNA phosphoester [23]. The divalent metal ions bind to specific sites and form a stable and active structure. Interestingly, in our experiments, all eight metal cations scarcely stimulated the enhancement of fluorescence signal (Figure 1b). A possible explanation is that Mg 2+ has unique physical and chemical properties, such as metal coordination number, ionic radius, electrophilicity, oxidation Scheme 1. The working principle of the probe.

Optimization of the Reaction Buffer Conditions
The reaction condition was mainly affected by the pH value of the reaction buffer as well as the type and concentration of cofactor ions. We optimized these conditions one by one. For all cleavage reactions conducted in this part, the reaction system contains 1× selection buffer, 0.5 pmol/µL RFD probe, 0.2 mM EGTA, and 30 µg/mL cell lysate. The mixture was incubated for 60 min at room temperature. The cleavage product was separated by 10% dPAGE. As shown in Figure 1a, the catalytic activity only occurred in the pH range from 7 to 10, which is an alkaline condition. Because more OH − are present in alkaline conditions, the inherent transesterification reaction may occur more easily between a phosphodiester linkage and the nearby 2-hydroxyl group [28]. Therefore, the catalytic ability of the DNAzyme was stronger under the selected alkaline conditions [42]. However, we did not test a reaction in pH greater than 10 considering that the probe was cleaved in higher pH conditions without the target [28]. According to the data, the maximum catalytic rate occurred at pH 8, which was similar to the bio-system and was chosen as the best pH condition in subsequent experiments.
In the RFD catalytic reaction, bivalent ions play an important and necessary role as metal cofactors [43][44][45][46][47][48][49]. In this work, we studied eight other common bivalent ions (Ba 2+ , Ca 2+ , Cu 2+ , Fe 2+ , Ni 2+ , Zn 2+ , Mn 2+ , Hg 2+ ) as replacements of Mg 2+ (10 mM) in the selection buffer. Generally speaking, Mg 2+ (or another divalent metal cation) as a cofactor can accelerate the cleavage of a particular RNA phosphoester [23]. The divalent metal ions bind to specific sites and form a stable and active structure. Interestingly, in our experiments, all eight metal cations scarcely stimulated the enhancement of fluorescence signal (Figure 1b). A possible explanation is that Mg 2+ has unique physical and chemical properties, such as metal coordination number, ionic radius, electrophilicity, oxidation state, and natural affinity with RNA [42,43]. Therefore, we selected Mg 2+ as the best cofactor for subsequent experiments. state, and natural affinity with RNA [42,43]. Therefore, we selected Mg 2+ as the best cofactor for subsequent experiments. We then studied the effect of the concentration of Mg 2+ on the activity of DNAzyme. We incubated different concentrations of Mg 2+ , from 0.1 mM to 10 mM, in the reaction system for 60 min at room temperature. The results are shown in Figure 1c. Considering that the rate of cleavage is largely dependent on the concentration of Mg 2+ , the cleavage efficiency was enhanced when the concentration of Mg 2+ increased and a stable cleaving value appeared at a concentration of 4 mM. The probe requires Mg 2+ as a cofactor to catalyze the RNA; thus, more Mg 2+ could accelerate the cleavage of the particular RNA phosphoester. Interestingly, we found that 4 mM is also the saturated concentration of Mg 2+ in the reaction, meaning that a higher concentration is not necessary. A lower concentration of the cofactor reduced the ionic effect in the nonspecific interaction, which is more applicable for practical implementation of the probe.

Sequence Minimization
In the original selection experiment, we established a library of 40-nt single-stranded DNA molecules, which acted as the active domain to recognize the target and cleave the substrate. The main reason for using a relatively long DNA sequence was that a longer DNA sequence may result in better catalytic efficiency [50][51][52]. Although DNA with long sequences have been isolated, the crucial domain pertains to only a shorter piece of the sequence, such as 8-17 [25] and 10-23 [25]. Therefore, we explored the active domain of the RFD probe, AAI2-5 (A0), for its ability to recognize MDA-MB-231 cells with relatively limited nucleotides via sequence minimization.
The obtained sequence of A0 in the previous work is shown in Figure 2. For the convenience of description, we divided the original A0 into four parts including P1, P2, catalytic domain and substrate. The parts of P1 and P2 were used as primers to amplify sequences in the selection process. The sequence signed in italic was catalytic domain which binds with target and cleaves the substrate. And the sequence modified with F, R, Q was substrate which could generate the fluorescence signal. Then, the catalytic domain was further divided into eight units with 5-nt each. In this strategy, the catalytic ability of the probes was scored by the cleavage efficiency after the reaction: a cleavage percentage less than 10% with a minus sign (−); 10%-20% with one plus sign (+); 20%-30% with two plus signs (++); and cleavage activity greater than 30% indicated that the deleted sequence had little We then studied the effect of the concentration of Mg 2+ on the activity of DNAzyme. We incubated different concentrations of Mg 2+ , from 0.1 mM to 10 mM, in the reaction system for 60 min at room temperature. The results are shown in Figure 1c. Considering that the rate of cleavage is largely dependent on the concentration of Mg 2+ , the cleavage efficiency was enhanced when the concentration of Mg 2+ increased and a stable cleaving value appeared at a concentration of 4 mM. The probe requires Mg 2+ as a cofactor to catalyze the RNA; thus, more Mg 2+ could accelerate the cleavage of the particular RNA phosphoester. Interestingly, we found that 4 mM is also the saturated concentration of Mg 2+ in the reaction, meaning that a higher concentration is not necessary. A lower concentration of the cofactor reduced the ionic effect in the nonspecific interaction, which is more applicable for practical implementation of the probe.

Sequence Minimization
In the original selection experiment, we established a library of 40-nt single-stranded DNA molecules, which acted as the active domain to recognize the target and cleave the substrate. The main reason for using a relatively long DNA sequence was that a longer DNA sequence may result in better catalytic efficiency [50][51][52]. Although DNA with long sequences have been isolated, the crucial domain pertains to only a shorter piece of the sequence, such as 8-17 [25] and 10-23 [25]. Therefore, we explored the active domain of the RFD probe, AAI2-5 (A0), for its ability to recognize MDA-MB-231 cells with relatively limited nucleotides via sequence minimization.
The obtained sequence of A0 in the previous work is shown in Figure 2. For the convenience of description, we divided the original A0 into four parts including P1, P2, catalytic domain and substrate. The parts of P1 and P2 were used as primers to amplify sequences in the selection process. The sequence signed in italic was catalytic domain which binds with target and cleaves the substrate. And the sequence modified with F, R, Q was substrate which could generate the fluorescence signal. Then, the catalytic domain was further divided into eight units with 5-nt each. In this strategy, the catalytic ability of the probes was scored by the cleavage efficiency after the reaction: a cleavage percentage less than 10% with a minus sign (−); 10%-20% with one plus sign (+); 20%-30% with two plus signs (++); and cleavage activity greater than 30% indicated that the deleted sequence had little effect on the probe activity and was represented by three plus signs (+++). All truncation sequences were examined using A0 as the reference deoxyribozyme. effect on the probe activity and was represented by three plus signs (+++). All truncation sequences were examined using A0 as the reference deoxyribozyme.
In the truncation design, the substrate domain remained constant, which ensured that the signal would be detected and other parts would be truncated to varying degrees. To first explore the function of the primes in the catalytic reaction, we truncated the P1 part-generating AP1 sequence and P2 part-generating AP2 sequence, respectively as: Figure 2. The original sequence of AAI2-5 (A0) which was divided into four parts including primer sequences P1 and P2, substrate, and active domain (italic letters) which was further divided into eight units. Figure 3 shows that only the original sequence, A0, exhibited a good cleavage percentage, while AP1 and AP2 presented apparent activity reduction with cleavage percent less than 10%. These results indicated that the two primer sequences are essential for maintaining the conformation of this RFD probe. In the next experiments, we reserved these two parts and then optimized the active domain using nucleotide truncation. We sequentially removed any 5-nt unit from the 5′ to 3′ end of the active domain and named the sequences A1-A8 ( Figure 4).   In the truncation design, the substrate domain remained constant, which ensured that the signal would be detected and other parts would be truncated to varying degrees. To first explore the function of the primes in the catalytic reaction, we truncated the P1 part-generating AP1 sequence and P2 part-generating AP2 sequence, respectively as: Figure 3 shows that only the original sequence, A0, exhibited a good cleavage percentage, while AP1 and AP2 presented apparent activity reduction with cleavage percent less than 10%. These results indicated that the two primer sequences are essential for maintaining the conformation of this RFD probe. In the next experiments, we reserved these two parts and then optimized the active domain using nucleotide truncation. We sequentially removed any 5-nt unit from the 5 to 3 end of the active domain and named the sequences A1-A8 (Figure 4). effect on the probe activity and was represented by three plus signs (+++). All truncation sequences were examined using A0 as the reference deoxyribozyme.
In the truncation design, the substrate domain remained constant, which ensured that the signal would be detected and other parts would be truncated to varying degrees. To first explore the function of the primes in the catalytic reaction, we truncated the P1 part-generating AP1 sequence and P2 part-generating AP2 sequence, respectively as: Figure 2. The original sequence of AAI2-5 (A0) which was divided into four parts including primer sequences P1 and P2, substrate, and active domain (italic letters) which was further divided into eight units. Figure 3 shows that only the original sequence, A0, exhibited a good cleavage percentage, while AP1 and AP2 presented apparent activity reduction with cleavage percent less than 10%. These results indicated that the two primer sequences are essential for maintaining the conformation of this RFD probe. In the next experiments, we reserved these two parts and then optimized the active domain using nucleotide truncation. We sequentially removed any 5-nt unit from the 5′ to 3′ end of the active domain and named the sequences A1-A8 (Figure 4).   effect on the probe activity and was represented by three plus signs (+++). All truncation sequences were examined using A0 as the reference deoxyribozyme.
In the truncation design, the substrate domain remained constant, which ensured that the signal would be detected and other parts would be truncated to varying degrees. To first explore the function of the primes in the catalytic reaction, we truncated the P1 part-generating AP1 sequence and P2 part-generating AP2 sequence, respectively as: Figure 2. The original sequence of AAI2-5 (A0) which was divided into four parts including primer sequences P1 and P2, substrate, and active domain (italic letters) which was further divided into eight units. Figure 3 shows that only the original sequence, A0, exhibited a good cleavage percentage, while AP1 and AP2 presented apparent activity reduction with cleavage percent less than 10%. These results indicated that the two primer sequences are essential for maintaining the conformation of this RFD probe. In the next experiments, we reserved these two parts and then optimized the active domain using nucleotide truncation. We sequentially removed any 5-nt unit from the 5′ to 3′ end of the active domain and named the sequences A1-A8 (Figure 4).   Interestingly, we found that all probes could cleave the substrate efficiently except for deletion unit 5 (A5, -) and unit 6 (A6, -) sequences. A5 and A6 presented poor catalytic activity, meaning that these 10 nucleotides are required for the catalytic function of deoxyribozyme, may engage in secondary interactions, and are responsible for substrate recognition. Therefore, we reserved these two units in the next minimization experiments. After the first sequence deletion round, we designed further truncation sequences with 10-nt and 15-nt deletion with constant units 5 and 6. Another six sequences named A9-A14 were obtained, which were found to exhibit medium (A9, A10, A11, A12, and A14) or good activity (A13). And the results were shown in Figure 5. It is important to mention that the A13 sequence showed a higher cleavage activity than A0. Interestingly, we found that all probes could cleave the substrate efficiently except for deletion unit 5 (A5, -) and unit 6 (A6, -) sequences. A5 and A6 presented poor catalytic activity, meaning that these 10 nucleotides are required for the catalytic function of deoxyribozyme, may engage in secondary interactions, and are responsible for substrate recognition. Therefore, we reserved these two units in the next minimization experiments. After the first sequence deletion round, we designed further truncation sequences with 10-nt and 15-nt deletion with constant units 5 and 6. Another six sequences named A9-A14 were obtained, which were found to exhibit medium (A9, A10, A11, A12, and A14) or good activity (A13). And the results were shown in Figure 5. It is important to mention that the A13 sequence showed a higher cleavage activity than A0. These results suggest that unit 1, unit 2, and unit 8 were dispensable in the catalytic reaction, while units 3-7 played an important role in maintaining the active structure. In this work, the active domain of the probe was reduced to 25 nt from 40 nt after optimization, which was synthesized more easily and economically.

Sensitivity of the Optimized Probe
According to the results, the original probe, A0, detected 0.5 μg/mL (nearly 5000 cell/mL) cell lysate proteins of MDA-MB-231 cells [29]. In this work, we tested the detection limit of A0 and A13 ( Figure 6). 0.5 pmol/μL probe A0 and A13 were incubated with 1× selection buffer, 0.2 mM EGTA, and different concentrations of MDA-MB-231 cell lysate, separately. To describe more conveniently, we used cell numbers to represent the concentration of the target. A0 showed 12.8% and 9.4% cleavage when the cell concentration was 6000 cell/mL and 4000 cell/mL, respectively. Compared to A0, A13 with 9.8% cleavage presented a better sensitivity when the cell concentration was 2000 cell/mL. In other words, through a series of optimization experiments, we obtained a more sensitive RFD probe for MDA-MB-231 cells. Probe A13 with a shorter sequence may fold a simpler, more preferable, and effective structure, while deleting some dispensable bases can also reduce the steric hindrance effect. The results suggest that our optimization enhanced the sensitivity of the assay and the feasibility of the probe's application. These results suggest that unit 1, unit 2, and unit 8 were dispensable in the catalytic reaction, while units 3-7 played an important role in maintaining the active structure. In this work, the active domain of the probe was reduced to 25 nt from 40 nt after optimization, which was synthesized more easily and economically.

Sensitivity of the Optimized Probe
According to the results, the original probe, A0, detected 0.5 µg/mL (nearly 5000 cell/mL) cell lysate proteins of MDA-MB-231 cells [29]. In this work, we tested the detection limit of A0 and A13 ( Figure 6). 0.5 pmol/µL probe A0 and A13 were incubated with 1× selection buffer, 0.2 mM EGTA, and different concentrations of MDA-MB-231 cell lysate, separately. To describe more conveniently, we used cell numbers to represent the concentration of the target. A0 showed 12.8% and 9.4% cleavage when the cell concentration was 6000 cell/mL and 4000 cell/mL, respectively. Compared to A0, A13 with 9.8% cleavage presented a better sensitivity when the cell concentration was 2000 cell/mL. In other words, through a series of optimization experiments, we obtained a more sensitive RFD probe for MDA-MB-231 cells. Probe A13 with a shorter sequence may fold a simpler, more preferable, and effective structure, while deleting some dispensable bases can also reduce the steric hindrance effect. The results suggest that our optimization enhanced the sensitivity of the assay and the feasibility of the probe's application.

Specificity of the New Probe
To investigate the specificity of the optimized probe, we subsequently tested the response of A13 to various cell lysates produced by different cell lines, including 3 normal cell lines (QSG 7701, HEB, MCF-10A) and 10 tumor cell lines (BT474, SK-BR-3, MCF-7, ZR-75-1, K562, HCT116, NCI-H69, CCRF-CEM, Hela, and U251), which contain breast cancer, liver cancer, small cell lung cancer, etc. All cell lysates were quantified using a BCA protein assay kit. In this experiment, 0.5 pmol/μL probe A13 was incubated with 60 μg/mL MDA-MB-231 cell lysate, 1× selection buffer, and 0.2 mM EGTA for 30 min and displayed significant cleavage efficiency, as shown in Figure 7. This probe is capable of specifically distinguishing MDA-MB-231 cells from normal cells and other tumor cells. These results strongly prove that the optimized probe A13 has excellent specificity like A0.

Specificity of the New Probe
To investigate the specificity of the optimized probe, we subsequently tested the response of A13 to various cell lysates produced by different cell lines, including 3 normal cell lines (QSG 7701, HEB, MCF-10A) and 10 tumor cell lines (BT474, SK-BR-3, MCF-7, ZR-75-1, K562, HCT116, NCI-H69, CCRF-CEM, Hela, and U251), which contain breast cancer, liver cancer, small cell lung cancer, etc. All cell lysates were quantified using a BCA protein assay kit. In this experiment, 0.5 pmol/µL probe A13 was incubated with 60 µg/mL MDA-MB-231 cell lysate, 1× selection buffer, and 0.2 mM EGTA for 30 min and displayed significant cleavage efficiency, as shown in Figure 7. This probe is capable of specifically distinguishing MDA-MB-231 cells from normal cells and other tumor cells. These results strongly prove that the optimized probe A13 has excellent specificity like A0.

Conclusions
In this work, we successfully optimized and characterized the RFD probe for MDA-MB-231 cells through series experiments, including reaction condition optimization and sequence minimization. We reduced the cofactor concentration of Mg 2+ to 4 mM, which weakened the ionic effect in the nonspecific interaction and is useful for practical application of the probe. After sequence truncation, we obtained a better probe, A13, which is 16-nt shorter than the initial probe, A0. Probe A13 presented excellent specificity and could distinguish MDA-MB-231 cells from 10 other tumor cell lines and 3 normal human cell lines. More importantly, A13 could detect as low as 2000 MDA-MB-231 tumor cells, which is more sensitive than probe A0. In summary, the shorter length of the probe displayed good sensitivity, excellent specificity, and more economical efficiency. These properties enhance the feasibility of the probe for further clinical application in breast cancer diagnosis. Based on the results, our optimization system may be efficiently used to develop a general strategy to establish an easyto-use DNAzyme-based assay for other targets.

Conclusions
In this work, we successfully optimized and characterized the RFD probe for MDA-MB-231 cells through series experiments, including reaction condition optimization and sequence minimization. We reduced the cofactor concentration of Mg 2+ to 4 mM, which weakened the ionic effect in the nonspecific interaction and is useful for practical application of the probe. After sequence truncation, we obtained a better probe, A13, which is 16-nt shorter than the initial probe, A0. Probe A13 presented excellent specificity and could distinguish MDA-MB-231 cells from 10 other tumor cell lines and 3 normal human cell lines. More importantly, A13 could detect as low as 2000 MDA-MB-231 tumor cells, which is more sensitive than probe A0. In summary, the shorter length of the probe displayed good sensitivity, excellent specificity, and more economical efficiency. These properties enhance the feasibility of the probe for further clinical application in breast cancer diagnosis. Based on the results, our optimization system may be efficiently used to develop a general strategy to establish an easy-to-use DNAzyme-based assay for other targets.