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.

Hierarchical assembly is more robust than egalitarian assembly in synthetic capsids.

Self-assembly of complex and functional materials remains a grand challenge in soft material science. Efficient assembly depends on a delicate balance between thermodynamic and kinetic effects, requiring fine-tuning affinities and concentrations of subunits. By contrast, we introduce an assembly paradigm that allows large error-tolerance in the subunit affinity and helps avoid kinetic traps. Our combined experimental and computational approach uses a model system of triangular subunits programmed to assemble into T = 3 icosahedral capsids comprising 60 units. The experimental platform uses DNA origami to create monodisperse colloids whose three-dimensional geometry is controlled to nanometer precision, with two distinct bonds whose affinities are controlled to kBT precision, quantified in situ by static light scattering. The computational model uses a coarse-grained representation of subunits, short-ranged potentials, and Langevin dynamics. Experimental observations and modeling reveal that when the bond affinities are unequal, two distinct hierarchical assembly pathways occur, in which the subunits first form dimers in one case and pentamers in another. These hierarchical pathways produce complete capsids faster and are more robust against affinity variation than egalitarian pathways, in which all binding sites have equal strengths. This finding suggests that hierarchical assembly may be a general engineering principle for optimizing self-assembly of complex target structures.

A machine learning approach to robustly determine director fields and analyze defects in active nematics.

Active nematics are dense systems of rodlike particles that consume energy to drive motion at the level of the individual particles. They exist in natural systems like biological tissues and artificial materials such as suspensions of self-propelled colloidal particles or synthetic microswimmers. Active nematics have attracted significant attention in recent years due to their spectacular nonequilibrium collective spatiotemporal dynamics, which may enable applications in fields such as robotics, drug delivery, and materials science. The director field, which measures the direction and degree of alignment of the local nematic orientation, is a crucial characteristic of active nematics and is essential for studying topological defects. However, determining the director field is a significant challenge in many experimental systems. Although director fields can be derived from images of active nematics using traditional imaging processing methods, the accuracy of such methods is highly sensitive to the settings of the algorithms. These settings must be tuned from image to image due to experimental noise, intrinsic noise of the imaging technology, and perturbations caused by changes in experimental conditions. This sensitivity currently limits automatic analysis of active nematics. To address this, we developed a machine learning model for extracting reliable director fields from raw experimental images, which enables accurate analysis of topological defects. Application of the algorithm to experimental data demonstrates that the approach is robust and highly generalizable to experimental settings that are different from those in the training data. It could be a promising tool for investigating active nematics and may be generalized to other active matter systems.

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.

From disks to channels: dynamics of active nematics confined to an annulus.

Confinement can be used to systematically tame turbulent dynamics occurring in active fluids. Although periodic channels are the simplest geometries to study confinement numerically, the corresponding experimental realizations require closed racetracks. Here, we computationally study 2D active nematics confined to such a geometry-an annulus. By systematically varying the annulus inner radius and channel width, we bridge the behaviors observed in the previously studied asymptotic limits of the annulus geometry: a disk and an infinite channel. We identify new steady-state behaviors, which reveal the influence of boundary curvature and its interplay with confinement. We also show that, below a threshold inner radius, the dynamics are insensitive to the presence of the inner hole. We explain this insensitivity through a simple scaling analysis. Our work sheds further light on design principles for using confinement to control the dynamics of active nematics.

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.

Engineering Metastability into a Virus-like Particle to Enable Triggered Dissociation.

For a virus-like particle (VLP) to serve as a delivery platform, the VLP must be able to release its cargo in response to a trigger. Here, we use a chemical biology approach to destabilize a self-assembling capsid for a subsequent triggered disassembly. We redesigned the dimeric hepatitis B virus (HBV) capsid protein (Cp) with two differentially addressable cysteines, C150 for reversibly crosslinking the capsid and C124 to react with a destabilizing moiety. The resulting construct, Cp150-V124C, assembles into icosahedral, 120-dimer VLPs that spontaneously crosslink via the C-terminal C150, leaving C124 buried at a dimer-dimer interface. The VLP is driven into a metastable state when C124 is reacted with the bulky fluorophore, maleimidyl BoDIPY-FL. The resulting VLP is stable until exposed to modest, physiologically relevant concentrations of reducing agent. We observe dissociation with FRET relaxation of polarization, size exclusion chromatography, and resistive-pulse sensing. Dissociation is slow, minutes to hours, with a characteristic lag phase. Mathematical modeling based on the presence of a nucleation step predicts disassembly dynamics that are consistent with experimental observations. VLPs transfected into hepatoma cells show similar dissociation behavior. These results suggest a generalizable strategy for designing a VLP that can release its contents in an environmentally responsive reaction.

Optimization of non-equilibrium self-assembly protocols using Markov state models.

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.

Geometrically programmed self-limited assembly of tubules using DNA origami colloids.

Self-assembly is one of the most promising strategies for making functional materials at the nanoscale, yet new design principles for making self-limiting architectures, rather than spatially unlimited periodic lattice structures, are needed. To address this challenge, we explore the tradeoffs between addressable assembly and self-closing assembly of a specific class of self-limiting structures: cylindrical tubules. We make triangular subunits using DNA origami that have specific, valence-limited interactions and designed binding angles, and we study their assembly into tubules that have a self-limited width that is much larger than the size of an individual subunit. In the simplest case, the tubules are assembled from a single component by geometrically programming the dihedral angles between neighboring subunits. We show that the tubules can reach many micrometers in length and that their average width can be prescribed through the dihedral angles. We find that there is a distribution in the width and the chirality of the tubules, which we rationalize by developing a model that considers the finite bending rigidity of the assembled structure as well as the mechanism of self-closure. Finally, we demonstrate that the distributions of tubules can be further sculpted by increasing the number of subunit species, thereby increasing the assembly complexity, and demonstrate that using two subunit species successfully reduces the number of available end states by half. These results help to shed light on the roles of assembly complexity and geometry in self-limited assembly and could be extended to other self-limiting architectures, such as shells, toroids, or triply periodic frameworks.

Multiscale Modeling of Hepatitis B Virus Capsid Assembly and Its Dimorphism.

Hepatitis B virus (HBV) is an endemic, chronic virus that leads to 800000 deaths per year. Central to the HBV lifecycle, the viral core has a protein capsid assembled from many copies of a single protein. The capsid protein adopts different (quasi-equivalent) conformations to form icosahedral capsids containing 180 or 240 proteins: T = 3 or T = 4, respectively, in Caspar-Klug nomenclature. HBV capsid assembly has become an important target for recently developed antivirals; nonetheless, the assembly pathways and mechanisms that control HBV dimorphism remain unclear. We describe computer simulations of the HBV assembly, using a coarse-grained model that has parameters learned from all-atom molecular dynamics simulations of a complete HBV capsid and yet is computationally tractable. Dynamical simulations with the resulting model reproduce experimental observations of HBV assembly pathways and products. By constructing Markov state models and employing transition path theory, we identify pathways leading to T = 3, T = 4, and other experimentally observed capsid morphologies. The analysis shows that capsid polymorphism is promoted by the low HBV capsid bending modulus, where the key factors controlling polymorphism are the conformational energy landscape and protein-protein binding affinities.

Polymorphic self-assembly of helical tubules is kinetically controlled.

In contrast to most self-assembling synthetic materials, which undergo unbounded growth, many biological self-assembly processes are self-limited. That is, the assembled structures have one or more finite dimensions that are much larger than the size scale of the individual monomers. In many such cases, the finite dimension is selected by a preferred curvature of the monomers, which leads to self-closure of the assembly. In this article, we study an example class of self-closing assemblies: cylindrical tubules that assemble from triangular monomers. By combining kinetic Monte Carlo simulations, free energy calculations, and simple theoretical models, we show that a range of programmable size scales can be targeted by controlling the intricate balance between the preferred curvature of the monomers and their interaction strengths. However, their assembly is kinetically controlled-the tubule morphology is essentially fixed shortly after closure, resulting in a distribution of tubule widths that is significantly broader than the equilibrium distribution. We develop a simple kinetic model based on this observation and the underlying free-energy landscape of assembling tubules that quantitatively describes the distributions. Our results are consistent with recent experimental observations of tubule assembly from triangular DNA origami monomers. The modeling framework elucidates design principles for assembling self-limited structures from synthetic components, such as artificial microtubules that have a desired width and chirality.

Microcompartment assembly around multicomponent fluid cargoes.

This article describes dynamical simulations of the assembly of an icosahedral protein shell around a bicomponent fluid cargo. Our simulations are motivated by bacterial microcompartments, which are protein shells found in bacteria that assemble around a complex of enzymes and other components involved in certain metabolic processes. The simulations demonstrate that the relative interaction strengths among the different cargo species play a key role in determining the amount of each species that is encapsulated, their spatial organization, and the nature of the shell assembly pathways. However, the shell protein-shell protein and shell protein-cargo component interactions that help drive assembly and encapsulation also influence cargo composition within certain parameter regimes. These behaviors are governed by a combination of thermodynamic and kinetic effects. In addition to elucidating how natural microcompartments encapsulate multiple components involved within reaction cascades, these results have implications for efforts in synthetic biology to colocalize alternative sets of molecules within microcompartments to accelerate specific reactions. More broadly, the results suggest that coupling between self-assembly and multicomponent liquid-liquid phase separation may play a role in the organization of the cellular cytoplasm.

Thermodynamic Size Control in Curvature-Frustrated Tubules: Self-Limitation with Open Boundaries.

We use computational modeling to investigate the assembly thermodynamics of a particle-based model for geometrically frustrated assembly, in which the local packing geometry of subunits is incompatible with uniform, strain-free large-scale assembly. The model considers discrete triangular subunits that drive assembly toward a closed, hexagonal-ordered tubule, but have geometries that locally favor negative Gaussian curvature. We use dynamical Monte Carlo simulations and enhanced sampling methods to compute the free energy landscape and corresponding self-assembly behavior as a function of experimentally accessible parameters that control assembly driving forces and the magnitude of frustration. The results determine the parameter range where finite-temperature self-limiting assembly occurs, in which the equilibrium assembly size distribution is sharply peaked around a well-defined finite size. The simulations also identify two mechanisms by which the system can escape frustration and assemble to unlimited size, and determine the particle-scale properties of subunits that suppress unbounded growth.

Tiling a tubule: how increasing complexity improves the yield of self-limited assembly.

The ability to design and synthesize ever more complicated colloidal particles opens the possibility of self-assembling a zoo of complex structures, including those with one or more self-limited length scales. An undesirable feature of systems with self-limited length scales is that thermal fluctuations can lead to the assembly of nearby, off-target states. We investigate strategies for limiting off-target assembly by using multiple types of subunits. Using simulations and energetics calculations, we explore this concept by considering the assembly of tubules built from triangular subunits that bind edge to edge. While in principle, a single type of triangle can assemble into tubules with a monodisperse width distribution, in practice, the finite bending rigidity of the binding sites leads to the formation of off-target structures. To increase the assembly specificity, we introduce tiling rules for assembling tubules from multiple species of triangles. We show that the selectivity of the target structure can be dramatically improved by using multiple species of subunits, and provide a prescription for choosing the minimum number of subunit species required for near-perfect yield. Our approach of increasing the system's complexity to reduce the accessibility of neighboring structures should be generalizable to other systems beyond the self-assembly of tubules.

Data-Driven Discovery of Active Nematic Hydrodynamics.

Active nematics can be modeled using phenomenological continuum theories that account for the dynamics of the nematic director and fluid velocity through partial differential equations (PDEs). While these models provide a statistical description of the experiments, the relevant terms in the PDEs and their parameters are usually identified indirectly. We adapt a recently developed method to automatically identify optimal continuum models for active nematics directly from spatiotemporal data, via sparse regression of the coarse-grained fields onto generic low order PDEs. After extensive benchmarking, we apply the method to experiments with microtubule-based active nematics, finding a surprisingly minimal description of the system. Our approach can be generalized to gain insights into active gels, microswimmers, and diverse other experimental active matter systems.

Vesicle shape transformations driven by confined active filaments.

In active matter systems, deformable boundaries provide a mechanism to organize internal active stresses. To study a minimal model of such a system, we perform particle-based simulations of an elastic vesicle containing a collection of polar active filaments. The interplay between the active stress organization due to interparticle interactions and that due to the deformability of the confinement leads to a variety of filament spatiotemporal organizations that have not been observed in bulk systems or under rigid confinement, including highly-aligned rings and caps. In turn, these filament assemblies drive dramatic and tunable transformations of the vesicle shape and its dynamics. We present simple scaling models that reveal the mechanisms underlying these emergent behaviors and yield design principles for engineering active materials with targeted shape dynamics.

Multiscale Microtubule Dynamics in Active Nematics.

In microtubule-based active nematics, motor-driven extensile motion of microtubule bundles powers chaotic large-scale dynamics. We quantify the interfilament sliding motion both in isolated bundles and in a dense active nematic. The extension speed of an isolated microtubule pair is comparable to the molecular motor stepping speed. In contrast, the net extension in dense 2D active nematics is significantly slower; the interfilament sliding speeds are widely distributed about the average and the filaments exhibit both contractile and extensile relative motion. These measurements highlight the challenge of connecting the extension rate of isolated bundles to the multimotor and multifilament interactions present in a dense 2D active nematic. They also provide quantitative data that is essential for building multiscale models.

Active liquid crystals powered by force-sensing DNA-motor clusters.

Cytoskeletal active nematics exhibit striking nonequilibrium dynamics that are powered by energy-consuming molecular motors. To gain insight into the structure and mechanics of these materials, we design programmable clusters in which kinesin motors are linked by a double-stranded DNA linker. The efficiency by which DNA-based clusters power active nematics depends on both the stepping dynamics of the kinesin motors and the chemical structure of the polymeric linker. Fluorescence anisotropy measurements reveal that the motor clusters, like filamentous microtubules, exhibit local nematic order. The properties of the DNA linker enable the design of force-sensing clusters. When the load across the linker exceeds a critical threshold, the clusters fall apart, ceasing to generate active stresses and slowing the system dynamics. Fluorescence readout reveals the fraction of bound clusters that generate interfilament sliding. In turn, this yields the average load experienced by the kinesin motors as they step along the microtubules. DNA-motor clusters provide a foundation for understanding the molecular mechanism by which nanoscale molecular motors collectively generate mesoscopic active stresses, which in turn power macroscale nonequilibrium dynamics of active nematics.

Avidity and surface mobility in multivalent ligand-receptor binding.

Targeted drug delivery relies on two physical processes: the selective binding of a therapeutic particle to receptors on a specific cell membrane, followed by transport of the particle across the membrane. In this article, we address some of the challenges in controlling the thermodynamics and dynamics of these two processes by combining a simple experimental system with a statistical mechanical model. Specifically, we characterize and model multivalent ligand-receptor binding between colloidal particles and fluid lipid bilayers, as well as the surface mobility of membrane-bound particles. We show that the mobility of the receptors within the fluid membrane is key to both the thermodynamics and dynamics of binding. First, we find that the particle-membrane binding free energy-or avidity-is a strongly nonlinear function of the ligand-receptor affinity. We attribute the nonlinearity to a combination of multivalency and recruitment of fluid receptors to the binding site. Our results also suggest that partial wrapping of the bound particles by the membrane enhances avidity further. Second, we demonstrate that the lateral mobility of membrane-bound particles is also strongly influenced by the recruitment of receptors. Specifically, we find that the lateral diffusion coefficient of a membrane-bound particle is dominated by the hydrodynamic drag against the aggregate of receptors within the membrane. These results provide one of the first direct validations of the working theoretical framework for multivalent interactions. They also highlight that the fluidity and elasticity of the membrane are as important as the ligand-receptor affinity in determining the binding and transport of small particles attached to membranes.

Programmable icosahedral shell system for virus trapping.

Broad-spectrum antiviral platforms that can decrease or inhibit viral infection would alleviate many threats to global public health. Nonetheless, effective technologies of this kind are still not available. Here, we describe a programmable icosahedral canvas for the self-assembly of icosahedral shells that have viral trapping and antiviral properties. Programmable triangular building blocks constructed from DNA assemble with high yield into various shell objects with user-defined geometries and apertures. We have created shells with molecular masses ranging from 43 to 925 MDa (8 to 180 subunits) and with internal cavity diameters of up to 280 nm. The shell interior can be functionalized with virus-specific moieties in a modular fashion. We demonstrate this virus-trapping concept by engulfing hepatitis B virus core particles and adeno-associated viruses. We demonstrate the inhibition of hepatitis B virus core interactions with surfaces in vitro and the neutralization of infectious adeno-associated viruses exposed to human cells.