Coronavirus hemagglutinin-esterase and spike proteins coevolve for functional balance and optimal virion avidity

Human coronaviruses OC43 and HKU1 are respiratory pathogens of zoonotic origin that have gained worldwide distribution. OC43 apparently emerged from a bovine coronavirus (BCoV) spillover. All three viruses attach to 9-O-acetylated sialoglycans via spike protein S with hemagglutinin-esterase (HE) acting as a receptor-destroying enzyme. In BCoV, an HE lectin domain promotes esterase activity toward clustered substrates. OC43 and HKU1, however, lost HE lectin function as an adaptation to humans. Replaying OC43 evolution, we knocked out BCoV HE lectin function and performed forced evolution-population dynamics analysis. Loss of HE receptor binding selected for second-site mutations in S, decreasing S binding affinity by orders of magnitude. Irreversible HE mutations led to cooperativity in virus swarms with low-affinity S minority variants sustaining propagation of high-affinity majority phenotypes. Salvageable HE mutations induced successive second-site substitutions in both S and HE. Apparently, S and HE are functionally interdependent and coevolve to optimize the balance between attachment and release. This mechanism of glycan-based receptor usage, entailing a concerted, fine-tuned activity of two envelope protein species, is unique among CoVs, but reminiscent of that of influenza A viruses. Apparently, general principles fundamental to virion-sialoglycan interactions prompted convergent evolution of two important groups of human and animal pathogens.