Thermodynamics of glycophorin A transmembrane helix dimerization in C14 betaine micelles

We have used sedimentation equilibrium analytical ultracentrifugation to measure the free energy change for the glycophorin A transmembrane helix-helix dimerization in C14 betaine micelles. By varying the amount of micellar C14 betaine, we show that the protein association reaction in the micellar C14 phase behaves as an ideal-dilute solution. In this hydrophobic environment, the mole-fraction standard state free energy change for self-association of the SNGpA99 glycophorin A construct is -5.7 (+/-0.3, N=5) kcal mol(-1) at 25 degrees C. Compared with previous results carried out in C(8)E(5) micellar solutions, the free energy of dimerization is 1.3 kcal mol(-1) less favorable in C14 betaine micelles. In contrast, when considered on a per-interface basis, the formation of the glycophorin A transmembrane dimer in C14 betaine micelles may be more favorable than the association of several designed transmembrane peptides.     

The net orientation of nicotinic receptor transmembrane alpha-helices in the resting and desensitized states

The net orientation of nicotinic acetylcholine receptor transmembrane alpha-helices has been probed in both the activatable resting and nonactivatable desensitized states using linear dichroism Fourier-transform infrared spectroscopy. Infrared spectra recorded from reconstituted nicotinic acetylcholine receptor membranes after 72 h exposure to (2)H2O exhibit an intense amide I component band near 1655 cm(-1) that is due predominantly to hydrogen-exchange-resistant transmembrane peptides in an alpha-helical conformation. The measured dichroism of this band is 2.37, suggesting a net tilt of the transmembrane alpha-helices of roughly 40 degrees from the bilayer normal, although this value overestimates the tilt angle because the measured dichroism at 1655 cm(-1) also reflects the dichroism of overlapping amide I component bands. Significantly, no change in the net orientation of the transmembrane alpha-helices is observed upon agonist binding. In fact, the main changes in structure and orientation detected upon desensitization involve highly solvent accessible regions of the polypeptide backbone. Our data are consistent with a capping of the ligand binding site by the solvent accessible C-loop with little change in the structure of the transmembrane domain in the desensitized state. Changes in structure at the interface between the ligand-binding and transmembrane domains may uncouple binding from gating.     

Proline in alpha-helix: stability and conformation studied by dynamics simulation

Free-energy simulations have been used to estimate the change in the conformational stability of short polyalanine alpha-helices when one of the alanines is replaced by a proline residue. For substituting proline in the middle of the helix the change in free energy of folding (delta delta G degrees) was calculated as 14 kJ/mol (3.4 kcal/mol), in excellent agreement with the one available experimental value. The helix containing proline was found to be strongly kinked; the free energy for reducing the angle of the kink from 40 degrees to 15 degrees was calculated, and found to be small. A tendency to alternate hydrogen bonding schemes was observed in the proline-containing helix. These observations for the oligopeptide agree well with the observation of a range of kink angles (18-35 degrees) and variety of hydrogen bonding schemes, in the rare instances where proline occurs in helices in globular proteins. For substituting proline at the N-terminus of the helix the change in free energy of folding (delta delta G degrees) was calculated as -4 kJ/mol in the first helical position (N1) and +6 kJ/mol in the second helical position (N2). The observed frequent occurrence of proline in position N1 in alpha-helices in proteins therefore has its origin in stability differences of secondary structure. The conclusion reached here that proline may be a better helix former in position N1 than (even) alanine, and thus be a helix initiator may be testable experimentally by measurements of fraction helical conformation of individual residues in oligopeptides of appropriate sequence. The relevance of these results in regards to the frequent occurrence of proline-containing helices in certain membrane proteins is discussed.     

Why is the biological hydrophobicity scale more accurate than earlier experimental hydrophobicity scales?

The recognition of transmembrane helices by the translocon is primarily guided by the average hydrophobicity of the potential transmembrane helix. However, the exact hydrophobicity of each amino acid can be identified in several different ways. The free energy of transfer for amino acid analogues between a hydrophobic media, for example, octanol and water can be measured or obtained from simulations, the hydrophobicity can also be estimated by statistical properties from known transmembrane segments and finally the contribution of each amino acid type for the probability of translocon recognition has recently been measured directly. Although these scales correlate quite well, there are clear differences between them and it is not well understood which scale represents neither the biology best nor what the differences are. Here, we try to provide some answers to this by studying the ability of different scales to recognize transmembrane helices and predict the topology of transmembrane proteins. From this analysis it is clear that the biological hydrophobicity scale as well scales created from statistical analysis of membrane helices perform better than earlier experimental scales that are mainly based on measurements of amino acid analogs and not directly on transmembrane helix recognition. Using these results we identified the properties of the scales that perform better than other scales. We find, for instance, that the better performing scales consider proline more hydrophilic. This shows that transmembrane recognition is not only governed by pure hydrophobicity but also by the helix preferences for amino acids, as proline is a strong helix breaker. Proteins 2014; 82:2190–2198. © 2014 Wiley Periodicals, Inc.     

Proline-induced distortions of transmembrane helices

Proline residues in the transmembrane (TM) alpha-helices of integral membrane proteins have long been suspected to play a key role for helix packing and signal transduction by inducing regions of helix distortion and/or dynamic flexibility (hinges). In this study we try to characterise the effect of proline on the geometric properties of TM alpha-helices. We have examined 199 transmembrane alpha-helices from polytopic membrane proteins of known structure. After examining the location of proline residues within the amino acid sequences of TM helices, we estimated the helix axes either side of a hinge and hence identified a hinge residue. This enabled us to calculate helix kink and swivel angles. The results of this analysis show that proline residues occur with a significant concentration in the centre of sequences of TM alpha-helices. In this location, they may induce formation of molecular hinges, located on average about four residues N-terminal to the proline residue. A superposition of proline-containing TM helices structures shows that the distortion induced is anisotropic and favours certain relative orientations (defined by helix kink and swivel angles) of the two helix segments.     

Calculating the free energy of association of transmembrane helices

A large number of experimental studies have been devoted to quantifying the interaction between transmembrane (TM) helices in detergent micelles and, more recently, in bilayers. Theoretical calculation of association free energy of TM helices would be useful for predicting the propensity of given sequences to oligomerize and for understanding the difference between association in micelles and in bilayers. In this article, the theoretical foundation for calculating the standard association free energy of TM helices is laid out and is applied to glycophorin A in both micelles and bilayers. The standard association free energy is decomposed into the effective energy, translational, rotational, and conformational entropy terms. The effective energy of association is obtained by molecular dynamics simulations in an implicit membrane model. The translational and rotational entropy of association is calculated from the probability distribution of the translational and rotational degrees of freedom obtained from the molecular dynamics simulations. The side-chain conformational entropy of association is estimated from the probability distribution obtained by rigid rotation of all side-chain dihedral angles. The calculated standard association free energy of glycophorin A in N-dodecylphosphocholine micelles is in good agreement with the experimental value. The translational entropy cost is larger, whereas the rotational entropy cost is smaller in bilayers than in micelles. The standard association free energy in 1,2-dimyristoyl-sn-glycero-3-phosphocholine bilayers is calculated to be approximately 1.3 kcal/mol more favorable than in N-dodecylphosphocholine micelles, consistent with available experimental data.     

Comparison of intramolecular and intermolecular reactions in protein folding

The intrinsic component of the standard free energy change for the formation of a disulfide bond in a protein molecule is compared to that for an analogous chemical reaction. The former reaction, which represents theintramolecular formation of a disulfide bond in a protein molecule from a cysteine group containing a mixed disulfide bond with glutathione, and a free cysteine residue, is a unimolecular reaction. In contrast, its chemical analogue is a bimolecular reaction, and corresponds to theintermolecular disulfide interchange between a mixed disulfide-bonded compound between a cysteine residue and glutathione, and a free cysteine molecule. The difference in the intrinsic free energy of the above two reactions is estimated by two different approaches. First, a theoretical estimate of the magnitude of the difference in free energy of the two reactions (for a standard state of 1 M) is obtained using a gas-phase statistical thermodynamic approach, which indicates that the intramolecular reaction is energetically favored over its intermolecular counterpart by as much as 15.6 kcal/mole. For comparison, an experimentally derived value is also obtained, using experimental data from a study by Konishi et al. of the regeneration of the protein ribonuclease A (RNase A) from its reduced form by reduced and oxidized glutathiones. The intrinsic component of the free energy change of the intramolecular reaction, as it occurs in the protein molecule, is obtained from such experimental data by accounting explicitly for the free energy change (assumed to be solely an entropy change) pertaining to the conformational changes (ring closure) that the protein molecule undergoes in the course of the reaction. On the basis of the value derived from such an experimental approach, the intramolecular reaction is also energetically more favorable as compared to its intermolecular analogue, but only by a difference of 2.3 kcal/mole (for a standard state of 1 M). The large apparent discrepancy between the two values estimated from the theoretical and experimental approaches is rationalized by the postulation of several additional factors not inherent in the gas-phase theoretical estimate, such as dehydration and intramolecular hydrogen-bonding effects, which can largely compensate for the otherwise favorable energetics of the intramolecular reaction.     

Coarse-grained molecular dynamics simulations of the energetics of helix insertion into a lipid bilayer

Experimental and computational studies have indicated that hydrophobicity plays a key role in driving the insertion of transmembrane alpha-helices into lipid bilayers. Molecular dynamics simulations allow exploration of the nature of the interactions of transmembrane alpha-helices with their lipid bilayer environment. In particular, coarse-grained simulations have considerable potential for studying many aspects of membrane proteins, ranging from their self-assembly to the relation between their structure and function. However, there is a need to evaluate the accuracy of coarse-grained estimates of the energetics of transmembrane helix insertion. Here, three levels of complexity of model system have been explored to enable such an evaluation. First, calculated free energies of partitioning of amino acid side chains between water and alkane yielded an excellent correlation with experiment. Second, free energy profiles for transfer of amino acid side chains along the normal to a phosphatidylcholine bilayer were in good agreement with experimental and atomistic simulation studies. Third, estimation of the free energy profile for transfer of an arginine residue, embedded within a hydrophobic alpha-helix, to the center of a lipid bilayer gave a barrier of approximately 15 kT. Hence, there is a substantial barrier to membrane insertion for charged amino acids, but the coarse-grained model still underestimates the corresponding free energy estimate (approximately 29 kT) from atomistic simulations (Dorairaj, S., and Allen, T. W. (2007) Proc. Natl. Acad. Sci. U.S.A. 104, 4943-4948). Coarse-grained simulations were then used to predict the free energy profile for transfer of a simple model transmembrane alpha-helix (WALP23) across a lipid bilayer. The results indicated that a transmembrane orientation was favored by about -70 kT.     

Interhelical hydrogen bonding drives strong interactions in membrane proteins

Polar residues in transmembrane alpha-helices may strongly influence the folding or association of integral membrane proteins. To test whether a motif that promotes helix association in a soluble protein could do the same within a membrane, we designed a model transmembrane helix based on the GCN4 leucine zipper. We found in both detergent micelles and biological membranes that helix association is driven strongly by asparagine, independent of the rest of the hydrophobic leucine and/or valine sequence. Hydrogen bonding between membrane helices gives stronger associations than the packing of surfaces in glycophorin A helices, creating an opportunity to stabilize structures, but also implying a danger that non-specific interactions might occur. Thus, membrane proteins may fold to avoid exposure of strongly hydrogen bonding groups at their lipid exposed surfaces.     

The stability of transmembrane helix interactions measured in a biological membrane

Despite some promising progress in the understanding of membrane protein folding and assembly, there is little experimental information regarding the thermodynamic stability of transmembrane helix interactions and even less on the stability of transmembrane helix-helix interactions in a biological membrane. Here we describe an approach that allows quantitative measurement of transmembrane helix interactions in a biological membrane, and calculation of changes in the interaction free energy resulting from substitution of single amino acids. Dimerization of several variants of the glycophorin A transmembrane domain are characterized and compared to the wild-type (wt) glycophorin A transmembrane helix dimerization. The calculated DeltaDeltaG(app) values are further compared with values found in the literature. In addition, we compare interactions between the wt glycophorin A transmembrane domain and helices in which critical glycine residues are replaced by alanine or serine, respectively. The data demonstrate that replacement of the glycine residues by serine is less destabilizing than replacement by alanine with a DeltaDeltaG(app) value of about 0.4 kcal/mol. Our study comprises the first measurement of a transmembrane helix interaction in a biological membrane, and we are optimistic that it can be further developed and applied.     

Motifs of two small residues can assist but are not sufficient to mediate transmembrane helix interactions

Both experimental and statistical searches for specific motifs that mediate transmembrane helix-helix interactions showed that two glycine residues separated by three intervening residues (GxxxG) provide a framework for specific interactions. Further work suggested that other motifs of small residues can mediate the interaction of transmembrane domains, so that the AxxxA-motif could also drive strong interactions of alpha-helices in soluble proteins. Thus, all these data indicate that a motif of two small residues in a distance of four might be enough to provide a framework for transmembrane helix-helix interaction. To test whether GxxxG is equivalent to (small)xxx(small), we investigated the effect of a substitution of either of the two Gly residues in the glycophorin A GxxxG-motif by Ala or Ser using the recently developed GALLEX system. The results of this mutational study demonstrate that, while a replacement of either of the two Gly by Ala strongly disrupts GpA homo-dimerization, the mutation to Ser partly stabilizes a dimeric structure. We suggest that the Ser residue can form a hydrogen bond with a backbone carbonyl group of the adjacent helix stabilizing a preformed homo-dimer. While (small)xxx(small) serves as a useful clue, the context of adjacent side-chains is essential for stable helix interaction, so each case must be tested.     

Effects of alanine substitutions in alpha-helices of sperm whale myoglobin on protein stability.

The peptide backbones in folded native proteins contain distinctive secondary structures, alpha-helices, beta-sheets, and turns, with significant frequency. One question that arises in folding is how the stability of this secondary structure relates to that of the protein as a whole. To address this question, we substituted the alpha-helix-stabilizing alanine side chain at 16 selected sites in the sequence of sperm whale myoglobin, 12 at helical sites on the surface of the protein, and 4 at obviously internal sites. Substitution of alanine for bulky side chains at internal sites destabilizes the protein, as expected if packing interactions are disrupted. Alanine substitutions do not uniformly stabilize the protein, either in capping positions near the ends of helices or at mid-helical sites near the surface of myoglobin. When corrected for the extent of exposure of each side chain replaced by alanine at a mid-helix position, alanine replacement still has no clear effect in stabilizing the native structure. Thus linkage between the stabilization of secondary structure and tertiary structure in myoglobin cannot be demonstrated, probably because of the relatively small free energy differences between side chains in stabilizing isolated helix. By contrast, about 80% of the variance in free energy observed can be accounted for by the loss in buried surface area of the native residue substituted by alanine. The differential free energy of helix stabilization does not account for any additional variation.     

Exploring peptide-membrane interactions with coarse-grained MD simulations

The interaction of α-helical peptides with lipid bilayers is central to our understanding of the physicochemical principles of biological membrane organization and stability. Mutations that alter the position or orientation of an α-helix within a membrane, or that change the probability that the α-helix will insert into the membrane, can alter a range of membrane protein functions. We describe a comparative coarse-grained molecular dynamics simulation methodology, based on self-assembly of a lipid bilayer in the presence of an α-helical peptide, which allows us to model membrane transmembrane helix insertion. We validate this methodology against available experimental data for synthetic model peptides (WALP23 and LS3). Simulation-based estimates of apparent free energies of insertion into a bilayer of cystic fibrosis transmembrane regulator-derived helices correlate well with published data for translocon-mediated insertion. Comparison of values of the apparent free energy of insertion from self-assembly simulations with those from coarse-grained molecular dynamics potentials of mean force for model peptides, and with translocon-mediated insertion of cystic fibrosis transmembrane regulator-derived peptides suggests a nonequilibrium model of helix insertion into bilayers.     

Changes in apparent free energy of helix-helix dimerization in a biological membrane due to point mutations

We present an implementation of the TOXCAT membrane protein self-association assay that measures the change in apparent free energy of transmembrane helix dimerization caused by point mutations. Quantifying the reporter gene expression from cells carrying wild-type and mutant constructs shows that single point mutations that disrupt dimerization of the transmembrane domain of glycophorin A reproducibly lower the TOXCAT signal more than 100-fold. Replicate cultures can show up to threefold changes in the level of expression of the membrane bound fusion construct, and correcting for these variations improves the precision of the calculated apparent free energy change. The remarkably good agreement between our TOXCAT apparent free energy scale and free energy differences from sedimentation equilibrium studies for point mutants of the glycophorin A transmembrane domain dimer indicate that sequence changes usually affect membrane helix-helix interactions quite similarly in these two very different environments. However, the effects of point mutations at threonine 87 suggest that intermonomer polar contacts by this side-chain contribute significantly to dimer stability in membranes but not in detergents. Our findings demonstrate that a comparison of quantitative measurements of helix-helix interactions in biological membranes and genuine thermodynamic data from biophysical measurements on purified proteins can elucidate how changes in the lipidic environment modulate membrane protein stability.     

Stabilizing nonpolar/polar side-chain interactions in the alpha-helix.

A simplistic, yet often used, view of protein stability is that amino acids attract other amino acids with similar polarity, whereas nonpolar and polar side chains repel. Here we show that nonpolar/polar interactions, namely Val or Ile bonding to Lys or Arg in alpha-helices, can in fact be stabilizing. Residues spaced i, i + 4 in alpha-helices are on the same face of the helix, with potential to favorably interact and stabilize the structure. We observe that the nonpolar/polar pairs Ile-Lys, Ile-Arg, and Val-Lys occur in protein helices more often than expected when spaced i, i + 4. Partially helical peptides containing pairs of nonpolar/polar residues were synthesized. Controls with i, i + 5 spacing have the residues on opposite faces of the helix and are less helical than the test peptides with the i, i + 4 interactions. Experimental circular dichroism results were analyzed with helix-coil theory to calculate the free energy for the interactions. All three stabilize the helix with DeltaG between -0.14 and -0.32 kcal x mol(-1). The interactions are hydrophobic with contacts between Val or Ile and the alkyl groups in Arg or Lys. Side chains such as Lys and Arg can thus interact favorably with both polar and nonpolar residues.     

Rotamer strain energy in protein helices - quantification of a major force opposing protein folding

It is widely believed that the dominant force opposing protein folding is the entropic cost of restricting internal rotations. The energetic changes from restricting side-chain torsional motion are more complex than simply a loss of conformational entropy, however. A second force opposing protein folding arises when a side-chain in the folded state is not in its lowest-energy rotamer, giving rotameric strain. chi strain energy results from a dihedral angle being shifted from the most stable conformation of a rotamer when a protein folds. We calculated the energy of a side-chain as a function of its dihedral angles in a poly(Ala) helix. Using these energy profiles, we quantify conformational entropy, rotameric strain energy and chi strain energy for all 17 amino acid residues with side-chains in alpha-helices. We can calculate these terms for any amino acid in a helix interior in a protein, as a function of its side-chain dihedral angles, and have implemented this algorithm on a web page. The mean change in rotameric strain energy on folding is 0.42 kcal mol-1 per residue and the mean chi strain energy is 0.64 kcal mol-1 per residue. Loss of conformational entropy opposes folding by a mean of 1.1 kcal mol-1 per residue, and the mean total force opposing restricting a side-chain into a helix is 2.2 kcal mol-1. Conformational entropy estimates alone therefore greatly underestimate the forces opposing protein folding. The introduction of strain when a protein folds should not be neglected when attempting to quantify the balance of forces affecting protein stability. Consideration of rotameric strain energy may help the use of rotamer libraries in protein design and rationalise the effects of mutations where side-chain conformations change.     

Evaluating tilt angles of membrane-associated helices: comparison of computational and NMR techniques

A computational method to calculate the orientation of membrane-associated alpha-helices with respect to a lipid bilayer has been developed. It is based on a previously derived implicit membrane representation, which was parameterized using the structures of 46 alpha-helical membrane proteins. The method is validated by comparison with an independent data set of six transmembrane and nine antimicrobial peptides of known structure and orientation. The minimum energy orientations of the transmembrane helices were found to be in good agreement with tilt and rotation angles known from solid-state NMR experiments. Analysis of the free-energy landscape found two types of minima for transmembrane peptides: i), Surface-bound configurations with the helix long axis parallel to the membrane, and ii), inserted configurations with the helix spanning the membrane in a perpendicular orientation. In all cases the inserted configuration also contained the global energy minimum. Repeating the calculations with a set of solution NMR structures showed that the membrane model correctly distinguishes native transmembrane from nonnative conformers. All antimicrobial peptides investigated were found to orient parallel and bind to the membrane surface, in agreement with experimental data. In all cases insertion into the membrane entailed a significant free-energy penalty. An analysis of the contributions of the individual residue types confirmed that hydrophobic residues are the main driving force behind membrane protein insertion, whereas polar, charged, and aromatic residues were found to be important for the correct orientation of the helix inside the membrane.     

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

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

Thermodynamic model of secondary structure for α-helical peptides and proteins

A thermodynamic model describing formation of α-helices by peptides and proteins in the absence of specific tertiary interactions has been developed. The model combines free energy terms defining α-helix stability in aqueous solution and terms describing immersion of every helix or fragment of coil into a micelle or a nonpolar droplet created by the rest of protein to calculate averaged or lowest energy partitioning of the peptide chain into helical and coil fragments. The α-helix energy in water was calculated with parameters derived from peptide substitution and protein engineering data and using estimates of nonpolar contact areas between side chains. The energy of nonspecific hydrophobic interactions was estimated considering each α-helix or fragment of coil as freely floating in the spherical micelle or droplet, and using water/cyclohexane (for micelles) or adjustable (for proteins) side-chain transfer energies. The model was verified for 96 and 36 peptides studied by ¹H-nmr spectroscopy in aqueous solution and in the presence of micelles, respectively ([set I] and [set 2]) and for 30 mostly α-helical globular proteins ([set 3]). For peptides, the experimental helix locations were identified from the published medium-range nuclear Overhauser effects detected by ¹H-nmr spectroscopy. For sets 1, 2, and 3, respectively, 93, 100, and 97% of helices were identified with average errors in calculation of helix boundaries of 1.3, 2.0, and 4.1 residues per helix and an average percentage of correctly calculated helix—coil states of 93, 89, and 81%, respectively. Analysis of adjustable parameters of the model (the entropy and enthalpy of the helix—coil transition, the transfer energy of the helix backbone, and parameters of the bound coil), determined by minimization of the average helix boundary deviation for each set of peptides or proteins, demonstrates that, unlike micelles, the interior of the effective protein droplet has solubility characteristics different from that for cyclohexane, does not bind fragments of coil, and lacks interfacial area. © 1997 John Wiley & Sons, Inc. Biopoly 42: 239–269, 1997     

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.