Single-molecule chemo-mechanical unfolding reveals multiple transition state barriers in a small single-domain protein

Ligand-dependent folding and unfolding dynamics and free energy landscapes of acylphosphatase

Acylphosphatase (AcP) is an enzyme which catalyses the hydrolysis of acylphosphate. The binding with the phosphate ion (Pi) assumes significance in preserving both the stability and enzymatic activity of AcP. While previous studies using single molecule force spectroscopy explored the mechanical properties of AcP, the influence of Pi on its folding and unfolding dynamic behaviors remains unexplored. In this work, using stable magnetic tweezers, we measured and compared the force-dependent folding and unfolding rates of AcP in the Tris buffer and phosphate buffer within a force range from 2 pN to 40 pN. We found that Pi exerts no discernible effect on the folding dynamics but consistently decreases the force-dependent unfolding rate of AcP by a constant ratio across the entire force spectrum. The free energy landscapes of AcP in the absence and presence of Pi are constructed. Our results reveal that Pi selectively binds to the native state of AcP, stabilizing it and suggesting the general properties of specific ligand-receptor interactions.

Pathfinder: Protein folding pathway prediction based on conformational sampling

The study of protein folding mechanism is a challenge in molecular biology, which is of great significance for revealing the movement rules of biological macromolecules, understanding the pathogenic mechanism of folding diseases, and designing protein engineering materials. Based on the hypothesis that the conformational sampling trajectory contain the information of folding pathway, we propose a protein folding pathway prediction algorithm named Pathfinder. Firstly, Pathfinder performs large-scale sampling of the conformational space and clusters the decoys obtained in the sampling. The heterogeneous conformations obtained by clustering are named seed states. Then, a resampling algorithm that is not constrained by the local energy basin is designed to obtain the transition probabilities of seed states. Finally, protein folding pathways are inferred from the maximum transition probabilities of seed states. The proposed Pathfinder is tested on our developed test set (34 proteins). For 11 widely studied proteins, we correctly predicted their folding pathways and specifically analyzed 5 of them. For 13 proteins, we predicted their folding pathways to be further verified by biological experiments. For 6 proteins, we analyzed the reasons for the low prediction accuracy. For the other 4 proteins without biological experiment results, potential folding pathways were predicted to provide new insights into protein folding mechanism. The results reveal that structural analogs may have different folding pathways to express different biological functions, homologous proteins may contain common folding pathways, and α-helices may be more prone to early protein folding than β-strands.

Free-energy landscape of two-state protein acylphosphatase with large contact order revealed by force-dependent folding and unfolding dynamics

Acylphosphatase (AcP) is a small protein with 98 amino acid residues that catalyzes the hydrolysis of carboxyl-phosphate bonds. AcP is a typical two-state protein with slow folding rate due to its relatively large contact order in the native structure. The mechanical properties and unfolding behavior of AcP has been studied by atomic force microscope. Here using stable magnetic tweezers, we measured the force-dependent folding rates within a force range 1-3 pN, and unfolding rates 15-40 pN. The obtained unfolding rates show different force sensitivities at forces below and above ∼27 pN, which determines a free-energy landscape with two energy barriers. Our results indicate that the free-energy landscape of small globule proteins have general Bactrian camel shape, and large contact order of the native state produces a high barrier dominate at low forces.

Molten globule-like transition state of protein barnase measured with calorimetric force spectroscopy

SignificanceUnderstanding the molecular forces driving the unfolded polypeptide chain to self-assemble into a functional native structure remains an open question. However, identifying the states visited during protein folding (e.g., the transition state between the unfolded and native states) is tricky due to their transient nature. Here, we introduce calorimetric force spectroscopy in a temperature jump optical trap to determine the enthalpy, entropy, and heat capacity of the transition state of protein barnase. We find that the transition state has the properties of a dry molten globule, that is, high free energy and low configurational entropy, being structurally similar to the native state. This experimental single-molecule study characterizes the thermodynamic properties of the transition state in funneled energy landscapes.

Observing the base-by-base search for native structure along transition paths during the folding of single nucleic acid hairpins

Biomolecular folding involves searching among myriad possibilities for the native conformation, but the elementary steps expected from theory for this search have never been detected directly. We probed the dynamics of folding at high resolution using optical tweezers, measuring individual trajectories as nucleic acid hairpins passed through the high-energy transition states that dominate kinetics and define folding mechanisms. We observed brief but ubiquitous pauses in the transition states, with a dwell time distribution that matched microscopic theories of folding quantitatively. The sequence dependence suggested that pauses were dominated by microbarriers from nonnative conformations during the search by each nucleotide residue for the native base-pairing conformation. Furthermore, the pauses were position dependent, revealing subtle local variations in energy-landscape roughness and allowing the diffusion coefficient describing the microscopic dynamics within the barrier to be found without reconstructing the shape of the energy landscape. These results show how high-resolution measurements can elucidate key microscopic events during folding to test fundamental theories of folding.

Protein Folding Dynamics as Diffusion on a Free Energy Surface: Rate Equation Terms, Transition Paths, and Analysis of Single-Molecule Photon Trajectories

The rates of protein (un)folding are often described as diffusion on the projection of a hyperdimensional energy landscape onto a few (ideally one) order parameters. Testing such an approximation by experiment requires resolving the reactive transition paths of individual molecules, which is now becoming feasible with advanced single-molecule spectroscopic techniques. This has also sparked the interest of theorists in better understanding reactive transition paths. Here we focus on these issues aiming to establish (i) practical guidelines for the mechanistic interpretation of transition path times (TPT) and (ii) methods to extract the free energy surface and protein dynamics from the maximum likelihood analysis of photon trajectories (MLA-PT). We represent the (un)folding rates as diffusion on a 1D free energy surface with the FRET efficiency as a reaction coordinate proxy. We then perform diffusive kinetic simulations on surfaces with two minima and a barrier, but with different shapes (curvatures, barrier height, and symmetry), coupled to stochastic simulations of photon emissions that reproduce current SM-FRET experiments. From the analysis of transition paths, we find that the TPT is inversely proportional to the barrier height (difference in free energy between minimum and barrier top) for any given surface shape, and that dividing the TPT into climb and descent segments provides key information about the barrier's symmetry. We also find that the original MLA-PT procedure used to determine the TPT from experiments underestimates its value, particularly for the cases with smaller barriers (e.g., fast folders), and we suggest a simple strategy to correct for this bias. Importantly, we also demonstrate that photon trajectories contain enough information to extract the 1D free energy surface's shape and dynamics (if TPT is >4-5-fold longer than the interphoton time) using the MLA-PT directly implemented with a diffusive free energy surface model. When dealing with real (unknown) experimental data, the comparison between the likelihoods of the free energy surface and discrete kinetic three-state models can be used to evaluate the statistical significance of the estimated free energy surface.

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.

Anisotropy in mechanical unfolding of protein upon partner-assisted pulling and handle-assisted pulling

Proteins as force-sensors respond to mechanical cues and regulate signaling in physiology. Proteins commonly connect the source and response points of mechanical cues in two conformations, independent proteins in end-to-end geometry and protein complexes in handshake geometry. The force-responsive property of independent proteins in end-to-end geometry is studied extensively using single-molecule force spectroscopy (SMFS). The physiological significance of the complex conformations in force-sensing is often disregarded as mere surge protectors. However, with the potential of force-steering, protein complexes possess a distinct mechano-responsive property over individual force-sensors. To decipher, we choose a force-sensing protein, cadherin-23, from tip-link complex and perform SMFS using end-to-end geometry and handshake complex geometry. We measure higher force-resilience of cadherin-23 with preferential shorter extensions in handshake mode of pulling over the direct mode. The handshake geometry drives the force-response of cadherin-23 through different potential-energy landscapes than direct pulling. Analysis of the dynamic network structure of cadherin-23 under tension indicates narrow force-distributions among residues in cadherin-23 in direct pulling, resulting in low force-dissipation paths and low resilience to force. Overall, the distinct and superior mechanical responses of cadherin-23 in handshake geometry than single protein geometry highlight a probable evolutionary drive of protein-protein complexes as force-conveyors over independent ones.

Energy Landscape of Ubiquitin Is Weakly Multidimensional

Single molecule pulling experiments report time-dependent changes in the extension (X) of a biomolecule as a function of the applied force (f). By fitting the data to one-dimensional analytical models of the energy landscape, we can extract the hopping rates between the folded and unfolded states in two-state folders as well as the height and the location of the transition state (TS). Although this approach is remarkably insightful, there are cases for which the energy landscape is multidimensional (catch bonds being the most prominent). To assess if the unfolding energy landscape in small single domain proteins could be one-dimensional, we simulated force-induced unfolding of ubiquitin (Ub) using the coarse-grained self-organized polymer-side chain (SOP-SC) model. Brownian dynamics simulations using the SOP-SC model reveal that the Ub energy landscape is weakly multidimensional (WMD), governed predominantly by a single barrier. The unfolding pathway is confined to a narrow reaction pathway that could be described as diffusion in a quasi-1D X-dependent free energy profile. However, a granular analysis using the P_fold analysis, which does not assume any form for the reaction coordinate, shows that X alone does not account for the height and, more importantly, the location of the TS. The f-dependent TS location moves toward the folded state as f increases, in accord with the Hammond postulate. Our study shows that, in addition to analyzing the f-dependent hopping rates, the transition state ensemble must also be determined without resorting to X as a reaction coordinate to describe the unfolding energy landscapes of single domain proteins, especially if they are only WMD.

Side-Chain Polarity Modulates the Intrinsic Conformational Landscape of Model Dipeptides

The intrinsic conformational preferences of small peptides may provide additional insight into the thermodynamics and kinetics of protein folding. In this study, we explore the underlying energy landscapes of two model peptides, namely, Ac-Ala-NH₂ and Ac-Ser-NH₂, using geometry-optimization-based tools developed within the context of energy landscape theory. We analyze not only how side-chain polarity influences the structural preferences of the dipeptides, but also other emergent properties of the landscape, including heat capacity profiles, and kinetics of conformational rearrangements. The contrasting topographies of the free energy landscape agree with recent results from Fourier transform microwave spectroscopy experiments, where Ac-Ala-NH₂ was found to exist as a mixture of two conformers, while Ac-Ser-NH₂ remained structurally locked, despite exhibiting an apparently rich conformational landscape.

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.

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.

Folding and Misfolding of a PDZ Tandem Repeat

Although the vast majority of the human proteome is represented by multi-domain proteins, the study of multi-domain folding and misfolding is a relatively poorly explored field. The protein Whirlin is a multi-domain scaffolding protein expressed in the inner ear. It is characterized by the presence of tandem repeats of PDZ domains. The first two PDZ domains of Whirlin (PDZ1 and PDZ2 - namely P1P2) are structurally close and separated by a disordered short linker. We recently described the folding mechanism of the P1P2 tandem. The difference in thermodynamic stability of the two domains allowed us to selectively unfold one or both PDZ domains and to pinpoint the accumulation of a misfolded intermediate, which we demonstrated to retain physiological binding activity. In this work, we provide an extensive characterization of the folding and unfolding of P1P2. Based on the observed data, we describe an integrated kinetic analysis that satisfactorily fits the experiments and provides a valuable model to interpret multi-domain folding. The experimental and analytical approaches described in this study may be of general interest for the interpretation of complex multi-domain protein folding kinetics.

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.

Snapshotting the transient conformations and tracing the multiple pathways of single peptide folding using a solid-state nanopore

A fundamental question relating to protein folding/unfolding is the time evolution of the folding of a protein into its precisely defined native structure. The proper identification of transition conformations is essential for accurately describing the dynamic protein folding/unfolding pathways. Owing to the rapid transitions and sub-nm conformation differences involved, the acquisition of the transient conformations and dynamics of proteins is difficult due to limited instrumental resolution. Using the electrochemical confinement effect of a solid-state nanopore, we were able to snapshot the transient conformations and trace the multiple transition pathways of a single peptide inside a nanopore. By combining the results with a Markov chain model, this new single-molecule technique is applied to clarify the transition pathways of the β-hairpin peptide, which shows nonequilibrium fluctuations among several blockage current stages. This method enables the high-throughput investigation of transition pathways experimentally to access previously obscure peptide dynamics, which is significant for understanding the folding/unfolding mechanisms and misfolding of peptides or proteins.

A solid-state nanopore based method is described for resolving protein-folding-related problems via snapshotting the folding intermediates and characterizing the kinetics of a single peptide.

The road less traveled in protein folding: evidence for multiple pathways

Free Energy Landscape theory of Protein Folding, introduced over 20 years ago, implies that a protein has many paths to the folded conformation with the lowest free energy. Despite the knowledge in principle, it has been remarkably hard to detect such pathways. The lack of such observations is primarily due to the fact that no one experimental technique can detect many parts of the protein simultaneously with the time resolution necessary to see such differences in paths. However, recent technical developments and employment of multiple experimental probes and folding prompts have illuminated multiple folding pathways in a number of proteins that had all previously been described with a single path.

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.

Disordered proteins follow diverse transition paths as they fold and bind to a partner

Transition paths of macromolecular conformational changes such as protein folding are predicted to be heterogeneous. However, experimental characterization of the diversity of transition paths is extremely challenging because it requires measuring more than one distance during individual transitions. In this work, we used fast three-color single-molecule Förster resonance energy transfer spectroscopy to obtain the distribution of binding transition paths of a disordered protein. About half of the transitions follow a path involving strong non-native electrostatic interactions, resulting in a transition time of 300 to 800 microseconds. The remaining half follow more diverse paths characterized by weaker electrostatic interactions and more than 10 times shorter transition path times. The chain flexibility and non-native interactions make diverse binding pathways possible, allowing disordered proteins to bind faster than folded proteins.

Synonymous codon substitutions perturb cotranslational protein folding in vivo and impair cell fitness

In the cell, proteins are synthesized from N to C terminus and begin to fold during translation. Cotranslational folding mechanisms are therefore linked to elongation rate, which varies as a function of synonymous codon usage. However, synonymous codon substitutions can affect many distinct cellular processes, which has complicated attempts to deconvolve the extent to which synonymous codon usage can promote or frustrate proper protein folding in vivo. Although previous studies have shown that some synonymous changes can lead to different final structures, other substitutions will likely be more subtle, perturbing predominantly the protein folding pathway without radically altering the final structure. Here we show that synonymous codon substitutions encoding a single essential enzyme lead to dramatically slower cell growth. These mutations do not prevent active enzyme formation; instead, they predominantly alter the protein folding mechanism, leading to enhanced degradation in vivo. These results support a model in which synonymous codon substitutions can impair cell fitness by significantly perturbing cotranslational protein folding mechanisms, despite the chaperoning provided by the cellular protein homeostasis network.

Urea-aromatic interactions in biology

Noncovalent interactions are key determinants in both chemical and biological processes. Among such processes, the hydrophobic interactions play an eminent role in folding of proteins, nucleic acids, formation of membranes, protein-ligand recognition, etc.. Though this interaction is mediated through the aqueous solvent, the stability of the above biomolecules can be highly sensitive to any small external perturbations, such as temperature, pressure, pH, or even cosolvent additives, like, urea-a highly soluble small organic molecule utilized by various living organisms to regulate osmotic pressure. A plethora of detailed studies exist covering both experimental and theoretical regimes, to understand how urea modulates the stability of biological macromolecules. While experimentalists have been primarily focusing on the thermodynamic and kinetic aspects, theoretical modeling predominantly involves mechanistic information at the molecular level, calculating atomistic details applying the force field approach to the high level electronic details using the quantum mechanical methods. The review focuses mainly on examples with biological relevance, such as (1) urea-assisted protein unfolding, (2) urea-assisted RNA unfolding, (3) urea lesion interaction within damaged DNA, (4) urea conduction through membrane proteins, and (5) protein-ligand interactions those explicitly address the vitality of hydrophobic interactions involving exclusively the urea-aromatic moiety.

Thermal versus mechanical unfolding in a model protein

Force spectroscopy techniques are often used to learn about the free energy landscape of single biomolecules, typically by recovering free energy quantities that, extrapolated to zero force, are compared to those measured in bulk experiments. However, it is not always clear how the information obtained from a mechanically perturbed system can be related to the information obtained using other denaturants since tensioned molecules unfold and refold along a reaction coordinate imposed by the force, which is not likely to be meaningful in its absence. Here, we explore this dichotomy by investigating the unfolding landscape of a model protein, which is unfolded first mechanically through typical force spectroscopy-like protocols and next thermally. When unfolded by nonequilibrium force extension and constant force protocols, we recover a simple two-barrier landscape as the protein reaches the extended conformation through a metastable intermediate. Interestingly, folding-unfolding equilibrium simulations at low forces suggested a totally different scenario, where this metastable state plays little role in the unfolding mechanism, and the protein unfolds through two competing pathways [R. Tapia-Rojo et al., J. Chem. Phys. 141, 135102 (2014)]. Finally, we use Markov state models to describe the configurational space of the unperturbed protein close to the critical temperature. The thermal dynamics is well understood by a one-dimensional landscape along an appropriate reaction coordinate, however it is very different from the mechanical picture. In this sense, the results of our protein model for the mechanical and thermal descriptions provide incompatible views of the folding/unfolding landscape of the system, and the estimated quantities to zero force result are hard to interpret.

Principles for enhancing virus capsid capacity and stability from a thermophilic virus capsid structure

The capsids of double-stranded DNA viruses protect the viral genome from the harsh extracellular environment, while maintaining stability against the high internal pressure of packaged DNA. To elucidate how capsids maintain stability in an extreme environment, we use cryoelectron microscopy to determine the capsid structure of thermostable phage P74-26 to 2.8-Å resolution. We find P74-26 capsids exhibit an overall architecture very similar to those of other tailed bacteriophages, allowing us to directly compare structures to derive the structural basis for enhanced stability. Our structure reveals lasso-like interactions that appear to function like catch bonds. This architecture allows the capsid to expand during genome packaging, yet maintain structural stability. The P74-26 capsid has T = 7 geometry despite being twice as large as mesophilic homologs. Capsid capacity is increased with a larger, flatter major capsid protein. Given these results, we predict decreased icosahedral complexity (i.e. T ≤ 7) leads to a more stable capsid assembly.

Viral capsids need to protect the genome against harsh environmental conditions and cope with high internal pressure from the packaged genome. Here, the authors determine the structure of the thermostable phage P74-26 capsid at 2.8-Å resolution and identify features underlying enhanced capsid capacity and structural stability.

Diverse Folding Pathways of HIV-1 Protease Monomer on a Rugged Energy Landscape

The modern energy landscape theory of protein folding predicts multiple folding pathways connecting a myriad of unfolded conformations and a well-defined folded state. However, direct experimental observation of heterogeneous folding pathways is difficult. Naturally evolved proteins typically exhibit a smooth folding energy landscape for fast and efficient folding by avoiding unfavorable kinetic traps. In this case, rapid fluctuations between unfolded conformations result in apparent two-state behavior and make different pathways indistinguishable. However, the landscape roughness can be different, depending on the selection pressures during evolution. Here, we characterize the unusually rugged folding energy landscape of human immunodeficiency virus-1 protease monomer using single-molecule Förster resonance energy transfer spectroscopy. Our data show that fluctuations between unfolded conformations are slow, which enables the experimental observation of heterogeneous folding pathways as predicted by the landscape theory. Although the landscape ruggedness is sensitive to the mutations and fluorophore locations, the folding rate is similar for various protease constructs. The natural evolution of the protease to have a rugged energy landscape likely results from intrinsic pressures to maintain robust folding when human immunodeficiency virus-1 mutates frequently, which is essential for its survival.

A Plasmonic Approach to Study Protein Interaction Kinetics through the Dimerization of Functionalized Ag Nanoparticles

Understanding the kinetics of protein interactions plays a key role in biology with significant implications for the design of analytical methods for disease monitoring and diagnosis in medical care, research and industrial applications. Herein, we introduce a novel plasmonic approach to study the binding kinetics of protein-ligand interactions following the formation of silver nanoparticles (Ag NPs) dimers by UV-Vis spectroscopy that can be used as probes for antigen detection and quantification. To illustrate and test the method, the kinetics of the prototype biotin-streptavidin (Biot-STV) pair interaction was studied. Controlled aggregates (dimers) of STV functionalized Ag NPs were produced by adding stoichiometric quantities of gliadin-specific biotinylated antibodies (IgG-Biot). The dimerization kinetics was studied in a systematic way as a function of Ag NPs size and at different concentrations of IgG-Biot. The kinetics data have shown to be consistent with a complex reaction mechanism in which only the Ag NPs attached to the IgG-Biot located in a specific STV site are able to form dimers. These results help in elucidating a complex reaction mechanism involved in the dimerization kinetics of functionalized Ag NPs, which can serve as probes in surface plasmon resonance-based bioassays for the detection and quantification of different biomarkers or analytes of interest.

Mechanical unfolding of spectrin reveals a super-exponential dependence of unfolding rate on force

We investigated the mechanical unfolding of single spectrin molecules over a broad range of loading rates and thus unfolding forces by combining magnetic tweezers with atomic force microscopy. We find that the mean unfolding force increases logarithmically with loading rate at low loading rates, but the increase slows at loading rates above 1pN/s. This behavior indicates an unfolding rate that increases exponentially with the applied force at low forces, as expected on the basis of one-dimensional models of protein unfolding. At higher forces, however, the increase of the unfolding rate with the force becomes faster than exponential, which may indicate anti-Hammond behavior where the structures of the folded and transition states become more different as their free energies become more similar. Such behavior is rarely observed and can be explained by either a change in the unfolding pathway or as a reflection of a multidimensional energy landscape of proteins under force.

Nonnative structure in a peptide model of the unfolded state of superoxide dismutase 1 (SOD1): Implications for ALS-linked aggregation

Dozens of mutations throughout the sequence of the gene encoding superoxide dismutase 1 (SOD1) have been linked to toxic protein aggregation in the neurodegenerative disease amyotrophic lateral sclerosis (ALS). A parsimonious explanation for numerous genotypes resulting in a common phenotype would be mutation-induced perturbation of the folding free-energy surface that increases the populations of high-energy states prone to aggregation. The absence of intermediates in the folding of monomeric SOD1 suggests that the unfolded ensemble is a potential source of aggregation. To test this hypothesis, here we dissected SOD1 into a set of peptides end-labeled with FRET probes to model the local behavior of the corresponding sequences in the unfolded ensemble. Using time-resolved FRET, we observed that the peptide corresponding to the Loop VII-β8 sequence at the SOD1 C terminus was uniquely sensitive to denaturant. Utilizing a two-dimensional form of maximum entropy modeling, we demonstrate that the sensitivity to denaturant is the surprising result of a two-state-like transition from a compact to an expanded state. Variations of the peptide sequence revealed that the compact state involves a nonnative interaction between the disordered N terminus and the hydrophobic C terminus of the peptide. This nonnative intramolecular structure could serve as a precursor for intermolecular association and result in aggregation associated with ALS. We propose that this precursor would provide a common molecular target for therapeutic intervention in the dozens of ALS-linked SOD1 mutations.

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.

Measuring the average shape of transition paths during the folding of a single biological molecule

Transition paths represent the parts of a reaction where the energy barrier separating products and reactants is crossed. They are essential to understanding reaction mechanisms, yet many of their properties remain unstudied. Here, we report measurements of the average shape of transition paths, studying the folding of DNA hairpins as a model system for folding reactions. Individual transition paths were detected in the folding trajectories of hairpins with different sequences held under tension in optical tweezers, and path shapes were computed by averaging all transitions in the time domain, 〈t(x)〉, or by averaging transitions of a given duration in the extension domain, 〈x(t|τ)〉 _τ Whereas 〈t(x)〉 was close to straight, with only a subtle curvature, 〈x(t|τ)〉 _τ had more pronounced curvature that fit well to theoretical expectations for the dominant transition path, returning diffusion coefficients similar to values obtained previously from independent methods. Simulations suggested that 〈t(x)〉 provided a less reliable representation of the path shape than 〈x(t|τ)〉 _τ , because it was far more sensitive to the effects of coupling the molecule to the experimental force probe. Intriguingly, the path shape variance was larger for some hairpins than others, indicating sequence-dependent changes in the diversity of transition paths reflective of differences in the character of the energy barriers, such as the width of the barrier saddle-point or the presence of parallel paths through multiple barriers between the folded and unfolded states. These studies of average path shapes point the way forward for probing the rich information contained in path shape fluctuations.

Nature's Shortcut to Protein Folding

This Feature Article presents a view of the protein folding transition based on the hypothesis that Nature has built features within the sequences that enable a Shortcut to efficient folding. Nature's Shortcut is proposed to be the early establishment of a set of nonlocal weak contacts, constituting protein loops that significantly constrain regions of the collapsed disordered protein into a native-like low-resolution fluctuating topology of major sections of the backbone. Nature's establishment of this scaffold of nonlocal contacts is claimed to bypass what would otherwise be a nearly hopeless unaided search for the final three-dimensional structure in proteins longer than ∼100 amino acids. To support this main contention of the Feature Article, the loop hypothesis (LH) description of early folding events is experimentally tested with time-resolved Förster resonance energy transfer techniques for adenylate kinase, and the data are shown to be consistent with theoretical predictions from the sequential collapse model (SCM). The experimentally based LH and the theoretically founded SCM are argued to provide a unified picture of the role of nonlocal contacts as constituting Nature's Shortcut to protein folding. Importantly, the SCM is shown to reliably predict key nonlocal contacts utilizing only primary sequence information. This view on Nature's Shortcut is open to the protein community for further detailed assessment, including its practical consequences, by suitable application of advanced experimental and computational techniques.

A small single-domain protein folds through the same pathway on and off the ribosome

In vivo, proteins fold and function in a complex environment subject to many stresses that can modulate a protein's energy landscape. One aspect of the environment pertinent to protein folding is the ribosome, since proteins have the opportunity to fold while still bound to the ribosome during translation. We use a combination of force and chemical denaturant (chemomechanical unfolding), as well as point mutations, to characterize the folding mechanism of the src SH3 domain both as a stalled ribosome nascent chain and free in solution. Our results indicate that src SH3 folds through the same pathway on and off the ribosome. Molecular simulations also indicate that the ribosome does not affect the folding pathway for this small protein. Taken together, we conclude that the ribosome does not alter the folding mechanism of this small protein. These results, if general, suggest the ribosome may exert a bigger influence on the folding of multidomain proteins or protein domains that can partially fold before the entire domain sequence is outside the ribosome exit tunnel.

Collapse Transitions of Proteins and the Interplay Among Backbone, Sidechain, and Solvent Interactions

Proteins can collapse into compact globules or form expanded, solvent-accessible, coil-like conformations. Additionally, they can fold into well-defined three-dimensional structures or remain partially or entirely disordered. Recent discoveries have shown that the tendency for proteins to collapse or remain expanded is not intrinsically coupled to their ability to fold. These observations suggest that proteins do not have to form compact globules in aqueous solutions. They can be intrinsically disordered, collapsed, or expanded, and even form well-folded, elongated structures. This ability to decouple collapse from folding is determined by the sequence details of proteins. In this review, we highlight insights gleaned from studies over the past decade. Using a polymer physics framework, we explain how the interplay among sidechains, backbone units, and solvent determines the driving forces for collapsed versus expanded states in aqueous solvents.

Extracting the Hidden Distributions Underlying the Mean Transition State Structures in Protein Folding

The inherent conflict between noncovalent interactions and the large conformational entropy of the polypeptide chain forces folding reactions and their mechanisms to deviate significantly from chemical reactions. Accordingly, measures of structure in the transition state ensemble (TSE) are strongly influenced by the underlying distributions of microscopic folding pathways that are challenging to discern experimentally. Here, we present a detailed analysis of 150,000 folding transition paths of five proteins at three different thermodynamic conditions from an experimentally consistent statistical mechanical model. We find that the underlying TSE structural distributions are rarely unimodal, and the average experimental measures arise from complex underlying distributions. Unfolding pathways also exhibit subtle differences from folding counterparts due to a combination of Hammond behavior and native-state movements. Local interactions and topological complexity, to a lesser extent, are found to determine pathway heterogeneity, underscoring the importance of the balance between local and nonlocal energetics in protein folding.

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.

Two Disease-Causing SNAP-25B Mutations Selectively Impair SNARE C-terminal Assembly

Synaptic exocytosis relies on assembly of three soluble N-ethylmaleimide-sensitive factor attachment receptor (SNARE) proteins into a parallel four-helix bundle to drive membrane fusion. SNARE assembly occurs by stepwise zippering of the vesicle-associated SNARE (v-SNARE) onto a binary SNARE complex on the target plasma membrane (t-SNARE). Zippering begins with slow N-terminal association followed by rapid C-terminal zippering, which serves as a power stroke to drive membrane fusion. SNARE mutations have been associated with numerous diseases, especially neurological disorders. It remains unclear how these mutations affect SNARE zippering, partly due to difficulties to quantify the energetics and kinetics of SNARE assembly. Here, we used single-molecule optical tweezers to measure the assembly energy and kinetics of SNARE complexes containing single mutations I67T/N in neuronal SNARE synaptosomal-associated protein of 25kDa (SNAP-25B), which disrupt neurotransmitter release and have been implicated in neurological disorders. We found that both mutations significantly reduced the energy of C-terminal zippering by ~10 k_BT, but did not affect N-terminal assembly. In addition, we observed that both mutations lead to unfolding of the C-terminal region in the t-SNARE complex. Our findings suggest that both SNAP-25B mutations impair synaptic exocytosis by destabilizing SNARE assembly, rather than stabilizing SNARE assembly as previously proposed. Therefore, our measurements provide insights into the molecular mechanism of the disease caused by SNARE mutations.

Exploring the Denatured State Ensemble by Single-Molecule Chemo-Mechanical Unfolding: The Effect of Force, Temperature, and Urea

While it is widely appreciated that the denatured state of a protein is a heterogeneous conformational ensemble, there is still debate over how this ensemble changes with environmental conditions. Here, we use single-molecule chemo-mechanical unfolding, which combines force and urea using the optical tweezers, together with traditional protein unfolding studies to explore how perturbants commonly used to unfold proteins (urea, force, and temperature) affect the denatured-state ensemble. We compare the urea m-values, which report on the change in solvent accessible surface area for unfolding, to probe the denatured state as a function of force, temperature, and urea. We find that while the urea- and force-induced denatured states expose similar amounts of surface area, the denatured state at high temperature and low urea concentration is more compact. To disentangle these two effects, we use destabilizing mutations that shift the T_m and C_m. We find that the compaction of the denatured state is related to changing temperature as the different variants of acyl-coenzyme A binding protein have similar m-values when they are at the same temperature but different urea concentration. These results have important implications for protein folding and stability under different environmental conditions.

Proteins Breaking Bad: A Free Energy Perspective

Protein aging may manifest as a mechanical disease that compromises tissue elasticity. As proved recently, while proteins respond to changes in force with an instantaneous elastic recoil followed by a folding contraction, aged proteins break bad, becoming unstructured polymers. Here, we explain this phenomenon in the context of a free energy model, predicting the changes in the folding landscape of proteins upon oxidative aging. Our findings validate that protein folding under force is constituted by two separable components, polymer properties and hydrophobic collapse, and demonstrate that the latter becomes irreversibly blocked by oxidative damage. We run Brownian dynamics simulations on the landscape of protein L octamer, reproducing all experimental observables, for a naive and damaged polyprotein. This work provides a unique tool to understand the evolving free energy landscape of elastic proteins upon physiological changes, opening new perspectives to predict age-related diseases in tissues.

Structure-Based Prediction of Protein-Folding Transition Paths

We propose a general theory to describe the distribution of protein-folding transition paths. We show that transition paths follow a predictable sequence of high-free-energy transient states that are separated by free-energy barriers. Each transient state corresponds to the assembly of one or more discrete, cooperative units, which are determined directly from the native structure. We show that the transition state on a folding pathway is reached when a small number of critical contacts are formed between a specific set of substructures, after which folding proceeds downhill in free energy. This approach suggests a natural resolution for distinguishing parallel folding pathways and provides a simple means to predict the rate-limiting step in a folding reaction. Our theory identifies a common folding mechanism for proteins with diverse native structures and establishes general principles for the self-assembly of polymers with specific interactions.

Kinetic Ductility and Force-Spike Resistance of Proteins from Single-Molecule Force Spectroscopy

Ductile materials can absorb spikes in mechanical force, whereas brittle ones fail catastrophically. Here we develop a theory to quantify the kinetic ductility of single molecules from force spectroscopy experiments, relating force-spike resistance to the differential responses of the intact protein and the unfolding transition state to an applied mechanical force. We introduce a class of unistable one-dimensional potential surfaces that encompass previous models as special cases and continuously cover the entire range from ductile to brittle. Compact analytic expressions for force-dependent rates and rupture-force distributions allow us to analyze force-clamp and force-ramp pulling experiments. We find that the force-transmitting protein domains of filamin and titin are kinetically ductile when pulled from their two termini, making them resistant to force spikes. For the mechanostable muscle protein titin, a highly ductile model reconciles data over 10 orders of magnitude in force loading rate from experiment and simulation.

Meandering Down the Energy Landscape of Protein Folding: Are We There Yet?

As judged by a single publication metric, the activity in the protein folding field has been declining over the past 5 years, after enjoying a decade-long growth. Does this development indicate that the field is sunsetting or is this decline only temporary? Upon surveying a small territory of its landscape, we find that the protein folding field is still quite active and many important findings have emerged from recent experimental studies. However, it is also clear that only continued development of new techniques and methods, especially those enabling dissection of the fine details and features of the protein folding energy landscape, will fuel this old field to move forward.