Probing the Energy Landscape of the Membrane Protein Bacteriorhodopsin

Mechanistic Insight into the Mechanical Unfolding of the Integral Membrane Diacylglycerol Kinase

The essential forces stabilizing membrane proteins and governing their folding and unfolding are difficult to decipher. Single-molecule atomic force spectroscopy mechanically unfolds individual membrane proteins and quantifies their dynamics and energetics. However, it remains challenging to structurally assign unfolding intermediates precisely and to deduce dominant interactions between specific residues that facilitate either the localized stabilization of these intermediates or the global assembly of membrane proteins. Here, we performed force spectroscopy experiments and multiscale molecular dynamics simulations to study the unfolding pathway of diacylglycerol kinase (DGK), a small trimeric multispan transmembrane enzyme. The remarkable agreement between experiments and simulations allowed precise structural assignment and interaction analysis of unfolding intermediates, bypassing existing limitations on structural mapping, and thus provided mechanistic explanations for the formation of these states. DGK unfolding was found to proceed with structural segments varying in size that do not correlate with its secondary structure. We identified intermolecular side-chain packing interactions as one of the major contributions to the stability of unfolding intermediates. Mutagenesis creating packing defects induced a dramatic decrease in the mechano-stability of corresponding intermediates and also in the thermo-stability of DGK trimer, in good agreement with predictions from simulations. Hence, the molecular determinants of the mechano- and thermo-stability of a membrane protein can be identified at residue resolution. The accurate structural assignment established and microscopic mechanism revealed in this work may substantially expand the scope of single-molecule studies of membrane proteins.

Single-Molecule Force Spectroscopy of Membrane Protein Folding

Single-molecule force spectroscopy is a unique method that can probe the structural changes of single proteins at a high spatiotemporal resolution while mechanically manipulating them over a wide force range. Here, we review the current understanding of membrane protein folding learned by using the force spectroscopy approach. Membrane protein folding in lipid bilayers is one of the most complex biological processes in which diverse lipid molecules and chaperone proteins are intricately involved. The approach of single protein forced unfolding in lipid bilayers has produced important findings and insights into membrane protein folding. This review provides an overview of the forced unfolding approach, including recent achievements and technical advances. Progress in the methods can reveal more interesting cases of membrane protein folding and clarify general mechanisms and principles.

Free-energy changes of bacteriorhodopsin point mutants measured by single-molecule force spectroscopy

Single amino acid mutations provide quantitative insight into the energetics that underlie the dynamics and folding of membrane proteins. Chemical denaturation is the most widely used assay and yields the change in unfolding free energy (ΔΔG). It has been applied to >80 different residues of bacteriorhodopsin (bR), a model membrane protein. However, such experiments have several key limitations: 1) a nonnative lipid environment, 2) a denatured state with significant secondary structure, 3) error introduced by extrapolation to zero denaturant, and 4) the requirement of globally reversible refolding. We overcame these limitations by reversibly unfolding local regions of an individual protein with mechanical force using an atomic-force-microscope assay optimized for 2 μs time resolution and 1 pN force stability. In this assay, bR was unfolded from its native bilayer into a well-defined, stretched state. To measure ΔΔG, we introduced two alanine point mutations into an 8-amino-acid region at the C-terminal end of bR's G helix. For each, we reversibly unfolded and refolded this region hundreds of times while the rest of the protein remained folded. Our single-molecule-derived ΔΔG for mutant L223A (-2.3 ± 0.6 kcal/mol) quantitatively agreed with past chemical denaturation results while our ΔΔG for mutant V217A was 2.2-fold larger (-2.4 ± 0.6 kcal/mol). We attribute the latter result, in part, to contact between Val²¹⁷ and a natively bound squalene lipid, highlighting the contribution of membrane protein-lipid contacts not present in chemical denaturation assays. More generally, we established a platform for determining ΔΔG for a fully folded membrane protein embedded in its native bilayer.

Phosphorylation of the Chaperone-Like HspB5 Rescues Trafficking and Function of F508del-CFTR

Cystic Fibrosis is a lethal monogenic autosomal recessive disease linked to mutations in Cystic Fibrosis Transmembrane Conductance Regulator (CFTR) protein. The most frequent mutation is the deletion of phenylalanine at position 508 of the protein. This F508del-CFTR mutation leads to misfolded protein that is detected by the quality control machinery within the endoplasmic reticulum and targeted for destruction by the proteasome. Modulating quality control proteins as molecular chaperones is a promising strategy for attenuating the degradation and stabilizing the mutant CFTR at the plasma membrane. Among the molecular chaperones, the small heat shock protein HspB1 and HspB4 were shown to promote degradation of F508del-CFTR. Here, we investigated the impact of HspB5 expression and phosphorylation on transport to the plasma membrane, function and stability of F508del-CFTR. We show that a phosphomimetic form of HspB5 increases the transport to the plasma membrane, function and stability of F508del-CFTR. These activities are further enhanced in presence of therapeutic drugs currently used for the treatment of cystic fibrosis (VX-770/Ivacaftor, VX-770+VX-809/Orkambi). Overall, this study highlights the beneficial effects of a phosphorylated form of HspB5 on F508del-CFTR rescue and its therapeutic potential in cystic fibrosis.

On the Interpretation of Force-Induced Unfolding Studies of Membrane Proteins Using Fast Simulations

Single-molecule force spectroscopy has proven extremely beneficial in elucidating folding pathways for membrane proteins. Here, we simulate these measurements, conducting hundreds of unfolding trajectories using our fast Upside algorithm for slow enough speeds to reproduce key experimental features that may be missed using all-atom methods. The speed also enables us to determine the logarithmic dependence of pulling velocities on the rupture levels to better compare to experimental values. For simulations of atomic force microscope measurements in which force is applied vertically to the C-terminus of bacteriorhodopsin, we reproduce the major experimental features including even the back-and-forth unfolding of single helical turns. When pulling laterally on GlpG to mimic the experiment, we observe quite different behavior depending on the stiffness of the spring. With a soft spring, as used in the experimental studies with magnetic tweezers, the force remains nearly constant after the initial unfolding event, and a few pathways and a high degree of cooperativity are observed in both the experiment and simulation. With a stiff spring, however, the force drops to near zero after each major unfolding event, and numerous intermediates are observed along a wide variety of pathways. Hence, the mode of force application significantly alters the perception of the folding landscape, including the number of intermediates and the degree of folding cooperativity, important issues that should be considered when designing experiments and interpreting unfolding data.

Forced Unfolding Mechanism of Bacteriorhodopsin as Revealed by Coarse-Grained Molecular Dynamics

Developments in atomic force microscopy have opened up a new path toward single-molecular phenomena; in particular, during the process of pulling a membrane protein out of a lipid bilayer. However, the characteristic features of the force-distance (F-D) curve of a bacteriorhodopsin in purple membrane, for instance, have not yet been fully elucidated in terms of physicochemical principles. To address the issue, we performed a computer simulation of bacteriorhodopsin with, to our knowledge, a novel coarse-grained (C-G) model. Peptide planes are represented as rigid spheres, while the surrounding environment consisting of water solvents and lipid bilayers is represented as an implicit continuum. Force-field parameters were determined on the basis of auxiliary simulations and experimental values of transfer free energy of each amino acid from water to membrane. According to Popot's two-stage model, we separated molecular interactions involving membrane proteins into two parts: I) affinity of each amino acid to the membrane and intrahelical hydrogen bonding between main chain peptide bonds; and II) interhelix interactions. Then, only part I was incorporated into the C-G model because we assumed that the part plays a dominant role in the forced unfolding process. As a result, the C-G simulation has successfully reproduced the key features, including peak positions, of the experimental F-D curves in the literature, indicating that the peak positions are essentially determined by the residue-lipid and intrahelix interactions. Furthermore, we investigated the relationships between the energy barrier formation on the forced unfolding pathways and the force peaks of the F-D curves.

Applications of Single-Molecule Methods to Membrane Protein Folding Studies

Protein folding is a fundamental life process with many implications throughout biology and medicine. Consequently, there have been enormous efforts to understand how proteins fold. Almost all of this effort has focused on water-soluble proteins, however, leaving membrane proteins largely wandering in the wilderness. The neglect has occurred not because membrane proteins are unimportant but rather because they present many theoretical and technical complications. Indeed, quantitative membrane protein folding studies are generally restricted to a handful of well-behaved proteins. Single-molecule methods may greatly alter this picture, however, because the ability to work at or near infinite dilution removes aggregation problems, one of the main technical challenges of membrane protein folding studies.

Single-molecule force spectroscopy of membrane proteins from membranes freely spanning across nanoscopic pores

Single-molecule force spectroscopy (SMFS) provides detailed insight into the mechanical (un)folding pathways and structural stability of membrane proteins. So far, SMFS could only be applied to membrane proteins embedded in native or synthetic membranes adsorbed to solid supports. This adsorption causes experimental limitations and raises the question to what extent the support influences the results obtained by SMFS. Therefore, we introduce here SMFS from native purple membrane freely spanning across nanopores. We show that correct analysis of the SMFS data requires extending the worm-like chain model, which describes the mechanical stretching of a polypeptide, by the cubic extension model, which describes the bending of a purple membrane exposed to mechanical stress. This new experimental and theoretical approach allows to characterize the stepwise (un)folding of the membrane protein bacteriorhodopsin and to assign the stability of single and grouped secondary structures. The (un)folding and stability of bacteriorhodopsin shows no significant difference between freely spanning and directly supported purple membranes. Importantly, the novel experimental SMFS setup opens an avenue to characterize any protein from freely spanning cellular or synthetic membranes.

The cytosolic chaperone α-crystallin B rescues folding and compartmentalization of misfolded multispan transmembrane proteins

The α-crystallin B chain (CRYAB or HspB5) is a cytosolic chaperone belonging to the small heat shock protein family, which is known to help in the folding of cytosolic proteins. Here we show that CRYAB binds the mutant form of at least two multispan transmembrane proteins (TMPs), exerting an anti-aggregation activity. It rescues the folding of mutant Frizzled4, which is responsible for a rare autosomal dominant form of familial exudative vitreoretinopathy (Fz4-FEVR), and the mutant ATP7B Cu transporter (ATP7B-H1069Q) associated with a common form of Wilson's disease. In the case of Fz4-FEVR, CRYAB prevents the formation of inter-chain disulfide bridges between the lumenal ectodomains of the aggregated mutant chains, which enables correct folding and promotes appropriate compartmentalization on the plasma membrane. ATP7B-H1069Q, with help from CRYAB, folds into the proper conformation, moves to the Golgi complex, and responds to copper overload in the same manner as wild-type ATP7B. These findings strongly suggest that CRYAB plays a pivotal role, previously undetected, in the folding of multispan TMPs and, from the cytosol, is able to orchestrate folding events that take place in the lumen of the ER. Our results contribute to the explanation of the complex scenario behind multispan TMP folding; additionally, they serve to expose interesting avenues for novel therapeutic approaches.

Atomic force microscopy and spectroscopy to probe single membrane proteins in lipid bilayers

The atomic force microscope (AFM) has opened vast avenues hitherto inaccessible to the biological scientist. The high temporal (millisecond) and spatial (nanometer) resolutions of the AFM are suited for studying many biological processes in their native conditions. The AFM cantilever stylus is aptly termed as a "lab on a tip" owing to its versatility as an imaging tool as well as a handle to manipulate single bonds and proteins. Recent examples assert that the AFM can be used to study the mechanical properties and monitor processes of single proteins and single cells, thus affording insight into important mechanistic details. This chapter specifically focuses on practical and analytical protocols of single-molecule AFM methodologies related to high-resolution imaging and single-molecule force spectroscopy of membrane proteins. Both these techniques are operator oriented, and require specialized working knowledge of the instrument, theoretical, and practical skills.

Cholesterol increases kinetic, energetic, and mechanical stability of the human β2-adrenergic receptor

The steroid cholesterol is an essential component of eukaryotic membranes, and it functionally modulates membrane proteins, including G protein-coupled receptors. To reveal insight into how cholesterol modulates G protein-coupled receptors, we have used dynamic single-molecule force spectroscopy to quantify the mechanical strength and flexibility, conformational variability, and kinetic and energetic stability of structural segments stabilizing the human β(2)-adrenergic receptor (β(2)AR) in the absence and presence of the cholesterol analog cholesteryl hemisuccinate (CHS). CHS considerably increased the kinetic, energetic, and mechanical stability of almost every structural segment at sufficient magnitude to alter the structure and functional relationship of β(2)AR. One exception was the structural core segment of β(2)AR, which establishes multiple ligand binding sites, and its properties were not significantly influenced by CHS.

Ligand-specific interactions modulate kinetic, energetic, and mechanical properties of the human β2 adrenergic receptor

G protein-coupled receptors (GPCRs) are a class of versatile proteins that transduce signals across membranes. Extracellular stimuli induce inter- and intramolecular interactions that change the functional state of GPCRs and activate intracellular messenger molecules. How these interactions are established and how they modulate the functional state of GPCRs remain to be understood. We used dynamic single-molecule force spectroscopy to investigate how ligand binding modulates the energy landscape of the human β2 adrenergic receptor (β2 AR). Five different ligands representing either agonists, inverse agonists or neutral antagonists established a complex network of interactions that tuned the kinetic, energetic, and mechanical properties of functionally important structural regions of β2 AR. These interactions were specific to the efficacy profile of the ligands investigated and suggest that the functional modulation of GPCRs follows structurally well-defined interaction patterns.

Single-molecule force spectroscopy from nanodiscs: an assay to quantify folding, stability, and interactions of native membrane proteins

Single-molecule force spectroscopy (SMFS) can quantify and localize inter- and intramolecular interactions that determine the folding, stability, and functional state of membrane proteins. To conduct SMFS the membranes embedding the membrane proteins must be imaged and localized in a rather time-consuming manner. Toward simplifying the investigation of membrane proteins by SMFS, we reconstituted the light-driven proton pump bacteriorhodopsin into lipid nanodiscs. The advantage of using nanodiscs is that membrane proteins can be handled like water-soluble proteins and characterized with similar ease. SMFS characterization of bacteriorhodopsin in native purple membranes and in nanodiscs reveals no significant alterations of structure, function, unfolding intermediates, and strengths of inter- and intramolecular interactions. This demonstrates that lipid nanodiscs provide a unique approach for in vitro studies of native membrane proteins using SMFS and open an avenue to characterize membrane proteins by a wide variety of SMFS approaches that have been established on water-soluble proteins.

Locating an extracellular K+-dependent interaction site that modulates betaine-binding of the Na+-coupled betaine symporter BetP

BetP, a trimeric Na(+)-coupled betaine symporter, senses hyperosmotic stress via its cytoplasmic C-terminal domain and regulates transport activity in dependence of the cytoplasmic K(+)-concentration. This transport regulation of BetP depends on a sophisticated interaction network. Using single-molecule force spectroscopy we structurally localize and quantify these interactions changing on K(+)-dependent transport activation and substrate-binding. K(+) significantly strengthened all interactions, modulated lifetimes of functionally important structural regions, and increased the mechanical rigidity of the symporter. Substrate-binding could modulate, but not establish most of these K(+)-dependent interactions. A pronounced effect triggered by K(+) was observed at the periplasmic helical loop EH2. Tryptophan quenching experiments revealed that elevated K(+)-concentrations akin to those BetP encounters during hyperosmotic stress trigger the formation of a periplasmic second betaine-binding (S2) site, which was found to be at a similar position reported previously for the BetP homologue CaiT. In BetP, the presence of the S2 site strengthened the interaction between EH2, transmembrane α-helix 12 and the K(+)-sensing C-terminal domain resulting in a K(+)-dependent cooperative betaine-binding.

Efficient unfolding pattern recognition in single molecule force spectroscopy data

Background

Single-molecule force spectroscopy (SMFS) is a technique that measures the force necessary to unfold a protein. SMFS experiments generate Force-Distance (F-D) curves. A statistical analysis of a set of F-D curves reveals different unfolding pathways. Information on protein structure, conformation, functional states, and inter- and intra-molecular interactions can be derived.

Results

In the present work, we propose a pattern recognition algorithm and apply our algorithm to datasets from SMFS experiments on the membrane protein bacterioRhodopsin (bR). We discuss the unfolding pathways found in bR, which are characterised by main peaks and side peaks. A main peak is the result of the pairwise unfolding of the transmembrane helices. In contrast, a side peak is an unfolding event in the alpha-helix or other secondary structural element. The algorithm is capable of detecting side peaks along with main peaks.

Therefore, we can detect the individual unfolding pathway as the sequence of events labeled with their occurrences and co-occurrences special to bR's unfolding pathway. We find that side peaks do not co-occur with one another in curves as frequently as main peaks do, which may imply a synergistic effect occurring between helices. While main peaks co-occur as pairs in at least 50% of curves, the side peaks co-occur with one another in less than 10% of curves. Moreover, the algorithm runtime scales well as the dataset size increases.

Conclusions

Our algorithm satisfies the requirements of an automated methodology that combines high accuracy with efficiency in analyzing SMFS datasets. The algorithm tackles the force spectroscopy analysis bottleneck leading to more consistent and reproducible results.

Velocity-dependent mechanical unfolding of bacteriorhodopsin is governed by a dynamic interaction network

Bacteriorhodopsin is a model system for membrane proteins. This seven transmembrane helical protein is embedded within a membrane structure called purple membrane. Its structural stability against mechanical stress was recently investigated by atomic force microscopy experiments, in which single proteins were extracted from the purple membrane. Here, we study this process by all-atom molecular dynamics simulations, in which single bacteriorhodopsin molecules were extracted and unfolded from an atomistic purple membrane model. In our simulations, key features from the experiments like force profiles and location of key residues that resist mechanical unfolding were reproduced. These key residues were seen to be stabilized by a dynamic network of intramolecular interactions. Further, the unfolding pathway was found to be velocity-dependent. Simulations in which the mechanical stress was released during unfolding revealed relaxation motions that allowed characterization of the nonequilibrium processes during fast extraction.

Adhesive properties of Staphylococcus epidermidis probed by atomic force microscopy

Mapping of the surface properties of Staphylococcus epidermidis and of biofilm forming bacteria in general is a key to understand their functions, particularly their adhesive properties. To gain a comprehensive view of the structural and chemical properties of S. epidermidis, four different strains (biofilm positive and biofilm negative strains) were analyzed using in situ atomic force microscopy (AFM). Force measurements performed using bare hydrophilic silicon nitride tips disclosed similar adhesive properties for each strain. However, use of hydrophobic tips showed that hydrophobic forces are not the driving forces for adhesion of the four strains. Rather, the observation of sawtooth force-distance patterns on the surface of biofilm positive strains documents the presence of modular proteins such as Aap that may mediate cell adhesion. Treatment of two biofilm positive strains with two chemical inhibitor compounds leads to a loss of adhesion, suggesting that AFM could be a valuable tool to screen for anti-adhesion molecules.

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.

Nanomechanics of full-length nebulin: an elastic strain gauge in the skeletal muscle sarcomere

Nebulin, a family of giant modular proteins (MW 700-800 kDa), acts as a F-actin thin filament ruler and calcium-linked regulator of actomyosin interaction. The nanomechanics of full length, native rabbit nebulin was investigated with an atomic force microscope by tethering, bracketing, and stretching full-length molecules via pairs of site-specific antibodies that were attached covalently, one to a protein resistant self-assembled monolayer of oligoethylene glycol and the other to the cantilever. Using this new nanomechanics platform that enables the identification of single molecule events via an unbiased analysis of detachment force and distance of all force curves, we showed that nebulin is elastic and extends to approximately 1 microm by external force up to an antibody detachment force of approximately 300-400 pN. Upon stretching, nebulin unravels and yields force spectra with craggy mountain range profiles with variable numbers and heights of force peaks. The peak spacings, analyzed by the model-independent, empirical Hilbert-Huang transform method, displayed underlying periodicities at approximately 15 and approximately 22 nm that may result from the unfolding of one or more nebulin modules between force peaks. Nebulin may act as an elastic strain gauge that interacts optimally with actin only under appropriate strain and stress. This stretch to match protein ruler may also exert a compressive force that stabilizes thin filaments against stress during contraction. We propose that the elasticity of nebulin is integral and essential in the muscle sarcomere.

Substrate binding tunes conformational flexibility and kinetic stability of an amino acid antiporter

We used single molecule dynamic force spectroscopy to unfold individual serine/threonine antiporters SteT from Bacillus subtilis. The unfolding force patterns revealed interactions and energy barriers that stabilized structural segments of SteT. Substrate binding did not establish strong localized interactions but appeared to be facilitated by the formation of weak interactions with several structural segments. Upon substrate binding, all energy barriers of the antiporter changed thereby describing the transition from brittle mechanical properties of SteT in the unbound state to structurally flexible conformations in the substrate-bound state. The lifetime of the unbound state was much shorter than that of the substrate-bound state. This leads to the conclusion that the unbound state of SteT shows a reduced conformational flexibility to facilitate specific substrate binding and a reduced kinetic stability to enable rapid switching to the bound state. In contrast, the bound state of SteT showed an increased conformational flexibility and kinetic stability such as required to enable transport of substrate across the cell membrane. This result supports the working model of antiporters in which alternate substrate access from one to the other membrane surface occurs in the substrate-bound state.

A new atomic force microscope force ramp technique using digital force feedback control reveals mechanically weak protein unfolding events

We have developed a new force ramp modification of the atomic force microscope (AFM) which can control multiple unfolding events of a multi-modular protein using software-based digital force feedback control. With this feedback the force loading rate can be kept constant regardless the length of soft elastic linkage or number of unfolded polypeptide domains. An unfolding event is detected as a sudden drop in force, immediately after which the feedback control reduces the applied force to a low value of a few pN by lowering the force set point. Hence the remaining folded domains can relax and the subsequent force ramp is applied to relaxed protein domains identically in each case. We have applied this technique to determine the kinetic parameters x(u), which is the distance between the native state and transition state, and α(0), which is the unfolding rate constant at zero force, for the mechanical unfolding of a pentamer of I27 domains of titin. In each force ramp the unfolding probability depends on the number of folded domains remaining in the system and we had to take account of this effect in the analysis of unfolding force data. We obtained values of x(u) and α(0) to be 0.28 nm and 1.02 × 10(-3) s(-1), which are in good agreement with those obtained from conventional constant velocity experiments. This method reveals unfolding data at low forces that are not seen in constant velocity experiments and corrects for the change in stiffness that occurs with most mechanical systems throughout the unfolding process to allow constant force ramp experiments to be carried out. In addition, a mechanically weak structure was detected, which formed from the fully extended polypeptide chain during a force quench. This indicates that the new technique will allow studies of the folding kinetics of previously hidden, mechanically weak species.

High-throughput single-molecule force spectroscopy for membrane proteins

Atomic force microscopy-based single-molecule force spectroscopy (SMFS) is a powerful tool for studying the mechanical properties, intermolecular and intramolecular interactions, unfolding pathways, and energy landscapes of membrane proteins. One limiting factor for the large-scale applicability of SMFS on membrane proteins is its low efficiency in data acquisition. We have developed a semi-automated high-throughput SMFS (HT-SMFS) procedure for efficient data acquisition. In addition, we present a coarse filter to efficiently extract protein unfolding events from large data sets. The HT-SMFS procedure and the coarse filter were validated using the proton pump bacteriorhodopsin (BR) from Halobacterium salinarum and the L-arginine/agmatine antiporter AdiC from the bacterium Escherichia coli. To screen for molecular interactions between AdiC and its substrates, we recorded data sets in the absence and in the presence of L-arginine, D-arginine, and agmatine. Altogether ∼400 000 force-distance curves were recorded. Application of coarse filtering to this wealth of data yielded six data sets with ∼200 (AdiC) and ∼400 (BR) force-distance spectra in each. Importantly, the raw data for most of these data sets were acquired in one to two days, opening new perspectives for HT-SMFS applications.

Molecular energy dissipation in nanoscale networks of Dentin Matrix Protein 1 is strongly dependent on ion valence

The fracture resistance of biomineralized tissues such as bone, dentin, and abalone is greatly enhanced through the nanoscale interactions of stiff inorganic mineral components with soft organic adhesive components. A proper understanding of the interactions that occur within the organic component, and between the organic and inorganic components, is therefore critical for a complete understanding of the mechanics of these tissues. In this paper, we use Atomic Force Microscope (AFM) force spectroscopy and dynamic force spectroscopy to explore the effect of ionic interactions within a nanoscale system consisting of networks of Dentin Matrix Protein 1 (DMP1) (a component of both bone and dentin organic matrix), a mica surface, and an AFM tip. We find that DMP1 is capable of dissipating large amounts of energy through an ion-mediated mechanism, and that the effectiveness increases with increasing ion valence.

Role of extracellular glutamic acids in the stability and energy landscape of bacteriorhodopsin

Bacteriorhodopsin (BR), a specialized nanomachine, converts light energy into a proton gradient to power Halobacterium salinarum. In this work, we analyze the mechanical stability of a BR triple mutant in which three key extracellular residues, Glu(9), Glu(194), and Glu(204), were mutated simultaneously to Gln. These three Glu residues are involved in a network of hydrogen bonds, in cation binding, and form part of the proton release pathway of BR. Changes in these features and the robust photocycle dynamics of wild-type (WT) BR are apparent when the three extracellular Glu residues are mutated to Gln. It is speculated that such functional changes of proteins go hand in hand with changes in their mechanical properties. Here, we apply single-molecule dynamic force spectroscopy to investigate how the Glu to Gln mutations change interactions, reaction pathways, and the energy barriers of the structural regions of WT BR. The altered heights and positions of individual energy barriers unravel the changes in the mechanical and the unfolding kinetic properties of the secondary structures of WT BR. These changes in the mechanical unfolding energy landscape cause the proton pump to choose unfolding pathways differently. We suggest that, in a similar manner, the changed mechanical properties of mutated BR alter the functional energy landscape favoring different reaction pathways in the light-induced proton pumping mechanism.

From valleys to ridges: exploring the dynamic energy landscape of single membrane proteins

Membrane proteins are involved in essential biological processes such as energy conversion, signal transduction, solute transport and secretion. All biological processes, also those involving membrane proteins, are steered by molecular interactions. Molecular interactions guide the folding and stability of membrane proteins, determine their assembly, switch their functional states or mediate signal transduction. The sequential steps of molecular interactions driving these processes can be described by dynamic energy landscapes. The conceptual energy landscape allows to follow the complex reaction pathways of membrane proteins while its modifications describe why and how pathways are changed. Single-molecule force spectroscopy (SMFS) detects, quantifies and locates interactions within and between membrane proteins. SMFS helps to determine how these interactions change with temperature, point mutations, oligomerization and the functional states of membrane proteins. Applied in different modes, SMFS explores the co-existence and population of reaction pathways in the energy landscape of the protein and thus reveals detailed insights into local mechanisms, determining its structural and functional relationships. Here we review how SMFS extracts the defining parameters of an energy landscape such as the barrier position, reaction kinetics and roughness with high precision.

Vertebrate membrane proteins: structure, function, and insights from biophysical approaches

Membrane proteins are key targets for pharmacological intervention because they are vital for cellular function. Here, we analyze recent progress made in the understanding of the structure and function of membrane proteins with a focus on rhodopsin and development of atomic force microscopy techniques to study biological membranes. Membrane proteins are compartmentalized to carry out extra- and intracellular processes. Biological membranes are densely populated with membrane proteins that occupy approximately 50% of their volume. In most cases membranes contain lipid rafts, protein patches, or paracrystalline formations that lack the higher-order symmetry that would allow them to be characterized by diffraction methods. Despite many technical difficulties, several crystal structures of membrane proteins that illustrate their internal structural organization have been determined. Moreover, high-resolution atomic force microscopy, near-field scanning optical microscopy, and other lower resolution techniques have been used to investigate these structures. Single-molecule force spectroscopy tracks interactions that stabilize membrane proteins and those that switch their functional state; this spectroscopy can be applied to locate a ligand-binding site. Recent development of this technique also reveals the energy landscape of a membrane protein, defining its folding, reaction pathways, and kinetics. Future development and application of novel approaches during the coming years should provide even greater insights to the understanding of biological membrane organization and function.

Mechanical properties of bovine rhodopsin and bacteriorhodopsin: possible roles in folding and function

Molecular interactions and mechanical properties that contribute to the stability and function of proteins are complex and of fundamental importance. In this study, we used single-molecule dynamic force spectroscopy (DFS) to explore the interactions and the unfolding energy landscape of bovine rhodopsin and bacteriorhodopsin. An analysis of the experimental data enabled the extraction of parameters that provided insights into the kinetic stability and mechanical properties of these membrane proteins. Individual structural segments of rhodopsin and bacteriorhodopsin have different properties. A core of rigid structural segments was observed in rhodopsin but not in bacteriorhodopsin. This core may reflect differences in mechanisms of protein folding between the two membrane proteins. The different structural rigidity of the two proteins may also reflect their adaptation to differing functions.

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.

Oriented purple membrane monolayers covalently attached to gold by multiple thiole linkages analyzed by single molecule force spectroscopy

Highly oriented monolayers of bacteriorhodopsin (BR) in purple membrane (PM) form are obtained by the reaction of BR-Q3C, where a cysteine was introduced into the N-terminal region, with a gold surface. Single molecule force spectroscopy was used to show that about 50% of the BRs are covalently bound to the surface. The linkage between the cysteine and the gold causes an additional characteristic peak in the force-distance curves to appear. Because several thousand cysteine-gold bonds exist between each PM patch and the surface, the PM is irreversibly bound. Such oriented PM monolayers may serve as an interface between metal surfaces and biomaterials, which may be linked to the PM surface chemically. Photoelectric applications of BR will benefit from the high degree of orientation obtained by this method.

Single-molecule force spectroscopy of mycobacterial adhesin-adhesin interactions

The heparin-binding hemagglutinin (HBHA) is one of the few virulence factors identified for Mycobacterium tuberculosis. It is a surface-associated adhesin that expresses a number of different activities, including mycobacterial adhesion to nonphagocytic cells and microbial aggregation. Previous evidence indicated that HBHA is likely to form homodimers or homopolymers via a predicted coiled-coil region located within the N-terminal portion of the molecule. Here, we used single-molecule atomic-force microscopy to measure individual homophilic HBHA-HBHA interaction forces. Force curves recorded between tips and supports derivatized with HBHA proteins exposing their N-terminal domains showed a bimodal distribution of binding forces reflecting the formation of dimers or multimers. Moreover, the binding peaks showed elongation forces that were consistent with the unfolding of alpha-helical coiled-coil structures. By contrast, force curves obtained for proteins exposing their lysine-rich C-terminal domains showed a broader distribution of binding events, suggesting that they originate primarily from intermolecular electrostatic bridges between cationic and anionic residues rather than from specific coiled-coil interactions. Notably, similar homophilic HBHA-HBHA interactions were demonstrated on live mycobacteria producing HBHA, while they were not observed on an HBHA-deficient mutant. Together with the fact that HBHA mediates bacterial aggregation, these observations suggest that the single homophilic HBHA interactions measured here reflect the formation of multimers that may promote mycobacterial aggregation.

Electrostatic and steric interactions determine bacteriorhodopsin single-molecule biomechanics

Bacteriorhodopsin (bR) is a haloarchaeal membrane protein that converts the energy of single photons into large structural changes to directionally pump protons across purple membrane. This is achieved by a complex combination of local dynamic interactions controlling bR biomechanics at the submolecular level, producing efficient amplification of the retinal photoisomerization. Using single molecule force spectroscopy at different salt concentrations, we show that tryptophan (Trp) residues use steric specific interactions to create a rigid scaffold in bR extracellular region and are responsible for the main unfolding barriers. This scaffold, which encloses the retinal, controls bR local mechanical properties and anchors the protein into the membrane. Furthermore, the stable Trp-based network allows ion binding to two specific sites on the extracellular loops (BC and FG), which are involved in proton release and lateral transport. In contrast, the cytoplasmic side of bR is mainly governed by relatively weak nonspecific electrostatic interactions that provide the flexibility necessary for large cytoplasmic structural rearrangements during the photocycle. The presence of an extracellular Trp-based network tightly enclosing the retinal seems common to most haloarchaeal rhodopsins, and could be relevant to their exceptional efficiency.

Single molecule techniques for the study of membrane proteins

Single molecule techniques promise novel information about the properties and behavior of individual particles, thus enabling access to molecular heterogeneities in biological systems. Their recent developments to accommodate membrane studies have significantly deepened the understanding of membrane proteins. In this short review, we will describe the basics of the three most common single-molecule techniques used on membrane proteins: fluorescence correlation spectroscopy, single particle tracking, and atomic force microscopy. We will discuss the most relevant findings made during the recent years and their contribution to the membrane protein field.

Exploring the molecular forces within and between CbsA S-layer proteins using single molecule force spectroscopy

We used single molecule atomic force microscopy (AFM) to gain insight into the molecular forces driving the folding and assembly of the S-layer protein CbsA. Force curves recorded between tips and supports modified with CbsA proteins showed sawtooth patterns with multiple force peaks of 58+/-26pN that we attribute to the unfolding of alpha-helices, in agreement with earlier secondary structure predictions. The average unfolding force increased with the pulling speed but was independent on the interaction time. Force curves obtained for CbsA peptides truncated in their C-terminal region showed similar periodic features, except that fewer force peaks were seen. Furthermore, the average unfolding force was 83+/-45pN, suggesting the domains were more stable. By contrast, cationic peptides truncated in their N-terminal region showed single force peaks of 366+/-149pN, presumably reflecting intermolecular electrostatic bridges rather than unfolding events. Interestingly, these large intermolecular forces increased not only with pulling speed but also with interaction time. We expect that the intra- and intermolecular forces measured here may play a significant role in controlling the stability and assembly of the CbsA protein.

Free energy of membrane protein unfolding derived from single-molecule force measurements

Mechanical single-molecule techniques offer exciting possibilities to investigate protein folding and stability in native environments at submolecular resolution. By applying a free-energy reconstruction procedure developed by Hummer and Szabo, which is based on a statistical theorem introduced by Jarzynski, we determined the unfolding free energy of the membrane proteins bacteriorhodopsin (BR), halorhodopsin, and the sodium-proton antiporter NhaA. The calculated energies ranged from 290.5 kcal/mol for BR to 485.5 kcal/mol for NhaA. For the remarkably stable BR, the equilibrium unfolding free energy was independent of pulling rate and temperature ranging between 18 and 42 degrees C. Our experiments also revealed heterogeneous energetic properties in individual transmembrane helices. In halorhodopsin, the stabilization of a short helical segment yielded a characteristic signature in the energy profile. In NhaA, a pronounced peak was observed at a functionally important site in the protein. Since a large variety of single- and multispan membrane proteins can be tackled in mechanical unfolding experiments, our approach provides a basis for systematically elucidating energetic properties of membrane proteins with the resolution of individual secondary-structure elements.

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.

Membrane proteins shape up: understanding in vitro folding

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

Sequential unfolding of individual helices of bacterioopsin observed in molecular dynamics simulations of extraction from the purple membrane

Multiple molecular dynamics simulations of bacterioopsin pulling from its C-terminus show that its alpha-helices unfold individually. In the first metastable state observed in the simulations, helix G is unfolded at its C-terminal segment while the rest of helix G (residues 200-216) is folded and opposes resistance because of a salt-bridge network consisting of Asp-212 and Lys-216 on helix G and Arg-82 and Asp-85 on helix C. Helix G unfolds inside the bundle because the external force is applied to its C-terminal end in a direction perpendicular to the surface of the membrane. Inversely, helix F has to flip by 180 degrees to exit from the membrane because the applied force and the helical N-C axis point in opposite directions. At the highest peak of the force, which cannot be interpreted in single-molecule force spectroscopy experiments, helix F has a pronounced kink at Pro-186. Mutation of Pro-186 and/or the charged side chains mentioned above, which are involved in very favorable electrostatic interactions in the low-dielectric region of the membrane, are expected to reduce the highest peak of the force. Helices E and D unfold in a similar way to helices G and F, respectively. Hence, the force-distance profile and sequence of events during forced unfolding of bacterioopsin are influenced by the up-and-down topology of the seven-helix bundle. The sequential extraction of individual helices from the membrane suggests that the spontaneous (un)folding of bacterioopsin proceeds through metastable bundles of fewer than seven helices. The metastable states observed in the simulations provide atomic level evidence that corroborates the interpretation of very recent force spectroscopy experiments of bacteriorhodopsin refolding.

Transmembrane helix-helix association: relative stabilities at low pH

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

Pulling single bacteriorhodopsin out of a membrane: Comparison of simulation and experiment

Mechanical unfolding of single bacteriorhodopsins from a membrane bilayer is studied using molecular dynamics simulations. The initial conformation of the lipid membrane is determined through all-atom simulations and then its coarse-grained representation is used in the studies of stretching. A Go-like model with a realistic contact map and with Lennard-Jones contact interactions is applied to model the protein-membrane system. The model qualitatively reproduces the experimentally observed differences between force-extension patterns obtained on bacteriorhodopsin at different temperatures and predicts a lack of symmetry in the choice of the terminus to pull by. It also illustrates the decisive role of the interactions of the protein with the membrane in determining the force pattern and thus the stability of transmembrane proteins.

Energetics of membrane protein folding and stability

The critical role of membrane proteins in a myriad of biological and physiological functions has spawned numerous investigations over the past several decades with the long-term goal of identifying the molecular origins and energetic forces that stabilize these proteins within the membrane. Parallel structural and thermodynamics studies on several systems have provided significant insight regarding the driving forces governing folding, assembly, insertion, and translocation of membrane proteins. The present review surveys families of membrane-associated proteins including alpha-helical and beta-barrel structures, viral surface receptors, and pore-forming toxins, citing representative proteins within each of these classes for further scrutiny in terms of structure-function relationships and global conformational stability. This overview presents seminal findings from pioneering studies on the energetics of membrane protein folding and stability to modern techniques that are exploiting the use of molecular genetics and single molecule studies. An overall consensus regarding the molecular origins of membrane protein stability is that a number of intrinsic properties resemble features of soluble proteins, yet there are distinct energetic differences arising from specific intra- and intermolecular interactions within the membrane. The combined efforts from structural, energetics, and dynamics approaches offer unique insights and improve our fundamental understanding of the driving forces dictating membrane protein folding and stability.

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

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

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.

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.

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

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

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.