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.     

Unfolding a linker between helical repeats

In many multi-repeat proteins, linkers between repeats have little secondary structure and place few constraints on folding or unfolding. However, the large family of spectrin-like proteins, including alpha-actinin, spectrin, and dystrophin, share three-helix bundle, spectrin repeats that appear in crystal structures to be linked by long helices. All of these proteins are regularly subjected to mechanical stress. Recent single molecule atomic force microscopy (AFM) experiments demonstrate not only forced unfolding but also simultaneous unfolding of tandem repeats at finite frequency, which suggests that the contiguous helix between spectrin repeats can propagate a cooperative helix-to-coil transition. Here, we address what happens atomistically to the linker under stress by steered molecular dynamics simulations of tandem spectrin repeats in explicit water. The results for alpha-actinin repeats reveal rate-dependent pathways, with one pathway showing that the linker between repeats unfolds, which may explain the single-repeat unfolding pathway observed in AFM experiments. A second pathway preserves the structural integrity of the linker, which explains the tandem-repeat unfolding event. Unfolding of the linker begins with a splay distortion of proximal loops away from hydrophobic contacts with the linker. This is followed by linker destabilization and unwinding with increased hydration of the backbone. The end result is an unfolded helix that mechanically decouples tandem repeats. Molecularly detailed insights obtained here aid in understanding the mechanical coupling of domain stability in spectrin family proteins.     

Spectrin domains lose cooperativity in forced unfolding

Spectrin is a multidomain cytoskeletal protein, the component three-helix bundle domains are expected to experience mechanical force in vivo. In thermodynamic and kinetic studies, neighboring domains of chicken brain alpha-spectrin R16 and R17 have been shown to behave cooperatively. Is this cooperativity maintained under force? The effect of force on these spectrin domains was investigated using atomic force microscopy. The response of the individual domains to force was compared to that of the tandem repeat R1617. Importantly, nonhelical linkers (all-beta immunoglobulin domains) were used to avoid formation of nonnative helical linkers. We show that, in contrast to previous studies on spectrin repeats, only 3% of R1617 unfolding events gave an increase in contour length consistent with cooperative two-domain unfolding events. Furthermore, the unfolding forces for R1617 were the same as those for the unfolding of R16 or R17 alone. This is a strong indication that the cooperative unfolding behavior observed in the stopped-flow studies is absent between these spectrin domains when force is acting as a denaturant. Our evidence suggests that the rare double unfolding events result from misfolding between adjacent repeats. We suggest that this switch from cooperative to independent behavior allows multidomain proteins to maintain integrity under applied force.     

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.     

Mutation of conserved histidines alters tertiary structure and nanomechanics of consensus ankyrin repeats

The conserved TPLH tetrapeptide motif of ankyrin repeats (ARs) plays an important role in stabilizing AR proteins, and histidine (TPLH)-to-arginine (TPLR) mutations in this motif have been associated with a hereditary human anemia, spherocytosis. Here, we used a combination of atomic force microscopy-based single-molecule force spectroscopy and molecular dynamics simulations to examine the mechanical effects of His → Arg substitutions in TPLH motifs in a model AR protein, NI6C. Our molecular dynamics results show that the mutant protein is less mechanically stable than the WT protein. Our atomic force microscopy results indicate that the mechanical energy input necessary to fully unfold the mutant protein is only half of that necessary to unfold the WT protein (53 versus 106 kcal/mol). In addition, the ability of the mutant to generate refolding forces is also reduced. Moreover, the mutant protein subjected to cyclic stretch-relax measurements displays mechanical fatigue, which is absent in the WT protein. Taken together, these results indicate that the His → Arg substitutions in TPLH motifs compromise mechanical properties of ARs and suggest that the origin of hereditary spherocytosis may be related to mechanical failure of ARs.     

The Molecular Mechanism Underlying Mechanical Anisotropy of the Protein GB1

Mechanical responses of elastic proteins are crucial for their biological function and nanotechnological use. Loading direction has been identified as one key determinant for the mechanical responses of proteins. However, it is not clear how a change in pulling direction changes the mechanical unfolding mechanism of the protein. Here, we combine protein engineering, single-molecule force spectroscopy, and steered molecular dynamics simulations to systematically investigate the mechanical response of a small globular protein GB1. Force versus extension profiles from both experiments and simulations reveal marked mechanical anisotropy of GB1. Using native contact analysis, we relate the mechanically robust shearing geometry with concurrent rupture of native contacts. This clearly contrasts the sequential rupture observed in simulations for the mechanically labile peeling geometry. Moreover, we identify multiple distinct mechanical unfolding pathways in two loading directions. Implications of such diverse unfolding mechanisms are discussed. Our results may also provide some insights for designing elastomeric proteins with tailored mechanical properties.     

Force-Regulated Refolding of the Mechanosensory Domain in the Platelet Glycoprotein Ib-IX Complex

In platelets, the glycoprotein (GP) Ib-IX receptor complex senses blood shear flow and transmits the mechanical signals into platelets. Recently, we have discovered a juxtamembrane mechanosensory domain (MSD) within the GPIbα subunit of GPIb-IX. Mechanical unfolding of the MSD activates GPIb-IX signaling into platelets, leading to their activation and clearance. Using optical tweezer-based single-molecule force measurement, we herein report a systematic biomechanical characterization of the MSD in its native, full-length receptor complex and a recombinant, unglycosylated MSD in isolation. The native MSD unfolds at a resting rate of 9 × 10^?3 s^?1. Upon exposure to pulling forces, MSD unfolding accelerates exponentially over a force scale of 2.0 pN. Importantly, the unfolded MSD can refold with or without applied forces. The unstressed refolding rate of MSD is ～17 s^?1 and slows exponentially over a force scale of 3.7 pN. Our measurements confirm that the MSD is relatively unstable, with a folding free energy of 7.5 k_BT. Because MSD refolding may turn off GPIb-IX’s mechanosensory signals, our results provide a mechanism for the requirement of a continuous pulling force of >15 pN to fully activate GPIb-IX.     

Forced Unfolding of Apocytochrome <Emphasis Type="Italic">b</Emphasis><Subscript>5</Subscript> by Steered Molecular Dynamics Simulation

Apocytochrome b5 (apocyt b5), a small b-type cytochrome with heme prosthetic group removal, has been subjected to steered molecular dynamics (SMD) simulations for investigating the consequences of mechanical force-induced unfolding. Both constant velocity (0.5 and 1.0 A/ps) and constant force (500, 750 and 1000 pN) stretching have been employed to model forced unfolding of apocyt b5. The results of SMD simulations elucidate that apocyt b5 is protected against external stress mainly through the interstrand hydrogen bonding between its beta1-beta2 and beta2-beta3 strands, highlighting the importance of hydrophobic core 2 in stabilization of apocyt b5. The existence of intermediate states manifested by current simulations in the forced unfolding pathway of apocyt b5 is different from the observations in pervious thermal or chemical unfolding studies in the absence of force. The present study could thus provide insights into the relationship between the two cooperative functional modules of apocyt b5 and also guide the rational molecular design of heme proteins.     

Nanomechanics of Ig-like domains of human contactin (BIG-2)

Contactins are modular extracellular cell matrix proteins that are present in the brain, and they are responsible for the proper development and functioning of neurons. They contain six immunoglobulin-like IgC2 domains and four fibronectin type III repeats. The interactions of contactin with other proteins are poorly understood. The mechanical properties of all IgC2 domains of human contactin 4 were studied using a steered molecular dynamics approach and CHARMM force field with an explicit TIP3P water environment on a 10-ns timescale. Force spectra of all domains were determined computationally and the nanomechanical unfolding process is described. The domains show different mechanical stabilities. The calculated maxima of the unfolding force are in the range of 900–1700 pN at a loading rate of 7 N/s. Our data indicate that critical regions of IgC2 domains 2 and 3, which are responsible for interactions with tyrosine phosphatases and are important in nervous system development, are affected by even weak mechanical stretching. Thus, tensions present in the cell may modulate cellular activities related to contactin function. The present data should facilitate the interpretation of atomic force microscope single-molecule spectra of numerous proteins with similar IgC2 motives.     

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.     

States and transitions during forced unfolding of a single spectrin repeat

Spectrin is a vital and abundant protein of the cytoskeleton. It has an elongated structure that is made by a chain of so-called spectrin repeats. Each repeat contains three antiparallel alpha-helices that form a coiled-coil structure. Spectrin forms an oligomeric structure that is able to cross-link actin filaments. In red cells, the spectrin/actin meshwork underlying cell membrane is thought to be responsible for special elastic properties of the cell. In order to determine mechanical unfolding properties of the spectrin repeat, we have used single molecule force spectroscopy to study the states of unfolding of an engineered polymeric protein consisting of identical spectrin domains. We demonstrate that the unfolding of spectrin domains can occur in a stepwise fashion during stretching. The force-extension patterns exhibit features that are compatible with the existence of at least one intermediate between the folded and the completely unfolded conformation. Only those polypeptides that still contain multiple intact repeats display intermediates, indicating a stabilisation effect. Precise force spectroscopy measurements on single molecules using engineered protein constructs reveal states and transitions during the mechanical unfolding of spectrin. Single molecule force spectroscopy appears to open a new window for the analysis of transition probabilities between different conformational states.     

Influence of lateral association on forced unfolding of antiparallel spectrin heterodimers

Protein extensibility appears to be based broadly on conformational changes that can in principle be modulated by protein-protein interactions. Spectrin family proteins, with their extensible three-helix folds, enable evaluation of dimerization effects at the single molecule level by atomic force microscopy. Although some spectrin family members function physiologically only as homodimers (e.g. alpha-actinin) or are strictly monomers (e.g. dystrophin), alpha- and beta-spectrins are stable as monomeric forms but occur physiologically as alpha,beta-heterodimers bound laterally lengthwise. For short constructs of alpha- and beta-spectrin, either as monomers or as alpha,beta-dimers, sawtooth patterns in atomic force microscopy-forced extension show that unfolding stochastically extends repeats approximately 4-5-fold greater in length than native conformations. For both dimers and monomers, distributions of unfolding lengths appear bimodal; major unfolding peaks reflect single repeats, and minor unfolding peaks at twice the length reflect tandem repeats. Cooperative unfolding thus propagates through helical linkers between serial repeats (1, 2). With lateral heterodimers, however, the force distribution is broad and shifted to higher forces. The associated chains in a dimer can stay together and unfold simultaneously in addition to unfolding independently. Weak lateral interactions do not inhibit unfolding, but strong lateral interactions facilitate simultaneous unfolding analogous to serial repeat coupling within spectrin family proteins.     

Full Reconstruction of a Vectorial Protein Folding Pathway by Atomic Force Microscopy and Molecular Dynamics Simulations

During co-translational folding, the nascent polypeptide chain is extruded sequentially from the ribosome exit tunnel and, under severe conformational constraints, is dictated by its one-dimensional geometry. How do such vectorial constraints impact the folding pathway? Here, we combine single-molecule atomic force spectroscopy and steered molecular dynamics simulations to examine protein folding in the presence of one-dimensional constraints that are similar to those imposed on the nascent polypeptide chain. The simulations exquisitely reproduced the experimental unfolding and refolding force extension relationships and led to the full reconstruction of the vectorial folding pathway of a large polypeptide, the 253-residue consensus ankyrin repeat protein, NI6C. We show that fully stretched and then relaxed NI6C starts folding by the formation of local secondary structures, followed by the nucleation of three N-terminal repeats. This rate-limiting step is then followed by the vectorial and sequential folding of the remaining repeats. However, after partial unfolding, when allowed to refold, the C-terminal repeats successively regain structures without any nucleation step by using the intact N-terminal repeats as a template. These results suggest a pathway for the co-translational folding of repeat proteins and have implications for mechanotransduction.     

The mechanical hierarchies of fibronectin observed with single-molecule AFM

Mechanically induced conformational changes in proteins such as fibronectin are thought to regulate the assembly of the extracellular matrix and underlie its elasticity and extensibility. Fibronectin contains a region of tandem repeats of up to 15 type III domains that play critical roles in cell binding and self-assembly. Here, we use single-molecule force spectroscopy to examine the mechanical properties of fibronectin (FN) and its individual FNIII domains. We found that fibronectin is highly extensible due to the unfolding of its FNIII domains. We found that the native FNIII region displays strong mechanical unfolding hierarchies requiring 80 pN of force to unfold the weakest domain and 200 pN for the most stable domain. In an effort to determine the identity of the weakest/strongest domain, we engineered polyproteins composed of an individual domain and measured their mechanical stability by single-protein atomic force microscopy (AFM) techniques. In contrast to chemical and thermal measurements of stability, we found that the tenth FNIII domain is mechanically the weakest and that the first and second FNIII domains are the strongest. Moreover, we found that the first FNIII domain can acquire multiple, partially folded conformations, and that their incidence is modulated strongly by its neighbor FNIII domain. The mechanical hierarchies of fibronectin demonstrated here may be important for the activation of fibrillogenesis and matrix assembly.     

Mechanical Characterization of Protein L in the Low-Force Regime by Electromagnetic Tweezers/Evanescent Nanometry

Mechanical manipulation at the single molecule level of proteins exhibiting mechanical stability poses a technical challenge that has been almost exclusively approached by atomic force microscopy (AFM) techniques. However, due to mechanical drift limitations, AFM techniques are restricted to experimental recordings that last less than a minute in the high-force regime. Here we demonstrate a novel combination of electromagnetic tweezers and evanescent nanometry that readily captures the forced unfolding trajectories of protein L at pulling forces as low as 10 ∼ 15 pN. Using this approach, we monitor unfolding and refolding cycles of the same polyprotein for a period of time longer than 30 min. From such long-lasting recordings, we obtain ensemble averages of unfolding step sizes and rates that are consistent with single-molecule AFM data obtained at higher stretching forces. The unfolding kinetics of protein L at low stretching forces confirms and extends the observations that the mechanical unfolding rate is exponentially dependent on the pulling force within a wide range of stretching forces spanning from 13 pN up to 120 pN. Our experiments demonstrate a novel approach for the mechanical manipulation of single proteins for extended periods of time in the low-force regime.     

Mechanical Design of the Third FnIII Domain of Tenascin-C

By combining single-molecule atomic force microscopy (AFM), proline mutagenesis and steered molecular dynamics (SMD) simulations, we investigated the mechanical unfolding dynamics and mechanical design of the third fibronectin type III domain of tenascin-C (TNfn3) in detail. We found that the mechanical stability of TNfn3 is similar to that of other constituting FnIII domains of tenascin-C, and the unfolding process of TNfn3 is an apparent two-state process. By employing proline mutagenesis to block the formation of backbone hydrogen bonds and introduce structural disruption in β sheet, we revealed that in addition to the important roles played by hydrophobic core packing, backbone hydrogen bonds in β hairpins are also responsible for the overall mechanical stability of TNfn3. Furthermore, proline mutagenesis revealed that the mechanical design of TNfn3 is robust and the mechanical stability of TNfn3 is very resistant to structural disruptions caused by proline substitutions in β sheets. Proline mutant F88P is one exception, as the proline mutation at position 88 reduced the mechanical stability of TNfn3 significantly and led to unfolding forces of < 20 pN. This result suggests that Phe88 is a weak point of the mechanical resistance for TNfn3. We used SMD simulations to understand the molecular details underlying the mechanical unfolding of TNfn3. The comparison between the AFM results and SMD simulations revealed similarities and discrepancies between the two. We compared the mechanical unfolding and design of TNfn3 and its structural homologue, the tenth FnIII domain from fibronectin. These results revealed the complexity underlying the mechanical design of FnIII domains and will serve as a starting point for systematically analyzing the mechanical architecture of other FnIII domains in tenascins-C, and will help to gain a better understanding of some of the complex features observed for the stretching of native tenascin-C.     

Stabilization provided by neighboring strands is critical for the mechanical stability of proteins

Single-molecule force spectroscopy studies and steered molecular dynamics simulations have revealed that protein topology and pulling geometry play important roles in determining the mechanical stability of proteins. Most studies have focused on local interactions that are associated with the force-bearing β-strands. Interactions mediated by neighboring strands are often overlooked. Here we use Top7 and barstar as model systems to illustrate the critical importance of the stabilization effect provided by neighboring β-strands on the mechanical stability. Using single-molecule atomic force microscopy, we showed that Top7 and barstar, which have similar topology in their force-bearing region, exhibit vastly different mechanical-stability characteristics. Top7 is mechanically stable and unfolds at ∼150 pN, whereas barstar is mechanically labile and unfolds largely below 50 pN. Steered molecular dynamics simulations revealed that stretching force peels one force-bearing strand away from barstar to trigger unfolding, whereas Top7 unfolds via a substructure-sliding mechanism. This previously overlooked stabilization effect from neighboring β-strands is likely to be a general mechanism in protein mechanics and can serve as a guideline for the de novo design of proteins with significant mechanical stability and novel protein topology.     

Single molecule force spectroscopy reveals a weakly populated microstate of the FnIII domains of tenascin

The native states of proteins exist as an ensemble of conformationally similar microstates. The fluctuations among different microstates are of great importance for the functions and structural stability of proteins. Here, we demonstrate that single molecule atomic force microscopy (AFM) can be used to directly probe the existence of multiple folded microstates. We used the AFM to repeatedly stretch and relax a recombinant tenascin fragment TNfnALL to allow the fibronectin type III (FnIII) domains to undergo repeated unfolding/refolding cycles. In addition to the native state, we discovered that some FnIII domains can refold from the unfolded state into a previously unrecognized microstate, N* state. This novel state is conformationally similar to the native state, but mechanically less stable. The native state unfolds at approximately 120 pN, while the N* state unfolds at approximately 50 pN. These two distinct populations of microstates constitute the ensemble of the folded states for some FnIII domains. An unfolded FnIII domain can fold into either one of the two microstates via two distinct folding routes. These results reveal the dynamic and heterogeneous picture of the folded ensemble for some FnIII domains of tenascin, which may carry important implications for the mechanical functions of tenascins in vivo.     

Determination of structural models of the complex between the cytoplasmic domain of erythrocyte band 3 and ankyrin-R repeats 13-24

The adaptor protein ankyrin-R interacts via its membrane binding domain with the cytoplasmic domain of the anion exchange protein (AE1) and via its spectrin binding domain with the spectrin-based membrane skeleton in human erythrocytes. This set of interactions provides a bridge between the lipid bilayer and the membrane skeleton, thereby stabilizing the membrane. Crystal structures for the dimeric cytoplasmic domain of AE1 (cdb3) and for a 12-ankyrin repeat segment (repeats 13-24) from the membrane binding domain of ankyrin-R (AnkD34) have been reported. However, structural data on how these proteins assemble to form a stable complex have not been reported. In the current studies, site-directed spin labeling, in combination with electron paramagnetic resonance (EPR) and double electron-electron resonance, has been utilized to map the binding interfaces of the two proteins in the complex and to obtain inter-protein distance constraints. These data have been utilized to construct a family of structural models that are consistent with the full range of experimental data. These models indicate that an extensive area on the peripheral domain of cdb3 binds to ankyrin repeats 18-20 on the top loop surface of AnkD34 primarily through hydrophobic interactions. This is a previously uncharacterized surface for binding of cdb3 to AnkD34. Because a second dimer of cdb3 is known to bind to ankyrin repeats 7-12 of the membrane binding domain of ankyrin-R, the current models have significant implications regarding the structural nature of a tetrameric form of AE1 that is hypothesized to be involved in binding to full-length ankyrin-R in the erythrocyte membrane.     

Examining the influence of linkers and tertiary structure in the forced unfolding of multiple-repeat spectrin molecules

The unfolding pathways of multiple-repeat spectrin molecules were examined using steered molecular dynamics (SMD) simulations to forcibly unfold double- and triple-repeat spectrin molecules. Although SMD has previously been used to study other repeating-domain proteins, spectrin offers a unique challenge in that the linker connecting repeat units has a definite secondary structure, that of an alpha-helix. Therefore, the boundary conditions imposed on a double- or triple-repeat spectrin must be carefully considered if any relationship to the real system is to be deduced. This was accomplished by imposing additional forces on the system which ensure that the terminal alpha-helices behave as if there were no free noncontiguous helical ends. The results of the SMD simulations highlight the importance of the rupture of the alpha-helical linker on the subsequent unfolding events. Rupture of the linker propagates unfolding in the adjacent repeat units by destabilizing the tertiary structure, ultimately resulting in complete unfolding of the affected repeat unit. Two dominant classes of unfolding pathways are observed after the initial rupture of a linker which involve either rupture of another linker (possibly adjacent) or rupture of the basic tertiary structure of a repeat unit. The relationship between the force response observed on simulation timescales and those of experiment or physiological conditions is also discussed.