Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.
Advertisement
Nature Communications volume 13, Article number: 5775 (2022)
1639
2
1
Metrics details
Elucidating how individual mutations affect the protein energy landscape is crucial for understanding how proteins evolve. However, predicting mutational effects remains challenging because of epistasis—the nonadditive interactions between mutations. Here, we investigate the biophysical mechanism of strain-specific epistasis in the nonstructural protein 1 (NS1) of influenza A viruses (IAVs). We integrate structural, kinetic, thermodynamic, and conformational dynamics analyses of four NS1s of influenza strains that emerged between 1918 and 2004. Although functionally near-neutral, strain-specific NS1 mutations exhibit long-range epistatic interactions with residues at the p85β-binding interface. We reveal that strain-specific mutations reshaped the NS1 energy landscape during evolution. Using NMR spin dynamics, we find that the strain-specific mutations altered the conformational dynamics of the hidden network of tightly packed residues, underlying the evolution of long-range epistasis. This work shows how near-neutral mutations silently alter the biophysical energy landscapes, resulting in diverse background effects during molecular evolution.
Addressing the relationship between mutations and their effects on the functional landscape can provide a detailed evolutionary trajectory of a protein1. However, nonadditive interactions between mutations make comprehensive mapping of this sequence-function space extremely difficult to achieve, if not impossible2.
For example, the effects of two single mutations at sites a and b, denoted as ΔGa and ΔGb, respectively, can cause the double mutant of the two sites to show either additive (ΔΔGa,b = ΔGa + ΔGb) or nonadditive behaviors (ΔΔGa,b ≠ ΔGa + ΔGb)3. Epistasis refers to the nonadditive effect between two mutations. Various patterns of epistatic interactions are possible, and the same mutation can result in different outcomes in different backgrounds3,4. Epistasis also allows neutral or near-neutral mutations to have an evolutionary impact by diversifying genetic backgrounds, which might serve as a pre-adapted platform in the new environment5,6,7. Alternatively, these mutations may have permissive roles in which subsequent mutations would otherwise be deleterious8,9.
Although many studies on epistasis have focused on the mutational effects on the functional landscape, mapping the mutational effect on the energy landscape is critical for understanding the mechanistic basis of epistasis10,11. In particular, because neutral mutations are functionally silent on their own, deep mechanistic studies on how they alter biophysical traits of a protein are essential6,8.
From the physical perspective, epistasis can be classified as a short- or long-range mutational effect. While a short-range epistatic interaction between directly contacting residues is structurally intuitive, the mechanistic basis of long-range epistatic interactions between spatially distant residues requires further refinement, although some interesting mechanisms are available12,13,14,15. Moreover, epistatic effects of neutral mutations often occur through long-range interactions8,16. Thus, the biophysical foundation of long-range epistasis can help understand how neutral mutations contribute to protein evolution9.
In this study, we trace evolutionary changes in the in vitro functional and biophysical traits of nonstructural protein 1 (NS1) of influenza A viruses (IAV). To this end, we selected four NS1s of human IAV strains that emerged between 1918 and 2004: 1918 H1N1, Puerto Rico 8 (PR8) H1N1, Udorn (Ud) H3N2, and Vietnam (VN) H5N1.
NS1 of an influenza virus is a multifunctional virulence factor that antagonizes the host’s innate immune response and is associated with many strain-specific functions in viral replication and host range17,18. Moreover, NS1 is one of the most frequently mutated19 and mutation-tolerant proteins20 in the IAV genome. Multiple strain-specific mutations in NS1 were identified as adaptive mutations with respect to viral virulence and replication21,22; however, many remain uncharacterized, mainly because of the lack of functional phenotypes. Especially, many mutations occur in the effector domain (ED) of NS1 that interacts with a number of host factors, including retinoic acid-inducible gene 1 (RIG-I), tripartite motif containing 25 (TRIM25), and phosphoinositide 3-kinase (PI3K), during the infection cycle23,24,25,26.
Outlining the present study, we first present that strain-specific mutations in NS1 EDs have near-neutral effects on the association kinetics to the p85β subunit of PI3K. We subsequently find that near-neutral, strain-specific mutations have long-range epistatic interactions with the residues on the p85β-binding interface of NS1. Using isothermal titration calorimetry (ITC), we then find that strain-specific mutations altered the thermodynamic energy landscape of the NS1:p85β interaction. To further address how strain-specific mutations altered the energy landscape of NS1s and enabled long-range epistasis during evolution, we conducted extensive biophysical analyses of NS1s. Structural analyses revealed conformational variations occurred during NS1 evolution. Especially, using NMR spin dynamics of backbones and side-chains, we reveal how strain-specific mutations influence long-range epistasis through dynamic reorganization of packing interactions in the hydrophobic core. Consequently, our study provides a high-resolution mechanism by which neutral mutations contribute to protein evolution by reshaping the energy landscape, despite the absence of functional changes.
The four NS1 EDs (referred to as NS1s) employed in this study contain 3–14 % of residues that are strain-specifically mutated in the protein sequence with respect to the 1918 NS1 (Supplementary Fig. 1a). Many of the mutations are conservative between similar types of amino acids located at buried or partially buried sites in the NS1 structure; the average relative solvent accessible surface area of mutated residues was 32.2% (Fig. 1a, b and Supplementary Fig. 1b).
a The sequence logo of four NS1s. Completely conserved residues are in gray. Mutated residues are color coded according to the relative solvent accessible surface area. Asterisks and black filled circles represent the core and rim interface residues, respectively. b (left panel) Positions of strain-specific mutations are shown as ball-and-sticks. Partially or completely buried residues are shown in blue. Solvent exposed residues are shown in light blue. (right panel) The p85β-binding interface is shown in surface model. The bound p85β is shown green. c BLI-based kinetic parameters for binding between NS1s and p85β. The uncertainties correspond to the standard deviation of three repeated measurements.
To study the functional effect of strain-specific mutations, we first confirmed that all four NS1s bind to the full-length PI3K using co-immunoprecipitation (Supplementary Fig. 2a). Next, we characterized the interaction between NS1 and p85β using BLI (Biolayer interferometry). Compared to the 1918 NS1, other NS1s dissociated from p85β with a smaller dissociation rate constant, resulting in a moderately higher binding affinity to p85β (Fig. 1c and Supplementary Figs. 2b–e).
However, there was no difference in the binding among the later strains of 1918 IAV; namely, PR8, Ud, and VN NS1s shared virtually identical binding (kon) and unbinding (koff) rate constants (Fig. 1c). Notably, the kon value was highly similar in all four NS1s, including 1918 NS1, showing that strain-specific mutations across the NS1s have a near-neutral effect on the p85β-binding rate constant. These results also suggest that kon is a biophysical constraint of NS1 evolution27,28.
We further tested whether the energetic contribution of the individual p85β-binding residues to kon is conserved across different NS1 backgrounds (Supplementary Fig. 3). The p85β-binding interface on NS1 consists of 20 residues based on the crystal structure of the complex23,29, while our previous study identified residues 89, 95, 98, 145, and 146 of NS1 as the major players determining the p85β-binding rate constant29; thus, these residues are referred to as core interface residues (Fig. 2a).
a The core and rim interface residues are shown as cyan sticks and wheat spheres, respectively. b Definition of ΔGon (WT–MT) in the binding reaction coordinate between NS1 and p85β. Negative ΔGon (WT–MT) indicates deleterious mutational effects on association kinetics (kon). c Strain-specific ΔGon of core interface residues. d Evolution of ΔGon of core interface residues as a function of years of the occurrence of IAV strains. e–g Epistatic interactions (ΔΔGepi) between strain-specific mutations (A → a) and mutations of individual core interface residues (B → b). ab* corresponds to the experimentally measured effects of double mutations, while ab corresponds to the calculated additive (non-epistatic) effects. Red and faint blue arrows correspond to the experimentally measured and calculated additive effects of mutations, respectively. The difference between ab* and ab defines the pattern of epistasis. Epistatic interactions of Y89F mutation between e 1918 and PR8, f PR8 and Ud, and g Ud and VN NS1 backgrounds. c–g All bars or filled symbols represent ΔGon or ΔΔGon values calculated using mean kon values. Error bars represent the propagated standard deviation of three repeated measurements of kon values. Source data are provided as a Source data file.
Using Ala-scanning mutagenesis30, we measured the extent to which individual core interface residues contribute to the p85β-binding rate constant in the four NS1 backgrounds. For Y89, however, we incorporated F89 instead of A89 because the Y89A mutation reduced binding too much to quantify the effect. We first calculated ΔGon (WT–MT) of the individual mutants with respect to their wild-type constructs (Fig. 2b); a larger negative number means a larger deleterious effect of a mutation. Notably, individual core interface residues have substantially different binding energetics in different NS1 backgrounds (Fig. 2c). Next, to follow the evolutionary trajectory of the binding energetics at the residue level, we plotted the ΔGon of the core interface residues along with the year of occurrence of IAV strains (Fig. 2d). Overall, the energetic contribution of core residues decreased from 1918 to 1972 (Ud) and then increased in 2004 (VN).
These background-dependent mutational effects suggest that the core binding residues are energetically coupled with strain-specific mutations of NS1, i.e., epistatic interactions. So, we measured the epistatic interactions using a thermodynamic cycle analysis (Fig. 2e–g and Supplementary Fig. 4). Note that we considered a different NS1 background as a single mutant, although multiple mutations are involved in different NS1s. Thus, energetic coupling (i.e., epistasis) in our study represents the interaction between mutations of individual core interface residues and the collection of strain-specific residues. Following the standard notation for epistasis4, here, we define that additive and negative epistasis corresponds to less (i.e., larger kon values) and more (i.e., smaller kon values) deleterious effects, respectively, by a double mutation than the sum of two single mutational effects. Sign epistasis means that the effect of a double mutation is beneficial with respect to one type of single mutation background and deleterious with respect to another single mutation background.
Indeed, we observed various types of epistasis, including additive, negative, and sign epistasis Supplementary Fig. 4). Moreover, the same residue exhibited different types of epistasis depending on the choice of NS1 background. For example, Y89F showed sign epistasis between the 1918 and PR8 NS1 backgrounds but negative epistasis between Ud and VN NS1 backgrounds (Fig. 2e–g).
Then, to examine the evolutionary trajectory of epistasis, we calculated how the epistatic interaction evolved in different NS1s relative to the 1918 NS1 (Supplementary Fig. 5). This trajectory is equivalent to cumulative epistatic interactions starting in 1918. Overall, core interface residues showed positive epistatic interactions with strain-specific mutations between 1918 and 1972, indicating that the binding interface mutations have less deleterious effects on kon values of NS1s of later IAV strains. This trend reversed between 1972 and 2004, with large negative epistasis. As a result, overall epistatic interactions seemed to fluctuate around null epistasis. This outcome indicates that strain-specific mutations have near-neutral effects on the overall association kinetics, while they altered the energetic contribution of individual residues during NS1 evolution. Structurally, however, the majority of strain-specific mutations are located remotely from the p85β-binding interface, indicating long-range epistasis (Fig. 1b).
Long-range epistasis—energetic coupling between distant residues—suggests that strain-specific mutations altered the protein energy landscape3,4. The protein energy landscape can be probed by examining biophysical properties such as binding thermodynamics31. Thus, we characterized the thermodynamic signature for the binding of NS1s to p85β using isothermal titration calorimetry (ITC). Here, changes in binding thermodynamic signature across NS1s represent the evolution of the binding energy landscape accompanying changes in enthalpic (ΔH) and entropic (TΔS) properties.
Overall, binding was driven by favorable enthalpy and entropy changes (Fig. 3a and b). However, Ud NS1 showed a striking difference from other NS1s in that the binding was not accompanied by measurable heat; that is, the binding was driven exclusively by a favorable entropy change (Fig. 3a, b). Although minor compared to Ud NS1, other NS1s also showed sizable variations in thermodynamic signature, indicating strain-specific variations in the binding energy landscape.
a ITC thermograms and binding isotherms for the interaction between NS1s and p85β. Red lines represent fit curves using 1:1 binding stoichiometry. b Estimated thermodynamic parameters for binding between NS1s and p85β. Uncertainties represent the standard deviation of three repeated measurements. N.D. represents not detected because of the absence of heat in the binding reaction. c Surface representation of NS1. Core interface residues are shown as sticks. Residues 95 and 98 of 1918 and Ud NS1s are shown in cyan and magenta, respectively. d The structure of p85β bound to NS1 is shown to illustrate the interaction with the core interface residues.
This result is surprising because the protein sequence of the p85β-binding interface is almost completely conserved across the NS1s. The only difference is that Ud NS1 contains mutations at two of the 20 interface residues compared to the other three NS1s. However, these mutations are isosteric: L to I and M to L for residues 95 and 98, respectively (Fig. 3c, d). Moreover, grafting the Ud sequence to 1918 NS1 did not affect its binding to p85β23.
To further confirm that the two mutations are not responsible for the drastic change in the thermodynamic signature, we incorporated I95L/L98M mutations into Ud NS1; i.e., this mutant contains the same sequence in the binding interface as other NS1s. ITC data showed that the binding of Ud-I95L/L98M to p85β is not accompanied by heat as in wild-type Ud NS1 (Fig. 3a, b). We also tested the effect of grafting the Ud sequence to the PR8 background (i.e., PR8 L95I/M98L) and found no significant change in the thermodynamic signature for binding to p85β compared to that of wild-type PR8 NS1 (Fig. 3a, b). These results demonstrated that the direct binding interface is not responsible for the variation in the thermodynamic signature for binding. Instead, our results suggested that strain-specific mutations reshaped the energy landscape of NS1 via long-range interactions, resulting in strain-dependent thermodynamic signatures for binding to p85β. However, its underlying mechanism remains to be determined. In particular, future studies on the NS1-p85β complex are warranted to reveal the mechanism.
These results are reminiscent of the recent finding that proteins explore diverse thermodynamic mechanisms by neutral or selective mutations during adaption to high temperatures, referred to as thermodynamic systems drift32. Thus, our results support the notion that proteins can acquire the same functional phenotype using completely different biophysical mechanisms32,33,34.
To reveal the mechanistic bases of long-range epistasis and how strain-specific mutations reshaped the energy landscape of NS1 during evolution, we first examined the structure of the NS1:p85β complexes. Noticeable changes in the complex structures could be a plausible explanation for the strain-specific epistasis and variation in the thermodynamic signature. The structure of the NS1:p85β complex was available for 1918 and PR8 NS1s23,35, and we determined a new crystal structure of VN NS1 bound to p85β in this study (Supplementary Table 1).
Interestingly, the conformations of NS1s in the complexes were almost identical, with RMSD < 0.5 Å (Fig. 4a, b). Although the complex structure of Ud NS1 is not available, it is reasonable to assume that it has a structure similar to other NS1s, considering the sequence identity in the binding interface. This result indicated that the complex structure is not responsible for strain-specific mutational effects on NS1.
a RMSD (n = 4) of backbone Cartesian coordinates of NS1s. Source data are provided as a Source data file. See panel e for PDB IDs of NS1 structures included in the calculation. b Superimposed structures of 1918 (cyan), PR8 (orchid), and VN (slate blue) in complex with p85β. The structure of p85β is not shown for clarity. c Superimposed structures of free NS1s (n = 13) from 1918, PR8, Ud, and VN strains. The conformationally most heterogeneous region is marked with a red dotted circle. d Schematic showing the binding reaction between NS1 and p85β. BI and BC conformers are shown in orchid and pale blue, respectively. e The result of hierarchical cluster analysis of free NS1 structures (n = 16). f Energy landscapes illustrate the conformational exchange between BI-like and BC-like states of NS1s. Exchanging rate constants (kex) and population of BI conformers were determined by NMR CPMG RD (see Supplementary Fig. 7). N.D. not detected.
Given the structural invariance of the complexes, we hypothesized that the structure and dynamics of NS1 in the unbound state are most likely responsible for the strain-specific epistasis and thermodynamic signature for binding. The following observations further support this hypothesis. The observed epistasis represents the change in the free-energy difference between the unbound and transition states (Fig. 2b), and the ITC result corresponds to the free-energy difference between the unbound and bound states. Therefore, both results are directly related to the structure and dynamics of free NS1.
To study the evolutionary influence on the conformation of free NS1, we first compared all the crystal structures of free NS1s (n = 16) of the four IAV strains available in PDB. The structural alignment showed substantial conformational heterogeneity among free NS1s (Fig. 4a, c). Previously, we showed that 1918 NS1 populates two conformations23: p85β binding-competent (BC) and -incompetent (BI) conformers (Fig. 4d). Indeed, hierarchical cluster analysis showed that all 16 crystal structures could be clustered into either BI-like or BC-like conformations (Fig. 4e).
However, which of the BI and BC states represents the major conformer in the solution remained uncertain. To identify the major conformer of free NS1s, we compared NMR NOESY cross-peaks specific to either BI-like or BC-like conformers. Based on the NOESY cross-peaks pattern, we concluded that the free NS1s of the four IAV strains populate the BC-like conformer as the major species (Supplementary Fig. 6).
Considering the presence of the two groups of crystal structures, it is likely that the BC state of NS1 undergoes conformational exchange with the BI state. We further characterized the BI-BC conformational dynamics using the NMR 15N CPMG relaxation dispersion36,37,38. Indeed, PR8 and VN NS1s underwent the BI-BC dynamics on a sub-ms timescale (Fig. 4f and Supplementary Fig. 7). Our previous study showed that 1918 NS1 also undergoes the BI-BC conformational dynamics on a similar timescale23. The analysis of slow-exchanging residues on the NMR timescale revealed that the population of the BC-like conformer was ~99% for PR8 and VN NS1s and ~90% for 1918 NS1 (Supplementary Fig. 7).
On the contrary, Ud NS1 was static and only populated the BC conformer (>99%) (Fig. 4f)23, which is consistent with the fact that all available Ud NS1 structures showed only a BC-like conformer. Although Ud NS1 could be dynamic in different timescales that are not detectable by NMR CPMG-RD, it was obvious that strain-specific mutations reshaped the energy landscape of free NS1s, as evidenced by altered conformational dynamics (Fig. 4f).
Given that the epistatic interactions were probed based on the mutational effects on binding rate constants, we expected that the energy landscapes of the BC state differ strain-specifically; the BC conformer is the major species and responsible for binding to p85β. Therefore, we examined the conformational variations in the BC state between NS1s.
To this end, we compared NMR chemical shifts (CS) of the backbone amide (1H and 15N), 13Cα and 13Cβ resonances in four NS1s (Fig. 5a and Supplementary Fig. 8). The root-mean-square deviation of chemical shift (RMSD-CS) among the NS1s reports the evolutionary variation in the conformation of the BC state. In other words, conformational fluctuations of NS1 during evolution can be mapped in a site-specific manner using RMSD-CS. Residues with an RMSD-CS larger than the third quartile Q3 were identified as structurally variable residues across NS1s (Fig. 5a, b). Here, we excluded strain-specific mutations from the calculation of Q3; note that the inclusion of mutated residues altered the selection of only five residues and did not affect our analysis.
a RMSD-CS of backbone amides (15NH) resonances of NS1s (n = 3 or 4). Green bars correspond to strain-specific mutational sites. Dotted line represents Q3. Orange arrows show positions of core interface residues. b Residues whose RMSD-CS (15NH) > Q3 are shown as light blue spheres. Core interface residues are shown as orange sticks. c Strain-specific mutation sites are shown as spheres. Mutations in the hydrophobic core are shown in blue spheres. d RMSD-CS of methyl (13CH3) resonances of NS1s. Note that x-axis represents approximate residue number due to space limitation. The same plot with explicit residue numbers is shown in Supplementary Fig. 10. a, d Source data are provided as a Source data file. e Residues whose RMSD-CS (13CH3) > Q3 are shown as purple spheres. The β-sheet shared by the hydrophobic core and p85-binding surface is shown in cyan. f Conformationally variable methyl-containing side-chains and β6-α2 loop are shown in magenta and blue spheres, respectively, in the complex with p85β (PDB ID: 6U28).
Many of the identified residues are spatially close to the mutation sites (Fig. 5a–c), indicating that the chemical shift changes represent mutational effects on the conformation of nearby residues. Moreover, the rim interface showed conformational variation across NS1s; for example, several residues in the β6-α2 loop in the rim interface showed conformational variation (Fig. 5b). This result is also consistent with the crystallographic analysis of BC conformers (Fig. 4a; blue line).
To further examine the side-chain conformational changes, we compared RMSD-CS of side-chain methyl resonances (1H and 13C) of Ile, Leu, Val, Thr, Met, and Ala (ILVTMA) in four NS1s (Fig. 5d and Supplementary Figs. 9 and 10, and Supplementary Tables 2–5). Notably, all conformationally variable side-chains (i.e., RMSD-CS > Q3) clustered in the hydrophobic core (Fig. 5e). The hydrophobic core is spatially located on the opposite side of the β-sheet harboring core p85β-binding residues (Fig. 5e, f). Moreover, the hydrophobic core contains multiple strain-specific mutation sites (Fig. 5c). This result indicates that the hydrophobic core might play a key role in transmitting mutational effects to the p85β-binding interface through a long-distance network of physical contacts. Notably, none of the core interface residues showed considerable conformational variation across NS1s, while the hydrophobic core and rim interface (i.e., β6-α2 loop) were sensitive to strain-specific mutations (Fig. 5f).
The conformational variations in the hydrophobic core and rim interface across NS1s can result from diverse biophysical mechanisms32,34. Thus, we further examined the mechanistic basis of the long-range mutational effects on the rim interface.
We first tested whether the evolutionary conformational variation of NS1 is the result of changes in the overall protein stability. Long-range epistasis can be mediated by differential mutational effects on protein stability4,11,39,40. We measured the protein stability of NS1s using H/D-exchange coupled with NMR and the thermal shift analysis. Intriguingly, all four NS1s showed similar global protein stability (Supplementary Fig. 11), indicating that the long-range epistasis is not due to the differential stability of NS1s. These results also suggest that the stability of NS1s was already optimized during natural evolution.
Protein motion is another commonly postulated mechanism of epistatic interactions14,15,41; however, there is a paucity of direct experimental evidence. To test the role of protein motion in long-range epistasis, we hypothesized that mutations affected the rim interface by altering the flexibility of hydrophobic core residues of NS1. However, it should be noted that our hypothesis does not assume a contiguous pathway between a mutational site to the rim interface. Instead, we expected that the flexibility of individual residues changes differently depending on the position and nature of the mutations. For example, a comparison between strain-1 and strain-3 (Fig. 6) can show the long-range interaction between residues A and F. However, this comparison does not reveal other hidden relationships between residues E and F, for which strain-1 and strain-2 need to be compared. Moreover, the same residue may become more flexible in one NS1 but more rigid in another (Fig. 6). Taken together, studying the residue-specific flexibility within a specific NS1 might be insufficient to identify all the residues responsible for long-range epistasis.
Interactions between residues A and F may not be detectable in individual NS1s, as strain-specific mutations (denoted by a, b, and c) differently perturb the flexibility of individual residues. However, residues with large changes in flexibility (red balls) during evolution can be identified by comparing the standard deviation of flexibility among NS1s. The identified residues might collectively form a network in the hydrophobic core and rim interface when they are projected onto the NS1 structure. Our hypothesis suggests that the networked residues interact with each other through direct or indirect contact during NS1 evolution.
However, a large statistical variation in the residue-specific flexibility among NS1s indicates that the residue underwent significant changes in flexibility during NS1 evolution. Therefore, we can accept the hypothesis if the residues with a large statistical variation in flexibility are clustered in the hydrophobic core and rim interface (Fig. 6).
To characterize mutational effects on the flexibility of individual residues of NS1s, we measured the NMR order parameter (S2) of the backbone amide bonds (1H-15N)42. Generally, NMR S2 characterizes the energy landscape shaped by the conformational flexibility of individual chemical bonds. Here, S2 probes the energy landscape basin corresponding to the BC-like conformer because it populates 90–99% of the entire population (Fig. 7a). For comparison, the NMR CPMG-RD probed the energy landscape shaped by conformational exchange in the slower (typically μs – ms) timescale (Fig. 7a).
a Schematic showing the energy landscape mapping by NMR CPMG-RD and S2 analyses. b Standard deviation (SD) of backbone S2NH values acquired for four NS1s. Red and blue bars correspond to core and rim interface residues, respectively. Dotted line represents Q3. c Residues whose SD-S2NH > Q3 are shown as blue spheres. Strain-specific mutational sites are shown in yellow spheres. Core interface residues are shown in red sticks. d SD of methyl axis order parameter (S2axis) acquired for four NS1s. b, d The exact sample size (n) for individual residues varies owing to peak overlap and variations in peak intensities. The sample sizes for individual residues are provided as a Source data file. Dotted line represents Q3. Note that x-axis corresponds to approximate residue number due to space limitation. The same plot with explicit residue number is presented in Supplementary Figs. 14 and 15. e NS1 structure with residues whose SD-S2axis > Q3 are shown as tan spheres. Methyl groups examined by NMR are shown in purple. β-sheet of NS1 is shown in cyan. f All dynamically variable residues are shown as spheres. Residues with SD-S2NH > Q3 and SD-S2axis > Q3 are shown in light blue and purple, respectively. Residues whose backbone and side-chain are both dynamically variable are shown in orange. g All dynamically variable residues are shown in spheres. Residue 171 is shown in green. Rim interface residues that are dynamically variable are shown in slate blue with residue numbers in purple. Core binding residues are shown in orange sticks. Source data are provided as a Source data file.
Residues with a large standard deviation of S2 values (i.e., SD-S2 > Q3) across NS1s were identified as dynamically variable residues; that is, the conformational flexibility of these residues fluctuated widely during the evolution of NS1 (Fig. 7b and Supplementary Fig. 12). Overall, the backbone amide order parameter (S2NH) analysis largely recapitulated the mutational effects on the structure: the core interface residues only showed slight variation in S2NH across NS1s, whereas residues with a large SD-S2NH were mainly in the hydrophobic core and rim interface (Fig. 7c). Noticeably, the dynamically variable residues are mainly located at or near the mutational sites and are mainly in the well-folded secondary structures, not in the flexible loops (Fig. 7c). This result suggests that the observed variation in dynamics reflects the bona fide mutational effects and not the variation of intrinsically flexible regions by subtle changes in experimental conditions. Thus, site-specific conformational dynamics most likely varied along with the evolution of NS1.
Interestingly, three hydrophobic core residues with a dynamically variable backbone also showed conformational variations in the side-chain based on methyl chemical shift changes (Supplementary Fig. 13). Generally, however, backbone and side-chain dynamics do not necessarily correlate, especially for residues with a long side-chain43. So, we hypothesized that the high packing density in the hydrophobic core enabled the facile propagation of mutational effects to a longer range, and as a result, both backbone and side-chain dynamics probe mutational effects in a complementary manner.
To test the idea, we examined side-chain dynamics by measuring the methyl axis S2 values based on deuterium relaxation rates in the 13CH2D moiety of the ILVTMA residues44 (Supplementary Fig. 14 and 15). The dynamically variable side-chains across NS1s (SD-S2axis > Q3) clustered closely in the hydrophobic core; in contrast, none of the core interface residues were dynamically variable across NS1s (Fig. 7d, e). Moreover, dynamically variable side-chains are near the residues whose backbone is dynamically variable among NS1s (Fig. 7f and Supplementary Fig. 12). For some residues, both backbone and side-chain were dynamically variable. These dynamically variable residues directly contact each other (Fig. 7f).
These observations support our hypothesis that a high packing density in the hydrophobic core enables the long-range relay of mutational effects. For example, residue 171 is the most frequently mutated position in NS1 of all human IAVs and directly contacts with residues whose side-chain or backbone showed dynamic variation during NS1 evolution (Fig. 7g). This example illustrates how strain-specific mutational effects are transmitted to the rim interface area. Although core residues make a significant energetic contribution to binding, rim residues also play a role in modulating the local environment around core residues45,46,47,48. Subsequently, it is likely that altered dynamics of the rim interface affect the sampling of productive conformations during the association process to p85β, resulting in strain-dependent epistasis. From the energy landscape perspective, strain-specific mutations altered the conformational energy landscape of NS1s by rearranging the packing interactions, as evidenced by the clustering of dynamically variable residues.
These results demonstrate that near-neutral mutations can silently alter the energy landscape without noticeable changes in functional phenotypes. We speculate that the cumulative effects of near-neutral mutations on the molecular energy landscape help proteins explore new functions during evolution. Therefore, understanding how mutations reshape the biophysical landscape can offer a new mechanistic basis for protein evolution and function14,49.
The evolutionary inference of our results should consider the following caveats. The strain-specific mutations in NS1s employed in the present study might not be neutral with respect to in vivo viral fitness. For example, although the interaction with p85β is highly conserved, NS1s showed some differences in their interactions with other host factors24,50. Thus, we cannot exclude that the mutations in NS1s could affect interactions with other host factors. Here, we cautiously define the “neutral mutation” as limited to the association kinetics to p85β, not concerning other cellular or in vivo functions.
In the present study, we found long-range, strain-specific epistasis in the NS1s of human IAVs and sought to determine its mechanistic basis. Our structural study revealed that strain-specific mutations have limited impacts on the conformation of core p85β-binding residues. In contrast, the hydrophobic core and rim interface showed strain-specific conformational variations among free NS1s. Further studies revealed the hidden network of residues showing varied conformational dynamics during NS1 evolution. These networked residues structurally connect the mutational sites and the rim interface through the hydrophobic core. Therefore, our findings illustrate the high-resolution mechanism whereby a network of closely packed residues underlies long-range epistasis.
Recent studies showed that epistasis is a pervasive phenomenon in protein evolution14,51. A tightly packed hydrophobic core is a common structural feature of well-folded proteins. We speculate that mutations altering the packing interactions might contribute to intramolecular epistatic interactions. Thus, our study provides mechanistic insights into why epistasis prevails. This result is also related to coupled conformational dynamics in the evolution of enzyme catalysis and protein allostery14,15,34,41,49,52,53.
It is noteworthy that the network view of dynamically variable residues only emerges by the statistical measure of fluctuations in conformational (RMSD-CS) and dynamical (SD-S2) data across multiple NS1s (Fig. 6). This result indicates that the conformational or dynamic responses of individual residues in the network differ depending on the NS1 background. Indeed, this statistical picture aligns with the fluctuating pattern of long-range epistasis during NS1 evolution (Supplementary Figs. 4 and 5); strain-specific mutations alter the pattern of epistasis because one mutation remodels the packing interaction differently from another.
Our study also indicated that the molecular energy landscape of NS1s responds more sensitively to mutations during evolution than the functional landscape. For example, strain-specific mutations drastically alter the thermodynamic signature for binding to p85β in addition to affecting its structure and dynamics. Despite drastic changes in the molecular properties, the functional phenotype (i.e., association to p85β) did not change accordingly. This result is reminiscent of thermodynamic system drift, indicating that diverse biophysical mechanisms can result in the same phenotype during protein evolution33. Because of the diverse hosts and high zoonotic potential of IAVs, the NS1s employed in this study are the outcome of diverse evolutionary adaptations54. We conjecture that strain-specific mutations allowed NS1s to explore different evolutionary pathways, resulting in various molecular properties, while the key functionality is conserved. Our results highlight the need for integrated biophysical approaches to reveal the hidden molecular foundations of protein evolution.
Genes encoding all NS1 EDs and p85β proteins were prepared by gene-synthesis service from Genscript: 1918 NS1 (residues 80–205), PR8 NS1 (residues 80–206), Ud NS1 (residues 83–205), VN NS1 (residues 79–206), and p85β (residues 435–599). All NS1 proteins contain W187A mutation to prevent protein dimerization and aggregation55. The incorporation of W187A mutation does not interfere with the binding to p85β23. All proteins were expressed with the N-terminal His6 and SUMO tags in BL21 (DE3) E. coli cells and purified by Ni2+ NTA column chromatography. Tags were removed by SUMO protease and further purified by additional Ni2+ NTA column and gel-filtration chromatography. All purified proteins were >95% pure; the purity of protein samples was assessed by sodium dodecyl sulfate polyacrylamide gel electrophoresis.
Isotopically labeled NS1 proteins for NMR studies were prepared by growing BL21 (DE3) E. coli cells in an M9 medium containing 15NH4Cl and 13C6 glucose as the sole nitrogen and carbon sources. For deuterium relaxation measurements, NS1 proteins were grown in an M9 medium containing 50% D2O/50% H2O with 15NH4Cl and 13C6 glucose. All NMR samples were prepared in a buffer containing 20 mM sodium phosphate (pH 7), 80 mM NaCl, 1 mM TCEP, 1 mM EDTA, and 10% D2O for PR8, Ud, and VN NS1s. NMR sample for 1918 NS1 was prepared in a buffer containing 20 mM sodium phosphate (pH 7), 80 mM NaCl, 1 mM TCEP, 1 mM EDTA, 5 mM ATP, 5 mM MgCl2, and 10% D2O. NS1 concentration in NMR samples was 150–200 μM.
The binding of surface-immobilized NS1 to p85β was measured at 25 °C using an Octet RED biolayer interferometry (Pall ForteBio). The N-terminal His6 and SUMO-tagged NS1 proteins were used for immobilization. The buffer was 20 mM sodium phosphate (pH 7), 150 mM NaCl, 1% bovine serum albumin, and 0.6 M sucrose56. All reported values are the average and standard deviation of three repeated measurements. Association and dissociation data were fit using single exponential growth and decay functions, respectively, using GraphPad Prism (ver. 9). KD values were calculated using measured kon and koff values (KD = koff/kon). For core interface mutants, only kon values were calculated by fitting a linear equation to the plot of kobs vs. [p85].
NMR experiments were conducted on Bruker 800- and 600-MHz spectrometers equipped with a cryogenic probe at the Biomolecular NMR facility (Texas A&M University). All NMR data were collected at 25 °C. NMR data were processed using NMRPipe57 and NMRFAM-SPARKY58. Sequential backbone assignments were performed using 3D triple resonance experiments, including HNCO, HN(CA)CO, HNCACB, CBCA(CO)NH, and HBHA(CO)NH. Side-chain methyl resonance assignments were performed using HCCH-TOCSY, H(CCO)NH, and C(CO)NH experiments. 13C-edited NOESY HSQC was performed with 100 ms mixing time. Chemical shifts for backbone atoms were deposited in the Biological Magnetic Resonance Bank under accession code 51403 for VN NS1 and 51404 for PR8 NS1.
The RMSD-CS values were calculated using the following equation.
where ωx represents the weighting factor accounting for the differences in the gyromagnetic ratios of different types of nuclei. CS represents a chemical shift. i represents different strains (e.g., 1918, PR8, UD, and VN), and N represents the number of NS1s included in the calculation.
Constant relaxation time 15N CPMG RD36,37,38,59 data were recorded at 25 °C on Bruker 600- and 800-MHz NMR spectrometers with CPMG frequencies, ranging from 50 to 1000 Hz, as pseudo-three-dimensional experiments. The uncertainty of R2,eff was estimated by duplicated experiments. The relaxation dispersion data were fit using the Carver-Richards equation60 for a two-state exchange model. Residues undergoing a slow-exchange process were identified based on the Akaike information criterion.
Three relaxation parameters were measured for backbone amides of NS1s at 25 °C on Bruker 600 and 800 MHz NMR spectrometers. For 15N R1 and R2 constants, five relaxation time points were taken with a recycle delay of 2 s, as previously described. For the heteronuclear nuclear Overhauser effect measurements, a recycle delay of 10 s was used in the reference experiment. The steady-state saturation of protons was performed by applying 180° pulses for 4 s61. 2H R1 and R1ρ relaxation parameters were measured at 25 °C on a Bruker 600 MHz NMR spectrometer; pulse sequences provided by Dr. José A. Caro were used after minor modifications. For 2H R1, six relaxation time points were taken between 0.05 and 50 ms. For 2H R1ρ, six relaxation time points were taken between 0.2 and 20 ms. Uncertainties of the relaxation parameters were estimated using duplicated measurements.
The Lipari-Szabo model-free formalism was used to calculate order parameters42. Overall correlation times and rotational diffusion tensors for the axially symmetric model were estimated using the program Quadric62. Backbone order parameters S2NH were calculated using Mathematica (version 12) following the model-selection protocol and the Akaike information criterion63. Side-chain methyl axis order parameters S2axis were calculated using the following equation64.
where S2 is an order parameter for the methyl group, τm is the overall correlation time, and τi−1 = τm−1 + τe,i−1 with τe,i is the effective correlation time of the internal motions for the 13C-2H bond vector. S2axis was calculated by S2 = 0.111S2axis, assuming tetrahedral geometry for the methyl group64.
S2axis values of Leu 13Cδ, whose chemical shifts are larger than 24.5 ppm, were eliminated to prevent artifacts arising from a strong 13C-13C coupling on S2axis65.
Buffer containing 150 μM 15N labeled NS1 was exchanged against 20 mM sodium phosphate (pD 7), 80 mM NaCl, 1 mM TCEP, and 1 mM EDTA in 99.9% deuterium oxide using a spin desalting column (Micro Bio-Spin 6, BioRad). 15N-HSQC spectra were measured every 40 min for 106 h at 25 °C on Bruker 600 MHz NMR spectrometer. H/D exchange rates of backbone amide resonances of NS1 were calculated from a single exponential fit. Protection factors were calculated from the H/D exchange rate divided by the predicted intrinsic exchange rate from the program SPHERE66.
The VN NS1 W187A:p85β complex was crystallized at 277 K by sitting-drop vapor diffusion in 20 mM Tris (pH 7.0) and 80 mM NaCl. The crystal was flash-frozen in liquid nitrogen in the reservoir solution containing 25% (v/v) glycerol. X-ray diffraction datasets were collected at 120 K using an R-AXIS IV++ image plate detector mounted on a Rigaku MicroMax 007HF X-ray generator. The data were processed using iMosflm in the CCP4 package67. The structure was remodeled and refined with Coot68 and the Phenix package (Supplementary Table 1). The crystal structure of VN NS1 W187A in a complex with p85β is deposited in the PDB (ID code: 7RCH). Protein structures were visualized using Pymol and ChimeraX 1.469.
All ITC samples were prepared in 20 mM sodium phosphate (pH 7), 100 mM NaCl, and 2 mM TCEP. Data were recorded at 25 °C using a Microcal-PEAQ-ITC calorimeter (Malvern Panalytical). 100 μM p85β was in the syringe, and 10 μM NS1 was in the cell. 19 consecutive 2 μl aliquots of p85β were titrated into the cell. Data were fit using a 1:1 binding model to obtain Ka, ΔH, and ΔS using software provided with the instrument. All reported parameters are the average and standard deviation of two repeated measurements.
Protein melting temperatures (Tm) were measured using thermal shift analysis. 0.1 μM NS1 was mixed with Sypro Orange (S5692, Sigma-Aldrich) in 20 mM sodium phosphate (pH 7), 100 mM NaCl, and 2 mM TCEP. Data were recorded from 20 °C to 80 °C using a Quantamaster 400 fluorometer (Photon Technology International). Tm was calculated by fitting data to Boltzman sigmoid equation (GraphPad Prism).
The NS1 sequences were aligned using Clustal Omega70 and the sequence logo was generated using WebLogo71.
Further information on research design is available in the Nature Research Reporting Summary linked to this article.
All data generated in this study are available within the Article and Supplementary information. The coordinate of VN NS1 in complex with p85β are available in the Protein Data Bank (PDB) under accession code 7RCH. NMR chemical shifts for backbone atoms are available in the Biological Magnetic Resonance Bank under accession code 51403 for VN NS1 and 51404 for PR8 NS1. Source data are provided with this paper.
Fowler, D. M. & Fields, S. Deep mutational scanning: a new style of protein science. Nat. Methods 11, 801–807 (2014).
Article PubMed PubMed Central CAS Google Scholar
Sailer, Z. R. & Harms, M. J. Molecular ensembles make evolution unpredictable. Proc. Natl Acad. Sci. USA 114, 11938–11943 (2017).
Article ADS PubMed PubMed Central CAS Google Scholar
Domingo, J., Baeza-Centurion, P. & Lehner, B. The causes and consequences of genetic interactions (Epistasis). Annu. Rev. Genomics Hum. Genet. 20, 433–460 (2019).
Article PubMed CAS Google Scholar
Starr, T. N. & Thornton, J. W. Epistasis in protein evolution. Protein Sci. 25, 1204–1218 (2016).
Article PubMed PubMed Central CAS Google Scholar
Zheng, J., Payne, J. L. & Wagner, A. Cryptic genetic variation accelerates evolution by opening access to diverse adaptive peaks. Science 365, 347–353 (2019).
Article ADS PubMed CAS Google Scholar
Baier, F. et al. Cryptic genetic variation shapes the adaptive evolutionary potential of enzymes. Elife 8, e40789 (2019).
Article PubMed PubMed Central Google Scholar
Raman, A. S., White, K. I. & Ranganathan, R. Origins of allostery and evolvability in proteins: a case study. Cell 166, 468–480 (2016).
Article PubMed CAS Google Scholar
Bloom, J. D., Gong, L. I. & Baltimore, D. Permissive secondary mutations enable the evolution of influenza oseltamivir resistance. Science 328, 1272–1275 (2010).
Article ADS PubMed PubMed Central CAS Google Scholar
Ortlund, E. A., Bridgham, J. T., Redinbo, M. R. & Thornton, J. W. Crystal structure of an ancient protein: evolution by conformational epistasis. Science 317, 1544–1548 (2007).
Article ADS PubMed PubMed Central CAS Google Scholar
Harms, M. J. & Thornton, J. W. Evolutionary biochemistry: revealing the historical and physical causes of protein properties. Nat. Rev. Genet. 14, 559–571 (2013).
Article PubMed PubMed Central CAS Google Scholar
Miton, C. M., Buda, K. & Tokuriki, N. Epistasis and intramolecular networks in protein evolution. Curr. Opin. Struct. Biol. 69, 160–168 (2021).
Article PubMed CAS Google Scholar
Nelson, E. D. & Grishin, N. V. Long-range epistasis mediated by structural change in a model of ligand binding proteins. PLoS ONE 11, e0166739 (2016).
Article PubMed PubMed Central CAS Google Scholar
Yu, H. & Dalby, P. A. Coupled molecular dynamics mediate long- and short-range epistasis between mutations that affect stability and aggregation kinetics. Proc. Natl Acad. Sci. USA 115, E11043–E11052 (2018).
Article ADS PubMed PubMed Central CAS Google Scholar
Acevedo-Rocha, C. G. et al. Pervasive cooperative mutational effects on multiple catalytic enzyme traits emerge via long-range conformational dynamics. Nat. Commun. 12, 1621 (2021).
Article ADS PubMed PubMed Central CAS Google Scholar
Campitelli, P., Modi, T., Kumar, S. & Ozkan, S. B. The role of conformational dynamics and allostery in modulating protein evolution. Annu Rev. Biophys. 49, 267–288 (2020).
Article PubMed CAS Google Scholar
Carr, J., Ives, J. & Roberts, N. Virological assessment in vitro and in vivo of an influenza H1N1 virus with a H274Y mutation in the neuraminidase [abstract 79]. In: Program and abstracts of the 13th International Conference on Antiviral Research (Baltimore)) (2000).
Carrillo, B. et al. The influenza A virus protein NS1 displays structural polymorphism. J. Virol. 88, 4113–4122 (2014).
Article PubMed PubMed Central CAS Google Scholar
Ayllon, J., García-Sastre, A. & Hale, B. G. Influenza A viruses and PI3K: are there time, place and manner restrictions? Virulence 3, 411–414 (2012).
Article PubMed PubMed Central Google Scholar
Obenauer, J. C. et al. Large-scale sequence analysis of avian influenza isolates. Science 311, 1576–1580 (2006).
Article ADS PubMed CAS Google Scholar
Heaton, N. S., Sachs, D., Chen, C. J., Hai, R. & Palese, P. Genome-wide mutagenesis of influenza virus reveals unique plasticity of the hemagglutinin and NS1 proteins. Proc. Natl Acad. Sci. USA 110, 20248–20253 (2013).
Article ADS PubMed PubMed Central CAS Google Scholar
Evseev, D. & Magor, K. E. Molecular evolution of the influenza A virus non-structural protein 1 in interspecies transmission and adaptation. Front. Microbiol. 12, 693204 (2021).
Article PubMed PubMed Central Google Scholar
Selman, M., Dankar, S. K., Forbes, N. E., Jia, J. J. & Brown, E. G. Adaptive mutation in influenza A virus non-structural gene is linked to host switching and induces a novel protein by alternative splicing. Emerg. Microbes Infect. 1, 1–10 (2012).
Article Google Scholar
Cho, J. H. et al. Molecular recognition of a host protein by NS1 of pandemic and seasonal influenza A viruses. Proc. Natl Acad. Sci. USA 117, 6550–6558 (2020).
Article ADS PubMed PubMed Central CAS Google Scholar
Jureka, A. S., Kleinpeter, A. B., Cornilescu, G., Cornilescu, C. C. & Petit, C. M. Structural basis for a novel interaction between the NS1 protein derived from the 1918 influenza virus and RIG-I. Structure 23, 2001–2010 (2015).
Article PubMed PubMed Central CAS Google Scholar
Koliopoulos, M. G. et al. Molecular mechanism of influenza A NS1-mediated TRIM25 recognition and inhibition. Nat. Commun. 9, 1820 (2018).
Article ADS PubMed PubMed Central CAS Google Scholar
Hale, B. G. & Randall, R. E. PI3K signalling during influenza A virus infections. Biochem. Soc. Trans. 35, 186–187 (2007).
Article PubMed CAS Google Scholar
Sikosek, T. & Chan, H. S. Biophysics of protein evolution and evolutionary protein biophysics. J. R. Soc. Interface 11, 20140419 (2014).
Article PubMed PubMed Central CAS Google Scholar
Echave, J. Beyond stability constraints: a biophysical model of enzyme evolution with selection on stability and activity. Mol. Biol. Evol. 36, 613–620 (2019).
Article PubMed CAS Google Scholar
Dubrow, A., Kim, I., Topo, E. & Cho, J. H. Understanding the binding transition state after the conformational selection step: the second half of the molecular recognition process between NS1 of the 1918 influenza virus and host p85β. Front. Mol. Biosci. 8, 716477 (2021).
Article PubMed PubMed Central CAS Google Scholar
Cunningham, B. C. & Wells, J. A. High-resolution epitope mapping of hGH-receptor interactions by alanine-scanning mutagenesis. Science 244, 1081–1085 (1989).
Article ADS PubMed CAS Google Scholar
Manhart, M. & Morozov, A. V. Protein folding and binding can emerge as evolutionary spandrels through structural coupling. Proc. Natl Acad. Sci. USA 112, 1797–1802 (2015).
Article ADS PubMed PubMed Central CAS Google Scholar
Hart, K. M. et al. Thermodynamic system drift in protein evolution. PloS Biol. 12, e1001994 (2014).
Article PubMed PubMed Central CAS Google Scholar
Morgunova, E. et al. Two distinct DNA sequences recognized by transcription factors represent enthalpy and entropy optima. Elife 7, e32963 (2018).
Article PubMed PubMed Central Google Scholar
Pinney, M. M. et al. Parallel molecular mechanisms for enzyme temperature adaptation. Science 371, eaay2784 (2021).
Article PubMed CAS Google Scholar
Hale, B. G. et al. Structural insights into phosphoinositide 3-kinase activation by the influenza A virus NS1 protein. Proc. Natl Acad. Sci. USA 107, 1954–1959 (2010).
Article ADS PubMed PubMed Central Google Scholar
Carr, H. Y. & Purcell, E. M. Effects of diffusion on free precession in nuclear magnetic resonance experiments. Phys. Rev. 94, 630–638 (1954).
Article ADS CAS Google Scholar
Meiboom, S. & Gill, D. Modified spin-echo method for measuring nuclear relaxation times. Rev. Sci. Instrum. 29, 688–691 (1958).
Article ADS CAS Google Scholar
Loria, J. P., Rance, M. & Palmer, A. G. A relaxation-compensated Carr-Purcell-Meiboom-Gill sequence for characterizing chemical exchange by NMR spectroscopy. J. Am. Chem. Soc. 121, 2331–2332 (1999).
Article CAS Google Scholar
Olson, C. A., Wu, N. C. & Sun, R. A comprehensive biophysical description of pairwise epistasis throughout an entire protein domain. Curr. Biol. 24, 2643–2651 (2014).
Article PubMed PubMed Central CAS Google Scholar
Kumar, A. et al. Stability-mediated epistasis restricts accessible mutational pathways in the functional evolution of avian hemoglobin. Mol. Biol. Evol. 34, 1240–1251 (2017).
Article PubMed PubMed Central CAS Google Scholar
González, M. M., Abriata, L. A., Tomatis, P. E. & Vila, A. J. Optimization of conformational dynamics in an epistatic evolutionary trajectory. Mol. Biol. Evol. 33, 1768–1776 (2016).
Article PubMed PubMed Central CAS Google Scholar
Lipari, G. & Szabo, A. Model-free approach to the interpretation of nuclear magnetic resonance relxation in macromolecules. 1. Theory and range of validity. J. Am. Chem. Soc. 104, 4546–4559 (1982).
Article CAS Google Scholar
Trbovic, N., Cho, J. H., Abel, R., Friesner, R. A. & Rance, M. Palmer AGr. Protein side-chain dynamics and residual conformational entropy. J. Am. Chem. Soc. 131, 615–622 (2009).
Article PubMed PubMed Central CAS Google Scholar
Farrow, N. A., Zhang, O., Szabo, A., Torchia, D. A. & Kay, L. E. Spectral density function mapping using 15N relaxation data exclusively. J. Biomol. NMR 6, 153–162 (1995).
Article PubMed CAS Google Scholar
Hendsch, Z. S. & Tidor, B. Do salt bridges stabilize proteins? A continuum electrostatic analysis. Protein Sci. 3, 211–226 (1994).
Article PubMed PubMed Central CAS Google Scholar
Joughin, B. A., Green, D. F. & Tidor, B. Action-at-a-distance interactions enhance protein binding affinity. Protein Sci. 14, 1363–1369 (2005).
Article PubMed PubMed Central CAS Google Scholar
Kiel, C., Selzer, T., Shaul, Y., Schreiber, G. & Herrmann, C. Electrostatically optimized Ras-binding Ral guanine dissociation stimulator mutants increase the rate of association by stabilizing the encounter complex. Proc. Natl Acad. Sci. USA 101, 9223–9228 (2004).
Article ADS PubMed PubMed Central CAS Google Scholar
Pace, C. N. et al. Contribution of hydrogen bonds to protein stability. Protein Sci. 23, 652–661 (2014).
Article PubMed PubMed Central CAS Google Scholar
Otten, R. et al. How directed evolution reshapes the energy landscape in an enzyme to boost catalysis. Science 370, 1442–1446 (2020).
Article ADS PubMed CAS Google Scholar
Kuo, R. L. et al. Role of N terminus-truncated NS1 proteins of influenza A virus in inhibiting IRF3 activation. J. Virol. 90, 4696–4705 (2016).
Article PubMed PubMed Central CAS Google Scholar
Breen, M. S., Kemena, C., Vlasov, P. K., Notredame, C. & Kondrashov, F. A. Epistasis as the primary factor in molecular evolution. Nature 490, 535–538 (2012).
Article ADS PubMed CAS Google Scholar
Otten, R. et al. Rescue of conformational dynamics in enzyme catalysis by directed evolution. Nat. Commun. 9, 1314 (2018).
Article ADS PubMed PubMed Central CAS Google Scholar
Hajredini, F., Piserchio, A. & Ghose, R. Long-range dynamic correlations regulate the catalytic activity of the bacterial tyrosine kinase Wzc. Sci. Adv. 6, eabd3718 (2020).
Article ADS PubMed CAS Google Scholar
Taubenberger, J. K. & Kash, J. C. Influenza virus evolution, host adaptation, and pandemic formation. Cell Host Microbe 7, 440–451 (2010).
Article PubMed PubMed Central CAS Google Scholar
Aramini, J. M. et al. Dimer interface of the effector domain of non-structural protein 1 from influenza A virus: an interface with multiple functions. J. Biol. Chem. 286, 26050–26060 (2011).
Article PubMed PubMed Central Google Scholar
Dubrow, A., Zuniga, B., Topo, E. & Cho, J. H. Suppressing nonspecific binding in biolayer interferometry experiments for weak ligand-analyte interactions. ACS Omega 7, 9206–9211 (2022).
Article PubMed PubMed Central CAS Google Scholar
Delaglio, F. et al. NMRPipe: a multidimensional spectral processing system based on UNIX pipes. J. Biomol. NMR 6, 277–293 (1995).
Article PubMed CAS Google Scholar
Lee, W., Tonelli, M. & Markley, J. L. NMRFAM-SPARKY: enhanced software for biomolecular NMR spectroscopy. Bioinformatics 31, 1325–1327 (2015).
Article PubMed Google Scholar
Mulder, F. A., Skrynnikov, N. R., Hon, B., Dahlquist, F. W. & Kay, L. E. Measurement of slow (micros-ms) time scale dynamics in protein side chains by (15)N relaxation dispersion NMR spectroscopy: application to Asn and Gln residues in a cavity mutant of T4 lysozyme. J. Am. Chem. Soc. 123, 967–975 (2001).
Article PubMed CAS Google Scholar
Carver, J. P. & Richards, R. E. A general state-site solution for the chemical exchange produced dependence of T2 upon the Carr-Purcell pulse separation. J. Magn. Reson. 6, 89–105 (1972).
ADS CAS Google Scholar
Ferrage, F., Reichel, A., Battacharya, S., Cowburn, D. & Ghose, R. On the measurement of 15N-{1H} nuclear Overhauser effects. 2. Effects of the saturation scheme and water signal suppression. J. Magn. Reson 207, 294–303 (2010).
Article ADS PubMed PubMed Central CAS Google Scholar
Lee, L. K., Rance, M., Chazin, W. J. & Palmer, A. G. Rotational diffusion anisotropy of proteins from simultaneous analysis of 15N and 13C alpha nuclear spin relaxation. J. Biomol. NMR 9, 287–298 (1997).
Article PubMed CAS Google Scholar
Mandel, A. M., Akke, M. & Palmer, A. G. Backbone dynamics of Escherichia coli ribonulecease H1: Correlations with structure and function in an active enzyme. J. Mol. Biol. 246, 144–163 (1995).
Article PubMed CAS Google Scholar
Muhandiram, D. R., Yamazaki, T., Sykes, B. D. & Kay, L. E. Measurement of 2H T1 and T1.rho. relaxation times in uniformly 13C-labeled and fractionally 2H-labeled proteins in solution. J. Am. Chem. Soc. 117, 11536–11544 (1995).
Article CAS Google Scholar
Hoffmann, F., Xue, M., Schäfer, L. V. & Mulder, F. A. A. Narrowing the gap between experimental and computational determination of methyl group dynamics in proteins. Phys. Chem. Chem. Phys. 20, 24577–24590 (2018).
Article PubMed CAS Google Scholar
Bai, Y., Milne, J. S., Mayne, L. & Englander, S. W. Primary structure effects on peptide group hydrogen exchange. Proteins 17, 75–86 (1993).
Article PubMed PubMed Central CAS Google Scholar
Adams, P. D. et al. PHENIX: a comprehensive Python-based system for macromolecular structure solution. Acta Crystallogr. D. Biol. Crystallogr. 66, 213–221 (2010).
Article PubMed PubMed Central CAS Google Scholar
Emsley, P., Lohkamp, B., Scott, W. G. & Cowtan, K. Features and development of Coot. Acta Crystallogr. D. Biol. Crystallogr. 66, 486–501 (2010).
Article PubMed PubMed Central CAS Google Scholar
Pettersen, E. F. et al. UCSF ChimeraX: structure visualization for researchers, educators, and developers. Protein Sci. 30, 70–82 (2021).
Article PubMed CAS Google Scholar
Madeira, F. et al. Search and sequence analysis tools services from EMBL-EBI in 2022. Nucleic Acids Res. 50, W276–W279 (2022).
Article PubMed Central Google Scholar
Crooks, G. E., Hon, G., Chandonia, J. M. & Brenner, S. E. WebLogo: a sequence logo generator. Genome Res. 14, 1188–1190 (2004).
Article PubMed PubMed Central CAS Google Scholar
Download references
We thank Prof. Wonmuk Hwang, Prof. Tatyana Igumenova, and James Gonzales for carefully reading the manuscript and providing constructive comments. Support from NIH grant R01GM127723 (J.H.C.), the Welch Foundation A-2028-20200401 (J.H.C.), and USDA National Institute of Food and Agriculture grant Hatch project 1020344 (J.H.C.) is acknowledged.
Department of Biochemistry and Biophysics, Texas A&M University, College Station, TX, 77843, USA
Iktae Kim, Alyssa Dubrow, Bryan Zuniga, Baoyu Zhao, Noah Sherer, Abhishek Bastiray, Pingwei Li & Jae-Hyun Cho
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
J.H.C. conceived the project. I.K. conducted NMR resonance assignments and ITC measurements. A.D. and B.Zu. acquired BLI data. B.Zh., I.K., P.L., and J.H.C. conducted protein crystallization and structure determination. N.S. and A.B. conducted co-IP.
Correspondence to Jae-Hyun Cho.
The authors declare no competing interests.
Nature Communications thanks the anonymous reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
Reprints and Permissions
Kim, I., Dubrow, A., Zuniga, B. et al. Energy landscape reshaped by strain-specific mutations underlies epistasis in NS1 evolution of influenza A virus. Nat Commun 13, 5775 (2022). https://doi.org/10.1038/s41467-022-33554-9
Download citation
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41467-022-33554-9
Anyone you share the following link with will be able to read this content:
Sorry, a shareable link is not currently available for this article.
Provided by the Springer Nature SharedIt content-sharing initiative
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.
Advertisement
© 2023 Springer Nature Limited
Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.