A simple two-state protein unfolds mechanically via multiple heterogeneous pathways at single-molecule resolution

Identical sequences, different behaviors: Protein diversity captured at the single-molecule level

The classical "one sequence, one structure, one function" paradigm has shaped much of our intuition of how proteins work inside the cell. Partially due to the insight provided by bulk biochemical assays, individual biomolecules are often assumed to behave as identical entities, and their characterization relies on ensemble averages that flatten any conformational diversity into a unique phenotype. While the emergence of single-molecule techniques opened the gates to interrogating individual molecules, technical shortcomings typically limit the duration of these measurements, which precludes a complete characterization of an individual protein and, hence, capturing the heterogeneity among molecular populations. Here, we introduce an ultrastable magnetic tweezers design, which enables us to measure the folding dynamics of a single protein during several uninterrupted days with high temporal and spatial resolution. Thanks to this instrumental development, we fully characterize the nanomechanics of two proteins with a very distinct force response, the talin R3IVVI domain and protein L. Days-long recordings on the same protein individual accumulate thousands of folding transitions with submicrosecond resolution, allowing us to reconstruct their free energy landscapes and describe how they evolve with force. By mapping the nanomechanical identity of many different protein individuals, we directly capture their molecular diversity as a quantifiable dispersion on their force response and folding kinetics. By significantly expanding the measurable timescales, our instrumental development offers a tool for profiling individual molecules, opening the gates to directly characterizing biomolecular heterogeneity.

Polyols, increasing global stability of cytochrome c, destabilize the thermal unfolding intermediate

Studies on the intermediate states of proteins provide essential information on folding pathway and energy landscape of proteins. Osmolytes, known to alter the stability of proteins, might also affect the structure and energy states of folding intermediates. This was examined using cytochrome c (Cyt) as a model protein which forms a spectroscopically detectable intermediate during thermal denaturation transition. Most of the secondary structure and the native heme-ligation were intact in the intermediate state of the protein. Denaturants, urea and guanidinium hydrochloride, and ionic salt destabilizes the intermediate and drive the protein to follow two-state transition. The effect of polyol class of osmolytes, glycol, glycerol, erythritol, xylitol and sorbitol (with OH-groups two to six), on the intermediate was studied using Soret absorbance and far-UV circular dichroism. With the increasing concentration of any of the polyols, the transition-midpoint temperature (T_m) and the enthalpy change (ΔH) for native to intermediate transition were decreased. This indicated that the intermediate was destabilized by the polyols. However, the polyols increased the overall stability of the protein by increasing T_m and ΔH for intermediate to unfolded transition, except for glycol which destabilized the protein. These results show that the polyols could alter the energy state of the intermediate, and the effect of lower and higher polyols might be different on the stability and folding pathway of the protein.Communicated by Ramaswamy H. Sarma.

Protein folding in vitro and in the cell: From a solitary journey to a team effort

Correct protein folding is essential for the health and function of living organisms. Yet, it is not well understood how unfolded proteins reach their native state and avoid aggregation, especially within the cellular milieu. Some proteins, especially small, single-domain and apparent two-state folders, successfully attain their native state upon dilution from denaturant. Yet, many more proteins undergo misfolding and aggregation during this process, in a concentration-dependent fashion. Once formed, native and aggregated states are often kinetically trapped relative to each other. Hence, the early stages of protein life are absolutely critical for proper kinetic channeling to the folded state and for long-term solubility and function. This review summarizes current knowledge on protein folding/aggregation mechanisms in buffered solution and within the bacterial cell, highlighting early stages. Remarkably, teamwork between nascent chain, ribosome, trigger factor and Hsp70 molecular chaperones enables all proteins to overcome aggregation propensities and reach a long-lived bioactive state.

Using steered molecular dynamic tension for assessing quality of computational protein structure models

The native structures of proteins, except for notable exceptions of intrinsically disordered proteins, in general take their most stable conformation in the physiological condition to maintain their structural framework so that their biological function can be properly carried out. Experimentally, the stability of a protein can be measured by several means, among which the pulling experiment using the atomic force microscope (AFM) stands as a unique method. AFM directly measures the resistance from unfolding, which can be quantified from the observed force-extension profile. It has been shown that key features observed in an AFM pulling experiment can be well reproduced by computational molecular dynamics simulations. Here, we applied computational pulling for estimating the accuracy of computational protein structure models under the hypothesis that the structural stability would positively correlated with the accuracy, i.e. the closeness to the native, of a model. We used in total 4929 structure models for 24 target proteins from the Critical Assessment of Techniques of Structure Prediction (CASP) and investigated if the magnitude of the break force, that is, the force required to rearrange the model's structure, from the force profile was sufficient information for selecting near-native models. We found that near-native models can be successfully selected by examining their break forces suggesting that high break force indeed indicates high stability of models. On the other hand, there were also near-native models that had relatively low peak forces. The mechanisms of the stability exhibited by the break forces were explored and discussed.

Molecular interpretation of single-molecule force spectroscopy experiments with computational approaches

Single molecule force-spectroscopy techniques have granted access to unprecedented molecular-scale details about biochemical and biological mechanisms. However, the interpretation of the experimental data is often challenging. Computational and simulation approaches (all-atom steered MD simulations in particular) are key to provide molecular details about the associated mechanisms, to help test different hypotheses and to predict experimental results. In this review, particular recent efforts directed towards the molecular interpretation of single-molecule force spectroscopy experiments on proteins and protein-related systems (often in close collaboration with experimental groups) will be presented. These results will be discussed in the broader context of the field, highlighting the recent achievements and the ongoing challenges for computational biophysicists and biochemists. In particular, I will focus on the input gained from molecular simulations approaches to rationalize the origin of the unfolded protein elasticity and the protein conformational behavior under force, to understand how force denaturation differs from chemical, thermal or shear unfolding, and to unravel the molecular details of unfolding events for a variety of systems. I will also discuss the use of models based on Langevin dynamics on a 1-D free-energy surface to understand the effect of protein segmentation on the work exerted by a force, or, at the other end of the spectrum of computational techniques, how quantum calculations can help to understand the reactivity of disulfide bridges exposed to force.

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.

Two energy barriers and a transient intermediate state determine the unfolding and folding dynamics of cold shock protein

Cold shock protein (Csp) is a typical two-state folding model protein which has been widely studied by biochemistry and single molecule techniques. Recently two-state property of Csp was confirmed by atomic force microscopy (AFM) through direct pulling measurement, while several long-lifetime intermediate states were found by force-clamp AFM. We systematically studied force-dependent folding and unfolding dynamics of Csp using magnetic tweezers with intrinsic constant force capability. Here we report that Csp mostly folds and unfolds with a single step over force range from 5 pN to 50 pN, and the unfolding rates show different force sensitivities at forces below and above ~8 pN, which determines a free energy landscape with two barriers and a transient intermediate state between them along one transition pathway. Our results provide a new insight on protein folding mechanism of two-state proteins.

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.

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.

Highly Dynamic Polynuclear Metal Cluster Revealed in a Single Metallothionein Molecule

Human metallothionein (MT) is a small-size yet efficient metal-binding protein, playing an essential role in metal homeostasis and heavy metal detoxification. MT contains two domains, each forming a polynuclear metal cluster with an exquisite hexatomic ring structure. The apoprotein is intrinsically disordered, which may strongly influence the clusters and the metal-thiolate (M-S) bonds, leading to a highly dynamic structure. However, these features are challenging to identify due to the transient nature of these species. The individual signal from dynamic conformations with different states of the cluster and M-S bond will be averaged and blurred in classic ensemble measurement. To circumvent these problems, we combined a single-molecule approach and multiscale molecular simulations to investigate the rupture mechanism and chemical stability of the metal cluster by a single MT molecule, focusing on the Zn₄S₁₁ cluster in the α domain upon unfolding. Unusual multiple unfolding pathways and intermediates are observed for both domains, corresponding to different combinations of M-S bond rupture. None of the pathways is clearly preferred suggesting that unfolding proceeds from the distribution of protein conformational substates with similar M-S bond strengths. Simulations indicate that the metal cluster may rearrange, forming and breaking metal-thiolate bonds even when MT is folded independently of large protein backbone reconfiguration. Thus, a highly dynamic polynuclear metal cluster with multiple conformational states is revealed in MT, responsible for the binding promiscuity and diverse cellular functions of this metal-carrier protein.

Mobile Protons Limit the Stability of Salt Bridges in the Gas Phase: Implications for the Structures of Electrosprayed Protein Ions

Electrosprayed protein ions can retain native-like conformations. The intramolecular contacts that stabilize these compact gas-phase structures remain poorly understood. Recent work has uncovered abundant salt bridges in electrosprayed proteins. Salt bridges are zwitterionic BH⁺/A^- contacts. The low dielectric constant in the vacuum strengthens electrostatic interactions, suggesting that salt bridges could be a key contributor to the retention of compact protein structures. A problem with this assertion is that H⁺ are mobile, such that H⁺ transfer can convert salt bridges into neutral B⁰/HA⁰ contacts. This possible salt bridge annihilation puts into question the role of zwitterionic motifs in the gas phase, and it calls for a detailed analysis of BH⁺/A^- versus B⁰/HA⁰ interactions. Here, we investigate this issue using molecular dynamics (MD) simulations and electrospray experiments. MD data for short model peptides revealed that salt bridges with static H⁺ have dissociation energies around 700 kJ mol^-1. The corresponding B⁰/HA⁰ contacts are 1 order of magnitude weaker. When considering the effects of mobile H⁺, BH⁺/A^- bond energies were found to be between these two extremes, confirming that H⁺ migration can significantly weaken salt bridges. Next, we examined the protein ubiquitin under collision-induced unfolding (CIU) conditions. CIU simulations were conducted using three different MD models: (i) Positive-only runs with static H⁺ did not allow for salt bridge formation and produced highly expanded CIU structures. (ii) Zwitterionic runs with static H⁺ resulted in abundant salt bridges, culminating in much more compact CIU structures. (iii) Mobile H⁺ simulations allowed for the dynamic formation/annihilation of salt bridges, generating CIU structures intermediate between scenarios (i) and (ii). Our results uncover that mobile H⁺ limit the stabilizing effects of salt bridges in the gas phase. Failure to consider the effects of mobile H⁺ in MD simulations will result in unrealistic outcomes under CIU conditions.

Cold-Shock Domains-Abundance, Structure, Properties, and Nucleic-Acid Binding

Simple Summary

Proteins are composed of compact domains, often of known three-dimensional structure, and natively unstructured polypeptide regions. The abundant cold-shock domain is among the set of canonical nucleic acid-binding domains and conserved from bacteria to man. Proteins containing cold-shock domains serve a large variety of biological functions, which are mostly linked to DNA or RNA binding. These functions include the regulation of transcription, RNA splicing, translation, stability and sequestration. Cold-shock domains have a simple architecture with a conserved surface ideally suited to bind single-stranded nucleic acids. Because the binding is mostly by non-specific molecular interactions which do not involve the sugar-phosphate backbone, cold-shock domains are not strictly sequence-specific and do not discriminate reliably between DNA and RNA. Many, but not all functions of cold shock-domain proteins in health and disease can be understood based of the physical and structural properties of their cold-shock domains.

Abstract

The cold-shock domain has a deceptively simple architecture but supports a complex biology. It is conserved from bacteria to man and has representatives in all kingdoms of life. Bacterial cold-shock proteins consist of a single cold-shock domain and some, but not all are induced by cold shock. Cold-shock domains in human proteins are often associated with natively unfolded protein segments and more rarely with other folded domains. Cold-shock proteins and domains share a five-stranded all-antiparallel β-barrel structure and a conserved surface that binds single-stranded nucleic acids, predominantly by stacking interactions between nucleobases and aromatic protein sidechains. This conserved binding mode explains the cold-shock domains’ ability to associate with both DNA and RNA strands and their limited sequence selectivity. The promiscuous DNA and RNA binding provides a rationale for the ability of cold-shock domain-containing proteins to function in transcription regulation and DNA-damage repair as well as in regulating splicing, translation, mRNA stability and RNA sequestration.

Lessons about Protein Folding and Binding from Archetypal Folds

The function of proteins as biological nanomachines relies on their ability to fold into complex 3D structures, bind selectively to partners, and undergo conformational changes on cue. The native functional structures, and the rates of interconversion between conformational states (folded-unfolded, bound-free), are all encoded in the physical chemistry of their amino acid sequence. However, despite extensive research over decades, this code has proven difficult to fully crack, in terms of both prediction and understanding the molecular mechanisms at play.Earlier work on single-domain proteins reported a commonality of slow rates (10^-2-10² s^-1) and simple behavior in both kinetic and thermodynamic unfolding experiments, which suggested the process was all-or-none and thereby analogous to a chemical reaction (e.g., A ⇄ B). In the absence of a first-principles pre-exponential factor for protein (un)folding dynamics, the rates could only be interpreted in relative terms, e.g., the changes induced by mutation, and hence, neither the height of nor the entropic contribution to the free energy barriers was known. The rates were also many orders of magnitude too slow for direct atomistic simulations, and the computational focus was on predicting rate changes induced by mutation via coarse grained simulations. However, even the effects of mutation proved to be strikingly homogeneous with all experimental data clustering at ∼1/3 of the free energy perturbation recovered on folding and ∼2/3 on unfolding.The implementation of ultrafast kinetic methods turned the field upside down because they allowed researchers to measure the time scales of elementary (un)folding motions, which set the pre-exponential factor for protein conformational transitions at ∼1 μs. In parallel, we and others set out to investigate the simplest possible protein structures capable of autonomous folding, which we defined as archetypal folds. The rationale was to recapitulate the hierarchical organization of protein structure, starting from the bottom up. The study of fold archetypes ended up opening new research avenues in protein (un)folding, but also making unexpected connections with the folding upon binding of intrinsically disordered proteins and suggesting their functioning as conformational rheostats.This Account describes our work on the kinetic, thermodynamic, mechanistic, and functional analysis of fold archetypes. We first discuss the kinetic studies, emphasizing their impact on our understanding of (un)folding rates, of barrierless (downhill) folding, and as benchmarks for atomistic simulations. We continue with the thermodynamic analysis, introducing the differential scanning calorimetry, multiprobe, and NMR approaches that we developed to dissect their gradual, minimally cooperative (un)folding transitions and to probe the underlying mechanisms with unprecedented detail. The last two sections cover single-molecule analyses and some recent, mostly computational, results on the exploration of possible biological and technological roles for the gradual conformational transitions of fold archetypes.

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.

Next Generation Methods for Single-Molecule Force Spectroscopy on Polyproteins and Receptor-Ligand Complexes

Single-molecule force spectroscopy with the atomic force microscope provides molecular level insights into protein function, allowing researchers to reconstruct energy landscapes and understand functional mechanisms in biology. With steadily advancing methods, this technique has greatly accelerated our understanding of force transduction, mechanical deformation, and mechanostability within single- and multi-domain polyproteins, and receptor-ligand complexes. In this focused review, we summarize the state of the art in terms of methodology and highlight recent methodological improvements for AFM-SMFS experiments, including developments in surface chemistry, considerations for protein engineering, as well as theory and algorithms for data analysis. We hope that by condensing and disseminating these methods, they can assist the community in improving data yield, reliability, and throughput and thereby enhance the information that researchers can extract from such experiments. These leading edge methods for AFM-SMFS will serve as a groundwork for researchers cognizant of its current limitations who seek to improve the technique in the future for in-depth studies of molecular biomechanics.

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.

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.

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

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

Multistep Protein Unfolding Scenarios from the Rupture of a Complex Metal Cluster Cd3S9

Protein (un)folding is a complex and essential process. With the rapid development of single-molecule techniques, we can detect multiple and transient proteins (un)folding pathways/intermediates. However, the observation of multiple multistep (>2) unfolding scenarios for a single protein domain remains limited. Here, we chose metalloprotein with relatively stable and multiple metal-ligand coordination bonds as a system for such a purpose. Using AFM-based single-molecule force spectroscopy (SMFS), we successfully demonstrated the complex and multistep protein unfolding scenarios of the β-domain of a human protein metallothionein-3 (MT). MT is a protein of ~60 amino acids (aa) in length with 20 cysteines for various metal binding, and the β-domain (βMT) is of ~30 aa with an M3S9 metal cluster. We detected four different types of three-step protein unfolding scenarios from the Cd-βMT, which can be possibly explained by the rupture of Cd-S bonds in the complex Cd3S9 metal cluster. In addition, complex unfolding scenarios with four rupture peaks were observed. The Cd-S bonds ruptured in both single bond and multiple bonds modes. Our results provide not only evidence for multistep protein unfolding phenomena but also reveal unique properties of metalloprotein system using single-molecule AFM.

Fodis: Software for Protein Unfolding Analysis

The folding dynamics of proteins at the single-molecule level has been studied with single-molecule force spectroscopy experiments for 20 years, but a common standardized method for the analysis of the collected data and for sharing among the scientific community members is still not available. We have developed a new open-source tool-Fodis-for the analysis of the force-distance curves obtained in single-molecule force spectroscopy experiments, providing almost automatic processing, analysis, and classification of the obtained data. Our method provides also a classification of the possible unfolding pathways and the structural heterogeneity present during the unfolding of proteins.

Three Weaknesses for Three Perturbations: Comparing Protein Unfolding Under Shear, Force, and Thermal Stresses

The perturbation of a protein conformation by a physiological fluid flow is crucial in various biological processes including blood clotting and bacterial adhesion to human tissues. Investigating such mechanisms by computer simulations is thus of great interest, but it requires development of ad hoc strategies to mimic the complex hydrodynamic interactions acting on the protein from the surrounding flow. In this study, we apply the Lattice Boltzmann Molecular Dynamics (LBMD) technique built on the implicit solvent coarse-grained model for protein Optimized Potential for Efficient peptide structure Prediction (OPEP) and a mesoscopic representation of the fluid solvent, to simulate the unfolding of a small globular cold-shock protein in shear flow and to compare it to the unfolding mechanisms caused either by mechanical or thermal perturbations. We show that each perturbation probes a specific weakness of the protein and causes the disruption of the native fold along different unfolding pathways. Notably, the shear flow and the thermal unfolding exhibit very similar pathways, while because of the directionality of the perturbation, the unfolding under force is quite different. For force and thermal disruption of the native state, the coarse-grained simulations are compared to all-atom simulations in explicit solvent, showing an excellent agreement in the explored unfolding mechanisms. These findings encourage the use of LBMD based on the OPEP model to investigate how a flow can affect the function of larger proteins, for example, in catch-bond systems.

Single glucose molecule transport process revealed by force tracing and molecular dynamics simulations

Transporting individual molecules across cell membranes is a fundamental process in cellular metabolism. Although the crystal diffraction technique has greatly contributed to our understanding of the structures of the involved transporters, a description of the dynamic transport mechanism at the single-molecule level has been extremely elusive. In this study, we applied atomic force microscopy (AFM)-based force tracing to directly monitor the transport of a single molecule, d-glucose, across living cell membranes. Our results show that the force to transport a single molecule of d-glucose across cell membranes is 37 ± 9 pN, and the corresponding transport interval is approximately 20 ms, while the average speed is approximately 0.3 μm s^-1. Furthermore, our calculated force profile from molecular dynamics simulations showed quantitatively good agreement with the force tracing observation and revealed detailed information regarding the glucose transport path, indicating that two salt bridges, K38/E299 and K300/E426, play critical roles during glucose transport across glucose transporter 1 (GLUT1). This role was further verified using biological experiments that disrupted these two bridges and measured the uptake of glucose into the cells. Our approaches led to the first unambiguous description of the glucose transport process across cell membranes at the single-molecule level and demonstrated the biological importance of the two salt bridges for transporting glucose across GLUT1.

Minichaperone (GroEL191-345) mediated folding of MalZ proceeds by binding and release of native and functional intermediates

The isolated apical domain of GroEL consisting of residues 191-345 (known as "minichaperone") binds and assists the folding of a wide variety of client proteins without GroES and ATP, but the mechanism of its action is still unknown. In order to probe into the matter, we have examined minichaperone-mediated folding of a large aggregation prone protein Maltodextrin-glucosidase (MalZ). The key objective was to identify whether MalZ exists free in solution, or remains bound to, or cycling on and off the minichaperone during the refolding process. When GroES was introduced during refolding process, production of the native MalZ was inhibited. We also observed the same findings with a trap mutant of GroEL, which stably captures a predominantly non-native MalZ released from minichaperone during refolding process, but does not release it. Tryptophan and ANS fluorescence measurements indicated that refolded MalZ has the same structure as the native MalZ, but that its structure when bound to minichaperone is different. Surface plasmon resonance measurements provide an estimate for the equilibrium dissociation constant KD for the MalZ-minichaperone complex of 0.21 ± 0.04 μM, which are significantly higher than for most GroEL clients. This showed that minichaperone interacts loosely with MalZ to allow the protein to change its conformation and fold while bound during the refolding process. These observations suggest that the minichaperone works by carrying out repeated cycles of binding aggregation-prone protein MalZ in a relatively compact conformation and in a partially folded but active state, and releasing them to attempt to fold in solution.

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.

Differences in the mechanical unfolding pathways of apo- and copper-bound azurins

Metalloproteins carry out diverse biological functions including metal transport, electron transfer, and catalysis. At present, the influence of metal cofactors on metalloprotein stability is not well understood. Here, we report the mechanical stability and unfolding pathway of azurin, a cupredoxin family protein with β-barrel topology and type I copper-binding centre. Single-molecule force spectroscopy (SMFS) experiments reveal 2-state and 3-state unfolding pathways for apo-azurin. The intermediate in the 3-state pathway occurs at an unfolding contour length of 7.5 nm from the native state. Steered molecular dynamics (SMD) simulations show that apo-azurin unfolds via a first transition state (TS) where β2Β–β8 and β7–β8 strand pairs rupture to form the intermediate, which subsequently unfolds by the collective rupture of remaining strands. SMFS experiments on holo-azurin exhibit an additional 4-state pathway besides the 2-state and 3-state pathways. The unfolding contour length leading to the first intermediate is 6.7 nm suggesting a sequestration of ~1 nm polypeptide chain length by the copper. SMD simulations reveal atomistic details of the copper sequestration and predict a combined β4–β7 pair and copper coordination sphere rupture to create the third TS in the 4-state pathway. Our systematic studies provide detailed mechanistic insights on modulation of protein mechanical properties by metal-cofactors.

Mechanics of Protein Adaptation to High Temperatures

Inspired by Somero's corresponding state principle that relates protein enhanced thermal stability with mechanical rigidity, we deployed state of the art computational techniques (based on atomistic steered molecular dynamics and Hamiltonian-replica exchange simulations) to study the in silico realization of mechanical and thermal unfolding of two homologous Csp proteins that have evolved to thrive in different thermal environments. By complementing recent single-molecule experiments, we unambiguously show that, for these homologues whose structures are very similar, the increased thermal resistance of the thermophilic variant is not associated with an increased mechanical stability. Our approach provides microscopic insights that are otherwise inaccessible to experimental techniques, and explains why the protein weak spots for thermal and mechanical denaturation are distinct.

Hydrophilicities of amylose and natural cellulose are regulated by the linkage between sugar rings

Comparative studies of single molecule force spectroscopy and molecular dynamics simulations indicate that natural cellulose is more hydrophobic than amylose at the single-chain level, implying that the hydrophobicities of these polymeric isomers are regulated by only one parameter in the chains, the linkage between the sugar rings.

Role of the hydrophobic core in cytoskeleton protein: cardiac myosin binding protein C

Cardiac myosin binding protein C is the object of analysis presented in this paper. The fuzzy oil drop model was applied to analyze the status of the hydrophobic core in two forms of this protein: WT and R502W mutant. The status of the mutant is revealed to be of lower stability than the WT form. The high order of the hydrophobic core is interpreted as the factor of stability of the tertiary structure. The muscle proteins, which undergo significant structural changes as the consequence of external stretching forces, are expected to return to initial structures after the release of an external force. The mutant R502W appears to represent lower stability; thus, the return to the initial structure may be of lower probability. The comparable analysis to other muscle domains (titin) and immunoglobulin domains suggests the very subtle relation to the biological activity of these proteins.

Toward a quantitative description of microscopic pathway heterogeneity in protein folding

Experimentally consistent statistical modeling of protein folding thermodynamics reveals unprecedented complexity with numerous parallel folding routes in five different proteins.

Reconciling Intermediates in Mechanical Unfolding Experiments with Two-State Protein Folding in Bulk

Most experimentally well-characterized single domain proteins of less than 100 residues have been found to be two-state folders. That is, only two distinct populations can explain both equilibrium and kinetic measurements. Results from single molecule force spectroscopy, where a protein is unfolded by applying a mechanical pulling force to its ends, have largely confirmed this description for proteins found to be two-state in ensemble experiments. Recently, however, stable intermediates have been reported in mechanical unfolding experiments on a cold-shock protein previously found to be a prototypical two-state folder. Here, we tackle this discrepancy using free energy landscapes and Markov state models derived from coarse-grained molecular simulations. We show that protein folding intermediates can be selectively stabilized by the pulling force and that the populations of these intermediates vary in a force-dependent manner. Our model qualitatively captures the experimental results and suggests a possible origin of the apparent discrepancy.

The Power of Force: Insights into the Protein Folding Process Using Single-Molecule Force Spectroscopy

One of the major challenges in modern biophysics is observing and understanding conformational changes during complex molecular processes, from the fundamental protein folding to the function of molecular machines. Single-molecule techniques have been one of the major driving forces of the huge progress attained in the last few years. Recent advances in resolution of the experimental setups, aided by theoretical developments and molecular dynamics simulations, have revealed a much higher degree of complexity inside these molecular processes than previously reported using traditional ensemble measurements. This review sums up the evolution of these developments and gives an outlook on prospective discoveries.