Quantifying the Native Energetics Stabilizing Bacteriorhodopsin by Single-Molecule Force Spectroscopy

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.

Quantifying a light-induced energetic change in bacteriorhodopsin by force spectroscopy

Ligand-induced conformational changes are critical to the function of many membrane proteins and arise from numerous intramolecular interactions. In the photocycle of the model membrane protein bacteriorhodopsin (bR), absorption of a photon by retinal triggers a conformational cascade that results in pumping a proton across the cell membrane. While decades of spectroscopy and structural studies have probed this photocycle in intricate detail, changes in intramolecular energetics that underlie protein motions have remained elusive to experimental quantification. Here, we measured these energetics on the millisecond time scale using atomic-force-microscopy-based single-molecule force spectroscopy. Precisely, timed light pulses triggered the bR photocycle while we measured the equilibrium unfolding and refolding of the terminal 8-amino-acid region of bR's G-helix. These dynamics changed when the EF-helix pair moved ~9 Å away from this end of the G helix during the "open" portion of bR's photocycle. In ~60% of the data, we observed abrupt light-induced destabilization of 3.4 ± 0.3 kcal/mol, lasting 38 ± 3 ms. The kinetics and pH-dependence of this destabilization were consistent with prior measurements of bR's open phase. The frequency of light-induced destabilization increased with the duration of illumination and was dramatically reduced in the triple mutant (D96G/F171C/F219L) thought to trap bR in its open phase. In the other ~40% of the data, photoexcitation unexpectedly stabilized a longer-lived putative misfolded state. Through this work, we establish a general single-molecule force spectroscopy approach for measuring ligand-induced energetics and lifetimes in membrane proteins.

Structural heterogeneity of the ion and lipid channel TMEM16F

Transmembrane protein 16 F (TMEM16F) is a Ca2+-activated homodimer which functions as an ion channel and a phospholipid scramblase. Despite the availability of several TMEM16F cryogenic electron microscopy (cryo-EM) structures, the mechanism of activation and substrate translocation remains controversial, possibly due to restrictions in the accessible protein conformational space. In this study, we use atomic force microscopy under physiological conditions to reveal a range of structurally and mechanically diverse TMEM16F assemblies, characterized by variable inter-subunit dimerization interfaces and protomer orientations, which have escaped prior cryo-EM studies. Furthermore, we find that Ca2+-induced activation is associated to stepwise changes in the pore region that affect the mechanical properties of transmembrane helices TM3, TM4 and TM6. Our direct observation of membrane remodelling in response to Ca2+ binding along with additional electrophysiological analysis, relate this structural multiplicity of TMEM16F to lipid and ion permeation processes. These results thus demonstrate how conformational heterogeneity of TMEM16F directly contributes to its diverse physiological functions.

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.

Single-Molecule Force Spectroscopy Reveals Cation-π Interactions in Aqueous Media Are Highly Affected by Cation Dehydration

Cation-π interactions underlie many important processes in biology and materials science. However, experimental investigations of cation-π interactions in aqueous media remain challenging. Here, we studied the cation-π binding strength and mechanism by pulling two hydrophobic polymers with distinct cation binding properties, i.e., poly-pentafluorostyrene and polystyrene, in aqueous media using single-molecule force spectroscopy and nuclear magnetic resonance measurement. We found that the interaction strengths linearly depend on the cation concentrations, following the order of Li^{+}<NH_{4}^{+}<Na^{+}<K^{+}. The binding energies are 0.03-0.23  kJ mol^{-1} M^{-1}. This order is distinct from the strength of cation-π interactions in gas phase and may be caused by the different dehydration ability of the cations. Taken together, our method provides a unique perspective to investigate cation-π interactions under physiologically relevant conditions.

How physical forces drive the process of helical membrane protein folding

Protein folding is a fundamental process of life with important implications throughout biology. Indeed, tens of thousands of mutations have been associated with diseases, and most of these mutations are believed to affect protein folding rather than function. Correct folding is also a key element of design. These factors have motivated decades of research on protein folding. Unfortunately, knowledge of membrane protein folding lags that of soluble proteins. This gap is partly caused by the greater technical challenges associated with membrane protein studies, but also because of additional complexities. While soluble proteins fold in a homogenous water environment, membrane proteins fold in a setting that ranges from bulk water to highly charged to apolar. Thus, the forces that drive folding vary in different regions of the protein, and this complexity needs to be incorporated into our understanding of the folding process. Here, we review our understanding of membrane protein folding biophysics. Despite the greater challenge, better model systems and new experimental techniques are starting to unravel the forces and pathways in membrane protein folding.

Transforming de novo protein α3D into a mechanically stable protein by zinc binding

α₃D is a de novo designed three-helix bundle protein. Like most naturally occurring helical proteins, it is mechanically labile with an unfolding force of <15 pN, revealed by atomic force microscopy-based single-molecule force spectroscopy (AFM-SMFS). This protein has been further designed with a tri-cysteine metal-binding site, named α₃DIV, which can bind heavy transition metals. Here, we demonstrate that incorporating such a metal-binding site can transform this mechanically labile protein into a stable one. We show that zinc binds to the tri-cysteine site and increases the unfolding force to ∼160 pN. This force is one order of magnitude higher than that of the apo-protein (<15 pN). Moreover, the unfolding mechanism of Zn-α₃DIV indicates the correct zinc binding with the tri-cysteine site, forming three mechanostable Zn-thiolate bonds. Thus, α₃DIV could be a potential α-helical structure-based building block for synthesizing biomaterials with tunable mechanical properties.

Single-Molecule Force Spectroscopy of Protein Folding

The use of force probes to induce unfolding and refolding of single molecules through the application of mechanical tension, known as single-molecule force spectroscopy (SMFS), has proven to be a powerful tool for studying the dynamics of protein folding. Here we provide an overview of what has been learned about protein folding using SMFS, from small, single-domain proteins to large, multi-domain proteins. We highlight the ability of SMFS to measure the energy landscapes underlying folding, to map complex pathways for native and non-native folding, to probe the mechanisms of chaperones that assist with native folding, to elucidate the effects of the ribosome on co-translational folding, and to monitor the folding of membrane proteins.

Computational Investigation of the pH Dependence of Stability of Melanosome Proteins: Implication for Melanosome formation and Disease

Intravesicular pH plays a crucial role in melanosome maturation and function. Melanosomal pH changes during maturation from very acidic in the early stages to neutral in late stages. Neutral pH is critical for providing optimal conditions for the rate-limiting, pH-sensitive melanin-synthesizing enzyme tyrosinase (TYR). This dramatic change in pH is thought to result from the activity of several proteins that control melanosomal pH. Here, we computationally investigated the pH-dependent stability of several melanosomal membrane proteins and compared them to the pH dependence of the stability of TYR. We confirmed that the pH optimum of TYR is neutral, and we also found that proteins that are negative regulators of melanosomal pH are predicted to function optimally at neutral pH. In contrast, positive pH regulators were predicted to have an acidic pH optimum. We propose a competitive mechanism among positive and negative regulators that results in pH equilibrium. Our findings are consistent with previous work that demonstrated a correlation between the pH optima of stability and activity, and they are consistent with the expected activity of positive and negative regulators of melanosomal pH. Furthermore, our data suggest that disease-causing variants impact the pH dependence of melanosomal proteins; this is particularly prominent for the OCA2 protein. In conclusion, melanosomal pH appears to affect the activity of multiple melanosomal proteins.

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.

Correcting molecular transition rates measured by single-molecule force spectroscopy for limited temporal resolution

Equilibrium free-energy-landscape parameters governing biomolecular folding can be determined from nonequilibrium force-induced unfolding by measuring the rates k for transitioning back and forth between states as a function of force F. However, bias in the observed forward and reverse rates is introduced by limited effective temporal resolution, which includes the mechanical response time of the force probe and any smoothing used to improve the signal-to-noise ratio. Here we use simulations to characterize this bias, which is most prevalent when the ratio of forward and reverse rates is far from unity. We find deviations in k(F) at high rates, due to unobserved transitions from short- to long-lived states, and at low rates, due to the corresponding unobserved transitions from long- to short-lived states. These missing events introduce erroneous curvature in log(k) vs F that leads to incorrect landscape parameter determination. To correct the measured k(F), we derive a pair of model-independent analytical formulas. The first correction accounts for unobserved transitions from short- to long-lived states, but does surprisingly little to correct the erroneous energy-landscape parameters. Only by subsequently applying the second formula, which corrects the corresponding reverse process, do we recover the expected k(F) and energy-landscape quantities. Going forward, these corrections should be applied to transition-rate data whenever the highest measured rate is not at least an order of magnitude slower than the effective temporal resolution.