vpu transmembrane peptide structure obtained by site-specific fourier transform infrared dichroism and global molecular dynamics searching.

The recently developed method of site-directed Fourier transform infrared dichroism for obtaining orientational constraints of oriented polymers is applied here to the transmembrane domain of the vpu protein from the human immunodeficiency virus type 1 (HIV-1). The infrared spectra of the 31-residue-long vpu peptide reconstituted in lipid vesicles reveal a predominantly alpha-helical structure. The infrared dichroism data of the (13)C-labeled peptide yielded a helix tilt beta = (6.5 +/- 1.7) degrees from the membrane normal. The rotational pitch angle omega, defined as zero for a residue located in the direction of the helix tilt, is omega = (283 +/- 11) degrees for the (13)C labels Val(13)/Val(20) and omega = (23 +/- 11) degrees for the (13)C labels Ala(14)/Val(21). A global molecular dynamics search protocol restraining the helix tilt to the experimental value was performed for oligomers of four, five, and six subunits. From 288 structures for each oligomer, a left-handed pentameric coiled coil was obtained, which best fits the experimental data. The structure reveals a pore occluded by Trp residues at the intracellular end of the transmembrane domain.     

The transmembrane domain of the oncogenic mutant ErbB-2 receptor: a structure obtained from site-specific infrared dichroism and molecular dynamics

ErbB-2 is a member of the family of epidermal growth factor receptors, which shows an oncogenic mutation in the rat gene neu, Val664Glu in the transmembrane domain that causes permanent dimerisation and subsequently leads to uncontrollable cell division and tumour formation. We have obtained the alpha-helical structure of the mutant transmembrane domain dimer experimentally with site-specific infrared dichroism (SSID) based on six transmembrane peptides with 13C18O carbonyl group-labelled residues. The derived orientational data indicate a local helix tilt ranging from 28(+/-6) degrees to 22(+/-4) degrees. Altogether using orientational constraints from SSID and experimental alpha-helical constraints while performing a systematic conformational search including molecular dynamics simulation in a lipid bilayer, we have obtained a unique experimentally defined atomic structure. The resulting structure consists of a right handed alpha-helical bundle with the residues Ile659, Val663, Leu667, Ile671, Val674 and Leu679 in the dimerisation interface. The right-handed bundle is in contrast to the left-handed structures obtained in previous modelling efforts. In order to facilitate tight helical packing, the spacious Glu664 residues do not interact directly but with water molecules that enter the bilayer.     

The influenza C virus CM2 protein can alter intracellular pH, and its transmembrane domain can substitute for that of the influenza A virus M2 protein and support infectious virus production

The influenza C virus CM2 protein and a chimeric influenza A virus M2 protein (MCM) containing the CM2 transmembrane domain were assessed for their ability to functionally replace the M2 protein. While all three proteins could alter cytosolic pH to various degrees when expressed from cDNA, only M2 and MCM could at least partially restore infectious virus production to M2-deficient influenza A viruses. The data suggest that while the CM2 ion channel activity is similar to that of M2, sequences in the extracellular and/or cytoplasmic domains play important roles in infectious virus production.     

Role of the CM2 protein in the influenza C virus replication cycle

CM2 is the second membrane protein of influenza C virus. Although its biochemical characteristics, coding strategy, and properties as an ion channel have been extensively studied, the role(s) of CM2 in the virus replication cycle remains to be clarified. In order to elucidate this role, in the present study we generated CM2-deficient influenza C virus-like particles (VLPs) and examined the VLP-producing 293T cells, VLPs, and VLP-infected HMV-II cells. Quantification of viral RNA (vRNA) in the VLPs by real-time PCR revealed that the CM2-deficient VLPs contain approximately one-third of the vRNA found in wild-type VLPs although no significant differences were detected in the expression levels of viral components in VLP-producing cells or in the number and morphology of the generated VLPs. This finding suggests that CM2 is involved in the genome packaging process into VLPs. Furthermore, HMV-II cells infected with CM2-deficient VLPs exhibited significantly reduced reporter gene expression. Although CM2-deficient VLPs could be internalized into HMV-II cells as efficiently as wild-type VLPs, a smaller amount of vRNA was detected in the nuclear fraction of CM2-deficient VLP-infected cells than in that of wild-type VLP-infected cells, suggesting that the uncoating process of the CM2-deficient VLPs in the infected cells did not proceed in an appropriate manner. Taken together, the data obtained in the present study indicate that CM2 has a potential role in the genome packaging and uncoating processes of the virus replication cycle.     

Modeling the structure of the respiratory syncytial virus small hydrophobic protein by silent-mutation analysis of global searching molecular dynamics

Human respiratory syncytial virus (RSV) encodes a small hydrophobic (SH) protein, whose function in the life cycle of the virus is unknown. Recent channel activity measurements of the protein suggest that like other viroporins, SH may assemble into a homo-oligomeric ion channel. To further our understanding of this potentially important protein, a new strategy was implemented in order to model the transmembrane oligomeric bundle of the protein. Global searching molecular dynamic simulations of SH proteins from eight different viral strains, each at different oligomeric states, as well as different lengths of the putative transmembrane domain, were undertaken. Taken together, a total of 45 different global molecular dynamic simulations pointed to a single pentameric structure for the protein that was found in all of the different variants. The model of the structure obtained is a channel-like homopentamer whose minimal transmembrane pore diameter is 1.46 A.     

Global fold and backbone dynamics of the hepatitis C virus E2 glycoprotein transmembrane domain determined by NMR

E1 and E2 are two hepatitis C viral envelope glycoproteins that assemble into a heterodimer that is essential for membrane fusion and penetration into the target cell. Both extracellular and transmembrane (TM) glycoprotein domains contribute to this interaction, but study of TM–TM interactions has been limited because synthesis and structural characterization of these highly hydrophobic segments present significant challenges. In this NMR study, by successful expression and purification of the E2 transmembrane domain as a fusion construct we have determined the global fold and characterized backbone motions for this peptide incorporated in phospholipid micelles. Backbone resonance frequencies, relaxation rates and solvent exposure measurements concur in showing this domain to adopt a helical conformation, with two helical segments spanning residues 717–726 and 732–746 connected by an unstructured linker containing the charged residues D728 and R730 involved in E1 binding. Although this linker exhibits increased local motions on the ps timescale, the dominating contribution to its relaxation is the global tumbling motion with an estimated correlation time of 12.3 ns. The positioning of the helix–linker–helix architecture within the mixed micelle was established by paramagnetic NMR spectroscopy and phospholipid-peptide cross relaxation measurements. These indicate that while the helices traverse the hydrophobic interior of the micelle, the linker lies closer to the micelle perimeter to accommodate its charged residues. These results lay the groundwork for structure determination of the E1/E2 complex and a molecular understanding of glycoprotein heterodimerization.     

Contribution of energy values to the analysis of global searching molecular dynamics simulations of transmembrane helical bundles

Molecular interactions between transmembrane alpha-helices can be explored using global searching molecular dynamics simulations (GSMDS), a method that produces a group of probable low energy structures. We have shown previously that the correct model in various homooligomers is always located at the bottom of one of various possible energy basins. Unfortunately, the correct model is not necessarily the one with the lowest energy according to the computational protocol, which has resulted in overlooking of this parameter in favor of experimental data. In an attempt to use energetic considerations in the aforementioned analysis, we used global searching molecular dynamics simulations on three homooligomers of different sizes, the structures of which are known. As expected, our results show that even when the conformational space searched includes the correct structure, taking together simulations using both left and right handedness, the correct model does not necessarily have the lowest energy. However, for the models derived from the simulation that uses the correct handedness, the lowest energy model is always at, or very close to, the correct orientation. We hypothesize that this should also be true when simulations are performed using homologous sequences, and consequently lowest energy models with the right handedness should produce a cluster around a certain orientation. In contrast, using the wrong handedness the lowest energy structures for each sequence should appear at many different orientations. The rationale behind this is that, although more than one energy basin may exist, basins that do not contain the correct model will shift or disappear because they will be destabilized by at least one conservative (i.e. silent) mutation, whereas the basin containing the correct model will remain. This not only allows one to point to the possible handedness of the bundle, but can be used to overcome ambiguities arising from the use of homologous sequences in the analysis of global searching molecular dynamics simulations. In addition, because clustering of lowest energy models arising from homologous sequences only happens when the estimation of the helix tilt is correct, it may provide a validation for the helix tilt estimate.     

Structure prediction of the dimeric neu/ErbB-2 transmembrane domain from multi-nanosecond molecular dynamics simulations

Dimerization of the neu/ErbB-2 receptor tyrosine kinase is a necessary but not a sufficient step for signaling. Despite the efforts expended to identify the molecular interactions responsible for receptor-receptor contacts and particularly those involving the transmembrane domain, structural details are still unknown. In this work, molecular dynamics simulations of the helical transmembrane domain (TM) of neu and ErbB-2 receptors are used to predict their dimer structure both in the wild and oncogenic forms. A global conformational search method, applied to define the best orientations of parallel helices, showed an energetically favorable configuration with the specific mutation site within the interface, common for both the nontransforming and the transforming neu/ErbB-2 TM dimers. Starting from this configuration, a total of 10 simulations, about 1.4 ns each, performed in vacuum, without any constraints, show that the two helices preferentially wrap in left-handed interactions with a packing angle at about 20°. The resulting structures are nonsymmetric and the hydrogen bond network analysis shows that helices experience π local distortions that facilitate inter-helix hydrogen bond interactions and may result in a change in the helix packing, leading to a symmetric interface. For the mutated sequences, we show that the Glu side chain interacts directly with its cognate or with carbonyl groups of the facing backbone. We show that the connectivity between interfacial residues conforms to the knobs-into-holes packing mode of transmembrane helices. The dimeric interface described in our models is discussed with respect to mutagenesis studies. Received: 12 March 1999 / Revised version: 23 August 1999 / Accepted: 23 August 1999     

Defining the transmembrane helix of M2 protein from influenza A by molecular dynamics simulations in a lipid bilayer

Integral membrane proteins containing at least one transmembrane (TM) alpha-helix are believed to account for between 20% and 30% of most genomes. There are several algorithms that accurately predict the number and position of TM helices within a membrane protein sequence. However, these methods tend to disagree over the beginning and end residues of TM helices, posing problems for subsequent modeling and simulation studies. Molecular dynamics (MD) simulations in an explicit lipid and water environment are used to help define the TM helix of the M2 protein from influenza A virus. Based on a comparison of the results of five different secondary structure prediction algorithms, three different helix lengths (an 18mer, a 26mer, and a 34mer) were simulated. Each simulation system contained 127 POPC molecules plus approximately 3500-4700 waters, giving a total of approximately 18,000-21,000 atoms. Two simulations, each of 2 ns duration, were run for the 18mer and 26mer, and five separate simulations were run for the 34mer, using different starting models generated by restrained in vacuo MD simulations. The total simulation time amounted to 11 ns. Analysis of the time-dependent secondary structure of the TM segments was used to define the regions that adopted a stable alpha-helical conformation throughout the simulation. This analysis indicates a core TM region of approximately 20 residues (from residue 22 to residue 43) that remained in an alpha-helical conformation. Analysis of atomic density profiles suggested that the 18mer helix revealed a local perturbation of the lipid bilayer. Polar side chains on either side of this region form relatively long-lived H-bonds to lipid headgroups and water molecules.     

Stability and membrane orientation of the fukutin transmembrane domain: a combined multiscale molecular dynamics and circular dichroism study

The N-terminal domain of fukutin-I has been implicated in the localization of the protein in the endoplasmic reticulum and Golgi Apparatus. It has been proposed to mediate this through its interaction with the thinner lipid bilayers found in these compartments. Here we have employed multiscale molecular dynamics simulations and circular dichroism spectroscopy to explore the structure, stability, and orientation of the short 36-residue N-terminus of fukutin-I (FK1TMD) in lipids with differing tail lengths. Our results show that FK1TMD adopts a stable helical conformation in phosphatidylcholine lipids when oriented with its principal axis perpendicular to the bilayer plane. The stability of the helix is largely insensitive to the lipid tail length, preventing hydrophobic mismatch by virtue of its mobility and ability to tilt within the lipid bilayers. This suggests that changes in FK1TMD tilt in response to bilayer properties may be implicated in the regulation of its trafficking. Coarse-grained simulations of the complex Golgi membrane suggest the N-terminal domain may induce the formation of microdomains in the surrounding membrane through its preferential interaction with 1,2-dipalmitoyl-sn-glycero-3-phosphatidylinositol 4,5-bisphosphate lipids.     

Characterization of a second protein (CM2) encoded by RNA segment 6 of influenza C virus.

The biochemical properties of a second protein (CM2) encoded by RNA segment 6 of influenza C virus were investigated. Three forms of CM2 with different electrophoretic mobilities (CM2(0), CM2a, and CM2b) were detected in infected cells by immunoprecipitation with antiserum to the glutathione S-transferase (GST)-CM2 fusion protein. Treatment of infected cells with tunicamycin and digestion of immunoprecipitated proteins with endoglycosidase H or peptide-N-glycosidase F suggested that a mannose-rich oligosaccharide core is added to unglycosylated CM2(0) (Mr, approximately 16,000) to form CM2a (Mr, approximately 18,000) and that the processing of the carbohydrate chain from the high-mannose type to the complex type converts CM2a into CM2b, which is heterogeneous in electrophoretic mobility (Mr, approximately 22,000 to 30,000). Labeling of infected cells with [3H]palmitic acid showed that CM2 is fatty acylated. The fatty acid bond was sensitive to treatment with hydroxylamine and mercaptoethanol, which indicates a labile thioester-type linkage. The CM2 protein was also found to form disulfide-linked dimers and tetramers on sodium dodecyl sulfate-polyacrylamide gels under nonreducing conditions. Trypsin treatment of infected cell surfaces as well as of microsome vesicles from infected cells followed by immunoprecipitation with antiserum to the GST fusion protein containing the 56 C-terminal amino acid residues of CM2 suggested that this C-terminal domain is intracellular and exposed to the cytoplasms of microsomes. Furthermore, evidence that a small amount of CM2 is incorporated into progeny virus particles was obtained by Western blot analysis. These results, altogether, suggest that CM2 is an integral membrane protein with biochemical properties similar to those of influenza A virus M2 and influenza B virus NB proteins.     

The effect of the cytoplasmic tail of influenza C virus CM2 protein on its biochemical properties and intracellular processing

CM2 is an integral membrane protein encoded by the influenza C virus M gene. To examine the effects of the cytoplasmic tail of CM2 on its biochemical properties, deletion and substitution mutations were introduced into CM2 cytoplasmic tail at residues 47–115, and the expressed CM2 mutants were investigated. Although the cytoplasmic tail is not essential for the oligomerization of CM2, it may affect the degree of oligomerization. The residues 47–48, 67–69, 73–90 and 113–115 were all required for the proper expression of CM2. Pulse-chase experiments suggest that residues 47–48, 67–69, 73–75 and 79–87 stabilize CM2, thereby affecting CM2 expression. The C-terminal region at residues 61–115 is not essential for CM2 transport to the cell surface, and a 14-amino-acid, but not an 11-amino-acid, cytoplasmic tail is sufficient for the cell surface expression of CM2. These results suggest that either certain amino acid sequences or the length of the CM2 cytoplasmic tail are necessary for the proper conformational maturation, stability, expression level and intracellular transport of CM2.     

Structure and functional dynamics characterization of the ion channel of the human respiratory syncytial virus (hRSV) small hydrophobic protein (SH) transmembrane domain by combining molecular dynamics with excited normal modes

The human respiratory syncytial virus (hRSV) is the major cause of lower respiratory tract infection in children and elderly people worldwide. Its genome encodes 11 proteins including SH protein, whose functions are not well known. Studies show that SH protein increases RSV virulence degree and permeability to small compounds, suggesting it is involved in the formation of ion channels. The knowledge of SH structure and function is fundamental for a better understanding of its infection mechanism. The aim of this study was to model, characterize, and analyze the structural behavior of SH protein in the phospholipids bilayer environment. Molecular modeling of SH pentameric structure was performed, followed by traditional molecular dynamics (MD) simulations of the protein immersed in the lipid bilayer. Molecular dynamics with excited normal modes (MDeNM) was applied in the resulting system in order to investigate long time scale pore dynamics. MD simulations support that SH protein is stable in its pentameric form. Simulations also showed the presence of water molecules within the bilayer by density distribution, thus confirming that SH protein is a viroporin. This water transport was also observed in MDeNM studies with histidine residues of five chains (His22 and His51), playing a key role in pore permeability. The combination of traditional MD and MDeNM was a very efficient protocol to investigate functional conformational changes of transmembrane proteins that act as molecular channels. This protocol can support future investigations of drug candidates by acting on SH protein to inhibit viral infection.

Open image in new window

Graphical Abstract The ion channel of the human respiratory syncytial virus (hRSV) small hydrophobic protein (SH) transmembrane domain?

     

Exploring conformational states and helical packings in the P2X receptor transmembrane domain by molecular dynamics simulation

The P2X receptor is a trimeric transmembrane protein that acts as an ATP-gated ion channel. Its transmembrane domain (TMD) contains only six helices and three of them, the M2 helices, line the ion conduction pathway. Here, using molecular dynamics simulation, I identify four conformational states of the TMD that are associated with four types of packing between M2 helices. Packing in the extracellular half of the M2 helix produces closed conformations, while packing in the intracellular half produces both open and closed conformations. State transition is observed and supports a mechanism where iris-like twisting of the M2 helices switches the location of helical packing between the extracellular and the intracellular halves of the helices. In addition, this twisting motion alters the position and orientation of residue side-chains relative to the pore and therefore influences the pore geometry and possibly ion permeation. Helical packing, on the other hand, may restrict the twisting motion and generate discrete conformational states.     

Structure determination of UL49.5 transmembrane protein from bovine herpesvirus 1 by NMR spectroscopy and molecular dynamics

The transporter associated with antigen processing (TAP) directly participates in the immune response as a key component of the cytosolic peptide to major histocompatibility complex (MHC) class I protein loading machinery. This makes TAP an important target for viruses avoiding recognition by CD8+ T lymphocytes. Its activity can be suppressed by the UL49.5 protein produced by bovine herpesvirus 1, although the mechanism of this inhibition has not been understood so far.Therefore, the main goal of our study was to investigate the 3D structure of bovine herpesvirus 1 - encoded UL49.5 protein. The final structure of the inhibitor was established using circular dichroism (CD), 2D nuclear magnetic resonance (NMR), and molecular dynamics (MD) in membrane mimetic environments. In NMR studies, UL49.5 was represented by two fragments: the extracellular region (residues 1–35) and the transmembrane-intracellular fragment (residues 36–75), displaying various functions during viral invasion. After the empirical structure determination, a molecular docking procedure was used to predict the complex of UL49.5 with the TAP heterodimer.Our results revealed that UL49.5 adopted a highly flexible membrane-proximal helical structure in the extracellular part. In the transmembrane region, we observed two short α-helices. Furthermore, the cytoplasmic part had an unordered structure. Finally, we propose three different orientations of UL49.5 in the complex with TAP. Our studies provide, for the first time, the experimental structural information on UL49.5 and structure-based insight in its mechanism of action which might be helpful in designing new drugs against viral infections.     

Structure of outer membrane protein A transmembrane domain by NMR spectroscopy

We have determined the three-dimensional fold of the 19 kDa (177 residues) transmembrane domain of the outer membrane protein A of Escherichia coli in dodecylphosphocholine (DPC) micelles in solution using heteronuclear NMR. The structure consists of an eight-stranded beta-barrel connected by tight turns on the periplasmic side and larger mobile loops on the extracellular side. The solution structure of the barrel in DPC micelles is similar to that in n-octyltetraoxyethylene (C(8)E(4)) micelles determined by X-ray diffraction. Moreover, data from NMR dynamic experiments reveal a gradient of conformational flexibility in the structure that may contribute to the membrane channel function of this protein.     

Intrinsically unstructured regions in the C domain of the influenza virus M1 protein

The M1 matrix protein of the influenza virus is one of the main structural components of the virion that performs several different functions in the infected cell. X-ray analysis (with 2.08 Å resolution) has been performed for the N-terminal part of the M1 protein (residues 2–158) but not for its C-terminal domain (159–252). In the present study, we analyzed the structure of the M1 protein of the influenza virus A/Puerto Rico/8/34 (H1N1) strain in acidic solution using tritium planigraphy. The incorporation of tritium label into the domains of the M1 protein were studied; the C domain and the interdomain loops are preferentially accessible to tritium. Analytical centrifugation and dynamic laser light scattering demonstrated anomalous hydrodynamic parameters and low structuredness of the M1 protein, which has also been confirmed by circular dichroism data. Bioinformatic analysis of the M1 protein sequence revealed intrinsically unstructured segments that were concentrated in the C domain and interdomain loops between the N-, M-, and C domains. We suggest that the multifunctionality of the M1 protein in a cell is determined by the plasticity of its tertiary structure, which is caused by the presence of intrinsically unstructured segments.     

Structural elucidation of transmembrane transporter protein bilitranslocase: Conformational analysis of the second transmembrane region TM2 by molecular dynamics and NMR spectroscopy

Membrane proteins represent about a third of the gene products in most organisms, as revealed by the genome sequencing projects. They account for up to two thirds of known drugable targets, which emphasizes their critical pharmaceutical importance. Here we present a study on bilitranslocase (BTL) (TCDB 2.A.65), a membrane protein primarily involved in the transport of bilirubin from blood to liver cells. Bilitranslocase has also been identified as a potential membrane transporter for cellular uptake of several drugs and due to its implication in drug uptake, it is extremely important to advance the knowledge about its 3D structure. However, at present, only a limited knowledge is available beyond the primary structure of BTL. It has been recently confirmed experimentally that one of the four computationally predicted transmembrane segments of bilitranslocase, TM3, has a helical structure with hydrophilic amino acid residues oriented towards one side, which is typical for transmembrane domains of membrane proteins. In this study we confirmed by the use of multidimensional NMR spectroscopy that the second transmembrane segment, TM2, also appears in a form of α-helix. The stability of this polypeptide chain was verified by molecular dynamics (MD) simulation in dipalmitoyl phosphatidyl choline (DPPC) and in sodium dodecyl sulfate (SDS) micelles. The two α-helices, TM2 corroborated in this study, and TM3 confirmed in our previous investigation, provide reasonable building blocks of a potential transmembrane channel for transport of bilirubin and small hydrophilic molecules, including pharmaceutically active compounds.