Unravelling the folding of bacteriorhodopsin

Helical reorganization in the context of membrane protein folding: Insights from simulations with bacteriorhodopsin (BR) fragments

Membrane protein folding is distinct from folding of soluble proteins. Conformational acquisition in major membrane protein subclasses can be delineated into insertion and folding processes. An exception to the "two stage" folding, later developed to "three stage" folding, is observed within the last two helices in bacteriorhodopsin (BR), a system that serves as a model membrane protein. We employ a reductionist approach to understand interplay of molecular factors underlying the apparent defiance. Leveraging available solution NMR structures, we construct, sample in silico, and analyze partially (PIn) and fully inserted (FIn) BR membrane states. The membrane lateral C-terminal helix (CH) in PIn is markedly prone to transient structural distortions over microsecond timescales; a disorder prone region (DPR) is thereby identified. While clear transmembrane propensities are not acquired, the distortions induce alterations in local membrane curvature and area per lipid. Importantly, energetic decompositions reveal that overall, the N-terminal helix (NH) is thermodynamically more stable in the PIn. Higher overall stability of the FIn arises from favorable interactions between the NH and the CH. Our results establish lack of spontaneous transition of the PIn to the FIn, and attributes their partitioning to barriers that exceed those accessible with thermal fluctuations. This work paves the way for further detailed studies aimed at determining the thermo-kinetic roles of the initial five helices, or complementary external factors, in complete helical folding and insertion in BR. We comment that complementing such efforts with the growing field of machine learning assisted energy landscape searches may offer unprecedented insights.

Folding speeds of helical membrane proteins

Membrane proteins play key roles in human health, contributing to cellular signaling, ATP synthesis, immunity, and metabolite transport. Protein folding is the pivotal early step for their proper functioning. Understanding how this class of proteins adopts their native folds could potentially aid in drug design and therapeutic interventions for misfolding diseases. It is an essential piece in the whole puzzle to untangle their kinetic complexities, such as how rapid membrane proteins fold, how their folding speeds are influenced by changing conditions, and what mechanisms are at play. This review explores the folding speed aspect of multipass α-helical membrane proteins, encompassing plausible folding scenarios based on the timing and stability of helix packing interactions, methods for characterizing the folding time scales, relevant folding steps and caveats for interpretation, and potential implications. The review also highlights the recent estimation of the so-called folding speed limit of helical membrane proteins and discusses its consequent impact on the current picture of folding energy landscapes.

Monitoring the Progression of Cell-Free Expression of Microbial Rhodopsins by Surface Enhanced IR Spectroscopy: Resolving a Branch Point for Successful/Unsuccessful Folding

The translocon-unassisted folding process of transmembrane domains of the microbial rhodopsins sensory rhodopsin I (HsSRI) and II (HsSRII), channelrhodopsin II (CrChR2), and bacteriorhodopsin (HsBR) during cell-free expression has been investigated by Surface-Enhanced Infrared Absorption Spectroscopy (SEIRAS). Up to now, only a limited number of rhodopsins have been expressed and folded into the functional holoprotein in cell free expression systems, while other microbial rhodopsins fail to properly bind the chromophore all-trans retinal as indicated by the missing visible absorption. SEIRAS experiments suggest that all investigated rhodopsins lead to the production of polypeptides, which are co-translationally inserted into a solid-supported lipid bilayer during the first hour after the in-vitro expression is initiated. Secondary structure analysis of the IR spectra revealed that the polypeptides form a comparable amount of α-helical structure during the initial phase of insertion into the lipid bilayer. As the process progressed (>1 h), only HsBR exhibited a further increase and association of α-helices to form a compact tertiary structure, while the helical contents of the other rhodopsins stagnated. This result suggests that the molecular reason for the unsuccessful cell-free expression of the two sensory rhodopsins and of CrChR2 is not due to the translation process, but rather to the folding process during the post-translational period. Taking our previous observation into account that HsBR fails to form a tertiary structure in the absence of its retinal, we infer that the chromophore retinal is an integral component of the compaction of the polypeptide into its tertiary structure and the formation of a fully functional protein.

The Role of the Membrane in Transporter Folding and Activity

The synthesis, folding, and function of membrane transport proteins are critical factors for defining cellular physiology. Since the stability of these proteins evolved amidst the lipid bilayer, it is no surprise that we are finding that many of these membrane proteins demonstrate coupling of their structure or activity in some way to the membrane. More and more transporter structures are being determined with some information about the surrounding membrane, and computational modeling is providing further molecular details about these solvation structures. Thus, the field is moving towards identifying which molecular mechanisms - lipid interactions, membrane perturbations, differential solvation, and bulk membrane effects - are involved in linking membrane energetics to transporter stability and function. In this review, we present an overview of these mechanisms and the growing evidence that the lipid bilayer is a major determinant of the fold, form, and function of membrane transport proteins in membranes.

Detergent-free Lipodisq Nanoparticles Facilitate High-Resolution Mass Spectrometry of Folded Integral Membrane Proteins

Integral membrane proteins pose considerable challenges to mass spectrometry (MS) owing to the complexity and diversity of the components in their native environment. Here, we use native MS to study the post-translational maturation of bacteriorhodopsin (bR) and archaerhodopsin-3 (AR3), using both octyl-glucoside detergent micelles and lipid-based nanoparticles. A lower collision energy was required to obtain well-resolved spectra for proteins in styrene-maleic acid copolymer (SMA) Lipodisqs than in membrane scaffold protein (MSP) Nanodiscs. By comparing spectra of membrane proteins prepared using the different membrane mimetics, we found that SMA may favor selective solubilization of correctly folded proteins and better preserve native lipid interactions than other membrane mimetics. Our spectra reveal the correlation between the post-translation modifications (PTMs), lipid-interactions, and protein-folding states of bR, providing insights into the process of maturation of the photoreceptor proteins.

Using Fluorescence Quenching Titration to Determine the Orientation of a Model Transmembrane Protein in Mimic Membranes

After synthesis of transmembrane proteins (TMPs), they are transferred and inserted into plasma membranes to play biological functions. Crucially, orientation of TMPs in membranes determines whether they have biological activities. In cellular environments, a number of cofactors, such as translocon, can assist TMPs to be inserted into membranes in defined orientations. During in vitro reconstitution of TMPs with mimic membranes, both insertion and orientation of TMPs are primarily determined by interactions with the membrane. Yet the knowledge is limited, hindering the in vitro applications of TMPs. Here, we take Bacteriorhodopsin (bR) as a model TMP, using fluorescence quenching titration experiment to identify orientation of bR in mimic membranes, examining effects of a number of factors, including lipid composition, pH value, ionic strength and membrane curvature. The most effective determinant is the lipid type, which modulates insertion and orientation of bR in membranes by changing the membrane surface charge and the membrane fluidity. Both the pH value and the ionic strength play secondary roles by tuning the nature of the electrostatic interaction. The membrane curvature was found to have a minor effect on orientation of bR in membranes. By comparing orientations of bR in folded and unfolded states, no obvious change was observed, informing that nascent proteins could be inserted into membranes in defined orientations before folding into the native state inside the membrane.

Recent advances in biophysical studies of rhodopsins - Oligomerization, folding, and structure

Retinal-binding proteins, mainly known as rhodopsins, function as photosensors and ion transporters in a wide range of organisms. From halobacterial light-driven proton pump, bacteriorhodopsin, to bovine photoreceptor, visual rhodopsin, they have served as prototypical α-helical membrane proteins in a large number of biophysical studies and aided in the development of many cutting-edge techniques of structural biology and biospectroscopy. In the last decade, microbial and animal rhodopsin families have expanded significantly, bringing into play a number of new interesting structures and functions. In this review, we will discuss recent advances in biophysical approaches to retinal-binding proteins, primarily microbial rhodopsins, including those in optical spectroscopy, X-ray crystallography, nuclear magnetic resonance, and electron paramagnetic resonance, as applied to such fundamental biological aspects as protein oligomerization, folding, and structure.

In-Situ Observation of Membrane Protein Folding during Cell-Free Expression

Proper insertion, folding and assembly of functional proteins in biological membranes are key processes to warrant activity of a living cell. Here, we present a novel approach to trace folding and insertion of a nascent membrane protein leaving the ribosome and penetrating the bilayer. Surface Enhanced IR Absorption Spectroscopy selectively monitored insertion and folding of membrane proteins during cell-free expression in a label-free and non-invasive manner. Protein synthesis was performed in an optical cell containing a prism covered with a thin gold film with nanodiscs on top, providing an artificial lipid bilayer for folding. In a pilot experiment, the folding pathway of bacteriorhodopsin via various secondary and tertiary structures was visualized. Thus, a methodology is established with which the folding reaction of other more complex membrane proteins can be observed during protein biosynthesis (in situ and in operando) at molecular resolution.

Measurements of the effect of membrane asymmetry on the mechanical properties of lipid bilayers

We detail an approach for constructing asymmetric membranes and characterising their mechanical properties, leading to the first measurement of the effect of asymmetry on lipid bilayer mechanics. Our results demonstrate that asymmetry induces a significant increase in rigidity compared to symmetric membranes. Given that all biological membranes are asymmetric our findings have profound implications for the role of this phenomenon in biology.

Heterogeneity in Regeneration of Bacteriorhodopsin from Bacterio-Opsin and All-transRetinal at High Temperatures: Implications for Dynamic Structural Fluctuations

Measurements of regeneration kinetics were performed in order to investigate the regeneration mechanisms of bacteriorhodopsin (bR) from thermally unfolded bacterio-opsin (bO) and all-trans retinal. Regeneration kinetics data were successfully fitted to a single exponential function when regeneration was performed at 25 degrees C after incubation at high temperatures. Conversely, the process of regeneration after the addition of retinal to bO at high temperatures occurred at two different rate constants. These findings strongly suggest that the slower regeneration of bR at high temperatures occurs as a result of dynamic structural fluctuation of bO, whereas the faster process corresponds to regeneration from bO, which retains a native structure capable of retinal binding.

Rab1a and Rab5a preferentially bind to binary lipid compositions with higher stored curvature elastic energy

Rab proteins are a large family of GTP-binding proteins that regulate cellular membrane traffic and organelle identity. Rab proteins cycle between association with membranes and binding to RabGDI. Bound on membranes, each Rab has a very specific cellular location and it is this remarkable degree of specificity with which Rab GTPases recognize distinct subsets of intracellular membranes that forms the basis of their ability to act as key cellular regulators, determining the recruitment of downstream effectors to the correct membrane at the correct time. The molecular mechanisms controlling Rab localization remain poorly understood. Here, we present a fluorescence-based assay to investigate Rab GTPase membrane extraction and delivery by RabGDI. Using EGFP-Rab fusion proteins the amount of Rab:GDI complex obtained by GDI extraction of Rab proteins from HEK293 membranes could be determined, enabling control of complex concentration. Subsequent partitioning of the Rab GTPases into vesicles made up of artificial binary lipid mixtures showed for the first time, that the composition of the target membrane plays a key role in the localization of Rab proteins by sensing the stored curvature elastic energy in the membrane.

An empirical test of convergent evolution in rhodopsins

Rhodopsins are photochemically reactive membrane proteins that covalently bind retinal chromophores. Type I rhodopsins are found in both prokaryotes and eukaryotic microbes, whereas type II rhodopsins function as photoactivated G-protein coupled receptors (GPCRs) in animal vision. Both rhodopsin families share the seven transmembrane α-helix GPCR fold and a Schiff base linkage from a conserved lysine to retinal in helix G. Nevertheless, rhodopsins are widely cited as a striking example of evolutionary convergence, largely because the two families lack detectable sequence similarity and differ in many structural and mechanistic details. Convergence entails that the shared rhodopsin fold is so especially suited to photosensitive function that proteins from separate origins were selected for this architecture twice. Here we show, however, that the rhodopsin fold is not required for photosensitive activity. We engineered functional bacteriorhodopsin variants with novel folds, including radical noncircular permutations of the α-helices, circular permutations of an eight-helix construct, and retinal linkages relocated to other helices. These results contradict a key prediction of convergence and thereby provide an experimental attack on one of the most intractable problems in molecular evolution: how to establish structural homology for proteins devoid of discernible sequence similarity.

Vibrational coupling between helices influences the amide I infrared absorption of proteins: application to bacteriorhodopsin and rhodopsin

The amide I spectrum of multimers of helical protein segments was simulated using transition dipole coupling (TDC) for long-range interactions between individual amide oscillators and DFT data from dipeptides (la Cour Jansen et al. J. Chem. Phys.2006, 125, 44312) for nearest neighbor interactions. Vibrational coupling between amide groups on different helices shift the helix absorption to higher wavenumbers. This effect is small for helix dimers (1 cm(-1)) at 10 Å distance and only moderately affected by changes in the relative orientation between the helices. However, the effect becomes considerable when several helices are bundled in membrane proteins. Particular examples are the 7-helix membrane proteins bacteriorhodopsin (BR) and rhodopsin, where the upshift is 4.3 and 5.3 cm(-1), respectively, due to interhelical coupling within a BR monomer. A further upshift of 4.0 cm(-1) occurs when BR monomers associate to trimers. We propose that interhelical vibrational coupling explains the experimentally observed unusually high wavenumber of the amide I band of BR.

Engineering de novo membrane-mediated protein-protein communication networks

Mechanical properties of biological membranes are known to regulate membrane protein function. Despite this, current models of protein communication typically feature only direct protein-protein or protein-small molecule interactions. Here we show for the first time that, by harnessing nanoscale mechanical energy within biological membranes, it is possible to promote controlled communication between proteins. By coupling lipid-protein modules and matching their response to the mechanical properties of the membrane, we have shown that the action of phospholipase A(2) on acyl-based phospholipids triggers the opening of the mechanosensitive channel, MscL, by generating membrane asymmetry. Our findings confirm that the global physical properties of biological membranes can act as information pathways between proteins, a novel mechanism of membrane-mediated protein-protein communication that has important implications for (i) the underlying structure of signaling pathways, (ii) our understanding of in vivo communication networks, and (iii) the generation of building blocks for artificial protein networks.

Progress in understanding the role of lipids in membrane protein folding

Detailed investigations of membrane protein folding present a number of serious technical challenges. Most studies addressing this subject have emphasized aspects of protein amino acid sequence and structure. While it is generally accepted that the interplay between proteins and lipids plays an important role in membrane protein folding, the role(s) played by membrane lipids in this process have only recently been explored in any detail. This review is intended to summarize recent studies in which particular lipids or membrane physical properties have been shown to play a role in the folding of intact, functionally competent integral membrane proteins. This article is part of a Special Issue entitled: Protein Folding in Membranes.

Stable folding core in the folding transition state of an alpha-helical integral membrane protein

Defining the structural features of a transition state is important in understanding a folding reaction. Here, we use Φ-value and double mutant analyses to probe the folding transition state of the membrane protein bacteriorhodopsin. We focus on the final C-terminal helix, helix G, of this seven transmembrane helical protein. Φ-values could be derived for 12 amino acid residues in helix G, most of which have low or intermediate values, suggesting that native structure is disrupted at these amino acid positions in the transition state. Notably, a cluster of residues between E204 and M209 all have Φ-values close to zero. Disruption of helix G is further confirmed by a low Φ-value of 0.2 between residues T170 on helix F and S226 on helix G, suggesting the absence of a native hydrogen bond between helices F and G. Φ-values for paired mutations involved in four interhelical hydrogen bonds revealed that all but one of these bonds is absent in the transition state. The unstructured helix G contrasts with Φ-values along helix B that are generally high, implying native structure in helix B in the transition state. Thus helix B seems to constitute part of a stable folding nucleus while the consolidation of helix G is a relatively late folding event. Polarization of secondary structure correlates with sequence position, with a structured helix B near the N terminus contrasting with an unstructured C-terminal helix G.

Dual energy landscape: the functional state of the β-barrel outer membrane protein G molds its unfolding energy landscape

We applied dynamic single-molecule force spectroscopy to quantify the parameters (free energy of activation and distance of the transition state from the folded state) characterizing the energy barriers in the unfolding energy landscape of the outer membrane protein G (OmpG) from Escherichia coli. The pH-dependent functional switching of OmpG directs the protein along different regions on the unfolding energy landscape. The two functional states of OmpG take the same unfolding pathway during the sequential unfolding of β-hairpins I-IV. After the initial unfolding events, the unfolding pathways diverge. In the open state, the unfolding of β-hairpin V in one step precedes the unfolding of β-hairpin VI. In the closed state, β-hairpin V and β-strand S11 with a part of extracellular loop L6 unfold cooperatively, and subsequently β-strand S12 unfolds with the remaining loop L6. These two unfolding pathways in the open and closed states join again in the last unfolding step of β-hairpin VII. Also, the conformational change from the open to the closed state witnesses a rigidified extracellular gating loop L6. Thus, a change in the conformational state of OmpG not only bifurcates its unfolding pathways but also tunes its mechanical properties for optimum function.

Towards a structural understanding of the smallest known oncoprotein: investigation of the bovine papillomavirus E5 protein using solution-state NMR

The homodimeric E5 protein from bovine papillomavirus activates the platelet-derived growth factor β receptor through transmembrane (TM) helix-helix interactions leading to uncontrolled cell growth. Detailed structural information for the E5 dimer is essential if we are to uncover its unique mechanism of action. In vivo mutagenesis has been used to identify residues in the TM domain critical for dimerization, and we previously reported that a truncated synthetic E5 peptide forms dimers via TM domain interactions. Here we extend this work with the first application of high-resolution solution-state NMR to the study of the E5 TM domain in SDS micelles. Using selectively 15N-labelled peptides, we first probe sample homogeneity revealing two predominate species, which we interpret to be monomer and dimer. The equilibrium between the two states is shown to be dependent on detergent concentration, revealed by intensity shifts between two sets of peaks in 15N-(1)H HSQC experiments, highlighting the importance of sample preparation when working with these types of proteins. This information is used to estimate a free energy of association (ΔGx°=-3.05 kcal mol(-1)) for the dimerization of E5 in SDS micelles. In addition, chemical shift changes have been observed that indicate a more pronounced change in chemical environment for those residues expected to be at the dimer interface in vivo versus those that are not. Thus we are able to demonstrate our in vitro dimer is comparable to that defined in vivo, validating the biological significance of our synthetic peptide and providing a solid foundation upon which to base further structural studies. Using detergent concentration to modulate oligomeric state and map interfacial residues by NMR could prove useful in the study of other homo-oligomeric transmembrane proteins.

In vitro dimerization of the bovine papillomavirus E5 protein transmembrane domain

The E5 protein from bovine papillomavirus is a type II membrane protein and the product of the smallest known oncogene. E5 causes cell transformation by binding and activating the platelet-derived growth factor beta receptor (PDGFbetaR). In order to productively interact with the receptor, it is thought that E5 binds as a dimer. However, wild-type E5 and various mutants have also been shown to form trimers, tetramers, and even higher order oligomers. The residues in E5 that drive and stabilize a dimeric state are also still in question. At present, two different models for the E5 dimer exist in the literature, one symmetric and one asymmetric. There is universal agreement, however, that the transmembrane (TM) domain plays a vital role in stabilizing the functional oligomer; indeed, mutation of various TM domain residues can abolish E5 function. In order to better resolve the role of the E5 TM domain in function, we have undertaken the first quantitative in vitro characterization of the E5 TM domain in detergent micelles and liposomes. Circular and linear dichroism analyses verify that the TM domain adopts a stable alpha-helical structure and is able to partition efficiently across lipid bilayers. SDS-PAGE and analytical ultracentrifugation demonstrate for the first time that the TM domain of E5 forms a strong dimer with a standard state free energy of dissociation of 5.0 kcal mol (-1). We have used our new results to interpret existing models of E5 dimer formation and provide a direct link between TM helix interactions and E5 function.

Combined kinetic and thermodynamic analysis of alpha-helical membrane protein unfolding

The analytical toolkit developed for investigations into water-soluble protein folding has yet to be applied in earnest to membrane proteins. A major problem is the difficulty in collecting kinetic data, which are crucial to understanding any reaction. Here, we combine kinetic and thermodynamic studies of the reversible unfolding of an alpha-helical membrane protein to provide a definitive value for the reaction free energy and a means to probe the transition state. Our analyses show that the major unfolding step in the SDS-induced denaturation of bacteriorhodopsin involves a reduction in alpha-helical structure and proceeds with a large free-energy change; both our equilibrium and kinetic measurements predict that the free energy of unfolding in the absence of denaturant is +20 kcal.mol(-1), with an associated m-value of 25 kcal.mol(-1). The rate of unfolding in the absence of denaturant, k(u)(H(2)O), is surprisingly very slow ( approximately 10(-15) s(-1)). The kinetics also give information on the transition state for this major unfolding step, with a value for beta (m(f)/[m(f) + m(u)]) of approximately 0.1, indicating that the transition state is close to the unfolded state. We thus present a basis for mapping the structural and energetic properties of membrane protein folding by mutagenesis and classical kinetics.

Deciphering molecular interactions of native membrane proteins by single-molecule force spectroscopy

Molecular interactions are the basic language of biological processes. They establish the forces interacting between the building blocks of proteins and other macromolecules, thus determining their functional roles. Because molecular interactions trigger virtually every biological process, approaches to decipher their language are needed. Single-molecule force spectroscopy (SMFS) has been used to detect and characterize different types of molecular interactions that occur between and within native membrane proteins. The first experiments detected and localized molecular interactions that stabilized membrane proteins, including how these interactions were established during folding of alpha-helical secondary structure elements into the native protein and how they changed with oligomerization, temperature, and mutations. SMFS also enables investigators to detect and locate molecular interactions established during ligand and inhibitor binding. These exciting applications provide opportunities for studying the molecular forces of life. Further developments will elucidate the origins of molecular interactions encoded in their lifetimes, interaction ranges, interplay, and dynamics characteristic of biological systems.

Assembly of the Cytochrome bo3 Complex

An understanding of the mechanisms that govern the assembly of macromolecular protein complexes is fundamental for studying their function and regulation. With this in mind, we have determined the assembly pathway for the membrane-embedded cytochrome bo(3) of Escherichia coli. We show that there is a preferred order of assembly, where subunits III and IV assemble first, followed by subunit I and finally subunit II. We also show that cofactor insertion catalyses assembly. These findings provide novel insights into the biogenesis of this model membrane protein complex.

An unfolding story of helical transmembrane proteins

Reversible unfolding of helical transmembrane proteins could provide valuable information about the free energy of interaction between transmembrane helices. Thermal unfolding experiments suggest that this process for integral membrane proteins is irreversible. Chemical unfolding has been accomplished with organic acids, but the unfolding or refolding pathways involve irreversible steps. Sodium dodecyl sulfate (SDS) has been used as a perturbant to study reversible unfolding and refolding kinetics. However, the interpretation of these experiments is not straightforward. It is shown that the results could be explained by SDS binding without substantial unfolding. Furthermore, the SDS-perturbed state is unlikely to include all of the entropy terms involved in an unfolding process. Alternative directions for future research are suggested: fluorinated alcohols in homogeneous solvent systems, inverse micelles, and fragment association studies.

Membrane proteins shape up: understanding in vitro folding

Progress in structural biology has begun to reveal the precise architecture of integral membrane proteins. However, the manner in which these complex structures are achieved remains unclear. Recent developments are starting to shed light on the unfolding and folding of a small but growing number of membrane proteins. Mechanistic details derived from kinetic and thermodynamic experiments now enable comparison of the folding of different membrane proteins and their water-soluble cousins. This work also has important implications for other structural and functional studies of membrane proteins in vitro.

Transmembrane helix-helix association: relative stabilities at low pH

We have previously studied the unfolding equilibrium of bacterioopsin in a single phase solvent, using Förster mechanism fluorescence resonance energy transfer (FRET) as a probe, from tryptophan donors to a dansyl acceptor. We observed an apparent unfolding transition in bacterioopsin perturbed by increasing ethanol concentrations [Nannepaga et al. (2004) Biochemistry 43, 50-59]. We have further investigated this transition and find that the unfolding is pH-dependent. We have now measured the apparent pK of acid-induced unfolding of bacterioopsin in 90% ethanol. When the acceptor is on helix B (Lys 41), the apparent pK for unfolding is 4.75; on the EF connecting loop (Cys 163), 5.15; and on helix G (Cys 222), 5.75. Five-helix proteolytic fragments are less stable. The apparent unfolding pKs are 5.46 for residues 72-248 (Cys 163) and 7.36 for residues 1-166 (Lys 41). When interpreted in terms of a simple equilibrium model for unfolding, the apparent pKs give relative free energies of unfolding in the range of -0.54 to -3.5 kcal/mol. The results suggest that the C-terminal helix of bacterioopsin is less stably folded than the N-terminal helices. We analyzed the pairwise helix-helix interaction surfaces of bacteriorhodopsin and three other seven-transmembrane-helix proteins on the basis of crystal structures. The results show that the interaction surfaces are smoother and the helix axis separations are closer in the amino-terminal two-thirds of the proteins compared with the carboxyl-terminal one-third. However, the F helix is important in stabilizing the folded structure, as shown by the instability of the 1-166 fragment. Considering the high-resolution crystal structure of bacteriorhodopsin, there are no obvious helix-helix interactions involving protein side chains which would be destabilized by protonation at the estimated pH of the unfolding transitions. However, a number of helix-bridging water molecules could become protonated, thereby weakening the helix-helix interactions.

Characterizing folding, structure, molecular interactions and ligand gated activation of single sodium/proton antiporters

Using the example of sodium/proton antiporter from Escherichia coli NhaA, we review the capabilities of single-molecule atomic force microscopy and force spectroscopy to observe structural and functional insights of a membrane protein, which are not attainable by other traditional methods. While atomic force microscopy provides high-resolution topographs of single membrane proteins, their oligomeric state and assembly, single-molecule force spectroscopy experiments detect molecular interactions of the protein. The sensitivity of this method makes it possible to detect and locate interactions that stabilize secondary structures such as transmembrane alpha-helices, polypeptide loops and segments within them. Controlled refolding experiments using single-molecule force spectroscopy observed individual secondary structure segments folding into the functional protein. Various folding pathways of NhaA were detected, each one exhibiting a certain probability to be taken. Time-lapse refolding experiments enabled determining the folding kinetics and hierarchy of individual secondary structural elements. Recent examples detected and located the ligand binding of an antiporter. Similarly, inhibitor binding and location can be detected which in future may guide towards comparative studies of agonist and antagonist altering the functional state of a membrane protein. We review current and future potentials of these approaches to characterize the action of pharmacological molecules on the antiporter function.

Bacteriorhodopsin Folds into the Membrane against an External Force

Despite their crucial importance for cellular function, little is known about the folding mechanisms of membrane proteins. Recently details of the folding energy landscape were elucidated by atomic force microscope (AFM)-based single molecule force spectroscopy. Upon unfolding and extraction of individual membrane proteins energy barriers in structural elements such as loops and helices were mapped and quantified with the precision of a few amino acids. Here we report on the next logical step: controlled refolding of single proteins into the membrane. First individual bacteriorhodopsin monomers were partially unfolded and extracted from the purple membrane by pulling at the C-terminal end with an AFM tip. Then by gradually lowering the tip, the protein was allowed to refold into the membrane while the folding force was recorded. We discovered that upon refolding certain helices are pulled into the membrane against a sizable external force of several tens of picoNewton. From the mechanical work, which the helix performs on the AFM cantilever, we derive an upper limit for the Gibbs free folding energy. Subsequent unfolding allowed us to analyze the pattern of unfolding barriers and corroborate that the protein had refolded into the native state.

Kinetics of an individual transmembrane helix during bacteriorhodopsin folding

The kinetics of an individual helix of bacteriorhodopsin have been monitored during folding of the protein into lipid bilayer vesicles. A fluorescence probe was introduced at individual sites throughout helix D of bacteriorhodopsin and the changes in the fluorescence of the label were time-resolved. Partially denatured, labelled bacteriorhodopsin in SDS was folded directly into phosphatidylcholine lipid vesicles. Stopped-flow mixing of the reactants allowed the folding kinetics to be monitored with millisecond time resolution by time-resolving changes in the label fluorescence, intrinsic protein fluorescence as well as in the absorption of the retinal chromophore. Monitoring specific positions on helix D showed that two kinetic phases were altered compared to those determined by monitoring the average protein behaviour. These two phases, of 6.7 s(-1) and 0.33 s(-1), were previously assigned to formation of a key apoprotein intermediate during bacteriorhodopsin folding. The faster 6.7s(-1) phase was missing when time-resolving fluorescence changes of labels attached to the middle of helix D. The amplitude of the 0.33 s(-1) phase increased along the helix, as single labels were attached in turn from the cytoplasmic to the extracellular side. An interpretation of these results is that the 6.7 s(-1) phase involves partitioning of helix D within the lipid headgroups of the bilayer vesicle, while the 0.33 s(-1) phase could reflect transmembrane insertion of this helix. In addition, a single site on helix G was monitored during folding. The results indicate that, unlike helix D, the insertion of helix G cannot be differentiated from the average protein behaviour. The data show that, while folding of bacteriorhodopsin from SDS into lipids is a co-operative process, it is nevertheless possible to obtain information on specific regions of a membrane protein during folding in vitro.

Observing Folding Pathways and Kinetics of a Single Sodium-proton Antiporter from Escherichia coli

Mechanisms of folding and misfolding of membrane proteins are of interest in cell biology. Recently, we have established single-molecule force spectroscopy to observe directly the stepwise folding of the Na+/H+ antiporter NhaA from Escherichia coli in vitro. Here, we improved this approach significantly to track the folding intermediates of a single NhaA polypeptide forming structural segments such as the Na+-binding site, transmembrane alpha-helices, and helical pairs. The folding rates of structural segments ranged from 0.31 s(-1) to 47 s(-1), providing detailed insight into a distinct folding hierarchy of an unfolded polypeptide into the native membrane protein structure. In some cases, however, the folding chain formed stable and kinetically trapped non-native structures, which could be assigned to misfolding events of the antiporter.

Defining the Structural Basis for Assembly of a Transmembrane Cytochrome

To define the structural basis for cofactor binding to membrane proteins, we introduce a manageable model system, which allows us, for the first time, to study the influence of individual transmembrane helices and of single amino acid residues on the assembly of a transmembrane cytochrome. In vivo as well as in vitro analyses indicate central roles of single amino acid residues for either interaction of the transmembrane helices or for binding of the cofactor. The results clearly show that interaction of the PsbF transmembrane helix is independent from binding of the heme cofactor. On the other hand, binding of the cofactor highly depends on helix-helix interactions. By site-directed mutagenesis critical amino acid residues were identified, which are involved in the assembly of a functional transmembrane cytochrome. Especially, a highly conserved glycine residue is critical for interaction of the transmembrane helices and assembly of the cytochrome. Based on the two-stage-model of alpha-helical membrane protein folding, the presented results clearly indicate a third stage of membrane protein folding, in which a cofactor binds to a pre-assembled transmembrane protein.

Transmembrane helix prediction: a comparative evaluation and analysis

The prediction of transmembrane (TM) helices plays an important role in the study of membrane proteins, given the relatively small number (approximately 0.5% of the PDB) of high-resolution structures for such proteins. We used two datasets (one redundant and one non-redundant) of high-resolution structures of membrane proteins to evaluate and analyse TM helix prediction. The redundant (non-redundant) dataset contains structure of 434 (268) TM helices, from 112 (73) polypeptide chains. Of the 434 helices in the dataset, 20 may be classified as 'half-TM' as they are too short to span a lipid bilayer. We compared 13 TM helix prediction methods, evaluating each method using per segment, per residue and termini scores. Four methods consistently performed well: SPLIT4, TMHMM2, HMMTOP2 and TMAP. However, even the best methods were in error by, on average, about two turns of helix at the TM helix termini. The best and worst case predictions for individual proteins were analysed. In particular, the performance of the various methods and of a consensus prediction method, were compared for a number of proteins (e.g. SecY, ClC, KvAP) containing half-TM helices. The difficulties of predicting half-TM helices suggests that current prediction methods successfully embody the two-state model of membrane protein folding, but do not accommodate a third stage in which, e.g., short helices and re-entrant loops fold within a bundle of stable TM helices.

Insertion kinetics of a denatured alpha helical membrane protein into phospholipid bilayer vesicles

Membrane protein folding has suffered from a lack of detailed kinetic studies, particularly with regard to the insertion of denatured protein into lipid bilayers. We present a detailed in vitro kinetic study of the association of a denatured, transmembrane alpha helical protein with lipid vesicles. The mechanism of folding of Escherichia coli diacylglycerol kinase from a partially denatured state in urea has been investigated. The protein associates with lipid vesicles to give a protein, vesicle complex with an apparent association constant of 2 x 10(6) M(-1) s(-1). This association rate approaches the diffusion limit of the protein, vesicle reaction. The association of the protein with lipid vesicles is followed by a slower process occurring at observed rate of 0.031 s(-1), involving insertion into the bilayer and generation of a functional oligomer of diacylglycerol kinase. Protein aggregation competes with vesicle insertion. The urea-denatured protein monomers begin to aggregate as soon as the urea is diluted. This aggregation is faster than the association of the protein with vesicles so that most protein aggregates before it inserts into a vesicle. Increasing the vesicle concentration favours insertion of protein monomers, but at high vesicle concentrations monomers are primarily in separate vesicles and do not associate to form functional oligomers. Irreversible aggregation limits the yield of functional protein, while the data also suggest that lipid vesicles can reverse another aggregation reaction, leading to the recovery of correctly folded protein.

Circular dichroism of heterochromophoric and partially regenerated purple membrane: Search for exciton coupling

In order to determine the origin of the bisignate CD spectra of native purple membrane, heterochromophoric analogues containing bacteriorhodopsin regenerated with native all-trans-retinal and retinal analogues were investigated. The data collected for the purple membrane samples containing two different chromophores suggest the additive character of the CD spectra. This conclusion was supported by a series of spectra using 5,6-dihydroretinal and 3-dehydroretinal and by using 33% regenerated PM in buffer and in presence of osmolytes. Our results support the idea of conformational heterogeneity of the chromophores in the bR in the trimer, suggesting that the three bR subunits in the trimer are not conformationally equal, and therefore, the bisignate CD spectrum of bR in the purple membrane occurs rather due to a superposition of the CD spectra from variously distorted bR subunits in the trimer than interchromophoric exciton-coupling interactions.

Controlling the folding efficiency of an integral membrane protein

Research into the folding mechanisms of integral membrane proteins lags far behind that of water-soluble proteins, to the extent that the term protein folding is synonymous with water-soluble proteins. Hydrophobic membrane proteins, and particularly those with transmembrane alpha-helical motifs, are frequently considered too difficult to work with. We show that the stored curvature elastic stress of lipid bilayers can be used to guide the design of efficient folding systems for these integral membrane proteins. The curvature elastic stress of synthetic phosphatidylcholine/phosphatidylethanolamine lipid bilayers can be used to control both the rate of folding and the yield of folded protein. The use of a physical bilayer property generalises this approach beyond the particular chemistry of the lipids involved.

Snorkeling Preferences Foster an Amino Acid Composition Bias in Transmembrane Helices

By analyzing transmembrane (TM) helices in known structures, we find that some polar amino acids are more frequent at the N terminus than at the C terminus. We propose the asymmetry occurs because most polar amino acids are better able to snorkel their polar atoms away from the membrane core at the N terminus than at the C terminus. Two findings lead us to this proposition: (1) side-chain conformations are influenced strongly by the N or C-terminal position of the amino acid in the bilayer, and (2) the favored snorkeling direction of an amino acid correlates well with its N to C-terminal composition bias. Our results suggest that TM helix predictions should incorporate an N to C-terminal composition bias, that rotamer preferences of TM side-chains are position-dependent, and that the ability to snorkel influences the evolutionary selection of amino acids for the helix N and C termini.

Structure and mechanism in prokaryotic mechanosensitive channels

Mechanosensitive channels function as electromechanical switches with the capability to sense the physical state of lipid bilayers. The X-ray crystal structures of MscL and MscS offer a unique opportunity to identify the types of protein motions associated with the opening and closing of these structurally unrelated channels, while providing the framework to address a mechanism of tension sensing that is defined by channel-lipid interactions. Recent functional, structural and dynamic data offer fresh insights into the molecular basis of gating for these membrane proteins.

A Monte Carlo study of the self-assembly of bacteriorhodopsin

Bacteriorhodopsin (BR) and specific lipid molecules self-assemble into a quasi two-dimensional lattice structure known as the purple membrane (PM). In the PM, BR molecules exist in a trimeric form with lipid molecules present in the space enclosed by each trimeric unit and in the inter-trimer space. These trimeric units, which have a roughly circular cross-section, are arranged in hexagonal patterns with long-ranged crystalline order. In this work, we investigate the self-assembly of BR in the PM via Monte Carlo simulations of a two-dimensional model of the membrane and proteins. The protein molecules are modeled as 120 degrees sectors of a circle and the lipid molecules enter into the model through effective protein-protein interactions. The sectors cannot overlap with each other, and in addition to this excluded volume interaction there are site-site attractive interactions between specific points of the proteins to mimic interactions between helices on the proteins and lipid-induced interactions. At low values of the attractive well depth, the proteins are found in the monomeric form at all concentrations. At moderate and high values of the attractive well depth, trimers are formed as the concentration increases, and with a further increase in concentration the trimers organize into a hexagonal lattice. The interactions between the proteins and those induced by the intra-trimer lipids play an equally important role in the formation of trimers and the lattice. The lipids in the inter-trimer space cause the trimers to orient in a specific direction in the hexagonal crystal lattice.

Membrane proteins: shaping up

Over recent years, much progress has been made in the identification and characterization of factors involved in the biosynthesis of integral membrane proteins of the helix-bundle type. In addition, our knowledge of membrane protein structure and the forces stabilizing helix-helix interactions in a lipid environment is expanding rapidly. However, it is still not clear how a membrane protein folds into its final form in vivo, nor what constraints there are on the folded structure that results from the mechanistic details of translocon-mediated assembly rather than simply from the thermodynamics of protein-lipid interactions.

Effects of cholesterol and model transmembrane proteins on drug partitioning into lipid bilayers as analysed by immobilized-liposome chromatography

We have analysed how cholesterol and transmembrane proteins in phospholipid bilayers modulate drug partitioning into the bilayers. For this purpose we determined the chromatographic retention of drugs on liposomes or proteoliposomes entrapped in gel beads. The drug retention per phospholipid amount (the capacity factor Ks) reflects the drug partitioning. Cholesterol in the bilayers decreased the Ks value and hence the partitioning into the membrane in proportion to the cholesterol fraction. On average this cholesterol effect decreased with increasing temperature. Model transmembrane proteins, the glucose transporter GLUT1 and bacteriorhodopsin, interacted electrostatically with charged drugs to increase or decrease the drug partitioning into the bilayers. Bacteriorhodopsin proteoliposomes containing cholesterol combined the effects of the protein and the cholesterol and approached the partitioning properties of red blood cell membranes. For positively charged drugs the correlation between calculated intestinal permeability and log Ks was fair for both liposomes and bacteriorhodopsin-cholesterol proteoliposomes. Detailed modeling of solute partitioning into biological membranes may require an extensive knowledge of their structures.