The value of electrostatic potentials of the spike receptor binding and N-terminal domains in addressing transmissibility and infectivity of SARS-CoV-2 variants of concern.

Dear Editor

Ibrahim and colleagues, in this Journal, published a predicted structure for the SARS-CoV-2 receptor binding site. We and others have focussed attention upon the electrostatic potential (EP) of the receptor binding (EPrbd) and N-terminal (EPntd) domains of SARS-CoV-2 Spike (S protein) that in all previously studied SARS-CoV-2 Variants of Concern (VOC) have a positive charge and strongly interact with negatively-charged patches of the ACE2 and other receptors.2, 3, 4, 5

SARS-CoV-2 Omicron variant is the last and most transmissible of all SARS-CoV-2 mutants. It has an unprecedented, high number of mutations in the receptor binding domain and a novel three amino acids (EPE) insertion in the N-terminal domain of the S Protein, which have determined a new SARS-CoV-2 pandemic trajectory. Initial observations of qualitative nature suggested increasingly positive EPrbd in the successive, worldwide-spreading VOC. We here particularly focus upon VOC which have been able to spread globally and overtake the previously dominant viral strain(s). We assign numerical values to the EPrbd and EPntd, and identify some peculiarities of the Omicron variant that could markedly impact virus transmissibility and pathogenesis.

VOC RBD and NTD, their modeling and energy minimization were as previously described. Net RBD and NTD charge at pH= 7.0 were predicted using the PROPKA3, a software allowing pK calculationa of protein ionizable residues, taking int account the influence of the local structural environment, inclusive of residue solvent exposure, presence of hydrogen bonds or interactions with other charged groups. The pKa values were used to predict the overall domain charge (EPrbd and EPntd). Surface electrostatic potential has been calculated for each variant RBD and NTD by means of the program APBS (Adaptive Poisson-Boltzmann Solver) within the PyMOL environment using default parameters. Electrostatic potential was mapped onto the domain accessibility surfaces and represented with a color gradient ranging from blue to red indicating positive and negative values, respectively.

The calculated pK were used to assign numerical values to the EPrbd and EPntd. Knowing the importance of the RBD-ACE2 binding strength for SARS-CoV-2 transmissibility, we also compared the EPrbd values to the estimated effective reproduction numbers (ERN) of each VOC. , Confirming color gradient visualization, the EPrbd was positive and linearly increased from the Wuhan strain (2.15) to successive variants, reaching its highest value (5.22) in Omicron VOC (Table 1 ). Other noticeable data in Table 1 are, first, that EPrbd values fall within the range of the estimated ERN of the Wuhan reference, Alpha and Delta but not Omicron variant which has a remarkably high ERN value range. The second is the equivalence of the EPrbd fold increase of Omicron relative the Wuhan strain (2.44) to the (around) 2.4 fold enhanced affinity of Omicron RBD binding to ACE2 relative to the Wuhan strain, as measured by surface plasmon resonance. The third is the negative EPntd of the Omicron variant, contrasting the EPntd positivity of all other SARS-CoV-2 strains. All this highlights the peculiarity of the Omicron VOC in such basic chemical and functionals parameters as the RBD and NTD electrostatic potentials.

Table 1

Predicted net charge of the Spike RBD and NTD variants in the folded state calculated with the PROPKA program, and compared with the ranges of the estimated effective reproduction number (ERN)a.

SARS-CoV-2	Pango Lineage	EPrbd	EPntd	Estimated  ERN (range)a
Wuhan	B.1	2.15	1.30	1.4–2.5
Alpha	B.1.1.7	3.18	1.69	1.8–4.7
Delta	B.1.617.2	4.15	1.28	3.4–8.6
Omicron	B.A.1	5.22	-1.10	9.7--13.0

Taken from Refs. 2,7; ERN: effective reproduction number.

Our data suggest that EPrbd could be a proxy of the binding affinity RBD/ACE2, then impacting on virus transmissibility. Hypothetically, the higher the EPrbd, the higher the binding affinity to ACE2, the more transmissible the VOC. This hypothesis appears to hold with the Wuhan reference strain, Alpha and Delta VOC as their EPrbd values do linearly increase and parallel the value range of their effective reproduction numbers. as estimated by epidemiological and clinical studies. , However, it does not hold for the Omicron variant since its EPrbd value is substantially lower than the estimated ERN, This VOC is characterized by 37 mutations in the sequence of the Spike (S protein), of which fifteen are just located in its RBD.7, 8, 9 While they appear to have preserved if not increased (as our EPrbd data and a previous report suggests; 9) the binding strength to ACE2 receptor, there is both laboratory and clinical evidence that this mutational array allows SARS-CoV-2 escaping from neutralizing antibodies of vaccinated and COVID-19 convalescent subjects, as well as the great majority of the monoclonal antibodies in clinical use. , Thus, enhanced transmissibility and immuno-escape are the two features of Omicron worldwide spread.

Mutations in other molecular constituents, both within and outside the S protein could also contribute to virus transmissibility. One is the electrostatic potential of S-protein N-terminal domain . The negative value of the Omicron EPntd is rather unusual, in clear contrast with the positivity of all other variants, as already reported. Rather than affecting RBD-ACE2 interaction, the S protein-NTD could play a role in virus binding to other host cell surface constituent(s). SARS-CoV-2 has been shown to use the NTD domain to recognize and bind to several host receptors or candidate-co-receptors , , If EP positivity is important for SARS-CoV-2 binding to the above receptors, as it is for RBD-ACE2 interaction, then the negative EPntd value would imply a less efficient receptor binding of the Omicron variant relative to that of Wuhan, Alpha and Delta variants. Given that at least one receptor is markedly expressed in lung and bronchial cells , the EPntd negative value could be one of the factors contributing to explain the supposed lower pathogenicity of the Omicron variant for the low respiratory tract.

We recognize the limitations of our in silico data, though they are tempered by the similarity between our ratio EPrbd Omicron variant /Wuhan strain to that measured by Cameroni et al. using a surface plasmon resonance technique. We also recognize the rather speculative nature of some of our suggestions. Experimental confirmations are needed. All above considered, our data highlight some special features of S protein RBD and NTD moieties possessed by the most transmissble and contagious Omicron variant. We invite to consider accurate EP determinations as a help in disentangling the complex mechanisms of SARS-CoV-2 transmissibility and infectivity.

Declaration of Competing Interest

We report no competing interests.