key: cord-0963866-vmyd7ur1 authors: Lobo, Vanessa R.; Warwicker, Jim title: Predicted pH-dependent stability of SARS-CoV-2 spike protein trimer from interfacial acidic groups date: 2021-09-02 journal: Comput Struct Biotechnol J DOI: 10.1016/j.csbj.2021.08.049 sha: 04273f28b1dc64e12df3730c3f6f5826ed577b97 doc_id: 963866 cord_uid: vmyd7ur1 Transition between receptor binding domain (RBD) up and down forms of the SARS-CoV-2 spike protein trimer is coupled to receptor binding and is one route by which variants can alter viral properties. It is becoming apparent that key roles in the transition are played by pH and a more compact closed form, termed locked. Calculations of pH-dependence are made for a large set of spike trimers, including locked form trimer structures that have recently become available. Several acidic sidechains become sufficiently buried in the locked form to give a predicted pH-dependence in the mild acidic range, with stabilisation of the locked form as pH reduces from 7.5 to 5, consistent with emerging characterisation by cryo-electron microscopy. The calculated pH effects in pre-fusion spike trimers are modulated mainly by aspartic acid residues, rather than the more familiar histidine role at mild acidic pH. These acidic sidechains are generally surface located and weakly interacting when not in a locked conformation. According to this model, their replacement (perhaps with asparagine) would remove the pH-dependent destabilisation of locked spike trimer conformations, and increase their recovery at neutral pH. This would provide an alternative or supplement to the insertion of disulphide linkages for stabilising spike protein trimers, with potential relevance for vaccine design. The spike glycoprotein (S protein) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has received much attention, as the primary point of attachment between virus and ACE2 receptor on a cell surface [1] , and as a target in vaccine design [2] . An original focus on the pre-fusion S protein trimer structure in closed and open states, with receptor binding domain (RBD) down and up, respectively, has expanded to include a further RBD down state, termed locked [3] , first observed with in a structure with linoleic acid bound [4] . Here, the terms up and down are used to delineate RBD position, so that the open form maps to RBD up, but both closed and locked forms map to RBD down. This 3 state pre-fusion S trimer picture will be expanded with greater characterisation of post-fusion trimers [5] , but the transitions between pre-fusion trimers are of key interest for the early steps in cell entry, and also for the export of fusion competent virus. Further, spike protein trimer expressed at the cell surface is able to mediate fusion between cells [6] . The role of pH has been studied for cell entry and exit of coronaviruses. Both plasma membrane fusion and endosomal mediated fusion pathways can be used [7] , with fusion primed by S1/S2 spike protein cleavage during the biosynthetic exit pathway [8] . S1/S2 and S2/S2' cleavage sites, and the fusion peptide, are critical locations in considering spike trimer conformational transitions [9] . Most cryo-electron microscopy (cryo-EM) structures of spike trimers have been resolved at pH 7.5 -8. Studies at mild acidic pH have found structural transitions [3, 10] , interpreted as relatively favouring closed over locked forms as pH is increased from mild acidic to neutral [3] . Other cryo-EM analysis shows changes in the fraction of RBD up and down conformations between pH values of 6.5 and 8 [11] . Cell biological investigation of the exit pathway for newly synthesised SARS-CoV-2 virus has revealed a novel mechanism where lysosomes are repurposed for viral transit to the cell surface, with a moderate de-acidification of their mean pH from 4.7 to 5.7 [12] . Other, theoretical, work has looked at the wider effects of virus-induced pH change in SARS-CoV-2 infection [13] . Biophysical and structural analysis of the D614G S protein mutation that emerged in Spring 2020 indicates that it leads to greater availability of the RBD for ACE2 binding, through altered structural transitions [14] . These transitions include packing changes that resemble closed to locked form differences [15] . A model has been proposed in which the locked form is prevalent in the virus export pathway, giving protection of newly synthesised spike proteins at acidic pH, with the D614G mutation coupling to transition between locked and closed identifiers) were SARS-CoV-2 (P0DTC2), SARS (P59594), MERS (K9N5Q8), HCoV-HKU1 (Q5MQD0), HCoV-OC43 (P36334), HCoV-NL63 (Q6Q1S2), and HCoV-229E (P15423). Solvent accessible surface area (abbreviated to SASA) was calculated for S protein monomer burial within a trimer. For each monomer within a trimer (after non-spike components have been removed), SASA was calculated for the monomer alone and for the monomer within the trimer, the difference (d-SASA) giving the extent of burial. The contribution of polar and non-polar atoms to this burial was recorded, following previous methodology [35] . To facilitate comparison between sub-groups, d-SASA values were averaged over monomers within a sub-group. A key feature is that all calculations were referenced to an alignment of amino acid sequences for all monomers (spike387 and additional monomers added in subsequent formation of sub-groups). This ensures that all comparisons and heat maps are correctly aligned, allowing for engineered monomers and residues present in a sequence but disordered in a structure. Amino acid numbering throughout is given as the equivalent location in UniProt entry P0DTC2. Missing coordinates are not modelled, structurebased calculations use only the ordered residues reported in each structure. From vectors of 1 (ordered) and 0 (disordered) running along the sequence of each monomer, similarity to chain A of the locked conformation 6zb4 was calculated by vector dot product. Monomers were then listed in heat maps (of order/disorder or d-SASA) according to similarity with 6zb4A. Electrostatics calculations used pK a predictions following the Finite Difference Poisson Boltzmann (FDPB) / Debye-Hückel (DH) hybrid method, termed FD/DH [26] . Ionisable groups that are not buried can sample both FDPB and DH schemes, and these calculations are combined with Boltzmann weighting. Since the latter is a simple water-dominated and highly damped estimate of electrostatic interactions, non-buried groups will not have high pK a deviations (pK a s) from normal values, unless those are stabilising. On the other hand, buried groups, shielded from water-dominated interactions, can have pK a values that relate to substantial destabilisation as well as stabilisation. The method has been implemented in a web tool, following testing with sites known to determine the pH-dependence of influenza virus hemagglutinin stability [27] . Here, the focus is on predicted stability change between pH 7.5 (extracellular, cytoplasmic pH) and the lower pH experienced in import and export pathways (pH 5 was used as the lower value). This pH-dependent energy can be derived using the relationship [36] [37] : G pH-dep = (2.303RTQ)dpH (1) where G pH-dep is the pH-dependent contribution to conformational stability (over the pH range of the integration, 7.5 to 5), differenced between two states. The states here are ionisable groups interacting in the FD/DH scheme, and the same set without interactions (null state with normal model compound pK a s). Q is the charge difference between those two states, R is the Universal Gas constant, and T is set at 298.15 K. Since the null state is uniform between structures, results for each structure can be compared to assess differential responses to pH change in the mild acidic range. Groups responsible for the predicted effects can be identified since Q separates into component ionisations. For 5 of the 6 structural sub-groups, monomers are either all RBD down or all RBD up, and in many cases will have been refined with imposed trimer symmetry. It is still useful to calculate on all monomers to assess the errors associated with numerical parts of the calculation (use of a finite difference grid and Monte Carlo sampling of protonation states). The set of predicted pH-dependent energies for each ionisable group reproduce closely, for example with correlation coefficients of 0.99 between monomers of 6vxx, indicating self-consistency of the electrostatics calculations. An additional method was used to calculate pK a s, in order to compare with FD/DH results. The empirical calculation tool PROPKA [38] implemented in PDB2PQR at the APBS server [39] , was used (with default parameters) for the 6 sub-groups of spike proteins (36 trimers in total). Predicted pK a s from PROPKA were averaged for various single or groups of amino acids, in each of the 6 spike protein sub-groups. Molecular viewers used to analyse spike proteins were NGL [40] , Swiss PDB Viewer [41] , and PyMOL. For representative trimer structures in the locked (6zp2), closed (6zp1) and open (7a98) sub-groups, d-SASA values were transferred to the B factor column (capped at 99 Å 2 ), facilitating visualisation of monomer burial within a trimer. The ratio of non-polar to polar contributions to d-SASA was calculated for the S1 and S2 subunits of each monomer (Fig. 2b) . For external protein surfaces not involved in proteinprotein interactions the most common value of this ratio, in a distribution over patches, is about 1.2, increasing for more hydrophobic, interacting surfaces [35] . For the open sub-group, with least monomer burial, non-polarity increases for both S1 and S2 contributions. Conversely, monomer interface within a trimer becomes successively more polar for sub-groups that have higher d-SASA, with polarity for S1 d-SASA in the locked form close to that seen for noninterfacial surface (Fig. 2b) . It can be concluded that transition from RBD up to RBD down, and further transition to the locked form, involves burial of surface that is not especially hydrophobic. Comparison of Fig. 3a and Fig. 3b reinforces the observation that d-SASA differences between S trimer sub-groups are more extensively spread over spike monomers than are order/disorder differences. The extent of d-SASA change between locked and closed forms, in comparison with that between closed and open forms, is also emphasised. In light of reported pH-dependent transition between S trimer conformations [3, 10] , pK a predictions were made for ionisable residues in trimer sub-groups. To directly address changes in the mild acidic pH range, the predicted contribution to change in conformational stability between pH 7.5 (cytoplasm, extracellular) and pH 5 (endosome) was calculated (summed and for each ionisable group). The results are added to Fig. 3 as a heat map of calculated per-residue pH-dependence, averaged over monomers within each sub-group (Fig. 3c) . Towards the C-terminal part of S2 are several histidines [27, 42] , for which burial and predicted resistance to protonation from pH 7.5 to 5 leads to destabilisation in all sub-groups. Of more interest for differences between sub-groups are amino acids with acidic sidechains that have elevated pK a s and are predicted to be stabilising from pH 7.5 to 5, since destabilisation is reduced. Destabilisation derives from a local environment that suppresses ionisation at a pH above the normal pK a . Several Asp and Glu sidechains are of interest, particularly those for which the effect is exhibited in locked forms but not others, as a consequence of the increased burial. Summed values of predicted pH-dependence give a stabilisation of locked form subgroups as pH falls from 7.5 to 5, but a destabilisation of other sub-groups over the same range ( Fig. 4) . The consequence of these predictions is a model in which the locked form is preferentially stabilised in low pH secretory pathway vesicles, whereas closed (unlocked) and open forms are stabilised at extracellular, neutral pH, as suggested from structural studies [3] . Looking at individual group contributions, 3 His in S2 (H1048, H1064, H1088) are buried and largely inaccessible to solvent throughout the sub-groups, so that they give a uniform destabilisation over the pH 7.5 to 5 range, as each His transits through the pK a (6.3) that would normally mark ionisation. These groups are not predicted to contribute differentially to pH-dependence of locked, closed, and open forms. Next, D398 was proposed to couple to pH-dependence of closed (mostly unlocked) to open forms [42] . This result is replicated in the current larger study (Fig. 4) , with a predicted D398 contribution to variation of pH-dependence between locked, closed, and open forms, but in the opposite sense to the overall effect. Two histidines have been suggested as potential mediators of locked trimer stabilisation at acidic pH [3] , but no clear differential contributions of H49 and H519 are seen between sub-groups of structures, with the FD/DH calculation of pH-dependent stability (Fig. 4 ). It is of interest to examine whether structure resolution is a factor in the results shown in Fig. 4 , in particular for the key predicted differences in pH-dependence between locked and closed sub-groups. A comparison of reported resolution and segment Manders' overlap coefficients for the 36 sub-group trimers ( Supplementary Fig. S2a ) shows the general decline of SMOC values for better resolution structures that has been observed previously [43] , and interpreted as the influence of map resolution on the level of detail required for a given fit of model to map. Given this rather complex relationship, a SMOC criterion was excluded from filtering for structure quality, leaving the reported resolution. Filtering at 3.0 Å resolution gives 5 of 8 locked form sub-group trimers and 4 of 8 closed form sub-group trimers. Predicted pHdependence for these resolution-filtered structures ( Supplementary Fig. S2b ) matches those given for the full locked and closed form sub-groups (Fig. 4) , thus giving confidence in those calculations. In order to provide context for this particular calculation method, the PROPKA empirical algorithm for pK a prediction [38] was employed, as implemented at the APBS/PDB2PQR server [39] . Here the pH-dependent stability contributions of individual groups are not available, so pK a values were used, which are positive for aspartic acid and glutamic acid destabilisation and negative for histidine destabilisation ( Supplementary Fig. S3 ). Over different selections of groups, similar predicted effects are seen to those for the FD/DH method (Fig. 4) . The 3 S2 subunit His (H1048, H1064, H1088) are predicted to be substantially destabilising in all S protein sub-groups. D398 and D614 are predicted to be destabilising, with similar trends between sub-groups to those in Fig. 4 . The subset of groups according to the PROPKA method. Interestingly, H519 is in the proximity of D979, one of the highlighted acidic groups in Fig. 3c , and the mutation H49Y increases cell entry (relative to wild type), in an S-pseudotyped lentiviral system (although much less than D614G) [44] . Given these observations, it would be reasonable to keep H49 and H519 in mind when considering potential sites of pH-sensing in the S protein trimer. A feature of predicted pH-dependence for which the picture changes when locked (RBD down) structures are separated from closed (RBD down), is the influence of D614. Contributions of D614 to the calculated pH-dependence are seen for locked and pHlocked subgroups, in a sense to stabilise those forms at pH 5 relative to pH 7.5 ( Fig. 4 ). It has been reported that the S protein D614G mutation increases virus infectivity [45] , and suggested that the mutation leads to a greater population of up relative to down RBDs due to loss of an intermonomer latch formed by the sidechain of D614 [14] . Recent structural studies with an S protein ectodomain engineered with D614G, report that the increase in infectivity results from stabilisation of the S trimer against dissociation [5] . One of the structures (7krq) from that D614G study has all RBDs down, and d-SASA for monomers in this structure ( Supplementary Fig. S4 ) are similar to those seen for locked forms (Fig. 2a) , consistent with trimer stabilisation. Generally, the D614Gset sub-group exhibits a wide distribution of monomer d-SASA values, but these structures mostly have mixed up and down RBD within each trimer ( Supplementary Fig. S4 ). D614 is predicted to be relatively destabilising at neutral pH in locked and pHlocked forms (Fig. 3c, Fig. 4 ). Even after incorporating stabilising D614 sidechain interaction with K854, the reduction in solvent exposure due to ordering of the loop around residue 840 leads to a predicted average destabilisation, at pH 7.5 compared to pH 5 and averaged over the locked sub-group, of 7.7 kJ/mol. This implies a mean D614 pK a shift to about 6.0, from the intrinsic value of 4.0, consistent with wide-spread D614 -K854 salt-bridge formation in locked structures at pH 7.5 or 8. There are examples of this salt-bridge in all sub-groups (other than D614Gset), although at lower incidence outside of the locked sub-group. Our interpretation is that D614 -K854 interaction is solvent exposed and relatively weak outside of locked foms. Then, ordering of the loop around residue 840 in the locked form decreases solvent exposure, putting a greater onus on salt-bridge formation to reduce destabilisation of D614. Summed over interactions though, destabilisation of D614 is predicted to result for the locked form, and will be greater at pH 7.5 than at pH 5. This model is consistent with stabilisation of the D614G S trimer against dissociation at neutral pH [46] . Referring to the 3 structures in the pHlocked sub-group, for two (6xlu, pH 4.0 and 7jwy, pH 4.5) the 840 loop is ordered. In both of these, K854 is outside of salt-bridge range, consistent with elevated pK a and protonation of D614 at acidic pH. In common with D614, other predicted sites of pH-dependence are coupled to a tightening of monomer-monomer interfaces in the locked form S trimer ( Fig. 5a ), so that solvent exposure is reduced and Asp/Glu sidechain pK a s are elevated, with predicted destabilisation of the locked form as pH increases through the mild acidic pH range. Presumably, pH-independent features, such as linoleic acid binding [4] , also couple to the energetic balance between locked and closed forms. For the pH-dependent component, it is notable that a tendency towards a locked form at acidic pH is still evident for the D614G S protein trimer [3] , consistent with the existence of other pH-dependent sites (Fig. 3c, Fig. 5a ). Whereas the D614 -K854 interaction is between monomers, neighbouring basic residues for other acidic sites of predicted pH-dependence, where they exist, tend to be from the same monomer. Nevertheless, the predicted effect is the same, it is the solvent exposure of the acidic group rather than the details of salt-bridging that is the key determinant of predicted pH-dependence. Several of the implicated acidic residues are present in only SARS-CoV-2 and SARS, of coronaviruses known to infect humans (Fig. 5b) . Display of d-SASA on S protein structure shows the expected difference between open ( Fig. 6a ) and closed ( Fig. 6b ) forms, and emphasises the difference (that is not evident by overall RBD location) between closed (Fig. 6b ) and locked (Fig. 6c) forms. General location ( Fig. 6d ) and detail (Fig. 6e ) are shown for residues that are present in the more compact interfaces of the locked form trimer, and which are predicted to contribute to pH-dependence over the mild acidic range. Some residues (D586, D614, D830, D843) flank regions of known importance for structural transitions (S1/S2 and S2/S2' cleavage sites, fusion peptide), and other residues (D405, D420, E465, D979, D994) are adjacent to a site commonly engineered to generate more stable trimers [3] . Additionally, sites of predicted pH-dependence in the RBD are proximal to the linoleic acid binding site [4] . It is therefore suggested that members of this set of acidic residues (Fig. 5b, Fig. 6e ), beyond D614, couple pH-dependence to conformational biases of functional relevance. In this model for pH-dependence of locked to closed form transition, a surface acidic group, whether or not paired with a basic group, loses substantial SASA on formation or tightening of an interface. The acidic pK a increases in the range pH 4 to 7. A bias against interface formation is introduced, which is relieved as pH reduced. Depending on the conformational details and interactions with neighbouring groups, a buried acidic group could be stabilising, but that is not predicted here. Although functional ionisation of carboxylate sidechain groups at mild acidic pH is most well-known for acid-base catalytic systems, such as the elevated pK a of E35 in hen egg white lysozyme [47] , many non-catalytic occurrences have been discussed. In viruses, coiledcoil formation in a variant influenza hemagglutinin is pH-dependent between pH 4.5 and 7.1, possibly mediated by E69 and E74 [48] . A model for pH-dependent association of some antigenic peptides and class II proteins of the major histocompatibility complex suggests that buried interfacial Asp or Glu residues mediate this process, with pK a s higher than 7 [49] . Interfacial residues D612 and E613 of PsaB in photosystem I are proposed to be responsible for mild acidic pH-dependence of electron transfer in complexes with plastocyanin or cytochrome c 6 [50] . The periplasmic chaperone HdeA of Escherichia coli acts on a shift from neutral to lower pH via dimer dissociation and partial unfolding, these changes mediated by two aspartic acid residues [51] . In a further example, aspartic acid residues have been proposed as determining neutral to low pH conformational variation in an antigen-binding region of a therapeutic monoclonal antibody [52] . An acidic pH-dependent stability has been engineered into staphylococcal nuclease with the mutation V66E, burying a glutamic acid (measured pK a 8.8) in the protein interior [53] . Here, the unfavourable stability consequence of glutamic acid burial is absorbed by the stability of the enzyme fold, with the addition of water molecules that penetrate to form a small channel around V66E [53] . For spike protein, it is hypothesised that unfavourable burial of acidic groups could also be coupled to other conformational factors (that are energetically favourable), including formation of the tighter locked form interfaces and linoleic acid binding. In general terms, there are some parallels with conformational transitions and pocket factor binding in picornaviruses [54] . Given the background of acidic sidechains mediating mild acidic pH-dependent effects, it is reasonable to investigate similar interactions when considering the pH-dependence of spike protein conformation, supplementing the assessment of histidine protonation in regard to pHinduced changes in viruses [55] . Analysis of spike trimer structures shows the expected large SASA change, for monomer burial, between closed and open forms, and an even larger change and interface tightening between closed and locked forms. A subset of acidic groups was identified in the locked form monomer-monomer interfaces that are predicted to mediate pH-dependent stability at mild acidic pH. These are located in regions that flank sites of known importance for stability engineering and/or function, such as cleavage sites and fusion peptide. D614 is amongst the groups identified and, in common with other acidic groups in the subset, is predicted to destabilise the locked form at neutral pH. This result is consistent with a report that the G614 mutation can exist in a form that resembles locked and stabilises the trimer at neutral pH, leading to suppression of S1 subunit loss and increased RBD -ACE2 binding [46] . The current work is also consistent with a model for pH-dependent conformational masking of RBDs within the S protein trimer, derived from the 3 structures of the pHlocked set [10] . That study highlighted D614 and D586 (in common with our calculations), and also D574. If the prediction of a wider set of acidic groups relevant for pH-dependence is correct, then it might be asked why, from this set, only mutation at D614 has been seen widely seen in SARS-CoV-2 genomes. Of note though, a study of recurring spike protein mutations highlights several regions that overlap with, or are adjacent to, segments identified in Fig. 6e [56] . These mutations could influence folding and packing around the sites predicted to mediate pHdependence. Since pH modulation of the D614G spike trimer conformation has been reported [3] , there must be groups beyond D614 responsible for some degree of pH-coupling. The recently discovered lysosomal exit pathway for newly synthesised SARs-CoV-2 virus provides a mild acidic environment, even including the observed moderate deacidification [12] . Reduced acidification is thought to reduce activity of lysosomal proteases, thereby protecting exposed viral proteins from proteolysis, beyond the S1/S2 priming cleavage of S protein that confers an advantage for SARS-CoV-2 infection [57] . Stabilising the locked form at acidic pH is proposed to contribute to protection of the spike protein along the secretory pathway [3] . The current results are consistent with this model, extending the set of groups that may be responsible for pH-dependent stabilisation (beyond D614). Importantly, the balance between priming cleavage and stabilisation along the biosynthetic pathway is the focus for ongoing viral adaptation, including the observation that P681R mutation in the B.1.617 lineage optimises the furin cleavage site and could enhance transmissibility and pathogenicity [58] . An implication of the current work is that it should be possible to engineer more stable locked forms at neutral pH through mutation of the acidic groups, (e.g. Asp to Asn), at sites indicated in Fig. 6e . In line with the observation of close packing around G614, Asp to Gly changes could be tested, although the D614N mutation also leads to a more compact spike trimer [59] . Such analysis would probe the balance of stabilising and destabilising contributions to spike trimer structure during trafficking, and it could provide further insight into production of specific spike protein conformational forms for vaccines. Support from the UK Engineering and Physical Sciences Research Council (award EP/N024796/1) is gratefully acknowledged. showing the ratio of non-polar to polar buried d-SASA for the S1 and S2 subunits, upon monomer burial in the trimer. Variation of this ratio is much greater for S1 than for S2 (seen with equal scaling of the axes), with higher values for the open form indicating that residual buried surface with RBD up is relatively non-polar. A-531-T 6ZGE-ABC-A-532-N 6ZGE-ABC-A-533-L 6ZGE-ABC-A-534-V 6ZGE-ABC-A-535-K 6ZGE-ABC-A-536-N 6ZGE-ABC-A-537-K 6ZGE-ABC-A-538-C 6ZGE-ABC-A-539-V 6ZGE-ABC-A-540-N 6ZGE-ABC-A-541-F 6ZGE-ABC-A-542-N 6ZGE-ABC-A-543-F 6ZGE-ABC-A-544-N 6ZGE-ABC-A-545-G 6ZGE-ABC-A-546-L 6ZGE-ABC-A-547-T 6ZGE-ABC-A-548-G 6ZGE-ABC-A-549-T 6ZGE-ABC-A-550-G 6ZGE-ABC-A-551-V 6ZGE-ABC-A-552-L 6ZGE-ABC-A-553-T 6ZGE-ABC-A-554-E 6ZGE-ABC-A-555-S 6ZGE-ABC-A-556-N 6ZGE-ABC-A-557-K 6ZGE-ABC-A-558-K 6ZGE-ABC-A-559-F 6ZGE-ABC-A-560-L 6ZGE-ABC-A-561-P 6ZGE-ABC-A-562-F 6ZGE-ABC-A-563-Q 6ZGE-ABC-A-564-Q 6ZGE-ABC-A-565-F 6ZGE-ABC-A-566-G 6ZGE-ABC-A-567-R 6ZGE-ABC-A-568-D 6ZGE-ABC-A-569-I 6ZGE-ABC-A-570-A 6ZGE-ABC-A-571-D 6ZGE-ABC-A-572-T 6ZGE-ABC-A-573-T 6ZGE-ABC-A-574-D 6ZGE-ABC-A-575-A 6ZGE-ABC-A-576-V 6ZGE-ABC-A-577-R 6ZGE-ABC-A-578-D 6ZGE-ABC-A-579-P 6ZGE-ABC-A-580-Q 6ZGE-ABC-A-581-T 6ZGE-ABC-A-582-L 6ZGE-ABC-A-583-E 6ZGE-ABC-A-584-I 6ZGE-ABC-A-585-L 6ZGE-ABC-A-586-D 6ZGE-ABC-A-587-I 6ZGE-ABC-A-588-T 6ZGE-ABC-A-589-P 6ZGE-ABC-A-590-C 6ZGE-ABC-A-591-S 6ZGE-ABC-A-592-F 6ZGE-ABC-A-593-G 6ZGE-ABC-A-594-G 6ZGE-ABC-A-595-V 6ZGE-ABC-A-596-S 6ZGE-ABC-A-597-V 6ZGE-ABC-A-598-I 6ZGE-ABC-A-599-T 6ZGE-ABC-A-600-P 6ZGE-ABC-A-601-G 6ZGE-ABC-A-602-T 6ZGE-ABC-A-603-N 6ZGE-ABC-A-604-T 6ZGE-ABC-A-605-S 6ZGE-ABC-A-606-N 6ZGE-ABC-A-607-Q 6ZGE-ABC-A-608-V 6ZGE-ABC-A-609-A 6ZGE-ABC-A-610-V 6ZGE-ABC-A-611-L 6ZGE-ABC-A-612-Y 6ZGE-ABC-A-613-Q 6ZGE-ABC-A-614-D 6ZGE-ABC-A-615-V 6ZGE-ABC-A-616-N 6ZGE-ABC-A-617-C 6ZGE-ABC-A-null-Z 6ZGE-ABC-A-null-Z 6ZGE-ABC-A-null-Z 6ZGE-ABC-A-null-Z 6ZGE-ABC-A-null-Z 6ZGE-ABC-A-null-Z 6ZGE-ABC-A-null-Z 6ZGE-ABC-A-null-Z 6ZGE-ABC-A-null-Z 6ZGE-ABC-A-null-Z 6ZGE-ABC-A-null-Z 6ZGE-ABC-A-null-Z 6ZGE-ABC-A-null-Z 6ZGE-ABC-A-null-Z 6ZGE-ABC-A-null-Z 6ZGE-ABC-A-633-W 6ZGE-ABC-A-634-R 6ZGE-ABC-A-635-V 6ZGE-ABC-A-636-Y 6ZGE-ABC-A-637-S 6ZGE-ABC-A-638-T 6ZGE-ABC-A-639-G 6ZGE-ABC-A-640-S 6ZGE-ABC-A-641-N 6ZGE-ABC-A-642-V 6ZGE-ABC-A-643-F 6ZGE-ABC-A-644-Q 6ZGE-ABC-A-645-T 6ZGE-ABC-A-646-R 6ZGE-ABC-A-647-A 6ZGE-ABC-A-648-G 6ZGE-ABC-A-649-C 6ZGE-ABC-A-650-L 6ZGE-ABC-A-651-I 6ZGE-ABC-A-652-G 6ZGE-ABC-A-653-A 6ZGE-ABC-A-654-E 6ZGE-ABC-A-655-H 6ZGE-ABC-A-656-V 6ZGE-ABC-A-657-N 6ZGE-ABC-A-658-N 6ZGE-ABC-A-659-S 6ZGE-ABC-A-660-Y 6ZGE-ABC-A-661-E 6ZGE-ABC-A-662-C 6ZGE-ABC-A-663-D 6ZGE-ABC-A-664-I 6ZGE-ABC-A-665-P 6ZGE-ABC-A-666-I 6ZGE-ABC-A-667-G 6ZGE-ABC-A-668-A 6ZGE-ABC-A-669-G 6ZGE-ABC-A-670-I 6ZGE-ABC-A-671-C 6ZGE-ABC-A-672-A 6ZGE-ABC-A-673-S 6ZGE-ABC-A-674-Y 6ZGE-ABC-A-675-Q 6ZGE-ABC-A-676-T 6ZGE-ABC-A-null-Z 6ZGE-ABC-A-null-Z 6ZGE-ABC-A-null-Z 6ZGE-ABC-A-null-Z 6ZGE-ABC-A-null-Z 6ZGE-ABC-A-null-Z 6ZGE-ABC-A-null-Z 6ZGE-ABC-A-null-Z 6ZGE-ABC-A-null-Z 6ZGE-ABC-A-null-Z 6ZGE-ABC-A-null-Z 6ZGE-ABC-A-null-Z 6ZGE-ABC-A-689-S 6ZGE-ABC-A-690-Q 6ZGE-ABC-A-691-S 6ZGE-ABC-A-692-I 6ZGE-ABC-A-693-I 6ZGE-ABC-A-694-A 6ZGE-ABC-A-695-Y 6ZGE-ABC-A-696-T 6ZGE-ABC-A-697-M 6ZGE-ABC-A-698-S 6ZGE-ABC-A-699-L 6ZGE-ABC-A-700-G 6ZGE-ABC-A-701-A 6ZGE-ABC-A-702-E 6ZGE-ABC-A-703-N 6ZGE-ABC-A-704-S 6ZGE-ABC-A-705-V 6ZGE-ABC-A-706-A 6ZGE-ABC-A-707-Y 6ZGE-ABC-A-708-S 6ZGE-ABC-A-709-N 6ZGE-ABC-A-710-N 6ZGE-ABC-A-711-S 6ZGE-ABC-A-712-I 6ZGE-ABC-A-713-A 6ZGE-ABC-A-714-I 6ZGE-ABC-A-715-P 6ZGE-ABC-A-716-T 6ZGE-ABC-A-717-N 6ZGE-ABC-A-718-F 6ZGE-ABC-A-719-T 6ZGE-ABC-A-720-I 6ZGE-ABC-A-721-S 6ZGE-ABC-A-722-V 6ZGE-ABC-A-723-T 6ZGE-ABC-A-724-T 6ZGE-ABC-A-725-E 6ZGE-ABC-A-726-I 6ZGE-ABC-A-727-L 6ZGE-ABC-A-728-P 6ZGE-ABC-A-729-V 6ZGE-ABC-A-730-S 6ZGE-ABC-A-731-M 6ZGE-ABC-A-732-T 6ZGE-ABC-A-733-K 6ZGE-ABC-A-734-T 6ZGE-ABC-A-735-S 6ZGE-ABC-A-736-V 6ZGE-ABC-A-737-D 6ZGE-ABC-A-738-C 6ZGE-ABC-A-739-T 6ZGE-ABC-A-740-M 6ZGE-ABC-A-741-Y 6ZGE-ABC-A-742-I 6ZGE-ABC-A-743-C 6ZGE-ABC-A-744-G 6ZGE-ABC-A-745-D 6ZGE-ABC-A-746-S 6ZGE-ABC-A-747-T 6ZGE-ABC-A-748-E 6ZGE-ABC-A-749-C 6ZGE-ABC-A-750-S 6ZGE-ABC-A-751-N 6ZGE-ABC-A-752-L 6ZGE-ABC-A-753-L 6ZGE-ABC-A-754-L 6ZGE-ABC-A-755-Q 6ZGE-ABC-A-756-Y 6ZGE-ABC-A-757-G 6ZGE-ABC-A-758-S 6ZGE-ABC-A-759-F 6ZGE-ABC-A-760-C 6ZGE-ABC-A-761-T 6ZGE-ABC-A-762-Q 6ZGE-ABC-A-763-L 6ZGE-ABC-A-764-N 6ZGE-ABC-A-765-R 6ZGE-ABC-A-766-A 6ZGE-ABC-A-767-L 6ZGE-ABC-A-768-T 6ZGE-ABC-A-769-G 6ZGE-ABC-A-770-I 6ZGE-ABC-A-771-A 6ZGE-ABC-A-772-V 6ZGE-ABC-A-773-E 6ZGE-ABC-A-774-Q 6ZGE-ABC-A-775-D 6ZGE-ABC-A-776-K 6ZGE-ABC-A-777-N 6ZGE-ABC-A-778-T 6ZGE-ABC-A-779-Q 6ZGE-ABC-A-780-E 6ZGE-ABC-A-781-V 6ZGE-ABC-A-782-F 6ZGE-ABC-A-783-A 6ZGE-ABC-A-784-Q 6ZGE-ABC-A-785-V 6ZGE-ABC-A-786-K 6ZGE-ABC-A-787-Q 6ZGE-ABC-A-788-I 6ZGE-ABC-A-789-Y 6ZGE-ABC-A-790-K 6ZGE-ABC-A-791-T 6ZGE-ABC-A-792-P 6ZGE-ABC-A-793-P 6ZGE-ABC-A-794-I 6ZGE-ABC-A-795-K 6ZGE-ABC-A-796-D 6ZGE-ABC-A-797-F 6ZGE-ABC-A-798-G 6ZGE-ABC-A-799-G 6ZGE-ABC-A-800-F 6ZGE-ABC-A-801-N 6ZGE-ABC-A-802-F 6ZGE-ABC-A-803-S 6ZGE-ABC-A-804-Q 6ZGE-ABC-A-805-I 6ZGE-ABC-A-806-L 6ZGE-ABC-A-807-P 6ZGE-ABC-A-808-D 6ZGE-ABC-A-809-P 6ZGE-ABC-A-810-S 6ZGE-ABC-A-811-K 6ZGE-ABC-A-812-P 6ZGE-ABC-A-813-S 6ZGE-ABC-A-814-K 6ZGE-ABC-A-815-R 6ZGE-ABC-A-816-S 6ZGE-ABC-A-817-F 6ZGE-ABC-A-818-I 6ZGE-ABC-A-819-E 6ZGE-ABC-A-820-D 6ZGE-ABC-A-821-L 6ZGE-ABC-A-822-L 6ZGE-ABC-A-823-F 6ZGE-ABC-A-824-N 6ZGE-ABC-A-825-K 6ZGE-ABC-A-826-V 6ZGE-ABC-A-827-T 6ZGE-ABC-A-828-L 6ZGE-ABC-A-829-A 6ZGE-ABC-A-830-D 6ZGE-ABC-A-831-A 6ZGE-ABC-A-832-G 6ZGE-ABC-A-833-F 6ZGE-ABC-A-834-I 6ZGE-ABC-A-835-K 6ZGE-ABC-A-836-Q 6ZGE-ABC-A-837-Y 6ZGE-ABC-A-838-G 6ZGE-ABC-A-839-D 6ZGE-ABC-A-840-C 6ZGE-ABC-A-841-L 6ZGE-ABC-A-842-G 6ZGE-ABC-A-843-D 6ZGE-ABC-A-844-I 6ZGE-ABC-A-845-A 6ZGE-ABC-A-846-A 6ZGE-ABC-A-847-R 6ZGE-ABC-A-848-D 6ZGE-ABC-A-849-L 6ZGE-ABC-A-850-I 6ZGE-ABC-A-851-C 6ZGE-ABC-A-852-A 6ZGE-ABC-A-853-Q 6ZGE-ABC-A-854-K 6ZGE-ABC-A-855-F 6ZGE-ABC-A-856-N 6ZGE-ABC-A-857-G 6ZGE-ABC-A-858-L 6ZGE-ABC-A-859-T 6ZGE-ABC-A-860-V 6ZGE-ABC-A-861-L 6ZGE-ABC-A-862-P 6ZGE-ABC-A-863-P 6ZGE-ABC-A-864-L 6ZGE-ABC-A-865-L 6ZGE-ABC-A-866-T 6ZGE-ABC-A-867-D 6ZGE-ABC-A-868-E 6ZGE-ABC-A-869-M 6ZGE-ABC-A-870-I 6ZGE- A-null-Z 6ZB4-A 6ZB4-B 6ZB4-C 6ZB5-A 6ZB5-B 6ZB5-C 6XR8-A 6XR8-B 6XR8-C 6ZGE-A 6ZGE-B 6ZGE-C 6ZGI-A 6ZGI-B 6ZGI-C 7JJI-A 7JJI-B 7JJI-C 7DF3-A 7DF3-B 7DF3-C 7DDD-A 7DDD-B 7DDD-C 6ZP2-A 6ZP2-B 6ZP2-C 6ZGH-C 6ZGH-B 7DDN-B 7DDN-C 7DF4-B 7DF4-D 7DF4-C 7DDN-A 6ZGG-B 7A95-C 6ZGG-A 7A95-A 6ZGG-C 7A93-A 7A93-C 7A95-B 7A97-B 7A98-A 7A98-B 7A98-C 7A97-C 7A97-A 7A96-A 7A96-B 7A96-C 7A93-B 7DK3-A 7DK3-B 7DCC-E 7DCC-I 7DCC-K 7DD2-C 7DD2-K 7DD8-C 7DCX-C 7DCX-K 7DK6-A 7DK6-B 7DK7-A 7DK7-B 7DK7-C 7DK5-A 7DD2-D 7DD8-D 7DD8-E 7DCX-D 7DK6-C 7DK5-B 7DK5-C 7DK4-A 7DK4-B 7DK4-C 7KSG-A 7KSG-B 6ZXN-A 7CWU-B 7CWS-O 7CWS-Q 7CWS-R 7CWU-A 7A25-A 7A29-A 7A29-B 7CWL-C 7A94-A 7CWL-A 7KSG-C 7CWM-C 7CWN-A 7CWN-C 7CWM-B 7A94-B 7CWN-B 7CWU-C 7CWL-B 7A94-C 6ZXN-B 6ZXN-C 7A25-B 7A25-C 7A29-C 7DK3-C 6ZP0-A 6ZP0-B 6ZP0-C 6ZOY-A 6ZOY-B 6ZOY-C 6ZP1-A 6ZP1-B 6ZP1-C 6ZOX-A 6ZOX-B 6ZOX-C 7K43-A 7K43-B 7K43-E 7CAI-A 7CAI-B 7CAI-C 7CAK-A 7CAK-B 7CAK-C 6XEY-C 6XEY-A 7CWM-A 7C2L-B 7C2L-C 6XEY-B 7KNI-A 7KNI-B 7KNI-C 7KNE-B 7KNH-A 7KNB-B 7KMZ-A 7KNH-B 7KMZ-B 7KMS-B 7KMS-A 7KMS-C 7CT5-A 7CT5-C 7CT5-B 7K8S-A 7K8S-B 7K8S-C 7K8Y-B 7K8Y-D 7K8Y-E 7K8X-A 7K8X-B 7K8X-C 7KKK-A 7KKK-C 7KKK-E 7KKL-A 7KKL-C 7KKL-D 6ZDH-A 6ZDH-B 6ZDH-C 6XF5-A 6XF5-B 6XF5-C 7K8T-A 7K8T-B 7K8T-C 7KL9-B 7K4N-A 7K4N-B 7K4N-E 7KNE-A 7KNE-C 7KNH-C 7KNB-A 7KNB-C 7KMZ-C 7KL9-A 7KL9-C 7JW0-A 7JW0-B 7JW0-E 7JVC-A 7JVC-B 7JVC-E 7K8U-C 6Z43-C 6XF6-C 6Z97-C 6ZHD-B 6Z43-B 6XF6-A 7K8W-G 6ZHD-C 7K8U-A 7JWB-A 7K8W-A 7JZN-C 6Z97-A 7CHH-A 6ZHD-A 6Z43-A 7CHH-C 7K8U-B 7JWB-B 6XCN-A 6XCN-C 6XCN-E 6XCM-B 6XCM-C 7JZL-C 7K90-A 7K90-B 7K90-C 6Z97-B 7BYR-B 7KJ3-A 7KJ3-B 7KJ2-B 7KJ4-A 7KJ4-B 7KJ4-C 7KJ5-B 7JZN-B 6XM0-C 6XM0-A 6XM4-A 6XM5-C 6XM3-A 7CHH-B 7JZN-A 7JZL-A 7JZL-B 7K8V-C 6ZWV-A 6ZWV-B 6ZWV-C 7K8V-A 7K8V-B 7LD1-A 7LD1-B 7JV6-A 7JV6-B 7JV6-E 7K8Z-C 7LCN-A 7LCN-C 7LCN-K 6X6P-A 6X6P-B 6X6P-C 6XM3-C 7KDI-A 7KDI-B 7KDI-C 7KDG-A 7KDG-B 7KDG-C 7L06-A 7L06-B 7L06-C 7KE6-A 7KE6-B 7KE6-C 6X2C-A 6X2C-B 6X2C-C 7KDK-A 7KDK-B 7KDK-C 7KE4-A 7KE4-B 7KE4-C 7L02-A 7L02-B 7L02-C 7L09-A 7L09-B 7L09-C 7KE8-A 7KE8-B 7KE8-C 6X29-A 6X29-B 6X29-C 7KE7-A 7KE7-B 7KE7-C 6VXX-A 6VXX-B 6VXX-C 7KDH-A 7K8Z-A 7KJ2-A 7KJ5-A 7KDJ-A 7KEB-A 7KE9-A 6VYB-A 6X2A-A 7KEC-A 7KEA-A 6XCM-A 7KDL-A 7KDH-C 6X2B-A 6X2A-C 7KJ3-C 7KJ2-C 7KJ5-C 7KDJ-C 7KEB-C 7KE9-C 6VYB-C 7KEC-C 7KEA-C 7KDL-C 6WPT-C 6WPS-A 6WPS-B 6WPS-E 7KDH-B 6X79-A 6X79-B 6X79-C 6XKL-C 6WPT-A 6XM4-C 7JWY-A 6XLU-A 7BYR-C 6VSB-C 6XKL-B 7K8W-B 7LD1-C 6X2B-B 6X2B-C 7KDJ-B 7KEB-B 7KE9-B 7K8Z-B 6VYB-B 7KEC-B 7KEA-B 7KDL-B 7JWB-C 6X2A-B 6VSB-B 7A4N-A 7A4N-B 7A4N-C 6WPT-B 7AD1-C 7AD1-A 6XF6-B 7JWY-C 6XLU-C 6XM5-A 7CN9-C 6XM3-B 6XM4-B 6ZOZ-A 6ZOZ-B 6ZOZ-C 7JWY-B 6XLU-B 6XM5-B 6XM0-B 6XKL-A 7AD1-B 7C2L-A 7BYR-A 6VSB-A 7CN9-B 7CN9-A 6XS6-C 6XS6-A 6XS6-B 6ZGH Structural basis for the recognition of SARS-CoV-2 by full-length human ACE2 Status Report on COVID-19 Vaccines Development Structures and function of locked conformations of SARS-CoV-2 spike Free fatty acid binding pocket in the locked structure of SARS-CoV-2 spike protein Distinct conformational states of SARS-CoV-2 spike protein Syncytia formation by SARS-CoV-2-infected cells Coronavirus biology and replication: implications for SARS-CoV-2 Coronavirus entry: how we arrived at SARS-CoV-2 Receptor binding and priming of the spike protein of SARS-CoV-2 for membrane fusion Cryo-EM Structures of SARS-CoV-2 Spike without and with ACE2 Reveal a pH-Dependent Switch to Mediate Endosomal Positioning of Receptor-Binding Domains Conformational flexibility and structural variability of SARS-CoV2 S protein beta-Coronaviruses Use Lysosomes for Egress Instead of the Biosynthetic Secretory Pathway Thermodynamics and SARS-CoV-2: neurological effects in post-Covid 19 syndrome Structural and Functional Analysis of the D614G SARS-CoV-2 Spike Protein Variant Structural impact on SARS-CoV-2 spike protein by D614G substitution Enhancing the Prefusion Conformational Stability of SARS-CoV-2 Spike Protein Through Structure-Guided Design Structural Basis of SARS-CoV-2 and SARS-CoV Antibody Interactions Biomolecular Simulations in the Time of COVID19, and After AI-Driven Multiscale Simulations Illuminate Mechanisms of SARS-CoV-2 Spike Dynamics Molecular Simulations suggest Vitamins, Retinoids and Steroids as Ligands of the Free Fatty Acid Pocket of the SARS-CoV-2 Spike Protein* Uncovering cryptic pockets in the SARS-CoV-2 spike glycoprotein Temperature effect on the SARS-CoV-2: A molecular dynamics study of the spike homotrimeric glycoprotein Site-specific steric control of SARS-CoV-2 spike glycosylation Assessment of proton-coupled conformational dynamics of SARS and MERS coronavirus papain-like proteases: Implication for designing broad-spectrum antiviral inhibitors Proton-Coupled Conformational Activation of SARS Coronavirus Main Proteases and Opportunity for Designing Small-Molecule Broad-Spectrum Targeted Covalent Inhibitors Improved pKa calculations through flexibility based sampling of a waterdominated interaction scheme protein-sol pKa: prediction of electrostatic frustration, with application to coronaviruses The worldwide Protein Data Bank (wwPDB): ensuring a single, uniform archive of PDB data UniProt: a worldwide hub of protein knowledge Basic local alignment search tool A thermostable, closed SARS-CoV-2 spike protein trimer Making the invisible enemy visible Fast, scalable generation of high-quality protein multiple sequence alignments using Clustal Omega The EMBL-EBI search and sequence analysis tools APIs in 2019 Web-based display of protein surface and pH-dependent properties for assessing the developability of biotherapeutics Prediction of pH-dependent properties of proteins Multiple-site titration and molecular modeling: two rapid methods for computing energies and forces for ionizable groups in proteins Improved Treatment of Ligands and Coupling Effects in Empirical Calculation and Rationalization of pKa Values Improvements to the APBS biomolecular solvation software suite NGL Viewer: Web-based molecular graphics for large complexes SWISS-MODEL and the Swiss-PdbViewer: an environment for comparative protein modeling A model for pH coupling of the SARS-CoV-2 spike protein open/closed equilibrium Cryo-EM model validation recommendations based on outcomes of the 2019 EMDataResource challenge SARS-CoV-2 D614G spike mutation increases entry efficiency with enhanced ACE2-binding affinity Tracking Changes in SARS-CoV-2 Spike: Evidence that D614G Increases Infectivity of the COVID-19 Virus Structural impact on SARS-CoV-2 spike protein by D614G substitution Measurement of the individual pKa values of acidic residues of hen and turkey lysozymes by two-dimensional 1H NMR Influence of a heptad repeat stutter on the pH-dependent conformational behavior of the central coiled-coil from influenza hemagglutinin HA2 pH stability of HLA-DR4 complexes with antigenic peptides Residues PsaB Asp612 and PsaB Glu613 of photosystem I confer pH-dependent binding of plastocyanin and cytochrome c(6) Chaperone activation by unfolding Investigation of anomalous charge variant profile reveals discrete pH-dependent conformations and conformationdependent charge states within the CDR3 loop of a therapeutic mAb High apparent dielectric constants in the interior of a protein reflect water penetration Viral cell recognition and entry The Role of histidine residues in low-pH-mediated viral membrane fusion SARS-CoV-2 convergent evolution as a guide to explore adaptive advantage The furin cleavage site in the SARS-CoV-2 spike protein is required for transmission in ferrets The SARS-CoV-2 variants associated with infections in India, B.1.617, show enhanced spike cleavage by furin Stabilizing the closed SARS-CoV-2 spike trimer A-272-P 6ZGE-ABC-A-273-R 6ZGE-ABC-A-274-T 6ZGE-ABC-A-275-F 6ZGE-ABC-A-276-L 6ZGE-ABC-A-277-L 6ZGE-ABC-A-278-K 6ZGE-ABC-A-279-Y 6ZGE-ABC-A-280-N 6ZGE-ABC-A-281-E 6ZGE-ABC-A-282-N 6ZGE-ABC-A-283-G 6ZGE-ABC-A-284-T 6ZGE-ABC-A-285-I 6ZGE-ABC-A-286-T 6ZGE-ABC-A-287-D 6ZGE-ABC-A-288-A 6ZGE-ABC-A-289-V 6ZGE-ABC-A-290-D 6ZGE-ABC-A-291-C 6ZGE-ABC-A-292-A 6ZGE-ABC-A-293-L 6ZGE-ABC-A-294-D 6ZGE-ABC-A-295-P 6ZGE-ABC-A-296-L 6ZGE-ABC-A-297-S 6ZGE-ABC-A-298-E 6ZGE-ABC-A-299-T 6ZGE-ABC-A-300-K 6ZGE-ABC-A-301-C 6ZGE-ABC-A-302-T 6ZGE-ABC-A-303-L 6ZGE-ABC-A-304-K 6ZGE-ABC-A-305-S 6ZGE-ABC-A-306-F 6ZGE-ABC-A-307-T 6ZGE-ABC-A-308-V 6ZGE-ABC-A-309-E 6ZGE-ABC-A-310-K 6ZGE-ABC-A-311-G 6ZGE-ABC-A-312-I 6ZGE-ABC-A-313-Y 6ZGE-ABC-A-314-Q 6ZGE-ABC-A-315-T 6ZGE-ABC-A-316-S 6ZGE-ABC-A-317-N 6ZGE-ABC-A-318-F 6ZGE-ABC-A-319-R 6ZGE-ABC-A-320-V 6ZGE-ABC-A-321-Q 6ZGE-ABC-A-322-P 6ZGE-ABC-A-323-T 6ZGE-ABC-A-324-E 6ZGE-ABC-A-325-S 6ZGE-ABC-A-326-I 6ZGE-ABC-A-327-V 6ZGE-ABC-A-328-R 6ZGE-ABC-A-329-F 6ZGE-ABC-A-330-P 6ZGE-ABC-A-331-N 6ZGE-ABC-A-332-I 6ZGE-ABC-A-333-T 6ZGE-ABC-A-334-N 6ZGE-ABC-A-335-L 6ZGE-ABC-A-336-C 6ZGE-ABC-A-337-P 6ZGE-ABC-A-338-F 6ZGE-ABC-A-339-G 6ZGE-ABC-A-340-E 6ZGE-ABC-A-341-V 6ZGE-ABC-A-342-F 6ZGE-ABC-A-343-N 6ZGE-ABC-A-344-A 6ZGE-ABC-A-345-T 6ZGE-ABC-A-346-R 6ZGE-ABC-A-347-F 6ZGE-ABC-A-348-A 6ZGE-ABC-A-349-S 6ZGE-ABC-A-350-V 6ZGE-ABC-A-351-Y 6ZGE-ABC-A-352-A 6ZGE-ABC-A-353-W 6ZGE-ABC-A-354-N 6ZGE-ABC-A-355-R 6ZGE-ABC-A-356-K 6ZGE-ABC-A-357-R 6ZGE-ABC-A-358-I 6ZGE-ABC-A-359-S 6ZGE-ABC-A-360-N 6ZGE-ABC-A-361-C 6ZGE-ABC-A-362-V 6ZGE-ABC-A-363-A 6ZGE-ABC-A-364-D 6ZGE-ABC-A-365-Y 6ZGE-ABC-A-366-S 6ZGE-ABC-A-367-V 6ZGE-ABC-A-368-L 6ZGE-ABC-A-369-Y 6ZGE-ABC-A-370-N 6ZGE-ABC-A-371-S 6ZGE-ABC-A-372-A 6ZGE-ABC-A-373-S 6ZGE-ABC-A-374-F 6ZGE-ABC-A-375-S 6ZGE-ABC-A-376-T 6ZGE-ABC-A-377-F 6ZGE-ABC-A-378-K 6ZGE-ABC-A-379-C 6ZGE-ABC-A-380-Y 6ZGE-ABC-A-381-G 6ZGE-ABC-A-382-V 6ZGE-ABC-A-383-S 6ZGE-ABC-A-384-P 6ZGE-ABC-A-385-T 6ZGE-ABC-A-386-K 6ZGE-ABC-A-387-L 6ZGE-ABC-A-388-N 6ZGE-ABC-A-389-D 6ZGE-ABC-A-390-L 6ZGE-ABC-A-391-C 6ZGE-ABC-A-392-F 6ZGE-ABC-A-393-T 6ZGE-ABC-A-394-N 6ZGE-ABC-A-395-V 6ZGE-ABC-A-396-Y 6ZGE-ABC-A-397-A 6ZGE-ABC-A-398-D 6ZGE-ABC-A-399-S 6ZGE-ABC-A-400-F 6ZGE-ABC-A-401-V 6ZGE-ABC-A-402-I 6ZGE-ABC-A-403-R 6ZGE-ABC-A-404-G 6ZGE-ABC-A-405-D 6ZGE-ABC-A-406-E 6ZGE-ABC-A-407-V 6ZGE-ABC-A-408-R 6ZGE-ABC-A-409-Q 6ZGE-ABC-A-410-I 6ZGE-ABC-A-411-A 6ZGE-ABC-A-412-P 6ZGE-ABC-A-413 Monomer-monomer contact area (d-SASA) in D614Gset spike trimers S2-d-SASA S1-d-SASA The authors thank Lorena Zuzic and Dr Agnel Joseph for discussions, and the computational shared facility at the University of Manchester for support and resources.