Self-assembly of multicomponent structures in and out of equilibrium

Mechanism of DNA origami folding elucidated by mesoscopic simulations

Many experimental and computational efforts have sought to understand DNA origami folding, but the time and length scales of this process pose significant challenges. Here, we present a mesoscopic model that uses a switchable force field to capture the behavior of single- and double-stranded DNA motifs and transitions between them, allowing us to simulate the folding of DNA origami up to several kilobases in size. Brownian dynamics simulations of small structures reveal a hierarchical folding process involving zipping into a partially folded precursor followed by crystallization into the final structure. We elucidate the effects of various design choices on folding order and kinetics. Larger structures are found to exhibit heterogeneous staple incorporation kinetics and frequent trapping in metastable states, as opposed to more accessible structures which exhibit first-order kinetics and virtually defect-free folding. This model opens an avenue to better understand and design DNA nanostructures for improved yield and folding performance.

Nonequilibrium design strategies for functional colloidal assemblies

We use a nonequilibrium variational principle to optimize the steady-state, shear-induced interconversion of self-assembled nanoclusters of DNA-coated colloids. Employing this principle within a stochastic optimization algorithm allows us to identify design strategies for functional materials. We find that far-from-equilibrium shear flow can significantly enhance the flux between specific colloidal states by decoupling trade-offs between stability and reactivity required by systems in equilibrium. For isolated nanoclusters, we find nonequilibrium strategies for amplifying transition rates by coupling a given reaction coordinate to the background shear flow. We also find that shear flow can be made to selectively break detailed balance and maximize probability currents by coupling orientational degrees of freedom to conformational transitions. For a microphase consisting of many nanoclusters, we study the flux of colloids hopping between clusters. We find that a shear flow can amplify the flux without a proportional compromise on the microphase structure. This approach provides a general means of uncovering design principles for nanoscale, autonomous, functional materials driven far from equilibrium.

Mechanistic Study of Amphiphilic-Assisted Self-Assembled Cadmium Sulfide Quantum Dots into 3D Superstructures

Self-assembling of nanoparticles into complex superstructures is very challenging, which usually depends on postorganizing techniques or pre-existing templates such as polypeptide chains or DNA or external stimulus. Such self-assembled processes typically lead to close-packed structures. Here, it has been demonstrated that under carefully template-free reaction conditions CdS quantum dots (QDs) could be synthesized and simultaneously self-assembled into complex superstructures without compromising individual QD properties. The superstructures of CdS QDs attained by the chemical-based method demonstrate Stokes-shifted photoluminescence (PL) from trap states. Remarkably, the PL decay of superstructures exhibits a single-exponential feature. This behavior is unusual for the synthesized superstructures, indicating that the trap states are restricted to a narrow range. The growth mechanism of these superstructures is explained through the formation of liquid crystal phases (LCPs) with the help of a small-angle X-ray scattering (SAXS) analysis.

A universal method for analyzing copolymer growth

Polymers consisting of more than one type of monomer, known as copolymers, are vital to both living and synthetic systems. Copolymerization has been studied theoretically in a number of contexts, often by considering a Markov process in which monomers are added or removed from the growing tip of a long copolymer. To date, the analysis of the most general models of this class has necessitated simulation. We present a general method for analyzing such processes without resorting to simulation. Our method can be applied to models with an arbitrary network of sub-steps prior to addition or removal of a monomer, including non-equilibrium kinetic proofreading cycles. Moreover, the approach allows for a dependency of addition and removal reactions on the neighboring site in the copolymer and thermodynamically self-consistent models in which all steps are assumed to be microscopically reversible. Using our approach, thermodynamic quantities such as chemical work; kinetic quantities such as time taken to grow; and statistical quantities such as the distribution of monomer types in the growing copolymer can be directly derived either analytically or numerically from the model definition.

Catalystlike role of impurities in speeding layer-by-layer growth

Molecular self-assembly is usually done at low supersaturation, leading to low rates of growth, in order to allow time for binding mistakes to anneal. However, such conditions can lead to prohibitively long assembly times where growth proceeds by the slow nucleation of successive layers. Here we use a lattice model of molecular self-assembly to show that growth in this regime can be sped up by impurities, which lower the free-energy cost of layer nucleation. Under certain conditions impurities behave almost as a catalyst in that they are present at high concentration at the surface of the assembling structure, but at low concentration in the bulk of the assembled structure. Extrapolation of our numerics using simple analytic arguments suggests that this mechanism can reduce growth times by orders of magnitude in parameter regimes applicable to molecular systems.

Typical Stochastic Paths in the Transient Assembly of Fibrous Materials

When chemically fueled, molecular self-assembly can sustain dynamic aggregates of polymeric fibers-hydrogels-with tunable properties. If the fuel supply is finite, the hydrogel is transient, as competing reactions switch molecular subunits between active and inactive states, drive fiber growth and collapse, and dissipate energy. Because the process is away from equilibrium, the structure and mechanical properties can reflect the history of preparation. As a result, the formation of these active materials is not readily susceptible to a statistical treatment in which the configuration and properties of the molecular building blocks specify the resulting material structure. Here, we illustrate a stochastic-thermodynamic and information-theoretic framework for this purpose and apply it to these self-annihilating materials. Among the possible paths, the framework variationally identifies those that are typical-loosely, the minimum number with the majority of the probability. We derive these paths from computer simulations of experimentally-informed stochastic chemical kinetics and a physical kinetics model for the growth of an active hydrogel. The model reproduces features observed by confocal microscopy, including the fiber length, lifetime, and abundance as well as the observation of fast fiber growth and stochastic fiber collapse. The typical mesoscopic paths we extract are less than 0.23% of those possible, but they accurately reproduce material properties such as mean fiber length.

Sequence-Regulated Supracolloidal Copolymers via Copolymerization-Like Coassembly of Binary Mixtures of Patchy Nanoparticles

Synthetic copolymers of molecular systems serve as an inspiration for creation of one-dimensional copolymer-like superstructures via coassembly of anisometric nanoparticles. In contrast to the covalent and molecular copolymers, the details of formation mechanisms of copolymer-like superstructures, as well as the factors determining their length and the sequences of arranged nanoparticles, are still poorly understood. Herein, we propose a joint theoretical-computational framework to probe into the coassembly mechanism and kinetics of binary mixtures of patchy nanoparticles. By applying the coarse-grained molecular dynamics simulations, it is demonstrated that the coassembly of patchy nanoparticles markedly resembles many aspects of molecular step-growth copolymerization, and the sequences of nanoparticles inside the copolymer-like superstructures can be finely regulated by the relative activity and the initial ingredient of patchy nanoparticles as well as the coassembly strategy. A quantitatively copolymerization-like model is developed to account for the coassembly kinetics of patchy nanoparticles and the sequence distribution of arranged nanoparticles, all governed by the elaborate design of lower-level building units. The jointly theoretical and simulated studies offer mechanistic insights into the copolymerization-like kinetics and the sequence prediction for the coassembly of binary mixtures of patchy nanoparticles, paving the way toward the rational design of copolymer-like superstructures with various sequences and functionalities.

Nonequilibrium correlations in minimal dynamical models of polymer copying

Living systems produce "persistent" copies of information-carrying polymers, in which template and copy sequences remain correlated after physically decoupling. We identify a general measure of the thermodynamic efficiency with which these nonequilibrium states are created and analyze the accuracy and efficiency of a family of dynamical models that produce persistent copies. For the weakest chemical driving, when polymer growth occurs in equilibrium, both the copy accuracy and, more surprisingly, the efficiency vanish. At higher driving strengths, accuracy and efficiency both increase, with efficiency showing one or more peaks at moderate driving. Correlations generated within the copy sequence, as well as between template and copy, store additional free energy in the copied polymer and limit the single-site accuracy for a given chemical work input. Our results provide insight into the design of natural self-replicating systems and can aid the design of synthetic replicators.

Topological localization in out-of-equilibrium dissipative systems

In this paper, we report that notions of topological protection can be applied to stationary configurations that are driven far from equilibrium by active, dissipative processes. We consider two physically disparate systems: stochastic networks governed by microscopic single-particle dynamics, and collections of driven interacting particles described by coarse-grained hydrodynamic theory. We derive our results by mapping to well-known electronic models and exploiting the resulting correspondence between a bulk topological number and the spectrum of dissipative modes localized at the boundary. For the Markov networks, we report a general procedure to uncover the topological properties in terms of the transition rates. For the active fluid on a substrate, we introduce a topological interpretation of fluid dissipative modes at the edge. In both cases, the presence of dissipative couplings to the environment that break time-reversal symmetry are crucial to ensuring topological protection. These examples constitute proof of principle that notions of topological protection do indeed extend to dissipative processes operating out of equilibrium. Such topologically robust boundary modes have implications for both biological and synthetic systems.

Theory of amyloid fibril nucleation from folded proteins

We present a theoretical model for the nucleation of amyloid fibrils. In our model we use helix-coil theory to describe the equilibrium between a soluble native state and an aggregation-prone unfolded state. We then extend the theory to include oligomers with β-sheet cores and calculate the free energy of these states using estimates for the energies of H-bonds, steric zipper interactions, and the conformational entropy cost of forming secondary structure. We find that states with fewer than ~10 β-strands are unstable relative to the dissociated state and three β-strands is the highest free energy state. We then use a modified version of Classical Nucleation Theory to compute the nucleation rate of fibrils from a supersaturated solution of monomers, dimers, and trimers. The nucleation rate has a non-monotonic dependence on denaturant concentration reflecting the competing effects of destabilizing the fibril and increasing the concentration of unfolded monomers. We estimate heterogeneous nucleation rates and discuss the application of our model to secondary nucleation.

Design principles for nonequilibrium self-assembly

We consider an important class of self-assembly problems, and using the formalism of stochastic thermodynamics, we derive a set of design principles for growing controlled assemblies far from equilibrium. The design principles constrain the set of configurations that can be obtained under nonequilibrium conditions. Our central result provides intuition for how equilibrium self-assembly landscapes are modified under finite nonequilibrium drive.

Self-organized architectures from assorted DNA-framed nanoparticles

The science of self-assembly has undergone a radical shift from asking questions about why individual components self-organize into ordered structures, to manipulating the resultant order. However, the quest for far-reaching nanomanufacturing requires addressing an even more challenging question: how to form nanoparticle (NP) structures with designed architectures without explicitly prescribing particle positions. Here we report an assembly concept in which building instructions are embedded into NPs via DNA frames. The integration of NPs and DNA origami frames enables the fabrication of NPs with designed anisotropic and selective interactions. Using a pre-defined set of different DNA-framed NPs, we show it is possible to design diverse planar architectures, which include periodic structures and shaped meso-objects that spontaneously emerge on mixing of the different topological types of NP. Even objects of non-trivial shapes, such as a nanoscale model of Leonardo da Vinci's Vitruvian Man, can be self-assembled successfully.

Landscape of kinetically trapped binary assemblies

For two-component assemblies, an inherent structure diagram (ISD) is the relationship between set inter-subunit energies and the types of kinetic traps (inherent structures) one may obtain from those energies. It has recently been shown that two-component ISDs are apportioned into regions or plateaux within which inherent structures display uniform features (e.g., stoichometries and morphologies). Interestingly, structures from one of the plateaux were also found to be robust outcomes of one type of non-equilibrium growth, which indicates the usefulness of the two-component ISD in predicting outcomes of some types of far-from-equilibrium growth. However, little is known as to how the ISD is apportioned into distinct plateaux. Also, while each plateau displays classes of structures that are morphologically distinct, little is known about the source of these distinct morphologies. This article outlines an analytic treatment of the two-component ISD and shows that the manner in which any ISD is apportioned arises from a single unitless order parameter. Additionally, the analytical framework allows for the characterization of local properties of the trapped structures within each ISD plateau. This work may prove to be useful in the design of novel classes of robust nonequilibrium assemblies.

Heterogeneity of functional groups in a metal-organic framework displays magic number ratios

Multiple organic functionalities can now be apportioned into nanoscale domains within a metal-coordinated framework, posing the following question: how do we control the resulting combination of "heterogeneity and order"? Here, we report the creation of a metal-organic framework, MOF-2000, whose two component types are incorporated in a 2:1 ratio, even when the ratio of component types in the starting solution is varied by an order of magnitude. Statistical mechanical modeling suggests that this robust 2:1 ratio has a nonequilibrium origin, resulting from kinetic trapping of component types during framework growth. Our simulations show how other "magic number" ratios of components can be obtained by modulating the topology of a framework and the noncovalent interactions between component types, a finding that may aid the rational design of functional multicomponent materials.

An assembly funnel makes biomolecular complex assembly efficient

Like protein folding and crystallization, the self-assembly of complexes is a fundamental form of biomolecular organization. While the number of methods for creating synthetic complexes is growing rapidly, most require empirical tuning of assembly conditions and/or produce low yields. We use coarse-grained simulations of the assembly kinetics of complexes to identify generic limitations on yields that arise because of the many simultaneous interactions allowed between the components and intermediates of a complex. Efficient assembly occurs when nucleation is fast and growth pathways are few, i.e. when there is an assembly "funnel". For typical complexes, an assembly funnel occurs in a narrow window of conditions whose location is highly complex specific. However, by redesigning the components this window can be drastically broadened, so that complexes can form quickly across many conditions. The generality of this approach suggests assembly funnel design as a foundational strategy for robust biomolecular complex synthesis.

Growth of equilibrium structures built from a large number of distinct component types

We use simple analytic arguments and lattice-based computer simulations to study the growth of structures made from a large number of distinct component types. Components possess 'designed' interactions, chosen to stabilize an equilibrium target structure in which each component type has a defined spatial position, as well as 'undesigned' interactions that allow components to bind in a compositionally-disordered way. We find that high-fidelity growth of the equilibrium target structure can happen in the presence of substantial attractive undesigned interactions, as long as the energy scale of the set of designed interactions is chosen appropriately. This observation may help explain why equilibrium DNA 'brick' structures self-assemble even if undesigned interactions are not suppressed [Ke et al. Science, 338, 1177, (2012)]. We also find that high-fidelity growth of the target structure is most probable when designed interactions are drawn from a distribution that is as narrow as possible. We use this result to suggest how to choose complementary DNA sequences in order to maximize the fidelity of multicomponent self-assembly mediated by DNA. We also comment on the prospect of growing macroscopic structures in this manner.

Thermodynamics versus Kinetics Dichotomy in the Linear Self-Assembly of Mixed Nanoblocks

We report classical and replica exchange molecular dynamics simulations that establish the mechanisms underpinning the growth kinetics of a binary mix of nanorings that form striped nanotubes via self-assembly. A step-growth coalescence model captures the growth process of the nanotubes, which suggests that high aspect ratio nanostructures can grow by obeying the universal laws of self-similar coarsening, contrary to systems that grow through nucleation and elongation. Notably, striped patterns do not depend on specific growth mechanisms, but are governed by tempering conditions that control the likelihood of depropagation and fragmentation.

Self-assembly at a nonequilibrium critical point

We use analytic theory and computer simulation to study patterns formed during the growth of two-component assemblies in two and three dimensions. We show that these patterns undergo a nonequilibrium phase transition, at a particular growth rate, between mixed and demixed arrangements of component types. This finding suggests that principles of nonequilibrium statistical mechanics can be used to predict the outcome of multicomponent self-assembly, and suggests an experimental route to the self-assembly of multicomponent structures of a qualitatively defined nature.

Kinetic theory of amyloid fibril templating

The growth of amyloid fibrils requires a disordered or partially unfolded protein to bind to the fibril and adapt the same conformation and alignment established by the fibril template. Since the H-bonds stabilizing the fibril are interchangeable, it is inevitable that H-bonds form between incorrect pairs of amino acids which are either incorporated into the fibril as defects or must be broken before the correct alignment can be found. This process is modeled by mapping the formation and breakage of H-bonds to a one-dimensional random walk. The resulting microscopic model of fibril growth is governed by two timescales: the diffusion time of the monomeric proteins, and the time required for incorrectly bound proteins to unbind from the fibril. The theory predicts that the Arrhenius behavior observed in experiments is due to off-pathway states rather than an on-pathway transition state. The predicted growth rates are in qualitative agreement with experiments on insulin fibril growth rates as a function of protein concentration, denaturant concentration, and temperature. These results suggest a templating mechanism where steric clashes due to a single mis-aligned molecule prevent the binding of additional molecules.