Molecular dynamics simulations of the Nip7 proteins from the marine deep- and shallow-water Pyrococcus species
BMC Structural Biology volume 14, Article number: 23 (2014)
The identification of the mechanisms of adaptation of protein structures to extreme environmental conditions is a challenging task of structural biology. We performed molecular dynamics (MD) simulations of the Nip7 protein involved in RNA processing from the shallow-water (P. furiosus) and the deep-water (P. abyssi) marine hyperthermophylic archaea at different temperatures (300 and 373 K) and pressures (0.1, 50 and 100 MPa). The aim was to disclose similarities and differences between the deep- and shallow-sea protein models at different temperatures and pressures.
The current results demonstrate that the 3D models of the two proteins at all the examined values of pressures and temperatures are compact, stable and similar to the known crystal structure of the P. abyssi Nip7. The structural deviations and fluctuations in the polypeptide chain during the MD simulations were the most pronounced in the loop regions, their magnitude being larger for the C-terminal domain in both proteins. A number of highly mobile segments the protein globule presumably involved in protein-protein interactions were identified. Regions of the polypeptide chain with significant difference in conformational dynamics between the deep- and shallow-water proteins were identified.
The results of our analysis demonstrated that in the examined ranges of temperatures and pressures, increase in temperature has a stronger effect on change in the dynamic properties of the protein globule than the increase in pressure. The conformational changes of both the deep- and shallow-sea protein models under increasing temperature and pressure are non-uniform. Our current results indicate that amino acid substitutions between shallow- and deep-water proteins only slightly affect overall stability of two proteins. Rather, they may affect the interactions of the Nip7 protein with its protein or RNA partners.
High temperatures and pressures cause damage to living cells. For humans and the best studied organisms, conditions with temperature close to 27°C (300 K) and atmospheric pressure around 0.1 MPa are optimal. However, there exist organisms, which colonize habitats extreme, life-incompatible for humans. Such conditions are near deep hot springs colonized by communities of organisms, the extremophiles . Their life is sustained under conditions with temperature as high as 100°C (373 K) and pressure above 20 MPa, exceeding the atmospheric by 200 times. The mechanisms by which cell survival is provided are not well understood. Their elucidation would provide answers to some fundamental questions on the origins of life and the early adaptation of microorganisms , also on adaptation to the conditions of diverse ecological niches . A timely challenge is the identification of the molecular mechanisms of the evolutionary adaptation of the genomes and proteomes of the living beings to the conditions of high temperature - and pressure -. Extremophiles were supposed to provide unprecedented opportunities for biotechnological explorations (single enzyme catalysis) -. Indeed, based on research on extremophiles, enzymes were developed for biotechnological applications.
To study the possible mechanisms of the influence of high temperature and pressure on the protein dynamics and protein adaptation to altered environmental pressure, we used here computer MD simulations of two homologous Nip7 proteins from hyperthermophilic (optimal growth temperature close to 100°C) archaea, shallow water P. furiosus (hydrostatic pressure close to atmospheric)  and deep-water, P. abyssi (hydrostatic pressure close to ~20 MPa) ,. The goal of our work was to compare the dynamics properties of the models we built at high and low temperatures, also at atmospheric and high pressures to identify their common characteristics, also their differences presumably related to the different depths of the organisms’ habitats.
Nip7 was initially identified in yeast as required for processing of the 27S pre-rRNA to form the mature 25S and 5.8S rRNAs . It localizes to the nucleolus but was also found to sediment in the region of free 60S subunits in sucrose density gradients , which is consistent with its presence in pre-60S complexes . Experimental evidence suggests that the P. abyssi Nip7 may be an exosome regulatory factor. It binds preferentially to U- and AU-rich RNAs and strongly inhibits the exosome due to its association with both the exosome complex and the substrate RNA .
The 3D structure of P. abyssi Nip7 protein is known (PDB ID 2P38; Figure 1A). The protein polypeptide chain is 166 amino acids long of which 155 residues are represented in the 3D structure. The protein consists of two α-β domains , Figure 1. The N-terminal domain (residues 1-90) is composed of five antiparallel β-strands surrounded by three α-helices and one 310 helix. There is an assumption that archaeal Nip7 may interact with exosome via its N-terminal domain, thereby controlling the exosome function . However, the molecular mechanism of this interaction is unknown.
The C-terminal domain is assigned to the PUA class. . It includes amino acid residues 91-155 and is comprised of a mixture of β-sheets, one α-helix, and one short 310 helix . This domain, named after pseudouridine synthases and archaeosine-specific transglycosylases , was initially described in tRNA modifying enzymes and in pseudo-uridine synthases from Archaea and eukaryotes , and it has been proposed to mediate protein-RNA interactions. Comparative structural analyses revealed that the residues involved in RNA contacts are conservative in the archaeal PUA domains . In the P. abyssi Nip7, these are residues K103, L107, D113, P115, E117, R151, R152, K155, L157, K158 (Figure 1), whose interactions with RNA have been confirmed by the results of mutation experiments . The PUA domain contacts the RNA molecule using a glycine-containing loop, which connects the fifth α-helix and the β9 strand via residues of the β12 strand (Figure 1). The structural alignment suggests that Nip7 may use a mechanism similar to ArcTGT for RNA interaction, which binds to the bottom of the tRNAVal acceptor stem through the major groove. In the yeast and human Nip7, however, some RNA-binding residues are replaced by a glycine and three residues with hydrophobic side chain .
The Nip7 from P. abyssi and its P. furiosus homolog share 70% identity in their sequences (Figure 1B). Comparative analysis identified the excess of radical versus conservative amino acid substitutions fixation rates in Nip7 P. furiosus after divergence from the deep-sea ancestor it shared with P. abyssi and P. horikoshii. It was suggested that Nip7 and some other proteins, which were concerned predominantly with “translation machinery“ and “ribosomal function“ evolved under positive Darwinian selection, resulted from their adaptation to altered conditions of elusive pressures. However, the molecular mechanisms underlying this mode of gene evolution remained unclear.
A deletion/insertion free alignment is an advantage for simulations of the P. furiosus 3D structure allowing to omit the reconstruction of loops. The Nip7 proteins are simulated at atmospheric pressures, also at 50 and 100 MPa. Given the fact that the temperatures optimal for these two organisms are close to 373 K, we consider simulations at room (300 K) and elevated (373 K) temperatures.
The results show that the models of the two proteins at all the examined pressures and temperatures are compact, stable and similar to the crystal Nip7 structure. Analysis of protein structure dynamics at different pressures and temperatures allowed us to disclose similarities and differences between the deep- and shallow-sea protein models at different temperatures and pressures.
Stability of models
Figure 2 shows changes in the root mean - square deviations (the RMSDs) of the protein Cα atoms relative to the starting structures during the MD simulations for the trajectories at 50 MPa. From Figure 2 it follows that the structures achieve equilibrium starting from about 10 ns of the trajectory periods. The deviations from the starting structure are about 2 Å for the NIP7-ABY model, they are about 2.25 Å for the NIP7-FUR model.
We compare the structures during simulations of the two protein models with the crystal structure of the P. abyssi Nip7 (2P38:A). Table 1 gives the mean RMSD values and the 95% confidence intervals at different temperatures and pressures calculated using five MD trajectories for each parameter set.
As shown in Table 1, the protein deviates from the crystal structure at higher temperatures; however, the mean RMSD is not greater than 2.41 Å for NIP7-FUR and 2.0 Å for NIP7-ABY, thereby supporting the inference that the structures are stable throughout the MD simulations at all the parameters under study. Table 1 also shows that the RMSD for both proteins tends to increase with rising temperature.
The radius of gyration (Rg) is another indicator of stability of a structure during simulation. The average Rgs of the protein structures for trajectories at different pressures and temperatures are compared in Figure 3. As shown in Figure 3, most differences between the Rg values for different trajectories do not exceed the standard errors of means, thereby providing evidence that our models are stable. It should be noted that this is associated with somewhat greater Rg for NIP7-FUR than NIP7-ABY. Two-way ANOVA of the Rg values for NIP7-FUR and NIP7-ABY demonstrated that alterations in pressure and temperature did not result in significant changes in the mean Rg values for both protein models (Additional file 1).
Analysis of solvent accessibility
We compute for each trajectory the average values of solvent accessibility of the residues for the entire protein (SASt), also separately for the polar (SASp) and hydrophobic (SASh) surface fractions. We calculate the mean and standard error of these parameters over five runs for each protein model and pressure-temperature values set. The results are given in Figure 4(A-C).
As seen in Figure 4, the total area is somewhat smaller for NIP7-FUR. Higher values for the solvent accessible polar portion of the residues are observed for the NIP7-ABY than for the NIP7-FUR model (compare the blue/red line values for NIP7-ABY with the brown/green values for NIP7-FUR in Figure 4B). The area of the hydrophobic part of the model is smaller (the same lines in Figure 4C). These observations concerned all the pressure and temperature parameters we studied.
The following trends in the values of the SAS parameters of the NIP7-FUR model are notable: (1) the rise in temperature at the same pressures increases the solvent accessibility of the polar portion of the protein (SASp values for NIP7-FUR at 373 K, brown, are higher than the for NIP7-FUR 300 K, green, with the exception of 100 MPa; Figure 4B) and decreases the area of the hydrophobic portion of the surface, this is associated with a small change in the total area (Figure 4C); (2) the increase in pressure systematically reduces the SAS values (brown and green lines, Figure 4B, C). However, most changes are in the range of the mean standard errors.
As for the NIP7-ABY model, the behavior of its SAS parameters was ambiguous. SASp at atmospheric pressure and 373 K is greater, than at 300 K, whereas at 100 MPa, on the contrary, is lower (compare the red/blue line Y values, Figure 4B). A significant increase in SASh at elevating temperature (red line is above blue one, Figure 4C; 100 MPa was an exception) was other characteristic of this model. At the constant high temperature, increase in pressure resulted in decrease in the SAS values of all the three types (red lines on Figure 4A, B, C). In contrast, at constant T = 300 K, the SASp values increased with increasing pressure, although within the limits of standard deviations (blue line, Figure 4B). Changes of SASh for the NIP7-ABY model demonstrates no systematic trend at increasing pressure (blue line, Figure 4C); as for SASh, it slightly changes at P = 50 MPa, as compared with P = 0.1 MPa, it rises slightly with pressure increasing to 100 MPa, becoming actually equal to SASh at P = 0.1 MPa.
We estimated the influence of temperature and pressure on change in the surface area of the models on the basis of two-way ANOVA (Additional file 2). Taken together, the data agree with the analysis shown in Figure 4. Temperature exerts a significant influence on the hydrophobic part of the NIP7-ABY surface area, thereby demonstrating that the influence of pressure and temperature on change in the total area of protein surface is not additive. This presumably reflects the difference in the SAS values between the high and low temperature NIP7-ABY trajectories. As for the NIP7-FUR model, pressure exerts a significant influence on both the hydrophobic part of the area and SASt.
Thus, pressure and temperature differently affect residue solvent accessibility for the two protein models. Pressure results in a decrease in the solvent accessible surface area. This effect is the most conspicuous at high temperatures.
Local structure of the polypeptide chain
The next step was to define how the conformation of the different regions of the polypeptide chain deviates from the X-ray Nip7 structural model. For different trajectories, we built graphs showing the dependencies of the local structure deviation (RMSDL) parameter (see Methods). These graphs concerned every residue; we estimated also the mean RMSDL error, which characterizes its fluctuation during the simulations. (Additional file 3: Figure S1). The graphs give prominence to the non-uniform local changes in the protein structure during the simulations. There are regions where deviations are quite small (~1 Å), while in the other regions changes in the local chain conformation relative to the crystal structure are considerable (~4 Å). In such cases, the RMSDLs of the polypeptide chain conformation are greater in those regions, which correspond to the loops connecting the secondary structure elements. Regions with great RMSDL values are, as a rule, located where amino acids substitutions are numerous (see Additional file 3: Figure S1A). As a result, the evolutionary and structurally variable regions are virtually the same (with the exception of positions 85-95 making up the segment between the α4-α5 helices, Additional file 3: Figure S1A). For this segment, the RMSDL values are low, whereas the number of differences between the P. abyssi and P. furiosus Nip7 sequences is high.
Other features concern the differences in the RMSDLs values between the N- and C-terminal domains . There are regions with small RMSDL values and peaks in the loop regions (positions 30-40, which correspond to the β3-β4 loop, positions 49-59 and 69-79) for the N-terminal domain. The graphs of structural variations in the C-terminal, the PUA-domain, also contain maxima at positions 105-110, 115-125, 143-148 (Additional file 3: Figure S1) and minima; however, the minima values are substantially higher than those for the N-terminal domain. From comparisons of the deviation graphs for the different trajectories, it is evident that the regions corresponding to the RMSDL maxima are by and large the same for the trajectories that correspond to both high pressures and temperatures (Additional file 3: Figure S2). This means that the local changes in the conformation of the protein chain in our models occur in the same region of the structure despite the different nature of the destabilizing factor (changes in pressure or temperature).
The effects of pressure and temperature on changes in the conformation of the polypeptide chain in the Nip7 models are now considered. The most conspicuous fact is that, at the same pressure, elevated temperature causes an increase in the RMSDLs for some protein residues (compare Additional file 3: Figure S1, panel pairs A-B, C-D, E-F). The changes in the RMSDLs of the polypeptide chain resulting from elevated temperature occur non-uniformly. For example, greater RMSDLs at elevated temperature are observed in the NIP7-ABY model in the region of the α2-α3 loop (at positions 49-59). Conversely, for some regions in this model, an increase in temperature results in their decrease, for example residues 109-127 of the C-terminal domain. A similar pattern is observed for this model at low (Additional file 3: Figure S1, panels A-B), moderate (C-D) and high (E-F) pressures. However, for some regions of the protein changes in the RMSDLs at elevated temperature are only slight.
The interesting features brought out by comparing the two proteins are the higher RMSDL values for the C-terminal domain in NIP7-FUR compared to NIP7-ABY at any pressure-temperature parameter (Additional file 3: Figure S1).
To make the local structure changes more prominent with respect to the temperature and pressure, we performed an analysis of the sensitivity of the residue RMSDL parameter to temperature (St i ) and pressure (Sp i ) expressed as the ratio of the RMSDL at high and low temperature and high and low pressure (see Methods). We also estimated the dependence of the RMSDL parameter for each residue on the temperature and pressure using two-way ANOVA (see Methods). The influence of temperature and pressure on the RMSDL parameter was significant, if the corresponding F-statistics was above the threshold at α = 0.05 (see Additional file 4).
Figure 5A illustrates that temperature sensitivity of the RMSDL parameter is different, being dependent on protein position in the two models. To begin with, there are protein regions for which temperature significantly increases the RMSDL ( >1) for both models. To such regions belongs, in the first place, the long segment extending from the C-terminal part of the β4 strand to the N-terminal end of the β5 strand (positions 41-65). The F-statistics for most of its residues exceeds by many times the 5% critical value for both models (Additional file 4). In the graphs given in Additional file 3: Figure S1, this region corresponds to the peak RMSDL values, which increase when temperature elevates. A small, yet significant, increase in the RMSDL is established in the N-terminal domain for both models at positions 26-27 (the β3 strand of the N-terminal domain). A significant increase in the RMSDL residues at positions 31-37 (β3-β4) under the influence of high temperature is also characteristic of the NIP7-FUR model. High values of the parameter are observed for the N-terminal regions of both models, however, the influence of temperature on change in the RMSDL for these residues proved to be insignificant.
An interesting pattern was observed for the C-terminal domain. <1 for the great majority of residues, the differences were significant for many of them. This meant that elevation in temperature made the conformation of the C-terminal domain to a great extent similar to that of the crystal structures. These changes were the most prominent for region 109-127 (α6-β10). This was associated with lower values characteristic of the NIP7-ABY model. Another region, which was also subject to conformational change at elevated temperature in the models, joined two protein domains (positions 80-93). A significant decrease in the RMSDL was characteristic of the C-terminal residues in both models, too.
The sensitivity of the RMSDL values to elevated pressure is shown in Figure 5B. We compared and with respect to the positions at which their proportions are significantly different for both models. We found that their number is considerably smaller for than . For the NIP7-FUR model, significant changes are observed in the N-terminal domain (positions 30-38). For the NIP7-ABY model, fluctuations in the values are more expressed, however, changes at positions 41-44 (a decrease in the RMSDL) and 85-89 (an increase in the RMSDL) are significant.
Analysis of the secondary structure
To describe in more detail the conformational changes, we analyzed the secondary structure of the Nip7 protein models at different MD trajectories. Additional file 3: Figure S2 provides evidence for the stability of the secondary structure during the simulations. For example, although the RMSDL values for the loop 69-79 (the β5-β6 region) are high (Additional file 3: Figure S1), the secondary structure pattern of the polypeptide chain demonstrates no changes dependent on pressure and temperature. This may suggest that this β-hairpin changes its conformation without breaking hydrogen bonds, just by bending or twisting.
Analysis of the data on the changes in the protein secondary structure (Additional file 3: Figure S2) shows also that a part of its elements is unstable when pressure and temperature change. For example, this is characteristic of the loop formed by residues 48-62 (the α2-α3 segment). At different pressures and temperatures, the conformation of these residues varies from the states of the α-helix (H), the 310 -helix (G), the turn (T), and to the bend (B). A local maximum of the deviations of the polypeptide chain from the crystal structure, the RMSDL (Additional file 3: Figure S1) are observed in the segment. Thus, this protein segment proves to be very mobile.
Other regions of smaller size with unstable secondary structures are the first N-terminal residues, the β10 terminal residues (transitions from the loop to the extended conformation) and the β11-β12 loop (transitions between the turn, the bend and the loop conformations). These regions show characteristic transitions from the loop to the bend, also the maximum RMSDL values (Additional file 3: Figure S1).
The above described conformational changes are supported by the trajectory snapshots providing new dimensions of the credibility of our models (Figure 6; Additional file 3: Figure S3). Thus, residues 69-75, which belong to the β-turn, have a stable secondary structure, but high RMSDL values. These residues bend towards the globular part of the N-terminal domain during the simulations. In so doing, they twist slightly the β-sheet. The sweep of the loop 69-75 bend depends on pressure and temperature. Similar conformational changes (the bend towards the centre of the N-terminal domain without significant impairment of the secondary structure) are also observed for the loop between the β3-β4 strands. It also shows characteristic RMSDL deviations (Additional file 3: Figure S1). As seen in Figure 6B, the range of the fluctuations widens at high temperatures in the NIP7-ABY model, a similar widening is true for the NIP7-FUR model (Figure 6, C, D). The structural fluctuations in the segment 45-55 and those in the loops between the β-strands of the terminal domain are conspicuous, too.
Fluctuations in the polypeptide chain
We estimated the structural flexibility of the two proteins at different temperatures and pressures. The graphs display the dependencies of the Root-Mean-Squared-Fluctuation (the RMSF) on the residue number (Additional file 3: Figure S4).
The peak regions of the RMSF values enclose predominantly loops in the sequence regions at positions 31-41 (the β3-β4 loop), 50-60 (the α3 region), 69-79 (the β5-β6 region), 98-110 (the β8-α5 region), 120-125 (the β9-β10 region), 132-135 (the β10-β11 region), 140-150 (the β11-β12 loop). These regions correspond to those where the RMSFs are the largest. In addition, the highest values for the local root-mean-square deviations of the polypeptide chain (the RMSDLs) are observed in them (compare with Additional file 3: Figure S1).
The other features are the high RMSF values for the residues at the amino (N) and carboxy (C)-termini of the proteins. Notably, the RMSF values are somewhat higher for the trajectories at high temperature than for those at 300 K. These differences were observed for both NIP7-ABY and NIP7-FUR (Additional file 3: Figure S4).
The sensitivity of the RMSF value to increase in temperature and pressure was analyzed for each and every amino acid residue in the same manner as was for the RMSDL parameter. The results are shown in Figure 7 (see Additional file 5 for the two-way ANOVA results). Analysis of the sensitivity profiles of the residues to temperature demonstrated that increase in temperature raised the value of the local polypeptide fluctuations for all the residues at high significance ( >1). However, the graph displays regions weakly sensitive to increase in temperature, even insensitive to it (these are residues 36, 121, 123 and 131 for the NIP7-ABY model and residue 146 for the NIP7-FUR model). The RMSF values are high in the residues of the loop already at 300 K, while increase in temperature does not cause any considerable increase in their mobility.
The results are different for fluctuations in the Nip7 models of the deep- and shallow-sea organisms at increased pressure (Figure 7B). The parameter is above unity for a small fraction of residues. For the N-terminal domain, the decrease in fluctuations under the effect of pressure is greater for the NIP7-ABY than the NIP7-FUR model. In this domain, pressure significantly influences residues located at positions 6-15 (the N-terminal part of the α1 helix), positions 23, 29 (the loop between β2 and β3), and 54-55 (the α2-α3 loop).
For the NIP7-FUR RMSF, a significant influence of pressure is observed in residues mostly in the C-terminal part of the protein. A wide region of these residues spans positions 97-110 (the α5 helix), the most sensitive region to increase in pressure and <1 (decrease in the RMSF). Several other residues for which RMSF decreased by pressure significantly are located the β10 region (positions 125-130), β11 (138-140), β12 (152-153) and at position 89. Several positions demonstrating significant decrease in the RMSF under pressure are located in the N-terminal domain (66, 72, 78-80). There are, however, residues for which the RMSF changes significantly when pressure increases. They are observed for NIP7-FUR only (54, 124, 132). Interestingly, the single position, which demonstrates the significant dependence of its RMSF value on the pressure for both models, 54, is located in the α2-α3 loop. However, the influence is inversed for this residue in the NIP7-FUR and NIP7-ABY models.
To analyze the influence of pressure and temperature on fluctuations in the models, we performed two-way ANOVA at different temperatures for both the entire models and separately for each domain. The mean RMSF values for the Cα atoms in the NIP7-ABY and NIP7-FUR models at different trajectories were compared (see Additional file 6). Statistical analysis demonstrates that, for the NIP7-ABY model at 300 K, change in pressure does not affect significantly fluctuations in the N- and C- domains. However, differences in fluctuations between these domains are significant. In general, the C-terminal domain demonstrates larger fluctuations compared with the N-terminal domain. At high temperature, the influence of pressure on the RMSF of two domains becomes significant. The fluctuations in both domains decreased at high pressure. The difference in the fluctuations between two domains at high temperature remains the same as at 300 K. In other words, the fluctuations are greater for the C-terminal domain.
For the NIP7-FUR model, pressure exerts a significant influence on the RMSF values in the two domains both at low and high temperatures (see Additional file 6). The fluctuations in both domains decrease with increase in pressure, like in the NIP7-ABY model. The differences in the RMSF values for the two domains are significant at 300 and 373 K.
Analysis of salt bridges
The crystal structure of P. abyssi contains eight pairs of residues forming salt bridges . It should be noted that the positions of the side chain atoms have not been resolved for some residues capable of forming salt bridges probably due to their decreased stability under crystallization conditions. These are GLU51, LYS76, for example. In our simulations it is feasible to reconstruct side chain atoms for these residues and to take them into account in salt-bridge stability analysis. For each and every possible salt bridge, the proportion of structures was estimated where such a bridge was formed (salt bridge persistence; see Methods). This was done for all the trajectories and models under study. Analysis of the NIP7-ABY and NIP7-FUR models revealed that there were stable salt bridges (occurring, on average, in 70% of the structures in the different trajectories of the two proteins), moderately stable (20-70% persistence) and a number of unstable salt bridges occurring in less than 20% of the structures (Additional file 7). The number of stable and moderately stable salt bridges was larger in NIP7-ABY, twenty, than in NIP7-FUR, thirteen. Interestingly, the residues with unresolved side chain atoms in the crystal structure formed moderately stable salt bridges only in the NIP7-FUR and NIP7-ABY models.
Six salt bridges (GLU10-ARG4, GLU75-ARG37, ARG37-GLU10, LYS20-GLU17, ARG148-ASP109, GLU33-ARG4) are the most stable in the NIP7-ABY model (Additional file 7). These form in more than 70% of the structures in both the NIP7-ABY and NIP7-FUR models. The GLU33-ARG4 and GLU131-ARG116 pairs (unresolved in the 2P38:A crystal structure) may be referred to the most conservative in NIP7-FUR (being moderately stable in NIP7-ABY). Four stable salt bridges from the N-terminal domain (GLU10-ARG4, GLU75-ARG37, ARG37-GLU10, GLU33-ARG4) form a network that links the alpha-helix and the beta-strands. The remaining salt bridges may be assigned to the moderately stable. Some stable or moderately stable salt bridges in NIP7-ABY are unstable in the NIP7-FUR (148-109, 147-113). There are also two salt bridges moderately stable in the NIP7-FUR and unstable in the NIP7-ABY models (146-112, 76-17) (Table 2).
Analysis of the functional role of the most stable salt bridges demonstrated (Table 2) that most of them stabilized the packaging of the elements of the protein secondary structure. However, a number of salt bridges in our models formed residues for which interactions with RNA has been demonstrated . There are four such pairs among stable bridges for the NIP7-ABY model and only one for NIP7-FUR model (148-109) is formed by a residue pair with both involved in the interaction with RNA, the other salt bridges are either unstable or missing in the NIP7-FUR model.
How high temperatures and pressures may affect the formation of salt bridges in the two different protein models? To answer this question, we introduced for each salt bridge its index of sensitivity to pressure (Sp) and temperature (St) in the same manner as we did for the RMSDL and the RMSF parameters (Methods). The higher were the values of these indices, the greater was the persistence of a bridge at high pressure or temperature. To estimate the significance of the influence of pressure and temperature on the stability of a salt bridge, we applied two-way ANOVA and obtained values of the statistics F(T), F(P) (see Methods). We assumed that high temperature or pressure exerted an influence on stability of a salt bridge, if the corresponding F(T) and F(P) statistics exceeded the 0.05 significance threshold. The results for the stable and moderately stable salt bridges common to both models are given in Table 2, Additional files 8 and 9. They demonstrate that for the NIP7-ABY model, of the 20 stable salt bridges 3 (15%) are affected by high temperatures (131-116, 76-70 and 88-23). One salt bridge is significantly subject to the influence of pressure (76-70) both also become stable at high temperature. This bridge becomes also stable at high temperatures in NIP7-FUR. Both models also contain salt bridges, whose temperature sensitivity considerably, yet insignificantly, deviates from 1 (St <0.8 or St >1.2). Most become stable when temperature rises (St >1). Thus, elevation in temperature for a number of salt bridges in our models changes their persistence mainly in favor of its increase.
A more detailed consideration of the salt bridge stability may reveal certain interesting details. Bridge 131-116, whose stability, in the case of the deep-water protein, is under the significant influence of temperature, and is formed by the side chain groups of residues ASP131 and ARG116 for the NIP7-ABY model. In the case of the shallow-water protein, aspartic acid is substituted by glutamic acid at position 131 and no significant influence of temperature is observed. Bridge 148-109 is formed by the side chain groups of residues ARG148-ASP109 for the NIP7-ABY model. In the case of the shallow-water protein, arginine is substituted by lysine at position 148 and, as a result, this bridge turns out to be unstable for the NIP7-FUR model. Therefore, the structure of the side chain radicals has a strong impact on the stability of these salt bridges.
Influence of protein type on the dynamics of protein structure
To estimate the influence amino acid substitutions exerted on the structural properties of the NIP7 proteins in comparison with factors such as temperature and pressure, we performed three-way ANOVA for each residue. The results are given in Additional file 10.
Analysis of the influence temperature exerts on the RMSDL in the two models shows that the local structure deviates significantly from the crystal 2P38:A for 77 of the 143 positions we analyzed (53%). Pressure exerts a significant influence on the RMSDL at only 9 positions of the protein (6%), while protein type exerts a significant influence on the conformational deviations of the polypeptide chain for 91 positions, 63% (Figure 8).
The most characteristic differences are observed for the C-terminal domain of the protein (Figure 8). Furthermore, significant differences in the RMSDL parameter depending on model type are observed for regions 16-28, 40-45, 67-77 and position 5. As seen, the deviations from the crystal structure are greater for the NIP7-FUR than the NIP7-ABY model.
A similar analysis was carried out for the RMSF parameters of the amino acid residues (Additional file 10). The results showed that a rise in temperature caused a significant increase in the RMSF for all, without exception, amino acid residues of the NIP7 protein. Pressure exerted a significant influence on change in 26 amino acid residues (16%). Type of model exerted an influence on change in fluctuations for 56 (36%) of the amino acid residues.
Significant differences between the RMSF values for the residues of the two models were observed both in the case when fluctuations were higher for the NIP7-ABY model and, conversely, for the NIP7-FUR model (Figure 9). Residues at positions of the N-terminal domain (9, 11, 12, 17, 20, 21, 23, 24-26) corresponding to the α1 helix and the β2 strand, also position 34 (the loop between β3 and β4), also the C-terminal domain (96-107, 115-118, 154-155) are referred to the first group. Positions 65-67, 78-85, 89, 120-128, 131-133, 138-139 are referred to the second group.
Effect of temperature and pressure on the structure and dynamics of Nip7
In our work we compare the dynamics of two Nip7 protein models, NIP7-ABY from P. abyssi, a deep water organism, and NIP7-FUR, from P. furiosus, a shallow water organism, at different temperatures and pressures.
The hallmark feature of our models is their stability throughout the simulation period (40 ns). It is well maintained at both elevated temperatures and pressures. The RMSDs from the crystal structure (2P38:A) are, on average, smaller for the NIP7-ABY than the NIP7-FUR model. It could be remembered, however, that the NIP7-FUR model is not crystallographic, but it has been rather produced by homology modeling, which could indeed justify a larger RMSD. The secondary structure of the model changes slightly, thereby providing evidence that the structural differences between the two proteins are related to change in loop conformation and in mutual disposition of the secondary structure elements.
Exposure to high temperature for the two proteins results in change in the equilibrated conformations of the polypeptide chain, on the one hand (seen when comparing the RMSDL profiles in Figure 5), and in increase in its fluctuations, on the other hand.
Conformational changes at increasing temperature are related to decrease in structural deviations from the crystal structure in the C-terminal domain, including the positions that form the interaction with RNA. The deviations increase in the regions of the helices α2-α3 in the N-terminal domain. It should be noted that the specific influence of elevated temperature on the conformation of the protein, on its active centers in particular, is one of the possible mechanisms providing optimal protein activity in this region. For example, the important role of temperature-dependent conformational transitions of protein, at its active site particularly, in the enzymatic activity of the NADH oxidase from Thermus thermophilus has been well demonstrated at high temperatures . Also it was demonstrated that some mutation may increase and some of them may decrease the thermostability of thioredoxins from Escherichia coli and Bacillus acidocaldarius; the thermostability of these proteins was revealed to depend on ionic interactions between the thermolabile regions . Such conformational changes may be important to the Nip7 functions as well.
The comparison of the protein models obtained at 373 and 300 K demonstrates that at high temperature the fluctuations in the polypeptide chain increase, although the compactness of the protein structure is preserved (Figure 3). The increase in the fluctuations of the polypeptide chain is non-uniform. The regions subject to the fluctuations are mainly loops and the terminal regions of the protein. Like in the case of the RMSDL, the region α2-α3 proved to be the most subject to increase in the RMSF among the internal regions of the chain. It is pertinent to note that a similar effect of high temperatures on the conformation of the polypeptide chain is characteristic of the high temperature protein models -.
The effect of increased pressure on protein structure is manifold. The structural deviations of the NIP7-ABY from the conformation of the Nip7 model based on the X-ray diffraction structural data at high temperatures decrease as pressure increases. This trend holds true for the NIP7-FUR model (Table 1). The increased pressure results in a decrease in the accessible surface area. This decrease is characteristic of NIP7-FUR for which the influence pressure exerts on the SAS value is significant (Additional file 2). A significant decrease in the SAS area is observed for the NIP7-ABY model, only at high temperature. In this way, increase in pressure makes the Nip7 structure more compact, an effect more expressed at high temperature.
The increase in pressure produces structural rearrangements of the protein. Thus, for the NIP7-ABY model, this increase results in a decrease in the local deviations of protein conformation from the crystal structure. This is in complete agreement with the conclusions that the pressure affects essentially slower motions which imply structural rearrangements of the protein globule . The pressurization effect is more pronounced for the NIP7-ABY model and it affects the conformation non-uniformly, increasing the RMSDL in some regions and decreasing it in others (Figure 5B). As for the NIP7-FUR model most conformational changes occur under high pressure in the N-terminal domain.
The MD simulations of protein structures subject to increased pressure have yielded abundant evidence of reduced fluctuations in the polypeptide chain due to pressurization -. They are convergent in demonstrating that, if the protein globule does not denaturate, pressure can stabilize the polypeptide chain. An increase in protein stability resulting from an increase in pressure has been demonstrated experimentally for the exemplary glutamate dehydrogenases from the hyperthermophilic archaea Pyrococcus furiosus and Thermococcus litoralis. The published data support ours: indeed, an increase in pressure stabilized significantly the protein globule, also the protein complex formed by the enzymes. The authors suggested that such a stabilization may arise through changes in the stability of the native states due to reduced fluctuations in the polypeptide chain at high pressures and temperatures . However, our results demonstrate that the magnitude of the effect can be different for protein parts and for proteins from different organisms.
Functional implications for Nip7
Our current results demonstrate that the C-terminal domain is subject to larger structural displacements and fluctuations than the N-terminal domain during the MD simulations. This domain is plastic possibly because it is small, with just ~60 amino acid residues, and is stabilized predominantly by the hydrophobic core (the few salt bridges in this domain coordinate the positions of the loops and the terminal helical regions). The plasticity of the DNA/RNA-binding domains makes non-specific nucleotide-binding feasible -. Such a flexible structure, in contrast to the strictly coordinated components, which obey the key-lock rule, provides the possible binding to the poly-U RNA and poly-AU RNA sequences with a weak secondary structure ,. Therefore, the plasticity that we currently observed for the Nip7 PUA-domain may presumably be its significant functional property.
The Nip7 PUA-domain has an interesting feature: the formation of salt bridges by the residues shaping the interaction with RNA. Possibly, this effect allows the conformation of chain side groups of these residues to stabilize so as to facilitate binding to RNA. It is encouraging that a similar effect of the stabilization of the side group amino acids at the expense of the formation of rigidifying salt bridges has been observed in the active centre of acylphosphatase from Pyrococcus horikoshii. This is pertinent to our current observations: such a mechanism is to a great extent characteristic of the model for the deep-water NIP7-ABY protein.
The N-terminal domain, on the whole, is subject to smaller fluctuations compared to the C-terminal domain (Additional file 6). In this connection, the structurally unstable segments of this domain (positions 30-40, 49-59, 69-79) are outstanding. Two of these segments are β-hairpins (positions 30-40, 69-79), the deviations from the crystal structure in these segments occur without significant changes in the secondary structure, so that they bend in the direction of the α1-helix in the two models. These conformational changes may result from interactions in the salt bridges network (GLU10-ARG4, ARG37-GLU10, GLU75-ARG37) bringing close together the β-sheet and the α1 helix of the N-terminal domain. This interaction network may, in turn, become stabilized through the voluminous hydrophobic nucleus formed by hydrophobic residues in the centre of the β1-β5 strands.
In terms of the putative functional role in the N-terminal domain, it appears worthwhile to consider the loop between the second and the third α-helices (positions 49-59). The loop is contiguous with a lengthy nonpolar region at the surface of the N-terminal domain, which extends from the α2 helix to the groove between the N- and C-terminal domains so that is comes to lie on the side opposite to the polar portion of the Nip7 surface (Figure 1) . The conformational changes in the α2-α3 region may be due to the high twisting tension of the main chain in it. It undergoes a sharp bend at the very end of the short α3-helix, then forms one helix turn, which is converted into the β5 strand structure (Figure 1A). The positions of the N- and C- terminal residues of this fragment are fixed by the globule of the N-terminal domain. Then, the loop itself comes to lie aside from the globular domain without imposing steric constraints on the conformational changes. Additional destabilizing factors for the conformation of this region are alternating polar and hydrophobic (positions 46, 49, 53, 55) residues, which, as a result of conformational changes in the main chain of the α2-α3 segment, can by turns face solvent. It may be assumed that this bend together with the lengthy region of the nonpolar surface may play an important role in the interaction of the Nip7 protein with the exosomal protein partner ,. Interestingly, prediction of the protein-protein interaction sites using 3D structure by the SPPIDER web-server demonstrated that most residues of this region could be involved in protein-protein interactions with high probability both for the NIP7-ABY and NIP7-FUR models (Additional file 11).
The structural lability of proteins may accomplish a role of consequence in molecular recognition. For example, the regions of the polypeptide chain of the monomers devoid of an orderly structure can provide the protein interaction  and adopt an orderly conformation in the process of binding . With respect to the Nip7 protein, intrinsic flexibility appears to be more likely  because it ensures the presence of a conformational ensemble required for the formation of RNA-protein and protein-protein interactions .
Comparison of the protein dynamics between the shallow- and deep-water organisms
The sequences of P. abyssi and P. furiosus Nip7 in the 3D structure differ by 47 substitutions (~30%; Figure 1B); 19 of them result in radical changes of the physico-chemical properties of amino acids. These substitutions may be assigned to categories I-V.
Substitutions resulting in the replacement of a polar residue in P. abyssi by a nonpolar in P. furiosus (7 substitutions): 2 (R → I), 21 (Y → F), 58 (Y → F), 59 (S → A), 66 (T → M), 80 (N → A), 119 (K → V).
P. abyssi → P. furiosus substitutions, which result in replacement of a charged amino acid by an uncharged one at retained polarity (5 substitutions): 23 (E → T), 27 (E → N), 44 (E → N), 144 (K → S), 154 (K → T).
Substitutions of nonpolar amino acids in P. abyssi by polar ones in P. furiosus (1 substitution): 140 (L → R).
P. abyssi → P. furiosus substitutions of the uncharged polar side group by the charged side group: 121 (Q → K).
Substitutions P. abyssi → P. furiosus resulting in the change of the charge sign to the opposite one (5 substitutions): 88 (K → D), 113 (E → K), 117 (K → Q), 147 (R → E).
As a result of such substitutions, the hydrophobic portion of SAS increased in the shallow-water organism compared with its deep-water counterpart. The question then was: How could this influence the stability, dynamics and functional properties of the proteins? True, both proteins, in all the conditions we examined, proved to be stable and had a similar conformation. However it was perplexing to single out a parameter in the dynamics that could be significantly influenced by these substitutions.
It will be remembered that the majority of substitutions that belong to categories I and II lie in the N-terminal domain of the protein, they are in contact with its surface region, which contains mobile residues (the β2-strand and the β2-β3 loop), also residues at positions α2-α3 (Figure 8) and it may presumably be involved in protein-protein interactions (Figure 10; Additional file 11).
The substitutions we observed in the N-terminal domain may, probably, somehow modify the interaction character of this region with its molecular partner. Substitutions in this stretch are associated with increase in hydrophobicity of residues in the shallow-water organisms. Changes of such kind can promote the stabilization of protein interaction of Nip7 and its partner. This is because hydrophobic interactions contribute considerably to the stability of protein-protein interfaces . Another fact will be remembered here. When pressure close to atmospheric and temperature rises, hydrophobic interactions favor closer interactions of nonpolar protein regions , which can additionally stabilize the Nip7-exosome protein complex in P. furiosus. If so, the observed substitutions may truly be of adaptive character. We have previously analyzed  the indices of asymmetry for amino acid substitutions in hundreds of Pyrococcus proteins starting from its common deep-sea ancestor P. furiosus, P. horikoshii and P. abyssi to the extant shallow-sea organism P. furiosus. According to the results, the values of the indices are positively correlated with the scales of amino acid hydrophobicity and negatively with the polarity scales indicating that polar residues were predominantly substituted by nonpolar when the shallow-sea habitats replaced the deep-sea ones. However, the issue as to whether this trend is a reflection of the stabilization of protein-protein interactions or of the entire structure requires further study.
Furthermore, the character of changes in fluctuations for the chain stretch in the α2-α3 region at increasing pressure is drastically different for the two proteins: they significantly increase in the shallow-water protein and decrease, conversely, in its deep-water counterpart (Figure 7А). A general decrease in fluctuations throughout the entire N-terminal domain at increasing pressure is also characteristic of the deep-sea protein. The magnitude of fluctuations for efficient binding must be within definite limits and their general decrease at reduction in pressure is also additive in the deep-water protein.
Substitutions exert a significant influence on the conformational parameters in proteins connecting two protein domains (the RMSDL and the RMSF parameters) in the β6-α4 region. This is the stretch that joins together two protein domains. Probably, mutations in it favor change in the mutual orientation of the domains and result in an increase in the structural deviations in the C-terminal domain of the protein in the NIP7-FUR as compared with the NIP7-ABY model (Figure 8).
As for the RNA-binding domain, radical changes that could alter residue charge are characteristic of the shallow-water protein. As a consequence, the formation of certain salt bridges, those between residues involved in RNA recognition in particular, is impaired. It is also worth to note, that the number of residues associated with significant differences in RMSDL/RMSF parameters between the shallow- and deep-water protein models are higher in the RNA-binding domain. The differences we observed in the course of MD simulations may reflect the different ways and means RNA interacts with the Nip7 proteins at high and low pressures.
The MD simulations we performed currently for the Nip7 protein from the shallow-water (P. furiosus) and deep-water (P. abyssi) marine hyperthermophylic archaea demonstrated that their structures are stable at different temperatures (300 and 373 K) and pressures (0.1-100 MPa). Increase in temperatures and pressures caused conformational changes in about the same protein regions, in the loops first of all. In the examined ranges of temperatures and pressures, increase in temperature had a more pronounced effect on changes in the dynamic properties of the protein globule than increase in pressure. A number of highly mobile segments of protein globule presumably involved in protein-protein interactions were identified. Substitutions of the polar residues in the deep-water organisms by hydrophobic in the shallow-water ones were observed in some of the regions. These replacements may be evidence of change in the interaction pattern of the Nip7 protein and the proteins of the exosome complex in the course of the divergence of the shallow-water organisms from the deep-water ones. Taken together, the results of such a complex comparative analysis of the MD of proteins from organisms belonging to different ecological niches may be useful in refining ideas of their functioning and adaptation to extreme environmental conditions.
The model of the P. abyssi Nip7 protein was retrieved from the PDB , identifier 2P38 . The atomic model of the P. abyssi Nip7 contains two independent monomers and 252 water molecules. The monomers are symmetric with respect to each other, with a head-to-tail orientation, they form a dimer without showing structural differences (the RMSD of the Cα atoms between two monomers is 0.11 Å). The contact surface between the monomers is about 700-Å2 per molecule, or 8% of its total surface. These estimates suggest that the dimeric structure might have resulted from crystallization and, if so, it may be biologically irrelevant . Hence, 2P38:A is the only monomer included in analysis. The first four and the last seven amino acid residues not located in the experiment were discarded. As a result, the total length of the P. abyssi Nip7 protein model is 155 amino acid residues. This model is referred to as NIP7-ABY, the amino acid residues are numbered according to their order numbers in the structure (1VAL-155ASP).
The P. furiosus Nip7 model (NIP7-FUR) was built using the method of reconstruction based on homology. It was retrieved from the ModBase databank , ID Q8TZP7. The P. abyssi and P. furiosus Nip7 alignments contain neither deletions nor insertions and, hence, NIP7-FUR reliability is ensured (Figure 1B). The number of residues of the NIP7-FUR model is 155, their numbering corresponds to the one for NIP7-ABY (1ILE-155ASP). Based on the Cα atoms, the RMSD of the initial structures of NIP7-FUR and NIP7-ABY is 0.1 Å.
Here we examine the MD simulations of the Nip7 protein models at various pressures (0.1, 50, and 100 MPa) and at temperatures 300 and 373 K.
Preparation of models for MD simulation
The MD of the protein 3D structure is simulated using the program package GROMACS  4.5.3. The algorithm for preparation of the models is as follows. To remove undesirable tensions, the structure energy is minimized in vacuum. Then, the molecule is placed in a cubic cell (edge 200 Å) containing water molecules in the SPCE model. Then, ion molecules are added to the cell to neutralize the charge. The NIP7-ABY model consists of 2 585 protein atoms, 32 900 water molecules, and 1 ion. The NIP7-FUR model is composed of 2 544 protein atoms, 32 900 water molecules, and 2 ions. Thereafter, the protein solvent system is minimized for 2 000 steps. There are two different techniques for further preparation of the models. In the case of simulation at 300 K, there are two preparatory steps of simulation at the given pressure-temperature values, 200-ps each; in that of simulation at 373 K, the number of preparatory simulations is nine, with a gradual increase in temperature by 8 K at each simulation. Pressure in the system is set at simulation onset.
MD simulations are performed using the LINCS  algorithm in the amber99sb-ildn force field . Electrostatic interactions are computed using the Partial Mesh Ewald (PME) algorithm. Pressure coupling was done using Parrinello-Rahman method . Temperature coupling was done using v-rescale method .
The simulation time was 40 ns. To take into account slower dynamics of proteins at high pressure due to increase in viscosity of water for each set of parameters we perform 5 independent trajectory runs. As a result, 30 experiments are considered for each protein model.
Data were saved every 80 ps. We collected statistics for protein conformation parameters for model structures after equilibration from 20 to 40 ns.
Structure similarity analysis
We characterized changes in the structures by the Cα atoms RMSDs and the solvent accessible surface area of the residues (the SAS); the secondary structure parameters based on the DSSP algorithm  were calculated using the GROMACS package with the help of the g_rms, g_sas and do_dssp programs, respectively.
To define the regions of the protein structure most susceptible to pressure/temperature impact, in a sliding window of the size q + 1, the local deviation of the structure for each residue, the RMSDL after superposition of two protein structures using all Cα atoms was calculated as:
where v, w are two structures superposed to minimize their Cα RMSD; i is the index of residue for which RMSDL is calculated; v jk , w jk (k = x, y, z) are x, y, z coordinates of the Cα atoms of the aligned residues j in each of the v, w structures; the q parameter was set to 10 yielding a window of 11 amino acids to calculate the RMSDL.
Secondary structure analysis
To analyze the stability of the secondary structure during the simulations, we analyze the 3D model at each step using the DSSP program . A hundred 3D structures obtained at the last steps of the MD simulation were used as input for the DSSP program. The occurrence frequencies of a secondary structure of a particular type at a concrete position in the secondary structure were determined. The total secondary structure for a given trajectory included the type of the secondary structure with the highest occurrence frequency. The procedure was performed on all the 24 trajectories.
To estimate the variability in the secondary structure types in the 12 trajectories for each of two models we used the Protein Variability Server .
Analysis of the fluctuation in the polypeptide chain
Fluctuations in the polypeptide chain were characterized by the Cα atom RMSFs computed using the GROMACS package program g_rmsf. The average Cα atoms RMSF value in Angstroms was computed for all the residues in the Nip7 model and N- and C- terminal domains. To determine the significance of the differences between the average Cα atoms RMSFs of the different trajectories or models, we used two-way ANOVA (see below).
Dependence of the models’ conformational parameters on temperature, pressure and protein sequence
To determine the influence of temperature, pressure and amino acid substitutions on the protein parameters, we used ANOVA of two types .
First, to estimate the influence of pressure and temperature on the protein structural characteristics (for example, Rg, SAS, salt bridge persistence or RMSF of the i-th residue) we used two-way ANOVA. In this analysis the influence of two factors was examined: temperature (T = 300, 373 K) and pressure (P = 0.1, 50, 100 MPa). In each set of modeling parameters, we estimated the means of this characteristic on the basis of the results of analysis of 5 runs. In this way, for the structural characteristic we analyzed the ANOVA table with r = 2 (temperature), c = 3 (pressure), the same number of observations for each parameter set, whose sum total was n = 30 (2 ± 3 ± 5). We tested the hypothesis that the means of the characteristic are independent of temperature and pressure. The critical F-values were set at α = 0.05 for F (P) = 3.4 (df = 2), temperature F (T) = 4.26 (df = 1), and for interaction of these factors F (PxT) = 3.4 (df = 24).
Second, to estimate the effect of increased temperature on the conformational parameters of the protein, their index of sensitivity (St) was computed. To this end, we calculated the average value of a model’s conformational parameter, X (for example, the RMSF of the i-th residue), for the trajectories obtained at high temperatures (X T=373 ) and for those obtained at low temperatures (X T=300 ). The index of sensitivity was expressed as . The St value close to unity meant that increased temperature had no considerable effect on the value of the conformational parameter compared to its value at low temperature. The St value considerably larger than unity implied that increased temperature increased the value of the structural parameter compared to its value at low temperature. The St value considerably smaller than unity imply that increase in temperature reduces the value of the structural parameter. A number of conformational parameters, such as the RMSDL i , the local structural deviations of the polypeptide chain for the i- th residue, the RMSF i , the root-mean-square fluctuations in the i- th residue, and salt bridge persistence (see below) were considered.
Similarly, the sensitivity of the structural parameters to increased pressure Sp was estimated. We determined the value of the parameters for the trajectories at high pressures (50 and 100 MPa), X P-High , atmospheric and moderate pressure (0.1 MPa), X P-low . The sensitivity to pressure was expressed as .
Third, in analysis of the structural characteristics of the amino acid residues the RMSDL, the RMSF, we estimated the influence of protein type (NIP7-ABY, NIP7-FUR); in such a way, we treated the three-way ANOVA model that considered temperature, pressure and protein type (r = 2, c = 3, l = 2; n = 60). The additional factor (protein type) reflected the influence amino acid substitutions exerted on the RMSDL and RMSF parameters for each of the residues. It should be stipulated, however, that such an analysis did not provide information about the influence of specific amino acid substitutions on changes in a conformational parameter of a residue. The critical F-values in this analysis were set at α = 0.05 for pressure F (P) = 3.4 (df = 2), for temperature F (T) = 4.26 (df = 1), and for interaction of these factors F (PxT) = 3.4 (df = 24).
Analysis of salt bridges
To identify salt bridges, we used the criterion according to which the distance between two oppositely changed atoms should not exceed 4 Å . Salt bridge persistence was expressed as the proportion of structures in the equilibrium region of trajectories (20-40 ns) for which the salt bridge was identified.
The sensitivity of the salt bridge to temperature (St) and pressure (Sp) was defined as the ration of the salt bridge persistence at high and low temperature and high and low pressure. The St (Sp) value close to unity meant that increased temperature (pressure) had no considerable effect on the persistence of the salt bridge. The St (Sp) value considerably larger than unity implied that increased temperature (pressure) increased the persistence of the salt bridge compared to its value at low temperature. The St (Sp) value considerably smaller than unity imply that increase in temperature (pressure) reduces the persistence of the salt bridge.
The dependence of salt bridge persistence on temperature and pressure was estimated using two-way ANOVA (see above).
Prediction of the protein-protein interactions
The Accelrys Discovery Studio Visualizer was applied to visualize the structure of the models (http://accelrys.com/products/discovery-studio/).
Statistical analysis of the data was performed using the Excel and Statistica 6.0 package.
Availability of supporting data
All the supporting data are included as additional files.
KEM and NAA performed the MD simulations of the Nip7 proteins. KEM, DAA and YNV performed the analysis of the MD results. EVB, NAK and DAA initiated the comparative molecular dynamics study of the Nip7 proteins and participated in its coordination. All authors read and approved the final manuscript.
Radius of gyration
Root mean squared deviation
Root mean squared fluctuation
Local root mean squared deviation
Solvent accessible surface
Model of the Nip7 protein of Pyrococcus abyssi
Model of the Nip7 protein of Pyrococcus furiosus
Rothschild LJ, Mancinelli RL: Life in extreme environments. Nature 2001, 409(6823):1092–1101. 10.1038/35059215
Daniel I, Oger P, Winter R: Origins of life and biochemistry under high-pressure conditions. Chem Society Rev 2006, 35(10):858–875. 10.1039/b517766a
Brooks AN, Turkarslan S, Beer KD, Yin Lo F, Baliga NS: Adaptation of cells to new environments. Wiley Interdiscip Rev Syst Biol Med 2011, 3(5):544–561. 10.1002/wsbm.136
Sterner R, Liebl W: Thermophilic adaptation of protein. Crit Rev Biochem Mol Biol 2001, 36: 39–106. 10.1080/20014091074174
Makarova KS, Wolf YI, Koonin EV: Potential genomic determinants of hyperthermophily. Trends Genet 2003, 19(4):172–176. 10.1016/S0168-9525(03)00047-7
Berezovsky IN, Shakhnovich EI: Physics and evolution of thermophilic adaptation. Proc Natl Acad Sci U S A 2005, 102(36):12742–12747. 10.1073/pnas.0503890102
McDonald JH: Temperature adaptation at homologous sites in proteins from nine thermophile-mesophile species pairs. Genome Biol Evol 2010, 2: 267. 10.1093/gbe/evq017
Jollivet D, Mary J, Gagniere N, Tanguy A, Fontanillas E, Boutet I, Hourdez S, Segurens B, Weissenbach J, Poch O, Lecompte O: Proteome adaptation to high temperatures in the ectothermic hydrothermal vent Pompeii worm. PLoS One 2012, 7(2):e31150. 10.1371/journal.pone.0031150
Di Giulio M: A comparison of proteins from Pyrococcus furiosus and Pyrococcus abyssi : barophily in the physicochemical properties of amino acids and in the genetic code. Gene 2005, 346: 1–6. 10.1016/j.gene.2004.10.008
Simonato F, Campanaro S, Lauro FM, Vezzi A, D’Angelo M, Vitulo N, Bartlett DH: Piezophilic adaptation: a genomic point of view. J Biotechnol 2006, 126(1):11–25. 10.1016/j.jbiotec.2006.03.038
Campanaro S, Treu L, Valle G: Protein evolution in deep sea bacteria: an analysis of amino acids substitution rates. BMC Evol Biol 2008, 8(1):313. 10.1186/1471-2148-8-313
Oger PM, Jebbar M: The many ways of coping with pressure. Res Microbiol 2010, 161(10):799–809. 10.1016/j.resmic.2010.09.017
Podar M, Reysenbach AL: New opportunities revealed by biotechnological explorations of extremophiles. Curr Opin Biotechnol 2006, 17(3):250–255. 10.1016/j.copbio.2006.05.002
Frock AD, Kelly RM: Extreme thermophiles: moving beyond single-enzyme biocatalysis. Curr Opin Chem Eng 2012, 1(4):363–372. 10.1016/j.coche.2012.07.003
Liszka MJ, Clark ME, Schneider E, Clark DS: Nature versus nurture: developing enzymes that function under extreme conditions. Annu Rev Chem Biomol Eng 2012, 3: 77–102. 10.1146/annurev-chembioeng-061010-114239
Singh OV: Extremophiles: Sustainable Resources and Biotechnological Implications. Wiley-Verlag, New York; 2012.
Fiala G, Stetter KO: Pyrococcus furiosus sp. nov. represents a novel genus of marine heterotrophic archaebacteria growing optimally at 100°C. Arch Microbiol 1986, 145: 56–61. 10.1007/BF00413027
Erauso G, Reysenbach A-L, Godfroy A, Meunier J-R, Crump B, Partensky F, Baross JA, Marteinsson V, Barbier G, Pace NR, Prieur D: Pyrococcus abyssi sp. nov., a new hyperthermophilic archaeon isolated from a deep-sea hydrothermal vent. Arch Microbiol 1993, 160: 338–349. 10.1007/BF00252219
González JM, Masuchi Y, Robb FT, Ammerman JW, Maeder DL, Yanagibayashi M, Tamaoka J, Kato C: Pyrococcus horikoshii sp. nov., a hyperthermophilic archaeon isolated from a hydrothermal vent at the Okinawa Trough. Extremophiles 1998, 2: 123–130. 10.1007/s007920050051
Zanchin NIT, Goldfarb DS: Nip7p interacts with Nop8p, an essential nucleolar protein required for 60S ribosome biogenesis, and the exosome subunit Rrp43p. Mol Cell Biol 1999, 19: 1518–1525.
Bassler J, Grandi P, Gadal O, Lessmann T, Petfalski E, Tollervey D, Lechner J, Hurt E: Identification of a 60S preribosomal particle that is closely linked to nuclear export. Mol Cell Biol 2001, 8: 517–529.
Luz JS, Ramos CR, Santos MC, Coltri PP, Palhano FL, Foguel D, Oliveira CC: Identification of archaeal proteins that affect the exosome function in vitro. BMC Biochem 2010, 11(1):22. 10.1186/1471-2091-11-22
Coltri PP, Guimaraez BG, Granato DC, Luz JS, Teixeira EC, Oliveira CC, Zanchin NIT: Structural insights into the interaction of the Nip7 PUA domain with polyuridine RNA. Biochemistry 2007, 46: 14177–14187. 10.1021/bi7015876
Pérez-Arellano I, Gallego J, Cervera J: The PUA domain? a structural and functional overview. FEBS J 2007, 274(19):4972–4984. 10.1111/j.1742-4658.2007.06031.x
Aravind L, Koonin E: Novel predicted RNA-binding domains associated with the translation machinery. J Mol Evol 1999, 48: 291–302. 10.1007/PL00006472
Hallberg BM, Ericsson UB, Johnson KA, Andersen NM, Douthwaite S, Nordlund P, Beuscher AE, Erlandsen H: The structure of the RNA m5C methyltransferase YebU from Escherichia coli reveals a C-terminal RNArecruiting PUA domain. J Mol Biol 2006, 360: 774–787. 10.1016/j.jmb.2006.05.047
Gunbin KV, Afonnikov DA, Kolchanov NA: Molecular evolution of the hyperthermophilic archaea of the Pyrococcus genus: analysis of adaptation to different environmental conditions. BMC Genomics 2009, 10(1):639. 10.1186/1471-2164-10-639
Barlow DJ, Thornton JM: Ion-pairs in proteins. J Mol Biol 1983, 168(4):867–885. 10.1016/S0022-2836(83)80079-5
Merkley ED, Daggett V, Parson WW: A temperature-dependent conformational change of NADH oxidase from Thermus thermophilus HB8. Proteins 2012, 80(2):546–555. 10.1002/prot.23219
Polyansky AA, Kosinsky YA, Efremov RG: Correlation of local changes in the temperature-dependent conformational flexibility of thioredoxins with their thermostability. Russ J Bioorg Chem 2004, 30: 421–430. 10.1023/B:RUBI.0000043784.51859.41
Martinez R, Schwaneberg U, Roccatano D: Temperature effects on structure and dynamics of the psychrophilic protease subtilisin S41 and its thermostable mutants in solution. Protein Eng Des Sel 2011, 24(7):533–544. 10.1093/protein/gzr014
Lee KJ: Molecular dynamics simulations of a hyperthermophilic and a mesophilic protein L30e. J Chem Inf Model 2011, 52(1):7–15. 10.1021/ci200184y
Priyakumar UD, Ramakrishna S, Nagarjuna KR, Reddy SK: Structural and energetic determinants of thermal stability and hierarchical unfolding pathways of hyperthermophilic proteins, Sac7d and Sso7d. J Phys Chem B 2010, 114(4):1707–1718. 10.1021/jp909122x
Tiberti M, Papaleo E: Dynamic properties of extremophilic subtilisin-like serine-proteases. J Struct Biol 2011, 174(1):69–83. 10.1016/j.jsb.2011.01.006
Calandrini V, Kneller GR: Influence of pressure on the slow and fast fractional relaxation dynamics in lysozyme: a simulation study. J Chem Phys 2008, 128: 065102. 10.1063/1.2828769
Capece L, Marti MA, Bidon-Chanal A, Nadra A, Luque FJ, Estrin DA: High pressure reveals structural determinants for globin hexacoordination: neuroglobin and myoglobin cases. Proteins 2009, 75(4):885–894. 10.1002/prot.22297
McCarthy AN, Grigera JR: Effect of pressure on the conformation of proteins. A molecular dynamics simulation of lysozyme. J Mol Graph Model 2006, 24(4):254–261. 10.1016/j.jmgm.2005.09.007
Laurent AD, Mironov VA, Chapagain PP, Nemukhin AV, Krylov AI: Exploring structural and optical properties of fluorescent proteins by squeezing: modeling high-pressure effects on the mStrawberry and mCherry red fluorescent proteins. J Phys Chem B 2012, 116(41):12426–12440. 10.1021/jp3060944
Sun MM, Tolliday N, Vertiani C, Robb FT, Clark DS: Pressure-induced thermostabilization of glutamate dehydrogenase from the hyperthermophile Pyrococcus furiosus. Protein Sci 1999, 8(5):1056–1063. 10.1110/ps.8.5.1056
Sun MM, Caillot R, Mark G, Robb FT, Clark DS: Mechanism of pressure-induced thermostabilization of proteins: studies of glutamate dehydrogenases from the hyperthermophile Thermococcus litoralis. Protein Sci 2001, 10(9):1750–1757. 10.1110/ps.4001
Kalodimos CG, Biris N, Bonvin AMJJ, Levandoski MM, Guennuegues M, Boelens R, Kaptein R: Structure and flexibility adaptation in nonspecific and specific protein-DNA complexes. Science 2004, 305(5682):386–389. 10.1126/science.1097064
Brown C, Campos-León K, Strickland M, Williams C, Fairweather V, Brady RL, Cramp MP, Gaston K: Protein flexibility directs DNA recognition by the papillomavirus E2 proteins. Nucleic Acids Res 2011, 39(7):2969–2980. 10.1093/nar/gkq1217
Boehr DD: Promiscuity in protein-RNA interactions: conformational ensembles facilitate molecular recognition in the spliceosome. Bioessays 2012, 34(3):174–180. 10.1002/bies.201100152
Lam SY, Yeung RC, Yu TH, Sze KH, Wong KB: A rigidifying salt-bridge favors the activity of thermophilic enzyme at high temperatures at the expense of low-temperature activity. PLoS Biol 2011, 9(3):e1001027. 10.1371/journal.pbio.1001027
Fong JH, Shoemaker BA, Garbuzynskiy SO, Lobanov MY, Galzitskaya OV, Panchenko AR: Intrinsic disorder in protein interactions: insights from a comprehensive structural analysis. PLoS Comput Biol 2009, 5(3):e1000316. 10.1371/journal.pcbi.1000316
Uversky VN, Dunker AK: The case for intrinsically disordered proteins playing contributory roles in molecular recognition without a stable 3D structure. F1000 Biol Rep 2013, 5: 1. 10.3410/B5-1
Janin J, Sternberg MJ: Protein flexibility, not disorder, is intrinsic to molecular recognition. F1000 Biol Rep 2013, 5: 2. 10.3410/B5-2
Boehr DD, Nussinov R, Wright PE: The role of dynamic conformational ensembles in biomolecular recognition. Nat Chem Biol 2009, 5(11):789–796. 10.1038/nchembio.232
Jones S, Thornton JM: Principles of protein-protein interactions. Proc Natl Acad Sci U S A 1996, 93(1):13–20. 10.1073/pnas.93.1.13
Ferrara CG, Ghara O, Grigera JR: Aggregation of nonpolar solutes in water at different pressures and temperatures: the role of hydrophobic interaction. J Chem Phys 2012, 137(13):135104. 10.1063/1.4755752
Afonnikov DA, Medvedev KE, Gunbin KV, Kolchanov NA: Important role of hydrophobic interactions in high-pressure adaptation of proteins. Dokl Biochem and Biophys 2011, 438(1):113–116. 10.1134/S160767291103001X
Berman HM, Westbrook J, Feng Z, Gilliland G, Bhat TN, Weissig H, Shindyalov IN, Bourne PE: The protein data bank. Nucleic Acids Res 2000, 28: 235–242. 10.1093/nar/28.1.235
Pieper U, Eswar N, Davis FP, Braberg H, Madhusudhan MS, Rossi A, Marti-Renom M, Karchin R, Webb BM, Eramian D, Shen MY, Kelly L, Melo F, Sali A: ModBase, a database of annotated comparative protein structure models, and associated resources. Nucleic Acids Res 2006, 34: D291-D295. 10.1093/nar/gkj059
Hess B, Kutzner C, Van Der Spoel D, Lindahl E: GROMACS 4: algorithms for highly efficient, load-balanced, and scalable molecular simulation. J Chem Theory Comput 2008, 4(3):435–447. 10.1021/ct700301q
Hess B, Bekker H, Berendsen HJC, Fraaije JGEM: LINCS: a linear constraint solver for molecular simulations. J Comp Chem 1997, 18(12):1463–1472. 10.1002/(SICI)1096-987X(199709)18:12<1463::AID-JCC4>3.0.CO;2-H
Lindorff-Larsen K, Piana S, Palmo K, Maragakis P, Klepeis JL, Dror RO, Shaw DE: Improved side-chain torsion potentials for the Amber ff99SB protein force field. Proteins 2010, 78(8):1950–1958.
Parrinello M, Rahman A: Polymorphic transitions in single crystals: a new molecular dynamics method. J Appl Phys 1981, 52(12):7182–7190. 10.1063/1.328693
Bussi G, Donadio D, Parrinello M: Canonical sampling through velocity rescaling. J Chem Phys 2006, 126(1):014101–014107. 10.1063/1.2408420
Kabsch W, Sander C: Dictionary of protein secondary structure: pattern recognition of hydrogen-bonded and geometrical features. Biopolymers 1983, 22(12):2577–2637. 10.1002/bip.360221211
Garcia-Boronat M, Diez-Rivero CM, Reinherz EL, Reche PA: PVS: a web server for protein sequence variability analysis tuned to facilitate conserved epitope discovery. Nucleic Acids Res 2008, 36: W35-W41. 10.1093/nar/gkn211
Kendall MG, Stuart A: The Advanced Theory of Statistics. Volume 3. Charles Griffin & Co, London; 1968.
Porollo A, Meller J: Prediction-based fingerprints of protein-protein interactions. Proteins 2007, 66: 630–645. 10.1002/prot.21248
The SB RAS Siberian Supercomputing Center and the Novosibirsk State University High-Performance Computing Center are gratefully acknowledged for providing computer facilities. This work was in part supported by the Ministry of Education and Science of the Russian Federation, RFBR (projects 11-04-01771 and 12-04-00135); Integration Projects № 130 and 39 from the SB RAS; Program № 6.6 and “Biosphere Origin and Evolution” from the RAS; EVB acknowledges the support from a grant of RSF 14-13-00834.
The authors are grateful to the referees for their helpful comments, to K.V. Gunbin for helpful discussions, to A.N. Fadeeva for translating the manuscript from Russian into English and to N.V. Kuchin for technical support.
The authors declare that they have no competing interests.
Electronic supplementary material
Additional file 3: Analysis of the RMSDL, the NIP7-ABY and NIP7-FUR models; Analysis of the secondary structure, the NIP7-ABY and NIP7-FUR models; Snapshots of the NIP7-ABY and NIP7-FUR models at different temperatures and pressures; Analysis of the Cα RMSF, the NIP7-ABY and NIP7-FUR models. (PDF 3 MB)
Additional file 6: 2-way ANOVA of the fluctuations in the C- and N-terminal domains of the NIP7-ABY and NIP7-FUR models. (XLS 52 KB)
Additional file 7: Analysis of salt bridge persistence during molecular dynamics simulation at different pressures and temperatures. (XLS 86 KB)
Additional file 11: Prediction of the protein-protein interaction sites for the NIP7-ABY and NIP7-FUR models by the SPPIDER web-server. (XLS 60 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Medvedev, K.E., Alemasov, N.A., Vorobjev, Y.N. et al. Molecular dynamics simulations of the Nip7 proteins from the marine deep- and shallow-water Pyrococcus species. BMC Struct Biol 14, 23 (2014). https://doi.org/10.1186/s12900-014-0023-z