SARS-CoV-2 sequence data were downloaded from Nextstrain (Hadfield et al. 2018) in January 2023 and 5.06 million high quality sequences were selected for analysis. The N-protein amino acid sequences exhibit ≈43 million instances of mutations distributed across ≈92% of its residues. We have previously characterized this dataset with regard to the amino acid mutational landscape of N-protein, and found mutation frequencies that are strongly dependent on position and largely time-invariant, except for the defining mutations arising in variants of concern, the latter comprising ≈36% Delta-variant and ≈49% Omicron-variant sequences (Schuck and Zhao 2023). A histogram of the number of different amino acids mutations that are found at each residue is shown in Figure 1B. It may be discerned that sequence plasticity is highest in the IDRs, with an average of 5.2 different amino acid mutations compared to 2.9 different mutations on average in the folded domains.

Exploiting the N-protein mutational landscape and sequence data, previous work in our laboratory has focused on local amino acid sequence properties such as mutation effects on transient structural features in the linker IDR (Zhao et al. 2023) and the creation of short linear motifs (Schuck and Zhao 2023). However, nonlocal biophysical properties may also be functionally critical and evolutionarily conserved despite amino acid sequence heterogeneity in IDRs (Zarin et al. 2017; Zarin et al. 2021). The sequence ensembles extracted from the genomic database allow us to ask whether physicochemical properties are constrained or can vary across viable sequences of the mutation spectrum.

To this end, genome data were sorted into unique groups with distinct N-protein amino acid sequences, each sequence carrying a set of distinct mutations that represent a viable N-protein species. For a robust analysis, each mutated sequence was required to be represented in at least 10 different genomes in the database. This led to 6,300 distinct full-length N-protein sequences (N-FL; 1-419). We similarly subdivided the N-protein in different regions (Figure 1A) and grouped unique sets of mutations in each region: For the folded domains we found 720 distinct NTD (N:45-179) and 399 distinct CTD (N:248-363) sequences, while for the IDRs there are 512 N-arm (N:1-44), 1039 linker (N:175-247), and 556 C-arm (N:364-419) sequences. Further subdividing the linker there are 349 distinct sequences for the SR-rich region (N:175-205) and 442 for the L-rich region (N:206-247), respectively. Finally, similarly subdividing the C-arm we obtained the 176 sequences for the N3 region (N:390-419) and 242 for the remainder of the C-arm (N:364-389).

We first examine polarity and hydrophobicity of N-protein and different regions based on their amino acid compositions. As shown in beehive plots of Figure 2, where each of the partially overlapping black dots represents one species from the cloud of mutant sequences, the index values of all N-FL sequences fall within a very narrow range (left column). Properties of the full-length protein may obscure significant differences on a smaller scale, in particular since the polarity and hydrophobicity indices are weighted-average properties. Focusing on folded N-protein modules, we find that hydrophobicity is uniformly high and polarity correspondingly low in the folded NTD and CTD domains, which is consistent with anticipated physicochemical requirements of burying residues in folded structures. By contrast, IDRs exhibit significantly higher polarity and lower hydrophobicity. In particular, the N-arm and C-arm are most polar: despite a very large dispersion across the mutant spectrum, their values do not overlap with those of the folded domains.

It is useful to subdivide the linker IDR further to distinguish the SR-rich region (N:175-205), which exhibits high polarity and low hydrophobicity, from the L-rich region (N:206-247), which exhibits opposite behavior and is among the sequence stretches with lowest polarity values and very high hydrophobicity (Figure 2, red arrows in magenta shaded columns). Despite significant spread across the mutant spectrum, there is no overlap in these properties, which suggests biophysical constraints require the distinct polar and non-polar properties of the SR-rich region and the L-rich region, respectively. Indeed, these regions in the linker IDR have been recognized to play distinct functional roles: The SR-region provides a major hub for phosphorylation, aids in NA-binding, and mediates NA-binding induced allosteric interactions between NTD and the L-rich region (Pontoriero et al. 2022; Yaron et al. 2022; Zhao et al. 2023). This is distinct from the L-rich region, which has a propensity for the formation of transient helices that interact with NSP3 (Bessa et al. 2022), and can assemble via hydrophobic interactions to from coiled-coiled oligomers that contribute to the architecture of RNPs in viral assembly (Adly et al. 2023; Zhao et al. 2023).

Similarly, the C-arm IDR can be subdivided in the N3 region (N:390-419) and the remainder (‘C-arm1’, N:364-389), which also have strikingly different properties (Figure 2, blue arrows in cyan shaded columns): Whereas the connecting C-arm portion is by far the most polar, the N-terminal N3-region is among the most hydrophobic regions of the entire protein. Interestingly, the N3 region contains a transient helix (Cubuk et al. 2021; Zhao et al. 2022; Zhao et al. 2023), which may be involved in recognition of the packaging signal and M-protein interactions localized here (Kuo et al. 2016; Masters 2019). Again, the difference in the physicochemical properties of these regions persists throughout the entire ensemble of sequences despite their significant spread and high mutation frequencies (Figure 1B).

The net charges of the different N-protein regions are displayed in Figure 3A. Similar to polarity and hydrophobicity, viable sequences can have significant spread of net charges among all the mutants, amounting to departures by ±(1-2) from the ancestral sequence. This is expected considering the replacement and introduction of charged residues in the mutational landscape, for example, including those from the defining substitutions of variants. The positive charge of the overall basic protein is shared similarly among all folded domains and IDRs. However, noteworthy is again the contrast arising from subdivision of the linker and C-arm, which displays uneven and non-overlapping distributions: despite the strongly basic character of the linker, its L-rich sequence is nearly neutral; similarly, the basic C-arm splits into an even more basic C-arm1 and an acidic N3 tail region. These differences are highly significant and persist throughout the mutant spectrum.

It is well-established that intracellular N-protein can be heavily phosphorylated (in contrast to N-protein in the virion) (Fung and Liu 2018; Carlson et al. 2020; Johnson et al. 2022; Yaron et al. 2022). As reviewed in (Yaron et al. 2022), most serine, threonine and tyrosine residues in the disordered regions (30 of 37) have been found phosphorylated in different proteomic analyses. Accordingly, we estimated the maximum charge when all of these residues in the IDRs are phosphorylated (Figure 3B). This leads to a negative charge for all IDRs. As might be expected, the largest impact was found in the SR-rich region of the linker, which carries the highest density of phosphorylation sites. Interestingly, despite the considerable spread of net charges within families of mutant sequences, the differences between the regions remain highly significant.

It is noteworthy that the defining mutations of the Delta- and Omicron-variant (denoted by cyan and magenta horizontal lines, respectively) do impact the hydrophobicity, polarity, and charges in all of the N-protein regions. However, their values do not stand out from the clouds of values across the mutant spectrum, which include more extreme values throughout.

The distinct physicochemical properties of the linker and C-arm sub-segments persist throughout the mutant spectrum, which suggests these constitute biophysical constraints for functional SARS-CoV-2 N-protein. Therefore, we asked whether this holds true for N-protein from related coronaviruses such as SARS-CoV-1 (P59595.1), Middle East respiratory syndrome coronavirus (MERS, YP_009047211.1), murine hepatitis virus (MHV, NP_045302.1), human coronavirus NL63 (Q6Q1R8.1), and the 229E-related bat coronavirus APD51511.1. To this end, we used their sequence alignment to SARS-CoV-2 N-protein (shown previously (Zhao et al. 2022)) to subdivide all N-proteins into equivalent regions (Supplemental File S1). As shown in Table 1, the resulting peptides present high sequence identity scores for the FL protein and the folded domains, but, with exception of SARS-CoV-1, have little to no sequence identity in the IDRs. This observation is consistent with the high mutation frequency of the IDRs.

The resulting peptides were subjected to the same analyses of physicochemical properties described above for SARS-CoV-2 N-protein. The results are displayed in Figures 2 and 3 as symbols. With regard to hydrophobicity (Figure 2B), the FL proteins and folded domains show values within the range of the SARS-CoV-2 mutant spectrum. By contrast, more significant spread is observed in most IDR peptides. Nonetheless, the pattern observed for SARS-CoV-2 of hydrophobicity and polarity values of IDRs relative to those of the folded domains, and the pattern comparing subdivisions of the IDRs is closely mirrored for SARS-CoV-1, MERS, and MHV (red symbols). Similar patterns, although with some divergence, are observed for the NL63 and APD51511.1 IDRs (grey pentagrams and diamonds, respectively) which have the least sequence identity to SARS-CoV-2.

Polarity values (Figure 2A) of all coronavirus linker peptides are higher than either FL, NTD, or CTD, and the subdivision of the linker in the peptides corresponding to SR-rich and L-rich regions of SARS-CoV-2 follow the same qualitative trend, with higher polarity in the equivalent SR-rich and lower polarity in the equivalent L-rich peptides for all coronaviruses studied. Similarly, the properties of the equivalent C-arm and subdivision of C-arm1 and N3 peptides for SARS-CoV-1, MERS, and MHV (red symbols) closely track the values from the mutant spectrum of SARS-CoV-2, although this is not the case for the more distant NL63 and APD51511.1 (grey symbols).

Charge properties of related coronaviruses follow a similar pattern of SARS-CoV-2 (Figure 3A), although with somewhat greater differences, particularly again for NL63 and APD51511.1. Peptides corresponding to L-rich regions exhibit low charge, distinctly below those of the SR-rich regions, and similarly, N3 peptides have lower charges than C-arm-1 peptides of the corresponding viral species, and are nearly all acidic. Even though it is unclear to what extent IDRs of other coronaviruses can be phosphorylated, their amino acid composition would provide similar potential as SARS-CoV-2, as the completely phosphorylated charges of all peptides follow closely those of SARS-CoV-2 (Figure 3B).

This suggests that the charge properties and phosphorylation, like polarity and hydrophobicity, of the equivalent IDR sub-regions are functional biophysical constraints maintained across related coronaviruses despite little sequence conservation.

Unfortunately, it is impossible to express and experimentally characterize biophysical properties of all mutant species. Therefore, to assess the range of phenotype variation we examine only six exemplary protein constructs related to variants of concern in comparison with the Wuhan-Hu-1 reference molecule, N_ref (Table 2): 1) N:R203K/G204R with a double mutation in the disordered linker that arose early in the Alpha-variant (B.1.1.7), but occurs also in the Gamma-variant (P.1), and all Omicron-variants (BA.1 through BA.5). It was found to modulate phosphorylation of cytosolic N-protein, enhance assembly in a VLP assay, and increase viral fitness (Johnson et al. 2022; Syed et al. 2022; Syed et al. 2023); 2) N:P13L/Δ31-33 carrying the mutation P13L and the deletion Δ31-33 that are part of the defining mutations of all Omicron variants, with P13L epidemiologically ranked as the most statistically significant N-protein mutation linked to increased fitness (Oulas et al. 2021; Obermeyer et al. 2022); 3) N_o is a combination of N:R203K/G204R and N:P13L/Δ31-33, carrying thereby the complete set of defining mutations of the BA.1 Omicron-variant; 4) N:G215C with a key mutation in the disordered linker that was associated with the rise of the 21J clade of the Delta-variant, and found to modulate a transient helix in the L-rich linker region (Zhao et al. 2022); 5) N:D63G containing another defining mutation of the Delta-variant, located in the NTD and epidemiologically ranked above G215C in increasing SARS-CoV-2 fitness (Obermeyer et al. 2022); and 6) N_δ carrying all four defining mutations D63G, R203M, G215C, D377Y of the Delta-variant. As detailed in Table 2, all of these species are found in the genomic database, and in combination with additional mutations occur in a high fraction of all genomes (exceeding the frequency of the ancestral Wuhan-Hu-1 N-protein by an order of magnitude). However, with the exception of N:G215C, none of the mutants has been studied in detail with regard to their macromolecular biophysical properties.

All mutations considered here are within the IDRs, except for N:D63G, a mutation characteristic of the Delta variant. The effect of the N:D63G mutation in the NTD is highlighted in the shift of the intrinsic fluorescence quantum yield of this mutant in comparison to N_ref (Figure 4A). This may be attributed to changes in the local environment of tryptophan W108, which is partially surface exposed and structurally near the aspartic acid D63, as indicated by AlphaFold structural predictions (Supplementary Figure S1). D63G ablates a negative surface charge near the nucleic acid (NA) binding site of the NTD, which poses the question whether this mutation alters NA binding affinity. We assessed this using sedimentation velocity analytical ultracentrifugation (SV-AUC) with the oligonucleotide T₁₀ as a NA probe. T₁₀ is comparable in length to the NTD binding canyon for NA but does not permit multi-valent binding (Dinesh et al. 2020; Zhao et al. 2021). No significant differences in the intrinsic binding affinity to T₁₀ was detected between N:D63G, other mutants, and the ancestral species (Supplementary Figure S2).

A parameter of great interest from an evolutionary perspective is the thermal stability of the folded domains. This property can be assessed experimentally by differential scanning fluorometry (DSF), which reports on temperature-driven changes in the environment of aromatic amino acids due to changes in solvent exposure (Eftink 2000). Such changes may occur during unfolding or as a result of other conformational changes. In the case of N-protein, conveniently all tryptophan and tyrosine residues of N-protein are located in the NTD and CTD, such that changes in the intrinsic fluorescence report exclusively on changes in the state of the folded domains. As shown in Figure 4B, a major transition is observed with an inflection point at T_i ≈ 49 °C. Compared to the reproducibility of transition temperatures of ±0.1 °C, significant shifts from the ancestral N-protein can be discerned: While Omicron mutations N_o, N:R203K/G204R, and N:P13L/Δ31-33 are neutral, those occurring in the Delta variant (N:D63G, N:G215C, and N_δ) are destabilizing, i.e., they lower the transition temperature. Interestingly, apparent destabilization of the folded domains occurs in N:G215C despite the absence of mutations in the folded domains – 215C being located in the middle of the linker IDR. This non-local mutation effect points to altered intra-molecular interactions between IDRs and the folded domains, and/or changes in contacts between folded domains mediated through an altered oligomeric state. (This is corroborated in non-natural point mutants N:L222P and N:L222P/R226P which abrogate linker helix oligomerization (Zhao et al. 2023) and exhibit T_i-values of ≈51 °C.) the observation of T_i of Furthermore, Figure 4B shows additional transitions occur at higher temperatures broadly in the range of 60 – 70 °C. While their origin is unclear, this signal may accompany the formation of higher-order structure. It is noteworthy that N:G215C is also distinctly different in this feature.

Secondary structure information from the entire molecule including the IDRs can be extracted from circular dichroism (CD) spectra. As may be observed from Figure 4C, significant variation occurs both in the magnitude of the negative ellipticity at ≈200 nm, which mainly reflects disordered residues, as well as the magnitude of the negative ellipticity at ≈220 nm, which reports on helical structure. Compared to the ancestral N_ref, significantly less disorder and greater helicity is observed for N:G215C (and to lesser extent also for N_δ), whereas slightly more disorder is indicated for N:R203K/G204R. Little difference to the ancestral molecule is observed for N_o, N:P13L/Δ31-33, and N:D63G. The absence of significant changes for N:D63G is consistent with this mutation having only a subtle, if any, impact on the NTD conformation. For N:G215C, increased helicity can be attributed to the stabilization of transient helices in the leucine-rich region of the central linker IDR, as shown previously (Zhao et al. 2022; Zhao et al. 2023). N_o consists of the R203K/G204R mutation combined those of N:P13L/Δ31-33. It is noteworthy that the increase in disorder from the R203K/G204R mutation is not exhibited by N_o, while secondary structure of N:P13L/Δ31-33 is very similar to N_o, suggesting some epistatic interaction between mutations of the linker and N-arm.

Tertiary and quaternary structure can be assessed by sedimentation velocity analytical ultracentrifugation (SV-AUC) (Figure 5A). As reported previously, the ancestral N-protein in the absence of NA is a tightly linked dimer sedimenting at ≈4 S (Zhao et al. 2021), but the G215C mutation promotes the formation of higher oligomers via stabilization of coiled-coil interactions of transient helices in the L-rich linker region (Zhao et al. 2022; Zhao et al. 2023). This is consistent with the enhanced helical content of this mutant (Figure 4C). Oligomerization beyond the dimeric N_ref is also observed for N_δ, which incorporates the 215C mutation, but less than for N:G215C. This is consistent with the intermediate helical content of N_δ observed in CD. This may be due to an additional linker mutation of N_δ, R203M, which may limit the stabilization of a helical linker conformation by G215C and thereby weaken the self-association to form higher oligomers.

N-protein has a propensity to form large particles and undergo LLPS, which can be promoted at higher temperatures (Iserman et al. 2020; Zhao et al. 2021). Figure 5B shows the z-average particle size measured by dynamic light scattering (DLS) as a function of temperature. Particle formation is governed by a combination of processes, including the hydrophobicity-driven stabilization of the linker helix and its self-association, ultra-weak interactions across the entire protein contributing to LLPS, and unfolding and aggregation processes. This complicates a comparison of the temperature transitions observed in DSF (Figure 4B) and DLS (Figure 5B) (and a further technical difficulty may be potential differences in temporal lag of conformational rearrangements versus particle assembly kinetics).

Nevertheless, several clear observations can be made. As reported previously, N_ref forms clusters and particles at >55 °C (Zhao et al. 2021), which is strongly enhanced and occurs at a lower temperature for N:G215C, due to the enhancement of the linker oligomerization (Figure 5B) (Zhao et al. 2023). Very similar behavior is observed for N_δ, which suggests that at higher temperatures any inhibitory role of the R203M mutation on self-association may be less relevant compared to G215C. It is interesting to note that, correspondingly, both show a lower T_i in DSF. More moderate enhancement of particle formation is observed for N:D63G, which shows an onset already at ≈50 °C and larger particle averages than the ancestral protein. This also correlates with its significantly lower T_i in DSF. Thus, even subtle structural changes (as shown in Supplementary Figure S1) can impact the assembly behavior.

The opposite effect, strong inhibition of particle formation, is observed for the N:R203K/G204R double mutant. Here, particles form only at temperatures > 70 °C, as a mixture of smaller clusters with some very large aggregates that adventitiously enter the light path in DLS and cause fluctuations in the z-average Stokes radius. Interestingly, although N_o comprises the R203K/G204R mutation, N_o does not share this behavior but instead exhibits slightly enhanced particle formation relative to the ancestral N_ref, comparable to N:D63G. This points to the role of additional mutations in N_o, which besides R203K/G204R features the N-arm mutations P13L and Δ31-33. Interestingly, by themselves in N:P13L/Δ31-33 the particle formation is also suppressed relative to N_ref, although less so than for N:R203K/G204R. This again points to non-additive effects, suggesting that the N-arm IDR mutations interact with the mutations of the linker IDR to promote particle formation of N_o.

To examine the role of the N-arm mutations of the Omicron variants in greater detail we restricted the protein to the N-arm peptide and compared solution behavior of the N-arm constructs N_ref:(1-43) with the Omicron N-arm N:P13L/Δ31-33(1:43), as well as the N-arm with individual mutation N:P13L(1:43) and deletion N:Δ31-33(1:43). Unexpectedly, solutions of N:P13L/Δ31-33(1:43) exhibited elevated viscosity after storage for several days at 4 °C in 20 mM HEPES, 150 mM NaCl, pH 7.5. Since this is a tell-tale of weak protein interactions, we carried out ColabFold structural predictions. Even though ColabFold is trained to predict folded structures, it has been found to be frequently successfully in predicting transient folds in IDRs (Alderson et al. 2022; Zhao et al. 2023). Indeed, it predicts that replacement of proline at position 13 by leucine allows for formation of parallel sheets symmetrically arranged in higher-order N-arm oligomers (Supplementary Figure S3). We proceeded to test oligomerization of the N-arm constructs experimentally in hydrodynamic studies. Figure 6A shows autocorrelation functions of all peptides. While the reference N-arm N_ref:(1-43) and the construct carrying the Δ31-33 deletion behave as expected for non-interacting peptides of this size, the N-arm constructs carrying the P13L mutation (in particular, the Omicron N-arm N:P13L/Δ31-33(1-43)) exhibit very large correlation times. This may be indicative of either formation of large particles or the presence of weak interaction networks as in gels. Similarly, in SV-AUC (Figure 6B) the ancestral reference and the Δ31-33 deletion mutant sediment as expected for non-interacting N-arm peptides (Zhao et al. 2023), whereas rapidly sedimenting, anomalously shaped boundaries with ≈100-fold larger sedimentation coefficient were observed for the Omicron N-arm and the construct carrying solely the P13L mutation. This unequivocally demonstrates the introduction of new protein self-association interfaces from the P13L mutation. They are weak and not apparent in studies of the full-length protein N:P13L/Δ31-33 at low micromolar concentrations, but oligomers can be populated at the ≈100-fold higher achievable concentrations of the peptides, which mirrors the concentration range for in vitro observation of interactions of the leucine-rich linker helices (Zhao et al. 2023).

The ability for N-protein to undergo LLPS in is thought to be crucial for several functions including interactions with stress granules, RNP assembly, and interactions with viral M-protein (Iserman et al. 2020; Savastano et al. 2020; Lu et al. 2021; Carlson et al. 2022; Cascarina and Ross 2022). Weak protein-protein interactions and cluster formation such as shown in Figure 5 and 6 can be coupled to LLPS, or alternatively LLPS may occur independent of clusters following Flory-Huggins theory (Kar et al. 2022). Therefore, we examined the impact of mutations on the propensity for LLPS. As shown in Figure 7 (top left), N_ref readily forms droplets in the presence of T₄₀ oligonucleotides. Under the same conditions N:R203K/G204R (bottom left) does not display droplets, but forms few large particles with fibrillar morphology. In stark contrast, N:P13L/Δ31-33 (bottom center) readily forms droplets that appear to be rapidly merging and growing. The combination of these mutations in N_o exhibits an intermediate propensity for LLPS with droplets in a dispersion of sizes. Mutations comprising the Delta variant (Figure 7 top) have a smaller impact on LLPS than those from the Omicron variant (Figure 7 bottom).

The Wuhan-Hu-1 isolate (GenBank QHD43423) (Wu et al. 2020) was used as the ancestral reference. Sequence data were based on consensus sequences of SARS-CoV-2 isolates submitted to the GISAID as previously described (Zhao et al. 2022; Schuck and Zhao 2023). Briefly, sequence data were downloaded on January 20, 2023 from Nextstrain (Hadfield et al. 2018) and 5.06 million high quality preprocessed sequences were included in the analysis. 746 sequences exhibiting insertions in the N-protein were omitted, as well as those with more than 10 deletions in N-protein and those represented in fewer than 10 genome instances.

The resulting sequence database was parsed for different unique sequences for N-proteins and different segments, using MATLAB (MathWorks, Natick, MA). Sequence hydrophobicity was calculated in RStudio (https://posit.co/) using the package PEPTIDES (Osorio et al. 2015) and polarity and charge using the package ALAKAZAM (Gupta et al. 2015). For maximally phosphorylated charge, -2 was added to the total charge for each serine, threonine, and tyrosine in the IDRs.

Alignment of SARS and related coronavirus sequences (SARS-CoV-1 P59595.1, MERS YP_009047211.1, MHV NP_045302.1, human coronavirus NL63 Q6Q1R8.1, and 229E-related bat coronavirus APD51511.1) was carried out with COBALT at NLM (Papadopoulos and Agarwala 2007), as shown in (Zhao et al. 2022). This alignment was used to dissect related viruses into regions corresponding to the SARS-CoV-2 regions (N-arm, NTD, linker, SR-rich, L-rich, CTD, Carm, Carm1, N3). The resulting segments of the related viruses were subjected to analysis of physicochemical properties as described above. Sequence similarity of the corresponding regions relative to the SARS-CoV-2 regions was calculated using BLAST blastp suite (Altschul et al. 1997), using an expectation threshold of 0.9, word size 2, and BLOSUM63 scoring matrix.

Structural predictions for NTD and N-arm were carried out using ColabFold (Mirdita et al. 2022) and graphics were generated using ChimeraX (Pettersen et al. 2021).

N:D63G and N:G215C were purchased from EXONBIO (catalog# 19CoV-N170 and 19CoV-N180, San Diego, CA), while N_ref, N:R203K/G204R, N:P13L/Δ31-33, N_o, and N_δ were expressed in house as described previously (Zhao et al. 2022; Zhao et al. 2023). Briefly, the full-length protein with an N-terminal Tobacco Etch Virus (TEV) cleavage site and 6His tag was cloned into the pET-29a(+) expression vector and transformed into One Shot BL21(DE3)pLysS E. coli (Thermo Fisher Scientific, Carlsbad, CA). After cell lysis, the protein was bound to a Ni-NTA column, and unfolded and refolded to remove residual protein-bound bacterial nucleic acid (Carlson et al. 2020). After elution the 6xHis tag was cleaved and the protein purified by size exclusion chromatography. Greater than 95% purity of the proteins was confirmed by SDS-PAGE, and the ratio of absorbance at 260 nm and 280 nm of ∼0.50-0.55 confirmed absence of nucleic acid. For a subset of mutants, the protein sequence and mass were tested and confirmed by LC-MS/MS and LC-MS mass spectrometry, respectively. Biophysical experiments were preceded by dialysis in either high-salt buffer (HS) consisting of 20mM HEPES, 150mM NaCl, pH 7.5, or low-salt buffer (LS) consisting of 10.1 mM Na₂PO₄, 1.8 mM KH₂PO₄, 2.7 mM KCl, 10 mM NaCl, pH 7.4.

The oligonucleotide T₄₀ was purchased from Integrated DNA Technologies (Skokie, IL), as purified by HPLC and lyophilized. N-arm peptides were purchased from ABI Scientific (Sterling, VA), as purified by HPLC, examined by MALDI for purity and identity, and lyophilized.

CD spectra were acquired in a Chirascan Q100 (Applied Photophysics, U.K.), using cuvettes of 1 mm pathlength, and data acquisition with 1 nm steps and 1 sec integration time. Results are averages of 3 acquisitions, corrected for buffer background. Protein concentration was 3 µM in buffer LS, except N_o in buffer HS.

For the acquisition of fluorescence spectra, protein samples at 1 µM were loaded into a quartz cuvette with 1.0 cm optical pathlength. Steady-state tryptophan (Trp) fluorescence emission spectra in the range from 305 nm to 500 nm were recorded in a spectrofluorimeter (QuantaMaster, Photon Technology) with excitation at 295 nm using a 1.0 nm increment. Each data point is an average of 10 accumulations.

DSF was carried out in a Tycho instrument (Nanotemper, Germany) as previously described (Zhao et al. 2021). Briefly, 10 µL samples were aspirated in capillaries (TY-C001, Nanotemper, Germany), and intrinsic fluorescence was measured at 350 nm and 330 nm while the temperature was ramped from 35 °C to 95 °C at a rate of 30 °C/min. The first derivative of the intensity ratio was calculated as a function of temperature. DSF experiments were carried out at protein concentrations of 2 µM in buffer LS, except for N:R203K/G204R which was measured in buffer HS. As a buffer control, the difference in T_i for N_ref in LS and HS buffer was measured and found to be within error of data acquisition.

SV-AUC experiments were carried out in a ProteomeLab XL-I analytical ultracentrifuge (Beckman Coulter, Indianapolis, IN) in standard configurations (Schuck et al. 2015). Briefly, 2 µM protein samples were filled in cell assemblies composed of charcoal-filled Epon double-sector centerpieces with sapphire windows, inserted in an 8-hole AN-50 TI rotor and temperature equilibrated. After acceleration to 50,000 rpm data acquisition commenced using the absorbance optical detector at 280 nm and the interference optical detector. Data were analyzed in SEDFIT (sedfitsedphat.nibib.nih.gov/software) in terms of a sedimentation coefficient distribution c(s) (Schuck 2016). Nucleic acid binding experiments were analyzed with isotherms of signal weighted-average sedimentation coefficients in SEDPHAT (Schuck and Zhao 2017). For studies of the N-arm peptide species, 400 µM peptide samples were studied by gravitational sweep sedimentation using rotor speed steps of 3,000 rpm, 10,000 rpm, 40,000 rpm, and 55,000 rpm (Ma et al. 2016) and analyzed with a model for apparent sedimentation coefficient distributions ls-g*(s) (Schuck 2016) as a qualitative representation of rapidly migrating boundaries of N:P13L(1:43) and N:P13L/Δ31-33(1:43), or with c(s) distributions for N_ref:(1:43) and N:Δ31-33(1:43).

Temperature-dependent DLS autocorrelation data of N-protein species were collected in a NanoStar instrument (Wyatt Technology, Santa Barbara, CA) equipped with a 658 nm laser and using a detection angle of 90°. 100 µL samples at 3 µL N-protein in LS buffer were inserted into a 1 µL quartz cuvette (WNQC01-00, Wyatt Instruments), with excess sample to prevent evaporation in the observation chamber. A temperature ramp rate of 1°C/min was applied with 5 sec data acquisitions and averaging 3 replicates for each temperature point. Data were collected and processed with the software Dynamics 7.4 (Wyatt Instruments).

DLS studies of N-arm peptides were carried out in a Prometheus Panta (Nanotemper, Germany) instrument at 20°C. The samples were loaded into a capillary (Nanotemper PR-AC002) and ACFs were acquired using the 405 nm laser at the detection angle of 140°.

Optical imaging of in vitro phase-separated condensates was carried out as described previously (Zhao et al. 2021). For reaction mixtures, N protein and T₄₀ in buffer LS were combined and mixed immediately prior to imaging. 20 µL samples were transferred onto a glass-bottom 35 mm dish (catalog # Part No: P35G-1.5-20-C, MatTek) for imaging at room temperature. Images were acquired on a Nikon Ti-E microscope equipped with a 100X 1.49 NA oil objective lens (LIDA light engine, Lumencor, Beaverton, OR) and recorded with a Prime 95B camera (Teledyne Photometrics) with a pixel size of 110 nanometers. Images were background subtracted and contrast enhanced using MATLAB (Mathworks, Natick, MA).