The spike protein of SARS-CoV-2 is known to interact
with the human
ACE2 protein via its receptor binding domain (RBD). We have investigated
the molecular nature of this interprotein interaction and the associated
free energy diagrams for the unbinding of the two proteins for SARS-CoV-2
and some of its known variants through all-atom simulations. The present
work involves generation and analysis of 2.5 μs of unbiased
and 4.2 μs of biased molecular dynamics trajectories in total
for five explicitly solvated RBD-ACE2 systems at full atomic level.
First, we have made a comparative analysis of the details of residue-wise
specific interactions of the spike protein with ACE2 for SARS-CoV-1
and SARS-CoV-2. It is found that the average numbers of both direct
interprotein and water-bridged hydrogen bonds between the RBD and
ACE2 are higher for SARS-CoV-2 than SARS-CoV-1. These higher hydrogen
bonded interactions are further aided by enhanced nonspecific electrostatic
attractions between the two protein surfaces for SARS-CoV-2. The free
energy calculations reveal that there is an increase in the free energy
barrier by ∼1.5 kcal/mol for the unbinding of RBD from ACE2
for SARS-CoV-2 compared to that for SARS-CoV-1. Subsequently, we considered
the RBDs of three variants of SARS-CoV-2, namely N501Y, E484Q/L452R,
and N440K. The free energy barrier of protein unbinding for the N501Y
variant is found to be ∼4 kcal/mol higher than the wild type
SARS-CoV-2 which can be attributed to additional specific interactions
involving Tyr501 of RBD and Lys353 and Tyr42 of ACE2 and also enhanced
nonspecific electrostatic interaction between the protein surfaces.
For the other two mutant variants of E484Q/L452R and N440K, the free
energy barrier for protein unbinding increases by ∼2 and ∼1
kcal/mol, respectively, compared with the wild type SARS-CoV-2, which
can be attributed to an increase in the number of interprotein hydrogen
bonds for the former and also to enhanced positive electrostatic potential
on the RBD surfaces for both of the variants. The successive breaking
of interprotein hydrogen bonds along the free energy pathway of the
unbinding process is also found out for all five systems studied here.