The Magnetic Bead Computing Model of the 0-1 Integer Programming Problem Based on DNA Cycle Hybridization

Magnetic beads and magnetic Raman technology substrates have good magnetic response ability and surface-enhanced Raman technology (SERS) activity.-erefore, magnetic beads exhibit high sensitivity in SERS detection. In this paper, DNA cycle hybridization and magnetic bead models are combined to solve 0-1 integer programming problems. First, the model maps the variables to DNA strands with hairpin structures and weights them by the number of hairpin DNA strands. -is result can be displayed by the specific binding of streptavidin and biotin. Second, the constraint condition of the 0-1 integer programming problem can be accomplished by detecting the signal intensity of the biological barcode to find the optimal solution. Finally, thismodel can be used to solve the general 0-1 integer programming problem and has more extensive applications than the previous DNA computing model.


Introduction
With the development of science and technology, traditional computing has been unable to meet people's requirements when dealing with massive data and information processing, and people have started to explore new fields of computing. Since Adleman proposed the use of DNA computing to solve the directed Hamilton path in 1994, DNA computing has received increasing attention from researchers [1]. In 2000, Head et al. proposed a new method of computing by using DNA plasmids and reported the NP-complete problem concerning the cardinality of the largest independent subset of the vertex set of the computing graph [2]. In 2011, Zhang et al. designed a DNA word set based on minimum free energy [3]. In 2017, Yin and Cui reported the integer programming problem based on the plasmid DNA computing model [4]. In 2018, Ramanamurthy introduced the basic structure of DNA and DNA processing tools [5]. In 2019, Tang established a dynamic NAND computing model using DNA origami [6]. In the same year, Yang et al. used DNA origami and hybridization chain reaction to solve a new computational model for solving the knapsack problem [7]. DNA cycle hybridization chain reaction is a process of alternating hybridization of two DNA molecules with different hairpin structures induced under the induction of a trigger strand.
is process is spontaneous and does not require the involvement of enzymes. With the development of science and technology, DNA cycle hybridization chain reaction has been applied to many fields, such as biosensing, biomedicine, proteins, and others. In 2004, Dirks first proposed the concept and indicated that DNA can be used as an amplified transducer for biosensing applications [8]. In 2016, Guo proposed a new chemical immunoassay method for signal amplification that can detect multiple tumor biomarkers simultaneously [9]. In the same year, Yang designed an Aptamer-Binding Directed DNA Origami Pattern for logic gates [10]. In 2018, Li et al. proposed a method for label-free lighting of fluorescent sensors using hybrid chain reaction and DNA triple-strand assembly [11]. In 2018, Xiao designed multiple chemiluminescence imaging and used it for sensitive screening and detection of protein biomarkers through the use of DNA microarray and hybridization chain reaction amplification integration induced by adjacent binding [12]. 0-1 integer programming is a special type of integer programming problem, and its variable values are only 0 and 1. It is widely used in a variety of problems, such as line design, backpack problems, and dispatch problems [13][14][15][16][17]. In 2006, Yin et al. designed a molecular beacon model for solving integer linear programming problems [13]. In 2010, Huang used the advantage of DNA tiles to build a molecular computing system to solve the 0-1 programming problem [14]. In 2017, Li devised a new DNA computing method to solve the 0-1 programming problem. e method is based on a self-assembled nanoparticle probe, which reduces the cost of the model and improves the sensitivity and accuracy of detection [17].
In this paper, based on DNA cycle hybridization chain reaction, a magnetic bead model is constructed to solve the general 0-1 integer programming problem. is magnetic bead model maps variables into the hairpin structure DNA strands and maps weights according to the number of hairpin structure DNA strands. is structure can be displayed according to the specific combination of streptavidin and biotin, and SERS detection has been used [18][19][20]. e general structure of this paper is as follows. First, the basic principle of the DNA cycle hybridization chain reaction and the general form of the 0-1 integer programming problem are given. Second, the magnetic bead model is constructed, and the algorithm steps of the model and a concrete example are given. Finally, this paper uses Visual DSD software to simulate and analyze the optimal solution of the 0-1 integer programming problem and then provides the conclusion. e resulting hairpin structure H 1 is opened, exposing the gap n − m * .
is exposed gap region n − m * happens to have complementary base pairing with n * − m from the H 2 region. Opening the structure of H 2 to expose the area l * − m * , the exposed area l * − m * from H 2 will continue to pair with the next base complement of H 1 . Opening the structure of H 1 to expose the area n − m * , this exposed gap region will continue to have complementary base pairing with n * − m from the H 2 region. Opening the structure of H 2 , H 1 and H 2 are turned on alternately in turn until the sum in the solution is used up. A double-stranded DNA product that hybridizes H 1 and alternately H 2 with a gap is generated.
e most commonly used method for solving the 0-1 integer programming problem is the exhaustive method (also known as the forced search method), which traverses the entire search space, but this method is time-and laborconsuming.
In this paper, the general 0-1 integer programming problem is solved based on the DNA cyclc hybridization chain reaction, which is a generalization of the assignment problem.
e algorithm steps of the general 0-1 integer programming problem are discussed as follows: Step 1. All possible solutions to the problem with variables equal to 0 or 1 are given Step 2. According to the Raman signal intensity, the nonfeasible solutions in the first constraint condition are eliminated, and the feasible solutions are retained Step 3. Repeat step 2 with the reserved feasible solutions; then, all nonfeasible solutions can be eliminated, and all feasible solutions of the problem can be obtained Step 4. e corresponding objective function values of each feasible solution are compared to obtain the optimal solution of the objective function.

Magnetic Bead Model of the 0-1 Integer
Programming Problem Based on DNA Cycle Hybridization

Building a Magnetic Bead Computing Model.
In this paper, we study a new magnetic bead computing model for 0-1 integer programming problems. e magnetic bead computing model effectively utilizes the specific binding effect of streptavidin and biotin through DNA cycle hybridization technology and Raman technology to detect the signal released by the biological barcode. Table 1 shows the coding sequences of the three DNA molecules. Figure 2 shows the DNA cycle hybridization process.
As shown in Figure 2, at the optimal experimental conditions of 37°C and a pH of 7.4, the concentration of the captured stranded DNA was 1.0 × 10 − 7 M [21]. e single strand of capture DNA was fixed on the magnetic bead by means of the amide bond between -COOH modified by the magnetic bead and -NH2 on the DNA strand. e single strand of capture DNA fixed on the magnetic bead was complementarily paired with the base at the sticky end of the strand of the hairpin DNA1 strand, thus opening the hairpin structure DNA1 strand. When the hairpin structure DNA1 is opened, it is complementary with the sticky end base of the hairpin structure DNA2, thus opening the hairpin structure of DNA2. After the hairpin structure DNA2 is opened, it continues to be complementary to the sticky end of the hairpin structure of DNA1. In this way, the hairpin structure DNA1 strand and DNA2 strand cycle hybridizes successively, forming "magnetic bead-capture DNA-DNA1-DNA2-DNA1-...-DNA2-DNA1-DNA2," a special doublestranded DNA molecule. Until DNA1 and DNA2 are consumed in the solution, the sticky end of DNA1 and the sticky end of DNA2 are both modified by biotin and bind to the strepavidin-modified nanobiotic barcode specifically to achieve signal release.
In summary, for the 0-1 integer programming problem with n variables (x 1 , x 2 , . . . , x n ) and m constraint equations, the specific algorithm of the general 0-1 integer programming computing model is as follows: Step 1. First, for n variables in each constraint condition, n magnetic beads with capture DNA were designed (magnetic beads with different radii represent different variables). Second, two types of hairpin structure DNA strands were designed, which were named DNA1 and DNA2. e capture DNA fixed on the magnetic bead can open the hairpin structure DNA1. A gap appeared after the hairpin structure DNA1 was opened, which could be further opened to design the hairpin structure DNA2. In this way, DNA1, DNA2, DNA1, DNA2... cyclically cross each other in turn. Until DNA1 and DNA2 in the solution are consumed (the specific coding sequence design of DNA1 and DNA2 is shown in Table 1), when x i � 1, the sticky ends of the DNA1 strand and the sticky ends of the DNA2 strand were modified with biotin; when x i � 0, the sticky ends of the DNA1 strand and the sticky ends of the DNA2 strand did not need to be modified with biotin, as shown in Figure 3.
Step 2. First, a proper number of biological barcodes were placed in the data pool and mixed evenly. Second, a set of test tubes was prepared for each constraint condition, and each set of test tubes had 2 n test tubes (where n represents the number of variables in the constraint condition). Finally, equal amounts of the solution were placed in the desired tubes.
Step 3. For the first constraint, according to the number of possible solutions k, take out the test tubes according to step 2 and group them and place a magnetic bead with capture DNA in each test tube of each group. After that, put equal amounts of DNA1 and DNA2 into the solution according to the weight coefficients of the variables in the constraint condition. at is, the total amounts of DNA1 and DNA2 are the same as the weight coefficients of the variables. At the same time, according to the characteristics of DNA cycle hybridization, we alternately put DNA1 and DNA2 into the solution every time and put the DNA1 strand first.
Step 4. When the biological barcode in the solution is combined with the sticky ends of the biotin of DNA1 and DNA2, the cycle hybridization signal will be amplified, and the feasibility solution will be judged by the intensity of amplification of the cycle hybridization signal. Here, it is stated that, when there is no biotin at the sticky ends of DNA1 and DNA2, the signal intensity is 0, when 1 biological barcode in the solution binds to 1 biotin, the signal intensity is 1, and so on, and when biological barcodes bind to biotin in the solution, the signal intensity is a, where a represents the coefficient in front of each variable, namely, the weight.

Mathematical Problems in Engineering
Step 5. Find a feasible solution satisfying the first constraint condition by detecting the number of biological barcodes.
Step 6. Repeat steps 4-6 above for the feasible solution obtained from the previous constraint condition, and we can obtain the feasible solution that satisfies all of the constraints.
Step 7. Calculate each feasible solution corresponding to the objective function value and, finally, judge the optimal integer programming solution.

Example Analysis.
A general 0-1 integer programming problem is discussed in detail as follows: Step 1. For the variables in each constraint condition, magnetic beads with capture DNA were designed. In this problem, 3, 2, and 2 magnetic beads with capture DNA were designed for the three constraints, with magnetic bead radii of 2 nm, 4 nm, and 6 nm, respectively. en, two types of hairpin structure DNA strands were designed, known as DNA1 and DNA2. When the value of variable x i is 1, the sticky ends of the DNA1 and DNA2 strands are modified with biotin. When the value of the variable x i is 0, the sticky ends of the DNA1 and DNA2 strands do not need to be modified with biotin.
Step 2. An appropriate amount of the biological barcode was placed in the solution and mixed evenly. ree sets of test tubes were prepared, and the number of each set of test tubes was 8, 4, and 4. e correct amount and equal amount of solution were placed into the test tubes.
ere are three test tubes in each set of test tubes, each of which is put into a magnetic bead with captured DNA, which are recorded as x 1 , x 2 , x 3 , respectively, and the radii of the magnetic beads are 2 nm, 4 nm, and 6 nm, Put DNA1 and DNA2 into the respective test tubes according to the x 1 , x 2 , x 3 coefficients in the constraint condition.
Step 4. e specific process is shown in Figure 4.
us, the feasible solutions satisfying the first two constraints are 4 (0,1,1) and 6 (1,0,1). Because the third constraint condition does not involve variable x 1 , we only need to consider the values x 2 and x 3 , and the values are 4 (1, 1) and 6 (0, 1). Continue to steps 4 and 5. e specific process is shown in Figure 6 below.
Step 7. Finally, group 6 of solutions (1, 0, 1) is a feasible solution that satisfies all constraints. Substituting the feasible solution into the objective function, the minimum objective function of the 0-1 integer programming problem can be obtained as 9.

Discussion
Visual DSD is a simulation software commonly used in DNA computing and hybridization chain reaction. is paper uses Visual DSD software to simulate and analyze the optimal solution of the 0-1 integer programming problem. e optimal solution of the example integer programming problem is (x 1 , x 2 , x 3 ) � (1, 0, 1). For variables x 1 � 1, add the hairpin structure DNA1 strand and DNA2 strand, and because the reaction is just started, the concentration of reactants is higher and the reaction speed is faster. e concentration of the hairpin structure DNA1 and DNA2 strands decreases rapidly in a short time and eventually gradually approaches 0. For sp5, the intermediate product of the reaction, because the cycle hybridization reaction is carried out step by step, the concentration of the strand first increases and then decreases before finally approaching 0. e concentration of the final product sp4 gradually increases and finally tends to be stable. e specific reaction process is shown in Figure 7. e simulation results show that the model is feasible and consistent with the expected results.
Previous models, such as the DNA origami base, circular logic gate, and others, cannot solve the weighted integer programming problem, which increases the understanding space virtually. e magnetic bead model proposed in this paper, which can solve the 0-1 integer programming problem with weight, can solve the general 0-1 integer programming problem, so it is more widely used. x 1 x 2 x 3 x 1 x 2 x 3 x 1 x 2 x 3 x 1 x 2 x 3 x 1 x 2 x 3 x 1 x 2 x 3 x 1 x 2 x 3 x 1 x 2 x 3 x 1 x 2 x 3 x 1 x 2 x 3 x 1 x 2 x 3 x 1 x 2 x 3 x 1 x 2 x 3 x 1 x 2 x 3 x 1 x 2 x 3 x 1 x 2 x 3 Mathematical Problems in Engineering

Conclusion
In this paper, a magnetic bead model for solving the 0-1 integer programming problem was established based on the DNA cycle hybridization chain reaction and the specific binding effect of streptavidin and biotin. Compared with the previous DNA computing model, this model has the following advantages. First, there is no requirement for enzymes in the operation process, which can reduce the experimental cost and improve the versatility of the model. Second, the intensity of the signal is used to judge the feasibility of the solution. is can improve the accuracy and practicability of the detection results. Finally, this model can be used to solve the general 0-1 integer programming problem and has more extensive applications than the previous DNA computing model. However, this method still has some shortcomings, such as a large number of steps and long operation time. erefore, these aspects still need to be studied further.

Data Availability
No data were used to support this study.

Conflicts of Interest
e authors declare that they have no conflicts of interest. x 2 x 3 x 2 x 3 x 2 x 3 x 2 x 3  x 1 x 2 x 1 x 2 x 1 x 2 x 1 x 2 x 1 x 2 x 1 x 2 Figure 5: Diagram of the solution under the second constraint. 6 Mathematical Problems in Engineering