Direct observation of multiple misfolding pathways in a single prion protein molecule

Probing the stability and interdomain interactions in the ABC transporter OpuA using single-molecule optical tweezers

Transmembrane transporter proteins are essential for maintaining cellular homeostasis and, as such, are key drug targets. Many transmembrane transporter proteins are known to undergo large structural rearrangements during their functional cycles. Despite the wealth of detailed structural and functional data available for these systems, our understanding of their dynamics and, consequently, how they function is generally limited. We introduce an innovative approach that enables us to directly measure the dynamics and stability of interdomain interactions of transmembrane proteins using optical tweezers. Focusing on the osmoregulatory ATP-binding cassette transporter OpuA from Lactococcus lactis, we examine the mechanical properties and potential interactions of its substrate-binding domains. Our measurements are performed in lipid nanodiscs, providing a native-mimicking environment for the transmembrane protein. The technique provides high spatial and temporal resolution and allows us to study the functionally relevant motions and interdomain interactions of individual transmembrane transporter proteins in real time in a lipid bilayer.

Mutations of evolutionarily conserved aromatic residues suggest that misfolding of the mouse prion protein may commence in multiple ways

The misfolding of the mammalian prion protein from its α-helix rich cellular isoform to its β-sheet rich infectious isoform is associated with several neurodegenerative diseases. The determination of the structural mechanism by which misfolding commences, still remains an unsolved problem. In the current study, native-state hydrogen exchange coupled with mass spectrometry has revealed that the N state of the mouse prion protein (moPrP) at pH 4 is in dynamic equilibrium with multiple partially unfolded forms (PUFs) capable of initiating misfolding. Mutation of three evolutionarily conserved aromatic residues, Tyr168, Phe174, and Tyr217 present at the interface of the β2-α2 loop and the C-terminal end of α3 in the structured C-terminal domain of moPrP significantly destabilize the native state (N) of the protein. They also reduce the free energy differences between the N state and two PUFs identified as PUF1 and PUF2**. It is shown that PUF2** in which the β2-α2 loop and the C-terminal end of α3 are disordered, has the same stability as the previously identified PUF2*, but to have a very different structure. Misfolding can commence from both PUF1 and PUF2**, as it can from PUF2*. Hence, misfolding can commence and proceed in multiple ways from structurally distinct precursor conformations. The increased extents to which PUF1 and PUF2** are populated at equilibrium in the case of the mutant variants, greatly accelerate their misfolding. The results suggest that the three aromatic residues may have been evolutionarily selected to impede the misfolding of moPrP.

Distinguishing features of fold-switching proteins

Though many folded proteins assume one stable structure that performs one function, a small-but-increasing number remodel their secondary and tertiary structures and change their functions in response to cellular stimuli. These fold-switching proteins regulate biological processes and are associated with autoimmune dysfunction, severe acute respiratory syndrome coronavirus-2 infection, and more. Despite their biological importance, it is difficult to computationally predict fold switching. With the aim of advancing computational prediction and experimental characterization of fold switchers, this review discusses several features that distinguish fold-switching proteins from their single-fold and intrinsically disordered counterparts. First, the isolated structures of fold switchers are less stable and more heterogeneous than single folders but more stable and less heterogeneous than intrinsically disordered proteins (IDPs). Second, the sequences of single fold, fold switching, and intrinsically disordered proteins can evolve at distinct rates. Third, proteins from these three classes are best predicted using different computational techniques. Finally, late-breaking results suggest that single folders, fold switchers, and IDPs have distinct patterns of residue-residue coevolution. The review closes by discussing high-throughput and medium-throughput experimental approaches that might be used to identify new fold-switching proteins.

Functional, pathogenic, and pharmacological roles of protein folding intermediates

Protein expression and function in eukaryotic cells are tightly harmonized processes modulated by the combination of different layers of regulation, including transcription, processing, stability, and translation of messenger RNA, as well as assembly, maturation, sorting, recycling, and degradation of polypeptides. Integrating all these pathways and the protein quality control machinery, deputed to avoid the production and accumulation of aberrantly folded proteins, determines protein homeostasis. Over the last decade, the combined development of accurate time-resolved experimental techniques and efficient computer simulations has opened the possibility of investigating biological mechanisms at atomic resolution with physics-based models. A meaningful example is the reconstruction of protein folding pathways at atomic resolution, which has enabled the characterization of the folding kinetics of biologically relevant globular proteins consisting of a few hundred amino acids. Combining these innovative computational technologies with rigorous experimental approaches reveals the existence of non-native metastable states transiently appearing along the folding process of such proteins. Here, we review the primary evidence indicating that these protein folding intermediates could play roles in disparate biological processes, from the posttranslational regulation of protein expression to disease-relevant protein misfolding mechanisms. Finally, we discuss how the information encoded into protein folding pathways could be exploited to design an entirely new generation of pharmacological agents capable of promoting the selective degradation of protein targets.

Enhanced statistical sampling reveals microscopic complexity in the talin mechanosensor folding energy landscape

Statistical mechanics can describe the major conformational ensembles determining the equilibrium free-energy landscape of a folding protein. The challenge is to capture the full repertoire of low-occurrence conformations separated by high kinetic barriers that define complex landscapes. Computationally, enhanced sampling methods accelerate the exploration of molecular rare events. However, accessing the entire protein's conformational space in equilibrium experiments requires technological developments to enable extended observation times. We developed single-molecule magnetic tweezers to capture over a million individual transitions as a single talin protein unfolds and refolds under force in equilibrium. When observed at classically-probed timescales, talin folds in an apparently uncomplicated two-state manner. As the sampling time extends from minutes to days, the underlying energy landscape exhibits gradually larger signatures of complexity, involving a finite number of well-defined rare conformations. A fluctuation analysis allows us to propose plausible structures of each low-probability conformational state. The physiological relevance of each distinct conformation can be connected to the binding of the cytoskeletal protein vinculin, suggesting an extra layer of complexity in talin-mediated mechanotransduction. More generally, our experiments directly test the fundamental notion that equilibrium dynamics depend on the observation timescale.

Evolutionarily Conserved Proline Residues Impede the Misfolding of the Mouse Prion Protein by Destabilizing an Aggregation-competent Partially Unfolded Form

The misfolding of the prion protein has been linked to several neurodegenerative diseases. Despite extensive studies, the mechanism of the misfolding process remains poorly understood. The present study structurally delineates the role of the conserved proline residues present in the structured C-terminal domain of the mouse prion protein (moPrP) in the misfolding process. It is shown that mutation of these Pro residues to Ala leads to destabilization of the native (N) state, and also to rapid misfolding. Using hydrogen-deuterium exchange (HDX) studies coupled with mass spectrometry (MS), it has been shown that the N state of moPrP is in rapid equilibrium with a partially unfolded form (PUF2*) at pH 4. It has been shown that the Pro to Ala mutations make PUF2* energetically more accessible from the N state by stabilizing it relative to the unfolded (U) state. The apparent rate constant of misfolding is found to be linearly proportional to the extent to which PUF2* is populated in equilibrium with the N state, strongly indicating that misfolding commences from PUF2*. It has also been shown that the Pro residues restrict the boundary of the structural core of the misfolded oligomers. Overall, this study highlights how the conserved proline residues control misfolding of the prion protein by modulating the stability of the partially unfolded form from which misfolding commences.

Probing the origin of prion protein misfolding via reconstruction of ancestral proteins

Prion diseases are fatal neurodegenerative diseases caused by pathogenic misfolding of the prion protein, PrP. They are transmissible between hosts, and sometimes between different species, as with transmission of bovine spongiform encephalopathy to humans. Although PrP is found in a wide range of vertebrates, prion diseases are seen only in certain mammals, suggesting that infectious misfolding was a recent evolutionary development. To explore when PrP acquired the ability to misfold infectiously, we reconstructed the sequences of ancestral versions of PrP from the last common primate, primate-rodent, artiodactyl, placental, bird, and amniote. Recombinant ancestral PrPs were then tested for their ability to form β-sheet aggregates, either spontaneously or when seeded with infectious prion strains from human, cervid, or rodent species. The ability to aggregate developed after the oldest ancestor (last common amniote), and aggregation capabilities diverged along evolutionary pathways consistent with modern-day susceptibilities. Ancestral bird PrP could not be seeded with modern-day prions, just as modern-day birds are resistant to prion disease. Computational modeling of structures suggested that differences in helix 2 could account for the resistance of ancestral bird PrP to seeding. Interestingly, ancestral primate PrP could be converted by all prion seeds, including both human and cervid prions, raising the possibility that species descended from an ancestral primate have retained the susceptibility to conversion by cervid prions. More generally, the results suggest that susceptibility to prion disease emerged prior to ~100 million years ago, with placental mammals possibly being generally susceptible to disease.

Evidence of Orientation-Dependent Early States of Prion Protein Misfolded Structures from Single Molecule Force Spectroscopy

Prion diseases are neurodegenerative disorders characterized by the presence of oligomers and amyloid fibrils. These are the result of protein aggregation processes of the cellular prion protein (PrP^C) into amyloidal forms denoted as prions or PrP^Sc. We employed atomic force microscopy (AFM) for single molecule pulling (single molecule force spectroscopy, SMFS) experiments on the recombinant truncated murine prion protein (PrP) domain to characterize its conformations and potential initial oligomerization processes. Our AFM-SMFS results point to a complex scenario of structural heterogeneity of PrP at the monomeric and dimer level, like other amyloid proteins involved in similar pathologies. By applying this technique, we revealed that the PrP C-terminal domain unfolds in a two-state process. We used two dimeric constructs with different PrP reciprocal orientations: one construct with two sequential PrP in the N- to C-terminal orientation (N-C dimer) and a second one in the C- to C-terminal orientation (C-C dimer). The analysis revealed that the different behavior in terms of unfolding force, whereby the dimer placed C-C dimer unfolds at a higher force compared to the N-C orientation. We propose that the C-C dimer orientation may represent a building block of amyloid fibril formation.

Tethering Complex Proteins and Protein Complexes for Optical Tweezers Experiments

Tethering proteins to force probes, typically micrometer-sized beads, is a prerequisite for dissecting their properties with optical tweezers. DNA handles serve as spacers between the tethered protein of interest and the bead surface. Attachment sites of the DNA handles to both the surface of beads and to the protein of interest must be mechanically stable for optical tweezers experiments. The most prominent method for attaching DNA handles to proteins utilizes thiol chemistry, linking modified DNA to engineered cysteines in the target protein. This method, although experimentally straightforward, is impractical for the large number of proteins that endogenously contain multiple or essential cysteines at undesired positions. Here, we describe two alternative approaches that take advantage of genetically encoded tag sequences in the target protein. The first method uses the enzymes Sfp and BirA, and the second uses the more recently described SpyTag-SpyCatcher system. We outline the process of generating the DNA handles themselves, as well as how to make the DNA-protein chimeras for carrying out optical tweezers experiments. These methods have robustly worked for several diverse and complex proteins, including ones that are difficult to produce or purify, and for protein-containing complexes such as the ribosome. They will be useful in cases where chemistry-based approaches are impractical or not feasible.

Effects of Protein Crowders and Charge on the Folding of Superoxide Dismutase 1 Variants: A Computational Study

The neurodegenerative disease amyotrophic lateral sclerosis (ALS) is associated with the misfolding and aggregation of the metalloenzyme protein superoxide dismutase 1 (SOD1) via mutations that destabilize the monomer-dimer interface. In a cellular environment, crowding and electrostatic screening play essential roles in the folding and aggregation of the SOD1 monomers. Despite numerous studies on the effects of mutations on SOD1 folding, a clear understanding of the interplay between crowding, folding, and aggregation in vivo remains lacking. Using a structure-based minimal model for molecular dynamics simulations, we investigate the role of self-crowding and charge on the folding stability of SOD1 and the G41D mutant where experimentalists were intrigued by an alteration of the folding mechanism by a single point mutation from glycine to charged aspartic acid. We show that unfolded SOD1 configurations are significantly affected by charge and crowding, a finding that would be extremely costly to achieve with all-atom simulations, while the native state is not significantly altered. The mutation at residue 41 alters the interactions between proteins in the unfolded states instead of those within a protein. This paper suggests electrostatics may play an important role in the folding pathway of SOD1 and modifying the charge via mutation and ion concentration may change the dominant interactions between proteins, with potential impacts for aggregation of the mutants. This work provides a plausible reason for the alteration of the unfolded states to address why the mutant G41D causes the changes to the folding mechanism of SOD1 that have intrigued experimentalists.

Heterogeneity in Protein Folding and Unfolding Reactions

Proteins have dynamic structures that undergo chain motions on time scales spanning from picoseconds to seconds. Resolving the resultant conformational heterogeneity is essential for gaining accurate insight into fundamental mechanistic aspects of the protein folding reaction. The use of high-resolution structural probes, sensitive to population distributions, has begun to enable the resolution of site-specific conformational heterogeneity at different stages of the folding reaction. Different states populated during protein folding, including the unfolded state, collapsed intermediate states, and even the native state, are found to possess significant conformational heterogeneity. Heterogeneity in protein folding and unfolding reactions originates from the reduced cooperativity of various kinds of physicochemical interactions between various structural elements of a protein, and between a protein and solvent. Heterogeneity may arise because of functional or evolutionary constraints. Conformational substates within the unfolded state and the collapsed intermediates that exchange at rates slower than the subsequent folding steps give rise to heterogeneity on the protein folding pathways. Multiple folding pathways are likely to represent distinct sequences of structure formation. Insight into the nature of the energy barriers separating different conformational states populated during (un)folding can also be obtained by resolving heterogeneity.

High-Resolution Optical Tweezers Combined with Multicolor Single-Molecule Microscopy

We present an instrument that combines high-resolution optical tweezers and multicolor confocal fluorescence spectroscopy. Biological macromolecules exhibit complex conformation and stoichiometry changes in coordination with their motion and activity. To further our understanding of the complex machinery of life, we need methods that can simultaneously probe more than one degree of freedom of single molecules and complexes. Fluorescence optical tweezers, or "fleezers," combine the capabilities of optical tweezers and single-molecule fluorescence microscopy into a single instrument. Here we present the latest generation of a high-resolution fleezers instrument integrated with multicolor fluorescence spectroscopy. The tweezers portion of the instrument can manipulate biological macromolecules with pN scale forces while measuring subnanometer distances. Simultaneous with tweezers measurements, the multicolor fluorescence capability allows the direct observation of multiple molecules or multiple degrees of freedom which allows, for example, the observation of multiple proteins simultaneously within a complex. The instrument incorporates three fluorescence excitation lasers, all sourced from a single-mode optical fiber allowing a reliable alignment scheme, that allows, for example, three independent fluorescent probes or fluorescence resonance energy transfer (FRET) measurements and also increases flexibility in the choice of fluorescent probes. To avoid photobleaching and improve tweezers stability, the instrument implements a timesharing (using a single trap laser to produce a pair of traps via rapid switching between two locations) and interlacing (turning the trapping beam off when the fluorescence excitation beams are on and vice versa) scheme using acousto-optic modulators (AOM) to rapidly and precisely modulate lasers. Our latest "random phase" trap AOM control method obliterates previous residual trap positioning and bead position measurement errors. Here we present the general design principles and detailed construction and testing protocols for the instrument.

Tweezepy: A Python package for calibrating forces in single-molecule video-tracking experiments

Single-molecule force spectroscopy (SMFS) instruments (e.g., magnetic and optical tweezers) often use video tracking to measure the three-dimensional position of micron-scale beads under an applied force. The force in these experiments is calibrated by comparing the bead trajectory to a thermal motion-based model with the drag coefficient, γ, and trap spring constant, κ, as parameters. Estimating accurate parameters is complicated by systematic biases from spectral distortions, the camera exposure time, parasitic noise, and least-squares fitting methods. However, while robust calibration methods exist that correct for these biases, they are not always used because they can be complex to implement computationally. To address this barrier, we present Tweezepy: a Python package for calibrating forces in SMFS video-tracking experiments. Tweezepy uses maximum likelihood estimation (MLE) to estimate parameters and their uncertainties from a single bead trajectory via the power spectral density (PSD) and Allan variance (AV). It is well-documented, fast, easy to use, and accounts for most common sources of biases in SMFS video-tracking experiments. Here, we provide a comprehensive overview of Tweezepy's calibration scheme, including a review of the theory underlying thermal motion-based parameter estimates, a discussion of the PSD, AV, and MLE, and an explanation of their implementation.

Identification and correction of miscalibration artifacts based on force noise for optical tweezers experiments

Single-molecule force spectroscopy using optical tweezers continues to provide detailed insights into the behavior of nanoscale systems. Obtaining precise measurements of their mechanical properties is highly dependent on accurate instrument calibration. Therefore, instrumental drift or inaccurate calibration may prevent reaching an accuracy at the theoretical limit and may lead to incorrect conclusions. Commonly encountered sources of error include inaccuracies in the detector sensitivity and trap stiffness and neglecting the non-harmonicity of an optical trap at higher forces. Here, we first quantify the impact of these artifacts on force-extension data and find that a small deviation of the calibration parameters can already have a significant downstream effect. We then develop a method to identify and remove said artifacts based on differences in the theoretical and measured noise of bead fluctuations. By applying our procedure to both simulated and experimental data, we can show how effects due to miscalibration and trap non-linearities can be successfully removed. Most importantly, this correction can be performed post-measurement and could be adapted for data acquired using any force spectroscopy technique.

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.

Slow Escape from a Helical Misfolded State of the Pore-Forming Toxin Cytolysin A

The pore-forming toxin cytolysin A (ClyA) is expressed as a large α-helical monomer that, upon interaction with membranes, undergoes a major conformational rearrangement into the protomer conformation, which then assembles into a cytolytic pore. Here, we investigate the folding kinetics of the ClyA monomer with single-molecule Förster resonance energy transfer spectroscopy in combination with microfluidic mixing, stopped-flow circular dichroism experiments, and molecular simulations. The complex folding process occurs over a broad range of time scales, from hundreds of nanoseconds to minutes. The very slow formation of the native state occurs from a rapidly formed and highly collapsed intermediate with large helical content and nonnative topology. Molecular dynamics simulations suggest pronounced non-native interactions as the origin of the slow escape from this deep trap in the free-energy surface, and a variational enhanced path-sampling approach enables a glimpse of the folding process that is supported by the experimental data.

Unfolded and intermediate states of PrP play a key role in the mechanism of action of an antiprion chaperone

Prion and prion-like diseases involve the propagation of misfolded protein conformers. Small-molecule pharmacological chaperones can inhibit propagated misfolding, but how they interact with disease-related proteins to prevent misfolding is often unclear. We investigated how pentosan polysulfate (PPS), a polyanion with antiprion activity in vitro and in vivo, interacts with mammalian prion protein (PrP) to alter its folding. Calorimetry showed that PPS binds two sites on natively folded PrP, but one PPS molecule can bind multiple PrP molecules. Force spectroscopy measurements of single PrP molecules showed PPS stabilizes not only the native fold of PrP but also many different partially folded intermediates that are not observed in the absence of PPS. PPS also bound tightly to unfolded segments of PrP, delaying refolding. These observations imply that PPS can act through multiple possible modes, inhibiting misfolding not only by stabilizing the native fold or sequestering natively folded PrP into aggregates, as proposed previously, but also by binding to partially or fully unfolded states that play key roles in mediating misfolding. These results underline the likely importance of unfolded states as critical intermediates on the prion conversion pathway.

Simple way to correct the drift in surface-coupled optical tweezers using the laser reflection pattern

The surface-coupled optical tweezers are widely used to resolve small units of motion in biology. However, such motions could readily be interfered by the drift between the trap and surface. We present a simple and low-cost method to correct the drift both actively and passively based on video tracking the distance between the laser reflection pattern and the reference bead. As a result, we achieved sub-nanometer resolution and stability for the stuck bead over a broad range of averaging time (0.002-100 s) as demonstrated by the Allan deviation analysis. The sub-nanometer resolution was further manifested with step measurement. Finally, in double-stranded DNA and DNA hairpin stretching experiments, an extension resolution of 1-2 nm with the stability over 120 s has been demonstrated under a constant force. This work thus provides an easy way to bring the benefit of nanometer resolution and long-term stability to the surface-coupled optical tweezers.

From folding to function: complex macromolecular reactions unraveled one-by-one with optical tweezers

Single-molecule manipulation with optical tweezers has uncovered macromolecular behaviour hidden to other experimental techniques. Recent instrumental improvements have made it possible to expand the range of systems accessible to optical tweezers. Beyond focusing on the folding and structural changes of isolated single molecules, optical tweezers studies have evolved into unraveling the basic principles of complex molecular processes such as co-translational folding on the ribosome, kinase activation dynamics, ligand-receptor binding, chaperone-assisted protein folding, and even dynamics of intrinsically disordered proteins (IDPs). In this mini-review, we illustrate the methodological principles of optical tweezers before highlighting recent advances in studying complex protein conformational dynamics - from protein synthesis to physiological function - as well as emerging future issues that are beginning to be addressed with novel approaches.

Optical tweezers in single-molecule biophysics

Optical tweezers have become the method of choice in single-molecule manipulation studies. In this Primer, we first review the physical principles of optical tweezers and the characteristics that make them a powerful tool to investigate single molecules. We then introduce the modifications of the method to extend the measurement of forces and displacements to torques and angles, and to develop optical tweezers with single-molecule fluorescence detection capabilities. We discuss force and torque calibration of these instruments, their various modes of operation and most common experimental geometries. We describe the type of data obtained in each experimental design and their analyses. This description is followed by a survey of applications of these methods to the studies of protein-nucleic acid interactions, protein/RNA folding and molecular motors. We also discuss data reproducibility, the factors that lead to the data variability among different laboratories and the need to develop field standards. We cover the current limitations of the methods and possible ways to optimize instrument operation, data extraction and analysis, before suggesting likely areas of future growth.

Mapping Distinct Sequences of Structure Formation Differentiating Multiple Folding Pathways of a Small Protein

To determine experimentally how the multiple folding pathways of a protein differ, in the order in which the structural parts are assembled, has been a long-standing challenge. To resolve whether structure formation during folding can progress in multiple ways, the complex folding landscape of monellin has been characterized, structurally and temporally, using the multisite time-resolved FRET methodology. After an initial heterogeneous polypeptide chain collapse, structure formation proceeds on parallel pathways. Kinetic analysis of the population evolution data across various protein segments provides a clear structural distinction between the parallel pathways. The analysis leads to a phenomenological model that describes how and when discrete segments acquire structure independently of each other in different subensembles of protein molecules. When averaged over all molecules, structure formation is seen to progress as α-helix formation, followed by core consolidation, then β-sheet formation, and last end-to-end distance compaction. Parts of the protein that are closer in the primary sequence acquire structure before parts separated by longer sequence.

Slow Transition Path Times Reveal a Complex Folding Barrier in a Designed Protein

De-novo designed proteins have received wide interest as potential platforms for nano-engineering and biomedicine. While much work is being done in the design of thermodynamically stable proteins, the folding process of artificially designed proteins is not well-studied. Here we used single-molecule force spectroscopy by optical tweezers to study the folding of ROSS, a de-novo designed 2x2 Rossmann fold. We measured a barrier crossing time in the millisecond range, much slower than what has been reported for other systems. While long transition times can be explained by barrier roughness or slow diffusion, we show that isotropic roughness cannot explain the measured transition path time distribution. Instead, this study shows that the slow barrier crossing of ROSS is caused by the population of three short-lived high-energy intermediates. In addition, we identify incomplete and off-pathway folding events with different barrier crossing dynamics. Our results hint at the presence of a complex transition barrier that may be a common feature of many artificially designed proteins.

Single-molecule studies of amyloid proteins: from biophysical properties to diagnostic perspectives

In neurodegenerative diseases, a wide range of amyloid proteins or peptides such as amyloid-beta and α-synuclein fail to keep native functional conformations, followed by misfolding and self-assembling into a diverse array of aggregates. The aggregates further exert toxicity leading to the dysfunction, degeneration and loss of cells in the affected organs. Due to the disordered structure of the amyloid proteins, endogenous molecules, such as lipids, are prone to interact with amyloid proteins at a low concentration and influence amyloid cytotoxicity. The heterogeneity of amyloid proteinscomplicates the understanding of the amyloid cytotoxicity when relying only on conventional bulk and ensemble techniques. As complementary tools, single-molecule techniques (SMTs) provide novel insights into the different subpopulations of a heterogeneous amyloid mixture as well as the cytotoxicity, in particular as involved in lipid membranes. This review focuses on the recent advances of a series of SMTs, including single-molecule fluorescence imaging, single-molecule force spectroscopy and single-nanopore electrical recording, for the understanding of the amyloid molecular mechanism. The working principles, benefits and limitations of each technique are discussed and compared in amyloid protein related studies.. We also discuss why SMTs show great potential and are worthy of further investigation with feasibility studies as diagnostic tools of neurodegenerative diseases and which limitations are to be addressed.

Effects of Ligand Binding on the Energy Landscape of Acyl-CoA-Binding Protein

Binding of ligands is often crucial for function yet the effects of ligand binding on the mechanical stability and energy landscape of proteins are incompletely understood. Here, we use a combination of single-molecule optical tweezers and MD simulations to investigate the effect of ligand binding on the energy landscape of acyl-coenzyme A (CoA)-binding protein (ACBP). ACBP is a topologically simple and highly conserved four-α-helix bundle protein that acts as an intracellular transporter and buffer for fatty-acyl-CoA and is active in membrane assembly. We have previously described the behavior of ACBP under tension, revealing a highly extended transition state (TS) located almost halfway between the unfolded and native states. Here, we performed force-ramp and force-jump experiments, in combination with advanced statistical analysis, to show that octanoyl-CoA binding increases the activation free energy for the unfolding reaction of ACBP without affecting the position of the transition state along the reaction coordinate. It follows that ligand binding enhances the mechanical resistance and thermodynamic stability of the protein, without changing its mechanical compliance. Steered molecular dynamics simulations allowed us to rationalize the results in terms of key interactions that octanoyl-CoA establishes with the four α-helices of ACBP and showed that the unfolding pathway is marginally affected by the ligand. The results show that ligand-induced mechanical stabilization effects can be complex and may prove useful for the rational design of stabilizing ligands.

Entropy-Based Strategies for Rapid Pre-Processing and Classification of Time Series Data from Single-Molecule Force Experiments

Recent advances in single-molecule science have revealed an astonishing number of details on the microscopic states of molecules, which in turn defined the need for simple, automated processing of numerous time-series data. In particular, large datasets of time series of single protein molecules have been obtained using laser optical tweezers. In this system, each molecular state has a separate time series with a relatively uneven composition from the point of view-point of local descriptive statistics. In the past, uncertain data quality and heterogeneity of molecular states were biased to the human experience. Because the data processing information is not directly transferable to the black-box-framework for an efficient classification, a rapid evaluation of a large number of time series samples simultaneously measured may constitute a serious obstacle. To solve this particular problem, we have implemented a supervised learning method that combines local entropic models with the global Lehmer average. We find that the methodological combination is suitable to perform a fast and simple categorization, which enables rapid pre-processing of the data with minimal optimization and user interventions.

Single-Molecule Studies of Protein Folding with Optical Tweezers

Manipulation of individual molecules with optical tweezers provides a powerful means of interrogating the structure and folding of proteins. Mechanical force is not only a relevant quantity in cellular protein folding and function, but also a convenient parameter for biophysical folding studies. Optical tweezers offer precise control in the force range relevant for protein folding and unfolding, from which single-molecule kinetic and thermodynamic information about these processes can be extracted. In this review, we describe both physical principles and practical aspects of optical tweezers measurements and discuss recent advances in the use of this technique for the study of protein folding. In particular, we describe the characterization of folding energy landscapes at high resolution, studies of structurally complex multidomain proteins, folding in the presence of chaperones, and the ability to investigate real-time cotranslational folding of a polypeptide.

Structural insight into conformational change in prion protein by breakage of electrostatic network around H187 due to its protonation

A conformational change from normal prion protein(PrP^C) to abnormal prion protein(PrP^SC) induces fatal neurodegenerative diseases. Acidic pH is well-known factors involved in the conformational change. Because the protonation of H187 is strongly linked to the change in PrP stability, we examined the charged residues R156, E196, and D202 around H187. Interestingly, there have been reports on pathological mutants, such as H187R, E196A, and D202N. In this study, we focused on how an acidic pH and pathological mutants disrupt this electrostatic network and how this broken network destabilizes PrP structure. To do so, we performed a temperature-based replica-exchange molecular dynamics (T-REMD) simulation using a cumulative 252 μs simulation time. We measured the distance between amino acids comprising four salt bridges (R156–E196/D202 and H187–E196/D202). Our results showed that the spatial configuration of the electrostatic network was significantly altered by an acidic pH and mutations. The structural alteration in the electrostatic network increased the RMSF value around the first helix (H1). Thus, the structural stability of H1, which is anchored to the H2–H3 bundle, was decreased. It induces separation of R156 from the electrostatic network. Analysis of the anchoring energy also shows that two salt-bridges (R156-E196/D202) are critical for PrP stability.

Energetic dependencies dictate folding mechanism in a complex protein

Large proteins composed of multiple domains are abundant in all proteomes, but their folding and structural dynamics remain poorly understood. Using single-molecule force spectroscopy, we have defined how stabilizing interfaces among the domains of elongation factor G (EF-G) shape its folding pathway. Contrary to the expectation that multidomain proteins fold sequentially as they emerge from the ribosome, we find that folding cannot be completed until the full protein has been synthesized. This posttranslational folding mechanism results in a propensity for misfolding. It is dictated by an energetic coupling among domains that enables conformational flexibility crucial for EF-G function. EF-G thus provides an example of how distinct biological ends—robust folding and functionally important flexibility—come into conflict during protein biogenesis.

Large proteins with multiple domains are thought to fold cotranslationally to minimize interdomain misfolding. Once folded, domains interact with each other through the formation of extensive interfaces that are important for protein stability and function. However, multidomain protein folding and the energetics of domain interactions remain poorly understood. In elongation factor G (EF-G), a highly conserved protein composed of 5 domains, the 2 N-terminal domains form a stably structured unit cotranslationally. Using single-molecule optical tweezers, we have defined the steps leading to fully folded EF-G. We find that the central domain III of EF-G is highly dynamic and does not fold upon emerging from the ribosome. Surprisingly, a large interface with the N-terminal domains does not contribute to the stability of domain III. Instead, it requires interactions with its folded C-terminal neighbors to be stably structured. Because of the directionality of protein synthesis, this energetic dependency of domain III on its C-terminal neighbors disrupts cotranslational folding and imposes a posttranslational mechanism on the folding of the C-terminal part of EF-G. As a consequence, unfolded domains accumulate during synthesis, leading to the extensive population of misfolded species that interfere with productive folding. Domain III flexibility enables large-scale conformational transitions that are part of the EF-G functional cycle during ribosome translocation. Our results suggest that energetic tuning of domain stabilities, which is likely crucial for EF-G function, complicates the folding of this large multidomain protein.

Characterizing the inhibition of α-synuclein oligomerization by a pharmacological chaperone that prevents prion formation by the protein PrP

Aggregation of the disordered protein α-synuclein into amyloid fibrils is a central feature of synucleinopathies, neurodegenerative disorders that include Parkinson's disease. Small, pre-fibrillar oligomers of misfolded α-synuclein are thought to be the key toxic entities, and α-synuclein misfolding can propagate in a prion-like way. We explored whether a compound with anti-prion activity that can bind to unfolded parts of the protein PrP, the cyclic tetrapyrrole Fe-TMPyP, was also active against α-synuclein aggregation. Observing the initial stages of aggregation via fluorescence cross-correlation spectroscopy, we found that Fe-TMPyP inhibited small oligomer formation in a dose-dependent manner. Fe-TMPyP also inhibited the formation of mature amyloid fibrils in vitro, as detected by thioflavin T fluorescence. Isothermal titration calorimetry indicated Fe-TMPyP bound to monomeric α-synuclein with a stoichiometry of 2, and two-dimensional heteronuclear single quantum coherence NMR spectra revealed significant interactions between Fe-TMPyP and the C-terminus of the protein. These results suggest commonalities among aggregation mechanisms for α-synuclein and the prion protein may exist that can be exploited as therapeutic targets.

Non-equilibrium dynamics of a nascent polypeptide during translation suppress its misfolding

Protein folding can begin co-translationally. Due to the difference in timescale between folding and synthesis, co-translational folding is thought to occur at equilibrium for fast-folding domains. In this scenario, the folding kinetics of stalled ribosome-bound nascent chains should match the folding of nascent chains in real time. To test if this assumption is true, we compare the folding of a ribosome-bound, multi-domain calcium-binding protein stalled at different points in translation with the nascent chain as is it being synthesized in real-time, via optical tweezers. On stalled ribosomes, a misfolded state forms rapidly (1.5 s). However, during translation, this state is only attained after a long delay (63 s), indicating that, unexpectedly, the growing polypeptide is not equilibrated with its ensemble of accessible conformations. Slow equilibration on the ribosome can delay premature folding until adequate sequence is available and/or allow time for chaperone binding, thus promoting productive folding.

Co-translational protein folding is thought to occur at equilibrium for fast-folding domains. Here authors use optical tweezers to show that the folding kinetics of stalled ribosome-bound nascent chains do not match the folding of nascent chains in real time.

Microsecond sub-domain motions and the folding and misfolding of the mouse prion protein

Protein aggregation appears to originate from partially unfolded conformations that are sampled through stochastic fluctuations of the native protein. It has been a challenge to characterize these fluctuations, under native like conditions. Here, the conformational dynamics of the full-length (23-231) mouse prion protein were studied under native conditions, using photoinduced electron transfer coupled to fluorescence correlation spectroscopy (PET-FCS). The slowest fluctuations could be associated with the folding of the unfolded state to an intermediate state, by the use of microsecond mixing experiments. The two faster fluctuations observed by PET-FCS, could be attributed to fluctuations within the native state ensemble. The addition of salt, which is known to initiate the aggregation of the protein, resulted in an enhancement in the time scale of fluctuations in the core of the protein. The results indicate the importance of native state dynamics in initiating the aggregation of proteins.

Using Single-Molecule Chemo-Mechanical Unfolding to Simultaneously Probe Multiple Structural Parameters in Protein Folding

While single-molecule force spectroscopy has greatly advanced the study of protein folding, there are limitations to what can be learned from studying the effect of force alone. We developed a novel technique, chemo-mechanical unfolding, that combines multiple perturbants—force and chemical denaturant—to more fully characterize the folding process by simultaneously probing multiple structural parameters—the change in end-to-end distance, and solvent accessible surface area. Here, we describe the theoretical background, experimental design, and data analysis for chemo-mechanical unfolding experiments probing protein folding thermodynamics and kinetics. This technique has been applied to characterize parallel protein folding pathways, the protein denatured state, protein folding on the ribosome, and protein folding intermediates.

The Ribosome Cooperates with a Chaperone to Guide Multi-domain Protein Folding

Multi-domain proteins, containing several structural units within a single polypeptide, constitute a large fraction of all proteomes. Co-translational folding is assumed to simplify the conformational search problem for large proteins, but the events leading to correctly folded, functional structures remain poorly characterized. Similarly, how the ribosome and molecular chaperones promote efficient folding remains obscure. Using optical tweezers, we have dissected early folding events of nascent elongation factor G, a multi-domain protein that requires chaperones for folding. The ribosome and the chaperone trigger factor reduce inter-domain misfolding, permitting folding of the N-terminal G-domain. Successful completion of this step is a crucial prerequisite for folding of the next domain. Unexpectedly, co-translational folding does not proceed unidirectionally; emerging unfolded polypeptide can denature an already-folded domain. Trigger factor, but not the ribosome, protects against denaturation. The chaperone thus serves a previously unappreciated function, helping multi-domain proteins overcome inherent challenges during co-translational folding.

Polarization induced control of optical trap potentials in binary liquids

We illustrate control of a polarized laser optical trapping potential landscape through the nonideal mixing of binary liquids. The inherent trapping potential asymmetry (ITPA) present in the trapping region results from the asymmetric intensity distribution in focal volume due to the high numerical aperture objective lens. Experimentally, we show that this ITPA effect can be modified and/or removed by the use of binary liquid mixtures. From our femtosecond optical tweezers experiments, we determine the topograph of the trapping potential base on the fluctuation-dissipation theorem. Additionally, the Brownian motion of the trapped bead is sensitive to the frictional force (FF) of the surroundings that is exerted by clusters of water and alcohol binary mixture through extended hydrogen bonding. Thus, using these two effects, ITPA and FF of the medium, we have shown that one can indeed modify the effective trapping potential landscape. Water-alcohol binary mixtures display a nonlinear dependence on the microrheological properties of the solvent composition as a result of rigid cluster formation. Volumetrically, at about 30% methanol in water binary mixture, the trapping asymmetry is minimal. In this particular binary mixture composition, the hydrophobic part of the methanol molecule is surrounded by ‘cages’ of water molecules. Enhanced H-bonding network of water molecules results in higher viscosity, which contributes to the higher frictional force. Increased viscosity decreases the degree of anisotropy due to hindered dipolar rotation. However, at higher methanol concentrations, the methanol molecules are no longer contained within the water cages and are free to move, which decrease their overall bulk viscosity. Thus, for pure solvents, experimentally measured anisotropy matches quite well with the theoretical prediction, but this fails in case of the binary mixtures due to the increased frictional force exerted by binary mixtures that result from the formation of cage-like structures.

Mechanobiology: protein refolding under force

The application of direct force to a protein enables to probe wide regions of its energy surface through conformational transitions as unfolding, extending, recoiling, collapsing, and refolding. While unfolding under force typically displayed a two-state behavior, refolding under force, from highly extended unfolded states, displayed a more complex behavior. The first recording of protein refolding at a force quench step displayed an initial rapid elastic recoil, followed by a plateau phase at some extension, concluding with a collapse to a final state, at which refolding occurred. These findings stirred a lively discussion, which led to further experimental and theoretical investigation of this behavior. It was demonstrated that the polymeric chain of the unfolded protein is required to fully collapse to a globular conformation for the maturation of native structure. This behavior was modeled using one-dimensional free energy landscape over the end-to-end length reaction coordinate, the collective measured variable. However, at low forces, conformational space is not well captured by such models, and using two-dimensional energy surfaces provides further insight into the dynamics of this process. This work reviews the main concepts of protein refolding under constant force, which is essential for understanding how mechanotransducing proteins operate in vivo.

The Complex Conformational Dynamics of Neuronal Calcium Sensor-1: A Single Molecule Perspective

The human neuronal calcium sensor-1 (NCS-1) is a multispecific two-domain EF-hand protein expressed predominantly in neurons and is a member of the NCS protein family. Structure-function relationships of NCS-1 have been extensively studied showing that conformational dynamics linked to diverse ion-binding is important to its function. NCS-1 transduces Ca²⁺ changes in neurons and is linked to a wide range of neuronal functions such as regulation of neurotransmitter release, voltage-gated Ca²⁺ channels and neuronal outgrowth. Defective NCS-1 can be deleterious to cells and has been linked to serious neuronal disorders like autism. Here, we review recent studies describing at the single molecule level the structural and mechanistic details of the folding and misfolding processes of the non-myristoylated NCS-1. By manipulating one molecule at a time with optical tweezers, the conformational equilibria of the Ca²⁺-bound, Mg²⁺-bound and apo states of NCS-1 were investigated revealing a complex folding mechanism underlain by a rugged and multidimensional energy landscape. The molecular rearrangements that NCS-1 undergoes to transit from one conformation to another and the energetics of these reactions are tightly regulated by the binding of divalent ions (Ca²⁺ and Mg²⁺) to its EF-hands. At pathologically high Ca²⁺ concentrations the protein sometimes follows non-productive misfolding pathways leading to kinetically trapped and potentially harmful misfolded conformations. We discuss the significance of these misfolding events as well as the role of inter-domain interactions in shaping the energy landscape and ultimately the biological function of NCS-1. The conformational equilibria of NCS-1 are also compared to those of calmodulin (CaM) and differences and similarities in the behavior of these proteins are rationalized in terms of structural properties.

Structural mechanisms of oligomer and amyloid fibril formation by the prion protein

Misfolding and aggregation of the prion protein is responsible for multiple neurodegenerative diseases. Works from several laboratories on folding of both the WT and multiple pathogenic mutant variants of the prion protein have identified several structurally dissimilar intermediates, which might be potential precursors to misfolding and aggregation. The misfolded aggregates themselves are morphologically distinct, critically dependent on the solution conditions under which they are prepared, but always β-sheet rich. Despite the lack of an atomic resolution structure of the infectious pathogenic agent in prion diseases, several low resolution models have identified the β-sheet rich core of the aggregates formed in vitro, to lie in the α2-α3 subdomain of the prion protein, albeit with local stabilities that vary with the type of aggregate. This feature article describes recent advances in the investigation of in vitro prion protein aggregation using multiple spectroscopic probes, with particular focus on (1) identifying aggregation-prone conformations of the monomeric protein, (2) conditions which trigger misfolding and oligomerization, (3) the mechanism of misfolding and aggregation, and (4) the structure of the misfolded intermediates and final aggregates.

Size and topology modulate the effects of frustration in protein folding

The presence of conflicting interactions, or frustration, determines how fast biomolecules can explore their configurational landscapes. Recent experiments have provided cases of systems with slow reconfiguration dynamics, perhaps arising from frustration. While it is well known that protein folding speed and mechanism are strongly affected by the protein native structure, it is still unknown how the response to frustration is modulated by the protein topology. We explore the effects of nonnative interactions in the reconfigurational and folding dynamics of proteins with different sizes and topologies. We find that structural correlations related to the folded state size and topology play an important role in determining the folding kinetics of proteins that otherwise have the same amount of nonnative interactions. In particular, we find that the reconfiguration dynamics of α-helical proteins are more susceptible to frustration than β-sheet proteins of the same size. Our results may explain recent experimental findings and suggest that attempts to measure the degree of frustration due to nonnative interactions might be more successful with α-helical proteins.

Communication: Coordinate-dependent diffusivity from single molecule trajectories

Single-molecule observations of biomolecular folding are commonly interpreted using the model of one-dimensional diffusion along a reaction coordinate, with a coordinate-independent diffusion coefficient. Recent analysis, however, suggests that more general models are required to account for single-molecule measurements performed with high temporal resolution. Here, we consider one such generalization: a model where the diffusion coefficient can be an arbitrary function of the reaction coordinate. Assuming Brownian dynamics along this coordinate, we derive an exact expression for the coordinate-dependent diffusivity in terms of the splitting probability within an arbitrarily chosen interval and the mean transition path time between the interval boundaries. This formula can be used to estimate the effective diffusion coefficient along a reaction coordinate directly from single-molecule trajectories.

A peptide-display protein scaffold to facilitate single molecule force studies of aggregation-prone peptides

Protein aggregation is linked with the onset of several neurodegenerative disorders, including Parkinson's disease (PD), which is associated with the aggregation of α‐synuclein (αSyn). The structural mechanistic details of protein aggregation, including the nature of the earliest protein–protein interactions, remain elusive. In this study, we have used single molecule force spectroscopy (SMFS) to probe the first dimerization events of the central aggregation‐prone region of αSyn (residues 71–82) that may initiate aggregation. This region has been shown to be necessary for the aggregation of full length αSyn and is capable of forming amyloid fibrils in isolation. We demonstrate that the interaction of αSyn_71‐82 peptides can be studied using SMFS when inserted into a loop of protein L, a mechanically strong and soluble scaffold protein that acts as a display system for SMFS studies. The corresponding fragment of the homolog protein γ‐synuclein (γSyn), which has a lower aggregation propensity, has also been studied here. The results from SMFS, together with native mass spectrometry and aggregation assays, demonstrate that the dimerization propensity of γSyn_71‐82 is lower than that of αSyn_71‐82, but that a mixed αSyn_71‐82: γSyn_71‐82 dimer forms with a similar propensity to the αSyn_71‐82 homodimer, slowing amyloid formation. This work demonstrates the utility of a novel display method for SMFS studies of aggregation‐prone peptides, which would otherwise be difficult to study.

Protein Tethering for Folding Studies

Optical tweezers allow the detection of unfolding and refolding transitions in individual proteins, and how interacting molecules such as chaperones affect these transitions. Typical methods that tether individual proteins are based on cysteine chemistry, which is less suitable for proteins with essential cysteines. Here we describe a cysteine-independent tethering protocol that can be performed in situ.

Partially native intermediates mediate misfolding of SOD1 in single-molecule folding trajectories

Prion-like misfolding of superoxide dismutase 1 (SOD1) is associated with the disease ALS, but the mechanism of misfolding remains unclear, partly because misfolding is difficult to observe directly. Here we study the most misfolding-prone form of SOD1, reduced un-metallated monomers, using optical tweezers to measure unfolding and refolding of single molecules. We find that the folding is more complex than suspected, resolving numerous previously undetected intermediate states consistent with the formation of individual β-strands in the native structure. We identify a stable core of the protein that unfolds last and refolds first, and directly observe several distinct misfolded states that branch off from the native folding pathways at specific points after the formation of the stable core. Partially folded intermediates thus play a crucial role mediating between native and non-native folding. These results suggest an explanation for SOD1's propensity for prion-like misfolding and point to possible targets for therapeutic intervention.

Force Spectroscopy with 9-μs Resolution and Sub-pN Stability by Tailoring AFM Cantilever Geometry

Atomic force microscopy (AFM)-based single-molecule force spectroscopy (SMFS) is a powerful yet accessible means to characterize the unfolding/refolding dynamics of individual molecules and resolve closely spaced, transiently occupied folding intermediates. On a modern commercial AFM, these applications and others are now limited by the mechanical properties of the cantilever. Specifically, AFM-based SMFS data quality is degraded by a commercial cantilever's limited combination of temporal resolution, force precision, and force stability. Recently, we modified commercial cantilevers with a focused ion beam to optimize their properties for SMFS. Here, we extend this capability by modifying a 40 × 18 μm² cantilever into one terminated with a gold-coated, 4 × 4 μm² reflective region connected to an uncoated 2-μm-wide central shaft. This "Warhammer" geometry achieved 8.5-μs resolution coupled with improved force precision and sub-pN stability over 100 s when measured on a commercial AFM. We highlighted this cantilever's biological utility by first resolving a calmodulin unfolding intermediate previously undetected by AFM and then measuring the stabilization of calmodulin by myosin light chain kinase at dramatically higher unfolding velocities than in previous AFM studies. More generally, enhancing data quality via an improved combination of time resolution, force precision, and force stability will broadly benefit biological applications of AFM.

Comparing the energy landscapes for native folding and aggregation of PrP

Protein sequences are evolved to encode generally one folded structure, out of a nearly infinite array of possible folds. Underlying this code is a funneled free energy landscape that guides folding to the native conformation. Protein misfolding and aggregation are also a manifestation of free-energy landscapes. The detailed mechanisms of these processes are poorly understood, but often involve rare, transient species and a variety of different pathways. The inherent complexity of misfolding has hampered efforts to measure aggregation pathways and the underlying energy landscape, especially using traditional methods where ensemble averaging obscures important rare and transient events. We recently studied the misfolding and aggregation of prion protein by examining 2 monomers tethered in close proximity as a dimer, showing how the steps leading to the formation of a stable aggregated state can be resolved in the single-molecule limit and the underlying energy landscape thereby reconstructed. This approach allows a more quantitative comparison of native folding versus misfolding, including fundamental differences in the dynamics for misfolding. By identifying key steps and interactions leading to misfolding, it should help to identify potential drug targets. Here we describe the importance of characterizing free-energy landscapes for aggregation and the challenges involved in doing so, and we discuss how single-molecule studies can help test proposed structural models for PrP aggregates.

Experimental measurement of binding energy, selectivity, and allostery using fluctuation theorems

Thermodynamic bulk measurements of binding reactions rely on the validity of the law of mass action and the assumption of a dilute solution. Yet, important biological systems such as allosteric ligand-receptor binding, macromolecular crowding, or misfolded molecules may not follow these assumptions and may require a particular reaction model. Here we introduce a fluctuation theorem for ligand binding and an experimental approach using single-molecule force spectroscopy to determine binding energies, selectivity, and allostery of nucleic acids and peptides in a model-independent fashion. A similar approach could be used for proteins. This work extends the use of fluctuation theorems beyond unimolecular folding reactions, bridging the thermodynamics of small systems and the basic laws of chemical equilibrium.

Conformational dynamics of the frameshift stimulatory structure in HIV-1

Programmed ribosomal frameshifting (PRF) in HIV-1 is thought to be stimulated by a hairpin in the mRNA, although a pseudoknot-like triplex has also been proposed. Because the conformational dynamics of the stimulatory structure under tension applied by the ribosomal helicase during translation may play an important role in PRF, we used optical tweezers to apply tension to the HIV stimulatory structure and monitor its unfolding and refolding dynamics. The folding and unfolding kinetics and energy landscape of the hairpin were measured by ramping the force on the hairpin up and down, providing a detailed biophysical characterization. Unexpectedly, whereas unfolding reflected the simple two-state behavior typical of many hairpins, refolding was more complex, displaying significant heterogeneity. Evidence was found for multiple refolding pathways as well as previously unsuspected, partially folded intermediates. Measuring a variant mRNA containing only the sequence required to form the proposed triplex, it behaved largely in the same way. Nonetheless, very rarely, high-force unfolding events characteristic of pseudoknot-like structures were observed. The rare occurrence of the triplex suggests that the hairpin is the functional stimulatory structure. The unusual heterogeneity of the hairpin dynamics under tension suggests a possible functional role in PRF similar to the dynamics of other stimulatory structures.