Next Article in Journal
Gal-3BP in Viral Infections: An Emerging Role in Severe Acute Respiratory Syndrome Coronavirus 2
Next Article in Special Issue
Theoretical Explanation for the Rarity of Antibody-Dependent Enhancement of Infection (ADE) in COVID-19
Previous Article in Journal
Inactivation of the Lateral Hypothalamus Attenuates Methamphetamine-Induced Conditioned Place Preference through Regulation of Kcnq3 Expression
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 23
Issue 13
10.3390/ijms23137304
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Electrostatic Map of the SARS-CoV-2 Virion Specifies Binding Sites of the Antiviral Cationic Photosensitizer

by
Vladimir Fedorov
1,*,
Ekaterina Kholina
1,
Sergei Khruschev
1,
Ilya Kovalenko
1,2,
Andrew Rubin
1 and
Marina Strakhovskaya
1,2
1
Faculty of Biology, Lomonosov Moscow State University, 119234 Moscow, Russia
2
Federal Scientific and Clinical Center of Specialized Types of Medical Care and Medical Technologies, Federal Medical and Biological Agency of Russia, 115682 Moscow, Russia
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2022, 23(13), 7304; https://doi.org/10.3390/ijms23137304
Submission received: 4 June 2022 / Revised: 26 June 2022 / Accepted: 28 June 2022 / Published: 30 June 2022
(This article belongs to the Special Issue State-of-the-Art Molecular Biophysics in Russia)
Browse Figures
Versions Notes

Abstract

:
Electrostatics is an important part of virus life. Understanding the detailed distribution of charges over the surface of a virus is important to predict its interactions with host cells, antibodies, drugs, and different materials. Using a coarse-grained model of the entire viral envelope developed by D. Korkin and S.-J. Marrink’s scientific groups, we created an electrostatic map of the external surface of SARS-CoV-2 and found a highly heterogeneous distribution of the electrostatic potential field of the viral envelope. Numerous negative patches originate mainly from negatively charged lipid domains in the viral membrane and negatively charged areas on the “stalks” of the spike (S) proteins. Membrane (M) and envelope (E) proteins with the total positive charge tend to colocalize with the negatively charged lipids. In the E protein pentamer exposed to the outer surface, negatively charged glutamate residues and surrounding lipids form a negative electrostatic potential ring around the channel entrance. We simulated the interaction of the antiviral octacationic photosensitizer octakis(cholinyl)zinc phthalocyanine with the surface structures of the entire model virion using the Brownian dynamics computational method implemented in ProKSim software (version r661). All mentioned negatively charged envelope components attracted the photosensitizer molecules and are thus potential targets for reactive oxygen generated in photosensitized reactions.

1. Introduction

In the past two decades, the attention of researchers from different fields of science has been steadily drawn to coronaviruses (CoVs), enveloped viruses with single-stranded, positive-sense RNA genomes that infect birds and mammals [1,2]. CoVs are of particular concern due to their zoonotic potential and high plasticity, ability to infect multiple hosts, quick adaptation to new hosts and changing environmental conditions [3]. This is facilitated by large genomes, jumping mutations, and recombination potential.
In SARS-CoV-2, RNA encodes 16 non-structural proteins, nine accessory proteins, and four major structural proteins [4]—the nucleocapsid (N) protein, the envelope (E) protein, the membrane (M) protein, and the spike (S) protein. The N protein binds to viral RNA to form an RNP assembly [5] and determines its packaging in the ~80 nm-diameter lumen [6]. Through interactions with the M protein, N incorporates viral RNA in the budding virions independently of other proteins [7].
When CoVs bud in the endoplasmic-reticulum Golgi intermediate compartment (ERGIC), the nucleocapsid is covered by a lipid bilayer, derived from host lipids with embedded transmembrane viral proteins (E, M, and S) [8]. It is known that ER membranes are enriched in phosphatidylcholine (PC) and phosphatidylinositol (PI) [9], and, to mimic ERGIC, the lipid bilayer should contain PC, PI, phosphatidylethanolamine (PE), phosphatidylserine (PS), and cholesterol (Chol) [10]. SARS-CoV-2 and other CoVs use the small (76–109 amino acids) E protein to take control of ERGIC [11,12], but in the viral envelope, the E protein is represented by only a few copies. E is an integral membrane protein, with a single hydrophobic domain [12] that forms a homopentameric channel, with mild cation selectivity in negatively charged lipid membranes, termed “viraporin” [13].
The most abundant M protein is responsible for sorting and the incorporation of viral components in the budding virions and maintaining the shape and size of CoVs [4,7,14]. The M protein (about 230 amino acids) has three transmembrane domains and usually exists in the form of functional dimers [14,15]. It forms the matrix of the envelope and interacts with all other structural proteins.
The highly glycosylated S proteins (each about 1300 amino acids) form homotrimeric spikes, which protrude for about 20 nm from the viral membrane and make up the “corona” appearance. Each virus particle of SARS-CoV-2 bears 26 ± 15 spikes [6]. Host cell proteases, such as furin, cleave the S protein into S1 and S2 subunits in the infected cell, and thus in some CoVs the S1 and S2 subunits are associated non-covalently. The receptor-binding domain (RBD), responsible for coronavirus binding to a host cell receptor, is located in the S1 subunit of the spike protein, while the fusion peptide is assigned to the S2 subunit. The S2 single transmembrane domain anchors the S protein in the envelope [16].
Biophysics and structural biology provide key information for drug design to fight pandemic viruses [17,18]. However, the presence of lipid bilayers creates difficulties for the complete reconstruction of enveloped viruses by traditional methods of structural biology [19]. The complexity of the viral supramolecular structures, and the interactions between viral components and active substances, can be successfully described using so-called “computational microscopy”, with a spatiotemporal resolution unmatched by other methods [20].
Among the advances in computational virology bringing us closer to creating realistic models of whole virions are: the study of the structural dynamics of viral capsids, by employing a coarse-graining (CG) molecular dynamics method [21]; atomic HIV-1 capsid model [22]; simulation of Herpes simplex virus type 2 B-capsid and chromatin models in GENeralized-Ensemble Simulation System (GENESIS) [23]; a dynamic and integrative computational model of an influenza A virion [24]; computational model of the envelope of the dengue virion [25]; Martini CG models of intact SARS-CoV and SARS-CoV-2 envelopes [26]; and multiscale CG models of the SARS-CoV-2 virion [27,28].
In this study, we utilized the CG model of SARS-CoV-2 developed in [28], and kindly provided by W. Pezeshkian, D. Korkin and S.J. Marrink, to simulate the interaction of the photosensitizer (PS) octakis(cholinyl)zinc phthalocyanine (ZnPcChol8+) with the surface structures of the whole virion. We have recently shown that ZnPcChol8+ is a potent antiviral PS against bovine CoV [29], H5N8 avian influenza virus [30], and SARS-CoV-2 [31] in vitro. Proteins and unsaturated lipids of the viral envelope are readily oxidized with reactive oxygen species and free radicals generated by photoactivated PS in photodynamic reactions. As the photodynamic impact increased, the successive morphological changes in the H5N8 avian influenza virus and bovine CoV were as follows: the disappearance of the spikes, the change in the shape and size of the viral particles, the destruction of the envelope, and the complete collapse of the viral structure [29,30]. We analyzed electrostatic potential and simulated the binding of ZnPcChol8+ to the S protein of SARS-CoV, MERS-CoV, and SARS-CoV-2 [31,32] by means of Brownian dynamics (BD). In all three studied S proteins, we found a major binding site for ZnPcChol8+ at the connection of the S protein stalk and the head, adjacent to the heptad repeat 2 (HR2) domain. Hence, we created the electrostatic map of the whole SARS-CoV-2 virion and extended PS interaction simulations to the entire viral envelope to create a model which predicts PS binding sites and relevant targets of the reactive oxygen.

2. Results

2.1. Electrostatic Potential Distribution of the SARS-CoV-2 Envelope

The lipid and protein composition of SARS-CoV-2’s viral envelope is represented in Figure 1a. The CG model envelope includes three protein types, namely M, E, and S proteins. The M protein dominates and is presented in 1003 copies, the S protein is the only one protruding from the viral membrane and is not so numerous, it presented in 25 copies. In contrast to the other proteins, the viral membrane contains only two E proteins. The lipids of the viral membrane are presented by neutral palmitoyl oleoyl phosphatidylethanolamine (POPE, 20%), palmitoyl oleoyl phosphatidylcholine (POPC, 59%), cholesterol (CHOL, 4.5%) molecules and by negatively charged palmitoyl oleoyl phosphatidylinositol (POPI, 10%), palmitoyl oleoyl phosphatidylserine (POPS, 2%) and cardiolipin (CDL2, 4.5%).
The SARS-CoV-2 virus envelope possesses a positive net charge of 8856 elementary charges, mainly due to the presence of more than a thousand positively charged M protein molecules (charge of each is +22). However, the SARS-CoV-2 virus envelope demonstrates a highly heterogeneous electrostatic potential field, with large alternating areas of positive and negative potential (Figure 1b). Negative patches of electrostatic potential are generated mainly by POPI (−1), POPS (−1), and CDL2 (−2) molecules and some negatively charged amino acids of S protein, as described in our previous papers [31,32]. Interestingly, negatively charged lipid molecules are prone to forming domains in the viral membrane (Figure 1 and Figure 2a), resulting in the existence of vast negatively charged areas on the surface of the SARS-CoV-2 virus envelope.

2.2. Interactions of Photosensitizer with the SARS-CoV-2 Envelope Revealed by Brownian Dynamics

To investigate which components of SARS-CoV-2 viral membrane attract ZnPcChol8+ molecule, we performed 40 thousand independent BD simulations of ZnPcChol8+ molecule relative to immobile viral envelope. Panels C and D of Figure 1 demonstrate the resulting 40,000 different electrostatic encounter complexes of photosensitizer with viral envelope, with electrostatic attraction energy exceeding 8 kT. Table 1 summarizes the statistics of contacts of the ZnPcChol8+ molecule with different components of the viral envelope.
About 43% of the ZnPcChol8+ molecules were found in electrostatic encounter complexes with SARS-CoV-2 membrane proteins. Among them, about 80% formed electrostatic contacts with S proteins; about 19%—with M proteins; and less than 1%—with E proteins. Most of the ZnPcChol8+ molecules bound to S proteins (91%) did not contact any membrane lipids or other proteins. On the contrary, almost all ZnPcChol8+ molecules bound to M proteins also interacted with negatively charged lipids (POPI, POPS, and CDL2); 70% of these molecules interacted with at least two charged lipid molecules. As for E proteins, about one-third of ZnPcChol8+ molecules bound to these proteins also formed contacts with negatively charged lipids, namely POPI (about 60%) and CDL2 (about 40%).
The remaining 57% of the ZnPcChol8+ molecules had no close contacts with membrane proteins but were bound to negatively charged lipids only; 63% of these molecules simultaneously interacted with at least two charged lipid molecules. Among the ZnPcChol8+ molecules interacting exclusively with charged lipid molecules (see Table 2), the most typical combination was the contact with one POPI molecule and one CDL2 molecule (19%). In 9% of cases, the ZnPcChol8+ molecule interacted with two POPI molecules, and in 7% of cases—with two CDL2 molecules. Another 7% of cases interacted with two CDL2 molecules and one POPI molecule.
From Table 2 it follows that in most cases, when ZnPcChol8+ interacts exclusively with lipids, the PS molecule forms contacts with several lipid molecules at once, due to the presence of patches of negatively charged lipids in the viral membrane. These patches are formed by either the same or different types of lipids (Figure 2a). The colocalization of negatively charged lipids POPI, POPS, and CDL2 results in formation of large regions of negative electrostatic potential (Figure 2b). As one can expect, these areas intensively attract ZnPcChol8+ molecules (Figure 2c,d).
On their own, the most abundant in the viral envelope are M proteins, with a total positive charge fail to attract PS molecules. However, M proteins due to their positive charge are colocalized with negatively charged lipids (Figure 2a), which attract ZnPcChol8+ molecules (Figure 2c). Therefore, in 19% of electrostatic complexes we observed contacts of PS with M proteins.
In the pentameric ion channel, each monomer, that is, the E protein, comprises three domains: the N-terminal, transmembrane and C-terminal domains. The N-terminal domain exposed outside the virion contains two negatively charged glutamate residues, E7 and E8. Despite the total positive charge (+10) of the E pentamer, these glutamates generate a ring area of negative electrostatic potential around the channel entrance (Figure 3a). In addition, the E protein is surrounded by negatively charged lipids, mainly CDL2, which further contribute to the negative surface potential around the channel. The resulting vast area of negative electrostatic potential intensively attracts ZnPcChol8+ molecules (Figure 3b).
Figure 4a demonstrates the section of the viral envelope with one of the 25 glycosylated S trimers embedded in the viral membrane. The S protein protruding from the viral membrane is one of the electrostatically attractive targets for ZnPcChol8+ [31,32]. As shown in these papers by means of all-atom molecular modeling, the electrostatic field is heterogeneously distributed on the surface of the “head” of the S protein. The “stalk” of the S protein is highly negatively charged and provides the major binding site to ZnPcChol8+ molecules (Figure 4c). In [31,32], we investigated the binding ability of non-glycosylated S protein. In the current CG model of SARS-CoV-2 envelope, the S protein trimer is glycosylated with 579 residues, among which 12 are negatively charged, that may influence ZnPcChol8+ binding. In accordance with our previous results, the major binding sites for ZnPcChol8+ are located in the linkers between FP and HR1, HR1 and HR2, as well as in the HR2 domain (Figure 4b,d). The presence of glycans on the surface of S proteins has little effect on ZnPcChol8+ binding. However, we found additional minor binding sites of ZnPcChol8+ on the surface of NTD and RBD domains of the S protein (Figure 4d).

3. Discussion

Electrostatic forces play a prominent role in the life cycle of coronaviruses, primarily by participating in the binding of the S protein to the host cell [33,34,35]. Mutations affecting the surface electrostatic potential may ultimately alter the infectivity, pathogenicity, and transmission rate of the virus [34].
In a polar aquatic environment, viruses have a pH-dependent surface charge due to the protonation/deprotonation of carboxyl and amino groups, and can be characterized by the so-called isoelectric point (IEP). IEP is considered a reliable quantitative characteristic that is useful in comparing the interactions of different viruses with the environment under different experimental conditions [36]. Michen and Graule [37] analyzed 152 IEP values from 104 viruses and found IEP to be in the range from 1.9 to 8.4, most often from 3.5 to 7, which may indicate an excess of acidic amino acids (aspartic acid and/or glutamic acid) exposed on the viral surface. The calculated IEP of SARS-CoV-2 ranges from 5.2 to 6.2, and when measured by means of chemical force microscopy in 20 mM salt is 5.2–5.3 [38]. The calculated IEP of the structural envelope S, E, and M proteins of SARS-CoV-2 are 6.24, 8.57, and 9.51, respectively [39]. In addition to proteins, phospholipids of the viral membrane, which originate from the membranes of the host cell, contribute to the total surface charge of enveloped viruses. Their diversity and variable ratio affect the total charge of virions, and complex lipid profiles cannot be predicted from the viral genome [36].
Another useful characteristic that reflects the total charge of a viral particle, its adsorption properties, and stability in dispersion is the zeta potential. The average zeta potential of the porcine coronavirus dispersed in double-distilled water was determined to be −25.675 mV [40]. However, zeta potential measurements require concentrated viral suspensions and are not applicable to intact pathogenic viruses.
To predict the precise binding of antiviral compounds, it is important to study the local electric charge distribution of the viral envelope components exposed to the surface. As we demonstrated, domains of S and M proteins protruding from the viral membrane have highly heterogeneous potential, while the exposed outside surface of the E protein is almost negatively charged. The total charge of proteins in the current model viral envelope is +21,261. The viral envelope is composed of a significant fraction of negatively charged lipids, contributing 12,405 elementary negative charges to the total +8856 net virion charge. As we demonstrated, negatively charged lipids form patches of condensed negative potential and provide more than a half of ZnPcChol8+ binding events. About 40% of ZnPcChol8+ binding events occur due to direct contact with envelope proteins, preferentially with S proteins. A high level of the S protein glycosylation almost does not influence the binding of polycationic PS. Glycan residues provide no steric barrier and leave enough space for PS binding. Additionally, predominantly neutral glycans do not shield the S protein surface charges, and so do not hinder electrostatic interactions with PS molecules.
The M protein with a total positive charge +22, the most abundant in the viral envelope, almost does not bind ZnPcChol8+ molecules. However, some binding events do happen due to the existence of local negative areas on the M protein surface. Despite the slight positive charge of the E protein, the N-terminal domain contains two glutamates, E7 and E8. These residues are exposed to the outside envelope surface and form a ring-like area of negative potential, contouring the entrance into the pentameric channel, attracting ZnPcChol8+ molecules.
We have recently shown that ZnPcChol8+ is a potent antiviral PS against bovine CoV [29], H5N8 avian influenza virus [30], and SARS-CoV-2 [31] in vitro. Using transmission electron microscopy with negative contrast, we have identified [29,30] a number of morphological changes common to the avian influenza virus and bovine coronavirus that occur during photodynamic inactivation: the disappearance of the spikes, the change in the shape and size of the viral particles, the destruction of the envelope, and the complete collapse of the viral structure. All of these lesions resulted in a loss of virus infectivity, which was determined by virus titration. In accordance with these morphological changes, we found that the most preferable binding site of PS is located on the S protein. Taking into account the small diffusion radius of singlet oxygen in the biological environment, the observed photosensitized damage to the viral membranes is in favor of close contacts of the PS with lipid membrane components. Actually, the negatively charged lipids in the viral model envelope provide multiple binding sites for PS molecules. Thus, the study of the detailed electrostatic map of the whole virion using a computer model provides unique opportunities to reveal the binding sites of charged molecules on the surface of the virus.

4. Materials and Methods

The interaction of PS with the SARS-CoV-2 envelope, both represented at the CG level, was studied by BD. The CG approach has the advantage of dealing with simplified molecular models, compared to the all-atom level. We used Martini CG models, based on the assumption that groups of several heavy atoms of particular chemical groups could be combined to represent one CG particle. The CG Martini force field proved to be an effective approach for studying different biomolecular systems [41].
In this study, we used the CG model based on Martini 3 force field of SARS-CoV-2 envelope developed in [28], and kindly provided to us by W. Pezeshkian, D. Korkin and S.-J. Marrink (8 February 2022 version). The lipids of the viral membrane were presented by 11817 POPE, 34860 POPC, 2658 CHOL, 5908 POPI, 1181 POPS, and 2658 CDL2 molecules. In total, 2 E proteins, 25 S proteins, and 1003 M proteins were embedded into the viral membrane. The CG model of ZnPcChol8+ was previously developed in our laboratory [42].
The electrostatic potential field and BD calculations were conducted using ProKSim (Protein Kinetics Simulator) software [43,44]. The Poisson–Boltzmann equation was used to calculate electrostatic potential field of SARS-CoV-2 virion and ZnPcChol8+ molecule. SARS-CoV-2 virion envelope was represented as a low dielectric area (ε = 2) with spatially fixed partial charges. The water solvent (dielectric constant ε = 80) with ions was described implicitly by the ionic strength 100 mM. Electrostatic cutoff radius was 3.5 nm.
In BD simulations, the SARS-CoV-2 virion was placed in the middle of a virtual reaction volume with mirror boundary conditions and dimensions of 200 × 200 × 200 nm. In each simulation, a ZnPcChol8+ molecule was randomly placed in the volume not occupied by the virion. The virion was immobile while the ZnPcChol8+ molecule was moving under the action of random and electrostatic forces. The calculation continued until the energy of electrostatic attraction of ZnPcChol8+ to the virion reached 8 kT. The final ZnPcChol8+ position with respect to the virion was saved for further analysis, and then the simulation was restarted from another random position of ZnPcChol8+. In total, we obtained 40,000 configurations of encounter complexes of ZnPcChol8+ with the SARS-CoV-2 virion.
We considered that contact between ZnPcChol8+ molecule and a specific component of the virion was formed if any CG bead of PS approached any CG bead of the virion within a distance of 0.6 nm. To characterize the interaction of PS with virion components throughout the obtained ensemble of structures, we used a custom Python script utilizing the MDAnalysis library [45], which identifies and counts particular components of viral envelope which are in contact with the PS molecule. Contact probability for particular amino acid residues of the S protein was calculated as the average fraction of contacts in PS electrostatic complexes with all 25 S protein molecules over the entire ensemble of structures. The contact probability for every amino acid residue was visualized on the S protein primary sequence by color. PyMOL software [46] was used to visualize molecular structures.

5. Conclusions

In conclusion, using a coarse-grained model of the SARS-CoV-2 envelope, we identified multiple heterogeneously distributed patches of negative potential at the envelope surface structures. These sites, formed mainly by negatively charged membrane phospholipids, as well as exposed exterior negatively charged amino acid residues on the S proteins, attract positively charged PS molecules and are consistent with the previously observed spike loss and membrane destruction, as a result of the photodynamic inactivation of the coronavirus with the same PS. The obtained results of model experiments explain how positively charged antivirals can bind to the coronavirus envelope, despite its net positive charge.

Author Contributions

Conceptualization, I.K., A.R. and M.S.; methodology, V.F., S.K. and I.K.; investigation, V.F., E.K. and S.K.; writing—original draft preparation, all authors; writing—review and editing, all authors; supervision, A.R. and M.S.; project administration, M.S.; funding acquisition, M.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by RFBR (Russian Foundation for Basic Research), project number 20-04-60084.

Institutional Review Board Statement

No ethical approval was required for this study.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We would like to thank W. Pezeshkian, D. Korkin and S.-J. Marrink for providing the CG model (8 February 2022 version) of the SARS-CoV-2 envelope.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Zmasek, C.M.; Lefkowitz, E.J.; Niewiadomska, A.; Scheuermann, R.H. Genomic evolution of the Coronaviridae family. Virology 2022, 570, 123–133. [Google Scholar] [CrossRef] [PubMed]
  2. Mousavizadeh, L.; Ghasemi, S. Genotype and phenotype of COVID-19: Their roles in pathogenesis. J. Microbiol. Immunol. Infect. 2021, 54, 159–163. [Google Scholar] [CrossRef] [PubMed]
  3. Ghai, R.R.; Carpenter, A.; Liew, A.Y.; Martin, K.B.; Herring, M.K.; Gerber, S.I.; Hall, A.J.; Sleeman, J.M.; VonDobschuetz, S.; Behravesh, C.B. Animal reservoirs and hosts for emerging alphacoronaviruses and betacoronaviruses. Emerg. Infect. Dis. 2021, 27, 1015–1022. [Google Scholar] [CrossRef] [PubMed]
  4. Yoshimoto, F.K. The proteins of Severe Acute Respiratory Syndrome Coronavirus-2 (SARS CoV-2 or n-COV19), the cause of COVID-19. Protein J. 2020, 39, 198–216. [Google Scholar] [CrossRef] [PubMed]
  5. Matsuo, T. Viewing SARS-CoV-2 nucleocapsid protein in terms of molecular flexibility. Biology 2021, 10, 454. [Google Scholar] [CrossRef]
  6. Yao, H.; Song, Y.; Chen, Y.; Wu, N.; Xu, J.; Sun, C.; Zhang, J.; Weng, T.; Zhang, Z.; Wu, Z.; et al. Molecular architecture of the SARS-CoV-2 Virus. Cell 2020, 183, 730–738.e13. [Google Scholar] [CrossRef]
  7. Siu, Y.L.; Teoh, K.T.; Lo, J.; Chan, C.M.; Kien, F.; Escriou, N.; Tsao, S.W.; Nicholls, J.M.; Altmeyer, R.; Peiris, J.S.M.; et al. The M, E, and N structural proteins of the severe acute respiratory syndrome coronavirus are required for efficient assembly, trafficking, and release of virus-like particles. J. Virol. 2008, 82, 11318–11330. [Google Scholar] [CrossRef] [Green Version]
  8. Klein, S.; Cortese, M.; Winter, S.L.; Wachsmuth-Melm, M.; Neufeldt, C.J.; Cerikan, B.; Stanifer, M.L.; Boulant, S.; Bartenschlager, R.; Chlanda, P. SARS-CoV-2 structure and replication characterized by in situ cryo-electron tomography. Nat. Commun. 2020, 11, 5885. [Google Scholar] [CrossRef]
  9. Van Meer, G.; de Kroon, A.I. Lipid map of the mammalian cell. J. Cell Sci. 2011, 124 Pt 1, 5–8. [Google Scholar] [CrossRef] [Green Version]
  10. Mandala, V.S.; McKay, M.J.; Shcherbakov, A.A.; Dregni, A.J.; Kolocouris, A.; Hong, M. Structure and drug binding of the SARS-CoV-2 envelope protein transmembrane domain in lipid bilayers. Nat. Struct. Mol. Biol. 2020, 27, 1202–1208. [Google Scholar] [CrossRef]
  11. Barrantes, F.J. Structural biology of coronavirus ion channels. Acta Crystallogr. D Struct. Biol. 2021, 77 Pt 4, 391–402. [Google Scholar] [CrossRef] [PubMed]
  12. Ruch, T.R.; Machamer, C.E. The coronavirus E protein: Assembly and beyond. Viruses 2012, 4, 363–382. [Google Scholar] [CrossRef] [Green Version]
  13. Verdiá-Báguena, C.; Nieto-Torres, J.L.; Alcaraz, A.; DeDiego, M.L.; Torres, J.; Aguilella, V.M.; Enjuanes, L. Coronavirus E protein forms ion channels with functionally and structurally-involved membrane lipids. Virology 2012, 432, 485–494. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Neuman, B.W.; Kiss, G.; Kunding, A.H.; Bhella, D.; Baksh, M.F.; Connelly, S.; Droese, B.; Klaus, J.P.; Makino, S.; Sawicki, S.G.; et al. A structural analysis of M protein in coronavirus assembly and morphology. J. Struct. Biol. 2011, 174, 11–22. [Google Scholar] [CrossRef] [PubMed]
  15. Cao, Y.; Yang, R.; Wang, W.; Jiang, S.; Yang, C.; Liu, N.; Dai, H.; Lee, I.; Meng, X.; Yuan, Z. Probing the formation, structure and free energy relationships of M protein dimers of SARS-CoV-2. Comput. Struct. Biotechnol. J. 2022, 20, 573–582. [Google Scholar] [CrossRef]
  16. Jackson, C.B.; Farzan, M.; Chen, B.; Choe, H. Mechanisms of SARS-CoV-2 entry into cells. Nat. Rev. Mol. Cell Biol. 2022, 23, 3–20. [Google Scholar] [CrossRef]
  17. Barrantes, F.J. The contribution of biophysics and structural biology to current advances in COVID-19. Annu. Rev. Biophys. 2021, 50, 493–523. [Google Scholar] [CrossRef]
  18. Bárcena, M.; Barnes, C.O.; Beck, M.; Bjorkman, P.J.; Canard, B.; Gao, G.F.; Gao, Y.; Hilgenfeld, R.; Hummer, G.; Patwardhan, A.; et al. Structural biology in the fight against COVID-19. Nat. Struct. Mol. Biol. 2021, 28, 2–7. [Google Scholar] [CrossRef] [PubMed]
  19. Li, S. Cryo-electron tomography of enveloped viruses. Trends Biochem. Sci. 2022, 47, 173–186. [Google Scholar] [CrossRef]
  20. Ingólfsson, H.I.; Arnarez, C.; Periole, X.; Marrink, S.J. Computational ‘microscopy’ of cellular membranes. J. Cell Sci. 2016, 129, 257–268. [Google Scholar] [CrossRef] [Green Version]
  21. Arkhipov, A.; Freddolino, P.L.; Schulten, K. Stability and dynamics of virus capsids described by coarse-grained modeling. Structure 2006, 14, 1767–1777. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Zhao, G.; Perilla, J.R.; Yufenyuy, E.L.; Meng, X.; Chen, B.; Ning, J.; Ahn, J.; Gronenborn, A.M.; Schulten, K.; Aiken, C.; et al. Mature HIV-1 capsid structure by cryo-electron microscopy and all-atom molecular dynamics. Nature 2013, 497, 643–646. [Google Scholar] [CrossRef] [PubMed]
  23. Tan, C.; Jung, J.; Kobayashi, C.; Torre, D.U.; Takada, S.; Sugita, Y. Implementation of residue-level coarse-grained models in GENESIS for large-scale molecular dynamics simulations. PLoS Comput. Biol. 2022, 18, e1009578. [Google Scholar] [CrossRef] [PubMed]
  24. Reddy, T.; Shorthouse, D.; Parton, D.L.; Jefferys, E.; Fowler, P.W.; Chavent, M.; Baaden, M.; Sansom, M.S.P. Nothing to sneeze at: A dynamic and integrative computational model of an influenza A virion. Structure 2015, 23, 584–597. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Reddy, T.; Sansom, M.S. The role of the membrane in the structure and biophysical robustness of the dengue virion envelope. Structure 2016, 24, 375–382. [Google Scholar] [CrossRef] [Green Version]
  26. Wang, B.; Zhong, C.; Tieleman, D.P. Supramolecular organization of SARS-CoV and SARS-CoV-2 virions revealed by coarse-grained models of intact virus envelopes. J. Chem. Inf. Model. 2022, 62, 176–186. [Google Scholar] [CrossRef]
  27. Yu, A.; Pak, A.J.; He, P.; Monje-Galvan, V.; Casalino, L.; Gaieb, Z.; Dommer, A.C.; Amaro, R.E.; Voth, G.A. A multiscale coarse-grained model of the SARS-CoV-2 virion. Biophys. J. 2021, 120, 1097–1104. [Google Scholar] [CrossRef]
  28. Pezeshkian, W.; Grünewald, F.; Narykov, O.; Lu, S.; Wassenaar, T.A.; Marrink, S.J.; Korkin, D. Molecular architecture of SARS-CoV-2 envelope by integrative modeling. bioRxiv 2021. [Google Scholar] [CrossRef]
  29. Zhukhovitsky, V.; Shevlyagina, N.; Zubasheva, M.; Russu, L.; Gushchin, V.; Meerovich, G.; Strakhovskaya, M. Infectivity and morphology of bovine coronavirus inactivated in vitro by cationic photosensitizers. Viruses 2022, 14, 1053. [Google Scholar] [CrossRef] [PubMed]
  30. Korneev, D.; Kurskaya, O.; Sharshov, K.; Eastwood, J.; Strakhovskaya, M. Ultrastructural aspects of photodynamic inactivation of highly pathogenic avian H5N8 influenza virus. Viruses 2019, 11, 955. [Google Scholar] [CrossRef] [Green Version]
  31. Sharshov, K.; Solomatina, M.; Kurskaya, O.; Kovalenko, I.; Kholina, E.; Fedorov, V.; Meerovich, G.; Rubin, A.; Strakhovskaya, M. The photosensitizer octakis(cholinyl)zinc phthalocyanine with ability to bind to a model spike protein leads to a loss of SARS-CoV-2 infectivity in vitro when exposed to far-red LED. Viruses 2021, 13, 643. [Google Scholar] [CrossRef]
  32. Fedorov, V.; Kholina, E.; Khruschev, S.; Kovalenko, I.; Rubin, A.; Strakhovskaya, M. What binds cationic photosensitizers better: Brownian dynamics reveals key interaction sites on spike proteins of SARS-CoV, MERS-CoV, and SARS-CoV-2. Viruses 2021, 13, 1615. [Google Scholar] [CrossRef] [PubMed]
  33. Xie, Y.; Du, D.; Karki, C.B.; Guo, W.; Lopez-Hernandez, A.E.; Sun, S.; Juarez, B.Y.; Li, H.; Wang, J.; Li, L. Revealing the mechanism of SARS-CoV-2 spike protein binding with ACE2. Comput. Sci. Eng. 2020, 22, 21–29. [Google Scholar] [CrossRef] [PubMed]
  34. Pascarella, S.; Ciccozzi, M.; Zella, D.; Bianchi, M.; Benedetti, F.; Benvenuto, D.; Broccolo, F.; Cauda, R.; Caruso, A.; Angeletti, S.; et al. SARS-CoV-2 B.1.617 Indian variants: Are electrostatic potential changes responsible for a higher transmission rate? J. Med. Virol. 2021, 93, 6551–6556. [Google Scholar] [CrossRef] [PubMed]
  35. Pascarella, S.; Ciccozzi, M.; Bianchi, M.; Benvenuto, D.; Cauda, R.; Cassone, A. 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. J. Infect. 2022, 84, e62–e63. [Google Scholar] [CrossRef] [PubMed]
  36. Heffron, J.; Mayer, B.K. Virus isoelectric point estimation: Theories and methods. Appl. Environ. Microbiol. 2021, 87, e02319-20. [Google Scholar] [CrossRef]
  37. Michen, B.; Graule, T. Isoelectric points of viruses. J. Appl. Microbiol. 2010, 109, 388–397. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Areo, O.; Joshi, P.U.; Obrenovich, M.; Tayahi, M.; Heldt, C.L. Single-particle characterization of SARS-CoV-2 isoelectric point and comparison to variants of interest. Microorganisms 2021, 9, 1606. [Google Scholar] [CrossRef]
  39. Scheller, C.; Krebs, F.; Minkner, R.; Astner, I.; Gil-Moles, M.; Wätzig, H. Physicochemical properties of SARS-CoV-2 for drug targeting, virus inactivation and attenuation, vaccine formulation and quality control. Electrophoresis 2020, 41, 1137–1151. [Google Scholar] [CrossRef]
  40. Ghatak, S.; Khona, D.K.; Sen, A.; Huang, K.; Jagdale, G.; Singh, K.; Gopalakrishnan, V.; Cornetta, K.G.; Roy, S.; Khanna, S.; et al. Electroceutical fabric lowers zeta potential and eradicates coronavirus infectivity upon contact. Sci. Rep. 2021, 11, 21723. [Google Scholar] [CrossRef]
  41. Marrink, S.J.; Tieleman, D.P. Perspective on the Martini model. Chem. Soc. Rev. 2013, 42, 6801–6822. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Orekhov, P.S.; Kholina, E.G.; Bozdaganyan, M.E.; Nesterenko, A.M.; Kovalenko, I.B.; Strakhovskaya, M.G. Molecular mechanism of uptake of cationic photoantimicrobial phthalocyanine across bacterial membranes revealed by molecular dynamics simulations. J. Phys. Chem. B 2018, 122, 3711–3722. [Google Scholar] [CrossRef] [PubMed]
  43. Kovalenko, I.B.; Khrushchev, S.S.; Fedorov, V.A.; Riznichenko, G.Y.; Rubin, A.B. The role of electrostatic interactions in the process of diffusional encounter and docking of electron transport proteins. Dokl. Biochem. Biophys. 2016, 468, 183–186. [Google Scholar] [CrossRef]
  44. Fedorov, V.A.; Kovalenko, I.B.; Khruschev, S.S.; Ustinin, D.M.; Antal, T.K.; Riznichenko, G.Y.; Rubin, A.B. Comparative analysis of plastocyanin- cytochrome f complex formation in higher plants, green algae and cyanobacteria. Physiol. Plant. 2019, 166, 320–335. [Google Scholar] [CrossRef] [PubMed]
  45. Gowers, R.J.; Linke, M.; Barnoud, J.; Reddy, T.J.E.; Melo, M.N.; Seyler, S.L.; Dotson, D.L.; Domanski, J.; Buchoux, S.; Kenney, I.M.; et al. MDAnalysis: A Python package for the rapid analysis of molecular dynamics simulations. In Proceedings of the 15th Python in Science Conference, Austin, TX, USA, 11–17 July 2016; Benthall, S., Rostrup, S., Eds.; SciPy: Austin, TX, USA, 2016; pp. 98–105. [Google Scholar] [CrossRef] [Green Version]
  46. Schrödinger, L.L.C.; DeLano, W. The PyMOL Molecular Graphics System, Version 2.5. Available online: https://pymol.org/ (accessed on 3 June 2022).
Ijms 23 07304 g001 550
Figure 1. Three-dimensional structure of CG molecular model and electrostatic surface potential of SARS-CoV-2 virus envelope. CG beads colored in accordance with envelope component types (a). Rectangular boxes pick up the areas shown in Figure 2, Figure 3 and Figure 4 (green, yellow and blue, respectively). Surface of the virion colored by electrostatic potential value from −50 mV (red) to +50 mV (blue) (b). Electrostatic encounter complexes of ZnPcChol8+ with SARS-CoV-2 virion obtained by BD (c,d). The virion is colored in accordance with the envelope component types (c) and electrostatic potential (d). CG beads of neutral lipids, namely POPC, POPE and CHOL, are shown in gray; CG beads of negatively charged lipids are shown in shades of red: POPI molecules—in brown; POPS—in magenta; CDL2—in red. M proteins are shown in pale yellow: E proteins—in white; S proteins—in mustard with black glycans. ZnPcChol8+ molecules are visualized as cyan spheres.
Figure 1. Three-dimensional structure of CG molecular model and electrostatic surface potential of SARS-CoV-2 virus envelope. CG beads colored in accordance with envelope component types (a). Rectangular boxes pick up the areas shown in Figure 2, Figure 3 and Figure 4 (green, yellow and blue, respectively). Surface of the virion colored by electrostatic potential value from −50 mV (red) to +50 mV (blue) (b). Electrostatic encounter complexes of ZnPcChol8+ with SARS-CoV-2 virion obtained by BD (c,d). The virion is colored in accordance with the envelope component types (c) and electrostatic potential (d). CG beads of neutral lipids, namely POPC, POPE and CHOL, are shown in gray; CG beads of negatively charged lipids are shown in shades of red: POPI molecules—in brown; POPS—in magenta; CDL2—in red. M proteins are shown in pale yellow: E proteins—in white; S proteins—in mustard with black glycans. ZnPcChol8+ molecules are visualized as cyan spheres.
Ijms 23 07304 g001
Ijms 23 07304 g002 550
Figure 2. A fragment of the envelope of the SARS-CoV-2 virus, with M proteins and patches of charged lipids, which is an enlarged view of the rectangular area highlighted in green in Figure 1a (a). Molecular surface distribution of electrostatic potential from −100 mV (red) to +100 mV (blue) (b). Positions of the Zn atom in electrostatic encounter complexes of ZnPcChol8+ with SARS-CoV-2 virion obtained by BD (c,d). Virus envelope is colored in accordance with component types (c) and electrostatic potential (d). Neutral lipids, namely POPC, POPE, and CHOL, are shown as gray spheres; POPI molecules—as brown spheres; POPS—as magenta spheres; CDL2—as red spheres. Positions of the Zn atom of ZnPcChol8+ are shown as cyan spheres.
Figure 2. A fragment of the envelope of the SARS-CoV-2 virus, with M proteins and patches of charged lipids, which is an enlarged view of the rectangular area highlighted in green in Figure 1a (a). Molecular surface distribution of electrostatic potential from −100 mV (red) to +100 mV (blue) (b). Positions of the Zn atom in electrostatic encounter complexes of ZnPcChol8+ with SARS-CoV-2 virion obtained by BD (c,d). Virus envelope is colored in accordance with component types (c) and electrostatic potential (d). Neutral lipids, namely POPC, POPE, and CHOL, are shown as gray spheres; POPI molecules—as brown spheres; POPS—as magenta spheres; CDL2—as red spheres. Positions of the Zn atom of ZnPcChol8+ are shown as cyan spheres.
Ijms 23 07304 g002
Ijms 23 07304 g003 550
Figure 3. Detailed view of the surface distribution of electrostatic potential from −100 mV (red) to +100 mV (blue) of the pentameric E protein, embedded into the viral membrane (a), and its complexes with ZnPcChol8+ (b). CG beads representing the Zn atom of ZnPcChol8+ molecules in complexes with the SARS-CoV-2 envelope are shown as cyan spheres. Green line contours the borders of the E protein. Negatively charged surface amino acids of the E protein and negatively charged surrounding lipid molecules are labeled in panel (a). The view corresponds to the yellow frame in Figure 1a.
Figure 3. Detailed view of the surface distribution of electrostatic potential from −100 mV (red) to +100 mV (blue) of the pentameric E protein, embedded into the viral membrane (a), and its complexes with ZnPcChol8+ (b). CG beads representing the Zn atom of ZnPcChol8+ molecules in complexes with the SARS-CoV-2 envelope are shown as cyan spheres. Green line contours the borders of the E protein. Negatively charged surface amino acids of the E protein and negatively charged surrounding lipid molecules are labeled in panel (a). The view corresponds to the yellow frame in Figure 1a.
Ijms 23 07304 g003
Ijms 23 07304 g004 550
Figure 4. Effects of SARS-CoV-2 S protein glycosylation on ZnPcChol8+ binding. The local region of SARS-CoV-2 envelope with the S protein embedded into viral membrane (a), and molecular surface distribution of electrostatic potential from −100 mV (red) to +100 mV (blue) (b). Lipids and M proteins are the same as in Figure 1a,d. The S protein is shown as a mustard surface with black neutrally charged glycans and fuchsia negatively charged glycans. CG beads representing the Zn atom of ZnPcChol8+ molecules in complexes with SARS-CoV-2 envelope are shown as cyan spheres (b). The view corresponds to the blue frame in Figure 1a. The binding ability of non-glycosylated (c) and glycosylated (d) S proteins in relation to ZnPcChol8+.
Figure 4. Effects of SARS-CoV-2 S protein glycosylation on ZnPcChol8+ binding. The local region of SARS-CoV-2 envelope with the S protein embedded into viral membrane (a), and molecular surface distribution of electrostatic potential from −100 mV (red) to +100 mV (blue) (b). Lipids and M proteins are the same as in Figure 1a,d. The S protein is shown as a mustard surface with black neutrally charged glycans and fuchsia negatively charged glycans. CG beads representing the Zn atom of ZnPcChol8+ molecules in complexes with SARS-CoV-2 envelope are shown as cyan spheres (b). The view corresponds to the blue frame in Figure 1a. The binding ability of non-glycosylated (c) and glycosylated (d) S proteins in relation to ZnPcChol8+.
Ijms 23 07304 g004
Table
Table 1. Statistics of contacts of ZnPcChol8+ molecule with different components of the viral envelope.
Table 1. Statistics of contacts of ZnPcChol8+ molecule with different components of the viral envelope.
Components of the Viral EnvelopeFraction, %
proteinsS33.9
M8.4
E0.6
lipidsnegatively chargedPOPI43.9
CDL240.6
POPS11.5
unchargedPOPC61.1
POPE39.4
CHOL0.5
Table
Table 2. Distribution of electrostatic complexes of ZnPcChol8+ with negatively charged lipids by the number of contacts with different lipid types. The cases with occurrence less than 2% are omitted.
Table 2. Distribution of electrostatic complexes of ZnPcChol8+ with negatively charged lipids by the number of contacts with different lipid types. The cases with occurrence less than 2% are omitted.
Number of Lipid Partners of ZnPcChol8+POPIPOPSCDL2Fraction, %
100118
10103
110013
20113
20027
22009
21104
210119
31113
31023
32102
32017
33003
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Fedorov, V.; Kholina, E.; Khruschev, S.; Kovalenko, I.; Rubin, A.; Strakhovskaya, M. Electrostatic Map of the SARS-CoV-2 Virion Specifies Binding Sites of the Antiviral Cationic Photosensitizer. Int. J. Mol. Sci. 2022, 23, 7304. https://doi.org/10.3390/ijms23137304

AMA Style

Fedorov V, Kholina E, Khruschev S, Kovalenko I, Rubin A, Strakhovskaya M. Electrostatic Map of the SARS-CoV-2 Virion Specifies Binding Sites of the Antiviral Cationic Photosensitizer. International Journal of Molecular Sciences. 2022; 23(13):7304. https://doi.org/10.3390/ijms23137304

Chicago/Turabian Style

Fedorov, Vladimir, Ekaterina Kholina, Sergei Khruschev, Ilya Kovalenko, Andrew Rubin, and Marina Strakhovskaya. 2022. "Electrostatic Map of the SARS-CoV-2 Virion Specifies Binding Sites of the Antiviral Cationic Photosensitizer" International Journal of Molecular Sciences 23, no. 13: 7304. https://doi.org/10.3390/ijms23137304

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop