This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

Self-Assembling DNA Complexes with a Wheel Graph Structure

Gabriel Lopez California State University, San Bernardino, San Bernardino, CA gabriel.lopez.x15@gmail.com Cory Johnson corrine.johnson@csusb.edu
Abstract

The Watson-Crick complementary properties of DNA make DNA a useful tool for the self-assembly of various target complexes. Concepts from graph theory can be used to model the self-assembling process in which the vertices of the graph represent kk-armed branched junction molecules, called tiles. We seek to determine the minimum number of tile and cohesive-end types necessary to create the desired self-assembled complex. Although results are known for a few infinite classes of graphs, many classes of graphs remain unsolved. We present results for the wheel graph within the restrictions of three different settings.

keywords:
graph theory , DNA self-assembly , flexible tile model , wheel graph
journal:

1 Introduction

1.1 Self Assembling DNA Complexes

Advancements in nanotechnology have allowed biologists to construct various complexes from genetic material. This is accomplished by utilizing the complementary properties of DNA molecules so that the molecules bond in a specific way. These complexes, while proven to be incredibly useful, can also be costly to synthesize in a lab. One of our objectives is to make the creation of these structures more efficient.

There are several applications of DNA self-assembly including targeted drug delivery as a treatment for various types of cancer, gene therapy, the formation of fine screen filters, biomolecular computing, and biosensors (see [1, 3, 4, 5, 6]). Our focus has been on complexes formed from flexible kk-armed branched junction DNA molecules. While some general results are known, there are many open problems in finding the most optimal construction. See [2] for overview of problem and summary of main results. The notion of optimizing the construction of these complexes comes from the fact that creating a complex will require a certain number of bond-edge types, as well as a certain number of tile types. We use concepts from graph theory to model this design process. We then optimize the construction of a target graph by providing an algorithm that will construct the graph with as few tile and bond-edge types as possible under the constraints of one of the three following scenarios:

  • 1.

    Scenario 1: A pot of tile types may realize other graphs of smaller order than the target graph (that is, a graph on fewer vertices) or nonisomorphic graphs of the same size as the target graph.

  • 2.

    Scenario 2: A pot of tile types may realize graphs which are of the same order as the target graph, but are nonisomorphic. Graphs of smaller order than the target graph cannot be realized by the pot.

  • 3.

    Scenario 3: A pot of tile types may not realize graphs of smaller order nor any nonisomorphic graphs of the same order as the target graph.

We introduce the following notation for the minimum number of tile types and bond-edge types needed to construct a complex.

Definition.

Let GG be a target graph. We denote the minimum number of bond-edge types needed to construct GG in Scenario ii by Bi(G),B_{i}(G), and the minimum number of tile types needed to construct GG in Scenario ii by Ti(G).T_{i}(G).

1.2 Concepts from Graph Theory

A DNA complex can be modeled using a graph, a mathematical object consisting of a nonempty set of nodes, called vertices, and a set of edges which connect pairs of distinct vertices.

Take as example the DNA complex in Figure 1 which takes the shape of a cube. This structure can be modeled by the graph in Figure 2.

Refer to caption
Figure 1: A DNA complex.
Refer to caption
Figure 2: The target graph associated with the complex in Figure 1.

The structure of a DNA complex becomes more easily studied as a graph rather than a complex made from genetic material. The kk-armed DNA molecules that comprise the complex are modeled by objects called tiles, which consist of a single vertex incident to kk half-edges. Each half-edge represents one arm of the molecule, with some paired sequence of DNA bases at the end. Since the end of each branch is unpaired, then in a lab environment each branch of the molecule seeks its complementary bond type to complete the sequence. We refer to a collection of tiles used to construct a graph as a pot of tiles that realizes the target graph. See an example of a branched molecule along with its corresponding tile in Figure 3.

Refer to caption
Refer to caption
Figure 3: A 3-armed branched junction molecule (left) with example tile representation (right)

Half edges are determined by their respective bond-edge types, which we denote with some letter, with hatted and unhatted labels to denote pairs of complementary bond-edge types. For example, a half edge labeled aa must bond with a half edge labeled a^\hat{a}. In certain situations, it is useful to think of the half edges of each tile bonding together to form directed edges. For the purpose of this paper, the edge will be directed from the vertex with the unhatted half edge to the vertex with the hatted half edge.

Constructed using tiles rather than a plain graph, the cube complex in Figure 1 could be represented by the graph on the left in Figure 4. As a directed graph, the same complex may be constructed as seen on the right side of Figure 4.

Refer to caption
Refer to caption
Figure 4: A cube structure built with tiles (left), with equivalent digraph representation (right).

The process of modeling a self-assembling DNA complex with a (directed) graph is at the core of the work presented here. We seek to find a labeling of a target graph with a minimum alphabet of bond-edge types and a minimum set of distinct tiles.

Shifting our focus from complexes to graphs, it is useful to have a few concepts from graph theory in mind. We use the following graph theory concepts throughout this paper.

Definition.

Let GG be a graph. We denote the number of different vertex degrees in GG by av(G)av(G), the number of distinct even vertex degrees in GG by ev(G)ev(G), and the number of distinct odd vertex degrees by ov(G).ov(G).

Definition.

A cycle in a graph GG that contains every vertex of GG is called a Hamilton cycle. Any graph which possesses a Hamilton cycle is called Hamiltonian.

1.3 The Construction Matrix

We are able to use tools from linear algebra to aid us in verifying that a given pot of tiles has the desired properties. Specifically, there is a linear system that governs the ratio of tile and bond-edge types needed when constructing a target graph. This is due to the fact that for every half edge on every tile, there must be a complementary half edge of the same bond-edge type somewhere else in the pot.

The following results are directly from [2] and are written here for later reference:

Definition.

Let P={t1,,tp}P=\{t_{1},\ldots,t_{p}\} be a pot of pp distinct tiles. We define Ai,jA_{i,j} to be the number of cohesive ends of type aia_{i} on tile tjt_{j}, and A^i,j\hat{A}_{i,j} to be the number of cohesive ends of type a^i\hat{a}_{i} on tile tjt_{j}. Let zi,jz_{i,j} be the net number of cohesive ends of type aia_{i} on tile type tjt_{j}, i.e. zi,j=Ai,jA^i,jz_{i,j}=A_{i,j}-\hat{A}_{i,j}. Define rir_{i} to be the proportion of tile type tit_{i} to be used in the assembly process.

Proposition 1.

Let P={t1,t2,tp}P=\{t_{1},t_{2},...t_{p}\} be a pot. Then the number of hatted cohesive ends of each bond type must equal the number of unhatted cohesive ends of the same type that appear in the construction of GG. That is, for all ii, jrjzi,j=0\sum_{j}r_{j}z_{i,j}=0.

The following system of equations captures the requirements outlined in Proposition 1.

z1,1r1+z1,2r2++z1,prp\displaystyle z_{1,1}r_{1}+z_{1,2}r_{2}+\dots+z_{1,p}r_{p} =0\displaystyle=0
\displaystyle\vdots
zm,1r1+zm,2r2++zm,prp\displaystyle z_{m,1}r_{1}+z_{m,2}r_{2}+\dots+z_{m,p}r_{p} =0\displaystyle=0
r1+r2++rp\displaystyle r_{1}+r_{2}+\dots+r_{p} =1\displaystyle=1

The construction matrix of PP, denoted M(P)M(P), is the corresponding augmented matrix:

[z1,1z1,20zm,1zm,20111].\begin{bmatrix}z_{1,1}&z_{1,2}&\dots&0\\ \vdots&\vdots&\ddots&\vdots\\ z_{m,1}&z_{m,2}&\dots&0\\ 1&1&\dots&1\end{bmatrix}. (1)
Definition.

Given a pot PP, we define C(P)C(P) to be the set of graphs that can be constructed from PP. The set of graphs of minimum order that may be constructed from PP is denoted Cmin(P)C_{min}(P). We write mPm_{P} for the order of a smallest graph that may be constructed from PP.

Proposition 2.

Let P={t1,t2,,tp}.P=\{t_{1},t_{2},\dots,t_{p}\}. Then:

  1. 1.

    If a graph GG of order nn may be constructed from PP, using RjR_{j} tiles of type tjt_{j}, then (1/n)R1,R2,,Rp(1/n)\langle R_{1},R_{2},\dots,R_{p}\rangle is a solution of the construction matrix M(P).M(P).

  2. 2.

    If r1,,rp\langle r_{1},\dots,r_{p}\rangle is a solution of the construction matrix M(P)M(P), and there is a positive integer nn such that nrj0nr_{j}\in\mathbb{Z}_{\geq 0} for all jj, then there is a graph of order nn that may be constructed using nrjnr_{j} tiles of type tjt_{j}.

  3. 3.

    mp=min{lcm{bj|rj0 and rj=aj/bj},where r1,,rp is a solution to M(P)}m_{p}=\textrm{min}\{lcm\{b_{j}|r_{j}\neq 0\textrm{ and }r_{j}=a_{j}/b_{j}\},\textrm{where }\langle r_{1},\dots,r_{p}\rangle\textrm{ is a solution to }M(P)\} where the minimum is taken over all solutions to M(P)M(P) such that rj0r_{j}\geq 0 and aj/bja_{j}/b_{j} is in reduced form for all jj.

This paper describes results for target complexes that can be modeled with a wheel graph.

1.4 Wheel Graphs

A wheel graph on nn vertices, denoted WnW_{n}, is a graph consisting of a cycle on n1n-1 vertices and an additional vertex (often called the hub) which is adjacent to every vertex on the cycle.

Refer to caption
Figure 5: The wheel graph W7W_{7}.

For the purpose of this paper, we refer to the n1n-1 cycle subgraph as the outer cycle. The edges that connect the hub to the vertices on the outer cycle are called spokes.

A property of wheel graphs that we use later on in Section 3.2 is that they are Hamiltonian. One possible Hamilton cycle is shown in Figure 6.

Refer to caption
Figure 6: The Hamilton cycle in a wheel graph visits all nn vertices of WnW_{n}.

In the case of a wheel graph on 4 vertices, we note that W4W_{4} is isomorphic to the complete graph on 4 vertices, K4.K_{4}. All results presented in this paper for the case where n=4n=4 are consistent with the results for K4K_{4} in [2]. Unless otherwise stated, this paper assumes any wheel graph has order greater than 4; that is, we assume n>4n>4.

2 Methods

Previous authors have provided a number of useful theorems that we use to prove new results on the construction of wheel graphs. The following are results from [2]:

Theorem 1.

av(G)T1(G)ev(G)+2ov(G)av(G)\leq T_{1}(G)\leq ev(G)+2ov(G)

Theorem 1 tells us that in Scenario 1, the number of distinct tile types needed to construct a target graph is bounded below by the number of distinct vertex degrees that are present in the graph, and bounded above by the sum of the number of distinct even degrees and twice the number of distinct odd degrees present in the graph. The lower bound holds for T2(G)T_{2}(G) and T3(G)T_{3}(G) as well.

We also take advantage of the following preposition regarding the hierarchy of tile type and bond-edge type minima between scenarios.

Proposition 3.

B1(G)B2(G)B3(G)B_{1}(G)\leq B_{2}(G)\leq B_{3}(G) and T1(G)T2(G)T3(G)T_{1}(G)\leq T_{2}(G)\leq T_{3}(G).

The next theorem gives the minimum number of tile and bond-edge types needed to construct the cycle graph on nn vertices, CnC_{n}, in Scenario 3. This result is useful since the cycle graph is a subgraph of the wheel graph.

Proposition 4.

B3(Cn)=n2B_{3}(C_{n})=\left\lceil\frac{n}{2}\right\rceil and T3(Cn)=n2+1T_{3}(C_{n})=\left\lceil\frac{n}{2}\right\rceil+1.

The following result is another lemma that we use in the proof of one of our theorems.

Lemma 1.

In Scenario 3, if PP is a pot such that {G}=Cmin(P)\{G\}=C_{min}(P), and GG has no loops, then no tile type tPt\in P may be used for two adjacent vertices in GG.

3 Results

3.1 Wheel Graphs in Scenario 1 and Scenario 2

We have developed an algorithm to construct wheel graphs in Scenarios 1 and 2 in the most efficient way possible. We use a pot of two tile types and a single bond-edge type to construct a wheel graph of any order.

The pot we propose is given below:

P={t1,t2},P=\{t_{1},t_{2}\}, (2)

where

t1={a,a^2},and t2={an1}.t_{1}=\{a,\hat{a}^{2}\},\textrm{and }t_{2}=\{a^{n-1}\}.

This leads us to our first theorem.

Theorem 2.

Let WnW_{n} be a wheel graph on nn vertices. Then T1(Wn)=2T_{1}(W_{n})=2 and B1(Wn)=1B_{1}(W_{n})=1.

Refer to caption
Figure 7: W7W_{7} constructed using the pot of tiles in Equation (2).
Proof.

Let WnW_{n} denote a wheel graph on nn vertices. From [2], we know that B1(Wn)=1B_{1}(W_{n})=1. We observe that av(Wn)=2av(W_{n})=2, so by Theorem 1, we know that 2T1(Wn)2\leq T_{1}(W_{n}). This lower bound is achieved by the pot proposed above. We may now conclude that T1(Wn)=2T_{1}(W_{n})=2. ∎

Corollary 1.

Let WnW_{n} be a wheel graph on nn vertices. Then T2(Wn)=2T_{2}(W_{n})=2 and B2(Wn)=1B_{2}(W_{n})=1.

Proof.

By Theorem 1 and Proposition 3, 2T1(Wn)T2(Wn)2\leq T_{1}(W_{n})\leq T_{2}(W_{n}). This lower bound is achieved by the pot proposed above. We will now use the construction matrix to verify that no smaller graphs may be constructed from the same pot. The pot of tiles in (2) gives us the following construction matrix:

[1n10111],\begin{bmatrix}-1&n-1&0\\ 1&1&1\\ \end{bmatrix},

which yields the unique solution n1n,1n\langle\frac{n-1}{n},\frac{1}{n}\rangle, thus by Proposition 2 a graph on nn vertices is the smallest graph realized by this pot. Therefore, T2(Wn)=2T_{2}(W_{n})=2 and B2(Wn)=1B_{2}(W_{n})=1. ∎

It is worth noting that while other graphs on nn vertices nonisomorphic to the wheel graph are able to be constructed from this same pot of tiles, the conditions of Scenarios 1 and 2 do not prohibit this.

Refer to caption
Figure 8: Nonisomorphic graph realized by the pot in Scenario 2

3.2 Wheel Graphs in Scenario 3

Wheel graphs in Scenario 3 require a different construction algorithm compared to the construction in Scenarios 1 and 2 due to the stricter condition that no nonisomorphic graphs of the same order may be realized by a pot of tiles. To begin, we introduce two lemmas that aid us in optimizing the construction of WnW_{n}.

Lemma 2.

If PP is a pot such that {Wn}=Cmin(P)\{W_{n}\}=C_{min}(P), then no bond-edge type used in the construction of WnW_{n} may be used both on the outer cycle and a nonincident spoke.

Proof.

We proceed by contradiction. Suppose there exists a pot PP such that {Wn}=Cmin(P)\{W_{n}\}=C_{min}(P) and some bond-edge type, say aa, appears on an edge on the outer cycle as well as a nonincident spoke edge. See the graphs on the left side of Figure 9.

Refer to caption
Figure 9: Graph with bond-edge type on the outer cycle and a nonincident spoke edge

Without loss of generality, we may assume the edge on the outer cycle is oriented counterclockwise, call this edge (u,v)(u,v). The spoke edge with the same bond-edge type may be oriented either towards the hub or towards the outer cycle, but in either case, it is possible that the half edges on the spoke may reattach with the half edges on the outer cycle. This will produce a multiple edge with the hub and vertex uu or vv, as seen on the right side in Figure 9. The resulting graph is nonisomorphic to WnW_{n}, thus contradicting the assumption that {Wn}=Cmin(P)\{W_{n}\}=C_{min}(P). ∎

Lemma 3.

If PP is a pot such that {Wn}=Cmin(P)\{W_{n}\}=C_{min}(P), then a bond-edge type can be used at most twice on the outer cycle in the construction of WnW_{n}.

Proof.

By way of contradiction, suppose {Wn}=Cmin(P)\{W_{n}\}=C_{min}(P) and suppose in the construction of WnW_{n} more than two edges along the outer cycle are labeled with the same bond-edge type. By the Pigeonhole Principle, at least two of these edges must be directed in the same orientation. Without loss of generality, suppose these two edges are oriented clockwise. See Figure 10 for reference.

Refer to caption
Figure 10: Graph with bond-edge type appearing at least three times on the outer cycle

This arrangement allows for the edges to detach and rejoin as in Figure 11 to form the complex HH. We will now show that this will necessarily yield a nonisomorphic graph.

Refer to caption
Figure 11: Resulting graph HH

Recall, WnW_{n} has a Hamilton cycle, as seen in Figure 6. If graph HH is isomorphic to the original wheel graph, then it must have a Hamilton cycle.

However, refer to Figure 11 to see that the set of vertices excluding the hub, hh, may be partitioned into two sets so that no two vertices in different sets are adjacent. Suppose without loss of generality that we start a walk at vertex hh. From there, a walk of distinct vertices will include each vertex, uiu_{i} for 1ik1\leq i\leq k. The only way to include vertex v1v_{1} in a cycle is for the walk to pass through vertex hh. Thus, no Hamilton cycle exists in the graph HH. Therefore, the pot PP realizes a graph nonisomorphic to WnW_{n}, contradicting the assumption that {Wn}=Cmin(P)\{W_{n}\}=C_{min}(P). ∎

With these new lemmas, we have developed an algorithm to construct a wheel graph that builds upon the construction of the cycle graph in Scenario 3 (see Proposition 4). In particular, to construct the wheel graph WnW_{n}, we begin with the construction of the cycle graph Cn1C_{n-1}, pictured in Figure 12.

Refer to caption
Figure 12: The optimized cycle graph C6C_{6}.

From the optimized construction of the cycle graph, we can append the hub vertex and all necessary edges to the vertices on the outer cycle with one new bond-edge type.

Refer to caption
Figure 13: The proposed optimization of the wheel graph.

To prove that the construction described above uses the minimum number of tile types and bond-edge types needed to construct the wheel graph on nn vertices, it suffices to show that the proposed labeling constructs the wheel graph as efficiently as possible, i.e., that it is not possible to construct WnW_{n} with fewer bond-edge types or tile types in Scenario 3.

Theorem 3.

B3(Wn)=n2+1B_{3}(W_{n})=\left\lfloor\frac{n}{2}\right\rfloor+1.

Proof.

In the wheel graph, WnW_{n}, there are n1n-1 edges on the outer cycle. By Lemma 3, no bond-edge type on the outer cycle may be used more than twice. Thus, B3(Wn)n12=n2.B_{3}(W_{n})\geq\left\lceil\frac{n-1}{2}\right\rceil=\left\lfloor\frac{n}{2}\right\rfloor.

If any bond-edge type from the outer cycle is used on more than one spoke edge, then that bond-edge type will appear on both the outer cycle and some nonincident spoke. By Lemma 2, this is not permitted. Since there are n1n-1 spoke edges, then it is necessary to introduce at least one new bond-edge type. Thus, B3(Wn)n2+1.B_{3}(W_{n})\geq\lfloor\frac{n}{2}\rfloor+1.

The following pots of tiles realize WnW_{n} using exactly n2+1\lfloor\frac{n}{2}\rfloor+1 bond-edge types. Note that PevenP_{even} and PoddP_{odd} are the pots when nn is even or odd, respectively.

Peven={\displaystyle P_{even}=\Big{\{} t1={a1n1},t2={a^1,a22},ti={a^1,a^i1,ai},for i=3,,n2+1,\displaystyle t_{1}=\{a_{1}^{n-1}\},t_{2}=\{\hat{a}_{1},a_{2}^{2}\},t_{i}=\{\hat{a}_{1},\hat{a}_{i-1},a_{i}\},\text{for }i=3,\ldots,\left\lfloor\frac{n}{2}\right\rfloor+1, (3)
tn2+2={a^1,a^n2,a^n2+1}}\displaystyle t_{\lfloor\frac{n}{2}\rfloor+2}=\{\hat{a}_{1},\hat{a}_{\lfloor\frac{n}{2}\rfloor},\hat{a}_{\lfloor\frac{n}{2}\rfloor+1}\}\Big{\}}
Podd={\displaystyle P_{odd}=\Big{\{} t1={a1n1},t2={a^1,a22},ti={a^1,a^i1,ai},for i=3,,n2+1,\displaystyle t_{1}=\{a_{1}^{n-1}\},t_{2}=\{\hat{a}_{1},a_{2}^{2}\},t_{i}=\{\hat{a}_{1},\hat{a}_{i-1},a_{i}\},\text{for }i=3,\ldots,\left\lfloor\frac{n}{2}\right\rfloor+1, (4)
tn2+2={a^1,a^n2+12}}\displaystyle t_{\lfloor\frac{n}{2}\rfloor+2}=\{\hat{a}_{1},\hat{a}_{\lfloor\frac{n}{2}\rfloor+1}^{2}\}\Big{\}}

The pot of tiles where nn is even produces the following construction matrix:

M(Peven)=[(n1)1111110021000000011000000011000000010000000011000000110111111111]M(P_{even})=\begin{bmatrix}(n-1)&-1&-1&-1&-1&\cdots&-1&-1&0\\ 0&2&-1&0&0&\cdots&0&0&0\\ 0&0&1&-1&0&\cdots&0&0&0\\ 0&0&0&1&-1&\cdots&0&0&0\\ 0&0&0&0&1&\cdots&0&0&0\\ \vdots&\vdots&\vdots&\vdots&\vdots&\ddots&\vdots&\vdots&\vdots\\ 0&0&0&0&0&\cdots&-1&-1&0\\ 0&0&0&0&0&\cdots&1&-1&0\\ 1&1&1&1&1&1&1&1&1\end{bmatrix} (5)

This construction matrix has the unique solution 1n,1n,2n,,2n,1n,1n\langle\frac{1}{n},\frac{1}{n},\frac{2}{n},\dots,\frac{2}{n},\frac{1}{n},\frac{1}{n}\rangle.

The pot of tiles where nn is even produces the following construction matrix:

M(Podd)=[(n1)1111110021000000011000000011000000010000000010000000120111111111]M(P_{odd})=\begin{bmatrix}(n-1)&-1&-1&-1&-1&\cdots&-1&-1&0\\ 0&2&-1&0&0&\cdots&0&0&0\\ 0&0&1&-1&0&\cdots&0&0&0\\ 0&0&0&1&-1&\cdots&0&0&0\\ 0&0&0&0&1&\cdots&0&0&0\\ \vdots&\vdots&\vdots&\vdots&\vdots&\ddots&\vdots&\vdots&\vdots\\ 0&0&0&0&0&\cdots&-1&0&0\\ 0&0&0&0&0&\cdots&1&-2&0\\ 1&1&1&1&1&1&1&1&1\end{bmatrix} (6)

This construction matrix has the unique solution 1n,1n,2n,,2n,1n\langle\frac{1}{n},\frac{1}{n},\frac{2}{n},\dots,\frac{2}{n},\frac{1}{n}\rangle. By Proposition 2, the smallest graph that can be realized by this pot is on nn vertices.

We must now show that the proposed pots will not realize graphs nonisomorphic to WnW_{n}. In both proposed pots, the solutions from (5) and (6) give the ratio of tile types to be used when constructing a graph of order nn. In both pots, only one copy of tile t1t_{1} is used, and since every other tile has one bond-edge type of type a1^\hat{a_{1}}, and t1t_{1} is only tile containing a bond-edge of type a1a_{1}, then t1t_{1} must bond with every other tile from PP in order to form a complete complex. In both even and odd cases, tile type t2t_{2} is used exactly once. The only tile type that is able to bond with the unattached edges of t2t_{2} is t3t_{3}, so the two copies of t3t_{3} must bond with t2t_{2}. Once the two copies of t3t_{3} bond to t2t_{2} and to t1t_{1}, then each of these two copies has a free cohesive end that can only be complemented by a cohesive end from t4t_{4}, so each copy of t4t_{4} must bond with one of the copies of t3t_{3}. In general, for 3in2+1,3\leq i\leq\lfloor\frac{n}{2}\rfloor+1, tile tit_{i} can only bond with tiles t1t_{1}, ti1,t_{i-1}, and ti+1t_{i+1}. For cases in which nn is even, tn2+2t_{\lfloor\frac{n}{2}\rfloor+2} must bond with tiles tn2t_{\lfloor\frac{n}{2}\rfloor} and tn2+1t_{\lfloor\frac{n}{2}\rfloor+1}. For cases in which nn is odd, tn2+2t_{\lfloor\frac{n}{2}\rfloor+2} must bond with two copies of tile tn2+1t_{\lfloor\frac{n}{2}\rfloor+1}. The resulting complex is isomorphic to WnW_{n}.

Thus, B3(Wn)=n2+1B_{3}(W_{n})=\left\lfloor\frac{n}{2}\right\rfloor+1. ∎

The corresponding result for minimum number of tile types needed to construct WnW_{n} follows from Theorem 3.

Theorem 4.

T3(Wn)=n2+2T_{3}(W_{n})=\lfloor\frac{n}{2}\rfloor+2.

Proof.

Since the only vertex of degree n1n-1 is the hub, this vertex must have its own tile type.

The pots in (3) and (4) suggest that the remaining vertices, all of which are of degree 3, can be labeled with n2+1\lfloor\frac{n}{2}\rfloor+1 tile types.

If it were possible to construct the outer cycle of WnW_{n} using exactly n2\lfloor\frac{n}{2}\rfloor tile types, then every tile on the outer cycle would appear twice if nn is odd, and if nn is even, exactly one tile type on the outer cycle would appear once. Suppose for sake of contradiction that n2\lfloor\frac{n}{2}\rfloor tile types may be used on the outer cycle. Label some pair of tiles on the outer cycle t1t_{1}, where t1={a,b,c}t_{1}=\{a,b,c\} (where a,b,ca,b,c are not necessarily distinct). By Lemma 1, these two copies must have a distance between them of at least two.

Refer to caption
Refer to caption
Figure 14: Two copies of t1t_{1} are placed on the outer cycle. Paths along the outer cycle must be of odd or even length.

By Lemma 2, a single bond type, say aa, must be used on both spoke edges. By the proof of Lemma 3, each of the remaining two bond-edge types of t1t_{1} must be oriented opposingly along the outer cycle (one bond-edge type, say bb, is colored in blue in Figure 14). Furthermore, since each bond-edge type has appeared twice, then by Lemma 3, the two vertices adjacent to those of t1t_{1} via the bond-edge type bb must be of the same tile type, say t2t_{2}. By the nature of the wheel graph, there must exist two paths between the two copies of t1t_{1} whose union is exactly the outer cycle.

If at least one such path has odd length, then we see on the right side of Figure 14 that applying Lemmas 2 and 3 to tile type t2t_{2} and each subsequent tile type, eventually some tile type will be adjacent to itself, which violates Lemma 1. This contradicts the assumption that the outer cycle of WnW_{n} could be constructed using n2\lfloor\frac{n}{2}\rfloor tile types.

If both paths along the outer cycle are of even length, then nn must be odd. We see on the left side of Figure 14 that repeating the reasoning above, some tile type will be used exactly once, which contradicts the assumption that every tile could be used twice.

Alternatively, if the outer cycle of WnW_{n} could be constructed in strictly less than n2\lfloor\frac{n}{2}\rfloor tile types, then at least one tile type would appear three or more times on a vertex on the outer cycle. Suppose without loss of generality that the tile type is t={a,b,c}t=\{a,b,c\} (where a,b,ca,b,c may not necessarily be distinct bond-edge types). By Lemma 1, tile type tt may not be adjacent to itself, so we know each instance of tt must be at least distance 2 from itself. To satisfy Lemma 2, some bond-edge type, say aa, must be used on all the spoke edges. This leaves three edges on the outer cycle to be labeled with the bond-edge types bb and cc, which violates Lemma 3. Thus, we may not use fewer than n2\lfloor\frac{n}{2}\rfloor tile types to construct WnW_{n}.

Since each case shows that a pot using fewer than n2+2\lfloor\frac{n}{2}\rfloor+2 tiles violates the conditions of Scenario 3, then T3(Wn)n2+2T_{3}(W_{n})\geq\lfloor\frac{n}{2}\rfloor+2. Since the pots given in Equations (3) and (4) achieve this lower bound, then T3(Wn)=n2+2T_{3}(W_{n})=\lfloor\frac{n}{2}\rfloor+2. ∎

4 Conclusion

In this paper we have presented the minimum number of distinct tile types and bond-edge types needed to construct a self-assembled DNA complex with a wheel graph structure. There are still several other infinite classes of graphs beyond wheel graphs to study, and we anticipate that the results presented here may aid future authors in finding similar results for graphs which contain WnW_{n} as a subgraph.

The minimum number of distinct tile and bond edge-types needed to construct the wheel graph in Scenarios 1 and 2 are not only equal, but also constant, i.e., they do not increase as a function of the number of vertices of the target graph. The minima in Scenario 3 do increase as a function of nn, which is to be expected, given the more stringent conditions that come with Scenario 3.

In Scenario 3, we established two lemmas that aided us in finding lower bounds for the minimum number of tile types and bond-edge types. We believe this method of finding restrictions on the construction of a given complex may be generalized or adapted to find more results for additional classes of graphs.

5 Acknowledgements

This work was supported in part by the CSU San Bernardino Office of Student Research.

References

References

  • [1] Leonard M Adleman. Molecular computation of solutions to combinatorial problems. Science, 266(5187):1021–1024, 1994.
  • [2] Joanna Ellis-Monaghan, Greta Pangborn, Laura Beaudin, David Miller, Nick Bruno, and Akie Hashimoto. Minimal tile and bond-edge types for self-assembling dna graphs. In Discrete and Topological Models in Molecular Biology, pages 241–270. Springer, 2014.
  • [3] Hongzhou Gu, Jie Chao, Shou-Jun Xiao, and Nadrian C Seeman. A proximity-based programmable dna nanoscale assembly line. Nature, 465(7295):202, 2010.
  • [4] Thom H LaBean and Hanying Li. Constructing novel materials with dna. Nano Today, 2(2):26–35, 2007.
  • [5] Andre V. Pinheiro, Dongran Han, William M. Shih, and Hao Yan. Challenges and opportunities for structural dna nanotechnology. Nature Nanotechnology, 6(12):763–773, 2011.
  • [6] Nadrian C Seeman. An overview of structural dna nanotechnology. Molecular biotechnology, 37(3):246, 2007.