Specificity, flexibility and valence of DNA bonds for guided emulsion architecture

Abstract

A method of forming an end product by self-assembly of a first component having a patch of a linker component, such as DNA strands, cadherins, adhesive proteins and nanoparticle linkers. Such emulsions can be used to process personal care products, skin cremes, foods and animal feedstocks.

Description

SEQUENCE LISTING

The instant application contains a Sequence Listing which has been submitted electronically in ASCII format and is hereby incorporated by reference in its entirety. Said ASCII copy, created on May 30, 2014, is named 046434-0453_SL.txt and is 1,779 bytes in size.

FIELD OF THE INVENTION

This invention is directed to the self-assembly and directed controlled assembly of complex particle architectures. More particularly this invention is directed to methods, system and product materials employing DNA component sequences, cadherins, other adhesive proteins on cell membranes and solid nanoparticles serving as droplet linkers to construct a wide range of particle assemblages and architectures.

BACKGROUND OF THE INVENTION

Self-assembly of particles is of great interest for the design of complex particulate architectures to create smart nano-materials with tunable optical, mechanical or electronic properties. Grafting linking components, onto liquid interfaces of emulsions can lead to new architectural possibilities. In one example of particle linking components, the specific and programmable interaction between complementary DNA components has been used to assemble colloidal molecules with specific symmetries imposed by the positions of the grafted DNA.

Consequently it is desirable to further expand the application of selected linking components to improve formulation of particle architectures and create new systems for various applications.

SUMMARY OF THE INVENTION

In one preferred embodiment DNA interactions have been developed to establish that the size and number of these adhesion patches (valency) can be controlled. Valence values of 2 lead to flexible polymers of emulsion droplets, while valence values above 4 lead to rigid droplet networks. In one example, a simple thermodynamic model quantitatively describes the increase in the patch size with droplet radii, DNA concentration and the stiffness of the tether to the sticky-end. The patches are formed between droplets with complementary DNA strands or alternatively with complementary colloidal nanoparticles to mediate DNA binding between droplets. This emulsion system opens the route to directed self-assembly of more complex structures through distinct DNA bonds with varying strengths and controlled valence and flexibility.

DNA strands can be grafted onto thermal oil-in-water emulsions. Mixing two emulsions with complementary DNA strands leads to their specific binding through strong yet reversible adhesion patches. Unlike solid colloidal particles, liquid droplets are able to rearrange within the packed structure once they are bound together. Moreover, the deformation of the emulsion droplets, i.e. the size of the adhesion patch, provides a direct probe of the free energy of binding. A thermodynamic model has been developed to relate the adhesion size to the binding free energy and discover that the entropy loss upon binding plays an important role. The validity of the model has been tested by varying the DNA surface density the stiffness of the tether and the droplet size. This system sheds light on the mechanisms of adhesion between contacting liquid surfaces. Emulsion self-assembly leads to segregated floppy networks, which are amorphous materials with advantageous rheological properties. In addition, it has been determined that colloidal nanoparticles can serve as mediators of the DNA interaction between droplets. Controlling their concentration determines the valence of the droplets and enables us to uniquely create linear emulsion strings or those that fold into compact clusters. Consequently, this system and method provides a highly advantageous tool for self-assembly and controlled or directed applications.

In other embodiments cadherins, other adhesive proteins on cell membranes and selected nanoparticles can also be used to establish controlled self-assembly linkages. A variety of commercial applications are possible for these embodiments, such as personal care products, skin creams, foods and feedstocks for animals.

The foregoing summary is illustrative only and is not intended to be in any way limiting. In addition to the illustrative aspects, embodiments, and features described above, further aspects, embodiments, and features will become apparent by reference to the following drawings and the detailed description. These and other objects, advantages and features of the invention and related claims herein, together with the organization and manner of operation therefore, will become apparent from the following detailed description where taken in conjunction with the accompanying drawings wherein like elements have like numerals throughout the several drawings described below.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1A illustrates an adhesive DNA-emulsion system wherein silicone oil droplets are stabilized with phospholipids, some of which are labeled with fluorescently dyed streptavidin that allow the grafting of S or S′ DNA strands; S functionalized emulsion is dyed in green (“g”) while the S′ one is dyed in red (“r”); complementary S and S′ sticky ends (SEQ ID NOS 4 and 5, respectively) then bind and form adhesion patches enriched in DNA tethers between the droplets; in FIG. 1B colloids are coated with the S DNA strand and dyed in green (“g”) which can stick to S′ functionalized droplets and eventually bridge to droplets together through a colloidal patch; FIGS. 1C(1) and 1C(2) are for a temperature below the DNA melting temperature Tm=50° C. and the complementary thermal droplets aggregate into clusters (FIG. 1C(1)) that are disrupted above the Tm (FIG. 1C(2)); and FIGS. 1D(1) and 1D(2) similarly show the S coated particles only stick to the complementary S′ emulsion below Tm (see FIG. 1D(1)) to form large composite structures of multiple droplets linked through particles and those clusters as well as individual droplet-particle interactions dissociate above Tm (see FIG. 1D(2));

FIGS. 2A(1)-A(3) illustrate entropy-dependent adhesion patches between complementary emulsions with confocal imaging of S (green or “g”) and S′ (red or “r”) complementary emulsions interaction reveal the formation of adhesion patches for double stranded (dsb); FIGS. 2B(1)-B(3) are a single stranded (ssb) construct; and FIG. 2C shows diameter of the patch as a function of the squared radius <R> for two interacting droplets of radii R₁ and R₂ for dsb (squares) and ssb DNA (circles), with

$\langle R\rangle =\sqrt{\frac{2R_1^2{R}_2^2}{R_1^2+R_2^2}};$ data are fitted by the model (dashed lines) on FIG. 2D; growth of the patch size (diameter/<R>) is plotted as a function of time for dsb (squares) and ssb (circles) DNA; and triangles record the growth of ssb patches after addition of complementary strands, which turn into dsb patches after ˜40 min;

FIGS. 3A(1)-3A(3) illustrate dependence of patch size and intensity on DNA coverage, showing confocal images of the adhesion patches between complementary droplets for increasing DNA coverages on the droplets surface (DNA coverage from left to right: FIG. 3A(1) shows ˜17% (middle); FIG. 3A(2) shows ˜34% (right); and FIG. 3A(3) shows 100% which is ˜1400 strands/μm²; FIG. 3B shows patch angle θ and FIG. 3C shows relative intensity plotted as a function of the relative DNA coverage on the droplets and the experimental data are well fitted by the model (solid lines);

FIGS. 4A(1)-4C(4) show diffusivity of complementary bonds, wherein FIGS. 4A(1)-4A(5) show green emulsion droplet diffusing on a complementary droplet along with the adhesion patch; FIGS. 4B(1)-4B(3) show two S coated particles bound to a S′ coated droplet diffuse on its surface; FIG. 4B(4) shows the mean square displacement of the particles with time reveals a diffusive behavior with a diffusion coefficient D=0.012 μm²/s on the surface of phospholipid stabilized emulsions; and FIGS. 4C(1)-4C(5) show an emulsion network connected by specific DNA bonds;

FIGS. 5A(1)-5D(4) illustrate self-assembly of the colloid-emulsion system, wherein FIGS. 5A(1)-A(4) in particular show micrographs of various configurations, FIG. 5A(5) shows a schematic chain, and FIGS. 5A(6)-5A(10) show layer formations of polymer chains; polymer chains of emulsion droplets are stabilized by two adhesion patches formed using complementary nanoparticles; and the divalent emulsion droplets result from a low nanoparticle/droplet ratio ˜5; FIGS. 5B(1)-5B(5) show that at higher nanoparticle/droplet ratios more than two adhesion patches can form over time, leading in this case to a trivalent 20 structure; FIG. 5C(1) shows a schematic and 5C(2)-5C(5) show large particle/droplet ratios, ˜100, can produce multivalent droplets and lead to compact rigid structures in which the beads (green or “g”) all assemble between the droplets (red or “r”) contacts to minimize the system's energy, and FIG. 5D(1) shows a schematic polymer and FIGS. 5D(2)-5D(4) show time lapse formation of the colloidal nanoparticles assembling as rings between droplet contacts to maximize the amount of droplet/particle adhesive area with little droplet deformation;

FIG. 6 illustrates a block schematic diagram of a computer system to carry out determination of various products and associated models for creating selected chemical architecture of the invention;

FIG. 7 illustrates confocal imaging of the adhesion patches between the droplets. The fluorescence intensity along the droplets surface is analyzed with the Oval Profile plugin in image 7. Once the patch is identified from the circular intensity profile (see next figure), the fluorescence intensity along two radii inside the outside of the patch are also measured for contrast measurements;

FIG. 8 illustrates oval profile along the circumference of a droplet. The central plateau has an average intensity of (can this be moved down to be level with the line?) I_patch+I_noise=88±9 AU. We define the width of the peak as the diameter of the patch, which is here of 20 pixels or d_patch=(please move down so level on line) 1.2 μm. The average intensity in the peripheral region is the background noise estimated from the droplets central intensity is evaluated to be I_noise=9.8±3 AU. Averaging this among different slices significantly reduces the error bar. Assuming that the concentration of streptavidin is linearly proportional to the measured fluorescence intensity, the streptavidin enrichment in the patch is

$\frac{C_{\mathrm{patch}}}{C_{\mathrm{emulsion}}}=\frac{I_{\mathrm{patch}}}{I_{\mathrm{emulsion}}}=7.5.$ This value as well as d_patch are used in the main text in FIGS. 8 and 9, using statistics from different slices, patches and droplets;

FIG. 9a illustrates the bonds between complementary red and green emulsions are mobile, even after the structure reaches the maximum droplet connectivity (left). The distribution of the number of red-green sticky contacts can be measured with many confocal scan (see FIG. 9b); and

FIG. 10 illustrates development of emulsion polymeric chains beginning with a dilute suspension of S′ coated emulsion droplets is mixed with the complementary S coated particles top a ratio of particle/droplet ˜5%. The capillary is first tilted horizontally to force the creaming of the droplets in the top corner. The droplets are then brought together by an additional slight vertical tilting that makes them slide together and adhere through the particles to form emulsion polymeric chains.

DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS

Various methods and chemical components were implemented and are described herein to illustrate details of preferred embodiments. A number of “valence” values were used to control size and architecture development for product materials. These advantageous features and attributes will be illustrated by the following description and examples.

A model emulsion system using DNA as an example linker system is shown in FIG. 1A and is stabilized with a mixture of two phospholipids, Egg-Phosphatidylcholine (EPC) and PEG-biotinylated lipids, and a small amount of a conventional SDS surfactant. The biotinylated lipids are saturated with fluorescent streptavidin, which in turn binds to a second biotin from the functional DNA strand. The green emulsion (Alexa Fluor 488® streptavidin) is coated with the S sequence, while the red emulsion (Alexa Fluor 633® streptavidin) is coated with the S′ sequence (see Example section). In addition to the complementary sticky ends, the binders have an identical backbone of either single stranded or double stranded DNA that serves as a tether to the binder. Alternatively, we study a hybrid system of interacting colloidal nanoparticles and emulsions, as illustrated in FIG. 1B. In this system the emulsion is only decorated with the S′ DNA sequence (red or “r”), which interacts with colloids coated with the complementary S sequence (green or “g”). As a result, a colloidal particle binds droplets together.

The thermal reversibility of DNA interactions in these systems allows switching the adhesion between droplets on and off by cycling the temperature above and below the DNA melting temperature of T_m=50° C. At room temperature, complementary emulsions are mixed together and diffuse to form adhesive clusters (see FIG. 1C(1)) that dissociate upon heating above T_m (see FIG. 1C(2)). Similarly, the emulsion-colloid system below T_m reveals multivalent particle-emulsion structures (see FIG. 1D(1)) that are separated upon heating (see FIG. 1D(2)).

When complementary droplets meet, their DNA strands hybridize to form double-stranded (ds) DNA. Therefore, the presence of green and red streptavidin, associated with each DNA strand, lead to yellow adhesion patches in regions of hybridization, as shown in FIGS. 2A(1)-2A(3) (dsb) and FIGS. 2B(1)-B(3) (ssb). The patches in emulsions with double stranded backbone (dsb) tethers are significantly larger than those with single stranded backbone (ssb) tethers. In both cases, the plot of the patch size as a function of the reduced radius of each pair of contacting droplets in FIG. 2C reveals a linear dependence. Thus, the measured ratio of dsb to ssb patch sizes of ˜1.6 is independent of droplet size. This result indicates that the strength of binding between droplets can be modulated in situ by adding the complementary strand to ssb functionalized emulsions. The patch size is observed to transition from the average size expected for ssb DNA interactions to the ˜1.6 times larger adhesions corresponding to dsb DNA in less than one hour, as shown in FIG. 2D.

Another way to increase the binding strength, as well as the number of patches per droplet, is to increase the DNA coverage on the droplets, C, as shown in the 3D perspective images in FIGS. 3A(1) and A(3). The patch size (θ) and the relative intensity, defined as the ratio of the intensity of fluorescence in the patch and that of the droplet surface, increase as a function of DNA coverage in FIGS. 3B and C. This capability allows one to tune the reversibility of droplet interactions and the temperature at which the structures melt.

In one preferred embodiment a model of the experimental observation is implemented by a statistical mechanic model. It is based on the assumption that binders are recruited into the contact area until the binding energy balances the energy cost upon droplet deformation and the entropy penalty due to the immobilization of the DNA tether in the patch. Consider two complementary emulsion droplets with the same radius Re and DNA surface density N₀/(4R_e²), where N₀ is the total number of DNA on the droplet. The two droplets interact to form an adhesive patch of radius r_p and deformation angle θ=r_p/R, as shown in FIG. 1A. N_β is the number of DNA binders inside the patch between two complementary droplets, and N_α the amount of free DNA strands remaining on the droplet surface. The global free energy change from the unbound to the bound state is then:ΔF−ΔE_DNAβ−2T(S_β+S_α)+E_deform−F_unbound  (1)where E_dna is the binding energy, T is the temperature, S is the entropy of binding and E_deform is the energy cost to deform the interface. Subscripts β and α indicate the binding patch region and the unbound surface, respectively. In the dilute case, S_β=−kN_βln [C_β/C₀] where C_β=N_β/πr_p²] is the surface density of DNA in the adhesive patch and C₀ is the reference concentration, which cancels out in the calculation. The binding free energy for the mobile DNA patch is estimated in mean field: ΔE_DNA,β=N_β[ΔG_DNA−TΔS_r−k_BT ln(A_wC_β)], where ΔG_DNA is the free energy of hybridization of free DNA in solution, ΔS_r is the entropy loss due to rotational constraints of hybridized DNA strands at low temperature, A_w the area in which two DNA strands can move relative to each other when hybridized and ln(A_wC_β) is the translational entropy penalty for two DNA strands bound in the patch. The deformation energy is given by E_deform=σπR_e²θ⁴/2, where σ is the surface tension. Since there are only two independent parameters in the problem, C_β and r_p, the global energy is minimized to obtain the profile of C_β(C_α, A_w) and θ(C_α, A_w):

$\begin{array}{cc}\theta \left(C_{\alpha },A_w\right)=r_p/R_e={\left[\frac{kT\left[-21n\left(1-C_{\beta }{A}_{\mathrm{strep}}\right)-C_{\beta }{A}_{\mathrm{strep}}\right]}{\sigma A_{\mathrm{strep}}}\right]}^{1/2}& \left(2\right)\end{array}$ assuming that the binding free energy ΔG_DNA, the surface tension σ, the streptavidin size A_strep and the temperature T are kept constant. Since the double stranded tether is much longer and can reach as far as ˜26 nm, further than that of the single stranded tether of ˜4.5 nm, the area of relative motion of bound DNA strands, A_w, is also much larger, estimated to be ˜2000 nm² compared to only ˜60 nm². Therefore, the dsb case loses less entropy upon binding which quantitatively explains the ˜1.6 fold larger average patch size, as shown in FIG. 2C. Moreover, the adhesion strength dependence is captured on the DNA surface coverage with only two fitting parameters: the rotational entropy loss ΔS_r=16R (where R is the gas constant) and the maximum DNA binding capability, N=12 pmol, for an emulsion sample of 30 μL. These fitting parameters are consistent with previously reported values. The agreement of the model with the experimental data for the three trends shown in FIG. 2C and FIGS. 3B and 3C illustrate validity of the model enabling further exploitation of that model for advantageous methods and systems.

The fluidity of the droplet surface enables rearrangements in bound structures and allows for the self-assembly of programmable geometries. Adhesion patches are free to diffuse despite the high binding energy of ˜20000 DNA connections in an average-sized patch with a 1.6 μm diameter. FIGS. 4A(1)-A(5) show the diffusion of droplets that are bound through DNA interactions, but free to rotate with respect to each other and thus enable exploring available configurations. A hybrid system of particles and emulsions is used to quantify the diffusion of adhesion patches (see FIGS. 4B(1)-B(4)). The beads serve as reporters for the lipid motion on the monolayer surface. To measure relative motion, two colloidal particles coated with the S DNA sequence are attached onto the surface of a S′ functionalized droplet through at least 200 DNA bonds. The mean square displacement of one bead with respect to the other (see FIG. 4B(4)) yields a diffusion constant of D˜0.012 μm²/s. This value is significantly smaller than both the diffusion of a single lipid of size ˜1 nm in a fluid model membrane with D˜1-10 μm²/s [27-30] and that of a 1 μm colloidal particle with D_particle˜0.5 μm²/s. This slow diffusion of the particle is due to the strong hydrodynamic drag of an adhesive lipid patch of radius ˜100 nm, which is expected to be two orders of magnitude lower than that of a single lipid.

Allowing the emulsions to cream to the surface assembles floppy networks of bound droplets that are organized by the specificity of the DNA bonds, as shown in FIGS. 4C(1)-C(5). Once the maximum droplet connectivity is achieved, no further rearrangements in the structure are observed. Nevertheless, the bonds continue to be mobile owing to the liquid interfaces. While the coordination number distribution of such networks can be tuned by the concentration of binders on the surface, their structure remains amorphous.

Alternatively, the complementary colloid-emulsion hybrid system, as shown in FIG. 1B, can be used as a versatile tool for self and directed-assembly. Combining the dilute emulsion (˜1000 droplet/mm²) with complementary nanoparticles at a low particle/droplet ratio of ˜5 enables the formation of linear emulsion polymers of different lengths, as shown in FIGS. 5A(1)-A(10). The linear arrangement of the droplets is induced by constraining them in a 1-D line at the edge of a tilted rectangular capillary. The binding colloids are recruited exclusively to the emulsion contacts after overnight incubation to form two patches per emulsion droplet and thus prevent branching. This leads to polymer chains with a valency of two that diffuse over time due to the mobility of the particulate joints between the droplets. They remain in a linear configuration because the particles are too small to bridge three droplets. On the other hand, at a higher particle/droplet ratio of ˜100 we observe multivalency and folding of the linear chains into compact structures over time, as shown in FIGS. 5B(1)-B(5). In that case, the binding energy of the multivalent particles and their connectivity is sufficient to arrest the resulting structure. The geometry of the final compact structure depends on the valency, the number of droplets in the cluster, as well as the kinetic pathway, and ranges from triangular lattices to flowers, as shown in FIG. 5C(1)-C(6). These structures lower their energy by maximizing the number of colloids that occupy emulsion contacts. Unlike droplet-droplet patches, shown in FIGS. 2A(1)-2b(3), 3a(1)-a(3) and 4a(1)-a(5) and C(1)-C(5), the solid polystyrene particles are large enough to bridge two emulsion droplets without reaching their contact point. This leads to a ring structure in the adhesion zone, as shown in the confocal images of clusters in FIGS. 5D(1)-D(4). The circular arrangements maximizes the number of particles per emulsion contact and thereby minimizes the global free energy.

Self-assembly of thermal emulsion polymer chains can be achieved with programmable droplet interactions using DNA interactions, cadherins and selected nanoparticles. Controlling the number of binders and the length of the chain one can obtain divalent, trivalent and multivalent structures. In addition, the mobility of adhesive patches within these structures allows them to evolve into geometries that are governed by the underlying free energy landscape. Furthermore, such interactions allow one to program the shape of the free energy landscape through the control of bond specificity, strength, flexibility and valency. This system promises to become a highly advantageous system and method for producing products by directed self-assembly because it has the potential of building intelligently designed materials, such as colloidal crystals or artificial self-replicating materials, with no external inputs. These materials and methods can be used for a variety of commercial applications including, without limitation, controlled formation of personal care products, food processing, skin creams, pharmaceutical products, foods and animal feedstocks.

Various embodiments are described in the general context of method steps, which may be implemented in one embodiment by a computer 100 having an embedded program in a non-transitory storage medium 200 and including computer-executable instructions, such as program code, executed by the computer 100 in networked environments or in the cloud 300. In FIG. 6 is shown a block diagram illustrating such a system. Generally, program modules include routines, programs, objects, components, data structures, etc. that perform particular tasks or implement particular abstract data types. Computer-executable instructions, associated data structures, and program modules represent examples of program code for executing steps of the methods disclosed herein. The particular sequence of such executable instructions or associated data structures represents examples of corresponding acts for implementing the functions described in such steps.

Software and web implementations of the present invention could be accomplished with standard programming techniques with rule based logic and other logic to accomplish the various database searching steps, correlation steps, comparison steps and decision steps. It should also be noted that the words “component” and “module,” as used herein and in the claims, are intended to encompass implementations using one or more lines of software code, and/or hardware implementations, and/or equipment for receiving manual inputs.

The following non-limiting Examples illustrate selected aspects of the invention.

EXAMPLE I

Synthesis of Biomimetic Emulsion:

The protocol for the emulsion preparation is a conventional, well known methodology. The oil droplets are stabilized with egg L-α-phosphatidylcholine (EPC) lipids and the DSPE-PEG(2000) biotinylated lipids from Avanti Polar Lipids at a molar ratio of 92:8, respectively. After cooling to room temperature the lipid containing oil (10 mL) can be emulsified in two different buffers to reach different droplet sizes. Athermal droplets are obtained when emulsified in a buffer containing 5 mMSDS and w_t=18% dextran and sheared at 22 rpm in the narrow gap coquette mixer. Smaller droplets sizes are obtained with a buffer containing 5 mMSDS and w_t=4.5% alginate and a shearing rate of 30 rpm. The resulting emulsions are stable over several weeks at 4° C.

DNA coated emulsion preparation: The emulsion is first coated with two different fluorescent streptavidins: Alexa Fluor 488® and Alexa Fluor 633® streptavidin (Molecular Probes). 100 μL of emulsion is mixed with 10 μL of 1 mg/mL streptavidin and 300 μL of buffer containing 2 mM Tris pH=7 and 1 mM SDS. The solution is incubated 1 h at 4° C. and then washed twice with the same buffer, before a final wash in the DNA buffer (1 mM SDS, 5 mM PBS, 4 mM MgCl2). The DNA can then be added to the streptavidin coated emulsion: 10 μL of 50 μM DNA is added to the solution and incubated 1 h at room temperature. The remaining unbound DNA is finally washed three times in the DNA buffer. One type of fluorescent streptavidin is specifically associated with one DNA strand in order to distinguish them by microscopy.

DNA Constructs:

The S strand has a sequence: 5′-BiotinTEG-49bp backbone-GGATGAAGATG-3′ (SEQ ID NO: 1);

The S′ strand has a sequence: 5′-BiotinTEG-49bp backbone-CATCTTCATCC-3″ (SEQ ID NO: 2);

The CS strand complementary to the backbone: 5′-TCG TAA TGA AAG GCA GGG CTC TCT GGA TTG ACT GTG CGAAGG GTA GCG AT-3″ (SEQ ID NO: 3)

TEG: Tetra-Ethylene Glycol

Confocal microscope: The samples are imaged using a fast scanning confocal

microscope (Leica TCS SPS 11).

Light microscope with a temperature stage: A Leica DMRXA microscope with Qimaging Retiga 1300 camera is used to obtain microscopic images. A temperature stage is built on the light microscope to provide fast in-situ temperature control. Briefly, 1000Ω ITO glass is placed on a 3 mm thick copper plate, two ends of which are connected to peltiers (2.5 cm by 2.5 cm) then to a thermal sink with constant temperature. With LakeShore DRC 93C Temperature Controller and LakeShore PT-111 temperature sensor, we are able to control and detect the temperature with <0.5° C. relative error.

EXAMPLE II

Consider two interacting droplets of the same radius R and coated with complementary strands of DN, when two DNA strands from opposite surfaces bind together, they gain binding energy but lose entropy due to the spatial constraint into the patch. The system also endures a deformation energy cost to allow the droplets to deform and form the adhesion patch area.

First we use the simple ‘lattice model’ (or ‘box model’) to get the entropy of molecules in a non-dilute solution. Given the total area S_total, the molecule (streptavidin) size A_strep, the number of streptavidins N and the total number of sites available on the droplet surface N₀=S_total A_strep where (move down to be on same line)

$N_0=\frac{S_{\mathrm{total}}}{A_{\mathrm{strep}}},$ the total number of configurations reads:

$\begin{array}{cc}\Omega =\frac{N_0!}{\left(N_0-N\right)!N!}& \left(3\right)\end{array}$

Using the Stirling's formula the entropy is approximated to:

$\begin{array}{cc}S=k\ln \left(\Omega \right)=k\left[N\ln \left(\frac{N_0-N}{N}\right)+N_0\ln \left(\frac{N_0}{N_0-N}\right)\right]& \left(4\right)\end{array}$

With the binding energy for a pair of DNA sticky ends (move down on same line) ΔG_DNA=ΔF_DNA−TΔS_DNA, the deformation energy of the droplet [2] and the entropy term derived in Equation (2), we can write down the global free energy difference between the bound state and the non-interacting droplets state as follows:ΔF=E_DNA,β−2(TS_β+TS_α)+E_deformation−F_unbound  (5)

The subscript β refers to the adhesive patch region while α refers to the non-interacting region on the rest of the droplet surface. Since two droplets interact to form a patch, the entropy term has to be taken into account twice which justifies the prefactor.

The energy terms in Equation (3) can be written:

$\begin{array}{cc}{E}_{\mathrm{DN}A,\beta }=-\mathrm{kT}\ln \left[{\left(1+\exp \left(-\frac{\Delta G_{\mathrm{DN}A}-T\Delta S_r-T\ln \left(A_w{C}_{\beta }\right)}{\mathrm{kT}}\right)\right)}^{N_a}-1\right]& \left(6\right)\end{array}$

The entropy and the deformation energy in Equation (3) can be written:

$\begin{array}{cc}{S}_{\beta }=k\left[N_{\beta }\ln \left(\frac{N_{\beta 0}-N_{\beta }}{N_{\beta }}\right)+N_{\mathrm{\beta 0}}\ln \left(\frac{N_{\mathrm{\beta 0}}}{N_{\mathrm{\beta 0}}-N_{\beta }}\right)\right]S_{\alpha }=k\left[N_{\alpha }\ln \left(\frac{N_{\mathrm{\alpha 0}}-N_{\alpha }}{N_{\alpha }}\right)+N_{\mathrm{\alpha 0}}\ln \left(\frac{N_{\mathrm{\alpha 0}}}{N_{\mathrm{\alpha 0}}-N_{\alpha }}\right)\right]E_{\mathrm{deformation}}=\frac{1}{2}\mathrm{\sigma \pi }{R}^2{\theta }^4=\frac{1}{2}\mathrm{\sigma \pi }\frac{r_P^4}{R^2}& \left(7\right)\end{array}$

Where r_p is the radius of the enriched patch; θ is defined as the deformation angle r_p/R; σ is the surface tension of the emulsion; N_+N_=N gives the total number of streptavidins, N_β of them being in the binding patch; (try to move all these inserts to be level with line) N_αθ+N_βθ=N_θ=4πR²/A_strep gives the total number of biotin sites on a emulsion, while N_βθ=πr_p²/A_strep is the number of sites available in the patch area; C_β=N_β/(πr_p²) is the concentration of streptavidin in the patch; A_w is the area over which two bound DNA strands could move relative to each other while remaining hybridized; ΔS_r and k ln(A_wC_β) are the configurational entropy lost due to rotational and translational confinement of hybridized DNA sticky ends, respectively.

We now minimize this global free energy ΔF with respect to two independent parameters in the equations: N_βθ and N_β. Note that we could conversely use the two independent parameters C_β and r_p instead, which would result in the same equations.

The first equation leads to the chemical potential equilibrium. In the strong binding case where ΔF_DNA−TΔS_DNA−TΔS_r at least a few kT, this first equation can be simplified as follows:

$\begin{array}{cc}\Delta G_{\mathrm{DNA}}-T\Delta S_r-\mathrm{kT}-\mathrm{kT}\ln \left(A_w{C}_{\beta }\right)-2T\left[k\ln \left(\frac{N_{\mathrm{\beta 0}}-N_{\beta }}{N_{\beta }}\right)-k\ln \left(\frac{N_{\mathrm{\alpha 0}}-N_{\alpha }}{N_{\alpha }}\right)\right]=0& \left(8\right)\end{array}$ The second equation reads:

$\begin{array}{cc}\mathrm{kT}\frac{N_{\beta }}{N_{\mathrm{\beta 0}}}-2\mathrm{kT}\left[\ln \left(\frac{N_{\mathrm{\beta 0}}}{N_{\mathrm{\beta 0}}-N_{\beta }}\right)-\ln \left(\frac{N_{\mathrm{\alpha 0}}}{N_{\mathrm{\alpha 0}}-N_{\alpha }}\right)\right]+\frac{\sigma N_{\mathrm{\beta 0}}{A}_{\mathrm{strep}}^2}{\pi R^2}=0& \left(9\right)\end{array}$

The resulting N_β and N_βθ, directly leading to values of r_p and C_β/C_α, can be solved numerically which allows the comparison with our experimental values for the patch size r_p and contrast C_β/C_α. These analytical solutions are obtained under the approximation that there is an infinite dilute reservoir with a constant supply N_α/N_αθ=Const=d. This approximation is reasonable for our experimental condition, since N_α/N_αθ<0.1 and the relative change in C_α is less than 10%, even with the most enrichment condition.

As a result, the approximate solution to Equation (8) is:

$\begin{array}{cc}\frac{N_{\beta }}{N_{\mathrm{\beta 0}}}=\frac{2{\mathrm{cd}}^2+e^{a\text{/}b}-e^{\frac{a}{2b}}\sqrt{4{\mathrm{cd}}^2+e^{a\text{/}b}}}{2{\mathrm{cd}}^2}& \left(10\right)\end{array}$

Where α=ΔG_DNA−TΔS_r−kT, b=kT, c=A_w/A_strep. This expression is directly linked to the measured patch intensity contrast

$C_{\beta }\text{/}{C}_{\alpha }=\frac{N_{\beta }}{N_{\mathrm{\beta 0}}}·\frac{N_{\mathrm{\alpha 0}}}{N_{\alpha }}=\frac{N_{\beta }}{d·N_{\mathrm{\beta 0}}}.$ Since we know the relation N_βθ=πr_p²/A_strep, Equation (9) directly gives us:

$\begin{array}{cc}{r}_p^2=R^2\frac{\mathrm{kT}\left[2\ln \left(\frac{N_{\mathrm{\beta 0}}}{N_{\mathrm{\beta 0}}-N_{\beta }}\right)-\frac{N_{\beta }}{N_{\mathrm{\beta 0}}}\right]}{\sigma A_{\mathrm{strep}}}& \left(11\right)\\ \theta =\sqrt{\frac{\mathrm{kT}\left[2\ln \left(\frac{N_{\mathrm{\beta 0}}}{N_{\mathrm{\beta 0}}-N_{\beta }}\right)-\frac{N_{\beta }}{N_{\mathrm{\beta 0}}}\right]}{\sigma A_{\mathrm{strep}}}}& \left(12\right)\end{array}$

We can now compare our experimental values to the ones found analytically here.

For the DNA sequence used in the experiments, ΔG_DNA=ΔF_DNA−TΔS_DNA is ≈−20 kT at room temperature, and the experimental value for TΔS_r are −14.6 kT for the double-stranded backbone DNA and ˜14.8 kT for the single-stranded backbone one. We therefore use the same fitting parameters for both the ssb and dsb case. A_w/A_strep≅34 dsb while it is only A_w/A_strep≅1 for ssb DNA. This discrepancy is due to the different rigidities of the DNA strands: double stranded DNA is more rigid and rod-like and can thus reach a large number of strands on the opposite surface, whereas single stranded DNA behaves as a very exible polymer in our buffer conditions, with a persistent length of ˜2 nm leading to a smaller end-to-end distance of ≈6 nm, σ≅15 mN/m for phospholipid emulsions co-stabilized with 1 mM SDS. With a streptavidin size of A_strep=60 nm² and the initial streptavidin surface concentration of 1400/μm², this leads to d_min≅0.09.

Experimentally we vary the DNA surface density d=N_DNA/N_max by changing the amount of DNA introduced in the system N_DNA˜1 pmol, 2 pmol, 4 pmol, 8 pmol, 20 pmol or 80 pmol. An ideal emulsion packing of 30 μL, as used in this experiment, can bind up to N_max˜30 pmol of DNA. Nevertheless the experiments require two washing steps of the emulsions before DNA conjugation, which is suspected to significantly reduce this number.

As a result, all the data in FIG. 8c and FIGS. 9b and 9c can be fitted with only two fitting parameters: ΔS_p=−16R and N_max=12 pmol.

Polydisperse Emulsion Droplets Interaction

In the approximation of infinite reservoir, the only radius dependent term in the above set of equations are from the deformation energy. The deformation energy of the emulsions should be corrected as:

$E_{\mathrm{deformation}}=\frac{1}{4}\mathrm{\sigma \pi }\frac{r_p^4}{R_1^2}+\frac{1}{4}\mathrm{\sigma \pi }\frac{r_p^4}{R_2^2}$ with different radius of contacting emulsions, R1 and R2, in the lowest order approximation. We define a square-averaged radius

$<R>=\sqrt{\frac{2R_1^2{R}_2^2}{R_1^2+R_2^2}},$ so that

$E_{\mathrm{deformation}}=\frac{1}{2}\mathrm{\sigma \pi }\frac{r_p^4}{<R{>}^2}$ and we can use all the equations in the previous section replacing R with <R>, as plotted or used in main text FIGS. 2c and d.

Linear Regression, Additional Geometry and Diffraction Limit

We fit d_p as a function of <R> with a simple linear regression relation rather than a line d_p=θR going through the origin as suggested by our model. The origin of this choice lies in geometrical arguments. Indeed the DNA constructs can be stretched, which leads to enrichment outside of the geometrically predicted adhesion patch. This additional area leads to a geometric factor

$\frac{\Delta L}{\theta }$ contributing to me paten size. Both double-stranded and single-stranded DNA can extend up to ΔL˜12 nm as estimated respectively from conventional teachings and a worm like chain model. This leads to an entropy loss of ˜2−3 kT which reduces DNA concentration by half

As a result the model is modified to include this additional term:

$\begin{array}{cc}{d}_p=2\theta R+\frac{\Delta L}{\theta }& \left(13\right)\end{array}$

The estimate for the respective intersections for dsb and ssb DNA give the values of ΔL/θ˜80 nm and ˜160 nm, which are smaller than the experimental values of ˜150 nm (dsb) and ˜210 nm (ssb) that are certainly fixed by our experimental diffraction limit of ˜150 nm.

Nevertheless, the fitting curves with either

$d_p=2\theta R+\frac{\Delta L}{\theta }\mathrm{or}{d}_p=2\theta R$ are similar with the data shown in FIG. 2C and the fitting parameter ΔS_p only changes by ˜5% to fit the data. Therefore the limits of our experimental accuracy do not allow us to discriminate between the two relationships.

The foregoing description of illustrative embodiments has been presented for purposes of illustration and of description. It is not intended to be exhaustive or limiting with respect to the precise form disclosed, and modifications and variations are possible in light of the above teachings or may be acquired from practice of the disclosed embodiments. It is intended that the scope of the invention be defined by the claims appended hereto and their equivalents.

Claims

1. A method of self-assembly to form end product, comprising the steps of: providing a first emulsion droplet for assembly;providing a second emulsion droplet for assembly;providing a first linker comprising a colloidal nanoparticle which can be coupled to the first emulsion droplet;forming a first patch, comprising the first linker, between the first emulsion droplet and the second emulsion droplet.

2. The method as defined in claim 1 wherein pseudo-bonding valency is established for the linker components.

3. The method as defined in claim 2 wherein a valency of 2 enables formation of flexible polymers of emulsion droplets.

4. The method as defined in claim 2 wherein a pseudo-valency of 4 enables formation of rigid polymer networks.

5. The method as defined in claim 1 wherein the first linker further comprises a plurality of single stranded DNA, attached to the colloidal nanoparticle, having a first sequence and further wherein the first emulsion droplet and the second emulsion droplet each have attached thereto a plurality of single stranded DNA having a second sequence, the second sequence and the first sequence being complementary.

6. The method as defined in claim 1 further including the step of mixing the end production with a personal care product matrix.

7. The method as defined in claim 1 further including the step of mixing the end product with a food product matrix.

8. The method as defined in claim 1 wherein the end production consists of an amorphous material having selectable rheological properties, thereby enabling processing of the end product to operate consumer products having desired properties.

9. The method as defined in claim 1 further including the step of cycling temperature of the end product.

10. The method as defined in claim 1 wherein the linker component comprises DNA strands further including the step of increased coverage of the DNA strands on the first component, thereby enabling tuning of reversibility of self-assembly.

11. A method of self-assembly to form end product, comprising the steps of: providing an emulsion having a plurality of emulsion droplets and a plurality of linkers, wherein each linker comprises a nanoparticle having attached thereto a first DNA strand and a second DNA strand and the concentration of the plurality of linkers is selected to provide a desired valency;providing a first linker and a second linker of the plurality of linkers, each of the first linker and the second linker engagable with the plurality of emulsion droplets and engageable with each other;engaging the first DNA strand with a first emulsion droplet of the plurality of droplets;engaging the second DNA strand with a second emulsion droplet of the plurality of droplets;forming a patch between the first emulsion droplet and the second emulsion droplet comprising the first DNA strand, the nanoparticle, and the second DNA strand;engaging the first emulsion droplet with additional emulsion droplets to achieve the desired valency andcreaming the emulsion forming a floppy network of the plurality of emulsion droplets interconnected by a plurality of patches.

12. The method of claim 11, wherein the emulsion has a ratio of 5 to 1 of the plurality of linkers to the plurality of droplets, whereby a linear assembly is formed.

13. The method of claim 11, wherein the emulsion has a ratio of 100 to 1 of the plurality of linkers to the plurality of emulsion droplets, whereby a folded assembly is formed.

14. The method of claim 11, further comprising selecting a valency for the floppy network.