Effect of the Lys62Ala Mutation on the Thermal Stability of BstHPr Protein by Molecular Dynamics

We analyzed the thermal stability of the BstHPr protein through the site-directed point mutation Lys62 replaced by Ala residue using molecular dynamics simulations at five different temperatures: 298, 333, 362, 400, and 450 K, for periods of 1 μs and in triplicate. The results from the mutant thermophilic BstHPrm protein were compared with those of the wild-type thermophilic BstHPr protein and the mesophilic BsHPr protein. Structural and molecular interaction analyses show that proteins lose stability as temperature increases. Mutant and wild-type proteins behave similarly up to 362 K. However, at 400 K the mutant protein shows greater structural instability, losing more buried hydrogen bonds and exposing more of its non-polar residues to the solvent. Therefore, in this study, we confirmed that the salt bridge network of the Glu3–Lys62–Glu36 triad, made up of the Glu3–Lys62 and Glu36–Lys62 ion pairs, provides thermal stability to the thermophilic BstHPr protein.


Introduction
Proteins are macromolecules made up of chains of amino acids and play an important role in the structure and function of cells.They have numerous and important functions in the organism, as they can carry out different processes such as enzymatic catalysis, molecular transport, mechanical resistance, protection against pathogens, metabolic regulation, etc. [1,2].For a protein to perform these functions, it must maintain a stable fold, which is achieved during the transcription of the mRNA information in the ribosome.The chemical nature of the amino acids and the environment in which the proteins are found are essential for suitable folding [3].
A large number of proteins have been characterized thanks to advances in the fields of molecular biology and biochemistry.Structural information and thermodynamic data have been the basis for the development of algorithms and methods for the explanation of protein-folding mechanisms [1,2].In addition, thermodynamic data of these biomolecules are essential to understanding and predicting the stability when mutations are made to improve their resistance to unfolding under the influence of external factors such as solvents, salts, temperature, etc. [4][5][6].
The study of proteins from extremophilic organisms has generated great interest in the field of protein folding, since the elucidation of this process contributes to better knowledge of the mechanisms that stabilize proteins, which is fundamental to improving their design, for example by developing more efficient enzymes that work at high temperatures [7][8][9].
Proteins from thermophilic and hyperthermophilic organisms show stable adaptations to high-temperature environments, e.g., hydrothermal vents, hot springs, mud volcanoes, etc. [10].In recent decades, studies have focused on analyzing the unusually high stability of these proteins [4][5][6][7][8][9].The knowledge gained has been used to improve proteins from other types of organisms, such as mesophiles.There are different reasons for this, the most common being that temperature is a variable used in most industrial processes, for example, in the food industry [7].Research has suggested that there are several mechanisms for increasing the thermal stability of proteins.However, none of them provide a universal solution, as it has been found that proteins activate different physicochemical pathways to withstand high temperatures and not lose their functionality within the ecosystems where their host organisms live [10].An alternative way to study and improve the thermal stability of proteins is to perform point mutations of sites, i.e., mutations can be performed by making a targeted exchange of amino acids [6,11,12].Experimental studies such as circular dichroism (CD), differential scanning calorimetry (DSC), and fluorescence spectroscopy are used to determine whether thermal resistance increases [5,13].
The thermostability analyses of the BstHPr protein have been reported in various experimental studies [12][13][14][15][16]. Their structural and thermodynamic data have been compared with those of homologous proteins from mesophilic organisms such as B. subtilis (Bs), E. coli (Ec), E. faecalis (Ef ), and M. capricolum (Mc); moderate thermophiles such S. thermophilus (St); thermophiles as T. tengcongensis (Tt); and haloalkaliphilic organisms such as B. halodurans (Bh) and O. iheyensis (Oi) [12,14,16].Stability measurements of the Bst, Bh, St, Bs, and Oi HPr proteins carried out by thermal and solvent denaturation, under different experimental conditions of pH, salinity, and temperature, have shown that BstHPr has the highest thermal stability.For example, Table 1 shows a comparison of thermodynamic data, such as free energy of stabilization (∆G S ), temperature of maximal stability (T S ), melting temperature (T m ), change in heat capacity (∆C p ), and change in enthalpy (∆H), between the BstHPr and BsHPr proteins at pH = 7.0 [12][13][14].In addition, it has been observed that the BstHPr protein forms a salt bridge between the Asp11 and Lys57 residues [12,16].The site-directed Lys57Thr mutation shows that the Asp11-Lys57 salt bridge plays an important role in the thermostability of the BstHPr protein [16].These facts have shed light on the structural behavior and the mechanisms involved in the thermal stability of the BstHPr protein but have not been conclusive, since the key molecular interactions that stabilize it at elevated temperatures are still not well understood.
Recently, our group performed molecular dynamics (MD) simulations of HPr proteins from Bacillus stearothermophilus (BstHPr) and Bacillus subtilis (BsHPr) organisms to elucidate the molecular mechanisms that provide thermal resistance to the thermophilic protein at elevated temperatures.Structural and molecular interaction results showed that the salt bridge network formed by the Glu3-Lys62-Glu36 triad is a key factor in its stability [20].To confirm this hypothesis, in this work, we present MD analysis of a mutant structure (BstHPrm) consisting of the replacement of the lysine residue at position 62 with an alanine residue, affecting the aforementioned triad.Moreover, the effects of temperature on the mutant structure of the HPr protein were analyzed and compared with those obtained from the thermophilic wild-type (BstHPr) and mesophilic (BsHPr) structures.

Structural Behavior
Figure 1 shows the root mean square deviation (RMSD), radius of gyration (Rg), and fraction of native contact (Q) trajectories for BstHPr and BstHPrm proteins at 298, 362, and 400 K.Both proteins exhibit almost the same structural behavior during the simulation time at 298 K (Figure 1a,d,g), indicating that their structures are stable.However, when the temperature increases at 362 K, the mutant protein shows slight structural increases in its overall fluctuations, expansions, and loss of native contacts, for example, in the range of 300 to 400 ns, compared with the wild-type protein (Figure 1b,e,h).These changes become more pronounced as the temperature increases to 400 K.It is observed that after 100 ns the mutant protein displays larger global fluctuations, undergoing drastic structural expansion and compaction states with significant losses of native contacts (Figure 1c,f,i).In particular, around 490 ns the BstHPrm protein shows larger fluctuation and structural expansion than the BstHPr protein.an alanine residue, affecting the aforementioned triad.Moreover, the effects of temperature on the mutant structure of the HPr protein were analyzed and compared with those obtained from the thermophilic wild-type (BstHPr) and mesophilic (BsHPr) structures.

Structural Behavior
Figure 1 shows the root mean square deviation (RMSD), radius of gyration (Rg), and fraction of native contact (Q) trajectories for BstHPr and BstHPrm proteins at 298, 362, and 400 K.Both proteins exhibit almost the same structural behavior during the simulation time at 298 K (Figure 1a,d,g), indicating that their structures are stable.However, when the temperature increases at 362 K, the mutant protein shows slight structural increases in its overall fluctuations, expansions, and loss of native contacts, for example, in the range of 300 to 400 ns, compared with the wild-type protein (Figure 1b,e,h).These changes become more pronounced as the temperature increases to 400 K.It is observed that after 100 ns the mutant protein displays larger global fluctuations, undergoing drastic structural expansion and compaction states with significant losses of native contacts (Figure 1c,f,i).In particular, around 490 ns the BstHPrm protein shows larger fluctuation and structural expansion than the BstHPr protein.Figure 2 shows the average values (avg) and their corresponding standard deviation (SD) of the RMSD (Figure 2a), Rg (Figure 2b), and Q (Figure 2c) of the BsHPr, BstHPr, and BstHPrm proteins.It should be noted that the BsHPr tag refers to the HPr protein from the mesophilic organism Bacillus subtilis, whose results are reported in reference [20].In this figure, it is observed that the mutant protein presents almost the same structural behavior as the mesophilic protein.In particular, (a) the average values of RMSD and Q are practically identical at a temperature higher than 362 K and (b) the mutant protein remains more compact than the mesophilic one in the unfolded state at 450 K, as measured by the Rg. Figure 2 shows the average values (avg) and their corresponding standard deviation (SD) of the RMSD (Figure 2a), Rg (Figure 2b), and Q (Figure 2c) of the BsHPr, BstHPr, and BstHPrm proteins.It should be noted that the BsHPr tag refers to the HPr protein from the mesophilic organism Bacillus subtilis, whose results are reported in reference [20].In this figure, it is observed that the mutant protein presents almost the same structural behavior as the mesophilic protein.In particular, (a) the average values of RMSD and Q are practically identical at a temperature higher than 362 K and (b) the mutant protein remains more compact than the mesophilic one in the unfolded state at 450 K, as measured by the Rg.S1-S3 of the Supplementary Materials show the corresponding numerical values.Figure S4 shows the boxplots for these analyses to provide a more detail examination of the statistical distribution of the data.
Figure 3 shows the behavior of the β-strand and α-helix native secondary structures (SS) of the three proteins.These proteins contain the same number of secondary structures at 300 K, i.e., the content of β-strand and α-helix is approximately 26 and 34%, respectively; thus, around 40% of the secondary structure is random coil.Moreover, it is observed that with increasing temperature, the BsHPr and BstHPrm proteins lose β-strand and α-helix structures compared to the BstHPr protein.The most affected structures are the α-helices since they show higher losses at lower temperatures than the β-strands, i.e., more α-helix structures are lost from 362 K, while the β-strands lose a higher structural percentage from 400 K onwards.Additionally, the average percentage of β-strands lost from the mutant protein is between those of the wild-type and mesophilic proteins at 400 and 450 K (Figure 3a).On the other hand, the average percentage of the α-helices lost by the mutant protein is almost identical to that of the mesophilic protein at 362 and 400 K (Figure 3b), ensuring that α-helix structures are more affected than the β-strand ones as the temperature increases.S1-S3 of the Supplementary Materials show the corresponding numerical values.Figure S4 shows the boxplots for these analyses to provide a more detail examination of the statistical distribution of the data.
Figure 3 shows the behavior of the β-strand and α-helix native secondary structures (SS) of the three proteins.These proteins contain the same number of secondary structures at 300 K, i.e., the content of β-strand and α-helix is approximately 26 and 34%, respectively; thus, around 40% of the secondary structure is random coil.Moreover, it is observed that with increasing temperature, the BsHPr and BstHPrm proteins lose β-strand and α-helix structures compared to the BstHPr protein.The most affected structures are the α-helices since they show higher losses at lower temperatures than the β-strands, i.e., more α-helix structures are lost from 362 K, while the β-strands lose a higher structural percentage from 400 K onwards.Additionally, the average percentage of β-strands lost from the mutant protein is between those of the wild-type and mesophilic proteins at 400 and 450 K (Figure 3a).On the other hand, the average percentage of the α-helices lost by the mutant protein is almost identical to that of the mesophilic protein at 362 and 400 K (Figure 3b), ensuring that α-helix structures are more affected than the β-strand ones as the temperature increases.
Considering only the β-strands of the mutant protein, the β 4 -strand shows greater structural loss, which is to be expected since the mutation has been made in this structure.As expected, the other structures affected are the β 1 -and β 2 -strands, which contribute with their Glu3 and Glu36 residues to form the Glu3-Lys62-Glu36 salt bridge network.Therefore, the β 3 -strand does not present structural modifications, i.e., the structural loss is similar to that of the wild-type protein.Considering only the β-strands of the mutant protein, the β4-strand shows greater structural loss, which is to be expected since the mutation has been made in this structure.As expected, the other structures affected are the β1and β2-strands, which contribute with their Glu3 and Glu36 residues to form the Glu3-Lys62-Glu36 salt bridge network.Therefore, the β3-strand does not present structural modifications, i.e., the structural loss is similar to that of the wild-type protein.
On the other hand, considering the α-helices of the same protein, the mutation of site 62 causes the α1-helix to be the most affected of these structures, while the α2-helix remains almost unchanged to that of the wild-type protein.The changes produced in the α1-helix may occur because the Glu84-Arg17 salt bridge decreases its frequency of formation, i.e., the mutation causes the weakening of the salt bridge; thus, the secondary structure is destabilized.Details of these facts will be discussed later.The individual behavior of each βstrand and α-helix is provided in Figures S6 and S7, while the secondary structure profiles are given in Figures S8-S10 in the Supplementary Materials.Those profiles were calculated using the "define secondary structures of proteins" (DSSP) algorithm [20].

Molecular Interactions
We first analyzed the hydrogen bonds (HBs) in both proteins, which are classified into two types: (1) bonds that form between residues of the secondary structures (labeled HBpp) and (2) interactions between protein residues with the solvent (labeled HBps).
Figure 4 shows the trajectories of the HBpp and HBps values of the BstHPr and BstH-Prm proteins at the temperatures of 298, 362, and 400 K.Both proteins show similar trends for the two types of hydrogen bonds.Only slight changes in these trends are observed; for example, HBpp decreases in the BstHPrm protein between 100 and 500 ns at 400 K (Figure 4c) and HBps decreases in the BstHPr protein at around 150-500 ns at 400 K (Figure 4f).On the other hand, considering the α-helices of the same protein, the mutation of site 62 causes the α 1 -helix to be the most affected of these structures, while the α 2 -helix remains almost unchanged to that of the wild-type protein.The changes produced in the α 1 -helix may occur because the Glu84-Arg17 salt bridge decreases its frequency of formation, i.e., the mutation causes the weakening of the salt bridge; thus, the secondary structure is destabilized.Details of these facts will be discussed later.The individual behavior of each β-strand and α-helix is provided in Figures S6 and S7, while the secondary structure profiles are given in Figures S8-S10 in the Supplementary Materials.Those profiles were calculated using the "define secondary structures of proteins" (DSSP) algorithm [20].

Molecular Interactions
We first analyzed the hydrogen bonds (HBs) in both proteins, which are classified into two types: (1) bonds that form between residues of the secondary structures (labeled HBpp) and ( 2) interactions between protein residues with the solvent (labeled HBps).
Figure 4 shows the trajectories of the HBpp and HBps values of the BstHPr and BstHPrm proteins at the temperatures of 298, 362, and 400 K.Both proteins show similar trends for the two types of hydrogen bonds.Only slight changes in these trends are observed; for example, HBpp decreases in the BstHPrm protein between 100 and 500 ns at 400 K (Figure 4c) and HBps decreases in the BstHPr protein at around 150-500 ns at 400 K (Figure 4f).
Figure 5 shows the average values and standard deviations of both hydrogen bonds for the BsHPr, BstHPr, and BstHPrm proteins at the five temperatures studied.The HBpp trends of the three proteins are similar, while the HBps trend is similar between BstHPr and BstHPrm proteins.Yet, these values in the mesophilic protein are higher at all temperatures; this occurs because the BsHPr protein possesses a higher number of polar residues that are exposed to the solvent [20].Figure 5 shows the average values and standard deviations of both hydrogen bonds for the BsHPr, BstHPr, and BstHPrm proteins at the five temperatures studied.The HBpp trends of the three proteins are similar, while the HBps trend is similar between BstHPr and BstHPrm proteins.Yet, these values in the mesophilic protein are higher at all temperatures; this occurs because the BsHPr protein possesses a higher number of polar residues that are exposed to the solvent [20].Figure 5 shows the average values and standard deviations of both hydrogen bonds for the BsHPr, BstHPr, and BstHPrm proteins at the five temperatures studied.The HBpp trends of the three proteins are similar, while the HBps trend is similar between BstHPr and BstHPrm proteins.Yet, these values in the mesophilic protein are higher at all temperatures; this occurs because the BsHPr protein possesses a higher number of polar residues that are exposed to the solvent [20].S4 of the Supplementary Materials shows the corresponding numerical values.The boxplots of these analyses are given in Figure S13 of the Supplementary Materials.
Taking as reference the average values of the HBpp and HBps of the BstHPr and BstHPrm protein structures at 298 K, it is possible to calculate the number and percentage of hydrogen bonds lost for each temperature, which is summarized in Tables 2 and 3.These analyses show that both proteins lose approximately the same number and percentage of HBpp and HBps from 298 to 450 K.However, the mutant protein loses more HBpp at 400 K (19.5%), indicating that the interior of its structure is weakened more than that of the wild-type protein at this point.
Secondly, we analyzed the solvent accessible surface area (SASA) behavior of both proteins for two residue types: polar (labeled SASAp) and non-polar (labeled SASAnp). Figure 6 shows the trajectories of the SASAp and SASAnp values of the BstHPr and BstHPrm proteins at 298, 362, and 400 K. From this figure, it can be seen that both trends increase for the two protein structures as the temperature increases.However, the SASAp trend of the BstHPrm protein is higher than that of the BstHPr one, indicating that the mutation causes greater exposure of the polar residues to the solvent.While the SASAnp of the BstHPrm protein is slightly lower than that of the BstHPr protein at 298 K, this behavior changes as the temperature increases, i.e., the SASAnp of the BstHPrm protein is higher compared to that of the BstHPr protein at 400 K. Therefore, at low temperatures (e.g., 298 K) the non-polar residues of the mutant protein are more buried inside its structure; however, these lose stability as the temperature increases (e.g., 400 K), causing greater contact with the solvent in comparison to the wild-type protein.
Considering the average SASAp and SASAnp values of the BstHPr and BstHPrm proteins, one can calculate the amount and percentage of area that increased or decreased in the thermophilic protein due to the mutation, which is shown in Table 4.  Considering the average SASAp and SASAnp values of the BstHPr and BstHPrm proteins, one can calculate the amount and percentage of area that increased or decreased in the thermophilic protein due to the mutation, which is shown in Table 4.The surface area of polar residues increased for all temperatures.It should be noted that the largest increase is observed at 400 K, whereas the surface area of the non-polar residues increases only at 400 K (values denoted with an asterisk), but this value decreases for the other temperatures.
These analyses confirm that the polar residues of the BstHPrm protein are more exposed to the solvent than those of the BstHPr protein.Furthermore, the non-polar residues of the BstHPrm protein are more buried inside the structure at low temperatures, e.g., at 298 K (see Figure 6d), but these residues become more exposed to the solvent as the temperature increases, e.g., at 400 K (see Figure 6f).The surface area of polar residues increased for all temperatures.It should be noted that the largest increase is observed at 400 K, whereas the surface area of the non-polar residues increases only at 400 K (values denoted with an asterisk), but this value decreases for the other temperatures.
These analyses confirm that the polar residues of the BstHPrm protein are more exposed to the solvent than those of the BstHPr protein.Furthermore, the non-polar residues of the BstHPrm protein are more buried inside the structure at low temperatures, e.g., at 298 K (see Figure 6d), but these residues become more exposed to the solvent as the temperature increases, e.g., at 400 K (see Figure 6f).
Figure 7 shows the average values and their corresponding standard deviations of the SASAp and SASAnp of the BsHPr, BstHPr, and BstHPrm proteins at the five temperatures analyzed.The SASAp behaviors of the BsHPr and BstHPrm proteins are similar, while their SASAnp trends are lower than the values of the BstHPr protein at 298 and 333 K.However, as stated in the previous paragraph, from 362 K onwards it is observed that the mutant protein changes its tendency concerning the wild-type protein.In addition, the same behavior is observed in the mesophilic protein.
In addition to calculating the behavior of the hydrophobic core residues through SASA, we performed the analysis of the ILV clusters of the mutant protein at 750 ns at the different temperatures studied.The ILV clusters for both proteins at 298, 362, and 400 K are given in Figure 8.These results show how these clusters change with the temperature.In particular, the mutant protein forms three clusters at 400 K, while the thermophilic protein forms one cluster at the same temperature.This fact provides insight into the stability of the hydrophobic cluster of the BstHPr protein because the point mutation Lys62Ala triggers structural instability, causing the hydrophobic core to be exposed to the solvent, as previously described with the SASA analyses.
Figure 7 shows the average values and their corresponding standard deviations of the SASAp and SASAnp of the BsHPr, BstHPr, and BstHPrm proteins at the five temperatures analyzed.The SASAp behaviors of the BsHPr and BstHPrm proteins are similar, while their SASAnp trends are lower than the values of the BstHPr protein at 298 and 333 K.However, as stated in the previous paragraph, from 362 K onwards it is observed that the mutant protein changes its tendency concerning the wild-type protein.In addition, the same behavior is observed in the mesophilic protein.S5 of the Supplementary Materials shows the corresponding numerical values.The boxplots of these analyses are given in Figure S16 of the Supplementary Materials.
In addition to calculating the behavior of the hydrophobic core residues through SASA, we performed the analysis of the ILV clusters of the mutant protein at 750 ns at the different temperatures studied.The ILV clusters for both proteins at 298, 362, and 400 K are given in Figure 8.These results show how these clusters change with the temperature.In particular, the mutant protein forms three clusters at 400 K, while the thermophilic protein forms one cluster at the same temperature.This fact provides insight into the stability of the hydrophobic cluster of the BstHPr protein because the point mutation Lys62Ala triggers structural instability, causing the hydrophobic core to be exposed to the solvent, as previously described with the SASA analyses.S5 of the Supplementary Materials shows the corresponding numerical values.The boxplots of these analyses are given in Figure S16 of the Supplementary Materials.
Additionally, we analyzed the formation of ionic pairs in both proteins.In our previous work [20], in which the thermal stability of the BstHPr protein was analyzed, the formation of five salt bridges (SBs) was observed: Asp79-Lys83, Glu84-Arg17, Asp11-Lys57, Glu3-Lys62, and Glu36-Lys62, where the latter two ionic pairs form the Glu3-Lys62-Glu36 triad.The residues of this triad of Glu3, Lys62, and Glu36 are located on the β 1 -, β 4 -, and β 2 -strands, respectively.Thus, the structural and physicochemical nature of the β 4 -strand is expected to be affected, since the Lys62 residue has been mutated by Ala62.Except for the Asp79-Lys83 salt bridge, which is of the intra-molecular type, the other SBs are formed between residues that are in different secondary structures, i.e., they are of the inter-molecular type.Table 5 shows these values.On the other hand, in this work, in which we analyzed the thermal stability of the BstHPr protein but performed the point mutation Lys62Ala, the formation of four salt bridges was observed in the BstHPrm protein: Asp79-Lys83, Glu84-Arg17, Asp11-Lys57, and Glu32-Lys45.The SB Glu32-Lys45 is the only new ionic interaction formed in the mutant protein.Table 6 shows the average frequency of each SB from the three simulations of the BstHPrm protein for the temperatures studied.(d-f) BstHPrm proteins at 298, 362, and 400 K. β-strands, α-helices, and random coils are colored in red, blue, and yellow, respectively.ILV clusters were calculated using the ProteinTools tool (https://proteintools.uni-bayreuth.de/(accessed on 24 March 2024)) [21] and are shown in green, cyan, and magenta.Each snapshot was made using PyMOL version 2.5.2 software for simulation 1 at 750 ns.The initial and final residues of both proteins are shown in gray.Details of these clusters can be found in Figure S17 of the Supplementary Materials.
Additionally, we analyzed the formation of ionic pairs in both proteins.In our previous work [20], in which the thermal stability of the BstHPr protein was analyzed, the formation of five salt bridges (SBs) was observed: Asp79-Lys83, Glu84-Arg17, Asp11-Lys57, Glu3-Lys62, and Glu36-Lys62, where the latter two ionic pairs form the Glu3-Lys62-Glu36 triad.The residues of this triad of Glu3, Lys62, and Glu36 are located on the β1-, β4-, and β2-strands, respectively.Thus, the structural and physicochemical nature of the β4strand is expected to be affected, since the Lys62 residue has been mutated by Ala62.Except for the Asp79-Lys83 salt bridge, which is of the intra-molecular type, the other SBs are formed between residues that are in different secondary structures, i.e., they are of the inter-molecular type.Table 5 shows these values.Considering the values in both Tables 5 and 6, it can be seen that the mutation causes the following changes in these molecular interactions with increasing temperature: (a) The Asp79-Lys83 salt bridge, located on the α 3 -helix, does not undergo meaningful changes; only its frequency is slightly higher at 333 K. (b) The Glu84-Arg17 salt bridge, formed between the α 3 -helix and α 1 -helix, undergoes significant changes, since it presents lower frequencies at all temperatures, causing the secondary structures to be less stable.
(c) The Asp11-Lys57 salt bridge, formed between two loops (that are in the β 1 -strand/α 1helix and the α 2 -helix/β 4 -strand structures), increases its frequency at 298 and 333 K but decreases for the other temperatures.In particular, it decreases drastically at 400 K.(d) The Glu32-Lys45 salt bridge is formed between the β 2 -strand and the first residue of the loop after the β 3 -strand.This interaction is very weak, as its frequency is less than 0.3 at 298 and 333 K.It reaches values slightly higher than 0.3 at 362 and 400 K, but it is almost lost at 450 K.
In addition to these frequency measurements, we calculated the distances of the salt bridges for both proteins.These values can be found in Tables S6 and S7 in the Supplementary Materials.The structural comparison between the residues of the Glu3-Lys62-Glu36 triad in BstHPr and the corresponding residues of the BstHPrm and BsHPr proteins are given in Figure S18 of the Supplementary Materials.
A salient feature to highlight is that due to the mutation, different regions of the protein lose structural strength.For example, the SB Glu84-Arg17 is weakened, as it presents lower frequencies compared to those observed in the wild-type BstHPr protein, indicating that the Lys62 residue forming the Glu3-Lys62-Glu36 triad is crucial in the structural stability of the protein.
Lastly, we performed a brief analysis of the flexibility of the active site in both proteins.The main function of the HPr proteins is to transfer a phosphate group within the PTS system in bacteria [12,16,19].The flexibility of two residues belonging to the active center of the HPr protein, His15 and Arg17, is essential for this transport.Depending on the spatial relationship between these residues, the HPr protein can adopt two main conformations: open state (OS) and closed state (CS).The OS is considered when the distance between the N δ1 of His15 and N of Arg17 is d Nδ1-N > 7.5 Å, and the CS when this distance is d Nδ1-N < 4.5 Å.It has been proposed that dynamic conversion between the OS and the CS is essential for the phosphorylation of this protein and its effect as a transcriptional regulator [19].Through the simulations performed in this study, it was possible to observe that BstHPr and BstHPrm proteins adopt both states frequently during the MD trajectory; hence, the ratio of OS/CS was calculated according to the total number of states achieved, as indicated in Figure 9. is essential for the phosphorylation of this protein and its effect as a transcriptional regulator [19].Through the simulations performed in this study, it was possible to observe that BstHPr and BstHPrm proteins adopt both states frequently during the MD trajectory; hence, the ratio of OS/CS was calculated according to the total number of states achieved, as indicated in Figure 9.

Discussion
BstHPr and BstHPrm proteins lose their structural hierarchy with increasing temperature, as the dominant arrangements at high temperatures no longer correspond to αhelices and β-strands but rather to less ordered structures such as loops and turns.The αhelix and β-strand secondary structures are stable at 298 K and almost all are lost at 450 K

Discussion
BstHPr and BstHPrm proteins lose their structural hierarchy with increasing temperature, as the dominant arrangements at high temperatures no longer correspond to α-helices and β-strands but rather to less ordered structures such as loops and turns.The α-helix and β-strand secondary structures are stable at 298 K and almost all are lost at 450 K (Figure 3), indicating that both proteins are unfolded.However, it can be claimed that the most structurally stable protein throughout the simulations is the wild-type BstHPr variant.This conclusion has been reached because, when comparing the secondary structures of the two proteins, the β 4 -, β 1 -, and β 2 -strands and the α 1 -and α 3 -helices in the BstHPrm protein are more destabilized up to 400 K.
RMSD analysis shows that both proteins maintain their structural stability and almost similar behavior between 298 and 362 K (Figures 1a,b and 2a), which agrees with the secondary structure profiles.At the temperature of 400 K, there is a difference in fluctuations between the proteins, with the BstHPrm protein exhibiting larger structural global fluctuations than the BstHPr one (Figure 1c).At the temperature of 450 K, drastic fluctuations are observed in the trajectories of the three simulations for both proteins (refer to Figure S1).This indicates a noticeable structural difference compared to the reference or native structures at 298 K, showing that both proteins have reached the unfolded state.
The radius of gyration of both proteins increases its average value and standard deviation between 298 and 450 K (Figure 2b), indicating that the expansion and compaction processes of the protein structures are meaningful as the temperature increases.In particular, this fact is noticeable at 450 K, since at this point structural unfolding occurs (as shown in Figure S2).Furthermore, the mutant variant experiences greater structural contractions and expansions than the wild-type variant at 400 K (Figure 1f), which is consistent with the RMSD behavior and secondary structure conformations.
The fraction of native contacts also decreases in both proteins (Figure 2c) from 298 to 450 K, i.e., topological interactions are lost regarding the initial structures (conformations at t = 0 ns).This statement is consistent with the aforementioned analyses, specifically the analysis of secondary structures (Figure 3), which indicates that the proteins lose a large number of their α-helix and β-strand structures at 450 K.However, it is clear that out of the two proteins, and as expected, the native contacts of the BstHPrm protein decay more than those of the BstHPr protein, which is illustrated in Figures 1h,i and 2c.In particular, these contacts fall more in the BstHPrm protein at 400 K.
It is known that each hydrogen bond in proteins contributes energetically on average 1 kcal/mol [22]; therefore, increasing the number of HBs improves the thermostability of these biomolecules.However, in this study, HBpp and HBps interactions decreased for both proteins with increasing temperature (Figures 4 and 5).This occurs because the energy yielded to the system causes an increase in the kinetic energy of the particles, promoting the rupture of HB.If the number of HB decreases, then secondary structures are drastically lost, global fluctuations and compaction/expansion states increase, and more native contacts are lost, as described in the previous paragraphs.
In this way, comparing the behavior of the two proteins, it is observed that the reductions in HBpp and HBps are almost similar up to 362 K. Nonetheless, when the temperature increases to 400 K, more HBpp of the BstHPrm protein is lost than in the BstHPr protein, e.g., in the range from 298 to 400 K 19.5 and 15.1% is lost, respectively, whereas between 362 and 400 K 14.0 and 9.2% is lost, respectively.The opposite case arises with HBps since it increases in the mutant variant and decreases in the wild-type protein, e.g., HBps decreases by 1.6% in the BstHPr and increases by 0.5% in the BstHPrm between 362 and 400 K.These facts occur because the mutant protein is destabilized and unfolds more at this temperature (400 K).The residues buried in the structure are exposed to the solvent and the HBpp interactions are broken; consequently, they tend to form more hydrogen bonds with the solvent.This agrees with the structural analyses of RMSD, Rg, Q, and SS profiles, since at 400 K the most meaningful changes between the two proteins are observed.Additionally, from the native to the unfolded state (298-450 K), HBpp decreases by 33.3 and 33.9% for BstHPr and BstHPrm proteins, respectively, and HBps shows losses of 6.6% for BstHPr and 7.0% for BstHPrm, i.e., both proteins lose almost the same number of HBs, yielding structural similarities at 450 K.
Hydrophobic interactions and packing of the hydrophobic core of proteins are calculated indirectly through SASA.This is achieved by calculating the area of the proteins in contact with the solvent.The measurement is made on the total or part of the residues exposed to the solvent and can be distributed into polar and non-polar residue contacts.SASA is another useful parameter that indicates the loss of stability and unfolding of proteins, since hydrophobic zones are exposed when increasing the temperature, which usually interact with each other through non-polar residues that hardly have contact with the solvent.Consequently, in this work, the hydrophobic zones remain stable in both BstHPr and BstHPrm proteins at 298 K; nevertheless, as the temperature increases up to 450 K, these zones are exposed to the medium, i.e., the interaction between the non-polar residues decreases, leading to the instability of the hydrophobic core.The hydrophobic core reaches its maximum exposure for the two proteins at 450 K (Figure 7).
The Lys62Ala mutation causes the polar area to increase by an average of 8.7% at 298 K, keeping this trend almost constant up to 450 K, where the average difference is equal to 5.9%; that is, although the temperature increases, the difference in polar area between the proteins does not change significantly.Conversely, the mutation causes the non-polar area to decrease slightly up to 362 K; however, this area increases at 400 K by 2.3%, i.e., the hydrophobic core of the BstHPrm protein is more exposed to the solvent than that of the BstHPr protein at this temperature.In other words, the hydrophobic interactions of BstHPrm are weaker at this point, losing more structural stability.In addition, SASAnp increases by an average of 14.3% for the BstHPrm protein and by an average of 8.3% for the BstHPr protein in the range of 298 to 400 K.The proteins reach their maximum nonpolar area exposure at 450 K; as such, the area increases by 39.0 and 39.7% for BstHPr and BstHPrm, respectively, between 298 and 450 K.In addition to these results, the ILV cluster analysis shows that the hydrophobic core of the mutant protein loses stability at 400 K as the number of these clusters increases.These facts are consistent with the structural and hydrogen bond analyses, as discussed in the aforementioned paragraphs.
Salt bridges play an important role in the stability of proteins, especially when they are exposed to high temperatures.In this work, the mutant protein conserves the salt bridges Asp79-Lys83, Glu84-Arg17, and Asp11-Lys57 from the wild-type protein.The Asp79-Lys83 salt bridge is the most stable one to temperature changes during MD simulations since it does not show any significant variations in its average formation frequency compared to the BstHPr protein, and according to our frequency criterion, its value is greater than 0.3 up to 450 K.The stability of this SB is due to its intramolecular nature, i.e., it is located on the α 3 -helix structure, which is slightly affected by temperature (see Figure S7c).On the other hand, the Glu84-Arg17 salt bridge is the least stable to temperature changes, as its average formation frequency decreases largely upon mutation (Table 6).The mutation Lys62Ala causes the formation of the Glu32-Lys45 salt bridge in the BstHPrm variant.This interaction is unlikely to contribute to its stability since the formation frequency is less than 0.3 at 298 K; that is, it presents an average frequency of 0.246.However, it may provide some structural stability when increasing to 362 and 400 K, since the average frequency increases to 0.334 and 0.351, respectively.The β 1 -and β 4 -strands are the secondary structures most affected by the mutation, as they lose more structure when the temperature reaches 400 K (see Figure S6a,d).Additionally, more β-strand fragments are observed in the BstHPr protein than in the BstHPrm one at 450 K (see Figure S6 and bottom panels of Figures S9 and S10), indicating that the mutation affects the structural behavior of these structures up to the unfolded state.Analyzing these facts and comparing them with the results of the parameters RMSD, Rg, Q, SS, HB, and SASA, it is gathered that the salt bridge network formed by the Glu3-Lys62-Glu36 triad provides greater thermal stability to the BstHPr protein.These parameters are less affected in the wild-type protein than in the mutant protein up to 400 K.
An interesting fact is that when comparing the behavior of the structural (RMSD, Rg, Q, and SS) and molecular interaction (HB and SASA) parameters of the BstHPr and BstHPrm proteins with their mesophilic BsHPr counterpart, the mutant protein shows identical trends to those of the mesophilic protein up to 400 K (the exception to this is the behavior of HBps), indicating that the mutation possibly causes a decrease in the melting point (Tm) of the thermophilic protein.Additionally, the net charge in the mutant protein was calculated to be −2 at 298 K, while the thermophilic (wild-type) and mesophilic proteins have net charges equal to −1 and −4, respectively.Therefore, the mutation increases the charge repulsion potential, destabilizing the structure.
It has not been possible to find experimental data on the formation of this Glu3-Lys62-Glu36 triad or on the site-directed point mutation Lys62Ala in the BstHPr protein.However, it has been reported that the HPr protein from the thermophilic organism Thermoanaerobacter tengcongensis (TtHPr) forms the Glu3-Lys62 salt bridge, which might compact the structure between β1 and β4 strands, contributing to its thermal stability.Therefore, the site-directed Lys62Thr was mutated in the TtHPr variant.The authors concluded that the measurements by CD spectra provide firm evidence that the Glu3-Lys62 salt bridge is key in the thermostability of the TtHPr protein [16].
In summary, the Glu3-Lys62-Glu36 triad, formed by the Glu3-Lys62 and Glu36-Lys62 salt bridges, plays a crucial role in the stability of the hydrophobic core of the thermophilic BstHPr protein, avoiding its exposure to the solvent and thereby preventing the breaking of buried hydrogen bonds and decreasing the electrostatic surface potential.Thus, this triad functions as a "natural molecular staple" that keeps the structure packaged and stable against energetic changes caused by the thermal increase.
Finally, the OS/CS ratio as a function of temperature may explain the functional mechanism of the HPr proteins regarding their structural stability.Although these proteins have a common folding (open-faced β-sandwich type), their geometry in the vicinity of the active site has marked differences.For example, the crystallographic structure of BsHPr exhibits a distance of d Nδ1-N = 3.9 Å, adopting the CS state, while the BstHPr crystal exhibits the OS state, with a distance equal to d Nδ1-N = 8.8 Å [19].Consequently, these differences complicate the comparison between two different proteins, even when performing MD simulations.These structural changes should directly impact the catalytic activity of HPr proteins.However, direct measurement of activity requires additional MD protocols from those used in this investigation.
On the other hand, the only difference between BstHPr and BstHPrm is the Lys62Ala point mutation; thus, a direct comparison can be made between these proteins.Figure 9 shows that the ratio of OS/CS approaches 1 at low temperatures.However, this equilibrium is considerably affected due to unfolding effects at elevated temperatures.There is a correlation between the structural changes in the region of the Lys62 residue and the behavior of the active center, since the mutant protein exhibits more OS than that of the wild-type one at 400 K. Obviously, both proteins reach OS/CS equilibrium at 450 K, as they are fully unfolded.
strands and three α-helices, resulting in a β1α1β2β3α2β4α3 arrangement (Figu [14,16,23].The BstHPr crystallographic structure was protonated using the PropK gram [24] to achieve neutral pH, and the point mutation Lys62Ala was created PyMOL version 2.5.2 (https://pymol.org(accessed on 24 March 2024)).To confirm t mutation is a destabilizing site in the BstHPr protein, we carried out virtual mu using 6 different methods [25].Details and comments are presented in Table S8 Supplementary Materials.[26,27] and the AMBER99SB force field [28] were used to perform the MD simulations.Wild-type BstHPr and mutant BstHPrm proteins were solvated using the SPCE water model [29] inside a dodecahedral box with periodic boundary conditions.The distance between the proteins and the edge of the box was equal to 1.0 nm.The systems' total charge was neutralized at pH = 7.0 by adding sufficient counterions employing the "genion" function of GROMACS [26].The systems BstHPr and BstHPrm consisted of 1290 and 1278 protein atoms, 15,084 and 13,131 solvent atoms, 1 and 0 Na + ions, and 2 and 2 Cl -ions, yielding a total of 16,377 and 14,411 atoms, respectively.

GROMACS v2020.3 program
The particle mesh Ewald (PME) algorithm was used to evaluate the long-range electrostatic interactions [30].Specifically, one Coulomb cutoff equal to 1.0 nm and Fourier spacing of 0.16 nm with a 4th-order cubic interpolation were used.In addition, one cutoff of 1.0 nm for the van der Waals interactions was set.The bond lengths were constrained using the LINCS method [31] with a 2 ps integration step.The "steepest descent" function [32] was used to minimize the energy for both proteins, and the processes converged at 624 and 446 steps (−2.68 × 10 5 and −2.32 × 10 5 kJ/mol) for the BstHPr and BstHPrm proteins, respectively.The velocity-rescaling algorithm was used to equilibrate the NVT ensemble [33]; the temperatures were fixed at 298, 333, 362, 400, and 450 K.The systems were simulated at temperatures equal to 400 K and 450 K to accelerate global and atomic fluctuations of structures to activate protein unfolding, since the experimental and simulated timescales for this process are very different [20].The NPT ensemble was equilibrated at 1.0 bar using the Parrinello-Rahman barostat [34].Both equilibrations were achieved during 100 ps.
In this work, we performed fifteen simulations of the mutant BstHPrm protein, which were then compared with the simulations obtained of the wild-type BstHPr and mesophilic BsHPr proteins given in reference [20].This way, the same time and number of configurations were simulated, i.e., 1 µs of simulation and 10,001 configurations.Moreover, we repeated simulations at each temperature in triplicate.

Simulation Analysis
As stated in previous work [20], here we carried out the following analyses: (a) root mean square deviation, (b) radius of gyration, (c) fraction of native contacts, (d) secondary structure profiles, (e) hydrogen bonds, (f) solvent accessible surface area, (g) cluster ILV, and (h) salt bridges.However, we briefly describe some noteworthy keys of these analyses.Eaton model [35].For this calculation, the first conformation of the simulations is defined as the native structure (t = 0 ns).The total number of contacts in the native structure is taken as Q = 1, and from this reference, the contacts for the remaining conformations of the trajectory are obtained.(d) Secondary structure profiles: SS assignment is performed using the define secondary structure of proteins (DSSP) algorithm [36].This algorithm considers 8 types of SS: α-helix, π-helix, 3 10 -helix, β-strand, β-bridge, random coil, bend, and turn.After this calculation, Micsonai et al. classified these structures in three different groups [37], that is, in the α-helix SS the 3 helix structures (α-helix, π-helix, and 3 10 -helix) are included, in the β-strand SS only the β-strand is considered, and in the random coil the remaining structures (β-bridge, random coil, bend, and turn) are included.Micsonai et al. proposed this classification from protein structure data and their respective circular dichroism spectra.Therefore, the set of the three classifications (α-helix, β-strand, and random coil) is considered 100% of the secondary structure.(e) Hydrogen bonds: For this calculation, the distance r and the angle θ between the mass centers of the acceptor (A) and donor (D) atoms of the proton (H) are considered (r AD ≤ 3.5 Å and θ AD ≤ 30 • ).(f) Solvent accessible surface area: This parameter is determined using the Lee and Richards approximation: one solvent sphere with a radius of 1.4 Å is used [38].(g) Cluster ILV: The contacts of structural units (CSU) algorithm is used to find the groupings of isoleucine, leucine, and valine residues within proteins [39].This methodology analyzes atoms as spheres with a van der Waals radius.The contact of two atoms, A and B, is considered, i.e., a test sphere on the surface of atom A must overlap at least 10 Å with the surface of the sphere of atom B. If this contact occurs between the atoms of the residues Ile, Leu, and Val, they are considered part of a cluster.Therefore, different ILV clusters can be expected in the proteins.(h) Salt bridges: The Barlow and Thornton criterion is taken to measure SB formation [40], i.e., r SB ≤ 0.4 nm.In addition, the ionic pairs were calculated using the GetContacts program (https://getcontacts.github.io/(accessed on 22 February 2024)), taking as a criterion of formation that the average frequency must be equal to 0.3 during the trajectories of the three replicas for each temperature.

Conclusions
By comparing and analyzing the structural and molecular interaction parameters between both proteins (BstHPr and BstHPrm), the salt bridge network formed with the Glu3-Ala62-Glu36 triad is essential for maintaining the thermal stability of the thermophilic protein, since without these interactions two phenomena take place, namely: (a) Global fluctuations increase, compaction/expansion processes increase, topological native contacts decrease, ordered secondary structures are lost, and disordered structures increase; (b) Buried hydrogen bonds decrease and those formed with the solvent increase, while the non-polar residues are more exposed to the solvent, causing a loss of hydrophobic interactions.

Figure 1 .
Figure 1.Time evolution of the: (a-c) RMSD; (d-f) Rg; and (g-i) Q for the BstHPr (red line) and BstHPrm (black line) proteins at 298 K, 362 K, and 400 K.The solid lines and shaded areas represent the average values and their standard deviation values from the three replicas of MD simulations, respectively.Details on these trajectories and the average values for the temperatures analyzed are given in Figures S1-S3 in the Supplementary Materials.

Figure 1 .
Figure 1.Time evolution of the: (a-c) RMSD; (d-f) Rg; and (g-i) Q for the BstHPr (red line) and BstHPrm (black line) proteins at 298 K, 362 K, and 400 K.The solid lines and shaded areas represent the average values and their standard deviation values from the three replicas of MD simulations, respectively.Details on these trajectories and the average values for the temperatures analyzed are given in Figures S1-S3 in the Supplementary Materials.

Figure 2 .
Figure 2. Average values and standard deviations of the (a) RMSD, (b) Rg, and (c) Q for the BsHPr, BstHPr, and BstHPrm proteins at 298, 333, 362, 400, and 450 K.The values have been calculated from the three replicas of MD simulations.Symbols represent the avg values and the SD values are indicated with bars.The dashed lines are only guides for the eye.TablesS1-S3of the Supplementary Materials show the corresponding numerical values.FigureS4shows the boxplots for these analyses to provide a more detail examination of the statistical distribution of the data.

Figure 2 .
Figure 2. Average values and standard deviations of the (a) RMSD, (b) Rg, and (c) Q for the BsHPr, BstHPr, and BstHPrm proteins at 298, 333, 362, 400, and 450 K.The values have been calculated from the three replicas of MD simulations.Symbols represent the avg values and the SD values are indicated with bars.The dashed lines are only guides for the eye.TablesS1-S3of the Supplementary Materials show the corresponding numerical values.FigureS4shows the boxplots for these analyses to provide a more detail examination of the statistical distribution of the data.

Figure 3 .
Figure 3. Average percentage values and standard deviations of the native structures: (a) β-strands and (b) α-helices for the BsHPr, BstHPr, and BstHPrm proteins at 298, 333, 362, 400, and 450 K.The values have been calculated from the three replicas of MD simulations.Symbols represent the avg values and the SD values are indicated with bars.Dashed lines are only guides for the eye.The boxplots of these analyses are given in Figure S5 of the Supplementary Materials.

Figure 3 .
Figure 3. Average percentage values and standard deviations of the native structures: (a) β-strands and (b) α-helices for the BsHPr, BstHPr, and BstHPrm proteins at 298, 333, 362, 400, and 450 K.The values have been calculated from the three replicas of MD simulations.Symbols represent the avg values and the SD values are indicated with bars.Dashed lines are only guides for the eye.The boxplots of these analyses are given in Figure S5 of the Supplementary Materials.

Figure 4 .
Figure 4. Time evolution of hydrogen bonds for the BstHPr (red line) and BstHPrm (black line) proteins at (a) 298 K, (b) 362 K, and (c) 400 K for HBpp and (d) 298 K, (e) 362 K, and (f) 400 K for HBps.The solid lines and shaded areas represent the average values and their standard deviation values from the three replicas of MD simulations, respectively.Details on these trajectories and the average values for the temperatures analyzed are given in Figures S11 and S12 of the Supplementary Materials.

Figure 4 .
Figure 4. Time evolution of hydrogen bonds for the BstHPr (red line) and BstHPrm (black line) proteins at (a) 298 K, (b) 362 K, and (c) 400 K for HBpp and (d) 298 K, (e) 362 K, and (f) 400 K for HBps.The solid lines and shaded areas represent the average values and their standard deviation values from the three replicas of MD simulations, respectively.Details on these trajectories and the average values for the temperatures analyzed are given in Figures S11 and S12 of the Supplementary Materials.

Figure 4 .
Figure 4. Time evolution of hydrogen bonds for the BstHPr (red line) and BstHPrm (black line) proteins at (a) 298 K, (b) 362 K, and (c) 400 K for HBpp and (d) 298 K, (e) 362 K, and (f) 400 K for HBps.The solid lines and shaded areas represent the average values and their standard deviation values from the three replicas of MD simulations, respectively.Details on these trajectories and the average values for the temperatures analyzed are given in Figures S11 and S12 of the Supplementary Materials.

Figure 5 .
Figure 5. Average values and standard deviation of the (a) HBpp and (b) HBps for the BsHPr, BstHPr, and BstHPrm proteins at 298, 333, 362, 400, and 450 K.The values have been calculated from the three replicas of the MD simulations.Symbols represent the avg values, and the SD values are indicated with bars.The dashed lines are only guides for the eye.TableS4of the Supplementary Materials shows the corresponding numerical values.The boxplots of these analyses are given in FigureS13of the Supplementary Materials.

Figure 6 .
Figure 6.Time evolution of solvent accessible surface area for the BstHPr (red line) and BstHPrm (black line) proteins at (a) 298 K, (b) 362 K, and (c) 400 K for SASAp and (d) 298 K, (e) 362 K, and (f) 400 K for SASAnp.The solid lines and shaded areas represent the average values and their standard deviation from the three replicas of MD simulations, respectively.Details on these trajectories and the average values for the temperatures analyzed are given in Figures S14 and S15 of the Supplementary Materials.

Figure 6 .
Figure 6.Time evolution of solvent accessible surface area for the BstHPr (red line) and BstHPrm (black line) proteins at (a) 298 K, (b) 362 K, and (c) 400 K for SASAp and (d) 298 K, (e) 362 K, and (f) 400 K for SASAnp.The solid lines and shaded areas represent the average values and their standard deviation from the three replicas of MD simulations, respectively.Details on these trajectories and the average values for the temperatures analyzed are given in Figures S14 and S15 of the Supplementary Materials.

Figure 7 .
Figure 7. Average values and standard deviation of the (a) SASAp and (b) SASAnp for the BsHPr, BstHPr, and BstHPrm proteins at 298, 333, 362, 400, and 450 K.The values have been calculated from the three replicas of the MD simulations.Symbols represent the avg values, and the SD values are indicated with bars.The dashed lines are only guides for the eye.TableS5of the Supplementary Materials shows the corresponding numerical values.The boxplots of these analyses are given in FigureS16of the Supplementary Materials.

Figure 7 .
Figure 7. Average values and standard deviation of the (a) SASAp and (b) SASAnp for the BsHPr, BstHPr, and BstHPrm proteins at 298, 333, 362, 400, and 450 K.The values have been calculated from the three replicas of the MD simulations.Symbols represent the avg values, and the SD values are indicated with bars.The dashed lines are only guides for the eye.TableS5of the Supplementary Materials shows the corresponding numerical values.The boxplots of these analyses are given in FigureS16of the Supplementary Materials.

Figure 8 .
Figure 8. ILV clusters for: (a-c) BstHPr and(d-f) BstHPrm proteins at 298, 362, and 400 K. β-strands, α-helices, and random coils are colored in red, blue, and yellow, respectively.ILV clusters were calculated using the ProteinTools tool (https://proteintools.uni-bayreuth.de/(accessed on 24 March 2024))[21] and are shown in green, cyan, and magenta.Each snapshot was made using PyMOL version 2.5.2 software for simulation 1 at 750 ns.The initial and final residues of both proteins are shown in gray.Details of these clusters can be found in FigureS17of the Supplementary Materials.

Figure 8 .
Figure 8. ILV clusters for: (a-c) BstHPr and(d-f) BstHPrm proteins at 298, 362, and 400 K. β-strands, α-helices, and random coils are colored in red, blue, and yellow, respectively.ILV clusters were calculated using the ProteinTools tool (https://proteintools.uni-bayreuth.de/(accessed on 24 March 2024))[21] and are shown in green, cyan, and magenta.Each snapshot was made using PyMOL version 2.5.2 software for simulation 1 at 750 ns.The initial and final residues of both proteins are shown in gray.Details of these clusters can be found in FigureS17of the Supplementary Materials.

Figure 9 .
Figure 9. Ratio of OS/CS for BstHPr and BstHPrm proteins at different temperatures of simulation.Symbols represent the avg values, and the SD values are indicated with bars.The dashed lines are only guides for the eye.

Figure 9 .
Figure 9. Ratio of OS/CS for BstHPr and BstHPrm proteins at different temperatures of simulation.Symbols represent the avg values, and the SD values are indicated with bars.The dashed lines are only guides for the eye.

Figure 10 .
Figure 10.Snapshot of the BstHPr protein at t = 0 ns of simulation at 298 K. Structure create PyMOL version 2.5.2 software.The color codes of the SS are the same as in Figure 8.

Figure 10 .
Figure 10.Snapshot of the BstHPr protein at t = 0 ns of simulation at 298 K. Structure created using PyMOL version 2.5.2 software.The color codes of the SS are the same as in Figure 8.
(a) Root mean square deviation: The first structural conformation of the DM simulation is used as the reference structure (t = 0 ns).(b) Radius of gyration: This parameter is calculated from the protein center of mass.(c) Fraction of native contacts: This indicator is determined using the Best-Hummer-

Table 1 .
Thermodynamic data of the BstHPr and BsHPr proteins.

Table 2 .
Approximate number and percentage of HBpp lost for BstHPr and BstHPrm.

Table 3 .
Approximate number and percentage of HBps lost for BstHPr and BstHPrm.

Table 4 .
Increased and decreased area in the BstHPr protein by the mutation.
* These values increase.

Table 4 .
Increased and decreased area in the BstHPr protein by the mutation.

Table 6 .
Average frequencies of salt bridges for the BstHPrm protein.