Viral self-assembly as a thermodynamic process

Nanomechanical resilience and thermal stability of RSJ2 phage

As the world moves toward a green economy and sustainable agriculture, bacterial viruses or bacteriophages (phages) become attractive biocontrol agents for controlling crop diseases. Effective utilization of phages in farms requires integrated knowledge of crops, pathogens, phages, and surroundings. Phages must encounter environmental fluctuations, including temperature, and must remain infectious for successful bacteria lysis. This work studied a soilborne RSJ2 phage discovered in Thailand, which can eliminate Ralstonia solanacearum, causing bacterial wilt disease in chili. We investigated how phage infectivity and nanomechanics responded to thermal changes. The plaque-based assay showed that the infectivity of the RSJ2 phage was stable within 24-40 °C, an average temperature fluctuation in tropical regions. The structural examination also showed that the phage remained intact. The nanomechanical property of the phage was inspected by the atomic force microscopy-based nanoindentation. The result revealed that the phage stiffness within 24-40 °C was statistically similar (0.05-0.06 N/m). Upon heating at 40 °C for 1, 5, and 10 h and resting at 25 °C, the stiffness of the phage particles increased to 0.09-0.11 N/m (54-83% increase). The stiffness results suggest structural adaptation of the protein subunits as a response to thermal alteration. The study exhibits that the phage structure is highly dynamic and can nanomechanically respond to varying temperatures. The phage stiffness may reveal insight into phage adaptation to environmental factors. Equipped with the knowledge of phage infectivity, structure, and nanomechanics, we can design practical guidelines for effective phage usage in farming and propelling green and safe agriculture.

Virtual indentation of the empty capsid of the minute virus of mice using a minimal coarse-grained model

A minimal coarse-grained model for T=1 viral capsids assembled from 20 protein rigid trimers has been designed by extending a previously proposed form of the interaction energy written as a sum of anisotropic pairwise interactions between the trimeric capsomers. The extension of the model has been performed to properly account for the coupling between two internal coordinates: the one that measures the intercapsomer distance and the other that gives the intercapsomer dihedral angle. The model has been able to fit with less than a 10% error the atomic force microscopy (AFM) indentation experimental data for the empty capsid of the minute virus of mice (MVM), providing in this way an admissible picture of the main mechanisms behind the capsid deformations. In this scenario, the bending of the intercapsomer dihedral angle is the angular internal coordinate that can support larger deformations away from its equilibrium values, determining important features of the AFM indentation experiments as the elastic constants along the three symmetry axes of the capsid and the critical indentations. From the value of one of the parameters of our model, we conclude that trimers in the MVM must be quite oblate tops, in excellent agreement with their known structure. The transition from the linear to the nonlinear regimes sampled in the indentation process appears to be an interesting topic for future research in physical virology.

Close packings of identical proteins in small spherical capsids and similar proteinaceous shells

Understanding the principles governing protein arrangement in viral capsids and structurally similar protein shells can enable the development of new antiviral strategies and the design of artificial protein cages for various applications. We study these principles within the context of the close packing problem, by analyzing dozens of small spherical shells assembled from a single type of protein. First, we use icosahedral spherical close packings containing 60T identical disks, where T ≤ 4, to rationalize the protein arrangement in twenty real icosahedral shells both satisfying and violating the paradigmatic Caspar-Klug model. We uncover a striking correspondence between the protein mass centers in the considered shells and the centers of disks in the close packings. To generalize the packing model, we consider proteins with a weak shape anisotropy and propose an interaction energy, minimization of which allows us to obtain spherical dense packings of slightly anisotropic structural units. In the case of strong anisotropy, we model the proteins as sequences of overlapping discs of different sizes, with minimum energy configuration not only resulting in packings, accurately reproducing locations and orientations of individual proteins, but also revealing that icosahedral packings that display the handedness of real capsids are energetically more favorable. Finally, by introducing effective disc charges, we rationalize the formation of inter-protein bonds in protein shells.

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.

Relaxational dynamics of the T-number conversion of virus capsids

We extend a recently proposed kinetic theory of virus capsid assembly based on Model A kinetics and study the dynamics of the interconversion of virus capsids of different sizes triggered by a quench, that is, by sudden changes in the solution conditions. The work is inspired by in vitro experiments on functionalized coat proteins of the plant virus cowpea chlorotic mottle virus, which undergo a reversible transition between two different shell sizes (T = 1 and T = 3) upon changing the acidity and salinity of the solution. We find that the relaxation dynamics are governed by two time scales that, in almost all cases, can be identified as two distinct processes. Initially, the monomers and one of the two types of capsids respond to the quench. Subsequently, the monomer concentration remains essentially constant, and the conversion between the two capsid species completes. In the intermediate stages, a long-lived metastable steady state may present itself, where the thermodynamically less stable species predominate. We conclude that a Model A based relaxational model can reasonably describe the early and intermediate stages of the conversion experiments. However, it fails to provide a good representation of the time evolution of the state of assembly of the coat proteins in the very late stages of equilibration when one of the two species disappears from the solution. It appears that explicitly incorporating the nucleation barriers to assembly and disassembly is crucial for an accurate description of the experimental findings, at least under conditions where these barriers are sufficiently large.

Packing and trimer-to-dimer protein reconstruction in icosahedral viral shells with a single type of symmetrical structural unit

Understanding the principles of protein packing and the mechanisms driving morphological transformations in virus shells (capsids) during their maturation can be pivotal for the development of new antiviral strategies. Here, we study how these principles and mechanisms manifest themselves in icosahedral viral capsids assembled from identical symmetric structural units (capsomeres). To rationalize such shells, we model capsomers as symmetrical groups of identical particles interacting with a short-range potential typical of the classic Tammes problem. The capsomere particles are assumed to retain their relative positions on the vertices of planar polygons placed on the spherical shell and to interact only with the particles from other capsomeres. Minimization of the interaction energy enforces equal distances between the nearest particles belonging to neighboring capsomeres and minimizes the number of different local environments. Thus, our model implements the Caspar and Klug quasi-equivalence principle and leads to packings strikingly similar to real capsids. We then study a reconstruction of protein trimers into dimers in a Flavivirus shell during its maturation, connecting the relevant structural changes with the modifications of the electrostatic charges of proteins, wrought by the oxidative switch in the bathing solution that is essential for the process. We highlight the key role of pr peptides in the shell reconstruction and show that the highly ordered arrangement of these subunits in the dimeric state is energetically favored at a low pH level. We also discuss the electrostatic mechanisms controlling the release of pr peptides in the last irreversible step of the maturation process.

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.

Inositol Hexakisphosphate (IP6) Accelerates Immature HIV-1 Gag Protein Assembly toward Kinetically Trapped Morphologies

During the late stages of the HIV-1 lifecycle, immature virions are produced by the concerted activity of Gag polyproteins, primarily mediated by the capsid (CA) and spacer peptide 1 (SP1) domains, which assemble into a spherical lattice, package viral genomic RNA, and deform the plasma membrane. Recently, inositol hexakisphosphate (IP6) has been identified as an essential assembly cofactor that efficiently produces both immature virions in vivo and immature virus-like particles in vitro. To date, however, several distinct mechanistic roles for IP6 have been proposed on the basis of independent functional, structural, and kinetic studies. In this work, we investigate the molecular influence of IP6 on the structural outcomes and dynamics of CA/SP1 assembly using coarse-grained (CG) molecular dynamics (MD) simulations and free energy calculations. Here, we derive a bottom-up, low-resolution, and implicit-solvent CG model of CA/SP1 and IP6, and simulate their assembly under conditions that emulate both in vitro and in vivo systems. Our analysis identifies IP6 as an assembly accelerant that promotes curvature generation and fissure-like defects throughout the lattice. Our findings suggest that IP6 induces kinetically trapped immature morphologies, which may be physiologically important for later stages of viral morphogenesis and potentially useful for virus-like particle technologies.

Local rules for the self-assembly of a non-quasi-equivalent viral capsid

The structures of many large bacteriophages, such as the P23-77 capsids, do not adhere strictly to the quasi-equivalence principle of viral architecture. Although the general architecture of the P23-77 capsids is classed as T=28d, it self-assembles from multiple copies of two types of coat protein subunits, and the resulting hexameric capsomers do not conform to the Caspar-Klug paradigm. There are two types of hexamers with distinct internal organization, that are located at specific positions in the capsid. It is an open problem which assembly mechanism can lead to such a complex capsid organization. Here we propose a simple set of local rules that can explain how such non-quasi-equivalent capsid structures can arise as a result of self-assembly.

Mechanical Capsid Maturation Facilitates the Resolution of Conflicting Requirements for Herpesvirus Assembly

Most viruses undergo a maturation process from a weakly self-assembled, noninfectious particle to a stable, infectious virion. For herpesviruses, this maturation process resolves several conflicting requirements: (i) assembly must be driven by weak, reversible interactions between viral particle subunits to reduce errors and minimize the energy of self-assembly, and (ii) the viral particle must be stable enough to withstand tens of atmospheres of DNA pressure resulting from its strong confinement in the capsid. With herpes simplex virus 1 (HSV-1) as a prototype of human herpesviruses, we demonstrated that this mechanical capsid maturation is mainly facilitated through capsid binding auxiliary protein UL25, orthologs of which are present in all herpesviruses. Through genetic manipulation of UL25 mutants of HSV-1 combined with the interrogation of capsid mechanics with atomic force microscopy nano-indentation, we suggested the mechanism of stepwise binding of distinct UL25 domains correlated with capsid maturation and DNA packaging. These findings demonstrate another paradigm of viruses as elegantly programmed nano-machines where an intimate relationship between mechanical and genetic information is preserved in UL25 architecture.

IMPORTANCE The minor capsid protein UL25 plays a critical role in the mechanical maturation of the HSV-1 capsid during virus assembly and is required for stable DNA packaging. We modulated the UL25 capsid interactions by genetically deleting different UL25 regions and quantifying the effect on mechanical capsid stability using an atomic force microscopy (AFM) nanoindentation approach. This approach revealed how UL25 regions reinforced the herpesvirus capsid to stably package and retain pressurized DNA. Our data suggest a mechanism of stepwise binding of two main UL25 domains timed with DNA packaging.

Phase transitions in virology

Viruses have established relationships with almost every other living organism on Earth and at all levels of biological organization: from other viruses up to entire ecosystems. In most cases, they peacefully coexist with their hosts, but in most relevant cases, they parasitize them and induce diseases and pandemics, such as the AIDS and the most recent avian influenza and COVID-19 pandemic events, causing a huge impact on health, society, and economy. Viruses play an essential role in shaping the eco-evolutionary dynamics of their hosts, and have been also involved in some of the major evolutionary innovations either by working as vectors of genetic information or by being themselves coopted by the host into their genomes. Viruses can be studied at different levels of biological organization, from the molecular mechanisms of genome replication, gene expression and encapsidation, to global pandemics. All these levels are different and yet connected through the presence of threshold conditions allowing for the formation of a capsid, the loss of genetic information or epidemic spreading. These thresholds, as occurs with temperature separating phases in a liquid, define sharp qualitative types of behaviour. Thesephase transitionsare very well known in physics. They have been studied by means of simple, but powerful models able to capture their essential properties, allowing us to better understand them. Can the physics of phase transitions be an inspiration for our understanding of viral dynamics at different scales? Here we review well-known mathematical models of transition phenomena in virology. We suggest that the advantages of abstract, simplified pictures used in physics are also the key to properly understanding the origins and evolution of complexity in viruses. By means of several examples, we explore this multilevel landscape and how minimal models provide deep insights into a diverse array of problems. The relevance of these transitions in connecting dynamical patterns across scales and their evolutionary and clinical implications are outlined.

Self-assembly of polyhedral bilayer vesicles from Piezo ion channels

Piezo ion channels underlie many forms of mechanosensation in vertebrates and have been found to bend the membrane into strongly curved dome shapes. We develop a methodology describing the self-assembly of lipids and Piezo proteins into polyhedral bilayer vesicles. We validate this methodology for bilayer vesicles formed from bacterial mechanosensitive channels of small conductance, for which experiments found a polyhedral arrangement of proteins with snub cube symmetry and a well-defined characteristic vesicle size. On this basis, we calculate the self-assembly diagram for polyhedral bilayer vesicles formed from Piezo proteins. We find that the radius of curvature of the Piezo dome provides a critical control parameter for the self-assembly of Piezo vesicles, with high abundances of Piezo vesicles with octahedral, icosahedral, and snub cube symmetry with increasing Piezo dome radius of curvature.

Minimal Design Principles for Icosahedral Virus Capsids

The geometrical structures of single- and multiple-shell icosahedral virus capsids are reproduced as the targets that minimize the cost corresponding to relatively simple design functions. Capsid subunits are first identified as building blocks at a given coarse-grained scale and then represented in these functions as point particles located on an appropriate number of concentric spherical surfaces. Minimal design cost is assigned to optimal spherical packings of the particles. The cost functions are inspired by the packings favored for the Thomson problem, which minimize the electrostatic potential energy between identical charged particles. In some cases, icosahedral symmetry constraints are incorporated as external fields acting on the particles. The simplest cost functions can be obtained by separating particles in disjoint nonequivalent sets with distinct interactions, or by introducing interacting holes (the absence of particles). These functions can be adapted to reproduce any capsid structure found in real viruses. Structures absent in Nature require significantly more complex designs. Measures of information content and complexity are assigned to both the cost functions and the capsid geometries. In terms of these measures, icosahedral structures and the corresponding cost functions are the simplest solutions.

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.

Buckling transitions and soft-phase invasion of two-component icosahedral shells

What is the optimal distribution of two types of crystalline phases on the surface of icosahedral shells, such as of many viral capsids? We here investigate the distribution of a thin layer of soft material on a crystalline convex icosahedral shell. We demonstrate how the shapes of spherical viruses can be understood from the perspective of elasticity theory of thin two-component shells. We develop a theory of shape transformations of an icosahedral shell upon addition of a softer, but still crystalline, material onto its surface. We show how the soft component "invades" the regions with the highest elastic energy and stress imposed by the 12 topological defects on the surface. We explore the phase diagram as a function of the surface fraction of the soft material, the shell size, and the incommensurability of the elastic moduli of the rigid and soft phases. We find that, as expected, progressive filling of the rigid shell by the soft phase starts from the most deformed regions of the icosahedron. With a progressively increasing soft-phase coverage, the spherical segments of domes are filled first (12 vertices of the shell), then the cylindrical segments connecting the domes (30 edges) are invaded, and, ultimately, the 20 flat faces of the icosahedral shell tend to be occupied by the soft material. We present a detailed theoretical investigation of the first two stages of this invasion process and develop a model of morphological changes of the cone structure that permits noncircular cross sections. In conclusion, we discuss the biological relevance of some structures predicted from our calculations, in particular for the shape of viral capsids.

Viruses as 'Truffle Hounds': Molecular Tools for Untangling Brain Cellular Pathology

The ability of viruses to evolve several orders of magnitude faster than their host cells has enabled them to exploit host cellular machinery by selectively recruiting multiprotein complexes (MPCs) for their catalyzed assembly and replication. This hijacking may depend on alternative, 'moonlighting' functions of host proteins that deviate from their canonical functions thereby inducing cellular pathology. Here, we posit that if virus-induced cellular pathology is similar to that of other, unknown (non-viral) causes, the identification and molecular characterization of the host proteins involved in virus-mediated cellular pathology can be leveraged to decipher the non-viral disease-relevant mechanisms. We focus on how virus-induced aberrant proteostasis and protein aggregation resemble the cellular pathology of sporadic neurodegenerative diseases (NDs) and how this can be exploited for drug discovery.

Effect of protein steric constraints on the symmetry of membrane protein polyhedra

Experiments have shown that, in an aqueous environment, lipids and membrane proteins can self-assemble into membrane protein polyhedral nanoparticles (MPPNs). MPPNs are closed, spherical vesicles composed of a lipid bilayer membrane and membrane proteins, with a polyhedral arrangement of membrane proteins. The observed symmetry and size of MPPNs can be understood from the interplay of protein-induced lipid bilayer deformations in MPPNs, topological defects in protein packing necessitated by the spherical shape of MPPNs, and thermal fluctuations in MPPN self-assembly. We explore here the effect of protein steric constraints on MPPN shape. The protein steric constraints considered here may arise from a well-defined shape of protein domains outside the membrane, entropic repulsion between membrane proteins with flexible domains outside the membrane, or binding of other molecules to membrane proteins. Calculating MPPN self-assembly diagrams under protein steric constraints we find that protein steric constraints can strongly affect MPPN self-assembly. Depending on the specific scenario considered, protein steric constraints can leave large portions of the MPPN self-assembly diagrams with no clearly defined MPPN symmetry or substantially expand the regions of MPPN self-assembly diagrams dominated by highly symmetric MPPN states, such as MPPNs with icosahedral or snub cube symmetry. Our results suggest that modification of protein steric constraints may allow the directed self-assembly of MPPNs with specified symmetry, size, and protein composition and may thus facilitate the further utilization of MPPNs for membrane protein structural analysis or targeted drug delivery.

Should Virus Capsids Assemble Perfectly? Theory and Observation of Defects

Although published structural models of viral capsids generally exhibit a high degree of regularity or symmetry, structural defects might be expected because of the fluctuating environment in which capsids assemble and the requirement of some capsids for disassembly before genome delivery. Defective structures are observed in computer simulations, and are evident in single-particle cryoelectron microscopy studies. Here, we quantify the conditions under which defects might be expected, using a statistical mechanics model allowing for ideal, defective, and vacant sites. The model displays a threshold in affinity parameters below which there is an appreciable population of defective capsids. Even when defective sites are not allowed, there is generally some population of vacancies. Analysis of single particles in cryoelectron microscopy micrographs yields a confirmatory ≳15% of defective particles. Our findings suggest structural heterogeneity in virus capsids may be under-appreciated, and also points to a nontraditional strategy for assembly inhibition.

Shape morphology of dipolar domains in planar and spherical monolayers

We present a continuum theory for predicting the equilibrium shape and size of dipolar domains formed during liquid-liquid phase coexistence in planar and spherical monolayers. Our main objective is to assess the impact of the monolayer surface curvature on domain morphology. Following previous investigators, we base our analysis around minimizing the free energy, with contributions from line tension and electrostatic dipolar repulsions. Assuming a monodisperse system of circularly symmetric domains, we calculate self-energies and interaction energies for planar and spherical monolayers and determine the equilibrium domain size from the energy minima. We subsequently evaluate the stability of the circularly symmetric domain shapes to an arbitrary, circumferential distortion of the perimeter via a linear stability analysis. We find that the surface curvature generally promotes the formation of smaller, circularly symmetric domains instead of larger, elongated domains. We rationalize these results by examining the effect of the curvature on the intra- and inter-domain dipolar repulsions. We then present a phase diagram of domain shape morphologies, parameterized in terms of the domain area fraction and the monolayer curvature. For typical domain dimensions of 1-30 µm, our theoretical results are relevant to monolayers (and possibly also bilayers) in liquid-liquid phase coexistence with radii of curvature of 1-100 µm.

Crystal-like order and defects in metazoan epithelia with spherical geometry

Since Robert Hooke studied cork cell patterns in 1665, scientists have been puzzled by why cells form such ordered structures. The laws underlying this type of organization are universal, and we study them comparing the living and non-living two-dimensional systems self-organizing at the spherical surface. Such-type physical systems often possess trigonal order with specific elongated defects, scars and pleats, where the 5-valence and 7-valence vertices alternate. In spite of the fact that the same physical and topological rules are involved in the structural organization of biological systems, such topological defects were never reported in epithelia. We have discovered them in the follicular spherical epithelium of ascidians that are emerging models in developmental biology. Surprisingly, the considered defects appear in the epithelium even when the number of cells in it is significantly less than the previously known threshold value. We explain this result by differences in the cell sizes and check our hypothesis considering the self-assembly of different random size particles on the spherical surface. Scars, pleats and other complex defects found in ascidian samples can play an unexpected and decisive role in the permanent renewal and reorganization of epithelia, which forms or lines many tissues and organs in metazoans.

A minimal coarse-grained model for the low-frequency normal mode analysis of icosahedral viral capsids

The main goal of this work is the design of a coarse-grained theoretical model of minimal resolution for the study of the physical properties of icosahedral virus capsids within the linear-response regime. In this model the capsid is represented as an interacting many-body system whose composing elements are capsid subunits (capsomers), which are treated as three-dimensional rigid bodies. The total interaction potential energy is written as a sum of pairwise capsomer-capsomer interactions. Based on previous work [Gomez Llorente et al., Soft Matter, 2014, 10, 3560], a minimal and complete anisotropic binary interaction that includes a full Hessian matrix of independent force constants is proposed. In this interaction model, capsomers have rotational symmetry around an axis of order n > 2. The full coarse-grained model is applied to analyse the low-frequency normal-mode spectrum of icosahedral T = 1 capsids. The model performance is evaluated by fitting its predicted spectrum to the full-atom results for the Satellite Tobacco Necrosis Virus (STNV) capsid [Dykeman and Sankey, Phys. Rev. Lett., 2008, 100, 028101]. Two capsomer choices that are compatible with the capsid icosahedral symmetry are checked, namely pentamers (n = 5) and trimers (n = 3). Both subunit types provide fair fits, from which the magnitude of the coarse-grained force constants for a real virus is obtained. The model is able to uncover latent instabilities whose analysis is fully consistent with the current knowledge about the STNV capsid, which does not self-assemble in the absence of RNA and is thermally unstable. The straightforward generalisability of the model beyond the linear regime and its completeness make it a promising tool to theoretically interpret many experimental data such as those provided by the atomic force microscopy or even to better understand processes far from equilibrium such as the capsid self-assembly.

Kinetics of empty viral capsid assembly in a minimal model

The efficient construction of a protective protein shell or capsid is one of the most crucial steps in the replication cycle of a virus. The formation of the simplest capsid typically proceeds by the spontaneous assembly of identical building blocks. This process can also be achieved in vitro even in the absence of genetic material, thus opening the door to the production of artificial viral cages for a myriad of applications. In this work, we analyze the efficiency and the kinetic peculiarities of this self-assembly process using Brownian Dynamics simulations. We use a minimal model that considers identical assembly units and is able to reproduce successfully the correct final architecture of spherical capsids. The selection of a specific size and structure is achieved by changing a single parameter that imposes an angular anisotropy on the interaction. We analyze how the geometrical constraints of the interaction affect the efficiency of the assembly. We find that the optimal conditions for an efficient assembly from a kinetic point of view strongly depart from the lowest capsid energy corresponding to the minimum of the potential energy landscape. Our work illustrates the important differences between the equilibrium and dynamic characteristics of viral self-assembly, and provides important insights on how to design specific interactions for a successful assembly of artificial viral cages.

Supramolecular copolymers predominated by alternating order: Theory and application

We investigate the copolymerization behavior of a two-component system into quasilinear self-assemblies under conditions that interspecies binding is favored over identical species binding. The theoretical framework is based on a coarse-grained self-assembled Ising model with nearest neighbor interactions. In Ising language, such conditions correspond to the antiferromagnetic case giving rise to copolymers with predominantly alternating configurations. In the strong coupling limit, we show that the maximum fraction of polymerized material and the average length of strictly alternating copolymers depend on the stoichiometric ratio and the activation free energy of the more abundant species. They are substantially reduced when the stoichiometric ratio noticeably differs from unity. Moreover, for stoichiometric ratios close to unity, the copolymerization critical concentration is remarkably lower than the homopolymerization critical concentration of either species. We further analyze the polymerization behavior for a finite and negative coupling constant and characterize the composition of supramolecular copolymers. Our theoretical insights rationalize experimental results of supramolecular polymerization of oppositely charged monomeric species in aqueous solutions.

RNA-Mediated Virus Assembly: Mechanisms and Consequences for Viral Evolution and Therapy

Viruses, entities composed of nucleic acids, proteins, and in some cases lipids lack the ability to replicate outside their target cells. Their components self-assemble at the nanoscale with exquisite precision-a key to their biological success in infection. Recent advances in structure determination and the development of biophysical tools such as single-molecule spectroscopy and noncovalent mass spectrometry allow unprecedented access to the detailed assembly mechanisms of simple virions. Coupling these techniques with mathematical modeling and bioinformatics has uncovered a previously unsuspected role for genomic RNA in regulating formation of viral capsids, revealing multiple, dispersed RNA sequence/structure motifs [packaging signals (PSs)] that bind cognate coat proteins cooperatively. The PS ensemble controls assembly efficiency and accounts for the packaging specificity seen in vivo. The precise modes of action of the PSs vary between viral families, but this common principle applies across many viral families, including major human pathogens. These insights open up the opportunity to block or repurpose PS function in assembly for both novel antiviral therapy and gene/drug/vaccine applications.

Packaging of Genomic RNA in Positive-Sense Single-Stranded RNA Viruses: A Complex Story

The packaging of genomic RNA in positive-sense single-stranded RNA viruses is a key part of the viral infectious cycle, yet this step is not fully understood. Unlike double-stranded DNA and RNA viruses, this process is coupled with nucleocapsid assembly. The specificity of RNA packaging depends on multiple factors: (i) one or more packaging signals, (ii) RNA replication, (iii) translation, (iv) viral factories, and (v) the physical properties of the RNA. The relative contribution of each of these factors to packaging specificity is different for every virus. In vitro and in vivo data show that there are different packaging mechanisms that control selective packaging of the genomic RNA during nucleocapsid assembly. The goals of this article are to explain some of the key experiments that support the contribution of these factors to packaging selectivity and to draw a general scenario that could help us move towards a better understanding of this step of the viral infectious cycle.

Self-assembling outside equilibrium: emergence of structures mediated by dissipation

Self-assembly under non-equilibrium conditions may give rise to the formation of structures not available at equilibrium.

A dimorphism shift of hepatitis B virus capsids in response to ionic conditions

The dimorphism of HBV capsids (coexistence of T = 3 and T = 4 capsids) was found to be regulatable by controlling the rate of capsid nucleation using cations such as K+ or Ca2+: a quick addition of highly concentrated monovalent and/or multivalent counter-cations resulted in a morphism transition from a thermodynamically more stable, T = 4 capsid-dominant state (>80% of total capsids) to a new state containing ∼1 : 1 amounts of T = 3 and T = 4 capsids. These results suggested that the salts with strong charge screening ability could narrow the difference in nucleation energy barriers between the two states, which were not inter-convertible once formed. The effect of salts was more significant than other factors such as pH or protein concentration in achieving such a dimorphism shift. The general mechanism of HBV capsid dimorphism described here provides a new perspective in understanding the virus assembly during infection and directing the design of non-infectious capsids for nanotechnology applications.

Electrostatics-Driven Inflation of Elastic Icosahedral Shells as a Model for Swelling of Viruses

We develop a clear theoretical description of radial swelling in virus-like particles that delineates the importance of electrostatic contributions to swelling in the absence of any conformational changes. The model couples the elastic parameters of the capsid-represented as a continuous elastic shell-to the electrostatic pressure acting on it. We show that different modifications of the electrostatic interactions brought about by, for instance, changes in pH or solution ionic strength are often sufficient to achieve the experimentally observed swelling (∼10% of the capsid radius). Additionally, we derive analytical expressions for the electrostatics-driven radial swelling of virus-like particles that enable one to quickly estimate the magnitudes of physical quantities involved.

A modelling paradigm for RNA virus assembly

Virus assembly, a key stage in any viral life cycle, had long been considered to be primarily driven by protein-protein interactions and nonspecific interactions between genomic RNA and capsid protein. We review here a modelling paradigm for RNA virus assembly that illustrates the crucial roles of multiple dispersed, specific interactions between viral genomes and coat proteins in capsid assembly. The model reveals how multiple sequence-structure motifs in the genomic RNA, termed packaging signals, with a shared coat protein recognition motif enable viruses to overcome a viral assembly-equivalent of Levinthal's Paradox in protein folding. The fitness advantages conferred by this mechanism suggest that it should be widespread in viruses, opening up new perspectives on viral evolution and anti-viral therapy.

The role of the encapsulated cargo in microcompartment assembly

Bacterial microcompartments are large, roughly icosahedral shells that assemble around enzymes and reactants involved in certain metabolic pathways in bacteria. Motivated by microcompartment assembly, we use coarse-grained computational and theoretical modeling to study the factors that control the size and morphology of a protein shell assembling around hundreds to thousands of molecules. We perform dynamical simulations of shell assembly in the presence and absence of cargo over a range of interaction strengths, subunit and cargo stoichiometries, and the shell spontaneous curvature. Depending on these parameters, we find that the presence of a cargo can either increase or decrease the size of a shell relative to its intrinsic spontaneous curvature, as seen in recent experiments. These features are controlled by a balance of kinetic and thermodynamic effects, and the shell size is assembly pathway dependent. We discuss implications of these results for synthetic biology efforts to target new enzymes to microcompartment interiors.

Compactness of viral genomes: effect of disperse and localized random mutations

Genomes of single-stranded RNA viruses have evolved to optimize several concurrent properties. One of them is the architecture of their genomic folds, which must not only feature precise structural elements at specific positions, but also allow for overall spatial compactness. The latter was shown to be disrupted by random synonymous mutations, a disruption which can consequently negatively affect genome encapsidation. In this study, we use three mutation schemes with different degrees of locality to mutate the genomes of phage MS2 and Brome Mosaic virus in order to understand the observed sensitivity of the global compactness of their folds. We find that mutating local stretches of their genomes' sequence or structure is less disruptive to their compactness compared to inducing randomly-distributed mutations. Our findings are indicative of a mechanism for the conservation of compactness acting on a global scale of the genomes, and have several implications for understanding the interplay between local and global architecture of viral RNA genomes.

Packaging DNA Origami into Viral Protein Cages

The DNA origami technique is a widely used method to create customized, complex, spatially well-defined two-dimensional (2D) and three-dimensional (3D) DNA nanostructures. These structures have huge potential to serve as smart drug-delivery vehicles and molecular devices in various nanomedical and biotechnological applications. However, so far only little is known about the behavior of these novel structures in living organisms or in cell culture/tissue models. Moreover, enhancing pharmacokinetic bioavailability and transfection properties of such structures still remains a challenge. One intriguing approach to overcome these issues is to coat DNA origami nanostructures with proteins or lipid membranes. Here, we show how cowpea chlorotic mottle virus (CCMV) capsid proteins (CPs) can be used for coating DNA origami nanostructures. We present a method for disassembling native CCMV particles and isolating the pure CP dimers, which can further bind and encapsulate a rectangular DNA origami shape. Owing to the highly programmable nature of DNA origami, packaging of DNA nanostructures into viral protein cages could find imminent uses in enhanced targeting and cellular delivery of various active nano-objects, such as enzymes and drug molecules.

Self-assembly of amphiphilic truncated cones to form hollow nanovesicles

To mimic the unique properties of capsid (protein shell of a virus), we performed Brownian dynamics simulations of the self-assembly of amphiphilic truncated cone particles with anisotropic interactions. The particle shape of a truncated cone in our simulations depended on the cone angle θ, truncated height h_c and particle type (A_xB_y and B_xA_yB_z). The hydrophobic A moieties and hydrophilic B moieties are responsible for attractive and repulsive interactions, respectively. By varying the particle shape, truncated cones can assemble into hollow and vesicle-like clusters with a specific cluster size N. To assemble into hollow vesicles, the truncated height h_c must be below a critical value. When h_c exceeds this critical value, malformation will occur. The dynamics shows that the vesicle formation occurs in three stages: initially the growth is slow, then rapid, and finally it slows down. The truncated height h_c has a stronger impact on the growth kinetics than the cone angle θ or the particle type. We explored how the cluster packing depended on the cooling rate and particle number as well as discussing the relationship between the cluster geometry and the interparticle interactions. Further, we also discuss possible methods to experimentally prepare the truncated cones. The results of our work deepen our understanding of the self-assembly behavior of truncated cones and our results will aid the effective design of particle building blocks for novel nanostructures.

To mimic the unique properties of capsid (protein shell of a virus), we performed Brownian dynamics simulations of the self-assembly of amphiphilic truncated cone particles with anisotropic interactions.

Investigating the thermal dissociation of viral capsid by lattice model

The dissociation of icosahedral viral capsids was investigated by a homogeneous and a heterogeneous lattice model. In thermal dissociation experiments with cowpea chlorotic mottle virus and probed by small-angle neutron scattering, we observed a slight shrinkage of viral capsids, which can be related to the strengthening of the hydrophobic interaction between subunits at increasing temperature. By considering the temperature dependence of hydrophobic interaction in the homogeneous lattice model, we were able to give a better estimate of the effective charge. In the heterogeneous lattice model, two sets of lattice sites represented different capsid subunits with asymmetric interaction strengths. In that case, the dissociation of capsids was found to shift from a sharp one-step transition to a gradual two-step transition by weakening the hydrophobic interaction between AB and CC subunits. We anticipate that such lattice models will shed further light on the statistical mechanics underlying virus assembly and disassembly.

Ground state instabilities of protein shells are eliminated by buckling

We propose a hybrid discrete-continuum model to study the ground state of protein shells. The model allows for shape transformation of the shell and buckling transitions as well as the competition between states with different symmetries that characterize discrete particle models with radial pair potentials. Our main results are as follows. For large Föppl-von Kármán (FvK) numbers the shells have stable isometric ground states. As the FvK number is reduced, shells undergo a buckling transition resembling that of thin-shell elasticity theory. When the width of the pair potential is reduced below a critical value, then buckling coincides with the onset of structural instability triggered by over-stretched pair potentials. Chiral shells are found to be more prone to structural instability than achiral shells. It is argued that the well-width appropriate for protein shells lies below the structural instability threshold. This means that the self-assembly of protein shells with a well-defined, stable structure is possible only if the bending energy of the shell is sufficiently low so that the FvK number of the assembled shell is above the buckling threshold.

CREIM: Coffee Ring Effect Imaging Model for Monitoring Protein Self-Assembly in Situ

Protein self-assembly is fundamental to nanotechnology. Self-assembling structures are produced under static in vitro conditions typically forming over hours. In contrast, hydrodynamic intracellular environments employ far shorter time scales to compartmentalize highly concentrated protein solutions. Herein, we exploit the radial capillary flow within a drying sessile droplet (the coffee ring effect) to emulate dynamic native environments and monitor an archetypal protein assembly in situ using high-speed super-resolution imaging. We demonstrate that the assembly can be empirically driven to completion within minutes to seconds without apparent changes in supramolecular morphology. The model offers a reliable tool for the diagnosis and engineering of self-assembling systems under nonequilibrium conditions.

Shape Morphologies of Icosahedral Two-Component Vesicles

What are the features of partitioning of crystalline materials on the surface of a two-component icosahedral vesicle? We model the response of the rigid hardly stretchable crystalline icosahedra upon addition of a softer component on its surface. We demonstrate how the soft phase "invades" the shell regions with the highest elastic energy density around 12 5-fold topological defects. We explore the phase diagram of these inhomogeneous shells as a function of the soft material fraction, shell radius, and elastic moduli of the two phases. The findings are compared with the recent computer simulation findings, and their biological relevance, for example, for the structure of icosahedral viruses, is also discussed.

Beyond icosahedral symmetry in packings of proteins in spherical shells

The formation of quasi-spherical cages from protein building blocks is a remarkable self-assembly process in many natural systems, where a small number of elementary building blocks are assembled to build a highly symmetric icosahedral cage. In turn, this has inspired synthetic biologists to design de novo protein cages. We use simple models, on multiple scales, to investigate the self-assembly of a spherical cage, focusing on the regularity of the packing of protein-like objects on the surface. Using building blocks, which are able to pack with icosahedral symmetry, we examine how stable these highly symmetric structures are to perturbations that may arise from the interplay between flexibility of the interacting blocks and entropic effects. We find that, in the presence of those perturbations, icosahedral packing is not the most stable arrangement for a wide range of parameters; rather disordered structures are found to be the most stable. Our results suggest that (i) many designed, or even natural, protein cages may not be regular in the presence of those perturbations and (ii) optimizing those flexibilities can be a possible design strategy to obtain regular synthetic cages with full control over their surface properties.

Illuminating the Reaction Pathways of Viromimetic Assembly

The coassembly of well-defined biological nanostructures relies on a delicate balance between attractive and repulsive interactions between biomolecular building blocks. Viral capsids are a prototypical example, where coat proteins exhibit not only self-interactions but also interact with the cargo they encapsulate. In nature, the balance between antagonistic and synergistic interactions has evolved to avoid kinetic trapping and polymorphism. To date, it has remained a major challenge to experimentally disentangle the complex kinetic reaction pathways that underlie successful coassembly of biomolecular building blocks in a noninvasive approach with high temporal resolution. Here we show how macromolecular force sensors, acting as a genome proxy, allow us to probe the pathways through which a viromimetic protein forms capsids. We uncover the complex multistage process of capsid assembly, which involves recruitment and complexation, followed by allosteric growth of the proteinaceous coat. Under certain conditions, the single-genome particles condense into capsids containing multiple copies of the template. Finally, we derive a theoretical model that quantitatively describes the kinetics of recruitment and growth. These results shed new light on the origins of the pathway complexity in biomolecular coassembly.

Imaging and Quantitation of a Succession of Transient Intermediates Reveal the Reversible Self-Assembly Pathway of a Simple Icosahedral Virus Capsid

Understanding the fundamental principles underlying supramolecular self-assembly may facilitate many developments, from novel antivirals to self-organized nanodevices. Icosahedral virus particles constitute paradigms to study self-assembly using a combination of theory and experiment. Unfortunately, assembly pathways of the structurally simplest virus capsids, those more accessible to detailed theoretical studies, have been difficult to study experimentally. We have enabled the in vitro self-assembly under close to physiological conditions of one of the simplest virus particles known, the minute virus of mice (MVM) capsid, and experimentally analyzed its pathways of assembly and disassembly. A combination of electron microscopy and high-resolution atomic force microscopy was used to structurally characterize and quantify a succession of transient assembly and disassembly intermediates. The results provided an experiment-based model for the reversible self-assembly pathway of a most simple (T = 1) icosahedral protein shell. During assembly, trimeric capsid building blocks are sequentially added to the growing capsid, with pentamers of building blocks and incomplete capsids missing one building block as conspicuous intermediates. This study provided experimental verification of many features of self-assembly of a simple T = 1 capsid predicted by molecular dynamics simulations. It also demonstrated atomic force microscopy imaging and automated analysis, in combination with electron microscopy, as a powerful single-particle approach to characterize at high resolution and quantify transient intermediates during supramolecular self-assembly/disassembly reactions. Finally, the efficient in vitro self-assembly achieved for the oncotropic, cell nucleus-targeted MVM capsid may facilitate its development as a drug-encapsidating nanoparticle for anticancer targeted drug delivery.

Virus Matryoshka: A Bacteriophage Particle-Guided Molecular Assembly Approach to a Monodisperse Model of the Immature Human Immunodeficiency Virus

Immature human immunodeficiency virus type 1 (HIV-1) is approximately spherical, but is constructed from a hexagonal lattice of the Gag protein. As a hexagonal lattice is necessarily flat, the local symmetry cannot be maintained throughout the structure. This geometrical frustration presumably results in bending stress. In natural particles, the stress is relieved by incorporation of packing defects, but the magnitude of this stress and its significance for the particles is not known. In order to control this stress, we have now assembled the Gag protein on a quasi-spherical template derived from bacteriophage P22. This template is monodisperse in size and electron-transparent, enabling the use of cryo-electron microscopy in structural studies. These templated assemblies are far less polydisperse than any previously described virus-like particles (and, while constructed according to the same lattice as natural particles, contain almost no packing defects). This system gives us the ability to study the relationship between packing defects, curvature and elastic energy, and thermodynamic stability. As Gag is bound to the P22 template by single-stranded DNA, treatment of the particles with DNase enabled us to determine the intrinsic radius of curvature of a Gag lattice, unconstrained by DNA or a template. We found that this intrinsic radius is far larger than that of a virion or P22-templated particle. We conclude that Gag is under elastic strain in a particle; this has important implications for the kinetics of shell growth, the stability of the shell, and the type of defects it will assume as it grows.

Hidden symmetry of small spherical viruses and organization principles in "anomalous" and double-shelled capsid nanoassemblies

We propose the principles of structural organization in spherical nanoassemblies with icosahedral symmetry constituted by asymmetric protein molecules. The approach modifies the paradigmatic geometrical Caspar and Klug (CK) model of icosahedral viral capsids and demonstrates the common origin of both the "anomalous" and conventional capsid structures. In contrast to all previous models of "anomalous" viral capsids the proposed modified model conserves the basic structural principles of the CK approach and reveals the common hidden symmetry underlying all small viral shells. We demonstrate the common genesis of the "anomalous" and conventional capsids and explain their structures in the same frame. The organization principles are derived from the group theory analysis of the positional order on the spherical surface. The relationship between the modified CK geometrical model and the theory of two-dimensional spherical crystallization is discussed. We also apply the proposed approach to complex double-shelled capsids and capsids with protruding knob-like proteins. The introduced notion of commensurability for the concentric nanoshells explains the peculiarities of their organization and helps to predict analogous, but yet undiscovered, double-shelled viral capsid nanostructures.

Symmetry and Size of Membrane Protein Polyhedral Nanoparticles

In recent experiments [T. Basta et al., Proc. Natl. Acad. Sci. U.S.A. 111, 670 (2014)] lipids and membrane proteins were observed to self-assemble into membrane protein polyhedral nanoparticles (MPPNs) with a well-defined polyhedral protein arrangement and characteristic size. We develop a model of MPPN self-assembly in which the preferred symmetry and size of MPPNs emerge from the interplay of protein-induced lipid bilayer deformations, topological defects in protein packing, and thermal effects. With all model parameters determined directly from experiments, our model correctly predicts the observed symmetry and size of MPPNs. Our model suggests how key lipid and protein properties can be modified to produce a range of MPPN symmetries and sizes in experiments.

Kinetic Transition Networks for the Thomson Problem and Smale's Seventh Problem

The Thomson problem, arrangement of identical charges on the surface of a sphere, has found many applications in physics, chemistry and biology. Here, we show that the energy landscape of the Thomson problem for N particles with N=132, 135, 138, 141, 144, 147, and 150 is single funneled, characteristic of a structure-seeking organization where the global minimum is easily accessible. Algorithmically, constructing starting points close to the global minimum of such a potential with spherical constraints is one of Smale's 18 unsolved problems in mathematics for the 21st century because it is important in the solution of univariate and bivariate random polynomial equations. By analyzing the kinetic transition networks, we show that a randomly chosen minimum is, in fact, always "close" to the global minimum in terms of the number of transition states that separate them, a characteristic of small world networks.

Physical Ingredients Controlling Stability and Structural Selection of Empty Viral Capsids

One of the crucial steps in the viral replication cycle is the self-assembly of its protein shell. Typically, each native virus adopts a unique architecture, but the coat proteins of many viruses have the capability to self-assemble in vitro into different structures by changing the assembly conditions. However, the mechanisms determining which of the possible capsid shapes and structures is selected by a virus are still not well-known. We present a coarse-grained model to analyze and understand the physical mechanisms controlling the size and structure selection in the assembly of empty viral capsids. Using this model and Monte Carlo simulations, we have characterized the phase diagram and stability of T = 1,3,4,7 and snub cube shells. In addition, we have studied the tolerance of different shells to changes in physical parameters related to ambient conditions, identifying possible strategies to induce misassembly or failure. Finally, we discuss the factors that select the shape of a capsid as spherical, faceted, elongated, or decapsidated. Our model sheds important light on the ingredients that control the assembly and stability of viral shells. This knowledge is essential to get capsids with well-defined size and structure that could be used for promising applications in medicine or bionanotechnology.

Energetically favoured defects in dense packings of particles on spherical surfaces

The dense packing of interacting particles on spheres has proved to be a useful model for virus capsids and colloidosomes. Indeed, icosahedral symmetry observed in virus capsids corresponds to potential energy minima that occur for magic numbers of, e.g., 12, 32 and 72 identical Lennard-Jones particles, for which the packing has exactly the minimum number of twelve five-fold defects. It is unclear, however, how stable these structures are against thermal agitation. We investigate this property by means of basin-hopping global optimisation and Langevin dynamics for particle numbers between ten and one hundred. An important measure is the number and type of point defects, that is, particles that do not have six nearest neighbours. We find that small icosahedral structures are the most robust against thermal fluctuations, exhibiting fewer excess defects and rearrangements for a wide temperature range. Furthermore, we provide evidence that excess defects appearing at low non-zero temperatures lower the potential energy at the expense of entropy. At higher temperatures defects are, as expected, thermally excited and thus entropically stabilised. If we replace the Lennard-Jones potential by a very short-ranged (Morse) potential, which is arguably more appropriate for colloids and virus capsid proteins, we find that the same particle numbers give a minimum in the potential energy, although for larger particle numbers these minima correspond to different packings. Furthermore, defects are more difficult to excite thermally for the short-ranged potential, suggesting that the short-ranged interaction further stabilises equilibrium structures.

The Robust Assembly of Small Symmetric Nanoshells

Highly symmetric nanoshells are found in many biological systems, such as clathrin cages and viral shells. Many studies have shown that symmetric shells appear in nature as a result of the free-energy minimization of a generic interaction between their constituent subunits. We examine the physical basis for the formation of symmetric shells, and by using a minimal model, demonstrate that these structures can readily grow from the irreversible addition of identical subunits. Our model of nanoshell assembly shows that the spontaneous curvature regulates the size of the shell while the mechanical properties of the subunit determine the symmetry of the assembled structure. Understanding the minimum requirements for the formation of closed nanoshells is a necessary step toward engineering of nanocontainers, which will have far-reaching impact in both material science and medicine.