Breaking a virus: Identifying molecular level failure modes of a viral capsid by multiscale modeling

Abstract

We use coarse-grained (CG) simulations to study the deformation of empty Cowpea Chlorotic Mottle Virus (CCMV) capsids under uniaxial compression, from the initial elastic response up to capsid breakage. Our CG model is based on the MARTINI force field and has been amended by a stabilizing elastic network, acting only within individual proteins, that was tuned to capture the fluctuation spectrum of capsid protein dimers, obtained from all atom simulations. We have previously shown that this model predicts force-compression curves that match AFM indentation experiments on empty CCMV capsids. Here we investigate details of the actual breaking events when the CCMV capsid finally fails. We present a symmetry classification of all relevant protein contacts and show that they differ significantly in terms of stability. Specifically, we show that interfaces which break readily are precisely those which are believed to form last during assembly, even though some of them might share the same contacts as other non-breaking interfaces. In particular, the interfaces that form pentamers of dimers never break, while the virtually identical interfaces within hexamers of dimers readily do. Since these units differ in the large-scale geometry and, most noticeably, the cone-angle at the center of the 5- or 6-fold vertex, we propose that the hexameric unit fails because it is pre-stressed. This not only suggests that hexamers of dimers form less frequently during the early stages of assembly; it also offers a natural explanation for the well-known β-barrel motif at the hexameric center as a post-aggregation stabilization mechanism. Finally, we identify those amino acid contacts within all key protein interfaces that are most persistent during compressive deformation of the capsid, thereby providing potential targets for mutation studies aiming to elucidate the key contacts upon which overall stability rests.