Cooperative unfolding of distinctive mechanoreceptor domains transduces force into signals
Abstract
How cells sense their mechanical environment and transduce forces into biochemical signals is a crucial yet unresolved question in mechanobiology. Platelets use receptor glycoprotein Ib (GPIb), specifically its α subunit (GPIbα), to signal as they tether and translocate on von Willebrand factor (VWF) of injured arterial surfaces against blood flow. Force elicits catch bonds to slow VWF–GPIbα dissociation and unfolds the GPIbα leucine-rich repeat domain (LRRD) and juxtamembrane mechanosensitive domain (MSD). How these mechanical processes trigger biochemical signals remains unknown. Here we analyze these extracellular events and the resulting intracellular Ca2+ on a single platelet in real time, revealing that LRRD unfolding intensifies Ca2+ signal whereas MSD unfolding affects the type of Ca2+ signal. Therefore, LRRD and MSD are analog and digital force transducers, respectively. The >30 nm macroglycopeptide separating the two domains transmits force on the VWF–GPIbα bond (whose lifetime is prolonged by LRRD unfolding) to the MSD to enhance its unfolding, resulting in unfolding cooperativity at an optimal force. These elements may provide design principles for a generic mechanosensory protein machine.
https://doi.org/10.7554/eLife.15447.001eLife digest
Platelets – the blood clotting cells – have the ability to detect, interpret and respond to mechanical forces, such as those generated by the flow of blood. The magnitude and duration of the forces detected by the platelets influences whether they form a blood clot. Understanding how the platelets respond to mechanical forces is therefore crucial for our knowledge of conditions such as thrombosis, where blood clots form inside vessels and block them. Clots that form within arteries are associated with heart attack and stroke, which account for around one third of all deaths worldwide.
Cells can sense external forces via individual proteins on their surface and transmit the mechanical information across the cell membrane. This triggers signals within the cell that influence how it responds. However, the molecular details of these “mechanosensory” processes remain poorly understood.
To patch up damaged blood vessels, platelets use a protein on their surface named glycoprotein Ibα (GPIbα) to bind to a plasma protein called von Willebrand factor that adheres to the vessel wall. This binding tethers the platelet to the blood vessel and activates it during clot formation. Previous studies suggested that mechanical force affects how this binding triggers the signals that activate platelets.
Ju, Chen et al. used a homebuilt nanotool to pull on platelet GPIbα while it was bound to von Willebrand factor. This revealed that two distinct domains of the GPIbα protein unfold to relay information about the force, such as its magnitude and duration, to the platelet to trigger biochemical signalling inside the cell. The unfolding of each GPIbα domain has a distinct role in determining the quantity and quality of the signals. The unfolding events work synergistically – they occur together to produce an effect that’s greater than the sum of their individual effects. However, pulling on GPIbα via a mutant form of von Willebrand factor eliminated the synergy between the two unfolding events, therefore hindering the effective conversion of mechanical forces into biochemical signals.
Notably, the two GPIbα domains unfolded by force exist in many protein families, including those involved in mediating cell adhesion and detecting signals. The biophysical tools developed by Ju, Chen et al. could be extended to analyze how mechanical cues are presented, received, transmitted and converted into biochemical signals in other cell types and biological systems. Furthermore, the structural insights gained from the platelet GPIbα system may help to design a generic mechanosensory protein machine.
https://doi.org/10.7554/eLife.15447.002Introduction
Platelets can serve as a natural model system for studying cell mechanosensing as they rapidly respond to changes in hydrodynamic forces and substrate stiffness due to vascular pathology (Jackson, 2011; Qiu et al., 2015). Previous studies have suggested the role of GPIbα as a mechanoreceptor, for force exerted on it via its ligand VWF induces platelet signaling (Ruggeri, 2015). Conceptually, this coupled mechanical-biochemical process (mechanosensing) can be broken down into four steps: 1) Mechanopresentation: the receptor binding domain A1 is exposed by structural changes in VWF induced by elongational flow and collagen immobilization (Ju et al., 2015a; Springer, 2014); 2) Mechanoreception: GPIbα LRRD receives the force signal via engaging VWF-A1 to tether the platelet against shear stress; 3) Mechanotransmission: force is propagated from the LRRD through the mucin-like macroglycopeptide (MP) stalk (cf. Figure 2A) (Fox et al., 1988) and the MSD across the membrane to adaptor and signaling molecules (e.g. 14-3-3ζ) inside the platelet (cf. Figure 7G); and 4) Mechanotransduction: force induces mechano-chemical changes to convert mechanical cues to biochemical signals. Some of these steps have been characterized separately. For example, GPIbα forms catch-slip bonds with wild-type (WT) A1 in >15 pN, such that the bond lifetime first increases with force, reaches a maximum at ~25 pN, and decreases thereafter; whereas it forms slip-only bonds with type 2B von Willebrand disease (VWD) mutant (e.g. A1R1450E), such that the bond lifetime decreases monotonically with force (Ju et al., 2013; Yago et al., 2008). As another example, force induces unfolding of the LRRD, which prolongs A1–GPIbα bond lifetime (Ju et al., 2015b), and of the MSD, which is hypothesized to play a role in platelet signaling (Zhang et al., 2015). However, how these inter-connected steps are orchestrated to enable the information encoded by force to be translated into biochemical signals is still poorly understood.
We used a biomembrane force probe (BFP) to recapitulate the above process in a single-cell and single-molecular bond level to address the following questions: 1) What molecular events would be induced in GPIbα and how these events are regulated mechanically? 2) Whether, and if so, how changes in presentation of force by VWF-A1 mutation would affect the force reception by GPIbα and its response to force? 3) What features of the force (waveforms) could be sensed by the platelet via GPIbα to initiate intraplatelet calcium fluxes? 4) What proximal events may be responsible for transducing force into a biochemical signal? By manipulating the mechanopresentation and mechanoreception steps then analyzing the resulting mechanotransmission and mechanotransduction steps, we gained insights into the inner workings of this GPIbα-mediated mechanosensory machine.
Results
In the BFP setup, a probe bead was functionalized with VWF-A1 or an anti-GPIb monoclonal antibody (mAb) to serve as a surrogate subendothelial surface (Figure 1A,B). It was attached to the apex of a micropipette-aspirated red blood cell (RBC) to form an ultrasensitive force transducer (Liu et al., 2014). A platelet was aspired by the target pipette to contact the bead in repetitive force-ramp or force-clamp cycles to mimic the sequential formation, force loading, and dissociation of VWF–GPIbα bonds during the translocation of a platelet on the sub-endothelium (Video 1; Figure 1C and Figure 1—figure supplement 1; Materials and methods). Adhesion frequencies from these cycles were kept low (<20%) by adjusting the ligand or antibody density, a condition required for the platelet to be pulled predominantly (>89%) by a single GPIbα bond (Chesla et al., 1998; Zhu et al., 2002). Control experiments using beads lacking ligand showed no binding, and blocking with mAb AK2 (epitope mapped to leucine-rich repeat 1–2 overlapping the A1 binding site, cf. Figure 2A) eliminated GPIbα binding to A1 but not to mAb WM23 (epitope mapped to the MP below LRRD, cf. Figure 2A) (Figure 1D) (Dong et al., 2001). This confirmed binding specificity and that the binding site of A1 is within LRRD but the binding site of WM23 is outside (Zhang et al., 2015).
Identification of LRRD and MSD unfolding
Using an optical trap, Zhang et al. observed force-induced MSD unfolding in purified recombinant full-length GPIb-IX and a GPIbα stalk region construct (Zhang et al., 2015). Using a BFP, we observed LRRD unfolding in glycocalicin (GC) (Ju et al., 2015b), the extracellular segment of GPIbα lacking the MSD (Liang et al., 2013) (Figure 2A–C). Here we pulled GPIbα on platelets via A1 and observed two unfolding signatures, one in the ramping and the other in the clamping phases of the force trace (Figure 1E). Unfolding that occurred in the ramping phase is termed ramped unfolding, which is featured by a sudden force kink at 5–20 pN as observed in previous studies of GPIbα unfolding (Ju et al., 2015b; Zhang et al., 2015). Similar to findings of protein unfolding studies (Kellermayer et al., 1997; Rief et al., 1997; Tskhovrebova et al., 1997; Zhang et al., 2009a, 2015), both the force-extension curves before and after unfolding were well fitted by the worm-like chain (WLC) model (Figure 1—figure supplement 2). Unfolding that occurred in the clamping phase is termed clamped unfolding, which is featured by an abrupt force drop (Figure 1F). Although not observed in the previous studies of GPIbα unfolding (Ju et al., 2015b; Zhang et al., 2015), this feature has been described in protein unfolding studies using force-clamp experiments (Oberhauser et al., 2001; Tskhovrebova et al., 1997).
Unfolding lengths derived from both signatures were measured from the probe bead position vs. time data (Figure 1E insert, 1F and Figure 1—figure supplement 2). The lengths of individual ramped unfolding events distributed tri-modally with three subpopulations (Figure 2D and Figure 2—figure supplement 1; Materials and methods). The first subpopulation coincides with the ramped unfolding length distribution from WM23 vs. platelet experiments (Figure 2E, white bars). WM23 binds the MP region below the LRRD (Figure 2A), hence could unfold MSD only. The average unfolding force vs. length data from the WM23 experiment was well fitted by the WLC model, yielding a contour length of 25.99 ± 0.85 nm (Figure 2H) that matches the previously reported MSD contour length (Zhang et al., 2015). The average unfolding force vs. length data from the A1 experiment overlaid well on the same WLC model fit (Figure 2H). These results identify the first subpopulation in Figure 2D as MSD unfolding.
The second subpopulation in Figure 2D matches the histogram of ramped unfolding lengths of GC pulled via A1 (Figure 2E, blue bars) that ranges from 18–56 nm and peaks at 36 nm (length of leucine-rich repeats 3–6). The average unfolding force vs. length plots derived from the A1 vs. platelet and A1 vs. GC experiments overlaid well on the same WLC model fit (Figure 2I). The best-fit contour length (70.29 ± 3.56 nm) matches the length of LRRD, calculated using a 4-Å contour length per residue (Ju et al., 2015b). These results identify the second subpopulation in Figure 2D as LRRD unfolding.
The third subpopulation in Figure 2D can be identified as concurrent unfolding of both MSD and LRRD that occurred within too short a time elapse to be distinguished by our BFP as two separate events, because its maximum unfolding length (85 nm) matches the sum of the observed maximum MSD and LRRD unfolding lengths.
Similar tri-modally distributed ramped unfolding lengths were obtained by using mAb AN51 (epitope mapped to the N-terminal flanking region above LRRD, cf. Figure 2A) instead of A1 to pull the platelet GPIbα (Figure 2F), and the second subpopulation also matches the ramped unfolding length distribution obtained using AN51 to pull GC (Figure 2G, blue bars). These results are expected because the unfolding lengths are determined by the respective primary structures of the LRRD and MSD, and as such should not depend on the 'grabbing handle' used to pull GPIbα. The consistence of the A1 and AN51 results imparts confidence in our identification of the three subpopulations as unfolding of MSD, LRRD, and both, respectively.
Interestingly, the two force waveforms induced unfolding of different GPIbα domains. Clamped forces unfolded only MSD as the lengths of clamped unfolding distribute as a single peak at 20 nm (Figure 2E,G, red bars), matching the first subpopulation in Figure 2D,F, respectively, regardless of whether platelet GPIbα was engaged by A1 or AN51. Furthermore, unfolding of LRRD in GC was induced only by ramped forces but not clamped forces (Figure 2J,K). By comparison, pulling platelet GPIbα via WM23 with both ramped and clamped forces induced MSD unfolding events with similar occurrence frequencies and unfolding lengths (Figure 2J,K). These results indicate that MSD can be unfolded by increasing forces as well as constant forces. By comparison, LRRD unfolding requires increasing forces. Some force-clamp cycles (Figure 1C; Video 2) generated two consecutive unfolding events, one in the ramping and the other in the clamping phase (Figure 1E). The respective unfolding lengths of the ramped and clamped unfolding events were 34–55 nm and 13–25 nm that totaled 47–80 nm, agreeing with those of the LRRD, MSD, and MSD+LRRD subpopulations in Figure 2D,F. Together, these results provide criteria to determine whether and which GPIbα domain(s) is unfolded (Figure 3—source data 1A).
Force- and ligand-dependent cooperativity between LRRD and MSD unfolding
To characterize the mechanical response of GPIbα, we measured the frequency, force and length of LRRD and MSD unfolding induced by a range of clamped forces exerted on platelet GPIbα or GC by A1WT or a type 2B VWD mutant A1R1450E. The ramped unfolding frequencies of both domains were extremely low at ≤10 pN but increased with the higher levels of clamped forces (Figure 3A,B). Interestingly, LRRD, but not MSD, unfolded more frequently when platelet GPIbα (Figure 3A) and GC (Figure 3—figure supplement 1) were pulled by A1WT than A1R1450E. The ramped unfolding forces of both domains increased with the clamped force and were indifferent to whether force was applied via WT or R1450E mutant of A1 (Figure 3C,D). In general, a higher force was required to unfold LRRD than MSD. Surprisingly, pulling platelet GPIbα via different ligands generated distinctive MSD clamped unfolding frequency vs. force plots: increasing initially and decreasing after reaching maximal at 25 pN when pulled by A1WT, but decreasing monotonically when pulled by A1R1450E (Figure 3E). These data suggest that the mechanoreceptor GPIbα may be able to interpret mechanical cues and discriminate ligands by responding to different force waveforms applied via different ligands with distinct LRRD and MSD unfolding frequencies. In addition, the distinctive force-dependences of two subpopulations of events that we deemed as respective LRRD and MSD unfolding provide further support for our criteria for their identification and classification.
-
Figure 3—source data 1
- https://doi.org/10.7554/eLife.15447.011
The spatial separation of LRRD and MSD by the >30 nm long MP stalk and the distinctive dependences of their unfolding on the force waveform would seem to favor these two GPIbα domains to unfold independently. This hypothesis predicts that the probability for LRRD and MSD to unfold concurrently should be equal to the product of the respective probabilities for LRRD and MSD to unfold separately. To test this hypothesis, we estimated these probabilities from the observed unfolding occurrence frequencies. At 25 pN, the 34.5% of BFP force-clamp cycles with unfolding events consist of 7.6, 17.2, 6.9, and 2.8% of unfolding of LRRD alone, MSD alone, LRRD and MSD sequentially, and concurrently (Figure 3—source data 1A). Significantly, the frequency of observing both LRRD and MSD unfolding in the same binding cycle, calculated by pooling together both cases of two domains unfolding sequentially and concurrently, P(MSD+LRRD), is much higher than the product of their respective occurrence frequencies, P(MSD)×P(LRRD), which is the joint probability for both to unfold assuming that they were independent (Figure 3—source data 1B).
These data suggest that the two GPIbα domains may unfold cooperatively, i.e., one domain unfolding may increase the likelihood for the other to unfold. To quantify the degree of such cooperativity, we defined a relative probability difference, ∆P/P = [P(MSD+LRRD) - P(MSD)×P(LRRD)]/[P(MSD)×P(LRRD)]. ∆P/P > 0 indicates positive cooperativity between LRRD and MSD unfolding. No cooperativity was observed at 10 pN because this force was insufficient to induce appreciable LRRD unfolding. Pulling with A1WT by a 25 pN clamped force generated high cooperativity, and further increase in force decreased cooperativity (Figure 3F). Remarkably, unfolding cooperativity was completely abolished at all forces when applied via the VWD mutant A1R1450E (Figure 3F).
We used χ2 test to determine if the hypothesis that MSD and LRRD unfolded independently should be rejected (Materials and methods). At 25 pN, LRRD unfolding significantly enhanced MSD unfolding (p = 3.09 × 10–4). The χ2 test results are depicted as negative log p-values vs. force plots in Figure 3G,H for A1WT and A1R1450E, respectively. Interestingly, significant (p = 0.05, dashed horizontal lines) unfolding cooperativity was observed only for A1WT at 25 and 40 pN. These data show that the cooperativity between LRRD and MSD unfolding is force- and ligand-dependent.
Model for cooperativity between LRRD and MSD unfolding
To elucidate the mechanism underlying the force- and ligand-dependent unfolding cooperativity, we note that when the MSD unfolding events were separately analyzed according to their occurrence in the ramping or clamping phase, MSD clamped, but not ramped, unfolding was significantly (p= 8.79 × 10–3 vs. 0.076 at 25 pN) enhanced by LRRD unfolding (Figure 3G), which occurred in the ramping phase only. This dominance of cooperativity by sequential rather than concurrent unfolding suggests a model for LRRD unfolding to impact MSD unfolding, which includes three ideas. The first idea has to do with the MSD time-to-unfold, tu (cf. Figure 1F). Our force-clamp measurements revealed similar tu values induced by A1WT or A1R1450E pulling (Figure 4A). The only exception is at 10 pN where a shorter tu was induced by A1WT than A1R1450E. This can be explained by their differential bond lifetimes (Figure 4B,C). Compared to A1R1450E, the much shorter lifetime of GPIbα bond with A1WT at 10 pN may underestimate tu because early dissociation of GPIbα would prevent observation of slow MSD unfolding events. This reasoning provides the second idea for our model: MSD clamped unfolding should occur before A1–GPIbα dissociation. The third idea comes from our previous observation (Ju et al., 2015b) that LRRD unfolding significantly prolongs GPIbα bond lifetime with A1WT (Figure 4B) but not A1R1450E (Figure 4C). Combining these three ideas, our model proposes that the A1–GPIbα bond lifetime, regulated by force and prolonged by LRRD unfolding in respective ligand-specific manners, determines the occurrence of MSD clamped unfolding, which, despite its ligand-independent unfolding kinetics, generates a cooperativity pattern that maximizes at the optimal force of 25 pN for A1WT but not for A1R1450E.
-
Figure 4—source data 1
- https://doi.org/10.7554/eLife.15447.014
To formulate the model mathematically, we multiplied the respective probability densities of the exponentially distributed MSD time-to-unfold (tu) (Figure 4—figure supplement 1A) and the dual-exponentially distributed lifetime (tb) of GPIbα bonds with A1WT or A1R1450E (Ju et al., 2013) to construct a joint probability density surface over the tu-tb plane (Figure 4D and Figure 4—figure supplement 1B). The predicted MSD clamped unfolding probability is the volume under this surface over the region 0<tu<tb<∞ (Materials and methods). When the model was tested against experiment, not only did the calculated force-dependent MSD unfolding frequency match the biphasic pattern for A1WT (Figure 4E) and the monophasic pattern for A1R1450E (Figure 4F), but it also compared well with the observed occurrence frequencies numerically at all forces. Remarkably, the model predicts both the quantitative enhancement of MSD unfolding by LRRD unfolding for A1WT and the lack of enhancement for A1R1450E without a single freely adjustable fitting parameter. The excellent agreement between theory and experiment has provided strong support for our model and explained the data in Figure 4E,F.
Platelet signaling induced by mechanoreception via a single GPIbα
Platelet translocation on VWF signals through GPIbα to induce Ca2+ fluxes (Mazzucato et al., 2002; Nesbitt et al., 2002). We optimized the fluorescence BFP (fBFP) method (Chen et al., 2015; Liu et al., 2014) for single-platelet calcium imaging and studied how platelet signaling was triggered by GPIbα mechanoreception via a sequence of intermittent single bonds under a range of clamped forces. The Ca2+ signals over the 200-s observation window of repeated platelet contact cycles were classified into three types (Figure 5A,B): i) null-type, featured by a basal trace with a maximum Ca2+ intensity increase (normalized by its initial value) ΔImax<0.05; ii) α-type, featured by an initial latent phase followed by a spike (mostly ΔImax>0.5) with a quick decay (Video 2); iii) β-type, featured by fluctuating signals around the baseline or gradually increasing signals to an intermediate level (mostly ΔImax<0.5) followed by a gradual decay to baseline (Video 3). The null type reflects the baseline with background noise, while the α- and β-types match the previous characterization of platelet internal Ca2+release triggered by VWF–GPIbα bonds measured in flow chamber experiments (Mazzucato et al., 2002). For each platelet, the calcium trace was overlaid with the sequential binding events, bond lifetimes, and their accumulation over the repeated platelet binding cycles (Figure 5B and Figure 5—figure supplement 1A,B).
Pulled by a 25 pN clamp force, A1WT–GPIbα bonds triggered much higher ΔImax than controls (Figure 5C), showing 28, 42, and 30% of null-, α- and β-types, respectively (Figure 5D). Control experiments that merely held aspirated platelets or contacted them by beads without coating any ligand showed null-type Ca2+ only (Figure 5C,D). The α-type Ca2+ could also be triggered by pulling GPIbα with AN51 but not with an anti-GPIbβ mAb (Figure 5C,D and Figure 5—figure supplement 1A,B), despite that GPIbβ is tightly connected to GPIbα within one GPIb complex and has been postulated to play a role in signaling through GPIb (Strassel et al., 2006). These data demonstrated the necessity of GPIbα engagement to trigger intraplatelet Ca2+ and agree with the previous report that α-type Ca2+ peaks occur when platelets are transiently arrested in the whole blood flow (Mazzucato et al., 2002;Nesbitt et al., 2002)
The concurrent measurements of A1–GPIbα binding kinetics and intraplatelet Ca2+ allowed us to determine the pre-Ca2+ bond lifetimes (Figure 5B), enabling single platelet correlative analysis of binding and signaling. Using the normalized maximum calcium intensity ΔImax to represent the Ca2+ level, we compared its correlations with three statistics of A1WT–GPIbα bond lifetimes occurred prior to calcium onset. We found that ΔImax correlates best with the pre-Ca2+ longest bond lifetime tmax (Figure 5E), similarly well with the pre-Ca2+ cumulative lifetime ∑ti (Figure 5—figure supplement 1C), but poorly with the pre-Ca2+ average lifetime <t> (Figure 5—figure supplement 1D). Careful examination of many overlaid calcium and bond lifetime traces, as exemplified in Figure 5B, revealed that the ∑ti values are generally dominated by the tmax values, which are usually much longer than the rest of the pre-Ca2+ bond lifetimes and are immediately followed by the calcium onset before observing additional shorter bond lifetimes. In other words, for each platelet usually ∑ti could be approximated by tmax but <t> is of a smaller and variable value. This observation explains why calcium correlates equally well with tmax and ∑ti but not with <t>. Importantly, these results also suggest that a single long-lived GPIbα bond is sufficient to trigger Ca2+ in a platelet. This assertion has been further supported by the parallel analysis of the data for A1R1450E. Although for R1450E the tmax value was significantly shorter and the α-type Ca2+ population was greatly reduced, the ΔImax still showed similar correlation with tmax (Figure 5F). Thus the pre-Ca2+ longest bond lifetime correlates the Ca2+ strength and type.
GPIbα discriminates ligands and shows different force-dependent calcium responses
We next asked whether the mechanoreceptor GPIbα is capable of sensing differences in the force waveform and discriminating the ligand through which force is applied. We first performed force-ramp experiment to generate a wide range of rupture forces using three ramping rates: 100, 1000 and 10,000 pN/s. However, only low levels of β-type Ca2+ were resulted (Figure 6A,B), showing no correlation with the largest rupture force prior to calcium onset (Figure 6C,D, right ordinate), regardless of whether platelets were tested by A1WT (Figure 6A,C) or A1R1450E (Figure 6B,D). In sharp contrast, much higher levels of Ca2+ of α- and β-types were induced by clamped forces applied to GPIbα via either A1WT (Figure 6E,G) or A1R1450E (Figure 6F,H) despite their much lower levels than the rupture forces seen in the force-ramp experiments (Figure 6C,D). Concurrently, the longest of A1–GPIbα bond lifetimes that occurred prior to Ca2+ onset was measured on each platelet and averaged over all platelets in each group. This pre-Ca2+ longest bond lifetime, tmax, exhibited catch-slip bond behavior for A1WT and slip-only bond behavior for A1R1450E (Figure 6G,H, right ordinate), just as the corresponding average bond lifetimes previously measured regardless of the intraplatelet calcium (Ju et al., 2013; Yago et al., 2008). Remarkably, the force-dependent pattern of calcium signals matched that of the pre-Ca2+ longest bond lifetimes for both A1WT and A1R1450E. The ligand-independent positive correlation of Ca2+ signal with tmax is consistent with a previously observed inverse correlation between the cytosolic Ca2+ level and the translocation velocity of platelets on immobilized VWF (Nesbitt et al., 2002). This is expected because the platelet translocation velocity is an inverse metric of VWF–GPIbα bond lifetime (Ju et al., 2013; Yago et al., 2008).
Interplay among GPIbα engagement duration, domain unfolding and signal initiation
The findings that durable force is important to both MSD unfolding and Ca2+ triggering prompted us to investigate the relationship between GPIbα domain unfolding and platelet signal initiation. We segregated the Ca2+ data generated by a 25 pN clamped force on A1WT–GPIbα bonds according to whether or not and, if so, which domain(s) was (were) unfolded prior to calcium onset. Platelets whose tests contained no unfolding event showed short tmax and low calcium of β- and null-types (Figure 7A). Platelets whose tests contained at least one pre-Ca2+ MSD unfolding event but no LRRD unfolding showed slightly longer tmax and higher Ca2+ of mostly α-type. By comparison, only β-type Ca2+ was observed in the rare (2.6%) cases where LRRD unfolded but MSD did not. Since in these cases the tmax values were much longer, this data excludes tmax to be the direct determining parameter for the Ca2+ type. Remarkably, the group with pre-Ca2+ unfolding of both LRRD and MSD exhibited long tmax and high Ca2+ of mostly α-type.
Consistent results were obtained using WM23 to pull GPIbα to bypass LRRD unfolding (cf. Figure 2A), showing significantly longer tmax and higher calcium for platelets with than without a MSD unfolding event and a clear α- vs. β-type signal distinction between them (Figure 7B). The higher α-type Ca2+ triggering efficiency of WM23 than A1 (compare MSD+ columns in Figure 7A and B) may be explained, at least in part, by the slower kinetics of GPIbα dissociation from WM23 than A1, which generated 70% more bond lifetime events by contacting a platelet for 200 s with WM23 than A1 (Figure 7—figure supplement 1A). The bonds were also 74% more durable (Figure 7—figure supplement 1B), resulting in a slightly higher (although not significant) MSD clamped unfolding probability per lifetime event for WM23 (Figure 7—figure supplement 1C). The expected number of MSD clamped unfolding over a 200-s experimental period, calculated as the product of the number of lifetime events and the unfolding probability per lifetime event, was significantly higher for WM23 than A1 (Figure 7—figure supplement 1D). The probability of platelet to have at least one MSD unfolding, calculated as 1– (1– unfolding probability per lifetime event) ^ (# lifetime events), is 40% and 60% for A1 and WM23, respectively (Figure 7—figure supplement 1E), close to the experimental results (51% vs. 69%, Figure 7A,B).
Interestingly, despite their high levels, ramped forces generated very few MSD unfolding events and triggered only null/β- but not α-type Ca2+ regardless of whether A1WT or WM23 was used to pull (Figure 7C,D). Together, these data indicate that both force-induced MSD unfolding and bond lifetime are necessary for inducing α-type Ca2+ signal.
Using a sensitivity-specificity analysis (Figure 5—figure supplement 2A,B) that slides a putative threshold through the tmax vs. Ca2+ type data, we found tmax>2 s to be the best predictor for A1WT to trigger α- rather than β-type Ca2+ (Figure 5E, dashed lines), which agrees with the fact that the 2 s threshold is much shorter than the average tmax of α-type Ca2+, but exceeds that of β-type Ca2+ and MSD time-to-unfold (Figure 5—figure supplement 2C). Thus, a longer-lived pre-Ca2+ bond favors MSD unfolding, thereby triggering α-type Ca2+; otherwise, it only triggers β-type Ca2+. Together, our data suggests separate roles of LRRD and MSD unfolding in GPIbα signaling, with the former intensifying the Ca2+ level and the latter determining the Ca2+ signal type .
Perturbing cytoplasmic association of GPIbα with 14-3-3ζ inhibits mechanotransduction
To understand GPIbα-mediated mechanosensing requires analysis of not only ligand binding and domain unfolding in the extracellular segment of GPIbα, but also events in its cytoplasmic region. 14-3-3ζ is a cytoplasmic protein that has direct association with both GPIbα and GPIbβ C-termini (Calverley et al., 1998) and regulates GPIb signal transduction(Dai et al., 2005). To investigate the role of 14-3-3ζ in GPIbα-mediated Ca2+ signaling, we perturbed the system with a myristoylated peptide (MPαC) that mimics the 14-3-3ζ binding sequence of GPIbα, thereby blocking the association of 14-3-3ζ with GPIbα cytoplasmic tail. Consistent with the previously reported signaling inhibition effect (Dai et al., 2005; Yin et al., 2013), MPαC reduced the fraction of α-type Ca2+ from 34 to 3% without affecting β-type Ca2+, whereas a scramble peptide MαCsc had no effect (Figure 7E). Similar results were obtained by pulling GPIbα via WM23 on platelets (Figure 7F). Thus, GPIb–14-3-3ζ association, a biochemical event, is crucial for the transduction of MSD unfolding, a mechanical event, into intracellular signals. These observations indicate that GPIb–14-3-3ζ serves, at least in part, as a mechanotransducer (Figure 7G).
Discussion
The mechanoreception of GPIbα has been supported by direct observations of transient intracellular Ca2+ spike (termed type α/β peak) upon platelet translocation on VWF under flow (Mazzucato et al., 2002; Nesbitt et al., 2002). However, many questions remain. Using fBFP real-time single-bond, single-platelet analysis of force-regulated ligand binding kinetics, receptor unfolding dynamics and intraplatelet calcium mobilization, we have: 1) identified, characterized and mathematically modeled the force- and ligand-dependent cooperativity between LRRD and MSD unfolding (Figure 7G); 2) defined an optimal magnitude and threshold duration of clamped force for platelet signal initiation via a single GPIbα bond (Figure 7—figure supplement 2); 3) uncovered a mechanopresentation defect in a type 2B VWD mutant A1R1450E; 4) delineated the interplay among ligand engagement, GPIbα domain unfolding and signal triggering; and 5) revealed inhibition of GPIbα mechanotransduction by perturbing its cytoplasmic association with 14-3-3ζ.
It is an interesting yet challenging problem to define the minimum mechanical stimulation for inducing signal transduction. We demonstrated that a single A1–GPIbα bond can induce calcium in a platelet without clustering GPIb by multimeric ligands. Specifically, in 83% of the cases where α-type Ca2+ was triggered, only a single MSD unfolding event was observed before the Ca2+ onset. Of these, 43% had only one pre-Ca2+ bond lifetime event. Thus, pulling a single GPIbα by a 25 pN force for >2 s to unfold MSD once is necessary and sufficient to induce α-type Ca2+ signals. By comparison, to trigger calcium in a naïve CD8+ T lymphocyte requires a sequence of intermittent bonds with a total of >10 s lifetimes under 10 pN accumulated in the first 60 s of contacts between T cell receptors and agonist peptide-major histocompatibility complex molecules (Liu et al., 2014). In both cases a threshold of force duration is required, as ramped forces without a clamp phase are unable to trigger appreciable levels of α-type Ca2+ regardless of its magnitude.
The binding defects of VWF-A1 with type 2B VWD mutations have long been recognized (Ruggeri, 2004). A recent study has shown a type 2B mutation, A1V1316M, causes additional signaling defects (Casari et al., 2013). We showed that another type 2B mutation, A1R1450E, also has signaling defect. Interestingly, the defect to induce calcium has the same root as the binding defect, namely, the conversion of the wild-type catch-slip bond to the mutant slip-only bond (Ju et al., 2013; Yago et al., 2008). Consequently, force exerted on GPIbα by A1R1450E is less able to unfold LRRD, lasts shorter at 25 pN to unfold MSD less frequently, does not generate unbinding cooperativity between the LRRD and MSD, and induces lower level and frequency of α-type Ca2+ at >10 pN. Thus, the mechanical requirements for signal induction manifest as force-dependencies of VWF–GPIbα bond lifetime, MSD unfolding frequency, unfolding cooperativity, and Ca2+ level/type that display similar patterns for the same A1 construct (WT or R1450E) but different patterns between A1WT and A1R1450E. These findings show that the GPIbα mechanoreceptor can discriminate ligands and shed light to the biophysical mechanisms of type 2B VWD.
Our new data on the interplay among VWF binding, GPIbα unfolding, and Ca2+ signaling have provided new insights into the inner workings of the A1–GPIb–14-3-3ζ molecular assembly (Figure 7; Video 1). By residing in the juxtamembrane stalk region, the MSD has been shown to be mechanosensitive, and hypothesized to play a role in activating platelets (Zhang et al., 2015). In the present work, we found that MSD unfolding is required to trigger α-type Ca2+, showing that this extracellular mechanical event is necessary for transduction of the information embedded in the force waveform into intracellular biochemical signals via the 14-3-3ζ connection (Figure 7G). By overlapping with the ligand-binding site, the LRRD can feel the structural variation in the A1 and respond with an altered unfolding frequency and changed bond lifetime. Importantly, LRRD unfolding prolongs A1–GPIbα bond lifetime to facilitate MSD unfolding, thereby increasing the frequency of α-type Ca2+ and its level (Figure 7—figure supplement 2). Thus, our study has elucidated part of a mechanosensor that includes three components: 1) a MSD in the juxtamembrane region whose conformational change results in a binary decision of Ca2+ type, 2) a LRRD in the ligand-binding region whose conformational change leads to continuous alterations in ligand-binding duration, signal level and fractions of different signal types, and 3) a MP stalk that transmits force over a distance and provides coupling between the two unfoldable domains. The differential unfolding behaviors of the LRRD and MSD in response to distinct force waveforms provide a simple mechanical mechanism for unfolding cooperativity, by setting the response order such that LRRD unfolds first during force ramp to give more time for MSD to unfold during force clamp. These principles may be helpful for design of a generic mechanosensor, e.g., using synthetic biology approaches.
In addition to an increased MSD unfolding frequency, the cooperativity between LRRD and MSD unfolding may manifest as an increased LRRD unfolding extent. This has been suggested by the observation that the values of the first two peaks (20 and 36 nm) do not add up to that of the third (65-70 nm) in the unfolding length histograms. We note that the observed maximum unfolding length from the GC experiments (56 nm) is smaller than the calculated contour length of LRRD (70 nm). Since LRRD consists of 8 leucine-rich repeats, this suggests that only some but not all of the repeats were unfolded in any given observation. Thus, the LRRD dataset contains mixed populations of partial unfolding events. By comparison, the MSD dataset may be a more uniform population as the observed maximum unfolding length from the WM23 experiments (27 nm) matches the calculated contour length of MSD (26 nm). These considerations suggest possible explanations for the observation that events in which both LRRD and MSD unfold generate more length than the sum of lengths generated from events in which either LRRD or MSD unfolds: LRRD unfolding events with a higher number of unfolded leucine-rich repeats may facilitate MSD unfolding more effectively than those with a lower number of unfolded leucine-rich repeats. Alternatively, MSD unfolding, once happens, may induce more leucine-rich repeats in LRRD to unfold. Note that these two mechanisms are not mutually exclusive.
Studies in mechanosensitive ion channels and enzymes have provided knowledge on how biomolecules respond to force and transduce mechanical stimulations into biochemical signals. For example, ion channels open and close in response to stress within the lipid bilayer or force within a protein link that can do work on the channel and stabilize its state (Sukharev and Sachs, 2012). Mechanosensitive enzymes or substrates, such as vinculin (del Rio et al., 2009; Grashoff et al., 2010) or A2 domain of VWF (Wu et al., 2010; Zhang et al., 2009b), change conformations in response to forces to expose a cryptic site to enable enzymatic reaction. By comparison, for the system studied herein, force signals are received by the receptor via ligand interaction, hence mediated by their binding kinetics. The process of transducing the extracellular mechanical events (i.e., LRRD and/or MSD unfolding) across the cell membrane is likely mechanical rather than chemical (i.e., ion influx).
The study of GPIb mechanosensing may help understand how mechanical force regulates platelet thrombotic functions. For example, in response to shear gradients resulting from flow perturbations, discoid platelets aggregate rapidly in a manner independent of agonist activation pathways (Jain et al., 2016; Nesbitt et al., 2009; Yong et al., 2011). This intriguing phenomenon has significant implication in atherothrombosis and medical device thrombotic fouling. Other than the requirement for VWF–GPIbα binding, the underlying mechanism of such purely force-induced platelet thrombosis remains elusive. Here LRRD unfolding may play a role because it requires an increasing force (resembles shear gradient) and strengthens VWF–GPIbα bonds (Ju et al., 2015b) (Figure 7—figure supplement 2C). In addition, our findings may have broad implications since LRRD is a common structure shared by many adhesion and signaling receptors, e.g., toll-like receptors (Bella et al., 2008).
Materials and methods
Proteins and peptides
Request a detailed protocolRecombinant monomeric VWF-A1 (residues 1238–1471) WT and type 2B mutant R1450E (Cruz et al., 2000; Morales et al., 2006) generated by E.coli were gifts of Miguel A. Cruz (Baylor College of Medicine, Houston, TX). Glycocalicin was purified from outdated platelets (Fox et al., 1988). Three anti-GPIbα mAbs were used: AK2 (Abcam, Cambridge, MA), AN51 (Millipore, Billerica, MA) and WM23 (a gift from Michael Berndt, Curtin University, WA, Australia and Renhao Li, Emory University, Atlanta, GA). Anti-GPIbβ mAb LS-B3174 was purchased (LifeSpan BioSciences, Seattle, WA). Anti-VWF-A1 mAb 6G1 was a gift from Michael Berndt. Myristoylated peptides (MPαC, C13H27CONH-SIRYSGHpSL) and myristoylated scrambled control peptide (MαCsc, C13H27CONH-LSISYGSHR) were produced as previously described (Dai et al., 2005; Yin et al., 2013).
Red blood cells (RBCs) and platelets
Request a detailed protocolHuman RBCs and platelets were collected abiding a protocol approved by the Institute Review Broad of Georgia Institute of Technology. RBCs were prepared as previously described (Ju et al., 2013a). To obtain fresh discoid human platelets, whole blood was drawn slowly from a vein of a healthy volunteer to fill in a 3 ml syringe preloaded with 0.5 ml ACD buffer (6.25 g sodium citrate, 3.1 g citric acid anhidrous, 3.4 g D-glucose in 250 ml deionized H2O, pH 6.7). The whole blood was centrifuged at 150 g for 15 min without brake. Platelet-rich plasma was extracted and centrifuged at 900 g for another 10 min. The platelet pellet was resuspended into Hepes-Tyrode buffer (134 mM NaCl, 12 mM NaHCO3, 2.9 mM KCl, 0.34 mM sodium phosphate monobasic, 5 mM HEPES, and 5 mM glucose, 1% bovine serum albumin (BSA), pH 7.4). For Ca2+ imaging experiments, isolated platelets were incubated with Fura-2-AM (Life Technologies, Grand Island, NY) at 30 μM for 30 min. For treatment with MPαC or MαCsc, the peptide pre-dissolved in DMSO was resuspended into the platelet suspension to reach a final concentration of 25 μM and incubated at 37°C for 30 min.
Functionalization of glass beads
Request a detailed protocolA1WT, A1R1450E and antibodies were pre-coupled covalently with maleimide-PEG3500-NHS (MW ~3500; JenKem, TX). As previously described (Ju et al., 2013a, 2015a), the modified proteins were then mixed with streptavidin (SA)-maleimide (Sigma-Aldrich, St. Louis, MO) in carbonate/bicarbonate buffer (pH 8.5) and together linked to silanized borosilicate beads (Thermo Fisher Scientific, Waltham, MA) in phosphate buffer (pH 6.8). Site densities of ligands on beads were measured using the previously described flow cytometry method (Ju et al., 2015a).
Fluorescence biomembrane force probe (fBFP)
Request a detailed protocolOur fBFP was developed to simultaneously measure the binding kinetics of single receptor–ligand bonds (Ju et al., 2015a, 2013, 2015c) and the mechanics of single protein conformational changes (Chen et al., 2012; Ju et al., 2015b), as did our original BFP, and receptor-initiated intracellular signaling with a concurrent fluorescent imaging module (Chen et al., 2015; Liu et al., 2014). Bond formation, force application, receptor conformational change, and bond dissociation were enabled and monitored in controlled BFP cycles of a few seconds each. Intraplatelet calcium fluxes were ratiometracally imaged as a signaling readout.
In a BFP cycle, the platelet was driven to approach, impinge and hold the probe with a 20–30 pN compressive force for a contact time of 2 s to allow for bond formation, and then retract (ramp) for bond detection (Figure 1C and Figure 1—figure supplement 1A i–iv). Displacement of the probe bead was tracked, which reflected the force exerted on it. During the ramping phase, a bond event was signified by a tensile force signal (Figure 1—figure supplement 1C,D), while no tensile force was detected in a no-bond event (Figure 1—figure supplement 1B). Bond and no-bond events were enumerated to calculate an adhesion frequency in 50 repeated cycles for each bead and platelet pair. At least 3 bead–platelet pairs were measured and their adhesion frequencies were used to calculate mean ± s.e.m. values. To define the minimum requirement for GPIb mechanoreception, adhesions were adjusted to be infrequent (<20%) by titrating the densities of randomly distributed A1 and mAb on the probe beads (Figure 1D). This ensured that most (>89%) platelet–bead binding events were mediated by non-clustered single-bonds (Chesla et al., 1998).
To quantify intraplatelet Ca2+ mobilization, we used ratiometric imaging with a light source that alternates two excitation wavelengths (340 ± 10 nm to excite Ca2+-bound Fura-2, and 380 ± 10 nm to excite Ca2+-free Fura-2). The emission light from the excited Fura-2 (with or without Ca2+ binding) was captured by a fluorescence camera. To maintain the physiological temperature (37°C) inside the cell chamber, a custom-designed temperature control system made in house was integrated into the fBFP. Details about the Ca2+ imaging analysis and temperature control have been previously described (Chen et al., 2015; Liu et al., 2014).
Force-clamp and force-ramp experiment modes
Request a detailed protocolIn the force-clamp mode, the target pipette was driven to repeatedly contact the probe bead for 2 s and retract at a constant speed (3.3 μm/s). Multiplying the BFP spring constant (0.3 pN/nm), this would translate to a linearly increasing force at a constant ramping rate (1000 pN/s). Upon detection of bond event, a feedback loop controls the retraction so that it would be paused at a desired clamped force (10, 25, 40 and 60 pN) to wait for bond dissociation (Figure 1—figure supplement 1C). After that the target pipette would return to the original position to complete the cycle (Figure 1—figure supplement 1Av–vii). Each platelet was interrogated for a continuous time of 200 s to generate a force spectroscopy trace exemplified in Figure 1C before changing to a new pair of BFP bead and platelet. Lifetimes were measured from the instant when the force reached the desired level to the instant of bond dissociation (Figure 1C) (Ju et al., 2015a, 2013). In the force-ramp mode, the force was loaded at different ramping rates (100, 1000 or 10,000 pN/s). The target was retracted continuously until bond rupture without holding at a constant force (Figure 1—figure supplement 1D).
Protein unfolding analysis
Request a detailed protocolUnfolding of GPIbα in the ramping phase was signified by a sudden force stagnation or drop (kink) as opposed to the linearly increasing force signals (Figure 1E). To determine the unfolding length, we derived a force vs. extension curve (Figure 1E inset) from the differential displacement between the BFP tracking system (probe position) and the piezoelectric actuator feedback system (target position) as previously described (Chen et al., 2012). The unfolding length was given by the sudden extension increase without force increase, the result of which was comparable to the differential contour length derived by fitting the prior- and post-unfolding force-extension curves with the worm-like chain (WLC) model (Bustamante et al., 1994) (Figure 1—figure supplement 2).
To reveal distinct populations, we used the nonparametric kernel density estimation to detect peaks in the data distribution of ensemble ramped unfolding lengths (Freedman and Diaconis, 1981) (Figure 2—figure supplement 1A,B) and used a reliable data-based bandwidth selection method (Sheather and Jones, 1991) to determine the optimal bin width as 5 nm for the data in Figure 2D,F. Using a 9.53 nm bin width determined by the Freedman-Diaconis formula and Freedman and Diaconis’ heuristic rule (Freedman and Diaconis, 1981) for histogram analysis also revealed three peaks, although the valley separating the first two peaks consists of a single low fraction bin only (Figure 2—figure supplement 1C). The first two peaks in Figure 2D were suggested as MSD and LRRD unfolding, respectively, based on their favorable comparisons to the respective WM23 vs. platelet and A1WT vs. GC data in Figure 2E, which allowed only MSD or LRRD to unfoldrespectively. To test these hypotheses, we analyzed the molecular mechanics by sorting the unfolding forces and lengths into subgroups from the respective WM23 vs. platelet and A1WT vs. GC experiments, plotted the average unfolding force vs. length data, and fitted the data to the WLC model (Zhang et al., 2009a, 2015). The best-fit curves were then served as standards to calibrate the average unfolding force vs. length data from the A1WT vs. platelet experiment (Figure 2H,I). The agreement between the data and the calibrated WLC curves rigorously verified the hypothetical identities of the first two peaks in Figure 2D.
Unfolding of GPIbα in the clamping phase was signified by a sudden force decrease (Figure 1E). The unfolding length was calculated from force change divided by BFP spring constant (Figure 1F), similar to the integrin extension length measurement from the previous distance-clamp analysis (Chen et al., 2012, 2016). The duration from the beginning of the clamping phase to the beginning of unfolding was the time-to-unfold, tu (Figure 1F).
Testing hypothesis for cooperative unfolding
Request a detailed protocolThe cooperativity between LRRD and MSD unfolding at a clamped force (e.g. 25 pN) was determined by testing the null hypothesis that the two domains unfolded independently. The frequencies of LRRD unfolding and MSD unfolding pulled by A1WT were calculated using data from Figure 3—source data 1A (take 25 pN for example):
MSD = '+' | MSD = '−' | Row total | |
---|---|---|---|
LRRD = '+' | 14 (9.7%) | 11 (7.6%) | 25 (17.2%) |
LRRD = '−' | 25 (17.2%) | 95 (65.5%) | 120 (82.8%) |
Column total | 39 (26.9%) | 106 (73.1%) | 145 (100%) |
The Pearson's χ2 test was used to test the null hypothesis (H0) that LRRD unfolding and MSD unfolding are independent. The χ2 statistic was calculated as follow:
where is the observed count and is the expected count under the null hypothesis. The subscripts i and j denote LRRD and MSD respectively, whose values 1 and 2 denote with (+) and without (-) unfolding respectively. The system has (2−1) × (2−1) = 1 degree of freedom. The small p-value (3.09 × 10–4 from the above χ2) requires that we reject the null hypothesis and accept the alternative hypothesis that LRRD and MSD unfolding are not independent. In other words, cooperativity exists between LRRD and MSD unfolding when GPIbα on platelets was pulled by A1WT. Similar statistical analyses were employed to assess cooperative unfolding between LRRD and MSD when GPIbα was pulled by A1WT, A1R1450E or AN51 at different clamped forces. The levels of significance were presented as –log10 (p-values) (Figure 3G,H).
Model for observing MSD clamped unfolding
Request a detailed protocolThe measured MSD time-to-unfold tu distributed as a single exponential decay: , where pu is the probability density and ku is the unfolding rate of MSD under a clamping force. By fitting the semi-log plotted experimental distribution with a straight line (Figure 4—figure supplement 1A), the unfolding rate at 25 pN was evaluated from the negative slope or the reciprocal average time-to-unfold, ku = 1/<tu> = 0.870s-1. Modeling the force-dependent MSD unfolding rate (Figure 4A) by the Bell equation (Bell, 1978), we found the zero-force unfolding rate and the width of the energy barrier to be 0.26 s−1 and 0.242 nm for pulling GPIbα via A1 on a live platelet. The first value is much larger and the second value is much smaller than the respective values previously obtained using an optical tweezer to measure ramped unfolding of MSD in purified GPIbα constructs (Zhang et al., 2015).
We previously reported that the A1–GPIbα bond lifetime tb distributed as a dual exponential decay with a fast- and a slow-dissociating off-rate (Ju et al., 2013). We also recently showed that unfolding of LRRD prolongs A1–GPIbα bond lifetime (Ju et al., 2015b). These results were also observed in this work, which comprise individual bond lifetime measurements that give rise to the averaged results in Figure 4B,C. Therefore, the probability densities for a A1–GPIbα bond, with and without prior LRRD unfolding in the ramping phase, to dissociate at time tb during the clamping phase are:
where kij and wij (wi1 + wi2 = 1) denote, respectively, off-rates and associated fractions of bonds under a clamped force, with the first subscript indicating without (Equation 1) or with (Equation 2) a prior LRRD unfolding event and the second subscript indicating the fast (Equation 1) or slow (Equation 2) dissociation pathway. By fitting the above model to the lifetime ensemble data, the parameters were calculated (Figure 4—source data 1).
Assuming that MSD unfolding and A1–GPIbα unbinding are independent events, the joint probability density for MSD unfolding at time tu and A1–GPIbα unbinding at time tb is where i = 1, 2 depending on whether LRRD unfolding occurs. This joint probability is depicted as a surface in Figure 4D and Figure 4—figure supplement 1B, using respective A1WT and A1R1450E data measured at 25 pN clamped force. The condition for observing MSD clamped unfolding is that the A1–GPIbα bond lifetime tb lasts longer than the time-to-unfold tu. Thus, the probability of observing MSD unfolding in the clamping phase Pui is the volume under the probability density surface over the region 0<tu<tb<∞, which is marked by the vertical red planes in Figure 4D and Figure 4—figure supplement 1B. For instance, in the absence of prior LRRD unfolding, the probability of observing MSD unfolding in the clamping phase of 25 pN is:
Similarly, in the presence of prior LRRD unfolding,
The model was applied to predict the MSD clamped unfolding probabilities under different clamped forces pulled by A1WT and A1R1450E (Figure 4E,F). For A1R1450E, the ensemble MSD clamped unfolding events were no longer segregated into LRRD- and LRRD+ groups, because few MSD clamped unfolding events occurred following LRRD unfolding due to the reduced bond lifetime.
Threshold and sensitivity-specificity analysis
Request a detailed protocolWe used the sensitivity-specificity analysis to solve the optimal threshold t0 for pre-Ca2+ longest lifetime tmax separating α and β Ca2+ types. There are 4 possible outcomes (Figure 5—figure supplement 2A): a false positive (FP) happens when tmax>t0 and a β-type Ca2+ was observed; a false negative (FN) happens when tmax≤t0 and an α-type Ca2+ was observed; a true positive (TP) happens when and an α-type Ca2+ was observed, and a true negative (TN) happens when tmax≤t0 and a β-type Ca2+ was observed. The sensitivity or true positive rate defines the fraction of true positive among all positive results, TP/(TP+FN), whereas the specificity or true negative rate defines the fraction of true negative among all negative results, TN/(TN+FP). The optimal threshold is solved by minimizing the total counts of false positive and false negative. To do that, a receiver operating characteristic (ROC) curve was created by plotting the TP rate (sensitivity) against the FP rate (1- specificity) at various tmax values from which the optimal threshold t0 that achieved the best sensitivity and specificity was identified (Figure 5—figure supplement 2B).
Statistical analysis
Request a detailed protocolTwo-tailed Students’ t-test was used to assess significance for group comparisons. Pearson correlation coefficient was used as a measure for linear dependency between two variables (Ca2+ level and binding kinetics).
To determine errors in classification of different Ca2+ types (null, β, α) and in identification of unfolding (no unfolding, LRRD, MSD, MSD+LRRD), we assume that observed counts follow a multinomial distribution with total counts and event probabilities . We then use the fraction of the i-th category, ni/n, as an estimate for the i-th event probability and as the associated standard error s.e.m.
References
-
The leucine-rich repeat structureCellular and Molecular Life Sciences 65:2307–2333.https://doi.org/10.1007/s00018-008-8019-0
-
Entropic elasticity of lambda-phage DNAScience 265:1599–1600.https://doi.org/10.1126/science.8079175
-
Human signaling protein 14-3-3zeta interacts with platelet glycoprotein Ib subunits Ibalpha and IbbetaBlood 91:1295–1303.
-
von Willebrand factor mutation promotes thrombocytopathy by inhibiting integrin αIIbβ3Journal of Clinical Investigation 123:5071–5081.https://doi.org/10.1172/JCI69458
-
Observing force-regulated conformational changes and ligand dissociation from a single integrin on cellsJournal of Cell Biology 199:497–512.https://doi.org/10.1083/jcb.201201091
-
Measuring two-dimensional receptor-ligand binding kinetics by micropipetteBiophysical Journal 75:1553–1572.https://doi.org/10.1016/S0006-3495(98)74074-3
-
Mapping the glycoprotein Ib-binding site in the von willebrand factor A1 domainJournal of Biological Chemistry 275:19098–19105.https://doi.org/10.1074/jbc.M002292200
-
Structure of the glycoprotein Ib.IX complex from platelet membranesJournal of Biological Chemistry 263:4882–4890.
-
On the histogram as a density estimator: L 2 theoryProbability Theory and Related Fields 57:453–476.https://doi.org/10.1007/BF01025868.
-
Arterial thrombosis—insidious, unpredictable and deadlyNature Medicine 17:1423–1436.https://doi.org/10.1038/nm.2515
-
The N-terminal flanking region of the A1 domain regulates the force-dependent binding of von Willebrand factor to platelet glycoprotein IbαJournal of Biological Chemistry 288:32289–32301.https://doi.org/10.1074/jbc.M113.504001
-
Transport regulation of two-dimensional receptor-ligand associationBiophysical Journal 108:1773–1784.https://doi.org/10.1016/j.bpj.2015.02.023
-
Specific inhibition of ectodomain shedding of glycoprotein Ibα by targeting its juxtamembrane shedding cleavage siteJournal of Thrombosis and Haemostasis 11:2155–2162.https://doi.org/10.1111/jth.12425
-
Distinct glycoprotein Ib/V/IX and integrin alpha IIbbeta 3-dependent calcium signals cooperatively regulate platelet adhesion under flowJournal of Biological Chemistry 277:2965–2972.https://doi.org/10.1074/jbc.M110070200
-
Type IIB von Willebrand disease: a paradox explains how von Willebrand factor worksJournal of Thrombosis and Haemostasis 2:2–6.https://doi.org/10.1111/j.1538-7836.2003.00523.x
-
A reliable data-based bandwidth selection method for kernel density estimationJournal of the Royal Statistical Society Series B 53:683–690.https://doi.org/10.2307/2345597
-
Molecular force transduction by ion channels: diversity and unifying principlesJournal of Cell Science 125:3075–3083.https://doi.org/10.1242/jcs.092353
-
Platelet glycoprotein Ibalpha forms catch bonds with human WT vWF but not with type 2B von Willebrand disease vWFThe Journal of Clinical Investigation 118:3195–3207.https://doi.org/10.1172/JCI35754
-
Role for platelet glycoprotein Ib-IX and effects of its inhibition in endotoxemia-induced thrombosis, thrombocytopenia, and mortalityArteriosclerosis, Thrombosis, and Vascular Biology 33:2529–2537.https://doi.org/10.1161/ATVBAHA.113.302339
-
Measuring receptor/ligand interaction at the single-bond level: experimental and interpretative issuesAnnals of Biomedical Engineering 30:305–314.https://doi.org/10.1114/1.1467923
Article and author information
Author details
Funding
Diabetes Australia (IRMA G179720)
- Lining Ju
University of Sydney (2016 Sydney Medical School ECR Kickstart Grant)
- Lining Ju
National Science Foundation (DMS-1505256)
- Lingzhou Xue
National Heart, Lung, and Blood Institute (HL062350)
- Xiaoping Du
National Heart, Lung, and Blood Institute (HL080264)
- Xiaoping Du
National Heart, Lung, and Blood Institute (HL125356)
- Xiaoping Du
National Heart, Lung, and Blood Institute (HL132019)
- Cheng Zhu
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank MA Cruz for providing VWF-A1 proteins; MC Berndt and R Li for providing the WM23 antibody and personal communications; W Lam lab for the blood draw; SP Jackson lab for constructive comments and financial support to prepare the manuscript; J Lou for discussion on GPIbα structure and WLC model fitting; F Zhou for generating 3D probability density surface plots. B Liu for sharing the Ca2+ analysis Matlab codes; W Chen, J Hong, Q Ji, and C Ge for fBFP instrumentation supports; J-F Dong, RP McEver, LV McIntire, R Andrews, H Lu and J Liao for helpful comments. This work was supported by NIH grant HL132019 (CZ), HL062350 (XD), HL080264 (XD) and HL125356 (XD), Diabetes Australia research grant G179720 (LJ), Sydney Medical School 2016 early-career researcher kickstart grant (LJ), and NSF grant DMS-1505256 (LX).
Ethics
Human subjects: Human RBCs and platelets for BFP experiments were collected abiding a protocol (#H12354) approved by the Institute Review Broad of Georgia Institute of Technology. Informed consent was obtained from each blood donor.
Copyright
© 2016, Ju et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 2,577
- views
-
- 578
- downloads
-
- 67
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Structural Biology and Molecular Biophysics
Ciliary rootlets are striated bundles of filaments that connect the base of cilia to internal cellular structures. Rootlets are critical for the sensory and motile functions of cilia. However, the mechanisms underlying these functions remain unknown, in part due to a lack of structural information of rootlet organization. In this study, we obtain 3D reconstructions of membrane-associated and purified rootlets from mouse retina using cryo-electron tomography. We show that flexible protrusions on the rootlet surface, which emanate from the cross-striations, connect to intracellular membranes. In purified rootlets, the striations were classified into amorphous (A)-bands, associated with accumulations on the rootlet surface, and discrete (D)-bands corresponding to punctate lines of density that run through the rootlet. These striations connect a flexible network of longitudinal filaments. Subtomogram averaging suggests the filaments consist of two intertwined coiled coils. The rootlet’s filamentous architecture, with frequent membrane-connecting cross-striations, lends itself well for anchoring large membranes in the cell.
-
- Structural Biology and Molecular Biophysics
Although the αC-β4 loop is a stable feature of all protein kinases, the importance of this motif as a conserved element of secondary structure, as well as its links to the hydrophobic architecture of the kinase core, has been underappreciated. We first review the motif and then describe how it is linked to the hydrophobic spine architecture of the kinase core, which we first discovered using a computational tool, local spatial Pattern (LSP) alignment. Based on NMR predictions that a mutation in this motif abolishes the synergistic high-affinity binding of ATP and a pseudo substrate inhibitor, we used LSP to interrogate the F100A mutant. This comparison highlights the importance of the αC-β4 loop and key residues at the interface between the N- and C-lobes. In addition, we delved more deeply into the structure of the apo C-subunit, which lacks ATP. While apo C-subunit showed no significant changes in backbone dynamics of the αC-β4 loop, we found significant differences in the side chain dynamics of K105. The LSP analysis suggests disruption of communication between the N- and C-lobes in the F100A mutant, which would be consistent with the structural changes predicted by the NMR spectroscopy.