Convergence of experimental, computational and evolutionary approaches predicts the presence of a tetrameric form for CD3-zeta

Use of Isotope-Edited FTIR to Derive a Backbone Structure of a Transmembrane Protein

Solving structures of membrane proteins has always been a formidable challenge, yet even upon success, the results are normally obtained in a mimetic environment that can be substantially different from a biological membrane. Herein, we use noninvasive isotope-edited FTIR spectroscopy to derive a structural model for the SARS coronavirus E protein transmembrane domain in lipid bilayers. Molecular-dynamics-based structural refinement, incorporating the IR-derived orientational restraints points to the formation of a helical hairpin structure. Disulfide cross-linking and X-ray reflectivity depth profiling provide independent support of the results. The unusually short helical hairpin structure of the protein might explain its ability to deform bilayers and is reminiscent of other peptides with membrane disrupting functionalities. Taken together, we show that isotope-edited FTIR is a powerful tool to analyze small membrane proteins in their native environment, enabling us to relate the unusual structure of the SARS E protein to its function.

Organization of the resting TCR in nanoscale oligomers

Despite the low affinity of the T-cell antigen receptor (TCR) for its peptide/major histocompatibility complex (pMHC) ligand, T cells are very sensitive to their antigens. This paradox can be resolved if we consider that the TCR may be organized into pre-existing oligomers or nanoclusters. Such structures could improve antigen recognition by increasing the functional affinity (avidity) of the TCR-pMHC interaction and by allowing cooperativity between individual TCRs. Up to approximately 20 TCRs become tightly apposed in these nanoclusters, often in a linear manner, and such structures could reflect a relatively generalized phenomenon: the non-random concentration of membrane receptors in specific areas of the plasma membrane known as protein islands. The association of TCRs into nanoclusters can explain the enhanced kinetics of the pMHC-TCR interaction in two dimensional versus three dimensional systems, but also their existence calls for a revision of the TCR triggering models based on pMHC-induced TCR clustering. Interestingly, the B-cell receptor and the FcεRI have also been shown to form nanoclusters, suggesting that the formation of pre-existing receptor oligomers could be widely used in the immune system.

Gaining insight into membrane protein structure using isotope-edited FTIR

FTIR spectroscopy has long been used as a tool used to gain average structural information on proteins. With the advent of stable isotope editing, FTIR can be used to derive accurate information on isolated amino acids. In particular, in an anisotropic sample such as membrane layers, it is possible to measure the orientation of the peptidic carbonyl groups. Herein, we review the theory that enables one to obtain accurate restraints from FTIR spectroscopy, alongside considerations for sample suitability and general applicability. We also propose approaches that may be used to generate structural models of simple membrane proteins based on FTIR orientational restraints. This article is part of a Special Issue entitled: FTIR in membrane proteins and peptide studies.

The integration of signaling and the spatial organization of the T cell synapse

Engagement of the T cell antigen receptor (TCR) triggers signaling pathways that lead to T cell selection, differentiation and clonal expansion. Superimposed onto the biochemical network is a spatial organization that describes individual receptor molecules, dimers, oligomers and higher order structures. Here we discuss recent findings and new concepts that may regulate TCR organization in naïve and memory T cells. A key question that has emerged is how antigen-TCR interactions encode spatial information to direct T cell activation and differentiation. Single molecule super-resolution microscopy may become an important tool in decoding receptor organization at the molecular level.

Environment Polarity in Proteins Mapped Noninvasively by FTIR Spectroscopy

The polarity pattern of a macromolecule is of utmost importance to its structure and function. For example, one of the main driving forces for protein folding is the burial of hydrophobic residues. Yet polarity remains a difficult property to measure experimentally, due in part to its non-uniformity in the protein interior. Herein, we show that FTIR linewidth analysis of noninvasive 1-(13)C=(18)O labels can be used to obtain a reliable measure of the local polarity, even in a highly multi-phasic system, such as a membrane protein. We show that in the Influenza M2 H(+) channel, residues that line the pore are located in an environment that is as polar as fully solvated residues, while residues that face the lipid acyl chains are located in an apolar environment. Taken together, FTIR linewidth analysis is a powerful, yet chemically non-perturbing approach to examine one of the most important properties in proteins - polarity.

Nanocluster formation: more with memory

The molecular basis of quicker and stronger responses by memory T cells is elusive. In this issue of Immunity, Kumar et al. (2011) demonstrate that more T cell receptor nanoclusters are present on memory cells than naive T cells before antigen stimulation, suggesting a basis for quick memory response.

Increased sensitivity of antigen-experienced T cells through the enrichment of oligomeric T cell receptor complexes

Although memory T cells respond more vigorously to stimulation and they are more sensitive to low doses of antigen than naive T cells, the molecular basis of this increased sensitivity remains unclear. We have previously shown that the T cell receptor (TCR) exists as different-sized oligomers on the surface of resting T cells and that large oligomers are preferentially activated in response to low antigen doses. Through biochemistry and electron microscopy, we now showed that previously stimulated and memory T cells have more and larger TCR oligomers at the cell surface than their naive counterparts. Reconstitution of cells and mice with a point mutant of the CD3ζ subunit, which impairs TCR oligomer formation, demonstrated that the increased size of TCR oligomers was directly responsible for the increased sensitivity of antigen-experienced T cells. Thus, we propose that an "avidity maturation" mechanism underlies T cell antigenic memory.

Computational study of the Na+/H + antiporter from Vibrio parahaemolyticus

Sodium proton antiporters are ubiquitous membrane proteins that catalyze the exchange of Na(+) for protons throughout the biological world. The Escherichia coli NhaA is the archetypal Na(+)/H(+) antiporter and is absolutely essential for survival in high salt concentrations under alkaline conditions. Its crystal structure, accompanied by extensive molecular dynamics simulations, have provided an atomically detailed model of its mechanism. In this study, we utilized a combination of computational methodologies in order to construct a structural model for the Na(+)/H(+) antiporter from the gram-negative bacterium Vibrio parahaemolyticus. We explored its overall architecture by computational means and validated its stability and robustness. This protein belongs to a novel group of NhaA proteins that transports not only Na(+) and Li(+) as substrate ions, but K(+) as well, and was also found to miss a β-hairpin segment prevalent in other homologs of the Bacteria domain. We propose, for the first time, a structure of a prototype model of a β-hairpin-less NhaA that is selective to K(+). Better understanding of the Vibrio parahaemolyticus NhaA structure-function may assist in studies on ion transport, pH regulation and designing selective blockers.

Empirical amide I vibrational frequency map: application to 2D-IR line shapes for isotope-edited membrane peptide bundles

The amide I vibrational mode, primarily associated with peptide-bond carbonyl stretches, has long been used to probe the structures and dynamics of peptides and proteins by infrared (IR) spectroscopy. A number of ab initio-based amide I vibrational frequency maps have been developed for calculating IR line shapes. In this paper, a new empirical amide I vibrational frequency map is developed. To evaluate its performance, we applied this map to a system of isotope-edited CD3-zeta membrane peptide bundles in aqueous solution. The calculated 2D-IR diagonal line widths vary from residue to residue and show an asymmetric pattern as a function of position in the membrane. The theoretical results are in fair agreement with experiments on the same system. Through analysis of the computed frequency time-correlation functions, it is found that the 2D-IR diagonal widths are dominated by contributions from the inhomogeneous frequency distributions, from which it follows that these widths are a good probe of the extent of local structural fluctuations. Thus, the asymmetric pattern of line widths follows from the asymmetric structure of the bundle in the membrane.

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.

Structural disorder of the CD3zeta transmembrane domain studied with 2D IR spectroscopy and molecular dynamics simulations

In a recently reported study [Mukherjee, et al. Proc. Natl. Acad. Sci. U.S.A. 2006, 103, 3528] we used 2D IR spectroscopy and 1-(13)C=(18)O isotope labeling to measure the vibrational dynamics of 11 amide I modes in the CD3zeta transmembrane domain. We found that the homogeneous line widths and population relaxation times were all nearly identical, but that the amount of inhomogeneous broadening correlated with the position of the amide group inside the membrane. In this study, we use molecular dynamics simulations to investigate the structural and dynamical origins of these experimental observations. We use two models to convert the simulations to frequency trajectories from which the mean frequencies, standard deviations, frequency correlation functions, and 2D IR spectra are calculated. Model 1 correlates the hydrogen-bond length to the amide I frequency, whereas model 2 uses an ab initio-based electrostatic model. We find that the structural distributions of the peptidic groups and their environment are reflected in the vibrational dynamics of the amide I modes. Environmental forces from the water and lipid headgroups partially denature the helices, shifting the infrared frequencies and creating larger inhomogeneous distributions for residues near the ends. The least inhomogeneously broadened residues are those located in the middle of the membrane where environmental electrostatic forces are weakest and the helices are most ordered. Comparison of the simulations to experiment confirms that the amide I modes near the C-terminal are larger than at the N-terminal because of the asymmetric structure of the peptide bundle in the membrane. The comparison also reveals that residues at a kink in the alpha-helices have broader line widths than more helical parts of the peptide because the peptide backbone at the kink exhibits a larger amount of structural disorder. Taken together, the simulations and experiments reveal that infrared line shapes are sensitive probes of membrane protein structural and environmental heterogeneity.

Isotope-edited IR spectroscopy for the study of membrane proteins

Fourier transform infrared (FTIR) spectroscopy has long been a powerful tool for structural analysis of membrane proteins. However, because of difficulties in resolving contributions from individual residues, most of the derived measurements tend to yield average properties for the system under study. Isotope editing, through its ability to resolve individual vibrations, establishes FTIR as a method that is capable of yielding accurate structural data on individual sites in a protein.

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 .

Model of a putative pore: the pentameric alpha-helical bundle of SARS coronavirus E protein in lipid bilayers

The coronavirus responsible for the severe acute respiratory syndrome contains a small envelope protein, E, with putative involvement in host apoptosis and virus morphogenesis. To perform these functions, it has been suggested that protein E can form a membrane destabilizing transmembrane (TM) hairpin, or homooligomerize to form a TM pore. Indeed, in a recent study we reported that the α-helical putative transmembrane domain of E protein (ETM) forms several SDS-resistant TM interactions: a dimer, a trimer, and two pentameric forms. Further, these interactions were found to be evolutionarily conserved. Herein, we have studied multiple isotopically labeled ETM peptides reconstituted in model lipid bilayers, using the orientational parameters derived from infrared dichroic data. We show that the topology of ETM is consistent with a regular TM α-helix. Further, the orientational parameters obtained unequivocally correspond to a homopentameric model, by comparison with previous predictions. We have independently confirmed that the full polypeptide of E protein can also aggregate as pentamers after expression in Escherichia coli. This interaction must be stabilized, at least partially, at the TM domain. The model we report for this pentameric α-helical bundle may explain some of the permabilizing properties of protein E, and should be the basis of mutagenesis efforts in future functional studies.

Picosecond dynamics of a membrane protein revealed by 2D IR

Fast protein dynamics can be missed with techniques that have relatively slow observation times. Using 2D IR spectroscopy and isotope labeling, we have probed the rapid, picosecond dynamics of a membrane protein in its native environment. By measuring the homogeneous and inhomogeneous IR linewidths of 11 amide I modes (backbone carbonyl stretch), we have captured the structural distributions and dynamics of the CD3zeta protein along its transmembrane segment that are lost with slower time-scale techniques. We find that the homogeneous lifetimes and population relaxation times are the same for almost all of the residues. In contrast, the inhomogeneous linewidths vary significantly with the largest inhomogeneous distribution occurring for residues near the N terminus and the narrowest near the center. This behavior is highly consistent with a recently reported experimental model of the protein and water accessibility as observed by molecular dynamics simulations. The data support the proposed CD3zeta peptide structure, and the simulations point to the structural disorder of water and lipid head-groups as the main source of inhomogeneous broadening. Taken together, this rigorous analysis of the vibrational dynamics of a membrane peptide provides experimental insight into a time regime of motions that has so far been largely unexplored.

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.

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

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

Site-specific dichroism analysis utilizing transmission FTIR

Infrared spectroscopy has long been used to examine the average secondary structure and orientation of membrane proteins. With the recent utilization of site-specific isotope labeling (e.g., peptidic 1-(13)C = (18)O) it is now possible to examine localized properties, rather than global averages. The technique of site-specific infrared dichroism (SSID) capitalized on this fact, and derives site-specific orientational restraints for the labeled amino acids. These restraints can then be used to solve the backbone structure of simple alpha-helical bundles, emphasizing the capabilities of this approach. So far SSID has been carried out in attenuated total internal reflection optical mode, with all of the respective caveats of attenuated total internal reflection. In this report we extend SSID through the use of transmission infrared spectroscopy tilt series. We develop the corresponding theory and demonstrate that accurate site-specific orientational restraints can be derived from a simple transmission experiment.

Oligomerization, Biogenesis, and Signaling Is Promoted by a Glycophorin A-like Dimerization Motif in Transmembrane Domain 1 of a Yeast G Protein-coupled Receptor

G protein-coupled receptors (GPCRs) can form dimeric or oligomeric complexes in vivo. However, the functions and mechanisms of oligomerization remain poorly understood for most GPCRs, including the alpha-factor receptor (STE2 gene product) of the yeast Saccharomyces cerevisiae. Here we provide evidence indicating that alpha-factor receptor oligomerization involves a GXXXG motif in the first transmembrane domain (TM1), similar to the transmembrane dimerization domain of glycophorin A. Results of fluorescence resonance energy transfer, fluorescence microscopy, endocytosis assays of receptor oligomerization in living cells, and agonist binding assays indicated that amino acid substitutions affecting the glycine residues of the GXXXG motif impaired alpha-factor receptor oligomerization and biogenesis in vivo but did not significantly impair agonist binding affinity. Mutant receptors exhibited signaling defects that were not due to impaired cell surface expression, indicating that oligomerization promotes alpha-factor receptor signal transduction. Structure-function studies suggested that the GXXXG motif in TM1 of the alpha-factor receptor promotes oligomerization by a mechanism similar to that used by the GXXXG dimerization motif of glycophorin A. In many mammalian GPCRs, motifs related to the GXXXG sequence are present in TM1 or other TM domains, suggesting that similar mechanisms are used by many GPCRs to form dimers or oligomeric arrays.

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.

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.

Multiple site-specific infrared dichroism of CD3-zeta, a transmembrane helix bundle

The structure of the transmembrane domain of CD3-zeta a component of the T-cell receptor involved in signal transduction, has been studied in its native state (a lipid bilayer) by multiple site-specific infrared dichroism. For the first time, the transmembrane domain has been labelled at multiple positions along the sequence, representing a total of 11 samples, each labelled at a different residue with an isotopically modified carbonyl group, (13)C [double bond] (18)O. A strategy is outlined that, based on the above data, can yield the rotational orientation and the local helix tilt for each labelled residue, giving a detailed description of helix geometry. The results obtained indicate that the transmembrane segment is in an alpha-helical conformation throughout, with an average helix tilt of 12 degrees. The N-terminal side of the helix is more tilted than the C-terminal. In an accompanying paper we describe the implementation of the infrared data in a model-building study of the CD3-zeta transmembrane complex. The model obtained is entirely consistent with results based on evolutionary conservation data. Taken together, this study represents the first step towards elucidation of the backbone structure of a transmembrane alpha-helical bundle by infrared spectroscopy.