Side-chain Contributions to Membrane Protein Structure and Stability

Lipid exposure prediction enhances the inference of rotational angles of transmembrane helices

Background

Since membrane protein structures are challenging to crystallize, computational approaches are essential for elucidating the sequence-to-structure relationships. Structural modeling of membrane proteins requires a multidimensional approach, and one critical geometric parameter is the rotational angle of transmembrane helices. Rotational angles of transmembrane helices are characterized by their folded structures and could be inferred by the hydrophobic moment; however, the folding mechanism of membrane proteins is not yet fully understood. The rotational angle of a transmembrane helix is related to the exposed surface of a transmembrane helix, since lipid exposure gives the degree of accessibility of each residue in lipid environment. To the best of our knowledge, there have been few advances in investigating whether an environment descriptor of lipid exposure could infer a geometric parameter of rotational angle.

Results

Here, we present an analysis of the relationship between rotational angles and lipid exposure and a support-vector-machine method, called TMexpo, for predicting both structural features from sequences. First, we observed from the development set of 89 protein chains that the lipid exposure, i.e., the relative accessible surface area (rASA) of residues in the lipid environment, generated from high-resolution protein structures could infer the rotational angles with a mean absolute angular error (MAAE) of 46.32˚. More importantly, the predicted rASA from TMexpo achieved an MAAE of 51.05˚, which is better than 71.47˚ obtained by the best of the compared hydrophobicity scales. Lastly, TMexpo outperformed the compared methods in rASA prediction on the independent test set of 21 protein chains and achieved an overall Matthew’s correlation coefficient, accuracy, sensitivity, specificity, and precision of 0.51, 75.26%, 81.30%, 69.15%, and 72.73%, respectively. TMexpo is publicly available at http://bio-cluster.iis.sinica.edu.tw/TMexpo.

Conclusions

TMexpo can better predict rASA and rotational angles than the compared methods. When rotational angles can be accurately predicted, free modeling of transmembrane protein structures in turn may benefit from a reduced complexity in ensembles with a significantly less number of packing arrangements. Furthermore, sequence-based prediction of both rotational angle and lipid exposure can provide essential information when high-resolution structures are unavailable and contribute to experimental design to elucidate transmembrane protein functions.

Reversible folding of human peripheral myelin protein 22, a tetraspan membrane protein

Misfolding of the α-helical membrane protein peripheral myelin protein 22 (PMP22) has been implicated in the pathogenesis of the common neurodegenerative disease known as Charcot-Marie-Tooth disease (CMTD) and also several other related peripheral neuropathies. Emerging evidence suggests that the propensity of PMP22 to misfold in the cell may be due to an intrinsic lack of conformational stability. Therefore, quantitative studies of the conformational equilibrium of PMP22 are needed to gain insight into the molecular basis of CMTD. In this work, we have investigated the folding and unfolding of wild type (WT) human PMP22 in mixed micelles. Both kinetic and thermodynamic measurements demonstrate that the denaturation of PMP22 by n-lauroyl sarcosine (LS) in dodecylphosphocholine (DPC) micelles is reversible. Assessment of the conformational equilibrium indicates that a significant fraction of unfolded PMP22 persists even in the absence of the denaturing detergent. However, we find the stability of PMP22 is increased by glycerol, which facilitates quantitation of thermodynamic parameters. To our knowledge, this work represents the first report of reversible unfolding of a eukaryotic multispan membrane protein. The results indicate that WT PMP22 possesses minimal conformational stability in micelles, which parallels its poor folding efficiency in the endoplasmic reticulum. Folding equilibrium measurements for PMP22 in micelles may provide an approach to assess the effects of cellular metabolites or potential therapeutic agents on its stability. Furthermore, these results pave the way for future investigation of the effects of pathogenic mutations on the conformational equilibrium of PMP22.

Stabilizing membrane proteins through protein engineering

Integral membrane proteins (IMPs) are crucial components of all cells but are difficult to study in vitro because they are generally unstable when removed from their native membranes using detergents. Despite the major biomedical relevance of IMPs, less than 1% of Protein Data Bank (PDB) entries are IMP structures, reflecting the technical gap between studies of soluble proteins compared to IMPs. Stability can be engineered into IMPs by inserting stabilizing mutations, thereby generating proteins that can be successfully applied to biochemical and structural studies when solubilized in detergent micelles. The identification of stabilizing mutations is not trivial, and this review will focus on the methods that have been used to identify stabilized membrane proteins, including alanine scanning and screening, directed evolution and computational design.

Membrane protein native state discrimination by implicit membrane models

Four implicit membrane models [IMM1, generalized Born (GB)-surface area-implicit membrane (GBSAIM), GB with a simple switching (GBSW), and heterogeneous dielectric GB (HDGB)] were tested for their ability to discriminate the native conformation of five membrane proteins from 450 decoys generated by the Rosetta-Membrane program. The energy ranking of the native state and Z-scores were used to assess the performance of the models. The effect of membrane thickness was examined and was found to be substantial. Quite satisfactory discrimination was achieved with the all-atom IMM1 and GBSW models at 25.4 Å thickness and with the HDGB model at 28.5 Å thickness. The energy components by themselves were not discriminative. Both van der Waals and electrostatic interactions contributed to native state discrimination, to a different extent in each model. Computational efficiency of the models decreased in the order: extended-atom IMM1 > all-atom IMM1 > GBSAIM > GBSW > HDGB. These results encourage the further development and use of implicit membrane models for membrane protein structure prediction.

Exploring the active site structure of a photoreceptor protein by Raman optical activity

We have developed a near-infrared excited Raman optical activity (ROA) spectrometer and report the first measurement of near-infrared ROA spectra of a light-driven proton pump, bacteriorhodopsin. Our results demonstrate that a near-infrared excitation enables us to measure the ROA spectra of the chromophore within a protein environment. Furthermore, the ROA spectra of the all-trans, 15-anti and 13-cis, 15-syn isomers differ significantly, indicating a high structural sensitivity of the ROA spectra. We therefore expect that future applications of the near-infrared ROA will allow the experimental elucidation of the active site structures in other proteins as well as reaction intermediates.

Prediction and design of outer membrane protein-protein interactions

Protein-protein interactions (PPI) play central roles in biological processes, motivating us to understand the structural basis underlying affinity and specificity. In this chapter, we focus on biochemical and computational design strategies of assessing and detecting PPIs of β-barrel outer membrane proteins (OMPs). A few case studies are presented highlighting biochemical techniques used to dissect the energetics of oligomerization and determine amino acids forming the key interactions of the PPI sites. Current computational strategies for detecting/predicting PPIs are introduced, and examples of computational and rational engineering strategies applied to OMPs are presented.

Shifting hydrogen bonds may produce flexible transmembrane helices

The intricate functions of membrane proteins would not be possible without bends or breaks that are remarkably common in transmembrane helices. The frequent helix distortions are nevertheless surprising because backbone hydrogen bonds should be strong in an apolar membrane, potentially rigidifying helices. It is therefore mysterious how distortions can be generated by the evolutionary currency of random point mutations. Here we show that we can engineer a transition between distinct distorted helix conformations in bacteriorhodopsin with a single-point mutation. Moreover, we estimate the energetic cost of the conformational transitions to be smaller than 1 kcal/mol. We propose that the low energy of distortion is explained in part by the shifting of backbone hydrogen bonding partners. Consistent with this view, extensive backbone hydrogen bond shifts occur during helix conformational changes that accompany functional cycles. Our results explain how evolution has been able to liberally exploit transmembrane helix bending for the optimization of membrane protein structure, function, and dynamics.

Identifying and measuring transmembrane helix-helix interactions by FRET

Specific interactions between helical transmembrane domains (TMs) play essential roles in the mechanisms governing the folding, stability and assembly of integral membrane proteins. Thus, it is appealing to identify helix-helix contacts and to seek the structural determinants of such interactions at the molecular level. Here, we provide a protocol for detecting and measuring specific helix-helix interactions in liposomes by Förster resonance energy transfer (FRET), using peptides corresponding to the TM domains of an integral membrane protein. We give a detailed procedure and practical guidelines on how to design, prepare, handle, and characterize fluorescently labeled TM peptides reconstituted in large unilamellar lipid vesicles. We also discuss some critical aspects of FRET measurements to ensure the correct analysis and interpretation of spectral data. Our method uses tryptophan/pyrene as the donor-acceptor FRET pair, but it can be easily adapted to other fluorescence pairs and to other membrane mimetic environments. The ability to identify crucial interhelical contacts is a valuable tool for the study of the stability, assembly, and function of the important and experimentally challenging helical membrane proteins.

Revisiting the folding kinetics of bacteriorhodopsin

The elucidation of the physical principles that govern the folding and stability of membrane proteins is one of the greatest challenges in protein science. Several insights into the folding of α-helical membrane proteins have come from the investigation of the conformational equilibrium of H. halobium bacteriorhodopsin (bR) in mixed micelles using SDS as a denaturant. In an effort to confirm that folded bR and SDS-denatured bR reach the same conformational equilibrium, we found that bR folding is significantly slower than has been previously known. Interrogation of the effect of the experimental variables on folding kinetics reveals that the rate of folding is dependent not only on the mole fraction of SDS but also on the molar concentrations of mixed micelle components, a variable that was not controlled in the previous study of bR folding kinetics. Moreover, when the molar concentrations of mixed micelle components are fixed at the concentrations commonly employed for bR equilibrium studies, conformational relaxation in the transition zone is slower than hydrolysis of the retinal Schiff base. As a result, the conformational equilibrium between folded bR and SDS-denatured bR cannot be achieved under the conventional condition. Our finding suggests that the molar concentrations of mixed micelle components are important experimental variables in the investigation of the kinetics and thermodynamics of bR folding and should be accounted for to ensure the accurate assessment of the conformational equilibrium of bR without the interference of retinal hydrolysis.

Computational studies of membrane proteins: models and predictions for biological understanding

We discuss recent progresses in computational studies of membrane proteins based on physical models with parameters derived from bioinformatics analysis. We describe computational identification of membrane proteins and prediction of their topology from sequence, discovery of sequence and spatial motifs, and implications of these discoveries. The detection of evolutionary signal for understanding the substitution pattern of residues in the TM segments and for sequence alignment is also discussed. We further discuss empirical potential functions for energetics of inserting residues in the TM domain, for interactions between TM helices or strands, and their applications in predicting lipid-facing surfaces of the TM domain. Recent progresses in structure predictions of membrane proteins are also reviewed, with further discussions on calculation of ensemble properties such as melting temperature based on simplified state space model. Additional topics include prediction of oligomerization state of membrane proteins, identification of the interfaces for protein-protein interactions, and design of membrane proteins. This article is part of a Special Issue entitled: Protein Folding in Membranes.

Advances in the mass spectrometry of membrane proteins: from individual proteins to intact complexes

Rapid advances in structural genomics and in large-scale proteomic projects have yielded vast amounts of data on soluble proteins and their complexes. Despite these advances, progress in studying membrane proteins using mass spectrometry (MS) has been slow. This is due in part to the inherent solubility and dynamic properties of these proteins, but also to their low abundance and the absence of polar side chains in amino acid residues. Considerable progress in overcoming these challenges is, however, now being made for all levels of structural characterization. This progress includes MS studies of the primary structure of membrane proteins, wherein sophisticated enrichment and trapping procedures are allowing multiple posttranslational modifications to be defined through to the secondary structure level in which proteins and peptides have been probed using hydrogen exchange, covalent, or radiolytic labeling methods. Exciting possibilities now exist to go beyond primary and secondary structure to reveal the tertiary and quaternary interactions of soluble and membrane subunits within intact assemblies of more than 700 kDa.

Interaction of a two-transmembrane-helix peptide with lipid bilayers and dodecyl sulfate micelles

To probe structural changes that occur when a membrane protein is transferred from lipid bilayers to SDS micelles, a fragment of bacteriorhodopsin containing transmembrane helical segments A and B was studied by fluorescence spectroscopy, molecular dynamics (MD) simulation, and stopped flow kinetics. In lipid bilayers, Förster resonance energy transfer (FRET) was observed between tyrosine 57 on helix B and tryptophans 10 and 12 on helix A. FRET efficiency decreased substantially when the peptide was transferred to SDS. MD simulation showed no evidence for significant disruption of helix-helix interactions in SDS micelles. However, a cluster of water molecules was observed to form a hydrogen-bonded network with the phenolic hydroxyl group of tyrosine 57, which probably causes the disappearance of tyrosine-to-tryptophan FRET in SDS. The tryptophan quantum yield decreased in SDS, and the change occurred at nearly the same rate as membrane solubilization. The results provide a clear example of the importance of corroborating distance changes inferred from FRET by using complementary methods.

Overcoming hysteresis to attain reversible equilibrium folding for outer membrane phospholipase A in phospholipid bilayers

The free energy of unfolding of a membrane protein from lipids into water (ΔG(o)(w,l)) describes its equilibrium thermodynamic stability. Knowing this parameter gives insight into a membrane protein's sequence-structure-energy relationships. However, there are few measures of membrane protein stability because of the technical difficulties associated with unfolded and partially folded states. Here, we describe the experimental process that allowed us to measure the ΔG(o)(w,l) of the outer membrane phospholipase A into large unilamellar vesicles (LUVs) of 1,2-dilauroyl-sn-glycero-3-phosphocholine. To arrive at this reversible folding condition, we screened a large number of experimental variables: temperature, incubation time, salt concentration, pH, lipid composition and liposome morphology. The principal challenge we encountered under most conditions was hysteresis between folding and unfolding titrations. A second factor that compromised reversible folding was the observation that a fraction of the protein population tended to aggregate. We found that hysteresis could be completely eliminated on a feasible timescale by conducting experiments at acidic pH, by the slow dilution of the protein in the initial titration setup and by utilizing a low concentration of a detergent as a temporary "holdase" to solubilize the protein upon its initial dilution into folding conditions. We confirmed that the detergent did not disrupt the LUVs using fluorescence emission of lipid-sensitive dyes and light scattering. The results of our parameter search should be generally useful for efforts to measure ΔG(o)(w,l) for other membrane proteins.

Thermodynamic stability of bacteriorhodopsin mutants measured relative to the bacterioopsin unfolded state

The stability of bacteriorhodopsin (bR) has often been assessed using SDS unfolding assays that monitor the transition of folded bR (bR(f)) to unfolded (bR(u)). While many criteria suggest that the unfolding curves reflect thermodynamic stability, slow retinal (RET) hydrolysis during refolding makes it impossible to perform the most rigorous test for equilibrium, i.e., superimposable unfolding and refolding curves. Here we made a new equilibrium test by asking whether the refolding rate in the transition zone is faster than RET hydrolysis. We find that under conditions we have used previously, refolding is in fact slower than hydrolysis, strongly suggesting that equilibrium is not achieved. Instead, the apparent free energy values reported previously are dominated by unfolding rates. To assess how different the true equilibrium values are, we employed an alternative method by measuring the transition of bR(f) to unfolded bacterioopsin (bO(u)), the RET-free form of unfolded protein. The bR(f)-to-bO(u) transition is fully reversible, particular when we add excess RET. We compared the difference in unfolding free energies for 13 bR mutants measured by both assays. For 12 of the 13 mutants with a wide range of stabilities, the results are essentially the same within experimental error. The congruence of the results is fortuitous and suggests the energetic effects of most mutations may be focused on the folded state. The bR(f)-to-bO(u) reaction is inconvenient because many days are required to reach equilibrium, but it is the preferable measure of thermodynamic stability. 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.

Simple screening method for improving membrane protein thermostability

Biochemical and biophysical analysis on integral membrane proteins often requires monodisperse and stable protein samples. Here we describe a method to characterize protein thermostability by measuring its melting temperature in detergent using analytical size-exclusion chromatography. This quantitative method can be used to screen for compounds and conditions that stabilize the protein. With this technique we were able to assess and improve the thermostability of several membrane proteins. These conditions were in turn used to assist purification, to identify protein ligand and to improve crystal quality.

Computational design of membrane proteins

This article reviews the recent successes of computational protein design techniques applied to integral membrane proteins. This emerging area is still handicapped by significant difficulties in the experimental characterization of the stability and structure of the designed proteins. Nevertheless, by focusing on oligomeric complexes of single-span transmembrane (TM) peptides with detectable activity, the computational design of membrane proteins has already produced very exciting results. The 'take-home message' is that optimization of van der Waals packing and hydrogen bonding (both 'canonical' and weak Cα-H⋯O bonds) can produce functional structures of remarkable stability and specificity in the membrane.

Contribution of hydrophobic interactions to protein stability

Our goal was to gain a better understanding of the contribution of hydrophobic interactions to protein stability. We measured the change in conformational stability, Δ(ΔG), for hydrophobic mutants of four proteins: villin headpiece subdomain (VHP) with 36 residues, a surface protein from Borrelia burgdorferi (VlsE) with 341 residues, and two proteins previously studied in our laboratory, ribonucleases Sa and T1. We compared our results with those of previous studies and reached the following conclusions: (1) Hydrophobic interactions contribute less to the stability of a small protein, VHP (0.6±0.3 kcal/mol per -CH(2)- group), than to the stability of a large protein, VlsE (1.6±0.3 kcal/mol per -CH(2)- group). (2) Hydrophobic interactions make the major contribution to the stability of VHP (40 kcal/mol) and the major contributors are (in kilocalories per mole) Phe18 (3.9), Met13 (3.1), Phe7 (2.9), Phe11 (2.7), and Leu21 (2.7). (3) Based on the Δ(ΔG) values for 148 hydrophobic mutants in 13 proteins, burying a -CH(2)- group on folding contributes, on average, 1.1±0.5 kcal/mol to protein stability. (4) The experimental Δ(ΔG) values for aliphatic side chains (Ala, Val, Ile, and Leu) are in good agreement with their ΔG(tr) values from water to cyclohexane. (5) For 22 proteins with 36 to 534 residues, hydrophobic interactions contribute 60±4% and hydrogen bonds contribute 40±4% to protein stability. (6) Conformational entropy contributes about 2.4 kcal/mol per residue to protein instability. The globular conformation of proteins is stabilized predominantly by hydrophobic interactions.

Mutations that probe the cooperative assembly of O⁶-alkylguanine-DNA alkyltransferase complexes

O(6)-Alkylguanine-DNA alkyltransferase (AGT) repairs mutagenic O(6)-alkylguanine and O(4)-alkylthymine adducts present in DNA that has been exposed to alkylating agents. AGT binds DNA cooperatively, and models of cooperative complexes predict that residues 1-7 of one protein molecule and residues 163-169 of a neighboring protein are closely juxtaposed. To test these models, we used directed mutagenesis to substitute triplets of alanine for triplets of native residues across these two sequences. Six of eight designed mutants expressed AGT at detectable levels. All mutant AGTs that were expressed were folded compactly, bound DNA with stoichiometries equivalent to that of the wild-type protein, and were able to protect Escherichia coli to varying degrees from the potent alkylating agent N-methyl-N'-nitro-N-nitrosoguanidine (MNNG). All mutations attenuated DNA binding cooperativity, but unexpectedly, they also reduced the affinity of AGT for DNA. This suggests that the protein-protein and protein-DNA interactions of AGT are strongly coupled. When normalized for differences in AGT expression, cells expressing mutants KDC(3-5)-AAA, DCE(4-6)-AAA, and KEW(165-167)-AAA were significantly more susceptible to MNNG than wild-type cells. This is the first evidence, to the best of our knowledge, of a role for residues at the protein-protein interface and, by implication, cooperative protein-protein interactions in the cell-protective mechanisms of AGT.

Viral channel forming proteins - modeling the target

The cellular and subcellular membranes encounter an important playground for the activity of membrane proteins encoded by viruses. Viral membrane proteins, similar to their host companions, can be integral or attached to the membrane. They are involved in directing the cellular and viral reproduction, the fusion and budding processes. This review focuses especially on those integral viral membrane proteins which form channels or pores, the classification to be so, modeling by in silico methods and potential drug candidates. The sequence of an isolate of Vpu from HIV-1 is aligned with host ion channels and a toxin. The focus is on the alignment of the transmembrane domains. The results of the alignment are mapped onto the 3D structures of the respective channels and toxin. The results of the mapping support the idea of a 'channel-pore dualism' for Vpu.

Residue-specific side-chain packing determines the backbone dynamics of transmembrane model helices

The transmembrane domains (TMDs) of membrane-fusogenic proteins contain an overabundance of β-branched residues. In a previous effort to systematically study the relation among valine content, fusogenicity, and helix dynamics, we developed model TMDs that we termed LV-peptides. The content and position of valine in LV-peptides determine their fusogenicity and backbone dynamics, as shown experimentally. Here, we analyze their conformational dynamics and the underlying molecular forces using molecular-dynamics simulations. Our study reveals that backbone dynamics is correlated with the efficiency of side-chain to side-chain van der Waals packing between consecutive turns of the helix. Leu side chains rapidly interconvert between two rotameric states, thus favoring contacts to its i±3 and i±4 neighbors. Stereochemical restraints acting on valine side chains in the α-helix force both β-substituents into an orientation where i,i±3 interactions are less favorable than i,i±4 interactions, thus inducing a local packing deficiency at VV3 motifs. We provide a quantitative molecular model to explain the relationship among chain connectivity, side-chain mobility, and backbone flexibility. We expect that this mechanism also defines the backbone flexibility of natural TMDs.

Probing membrane protein unfolding with pulse proteolysis

Technical challenges have greatly impeded the investigation of membrane protein folding and unfolding. To develop a new tool that facilitates the study of membrane proteins, we tested pulse proteolysis as a probe for membrane protein unfolding. Pulse proteolysis is a method to monitor protein folding and unfolding, which exploits the significant difference in proteolytic susceptibility between folded and unfolded proteins. This method requires only a small amount of protein and, in many cases, may be used with unpurified proteins in cell lysates. To evaluate the effectiveness of pulse proteolysis as a probe for membrane protein unfolding, we chose Halobacterium halobium bacteriorhodopsin (bR) as a model system. The denaturation of bR in SDS has been investigated extensively by monitoring the change in the absorbance at 560 nm (A(560)). In this work, we demonstrate that denaturation of bR by SDS results in a significant increase in its susceptibility to proteolysis by subtilisin. When pulse proteolysis was applied to bR incubated in varying concentrations of SDS, the remaining intact protein determined by electrophoresis shows a cooperative transition. The midpoint of the cooperative transition (C(m)) shows excellent agreement with that determined by A(560). The C(m) values determined by pulse proteolysis for M56A and Y57A bRs are also consistent with the measurements made by A(560). Our results suggest that pulse proteolysis is a quantitative tool to probe membrane protein unfolding. Combining pulse proteolysis with Western blotting may allow the investigation of membrane protein unfolding in situ without overexpression or purification.

Structural origins of nitroxide side chain dynamics on membrane protein α-helical sites

Understanding the structure and dynamics of membrane proteins in their native, hydrophobic environment is important to understanding how these proteins function. EPR spectroscopy in combination with site-directed spin labeling (SDSL) can measure dynamics and structure of membrane proteins in their native lipid environment; however, until now the dynamics measured have been qualitative due to limited knowledge of the nitroxide spin label's intramolecular motion in the hydrophobic environment. Although several studies have elucidated the structural origins of EPR line shapes of water-soluble proteins, EPR spectra of nitroxide spin-labeled proteins in detergents or lipids have characteristic differences from their water-soluble counterparts, suggesting significant differences in the underlying molecular motion of the spin label between the two environments. To elucidate these differences, membrane-exposed α-helical sites of the leucine transporter, LeuT, from Aquifex aeolicus, were investigated using X-ray crystallography, mutational analysis, nitroxide side chain derivatives, and spectral simulations in order to obtain a motional model of the nitroxide. For each crystal structure, the nitroxide ring of a disulfide-linked spin label side chain (R1) is resolved and makes contacts with hydrophobic residues on the protein surface. The spin label at site I204 on LeuT makes a nontraditional hydrogen bond with the ortho-hydrogen on its nearest neighbor F208, whereas the spin label at site F177 makes multiple van der Waals contacts with a hydrophobic pocket formed with an adjacent helix. These results coupled with the spectral effect of mutating the i ± 3, 4 residues suggest that the spin label has a greater affinity for its local protein environment in the low dielectric than on a water-soluble protein surface. The simulations of the EPR spectra presented here suggest the spin label oscillates about the terminal bond nearest the ring while maintaining weak contact with the protein surface. Combined, the results provide a starting point for determining a motional model for R1 on membrane proteins, allowing quantification of nitroxide dynamics in the aliphatic environment of detergent and lipids. In addition, initial contributions to a rotamer library of R1 on membrane proteins are provided, which will assist in reliably modeling the R1 conformational space for pulsed dipolar EPR and NMR paramagnetic relaxation enhancement distance determination.

Unfolding free energy of a two-domain transmembrane sugar transport protein

Understanding how an amino acid sequence folds into a functional, three-dimensional structure has proved to be a formidable challenge in biological research, especially for transmembrane proteins with multiple alpha helical domains. Mechanistic folding studies on helical membrane proteins have been limited to unusually stable, single domain proteins such as bacteriorhodopsin. Here, we extend such work to flexible, multidomain proteins and one of the most widespread membrane transporter families, the major facilitator superfamily, thus showing that more complex membrane proteins can be successfully refolded to recover native substrate binding. We determine the unfolding free energy of the two-domain, Escherichia coli galactose transporter, GalP; a bacterial homologue of human glucose transporters. GalP is reversibly unfolded by urea. Urea causes loss of substrate binding and a significant reduction in alpha helical content. Full recovery of helical structure and substrate binding occurs in dodecylmaltoside micelles, and the unfolding free energy can be determined. A linear dependence of this free energy on urea concentration allows the free energy of unfolding in the absence of urea to be determined as +2.5 kcal·mol(-1). Urea has often been found to be a poor denaturant for transmembrane helical structures. We attribute the denaturation of GalP helices by urea to the dynamic nature of the transporter structure allowing denaturant access via the substrate binding pocket, as well as to helical structure that extends beyond the membrane. This study gives insight into the final, critical folding step involving recovery of ligand binding for a multidomain membrane transporter.

The contribution of a covalently bound cofactor to the folding and thermodynamic stability of an integral membrane protein

The factors controlling the stability, folding, and dynamics of integral membrane proteins are not fully understood. The high stability of the membrane protein bacteriorhodopsin (bR), an archetypal member of the rhodopsin photoreceptor family, has been ascribed to its covalently bound retinal cofactor. We investigate here the role of this cofactor in the thermodynamic stability and folding kinetics of bR. Multiple spectroscopic probes were used to determine the kinetics and energetics of protein folding in mixed lipid/detergent micelles in the presence and absence of retinal. The presence of retinal increases extrapolated values for the overall unfolding free energy from 6.3 ± 0.4 kcal mol(-1) to 23.4 ± 1.5 kcal mol(-1) at zero denaturant, suggesting that the cofactor contributes 17.1 kcal mol(-1) towards the overall stability of bR. In addition, the cooperativity of equilibrium unfolding curves is markedly reduced in the absence of retinal with overall m-values decreasing from 31.0 ± 2.0 kcal mol(-1) to 10.9 ± 1.0 kcal mol(-1), indicating that the folded state of the apoprotein is less compact than the equivalent for the holoprotein. This change in the denaturant response means that the difference in the unfolding free energy at a denaturant concentration midway between the two unfolding curves is only ca 3-6 kcal mol(-1). Kinetic data show that the decrease in stability upon removal of retinal is associated with an increase in the apparent intrinsic rate constant of unfolding, k(u)(H2O), from ~1 × 10(-16) s(-1) to ~1 × 10(-4) s(-1) at 25 °C. This correlates with a decrease in the unfolding activation energy by 16.3 kcal mol(-1) in the apoprotein, extrapolated to zero SDS. These results suggest that changes in bR stability induced by retinal binding are mediated solely by changes in the activation barrier for unfolding. The results are consistent with a model in which bR is kinetically stabilized via a very slow rate of unfolding arising from protein-retinal interactions that increase the rigidity and compactness of the polypeptide chain.

ORF8a of SARS-CoV forms an ion channel: experiments and molecular dynamics simulations

ORF8a protein is 39 residues long and contains a single transmembrane domain. The protein is synthesized using solid phase peptide synthesis and reconstituted into artificial lipid bilayers that forms cation-selective ion channels with a main conductance level of 8.9±0.8pS at elevated temperature (38.5°C). Computational modeling studies including multi nanosecond molecular dynamics simulations in a hydrated POPC lipid bilayer are done with a 22 amino acid transmembrane helix to predict a putative homooligomeric helical bundle model. A structural model of a pentameric bundle is proposed with cysteines, serines and threonines facing the pore.

Protein folding in membranes

Separation of cells and organelles by bilayer membranes is a fundamental principle of life. Cellular membranes contain a baffling variety of proteins, which fulfil vital functions as receptors and signal transducers, channels and transporters, motors and anchors. The vast majority of membrane-bound proteins contain bundles of alpha-helical transmembrane domains. Understanding how these proteins adopt their native, biologically active structures in the complex milieu of a membrane is therefore a major challenge in today's life sciences. Here, we review recent progress in the folding, unfolding and refolding of alpha-helical membrane proteins and compare the molecular interactions that stabilise proteins in lipid bilayers. We also provide a critical discussion of a detergent denaturation assay that is increasingly used to determine membrane-protein stability but is not devoid of conceptual difficulties.

Practical considerations of membrane protein instability during purification and crystallisation

Crystallisation of integral membranes requires milligrams of purified protein in a homogeneous, monodisperse state, and crucially, the membrane protein must also be fully functional and stable. The stability of membrane proteins in solution is dependent on the type of detergents used, but unfortunately the use of the most stabilising detergent can often decrease the probability of obtaining crystals that diffract to high resolution, especially of small membrane proteins. A number of strategies have been developed to facilitate the purification of membrane proteins in a functional form, which have led to new possibilities for structure determination.

Similar energetic contributions of packing in the core of membrane and water-soluble proteins

A major driving force for water-soluble protein folding is the hydrophobic effect, but membrane proteins cannot make use of this stabilizing contribution in the apolar core of the bilayer. It has been proposed that membrane proteins compensate by packing more efficiently. We therefore investigated packing contributions experimentally by observing the energetic and structural consequences of cavity creating mutations in the core of a membrane protein. We observed little difference in the packing energetics of water and membrane soluble proteins. Our results imply that other mechanisms are employed to stabilize the structure of membrane proteins.

Sequential dissociation of subunits from bovine heart cytochrome C oxidase by urea

The quaternary stability of purified, detergent-solubilized, cytochrome c oxidase (CcO) was probed using two chemical denaturants, urea and guanidinium chloride (GdmCl). Each chaotrope induces dissociation of five subunits in a concentration-dependent manner. These five subunits are not scattered over the surface of CcO but are clustered together in close contact at the dimer interface. Increasing the concentration of urea selectively dissociates subunits from CcO in the following order: VIa and VIb, followed by III and VIIa, and finally Vb. After incubation in urea for 10 min at room temperature, the sigmoidal dissociation transitions were centered at 3.7, 4.6, and 7.0 M urea, respectively. The secondary structure of CcO was only minimally perturbed, indicating that urea causes disruption of subunit interactions without urea-induced conformational changes. Incubation of CcO in urea for 120 min produced similar results but shifted the sigmoidal dissociation curves to lower urea concentrations. Incubation of CcO with increasing concentrations of GdmCl produces an analogous effect; however, the GdmCl-induced dissociation of subunits occurs at lower concentrations and with a narrower concentration range. Thermodynamic parameters for each subunit dissociation were evaluated from the sigmoidal dissociation data by assuming a single transition from bound to dissociated subunit. The free energy change accompanying urea-induced dissociation of each subunit ranged from 18.0 to 29.7 kJ/mol, which corresponds to 0.32-0.59 kJ/mol per 100 A(2) of newly exposed solvent-accessible surface area. These values are 30-50-fold smaller than previously reported for the unfolding of soluble or membrane proteins.

Calcium binding to the purple membrane: A molecular dynamics study

The purple membrane (PM) is a specialized membrane patch found in halophilic archaea, containing the photoreceptor bacteriorhodopsin (bR). It is long known that calcium ions bind to the PM, but their position and role remain elusive to date. Molecular dynamics simulations in conjunction with a highly detailed model of the PM have been used to investigate the stability of calcium ions placed at three proposed cation binding sites within bR, one near the Schiff base, one in the region of the proton release group, and one near Glu9. The simulations suggest that, of the sites investigated, the binding of calcium ions was most likely at the proton release group. Binding in the region of the Schiff base, while possible, was associated with significant changes in local geometry. Calcium ions placed near Glu9 in the interior of bR (simultaneously to a Ca(2+) near the Schiff base and another one near the Glu194-Glu204 site) were not stable. The results obtained are discussed in relation to recent experimental observations and theoretical considerations.

Folding scene investigation: membrane proteins

Investigations into protein folding have concentrated on experimentally tractable proteins with the result that membrane protein folding remains unsolved. New evidence is providing insight into the nature of the interactions stabilising the folded state of α-helical membrane proteins as well as giving hints on the character of the folding transition state. These developments show that classical methods used for water-soluble proteins can be successfully adapted for membrane proteins. The advances, coupled with increasing numbers of solved crystal structures, augur well for future research into the mechanisms of membrane protein folding.

The transition state for integral membrane protein folding

Biology relies on the precise self-assembly of its molecular components. Generic principles of protein folding have emerged from extensive studies on small, water-soluble proteins, but it is unclear how these ideas are translated into more complex situations. In particular, the one-third of cellular proteins that reside in biological membranes will not fold like water-soluble proteins because membrane proteins need to expose, not hide, their hydrophobic surfaces. Here, we apply the powerful protein engineering method of Phi-value analysis to investigate the folding transition state of the alpha-helical membrane protein, bacteriorhodopsin, from a partially unfolded state. Our results imply that much of helix B of the seven-transmembrane helical protein is structured in the transition state with single-point alanine mutations in helix B giving Phi values >0.8. However, residues Y43 and T46 give lower Phi values of 0.3 and 0.5, respectively, suggesting a possible reduction in native structure in this region of the helix. Destabilizing mutations also increase the activation energy of folding, which is accompanied by an apparent movement of the transition state toward the partially unfolded state. This apparent transition state movement is most likely due to destabilization of the structured, unfolded state. These results contrast with the Hammond effect seen for several water-soluble proteins in which destabilizing mutations cause the transition state to move toward, and become closer in energy to, the folded state. We thus introduce a classic folding analysis method to membrane proteins, providing critical insight into the folding transition state.