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

A conserved tetrameric interaction of cry toxin helix α3 suggests a functional role for toxin oligomerization

Crystal (Cry) toxins are widely used for insect control, but their mechanism of toxicity is still uncertain. These toxins can form lytic pores in vitro, and water soluble tetrameric pre-pore intermediates have been reported. Even the precise oligomeric state of the toxin in membranes, trimeric or tetrameric, is still a debated issue. Based on previous reports, we have assumed that interactions between toxin monomers in solution are at least partly mediated by domain I, and we have analyzed in silico the homo-oligomerization tendencies of the domain I α-helices individually. Using many homologous sequences for each α-helix, our strategy allows selection of evolutionarily conserved interactions. These interactions appeared only in helices α3 and α5, but only α3 produced a suitably oriented or α-helical sample in lipid bilayers, forming homotetramers in C14-betaine, and allowing determination of its rotational orientation in lipid bilayers using site-specific infrared dichroism (SSID). The determined orientation in the tetrameric model is in agreement with only one of the evolutionarily conserved models. In addition mutation R99E, which was found to inhibit oligomerization experimentally, greatly destabilized the tetramer in molecular dynamic simulations. In this model, helix 3 is able to form inter-monomer interactions without significant rearrangements of domain I, which is compatible with the available crystal structure of Cry toxins in solution. The model presented here at least partially explains the reported tetrameric oligomerization of Cry toxins in solution and the inhibition of this oligomerization by a synthetic α3 peptide.

STRUCTURE PREDICTION OF HELICAL TRANSMEMBRANE PROTEINS AT TWO LENGTH SCALES

As the first step toward a multi-scale, hierarchical computational approach for membrane protein structure prediction, the packing of transmembrane helices was modeled at the residue and atom levels, respectively. For predictions at the residue level, the helix-helix and helix-membrane interactions were described by a set of knowledge-based energy functions. For predictions at the atom level, CHARMM19 force field was used. To facilitate the system to overcome energy barriers, the Wang–Landau method was employed, where a random walk is performed in the energy space with a uniform probability. Native-like structures were predicted at both levels for two model systems, each of which consists of two transmembrane helices. Interestingly, consistent results were obtained from simulations at the residue and atom levels for the same system, strongly suggesting the feasibility of a hierarchical approach for membrane protein structure predictions.

Structure and inhibition of the SARS coronavirus envelope protein ion channel

The envelope (E) protein from coronaviruses is a small polypeptide that contains at least one α-helical transmembrane domain. Absence, or inactivation, of E protein results in attenuated viruses, due to alterations in either virion morphology or tropism. Apart from its morphogenetic properties, protein E has been reported to have membrane permeabilizing activity. Further, the drug hexamethylene amiloride (HMA), but not amiloride, inhibited in vitro ion channel activity of some synthetic coronavirus E proteins, and also viral replication. We have previously shown for the coronavirus species responsible for severe acute respiratory syndrome (SARS-CoV) that the transmembrane domain of E protein (ETM) forms pentameric α-helical bundles that are likely responsible for the observed channel activity. Herein, using solution NMR in dodecylphosphatidylcholine micelles and energy minimization, we have obtained a model of this channel which features regular α-helices that form a pentameric left-handed parallel bundle. The drug HMA was found to bind inside the lumen of the channel, at both the C-terminal and the N-terminal openings, and, in contrast to amiloride, induced additional chemical shifts in ETM. Full length SARS-CoV E displayed channel activity when transiently expressed in human embryonic kidney 293 (HEK-293) cells in a whole-cell patch clamp set-up. This activity was significantly reduced by hexamethylene amiloride (HMA), but not by amiloride. The channel structure presented herein provides a possible rationale for inhibition, and a platform for future structure-based drug design of this potential pharmacological target.

Coronaviruses are viral pathogens that cause a variety of lethal diseases in birds and mammals, and common colds in humans. In 2003, however, an animal coronavirus was able to infect humans and produced severe acute respiratory syndrome (SARS), causing a near pandemic. Such events are likely to reoccur in the future, and new antiviral strategies are necessary. A small coronavirus protein called ‘envelope’ is important for pathogenesis, affecting the formation of the viral envelope and the distribution of the virus in the body. In vitro studies have shown that synthetic coronavirus envelope proteins have channel activity that in some cases has been inhibited by the drug hexamethylene amiloride, but not by amiloride. In the present paper, we have characterized the structure responsible for this channel activity. We have also determined the binding site of the drug hexamethylene amiloride in the channel, and shown that amiloride has only a mild effect on the NMR signals from the protein. The validity of these results is supported using mammalian cells expressing full length SARS-CoV E, where channel activity was inhibited by hexamethylene amiloride, but only mildly by amiloride. The structural model described for this channel provides a valuable insight into coronavirus envelope protein ion channel activity, and could serve as a platform for the development of novel anti-viral drugs.

Structural flexibility of the pentameric SARS coronavirus envelope protein ion channel

Coronaviruses contain a small envelope membrane protein with cation-selective ion channel activity mediated by its transmembrane domain (ETM). In a computational study, we proposed that ion channel activity can be explained by either of two similar ETM homopentameric transmembrane α-helical bundles, related by a ∼50° rotation of the helices. Later, we tested this prediction, using site-specific infrared dichroism of a lysine-flanked isotopically labeled ETM peptide from the virus responsible for the severe acute respiratory syndrome, SARS, reconstituted in lipid bilayers. However, the data were consistent with the presence of a kink at the center of the ETM α-helix, and it did not fit completely either computational model. Herein, we have used native ETM, without flanking lysines, and show that the helix orientation is now consistent with one of the predicted models. ETM only produced one oligomeric form, pentamers, in the lipid-mimic detergent dodecylphosphocholine and in perfluorooctanoic acid. We thus report the correct backbone model for the pentameric α-helical bundle of ETM. The disruptive effects caused by terminal lysines probably highlight the conformational flexibility required during ion channel function.

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.

A trimeric building block model for Cry toxins in vitro ion channel formation

The crystal (Cry) insecticidal toxins, or delta-endotoxins, are lethal to a wide variety of insect larvae, and are therefore very important in insect control. Toxicity has been explained by formation of transmembrane oligomeric pores or ion channels and, more recently, by the ability of the monomeric toxin to subvert cellular signaling pathways. The structure, topology, and precise role of the putative pore in toxicity are not known. However, in vitro biophysical studies suggest that helices alpha4 and alpha5 in domain I insert into the lipid bilayer as an alpha-helical hairpin. Mutagenesis studies have assigned an important role to alpha5 in maintaining oligomerization, and to alpha4 in channel formation. To detect the possible homo-oligomerizing tendencies of these two helices, we have used the evolutionary conservation data contained in sixteen Cry homologs in order to filter non-native interactions found during a global conformational search. No conserved homo-oligomer was found for alpha4, but a right handed trimeric alpha5 model was present in the simulations of all Cry sequences. We propose a model for Cry toxin oligomerization based on sequence analysis and available mutagenesis data.

Unambiguous prediction of human integrin transmembrane heterodimer interactions using only homologous sequences

Part of the interaction between the alpha- and beta-subunits of integrins is known to take place at the transmembrane (TM) domain, where both heteromeric and homomeric aggregates have been reported in vivo and in vitro. In a recent computational study, totally independent from biochemical or biophysical data, we explored the plausibility of various TM homo-oligomers using evolutionary conservation data as a filter for non-native interactions. We showed that several homodimeric and homotrimeric interactions for alpha- and beta-chains are evolutionarily conserved. We report herein the results of the application of the same exhaustive approach to the integrin heterodimer. We have studied all known human TM integrin alphabeta pairs, and we show unambiguously that two models of interaction are evolutionarily conserved. These two models are consistent with those proposed previously based on mutagenesis and crosslinking. Comparison with previous experimental data strongly supports that a glycophorin A-like model is an intermediate form of interaction between the resting state and the active form, where chain separation occurs. Surprisingly, these two models are also conserved when considering most of the possible alphabeta pair combinations, suggesting that specific pairing of integrins is not determined by the TM domain, which has remained unchanged in spite of the variety of known integrin functions. This fact highlights a common ancestral mechanism for signal transduction that has remained through evolution. In a broader context, our results show that it is possible to obtain correct and detailed interactions of alpha-helical heterodimers with total independence of experimental data.

Transmembrane Helix Packing of ErbB/Neu Receptor in Membrane Environment: A Molecular Dynamics Study

Dimerization or oligomerization of the ErbB/Neu receptors are necessary but not sufficient for initiation of receptor signaling. The two intracellular domains must be properly oriented for the juxtaposition of the kinase domains allowing trans-phosphorylation. This suggests that the transmembrane (TM) domain acts as a guide for defining the proper orientation of the intracellular domains. Two structural models, with the two helices either in left-handed or in right-handed coiling have been proposed as the TM domain structure of the active receptor. Because experimental data do not distinguish clearly helix-helix packing, molecular dynamics (MD) simulations are used to investigate the energetic factors that drive Neu TM-TM interactions of the wild and the oncogenic receptor (Val664/Glu mutation) in DMPC or in POPC environments. MD results indicate that helix-lipid interactions in the bilayer core are extremely similar in the two environments and raise the role of the juxtamembrane residues in helix insertion and helix-helix packing. The TM domain shows a greater propensity to adopt a left-handed structure in DMPC, with helices in optimal position for strong inter-helical Hbonds induced by the Glu mutation. In POPC, the right-handed structure is preferentially formed with the participation of water in inter-helical Hbonds. The two structural arrangements of the Neu(TM) helices both with GG4 residue motif in close contact at the interface are permissible in the membrane environment. According to the hypothesis of a monomer-dimer equilibrium of the proteins it is likely that the bilayer imposes structural constraints that favor dimerization-competent structure responsible of the proper topology necessary for receptor activation.

Membrane assembly of simple helix homo-oligomers studied via molecular dynamics simulations

The assembly of simple transmembrane helix homo-oligomers is studied by combining a generalized Born implicit membrane model with replica exchange molecular dynamics simulations to sample the conformational space of various oligomerization states and the native oligomeric conformation. Our approach is applied to predict the structures of transmembrane helices of three proteins--glycophorin A, the M2 proton channel, and phospholamban--using only peptide sequence and the native oligomerization state information. In every case, the methodology reproduces native conformations that are in good agreement with available experimental structural data. Thus, our method should be useful in the prediction of native structures of transmembrane domains of other peptides. When we ignore the experimental constraint on the native oligomerization state and attempt de novo prediction of the structure and oligomerization state based only on sequence and simple energetic considerations, we identify the pentamer as the most stable oligomer for phospholamban. However, for the glycophorin A and the M2 proton channels, we tend to predict higher oligomers as more stable. Our studies demonstrate that reliable predictions of the structure of transmembrane helical oligomers can be achieved when the observed oligomerization state is imposed as a constraint, but that further efforts are needed for the de novo prediction of both structure and oligomeric state.

Sequence dependence of BNIP3 transmembrane domain dimerization implicates side-chain hydrogen bonding and a tandem GxxxG motif in specific helix-helix interactions

The transmembrane domain of the pro-apoptotic protein BNIP3 self-associates strongly in membranes and in detergents. We have used site-directed mutagenesis to analyze the sequence dependence of BNIP3 transmembrane domain dimerization, from which we infer the physical basis for strong and specific helix-helix interactions in this system. Hydrophobic substitutions identify six residues as critical to dimerization, and the pattern of sensitive residues suggests that the BNIP3 helices interact at a right-handed crossing angle. Based on the dimerization propensities of single point mutants, we propose that: polar residues His173 and Ser172 make inter-monomer hydrogen bonds to one another through their side-chains; Ala176, Gly180, and Gly184 form a tandem GxxxG motif that allows close approach of the helices; and Ile183 makes inter-monomer van der Waals contacts. Since neither the tandem GxxxG motif nor the hydrogen bonding pair is sufficient to drive dimerization, our results demonstrate the importance of sequence context for either hydrogen bonding or GxxxG motif involvement in BNIP3 transmembrane helix-helix interactions. In this study, hydrophobic substitutions away from the six interfacial positions have almost no effect on dimerization, confirming the expectation that hydrophobic replacements affect helix-helix interactions only if they interfere with packing or hydrogen bonding by interfacial residues. However, changes to slightly polar residues are somewhat disruptive even when located away from the interface, and the degree of disruption correlates with the decrease in hydrophobicity. Changing the hydrophobicity of the BNIP3 transmembrane domain alters its helicity and protection of its backbone amides. We suggest that polar substitutions decrease the fraction of dimer by stabilizing an unfolded monomeric state of the transmembrane span, rather than by affecting helix-helix interactions. This result has broad implications for interpreting the sequence dependence of membrane protein stability in detergents.

Prediction of transmembrane helix orientation in polytopic membrane proteins

Background

Membrane proteins compose up to 30% of coding sequences within genomes. However, their structure determination is lagging behind compared with soluble proteins due to the experimental difficulties. Therefore, it is important to develop reliable computational methods to predict structures of membrane proteins.

Results

We present a method for prediction of the TM helix orientation, which is an essential step in ab initio modeling of membrane proteins. Our method is based on a canonical model of the heptad repeat originally developed for coiled coils. We identify the helical surface patches that interface with lipid molecules at an accuracy of about 88% from the sequence information alone, using an empirical scoring function LIPS (LIPid-facing Surface), which combines lipophilicity and conservation of residues in the helix. We test and discuss results of prediction of helix-lipid interfaces on 162 transmembrane helices from 18 polytopic membrane proteins and present predicted orientations of TM helices in TRPV1 channel. We also apply our method to two structures of homologous cytochrome b₆f complexes and find discrepancy in the assignment of TM helices from subunits PetG, PetN and PetL. The results of LIPS calculations and analysis of packing and H-bonding interactions support the helix assignment found in the cytochrome b₆f structure from green alga but not the assignment of TM helices in the cyanobacterium b₆f structure.

Conclusion

LIPS calculations can be used for the prediction of helix orientation in ab initio modeling of polytopic membrane proteins. We also show with the example of two cytochrome b₆f structures that our method can identify questionable helix assignments in membrane proteins. The LIPS server is available online at .

Two types of transmembrane homomeric interactions in the integrin receptor family are evolutionarily conserved

Integrins are heterodimers, but recent in vitro and in vivo experiments suggest that they are also able to associate through their transmembrane domains to form homomeric interactions. Two fundamental questions are the biological relevance of these aggregates and their form of interaction in the membrane domain. Although in vitro experiments have shown the involvement of a GxxxG‐like motif, several crosslinking in vivo data are consistent with an almost opposite form of interaction between the transmembrane α‐helices. In the present work, we have explored these two questions using molecular dynamics simulations for all available integrin types. We have tested the hypothesis that homomeric interactions are evolutionary conserved, and essential for the cell, using conservative substitutions to filter out nonnative interactions. Our results show that two models, one involving a GxxxG‐like motif (model I) and an almost opposite form of interaction (model II) are conserved across all α and β integrin types, both in homodimers and homotrimers, with different specificities. No conserved interaction was found for homotetramers. Our results are completely independent from experimental data, both during molecular dynamics simulations and in the selection of the correct models. We rationalize previous seemingly conflicting findings regarding the nature of integrin interhelical homomeric interactions. Proteins 2006. © 2006 Wiley‐Liss, Inc.

Energetics of the native and non-native states of the glycophorin transmembrane helix dimer

Using an implicit membrane model (IMM1), we examine whether the structure of the transmembrane domain of Glycophorin A (GpA) could be predicted based on energetic considerations alone. The energetics of native GpA shows that van der Waals interactions make the largest contribution to stability. Although specific electrostatic interactions are stabilizing, the overall electrostatic contribution is close to zero. The GXXXG motif contributes significantly to stability, but residues outside this motif contribute almost twice as much. To generate non-native states a global conformational search was done on two segments of GpA: an 18-residue peptide (GpA74-91) that is embedded in the membrane and a 29-residue peptide (GpA70-98) that has additional polar residues flanking the transmembrane region. Simulated annealing was done on a large number of conformations generated from parallel, antiparallel, left- and right-handed starting structures by rotating each helix at 20 degrees intervals around its helical axis. Several crossing points along the helix dimer were considered. For 18-residue parallel topology, an ensemble of native-like structures was found at the lowest effective energy region; the effective energy is lowest for a right-handed structure with an RMSD of 1.0 A from the solid-state NMR structure with correct orientation of the helices. For the 29-residue peptide, the effective energies of several left-handed structures were lower than that of the native, right-handed structure. This could be due to deficiencies in modeling the interactions between charged sidechains and/or omission of the sidechain entropy contribution to the free energy. For 18-residue antiparallel topology, both IMM1 and a Generalized Born model give effective energies that are lower than that of the native structure. In contrast, the Poisson-Boltzmann solvation model gives lower effective energy for the parallel topology, largely because the electrostatic solvation energy is more favorable for the parallel structure. IMM1 seems to underestimate the solvation free energy advantage when the CO and NH dipoles just outside the membrane are parallel. This highlights the importance of electrostatic interactions even when these are not obvious by looking at the structures.

Energetics and stability of transmembrane helix packing: A replica-exchange simulation with a knowledge-based membrane potential

The energetics and stability of the packing of transmembrane helices were investigated by Monte Carlo simulations with the replica-exchange method. The helices were modeled with a united atom representation, and the CHARMM19 force field was employed. Based on known experimental structures of membrane proteins, an implicit knowledge-based potential was developed to describe the helix-membrane interactions at the residue level, whose validity was tested through prediction of the orientations when single helices were inserted into a membrane. Two systems were studied in this article, namely the glycophorin A dimer, and helices A and B of Bacteriorhodopsin. For the glycophorin A dimer, the most stable structure (0.5 A away from the experimental structure) is mainly stabilized by the favorable helix-helix interactions, and has the most population regardless of the helix-membrane interaction. However, for helices A and B of Bacteriorhodopsin, it was found that the packing determined by helix-helix interactions is nonspecific, and a native-like structure (0.2 A from the experimental one) can be identified from several structural analogs as the most stable one only after applying the membrane potential. Our results suggest that the contribution from the helix-membrane interaction could be critical in the correct packing of transmembrane helices in the membrane.

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.

Computer simulations of membrane proteins

Computer simulations are rapidly becoming a standard tool to study the structure and dynamics of lipids and membrane proteins. Increasing computer capacity allows unbiased simulations of lipid and membrane-active peptides. With the increasing number of high-resolution structures of membrane proteins, which also enables homology modelling of more structures, a wide range of membrane proteins can now be simulated over time spans that capture essential biological processes. Longer time scales are accessible by special computational methods. We review recent progress in simulations of membrane proteins.

Structure prediction of small transmembrane helix bundles

In this work, we will introduce a novel computational approach to predict the structures of small helical hetero-oligomeric transmembrane bundles. The approach is based on the generation and evaluation of a large library of randomly generated helix bundle conformations. This library will be evaluated by energy-dependent distributions of the structural parameters of the conformations. The approach enables us to model a subunit of cytochrome c oxidase (occ), consisting of four TM helices, to an accuracy of 1.7A and the transducer protein of the sensory Rhodopsin II-transducer complex to an accuracy of 2.3A when including two transducer-contacting Rhodopsin helices. As the approach does not afford a unique solution for each protein, experimental data would be needed to discriminate the possible models. In addition to predicting the structure of helix bundles, one can also gain insight into possible higher-energy conformations or flexible regions of the protein.

An implicit membrane generalized born theory for the study of structure, stability, and interactions of membrane proteins

Exploiting recent developments in generalized Born (GB) electrostatics theory, we have reformulated the calculation of the self-electrostatic solvation energy to account for the influence of biological membranes. Consistent with continuum Poisson-Boltzmann (PB) electrostatics, the membrane is approximated as an solvent-inaccessible infinite planar low-dielectric slab. The present membrane GB model closely reproduces the PB electrostatic solvation energy profile across the membrane. The nonpolar contribution to the solvation energy is taken to be proportional to the solvent-exposed surface area (SA) with a phenomenological surface tension coefficient. The proposed membrane GB/SA model requires minor modifications of the pre-existing GB model and appears to be quite efficient. By combining this implicit model for the solvent/bilayer environment with advanced computational sampling methods, like replica-exchange molecular dynamics, we are able to fold and assemble helical membrane peptides. We examine the reliability of this model and approach by applications to three membrane peptides: melittin from bee venom, the transmembrane domain of the M2 protein from Influenza A (M2-TMP), and the transmembrane domain of glycophorin A (GpA). In the context of these proteins, we explore the role of biological membranes (represented as a low-dielectric medium) in affecting the conformational changes in melittin, the tilt of transmembrane peptides with respect to the membrane normal (M2-TMP), helix-to-helix interactions in membranes (GpA), and the prediction of the configuration of transmembrane helical bundles (GpA). The present method is found to perform well in each of these cases and is anticipated to be useful in the study of folding and assembly of membrane proteins as well as in structure refinement and modeling of membrane proteins where a limited number of experimental observables are available.

Effective energy function for proteins in lipid membranes

A simple extension of the EEF1 energy function to heterogeneous membrane-aqueous media is proposed. The extension consists of (a) development of solvation parameters for a nonpolar phase using experimental data for the transfer of amino acid side-chains from water to cyclohexane, (b) introduction of a heterogeneous membrane-aqueous system by making the reference solvation free energy of each atom dependent on the vertical coordinate, (c) a modification of the distance-dependent dielectric model to account for reduced screening of electrostatic interactions in the membrane, and (d) an adjustment of the EEF1 aqueous model in light of recent calculations of the potential of mean force between amino acid side-chains in water. The electrostatic model is adjusted to match experimental observations for polyalanine, polyleucine, and the glycophorin A dimer. The resulting energy function (IMM1) reproduces the preference of Trp and Tyr for the membrane interface, gives reasonable energies of insertion into or adsorption onto a membrane, and allows stable 1-ns MD simulations of the glycophorin A dimer. We find that the lowest-energy orientation of melittin in bilayers varies, depending on the thickness of the hydrocarbon layer.

Membrane proteins: the 'Wild West' of structural biology

Historically, the task of determining the structure of membrane proteins has been hindered by experimental difficulties associated with their lipid-embedded domains. Here, we provide an overview of recently developed experimental and predictive tools that are changing our view of this largely unexplored territory - the 'Wild West' of structural biology. Crystallography, single-particle methods and atomic force microscopy are being used to study huge membrane proteins with increasing detail. Solid-state nuclear magnetic resonance strategies provide orientational constraints for structure determination of transmembrane (TM) alpha-helices and accurate measurements of intramolecular distances, even in very complex systems. Longer distance constraints are determined by site-directed spin-labelling electron paramagnetic resonance, but current labelling strategies still constitute some limitation. Other methods, such as site-specific infrared dichroism, enable orientational analysis of TM alpha-helices in aligned bilayers and, combined with novel computational and predictive tools that use evolutionary conservation data, are being used to analyze TM alpha-helical bundles.