The promise of self-assembly to enable the bottom-up formation of materials with prescribed architectures and functions has driven intensive efforts to uncover rational design principles for maximizing the yield of a target structure. Yet, despite many successful examples of self-assembly, ensuring kinetic accessibility of the target structure remains an unsolved problem in many systems. In particular, long-lived kinetic traps can result in assembly times that vastly exceed experimentally accessible timescales. One proposed solution is to design non-equilibrium assembly protocols in which system parameters change over time to avoid such kinetic traps. Here, we develop a framework to combine Markov state model (MSM) analysis with optimal control theory to compute a time-dependent protocol that maximizes the yield of the target structure at a finite time. We present an adjoint-based gradient descent method that, in conjunction with MSMs for a system as a function of its control parameters, enables efficiently optimizing the assembly protocol. We also describe an interpolation approach to significantly reduce the number of simulations required to construct the MSMs. We demonstrate our approach with two examples; a simple semi-analytic model for the folding of a polymer of colloidal particles, and a more complex model for capsid assembly. Our results show that optimizing time-dependent protocols can achieve significant improvements in the yields of selected structures, including equilibrium free energy minima, long-lived metastable structures, and transient states.
The considerable interest in two-dimensional (2D) materials and complex molecular topologies calls for a robust experimental system for single-molecule studies. In this work, we study the equilibrium properties and deformation response of a complex DNA structure called a kinetoplast, a 2D network of thousands of linked rings akin to molecular chainmail. Examined in good solvent conditions, kinetoplasts appear as a wrinkled hemispherical sheet. The conformation of each kinetoplast is dictated by its network topology, giving it a unique shape, which undergoes small-amplitude thermal fluctuations at subsecond timescales, with a wide separation between fluctuation and diffusion timescales. They deform elastically when weakly confined and swell to their equilibrium dimensions when the confinement is released. We hope that, in the same way that linear DNA became a canonical model system on the first investigations of its polymer-like behavior, kinetoplasts can serve that role for 2D and catenated polymer systems.
more » « less- Award ID(s):
- 1936696
- NSF-PAR ID:
- 10126567
- Publisher / Repository:
- Proceedings of the National Academy of Sciences
- Date Published:
- Journal Name:
- Proceedings of the National Academy of Sciences
- Volume:
- 117
- Issue:
- 1
- ISSN:
- 0027-8424
- Page Range / eLocation ID:
- p. 121-127
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
-
more » « less
The classic picture of soft material mechanics is that of rubber elasticity, in which material modulus is related to the entropic elasticity of flexible polymeric linkers. The rubber model, however, largely ignores the role of valence (i.e., the number of network chains emanating from a junction). Recent work predicts that valence, and particularly the Maxwell isostatic point, plays a key role in determining the mechanics of semiflexible polymer networks. Here, we report a series of experiments confirming the prominent role of valence in determining the mechanics of a model system. The system is based on DNA nanostars (DNAns): multiarmed, self-assembled nanostructures that form thermoreversible equilibrium gels through base pair-controlled cross-linking. We measure the linear and nonlinear elastic properties of these gels as a function of DNAns arm number, f, and concentration [DNAns]. We find that, as f increases from three to six, the gel’s high-frequency plateau modulus strongly increases, and its dependence on [DNAns] transitions from nonlinear to linear. Additionally, higher-valence gels exhibit less strain hardening, indicating that they have less configurational freedom. Minimal strain hardening and linear dependence of shear modulus on concentration at high f are consistent with predictions for isostatic systems. Evident strain hardening and nonlinear concentration dependence of shear modulus suggest that the low-f networks are subisostatic and have a transient, potentially fractal percolated structure. Overall, our observations indicate that network elasticity is sensitive both to entropic elasticity of network chains and to junction valence, with an apparent isostatic point
in agreement with the Maxwell prediction. -
Abstract Nanofluidic structures have over the last two decades emerged as a powerful platform for detailed analysis of DNA on the kilobase pair length scale. When DNA is confined to a nanochannel, the combination of excluded volume and DNA stiffness leads to the DNA being stretched to near its full contour length. Importantly, this stretching takes place at equilibrium, without any chemical modifications to the DNA. As a result, any DNA can be analyzed, such as DNA extracted from cells or circular DNA, and it is straight-forward to study reactions on the ends of linear DNA. In this comprehensive review, we first give a thorough description of the current understanding of the polymer physics of DNA and how that leads to stretching in nanochannels. We then describe how the versatility of nanofabrication can be used to design devices specifically tailored for the problem at hand, either by controlling the degree of confinement or enabling facile exchange of reagents to measure DNA–protein reaction kinetics. The remainder of the review focuses on two important applications of confining DNA in nanochannels. The first is optical DNA mapping, which provides the genomic sequence of intact DNA molecules in excess of 100 kilobase pairs in size, with kilobase pair resolution, through labeling strategies that are suitable for fluorescence microscopy. In this section, we highlight solutions to the technical aspects of genomic mapping, including the use of enzyme-based labeling and affinity-based labeling to produce the genomic maps, rather than recent applications in human genetics. The second is DNA–protein interactions, and several recent examples of such studies on DNA compaction, filamentous protein complexes, and reactions with DNA ends are presented. Taken together, these two applications demonstrate the power of DNA confinement and nanofluidics in genomics, molecular biology, and biophysics.more » « less
-
Abstract. Evidence has accumulated that secondary organic aerosols (SOAs) exhibit complex morphologies with multiple phases that can adopt amorphous semisolid or glassy phase states. However, experimental analysis and numerical modeling on the formation and evolution of SOA still often employ equilibrium partitioning with an ideal mixing assumption in the particle phase. Here we apply the kinetic multilayer model of gas–particle partitioning (KM-GAP) to simulate condensation of semi-volatile species into a core–shell phase-separated particle to evaluate equilibration timescales of SOA partitioning. By varying bulk diffusivity and the activity coefficient of the condensing species in the shell, we probe the complex interplay of mass transfer kinetics and the thermodynamics of partitioning. We found that the interplay of non-ideality and phase state can impact SOA partitioning kinetics significantly. The effect of non-ideality on SOA partitioning is slight for liquid particles but becomes prominent in semisolid or solid particles. If the condensing species is miscible with a low activity coefficient in the viscous shell phase, the particle can reach equilibrium with the gas phase long before the dissolution of concentration gradients in the particle bulk. For the condensation of immiscible species with a high activity coefficient in the semisolid shell, the mass concentration in the shell may become higher or overshoot its equilibrium concentration due to slow bulk diffusion through the viscous shell for excess mass to be transferred to the core phase. Equilibration timescales are shorter for the condensation of lower-volatility species into semisolid shell; as the volatility increases, re-evaporation becomes significant as desorption is faster for volatile species than bulk diffusion in a semisolid matrix, leading to an increase in equilibration timescale. We also show that the equilibration timescale is longer in an open system relative to a closed system especially for partitioning of miscible species; hence, caution should be exercised when interpreting and extrapolating closed-system chamber experimental results to atmosphere conditions. Our results provide a possible explanation for discrepancies between experimental observations of fast particle–particle mixing and predictions of long mixing timescales in viscous particles and provide useful insights into description and treatment of SOA in aerosol models.more » « less
-
more » « less
Abstract In this paper, we demonstrate a molecular system for the first
active self-assemblylinear DNA polymer that exhibits programmable molecular exponential growth in real time, also the first to implement “internal” parallel insertion that does not rely on adding successive layers to “external” edges for growth. Approaches like this can produce enhanced exponential growth behavior that is less limited by volume and external surface interference, for an early step toward efficiently building two and three dimensional shapes in logarithmic time. We experimentally demonstrate the division of these polymers via the addition of a single DNA complex that competes with the insertion mechanism and results in the exponential growth of a population of polymers per unit time. In the supplementary material, we note that an “extension” beyond conventional Turing machine theory is needed to theoretically analyze exponential growth itself in programmable physical systems. Sequential physical Turing Machines that run a roughly constant number of Turing steps per unit time cannot achieve an exponential growth of structure per time. In contrast, the “active” self-assembly model in this paper, computationally equivalent to a Push-Down Automaton, is exponentially fast when implemented in molecules, but istaxonomically less powerful than a Turing machine. In this sense, a physical Push-Down Automaton can bemore powerful than a sequential physical Turing Machine, even though the Turing Machine can compute any computable function. A need for an “extended” computational/physical theory arises, described in the supplementary material section S1.
