1
|
Niblo JK, Swartley JR, Zhang Z, DuBay KH. 2D capsid formation within an oscillatory energy landscape: orderly self-assembly depends on the interplay between a dynamic potential and intrinsic relaxation times. SOFT MATTER 2024. [PMID: 39046256 DOI: 10.1039/d4sm00455h] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 07/25/2024]
Abstract
Multiple dissipative self-assembly protocols designed to create novel structures or to reduce kinetic traps have recently emerged. Specifically, temporal oscillations of particle interactions have been shown effective at both aims, but investigations thus far have focused on systems of simple colloids or their binary mixtures. In this work, we expand our understanding of the effect of temporally oscillating interactions to a two-dimensional coarse-grained viral capsid-like model that undergoes a self-limited assembly. This model includes multiple intrinsic relaxation times due to the internal structure of the capsid subunits and, under certain interaction regimes, proceeds via a two-step nucleation mechanism. We find that oscillations much faster than the local intrinsic relaxation times can be described via a time averaged inter-particle potential across a wide range of interaction strengths, while oscillations much slower than these relaxation times result in structures that adapt to the attraction strength of the current half-cycle. Interestingly, oscillation periods similar to these relaxation times shift the interaction window over which orderly assembly occurs by enabling error correction during the half-cycles with weaker attractions. Our results provide fundamental insights to non-equilibrium self-assembly on temporally variant energy landscapes.
Collapse
Affiliation(s)
- Jessica K Niblo
- Department of Chemistry, University of Virginia, McCormick Road, PO Box 400319, Charlottesville, VA 22903-4319, USA.
| | - Jacob R Swartley
- Department of Chemistry, University of Virginia, McCormick Road, PO Box 400319, Charlottesville, VA 22903-4319, USA.
| | - Zhongmin Zhang
- Department of Chemistry, University of North Carolina at Chapel Hill, Campus Box 3290, Chapel Hill, NC 27599-3290, USA
| | - Kateri H DuBay
- Department of Chemistry, University of Virginia, McCormick Road, PO Box 400319, Charlottesville, VA 22903-4319, USA.
| |
Collapse
|
2
|
Osat S, Metson J, Kardar M, Golestanian R. Escaping Kinetic Traps Using Nonreciprocal Interactions. PHYSICAL REVIEW LETTERS 2024; 133:028301. [PMID: 39073937 DOI: 10.1103/physrevlett.133.028301] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/01/2023] [Accepted: 06/10/2024] [Indexed: 07/31/2024]
Abstract
Kinetic traps are a notorious problem in equilibrium statistical mechanics, where temperature quenches ultimately fail to bring the system to low energy configurations. Using multifarious self-assembly as a model system, we introduce a mechanism to escape kinetic traps by utilizing nonreciprocal interactions between components. Introducing nonequilibrium effects offered by broken action-reaction symmetry in the system pushes the trajectory of the system out of arrested dynamics. The dynamics of the model is studied using tools from the physics of interfaces and defects. Our proposal can find applications in self-assembly, glassy systems, and systems with arrested dynamics to facilitate escape from local minima in rough energy landscapes.
Collapse
|
3
|
Pimonova Y, Carpenter JE, Gruenwald M. Thermodynamic Stability Is a Poor Indicator of Cocrystallization in Models of Organic Molecules. J Am Chem Soc 2024; 146:2805-2815. [PMID: 38241026 DOI: 10.1021/jacs.3c13030] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2024]
Abstract
Cocrystallizing a given molecule with another can be useful for adjusting the physical properties of molecules in the solid state. However, most combinations of molecules do not readily cocrystallize but form either one-component crystals or amorphous solids. Computational methods of crystal structure prediction can, in principle, identify the thermodynamically stable cocrystal and thus predict if molecules will cocrystallize or not. However, the pronounced polymorphism and tendency of many organic molecules to form disordered solids suggest that kinetic factors can play an important role in cocrystallization. The question remains: if a binary system of molecules has a thermodynamically stable cocrystal, will it indeed cocrystallize? To address this question, we simulate the crystallization of more than 2600 distinct pairs of chiral model molecules of similar size in 2D and calculate accurate crystal energy landscapes for all of them. Our analysis shows that thermodynamic criteria alone are unreliable in the prediction of cocrystallization. While the vast majority of cocrystals that form in our simulations are thermodynamically favorable, most coformer systems that have a thermodynamically stable cocrystal do not cocrystallize. We furthermore show that cocrystallization rates increase 3-fold when coformers are used that do not form well-ordered single-component crystals. Our results suggest that kinetic factors of cocrystallization are decisive in many cases.
Collapse
Affiliation(s)
- Yulia Pimonova
- Department of Chemistry, University of Utah, Salt Lake City, Utah 84112, United States
| | - John E Carpenter
- Department of Chemistry, University of Utah, Salt Lake City, Utah 84112, United States
| | - Michael Gruenwald
- Department of Chemistry, University of Utah, Salt Lake City, Utah 84112, United States
| |
Collapse
|
4
|
Curatolo AI, Kimchi O, Goodrich CP, Krueger RK, Brenner MP. A computational toolbox for the assembly yield of complex and heterogeneous structures. Nat Commun 2023; 14:8328. [PMID: 38097568 PMCID: PMC10721878 DOI: 10.1038/s41467-023-43168-4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/31/2022] [Accepted: 11/02/2023] [Indexed: 12/17/2023] Open
Abstract
The self-assembly of complex structures from a set of non-identical building blocks is a hallmark of soft matter and biological systems, including protein complexes, colloidal clusters, and DNA-based assemblies. Predicting the dependence of the equilibrium assembly yield on the concentrations and interaction energies of building blocks is highly challenging, owing to the difficulty of computing the entropic contributions to the free energy of the many structures that compete with the ground state configuration. While these calculations yield well known results for spherically symmetric building blocks, they do not hold when the building blocks have internal rotational degrees of freedom. Here we present an approach for solving this problem that works with arbitrary building blocks, including proteins with known structure and complex colloidal building blocks. Our algorithm combines classical statistical mechanics with recently developed computational tools for automatic differentiation. Automatic differentiation allows efficient evaluation of equilibrium averages over configurations that would otherwise be intractable. We demonstrate the validity of our framework by comparison to molecular dynamics simulations of simple examples, and apply it to calculate the yield curves for known protein complexes and for the assembly of colloidal shells.
Collapse
Affiliation(s)
- Agnese I Curatolo
- School of Engineering and Applied Sciences, Harvard University, Cambridge, MA, 02138, USA
| | - Ofer Kimchi
- Lewis-Sigler Institute, Princeton University, Princeton, NJ, 08544, USA
| | - Carl P Goodrich
- Institute of Science and Technology Austria, A-3400, Klosterneuburg, Austria
| | - Ryan K Krueger
- School of Engineering and Applied Sciences, Harvard University, Cambridge, MA, 02138, USA
| | - Michael P Brenner
- School of Engineering and Applied Sciences, Harvard University, Cambridge, MA, 02138, USA.
- Department of Physics, Harvard University, Cambridge, MA, 02138, USA.
| |
Collapse
|
5
|
Das A, Limmer DT. Nonequilibrium design strategies for functional colloidal assemblies. Proc Natl Acad Sci U S A 2023; 120:e2217242120. [PMID: 37748070 PMCID: PMC10556551 DOI: 10.1073/pnas.2217242120] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/08/2022] [Accepted: 08/17/2023] [Indexed: 09/27/2023] Open
Abstract
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.
Collapse
Affiliation(s)
- Avishek Das
- Department of Chemistry, University of California, Berkeley, CA94720
| | - David T. Limmer
- Department of Chemistry, University of California, Berkeley, CA94720
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, CA94720
- Material Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, CA94720
- Kavli Energy NanoSciences Institute, University of California, Berkeley, CA94720
| |
Collapse
|
6
|
Lochab V, Ewim ED, Prakash S. Continuous flow microfluidics for colloidal particle assembly on porous substrates. SOFT MATTER 2023; 19:2564-2569. [PMID: 36942885 DOI: 10.1039/d2sm01414a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/18/2023]
Abstract
Self-assembly of colloidal particles for 'bottom-up' fabrication of various patterns and structures is critical for a range of applications including, but not limited to, energy generation and storage, nanomaterial structures, biomimetics, and biosensing. Multiple self-assembly techniques, such as substrate templating-via topological or chemical patterning-and solvent evaporation were discussed in our previous papers and have been developed for the deposition of patterned self-assembled structures, such as bands of colloidal particles, on various substrates. While the templating techniques are limited in applications due to the requirements for pattern-specific prior substrate engineering to fabricate the desired structure, solvent evaporation requires longer assembly times and precise control over environmental conditions. In this paper, a template-free, continuous flow process, which is facilitated by continuous solvent drainage through porous substrates, is demonstrated for the self-assembly of colloidal particles into high-aspect ratio (>103, length to width) structures, such as linear arrays or grid structures. Colloidal particles were assembled both on polymeric and metallic porous membranes, with rapid assembly times.
Collapse
Affiliation(s)
- Varun Lochab
- Department of Mechanical and Aerospace Engineering, The Ohio State University, Columbus, 43210, OH, USA.
| | - E Daniel Ewim
- Department of Mechanical and Aerospace Engineering, The Ohio State University, Columbus, 43210, OH, USA.
| | - Shaurya Prakash
- Department of Mechanical and Aerospace Engineering, The Ohio State University, Columbus, 43210, OH, USA.
| |
Collapse
|
7
|
Osat S, Golestanian R. Non-reciprocal multifarious self-organization. NATURE NANOTECHNOLOGY 2023; 18:79-85. [PMID: 36509920 PMCID: PMC9879770 DOI: 10.1038/s41565-022-01258-2] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/14/2022] [Accepted: 10/06/2022] [Indexed: 05/14/2023]
Abstract
A hallmark of living systems is the ability to employ a common set of building blocks that can self-organize into a multitude of different structures. This capability can only be afforded in non-equilibrium conditions, as evident from the energy-consuming nature of the plethora of such dynamical processes. To achieve automated dynamical control of such self-assembled structures and transitions between them, we need to identify the fundamental aspects of non-equilibrium dynamics that can enable such processes. Here we identify programmable non-reciprocal interactions as a tool to achieve such functionalities. The design rule is composed of reciprocal interactions that lead to the equilibrium assembly of the different structures, through a process denoted as multifarious self-assembly, and non-reciprocal interactions that give rise to non-equilibrium dynamical transitions between the structures. The design of such self-organized shape-shifting structures can be implemented at different scales, from nucleic acids and peptides to proteins and colloids.
Collapse
Affiliation(s)
- Saeed Osat
- Max Planck Institute for Dynamics and Self-Organization (MPIDS), Göttingen, Germany
| | - Ramin Golestanian
- Max Planck Institute for Dynamics and Self-Organization (MPIDS), Göttingen, Germany.
- Rudolf Peierls Centre for Theoretical Physics, University of Oxford, Oxford, UK.
| |
Collapse
|
8
|
Trubiano A, Hagan MF. Optimization of non-equilibrium self-assembly protocols using Markov state models. J Chem Phys 2022; 157:244901. [PMID: 36586982 PMCID: PMC9788858 DOI: 10.1063/5.0130407] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/11/2022] [Accepted: 12/04/2022] [Indexed: 12/12/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Anthony Trubiano
- Martin Fisher School of Physics, Brandeis University, Waltham, Massachusetts 02454, USA
| | - Michael F. Hagan
- Martin Fisher School of Physics, Brandeis University, Waltham, Massachusetts 02454, USA
| |
Collapse
|
9
|
Bone RA, Green JR. Optimizing dynamical functions for speed with stochastic paths. J Chem Phys 2022; 157:224101. [PMID: 36546817 DOI: 10.1063/5.0125479] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
Abstract
Living systems are built from microscopic components that function dynamically; they generate work with molecular motors, assemble and disassemble structures such as microtubules, keep time with circadian clocks, and catalyze the replication of DNA. How do we implement these functions in synthetic nanostructured materials to execute them before the onset of dissipative losses? Answering this question requires a quantitative understanding of when we can improve performance and speed while minimizing the dissipative losses associated with operating in a fluctuating environment. Here, we show that there are four modalities for optimizing dynamical functions that can guide the design of nanoscale systems. We analyze Markov models that span the design space: a clock, ratchet, replicator, and self-assembling system. Using stochastic thermodynamics and an exact expression for path probabilities, we classify these models of dynamical functions based on the correlation of speed with dissipation and with the chosen performance metric. We also analyze random networks to identify the model features that affect their classification and the optimization of their functionality. Overall, our results show that the possible nonequilibrium paths can determine our ability to optimize the performance of dynamical functions, despite ever-present dissipation, when there is a need for speed.
Collapse
Affiliation(s)
- Rebecca A Bone
- Department of Chemistry, University of Massachusetts Boston, Boston, Massachusetts 02125, USA
| | - Jason R Green
- Department of Chemistry, University of Massachusetts Boston, Boston, Massachusetts 02125, USA
| |
Collapse
|
10
|
Bupathy A, Frenkel D, Sastry S. Temperature protocols to guide selective self-assembly of competing structures. Proc Natl Acad Sci U S A 2022; 119:2119315119. [PMID: 35165184 PMCID: PMC8872760 DOI: 10.1073/pnas.2119315119] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 01/17/2022] [Indexed: 11/18/2022] Open
Abstract
Multicomponent self-assembly mixtures offer the possibility of encoding multiple target structures with the same set of interacting components. Selective retrieval of one of the stored structures has been attempted by preparing an initial state that favors the assembly of the required target, through seeding, concentration patterning, or specific choices of interaction strengths. This may not be possible in an experiment where on-the-fly reconfiguration of the building blocks to switch functionality may be required. In this paper, we explore principles of inverse design of a multicomponent, self-assembly mixture capable of encoding two competing structures that can be selected through simple temperature protocols. We design the target structures to realize the generic situation in which one of the targets has the lower nucleation barrier, while the other is globally more stable. We observe that, to avoid the formation of spurious or chimeric aggregates, the number of neighboring component pairs that occur in both structures should be minimal. Our design also requires the inclusion of components that are part of only one of the target structures. We observe, however, that to maximize the selectivity of retrieval, the component library itself should be maximally shared by the two targets, within such a constraint. We demonstrate that temperature protocols can be designed that lead to the formation of either one of the target structures with high selectivity. We discuss the important role played by secondary aggregation products in improving selectivity, which we term "vestigial aggregates."
Collapse
Affiliation(s)
- Arunkumar Bupathy
- Theoretical Sciences Unit, Jawaharlal Nehru Centre for Advanced Scientific Research, Bangalore 560064, India
| | - Daan Frenkel
- Centre for Computational Chemistry, Department of Chemistry, University of Cambridge, Cambridge CB2 1EW, United Kingdom
| | - Srikanth Sastry
- Theoretical Sciences Unit, Jawaharlal Nehru Centre for Advanced Scientific Research, Bangalore 560064, India;
| |
Collapse
|
11
|
Yuan M, Tang X, Han W. Anatomy and Formation Mechanisms of Early Amyloid-β Oligomers with Lateral Branching: Graph Network Analysis on Large-Scale Simulations. Chem Sci 2022; 13:2649-2660. [PMID: 35356670 PMCID: PMC8890322 DOI: 10.1039/d1sc06337e] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/15/2021] [Accepted: 02/08/2022] [Indexed: 11/29/2022] Open
Abstract
Oligomeric amyloid-β aggregates (AβOs) effectively trigger Alzheimer's disease-related toxicity, generating great interest in understanding their structures and formation mechanisms. However, AβOs are heterogeneous and transient, making their structure and formation difficult to study. Here, we performed graph network analysis of tens of microsecond massive simulations of early amyloid-β (Aβ) aggregations at near-atomic resolution to characterize AβO structures with sizes up to 20-mers. We found that AβOs exhibit highly curvilinear, irregular shapes with occasional lateral branches, consistent with recent cryo-electron tomography experiments. We also found that Aβ40 oligomers were more likely to develop branches than Aβ42 oligomers, explaining an experimental observation that only Aβ40 was trapped in network-like aggregates and exhibited slower fibrillization kinetics. Moreover, AβO architecture dissection revealed that their curvilinear appearance is related to the local packing geometries of neighboring peptides and that Aβ40's greater branching ability originates from specific C-terminal interactions at branching interfaces. Finally, we demonstrate that whether Aβ oligomerization causes oligomers to elongate or to branch depends on the sizes and shapes of colliding aggregates. Collectively, this study provides bottom-up structural information for understanding early Aβ aggregation and AβO toxicity. Graph network analysis on large-scale simulations uncovers the differential branching behaviours of large Aβ40 and Aβ42 oligomers.![]()
Collapse
Affiliation(s)
- Miao Yuan
- State Key Laboratory of Chemical Oncogenomics, Guangdong Provincial Key Laboratory of Chemical Genomics, School of Chemical Biology and Biotechnology, Peking University Shenzhen Graduate School Shenzhen 518055 China
| | - Xuan Tang
- State Key Laboratory of Chemical Oncogenomics, Guangdong Provincial Key Laboratory of Chemical Genomics, School of Chemical Biology and Biotechnology, Peking University Shenzhen Graduate School Shenzhen 518055 China
| | - Wei Han
- State Key Laboratory of Chemical Oncogenomics, Guangdong Provincial Key Laboratory of Chemical Genomics, School of Chemical Biology and Biotechnology, Peking University Shenzhen Graduate School Shenzhen 518055 China
| |
Collapse
|
12
|
Abstract
The inverse problem of designing component interactions to target emergent structure is fundamental to numerous applications in biotechnology, materials science, and statistical physics. Equally important is the inverse problem of designing emergent kinetics, but this has received considerably less attention. Using recent advances in automatic differentiation, we show how kinetic pathways can be precisely designed by directly differentiating through statistical physics models, namely free energy calculations and molecular dynamics simulations. We consider two systems that are crucial to our understanding of structural self-assembly: bulk crystallization and small nanoclusters. In each case, we are able to assemble precise dynamical features. Using gradient information, we manipulate interactions among constituent particles to tune the rate at which these systems yield specific structures of interest. Moreover, we use this approach to learn nontrivial features about the high-dimensional design space, allowing us to accurately predict when multiple kinetic features can be simultaneously and independently controlled. These results provide a concrete and generalizable foundation for studying nonstructural self-assembly, including kinetic properties as well as other complex emergent properties, in a vast array of systems.
Collapse
|
13
|
Trubiano A, Holmes-Cerfon M. Thermodynamic stability versus kinetic accessibility: Pareto fronts for programmable self-assembly. SOFT MATTER 2021; 17:6797-6807. [PMID: 34223604 DOI: 10.1039/d1sm00681a] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
A challenge in designing self-assembling building blocks is to ensure the target state is both thermodynamically stable and kinetically accessible. These two objectives are known to be typically in competition, but it is not known how to simultaneously optimize them. We consider this problem through the lens of multi-objective optimization theory: we develop a genetic algorithm to compute the Pareto fronts characterizing the tradeoff between equilibrium probability and folding rate, for a model system of small polymers of colloids with tunable short-ranged interaction energies. We use a coarse-grained model for the particles' dynamics that allows us to efficiently search over parameters, for systems small enough to be enumerated. For most target states there is a tradeoff when the number of types of particles is small, with medium-weak bonds favouring fast folding, and strong bonds favouring high equilibrium probability. The tradeoff disappears when the number of particle types reaches a value m*, that is usually much less than the total number of particles. This general approach of computing Pareto fronts allows one to identify the minimum number of design parameters to avoid a thermodynamic-kinetic tradeoff. However, we argue, by contrasting our coarse-grained model's predictions with those of Brownian dynamics simulations, that particles with short-ranged isotropic interactions should generically have a tradeoff, and avoiding it in larger systems will require orientation-dependent interactions.
Collapse
Affiliation(s)
- Anthony Trubiano
- Courant Institute of Mathematical Sciences, New York University, New York, New York 10012, USA.
| | - Miranda Holmes-Cerfon
- Courant Institute of Mathematical Sciences, New York University, New York, New York 10012, USA.
| |
Collapse
|
14
|
Hagan MF, Grason GM. Equilibrium mechanisms of self-limiting assembly. REVIEWS OF MODERN PHYSICS 2021; 93:025008. [PMID: 35221384 PMCID: PMC8880259 DOI: 10.1103/revmodphys.93.025008] [Citation(s) in RCA: 38] [Impact Index Per Article: 12.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/08/2023]
Abstract
Self-assembly is a ubiquitous process in synthetic and biological systems, broadly defined as the spontaneous organization of multiple subunits (e.g. macromolecules, particles) into ordered multi-unit structures. The vast majority of equilibrium assembly processes give rise to two states: one consisting of dispersed disassociated subunits, and the other, a bulk-condensed state of unlimited size. This review focuses on the more specialized class of self-limiting assembly, which describes equilibrium assembly processes resulting in finite-size structures. These systems pose a generic and basic question, how do thermodynamic processes involving non-covalent interactions between identical subunits "measure" and select the size of assembled structures? In this review, we begin with an introduction to the basic statistical mechanical framework for assembly thermodynamics, and use this to highlight the key physical ingredients that ensure equilibrium assembly will terminate at finite dimensions. Then, we introduce examples of self-limiting assembly systems, and classify them within this framework based on two broad categories: self-closing assemblies and open-boundary assemblies. These include well-known cases in biology and synthetic soft matter - micellization of amphiphiles and shell/tubule formation of tapered subunits - as well as less widely known classes of assemblies, such as short-range attractive/long-range repulsive systems and geometrically-frustrated assemblies. For each of these self-limiting mechanisms, we describe the physical mechanisms that select equilibrium assembly size, as well as potential limitations of finite-size selection. Finally, we discuss alternative mechanisms for finite-size assemblies, and draw contrasts with the size-control that these can achieve relative to self-limitation in equilibrium, single-species assemblies.
Collapse
Affiliation(s)
- Michael F Hagan
- Martin Fisher School of Physics, Brandeis University, Waltham, MA 02454, USA
| | - Gregory M Grason
- Department of Polymer Science and Engineering, University of Massachusetts, Amherst, MA 01003, USA
| |
Collapse
|
15
|
Giammona J, Campàs O. Physical constraints on early blastomere packings. PLoS Comput Biol 2021; 17:e1007994. [PMID: 33497383 PMCID: PMC7864451 DOI: 10.1371/journal.pcbi.1007994] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/25/2020] [Revised: 02/05/2021] [Accepted: 12/15/2020] [Indexed: 11/24/2022] Open
Abstract
At very early embryonic stages, when embryos are composed of just a few cells, establishing the correct packing arrangements (contacts) between cells is essential for the proper development of the organism. As early as the 4-cell stage, the observed cellular packings in different species are distinct and, in many cases, differ from the equilibrium packings expected for simple adherent and deformable particles. It is unclear what are the specific roles that different physical parameters, such as the forces between blastomeres, their division times, orientation of cell division and embryonic confinement, play in the control of these packing configurations. Here we simulate the non-equilibrium dynamics of cells in early embryos and systematically study how these different parameters affect embryonic packings at the 4-cell stage. In the absence of embryo confinement, we find that cellular packings are not robust, with multiple packing configurations simultaneously possible and very sensitive to parameter changes. Our results indicate that the geometry of the embryo confinement determines the packing configurations at the 4-cell stage, removing degeneracy in the possible packing configurations and overriding division rules in most cases. Overall, these results indicate that physical confinement of the embryo is essential to robustly specify proper cellular arrangements at very early developmental stages. At the initial stages of embryogenesis, the precise arrangement of cells in the embryo is critical to ensure that each cell gets the right chemical and physical signals to guide the formation of the organism. Even when the embryo is made of only four cells, different species feature varying cellular arrangements: cells in mouse embryos arrange as a tetrahedron, in the nematode worm C. elegans cells make a diamond and in sea urchins cells arrange in a square configuration. How do cells in embryos of different species control their arrangements? Using computer simulations, we studied how cell divisions, physical contacts between cells and the confinement of the embryo by an eggshell affect the arrangements of cells when the embryos have only 4 cells. We find that the shape of the confining eggshell plays a key role in controlling the cell arrangements, removing unwanted arrangements and robustly specifying the proper contacts between cells. Our results highlight the important roles of embryonic confinement in establishing the proper cell-cell contacts as the embryo starts to develop.
Collapse
Affiliation(s)
- James Giammona
- California NanoSystems Institute, University of California, Santa Barbara, California, United States of America
- Department of Physics, University of California, Santa Barbara, California, United States of America
| | - Otger Campàs
- California NanoSystems Institute, University of California, Santa Barbara, California, United States of America
- Department of Mechanical Engineering, University of California, Santa Barbara, California, United States of America
- Department of Molecular, Cell and Developmental Biology, University of California, Santa Barbara, California, United States of America
- Center for Bioengineering, University of California, Santa Barbara, California, United States of America
- Cluster of Excellence Physics of Life, TU Dresden, Dresden, Germany
- * E-mail:
| |
Collapse
|
16
|
Das A, Limmer DT. Variational design principles for nonequilibrium colloidal assembly. J Chem Phys 2021; 154:014107. [DOI: 10.1063/5.0038652] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
Affiliation(s)
- Avishek Das
- Department of Chemistry, University of California, Berkeley, California 94609, USA
| | - David T. Limmer
- Department of Chemistry, University of California, Berkeley, California 94609, USA
- Kavli Energy NanoScience Institute, Berkeley, California 94609, USA
- Materials Science Division, Lawrence Berkeley National Laboratory, Berkeley, California 94609, USA
- Chemical Science Division, Lawrence Berkeley National Laboratory, Berkeley, California 94609, USA
| |
Collapse
|
17
|
Holmes-Cerfon M. Simulating sticky particles: A Monte Carlo method to sample a stratification. J Chem Phys 2020; 153:164112. [PMID: 33138386 DOI: 10.1063/5.0019550] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023] Open
Abstract
Many problems in materials science and biology involve particles interacting with strong, short-ranged bonds that can break and form on experimental timescales. Treating such bonds as constraints can significantly speed up sampling their equilibrium distribution, and there are several methods to sample probability distributions subject to fixed constraints. We introduce a Monte Carlo method to handle the case when constraints can break and form. More generally, the method samples a probability distribution on a stratification: a collection of manifolds of different dimensions, where the lower-dimensional manifolds lie on the boundaries of the higher-dimensional manifolds. We show several applications of the method in polymer physics, self-assembly of colloids, and volume calculation in high dimensions.
Collapse
Affiliation(s)
- Miranda Holmes-Cerfon
- Courant Institute of Mathematical Sciences, New York University, New York, New York 10012, USA
| |
Collapse
|
18
|
Rogers WB. A mean-field model of linker-mediated colloidal interactions. J Chem Phys 2020; 153:124901. [DOI: 10.1063/5.0020578] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- W. Benjamin Rogers
- Martin A. Fisher School of Physics, Brandeis University, Waltham, Massachusetts, 02453, USA
| |
Collapse
|
19
|
Lowensohn J, Hensley A, Perlow-Zelman M, Rogers WB. Self-Assembly and Crystallization of DNA-Coated Colloids via Linker-Encoded Interactions. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2020; 36:7100-7108. [PMID: 32013444 DOI: 10.1021/acs.langmuir.9b03391] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
Coating colloidal particles with DNA is a promising strategy to make functional nanoscale materials because the particles can be programmed to spontaneously self-assemble into complex, ordered structures. In this Article, we explore the phase behavior and types of structures that can be formed when interactions between DNA-coated colloids are specified by linker DNA strands dispersed in solution. We show that linker-mediated interactions direct the self-assembly of colloids into equilibrium crystal structures. Furthermore, we demonstrate how different linker sequences and concentrations produce different crystal lattices, whose symmetry and compositional order are encoded exclusively by the linker-mediated interactions. These results illustrate how linkers can be used to separate the assembly instructions, encoded in the linkers, from the colloids themselves. We also examine the phase behavior of asymmetric linkers, which bind more strongly to one colloidal species than the other. We find that asymmetry strongly influences the concentration dependence of the colloidal interactions, which we explain using a mean-field model. We also find evidence that asymmetric linkers might help to reduce kinetic bottlenecks to colloidal crystallization. Together, our findings expand the design rules of linker-mediated self-assembly and make connections between the various schemes for programming assembly of DNA-coated colloids reported in the literature.
Collapse
Affiliation(s)
- Janna Lowensohn
- Martin A. Fisher School of Physics, Brandeis University, Waltham, Massachusetts 02453, United States
| | - Alexander Hensley
- Martin A. Fisher School of Physics, Brandeis University, Waltham, Massachusetts 02453, United States
| | - Michael Perlow-Zelman
- Martin A. Fisher School of Physics, Brandeis University, Waltham, Massachusetts 02453, United States
| | - W Benjamin Rogers
- Martin A. Fisher School of Physics, Brandeis University, Waltham, Massachusetts 02453, United States
| |
Collapse
|
20
|
Carpenter JE, Grünwald M. Heterogeneous Interactions Promote Crystallization and Spontaneous Resolution of Chiral Molecules. J Am Chem Soc 2020; 142:10755-10768. [DOI: 10.1021/jacs.0c02097] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
Affiliation(s)
- John E. Carpenter
- Department of Chemistry, University of Utah, Salt Lake City, Utah 84112, United States
| | - Michael Grünwald
- Department of Chemistry, University of Utah, Salt Lake City, Utah 84112, United States
| |
Collapse
|
21
|
Kumar R, Coli GM, Dijkstra M, Sastry S. Inverse design of charged colloidal particle interactions for self assembly into specified crystal structures. J Chem Phys 2019; 151:084109. [DOI: 10.1063/1.5111492] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Rajneesh Kumar
- Theoretical Sciences Unit, Jawaharlal Nehru Centre for Advanced Scientific Research, Bengaluru 560064, India
| | - Gabriele M. Coli
- Soft Condensed Matter, Debye Institute for Nanomaterials Science, Utrecht University, Princetonplein 5, 3584 CC Utrecht, The Netherlands
| | - Marjolein Dijkstra
- Soft Condensed Matter, Debye Institute for Nanomaterials Science, Utrecht University, Princetonplein 5, 3584 CC Utrecht, The Netherlands
| | - Srikanth Sastry
- Theoretical Sciences Unit, Jawaharlal Nehru Centre for Advanced Scientific Research, Bengaluru 560064, India
| |
Collapse
|
22
|
Nicholson SB, Bone RA, Green JR. Typical Stochastic Paths in the Transient Assembly of Fibrous Materials. J Phys Chem B 2019; 123:4792-4802. [PMID: 31063371 DOI: 10.1021/acs.jpcb.9b02811] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
Abstract
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.
Collapse
Affiliation(s)
- Schuyler B Nicholson
- Department of Chemistry , University of Massachusetts Boston , Boston , Massachusetts 02125 , United States
| | - Rebecca A Bone
- Department of Chemistry , University of Massachusetts Boston , Boston , Massachusetts 02125 , United States
| | - Jason R Green
- Department of Chemistry , University of Massachusetts Boston , Boston , Massachusetts 02125 , United States.,Department of Physics , University of Massachusetts Boston , Boston , Massachusetts 02125 , United States.,Center for Quantum and Nonequilibrium Systems , University of Massachusetts Boston , Boston , Massachusetts 02125 , United States
| |
Collapse
|
23
|
Abstract
Structural hierarchy, in which materials possess distinct features on multiple length scales, is ubiquitous in nature. Diverse biological materials, such as bone, cellulose, and muscle, have as many as 10 hierarchical levels. Structural hierarchy confers many mechanical advantages, including improved toughness and economy of material. However, it also presents a problem: Each hierarchical level adds a new source of assembly errors and substantially increases the information required for proper assembly. This seems to conflict with the prevalence of naturally occurring hierarchical structures, suggesting that a common mechanical source of hierarchical robustness may exist. However, our ability to identify such a unifying phenomenon is limited by the lack of a general mechanical framework for structures exhibiting organization on disparate length scales. Here, we use simulations to substantiate a generalized model for the tensile stiffness of hierarchical filamentous networks with a nested, dilute triangular lattice structure. Following seminal work by Maxwell and others on criteria for stiff frames, we extend the concept of connectivity in network mechanics and find a similar dependence of material stiffness upon each hierarchical level. Using this model, we find that stiffness becomes less sensitive to errors in assembly with additional levels of hierarchy; although surprising, we show that this result is analytically predictable from first principles and thus potentially model independent. More broadly, this work helps account for the success of hierarchical, filamentous materials in biology and materials design and offers a heuristic for ensuring that desired material properties are achieved within the required tolerance.
Collapse
|
24
|
Guan Q, Yuan L, Gu A, Liang G. Fabrication of In Situ Nanofiber-Reinforced Molecular Composites by Nonequilibrium Self-Assembly. ACS APPLIED MATERIALS & INTERFACES 2018; 10:39293-39306. [PMID: 30338977 DOI: 10.1021/acsami.8b15037] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
Although the concept of molecular composites (MCs) is very promising, there are major obstacles arising from the immiscibility of the rigid-rod with the random-coil polymers. Here, we developed a novel method for fabricating an in situ reinforced MC system with nonequilibrium self-assembled nanofibrous structures based on bisphenol A epoxy resin, 4,4'-diaminodiphenylsulfone, bismaleimide, and a polyphenylene ether (PPO) oligomer. A variety of spectroscopic and morphological techniques were used to probe the structural evolution from the emergence of nanofibrils, to growth and aggregation of nanofibers, and then to the formation of in situ reinforced MC with strong interfacial interactions. The in situ nanofibers within the polymer matrix could be formed by the polymerization force extruding the PPO phase through the interspaces within the simultaneous interpenetrating network polymers during the cure process of the thermosetting resin system. Compared to the control sample, the in situ nanofiber-reinforced MC exhibited better thermal properties and flame retardancy. In particular, the obtained MC showed a significant improvement in glass transition temperature and mechanical properties, which were mainly attributed to the restriction of high thermal stability of PPO on the segmental motion of polymer chains, the toughening and reinforcement behaviors of PPO nanofibers on the matrix, and the chemical interaction at the PPO/matrix interface.
Collapse
Affiliation(s)
- Qingbao Guan
- Jiangsu Key Laboratory of Advanced Functional Polymer Design and Application, Department of Materials Science and Engineering, College of Chemistry, Chemical Engineering and Materials Science , Soochow University , Suzhou 215123 , P. R. China
| | - Li Yuan
- Jiangsu Key Laboratory of Advanced Functional Polymer Design and Application, Department of Materials Science and Engineering, College of Chemistry, Chemical Engineering and Materials Science , Soochow University , Suzhou 215123 , P. R. China
| | - Aijuan Gu
- Jiangsu Key Laboratory of Advanced Functional Polymer Design and Application, Department of Materials Science and Engineering, College of Chemistry, Chemical Engineering and Materials Science , Soochow University , Suzhou 215123 , P. R. China
| | - Guozheng Liang
- Jiangsu Key Laboratory of Advanced Functional Polymer Design and Application, Department of Materials Science and Engineering, College of Chemistry, Chemical Engineering and Materials Science , Soochow University , Suzhou 215123 , P. R. China
| |
Collapse
|
25
|
Fonseca P, Romano F, Schreck JS, Ouldridge TE, Doye JPK, Louis AA. Multi-scale coarse-graining for the study of assembly pathways in DNA-brick self-assembly. J Chem Phys 2018; 148:134910. [PMID: 29626893 DOI: 10.1063/1.5019344] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023] Open
Abstract
Inspired by recent successes using single-stranded DNA tiles to produce complex structures, we develop a two-step coarse-graining approach that uses detailed thermodynamic calculations with oxDNA, a nucleotide-based model of DNA, to parametrize a coarser kinetic model that can reach the time and length scales needed to study the assembly mechanisms of these structures. We test the model by performing a detailed study of the assembly pathways for a two-dimensional target structure made up of 334 unique strands each of which are 42 nucleotides long. Without adjustable parameters, the model reproduces a critical temperature for the formation of the assembly that is close to the temperature at which assembly first occurs in experiments. Furthermore, the model allows us to investigate in detail the nucleation barriers and the distribution of critical nucleus shapes for the assembly of a single target structure. The assembly intermediates are compact and highly connected (although not maximally so), and classical nucleation theory provides a good fit to the height and shape of the nucleation barrier at temperatures close to where assembly first occurs.
Collapse
Affiliation(s)
- Pedro Fonseca
- Rudolf Peierls Centre for Theoretical Physics, University of Oxford, 1 Keble Road, Oxford OX1 3NP, United Kingdom
| | - Flavio Romano
- Dipartimento di Scienze Molecolari e Nanosistemi, Università Ca' Foscari, Via Torino 155, 30172 Venezia Mestre, Italy
| | - John S Schreck
- Department of Chemical Engineering, Columbia University, 500 W 120th St., New York, New York 10027, USA
| | - Thomas E Ouldridge
- Department of Bioengineering, Imperial College London, London SW7 2AZ, United Kingdom
| | - Jonathan P K Doye
- Physical and Theoretical Chemistry Laboratory, Department of Chemistry, University of Oxford, South Parks Road, Oxford OX1 3QZ, United Kingdom
| | - Ard A Louis
- Rudolf Peierls Centre for Theoretical Physics, University of Oxford, 1 Keble Road, Oxford OX1 3NP, United Kingdom
| |
Collapse
|
26
|
Abstract
A general formulation for constructing addressable atomic clusters is introduced, based on one or more reference structures. By modifying the well depths in a given interatomic potential in favour of nearest-neighbour interactions that are defined in the reference(s), the potential energy landscape can be biased to make a particular permutational isomer the global minimum. The magnitude of the bias changes the resulting potential energy landscape systematically, providing a framework to produce clusters that should self-organise efficiently into the target structure. These features are illustrated for small systems, where all the relevant local minima and transition states can be identified, and for the low-energy regions of the landscape for larger clusters. For a 55-particle cluster, it is possible to design a target structure from a transition state of the original potential and to retain this structure in a doubly addressable landscape. Disconnectivity graphs based on local minima that have no direct connections to a lower minimum provide a helpful way to visualise the larger databases. These minima correspond to the termini of monotonic sequences, which always proceed downhill in terms of potential energy, and we identify them as a class of biminimum. Multiple copies of the target cluster are treated by adding a repulsive term between particles with the same address to maintain distinguishable targets upon aggregation. By tuning the magnitude of this term, it is possible to create assemblies of the target cluster corresponding to a variety of structures, including rings and chains.
Collapse
Affiliation(s)
- David J Wales
- University Chemical Laboratories, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| |
Collapse
|
27
|
Abstract
Contemporary chemical and material engineering often takes inspiration from nature, targeting for example strong yet light materials and materials composed of highly ordered domains at multiple different lengthscales for fundamental science and applications in e.g. sensing, catalysis, coating technology, and delivery. The preparation of such hierarchically structured functional materials through guided bottom-up assembly of synthetic building blocks requires a high level of control over their synthesis, interactions and assembly pathways. In this perspective we showcase recent work demonstrating how molecular control can be exploited to direct colloidal assembly into responsive materials with mechanical, optical or electrical properties that can be adjusted post-synthesis with external cues.
Collapse
Affiliation(s)
- M Gerth
- Laboratory of Physical Chemistry, and Institute of Complex Molecular Systems, Eindhoven University of Technology, P.O. Box 513, 5600 MD, Eindhoven, The Netherlands
| | | |
Collapse
|
28
|
Fejer SN, Mantell RG, Wales DJ. Designing hierarchical molecular complexity: icosahedra of addressable icosahedra. Mol Phys 2018. [DOI: 10.1080/00268976.2018.1439190] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/17/2022]
Affiliation(s)
| | | | - David J. Wales
- University Chemical Laboratories, University of Cambridge , Cambridge, UK
| |
Collapse
|
29
|
Pestana LR, Minnetian N, Lammers LN, Head-Gordon T. Dynamical inversion of the energy landscape promotes non-equilibrium self-assembly of binary mixtures. Chem Sci 2018; 9:1640-1646. [PMID: 29675210 PMCID: PMC5887813 DOI: 10.1039/c7sc03524a] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/11/2017] [Accepted: 12/26/2017] [Indexed: 11/21/2022] Open
Abstract
When driven out of equilibrium, many diverse systems can form complex spatial and dynamical patterns, even in the absence of attractive interactions.
When driven out of equilibrium, many diverse systems can form complex spatial and dynamical patterns, even in the absence of attractive interactions. Using kinetic Monte Carlo simulations, we investigate the phase behavior of a binary system of particles of dissimilar size confined between semiflexible planar surfaces, in which the nanoconfinement introduces a non-local coupling between particles, which we model as an activation energy barrier to diffusion that decreases with the local fraction of the larger particle. The system autonomously reaches a cyclical non-equilibrium state characterized by the formation and dissolution of metastable micelle-like clusters with the small particles in the core and the large ones in the surrounding corona. The power spectrum of the fluctuations in the aggregation number exhibits 1/f noise reminiscent of self-organized critical systems. We suggest that the dynamical metastability of the micellar structures arises from an inversion of the energy landscape, in which the relaxation dynamics of one of the species induces a metastable phase for the other species.
Collapse
Affiliation(s)
- Luis Ruiz Pestana
- Chemical Sciences Division , Lawrence Berkeley National Laboratory , University of California , Berkeley , USA .
| | | | - Laura Nielsen Lammers
- Earth and Environmental Science Area , Lawrence Berkeley National Laboratory , University of California , Berkeley , USA.,Department of Environmental Science, Policy, and Management , University of California , Berkeley , USA
| | - Teresa Head-Gordon
- Chemical Sciences Division , Lawrence Berkeley National Laboratory , University of California , Berkeley , USA . .,Department of Chemistry , University of California , Berkeley , USA.,Department of Bioengineering , University of California , Berkeley , USA.,Department of Chemical and Biomolecular Engineering , University of California , Berkeley , USA
| |
Collapse
|
30
|
Halverson JD, Tkachenko AV. Communication: Programmable self-assembly of thin-shell mesostructures. J Chem Phys 2017; 147:141103. [DOI: 10.1063/1.4999654] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023] Open
Affiliation(s)
- Jonathan D. Halverson
- Center for Functional Nanomaterials, Brookhaven National Laboratory, Upton, New York 11973, USA
| | - Alexei V. Tkachenko
- Center for Functional Nanomaterials, Brookhaven National Laboratory, Upton, New York 11973, USA
| |
Collapse
|
31
|
Pietropaolo A, Tang S, Raymo FM. Free-energy predictions and absorption spectra calculations for supramolecular nanocarriers and their photoactive cargo. NANOSCALE 2017; 9:4989-4994. [PMID: 28383081 DOI: 10.1039/c7nr00839b] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
We reconstructed the free-energy landscape for supramolecular nanoparticles of amphiphilic methacrylated-based co-polymers. Their self-assembly in aqueous solution and encapsulation of borondipyrromethene (BODIPY) derivatives were enforced through atomistic free-energy simulations. The BODIPY binding modes detected in each of the free-energy basins were validated through a comparison of theoretical absorption spectra, calculated at the TD-DFT level, to their experimental counterparts. The nanoparticle distribution is controlled within a thermodynamic regime, with free-energy barriers approaching 8 kcal mol-1, enabling the existence of different-sized nanoparticles in aqueous solution at room temperature. Two types of supramolecular morphologies were identified. One is compact and spherical in shape and the other is large and donut-like, with the former more stable than the latter by 4 kcal mol-1. The morphology of the supramolecular host affects the binding mode of the BODIPY guests. Stacked BODIPY aggregates are encapsulated in the spherical nanocarriers, whereas isolated chromophores associate with the donut-shaped assemblies.
Collapse
Affiliation(s)
- Adriana Pietropaolo
- Dipartimento di Scienze della Salute, Università di Catanzaro, Viale Europa 88100, Catanzaro, Italy.
| | | | | |
Collapse
|
32
|
Abstract
Colloidal particles endowed with specific time-dependent interactions are a promising route for realizing artificial materials that have the properties of living ones. Previous work has demonstrated how this system can give rise to self-replication. Here, we introduce the process of colloidal catalysis, in which clusters of particles catalyze the creation of other clusters through templating reactions. Surprisingly, we find that simple templating rules generically lead to the production of huge numbers of clusters. The templating reactions among this sea of clusters give rise to an exponentially growing catalytic cycle, a specific realization of Dyson's notion of an exponentially growing metabolism. We demonstrate this behavior with a fixed set of interactions between particles chosen to allow a catalysis of a specific six-particle cluster from a specific seven-particle cluster, yet giving rise to the catalytic production of a sea of clusters of sizes between 2 and 11 particles. The fact that an exponentially growing cycle emerges naturally from such a simple scheme demonstrates that the emergence of exponentially growing metabolisms could be simpler than previously imagined.
Collapse
|
33
|
Kallus Y, Holmes-Cerfon M. Free energy of singular sticky-sphere clusters. Phys Rev E 2017; 95:022130. [PMID: 28297917 DOI: 10.1103/physreve.95.022130] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/16/2016] [Indexed: 06/06/2023]
Abstract
Networks of particles connected by springs model many condensed-matter systems, from colloids interacting with a short-range potential and complex fluids near jamming, to self-assembled lattices and various metamaterials. Under small thermal fluctuations the vibrational entropy of a ground state is given by the harmonic approximation if it has no zero-frequency vibrational modes, yet such singular modes are at the epicenter of many interesting behaviors in the systems above. We consider a system of N spherical particles, and directly account for the singularities that arise in the sticky limit where the pairwise interaction is strong and short ranged. Although the contribution to the partition function from singular clusters diverges in the limit, its asymptotic value can be calculated and depends on only two parameters, characterizing the depth and range of the potential. The result holds for systems that are second-order rigid, a geometric characterization that describes all known ground-state (rigid) sticky clusters. To illustrate the applications of our theory we address the question of emergence: how does crystalline order arise in large systems when it is strongly disfavored in small ones? We calculate the partition functions of all known rigid clusters up to N≤21 and show the cluster landscape is dominated by hyperstatic clusters (those with more than 3N-6 contacts); singular and isostatic clusters are far less frequent, despite their extra vibrational and configurational entropies. Since the most hyperstatic clusters are close to fragments of a close-packed lattice, this underlies the emergence of order in sticky-sphere systems, even those as small as N=10.
Collapse
Affiliation(s)
- Yoav Kallus
- Santa Fe Institute, 1399 Hyde Park Road, Santa Fe, New Mexico 87501, USA
| | - Miranda Holmes-Cerfon
- Courant Institute of Mathematical Sciences, New York University, New York, New York 10012, USA
| |
Collapse
|
34
|
Using active colloids as machines to weave and braid on the micrometer scale. Proc Natl Acad Sci U S A 2017; 114:257-262. [PMID: 28034922 DOI: 10.1073/pnas.1608838114] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Controlling motion at the microscopic scale is a fundamental goal in the development of biologically inspired systems. We show that the motion of active, self-propelled colloids can be sufficiently controlled for use as a tool to assemble complex structures such as braids and weaves out of microscopic filaments. Unlike typical self-assembly paradigms, these structures are held together by geometric constraints rather than adhesive bonds. The out-of-equilibrium assembly that we propose involves precisely controlling the 2D motion of active colloids so that their path has a nontrivial topology. We demonstrate with proof-of-principle Brownian dynamics simulations that, when the colloids are attached to long semiflexible filaments, this motion causes the filaments to braid. The ability of the active particles to provide sufficient force necessary to bend the filaments into a braid depends on a number of factors, including the self-propulsion mechanism, the properties of the filament, and the maximum curvature in the braid. Our work demonstrates that nonequilibrium assembly pathways can be designed using active particles.
Collapse
|
35
|
Abstract
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.
Collapse
|
36
|
Boles MA, Engel M, Talapin DV. Self-Assembly of Colloidal Nanocrystals: From Intricate Structures to Functional Materials. Chem Rev 2016; 116:11220-89. [PMID: 27552640 DOI: 10.1021/acs.chemrev.6b00196] [Citation(s) in RCA: 1049] [Impact Index Per Article: 131.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
Chemical methods developed over the past two decades enable preparation of colloidal nanocrystals with uniform size and shape. These Brownian objects readily order into superlattices. Recently, the range of accessible inorganic cores and tunable surface chemistries dramatically increased, expanding the set of nanocrystal arrangements experimentally attainable. In this review, we discuss efforts to create next-generation materials via bottom-up organization of nanocrystals with preprogrammed functionality and self-assembly instructions. This process is often driven by both interparticle interactions and the influence of the assembly environment. The introduction provides the reader with a practical overview of nanocrystal synthesis, self-assembly, and superlattice characterization. We then summarize the theory of nanocrystal interactions and examine fundamental principles governing nanocrystal self-assembly from hard and soft particle perspectives borrowed from the comparatively established fields of micrometer colloids and block copolymer assembly. We outline the extensive catalog of superlattices prepared to date using hydrocarbon-capped nanocrystals with spherical, polyhedral, rod, plate, and branched inorganic core shapes, as well as those obtained by mixing combinations thereof. We also provide an overview of structural defects in nanocrystal superlattices. We then explore the unique possibilities offered by leveraging nontraditional surface chemistries and assembly environments to control superlattice structure and produce nonbulk assemblies. We end with a discussion of the unique optical, magnetic, electronic, and catalytic properties of ordered nanocrystal superlattices, and the coming advances required to make use of this new class of solids.
Collapse
Affiliation(s)
- Michael A Boles
- Department of Chemistry and James Franck Institute, University of Chicago , Chicago, Illinois 60637, United States
| | - Michael Engel
- Institute for Multiscale Simulation, Friedrich-Alexander University Erlangen-Nürnberg , 91052 Erlangen, Germany.,Department of Chemical Engineering, University of Michigan , Ann Arbor, Michigan 48109, United States
| | - Dmitri V Talapin
- Department of Chemistry and James Franck Institute, University of Chicago , Chicago, Illinois 60637, United States.,Center for Nanoscale Materials, Argonne National Lab , Argonne, Illinois 60439, United States
| |
Collapse
|
37
|
Halverson JD, Tkachenko AV. Sequential programmable self-assembly: Role of cooperative interactions. J Chem Phys 2016; 144:094903. [DOI: 10.1063/1.4942615] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022] Open
Affiliation(s)
- Jonathan D. Halverson
- Center for Functional Nanomaterials, Brookhaven National Laboratory, Upton, New York 11973, USA
| | - Alexei V. Tkachenko
- Center for Functional Nanomaterials, Brookhaven National Laboratory, Upton, New York 11973, USA
| |
Collapse
|
38
|
Jacobs WM, Frenkel D. Self-Assembly of Structures with Addressable Complexity. J Am Chem Soc 2016; 138:2457-67. [DOI: 10.1021/jacs.5b11918] [Citation(s) in RCA: 66] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Affiliation(s)
- William M. Jacobs
- Department
of Chemistry and Chemical Biology, Harvard University, 12 Oxford
Street, Cambridge, Massachusetts 02138, United States
| | - Daan Frenkel
- Department
of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| |
Collapse
|
39
|
Miskin MZ, Khaira G, de Pablo JJ, Jaeger HM. Turning statistical physics models into materials design engines. Proc Natl Acad Sci U S A 2016; 113:34-9. [PMID: 26684770 PMCID: PMC4711868 DOI: 10.1073/pnas.1509316112] [Citation(s) in RCA: 46] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Despite the success statistical physics has enjoyed at predicting the properties of materials for given parameters, the inverse problem, identifying which material parameters produce given, desired properties, is only beginning to be addressed. Recently, several methods have emerged across disciplines that draw upon optimization and simulation to create computer programs that tailor material responses to specified behaviors. However, so far the methods developed either involve black-box techniques, in which the optimizer operates without explicit knowledge of the material's configuration space, or require carefully tuned algorithms with applicability limited to a narrow subclass of materials. Here we introduce a formalism that can generate optimizers automatically by extending statistical mechanics into the realm of design. The strength of this approach lies in its capability to transform statistical models that describe materials into optimizers to tailor them. By comparing against standard black-box optimization methods, we demonstrate how optimizers generated by this formalism can be faster and more effective, while remaining straightforward to implement. The scope of our approach includes possibilities for solving a variety of complex optimization and design problems concerning materials both in and out of equilibrium.
Collapse
Affiliation(s)
- Marc Z Miskin
- James Franck Institute and Department of Physics, The University of Chicago, Chicago, IL 60637;
| | - Gurdaman Khaira
- Institute for Molecular Engineering, The University of Chicago, Chicago, IL 60637
| | - Juan J de Pablo
- Institute for Molecular Engineering, The University of Chicago, Chicago, IL 60637; Institute for Molecular Engineering, Argonne National Laboratory, Lemont, IL 60439
| | - Heinrich M Jaeger
- James Franck Institute and Department of Physics, The University of Chicago, Chicago, IL 60637
| |
Collapse
|
40
|
Madge J, Miller MA. Design strategies for self-assembly of discrete targets. J Chem Phys 2015; 143:044905. [PMID: 26233162 DOI: 10.1063/1.4927671] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023] Open
Abstract
Both biological and artificial self-assembly processes can take place by a range of different schemes, from the successive addition of identical building blocks to hierarchical sequences of intermediates, all the way to the fully addressable limit in which each component is unique. In this paper, we introduce an idealized model of cubic particles with patterned faces that allows self-assembly strategies to be compared and tested. We consider a simple octameric target, starting with the minimal requirements for successful self-assembly and comparing the benefits and limitations of more sophisticated hierarchical and addressable schemes. Simulations are performed using a hybrid dynamical Monte Carlo protocol that allows self-assembling clusters to rearrange internally while still providing Stokes-Einstein-like diffusion of aggregates of different sizes. Our simulations explicitly capture the thermodynamic, dynamic, and steric challenges typically faced by self-assembly processes, including competition between multiple partially completed structures. Self-assembly pathways are extracted from the simulation trajectories by a fully extendable scheme for identifying structural fragments, which are then assembled into history diagrams for successfully completed target structures. For the simple target, a one-component assembly scheme is most efficient and robust overall, but hierarchical and addressable strategies can have an advantage under some conditions if high yield is a priority.
Collapse
Affiliation(s)
- Jim Madge
- Department of Chemistry, Durham University, South Road, Durham DH1 3LE, United Kingdom
| | - Mark A Miller
- Department of Chemistry, Durham University, South Road, Durham DH1 3LE, United Kingdom
| |
Collapse
|
41
|
Kaiser A, Popowa K, Löwen H. Active dipole clusters: From helical motion to fission. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2015; 92:012301. [PMID: 26274156 DOI: 10.1103/physreve.92.012301] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/11/2015] [Indexed: 06/04/2023]
Abstract
The structure of a finite particle cluster is typically determined by total energy minimization. Here we consider the case where a cluster of soft-sphere dipoles becomes active, i.e., when the individual particles exhibit an additional self-propulsion along their dipole moments. We numerically solve the overdamped equations of motion for soft-sphere dipoles in a solvent. Starting from an initial metastable dipolar cluster, the self-propulsion generates a complex cluster dynamics. The final cluster state has in general a structure widely different to the initial one, the details depend on the model parameters and on the protocol of how the self-propulsion is turned on. The center of mass of the cluster moves on a helical path, the details of which are governed by the initial cluster magnetization. An instantaneous switch to a high self-propulsion leads to fission of the cluster. However, fission does not occur if the self-propulsion is increased slowly to high strengths. Our predictions can be verified through experiments with self-phoretic colloidal Janus particles and for macroscopic self-propelled dipoles in a highly viscous solvent.
Collapse
Affiliation(s)
- Andreas Kaiser
- Institut für Theoretische Physik II: Weiche Materie, Heinrich-Heine-Universität Düsseldorf, Universitätsstraße 1, 40225 Düsseldorf, Germany
| | - Katarina Popowa
- Institut für Theoretische Physik II: Weiche Materie, Heinrich-Heine-Universität Düsseldorf, Universitätsstraße 1, 40225 Düsseldorf, Germany
| | - Hartmut Löwen
- Institut für Theoretische Physik II: Weiche Materie, Heinrich-Heine-Universität Düsseldorf, Universitätsstraße 1, 40225 Düsseldorf, Germany
| |
Collapse
|
42
|
Jacobs WM, Reinhardt A, Frenkel D. Communication: theoretical prediction of free-energy landscapes for complex self-assembly. J Chem Phys 2015; 142:021101. [PMID: 25591330 DOI: 10.1063/1.4905670] [Citation(s) in RCA: 31] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present a technique for calculating free-energy profiles for the nucleation of multicomponent structures that contain as many species as building blocks. We find that a key factor is the topology of the graph describing the connectivity of the target assembly. By considering the designed interactions separately from weaker, incidental interactions, our approach yields predictions for the equilibrium yield and nucleation barriers. These predictions are in good agreement with corresponding Monte Carlo simulations. We show that a few fundamental properties of the connectivity graph determine the most prominent features of the assembly thermodynamics. Surprisingly, we find that polydispersity in the strengths of the designed interactions stabilizes intermediate structures and can be used to sculpt the free-energy landscape for self-assembly. Finally, we demonstrate that weak incidental interactions can preclude assembly at equilibrium due to the combinatorial possibilities for incorrect association.
Collapse
Affiliation(s)
- William M Jacobs
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Aleks Reinhardt
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Daan Frenkel
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| |
Collapse
|
43
|
Rational design of self-assembly pathways for complex multicomponent structures. Proc Natl Acad Sci U S A 2015; 112:6313-8. [PMID: 25941388 DOI: 10.1073/pnas.1502210112] [Citation(s) in RCA: 74] [Impact Index Per Article: 8.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
The field of complex self-assembly is moving toward the design of multiparticle structures consisting of thousands of distinct building blocks. To exploit the potential benefits of structures with such "addressable complexity," we need to understand the factors that optimize the yield and the kinetics of self-assembly. Here we use a simple theoretical method to explain the key features responsible for the unexpected success of DNA-brick experiments, which are currently the only demonstration of reliable self-assembly with such a large number of components. Simulations confirm that our theory accurately predicts the narrow temperature window in which error-free assembly can occur. Even more strikingly, our theory predicts that correct assembly of the complete structure may require a time-dependent experimental protocol. Furthermore, we predict that low coordination numbers result in nonclassical nucleation behavior, which we find to be essential for achieving optimal nucleation kinetics under mild growth conditions. We also show that, rather surprisingly, the use of heterogeneous bond energies improves the nucleation kinetics and in fact appears to be necessary for assembling certain intricate 3D structures. This observation makes it possible to sculpt nucleation pathways by tuning the distribution of interaction strengths. These insights not only suggest how to improve the design of structures based on DNA bricks, but also point the way toward the creation of a much wider class of chemical or colloidal structures with addressable complexity.
Collapse
|
44
|
Walker DM, Tordesillas A, Brodu N, Dijksman JA, Behringer RP, Froyland G. Self-assembly in a near-frictionless granular material: conformational structures and transitions in uniaxial cyclic compression of hydrogel spheres. SOFT MATTER 2015; 11:2157-2173. [PMID: 25634109 DOI: 10.1039/c4sm02384f] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/04/2023]
Abstract
We use a Markov transition matrix-based analysis to explore the structures and structural transitions in a three-dimensional assembly of hydrogel spheres under cyclic uniaxial compression. We apply these methods on experimental data obtained from a packing of nearly frictionless hydrogel balls. This allows an exploration of the emergence and evolution of mesoscale internal structures - a key micromechanical property that governs self-assembly and self-organization in dense granular media. To probe the mesoscopic force network structure, we consider two structural state spaces: (i) a particle and its contacting neighbours, and (ii) a particle's local minimal cycle topology summarized by a cycle vector. In both spaces, our analysis of the transition dynamics reveals which structures and which sets of structures are most prevalent and most likely to transform into each other during the compression/decompression of the material. In compressed states, structures rich in 3-cycle or triangle topologies form in abundance. In contrast, in uncompressed states, transitions comprising poorly connected structures are dominant. An almost-invariant transition set within the cycle vector space is discovered that identifies an intermediate set of structures crucial to the material's transition from weakly jammed to strongly jammed, and vice versa. Preferred transition pathways are also highlighted and discussed with respect to thermo-micro-mechanical constitutive formulations.
Collapse
Affiliation(s)
- David M Walker
- Department of Mathematics and Statistics, University of Melbourne, Parkville, VIC 3010, Australia.
| | | | | | | | | | | |
Collapse
|
45
|
Undesired usage and the robust self-assembly of heterogeneous structures. Nat Commun 2015; 6:6203. [DOI: 10.1038/ncomms7203] [Citation(s) in RCA: 41] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/07/2014] [Accepted: 01/05/2015] [Indexed: 01/04/2023] Open
|
46
|
Zenk J, Schulman R. An assembly funnel makes biomolecular complex assembly efficient. PLoS One 2014; 9:e111233. [PMID: 25360818 PMCID: PMC4215988 DOI: 10.1371/journal.pone.0111233] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/28/2014] [Accepted: 09/30/2014] [Indexed: 11/18/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- John Zenk
- Chemical and Biomolecular Engineering, Johns Hopkins University, Baltimore, Maryland, United States of America
| | - Rebecca Schulman
- Chemical and Biomolecular Engineering, Johns Hopkins University, Baltimore, Maryland, United States of America
- Computer Science, Johns Hopkins University, Baltimore, Maryland, United States of America
- * E-mail:
| |
Collapse
|
47
|
Size limits of self-assembled colloidal structures made using specific interactions. Proc Natl Acad Sci U S A 2014; 111:15918-23. [PMID: 25349380 DOI: 10.1073/pnas.1411765111] [Citation(s) in RCA: 58] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022] Open
Abstract
We establish size limitations for assembling structures of controlled size and shape out of colloidal particles with short-ranged interactions. Through simulations we show that structures with highly variable shapes made out of dozens of particles can form with high yield, as long as each particle in the structure binds only to the particles in their local environment. To understand this, we identify the excited states that compete with the ground-state structure and demonstrate that these excited states have a completely topological characterization, valid when the interparticle interactions are short-ranged. This allows complete enumeration of the energy landscape and gives bounds on how large a colloidal structure can assemble with high yield. For large structures the yield can be significant, even with hundreds of particles.
Collapse
|
48
|
Grünwald M, Geissler PL. Patterns without patches: hierarchical self-assembly of complex structures from simple building blocks. ACS NANO 2014; 8:5891-7. [PMID: 24816138 PMCID: PMC4071109 DOI: 10.1021/nn500978p] [Citation(s) in RCA: 35] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/18/2014] [Accepted: 05/12/2014] [Indexed: 05/26/2023]
Abstract
Nanoparticles with "sticky patches" have long been proposed as building blocks for the self-assembly of complex structures. The synthetic realizability of such patchy particles, however, greatly lags behind predictions of patterns they could form. Using computer simulations, we show that structures of the same genre can be obtained from a solution of simple isotropic spheres, with control only over their sizes and a small number of binding affinities. In a first step, finite clusters of well-defined structure and composition emerge from natural dynamics with high yield. In effect a kind of patchy particle, these clusters can further assemble into a variety of complex superstructures, including filamentous networks, ordered sheets, and highly porous crystals.
Collapse
Affiliation(s)
- Michael Grünwald
- Computational Physics, University of Vienna, Sensengasse 8, 1090 Vienna, Austria
| | - Phillip L. Geissler
- Department of Chemistry, University of California, Berkeley, California 94720, United States
| |
Collapse
|
49
|
Long AW, Ferguson AL. Nonlinear Machine Learning of Patchy Colloid Self-Assembly Pathways and Mechanisms. J Phys Chem B 2014; 118:4228-44. [DOI: 10.1021/jp500350b] [Citation(s) in RCA: 59] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
Affiliation(s)
- Andrew W. Long
- Department of Materials Science
and Engineering, University of Illinois at Urbana−Champaign, Urbana, Illinois 61801, United States
| | - Andrew L. Ferguson
- Department of Materials Science
and Engineering, University of Illinois at Urbana−Champaign, Urbana, Illinois 61801, United States
| |
Collapse
|
50
|
Ruiz L, Keten S. Directing the self-assembly of supra-biomolecular nanotubes using entropic forces. SOFT MATTER 2014; 10:851-861. [PMID: 24652037 DOI: 10.1039/c3sm52600c] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/03/2023]
Abstract
Peptide self-assembly, ubiquitous in biology, is one of the most promising 'bottom-up' approaches for the generation of synthetic supramolecular architectures. However, directing the self-assembly of functional peptides into predictable ordered structures most often requires precise tuning of weak intermolecular forces. Existing strategies are generally based on specific interactions between molecular mediators that require complex chemical synthesis pathways and elaborated design rules. Here we establish a theoretical framework that delineates a generic route towards directing the self-assembly of small peptides by simply using entropic forces generated by the polymer chains attached to the peptides. We demonstrate the viability of this concept for polymer-conjugated peptide nanotubes using coarse-grained molecular dynamics (CGMD) simulations combined with theoretical calculations. We show that conjugated polymer chains create an entropic penalty due to chain confinement upon assembly, and illustrate that the self-assembly process can be directed by merely varying the degree of polymer conjugation. Specifically, the entropic penalty, and consequently, the binding energy between peptides can be greatly varied by changing the length and the number of conjugated polymers. Extending this concept for peptides with different degrees of conjugation reveals a path towards controlling the stacking sequence of binary mixtures. Remarkably, we find that a large disparity in the conjugation degree of the two peptides results in a preference towards alternating mixed sequences that minimize the entropic penalty of confinement in the thermodynamic limit. Our study explains recent experiments on polymer-peptide conjugates and sets the stage for utilizing entropic forces to guide the stacking sequence of functional macrocycles in tubular assemblies.
Collapse
Affiliation(s)
- Luis Ruiz
- Department of Civil & Environmental Engineering and Mechanical Engineering, Northwestern University, Evanston, Illinois 60208-3111, USA.
| | | |
Collapse
|