Continuum theory of retroviral capsids

Hexagonal Lattices of HIV Capsid Proteins Explored by Simulations Based on a Thermodynamically Consistent Model

HIV capsid proteins (CAs) may self-assemble into a variety of shapes under in vivo and in vitro conditions. Here, we employed simulations based on a residue-level coarse-grained (CG) model with full conformational flexibility to investigate hexagonal lattices, which are the underlying structural pattern for CA aggregations. Facilitated by enhanced sampling simulations to rigorously calculate CA dimerization and polymerization affinities, we calibrated our model to reproduce the experimentally measured affinities. Using the calibrated model, we performed unbiased simulations on several large systems consisting of 1512 CA subunits, allowing reversible binding and unbinding of the CAs in a thermodynamically consistent manner. In one simulation, a preassembled hexagonal CA sheet developed spontaneous curvatures reminiscent of those observed in experiments, and the edges of the sheet exhibited local curvatures larger than those of the interior. In other simulations starting with randomly distributed CAs at different concentrations, existing CA assemblies grew by binding free capsomeres to the edges and by merging with other assemblies. At high CA concentrations, rapid establishment of predominant aggregates was followed by much slower adjustments toward more regular hexagonal lattices, with increasing numbers of intact CA hexamers and pentamers being formed. Our approach of adapting a general CG model to specific systems by using experimental binding data represents a practical and effective strategy for simulating and elucidating intricate protein aggregations.

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.

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.

The equilibrium structure of self-assembled protein nano-cages

Understanding how highly symmetric, robust, monodisperse protein nano-cages self-assemble can have major applications in various areas of bio-nanotechnology, such as drug delivery, biomedical imaging and gene therapy. We develop a model to investigate the assembly of protein subunits into the structures with different sizes and symmetries. Using Monte Carlo simulation, we obtain global minimum energy structures. Our results suggest that the physical properties including the spontaneous curvature, flexibility and bending rigidity of coat proteins are sufficient to predict the size, symmetry and shape selectivity of the assembly products. Further, on a thermodynamic basis, we discuss the polymorphism of nano-cages observed in assembly experiments.

Why large icosahedral viruses need scaffolding proteins

While small single-stranded viral shells encapsidate their genome spontaneously, many large viruses, such as the herpes simplex virus or infectious bursal disease virus (IBDV), typically require a template, consisting of either scaffolding proteins or an inner core. Despite the proliferation of large viruses in nature, the mechanisms by which hundreds or thousands of proteins assemble to form structures with icosahedral order (IO) is completely unknown. Using continuum elasticity theory, we study the growth of large viral shells (capsids) and show that a nonspecific template not only selects the radius of the capsid, but also leads to the error-free assembly of protein subunits into capsids with universal IO. We prove that as a spherical cap grows, there is a deep potential well at the locations of disclinations that later in the assembly process will become the vertices of an icosahedron. Furthermore, we introduce a minimal model and simulate the assembly of a viral shell around a template under nonequilibrium conditions and find a perfect match between the results of continuum elasticity theory and the numerical simulations. Besides explaining available experimental results, we provide a number of predictions. Implications for other problems in spherical crystals are also discussed.

Defects and Chirality in the Nanoparticle-Directed Assembly of Spherocylindrical Shells of Virus Coat Proteins

Virus coat proteins of small isometric plant viruses readily assemble into symmetric, icosahedral cages encapsulating noncognate cargo, provided the cargo meets a minimal set of chemical and physical requirements. While this capability has been intensely explored for certain virus-enabled nanotechnologies, additional applications require lower symmetry than that of an icosahedron. Here, we show that the coat proteins of an icosahedral virus can efficiently assemble around metal nanorods into spherocylindrical closed shells with hexagonally close-packed bodies and icosahedral caps. Comparison of chiral angles and packing defects observed by in situ atomic force microscopy with those obtained from molecular dynamics models offers insight into the mechanism of growth, and the influence of stresses associated with intrinsic curvature and assembly pathways.

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.

Physical properties of the HIV-1 capsid from all-atom molecular dynamics simulations

Human immunodeficiency virus type 1 (HIV-1) infection is highly dependent on its capsid. The capsid is a large container, made of ∼1,300 proteins with altogether 4 million atoms. Although the capsid proteins are all identical, they nevertheless arrange themselves into a largely asymmetric structure made of hexamers and pentamers. The large number of degrees of freedom and lack of symmetry pose a challenge to studying the chemical details of the HIV capsid. Simulations of over 64 million atoms for over 1 μs allow us to conduct a comprehensive study of the chemical–physical properties of an empty HIV-1 capsid, including its electrostatics, vibrational and acoustic properties, and the effects of solvent (ions and water) on the capsid. The simulations reveal critical details about the capsid with implications to biological function.

The large and complex HIV-1 capsid modulates several molecular events during HIV-1’s infective cycle. Here the authors use all-atom molecular dynamic simulations to probe the biophysical properties of the genome-free HIV-1 capsid.

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.

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.

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.

Molecular Architecture of the Retroviral Capsid

Retroviral capsid cores are proteinaceous containers that self-assemble to encase the viral genome and a handful of proteins that promote infection. Their function is to protect and aid in the delivery of viral genes to the nucleus of the host, and, in many cases, infection pathways are influenced by capsid-cellular interactions. From a mathematical perspective, capsid cores are polyhedral cages and, as such, follow well-defined geometric rules. However, marked morphological differences in shapes exist, depending on virus type. Given the specific roles of capsid in the viral life cycle, the availability of detailed molecular structures, particularly at assembly interfaces, opens novel avenues for targeted drug development against these pathogens. Here, we summarize recent advances in the structure and understanding of retroviral capsid, with particular emphasis on assemblies and the capsid cores.

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.

Applying molecular crowding models to simulations of virus capsid assembly in vitro

Virus capsid assembly has been widely studied as a biophysical system, both for its biological and medical significance and as an important model for complex self-assembly processes. No current technology can monitor assembly in detail and what information we have on assembly kinetics comes exclusively from in vitro studies. There are many differences between the intracellular environment and that of an in vitro assembly assay, however, that might be expected to alter assembly pathways. Here, we explore one specific feature characteristic of the intracellular environment and known to have large effects on macromolecular assembly processes: molecular crowding. We combine prior particle simulation methods for estimating crowding effects with coarse-grained stochastic models of capsid assembly, using the crowding models to adjust kinetics of capsid simulations to examine possible effects of crowding on assembly pathways. Simulations suggest a striking difference depending on whether or not a system uses nucleation-limited assembly, with crowding tending to promote off-pathway growth in a nonnucleation-limited model but often enhancing assembly efficiency at high crowding levels even while impeding it at lower crowding levels in a nucleation-limited model. These models may help us understand how complicated assembly systems may have evolved to function with high efficiency and fidelity in the densely crowded environment of the cell.

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.

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.

Determinants of the HIV-1 core assembly pathway

Based on structural information, we have analyzed the mechanism of mature HIV-1 core assembly and the contributions of structural elements to the assembly process. Through the use of several in vitro assembly assay systems, we have examined details of how capsid (CA) protein helix 1, ß-hairpin and cyclophilin loop elements impact assembly-dependent protein interactions, and we present evidence for a contribution of CA helix 6 to the mature assembly-competent conformation of CA. Additional experiments with mixtures of proteins in assembly reactions provide novel analyses of the mature core assembly mechanism. Our results support a model in which initial assembly products serve as scaffolds for further assembly by converting incoming subunits to assembly proficient conformations, while mutant subunits increase the probability of assembly termination events.

Nematic ordering of polymers in confined geometry applied to DNA packaging in viral capsids

A density functional theory of the spatial distribution and biaxial nematic order of polymers of arbitrary length and rigidity inside a spherical cavity is proposed. The local order of different chain segments is considered as an alignment to a spatially varying director field of cylindrical symmetry. The steric interactions are taken into account in the second virial approximation. Polymer density and orientational order distributions inside the spherically cavity are the principal results. It was found that short and flexible polymers were located at the center of the sphere and were orientationaly disordered. Upon increasing polymer length and/or polymer rigidity, the location of the polymer was continuously shifted toward the surface of the spherical cavity and the polymer segments became gradually more aligned. Parameters have been selected to model the behavior of genomes in spherical viral capsids.

Normal mode analysis and applications in biological physics

Normal mode analysis has become a popular and often used theoretical tool in the study of functional motions in enzymes, viruses, and large protein assemblies. The use of normal modes in the study of these motions is often extremely fruitful since many of the functional motions of large proteins can be described using just a few normal modes which are intimately related to the overall structure of the protein. In this review, we present a broad overview of several popular methods used in the study of normal modes in biological physics including continuum elastic theory, the elastic network model, and a new all-atom method, recently developed, which is capable of computing a subset of the low frequency vibrational modes exactly. After a review of the various methods, we present several examples of applications of normal modes in the study of functional motions, with an emphasis on viral capsids.

Role of protein interactions in defining HIV-1 viral capsid shape and stability: a coarse-grained analysis

Coarse-grained models of the HIV-1 CA dimer are constructed based on all-atom molecular dynamics simulations. Coarse-grained representations of the capsid shell, which is composed of approximately 1500 copies of CA proteins, are constructed and their stability is examined. A key interaction between carboxyl and hexameric amino terminal domains is shown to generate the curvature of the capsid shell. It is demonstrated that variation of the strength of this interaction for different subunits in the lattice can cause formation of asymmetric, conical-shaped closed capsid shells, and it is proposed that variations, in the structure of the additional carboxyl-amino terminal binding interface during self-assembly, are important aspects of capsid cone formation. These results are in agreement with recent structural studies of the capsid hexamer subunit, which suggest that variability in the binding interface is a cause of the differences in subunit environments that exist in a conical structure.

Emergent complexity from simple anisotropic building blocks: shells, tubes, and spirals

We describe a remarkably simple, generic, coarse-grained model involving anisotropic interactions, and characterize the global minima for clusters as a function of various parameters. Appropriate choices for the anisotropic interactions can reproduce a wide variety of complex morphologies as global minima, including spheroidal shells, tubular, helical and even head-tail morphologies, elucidating the physical principles that drive the assembly of these mesoscopic structures. Our model captures several experimental observations, such as the existence of competing morphologies, capsid polymorphism, and the effect of scaffolding proteins on capsid assembly.

Invariant polymorphism in virus capsid assembly

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

Fluorescence fluctuation spectroscopy on viral-like particles reveals variable gag stoichiometry

Fluorescence fluctuation spectroscopy determines the brightness, size, and concentration of fluorescent particles from the intensity bursts generated by individual particles passing through a small observation volume. Brightness provides a measure of the number of fluorescently labeled proteins within a complex and has been used previously to determine the stoichiometry of small oligomers in cells. We extend brightness analysis to large macromolecular protein complexes containing thousands of proteins and determine their stoichiometry. This study investigates viral-like particles (VLP) formed from human immunodeficiency virus type 1 (HIV-1) Gag protein expressed in COS-1 cells using fluorescence fluctuation spectroscopy to determine the stoichiometry of HIV-1 Gag within the particles. Control experiments establish that the stoichiometry and size of VLPs are not influenced by labeling of HIV-1 Gag with a fluorescent protein. The experiments further show that the brightness scales linearly with the amount of labeled Gag within the particle. Brightness analysis shows that the Gag stoichiometry of VLPs formed in COS-1 cells is not constant, but varies with the amount of transfected DNA plasmid. We observed HIV-1 Gag stoichiometries ranging from approximately 750 to approximately 2500, whereas the size of the VLPs remains unchanged. This result indicates that large areas of the VLP membrane are void of Gag protein. Therefore, a closed layer of HIV-1 Gag at the membrane is not required for VLP production. This study shows that brightness analysis has the potential to become an important tool for investigating large molecular complexes by providing quantitative information about their size and composition.

Nonequilibirum assembly, retroviruses, and conical structures

We study the spontaneous assembly of viral shells composed of several identical subunits under nonequilibrium conditions. We find that within the basic continuum elasticity framework, the nonequilibrium assembly process is able to predict the formation of structures pertinent to retroviruses. Our minimal model of assembly yields a unified one-dimensional phase diagram in which the appearance of spherical, irregular, conical and cylindrical structures of retroviruses is seen to be governed by the spontaneous curvature of protein subunits.

Retroviral capsid assembly: a role for the CA dimer in initiation

In maturing retroviral virions, CA protein assembles to form a capsid shell that is essential for infectivity. The structure of the two folded domains [N-terminal domain (NTD) and C-terminal domain (CTD)] of CA is highly conserved among various retroviruses, and the capsid assembly pathway, although poorly understood, is thought to be conserved as well. In vitro assembly reactions with purified CA proteins of the Rous sarcoma virus (RSV) were used to define factors that influence the kinetics of capsid assembly and provide insights into underlying mechanisms. CA multimerization was triggered by multivalent anions providing evidence that in vitro assembly is an electrostatically controlled process. In the case of RSV, in vitro assembly was a well-behaved nucleation-driven process that led to the formation of structures with morphologies similar to those found in virions. Isolated RSV dimers, when mixed with monomeric protein, acted as efficient seeds for assembly, eliminating the lag phase characteristic of a monomer-only reaction. This demonstrates for the first time the purification of an intermediate on the assembly pathway. Differences in the intrinsic tryptophan fluorescence of monomeric protein and the assembly-competent dimer fraction suggest the involvement of the NTD in the formation of the functional dimer. Furthermore, in vitro analysis of well-characterized CTD mutants provides evidence for assembly dependence on the second domain and suggests that the establishment of an NTD-CTD interface is a critical step in capsid assembly initiation. Overall, the data provide clear support for a model whereby capsid assembly within the maturing virion is dependent on the formation of a specific nucleating complex that involves a CA dimer and is directed by additional virion constituents.

Energy landscapes for shells assembled from pentagonal and hexagonal pyramids

We present new rigid body potentials that should favour efficient self-assembly of pentagonal and hexagonal pyramids into icosahedral shells over a wide range of temperature. By adding an extra repulsive site opposite the existing apex sites of the pyramids considered in a previously published model, frustrated energy landscapes are transformed into systems identified with self-assembling properties. The extra interaction may be considered analogous to a hydrophobic-hydrophilic repulsion, as in micelle formation.

Generalized structural polymorphism in self-assembled viral particles

The protein shells, called capsids, of nearly all spherical viruses adopt icosahedral symmetry; however, self-assembly of such empty structures often occurs with multiple misassembly steps resulting in the formation of aberrant structures. Using simple models that represent the coat proteins preassembled in the two different predetermined species that are common motifs of viral capsids (i.e., pentameric and hexameric capsomers), we perform molecular dynamics simulations of the spontaneous self-assembly of viral capsids of different sizes containing T = 1,3,4,7,9,12,13,16, and 19 proteins in their icosahedral repeating unit. We observe, in addition to icosahedral capsids, a variety of nonicosahedral yet highly ordered and enclosed capsules. Such structural polymorphism is demonstrated to be an inherent property of the coat proteins, independent of the capsid complexity and the elementary kinetic mechanisms. Moreover, there exist two distinctive classes of polymorphic structures: aberrant capsules that are larger than their respective icosahedral capsids, in T = 1-7 systems; and capsules that are smaller than their respective icosahedral capsids when T = 7-19. Different kinetic mechanisms responsible for self-assembly of those classes of aberrant structures are deciphered, providing insights into the control of the self-assembly of icosahedral capsids.

Elastic properties and mechanical stability of chiral and filled viral capsids

The elasticity and mechanical stability of empty and filled viral capsids under external force loading are studied in a combined analytical and numerical approach. We analyze the influence of capsid structure and chirality on the mechanical properties. We find that generally skew shells have lower stretching energy. For large Föppl-von Kármán numbers gamma (gamma approximately 10(5)), skew structures are stiffer in their elastic response than nonchiral ones. The discrete structure of the capsules not only leads to buckling for large gamma but also influences the breakage behavior of capsules below the buckling threshold: the rupture force shows a gamma1/4 scaling rather than a gamma1/2 scaling as expected from our analytical results for continuous shells. Filled viral capsids are exposed to internal anisotropic pressure distributions arising from regularly packaged DNA coils. We analyze their influence on the elastic properties and rupture behavior and we discuss possible experimental consequences. Finally, we numerically investigate specific sets of parameters corresponding to specific phages such as phi29 and cowpea chlorotic mottle virus (CCMV). From the experimentally measured spring constants we make predictions about specific material parameters (such as bending rigidity and Young's modulus) for both empty and filled capsids.

Critical role of conserved hydrophobic residues within the major homology region in mature retroviral capsid assembly

During retroviral maturation, the CA protein undergoes dramatic structural changes and establishes unique intermolecular interfaces in the mature capsid shell that are different from those that existed in the immature precursor. The most conserved region of CA, the major homology region (MHR), has been implicated in both immature and mature assembly, although the precise contribution of the MHR residues to each event has been largely undefined. To test the roles of specific MHR residues in mature capsid assembly, an in vitro system was developed that allowed for the first-time formation of Rous sarcoma virus CA into structures resembling authentic capsids. The ability of CA to assemble organized structures was destroyed by substitutions of two conserved hydrophobic MHR residues and restored by second-site suppressors, demonstrating that these MHR residues are required for the proper assembly of mature capsids in addition to any role that these amino acids may play in immature particle assembly. The defect caused by the MHR mutations was identified as an early step in the capsid assembly process. The results provide strong evidence for a model in which the hydrophobic residues of the MHR control a conformational reorganization of CA that is needed to initiate capsid assembly and suggest that the formation of an interdomain interaction occurs early during maturation.

Radial distribution of RNA genomes packaged inside spherical viruses

The problem of RNA genomes packaged inside spherical viruses is studied. The RNA-capsid attraction is assumed to be nonspecific and occurs at the inner capsid surface only. For weak attraction, RNA concentration is maximum at the center of the capsid to maximize their configurational entropy. For stronger attraction, RNA concentration peaks near the capsid surface. In the latter case, the competition between the branching of RNA secondary structure and its adsorption to the inner capsid results in the formation of a dense layer of RNA near capsid surface. The layer thickness is a slowly varying (logarithmic) function of the capsid inner radius. Consequently, the amount of RNA packaged is proportional to the capsid area (or the number of proteins) instead of its volume. The numerical profiles describe reasonably well the observed RNA concentration profiles of various viruses.

Self-assembly approaches to nanomaterial encapsulation in viral protein cages

A perspective on abiotic material encapsulation inside virus capsids is provided. The emphasis is on the physical principles of virus assembly relevant to packaging, strategies for encapsulation and capsid modification, and on emerging applications.

Role of electrostatic interactions in the assembly of empty spherical viral capsids

We examine the role of electrostatic interactions in the assembly of empty spherical viral capsids. The charges on the protein subunits that make the viral capsid mutually interact and are expected to yield electrostatic repulsion acting against the assembly of capsids. Thus, attractive protein-protein interactions of nonelectrostatic origin must act to enable the capsid formation. We investigate whether the interplay of repulsive electrostatic and attractive interactions between the protein subunits can result in the formation of spherical viral capsids of a preferred radius. For this to be the case, we find that the attractive interactions must depend on the angle between the neighboring protein subunits (i.e., on the mean curvature of the viral capsid) so that a particular angle(s) is (are) preferred energywise. Our results for the electrostatic contributions to energetics of viral capsids nicely correlate with recent experimental determinations of the energetics of protein-protein contacts in the hepatitis B virus [P. Ceres A. Zlotnick, Biochemistry 41, 11525 (2002)].

Viral capsids: mechanical characteristics, genome packaging and delivery mechanisms

The main functions of viral capsids are to protect, transport and deliver their genome. The mechanical properties of capsids are supposed to be adapted to these tasks. Bacteriophage capsids also need to withstand the high pressures the DNA is exerting onto it as a result of the DNA packaging and its consequent confinement within the capsid. It is proposed that this pressure helps driving the genome into the host, but other mechanisms also seem to play an important role in ejection. DNA packaging and ejection strategies are obviously dependent on the mechanical properties of the capsid. This review focuses on the mechanical properties of viral capsids in general and the elucidation of the biophysical aspects of genome packaging mechanisms and genome delivery processes of double-stranded DNA bacteriophages in particular.

Defect scars on flexible surfaces with crystalline order

The crystallography of two-dimensional particle packings on flexible surfaces of spherical topology is investigated. Examples are viral capsids and crystalline vesicles. Computer simulations of dynamically triangulated surfaces are employed to study the shape and structure of lattice defects as a function of the Föppl-von Kármán number gamma. We find that grain-boundary scars become much more fuzzy with increasing temperature, that the size of grain-boundary scars saturates with increasing vesicle radius, and that the buckling transition shifts to higher values of gamma due to the presence of scars.

A precise packing sequence for self-assembled convex structures

Molecular simulations of the self-assembly of cone-shaped particles with specific, attractive interactions are performed. Upon cooling from random initial conditions, we find that the cones self-assemble into clusters and that clusters comprised of particular numbers of cones (e.g., 4-17, 20, 27, 32, and 42) have a unique and precisely packed structure that is robust over a range of cone angles. These precise clusters form a sequence of structures at specific cluster sizes (a "precise packing sequence") that for small sizes is identical to that observed in evaporation-driven assembly of colloidal spheres. We further show that this sequence is reproduced and extended in simulations of two simple models of spheres self-assembling from random initial conditions subject to convexity constraints, including an initial spherical convexity constraint for moderate- and large-sized clusters. This sequence contains six of the most common virus capsid structures obtained in vivo, including large chiral clusters and a cluster that may correspond to several non-icosahedral, spherical virus capsids obtained in vivo. Our findings suggest that this precise packing sequence results from free energy minimization subject to convexity constraints and is applicable to a broad range of assembly processes.

Irreversible growth model for virus capsid assembly

We model the spontaneous assembly of a capsid (a virus' closed outer shell) from many copies of identical units, using entirely irreversible steps and only information local to the growing edge. Our model is formulated in terms of (i) an elastic Hamiltonian with stretching and bending stiffness and a spontaneous curvature, and (ii) a set of rate constants for the addition of new units or bonds. An ensemble of highly irregular capsids is generated, unlike the well-known icosahedrally symmetric viruses, but (we argue) plausible as a way to model the irregular capsids of retroviruses such as HIV. We found that (i) the probability of successful capsid completion decays exponentially with capsid size; (ii) capsid size depends strongly on spontaneous curvature and weakly on the ratio of the bending and stretching elastic stiffnesses of the shell; (iii) the degree of localization of Gaussian curvature (a measure of facetedness) depends heavily on the ratio of elastic stiffnesses.