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.     

Characterizing molecular interactions in different bacteriorhodopsin assemblies by single-molecule force spectroscopy

Using single-molecule force spectroscopy we characterized inter- and intramolecular interactions stabilizing structural segments of individual bacteriorhodopsin (BR) molecules assembled into trimers and dimers, and monomers. While the assembly of BR did not vary the location of these structural segments, their intrinsic stability could change up to 70% increasing from monomer to dimer to trimer. Since each stable structural segment established one unfolding barrier, we conclude that the locations of unfolding barriers were determined by intramolecular interactions but that their strengths were strongly influenced by intermolecular interactions. Subtracting the unfolding forces of the BR trimer from that of monomer allowed us to calculate the contribution of inter- and intramolecular interactions to the membrane protein stabilization. Statistical analyses showed that the unfolding pathways of differently assembled BR molecules did not differ in their appearance but in their population. This suggests that in our experiments the membrane protein assembly does not necessarily change the location of unfolding barriers within the protein, but certainly their strengths, and thus alters the probability of a protein to choose certain unfolding pathways.     

Probing force-induced unfolding intermediates of a single staphylococcal nuclease molecule and the effect of ligand binding

Single-molecule manipulation techniques have given experimental access to unfolding intermediates of proteins that are inaccessible in conventional experiments. A detailed characterization of the intermediates is a challenging problem that provides new possibilities for directly probing the energy landscape of proteins. We investigated single-molecule mechanical unfolding of a small globular protein, staphylococcal nuclease (SNase), using atomic force microscopy. The unfolding trajectories of the protein displayed sub-molecular and stochastic behavior with typical lengths corresponding to the size of the unfolded substructures. Our results support the view that the single protein unfolds along multiple pathways as suggested in recent theoretical studies. Moreover, we found the drastic change, caused by the ligand and inhibitor bindings, in the mechanical unfolding dynamics.     

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.     

Probing the energy landscape of the membrane protein bacteriorhodopsin

The folding and stability of transmembrane proteins is a fundamental and unsolved biological problem. Here, single bacteriorhodopsin molecules were mechanically unfolded from native purple membranes using atomic force microscopy and force spectroscopy. The energy landscape of individual transmembrane alpha helices and polypeptide loops was mapped by monitoring the pulling speed dependence of the unfolding forces and applying Monte Carlo simulations. Single helices formed independently stable units stabilized by a single potential barrier. Mechanical unfolding of the helices was triggered by 3.9-7.7 A extension, while natural unfolding rates were of the order of 10(-3) s(-1). Besides acting as individually stable units, helices associated pairwise, establishing a collective potential barrier. The unfolding pathways of individual proteins reflect distinct pulling speed-dependent unfolding routes in their energy landscapes. These observations support the two-stage model of membrane protein folding in which alpha helices insert into the membrane as stable units and then assemble into the functional protein.     

Complex stability of single proteins explored by forced unfolding experiments

In the last decade atomic force microscopy has been used to measure the mechanical stability of single proteins. These force spectroscopy experiments have shown that many water-soluble and membrane proteins unfold via one or more intermediates. Recently, Li and co-workers found a linear correlation between the unfolding force of the native state and the intermediate in fibronectin, which they suggested indicated the presence of a molecular memory or multiple unfolding pathways (1). Here, we apply two independent methods in combination with Monte Carlo simulations to analyze the unfolding of alpha-helices E and D of bacteriorhodopsin (BR). We show that correlation analysis of unfolding forces is very sensitive to errors in force calibration of the instrument. In contrast, a comparison of relative forces provides a robust measure for the stability of unfolding intermediates. The proposed approach detects three energetically different states of alpha-helices E and D in trimeric BR. These states are not observed for monomeric BR and indicate that substantial information is hidden in forced unfolding experiments of single proteins.     

Complex Stability of Single Proteins Explored by Forced Unfolding Experiments

In the last decade atomic force microscopy has been used to measure the mechanical stability of single proteins. These force spectroscopy experiments have shown that many water-soluble and membrane proteins unfold via one or more intermediates. Recently, Li and co-workers found a linear correlation between the unfolding force of the native state and the intermediate in fibronectin, which they suggested indicated the presence of a molecular memory or multiple unfolding pathways (1). Here, we apply two independent methods in combination with Monte Carlo simulations to analyze the unfolding of α-helices E and D of bacteriorhodopsin (BR). We show that correlation analysis of unfolding forces is very sensitive to errors in force calibration of the instrument. In contrast, a comparison of relative forces provides a robust measure for the stability of unfolding intermediates. The proposed approach detects three energetically different states of α-helices E and D in trimeric BR. These states are not observed for monomeric BR and indicate that substantial information is hidden in forced unfolding experiments of single proteins.     

Observing the osmophobic effect in action at the single molecule level

Protecting osmolytes are widespread small organic molecules able to stabilize the folded state of most proteins against various denaturing stresses in vivo. The osmophobic model explains thermodynamically their action through a preferential exclusion of the osmolyte molecules from the protein surface, thus favoring the formation of intrapeptide hydrogen bonds. Few works addressed the influence of protecting osmolytes on the protein unfolding transition state and kinetics. Among those, previous single molecule force spectroscopy experiments evidenced a complexation of the protecting osmolyte molecules at the unfolding transition state of the protein, in apparent contradiction with the osmophobic nature of the protein backbone. We present single-molecule evidence that glycerol, which is a ubiquitous protecting osmolyte, stabilizes a globular protein against mechanical unfolding without binding into its unfolding transition state structure. We show experimentally that glycerol does not change the position of the unfolding transition state as projected onto the mechanical reaction coordinate. Moreover, we compute theoretically the projection of the unfolding transition state onto two other common reaction coordinates, that is, the number of native peptide bonds and the weighted number of native contacts. To that end, we augment an analytic Ising-like protein model with support for group-transfer free energies. Using this model, we find again that the position of the unfolding transition state does not change in the presence of glycerol, giving further support to the conclusions based on the single-molecule experiments.     

Predicting repeat protein folding kinetics from an experimentally determined folding energy landscape

The Notch ankyrin domain is a repeat protein whose folding has been characterized through equilibrium and kinetic measurements. In previous work, equilibrium folding free energies of truncated constructs were used to generate an experimentally determined folding energy landscape (Mello and Barrick, Proc Natl Acad Sci USA 2004;101:14102–14107). Here, this folding energy landscape is used to parameterize a kinetic model in which local transition probabilities between partly folded states are based on energy values from the landscape. The landscape‐based model correctly predicts highly diverse experimentally determined folding kinetics of the Notch ankyrin domain and sequence variants. These predictions include monophasic folding and biphasic unfolding, curvature in the unfolding limb of the chevron plot, population of a transient unfolding intermediate, relative folding rates of 19 variants spanning three orders of magnitude, and a change in the folding pathway that results from C‐terminal stabilization. These findings indicate that the folding pathway(s) of the Notch ankyrin domain are thermodynamically selected: the primary determinants of kinetic behavior can be simply deduced from the local stability of individual repeats.     

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

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

Mechanical unfoldons as building blocks of maltose-binding protein

Identifying independently folding cores or substructures is important for understanding and assaying the structure, function and assembly of large proteins. Here, we suggest mechanical stability as a criterion to identify building blocks of the 366 amino acid maltose-binding protein (MBP). We find that MBP, when pulled at its termini, unfolds via three (meta-) stable unfolding intermediates. Consequently, the MBP structure consists of four structural blocks (unfoldons) that detach sequentially from the folded structure upon force application. We used cysteine cross-link mutations to characterize the four unfoldons structurally. We showed that many MBP constructs composed of those building blocks indeed form stably folded structures in solution. Mechanical unfoldons may provide a new tool for a systematic search for stable substructures of large proteins.     

Mechanically unfolding protein L using a laser-feedback-controlled cantilever

Force spectroscopy using the atomic force microscope (AFM) can yield important information on the strength and lifetimes of the folded states of single proteins and their complexes when they are loaded with force. For example, by mechanically unfolding concatenated proteins at different velocities, a dynamic force spectrum can be built up that allows reconstruction of the energy landscape that the protein traverses during unfolding. To characterize fully the unfolding landscape, however, it is necessary both to explore the entire force spectrum and to characterize each species populated during unfolding. In the conventional AFM apparatus, force is applied to the protein construct through a compliant cantilever. This limits the dynamic range of the force spectrum that can be probed, and the cantilever recoil after unfolding may mask the presence of metastable intermediates. Here, we describe to our knowledge a new technique—constant-deflection AFM—in which the compliance of the AFM cantilever is removed. Using this technique, we show that protein L exhibits a more complex unfolding energy landscape than previously detected using the conventional technique. This technique is also able to detect the presence of a refolding intermediate whose formation is otherwise prevented by cantilever recoil.     

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.     

Evidence for sequential barriers and obligatory intermediates in apparent two-state protein folding

Many small proteins fold fast and without detectable intermediates. This is frequently taken as evidence against the importance of partially folded states, which often transiently accumulate during folding of larger proteins. To get insight into the properties of free energy barriers in protein folding we analyzed experimental data from 23 proteins that were reported to show non-linear activation free-energy relationships. These non-linearities are generally interpreted in terms of broad transition barrier regions with a large number of energetically similar states. Our results argue against the presence of a single broad barrier region. They rather indicate that the non-linearities are caused by sequential folding pathways with consecutive distinct barriers and a few obligatory high-energy intermediates. In contrast to a broad barrier model the sequential model gives a consistent picture of the folding barriers for different variants of the same protein and when folding of a single protein is analyzed under different solvent conditions. The sequential model is also able to explain changes from linear to non-linear free energy relationships and from apparent two-state folding to folding through populated intermediates upon single point mutations or changes in the experimental conditions. These results suggest that the apparent discrepancy between two-state and multi-state folding originates in the relative stability of the intermediates, which argues for the importance of partially folded states in protein folding.     

Structure-Based Derivation of Protein Folding Intermediates and Energies from Optical Tweezers

Optical tweezers (OTs) measure the force-dependent time-resolved extension of a single macromolecule tethered between two trapped beads. From this measurement, it is possible to determine the folding intermediates, energies, and kinetics of the macromolecule. Previous data analysis generally has used the extension as a reaction coordinate to characterize the observed folding transitions. Despite its convenience, the extension poorly describes folding in the absence of force. Here, we chose the contour length of the unfolded polypeptide as a reaction coordinate and modeled the extensions of protein structures along their predicted folding pathways based on high-resolution structures of the proteins in their native states. We included the extension in our model to calculate the total extensions, energies, and transition rates of the proteins as a function of force. We fit these calculations to the corresponding experimental measurements and obtained the best-fit conformations and energies of proteins in different folding states. We applied our method to analyze single-molecule trajectories of two representative protein complexes responsible for membrane fusion, the HIV-1 glycoprotein 41 and the synaptic SNARE proteins, which involved transitions between two and five states, respectively. Nonlinear fitting of the model to the experimental data revealed the structures of folding intermediates and transition states and their associated energies. Our results demonstrate that the contour length is a useful reaction coordinate to characterize protein folding and that intrinsic extensions of protein structures should be taken into account to properly derive the conformations and energies of protein folding intermediates from single-molecule manipulation experiments.     

How cooperative are protein folding and unfolding transitions?

A thermodynamically and kinetically simple picture of protein folding envisages only two states, native (N) and unfolded (U), separated by a single activation free energy barrier, and interconverting by cooperative two‐state transitions. The folding/unfolding transitions of many proteins occur, however, in multiple discrete steps associated with the formation of intermediates, which is indicative of reduced cooperativity. Furthermore, much advancement in experimental and computational approaches has demonstrated entirely non‐cooperative (gradual) transitions via a continuum of states and a multitude of small energetic barriers between the N and U states of some proteins. These findings have been instrumental towards providing a structural rationale for cooperative versus noncooperative transitions, based on the coupling between interaction networks in proteins. The cooperativity inherent in a folding/unfolding reaction appears to be context dependent, and can be tuned via experimental conditions which change the stabilities of N and U. The evolution of cooperativity in protein folding transitions is linked closely to the evolution of function as well as the aggregation propensity of the protein. A large activation energy barrier in a fully cooperative transition can provide the kinetic control required to prevent the accumulation of partially unfolded forms, which may promote aggregation. Nevertheless, increasing evidence for barrier‐less “downhill” folding, as well as for continuous “uphill” unfolding transitions, indicate that gradual non‐cooperative processes may be ubiquitous features on the free energy landscape of protein folding.     

Mechanical unfolding of RNA: from hairpins to structures with internal multiloops

Mechanical unfolding of RNA structures, ranging from hairpins to ribozymes, using laser optical tweezer experiments have begun to reveal the features of the energy landscape that cannot be easily explored using conventional experiments. Upon application of constant force (f), RNA hairpins undergo cooperative transitions from folded to unfolded states whereas subdomains of ribozymes unravel one at a time. Here, we use a self-organized polymer model and Brownian dynamics simulations to probe mechanical unfolding at constant force and constant-loading rate of four RNA structures of varying complexity. For simple hairpins, such as P5GA, application of constant force or constant loading rate results in bistable cooperative transitions between folded and unfolded states without populating any intermediates. The transition state location (DeltaxFTS) changes dramatically as the loading rate is varied. At loading rates comparable to those used in laser optical tweezer experiments, the hairpin is plastic, with DeltaxFTS being midway between folded and unfolded states; whereas at high loading rates, DeltaxFTS moves close to the folded state, i.e., RNA is brittle. For the 29-nucleotide TAR RNA with the three-nucleotide bulge, unfolding occurs in a nearly two-state manner with an occasional pause in a high free energy metastable state. Forced unfolding of the 55 nucleotides of the Hepatitis IRES domain IIa, which has a distorted L-shaped structure, results in well-populated stable intermediates. The most stable force-stabilized intermediate represents straightening of the L-shaped structure. For these structures, the unfolding pathways can be predicted using the contact map of the native structures. Unfolding of a RNA motif with internal multiloop, namely, the 109-nucleotide prohead RNA that is part of the 29 DNA packaging motor, at constant value of rf occurs with three distinct rips that represent unraveling of the paired helices. The rips represent kinetic barriers to unfolding. Our work shows 1), the response of RNA to force is largely determined by the native structure; and 2), only by probing mechanical unfolding over a wide range of forces can the underlying energy landscape be fully explored.     

Protein conformational changes in the bacteriorhodopsin photocycle

We report a comprehensive electron crystallographic analysis of conformational changes in the photocycle of wild-type bacteriorhodopsin and in a variety of mutant proteins with kinetic defects in the photocycle. Specific intermediates that accumulate in the late stages of the photocycle of wild-type bacteriorhodopsin, the single mutants D38R, D96N, D96G, T46V, L93A and F219L, and the triple mutant D96G/F171C/F219L were trapped by freezing two-dimensional crystals in liquid ethane at varying times after illumination with a light flash. Electron diffraction patterns recorded from these crystals were used to construct projection difference Fourier maps at 3.5 A resolution to define light-driven changes in protein conformation.Our experiments demonstrate that in wild-type bacteriorhodopsin, a large protein conformational change occurs within approximately 1 ms after illumination. Analysis of structural changes in wild-type and mutant bacteriorhodopsins under conditions when either the M or the N intermediate is preferentially accumulated reveals that there are only small differences in structure between M and N intermediates trapped in the same protein. However, a considerably larger variation is observed when the same optical intermediate is trapped in different mutants. In some of the mutants, a partial conformational change is present even prior to illumination, with additional changes occurring upon illumination. Selected mutations, such as those in the D96G/F171C/F219L triple mutant, can sufficiently destabilize the wild-type structure to generate almost the full extent of the conformational change in the dark, with minimal additional light-induced changes. We conclude that the differences in structural changes observed in mutants that display long-lived M, N or O intermediates are best described as variations of one fundamental type of conformational change, rather than representing structural changes that are unique to the optical intermediate that is accumulated. Our observations thus support a simplified view of the photocycle of wild-type bacteriorhodopsin in which the structures of the initial state and the early intermediates (K, L and M1) are well approximated by one protein conformation, while the structures of the later intermediates (M2, N and O) are well approximated by the other protein conformation. We propose that in wild-type bacteriorhodopsin and in most mutants, this conformational change between the M1 and M2 states is likely to make an important contribution towards efficiently switching proton accessibility of the Schiff base from the extracellular side to the cytoplasmic side of the membrane.     

Using circular dichroism collected as a function of temperature to determine the thermodynamics of protein unfolding and binding interactions

Circular dichroism (CD) is an excellent spectroscopic technique for following the unfolding and folding of proteins as a function of temperature. One of its principal applications is to determine the effects of mutations and ligands on protein and polypeptide stability. If the change in CD as a function of temperature is reversible, analysis of the data may be used to determined the van't Hoff enthalpy and entropy of unfolding, the midpoint of the unfolding transition and the free energy of unfolding. Binding constants of protein-protein and protein-ligand interactions may also be estimated from the unfolding curves. Analysis of CD spectra obtained as a function of temperature is also useful to determine whether a protein has unfolding intermediates. Measurement of the spectra of five folded proteins and their unfolding curves at a single wavelength requires approximately 8 h.     

Shape of the free energy barriers for protein folding probed by multiple perturbation analysis

The characterization of the free energy barriers has been a major goal in studies on the mechanism of protein folding. Testing the effect of mutations or denaturants on protein folding reactions revealed that transition state movement is rare, suggesting that folding barriers are robust and narrow maxima on the free energy landscape. Here we demonstrate that the application of multiple perturbations allows the observation of small transition state movements that escape detection in single perturbation experiments. We used tendamistat as a model protein to test the broadness of the free energy barriers. Tendamistat folds over two consecutive transition states and through a high-energy intermediate. Measuring the combined effect of temperature and denaturant on the position of the transition state in the wild-type protein and in several mutants revealed that the early transition state shows significant transition state movement. Its accessible surface area state becomes more native-like with destabilization of the native state by temperature. To the same extent, the entropy of the early transition state becomes more native-like with increasing denaturant concentration, in accordance with Hammond behavior. The position of the late transition state, in contrast, is much less sensitive to the applied perturbations. These results suggest that the barriers in protein folding become increasingly narrow as the folding polypeptide chain approaches the native state.