Minimal stabilizing bond sets for Caspar-Klug capsid spring networks: a systematic three-dimensional choice achieving exactly six zero modes

In the coarse-grained normal-mode model of Englert, Peeters and Taormina [Phys. Rev. E 78, 031908 (2008)], an icosahedral Caspar-Klug viral capsid with triangulation number $T=n$ is represented by $60n$ point masses at the coat-protein centres of mass, joined by springs wherever the VIPER database of structural virology reports an inter-protein association energy. The authors showed that stable capsids in this model always exhibit a low-frequency plateau of 24 near-zero modes with universal icosahedral representation content, and they left open the mathematical question of how to choose and place, in three dimensions, the minimal set of stabilizing inter-protein bonds: counting requires at least $60(3n-q)-6$ independent bonds beyond the on-face constraints, but no systematic construction or characterization is known.

Background

Most viruses package their genome in a protein shell, the capsid, which for a large class of viruses assembles with icosahedral symmetry. Caspar and Klug (1962) classified such shells by a triangulation number $T=n$ with $n=h^2+hk+k^2$ for non-negative integers $h,k$ not both zero: a $T=n$ capsid contains $60n$ copies of the coat protein, organised as twelve pentamers (clusters of five subunits) at the vertices of an icosahedron and $10(n-1)$ hexamers (clusters of six) at the global three-fold and local six-fold symmetry axes.

Large, slow conformational changes of viruses — maturation and genome release — are widely believed to be dominated by a few collective, low-frequency vibrational modes, which motivates normal mode analysis (NMA): diagonalising the harmonic, i.e. quadratic, approximation of the potential energy about an equilibrium structure. The standard coarse-graining is the elastic network model introduced by Tirion (1996), in which each residue or subunit becomes a point mass and springs connect all pairs of masses closer than a distance cutoff; it uses a single spring constant and therefore cannot distinguish strong from weak inter-protein bonds. Englert, Peeters and Taormina (2008) instead used the VIPERdb database of structural virology (Shepherd et al. 2006), which lists pairwise inter-protein association energies, to decide which springs exist and to weight them relatively: spring constants $\kappa_{mn}=\rho_{mn}\kappa$, where $\rho_{mn}$ is the association energy of the protein pair $(m,n)$ divided by the largest listed value and $\kappa$ is a single free scale, reflecting that the absolute values of the database energies are not trusted. With each protein a unit point mass at its centre of mass (from crystallographic coordinates), a $T=n$ capsid has $180n$ degrees of freedom and a stiffness (force) matrix — the Hessian of the quadratic spring energy — whose null space always contains the six trivial zero modes, i.e. the three global translations and three global rotations. In this model a capsid is called stable precisely when its stiffness matrix has exactly six zero eigenvalues.

Their central result is a symmetry-based explanation of a robust spectral feature: in every stable capsid they studied, the low-frequency spectrum contains a plateau of 24 near-zero modes. The cleanest idealisation merges the $3n$ proteins of each of the twenty icosahedral faces into a single mass per face; the masses sit at the vertices of the dodecahedron (the polyhedron with 20 vertices, 30 edges and 12 pentagonal faces, dual to the icosahedron), and nearest-neighbour springs form a cage with 30 edges. This cage has $3\times 20-6=54$ genuine degrees of freedom against only 30 constraints, hence 24 non-trivial zero modes; a group-theoretical comparison of displacement representations shows that these transform as $\Gamma_{3'}^{+}+\Gamma_{3'}^{-}+\Gamma_{4}^{+}+\Gamma_{4}^{-}+\Gamma_{5}^{+}+\Gamma_{5}^{-}$, a sum of non-singlet irreducible representations of the full icosahedral group $I_h$ (the 120-element group obtained by adjoining inversion to the 60 proper rotations; $\Gamma_{d}^{\pm}$ denotes a $d$-dimensional irreducible representation of parity $\pm$, with the prime distinguishing an inequivalent representation of the same dimension). Restoring the full, realistic bond structure lifts these zero modes only slightly, producing the plateau, while the first fully symmetric mode always lies higher. Related continuum analyses of soft modes of icosahedral shells (Widom, Lidmar and Nelson 2007) predate this work and address the buckling regime rather than the discrete bond network.

For general $T=n$ the same reasoning is combinatorial. If the $9n$ degrees of freedom within one face are constrained by $3q$ on-face links ($q$ an integer fixed by the capsid's bond structure; e.g. $q=4$ for Rice Yellow Mottle Virus after discarding its two weakest bonds, $q=12$ for phage HK97), then without inter-face bonds the capsid has $N_0=(9n-3q)\times 20=60(3n-q)$ zero modes, so a stable capsid requires at least $60(3n-q)-6$ further independent bonds. The authors state that how these bonds must be chosen and placed in three dimensions 'is a mathematical question which requires further investigation', and their companion paper (Peeters and Taormina 2009) calls the minimisation of the number of bonds needed for stability of a three-dimensional spring-mass model 'a very interesting and deep mathematical question'.

This is a rigidity-theory question. A network of point masses joined by fixed-length bars is isostatic, or minimally rigid, when it is rigid but loses rigidity upon removal of any single bar. In two dimensions, Laman's theorem (1970) characterises such networks combinatorially (every subset of $v$ vertices spans at most $2v-3$ bars, with equality overall), which underlies the linear-time pebble game algorithm of Jacobs and Hendrickson (1997) that decides two-dimensional network stability. No analogue of either result is known in three dimensions — a famous gap, already noted in the companion paper — and the capsid version adds a symmetry constraint: admissible bond sets are unions of orbits of the icosahedral group, so bond counts come in multiples of 30 or 60, and symmetry can force additional degeneracies. Representation-level necessary conditions for symmetric isostatic frameworks exist (Connelly, Fowler, Guest, Schulze and Whiteley 2009) but have not been applied to capsid geometries. Pebble-game rigidity analysis has been applied to one capsid, cowpea chlorotic mottle virus (Hespenheide, Jacobs and Thorpe 2004), but as an analysis of a given network, not a construction of minimal stabilising ones. The unresolved target below asks for exactly this missing systematic construction and its consequences for the universality of the 24-mode plateau.

Problem Statement

The problem is posed within the harmonic point-mass spring-network model of Englert, Peeters and Taormina (2008): $60n$ unit point masses at the coat-protein centres of mass of a Caspar-Klug capsid with triangulation number $T=n$ ($n=h^2+hk+k^2$), springs assigned according to the capsid's inter-protein association-energy data, stiffness matrix the Hessian of the corresponding quadratic potential, and stability defined as the stiffness matrix having exactly six zero eigenvalues (the three global translations and three global rotations).

Primary question. Given the icosahedrally symmetric mass configuration and the on-face constraint structure of a $T=n$ capsid ($3q$ independent on-face bonds per icosahedral face, $q$ fixed by the bond data), give a systematic three-dimensional choice and placement of a minimal set of stabilizing inter-protein bonds for general $T=n$ such that the stiffness matrix has exactly six zero modes. Since the on-face constraints alone leave

$$N_0=(9n-3q)\times 20=60(3n-q)$$

zero modes. The target consists of three parts; (i) fixes the admissible class of bond sets:

(i) admissible bond sets are unions of orbits of the icosahedral group (hence of sizes in multiples of 30 or 60); such a set is stabilizing when its addition makes the null space of the stiffness matrix exactly six-dimensional; (ii) establish minimality — determine, for each $T=n$, the minimum orbit-union set of inter-face bonds whose addition leaves the stiffness matrix with exactly six zero modes, giving the construction where such a set exists and a proof of non-existence where it does not; the count $60(3n-q)-6$ is a stated lower bound on the number of independent inter-face constraints, not itself the target of the construction, since admissible bond sets are unions of icosahedral orbits of size a multiple of 30 or 60 and the realized minimum may exceed it; (iii) hold for arbitrary $T=n$ rather than only for specific capsids, and explain the relationship between minimal stabilization and the source's universal observation that VIPER-guided near-minimal bond sets always leave a low-frequency plateau of exactly 24 modes with icosahedral content $\Gamma_{3'}^{+}+\Gamma_{3'}^{-}+\Gamma_{4}^{+}+\Gamma_{4}^{-}+\Gamma_{5}^{+}+\Gamma_{5}^{-}$. The authors' heuristic — keep the strongest $2(3n-q)-1$ edge-crossing bonds per icosahedral edge (bonds joining proteins of the two faces that meet along that edge) — has been verified only numerically on specific capsids (satellite tobacco mosaic virus (STMV), Rice Yellow Mottle Virus (RYMV), bacteriophage HK97 and others) and is not a solution.

Secondary question (same objective, raised in the companion paper Peeters and Taormina 2009). Provide a mathematical understanding of the anomalous low-frequency plateau of the all-pentamer capsid of simian virus 40 (SV40), in particular whether its 30 near-zero modes are a generic signature of all-pentamer capsids rather than the 24-mode structure. The secondary question may be resolved independently of the primary question.

Both parts serve one scientific objective: a mathematical theory of minimal stabilization and zero-mode structure of icosahedrally symmetric capsid spring networks within the stated model class.

References

1. François Englert, Kasper Peeters, Anne Taormina, 'Twenty-four near-instabilities of Caspar-Klug viruses', Physical Review E 78, 031908 (2008). DOI: 10.1103/PhysRevE.78.031908; arXiv:0804.4275. https://doi.org/10.1103/PhysRevE.78.031908
2. Kasper Peeters, Anne Taormina, 'Group theory of icosahedral virus capsid vibrations: A top-down approach', Journal of Theoretical Biology 256, 607-624 (2009). DOI: 10.1016/j.jtbi.2008.10.019; arXiv:0806.1029. https://doi.org/10.1016/j.jtbi.2008.10.019
3. D. L. D. Caspar, A. Klug, 'Physical Principles in the Construction of Regular Viruses', Cold Spring Harbor Symposia on Quantitative Biology 27 (1962). DOI: 10.1101/sqb.1962.027.001.005. https://doi.org/10.1101/sqb.1962.027.001.005
4. Monique M. Tirion, 'Large Amplitude Elastic Motions in Proteins from a Single-Parameter, Atomic Analysis', Physical Review Letters 77, 1905-1908 (1996). DOI: 10.1103/PhysRevLett.77.1905. https://doi.org/10.1103/PhysRevLett.77.1905
5. Christopher M. Shepherd, Ivanka A. Borelli, Gaurav Lander, Preethi Natarajan, Vijay Siddavanahalli, Chandrajit Bajaj, John E. Johnson, Charles L. Brooks III, V. S. Reddy, 'VIPERdb: a relational database for structural virology', Nucleic Acids Research 34, D386-D389 (2006). DOI: 10.1093/nar/gkj032. https://doi.org/10.1093/nar/gkj032
6. M. Widom, J. Lidmar, David R. Nelson, 'Soft modes near the buckling transition of icosahedral shells', Physical Review E 76, 031911 (2007). DOI: 10.1103/PhysRevE.76.031911. https://doi.org/10.1103/PhysRevE.76.031911
7. G. Laman, 'On graphs and rigidity of plane skeletal structures', Journal of Engineering Mathematics 4, 331-340 (1970). DOI: 10.1007/BF01534980. https://doi.org/10.1007/BF01534980
8. Donald J. Jacobs, Bruce Hendrickson, 'An Algorithm for Two-Dimensional Rigidity Percolation: The Pebble Game', Journal of Computational Physics 137, 346-365 (1997). DOI: 10.1006/jcph.1997.5809. https://doi.org/10.1006/jcph.1997.5809
9. R. Connelly, P. W. Fowler, S. D. Guest, B. Schulze, W. J. Whiteley, 'When is a symmetric pin-jointed framework isostatic?', International Journal of Solids and Structures 46, 762-773 (2009). DOI: 10.1016/j.ijsolstr.2008.09.023. https://doi.org/10.1016/j.ijsolstr.2008.09.023
10. B. M. Hespenheide, D. J. Jacobs, M. F. Thorpe, 'Structural rigidity in the capsid assembly of cowpea chlorotic mottle virus', Journal of Physics: Condensed Matter 16, S5055-S5064 (2004). DOI: 10.1088/0953-8984/16/44/003. https://doi.org/10.1088/0953-8984/16/44/003
11. Eric C. Dykeman, Otto F. Sankey, 'Low Frequency Mechanical Modes of Viral Capsids: An Atomistic Approach', Physical Review Letters 100, 028101 (2008). DOI: 10.1103/PhysRevLett.100.028101. https://doi.org/10.1103/PhysRevLett.100.028101
12. Eric C. Dykeman, Otto F. Sankey, 'Atomistic modeling of the low-frequency mechanical modes and Raman spectra of icosahedral virus capsids', Physical Review E 81, 021918 (2010). DOI: 10.1103/PhysRevE.81.021918. https://doi.org/10.1103/PhysRevE.81.021918
13. Zheng Yang, Ivet Bahar, Michael Widom, 'Vibrational Dynamics of Icosahedrally Symmetric Biomolecular Assemblies Compared with Predictions Based on Continuum Elasticity', Biophysical Journal 96, 4438-4448 (2009). DOI: 10.1016/j.bpj.2009.03.016. https://doi.org/10.1016/j.bpj.2009.03.016
14. M. Martín-Bravo, J. M. Gomez Llorente, J. Hernández-Rojas, 'A minimal coarse-grained model for the low-frequency normal mode analysis of icosahedral viral capsids', Soft Matter 16, 3443-3455 (2020). DOI: 10.1039/d0sm00299b. https://doi.org/10.1039/d0sm00299b
15. Yin-Chen Hsieh, Frédéric Poitevin, Marc Delarue, Patrice Koehl, 'Comparative Normal Mode Analysis of the Dynamics of DENV and ZIKV Capsids', Frontiers in Molecular Biosciences 3, 85 (2016). DOI: 10.3389/fmolb.2016.00085. https://doi.org/10.3389/fmolb.2016.00085
16. Lei Zhang, C. Q. Ru, 'Free vibration of biopolymer spherical shells of high structural heterogeneity', AIP Advances 8, 075006 (2018). DOI: 10.1063/1.5036672. https://doi.org/10.1063/1.5036672
17. Lei Zhang, Yang Luo, Hao Wang, Yuanxin Luo, Jian Wu, 'The role of DNA/RNA core on free vibration characteristics of viral capsids', Journal of Applied Physics 138, 084703 (2025). DOI: 10.1063/5.0272224. https://doi.org/10.1063/5.0272224
18. Diane L. Lynch, Anna Pavlova, Zixing Fan, James C. Gumbart, 'Understanding Virus Structure and Dynamics through Molecular Simulations', Journal of Chemical Theory and Computation 19, 3025-3036 (2023). DOI: 10.1021/acs.jctc.3c00116. https://doi.org/10.1021/acs.jctc.3c00116
19. K. T. Tsen, Eric C. Dykeman, Otto F. Sankey, Nien-Tsung Lin, Shaw-Wei D. Tsen, Juliann G. Kiang, 'Observation of the low frequency vibrational modes of bacteriophage M13 in water by Raman spectroscopy', Virology Journal 3, 79 (2006). DOI: 10.1186/1743-422X-3-79. https://doi.org/10.1186/1743-422X-3-79
20. K. T. Tsen, Shaw-Wei D. Tsen, Otto F. Sankey, Juliann G. Kiang, 'Selective inactivation of micro-organisms with near-infrared femtosecond laser pulses', Journal of Physics: Condensed Matter 19, 472201 (2007). DOI: 10.1088/0953-8984/19/47/472201. https://doi.org/10.1088/0953-8984/19/47/472201
21. Eric C. Dykeman, Daryn Benson, K.-T. Tsen, Otto F. Sankey, 'Simulations of impulsive laser scattering of biological protein assemblies: Application to M13 bacteriophage', Physical Review E 80, 041909 (2009). DOI: 10.1103/PhysRevE.80.041909. https://doi.org/10.1103/PhysRevE.80.041909
22. Jeffrey C. Wigle, Eric A. Holwitt, Larry E. Estlack, Gary D. Noojin, Katharine E. Saunders, Valdislav V. Yakovlev, Benjamin A. Rockwell, 'No effect of femtosecond laser pulses on M13, E. coli, DNA, or protein', Journal of Biomedical Optics 19, 015008 (2014). DOI: 10.1117/1.JBO.19.1.015008. https://doi.org/10.1117/1.JBO.19.1.015008
23. Szu-Chi Yang, Huan-Chun Lin, Tzu-Ming Liu, Jen-Tang Lu, Wan-Ting Hung, Yuru Huang, Yi-Chun Tsai, Chuan-Liang Kao, Shih-Yuan Chen, Chi-Kuang Sun, 'Efficient Structure Resonance Energy Transfer from Microwaves to Confined Acoustic Vibrations in Viruses', Scientific Reports 5, 18030 (2015). DOI: 10.1038/srep18030. https://doi.org/10.1038/srep18030
24. S. Sirotkin, A. Mermet, M. Bergoin, V. Ward, J. L. Van Etten, 'Viruses as nanoparticles: Structure versus collective dynamics', Physical Review E 90, 022718 (2014). DOI: 10.1103/PhysRevE.90.022718. https://doi.org/10.1103/PhysRevE.90.022718
25. Jeffrey Burkhartsmeyer, Yanhong Wang, Kam Sing Wong, Reuven Gordon, 'Optical Trapping, Sizing, and Probing Acoustic Modes of a Small Virus', Applied Sciences 10, 394 (2020). DOI: 10.3390/app10010394. https://doi.org/10.3390/app10010394
26. Shui-Jing Tang, Mingjie Zhang, Jialve Sun, Jia-Wei Meng, Xiao Xiong, Qihuang Gong, Dayong Jin, Qi-Fan Yang, Yun-Feng Xiao, 'Single-particle photoacoustic vibrational spectroscopy using optical microresonators', Nature Photonics 17, 951-956 (2023). DOI: 10.1038/s41566-023-01264-3. https://doi.org/10.1038/s41566-023-01264-3
27. Yaqing Zhang, Rihan Wu, Md. Shahjahan, Canchai Yang, Dohun Pyeon, Elad Harel, 'Nanoscopic acoustic vibrational dynamics of a single virus captured by ultrafast spectroscopy', Proceedings of the National Academy of Sciences 122, e2420428122 (2025). DOI: 10.1073/pnas.2420428122. https://doi.org/10.1073/pnas.2420428122
28. L. E. Perotti, J. Rudnick, R. F. Bruinsma, W. S. Klug, 'Statistical Physics of Viral Capsids with Broken Symmetry', Physical Review Letters 115, 058101 (2015). DOI: 10.1103/PhysRevLett.115.058101. https://doi.org/10.1103/PhysRevLett.115.058101
29. Emilio Zappa, Giuliana Indelicato, Alberto Albano, Paolo Cermelli, 'A Ginzburg-Landau model for the expansion of a dodecahedral viral capsid' (2013). arXiv:1301.4036. https://arxiv.org/abs/1301.4036
30. Pierre-Philippe Dechant, 'Root systems & Clifford algebras: from symmetries of viruses to E8 & ADE correspondences' (2022). arXiv:2204.05718. https://arxiv.org/abs/2204.05718
31. Brian Bothner, Jonathan K. Hilmer, 'Probing Viral Capsids in Solution', in Structural Virology (RSC Biomolecular Sciences, 2010). DOI: 10.1039/9781849732239-00041. https://doi.org/10.1039/9781849732239-00041
32. Hon-Wai Leong, Lock Yue Chew, Kerson Huang, 'Normal modes and phase transition of the protein chain based on the Hamiltonian formalism', Physical Review E 82, 011915 (2010). DOI: 10.1103/PhysRevE.82.011915. https://doi.org/10.1103/PhysRevE.82.011915
33. Colin Brown, Anuradha Agarwal, Antoni Luque, 'pyCapsid: identifying dominant dynamics and quasi-rigid mechanical units in protein shells', Bioinformatics 40, btad761 (2023). DOI: 10.1093/bioinformatics/btad761. https://doi.org/10.1093/bioinformatics/btad761