The total charge on spike protein, main protease and the two in complex with rutin were neutralized by adding –6, –4, –6, and –4 charges, respectively, and were incorporated into the simulation system by compensating the water molecules in the arbitrary locations inside the simulation box.