- 1Division of Medicinal Chemistry and Pharmacognosy, College of Pharmacy, The Ohio State University, Columbus, OH, United States
- 2Department of Chemical and Biomolecular Engineering, College of Engineering, The Ohio State University, Columbus, OH, United States
- 3Department of Chemical and Environmental Engineering, The University of Cincinnati, Cincinnati, OH, United States
- 4Translational Data Analytics Institute (TDAI), The Ohio State University, Columbus, OH, United States
COVID-19 has become an unprecedented threat to human health. The SARS-CoV-2 envelope (E) protein plays a critical role in the viral maturation process and pathogenesis. Despite intensive investigation, its structure in physiological conditions remains mysterious: no high-resolution full-length structure is available and only an NMR structure of the transmembrane (TM) region has been determined. Here, we present a refined E protein structure, using molecular dynamics (MD) simulations to investigate its structure and dynamics in a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) bilayer system. Our initial homology model based upon the SARS-CoV E protein structure is shown to be unstable in the lipid bilayer, and the H3 helices tend to move away from the membrane center to the membrane-water interface. A more stable model was developed by replacing all H3 helices with the fully equilibrated H3 structure sampled in the MD simulations. This refined model exhibited more favorable contacts with lipids and water than the original homology model and induced local membrane curvature, decreasing local lipid order. Interestingly, the pore radius profiles showed that the channel in both homology and refined models remained in a closed state throughout the simulations. We also demonstrated the utility of this structure to develop anti-SARS-CoV-2 drugs by docking a library of FDA-approved, investigational, and experimental drugs to the refined E protein structure, identifying 20 potential channel blockers. This highlights the power of MD simulations to refine low-resolution structures of membrane proteins in a native-like membrane environment, shedding light on the structural features of the E protein and providing a platform for the development of novel antiviral treatments.
Introduction
Over 2 years have passed since the outbreak of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), which is also referred to as the novel coronavirus disease (COVID-19). As of January 2022, there have been over 300 million confirmed cases of COVID-19 including over five million deaths. Thanks to the strict social distancing in many countries and the rollout of vaccines from Pfizer/BioNtech, Moderna, and Johnson & Johnson, etc., the spread of the virus has been tremendously slowed down. However, the emergence of several SARS-CoV-2 variants, such as the Delta and Omicron variants has caused renewed COVID outbreaks, increasing acute concerns and reminding the world to maintain efforts to find new and better therapeutics and antiviral agents.
The genome of SARS-CoV-2 encodes three membrane proteins: the spike (S) protein, the membrane (M) protein, and the envelope (E) protein. S protein binds to the receptor on host cell surface and mediates viral entry (Zhou et al., 2020; Jackson et al., 2022). M protein plays a critical role in virus assembly and budding processes (Weiss and Navas-Martin, 2005; Plescia et al., 2021). E protein is composed of 75 amino acids and interacts with M protein, which is important for viral formation and release (Lim and Liu, 2001; Hsieh et al., 2008). The E protein is an integral membrane protein. It forms a pH-sensitive pentameric cation channel on the ERGIC/Golgi membrane, which is permeable to Ca2+, Na+, and K+ (Verdiá-Báguena et al., 2012; Nieto-Torres et al., 2015; Cabrera-Garcia et al., 2021; Verdiá-Báguena et al., 2021; Xia et al., 2021). Deletion of the E protein in recombinant coronaviruses resulted in significantly reduced viral titers and viral maturation, while mutations of the E protein affected post-translational modification and impaired ion conductivity, leading to reduced viral formation and pathogenicity (DeDiego et al., 2007; Boscarino et al., 2008; Lopez et al., 2008; Verdiá-Báguena et al., 2012; DeDiego et al., 2014; Castaño-Rodriguez et al., 2018). The E protein, together with the M protein, modulates the N-glycosylation of the S protein and is essential for the optimal assembly of viral particles (Siu et al., 2008; Boson et al., 2021). The E protein alone was found to cause cell death in vitro and trigger strong immune responses in vivo (Xia et al., 2021). A recent study demonstrated that the SARS-CoV-2 E protein stimulates the production of an inflammatory chemokine by binding and activating the TLR2 signaling cascade (Planès et al., 2022). Thus, it is hypothesized that inhibition of the E protein’s ion channel activity could disrupt viral assembly, decrease viral maturation, and eventually reduce viral infection. Several studies have identified potential E protein viroporin blockers, some of which exhibited protective effects in vitro against SARS-CoV/SARS-CoV-2, such as hexamethylene amiloride, amantadine, glicazide, and tretinoin, etc (Pervushin et al., 2009; Dey et al., 2020; Mandala et al., 2020; Singh Tomar and Arkin, 2020; Das et al., 2021; Park et al., 2021; Xia et al., 2021). Compared to the heavily studied S protein which has exhibited rapid mutation that often correlates with enhanced infection and death rates, the E protein appears to be more stable and conserved among detected variants. The E protein is also less thoroughly investigated as a therapeutic target for vaccine and drug development against coronavirus (Alam et al., 2020; Tilocca et al., 2020; Das and Roy, 2021; Troyano-Hernáez et al., 2021). The stability of the E protein may point to it being a better long-term target for therapeutic intervention against SARS-CoV-2.
The structure of the E protein can be three regions: the N-terminal domain (NTD), the transmembrane domain (TMD), and the C-terminal domain (CTD) (Figure 1A) (Schoeman and Fielding, 2019). The TM regions oligomerize into a homopentameric channel (Pervushin et al., 2009; Surya et al., 2018; Mandala et al., 2020). Mutations in the hydrophobic TMD, such as N15A and V25F, led to the abrogation of ion conductance (Verdiá-Báguena et al., 2012). The CTD, facing the cytoplasm, contains a PDZ-binding motif (final four amino acid residues DLLV), which can potentially bind to over 400 host cellular proteins that carry a PDZ domain, like PALS1 or ZO1 (Teoh et al., 2010; Münz et al., 2012; Castaño-Rodriguez et al., 2018; Chai et al., 2021; Shepley-McTaggart et al., 2021). Although the SARS-CoV-2 E protein has been intensively studied since the outbreak of COVID-19, its structure remained largely unknown until a high-resolution solid-state NMR structure of its TM region was reported in late 2020 (PDB code: 7K3G) (Mandala et al., 2020). A full-length (PDB code: 2MM4) and a near full-length (residues 8-65) (PDB code: 5X29) solution NMR structures of its close analog, the SARS-CoV E protein, were determined in sodium dodecyl sulfate (SDS) and 1-myristoyl-2-hydroxy-sn-glycero-3-phospho-(1′-rac-glycerol) (LMPG) micelles, respectively, which could serve as good templates for homology modeling of the SARS-CoV-2 E protein (Li et al., 2014; Surya et al., 2018).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1027223/fmolb-09-1027223-HTML/image_m/fmolb-09-1027223-g001.jpg)
FIGURE 1. (A) Sequence alignment of SARS-CoV-2 and SARS-CoV E protein. Non-conserved residues are highlighted in red. (B) Superposition of the homology model of SARS-CoV-2 E protein (green) and the NMR structure (PDB code: 5X29) of SARS-CoV E protein (cyan). (C) Assembled pentameric model of SARS-CoV-2 E protein. (D) Top view and (E) side view of the SARS-CoV-2 E protein embedded in a POPC bilayer. Chain A-E are colored with red, blue, green, magenta and orange, respectively. Lipids are shown as blue lines, with the phosphorus atoms of the lipid shown as olive spheres.
Here, we built a homology model of SARS-CoV-2 E protein and applied microsecond-long all-atom molecular dynamics (AA-MD) simulations to study the structural and dynamic properties of the SARS-CoV-2 E protein in a model 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) bilayer. Our study indicated that MD simulations can be employed to refine NMR structures determined from non-native environments. The refined E protein structure allowed us to explore the interplay between the E protein and the membrane in which it natively exists. The refined structure also provided a structural basis to start identifying potential channel blockers by molecular docking of a library of FDA-approved, investigational and experimental drugs against the E protein.
Results and discussion
SARS-CoV-2 E protein was constructed through homology modeling and assembled into a lipid bilayer
The full-length structure of SARS-CoV-2 E protein remains ambiguously defined despite numerous investigations. A partial solid-state NMR structure of the protein, the TM domain, was determined by Hong’s group (PDB code: 7K3G) (Mandala et al., 2020), while Surya et al. reported a solution NMR structure of the SARS-CoV E protein in LMPG micelles (PDB code: 5X29), which covered residues 8-65 (Surya et al., 2018). Since the overall and TMD of the SARS-CoV and SARS-CoV-2 E proteins share 94.7% and 100% sequence identity, respectively (Figure 1A), we built a homology model of the near full-length SARS-CoV-2 E protein (residues 8-65) with MODELLER (Šali and Blundell, 1993) using the 5X29 structure as the template (Figure 1B). The homology model has the same topology as the template: a transmembrane α-helix (residue 8-38) followed by two short α-helices (H2 and H3) connected by a flexible linker in the CTD (residue 39-65). The stereochemical quality of the homology model was validated by SAVES v6.0 developed by DOE lab at UCLA (https://saves.mbi.ucla.edu/) and Molprobity developed by Williams et al. at Duke University (http://molprobity.biochem.duke.edu/) (Williams et al., 2018) and the results were summarized in Supplementary Table S1. All residues in our homology model fall in the allowed regions and the ERRAT overall quality factor is 88 with a 3.01 clash score, indicating that the homology model has an acceptable quality. Interestingly, the two NMR structures have a relatively low ERRAT overall quality factor and 3.6% of 7K3G’s residues even fall in the disallowed regions. The other three homology models all have comparable or slightly better qualities than our homology model. Thus, we think that our homology model can be reliably used as the starting structure for MD simulations.
The homology model was assembled into a pentameric complex (Figure 1C) and subsequently inserted into a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) bilayer using the CHARMM GUI webserver (Figures 1D,E) (Jo et al., 2008; Wu et al., 2014). However, visual inspection of the assembled system suggests that some polar and charged residues in the CTDs of the E protein are buried in the lipid bilayer and thus have unfavorable interactions with the hydrophobic lipid tails. This unphysical arrangement of the CTDs might result from the non-native condition in which the SARS-CoV E protein was determined - detergent micelles - instead of more native-like phospholipid bicelles or nanodiscs (Lau et al., 2009; Suk et al., 2012; Surya et al., 2013). Thus, we aimed to refine the homology model of SARS-CoV-2 E protein by all-atom molecular dynamics (AA-MD) simulations. Three separate simulations were performed for the E protein homology model embedded in a POPC bilayer at 303K for 2 µs each. An additional simulation of the same system was performed at 323 K for 1 µs in order to speed up the conformational sampling.
SARS-CoV-2 E protein displays significant structural fluctuations in AA-MD simulations
During the simulations, the E protein (referred to as the homology model below) displays significant structural fluctuations as indicated by the large root mean square deviation (RMSD) values of the heavy atoms, quickly rising to ∼7 Å within the first 50 ns (Figure 2A). The CTDs contribute more to the structural fluctuations as indicated by the larger RMSDs (>5 Å) than the TMs (∼3Å) throughout the simulations (Supplementary Figures S1A, S1C). A similar fluctuation pattern is observed in the root mean square fluctuations (RMSFs) of the Cα atoms, where the RMSFs of CTDs are much higher than those of the TM regions (Supplementary Figure S2A). Within the CTDs, the H3 helices fluctuate more than the H2 helices (Figures 2B,C). During the simulations, all five H2s stay at the membrane-water interface, but the H3s of four monomers (chain A, B, C and E) move away from the hydrophobic core toward the lipid headgroups while only one H3 (chain D) does not display any significant upward or downward movement and remains at the lipid tailgroup level (Figures 2D,E). This movement is likely due to the incompatibility between the polar and charged residues in the CTDs and the hydrophobic lipid core, suggesting that the NMR structure determined in detergent micelles does not represent a functionally relevant conformation of the E protein in its native environment.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1027223/fmolb-09-1027223-HTML/image_m/fmolb-09-1027223-g002.jpg)
FIGURE 2. (A) RMSDs of the homology model and the refined model. The lines are the mean values averaged from the three parallel simulations. The shade area is the standard deviation. RMSDs of (B) the H2s and (C) the H3s of the homology model. (D) The density profile of the centers of mass (COMs) of H3s along the membrane normal in the simulations of the homology model. The average position of the phosphorus atoms on the upper leaflet is shown as a dashed line for reference. (E) The initial (left) and last frame (right) of the homology model in a POPC bilayer from the AA-MD simulations. H3 helices of chain A-E are colored with red, blue, green, magenta and orange, respectively. Lipids are shown as blue lines and the phosphorus atoms are shown as olive spheres.
A refined SARS-CoV-2 E protein model shows smaller structural fluctuations
Simulation of the homology model of the SARS-CoV-2 E protein showed that the homology modeling approach did not result in a realistic conformation in the membrane environment. Thus, we built a new structural model of the SARS-CoV-2 E protein by taking a structural frame at ∼500 ns from the 323K simulation trajectory, in which the H3 of its chain C has completely moved to the membrane-water interface. We replicated this chain C model four times and assembled the resulting five protomers into a pentamer using the initial homology model as a template. The new pentameric model (referred to as the refined model below) was subsequently inserted into the POPC bilayer (Figure 3A). In the refined model, all H3s were located at the membrane-water interface, with H2s partially unfolded. Three MD simulations of the refined model were performed at 303K for 2 µs each. The overall RMSD values of the refined model are smaller than those of the homology model (Figure 2A). Particularly, compared to the original homology model, the refined model displays smaller RMSDs and RMSFs for both TMDs and CTDs (both H2s and H3s) (Figures 3B, 3C; Supplementary Figures S1B, S1D, S2B), indicating that the new model is more stable. This is likely due to the elimination of the unfavorable interactions between CTDs and lipid tails in the homology model. H2s stay in the water layer while H3s remain at the membrane-water interface for most of the time throughout the simulations (Figures 3D,E), only occasionally dissociating and re-associating from the interface area, indicated by the sudden increase of RMSD and the movement of the centers of mass (COMs) along the membrane normal (Figure 3C; Supplementary Figure S3B). Interestingly, the TM helices in the refined model are found to tilt slightly more than those in the homology model, which might also facilitate the H3s approaching the membrane headgroup region (Supplementary Figure S4).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1027223/fmolb-09-1027223-HTML/image_m/fmolb-09-1027223-g003.jpg)
FIGURE 3. (A) The refined model assembled in a POPC bilayer (side view). RMSDs of (B) the H2s and (C) the H3s of the refined model. (D) The density profile of the H3s COMs along the membrane normal in the simulations of the refined model. The average position of the phosphorus atoms on the upper leaflet is shown as a dashed line for reference. (E) The initial (left) and last frame (right) of the refined model in a POPC bilayer from the AA-MD simulations. The representation and color scheme are the same as Figure 2E.
We next calculated the α-helical probability for each residue from the AA-MD simulations. The secondary structure content of the TM domains is largely maintained in both models during the simulations, while the TM regions in the refined model show a slightly higher α-helical content than those in the homology model (∼80% vs. ∼70%) (Supplementary Figure S5). In contrast, the H2 helices in the refined model are found to be less structured than those in the homology model. This is probably because all H2 helices in the refined model are fully immersed in the water layer, which caused them to be at fast unfolding/refolding equilibrium, while those helices in the homology model are only partially above the membrane surface and remain largely folded in an amphiphilic environment (Figures 4A,B). Unlike the H2 helices, the H3 helices in the refined model show a slight increase of the α-helix content when comparted to the homology model. This may be because these helices can make more favorable interactions with lipid headgroups in the refined model than in the homology model where the H3s are trapped in a hydrophobic environment and are thus energetically less stable (Figures 4C,D). These results are in good agreement with a recently published study where the authors demonstrated the secondary structure of the TM regions is mostly maintained while the H2s and H3s are more disordered throughout the atomistic MD simulations (Kuzmin et al., 2022).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1027223/fmolb-09-1027223-HTML/image_m/fmolb-09-1027223-g004.jpg)
FIGURE 4. Time-evolution of α-helical content of H2s and H3s in (A) the homology and (B) the refined model. Average probability of the α-helical structure for each residue from simulations of (C) the homology model and (D) the refined model. Cyan: TM regions; green: H2 helices; magenta: H3 helices.
The refined model is stabilized by more H bonds and contacts with lipids/water
To understand how the refined model is stabilized relative to the homology model, we examined the molecular interactions between the protein complex and lipid or water molecules in the local environment. In both models, the CTDs establish more H bonds with lipids and water than the TM regions (Figures 5A–D). The CTDs in the refined model form more H bonds (∼30) with lipids than in the homology model (∼20) because the CTDs in the refined model were in closer proximity to lipid headgroups. Figure 5E depicts the number of contacts with lipids for individual residues of the E protein, which confirms that the majority of residues involved in the protein-lipid interactions is located in the CTD regions and the H3s form more contacts with lipids in the refined model than in the homology model. Strikingly, the TM domains in the refined model form fewer H bonds with water than those in the homology model, while the number of H bonds between the TM regions and lipids is approximately equal for both models. Close inspection reveals a substantial amount of water molecules in the homology model moving into the lipid bilayer and staying in the vicinity of the protein. This appears to be due to the burial of the polar H3s in the membrane hydrophobic core which greatly destabilizes the lipid bilayer. In contrast, the lipid bilayer remains intact with no water molecules trapped inside the membrane throughout the simulations of the refined model (Supplementary Figure S6). These data show that the integrity of the membrane is dependent on the presence of the different protein models, substantiating that the refined model is likely a more physiologically-relevant conformation of the E protein in the membrane environment than the original homology model.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1027223/fmolb-09-1027223-HTML/image_m/fmolb-09-1027223-g005.jpg)
FIGURE 5. (A–D) H bonds formed between lipids/water and the homology and refined models. (E) Total number of contacts with lipids per residue (within a 4 Å cut-off) in both models.
E protein induces local membrane curvature
We next investigated how the presence of the E protein affects various membrane properties. We first calculated the bilayer thickness as a function of the distance from the COM of the E protein (Figure 6A). Compared to a pure bilayer system, the bilayer thickness decreases within ∼3 nm of the protein COM in both models, while the presence of the protein seems to have little or no effects on the bilayer thickness between 4 and 6 nm from the protein COM (the largest radius of the protein is ∼2.5 nm), suggesting that the bilayer around the E protein bends toward the membrane core. It is noteworthy that, along with the movement of the E protein, especially H3s, this membrane bending motion also contributes to bringing the H3 helices close to the membrane-water interfacial region. These results are in agreement with Kuzmin et al.‘s work, where the authors found that the induction of membrane curvature by the E protein led to thinner membrane close to the TM regions but thicker around the CTDs (Kuzmin et al., 2022). Yet, another study by Collins et al. found that the E protein was not able to facilitate membrane curvature (Collins et al., 2021). These inconsistent results might arise from the different protein structural models (Heo and Feig’s predicted model (Heo and Feig, 2020) vs. our refined model generated from MD simulations) or the different membrane systems used in the studies (ERGIC mimic vs. pure POPC).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1027223/fmolb-09-1027223-HTML/image_m/fmolb-09-1027223-g006.jpg)
FIGURE 6. (A) Membrane thickness and (B) order parameter as a function of the distance from the COM of the E protein. The dashed lines indicate the corresponding values calculated from the AA-MD simulation (200 ns) of a pure POPC bilayer. (C) STilt and (D) SSplay as a function of the distance from the COM of the E protein. STilt is calculated as the average tilt angle of a lipid with respect to the bilayer normal. SSplay reflects the divergence of local lipids tilt.
In parallel to bilayer thinning, a decrease in the order parameter (SCH) of lipid acyl chains is found in both models when compared to the pure POPC bilayer, especially for lipids within 1-4 nm of the protein. This indicates that the protein makes the surrounding lipids more flexible and less ordered (Figure 6B). We also observed that both models decrease the tilt (STilt) and splay (SSplay) parameters of surrounding lipids (Figures 6C,D). Since lipid tilt and splay are known to be correlated with membrane local rigidity (Watson et al., 2012; Khelashvili et al., 2013), these results suggest that the E protein can soften local lipid bilayer.
Taken together, our simulations show interplay between SARS-CoV-2 E protein and the membrane, which rearranges the position/orientation of the H3 helices relative to the membrane, eventually leading to more favorable accommodation of the protein in the membrane. This rearrangement is achieved through a mutual adaptation mechanism. The E protein relieves frustration due to the burial of its H3s in the membrane, through the movement of the H3 helices and the shift/tilting of the TM helices, while the membrane adapts to the protein by changing its thickness and lipid order. The E protein is known to play a critical role in the viral maturation and budding process of SARS-CoV-2 infection. However, the exact mechanism remains poorly understood. Our results reveal that the presence of the protein induces a local thinning and softening of the surrounding lipid bilayer, which in turn facilitates accommodating the TM helices comprising polar residues on both termini that need to be solvent-exposed. This window into the complexity of the E protein-membrane interaction provides an additional basis for understanding how the E protein contributes to the viral budding process.
The viroporin channel remains in a closed state in both models
The ion channel of the E protein is lined primarily with hydrophobic residues (e.g., LEU18, ALA22, PHE26, VAL29, ILE33 and LEU37) and two hydrophilic residues (i.e., GLU8 and ASN15). We calculated the channel radius profiles of both models (Figures 7A,B). Since the TM domains are relatively stable during the simulations, the pore radius profiles of both models are similar to each other. PHE26 is identified as the bottleneck (narrowest) residue in the middle of the pore, which could mediate the opening and closure of the channel, in line with several previous studies (Figures 7C,D) (Sarkar and Saha, 2020; Monje-Galvan and Voth, 2021). Both models remain in the closed conformation, and the transition from the closed to the open state is not observed throughout the simulations. While the channel is not well-hydrated in the refined model, the hydration profile of the homology model seems to indicate a well-hydrated channel. However, as discussed above, the hydration of the homology model is not due to water molecules passing through the channel but rather being trapped inside the membrane (Figures 7E,F).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1027223/fmolb-09-1027223-HTML/image_m/fmolb-09-1027223-g007.jpg)
FIGURE 7. Time-averaged pore radius profiles of (A) the homology model and (B) the refined model. The grey shaded areas correspond to the one-sigma (dark gray) and two-sigma (light gray) confidence intervals. The pore lining residues (facing the channel more than 50% of the time) are indicated as dots colored by their hydrophobicity. Top view of the bottleneck residue PHE26 in the channel of (C) the homology model and (D) the refined model. Pore hydration profiles of (E) the homology model and (F) the refined model. Coordinates of water molecules along the z axis were plotted as a function of time. Each water molecule was represented by a blue dot.
Repurposing drugs targeting the E protein through virtual screening
The ion channel formed by the pentameric E protein is critical to the viral assembly and maturation (Verdiá-Báguena et al., 2012; Nieto-Torres et al., 2015; Cabrera-Garcia et al., 2021). Impaired ion conductivity caused decreased virulence and reduced viral infection (Boscarino et al., 2008; Castaño-Rodriguez et al., 2018; Xia et al., 2021). Several studies sought to identify the E protein channel blockers by in vitro cell-based assays and in silico screening (Dey et al., 2020; Das et al., 2021; Mukherjee et al., 2021; Xia et al., 2021). Given a better equilibrated and likely more physiologically relevant structural model of the SARS-CoV-2 E protein, we conducted virtual screening of a library of FDA-approved, investigational and experimental drugs to discover potential E protein ion channel inhibitors. The drug library containing 9,213 compounds was obtained from the DrugBank and screened in silico by Schrodinger Glide on the Ohio Supercomputer Center (Ohio-Supercomputer-Center, 1987; Halgren et al., 2004; Wishart et al., 2017). The docking results were ranked by Glide score (gscore) and the top 100 results were visually inspected. The entire channel can be divided by the bottleneck residue PHE26 into an upper and a lower channel. All the docking poses are located in the lower channel. Many of the top-ranking compounds form H bonds with ASN15 side chain at the bottom of the channel. The mutation of ASN15, such as N15A, was found to reduce ion conductance (Verdiá-Báguena et al., 2012). Additionally, some top-ranking compounds form H bonds with the backbone of LEU18. Only one molecule (DB11262) establishes π-π stacking interaction with PHE26. 20 compounds were selected from the top-ranking results based on the key interactions formed with the E protein (Table 1) and their physiochemical properties and mode of action/indications were summarized in Supplementary Tables S2, S3. The representative binding poses of the selected compounds are shown in Figure 8; Supplementary Figure S8. We found that most of these compounds are beyond the rule-of-five scope, suggesting that in order to fully block the channel, large compound size seems necessary. Interestingly, among these twenty compounds, two molecules (DB12138 and DB11871) might also inhibit the SARS-CoV-2 main protease according to the screening work by Sharma et al. (Sharma et al., 2021). These results might provide an alternative starting point for the development of therapeutics against COVID-19.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1027223/fmolb-09-1027223-HTML/image_m/fmolb-09-1027223-t001.jpg)
TABLE 1. List of 20 top-ranked compounds and their Glide score (gscore) and the key interactions formed with the E protein.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1027223/fmolb-09-1027223-HTML/image_m/fmolb-09-1027223-g008.jpg)
FIGURE 8. Docking poses of (A) DB03005, (B) DB06942, (C) DB04042 and (D) DB11262. For clarity, one of the chains of the E protein is not shown.
Conclusion
The SARS-CoV-2 envelope (E) protein is thought to play essential roles in the viral maturation and pathogenesis, but its structure in physiological conditions is not well understood. This has greatly hampered our fundamental understanding of E protein-mediated membrane fusion and viral assembly. This study began with a homology model of the E protein and investigated its structure and dynamics in a POPC lipid bilayer with molecular dynamic (MD) simulations. Our results show that the homology model, especially the C-terminal domain, is not stable when embedded in the lipid bilayer. The H3 helices of four monomers are found to move to the lipid headgroup region in the simulations. Based on this observation, we built a refined E protein structure by replacing the H3 structure with the one fully equilibrated in the MD simulations. Not surprisingly, the refined model displays smaller structural fluctuations. The stability of the refined model likely arises from the more H bonds and favorable interactions of the H3 helices with both lipids and water.
The presence of the SARS-CoV-2 E protein induces a local membrane curvature - thinning and softening the surrounding lipid bilayer. This could partially facilitate the accommodation of the unstable H3 helices in the original homology model. Furthermore, our simulations indicate that both models represent a closed channel with a hydrophobically occluded pore. Finally, given a better E protein structural model, we conducted virtual screening against a library of 9,213 FDA drugs and discovered 20 compounds that could potentially function as E protein channel blockers.
The refined structural model for the SARS-CoV-2 E protein is offered as a platform to aide in the interpretation of experimental and computational results as the community continues to investigate the role of the E protein viral maturation and pathogenesis. As we demonstrate, it is also expected that this refined model will facilitate the development anti-COVID-19 drugs targeting the E protein.
Methods
Homology modeling and system preparation
The primary sequences of SARS-CoV E protein (Accession number: AYV99820.1) and SARS-CoV-2 E protein (Accession number: QII57162.1) were retrieved from NCBI protein database at www.ncbi.nlm.nih.gov/protein/. These two proteins have high similarity (94.74% identity) in sequence, which were compared at http://blast.ncbi.nlm.nih.gov. The structure of SARS-CoV E protein (PDB 5X29) was used as the template structure for homology modeling by Modeller 9.1 (Šali and Blundell, 1993). The sequence of the SARS-CoV-2 E protein was aligned with the structure of the SARS-CoV E protein using the align2d command in MODELLER. 10,000 3D models of the SARS-CoV-2 E protein were generated using the AutoModel class. The results were plotted by their DOPE score and molpdf values. Three structures with the lowest DOPE score and molpdf values (query.B99990502, query.B99990876, and query.B99990918) were visually inspected and query.B99990502 was selected as the final model due to its smallest structural deviation from the template. The model was submitted to SAVES v6.0 developed by the DOE lab at UCLA (https://saves.mbi.ucla.edu/) and Molprobity developed by Williams et al. at Duke University (http://molprobity.biochem.duke.edu/) (Williams et al., 2018) for stereochemical quality check. The final pentameric structure of SARS-CoV-2 E protein was built with the homology modeling structure by GalaxyHomomer (https://galaxy.seoklab.org/cgi-bin/submit.cgi?type=HOMOMER). Each pentameric E protein structure was embedded in a membrane bilayer containing 150 1-palmitoyl-2-oleoyl-glycero-3-phosphocholine (POPC) molecules on each leaflet using CHARMM-GUI Membrane Builder (Jo et al., 2008; Wu et al., 2014). The bury depth was determined by PPM (one of the orientation options on CHARMM-GUI), which predicts the membrane embedding orientation and depth for a given input protein structure at https://opm.phar.umich.edu/ppm_server2_cgopm. Approximately 20,000 TIP3P water molecules and 0.15 M NaCl were added to each system to ensure electric neutrality.
AA-MD simulations
All simulations were performed using AMBER following the 6-step protocol provided on the CHARMM-GUI web-server (D.A. Case, 2018). The Lipid14 force field was used for lipids (Dickson et al., 2014). The ff14SB was used for the protein (Maier et al., 2015). The TIP3P was used for water molecules (Price and Brooks, 2004). The AMBER-99 force field was used for ions (Chen and Pappu, 2007). The non-bonded pair list was updated every 1,000 steps with the cutoff of 12 Å. Long-range electrostatic interactions were calculated using particle mesh Ewald (PME) method (Darden et al., 1993). Bonds involving hydrogen were constrained with the SHAKE algorithm (Ryckaert et al., 1977). The temperature was maintained at 303.15 K with a friction coefficient of 1.0 ps−1 using the Langevin dynamics algorithm and the pressure was maintained at 1.0 bar using the Berendsen pressure coupling algorithm (Berendsen et al., 1984). The pressure was coupled semi-isotropically, in which the x and y directions were coupled together, and the z direction fluctuates independently. Three independent 2 μs simulations were run for the homology and the refined models, respectively. One simulation of the homology model was run at 323.15 K for 1 μs. One simulation of a pure lipid bilayer consisting of 300 POPC molecules (150 on each leaflet) was run at 303.15K for 200 ns.
Virtual screening
A compound library containing 9,213 FDA-approved, investigational and experimental drugs was obtained from the DrugBank (Wishart et al., 2017). Each of these compounds was docked into the E protein channel by Glide XP (Madhavi Sastry et al., 2013) from the Schrodinger suite installed on the Ohio Supercomputer Center (Ohio-Supercomputer-Center, 1987). Ligand structures were prepared by Ligprep (Schrödinger, 2021a). The protein structure was prepared by Protein Preparation Wizard (Madhavi Sastry et al., 2013). The docking results were analyzed based on docking scores and visual inspection. The images of the docking poses were rendered by Maestro (Schrödinger, 2021b).
Analysis
All the images, structures and trajectories were visualized and rendered on the Visual Molecular Dynamics (VMD) software package (Humphrey et al., 1996) and Pymol (Schrodinger, 2015) except Figures 7A,B which were generated by The Channel Annotation Package (Klesse et al., 2019). All figures were plotted with OriginPro (OriginLab Corporation, Northampton, MA, United States). All analyses were performed with CPPTRAJ or Python using in-house scripts (Roe and Cheatham, 2013).
Membrane thickness is calculated using the positions of phosphate atoms of the lipid head. The upper and lower surface are represented by the two-dimensional Fourier series to the order of 3 in both x and y direction. Coefficients are determined by fitting the series to the phosphate atoms positions for every frame (the first 500ns are removed).
For each phosphate atom, the thickness at that position is calculated by
The order parameter is defined as:
where α is the angle between the C-H bond and the bilayer normal.
Tilt is defined as the angle between the bilayer normal and a vector that connects the COM of the phosphate and glyceride C2 atom with the COM of the three terminal carbon atoms on the lipid tail of POPC (Khelashvili et al., 2013). The splay angle of a lipid is the average angle of its head-to-tail vector with that of surrounding lipids within 10 Å. The STilt and SSplay are derived analogously from Eq. 2 for the tilt and splay angles, respectively.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation. The homology model of the E protein, the refined model of the E protein, and computational models of the E protein bound with drug molecules are available at https://u.osu.edu/chenglab/research/.
Author contributions
XC and JN conceived the ideas. RY performed the simulations and did most of the analyses. SW assisted in constructing the homology model. SW and GR did some of the analyses. RY, XC, and JN wrote the manuscript and all authors reviewed and approved the manuscript for publication.
Funding
This work was partially supported by the start-up fund from the College of Pharmacy at The Ohio State University (to XC). This work has also used resources from the Ohio Supercomputer Center (OSC) in support of COVID-19-related research.
Acknowledgments
This research was supported by the startup fund from the College of Pharmacy and the Discovery Themes at The Ohio State University (OSU) and OSU’s Translational Data Analytics Institute (TDAI) Interdisciplinary Research Pilot Award. This research used the resources of Ohio Supercomputer Center (OSC).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmolb.2022.1027223/full#supplementary-material
References
Alam, I., Kamau, A. A., Kulmanov, M., Jaremko, Ł., Arold, S. T., Pain, A., et al. (2020). Functional pangenome analysis shows key features of E protein are preserved in SARS and SARS-CoV-2. Front. Cell. Infect. Microbiol. 10, 405. doi:10.3389/fcimb.2020.00405
Berendsen, H. J. C., Postma, J. P. M., Gunsteren, W. F. v., DiNola, A., and Haak, J. R. (1984). Molecular dynamics with coupling to an external bath. J. Chem. Phys. 81 (8), 3684–3690. doi:10.1063/1.448118
Boscarino, J. A., Logan, H. L., Lacny, J. J., and Gallagher, T. M. (2008). Envelope protein palmitoylations are crucial for murine coronavirus assembly. J. Virol. 82 (6), 2989–2999. doi:10.1128/JVI.01906-07
Boson, B., Legros, V., Zhou, B., Siret, E., Mathieu, C., Cosset, F.-L., et al. (2021). The SARS-CoV-2 envelope and membrane proteins modulate maturation and retention of the spike protein, allowing assembly of virus-like particles. J. Biol. Chem. 296, 100111. doi:10.1074/jbc.RA120.016175
Cabrera-Garcia, D., Bekdash, R., Abbott, G. W., Yazawa, M., and Harrison, N. L. (2021). The envelope protein of SARS-CoV-2 increases intra-Golgi pH and forms a cation channel that is regulated by pH. J. Physiol. 599 (11), 2851–2868. doi:10.1113/JP281037
Case, D. A., Brozell, S. R., Cerutti, D. S., Cheatham, T. E., Cruzeiro, V. W. D., et al. (2018). Amber 2018. San Francisco: University of California.
Castaño-Rodriguez, C., Honrubia, J. M., Gutiérrez-Álvarez, J., DeDiego, M. L., Nieto-Torres, J. L., Jimenez-Guardeño, J. M., et al. (2018). Role of severe acute respiratory syndrome coronavirus viroporins E, 3a, and 8a in replication and pathogenesis. mBio 9 (3), e02325–e02317. doi:10.1128/mBio.02325-17
Chai, J., Cai, Y., Pang, C., Wang, L., McSweeney, S., Shanklin, J., et al. (2021). Structural basis for SARS-CoV-2 envelope protein recognition of human cell junction protein PALS1. Nat. Commun. 12 (1), 3433. doi:10.1038/s41467-021-23533-x
Chen, A. A., and Pappu, R. V. (2007). Parameters of monovalent ions in the AMBER-99 forcefield: Assessment of inaccuracies and proposed improvements. J. Phys. Chem. B 111 (41), 11884–11887. doi:10.1021/jp0765392
Collins, L. T., Elkholy, T., Mubin, S., Hill, D., Williams, R., Ezike, K., et al. (2021). Elucidation of SARS-cov-2 budding mechanisms through molecular dynamics simulations of M and E protein complexes. J. Phys. Chem. Lett. 12 (51), 12249–12255. doi:10.1021/acs.jpclett.1c02955
Darden, T., York, D., and Pedersen, L. (1993). Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems. J. Chem. Phys. 98 (12), 10089–10092. doi:10.1063/1.464397
Das, G., Das, T., Chowdhury, N., Chatterjee, D., Bagchi, A., and Ghosh, Z. (2021). Repurposed drugs and nutraceuticals targeting envelope protein: A possible therapeutic strategy against COVID-19. Genomics 113 (2), 1129–1140. doi:10.1016/j.ygeno.2020.11.009
Das, J. K., and Roy, S. (2021). A study on non-synonymous mutational patterns in structural proteins of SARS-CoV-2. Genome 64 (7), 665–678. doi:10.1139/gen-2020-0157
DeDiego, M. L., Alvarez, E., Almazán, F., Rejas, M. T., Lamirande, E., Roberts, A., et al. (2007). A severe acute respiratory syndrome coronavirus that lacks the E gene is attenuated in vitro and in vivo. J. Virol. 81 (4), 1701–1713. doi:10.1128/JVI.01467-06
DeDiego, M. L., Nieto-Torres, J. L., Jimenez-Guardeño, J. M., Regla-Nava, J. A., Castaño-Rodriguez, C., Fernandez-Delgado, R., et al. (2014). Coronavirus virulence genes with main focus on SARS-CoV envelope gene. Virus Res. 194, 124–137. doi:10.1016/j.virusres.2014.07.024
Dey, D., Borkotoky, S., and Banerjee, M. (2020). In silico identification of Tretinoin as a SARS-CoV-2 envelope (E) protein ion channel inhibitor. Comput. Biol. Med. 127, 104063. doi:10.1016/j.compbiomed.2020.104063
Dickson, C. J., Madej, B. D., Skjevik, Å. A., Betz, R. M., Teigen, K., Gould, I. R., et al. (2014). Lipid14: The amber lipid force field. J. Chem. Theory Comput. 10 (2), 865–879. doi:10.1021/ct4010307
Halgren, T. A., Murphy, R. B., Friesner, R. A., Beard, H. S., Frye, L. L., Pollard, W. T., et al. (2004). Glide: A new approach for rapid, accurate docking and scoring. 2. Enrichment factors in database screening. J. Med. Chem. 47 (7), 1750–1759. doi:10.1021/jm030644s
Heo, L., and Feig, M. (2020). Modeling of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) proteins by machine learning and physics-based refinement. bioRxiv. 2003, 2020.03.25.008904. doi:10.1101/2020.03.25.008904
Hsieh, Y.-C., Li, H.-C., Chen, S.-C., and Lo, S.-Y. (2008). Interactions between M protein and other structural proteins of severe, acute respiratory syndrome-associated coronavirus. J. Biomed. Sci. 15 (6), 707–717. doi:10.1007/s11373-008-9278-3
Humphrey, W., Dalke, A., and Schulten, K. (1996). Vmd: Visual molecular dynamics. J. Mol. Graph. 14 (1), 33–38. doi:10.1016/0263-7855(96)00018-5
Jackson, C. B., Farzan, M., Chen, B., and Choe, H. (2022). Mechanisms of SARS-CoV-2 entry into cells. Nat. Rev. Mol. Cell Biol. 23 (1), 3–20. doi:10.1038/s41580-021-00418-x
Jo, S., Kim, T., Iyer, V. G., and Im, W. (2008). CHARMM-GUI: A web-based graphical user interface for CHARMM. J. Comput. Chem. 29 (11), 1859–1865. doi:10.1002/jcc.20945
Khelashvili, G., Kollmitzer, B., Heftberger, P., Pabst, G., and Harries, D. (2013). Calculating the bending modulus for multicomponent lipid membranes in different thermodynamic phases. J. Chem. Theory Comput. 9 (9), 3866–3871. doi:10.1021/ct400492e
Klesse, G., Rao, S., Sansom, M. S. P., and Tucker, S. J. (2019). Chap: A versatile tool for the structural and functional annotation of ion channel pores. J. Mol. Biol. 431 (17), 3353–3365. doi:10.1016/j.jmb.2019.06.003
Kuzmin, A., Orekhov, P., Astashkin, R., Gordeliy, V., and Gushchin, I. (2022). Structure and dynamics of the SARS-CoV-2 envelope protein monomer. Proteins 90 (5), 1102–1114. doi:10.1002/prot.26317
Lau, T.-L., Kim, C., Ginsberg, M. H., and Ulmer, T. S. (2009). The structure of the integrin alphaIIbbeta3 transmembrane complex explains integrin transmembrane signalling. EMBO J. 28 (9), 1351–1361. doi:10.1038/emboj.2009.63
Li, Y., Surya, W., Claudine, S., and Torres, J. (2014). Structure of a conserved golgi complex-targeting signal in coronavirus envelope proteins. J. Biol. Chem. 289 (18), 12535–12549. doi:10.1074/jbc.M114.560094
Lim, K. P., and Liu, D. X. (2001). 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. J. Biol. Chem. 276 (20), 17515–17523. doi:10.1074/jbc.M009731200
Lopez, L. A., Riffle, A. J., Pike, S. L., Gardner, D., and Hogue, B. G. (2008). Importance of conserved cysteine residues in the coronavirus envelope protein. J. Virol. 82 (6), 3000–3010. doi:10.1128/JVI.01914-07
Madhavi Sastry, G., Adzhigirey, M., Day, T., Annabhimoju, R., and Sherman, W. (2013). Protein and ligand preparation: Parameters, protocols, and influence on virtual screening enrichments. J. Comput. Aided. Mol. Des. 27 (3), 221–234. doi:10.1007/s10822-013-9644-8
Maier, J. A., Martinez, C., Kasavajhala, K., Wickstrom, L., Hauser, K. E., and Simmerling, C. (2015). ff14SB: Improving the accuracy of protein side chain and backbone parameters from ff99SB. J. Chem. Theory Comput. 11 (8), 3696–3713. doi:10.1021/acs.jctc.5b00255
Mandala, V. S., McKay, M. J., Shcherbakov, A. A., Dregni, A. J., Kolocouris, A., and Hong, M. (2020). Structure and drug binding of the SARS-CoV-2 envelope protein transmembrane domain in lipid bilayers. Nat. Struct. Mol. Biol. 27 (12), 1202–1208. doi:10.1038/s41594-020-00536-8
Monje-Galvan, V., and Voth, G. A. (2021). Molecular interactions of the M and E integral membrane proteins of SARS-CoV-2. bioRxiv. 2004, 2021.04.29.442018. doi:10.1101/2021.04.29.442018
Mukherjee, S., Harikishore, A., and Bhunia, A. (2021). Targeting C-terminal helical bundle of NCOVID19 envelope (E) protein. Int. J. Biol. Macromol. 175, 131–139. doi:10.1016/j.ijbiomac.2021.02.011
Münz, M., Hein, J., and Biggin, P. C. (2012). The role of flexibility and conformational selection in the binding promiscuity of PDZ domains. PLoS Comput. Biol. 8 (11), e1002749. doi:10.1371/journal.pcbi.1002749
Nieto-Torres, J. L., Verdiá-Báguena, C., Jimenez-Guardeño, J. M., Regla-Nava, J. A., Castaño-Rodriguez, C., Fernandez-Delgado, R., et al. (2015). Severe acute respiratory syndrome coronavirus E protein transports calcium ions and activates the NLRP3 inflammasome. Virology 485, 330–339. doi:10.1016/j.virol.2015.08.010
Ohio-Supercomputer-Center (1987. Ohio Supercomputer Center. Columbus OH: Ohio Supercomputer Center. (Available at: http://osc.edu/ark:/19495/f5s1ph73).
Park, S. H., Siddiqi, H., Castro, D. V., De Angelis, A. A., Oom, A. L., Stoneham, C. A., et al. (2021). Interactions of SARS-CoV-2 envelope protein with amilorides correlate with antiviral activity. PLoS Pathog. 17 (5), e1009519. doi:10.1371/journal.ppat.1009519
Pervushin, K., Tan, E., Parthasarathy, K., Lin, X., Jiang, F. L., Yu, D., et al. (2009). Structure and inhibition of the SARS coronavirus envelope protein ion channel. PLoS Pathog. 5 (7), e1000511. doi:10.1371/journal.ppat.1000511
Planès, R., Bert, J.-B., Tairi, S., BenMohamed, L., and Bahraoui, E. (2022). SARS-CoV-2 envelope (E) protein binds and activates TLR2 pathway: A novel molecular target for COVID-19 interventions. Viruses 14 (5), 999. doi:10.3390/v14050999
Plescia, C. B., David, E. A., Patra, D., Sengupta, R., Amiar, S., Su, Y., et al. (2021). SARS-CoV-2 viral budding and entry can be modeled using BSL-2 level virus-like particles. J. Biol. Chem. 296, 100103. doi:10.1074/jbc.RA120.016148
Price, D. J., and Brooks, C. L. (2004). A modified TIP3P water potential for simulation with Ewald summation. J. Chem. Phys. 121 (20), 10096–10103. doi:10.1063/1.1808117
Roe, D. R., and Cheatham, T. E. (2013). PTRAJ and CPPTRAJ: Software for processing and analysis of molecular dynamics trajectory data. J. Chem. Theory Comput. 9 (7), 3084–3095. doi:10.1021/ct400341p
Ryckaert, J.-P., Ciccotti, G., and Berendsen, H. J. C. (1977). Numerical integration of the cartesian equations of motion of a system with constraints: Molecular dynamics of n-alkanes. J. Comput. Phys. 23 (3), 327–341. doi:10.1016/0021-9991(77)90098-5
Šali, A., Blundell, T. L., and Sali, A. (1993). Comparative protein modelling by satisfaction of spatial restraints. J. Mol. Biol. 234 (3), 779–815. doi:10.1006/jmbi.1993.1626
Sarkar, M., and Saha, S. (2020). Structural insight into the role of novel SARS-CoV-2 E protein: A potential target for vaccine development and other therapeutic strategies. PLOS ONE 15 (8), e0237300. doi:10.1371/journal.pone.0237300
Schoeman, D., and Fielding, B. C. (2019). Coronavirus envelope protein: Current knowledge. Virol. J. 16 (1), 69. doi:10.1186/s12985-019-1182-0
Schrödinger, L. L. C. (2021a). Schrödinger release. New York, NY: LigPrep, Schrödinger, LLC. 2022-2.
Schrödinger, L. L. C. (2021b). Schrödinger release. New York, NY: Maestro, Schrödinger, LLC. 2022-2.
Schrodinger, L. L. C. (2015). The AxPyMOL molecular graphics plugin for microsoft PowerPoint. Version 1.8
Sharma, T., Abohashrh, M., Baig, M. H., Dong, J.-J., Alam, M. M., Ahmad, I., et al. (2021). Screening of drug databank against WT and mutant main protease of SARS-CoV-2: Towards finding potential compound for repurposing against COVID-19. Saudi J. Biol. Sci. 28 (5), 3152–3159. doi:10.1016/j.sjbs.2021.02.059
Shepley-McTaggart, A., Sagum, C. A., Oliva, I., Rybakovsky, E., DiGuilio, K., Liang, J., et al. (2021). SARS-CoV-2 Envelope (E) protein interacts with PDZ-domain-2 of host tight junction protein ZO1. PLOS ONE 16 (6), e0251955. doi:10.1371/journal.pone.0251955
Singh Tomar, P. P., and Arkin, I. T. (2020). SARS-CoV-2 E protein is a potential ion channel that can be inhibited by Gliclazide and Memantine. Biochem. Biophys. Res. Commun. 530 (1), 10–14. doi:10.1016/j.bbrc.2020.05.206
Siu, Y. L., Teoh, K. T., Lo, J., Chan, C. M., Kien, F., Escriou, N., et al. (2008). 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. 82(22), 11318–11330. doi:10.1128/JVI.01052-08
Suk, J.-E., Situ, A. J., and Ulmer, T. S. (2012). Construction of covalent membrane protein complexes and high-throughput selection of membrane mimics. J. Am. Chem. Soc. 134 (22), 9030–9033. doi:10.1021/ja304247f
Surya, W., Li, Y., Millet, O., Diercks, T., and Torres, J. (2013). Transmembrane and juxtamembrane structure of αL integrin in bicelles. PLOS ONE 8 (9), e74281. doi:10.1371/journal.pone.0074281
Surya, W., Li, Y., and Torres, J. (2018). Structural model of the SARS coronavirus E channel in LMPG micelles. Biochim. Biophys. Acta. Biomembr. 1860 (6), 1309–1317. doi:10.1016/j.bbamem.2018.02.017
Teoh, K.-T., Siu, Y.-L., Chan, W.-L., Schlüter, M. A., Liu, C.-J., Peiris, J. S. M., et al. (2010). The SARS coronavirus E protein interacts with PALS1 and alters tight junction formation and epithelial morphogenesis. Mol. Biol. Cell 21 (22), 3838–3852. doi:10.1091/mbc.e10-04-0338
Tilocca, B., Soggiu, A., Sanguinetti, M., Babini, G., De Maio, F., Britti, D., et al. (2020). Immunoinformatic analysis of the SARS-CoV-2 envelope protein as a strategy to assess cross-protection against COVID-19. Microbes Infect. 22 (4), 182–187. doi:10.1016/j.micinf.2020.05.013
Troyano-Hernáez, P., Reinosa, R., and Holguín, Á. (2021). 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. Viruses 13 (2), 243. doi:10.3390/v13020243
Verdiá-Báguena, C., Aguilella, V. M., Queralt-Martín, M., and Alcaraz, A. (2021). Transport mechanisms of SARS-CoV-E viroporin in calcium solutions: Lipid-dependent Anomalous Mole Fraction Effect and regulation of pore conductance. Biochim. Biophys. Acta. Biomembr. 1863 (6), 183590. doi:10.1016/j.bbamem.2021.183590
Verdiá-Báguena, C., Nieto-Torres, J. L., Alcaraz, A., DeDiego, M. L., Torres, J., Aguilella, V. M., et al. (2012). Coronavirus E protein forms ion channels with functionally and structurally-involved membrane lipids. Virology 432 (2), 485–494. doi:10.1016/j.virol.2012.07.005
Watson, M. C., Brandt, E. G., Welch, P. M., and Brown, F. L. H. (2012). Determining biomembrane bending rigidities from simulations of modest size. Phys. Rev. Lett. 109 (2), 028102. doi:10.1103/PhysRevLett.109.028102
Weiss, S. R., and Navas-Martin, S. (2005). Coronavirus pathogenesis and the emerging pathogen severe acute respiratory syndrome coronavirus. Microbiol. Mol. Biol. Rev. 69(4), 635–664. doi:10.1128/MMBR.69.4.635-664.2005
Williams, C. J., Headd, J. J., Moriarty, N. W., Prisant, M. G., Videau, L. L., Deis, L. N., et al. (2018). MolProbity: More and better reference data for improved all-atom structure validation. Protein Sci. 27 (1), 293–315. doi:10.1002/pro.3330
Wishart, D. S., Feunang, Y. D., Guo, A. C., Lo, E. J., Marcu, A., Grant, J. R., et al. (2017). DrugBank 5.0: A major update to the DrugBank database for 2018. Nucleic Acids Res. 46 (D1), D1074–D1082. doi:10.1093/nar/gkx1037
World Health Organization (2022). [Online]Availble At: https://covid19.who.int/(Accessed Janurary, 2022).
Wu, E. L., Cheng, X., Jo, S., Rui, H., Song, K. C., Dávila-Contreras, E. M., et al. (2014). CHARMM-GUI Membrane Builder toward realistic biological membrane simulations. J. Comput. Chem. 35 (27), 1997–2004. doi:10.1002/jcc.23702
Xia, B., Shen, X., He, Y., Pan, X., Liu, F.-L., Wang, Y., et al. (2021). SARS-CoV-2 envelope protein causes acute respiratory distress syndrome (ARDS)-like pathological damages and constitutes an antiviral target. Cell Res. 31, 847–860. doi:10.1038/s41422-021-00519-4
Keywords: COVID-19, coronavirus, envelope (E) protein, SARS-CoV-2, E protein inhibitors
Citation: Yang R, Wu S, Wang S, Rubino G, Nickels JD and Cheng X (2022) Refinement of SARS-CoV-2 envelope protein structure in a native-like environment by molecular dynamics simulations. Front. Mol. Biosci. 9:1027223. doi: 10.3389/fmolb.2022.1027223
Received: 24 August 2022; Accepted: 21 September 2022;
Published: 10 October 2022.
Edited by:
Jian Wang, The Pennsylvania State University, United StatesReviewed by:
Yi He, University of New Mexico, United StatesYunjie Zhao, Central China Normal University, China
Copyright © 2022 Yang, Wu, Wang, Rubino, Nickels and Cheng. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Sijin Wu, wusj@dicp.ac.cn; Jonathan D. Nickels, Jonathan.Nickels@uc.edu; Xiaolin Cheng, cheng.1302@osu.edu