key: cord-0776656-g5zhoyo9
authors: Mukherjee, Shruti; Bhattacharyya, Dipita; Bhunia, Anirban
title: Host-membrane interacting interface of the SARS coronavirus envelope protein: Immense functional potential of C-terminal domain
date: 2020-08-11
journal: Biophys Chem
DOI: 10.1016/j.bpc.2020.106452
sha: 1b468fd8198807250b4cc621a532bddab11b406b
doc_id: 776656
cord_uid: g5zhoyo9

The Envelope (E) protein in SARS Coronavirus (CoV) is a small structural protein, incorporated as part of the envelope. A major fraction of the protein has been known to be associated with the host membranes, particularly organelles related to intracellular trafficking, prompting CoV packaging and propagation. Studies have elucidated the central hydrophobic transmembrane domain of the E protein being responsible for much of the viroporin activity in favor of the virus. However, newer insights into the organizational principles at the membranous compartments within the host cells suggest further complexity of the system. The lesser hydrophobic Carboxylic-terminal of the protein harbors interesting amino acid sequences- suggesting at the prevalence of membrane-directed amyloidogenic properties that remains mostly elusive. These highly conserved segments indicate at several potential membrane-associated functional roles that can redefine our comprehensive understanding of the protein. This should prompt further studies in designing and characterizing of effective targeted therapeutic measures.

induce global metabolic changes on infected cells, leading to the rearrangement of the lipid dynamics to facilitate viral multiplication. (27) In some cases, these alterations produce the reorganization of intracellular membranes of the host cell, building an adequate microenvironment for viral replication. All these findings highlight the intimate connections between viruses and host lipid interactions.

Typical of the genus, the initial steps of SARS-CoVs entry include the attachment of the virus particle to a cellular receptor, angiotensin-converting enzyme 2 (ACE2-specific for SARS-CoV) located on the host plasma membrane surface. (28) The viral genome has to take entry into the host cell to reach the replication sites. Different lipids found either on the plasma membrane and/or the endosomal membranes can contribute to these processes by enabling receptor clustering, virus internalization, or membrane fusion.

The host spectrum/tropism of individual CoVs have been known to be primarily determined by the S protein. (29) The functional aspects of the S protein have been known to mediate receptor binding and virus-cell fusion that occur at either the interface of the plasma membrane or within the endosomal vesicles. (30) The mechanism of membrane fusion of SARS-CoV-1 mostly belongs to the type I fusion system, where S protein of SARS-CoV-1 comprises of a number of membrane binding regions in the S2 domain identified as fusogenic peptide or FPs by several researchers. (31) (32) (33) Recently Bhattacharjya et. al. have shown that one of these fusion peptides adopts β-sheet conformation either in the presence of phosphatidylcholine or phosphatidylcholine/cholesterol comprising membrane models with the help of circular dichroism spectroscopy. (34) Moreover, it has also been shown that the peptide undergoes oligomerization with increasing concentration of cholesterol in the membrane. In addition to receptor binding, the action of host cell-specific proteases can likely cleave the S protein, in turn, regulating the CoV infection and host-or tissue-tropism. (35) Once inside the host cellular environment, positive-stranded RNA viruses such as CoV compartmentalize their replication machinery efficiently within the host membranous compartments to evade detection by host pattern recognition receptors. (36) Such compartments are formed for the benefit of the virus by taking advantage of cellular pathways and lipid-modifying enzymes. (27) After transcription within the rough Endoplasmic Reticulum (ER), the viral membrane proteins associate with the Endoplasmic Reticulum-Golgi intermediate compartment (or the ERGIC) before being folded and promoting the budding off of the viral particle by hijacking the cellular secretory mechanism.(37) A major portion of the viral E protein has been known to remain localized at the site of intracellular trafficking, viz. the ERGIC system, participating in CoV assembly and budding. (38, 39) Thus, much of the lipid components of the virus particles are comparable to these host organellar membranes being rich in phosphatidylcholines, phosphatidylethanolamines and a small fraction of phosphatidylinositols. (40) This has provided direct evidence that the host membranes can serve as crucial functional domains that are indispensable for viral propagation. Thus, mechanistic insight into the functional interface of host-membranes is particularly essential to obtain a comprehensive understanding of the target conformers for effective therapeutic strategies.

E protein of CoV is an integral membrane protein about 76-109 residues long and has been known to be most dynamic at the membranous functional interface. (41) The primary structural elements of the protein make it particularly interesting that delve into its possible role in viral propagation. Interestingly, despite the everevolving nature of the virion particles, the primary sequence of the protein remains mostly conserved. The SARS-CoV-2 E protein bears about 97% sequence similarity with that of the SARS-CoV. This alternatively directs into understanding the uniqueness of the sequence that is conserved to mediate some crucial functional roles of the protein (Figure1).

The primary amino acid sequence reveals a very short N-terminal hydrophilic sequence of about ten amino acids. (42, 43) This is followed by a stretch of hydrophobic sequence that serves as the transmembrane domain (TMD). This twenty-eight residue long stretch has been known for the membrane-induced topologies of the protein that have been extensively defined and studied. (44) The absence of a canonical cleaved signal sequence suggests equal likelihood for the protein to be a type II (with its C-terminal targeted to the ER lumen) or type III (with its N-terminal targeted to the ER lumen) membrane protein. Therefore, the topology of the CoV E protein is highly debated in the literature and remains largely inconclusive. Table 1 summarizes the contradictory topology predicted for some of the widely studied CoV strains. The predicted topologies of the SARS CoV-2 E protein using prediction servers have been represented in Table 2 . However, further complicating matters, the topology models devised using prediction programs are mostly inconsistent with the previous experimental observations. Two membrane-associated topologies have been demonstrated for E protein (Figure2)-hairpin or transmembrane. (24, 45) In a report by Arbely et al. (2004) , the protein was shown to adopt a highly unusual topology, comprising of a short transmembrane helical hairpin conformation that inverses about a previously unidentified pseudo-center of symmetry. (46) This hairpin structure could deform lipid bilayers and cause fragmentation. (46, 47) It has been suggested that post-translational modifications and several intermediate stages of maturation in the ER and Golgi environment can modulate the overall conformation of the molecule in response to the physiochemical properties of the membrane. (48) (49) (50) The specific lipid composition of the particular lipid bilayers of the ER-Golgi protein synthesis-packing machinery of the cell defines divergent physicochemical properties at the membranous interface. The lipid packing density has been known to underlie significant features of the resulting membrane-inducing characteristic local curvatures, surface charge, phase behavior, elasticity, hydrophobicity, degree of hydration, etc.(51) Membrane phase behavior in a cell has been known to permit the transient concentration of specific proteins and lipids into dynamic nanoscopic domains. (52) The level of palmitoylation may moderate the relative proportion of the membrane-associated conformers of the protein that again, in turn, can modulate the local membrane curvatures, facilitating further favorable association of the protein conformers. (53) Substitutions of the hydrophobic amino acid residues in the TMD of the E protein with charged amino acids have been shown to alter the migrating properties of the E protein in SDS-PAGE. (44) This, therefore, suggested changes in the overall conformation and membrane association of the mutant forms in comparison to the wild type E protein.

The TMD is followed by a lesser hydrophobic C-terminus with a fraction incorporated as part of the envelope in the virion particle ( Figure 1 ). The C-terminus, in particular, has been known to be associated with the redirection to the host ER and Golgi complex membranes that underlies some of the major functional domains of the protein. (54, 55) Comprising of more than 50% of the entire protein sequence, this segment has been known to influence the topology in the host cell. There are several conflicting theories that either suggests a luminal localization of the C-terminal when associated with the ERGIC system or directed towards the cytoplasm. (38, 45) The several experimental and in-silico prediction studies have been mostly inconclusive with respect to its localization and hence functional domain in the different strains of CoV. Nevertheless, the amino acid sequence of the C-terminal harbors some of the answers to the many known functional attribute of the protein sequence. Interestingly, it comprises of the conserved "FYxY" peptide sequence (where x= any amino acid) that can be correlated to its high amyloidogenic propensity. (15, 25) Also, another highly conserved Cysteine containing motif-"CxxC" can enable disulfide isomerization to prompt membrane-directed conformational changes. (56) Apart from these highly conserved sequences throughout the genus, there are distinct potent glycosylation sites along the stretch that can serve as chaperone interacting motifs to help in the protein folding and/or aid in J o u r n a l P r e -p r o o f Journal Pre-proof trafficking along with the cellular machinery.(57) Glycosylation of particular asparagine residues (Asn 45, Asn 48, Asn 64, and Asn 68) in the SARS-CoV has been shown to be crucial in maintaining the proteinoligomerization events associated with the host membranes. (58) Conversely, maintaining the glycosylation events have been suggested to be critical to maintaining the broad functional roles of the protein. (59, 60) Further terminally, there is a distinct phosphorylation domain comprising of the "S 67 SR 69 " motif. Additionally, a PDZbinding domain (PBM) in the C-terminal plays an essential role in host cell modifications necessary for viral infection and pathogenesis (Figure 1 ). (61, 62) Mutation studies with the PBM have established its importance for viral infection, making it a significant factor of research for therapeutic and vaccine designing. (61, 62) 2. The E protein as a functional viroporin.

The primary structural features of the E protein have striking similarities with the viroporin class of proteins. (63) (64) (65) Viroporins ideally comprise of about 60-120 amino acids, with a distinct hydrophobic transmembrane domain. They are known for their interactions with membrane surfaces, often resulting in the expansion of the lipid bilayer. (65) Viroporins participate in several viral functions, including the prompting of viral particle release from host cells. (66) Although not essential for viral replication, viroporins modulate the host cellular machinery to make it favorable for viral propagation. They affect host cellular functions, including membrane vesicularization, glycoprotein tracking, and membrane permeability. The transmembrane domain could form hydrophilic pores in the membranes of virus-infected cells by oligomerization. (63, 67) These hydrophilic channels would allow low molecular weight hydrophilic molecules to cross the membrane barrier, leading to the disruption of membrane potential, collapse of ionic gradients, and release of essential compounds from the cell Several studies with SARS CoV E protein have demonstrated the structural and functional similarities with viroporins. It has been demonstrated that the SARS-CoV E protein could enhance the membrane permeability of bacterial cells to o-nitrophenyl-β-d-galactopyranoside and hygromycin B, suggesting that the protein may function as a viroporin. (68) Similar observations were reported on mouse hepatitis virus E protein. (64) The E protein, as a viroporin, is present in low copy numbers in the virus particle, but has been implicated in membrane scission being mostly associated with the ERGIC and cis Golgi-consistent with a predominant role as a mediator of virus assembly and release at this location. (69) The lipid content at these locations may also enhance virus budding. (70) compartments are formed by manipulating the protein-synthesis, packaging, and distribution system of the cell involving the ER-ERGIC-Golgi complex. This helps the viral replication machinery to evade detection by host pattern recognition receptors ( Figure 2 ). Both viral and hijacked host proteins are used in this process, taking advantage of cellular pathways and lipid-modifying enzymes to the benefit of the virus, sequestering newly formed RNAs away from host immune sensors. (75, 76) This eventually allows the virus to manipulate the normal secretory pathway of the host cellular machinery to transport and deliver the virus protein cargo at the site of final packaging, and eventual budding. (69) The virulence of the strain is, in fact, determined by the efficiency in their mechanism to concentrate the RNA synthesis machinery and subsequent packaging within these specialized compartments that are, in turn, modulated by the associated structural proteins.(77) Among the structural proteins, the E protein has been particularly intriguing either alone or in collaboration with the M protein. (71) Interestingly, studies have shown that overexpression of the E protein is sufficient to induce the packaging and assembly of the CoV particles upon efficient manipulation of the host cellular machinery. (78) Rottier et al. suggested that the E protein has no genuine structural function in the virion envelope itself where it occupies frequent, regular positions within the lattice built by M protein. (79) However, its primary role lies in the functional interface of the host packaging and assembling membranes ( Figure 3 ). The frequent but strategic positions within the lattice hint at a possible morphogenetic function to generate the required membrane curvature for the final viral particle formation ( Figure  3 , inset).

Studies have endowed the C-terminal to be associated with the redirection of the E protein to the ER, and Golgi complex for assembly, budding, and intracellular trafficking of infectious virions.(55) Interestingly, sequence analyses of the C-terminal residues suggest the prevalence of some degenerate ER export signal sequences. These are closely comparable to the most common FxxxFxxxF(80) and the [R/K](x)[R/K](81) dibasic motifs that have been known to be responsible for the export of glycosyltransferases. The last 6-residue "R 103 DKLYS 108 " stretch from the IBV E protein terminal has been defined to serve as the ER retention signal. (72) The site-directed mutagenesis of this composite lysine residue (K105) to glutamine resulted in the accumulation of E in the Golgi apparatus. [56] Maturing through the Golgi complex, the E protein is expected to impart some very crucial functions with respect to vesicularization events, by inducing local curvatures at the site of membrane-directed accumulation ( Figure 3 ). Two principle mechanisms have been described for moving and delivering cargo proteins through the secretory pathway of the Golgi complex, the cisternal maturation, and the formation of mega-vesicles. (82, 83) Despite the extensive studies performed, experimental evidence for both the systems remains inadequate and far from being elusive. During CoV infection, virions have been observed in large vesicle depots resembling megavesicles derived from Golgi/ERGIC membranes, indicating that remodeling of the Golgi complex may be crucial for virion trafficking. (84) Understanding the structural transitions that the E protein undergoes, from after its synthesis to its association with the ERGIC membrane complex, would greatly facilitate the characterization of potential targets. We expect the protein to undergo stepwise dynamic structural transitions with specific roles bound to the membranous component. This gradual transition is dependent equally on both the respective membrane composition along with the conformational changes in the protein as it is translated, modified, and matures along with the ER-Golgi network system. Previous studies performed by Knoops as membrane fusion and fission to prompt viral propagation by hijacking the entire cellular machinery. However, progress in our in-depth understanding of this whole machinery has been stymied by the inability to track the functional interface of the protein at the atomic-level resolution in real-time. Identifying these intermediate functional conformers could broaden our chances of finding the best target motif for therapeutic design.

The CoV E protein has been known for its possible pro-apoptotic viroporin activity majorly by mobilizing calcium ions, prompting an overall change in the ionic environment of the homeostatic cell. (86, 87) The virus particles, generally, hijack the cellular homeostatic machinery of the host cell by compromising the membrane integrity that facilitates its regulatory mechanisms within the environment of the host cells. It has been suggested that the E protein can directly disrupt membranes through the formation of ion channels upon membrane-directed oligomerization events mediated upon the insertion of the TMD. In this context, the TMD has received much interest owing to their membrane-spanning attribute. Molecular dynamic simulation studies with synthetic peptides have revealed that the TMD could form ion channels by homo-oligomerization into stable dimers, trimers, and pentamers. 90) Additionally, in silico prediction studies based on high-resolution NMR spectroscopic studies provided useful insight into understanding the probable conformers that define the TMD induced ion channels. Hence, these studies directed the defining of therapeutics targeting the resultant ion channel pores. (91) The hollow structure of the oligomers allows the unregulated efflux of ions from the cellular compartments across the hydrophilic interior of the pore, modulating the membrane potential, altering the ionic environment of the lumen. Alterations in ion concentration would promote translation of viral mRNAs, given the translation of mRNAs from many cytolytic animal viruses is reasonably resistant to high sodium concentrations. (63, 67) Since the secretory pathway of the host cellular machinery is highly sensitive to perturbations in ionic strength, viral propagation can trigger a detrimental apoptotic cascade.

Unregulated ion channel formation upon E protein oligomerization in the host membrane is believed to be one of the primary processes underlying the viral pathophysiology.(92) As discussed previously, much of the studies have been focusing on the TMD for the high molecular weight conformers (multimers) of the protein upon host membrane interactions. (93, 94) However, interestingly, the immediate adjacent segment (towards the C-terminal) hint, at least in part, at a possible role in the membrane directed oligomerization propensity of the molecule. Upon inspection into the C-terminal segment of the protein, we can see that it is marked by the presence of unique motifs that can prompt much of the membrane-directed conformational changes of the protein molecule. The primary sequence of the C-terminal region reveals possible functional characteristics of the protein that are lesser known in literature. Interestingly, the sequence is distinctive of amyloidogenic proteins and peptides that have been studied in the literature for their membrane-directed indispensable functions. Highresolution structural models of amyloid oligomers have been obtained in the context of membrane interactions in neurodegenerative diseases.(95) Membrane-induced oligomerization prompts membrane compromise, often leading to vesicularization into smaller lipid aggregates. (96) (97) (98) Studies with amyloid proteins have hinted at the possibilities of a change in the local curvature being induced by the membrane-associated high molecular weight conformers of the proteins (Figure 3 , inset).(98) Several reports have suggested specific lipid domains to serve as templates for protein anchorage. This is followed by rearrangement of the protein conformers resulting in aggregates that form "wedges" within the bilayer to facilitate curvature change. (99, 100) Such protein aggregates are driven mostly by the hydrophobic interactions that relocate within the acyl chain region of the lipid bilayer. These interactions result in a change in the lipid packing density and hence the overall membrane integrity. (101) The structural similarity of the C-terminal segment of the E protein to amyloid proteomics indicates at analogous purposes that can be complementary to or function in parallel to its viroporin activities. The oligomerization properties of E protein have been defined previously.(25) E protein expressed in bacterial and mammalian cell systems under reducing conditions was shown to exist as monomers. However, upon change to non-reducing conditions, they formed homodimers and homotrimers. (68) The existing literature has focused on the phenomenon of TMD multimerization associated with the ion channel formation that lacks mechanistic details on the initiation. However, site-directed mutagenesis studies revealed that the two-cysteine residues, immediately adjacent to the TMD in the C-terminal, were essential for the oligomerization, leading to the induction of membrane permeability. The TMD is immediately followed by a "CxxC" redox motif, which is highly conserved among the coronaviridae family of viruses. This motif is very crucial for several reasons. The oxidized CxxC motif can maintain the structural topology of the transmembrane region as well as that of its contiguous cytoplasmic domain, inclusive of glycosylation sites, involved in signaling and protein-protein interaction.(102) A probable activation of a thiol in this motif can trigger the participation of the other Cys residues in the formation of the inter-subunit disulfide bond. The activated thiol will then attack the disulfide and cause its isomerization into a disulfide isomer within the motif. This may lead to the refolding of the transmembrane region and may activate its fusogenic potential.(103, 104) Interestingly, a similar CxxC motif is present as a part of the C-terminal of S protein sequence. Bioinformatic studies have suggested at the possibilities of forming inter-protein disulphide bridges that might indicate at a crossover or cooperation between these two structural proteins in the viral membrane interactions. (41) The formation of a disulfide bond may also play a crucial role in the oligomerization of the E protein, forming stable dimers, trimers, and pentamers depending on its functional requirement.(105) Thus even though the TMD spans the lipid bilayer, the CxxC motif could serve as an essential key to defining the membrane-associated oligomerization events-providing newer targets for preemptive therapeutic intervention.

Apart from the consequential involvement of this segment in the viroporin function, the amyloidogenic segment of the E proteins can function independently at the membranous interface. Further downstream, the Cterminal harbors the "FYxY" composite peptide sequence ( Figure 4 ) that has been known for a high propensity of amyloid formation. (106, 107) Interestingly, a self-assembling short peptide segment with the composite sequence has been studied to have a membrane-mediated function. Bhunia and co-workers studied a 9-residue long peptide stretch (TK9, T 55 VYVYSRVK 63 ) in SARS CoV E protein (Figure 4 ) for its physical changes in the presence of both zwitterionic and negatively charged model membrane micelles.(15) The TK9 peptide segment was shown to adopt amyloid fibrillary structures in solution ( Figure 4 ). (25) Further studies with TK9 provided evidence for its aggregation propensity that closely resembles amyloid proteins. Intriguingly, the composite 5-residue peptide motif-T 55 VYV 58 in TK9 was found to serve as the critical sequence motif that may be critical for the beta-sheet formation. (25) However, in vitro measures demonstrated the differential orientation of the TK9 peptide with respect to the membranous environment: correlating well with the structure and function at the membrane-protein interactions interface. This atomic-resolution NMR based structural elucidation in the presence of membrane mimics have improved understanding of the molecular mechanism of CoV infection under physiological conditions. Studies demonstrated the helical propensity of this particular segment from the C-terminal region (Figure 4 ) of the protein to impart host specificity to the composite virus. The peptide's affinity was further J o u r n a l P r e -p r o o f Journal Pre-proof manifested by its pronounced membrane-integrity disruption ability towards the mammalian compared to the bacterial membrane mimic-implicated in the viral pathogenesis. TK9 penetrated deeper into the acyl region of the neutral micelles mimicking the mammalian membrane and its intracellular organelles as opposed to the superficial conformation when in association with the bacterial model membrane mimics (Figure 4 ). This provided evidence for the host membrane-associated conformational change in the protein segments that is specific to a particular lipid composition (species specificity). This study alternatively offered useful insight into the fact that the composite lipid molecules play particularly important roles in the interaction dynamics and conformational uptake. These segment-based studies highlighted the exclusive function of the amyloidogenic segment of the E protein that acts independently of the TMD-associated oligomerization events. Similarly, another short 9-residue peptide, NK9 (N 45 IVNVSLVK 53 ) further upstream, was identified and characterized, being marked by sheet breaker-hydrophobic residues, like Ile, Val, and Leu. (25) The experimental intervention into these segments can provide useful insight into the topology of the protein, its membrane-bound orientation, and eventual functional conformation uptake.

These preliminary studies have provided useful links to the immense possibilities of the functional membrane-directed conformers of the E protein, which remains largely unexplored. Nevertheless, epitope mapping for these interfaces can serve as critical targets for effective therapeutic intervention. Preliminary sequence homology analyses performed to compare the E protein primary sequence of SARS CoV to that of the novel CoV-2, revealed a minimal difference with the "FYxY" motif remaining conserved (Figure 1 ). The first Threonine and the second Valine from the TK9 sequence have only been substituted by a Serine and Phenylalanine residues, respectively ( Figure 4 ). This might indicate at an unperturbed oligomerization potency of the segment. This conserved amyloidogenic sequence indicates a very crucial functional role of the segment that contributes significantly to the virulence. The fact that hydrophobic interactions can provide stability in localization within the acyl-chain layer of the lipid bilayers leaves us with some obvious open questions. Does the C-terminal induce the initial docking of the protein to the host membranes? Do these hydrophobic interactions prompt the downstream oligomerization events to the formation of the membrane-associated multimers? Alternatively, does the lipid-composition specificity of the C-terminal determine the conformational uptake and the eventual functional topology?

These questions are vital for a comprehensive understanding of viral propagation. Interestingly, hydrophobic protection of these segments within the acyl-chain region can serve as open templates for downstream aggregation. This can result in membrane-associated aggregates that, in analogy to the amyloidogenic proteins, can modulate the overall membrane integrity. Thus, it is essential to study the interface of interaction between the E protein and the host membrane, either as whole or directed to particular functional peptide segments.

One of the major functional aspects of CoV E protein has been associated with virus assembly and its subsequent release for propagation in the host physiology. Indirect immunofluorescence microscopy showed that E protein is localized to the Golgi complex in cells transiently expressing IBV E.(108) When co-expressed with IBV M, both from cDNA and in IBV infection, the two proteins are colocalized in Golgi membranes, near the CoV budding site.(109) Thus, even though IBV E is present at low levels in virions, it is apparently expressed at high levels in infected cells near the site of virus assembly. The subcellular localization of the SARS CoV E protein was analyzed using the sera from immunized mice by ELISA and immunofluorescence using cells J o u r n a l P r e -p r o o f infected with recombinant (rSARS-CoV) and viruses lacking the E gene (rSARS-CoV-∆E) as a negative control. (38) Several parallel studies have endowed the C-terminal to be associated with the redirection of the E protein to the ER, and Golgi complex for the intracellular trafficking of infectious virions. (55) Experimental studies have shown that the C-terminal domain of the E protein, in fact, plays a crucial role in virus budding. (110, 111) The deletion of the domain resulted in the free distribution of the mutant protein and a dysfunctional viral assembly. Alternatively, mutations introduced into the cytoplasmic tail of MHV E protein by targeted RNA recombination resulted in elongated virions, further indicating at the critical involvement of the domain in the budding mechanism. (72) The difference in the lipid composition along the secretory pathway might, in fact, have a direct role on the E protein topology. Moving from the ER to the Golgi and finally the plasma membrane, there is a definite gradient in the concentration of specific lipid components that induce a unique physical property of the lipidbilayers in terms of thickness and rigidity. (40) The ER presents a lower fraction of Cholesterol or Sphingolipids that allows maintenance of a relatively fluid lipid compartments for the remodeling associated with the protein association, sorting and accumulation. (112) This fluidity can have a direct role in the gradual compartmentalization of the ER-associated membranes prompting viral particle formations.

The host membrane-mediated oligomerization of the E protein may be responsible for much of the membrane compromise in favor of the viral propagation. The membrane-associated aggregates can compromise the overall lipid bilayer integrity. (113, 114) The hydrophobic interactions involving the lipid acyl regions can result in increased fluidity of the membrane surface, prompting phase separation of the lipid domains. Biomolecular phase separation is suggested to drive the organization of cellular organelles involved in cellular packaging and compartmentalization.(115) Differential distribution of the membrane proteins and their site of interaction with specific lipid rafts, help determine the local curvature and hence characterize propensity to bud off from the surface. (112, 116) The membranous compartments of ER and the Golgi complex presents with a fair fraction of zwitterionic and neutral lipids-i.e., phosphatidylcholine and phosphatidylethanolamines. (112) In vitro studies by Khattari et al have demonstrated that the membrane-directed topology of the SARS CoV E protein imposes a direct constraint on the lipid-bilayer thickness and the acyl-chain ordering. (117) Their studies suggested that the increasing protein concentrations in the organellar membrane along the secretory pathway (greater protein/lipid concentrations), the TMD of the SARS CoV E protein induces bilayer lipid arrangements. Specific inter-molecular interactions with the lipid molecules and intra-chain hydrogen bonds and salt-bridge interactions allow the uptake of a topology that traverses across the membrane bilayer. However, the exact transmembrane topology of the protein was highly debated and remains elusive.(23) Nevertheless, these intra-and intermolecular interactions can be considered to have the ability to rearrange and induce direct morphological changes in the host membranes. The final viral packaging requires assembling of the replicated RNA strands and the other viral proteins before the final membrane budding.

The putative transmembrane domains of the E protein can also serve a 'catalytic' function in the membrane packaging. (118) Cross-talk between distant protein molecules can, in fact, serve as a "zipper" to close the 'neck' of the viral particle as it pinches off from the membrane in the terminal phase of budding. (8, 119) Despite the immense possibilities, much of the progress in our understanding of the actual mechanism of action is still hypothetical and far from being elusive. It is, therefore, essential to study the interface of interaction between the E protein and the replication complex. Particularly, the membrane-mediated aggregation, that compromises the local curvature, inducing the final scission.

CoV gained particular notoriety when the SARS CoV outbreak shook the world around 2002-2003. The aftermath leads to the identification of many newer family members. Soon after that, studies were in full swing all around the globe to determine the mechanism of viral propagation. Several parallel studies provided useful insight into the understanding of the structural and functional uniqueness of the viral proteins and RNA in their search for finding the ideal target for therapeutic intervention. Almost immediately after SARS, the endemic spread of the MERS-CoV in 2012 reignited the urge to gain in-depth insights into the system all over the world for designing effective therapeutic measures. The recent outbreak of the novel COVID-19 soon turned into a pandemic, threatening health standards globally. Lu et al. first published the genomic characterization of COVID-19 and also reported it to be significantly divergent from SARS-CoV, claiming it to be a new human-infecting βcoronavirus. (120) A comparative analysis of small membrane/ envelope protein of different coronaviruses could provide an interpretation of the molecular events that could advance these viruses from developing acute infections to one producing a pandemic. An analysis of the sequence conservation can also help in identifying the regions required in the typical functions of membrane interaction, oligomerization, localization in the host cell, and viral infection. Such information is of immediate significance and would contribute to vaccine designing and facilitate the evaluation of vaccine candidate immunogenicity. Additionally, this would help in reflecting the potential effects of mutational events as the virus is transmitted through the human population. Recent studies have revealed that the homology with a CoV strain isolated from pangolin was ~99%. This indicated that while SARS-CoV-2 have been known to evolve from the bat CoV, pangolins might have helped as an intermediate host. (121) Nevertheless, further research is required to justify these theories.

We performed a BLAST (https://blast.ncbi.nlm.nih.gov/Blast.cgi) search for the E protein and selected the sequences from the different genus of the coronaviridae family. Homology analyses and sequence alignment were conducted using the MEGA software. The Muscle program of the MEGA software was used to perform multiple sequence alignment, and a phylogenetic tree was created by using a maximum likelihood approach ( Figure  5 ). (122) Predictive tools such as TMHMM (123), Memsat(124), Phobius (125) , and the hydrophobic moment plot method(126) of Eisenberg and co-workers allowed the calculation of the transmembrane regions. Much has been defined to correlate the structural conformation of the TMD with its ion channel formation. However, sequence analyses of the C-terminal have revealed crucial information that has been out of the limelight for several years. Despite the fact that the C-terminal has shown the potential to underlie much of the protein's dynamic functionality associated with membrane budding and scission, it has received much lesser attention.

Several secondary structure prediction programs predicted two β-strands within the SARS-CoV E protein tail sequence, including the N 45 IVNVSLVK 53 and the T 55 VYVYSRVK 63 . The predicted β-strands fit the criteria for forming a β-hairpin, which is a simple structural motif with two β-strands linked by a short loop of two to five amino acid residues. Interestingly, a highly conserved proline residue, Pro 54, resides between the predicted βstrands-responsible for the β-coil-β motif. (54) The antiparallel β-strands form hydrogen bonds to stabilize the hairpin structure. Previous reports with the synthetic peptide fragment of E (i.e., I 46 -S 60 ), encompassing the predicted β-hairpin segment, was found to produce ∼100% β-structure and was completely resistant to Hydrogen-deuterium exchange in Fourier transform experiments. Alternatively, the study showed that titration with the drug molecule, HMA (5-(N, N-hexamethylene) amiloride) resulted in large perturbations of the Val 49 and Leu 65 residues, which are far apart in the sequence. The structural model proposed by Jaume Torres's group, therefore, suggested that these two residues were coming spatially close in the membrane-associated pentamer and might belong to different monomers. Comparison with mutation-based data hypothesized this β structure to be in dynamic equilibrium with an α-helical intermediate form. In fact, a delicate balance between J o u r n a l P r e -p r o o f these two forms may alter the membrane-dynamics and subsequent processes in the infected cell (e.g., membrane scission, binding to protein partners, or E protein localization). Disturbing this predicted β-strand region along with the substitution of the conserved proline residues had resulted in the disruption of the cytoplasmic Golgi complex signal, consequently changing the subcellular localization of E protein. (55) Taken together, the results support the hypothesis that the three-dimensional structure of E protein and not the primary sequence dictates much of its function. The sum of amino acid substitutions per site from between sequences is presented in Figure 5 . Analyses were performed using the Poisson correction model. The estimated value of the shape parameter for the discrete Gamma Distribution is ~15.1, suggesting that the occurrence of the conserved amino acid residues among the coronavirus family E protein is not entirely random. However, upon close inspection with other E protein sequences of different species of the same genus, the value decreased to ~6.8. Substitution patterns and rates were estimated under the Jones-Taylor-Thornton (1992) model (+G). (127, 128) The model evolutionary rate differences among sites were obtained using a discrete Gamma distribution (10 categories, [+G]). All the ambiguous positions were eliminated for each of the sequence pairs (pairwise deletion option). Evolutionary analyses were conducted in MEGA X. (129) Collectively, the data reveals that sites observed to be constant among sequences might have crucial biological significance correlated to their function. Besides, deleterious mutations are much more likely to be found in population-level data. This indicates that sites that would predictably be considered as 'invariable' at the phylogenetic level might have transient polymorphisms at the population level. In the analyses of the E protein sequences from the cornaviridae family, altering the number of gamma categories did not have any noticeable effect on the estimate of the substitution rate.

CoV E protein has received lesser attention when compared to the other structural proteins (like S and M) majorly owing to a very low copy number in the viral envelope. (42, 130, 131) It is the smallest and yet has been the most puzzling of the major structural proteins. Earlier deletion based studies had proven that the viral life cycle endures the absence of E protein, implying that other viral genes could counterbalance for its loss. (132) However, recent evidence collected with recombinant CoVs missing E protein, exhibit considerably lowered viral titers, crippled viral development, or produce progeny incapable of further propagation. (10, 11) Several parallel studies have highlighted the crucial role played by the E and M proteins in the propagation of the viral genome into the host cell and virulence. Remarkably, studies have reported that viruses produced by the deletion of the SARS-CoV E gene can be attenuated in at least three different animal models that conferred protection against SARS-CoV pathogenesis. (133, 134) But, the lack of complete information and a limited number of findings have prevented an understanding of the exact mechanism underlying the definite functional role of E protein in viral infection. Nevertheless, these have provided the basis of understanding that the CoV E might be involved in several aspects of the viral replication cycle, including host cell responses such as apoptosis, inflammation, and even autophagy in association with other nonstructural viral proteins.(135) Soon after the identification, studies have been in full swing all around the globe to determine the mechanism of the viroporin action. Several parallel studies from all over provided useful insight into understanding the structural and functional uniqueness of the viral protein in their search for finding the absolute target for therapeutic intervention.

Extensive studies performed in Prof. Torres's laboratory provide useful insight into the structural conformation of the ion channel pores created by E protein that seems to be a very attractive target for therapeutic intervention. (91) Their studies have prompted the designing and characterization of specific inhibitors for these channels that would perturb the primary functional attribute of the protein. Hexamethylene amiloride has been studied to have blocked this E protein-associated ion channel activity in the mammalian cells expressing SARS-J o u r n a l P r e -p r o o f CoV E protein. (136) However, viroporin proteins, like the SARS-CoV E protein, can exhibit a multitude of diverse functions unrelated to their ion-channel properties that are yet to be defined. The dynamic functional interfaces of the protein have been underestimated for its functionality in the propagation. This study can be correlated with earlier docking screening by Gupta et al., where more than 4000 phytochemicals were evaluated on the same protein, and three of those, i.e., belachinal, macaflavanone E, and vibsanol B were found to be particularly promising. (130) Based on the experimental evidences, it is worth mentioning that the two hits, i.e., glycyrrhizic acid and cepharanthine, are supporting their significant activity against the SARS viruses. (130) More recently, Anatoly Chernyshev has established the potential of the SARS-CoV-2 E protein as a pharmaceutical target. (94, 137) These studies have defined 36 approved drugs, which theoretically could block the E protein-induced ion channel and eventually hinder the virus's life cycle. (138) However, these studies have been based on in silico approaches, and hence much of the drug candidates have to be subjected to further experimental investigations.

Considering the lengthy procedure of new drug development, the existing strategy of drug repurposing has turned into one of the preferred solutions for the immediate treatment of SARS-CoV-2 affected individuals. (139) Long-term drug development objectives for the pharmaceutical industry incorporate the identification of inhibitors targeted at the replication or infection routes associated with SARS-CoV-2 or other allied coronavirus infections. (140) 

Knowledge of the coronavirus structural proteins and their mechanisms of viral action can be crucial in defining alternative therapeutic targets. Despite the extensive research efforts worldwide, a definite cure to target CoV propagation is still elusive. Effective therapeutic intervention into viral systems without compromising the host cell integrity requires an in-depth understanding of the specific functional interface. Especially, with the ever-evolving viral components, as seen in the COVID-19 specimen, immense multidimensional research is the need for the hour. Among the structural proteins, the E protein have not been a convincing choice for therapeutic targeting, despite the evidence for a crucial viroporin activity. Much of the knowledge has remained restricted to only its ion channel formation ability. Although this might be serving as a significant functionality, it is definitely not the comprehensive representation of the immense potential that the protein is possibly capable of undertaking. We believe, in-depth understanding of its functional interface with the host-membranes can undermine newer possibilities and functional attributes that have not received the deserved attention so far. Attempts to determine the crucial protein domains associated with the host membranes are innovative in terms of targeted therapeutic designing. Several studies have been gradually exemplifying the protein's essential role in viral packaging and propagation. The conserved sequences in the C-terminal region directly prompt the understanding of evolutionarily crucial functional roles of this segment that is yet to be conclusively defined. The knowledge is incomplete, essential insight into the protein's functional interface should provide newer insight into understanding and structure-function correlation. This should enable the effective designing of targeted therapeutics. The entire portrayals of biologics with characters suitable for inhibiting several key CoV proteins could serve as a basis for drug development. J o u r n a l P r e -p r o o f infected cell Table 2 : List of predicted topologies of SARS CoV-2 Envelope protein using in silico prediction servers. Alternatively, the segment has been studied for the differential membrane-directed functioning. The peptide undergoes a helical conformation when in association with membranous environments. Comparison between bacterial and mammalian model membrane mimicking systems showed the different orientation of the peptide. This suggested the specific functional role of this peptide segment and its membrane-directed structural change. Preliminary studies with the TY5 peptide segment (underlined) had also shown the significance of the "FYxY" sequence, characteristic of amyloid proteins, and peptides. (15, 25) The figure was prepared with https://biorender.com.

J o u r n a l P r e -p r o o f Highlights:

 The SARS CoV Envelope protein is a small structural protein of the virus that has been suggested to have significant viroporin like activity.

 Majority of its function is mediated at the interface of host-membrane interactions J o u r n a l P r e -p r o o f  Focus at the membrane-directed features of the protein provide useful insight into gaining mechanistic insight into its viroporin functions.

 Studies have elaborated the central hydrophobic transmembrane domain of E protein, known to affect ion-channel formation.

 The C-terminal region of the protein show further potential host-membrane directed functional roles.

 The highly conserved amyloidogenic amino acid stretches of the C-terminal suggest at significant contribution to CoV propagation. 

SARS-CoV-2: an Emerging Coronavirus that Causes a Global Threat

World Health Organization declares global emergency: A review of the 2019 novel coronavirus (COVID-19)

A major outbreak of severe acute respiratory syndrome in Hong Kong

Isolation of a novel coronavirus from a man with pneumonia in Saudi Arabia

Characterizing super-spreading in microblog: An epidemicbased information propagation model

Emerging technologies for the detection of viral infections

SARS coronavirus replicase proteins in pathogenesis

The molecular biology of coronaviruses

Efficient assembly and release of SARS coronavirus-like particles by a heterologous expression system

Coronaviruses: an overview of their replication and pathogenesis

Research and Development on Therapeutic Agents and Vaccines for COVID-19 and Related Human Coronavirus Diseases

Identification of novel cellular targets for therapeutic intervention against Ebola virus infection by siRNA screening

Membrane dynamics associated with viral infection

Structural model of the SARS coronavirus E channel in LMPG micelles

Structural insights of a self-assembling 9-residue peptide from the Cterminal tail of the SARS corona virus E-protein in DPC and SDS micelles: A combined high and low resolution spectroscopic study

Targeting viral entry as a strategy for broad-spectrum antivirals

Analysis of therapeutic targets for SARS-CoV-2 and discovery of potential drugs by computational methods

Host cell proteases: Critical determinants of coronavirus tropism and pathogenesis

Absence of E protein arrests transmissible gastroenteritis coronavirus maturation in the secretory pathway

Generation of a replication-competent, propagation-deficient virus vector based on the transmissible gastroenteritis coronavirus genome

Drug targets for corona virus: A systematic review

SARS coronavirus E protein forms cation-selective ion channels

The coronavirus E protein: assembly and beyond

A single polar residue and distinct membrane topologies impact the function of the infectious bronchitis coronavirus E protein

Self-assembly of a nine-residue amyloid-forming peptide fragment of SARS corona virus E-protein: mechanism of self aggregation and amyloidinhibition of hIAPP

Viruses and lipids

Role of lipids in virus replication

Angiotensin-converting enzyme 2 (ACE2) as a SARS-CoV-2 receptor: molecular mechanisms and potential therapeutic target

Coronavirus Spike Protein and Tropism Changes

Mechanisms of coronavirus cell entry mediated by the viral spike protein

NMR structures and localization of the potential fusion peptides and the pretransmembrane region of SARS-CoV: Implications in membrane fusion

NMR structure and localization of a large fragment of the SARS-CoV fusion protein: Implications in viral cell fusion

Mechanistic insights of host cell fusion of SARS-CoV-1 and SARS-CoV-2 from atomic resolution structure and membrane dynamics

Membrane Cholesterol Modulates Oligomeric Status and Peptide-Membrane Interaction of Severe Acute Respiratory Syndrome Coronavirus Fusion Peptide

Structure, Function, and Evolution of Coronavirus Spike Proteins

SARS coronavirus and innate immunity

Endoplasmic Reticulum: The Favorite Intracellular Niche for Viral Replication and Assembly

Subcellular location and topology of severe acute respiratory syndrome coronavirus envelope protein

Coronavirus envelope (E) protein remains at the site of assembly

Membrane lipid therapy: Modulation of the cell membrane composition and structure as a molecular base for drug discovery and new disease treatment

The E protein is a multifunctional membrane protein of SARS-CoV

Exceptional flexibility in the sequence requirements for coronavirus small envelope protein function

Analysis of SARS-CoV E protein ion channel activity by tuning the protein and lipid charge

Biochemical and functional characterization of the membrane association and membrane permeabilizing activity of the severe acute respiratory syndrome coronavirus envelope protein

Biochemical evidence for the presence of mixed membrane topologies of the severe acute respiratory syndrome coronavirus envelope protein expressed in mammalian cells

A highly unusual palindromic transmembrane helical hairpin formed by SARS coronavirus E protein

A balance between membrane elasticity and polymerization energy sets the shape of spherical clathrin coats

Golgi structure formation, function, and post-translational modifications in mammalian cells

Protein folding in the endoplasmic reticulum

Co-and Post-Translational Protein Folding in the ER

The 2018 biomembrane curvature and remodeling roadmap

ER membranes exhibit phase behavior at sites of organelle contact

Membranes: a meeting point for lipids, proteins and therapies

Structure of a conserved Golgi complex-targeting signal in coronavirus envelope proteins

Identification of a Golgi complex-targeting signal in the cytoplasmic tail of the severe acute respiratory syndrome coronavirus envelope protein

The CXXC motif: imperatives for the formation of native disulfide bonds in the cell

Importance of N-linked glycosylation in the functional expression of murine CD1d1

A single asparaginelinked glycosylation site of the severe acute respiratory syndrome coronavirus spike glycoprotein facilitates inhibition by mannose-binding lectin through multiple mechanisms

A Coronavirus E Protein Is Present in Two Distinct Pools with Different Effects on Assembly and the Secretory Pathway

Coronavirus envelope protein: current knowledge

The SARS coronavirus E protein interacts with PALS1 and alters tight junction formation and epithelial morphogenesis

Emerging theme: cellular PDZ proteins as common targets of pathogenic viruses

Viroporin activity of murine hepatitis virus E protein

Viroporins: structure, function and potential as antiviral targets

Relevance of Viroporin Ion Channel Activity on Viral Replication and Pathogenesis. Viruses

Modification of membrane permeability induced by animal viruses early in infection

Expression of SARS-coronavirus envelope protein in Escherichia coli cells alters membrane permeability

Membrane binding proteins of coronaviruses

Multifaceted roles for lipids in viral infection

The cytoplasmic tails of infectious bronchitis virus E and M proteins mediate their interaction

The missing link in coronavirus assembly. Retention of the avian coronavirus infectious bronchitis virus envelope protein in the pre-Golgi compartments and physical interaction between the envelope and membrane proteins

Generation of bovine respiratory syncytial virus (BRSV) from cDNA: BRSV NS2 is not essential for virus replication in tissue culture, and the human RSV leader region acts as a functional BRSV genome promoter

The Virus-Host Interplay: Biogenesis of +RNA Replication Complexes

Biogenesis and dynamics of the coronavirus replicative structures

Modification of intracellular membrane structures for virus replication

Translational control of viral gene expression in eukaryotes

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

Nucleocapsid-independent assembly of coronavirus-like particles by co-expression of viral envelope protein genes

Regulation of transport of the dopamine D1 receptor by a new membrane-associated ER protein

Endoplasmic reticulum export of glycosyltransferases depends on interaction of a cytoplasmic dibasic motif with Sar1

Megavesicles implicated in the rapid transport of intracisternal aggregates across the Golgi stack

Small cargo proteins and large aggregates can traverse the Golgi by a common mechanism without leaving the lumen of cisternae

Qualitative and quantitative ultrastructural analysis of the membrane rearrangements induced by coronavirus

SARS-coronavirus replication is supported by a reticulovesicular network of modified endoplasmic reticulum

Ion channel activity of influenza A virus M2 protein: characterization of the amantadine block

Influenza B virus BM2 protein has ion channel activity that conducts protons across membranes

The transmembrane oligomers of coronavirus protein E

Conductance and amantadine binding of a pore formed by a lysine-flanked transmembrane domain of SARS coronavirus envelope protein

Importance of conserved cysteine residues in the coronavirus envelope protein

Structure and inhibition of the SARS coronavirus envelope protein ion channel

Coronavirus E protein forms ion channels with functionally and structurally-involved membrane lipids

Channel-Inactivating Mutations and Their Revertant Mutants in the Envelope Protein of Infectious Bronchitis Virus

Pharmaceutical Targeting the Envelope Protein of SARS-CoV-2: the Screening for Inhibitors in Approved Drugs

Multitude NMR studies of α-synuclein familial mutants: probing their differential aggregation propensities

Monitoring peripheral protein oligomerization on biological membranes

Membrane domain modulation of Aβ

Biophysical insights into how surfaces, including lipid membranes, modulate protein aggregation related to neurodegeneration

Architectural and mechanistic insights into an EHD ATPase involved in membrane remodelling

A hinge in the distal end of the PACSIN 2 F-BAR domain may contribute to membrane-curvature sensing

Membrane curvature at a glance

Constitutively oxidized CXXC motifs within the CD3 heterodimeric ectodomains of the T cell receptor complex enforce the conformation of juxtaposed segments

Furin cleavage potentiates the membrane fusion-controlling intersubunit disulfide bond isomerization activity of leukemia virus Env

Intersubunit disulfide isomerization controls membrane fusion of human T-cell leukemia virus Env

Formation of transitory intrachain and interchain disulfide bonds accompanies the folding and oligomerization of simian virus 40 Vp1 in the cytoplasm

Assessing the role of aromatic residues in the amyloid aggregation of human muscle acylphosphatase

A possible role for pi-stacking in the self-assembly of amyloid fibrils

Infectious bronchitis virus E protein is targeted to the Golgi complex and directs release of virus-like particles

Physical interaction between the membrane (M) and envelope (E) proteins of the coronavirus avian infectious bronchitis virus (IBV)

Analysis of constructed E gene mutants of mouse hepatitis virus confirms a pivotal role for E protein in coronavirus assembly

Coronavirus pseudoparticles formed with recombinant M and E proteins induce alpha interferon synthesis by leukocytes

Membrane Lipid Composition: Effect on Membrane and Organelle Structure, Function and Compartmentalization and Therapeutic Avenues

Protein aggregation in a membrane environment

Protein aggregation and membrane lipid modifications under lactic acid stress in wild type and OPI1 deleted Saccharomyces cerevisiae strains

Biological phase separation: cell biology meets biophysics

Membrane lipids: where they are and how they behave

SARS coronavirus E protein in phospholipid bilayers: an x-ray study

Membrane topology and insertion of membrane proteins: search for topogenic signals

Another triple-spanning envelope protein among intracellularly budding RNA viruses: the torovirus E protein

Genomic characterisation and epidemiology of 2019 novel coronavirus: implications for virus origins and receptor binding

Identifying SARS-CoV-2

Coevolution in defining the functional specificity

Predicting transmembrane protein topology with a hidden Markov model: application to complete genomes

Improving the accuracy of transmembrane protein topology prediction using evolutionary information

Advantages of combined transmembrane topology and signal peptide prediction--the Phobius web server

Analysis of membrane and surface protein sequences with the hydrophobic moment plot

Trends in substitution models of molecular evolution

The rapid generation of mutation data matrices from protein sequences

MEGA X: Molecular Evolutionary Genetics Analysis across Computing Platforms

Membrane topology of coronavirus E protein

Characterization of the coronavirus mouse hepatitis virus strain A59 small membrane protein E

Two deletion variants of Middle East respiratory syndrome coronavirus found in a patient with characteristic symptoms

Coronavirus pathogenesis and the emerging pathogen severe acute respiratory syndrome coronavirus

Coronavirus virulence genes with main focus on SARS-CoV envelope gene

Genome-wide analysis of protein-protein interactions and involvement of viral proteins in SARS-CoV replication

Hexamethylene amiloride blocks E protein ion channels and inhibits coronavirus replication

2020. approaches to detect inhibitors of the human severe acute respiratory syndrome coronavirus envelope protein ion channel

Antiviral activity of cepharanthine against severe acute respiratory syndrome coronavirus in vitro

Drug repurposing strategies for COVID-19

Recent progress and challenges in drug development against COVID-19 coronavirus (SARS-CoV-2) -an update on the status

TGEV corona virus ORF4 encodes a membrane protein that is incorporated into virions