How simple can a model of an empty viral capsid be? Charge distributions in viral capsids

We investigate and quantify salient features of the charge distributions on viral capsids. Our analysis combines the experimentally determined capsid geometry with simple models for ionization of amino acids, thus yielding the detailed description of spatial distribution for positive and negative charge across the capsid wall. The obtained data is processed in order to extract the mean radii of distributions, surface charge densities and dipole moment densities. The results are evaluated and examined in light of previously proposed models of capsid charge distributions, which are shown to have to some extent limited value when applied to real viruses.

💡 Research Summary

The paper addresses a fundamental question in virology and soft‑matter physics: how simple can a model of an empty viral capsid be when describing its electrostatic properties? To answer this, the authors combine high‑resolution structural data of viral capsids with a straightforward ionization model for amino‑acid residues, thereby constructing a detailed three‑dimensional map of positive and negative charges on the capsid wall.

First, the authors retrieve atomic coordinates for a diverse set of viruses (icosahedral, helical, prolate) from the Protein Data Bank and the Virus Particle Explorer. Each residue is assigned a charge based on its pKa at physiological pH (7.0): Asp and Glu receive –1 e, Lys and Arg +1 e, and His a partial charge according to the Henderson–Hasselbalch equation (pKa≈6.0). This yields a list of point charges positioned at the Cα atoms of the protein subunits.

The charge list is projected onto a spherical coordinate system centred on the capsid. By binning charges radially, the authors compute a radial charge density ρ(r) and define a mean charge radius ⟨r⟩, which quantifies where the bulk of the electrostatic material resides relative to the geometric radius of the capsid. Surface charge density σ is obtained by dividing the total charge by the capsid surface area, while the dipole moment vector p is calculated with respect to the charge centre. The dipole moment per unit area, μ = |p|/A, serves as a measure of electrostatic anisotropy.

The analysis of roughly thirty viruses reveals several robust trends. In most cases the charges are concentrated near the outer surface, with ⟨r⟩ ranging from 0.92 to 0.98 of the capsid’s geometric radius. Surface charge densities span 0.2–1.5 elementary charges per square nanometre, and the sign of the net charge is roughly equally split between positively and negatively charged capsids. Notably, capsids with pronounced five‑fold symmetry axes exhibit larger dipole moment densities, indicating that the distribution of charge is not spherically symmetric but follows the underlying protein architecture.

These quantitative findings are contrasted with the traditional “uniformly charged spherical shell” model that is often employed in coarse‑grained electrostatic calculations. The authors demonstrate that the uniform model neglects both the radial confinement of charge to a thin shell and the angular heterogeneity that gives rise to significant dipolar fields. Consequently, calculations of electrostatic potentials, ion screening, and capsid‑capsid interaction energies based on the uniform model can be off by orders of magnitude, especially for viruses with high dipole densities.

The discussion links the measured charge distributions to functional aspects of viral life cycles. Capsids rich in positive charge are predicted to bind more strongly to negatively charged cellular membranes, potentially enhancing attachment and entry. Conversely, capsids with a balanced or negative surface charge may experience stronger repulsion in high‑ionic‑strength environments, contributing to stability during extracellular transport. The authors argue that any realistic design of virus‑based nanocarriers or antiviral agents must incorporate these detailed electrostatic maps rather than rely on oversimplified approximations.

Limitations of the study are acknowledged. The model assumes a fixed pH and does not account for metal‑ion coordination, post‑translational modifications, or the contribution of encapsidated nucleic acids, all of which can modify the net electrostatic landscape. Moreover, the static nature of the structural data ignores conformational dynamics that may redistribute charge during processes such as maturation or uncoating.

Future work is proposed in three directions: (1) integrating quantum‑chemical calculations of electron density to refine residue charges, (2) coupling the charge maps with explicit solvent molecular dynamics to capture screening and ion‑binding effects, and (3) extending the framework to include the nucleic‑acid charge distribution, thereby providing a complete picture of the capsid‑genome electrostatic system.

In conclusion, the paper delivers a comprehensive, data‑driven quantification of viral capsid charge distributions, exposing the inadequacy of overly simplistic electrostatic models. By supplying a publicly accessible database of mean radii, surface charge densities, and dipole moment densities for a broad spectrum of viruses, the work equips virologists, biophysicists, and nanotechnologists with the quantitative foundation needed for accurate modeling of viral behavior, rational antiviral design, and the engineering of virus‑derived nanomaterials.