key: cord-0741072-rnaktik9
authors: Surya, Wahyu; Li, Yan; Torres, Jaume
title: Structural model of the SARS coronavirus E channel in LMPG micelles
date: 2018-02-21
journal: Biochim Biophys Acta Biomembr
DOI: 10.1016/j.bbamem.2018.02.017
sha: 2c01b65a4468d43dad054457c80d05d5a2715ac6
doc_id: 741072
cord_uid: rnaktik9

Coronaviruses (CoV) cause common colds in humans, but are also responsible for the recent Severe Acute, and Middle East, respiratory syndromes (SARS and MERS, respectively). A promising approach for prevention are live attenuated vaccines (LAVs), some of which target the envelope (E) protein, which is a small membrane protein that forms ion channels. Unfortunately, detailed structural information is still limited for SARS-CoV E, and non-existent for other CoV E proteins. Herein, we report a structural model of a SARS-CoV E construct in LMPG micelles with, for the first time, unequivocal intermolecular NOEs. The model corresponding to the detergent-embedded region is consistent with previously obtained orientational restraints obtained in lipid bilayers and in vivo escape mutants. The C-terminal domain is mostly α-helical, and extramembrane intermolecular NOEs suggest interactions that may affect the TM channel conformation.

Coronaviruses (CoV) typically affect the respiratory tract and gut of mammals and birds. Approximately 30% of common colds are caused by two human coronaviruses -OC43 and 229E. Of particular interest are the viruses responsible for the severe acute respiratory syndrome (SARS), which produced a near pandemic in 2003 [1] , and the recent Middle East respiratory syndrome coronavirus (MERS-CoV) [2] .

No effective licensed treatments exist against coronavirus infections [3] [4] [5] , but live attenuated vaccines (LAVs) [6] [7] [8] [9] [10] and fusion inhibitors [11] are promising strategies. One CoV component critical for pathogenesis is the envelope (E) protein, as reported in several coronaviruses, e.g., MERS and SARS-CoVs [12] [13] [14] . The CoV envelope (E) proteins are short polypeptides (76-109 amino acids) with a single α-helical transmembrane (TM) domain [15] [16] [17] [18] [19] [20] [21] that form homopentameric ion channels (IC) with poor ion selectivity [22, 23] . CoV E proteins are mostly found in the endoplasmic reticulum-Golgi intermediate compartment (ERGIC) [24] [25] [26] [27] [28] [29] . In animal models, deletion of SARS-CoV E protein reduced pathogenicity and mortality [30] , whereas cellular models displayed up-and down-regulation of stress response and inflammation host genes, respectively [31] . The importance of E protein in pathogenesis has led to the development of LAVs based on deletion of E protein in SARS-and MERS-CoVs, although this led to compensatory mechanisms that recover virulence [32, 33] .

Specific critical features in the SARS-CoV E protein sequence have been identified that determine virulence, e.g., at the C-terminal tail [34] or in the TM domain [30] , and precise structural characterization of these regions could help in the design of E protein-based CoV LAVs. However, detailed structural knowledge is still very limited in the case of SARS-CoV E, and non-existent for other CoV E proteins.

A pentameric model for SARS-CoV E was initially proposed by the authors after an in silico conformational search [15] of TM domain oligomers. In that report, two pentameric models (termed 'A' and 'B') that were separated by a~50°rotation of their α-helices were selected. In model A, V25 adopts a more lumenal position, whereas in model B, the position of this residue is clearly interhelical (Fig. 1) . The pentameric organization of SARS-CoV E has been confirmed experimentally in various detergents: PFO, DPC or C-14 betaine [17, 18] , not only for synthetic TM (E TM ), but also for an 8-65 (E TR ) construct and for full length E protein (E FL ).

To confirm experimentally the orientation of the α-helices in the pentameric model, site specific infrared dichroism (SSID) measurements [35] were obtained in hydrated lipid bilayers, with 13 C = 18 O isotopically labeled synthetic E TM . However, the orientation of the α-helices turned out to be strongly dependent on the presence of 2 flanking lysine residues at each end of the peptides [16] : with flanking lysine residues, the orientation was a hybrid between models A and B (residues 17-24 were oriented consistent with model B, but from residue 24 onwards, orientation was as expected for model A), consistent with a 'bend of the α-helices around residues 25-27' [16] . Without terminal lysines, however, the orientation of the central five labeled consecutive residues, L21 to V25, was entirely consistent with model A [17] .

These initial results suggested that the conformation of the E TM pentamer may be very sensitive of the presence of extra residues and probably also, extramembrane domains. An NMR study was performed on a synthetic E TM (residues in DPC detergent micelles, where E TM was selectively labeled [20] . E TM was 15 N-labeled at A22, V24, V25, and 13 C, 15 Orientation of computational models A (orange) and B (cyan) [15] , where the side chain of V25 (F26 is only used to guide the eye) is indicated. The 'A-like' model obtained by NMR [20] is shown in red. In model B, the position of V25 is clearly interhelical. intermolecular NOEs were insufficient to distinguish between models A and B, and the monomer structure was fit to a model A template. More recently, recombinant SARS-CoV escape mutants were recovered after introducing a V25F channel-inactivating mutation in the E protein, [36] , that led to attenuation in a mouse model [30] . Revertant mutants regained fitness and pathogenicity whereas mutated E protein regained channel activity [30] . Surprisingly, escape mutations in E protein clustered along the helix interface opposite to residue V25, consistent with an interhelical orientation of this residue, as found in model B (Fig. 1, cyan) .

In the present paper, we report a more accurate model of the SARS-CoV E protein pentamer, in LMPG micelles. The construct we have used prolongs the TM domain with another 27 residues in the C-terminal domain (residues 8-65). Following established protocols [37] , two types of monomers were mixed, bearing different isotopical labels, that allowed unambiguous identification of ten intermonomeric NOEs. In a nutshell, the results are consistent with a TM model that appears to be a hybrid between models A and B: while overall being closer to model A, residue V25 has a clear 'model B-like' interhelical orientation, consistent with the revertant mutants that appeared in vivo.

The expression and purification methods for the truncated SARS-CoV E construct corresponding to residues 8-65 (E TR ) have been described previously [19] . This construct does not have cysteines, as these are not required for oligomerization [18, 19, 28, 38] . In the present work, M9 media was supplemented with an appropriate combination of 15 

Blue-native PAGE (BN-PAGE) was performed as described previously [41] . Lyophilized E TR protein was solubilized (0.1 mM) in sample buffer containing LMPG (lyso-myristoyl phosphatidylglycerol, Anatrace) at the indicated concentrations.

For α-helical bundle models, the rotational pitch angle of a residue, ω, defined arbitrarily as 0°or 180°when transition dipole moment, helix director, and the z-axis all reside in a single plane, was calculated as described elsewhere [42] . The final result is the average of the ω values calculated in each monomer. For a canonical α-helix, it is expected that Δω between two consecutive residues is~100°.

Lyophilized E TR protein (0.67 mM) was solubilized in 20 mM sodium phosphate pH 5.5, 50 mM NaCl, and 200 mM LMPG, i.e., a protein:detergent (P/D) molar ratio of 1:300. The same protein concentration and P/D ratio was used for the mixture of 15 N-D and 13 Clabeled samples. The solution was vortexed and sonicated several times until a clear solution was obtained, indicating protein reconstitution into detergent micelles.

NMR experiments were performed at 308 K using an Avance-II 700 NMR spectrometer with cryogenic probes. Sodium 2,2-dimethyl-2-silapentane-5-sulfonate (DSS) was used as the internal reference for 1 H nuclei. The chemical shifts of 13 C and 15 N nuclei were calculated from the 1 H chemical shifts. The NMR data were processed using TopSpin 3.1 (www.bruker-biospin.com) and analyzed using CARA (www.nmr.ch). Sequence-specific assignment of backbone 1 H N , 15 N, 13 C′ and 13 C α was achieved by using 2D [ 1 H- 15 

Intra-monomeric NOE distance restraints were obtained from 15 N-NOESY-HSQC and 13 C-NOESY-HSQC spectra (both with a mixing time of 120 ms). Backbone dihedral angle restraints (φ and ψ) were derived from 13 C′, 13 C α , 13 C β , 1 H α and 1 H β chemical shift values using TALOS+ [43] . Short-range and medium range NOE connectivities were used to establish sequence-specific 1 H NMR assignments and to identify elements of the regular secondary structure. Hydrogen bonds were derived from the NOE connectivity, and supported by the H/D exchange data. Monomer structure calculations were performed using CYANA 3.0 [44, 45] and visualized using PyMOL (Delano Scientific). All of the restraints used in the calculations to obtain a total of 10 monomer structures, and all the structure statistics, are summarized in Supplementary Tables S1 and S2.

Inter-monomeric NOE restraints were obtained from 3D 15 N-resolved NOESY-HSQC (250 ms mixing time) of two sets of asymmetrically deuterated samples: (1) 15 N/ 2 H-labeled E TR sample (ND), and (2) an equimolar mixture of 15 N/ 2 H-labeled and a non-deuterated 13 Clabeled E TR sample (ND + C). NOE cross-peaks appearing in sample ND + C but not in sample ND were assigned to inter-monomeric contacts. Conversely, resonances also appearing in the ND sample were attributed to incomplete deuteration, and were assigned to intramonomeric NOEs.

The pentamer structure was calculated using HADDOCK 2.2 [46] according to standard protocols. Ten inter-monomeric NOE restraints (defined as above) were described as ambiguous and unambiguous 5.0 Å distance restraints. Two segments were described as fully flexible: residues 37-47 and 40-54. A C5 symmetry restraint between all 5 subunits and pairwise non-crystallographic symmetry restraints between neighbouring subunits were applied. Initial rigid-body docking yielded 1000 structures, out of which 200 top-scoring structures (i.e., based on HADDOCK target function score) were selected for refinement by semi-flexible simulated annealing. These were then clustered based on RMSD, and the top-scoring cluster was selected (all 16 structures within the said cluster were grouped to form an ensemble).

Despite phospholipid isotropic bicelles may have been more membrane-like than detergent micelles, in our hands, phospholipid bicelles did not produce suitable spectra of E TR (not shown). Examples of significant differences observed in bicelles vs micelles have been reported, e.g., in the study of the integrin TM heterodimers [47] [48] [49] [50] [51] [52] or in viral channels [53] .

Nevertheless, we have shown previously that E TR is pentameric in various detergents [17, 18] , although none of them was suitable for NMR studies of E TR or E FL (not shown). E TR only produced reasonably good NMR spectra in DPC when SDS was also present [19] , but since SDS disrupts E TR oligomerization, we searched for other micellar environments. Lipid-like LMPG was found to produce good NMR spectra for E TR , although not for E FL . Therefore, E TR in LMPG was used in subsequent experiments. The use of the E TR construct instead of the fulllength E protein (E FL ) is justified since the 13 Cα chemical shifts of E TR and E FL protein in SDS or SDS/DPC were almost identical for residues 8-65 [19] . In addition, the secondary structure, obtained by CD/FTIR [18] , of E TR and E FL is similar and predominantly α-helical, whether in DPC, SDS, mixed (1:2 M ratio) SDS/DPC micelles or DMPC synthetic membranes [18, 19] .

Comparison of the HSQC spectrum of E TR /LMPG before and after exposure to D 2 O ( Fig. 2A) shows that only 20 residues are protected from hydrogen/deuterium (H/D) exchange. The protected residues correspond to the stretch L18-L37, unequivocally indicating the presence of a single TM domain in SARS-CoV E. This result is consistent with the stretch L18-L39 found to be protected in SDS micelles [19] . The chemical shift index (CSI)-based secondary structure of E TR (calculated by using TALOS+) obtained in LMPG (Fig. 2B) , has significantly higher helicity in C-terminal residues 52-55, when compared with the data obtained SDS or with a mixture SDS/DPC [19] .

The structure of E TR was calculated from 10 E TR monomer structures (Fig. 3A ) and the structure statistics are summarized in Supplementary  Table S1 . The E TR monomer in LMPG consists of three helical segments: the one encompassing the TM domain (H1, residues 12-37), a juxtamembrane middle helical segment (H2, residues [39] [40] [41] [42] [43] [44] [45] [46] [47] , and a Cterminal helix (H3, residues 52-65) (Fig. 3B) . In contrast, E TR in DPC/ SDS [19] was formed by only two helical segments separated by a long flexible link (Fig. 3C) . Compared to the results in SDS or SDS/DPC [19] , in LMPG helix H3 is extended by 3 residues on its N-terminal side, whereas a new helical segment, H2, is formed.

To assess the oligomerization of E TR in LMPG micelles, its migration in a BN-PAGE gel was analyzed at various protein-to-detergent (P/D) ratios (Fig. 4) . At the lowest P/D molar ratio (1:1000), E TR migrates as a ladder of increasingly larger oligomers where the fastest migrating band is assumed to correspond to monomers (lower star),~8 kDa, whereas at a high P/D ratio (1:125), E TR migrates with an apparent molecular weight of~150 kDa. These results are almost identical to those obtained previously for MERS-CoV E, for which a pentameric oligomer was determined using analytical ultracentrifugation in C-14 betaine. In that case, migration in BN-PAGE gels was also observed as a single~150 kDa band in detergents DPC, DHPC and LMPG [21] , and the ladder observed at higher detergent concentration conveniently provided an internal reference that served as a oligomeric size marker. Similar to E TR , by comparison with that ladder, we confidently assigned the single band observed for MERS-CoV E to pentameric oligomers. It should be noted that in BN-PAGE gels of membrane proteins, molecular weights can appear up to 80% higher due to a contribution of the dye [54] . We have shown this for tetrameric AQPZ, which migrated at 170 kDa instead of the expected~100 kDa, and with a viroporin, the SH protein pentamer [41] , which migrated as~66 kDa instead of 40 kDa. In the case of envelope E proteins, the effect is even more pronounced. In both SARS-CoV E TR and MERS-CoV E, the pentameric form appears at~150 kDa, therefore the monomer should appear at > 30 kDa. This is consistent with its migration above the AQPZ monomer (~25 kDa). The ladder ends with a pentamer, which is the predominant band at high P/D ratios. The proportion of large oligomers naturally decrease at low P/D ratios, but a significant amount of pentamer species is still present even at the 1:1000 P/D ratio. The NMR data was collected at a P/D molar ratio of 1:300, which should mostly be formed by pentamers. 

In a previous paper, we showed that monomeric E TR in SDS micelles was not affected by addition of the drug HMA [19] . However, after addition of DPC to SDS, to a SDS/DPC 1:4 M ratio, HMA induced clear chemical shift perturbations (CSPs), concomitant with E TR oligomerization. The oligomerization in DPC/SDS was not homogeneous, which precluded a more detailed study, whereas in LMPG a predominant oligomeric size is observed at a high protein-detergent ratio (Fig. 4) . Therefore, in LMPG the changes observed after HMA addition should more reliably represent the binding of HMA to E TR . HMA-induced CSPs were detected herein after addition of 7.75 mM HMA to 0.25 mM E TR in 200 mM LMPG micelles (P/D molar ratio 1:800) (Fig. 5) .

The average CSP value was 0.019 ppm, and several residues showed CSP > 1 S.D. from the average value, notably Thr-9, Leu-12, Ile-13, Ala-36 and Val-47. These results suggest the presence of two binding sites located at both ends of the TM domain. Given the long distance between Ala-36 and Val-47, the two HMA-interacting residues may be located in different monomers.

A pentameric model was obtained by docking the monomeric form of E TR using HADDOCK 2.2 [46] , which incorporated 10 inter-monomeric NOE restraints (Fig. 6A) . We note that 2 inter-monomeric NOEs are located at the extramembrane C-terminal tail: L39 HN -Y57 HB and V47 HN -N64 HN. The same figure shows a representative example of NOE E TR inter-monomer connectivity (Fig. 6B) . The remaining plots of inter-monomeric NOEs are shown in Fig. S1 . Structure statistics are summarized in Supplementary Table S2 .

The E TR pentamer is a right handed α-helical bundle where the Cterminal tails coil around each other (Fig. 7A ) likely owing to the 2 inter-monomeric restraints between the two C-terminal helices. Each pentamer subunit (Fig. 7B) has better defined structure compared to the monomer alone (Fig. 3A) . This is mainly due to decreased flexibility at the inter-helical segments, which were kept flexible during the docking, as the two C-terminal helices now adopt a relatively fixed conformation. This is also apparent from the RMSD values; the pentamer subunit RMSD values are significantly reduced as compared to the monomer (Fig. 7C) .

Notably, in this pentameric model, the location of V25 is interhelical ( Fig. 8B-C) , whereas in the previously proposed model it was closer to a lumenal orientation [16] . The rotational pitch of the residues in the TM domain of this pentameric model were measured individually [35] and compared to those from the computational models A and B [15] (Fig. 8D) . While values for residues 25-27 are closer to model B, the rest of the sequence is similar to model A, except at residue 28 which deviates from both models. For comparison, the rotational pitch close to model A for residues in E TM obtained previously by NMR in DPC micelles [20] is also shown. The present model has been constructed independently from A and B model templates, and the result appears to be a hybrid between the two [15] . This is not surprising since the in silico study assumed a certain rigidity in the TM α-helices [15] . Most of the residues in the model we report have an orientation consistent with model A. This is not surprising, since model A had the lowest energy value for each individual E protein homologs [15] . However, the model gets closer to model B in the turn that contains V25 (Fig. 8D ). This enables V25 to adopt a more interhelical orientation consistent with the revertant mutants that appeared in vivo [30] . Additionally, the helix kink region suggested by infrared dichroism data in lipid bilayers [16] is also observed, which supports the validity of the membrane-mimic environment used herein.

Finally, in LMPG micelles, the C-terminal tail of SARS-CoV E protein is α-helical, more so than observed in mixed DPC/SDS micelles [19] , and the presence of extramembrane NOEs suggest interactions between the C-terminal domains that may affect the pentameric conformation.

The atomic coordinates have been deposited in the Protein Data Fig. 6 . Inter-monomeric NOEs in E TR pentamer. (A) List of intermonomeric NOE contacts, with those located in the extramembrane C-terminal region in bold; (B) a representative example of NOE E TR inter-monomer connectivity (green lines). Selected strips correspond to a 15 N-NOESY-HSQC spectrum and NH protons of V14 for samples 15 N/ 2 H-labeled (ND), 15 N/ 2 Hlabeled + 13 C-labeled (ND + C), and 15 N/ 13 C-labeled (NC). The NOE strips from the NC sample are shown as reference, as they contain both intra and inter-monomer contacts. Strips corresponding to the remaining NOE connectivity are shown in Fig.  S1 . (Layout note: 1 column). Fig. 3A , black), structured helical segments of the monomer (see Fig. 3B , blue), and the pentamer ensemble (Fig. 7A, red) . The average RMSD value of the monomer (dashed line) and ± 1 S.D. values (grey band) are indicated. (Layout note: 1.5 column). The N-and C-terminus of one monomeric unit is indicated; (C) top view of a monomer-monomer TM interaction, showing the distances between the side chain of V25 and those of residues appearing in SARS-CoV E V25F revertant mutants [30] ; (D) differences in TM residue rotational orientation, ω, between the experimental model proposed here (LMPG) versus computational models A and B [15] and that of E TM obtained by NMR in DPC micelles [20] . The region with larger differences between the present model (LMPG) and model A (residues 25-28) is highlighted. (Layout note: 1.5 columns). Bank (PDB ID: 5X29). Assigned chemical shifts have been deposited at the Biological Magnetic Resonance Bank (BMRB ID: 36049).

The http://dx.doi.org/10.1016/j.bbamem.2018.02.017 associated with this article can be found, in online version.

SARS coronavirus: a new challenge for prevention and therapy

MERS: emergence of a novel human coronavirus

Current progress in antiviral strategies

Cell-based antiviral screening against coronaviruses: developing virus-specific and broad-spectrum inhibitors

Assessing activity and inhibition of Middle East respiratory syndrome coronavirus papain-like and 3C-like proteases using luciferase-based biosensors

Development and characterization of a recombinant infectious bronchitis virus expressing the ectodomain region of S1 gene of H120 strain

Production and immunogenicity of chimeric virus-like particles containing the spike glycoprotein of infectious bronchitis virus

A live, impaired-fidelity coronavirus vaccine protects in an aged, immunocompromised mouse model of lethal disease

Recombinant live vaccines to protect against the severe acute respiratory syndrome coronavirus

Review of infectious bronchitis virus around the world

Ready, set, fuse! The coronavirus spike protein and acquisition of fusion competence

A live attenuated severe acute respiratory syndrome coronavirus is immunogenic and efficacious in golden Syrian hamsters

Immunization with an attenuated severe acute respiratory syndrome coronavirus deleted in E protein protects against lethal respiratory disease

Engineering a replication-competent, propagation-defective Middle East respiratory syndrome coronavirus as a vaccine candidate

The transmembrane oligomers of coronavirus protein E

Model of a putative pore: the pentameric α-helical bundle of SARS coronavirus E protein in lipid bilayers

Structural flexibility of the pentameric SARS coronavirus envelope protein ion channel

Expression and purification of coronavirus envelope proteins using a modified betabarrel construct

Structure of a conserved Golgi complextargeting signal in coronavirus envelope proteins

Structure and inhibition of the SARS coronavirus envelope protein ion channel

MERS coronavirus envelope protein has a single transmembrane domain that forms pentameric ion channels

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

Hexamethylene amiloride blocks E protein ion channels and inhibits coronavirus replication

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

The 9-kDa hydrophobic protein encoded at the 3′ end of the porcine transmissible gastroenteritis coronavirus genome is membrane-associated

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

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

Importance of conserved cysteine residues in the coronavirus envelope protein

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

Inhibition of NF-κB-mediated inflammation in severe acute respiratory syndrome coronavirus-infected mice increases survival

Severe acute respiratory syndrome coronavirus envelope protein regulates cell stress response and apoptosis

Identification of the mechanisms causing reversion to virulence in an attenuated SARS-CoV for the design of a genetically stable vaccine

The PDZ-binding motif of severe acute respiratory syndrome coronavirus envelope protein is a determinant of viral pathogenesis

Severe acute respiratory syndrome coronaviruses with mutations in the E protein are attenuated and promising vaccine candidates

Multiple site-specific infrared dichroism of CD3-zeta, a transmembrane helix bundle

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

Unusual architecture of the p7 channel from hepatitis C virus

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

Practical protocols for production of very high yields of recombinant proteins using Escherichia coli

Isotope labeling strategies for the study of highmolecular-weight proteins by solution NMR spectroscopy

A cost-effective method for simultaneous homo-oligomeric size determination and monodispersity conditions for membrane proteins

Site-directed dichroism as a method for obtaining rotational and orientational constraints for oriented polymers

Protein backbone angle restraints from searching a database for chemical shift and sequence homology

Torsion angle dynamics for NMR structure calculation with the new program DYANA

Protein NMR structure determination with automated NOE assignment using the new software CANDID and the torsion angle dynamics algorithm DYANA

HADDOCK: a protein−protein docking approach based on biochemical or biophysical information

The structure of the integrin alphaIIbbeta3 transmembrane complex explains integrin transmembrane signalling

Transmembrane and juxtamembrane structure of alphaL integrin in bicelles

Structure of the integrin alphaIIb transmembrane segment

Oligomerization of the integrin αIIbβ3: roles of the transmembrane and cytoplasmic domains

Construction of covalent membrane protein complexes and high-throughput selection of membrane mimics

Integrin α1 has a long helix, extending from the transmembrane region to the cytoplasmic tail in detergent micelles

The small hydrophobic protein of the human respiratory syncytial virus forms pentameric ion channels

Oligomeric state of membrane transport proteins analyzed with blue native electrophoresis and analytical ultracentrifugation

Acknowledgements J.T. acknowledges the funding from Singapore MOE Tier 1 grant RG 51/13. The authors declare no conflict of interest.

Supplementary data to this article can be found online at https:// doi.org/10.1016/j.bbamem.2018.02.017.