OBSERVACIÓN DIRECTA POR CRIOTOMOGRAFÍA MICROELECTRÓNICA DE DEFORMACIÓN PLAQUETARIA INDUCIDA POR LA PROTEINA SPIKE (COVID-19)

 

• Christopher Cyrus Kuhn,
• Nirakar Basnet,
• Satish Bodakuntla,
• Pelayo Alvarez-Brecht,
• Scott Nichols,
• Antonio Martinez-Sanchez,
• Lorenzo Agostini,
• Young-Min Soh,
• Junichi Takagi,
• Christian Biertümpfel &
• Naoko Mizuno 

Nature Communications volume 14, Article number: 620 (2023) Cite this article

• 8627 Accesses

• 6 Citations

• 184 Altmetric

• Metrics details

Abstract

SARS-CoV-2 is a novel coronavirus responsible for the COVID-19 pandemic. Its high pathogenicity is due to SARS-CoV-2 spike protein (S protein) contacting host-cell receptors. A critical hallmark of COVID-19 is the occurrence of coagulopathies. Here, we report the direct observation of the interactions between S protein and platelets. Live imaging shows that the S protein triggers platelets to deform dynamically, in some cases, leading to their irreversible activation. Cellular cryo-electron tomography reveals dense decorations of S protein on the platelet surface, inducing filopodia formation. Hypothesizing that S protein binds to filopodia-inducing integrin receptors, we tested the binding to RGD motif-recognizing platelet integrins and find that S protein recognizes integrin α_vβ₃. Our results infer that the stochastic activation of platelets is due to weak interactions of S protein with integrin, which can attribute to the pathogenesis of COVID-19 and the occurrence of rare but severe coagulopathies.

Introduction

In 2019 a novel member of the Coronaviridae family was identified to cause a respiratory illness associated with an outbreak to a global extent^1,2. The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) likely emerged from a zoonotic transmission similar to previous epidemic pathogens SARS-CoV and MERS-CoV, and is the origin of the Corona virus disease 2019 (COVID-19) pandemic³. SARS-CoV-2 is an enveloped, positive-sense single-stranded RNA virus and belongs to the genus of betacoronaviruses. It is closely related to the bat coronavirus RatG and shows 79% homology to its predecessor SARS-CoV^4,5. Common symptoms are cough, pneumonia, and dyspnea, usually associated with a mild or moderate infection^6,7. A severe course of the disease can have dramatic outcomes like cardiovascular complications, respiratory failure, systemic shock, and multiple organ failure, leading to life-threatening conditions and potential death^8,9.

COVID-19 is associated with abnormalities in blood coagulation in severe cases. SARS-CoV-2 is detected in the blood samples collected from COVID-19 patients¹⁰. Although detected viral load is generally low, the amount of virus present in the plasma correlates with the severity of COVID-19^11,12. In a study observing the clinical aspects of COVID-19, 59.6% of the COVID-19 patients had viral loads in their blood. Particularly, in critical patients, a constant high amount of viral load (176 copies/ml) was observed, in contrast to the patients with mild cases (81.7 copies/ml)¹³. A low count of platelets (thrombocytopenia) together with the development of disseminated intravascular coagulation, myocardial infarction and non-vessel thrombotic complications are commonly observed in COVID-19 patients^8,14,15. Platelets isolated from COVID-19 patients have also shown abnormalities such as hyperactivity and an increase in their spreading behavior¹⁶. The causes of these abnormalities have been hypothesized as cytokines, antiphospholipid antibodies, interactions with other immune cells, and direct interaction between SARS-CoV-2 and platelets^{17,18,19,20,21,22}. Furthermore, isolated platelets from healthy donors mixed with SARS-CoV-2 or the SARS-CoV-2 spike (S) protein show a faster thrombin-dependent clot retraction and activate platelets independent of thrombin with upregulation of signaling factors²¹. A recent study further suggests the involvement of thrombin and tissue factor (TF) in the hyperactivation of platelets²³.

Owing to its relevance for the pathogenesis of SARS-CoV-2, S protein is central to the understanding of the molecular mechanisms of the SARS-CoV-2 infection. The petal shaped S protein forms a trimer that protrudes out of the viral membrane surface^24,25, and it is poised to engage with host cell receptors. SARS-CoV-2 and SARS-CoV S proteins show an overall similarity of 76%, although they have specific differences that impact their functions such as a furin cleavage site (PRRAR) unique to SARS-CoV-2, leading to its increased pathogenicity^26,27,28. Since the beginning of the pandemic, continuous mutations have been accumulated in S protein, making it challenging to battle infections and contain disease outbreaks.

The SARS-CoV-2 S protein consists of a S1 (residues 14–685) and a S2 (residues 686–1273) subunit, which are separated by host cell proteases²⁹. The receptor binding domain (RBD) in the S1 subunit is responsible for the attachment to host cells. Earlier structural studies by cryo-electron microscopy (cryo-EM) revealed open and closed states of S protein in the trimer^25,30. In the open state, in which RBDs are uplifted thereby revealing the receptor binding motif (RBM), S protein captures its major receptor angiotensin-converting enzyme 2 (ACE2)^25,30. In addition to ACE2, S protein is suggested to interact with several other host receptors including neuropilin-1 (NRP1)³¹ and CD147^32,33. Interestingly, SARS-CoV-2 is the only betacoronavirus containing an RGD (Arg-Gly-Asp) tripeptide motif in the RBD, which is typically recognized by several members of the integrin membrane receptor family^34,35. Bat-SL-CoVZC45 contains a RGD motif within the S protein but not within the RBD³⁶. Initial studies observed the involvement of integrin in the SARS-CoV-2 entry^37,38 and the binding of integrins to S protein^39,40. However, these findings await further validation and in-depth analysis. While there is fragmented information from ultrastructural, pathological, and patient studies connecting the SARS-CoV-2 severity to the impact of its spike protein on platelets, very little is known what causes the coagulation of platelets in the presence of SARS-CoV-2. Moreover, it is still under debate if ACE2 is present on the platelet surface, and therefore, it is not clear what the direct effect of SARS-CoV-2 on platelets is^16,19,21.

In this study, we probe the direct interaction of the SARS-CoV-2 S protein with platelets and visualize its effect on platelet morphology using live imaging and cryo-electron tomography (cryo-ET). In the presence of S protein, extensive elongation and increased spreading of platelets is observed. Notably, a population of abnormally shaped platelets resembling a proplatelet-like appearance is found, indicating an impact on the cytoskeleton-dependent platelet maturation process. Cryo-ET observations reveal actin-rich filopodia formation at the end of the elongated platelet and unexpectedly, there is a dense decoration of S protein on the filopodia surface. An orientation analysis reveals that S protein binds to the membrane surface with various angular distributions. Furthermore, based on the correlation of S protein binding and the filopodia formation, the interaction of platelet surface receptors and S protein are assessed in vitro. We find a weak but direct interaction of platelet residing RGD ligand integrin receptors α_vβ₃ and α₅β₁ with S protein but not with integrin α_IIbβ₃. Our results shed light on the abnormal behavior of platelets leading to coagulopathic events and micro-thrombosis caused by SARS-CoV-2 infection.

Results

Platelet deformation in the presence of SARS-CoV-2 S protein

One of the major pathological symptoms in COVID-19 patients is an abnormal platelet behavior. COVID-19 patients exhibit conditions such as thrombocytopenia, microvascular thrombosis, and coagulation, leading to the hypothesis that SARS-CoV-2 may directly cause platelet malfunctions^8,14,15. To assess the direct effect of SARS-CoV-2 S protein, we isolated platelets from healthy de-identified human blood donors and tested their morphological changes in the presence of S protein. Under several extracellular matrix proteins (ECMs; i.e., fibronectin and collagen I) or poly-L-lysine, platelets occasionally adhered to ECMs without S protein (Fig. 1, Supplementary Movie 1). Discoid-shaped platelets were predominant in these control conditions without S protein (Fig. 1A–C; “control”). However, in the presence of S protein, we observed the deformation of platelets to elongated morphologies (Fig. 1A–C, panel “Spike”). The deformation of platelets was quantified by the aspect ratio (ratio of the lengths of the longer to shorter axis, Fig. 1D), showing a median value of 1.960, 1.565, and 1.786 with collagen I, poly-L-lysine, and fibronectin respectively, while 1.687, 1.491 and 1.756 under control conditions without the addition of S protein. The corresponding circularities in the presence of S protein decreased to 0.559, 0.649, and 0.613 with collagen I, poly-L-lysine and fibronectin respectively, compared to 0.686, 0.668, and 0.650 under control conditions (Fig. 1E). S protein-induced elongation of platelets was observed under all coating backgrounds, indicating that the morphology change of platelets is a direct effect of S protein. The deformation of platelets reached to an extent of 32 µm in the long axis and 1.4 µm in the short axis in rare cases (Supplementary Fig. 1), though precise boundaries of platelets cannot be measured due to the limitations of the light microscopic resolution. The elongation of platelets (Supplementary Fig. 1) indicates the presence of proplatelets, a platelet precursor present in the circulatory system^41,42. Typically, proplatelets are processed by cytoskeletal remodeling and abscise into smaller platelets within the circulatory system^41,42,43. Our observation suggests that S protein may influence the maturation of platelets by acting on the cytoskeleton-based process necessary for proplatelet division. Occasionally, we observed a proplatelet shape, (Supplementary Fig. 1, white arrowhead), resulting in the formation of a wide hollow, ring with one or more bulges like gemstones on a ring. Interestingly, the effect of platelet deformation is less prominent in the poly-L-lysine background (Fig. 1B). We also observed less effect in the activation of platelets with S protein (Figs. 1F, 2D) with poly-L-lysine. As poly-L-lysine is not a platelet-specific adherence factor, it suggests that the deformation of platelets may be further aided by the binding to extracellular agents. To test the effects of S protein in the activation of platelet, a solid-phase sandwich ELISA assay and western blotting were performed. In the ELISA assay, platelet factor 4 (PF4) was measured to test the secretion of alpha-granules, a marker for the activation of platelets⁴⁴. The result showed an increase in secreted PF4 in platelet samples in the presence of S protein (Fig. 1G). Western blot analysis showed an increase in phosphorylated focal adhesion kinase (pFAK) as well, which responds to the signaling cascade of focal adhesion pathway⁴⁵. Interestingly, while pFAK was detected from platelets that adhered to the plate surface (1.45 ± 0.51 times higher than control without S protein), floating platelets did not show a detectable level of pFAK (Supplementary Fig. 1B–D). This indicates that the deformation of platelets is a reversible process that has little or only local influence.

Fig. 1: Comparison of platelet morphology with and without SARS-CoV-2 S protein.

A DIC images of platelets without (Control) and pre-incubated with 20 µg/ml S protein (Spike) on a collagen I support. * indicates collagen I fibers. B DIC images of platelets without and pre-incubated with S on a poly-L-lysine support. C DIC images of platelets without and pre-incubated with S protein on a fibronectin support. Scale bar: 5 µm (A–C). D Quantification of the aspect ratio of platelets (major axis/minor axis) on different coated surfaces, without and in the presence of S protein. The median aspect ratio is shown below the corresponding violin plot. The significance was determined by two-tailed Mann–Whitney U test. Number of platelets (collagen control: 442, collagen spike: 391, Poly-L-lysine control: 343, Poly-L-lysine spike: 430, Fibronectin control: 505, Fibronectin spike: 446) examined over at least three independent experiments from different donors. Each dot represents the value from a single platelet. E Quantification of the circularity of platelets on different coated surfaces, without and in the presence of S protein. The median circularity is shown below the corresponding violin plot. The significance was determined by two-tailed Mann–Whitney U test. Number of platelets (collagen control: 442, collagen spike: 391, Poly-L-lysine control: 343, Poly-L-lysine spike: 430, Fibronectin control: 505, Fibronectin spike: 446) examined over at least 3 independent experiments from different donors. Each dot represents the value from a single platelet. F Comparison of platelet activation, incubated with and without S protein, on Poly-L-Lysine (left), and on Fibronectin (right) Platelets with amoeba-like morphologies were defined as activated platelets. The median percentage of the activation is shown below the corresponding dot plot. The significance was determined by two-tailed Mann–Whitney U test. Number of platelets (collagen control: 442, collagen spike: 391, Poly-L-lysine control: 343, Poly-L-lysine spike: 430, Fibronectin control: 505, Fibronectin spike: 446) examined over at least three independent experiments from different donors. Each dot represents the value from a series of experiment. G Sandwich ELISA assay detecting PF4 release in the absence and presence of S protein. The mean ± SD is shown below the corresponding dot plot. The significance was determined by two-tailed Mann–Whitney U test. n = 3 independent experiments from n = 3 different donors. Source data are available as a Source Data file.

Full size image

Fig. 2: Platelet morphology depending on SARS-CoV-2 S protein concentration.

A DIC images of platelets in the presence of different amounts of S protein plated onto collagen I-coated surfaces. Scale bar: 5 µm. B Quantification of the aspect ratio of platelets on collagen I-coated surfaces, without and in the presence of different S protein concentrations. The median aspect ratio is shown below the corresponding violin plot. The significance was determined by two-tailed Mann–Whitney U test. n = 442 (control), n = 359 (0.2 µg/ml S protein), n = 443 (2 µg/ml S protein), n = 391 (20 µg/ml S protein) platelets examined over at least three independent experiments from different donors. Each dot represents the value from a single platelet. C Quantification of the circularity of platelets on collagen I-coated surfaces, without and in the presence of different S protein concentrations. The median circularity is shown below the corresponding violin plot. The significance was determined by two-tailed Mann–Whitney U test. n = 442 (control), n = 359 (0.2 µg/ml S protein), n = 443 (2 µg/ml S protein), n = 391 (20 µg/ml S protein) platelets examined over at least three independent experiments from different donors. Each dot represents the value from a single platelet. The plots for control and in the presence of 20 µg/ml spike protein in B and C are same as those in Fig. 1D and E, respectively. D Quantification of platelet activation on Collagen I depending on S protein concentration. The median percentage of the activation is shown below the corresponding dot plot. The significance was determined by two-tailed Mann–Whitney U test. n = 442 (control), n = 359 (0.2 µg/ml S protein), n = 443 (2 µg/ml S protein), n = 391 (20 µg/ml S protein) platelets examined over at least three independent experiments from different donors. Each dot represents the value from a single platelet. Each dot represents the value from a series of experiment. Source data are available as a Source Data file.

Full size image

On the surface of SARS-CoV-2, S protein is arranged with an average distance of 35 nm to its neighboring S protein with an estimated number of 24 S proteins per virus particle, though the reported distribution has no apparent geometrical order⁴⁶. In this way, S proteins are effectively locally concentrated on the viral surface. To reflect the local concentration effects and test the influence of the concentration of S protein on the deformation of platelets, various concentrations of S protein were added to platelets in tenfold steps from 0.2 µg/ml to 20 µg/ml (0.47 nM to 47 nM) and the morphology changes of platelets were assessed. We observed an increase in the elongation as well as activation of platelets in accordance with the concentration of S protein (Fig. 2). While S protein was able to cause the deformation of the platelets even at the lowest concentration (0.2 µg/ml, or 0.47 nM), more dominant effects were observed under high concentrations. Our observation suggests that higher concentrations of S protein have a more pronounced impact on platelets, which could mimic locally concentrated S protein on the viral surface.

Cryo-ET analysis shows S protein densely decorated on the platelet surface

To gain molecular insights into the morphological changes of platelets in the presence of S protein, we performed a cryo-ET analysis of platelets under a collagen type I background, both in the presence and absence of S protein (Supplementary Movie 2). Platelets incubated with S protein revealed extensive deformations (Fig. 3A–C, highlighted in purple, Supplementary Fig. 2) consistent with the light microscopic observations, while intact platelets exhibit their typical disc-like morphology (Fig. 3D–F, highlighted in yellow). High magnification observation (33,000×) revealed that the elongated morphology of platelets was facilitated by the remodeling of actin, forming filopodia-like architectures as narrow as 43 nm (Supplementary Fig. 2). Within the filopodia that we analyzed (Fig. 3G, I), we found actin filaments assembling into tightly packed bundles, comparable to those seen in pseudopodia of untreated platelets⁴⁷. Accompanied by longer actin filaments running parallel to the axis of the filopodia, shorter ones were bridging between the longer filaments and the membrane surface (Fig. 3H–J). The angular distribution of actin segments shows two peaks, namely at 10 and 80 degrees (Fig. 3H, arrows), those of 10° corresponding to the actin filaments that are running along the filopodia axis, while those of 80° connecting between the actin bundles and membranes. The difference of 70° between peaks reinforces the notion that the connections are made by the Arp2/3 complex, which is known to mediate actin branching at 70°⁴⁸. Control platelets showed filopodia-formation exclusively when contacting collagen I fibrils (Fig. 3D–F). Interestingly, in the presence of S protein, the morphological changes of platelets and their filopodia-like formation did not require the attachment to collagen I (Fig. 3A–C). This indicates that S protein itself has a contribution to the surface activation of platelets. Unexpectedly, we identified dense S protein densities decorating the external membrane surface of the filopodia protrusions, indicative of S protein (Fig. 3K, “+”, Supplementary Movie 2). In contrast, we only observed shorter faint densities of endogenous membrane receptors on the surface of intact platelets (Fig. 3K, “-”).

Fig. 3: Cryo-electron tomograms of platelets alone and in the presence of SARS-CoV-2 S protein on a collagen I support.

A, B Low magnification views of a platelet in the presence of S protein on collagen I. The dashed box in B represents the area of tomographic data collection in C. C Tomographic slice of the platelet protrusion in the presence of S protein. D, E Low magnification views of a platelet on collagen I. The dashed box in E represents the area of tomographic data collection in F. F Representative slice of the reconstructed tomogram of platelet plasma membrane. G Magnification of the filopodial structure from C with actin filaments running along the protrusion. H Angular arrangement of actin filaments along the platelet protrusion. I Length of the traced actin filaments of the tomographic reconstruction in G. J Length distribution of traced actin filaments depicted in I. Actin filaments ≥100 nm (12 in total) are not represented in the graph. K Magnified views on the platelet plasma membrane without and in the presence of S protein (E - extracellular, I –intracellular). Scale bars: A, D = 1 µm; B, E = 0.5 µm; C, F, G = 200 nm; K = 20 nm. Source data are available as a Source Data file.

Full size image

SARS-CoV-2 S protein binds to the platelet membrane surface flexibly

To further analyze the S protein densities on the platelet membrane surface, we manually selected and extracted the densities from 8 tomograms and analyzed them using subtomogram averaging approaches (Fig. 4A). To facilitate a focused alignment of the decorating protein without the influence of the membrane density, the membrane signal was subtracted using PySeg⁴⁹. We obtained a 3D-averaged density at a resolution of 13.8 Å (Fig. 4A, Supplementary Fig. 3) showing a characteristic trimeric shape of 15 nm in size. The obtained structure agreed well with our near-atomic resolution structure of S protein using the same protein batch for single-particle analysis (Fig. 4A, lower panel and Supplementary Fig. 3) as well as previously published structures^25,30,46,50, (Fig. 4A, fitted PDB 6vxx), validating the identity of S protein decorations. Furthermore, the analysis of the particles without application of C3 symmetry showed an asymmetric uplift of the tip of the S1 surface that connects to the extra density (Fig. 4A, right, circled in magenta). This shows that one of the three RBD domains of S protein is lifted up upon its binding to the host cell receptors. The extra density connected to S protein represents the density from the host cell, likely the platelet surface receptor recognized by S protein (Fig. 4I).

Fig. 4: SARS-CoV-2S protein reconstruction and membrane decoration analysis.

A Top-Left: Structure of S protein with closed conformation fitted in the subtomogram reconstruction. Top-Right: Structure of S protein calculated without C3 symmetry, revealing the uplifted RBD domain connected to additional densities from the host platelets. The additional densities connected to the open RBD domain is circled in magenta. Bottom-left: SPA-based structure of S protein in the closed conformation. Bottom-right: SPA-based structure of S protein in the open conformation. B Nearest neighbor distance distribution of S protein densities on the platelet surface membrane. The distances are calculated using the originally manually picked coordinates. The median distance between two S protein is 27.3 nm C Densities of the reconstructed S protein back-projected to the segmented platelet plasma membrane. The tomogram lacks top and bottom due to the missing wedge effect of tomographic data collection. D Orientation of S protein on the membrane surface. The scheme depicts the angle determination of S protein C3 axis and the normal of the platelet plasma membrane. The range from 45 to 135° was observed to be favorable for S protein interaction with the platelet surface. E Schematic depiction of S protein orientation in different angles towards platelet plasma membrane. F Distance of S protein from the membrane. The median distance from the center of S protein to the membrane is 16 nm. The box plot represents 25 and 75 percentiles (8.6 and 27 nm, whiskers). n = 539 S protein was assessed. G Visualization of the platelet plasma membrane curvedness and the position of S protein on the platelet surface. The top and bottom edge of the segmented membrane was excluded from the estimation. H Membrane curvature comparison of S protein bound and surface protein free areas on the platelet plasma membrane. The box plots represent 25, median and 75 percentiles. Control: 25% 0.043, median 0.091 and 75% 0.15. + S protein: 25% 0.11, median 0.14 and 75% 0.17. n = 2237 (control) and n = 539 (+Spike) points were assessed. I Additional densities (red allow-heads) between picked S protein and platelet plasma membrane. Scale bar: I: 20 nm. Source Data are available as a Source Data file.

Full size image

To assess how S protein recognizes the platelet surface, the alignment parameters of the individually analyzed S protein densities were applied to the 3D average and plotted back to the original tomograms (Fig. 4B–H). The distribution analysis of neighboring S protein showed that the peak population had a distance of 27.3 nm (median) apart, but some molecules were also more sparsely distributed (Fig. 4B). This measurement corresponds to a density of one S protein on a surface area of up to 585 nm² (a radial surface of 585 nm² is covered by one S protein), although no apparent periodical distribution was detected. Judging from the diameter of S protein (~17 nm), neighboring S protein closely located next to each other. S protein bound to the platelet membrane at a distance of 16 nm between the center of S protein to the membrane surface (Fig. 4F) and interestingly, with a wide range of angular distribution (Fig. 4D–E) with respect to the membrane surface, indicating its flexible attachment to the platelet surface. In addition, S protein binds to a slightly more curved membrane surface (Fig. 4G, H). This may be reflected by the fact that the binding of S protein induces filopodia formation with a membrane protrusion. Taken together, these results indicate that S protein approaches the platelet surface from various geometrical orientations to accommodate and enhance the docking to the membrane surface. Similarly, a broad angular distribution of S protein has been observed from the viral surface, due to several kinked points in the stalk region^46,50. Together with our observation, it suggests the orientational adjustments from both sides of S protein, namely the receptor binding S1 subunit and the stem side at the root on the virus, maximize the efficiency of the attachment of S protein to the host cell receptors. Consistent with the observed additional density on the lifted RBD domain (Fig. 4A, right), some of the tomograms showed extra densities bridging between plasma membrane and S protein (Fig. 4I, red arrows), presumably those of platelet receptors recognizing S protein. However, subtomogram analysis only yielded a faint density (Fig. 4A, right) without features, also suggesting a flexible attachment of S protein to its receptor on the membrane surface.

Platelet deformation in the presence of pseudotyped viral particles

After characterizing the effects of S protein on platelets, we hypothesized that locally concentrated S protein on a globular viral surface would be advantageous for increasing the local concentration and evaluated the influence of SARS-CoV-2 pseudo virus-like-particles (VLPs) on platelet deformation. We either generated or obtained SARS-CoV-2 pseudotyped VLPs that are fully intact vesicle-like entities as validated by negative staining EM (Supplementary Fig. 4A) and by their ability to infect HEK-293T-hACE2 cells (HEK-293T cells constitutively expressing ACE2 receptor, Supplementary Fig. 4B and C). The viral titer determined by flow cytometry was ~10⁴–10⁶ particles/ml, comparable to the reported preparation of SARS-CoV-2 pseudo VLPs⁵¹, however, it was low compared to VSV-G based lentiviruses. This low titer did not allow us to readily detect changes in platelet morphologies by live platelet imaging. However, we were able to find an example of a particle located in close proximity to a platelet filopodium (Supplementary Fig. 5B–F). Cryo-electron tomography revealed that the closest distance between this particle and the membrane surface of the filopodium was 20 nm (Supplementary Fig. 5G, H), similar to that measured for S protein alone (Fig. 4F 16 nm, from the center of S protein to membrane). The cross-section views of the particle showed decorations of proteins on the membrane surface (Supplementary Fig. 5I, red arrowheads), altogether suggesting that this vesicle may be a bound VLP. In comparison, we also found examples of extracellular vesicles with a similar shape and size (Supplementary Fig. 5A) that appeared to originate from intracellular vesicles, i.e. exosomes, containing alpha granules (Supplementary Fig. 5A, left) or budding out from the concave surface of plasma membrane, instead of filopodia and indicative of vesicle release through fusion of lysosome and plasma membrane. These results corroborate our in vitro data of purified S protein inducing morphological changes in platelets.

Integrin receptors recognize SARS-CoV-2 S protein

Several cell receptors were reported to recognize S protein. However, the presence of ACE2, the major S protein receptor, on the platelet surface is still inconclusive^16,19,21.Therefore, the relevance of the ACE2 receptor for platelet malfunction is still an open question. In contrast, integrin receptors are the major class of receptors expressed in platelets. Considering our structural analysis (Fig. 4) and the possibility that ACE2 is not abundantly expressed on platelets, we hypothesized that S protein may directly recognize integrin receptors. Interestingly, the RBD domain of S protein contains a stretch with an RGD-motif, which is a common motif among integrin ligands³⁵ and a direct interaction of tissue integrin α₅β₁ and SARS-CoV-2 S protein has been shown³⁹. We therefore tested the binding of S protein to known platelet integrin receptors α_IIbβ₃, α_vβ₃, and α₅β₁, enriched in the tissue but also expressed on platelets, all recognizing the RGD ligand motif. We used ELISA-like solid-phase equilibrium binding assays to detect the interaction of S protein with integrins (Fig. 5A). We detected the binding of integrin α₅β₁ and α_vβ₃ to S protein, while integrin α_IIbβ₃ does not have an apparent interaction with it (Fig. 5B). The extent of binding is most prominent with integrin α_vβ₃, while integrin α₅β₁ showed only a weak interaction. However, it should be noted that the observed binding of tested integrins was much weaker (<10-fold) compared to those for physiological integrin ligands: vitronectin for α_vβ₃, fibrinogen for α_IIbβ₃ and fibronectin for α₅β₁. Encouraged by our results, we tested the effect of platelet activation in the presence of cilengitide, a cyclic RGD pentapeptide⁵² that blocks the binding of integrin to RGD motif-containing extracellular ligands, and indeed the activation was reduced (Fig. 5C). These observations generally agree with a recent discussion of the relevance of integrin recognition by SARS-CoV-2 for vascular dysregulation³⁸.

Fig. 5: Interaction of integrin receptors with SARS-CoV-2 S protein and various ECM proteins.

A Scheme of the experimental setup. Immunoplates were coated with either S protein or ECM proteins. The ligands were incubated with various integrin-velcro constructs. Biotinylated anti-velcro polyclonal antibody, subsequently coupled to Streptavidin-HRP, was used to label ligand-bound integrins. Detection of the binding was measured at 405 nm, 10 min after addition of ABTS. The scheme was created with Biorender.com. B Binding of integrins α_vβ₃, α_IIbβ₃, and α₅β₁ to S protein and their physiological ECM ligands: α_vβ₃ - vitronectin, α_IIbβ₃ - fibrinogen, α₅β₁ - fibronectin. Data are from a representative experiment out of three independent ones, and shown as mean ± SD. The significance was determined by an two-tailed unpaired t test. C Quantification of platelet activation on Collagen I depending on S protein concentration and integrin inhibitor cilengitide. The significance was determined by a two-tailed Mann–Whitney U test. n = 363 (control), n = 278 (control + CTD), n = 339 (+Spike), n = 358 (+Spike, +CTD) platelets examined over at least three independent experiments from different donors. Source data are available as a Source Data file.

Full size image

Discussion

SARS-CoV-2 has shown unique pathological symptoms that can lead to a wide range of coagulopathic events in severe cases. In our study, we probed the direct effect of S protein to the change in morphology of platelets at a molecular level, and we directly visualized the binding of S protein to the platelet surface (summarized in Fig. 6). We hypothesized that the binding of the SARS-CoV-2 is mediated by integrin receptors based on the following reasons; (1) the activation of platelets is governed by filopodia formation, (2) filopodia formation is initiated by integrin receptors, (3) the major receptors on the platelets are integrin receptors and (4) SARS-CoV-2 S protein contains a RGD sequence in the RBD, which is recognized by a subtype of integrin, and therefore we tested the interaction of platelet-expressed integrins with S protein. Our integrin inhibition experiment using cilengitide and in vitro solid-phase binding assays support this hypothesis, particularly with the possibility that S protein recognizes integrin α_vβ₃. The binding of S protein to integrin was much lower compared to the interaction of integrins with their physiological ligands, and interestingly, we did not detect the binding to the major platelet integrin α_IIbβ₃. Previously, an increased binding of the activated integrin α_IIbβ₃ antibody PAC-1 to platelets was observed in the presence of S protein²¹. This may be due to an inside-out effect, in which the outside-in signaling is activated by the direct binding of S protein to integrin α_vβ₃ and in turn, α_IIbβ₃ would get activated through the intracellular signaling (inside-out). We surmise that the weak affinity of S protein to platelet integrin receptors and the reversible binding, may reflect the fact that blood clotting defects observed in patients are rare complications and occur in severe cases of COVID-19. However, here we should also note that there are other receptors on platelets that may also be accountable for the interaction with S protein^33,53 and combinatory effects of the binding of S protein to multiple receptors may also occur.

Fig. 6: Schematic representation of potential SARS-CoV-2 S platelet interaction.

First, S protein binds to receptors on the platelet surface, causing the deformation and priming the activation. Protrusions are forming as a consequence of actin remodeling. This leads to the activation of platelets by the formation of filopodia and the stabilization of the cytoskeleton network. The scheme was created with Biorender.com.

Full size image

SARS-CoV-2 is found in the blood stream of COVID-19 patients¹⁰, and an open question is how it can lead to rare but severe coagulation defects. We showed that the deformation of platelets itself does not always alter their intracellular signaling (Supplementary Fig. 1), or induces activation. It rather appears that platelets exposed to S protein are primed for the activation upon further stimuli, such as the attachment to an adhesion surface. Based on this observation, we speculate that the combination of the direct binding of S protein to platelets and other identified coagulation factors may induce a synergistic and irreversible activation of platelets, leading to coagulation. During SARS-CoV-2 infection, several other procoagulant players are active, for example the formation of neutrophil extracellular traps¹⁸, the release of TF²³, elevated fibrinogen levels⁵⁴ and dysregulated release of cytokines⁵⁵, creating a hypercoagulative environment in the context of COVID-19.

In our study, we visualized the adaptable attachment of S protein to the platelet plasma membrane with a high degree of flexibility for the engagement to continuously curved membrane surfaces (Fig. 4D, E). Similarly, it has been reported that the stalk domain of S protein proximal to the viral membrane surface contains three hinges, presumably allowing the flexible motion of individual S protein on the viral surface to adapt to curved host cell surfaces⁵⁰. This dual flexibility likely increases the probability for S protein to attach to a host cell receptor, thus, allowing an efficient action of S protein to the membrane surface.

Methods

Platelet Isolation

Human platelets were prepared from the blood of de-identified healthy donors according to the available protocol²¹, with minor modifications as described below. Immediately after blood was drawn from the donors, it was centrifuged for 20 min at 200 × g at RT. The top half of the platelet-rich plasma (PRP) was transferred to a fresh tube and gently mixed with an equal amount of HEP buffer (14 mM NaCl, 2.7 mM KCl, 3.8 mM HEPES, 5 mM EGTA, 1 μM Prostaglandin E1, pH 7.4.). To remove remaining cells other than platelets, the PRP solution was centrifuged for 20 min at 100 × g. Three fourth of the supernatant was carefully transferred to a fresh tube, and platelets were pelleted for 20 min at 800 × g. The supernatant was discarded, the pellet was then washed with a solution containing 10 mM sodium citrate, 150 mM NaCl, 1 mM EDTA, 1% (w/v) dextrose, pH 7.4 and resuspended in Tyrode’s buffer (134 mM NaCl, 12 mM NaHCO₃, 2.9 mM KCl, 0.34 mM Na₂HPO₄, 1 mM MgCl₂, 10 mM HEPES, 5 mM Glucose, 3 mg/ml BSA, pH 7.4). The isolated platelets were rested for 45–60 min. The platelet concentration was estimated using a hemacytometer.

Preparation of Coating

Poly-L-Lysine (Sigma #P2636) or ECM proteins fibronectin (R&D systems #3420-001-01) and collagen-I (Chrono-Log #385) were diluted to working concentration 0.01% for Poly-L-Lysine, 10 µg/ml for Fibronectin and 25 µg/ml for Collagen I respectively. Poly-L-Lysine and fibronectin were applied to imaging chambers (μ-Slide 8 well glass bottom, ibidi #80827-90) over night at 37 °C. Collagen I coating was applied for 30 min at 37 °C. Coated imaging chambers were washed with PBS and afterwards blocked with 1% BSA in PBS.

Platelet and S protein incubation

S protein (Cube Biotech #28703) and isolated platelets were mixed at final concentrations of 0.2 µg/ml, 2 µg/ml and 20 µg/ml of S protein and 1.5 × 10⁷ platelets/ml in Tyrode’s buffer. The samples were incubated for 4 h at 37 °C, and further diluted to 0.75 × 10⁷ platelets/ml. Controls were prepared by adding Cube Biotech’s S protein buffer (20 mM HEPES, 150 mM NaCl, 0.01% LMNG, pH 7.5) instead of S protein.

Differential inference contrast (DIC) microscopy and analysis

Prior to the imaging, the coated imaging chambers were fixed inside the live-cell imaging chamber and the chamber and water bath temperature were set to 37 °C. The platelet and S protein mixture was transferred to the coated imaging chambers using wide orifice tips. Platelets were settled for 15 min, and subsequently live imaging was performed for 60 min, with a frame time of 2 min, except one series of experiment recorded for 58 min. Acquired frames were analyzed using FIJI⁵⁶. The platelet shape was manually tracked and segmented every three frames using ‘freehand selection’ tool over the imaging period. The overall circularity and aspect ratio of the platelets were quantified using the FIJI measurement plugin with options of ‘circularity’ and ‘aspect ratio’. Circularity is calculated by 4π(area/perimeter²). A circularity value of 1.0 indicates a perfect circle and lower values indicate an elongated polygon. Aspect ratio is defined as the ratio of the major axis and minor axis of the fitted ellipse. The quantified data was analyzed using GraphPad PRISM. The measurements were done at least three times using platelets from at least 3 different donors.

Platelet activation in the presence of integrin inhibitor cilengitide

Isolated platelets were pre-incubated at 111 µg/ml cilengitide (Millipore-Sigma #ML1594) for 20 min at 37 °C. Subsequently, S-protein or vehicle control was added to the platelets so that the final concentration of cilengitide becomes 100 µg/ml. The platelet activation was assessed by live imaging in the same way as the DIC assay without cilengitide.

PF4 sandwich ELISA assay

Supernatant from platelets seeded on Collagen I either incubated with or without S-protein was collected. The PF4 concentration was determined using the Human PF4 CatchPoint® SimpleStep ELISA® Kit (Abcam #ab278096), according to the manufacturer’s protocol. Prepared samples, standards and antibody cocktail were added to appropriate wells and incubated for 1 h at RT while gently shaking. Subsequently, the solution was removed and wells were washed with 1× wash buffer PT. CatchPoint HRP Development Solution was added and incubated for 10 min in the dark while gently shaking. The fluorescence was measured using BioTek Synergy H1 plate reader at excitation at 530 nm and emission at 590 nm.

Immunoblot

Platelets were seeded on 8-well chamber cover glass slides coated with Collagen I. After exposing the platelets with 20 µg/ml spike protein or its vehicle control, the floating and adherent platelets were separately lysed in NP‑40 lysis buffer (50 mM Tris‑HCl pH 7.4, 150 mM NaCl, 1 mM NaF, 1 mM Na3VO4, 1% IGEPAL® CA‑630, 1 mM EDTA, protease inhibitors). The total protein concentration was determined using the Pierce™ Rapid Gold BCA Protein Assay Kit (Thermo Fisher #53226). 30 µg of total protein sample was mixed with Laemmli buffer (180 mM DTT (Sigma #D9779), 4% SDS (VWR #442444H), 160 mM Tris-HCl pH 6.8, 20% glycerol (VWR #24388.295), bromophenol blue). The samples were then heated at 95 °C for 5 min, spun down and were loaded on a BoltTM 4–12% Bis‑tris (Invitrogen #NW04120BOX), and SDS-PAGE was run at 180 V for 40 min in 1× MES buffer and subsequently transferred to a PVDF membrane using the Bio-Rad Trans-Blot® Turbo system according to the manufacturer’s instructions. The membrane was incubated with 4% ECL Blocking Agent (Fisher Scientific # 45001197) prepared in 1× TBST buffer (20 mM Tris‑HCl, 150 mM NaCl, 0.1% Tween® 20) at 4 °C for 3 h. When blotting against phospho-antibodies, the membrane was incubated with 4% BLOCK ACE (Bio-Rad # BUF029). The primary antibodies (Phospho-FAK (Tyr397) Polyclonal Antibody (Thermo Fisher, 44–624 G; 1:1000), Anti-GAPDH antibody Mouse monoclonal (Sigma #8795; 1:1000) for control were incubated at 4 °C overnight. After three washes with TBST, the membranes were incubated with secondary horseradish peroxidase-coupled goat anti–rabbit (Bethyl Laboratories, A120-201P) or anti-mouse (Bethyl Laboratories, A90-516P) antibody (1:10,000) for 60 min. The membranes were washed in TBST and the chemiluminescence signal was revealed by incubating with ECL substrate (Bio-Rad) for 2 min and imaged using Amersham™ Imager 600. The signals were quantified using FIJI and they were normalized using GAPDH signals. Anti-GAPDH antibody was used as a control as regulation of cytoskeleton components such as actin and tubulin may occur upon activation of platelet.

Cryo-electron tomography sample preparation and data acquisition

Quantifoil grids (MultiA Au200 & SiO₂ R1/4 Au200) were glow discharged using a Pelco easiGlow^TM at negative discharge, 15 mA plasma current and 0.38 mbar residual air pressure for 45 s. The glow discharged grids were placed in a cell culture dish (greiner #627170), covered in diluted coating solution, and incubated. After coating, the grids were washed in PBS and blocked with 1% BSA in PBS. After blocking, grids were covered in Tyrode’s buffer. Platelets, preincubated with 9.75 µg/ml S protein or control buffer (20 mM HEPES, 150 mM NaCl, 0.01% LMNG, pH 7.5) for 4 h, were added to cover the grids using a wide orifice and adhered for 1 h. 3 μl Tyrode’s buffer was added to grids with adhered platelets and vitrified in liquid ethane using a Thermo Fisher Scientific Vitrobot MarkVI, conditioned at 37 °C and 100% humidity.

The data was collected using a Titan Krios (Thermo Fischer Scientific), equipped with a Gatan Quantum 967 LS and K3 Summit direct detector at an acceleration voltage of 300 kV. Tilt-series were collected from −60° to 60° with 2° angular increment with a defocus range between −3 µm to −5  µm using a dose-symmetric acquisition scheme in SERIAL-EM software⁵⁷. At a nominal magnification of 33,000×, corresponding to a final pixel size of 2.76 Å. The total accumulated electron dose was 123 e⁻/Ų. Images were acquired as six-frame movies in super-resolution mode. A total of 8 tilt-series were assessed for this study.

Cryo-electron tomography reconstruction and segmentation

Images were motion-corrected and filtered according to their cumulative dose using the software MotionCor2⁵⁸. The tilt-series was aligned using the IMOD ETOMO package⁵⁹. Tomograms were reconstructed, unbinned and 4x binned, from aligned stacks as weighted back-projection in IMOD. The contrast of the tomograms was increased by applying Matlab based deconvolution filter⁶⁰. The CTF was estimated using GCTF⁶¹ and then the IMOD ctfphaseflip implementation was used for phase correction⁶².

Tomograms were manually segmented using AMIRA (Thermo Fisher Scientific) or 3DMOD (IMOD). The data was plotted using GraphPad PRISM.

Actin analysis

For the actin analysis, we analyzed 282 actin filaments. Actin was segmented manually using IMOD software⁵⁹. The IMOD model files were converted into coordinate files using IMOD (model2point). Each actin filament was segmented into 3 nm and 10 nm spaced segments for further analysis. For each actin filament the following two parameters were calculated: (a) angle of individual 10 nm actin segment against the longitudinal axis of platelets protrusion and (b) length of each actin filament. For determination of angle, a vector pointing towards the longitudinal axis of platelets is taken as reference vector and angle between this vector and each actin segment was calculated using Python3 NumPy library. Similarly, the length of each actin filament was calculated by adding the distance of all the segment present in each filament by using Python NumPy and SciPy libraries. The data was shown as histogram plots using GraphPad PRISM.

Subtomogram averaging of S protein and distance analysis

4167 S protein particles were manually picked from 4-times binned tomograms using IMOD softwared (3dmod). The coordinates of the picked particles were then transferred to RELION-3⁶³, and particles were extracted to a box with the size of 120 pixel from the original, unbinned electron tomograms. The extracted particles are located on the membrane surface and the membrane densities can interfere with the alignment process. To computationally supress the membrane density, Pyseg⁴⁹ scheme was used. Pyseg uses a basis of discrete Morse theory^64,65. Initially, all subvolumes where aligned with respect to the membrane using RELION3 as its density is stronger than the density of the protein. Afterwards, membrane densities for each subvolume are suppressed by assigning random background values to membrane voxels, membrane and background voxels are identified with an input mask. The initial template used for the alignment was our SPA 3D reconstruction of S protein but low pass filtered to 60 Å. At this resolution, only general shape and size of S protein was visible. The initial alignment was done using 3D auto-refine and 3D classification schemes from RELION-3. The particles were divided into 4 classes by 3D classification and the classes that showed the remaining membrane densities and noises were discarded for the further reconstruction. The class that showed the most features was selected and the final reconstruction was performed using 976 particles with C3 symmetry and with the mask that is created from the reconstruction of the previous run. The resolution was estimated by comparing the FSC of two separately computed averages from odd and even half-sets from the final refinement, the standard procedure available from RELION-3. The final resolution was estimated to be 13.8 Å with the FSC 0.143 criterion. The refinement without C3 symmetry was also performed, revealing the resolution of 20 Å.

Distances between S protein particles were calculated from the coordinates using Python3 (numpy and scipy libraries)⁶⁶. Positions of S protein were defined by manual picking. For each particle, the closest neighboring distance was plotted into the distance distribution histogram using GraphPad PRISM.

Membrane curvature analysis and determination of S protein orientation and distance from the membrane

The segmentation of the membrane was done manually using AMIRA (Thermo Fisher) software. The membrane curvature was then determined using python-based software PyCurv⁶⁷ using the standard workflow (https://github.com/kalemaria/pycurv). The segmented membrane from the binned tomogram was used as the input. The software then converts this segmentation, i.e., a set of voxels, into a surface, mesh of triangles. This surface of triangular mesh is converted into a surface graph, normal vectors and local curvature was then computed for every triangle center.

In order to measure the distance between the density corresponding to S protein and the membrane, the Euclidean distance between the refined coordinates of S protein particle to all the triangles on the membrane surface was calculated and the smallest distance was considered. Similarly, to determine the orientation of S protein with respect to the membrane, the angle between the vector pointing towards the longest axis of S protein and the normal vector of the closest triangle was calculated (Fig. 4D–E). Both distance and angular orientation was calculated using Python3 (numpy and scipy libraries). The position and orientation of S protein was visually assessed by the Place Object plug-in⁶⁸ in Chimera⁶⁹.

Pseudotyped SARS-CoV-2 S lentiviral particle production and transduction

Pseudotyped S lentivirus particles were either purchased (BPI #79981-1), obtained from BEI resources (BEI #NR-53818), NIAID, NIH or produced by adapting the published protocols^51,70. Necessary reagents for the production were obtained from BEI resources, NIAID, NIH. Briefly, X-Lenti 293T cells (Takara #632180) plated in six-well plates were co-transfected with 1 µg of lentiviral backbone ZsGreen (BEI #NR-52516); 0.22 µg each of the helper plasmids HDM-Hgpm2 (BEI #NR-52517), pRC-CMV-Rev1b (BEI #NR-52519), HDM-tat1b (BEI #NR-52518); and 0.34 µg of viral entry protein (SARS-CoV-2 Spike, BEI #NR-53742) or pCMV-VSVG (a gift from B. Weinberg, Addgene plasmid #8454) using 8 µl of TransIT-293 (Mirus Bio #MIR 2705) transfection reagent per well. On the next day, media was exchanged and collected after 60 h post transfection. The virus-containing supernatants were passed through a 0.45 µm filter and stored at −80 °C. To test the transducing ability of the viral particles, different volumes of the viral supernatant were added to HEK-293T-hACE2 (BEI #NR-52511, HEK-293T cells constitutively expressing ACE2) for 60 h and were analyzed for GFP-positive cells using a BD LSRFortessa Cell Analyzer. For negative stain experiments, commercially obtained or home-made virus concentrated using ultracentrifugation on a 20% sucrose layer⁷¹ was used.

Negative staining of pseudotyped SARS-CoV-2 S lentiviral particle

10 µl of concentrated Pseudotyped SARS-CoV-2 S lentiviral particles were incubated for 30 min on carbon coated grids (EMS #CF200-CU). Excess solution was removed by back blotting, the sample was washed three times in a drop of 1×PBS and H₂O, subsequently stained with 5 µl of 2% (w/v) uranyl acetate solution. The grids were imaged using a Tecnai T12 transmission electron microscope (FEI) operated at 120 keV.

Single particle analysis of SARS-Cov-2 S protein, data collection and image analysis

SARS-CoV-2 S protein (Cube Biotech) was recorded at 0.3 mg/ml. 3 µl of sample was applied on glow discharged (Pelco easiGlow^TM; t = 20 s; I = 20 mA) 200-mesh R1.2/1.3 Quantifoil girds. After blotting for 3.5 s at 4 °C and 100% humidity, the sample was vitrified in liquid ethane using a Vitrobot Mark IV.

The data was acquired on a Glacios (Thermo Fisher Scientific) operated at 200 keV and equipped with a Falcon 4 direct electron detector. Images were collected by EPU software (Thermo Fisher Scientific) with a pixel size of 0.93 Å with a defocus range from −0.8 to −2.4 µm. In total, 3060 movies, divided in 40 frames, were collected with a total dose of 49 e⁻/Ų.

All data processing steps were performed in cryoSAPRC v.3.3.1⁷². Motion correction in patch mode, CTF estimation in patch mode and subsequent blob picking were performed. An initial set of obtained particles was used for training a Topaz which was optimized over several rounds to extract the final set of 11,549 particles⁷³. Heterogenous refinement was performed to separate open and closed states of S proteins. The closed-state S protein class was further refined using C3 symmetry. The final map was reconstructed by non-uniform refinement with per particle CTF estimation and aberration correction from 7618 particles^74,75 at a resolution of 3.56 Å with the FSC 0.143 criterion. The final map was sharpened using DeepEMhancer⁷⁶ and visualized in ChimeraX⁷⁷. For the reconstruction of S protein in the open conformation, the same analysis scheme was applied using C1 symmetry on a final set of 3931 particles. The final resolution was estimated to be 7.44 Å with the FSC 0.143 criterion.

Ligand binding ELISA assay

Ligand binding assays were performed using a published protocol⁷⁸ with modifications described below. Solutions of human plasma vitronectin (5 µg/ml), human fibrinogen (20 µg/ml), bovine plasma fibronectin (10 µg/ml), S protein (residues 1–1212 with stabilizing mutations R684G, R685S, R678G, K998P, and V999P) with a C-terminal foldon motif and a His6 tag, 20 µg/ml)⁷⁹ and SARS-CoV-2 S1(residues 1-696 with a C-terminal foldon motif and a His6 tag, 20 µg/ml) in TBS were used to coat 96-well polyvinylchloride microtiter plates (Nunc Maxisorp #44-2404-21) for 6 h at RT. Non-coated well was used to determine unspecific binding background values. After blocking overnight at 4 °C (Blocking One, Nacalai #03953-95), velcro-tagged integrins α_vβ₃, α_IIbβ₃, and α₅β₁ diluted at 10 µg/ml in 20 mM Hepes, 150 mM NaCl, 1 mM MnCl₂, pH 7.2 were incubated for 2 h at RT to allow binding to the immobilized ligands. Bound integrins were quantified by an enzyme-linked immunosorbent (ELISA)-like solid-phase assay using biotinylated rabbit anti-velcro (against ACIS/BASE coiled-coil) antibody and HRP-conjugated streptavidin (VECTOR Laboratories #SA-5004). After addition of ABTS, readout was performed at 405 nm.

Statistics and reproducibility

The differential inference contrast microscopy, resulting in Figs. 1A–C and 2A were repeated three times for fibronectin/poly-l-lysine and four times for collagen I. The proplatelet-shapes in Supplementary Fig. 1A were observed in all spike-induced DIC experiments. The platelet shape under cryo-EM conditions in Fig. 3-G and Supplementary Fig. 2A–B were compared over a whole grid. The particles shown in Fig. 4I originate from the tomogram used for the membrane curvature analysis. The immunoblots in Supplementary Fig. 1D was performed using platelets from three different donors. S-protein pseudotyped viral particles were checked for their morphology on three different grids using negative stain. The S-protein pseudotyped viral particle in cryo-EM was observed for one grid.

Reporting summary

Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.

Data availability

The data that support this study are available from the corresponding authors upon reasonable request. Tomograms of platelets in the presence of SARS-CoV-2 S protein used in the figures were deposited to the Electron Microscopy Data Bank (EMDB) with accession codes EMD-26794 (platelet protrusion shown in Fig. 3) and EMD-26796 (platelet protrusion shown in Fig. 4). Additional tomograms used for subtomogram averaging were deposited to the Electron Microscopy Public Image Archive (EMPIAR) with accession code EMPIAR-11038 [10.6019/EMPIAR-10921]. The 3D map of the single particle reconstructed S protein has been deposited with the accession code EMD-26798. The 3D map of the 7.44 Å dataset is available upon reasonable request. The atomic model of the SARS-CoV-2 spike glycoprotein in the closed state 6VXX [https://doi.org/10.2210/pdb6VXX/pdb] and open state 6VYB [https://doi.org/10.2210/pdb6VYB/pdb] were used for comparison. Data underlying Figs. 1–5, and Supplementary Fig. 1 are available as a Source Data file. Source data are provided with this paper.

References

1. Zhu, N. et al. A Novel Coronavirus from patients with pneumonia in China, 2019. N. Engl. J. Med. 382, 727–733 (2020).

   Article  CAS  Google Scholar 

2. Wu, F. et al. A new coronavirus associated with human respiratory disease in China. Nature 579, 265–269 (2020).

   Article  ADS  CAS  Google Scholar 

3. Coronaviridae Study Group of the International Committee on Taxonomy of, V. The species Severe acute respiratory syndrome-related coronavirus: classifying 2019-nCoV and naming it SARS-CoV-2. Nat. Microbiol. 5, 536–544 (2020).

   Article  Google Scholar 

4. Lu, R. et al. Genomic characterisation and epidemiology of 2019 novel coronavirus: implications for virus origins and receptor binding. Lancet 395, 565–574 (2020).

   Article  CAS  Google Scholar 

5. Zhou, P. et al. A pneumonia outbreak associated with a new coronavirus of probable bat origin. Nature 579, 270–273 (2020).

   Article  ADS  CAS  Google Scholar 

6. Guo, Y. R. et al. The origin, transmission and clinical therapies on coronavirus disease 2019 (COVID-19) outbreak - an update on the status. Mil. Med. Res. 7, 11 (2020).

   CAS  Google Scholar 

7. Huang, C. et al. Clinical features of patients infected with 2019 novel coronavirus in Wuhan, China. Lancet 395, 497–506 (2020).

   Article  CAS  Google Scholar 

8. Al-Samkari, H. et al. COVID-19 and coagulation: bleeding and thrombotic manifestations of SARS-CoV-2 infection. Blood 136, 489–500 (2020).

   Article  CAS  Google Scholar 

9. Chen, N. et al. Epidemiological and clinical characteristics of 99 cases of 2019 novel coronavirus pneumonia in Wuhan, China: a descriptive study. Lancet 395, 507–513 (2020).

   Article  CAS  Google Scholar 

10. Li, L. et al. Analysis of viral load in different specimen types and serum antibody levels of COVID-19 patients. J. Transl. Med. 19, 30 (2021).

    Article  CAS  Google Scholar 

11. Fajnzylber, J. et al. SARS-CoV-2 viral load is associated with increased disease severity and mortality. Nat. Commun. 11, 5493 (2020).

    Article  ADS  CAS  Google Scholar 

12. Brasen, C. L. et al. Daily monitoring of viral load measured as SARS-CoV-2 antigen and RNA in blood, IL-6, CRP and complement C3d predicts outcome in patients hospitalized with COVID-19. Clin. Chem. Lab Med. 59, 1988–1997 (2021).

    Article  CAS  Google Scholar 

13. Chen, L. et al. Dynamics of blood viral load is strongly associated with clinical outcomes in coronavirus disease 2019 (COVID-19) patients: a prospective cohort study. J. Mol. Diagn. 23, 10–18 (2021).

    Article  Google Scholar 

14. Lippi, G., Plebani, M. & Henry, B. M. Thrombocytopenia is associated with severe coronavirus disease 2019 (COVID-19) infections: a meta-analysis. Clin. Chim. Acta 506, 145–148 (2020).

    Article  CAS  Google Scholar 

15. Tang, N., Li, D., Wang, X. & Sun, Z. Abnormal coagulation parameters are associated with poor prognosis in patients with novel coronavirus pneumonia. J. Thromb. Haemost. 18, 844–847 (2020).

    Article  CAS  Google Scholar 

16. Manne, B. K. et al. Platelet gene expression and function in patients with COVID-19. Blood 136, 1317–1329 (2020).

    Article  CAS  Google Scholar 

17. Lodigiani, C. et al. Venous and arterial thromboembolic complications in COVID-19 patients admitted to an academic hospital in Milan, Italy. Thromb. Res. 191, 9–14 (2020).

    Article  CAS  Google Scholar 

18. Middleton, E. A. et al. Neutrophil extracellular traps contribute to immunothrombosis in COVID-19 acute respiratory distress syndrome. Blood 136, 1169–1179 (2020).

    Article  CAS  Google Scholar 

19. Zaid, Y. et al. Platelets can associate with SARS-Cov-2 RNA and are hyperactivated in COVID-19. Circ. Res. 127, 1404–1418 (2020).

    Article  CAS  Google Scholar 

20. Zuo, Y. et al. Prothrombotic autoantibodies in serum from patients hospitalized with COVID-19. Sci. Transl. Med. 12, eabd3876 (2020).

    Article  CAS  Google Scholar 

21. Zhang, S. et al. SARS-CoV-2 binds platelet ACE2 to enhance thrombosis in COVID-19. J. Hematol. Oncol. 13, 120 (2020).

    Article  Google Scholar 

22. Koupenova, M. et al. SARS-CoV-2 initiates programmed cell death in platelets. Circ. Res. 129, 631–646 (2021).

    Article  CAS  Google Scholar 

23. Puhm, F. et al. Platelet activation by SARS-CoV-2 implicates the release of active tissue factor by infected cells. Blood Adv. 6, 3593–3605 (2022).

    Article  CAS  Google Scholar 

24. Weissenhorn, W. et al. Structural basis for membrane fusion by enveloped viruses. Mol. Membr. Biol. 16, 3–9 (1999).

    Article  CAS  Google Scholar 

25. Wrapp, D. et al. Cryo-EM structure of the 2019-nCoV spike in the prefusion conformation. Science 367, 1260–1263 (2020).

    Article  ADS  CAS  Google Scholar 

26. Coutard, B. et al. The spike glycoprotein of the new coronavirus 2019-nCoV contains a furin-like cleavage site absent in CoV of the same clade. Antivir. Res. 176, 104742 (2020).

    Article  CAS  Google Scholar 

27. Johnson, B. A. et al. Loss of furin cleavage site attenuates SARS-CoV-2 pathogenesis. Nature 591, 293–299 (2021).

    Article  ADS  CAS  Google Scholar 

28. Wan, Y., Shang, J., Graham, R., Baric, R. S. & Li, F. Receptor recognition by the novel coronavirus from wuhan: an analysis based on decade-long structural studies of SARS coronavirus. J. Virol. 94, e00127–20 (2020).

    Article  Google Scholar 

29. Hoffmann, M. et al. SARS-CoV-2 cell entry depends on ACE2 and TMPRSS2 and is blocked by a clinically proven protease inhibitor. Cell 181, 271–280.e278 (2020).

    Article  CAS  Google Scholar 

30. Walls, A. C. et al. Structure, function, and antigenicity of the SARS-CoV-2 spike glycoprotein. Cell 181, 281–292.e286 (2020).

    Article  CAS  Google Scholar 

31. Daly, J. L. et al. Neuropilin-1 is a host factor for SARS-CoV-2 infection. Science 370, 861–865 (2020).

    Article  ADS  CAS  Google Scholar 

32. Chen, Z. et al. Function of HAb18G/CD147 in invasion of host cells by severe acute respiratory syndrome coronavirus. J. Infect. Dis. 191, 755–760 (2005).

    Article  CAS  Google Scholar 

33. Wang, K. et al. CD147-spike protein is a novel route for SARS-CoV-2 infection to host cells. Signal Transduct. Target Ther. 5, 283 (2020).

    Article  CAS  Google Scholar 

34. Schumacher, S. et al. Structural insights into integrin alpha5beta1 opening by fibronectin ligand. Sci. Adv. 7, eabe9716 (2021).

    Article  ADS  CAS  Google Scholar 

35. Ruoslahti, E. RGD and other recognition sequences for integrins. Annu Rev. Cell Dev. Biol. 12, 697–715 (1996).

    Article  CAS  Google Scholar 

36. Hu, D. et al. Genomic characterization and infectivity of a novel SARS-like coronavirus in Chinese bats. Emerg. Microbes Infect. 7, 154 (2018).

    Article  Google Scholar 

37. Simons, P. et al. Integrin activation is an essential component of SARS-CoV-2 infection. Sci. Rep. 11, 20398 (2021).

    Article  ADS  CAS  Google Scholar 

38. Nader, D., Fletcher, N., Curley, G. F. & Kerrigan, S. W. SARS-CoV-2 uses major endothelial integrin alphavbeta3 to cause vascular dysregulation in-vitro during COVID-19. PLoS ONE 16, e0253347 (2021).

    Article  CAS  Google Scholar 

39. Beddingfield, B. J. et al. The integrin binding peptide, ATN-161, as a novel therapy for SARS-CoV-2 infection. JACC Basic Transl. Sci. 6, 1–8 (2021).

    Article  Google Scholar 

40. Park, E. J. et al. The spike glycoprotein of SARS-CoV-2 binds to beta1 integrins expressed on the surface of lung epithelial cells. Viruses 13, 645 (2021).

    Article  CAS  Google Scholar 

41. Schwertz, H. et al. Anucleate platelets generate progeny. Blood 115, 3801–3809 (2010).

    Article  CAS  Google Scholar 

42. Thon, J. N. et al. Cytoskeletal mechanics of proplatelet maturation and platelet release. J. Cell Biol. 191, 861–874 (2010).

    Article  CAS  Google Scholar 

43. Thon, J. N. et al. Microtubule and cortical forces determine platelet size during vascular platelet production. Nat. Commun. 3, 852 (2012).

    Article  ADS  Google Scholar 

44. Smith, C. W. Release of alpha-granule contents during platelet activation. Platelets 33, 491–502 (2022).

    Article  CAS  Google Scholar 

45. Guidetti, G. F., Torti, M. & Canobbio, I. Focal adhesion kinases in platelet function and thrombosis. Arterioscler Thromb. Vasc. Biol. 39, 857–868 (2019).

    Article  CAS  Google Scholar 

46. Ke, Z. et al. Structures and distributions of SARS-CoV-2 spike proteins on intact virions. Nature 588, 498–502 (2020).

    Article  ADS  CAS  Google Scholar 

47. Sorrentino, S. et al. Structural analysis of receptors and actin polarity in platelet protrusions. Proc. Natl Acad. Sci. USA 118, e2105004118 (2021).

    Article  CAS  Google Scholar 

48. Mullins, R. D., Heuser, J. A. & Pollard, T. D. The interaction of Arp2/3 complex with actin: nucleation, high affinity pointed end capping, and formation of branching networks of filaments. Proc. Natl Acad. Sci. USA 95, 6181–6186 (1998).

    Article  ADS  CAS  Google Scholar 

49. Martinez-Sanchez, A. et al. Template-free detection and classification of membrane-bound complexes in cryo-electron tomograms. Nat. Methods 17, 209–216 (2020).

    Article  CAS  Google Scholar 

50. Turonova, B. et al. In situ structural analysis of SARS-CoV-2 spike reveals flexibility mediated by three hinges. Science 370, 203–208 (2020).

    Article  ADS  CAS  Google Scholar 

51. Crawford, K. H. D. et al. Protocol and reagents for pseudotyping lentiviral particles with SARS-CoV-2 spike protein for neutralization assays. Viruses 12, 513 (2020).

    Article  CAS  Google Scholar 

52. Reardon, D. A., Nabors, L. B., Stupp, R. & Mikkelsen, T. Cilengitide: an integrin-targeting arginine-glycine-aspartic acid peptide with promising activity for glioblastoma multiforme. Expert Opin. Investig. Drugs 17, 1225–1235 (2008).

    Article  CAS  Google Scholar 

53. Burkhart, J. M. et al. The first comprehensive and quantitative analysis of human platelet protein composition allows the comparative analysis of structural and functional pathways. Blood 120, e73–e82 (2012).

    Article  CAS  Google Scholar 

54. Hadid, T., Kafri, Z. & Al-Katib, A. Coagulation and anticoagulation in COVID-19. Blood Rev. 47, 100761 (2021).

    Article  CAS  Google Scholar 

55. Fajgenbaum, D. C. & June, C. H. Cytokine storm. N. Engl. J. Med. 383, 2255–2273 (2020).

    Article  CAS  Google Scholar 

56. Schindelin, J. et al. Fiji: an open-source platform for biological-image analysis. Nat. Methods 9, 676–682 (2012).

    Article  CAS  Google Scholar 

57. Hagen, W. J. H., Wan, W. & Briggs, J. A. G. Implementation of a cryo-electron tomography tilt-scheme optimized for high resolution subtomogram averaging. J. Struct. Biol. 197, 191–198 (2017).

    Article  Google Scholar 

58. Zheng, S. Q. et al. MotionCor2: anisotropic correction of beam-induced motion for improved cryo-electron microscopy. Nat. Methods 14, 331–332 (2017).

    Article  CAS  Google Scholar 

59. Kremer, J. R., Mastronarde, D. N. & McIntosh, J. R. Computer visualization of three-dimensional image data using IMOD. J. Struct. Biol. 116, 71–76 (1996).

    Article  CAS  Google Scholar 

60. Tegunov, D. & Cramer, P. Real-time cryo-electron microscopy data preprocessing with Warp. Nat. Methods 16, 1146–1152 (2019).

    Article  CAS  Google Scholar 

61. Zhang, K. Gctf: Real-time CTF determination and correction. J. Struct. Biol. 193, 1–12 (2016).

    Article  ADS  CAS  Google Scholar 

62. Xiong, Q., Morphew, M. K., Schwartz, C. L., Hoenger, A. H. & Mastronarde, D. N. CTF determination and correction for low dose tomographic tilt series. J. Struct. Biol. 168, 378–387 (2009).

    Article  Google Scholar 

63. Zivanov, J. et al. New tools for automated high-resolution cryo-EM structure determination in RELION-3. Elife 7, e42166 (2018).

    Article  Google Scholar 

64. Sousbie, T. The persistent cosmic web and its filamentary structure – I. Theory and implementation. Monthly Not. R. Astronomical Soc. 414, 350–383 (2011).

    Article  ADS  Google Scholar 

65. Forman, R. A user’s guide to discrete Morse theory. Séminaire Lotharingien de. Combinatoire [Electron. only] 48, B48c (2002).

    MathSciNet  MATH  Google Scholar 

66. Nedozralova, H. et al. In situ cryo-electron tomography reveals local cellular machineries for axon branch development. J. Cell Biol. 221, e202106086 (2022).

    Article  CAS  Google Scholar 

67. Salfer, M., Collado, J. F., Baumeister, W., Fernandez-Busnadiego, R. & Martinez-Sanchez, A. Reliable estimation of membrane curvature for cryo-electron tomography. PLoS Comput. Biol. 16, e1007962 (2020).

    Article  ADS  CAS  Google Scholar 

68. Qu, K. et al. Structure and architecture of immature and mature murine leukemia virus capsids. Proc. Natl Acad. Sci. USA 115, E11751–E11760 (2018).

    Article  ADS  CAS  Google Scholar 

69. Pettersen, E. F. et al. UCSF Chimera-a visualization system for exploratory research and analysis. J. Comput. Chem. 25, 1605–1612 (2004).

    Article  CAS  Google Scholar 

70. Bodakuntla, S., Janke, C. & Magiera, M. M. Knocking out multiple genes in cultured primary neurons to study tubulin posttranslational modifications. Methods Mol. Biol. 2101, 327–351 (2020).

    Article  CAS  Google Scholar 

71. Tiscornia, G., Singer, O. & Verma, I. M. Production and purification of lentiviral vectors. Nat. Protoc. 1, 241–245 (2006).

    Article  CAS  Google Scholar 

72. Punjani, A., Rubinstein, J. L., Fleet, D. J. & Brubaker, M. A. cryoSPARC: algorithms for rapid unsupervised cryo-EM structure determination. Nat. Methods 14, 290–296 (2017).

    Article  CAS  Google Scholar 

73. Bepler, T. et al. Positive-unlabeled convolutional neural networks for particle picking in cryo-electron micrographs. Nat. Methods 16, 1153–1160 (2019).

    Article  CAS  Google Scholar 

74. Punjani, A., Zhang, H. & Fleet, D. J. Non-uniform refinement: adaptive regularization improves single-particle cryo-EM reconstruction. Nat. Methods 17, 1214–1221 (2020).

    Article  CAS  Google Scholar 

75. Zivanov, J., Nakane, T. & Scheres, S. H. W. Estimation of high-order aberrations and anisotropic magnification from cryo-EM data sets in RELION-3.1. IUCrJ 7, 253–267 (2020).

    Article  CAS  Google Scholar 

76. Sanchez-Garcia, R. et al. DeepEMhancer: a deep learning solution for cryo-EM volume post-processing. Commun. Biol. 4, 874 (2021).

    Article  Google Scholar 

77. Pettersen, E. F. et al. UCSF ChimeraX: Structure visualization for researchers, educators, and developers. Protein Sci. 30, 70–82 (2021).

    Article  CAS  Google Scholar 

78. Miyazaki, N., Iwasaki, K. & Takagi, J. A systematic survey of conformational states in beta1 and beta4 integrins using negative-stain electron microscopy. J. Cell Sci. 131, jcs216754 (2018).

    Article  Google Scholar 

79. Higuchi, Y. et al. Engineered ACE2 receptor therapy overcomes mutational escape of SARS-CoV-2. Nat. Commun. 12, 3802 (2021).

    Article  ADS  CAS  Google Scholar 

Download references

Acknowledgements

We thank the members of Mizuno Lab for discussions and help for this research. We thank William Wan for their insightful suggestions and advice for the tomographic analysis and membrane curvature quantification. We thank Ana Pasapera for the technical advice of signaling assays. The light microscopic data was collected at the NHLBI light microscopy core facility, the flow analysis was performed at the NHLBI flow cytometry core and the cryo-ET data was collected at the MICEF at the National Institutes of Health, USA and at the National Cancer Institute’s National Cryo-EM Facility at the Frederick National Laboratory for Cancer Research under contract HSSN261200800001E. This work used the computational resources of the NIH HPC Biowulf cluster (http://hpc.nih.gov). The following reagents were obtained through BEI Resources, NIAID, NIH: SARS-Related Coronavirus 2, Wuhan-Hu-1 Spike-Pseudotyped Lentivirus, Luc2/ZsGreen, NR-53818. N.M. acknowledges the Intramural Research Program of the National Heart Lung and Blood Institute (1ZIAHL006264 and 1ZIAHL006265), and the National Institute of Arthritis and Musculoskeletal and Skin Diseases of National Institutes of Health, USA, for funding.

Funding

Open Access funding provided by the National Institutes of Health (NIH).

Author information

Author notes

1. These authors contributed equally: Nirakar Basnet, Satish Bodakuntla, Pelayo Alvarez-Brecht.

Authors and Affiliations

1. Laboratory of Structural Cell Biology, National Heart, Lung, and Blood Institute, National Institutes of Health, 50 South Dr., Bethesda, MD, 20892, USA

   Christopher Cyrus Kuhn, Nirakar Basnet, Satish Bodakuntla, Pelayo Alvarez-Brecht, Scott Nichols, Lorenzo Agostini, Young-Min Soh, Christian Biertümpfel & Naoko Mizuno

2. Department of Computer Sciences, Faculty of Sciences - Campus Llamaquique, University of Oviedo, Oviedo, 33007, Spain

   Pelayo Alvarez-Brecht & Antonio Martinez-Sanchez

3. Health Research Institute of Asturias (ISPA), Avenida Hospital Universitario s/n, 33011, Oviedo, Spain

   Antonio Martinez-Sanchez

4. Osaka University Institute for Protein Research, Osaka University, 3-2 Yamadaoka, Suita, Osaka, 565-0871, Japan

   Junichi Takagi

5. National Institute of Arthritis and Musculoskeletal and Skin Diseases, National Institutes of Health, 50 South Dr., Bethesda, MD, 20892, USA

   Naoko Mizuno

Contributions

C.C.K., C.B., and N.M. conceived and supervised the project. C.C.K. isolated the platelets and prepared samples for the experiments. C.C.K. and N.B. collected the cryo-EM data. C.C.K., N.B., P.A.B., and A.M.S. analyzed the cryo-EM data with the assist of L.A. C.C.K. and S.N., with the help of S.B., analyzed the light microscopy data. S.B., with the help of C.C.K., produced pseudo-typed viral particles and performed the immunoblots. C.C.K., J.T., and C.B. performed and analyzed the ligand binding assays with the initial assessments by Y.M.S. S.B. and C.C.K. performed the negative staining of lentiviral particles with the help of N.M. S.B. performed the flow cytometry analysis. C.C.K., S.B., C.B., and N.M. wrote the initial draft, all authors discussed and commented on the manuscript.

Corresponding author

Correspondence to Naoko Mizuno.

Ethics declarations

Competing interests

The authors declare no competing interests.

Peer review

Peer review information

Nature Communications thanks Tione Buranda and the other, anonymous, reviewer(s) for their contribution to the peer review of this work.

Additional information

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Supplementary Information

Description of Additional Supplementary Files

Supplementary Movie 1

Supplementary Movie 2

Reporting Summary

Source data

Source Data

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

About this article

Cite this article

Kuhn, C.C., Basnet, N., Bodakuntla, S. et al. Direct Cryo-ET observation of platelet deformation induced by SARS-CoV-2 spike protein. Nat Commun 14, 620 (2023). https://doi.org/10.1038/s41467-023-36279-5

Download citation

• Received24 May 2022

• Accepted23 January 2023

• Published04 February 2023

• DOIhttps://doi.org/10.1038/s41467-023-36279-5

Share this article

Anyone you share the following link with will be able to read this content:

Provided by the Springer Nature SharedIt content-sharing initiative

Subjects

• Cell polarity
• Cell signalling
• Cryoelectron tomography
• Proteins