Unfolding mechanics of holo- and apocalmodulin studied by the atomic force microscope

The structural stability of calmodulin (CaM) has been investigated previously by chemical and thermal methods. The calcium-loaded form of CaM has been found to be exceptionally stable, because it can be exposed to temperatures >90 degrees C or to a 9 M urea solution without a marked change in its tertiary structure, and is therefore not experimentally accessible for unfolding studies using conventional analytical methods. In this study, we have developed a system for measuring the force for mechanically unfolding CaM using an atomic force microscope (AFM) by stretching the protein from its N- and C-terminal residues; we have been successful in obtaining force versus extension (F-E) curves for both apo and holo forms of CaM. In our experiment, distinguishable F-E curves have been obtained upon stretching of apoCaM and holoCaM to their full extensions. A very low force observed upon stretching of apoCaM indicated a relatively high flexibility of the apo form. On the contrary, a relatively high unfolding force and the appearance of a characteristic force peak were noted during full stretching of holoCaM. The F-E curve of the latter form of CaM most likely reflects a more rigid and probably more organized conformation of holoCaM than that of apoCaM. These experiments confirmed that the AFM is able to clearly distinguish two functionally distinct forms of CaM in terms of their mechanical properties.     

Extending a spectrin repeat unit. II: rupture behavior

A spectrin repeat unit was subject to extension using cyclic expansion nonequilibrium molecular dynamics. Periodic boundary conditions were used to examine the effects of the contiguous alpha-helical linker on the force response. The measured force-extension curve shows a linear increase in the force response when the spectrin repeat unit is extended by approximately 0.4 nm. After that point, the force response peaks and subsequently declines. The peak in the force response marks the point where the spectrin repeat unit undergoes a change in its material properties from a strongly elastic material to a mostly viscous one, on the timescales of the simulations. The force peak is also correlated with rupture of the alpha-helical linker, and is likely the event responsible for the peaks in the sawtooth-pattern force-extension curves measured by atomic force microscopy experiments. Rupture of the linker involves simultaneously breaking approximately four hydrogen bonds that maintain the alpha-helical linker. After this initial rupture, the linker undergoes simple helix-to-coil transitions as the spectrin repeat unit continues to be extended. The implications of linker rupture in the interpretation of unfolding and atomic force microscopy experiments are also discussed.     

Unfolding of titin immunoglobulin domains by steered molecular dynamics simulation.

Titin, a 1-microm-long protein found in striated muscle myofibrils, possesses unique elastic and extensibility properties in its I-band region, which is largely composed of a PEVK region (70% proline, glutamic acid, valine, and lysine residue) and seven-strand beta-sandwich immunoglobulin-like (Ig) domains. The behavior of titin as a multistage entropic spring has been shown in atomic force microscope and optical tweezer experiments to partially depend on the reversible unfolding of individual Ig domains. We performed steered molecular dynamics simulations to stretch single titin Ig domains in solution with pulling speeds of 0.5 and 1.0 A/ps. Resulting force-extension profiles exhibit a single dominant peak for each Ig domain unfolding, consistent with the experimentally observed sequential, as opposed to concerted, unfolding of Ig domains under external stretching forces. This force peak can be attributed to an initial burst of backbone hydrogen bonds, which takes place between antiparallel beta-strands A and B and between parallel beta-strands A' and G. Additional features of the simulations, including the position of the force peak and relative unfolding resistance of different Ig domains, can be related to experimental observations.     

Nanotechnology and Protein Mechanics

The atomic force microscope is currently used in our and many other laboratories to measure the mechanical response of polypeptide and proteins against tensile forces applied to well defined positions in their chemical structures. The resulting force vs. extension (F-E) curves are analyzed in relation to their known conformations under various conditions. The method can be extended to study the mechanical responses of other, often much larger biological structures, and extract the component proteins and DNAs from cell membranes and chromosomes.     

Thermal unfolding simulations of apo-calmodulin using leap-dynamics

The simulation method leap-dynamics (LD) has been applied to protein thermal unfolding simulations to investigate domain-specific unfolding behavior. Thermal unfolding simulations of the 148-residue protein apo-calmodulin with implicit solvent were performed at temperatures 290 K, 325 K, and 360 K and compared with the corresponding molecular dynamics trajectories in terms of a number of calculated conformational parameters. The main experimental results of unfolding are reproduced in showing the lower stability of the C-domain: at 290 K, both the N- and C-domains are essentially stable; at 325 K, the C-domain unfolds, whereas the N-domain remains folded; and at 360 K, both domains unfold extensively. This behavior could not be reproduced by molecular dynamics simulations alone under the same conditions. These results show an encouraging degree of convergence between experiment and LD simulation. The simulations are able to describe the overall plasticity of the apo-calmodulin structure and to reveal details such as reversible folding/unfolding events within single helices. The results show that by using the combined application of a fast and efficient sampling routine with a detailed molecular dynamics force field, unfolding simulations of proteins at atomic resolution are within the scope of current computational power.     

Force-induced unfolding of human telomeric G-quadruplex: A steered molecular dynamics simulation study

We study the unfolding of a parallel G-quadruplex from human telomeric DNA by mechanical stretching using steered molecular dynamics (MD) simulation. We find that the force curves and unfolding processes strongly depend on the pulling sites. With pulling sites located on the sugar-phosphate backbone, the force-extension curve shows a single peak and the unfolding proceeds sequentially. Pulling sites located on the terminal nucleobases lead to a force-extension curve with two peaks and the unfolding is more cooperative. Simulations of the refolding of partially unfolded quadruplexes show very different behavior for the two different pulling modalities. In particular, starting from an unfolded state prepared by nucleobase pulling leads to a long-lived intermediate state whose existence is also corroborated by the free energy profile computed with the Jarzynski equation. Based on this observation, we propose a novel folding pathway for parallel G-quadruplexes with the human telomere sequence.     

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.     

Cooperativity in forced unfolding of tandem spectrin repeats

Force-driven conformational changes provide a broad basis for protein extensibility, and multidomain proteins broaden the possibilities further by allowing for a multiplicity of forcibly extended states. Red cell spectrin is prototypical in being an extensible, multidomain protein widely recognized for its contribution to erythrocyte flexibility. Atomic force microscopy has already shown that single repeats of various spectrin family proteins can be forced to unfold reversibly under extension. Recent structural data indicates, however, that the linker between triple-helical spectrin repeats is often a contiguous helix, thus raising questions as to what the linker contributes and what defines a domain mechanically. We have examined the extensible unfolding of red cell spectrins as monomeric constructs of just two, three, or four repeats from the actin-binding ends of both alpha- and beta-chains, i.e., alpha(18-21) and beta(1-4) or their subfragments. In addition to single repeat unfolding evident in sawtooth patterns peaked at relatively low forces (<50 pN at 1 nm/ms extension rates), tandem repeat unfolding is also demonstrated in ensemble-scale analyses of thousands of atomic force microscopy contacts. Evidence for extending two chains and loops is provided by force versus length scatterplots which also indicate that tandem repeat unfolding occurs at a significant frequency relative to single repeat unfolding. Cooperativity in forced unfolding of spectrin is also clearly demonstrated by a common force scale for the unfolding of both single and tandem repeats.     

A kinetic molecular model of the reversible unfolding and refolding of titin under force extension.

Molecular elasticity is a physicomechanical property that is associated with a select number of polypeptides and proteins, such as the giant muscle protein, titin, and the extracellular matrix protein, tenascin. Both proteins have been the subject of atomic force microscopy (AFM), laser tweezer, and other in vitro methods for examining the effects of force extension on the globular (FNIII/Ig-like) domains that comprise each protein. In this report we present a time-dependent method for simulating AFM force extension and its effect on FNIII/Ig domain unfolding and refolding. This method treats the unfolding and refolding process as a standard three-state protein folding model (U right arrow over left arrow T right arrow over left arrow F, where U is the unfolded state, T is the transition or intermediate state, and F is the fully folded state), and integrates this approach within the wormlike chain (WLC) concept. We simulated the effect of AFM tip extension on a hypothetical titin molecule comprised of 30 globular domains (Ig or FNIII) and 25% Pro-Glu-Val-Lys (PEVK) content, and analyzed the unfolding and refolding processes as a function of AFM tip extension, extension rate, and variation in PEVK content. In general, we find that the use of a three-state protein-folding kinetic-based model and the implicit inclusion of PEVK domains can accurately reproduce the experimental force-extension curves observed for both titin and tenascin proteins. Furthermore, our simulation data indicate that PEVK domains exhibit extensibility behavior, assist in the unfolding and refolding of FNIII/Ig domains in the titin molecule, and act as a force "buffer" for the FNIII/Ig domains, particularly at low and moderate extension forces.     

Remarkable disparity in mechanical response among the extracellular domains of type I and II cadherins

Cadherins, a large family of calcium-dependent adhesion molecules, are critical for intercellular adhesion. While crystallographic structures for several cadherins show clear structural similarities, their relevant adhesive strengths vary and their mechanisms of adhesion between types I and II cadherin subfamilies are still unclear. Here, stretching of cadherins was explored experimentally by atomic force microscopy and computationally by steered molecular dynamics (SMD) simulations, where partial unfolding of the E-cadherin ectodomains was observed. The SMD simulations on strand-swapping cadherin dimers displayed similarity in binding strength, suggesting contributions of other mechanisms to explain the strength differences of cell adhesion in vivo. Systematic simulations on the unfolding of the extracellular domains of type I and II cadherins revealed diverse pathways. However, at the earliest stage, a remarkable similarity in unfolding was observed for the various type I cadherins that was distinct from that for type II cadherins. This likely correlates positively with their distinct adhesive properties, suggesting that the initial forced deformation in type I cadherins may be involved in cadherin-mediated adhesion.

An animated Interactive 3D Complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:JBSD:25     

Mechanical anisotropy of ankyrin repeats

Red blood cells are frequently deformed and their cytoskeletal proteins such as spectrin and ankyrin-R are repeatedly subjected to mechanical forces. While the mechanics of spectrin was thoroughly investigated in vitro and in vivo, little is known about the mechanical behavior of ankyrin-R. In this study, we combine coarse-grained steered molecular dynamics simulations and atomic force spectroscopy to examine the mechanical response of ankyrin repeats (ARs) in a model synthetic AR protein NI6C, and in the D34 fragment of native ankyrin-R when these proteins are subjected to various stretching geometry conditions. Our steered molecular dynamics results, supported by AFM measurements, reveal an unusual mechanical anisotropy of ARs: their mechanical stability is greater when their unfolding is forced to propagate from the N-terminus toward the C-terminus (repeats unfold at ~60 pN), as compared to the unfolding in the opposite direction (unfolding force ～ 30 pN). This anisotropy is also reflected in the complex refolding behavior of ARs. The origin of this unfolding and refolding anisotropy is in the various numbers of native contacts that are broken and formed at the interfaces between neighboring repeats depending on the unfolding/refolding propagation directions. Finally, we discuss how these complex mechanical properties of ARs in D34 may affect its behavior in vivo.     

Mechanical unfolding of single filamin A (ABP-280) molecules detected by atomic force microscopy

Filamin A (ABP-280), which is an actin-binding protein of 560 kDa as a dimer, can, together with actin filaments, produce an isotropic cross-linked three-dimensional network (actin/filamin A gel) that plays an important role in mechanical responses of cells in processes such as maintenance of membrane stability and translational locomotion. In this study, we investigated the mechanical properties of single filamin A molecules using atomic force microscopy. In force-extension curves, we observed sawtooth patterns corresponding to the unfolding of individual immunoglobulin (Ig)-fold domains of filamin A. At a pulling speed of 0.37 microm/s, the unfolding interval was sharply distributed around 30 nm, while the unfolding force ranged from 50 to 220 pN. This wide distribution of the unfolding force can be explained by variation in values of activation energy and the width of activation barrier of 24 Ig-fold domains of the filamin A at the unfolding transition. This unfolding can endow filamin A with great extensibility. The refolding of the unfolded chain of filamin A occurred when the force applied to the protein was reduced to near zero, indicating that its unfolding is reversible. Based on these results, we discuss here the physiological implications of the mechanical properties of single filamin A molecules.     

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.     

Pathways and intermediates in forced unfolding of spectrin repeats

Spectrin repeats are triple-helical coiled-coil domains found in many proteins that are regularly subjected to mechanical stress. We used atomic force microscopy technique and steered molecular dynamics simulations to study the behavior of a wild-type spectrin repeat and two mutants. The experiments indicate that spectrin repeats can form stable unfolding intermediates when subjected to external forces. In the simulations the unfolding proceeded via a variety of pathways. Stable intermediates were associated to kinking of the central helix close to a proline residue. A mutant stabilizing the central helix showed no intermediates in experiments, in agreement with simulation. Spectrin repeats may thus function as elastic elements, extendable to intermediate states at various lengths.     

Mechanical unfolding of a titin Ig domain: structure of transition state revealed by combining atomic force microscopy,protein engineering and molecular dynamics simulations

Titin I27 shows a high resistance to unfolding when subject to external force. To investigate the molecular basis of this mechanical stability, protein engineering Phi-value analysis has been combined with atomic force microscopy to investigate the structure of the barrier to forced unfolding. The results indicate that the transition state for forced unfolding is significantly structured, since highly destabilising mutations in the core do not affect the force required to unfold the protein. As has been shown before, mechanical strength lies in the region of the A' and G-strands but, contrary to previous suggestions, the results indicate clearly that side-chain interactions play a significant role in maintaining mechanical stability. Since Phi-values calculated from molecular dynamics simulations are the same as those determined experimentally, we can, with confidence, use the molecular dynamics simulations to analyse the structure of the transition state in detail, and are able to show loss of interactions between the A' and G-strands with associated A-B and E-F loops in the transition state. The key event is not a simple case of loss of hydrogen bonding interactions between the A' and G-strands alone. Comparison with Phi-values from traditional folding studies shows differences between the force and "no-force" transition states but, nevertheless, the region important for kinetic stability is the same in both cases. This explains the correspondence between hierarchy of kinetic stability (measured in stopped-flow denaturant studies) and mechanical strength in these titin domains.     

Detecting molecular interactions that stabilize native bovine rhodopsin

Using single-molecule force spectroscopy we probed molecular interactions within native bovine rhodopsin and discovered structural segments of well-defined mechanical stability. Highly conserved residues among G protein-coupled receptors were located at the interior of individual structural segments, suggesting a dual role for these segments in rhodopsin. Firstly, structural segments stabilize secondary structure elements of the native protein, and secondly, they position and hold the highly conserved residues at functionally important environments. Two main classes of force curves were observed. One class corresponded to the unfolding of rhodopsin with the highly conserved Cys110-Cys187 disulfide bond remaining intact and the other class corresponded to the unfolding of the entire rhodopsin polypeptide chain. In the absence of the Cys110-Cys187 bond, the nature of certain molecular interactions within folded rhodopsin was altered. These changes highlight the structural importance of this disulfide bond and may form the basis of dysfunctions associated with its absence.     

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.     

Chemistry on a single protein, vascular cell adhesion molecule-1, during forced unfolding

Proteins of many types experience tensile forces in their normal function, and vascular cell adhesion molecule-1 (VCAM-1) is typical in this. VCAM has seven Ig domains, and each has a disulfide bond (-S-S-) buried in its core that covalently stabilizes about half of each domain against unfolding. VCAM is extended here by single molecule atomic force microscopy in the presence or absence of reducing agents. In the absence of reducing agent, a sawtooth pattern of forced unfolding reveals an average period and total length consistent with disulfide locations in VCAM. With increasing reducing agent, accessible disulfides are specifically reduced (to SH); the average period for unfolding increases up to saturation together with additional metrics of unfolding. Steered molecular dynamics simulations of unfolding indeed show that the core disulfide bond is solvent-exposed in the very earliest stages of protein extension. Michaelis-Menten kinetics emerge with reduction catalyzed by force (tau(reduction) approximately 10(-4) s). The results establish single molecule reduction, one bond at a time, and show that mechanical forces can play a key role in modulating the redox state of cell adhesion proteins that are invariably stressed in cell adhesion.     

Biophysical studies on the RNA cores of satellite tobacco mosaic virus

Satellite tobacco mosaic virus (STMV) was probed using a variety of proteases. Consequences of the degradation were analyzed using gel electrophoresis, quasi-elastic light scattering (QELS), and atomic force microscopy (AFM). Proteolysis rates of 30 minutes for complete degradation of the protein capsid, up to many hours, were investigated. With each protease, degradation of virions 17 nm in diameter was shown by QELS to result in particles of 10 nm diameter, which is that of the RNA core observed in the virion by x-ray diffraction analysis. This was verified by direct visualization with atomic force microscopy. Using QELS, it was further shown that freshly prepared RNA cores remain as individual, stable, 10-nm condensed particles for 12 to 24 h. Clusters of particles then formed, followed by very large aggregates of 500 to 1000 nm diameter. AFM showed that the aggregates were composed of groups of the condensed RNA cores and were not due to unfolding of the nucleic acid. No unfolding of the core particles into extended conformation was seen by AFM until the samples were heated well beyond 90 degrees C. Mass spectrometry of RNA core particles revealed the presence of a major polypeptide whose amino acid sequence corresponded to residues 2 through 25 of the coat protein. Amino acids 13 through 25 were previously observed to be in direct contact with the RNA and are presumably protected from protease digestion. Low resolution difference Fourier analyses indicated the courses of the remainders of the amino terminal strands (amino acids 2-12) in intact virions. Any individual strand appears to have several choices of path, which accounts for the observed disorder at high resolution. These positively charged strands, serving as virtual polyamines, engage the helical segments of RNA. The intimate association of amino acid residues 2 through 25 with RNA likely contributes to the stability of the condensed conformation of the nucleic acid cores.     

Exact low-force kinetics from high-force single-molecule unfolding events

Mechanical forces play a key role in crucial cellular processes involving force-bearing biomolecules, as well as in novel single-molecule pulling experiments. We present an exact method that enables one to extrapolate, to low (or zero) forces, entire time-correlation functions and kinetic rate constants from the conformational dynamics either simulated numerically or measured experimentally at a single, relatively higher, external force. The method has twofold relevance: 1), to extrapolate the kinetics at physiological force conditions from molecular dynamics trajectories generated at higher forces that accelerate conformational transitions; and 2), to extrapolate unfolding rates from experimental force-extension single-molecule curves. The theoretical formalism, based on stochastic path integral weights of Langevin trajectories, is presented for the constant-force, constant loading rate, and constant-velocity modes of the pulling experiments. For the first relevance, applications are described for simulating the conformational isomerization of alanine dipeptide; and for the second relevance, the single-molecule pulling of RNA is considered. The ability to assign a weight to each trace in the single-molecule data also suggests a means to quantitatively compare unfolding pathways under different conditions.