Fluctuation-dissipation ratios in the dynamics of self-assembly

Limits of economy and fidelity for programmable assembly of size-controlled triply periodic polyhedra

We propose and investigate an extension of the Caspar-Klug symmetry principles for viral capsid assembly to the programmable assembly of size-controlled triply periodic polyhedra, discrete variants of the Primitive, Diamond, and Gyroid cubic minimal surfaces. Inspired by a recent class of programmable DNA origami colloids, we demonstrate that the economy of design in these crystalline assemblies-in terms of the growth of the number of distinct particle species required with the increased size-scale (e.g., periodicity)-is comparable to viral shells. We further test the role of geometric specificity in these assemblies via dynamical assembly simulations, which show that conditions for simultaneously efficient and high-fidelity assembly require an intermediate degree of flexibility of local angles and lengths in programmed assembly. Off-target misassembly occurs via incorporation of a variant of disclination defects, generalized to the case of hyperbolic crystals. The possibility of these topological defects is a direct consequence of the very same symmetry principles that underlie the economical design, exposing a basic tradeoff between design economy and fidelity of programmable, size controlled assembly.

Effect of ionic strength on the assembly of simian vacuolating virus capsid protein around poly(styrene sulfonate)

Virus-like particles (VLPs) are noninfectious nanocapsules that can be used for drug delivery or vaccine applications. VLPs can be assembled from virus capsid proteins around a condensing agent, such as RNA, DNA, or a charged polymer. Electrostatic interactions play an important role in the assembly reaction. VLPs assemble from many copies of capsid protein, with a combinatorial number of intermediates. Hence, the mechanism of the reaction is poorly understood. In this paper, we combined solution small-angle X-ray scattering (SAXS), cryo-transmission electron microscopy (TEM), and computational modeling to determine the effect of ionic strength on the assembly of Simian Vacuolating Virus 40 (SV40)-like particles. We mixed poly(styrene sulfonate) with SV40 capsid protein pentamers at different ionic strengths. We then characterized the assembly product by SAXS and cryo-TEM. To analyze the data, we performed Langevin dynamics simulations using a coarse-grained model that revealed incomplete, asymmetric VLP structures consistent with the experimental data. We found that close to physiological ionic strength, [Formula: see text] VLPs coexisted with VP1 pentamers. At lower or higher ionic strengths, incomplete particles coexisted with pentamers and [Formula: see text] particles. Including the simulated structures was essential to explain the SAXS data in a manner that is consistent with the cryo-TEM images.

Nonequilibrium design strategies for functional colloidal assemblies

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

Self-assembly coupled to liquid-liquid phase separation

Liquid condensate droplets with distinct compositions of proteins and nucleic acids are widespread in biological cells. While it is known that such droplets, or compartments, can regulate irreversible protein aggregation, their effect on reversible self-assembly remains largely unexplored. In this article, we use kinetic theory and solution thermodynamics to investigate the effect of liquid-liquid phase separation on the reversible self-assembly of structures with well-defined sizes and architectures. We find that, when assembling subunits preferentially partition into liquid compartments, robustness against kinetic traps and maximum achievable assembly rates can be significantly increased. In particular, both the range of solution conditions leading to productive assembly and the corresponding assembly rates can increase by orders of magnitude. We analyze the rate equation predictions using simple scaling estimates to identify effects of liquid-liquid phase separation as a function of relevant control parameters. These results may elucidate self-assembly processes that underlie normal cellular functions or pathogenesis, and suggest strategies for designing efficient bottom-up assembly for nanomaterials applications.

Applications of Charge Detection Mass Spectrometry in Molecular Biology and Biotechnology

Charge detection mass spectrometry (CDMS) is a single-particle technique where the masses of individual ions are determined from simultaneous measurement of their mass-to-charge ratio (m/z) and charge. Masses are determined for thousands of individual ions, and then the results are binned to give a mass spectrum. Using this approach, accurate mass distributions can be measured for heterogeneous and high-molecular-weight samples that are usually not amenable to analysis by conventional mass spectrometry. Recent applications include heavily glycosylated proteins, protein complexes, protein aggregates such as amyloid fibers, infectious viruses, gene therapies, vaccines, and vesicles such as exosomes.

Equilibrium mechanisms of self-limiting assembly

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.

Mechanisms of Scaffold-Mediated Microcompartment Assembly and Size Control

This article describes a theoretical and computational study of the dynamical assembly of a protein shell around a complex consisting of many cargo molecules and long, flexible scaffold molecules. Our study is motivated by bacterial microcompartments, which are proteinaceous organelles that assemble around a condensed droplet of enzymes and reactants. As in many examples of cytoplasmic liquid-liquid phase separation, condensation of the microcompartment interior cargo is driven by flexible scaffold proteins that have weak multivalent interactions with the cargo. Our results predict that the shell size, amount of encapsulated cargo, and assembly pathways depend sensitively on properties of the scaffold, including its length and valency of scaffold-cargo interactions. Moreover, the ability of self-assembling protein shells to change their size to accommodate scaffold molecules of different lengths depends crucially on whether the spontaneous curvature radius of the protein shell is smaller or larger than a characteristic elastic length scale of the shell. Beyond natural microcompartments, these results have important implications for synthetic biology efforts to target alternative molecules for encapsulation by microcompartments or viral shells. More broadly, the results elucidate how cells exploit coupling between self-assembly and liquid-liquid phase separation to organize their interiors.

Effects of non-pairwise repulsion on nanoparticle assembly

Electrostatic interactions provide a convenient way to modulate interactions between nanoparticles, colloids, and biomolecules because they can be adjusted by the solution pH or salt concentration. While the presence of salt provides an easy method to control the net interparticle interaction, the nonlinearities arising from electrostatic screening make it difficult to quantify the strength of the interaction. In particular, when charged particles assemble into clusters or aggregates, nonlinear effects render the interactions strongly non-pairwise. Here, we report Brownian dynamics simulations to investigate the effect that the non-pairwise nature of electrostatic interactions has on nanoparticle assembly. We compare these simulations to a system in which the electrostatics are modeled by a strictly pairwise Yukawa potential. We find that both systems show a narrow range in parameter space where the particles form well-ordered crystals. Bordering this range are regions where the net interactions are too weak to stabilize aggregated structures or strong enough that the system becomes kinetically trapped in a gel. The non-pairwise potential differs from the pairwise system in the appearance of an amorphous state for strongly charged particles. This state appears because the many-body electrostatic interactions limit the maximum density achievable in an assembly.

Experimental Evidence for a Structural-Dynamical Transition in Trajectory Space

Among the key insights into the glass transition has been the identification of a nonequilibrium phase transition in trajectory space which reveals phase coexistence between the normal supercooled liquid (active phase) and a glassy state (inactive phase). Here, we present evidence that such a transition occurs in experiments. In colloidal hard spheres, we find a non-Gaussian distribution of trajectories leaning towards those rich in locally favored structures (LFSs), associated with the emergence of slow dynamics. This we interpret as evidence for a nonequilibrium transition to an inactive LFS-rich phase. Reweighting trajectories reveals a first-order phase transition in trajectory space between a normal liquid and a LFS-rich phase. We also find evidence for a purely dynamical transition in trajectory space.

Kinetics of Vapor-Solid Phase Transitions: Structure, Growth, and Mechanism

The kinetics of the separation between low and high density phases in a single component Lennard-Jones model is studied via molecular dynamics simulations, at very low temperatures, in the space dimension d=2. For densities close to the vapor branch of the coexistence curve, disconnected nanoscale clusters of the high density phase exhibit essentially ballistic motion. Starting from nearly circular shapes, at the time of nucleation, these clusters grow via sticky collisions, gaining filamentlike nonequilibrium structure at a later time, with a very low fractal dimensionality. The origin of the latter is shown to lie in the low mobility of the constituent particles, in the corresponding cluster reference frame, due to the (quasi-long-range) crystalline order. Standard self-similarity in the domain pattern, typically observed in the kinetics of phase transitions, is found to be absent. This invalidates the common method, that provides a growth law comparable to that in solid mixtures, of quantifying growth. An appropriate alternative approach, involving the fractality, quantifies the growth of the characteristic "length" to be a power law with time, the exponent being strongly temperature dependent. The observed growth law is in agreement with the outcome of a nonequilibrium kinetic theory.

The inherent dynamics of isotropic- and nematic-phase liquid crystals

The geodesic (shortest) pathways through the potential energy landscape of a liquid can be thought of as defining what its dynamics would be if thermal noise were removed, revealing what we have called the "inherent dynamics" of the liquid. We show how these inherent paths can be located for a model liquid crystal former, showing, in the process, how the molecular mechanisms of translation and reorientation compare in the isotropic and nematic phases of these systems. These mechanisms turn out to favor the preservation of local orientational order even under macroscopically isotropic conditions (a finding consistent with the experimental observation of pseudonematic domains in these cases), but disfavor the maintenance of macroscopic orientational order, even in the nematic phase. While the most efficient nematic pathways that maintain nematic order are indeed shorter than those that do not, it is apparently difficult for the system to locate these paths, suggesting that molecular motion in liquid-crystal formers is dynamically frustrated, and reinforcing the sense that there are strong analogies between liquid crystals and supercooled liquids.

Hierarchical assembly may be a way to make large information-rich structures

Self-assembly in the laboratory can now yield 'information-rich' nanostructures in which each component is of a distinct type and has a defined spatial position. Ensuring the thermodynamic stability of such structures requires inter-component interaction energies to increase logarithmically with structure size, in order to counter the entropy gained upon mixing component types in solution. However, self-assembly in the presence of strong interactions results in general in kinetic trapping, so suggesting a limit to the size of an (equilibrium) structure that can be self-assembled from distinguishable components. Here we study numerically a two-dimensional hierarchical assembly scheme already considered in experiment. We show that this scheme is immune to the kinetic traps associated with strong 'native' interactions (interactions designed to stabilize the intended structure), and so, in principle, offers a way to make large information-rich structures. In this scheme the size of an assembled structure scales exponentially with the stage of assembly, and assembly can continue as long as random motion is able to bring structures into contact. The resulting superstructure could provide a template for building in the third dimension. The chief drawback of this scheme is that it is particularly susceptible to kinetic traps that result from 'non-native' interactions (interactions not required to stabilize the intended structure); the scale on which such a scheme can be realized therefore depends upon how effectively this latter kind of interaction can be suppressed.

Self-assembly and crystallisation of indented colloids at a planar wall

We report experimental and simulation studies of the structure of a monolayer of indented ("lock and key") colloids, on a planar surface. On adding a non-absorbing polymer with prescribed radius and volume fraction, depletion interactions are induced between the colloids, with controlled range and strength. For spherical particles, this leads to crystallisation, but the indented colloids crystallise less easily than spheres, in both simulation and experiment. Nevertheless, simulations show that indented colloids do form plastic (rotator) crystals. We discuss the conditions under which this occurs, and the possibilities of lower-symmetry crystal states. We also comment on the kinetic accessibility of these states.

Scaling up self-assembly: bottom-up approaches to macroscopic particle organization

This review presents an overview of recent work in the field of non-Brownian particle self-assembly. Compared to nanoparticles that naturally self-assemble due to Brownian motion, larger, non-Brownian particles (d > 6 μm) are less prone to autonomously organize into crystalline arrays. The tendency for particle systems to experience immobilization and kinetic arrest grows with particle radius. In order to overcome this kinetic limitation, some type of external driver must be applied to act as an artificial "thermalizing force" upon non-Brownian particles, inducing particle motion and subsequent crystallization. Many groups have explored the use of various agitation methods to overcome the natural barriers preventing self-assembly to which non-Brownian particles are susceptible. The ability to create materials from a bottom-up approach with these characteristics would allow for precise control over their pore structure (size and distribution) and surface properties (topography, functionalization and area), resulting in improved regulation of key characteristics such as mechanical strength, diffusive properties, and possibly even photonic properties. This review will highlight these approaches, as well as discuss the potential impact of bottom-up macroscale particle assembly. The applications of such technology range from customizable and autonomously self-assembled niche microenvironments for drug delivery and tissue engineering to new acoustic dampening, battery, and filtration materials, among others. Additionally, crystals made from non-Brownian particles resemble naturally derived materials such as opals, zeolites, and biological tissue (i.e. bone, cartilage and lung), due to their high surface area, pore distribution, and tunable (multilevel) hierarchy.

The Role of Packaging Sites in Efficient and Specific Virus Assembly

During the life cycle of many single-stranded RNA viruses, including many human pathogens, a protein shell called the capsid spontaneously assembles around the viral genome. Understanding the mechanisms by which capsid proteins selectively assemble around the viral RNA amidst diverse host RNAs is a key question in virology. In one proposed mechanism, short sequences (packaging sites) within the genomic RNA promote rapid and efficient assembly through specific interactions with the capsid proteins. In this work, we develop a coarse-grained particle-based computational model for capsid proteins and RNA that represents protein-RNA interactions arising both from nonspecific electrostatics and from specific packaging site interactions. Using Brownian dynamics simulations, we explore how the efficiency and specificity of assembly depend on solution conditions (which control protein-protein and nonspecific protein-RNA interactions) and the strength and number of packaging sites. We identify distinct regions in parameter space in which packaging sites lead to highly specific assembly via different mechanisms and others in which packaging sites lead to kinetic traps. We relate these computational predictions to in vitro assays for specificity in which cognate viral RNAs compete against non-cognate RNAs for assembly by capsid proteins.

Modeling Viral Capsid Assembly

I present a review of the theoretical and computational methodologies that have been used to model the assembly of viral capsids. I discuss the capabilities and limitations of approaches ranging from equilibrium continuum theories to molecular dynamics simulations, and I give an overview of some of the important conclusions about virus assembly that have resulted from these modeling efforts. Topics include the assembly of empty viral shells, assembly around single-stranded nucleic acids to form viral particles, and assembly around synthetic polymers or charged nanoparticles for nanotechnology or biomedical applications. I present some examples in which modeling efforts have promoted experimental breakthroughs, as well as directions in which the connection between modeling and experiment can be strengthened.

Biomolecular engineering of virus-like particles aided by computational chemistry methods

Multi-scale investigation of VLP self-assembly aided by computational methods is facilitating the design, redesign, and modification of functionalized VLPs.

Patterns without patches: hierarchical self-assembly of complex structures from simple building blocks

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.

Detection of late intermediates in virus capsid assembly by charge detection mass spectrometry

The assembly of hundreds of identical proteins into an icosahedral virus capsid is a remarkable feat of molecular engineering. How this occurs is poorly understood. Key intermediates have been anticipated at the end of the assembly reaction, but it has not been possible to detect them. In this work we have used charge detection mass spectrometry to identify trapped intermediates from late in the assembly of the hepatitis B virus T = 4 capsid, a complex of 120 protein dimers. Prominent intermediates are found with 104/105, 110/111, and 117/118 dimers. Cryo-EM observations indicate the intermediates are incomplete capsids and, hence, on the assembly pathway. On the basis of their stability and kinetic accessibility we have proposed plausible structures. The prominent trapped intermediate with 104 dimers is attributed to an icosahedron missing two neighboring facets, the 111-dimer species is assigned to an icosahedron missing a single facet, and the intermediate with 117 dimers is assigned to a capsid missing a ring of three dimers in the center of a facet.

Collective translational and rotational Monte Carlo moves for attractive particles

Virtual move Monte Carlo is a Monte Carlo (MC) cluster algorithm forming clusters via local energy gradients and approximating the collective kinetic or dynamic motion of attractive colloidal particles. We carefully describe, analyze, and test the algorithm. To formally validate the algorithm through highlighting its symmetries, we present alternative and compact ways of selecting and accepting clusters which illustrate the formal use of abstract concepts in the design of biased MC techniques: the superdetailed balance and the early rejection scheme. A brief and comprehensive summary of the algorithms is presented, which makes them accessible without needing to understand the details of the derivation.

Simulations of HIV capsid protein dimerization reveal the effect of chemistry and topography on the mechanism of hydrophobic protein association

Recent work has shown that the hydrophobic protein surfaces in aqueous solution sit near a drying transition. The tendency for these surfaces to expel water from their vicinity leads to self-assembly of macromolecular complexes. In this article, we show with a realistic model for a biologically pertinent system how this phenomenon appears at the molecular level. We focus on the association of the C-terminal domain (CA-C) of the human immunodeficiency virus capsid protein. By combining all-atom simulations with specialized sampling techniques, we measure the water density distribution during the approach of two CA-C proteins as a function of separation and amino acid sequence in the interfacial region. The simulations demonstrate that CA-C protein-protein interactions sit at the edge of a dewetting transition and that this mesoscopic manifestation of the underlying liquid-vapor phase transition can be readily manipulated by biology or protein engineering to significantly affect association behavior. Although the wild-type protein remains wet until contact, we identify a set of in silico mutations, in which three hydrophilic amino acids are replaced with nonpolar residues, that leads to dewetting before association. The existence of dewetting depends on the size and relative locations of substituted residues separated by nanometer length scales, indicating long-range cooperativity and a sensitivity to surface topography. These observations identify important details that are missing from descriptions of protein association based on buried hydrophobic surface area.

Analyzing mechanisms and microscopic reversibility of self-assembly

We use computer simulations to investigate self-assembly in a system of model chaperonin proteins, and in an Ising lattice gas. We discuss the mechanisms responsible for rapid and efficient assembly in these systems, and we use measurements of dynamical activity and assembly progress to compare their propensities for kinetic trapping. We use the analytic solution of a simple minimal model to illustrate the key features associated with such trapping, paying particular attention to the number of ways that particles can misbind. We discuss the relevance of our results for the design and control of self-assembly in general.

Quantifying reversibility in a phase-separating lattice gas: an analogy with self-assembly

We present dynamic measurements of a lattice gas during phase separation, which we use as an analogy for self-assembly of equilibrium ordered structures. We use two approaches to quantify the degree of reversibility of this process: First, we count events in which bonds are made and broken; second, we use correlation-response measurements and fluctuation-dissipation ratios to probe reversibility during different time intervals. We show how correlation and response functions can be related directly to microscopic (ir)reversibility and we discuss the time dependence and observable dependence of these measurements, including the role of fast and slow degrees of freedom during assembly.

Effective confinement as origin of the equivalence of kinetic temperature and fluctuation-dissipation ratio in a dense shear-driven suspension

We study response and velocity autocorrelation functions for a tagged particle in a shear driven suspension governed by underdamped stochastic dynamics. We follow the idea of an effective confinement in dense suspensions and exploit a time scale separation between particle reorganization and vibrational motion. This allows us to approximately derive the fluctuation-dissipation theorem in a "hybrid" form involving the kinetic temperature as an effective temperature and an additive correction term. We show numerically that even in a moderately dense suspension the latter is negligible. We discuss similarities and differences with a simple toy model, a single trapped particle in shear flow.

Mechanisms of kinetic trapping in self-assembly and phase transformation

In self-assembly processes, kinetic trapping effects often hinder the formation of thermodynamically stable ordered states. In a model of viral capsid assembly and in the phase transformation of a lattice gas, we show how simulations in a self-assembling steady state can be used to identify two distinct mechanisms of kinetic trapping. We argue that one of these mechanisms can be adequately captured by kinetic rate equations, while the other involves a breakdown of theories that rely on cluster size as a reaction coordinate. We discuss how these observations might be useful in designing and optimising self-assembly reactions.

Encapsulation of a polymer by an icosahedral virus

The coat proteins of many viruses spontaneously form icosahedral capsids around nucleic acids or other polymers. Elucidating the role of the packaged polymer in capsid formation could promote biomedical efforts to block viral replication and enable use of capsids in nanomaterials applications. To this end, we perform Brownian dynamics on a coarse-grained model that describes the dynamics of icosahedral capsid assembly around a flexible polymer. We identify several mechanisms by which the polymer plays an active role in its encapsulation, including cooperative polymer-protein motions. These mechanisms are related to experimentally controllable parameters such as polymer length, protein concentration and solution conditions. Furthermore, the simulations demonstrate that assembly mechanisms are correlated with encapsulation efficiency, and we present a phase diagram that predicts assembly outcomes as a function of experimental parameters. We anticipate that our simulation results will provide a framework for designing in vitro assembly experiments on single-stranded RNA virus capsids.

Mechanisms of capsid assembly around a polymer

Capsids of many viruses assemble around nucleic acids or other polymers. Understanding how the properties of the packaged polymer affect the assembly process could promote biomedical efforts to prevent viral assembly or nanomaterials applications that exploit assembly. To this end, we simulate on a lattice the dynamical assembly of closed, hollow shells composed of several hundred to 1000 subunits, around a flexible polymer. We find that assembly is most efficient at an optimum polymer length that scales with the surface area of the capsid; polymers that are significantly longer than optimal often lead to partial-capsids with unpackaged polymer "tails" or a competition between multiple partial-capsids attached to a single polymer. These predictions can be tested with bulk experiments in which capsid proteins assemble around homopolymeric RNA or synthetic polyelectrolytes. We also find that the polymer can increase the net rate of subunit accretion to a growing capsid both by stabilizing the addition of new subunits and by enhancing the incoming flux of subunits; the effects of these processes may be distinguishable with experiments that monitor the assembly of individual capsids.

Control of pathways and yields of protein crystallization through the interplay of nonspecific and specific attractions

We use computer simulation to study crystal-forming model proteins equipped with interactions that are both orientationally specific and nonspecific. Distinct dynamical pathways of crystal formation can be selected by tuning the strengths of these interactions. When the nonspecific interaction is strong, liquidlike clustering can precede crystallization; when it is weak, growth can proceed via ordered nuclei. Crystal yields are in certain parameter regimes enhanced by the nonspecific interaction, even though it promotes association without local crystalline order. Our results suggest that equipping nanoscale components with weak nonspecific interactions (such as depletion attractions) can alter both their dynamical pathway of assembly and optimize the yield of the resulting material.

How do self-assembling polymers and gels age compared to glasses?

Experiments on gels have provided contradictory results concerning the relation between correlation and response functions during aging. To clarify this puzzle, we numerically investigate the fluctuation-dissipation plot in equilibrium polymers and in network forming gels employing two distinct observables, (i) the density Fourier transform and (ii) the single-particle potential energy, to probe (i) diffusional processes and (ii) the development of a bond network. The plot behaves very differently for the two cases. Violation from the equilibrium behavior is found only for the second observable. The experimental implications of the discovered sensitivity to the choice of the probe are discussed, in particular, with respect to the existing experimental results.

Predicting the effective temperature of a glass

We explain the findings by Di Leonardo et al. [Phys. Rev. Lett. 84, 6054 (2000)10.1103/PhysRevLett.84.6054] that the effective temperature of a Lennard-Jones glass depends only on the final density in the volume and/or temperature jump that produces the glass. This is not only a property of the Lennard-Jones liquid, but a feature of all strongly correlating liquids. For such liquids data from a single quench simulation provide enough information to predict the effective temperature of any glass produced by jumping from an equilibrium state. This prediction is validated by simulations of the Kob-Andersen binary Lennard-Jones liquid and shown not to apply for the nonstrongly correlating monatomic Lennard-Jones Gaussian liquid.

Self-assembly of monodisperse clusters: Dependence on target geometry

We apply a simple model system of patchy particles to study monodisperse self-assembly using the Platonic solids as target structures. We find marked differences between the assembly behaviors of the different systems. Tetrahedra, octahedral, and icosahedra assemble easily, while cubes are more challenging and dodecahedra do not assemble. We relate these differences to the kinetics and thermodynamics of assembly, with the formation of large disordered aggregates a particular important competitor to correct assembly. In particular, the free energy landscapes of those targets that are easy to assemble are funnel-like, whereas for the dodecahedral system the landscape is relatively flat with little driving force to facilitate escape from disordered aggregates.

Geometric frustration in small colloidal clusters

We study the structure of clusters in a model colloidal system with competing interactions using Brownian dynamics simulations. A short-ranged attraction drives clustering, while a weak, long-ranged repulsion is used to model electrostatic charging in experimental systems. The former is treated with a short-ranged Morse attractive interaction, the latter with a repulsive Yukawa interaction. We consider the yield of clusters of specific structure as a function of the strength of the interactions, for clusters with m = 3,4,5,6,7,10 and 13 colloids. At sufficient strengths of the attractive interaction (around 10k(B)T), the average bond lifetime approaches the simulation timescale and the system becomes nonergodic. For small clusters, m≤5, where geometric frustration is not relevant, despite nonergodicity, for sufficient strengths of the attractive interaction the yield of clusters which maximize the number of bonds approaches 100%. However for m = 7 and higher, in the nonergodic regime we find a lower yield of these structures where we argue geometric frustration plays a significant role. m = 6 is a special case, where two structures, of octahedral and C(2v) symmetry, compete, with the latter being favoured by entropic contributions in the ergodic regime and by kinetic trapping in the nonergodic regime. We believe that our results should be valid as long as the one-component description of the interaction potential is valid. A system with competing electrostatic repulsions and van der Waals attractions may be such an example. However, in some cases, the one-component description of the interaction potential may not be appropriate.

Invariant polymorphism in virus capsid assembly

Directed self-assembly of designed viral capsids holds significant potential for applications in materials science and medicine. However, the complexity of preparing these systems for assembly and the difficulty of quantitative experimental measurements on the assembly process have limited access to critical mechanistic questions that dictate the final product yields and isomorphic forms. Molecular simulations provide a means of elucidating self-assembly of viral proteins into icosahedral capsids and are the focus of the present study. Using geometrically realistic coarse-grained models with specialized molecular dynamics methods, we delineate conditions of temperature and coat protein concentration that lead to the spontaneous self-assembly of T = 1 and T = 3 icosahedral capsids. In addition to the primary product of icosahedral capsids, we observe a ubiquitous presence of nonicosahedral yet highly symmetric and enclosed aberrant capsules in both T = 1 and T = 3 systems. This polymorphism in assembly products recapitulates the scope and morphology of particle types that have been observed in mis-assembly experiments of virus capsids. Moreover, we find that this structural polymorphism in the end point structures is an inherent property of the coat proteins and arises from condition-dependent kinetic mechanisms that are independent of the elemental mechanisms of capsid growth (as long as the building blocks of the coat proteins are all monomeric, dimeric, or trimeric) and the capsid T number. The kinetic mechanisms responsible for self-assembly of icosahedral capsids and aberrant capsules are deciphered; the self-assembly of icosahedral capsids requires a high level of assembly fidelity, whereas self-assembly of nonicosahedral capsules is a consequence of an off-pathway mechanism that is prevalent under nonoptimal conditions of temperature or protein concentration during assembly. The latter case involves kinetically trapped dislocations of pentamer-templated proteins with hexameric organization. These findings provide insights into the complex processes that govern viral capsid assembly and suggest some features of the assembly process that can be exploited to control the assembly of icosahedral capsids and nonicosahedral capsules.

The role of collective motion in examples of coarsening and self-assembly

The simplest prescription for building a patterned structure from its constituents is to add particles, one at a time, to an appropriate template. However, self-organizing molecular and colloidal systems in nature can evolve in much more hierarchical ways. Specifically, constituents (or clusters of constituents) may aggregate to form clusters (or clusters of clusters) that serve as building blocks for later stages of assembly. Here we evaluate the character and consequences of such collective motion in a set of prototypical assembly processes. We do so using computer simulations in which a system's capacity for hierarchical dynamics can be controlled systematically. By explicitly allowing or suppressing collective motion, we quantify its effects. We find that coarsening within a two dimensional attractive lattice gas (and an analogous off-lattice model in three dimensions) is naturally dominated by collective motion over a broad range of temperatures and densities. Under such circumstances, cluster mobility inhibits the development of uniform coexisting phases, especially when macroscopic segregation is strongly favored by thermodynamics. By contrast, the assembly of model viral capsids is not frustrated but is instead facilitated by collective moves, which promote the orderly binding of intermediates consisting of several monomers.

Mechanisms of size control and polymorphism in viral capsid assembly

We simulate the assembly dynamics of icosahedral capsids from subunits that interconvert between different conformations (or quasi-equivalent states). The simulations identify mechanisms by which subunits form empty capsids with only one morphology but adaptively assemble into different icosahedral morphologies around nanoparticle cargoes with varying sizes, as seen in recent experiments with brome mosaic virus (BMV) capsid proteins. Adaptive cargo encapsidation requires moderate cargo-subunit interaction strengths; stronger interactions frustrate assembly by stabilizing intermediates with incommensurate curvature. We compare simulation results to experiments with cowpea chlorotic mottle virus empty capsids and BMV capsids assembled on functionalized nanoparticles and suggest new cargo encapsidation experiments. Finally, we find that both empty and templated capsids maintain the precise spatial ordering of subunit conformations seen in the crystal structure even if interactions that preserve this arrangement are favored by as little as the thermal energy, consistent with experimental observations that different subunit conformations are highly similar.

Controlling viral capsid assembly with templating

We develop coarse-grained models that describe the dynamic encapsidation of functionalized nanoparticles by viral capsid proteins. We find that some forms of cooperative interactions between protein subunits and nanoparticles can dramatically enhance rates and robustness of assembly, as compared to the spontaneous assembly of subunits into empty capsids. For large core-subunit interactions, subunits adsorb onto core surfaces en masse in a disordered manner, and then undergo a cooperative rearrangement into an ordered capsid structure. These assembly pathways are unlike any identified for empty capsid formation. Our models can be directly applied to recent experiments in which viral capsid proteins assemble around functionalized inorganic nanoparticles [Sun, Proc. Natl. Acad. Sci. U.S.A. 104, 1354 (2007)]. In addition, we discuss broader implications for understanding the dynamic encapsidation of single-stranded genomic molecules during viral replication and for developing multicomponent nanostructured materials.