Membrane Activity and Viroporin Assembly for the SARS-CoV-2 E Protein Are Regulated by Cholesterol.

The SARS-CoV-2 E protein is an enigmatic viral structural protein with reported viroporin activity associated with the acute respiratory symptoms of COVID-19, as well as the ability to deform cell membranes for viral budding. Like many viroporins, the E protein is thought to oligomerize with a well-defined stoichiometry. However, attempts to determine the structure of the protein complex have yielded inconclusive results, suggesting several possible oligomers, ranging from dimers to pentamers. Here, we combined patch-clamp, confocal fluorescence microscopy on giant unilamellar vesicles, and atomic force microscopy to show that E protein can exhibit two modes of membrane activity depending on membrane lipid composition. In the absence or the presence of a low content of cholesterol, the protein forms short-living transient pores, which are seen as semi-transmembrane defects in a membrane by atomic force microscopy. Approximately 30 mol% cholesterol is a threshold for the transition to the second mode of conductance, which could be a stable pentameric channel penetrating the entire lipid bilayer. Therefore, the E-protein has at least two different types of activity on membrane permeabilization, which are regulated by the amount of cholesterol in the membrane lipid composition and could be associated with different types of protein oligomers.

Oligomeric State and Drug Binding of the SARS-CoV-2 Envelope Protein Are Sensitive to the Ectodomain.

The envelope (E) protein of SARS-CoV-2 is the smallest of the three structural membrane proteins of the virus. E mediates budding of the progeny virus in the endoplasmic reticulum Golgi intermediate compartment of the cell. It also conducts ions, and this channel activity is associated with the pathogenicity of SARS-CoV-2. The structural basis for these functions is still poorly understood. Biochemical studies of E in detergent micelles found a variety of oligomeric states, but recent 19F solid-state NMR data indicated that the transmembrane domain (ETM, residues 8-38) forms pentamers in lipid bilayers. Hexamethylene amiloride (HMA), an E inhibitor, binds the pentameric ETM at the lipid-exposed helix-helix interface. Here, we investigate the oligomeric structure and drug interaction of an ectodomain-containing E construct, ENTM (residues 1-41). Unexpectedly, 19F spin diffusion NMR data reveal that ENTM adopts an average oligomeric state of dimers instead of pentamers in lipid bilayers. A new amiloride inhibitor, AV-352, shows stronger inhibitory activity than HMA in virus-like particle assays. Distance measurements between 13C-labeled protein and a trifluoromethyl group of AV-352 indicate that the drug binds ENTM with a higher stoichiometry than ETM. We measured protein-drug contacts using a sensitivity-enhanced two-dimensional 13C-19F distance NMR technique. The results indicate that AV-352 binds the C-terminal half of the TM domain, similar to the binding region of HMA. These data provide evidence for the existence of multiple oligomeric states of E in lipid bilayers, which may carry out distinct functions and may be differentially targeted by antiviral drugs.

The mutual lipid-mediated effect of the transmembrane domain of SARS-CoV-2 E-protein and glycyrrhizin nicotinate derivatives on the localization in the lipid bilayer.

Glycyrrhizinic acid (GA) is one of the active substances in licorice root. It exhibits antiviral activity against various enveloped viruses, for example, SARS-CoV-2. GA derivatives are promising biologically active compounds from perspective of developing broad-spectrum antiviral agents. Given that GA nicotinate derivatives (Glycyvir) demonstrate activity against various DNA- and RNA-viruses, a search for a possible mechanism of action of these compounds is required. In the present paper, the interaction of Glycyvir with the transmembrane domain of the SARS-CoV-2 E-protein (ETM) in a model lipid membrane was investigated by NMR spectroscopy and molecular dynamics simulation. The lipid-mediated influence on localization of the SARS-CoV-2 E-protein by Glycyvir was observed. The presence of Glycyvir leads to deeper immersion of the ETM in lipid bilayer. Taking into account that E-protein plays a significant role in virus production and takes part in virion assembly and budding, the data on the effect of potential antiviral agents on ETM localization and structure in the lipid environment may provide a basis for further studies of potential coronavirus E-protein inhibitors.

PDZ2-conjugated-PLGA nanoparticles are tiny heroes in the battle against SARS-CoV-2.

The COVID-19 pandemic caused by SARS-CoV-2 has highlighted the urgent need for innovative antiviral strategies to fight viral infections. Although a substantial part of the overall effort has been directed at the Spike protein to create an effective global vaccination strategy, other proteins have also been examined and identified as possible therapeutic targets. Among them, although initially underestimated, there is the SARS-CoV-2 E-protein, which turned out to be a key factor in viral pathogenesis due to its role in virus budding, assembly and spreading. The C-terminus of E-protein contains a PDZ-binding motif (PBM) that plays a key role in SARS-CoV-2 virulence as it is recognized and bound by the PDZ2 domain of the human tight junction protein ZO-1. The binding between the PDZ2 domain of ZO-1 and the C-terminal portion of SARS-CoV-2 E-protein has been extensively characterized. Our results prompted us to develop a possible adjuvant therapeutic strategy aimed at slowing down or inhibiting virus-mediated pathogenesis. Such innovation consists in the design and synthesis of externally PDZ2-ZO1 functionalized PLGA-based nanoparticles to be used as intracellular decoy. Contrary to conventional strategies, this innovative approach aims to capitalize on the E protein-PDZ2 interaction to prevent virus assembly and replication. In fact, the conjugation of the PDZ2 domain to polymeric nanoparticles increases the affinity toward the E protein effectively creating a "molecular sponge" able to sequester E proteins within the intracellular environment of infected cells. Our in vitro studies on selected cellular models, show that these nanodevices significantly reduce SARS-CoV-2-mediated virulence, emphasizing the importance of exploiting viral-host interactions for therapeutic benefit.

A short sequence in the tail of SARS-CoV-2 envelope protein controls accessibility of its PDZ-binding motif to the cytoplasm.

The carboxy-terminal tail of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) envelope protein (E) contains a PDZ-binding motif (PBM) which is crucial for coronavirus pathogenicity. During SARS-CoV-2 infection, the viral E protein is expressed within the Golgi apparatus membrane of host cells with its PBM facing the cytoplasm. In this work, we study the molecular mechanisms controlling the presentation of the PBM to host PDZ (PSD-95/Dlg/ZO-1) domain-containing proteins. We show that at the level of the Golgi apparatus, the PDZ-binding motif of the E protein is not detected by E C-terminal specific antibodies nor by the PDZ domain-containing protein-binding partner. Four alanine substitutions upstream of the PBM in the central region of the E protein tail is sufficient to generate immunodetection by anti-E antibodies and trigger robust recruitment of the PDZ domain-containing protein into the Golgi organelle. Overall, this work suggests that the presentation of the PBM to the cytoplasm is under conformational regulation mediated by the central region of the E protein tail and that PBM presentation probably does not occur at the surface of Golgi cisternae but likely at post-Golgi stages of the viral cycle.

SARS-CoV-2 Envelope Protein Forms Clustered Pentamers in Lipid Bilayers.

The SARS-CoV-2 envelope (E) protein is a viroporin associated with the acute respiratory symptoms of COVID-19. E forms cation-selective ion channels that assemble in the lipid membrane of the endoplasmic reticulum Golgi intermediate compartment. The channel activity of E is linked to the inflammatory response of the host cell to the virus. Like many viroporins, E is thought to oligomerize with a well-defined stoichiometry. However, attempts to determine the E stoichiometry have led to inconclusive results and suggested mixtures of oligomers whose exact nature might vary with the detergent used. Here, we employ 19F solid-state nuclear magnetic resonance and the centerband-only detection of exchange (CODEX) technique to determine the oligomeric number of E's transmembrane domain (ETM) in lipid bilayers. The CODEX equilibrium value, which corresponds to the inverse of the oligomeric number, indicates that ETM assembles into pentamers in lipid bilayers, without any detectable fraction of low-molecular-weight oligomers. Unexpectedly, at high peptide concentrations and in the presence of the lipid phosphatidylinositol, the CODEX data indicate that more than five 19F spins are within a detectable distance of about 2 nm, suggesting that the ETM pentamers cluster in the lipid bilayer. Monte Carlo simulations that take into account peptide-peptide and peptide-lipid interactions yielded pentamer clusters that reproduced the CODEX data. This supramolecular organization is likely important for E-mediated virus assembly and budding and for the channel function of the protein.

Effect of an Amyloidogenic SARS-COV-2 Protein Fragment on α-Synuclein Monomers and Fibrils.

Aggregates of α-synuclein are thought to be the disease-causing agent in Parkinson's disease. Various case studies have hinted at a correlation between COVID-19 and the onset of Parkinson's disease. For this reason, we use molecular dynamics simulations to study whether amyloidogenic regions in SARS-COV-2 proteins can initiate and modulate aggregation of α-synuclein. As an example, we choose the nine-residue fragment SFYVYSRVK (SK9), located on the C-terminal of the envelope protein of SARS-COV-2. We probe how the presence of SK9 affects the conformational ensemble of α-synuclein monomers and the stability of two resolved fibril polymorphs. We find that the viral protein fragment SK9 may alter α-synuclein amyloid formation by shifting the ensemble toward aggregation-prone and preferentially rod-like fibril seeding conformations. However, SK9 has only a small effect on the stability of pre-existing or newly formed fibrils. A potential mechanism and key residues for potential virus-induced amyloid formation are described.

Structure and dynamics of the SARS-CoV-2 envelope protein monomer.

Coronaviruses, especially severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), present an ongoing threat to human wellbeing. Consequently, elucidation of molecular determinants of their function and interaction with the host is an important task. Whereas some of the coronaviral proteins are extensively characterized, others remain understudied. Here, we use molecular dynamics simulations to analyze the structure and dynamics of the SARS-CoV-2 envelope (E) protein (a viroporin) in the monomeric form. The protein consists of the hydrophobic α-helical transmembrane domain (TMD) and amphiphilic α-helices H2 and H3, connected by flexible linkers. We show that TMD has a preferable orientation in the membrane, while H2 and H3 reside at the membrane surface. Orientation of H2 is strongly influenced by palmitoylation of cysteines Cys40, Cys43, and Cys44. Glycosylation of Asn66 affects the orientation of H3. We also observe that the monomeric E protein both generates and senses the membrane curvature, preferably localizing with the C-terminus at the convex regions of the membrane; the protein in the pentameric form displays these properties as well. Localization to curved regions may be favorable for assembly of the E protein oligomers, whereas induction of curvature may facilitate the budding of the viral particles. The presented results may be helpful for a better understanding of the function of the coronaviral E protein and viroporins in general, and for overcoming the ongoing SARS-CoV-2 pandemic.

Elucidation of SARS-Cov-2 Budding Mechanisms through Molecular Dynamics Simulations of M and E Protein Complexes.

SARS-CoV-2 and other coronaviruses pose major threats to global health, yet computational efforts to understand them have largely overlooked the process of budding, a key part of the coronavirus life cycle. When expressed together, coronavirus M and E proteins are sufficient to facilitate budding into the ER-Golgi intermediate compartment (ERGIC). To help elucidate budding, we ran atomistic molecular dynamics (MD) simulations using the Feig laboratory's refined structural models of the SARS-CoV-2 M protein dimer and E protein pentamer. Our MD simulations consisted of M protein dimers and E protein pentamers in patches of membrane. By examining where these proteins induced membrane curvature in silico, we obtained insights around how the budding process may occur. Multiple M protein dimers acted together to induce global membrane curvature through protein-lipid interactions while E protein pentamers kept the membrane planar. These results could eventually help guide development of antiviral therapeutics that inhibit coronavirus budding.

Immunodominant linear B cell epitopes in the spike and membrane proteins of SARS-CoV-2 identified by immunoinformatics prediction and immunoassay.

SARS-CoV-2 continues to infect an ever-expanding number of people, resulting in an increase in the number of deaths globally. With the emergence of new variants, there is a corresponding decrease in the currently available vaccine efficacy, highlighting the need for greater insights into the viral epitope profile for both vaccine design and assessment. In this study, three immunodominant linear B cell epitopes in the SARS-CoV-2 spike receptor-binding domain (RBD) were identified by immunoinformatics prediction, and confirmed by ELISA with sera from Macaca fascicularis vaccinated with a SARS-CoV-2 RBD subunit vaccine. Further immunoinformatics analyses of these three epitopes gave rise to a method of linear B cell epitope prediction and selection. B cell epitopes in the spike (S), membrane (M), and envelope (E) proteins were subsequently predicted and confirmed using convalescent sera from COVID-19 infected patients. Immunodominant epitopes were identified in three regions of the S2 domain, one region at the S1/S2 cleavage site and one region at the C-terminus of the M protein. Epitope mapping revealed that most of the amino acid changes found in variants of concern are located within B cell epitopes in the NTD, RBD, and S1/S2 cleavage site. This work provides insights into B cell epitopes of SARS-CoV-2 as well as immunoinformatics methods for B cell epitope prediction, which will improve and enhance SARS-CoV-2 vaccine development against emergent variants.

SARS-CoV-2 structural coverage map reveals viral protein assembly, mimicry, and hijacking mechanisms.

We modeled 3D structures of all SARS-CoV-2 proteins, generating 2,060 models that span 69% of the viral proteome and provide details not available elsewhere. We found that ˜6% of the proteome mimicked human proteins, while ˜7% was implicated in hijacking mechanisms that reverse post-translational modifications, block host translation, and disable host defenses; a further ˜29% self-assembled into heteromeric states that provided insight into how the viral replication and translation complex forms. To make these 3D models more accessible, we devised a structural coverage map, a novel visualization method to show what is-and is not-known about the 3D structure of the viral proteome. We integrated the coverage map into an accompanying online resource (https://aquaria.ws/covid) that can be used to find and explore models corresponding to the 79 structural states identified in this work. The resulting Aquaria-COVID resource helps scientists use emerging structural data to understand the mechanisms underlying coronavirus infection and draws attention to the 31% of the viral proteome that remains structurally unknown or dark.

Computational Study on the Function of Palmitoylation on the Envelope Protein in SARS-CoV-2.

SARS-CoV-2 that caused COVID-19 has spread since the end of 2019. Its major effects resulted in over four million deaths around the whole world by August 2021. Therefore, understanding virulence mechanisms is important to prevent future outbreaks and for COVID-19 drug development. The envelope (E) protein is an important structural protein, affecting virus assembly and budding. The E protein pentamer is a viroporin, serving as an ion transferring channel in cells. In this work, we applied molecular dynamic simulations and topological and electrostatic analyses to study the effects of palmitoylation on the E protein pentamer. The results indicate that the cation transferring direction is more from the lumen to the cytosol. The structure of the palmitoylated E protein pentamer is more stable while the loss of palmitoylation caused the pore radius to reduce and even collapse. The electrostatic forces on the two sides of the palmitoylated E protein pentamer are more beneficial to attract cations in the lumen and to release cations into the cytosol. The results indicate the importance of palmitoylation, which can help the drug design for the treatment of COVID-19.

Structural basis for SARS-CoV-2 envelope protein recognition of human cell junction protein PALS1.

The COVID-19 pandemic, caused by the SARS-CoV-2 virus, has created global health and economic emergencies. SARS-CoV-2 viruses promote their own spread and virulence by hijacking human proteins, which occurs through viral protein recognition of human targets. To understand the structural basis for SARS-CoV-2 viral-host protein recognition, here we use cryo-electron microscopy (cryo-EM) to determine a complex structure of the human cell junction protein PALS1 and SARS-CoV-2 viral envelope (E) protein. Our reported structure shows that the E protein C-terminal DLLV motif recognizes a pocket formed exclusively by hydrophobic residues from the PDZ and SH3 domains of PALS1. Our structural analysis provides an explanation for the observation that the viral E protein recruits PALS1 from lung epithelial cell junctions. In addition, our structure provides novel targets for peptide- and small-molecule inhibitors that could block the PALS1-E interactions to reduce E-mediated virulence.

Structural basis of coronavirus E protein interactions with human PALS1 PDZ domain.

SARS-CoV-2 infection leads to coronavirus disease 2019 (COVID-19), which is associated with severe and life-threatening pneumonia and respiratory failure. However, the molecular basis of these symptoms remains unclear. SARS-CoV-1 E protein interferes with control of cell polarity and cell-cell junction integrity in human epithelial cells by binding to the PALS1 PDZ domain, a key component of the Crumbs polarity complex. We show that C-terminal PDZ binding motifs of SARS-CoV-1 and SARS-CoV-2 E proteins bind the PALS1 PDZ domain with 29.6 and 22.8 µM affinity, whereas the related sequence from MERS-CoV did not bind. We then determined crystal structures of PALS1 PDZ domain bound to both SARS-CoV-1 and SARS-CoV-2 E protein PDZ binding motifs. Our findings establish the structural basis for SARS-CoV-1/2 mediated subversion of Crumbs polarity signalling and serve as a platform for the development of small molecule inhibitors to suppress SARS-CoV-1/2 mediated disruption of polarity signalling in epithelial cells.

Quasispecies of SARS-CoV-2 revealed by single nucleotide polymorphisms (SNPs) analysis.

New SARS-CoV-2 mutants have been continuously indentified with enhanced transmission ever since its outbreak in early 2020. As an RNA virus, SARS-CoV-2 has a high mutation rate due to the low fidelity of RNA polymerase. To study the single nucleotide polymorphisms (SNPs) dynamics of SARS-CoV-2, 158 SNPs with high confidence were identified by deep meta-transcriptomic sequencing, and the most common SNP type was C > T. Analyses of intra-host population diversity revealed that intra-host quasispecies' composition varies with time during the early onset of symptoms, which implicates viral evolution during infection. Network analysis of co-occurring SNPs revealed the most abundant non-synonymous SNP 22,638 in the S glycoprotein RBD region and 28,144 in the ORF8 region. Furthermore, SARS-CoV-2 variations differ in an individual's respiratory tissue (nose, throat, BALF, or sputum), suggesting independent compartmentalization of SARS-CoV-2 populations in patients. The positive selection analysis of the SARS-CoV-2 genome uncovered the positive selected amino acid G251V on ORF3a. Alternative allele frequency spectrum (AAFS) of all variants revealed that ORF8 could bear alternate alleles with high frequency. Overall, the results show the quasispecies' profile of SARS-CoV-2 in the respiratory tract in the first two months after the outbreak.

Interactions of SARS-CoV-2 envelope protein with amilorides correlate with antiviral activity.

SARS-CoV-2 is the novel coronavirus that is the causative agent of COVID-19, a sometimes-lethal respiratory infection responsible for a world-wide pandemic. The envelope (E) protein, one of four structural proteins encoded in the viral genome, is a 75-residue integral membrane protein whose transmembrane domain exhibits ion channel activity and whose cytoplasmic domain participates in protein-protein interactions. These activities contribute to several aspects of the viral replication-cycle, including virion assembly, budding, release, and pathogenesis. Here, we describe the structure and dynamics of full-length SARS-CoV-2 E protein in hexadecylphosphocholine micelles by NMR spectroscopy. We also characterized its interactions with four putative ion channel inhibitors. The chemical shift index and dipolar wave plots establish that E protein consists of a long transmembrane helix (residues 8-43) and a short cytoplasmic helix (residues 53-60) connected by a complex linker that exhibits some internal mobility. The conformations of the N-terminal transmembrane domain and the C-terminal cytoplasmic domain are unaffected by truncation from the intact protein. The chemical shift perturbations of E protein spectra induced by the addition of the inhibitors demonstrate that the N-terminal region (residues 6-18) is the principal binding site. The binding affinity of the inhibitors to E protein in micelles correlates with their antiviral potency in Vero E6 cells: HMA ≈ EIPA > DMA >> Amiloride, suggesting that bulky hydrophobic groups in the 5' position of the amiloride pyrazine ring play essential roles in binding to E protein and in antiviral activity. An N15A mutation increased the production of virus-like particles, induced significant chemical shift changes from residues in the inhibitor binding site, and abolished HMA binding, suggesting that Asn15 plays a key role in maintaining the protein conformation near the binding site. These studies provide the foundation for complete structure determination of E protein and for structure-based drug discovery targeting this protein.

Characterization of the SARS-CoV-2 E Protein: Sequence, Structure, Viroporin, and Inhibitors.

The COVID-19 epidemic is one of the most influential epidemics in history. Understanding the impact of coronaviruses (CoVs) on host cells is very important for disease treatment. The SARS-CoV-2 envelope (E) protein is a small structural protein involved in many aspects of the viral life cycle. The E protein promotes the packaging and reproduction of the virus, and deletion of this protein weakens or even abolishes the virulence. This review aims to establish new knowledge by combining recent advances in the study of the SARS-CoV-2 E protein and by comparing it with the SARS-CoV E protein. The E protein amino acid sequence, structure, self-assembly characteristics, viroporin mechanisms and inhibitors are summarized and analyzed herein. Although the mechanisms of the SARS-CoV-2 and SARS-CoV E proteins are similar in many respects, specific studies on the SARS-CoV-2 E protein, for both monomers and oligomers, are still lacking. A comprehensive understanding of this protein should prompt further studies on the design and characterization of effective targeted therapeutic measures.

A study on non-synonymous mutational patterns in structural proteins of SARS-CoV-2.

SARS-CoV-2 is mutating and creating divergent variants across the world. An in-depth investigation of the amino acid substitutions in the genomic signature of SARS-CoV-2 proteins is highly essential for understanding its host adaptation and infection biology. A total of 9587 SARS-CoV-2 structural protein sequences collected from 49 different countries are used to characterize protein-wise variants, substitution patterns (type and location), and major substitution changes. The majority of the substitutions are distinct, mostly in a particular location, and lead to a change in an amino acid's biochemical properties. In terms of mutational changes, envelope (E) and membrane (M) proteins are relatively more stable than nucleocapsid (N) and spike (S) proteins. Several co-occurrence substitutions are observed, particularly in S and N proteins. Substitution specific to active sub-domains reveals that heptapeptide repeat, fusion peptides, transmembrane in S protein, and N-terminal and C-terminal domains in the N protein are remarkably mutated. We also observe a few deleterious mutations in the above domains. The overall study on non-synonymous mutation in structural proteins of SARS-CoV-2 at the start of the pandemic indicates a diversity amongst virus sequences.

Evolution of SARS-CoV-2 Envelope, Membrane, Nucleocapsid, and Spike Structural Proteins from the Beginning of the Pandemic to September 2020: A Global and Regional Approach by Epidemiological Week.

Monitoring acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genetic diversity and emerging mutations in this ongoing pandemic is crucial for understanding its evolution and assuring the performance of diagnostic tests, vaccines, and therapies against coronavirus disease (COVID-19). This study reports on the amino acid (aa) conservation degree and the global and regional temporal evolution by epidemiological week for each residue of the following four structural SARS-CoV-2 proteins: spike, envelope, membrane, and nucleocapsid. All, 105,276 worldwide SARS-CoV-2 complete and partial sequences from 117 countries available in the Global Initiative on Sharing All Influenza Data (GISAID) from 29 December 2019 to 12 September 2020 were downloaded and processed using an in-house bioinformatics tool. Despite the extremely high conservation of SARS-CoV-2 structural proteins (>99%), all presented aa changes, i.e., 142 aa changes in 65 of the 75 envelope aa, 291 aa changes in 165 of the 222 membrane aa, 890 aa changes in 359 of the 419 nucleocapsid aa, and 2671 changes in 1132 of the 1273 spike aa. Mutations evolution differed across geographic regions and epidemiological weeks (epiweeks). The most prevalent aa changes were D614G (81.5%) in the spike protein, followed by the R203K and G204R combination (37%) in the nucleocapsid protein. The presented data provide insight into the genetic variability of SARS-CoV-2 structural proteins during the pandemic and highlights local and worldwide emerging aa changes of interest for further SARS-CoV-2 structural and functional analysis.