Abstract
The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) outbreak is a public health emergency of international concern. The spike glycoprotein (S protein) of SARS-CoV-2 is a key target of antiviral drugs. Focusing on the existing S protein structure, molecular docking was used in this study to calculate the binding energy and interaction sites between 14 antiviral molecules with different structures and the SARS-CoV-2 S protein, and the potential drug candidates targeting the SARS-CoV-2 S protein were analyzed. Tizoxanide, dolutegravir, bictegravir, and arbidol were found to have high binding energies, and they effectively bind key sites of the S1 and S2 subunits, inhibiting the virus by causing conformational changes in S1 and S2 during the fusion of the S protein with host cells. Based on the interactions among the drug molecules, the S protein and the amino acid environment around the binding sites, rational structure-based optimization was performed using the molecular connection method and bioisosterism strategy to obtain Ti-2, BD-2, and Ar-3, which have much stronger binding ability to the S protein than the original molecules. This study provides valuable clues for identifying S protein inhibitor binding sites and the mechanism of the anti-SARS-CoV-2 effect as well as useful inspiration and help for the discovery and optimization of small molecule S protein inhibitors.