Mapping the energy surface of transmembrane helix-helix interactions

Viroporins in the Influenza Virus

Influenza is a highly contagious virus that causes seasonal epidemics and unpredictable pandemics. Four influenza virus types have been identified to date: A, B, C and D, with only A–C known to infect humans. Influenza A and B viruses are responsible for seasonal influenza epidemics in humans and are responsible for up to a billion flu infections annually. The M2 protein is present in all influenza types and belongs to the class of viroporins, i.e., small proteins that form ion channels that increase membrane permeability in virus-infected cells. In influenza A and B, AM2 and BM2 are predominantly proton channels, although they also show some permeability to monovalent cations. By contrast, M2 proteins in influenza C and D, CM2 and DM2, appear to be especially selective for chloride ions, with possibly some permeability to protons. These differences point to different biological roles for M2 in types A and B versus C and D, which is also reflected in their sequences. AM2 is by far the best characterized viroporin, where mechanistic details and rationale of its acid activation, proton selectivity, unidirectionality, and relative low conductance are beginning to be understood. The present review summarizes the biochemical and structural aspects of influenza viroporins and discusses the most relevant aspects of function, inhibition, and interaction with the host.

Molecular dynamics simulations of homo-oligomeric bundles embedded within a lipid bilayer

Using molecular dynamics simulations, we studied the structure, interhelix interactions, and dynamics of transmembrane proteins. Specifically, we investigated homooligomeric helical bundle systems consisting of synthetic α-helices with either the sequence Ac-(LSLLLSL)3-NH2 (LS2) or Ac-(LSSLLSL)3-NH2 (LS3). The LS2 and LS3 helical peptides are designed to have amphipathic characteristics that form ion channels in membrane. We simulated bundles containing one to six peptides that were embedded in palmitoyl-oleoyl-phosphatidylcholine (POPC) lipid bilayer and placed between two lamellae of water. We aim to provide a fundamental understanding of how amphipathic helical peptides interact with each other and their dynamical behaviors in different homooligomeric states. To understand structural properties, we examined the helix lengths, tilt angles of individual helices and the entire bundle, interhelix distances, interhelix cross-angles, helix hydrophobic-to-hydrophilic vector projections, and the average number of interhelix hydrophilic (serine-serine) contacts lining the pore of the transmembrane channel. To analyze dynamical properties, we calculated the rotational autocorrelation function of each helix and the cross-correlation of the rotational velocity between adjacent helices. The observed structural and dynamical characteristics show that higher order bundles containing four to six peptides are composed of multiple lower order bundles of one to three peptides. For example, the LS2 channel was found to be stable in a tetrameric bundle composed of a "dimer of dimers." In addition, we observed that there is a minimum of two strong hydrophilic contacts between a pair of adjacent helices in the dimer to tetramer systems and only one strong hydrophilic interhelix contact in helix pairs of the pentamer and hexamer systems. We believe these results are general and can be applied to more complex ion channels, providing insight into ion channel stability and assembly.

Effect of cysteine mutations in the extracellular domain of CM2 on the influenza C virus replication

CM2 is the second membrane protein of influenza C virus and possesses three conserved cysteines at residue 1, 6 and 20 in its extracellular domain, all of which are involved in the formation of disulfide-linked oligomers of the molecule. In the present study, to examine the effect of CM2 oligomerization on virus replication, we generated a mutant recombinant virus, rC1620A, in which all three cysteines on CM2 were substituted to alanines. The rC1620A virus was more attenuated than the recombinant wild-type (rWT) virus in cultured cells. The CM2 protein synthesized in rC1620A-infected cells could not apparently be detected as a tetramer and was transported to the cell surface less efficiently than was authentic CM2. The amount of CM2 protein incorporated into the rC1620A virions was comparable to that into the rWT virions, although the main CM2 species in the rC1620A virions was in the form of a dimer. Analyses of one-step grown virions and virus-infected cells could not provide evidence for any difference in growth between rC1620A and rWT. On the other hand, the amount of genome present in VLPs possessing the mutant CM2 (C1620A-VLPs) was approximately 31% of that in VLPs possessing wild-type CM2 (WT-VLPs). The incoming genome from VLPs was less efficiently transported to the nucleus in the C1620A-VLP-infected cells than in WT-VLP-infected cells, leading to reduced reporter gene expression in the C1620A-VLP-infected cells. Taken together, these findings demonstrate that CM2 oligomerization affects the packaging and uncoating processes. Thus, we concluded that disulfide-linked CM2 oligomers facilitate virus growth by affecting the replication processes.

ORF8a of SARS-CoV forms an ion channel: experiments and molecular dynamics simulations

ORF8a protein is 39 residues long and contains a single transmembrane domain. The protein is synthesized using solid phase peptide synthesis and reconstituted into artificial lipid bilayers that forms cation-selective ion channels with a main conductance level of 8.9±0.8pS at elevated temperature (38.5°C). Computational modeling studies including multi nanosecond molecular dynamics simulations in a hydrated POPC lipid bilayer are done with a 22 amino acid transmembrane helix to predict a putative homooligomeric helical bundle model. A structural model of a pentameric bundle is proposed with cysteines, serines and threonines facing the pore.

Coarse-grained molecular dynamics of tetrameric transmembrane peptide bundles within a lipid bilayer

The conformations of model transmembrane peptides are studied to understand the structural and dynamical aspects of tetrameric bundles using a series of coarse grain (CG) molecular dynamics (MD) simulations since membrane proteins play a crucial role in cell function. In this work, two different amphipathic models have been constructed using similar hydrophobic/hydrophilic characteristics with two structurally distinct morphologies to evaluate the effect of roughness and hydrophilic topology on the structure of tetrameric bundles, one class that forms an ion-channel and one class that does not. Free energy calculations of typical amphipathic peptide topologies show that using a relatively smooth surface morphology allows for a stable conformation of the tetramer bundle in a diamond formation. However, the model with side chains attached to the core in order to roughen the surface has a stable square tetramer bundle which is consistent with experimental data and all-atom (AA) MD simulations. Comparisons of the CG simulations with AA MD simulations are in reasonable agreement with the formation of tetrameric homo-oligomers, partitioning within the lipid bilayer and tilt angle with respect to the bilayer normal. We concluded that a square or diamond shape tetrameric homo-oligomers could be stabilized by rational design of the peptide morphology and topology of the surface, thus allowing us to tune the permeability of the bundle or channel.

Structure and ion channel activity of the human respiratory syncytial virus (hRSV) small hydrophobic protein transmembrane domain

The small hydrophobic (SH) protein from the human respiratory syncytial virus (hRSV) is a glycoprotein of approximately 64 amino acids with one putative alpha-helical transmembrane domain. Although SH protein is important for viral infectivity, its exact role during viral infection is not clear. Herein, we have studied the secondary structure, orientation, and oligomerization of the transmembrane domain of SH (SH-TM) in the presence of lipid bilayers. Only one oligomer, a pentamer, was observed in PFO-PAGE. Using polarized attenuated total reflection-Fourier transform infrared (PATR-FTIR) spectroscopy, we show that the SH-TM is alpha-helical. The rotational orientation of SH-TM was determined by site-specific infrared dichroism (SSID) at two consecutive isotopically labeled residues. This orientation is consistent with that of an evolutionary conserved pentameric model obtained from a global search protocol using 13 homologous sequences of RSV. Conductance studies of SH-TM indicate ion channel activity, which is cation selective, and inactive below the predicted pK(a) of histidine. Thus, our results provide experimental evidence that the transmembrane domain of SH protein forms pentameric alpha-helical bundles that form cation-selective ion channels in planar lipid bilayers. We provide a model for this pore, which should be useful in mutagenesis studies to elucidate its role during the virus cycle.

Amantadine-induced conformational and dynamical changes of the influenza M2 transmembrane proton channel

The M2 protein of influenza A virus forms a transmembrane proton channel important for viral infection and replication. Amantadine blocks this channel, thus inhibiting viral replication. Elucidating the high-resolution structure of the M2 protein and its change upon amantadine binding is crucial for designing antiviral drugs to combat the growing resistance of influenza A viruses against amantadine. We used magic-angle-spinning solid-state NMR to determine the conformation and dynamics of the transmembrane domain of the protein M2TMP in the apo- and amantadine-bound states in lipid bilayers. (13)C chemical shifts and torsion angles of the protein in 1,2-dilauroyl-sn-glycero-3-phosphatidylcholine (DLPC) bilayers indicate that M2TMP is alpha-helical in both states, but the average conformation differs subtly, especially at the G34-I35 linkage and V27 side chain. In the liquid-crystalline membrane, the complexed M2TMP shows dramatically narrower lines than the apo peptide. Analysis of the homogeneous and inhomogeneous line widths indicates that the apo-M2TMP undergoes significant microsecond-time scale motion, and amantadine binding alters the motional rates, causing line-narrowing. Amantadine also reduces the conformational heterogeneity of specific residues, including the G34/I35 pair and several side chains. Finally, amantadine causes the helical segment N-terminal to G34 to increase its tilt angle by 3 degrees , and the G34-I35 torsion angles cause a kink of 5 degrees in the amantadine-bound helix. These data indicate that amantadine affects the M2 proton channel mainly by changing the distribution and exchange rates among multiple low-energy conformations and only subtly alters the average conformation and orientation. Amantadine-resistant mutations thus may arise from binding-incompetent changes in the conformational equilibrium.

Side-chain conformation of the M2 transmembrane peptide proton channel of influenza a virus from 19F solid-state NMR

The M2 transmembrane peptide (M2TMP) of the influenza A virus forms a tetrameric helical bundle that acts as a proton-selective channel important in the viral life cycle. The side-chain conformation of the peptide is largely unknown and is important for elucidating the proton-conducting mechanism and the channel stability. Using a 19F spin diffusion NMR technique called CODEX, we have measured the oligomeric states and interhelical side chain-side chain 19F-19F distances at several residues using singly fluorinated M2TMP bound to DMPC bilayers. 19F CODEX data at a key residue of the proton channel, Trp41, confirm the tetrameric state of the peptide and yield a nearest-neighbor interhelical distance of approximately 11 A under both neutral and acidic pH. Since the helix orientation is precisely known from previous 15N NMR experiments and the backbone channel diameter has a narrow allowed range, this 19F distance constrains the Trp41 side-chain conformation to t90 (chi1 approximately 180 degrees , chi2 approximately 90 degrees ). This Trp41 rotamer, combined with a previously measured 15N-13C distance between His37 and Trp411, suggests that the His37 rotamer is t-160. The implication of the proposed (His37, Trp41) rotamers to the gating mechanism of the M2 proton channel is discussed. Binding of the antiviral drug amantadine to the peptide does not affect the F-F distance at Trp41. Interhelical 19F-19F distances are also measured at residues 27 and 38, each mutated to 4-19F-Phe. For V27F-M2TMP, the 19F-19F distances suggest a mixture of dimers and tetramers, whereas the L38F-M2TMP data indicate two tetramers of different sizes, suggesting side chain conformational heterogeneity at this lipid-facing residue. This work shows that 19F spin diffusion NMR is a valuable tool for determining long-range intermolecular distances that shed light on the mechanism of action and conformational heterogeneity of membrane protein oligomers.

Molecular dynamics simulation approach for the prediction of transmembrane helix–helix heterodimers assembly

Computational methods are useful to identify favorable structures of transmembrane (TM) helix oligomers when experimental data are not available or when they cannot help to interpret helix-helix association. We report here a global search method using molecular dynamics (MD) simulations to predict the structures of transmembrane homo and heterodimers. The present approach is based only on sequence information without any experimental data and is first applied to glycophorin A to validate the protocol and to the HER2-HER3 heterodimer receptor. The method successfully reproduces the experimental structures of the TM domain of glycophorin A (GpA(TM)) with a root mean square deviation of 1.5 A. The search protocol identifies three energetically stable models of the TM domain of HER2-HER3 receptor with favorable helix-helix arrangement, including right-handed and left-handed coiled-coils. The predicted TM structures exhibit the GxxxG-like motif at the dimer interface which is presumed to drive receptor oligomerization. We demonstrate that native structures of TM domain can be predicted without quantitative experimental data. This search protocol could help to predict structures of the TM domain of HER heterodimer family.

Mutagenesis data in the automated prediction of transmembrane helix dimers

We present a molecular modeling protocol that selects modeled protein structures based on experimental mutagenesis results. The computed effect of a point mutation should be consistent with its experimental effect for correct models; mutations that do not affect protein stability and function should not affect the computed energy of a correct model while destabilizing mutations should have unfavorable computed energies. On the other hand, an incorrect model will likely display computed energies that are inconsistent with experimental results. We added terms to our energy function which penalize models that are inconsistent with experimental results. This creates a selective advantage for models that are consistent with experimental results in the Monte Carlo simulated annealing protocol we use to search conformational space. We calibrated our protocol to predict the structure of transmembrane helix dimers using glycophorin A as a model system. Inclusion of mutational data in this protocol compensates for the limitations of our force field and the limitations of our conformational search. We demonstrate an application of this structure prediction protocol by modeling the transmembrane region of the BNIP3 apoptosis factor.

Computational design of heterochiral peptides against a helical target

Polypeptides incorporating D-amino acids occasionally occur in nature and are an important class of pharmaceutical molecules. With the use of heterochiral Monte Carlo (HCMC), a method inspired by the de novo design of proteins, we develop peptide scaffolds for interacting with a molecular target, a left-handed alpha-helix. The HCMC approach concurrently seeks to optimize a peptide sequence, its internal conformation, and its docked conformation with a target surface. Several major classes of interactions are observed: (1) homochiral interactions between two alphaL helices, (2) heterochiral interactions between an alphaL and an alphaR helix, and (3) heterochiral interactions between the alphaL target and novel nonhelical structures. We explore the application of HCMC to simulating the preferential enantioselectivity of heterochiral complexes. Implications for biomimetic design in molecular recognition are discussed.

Automated use of mutagenesis data in structure prediction

In the absence of experimental structural determination, numerous methods are available to indirectly predict or probe the structure of a target molecule. Genetic modification of a protein sequence is a powerful tool for identifying key residues involved in binding reactions or protein stability. Mutagenesis data is usually incorporated into the modeling process either through manual inspection of model compatibility with empirical data, or through the generation of geometric constraints linking sensitive residues to a binding interface. We present an approach derived from statistical studies of lattice models for introducing mutation information directly into the fitness score. The approach takes into account the phenotype of mutation (neutral or disruptive) and calculates the energy for a given structure over an ensemble of sequences. The structure prediction procedure searches for the optimal conformation where neutral sequences either have no impact or improve stability and disruptive sequences reduce stability relative to wild type. We examine three types of sequence ensembles: information from saturation mutagenesis, scanning mutagenesis, and homologous proteins. Incorporating multiple sequences into a statistical ensemble serves to energetically separate the native state and misfolded structures. As a result, the prediction of structure with a poor force field is sufficiently enhanced by mutational information to improve accuracy. Furthermore, by separating misfolded conformations from the target score, the ensemble energy serves to speed up conformational search algorithms such as Monte Carlo-based methods.

A systematic search method for the identification of tightly packed transmembrane parallel alpha-helices

Membrane proteins play a major role in number of biological processes such as signaling pathways. The determination of the three-dimensional structure of these proteins is increasingly important for our understanding of their structure-function relationships. Due to the difficulty in isolating membrane proteins for X-ray diffraction studies, computational techniques are being developed to generate the 3D structures of TM domains. Here, we present a systematic search method for the identification of energetically favorable and tightly packed transmembrane parallel alpha-helices. The first step in our systematic search method is the generation of 3D models for pairs of parallel helix bundles with all possible orientations followed by an energy-based filter to eliminate structures with severe non-bonded contacts. Then, a RMS-based filter was used to cluster these structures into families. Furthermore, these dimers were energy minimized using molecular mechanics force field. Finally, we identified the tightly packed parallel alpha-helices by using an interface surface area. To validate our search method, we compared our predicted GlycophorinA dimer structures with the reported NMR structures. With our search method, we are able to reproduce NMR structures of GPA with 0.9A RMSD. In addition, by considering the reported mutational data on GxxxG motif interactions, twenty percent of our predicted dimers are within in the 2.0A RMSD. The dimers obtained from our method were used to generate parallel trimeric and tetramer TM structures of GPA and found that the structure of GPA might exist only in a dimer form as reported earlier.

The transmembrane oligomers of coronavirus protein E

We have tested the hypothesis that severe acute respiratory syndrome (SARS) coronavirus protein E (SCoVE) and its homologs in other coronaviruses associate through their putative transmembrane domain to form homooligomeric α-helical bundles in vivo. For this purpose, we have analyzed the results of molecular dynamics simulations where all possible conformational and aggregational space was systematically explored. Two main assumptions were considered; the first is that protein E contains one transmembrane α-helical domain, with its N- and C-termini located in opposite faces of the lipid bilayer. The second is that protein E forms the same type of transmembrane oligomer and with identical backbone structure in different coronaviruses. The models arising from the molecular dynamics simulations were tested for evolutionary conservation using 13 coronavirus protein E homologous sequences. It is extremely unlikely that if any of our assumptions were not correct we would find a persistent structure for all the sequences tested. We show that a low energy dimeric, trimeric and two pentameric models appear to be conserved through evolution, and are therefore likely to be present in vivo. In support of this, we have observed only dimeric, trimeric, and pentameric aggregates for the synthetic transmembrane domain of SARS protein E in SDS. The models obtained point to residues essential for protein E oligomerization in the life cycle of the SARS virus, specifically N15. In addition, these results strongly support a general model where transmembrane domains transiently adopt many aggregation states necessary for function.

Lipid bilayer perturbations around a transmembrane nanotube: a coarse grain molecular dynamics study

The perturbations induced in a lipid bilayer by the presence of a transmembrane nanotube are investigated using coarse grained molecular dynamics. Meniscus formation by the lipids and tilting of the nanotube occur in response to hydrophobic mismatch, although these two effects do not compensate completely for the total mismatch. The lipid head-to-tail vector field is examined and shows strong ordering in the membrane plane regardless of the nanotube length. Molecular layering at the lipid-nanotube interface is reported. This study extends previous theoretical approaches to a more realistic setting.

A structural model of EmrE, a multi-drug transporter from Escherichia coli

Using a recently reported computational method, we describe an approach to model the structure of EmrE, a proton coupled multi-drug transporter of Escherichia coli. EmrE is the smallest ion-coupled transporter known; it functions as an oligomer and each monomer comprises four transmembrane segments. Because of its size, EmrE provides a unique experimental paradigm. The computational method does not afford a unique solution for the monomer. The experimental constraints available were used to select the most likely structure and to dock two monomers together to yield a dimer. The model is further validated by modeling of Hsmr, an EmrE homolog with a remarkable amino acid composition with over 40% of Ala and Val. The Hsmr model is similar to that of EmrE, with the majority of the Ala or Val residues facing the lipid. In addition, the model of EmrE features a putative substrate-binding site very similar to that observed in BmrR, a transcription activator of multi-drug transporters, with a similar substrate profile. The two crucial residues that couple proton fluxes with substrate binding in the homo-dimer of EmrE, Glu-14, have a spatial arrangement that agrees with proposed molecular mechanisms of transport.

Quantification of helix-helix binding affinities in micelles and lipid bilayers

A theoretical approach for estimating association free energies of alpha-helices in nonpolar media has been developed. The parameters of energy functions have been derived from DeltaDeltaG values of mutants in water-soluble proteins and partitioning of organic solutes between water and nonpolar solvents. The proposed approach was verified successfully against three sets of published data: (1) dissociation constants of alpha-helical oligomers formed by 27 hydrophobic peptides; (2) stabilities of 22 bacteriorhodopsin mutants, and (3) protein-ligand binding affinities in aqueous solution. It has been found that coalescence of helices is driven exclusively by van der Waals interactions and H-bonds, whereas the principal destabilizing contributions are represented by side-chain conformational entropy and transfer energy of atoms from a detergent or lipid to the protein interior. Electrostatic interactions of alpha-helices were relatively weak but important for reproducing the experimental data. Immobilization free energy, which originates from restricting rotational and translational rigid-body movements of molecules during their association, was found to be less than 1 kcal/mole. The energetics of amino acid substitutions in bacteriorhodopsin was complicated by specific binding of lipid and water molecules to cavities created in certain mutants.

Three-dimensional structure of the channel-forming trans-membrane domain of virus protein "u" (Vpu) from HIV-1

The three-dimensional structure of the channel-forming trans-membrane domain of virus protein "u" (Vpu) of HIV-1 was determined by NMR spectroscopy in micelle and bilayer samples. Vpu(2-30+) is a 36-residue polypeptide that consists of residues 2-30 from the N terminus of Vpu and a six-residue "solubility tag" at its C terminus that facilitates the isolation, purification, and sample preparation of this highly hydrophobic minimal channel-forming domain. Nearly all of the resonances in the two-dimensional 1H/15N HSQC spectrum of uniformly 15N labeled Vpu(2-30+) in micelles are superimposable on those from the corresponding residues in the spectrum of full-length Vpu, which indicates that the structure of the trans-membrane domain is not strongly affected by the presence of the cytoplasmic domain at its C terminus. The two-dimensional 1H/15N PISEMA spectrum of Vpu(2-30+) in lipid bilayers aligned between glass plates has been fully resolved and assigned. The "wheel-like" pattern of resonances in the spectrum is characteristic of a slightly tilted membrane-spanning helix. Experiments were also performed on weakly aligned micelle samples to measure residual dipolar couplings and chemical shift anisotropies. The analysis of the PISA wheels and Dipolar Waves obtained from both weakly and completely aligned samples show that Vpu(2-30+) has a trans-membrane alpha-helix spanning residues 8-25 with an average tilt of 13 degrees. The helix is kinked slightly at Ile17, which results in tilts of 12 degrees for residues 8-16 and 15 degrees for residues 17-25. A structural fit to the experimental solid-state NMR data results in a three-dimensional structure with precision equivalent to an RMSD of 0.4 A. Vpu(2-30+) exists mainly as an oligomer on PFO-PAGE and forms ion-channels, a most frequent conductance of 96(+/- 6) pS in lipid bilayers. The structural features of the trans-membrane domain are determinants of the ion-channel activity that may be associated with the protein's role in facilitating the budding of new virus particles from infected cells.

The closed state of a H+ channel helical bundle combining precise orientational and distance restraints from solid state NMR

An interhelical distance has been precisely measured by REDOR solid-state NMR spectroscopy in the transmembrane tetrameric bundle of M2-TMP, from the M2 proton channel of the influenza A viral coat. The high-resolution structure of the helical backbone has been determined using orientational restraints from uniformly aligned peptide preparations in hydrated dimyristoylphosphatidylcholine bilayers. Here, the distance between (15)N(pi) labeled His37 and (13)C(gamma) labeled Trp41 is determined to be less than 3.9 A. Such a short distance, in combination with the known tilt and rotational orientation of the individual helices, permits not only a determination of which specific side chain pairings give rise to the interaction, but also the side chain torsion angles and restraints for the tetrameric bundle can also be characterized. The resulting proton channel structure is validated in a variety of ways. Both histidine and tryptophan side chains are oriented in toward the pore where they can play a significant functional role. The channel appears to be closed by the proximity of the four indoles consistent with electrophysiology and mutagenesis studies of the intact protein at pH 7.0 and above. The pore maintains its integrity to the N terminal side of the membrane, and at the same time, a cavity is generated that appears adequate for binding amantadine. Finally, the observation of a 2 kHz coupling in the PISEMA spectrum of (15)N(pi)His37 validates the orientation of the His37 side chain based on the observed REDOR distance.

Structural aspects of oligomerization taking place between the transmembrane alpha-helices of bitopic membrane proteins

Recent advances in biophysical methods have been able to shed more light on the structures of helical bundles formed by the transmembrane segments of bitopic membrane proteins. In this manuscript, I attempt to review the biological importance and diversity of these interactions, the energetics of bundle formation, motifs capable of inducing oligomerization and methods capable of detecting, solving and predicting the structures of these oligomeric bundles. Finally, the structures of the best characterized instances of transmembrane alpha-helical bundles formed by bitopic membrane proteins are described in detail.

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.

Expression, purification, and activities of full-length and truncated versions of the integral membrane protein Vpu from HIV-1

Vpu is an 81-residue accessory protein of HIV-1. Because it is a membrane protein, it presents substantial technical challenges for the characterization of its structure and function, which are of considerable interest because the protein enhances the release of new virus particles from cells infected with HIV-1 and induces the intracellular degradation of the CD4 receptor protein. The Vpu-mediated enhancement of the virus release rate from HIV-1-infected cells is correlated with the expression of an ion channel activity associated with the transmembrane hydrophobic helical domain. Vpu-induced CD4 degradation and, to a lesser extent, enhancement of particle release are both dependent on the phosphorylation of two highly conserved serine residues in the cytoplasmic domain of Vpu. To define the minimal folding units of Vpu and to identify their activities, we prepared three truncated forms of Vpu and compared their structural and functional properties to those of full-length Vpu (residues 2-81). Vpu(2-37) encompasses the N-terminal transmembrane alpha-helix; Vpu(2-51) spans the N-terminal transmembrane helix and the first cytoplasmic alpha-helix; Vpu(28-81) includes the entire cytoplasmic domain containing the two C-terminal amphipathic alpha-helices without the transmembrane helix. Uniformly isotopically labeled samples of the polypeptides derived from Vpu were prepared by expression of fusion proteins in E. coli and were studied in the model membrane environments of lipid micelles by solution NMR spectroscopy and oriented lipid bilayers by solid-state NMR spectroscopy. The assignment of backbone resonances enabled the secondary structure of the constructs corresponding to the transmembrane and the cytoplasmic domains of Vpu to be defined in micelle samples by solution NMR spectroscopy. Solid-state NMR spectra of the polypeptides in oriented lipid bilayers demonstrated that the topology of the domains is retained in the truncated polypeptides. The biological activities of the constructs of Vpu were evaluated. The ion channel activity is confined to the transmembrane alpha-helix. The C-terminal alpha-helices modulate or promote the oligomerization of Vpu in the membrane and stabilize the conductive state of the channel, in addition to their involvement in CD4 degradation.