Skip to main content

Exploring the impact of taurine on the biochemical properties of urate oxidase: response surface methodology and molecular dynamics simulation


This paper investigates the impact of taurine as an additive on the structural and functional stability of urate oxidase. First, the effect of the processing parameters for the stabilization of Urate Oxidase (UOX) using taurine was examined using the response surface methodology (RSM) and the central composite design (CCD) model. Also, the study examines thermodynamic and kinetic parameters as well as structural changes of urate oxidase with and without taurine. Fluorescence intensity changes indicated static quenching during taurine binding. The obtained result indicates that taurine has the ability to preserve the native structural conformation of UOX. Furthermore, molecular dynamics simulation is conducted in order to get insights into the alterations in the structure of urate oxidase in the absence and presence of taurine under optimal conditions. The molecular dynamics simulation section investigated the formation of hydrogen bonds (H-bonds) between different components as well as analysis of root mean square deviation (RMSD), root mean square fluctuations (RMSF) and secondary structure. Lower Cα-RMSD and RMSF values indicate greater stabilization of the taurine-treated UOX structure compared to the free enzyme. The results of molecular docking indicate that the binding of taurine to the UOX enzyme through hydrophobic interactions is associated with a negative value for the Gibbs free energy.


Urate oxidase (UOX) (EC, is an enzyme that binds to copper and catalyzes the chemical reaction of uric acid oxidation, resulting in the production of hydrogen peroxide (H2O2) and 5-hydroxyisourate. Uricase has the potential to serve as a therapeutic enzyme for the purpose of reducing the concentration of uric acid in the plasma [1]. In marine animals, it is metabolized further to urea and glyoxylic acid and is thus not excreted [2,3,4]. Urate oxidase is frequently observed in bacterial and fungal organisms. Nevertheless, it is important to note that humans lack uricase enzymes as a result of genetic abnormalities in the urate oxidase gene [5, 6]. The utilization of recombinant uricase enzyme has proven to be effective in the management of hyperuricemia, exhibiting a minimal occurrence of hypersensitivity reactions and tophaceous gout [6,7,8]. Furthermore, the uricase enzyme serves as a reagent in clinical diagnostic kits utilized for the enzymatic quantification of uric acid [6, 7]. One of the significant challenges associated with the utilization of therapeutic enzymes is to their limited stability. One approach to improve protein stability is the utilization of additives [9]. Living organisms, including both prokaryotic and eukaryotic cells, when exposed to adverse environmental conditions, such as high salt concentrations, extreme temperatures, or low temperatures, employ a similar mechanism to safeguard their proteins. This mechanism involves the synthesis of organic compounds known as osmolytes. Stabilizing osmolytes are organic compounds with a low molecular mass that raise the midpoint of thermal denaturation. They are divided into amino acids and their derivatives, polyhydric alcohols, sugars, and methylammonium derivatives [10,11,12]. Osmolytes can be classified into two distinct classes, namely compatible and counteracting, based on their functional mechanisms in relation to protein activity [13]. The mechanistic effect of each osmolyte on the performance and stability of the protein is unique and special. Note that sometimes the proteins that have been exposed to osmolytes and folded are more active than in their natural state [13, 14]. Taurine osmolyte is used to stabilize the structure of proteins. Taurine osmolyte is an organic acid with a molar mass of 125.14 g/mol found in many natural food sources. Taurine osmolyte is derived from cysteine and biosynthesized in the body; it is also produced via chemical synthesis for commercial purposes. Proteins that have high thermal stability have many applications in the field of biotechnology. The subject discussed, which is of great importance, is increasing the resistance of proteins in acute conditions. The comprehensive assessment of the stability and functionality of the uricase enzyme in the presence of osmolytes has yet to be fully elucidated. In this study, the correlation between protein stability and activity alterations when exposed to taurine osmolyte was examined. This study also looked at how taurine osmolyte affects the activity and structural stability of the recombinant urate oxidase enzyme. The investigation focused on assessing the impact of taurine on the enzymatic activity of recombinant uricase through the utilization of the response surface methodology (RSM) in order to determine the best conditions. The response surface method was used to achieve recommended optimal conditions, including the concentration of taurine osmolyte, incubation time, and temperature. Thermodynamic and kinetic parameters, and structural studies including intrinsic, and ANS fluorescence of pristine urate oxidase and taurine-treated uricase enzyme were examined. The obtained results indicated that the thermal stability of the untreated recombinant uricase enzyme was comparatively poorer in comparison to the uricase enzyme treated with taurine osmolyte. The study employed atomistic molecular dynamics (MD) simulations to get structural insights into the experimental findings and to offer further information regarding the impact of taurine. Also, MD simulation was applied to check the stability and flexibility of the formed complexes [15]. The findings indicate that the presence of taurine led to an increase in the proportion of α-helical and β-sheet structures in the recombination uricase enzyme.

Materials and methods


Boric acid (CAS No. 10043-35-3), Taurine (CAS No. 107-35-7; purity ≥ 99%), and uric acid (CAS No. 69-93-2, ≥ 99%) were from Sigma-Aldrich. The compounds KH2PO4 (CAS No. 7778-77-0) and K2HPO4 (CAS No. 11-4-7758) were procured from the Merck company. Affinity chromatography column of Ni-NTA resin for His6-tagged proteins was purchased from Qiagen. NaCl, yeast extract, Tris-HCl (CAS No. 1185-53-1), tryptone (CAS No. 91079-40-2). The agar medium was used for both liquid and solid culturing purposes. Bradford (Bio-Rad Protein) solution for the quantification of protein concentration, agarose, Isopropyl β-D-1-thiogalactopyranoside (IPTG) (CAS No. 367-93-1), sodium acetate, kanamycin antibiotic, and imidazole were obtained from Sina-clon Co, Iran. The plasmid extraction kit was purchased from Takara company, Japan. Acrylamide (CAS No. 79-06-1), sodium chloride (CAS No. 7647-14-5), glycine, and all SDS-PAGE chemicals were obtained from Merck company.

Microorganisms and media composition

The urate oxidase gene was expressed using the pET method developed by Novagen. The host organism for plasmid clones in this study was Escherichia coli, commonly referred to as E. coli (Invitrogen, Carlsbad, CA, USA). The strain of E. coli used in this study for protein expression is BL21 (DE3). The plasmid pET-28a was employed in this study to express the recombinant uricase enzyme. The expression was regulated by the T7-promoter, which was induced by isopropyl-β-D-1-thiogalactopyranoside (IPTG) obtained from Sigma Aldrich, USA. The bacterial strains were cultured using Luria-Bertani (LB) medium, which consisted of 1% (w/v) tryptone, 0.5% (w/v) yeast extract, and 1% (w/v) NaCl. For LB agar, 2% (w/v) agar was also used in the medium. The pH of the medium was adjusted to 7, and kanamycin was added at a concentration of 100 µg/ml as required.

Purification of the urate oxidase enzyme expressed in E. Coli

For recombinant uricase enzyme synthesis, the plasmid vector pET-28a was used as an expression vector, and the Escherichia coli BL21 strain was employed as the expression host. The cells obtained from a 100 ml culture were subjected to lysis and afterwards resuspended in 1–3 ml of Tris-HCL buffer solution (pH 8) with a concentration of 50 mmol. The process of sonication was performed on ice for a duration of 3 min, with a cycle of 70 and intervals of 0.5 s. The lysate was subjected to centrifugation at a speed of 13,000 g, at a temperature of 4 ˚C, for a duration of 20 min. The resulting supernatants after centrifuging were collected for the purpose of purification. The liquid portion of the sample (supernatant) was loaded onto chromatography columns containing Ni-Sepharose, and subsequent purification steps were carried out. Initially, the columns were washed using a wash solution with a concentration of 50 mM (pH 8). This wash buffer consisted of 2.5 M NaCl, 200 mM NaH2PO4, and 20 mM imidazole. The purpose of this washing step was to remove any nonspecific protein bonds that may have occurred within the columns. Following this, purification was performed utilizing an elution buffer comprising NaCl and NaH2PO4, along with different concentrations of imidazole ranging from 50 to 300 mM [16]. The protein concentration was determined using the Bradford test (Bio-Rad Protein, 1976), with bovine serum albumin (BSA) serving as the standard protein. In the present investigation, the SDS-PAGE technique, utilizing Laemmli’s buffer system, was utilized to determine the purity of the purified enzyme.

Uricase assay

The blank reaction mixture is composed of a Tris-HCl buffer with a concentration of 0.1 M and a pH of 9.5. Additionally, the mixture contains 1 mM of EDTA. The reaction mixture comprised 0.9 ml of air-saturated 0.1 M Tris–HCl buffer (pH 9.5) containing 1 millimolar EDTA, 0.1 millimolar uric acid, and 20 µl of urate oxidase enzyme solution. The activity of the uricase enzyme was assessed by measuring the decrease in absorbance at a wavelength of 293 nm. The one unit of enzyme activity (U) was defined by determining the amount of enzyme necessary to catalyze the conversion of 1 µmol of uric acid to allantoin within a duration of 1 min (definition A) under room temperature and pH 9.5 conditions [17]. The following equation can express the activity of recombinant uricase enzyme:

$${enzyme}_{\left(\raisebox{1ex}{unite}\!\left/ \!\raisebox{-1ex}{ml}\right.\right)}=\frac{({\varDelta A}_{293 nm}/minTest-{\varDelta A}_{293 nm}/minBlank)\left(\text{t}\text{o}\text{t}\text{a}\text{l}\,\text{v}\text{o}\text{l}\text{u}\text{m}\text{e}\right)\left(\text{d}\text{f}\right)}{\left(12.6\right) 0.02}$$

The symbol “df” represents the dilution factor, the value 12.6 corresponds to the millimolar (mM) extinction coefficient of uric acid at a wavelength of 293 nm, and the quantity 0.02 (ml) represents the volume of the uricase enzyme.

RSM methodology

The influence of the variables affecting the activity of the urate oxidase enzyme in the presence of taurine osmolyte was investigated with a full fraction design. Three operating variables including temperatures, incubation time, and the concentration of taurine were examined. In order to optimize the three most prominent factors, 18 full factorial central composite designs were implemented using the response surface methodology. This was accomplished through the utilization of Design-Expert software version 11.

Optimum temperature, incubation time, and taurine concentration by RSM

In order to optimize the stability process of the uricase enzyme, central composite design (CCD) and the response surface methodology (RSM) were employed. This approach aimed to attain the utmost efficacy and maximum level of stability for the urate oxidase enzyme. Based on the first experimental findings, three factors were identified as the independent variables: incubation time (x1), taurine concentration (x2), and temperatures (x3). The response addressed in this study was the activity of the urate oxidase enzyme after stabilization, measured in units per milliliter (U/ml). The investigation focused on determining the optimal estimated temperature for the activity of the recombinant uricase enzyme. This was achieved by assessing the enzyme’s activity subsequent to pre-incubation at different temperatures, spanning from 16 to 40 ˚C, with intervals of 5 ˚C, for a duration of 6 min. In order to calculate the duration of the optimum incubation, the urate oxidase enzyme and taurine were subjected to incubation for a range of 5 to 35 min. The present study aimed to investigate and compare the catalytic activity of uricase enzymes in the presence and absence of taurine osmolyte under different temperature conditions and incubation times. To assess the impact of the optimal taurine concentration on the thermal stability of the recombinant urate oxidase enzyme, several concentrations of taurine osmolyte (0, 50, 150, 250, 350, and 450 mM) were examined with regard to the activity of the urate oxidase enzyme.

Thermal inactivation

In this study, two sets of experiments were conducted to investigate the thermostability of the urate oxidase enzyme treated with taurine. In order to investigate the process of thermal inactivation, the recombinant uricase enzyme was subjected to various temperatures, specifically 40, 55, and 60 ˚C, for a duration of 1 h. During this incubation period, measurements were taken at 5-minute intervals. It is important to note that the uricase enzyme was incubated without the presence of its substrate, uric acid. Subsequently, the samples were subjected to a cooling process after being placed on ice for a duration ranging from 30 to 60 min, with the aim of reinstating their native folding structure. Then, the residual activities of the drawn enzyme samples were measured. The residual activities of the drawn enzyme samples vs time plots were employed to determine the melting temperature of the enzyme. The melting temperature is defined as the temperature at which half of the native enzyme molecules have been denatured to an inactive conformation. The experiment was continued until no further decrease in enzyme activity (A/A0) was observed, indicating that the enzyme had reached its minimal level of activity, Amin. The Amin value refers to the minimal level of activity exhibited by an enzyme at a specific temperature, indicating the enzyme’s equilibrium condition and representing the overall quantity of active enzyme after cooling [18,19,20]. The experiment was conducted in triplicate, with each trial consisting of three measurements. The result presented is the average of these triple measurements.

Kinetic parameters

A more comprehensive comprehension of the stability of the recombinant urate oxidase enzyme through the utilization of taurine osmolyte can be achieved by evaluating the kinetic parameters, namely the maximum reaction velocity (Vmax) and the Michaelis Menten constant (Km). In order to determine the Michaelis-Menten constant (Km) and maximum reaction rate (Vmax) of uricase, experimental solutions were prepared by combining 20 mmol of boric acid with varying amounts of uric acid as the substrate. The determination of the Michaelis-Menten constant and maximum reaction velocity values was conducted using the Lineweaver-Burk plot, as described by Eqs. 8 and 9, respectively [21]. The turnover number value (Kcat) was obtained by dividing the maximum reaction velocity value by the final molar concentration of the recombinant urate oxidase enzyme. The experimental measurements for the untreated uricase enzyme and the uricase enzyme treated with taurine osmolyte were each replicated three times.


In the Lineweaver-Burk equation, \(\frac{{k}_{m}}{{V}_{max}}\) is the slope and \(\frac{1}{{V}_{max}}\) is the y-intercept

Fluorescence measurements

Intrinsic fluorescence

The experiment commenced by employing Amicon Ultra Centrifugal Filter Units to eliminate imidazole present in solutions containing purified uricase enzyme. Subsequently, desalting protein solutions were utilized to carry out the subsequent stages of the experiment. A Cary Eclipse fluorescence spectrophotometer (Varian) was used to assess the steady-state fluorescence of recombinant urate oxidase in the presence and absence of taurine osmolyte. The fluorescence spectra were obtained by exposing a sample to an excitation wavelength of 290 nm and measuring the emission wavelengths between 300 and 400 nm. This was done under varying concentrations of taurine at temperatures of 40, 55, and 60 °C. The temperature within the quartz cells was regulated using a thermostatically controlled circulating-water pump. The experiment involved the addition of the recombinant uricase enzyme at a concentration of 0.5 mg/mL to a Tris-HCl buffer solution with a volume of 170 mL. This experimental procedure was conducted over a concentration range spanning from 0 to 750 mM of taurine [22].

ANS fluorescence

The compound known as 8-anilino-1-naphthalenesulfonic acid (ANS) is a commonly employed fluorescent probe in the field of protein research, specifically for the purpose of characterizing protein binding sites. Initially, stock solutions of 8-anilino-1-naphthalenesulfonic acid (ANS) were produced using phosphate-buffered saline (PBS). Subsequently, the interactions between ANS and both treated and untreated enzymes were observed using a Cary Eclipse fluorescence spectrophotometer. All measurements were made at a temperature of 28 °C in the cuvette. Following the addition of ANS to the enzyme solution at a ratio of 1:30, the samples were thoroughly mixed and subsequently incubated in a light-restricted environment for a duration of 10 min before fluorescence measurements were taken. The excitation wavelength utilized for ANS fluorescence spectra was held constant at 370 nm, while emission wavelength scans were conducted within the range of 400 to 600 nm [23, 24].

Computational methods

Molecular dynamics simulation

The MD simulations were done on the A. flavus UOX crystallographic structure resolved at 2.30 Å (PDB ID: 1R56 from Protein Data Bank) [25]. Recent structures indicate that the MD simulation can help describe the stabilization or destabilization of the protein structure in the presence of osmolyte [9, 26,27,28,29]. The MD calculations were done on the native and taurine-treated enzyme using the GROMACS software under the optimum conditions (4.50 M taurine and 28 °C) obtained from the RSM method. The velocity rescale thermostat and Berendsen algorithm were used to maintain the temperature and pressure of the system (300 K and 1 bar, respectively) [30, 31]. The particle-mesh Ewald (PME) algorithm was used for long-range electrostatics [32]. Van der Waal interactions were modeled using Lennard-Jones 6–12 potentials with 14 Å cut-off. Bond lengths were constrained using the Linear Constraint Solver (LINCS) algorithm [33]. The system was energy minimized using the steepest descent; after an equilibration of 500 ps at constant pressure and temperature, the system was subjected to 50 ns of molecular dynamics (MD). Trajectories were saved at 2 ps intervals, and trajectory analyses were performed using the programs of GROMACS suite. The leap-frog algorithm was used as an integrator with a 2 fs time step to integrate Newton’s equations of motion. The gmx hbond utility in the GROMACS was applied to measure the number of hydrogen bonds (H-bonds). Analysis of the H-bonds employed a donor-acceptor cutoff distance of 3.5 Å and acceptor-donor-hydrogen cutoff angle of 30°. The root mean square deviation (RMSD) was calculated by comparing the simulated protein structure during simulation with the reference structure [34]. This study applied Cα-RMSD values to estimate protein structure similarity. To assess the residual mobility of the protein, the root-mean-square fluctuation (RMSF) was calculated for the Cα atoms of each residue. Therefore, the RMSF of residues indicates the flexibility or rigidity regions within the protein structure where those residues are located [35], a higher RMSF value suggests lower stability [36]. The solvent-accessible surface area (SASA) was calculated, which indicates the solvent-accessible surface of protein residues. An increase or decrease in the SASA values demonstrates the changes in the exposed amino acid residues and could influence the protein tertiary structure. The SASA values were calculated for the 50 ns trajectory of molecular dynamics simulation for free and treated-UOX enzymes. Also, the secondary structure analysis was done with the DSSP program [37]. In this protocol, hydrogen bonding and the other geometrical features were used to determine the secondary structure of the enzyme. To ensure the reproducibility of our results, we conducted the MD simulation three times with varying initial structures. The results of three runs showed a high degree of similarity (the results can be seen in the electronic Supplementary Material). This work presents the average results obtained from three simulations.

Molecular docking

We used MGLTools 1.5.6 AutoDockTools (ADTs) along with AutoGrid 4.2 and AutoDock 4.2 to prepare and perform molecular docking calculations between UOX enzyme and taurine molecule. The Lamarckian Genetic Algorithm (LGA) [38] was utilized to perform docking calculations. The target protein file was prepared using AutoDock 4.2 tool by excluding heteroatoms, structural water molecules, and co-factors, while retaining the protein-associated residues. This involved several steps such as assigning AD4 type atoms, calculating Gasteiger charges for each atom in the macromolecule, and adding polar hydrogens. Adding polar hydrogens is crucial to ensure accurate calculation of partial charge, while keeping all other values as default. The LIGPLOT program was utilized to illustrate the strengths of intermolecular interactions, including hydrogen bonds and hydrophobic interactions [39]. The energy evaluations were set to 2,500,000 and 27,000 generations, with a population size of 150. The gene mutation rate was set at 0.02, while the gene crossover rate was set at 0.8.

Results and discussion

Experimental results

Protein purification

The recombinant uricase enzyme was purified using previous methods, and the SDS-polyacrylamide gel electrophoresis technique was used to analyze the proteins eluted from the purification [40]. The observation of a single band with a molecular weight in the region of 35 kDa served as conclusive evidence for the expression of recombinant urate oxidase enzyme (Fig. 1) [41].

Fig. 1
figure 1

Recombinant urate oxidase purification. proteins eluted collected after Ni-NTA chromatography were visualized on SDS-PAGE.

Optimum temperature and incubation time plus taurine concentration by RSM

In order optimize the three most prominent factors (temperatures, incubation time, and the concentration of taurine), a total of 18 full factorial central composite designs (CCDs) were used in this study, using Design-Expert software version 11. The predicted and experimental results are presented in Table 1.

Table 1 Design matrix of experiments proposed by response surface methodology (RSM) and the experimental units acquired for each experiment

The results shown in Table 1 demonstrate a high degree of concordance between the predicted response values and the corresponding experimental data. The RSM model suggested a quadratic polynomial equation. The ultimate equation, expressed in terms of actual The RSM model suggested a quadratic polynomial equation. The final equation, in terms of the actual factor, can be stated as follows:

$$\mathbf{Activity}=+5.20717-0.049496\,\text{Incubation time}-0.006595\,\text{Taurine concentration}-0.151224\,\text{Temperature}-0.000152\,\text{Incubation time}\,\ast\,\text{Taurine concentration}+0.002154\,\text{Incubation time}\,\ast\,\text{Temperature}+0.000293\,\text{Taurine concentration}\,\ast\,\text{Temperature}+0.000557\,\text{Incubation time}^2+\text{7.12333E-06}\,\text{Taurine concentration}^2-\,\text{4.62963E-07}\,\text{Temperature}^2$$

The presence of a positive coefficient indicates a synergistic effect on the response, while a negative coefficient signifies an antagonistic effect. The RSM plot in Fig. 2 illustrates the impact of taurine concentration, incubation time, and temperature on the activity of the recombinant uricase enzyme.

Fig. 2
figure 2

The planned set of contour plots aims to illustrate the impact of three independent variables and their effects on response (activity): (a1/2) The influence of temperature, taurine concentration, and their interaction on activity, while keeping the incubation time constant at 25 min (b1/2) The impact of temperature, incubation time, and their interaction on activity, while maintaining a constant taurine concentration of 250 mM. (c1/2) The effects of taurine concentration, incubation time, and their interaction on activity, while holding the temperature constant at 28 °C

Based on the results shown in Fig. 2a1/2, it can be observed that at high temperatures with increasing taurine osmolyte concentration and at low temperatures with lower taurine osmolyte concentration, the enzyme achieves its highest possible level of activity. Moreover, it is evident from the data presented in Fig. 2b1/2 that the enzymatic activity increases at higher temperatures with increases in the incubation time of treated uricase enzyme with a constant concentration of taurine osmolyte. Moreover, Fig. 2c1/2 indicates that both a low taurine concentration and a long incubation time had a positive effect on the activity of the recombinant uricase enzyme, and a high taurine concentration and a short incubation time had the same effect. These results suggest that the effects of all three variables and their interactions with each other on the recombinant uricase enzyme activity have been significant. The optimum conditions for achieving maximum activity of the recombinant uricase enzyme were determined to be a temperature of 28 °C, a taurine concentration of 450 mM, and an incubation time of 25 min. The recombinant uricase enzyme activity of 2.05 U/ml was generally achieved under optimal conditions, as indicated in Table 1.

Thermal inactivation and stability

To investigate the mechanism of protein inactivation, it is crucial to assess kinetic and thermodynamic parameters, which are fundamental in thermal processes. The study investigated the irreversible thermal inactivation of enzymes that were treated with taurine and enzymes that were not treated with taurine, at various temperatures. In general, a population of protein molecules undergo irreversible structural changes, which can be seen from the reduction of the residual activity of the pristine enzyme at representative temperatures. Thus, a resetting of the operating temperature does not lead to the reestablishment of the initial enzyme activity. The time and temperature of incubation, as well as the specific protein type, determine the degree of reversibility for each protein. The examination of the thermo-inactivation process following a duration of 60 min indicated that the untreated enzyme exhibited lower stability compared to the enzyme subjected to taurine treatment across the temperature range of 40, 55, and 60 ˚C, particularly at elevated temperatures. The enzyme’s activity was shown to retain 40% and 20% of its initial activity after being incubated for 60 min at 55 °C with and without the presence of taurine, respectively (Fig. 3). After incubating the recombinant uricase enzyme with and without taurine osmolyte at representative temperatures within 40, 55, and 60 ˚C, irreversible thermal inactivation experiments were performed for 1 h for detailed analysis. The results indicated that during a 60 min incubation at 40 ˚C, the enzymes treated with taurine exhibited a retention of over 70% of their activity, as depicted in Fig. 4a, b, and c. The variation in melting points of proteins can be attributed to the distinctive primary sequence of amino acids that are characteristic of certain protein types. Various factors, such as pH and salt concentration, as well as post-translational modifications, can significantly impact the stability of protein structure and hence affect the melting temperature. Protein stability increases proportionally with greater values of the thermal melting temperature (Tm). In order to achieve the minimum level of activity, Amin, the data of treated and untreated uricase enzyme were plotted until no further drop in enzyme activity (ARes/A0) was observed (Fig. 5a/b). Then, according to Fig. 5c/d, the first parameter of thermal stability, i.e. the melting temperature of the treated and untreated uricase enzyme (Tm) was estimated.

Fig. 3
figure 3

Irreversible thermal inactivation of taurine treated enzyme (white) and pristine enzyme (black) at different temperatures for 60 min

Fig. 4
figure 4

Irreversible thermal inactivation of taurine treated enzyme 40 °C, b 55 °C, and c 60 °C for different time intervals

Fig. 5
figure 5

Determination of thermostability parameters: a and b, Reduction in the proportion of ARes compared with native activity (A0) over time in treated and untreated enzyme. Amin is the minimum activity at a given temperature that corresponds to the enzyme in the equilibrium state which represents the total amount of active enzyme upon cooling a the treated enzyme b the untreated enzyme. b and c, Estimation of melting temperature (Tm) that corresponds to the temperature when Amin drops down to 50% of the native activity (A0) in treated and untreated enzymes respectively.

Kinetic parameters

Based on previous kinetic studies, the Michaelis constant (Km) values of the urate oxidase enzyme were investigated under two conditions: in the absence and presence of taurine. Table 2 presents a comprehensive summary of the kinetic parameters acquired in the present study. The taurine-treated enzyme exhibited a twofold augmentation in Km compared to the untreated uricase enzyme. Based on the results, it can be shown that the enzyme treated with taurine osmolyte exhibited a greater affinity towards the substrate in comparison to the untreated uricase enzyme (as depicted in Fig. 6a, b, and Table 2). A high Km value indicates that a substantial amount of substrate is required to achieve enzyme saturation, implying that the enzyme exhibits a low affinity for the substrate. Conversely, a low Km value signifies a minimal quantity of substrate required to achieve enzyme saturation, indicating a higher affinity for the substrate. Hence, the data indicate that there is a change in the rate of conformational change in the presence and absence of taurine on the way to forming the Michaelis complex. An increase in the Km value of the enzyme treated with taurine may indicate that the enzyme’s structure underwent alterations during the stabilization process or that taurine introduces a spatial barrier that hinders the entry of the substrate (uric acid) into the enzyme’s active site [21, 24, 42, 43]. Furthermore, the enzyme’s catalytic efficiency or rate constants for the catalytic conversion of substrate into product (Kcat) and the ratio of Kcat to Km, which serves as an indicator of the specificity of the enzyme to the substrate, were examined in both the absence and presence of taurine osmolyte. The observed reduction in the value of this parameter for the enzyme in the presence of taurine osmolyte is probably attributable to the spatial barrier.

Fig. 6
figure 6

The kinetic parameters of uricase in the absence (a) and presence of taurine (b)

Table 2 The biochemical characteristics of uricase in the absence and presence of taurine

Fluorescence measurements

Fluorescence spectroscopy

The utilization of fluorescence spectroscopy has been widely employed to investigate the interactions between ligands and proteins, yielding valuable insights into the quenching mechanism, binding constants, and binding sites [44, 45]. Emission spectra alterations can yield insights about the structure and dynamics of a molecule [46, 47]. The main source of fluorescence in the uricase enzyme is attributed to the presence of tryptophan (Trp) and tyrosine (Tyr) residues. The Trp and Tyr residues of proteins have maximal fluorescence intensities of approximately 340 nm and 300 nm, respectively, when stimulated at 280 nm. Figure 7a, b, and c illustrates the fluorescence spectra of the uricase enzyme at different taurine concentrations. The fluorescence intensity of uricase showed a consistent reduction as the concentration of taurine increased, suggesting that taurine interacted with uricase and influenced its structure [48, 49].

Fig. 7
figure 7

Emission spectra of uricase in the presence of various taurine concentrations of (λexc = 280 nm) at a: 40 ˚C, b: 55 ˚C and c: 60 ˚C

Fluorescence quenching mechanism

Fluorescence quenching refers to the reduction in the efficiency of fluorescence emission by a fluorophore due to various chemical interactions, including excited-state reactions, formation of ground-state complexes, energy transfer, and collisional quenching. The methods of quenching can be categorized as either dynamic quenching or static quenching. These categories are differentiated by their respective dependencies on temperature and viscosity, or preferable by measuring their lifetimes. As the temperature increases, the diffusion coefficients also increase. Consequently, the dynamic quenching constants will also increase with rising temperature. In contrast, the rise in temperature is expected to lead to a reduction in the stability of complexes. Therefore, it is anticipated that the values of the static quenching constants will be lower [50]. The well-known Stern–Volmer equation was used to confirm the mechanism as follows [51]:

$$\frac{{F}_{0}}{F}=1+{K}_{q}{{\uptau }}_{0}\left[Q\right]=1+{K}_{SV}\left[Q\right]$$

F0 and F denote the fluorescence intensity when there is no quencher and when the quencher is present, respectively; [52] represents the concentration of the quencher, while τ0 denotes the fluorescence lifetime when there is no quencher present. The value of τ0 is always 10−8 seconds [53]. Kq represents the rate at which the biological macromolecule is quenched, whereas KSV denotes the Stern-Volmer quenching constant. KSV is also denoted by the following equation:

$${K}_{SV}={K}_{q}{{\uptau }}_{0}$$

This study aimed to investigate the quenching mechanism of uricase by obtaining fluorescence quenching spectra at three distinct temperatures (40, 55, and 60 °C) in the presence of various taurine concentrations. Figure 8 displays the Stern-Volmer plots depicting the quenching of uricase fluorescence by taurine. The computed KSV and Kq values are presented in Table 3. A linear Stern-Volmer plot typically suggests the presence of a single class of fluorophore, which is equally susceptible to quenching by the quencher. The assessment of static and dynamic quenching can be determined by observing the impact of temperature and viscosity. Another, more preferable approach is to conduct studies on fluorescence lifetime [22]. Increasing temperatures cause diffusion to occur more quickly, leading to an increase in collisional quenching. Additionally, higher temperatures typically cause loosely bound complexes to separate, resulting in a decrease in static quenching. Dynamic quenching is indicated by the rise in KSV value with increasing temperature. The Kq value was < 2.0 × 1010 L mol−1 s−1. Furthermore, the observation suggests that the decrease in uricase fluorescence caused by taurine is a form of static quenching, as seen in Fig. 7a, b, and c. The impact of dynamic quenching is widely acknowledged to affect only the excited state of the fluorophore, whereas static quenching alters the absorption spectrum of the fluorophore. This study showed that the quenching process occurs through static quenching, which is triggered by the formation of the uricase-taurine complex in its ground state. Osmolytes tend to increase the stability of proteins’ natural structure. In conclusion, these findings suggest that the interaction between the protein and the quencher is a complex one, which is more consistent with the static quenching process rather than dynamic collision quenching.

Fig. 8
figure 8

Stern–Volmer plots for the quenching of uricase by taurine at 40, 55, and 60 ˚C

Table 3 Stern-Volmer constants for the interaction of uricase with taurine at different temperatures

Calculation of binding parameters

The equation provided can be used to compute the binding constant (Kb) and the number of binding sites (n) in the static quenching process, where small molecules independently bind to a group of identical sites on a macromolecule.


where [52] is the concentration of the quencher; F0 and F are the fluorescence intensity in the absence and presence of the quencher, respectively; Ka is the binding constant, and n is the number of binding sites per uricase molecule. By plotting the logarithm of (F0 − F)/F against the logarithm of [52] (as depicted in Fig. 9), one may determine the values of n and Ka [54]. The results are succinctly presented in Table 4.

Fig. 9
figure 9

Double-log plots of the taurine quenching effect on uricase enzyme fluorescence at different temperatures

Significantly, the values of n at the experimental temperatures were nearly identical to 1, suggesting that there is just one binding site in uricase for taurine. Furthermore, the increase in temperature resulted in a decrease in the Ka value, indicating a decrease in the stability of the uricase-taurine complex. Therefore, the binding process was characterized as an exothermic reaction.

Table 4 Binging and thermodynamic parameters of UOX-Taurine interaction at three different temperatures

Thermodynamic parameters determination

The binding constant is influenced by temperature, indicating that the thermodynamic process is responsible for the creation of the uricase-taurine complex. An analysis of this dependency was made to reveal the interacting forces between taurine and uricase. Four non-covalent-type interactions attribute a major part of the bonding between ligands to themselves. The hydrophobic force, van der Waals forces, hydrogen bonding, and electrostatic forces are some of these interactions [22]. The major forces were determined based on the sign and magnitude of the thermodynamic parameters. The Vant Hoff equation can be used to approximate the value of ΔH and ΔS if they are only slightly variable and do not change much over the studied temperature range.

$$\begin{array}{ll}Ln{K}_{a}=-\frac{\varDelta H}{RT}+\frac{\varDelta S}{R}\\\varDelta G=-RTLn{k}_{a}=\varDelta H-T\varDelta S\end{array}$$

where Ka and T refer to the binding constant and temperature, respectively, and R is the gas constant (8.314 J mol−1 K−1). The ΔG values were determined and the obtained results were shown in Table 4. Also, the ΔH and ΔS values were calculated by graphing \(Ln{k}_{a}\) vs. \(\frac{1}{T}\) (Fig. 10).

Fig. 10
figure 10

plot lnKa against 1/T for interaction taurine with uricase at 313, 328, and 333 K

Ross and Subramanian have provided a comprehensive summary of the forces that govern the interactions between proteins and ligands, allowing for the identification of many forms of binding with distinct interactions. If the enthalpy change (ΔH) is negative and the entropy change (ΔS) is negative, then the binding reaction is mostly governed by van der Waals forces and hydrogen bond interactions. When the change in enthalpy (ΔH) is greater than zero and the change in entropy (ΔS) is also greater than zero, it indicates that hydrophobic interactions are the prevailing factor. If the change in enthalpy (ΔH) is negative and the change in entropy (ΔS) is positive, then the primary driving factors are electrostatic interactions. Enthalpy change (ΔH) can be considered constant when the temperature range is minimal. The utilization of enthalpy change (ΔH) and entropy change (ΔS) can serve to validate the binding modalities [55]. The interaction and forces governing the relationship between taurine and uricase can be described as follows: The negative values of ΔS and ΔH indicate that the van der Waals force and the hydrogen bonding interactions were the main factors influencing in the interaction of taurine with UOX. Furthermore, the negative ΔG values suggest that the binding process occurred spontaneously.

ANS fluorescence

The fluorescent probe ANS is employed for the purpose of identifying protein binding sites. The present investigation involved the utilization of spectroscopy techniques under controlled conditions, specifically at a temperature of 28 °C. The excitation wavelength employed was 370 nm, whereas the emission wavelength range spanned from 400 to 600 nm. Figure 11 displays the fluorescence emission curves pertaining to the interaction between ANS and both treated and untreated uricase enzymes. After the binding of ANS to the hydrophobic surface of the treated enzyme, there was a significant increase in fluorescence intensity. The results suggest that there has been a modification in the composition and conformation of the enzyme’s tertiary and quaternary structures, which has led to the placement of hydrophobic residues on the enzyme surface. This may allow for the establishment and arrangement of more stable hydrophobic regions or clefts by increasingly shielding them from a polar environment, thereby driving equilibrium towards the formation of a more stable protein complex [23, 56,57,58,59,60].

Fig. 11
figure 11

Fluorescence study of uricase in the absence and presence of taurine. Extrinsic fluorescence measurement of uricase using ANS. ANS (dashed line), ANS along with free uricase (dotted line), along with taurine treated uricase (solid line); The spectrums were taken by excitation at 280 nm

Molecular dynamics simulation results

MD simulations can be a valuable mean to gaining insight into the conformational properties of the enzyme [40, 61]. Our obtained theoretical results are valuable to improve the stability and performance development of UOX. MD simulations for native and taurine-treated enzymes were done through MD simulation package GROMACS 5.1.4 which adopts the CHARRM27 force field parameter for energy. To understand the enzyme activity, it is necessary to examine the protein structure in the presence of taurine.

RMSD and secondary structure analysis

The root mean square deviation (RMSD) and the root mean square fluctuation (RMSF) of Cα atoms of two simulated enzymes were determined and illustrated in Fig. 12a and b.

Fig. 12
figure 12

The root mean square deviation of Cα atoms (Cα-RMSD) and the root mean square fluctuation of Cα atoms (Cα-RMSF) of uricase in the absence (blue line) and presence of taurine (pink line)

The RMSD values reached saturation very fast, after about 2.3 ns, after which the variations reached the minimum values (see Fig. 12a). This evidence suggests that each system is simulated for a sufficient time to reach equilibrium and that the simulation time is enough to cover all possible states. At the end of the simulation, the RMSD values of naked and treated enzymes reached ~ 0.310 nm and ~ 0.144 nm, respectively. The lower Cα-RMSD value of the enzyme in the presence of taurine compared to the free enzyme indicates the stability of the UOX structure in the binary system.

The structural flexibility of both the naked and taurine-treated enzymes was assessed by calculating the Cα-RMSF (refer to Fig. 12b). In general, the pure water system exhibited greater fluctuations in residues compared to the binary system. Regions with higher flexibility were characterized by higher RMSF values. Consequently, the lower RMSF values observed in the binary system confirm that treating UOX with taurine enhances the stability of the enzyme.

The information about the local secondary structure changes has been obtained from the MD simulation. The UOX secondary structure during the 40 ns of MD simulations was investigated to examine the changes in the content of α-helix and β-sheet structures in the enzyme (See Fig. 13a and b).

Fig. 13
figure 13

Secondary structure assignment of the protein as a function of time for a: UOX and b: UOX/taurine for 40 ns

According to these results, the higher β-sheet and α-helix content of taurine-treated enzyme suggest that taurine may be responsible for the stability and higher activity of the enzyme at elevated temperatures and may preserve the enzyme’s active site integrity.

Table 5 reports the percentage of the secondary structures for the free and treated enzymes during the MD simulation.

Table 5 Secondary structure of UOX in the two simulated systems

The α-helix and β-sheets contents of the taurine-UOX enzyme (25.12% and 37.46%, respectively) increased compared with those of the free enzyme (23.90% and 32.61%, respectively). Further, the obtained results indicate that the β-bridge percentage in the treated enzyme has been higher than that of the naked enzyme (see Table 5). The obtained MD simulation results provide more details about the taurine effects on the enzyme structure.

The enzyme stability is estimated by a number of interactions including hydrogen bonding, electrostatic, hydrophobic interactions, and van der Waals forces [62, 63]. Thus, the taurine-UOX interaction changes the secondary structure composition of the enzyme. The hydrogen bond between the carbonyl oxygen and the amide hydrogen atoms forms the α-helix and β-sheet structures; thus, changes in the secondary structure composition of UOX may enhances its intramolecular H-bonds so that the content of the α-helix and β-sheet increases. In general, the results showed that taurine causes changes in the composition of the enzyme and thus boosts the enzyme activity, which is consistent with the results of the experimental section. The research findings confirmed the significant impact of taurine on the UOX stability.

In Fig. 14, it can be observed that the enzyme’s conformation was changed from a random coil to a β-sheet and α-helix in the binary system. This change in conformation could potentially maintain the integrity of the enzyme’s active site, contributing to its increased activity and stability at elevated temperatures.

Fig. 14
figure 14

The structures of a: free UOX and b: treated-taurine UOX at the end of 50 ns simulation (active site residues have been indicated with red color)

SASA and hydrogen bonding

Solvent accessible surface area (SASA) is a significant parameter to analyze the interactions between protein and solvent [64]. In other words, the SASA values show the conformational stability of the enzyme [65, 66]. The SASA values for the 50 ns simulation trajectory have been plotted in Fig. 15.

Fig. 15
figure 15

The SASA value of enzyme in the presence and absence of taurine

It is shown in Fig. 15 that the SASA value of UOX-taurine has diminished, which indicates that in the presence of taurine, the intramolecular hydrogen bonding of UOX grows, thus its accessible surface area by solvent molecules drops. The lower SASA value shows the higher thermodynamic stability of the protein. Thus, it can be stated that the lower SASA value in the UOX-taurine system shows a more compact enzyme structure as well as higher thermodynamic stability of UOX in the binary system (see Fig. 15).

Figure 16a reveals that the number of H-bonds between enzyme and water molecules in the UOX-taurine system has been lower than that of the UOX system, which is related to the enhancement in the number of intra-protein hydrogen bonds in the binary system (see Fig. 16b). Comparison of Figs. 15 and 16a shows that the SASA value corresponds well with the number of H-bonds between the UOX enzyme and water molecules.

Fig. 16
figure 16

The number of hydrogen bonds a: between UOX and water molecules and b: intra-protein and c: between UOX and taurine molecules

The taurine molecule, due to having hydrogen and hydroxyl oxygen atoms, can participate in the formation of hydrogen bonds with the UOX side chains (Fig. 16c). Thus, another reason for reduced number of H-bonds between UOX and water molecules may be related to the formation of hydrogen bonds between the enzyme and taurine molecules. In general, these results imply that the main driving forces for stabilizing enzyme structure are sharing hydrogen bonds between enzyme- enzyme and enzyme- taurine.

Molecular docking results

Molecular docking is a method used to predict the binding behavior and affinity between specific molecules and their template proteins in noncovalent binding [67, 68]. The target-based approach like docking can only be applied to proteins whose 3D structures have been identified [52]. In this study, docking simulations were conducted to analyze the interaction between taurine and the active site of the UOX enzyme. A total of 100 different conformations of taurine were examined during the docking process, and the conformation with the lowest binding energy was selected for further analysis. The results, shown in Table 6, include information on hydrogen bond and van der Waals energy, approximated free energy of binding, final intermolecular energy, final total internal energy, electrostatic energy, estimated inhibition constant, and torsional free energy.

Table 6 Docking results with the interacting residues

The presence of polar and hydrophobic residues around taurine suggests that the binding process is primarily driven by electrostatic and hydrophobic interactions. The negative value of the free binding energy confirms that UOX automatically binds to taurine. The best-docked orientation of taurine with the UOX active site was illustrated in Fig. 17.

Fig. 17
figure 17

Molecular docking analysis results. Docking simulation shows that A: Taurine molecule is located within the hydrophobic cavity of UOX active site. The orange dotted lines indicate the formation of the hydrogen bond

Additionally, the molecular docking data predicted the formation of hydrogen bonds between taurine and the binding pocket (Fig. 17). Therefore, these docking results confirmed the observations from the quenching fluorescence. Additionally, the quenching fluorescence findings suggested that the hydrogen bonding network and van der Waals forces were the primary binding forces contributing to the UOX-taurine interaction.


The obtained results suggest that estimating the optimal conditions to obtain the maximum activity using RSM-based CCD can save cost and time. The best conditions for UOX activity were obtained at a temperature of 28 °C, a taurine concentration of 450 mM, and an incubation time of 25 min. In general, the uricase activity of 2.05 U/ml was obtained under optimal conditions. The investigation of kinetic and thermodynamic parameters was conducted to study the mechanism of protein inactivation, which is an essential parameter in thermal processes. The results of the experiment indicated that the untreated urate oxidase enzyme exhibited lower stability under elevated temperatures in comparison to the enzyme treated with taurine. The Km values for the urate oxidase enzyme were investigated in the presence and absence of taurine osmolyte, and a two-fold increase in Km for the taurine-treated enzyme was observed, suggesting that it had a lower affinity for the substrate compared to the untreated enzyme. The other word, a low Km value signifies a minimal quantity of substrate required to achieve enzyme saturation, indicating a higher affinity for the substrate. Furthermore, MD simulation results showed that changes in the secondary structure composition of uricase enzyme enhanced its intramolecular hydrogen bonds, increasing the contents of α-helix and β-sheet. The presence of hydrogen bonds between the enzyme and taurine molecules also decreased the dynamics of water molecules, which had a significant role in the protection of protein structure. The ANS fluorescence of the enzyme was investigated, and the increased fluorescence intensity of the taurine-uricase-ANS complex compared to uricase-ANS indicated greater stabilization of the enzyme structure in the presence of taurine. Finally, the molecular docking results with the negative free binding energy value confirm that UOX automatically binds to taurine.

Availability of data ang materials

Not applicable.



urate oxidase


central composite designs


Molecular dynamics


Root-mean-square deviation


Solvent-accessible surface area


  1. Badoei-Dalfard A, Sohrabi N, Karami Z, Sargazi G. Fabrication of an efficient and sensitive colorimetric biosensor based on Uricase/Th-MOF for uric acid sensing in biological samples. Biosens Bioelectron. 2019;141: 111420.

    Article  Google Scholar 

  2. Shaaban MI, Abdelmegeed E, Ali YM. Cloning, expression, and purification of recombinant uricase enzyme from Pseudomonas aeruginosa Ps43 using Escherichia coli. J Microbiol Biotechnol. 2015;25(6):887–92.

    Article  Google Scholar 

  3. Feig DI, Kang D-H, Johnson RJ. Uric acid and cardiovascular risk. N Engl J Med. 2008;359(17):1811–21.

    Article  Google Scholar 

  4. Oda M, Satta Y, Takenaka O, Takahata N. Loss of urate oxidase activity in hominoids and its evolutionary implications. Mol Biol Evol. 2002;19(5):640–53.

    Article  Google Scholar 

  5. Pui C-H, Mahmoud HH, Wiley JM, Woods GM, Leverger G, Camitta B, Hastings C, Blaney SM, Relling MV, Reaman GH. Recombinant urate oxidase for the prophylaxis or treatment of hyperuricemia in patients with Leukemia or Lymphoma. J Clin Oncol. 2001;19(3):697–704.

    Article  Google Scholar 

  6. Richette P, Perez-Ruiz F, Doherty M, Jansen TL, Nuki G, Pascual E, Punzi L, So AK, Bardin T. Improving cardiovascular and renal outcomes in gout: what should we target? Nat Rev Rheumatol. 2014;10(11):654–61.

    Article  Google Scholar 

  7. Colloc’h N, Hajji ME, Bachet B, l’Hermite G, Schiltz M, Prangé T, Castro B, Mornon J-P. Crystal structure of the protein drug urate oxidase-inhibitor complex at 2.05 Å resolution. Nat Struct Biol. 1997;4(11):947–52.

    Article  Google Scholar 

  8. Vogt B. Urate oxidase (rasburicase) for treatment of severe tophaceous gout. Nephrol Dialysis Transpl. 2005;20(2):431–3.

    Article  Google Scholar 

  9. Taherimehr Z, Zaboli M, Torkzadeh-Mahani M. New insight into the molecular mechanism of the trehalose effect on urate oxidase stability. J Biomol Struct Dynamics. 2022;40(4):1461–71.

    Article  Google Scholar 

  10. Borowitzka LJ. Glycerol and Other Carbohydrate Osmotic Effectors. In: Gilles R, Gilles-Baillien M, (eds) Transport Processes, Iono- and Osmoregulation. Proceedings in Life Sciences. Berlin, Heidelberg: Springer; 1985.

  11. Jamal S, Poddar NK, Singh LR, Dar TA, Rishi V, Ahmad F. Relationship between functional activity and protein stability in the presence of all classes of stabilizing osmolytes. FEBS J. 2009;276(20):6024–32.

    Article  Google Scholar 

  12. Yancey PH, Clark ME, Hand SC, Bowlus RD, Somero GN. Living with water stress: evolution of osmolyte systems. Science. 1982;217(4566):1214–22.

    Article  Google Scholar 

  13. Rani A, Venkatesu P. Changing relations between proteins and osmolytes: a choice of nature. Phys Chem Chem Phys. 2018;20(31):20315–33.

    Article  Google Scholar 

  14. Kumar R. Role of naturally occurring osmolytes in protein folding and stability. Arch Biochem Biophys. 2009;491(1–2):1–6.

    Article  Google Scholar 

  15. Khan T, Abdullah M, Toor TF, Almajhdi FN, Suleman M, Iqbal A, Ali L, Khan A, Waheed Y, Wei D-Q. Evaluation of the whole proteome of Achromobacter xylosoxidans to identify vaccine targets for mRNA and peptides-based Vaccine Designing against the emerging respiratory and Lung Cancer-causing Bacteria. Front Med. 2022;8: 825876.

    Article  Google Scholar 

  16. Fattahian Y, Riahi-Madvar A, Mirzaee R, Torkzadeh-Mahani M, Asadikaram G. Heterologous expression, purification and characterization of a peroxidase isolated from Lepidium draba. Protein J. 2017;36(6):461–71.

    Article  Google Scholar 

  17. Hongoh Y, Sasaki T, Ishikawa H. Cloning, sequence analysis and expression in Escherichia coli of the gene encoding a uricase from the yeast-like symbiont of the brown planthopper, Nilaparvata lugens. Insect Biochem Mol Biol. 2000;30(2):173–82.

    Article  Google Scholar 

  18. Lohrasbi-Nejad A, Torkzadeh‐Mahani M, Hosseinkhani S. Hydrophobin‐1 promotes thermostability of firefly luciferase. FEBS J. 2016;283(13):2494–507.

    Article  Google Scholar 

  19. Nemati R, Molakarimi M, Mohseni A, Taghdir M, Khalifeh K, Sajedi H. Thermostability of Ctenophore and Coelenterate Ca2+-Regulated Apo-photoproteins: a comparative study. ACS Chem Biol. 2021;16(8):1538–45.

    Article  Google Scholar 

  20. Saqib AA, Siddiqui KS. How to calculate thermostability of enzymes using a simple approach. Biochem Mol Biol Educ. 2018;46(4):398–402.

    Article  Google Scholar 

  21. Lineweaver H, Burk D. The determination of enzyme dissociation constants. J Am Chem Soc. 1934;56(3):658–66.

    Article  Google Scholar 

  22. Farhadian S, Shareghi B, Saboury AA, Momeni L. Counteraction of lactose on the thermal stability and activity of α-chymotrypsin: thermodynamic, kinetic and docking studies. RSC Adv. 2016;6(76):72201–12.

    Article  Google Scholar 

  23. Deshpande M, Sathe SK. Interactions with 8-Anilinonaphthalene‐1‐sulfonic acid (ANS) and Surface Hydrophobicity of Black Gram (Vigna mungo) Phaseolin. J Food Sci. 2018;83(7):1847–55.

    Article  Google Scholar 

  24. Antony N, Balachandran S, Mohanan P. Immobilization of diastase α-amylase on nano zinc oxide. Food Chem. 2016;211:624–30.

    Article  Google Scholar 

  25. Retailleau P, Colloc’h N, Vivarès D, Bonnete F, Castro B, El Hajji M, Mornon J-P, Monard G, Prangé T. Complexed and ligand-free high-resolution structures of urate oxidase (Uox) from aspergillus flavus: a reassignment of the active-site binding mode. Acta Crystallogr Sect D: Biol Crystallogr. 2004;60(3):453–62.

    Article  Google Scholar 

  26. Bahrami H, Zahedi M. Comparison of the effects of sucrose molecules on alcohol dehydrogenase folding with those of sorbitol molecules on alcohol dehydrogenase folding using molecular dynamics simulation. J Iran Chem Soc. 2015;12(11):1973–82.

    Article  Google Scholar 

  27. Liu F-F, Dong X-Y, Sun Y. Molecular mechanism for the effects of trehalose on β-hairpin folding revealed by molecular dynamics simulation. J Mol Graph Model. 2008;27(4):421–9.

    Article  Google Scholar 

  28. Mehrnejad F, Ghahremanpour MM, Khadem-Maaref M, Doustdar F. Effects of osmolytes on the helical conformation of model peptide: molecular dynamics simulation. J Chem Phys. 2011;134(3):01B620.

    Article  Google Scholar 

  29. Zaboli M, Saeidnia F, Zaboli M, Torkzadeh-Mahani M. Stabilization of recombinant d-Lactate dehydrogenase enzyme with trehalose: response surface methodology and molecular dynamics simulation study. Process Biochem. 2021;101:26–35.

    Article  Google Scholar 

  30. Berendsen HJ, van Postma Jv WF, DiNola A, Haak J. Molecular dynamics with coupling to an external bath. J Chem Phys. 1984;81(8):3684–90.

    Article  Google Scholar 

  31. Bussi G, Donadio D, Parrinello M. Canonical sampling through velocity rescaling. J Chem Phys. 2007;126(1):014101.

    Article  Google Scholar 

  32. 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(12):10089–92.

    Article  Google Scholar 

  33. Hess B, Bekker H, Berendsen HJ, Fraaije JG. LINCS: a linear constraint solver for molecular simulations. J Comput Chem. 1997;18(12):1463–72.

    Article  Google Scholar 

  34. Monroe JI, Shirts MR. Converging free energies of binding in cucurbit [7] uril and octa-acid host–guest systems from SAMPL4 using expanded ensemble simulations. J Comput Aided Mol Des. 2014;28(4):401–15.

    Article  Google Scholar 

  35. Nurbaiti S, Martoprawiro MA, Akhmaloka, Hertadi R. The role of electrostatic interactions on klentaq1 insight for domain separation. Bioinform Biol Insights. 2012;6:BBI (S9390).

    Article  Google Scholar 

  36. Gao X, Liu Z, Cui W, Zhou L, Tian Y, Zhou Z. Enhanced thermal stability and hydrolytic ability of Bacillus subtilis aminopeptidase by removing the thermal sensitive domain in the non-catalytic region. PLoS ONE. 2014;9(3): e92357.

    Article  Google Scholar 

  37. Kabsch W, Sander C. Dictionary of protein secondary structure: pattern recognition of hydrogen-bonded and geometrical features. Biopolymers. 1983;22(12):2577–637.

    Article  Google Scholar 

  38. Morris GM, Goodsell DS, Halliday RS, Huey R, Hart WE, Belew RK, Olson AJ. Automated docking using a lamarckian genetic algorithm and an empirical binding free energy function. J Comput Chem. 1998;19(14):1639–62.

    Article  Google Scholar 

  39. Wallace AC, Laskowski RA, Thornton JM. LIGPLOT: a program to generate schematic diagrams of protein-ligand interactions. Protein Eng Des Selection. 1995;8(2):127–34.

    Article  Google Scholar 

  40. Torkzadeh-Mahani M, Zaboli M, Barani M, Torkzadeh‐Mahani M. A combined theoretical and experimental study to improve the thermal stability of recombinant D‐lactate dehydrogenase immobilized on a novel superparamagnetic Fe3O4NPs@ metal–organic framework. Appl Organomet Chem. 2020;34(5):e5581.

    Article  Google Scholar 

  41. Yamamoto K, Kojima Y, Kikuchi T, Shigyo T, Sugihara K, Takashio M, Emi S. Nucleotide sequence of the uricase gene from Bacillus sp. TB-90. J Biochem. 1996;119(1):80–4.

    Article  Google Scholar 

  42. Najavand S, Habibnejad M, Amani-Ghadim A, Rahimizadeh P, Pazhang M. Optimized immobilization of endoglucanase Cel9A onto glutaraldehyde activated chitosan nanoparticles by response surface methodology: the study of kinetic behaviors. Biotechnol Prog. 2020;36(3):e2960.

    Article  Google Scholar 

  43. Sant’Anna V, Cladera-Olivera F, Brandelli A. Kinetic and thermodynamic study of thermal inactivation of the antimicrobial peptide P34 in milk. Food Chem. 2012;130(1):84–9.

    Article  Google Scholar 

  44. Asgharzadeh S, Shareghi B, Farhadian S, Tirgir F. Effect of free L-cysteine on the structure and function of α-chymotrypsin. J Mol Liq. 2019;280:79–86.

    Article  Google Scholar 

  45. Raeessi-Babaheydari E, Farhadian S, Shareghi B. Evaluation of interaction between citrus flavonoid, naringenin, and pepsin using spectroscopic analysis and docking simulation. J Mol Liq. 2021;339: 116763.

    Article  Google Scholar 

  46. Teng Y, Liu R, Li C, Xia Q, Zhang P. The interaction between 4-aminoantipyrine and bovine serum albumin: multiple spectroscopic and molecular docking investigations. J Hazard Mater. 2011;190(1–3):574–81.

    Article  Google Scholar 

  47. Hu Y-J, Liu Y, Zhang L-X, Zhao R-M, Qu S-S. Studies of interaction between colchicine and bovine serum albumin by fluorescence quenching method. J Mol Struct. 2005;750(1–3):174–8.

    Article  Google Scholar 

  48. Bi S, Ding L, Tian Y, Song D, Zhou X, Liu X, Zhang H. Investigation of the interaction between flavonoids and human serum albumin. J Mol Struct. 2004;703(1–3):37–45.

    Article  Google Scholar 

  49. Zeng H-j, You J, Liang H-l, Qi T, Yang R, Qu L. -b. investigation on the binding interaction between silybin and pepsin by spectral and molecular docking. Int J Biol Macromol. 2014;67:105–11.

    Article  Google Scholar 

  50. Zhang Y-Z, Zhou B, Zhang X-P, Huang P, Li C-H, Liu Y. Interaction of malachite green with bovine serum albumin: determination of the binding mechanism and binding site by spectroscopic methods. J Hazard Mater. 2009;163(2–3):1345–52.

    Article  Google Scholar 

  51. Soares S, Mateus N, De Freitas V. Interaction of different polyphenols with bovine serum albumin (BSA) and human salivary α-amylase (HSA) by fluorescence quenching. J Agric Food Chem. 2007;55(16):6726–35.

    Article  Google Scholar 

  52. Zhao M, Chang H-T, Zhou Q, Zeng T, Shih C-S, Liu Z-P, Chen L, Wei D-Q. Predicting protein-ligand interactions based on chemical preference features with its application to new D-amino acid oxidase inhibitor discovery. Curr Pharm Design. 2014;20(32):5202–11.

    Article  Google Scholar 

  53. Lakowicz JR, Weber G. Quenching of fluorescence by oxygen. Probe for structural fluctuations in macromolecules. Biochemistry. 1973;12(21):4161–70.

    Article  Google Scholar 

  54. Ghalandari B, Divsalar A, Saboury AA, Haertlé T, Parivar K, Bazl R, Eslami-Moghadam M, Amanlou M. Spectroscopic and theoretical investigation of oxali–palladium interactions with β-lactoglobulin. Spectrochim Acta Part A Mol Biomol Spectrosc. 2014;118:1038–46.

    Article  Google Scholar 

  55. Soong Y-Y, Barlow PJ. Quantification of gallic acid and ellagic acid from longan (Dimocarpus longan Lour.) seed and mango (Mangifera indica L.) kernel and their effects on antioxidant activity. Food Chem. 2006;97(3):524–30.

    Article  Google Scholar 

  56. Gasymov OK, Glasgow BJ. ANS fluorescence: potential to augment the identification of the external binding sites of proteins. Biochimica et Biophysica Acta (BBA)-Proteins and proteomics. 2007;1774(3):403–11.

    Article  Google Scholar 

  57. Mortazavi M, Hosseinkhani S, Khajeh K, Ranjbar B, Emamzadeh AR. Spectroscopic and functional characterization of Lampyris turkestanicus luciferase: a comparative study. Acta Biochim Biophys Sin. 2008;40(5):365–74.

    Article  Google Scholar 

  58. Nakai S. Measurement of protein hydrophobicity. Curr Protocols Food Analytical Chemist. 2003;9(1):B5-2.

    Google Scholar 

  59. Niesen FH, Berglund H, Vedadi M. The use of differential scanning fluorimetry to detect ligand interactions that promote protein stability. Nat Protoc. 2007;2(9):2212–21.

    Article  Google Scholar 

  60. Poklar N, Lah J, Salobir M, Macek P, Vesnaver G. pH and temperature-induced molten globule-like denatured states of equinatoxin II: a study by UV-melting, DSC, far-and near-UV CD spectroscopy, and ANS fluorescence. Biochemistry. 1997;36(47):14345–52.

    Article  Google Scholar 

  61. Ghanbari-Ardestani S, Khojasteh-Band S, Zaboli M, Hassani Z, Mortezavi M, Mahani M, Torkzadeh-Mahani M. The effect of different percentages of triethanolammonium butyrate ionic liquid on the structure and activity of urate oxidase: molecular docking, molecular dynamics simulation, and experimental study. J Mol Liq. 2019;292: 111318.

    Article  Google Scholar 

  62. Pikkemaat MG, Linssen AB, Berendsen HJ, Janssen DB. Molecular dynamics simulations as a tool for improving protein stability. Protein Eng. 2002;15(3):185–92.

    Article  Google Scholar 

  63. Zaboli M, Raissi H, Zaboli M, Farzad F, Torkzadeh-Mahani M. Stabilization of D-lactate dehydrogenase diagnostic enzyme via immobilization on pristine and carboxyl-functionalized carbon nanotubes, a combined experimental and molecular dynamics simulation study. Arch Biochem Biophys. 2019;661:178–86.

    Article  Google Scholar 

  64. Saxena R. Effect of processing on bovine milk quality attributes for human consumption. Canada: McGill University; 2018.

    Google Scholar 

  65. Chan MK, Mukund S, Kletzin A, Adams M, Rees DC. Structure of a hyperthermophilic tungstopterin enzyme, aldehyde ferredoxin oxidoreductase. Science. 1995;267(5203):1463–9.

    Article  Google Scholar 

  66. Lee Y, Lee J-j, Kim S, Lee S-C, Han J, Heu W, Park K, Kim HJ, Cheong H-K, Kim D. Dissecting the critical factors for thermodynamic stability of modular proteins using molecular modeling approach. PLoS ONE. 2014;9(5): e98243.

    Article  Google Scholar 

  67. Asemi-Esfahani Z, Shareghi B, Farhadian S, Momeni L. Effect of Naphthol yellow S as a food dye on the lysozyme structure and its mechanisms of action. J Mol Liq. 2021;332: 115846.

    Article  Google Scholar 

  68. Eslami-Farsani R, Shareghi B, Farhadian S, Momeni L. Insight into the binding of glycerol with myoglobin spectroscopic and MD simulation approach. Int J Biol Macromol. 2020;159:433–43.

    Article  Google Scholar 

Download references


We would like to thank the Graduate University of Advanced Technology for the support of this work.


This research received no external funding.

Author information

Authors and Affiliations



P. Sh. and M. Z. wrote the main manuscript text and contributed equally to this work. All authors helped in performing and drafting the manuscript. The authors read and approved the final manuscript.

Corresponding author

Correspondence to Masoud Torkzadeh-Mahani.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Figure S1.

The root mean square deviation of Cα atoms (Cα-RMSD) of uricase in the absence (A) and presence of taurine (B) during the simulation (first run, second run and third run). Figure S2. The SASA value of uricase in the absence (A) and presence of taurine (B) during the simulation (first run, second run and third run). Figure S3. The number of hydrogen bonds between water molecules and uricase in the absence (A) and presence of taurine (B) during the simulation (first run, second run and third run). Figure S4. The number of hydrogen bonds within the protein in the absence (A) and presence of taurine (B) during the simulation (first run, second run and third run). Figure S5. The number of hydrogen bonds between UOX and taurine molecules during the simulation (first run, second run and third run).

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Shahmoradipour, P., Zaboli, M. & Torkzadeh-Mahani, M. Exploring the impact of taurine on the biochemical properties of urate oxidase: response surface methodology and molecular dynamics simulation. J Biol Eng 18, 10 (2024).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: