Abstract
Developing construction methods of materials tailored for given applications with absolute control over building block placement poses an immense challenge. DNA-coated colloids offer the possibility of realising programmable self-assembly, which, in principle, can assemble almost any structure in equilibrium, but remains challenging experimentally. Here, we propose an innovative system of linker-mediated mobile DNA-coated colloids (mDNACCs), in which mDNACCs are bridged by the free DNA linkers in solution, whose two single-stranded DNA tails can bind with specific single-stranded DNA receptors of complementary sequence coated on colloids. We formulate a mean-field theory efficiently calculating the effective interaction between mDNACCs, where the entropy of DNA linkers plays a nontrivial role. Particularly, when the binding between free DNA linkers in solution and the corresponding receptors on mDNACCs is strong, the linker-mediated colloidal interaction is determined by the linker entropy depending on the linker concentration.
INTRODUCTION
The ultimate goal of self-assembly is to program many distinct building blocks, and each of them occupies a specific location within a self-assembled structure (1). The recent development of DNA nanotechnology offers possibilities of programmable self-assembly using the specific hybridization between single-stranded DNAs (ssDNAs) (2–5). This works very well in programmable self-assembly of DNA bricks, in which a variety of designed superstructures consisting of thousands of preprogrammed DNA bricks were fabricated (6–9). However, similar ideas were not well applied in the designed self-assembly of DNA-coated colloids (DNACCs). One of the major challenges is that typically each colloid is coated with many DNA linkers, and the effective colloidal interaction mediated by DNA hybridization changes abruptly with temperature, which makes the system difficult reach the equilibrium ordered state (1). Thus, only a few groups were able to obtain three-dimensional (3D) crystals of DNACCs (10–17).
This is particularly detrimental for designed self-assembly of colloidal superstructures, where the temperature window for high-yield self-assembly narrows down quickly with the increasing structure size (18). In this work, we propose a new system of linker-mediated mobile DNACCs (linker-mediated mDNACCs). Experimentally, one can fabricate mDNACCs by grafting DNA linkers onto the lipid bilayer coated on the colloids (19, 20), which makes the grafted linkers mobile on colloidal surface. It was found that the phase diagram of conventional mDNACCs is not qualitatively different from the corresponding immobile DNACCs, where the freezing colloidal density drops to zero quickly with increasing the binding strength between ssDNAs (21). Here, we find that if the ssDNAs grafted on different colloids do not bind with each other directly, but rather through the bridging of free DNA linkers in solution, then the effective colloidal attraction does not diverge at the strong binding limit, i.e., low-temperature limit. The reason is due to the special entropic effect in the strong binding limit making the linker-mediated attraction between mDNACCs finite and solely depending on the concentration of free DNA linkers in solution, which can be well controlled over orders of magnitude experimentally.
RESULTS
Model and mean-field theory
We consider an equimolar AB-type binary system of volume V consisting of N mDNACCs (hard spheres) of radius R in the solution containing free DNA linkers of chemical potential μ. Each mDNACC i is coated with ni A- or B-type mobile ssDNA receptors of length rc, which can bind to the ssDNA tail of complementary sequence on a free DNA linker in solution with the binding free energy ΔGbind (Fig. 1A). Free DNA linkers are modeled as infinitely thin hard rods of length l, with two ssDNA ends of length rc (Fig. 1B). Here, we assume R ≫ l ≫ rc, and the free DNA linkers in solution and the mobile ssDNA receptors move much faster than colloids. During the motion of colloids, all linkers and receptors reach equilibrium quickly. One can write down the partition function and calculate the free energy of the linker system using the saddle point approximation. This can be used as the effective interaction between mDNACCs mediated by linkers (22–24), which has the contributions from both bonded linkers on mDNACCs and unbound free linkers in solution. If we assume that the excluded volume effect between the grafting sites of ssDNA receptors on the same particle can be approximated using the hard disk–like repulsion on a 2D surface, when the area fraction of the grafted ssDNA receptors on the particle is less than 5%, then the compressibility factor is less than 1.1 (25), and the interaction between the grafting sites can be approximated as an ideal gas. Moreover, we assume that the concentration of DNA linkers in solution is low, and except the binding between complementary ssDNAs, the interaction between them is negligible. These assumptions are typical in experiments. Thus, the grand canonical partition function for the bonded DNA linkers on mDNACCs is
(A) Schematic representation of binary mNDACCs, in which the colloidal interaction is mediated via the bridging of free DNA linkers that bind with the mobile ssDNA receptors on mDNACCs. (B) A magnification of the yellow region in (A). Here, mobile ssDNA receptors (red or blue) coated on mDNACCs can bind with the ssDNA tails (red or blue spheres) of free linkers of the same color. (C) The linker-mediated effective interaction βUeff(r) as a function of the center-to-center distance between two mDNACCs r for different βμ with βΔGbind = −3.0, ni = 50, l = 0.2R, and rc = 0.01R. Solid lines are from Eq. 7, and symbols are from the direct MC simulations with explicit linkers.
Using the free energy of the bonded linkers ℱ({mi, qij}), we can rewrite Eq. 1 into Z = ∑{mi, qij} exp (−βℱ({mi, qij})), and with the saddle point approximation ∂ℱ({mi, qij})/∂{mi, qij} = 0, we obtain
Combining Eqs. 3 and 4, we obtain a set of self-consistent equations
Then, the free energy of bonded DNA linkers at the saddle point is
Besides the free energy contribution from the linkers bonded on mDNACCs, the unbound free linkers in solution also contributes a depletion effect Udep (26, 27) (see section S2), and the resulting linker-mediated effective interaction between mDNACCs can be written as
Numerical verification of the mean-field theory
We perform MC simulations with explicit DNA linkers to verify Eq. 7. As the system spans multiple length scales, direct simulations of many mDNACCs with explicit DNA linkers are prohibitively expensive, and we choose to perform thermodynamic integrations (see Methods) to calculate the linker-mediated effective interaction between two different mDNACCs (like in Fig. 2C). We plot the calculated linker-mediated effective interaction βUeff(r) between two colloids coated with the same number but different types of mobile ssDNA receptors compared with the theoretical prediction of Eq. 7 in Fig. 1C for various μ. One can see that the numerically calculated effective interaction between two mDNACCs agrees well with the theoretical prediction, and Eq. 7 generally predicts a slightly stronger attraction than the one measured in direct simulations. The reason causing the discrepancy is twofold. First, it is known that the mean-field approach developed here generally predicts a slightly stronger attraction between mDNACCs compared with direct simulations with explicit linkers, as the mean-field theory assumes that the probability for any ssDNA receptor coated on colloids to be unbound is uncorrelated to the probability of an ssDNA tail on a free linker to be unbound (24, 30). The difference is visible when the number of ssDNA receptors or the free DNA linkers in the system is very small, which is not the case in our simulations. Second, the analytical form of Fcnf is only exact at rc → 0, while for finite small rc, it slightly overestimates the conformational entropy for bridging linkers. We believe that the overestimation of Fcnf is the major (but not the only) reason responsible for the overestimation on the amount of bridging linkers and stronger attraction (see section S2).
(A and B) βUeff(2R + 2rc) (A) and
As shown in Fig. 1C, βUeff(r) is negative and attractive at 2R < r < 2R + l with the minimum located around 2R + 2rc, whose magnitude indicates the strength of attraction between two mDNACCs. As shown in Fig. 2A, with increasing βμ from −10 to about 0, the attraction between two mDNACCs first becomes stronger, i.e., βUeff(2R + 2rc) becomes more negative, while further increasing βμ makes the attraction weaker. Simultaneously, the probability of forming bridges between two colloids
When μ → ∞, all ssDNA receptors on mDNACCs are bonded, i.e.,
(A) Singlet fraction in binary linker-mediated mDNACCs (packing fraction η = 0.1) as a function of βμ for different βΔGbind (with lines as a guide to the eye). (B to D) The typical snapshots of the system at βμ = −6.5 (B), −2.5 (C), and 2.5 (D), where nonsinglets are drawn 10 times smaller. Here, ni = 50 and βΔGbind = −4.0 at the packing fraction η = 0.1.
Entropic effects in the strong binding limit
Furthermore, we investigate the effective interaction between mDNACCs with changing the binding strength between DNA linkers and corresponding receptors, i.e., decreasing ΔGbind. One might think that the stronger binding between DNA linkers and receptors would naturally make the attraction between mDNACCs stronger and eventually diverge at ΔGbind → −∞. Intriguingly, however, as shown in Fig. 4A, at ΔGbind → −∞, βUeff(2R + 2rc) does not diverge but reaches a plateau depending on μ. The MC simulations for systems of many mDNACCs also show that the singlet fraction drops with decreasing ΔGbind and reach a plateau at the strong binding limit (Fig. 4B). To understand this counterintuitive phenomenon, we modify our mean-field theory as follows (see section S4). When ΔGbind → − ∞, all receptors on mDNACCs are occupied, i.e.,
(A) Effective pair interaction βUeff(2R + 2rc) as a function of βΔGbind for various βμ predicted by Eqs. 7 (solid lines) and 11 (dashed lines). (B) Singlet fraction in binary linker-mediated mDNACCs at η = 0.1 as a function of βΔGbind for different βμ (with lines as a guide to the eye). Here, ni = 50. (C) Effective pair interaction
Here, the “enthalpy” term ∑iniβΔGbind does not depend on the colloidal configuration and can be neglected, and ξaeβ(ΔGbind + μ) = Va′eβ(μ − Frep) is the entropy part of the partition function and does not depend on ΔGbind. Therefore, the resulting effective colloidal interaction at the strong binding limit is
It is known that for systems of short-range attractive particles, at the infinitely strong attraction limit, vibration entropy stabilizes the floppy crystals (21), which first nucleate from the dilute fluid. Using Eqs. 7 and 11, we perform isothermal-isobaric (NPT) MC simulations for the binary system of linker-mediated mDNACCs at βμ = 1.6 at the strong binding, i.e., βΔGbind = −6, and the calculated equation of state (EOS) is shown in Fig. 5. Compared with the EOS at βΔGbind → −∞, one can see that the random fluid first nucleates into CsCl at the pressure PR3/kBT ≃ 1.1, and this is qualitatively different from the situation of conventional mDNACC systems, in which the crystallization pressure of the system approaches zero at the strong binding limit (21).
EOS of the linker-mediated mDNACC systems at βμ = 1.6 with different βΔGbind, where the filled and open symbols are the results from NPT simulations expanding from perfect CsCl crystals and compressing from random fluids, respectively. The insets are the snapshots of fluid and nucleated CsCl crystal, respectively.
Generalization to multicomponent systems
In this section, we generalize our mean-field theory to calculate the effective interaction in a multicomponent system consisting of Nc types of mDNACCs, in which the total number of mDNACCs is N. mDNACC i is coated with only one type of ssDNA receptors, e.g., type I, and we regard the mDNACCs coated with type I ssDNA receptors as the type I particles. The number of ssDNA receptors coated on mDNACC i is ni. There are multiple types of free linkers in solution, which are described by a Nc × Nc connectivity matrix M, and if the free linkers connecting mDNACCs of type I and J exist in the solution, then MIJ = 1; otherwise, MIJ = 0. We assume that the lengths of all free linkers in solution and all ssDNA receptors coated on mDNACCs are l and rc, respectively. Then, the grand canonical partition function of the bonded linker is
With these, we further obtain a set of self-consistent equations
Similarly, when βΔGbind, I → −∞ ,1 ≤ I ≤ Nc, one can obtain the free energy of the bonded DNA linkers as
DISCUSSION
In conclusion, we have proposed a linker-mediated mDNACC system, in which the interaction between mDNACCs is bridged by the free DNA linkers in solution. We formulate a mean-field theory to calculate the effective interaction between mDNACCs, which well agrees with numerical simulations with explicit DNA linkers. The mean-field theory can be further used to construct MC simulations for efficiently simulating the collective self-assembly of linker-mediated mDNACCs. Furthermore, combining analytic theories with numerical simulations, we find previously unknown entropic effects in linker-mediated mDNACC systems. With increasing the concentration of free DNA linkers from zero, the linker-mediated effective interaction between mDNACCs changes nonmonotonically, and the strongest colloidal interaction appears at some intermediate free linker concentration, which induces a re-entrant melting in linker-mediated mDNACCs. Moreover, at fixed free linker concentration, with increasing the binding strength between free linkers in solution and the receptors on mDNACCs, we find that the linker-mediated attraction between mDNACCs becomes stronger and reaches a plateau at the strong binding limit, i.e., ΔGbind → −∞. This is due to the fact that at the strong binding limit, all receptors on mDNACCs are bonded. Therefore, whether forming bridges between mDNACCs does not change the enthalpy of the system, and the linker-mediated interaction between mDNACCs is dominated by the entropy of DNA linkers, which depends on the concentration of free linkers in solution. Here, we model the DNA linkers as stiff rods, and to model realistic DNA linkers that are essentially semiflexible polymers, one just needs to numerically calculate the partition functions ξa and ξb, while the theoretical framework and physics remain the same. Moreover, we anticipate that the mobility of the grafted ssDNA receptors does not change the qualitative physics. The entropy-driven effective colloidal interaction at the strong binding limit should also exist in the recent experimental system of linker-mediated immobile DNACCs, in which the re-entrant melting transition with increasing the concentration of free DNA linkers was already observed (31). As the concentration of free DNA linkers can be well controlled in experiments, this suggests a new way to precisely tune the colloidal interaction for addressable assembly of DNACCs to avoid the abrupt change of colloidal interaction with temperature in conventional systems of DNACCs. Similar types of entropy-driven colloidal interaction have also been observed in other systems of DNACCs, in which the effective interaction between DNACCs is induced by the direct hybridization of ssDNAs coated on colloids, and the general feature of these nondiverging colloidal interactions at the low temperature limit is that when the binding strength between ssDNAs is infinitely strong, the colloidal interaction is a result of the competition between different types of entropies, while the enthalpy does not change with the colloidal configuration (33, 34). Compared with these existing entropy-driven methods for addressable assembly of DNACCs (33, 34), the advantage of using linker-mediated mDNACCs is that to encode all possible specific interactions between N, different mDNACCs would need only N-distinct grafted sequences (31) instead of N(N − 1)/2 in the other systems, which is particularly important for large N (35), and each specific interaction can be easily switched on or off in situ by introducing or (using strong absorbers) removing the corresponding free DNA linkers, which offers an extra degree of freedom in controlling the pathway of designed assembly. Moreover, the mobile feature of the receptors coated on colloids makes it possible to derive the closed form for effective interactions between colloids (like Eqs. 7 and 11) to efficiently simulate, investigate, and design the collective self-assembly, for which direct simulations with explicit linkers are prohibitively expensive.
METHODS
Numerical calculation of the effective interaction between two mDNACCs with explicit DNA linkers
In the simulation box, we place two mDNACCs (hard spheres of radius R), and each of them is coated with ni different types (A or B) of ssDNA receptors of length rc, which can specifically bind with one end of the free DNA linkers of complementary sequences with length l. We consider that there is an absorption layer with thickness rc on particles, and the linker ends of complementary sequence in the layer can be bonded with a certain probability. Note that the linker ends can move in the absorption layer freely without any energy cost. The infinitely thin stiff linkers cannot overlap with the particles, and periodic boundary conditions are applied in all directions. As typically there are many linkers bonded onto each mDNACC, most of the trial moves of mDNACCs would be rejected. Therefore, in each of our MC simulations, we fix the position of two mDNACCs in the system with the center-to-center distance r, and the free energy of the system can be calculated through the thermodynamic integration
As our theory is applicable at R ≫ l ≫ rc, we simulate the system with l = 0.2R and rc = 0.01R.
To evaluate Eq. 19, we perform grand canonical (μVT) MC simulations, in which we fix the chemical potential of linkers and position of two particles in the system, and the number of DNA linkers changes under the control of chemical potential. Besides the conventional grand canonical MC moves used to add/remove linkers in the system, we also perform MC trial moves to randomly translate and rotate the free linkers. For bonded linkers with one free end, we perform two kinds of trial moves to speed up the equilibration: randomly rotate around the bonded end or rotate around the center of the bonded particle. For the binding/unbinding moves, we use the approach described in (36, 37). For these MC trial moves, we first randomly select a linker end in the adsorption layer, and if the end is in the unbound state, then a binding trial move is proposed; if the end is in the bound state, then an unbinding trial move is proposed. The formation of a new linker end-ssDNA bond is accepted with a probability
SUPPLEMENTARY MATERIALS
Supplementary material for this article is available at http://advances.sciencemag.org/cgi/content/full/6/21/eaaz6921/DC1
This is an open-access article distributed under the terms of the Creative Commons Attribution-NonCommercial license, which permits use, distribution, and reproduction in any medium, so long as the resultant use is not for commercial advantage and provided the original work is properly cited.
REFERENCES AND NOTES
- Copyright © 2020 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works. Distributed under a Creative Commons Attribution NonCommercial License 4.0 (CC BY-NC).