Von Willebrand factor A1 domain stability and affinity for GPIbα are differentially regulated by its O-glycosylated N- and C-linker

  1. Klaus Bonazza
  2. Roxana E Iacob
  3. Nathan E Hudson
  4. Jing Li
  5. Chafen Lu
  6. John R Engen
  7. Timothy A Springer  Is a corresponding author
  1. Program in Cellular and Molecular Medicine, Boston Children’s Hospital and Department of Biological Chemistry and Molecular Pharmacology, Harvard Medical School, United States
  2. Department of Chemistry and Chemical Biology, Northeastern University, United States
  3. Department of Physics, East Carolina University, United States

Abstract

Hemostasis in the arterial circulation is mediated by binding of the A1 domain of the ultralong protein von Willebrand factor (VWF) to GPIbα on platelets to form a platelet plug. A1 is activated by tensile force on VWF concatemers imparted by hydrodynamic drag force. The A1 core is protected from force-induced unfolding by a long-range disulfide that links cysteines near its N- and C-termini. The O-glycosylated linkers between A1 and its neighboring domains, which transmit tensile force to A1, are reported to regulate A1 activation for binding to GPIb, but the mechanism is controversial and incompletely defined. Here, we study how these linkers, and their polypeptide and O-glycan moieties, regulate A1 affinity by measuring affinity, kinetics, thermodynamics, hydrogen deuterium exchange (HDX), and unfolding by temperature and urea. The N-linker lowers A1 affinity 40-fold with a stronger contribution from its O-glycan than polypeptide moiety. The N-linker also decreases HDX in specific regions of A1 and increases thermal stability and the energy gap between its native state and an intermediate state, which is observed in urea-induced unfolding. The C-linker also decreases affinity of A1 for GPIbα, but in contrast to the N-linker, has no significant effect on HDX or A1 stability. Among different models for A1 activation, our data are consistent with the model that the intermediate state has high affinity for GPIbα, which is induced by tensile force physiologically and regulated allosterically by the N-linker.

Editor's evaluation

Through a series of elegant experiments, the authors concluded that the N-terminal linker together with its O-linked glycans stabilize the A1 in an inactive, native conformation. Removal of the N-linker or its O-linked glycans destabilizes the A1 structure and may activate the A1 into an intermediate, higher affinity state for GPIbα binding. The authors also concluded that the C-terminal linker and its O-linked glycans do not contribute to stabilizing the A1 structure, though they prevent GPIbα binding through steric hindrance and electrostatic repulsion. The conclusions are consistently supported by the data from multiple experimental techniques. The reported results are highly interesting for the hematology and biophysics communities.

https://doi.org/10.7554/eLife.75760.sa0

Introduction

The ultra-long length of the blood plasma protein von Willebrand factor (VWF) enables its activation by hydrodynamic forces at sites of hemorrhage (Springer, 2014; Figure 1A, B). VWF, a long concatemer of monomers linked head-to-head and tail-to-tail, undergoes a change of its multimeric superstructure from a random coil form (Parker and Lollar, 2021) to a thread-like, extended form when exposed to elongational flow (Bonazza et al., 2015; Fu et al., 2017). Hydrodynamic flow activates VWF to bind to glycoprotein Ib (GPIb) on platelets, the essential interaction required for hemostasis and thrombosis in the arteriolar circulation. Fluorescent imaging of single VWF concatemers tethered in a flow cell showed that binding of GPIbα not only requires elongation of VWF but also a half-maximal tensile force of ~20 pN transmitted through the backbone of VWF concatemers (Fu et al., 2017). Activation could be well fit to a model of transition from a low-affinity state 1 to a high-affinity state 2 (Figure 1B). Measuring A1–GPIbα binding and unbinding with laser tweezers also revealed two states, with a force dependent switch between two states with different on- and off-rates (Kim et al., 2015; Kim et al., 2010). In addition to such flex-bond measurements, A1 and GPIbα are also reported to form catch bonds, although these measurements lack single molecule fiduciary markers or fits to models from which on- and off-rates could be extracted (Ju et al., 2013). Despite multiple crystal structures of the A1 domain bound to GPIb, distinct conformational states of A1 that explain force-dependent regulation of A1 affinity for GPIbα have yet to emerge (Blenner et al., 2014; Dumas et al., 2004).

von Willebrand factor (VWF), hydrodynamic flow, A1 activation, and A1 linkers.

(A) N- and C-terminal linkage and domain organization of a VWF monomer in a VWF concatemer. Monomers are connected head-to-head and tail-to-tail. (B) When exposed to elongational flow and tethered (black triangle) on a vessel wall, VWF concatemers extend to a linear shape. At higher elongational flow, which exerts higher mechanical tension along the length of the concatemer (highest at the tether point and zero at the downstream end), A1 undergoes transition from a low- to a high-affinity state as a consequence of tensile forces transmitted to it through its linkers. The disulfide bond between Cys1272 and Cys1458 is indicated with red SS. (C) A1 protein constructs studied here that differ in length of linkers and were expressed in E. coli (not glycosylated) or in mammalian cells (O-glycosylated) (O-glycans shown as lollipops). The long-range disulfide bond is schematized in red. (D) A1 constructs were subjected to sodium dodecyl–sulfate polyacrylamide gel electrophoresis (SDS–PAGE) and stained with Coomassie Blue. Amount of protein added per lane differed between protein constructs but was identical for each construct on reducing and nonreducing PAGE.

The force-responsive A1 domain is unusual among protein domains in being linked to neighboring domains through O-glycosylated, mucin-like segments (Figure 1A, blue hashed lines). Tensile force exerted on A1 is transmitted through these linkers. It has been suggested that A1-flanking sequences (Auton et al., 2012; Interlandi et al., 2017; Ju et al., 2013; Miyata and Ruggeri, 1999; Sugimoto et al., 1993) or domains (Aponte-Santamaría et al., 2015; Ulrichts et al., 2006) may shield the GPIbα-binding epitope and thereby inhibit otherwise tight binding to GPIb. It was also proposed that the N- and C-linkers interact with each other and form an autoinhibitory module that masks A1 (Deng et al., 2017). A force-dependent signature for breakage of this module has been reported that correlates with the combined expected extension of the N- and C-linkers (Arce et al., 2021). An alternative mechanism has also been proposed by which the N-linker could regulate A1 affinity, that is, by regulating the relative stability of the A1 native and intermediate states. A1 has been found to have three states, native (N), intermediate (I), and denatured (D), while the homologous A2 and A3 domains in VWF have only two states, native and denatured (Auton et al., 2007). Lower stability of the A1 native state correlates with higher affinity (Auton et al., 2012; Auton et al., 2010a). This anticorrelation between A1 affinity and stability is also seen with von Willebrand disease type 2B gain-of-function mutations, which increase A1 affinity and decrease stability (Auton et al., 2010b; Tischer et al., 2017). Conversely, the G1342S decrease-in-function type 2M von Willebrand disease mutation decreases A1 affinity for GpIbα and increases A1 stability (Auton et al., 2009).

Mucin-like regions in proteins have unique characteristics distinct from both intrinsically disordered polypeptide segments and folded domains. The bulky O-linked glycans which include sialic acid that are attached to threonine and serine residues in mucins are highly solvated, repel one another, and make mucins extended, with an average length of ~2 Å/residue (Clemetson, 1983; Woollett et al., 1985). Electron microscopy showed that D3 and A1 in VWF, which are connected by the N-linker, have many possible orientations, and an interdomain distance as long as ~7 nm (Zhou et al., 2011). A1 is mechanically stabilized by a long-range disulfide bond between Cys1272 and Cys1458 (Figure 1B, C). Under tension, the core of A1 between these two cysteines, residues 1273–1457, is protected from complete unfolding. Only the portions of A1 external to this disulfide bond, residues 1238–1272 and 1458–1489, are directly subjected to the mechanical force transmitted through the spine of the VWF concatemer. However, multiple internal residues are noncovalently associated with folded external residues, and thus the stability of internal residues is expected to be affected by elongation of external residues by tensile force.

Here, to understand how binding of the VWF A1 domain to platelet GPIbα is regulated by its N- and C-terminal mucin-like linkers, we have measured the effects of both the glycan and polypeptide moieties of these linkers on stability, thermodynamics, and ligand-binding affinity and kinetics of A1. These linkers include the multiple serine and threonine residues that were found to be O-glycosylated during protein sequencing of native plasma VWF, including Thr-1248, Thr-1255, Thr-1256, and Ser-1263 in the N-linker (Titani et al., 1986). Mass spectrometry verified these assignments and showed a total of eight sialic acids over the four sites (Solecka et al., 2016). We show that both the N- and C-linkers affect affinity, that both their polypeptide and O-glycan moieties are important, and that the N-linker affects stability and hydrogen deuterium exchange (HDX) of A1 including in specific regions internal to its disulfide.

Results

Proteins

We compared the effect of the linkers N- and C-terminal to A1, which separate A1 from D3 and A2, respectively, and which contain O-glycosylation sites that have been chemically identified in native VWF (Solecka et al., 2016; Titani et al., 1986; Figure 1). Because the linkers are O-glycosylated in mammalian Expi293 cells and not in E. coli, we expressed A1 protein constructs in both cell types in order to test the effect in some assays of both polypeptide and O-glycan components of the linkers (Figure 1C). The N- and C-terminal boundaries of the constructs were based on allowing comparability to previous studies; references to the previous literature and a structural rationale for defining the minimal size of A1 as residues 1270–1463 and the maximal size of A1 plus linkers as residues 1240–1494 are described in Methods. To obtain A1 disulfide bond formation in bacteria, thioredoxin–A1 fusion proteins were expressed in E. coli SHuffle cells. The C-terminal boundary of the A1+N+C bacterial construct at Ser-1486 that was already on hand was extended by three residues in the mammalian A1+OGly+N+C construct to increase the likelihood of O-glycosylation of residue Ser-1486 and to include O-glycosylated residue Thr-1487 (Solecka et al., 2016; Titani et al., 1986). Constructs were purified using His tag affinity and by Superdex S200 size exclusion, and for the bacterial constructs, also with an intervening heparin affinity chromatography step. As expected based on O-glycosylation, the mammalian constructs migrated more slowly and showed increased heterogeneity in sodium dodecyl–sulfate polyacrylamide gel electrophoresis (SDS–PAGE) compared to the corresponding bacterial A1 constructs (Figure 1D). Furthermore, faster migration of all seven A1 constructs in nonreducing than reducing SDS–PAGE confirmed that all contained the long-range disulfide bond.

A1–GPIbα-binding kinetics and thermodynamics

Using biolayer interferometry (BLI), we measured A1 association to and dissociation from GPIbα immobilized through a biotinylated C-terminal avitag to BLI sensor tips. A1 preparations were subjected to S200 gel filtration the day of each measurement to remove aggregates. Global fits at all analyte concentrations to a single on- and off-rate for each A1 construct were good (Figure 2). In contrast, when gel filtration was omitted, data could not be fit to a single on- and off-rate. Differences among all seven A1 constructs in 150 mM NaCl showed that the linkers were of great importance in regulating A1 affinity for GPIbα (Figure 2). Comparisons to A1 short showed that nonglycosylated linkers in A1+N+C lowered affinity by 10-fold while the presence of glycosylated linkers in A1+OGly+N+C lowered affinity by 50-fold. The N-linker in A1+N lowered affinity by 4.3-fold while the C-linker in A1+C lowered affinity by 2.5-fold. Glycosylation of the N-linker was also more important than the C-linker; the affinity of A1+OGly+N was 10-fold lower than A1+N while the affinity of A1+OGly+C was only 2-fold lower than A1+C. O-Glycosylated linkers consistently lowered on-rates (Figure 2B). Overall, the results showed that both the polypeptide linker moiety and the O-glycan moiety of A1 linkers contributed to lowering A1 affinity for GPIbα, that the N-linker was more important than the C-linker in lowering affinity, and that the combined effect was very large at 50-fold.

Effect of A1 linkers, their glycosylation, and ionic strength on GPIbα-binding kinetics.

(A) Biolayer interferometry (BLI) traces in gray are shown for the indicated concentrations of the indicated A1 construct at the indicated NaCl concentration. In each panel, all A1 concentrations were fit globally to a one-site-binding model to obtain a single kon and koff. Thin black lines show fits. Fitting errors ranged from 1% to 5% of the fit values. Vertical dashed lines mark the beginning of the dissociation phase. (B) Kinetic constants and derived KD. Errors show difference from mean of n = 2 independent experiments.

To test the electrostatic contribution to binding affinity and on-rate, binding of A1 short and the three glycosylated constructs was additionally measured in 30 and 10 mM NaCl (Figure 2). Lowering NaCl concentration from 150 to 10 mM increased on-rates by an average of 17-fold for the four constructs. The effect of NaCl was predominantly on on-rate, since on- and off-rate together resulted in only a somewhat greater average increase in affinity of 23-fold (Figure 2B). Our results demonstrate that electrostatic complementarity between the binding interfaces in the basic A1 (pI 9.2) and acidic GPIbα (pI 5.5) domains (Huizinga et al., 2002) gives rise to electrostatic steering that speeds association kinetics.

To independently test affinity, and obtain binding thermodynamics, we used isothermal titration calorimetry (ITC). Although the GPIbα protein used in ITC lacked the avitag and its biotinylation, KD measurements in BLI and ITC were on average within 1.3-fold of one another and showed the same trends (Figure 3). A1 short showed a 39-fold increase in affinity compared to A1+OGly+N+C, the N-linker was more important than the C-linker, and the O-glycan moiety was more important than the polypeptide moiety in regulating affinity. Binding of all A1 protein constructs to GPIbα absorbed heats that showed the reaction was endothermic, that is, entropically driven (Figure 3B). The thermodynamics of A1–GPIbα binding showed enthalpy–entropy compensation (Liu and Guo, 2001; Starikov and Nordén, 2007): the higher the enthalpic contribution to binding, the higher the entropic cost, with all seven constructs following the same rank order for the enthalpic and entropic terms. Thus, the linkers affected both the enthalpic and entropic components of GPIbα binding.

Effect of A1 linkers and their glycosylation on GPIbα-binding thermodynamics.

(A) Panels show ITC traces of heat absorption after each injection (above) and fit of enthalpy to the binding isotherm (below) in 150 mM NaCl at 22°C. (B) Table summarizing the reaction Gibbs free energies (∆G) and the enthalpic (∆H) and entropic (TS) contributions. Errors show difference from mean of n = 2 independent experiments, except for A1+N+C which shows the fitting error from one experiment.

A1 stability and the thermodynamics of the intermediate state

The relative stabilities of A1 short and the three O-glycosylated A1 constructs to denaturation by heat or urea were measured by tryptophan fluorescence, whose emission maximum shifts to higher wavelength upon exposure to solvent during unfolding. The fluorescence intensity ratio (FIR), I350/I330, showed that constructs had three states, that is, exhibited two unfolding transitions (Figure 4A, B, D). Constructs containing the N-linker, that is A1+OGly+N+C and A1+OGly+N, were more stable to heat (Tm of 53.2 and 53.3°C for the first transition, respectively) than constructs lacking the N-linker, that is A1+OGly+C and A1 short (Tm of 49.9 and 50.7°C for the first transition, respectively) (Figure 4A, B). These results suggest that the N-linker stabilizes the native state of A1 relative to the intermediate state. Furthermore, all four constructs showed a second inflection point between 60° and 65° (second peak or shoulder of the derivative curve, Figure 4B). These results suggested two-step unfolding behavior and thus that A1 has three states, native (N), intermediate (I), and D (denatured) as previously observed (Auton et al., 2007).

Stability of A1 constructs to heat and urea and free energy difference between the native and intermediate A1 states.

(A–C) Effect of heating from 30 to 85°C on A1 stability and aggregation. (A) Unfolding was measured by excitation of tryptophan at 275 nm and measuring the ratio of the fluorescence intensity at 350 and 330 nm (fluorescence intensity ratio [FIR] = I350/I330). (B) The derivative of data in A (dFIR/dT) shows two unfolding transitions for each A1 construct and provides evidence for three states of A1: N, I, and U as described in Results. (C) Backscattered light intensity as a measure of A1 aggregation. (D) Unfolding of A1 by urea measured at 20°C by FIR as in panel A. Lines show fits to the three-state model of A1 unfolding as described in Methods and Results (see also Source Data). (E) Tabulation of the melting temperature of the native state (Tm) from panels A and B and the free energy difference between the native and intermediate states of A1 (ΔGI–N), the concentration of urea at which the I and N states are equally populated (Curea I–N), and m, from the fits in panel D.

A1 aggregation by heat was measured by backscattered light intensity. After the inflection point for the intermediate state, A1 short began to aggregate (Figure 4C compared to Figure 4B). A1+OGlyc+C required higher temperature for aggregation, which began to occur after the inflection point for the denatured state (Figure 4C compared to Figure 4B). Strikingly, the O-glycosylated N-linker completely protected A1+OGly+N and A1+OGly+N+C from aggregation at temperatures up to 85°C.

Denaturation in urea confirmed the existence of an intermediate state and was used to estimate the free energy of the intermediate state relative to the native state, that is, ΔGI–N. The A1 constructs showed two separate unfolding transitions, one between 1 and 4 M urea, and another above 5 M urea (Figure 4D). All constructs showed a plateau in between, or in the case of A1+OGly+C, followed a similar trend, but with a decrease in slope instead of a clear plateau. Similar three-state unfolding curves were previously seen for a construct similar to A1+N+C using circular dichroism (Auton et al., 2007). We fit the unfolding data to a three-state model. We assumed that all four constructs had the same m-value (a parameter that defines the cooperativity of transition) and that each state had a unique baseline FIR slope as a function of urea concentration that was shared among the four constructs (Methods). We thus determined the difference in energy between the I and N states (ΔGI–N), the concentration of urea at which the I and N states are equally populated (Curea I–N), and m, a value that is related to the size of folded domains (Ghosh and Dill, 2009; Santoro and Bolen, 1988; Figure 4E). The results showed that the O-glycosylated N-linker markedly increased the stability of the native state and thus the gap in energy between the native and intermediate state, with higher ΔGI–N values and transitions at higher urea concentrations for A1+OGly+N and A1+OGly+N+C than for A1 short and A1+OGly+C (Figure 4E).

Changes in HDX dynamics among A1 constructs correlate with affinity and stability differences

HDX MS measures solvent accessibility and H-bonding and reports on the dynamics and stability of proteins (Wales and Engen, 2006). The O-glycosylated A1 constructs and A1 short were allowed to exchange their backbone amide hydrogens for deuterium in D2O for varying periods of time, digested with pepsin, and the kinetics of deuterium exchange were measured for 108 A1 peptides (91% coverage, 5.8× redundancy, Figure 5 and Figure 5—figure supplements 15). HDX as a function of secondary structure and sequence position trended similarly for all four constructs studied (shown for A1 short in Figure 5A and for all in Figure 5—figure supplements 1 and 2). Figure 5A shows exchange for all 108 peptides at all time points, from 10 s to 4 hr. Over the A1 sequence from N- to C-terminal, the least exchange was seen for the β1-strand, the α1-helix, the β2- and β3-strands, the α3-helix, and the α6-helix, which all had at least one peptide with less than 20% deuteration at 4 hr. These slowly exchanging secondary structural elements neighbor one another and the GPIbα-binding site on A1 (Figure 5D).

Figure 5 with 5 supplements see all
The N-linker decreases A1 dynamics measured by hydrogen deuterium exchange (HDX).

(A) Relative deuterium exchange at all time points for A1 short as % of the available amide backbone H atoms in each peptide, colored according to the key. Figure 5—figure supplements 1 and 2 show data for all constructs. (B) Difference (Δ) in HDX at all time points of A1 short minus HDX for the other three constructs as indicated. In (A) and (B), residues with VWD type 2B mutations are shown below sequence numbers as red circles. (C) Full HDX kinetics for selected peptides. Figure 5—figure supplements 35 show data for all peptides. (D) Structure of A1 bound to GPIbα (PDB 1SQ0). A1 is colored rainbow from N- (blue) to C-terminus (red). Residues with VWD type 2B mutations are shown as Cα atom spheres. Labeled residues are shown in stick. GPIbα is shown in silver, from the β-finger to the β-switch.

Consistent differences in the amount of HDX were found in specific regions among the four A1 constructs (Figure 5B, C and Figure 5—figure supplements 35). The A1+OGly+N+C and A1+OGly+N constructs showed significantly less deuterium exchange (>0.7 Da) at almost all time points in peptides that were centered on the α1–β2 loop and the β3–α2 loop. These loops locate near to GPIbα (Figure 5B–D). We also observed reduced exchange in the α2–α3 loop, but the most meaningful difference was only seen at the 4 hr time point. The α1–β2 loop is nearby the long-range disulfide and in a region with multiple gain-of-function VWD type 2B mutations (Figure 5D). The β3–α2 loop is near to GPIbα and has four basic residues (labeled in Figure 5D) but they are not close enough to directly interact with acidic residues in the GPIbα leucine-rich repeats.

The HDX data showed no significant HDX differences between A1 short and A1+OGly+C. Further, the magnitude of HDX differences between A1 short and both A1+OGly+N+C and A1+OGly+N was very similar. However, while A1+OGly+N+C and A1+OGly+N were similar to one another in HDX, they both exchanged less deuterium than A1 short and A1+OGly+C, indicating that both A1 short and A1OGly+C were more dynamic than A1+OGly+N+C and A1+OGly+N. These HDX findings resembled the measurements of the stability of the native states relative to the intermediate states (Figure 4), which showed that A1+OGly+N+C and A1+OGly+N were similarly stable and more stable than A1 short and A1+OGly+ C, which were also similar in stability.

Discussion

Among domains in proteins, the A1 domain of VWF is highly unusual in being separated from neighboring domains by mucin-like linkers and being activated by tensile force transmitted through those linkers that originates from hydrodynamic force on VWF concatemers and bound platelets. Here, we have investigated the effects of these linkers on A1 thermodynamics and function in binding to its ligand, GPIbα. Previous studies have examined the effect of the N-terminal linker or both linkers on various types of quantitative or nonquantitative binding assays and thermodynamics and found inhibitory effects on GPIbα binding from adding the N-linker or both linkers to A1 (Arce et al., 2021; Auton et al., 2012; Deng et al., 2018; Deng et al., 2017; Interlandi et al., 2017; Ju et al., 2013; Miyata and Ruggeri, 1999; Nakayama et al., 2002; Tischer et al., 2017; Tischer et al., 2014). We found that both the N- and C-linker decreased A1 affinity for GPIbα and that for both the N- and C-linker, O-glycosylated linkers decreased affinity more than polypeptide linkers. For the C-linker, both the O-glycan and polypeptide moieties were similar in importance, giving ~2.2-fold decreases in affinity each, while for the N-linker, the O-glycan component was more influential (9.5-fold) than the polypeptide alone (4.3-fold). While adding O-glycosylated linkers on the N- and C-terminal ends of A1 lowered affinity substantially (51- and 6-fold, respectively), the presence of both O-glycosylated linkers lowered binding little compared to the presence of the N-linker alone (1.2-fold). The separate contributions of the polypeptide and glycan moieties of the A1 N- and C-linkers had not previously been examined; however, O-glycosylation of both linkers together was reported to decrease affinity from 3.9 to 32 µM (Tischer et al., 2020). While these reported affinities differ from ours by two- to threefold, they agree in terms of the affinity difference caused by O-glycosylation. Kinetic measurements showed that the effect of both the polypeptide and O-glycan components of the linkers was primarily on-rate driven. Overall, among the seven constructs, on-rates varied by up to 80-fold (A1 short versus A1+OGly+N+C) while off-rates varied up to 10-fold (A1+OGly+C versus A1+N).

Rapid binding is especially important for proteins that function in hemostasis; electrostatics can help orient proteins for proper binding before they contact one another, as previously demonstrated using the ionic strength dependence of binding of soluble GPIbα to VWF tethered to the wall of a flow chamber (Fu et al., 2017). Our kinetic measurements in 10 mM compared to 150 mM NaCl showed a 20-fold increase in affinity for A1+OGly+N+C largely driven by a sevenfold increase in kon, while VWF tethered to a flow chamber showed a 130-fold increase in affinity driven by an 11-fold increase in kon.

In 150 mM NaCl, the high-affinity state of tethered, mechanically tensioned VWF had an affinity of 80 nM for GPIbα, compared to the Kd values of 180 nM for A1 short and 9200 nM for A1+OGly+N+C measured here. Introduction of two VWD mutations into A1 was shown to increase affinity for GPIbα by 10-fold (Blenner et al., 2014). However, even though such a mutant would be expected to increase the affinity of the A1 short construct to a higher level than found for mechanically tensioned A1 in intact VWF, a crystal structure of A1 short with the two VWD mutations bound to GPIbα revealed only minor differences from a WT-GPIbα crystal structure (Blenner et al., 2014). Mechanical tension on the one hand and introduction of VWD mutations in the absence of tension and removal of the N- and C-linkers on the other hand may have different effects on the structure of the A1 domain, consistent with the hypothesis that A1 short with two VWD mutations bound to GPIbα did not reveal an A1 high-affinity state, despite an affinity for GPIbα comparable to that seen with mechanically tensioned VWF (Blenner et al., 2014).

We buttressed BLI affinity measurements here with ITC measurements of affinity and binding energies. While affinity measurements in the literature on essentially identical A1 fragments can differ up to 15-fold from those reported here (Huizinga et al., 2002), affinity measurements here on seven different A1 constructs by BLI and ITC were all within 1.6- to 1.04-fold of one another. We found it was necessary to remove aggregates from A1 preparations by gel filtration immediately prior to BLI measurements to avoid heterogeneous kinetics that required two on- and off-rates to fit, which appear in some of the literature on A1 binding to GPIbα (Arce et al., 2021; Tischer et al., 2017). In general, measurements in the liquid phase are more accurate than those in solid/liquid phases such as SPR and BLI, where a small amount of an avid aggregate can dominate binding kinetics.

We are unaware of previous thermodynamic measurements of A1 and GPIbα binding. ITC results show that the interaction of GPIbα with all seven A1 constructs was endothermic with a large entropic term and ΔG values that ranged from −8.2 to −14.7 kcal/mol at 22°C. Among all seven A1 constructs, the enthalpy and entropy terms were highly correlated, so that their rank orders were identical. Entropy-driven binding of proteins is usually attributed to the increase in water entropy when waters are released from hydrophobic binding interfaces and is often counterbalanced by an increase in enthalpy (Liu and Guo, 2001; Richards, 1977; Starikov and Nordén, 2007).

Our HDX studies revealed that the O-glycosylated N-linker has a significant effect in stabilizing the α1–β2, β3–α2, and the α2–α3 loops, as shown by less deuterium incorporation. Previous HDX studies have shown that destabilizing A1 by reducing and alkylating its long-range disulfide, introducing VWD type 2B mutations, or omitting N-linker residues 1238–1260 increase HDX in the same regions (Deng et al., 2018; Deng et al., 2017; Tischer et al., 2017; Tischer et al., 2014). The finding that stands out from previous reports, which did not study the C-linker individually in HDX, is that in contrast to the N-linker, the C-linker did not significantly affect HDX in any of these regions.

We directly measured the stability of A1 by measuring its unfolding by heat and urea. We found evidence for three states by both techniques, consistent with a previous finding of an intermediate state that is hypothesized to represent a state with high affinity for GPIbα (Auton et al., 2010a; Tischer et al., 2017; Tischer et al., 2014). A plateau between the first and second transition clearly demonstrated the existence of the intermediate state for three of our four constructs, in contrast to the decrease in slope seen previously. This difference was likely due to measurement here of the fluorescence of the single tryptophan residue, Trp-1313 in A1, as opposed to circular dichroism, which measures contributions from all A1 secondary structures. Trp-1313 is in the α1–β2 loop that has HDX dynamics that are sensitive to the presence of the N-linker, locates near the GPIbα β-finger, is mutated in VWD type 2B, and although it may report denaturation of only a portion of A1, Trp-1313 discriminates between the native, intermediate, and denatured states of A1 similarly to circular dichroism.

Heat and urea denaturation had similar linker-dependent effects on stability of A1. A1+OGly+N and A1+OGly+N+C were each ~3°C more stable than A1 short and A1+OGly+C. Likewise, A1+OGly+N and A1+OGly+N+C each had native (N) states that were relatively more stable than their intermediate (I) states, by ~2 kcal/mol, than found for A1 short and A1+OGly+C. Additionally, the N-linker was sufficient to completely protect A1 against heat-induced aggregation up to 85°C. Since physiological tensile force on A1 may induce the I state, which could make A1 and with it an entire VWF concatemer susceptible to aggregation, one of the several important physiological roles of the N-linker and its O-glycosylation may be to protect VWF from aggregation. The neighboring A2 domain in VWF can be completely and reversibly unfolded by tensile force and can also be reversibly heated to 90°C without aggregation in the presence of its O-linked N-linker (the C-linker of A1) and its C-linker (Xu and Springer, 2013; Zhang et al., 2009), which may have a similar protective role physiologically to prevent aggregation during physiological partial (A1) and full (A2) unfolding.

Our results on the effects of the O-linked N- and C-linkers on affinity and stability differ. The C-linker decreased A1 affinity for GPIbα in two different measures, BLI and ITC, but had no effect on A1 stability as measured by HDX, melting temperature, and the relative stability of the A1 N and I states in urea denaturation experiments. These results suggest that the C-linker decreases affinity by direct mechanisms such as by steric or electrostatic repulsion. O-Glycans are decorated with sialic acid, which would repel negatively charged GPIbα. In contrast to the C-linker, the N-linker not only decreased affinity but also increased A1 stability as shown by significantly decreased HDX, increased melting temperature, and increased stability of the N state relative to the I state. The ~2 kcal/mol difference in ΔG between the N and I state for A1+OGly+N compared to A1 short corresponds to a ~30-fold difference in population of the I state. Thus, if the I state was the high-affinity state, much of the 40-fold difference in affinity between the I and N states could be accounted for by the shift in equilibrium toward the native state caused by the N-linker, and the remainder could be caused by steric and electrostatic repulsion of GPIbα. In VWF concatemers, the ΔG required for activation by mechanical tension of the high-affinity state was measured as 1.9 kcal/mol (Fu et al., 2017), considerably lower than the difference in energy between the I and N states measured for A1+OGly+N+C here of 4.7 kcal/mol. These differences suggest that the directional nature of energy input by mechanical tension allows for more efficient activation of A1 than urea denaturation and may induce a distinct and less disordered intermediate state.

Three different mechanisms have been proposed for how the linkers affect A1 binding to GPIbα. The proposal that the N- and C-linker associate with one another to form an autoinhibitory module (Deng et al., 2018; Deng et al., 2017) is not supported by our data that show that the N-linker alone is sufficient to inhibit binding of A1 to GPIbα. Furthermore, we find that the N- and C-linker are not synergistic and to the contrary, have much less than additive effects on lowering binding affinity. It also has been proposed that residue Asp-1261 in the N-linker interacts with basic residues in the A1 domain to stabilize a low-affinity state of A1 (Interlandi et al., 2017). While our results do not bear on this mechanism, they do demonstrate that the O-linked glycans within the N-linker are sufficient to strongly inhibit binding of A1 to GPIbα. Finally, Auton et al. have described an intermediate state of A1 that is proposed to correspond to an activated high-affinity state of A1 (Auton et al., 2007; Auton et al., 2012; Tischer et al., 2014). Our data are consistent with and provide further evidence in support of this hypothesis.

The association of the N-linker but not the C-linker with A1 stability is consistent with structural analysis of A1 complexes with GPIbα that demonstrate that their association strains the N-terminal portion of A1 (reviewed in Blenner et al., 2014). The main site of GPIbα association is at the β-switch in its C-terminal cap (Figure 5D). The flexible β-finger in the N-terminal cap of GPIbα associates with A1 over a smaller region, including the α1–β2 loop, with few specific contacts. Furthermore, these contacts vary among GPIbα complexes. Strain in A1 upon binding GPIbα causes shifts in some residues in the α1–β2 loop and Cys-1272 away from GPIbα. Furthermore, the hydrogen bond between the backbone of Cys-1272 and the backbone of Arg-1308 (a residue often mutated in VWD), which is conserved in all structures of the isolated A1 domain, is lost in all A1–GPIbα complex structures. This loss admits a water molecule to the hydrophobic core of A1. These structural observations begin to define a pathway for conformational communication between the N-linker and regions of A1 that change upon GPIbα binding, including the α1–β2 loop, which showed increased HDX dynamics in the absence of the N-linker. Crystal structures of A1 diverge prior to residue 1270, suggesting that the O-glycosylated N-linker of A1 (defined in our experiments as residues 1238–1263) is natively disordered. It is now well established that natively disordered regions in proteins frequently have important roles in regulating protein allostery (Motlagh et al., 2014). We speculate that in the intermediate state, A1 would reshape to enable closer approach to the concave, acidic surface of GPIbα of more A1 basic residues, such as those in its β3–α2 loop shown in stick in Figure 5D.

In summary, our results show that the A1 N-linker, but not the C-linker, increases the stability of the A1 native state relative to the intermediate state, thus regulating conformational change. The previous paragraph discussed structural evidence for a pathway for communication from the N-linker to regions that are known to change structurally upon GPIbα binding and also have backbone motions that are influenced by the N-linker as shown by HDX. Additionally, our results suggest that the C-linker influences affinity for GPIbα by steric and possibly charge repulsion. The N-linker may repel GPIbα as well, in addition to regulating A1 conformation.

Materials and methods

Boundaries of A1 constructs

Request a detailed protocol

Boundaries for truncations of A1 used here were taken from what had been used in the literature and structural considerations. Constructs beginning here at residues 1238 and 1264 had previously been used in related papers in Auton et al., 2007; Auton et al., 2012; Auton et al., 2010a; Auton et al., 2010b and (Blenner et al., 2014), respectively. Similarly, a construct ending at residue 1469 was previously used in Blenner et al., 2014. Longer constructs used here ending at 1486 or 1489 ended at positions intermediate between longer constructs ending at 1480 (Madabhushi et al., 2014; Tischer et al., 2020) and 1493 (Deng et al., 2018). Twelve examples of A1 crystallized alone or in complexes (PDB codes 1m10, 1sq0, 1u0n, 1auq, 1u0o, 1ijk, 1ijb, 1fns, 1uex, 3hxo, 3hxq, and 4c2a) were superimposed to define boundaries between residues that were largely invariant in structure (part of the A1 domain) or variable among structures (part of linkers). Superimposition shows that the A1 backbone takes the same path beginning at residue F1270 and ending at residue E1463; the A1 disulfide links C1272 and C1458. The paths of two D3 assembly monomers (PDB code 6n29) are similar up to residue E1239 and the paths of five A2 monomers in PDB codes 3gxb and 3zqk are similar beginning at residues 1495. Thus, the minimum length of the A1 domain might be considered to be between residues 1270 and 1463 and the A1 linkers might be considered to extend up to the boundaries of the D3 and A2 domains, that is from residue 1240 to 1494. We are not aware of tests of the feasibility of obtaining good expression of constructs of these lengths and chose lengths that would allow comparison to previous results as described above.

Glycosylated proteins

Request a detailed protocol

Human VWF A1 domain constructs as shown in Figure 1C (pre-pro-VWF amino acid residue numbering), beginning with residue 1238 or 1264 and ending with residue 1469 or 1489 and followed with a C-terminal 6xHis tag were cloned into a bicistronic IRES-GFP expression vector, ET8 (Zhou and Springer, 2014) and transfected into Expi293F cells using lipofectamine reagent (Invitrogen). Stable transfectants were selected in the presence of 500 μg/ml G418 (Geneticin) and subjected to two rounds of fluorescence-activated cell sorting. In each sort, the 8% most fluorescent cells were expanded in Expi293 expression medium. Culture supernatants were harvested after 5 days, supplemented with 20 mM HEPES (N-2-hydroxyethylpiperazine-N'-2-ethanesulfonic acid), pH 7.4 and 0.2 mM NiCl2 (final concentration), and loaded onto a Ni-NTA agarose column by gravity (2 ml Ni-NTA agarose beads per 500 ml of culture supernatant). The column was pre-equilibrated in Washing Buffer (20 mM HEPES, pH 7.4, 150 mM NaCl, 0.2 mM NiCl2). After loading the supernatant, the column was washed with 5 column volumes of Washing Buffer followed by 10 column volumes of 20 mM HEPES, 1 M NaCl, 16 mM imidazole, 0.2 mM NiCl2. Protein was eluted with 20 mM HEPES, 150 mM NaCl, and 300 mM imidazole, and further purified by Superdex 200 (GE Healthcare) size-exclusion chromatography in 20 mM HEPES, 150 mM NaCl, pH 7.5.

Wild-type GPIbα, residues 1–290 of the mature protein with a C-terminal His6 tag, with or without an intervening Avi-tag sequence (GLNDIFEAQKIEWHE), was purified from culture supernatant of HEK293 stable transfectants (Blenner et al., 2014) by Ni-NTA affinity chromatography and Superdex 200 size exclusion as above.

Nonglycosylated proteins

Request a detailed protocol

The cDNA sequence encoding A1 beginning with residue 1238 or 1264 and ending with residue 1469 or 1486 (Figure 1C) was cloned into the pET32a vector with a thioredoxin fusion protein (LaVallie et al., 1993) and a TEV cleavage site at its N-terminus and a 6xHis tag at its C-terminus. Proteins were expressed in SHuffle cells (Lobstein et al., 2012), which were grown at 37°C in LB medium; at an A600 of 0.8, 1 mM IPTG was added and the temperature was shifted to 25°C. After 20 hr, cells were collected by centrifugation and lysed in a French press at 4°C in 20 mM HEPES (pH 7.4), 150 mM NaCl. The lysate was centrifuged at 17,000 × g at 4°C for 40 min. Ni-NTA agarose beads were added to the supernatant containing soluble A1 and shaken overnight at 4°C in 40 ml Falcon tubes (1.5 ml beads/40 ml supernatant). Beads were centrifuged at 2500 RPM for 5 min and resuspended in 10 ml buffer (20 mM HEPES, 150 mM NaCl, pH 7.4) per tube and vortexed for 2 min (first washing step). After two more washing steps, five cleaning steps were performed following the same procedure but with 20 mM HEPES, 1 M NaCl, and 16 mM imidazole. The beads were then packed into a column and eluted with 20 mM HEPES, 150 mM NaCl, 300 mM imidazole by gravity. The eluate was mixed with Tobacco Etch Virus (TEV) protease (A1:TEV mass ratio of 10:1) and dialysed against 20 mM HEPES, 150 mM NaCl overnight at 4°C. The thioredoxin-cleaved A1 was centrifuged for 10 min at 13,000 rpm to remove precipitate. The supernatant was further purified by heparin column (HiTrap Heparin HP, GE Healthcare), washed with 20 mM HEPES, 350 mM NaCl, and eluted with 20 mM HEPES, 800 mM NaCl. Finally, the protein was subjected to Superdex 200 size-exclusion chromatography as above.

Storage and repurification

Request a detailed protocol

Protein concentrations were determined from their A280 using the extinction coefficient calculated by ProtParam tool at the Expasy website. All proteins were stored at −80°C and subjected to a second round of size-exclusion chromatography on the day of measurements. This was essential to prevent protein aggregates from contributing to affinity and kinetics measurements, that is, to obtain data that can be reliably fitted with a 1:1 binding model.

Kinetic binding measurements

Request a detailed protocol

The GPIbα used for kinetic measurements contained an N-terminal Avi-tag which was biotinylated using a BirA biotin-protein ligase kit (Cat #BirA500, Avidity, Aurora, CO). BLI used Octet RED384 instrument and software (ForteBio). Streptavidin-functionalized sensors were dipped in biotinylated GPIbα (1 μM) for 60 s and quenched with 10 μg/ml biotinyl-lysine for 40 s. Binding and dissociation responses were recorded at varying A1 concentrations in 20 mM HEPES, pH 7.5, 0.02% Tween-20, and either 10, 30, or 150 mM NaCl as indicated. Response in a buffer-only reference well was subtracted (internal reference). The experiment was also carried out using the same concentrations of A1 and the internal reference but with sensors that had no absorbed GPIbα (parallel reference). Parallel reference responses were subtracted from the experimental responses with Octet software. Data at different A1 concentrations were fit globally to single kon and koff values for each experimental condition and plotted using Originlab (Origin, Northampton, MA). The ionic strength of 20 mM HEPES, pH 7.5 was calculated using https://www.liverpool.ac.uk/pfg/Research/Tools/BuffferCalc/Buffer.html.

Isothermal calorimetry

Request a detailed protocol

Proteins were dialyzed overnight against 150 mM NaCl, 20 mM Tris–HCl, pH 7.5, degassed, and centrifuged at 20,000 × g for 10 min. GPIbα (250 μM except 550 μM with A1+OGly+N+C) was titrated into an A1 protein solution (22 μM except 50 μM with A1+OGly+N+C) in a MicroCal iTC200 (GE Healthcare Life Sciences). A priming injection of 0.4 μl (not included in data analysis) was followed by 2 μl injections every 180 s. Data averaged over 2 s windows were analyzed using Originlab 7.

Stability to denaturant and heat

Request a detailed protocol

For urea denaturation, A1 (0.25–0.5 mg/ml) in 10 μl of 20 mM HEPES, pH 7.5, 150 mM NaCl was mixed with 10 μl of 200 mM HEPES, pH 7.5, 1.5 M NaCl, ×μl of 7.5 M urea, and (80-x) μl of water, where x gave final urea concentrations of 0.135–6 M in 24 concentration steps in 100 μl, and equilibrated for 24 hr at 20°C. Intrinsic Trp fluorescence was measured at 330 and 350 nm upon excitation at 275 nm with a Prometheus NT48 (NanoTemper) with automatic adjustment of excitation power.

FIR, emission intensities at 350 nm over 330 nm, was fit as a function of urea concentration to a three-state transition model using the linear extrapolation method (Santoro and Bolen, 1988). FIR at each urea concentration (Curea) is modeled as the population weighted contribution from the native state (N), intermediate state (I), and the denatured state (D): FIR (Curea)=FIRNPN+FIRIPI+FIRDPD where FIRN , FIRI, and FIRD represent the linear baseline for the N, I, and D states, respectively, which can be expressed as FIRj=bj+sjCUrea , with bj and sj as intercept and slope of the j state. PN , PI, and PD represent the population of N, I, and D states at each urea concentration, respectively. The population of state j (Pj) in the system can be related to the free energy of each state (N, I, D) based on the Boltzmann distribution, Pj=Exp(GjR.T)Exp(GNR.T)+Exp(GIR.T)+Exp(GDR.T) , where R is the gas constant, and T is absolute temperature of the system. With the N state as reference state (ΔGN=0), the free energy of the I and D states at each urea concentration can be expressed as ΔGI(CUrea)=ΔGI0+mI.CUrea and ΔGD(CUrea)=ΔGD0+mD.CUrea, respectively, where ΔGI0 and ΔGU0 are the free energy of the I and D state in absence of urea, and mI and mD are the cooperativity of transition from the N to I state, and from N to D state, respectively. For a cooperative two-state transition induced by a particular denaturant, the m-value is proportional to the surface area exposed upon protein denaturation and is generally proportional to the protein size. The linear relationship between m-values and globular protein chain lengths have been reported for denaturant-induced protein unfolding (Ghosh and Dill, 2009). Thus, the final fitting function for the three-state transition is as follows, where b and s are baseline intercept and slope values, respectively, for N or I states:

(1) FIR(Curea)=(bN+sNCUrea)+(bI+sICUrea)Exp(ΔGI0+mICUreaRT)+(bU+sUCUrea)Exp(ΔGD0+mDCUreaRT)1+Exp(ΔGI0+mICUreaRT)+Exp(ΔGD0+mDCUreaRT)

Nonlinear least square fit to Equation 1 was employed to fit the urea denaturation data for four A1 constructs. The data were globally fit with shared sN , sI , sD , mI, and mD , and individual bN , bI , bD , ΔGI0, and ΔGD0 for each construct. The rationale for fitting the baseline slope for each state to a shared value between the four constructs is that the N- and C-terminal linker does not contain any tryptophan. The rationale for the four constructs to share the same m-values for the I and D states is that the four constructs only differ in the O-glycosylated linkers which are not folded.

For heat denaturation, A1 (0.25–0.5 mg/ml) was in 20 mM HEPES, pH 7.5, 150 mM NaCl. FIR was measured in the Prometheus NT48 while temperature was ramped from 30 to 85°C at a rate of 1 °C/min using an excitation power of 50%. Additionally, the intensity of backscattered light was recorded as a measure of aggregation.

HDX mass spectrometry

Request a detailed protocol

HDX experiments were essentially as reported previously (Iacob et al., 2013; Le et al., 2018). Comprehensive experimental details and parameters are provided in Table 1, in the recommended (Masson et al., 2019) tabular format. All HDX MS data have been deposited to the ProteomeXchange Consortium via the PRIDE (Perez-Riverol et al., 2019) partner repository with dataset identifier PXD029942.

Table 1
Hydrogen deuterium exchange (HDX) MS data summary and list of experimental parameters.
DatasetA1 shortA1NA1CA1NC
HDX reaction details*15-Fold dilution with labeling buffer at 20°C, final D2O level = 93.3%, pHread = 7.10; equal volume quench buffer at 0°C, pHread = 2.5
HDX time course10 s, 1 m, 10 m, 1 hr, 4 hr
HDX controls3 undeuterated for each condition
Back exchange30–35%
Number of peptides102103101108
Sequence coverage93.2%92.5%91.8%90.9%
Avg. peptide length (a.a.)12.312.212.312.3
Redundancy5.45.144.895.77
Replicates2 technical for each condition
HDXrepeatability±0.15 relative Da
  1. *

    Labeling buffer: 20 mM HEPES, 150 mM NaCl, 99% D2O, pD 7.5; quench buffer: 300 mM sodium phosphate, 0.25 M Tris (2-carboxyethyl)phosphine hydrochloride (TCEP-HCl), H2O pH 2.5.

  2. All reported values are the average relative deuterium level as given by the DynamX software, which in some cases is an average of more than one charge state, across all peptides in both technical replicates. No statistical tests were applied to the HDX MS measurements.

Briefly, A1 short and each O-glycosylated protein construct (28 μM in 20 mM HEPES, 150 mM NaCl) was diluted 15-fold into 20 mM HEPES, 150 mM NaCl, 99% D2O (pD 7.5) at room temperature. At deuterium exchange time points from 10 s to 240 min, an aliquot was quenched by adjusting the pH to 2.5 with an equal volume of 300 mM sodium phosphate, 0.25 M tris (2-carboxyethyl)phosphine hydrochloride (TCEP-HCl), H2O. Samples were analyzed as previously described (Iacob et al., 2013; Wales et al., 2008). They were digested offline with 10 mg/ml pepsin in water for 5 min on ice, then injected into a custom Waters nanoACQUITY UPLC HDX Manager. All mass spectra were acquired using a Waters Synapt G2-Si HDMSE mass spectrometer. Comparison experiments were done under identical experimental conditions such that deuterium levels were not corrected for back exchange and are therefore reported as relative (Wales and Engen, 2006). All experiments were performed in duplicate. The error of measuring the mass of each peptide was ±0.15 Da in this experimental setup. The peptides were identified using PLGS 3.0.1 software and the HDX MS data were processed using DynamX 3.0 (Waters Corp., USA). Peptides common to A1+OGly+N+C and the shorter constructs were followed with HDX, with an overall sequence coverage of 90.9%. In total, 108 peptic peptides were followed with HDX uptake plots (Figure 5—figure supplements 15).

Data availability

Figure 4—source data 1–4 contain the numerical data used to generate the figures. All HDX MS data have been deposited to the ProteomeXchange Consortium via the PRIDE (Perez-Riverol et al., 2019) partner repository with dataset identifier PXD029942.

The following data sets were generated

References

    1. Clemetson KJ
    (1983)
    Platelet membrane glycoprotein I: structure and function The domain of glycoprotein I involved in the von Willebrand receptor
    Blood Cells 9:319–329.
    1. Sugimoto M
    2. Dent J
    3. McClintock R
    4. Ware J
    5. Ruggeri ZM
    (1993)
    Analysis of structure-function relationships in the platelet membrane glycoprotein Ib-binding domain of von Willebrand’s factor by expression of deletion mutants
    The Journal of Biological Chemistry 268:12185–12192.
    1. Woollett GR
    2. Williams AF
    3. Shotton DM
    (1985)
    Visualisation by low-angle shadowing of the leucocyte-common antigen A major cell surface glycoprotein of lymphocytes
    The EMBO Journal 4:2827–2830.

Decision letter

  1. Yelena Ginzburg
    Reviewing Editor; Icahn School of Medicine at Mount Sinai, United States
  2. Mone Zaidi
    Senior Editor; Icahn School of Medicine at Mount Sinai, United States
  3. P Clint Spiegel
    Reviewer; Western Washington University, United States

In the interests of transparency, eLife publishes the most substantive revision requests and the accompanying author responses.

Decision letter after peer review:

Thank you for submitting your article "Von Willebrand factor A1 domain affinity for GPIbα and stability are differentially regulated by its O-glycosylated N-linker and C-linker" for consideration by eLife. Your article has been reviewed by 3 peer reviewers, and the evaluation has been overseen by a Reviewing Editor and Mone Zaidi as the Senior Editor. The following individual involved in review of your submission has agreed to reveal their identity: P. Clint Spiegel (Reviewer #2).

The reviewers have discussed their reviews with one another, and the Reviewing Editor has drafted this to help you prepare a revised submission.

Essential revisions:

In this paper, the authors use several biophysical tools, including biolayer interferometry, isothermal calorimetry, protein denaturation, and hydrogen-deuterium exchange mass spectrometry, to systemically dissect the influence of the N- and C-terminal polypeptide linkers and O-glycan moieties of the VWF A1 domain on A1's activation state and binding affinity to platelet receptor GPIbα.

Through a series of elegant experiments, the authors concluded that the N-terminal linker together with its O-linked glycans stabilize the A1 in an inactive, native conformation. Removal of the N-linker or its O-linked glycans destabilizes the A1 structure and may activate the A1 into an intermediate, higher affinity state for GPIbα binding. The authors also concluded that the C-terminal linker and its O-linked glycans do not contribute to stabilizing the A1 structure, though they prevent GPIbα binding through steric hindrance and electrostatic repulsion. The data and the figure presented in the manuscript are overall clear, and the conclusions are consistently supported by the data from multiple experimental techniques. The reported results are highly interesting for the hematology and biophysics communities.

1) Define the A1 domain in terms of the sequence that delimits it and the sequences that define the N- and C- linkers and provide justification for the selection.

2) Show the standard deviation of the fit instead of showing the average of two experiments in figure 2 and provide additional description for the methods used to purify protein constructs in the gel images as the bands appear differentially pure, in different amounts in figure 1.

Reviewer #1 (Recommendations for the authors):

In this paper, the authors use several biophysical tools, including biolayer interferometry, isothermal calorimetry, protein denaturation, and hydrogen-deuterium exchange mass spectrometry, to systemically dissect the influence of the N- and C-terminal polypeptide linkers and O-glycan moieties of the VWF A1 domain on A1's activation state and binding affinity to platelet receptor GPIbα.

Through a series of elegant experiments, the authors concluded that the N-terminal linker together with its O-linked glycans stabilize the A1 in an inactive, native conformation. Removal of the N-linker or its O-linked glycans destabilizes the A1 structure and may activate the A1 into an intermediate, higher affinity state for GPIbα binding. The authors also concluded that the C-terminal linker and its O-linked glycans do not contribute to stabilizing the A1 structure, though they prevent GPIbα binding through steric hindrance and electrostatic repulsion. The data and the figure presented in the manuscript are overall clear, and the conclusions are consistently supported by the data from multiple experimental techniques. The reported results are highly interesting for the hematology and biophysics communities

I have only a few suggestions for revision and improvement, detailed below.

1) The length of the long C-linker is inconsistent among the constructs. While the constructs expressed from E. coli end at residue 1486, the mammalian constructs end at residue 1489. No further justification was given by the authors. Moreover, since the A2 domain starts at residue 1494, why didn't the authors include residues 1490-1493?

2) In regard to Figure 2, some of the global fits are not of high quality (e.g., A1+C with 150 mM NaCl). The authors should show the standard deviation of the fit instead of showing the average of two experiments.

3) It would be worth discussing the difference between the intermediate state identified by the thermal/chemical denaturation experiments and the more physiological mechanically tensioned structure. The intermediate state may or may not be the same high-affinity conformation induced by tensile pulling, since heat/urea destabilizes the entire A1 structure and tensile force/tension is not expected to spread into the structures protected by the long-ranged disulfide.

4. Lines 223-224: It would be nice to specify the exact constructs and conditions that account for the 80-fold and 10-fold differences.

Reviewer #2 (Recommendations for the authors):

Bonazza, Klaus et al. investigate the molecular details for how von Willebrand factor (VWF) binds to Glycoprotein Ib (GPIb) on platelets. Previous studies have demonstrated that VWF binds to GPIb only under conditions of high tensile force through an unfolding event that occurs for the A1 domain of VWF. This study emphasizes the role(s) of the unstructured linkers, along with their associated glycosylation states, that are N- and C-terminal to the A1 domain. Their data demonstrate that both linkers decrease the affinity of VWF A1 for GPIb; the N-terminal linker decreases the affinity more strongly whilst stabilizing the A1 fold, and glycosylation of the linkers contributes more to the decreased affinity than the peptide component. With biolayer interferometry (BLI), they present association and dissociation curves for this intermolecular interaction that fit well to a one-site binding model to calculate on- and off-rates with high confidence, showing a decrease in affinity for various constructs containing the N- and/or C- terminal linkers (with and without glycosylation). Additionally, they use BLI to show that the interaction is largely electrostatic by calculating higher affinities at lower ionic strength conditions. They further support the BLI data with measuring equilibrium binding thermodynamics with isothermal titration calorimetry (ITC). This experimental approach allows them to demonstrate similar affinity values to the BLI results, but also provides a deeper understanding for the thermodynamic basis for the VWF/GPIb interaction; the binding event is endothermic (entropically driven). In order to investigate how the N- and C-linkers affect protein stability, they employed both thermal and chemical denaturation measured with intrinsic tryptophan fluorescence. These data suggest that the A1 domain adopts an intermediate, partially unfolded state, which has differential stability dependent on the presence of the linkers; the N-terminal linker stabilizes the 'native' state, thus restricting the partially unfolded state that is hypothesized to be the GPIb binding state. Following back scattering intensity as a function of temperature in the same experiment, they also show that the N-linker prevents aggregation of the A1 domain. Lastly, the authors employ H/D exchange measured by mass spectrometry (HDX-MS) to show that the presence of the N-linker specifically increases protection patterns (decreased H/D exchange) for regions of the A1 domain that are adjacent to the GPIb binding interface and the internal disulfide bond. Taken together, what emerges from this study is a model for VWF A1/GPIb association that is regulated by glycosylated linkers whereby the N-linker stabilizes a 'native' conformation, which is only perturbed to an intermediate, partially unfolded state under high tensile force, leading to the activation of a high affinity binding state for GPIb. The data presented in this study are clearly communicated and appropriately controlled and interpreted.

Some perceived weaknesses are as follows:

1. The authors appropriately report a gel image of their purified protein constructs. The bands on the gel appear differentially pure, in different amounts, and there is no description in the manuscript about the relative 'diffuse' nature of the bands for glycosylated protein constructs.

2. The study describes the rigorous treatment of their thermal and chemical denaturation curves, and they briefly address the difference between measuring protein unfolding via intrinsic tryptophan fluorescence vs. circular dichroism. None of the curves presented illustrate error bars, but they do provide a measure of error in the associated table. Their measurements in this experiment are exclusively dependent on the environment of a single tryptophan residue in the A1 domain structure; their interpretations of the unfolding curves are warranted, but not without the caveat that it is only measuring the local environment of the tryptophan and not the entire protein structure.

My recommendations for the authors concern the gel image in figure 1 and the Materials and methods section. For figure 1, the protein bands are underwhelming. The reader would have more confidence in the study if this gel demonstrated a more uniform concentration of the protein constructs used in the study. Some of the bands are quite large and dominate the figure while some are comparatively weaker (+N+C in particular). For the Materials and methods section, I question the accuracy of a couple statements: it is stated that you have NiCl in the purification buffers. Ni is divalent, so this should be more accurately reported as NiCl2. In the section on non-glycosylated proteins, the authors state that SHuffle cells are grown at 37C, but these cells are typically grown at 30C due to their thermolability (if this is correct, that is fine, but it diverges from the manufacturer's protocol). Lastly, please define for Equation 1 what b and s are.

Reviewer #3 (Recommendations for the authors):

The authors examined the effect of N- and C-linkers of the von Willebrand factor (VWF) A1 domain on binding to platelet glycoprotein Ibalpha (GPIbalpha). They made recombinant A1 proteins with and without the linkers and with and without O-glycans. They measured the interaction of the proteins by biolayer interferometry (BLI), isothermal titration calorimetry (ITC). Additionally, they measured the thermal and urea-induced denaturation and the hydrogen deuterium exchange (HDX) kinetics of the proteins. They found that both the N-and C-linkers lower A1 affinity, the effects of which are increased by O-glycosylation. The N-linker, but not the C-linker, stabilized the A1 domain and affected HDX exchange. They conclude that the A1 domain is regulated allosterically by the N-linker.

Overall, the conclusions of the study appear to be supported by carefully conceived and conducted experiments. The manuscript is well-written and nicely illustrated. The dissociation constants measured by BLI differ significantly from those measured by ITC, indicating systematic errors in one or both of the methods. However, the linker and glycan effects are qualitatively the same. Thus, the authors' conclusions are supported by the results of two independent physical methods.

The findings support and extend previously published findings by several groups that N- and C-linkers of the A1 domain regulate its interaction with GPIbalpha. However, the authors choose to interpret their results as being inconsistent with the proposal that N-linker and C-linker associate with one another to form an autoinhibitory module because they find that the N-linker alone is sufficient to inhibit binding of A1 to GPIbalpha. However, it is possible that the C-linker also contributes to the regulation of A1 binding to GPIbalpha, which, in fact, is supported by the authors' own the BLI and ITC results.

It would be useful for the authors to define the A1 domain in terms of the sequence that delimits it and the sequences that define the N- and C- linkers. Is the A1 delimited by residues 1272 and 1458 or by residues 1264 – 1469 (Figure 1C)?

The organization at the end of the Introduction perhaps could be improved. The last paragraph begins with "Here, to understand … we have made fundamental measurements …" preparing the reader for what the authors did. However, the authors proceed in the next few sentences to continue to cite observations from the literature.

The authors examined the effect of N- and C-linkers of the von Willebrand factor (VWF) A1 domain on binding to platelet glycoprotein Ibalpha (GPIbalpha). They made recombinant A1 proteins with and without the linkers and with and without O-glycans. They measured the interaction of the proteins by biolayer interferometry (BLI), isothermal titration calorimetry (ITC). Additionally, they measured the thermal and urea-induced denaturation and the hydrogen deuterium exchange (HDX) kinetics of the proteins. They found that both the N-and C-linkers lower A1 affinity, the effects of which are increased by O-glycosylation. The N-linker, but not the C-linker, stabilized the A1 domain and affected HDX exchange. They conclude that the A1 domain is regulated allosterically by the N-linker.

Overall, the conclusions of the study appear to be supported by carefully conceived and conducted experiments. The manuscript is well-written and nicely illustrated. The dissociation constants measured by BLI differ significantly from those measured by ITC, indicating systematic errors in one or both of the methods. However, the linker and glycan effects are qualitatively the same. Thus, the authors' conclusions are supported by the results of two independent physical methods.

The findings support and extend previously published findings by several groups that N- and C-linkers of the A1 domain regulate its interaction with GPIbalpha. However, the authors choose to interpret their results as being inconsistent with the proposal that N-linker and C-linker associate with one another to form an autoinhibitory module because they find that the N-linker alone is sufficient to inhibit binding of A1 to GPIbalpha. However, it is possible that the C-linker also contributes to the regulation of A1 binding to GPIbalpha, which, in fact, is supported by the authors' own the BLI and ITC results.

It would be useful for the authors to define the A1 domain in terms of the sequence that delimits it and the sequences that define the N- and C- linkers. Is the A1 delimited by residues 1272 and 1458 or by residues 1264 – 1469 (Figure 1C)?

The organization at the end of the Introduction perhaps could be improved. The last paragraph begins with "Here, to understand … we have made fundamental measurements …" preparing the reader for what the authors did. However, the authors proceed in the next few sentences to continue to cite observations from the literature.

https://doi.org/10.7554/eLife.75760.sa1

Author response

Essential revisions:

In this paper, the authors use several biophysical tools, including biolayer interferometry, isothermal calorimetry, protein denaturation, and hydrogen-deuterium exchange mass spectrometry, to systemically dissect the influence of the N- and C-terminal polypeptide linkers and O-glycan moieties of the VWF A1 domain on A1's activation state and binding affinity to platelet receptor GPIbα.

Through a series of elegant experiments, the authors concluded that the N-terminal linker together with its O-linked glycans stabilize the A1 in an inactive, native conformation. Removal of the N-linker or its O-linked glycans destabilizes the A1 structure and may activate the A1 into an intermediate, higher affinity state for GPIbα binding. The authors also concluded that the C-terminal linker and its O-linked glycans do not contribute to stabilizing the A1 structure, though they prevent GPIbα binding through steric hindrance and electrostatic repulsion. The data and the figure presented in the manuscript are overall clear, and the conclusions are consistently supported by the data from multiple experimental techniques. The reported results are highly interesting for the hematology and biophysics communities.

1) Define the A1 domain in terms of the sequence that delimits it and the sequences that define the N- and C- linkers and provide justification for the selection.

The justification is summarized in Results in lines 107-110: "The N and C-terminal boundaries of the constructs were based on allowing comparability to previous studies; references to the previous literature and a structural rationale for defining the minimal size of A1 as residues 1270-1463 and the maximal size of A1 plus linkers as residues 1240 to 1494 are described in Methods.” A new section in Methods, "Boundaries of A1 constructs" appears on lines 341-357.

2) Show the standard deviation of the fit instead of showing the average of two experiments in figure 2 and provide additional description for the methods used to purify protein constructs in the gel images as the bands appear differentially pure, in different amounts in figure 1.

Global non-linear least square fitting the traces at all analyte concentrations, we obtained standard fitting errors ranging from 1% to 5% of fit values. These fitting errors were smaller than the difference from the mean of 2 independent experiments shown in Figure 2; therefore, we have reported these results in the Figure 2 legend.

The proteins may have varied from 90 to 99% pure but that would not affect the measured quantities. The proteins were at different concentrations and were not adjusted to identical amounts for SDS-PAGE, which is explained in the legend. More details are given in Methods on purification and a summary has been added to the Results.

Reviewer #1 (Recommendations for the authors):

In this paper, the authors use several biophysical tools, including biolayer interferometry, isothermal calorimetry, protein denaturation, and hydrogen-deuterium exchange mass spectrometry, to systemically dissect the influence of the N- and C-terminal polypeptide linkers and O-glycan moieties of the VWF A1 domain on A1's activation state and binding affinity to platelet receptor GPIbα.

Through a series of elegant experiments, the authors concluded that the N-terminal linker together with its O-linked glycans stabilize the A1 in an inactive, native conformation. Removal of the N-linker or its O-linked glycans destabilizes the A1 structure and may activate the A1 into an intermediate, higher affinity state for GPIbα binding. The authors also concluded that the C-terminal linker and its O-linked glycans do not contribute to stabilizing the A1 structure, though they prevent GPIbα binding through steric hindrance and electrostatic repulsion. The data and the figure presented in the manuscript are overall clear, and the conclusions are consistently supported by the data from multiple experimental techniques. The reported results are highly interesting for the hematology and biophysics communities

I have only a few suggestions for revision and improvement, detailed below.

1) The length of the long C-linker is inconsistent among the constructs. While the constructs expressed from E. coli end at residue 1486, the mammalian constructs end at residue 1489. No further justification was given by the authors. Moreover, since the A2 domain starts at residue 1494, why didn't the authors include residues 1490-1493?

We have explained this above.

2) In regard to Figure 2, some of the global fits are not of high quality (e.g., A1+C with 150 mM NaCl). The authors should show the standard deviation of the fit instead of showing the average of two experiments.

Please see above.

3) It would be worth discussing the difference between the intermediate state identified by the thermal/chemical denaturation experiments and the more physiological mechanically tensioned structure. The intermediate state may or may not be the same high-affinity conformation induced by tensile pulling, since heat/urea destabilizes the entire A1 structure and tensile force/tension is not expected to spread into the structures protected by the long-ranged disulfide.

This was included in the original manusucript,: "In VWF concatemers, the ΔG required for activation by mechanical tension of the high affinity state was measured as 1.9 kcal/mol {Fu, 2017 #24942}, considerably lower than the difference in energy between the I and N states measured for A1+OGly+N+C here of 4.7 kcal/mol. These differences suggest that the directional nature of energy input by mechanical tension allows for more efficient activation of A1 than urea denaturation and may induce a distinct and less disordered intermediate state."

There are many interactions between residues internal and external to the disulfide bond, including hydrogen bonds, van der Waals contacts, burial from water. Therefore, we expect that force will have a significant indirect effect on the structure and stability of regions internal to the disulfide. This is mentioned in the Discussion.

4. Lines 223-224: It would be nice to specify the exact constructs and conditions that account for the 80-fold and 10-fold differences.

Done.

Reviewer #2 (Recommendations for the authors):

[…]

Some perceived weaknesses are as follows:

1. The authors appropriately report a gel image of their purified protein constructs. The bands on the gel appear differentially pure, in different amounts, and there is no description in the manuscript about the relative 'diffuse' nature of the bands for glycosylated protein constructs.

This information has been added to Results.

2. The study describes the rigorous treatment of their thermal and chemical denaturation curves, and they briefly address the difference between measuring protein unfolding via intrinsic tryptophan fluorescence vs. circular dichroism. None of the curves presented illustrate error bars, but they do provide a measure of error in the associated table. Their measurements in this experiment are exclusively dependent on the environment of a single tryptophan residue in the A1 domain structure; their interpretations of the unfolding curves are warranted, but not without the caveat that it is only measuring the local environment of the tryptophan and not the entire protein structure.

We have revised " Trp-1313 is in the α1-β2 loop that has HDX dynamics that are sensitive to the presence of the N-linker, locates near the GPIbα β-finger, is mutated in VWD type 2B, and appears to be a more discrete reporter of A1 conformational state than circular dichroism."

to read

"Trp-1313 is in the α1-β2 loop that has HDX dynamics that are sensitive to the presence of the N-linker, locates near the GPIbα β-finger, is mutated in VWD type 2B, and although it may report denaturation of only a portion of A1, Trp-1313 discriminates between the native, intermediate, and denatured states of A1 similarly to circular dichroism."

My recommendations for the authors concern the gel image in figure 1 and the Materials and methods section. For figure 1, the protein bands are underwhelming. The reader would have more confidence in the study if this gel demonstrated a more uniform concentration of the protein constructs used in the study. Some of the bands are quite large and dominate the figure while some are comparatively weaker (+N+C in particular).

This has been addressed as discussed above.

For the Materials and methods section, I question the accuracy of a couple statements: it is stated that you have NiCl in the purification buffers. Ni is divalent, so this should be more accurately reported as NiCl2.

Thanks. Done.

In the section on non-glycosylated proteins, the authors state that SHuffle cells are grown at 37C, but these cells are typically grown at 30C due to their thermolability (if this is correct, that is fine, but it diverges from the manufacturer's protocol).

"SHuffle, a novel Escherichia coli protein expression strain capable of correctly folding disulfide bonded proteins in its cytoplasm Julie Lobstein, Charlie A Emrich, Chris Jeans, Melinda Faulkner, Paul Riggs and Mehmet Berkmen" describes multiple growth temperatures. We used 30°C after induction.

Lastly, please define for Equation 1 what b and s are.

Thanks, they were defined many lines above the equation and are now defined again immediately before Equation 1.

Reviewer #3 (Recommendations for the authors):

[…]

It would be useful for the authors to define the A1 domain in terms of the sequence that delimits it and the sequences that define the N- and C- linkers. Is the A1 delimited by residues 1272 and 1458 or by residues 1264 – 1469 (Figure 1C)?

Sections on this have been added as described above.

The organization at the end of the Introduction perhaps could be improved. The last paragraph begins with "Here, to understand … we have made fundamental measurements …" preparing the reader for what the authors did. However, the authors proceed in the next few sentences to continue to cite observations from the literature.

Good point. We have broken this up into two paragraphs.

https://doi.org/10.7554/eLife.75760.sa2

Article and author information

Author details

  1. Klaus Bonazza

    Program in Cellular and Molecular Medicine, Boston Children’s Hospital and Department of Biological Chemistry and Molecular Pharmacology, Harvard Medical School, Boston, United States
    Contribution
    Formal analysis, Funding acquisition, Investigation, Writing – original draft, Writing – review and editing
    Competing interests
    No competing interests declared
  2. Roxana E Iacob

    Department of Chemistry and Chemical Biology, Northeastern University, Boston, United States
    Contribution
    Formal analysis, Investigation
    Competing interests
    No competing interests declared
  3. Nathan E Hudson

    Department of Physics, East Carolina University, Greenville, United States
    Contribution
    Funding acquisition, Investigation
    Competing interests
    No competing interests declared
  4. Jing Li

    Program in Cellular and Molecular Medicine, Boston Children’s Hospital and Department of Biological Chemistry and Molecular Pharmacology, Harvard Medical School, Boston, United States
    Contribution
    Formal analysis, Writing – review and editing
    Competing interests
    No competing interests declared
  5. Chafen Lu

    Program in Cellular and Molecular Medicine, Boston Children’s Hospital and Department of Biological Chemistry and Molecular Pharmacology, Harvard Medical School, Boston, United States
    Contribution
    Methodology, Writing – review and editing
    Competing interests
    No competing interests declared
  6. John R Engen

    Department of Chemistry and Chemical Biology, Northeastern University, Boston, United States
    Contribution
    Formal analysis, Funding acquisition, Methodology, Writing – review and editing
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-6918-9476
  7. Timothy A Springer

    Program in Cellular and Molecular Medicine, Boston Children’s Hospital and Department of Biological Chemistry and Molecular Pharmacology, Harvard Medical School, Boston, United States
    Contribution
    Conceptualization, Funding acquisition, Methodology, Supervision, Writing – original draft, Writing – review and editing
    For correspondence
    springer@crystal.harvard.edu
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-6627-2904

Funding

National Hemophilia Foundation (Postdoctoral fellowship grant)

  • Klaus Bonazza

Austrian Science Fund (Erwin Schrödinger Fellowship J4081-B21)

  • Klaus Bonazza

National Institutes of Health (F32HL126386)

  • Nathan E Hudson

National Institutes of Health (R01-HL148755)

  • Timothy A Springer

The funders had no role in study design, data collection, and interpretation, or the decision to submit the work for publication.

Acknowledgements

We would like to thank Prof. Thomas Wales for helpful discussions and a research collaboration with the Waters Corporation (JRE). We thank Sha Wang and Amy Xu for early work. We thank Yang Su for writing and using code to display HDX MS data for figure preparation. We acknowledge support from the National Hemophilia Foundation and from the Austrian Science Fund (FWF), Erwin Schrödinger Fellowship number J4081-B21, NIH grant F32HL126386, and NIH grant R01-HL148755.

Senior Editor

  1. Mone Zaidi, Icahn School of Medicine at Mount Sinai, United States

Reviewing Editor

  1. Yelena Ginzburg, Icahn School of Medicine at Mount Sinai, United States

Reviewer

  1. P Clint Spiegel, Western Washington University, United States

Publication history

  1. Received: November 22, 2021
  2. Preprint posted: January 13, 2022 (view preprint)
  3. Accepted: April 6, 2022
  4. Version of Record published: May 9, 2022 (version 1)

Copyright

© 2022, Bonazza 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

  • 476
    Page views
  • 73
    Downloads
  • 2
    Citations

Article citation count generated by polling the highest count across the following sources: Crossref, PubMed Central, Scopus.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

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)

  1. Klaus Bonazza
  2. Roxana E Iacob
  3. Nathan E Hudson
  4. Jing Li
  5. Chafen Lu
  6. John R Engen
  7. Timothy A Springer
(2022)
Von Willebrand factor A1 domain stability and affinity for GPIbα are differentially regulated by its O-glycosylated N- and C-linker
eLife 11:e75760.
https://doi.org/10.7554/eLife.75760

Further reading

    1. Biochemistry and Chemical Biology
    2. Structural Biology and Molecular Biophysics
    Ziva Vuckovic, Jinan Wang ... David M Thal
    Research Article

    Allosteric modulation of G protein-coupled receptors (GPCRs) is a major paradigm in drug discovery. Despite decades of research, a molecular-level understanding of the general principles that govern the myriad pharmacological effects exerted by GPCR allosteric modulators remains limited. The M4 muscarinic acetylcholine receptor (M4 mAChR) is a validated and clinically relevant allosteric drug target for several major psychiatric and cognitive disorders. In this study, we rigorously quantified the affinity, efficacy, and magnitude of modulation of two different positive allosteric modulators, LY2033298 (LY298) and VU0467154 (VU154), combined with the endogenous agonist acetylcholine (ACh) or the high-affinity agonist iperoxo (Ipx), at the human M4 mAChR. By determining the cryo-electron microscopy structures of the M4 mAChR, bound to a cognate Gi1 protein and in complex with ACh, Ipx, LY298-Ipx, and VU154-Ipx, and applying molecular dynamics simulations, we determine key molecular mechanisms underlying allosteric pharmacology. In addition to delineating the contribution of spatially distinct binding sites on observed pharmacology, our findings also revealed a vital role for orthosteric and allosteric ligand–receptor–transducer complex stability, mediated by conformational dynamics between these sites, in the ultimate determination of affinity, efficacy, cooperativity, probe dependence, and species variability. There results provide a holistic framework for further GPCR mechanistic studies and can aid in the discovery and design of future allosteric drugs.

    1. Structural Biology and Molecular Biophysics
    Seoyoon Kim, Daehyo Lee ... Duyoung Min
    Tools and Resources

    Single-molecule tweezers, such as magnetic tweezers, are powerful tools for probing nm-scale structural changes in single membrane proteins under force. However, the weak molecular tethers used for the membrane protein studies have limited the observation of long-time, repetitive molecular transitions due to force-induced bond breakage. The prolonged observation of numerous transitions is critical in reliable characterizations of structural states, kinetics, and energy barrier properties. Here, we present a robust single-molecule tweezer method that uses dibenzocyclooctyne (DBCO) cycloaddition and traptavidin binding, enabling the estimation of the folding 'speed limit' of helical membrane proteins. This method is >100 times more stable than a conventional linkage system regarding the lifetime, allowing for the survival for ~12 h at 50 pN and ~1000 pulling cycle experiments. By using this method, we were able to observe numerous structural transitions of a designer single-chained transmembrane (TM) homodimer for 9 h at 12 pN, and reveal its folding pathway including the hidden dynamics of helix-coil transitions. We characterized the energy barrier heights and folding times for the transitions using a model-independent deconvolution method and the hidden Markov modeling (HMM) analysis, respectively. The Kramers rate framework yields a considerably low speed limit of 21 ms for a helical hairpin formation in lipid bilayers, compared to μs scale for soluble protein folding. This large discrepancy is likely due to the highly viscous nature of lipid membranes, retarding the helix-helix interactions. Our results offer a more valid guideline for relating the kinetics and free energies of membrane protein folding.