The effect of core destabilization on the mechanical resistance of I27

It is still unclear whether mechanical unfolding probes the same pathways as chemical denaturation. To address this point, we have constructed a concatamer of five mutant I27 domains (denoted (I27)(5)*) and used it for mechanical unfolding studies. This protein consists of four copies of the mutant C47S, C63S I27 and a single copy of C63S I27. These mutations severely destabilize I27 (DeltaDeltaG(UN) = 8.7 and 17.9 kJ mol(-1) for C63S I27 and C47S, C63S I27, respectively). Both mutations maintain the hydrogen bond network between the A' and G strands postulated to be the major region of mechanical resistance for I27. Measuring the speed dependence of the force required to unfold (I27)(5)* in triplicate using the atomic force microscope allowed a reliable assessment of the intrinsic unfolding rate constant of the protein to be obtained (2.0 x 10(-3) s(-1)). The rate constant of unfolding measured by chemical denaturation is over fivefold faster (1.1 x 10(-2) s(-1)), suggesting that these techniques probe different unfolding pathways. Also, by comparing the parameters obtained from the mechanical unfolding of a wild-type I27 concatamer with that of (I27)(5)*, we show that although the observed forces are considerably lower, core destabilization has little effect on determining the mechanical sensitivity of this domain.     

Steered molecular dynamics studies of titin I1 domain unfolding

The cardiac muscle protein titin, responsible for developing passive elasticity and extensibility of muscle, possesses about 40 immunoglobulin-like (Ig) domains in its I-band region. Atomic force microscopy (AFM) and steered molecular dynamics (SMD) have been successfully combined to investigate the reversible unfolding of individual Ig domains. However, previous SMD studies of titin I-band modules have been restricted to I27, the only structurally known Ig domain from the distal region of the titin I-band. In this paper we report SMD simulations unfolding I1, the first structurally available Ig domain from the proximal region of the titin I-band. The simulations are carried out with a view toward upcoming atomic force microscopy experiments. Both constant velocity and constant force stretching have been employed to model mechanical unfolding of oxidized I1, which has a disulfide bond bridging beta-strands C and E, as well as reduced I1, in which the disulfide bridge is absent. The simulations reveal that I1 is protected against external stress mainly through six interstrand hydrogen bonds between its A and B beta-strands. The disulfide bond enhances the mechanical stability of oxidized I1 domains by restricting the rupture of backbone hydrogen bonds between the A'- and G-strands. The disulfide bond also limits the maximum extension of I1 to approximately 220 A. Comparison of the unfolding pathways of I1 and I27 are provided and implications to AFM experiments are discussed.     

Computational investigation of the effect of thermal perturbation on the mechanical unfolding of titin I27

The emergence of single-molecule force measurement experiments has facilitated a better understanding of protein folding pathways and the thermodynamics involved. Computational methods such as steered molecular dynamics (SMD) simulations are helpful in providing atomistic level information on the unfolding pathways. Recent experimental studies have showed that combinations of single-molecule experiments with traditional methods such as chemical and/or thermal denaturation yield additional insights into the folding phenomenon. In this study, we report results from extensive computations (a total of about 60 SMD simulations with a total length of about 0.4 μs) that address the effect of thermal perturbation on the mechanical stability of the I27 domain of the protein titin. A wide range of temperatures (280-340 K) were considered for the pulling, which was done at both constant velocity and constant force using SMD simulations. Good agreement with experimental data, such as for the trends in changes in average force and the maximum force with respect to the temperature, was obtained. This study identifies two competing pathways for the mechanical unfolding of I27, and illustrates the significance of combining various techniques to examine protein folding.     

Dwell-time distribution analysis of polyprotein unfolding using force-clamp spectroscopy

Using the recently developed single molecule force-clamp technique we quantitatively measure the kinetics of conformational changes of polyprotein molecules at a constant force. In response to an applied force of 110 pN, we measure the dwell times of 1647 unfolding events of individual ubiquitin modules within each protein chain. We then establish a rigorous method for analyzing force-clamp data using order statistics. This allows us to test the success of a history-independent, two-state model in describing the kinetics of the unfolding process. We find that the average unfolding trajectory is independent of the number of protein modules N in each trajectory, which varies between 3 and 12 (the engineered protein length), suggesting that the unfolding events in each chain are uncorrelated. We then derive a binomial distribution of dwell times to describe the stochastic dynamics of protein unfolding. This distribution successfully describes 81% of the data with a single rate constant of alpha = 0.6 s(-1) for all N. The remainder of the data that cannot be accounted for suggests alternative unfolding barriers in the energy landscape of the protein. This method investigates the statistical features of unfolding beyond the average measurement of a single rate constant, thus providing an attractive alternative for measuring kinetics by force-clamp spectroscopy.     

Frequency modulation atomic force microscopy reveals individual intermediates associated with each unfolded I27 titin domain

In this study, we apply a dynamic atomic force microscopy (AFM) technique, frequency modulation (FM) detection, to the mechanical unfolding of single titin I27 domains and make comparisons with measurements made using the AFM contact or static mode method. Static mode measurements revealed the well-known force transition occurring at 100-120 pN in the first unfolding peak, which was less clear, or more often absent, in the subsequent unfolding peaks. In contrast, some FM-AFM curves clearly resolved a force transition associated with each of the unfolding peaks irrespective of the number of observed unfolded domains. As expected for FM-AFM, the frequency shift response of the main unfolding peaks and their intermediates could only be detected when the oscillation amplitudes used were smaller than the interaction lengths being measured. It was also shown that the forces measured for the dynamical interaction of the FM-AFM technique were significantly lower than those measured using the static mode. This study highlights the potential for using dynamic AFM for investigating biological interactions, including protein unfolding and the detection of novel unfolding intermediates.     

Probing the mechanical folding kinetics of TAR RNA by hopping, force-jump, and force-ramp methods

Mechanical unfolding and refolding of single RNA molecules have previously been observed in optical traps as sudden changes in molecular extension. Two methods have been traditionally used: "force-ramp", with the applied force continuously changing, and "hopping". In hopping experiments the force is held constant and the molecule jumps spontaneously between two different states. Unfolding/refolding rates are measured directly, but only over a very narrow range of forces. We have now developed a force-jump method to measure the unfolding and refolding rates independently over a wider range of forces. In this method, the applied force is rapidly stepped to a new value and either the unfolding or refolding event is monitored through changes in the molecular extension. The force-jump technique is compared to the force-ramp and hopping methods by using a 52-nucleotide RNA hairpin with a three-nucleotide bulge, i.e., the transactivation response region RNA from the human immunodeficiency virus. We find the unfolding kinetics and Gibbs free energies obtained from all three methods to be in good agreement. The transactivation response region RNA hairpin unfolds in an all-or-none two-state reaction at any loading rate with the force-ramp method. The unfolding reaction is reversible at small loading rates, but shows hysteresis at higher loading rates. Although the RNA unfolds and refolds without detectable intermediates in constant-force conditions (hopping and force-jump), it shows partially folded intermediates in force-ramp experiments at higher unloading rates. Thus, we find that folding of RNA hairpins can be more complex than a simple single-step reaction, and that application of several methods can improve understanding of reaction mechanisms.     

Pathway shifts and thermal softening in temperature-coupled forced unfolding of spectrin domains

Pathways of unfolding a protein depend in principle on the perturbation-whether it is temperature, denaturant, or even forced extension. Widely-shared, helical-bundle spectrin repeats are known to melt at temperatures as low as 40-45 degrees C and are also known to unfold via multiple pathways as single molecules in atomic force microscopy. Given the varied roles of spectrin family proteins in cell deformability, we sought to determine the coupled effects of temperature on forced unfolding. Bimodal distributions of unfolding intervals are seen at all temperatures for the four-repeat beta(1-4) spectrin-an alpha-actinin homolog. The major unfolding length corresponds to unfolding of a single repeat, and a minor peak at twice the length corresponds to tandem repeats. Increasing temperature shows fewer tandem events but has no effect on unfolding intervals. As T approaches T(m), however, mean unfolding forces in atomic force microscopy also decrease; and circular dichroism studies demonstrate a nearly proportional decrease of helical content in solution. The results imply a thermal softening of a helical linker between repeats which otherwise propagates a helix-to-coil transition to adjacent repeats. In sum, structural changes with temperature correlate with both single-molecule unfolding forces and shifts in unfolding pathways.     

The key event in force-induced unfolding of Titin's immunoglobulin domains

Steered molecular dynamics simulation of force-induced titin immunoglobulin domain I27 unfolding led to the discovery of a significant potential energy barrier at an extension of approximately 14 A on the unfolding pathway that protects the domain against stretching. Previous simulations showed that this barrier is due to the concurrent breaking of six interstrand hydrogen bonds (H-bonds) between beta-strands A' and G that is preceded by the breaking of two to three hydrogen bonds between strands A and B, the latter leading to an unfolding intermediate. The simulation results are supported by Angstrom-resolution atomic force microscopy data. Here we perform a structural and energetic analysis of the H-bonds breaking. It is confirmed that H-bonds between strands A and B break rapidly. However, the breaking of the H-bond between strands A' and G needs to be assisted by fluctuations of water molecules. In nanosecond simulations, water molecules are found to repeatedly interact with the protein backbone atoms, weakening individual interstrand H-bonds until all six A'-G H-bonds break simultaneously under the influence of external stretching forces. Only when those bonds are broken can the generic unfolding take place, which involves hydrophobic interactions of the protein core and exerts weaker resistance against stretching than the key event.     

Force-clamp spectroscopy of single-protein monomers reveals the individual unfolding and folding pathways of I27 and ubiquitin

Single-protein force experiments have relied on a molecular fingerprint based on tethering multiple single-protein domains in a polyprotein chain. However, correlations between these domains remain an issue in interpreting force spectroscopy data, particularly during protein folding. Here we first show that force-clamp spectroscopy is a sensitive technique that provides a molecular fingerprint based on the unfolding step size of four single-monomer proteins. We then measure the force-dependent unfolding rate kinetics of ubiquitin and I27 monomers and find a good agreement with the data obtained for the respective polyproteins over a wide range of forces, in support of the Markovian hypothesis. Moreover, with a large statistical ensemble at a single force, we show that ubiquitin monomers also exhibit a broad distribution of unfolding times as a signature of disorder in the folded protein landscape. Furthermore, we readily capture the folding trajectories of monomers that exhibit the same stages in folding observed for polyproteins, thus eliminating the possibility of entropic masking by other unfolded modules in the chain or domain-domain interactions. On average, the time to reach the I27 folded length increases with increasing quenching force at a rate similar to that of the polyproteins. Force-clamp spectroscopy at the single-monomer level reproduces the kinetics of unfolding and refolding measured using polyproteins, which proves that there is no mechanical effect of tethering proteins to one another in the case of ubiquitin and I27.     

Molecular force modulation spectroscopy revealing the dynamic response of single bacteriorhodopsins

Recent advances in atomic force microscopy allowed globular and membrane proteins to be mechanically unfolded on a single-molecule level. Presented is an extension to the existing force spectroscopy experiments. While unfolding single bacteriorhodopsins from native purple membranes, small oscillation amplitudes (6-9 nm) were supplied to the vertical displacement of the cantilever at a frequency of 3 kHz. The phase and amplitude response of the cantilever-protein system was converted to reveal the elastic (conservative) and viscous (dissipative) contributions to the unfolding process. The elastic response (stiffness) of the extended parts of the protein were in the range of a few tens pN/nm and could be well described by the derivative of the wormlike chain model. Discrete events in the viscous response coincided with the unfolding of single secondary structure elements and were in the range of 1 microNs/m. In addition, these force modulation spectroscopy experiments revealed novel mechanical unfolding intermediates of bacteriorhodopsin. We found that kinks result in a loss of unfolding cooperativity in transmembrane helices. Reconstructing force-distance spectra by the integration of amplitude-distance spectra verified their position, offering a novel approach to detect intermediates during the forced unfolding of single proteins.     

Monte Carlo simulation for single RNA unfolding by force

Using polymer elastic theory and known RNA free energies, we construct a Monte Carlo algorithm to simulate the single RNA folding and unfolding by mechanical force on the secondary structure level. For the constant force ensemble, we simulate the force-extension curves of the P5ab, P5abc deltaA, and P5abc molecules in equilibrium. For the constant extension ensemble, we focus on the mechanical behaviors of the RNA P5ab molecule, which include the unfolding force dependence on the pulling speed, the force-hysteresis phenomenon, and the coincidence of stretching-relaxing force-curves in thermal equilibrium. We particularly simulate the time traces of the end-to-end distance of the P5ab under the constant force in equilibrium, which also have been recorded in the recent experiment. The reaction rate constants for the folding and unfolding are calculated. Our results show that the agreement between the simulation and the experimental measurements is satisfactory.     

Mechanical unfolding of TNfn3: the unfolding pathway of a fnIII domain probed by protein engineering, AFM and MD simulation

Protein engineering Phi-value analysis combined with single molecule atomic force microscopy (AFM) was used to probe the molecular basis for the mechanical stability of TNfn3, the third fibronectin type III domain from human tenascin. This approach has been adopted previously to solve the forced unfolding pathway of a titin immunoglobulin domain, TI I27. TNfn3 and TI I27 are members of different protein superfamilies and have no sequence identity but they have the same beta-sandwich structure consisting of two antiparallel beta-sheets. TNfn3, however, unfolds at significantly lower forces than TI I27. We compare the response of these proteins to mechanical force. Mutational analysis shows that, as is the case with TI I27, TNfn3 unfolds via a force-stabilised intermediate. The key event in forced unfolding in TI I27 is largely the breaking of hydrogen bonds and hydrophobic interactions between the A' and G-strands. The mechanical Phi-value analysis and molecular dynamics simulations reported here reveal that significantly more of the TNfn3 molecule contributes to its resistance to force. Both AFM experimental data and molecular dynamics simulations suggest that the rate-limiting step of TNfn3 forced unfolding reflects a transition from the extended early intermediate to an aligned intermediate state. As well as losses of interactions of the A and G-strands and associated loops there are rearrangements throughout the core. As was the case for TI I27, the forced unfolding pathway of TNfn3 is different from that observed in denaturant studies in the absence of force.     

A Theoretical Model for the Mechanical Unfolding of Repeat Proteins

We consider the mechanical stretching of a polypeptide chain formed by multiple interacting repeats. The folding thermodynamics and the interactions among the repeats are described by the Ising model. Unfolded repeats act as soft entropic springs, whereas folded repeats respond to a force as stiffer springs. We show that the resulting force-extension curve may exhibit a pronounced force maximum corresponding to the unfolding of the first repeat. This event is followed by the unfolding of the remaining repeats, which takes place at a lower force. As the protein extension is increased, the force-extension curve of a sufficiently long repeat protein displays a plateau, where the force remains nearly constant and the protein unfolds sequentially so that the number of unfolded repeats is proportional to the extension. Such a sequential mechanical unfolding mechanism is displayed even by the repeat proteins whose thermal denaturation is highly cooperative, provided that they are long enough. By contrast, the unfolding of short repeat progressions can be cooperative.     

40 Probing single molecule RNA folding using temperature and force

Nucleic acids can be unfolded either by temperature, such as in UV melting, or by mechanical force using optical tweezers. In UV melting experiments, the folding free energy of nucleic acids at mesophilic temperatures are extrapolated from unfolding occurring at elevated temperatures. Additionally, single molecule unfolding experiments are typically performed only at room temperature, preventing calculation of changes in enthalpy and entropy. Here, we present temperature-controlled optical tweezers suitable for studying folding of single RNA molecules at physiological temperatures. Constant temperatures between 22 and 37?°C are maintained with an accuracy of 0.1?°C, whereas the optical tweezers display a spatial resolution of ～1?nm over the temperature range. Using this instrument, we measured the folding thermodynamics and kinetics of a 20-base-pair RNA hairpin by force-ramp and constant force experiments. Between 22 and 37?°C, the hairpin unfolds and refolds in a single step. Increasing temperature decreases the stability of the hairpin and thus decreases the force required to unfold it. The equilibrium force, at which unfolding and refolding rates are equal, drops ～1?pN as temperature increases every 5?°C. At each temperature, the folding energy can be quantified by reversible work done to unfold the RNA and from the equilibrium constant at constant forces. Over the experimental temperature range, the folding free energy of the hairpin depends linearly on temperature, indicating that ΔH is constant. The measured folding thermodynamics are further compared with the nearest neighbor calculations using Turner’s parameters of nucleic acid folding energetics.     

Dynamic measurement of single protein's mechanical properties

Dimerized (tandemly repeated) protein was constructed, and the stretching force during the unfolding of the single protein molecule was measured using an atomic force microscope. In quasistatic measurements using normal force-distance curve measurements, each monomer unit was unfolded step by step. To elucidate the conformational state at each extension length, we measured the relax-stress response of the protein using short stroke sinusoidal movements of the sample stage. This allowed us to investigate the dynamic response of the protein repeatedly without full stretching or rupturing. Although the protein molecule responded in-phase to the applied movement in most cases, we found a novel out-of-phase response around the stretching length where the second monomer unit unfolded. Applying the spring constant measured in the quasistatic experiment, the out-of-phase response was reproduced in the simple calculation, which suggested the folding and the unfolding at the second monomer unit were taking place repeatedly during the relax-stress response measurement.     

Changing the Mechanical Unfolding Pathway of FnIII10 by Tuning the Pulling Strength

We investigate the mechanical unfolding of the tenth type III domain from fibronectin (FnIII₁₀) both at constant force and at constant pulling velocity, by all-atom Monte Carlo simulations. We observe both apparent two-state unfolding and several unfolding pathways involving one of three major, mutually exclusive intermediate states. All three major intermediates lack two of seven native β-strands, and share a quite similar extension. The unfolding behavior is found to depend strongly on the pulling conditions. In particular, we observe large variations in the relative frequencies of occurrence for the intermediates. At low constant force or low constant velocity, all three major intermediates occur with a significant frequency. At high constant force or high constant velocity, one of them, with the N- and C-terminal β-strands detached, dominates over the other two. Using the extended Jarzynski equality, we also estimate the equilibrium free-energy landscape, calculated as a function of chain extension. The application of a constant pulling force leads to a free-energy profile with three major local minima. Two of these correspond to the native and fully unfolded states, respectively, whereas the third one can be associated with the major unfolding intermediates.     

Can non-mechanical proteins withstand force? Stretching barnase by atomic force microscopy and molecular dynamics simulation

Atomic force microscopy (AFM) experiments have provided intriguing insights into the mechanical unfolding of proteins such as titin I27 from muscle, but will the same be possible for proteins that are not physiologically required to resist force? We report the results of AFM experiments on the forced unfolding of barnase in a chimeric construct with I27. Both modules are independently folded and stable in this construct and have the same thermodynamic and kinetic properties as the isolated proteins. I27 can be identified in the AFM traces based on its previous characterization, and distinct, irregular low-force peaks are observed for barnase. Molecular dynamics simulations of barnase unfolding also show that it unfolds at lower forces than proteins with mechanical function. The unfolding pathway involves the unraveling of the protein from the termini, with much more native-like secondary and tertiary structure being retained in the transition state than is observed in simulations of thermal unfolding or experimentally, using chemical denaturant. Our results suggest that proteins that are not selected for tensile strength may not resist force in the same way as those that are, and that proteins with similar unfolding rates in solution need not have comparable unfolding properties under force.     

Comparing proteins by their unfolding pattern

Single molecule force spectroscopy has evolved into an important and extremely powerful technique for investigating the folding potentials of biomolecules. Mechanical tension is applied to individual molecules, and the subsequent, often stepwise unfolding is recorded in force extension traces. However, because the energy barriers of the folding potentials are often close to the thermal energy, both the extensions and the forces at which these barriers are overcome are subject to marked fluctuations. Therefore, force extension traces are an inadequate representation despite widespread use particularly when large populations of proteins need to be compared and analyzed. We show in this article that contour length, which is independent of fluctuations and alterable experimental parameters, is a more appropriate variable than extension. By transforming force extension traces into contour length space, histograms are obtained that directly represent the energy barriers. In contrast to force extension traces, such barrier position histograms can be averaged to investigate details of the unfolding potential. The cross-superposition of barrier position histograms allows us to detect and visualize the order of unfolding events. We show with this approach that in contrast to the sequential unfolding of bacteriorhodopsin, two main steps in the unfolding of the enzyme titin kinase are independent of each other. The potential of this new method for accurate and automated analysis of force spectroscopy data and for novel automated screening techniques is shown with bacteriorhodopsin and with protein constructs containing GFP and titin kinase.