Shape selection and mis-assembly in viral capsid formation by elastic frustration
Figures

Sketch of the different structures considered in this study.
(a) A hexagonally-ordered spherical cap of radius and geodesic radius without defects; (b) a spherical cap with a single disclination at the center (as shown) or multiple disclinations; (c) a rectangular ribbon with length and width , that it is called belt when the length becomes ; and (d) a cylindrical patch with size , that eventually becomes a cylinder of radius . In the bending-dominated regime, .

Comparison of free energy landscapes for different structures.
Free energy of formation versus the radius of the patch x in the bending-dominated regime for a defectless spherical shell (blue line), a spherical shell with defects (green), and a ribbon/belt (orange) for and . The optimal structure is the one with the minimum energy, which is the belt in this case. The dashed line represents the unfrustrated decrease of energy expected by the classical nucleation picture for the defectless spherical cap in the absence of elastic stresses. The inset zooms the nucleation barrier located at small patch sizes.

Assembly phase diagrams.
Phase diagrams of the most stable structures in terms of the scaled chemical potential and the scaled line tension for different values of the FvK number: a) , corresponding to the bending-dominated regime, (b) , (c) , and d) . Three possible equilibrium regions are present: belts (i.e. closed ribbons, in orange), frustrated capsids with one defect (red), closed shells with defects (green), and cylinders (purple). Additionally, a region corresponding to a metastable spherical cap without defects (blue) is shown only in (a). In the white region, the equilibrium state corresponds to disaggregated individual capsomers.

Comparison of simulation results with the theoretical phase diagram.
(a) Phase diagram of the most stable structure in terms of the scaled chemical potential and the scaled line tension for . Snapshots of the final outcome of the simulation for different initial capsomer concentrations for: b) , obtained with and a , potential. By increasing the concentration of capsomers one goes from a dissassembled state to the formation of a ribbon to a spherical shell with defects. (The ribbon and spherical shell are not closed in the snapshots due to the limited number of capsomers and finite simulation time). (c) , obtained with and a , potential. As concentration increases, one goes from a shell with one defect to a complete shell with many defects. By increasing the FvK number to (by setting ), the simulations clearly forms cylindrical tubes. (d) , obtained with and a , potential. In this case, the sequence is: disagreggated, metastable defectless shell, and closed spherical shell with many defects (a partial shell is shown). The simulation results agree qualitatively with the predictions from the theoretical phase diagram.

Free energy of formation of a spherical cap without defects, Equation A4, versus the radius of the patch for and different values of the scaled chemical potential , illustrating the situations in which no assembly is possible (short dashed line, ), a geometrically frustrated metastastable cap (dashed line, ) or stable (solid line, ) finite shell.

Scaled optimal radius of a spherical cap without defects as a function of the patch size for a FvK number .
The solid line is the exact solution, Equation A29, the short dashed line is the approximation for large shells or FvK, Equation A31, and the dashed line is the approximation for small shells or FvK, Equation A30.

Other structures obtained in the simulations.
(a) Bullet-shaped shell and (b) conical shell obtained in two different repetitions for , , with , , , and ; (c) T = 13 icosahedral shell obtained for , , with , , , and ; (d) coexistence between a cylinder and a partial spherocylindrical shell obtained for , , with , , , and ; (e) branched ribbon-like structure (Köhler et al., 2016) obtained at low values of the scaled line tension for , , with , , , and .
Tables
Estimates of the main geometric and elastic properties of different non-enveloped empty viral capsids.
The Young’s modulus has been evaluated from AFM nanoindentation experiments (Mateu, 2012; Michel et al., 2006; Ivanovska et al., 2007; Sae-Ueng et al., 2014) and, for SV40, from the experimental spring constant (van Rosmalen et al., 2018) using the standard thin shell formula (Ivanovska et al., 2004). The 2D Young’s Modulus was calculated as ; the effective diameter of the capsomers as (Santolaria, 2011) , where is the triangulation number; the line tension as (Luque et al., 2012) considering a typical binding energy ; and the FvK number as , with (Buenemann and Lenz, 2008).
Virus | T-number | Diameter (nm) | Thickness h (nm) | E (Gpa) | Y (N/m) | σ (nm) | Scaled line tension λ | Föppl-von Karman γ |
---|---|---|---|---|---|---|---|---|
CCMV | 3 | 28 | 3.8 | 0.14 | 0.53 | 5.9 | 0.00107 | 148 |
λ Procapsid | 7 | 50 | 4.0 | 0.16 | 0.64 | 6.8 | 0.000436 | 427 |
λ Capsid | 7 | 63 | 1.8 | 1.0 | 1.8 | 8.6 | 0.0000976 | 3344 |
SV40 | 7 | 45 | 6.0 | 0.033 | 0.2 | 6.1 | 0.00174 | 152 |