Molecular Dynamics Simulated Unfolding of von Willebrand Factor A Domains by Force

Structural hierarchy of mechanical extensibility in human von Willebrand factor multimers

The von Willebrand factor (VWF) is a multimeric glycoprotein composed of 80- to 120-nm-long protomeric units and plays a fundamental role in mediating platelet function at high shear. The exact nature of the shear-induced structural transitions have remained elusive; uncovering them requires the high-resolution quantitative analysis of gradually extended VWF. Here, we stretched human blood-plasma-derived VWF with molecular combing and analyzed the axial structure of the elongated multimers with atomic force microscopy. Protomers extended through structural intermediates that could be grouped into seven distinct topographical classes. Protomer extension thus progresses through the uncoiling of the C_1-6 domain segment, rearrangements among the N-terminal VWF domains, and unfolding and elastic extension of the A₂ domain. The least and most extended protomer conformations were localized at the ends and the middle of the multimer, respectively, revealing an apparent necking phenomenon characteristic of plastic-material behavior. The structural hierarchy uncovered here is likely to provide a spatial control mechanism to the complex functions of VWF.

Exposure of Von Willebrand Factor Cleavage Site in A1A2A3-Fragment under Extreme Hydrodynamic Shear

Von Willebrand Factor (vWf) is a giant multimeric extracellular blood plasma involved in hemostasis. In this work we present multi-scale simulations of its three-domains fragment A1A2A3. These three domains are essential for the functional regulation of vWf. Namely the A2 domain hosts the site where the protease ADAMTS13 cleavages the multimeric vWf allowing for its length control that prevents thrombotic conditions. The exposure of the cleavage site follows the elongation/unfolding of the domain that is caused by an increased shear stress in blood. By deploying Lattice Boltzmann molecular dynamics simulations based on the OPEP coarse-grained model for proteins, we investigated at molecular level the unfolding of the A2 domain under the action of a perturbing shear flow. We described the structural steps of this unfolding that mainly concerns the β-strand structures of the domain, and we compared the process occurring under shear with that produced by the action of a directional pulling force, a typical condition of single molecule experiments. We observe, that under the action of shear flow, the competition among the elongational and rotational components of the fluid field leads to a complex behaviour of the domain, where elongated structures can be followed by partially collapsed melted globule structures with a very different degree of exposure of the cleavage site. Our simulations pose the base for the development of a multi-scale in-silico description of vWf dynamics and functionality in physiological conditions, including high resolution details for molecular relevant events, e.g., the binding to platelets and collagen during coagulation or thrombosis.

Intradimer forces and their implication for conformations of von Willebrand factor multimers

The largest blood glycoprotein von Willebrand factor (VWF) responds to hydrodynamic stresses in the bloodstream with abrupt conformation changes, thus increasing its adhesivity to platelets and collagen. Arterial and microvascular hemostasis relies on mechanical and physicochemical properties of this macromolecule. Recently, it was discovered that the mechanical properties of VWF are controlled by multiple pH-dependent interactions with opposite trends within dimeric subunits. In this work, computer simulations reveal the effect of these intradimer forces on the conformation of VWF multimers in various hydrodynamic conditions. A coarse-grained computer model of VWF has been proposed and parameterized to give a good agreement with experimental data. The simulations suggest that strong attraction between VWF D4 domains increases the resistance to elongation under shear stress, whereas even intermediate attraction between VWF C domains contributes to VWF compaction in nonsheared fluid. It is hypothesized that the detailed subdimer dynamics of VWF concatamers may be one of the biophysical regulators of initial hemostasis and arterial thrombosis.

A mechano-reactive coarse-grained model of the blood-clotting agent von Willebrand factor

The von Willebrand Factor (vWF) is a large blood glycoprotein that aids in hemostasis. Within each vWF monomer, the A2 domain hosts a cleavage site for enzyme ADAMTS13, which regulates the size of vWF multimers. This cleavage site can only be exposed when an A2 domain unfolds, and the unfolding reaction energy landscape is highly sensitive to the force conditions on the domain. Based on previous optical tweezer experimental results, we advance here a new activated A2 monomer model (AA2MM) for coarse-grained modeling of vWF that accurately represents the force-based probabilistic change between the unfolded/refolded states. A system of springs is employed to mimic the complex mechanical response of vWF monomers subject to pulling forces. AA2MM was validated by comparing monomer scale simulation results to data from prior pulling experiments on vWF monomer fragments. The model was further validated by comparing multimer scale Brownian dynamics simulation results to experiments using microfluidic chamber microscopy to visualize tethered vWF proteins subject to flow. The A2 domain unfolding reaction was studied in bulk flow simulations (pure shear and elongation flow), giving evidence that elongational flow drives the vWF size regulation process in blood. The mechanoreactive, coarse-grained AA2MM accurately describes the complex mechanical coupling between human blood flow conditions and vWF protein reactivity.

The combined effects of the microcirculatory status and cardiopulmonary bypass on platelet count and function during cardiac surgery

BACKGROUND

Cardiac surgery with cardiopulmonary bypass is associated with important changes in the microcirculation, usually attributed to endothelial dysfunction. Another common finding of cardiac surgery is postoperative thrombocytopenia and platelet loss of function.

OBJECTIVE

To investigate the association between microvascular flow pattern and postoperative changes in platelet count and function in cardiac surgery patients.

METHODS

Twelve adult cardiac surgery patients received microvascular circulation (sidestream darkfield sublingual mucosa analysis) and platelet count and function (multiple electrode aggregometry ADPtest and TRAPtest) assessment before and after cardiopulmonary bypass.

RESULTS

After cardiopulmonary bypass, sublingual microcirculation showed a significantly (P = 0.001) decreased microvascular flow index and increased heterogeneity index (P = 0.006). Platelet function significantly decrease after cardiopulmonary bypass both at ADPtest (P = 0.011) and TRAPtest (P = 0.002). Preoperative patterns of poor microvascular perfusion (low perfused vessels density and total vessels density) were significantly associated with lower values of post-cardiopulmonary bypass platelet function (ADPtest, P = 0.009, TRAPtest, P = 0.031) and count (P = 0.048).

CONCLUSIONS

A preoperative disturbance of the microcirculation is associated with a greater postoperative platelet dysfunction. Endothelial damage, chemical and mechanical stimuli are the possible link between the two patterns.

Long-ranged Protein-glycan Interactions Stabilize von Willebrand Factor A2 Domain from Mechanical Unfolding

von Willebrand Factor (vWF) is a large multimeric protein that binds to platelets and collagen in blood clotting. vWF A2 domain hosts a proteolytic site for ADAMTS13 (A Disintegrin and Metalloprotease with a ThromboSpondin type 1 motif, member 13) to regulate the size of vWF multimers. This regulation process is highly sensitive to force conditions and protein-glycan interactions as the process occurs in flowing blood. There are two sites on A2 domain (N1515 and N1574) bearing various N-linked glycan structures. In this study, we used molecular dynamics (MD) simulation to study the force-induced unfolding of A2 domain with and without a single N-linked glycan type on each site. The sequential pullout of β-strands was used to represent a characteristic unfolding sequence of A2. This unfolding sequence varied due to protein-glycan interactions. The force-extension and total energy-extension profiles also show differences in magnitude but similar characteristic shapes between the systems with and without glycans. Systems with N-linked glycans encountered higher energy barriers for full unfolding and even for unfolding up to the point of ADAMTS13 cleavage site exposure. Interestingly, there is not much difference observed for A2 domain structure itself with and without glycans from standard MD simulations, suggesting roles of N-glycans in A2 unfolding through long-ranged protein-glycan interactions.

A common mechanism by which type 2A von Willebrand disease mutations enhance ADAMTS13 proteolysis revealed with a von Willebrand factor A2 domain FRET construct

Rheological forces in the blood trigger the unfolding of von Willebrand factor (VWF) and its A2 domain, exposing the scissile bond for proteolysis by ADAMTS13. Under quiescent conditions, the scissile bond is hidden by the folded structure due to the stabilisation provided by the structural specialisations of the VWF A2 domain, a vicinal disulphide bond, a calcium binding site and a N1574-glycan.The reduced circulating high MW multimers of VWF in patients with type 2A von Willebrand disease (VWD) may be associated with mutations within the VWF A2 domain and this is attributed to enhanced ADAMTS13 proteolysis. We investigated 11 VWF A2 domain variants identified in patients with type 2A VWD. In recombinant full-length VWF, enhanced ADAMTS13 proteolysis was detected for all of the expressed variants in the presence of urea-induced denaturation. A subset of the FLVWF variants displayed enhanced proteolysis in the absence of urea. The mechanism of enhancement was investigated using a novel VWF A2 domain FRET construct. In the absence of induced unfolding, 7/8 of the expressed mutants exhibited a disrupted domain fold, causing spatial separation of the N- and C- termini. Three of the type 2A mutants were not secreted when studied within the VWF A2 domain FRET construct. Urea denaturation revealed for all 8 secreted mutants reduced unfolding cooperativity and stability of the VWF A2 domain. As folding stability was progressively disrupted, proteolysis by ADAMTS13 increased. Due to the range of folding stabilities and wide distribution of VWF A2 domain mutations studied, we conclude that these mutations disrupt regulated folding of the VWF A2 domain. They enhance unfolding by inducing separation of N- and C-termini, thereby promoting a more open conformation that reveals its binding sites for ADAMTS13 and the scissile bond.

Mutual A domain interactions in the force sensing protein von Willebrand factor

The von Willebrand factor (VWF) is a glycoprotein in the blood that plays a central role in hemostasis. Among other functions, VWF is responsible for platelet adhesion at sites of injury via its A1 domain. Its adjacent VWF domain A2 exposes a cleavage site under shear to degrade long VWF fibers in order to prevent thrombosis. Recently, it has been shown that VWF A1/A2 interactions inhibit the binding of platelets to VWF domain A1 in a force-dependent manner prior to A2 cleavage. However, whether and how this interaction also takes place in longer VWF fragments as well as the strength of this interaction in the light of typical elongation forces imposed by the shear flow of blood remained elusive. Here, we addressed these questions by using single molecule force spectroscopy (SMFS), Brownian dynamics (BD), and molecular dynamics (MD) simulations. Our SMFS measurements demonstrate that the A2 domain has the ability to bind not only to single A1 domains but also to VWF A1A2 fragments. SMFS experiments of a mutant [A2] domain, containing a disulfide bond which stabilizes the domain against unfolding, enhanced A1 binding. This observation suggests that the mutant adopts a more stable conformation for binding to A1. We found intermolecular A1/A2 interactions to be preferred over intramolecular A1/A2 interactions. Our data are also consistent with the existence of two cooperatively acting binding sites for A2 in the A1 domain. Our SMFS measurements revealed a slip-bond behavior for the A1/A2 interaction and their lifetimes were estimated for forces acting on VWF multimers at physiological shear rates using BD simulations. Complementary fitting of AFM rupture forces in the MD simulation range adequately reproduced the force response of the A1/A2 complex spanning a wide range of loading rates. In conclusion, we here characterized the auto-inhibitory mechanism of the intramolecular A1/A2 bond as a shear dependent safeguard of VWF, which prevents the interaction of VWF with platelets.

Misfolding of vWF to pathologically disordered conformations impacts the severity of von Willebrand disease

The primary hemostatic von Willebrand factor (vWF) functions to sequester platelets from rheological blood flow and mediates their adhesion to damaged subendothelium at sites of vascular injury. We have surveyed the effect of 16 disease-causing mutations identified in patients diagnosed with the bleeding diathesis disorder, von Willebrand disease (vWD), on the structure and rheology of vWF A1 domain adhesiveness to the platelet GPIbα receptor. These mutations have a dynamic phenotypical range of bleeding from lack of platelet adhesion to severe thrombocytopenia. Using new rheological tools in combination with classical thermodynamic, biophysical, and spectroscopic metrics, we establish a high propensity of the A1 domain to misfold to pathological molten globule conformations that differentially alter the strength of platelet adhesion under shear flow. Rheodynamic analysis establishes a quantitative rank order between shear-rate-dependent platelet-translocation pause times that linearly correlate with clinically reported measures of patient platelet counts and the severity of thrombocytopenia. These results suggest that specific secondary structure elements remaining in these pathological conformations of the A1 domain regulate GPIbα binding and the strength of vWF-platelet interactions, which affects the vWD functional phenotype and the severity of thrombocytopenia.

Force-sensitive autoinhibition of the von Willebrand factor is mediated by interdomain interactions

Von Willebrand factor (VWF) plays a central role in hemostasis. Triggered by shear-stress, it adheres to platelets at sites of vascular injury. Inactivation of VWF has been associated to the shielding of its adhesion sites and proteolytic cleavage. However, the molecular nature of this shielding and its coupling to cleavage under shear-forces in flowing blood remain unknown. In this study, we describe, to our knowledge, a new force-sensory mechanism for VWF-platelet binding, which addresses these questions, based on a combination of molecular dynamics (MD) simulations, atomic force microscopy (AFM), and microfluidic experiments. Our MD simulations demonstrate that the VWF A2 domain targets a specific region at the VWF A1 domain, corresponding to the binding site of the platelet glycoprotein Ibα (GPIbα) receptor, thereby causing its blockage. This implies autoinhibition of the VWF for the binding of platelets mediated by the A1-A2 protein-protein interaction. During force-probe MD simulations, a stretching force dissociated the A1A2 complex, thereby unblocking the GPIbα binding site. Dissociation was found to be coupled to the unfolding of the A2 domain, with dissociation predominantly occurring before exposure of the cleavage site in A2, an observation that is supported by our AFM experiments. This suggests that the A2 domain prevents platelet binding in a force-dependent manner, ensuring that VWF initiates hemostasis before inactivation by proteolytic cleavage. Microfluidic experiments with an A2-deletion VWF mutant resulted in increased platelet binding, corroborating the key autoinhibitory role of the A2 domain within VWF multimers. Overall, autoinhibition of VWF mediated by force-dependent interdomain interactions offers the molecular basis for the shear-sensitive growth of VWF-platelet aggregates, and might be similarly involved in shear-induced VWF self-aggregation and other force-sensing functions in hemostasis.

Protein unfolding under force: crack propagation in a network

The mechanical unfolding of a set of 12 proteins with diverse topologies is investigated using an all-atom constraint-based model. Proteins are represented as polypeptides cross-linked by hydrogen bonds, salt bridges, and hydrophobic contacts, each modeled as a harmonic inequality constraint capable of supporting a finite load before breaking. Stereochemically acceptable unfolding pathways are generated by minimally overloading the network in an iterative fashion, analogous to crack propagation in solids. By comparing the pathways to those from molecular dynamics simulations and intermediates identified from experiment, it is demonstrated that the dominant unfolding pathways for 9 of the 12 proteins studied are well described by crack propagation in a network.

Transient-state fluctuationlike relation for the driving force on a biomolecule

In experiments and simulations the force acting on a single biomolecular system has been observed as a fluctuating quantity if the system is driven under constant velocity. We ask the question that is analogous to transient state entropy production and work fluctuation relations whether the force fluctuations observed in the single biomolecular system satisfy a transient state fluctuationlike relation, and the answer is in the affirmative. Using a constant velocity pulling steered molecular dynamics simulation study for protein unfolding, we confirm that the force fluctuations of this single biomolecular system satisfy a transient-state fluctuationlike relation 1/γ(T,v) ln[P(v)(+f)/P(v)(-f)] = f. P(v)(±f) is the probability of positive and negative values of forces f = f · for a given unfolding velocity of magnitude v and the pulling direction n, nis the unit vector of n, and γ(T,v) is a factor that depends on initial equilibrium temperature T and the unfolding velocity. For different unfolding velocities we find that the system in the nonequilibrium pulling region displays substantial negative fluctuation in its unfolding force when velocity decreases. A negative value of force may indicate the emergence of refolding behavior during protein unfolding. We also find that γ(T,v) ~ T(-δ)v(α) and the system relaxation time τ(T,v) ~ T(δ)v(-(1+α), where α and δ are scaling exponents.

The importance of vicinal cysteines, C1669 and C1670, for von Willebrand factor A2 domain function

The von Willebrand factor (VWF) A2 crystal structure has revealed the presence of a rare vicinal disulfide bond between C1669 and C1670, predicted to influence domain unfolding required for proteolysis by ADAMTS13. We prepared VWF A2 domain fragments with (A2-VicCC, residues 1473-1670) and without the vicinal disulfide bond (A2-DeltaCC, residues 1473-1668). Compared with A2-DeltaCC, A2-VicCC exhibited impaired proteolysis and also reduced binding to ADAMTS13. Circular dichroism studies revealed that A2-VicCC was more resistant to thermal unfolding than A2-DeltaCC. Mutagenesis of C1669/C1670 in full-length VWF resulted in markedly increased susceptibility to cleavage by ADAMTS13, confirming the important role of the paired vicinal cysteines in VWF A2 domain stabilization.

Force-induced cleavage of single VWFA1A2A3 tridomains by ADAMTS-13

A disintegrin and metalloprotease with a thrombospondin type 1 motifs 13 (ADAMTS-13) regulates hemostasis by cleaving the folded A2 domain of von Willebrand factor (VWF). The cleavage is regulated by forces as it occurs in flowing blood. We tested the hypothesis that force-induced A2 domain unfolding facilitates cleavage using atomic force microscopy to pull single VWF A1A2A3 tridomain polypeptides by platelet glycoprotein Ibalpha or antibodies to measure time, distance, and force. Structural destabilization of A1A2A3 was induced by 5- to 80-pN forces, manifesting as an abrupt molecular length increase distributed around 20 and 50 nm, probably because of uncoupling A1A2A3 (or partially unfolding A2) and fully unfolding A2, respectively. Time required to destabilize A1A2A3 first increased (catch), reaching a maximum of 0.2 seconds at 20pN, then decreased (slip) with increasing force, independent of ADAMTS-13. The time required to rupture A1A2A3 exhibited a similar catch-slip behavior when pulled by glycoprotein Ibalpha but only slip behavior when pulled by antibody, which was progressively shortened by increasing concentration of ADAMTS-13 after (but not before) structural destabilization, indicating that cleavage of A2 requires the force-induced A2 unfolding. Analysis with a model for single-substrate trimolecular enzymatic kinetics estimated a cleavage rate k(cat) of 2.9 (+/- 59) seconds and a K(d) of 5.6 (+/- 3.4) nM for ADAMTS-13/A1A2A3 binding. These findings quantify the mechanical regulation of VWF cleavage by ADAMTS-13 at the level of single A1A2A3 tridomain.