Economical routes to size-specific assembly of self-closing structures

Abstract

Self-assembly is one of the prevalent strategies used by living systems to fabricate ensembles of precision nanometer-scale structures and devices. The push for analogous approaches to create synthetic nanomaterials has led to the development of a large class of programmable crystalline structures. However, many applications require ‘self-limiting’ assemblies, which autonomously terminate growth at a well-defined size and geometry. For example, curved architectures such as tubules, vesicles, or capsids can be designed to self-close at a particular size, symmetry, and topology. But developing synthetic strategies for self-closing assembly has been challenging, in part because such structures are prone to polymorphism that arises from thermal fluctuations of their local curvature, a problem that worsens with increased target size. Here we demonstrate a strategy to eliminate this source of polymorphism in self-closing assembly of tubules by increasing the assembly complexity. In the limit of single-component assembly, we find that thermal fluctuations allow the system to assemble nearby, off-target structures with varying widths, helicities, and chirality. By increasing the number of distinct components, we reduce the density of off-target states, thereby increasing the selectivity of a user-specified target structure to nearly 100%. We further show that by reducing the design constraints by targeting either the pitch or the width of tubules, fewer components are needed to reach complete selectivity. Combining experiments with theory, our results reveal an economical limit, which determines the minimum number of components that are required to create arbitrary assembly sizes with full selectivity. In the future, this approach could be extended to more complex self-limited structures, such as shells or triply periodic surfaces.

Introduction

The design and control of self-assembly pathways is a promising route for creating complex, functional nanomaterials [1]. Recent successes in colloidal self-assembly focus primarily on synthesizing spatially unbounded, dense crystalline materials, or spatially-limited architectures, like clusters, membranes, and filaments, whose dimensions are specified by the building-block sizes [2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15]. However, Nature is brimming with examples from a different class of structures that have self-regulated cavity sizes that are much larger than the size of the individual building blocks. Importantly, these self-limiting cavity sizes are essential to the functionality of various biological devices and materials, such as responsive containers that can selectively package specific genetic material, as in viral capsids [16], or photonic nanostructures, like those found in some butterfly wings that produce their structural coloration [17]. Inspired by these examples, synthetic colloidal assembly has recently taken a large leap forward by developing new colloidal building blocks with directional binding, complex geometries, and specific interactions that can target assembly of similar self-limiting architectures, such as icosahedral shells and cylindrical tubules [18, 19, 20].

A common strategy used by Nature to assemble self-limiting structures exploits self-closure [21], in which accumulated curvature between bound subunits allows the assembly to close upon itself during growth, thereby terminating assembly at a finite size. While translating this principle of self-closure could open new doors in synthetic colloidal assembly, self-closure has an associated fundamental challenge that must be solved to assemble precise, self-limiting architectures: Thermal fluctuations give rise to variations in the curvature of the growing assembly, leading to a distribution of the cavity sizes that result in the final structures [22]. This problem relates to the fact that the programmed angles between subunits lead to closed loops around the cavities, as seen for different discrete assemblies in Fig. 1A. Though a certain average curvature may be targeted, the finite bending rigidity of the assemblies leads to a spread of sizes that these loops can form [23, 21, 24, 25]. This aspect of assembly is easily seen in the context of forming a ring, Fig. 1A, in which the collective fluctuations of many directional bonds can admit larger or smaller rings [18]. The larger the number of subunits in the ring, the broader the distribution of states that can be accessed, and thus the ability to selectively assemble the target ring becomes more difficult with increasing ring size. Because the self-limiting cavity sizes are fundamental to how these types of structures function, developing generally applicable solutions to this challenge is essential to realizing synthetic self-limiting architectures with functionalities that rival their natural counterparts.

Figure 1: Improving selectivity in self-closing assemblies. (A) Examples of self-closing architectures—the unit cell of a P-surface (right), a cylinder (middle), and an icosahedral shell (left)—and closed loops of monomers (denoted in red). Fluctuations in the angles between subunits lead to a distribution of closure sizes in loops. (B) Schematic of closure fluctuations for a tubule. Tubules form from wrapping a sheet while satisfying a periodic boundary condition. Any discrete assembly can be represented by a unique vector between two vertices on the triangulated plane, corresponding to the circumference of that tubule. The pair of numbers ($m,n$) represents the numbers of steps in the $\mathbf{m}$ and $\mathbf{n}$ directions of the lattice to create this circumference vector. States accessible to thermal fluctuations around the target state lie within the dotted circle. Examples of target and off-target assemblies are shown. (C) Scheme for removing off-target states. Accessible states are shown as purple and states allowed by the coloring pattern are white. As the number of colors increases, the distance between similar vertices grows and eventually removes all allowed off-target states from the fluctuation area (red circle). (D) Experimental scheme for realizing colorings of tubules. Lines on the monomer edges denote ssDNA locations. Arrows between colored monomers show binding rules. (E) Example TEM micrograph for the 16-color tubule. The bottom shows a single layer of the tubule. 10-nm-diameter gold-nanoparticles identify the red component. Scale bar is 200 nm.

In this article, we revisit the initial promise of programmable assembly, and demonstrate a general strategy for eliminating polymorphism in self-closing assembly by increasing the assembly complexity. By combining theory and experiment we reveal an economical limit, in which the selectivity of a single target structure can be increased to 100% using a minimal number of components. Our approach extends symmetry-based theories, such as that from Caspar and Klug [16], which identify the minimum number of components required to form a given structure. In particular, we demonstrate our strategy for the specific case of assembling cylindrical tubules from triangular DNA origami subunits (Fig. 1B–E). We increase the complexity of our assemblies by incorporating a multi-component, periodic coloring pattern that specifies the arrangement of various components within the tubule. By using colorings with more types of components, fewer tubule geometries are commensurate with both the geometry of the tubule and the periodicity of the coloring, which reduces the number of off-target states that are thermally accessible. By borrowing ideas from two-dimensional addressable assembly [13, 26, 27], we decouple the interactions between subunits from the geometry of the monomer to realize assemblies built from up to 16 unique components, demonstrating that our strategy allows assembly with essentially arbitrarily high complexity. Ultimately, we identify an economical threshold—which grows as the size of the cavity squared—where a single assembly state can be selected with a minimal amount of coloring complexity and rationalize this limit using a simple Helfrich energy model. We conclude by showing that one can target lower dimensional properties, like the width or the pitch, which reduces the constraints on assembly and thus provides a more economical scaling that is proportional to the self-limiting size.

Improving specificity using multiple components

Our strategy works as follows. A cylindrical tubule can be conceptualized as a sheet that closes upon itself. Because you can tile a sheet with identical equilateral triangles, only one component is required [20] and any tubule state can be identified with a unique pair of numbers ($m,n$) that corresponds to the shortest closed path around the tubule (directions $m$ and $n$ in Fig. 1B). But having only a single component also allows the sheet to close on itself in many different ways since all vertices of the triangulation are identical, Fig. 1B. While a single tubule state can be preferred, as specified by the assembly’s curvature, the finite bending rigidity of the sheet could admit many neighboring tubules, with similar curvatures and therefore similar bending energies. These accessible geometries can be thought of as off-target states that occupy an area of vertices around some ground-state vertex. The challenge of eliminating polymorphism in the final structures then boils down to removing all undesired states from within this area.

To accomplish this goal, we color triangles in a periodic way using an increasing number of colors, which removes translational symmetries of the sheet and reduces the number of allowed closure states, Fig. 1C. As additional complexity of coloring is used, the allowed closure states are pushed further apart and eventually there are no allowed vertices within the area of fluctuations. At this point, only a single accessible assembly state is permitted by the matching rules.

Figure 2: Increased selectivity with increased complexity. (A,B,C) (i) Probability distributions of different tubule types for 1-, 4-, and 14-color assemblies with the (6,0) monomer. Each circle denotes a single ($m,n$) tubule type and the area of each point denotes the relative probability. Illustrations of the most probable tubule for each distribution are shown with each plot. (ii) Representative TEM images of (6,0) tubule assemblies with the 1-, 4-, and 14-color designs. (D,E,F) (i) Probability distributions for 1-, 4-, and 16-color assemblies with the (10,0) monomer. (ii) Representative TEM images of (10,0) tubule assemblies with the 1-, 4-, and 16-color designs. Tubules in B(ii),C(ii),E(ii), and F(ii) have the corresponding tiling overlayed with a single color labeled with 10-nm-diameter gold nanoparticles. Only positive $n$ points are shown since single TEM images cannot differentiate between left and right-handed structures. All scale bars are 200 nm.

We develop an experimental system to test this concept, comprised of triangular subunits made via DNA origami, that targets a user-specified cylindrical tubule by controlling the interaction specificity and the dihedral angles between neighboring subunits [20]. To encode the dihedral angles, we bevel each side of the triangle. The monomers also need specific interactions to preserve their orientation with respect to the tubule axis. To realize specific interactions, we place six single-stranded DNA (ssDNA) segments with six-base-long binding domains along each edge (Fig. 1D). Using sticky-end hybridization allows us to encode many low cross-talk interactions [28, 29, 30, 27, 15] without changing any internal routing of the DNA origami, which could otherwise have unintended effects on the monomer structure [31, 19, 20]. Finally, we assemble tubules at a constant temperature chosen such that the intersubunit interactions are weak and reversible.

We find that a single component assembles a distribution of tubules with varying width and helcity. We first design a triangular monomer that targets a (6,0) tubule and classify the tubules that assemble with transmission electron microscopy (TEM). In particular, we measure both the width and pitch of each tubule to identify its ($m,n$) type, and then construct a distribution of states, Fig. 2A. The distribution reveals that the monomer prefers to form (7,0) tubules—close to our target of (6,0)—but that over half of the observed tubules form other assembly states distributed near the preferred state, including both achiral and chiral tubules.

To circumvent the formation of off-target states, we engineer the free-energy landscape around the preferred state by building the tubule from a periodic arrangement of a larger number of components. More specifically, we generate a matrix of pairwise subunit interactions by finding periodic colorings of the plane for varying numbers of colors, Fig. 1C (Suppl. Sec. II). The set of adjacent colors specifies which unique, specific interactions are required to assemble that pattern, which we implement by designing the sticky-end sequences of our origami subunits. Each color then corresponds to a subunit with a unique set of intersubunit interactions. Going forward, we refer to a component as a ‘color’.

By choosing colorings that permit the preferred (7,0) state, we aim to reduce the density of nearby, undesired assemblies. Figures 2B and C show the distributions of 4- and 14-color assemblies. As the number of colors, $N_{\mathrm{colors}}$, increases we find that the density of available states decreases and the probability of the most likely tubule state, which we call ‘selectivity’, increases. Importantly, images of individual tubules show that the quality of assemblies does not diminish as $N_{\mathrm{colors}}$ increases, Fig. 2B(ii) and C(ii). We also confirm the specificity of the intersubunit interactions by labeling a single color with 10-nm-diameter gold nanoparticles and performing TEM tomography. While we cannot unambiguously infer the colors of all of the components, we can conclude that the pattern of gold nanoparticles is consistent with our designed coloring in all cases.

Figure 3: Scaling for increased selectivity. (A) Probability of the most common assembly state, called selectivity, against $N_{\mathrm{colors}}$. (B) Selectivity rescaled by the fluctuation area $N_{\mathrm{colors}}B/C^{2}$. $C$ values are taken from the most probable tubule state. Fitting the single-component assemblies to the Helfrich model provides estimated values of the bending rigidity of 10.0 k${}_{\mathrm{B}}T$  and 6.7 k${}_{\mathrm{B}}T$ for the (6,0) and (10,0) monomers, respectively. See Suppl. Sec. VI and Reference [24] for a description of the theory.

Because fluctuations of the dihedral angles between adjacent subunits give rise to the breadth of states, we expect that targeting a larger diameter tubule will result in a broader distribution of off-target structures and thus will require a greater degree of complexity to achieve full selectivity. To test this hypothesis we design a second DNA origami monomer with bevel angles targeting a (10,0) tubule. Figure 2D shows the distribution of states for this monomer and an example of an assembled tubule. As anticipated, the number of off-target states increases in comparison to the (6,0) tubule. Again, we can improve selectivity through increased assembly complexity. Figures 2E and F show 4- and 16-color assemblies and a corresponding increase in the assembly selectivity. As before, we see that the quality of tubules is preserved with increasing complexity and that the interaction specificity is consistent with our designed colorings, as seen in the TEM images.

We use a simple Helfrich model to relate the distribution of available states to the physical properties of the assembly [32, 24, 33]. The model assumes that tubule assembly begins by the formation of a circular patch, and thus the area of states that can form around the target due to thermal fluctuations is $C^{2}/B$, where $C$ is the preferred circumference non-dimensionalized by the monomer edge length and $B$ is the bending rigidity non-dimensionalized by k${}_{\mathrm{B}}T$ (Suppl. Sec. VI). This area is independent of the assembly complexity. Our coloring patterns introduce another lengthscale: the distance between similar vertices. This length grows with the size of the primitive cell of the coloring, whose area increases as $N_{\mathrm{colors}}$ (Fig. 1C). When these two areas are comparable, we expect that the nearest allowed off-target assembly is outside the area accessible by thermal fluctuations. Thus, we expect a crossover to 100% selectivity when

$$N_{\mathrm{colors}}B/C^{2}\gtrsim 1.$$ (1)

Our data of the selectivity collapse when replotted according to our Helfrich model. In particular, Fig. 3 shows that the selectivity increases linearly with an increasing number of colors and plateaus to full selectivity when $N_{\mathrm{colors}}B/C^{2}\approx 1$. This crossover point corresponds to the most economical way to target a single assembly state, after which additional complexity provides no further benefit.

Limiting lower-dimensional properties of assemblies

While the promise of complete selectively is tantalizing, the scaling of the number of components with the square of the self-limited length becomes daunting for larger structures. Therefore, we explore the possibility that targeting a lower-dimensional quantity, like the width or pitch, could be accomplished using fewer components. For each tubule geometry, i.e., each ($m,n$), there is a unique width and pitch of the assembly. But there also exist families of tubules that have nearly constant width or pitch. Figure 4A shows iso-contours of these properties in ($m,n$) space. We note that lines of constant $n$ or $m+n$ approximate lines of constant pitch or width, respectively. Therefore, we are able to target families of nearly constant width or pitch using linear colorings, in which seams of the same color lie along a single lattice direction, by aligning them along either $\mathbf{n}$ or $\mathbf{m-n}$ directions, Fig. 4B.

Figure 4: Limiting single degrees of freedom in assembly. (A) Iso-contours of the pitch and width for tubules for varying $m$ and $n$. (B) Schematic of allowed states for linear colorings for pitch or width alignment and (C) TEM images of 4-color width- and pitch-controlled tubules with the (10,0) monomer. Lines show seams of labeled particles (solid and dashed lines denote the top and bottom layers of the tubule). (D) Probability distributions for the pitch and width of multi-component tubules compared to one-color tubules. Units for the pitch and width are the edge length of a monomer. Complete ($m,n$) distributions of these assemblies are shown in (A). (E) Selectivity of tubules that have the same $m+n$ (width-control, circles) or the same $n$ (pitch-control, squares). The dashed lines are theoretical predictions from Fig. S5. (F) Selectivity of the unconstrained dimension versus $N_{\mathrm{colors}}$.

We validate this scheme in experiment and confirm that it yields tubule distributions that are highly selective in their width or pitch. By labeling a single color with gold nanoparticles we can see that lines of the same color correspond to lattice directions that are closest to the circumferential or axial directions for the pitch- and width-controlled tubules, respectively (Fig. 4C). We also compare the distributions of tubule widths and pitch for two of our (10,0) tubule experiments against a single-color experiment, and find that our linear colorings create assemblies with tightly peaked distributions for the desired property, Fig. 4D (Suppl. Fig. S9).

Using a similar argument from the Helfrich model, we find that placing reduced constraints on the assembly leads to an economical threshold that scales linearly with the self-limited length scale, rather than quadratically as before. Since we target seams of similar states, we are now concerned not with the area of thermal fluctuations, but with their linear extent, $C/\sqrt{B}$. For linear colorings, the separation between the same component type grows as $N_{\mathrm{colors}}$.

To compare our experimental data to the Helfrich model, we plot the selectivity of tubule states with the same $n$ or $m+n$ against $N_{\mathrm{colors}}\sqrt{B}/C$. Figure 4E shows that our data collapse in agreement with the Helfrich model, with the economical point occurring when $N_{\mathrm{colors}}\sqrt{B}/C\approx 1$. However, we find that reducing the design constraints comes at the cost that the unconstrained dimension (e.g. the pitch when one targets the width) does not see any improvement in selectivity for increasing complexity, Fig. 4F.

Whereas we have so far focused on constraining the ($m,n$) states that can form, we conclude by showing how assembly complexity can also be used to constrain a different aspect of the tubule geometry: their length. We use the same class of linear colorings as above, but now we passivate a single edge of one color to prevent the further addition of monomers along that seam, Fig. 5. In this way, the number of colors constrains the tubule length. In Fig. 5 we show representative TEM micrographs of our (6,0) monomer with 3-, 7-, and 19-color assemblies. In all cases, we find that the system reaches an equilibrium state characterized by a mixture of monomers, small clusters, and finite, $n$-length assemblies (Suppl. Sec. VII). The spirit of this type of constraint is different than those mentioned previously since we are not increasing the complexity to reduce the number of states accessible by thermal fluctuations, but rather to cap the growth of assemblies, highlighting the various ways that assembly complexity can be used to engineer self-limitation.

Figure 5: Tubule length limitation. For achiral tubules (such as those made by the (6,0) monomer), a linear coloring constrains tubules to length $n$ by neutering one edge of the $n$-th color. Dashed and solid red lines denote the neutered interactions for $n=$3 and 7. TEM images for 3-, 7-, and 19-color assemblies are shown. All scale bars are 100 nm.

Conclusions

In this work, we have demonstrated a scheme for preventing off-target assemblies by increasing the complexity of the initial assembly mixture. By imposing complex colorings, implemented through matrices of specific intersubunit attractions, assembly states that were accessible due to fluctuations in the curvature become disallowed, leading to full selectivity of the target geometry. By imposing different numbers of constraints on the assembly, we were able to either achieve full selectivity for both the pitch and width of tubules (Fig. 2 and  3) or either one separately (Fig. 4). The benefit of selecting a single dimension of the structure, say the width, is that it can be done using significantly fewer colors, which becomes essential as the self-limiting dimension becomes arbitrarily large relative to the subunit size. Unfortunately, targeting the pitch or width alone comes with a lack of selectivity for the unconstrained dimension, e.g. the pitch. We show examples of this trade-off in Suppl. Fig. S8. While in the limit of large target sizes selectivity eventually falls off, we showed that increasing the complexity mitigates the marginal losses of selectivity, enhanced by a factor of $N_{\mathrm{colors}}$ over the standard Helfrich expectation.

We anticipate that our strategy is generalizable to other self-closing architectures. A useful point of comparison is the Caspar-Klug (CK) construction for icosahedral capsids [16]. By utilizing the symmetry of icosahedral shells, the CK rules provide a method for determining the minimal number of components needed to form a shell of arbitrary size. However, for too low of a bending modulus, capsids can assemble asymmetric, defective structures. These defects often come from disclinations that form at 5-fold or 3-fold points. By increasing the number of colors around certain symmetry points in the structure, beyond the minimum prescribed by CK theory, we anticipate that strategies like the one we present here could reduce the propensity for these types of systems to form defects.

Despite the great improvement in specificity that we see using multi-color assemblies, there are still challenges that need to be overcome. The foremost is the kinetics of growth. The timescale for assembly increases with increasing complexity since the chance of the correct color monomer finding a binding site at the growth front decreases with $N_{\mathrm{colors}}$ (Suppl. Sec. VIII). This decrease in assembly rate could be compensated for by increasing the monomer concentration, but this strategy could become costly. Instead, future strategies may exploit hierarchical pathways to assembly both for accelerating assembly and possibly reducing the desired complexity if certain kinetic paths can be found that disallow off-target structures. An important aspect of the kinetics that this work reveals is that assembly near the economical point is favorable. In the 14-color assembly with the (6,0) monomer, even though all closed structures formed the same ($m,n$) state, a small fraction of assemblies missed the point of closure and overgrew as a sheet that wrapped around itself like a scroll. We hypothesize that this type of misassembly will become even more prevalent if the assembly complexity passes the point of optimal economy, $N_{\mathrm{colors}}B/C^{2}\gg 1$.

An important consideration with all types of multi-component assemblies is the possibility of crosstalk between different components that can stabilize unintended contacts [26, 34]. In this context, we can ask why our designs worked so well and do not seem to suffer from the inevitability of cross-talk between a limited set of DNA sequences [29]? Here we highlight that the site addressability of DNA origami allows for the inclusion of interaction geometry as an additional knob for tuning the specificity of interactions [27]. In fact, for some of our interaction patterns, we reused strands from other interactions but altered their arrangement on the faces of the interacting edges. This combination of sequence specificity and spatial commensurability greatly increases the number of unique interactions that can be specified by a finite library of DNA sequences. Furthermore, increasing the number of sticky ends per side could be used to increase this combinatoric aspect further, allowing for many more components than have been explored here.

Coloring, the process of finding allowed, distinguishable identities for particles in a structure, is an attractive way to increase the complexity of assemblies, though the inverse design challenge of finding an appropriate coloring for a desired outcome is not always straightforward. Here we were able to exploit the symmetries and periodicity of the tubule to generate any desired coloring. More complex geometries, such as capsids or gyroids, require considering local point symmetries of the structure [16, 19, 35, 36] to encode the correct interactions, while fully addressable clusters often reveal their allowed colorings through detailed searches of all possible interacitons [37]. Recent efforts have also shown that colorings can be used to encode hierarchical structures into crystalline self-assemblies [15]. As assembly structures become more complex, finding the right scheme to encode large libraries of interactions may require new methods, such as SAT-assembly [38, 39]. Additionally, it will be equally important to be able to predict and thus program the right binding free energies to facilitate robust assembly [40, 41, 42].

Going forward, our scheme could be expanded to systems with subunits with both unique interactions and unique geometries to target a wider range of self-limiting architectures. Whereas we focused on a base monomer with fixed geometry, other target structures have the need for varying subunit geometries, as with icosahedral shells [19] or for surfaces with varying Gaussian curvature, such as helical cylinders, toroids, or open crystalline structures [36]. Moreover, many studies have recently demonstrated that DNA origami can be dynamically reconfigured [43, 44, 45, 46], so it is within sight to imagine a base monomer with adjustable edge lengths or bevel angles that could be used to produce new nanoscale devices.