Statistical analysis of predicted transmembrane alpha-helices.

Statistical analyses were undertaken for putative transmembrane alpha-helices obtained from a database representing the subset of membrane proteins available in Swiss-Prot. The average length of a transmembrane alpha-helix was found to be 22-21 amino acids with a large variation around the mean. The transfer free energy from water to oil of a transmembrane alpha-helix in bitopic proteins, -48 kcal/mol, is higher than that in polytopic proteins, -39 kcal/mol, and is nearly identical to that obtained by assuming a random distribution of solely hydrophobic amino acids in the alpha-helix. The amino acid composition of hydrophobic residues is similar in bitopic and polytopic proteins. In contrast, the more polar the amino acids are, the less likely they are to be found in bitopic proteins compared to polytopic ones. This most likely reflects the ability of alpha-helical bundles to shield the polarity of residues from the hydrophobic bilayer. One half of all amino acids were distributed nonrandomly in both bitopic and polytopic proteins. A preference was found for tyrosine and tryptophan residues to be at the ends of transmembrane alpha-helices. Correlated distribution analysis of amino acid pairs indicated that most amino acids are independently distributed in each helix. Exceptions are cysteine, tyrosine, and tryptophan which appear to cluster closely to one another and glycines which are preferentially found on the same side of alpha-helices.

The structure of the HIV-1 Vpu ion channel: modelling and simulation studies.

Vpu is an 81 amino acid auxiliary protein in HIV-1 which exhibits channel activity. We used two homo-pentameric bundles with the helical transmembrane segments derived from FTIR spectroscopy in combination with a global molecular dynamics search protocol: (i) tryptophans (W) pointing into the pore, and (ii) W facing the lipids. Two equivalent bundles have been generated using a simulated annealing via a restrained molecular dynamics simulations (SA/MD) protocol. A fifth model was generated via SA/MD with all serines facing the pore. The latter model adopts a very stable structure during the 2 ns of simulation. The stability of the models with W facing the pore depends on the starting structure. A possible gating mechanism is outlined.

Helicity, membrane incorporation, orientation and thermal stability of the large conductance mechanosensitive ion channel from E. coli.

In this report, we present structural studies on the large conductance mechanosensitive ion channel (MscL) from E. coli in detergent micelles and lipid vesicles. Both transmission Fourier transform infrared spectroscopy and circular dichroism (CD) spectra indicate that the protein is highly helical in detergents as well as liposomes. The secondary structure of the proteins was shown to be highly resistant towards denaturation (25-95 degrees C) based on an ellipticity thermal profile. Amide H+/D+ exchange was shown to be extensive (ca. 66%), implying that two thirds of the protein are water accessible. MscL, reconstituted in oriented lipid bilayers, was shown to possess a net bilayer orientation using dichroic ratios measured by attenuated total-reflection Fourier transform infrared spectroscopy. Here, we present and discuss this initial set of structural data on this new family of ion-channel proteins.

Use of a new label, (13)==(18)O, in the determination of a structural model of phospholamban in a lipid bilayer. Spatial restraints resolve the ambiguity arising from interpretations of mutagenesis data.

A structural model of pentameric phospholamban (Plb) in a lipid bilayer has been derived using a combination of experimental data, obtained from ATR-FTIR site-directed dichroism, and the implementation of the resulting restraints during a molecular dynamics simulation. Plb (residues 24-52) has been synthesised incorporating a new label, 1-(13)C==(18)O, at residues 42 and 43. We have not only determined the tilt of the helices, 10(+/-6) degrees, but also the relative orientation of the transmembrane segments, with an omega angle of -32(+/-10) degrees for L42. This angle is taken as zero in the direction of the helix tilt. Plb is a simple test case where site-directed dichroism has been applied to resolve the indeterminacy arising from the mutagenesis data available. The results presented point specifically to a single structural model for Plb.

Structural model of the phospholamban ion channel complex in phospholipid membranes.

Phospholamban is a 52 amino acid residue membrane protein involved with the regulation of calcium levels across sarcoplasmic reticulum membranes in cardiac muscle cells. The N-terminal 30 amino acid residues of the protein are largely hydrophilic and include two sites whose phosphorylation is thought to dissociate an inhibitory complex between phospholamban and Ca2+ ATPase. The C-terminal 22 amino acid residues are largely hydrophobic, anchor the protein in the membrane and are responsible for Ca2+ selective ion conductance. Specific interactions between the transmembrane domains stabilize a pentameric protein complex. We have obtained circular dichroism (CD), transmission Fourier transform infrared (FTIR) and attenuated total reflection Fourier transform infrared (ATR-FTIR) spectra of the full-length protein and have compared these results to those from a 28 residue peptide that includes the transmembrane domain. Both proteins reconstituted into phospholipid membranes are largely alpha-helical by CD and FTIR. Polarized ATR-FTIR measurements show that both the cytosolic and transmembrane helices are oriented perpendicular to the membrane plane with a tilt of 28 (+/- 6) degrees with respect to the membrane normal. This tilt angle is in close agreement to that calculated from a model for the transmembrane domain of phospholamban suggested by mutagenesis and molecular modeling. Phosphorylation does not significantly change the secondary structure or orientation of the protein. The pentameric complex is modeled as a left-handed coiled-coil of five long helices (40 (+/- 3) residues) that extend across the membrane from the lumenal carboxy terminus to the phosphorylation site in the cytoplasm. The helix bundle forms a perpendicular ion pore that may begin at a distance (17 to 29 A) from the membrane surface. Based on the above, we propose a mechanism by which phospholamban regulates Ca2+ levels across membranes that takes into account both its selective ion conductance and inhibitory association with the Ca2+ pump.

SARS E protein in phospholipid bilayers: an anomalous X-ray reflectivity study.

We report on an anomalous X-ray reflectivity study to locate a labelled residue of a membrane protein with respect to the lipid bilayer. From such experiments, important constraints on the protein or peptide conformation can be derived. Specifically, our aim is to localize an iodine-labelled phenylalanine in the SARS E protein, incorporated in DMPC phospholipid bilayers, which are deposited in the form of thick multilamellar stacks on silicon surfaces. Here, we discuss the experimental aspects and the difficulties associated with the Fourier synthesis analysis that gives the electron density profile of the membranes.

Substitution rates in alpha-helical transmembrane proteins.

It has been shown previously that some membrane proteins have a conserved core of amino acid residues. This idea not only serves to orient helices during model building exercises but may also provide insight into the structural role of residues mediating helix-helix interactions. Using experimentally determined high-resolution structures of alpha-helical transmembrane proteins we show that, of the residues within the hydrophobic transmembrane spans, the residues at lipid and subunit interfaces are more evolutionarily variable than those within the lipid-inaccessible core of a polypeptide's transmembrane domain. This supports the idea that helix-helix interactions within the same polypeptide chain and those at the interface between different polypeptide chains may arise in distinct ways. To show this, we use a new method to estimate the substitution rate of an amino acid residue given an alignment and phylogenetic tree of closely related proteins. This method gives better sensitivity in the otherwise-conserved transmembrane domains than a conventional similarity analysis and is relatively insensitive to the sequences used.

The effect of nucleotide bias upon the composition and prediction of transmembrane helices.

Transmembrane helices are the most readily predictable secondary structure components of proteins. They can be predicted to a high degree of accuracy in a variety of ways. Many of these methods compare new sequence data with the sequence characteristics of known transmembrane domains. However, the known transmembrane sequences are not necessarily representative of a particular organism. We attempt to demonstrate that parameters optimized for the known transmembrane domains are far from optimal when predicting transmembrane regions in a given genome. In particular, we have tested the effect of nucleotide bias upon the composition and hence the prediction characteristics of transmembrane helices. Our analysis shows that nucleotide bias of a genome has a strong and predictable influence upon the occurrences of several of the most important hydrophobic amino acids found within transmembrane helices. Thus, we show that nucleotide bias should be taken into account when determining putative transmembrane domains from sequence data.

The orientation of the antibiotic peptide maculatin 1.1 in DMPG and DMPC lipid bilayers. Support for a pore-forming mechanism.

Maculatin 1.1 is an antimicrobial peptide isolated from the Australian tree frog Litoria genimaculata that adopts an amphipathic, alpha-helical structure in solution. Its orientation and conformation when incorporated to pre-formed DMPG (1,2-dimyristoyl-sn-glycero-3-phosphoglycerol) and DMPC (1,2-dimyristoyl-sn-glycero-3-phosphocholine) vesicles was determined using polarised Fourier transform infrared-attenuated total reflection infrared and deuterium exchange experiments. For DMPG membranes, our results show insertion of 70% of the maculatin 1.1 molecules, with an angle of insertion of approximately 35 degrees to the membrane normal and with a predominant alpha-helical structure. These results suggest that maculatin 1.1 acts through a pore-forming mechanism to lyse bacterial membranes. A similar degree of insertion in DMPG (65%) and alpha-helical structure was observed for a biologically inactive, less amphipathic maculatin 1.1 analogue, P15A, although the helix tilt was found to be greater (46 degrees) than for maculatin 1.1. Similar experiments performed using DMPC liposomes showed poor insertion, less than 5%, for both maculatin 1.1 and its analogue. In addition, the shape of the amide I band in these samples is consistent with alpha-helix, beta-structure and disordered structures being present in similar proportion. These results clearly show that maculatin 1.1 inserts preferentially in negatively charged membranes (DMPG) which mimic the negatively charged membrane of Gram-positive bacteria. We attribute the high percentage of insertion of the biologically inactive analogue in DMPG to the fact that its concentration on the membrane surface in our experiments is likely to be much higher than that found in physiological conditions.

Viral ion channel proteins in model membranes: a comparative study by X-ray reflectivity.

We have investigated the effect of the transmembrane domain of three viral ion channel proteins on the lipid bilayer structure by X-ray reflectivity and scattering from oriented planar bilayers. The proteins show a similar effect on the lipid bilayer structural parameters: an increase in the lipid bilayer hydrophobic core, a decrease in the amplitude of the vertical density profile and a systematic change in the ordering of the acyl chains as a function of protein-to-lipid ratio. These results are discussed in a comparative view.

SARS coronavirus E protein in phospholipid bilayers: an x-ray study.

We investigated the structure of the hydrophobic domain of the severe acute respiratory syndrome E protein in model lipid membranes by x-ray reflectivity and x-ray scattering. In particular, we used x-ray reflectivity to study the location of an iodine-labeled residue within the lipid bilayer. The label imposes spatial constraints on the protein topology. Experimental data taken as a function of protein/lipid ratio P/L and different swelling states support the hairpin conformation of severe acute respiratory syndrome E protein reported previously. Changes in the bilayer thickness and acyl-chain ordering are presented as a function of P/L, and discussed in view of different structural models.

Structure of the influenza C virus CM2 protein transmembrane domain obtained by site-specific infrared dichroism and global molecular dynamics searching.

The 115-residue protein CM2 from Influenza C virus has been recently characterized as a tetrameric integral membrane glycoprotein. Infrared spectroscopy and site-directed infrared dichroism were utilized here to determine its transmembrane structure. The transmembrane domain of CM2 is alpha-helical, and the helices are tilted by beta = (14.6 +/- 3.0) degrees from the membrane normal. The rotational pitch angle about the helix axis omega for the 1-(13)C-labeled residues Gly(59) and Leu(66) is omega = (218 +/- 17) degrees, where omega is defined as zero for a residue pointing in the direction of the helix tilt. A detailed structure was obtained from a global molecular dynamics search utilizing the orientational data as an energy refinement term. The structure consists of a left-handed coiled-coil with a helix crossing angle of Omega = 16 degrees. The putative transmembrane pore is occluded by the residue Met(65). In addition hydrogen/deuterium exchange experiments show that the core is not accessible to water.

Recursive use of evolutionary conservation data in molecular modeling of membrane proteins A model of the multidrug H+ antiporter emrE.

Membrane proteins are currently the most biomedically important family of proteins, serving as targets for the majority of pharmaceutical agents. It is also clear that they are invariably abundant in all of the genomes sequence so far, representing up to a third of all open reading frames. Finally, and regrettably, it is clear that they are highly resistant to structural elucidation, representing less than 0.2% of the Protein Data Bank. Recent accomplishments in genome sequencing efforts, however, may help offset this imbalance through the availability of evolutionary conservation data. Herein, we develop a novel approach, utilizing a combination of evolutionary conservation data and global searching molecular dynamics simulations to model membrane proteins, deriving a model for the multidrug H+ antiporter EmrE, a transmembrane four-helix bundle. Structures resulting from an extensive, rotational molecular dynamics search, were evaluated by comparing the residue specific interaction energy and the evolutionary conservation data. Subsequent rounds of molecular dynamics, in which confinement of the search space was undertaken in order to achieve a self consistent result, point to a structure that best satisfies the evolutionary conservation data. As the conservation patterns calculated for each of the helices suggested that the different conservation pattern for helix 3 (as well as being the most conserved) might be due to the oligomeric nature of EmrE, a dodecamer of helices was constructed based on the result of a search of helix 3 as a trimer. The resulting interaction energy per residue in the final model is in reasonable agreement with the evolutionary data and consistent with recent site directed mutagenesis experiments, pointing to the strength of this method as a general tool.

Do more complex organisms have a greater proportion of membrane proteins in their genomes?

One may speculate that higher organisms require a proportionately greater abundance of membrane proteins within their genomes in order to furnish the requirements of differentiated cell types, compartmentalization, and intercellular signalling. With the recent availability of several complete prokaryotic genome sequences and sufficient progress in many eukaryotic genome sequencing projects, we seek to test this hypothesis. Using optimized hydropathy analysis of proteins in several, diverse proteomes, we show that organisms of the three domains of life-Eukarya, Eubacteria, and Archaea-have similar proportions of alpha-helical membrane proteins within their genomes and that these are matched by the complexity of the aqueous components.

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.

Are membrane proteins "inside-out" proteins?

One of the central paradigms of structural biology is that membrane proteins are "inside-out" proteins, in that they have a core of polar residues surrounded by apolar residues. This is the reverse of the characteristics found in water-soluble proteins. We have decided to test this paradigm, now that sufficient numbers of transmembrane alpha-helical structures are accessible to statistical analysis. We have analyzed the correlation between accessibility and hydrophobicity of both individual residues and complete helices. Our analyses reveal that hydrophobicity of residues in a transmembrane helical bundle does not correlate with any preferred location and that the hydrophilic vector of a helix is a poor indicator of the solvent exposed face of a helix. Neither polar nor hydrophobic residues show any bias for the exterior or the interior of a transmembrane domain. As a control, analysis of water-soluble helical bundles performed in a similar manner has yielded clear correlations between hydrophobicity and accessibility. We therefore conclude that, based on the data set used, membrane proteins as "inside-out" proteins is an unfounded notion, suggesting that packing of alpha-helices in membranes is better understood by maximization of van der Waal's forces, rather than by a general segregation of hydrophobicities driven by lipid exclusion.

Mapping the energy surface of transmembrane helix-helix interactions.

Transmembrane helices are no longer believed to be just hydrophobic segments that exist solely to anchor proteins to a lipid bilayer, but rather they appear to have the capacity to specify function and structure. Specific interactions take place between hydrophobic segments within the lipid bilayer whereby subtle mutations that normally would be considered innocuous can result in dramatic structural differences. That such specificity takes place within the lipid bilayer implies that it may be possible to identify the most favorable interaction surface of transmembrane alpha-helices based on computational methods alone, as shown in this study. Herein, an attempt is made to map the energy surface of several transmembrane helix-helix interactions for several homo-oligomerizing proteins, where experimental data regarding their structure exist (glycophorin A, phospholamban, Influenza virus A M2, Influenza virus C CM2, and HIV vpu). It is shown that due to symmetry constraints in homo-oligomers the computational problem can be simplified. The results obtained are mostly consistent with known structural data and may additionally provide a view of possible alternate and intermediate configurations.

Use of a single glycine residue to determine the tilt and orientation of a transmembrane helix. A new structural label for infrared spectroscopy.

Site-directed dichroism is an emerging technique for the determination of membrane protein structure. However, due to a number of factors, among which is the high natural abundance of (13)C, the use of this technique has been restricted to the study of small peptides. We have overcome these problems through the use of a double C-deuterated glycine as a label. The modification of a single residue (Gly) in the transmembrane segment of M2, a protein from the Influenza A virus that forms H(+)-selective ion channels, has allowed us to determine its helix tilt and rotational orientation. Double C-deuteration shifts the antisymmetric and symmetric stretching vibrations of the CD(2) group in glycine to a transparent region of the infrared spectrum where the dichroic ratio of these bands can be measured. The two dichroisms, along with the helix amide I dichroic ratio, have been used to determine the helix tilt and rotational orientation of M2. The results are entirely consistent with previous site-directed dichroism and solid-state NMR experiments, validating C-deuterated glycine (GlyCD(2)) as a structural probe that can now be used in the study of polytopic membrane proteins.

A new method to model membrane protein structure based on silent amino acid substitutions.

The importance of accurately modeling membrane proteins cannot be overstated, in lieu of the difficulties in solving their structures experimentally. Often, however, modeling procedures (e.g., global searching molecular dynamics) generate several possible candidates rather then pointing to a single model. Herein we present a new approach to select among candidate models based on the general hypothesis that silent amino acid substitutions, present in variants identified from evolutionary conservation data or mutagenesis analysis, do not affect the stability of a native structure but may destabilize the non-native structures also found. The proof of this hypothesis has been tested on the alpha-helical transmembrane domains of two homodimers, human glycophorin A and human CD3-zeta, a component of the T-cell receptor. For both proteins, only one structure was identified using all the variants. For glycophorin A, this structure is virtually identical to the structure determined experimentally by NMR. We present a model for the transmembrane domain of CD3-zeta that is consistent with predictions based on mutagenesis, homology modeling, and the presence of a disulfide bond. Our experiments suggest that this method allows the prediction of transmembrane domain structure based only on widely available evolutionary conservation data.

Exploring models of the influenza A M2 channel: MD simulations in a phospholipid bilayer.

The M2 protein of influenza A virus forms homotetrameric helix bundles, which function as proton-selective channels. The native form of the protein is 97 residues long, although peptides representing the transmembrane section display ion channel activity, which (like the native channel) is blocked by the antiviral drug amantadine. As a small ion channel, M2 may provide useful insights into more complex channel systems. Models of tetrameric bundles of helices containing either 18 or 22 residues have been simulated while embedded in a fully hydrated 1-palmitoyl-2-oleoyl-sn-glycerol-3-phosphatidylcholine bilayer. Several different starting models have been used. These suggest that the simulation results, at least on a nanosecond time scale, are sensitive to the exact starting structure. Electrostatics calculations carried out on a ring of four ionizable aspartate residues at the N-terminal mouth of the channel suggest that at any one time, only one will be in a charged state. Helix bundle models were mostly stable over the duration of the simulation, and their helices remained tilted relative to the bilayer normal. The M2 helix bundles form closed channels that undergo breathing motions, alternating between a tetramer and a dimer-of-dimers structure. Under these conditions either the channel forms a pocket of trapped waters or it contains a column of waters broken predominantly at the C-terminal mouth of the pore. These waters exhibit restricted motion in the pore and are effectively "frozen" in a way similar to those seen in previous simulations of a proton channel formed by a four-helix bundle of a synthetic leucine-serine peptide (, Biophys. J. 77:2400-2410).