Mutational landscape and in silico structure models of SARS-CoV-2 Spike Receptor Binding Domain reveal key molecular determinants for virus-host interaction

AbstractProtein-protein interactions between virus and host are crucial for infection. SARS-CoV-2, the causative agent of COVID-19 pandemic is an RNA virus prone to mutations. Formation of a stable binding interface between the Spike (S) protein Receptor Binding Domain (RBD) of SARS-CoV-2 and Angiotensin-Converting Enzyme 2 (ACE2) of host actuates viral entry. Yet, how this binding interface evolves as virus acquires mutations during pandemic remains elusive. Here, using a high fidelity bioinformatics pipeline, we analysed 31,403 SARS-CoV-2 genomes across the globe, and identified 444 non-synonymous mutations that cause 49 distinct amino acid substitutions in the RBD. Molecular phylogenetic analysis suggested independent emergence of these RBD mutants during pandemic. In silico structure modelling of interfaces induced by mutations on residues which directly engage ACE2 or lie in the near vicinity revealed molecular rearrangements and binding energies unique to each RBD mutant. Comparative structure analysis using binding interface from mouse that prevents SARS-CoV-2 entry uncovered minimal molecular determinants in RBD necessary for the formation of stable interface. We identified that interfacial interaction involving amino acid residues N487 and G496 on either ends of the binding scaffold are indispensable to anchor RBD and are well conserved in all SARS-like corona viruses. All other interactions appear to be required to locally remodel binding interface with varying affinities and thus may decide extent of viral transmission and disease outcome. Together, our findings propose the modalities and variations in RBD-ACE2 interface formation exploited by SARS-CoV-2 for endurance.ImportanceCOVID-19, so far the worst hit pandemic to mankind, started in January 2020 and is still prevailing globally. Our study identified key molecular arrangements in RBD-ACE2 interface that help virus to tolerate mutations and prevail. In addition, RBD mutations identified in this study can serve as a molecular directory for experimental biologists to perform functional validation experiments. The minimal molecular requirements for the formation of RBD-ACE2 interface predicted using in silico structure models may help precisely design neutralizing antibodies, vaccines and therapeutics. Our study also proposes the significance of understanding evolution of protein interfaces during pandemic.