Four-Analog Computation Based on DNA Strand Displacement

DNA strand displacement plays an important role in biological computations. The inherent advantages of parallelism, high storability, and cascading ha...
0 downloads 0 Views 3MB Size
This is an open access article published under an ACS AuthorChoice License, which permits copying and redistribution of the article or any adaptations for non-commercial purposes.

Article http://pubs.acs.org/journal/acsodf

Four-Analog Computation Based on DNA Strand Displacement Chengye Zou,† Xiaopeng Wei,*,† Qiang Zhang,*,† Chanjuan Liu,† Changjun Zhou,‡ and Yuan Liu† †

Faculty of Electronic Information and Electrical Engineering, Dalian University of Technology, Dalian 116024, China Key Laboratory of Advanced and Intelligent Computing, Dalian University, Ministry of Education, Dalian 116622, China



ABSTRACT: DNA strand displacement plays an important role in biological computations. The inherent advantages of parallelism, high storability, and cascading have resulted in increased functional circuit realization of DNA strand displacement on the nanoscale. Herein, we propose an analog computation with minus based on DNA strand displacement. The addition, subtraction, multiplication, and division gates as elementary gates could realize analog computation with minus. The advantages of this proposal are the analog computation with negative value and division computation. In this article, we provide the designs and principles of these elementary gates and demonstrate gate performance by simulation. Furthermore, to show the cascade property of gates, we computed a polynomial as an example by these gates.

1. INTRODUCTION The advancement of science and technology has necessitated higher requirements for computing speed and size of computer; however, these requirements may be limited by budget speed and computer size that are at maximal capacity. Therefore, new computer types are urgently needed to break through the bottleneck. Biological computers may satisfy the need in the future due to their parallelism and high storability. At present, biological calculation and DNA computation have attracted research interest, and increasing DNA circuits are being designed for both digital and analog computations.1−16 The branch migration progress of DNA strand displacement is suitable for the production of digital signal, and many digital DNA circuits are built by DNA strand displacement, even appearing large scale integrated DNA digital circuit.17,18 Compared with analog circuits, digital circuits can easily reduce the deviation; however, analog circuits have certain advantages. First, analog computation circuits consume less resource than digital computation circuits,1,2,5−7 for example the add computation of 31 + 7 for digital DNA circuit consumes 8 DNA strands to represent inputs, while an analog DNA circuit consumes only 2 DNA strands to represent inputs, thus it may have significance for experimental applications. Second, analog DNA circuits can dynamically reflect the concentration of the DNA strand, which is a function of the performance of the DNA circuit and the mutual relationships between multiple analog gates, which will be applied for analog control,19−21 for example DNA doctor22 and molecular robots.23 Third, there are binary bits in the DNA digital circuit, whose values are decided by a threshold; analog circuit is without threshold and analog circuit is more robust than digital circuit in some applications. © 2017 American Chemical Society

Song et al. have researched analog computation by DNA and strand displacement circuits, and they proposed three elementary arithmetic operations: addition, subtraction, and multiplication.5 Analog circuits are suitable for positive analog computation because the concentration of DNA strands is the value of the input and the output. The concentration of DNA strands can only be positive; thus, in their article, they have implemented only positive analog computation. The subtraction gate in particular permits the use of larger numbers and allows smaller values to be removed. In addition, their element gates are without the division gate; as a result, they utilize Newton iteration to realize the division computation, but the error is large. In this article, we proposed four new elementary arithmetic operations based on their three element gates, addition, subtraction, multiplication, and division, which can carry out analog computation with negative and positive values through a dual-rail system and compared the properties of our DNA analog circle with their design. The chemical reaction networks (CRNs) of the four element gates are originated from ref 5. We have simplified the CRNs of their three element gates and reduce reversible reactions.

2. RESULTS 2.1. Abstractions of the Gates. Every gate has four inputs and two outputs, the high corner marks + and to distinguish − here indicate positive and negative values. For example, when input values a1 and a2 are positive numbers, a+1 and a+2 are nonzero but a−1 and a−2 are zero; when input values a1 and a2 are negative numbers, a+1 and a+2 are zero but a−1 and a−2 are Received: May 8, 2017 Accepted: July 5, 2017 Published: August 2, 2017 4143

DOI: 10.1021/acsomega.7b00572 ACS Omega 2017, 2, 4143−4160

ACS Omega

Article

and Ca3 are the chemical species; we define the initial concentrations as [Gai]0 ≥ 2ra and [Ca3]0 ≥ 2ra. Ma1 and Sp2 are the output chemical species of this gate, and their concentrations at equilibrium, [Ma1]∞ and [Sp2]∞, represent P+a and P−a , respectively. When Ia+1 , Ia+2 , Ia−1 , and Ia−2 are reacted eventually, addition gate computers can be described as follows Pa = a1 + a 2 =[Ia1+]0 − [Ia1−]0 + [Ia +2 ]0 − [Ia −2 ]0 =[Ia1+]0 + [Ia +2 ]0 − ([Ia1−]0 + [Ia −2 ]0 ) =[Ma1]∞ − [Sp2 ]∞

Figure 1. Abstraction of analog DNA gates: addition gate, subtraction gate, multiplication gate, and division gate.

The output of the addition gate can be described as follows

nonzero; when output value Pa is positive, P+a is nonzero but P−a is zero; and when output value Pa is negative, P+a is zero but P−a is nonzero (Figure 1). 2.2. Addition Gate. CRNs 1a−1f5 of the addition gate can be described as follows

Ia1+

|Pa| = |a1 + a 2| =|[Ma1]∞ − [Sp2 ]∞ | ⎧ if [Ia1+]0 + [Ia +2 ]0 > ⎪ P+ = [Ma ] − [Sp ] a 1∞ 2 ∞ ⎪ [Ia1−]0 + [Ia−2 ]0 ⎪ =⎨ ⎪ if [Ia1+]0 + [Ia +2 ]0 < ⎪ Pa− = [Sp2]∞ − [Ma1]∞ ⎪ [Ia1−]0 + [Ia−2 ]0 ⎩

k1

+ Ga1 → Sp1

(1a)

k2

Ia +2 + Ga1 → Sp1

(1b)

k3

Ia1− + Ga 2 → Sp2 Ia−2

(1c) (1d)

k5

Sp1 +Ca3 → Ma1 + Ma 2

(1e)

k6

Sp2 +Ma1 → ⌀ Ia+1 ,

Ia+2 ,

Ia−1 ,

(3)

When [Ia+1 ]0 + [Ia+2 ]0 > [Ia−1 ]0 + [Ia−2 ]0, [Ma1]∞ is nonzero and [Sp2]∞ is zero; otherwise, [Ma1]∞ is zero and [Sp2]∞ is nonzero. There are no constraints on rate constants; ki (i = 1, ..., 6). The DNA reactions in the addition gate are shown in Figures 2 and 3. 2.3. Subtraction Gate. The subtraction gate is inspired by CRNs 5a−5f5

k4

+ Ga 2 → Sp2

(2)

(1f)

Ia−2

and are the input chemical species to the addition gate, where their initial concentrations [Ia+1 ]0, [Ia+2 ]0, [Ia−1 ]0, and [Ia−2 ]0 represent the four inputs Ia+1 , Ia+2 , Ia−1 , and Ia−2 , respectively. Therefore, a+1 = [Ia+1 ]0, a+2 = [Ia+2 ]0, a−1 = [Ia−1 ]0, and a−2 = [Ia−2 ]0, where a+1 , a+2 , a−1 , and a−2 ∈ (0, ra) and Gai (i = 1, 2)

k1

Is1+ + Ga1 → Sp1

(5a)

k2

Is−2 + Ga1 → Sp1

(5b)

Figure 2. Diagram of DNA reactions in the addition gate. The length of branch migration domains represented by gray domains is 20 nucleotides, the length of toeholds represented by colored domains is 5 nucleotides, and the arrow in a DNA strand indicates the 3′ to 5′ direction. The optimum temperature for the reactions is set at 25 °C.17,18 The forward and backward reaction rates are 2 × 10−3/nM/s. 4144

DOI: 10.1021/acsomega.7b00572 ACS Omega 2017, 2, 4143−4160

ACS Omega

Article

Figure 3. List of the DNA reactions in the addition gate; CRN 4 is adapted from ref 5.

where Is+1 , Is+2 , Is−1 , and Is−2 are the inputs of the subtraction gate; when they are reacted eventually, subtraction gate computers can be described as follows

k3

Is1− + Ga 2 → Sp2

(5c)

k4

Is+2 + Ga 2 → Sp2

Ps = s1 − s2

(5d)

=[Is1+]0 − [Is1−]0 − ([Is+2 ]0 − [Is−2 ]0 ) k5

Sp1 +Ca3 → Ma1 + Ma 2

=[Is1+]0 + [Is−2 ]0 − ([Is1−]0 + [Is+2 ]0 )

(5e)

=[Ma1]∞ − [Sp2 ]∞

k6

Sp2 +Ma1 → ⌀

(5f)

(6)

The output of the addition gate can be described as follows 4145

DOI: 10.1021/acsomega.7b00572 ACS Omega 2017, 2, 4143−4160

ACS Omega

Article k5

|Ps| = |s1 − s2|

G1 + Sy1 → Gy1 + Ng1

=|[Ma1]∞ − [Sp2 ]∞ | ⎧ if [Is1+]0 + [Is−2 ]0 > ⎪ P+ = [Ma ] − [Sp ] 1∞ 2 ∞ ⎪ s [Is1−]0 + [Is+2 ]0 ⎪ =⎨ ⎪ if [Is1+]0 + [Is−2 ]0 < − ⎪ P = [Sp ]∞ − [Ma1]∞ 2 ⎪ s [Is1−]0 + [Is+2 ]0 ⎩

(10a)

k6

Gd1 + Gy1 → Idn6 + Sg1

(10b)

k7

Nd1 + Ng1 → ⌀

(10c)

k8

Nd1 + Gd1 → ⌀

(10d)

k9

(7)

G2 + Sy2 → Gy2 + Ng 2

When [Is+1 ]0 + [Is−2 ]0 > [Is−1 ]0 + [Is+2 ]0, [Ma1]∞ is nonzero and [Sp2]∞ is zero; otherwise, [Ma1]∞ is zero and [Sp2]∞ is nonzero. There are no constraints on rate constants, ki (i = 1, ..., 6). The DNA reactions in the subtraction gate are similar to those in the addition gate; therefore, we exclude the diagram and the list of DNA reactions in the subtraction gate. 2.4. Multiplication Gate. Considering the symbols + and −, multiplication computation is classified as follows

(10e)

k10

Gd 2 + Gy2 → Idn6 + Sg 2

(10f)

k11

Nd 2 + Ng 2 → ⌀

(10g)

k12

Nd 2 + Gd 2 → ⌀

(10h)

k13

I1 + Sx1 → Gx1 + Nh1

Pm = m1m2 ⎧ + +⎫ ⎪ m1 m2 ⎪ ⎬ = Pm+ ⎪ m −m − ⎪ ⎪ 1 2⎭ =⎨ ⎫ ⎪ m1+m2− ⎪ ⎪ ⎬ = Pm− ⎪ m1−m2+ ⎪ ⎭ ⎩ ⎧ ⎧ m +m + ⎫ ⎪ ⎪ ̃1 ̃ 2 ⎪ ⎬ = Pm+ ⎪ ⎪ m̃ −m̃ − ⎪ ⎪ 1 2⎭ ⎪ ⎨ (8a) ⎪ ⎫ ⎪ m̃ 1+m̃ 2− ⎪ ⎪ ⎪ ⎬=0 ⎪ ⎪ − +⎪ ⎪ ⎩ m̃ 1 m̃ 2 ⎭ ⎪ =⎨ ⎧ m̂ +m̂ + ⎫ ⎪ ⎪ 1 2⎪ ⎪ ⎬=0 ⎪ m̂ −m̂ − ⎪ ⎪ 1 2 ⎭ ⎪ ⎪(8b)⎨ ⎪ ⎫ ⎪ m̂ 1+m̂ 2− ⎪ ⎪ ⎪ ⎬ = Pm− ⎪ ⎪ m̂ 1−m̂ 2+ ⎪ ⎪ ⎭ ⎩ ⎩

k14

Ht1 + Gx1 → mt + Sm1 k15

k16

k17

(11e)

k18

Ht 2 + Gx2 → mt + Sm2

(11f)

k19

Nt 2 + Nh 2 → ⌀

(11g)

k 20

Nt 2 + Ht 2 → ⌀

(11h) k 21

Idn6 + amplifier1 →  At + At + ··· + At   rm At

(12a)

k 22

mt + amplifier2 →  Bt + Bt + ··· +Bt   rm Bt

(8)

Im+1 ,

Im−1 ,

Im+2 ,

(12b)

Im−2

where and are the input chemical species to the gate; At and Bt are the output chemical species to the gate; and Gd1, Gd2, Ht1, Ht2, G1, G2, I1, and I2 are the breakup products of the input chemical species. [Im+1 ]0, [Im−1 ]0, [Im+2 ]0, and [Im−2 ]0 are the initial concentrations of the input chemical species, so that, [Im+1 ]0 = m+1 , [Im−1 ]0 = m−1 , [Im+2 ]0 = m+2 , and [Im−2 ]0 = m−2 . Sy1, Sy2, Nd1, Nd2, Sx1, Sx2, Nt1, Nt2, amplifier1, amplifier2, splitter1, splitter2, splitter3, and splitter4 are the composed chemical species of the multiplication gate, and (0, rm) is the input range; this implies that m+1 , m−1 , m+2 , m−2 ∈ (0, rm), and initial concentrations of the other composed chemical species Sy1, Sy2, Nd1, Nd2, Sx1, Sx2, Nt1, Nt2, splitter1, splitter2, splitter3, splitter4, amplifier1, and amplifier2 are set to rm, where rm = 10N (N = 0, 1, 2, 3, ...). When Im+1 , Im+2 , Im−1 , and Im−2 are reacted at equilibrium, [At]∞ and [Bt]∞ represent P+m and P−m, respectively. All other chemical species are intermediate products. According to eq 8, eqs 8a and 8b cannot be satisfied simultaneously. When the inputs are Im+1 and Im+2 or Im−1 and

k1

(9a)

k2

(9b)

k3

(9c)

k4

Im−2 + splitter4 → G2 + I 2

(11d)

I 2 + Sx2 → Gx2 + Nh 2

CRNs 9a−9d, 10a−10h,5 11a−11h,5 and 12a,12b of the multiplication gate are given by

Im+2 + splitter3 → G1 + I1

(11c)

Nt1 + Ht1 → ⌀

Figure 4. Splitting of four inputs in the multiplication gate.

Im1− + splitter2 → Gd 2 + Ht1

(11b)

Nt1 + Nh1 → ⌀

According to eq 8, we should split four inputs into eight inputs, as shown in Figure 4.

Im1+ + splitter1 → Gd1 + Ht 2

(11a)

(9d) 4146

DOI: 10.1021/acsomega.7b00572 ACS Omega 2017, 2, 4143−4160

ACS Omega

Article

Figure 5. (a) Diagram of DNA reactions in splitter1 and splitter2. CRNs of splitter1 and splitter2 are adapted from ref 4. (b) Diagram of DNA reactions in splitter3 and splitter4. CRNs of splitter3 and splitter4 are adapted from ref 4.

Im−2 , the products of splitter are Gd1, G1, Ht2, and I1 or Gd2, G2, Ht1, and I2, respectively, which can produce Idn6. This means that [At]∞ is nonzero, but [Bt]∞ is zero; otherwise, [At]∞ is zero, but [Bt]∞ is nonzero. The reaction rate constant must meet the following requirements

⎧ k1 , ⎪ ⎪ k1 , ⎪ ⎨ k1 , ⎪ ⎪ k1 , ⎪ ⎩ k1 ,

Because k6 = k7, Gd1 will be consumed by Gy1 and Nd1 at the same time and the concentration of Idn6 at equilibrium is [Idn6]∞ = m1+

k 2 ≪ k5 = k 7 [At]∞ = rm

k 2 ≪ k 9 = k11 k 2 ≪ k13 = k15

[Nd1]

=

[Im+2 ]0 m2+ = rm − [Im+2 ]0 rm − m2+

(15)

m1+m2+ = [Im1+]0 [Im+2 ]0 rm

(16)

Similarly, we have ⎧ m1−m2− = [Im1−]0 [Im−2 ]0 ⎪[At]∞ = rm r m ⎪ ⎪ ⎪ m +m − ⎨[Bt]∞ = rm 1 2 = [Im1+]0 [Im−2 ]0 rm ⎪ ⎪ m1−m2+ ⎪ = [Im1−]0 [Im+2 ]0 ⎪[Bt]∞ = rm r ⎩ m

(13)

Therefore, reactions 10a, 10c, 10e, 10g, 11a, 11c, 11e, 11e, and 11g must be completed before reactions 9a and 9b can start. Because reactions 10a and 10c are completed before reaction 9a, the concentration ratio of Gy1 and Nd1 satisfied [Gy1]

m2+ m1+m2+ + = + (rm − m2 ) rm

Then, to take advantage of amplification reaction 38, the concentration of At at equilibrium will be

k 2 ≪ k 3 , k4

k 2 ≪ k17 = k19

m2+

(14) 4147

(17) DOI: 10.1021/acsomega.7b00572 ACS Omega 2017, 2, 4143−4160

ACS Omega

Article

Figure 6. Diagram of DNA reactions in the 2× amplifier. CRNs of the amplifier are adapted from ref 4. k2

The diagrams of the DNA reactions of the splitter and amplifier are shown in Figure 5 and Figure 6, respectively. Figure 6 shows the reaction diagrams of the two-fold amplifier; we can increase the length of DNA strands of Pm1, Pm2, Qm1, Qm2, Qd1, and Qd2 to achieve multiple amplification. The DNA reactions in the multiplication gate are shown in Figure 7. 2.5. Division Gate. 2.5.1. Basic Division Gate. For the division computation with minus, division computations are classified as follows

Id1− + split 2 → Id 2 + Ih1 k3

Id +2 + split3 → Xb2 + Xm2

(21c)

k4

Id−2 + split4 → Md 2 + Bm2

(21d)

k5

Id1 + ID1 → Xg1 + Ng1 + Xb1

⎧ d +/ d + ⎫ ⎪ ⎪ 1 2⎬ = Pd+ ⎪ ⎪ d −/d − ⎭ ⎪ 1 2 Pd = d1d 2 = ⎨ ⎪ d1+/d 2− ⎫ ⎪ ⎪ ⎬ = Pd− ⎪ − + ⎪ d1 /d 2 ⎭ ⎩

(22a)

k6

Gi1 + Xb1 → Bi1 + Ig1 k7

Xb2 + Bi1 → ⌀

(22b) (22c)

k8

IM1 + Ng1 → ⌀

(22d)

k8

⎧ ⎧ ̃ + ̃ +⎫ ⎪ ⎪ ⎪ d1 /d 2 ⎬ + ⎪ ⎪ ̃ − ̃ − ⎪ = Pd ⎪ ⎪ d1 /d 2 ⎭ ⎪(20a)⎨ ⎪ ⎪ d ̃ +/ d ̃ − ⎫ ⎪ ⎪ ⎪ 1 2⎬=0 ⎪ ⎪ d ̃ −/d ̃ + ⎪ ⎩ 1 2⎭ ⎪ =⎨ ⎧ ̂ + ̂ +⎫ ⎪ ⎪ d1 /d 2 ⎪ ⎪ =0 ⎪ − −⎬ ⎪ ⎪ ̂ ̂ / d d ⎪ 1 2⎭ ⎪ ⎪(20b)⎨ + − ⎪ d ̂ /d ̂ ⎫ ⎪ ⎪ 1 2⎪ ⎪ ⎬ = Pd− ⎪ ̂ − ̂ +⎪ ⎪ ⎩ d1 /d 2 ⎭ ⎩

IM1 + Xb2 → Pg1 + I 2d1

(22e)

k9

Id 2 + ID2 → Xg 2 + Ng 2 + Md1

(22f)

k10

Gi 2 + Md1 → Bi 2 + Ig 2 k11

Md 2 + Bi 2 → ⌀

(22g) (22h)

k12

IM 2 + Ng 2 → ⌀

(22i)

k12

IM 2 + Md 2 → Pg 2 + I 2d1 k13

Ih1 + IH1 → Xh1 + Nh1 + Xm1 (20)

k14

Xm1 + Hm1 → Bh1 + Mh1

where |d2| ≥ |d1|, |d2| ≥ 1. According to eq 20, we should split four inputs into eight inputs, as shown in Figure 8. The following CRNs, 21a−21d, 22a−22j,5 23a−23j,5 and 24a−24d, are given for the division gate

k15

Xm2 + Bh1 → ⌀ k16

IM3 + Nh1 → ⌀

k1

Id1+ + split1 → Id1 + Ih 2

(21b)

k16

IM3 + Xm2 → Qb1 + M 2t1

(21a) 4148

(22j) (23a) (23b) (23c) (23d) (23e) DOI: 10.1021/acsomega.7b00572 ACS Omega 2017, 2, 4143−4160

ACS Omega

Article k17

Ih 2 + IH 2 → Xh 2 + Nh 2 + Bm1 k18

Bm1 + Hm2 → Bh 2 + Mh 2

k19

(23f)

Bm2 + Bh 2 → ⌀

(23g)

IM4 + Nh 2 → ⌀

k 20

(23h) (23i)

Figure 7. continued 4149

DOI: 10.1021/acsomega.7b00572 ACS Omega 2017, 2, 4143−4160

ACS Omega

Article

Figure 7. List of the DNA reactions in the multiplication gate; CRNs 19 and 20 are adapted from ref 5. k 20

IM4 + Bm2 → Qb2 + M 2t1 k 21

An + I 2d1 → Pn + N2a

Figure 8. Splitting of four inputs in the division gate. 4150

(23j)

(24a) DOI: 10.1021/acsomega.7b00572 ACS Omega 2017, 2, 4143−4160

ACS Omega

Article

Figure 9. continued

4151

DOI: 10.1021/acsomega.7b00572 ACS Omega 2017, 2, 4143−4160

ACS Omega

Article

Figure 9. List of the DNA reactions in the basic division gate; CRNs 30−32 are adapted from refs 4, 5. k 22

I 2d 2 + Pn → ⌀ k 23

Bn + M 2t1 → Qn + N1b k 24

M 2t 2 + Qn ⎯→ ⎯ ⌀

where Id+1 , Id−1 , Id+1 , and Id−2 are the input chemical species to the gate; I2d2 and M2t2 are the output chemical species to the gate; and Id1, Ih2, Id2, Ih1, Xb2, Xm2, Md2, and Bm2 are the breakup products of the input chemical species. [Id+1 ]0, [Id−1 ]0, [Id+2 ]0, and [Id−2 ]0 are the initial concentrations of the input chemical species, so that, [Id+1 ]0 = d+1 , [Id−1 ]0 = d−1 , [Id+2 ]0 = d+2 , and [Id−2 ]0 = d−2 , where ID1, ID2, Gi1, Gi2, IH1, IH2, Hm1, Hm2,

(24b)

(24c)

(24d) 4152

DOI: 10.1021/acsomega.7b00572 ACS Omega 2017, 2, 4143−4160

ACS Omega

Article

Figure 10. Main leak reactions according to ref 5 in the four analog computation gates: leak reactions 1a and 2 in the addition and subtraction gates; leak reactions 3−10a in the multiplication gate; and leak reactions 11a−14 in the basic division gate. The forward and backward leak reaction rates are 5 × 10−9/nM/s.

the input range, meaning that d+1 , d−1 , d+2 , d−2 ∈ (0, rm), and initial concentrations of the other composed chemical species are

IM1, IM2, IM3, IM4, I2d2, M2t2, split1, split2, split3, and split4 are composed of the chemical species of the division gate, (0, rm) is 4153

DOI: 10.1021/acsomega.7b00572 ACS Omega 2017, 2, 4143−4160

ACS Omega

Article

defined as [ID1]0 = [ID2]0 = [Gi1]0 =[ Gi2]0 = [IH1]0 = [IH2]0 = [Hm1]0 = [Hm2]0 = [An]0 = [Bn]0 = rm, [IM1]0 = [IM3]0 = [IM4]0 = [I2d2]0 = [M2t2]0 = 1. When Id+1 , Id+2 , Id−1 , and Id−2 are reacted at equilibrium, [I2d2]∞ and [M2t2]∞ represent P+d and P−d , respectively. All other chemical species are intermediate products. According to eq 20, eqs 20a and 20b cannot be satisfied simultaneously. When the inputs are Id+1 and Id+2 , the products of the splitter are Id1, Ih2, Xb2, and Xm2 and the concentration of M2t1 at equilibrium will be [M2t1]∞ = 1 corresponding to reaction 23e, where the reaction rate constant must meet the following requirements ⎧ k 8 ≪ k1 , k 3 , k 5 , k 6 , k 7 ⎪ ⎪ k12 ≪ k 2 , k4 , k 9 , k10 , k11 ⎨ ⎪ k16 ≪ k1 , k4 , k13 , k14 , k15 ⎪ ⎩ k 20 ≪ k 2 , k 3 , k17 , k18 , k19

The DNA reactions in the basic division gate are shown in Figure 9. 2.5.2. Improved Division Gate. In consideration of the limitations (|d2| ≥ |d1|, |d2| ≥ 1), we have taken the following actions to improve the analog computation of division 1. To satisfy 10N|d2| ≥ |d1| and 10N|d2| ≥ 1 (N = 0, 1, 2, 3, ...), we use a 10N-times amplifier to amplify the concentration of the input of Id+2 or Id−2 ; when N = 0, the improved division gate will be reduced to the basic division gate. 2. Amplified Id+2 or Id−2 is involved in the reactions of the splitter and basic division gate; thus, the concentration of I2d1 at equilibrium is [I 2d1]∞ = 1

⎧ − ⎞ − ⎛ ⎪[I 2d 2]∞ = 1 − ⎜1 − d1 ⎟ = d1 ⎪ 10 N d 2− ⎠ 10 N d 2− ⎝ ⎪ ⎛ ⎪ d+ ⎞ ⎪[M 2t 2]∞ = 1 − ⎜1 − N1 − ⎟ 10 d 2 ⎠ ⎝ ⎪ + ⎪ d ⎨ = N1 − ⎪ 10 d 2 ⎪ ⎛ ⎪ d1− ⎞ ⎪[M 2t 2]∞ = 1 − ⎜1 − N + ⎟ 10 d 2 ⎠ ⎝ ⎪ d1− ⎪ = N + ⎪ 10 d 2 ⎩

(26)

(27)

(35)

3. Finally, we use a 10N-times amplifier to amplify the output of the basic division gate. Therefore, we can extend the division analog computation to a real number range through the above actions.

(28)

Similarly, we have ⎧ ⎛ d− ⎞ ⎪[I 2d 2]∞ = 1 − ⎜1 − 1− ⎟ d2 ⎠ ⎪ ⎝ ⎪ d1− = − , [M 2t 2]∞ ⎪ d2 ⎪ ⎪ =0 ⎪ ⎛ d1+ ⎞ ⎪ = − − [M t ] 1 1 ⎟ ⎜ ⎪ 22∞ d 2− ⎠ ⎝ ⎪ ⎨ d+ = 1− , [I 2d 2]∞ ⎪ d2 ⎪ ⎪ =0 ⎪ ⎛ d1− ⎞ ⎪ = − − [M t ] 1 1 ⎟ ⎜ ⎪ 22∞ d 2+ ⎠ ⎝ ⎪ d− ⎪ = 1+ , [I 2d 2]∞ ⎪ d2 ⎪ ⎩ =0

(34)

Similarly, we have

Next, to take advantage of amplification reactions 24a−24d, the concentrations of I2d2 and M2t2 at equilibrium will be ⎧ ⎛ d1+ ⎞ d1+ ⎪ ⎪[I 2d 2]∞ = 1 − ⎜1 − + ⎟ = + d2 ⎠ d2 ⎝ ⎨ ⎪ ⎪[M t ] = 1 − 1 = 0 ⎩ 22∞

10 N d 2+

⎛ d+ ⎞ d+ [I 2d 2]∞ = 1 − ⎜1 − N1 + ⎟ = N1 + 10 d 2 ⎠ 10 d 2 ⎝

Owing to the rate of reactions 22d and 22e being identical, IM1 will be consumed by Ng1 and Xb2 at the same time and the concentration of I2d1 at equilibrium is d 2+ − d1+ d1+ = − 1 d1+ + (d 2+ − d1+) d 2+

d1+

Then, the concentration of I2d2 at equilibrium will be

[Ng1]

[I 2d1]∞ = 1

d1+ + (10 N d 2+ − d1+)

=1−

(33) (25)

Therefore, reactions 21a, 21c, 22a, 22b, and 22c must be completed before reactions 22d and 22e can start. Because reactions 22b and 22c are completed before reactions 22d and 22e, the concentration of Xb2 approaches [Id+2 ]0 − [Id+1 ]0 and then the concentration ratio of Ng1 and Xb2 is satisfied [Id1+]0 d1+ = = [Xb2 ] [Id +2 ]0 − [Id1+]0 d 2+ − d1+

10 N d 2+ − d1+

3. LEAK REACTIONS Figure 10 depicts the main leak reactions in the four analog computation gates. The reason of the main leak is that the base pairs in the circular portion of DNA strand can be temporarily broken and create a toehold for an invaded strand. For example, Ma2 can invade Sx1 and displace the q1 domain in the addition and subtraction gates, which is a typical reaction existing in the multiplication and division gates without splitter and amplifier. We have analyzed the leak reactions of the four gates to evaluate the affection of the leak reactions for the output of the gates as follows: (a) Affection of leak reactions in addition gates: (I) When the symbols of a1 and a2 are different, leak reaction 1a in the addition gate is absent. The reaction rate of leak reaction 2 (reaction 37) is reduced compared to that of reaction 1f; therefore, affection of leak reaction 2 can be ignored.

(29) 4154

DOI: 10.1021/acsomega.7b00572 ACS Omega 2017, 2, 4143−4160

ACS Omega

Article

Figure 11. Performance of the addition gate.

gate, multiplication gate, and division gate, respectively, where 0.02 ≤ r ≤ 0.04 nM.

(II) When the symbols of a1 and a2 are the same, leak reaction 1a (reaction 38) will delay the concentration of the output reaching to valid range. slow ⎧ ⎪ Sx + Ma 2 ⎯⎯⎯→ Aq1 + Sp1 ⎨ 1 ⎪ Aq + Sp → Sx1 + Ma 2 ⎩ 1 1

4. SIMULATION RESULT OF THE GATES To test the effectiveness of these four gates, we simulated three input ranges (0, 1), (0, 2), and (0, 4), where the ranges between the pink and green lines were the valid output ranges. The simulation performance for the addition gate is given in Figure 11. When the symbols of a1 and a2 are different, the output stays in the valid range for a longer time because the inputs are larger and the effect of leak is relatively smaller. When the symbols of a1 and a2 are the same, the period for outputs in the valid range is constant. Figure 12 shows the performance of the simulation for the subtraction gate. When the symbols of s1 and s2 are the same, the output stays in the valid range for a longer time because the inputs are larger and the effect of leak is relatively smaller. When the symbols of s1 and s2 are different, the period for outputs to stay in the valid range is constant. The simulation performance for the multiplication gate is shown by Figure 13; because the period for outputs to stay in the valid range is constant, the influence of leak for multiplication gates can be ignored. Figure 14 shows the simulation performance for the division gate; the period for output to stay in the valid range increases with increasing concentrations of inputs because the influence of leak decreases with increasing inputs.

(36)

slow

Ca3 + Sp2 ←→ Aq1 + q2

(37)

(b) Effect of leak reactions in subtraction gates: (I) When the symbols of s1 and s2 are the same, the affection of leak reaction can be ignored. (II) When the symbols of s1 and s2 are different, leak reaction 1a will delay the concentration of the output strand at equilibrium without affecting the value of output. (c) Effect of leak reactions in multiplication gates: Although the production of leak reactions 3, 5a, 7, and 9a will increase the value of output, we can neglect the affection of leak reactions in multiplication gates because these leak reactions have much less reaction rates compared to those of reactions 10b, 10f, 11b, and 11f. (d) Effect of leak reactions in division gates: In view of surplus of Xb2, Md2, Xm2, and Bm2 in reactions 22c and 22e; reactions 22h and 22j; reactions 23c and 23e; reactions 23h and 23j, respectively, the production of leak reactions 11a−14 will increase the concentrations of [I2d1] and [M2t1] slowly, which will reduce the value of output gradually. Above all, we use a valid range to show the performance of a gate under particular inputs; thus, pa − r ≤ output value ≤ pa + r, ps − r ≤ output value ≤ ps + r, pm − r ≤ output value ≤ pm + r, and pd − r ≤ output value ≤ pd + r are fixed to define the valid output range of the addition gate, subtraction

5. ANALOG DNA CIRCUIT TO COMPUTE xy g ( x , y) = x 2 + y 2 5.1. Principle of the Analog Circuit. The input and output strands have same properties, so our four analog gates are modular; therefore, we can build DNA circuits by the four analog gates. For the addition and subtraction gates, early arrivals will wait for latecomers. When the cancellation between Sp1 and Sp2 is finished, the remaining Sp1 or Sp2 will be the 4155

DOI: 10.1021/acsomega.7b00572 ACS Omega 2017, 2, 4143−4160

ACS Omega

Article

Figure 12. Performance of the subtraction gate.

Figure 13. Performance of the multiplication gate.

reactions, such that they can wait for others. It should be noted that the division gate should be installed at the end of the DNA circuit because I2d2 and M2t2 are statically prepared in advance and they will react with other DNA strands of below gates. Figure 15 shows an analog DNA circuit to compute xy g (x , y) = 2 2 , for x ∈ R, |y| ≥ 1.

output strands. For the multiplication gate, the concentration ratios of Gy1 and Nd1, Gy2 and Nd2, Gx1 and Nt1, and Gx2 and Nt2 are obtained as early as possible because the ratios of reactions 18 and 19 are much smaller than those of the other reactions. Furthermore, Im+1 and Im−1 or Im+2 and Im−2 can be freely chosen for preparation in a “dynamic” fashion by other gates. For the division gate, inputs IM1, IM2, IM3, and IM4 can be prepared in a “static” fashion and ratios of reactions with IM1, IM2, IM3, and IM4 are much smaller than those of other

x +y

5.2. Simulation of the Circuit to Compute g(x, y). In the simulation of g(x, y), x ∈ {−2, −1.75, −1.5, −1.25, −1, 4156

DOI: 10.1021/acsomega.7b00572 ACS Omega 2017, 2, 4143−4160

ACS Omega

Article

Figure 14. Performance of the basic division gate.

circuit, we fixed the valid range between 0.95g(x, y) and 1.05g(x, y) during the 106 s. Figure 16a gives an example of the computation of g(x, y), where x = 1 and y = −2. The figure of evolution of g(x, y) is similar to division gate, but the outputs of g(x, y) reach the valid range slower than division gate ([Id+1 ]0 = 1, [Id−1 ]0 = [Id+2 ]0 = 0, [Id−2 ]0 = 2), which is the result of operations of add and multiplication gates. In simulation, when y ∈ {1, 1.5, 2, 2.5, 3}, the outputs stay in the valid range for a longer time, as shown in Figure 16b, which is irrelevant to the value of outputs because the affection of leak reaction is reduced with enlargement of inputs. Consequently, time for outputs to stay in the valid range has positive and negative symmetries, which means that time for outputs to stay in the valid range is the same for identical absolute values of inputs. When y∈ {4, 5}, the outputs cannot stay in the valid range, meaning that the enlargement of error for outputs, rather than decrease of time for outputs, stays in the valid range.

6. COMPARISON AND ANALYSIS

Figure 15. Circuit to compute g (x , y) =

xy x2 + y2

In this article, we not only added the division gate but also realized analog computation by DNA strand displacement circuits with minus; furthermore, we made some improvements in DNA analog computation. 6.1. Reduce the Influence of Leak Reactions of the Amplifier. As Figure 18 shows, Om is a production of main leak reactions in the amplifier designed by Song et al. and will increase the output of the multiplication gate. The larger [A1]0 is more influential to the output, as shown in Figure 17, when [Im1]0 = 1.5 and [Im2]0 = 2. The leaks shown in Figure 18 were not evident in our amplifier, indicating that our multiplication gate was more stable and easier to operate. 6.2. Accuracy of Analog Division Computation. Although there was no division gate in the DNA analog com-

.

−0.75, −0.5, −0.25, 0, 0.25, 0.5, 0.75, 1, 1.25, 1.5, 1.75, 2} and y ∈ {1, 1.5, 2, 2.5, 3, 4, 5}. To quantify the performance of the 4157

DOI: 10.1021/acsomega.7b00572 ACS Omega 2017, 2, 4143−4160

ACS Omega

Article

putation gates designed by Song et al., they achieved division 1 computation r(x) = x , (0.5 < x < 1) by Newton iteration Yn + 1 = 2Yn − Yn2x

where lim Yn = n →∞

(38)

1 . x

Execution of the circuit to compute g(x, y) = 2y − y2x is shown in Figure 19, when x = 0.5 and y = 1.

Figure 19. Execution of the circuit to compute g(x, y).

Considering that 101x > 1 in the computation of r(x), we can use the improved division gate to compute r(x). The performance of the improved division gate is shown in Figure 20a, when [Id+1 ]0 = 1 and [Id+2 ]0 = 0.5. Then, we chose a 10× amplifier to amplify the output of the improved division gate when the output of the improved division gate reached a valid range; otherwise, I2d2 and M2t2 will be amplified first because I2d2 and M2t2 are statically prepared in advance. The process of amplification is shown in Figure 20b, where [I2d2]∞ = 0.1977 and I2d3 is the output of the 10× amplifier. On the basis of these principles, the results of our strategy were more accurate than those of the technique by Song et al.

Figure 16. (a) Evolution of DNA circuit to compute g(x, y) when x = 1 and y = −2. (b) Performance of the DNA circuit to compute g(x, y).

7. CONCLUSIONS We proposed four DNA analog computation gates and extended the computation range to real numbers. On the basis of the same properties of these gates, we constructed a DNA circuit to compute the polynomial function of inputs. Simulations showed that the time for outputs of circuit to reach a valid range was similar to that for a single basic division gate because DNA reactions are simultaneous; therefore, the

Figure 17. Evolution of the output with the concentration of DNA strand [A1]0, when [Im1]0 = 1.5 and [Im2]0 = 2.

Figure 18. Main leak reactions in the amplifier designed by Song et al.5 4158

DOI: 10.1021/acsomega.7b00572 ACS Omega 2017, 2, 4143−4160

ACS Omega



ACKNOWLEDGMENTS



REFERENCES

Article

This work is supported by the National Natural Science Foundation of China (Nos. 61425002, 61672121, 61572093, 61402066, 61402067, 61370005, and 31370778), the Program for Changjiang Scholars and Innovative Research Team in University (No. IRT_15R07), the Program for Liaoning Innovative Research Team in University (No. LT2015002), the Basic Research Program of the Key Lab in Liaoning Province Educational Department (Nos. LZ2014049 and LZ2015004), the Scientific Research Fund of Liaoning Provincial Education (Nos. L2015015 and L2014499), and the Program for Liaoning Key Lab of Intelligent Information Processing and Network Technology in University.

(1) Li, W.; Zhang, F.; Yan, H.; Liu, Y. DNA based arithmetic function: a half adder based on DNA strand displacement. Nanoscale 2016, 8, 3775−3784. (2) Sun, J. W.; Li, X.; Cui, G. Z.; Wang, Y. F. One-Bit Half AdderHalf Subtractor Logical Operation Based on the DNA Strand Displacement. J. Nanoelectron. Optoelectron. 2017, 12, 375−380. (3) Fern, J.; Scalise, D.; Cangialosi, A.; Howie, D.; Potters, L.; Schulman, R. DNA Strand-Displacement Timer Circuits. ACS Synth. Biol. 2017, 6, 190−193. (4) Lakin, M. R.; Stefanovic, D. Supervised Learning in Adaptive DNA Strand Displacement Networks. ACS Synth. Biol. 2016, 5, 885− 897. (5) Song, T.; Garg, S.; Mokhtar, R.; Bui, H.; Reif, J. Analog Computation by DNA Strand Displacement Circuits. ACS Synth. Biol. 2016, 5, 898−912. (6) Li, M. H.; Liu, F. F.; Song, M. In A Half-Subtracter Calculation Model Based on Stand Displacement Technology, 10th International Conference on Bio-Inspired Computing - Theories and Applications, Hefei, China, Sept 25−28, 2015. (7) Sun, J. W.; Li, X.; Huang, C.; Cui, G. Z.; Wang, Y. F. Two-Digit Full Subtractor Logical Operation Based on DNA Strand Displacement. In Bio-Inspired Computing-Theories and Applications. Pan, L. Q., Păun, G., Pérez-Jiménez, M. J., Song, T., Eds.; Academic Press, 2016; pp 21−29. (8) Chen, X. X.; Dong, Y. F.; Xiao, S. Y.; Liang, H. J. DNA and DNA computation based on toehold-mediated strand-displacement reactions. Acta Phys. Sin. 2016, 65, No. 178106. (9) Zhang, Z.; Fan, T. W.; Hsing, I. M. Integrating DNA strand displacement circuitry to the nonlinear hybridization chain reaction. Nanoscale 2017, 9, 2748−2754. (10) Cui, G.; Zhang, J.; Cui, Y.; Zhao, T.; Wang, Y. DNA StrandDisplacement Digital Logic Circuit with Fluorescence Resonance Energy Transfer Detection. J. Comput. Theor. Nanosci. 2015, 12, 2095−2100. (11) Deng, W.; Xu, H. G.; Ding, W.; Liang, H. J. DNA Logic Gate Based on Metallo-Toehold Strand Displacement. PLoS One 2014, 9, No. e111650. (12) Dou, B.; Yang, J.; Shi, K.; Yuan, R.; Xiang, Y. DNA-mediated strand displacement facilitates sensitive electronic detection of antibodies in human serums. Biosens. Bioelectron. 2016, 83, 156−161. (13) Giuffrida, M. C.; Zanoli, L. M.; D’Agata, D.; Finotti, A.; Gambari, R.; Spoto, G. Isothermal circular-strand-displacement polymerization of DNA and microRNA in digital microfluidic devices. Anal. Bioanal. Chem. 2015, 407, 1533−1543. (14) Rogers, W. B.; Manoharan, V. N. Programming colloidal phase transitions with DNA strand displacement. Science 2015, 347, 639− 642. (15) Song, M.; Wang, Y. C.; Li, M. H.; Dong, Y. F. Metal-Mediated Logic Computing Model with DNA Strand Displacement. J. Comput. Theor. Nanosci. 2015, 12, 2318−2321.

Figure 20. (a) Performance of the improved division gate. (b) Output of the basic division gate is amplified by the 10× amplifier.

computations in the DNA circuit are parallel, which we aimed to achieve in this study. Leak reaction is a common issue in a DNA strand displacement circuit,24 which was reduced in our work. In our amplifier and splitter, we eliminated the main leak reactions, similar to that of the main leak reactions in the amplifier designed by Song et al., which improved the computation of multiplication.

8. METHODS The performances of the four element gates, the improved division gate, and the computation of a polynomial are simulated by Language for Biochemical Systems (LBS). To speed up the simulation, we used MATLAB to run the code produced from LBS by Visual GEC. All DNA figures are drawn by Visual DSD. The unit of the values of the code in the supporting information is nM.



AUTHOR INFORMATION

Corresponding Authors

*E-mail: [email protected] (X.W.). *E-mail: [email protected] (Q.Z.). ORCID

Qiang Zhang: 0000-0003-3776-9799 Author Contributions

The study was carried out and the manuscript was written with contributions of all authors. All authors have approved the final version of the manuscript. Notes

The authors declare no competing financial interest. 4159

DOI: 10.1021/acsomega.7b00572 ACS Omega 2017, 2, 4143−4160

ACS Omega

Article

(16) Zhang, X. C.; Zhang, W. W.; Zhao, T. T.; Wang, Y. F.; Cui, G. Z. Design of Logic Circuits Based on Combinatorial Displacement of DNA Strands. J. Comput. Theor. Nanosci. 2015, 12, 1161−1164. (17) Yordanov, B.; Kim, J.; Petersen, R. T.; Shudy, A.; Kulkarni, V. V.; Phillips, A. Computational Design of Nucleic Acid Feedback Control Circuits. ACS Synth. Biol. 2014, 3, 600−616. (18) Qian, L.; Winfree, E. Scaling Up Digital Circuit Computation with DNA Strand Displacement Cascades. Science 2011, 332, 1196− 1201. (19) Qian, L.; Winfree, E.; Bruck, J. Neural network computation with DNA strand displacement cascades. Nature 2011, 475, 368−372. (20) Oishi, K.; Klavins, E. Biomolecular implementation of linear I/O systems. IET Syst. Biol. 2011, 5, 252−260. (21) Sawlekar, R.; Montefusco, F.; Kulkarni, V. V.; Bates, D. G. Implementing Nonlinear Feedback Controllers Using DNA Strand Displacement Reactions. IEEE Trans. Nanobiosci. 2016, 15, 443−454. (22) Benenson, Y.; Gil, B.; Ben-Dor, U.; Adar, R.; Shapiro, E. An autonomous molecular computer for logical control of gene expression. Nature 2004, 429, 423−429. (23) Sarpeshkar, R. Analog versus digital: Extrapolating from electronics to neurobiology. Neural Comput. 1998, 10, 1601−1638. (24) Thachuk, C.; Winfree, E.; Soloveichik, D. Leakless DNA Strand Displacement Systems. DNA Computing and Molecular Programming, Proceedings of the 21st International Conference, DNA 21, 2015; pp 133−153.

4160

DOI: 10.1021/acsomega.7b00572 ACS Omega 2017, 2, 4143−4160