Previous Article | Next Article ![]()
Journal of Virology, February 2002, p. 1944-1958, Vol. 76, No. 4
0022-538X/01/$04.00+0 DOI: 10.1128'JVI.76.4.1944-1958.2002
Copyright © 2002, American Society for Microbiology. All Rights Reserved.
Centre de Biophysique Moléculaire Numérique, Faculté Universitaire des Sciences Agronomiques de Gembloux, B-5030 Gembloux,1 Innogenetics, B-9052 Ghent, Belgium2
Received 18 June 2001/ Accepted 18 October 2001
|
|
|---|
-helix followed by a pair of transmembrane ß-strands. Several three-dimensional (3-D) models were generated. After energy minimization, their ability to interact with membranes was studied using the molecular hydrophobicity potentials calculation and the IMPALA procedure. The latter simulates interactions with a membrane by a Monte Carlo minimization of energy. These methods suggest that the ß-hairpins could anchor the glycoproteins in the ER membrane at least transiently. Anchoring could be stabilized by the adsorption of the nearby amphipathic
-helices at the membrane surface. The 3-D models correlate with experimental results which indicate that the E1-E2 transmembrane domains are involved in the heterodimerization and have ER retention properties. |
|
|---|
The E1 and E2 proteins have hydrophobic carboxy-terminal (C-terminal) regions that presumably act as membrane anchors (51). Indeed, deletion of this region of E2 leads to its secretion (41, 53). The E2 membrane anchor should begin at amino acid 718 of the polyprotein (42). The situation is not as clear for E1, since a truncated form ending at amino acid 340 is secreted only if it contains an internal deletion between amino acids 262 and 290, suggesting the presence of a second membrane anchor (39). However, for Hussy et al. (26) and Michalak et al. (41), truncated forms ending at amino acid 311 and 334 with no additional deletion can be secreted. E1, like E2, is therefore probably anchored only by its C-terminal hydrophobic domain. The limits of this transmembrane domain have not yet been clearly determined.
Some data suggest that the E1 and E2 truncated proteins are secreted as complexes in the bathing medium of cell cultures (39), while others indicate that the C-terminal part of E2 is required to form the E1-E2 complexes (53). These contradictions were reanalyzed by Michalak et al. (41). Using immunoprecipitation with anti-E1 or anti-E2 monoclonal antibodies, they detected no E1-E2 dimer when deletions were introduced into the C terminus of E2. They therefore confirmed the observations by Selby et al. (53) and pointed out that the C-terminal domain of E2 should be important for heterodimerization (41). More recently, the crucial role of E1 and E2 C-terminal domains in heterodimerization was demonstrated through mutagenesis studies (11, 45).
Since no efficient cell culture replication system is available, the assembly and release of HCV particles have not been directly studied. Some observations such as the lack of modification of the glycans by the Golgi enzymes, the location of the HCV glycoproteins in the endoplasmic reticulum (ER), and their absence at the cell surface suggest that the mature E1-E2 heterodimer does not leave the ER and that E1 and/or E2 contains an ER retention signal (15, 16, 48, 55). The ER retention is static rather than due to a recycling from the cis-Golgi region (19). In 1998, Cocquerel et al. (10) showed that the E2 glycoprotein is retained in the ER as the E1-E2 heterodimer because of a specific signal in the C-terminal 29 amino acids. Replacement of this domain with the anchor domain of a protein normally targeted to the plasma membrane (the CD4 protein) is sufficient to drive E2 at the cell surface. Moreover, when the ectodomain of CD4 is fused to the transmembrane domain of E2, the chimeric protein is stored in the ER (10). The same authors (9) looked for an ER retention signal in E1 by making chimeric proteins with the ectodomain of CD4 or CD8 and the C-terminal hydrophobic sequence of E1 (C-terminal 31 amino acids). They showed that the transmembrane domain of E1 retains the chimera in the ER. In addition, the glycans of the chimeric proteins were not processed by the Golgi enzymes. This indicates that the transmembrane domain of E1 is responsible for a true retention in the ER. Charged residues of the E1 and E2 C-terminal domains have been shown to be crucial for the ER retention (11).
We studied the E1 and E2 C-terminal anchor domains in silico using molecular modeling. Analysis of 25 HCV strains, as well as calculation and characterization of three-dimensional (3-D) models, suggests a peculiar conformation of the E1 and E2 anchor regions. Both C-terminal domains should be one amphipathic
-helix followed by a ß-hairpin. The ß-strands could insert in the membrane and stay there, at least transiently. The expected properties of the 3-D models correlate with the heterodimerization and ER retention properties of E1 and E2.
|
|
|---|
HCV sequences. The E1 and E2 sequences of the following 25 HCV strains were extracted from PIR, SwissProt, EMBL, and GenBank sequence databases: HCV-BEINNX11 (GenBank accession no. A48711), HCV-BK (SwissProt accession no. P26663), HCV-J (SwissProt accession no. P26662), HCV-J33 (GenBank accession no. D14484), HCV-J4/83 (GenBank accession no. D13558), HCV-J4/91 (GenBank accession no. D10750), HCV-JK1 (PIR accession no. S18030), HCV-JT (GenBank accession no. D11168), HCV-JT" (GenBank accession no. D11355), HCV-L2 (GenBank accession no. U01214), HCV-N (GenBank accession no. S62220), HCV-T (SwissProt accession no. P29846), HCV-Unkcds (EMBL accession no. M96362), HCV-K1-R1 (GenBank accession no. D50480), HCV-K1-R2 (GenBank accession no. D50481), HCV-K1-R3 (GenBank accession no. D50482), HCV-K1-S1 (GenBank accession no. D50483), HCV-K1-S2 (GenBank accession no. D50485), HCV-K1-S3 (GenBank accession no. D50484), HCV-1 (SwissProt accession no. P26664), HCV-H (SwissProt accession no. P27958), HCV-J1 (EMBL accession no. D10749), HCV-G9 (GenBank accession no. D14853), HCV-J6 (SwissProt accession no. P26660), and HCV-J8 (SwissProt accession no. P26661). Note that amino acid numbers always refer to the mature HCV-BEINNX11 glycoproteins. E1 extends from amino acid 1 (= amino acid 192 in the polyprotein) to amino acid 192 (= amino acid 383 in the polyprotein), and E2 extends from amino acid 1 (= amino acid 384 in the polyprotein) to amino acid 363 (= amino acid 746 in the polyprotein).
Kyte and Doolittle's method. The mean hydrophobicity of a stretch of sequence is calculated by moving a seven-residue window along the sequence (30). The value is attributed to the central amino acid of the window, and the hydrophobicity of each amino acid is defined in the original work (30).
HCA.
The hydrophobic cluster analysis (HCA) plot (21) relies upon the 2-D helical representation of the sequence: the sequence is written on a classical
-helix (3.6 residues per turn) smoothed on a cylinder. To make the view easier to handle, the cylinder is then cut parallel to its axis and unrolled. As some adjacent amino acids are separated by the flattening of the cylinder, the plot is duplicated to restore all the connections of residues. The hydrophobic residues (F, V, L, I, M, Y, and W) are circled and hatched; proline is represented by a star, glycine is represented by a diamond, serine is represented by a square with a dot in it, and threonine is represented by an empty square. This method allows ß-strands to be distinguished from
-helices, as the shape and morphology of the hydrophobic clusters are related to these secondary structures: long horizontal hydrophobic clusters denote hydrophobic or amphipathic
-helices, while short vertical or "mosaic" clusters are indicative of ß-strands. This representation is also useful in predicting membrane or transmembrane domains.
Peptide structure and conformational search.
Helical peptides were constructed with standard values of
and
angles (-60° and -40°) and subjected to a Simplex minimization (44). For the study of the putative ß-pins, conformational searches of the turn were performed using the stereoalphabet procedure (14, 36). The analysis was carried out on 8 residues and in two steps. During the first step, the (
,
) couples of the stereoalphabet were systematically attributed to the first five amino acids of the putative turn while the other residues were fixed as the ß-strand. The couples of
and
were as follows: -60°, -40° (right
-helix); -160°, 160° (ß-sheet); -140°, 80° (intermediate structure between ß-sheet and 310 helix); -80°, 160° (helix of the polyproline); -80°, 80° (310 helix); and 60°, 60° (left
-helix). From the 7,776 conformations tested, the 20 lowest-energy structures were subjected to a second analysis: the
and
angles of the last five amino acids of the putative turn were tested using all couple-of-angle values of the stereoalphabet (the conformation of the fourth and the fifth amino acids were thus studied twice). For each of the 20 first-step structures, the 10 lowest-energy structures were conserved in the second step (200 structures out of a total of 155,520 calculated conformations). The energy of these 200 structures was minimized using the Simplex method (44). The structure with the lowest energy was considered to be the most-probable structure. The energy of all structures was calculated as the sum of the torsion potential and the Van der Waals, electrostatic, and hydrophobic interactions (35). All calculations take into account the presence of a hydrophobic/hydrophilic interface with the concomitant variation of the dielectric constant (5).
MHP.
The molecular hydrophobicity potentials (MHP) were used to describe the hydrophobic or hydrophilic environment of peptides as previously described (6). For the MHP calculation, the assumption is made that the hydrophobic potential decreases exponentially with the distance as follows:
![]() |
Orientation towards a hydrophobic-hydrophilic interface.
A hydrophilic center is calculated using the following equation (5):
![]() |
i and Etri are the coordinate and the transfer energy of atom i when only the hydrophilic atoms of the molecule are taken into consideration (phi). A hydrophobic center (
pho) is calculated using a similar equation taking into account only the hydrophobic atoms. The
=
-
vector is stated as perpendicular to the interface. The interface position (
) is given by the equation:
![]() |
IMPALA.
In the IMPALA procedure (18), an empirical function, C, is used to describe the water-lipid interface. Considering that the properties of membranes are constant in the plane of the bilayer, C varies only along the z axis, perpendicular to that plane. The z axis has its origin at the center of the membrane, and the empirical function C is described by the following equation:
![]() |
and z0 are defined so that C
1 for z < -18.5 and z > 18.5 and C
0 for -13.5 < z < 13.5. The interactions between the peptide and the bilayer are described by the sum of two restraint functions as follows:
![]() |
![]() |
|
|
|---|
![]() View larger version (63K): [in a new window] |
FIG. 1. (A) HCA plot of the HCV-BEINNX11 E1 glycoprotein. (B) Minimal (dashed line), mean (solid line), and maximal (dotted line) values of the Kyte and Doolittle hydrophobicity along the E1 sequences of 25 HCV strains obtained by using a seven-residue window. Amino acid numbering refers to the position of the amino acids in the mature proteins.
|
![]() View larger version (69K): [in a new window] |
FIG. 2. (A) HCA plot of the HCV-BEINNX11 E2 glycoprotein. (B) Minimal (dashed line), mean (solid line), and maximal (dotted line) values of the Kyte and Doolittle hydrophobicity along the E2 sequences of 25 HCV strains obtained by using a seven-residue window. Amino acid numbering refers to the position of the amino acids in the mature proteins.
|
In the E1 C-terminal domain, several stretches with conserved properties were identified from the C- to the N-terminal extremities (Fig. 3A) as follows. (i) The Val180-Ala188 stretch contains only hydrophobic residues and ends with an aromatic residue (Phe187). In addition, hydrophilic, charged, and coil-favorable residues are located on both sides (Gly175, Asn176, Lys179, Gly189, Asp191, and Gly192). (ii) The Trp162-Tyr171 region has no or few charged and polar residues (a lysine is present in one sequence). It is edged by aromatic residues on both sides (Trp162, Tyr170, and Tyr171) and preceded by a hydrophilic stretch containing a coil-favorable residue (Gly159). (iii) The Ala140-Ala158 region is rather amphiphilic. It contains hydrophobic and charged or polar residues and is preceded by a stretch containing polar residues that are often present in coil or turn structures (Pro137-Thr139).
![]() View larger version (68K): [in a new window] |
FIG.3. Sequence alignment of E1 (A) and E2 (B) C-terminal domains from 25 HCV strains. A dash indicates an amino acid identical to the amino acid of HCV-BEINNX11 strain in the alignment. Amino acids are indicated according to the following characteristics (13, 60). Hydrophobic amino acids (Val, Leu, Ile, Met, Phe, Tyr, and Trp) which are preferentially implicated in transmembrane stretches are hatched. Among these, the aromatic residues (Phe, Tyr, and Trp) prefer the water-membrane interface. Charged amino acids (Asp, Glu, Arg, and Lys) which are avoided in transmembrane stretches are highlighted on a black background. His, Asn, and Gln are polar amino acids and are generally located in the water phase or at the hydrophobic-hydrophilic interface. Pro, the well-known helix breaker, and Gly, Ser, and Thr are frequent in coil and turn structures. They are indicated in italic type. Ala can be both in the membrane and the water phase. Domains 1, 2, and 3 for E1 and domains 1, 2, 3, and 4 for E2 are discussed in the text. Positions of alanine inserted by Op De Beeck et al. (45) are indicated by arrows on the top of each alignment. Stars denote the positions of alanine which have a significant effect on heterodimerization.
|
Based on these analyses, we suggest that the C-terminal domains of the E1 and E2 glycoproteins share a similar pattern (Fig. 3 and 4): a long amphiphilic region (18 or 19 residues) followed by two short hydrophobic stretches (7 to 10 residues) surrounded by hydrophilic residues.
![]() View larger version (39K): [in a new window] |
FIG. 4. Comparison of the HCA plots of the C-terminal domains of HCV-BEINNX11 E1 glycoprotein (A) and HCV-BEINNX11 E2 glycoprotein (B). (C) Schematic representation of the topological model. Putative transmembrane ß-strands (1 and 2) and a putative amphipathic -helix (3) are indicated.
|
-helices. The hydrophobic stretches of E1 and E2 are 7 to 10 residues long. This is inappropriate for classical transmembrane
-helices. In other words, since the mean distance between two consecutive residues of an
-helix is equal to 1.5 Å in the direction of the helix axis (13), these stretches should be 10.5 to 15 Å as
-helices. This is only half the width of a lipid bilayer. In a ß-strand, the distance between two adjacent amino acids is 3.4 Å (13). Therefore, the E1 and E2 stretches would be 23.8 Å long for the shortest strand and 34 Å for the longest strand. This type of configuration seems, at least arithmetically, to be the most likely to cross the membrane (Fig. 4C).
In the amphiphilic region before the putative transmembrane ß-strands, the HCA profile shows long horizontal hydrophobic clusters suggestive of amphipathic
-helices (21) (Fig. 4).
In order to test the hypotheses, several 3-D models were built for the HCV-BEINNX11 strain. These models were tested for membrane interaction.
3-D model of amphipathic
-helices.
Helices were built for every stretch of a minimal length of 13 amino acids between Ala140 and Ala158 for E1 and between Ile313 and Ile331 for E2. The energy of each
-helix was minimized using the Simplex method. Each structure was then oriented towards a hydrophobic-hydrophilic interface. Helices eliciting a tilt of less than 30° between their axis and this interface were considered to be amphipathic and were selected for a thorough study using the IMPALA procedure. The 10 IMPALA Monte Carlo simulations gave similar results for each peptide. The MHP were also calculated to check the hydropathy of each structure. The best candidate as amphipathic
-helix was selected as the stretch with a maximal length, an insertion angle close to 0°, and a maximal depth of penetration in the bilayer.
Table 1
summarizes the screening for putative amphipathic
-helices in E1. Most peptides behave as amphipathic
-helices, and Val142-Ala158 was selected as the best candidate (Fig. 5). This 17-residues stretch makes a 10° angle with respect to its hydrophobic-hydrophilic interface and with the IMPALA membrane surface (Fig. 5C). In this position, the mass center of the structure was located at 13.5 Å from the membrane center, indicating that the peptide is able to adsorb at the membrane surface by inserting between the phospholipid heads. The MHP confirm the amphipathy of the helix since the hydrophobic and hydrophilic isopotential envelopes are well segregated at the helix surface (Fig. 5B).
|
View this table: [in a new window] |
TABLE 1. Characterization of helical peptides found between Ala140 and Ala158a
|
![]() View larger version (53K): [in a new window] |
FIG. 5. HCV-BEINNX11 E1 putative amphipathic -helix (Val142 to Ala158). (A) Minimized 3-D structure. (B) MHP around the peptide (orange = hydrophobic; green = hydrophilic). (C) Best configuration found by the IMPALA procedure. The center of the phospholipid bilayer is indicated by an orange line, the limit between the phospholipid acyl chains and polar heads is indicated by a yellow line, and the limit between the water phase and the phospholipids is indicated by a green line.
|
-helices are amphipathic. The 16-amino-acid Val316-Ile331 peptide was selected as the best candidate, as its angle with respect to its hydrophobic-hydrophilic interface or with the IMPALA modeled lipid bilayer surface is close to 0° (Fig. 6). The mass center of the helix in the membrane is at 12.5 Å from the membrane center. The MHP confirm the segregation of the hydrophobic and hydrophilic faces separately (Fig. 6B). |
View this table: [in a new window] |
TABLE 2. Characterization of helical peptides found between Ile313 and Ile331a
|
![]() View larger version (58K): [in a new window] |
FIG. 6. HCV-BEINNX11 E2 putative amphipathic -helix (Val316 to Ile331). (A) Minimized 3-D structure. (B) MHP around the peptide (orange = hydrophobic; green = hydrophilic). (C) Best configuration found by the IMPALA procedure. The center of the phospholipid bilayer is indicated by an orange line, the limit between the phospholipid acyl chains and polar heads is indicated by a yellow line, and the limit between the water phase and the phospholipids is indicated by a green line.
|
The lowest energy structure obtained for E1 is shown in Fig. 7A. The analysis shows that the best conformation of the loop corresponds to an antiparallel pair of ß-strands (Trp162 to Met173 and Trp177 to Ala188). Interstrand hydrogen bonds are implicated in the energy stabilization of the structure. MHP indicate that the final structure is highly hydrophobic except at the backbone contact (due to accessible main-chain C
O and N-H groups), in the turn, and at the peptide N- and C- terminal extremities (Fig. 7B). The hydrophobic part is longer than 35 Å, and this should be sufficient to cross a membrane. The results of the 10 IMPALA simulations confirmed this hypothesis. The most favorable position of this structure is transmembranous (Fig. 7C), with the soluble and adsorbed states being unlikely.
![]() View larger version (80K): [in a new window] |
FIG. 7. HCV-BEINNX11 E1 putative transmembrane ß-strands (Trp162 to Ala188). (A) Best minimized 3-D structure obtained by the stereoalphabet procedure. (B) MHP around the peptide (orange = hydrophobic; green = hydrophilic). (C) Best configuration found by the IMPALA procedure. The center of the phospholipid bilayer is indicated by an orange line, the limit between the phospholipid acyl chains and polar heads is indicated by a yellow line, and the limit between the water phase and the phospholipids is indicated by a green line.
|
![]() View larger version (59K): [in a new window] |
FIG. 8. HCV-BEINNX11 E2 putative transmembrane ß-strands (Tyr335 to Ala359). (A) Best minimized 3-D structure obtained by the stereoalphabet procedure. (B) MHP around the peptide (orange = hydrophobic; green = hydrophilic). (C) Best configuration found by the IMPALA procedure. The center of the phospholipid bilayer is indicated by an orange line, the limit between the phospholipid acyl chains and polar heads is indicated by a yellow line, and the limit between the water phase and the phospholipids is indicated by a green line.
|
|
|
|---|
-helices. Even if this method has already been tried and tested, it has limited utility for the study of nonhelical transmembrane domains or of membrane spans with unusual properties. This is the case for HCV glycoproteins which contain hydrophobic stretches shorter than 18 residues.
In the first part of our investigation, the analysis of E1 and E2 C-terminal sequences of 25 HCV strains showed that both envelope proteins share a similar pattern: two short hydrophobic stretches separated by charged or polar amino acids and, in the case of E1, by coil-favorable residues. These stretches follow a region likely to be an amphipathic helix. We postulate first that the E1 and E2 proteins could be anchored in the ER membrane by the two short hydrophobic C-terminal stretches. Stereoalphabet analysis, MHP calculation, and IMPALA simulations indicate that these stretches could form transmembrane ß-pins. Second, these ß-pins could interact more easily with the lipid bilayer because of the adsorption of the amphipathic
-helix. These two hypotheses are discussed below.
The E1 and E2 glycoproteins are cleaved from the HCV polyprotein by host signal peptidase (16, 23, 31, 48). Both extremities of the proteins are thus in the ER lumen, at least transiently, and this is the reason why an even number of transmembrane stretches is expected. Our first hypothesis is that the two required transmembrane spans are the short hydrophobic C-terminal segments of E1 and E2. Another configuration would present a significant number of hydrophobic residues in the water phase and charged and polar residues in the membrane. Moreover, the proposed configuration favorably places the aromatic residues at the water-phospholipid interface.
The E1 and E2 C-terminal domains signal for static retention in the ER (9, 10, 20). CD4 and CD8 are usually targeted to the plasma membrane. If the last 31 amino acids of E1 are coupled to the CD4 or CD8 ectodomain, the chimeric protein remains in the ER compartment (9). Cocquerel et al. (10) have also shown that E2 is retained in the ER by its transmembrane domain (the 29 C-terminal residues) and that this segment is able to confer the ER-retention property to the CD4 ectodomain. These "cut-and-paste" experiments prove the ER retention properties of the E1 31 C-terminal and the E2 29 C-terminal residues. Moreover, they demonstrate that these fragments are independent entities with native folding, which supports the idea that they are the E1 and E2 transmembrane segments, as proposed in our model.
It is generally admitted that the hydrocarbon region of a fluid bilayer is typically 30 Å thick (34, 58). Analyzing the cytochrome c oxidase from bovine mitochondria, Wallin et al. (57), however, showed that this bundle of 28 transmembrane
-helices extends about ±20 Å from the center of the membrane. They further emphasized that the hydrophobic residues are within ±10 Å from the bilayer center and that this slice is symmetrically flanked by regions rich in aromatic residues. In 1999, Monne et al. (43) determined that the minimum length of a poly-Leu helical hairpin is 31 residues. The shortest helical hairpin would thus be made by two
-helices of 14 amino acids separated by a three-residue turn. Since the standard translation per residue in an
-helix is 1.5 Å (13), the hydrophobic portion of these helices would be 21 Å. Taking these results into consideration, the hydrophobic core of a membrane would thus be at least
20 Å thick. As specified earlier, the E1 and E2 hydrophobic stretches are 7 to 10 residues long (aromatic residues included). An
-helix structure thus appears to be inappropriate to cross the membrane. In a ß-strand conformation, the length of these stretches would be 23.8 or 34 Å. These lengths are adapted to a transmembrane topology and are consistent with the size of the strands of the porin ß-barrels, which are 6 to 17 amino acids long. MHP confirm this hypothesis, since the transmembrane ß-pins exhibit a hydrophobic span at least equal to 35 Å.
Extended structures are thought to be energetically unfavorable in a membrane environment because of the accessibility of main-chain polar groups (60). However, some peptides form ß-sheet structures when associated to a lipid bilayer (1, 28, 32, 37). Hexapeptide made of tryptophan and leucine (acetyl-Trp-Leu5 peptides) has been shown to assemble into oligomers of 10 to 20 antiparallel ß-strands in membranes (59). PATIR-FTIR measurements indicate that interstrand hydrogen bonds of this sheet are parallel to the membrane surface. The value of the dichroic ratio suggests that the ß-sheet has a perpendicular orientation. When the length of the Ac-Trp-Leun peptides is increased to nine residues, ß-sheets are still observed in POPC bilayers. The critical number of leucine residues for helix formation in a membrane should be between 10 and 14. These results reveal that the formation of ß-sheet in a membrane is easier than expected and thus support our prediction.
The calculation of different stretches from the amphipathic region of E1 and E2 supports our second hypothesis. Indeed, most helices are amphipathic, with a clear segregation of hydrophobic and hydrophilic sides, and adsorb onto the IMPALA bilayer. The Val142-to-Ala158 fragment of E1, and the Val316-to-Ile331 fragment of E2 are the best candidates. Since these amphipathic
-helices are very close to the transmembrane domains, they must be close to the viral membrane as well and could thus probably adsorb to it. We suggest that the membrane adsorption of these amphipathic
-helices stabilizes the membrane insertion of the E1 and E2 ß-pins. Characterization of truncated forms of E1 and E2 indicate that the truncation of the putative transmembrane domains is not sufficient to observe efficient secretion of products with native folding (26, 39, 41, 42, 53). This could be correlated to the presence of other membrane-interacting structures, such as adsorbed
-helices.
Our model (Fig. 9) is compatible with available experimental data. Several features explain the retention of proteins in the ER (reviewed in reference 22). "KDEL," di-lysine, or di-arginine patterns have been identified. Another signal for ER retention is the occurrence of one or more hydrophilic residues in the middle of a hydrophobic transmembrane domain (4, 33, 50). Recent studies suggest that more than one signal is required. The participation of both soluble and transmembrane domains has been illustrated with cytochrome P450 (56), cytochrome b5 (25), Sec12p (52), and the rubella virus E1 glycoprotein (24). Honsho et al. (25) showed that mutation of all or part of the transmembrane domain of cytochrome b5 with hydrophobic residues has almost no effect on the intracellular distribution. However, an increase in the length of the transmembrane domain over 22 amino acids suppresses the retention. Thus, the transmembrane domain length directly functions as a retention signal of cytochrome b5. These results are consistent with those of Pedrazzini et al. (46), who showed that the insertion of five amino acids into the transmembrane domain of cytochrome b5 resulted in cell surface localization. The length of the transmembrane domain is critical for retention of other proteins in the ER (49, 62) and in the Golgi apparatus (reviewed in references 12 and 22). According to these results, our model suggests that E1 and E2 anchor domains contain a bipartite signal responsible for ER retention. The length of E1 and E2 transmembrane stretches should be the first critical determinant for retention. Since there is a relation between the length and the structure of transmembrane spans (see above), we further propose that a ß structure could be correlated to the ER retention properties of these anchor domains. This would also be consistent with Flint and McKeating's results (20), which suggest that the HCV E1 and E2 glycoproteins are located in the ER because of the structural properties of their transmembrane segments rather than of specific patterns of sequences. The other part of the retention signal should be in the turn between the two ß-strands. Indeed, these stretches contain charged residues which are crucial for ER retention (11).
![]() View larger version (21K): [in a new window] |
FIG. 9. Topological models of the E1 and E2 C-terminal anchor domains.
|
The final topological models of E1 and E2 anchor domains thus consist of an adsorbed
-helix followed by two interacting transmembrane ß-strands (Fig. 9). This partly contradicts the results from Op De Beeck et al. (45), who used alanine insertion to demonstrate that the E1 and E2 C-terminal domains are implicated in E1-E2 dimerization. They further analyzed the 3-D structure of the half of the E1 putative transmembrane domain in membrane mimetic solvents by circular dichroism (CD) and nuclear magnetic resonance (NMR). CD spectra in 2,2,2-trifluoroethanol (TFE), sodium dodecyl sulfate, and lysophosphatidylcholine indicated an
-helical content equal to 48, 44, and 52%, respectively. In TFE, the 21-residue-long peptide (vol/vol) is more frequently helical as determined by NMR analysis.
If the alanine insertion (Fig. 3) clearly proves that the fragment is involved in heterodimerization, it does not prove which structure is required for that dimerization. Indeed, the insertion of one residue in the middle of a secondary-structure element would perturb the spatial distribution of the other residues in an
-helix as well as in a ß-strand. In both cases, the perturbation could prevent the dimerization. According to our model, alanine insertion should prevent heterodimerization if the insertion modifies the structure of a ß-strand or if it generates unfavorable cross-strand pairs of amino acids. Mutant proteins unable to form heterodimers could also occur if an alanine modifies the structure of the turn of one ß-pin or the interactions between the turns. The addition of one alanine in a ß-strand should not modify its secondary structure since the length of the transmembrane stretch should not be sufficiently increased to form an
-helix (see above). The insertion of alanine in a ß-strand should thus only affect favorable or unfavorable cross-strand pairs of residues. Only three of the alanines inserted in the ß-strands have an effect on heterodimerization which is reduced to less than 10% of the wild-type level (45). They are located in the first ß-strand of E1 (Fig. 3). Alanines inserted in the other putative ß-strands do not have such a drastic effect. Interestingly, the first ß-strand of E1 has fewer hydrophobic residues than the others (Fig. 3). It contains two glycines, both on the same side of the strand. This peculiar composition could explain its sensitivity to alanine insertion. Indeed, glycine has particular requirements related to cross-strands pairs (40). On the other hand, the introduction of one alanine in the other strands should not drastically modify the interactions between strands since they contain exclusively hydrophobic amino acids (Fig. 3). In the case of alanine insertion close to the turn between E1 and E2 strands, three mutants are impaired in heterodimerization (Fig. 3). The effect is less dramatic for E2, since the heterodimerization is reduced to about 35% of the wild-type level. This could be due to a modification of the turn structure, less compatible with the interaction of the two ß-pins.
CD and NMR analyses by Op de Beeck et al. tend to prove that the E1 transmembrane domain is one
-helix (45). These results must be considered with caution since they were not obtained with the complete E1 putative transmembrane domain. Furthermore, the use of membrane mimetic solvents could have induced an artifact, as TFE is a well-known helix-inducing solvent (3, 27, 37, 54). Finally, the model of Op De Beeck et al. does not take into account the requirement of crossing the membrane twice. They indeed suggest that the E1 and E2 transmembrane domains would first form hairpin structures and then reorganize into single transmembrane helices after the processing of E1 and E2 proteins. We support the argument that the E1 and E2 transmembrane domains adopt a ß-hairpin configuration and that this configuration could be more stable than expected.
Conclusions.
In addition to their role in membrane association, the E1 and E2 transmembrane domains (the 31 and 29 last amino acids, respectively) are responsible for heterodimerization and retention of the two HCV envelope glycoproteins in the ER compartment. This is crucial for the HCV life cycle. We put forward the hypothesis that the E1 and E2 transmembrane domains adopt, at least transiently, a transmembrane ß-pin conformation downstream of an amphipathic
-helix. The ß-structure, the multimerization, and the ER retention properties of the transmembrane stretches would be correlated to their length. To our knowledge, this model is in agreement with results from the literature. Further experimental results will be needed to validate or invalidate our model.
We thank the National Funds for Scientific Research of Belgium (FNRS), where B.C. is research fellows, L.L. is research associate, and R.B. is research director.
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»