Ultrafast Protein Folding in Membrane-Mimetic Environments

Proteins fold on timescales from hours to microseconds. In addition to protein size, sequence, and topology, the environment represents an equally important factor in determining folding speed. This is particularly relevant for proteins that require a lipid membrane or a membrane mimic to fold. However, only little is known about how properties of such a hydrophilic/hydrophobic interface modulate the folding landscape of membrane-interacting proteins. Here, we studied the influence of different membrane-mimetic micellar environments on the folding and unfolding kinetics of the helical-bundle protein Mistic. Devising a single-molecule fluorescence spectroscopy approach, we extracted folding and unfolding rates under equilibrium conditions and dissected the contributions from different detergent moieties to the free-energy landscape. While both polar and nonpolar moieties contribute to stability, they exert differential effects on the free-energy barrier: Hydrophobic burial stabilizes the folded state but not the transition state in reference to a purely aqueous environment; by contrast, zwitterionic headgroup moieties stabilize the folded state and, additionally, lower the free-energy barrier to accelerate the folding of Mistic to achieve ultrafast folding times down to 35 μs.     

Thermodynamics of unfolding of an integral membrane protein in mixed micelles

Quantitative studies of membrane protein folding and unfolding can be difficult because of difficulties with efficient refolding as well as a pronounced propensity to aggregate. However, mixed micelles, consisting of the anionic detergent sodium dodecyl sulfate and the nonionic detergent dodecyl maltoside facilitate reversible and quantitative unfolding and refolding. The 4-transmembrane helix protein DsbB from the inner membrane of Escherichia coli unfolds in mixed micelles according to a three-state mechanism involving an unfolding intermediate I. The temperature dependence of the kinetics of this reaction between 15 degrees and 45 degrees C supports that unfolding from I to the denatured state D is accompanied by a significant decrease in heat capacity. For water-soluble proteins, the heat capacity increases upon unfolding, and this is generally interpreted as the increased binding of water to the protein as it unfolds, exposing more surface area. The decrease in DsbB's heat capacity upon unfolding is confirmed by independent thermal scans. The decrease in heat capacity is not an artifact of the use of mixed micelles, since the water soluble protein S6 shows conventional heat-capacity changes in detergent. We speculate that it reflects the binding of SDS to parts of DsbB that are solvent-exposed in the native DM-bound state. This implies that the periplasmic loops of DsbB are relatively unstructured. This anomalous thermodynamic behavior has not been observed for beta-barrel membrane proteins, probably because they do not bind SDS so extensively. Thus the thermodynamic behavior of membrane proteins appears to be intimately connected to their detergent-binding properties.     

Transmembrane helix-helix interactions: comparative simulations of the glycophorin a dimer

The glycophorin helix dimer is a paradigm for the exploration of helix-helix interactions in integral membrane proteins. Two NMR structures of the dimer are known, one in a detergent micelle and one in a lipid bilayer. Multiple (4 x 50 ns) molecular dynamics simulations starting from each of the two NMR structures, with each structure in either a dodecyl phosphocholine (DPC) micelle or a dimyristoyl phosphatidylcholine (DMPC) bilayer, have been used to explore the conformational dynamics of the helix dimer. Analysis of the helix-helix interaction, mediated by the GxxxG sequence motif, suggests convergence of the simulations to a common model. This is closer to the NMR structure determined in a bilayer than to micelle structure. The stable dimer interface in the final simulation model is characterized by (i) Gly/Gly packing and (ii) Thr/Thr interhelix H-bonds. These results demonstrate the ability of extended molecular dynamics simulations in a lipid bilayer environment to refine membrane protein structures or models derived from experimental data obtained in protein/detergent micelles.     

Expression, purification, and characterization of Thermotoga maritima membrane proteins for structure determination

Structural studies of integral membrane proteins typically rely upon detergent micelles as faithful mimics of the native lipid bilayer. Therefore, membrane protein structure determination would be greatly facilitated by biophysical techniques that are capable of evaluating and assessing the fold and oligomeric state of these proteins solubilized in detergent micelles. In this study, an approach to the characterization of detergent-solubilized integral membrane proteins is presented. Eight Thermotoga maritima membrane proteins were screened for solubility in 11 detergents, and the resulting soluble protein-detergent complexes were characterized with small angle X-ray scattering (SAXS), nuclear magnetic resonance (NMR) spectroscopy, circular dichroism (CD) spectroscopy, and chemical cross-linking to evaluate the homogeneity, oligomeric state, radius of gyration, and overall fold. A new application of SAXS is presented, which does not require density matching, and NMR methods, typically used to evaluate soluble proteins, are successfully applied to detergent-solubilized membrane proteins. Although detergents with longer alkyl chains solubilized the most proteins, further characterization indicates that some of these protein-detergent complexes are not well suited for NMR structure determination due to conformational exchange and protein oligomerization. These results emphasize the need to screen several different detergents and to characterize the protein-detergent complex in order to pursue structural studies. Finally, the physical characterization of the protein-detergent complexes indicates optimal solution conditions for further structural studies for three of the eight overexpressed membrane proteins.     

Unfolding a transmembrane helix dimer: A FRET study in mixed micelles

The exact nature of membrane protein folding and assembly is not understood in detail yet. Addition of SDS to a membrane protein dissolved in mild, non-polar detergent results in formation of mixed micelles and in subsequent denaturation of higher ordered membrane protein structures. The exact nature of this denaturation event is, however, enigmatic, and separation of an individual helix pair in mixed micelles has also not been reported yet. Here we followed unfolding of the human glycophorin A transmembrane helix dimer in mixed micelles by fluorescence spectroscopy. Energy transfer between differently labelled glycophorin A transmembrane helices decreased with increasing SDS mole fractions albeit without modifying the helicity of the peptides. The energetics and kinetics of the dimer dissociation in mixed micelles is analyzed and discussed, and the experimental data demonstrate that mixed micelles can be used as a general method to investigate unfolding of α-helical membrane proteins.     

Protein structures in SDS micelle-protein complexes.

Sodium dodecyl sulfate (SDS) is used more often than any other detergent as an excellent denaturing or "unfolding" detergent. However, formation of ordered structure (alpha-helix or beta-sheet) in certain peptides is known to be induced by interaction with SDS micelles. The SDS-induced structures formed by these peptides are amphiphilic, having both a hydrophobic and a hydrophilic face. Previous work in this area has revealed that SDS induces helical folding in a wide variety of non-helical proteins. Here, we describe the interaction of several structurally unrelated proteins with SDS micelles and the correlation of these structures to helical amphiphilic regions present in the primary sequence. It is likely that the ability of native nonordered protein structures to form induced amphiphilic ordered structures is rather common.     

Molecular dynamics simulations of a bacterial autotransporter: NalP from Neisseria meningitidis

NalP is an autotransporter secretory protein found in the outer membrane of Neisseria meningitidis. The crystal structure of the NalP translocator domain revealed a transmembrane beta-barrel containing a central alpha-helix. The role of this alpha-helix, and of the conformational dynamics of the beta-barrel pore have been studied via atomistic molecular dynamics simulations. Three simulations, each of 10 ns duration, of NalP embedded within a solvated DMPC bilayer were performed. The helix was removed from the barrel interior in one simulation. The conformational stability of the protein is similar to that of other outer membrane proteins, e.g., OmpA, in comparable simulations. The transmembrane beta-barrel is stable even in the absence of the alpha-helix. Removal of the helix results in an influx of water into the pore region, suggesting the helix acts as a 'plug'. Water molecules entering the resultant pore form hydrogen bonds with the barrel lining that compensate for the loss of helix-barrel hydrogen bonds. The dimensions of the pore fluctuate over the course of the simulation revealing it to be flexible, but only wide enough to allow transport of the passenger domain in an unfolded or extended conformation. The simulations help us to understand the role of the central helix in plugging the pore and in maintaining the width of the barrel, and show that the NalP monomer is sufficient for the transport of the passenger domain in an unfolded or extended conformation.     

Rapid exploration of the folding topology of helical membrane proteins using paramagnetic perturbation

An understanding of the folding states of α-helical membrane proteins in detergent systems is important for functional and structural studies of these proteins. Here, we present a rapid and simple method for identification of the folding topology and assembly of transmembrane helices using paramagnetic perturbation in nuclear magnetic resonance spectroscopy. By monitoring the perturbation of signals from glycine residues located at specific sites, the folding topology and the assembly of transmembrane helices of membrane proteins were easily identified without time-consuming backbone assignment. This method is validated with Mistic (membrane-integrating sequence for translation of integral membrane protein constructs) of known structure as a reference protein. The folding topologies of two bacterial histidine kinase membrane proteins (SCO3062 and YbdK) were investigated by this method in dodecyl phosphocholine (DPC) micelles. Combing with analytical ultracentrifugation, we identified that the transmembrane domain of YbdK is present as a parallel dimer in DPC micelle. In contrast, the interaction of transmembrane domain of SCO3062 is not maintained in DPC micelle due to disruption of native structure of the periplasmic domain by DPC micelle.     

Thermodynamic and kinetic characterization of a beta-hairpin peptide in solution: an extended phase space sampling by molecular dynamics simulations in explicit water

The folding of the amyloidogenic H1 peptide MKHMAGAAAAGAVV taken from the syrian hamster prion protein is explored in explicit aqueous solution at 300 K using long time scale all-atom molecular dynamics simulations for a total simulation time of 1.1 mus. The system, initially modeled as an alpha-helix, preferentially adopts a beta-hairpin structure and several unfolding/refolding events are observed, yielding a very short average beta-hairpin folding time of approximately 200 ns. The long time scale accessed by our simulations and the reversibility of the folding allow to properly explore the configurational space of the peptide in solution. The free energy profile, as a function of the principal components (essential eigenvectors) of motion, describing the main conformational transitions, shows the characteristic features of a funneled landscape, with a downhill surface toward the beta-hairpin folded basin. However, the analysis of the peptide thermodynamic stability, reveals that the beta-hairpin in solution is rather unstable. These results are in good agreement with several experimental evidences, according to which the isolated H1 peptide adopts very rapidly in water beta-sheet structure, leading to amyloid fibril precipitates [Nguyen et al., Biochemistry 1995;34:4186-4192; Inouye et al., J Struct Biol 1998;122:247-255]. Moreover, in this article we also characterize the diffusion behavior in conformational space, investigating its relations with folding/unfolding conditions.     

Bacterial expression of a eukaryotic membrane protein in fusion to various Mistic orthologs

Mistic, a bacterial membrane-associating protein family, uniquely found in Bacillus species. It enhances expression of eukaryotic membrane proteins at the bacterial membrane. Mistic from B. subtilis (M110), expresses at the Escherichia coli membrane, however its shorter orthologs have been recently shown to be mainly cytoplasmic with varying membrane affinities. Based on that, we hypothesized that the expression level of membrane proteins fused to Mistic is correlated with the degree of membrane association of the particular Mistic protein. We compared expression levels by various Mistic proteins as fusion partners for the Aplysia californica Kv1.1 (aKv1.1) channel as a cargo membrane protein. Mistic from B. atrophaeus (M4), which has the highest membrane association among the shorter orthologs, enhanced expression of the transmembrane domain of aKv1.1 to the highest extent. In contrast, M1, which consists of the 84 C-terminal amino acids of M110 is the most soluble protein and showed the least capacity to express the channel. A chimeric Mistic, constructed with the first α-helix (H1) of M110 N-terminally fused to M4, did not increase the level of expression of aKv1.1 beyond those of either the M110 or the M4 fusions. The channel fused to M110, M4 or the aforementioned H1-M4 chimera, expresses in the highest quantity and quality among Mistic proteins, providing suitable sample for structural studies. Our data support the concept that expression levels of ‘Misticated’ membrane proteins are related to the independent chaperoning character of Mistic via direct membrane association, rather than related to specific sequence-dependent interaction with the E. coli translocon machinery.     

Unfolding of beta-sheet proteins in SDS

Beta-sheet proteins are particularly resistant to denaturation by sodium dodecyl sulfate (SDS). Here we compare unfolding of two beta-sandwich proteins TNfn3 and TII27 in SDS. The two proteins show different surface electrostatic potential. Correspondingly, TII27 unfolds below the critical micelle concentration via the formation of hemimicelles on the protein surface, whereas TNfn3 only unfolds around the critical micelle concentration. Isothermal titration calorimetry confirms that unfolding of TII27 sets in at lower SDS concentrations, although the total number of bound SDS molecules is similar at the end of unfolding. In mixed micelles with the nonionic detergent dodecyl maltoside, where the concentration of monomeric SDS is insignificant, the behavior of the two proteins converges. TII27 unfolds more slowly than TNfn3 in SDS and follows a two-mode behavior. Additionally TNfn3 shows inhibition of SDS unfolding at intermediate SDS concentrations. Mutagenic analysis suggests that the overall unfolding mechanism is similar to that observed in denaturant for both proteins. Our data confirm the kinetic robustness of beta-sheet proteins toward SDS. We suggest this is related to the inability of SDS to induce significant amounts of alpha-helix structure in these proteins as part of the denaturation process, forcing the protein to denature by global rather than local unfolding.     

Protein unfolding in detergents: effect of micelle structure,ionic strength,pH, and temperature

The 101-residue monomeric protein S6 unfolds in the anionic detergent sodium dodecyl sulfate (SDS) above the critical micelle concentration, with unfolding rates varying according to two different modes. Our group has proposed that spherical micelles lead to saturation kinetics in unfolding (mode 1), while cylindrical micelles prevalent at higher SDS concentrations induce a power-law dependent increase in the unfolding rate (mode 2). Here I investigate in more detail how micellar properties affect protein unfolding. High NaCl concentrations, which induce cylindrical micelles, favor mode 2. This is consistent with our model, though other effects such as electrostatic screening cannot be discounted. Furthermore, unfolding does not occur in mode 2 in the cationic detergent LTAB, which is unable to form cylindrical micelles. A strong retardation of unfolding occurs at higher LTAB concentrations, possibly due to the formation of dead-end protein-detergent complexes. A similar, albeit much weaker, effect is seen in SDS in the absence of salt. Chymotrypsin inhibitor 2 exhibits the same modes of unfolding in SDS as S6, indicating that this type of protein unfolding is not specific for S6. The unfolding process in mode 1 has an activation barrier similar in magnitude to that in water, while the activation barrier in mode 2 is strongly concentration-dependent. The strong pH-dependence of unfolding in SDS and LTAB suggests that the rate of unfolding in anionic detergent is modulated by repulsion between detergent headgroups and anionic side chains, while cationic side chains modulate unfolding rates in cationic detergents.     

Membrane proteins, lipids and detergents: not just a soap opera

Studying membrane proteins represents a major challenge in protein biochemistry, with one of the major difficulties being the problems encountered when working outside the natural lipid environment. In vitro studies such as crystallization are reliant on the successful solubilization or reconstitution of membrane proteins, which generally involves the careful selection of solubilizing detergents and mixed lipid/detergent systems. This review will concentrate on the methods currently available for efficient reconstitution and solubilization of membrane proteins through the use of detergent micelles, mixed lipid/detergent micelles and bicelles or liposomes. We focus on the relevant molecular properties of the detergents and lipids that aid understanding of these processes. A significant barrier to membrane protein research is retaining the stability and function of the protein during solubilization, reconstitution and crystallization. We highlight some of the lessons learnt from studies of membrane protein folding in vitro and give an overview of the role that lipids can play in stabilizing the proteins.     

Correct folding of the beta-barrel of the human membrane protein VDAC requires a lipid bilayer

Spontaneous membrane insertion and folding of beta-barrel membrane proteins from an unfolded state into lipid bilayers has been shown previously only for few outer membrane proteins of Gram-negative bacteria. Here we investigated membrane insertion and folding of a human membrane protein, the isoform 1 of the voltage-dependent anion-selective channel (hVDAC1) of mitochondrial outer membranes. Two classes of transmembrane proteins with either alpha-helical or beta-barrel membrane domains are known from the solved high-resolution structures. VDAC forms a transmembrane beta-barrel with an additional N-terminal alpha-helix. We demonstrate that similar to bacterial OmpA, urea-unfolded hVDAC1 spontaneously inserts and folds into lipid bilayers upon denaturant dilution in the absence of folding assistants or energy sources like ATP. Recordings of the voltage-dependence of the single channel conductance confirmed folding of hVDAC1 to its active form. hVDAC1 developed first beta-sheet secondary structure in aqueous solution, while the alpha-helical structure was formed in the presence of lipid or detergent. In stark contrast to bacterial beta-barrel membrane proteins, hVDAC1 formed different structures in detergent micelles and phospholipid bilayers, with higher content of beta-sheet and lower content of alpha-helix when inserted and folded into lipid bilayers. Experiments with mixtures of lipid and detergent indicated that the content of beta-sheet secondary structure in hVDAC1 decreased at increased detergent content. Unlike bacterial beta-barrel membrane proteins, hVDAC1 was not stable even in mild detergents such as LDAO or dodecylmaltoside. Spontaneous folding of outer membrane proteins into lipid bilayers indicates that in cells, the main purpose of membrane-inserted or associated assembly factors may be to select and target beta-barrel membrane proteins towards the outer membrane instead of actively assembling them under consumption of energy as described for the translocons of cytoplasmic membranes.     

Expression of a chloroplast ATP/ADP transporter in E. coli membranes: behind the Mistic strategy

Eukaryotic membrane protein expression is still a major bottleneck for structural studies. Production in E. coli often leads to low expression level and/or aggregated proteins. In the last decade, strategies relying on new fusion protein expression revealed promising results. Fusion with the amphipatic Mistic protein has been described to favor expression in E. coli membranes. Although, this approach has already been reported for a few membrane proteins, little is known about the activity of the fused proteins. We used this strategy and obtained high expression levels of a chloroplast ATP/ADP transporter from A. thaliana (NTT1) and characterized its transport properties. NTT1 fused to Mistic has a very low transport activity which can be recovered after in vivo Mistic fusion cleavage. Moreover, detailed molecular characterization of purified NTT1 mature form, NTT1 fused to Mistic or NTT1 cleaved-off from this fusion highlights the correct fold of the latter one. Therefore, considering the higher quantity of purified NTT1 mature form obtained via the Mistic fusion approach, this is a valuable strategy for obtaining quantities of pure and active proteins that are adequate for structural studies.     

How Amphipols Embed Membrane Proteins: Global Solvent Accessibility and Interaction with a Flexible Protein Terminus

Amphipathic polymers called amphipols provide a valuable alternative to detergents for keeping integral membrane proteins soluble in aqueous buffers. Here, we characterize spatial contacts of amphipol A8-35 with membrane proteins from two architectural classes: The 8-stranded β-barrel outer membrane protein OmpX and the α-helical protein bacteriorhodopsin. OmpX is well structured in A8-35, with its barrel adopting a fold closely similar to that in dihexanoylphosphocholine micelles. The accessibility of A8-35-trapped OmpX by a water-soluble paramagnetic molecule is highly similar to that in detergent micelles and resembles the accessibility in the natural membrane. For the α-helical protein bacteriorhodopsin, previously shown to keep its fold and function in amphipols, NMR data show that the imidazole protons of a polyhistidine tag at the N-terminus of the protein are exchange protected in the presence of detergent and lipid bilayer nanodiscs, but not in amphipols, indicating the absence of an interaction in the latter case. Overall, A8-35 exhibits protein interaction properties somewhat different from detergents and lipid bilayer nanodiscs, while maintaining the structure of solubilized integral membrane proteins.     

A Comprehensive Multidimensional-Embedded, One-Dimensional Reaction Coordinate for Protein Unfolding/Folding

The goal of the Dynameomics project is to perform, store, and analyze molecular dynamics simulations of representative proteins, of all known globular folds, in their native state and along their unfolding pathways. To analyze unfolding simulations, the location of the protein along the unfolding reaction coordinate (RXN) must be determined. Properties such as the fraction of native contacts and radius of gyration are often used; however, there is an issue regarding degeneracy with these properties, as native and nonnative species can overlap. Here, we used 15 physical properties of the protein to construct a multidimensional-embedded, one-dimensional RXN coordinate that faithfully captures the complex nature of unfolding. The unfolding RXN coordinates for 188 proteins (1534 simulations and 22.9 μs in explicit water) were calculated. Native, transition, intermediate, and denatured states were readily identified with the use of this RXN coordinate. A global native ensemble based on the native-state properties of the 188 proteins was created. This ensemble was shown to be effective for calculating RXN coordinates for folds outside the initial 188 targets. These RXN coordinates enable, high-throughput assignment of conformational states, which represents an important step in comparing protein properties across fold space as well as characterizing the unfolding of individual proteins.     

Conformation and environment of channel-forming peptides: a simulation study

Ion channel-forming peptides enable us to study the conformational dynamics of a transmembrane helix as a function of sequence and environment. Molecular dynamics simulations are used to study the conformation and dynamics of three 22-residue peptides derived from the second transmembrane domain of the glycine receptor (NK4-M2GlyR-p22). Simulations are performed on the peptide in four different environments: trifluoroethanol/water; SDS micelles; DPC micelles; and a DMPC bilayer. A hierarchy of alpha-helix stabilization between the different environments is observed such that TFE/water < micelles < bilayers. Local clustering of trifluoroethanol molecules around the peptide appears to help stabilize an alpha-helical conformation. Single (S22W) and double (S22W,T19R) substitutions at the C-terminus of NK4-M2GlyR-p22 help to stabilize a helical conformation in the micelle and bilayer environments. This correlates with the ability of the W22 and R19 side chains to form H-bonds with the headgroups of lipid or detergent molecules. This study provides a first atomic resolution comparison of the structure and dynamics of NK4-M2GlyR-p22 peptides in membrane and membrane-mimetic environments, paralleling NMR and functional studies of these peptides.     

A multidomain outer membrane protein from Pasteurella multocida: modelling and simulation studies of PmOmpA

PmOmpA is a two-domain outer membrane protein from Pasteurella multocida. The N-terminal domain of PmOmpA is a homologue of the transmembrane beta-barrel domain of OmpA from Escherichia coli, whilst the C-terminal domain of PmOmpA is a homologue of the extra-membrane Neisseria meningitidis RmpM C-terminal domain. This enables a model of a complete two domain PmOmpA to be constructed and its conformational dynamics explored via MD simulations of the protein embedded within two different phospholipid bilayers (DMPC and DMPE). The conformational stability of the transmembrane beta-barrel is similar to that of a homology model of OprF from Pseudomonas aeruginosa in bilayer simulations. There is a degree of water penetration into the interior of the beta-barrel, suggestive of a possible transmembrane pore. Although the PmOmpA model is stable over 20 ns simulations, retaining its secondary structure and fold integrity throughout, substantial flexibility is observed in a short linker region between the N- and the C-terminal domains. At low ionic strength, the C-terminal domain moves to interact electrostatically with the lipid bilayer headgroups. This study demonstrates that computational approaches may be applied to more complex, multi-domain outer membrane proteins, rather than just to transmembrane beta-barrels, opening the possibility of in silico proteomics approaches to such proteins.     

Solution NMR studies of the integral membrane proteins OmpX and OmpA from Escherichia coli

Membrane proteins are usually solubilized in polar solvents by incorporation into micelles. Even for small membrane proteins these mixed micelles have rather large molecular masses, typically beyond 50000 Da. The NMR technique TROSY (transverse relaxation-optimized spectroscopy) has been developed for studies of structures of this size in solution. In this paper, strategies for the use of TROSY-based NMR experiments with membrane proteins are discussed and illustrated with results obtained with the Escherichia coli integral membrane proteins OmpX and OmpA in mixed micelles with the detergent dihexanoylphosphatidylcholine (DHPC). For OmpX, complete sequence-specific NMR assignments have been obtained for the polypeptide backbone. The 13C chemical shifts and nuclear Overhauser effect data then resulted in the identification of the regular secondary structure elements of OmpX/DHPC in solution, and in the collection of an input of conformational constraints for the computation of the global fold of the protein. For OmpA, the NMR assignments are so far limited to about 80% of the polypeptide chain, indicating different dynamic properties of the reconstituted OmpA beta-barrel from those of OmpX. Overall, the present data demonstrate that relaxation-optimized NMR techniques open novel avenues for studies of structure, function and dynamics of integral membrane proteins.