Homologous modeling and binding ability analysis of Spike protein after point mutation of severe acute respiratory syndrome coronavirus 2 to receptor proteins and potential antiviral drugs

Abstract

Objective To explore the natural mutations in Spike protein (S protein) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and the changes of affinity between virus and associated receptors or drug molecules before and after the mutation based on whole length sequencing results. Methods In the study, the bioinformatics analysis of all the published sequences of SARS-CoV-2 was conducted and thus the high frequency mutation sites were affirmed. Taking advantages of PolyPhen-2, the functional influence of each mutation in S protein was prospected. The 3D homologous modelling was performed by SWISS-MODEL to establish mutated S protein structural model, in which the protein-docking was then implemented with angiotensin-converting enzyme 2 (ACE2), dipeptidyl peptidase-4 (DPP4) and aminopeptidase N (APN) by ZDOCK, and the combining capacity of each mutated S protein evaluated by FiPD. Finally, the binding ability between mutated S proteins and anti-virus drugs were prospected and evaluated through AutoDock-Chimera 1.14. Results The mutations in specific region of S protein had greater tendency to destroy the S protein function by analysis of mutated S protein structure. Protein-receptor docking analysis between naturally mutated S protein and host receptors showed that, in the case of spontaneous mutation, the binding ability of S protein to ACE2 tended to be weakened, while the binding ability of DPP4 tended to be enhanced, and there was no significant change in the binding ability of APN. According to the computational simulation results of affinity binding between small molecular drugs and S protein, the affinity of aplaviroc with S protein was significantly higher than that of other small molecule drug candidates. Conclusion The region from 400-1 100 amino acid in S protein of SARS-CoV-2 is the mutation sensitive part during natural state, which was more potential to mutate than other part in S protein during natural state. The mutated SARS-CoV-2 might tend to target human cells with DPP4 as a new receptor rather than keep ACE2 as its unique receptor for human infection. At the same time, aplaviroc, which was used for the treatment of human immunodeficiency virus (HIV) infection, may become a new promising treatment for SARS-CoV-2 and could be a potential choice for the development of SARS-CoV-2 drugs.