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.

The spike protein of SARS‐CoV‐2 exists in two major conformational states, namely the ‘open’ and ‘closed’ states which are also known as the ‘up’ and ‘down’ states, respectively. In its open state, the receptor binding domain (RBD) of the protein is exposed for binding with ACE2, whereas the spike RBD is inaccessible to ACE2 in the closed state of the protein. In the current work, we have performed all‐atom microsecond simulations of the full‐length trimeric spike protein solvated in explicit aqueous medium with an average system size of ~0.7 million atoms to understand the molecular nature of intra‐ and inter‐chain interactions, water‐bridged interactions between different residues that contribute to the stability of the open and closed states of the protein, and also the free energy landscape for transition between the open and closed states of the protein. We have also examined the changes of such interactions that are associated with switching from one state to the other through both unbiased and biased simulations at all‐atom level with total run length of 4 μs. Interestingly, after about 0.8 μs of unbiased molecular dynamics run of the spike system in the open state, we observed a gradual transition of the monomeric chain (B) from open to its partially closed or down state. Initially the residues at the interface of chain B RBD in the open state spike protein were at non‐hydrogen‐bonding distances from the residues of chain C RBD. However, the two RBDs gradually came closer and finally the residue S459 of the RBD of chain B made a hydrogen bond with F374 of chain C in the last 200 ns of the simulation along with formation of a few more hydrogen bonds involving other residues. Since no transition from closed to the open state of the protein is observed in the present 1 μs unbiased simulation of the closed state protein, the current study seems to suggest that the closed conformational state is preferred for the spike protein of SARS‐CoV‐2 in aqueous medium. Furthermore, calculations of the free energy surface of the conformational transition from open (up) to the closed (down) state using a biased simulation method reveal a free energy barrier of ~3.20 kcal/mol for the transition of RBD from open to the closed state, whereas the barrier for the reverse process is found to be significantly higher.

No abstract