Accessibility of introduced cysteines in chemoreceptor transmembrane helices reveals boundaries interior to bracketing charged residues

Two hydrophobic sequences, 24 and 30 residues long, identify the membrane-spanning segments of chemoreceptor Trg from Escherichia coli. As in other related chemoreceptors, these helical sequences are longer than the minimum necessary for an alpha-helix to span the hydrocarbon region of a biological membrane. Thus, the specific positioning of the segments relative to the hydrophobic part of the membrane cannot be deduced from sequence alone. With the aim of defining the positioning for Trg experimentally, we determined accessibility of a hydrophilic sulfhydryl reagent to cysteines introduced at each position within and immediately outside the two hydrophobic sequences. For both sequences, there was a specific region of uniformly low accessibility, bracketed by regions of substantial accessibility. The two low-accessibility regions were each 19 residues long and were in register in the three-dimensional organization of the transmembrane domain deduced from independent data. None of the four hydrophobic-hydrophilic boundaries for these two membrane-embedded sequences occurred at a charged residue. Instead, they were displaced one to seven residues internal to the charged side chains bracketing the extended hydrophobic sequences. Many hydrophobic sequences, known or predicted to be membrane-spanning, are longer than the minimum necessary helical length, but precise membrane boundaries are known for very few. The cysteine-accessibility approach provides an experimental strategy for determining those boundaries that could be widely applicable.     

Sequence context strongly modulates association of polar residues in transmembrane helices

Polar residues are capable of mediating the association of membrane-embedded helices through the formation of side-chain/side-chain inter-helical hydrogen bonds. However, the extent to which native van der Waals packing of the residues surrounding the polar locus can enhance, or interfere with, the interaction of polar residues has not yet been studied. We examined the propensities of four polar residues (aspartic acid, asparagine, glutamic acid, and glutamine) to promote self-association of transmembrane (TM) domains in several biologically derived sequence environments, including (i). four naturally occurring TM domains that contain a Glu or Gln residue (Tnf5/CD40 ligand, C79a/Ig-alpha, C79b/Ig-beta, and Fut3/alpha-fucosyltransferase); and (ii). variants of bacteriophage M13 major coat protein TM segment with Asp and Asn at interfacial and non-interfacial positions. Self-association was quantified by the TOXCAT assay, which measures TM helix self-oligomerization in the Escherichia coli inner membrane. While an appropriately placed polar residue was found in several cases to significantly stabilize TM helix-helix interactions through the formation of an interhelical hydrogen bond, in other cases the strongly polar residues did not enhance the association of the two helices. Overall, these results suggest that an innate structural mechanism may operate to control non-specific association of membrane-embedded polar residues.     

Proline residues in transmembrane alpha helices affect the folding of bacteriorhodopsin

Proline residues occur frequently in transmembrane alpha helices, which contrasts with their behaviour as helix-breakers in water-soluble proteins. The three membrane-embedded proline residues of bacteriorhodopsin have been replaced individually by alanine and glycine to give P50A, or P50G on helix B, P91A, or P91G on helix C, and P186A or P186G on helix F, and the effect on the protein folding kinetics has been investigated. The rate-limiting apoprotein folding step, which results in formation of a seven transmembrane, alpha helical state, was slower than wild-type protein for the Pro50 and Pro91 mutants, regardless of whether they were mutated to Ala or Gly. These proline residues give rise to several inter-helix contacts, which are therefore important in folding to the seven transmembrane helix state. No evidence for cis-trans isomerisations of the peptidyl prolyl bonds was found during this rate-limiting apoprotein folding step. Mutations of all three membrane-embedded proline residues affected the subsequent retinal binding and final folding to bacteriorhodopsin, suggesting that these proline residues contribute to formation of the retinal binding pocket within the helix bundle, again via helix/helix interactions. These results point to proline residues in transmembrane alpha helices being important in the folding of integral membrane proteins. The helix/helix interactions and hydrogen bonds that arise from the presence of proline residues in transmembrane alpha helices can affect the formation of transmembrane alpha helix bundles as well as cofactor binding pockets.     

Detection of drug‐induced conformational change of a transmembrane protein in lipid bilayers using site‐directed spin labeling

As a target of antiviral drugs, the influenza A M2 protein has been the focus of numerous structural studies and has been extensively explored as a model ion channel. In this study, we capitalize on the expanding body of high‐resolution structural data available for the M2 protein to design and interpret site‐directed spin‐labeling electron paramagnetic resonance spectroscopy experiments on drug‐induced conformational changes of the M2 protein embedded in lipid bilayers. We obtained data in the presence of adamantane drugs for two different M2 constructs (M2TM 22–46 and M2TMC 23–60). M2TM peptides were spin labeled at the N‐terminal end of the transmembrane domain. M2TMC peptides were spin labeled site specifically at cysteine residues substituted for amino acids within the transmembrane domain (L36, I39, I42, and L43) and the C‐terminal amphipathic helix (L46, F47, F48, C50, I51, Y52, R53, F54, F55, and E56). Addition of adamantane drugs brought about significant changes in measured electron paramagnetic resonance spectroscopy environmental parameters consistent with narrowing of the transmembrane channel pore and closer packing of the C‐terminal amphipathic helices.     

Chemical shift assignment of the transmembrane helices of DsbB, a 20-kDa integral membrane enzyme, by 3D magic-angle spinning NMR spectroscopy

The Escherichia coli inner membrane enzyme DsbB catalyzes disulfide bond formation in periplasmic proteins, by transferring electrons to ubiquinone from DsbA, which in turn directly oxidizes cysteines in substrate proteins. We have previously shown that DsbB can be prepared in a state that gives highly resolved magic-angle spinning (MAS) NMR spectra. Here we report sequential 13C and 15N chemical shift assignments for the majority of the residues in the transmembrane helices, achieved by three-dimensional (3D) correlation experiments on a uniformly 13C, 15N-labeled sample at 750-MHz 1H frequency. We also present a four-dimensional (4D) correlation spectrum, which confirms assignments in some highly congested regions of the 3D spectra. Overall, our results show the potential to assign larger membrane proteins using 3D and 4D correlation experiments and form the basis of further structural and dynamical studies of DsbB by MAS NMR.     

Initial structural and dynamic characterization of the M2 protein transmembrane and amphipathic helices in lipid bilayers

Amphipathic helices in membrane proteins that interact with the hydrophobic/hydrophilic interface of the lipid bilayer have been difficult to structurally characterize. Here, the backbone structure and orientation of an amphipathic helix in the full-length M2 protein from influenza A virus has been characterized. The protein has been studied in hydrated DMPC/DMPG lipid bilayers above the gel to liquid-crystalline phase transition temperature by solid-state NMR spectroscopy. Characteristic PISA (Polar Index Slant Angle) wheels reflecting helical wheels have been observed in uniformly aligned bilayer preparations of both uniformly 15N labeled and amino acid specific labeled M2 samples. Hydrogen/deuterium exchange studies have shown the very slow exchange of some residues in the amphipathic helix and more rapid exchange for the transmembrane helix. These latter results clearly suggest the presence of an aqueous pore. A variation in exchange rate about the transmembrane helical axis provides additional support for this claim and suggests that motions occur about the helical axes in this tetramer to expose the entire backbone to the pore.     

The role of proline residues in the dynamics of transmembrane helices: the case of bacteriorhodopsin

Proline residues in transmembrane helices have been found to have important roles in the functioning of membrane proteins. Moreover, Pro residues occur with high frequency in transmembrane α-helices, as compared to α-helices for soluble proteins. Here, we report several properties of the bacteriorhodopsin mutants P50A (helix B), P91A (helix C) and P186A (helix F). Compared to wild type, strongly perturbed behaviour has been found for these mutants. In the resting state, increased hydroxylamine accessibility and altered Asp-85 pK_a and light-dark adaptation were observed. On light activation, hydroxylamine accessibility was increased and proton transport activity, M formation kinetics and FTIR difference spectra of M and N intermediates showed clear distortions. On the basis of these alterations and the near identity of the crystalline structures of mutants with that of wild type, we conclude that the transmembrane proline residues of bacteriorhodopsin fulfil a dynamic role in both the resting and the light-activated states. Our results are consistent with the notion that mutation of Pro to Ala allows the helix to increase its flexibility towards the direction originally hindered by the steric clash between the ring Cγ and the carbonyl O of the i-4 residue, at the same time decreasing the mobility towards the opposite direction. Due to their properties, transmembrane Pro residues may serve as transmission elements of conformational changes during the transport process. We propose that these concepts can be extended to other transmembrane proteins.     

The strong dimerization of the transmembrane domain of the fibroblast growth factor receptor (FGFR) is modulated by C‐terminal juxtamembrane residues

The fibroblast growth factor receptor 3 (FGFR3) is a member of the FGFR subfamily of the receptor tyrosine kinases (RTKs) involved in signaling across the plasma membrane. Generally, ligand binding leads to receptor dimerization and activation. Dimerization involves the transmembrane (TM) domain, where mutations can lead to constitutive activation in certain cancer types and also in skeletal malformations. Thus, it has been postulated that FGFR homodimerization must be inherently weak to allow regulation, a feature reminiscent of α and β integrin TM interactions. However, we show herein that in FGFR3‐TM, four C‐terminal residues, CRLR, have a profound destabilizing effect in an otherwise strongly dimerizing TM peptide. In the absence of these four residues, the dimerizing propensity of FGFR3‐TM is comparable to glycophorin, as shown using various detergents. In addition, the expected enhanced dimerization induced by the mutation associated to the Crouzon syndrome A391E, was observed only when these four C‐terminal residues were present. In the absence of these four residues, A391E was dimer‐destabilizing. Finally, using site specific infrared dichroism and convergence with evolutionary conservation data, we have determined the backbone model of the FGFR3‐TM homodimer in model lipid bilayers. This model is consistent with, and correlates with the effects of, most known pathological mutations found in FGFR‐TM.     

The role of proline residues in the dynamics of transmembrane helices: the case of bacteriorhodopsin

Proline residues in transmembrane helices have been found to have important roles in the functioning of membrane proteins. Moreover, Pro residues occur with high frequency in transmembrane alpha-helices, as compared to alpha-helices for soluble proteins. Here, we report several properties of the bacteriorhodopsin mutants P50A (helix B), P91A (helix C) and P186A (helix F). Compared to wild type, strongly perturbed behaviour has been found for these mutants. In the resting state, increased hydroxylamine accessibility and altered Asp-85 pKa and light-dark adaptation were observed. On light activation, hydroxylamine accessibility was increased and proton transport activity, M formation kinetics and FTIR difference spectra of M and N intermediates showed clear distortions. On the basis of these alterations and the near identity of the crystalline structures of mutants with that of wild type, we conclude that the transmembrane proline residues of bacteriorhodopsin fulfil a dynamic role in both the resting and the light-activated states. Our results are consistent with the notion that mutation of Pro to Ala allows the helix to increase its flexibility towards the direction originally hindered by the steric clash between the ring Cgamma and the carbonyl O of the i-4 residue, at the same time decreasing the mobility towards the opposite direction. Due to their properties, transmembrane Pro residues may serve as transmission elements of conformational changes during the transport process. We propose that these concepts can be extended to other transmembrane proteins.     

Beta-branched residues adjacent to GG4 motifs promote the efficient association of glycophorin A transmembrane helices

Protein transmenembrane (TM) segments participating in helix-helix packing commonly contain small residue patterns (termed GG4 or "small-xxx-small" motifs) at i and i + 4 positions. Within many TM segments - such as the glycophorin A (GpA) sequence L75IxxGVxxGVxxT87- the G17y-xxx-Gly83 motif often occurs in combination with large, usually beta3-branched aliphatic residues at adjacent positions, typified here by Val30 and Val84 residues. To explore the importance of local P-branched character on GpA dimerization, we made systematic replacements to all 16 combinations of single or double Ile, Leu, and AIa residues at GpA TM Val/Val positions 80 and 84. Using the TOXCAT system to assay self-oligomerization in the Escherichia coli inner membrane--we observed that (i) combinations of Val and lie residues maintained, or improved dimerization levels; (ii) single Ala or Leu mutant combinations with Val or Ile maintained near-wild type dimerization affinities; and (iii) in the absence of beta-branching, i.e., Leu/Leu, Ala/Ala and Ala/Leu combinations, GpA dimerization was significantly diminished. An apparent capacity of lle-containing mutants to increase GpA dimerization versus WT likely arises from improved van der Waals packing (vs. Val) within the locus of helix contact, consistent with correlations we noted in lipid accessibility measurements. Examination of several synthetic peptides with sequences corresponding to selected GpA mutants (VV VI, IV II, and LL) confirmed their dimerization on sodium dodecylsulfate polyacrylamide gel electrophoresis (SDS-PAGE). The overall results reinforce the importance of a beta-branch-containing "ridge" residue to complement a "small-xxx-small groove" in promotion of TM-TM interactions.     

Effect of temperature on the NMR spectra of sonicated dispersions of isomers of dipalmitoylphosphatidylcholine

Sonicated dispersions of 1,2-dipalmitoylsn-glycero-3-phosphorylcholine and of 1,3-dipalmitoylglycero-2-phosphorylcholine were examined by proton nuclear magnetic resonance (NMR) as a function of temperature. The —(CH₂)_n)— peak in the spectrum of the sn-3-isomer of dipalmitoylphosphatidylcholine showed the characteristic dramatic changes in the peak intensity and width associated with the phase transition between the liquid crystalline and gel states of the phospholipid. This occurred over a 2–3°C temperature range with the midpoint of the transition at 38.5°C. With the 2-isomer the change in phase took place over a similar temperature range but the midpoint was at 33.8°C. This lower phase transition temperature is presumably the result of increased acyl chain mobility caused by the increased separation of the two acyl chains by the centre carbon of the glycerol backbone. The effect of sonication of the broadening of the range and lowering of the midpoint temperature of the phase transition from that of the corresponding unsonicated dispersions was similar with each isomer. This suggests that the overall geometry of the sonicated vesicles of the isomers is similar.     

Interactions of two transmembrane peptides in supported lipid bilayers studied by a P and N MAOSS NMR strategy

³¹P and ¹⁵N solid-state NMR with the magic angle-oriented sample spinning (MAOSS) strategy was used to investigate the effect of two model peptides on phospholipid bilayers mimicking biological membrane. One of the peptides, alamethicin, used as a reference of transmembrane alignment, has been shown to disrupt the lipid bilayer organisation, affecting the DMPC packaging. On the other hand, a α-helix alanine-rich peptide, K₃A₁₈K₃, with a ¹⁵N labelled alanine, did not present any effect in the DMPC bilayer organisation. The mean orientation of this peptide in the bilayer gave a transmembrane alignment of about 80%.     

Orientational constraints as three-dimensional structural constraints from chemical shift anisotropy: the polypeptide backbone of gramicidin A in a lipid bilayer.

Chemical shifts observed from samples that are uniformly aligned with respect to the magnetic field can be used as very high-resolution structural constraints. This constraint takes the form of an orientational constraint rather than the more familiar distance constraint. The accuracy of these constraints is dependent upon the quality of the tensor characterization. Both tensor element magnitudes and tensor orientations with respect to the molecular frame need to be considered. Here these constraints have been used to evaluate models for the channel conformation of gramicidin A. Of the three models used, the one experimentally derived model of gramicidin in sodium dodecyl sulfate micelles fits the data least well.     

The stability of transmembrane helices: A molecular dynamics study on the isolated helices of bacteriorhodopsin

Bacteriorhodopsin (bR) continues to be a proven testing ground for the study of integral membrane proteins (IMPs). It is important to study the stability of the individual helices of bR, as they are postulated to exist as independently stable transmembrane helices (TMHs) and also for their utility as templates for modeling other IMPs with the postulated seven-helix bundle topology. Toward this purpose, the seven helices of bR have been studied by molecular dynamics simulation in this study. The suitability of using the backbone-dependent rotamer library of side-chain conformations arrived at from the data base of globular protein structures in the case TMHs has been tested by another set of 7 helix simulations with the side-chain orientations taken from this library. The influence of the residue's net charge on the helix stability was examined by simulating the helices III, IV, and VI (from both of the above sets of helices) with zero net charge on the side chains. The results of these 20 simulations demonstrate in general the stability of the isolated helices of bR in conformity with the two-stage hypothesis of IMP folding. However, the helices I, II, V, and VII are more stable than the other three helices. The helical nature of certain regions of III, IV, and VI are influenced by factors such as the net charge and orientation of several residues. It is seen that the residues Arg, Lys, Asp, and Glu (charged residues), and Ser, Thr, Gly, and Pro, play a crucial role in the stability of the helices of bR. The backbone-dependent rotamer library for the side chains is found to be suitable for the study of TMHs in IMP. © 1996 John Wiley & Sons, Inc.     

Identification of histidine residues important in the catalysis and structure of aspartyl aminopeptidase

Aspartyl aminopeptidase (DAP), a widely distributed and abundant cytosolic enzyme, removes glutamyl or aspartyl residues from N-terminal acidic amino acid-containing peptides. DAP is a member of the M18 family of the MH clan of cocatalytic metallopeptidases. The human and mouse enzymes have been cloned. We have identified 8 highly homologous eukaryotic sequences that are probable aspartyl aminopeptidases. Eight histidine residues of human DAP were sequentially mutated to phenylalanine. Mutation of His94, His170, and His440 abolished enzymatic activity. His94 and His440 are postulated to be involved in binding cocatalytic zinc atoms by homology with other members of the MH clan. Mutation of His352 dramatically reduced enzyme activity. Gel-filtration analysis of the His352 mutant revealed destabilization of the quaternary structure and dissociation of the native 440-kDa enzyme. Mutation of His33 and of histidines residing in a cluster at residues 349, 359, and 363 all decreased k(cat). These studies reveal an important role for histidine residues both in catalysis and in the structural integrity of DAP.     

Solid-state NMR characterization of conformational plasticity within the transmembrane domain of the influenza A M2 proton channel

Membrane protein function within the membrane interstices is achieved by mechanisms that are not typically available to water-soluble proteins. The whole balance of molecular interactions that stabilize three-dimensional structure in the membrane environment is different from that in an aqueous environment. As a result interhelical interactions are often dominated by non-specific van der Waals interactions permitting dynamics and conformational heterogeneity in these interfaces. Here, solid-state NMR data of the transmembrane domain of the M2 protein from influenza A virus are used to exemplify such conformational plasticity in a tetrameric helical bundle. Such data lead to very high resolution structural restraints that can identify both subtle and substantial structural differences associated with various states of the protein. Spectra from samples using two different preparation protocols, samples prepared in the presence and absence of amantadine, and spectra as a function of pH are used to illustrate conformational plasticity.     

Estimating the length of transmembrane helices using Z-coordinate predictions

Zpred2 is an improved version of ZPRED, a predictor for the Z-coordinates of alpha-helical membrane proteins, that is, the distance of the residues from the center of the membrane. Using principal component analysis and a set of neural networks, Zpred2 analyzes data extracted from the amino acid sequence, the predicted topology, and evolutionary profiles. Zpred2 achieves an average accuracy error of 2.18 A (2.17 A when an independent test set is used), an improvement by 15% compared to the previous version. We show that this accuracy is sufficient to enable the predictions of helix lengths with a correlation coefficient of 0.41. As a comparison, two state-of-the-art HMM-based topology prediction methods manage to predict the helix lengths with a correlation coefficient of less than 0.1. In addition, we applied Zpred2 to two other problems, the re-entrant region identification and model validation. Re-entrants were able to be detected with a certain consistency, but not better than with previous approaches, while incorrect models as well as mispredicted helices of transmembrane proteins could be distinguished based on the Z-coordinate predictions.     

Differential repositioning of the second transmembrane helices from E. coli Tar and EnvZ upon moving the flanking aromatic residues

Aromatic tuning, i.e. repositioning aromatic residues found at the cytoplasmic end of transmembrane (TM) domains within bacterial receptors, has been previously shown to modulate signal output from the aspartate chemoreceptor (Tar) and the major osmosensor EnvZ of Escherichia coli. In the case of Tar, changes in signal output consistent with the vertical position of the native Trp-Tyr aromatic tandem within TM2 were observed. In contrast, within EnvZ, where a Trp-Leu-Phe aromatic triplet was repositioned, the surface that the triplet resided upon was the major determinant governing signal output. However, these studies failed to determine whether moving the aromatic residues was sufficient to physically reposition the TM helix within a membrane. Recent coarse-grained molecular dynamics (CG-MD) simulations predicted displacement of Tar TM2 upon moving the aromatic residues at the cytoplasmic end of the helix. Here, we demonstrate that repositioning the Trp-Tyr tandem within Tar TM2 displaces the C-terminal boundary of the helix relative to the membrane. In a similar analysis of EnvZ, an abrupt initial displacement of TM2 was observed but no subsequent movement was seen, suggesting that the vertical position of TM2 is not governed by the location of the Trp-Leu-Phe triplet. Our results also provide another set of experimental data, i.e. the resistance of EnvZ TM2 to being displaced upon aromatic tuning, which could be useful for subsequent refinement of the initial CG-MD simulations. Finally, we discuss the limitations of these methodologies, how moving flanking aromatic residues might impact steady-state signal output and the potential to employ aromatic tuning in other bacterial membrane-spanning receptors.     

Influence of assignment on the prediction of transmembrane helices in protein structures

α-Helical transmembrane proteins (TMP_α) are composed of a series of helices embedded in the lipid bilayer. Due to technical difficulties, few 3D structures are available. Therefore, the design of structural models of TMP_α is of major interest. We study the secondary structures of TMP_α by analyzing the influence of secondary structures assignment methods (SSAMs). For this purpose, a published and updated benchmark databank of TMP_α is used and several SSAMs (9) are evaluated. The analysis of the results points to significant differences in SSA depending on the methods used. Pairwise comparisons between SSAMs led to more than 10% of disagreement. Helical regions corresponding to transmembrane zones are often correctly characterized. The study of the sequence–structure relationship shows very limited differences with regard to the structural disagreement. Secondary structure prediction based on Bayes’ rule and using only a single sequence give correct prediction rates ranging from 78 to 81%. A structural alphabet approach gives a slightly better prediction, i.e., only 2% less than the best equivalent approach, whereas the prediction rate with a very different assignment bypasses 86%. This last result highlights the importance of the correct assignment choice to evaluate the prediction assessment.