Next Article in Journal
The Impact of Plant Phytochemicals on the Gut Microbiota of Humans for a Balanced Life
Previous Article in Journal
Impact of Multiple Sclerosis Risk Polymorphism rs7665090 on MANBA Activity, Lysosomal Endocytosis, and Lymphocyte Activation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Structure-Function Relationships in Temperature Effects on Bacterial Luciferases: Nothing Is Perfect

by
Anna A. Deeva
1,*,
Albert E. Lisitsa
1,
Lev A. Sukovatyi
1,
Tatiana N. Melnik
2,
Valentina A. Kratasyuk
1,3 and
Elena V. Nemtseva
1,3
1
Biophysics Department, Siberian Federal University, 660041 Krasnoyarsk, Russia
2
Institute of Protein Research, Russian Academy of Sciences, 142290 Pushchino, Russia
3
Institute of Biophysics, Siberian Branch of Russian Academy of Sciences, 660036 Krasnoyarsk, Russia
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2022, 23(15), 8119; https://doi.org/10.3390/ijms23158119
Submission received: 6 July 2022 / Revised: 20 July 2022 / Accepted: 21 July 2022 / Published: 23 July 2022
(This article belongs to the Section Molecular Biophysics)

Abstract

:
The evaluation of temperature effects on the structure and function of enzymes is necessary to understand the mechanisms underlying their adaptation to a constantly changing environment. In the current study, we investigated the influence of temperature variation on the activity, structural dynamics, thermal inactivation and denaturation of Photobacterium leiognathi and Vibrio harveyi luciferases belonging to different subfamilies, as well as the role of sucrose in maintaining the enzymes functioning and stability. We used the stopped-flow technique, differential scanning calorimetry and molecular dynamics to study the activity, inactivation rate, denaturation and structural features of the enzymes under various temperatures. It was found that P. leiognathi luciferase resembles the properties of cold-adapted enzymes with high activity in a narrow temperature range and slightly lower thermal stability than V. harveyi luciferase, which is less active, but more thermostable. Differences in activity at the studied temperatures can be associated with the peculiarities of the mobile loop conformational changes. The presence of sucrose does not provide an advantage in activity but increases the stability of the enzymes. Differential scanning calorimetry experiments showed that luciferases probably follow different denaturation schemes.

1. Introduction

Living systems have highly sensitive regulatory mechanisms to perceive a wide range of environmental conditions and quickly adapt to even the slightest changes. In particular, microorganisms can tune the salt concentration, density, viscosity and other characteristics of the cell matrix in response to acute stress or during self-regulating processes that enable them to maintain homeostasis [1,2,3,4]. One of the most important parameters that affects microbial metabolism and proliferation is temperature. It is generally accepted to distinguish microorganisms into psychrophiles, mesophiles, and thermophiles according to the range of temperatures which they can inhabit. The accumulation of water-soluble organic compounds either by de novo synthesis or by consumption from the surrounding environment is a common mechanism of osmoregulation, UV-protection, and thermal adaptation among extremophiles and mesophiles as well (including Bacteria, some Eukaryota and a few methanogenic Archaea species) [5]. These compounds, so called extremolytes, are of different chemical structure, comprising saccharides, polyols, heterosides, amino acids, and their derivatives [6]. The adjustments occurring in the cell are aimed at maintaining various parameters including viscosity, which is crucial for diffusion-controlled reactions [3], as well as the structure and functions of cytosol components and mainly proteins [7,8]. Furthermore, proteins themselves could be adapted to temperature shifts that affect their folding, structure, dynamics and functional features [9]. The current work focuses on temperature effects on the properties of bacterial luciferases.
Luciferases are responsible for light emitting reactions in luminous organisms, among which bacteria are the most ubiquitous. There are more than 25 species of bioluminescent bacteria found within five genera: Aliivibrio, Vibrio, and Photobacterium within the Vibrionaceae family, Photorhabdus within the Enterobacteriaceae family, and Shewanella within the Shewanellaceae family [10]. They can be successfully isolated in different seasons and at various depths from the warm coastal and pelagic waters, as well as from the northern seas of the global ocean [11,12]. There are also freshwater and terrestrial luminous bacteria, and the latter group consists of only Photorhabdus species that are parasites infecting human skin lesions or insects with nematodes as the intermediate host [13]. Cold adapted representatives can be found in all known genera except the mentioned Photorhabdus. Recent studies also revealed several dozens of psychrophilic luminous strains, including first luminous Kosakonia cowanii, associated with the digestive system of fishes from the White, Okhotsk, and Bering Seas [14,15]. Thus, luminous bacteria have adapted to a wide range of temperatures, from near zero to about 35 °C, although thermophilic species have not been identified.
Bacterial luciferase is a heterodimeric enzyme that catalyzes the oxidation of reduced flavin mononucleotide (FMNH), the monooxygenation of a long chain aliphatic aldehyde and the reduction of molecular oxygen resulting in the formation of the electronically excited intermediate that emits light [16]. The stability of luciferase to thermal inactivation in vitro could vary depending on the bacterial species from which it was isolated. For example, it was obtained that the remaining activity is higher for Vibrio (Beneckea) harveyi luciferase than for Photobacterium phosphoreum and Aliivibrio (Photobacterium) fischeri [17]. However, V. harveyi luciferase is less thermostable than its close homolog from Vibrio campbellii that differs by only eight amino acid residues (98.8% sequence identity) [18]. Also, Photorhabdus luminescens luciferase is remarkably more thermostable compared to V. harveyi, while their α-subunits share over 85% amino acid sequence identity [19]. A recent phylogenetic analysis of luciferase amino acid sequences showed that they fall into two highly conserved groups with a specific structure of the active site [20]. Each group comprises luciferases with the similar type of decay kinetics in reaction with dodecanal: “fast” type from the Aliivibrio and Photobacterium species, and “slow”—from Vibrio and Photorhabdus [21,22,23]. The majority of known psychrophilic luciferases belong to the group with the fast decay kinetics, however a few species with “slow” luciferase were also isolated from the northern latitudes, for example, Vibrio splendidus [15]. The presence of both types of luciferases among psychrophilic strains suggests that either luciferases acquired the properties of cold-active enzymes, regardless of the structural features of the active center, or the adaptation of luminous bacteria passed through a change in the other cellular components and did not affect the bioluminescent system.
This research aimed to reveal the effect of temperature on the reaction kinetics and structure of two types of bacterial luciferases and the role of extremolytes in maintaining the activity of these enzymes under unfavorable temperature conditions. We studied Vibrio harveyi luciferase characterized by the slow decay kinetics and Photobacterium leiognathi defined as the “fast” enzyme. Sucrose was chosen as a model extremolyte, since its buffer solutions proved to have no inhibitory effect on bacterial luciferase activity in contrast to glycerol solutions [24].
Our data showed that “fast” P. leiognathi luciferase remains active at relatively low temperatures, while the temperature range of “slow” V. harveyi luciferase functioning is shifted to higher temperatures. Molecular dynamics simulations revealed structural motifs, which could be critical for the activity maintenance at the corresponding temperatures: a mobile loop forming the active center and a hairpin loop located near the α/β-interface. We also found that the addition of sucrose significantly reduces the mobile loop fluctuations in the P. leiognathi structure at high temperature. On the contrary, sucrose was less effective in stabilizing both flexible motifs of V. harveyi luciferase. The analysis of thermal inactivation rate and denaturation thermodynamics of the luciferases revealed that V. harveyi enzyme is more thermostable. Sucrose was found to have a stabilizing effect against thermal inactivation and the denaturation of both luciferases. Thus, we carried out a detailed comparative study of the temperature effects on two types of bacterial luciferase that have never been performed before.

2. Results

2.1. Activity of V. harveyi and P. leiognathi Luciferases at Different Temperatures

The kinetics of light emission in reactions catalyzed by V. harveyi or P. leiognathi luciferase was studied in the buffer and 30 (w/w)% sucrose at temperatures in the range of 5–45 °C. After mixing with substrates, bacterial luciferase can only make a single turnover due to the quick autoxidation of unbound FMNH. Hence, the reaction kinetics is flash-like with a rapid increase succeeded by a relatively slow decay of light intensity (Figure 1). The following empirical kinetic parameters are estimated: the peak intensity (Imax), the total quantum yield (Q*), and the decay constant (kdecay). Their dependences on temperature are shown in Figure 2 and Figure S1 (in the Supplementary Material).
Comparing luciferases in the buffer, one can see that according to Imax V. harveyi luciferase demonstrates its highest level of activity over a wider temperature range than P. leiognathi enzyme (20–35 °C vs. 20–25 °C, Figure 2a,b). However, the amount of emitted light in a single turnover, Q*, in the range 10–25 °C is about twofold less for V. harveyi enzyme than for P. leiognathi one, but becomes almost equal in value for both luciferases at higher temperatures (Figure S1a,b). The decay rate, kdecay, increases under higher temperatures (≥35 °C) more sharply for P. leiognathi luciferase (Figure 2d). Thus, in spite of high homology, P. leiognathi luciferase exhibits the properties of a more thermolabile enzyme than the V. harveyi one.
To mimic a viscous environment and the extremolyte effects on bioluminescent reactions from two bacterial species, we studied reaction kinetics in sucrose solution (Figure 1). A concentration of 30% was used, which corresponded to 0.99 M of sucrose and gave more than a threefold higher viscosity as compared with the buffer. From the kinetic curves in Figure 1, one can see that the reaction noticeably slows down in a more viscous medium for both luciferases, which is due to the diffusion control of some reaction steps [24]. The change of kinetic parameters in sucrose solution over the entire temperature range is shown in Figure 2 and Figure S1. We observed that in the presence of sucrose the peak intensity decreases at low temperatures (for V. harveyi luciferase—in the range 5–20 °C, for P. leiognathi one—in the range 5–15 °C) (Figure 2a,b), while the effect on the decay constant appears only at high temperatures (Figure 2c,d). All of this leads to a decrease of the reaction quantum yield in sucrose solution at lower temperatures for both luciferases (Figure S1). Interestingly, at temperatures >30 °C in the presence of sucrose, kdecay is higher for V. harveyi luciferase and lower for the P. leiognathi one, as compared with the buffer.
The observed temperature effects on the activity of the enzymes indicate that (i) P. leiognathi luciferase is more sensitive to temperature change than the V. harveyi enzyme and (ii) viscous sucrose solution gives no advantages for bacterial luciferases’ functioning either under “hot” (35–45 °C) or in “cold” (5–15 °C) conditions. These facts will be discussed in Section 3.

2.2. Thermal Inactivation of V. harveyi and P. leiognathi Luciferases

We measured the remaining activity of V. harveyi and P. leiognathi luciferases at 20 °C after incubating the enzymes for different time intervals at various temperatures in the absence and presence of sucrose. The remaining activity at a certain time was expressed as percent of activity in the initial moment (i.e., without incubation), which is described in detail in Section 4.2. The time course of thermal inactivation of the enzymes at different temperatures is shown in Figure 3.
It is worth noting that for both enzymes we observed an induction phase (from 0.5 to 8 min, depending on temperature) (Figure 3). Incubation during this time resulted in the remaining activity of about 100%, indicating that the most thermally inactivated molecules recovered after the immediate cooling of incubated enzymes. The longer heat treatment resulted in a consistent decrease of the remaining activity, which was well described by a single exponential function. The corresponding rate constants of thermal inactivation of the bacterial luciferases are shown in Table 1. In a few cases an additional slow phase of inactivation was observed for both enzymes on further heating, namely at 50 °C (Figure 3a,c). Such an effect is usually explained by the aggregation of the protein molecules that prevents their refolding [25]. These data points were not included in the estimation of the inactivation rate constants as well as the points of the induction phase.
Comparing thermal inactivation of two luciferases, we found that the V. harveyi enzyme is less sensitive to heating, since its remaining activity was always greater than that of the P. leiognathi enzyme under the same conditions. The presence of 30% sucrose appeared to reduce the rate of thermal inactivation of both luciferases (Table 1 and Figure 3c,d). The effect was the most pronounced for P. leiognathi luciferase at 50 °C. In this case, enzyme inactivation slowed down more than four times.
Using the rate constants, we performed the Arrhenius analysis of thermal inactivation of the two luciferases in the buffer and sucrose solution (Figure 4).
In spite of the protective effect of sucrose, evident from the decrease of the inactivation rate constants (Table 1), it was revealed that the activation energy of thermal inactivation, Ea, does not change in viscous medium. In particular, Ea in the buffer and sucrose solution was found to be 237 ± 30 and 224 ± 7 kJ/mol for V. harveyi luciferase and 255 ± 27 and 243 ± 47 kJ/mol for P. leiognathi luciferase, respectively. Such estimations do not allow for the distinguishing of the degree of thermal stability of the two luciferases as well.

2.3. Thermal Stability of Luciferases as Revealed by Differential Scanning Calorimetry

To compare thermal stability of the luciferases, their denaturation was studied by differential scanning calorimetry. Figure 5 shows the temperature dependencies of excess heat capacity functions (Cpexp) obtained for bacterial luciferases at pH 6.9 in the buffer and in sucrose solution. It can be seen that the curves demonstrate peaks of cooperative heat absorption corresponding to protein denaturation. In Table 2 the thermodynamic parameters of the thermal denaturation of the luciferases are summarized. All samples are characterized by the lack of calorimetric reversibility and an asymmetry of the heat absorption peaks. The degree of peak asymmetry (ΔQ/ΔQ+) is shown in Table 2. Generally, the asymmetry of the heat absorption peaks can be associated either with the non-equilibrium of the denaturation process or with the accumulation during denaturation of one or more intermediate states of different stability. Nevertheless, it can be concluded that V. harveyi luciferase has slightly higher temperature stability than the P. leiognathi one, both in buffer and in 30% sucrose. However, the enthalpy of thermal denaturation is greater in P. leiognathi luciferase than in V. harveyi. This indicates that the amount of the interactions (hydrophobic, van der Waals) that disrupt with increasing temperature is greater in P. leiognathi luciferase than in the V. harveyi one. The thermal denaturation of two luciferases also differs in the degree of cooperativeness: P. leiognathi luciferase has a higher cooperativeness (ΔT1/2 is less), i.e., the destruction of the native structure occurs in a narrower temperature range.
The addition of 30% sucrose to the buffer stabilized the luciferases: their Tm shifted by 4.7 K (P. leiognathi) and 5.0 K (V. harveyi). However, the sucrose effects on the process of thermal denaturation of two bacterial luciferases are different. In P. leiognathi enzyme, the addition of sucrose led to a broadening of the peak and an increase of its asymmetry. Such an effect of organic substances on the heat absorption curve is characteristic; it was observed on different proteins and explained by a slowdown of the denaturation kinetics process [26]. In the case of V. harveyi luciferase, the addition of 30% sucrose did not change either ΔT1/2 or ΔQ/ΔQ+. Thus, our DSC experiments indicate that the protein denaturation schemes could be different for two studied luciferases.

2.4. Temperature Dependent Molecular Dynamics of V. harveyi and P. leiognathi Luciferases

To evaluate the role of structural dynamics of the proteins in experimentally obtained patterns, we conducted a molecular dynamics (MD) computation of two bacterial luciferases surrounded by water molecules at low (5, 15 °C), medium (27 °C) and high (45, 60 °C) temperatures.
Firstly, based on the obtained MD trajectories, we analyzed three parameters: the root-mean-square deviation of the backbone atoms (RMSD), the radius of gyration (Rg), and the solvent accessible surface area (SASA). It allowed assessing the protein stability in the course of molecular dynamics and overall structure change. The results for the last 10 ns of the simulations are summarized in Table S1. No essential change of RMSD and Rg values for both enzymes was detected with increasing temperature. The RMSD values of V. harveyi luciferase were found to be lower than that of P. leiognathi luciferase (1.6–2.3 Å vs. 1.9–2.6 Å, respectively) and at 60 °C the standard deviation of the RMSD of V. harveyi luciferase appeared to be less than that of the P. leiognathi enzyme. It could reflect less deviation of the V. harveyi luciferase structure from the initial one at higher temperatures. The SASA values were obtained in the ranges of 284–290 and 293–301 × 102 Å2 for the V. harveyi and P. leiognathi luciferases, respectively. It indicates that the last enzyme could be less compact.
Additionally, the molecular dynamics of the two luciferases surrounded by the mixture of water and sucrose molecules (corresponding to concentration 30%) was analyzed at 27 °C and 60 °C. It is worth noting that at high temperature the presence of sucrose reduced the deviation of SASA from the average value for V. harveyi luciferase (Table S1), while for the P. leiognathi enzyme there was no effect.
Secondly, we calculated the average root-mean-square fluctuations (RMSF) of Cα atoms of the luciferases in order to examine the structural flexibility of protein segments at different temperatures. For each temperature the change of flexibility was estimated as a difference according to ΔRMSF = RMSFT − RMSF27, where RMSFT and RMSF27 are parameters for temperature T and 27 °C, respectively. In Figures S2 and S3 for two luciferases ΔRMSF of an α-subunit, which bears the enzyme active center, and a β-subunit, which stabilizes active enzyme, are shown [27]. The temperature effects on flexibility of a functional mobile loop in α-subunits were found for both enzymes (Figure 6a,b). The increased fluctuation of the region near αHis150 in V. harveyi luciferase was also observed (Figure S2). This amino acid residue is located at the apex of the hairpin loop structure (142–160 a. r.), which passes along the periphery of an α/β-interface of the luciferase [28]. Significant fluctuations in this region at 45 °C could expose a part of the α/β-interface to the solvent molecules, which could affect the enzyme activity.
The mobile loop is a particularly important functional region adjacent to the active site of the luciferase. It undergoes conformational changes upon enzyme-substrate complex formation [29]. We observed that at 5 °C two segments of the mobile loop change their flexibility in different manners as compared to 27 °C. The loop segment about 260–275 a. r. became less mobile in V. harveyi luciferase (Figure 6a), while its rigidity in the P. leiognathi enzyme did not change (Figure 6b). The loop segment at about 275–290 a. r., on the contrary, had a higher mobility in V. harveyi luciferase (Figure 6a) and a lower mobility in P. leiognathi enzyme under low temperatures (Figure 6b). Hence, we revealed that the mobile loop flexibility differs between the two luciferases at 5 °C, and there is moderate negative correlation between the ΔRMSF profiles of P. leiognathi and V. harveyi mobile loop segments, with a correlation coefficient of −0.52 (Table S2).
Inspecting the possible influence of sucrose on the mobile loop of the luciferase, we compared ΔRMSF of the proteins surrounded by water and sucrose-water mixture at 27 and 60 °C (Figure 6c,d, Figures S4 and S5). The obtained data indicate that the flexibility of the mobile loop was reduced in the presence of sucrose. The effect was more pronounced for P. leiognathi luciferase: fluctuations of its mobile loop at 60 °C in the presence of sucrose became close to the level observed in water at 27 °C.
The presence of sucrose molecules in the vicinity of amino acids forming the mobile loops of the luciferases during MD modeling was estimated using the minimum-distance distribution function, MDDF [30]. Based on MDDF we constructed the density maps of sucrose appearance at a distance of 1.5–3.5 Å from amino acids of the mobile loops (Figure 7a,b and Figure S6). It turned out that the sucrose atoms located most often within the first solvation shell (in 1.5–2.0 Å) of aspartate negatively charged side chains. Moreover, the density of sucrose molecules near the mobile loop residues of P. leiognathi luciferase increases at high temperatures (Figure 7b and Figure S6).

3. Discussion

We experimentally compared the temperature effects on two bacterial luciferases in the following directions: (i) the enzymes functioning under various temperatures, (ii) the rates of enzymes’ thermal inactivation, (iii) the thermal denaturation of luciferases. We aimed to reveal signs of protein adaptation to temperature variation for any of two enzymes and to connect them with structural properties of the proteins. However, as thermal adaptation of luminous bacteria could have occurred due to other mechanisms, such as, for example, the accumulation of extremolytes, we tested the influence of sucrose on the temperature dependencies, both functional and structural, of the bacterial luciferases.
The study of reaction kinetics at 5–45 °C revealed the different sensitivity of two luciferases to the temperature of solution (Figure 2 and Figure S1). In particular, P. leiognathi luciferase responds by the pronounced activity change for each temperature shift of 5 °C, while V. harveyi luciferase provides about the same peak intensity within a wide range of 20–35 °C. The observed thermolability of P. leiognathi luciferase and its higher activity under optimum conditions (at 20–25 °C) as compared with V. harveyi enzyme are consistent with the conception about the structural and functional features of cold-adapted proteins [31]. Indeed, the activity of P. leiognathi luciferase at 10–15 °C, although not maximal, is significantly higher than that of V. harveyi. Our molecular dynamics simulations shed light on the structural reasons for this by showing that there are two flexible motifs within the catalytic α-subunit that exhibit distinct mobility in the two luciferases at low temperatures (Figure 6a,b). Apparently, the observed change in structural movements does not provide the necessary mobility of the active site residues of V. harveyi enzyme for effective catalysis under cold conditions. Generally, the structure of V. harveyi luciferase was more stable at 60 °C as indicated by a smaller standard deviation of RMSD (Table S1). Such behavior is characteristic of proteins capable of functioning at higher temperatures, and not at lower ones. It is worth noting that the majority of known psychrophilic strains of luminous bacteria do belong to the Photobacterium and Aliivibrio genus, which have “fast” luciferases [15]. Thus, the data on the activity of the studied luciferases at different temperatures are consistent with their phylogeny and molecular dynamics.
Unexpectedly, the empiric kinetic parameter—the decay constant (kdecay), turned out to be more sensitive to temperature change for the reaction of V. harveyi luciferase than of P. leiognathi (Figure 2c,d). In the range of 10–30 °C, kdecay increases 4.4-fold for the first enzyme and 3.2-fold for the latter. The decay constant is known to depend on rates of three elementary steps of the reaction: the dark decay of peroxyflavin intermediate (kdd), aldehyde binding, and the formation of an electronically excited intermediate [32]. The contribution of kdd into kdecay substantially differs for two luciferases. In reaction with decanal at 20 °C, the kdecay values are 0.32 and 0.21 s−1 for the V. harveyi and P. leiognathi enzymes, whereas the corresponding kdd are 0.06 and 0.15 s−1 [33,34]. Thus, distinct temperature dependences of kdecay for two luciferases could be caused by two reasons: (i) by a stronger destabilizing effect on peroxyflavin intermediate in V. harveyi reaction and/or (ii) by temperature influence on the others than kdd rate constants which could be more pronounced in the V. harveyi reaction due to their high contribution into kdecay. Earlier, we demonstrated that for the reaction of P. leiognathi luciferase in the presence of sucrose the rate of the dark decay of peroxyflavin intermediate becomes lower, while the rate of the formation of an electronically excited intermediate increases as compared with the buffer [24]. This could underlie the contrasting changes of kdecay observed for two luciferases at 40–45 °C in a sucrose solution (Figure 2c,d).
The possible mechanism of sucrose’s influence on the stability of the reaction intermediates at high temperatures was revealed by molecular dynamics simulations. We observed that at 60 °C the sucrose significantly reduces the fluctuations of the functionally important mobile loop near the active site of P. leiognathi (Figure 6d), most likely through direct temporary contacts with aspartate side chains (Figure 7b). For V. harveyi luciferase stabilization of the mobile loop by sucrose is less pronounced (Figure 6c), which can cause the increase of the kdecay at higher temperatures. Moreover, it was found that sucrose promotes higher fluctuations of αHis150 within the hairpin loop of V. harveyi luciferase, which could also influence kdecay as well (Figure S4).
Based on the observed temperature effects on the activity of the enzymes in the presence of sucrose, we can conclude that this cosolvent provides no advantages for the functioning of bacterial luciferases either at high (35–45 °C) or at low (5–15 °C) temperatures. Furthermore, we observed that, at low temperatures, sucrose suppresses luciferase activity, probably due to high viscosity which leads to increased diffusion control of the reaction. However, along with activity, the structural stability of proteins also plays an important role in maintenance of their function under extreme conditions, and saccharides are known to be one of the best protein protectors against denaturation [35]. Hence, we further examined the kinetics of thermal inactivation of V. harveyi and P. leiognathi luciferases in the buffer and in the presence of sucrose. In addition, heat-induced denaturation of the enzymes was investigated using DSC. We observed that under the same conditions the thermal inactivation rate is always smaller for V. harveyi luciferase than for P. leiognathi (Table 1). Sucrose slows down the thermal inactivation of the enzymes by a factor of two to four times, but without a significant reduction in the activation energy of the process. This sucrose effect could also be a result of a decrease in mobile loop fluctuations at high temperatures, since the enhanced mobility of the protein segments (especially near the active center) could trigger the enzyme transition to an inactive state under heating. The thermal denaturation of the luciferases studied with DSC indicated that the temperature of the maximum of the heat absorption peak (Tm) is slightly higher for the V. harveyi protein than for the P. leiognathi one and it increases in sucrose solution. Thus, both approaches showed that V. harveyi luciferase has higher temperature stability than P. leiognathi and that sucrose is able to protect the proteins against thermal inactivation and denaturation.
In addition, we also found significant differences in ΔHcal between the two luciferases. Many interactions, such as hydration, hydrophobic effect, van der Waals interactions, and hydrogen bonding (between polar surfaces in protein and in water) contribute to the enthalpy of denaturation [36]. Since the spatial structures of P. leiognathi and V. harveyi luciferases are very similar, the difference in ΔHcal can be explained by the presence of a slightly larger number of unstructured parts in V. harveyi luciferase compared to P. leiognathi luciferase. This was confirmed by our early data on CD spectroscopy: V. harveyi luciferase has slightly less secondary structure compared to P. leiognathi luciferase [37]. According to our MD data, we also see that the hairpin loop within the V. harveyi structure was more flexible than in the P. leiognathi structure at high temperatures. It is worthy of note that our previous investigation of the unfolding of V. harveyi and P. leiognathi luciferases with urea showed that the latter has a more stable tertiary and secondary structure [37]. In agreement with these results, our current data on the enthalpy of thermal denaturation also indicate that P. leiognathi luciferase has a greater number of intramolecular interactions, which are disrupted with increasing temperature (parameter ΔHcal in Table 2).
The chemical unfolding of the bacterial luciferase can be considered much more studied than thermal denaturation [38]. A number of investigations have proven a four-state scheme for denaturation of V. harveyi luciferase with urea [39,40], where the unfolding of the C-terminal domain of the α-subunit is followed by the dissociation of the two subunits and the subsequent denaturation of them. In our previous work [37], we demonstrated that P. leiognathi luciferase in urea solutions follows the same unfolding pathway. However, the results presented in this work suggest that the mechanism of thermal denaturation of the two studied enzymes could differ. This is not very surprising, since chemical and thermal denaturation processes do not have to correlate, as the denaturant interacts specifically with the protein, while temperature generally affects intramolecular contacts indiscriminately [41].
Lack of knowledge about thermostability and the activity of bacterial luciferase at different temperatures is one of the issues when it is implemented as a reporter for the real-time monitoring of various processes in microbial or eukaryotic cells [42,43,44]. Previous studies demonstrated that in order to obtain a reproducible luminescent response and to understand this response, a detailed analysis of the characteristics of the proteins themselves is necessary, as well as the influence of the environment on their structure and function [43,45]. Therefore, our data has significant implications for the use of two types of luciferase with distinct sensitivity to temperature in terms of both activity and stability as reporters.

4. Materials and Methods

4.1. Materials

Lyophilized recombinant luciferase P. leiognathi (99% purity) was purchased from Biolumdiagnostika Ltd. (Krasnoyarsk, Russia). Luciferase V. harveyi was lyophilized after expression in the E. coli strain BL21 (DE3), and purification was as described in [37] in the Photobiology laboratory of the Institute of Biophysics SB RAS (Krasnoyarsk, Russia). Flavin mononucleotide (FMN, Sigma-Aldrich, Burlington, MA, USA) and decanal (Acros Organics, Fair Lawn, NJ, USA) were used as luciferase substrates. Ethylenediaminetetraacetic acid (EDTA, ROTH) was applied as an electron donor during FMN photoreduction. Sucrose (Gerbu, Heidelberg, Germany) in a concentration of 30 wt% was used as a model extermolyte. Stock solutions of the reactants were prepared in a potassium phosphate buffer (0.05 M, pH 6.9), except for decanal, which was dissolved in ethanol (2 × 10−3 M). The concentrations of FMN and luciferase were determined spectrophotometrically with the extinction coefficients of ε445 = 12,400 M−1 cm−1 and ε280 = 80,000 M−1 cm−1, respectively.

4.2. Activity and Thermal Inactivation of Bacterial Luciferases

The reaction kinetics for the luciferases was measured at 5–45 °C in a single-turnover assay using a SX20 stopped-flow spectrometer (Applied Photophysics, Leatherhead, UK) as described elsewhere [24]. Decanal and photoreduced flavin mononucleotide in the buffer or in the sucrose solution were loaded into one drive syringe of the spectrometer and pre-incubated there for 5 min. Similarly, bacterial luciferase in the buffer or sucrose solution was pre-incubated for 5 min in another drive syringe. The enzyme and substrates were then mixed in the measuring cell and the kinetics of light emission was registered for 10–30 s. The cell and drive syringes were thermostated at the required temperature. The concentrations in the reaction mixtures were 0.5 µM for luciferase, 15 µM for flavin and 50 µM for decanal.
The thermal inactivation rate of the luciferases was estimated by measuring the remaining activity of the enzymes after their incubation for various times under the required temperature in the range 40–55 °C. Solid-state thermostat Gnom (DNA technology, Moscow, Russia) was used for enzyme incubation. Reaction kinetics were measured at 20 °C in a single-turnover assay with a SX20 stopped-flow spectrometer as described above, but without incubation in drive syringes. The remaining activity (R) of luciferase was then calculated by percentage according to R = (Q*/Q0*)·100, where Q* and Q0* are areas under the kinetic curves obtained with thermally treated and thermally untreated enzymes, respectively. The rate constant of thermal inactivation, k, was determined by fitting the dependence of R on incubation time, t, with the function R = A∙ekt using Origin 8.0 software (OriginLab, Northampton, MA, USA). The coefficient of determination, R2, was >0.96 for all of the approximations. The activation energy of thermal inactivation (Ea) was calculated according to the Arrhenius equation from the slope of dependence lnk (1/T).

4.3. Differential Scanning Calorimetry

Calorimetric measurements were made using a precision scanning microcalorimeter SCAL-1 (Scal Co. Ltd., Pushchino, Russia) with 0.33-mL glass cells at a scanning rate of 1 K/min and under 2.5 atm pressure [46]. The protein concentrations ranged from 1.0 to 4.0 mg/mL in 50 mM potassium-sodium phosphate buffer with a pH of 6.9. In this concentration range, the thermodynamic parameters remained stable. The experimental calorimetric profiles were corrected for the calorimetric baseline, and the molar partial heat capacity functions were calculated in a standard manner. The excess heat capacity, Cpexp, was evaluated by subtracting the linearly extrapolated initial and final heat capacity functions with correction for the difference of these functions by using a sigmoid baseline [47]. A typical value for the partial specific volume for globular proteins (0.73 cm3/g) was accepted arbitrarily, since it does not influence the calculated excess heat capacity. In the calculation of molar thermodynamic quantities, the molecular weight used was 80,000 Da, which corresponds to the dimer state of the luciferases.

4.4. MD Simulation

The previously prepared structures of the V. harveyi and P. leiognathi luciferases [37] were taken for molecular dynamics simulations (MD) at different temperatures (5, 15, 27, 45, 60 °C) using the GROMACS 2020.4 software package and CHARMM36 force field [48,49]. The enzymes were solvated in a cubic box with periodic boundary conditions and at least 12 Å away from the protein to each of the box edges. Two strategies were applied for the solvation procedure: (1) TIP3P explicit water [50] was used for simulation at all studied temperatures; and (2) a mixture with sucrose molecules [51] adjusted to simulate the concentration at 30 wt% was exploited for molecular dynamics at 300 and 333 K. The net charge of the systems was neutralized with sodium ions. The minimization of each system was performed with the steepest descent method (maximum force of 1000.0 kJ/mol). The system was then heated, and sequential execution of isochoric and isobaric equilibrations was performed for 10 ns each with a 2 fs time step. A V-rescale thermostat at the required temperature and a Parrinello–Rahman barostat at 1 bar of pressure was used. Protein heavy atoms were restrained during equilibration [52,53]. The cutoff distance for the short-range non-bonded interactions was 10 Å. Electrostatic interactions were treated using the particle-mesh Ewald method [54]. An integration step of 2.0 fs was used, and bonds were constrained with the LINCS algorithm [55]. We conducted 100-ns molecular dynamics simulations for each system. Each MD run was repeated three times to obtain the average value of the structural parameters for the studied proteins.

5. Conclusions

In this study we aimed at comparing the effect of different temperatures on the structure and function of two bacterial luciferases: the Vibrio harveyi and Photobacterium leiognathi species. Despite high homology, these enzymes belong to different subfamilies with specific structural features of the active center and distinct reaction kinetics [20]. Thus, V. harveyi luciferase is classified as a “slow” enzyme, whereas P. leiognathi luciferase is a “fast” one. To the best of our knowledge, systematic studies of the thermolability of “slow” and “fast” bacterial luciferases has yet to be carried out.
The results of our study indicate that the “fast” P. leiognathi luciferase is thermolabile and has maximal activity in a narrower temperature range than does the “slow” V. harveyi luciferase. The latter is weakly sensitive to temperature change and remains active up to 35 °C. However, P. leiognathi luciferase provides twice the total quantum yield per single turnover under optimal conditions. At lower temperatures, the activity of the “fast” enzyme is significantly higher compared to the “slow” one as well. This is probably due to the peculiarities of its functional mobile loop dynamics, resulting in the ability to more effectively stabilize the reaction intermediates reflected in the smaller decay constant of the reaction. Moreover, our data revealed an interesting pattern that despite the greater number of intramolecular contacts disrupted during the unfolding of P. leognathi luciferase, it has slightly lower temperature stability. Thus, luciferase with the fast decay kinetics demonstrates the properties of the cold-adapted enzymes, while luciferase with the slow kinetics is more resistant to high temperatures. All of this reflects the well-known interplay between the activity and structural stability of enzymes, implying that more stable proteins could not provide a high rate of catalysis due to the rigidity of the structure.
In addition, we studied the effect of sucrose on the stability and activity of two bacterial luciferases. Our investigation has revealed that sucrose does not stabilize both enzymes in terms of activity, but it can protect the proteins from thermal denaturation. It is noteworthy that the influence of sucrose on the denaturation mechanisms of two luciferases seems to be different and requires more detailed examination.
Thus, after experimental research and molecular dynamics modeling of temperature effects. we have found that, despite high homology, the V. harveyi and P. leiognathi luciferases exhibit distinct features in response to heating. However, to answer the question of whether the found structural and functional peculiarities are a sign of the entire bacterial luciferase subfamily (“slow” or “fast”) or is an individual feature of only the studied enzymes, further research is required.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms23158119/s1.

Author Contributions

Conceptualization, E.V.N. and V.A.K.; methodology, E.V.N. and A.A.D.; software, L.A.S.; validation, A.A.D. and E.V.N.; formal analysis, A.E.L., T.N.M. and L.A.S.; investigation, A.E.L., L.A.S. and T.N.M.; resources, A.A.D. and E.V.N.; data curation, A.A.D.; writing—original draft preparation, A.A.D., A.E.L., T.N.M. and L.A.S.; writing—review and editing, E.V.N., A.A.D. and V.A.K.; visualization, A.E.L. and L.A.S.; supervision, E.V.N.; project administration, A.A.D.; funding acquisition, A.A.D. All authors have read and agreed to the published version of the manuscript.

Funding

The research was partially funded by the Ministry of Science and Higher Education of the Russian Federation (projects No. FSRZ-2020-0006); by RFBR, Krasnoyarsk Territory and Krasnoyarsk Regional Fund of Science (projects No. 20-44-243002); by RFBR (project No. 20-34-90118).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

The authors are grateful to Ludmila Burakova (Institute of Biophysics SB RAS) for purification and lyophilization of V. harveyi luciferase.

Conflicts of Interest

The authors declare that they have no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Gunde-Cimerman, N.; Plemenitaš, A.; Oren, A. Strategies of adaptation of microorganisms of the three domains of life to high salt concentrations. FEMS Microbiol. Rev. 2018, 42, 353–375. [Google Scholar] [CrossRef] [PubMed]
  2. Odermatt, P.D.; Miettinen, T.P.; Lemiere, J.; Kang, J.H.; Bostan, E.; Manalis, S.R.; Huang, K.C.; Chang, F. Variations of intracellular density during the cell cycle arise from tip-growth regulation in fission yeast. Elife 2021, 10, e64901. [Google Scholar] [CrossRef] [PubMed]
  3. Persson, L.B.; Ambati, V.S.; Brandman, O. Cellular control of viscosity counters changes in temperature and energy availability. Cell 2020, 183, 1572–1585. [Google Scholar] [CrossRef]
  4. Van Den Berg, J.; Boersma, A.J.; Poolman, B. Microorganisms maintain crowding homeostasis. Nat. Rev. Microbiol. 2017, 15, 309–318. [Google Scholar] [CrossRef]
  5. Kunte, H.J. Osmoregulation in bacteria: Compatible solute accumulation and osmosensing. Environ. Chem. 2006, 3, 94–99. [Google Scholar] [CrossRef]
  6. Becker, J.; Wittmann, C. Microbial production of extremolytes—High-value active ingredients for nutrition, health care, and well-being. Curr. Opin. Biotechnol. 2020, 65, 118–128. [Google Scholar] [CrossRef]
  7. Arakawa, T.; Timasheff, S. The stabilization of proteins by osmolytes. Biophys. J. 1985, 47, 411–414. [Google Scholar] [CrossRef]
  8. Bolen, D.W. Protein stabilization by naturally occurring osmolytes. Methods Mol. Biol. 2001, 168, 17–36. [Google Scholar]
  9. Åqvist, J.; Isaksen, G.V.; Brandsdal, B.O. Computation of enzyme cold adaptation. Nat. Rev. Chem. 2017, 1, 1–14. [Google Scholar] [CrossRef]
  10. Dunlap, P. Biochemistry and genetics of bacterial bioluminescence. In Bioluminescence: Fundamentals and Applications in Biotechnology—Volume 1. Advances in Biochemical Engineering/Biotechnology; Thouand, G., Marks, R., Eds.; Springer: Berlin/Heidelberg, Germany, 2014; Volume 144, pp. 37–64. [Google Scholar]
  11. Vannier, T.; Hingamp, P.; Turrel, F.; Tanet, L.; Lescot, M.; Timsit, Y. Diversity and evolution of bacterial bioluminescence genes in the global ocean. NAR Genom. Bioinform. 2020, 2, lqaa018. [Google Scholar] [CrossRef] [Green Version]
  12. Urbanczyk, H.; Ast, J.C.; Kaeding, A.J.; Oliver, J.D.; Dunlap, P.V. Phylogenetic analysis of the incidence of lux gene horizontal transfer in Vibrionaceae. J. Bacteriol. 2008, 190, 3494–3504. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Colepicolo, P.; Cho, K.; Poinar, G.; Hastings, J. Growth and luminescence of the bacterium Xenorhabdus luminescens from a human wound. Appl. Environ. Microbiol. 1989, 55, 2601–2606. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Bazhenov, S.V.; Khrulnova, S.A.; Konopleva, M.N.; Manukhov, I.V. Seasonal changes in luminescent intestinal microflora of the fish inhabiting the Bering and Okhotsk seas. FEMS Microbiol. Lett. 2019, 366, fnz040. [Google Scholar] [CrossRef]
  15. Burtseva, O.; Kublanovskaya, A.; Baulina, O.; Fedorenko, T.; Lobakova, E.; Chekanov, K. The strains of bioluminescent bacteria isolated from the White Sea finfishes: Genera Photobacterium, Aliivibrio, Vibrio, Shewanella, and first luminous Kosakonia. J. Photochem. Photobiol. B 2020, 208, 111895. [Google Scholar] [CrossRef] [PubMed]
  16. Brodl, E.; Winkler, A.; Macheroux, P. Molecular mechanisms of bacterial bioluminescence. Comput. Struct. Biotechnol. J. 2018, 16, 551–564. [Google Scholar] [CrossRef]
  17. Holzman, T.F.; Baldwin, T.O. The effects of phosphate on the structure and stability of the luciferases from Beneckea, harveyi, Photobacterium, fischeri, and, Photobacterium phosphoreum. Biochem. Biophys. Res. Commun. 1980, 94, 1199–1206. [Google Scholar] [CrossRef]
  18. Suadee, C.; Nijvipakul, S.; Svasti, J.; Entsch, B.; Ballou, D.P.; Chaiyen, P. Luciferase from Vibrio campbellii is more thermostable and binds reduced FMN better than its homologues. J. Biochem. 2007, 142, 539–552. [Google Scholar] [CrossRef]
  19. Szittner, R.; Meighen, E. Nucleotide sequence, expression, and properties of luciferase coded by lux genes from a terrestrial bacterium. J. Biol. Chem. 1990, 265, 16581–16587. [Google Scholar] [CrossRef]
  20. Deeva, A.A.; Temlyakova, E.A.; Sorokin, A.A.; Nemtseva, E.V.; Kratasyuk, V.A. Structural distinctions of fast and slow bacterial luciferases revealed by phylogenetic analysis. Bioinformatics 2016, 32, 3053–3057. [Google Scholar] [CrossRef]
  21. Nealson, K.H.; Hastings, J.W. Bacterial bioluminescence: Its control and ecological significance. Microbiol. Rev. 1979, 43, 496–518. [Google Scholar] [CrossRef]
  22. Makemson, J.C.; Fulayfil, N.R.; Landry, W.; Van Ert, L.M.; Wimpee, C.F.; Widder, E.A.; Case, J.F. Shewanella woodyi sp. nov., an exclusively respiratory luminous bacterium isolated from the Alboran Sea. Int. J. Syst. Evol. Microbiol. 1997, 47, 1034–1039. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Schmidt, T.; Kopecky, K.; Nealson, K. Bioluminescence of the insect pathogen Xenorhabdus luminescens. Appl. Environ. Microbiol. 1989, 55, 2607–2612. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Lisitsa, A.E.; Sukovatyi, L.A.; Bartsev, S.I.; Deeva, A.A.; Kratasyuk, V.A.; Nemtseva, E.V. Mechanisms of viscous media effects on elementary steps of bacterial bioluminescent reaction. Int. J. Mol. Sci. 2021, 22, 8827. [Google Scholar] [CrossRef] [PubMed]
  25. Brown, I.; Dafforn, T.; Fryer, P.; Cox, P. Kinetic study of the thermal denaturation of a hyperthermostable extracellular α-amylase from Pyrococcus furiosus. Biochim. Biophys. Acta Proteins Proteom. 2013, 1834, 2600–2605. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Griebenow, K.; Klibanov, A.M. On protein denaturation in aqueous-organic mixtures but not in pure organic solvents. J. Am. Chem. Soc. 1996, 118, 11695–11700. [Google Scholar] [CrossRef]
  27. Campbell, Z.T.; Weichsel, A.; Montfort, W.R.; Baldwin, T.O. Crystal structure of the bacterial luciferase/flavin complex provides insight into the function of the β subunit. Biochemistry 2009, 48, 6085–6094. [Google Scholar] [CrossRef]
  28. Fisher, A.J.; Thompson, T.B.; Thoden, J.B.; Baldwin, T.O.; Rayment, I. The 1.5-A resolution crystal structure of bacterial luciferase in low salt conditions. J. Biol. Chem. 1996, 271, 21956–21968. [Google Scholar] [CrossRef] [Green Version]
  29. Campbell, Z.T.; Baldwin, T.O.; Miyashita, O. Analysis of the bacterial luciferase mobile loop by replica-exchange molecular dynamics. Biophys. J. 2010, 99, 4012–4019. [Google Scholar] [CrossRef] [Green Version]
  30. Martínez, L.; Shimizu, S. Molecular interpretation of preferential interactions in protein solvation: A solvent-shell perspective by means of minimum-distance distribution functions. J. Chem. Theory Comput. 2017, 13, 6358–6372. [Google Scholar] [CrossRef] [Green Version]
  31. Siddiqui, K.S. Defying the activity–stability trade-off in enzymes: Taking advantage of entropy to enhance activity and thermostability. Crit. Rev. Biotechnol. 2017, 37, 309–322. [Google Scholar] [CrossRef]
  32. Cline, T.W.; Hastings, J. Mutationally altered bacterial luciferase. Implications for subunit functions. Biochemistry 1972, 11, 3359–3370. [Google Scholar] [CrossRef] [PubMed]
  33. Lisitsa, A.; Sukovatyi, L.; Kratasyuk, V.; Nemtseva, E. Viscous media slow down the decay of the key intermediate in bacterial bioluminescent reaction. Dokl. Biochem. Biophys. 2020, 492, 162–165. [Google Scholar] [CrossRef] [PubMed]
  34. Meighen, E.A.; Bartlet, I. Complementation of subunits from different bacterial luciferases. Evidence for the role of the beta subunit in the bioluminescent mechanism. J. Biol. Chem. 1980, 255, 11181–11187. [Google Scholar] [CrossRef]
  35. Yancey, P.H.; Clark, M.E.; Hand, S.C.; Bowlus, R.D.; Somero, G.N. Living with water stress: Evolution of osmolyte systems. Science 1982, 217, 1214–1222. [Google Scholar] [CrossRef] [PubMed]
  36. Makhatadze, G.I.; Privalov, P.L. Energetics of protein structure. Adv. Protein Chem. 1995, 47, 307–425. [Google Scholar]
  37. Nemtseva, E.V.; Gulnov, D.V.; Gerasimova, M.A.; Sukovatyi, L.A.; Burakova, L.P.; Karuzina, N.E.; Melnik, B.S.; Kratasyuk, V.A. Bacterial luciferases from Vibrio harveyi and Photobacterium leiognathi demonstrate different conformational stability as detected by time-resolved fluorescence spectroscopy. Int. J. Mol. Sci. 2021, 22, 10449. [Google Scholar] [CrossRef]
  38. Petushkov, V.; Kratasyuk, G.; Kratasyuk, V.; Belobrov, P. Thermal inactivation of bacterial luciferase. Biochemistry 1982, 47, 1504–1508. [Google Scholar]
  39. Inlow, J.K.; Baldwin, T.O. Mutational analysis of the subunit interface of Vibrio harveyi bacterial luciferase. Biochemistry 2002, 41, 3906–3915. [Google Scholar] [CrossRef]
  40. Noland, B.W.; Baldwin, T.O. Demonstration of two independently folding domains in the α-subunit of bacterial luciferase by preferential ligand binding-induced stabilization. Biochemistry 2003, 42, 3105–3112. [Google Scholar] [CrossRef]
  41. Narayan, A.; Bhattacharjee, K.; Naganathan, A.N. Thermally versus chemically denatured protein states. Biochemistry 2019, 58, 2519–2523. [Google Scholar] [CrossRef] [Green Version]
  42. Dorn, J.G.; Frye, R.J.; Maier, R.M. Effect of temperature, pH, and initial cell number on luxCDABE and nah gene expression during naphthalene and salicylate catabolism in the bioreporter organism Pseudomonas putida RB1353. Appl. Environ. Microbiol. 2003, 69, 2209–2216. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Cui, B.; Zhang, L.; Song, Y.; Wei, J.; Li, C.; Wang, T.; Wang, Y.; Zhao, T.; Shen, X. Engineering an enhanced, thermostable, monomeric bacterial luciferase gene as a reporter in plant protoplasts. PLoS ONE 2014, 9, e107885. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Close, D.; Xu, T.; Smartt, A.; Rogers, A.; Crossley, R.; Price, S.; Ripp, S.; Sayler, G. The evolution of the bacterial luciferase gene cassette (lux) as a real-time bioreporter. Sensors 2012, 12, 732–752. [Google Scholar] [CrossRef] [Green Version]
  45. Mackey, B.; Cross, D.; Park, S. Thermostability of bacterial luciferase expressed in different microbes. J. Appl. Bacteriol. 1994, 77, 149–154. [Google Scholar] [CrossRef] [PubMed]
  46. Senin, A.A.; Potekhin, S.; Tiktopulo, E.; Filomonov, V. Differential scanning microcalorimeter SCAL-1. J. Therm. Anal. Calorim. 2000, 62, 153–160. [Google Scholar] [CrossRef]
  47. Privalov, P.L.; Potekhin, S.A. Scanning microcalorimetry in studying temperature-induced changes in proteins. Methods Enzymol. 1986, 131, 4–51. [Google Scholar] [PubMed]
  48. Van Der Spoel, D.; Lindahl, E.; Hess, B.; Groenhof, G.; Mark, A.E.; Berendsen, H.J. GROMACS: Fast, flexible, and free. J. Comput. Chem. 2005, 26, 1701–1718. [Google Scholar] [CrossRef]
  49. Best, R.B.; Zhu, X.; Shim, J.; Lopes, P.E.; Mittal, J.; Feig, M.; MacKerell, A.D., Jr. Optimization of the additive CHARMM all-atom protein force field targeting improved sampling of the backbone ϕ, ψ and side-chain χ1 and χ2 dihedral angles. J. Chem. Theory Comput. 2012, 8, 3257–3273. [Google Scholar] [CrossRef] [Green Version]
  50. Jorgensen, W.L.; Chandrasekhar, J.; Madura, J.D.; Impey, R.W.; Klein, M.L. Comparison of simple potential functions for simulating liquid water. J. Chem. Phys. 1983, 79, 926–935. [Google Scholar] [CrossRef]
  51. Guvench, O.; Mallajosyula, S.S.; Raman, E.P.; Hatcher, E.; Vanommeslaeghe, K.; Foster, T.J.; MacKerell, A.D., Jr. CHARMM additive all-atom force field for carbohydrate derivatives and its utility in polysaccharide and carbohydrate—Protein modeling. J. Chem. Theory Comput. 2011, 7, 3162–3180. [Google Scholar] [CrossRef] [Green Version]
  52. Bussi, G.; Donadio, D.; Parrinello, M. Canonical sampling through velocity rescaling. J. Chem. Phys. 2007, 126, 014101. [Google Scholar] [CrossRef] [Green Version]
  53. Parrinello, M.; Rahman, A. Polymorphic transitions in single crystals: A new molecular dynamics method. J. Appl. Pys. 1981, 52, 7182–7190. [Google Scholar] [CrossRef]
  54. Darden, T.; York, D.; Pedersen, L. Particle mesh Ewald: An N log (N) method for Ewald sums in large systems. J. Chem. Phys. 1993, 98, 10089–10092. [Google Scholar] [CrossRef] [Green Version]
  55. Hess, B.; Bekker, H.; Berendsen, H.J.; Fraaije, J.G. LINCS: A linear constraint solver for molecular simulations. J. Comput. Chem. 1997, 18, 1463–1472. [Google Scholar] [CrossRef]
Figure 1. The kinetics of bioluminescence emission in reactions, catalyzed by luciferases from V. harveyi (a) and P. leiognathi (b) at low (5 °C, blue lines), medium (25 °C, green lines) and high (40 °C, red lines) temperature in the buffer (solid lines) and in sucrose solution (dashed lines). Concentrations in the reaction mixture were: luciferase—0.5 µM, flavin—15 µM, decanal—50 µM.
Figure 1. The kinetics of bioluminescence emission in reactions, catalyzed by luciferases from V. harveyi (a) and P. leiognathi (b) at low (5 °C, blue lines), medium (25 °C, green lines) and high (40 °C, red lines) temperature in the buffer (solid lines) and in sucrose solution (dashed lines). Concentrations in the reaction mixture were: luciferase—0.5 µM, flavin—15 µM, decanal—50 µM.
Ijms 23 08119 g001
Figure 2. Dependence of kinetic parameters of luciferases from V. harveyi (a,c) and P. leiognathi (b,d) on temperature in the buffer (empty squares) and sucrose solution (filled squares): the peak intensity, Imax, (a,b) and the decay constant, kdecay, (c,d). The dashed lines are a guide for the eye. Concentrations in the reaction mixture were: luciferase—0.5 µM, flavin—15 µM, decanal—50 µM.
Figure 2. Dependence of kinetic parameters of luciferases from V. harveyi (a,c) and P. leiognathi (b,d) on temperature in the buffer (empty squares) and sucrose solution (filled squares): the peak intensity, Imax, (a,b) and the decay constant, kdecay, (c,d). The dashed lines are a guide for the eye. Concentrations in the reaction mixture were: luciferase—0.5 µM, flavin—15 µM, decanal—50 µM.
Ijms 23 08119 g002
Figure 3. The dependence of the remaining activity of V. harveyi (a,c) and P. leiognathi (b,d) luciferases in the buffer (a,b, empty markers) and sucrose solution (c,d, filled markers) on the incubation time at selected temperatures (45 °C—circles, 50 °C—squares, 55°C—triangles) in semi-logarithmic scales. Lines refer to approximations with exponential function. Concentrations in the reaction mixture were: luciferase—0.5 µM, flavin—15 µM, decanal—50 µM. The activity was measured at 20 °C after the immediate cooling of the incubated enzyme.
Figure 3. The dependence of the remaining activity of V. harveyi (a,c) and P. leiognathi (b,d) luciferases in the buffer (a,b, empty markers) and sucrose solution (c,d, filled markers) on the incubation time at selected temperatures (45 °C—circles, 50 °C—squares, 55°C—triangles) in semi-logarithmic scales. Lines refer to approximations with exponential function. Concentrations in the reaction mixture were: luciferase—0.5 µM, flavin—15 µM, decanal—50 µM. The activity was measured at 20 °C after the immediate cooling of the incubated enzyme.
Ijms 23 08119 g003aIjms 23 08119 g003b
Figure 4. Arrhenius plot of the thermal inactivation of V. harveyi (a) and P. leiognathi (b) luciferases in the buffer (empty circles) and sucrose solution (filled circles). The k is the enzyme thermal inactivation rate constant at each respective temperature, taken from Table 1. Dashed lines refer to a linear approximation.
Figure 4. Arrhenius plot of the thermal inactivation of V. harveyi (a) and P. leiognathi (b) luciferases in the buffer (empty circles) and sucrose solution (filled circles). The k is the enzyme thermal inactivation rate constant at each respective temperature, taken from Table 1. Dashed lines refer to a linear approximation.
Ijms 23 08119 g004
Figure 5. Temperature dependence of excess heat capacity functions for bacterial luciferases from V. harveyi (green lines) and P. leiognathi (blue lines) in the buffer (solid lines) and in 30% sucrose (dashed lines) at pH 6.9.
Figure 5. Temperature dependence of excess heat capacity functions for bacterial luciferases from V. harveyi (green lines) and P. leiognathi (blue lines) in the buffer (solid lines) and in 30% sucrose (dashed lines) at pH 6.9.
Ijms 23 08119 g005
Figure 6. ΔRMSF of Cα atoms of V. harveyi luciferase mobile loop (262–291 a. r.) in water (a) and water-sucrose (c) and P. leiognathi luciferase mobile loop (261–290 a. r.) in water (b) and water-sucrose (d) at different temperatures. The positive value of ΔRMSF corresponds to a more flexible segment as compared with the structure in water at 27 °C, while the negative ΔRMSF corresponds to a more rigid segment. Correlation coefficients between the ΔRMSF profiles of V. harveyi and P. leiognathi mobile loops are presented in Table S2.
Figure 6. ΔRMSF of Cα atoms of V. harveyi luciferase mobile loop (262–291 a. r.) in water (a) and water-sucrose (c) and P. leiognathi luciferase mobile loop (261–290 a. r.) in water (b) and water-sucrose (d) at different temperatures. The positive value of ΔRMSF corresponds to a more flexible segment as compared with the structure in water at 27 °C, while the negative ΔRMSF corresponds to a more rigid segment. Correlation coefficients between the ΔRMSF profiles of V. harveyi and P. leiognathi mobile loops are presented in Table S2.
Ijms 23 08119 g006aIjms 23 08119 g006b
Figure 7. The density map of sucrose molecules near the mobile loop residues of V. harveyi (a) and P. leiognathi (b) luciferase at 60 °C during simulation time. The intensity of the green color indicates the probability of sucrose appearance in the distance r from the residue.
Figure 7. The density map of sucrose molecules near the mobile loop residues of V. harveyi (a) and P. leiognathi (b) luciferase at 60 °C during simulation time. The intensity of the green color indicates the probability of sucrose appearance in the distance r from the residue.
Ijms 23 08119 g007aIjms 23 08119 g007b
Table 1. The thermal inactivation rate constants for the P. leiognathi and V. harveyi luciferases in the buffer and sucrose solutions (30%) at different temperatures.
Table 1. The thermal inactivation rate constants for the P. leiognathi and V. harveyi luciferases in the buffer and sucrose solutions (30%) at different temperatures.
V. harveyiP. leiognathi
T, °Ck, min−1T, °Ck, min−1
BufferSucroseBufferSucrose
450.15 ± 0.010.07 ± 0.01400.02 ± 0.0010.01 ± 0.001
480.71 ± 0.110.17 ± 0.01
501.36 ± 0.100.39 ± 0.02450.12 ± 0.010.05 ± 0.01
521.63 ± 0.240.80 ± 0.06500.52 ± 0.010.15 ± 0.01
552.69 ± 0.321.14 ± 0.05551.10 ± 0.030.56 ± 0.05
Table 2. The thermodynamic parameters of the thermal denaturation of the bacterial luciferases *.
Table 2. The thermodynamic parameters of the thermal denaturation of the bacterial luciferases *.
LuciferaseSolventTm, ±0.1 K (°C)ΔT1/2, ±0.1 KΔHcal, kJ/(K⋅mol)ΔQ/ΔQ+
P. leiognathiBuffer318.5 (45.4)4.2997 ± 801.3
30% sucrose323.2 (50.1)6.3850 ± 501.6
V. harveyiBuffer320.1 (47.0)7.4552 ± 451.3
30% sucrose325.1 (52.0)7.5252 ± 301.3
* Tm—the temperature of the maximum of heat absorption peak; ΔT1/2—half-width of the heat absorption peak; ΔHcal—calorimetric enthalpy; ΔQ/ΔQ+—degree of asymmetry, where ΔQ—heat released when warming up to temperature Tm, ΔQ+—heat released after Tm.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Deeva, A.A.; Lisitsa, A.E.; Sukovatyi, L.A.; Melnik, T.N.; Kratasyuk, V.A.; Nemtseva, E.V. Structure-Function Relationships in Temperature Effects on Bacterial Luciferases: Nothing Is Perfect. Int. J. Mol. Sci. 2022, 23, 8119. https://doi.org/10.3390/ijms23158119

AMA Style

Deeva AA, Lisitsa AE, Sukovatyi LA, Melnik TN, Kratasyuk VA, Nemtseva EV. Structure-Function Relationships in Temperature Effects on Bacterial Luciferases: Nothing Is Perfect. International Journal of Molecular Sciences. 2022; 23(15):8119. https://doi.org/10.3390/ijms23158119

Chicago/Turabian Style

Deeva, Anna A., Albert E. Lisitsa, Lev A. Sukovatyi, Tatiana N. Melnik, Valentina A. Kratasyuk, and Elena V. Nemtseva. 2022. "Structure-Function Relationships in Temperature Effects on Bacterial Luciferases: Nothing Is Perfect" International Journal of Molecular Sciences 23, no. 15: 8119. https://doi.org/10.3390/ijms23158119

APA Style

Deeva, A. A., Lisitsa, A. E., Sukovatyi, L. A., Melnik, T. N., Kratasyuk, V. A., & Nemtseva, E. V. (2022). Structure-Function Relationships in Temperature Effects on Bacterial Luciferases: Nothing Is Perfect. International Journal of Molecular Sciences, 23(15), 8119. https://doi.org/10.3390/ijms23158119

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop