Molecular Dynamics Simulations of Deformable Viral Capsomers

Most coarse-grained models of individual capsomers associated with viruses employ rigid building blocks that do not exhibit shape adaptation during self-assembly. We develop a coarse-grained general model of viral capsomers that incorporates their stretching and bending energies while retaining many features of the rigid-body models, including an overall trapezoidal shape with attractive interaction sites embedded in the lateral walls to favor icosahedral capsid assembly. Molecular dynamics simulations of deformable capsomers reproduce the rich self-assembly behavior associated with a general T=1 icosahedral virus system in the absence of a genome. Transitions from non-assembled configurations to icosahedral capsids to kinetically-trapped malformed structures are observed as the steric attraction between capsomers is increased. An assembly diagram in the space of capsomer-capsomer steric attraction and capsomer deformability reveals that assembling capsomers of higher deformability into capsids requires increasingly large steric attraction between capsomers. Increasing capsomer deformability can reverse incorrect capsomer-capsomer binding, facilitating transitions from malformed structures to symmetric capsids; however, making capsomers too soft inhibits assembly and yields fluid-like structures.