key: cord-0929328-vjkied8h
authors: Jana, Asis K.; Greenwood, Augustus B.; Hansmann, Ulrich H. E.
title: Presence of a SARS-CoV-2 Protein Enhances Amyloid Formation of Serum Amyloid A
date: 2021-08-09
journal: J Phys Chem B
DOI: 10.1021/acs.jpcb.1c04871
sha: a7c1233cf08e7db310978af3d32bf08e952b5a5f
doc_id: 929328
cord_uid: vjkied8h

[Image: see text] A marker for the severeness and disease progress of COVID-19 is overexpression of serum amyloid A (SAA) to levels that in other diseases are associated with a risk for SAA amyloidosis. To understand whether SAA amyloidosis could also be a long-term risk of SARS-CoV-2 infections, we have used long all-atom molecular dynamic simulations to study the effect of a SARS-CoV-2 protein segment on SAA amyloid formation. Sampling over 40 μs, we find that the presence of the nine-residue segment SK9, located at the C-terminus of the envelope protein, increases the propensity for SAA fibril formation by three mechanisms: it reduces the stability of the lipid-transporting hexamer shifting the equilibrium toward monomers, it increases the frequency of aggregation-prone configurations in the resulting chains, and it raises the stability of SAA fibrils. Our results therefore suggest that SAA amyloidosis and related pathologies may be a long-term risk of SARS-CoV-2 infections.

While even after serious complications, most COVID-19 survivors appear to recover completely, only little is known about the long-term effects of infections by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Diseaseassociated symptoms such as inflammation of blood vessels and overreaction of the immune system 1−3 are connected with spikes in the concentration of the human serum amyloid A (SAA) protein, with the level increasing as the disease progresses from mild to critical. 4−6 The more than thousand times higher blood concentrations of SAA in acute COVID-19 patients 4,5 are comparable to the ones seen in patients with various cancers or inflammatory diseases 7 where the overexpression of SAA is associated with systemic amyloidosis as a secondary illness. SAA amyloidosis is characterized by the formation and deposition of SAA amyloids in the blood vessels, causing inflammation, thrombosis, and eventually organ damage. Common complications of SAA amyloidosis such as kidney failure or high incidents of thrombosis are also frequently observed in COVID-19 patients. 8, 9 The similarity of symptoms suggests that SAA amyloidosis may exacerbate COVID-19 symptoms, 10 or that it is a long-term risk in COVID-19 survivors causing, for instance, the broad spectrum of symptoms in the multisystem inflammatory syndrome first reported in children and adolescents (MIS-C), 11 but also observed in adults (MIS-A). This hypothesis is the motivation for the present study where we use molecular dynamics simulations to probe how the presence of a SARS-CoV-2 protein fragment modulates the formation and stability of SAA amyloids. Such SARS-CoV-2-triggered amyloid formation has been observed in vitro for αSynuclein, 12 but not yet demonstrated for SAA.

The overexpression of SAA in some cancers or inflammatory diseases leads not in all cases to amyloidosis. Usually, after a spike, concentration levels decrease rapidly in a process that involves dissociation of SAA hexamers followed by cleavage of the released chains. In our previous work, 13 we proposed that the cleavage happens in part because fragments have a lower probability to reassemble into functional hexamers than the complete SAA 1−104 proteins. We also observed that, unlike other fragments, the most commonly found SAA 1−76 can switch between two structural motifs. The first one is easy to proteolyze (allowing to lower rapidly the SAA concentration) but vulnerable for aggregation, while the opposite is the case for the second motif. If amyloid formation takes longer than proteolysis, the aggregation-prone species dominates. However, if environmental conditions such as low pH encourage amyloid formation, the configurational ensemble shifts toward the more protected form. In this picture, amyloidosis happens when this mechanism for downregulating SAA concentration becomes overwhelmed or otherwise fails. In COVID-19 patients, this could happen in three ways: first, the presence of SARS-CoV-2 could lower the stability of the functional hexamers; second, it could increase the probability for the association of SAA fragments after cleavage; and third, it could enhance the stability of the resulting SAA fibrils; each possibility is shifting the equilibrium toward SAA fibril formation.

In this computational study, we probe all three possibilities but focus on the effect of viral proteins, i.e., neglecting the potential roles played by viral RNA. To reduce computational costs, we restrict ourselves to short amyloidogenic regions on viral proteins that are most likely to interact with SAA. An example is the nine-residue-segment S 55 FYVYSRVK 63 (SK9) on the C-terminal tail of the SARS-COV2-Envelope protein, see Figure 1 . While most of the 75-residue Envelope proteins are transmembrane or intracellular, this segment is located on the extracellular C-terminal tail. As its location makes it likely to interact with extracellular SAA proteins, and as its homolog in SARS-CoV-1 is known to form amyloids in solution, 15 we investigate in our simulations the interaction of the SK9 segment with the SAA hexamers, monomeric SAA 1−76 -fragments, and the SAA fibrils as shown in Figure 2 . While the use of such a small segment may lead to a different mechanism than one would see for the full protein, the danger seems minimal in our cases, as most of the Envelope protein can likely not interact with the extracellular SAA.

Our results show that the presence of viral SK9 fragment raises the risk for SAA fibril formation at all three stages: it reduces the stability of the lipid-transporting hexamer, by increasing the chance for the association of the SAA fragments after cleavage, and by enhancing the stability of the resulting SAA fibrils, in each case, shifting the equilibrium toward SAA fibril formation. Our results therefore suggest that SAA amyloidosis and related pathologies may be a long-term risk of SARS-CoV-2 infections.

System Preparation. To evaluate the effect of the nineresidue SK9 segment S 55 FYVYSRVK 63 (located in the Cterminal tail of the Envelope protein of the SARS-CoV-2 virus) on serum amyloid A (SAA) amyloid formation, we monitor the change in stability of the known hexamer, monomer, and fibril models upon binding, comparing the complex of SK9 and SAA to the corresponding "pure" SAA models.

We choose as the initial configuration for the SAA 1−104 hexamer the X-ray-resolved crystal structure, deposited in the Protein Data Bank (PDB) under identifier 4IP8 16 and shown by us in Figure 2a ,b. Note that we add here and in the following cases a NH 3 + -group at the N-terminus and a COO − -end group at the C-terminus. These end groups are chosen to make our simulations consistent with our previous work. 13 As in vivo SAA monomers are cleaved enzymatically, with SAA 1−76 the most common resulting fragment, we consider the following three models for our simulation of monomers. The first one, shown in Figure 2c , is generated by removing the residues 77−104 from the crystal structure of the full-length SAA monomer (PDB ID: 4IP9), 16 while the other two configurations were derived by us in ref 13 as typical motifs found in long-time simulations of the fragment, and named by us "helix-weakened" (Figure 2d ) and "helix-broken" (Figure 2e ). Finally, for simulation of SAA fibrils, we have used tetramers made of two folds (protofibrils) and two layers since we have identified in previous work 17 such 2F2L tetramers as the smallest stable fibril fragments. Our 2F2L model is shown in Figure 2f and is derived from the cryo-EM structure deposited in the PDB database under identifier 6MST. 18 Note that the fibril model is made of SAA fragments with residues 2− 55 since no fibril model for human SAA 1−76 is available. Initially, we also built a second fibril model, where we added the likely disordered missing residues 56−76 in a configuration predicted by homology modeling. However, we discarded this model for reasons discussed in the Results section.

Simulations starting from the above-described SAA models serve as a control against which we compare our simulations of the various SAA models interacting with SK9 segments S 55 FYVYSRVK 63 , which is located in the C-terminal tail of the Envelope protein of the SARS-CoV-2 virus. As with the exception of the transmembrane domain of residues 8−38 19 the Envelope protein has not been resolved, we have generated the initial configuration for the SK9 fragment from a model derived by a machine-learning approach and subsequent refinement by molecular dynamics. 20 For this purpose, we have removed from this model the residues 1−54 and 64−75, afterward capping the remaining nine-residue segment with a NH 3 + -group at the N-terminus and CONH 2 as the C-terminus. These end groups were chosen to avoid strong electrostatic interactions between the oppositely charged terminal residues.

Using the AutoDock Vina software, 21 we have generated start configurations for our simulations by docking SK9 segments in a ratio 1:1 with SAA chains in our models. The resulting binding positions of SK9 segments in the start configurations of either SAA hexamer, monomers, or fibrils are also shown in Figure 2 . In the case of the hexamer, binding positions are obtained from a global search, and the configuration with the lowest energy is chosen as the start configuration. On the other hand, we have performed three independent searches for each SAA 1−76 monomer, focusing on either the N-terminus, the helix-I− helix-II linker, or the C-terminal region and chose the respective lowest energy binding pose as the start configuration of the , and the C-terminal helix-IV (residues 73−88). N-and C-terminal residues are here and in all other subfigures represented by blue and red spheres, respectively. In the start configuration, the virus protein segment SK9 binds with SAA 1−104 hexamer at the N-terminal helix-I (green) or residues in the helix-III (yellow). The monomer structure of the SAA 1−76 fragment in (c) is derived from the X-ray crystal structure (PDB ID: 4IP9), while the representative configurations for the helix weakened (d) and helix broken (e) SAA 1−76 monomers are taken from our previous work. 13 The initial binding positions of SK9 at Nterminal helix-I (green), helix-I-helix-II linker (orange), and disordered C-terminus (yellow) are also shown. Finally, we show in (f) the fibril fragment 2F2L as extracted from the cryo-EM structure of the human SAA fibril (PDB ID: 6MST). In the start configuration, the SK9 segment binds to the SAA 2−55 fibril at either the N-terminus (green), C-terminal cavity (violet), or the packing interface (orange). The Journal of Physical Chemistry B pubs.acs.org/JPCB Article Simulation Protocol. All simulations are carried out using the GROMACS 2018 package 23 and are employing the CHARMM 36m all-atom force field 24 and TIP3P water. 22 Number and length of runs are also listed in Table 1 . For each system, we have taken the configurations generated above and minimized their energy using the steepest-descent algorithm. This step is followed by 200 ps of molecular dynamic simulations in the NVT ensemble (keeping the volume constant) at 310 K, and a subsequent run of 200 ps in the NPT ensemble keeping the pressure at 1 atm. During the NVT and NPT equilibrations, the nonhydrogen (heavy) atoms of protein are restrained with a force constant of 1000 kJ mol −1 nm −2 .

The so-equilibrated configurations are the start point of our production runs where the respective systems evolve at a constant temperature of 310 K and a constant pressure of 1 atm. For each setup, we followed three independent trajectories starting from different initial velocity distributions. Note that in these simulations the position of the SK9 segment is not restrained to the original docking sites. Instead, it can move and even detach over the course of the simulation, no longer interacting with SAA. If this happens, the behavior of the system will be the same as for the control (where SK9 is absent). Hence, to save computational resources, we stopped a simulation if the SK9 segment did not reattach within 100 ns. The length of all trajectories is listed in Table 1 .

The temperature is controlled during the simulation with a vrescale thermostat, 25 and the pressure with the Parrinello− Rahman barostat, 26 using a coupling constant of 2 ps. Keeping the water geometry fixed with SETTLE algorithm 27 and constraining non-water bonds including hydrogen atoms with the LINCS algorithm 28 allowed us to use a timestep of 2 fs for integrating the equations of motion. As we use threedimensional orthorhombic periodic boundary conditions, we have to use the particle-mesh Ewald (PME) 29 method for calculating electrostatic interactions. This is done with a realspace cutoff of 12 Å; a value also used as cutoff for Van der Waal interactions, where smoothing started at 10.5 Å.

Trajectory Analysis. Most of our analysis relies on GROMACS tools such as gmx_rms for calculating the rootmean-square deviation (RMSD) with respect to the initial configuration. Another example is the do_dssp tool that implements the dictionary of secondary structure in proteins (DSSP) 30 and allows calculation of residue-wise secondary structure propensities. For visualization and for calculating the solvent-accessible surface area (SASA), we use VMD software. 31 Quantities such as the cavity diameter ⟨d cavity ⟩ in the hexamer are calculated by in-house programs, averaging over the center-ofmass distances between the N-terminal helix-I regions of adjacent units of both layers. Another example is the fraction of native contacts, which is calculated using a soft cutoff algorithm and which is defined as 32 Here, r is the distance between two heavy atoms at which a contact is formed in the native state and r 0 denotes this distance in the native state. β denotes the smoothing parameter taken to be 5 Å, while λ represents the fluctuation when contact is formed, taken to be 1.8. Residue-wise binding affinities of the SK9 segment to SAA chains are estimated by calculating binding probabilities instead of free energies. This is because exact methods such as thermodynamic integration would have been too costly for calculating the later and approximate approaches such as MM/ PBSA or MM/GBSA perform poorly when, as in our case, electrostatic interactions between charged (R61 and K63 in the SK9 segment) and polar residues dominate. 33 Here, we define a binding site as the closest residue that has at least one nonhydrogen atom within 4.5 Å from the SK9 segment.

Effect of SK9 Segments on SAA Hexamers. One of the symptoms of COVID-19 is an increase in the concentration of SAA to levels that in a number of cancers and inflammatory diseases often, but not always, leads to amyloidosis as a secondary illness. The response to the overexpression of SAA includes in these diseases a dissociation of hexamers. The subsequent cleavage of the resulting SAA monomers into fragments is likely because these fragments have a lower probability than the full SAA 1−104 protein to reassemble as hexamer. 13 Hence, we start our investigation into the effect of the SK9 segment on SAA amyloidosis by probing how the presence of SK9 alters the equilibrium between hexamers and monomers.

As the monomer−hexamer equilibrium depends on the stability of the SAA hexamer, we explore first the change in stability of the SAA hexamer upon binding of SK9 segments. The loss of stability can be monitored by measuring over the length of the trajectory the root-mean-square deviation (RMSD) to the start configuration. If this RMSD is calculated over all heavy atoms in the hexamer, we call this a global RMSD. On the other hand, if the RMSD is calculated separately for each chain, and averaged over all six chains, we talk of a chain RMSD. Hence, the global RMSD measures the structural deviation of the entire hexamer, whereas the chain RMSD represents the structural distortion of each chain in the hexamer. Note that we evaluate in both cases the RMSD only for residues 1−76 to stay consistent with our previous work. 13 In Figure 3a ,b, we plot the time evolution of the global RMSD, averaged over three independent trajectories, comparing the case of SAA bound with the SK9 segment to our control, the hexamer in the absence of SK9. The global RMSD for the SAA hexamer rises in the control simulation initially by nearly 3.0 Å, but quickly reaches a plateau within the first 100 ns, and over the next 400 ns increases only gradually. On the other hand, when SK9 segments are bound to SAA chains, the global RMSD rises in the first 100 ns also by about 3.0 Å, but instead of approaching a plateau, it continues to rise in a stepwise fashion to a RMSD of 5.0 Å and higher. The differences are much smaller for the chain RMSD, indicating that the interaction with SK9 segments disturbs the association of SAA chains but not the configuration of individual chains.

This can also be seen by comparing the final configurations of both systems in Figure 3c ,d, where we show the binding of the SK9 segments to the interfacial regions and the resulting distortion of the hexamer. To quantify this distortion, we plot in Figure 3e ,f the time evolution of the cavity diameter ⟨d cavity ⟩ for both systems. Strongly correlated with the global RMSD, we find for the hexamer in the absence of SK9 that the cavity diameter approaches quickly a plateau, while for the hexamer in complex with SK9 segments, the values again increase stepwise, with a pronounced step in ⟨d cavity ⟩ at 250 ns.

To understand in more detail how SK9 affects the association of the SAA chains in the hexamer, we first calculate the residuewise binding probability of the SK9 segment toward SAA hexamer. Data are averaged over the final 500 ns of each trajectory and shown in Figure 4 . The binding probability map demonstrates that the SK9 segment preferentially binds with interfacial aromatic and hydrophobic residues (F3, F4, W18, and I65) and interchain salt-bridge forming residues (D12 and R25), thereby disrupting interfacial hydrophobic and stacking contacts, as well as interchain salt bridges. As a consequence, the solvent exposure of the hydrophobic residues in the hexamer increases. This can be seen in Figure 5a where we show the solvent-accessible surface area (SASA) of exposed hydrophobic residues, measured in each of the three trajectories over the last 500 ns using VMD with a spherical probe of 1.4 Å radius. We find that in the presence of SK9 the SASA distribution is shifted toward higher values indicating solvent exposure of interfacial residues. The mean SASA value for hydrophobic residues in the presence of SK9 is (10 320 ± 418) Å 2 , and (10 426 ± 260) Å 2 in absence of SK9. Thus, the hydrophobic solvent-accessible surface increases by about 9% in the presence of SK9 segments, leading to the lower stability of the hexamer indicated by higher global RMSD values in Figure 3a .

The loss in contacts resulting from binding to the SK9 segment can be seen in Figure 5b , where we plot the time evolution of the number of native interchain contacts, i.e., contacts between SAA chains in the hexamer that also exists in the start configuration. We define contacts by a distance cutoff of 4.5 Å and calculate the fraction of interchain native contacts as described in the Methods section. In both systems, this fraction decreases within the first 50 ns. However, it reaches a stable value with little fluctuation after 100 ns in the control simulation, whereas in the presence of SK9, the loss of contacts continues and is clearly more pronounced after about 250 ns. For a more fine-grained picture, we have calculated the inter-residue contact probabilities between monomers, defining two residues whose heavy atoms lie within 7 Å of each other as a contact pair. Analyzing these contacts, we find multiple π−π stacking and hydrophobic interactions involving the residues F 3 , F 4 , W 18 , F 68 , However, when present, SK9 segments bind with the abovelisted interfacial residues, disrupting the interchain network in the hexamer, and the number of interfacial contacts decreases to an average of (704 ± 81) contacts. As an example, we show in Figure 5c ,d representative snapshots of the A and C chains from the control simulation and from the hexamer binding with SK9.

In the control simulation (Figure 5c ), residues F 3 and F 4 from chain A are in contact with residues F 68 and F 69 from chain C, while in Figure 5d , tyrosine and valine residues from a SK9 segment bind with the residues F 3 and F 4 from chain A, thereby preventing π−π stacking with residues F 68 and F 69 from chain C. We have quantified this effect for the salt-bridge forming residues D 12 and R 25 . For this pair, we plot in Figure 5e the probability distribution of the center-of-mass distance for both the control, and for the hexamer in the presence of SK9 segments. Measurements for both systems are taken over the last 500 ns of all three respective trajectories. The shift in the centerof-mass distance toward larger values in the presence of SK9 demonstrates the disruption of the interfacial salt bridge formed by the two residues D 12 and R 25 in the control. The loss of this salt bridge in the presence of SK9 destabilizes the SAA hexamer. This can be seen in Figure 5f where for a representative trajectory of the SAA hexamer in the presence of SK9 segments we compare the time evolution of the center-of-mass distance between D 12 of chain A and R 25 of chain C with the time evolution of the global RMSD. The strong correlation between both quantities is especially visible at 250 ns, where the loss of the salt bridge (once the distance between the two residues is larger than 4.0 Å) is mirrored by a jump in RMSD. Hence, our first result is that the binding of the SK9 segment with the SAA hexamer competes with the interchain binding of SAA proteins, reducing the stability of the hexamer. Our computational resources did not allow us to observe disassembly of the hexamer, but the loss of stability indicates that in the presence of SK9 segments the equilibrium is shifted away from the hexamer and toward monomers. While this shift downregulates the activity of SAA, it also increases the chance for aggregation.

Effect of SK9 Segments on SAA 1−76 Monomers. We argue in ref 13 that after dissociation of the hexamer monomers are cleaved into fragments because this eases their proteolysis. The most common fragment SAA 1−76 is special in that it evolves into an ensemble of configurations dominated by two forms. The first one is easy to proteolyze (allowing for a quick decrease in SAA concentration) but vulnerable to aggregation, while the situation is reversed for the second motif. If amyloid formation takes longer than proteolysis, the aggregation-prone species prevails. However, if external factors encourage amyloid formation, the frequency of the more protected motif increases. Hence, in the second part of our investigation, we study how the interaction with SK9 alters this mechanism.

For this purpose, we have simulated three systems. The first one is a SAA 1−76 -fragment in the same configuration as seen in the full-sized free monomer (PDB ID: 4IP9). This is presumably The Journal of Physical Chemistry B pubs.acs.org/JPCB Article the structure of the fragment right after cleavage. We follow the time evolution of this fragment in the presence of SK9 segments in three molecular dynamics simulations, with the SK9 segment initially docked to either the helix-I−helix-II linker, the Nterminus, or to the C-terminal region. We then compare our results from these three simulations with control simulations, where the SK9 segment is absent. When initially binding to the helix-I−helix-II linker, the SK9 segment separates within 500 ns from the SAA monomer and moves away, easily monitored by the time evolution of the number of contacts between SK9 and the SAA monomer, a quantity that at separation drops from fluctuating between 100 and 300 to zero. As after separation, the time evolution of the SAA chain is similar to the control simulation; we discuss in the following only the other two cases. For instance, we show in Figure 6 the time evolution of the secondary structure for the trajectory where the SK9 segment initially binds with the disordered C-terminal region. Over the course of the simulation, the SK9 segment shifts its binding partners to hydrophobic residues in N-terminal helix-I (L 7 , A 10 , and F 11 ) and some residues encompassing the helix-III and disordered C-terminal region (W 53 , A 54 , E 56 , A 57 , D 60 , N 64 , I 65 , R 67 , F 68 , G 70 , and E 74 ), see also the corresponding snapshot in Figure 6 . As a consequence of this binding pattern, we observe formation of a β-hairpin, involving residues 28−30 and 33−35 in the helix-I−helix-II linker and in the beginning of helix-III. Helix-II (residues 32−47) unfolds completely over the course of the simulation. This β-hairpin formation was also observed in a previous study by Gursky et al., 34 but is on the time scale of our simulations not seen in the control simulations. Both the unfolding of helix-I and helix-II and the β-hairpin formation in the helix-I−helix-II linker region are also not observed, when the SK9 segment binds initially to the N-terminal helix-I, but the two trajectories are otherwise similar. Configurations in both trajectories resemble the easy to proteolyze, but aggregationprone, helix-weakened configurations of Figure 2d , sharing a similar reduced helicity of residues 63−69 in helix-III (the mean helicity percentages are (43 ± 41) and (69 ± 31)%, respectively) and a comparable solvent-accessible surface area for the first 11 N-terminal residues (with mean SASA values of (524 ± 90) and (552 ± 43) Å 2 respectively).

Without the viral amyloidogenic segment present, the SAA 1−76 -fragment can also evolve into the helix-broken form of Figure 2e , which is more difficult to proteolyze but also less aggregation-prone. Our simulations indicate that the presence of SK9 alters the distribution of the two forms, shifting it to the more aggregation-prone helix-weakened form. Depending on the effect that the presence of SK9 has on the stability of the two forms, this would imply a raised probability for SAA-fibril formation. Hence, we have also performed molecular dynamic simulations of helix-broken and helix-weakened SAA 1−76 fragments interacting with a SK9 segment initially docked to either the N-terminus, the helix-I−helix-I -linker, or the Cterminal region. The Journal of Physical Chemistry B pubs.acs.org/JPCB Article Independent on where the SK9 segment is initially docked, we observe that in all our trajectories, that start from SAA 1−76 in the helix-broken configuration, the virus segment disengages and moves away. Consequently, no noticeable differences to the control simulation are seen. This observation suggests that once the SAA fragment assumes a helix-broken configuration, it will not be affected by the presence of the viral segment and will stay in this less aggregation-prone motif. On the other hand, when the SK9 segment is initially bound to SAA 1−76 in the helixweakened form, we find after about 800 ns a clear signal for forming a N-terminal β-strand involving the first 11 residues. This observation is independent of where the SK9 segment initially binds to the SAA chain, see for instance, the snapshots of representative configurations in Figure 7a . The appearance of a N-terminal β-strand is important as this region is known to be crucial for amyloid formation, 35−37 and we have shown in ref 13 that fibril assembly starts with this region.

As in the earlier discussed simulations starting from a nativelike configuration for the SAA 1−76 fragment, we also observe the β-sheet formation in the helix-I−helix-II linker region or the disordered C-terminal region, see the corresponding snapshots in Figure 7b ,c. Note that while the initial binding sites are obtained from docking calculations, visual inspection of the trajectories shows that the SK9 segment does not stay bound at the initial position.

For this reason, we have also calculated the residue-wise binding probabilities of the SK9 segment toward the helixweakened SAA 1−76 conformations. Data are averaged over the final 3.0 μs of each trajectory and shown in Figure 8 . These binding probabilities indicate that the SK9 segment binds preferentially with hydrophobic (A 10 and A 14 ) and aromatic residues (F 11 and W 18 ) in the N-terminal helix-I region, along with binding to helix-III and the disordered C-terminal region through hydrophobic (A 54 , A 57 , I 58 , and A 61 ), aromatic (W 53 , F 68 , and F 69 ) and electrostatic interaction (E 56 , D 60 , E 74 , and D 75 ). Our data indicate that the binding affinity of the SK9 segment toward the disordered C-terminus is stronger than to the Nterminus. Note that the binding probability of the SK9 segment toward the helix-I−helix-II linker region is relatively low in our simulations, but when binding occurs, the segment forms a βhairpin with the linker region (shown in Figure 7b ).

The increased propensity for β-strand formation is confirmed by Figure 9 , where we plot the secondary structure as a function of time. No β-strands are formed in the control simulations (Figure 7a) , while there is a clear signal for it in Figure 7b where we show the same quantity for a trajectory where the SAA monomer binds initially with the SK9 segment at the helix-I− helix-II linker region.

Hence, a second effect by which SK9 can increase SAA amyloid formation is by shifting the equilibrium toward helix- The Journal of Physical Chemistry B pubs.acs.org/JPCB Article weakened SAA 1−76 fragments and by initiating the β-sheet formation in these fragments. The higher probability for forming an N-terminal β-strand (residues S 2 -S 5 ), known to be the start point for fibril formation, 35, 36, 38 results from the interaction between the N-terminus (residues 1−6) and C-terminal region (residues 65−72), more precisely π-π stacking (F 3 −F 68 and F 4 − F 69 ) and hydrophobic (F 6 −I 65 ) interactions. To specify these relations, we show in Figure 10 the difference between the residue−residue contact probabilities, which measured in simulations of the helix-weakened SAA 1−76 in the presence of a SK9 segment, and the ones measured in the control simulations where no SK9 segment is present. The contact propensities are averaged over the last 3.0 μs of all three trajectories. These relative contact probabilities show that hydrophobic interaction between helix-I and helix-III involving residues A 10 Effect of SK9 Segments on the Stability of SAA Fibrils. In the previous two sections, we have shown that SK9 segments reduce the stability of the hexamer, making it easier to dissolve the hexamer, and that they raise the aggregation propensity of monomer fragments by encouraging N-terminal residues to assume a β-strand configuration that is known to be crucial for fibril formation. 35−38 In our third set of simulations, we finally study the effect of the amyloidogenic segment on the stability of SAA fibrils. This is because the stabilization of the fibril will shift the equilibrium between hexamers, monomers, and fibrils toward the fibrils. Having identified in previous work 17 a twofold-two-layer (2F2L) tetramer as the critical size for fibril stability, we have simulated the effect of SK9 on such tetramers. The creation of the start configuration and the setup of simulations are described in the Methods section. Note that our fibril model is made of segments SAA 2−55 as no fibril model for human SAA 1−76 is available. A second fibril model with the missing residues 56−76 added by homology modeling was no longer considered by us after we realized that the SK9 segment does not bind to the disordered region of residues 56−76. This is because we had seen in previous work 17 that the added presence Following the two systems over 300 ns, we compare the change in stability along the trajectories between the two cases. Representative final configurations are shown together with the corresponding start configurations in Figure 11 . Visual inspection of these configurations suggests that the fibril is stabilized by the SK9 segment.

The above visual inspection is again quantified by following the time evolution of global RMSD (measured over the whole fibril) and chain RMSD in Figure 12 , where the latter one is the average over RMSD measurements of individual chains. Similar to the hexamer, we observe for both systems an initial increase in both RMSD values, which plateaus for the SK9-bound fibril after about 100 ns, while the RMSD continues to rise for the fibril in the absence of the viral segment until reaching a global RMSD of 8.0 Å after 200 ns. The stabilization of the SAA fibril in the presence of SK9 is also supported by the smaller root-meansquare fluctuation, also shown in Figure 12 , and calculated over the last 200 ns, as they indicate that the shape of individual chains changes less than in the control.

The observed stabilization of the chain architecture is interesting as the amyloidogenic segment is binding to the outside of the chains, i.e., the stabilization has to be indirect by either enhancing stacking or packing of chains. Analyzing the contact pattern, we find that the presence of SK9 leads to an increase in hydrogen bonds and hydrophobic contacts involved in the stacking of chains. While at the same time the SK9 segment reduces the contacts and hydrogen bonds involved in the packing of two folds, this loss of contacts is compensated by the new interactions with the SK9 segment. The binding probabilities of SK9 segments toward the SAA fibril shown in Figure 13 indicate that the SK9 segment binds in the N-terminal and the C-terminal region to hydrophobic and aromatic residues. This binding pattern stabilizes residue−residue stacking contacts but also reduces the interstrand packing distance (shown in Table 2 ) through forming contacts with certain charged (E 26 and D 33 ), polar residues (N 27 ), and hydrophobic residues (M 24 and A 30 ) at the packing interface. Here, we define the interstrand packing distance by the centerof-mass distance between two folds at their packing interface (residues 28−31). Note that the data for the control differ slightly from the one listed in ref 17 as our simulations are with 300 ns longer than the 100 ns of the simulations of ref 17. In Figure 14 , we plot the residue−residue stacking contact probability maps of the SAA fibril in the presence (a) or absence The Journal of Physical Chemistry B pubs.acs.org/JPCB Article (b) of the SK9 segment. The difference between the two cases is shown in Figure 14c . Our data indicate that the residue−residue contact probability of the first 21 residues in the N-terminal region, the most hydrophobic and amyloidogenic segment of the protein, increases in the presence of SK9. Especially enhanced are the π−π stacking and the hydrophobic interaction involving residues F 3 , F 4 , L 7 , F 6 , F 11 

The concentration of human serum amyloid A (SAA) in acute COVID-19 patients can grow to levels that in patients with certain cancers or inflammatory diseases may cause systemic amyloidosis as a secondary illness. Hence, SARS-CoV-2 infections may also increase the risk for SAA amyloid formation and subsequent pathologies. However, overexpression of SAA does not always lead to systemic amyloidosis, and mechanisms exist for downregulating SAA concentration and minimizing the risk for amyloidosis. In the present paper, we use molecular dynamics to study how the presence of SARS-CoV-2 proteins may interfere with these protection mechanisms by changing the propensity for forming SAA amyloids. To reduce the computational cost, we have restricted ourselves the nine-residuesegment S 55 FYVYSRVK 63 (SK9) on the C-terminal tail of the SARS-CoV-2-Envelope protein whose location makes it likely to interact with SAA proteins. Our simulations show that SARS-CoV-2 proteins can increase the risk for SAA fibril formation by three mechanisms. First, binding of the SK9 reduces the stability of the biologically active SAA hexamer in which SAA transports lipids during inflammation, shifting the equilibrium toward monomers. Monomers are SAA proteins subject to enzymatic cleavage into smaller fragments, and only these fragments are found in SAA fibrils. Hence, by shifting the equilibrium toward the monomers, the presence of the viral protein segment SK9 increases the risk for fibril formation. This risk is further enhanced by the interaction of SK9 with SAA fragments, which increases the frequency of the aggregation-prone form (called by us helix-weakened) and for this motif raises the propensity to form β-strands, especially for the first 11 residues known to be crucial for fibril formation. Finally, the presence of the amyloidogenic segment SK9 also stabilizes SAA fibrils, moving further the equilibrium toward the fibril, and therefore enhancing the probability for amyloid formation.

Hence, our simulations strengthen our hypothesis that SARS-CoV-2 infections raise the risk for SAA amyloidosis during or after COVID-19. As SAA amyloidosis is characterized by formation and deposition of SAA amyloids in the blood vessels, causing inflammation and thrombosis, it may be behind the broad spectrum of severe and at times life-threatening cardiovascular, gastrointestinal, dermatologic, and neurological symptoms, commonly summarized as multisystem inflammatory syndrome (MIS-C) sometimes observed in COVID-19 survivors. 11 It is interesting to speculate why amyloidogenic regions such as SK9 on SARS-CoV-2 proteins seem to have such a pronounced effect on SAA amyloid formation. One possibility would be that fibril formation is part of the immune response and serves as a way to entrap and neutralize the virus. Such a microbial protection hypothesis 39, 40 has been suggested in the context of Herpes Simplex I infections and the development of Alzheimer's disease. Amyloid formation by SAA may serve a similar role, with SAA amyloidosis would be a consequence of this mechanism becoming overwhelmed. Further work will be needed to test this hypothesis. The corresponding standard deviations of the means are listed in brackets. Figure 14 . Residue-wise interlayer contact probabilities measure in simulations of the two-fold two-layer SAA fibrils in the presence (a) and absence of (b) SK9 segment. A pair of residues is in contact if the center-of-mass distance is less than 7 Å. Data are averaged over the final 200 ns of each trajectory. The difference between the two stacking contact probabilities is shown in (c).

The Journal of Physical Chemistry B pubs.acs.org/JPCB Article

Microvascular Injury in the Brains of Patients with Covid-19

T-Cell Hyperactivation and Paralysis in Severe COVID-19 Infection Revealed by Single-Cell Analysis

The Pathogenesis and Treatment of the ′Cytokine Storm' in COVID-19

The Clinical Course and Prognostic Factors of Severe COVID-19 in Wuhan, China: A Retrospective Case-Control Study

Serum Amyloid A Is a Biomarker of Severe Coronavirus Disease and Poor Prognosis

Laboratory Characteristics of Patients Infected with the Novel SARS-CoV-2 Virus

Amyloid Load and Clinical Outcome in AA Amyloidosis in Relation to Circulating Concentration of Serum Amyloid A Protein

Thrombosis in COVID-19

COVID-19 and Acute Kidney Injury: A Systematic Review and Meta-Analysis. Pathogens

Hypothesis: AA Amyloidosis Is a Factor Causing Systemic Complications after Coronavirus Disease

DorsinviCase Series of Multisystem Inflammatory Syndrome in Adults Associated with SARS-CoV-2 Infection -United Kingdom and United States

Interactions between SARS-CoV-2 N-Protein and α-Synuclein Accelerate Amyloid Formation. bioRxiv 2021

Downregulation, and Aggregation of Serum Amyloid A

A Multiscale Coarse-Grained Model of the SARS-CoV-2 Virion

Self-Assembly of a Nine-Residue Amyloid-Forming Peptide Fragment of SARS Corona Virus E-Protein: Mechanism of Self Aggregation and Amyloid-Inhibition of HIAPP

Structural Mechanism of Serum Amyloid A-Mediated Inflammatory Amyloidosis. Proc. Natl. Acad. Sci

Stability of Human Serum Amyloid A Fibrils

Cryo-EM Fibril Structures from Systemic AA Amyloidosis Reveal the Species Complementarity of Pathological Amyloids

Structure and Drug Binding of the SARS-CoV-2 Envelope Protein in Phospholipid Bilayers

Modeling of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) Proteins by Machine Learning and Physics-Based Refinement

AutoDock Vina: Improving the Speed and Accuracy of Docking with a New Scoring Function, Efficient Optimization, and Multithreading

Comparison of Simple Potential Functions for Simulating Liquid Water

High Performance Molecular Simulations through Multi-Level Parallelism from Laptops to Supercomputers

CHARMM36m: An Improved Force Field for Folded and Intrinsically Disordered Proteins

Polymorphic Transitions in Single Crystals: A New Molecular Dynamics Method

An Analytical Version of the SHAKE and RATTLE Algorithm for Rigid Water Models

LINCS: A Linear Constraint Solver for Molecular Simulations

A Smooth Particle Mesh Ewald Method

Dictionary of Protein Secondary Structure: Pattern Recognition of Hydrogen-Bonded and Geometrical Features

Native Contacts Determine Protein Folding Mechanisms in Atomistic Simulations

End-Point Binding Free Energy Calculation with MM/PBSA and MM/GBSA: Strategies and Applications in Drug Design

Structural Basis for Lipid Binding and Function by an Evolutionarily Conserved Protein

Stability of the N-Terminal Helix and Its Role in

Identification of Regions Responsible for Heparin-Induced Amyloido-The Journal of Physical Chemistry B pubs

Electrostatic Interactions at N-and C-Termini Determine Fibril Polymorphism in Serum Amyloid A Fragments

Herpes Simplex Virus Type I Induces the Accumulation of Intracellular β-Amyloid in Autophagic Compartments and the Inhibition of the Non-Amyloidogenic Pathway in Human Neuroblastoma Cells

Corroboration of a Major Role for Herpes Simplex Virus Type 1 in Alzheimer's Disease. Front

Complete contact information is available at: https://pubs.acs.org/10.1021/acs.jpcb.1c04871

The authors declare no competing financial interest.

The simulations in this work were done using the SCHOONER cluster of the University of Oklahoma, FRONTERA on TACC (under grant MCB20016), or XSEDE resources allocated under grant MCB160005 (National Science Foundation). The authors acknowledge financial support from the National Institutes of Health under research grant GM120634. The authors thank Alan J. Ray for help with the figures.