Abstract
The changes in charge distribution caused by mutations in the spike protein may play a crucial role in balancing infectivity and immune evasion during the evolution of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). To explore how charge increments in spike protein variants influence viral evolution, a statistical analysis was conducted on 57 SARS-CoV-2 variants, examining relationships between charge distribution, lineage divergence, angiotensin-converting enzyme 2 (ACE2) affinity, immune evasion, and receptor-binding domain (RBD) expression. A phylogenetic tree was also reconstructed using only the charge properties of mutation sites. Results indicated that with increasing lineage divergence, overall positive charge initially rose sharply and then more gradually. Partitioning the spike protein into three domains-the RBD, the N-terminal flanking region (B-RBD), and the C-terminal flanking region (A-RBD)-revealed distinct patterns: positive charge increased in the RBD and A-RBD, whereas the B-RBD accumulated negative charge. Charge increments were negatively associated with ACE2 affinity and RBD expression but positively correlated with immune evasion. The k-mer-based tree derived from charge-reduced sequences showed a topology consistent with the whole-genome tree. These findings suggest that charge distribution in spike proteins is closely linked to viral evolution, with the opposing trends in the RBD and B-RBD potentially reflecting a balance between infectivity and immune escape.