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.     

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.     

Contour length and refolding rate of a small protein controlled by engineered disulfide bonds

The introduction of disulfide bonds into proteins creates additional mechanical barriers and limits the unfolded contour length (i.e., the maximal extension) measured by single-molecule force spectroscopy. Here, we engineer single disulfide bonds into four different locations of the human cardiac titin module (I27) to control the contour length while keeping the distance to the transition state unchanged. This enables the study of several biologically important parameters. First, we are able to precisely determine the end-to-end length of the transition state before unfolding (53 Angstrom), which is longer than the end-to-end length of the protein obtained from NMR spectroscopy (43 Angstrom). Second, the measured contour length per amino acid from five different methods (4.0 +/- 0.2 Angstrom) is longer than the end-to-end length obtained from the crystal structure (3.6 Angstrom). Our measurement of the contour length takes into account all the internal degrees of freedom of the polypeptide chain, whereas crystallography measures the end-to-end length within the "frozen" protein structure. Furthermore, the control of contour length and therefore the number of amino acids unraveled before reaching the disulfide bond (n) facilitates the test of the chain length dependence on the folding time (tau(F)). We find that both a power law scaling tau(F) lambda n(lambda) with lambda = 4.4, and an exponential scaling with n(0.6) fit the data range, in support of different protein-folding scenarios.     

Point mutations in membrane proteins reshape energy landscape and populate different unfolding pathways

Using single-molecule force spectroscopy, we investigated the effect of single point mutations on the energy landscape and unfolding pathways of the transmembrane protein bacteriorhodopsin. We show that the unfolding energy barriers in the energy landscape of the membrane protein followed a simple two-state behavior and represent a manifestation of many converging unfolding pathways. Although the unfolding pathways of wild-type and mutant bacteriorhodopsin did not change, indicating the presence of same ensemble of structural unfolding intermediates, the free energies of the rate-limiting transition states of the bacteriorhodopsin mutants decreased as the distance of those transition states to the folded intermediate states decreased. Thus, all mutants exhibited Hammond behavior and a change in the free energies of the intermediates along the unfolding reaction coordinate and, consequently, their relative occupancies. This is the first experimental proof showing that point mutations can reshape the free energy landscape of a membrane protein and force single proteins to populate certain unfolding pathways over others.     

Mechanics and structure of titin oligomers explored with atomic force microscopy

Titin is a giant polypeptide that spans half of the striated muscle sarcomere and generates passive force upon stretch. To explore the elastic response and structure of single molecules and oligomers of titin, we carried out molecular force spectroscopy and atomic force microscopy (AFM) on purified full-length skeletal-muscle titin. From the force data, apparent persistence lengths as long as approximately 1.5 nm were obtained for the single, unfolded titin molecule. Furthermore, data suggest that titin molecules may globally associate into oligomers which mechanically behave as independent wormlike chains (WLCs). Consistent with this, AFM of surface-adsorbed titin molecules revealed the presence of oligomers. Although oligomers may form globally via head-to-head association of titin, the constituent molecules otherwise appear independent from each other along their contour. Based on the global association but local independence of titin molecules, we discuss a mechanical model of the sarcomere in which titin molecules with different contour lengths, corresponding to different isoforms, are held in a lattice. The net force response of aligned titin molecules is determined by the persistence length of the tandemly arranged, different WLC components of the individual molecules, the ratio of their overall contour lengths, and by domain unfolding events. Biased domain unfolding in mechanically selected constituent molecules may serve as a compensatory mechanism for contour- and persistence-length differences. Variation in the ratio and contour length of the component chains may provide mechanisms for the fine-tuning of the sarcomeric passive force response.     

Mechanics and structure of titin oligomers explored with atomic force microscopy

Titin is a giant polypeptide that spans half of the striated muscle sarcomere and generates passive force upon stretch. To explore the elastic response and structure of single molecules and oligomers of titin, we carried out molecular force spectroscopy and atomic force microscopy (AFM) on purified full-length skeletal-muscle titin. From the force data, apparent persistence lengths as long as ∼1.5 nm were obtained for the single, unfolded titin molecule. Furthermore, data suggest that titin molecules may globally associate into oligomers which mechanically behave as independent wormlike chains (WLCs). Consistent with this, AFM of surface-adsorbed titin molecules revealed the presence of oligomers. Although oligomers may form globally via head-to-head association of titin, the constituent molecules otherwise appear independent from each other along their contour. Based on the global association but local independence of titin molecules, we discuss a mechanical model of the sarcomere in which titin molecules with different contour lengths, corresponding to different isoforms, are held in a lattice. The net force response of aligned titin molecules is determined by the persistence length of the tandemly arranged, different WLC components of the individual molecules, the ratio of their overall contour lengths, and by domain unfolding events. Biased domain unfolding in mechanically selected constituent molecules may serve as a compensatory mechanism for contour- and persistence-length differences. Variation in the ratio and contour length of the component chains may provide mechanisms for the fine-tuning of the sarcomeric passive force response.     

Solvation and desolvation effects in protein folding: native flexibility, kinetic cooperativity and enthalpic barriers under isostability conditions

As different parts of a protein chain approach one another during folding, they are expected to encounter desolvation barriers before optimal packing is achieved. This impediment originates from the water molecule's finite size, which entails a net energetic cost for water exclusion when the formation of compensating close intraprotein contacts is not yet complete. Based on recent advances, we extend our exploration of these microscopic elementary desolvation barriers' roles in the emergence of generic properties of protein folding. Using continuum Gō-like C(alpha) chain models of chymotrypsin inhibitor 2 (CI2) and barnase as examples, we underscore that elementary desolvation barriers between a protein's constituent groups can significantly reduce native conformational fluctuations relative to model predictions that neglected these barriers. An increasing height of elementary desolvation barriers leads to thermodynamically more cooperative folding/unfolding transitions (i.e., higher overall empirical folding barriers) and higher degrees of kinetic cooperativity as manifested by more linear rate-stability relationships under constant temperature. Applying a spatially non-uniform thermodynamic parametrization we recently introduced for the pairwise C(alpha) potentials of mean force, the present barnase model further illustrates that desolvation is a probable physical underpinning for the experimentally observed high intrinsic enthalpic folding barrier under isostability conditions.     

A novel pattern recognition algorithm to classify membrane protein unfolding pathways with high-throughput single-molecule force spectroscopy

MOTIVATION: Misfolding of membrane proteins plays an important role in many human diseases such as retinitis pigmentosa, hereditary deafness and diabetes insipidus. Little is known about membrane proteins as there are only very few high-resolution structures. Single-molecule force spectroscopy is a novel technique, which measures the force necessary to pull a protein out of a membrane. Such force curves contain valuable information on the protein structure, conformation, and inter- and intra-molecular forces. High-throughput force spectroscopy experiments generate hundreds of force curves including spurious ones and good curves, which correspond to different unfolding pathways. Manual analysis of these data is a bottleneck and source of inconsistent and subjective annotation. RESULTS: We propose a novel algorithm for the identification of spurious curves and curves representing different unfolding pathways. Our algorithm proceeds in three stages: first, we reduce noise in the curves by applying dimension reduction; second, we align the curves with dynamic programming and compute pairwise distances and third, we cluster the curves based on these distances. We apply our method to a hand-curated dataset of 135 force curves of bacteriorhodopsin mutant P50A. Our algorithm achieves a success rate of 81% distinguishing spurious from good curves and a success rate of 76% classifying unfolding pathways. As a result, we discuss five different unfolding pathways of bacteriorhodopsin including three main unfolding events and several minor ones. Finally, we link folding barriers to the degree of conservation of residues. Overall, the algorithm tackles the force spectroscopy bottleneck and leads to more consistent and reproducible results paving the way for high-throughput analysis of structural features of membrane proteins.     

Single molecule measurements of titin elasticity

Titin, with a massive single chain of 3--4MDa and multiple modular motifs, spans the half-sarcomere of skeletal and cardiac muscles and serves important, multifaceted functions. In recent years, titin has become a favored subject of single molecule observations by atomic force microscopy (AFM) and laser optical trap (LOT). Here we review these single titin molecule extension studies with an emphasis on understanding their relevance to titin elasticity in muscle function. Some fundamental aspects of the methods for single titin molecule investigations, including the application of dynamic force, the elasticity models for filamentous titin motifs, the technical foundations and calibrations of AFM and LOT, and titin sample preparations are provided. A chronological review of major publications on recent single titin extension observations is presented. This is followed by summary evaluations of titin domain folding/unfolding results and of elastic properties of filamentous titin motifs. Implications of these single titin measurements for muscle physiology/pathology are discussed and forthcoming advances in single titin studies are anticipated.     

The key event in force-induced unfolding of Titin's immunoglobulin domains

Steered molecular dynamics simulation of force-induced titin immunoglobulin domain I27 unfolding led to the discovery of a significant potential energy barrier at an extension of approximately 14 A on the unfolding pathway that protects the domain against stretching. Previous simulations showed that this barrier is due to the concurrent breaking of six interstrand hydrogen bonds (H-bonds) between beta-strands A' and G that is preceded by the breaking of two to three hydrogen bonds between strands A and B, the latter leading to an unfolding intermediate. The simulation results are supported by Angstrom-resolution atomic force microscopy data. Here we perform a structural and energetic analysis of the H-bonds breaking. It is confirmed that H-bonds between strands A and B break rapidly. However, the breaking of the H-bond between strands A' and G needs to be assisted by fluctuations of water molecules. In nanosecond simulations, water molecules are found to repeatedly interact with the protein backbone atoms, weakening individual interstrand H-bonds until all six A'-G H-bonds break simultaneously under the influence of external stretching forces. Only when those bonds are broken can the generic unfolding take place, which involves hydrophobic interactions of the protein core and exerts weaker resistance against stretching than the key event.     

Reference-free alignment and sorting of single-molecule force spectroscopy data

Single-molecule force spectroscopy has become a versatile tool for investigating the (un)folding of proteins and other polymeric molecules. Like other single-molecule techniques, single-molecule force spectroscopy requires recording and analysis of large data sets to extract statistically meaningful conclusions. Here, we present a data analysis tool that provides efficient filtering of heterogeneous data sets, brings spectra into register based on a reference-free alignment algorithm, and determines automatically the location of unfolding barriers. Furthermore, it groups spectra according to the number of unfolding events, subclassifies the spectra using cross correlation-based sorting, and extracts unfolding pathways by principal component analysis and clustering methods to extracted peak positions. Our approach has been tested on a data set obtained through mechanical unfolding of bacteriorhodopsin (bR), which contained a significant number of spectra that did not show the well-known bR fingerprint. In addition, we have tested the performance of the data analysis tool on unfolding data of the soluble multidomain (Ig27)(8) protein.     

A kinetic molecular model of the reversible unfolding and refolding of titin under force extension.

Molecular elasticity is a physicomechanical property that is associated with a select number of polypeptides and proteins, such as the giant muscle protein, titin, and the extracellular matrix protein, tenascin. Both proteins have been the subject of atomic force microscopy (AFM), laser tweezer, and other in vitro methods for examining the effects of force extension on the globular (FNIII/Ig-like) domains that comprise each protein. In this report we present a time-dependent method for simulating AFM force extension and its effect on FNIII/Ig domain unfolding and refolding. This method treats the unfolding and refolding process as a standard three-state protein folding model (U right arrow over left arrow T right arrow over left arrow F, where U is the unfolded state, T is the transition or intermediate state, and F is the fully folded state), and integrates this approach within the wormlike chain (WLC) concept. We simulated the effect of AFM tip extension on a hypothetical titin molecule comprised of 30 globular domains (Ig or FNIII) and 25% Pro-Glu-Val-Lys (PEVK) content, and analyzed the unfolding and refolding processes as a function of AFM tip extension, extension rate, and variation in PEVK content. In general, we find that the use of a three-state protein-folding kinetic-based model and the implicit inclusion of PEVK domains can accurately reproduce the experimental force-extension curves observed for both titin and tenascin proteins. Furthermore, our simulation data indicate that PEVK domains exhibit extensibility behavior, assist in the unfolding and refolding of FNIII/Ig domains in the titin molecule, and act as a force "buffer" for the FNIII/Ig domains, particularly at low and moderate extension forces.     

Unfolding barriers in bacteriorhodopsin probed from the cytoplasmic and the extracellular side by AFM

Selecting an individual membrane protein and probing its mechanical properties has become possible by AFM-based single-molecule force spectroscopy. In contrast to earlier studies, we extracted and unfolded bacteriorhodopsin monomers from the purple membrane not only from the cytoplasmic side, but also from the extracellular side, and recorded the force extension profiles. This way different pathways through the potential landscape are explored. A map of the 21 most dominant barriers with their positions relative to the amino acid sequences is given at an accuracy of +/-3 aa. Most barriers were found to provide resistance to forced unfolding only when extracted toward one of the sides. However, certain barriers have identical positions to within a few amino acids when probed from either of the sides, which typifies them as structural traps.     

Individual Globular Domains and Domain Unfolding Visualized in Overstretched Titin Molecules with Atomic Force Microscopy

Titin is a giant elastomeric protein responsible for the generation of passive muscle force. Mechanical force unfolds titin’s globular domains, but the exact structure of the overstretched titin molecule is not known. Here we analyzed, by using high-resolution atomic force microscopy, the structure of titin molecules overstretched with receding meniscus. The axial contour of the molecules was interrupted by topographical gaps with a mean width of 27.7 nm that corresponds well to the length of an unfolded globular (immunoglobulin and fibronectin) domain. The wide gap-width distribution suggests, however, that additional mechanisms such as partial domain unfolding and the unfolding of neighboring domain multimers may also be present. In the folded regions we resolved globules with an average spacing of 5.9 nm, which is consistent with a titin chain composed globular domains with extended interdomain linker regions. Topographical analysis allowed us to allocate the most distal unfolded titin region to the kinase domain, suggesting that this domain systematically unfolds when the molecule is exposed to overstretching forces. The observations support the prediction that upon the action of stretching forces the N-terminal ß-sheet of the titin kinase unfolds, thus exposing the enzyme’s ATP-binding site and hence contributing to the molecule’s mechanosensory function.     

Probing origins of molecular interactions stabilizing the membrane proteins halorhodopsin and bacteriorhodopsin

Single-molecule atomic force microscopy and spectroscopy were applied to detect molecular interactions stabilizing the structure of halorhodopsin (HR), a light-driven chloride pump from Halobacterium salinarum. Because of the high structural and sequence similarities between HR and bacteriorhodopsin, we compared their unfolding pathways and polypeptide regions that established structurally stable segments against unfolding. Unfolding pathways and structural segments stabilizing the proteins both exhibited a remarkably high similarity. This suggests that different amino acid compositions can establish structurally indistinguishable energetic barriers. These stabilizing domains rather result from comprehensive interactions of all amino acids within a structural region than from specific interactions. However, one additional unfolding barrier located within a short segment of helix E was detected for HR. This barrier correlated with a Pi-bulk interaction, which locally disrupts helix E and divides a structural stabilizing segment.     

Hopping around an entropic barrier created by force

We use Langevin dynamics to investigate the role played by the recently discovered force-induced entropic energy barrier on the two-state hopping phenomena that has been observed in single RNA, DNA and protein molecules placed under a stretching force. Simple considerations about the free energy of a molecule readily show that the application of force introduces an entropic barrier separating the collapsed state of the molecule, from a force-driven extended conformation. A notable characteristic of the force induced barrier is its long distances to transition state, up to tens of nanometers, which renders the kinetics of crossing this barrier highly sensitive to an applied force. Langevin dynamics across such force induced barriers readily demonstrates the hopping behavior observed for a variety of single molecules placed under force. Such hopping is frequently interpreted as a manifestation of two-state folding/unfolding reactions observed in bulk experiments. However, given that such barriers do not exist at zero force these reactions do not take place at all in bulk.     

Single-molecule force spectroscopy reveals a stepwise unfolding of Caenorhabditis elegans giant protein kinase domains

Myofibril assembly and disassembly are complex processes that regulate overall muscle mass. Titin kinase has been implicated as an initiating catalyst in signaling pathways that ultimately result in myofibril growth. In titin, the kinase domain is in an ideal position to sense mechanical strain that occurs during muscle activity. The enzyme is negatively regulated by intramolecular interactions occurring between the kinase catalytic core and autoinhibitory/regulatory region. Molecular dynamics simulations suggest that human titin kinase acts as a force sensor. However, the precise mechanism(s) resulting in the conformational changes that relieve the kinase of this autoinhibition are unknown. Here we measured the mechanical properties of the kinase domain and flanking Ig/Fn domains of the Caenorhabditis elegans titin-like proteins twitchin and TTN-1 using single-molecule atomic force microscopy. Our results show that these kinase domains have significant mechanical resistance, unfolding at forces similar to those for Ig/Fn β-sandwich domains (30-150 pN). Further, our atomic force microscopy data is consistent with molecular dynamic simulations, which show that these kinases unfold in a stepwise fashion, first an unwinding of the autoinhibitory region, followed by a two-step unfolding of the catalytic core. These data support the hypothesis that titin kinase may function as an effective force sensor.     

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 Position-Dependent Diffusion in Folding Reactions Using Single-Molecule Force Spectroscopy

Folding of proteins and nucleic acids involves a diffusive search over a multidimensional conformational energy landscape for the minimal-energy structure. When examining the projection of conformational motions onto a one-dimensional reaction coordinate, as done in most experiments, the diffusion coefficient D is generally position dependent. However, it has proven challenging to measure such position-dependence experimentally. We investigated the position-dependence of D in the folding of DNA hairpins as a simple model system in two ways: first, by analyzing the round-trip time to return to a given extension in constant-force extension trajectories measured by force spectroscopy, and second, by analyzing the fall time required to reach a given extension in force jump measurements. These methods yielded conflicting results: the fall time implied a fairly constant D, but the round-trip time implied variations of over an order of magnitude. Comparison of experiments with computational simulations revealed that both methods were strongly affected by experimental artifacts inherent to force spectroscopy measurements, which obscured the intrinsic position-dependence of D. Lastly, we applied Kramers’s theory to the kinetics of hairpins with energy barriers located at different positions along the hairpin stem, as a crude probe of D at different stem positions, and we found that D did not vary much as the barrier was moved along the reaction coordinate. This work underlines the difficulties faced when trying to deduce position-dependent diffusion coefficients from experimental folding trajectories.