Coronavirus rotational diffusivity

Abstract

Just 11 weeks after the confirmation of first infection, one team had already discovered and published [D. Wrapp et al., "Cryo-EM structure of the 2019-nCoV spike in the prefusion conformation,"Science 367(6483), 1260-1263 (2020)] in exquisite detail about the new coronavirus, along with how it differs from previous viruses. We call the virus particle causing the COVID-19 disease SARS-CoV-2, a spherical capsid covered with spikes termed peplomers. Since the virus is not motile, it relies on its own random thermal motion, specifically the rotational component of this thermal motion, to align its peplomers with targets. The governing transport property for the virus to attack successfully is thus the rotational diffusivity. Too little rotational diffusivity and too few alignments are produced to properly infect. Too much, and the alignment intervals will be too short to properly infect, and the peplomer is wasted. In this paper, we calculate the rotational diffusivity along with the complex viscosity of four classes of virus particles of ascending geometric complexity: tobacco mosaic, gemini, adeno, and corona. The gemini and adeno viruses share icosahedral bead arrangements, and for the corona virus, we use polyhedral solutions to the Thomson problem to arrange its peplomers. We employ general rigid bead-rod theory to calculate complex viscosities and rotational diffusivities, from first principles, of the virus suspensions. We find that our ab initio calculations agree with the observed complex viscosity of the tobacco mosaic virus suspension. From our analysis of the gemini virus suspension, we learn that the fine detail of the virus structure governs its rotational diffusivity. We find the characteristic time for the adenovirus from general rigid bead-rod theory. Finally, from our analysis of the coronavirus suspension, we learn that its rotational diffusivity descends monotonically with its number of peplomers.