Decoding molecular factors shaping human angiotensin converting enzyme 2 receptor usage by spike glycoprotein in lineage B beta-coronaviruses

Acquiring the human ACE2 receptor usage trait enables the coronaviruses to spill over to humans. However, the origin of the ACE2 usage trait in coronaviruses is poorly understood. Using a multi-disciplinary approach combining evolutionary bioinformatics and molecular dynamics simulation, we decode the principal driving force behind human ACE2 receptor recognition in coronaviruses. Genomic content, evolutionary divergence, and codon usage bias analysis reveal that SARS-CoV2 is evolutionarily divergent from other human ACE2-user CoVs, indicating that SARS-CoV2 originates from a different lineage. Sequence analysis shows that all the human ACE2-user CoVs contain two insertions in the receptor-binding motif (RBM) that directly interact with ACE2. However, the insertion sequences in SARS-CoV2 are divergent from other ACE2-user CoVs, implicating their different recombination origins. The potential of mean force calculations reveals that the high binding affinity of SARS-CoV2 RBD to human ACE2 is primarily attributed to its ability to form a higher number of hydrogen bonds than the other ACE2-user CoVs. The adaptive branch-site random effects likelihood method identifies positive selection bias across the ACE2 user CoVs lineages. Recombination and selection forces shape the spike evolution in human ACE2-using beta-CoVs to optimize the interfacial hydrogen bonds between RBD and ACE2. However, these evolutionary forces work within the constraints of nucleotide composition, ensuring optimum codon adaptation of the spike (S) gene within the host cell.