Elasticity theory and shape transitions of viral shells

Mechanical fatigue testing in silico: Dynamic evolution of material properties of nanoscale biological particles

Biological particles have evolved to possess mechanical characteristics necessary to carry out their functions. We developed a computational approach to "fatigue testing in silico", in which constant-amplitude cyclic loading is applied to a particle to explore its mechanobiology. We used this approach to describe dynamic evolution of nanomaterial properties and low-cycle fatigue in the thin spherical encapsulin shell, thick spherical Cowpea Chlorotic Mottle Virus (CCMV) capsid, and thick cylindrical microtubule (MT) fragment over 20 cycles of deformation. Changing structures and force-deformation curves enabled us to describe their damage-dependent biomechanics (strength, deformability, stiffness), thermodynamics (released and dissipated energies, enthalpy, and entropy) and material properties (toughness). Thick CCMV and MT particles experience material fatigue due to slow recovery and damage accumulation over 3-5 loading cycles; thin encapsulin shells show little fatigue due to rapid remodeling and limited damage. The results obtained challenge the existing paradigm: damage in biological particles is partially reversible owing to particle's partial recovery; fatigue crack may or may not grow with each loading cycle and may heal; and particles adapt to deformation amplitude and frequency to minimize the energy dissipated. Using crack size to quantitate damage is problematic as several cracks might form simultaneously in a particle. Dynamic evolution of strength, deformability, and stiffness, can be predicted by analyzing the cycle number (N) dependent damage, [Formula: see text] , where α is a power law and Nf is fatigue life. Fatigue testing in silico can now be used to explore damage-induced changes in the material properties of other biological particles. STATEMENT OF SIGNIFICANCE: Biological particles possess mechanical characteristics necessary to perform their functions. We developed "fatigue testing in silico" approach, which employes Langevin Dynamics simulations of constant-amplitude cyclic loading of nanoscale biological particles, to explore dynamic evolution of the mechanical, energetic, and material properties of the thin and thick spherical particles of encapsulin and Cowpea Chlorotic Mottle Virus, and the microtubule filament fragment. Our study of damage growth and fatigue development challenge the existing paradigm. Damage in biological particles is partially reversible as fatigue crack might heal with each loading cycle. Particles adapt to deformation amplitude and frequency to minimize energy dissipation. The evolution of strength, deformability, and stiffness, can be accurately predicted by analyzing the damage growth in particle structure.

Dynamics of fluid bilayer vesicles: Soft meshes and robust curvature energy discretization

Continuum models like the Helfrich Hamiltonian are widely used to describe fluid bilayer vesicles. Here we study the molecular dynamics compatible dynamics of the vertices of two-dimensional meshes representing the bilayer, whose in-plane motion is only weakly constrained. We show (i) that Jülicher's discretization of the curvature energy offers vastly superior robustness for soft meshes compared to the commonly employed expression by Gommper and Kroll and (ii) that for sufficiently soft meshes, the typical behavior of fluid bilayer vesicles can emerge even if the mesh connectivity remains fixed throughout the simulations. In particular, soft meshes can accommodate large shape transformations, and the model can generate the typical ℓ^{-4} signal for the amplitude of surface undulation modes of nearly spherical vesicles all the way up to the longest wavelength modes. Furthermore, we compare results for Newtonian, Langevin, and Brownian dynamics simulations of the mesh vertices to demonstrate that the internal friction of the membrane model is negligible, making it suitable for studying the internal dynamics of vesicles via coupling to hydrodynamic solvers or particle-based solvent models.

Effect of Mechanical Properties on Multicomponent Shell Patterning

Self-organized shells are fundamental in biological compartmentalization. They protect genomic material or enclose enzymes to aid the metabolic process. Studies of crystalline shells have shown the importance of the mechanical properties of building units in the shell morphology. However, the mechanism underlying the morphology of multicomponent assemblies is still poorly understood. Here, we analyze multicomponent closed shells that have different mechanical properties. By minimizing elastic energy, we show that heterogeneous bending rigidities regulate the surface pattern into circular, spikes, and ridge shapes. Interestingly, our continuum elasticity model recovers the patterns that have been proposed in bacterial microcompartments (BMCs), which are self-organized protein shells that aid the breakdown of complex molecules and allow bacteria to survive in hostile environments. In addition, our work elucidates the principles of pattern formation that can be used to design and engineer multicomponent microcompartments with a specific surface distribution of the components.

Fluctuating nonlinear spring theory: Strength, deformability, and toughness of biological nanoparticles from theoretical reconstruction of force-deformation spectra

We developed the Fluctuating Nonlinear Spring (FNS) model to describe the dynamics of mechanical deformation of biological particles, such as virus capsids. The theory interprets the force-deformation spectra in terms of the "Hertzian stiffness" (non-linear regime of a particle's small-amplitude deformations), elastic constant (large-amplitude elastic deformations), and force range in which the particle's fracture occurs. The FNS theory enables one to quantify the particles' elasticity (Young's moduli for Hertzian and bending deformations), and the limits of their strength (critical forces, fracture toughness) and deformability (critical deformations) as well as the probability distributions of these properties, and to calculate the free energy changes for the particle's Hertzian, elastic, and plastic deformations, and eventual fracture. We applied the FNS theory to describe the protein capsids of bacteriophage P22, Human Adenovirus, and Herpes Simplex virus characterized by deformations before fracture that did not exceed 10-19% of their size. These nanoshells are soft (~1-10-GPa elastic modulus), with low ~50-480-kPa toughness - a regime of material behavior that is not well understood, and with the strength increasing while toughness decreases with their size. The particles' fracture is stochastic, with the average values of critical forces, critical deformations, and fracture toughness comparable with their standard deviations. The FNS theory predicts 0.7-MJ/mol free energy for P22 capsid maturation, and it could be extended to describe uniaxial deformation of cylindrical microtubules and ellipsoidal cellular organelles.

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.

Irreversible and reversible morphological changes in the φ6 capsid and similar viral shells: symmetry and micromechanics

Understanding the physicochemical processes occurring in viruses during their maturation is of fundamental importance since only mature viruses can infect host cells. Here we consider the irreversible and reversible morphological changes that occur with the dodecahedral φ6 procapsid during the sequential packaging of 3 RNA segments forming the viral genome. It is shown that the dodecahedral shape of all the four observed capsid states is perfectly reproduced by a sphere radially deformed by only two irreducible spherical harmonics with icosahedral symmetry and wave numbers l = 6 and l = 10. The rotation of proteins around the 3-fold axes at the Procapsid → Intermediate 1 irreversible transformation is in fact also well described with the shear field containing only two irreducible harmonics with the same two wave numbers. The high stability of the Intermediate 1 state is discussed and the shapes of the Intermediate 2 state and Capsid (reversibly transforming back to the Intermediate 1 state) are shown to be mainly due to the isotropic pressure that the encapsidated RNA segments exert on the shell walls. The hidden symmetry of the capsid and the physicochemical features of the in vitro genome extraction from the viral shell are also elucidated.

Dislocation screening in crystals with spherical topology

Whereas disclination defects are energetically prohibitive in two-dimensional flat crystals, their existence is necessary in crystals with spherical topology, such as viral capsids, colloidosomes, or fullerenes. Such a geometrical frustration gives rise to large elastic stresses, which render the crystal unstable when its size is significantly larger than the typical lattice spacing. Depending on the compliance of the crystal with respect to stretching and bending deformations, these stresses are alleviated either by a local increase of the intrinsic curvature in proximity of the disclinations or by the proliferation of excess dislocations, often organized in the form of one-dimensional chains known as "scars." The associated strain field of the scars is such as to counterbalance the one resulting from the isolated disclinations. Here we develop a continuum theory of dislocation screening in two-dimensional closed crystals with genus one. Upon modeling the flux of scars emanating from a given disclination as an independent scalar field, we demonstrate that the elastic energy of closed two-dimensional crystals with various degrees of asphericity can be expressed as a simple quadratic function of the screened topological charge of the disclinations, at both zero and finite temperature. This allows us to predict the optimal density of the excess dislocations as well as the minimal stretching energy attained by the crystal.

Finite Temperature Phase Behavior of Viral Capsids as Oriented Particle Shells

A general phase plot is proposed for discrete particle shells that allows for thermal fluctuations of the shell geometry and of the inter-particle connectivities. The phase plot contains a first-order melting transition, a buckling transition, and a collapse transition and is used to interpret the thermodynamics of microbiological shells.

Icosadeltahedral Geometry of Geodesic Domes, Fullerenes and Viruses: A Tutorial on the T-Number

The Caspar–Klug (CK) classification of viruses is discussed by parallel examination of geometry of icosahedral geodesic domes, fullerenes, and viruses. The underlying symmetry of all structures is explained and thoroughly visually represented. Euler’s theorem on polyhedra is used to calculate the number of vertices, edges, and faces in domes, number of atoms, bonds, and pentagonal and hexagonal rings in fullerenes, and number of proteins and protein–protein contacts in viruses. The T-number, the characteristic for the CK classification, is defined and discussed. The superposition of fullerene and dome designs is used to obtain a representation of a CK virus with all the proteins indicated. Some modifications of the CK classifications are sketched, including elongation of the CK blueprint, fusion of two CK blueprints, dodecahedral view of the CK shapes, and generalized CK designs without a clearly visible geometry of the icosahedron. These are compared to cases of existing viruses.

Osmotic stress and pore nucleation in charged biological nanoshells and capsids

A model system is proposed to investigate the chemical equilibrium and mechanical stability of biological spherical-like nanoshells in contact with an aqueous solution with added dissociated electrolyte of a given concentration. The ionic chemical equilibrium across the permeable shell is investigated in the framework of an accurate Density Functional Theory (DFT) that incorporates electrostatic and hardcore correlations beyond the traditional mean-field (e.g., Poisson-Boltzmann) limit. The accuracy of the theory is tested by a direct comparison with Monte Carlo (MC) simulations. A simple analytical expression is then deduced which clearly highlights the entropic, electrostatic, and self-energy contributions to the osmotic stress over the shell in terms of the calculated ionic profiles. By invoking a continuum mean-field elastic approach to account for the shell surface stress upon osmotic stretching, the mechanical equilibrium properties of the shell under a wide variety of ionic strengths and surface charges are investigated. The model is further coupled to a continuum mechanical approach similar in structure to a Classical Nucleation Theory (CNT) to address the question of mechanical stability of the shells against a pore nucleation. This allows us to construct a phase diagram which delimits the mechanical stability of capsids for different ionic strengths and shell surface charges.

Actin shells control buckling and wrinkling of biomembranes

Global changes of cell shape under mechanical or osmotic external stresses are mostly controlled by the mechanics of the cortical actin cytoskeleton underlying the cell membrane. Some aspects of this process can be recapitulated in vitro on reconstituted actin-and-membrane systems. In this paper, we investigate how the mechanical properties of a branched actin network shell, polymerized at the surface of a liposome, control membrane shape when the volume is reduced. We observe a variety of membrane shapes depending on the actin thickness. Thin shells undergo buckling, characterized by a cup-shape deformation of the membrane that coincides with the one of the actin network. Thick shells produce membrane wrinkles, but do not deform their outer layer. For intermediate micrometer-thick shells, wrinkling of the membrane is observed, and the actin layer is slightly deformed. Confronting our experimental results with a theoretical description, we determine the transition between buckling and wrinkling, which depends on the thickness of the actin shell and the size of the liposome. We thus unveil the generic mechanism by which biomembranes are able to accommodate their shape against mechanical compression, through thickness adaptation of their cortical cytoskeleton.

pH-induced morphological changes of proteinaceous viral shells

Changes in environmental pH can induce morphological changes in empty proteinaceous shells of bacteriophages in vitro that are very similar to changes occurring in viral capsids in vivo after encapsidation of DNA. These changes in capsid shape and size cannot be explained with a simple elastic model alone. We propose a new theoretical framework that combines the elasticity of thin icosahedral shells with the pH dependence of capsid charge distribution. Minimization of the sum of elastic and electrostatic free energies leads to equilibrium shapes of viral shells that depend on a single elastic parameter and the detailed configuration of the imbedded protein charges. Based on the in vitro shell reconstructions of bacteriophage HK97 we elucidate the details of how the reversible transition between Prohead II and Expansion Intermediate II states of the HK97 procapsid is induced by pH changes, as well as some other features of the bacteriophage maturation.

Kirigami and the Caspar-Klug construction for viral shells with negative Gauss curvature

In this work we extend the Caspar-Klug construction to the archaeal viruses, which in recent years have captured the attention of many researchers for their ability to thrive in extreme environments. We assume that the shells of archaeal viruses are composed of hexamers and pentamers-as is true for icosahedral viruses-together with heptamers, necessary to introduce negative Gauss curvature. Following the original work of Caspar and Klug, we first construct models capable of reproducing the shape observed in electron microscopy images of archaeal viruses. Next, using the technique of kirigami, we present a systematic way to formulate archaeal virus templates from regular hexagonal lattices. Finally, we utilize the presented techniques to build finite element models of archaeal virus geometries and investigate their shapes as a function of material properties. In particular, using thin-shell elasticity theory, we describe a buckling transition as a function of a modified Föppl-von Kármán number γ^{★} and we show how changes in γ^{★} may initiate the tail formation in the Acidianus two-tailed archaeal virus.

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.

Shell buckling: from morphogenesis of soft matter to prospective applications

Being one of the commonest deformation modes for soft matter, shell buckling is the primary reason for the growth and nastic movement of many plants, as well as the formation of complex natural morphology. On-demand regulation of buckling-induced deformation associated with wrinkling, ruffling, folding, creasing and delaminating has profound implications for diverse scopes, which can be seen in its broad applications in microfabrication, 4D printing, actuator and drug delivery. This paper reviews the recent remarkable developments in the shell buckling of soft matter to explain the most representative natural morphogenesis from the perspectives of theoretical analysis in continuum mechanics, finite element analysis, and experimental validations. Imitation of buckling-induced shape transformation and its applications are also discussed for the innovations of sophisticated materials and devices in future.

Viral nanomechanics with a virtual atomic force microscope

One of the most important components of a virus is the protein shell or capsid that encloses its genetic material. The main role of the capsid is to protect the viral genome against external aggressions, facilitating its safe and efficient encapsulation and delivery. As a consequence, viral capsids have developed astonishing mechanical properties that are crucial for viral function. These remarkable properties have started to be unveiled in single-virus nanoindentation experiments, and are opening the door to the use of viral-derived artificial nanocages for promising bio- and nano-technological applications. However, the interpretation of nanoindentation experiments is often difficult, requiring the support of theoretical and simulation analysis. Here we present a 'Virtual AFM' (VAFM), a Brownian Dynamics simulation of a coarse-grained model of virus aimed to mimic the standard setup of atomic force microscopy (AFM) nanoindentation experiments. Despite the heavy level of coarse-graining, these simulations provide valuable information which is not accessible in experiments. Rather than focusing on a specific virus, the VAFM will be used to analyze how the mechanical response and breaking of viruses depend on different parameters controlling the effective interactions between capsid's structural units. In particular, we will discuss the influence of adsorption, the tip radius, and the rigidity and shape of the shell on its mechanical response.

Extension of Caspar-Klug theory to higher order pentagonal polyhedra

Many viral capsids follow an icosahedral fullerene-like structure, creating a caged polyhedral arrangement built entirely from hexagons and pentagons. Viral capsids consist of capsid proteins,which group into clusters of six (hexamers) or five (pentamers). Although the number of hexamers per capsid varies depending on the capsid size, Caspar-Klug Theory dictates there are exactly twelve pentamers needed to form a closed capsid.However, for a significant number of viruses, including viruses of the Papovaviridae family, the theory doesn’t apply. The anomaly of the Caspar-Klug Theory has raised a new question:“For which Caspar and Klug models can each hexamer be replaced with a pentamer while still following icosahedral symmetry?” This paper proposes an answer to this question by examining icosahedral viral capsid-like structures composed only of pentamers, called pentagonal polyhedra. The analysis shows that pentagonal polyhedra fall in a subclass of T, defined by P ≥ 7 and T = 1( mod 3).

Determination of prestress and elastic properties of virus capsids

Virus capsids are protein shells that protect the virus genome, and determination of their mechanical properties has been a topic of interest because of their potential use in nanotechnology and therapeutics. It has been demonstrated that stresses exist in virus capsids, even in their equilibrium state, due to their construction. These stresses, termed "prestresses" in this study, closely affect the capsid's mechanical behavior. Three methods-shape-based metric, atomic force microscope indentation, and molecular dynamics-have been proposed to determine the capsid elastic properties without fully accounting for prestresses. In this paper, we theoretically analyze the three methods used for mechanical characterization of virus capsids and numerically investigate how prestresses affect the capsid's mechanical properties. We consolidate all the results and propose that by using these techniques collectively, it is possible to accurately determine both the mechanical properties and prestresses in capsids.

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.

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.

Orientational phase transitions and the assembly of viral capsids

We present a Landau theory for large-l orientational phase transitions and apply it to the assembly of icosahedral viral capsids. The theory predicts two distinct types of ordering transitions. Transitions dominated by the l=6,10,12, and 18 icosahedral spherical harmonics resemble robust first-order phase transitions that are not significantly affected by chirality. The remaining transitions depend essentially on including mixed l states denoted as l=15+16 corresponding to a mixture of l=15 and l=16 spherical harmonics. The l=15+16 transition is either continuous or weakly first-order and it is strongly influenced by chirality, which suppresses spontaneous chiral symmetry breaking. The icosahedral state is in close competition with states that have tetrahedral, D_{5}, and octahedral symmetries. We present a group-theoretic method to analyze the competition between the different symmetries. The theory is applied to a variety of viral shells.

Viral self-assembly pathway and mechanical stress relaxation

The final shape of a virus is dictated by the self-assembly pathway of its constituents. Using standard thin-shell elasticity, we highlight the prominent role of the viral shell's spontaneous curvature in determining the assembly pathway. In particular, we demonstrate that the mechanical stress inherent to the growth of a curved surface can be relaxed in two different ways in the early steps of assembly, depending on the value of the spontaneous curvature of the surface. This important result explains why most viral shells have either a compact shape with icosahedral symmetry or an elongated shape lacking this symmetry.

In vitro protease cleavage and computer simulations reveal the HIV-1 capsid maturation pathway

HIV-1 virions assemble as immature particles containing Gag polyproteins that are processed by the viral protease into individual components, resulting in the formation of mature infectious particles. There are two competing models for the process of forming the mature HIV-1 core: the disassembly and de novo reassembly model and the non-diffusional displacive model. To study the maturation pathway, we simulate HIV-1 maturation in vitro by digesting immature particles and assembled virus-like particles with recombinant HIV-1 protease and monitor the process with biochemical assays and cryoEM structural analysis in parallel. Processing of Gag in vitro is accurate and efficient and results in both soluble capsid protein and conical or tubular capsid assemblies, seemingly converted from immature Gag particles. Computer simulations further reveal probable assembly pathways of HIV-1 capsid formation. Combining the experimental data and computer simulations, our results suggest a sequential combination of both displacive and disassembly/reassembly processes for HIV-1 maturation.

Two competing models—disassembly/reassembly and displacive—have been proposed for how immature spherical HIV virions transform into mature particles with conical cores. Here the authors provide evidence that both disassembly/reassembly and displacive processes occur sequentially during the maturation process.

Reaction-diffusion basis of retroviral infectivity

Retrovirus particle (virion) infectivity requires diffusion and clustering of multiple transmembrane envelope proteins (Env₃) on the virion exterior, yet is triggered by protease-dependent degradation of a partially occluding, membrane-bound Gag polyprotein lattice on the virion interior. The physical mechanism underlying such coupling is unclear and only indirectly accessible via experiment. Modelling stands to provide insight but the required spatio-temporal range far exceeds current accessibility by all-atom or even coarse-grained molecular dynamics simulations. Nor do such approaches account for chemical reactions, while conversely, reaction kinetics approaches handle neither diffusion nor clustering. Here, a recently developed multiscale approach is considered that applies an ultra-coarse-graining scheme to treat entire proteins at near-single particle resolution, but which also couples chemical reactions with diffusion and interactions. A model is developed of Env₃ molecules embedded in a truncated Gag lattice composed of membrane-bound matrix proteins linked to capsid subunits, with freely diffusing protease molecules. Simulations suggest that in the presence of Gag but in the absence of lateral lattice-forming interactions, Env₃ diffuses comparably to Gag-absent Env₃ Initial immobility of Env₃ is conferred through lateral caging by matrix trimers vertically coupled to the underlying hexameric capsid layer. Gag cleavage by protease vertically decouples the matrix and capsid layers, induces both matrix and Env₃ diffusion, and permits Env₃ clustering. Spreading across the entire membrane surface reduces crowding, in turn, enhancing the effect and promoting infectivity.This article is part of the themed issue 'Multiscale modelling at the physics-chemistry-biology interface'.

Assembly of the most topologically regular two-dimensional micro and nanocrystals with spherical, conical, and tubular shapes

We consider how to control the extension of curvature-induced defects in the hexagonal order covering different curved surfaces. In these frames we propose a physical mechanism for improving structures of two-dimensional spherical colloidal crystals (SCCs). For any SCC comprising of about 300 or less particles the mechanism transforms all extended topological defects (ETDs) in the hexagonal order into the point disclinations. Perfecting the structure is carried out by successive cycles of the particle implantation and subsequent relaxation of the crystal. The mechanism is potentially suitable for obtaining colloidosomes with better selective permeability. Our approach enables modeling the most topologically regular tubular and conical two-dimensional nanocrystals including various possible polymorphic forms of the HIV viral capsid. Different HIV-like shells with an arbitrary number of structural units (SUs) and desired geometrical parameters are easily formed. Faceting of the obtained structures is performed by minimizing the suggested elastic energy.

Imperfection sensitivity of pressured buckling of biopolymer spherical shells

Imperfection sensitivity is essential for mechanical behavior of biopolymer shells [such as ultrasound contrast agents (UCAs) and spherical viruses] characterized by high geometric heterogeneity. In this work, an imperfection sensitivity analysis is conducted based on a refined shell model recently developed for spherical biopolymer shells of high structural heterogeneity and thickness nonuniformity. The influence of related parameters (including the ratio of radius to average shell thickness, the ratio of transverse shear modulus to in-plane shear modulus, and the ratio of effective bending thickness to average shell thickness) on imperfection sensitivity is examined for pressured buckling. Our results show that the ratio of effective bending thickness to average shell thickness has a major effect on the imperfection sensitivity, while the effect of the ratio of transverse shear modulus to in-plane shear modulus is usually negligible. For example, with physically realistic parameters for typical imperfect spherical biopolymer shells, the present model predicts that actual maximum external pressure could be reduced to as low as 60% of that of a perfect UCA spherical shell or 55%-65% of that of a perfect spherical virus shell, respectively. The moderate imperfection sensitivity of spherical biopolymer shells with physically realistic imperfection is largely attributed to the fact that biopolymer shells are relatively thicker (defined by smaller radius-to-thickness ratio) and therefore practically realistic imperfection amplitude normalized by thickness is very small as compared to that of classical elastic thin shells which have much larger radius-to-thickness ratio.

Useful scars: Physics of the capsids of archaeal viruses

We propose a physical model for the capsids of tailed archaeal viruses as viscoelastic membranes under tension. The fluidity is generated by thermal motion of scarlike structures that are an intrinsic feature of the ground state of large particle arrays covering surfaces with nonzero Gauss curvature. The tension is generated by a combination of the osmotic pressure of the enclosed genome and an extension force generated by filamentous structure formation that drives the formation of the tails. In continuum theory, the capsid has the shape of a surface of constant mean curvature: an unduloid. Particle arrays covering unduloids are shown to exhibit pronounced subdiffusive and diffusive single-particle transport at temperatures that are well below the melting temperature of defect-free particle arrays on a surface with zero Gauss curvature.

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.

Toward Virus-Like Surface Plasmon Strain Sensors

The strong configuration dependence of collective surface plasmon resonances in an array of metal nanoparticles provides an opportunity to develop a bioinspired tool for sensing mechanical deformations in soft matter at the nanoscale. We study the feasibility of a strain sensor based on an icosahedral array of nanoparticles encapsulated by a virus capsid. When the system undergoes deformation, the optical scattering cross-section spectra as well as the induced electric field profile change. By numerical simulations, we examine how these changes depend on the symmetry and extent of the deformation and on both the propagation direction and polarization of the incident radiation. Such a sensor could prove useful in studies of the mechanisms of nanoparticle or virus translocation in the confines of a host cell.

Role of Genome in the Formation of Conical Retroviral Shells

Human immunodeficiency virus (HIV) capsid proteins spontaneously assemble around the genome into a protective protein shell called the capsid, which can take on a variety of shapes broadly classified as conical, cylindrical, and irregular. The majority of capsids seen in in vivo studies are conical in shape, while in vitro experiments have shown a preference for cylindrical capsids. The factors involved in the selection of the unique shape of HIV capsids are not well understood, and in particular the impact of RNA on the formation of the capsid is not known. In this work, we study the role of the genome and its interaction with the capsid protein by modeling the genomic RNA through a mean-field theory. Our results show that the confinement free energy for a homopolymeric model genome confined in a conical capsid is lower than that in a cylindrical capsid, at least when the genome does not interact with the capsid, which seems to be the case in in vivo experiments. Conversely, the confinement free energy for the cylinder is lower than that for a conical capsid if the genome is attracted to the capsid proteins as the in vitro experiments. Understanding the factors that contribute to the formation of conical capsids may shed light on the infectivity of HIV particles.

Effects of long-range interactions on curvature energies of viral shells

We formulate a theory of the effects of long-range interactions on the surface tension and spontaneous curvature of proteinaceous shells based on the general Deryaguin-Landau-Verwey-Overbeek mesoscale approach to colloid stability. We derive the full renormalization formulas for the elastic properties of the shell and consider in detail the renormalization of the spontaneous curvature as a function of the corresponding Hamaker coefficient, inner and outer capsid charges, and bathing solution properties. The renormalized spontaneous curvature is found to be a nonmonotonic function of several parameters describing the system.

HIV Capsid Assembly, Mechanism, and Structure

The HIV genome materials are encaged by a proteinaceous shell called the capsid, constructed from ∼1000-1500 copies of the capsid proteins. Because its stability and integrity are critical to the normal life cycle and infectivity of the virus, the HIV capsid is a promising antiviral drug target. In this paper, we review the studies shaping our understanding of the structure and dynamics of the capsid proteins and various forms of their assemblies, as well as the assembly mechanism.

Influence of material stretchability on the equilibrium shape of a Möbius band

We use a two-dimensional discrete, lattice-based model to show that Möbius bands made with stretchable materials are less likely to crease or tear. This stems from a delocalization of twisting strain that occurs if stretching is allowed. The associated low-energy configurations provide strategic target shapes for the guided assembly of nanometer and micron scale Möbius bands. To predict macroscopic band shapes for a given material, we establish a connection between stretchability and relevant continuum moduli, leading to insight regarding the practical feasibility of synthesizing Möbius bands from materials with continuum parameters that can be measured experimentally or estimated by upscale averaging.

Kinetics versus Thermodynamics in Virus Capsid Polymorphism

Virus coat proteins spontaneously self-assemble into empty shells in aqueous solution under the appropriate physicochemical conditions, driven by an interaction free energy per bond on the order of 2-5 times the thermal energy kBT. For this seemingly modest interaction strength, each protein building block nonetheless gains a very large binding free energy, between 10 and 20 kBT. Because of this, there is debate about whether the assembly process is reversible or irreversible. Here we discuss capsid polymorphism observed in in vitro experiments from the perspective of nucleation theory and of the thermodynamics of mass action. We specifically consider the potential contribution of a curvature free energy term to the effective interaction potential between the proteins. From these models, we propose experiments that may conclusively reveal whether virus capsid assembly into a mixture of polymorphs is a reversible or an irreversible process.

Mechanism of polymorphism and curvature of HIV capsid assemblies probed by 3D simulations with a novel coarse grain model

BACKGROUND

During the maturation process, HIV capsid proteins self-assemble into polymorphic capsids. The strong polymorphism precludes high resolution structural characterization under in vivo conditions. In spite of the determination of structural models for various in vitro assemblies of HIV capsid proteins, the assembly mechanism is still not well-understood.

METHODS

We report 3D simulations of HIV capsid proteins by a novel coarse grain model that captures the backbone of the rigid segments in the protein accurately. The effects of protein dynamics on assembly are emulated by a static ensemble of subunits in conformations derived from molecular dynamics simulation.

RESULTS

We show that HIV capsid proteins robustly assemble into hexameric lattices in a range of conditions where trimers of dimeric subunits are the dominant oligomeric intermediates. Variations of hexameric lattice curvatures are observed in simulations with subunits of variable inter-domain orientations mimicking the conformation distribution in solution. Simulations with subunits based on pentameric structural models lead to assembly of sharp curved structures resembling the tips of authentic HIV capsids, along a distinct pathway populated by tetramers and pentamers with the characteristic quasi-equivalency of viral capsids.

CONCLUSIONS

Our results suggest that the polymorphism assembly is triggered by the inter-domain dynamics of HIV capsid proteins in solution. The assembly of highly curved structures arises from proteins in conformation with a highly specific inter-domain orientation.

SIGNIFICANCE

Our work proposes a mechanism of HIV capsid assembly based on available structural data, which can be readily verified. Our model can be applied to other large biomolecular assemblies.

Modeling Viral Capsid Assembly

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

DNA buckling in bacteriophage cavities as a mechanism to aid virus assembly

While relatively simple biologically, bacteriophages are sophisticated biochemical machines that execute a precise sequence of events during virus assembly, DNA packaging, and ejection. These stages of the viral life cycle require intricate coordination of viral components whose structures are being revealed by single molecule experiments and high resolution (cryo-electron microscopy) reconstructions. For example, during packaging, bacteriophages employ some of the strongest known molecular motors to package DNA against increasing pressure within the viral capsid shell. Located upstream of the motor is an elaborate portal system through which DNA is threaded. A high resolution reconstruction of the portal system for bacteriophage ϕ29 reveals that DNA buckles inside a small cavity under large compressive forces. In this study, we demonstrate that DNA can also buckle in other bacteriophages including T7 and P22. Using a computational rod model for DNA, we demonstrate that a DNA buckle can initiate and grow within the small confines of a cavity under biologically-attainable force levels. The forces of DNA-cavity contact and DNA-DNA electrostatic repulsion ultimately limit cavity filling. Despite conforming to very different cavity geometries, the buckled DNA within T7 and P22 exhibits near equal volumetric energy density (∼1kT/nm(3)) and energetic cost of packaging (∼22kT). We hypothesize that a DNA buckle creates large forces on the cavity interior to signal the conformational changes to end packaging. In addition, a DNA buckle may help retain the genome prior to tail assembly through significantly increased contact area with the portal.

Nonuniform growth and topological defects in the shaping of elastic sheets

We demonstrate that shapes with zero Gaussian curvature, except at singularities, produced by the growth-induced buckling of a thin elastic sheet are the same as those produced by the Volterra construction of topological defects in which edges of an intrinsically flat surface are identified. With this connection, we study the problem of choosing an optimal pattern of growth for a prescribed developable surface, finding a fundamental trade-off between optimal design and the accuracy of the resulting shape which can be quantified by the length along which an edge should be identified.

Microcapsule mechanics: from stability to function

Microcapsules are reviewed with special emphasis on the relevance of controlled mechanical properties for functional aspects. At first, assembly strategies are presented that allow control over the decisive geometrical parameters, diameter and wall thickness, which both influence the capsule's mechanical performance. As one of the most powerful approaches the layer-by-layer technique is identified. Subsequently, ensemble and, in particular, single-capsule deformation techniques are discussed. The latter generally provide more in-depth information and cover the complete range of applicable forces from smaller than pN to N. In a theory chapter, we illustrate the physics of capsule deformation. The main focus is on thin shell theory, which provides a useful approximation for many deformation scenarios. Finally, we give an overview of applications and future perspectives where the specific design of mechanical properties turns microcapsules into (multi-)functional devices, enriching especially life sciences and material sciences.

Recent Developments in Molecular Simulation Approaches to Study Spherical Virus Capsids

Viruses are a particularly challenging systems to study via molecular simulation methods. Virus capsids typically consist of over 100 subunit proteins and reach dimensions of over 100 nm; solvated viruses capsid systems can be over 1 million atoms in size. In this review, I will present recent developments which have attempted to overcome the significant computational expense to perform simulations which can inform experimental studies, make useful predictions about biological phenomena and calculate material properties relevant to nanotechnology design efforts.

Impact of the topology of viral RNAs on their encapsulation by virus coat proteins

Single-stranded RNAs of simple viruses seem to be topologically more compact than other types of single-stranded RNA. It has been suggested that this has an evolutionary purpose: more compact structures are more easily encapsulated in the limited space that the cavity of the virus capsid offers. We employ a simple Flory theory to calculate the optimal amount of polymers confined in a viral shell. We find that the free energy gain or more specifically the efficiency of RNA encapsidation increases substantially with topological compactness. We also find that the optimal length of RNA encapsidated in a capsid increases with the degree of branching of the genome even though this effect is very weak. Further, we show that if the structure of the branching of the polymer is allowed to anneal, the optimal loading increases substantially.

Modeling and simulation of the mechanical response from nanoindentation test of DNA-filled viral capsids

Viruses can be described as biological objects composed mainly of two parts: a stiff protein shell called a capsid, and a core inside the capsid containing the nucleic acid and liquid. In many double-stranded DNA bacterial viruses (aka phage), the volume ratio between the liquid and the encapsidated DNA is approximately 1:1. Due to the dominant DNA hydration force, water strongly mediates the interaction between the packaged DNA strands. Therefore, water that hydrates the DNA plays an important role in nanoindentation experiments of DNA-filled viral capsids. Nanoindentation measurements allow us to gain further insight into the nature of the hydration and electrostatic interactions between the DNA strands. With this motivation, a continuum-based numerical model for simulating the nanoindentation response of DNA-filled viral capsids is proposed here. The viral capsid is modeled as large- strain isotropic hyper-elastic material, whereas porous elasticity is adopted to capture the mechanical response of the filled viral capsid. The voids inside the viral capsid are assumed to be filled with liquid, which is modeled as a homogenous incompressible fluid. The motion of a fluid flowing through the porous medium upon capsid indentation is modeled using Darcy's law, describing the flow of fluid through a porous medium. The nanoindentation response is simulated using three-dimensional finite element analysis and the simulations are performed using the finite element code Abaqus. Force-indentation curves for empty, partially and completely DNA-filled capsids are directly compared to the experimental data for bacteriophage λ. Material parameters such as Young's modulus, shear modulus, and bulk modulus are determined by comparing computed force-indentation curves to the data from the atomic force microscopy (AFM) experiments. Predictions are made for pressure distribution inside the capsid, as well as the fluid volume ratio variation during the indentation test.

Theoretical studies on assembly, physical stability and dynamics of viruses

All matter has to obey the general laws of physics and living matter is not an exception. Viruses have not only learnt how to cope with them, but have managed to use them for their own survival. In this chapter we will review some of the exciting physics behind viruses and discuss simple physical models that can shed some light on different aspects of the viral life cycle and viral properties. In particular, we will focus on how the structure and shape of the capsid, its assembly and stability, and the entry and exit of viral particles and their genomes can be understood using fundamental physics theories.

Fluctuating shells under pressure

Thermal fluctuations strongly modify the large length-scale elastic behavior of cross-linked membranes, giving rise to scale-dependent elastic moduli. Whereas thermal effects in flat membranes are well understood, many natural and artificial microstructures are modeled as thin elastic shells. Shells are distinguished from flat membranes by their nonzero curvature, which provides a size-dependent coupling between the in-plane stretching modes and the out-of-plane undulations. In addition, a shell can support a pressure difference between its interior and its exterior. Little is known about the effect of thermal fluctuations on the elastic properties of shells. Here, we study the statistical mechanics of shape fluctuations in a pressurized spherical shell, using perturbation theory and Monte Carlo computer simulations, explicitly including the effects of curvature and an inward pressure. We predict novel properties of fluctuating thin shells under point indentations and pressure-induced deformations. The contribution due to thermal fluctuations increases with increasing ratio of shell radius to thickness and dominates the response when the product of this ratio and the thermal energy becomes large compared with the bending rigidity of the shell. Thermal effects are enhanced when a large uniform inward pressure acts on the shell and diverge as this pressure approaches the classical buckling transition of the shell. Our results are relevant for the elasticity and osmotic collapse of microcapsules.