<?xml-model href='http://www.tei-c.org/release/xml/tei/custom/schema/relaxng/tei_all.rng' schematypens='http://relaxng.org/ns/structure/1.0'?><TEI xmlns="http://www.tei-c.org/ns/1.0">
	<teiHeader>
		<fileDesc>
			<titleStmt><title level='a'>Simplifying Chemical Reaction Network Implementations with Two-Stranded DNA Building Blocks</title></titleStmt>
			<publicationStmt>
				<publisher></publisher>
				<date>09/04/2020</date>
			</publicationStmt>
			<sourceDesc>
				<bibl> 
					<idno type="par_id">10198586</idno>
					<idno type="doi">10.4230/LIPIcs.DNA.2020.2</idno>
					<title level='j'>Leibniz international proceedings in informatics</title>
<idno>1868-8969</idno>
<biblScope unit="volume">174</biblScope>
<biblScope unit="issue"></biblScope>					

					<author>Robert F Johnson</author><author>Lulu Qian</author>
				</bibl>
			</sourceDesc>
		</fileDesc>
		<profileDesc>
			<abstract><ab><![CDATA[In molecular programming, the Chemical Reaction Network model is often used to describe real or hypothetical systems. Often, an interesting computational task can be done with a known hypothetical Chemical Reaction Network, but often such networks have no known physical implementation. One of the important breakthroughs in the field was that any Chemical Reaction Network can be physically implemented, approximately, using DNA strand displacement mechanisms. This allows us to treat the Chemical Reaction Network model as a programming language and the implementation schemes as its compiler. This also suggests that it would be useful to optimize the result of such a compilation, and in general to find effective ways to design better DNA strand displacement systems. We discuss DNA strand displacement systems in terms of "motifs", short sequences of elementary DNA strand displacement reactions. We argue that describing such motifs in terms of their inputs and outputs, then building larger systems out of the abstracted motifs, can be an efficient way of designing DNA strand displacement systems. We discuss four previously studied motifs in this abstracted way, and present a new motif based on cooperative 4-way strand exchange. We then show how Chemical Reaction Network implementations can be built out of abstracted motifs, discussing existing implementations as well as presenting two new implementations based on 4-way strand exchange, one of which uses the new cooperative motif. The new implementations both have two desirable properties not found in existing implementations, namely both use only at most 2-stranded DNA complexes for signal and fuel complexes and both are physically reversible. There are reasons to believe that those properties may make them more robust and energy-efficient, but at the expense of using more fuel complexes than existing implementation schemes.]]></ab></abstract>
		</profileDesc>
	</teiHeader>
	<text><body xmlns="http://www.tei-c.org/ns/1.0" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xmlns:xlink="http://www.w3.org/1999/xlink">
<div xmlns="http://www.tei-c.org/ns/1.0"><head n="1">Introduction</head><p>What does it mean to optimize a molecular system? One particular field in molecular programming is currently faced with that question. The Chemical Reaction Network (CRN) model is often used to describe systems of interacting molecules. The model can either describe real systems, to analyze their behavior and computational function, or describe hypothetical systems, with known computational function but perhaps no known physical</p><p>In terms of robust DSD systems and their fuels, we can take a lesson from experiments with seesaw gates <ref type="bibr">[28,</ref><ref type="bibr">40]</ref>. For a two-reactant two-product reaction, the Soloveichik et al. translation scheme uses 3-stranded fuels <ref type="bibr">[34]</ref>, the Cardelli scheme 4-stranded fuels <ref type="bibr">[4]</ref>, and the Qian et al. scheme <ref type="bibr">[27]</ref> (in the corrected version) a 5-stranded or a 7-stranded fuel. The seesaw gates compute logic gates which are less complex than chemical reactions, but they do so with only single strands and 2-stranded complexes <ref type="bibr">[28]</ref>. Possibly because of this, they have been used to build larger circuits and to be robust to experimental imperfections, such as unpurified strands <ref type="bibr">[40]</ref>.</p><p>For this purpose, we have been investigating implementing CRNs using only 2-stranded fuels. Simple DSD systems, such as detecting a desired sequence <ref type="bibr">[5]</ref> or AND gates <ref type="bibr">[16]</ref>, are often 2-stranded, in addition to the seesaw gates mentioned above. There is even a class of hairpin-based systems that construct larger structures from single-stranded initial complexes <ref type="bibr">[44]</ref>, including the Hybridization Chain Reaction often used in imaging <ref type="bibr">[11]</ref>, and a design for hairpin-based logic circuits <ref type="bibr">[12]</ref>. However, none of these are a full Chemical Reaction Network implementation, or even an equivalently powerful dynamical system -while logic gates are universal for computing functions, CRNs have a dynamical behavior that logic gates in general do not.</p><p>We focus in this work on DSD systems using only 2-stranded fuels and where all mechanisms are physically reversible. We focus on 2-stranded fuels for the robustness concerns above, as well as the theoretical question of whether 2-stranded complexes are sufficient for complex behavior (as discussed further in <ref type="bibr">[18]</ref>). We focus on physical reversiblility because it reduces the quantity of fuel consumed by reversible reactions. Many interesting computations and dynamical behaviors require reversible reactions. For example, logically reversible operations allow computation with arbitrarily low energy if they are implemented with physically reversible reactions <ref type="bibr">[2,</ref><ref type="bibr">3]</ref>, such as DSD implementations of stack machines <ref type="bibr">[27]</ref>, Gray code counters <ref type="bibr">[9]</ref>, and space-bounded computations <ref type="bibr">[37]</ref>. DNA buffers <ref type="bibr">[29]</ref> use reversible reactions to maintain stable <ref type="bibr">[30]</ref> and dynamical <ref type="bibr">[31]</ref> spatial patterns. DNA circuits can be reset to process new input signals when reversible reactions are used for restoring fuel molecules in response to reset signals <ref type="bibr">[14,</ref><ref type="bibr">13,</ref><ref type="bibr">12]</ref>. (Existing implementations often are or can be made physically reversible; Qian et al. <ref type="bibr">[27]</ref> demonstrate it explicitly, while simple methods to make other existing schemes <ref type="bibr">[34,</ref><ref type="bibr">4]</ref> physically reversible is an exercise for the interested reader.)</p><p>In this work, we discuss ways of implementing CRNs using only 2-stranded fuels and where all mechanisms are physically reversible. We discuss four known 2-stranded DSD motifs that can serve as building blocks for such implementations, and we present a new cooperative 4-way strand exchange motif that starts with 2-stranded complexes. We discuss two ways of implementing general CRNs with these motifs, and tradeoffs between the two schemes. Finally, we show how, using CRN bisimulation, these schemes can be proven correct assuming the assumptions of the formal DSD model reflect real DSD systems.</p><p>We believe that having abstract descriptions of simple motifs will help the design of complex DSD systems. Whatever complex behavior is desired, it may be easier to implement by combining the simple logical operations of known motifs. To demonstrate this, we first discuss the 5 motifs and their behavior on an abstract level, then show how various CRN implementations can be constructed and comprehended by combining those abstract behaviors.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head n="2">Two-stranded motifs</head><p>We identify five "motifs", or simple condensed reactions, out of which we build two-stranded CRN implementations. Four of these motifs have been previously studied, while one is new. We discuss the properties of each motif in itself, while in Section 3 we will discuss how</p><p>2:4</p><p>Simplifying CRN Implementations with Two-Stranded DNA Building Blocks those properties interact when building larger circuits. For building two-stranded CRNs, key questions about a given motif are what logical operation it represents, whether its outputs have the form of its inputs and/or the inputs of the other motifs, and whether its outputs and reverse gates are 2-stranded.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head>Toehold Exchange</head><p>A reversible 3-way strand displacement exchanges which of two strands is bound to a gate (Figure <ref type="figure">1</ref> (a)).  <ref type="bibr">[34]</ref>, seesaw gates <ref type="bibr">[28]</ref>, and various others <ref type="bibr">[47]</ref>.</p><p>Cooperative Hybridization (symmetric)  <ref type="bibr">[8]</ref>. This mechanism, like the two other cooperative motifs, is "cooperative" in the sense that it requires two inputs to simultaneously, "cooperatively", displace parts of the gate complexes for a productive reaction to happen.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head>Cooperative Hybridization (asymmetric)</head><p>Two 3-way strand displacement reactions occur simultaneously on either side of a gate complex, meeting in the middle and releasing an output strand only if both inputs are present (Figure <ref type="figure">1</ref> (c)). The input strands are unbound toehold-long domain combinations, while the output strand has those two long domains in combination with a different toehold; but with only one toehold, barring complex mechanisms either one but only one of them can react. However, even if both inputs are single strands the reverse gate is a 3-stranded complex, so this motif is not "reversible with 2-stranded fuels". Introduced and tested by Zhang <ref type="bibr">[46]</ref>. Experimentally tested by Dabby <ref type="bibr">[10]</ref>. Various mechanisms, simple and complex, based on 4-way strand exchange have been used experimentally in a number of devices <ref type="bibr">[41,</ref><ref type="bibr">24,</ref><ref type="bibr">5,</ref><ref type="bibr">16]</ref>.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head>4-way Strand Exchange</head></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head>4-way Cooperative Hybridization</head><p>Two 4-way branch migrations happen on either side of a gate, meeting in the middle and separating into two intermediate complexes (Figure <ref type="figure">2</ref>). Observe that the "top" toeholds (t and t) on the initial X and Y complexes end up on one of the two products, while the "bottom" toeholds (s * and v * ) end up on another. That is, each of the two products carries only half of the information of the original reactants, and products of different instances of this reaction can interact in the reverse reaction. If for example the (t, t) top half of this reaction interacted with a (v * , s * ) bottom half from a different instance, while the (s * , v * ) bottom half interacted with an (a, a) top half, the result would be X and Y complexes with the same form as the original reactants but different toehold combinations. The effect of such a quadruplet of reactions is strand exchange between one pair of complexes coupled to strand exchange between the other, simultaneously changing the open toehold combinations on distinct long domains. This is important because affecting distinct long domains in a coupled manner was the one thing that, under a set of additional restrictions that this mechanism satisfies, our previous work <ref type="bibr">[18]</ref> showed that uncooperative 4-way strand exchange could not do.</p><p>While the other four mechanisms discussed have been experimentally demonstrated to work, cooperative 4-way branch migration has not yet been tested. In particular, the final dissociation step requires 3 toeholds separated by two 4-way junctions to dissociate. We think this is plausible, based on Dabby's observation that 2 toeholds separated by one 4-way junction can dissociate <ref type="bibr">[10]</ref>; or, if this is not the case, that there is some 0 &lt; Length(l) &#8804; 6 for which that dissociation is possible and reversible. It is possible that Length(l) = 0 (i.e. no third toehold) will give the desired behavior, but from Dabby's results, "closed" (both toehold lengths at least 2) 4-way branch migration seems to proceed much faster than "open" 4-way branch migration. Thus we suspect that Length(l) = 0, and in particular Length(l) &#8805; 2, will give the desired fast and reversible reaction kinetics.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head>An abstraction for 4-way-based mechanisms</head><p>Common to both uncooperative and cooperative 4-way strand exchange is a basic signal complex: two strands, one long domain bound to its complement flanked by one bound pair of complementary toeholds and one open pair of non-complementary toeholds, as seen repeatedly in Figures <ref type="figure">1 (d</ref>) and 2. As both types of 4-way strand exchange transform complexes of this form into complexes of the same form with different domain combinations, we find an abstract description of this type of molecule useful. For example, we write the molecule with long domain X, open 3' (end of the DNA) toehold t, open 5' toehold s * , and bound toehold m as X(t, s; m). Note that the semicolon distinguishes open toeholds t, s * available for interaction from the closed (m, m * ) toehold pair that cannot interact with other complexes, but can be opened for interaction by a reaction. When the long domain is unimportant or universal, such as a system composed entirely of uncooperative 4-way strand exchange, we omit it and write simply (t, s; m). For experimental reasons we prefer to have strands made up of only non- * or only * domains, and design non- * and * domains to have distinct sequence properties (for example, using a three-letter code <ref type="bibr">[28]</ref>). Then X(t, s; m) unambiguously describes the top reactant of Figure <ref type="figure">1 (d)</ref>, with s understood to mean an open s * toehold. With that assumption, the top product in Figure <ref type="figure">1</ref> (d) would be X(m, n; s) * , with the first toehold listed still being on the 3' end of its strand, but now understood to mean an open m * toehold. Without that assumption, we might use a more general notation where those molecules are X(t, s * ; m) and X * (m * , n; s * ) respectively. The circle abstraction shown in said figures is also useful to illustrate strand exchange reactions. Each circle represents a strand with one long domain and two toeholds, where half-faded circles represent strands made of * domains. Thin connections (both figures) represent strands bonded directly, requiring matching domains; thick connections labelled with a toehold domain (horizontal in Figure <ref type="figure">2</ref>) represent strands connected by gate strands from a cooperative 4-way strand exchange reaction, which can be between any domains so long as the appropriate gate exists. Figure <ref type="figure">2</ref> A cooperative 4-way branch migration mechanism. Initial X and Y complexes combine with a gate that matches their open toehold combinations, producing two 3-stranded complexes each with one of the strands of X and one of the strands of Y . These complexes can recombine with each other or with the corresponding products of a similar reaction, which in the latter case will produce X and Y complexes with different toehold combinations. On the right, this reaction is shown in abstracted form. The cooperative 4-way CRN is based on groups of four of these reactions, two in the reverse of the direction shown, where in the reverse reactions each product of one forward reaction interacts with the corresponding product of the other forward reaction. Complexes are labeled with names in the abstract notation if applicable, and their role in the cooperative 4-way CRN implementation scheme. "Signal" and "fuel" complexes have 2 strands as desired; stable "intermediate" complexes can have any number of strands; and "transient" complexes will quickly decay to one side or the other of the reaction. The marking of X(t, s, m) as signal and Y (v, t, n) * as fuel is based on the CRN implementation scheme presented in Section 3, but in general the two can be any combination of signal and fuel, or could be intermediates of a more complex pathway.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head>D N</head><p>In Figure <ref type="figure">2</ref> we introduce a similar notation for the "intermediate" products of a cooperative 4-way strand exchange reaction, in that case (m, l * , n; X : t, Y * : t * ) and (n * , l, m * ; Y * : v * , X : s). Again the semicolon distinguishes the three open toeholds, listed from 5' to 3' end, from the bound long domain-toehold pairs; each of those pairs is listed as the domains that appear first in 5' to 3' order. Thus the full reaction is</p><p>assuming the appropriate fuel (top center), which we do not give a notation to and omit from the reaction, is present.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head n="3">Chemical Reaction Network implementations</head><p>The above motifs can be combined in various ways to construct implementations of arbitrary Chemical Reaction Networks. To implement arbitrary CRNs, the reaction <ref type="bibr">[26]</ref>. From a logical perspective, "join" and "fork" operations are sufficient; the above reactions represent those logics.</p><p>We take modular CRN bisimulation <ref type="bibr">[19]</ref> as the definition of a "correct" CRN implementation scheme. Given that a scheme is correct, there are a number of other conditions that would be useful to satisfy for various reasons, theoretical and practical. CRN implementations typically have signal complexes that are the primary form of a given formal species, and fuel complexes that are assumed to be always present and drive the reactions. For a CRN to have "only 2-stranded inputs", as desired in this work, means that all signal complexes and fuel complexes are single strands or 2-stranded. We implicitly assume that we are discussing systematic CRN implementations, where we give a template for a generic reaction and construct larger CRNs by combining independent copies of the template with different domain identities. In such a case we can ask how the number of toehold domains scales, i.e. whether different reactions can use the same toeholds or have to create new ones; as toeholds are limited in length by thermodynamics, a system with O(n) toeholds may be able to implement small CRNs but a system with O(1) toeholds is better if possible. Whether a scheme requires cooperative mechanisms is worth noting. Finally, it is desirable for reversible reactions (A + B C + D) to be implemented with physically reversible mechanisms, so that going forward and backward multiple times does not consume fuel; to be truly reversible, the 2-stranded fuel criterion should include the reverse fuels as well. For further discussion and formal definitions of these criteria, see <ref type="bibr">[18]</ref>, which also contains a proof that no CRN implementation scheme using only 4-way branch migration can satisfy all of them.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head>Toehold Exchange-based CRNs</head><p>Existing CRN implementations <ref type="bibr">[34,</ref><ref type="bibr">27,</ref><ref type="bibr">4]</ref> are often based on toehold exchange mechanisms where e.g. A + B &#8594; C is implemented by a toehold exchange reaction with A opening a toehold on the gate for a reaction involving B. These schemes can be understood in light of the motifs previously discussed: the property of toehold exchange that a different toehold on the gate is opened allows join and fork logic. The property that the released strand has a different long domain/toehold combination is used to pass signals between gates. The same shared-toehold logic could also be used with 4-way branch migration instead of toehold exchange, similar to the 4-way-based AND gate <ref type="bibr">[16]</ref> (although that gate itself uses a toehold hidden in a loop rather than a toehold shared between adjacent long domains, which is a line of investigation to be explored elsewhere).</p><p>Such a shared-toehold mechanism seems to require a 3-stranded complex for the gate molecule to achieve join logic, so it does not meet the goal of this paper, but is worth mentioning as the current state of the art. Another relevant mechanism using toehold exchange is the seesaw gate <ref type="bibr">[28]</ref>, where transduction logic combines with threshold logic to check whether the total amount of signal is more than either A or B can produce by itself. This achieves join logic for macroscopic signals but cannot satisfy criteria such as CRN bisimulation for individual molecules.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head>3-way Cooperative CRNs</head><p>The symmetric cooperative hybridization is A + B C + D logic, if we consider the same long domain in a different toehold context to be a different signal. Since toehold exchange reactions depend on the combination of long domain and toehold, this is valid. Thachuk et al. use a combination of symmetric cooperative hybridization and toehold exchange to implement leakless A + B &#8594; C + D reactions in exactly this manner <ref type="bibr">[38,</ref><ref type="bibr">39,</ref><ref type="bibr">42]</ref>).</p><p>From our perspective, the only problem is that symmetric cooperative hybridization with 1-stranded inputs produces 2-stranded products, and toehold exchange with a 2-stranded input signal produces a 3-stranded reverse gate. For physically reversible reactions, this  3-stranded gate would be considered a reverse fuel, and the system would not be made with entirely 2-stranded fuels. Thus this mechanism meets all our criteria for irreversible CRNs, but not reversible CRNs.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head>4-way-based CRNs with O(n) toeholds</head><p>The two-toehold-mediated 4-way strand exchange mechanism effectively exchanges toeholds on a common long domain; note that while the inputs both have t and s toeholds, the outputs have one with only t and one with only s. When a signal complex goes through multiple copies of this reaction with different fuels, it can turn any combination of toeholds into any other combination. When two signals with complementary pairs of toeholds meet in this reaction, it produces two signals with different combinations in A + B C + D logic. So for example, we can turn (a 1 , a 2 ; a 3 ) into (r 1 , r 2 ; r 3 ) and (b 1 , b 2 ; b 3 ) into (r 2 , r 1 ; r 4 ), which will react and produce (r 3 , r 4 ; r 2 ) * and (r 4 , r 3 ; r 1 ) * , which can be turned into (c 1 , c 2 ; c 3 ) and (d 1 , d 2 ; d 3 ) respectively. Thus two-toehold-mediated 4-way strand exchange alone can implement arbitrary reversible CRNs if we allow O(n) toeholds.</p><p>A list of all species involved is given in Table <ref type="table">1</ref>. Note that fuels (r 2 , r 1 ; r 5 ) and (r 1 , r 2 ; r 6 ) can interact, but the products can do nothing but reverse the reaction, and the same is true for (r 4 , r 3 ; c 2 ) * with (r 3 , r 4 ; d 2 ) * .</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head>4-way Cooperative CRNs</head><p>The cooperative 4-way strand exchange motif in Figure <ref type="figure">2</ref>, when its products recombine with products of a different instance of the reaction, simultaneously exchanges the toehold combinations on a complex with long domain X and a complex with long domain Y . If A(t, s; m) is the signal molecule for A, then simultaneously breaking the (t, s) combination on A and putting together a (u, v) combination on some long domain R is effectively converting A(t, s; m) R(v, u; n) * if all other molecules involved are considered fuels. Where R is unique to the reaction A + B C + D, we can convert the four signal species from their own long domains to the R domain, then use a two-toehold-mediated 4-way strand exchange reaction to implement the reaction itself. In contrast to the previous implementation scheme, that each reaction has a different long domain allows the toeholds (u, v, etc.) to be universal, using O(1) toeholds at the expense of requiring cooperative hybridization. In the notation used in Figure <ref type="figure">2</ref>, this quadruplet of reactions (with the appropriate top-center fuels assumed present but not written) is As this scheme is based on the O(n)-toehold scheme, we reuse the mechanism from Table <ref type="table">1</ref>. Assume all complexes in that list have long domain R, unique to the reaction A + B C + D. To the toeholds listed, add toeholds t, s, m, n, l, and let a 3 = b 3 = c 3 = d 3 = n * , with u and v in the above quadruplet renamed appropriately. Then use cooperative 4-way strand exchange to convert A(t, s; m) (R * (a * 1 , a * 2 ; n)) * = R(a 1 , a 2 ; n * ) (the fuel will have R * on the "top" strand with A), B(t, s; m) R(b 1 , b 2 ; n * ), C(t, s; m) R(c 1 , c 2 ; n * ), and D(t, s; m) R(d 1 , d 2 ; n * ). This gives a mechanism with one long domain per species, one long domain per reaction, and a total of 19 toeholds. Because the long domains now indicate species/reaction identity, the toeholds can be shared between all species and reactions without crosstalk.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head n="4">Correctness of the schemes</head><p>The correctness of the schemes can be verified by CRN bisimulation, a formal definition of correctness of a CRN implementation that implies several desirable properties <ref type="bibr">[19]</ref>. Below we give an intuitive explanation of why the schemes are correct that parallels the definition of CRN bisimulation; readers familiar with CRN bisimulation can fill in the details of the formal proof. Intuitively, CRN bisimulation consists of interpreting each DNA complex as zero or more formal species, then confirming that the behavior of the formal system and the interpreted DSD system are the same from any initial state. That is to say, any reaction of the DNA complexes should be interpreted as a reaction of formal species that is either valid or trivial ("anything that can happen, should"), and any reaction of the formal interpretation of a set of DNA complexes should be possible, perhaps after some trivial reactions, starting from that set of DNA complexes ("anything that should happen, can"). Table <ref type="table">1</ref> is effectively a proof of the correctness of the O(n)-toehold 4-way-based scheme according to CRN bisimulation <ref type="bibr">[19]</ref>. For each A + B C + D reaction, construct a copy of this mechanism with unique r i domains, but any a i domains in common with other reactions using the same formal species; reactions with fewer reactants or products can have one of A, B, C, or D as a fuel; reactions with more reactants or products should D N A 2 6</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head>2:10</head><p>Simplifying CRN Implementations with Two-Stranded DNA Building Blocks be broken into steps with at most 2 of each <ref type="bibr">[26]</ref>. DNA complexes in columns labeled A, B, C, or D are interpreted as one copy of the corresponding species, while complexes in columns labeled &#8709; are fuels. Formally, fuels are assumed always present and removed from the enumerated implementation CRN before bisimulation verification; so for example the physical pathway (r 2 , a 2 ; r 3 ) + (a 2 , r 2 ; r 5 ) (r 5 , r 3 ; r 2 ) * + (r 3 , r 5 ; a 2 ) * would be represented as (r 2 , a 2 ; r 3 ) (r 5 , r 3 ; r 2 ) * , and then interpreted as the trivial reaction A A. Using the abstraction for 4-way strand exchange notation, the table is structured such that each non-fuel species can interact with the (usually two) fuel species in the same row, producing the corresponding fuel+non-fuel pair above or below it; that the final A + B forms react to produce the final C + D forms, while their fuels also have a spurious-but-harmless reaction with each other; and that, given the uniqueness of the domains, no other intramodule or inter-module reactions exist. In CRN bisimulation, we say that a reaction interpreted as, for example, A</p><p>A is "trivial", and in this case all reactions are trivial except (r 1 , r 2 ; r 3 ) + (r 2 , r 1 ; r 4 ) (r 3 , r 4 ; r 2 ) * + (r 4 , r 3 ; r 1 ) * which is interpreted as the desired reaction A + B C + D. With (a 1 , a 2 ; a 3 ) etc. as the signal species, one can see that the signal species can implement the formal reaction, and any intermediate species can turn into the common species with the same interpretation by interacting with only fuels. Intuitively this is a good argument for correctness, and readers familiar with CRN bisimulation will recognize the above as a sufficient condition for modular CRN bisimulation with respect to the signal species as common species.</p><p>For the cooperative 4-way scheme, the same bisimulation logic applies. In the notation used in Figure <ref type="figure">2</ref> and Section 3, in e.g. A(t, s; m) R(a 1 , a 2 , n * ) the signal complex A(t, s; m), output complex R(a 1 , a 2 , n * ), and intermediate (m, l * , n; A : t, R : a 2 ) all interpreted as A, while the other three intermediates and all the fuels will each be interpreted as nothing. From there the bisimulation proof follows the O(n)-toeholds case. In this case the lack of crosstalk between modules is assured by the distinct long domains; even if toehold combinations are identical, different long domains will make the reaction unproductive. The remaining caveat is with the cooperative 4-way mechanism itself. We designed the system so that the toeholds along the cooperative reaction are always m, l, n. Thus, we assume that intermediates of the cooperative pathway will all have the matching m, l, n toeholds, and all three toeholds will bind and dissociate as a unit. Whether this is actually true or not will be determined experimentally; if not, there may be problematic crosstalk between, for example, an (A, R 1 ) and (A, R 2 ) pair of long domains which leads to temporarily duplicated signals. If it is true, however, then the result of such a crosstalk will be a release of one side with the other suspended, one of which carries the signal, and the system will be correct according to bisimulation.</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head n="5">Discussion</head><p>We discussed the use of DNA Strand Displacement to implement Chemical Reaction Networks, and the desire to create larger, more robust DSD CRN implementations. We then presented 2-stranded DSD motifs which we used to build 2-stranded CRN implementations, in the hope that they would be more robust than those which rely on 3-or-more-stranded complexes.</p><p>There is some indication that 2-stranded DSD systems in general are more robust (as we briefly reviewed in the introduction), but whether these particular systems are more robust than the current state-of-the-art CRN implementations is an open question. We can compare Soloveichik et al.'s original CRN scheme <ref type="bibr">[34,</ref><ref type="bibr">36]</ref> (which is reasonably representative of other toehold exchange schemes), our O(n)-toehold 4-way strand exchange scheme, and our (O(1)-toehold) cooperative 4-way strand exchange scheme. While 3-and 4-stranded complexes may be less robust, in other aspects the toehold exchange scheme is simpler than our two schemes: it uses one long domain per formal species, one long domain per reaction, and can be done with a single, universal toehold. To go from reactant signal species to product signal species in the toehold exchange scheme (as implemented experimentally <ref type="bibr">[36]</ref>) takes 4 toehold exchange steps in an A + B &#8594; C + D reaction, and generalizes naturally to n + m steps in an n-reactant m-product reaction. In contrast, while the cooperative 4-way scheme also uses one long domain per formal species and reaction, as described above it uses 19 universal toeholds and takes 30 reactions for A + B &#8594; C + D. (By "reaction" we mean roughly one condensed reaction as described in Peppercorn <ref type="bibr">[17]</ref>, generalized to include trimolecular reactions. So one toehold exchange or one 2-toehold-mediated 4-way strand exchange is one reaction, as is the cooperative 4-way strand exchange shown in Figure <ref type="figure">2</ref>; note that using that mechanism to exchange e.g. A(t, s; m) R(a 1 , a 2 ; n * ) takes 4 such reactions.) The O(n)-toeholds scheme takes only 14 reactions for A + B &#8594; C + D, but with one universal long domain it takes 3 toeholds per species and 6 per reaction, which may run out of design space for large CRNs. Also, 14 reactions is still much more than 4. These pathways are not provably optimal; we suspect they can be reduced to less than 14 and 30, but still more than 4.</p><p>The increase in number of reactions to implement A + B &#8594; C + D may just be a cost of using 2-stranded complexes. The fundamental question is, given a complex of a certain size, how much information can it store? How can complexes meant to represent A, C, and an E from another reaction all present different enough open and bound domains that none can undergo a reaction meant for a different one? With 3-stranded complexes and toehold exchange, the long domain identity and open toehold does this very efficiently. With 2-stranded complexes and 4-way strand exchange, we use pairs of toehold identity to represent signal identity, which means we need extra reactions to (a) change the toehold identity one strand at a time, and (b) ensure that intermediates of different pathways don't try to pass through the same toehold combination. This question, then, connects to another work of ours. The final result of that work was a proof that a systematic CRN implementation that satisfies certain desirable conditions, including using only 2-stranded inputs and the other conditions discussed at the beginning of Section 3, cannot be done with DSD using only 4-way branch migration <ref type="bibr">[18]</ref>. The steps taken to prove that result involve questions of what sort of transformations are possible with DSD reactions, and how and whether the possibility of certain transformations can depend on the features of the strands. This "dependence" is in the sense that the release of a strand in toehold exchange "depends on" the incoming strand having the correct toehold and long domain identities, or the way we have to structure our CRN implementations so that production of the output species depends on the inputs having the correct toehold identity pairs. Thus, further exploration of that line of investigation might help answer some of the questions suggested by the mechanisms in this paper, of whether 2-stranded complex based CRN implementations inherently require longer pathways, and quantitatively how much longer. Moreover, the investigation could be expanded to include other CRN implementations involving enzymes. For example, transcriptional circuits <ref type="bibr">[20,</ref><ref type="bibr">21]</ref>, PEN-DNA toolbox <ref type="bibr">[23,</ref><ref type="bibr">1]</ref>, primer exchange reaction cascades <ref type="bibr">[22]</ref>, and strand-displacing polymerase systems <ref type="bibr">[35,</ref><ref type="bibr">32,</ref><ref type="bibr">33]</ref> all have elementary reactions that can be abstracted as motifs and are candidates for formal analysis. In these systems, it is possible to start with fewer and simpler fuel molecules (e.g. single strands only) while more complex molecules can be generated by DNA polymerase to carry out desired reactions. In addition to 3-way and 4-way strand displacement with standard toeholds, other mechanisms could also be investigated, D N A 2 6</p></div>
<div xmlns="http://www.tei-c.org/ns/1.0"><head>2:12</head><p>Simplifying CRN Implementations with Two-Stranded DNA Building Blocks including remote <ref type="bibr">[15]</ref>, associative <ref type="bibr">[6]</ref>, and allosteric <ref type="bibr">[43]</ref> toeholds. These mechanisms may allow further simplification of the implementations as they enrich the design space with alternative representations of signals.</p><p>It is also worth discussing how we discovered the cooperative 4-way strand exchange motif and associated CRN implementation in the process of working out the impossibility proof in <ref type="bibr">[18]</ref>. We give an intuitive list of those conditions at the beginning of Section 3, but readers desiring a formal list of conditions should see <ref type="bibr">[18]</ref>. Two of the conditions are using only O(1) toeholds and not using cooperative mechanisms, so both the O(n) toeholds uncooperative 4-way strand exchange based scheme and the O(1) toeholds cooperative 4-way strand exchange based scheme satisfy all but one of the conditions, each failing to satisfy a different one. Thus in some sense this paper is the positive counterpart to the previous negative result, forming a tight upper and lower bound on the complexity of DSD implementations of CRNs. But this pair of results also has implications for design of DSD systems. The cooperative 4-way strand exchange motif and the process by which we came up with it is potentially a proof of concept that, in systematically eliminating possibilities in DSD systems, we can find new motifs in whatever remains. How exactly this can be generalized we do not know, but if it can be, it may make the process of designing DSD systems faster and more systematic.</p><p>Another aspect worth mentioning is the focus on motifs before building up CRN implementations. We argued that each of the 5 motifs has certain abstract behaviors, and that larger systems such as CRN implementations can be thought of in terms of those behaviors. When building large systems, it is much easier if one can build mid-sized building blocks out of the fundamental units, then build larger systems out of the mid-sized building blocks. Motifs take that role between fundamental DSD steps (bind, unbind, 3-way branch migration, 4-way branch migration) and systems on the scale of CRN implementations. To the extent that we were able to describe our CRN implementations in terms of the motifs rather than in terms of the underlying DSD steps, this approach should be considered for future DSD system design.</p></div></body>
		</text>
</TEI>
