Gurprit Sekhon1, Ranvir Singh1. 1. Department cum National Center for Human Genome Studies and Research, Panjab University, Chandigarh, 160014, India.
Abstract
Background: Human aldose reductase (hAR) is the first and rate-limiting enzyme of the polyol pathway. For the development of secondary complications of diabetes in chronic hyperglycemic conditions, one of the critical factors is the increased flux of glucose through the polyol pathway. Due to this clinical implication, hAR attracted considerable attention from the drug discovery perspective. In spite of extensive characterization in the context of biochemical and structural aspects, we know very little about the unfolding behavior of hAR. This study reports equilibrium unfolding studies of hAR. Methods: We carried out thermal denaturation and chemical-induced equilibrium unfolding studies of hAR monitored by circular dichroism and fluorescence spectroscopy. Results: Thermal denaturation studies presented a classical picture of two-state unfolding from native to the denatured state. The data was used to derive thermodynamic parameters and study the thermostability of hAR. Chemical induced equilibrium unfolding studies led us to discover an intermediate state, which gets populated at 3.5-4.0 M and 0.7-2.0 M of urea and GuHCl, respectively. Thermodynamic parameters derived from chemical-induced unfolding are in agreement with those obtained from thermal denaturation of hAR. Conclusion: This study revealed that aldose reductase unfolds from native to the unfolded state via an intermediate. Assessment of the thermodynamic stability of native, intermediate, and unfolded states shows that significant energy barriers separate these states, which ensures the cooperativity of unfolding. As hAR functions in cells that are under osmotic and oxidative stress, these in vitro findings may have implications for its native conformation under the physiological state. Copyright:
Background: Humanaldose reductase (hAR) is the first and rate-limiting enzyme of the polyol pathway. For the development of secondary complications of diabetes in chronic hyperglycemic conditions, one of the critical factors is the increased flux of glucose through the polyol pathway. Due to this clinical implication, hAR attracted considerable attention from the drug discovery perspective. In spite of extensive characterization in the context of biochemical and structural aspects, we know very little about the unfolding behavior of hAR. This study reports equilibrium unfolding studies of hAR. Methods: We carried out thermal denaturation and chemical-induced equilibrium unfolding studies of hAR monitored by circular dichroism and fluorescence spectroscopy. Results: Thermal denaturation studies presented a classical picture of two-state unfolding from native to the denatured state. The data was used to derive thermodynamic parameters and study the thermostability of hAR. Chemical induced equilibrium unfolding studies led us to discover an intermediate state, which gets populated at 3.5-4.0 M and 0.7-2.0 M of urea and GuHCl, respectively. Thermodynamic parameters derived from chemical-induced unfolding are in agreement with those obtained from thermal denaturation of hAR. Conclusion: This study revealed that aldose reductase unfolds from native to the unfolded state via an intermediate. Assessment of the thermodynamic stability of native, intermediate, and unfolded states shows that significant energy barriers separate these states, which ensures the cooperativity of unfolding. As hAR functions in cells that are under osmotic and oxidative stress, these in vitro findings may have implications for its native conformation under the physiological state. Copyright:
Humanaldose reductase (hAR) (EC 1.1.1.21) is an NADPH-dependent oxidoreductase that belongs to the superfamily of aldo-keto reductases
[1]. hAR is the first and rate-limiting enzyme of the polyol pathway and converts glucose to sorbitol
[2]. The polyol pathway is up-regulated under hyperglycemic conditions, and a significant proportion of glucose gets fluxed through this pathway, which leads to the accumulation of sorbitol, consumption of NADPH, and redox imbalance of NADPH/NADP
+ ratio. All these factors have been linked with various tissue-based pathologies associated with secondary complications of diabetes mellitus
[3]. Due to its clinical importance from the perspective of developing potent inhibitors to prevent or delay the onset of secondary diabetic complications, hAR has been widely studied in almost every aspect
[4].Extensive information is available in the literature about the structure and function of hAR, mainly related to active site of hAR from high-resolution crystal structures with a number of potential inhibitors
[5], flexibility in the hAR binding site pocket
[6] and the thermodynamics of closing/opening of the specificity pocket within binding site pocket of hAR
[7]. There is little investigation related to the folding/unfolding mechanism of hAR.Protein folding is a fundamental process in all living systems. The linear sequence of amino acids encodes the information required for a polypeptide to fold. Despite extensive research and progress made over the years to understand this fundamental process, a complete understanding of protein folding mechanism(s) and the folding code remains elusive
[8]. Various models have been proposed to explain the mechanism underlying protein folding reaction from time to time
[9]. ‘New view folding model’ based on energy landscape theory (ELT) and ‘defined pathways model’ based on the concept of foldons are current alternatives to explain the protein folding mechanism.
[10].Among several hallmarks of protein folding are its spontaneity and cooperativity, along with fast reaction timescale. Usually, a small free energy change is required to shift the equilibrium between the folded and unfolded state ensemble. Complex cellular solvation environment plays a vital role in regulating the equilibrium between folded and unfolded ensemble in vivo, which are separated by small free energy barrier (~5 - 10 Kcal mol-1)
[11].Under physiological conditions, protein structure fluctuates among different native conformations separated by close free energy barriers
[12]. Since hAR activity leads to sorbitol accumulation, leading to osmotic stress, it seems to function under stress conditions which might perturb its native conformation ensemble. Here we report on thermal denaturation and chemical induced unfolding studies of hAR. Thermal denaturation revealed a simple two-state transition, whereas chemical induced unfolding led us to discover an intermediate state during hAR unfolding.
Methods
Materials
All chemicals were reagent grade and purchased from Sigma-Aldrich.
Protein purification
The hAR cDNA cloned into expression vector pET-15b (Novagen) was a kind gift from Dr. Alberto Podjarny (Department of Integrated Structural Biology, Institut de Génétique et de Biologie Moléculaire et Cellulaire, CNRS, INSERM, France). The plasmid, coding for a hexahistidine-tagged protein, was expressed into
E. coli strain BL21 (DE3) (Novagen). The cells with recombinant plasmid were grown with 100 μM ampicillin at 37°C to an OD
600 nm value of 0.7, and protein expression was induced by adding 1 mM IPTG. Cells were grown for a further 3 hours at 37°C. All further operations were carried out at 4°C unless otherwise stated. Cells were centrifuged, resuspended, and lysed by sonication. A Ni-NTA affinity column (GE Healthcare) was used for protein purification. The material used for the stationary phase for the column was Ni-Sepharose, and the flow rate of the column was adjusted to 0.5 ml min
-1. Imidazole and other salts were removed by repeated dialysis in 50 mM potassium phosphate, pH-7 buffer containing 50 mM NaCl. Protein concentration was estimated using the molar extinction coefficient and absorbance reading at 280 nm. The histidine tag from recombinant protein was removed by thrombin (4 units of thrombin per mg of recombinant protein at room temperature for 3 hours). The cleaved protein was passed through the Ni-NTA column, and purified protein without tag was collected as flow-through. Enzyme activity was checked as per standard assay
[13]. We analyzed the homogeneity and molecular weight of hAR with and without histidine tags under denaturing conditions on 15% SDS-PAGE. Purified hAR was stored at -20°C for further studies.
Thermal denaturation
Thermal denaturation was carried out at a final concentration of 2.8 µM protein in 50 mM potassium phosphate buffer, pH 7.0, containing 50 mM NaCl and 0.1mM TCEP. The transition between 20–70 °C was followed using a far-UV circular dichroism (CD) signal at 222 nm by using a 0.1 cm path-length cuvette at a sampling rate of 1.0 °C min
-1 in a Jasco J-810 spectropolarimeter. After subtracting buffer blank, we have reported the change in ellipticity (millidegree) at 222 nm.
Chemical induced unfolding
To follow chemical induced unfolding, we prepared samples with 1.4 μM protein concentration in phosphate buffer (described earlier), containing different concentrations of GuHCl/urea. Samples were incubated for 12 hours to reach equilibrium at 25°C, after which no change in signal occurred either in fluorescence or CD spectra. Trp fluorescence (excitation at 295 nm and emission recorded between 300 nm to 400 nm) and ANS fluorescence (excitation at 370 nm and emission recorded between 400 nm to 600 nm) measurements were performed using a Hitachi F-7000 fluorescence spectrophotometer for GuHCl samples and Jasco J-815 spectropolarimeter for urea samples. Far-UV CD measurements were performed using Jasco J-810 spectropolarimeter for GuHCl samples and Jasco J-815 spectropolarimeter for urea samples. Quartz cuvette of 1 cm and 0.5 cm path-length were used for fluorescence and CD measurements, respectively. All measurements were done at 25 °C. We have reported spectra as ellipticity (millidegree) after baseline correction.
Heat capacity change (Δ
C
) calculations
We have calculated the value of Δ
C
for the unfolding of hAR from the change in accessible surface area (ΔASA) according to
Equation 1
[14]ProtSA webserver was used to calculate the change in accessible surface area from native to the unfolded conformation of hAR
[15].
Data analysis
GraphPad Prism version 7.04 for Windows (GraphPad Software, La Jolla, California) was used for the analysis of thermal denaturation and chemical-induced unfolding data based on two-state and three-state models respectively, as described in following sections.The Least-square analysis was used to fit the data to
Equation 2
[16].Where
A
and
A
are native and unfolded state baseline intercepts, respectively, and
b
and
b
are native and unfolded baseline slopes, respectively. Δ
H
is enthalpy change at melting temperature (
T
).
T is absolute temperature, and R is the gas constant.
Calculation of
ΔG values for the transition region
Signal for native (
Y
×
T) and unfolded baseline (
Y
×
T) for every point in transition region was calculated from
Equation 2. If
Y is signal for a particular point in transition region, then the fraction of unfolded protein (
F
) at this point is given by
Equation 3.The equilibrium constant (k) was calculated from the relative population of species using
Equation 4.Δ
G was calculated as a function of temperature using
Equation 5.
Thermal stability curve
The thermal stability curve of hAR was constructed based on
Equation 6–8
[17].Where Δ
H
and Δ
S
are enthalpy change and entropy change, respectively, at temperature
T with reference to
T
.
T
,
T
and
T
are the temperatures at which Δ
H, Δ
S, and Δ
G are zero, respectively. Δ
G
is the stabilization free energy of the native state relative to the unfolded state.The least-square analysis was used to fit the data to
Equation 9
[16].Where
A
,
A
, and
A
are the native, unfolded, and intermediate baseline intercepts, respectively, and b
N,
b
, and
b
are the native, unfolded, and intermediate baseline slopes, respectively. [
D] is denaturant concentration in a molar unit.
m
and
m
are denaturant gradient for native to intermediate and intermediate to unfolded state, respectively. Δ
G
and Δ
G
are stabilization free energy of intermediate state relative to the native and unfolded state, respectively.
Results
Thermal denaturation monitored by far-UV CD
Change in ellipticity at 222 nm fitted well based on a two-state model (
Figure 1A). This analysis gave values for Δ
H
and
T
, which along with the Δ
C
value calculated from
Equation 1 were used for non-linear regression of the transition region (±5 kJ mol
-1) to
Equation 8 (
Figure 1C). Values of Δ
H and Δ
S were calculated over an extended range of temperatures by using
Equation 7 and
Equation 8, respectively (
Figure 1D). The thermal stability curve is an extrapolation of transition region, assuming constant Δ
C
during the unfolding transition (
Figure 1E). The relationship between
T
,
T
and Δ
G (T
Δ
G
/Δ
C
is presented in
Figure 1F. Thermodynamic parameters obtained from the analysis of thermal denaturation data are listed in
Table 1. All raw data are available as
Underlying data
[18].
Figure 1.
Thermal denaturation studies of human aldose reductase (hAR) monitored by far-UV circular dichroism.
(
A) Change in ellipticity at 222 nm plotted as a function of temperature. (
B) The fraction of protein folded (black dots) and unfolded (cyan dots) plotted against temperature. (
C) The portion of the transition curve used in van't Hoff analysis. (
D) Plots of ΔH and ΔS as a function of temperature. (
E) Thermal stability curve of hAR. (
F) Triangular relationship among T
h, T
s, and ∆G
s. The explanation for T
g, T
g’, T
h, and T
s is given in the text. Dashed lines are extrapolations.The solid lines represent curves fitted to the unfolding transition; filled symbols represent data points from unfolding experiments; and red symbols represent outliers.
Table 1.
Thermodynamic parameters derived from equilibrium unfolding transitions of hAR during thermal denaturation and Chemical induced unfolding studies using different probes.
Chemical
unfolding
Probe
ΔG
(N-I)(KJ mol
-1)
ΔG
(I-U)(KJ mol
-1)
m
(N-I)(KJ mol
-1 M
-1)
m
(I-U)(KJ mol
-1 M
-1)
C
m(N-I)(M)
C
m(I-U)(M)
GuHCl [F
314]
9.47 ± 1.69
58.34 ± 1.14
48.34 ± 7.69
23.56 ± 0.46
0.16
2.48
GuHCl-ANS [F
480]
16.48 ± 0.54
57.26 ± 3.33
38.99 ± 1.15
22.38 ± 1.29
0.42
2.56
Urea [F
314]
30.88 ± 4.47
35.04 ± 4.16
13.11 ± 2.01
7.99 ± 0.84
2.35
4.49
Urea-ANS [F
480]
28.61 ± 2.39
38.6 ± 2.14
12.1 ± 1.13
8.95 ± 0.46
2.36
4.32
GuHCl [Θ
219]
13.68 ± 5.78
34.13 ± 3.75
32.17 ± 11.63
14.9 ± 1.58
0.43
2.29
Urea [Θ
221.8]
30.53 ± 6.16
24.43 ± 6.25
14.12 ± 3.03
5.99 ± 1.30
2.16
4.08
Thermal
denaturation
Probe
ΔG
u(KJ mol
-1)
T
g(K)
ΔC
p(KJ mol
-1 K
-1)
ΔH
g(KJ mol
-1)
ΔS
m(KJ mol
-1)
T
g`(K)
T
h(K)
T
s(K)
Temperature [Θ
222]
50.25
329 ± 0.01
17.57
779.20 ± 15.36
2.37
248
284.65
287.51
Thermal denaturation studies of human aldose reductase (hAR) monitored by far-UV circular dichroism.
(
A) Change in ellipticity at 222 nm plotted as a function of temperature. (
B) The fraction of protein folded (black dots) and unfolded (cyan dots) plotted against temperature. (
C) The portion of the transition curve used in van't Hoff analysis. (
D) Plots of ΔH and ΔS as a function of temperature. (
E) Thermal stability curve of hAR. (
F) Triangular relationship among T
h, T
s, and ∆G
s. The explanation for T
g, T
g’, T
h, and T
s is given in the text. Dashed lines are extrapolations.The solid lines represent curves fitted to the unfolding transition; filled symbols represent data points from unfolding experiments; and red symbols represent outliers.
Chemical induced unfolding monitored by fluorescence
There are six Trp residues in hAR, out of which four are part of the hydrophobic active site pocket in the core of the β-barrel, and two are buried in alpha-helices surrounding the barrel. Their fluorescence provided a global signal of change in the tertiary structure. ANS has been extensively used as a probe for non-native, partially unfolded conformations of the protein. The binding of ANS to hydrophobic regions results in a significant enhancement of ANS fluorescence and a pronounced blue-shift of the
λ
max
[19].Fluorescence emission profiles of hAR equilibrated with different concentrations of denaturants are presented in
Figure 2 (
Figure 2A and
Figure 3A for urea and GuHCl, respectively). A plot of
λ
max against denaturant concentration indicated a cooperative transition from native to unfolded state (
Figure 2B and
Figure 3B for urea and GuHCl, respectively). In the case of ANS fluorescence, a significant blue-shift of around 20 nm and 10 nm from the native to the intermediate state was observed for urea and GuHCl, respectively (
Figure 2B2 and
Figure 3B2 for urea and GuHCl, respectively).
Figure 2.
Urea induced unfolding studies of human aldose reductase (hAR)monitored by fluorescence.
(
A1) Trp fluorescence scans and (
A2) ANS fluorescence scans for all the samples. (
B1) λ
max (Trp fluorescence) and (
B2) λ
max (ANS fluorescence) against [urea]. (
C1) I
max (Trp fluorescence) and (
C2) I
max (ANS fluorescence) against [urea]. (
D1) I
295/314 (Trp fluorescence) and (
D2) I
370/480 (ANS fluorescence) against [urea]. (
E1) Trp fluorescence and (
E2) ANS fluorescence of samples in native (black), intermediate (cyan), and unfolded state (red). Solid lines represent curves fitted to the unfolding transitions; filled symbols represent data points from unfolding experiments; red symbols represent outliers.
Figure 3.
GuHCl-induced unfolding studies of human aldose reductase (hAR)monitored by fluorescence.
(
A1) Trp fluorescence scans and (
A2) ANS fluorescence scans. (
B1) λ
max (Trp fluorescence) and (
B2) λ
max (ANS fluorescence) against [GuHCl]. (
C1) I
max (Trp fluorescence) and (
C2) I
max (ANS fluorescence) against [GuHCl]. (
D1) I
295/314 (Trp fluorescence) and (
D2) I
370/480 (ANS fluorescence) against [GuHCl]. (
E1) Trp fluorescence and (
E2) ANS fluorescence of samples in native (black), intermediate (cyan), and unfolded state (red). Solid lines represent curves fitted to the unfolding transitions; filled symbols represent data points from unfolding experiments; red symbols represent outliers.
Urea induced unfolding studies of human aldose reductase (hAR)monitored by fluorescence.
(
A1) Trp fluorescence scans and (
A2) ANS fluorescence scans for all the samples. (
B1) λ
max (Trp fluorescence) and (
B2) λ
max (ANS fluorescence) against [urea]. (
C1) I
max (Trp fluorescence) and (
C2) I
max (ANS fluorescence) against [urea]. (
D1) I
295/314 (Trp fluorescence) and (
D2) I
370/480 (ANS fluorescence) against [urea]. (
E1) Trp fluorescence and (
E2) ANS fluorescence of samples in native (black), intermediate (cyan), and unfolded state (red). Solid lines represent curves fitted to the unfolding transitions; filled symbols represent data points from unfolding experiments; red symbols represent outliers.
GuHCl-induced unfolding studies of human aldose reductase (hAR)monitored by fluorescence.
(
A1) Trp fluorescence scans and (
A2) ANS fluorescence scans. (
B1) λ
max (Trp fluorescence) and (
B2) λ
max (ANS fluorescence) against [GuHCl]. (
C1) I
max (Trp fluorescence) and (
C2) I
max (ANS fluorescence) against [GuHCl]. (
D1) I
295/314 (Trp fluorescence) and (
D2) I
370/480 (ANS fluorescence) against [GuHCl]. (
E1) Trp fluorescence and (
E2) ANS fluorescence of samples in native (black), intermediate (cyan), and unfolded state (red). Solid lines represent curves fitted to the unfolding transitions; filled symbols represent data points from unfolding experiments; red symbols represent outliers.A plot of
I
max against denaturant concentration indicated the presence of an intermediate during unfolding transition (
Figure 2C and
Figure 3C for urea and GuHCl, respectively).
I
max in case of ANS fluorescence fits satisfactorily based on the three-state model (
Figure 2C2 and
Figure 3C2 for urea and GuHCl, respectively). For both urea and GuHCl induced unfolding, Trp fluorescence emission intensity at 314 nm fits satisfactorily to the three-state model. In the case of ANS fluorescence, both
I
and fluorescence emission intensity at 480 nm fit equally well based on three-state model. Thus, Trp fluorescence emission intensity at 314 nm and ANS fluorescence emission intensity at 480 nm were analyzed based on the three-state model to evaluate the thermodynamic stability of hAR (
Figure 2D and
Figure 3D for urea and GuHCl, respectively). The thermodynamic parameters obtained from the fittings are listed in
Table 1. Trp and ANS fluorescence demonstrate the presence of an intermediate state populated at 3.5-4.0 M and 0.7-2 M urea and GuHCl concentration, respectively, apart from the native and unfolded states (
Figure 2E and
Figure 3E for urea and GuHCl, respectively).
Chemical induced unfolding monitored by far-UV CD
Unfolding profiles of hAR equilibrated at different denaturant concentrations monitored by far-UV CD are presented in
Figure 4A1 and 4A2 for GuHCl and urea, respectively. The thermodynamic stability of hAR was determined based on the three-state model by plotting the change in ellipticity at 219/222 nm as a function of denaturant concentration (
Figure 4B1 and 4B2 for GuHCl and urea, respectively). The transition determined by far-UV CD detected intermediate state at similar concentrations of denaturant as interrogated by fluorescence spectroscopy. Far-UV CD profiles clearly distinguish three states (
Figure 4C1 and 4C2 for GuHCl and urea, respectively). Thermodynamic parameters derived from far-UV CD data are listed in
Table 1.
Figure 4.
GuHCl/urea induced unfolding studies of human aldose reductase (hAR)monitored by far-UV circular dichroism (CD).
(
A1) far-UV CD scans recorded for GuHCl and (
A2) far-UV CD scans recorded for urea for all the samples. (
B1) change in ellipticity at 219 nm against [GuHCl] and (
B2) change in ellipticity at 221.8 nm against [urea]. (
C) CD spectra of samples in native (black), intermediate (cyan), and unfolded state (red) for (
C1) GuHCl and (
C2) urea, respectively. The solid lines represent curves fitted to the unfolding transitions; filled symbols represent data points from unfolding experiments; red symbols represent outliers in data fitting.
GuHCl/urea induced unfolding studies of human aldose reductase (hAR)monitored by far-UV circular dichroism (CD).
(
A1) far-UV CD scans recorded for GuHCl and (
A2) far-UV CD scans recorded for urea for all the samples. (
B1) change in ellipticity at 219 nm against [GuHCl] and (
B2) change in ellipticity at 221.8 nm against [urea]. (
C) CD spectra of samples in native (black), intermediate (cyan), and unfolded state (red) for (
C1) GuHCl and (
C2) urea, respectively. The solid lines represent curves fitted to the unfolding transitions; filled symbols represent data points from unfolding experiments; red symbols represent outliers in data fitting.
Discussion
The intermediate state with enhanced ANS fluorescence and significant blue shift of λ
max pointed to an intermediate state with ‘molten like’ nature during the chemical-induced unfolding of hAR. Far-UV CD studies strongly suggest that the intermediate state retains significant secondary structure during urea- and GuHCl-induced unfolding.During chemical induced unfolding, hAR unfolds through an intermediate state, which is absent during thermal denaturation. Moderate concentration of denaturant is known to stabilize native or intermediate state
[20]. The absence of such a stabilizing agent may be the reason that we could not detect the intermediate state during thermal denaturation.In all three probes used in studying unfolding, the value of Δ
G
(
obtained is ~30 kJ mol
-1 and ~15 kJ mol
-1 for urea- and GuHCl-induced unfolding, respectively while a Δ
G
of ~70 kJ mol
-1 is almost same for both denaturants (
Table 1). Thus, while urea seems to stabilize the native state with respect to the intermediate state, GuHCl seems to stabilize the intermediate state with respect to the native state.Small molecules change the free energy landscape of protein upon binding by selectively stabilizing native or intermediate/unfolded state
[21]. The difference in the value of Δ
G
obtained from thermal denaturation and chemical-induced unfolding is ~20 kJ mol
-1 (
Table 1), which is most likely due to free energy of stabilization and destabilization by urea and GuHCl, respectively.Values of Δ
G
obtained from the analysis of ANS fluorescence data are 16.48 and 28.61 kJ mol
-1 for GuHCl and urea, respectively (
Table 1), which indicates that a steep energy barrier does not separate intermediate state from the native state. Values for Δ
G
obtained from ANS fluorescence are 57.26 and 38.6 kJ mol
-1 for GuHCl and urea, respectively (
Table 1), which indicates that the intermediate state is separated from unfolded state by a high energy barrier. Thus, the intermediate state of hAR is close to its native state, which makes it functionally more relevant. These studies are essential to access the effect of other physiological relevant molecules on hAR stability
[22]. Protein stability of hAR in the cellular environment under hyperglycemia-induced cellular stress conditions has to be studied to aid rational drug design
[23].In summary, equilibrium unfolding studies of hAR have led us to discover that hAR unfolds through an intermediate state, which is close to the native state, and might have physiological relevance under hyperglycemic conditions in diabetes.
Data availability
Underlying data
Figshare: data_f1000_hAR_unfolding.zip.
https://doi.org/10.6084/m9.figshare.8001998.v1
[18].This project contains raw data for thermal denaturation and chemical induced unfolding studies on humanaldose reductase.Data are available under the terms of the
Creative Commons Attribution 4.0 International license (CC-BY 4.0).I am sorry for that I can not approve the revised versioon since the authors did not address most of my concerns raised during last round of review, particular of the change in enzymatic activity, the validity of the calculations and the effect of pathological stresses.I confirm that I have read this submission and believe that I have an appropriate level of expertise to confirm that it is of an acceptable scientific standard, however I have significant reservations, as outlined above.The authors seem to have carefully addressed all points raised.I confirm that I have read this submission and believe that I have an appropriate level of expertise to confirm that it is of an acceptable scientific standard.In this research, the authors studied the equilibrium unfolding pathway of humanaldose reductase (hAR) by spectroscopic methods. Although this pure unfolding study may provide information of how hAR unfolds, the underlying physiological relevance is unclear. Furthermore, the present study needs to be improved in the quality of the data and interpretation.During my preparation of the comments, Professor Uversky and Professor Mueller have published their excellent comments and suggestions. I agree with those points raised in their reports. Below are some additional comments that could be considered by the authors to improve their work.Although the authors have mentioned the importance of hAR in the Abstract and Introduction, especially for that it is an enzyme functions under osmotic stresses, it is unclear for the physiological relevance of the information obtained here using heat and chemical denaturants. I suggest the authors considered the physiologically relevant factors such as osmotic stress in their revised or future work.hAR is an enzyme. I suggest the authors monitor the activity of the enzyme during unfolding to provide information regarding the local structural changes around the active site.I am not convinced by the calculations of the thermodynamic parameters. The equations listed in the manuscript are derived from reverse folding processes. However, the authors did not show that the folding of hAR was reversible or not.The authors concluded that the thermal unfolding of hAR was a two-state process by far-UV CD spectra analysis. However, only one method/probe is difficult for the definition of a two-state transition. It’s better to study the transition with two or more probes and pure two-state transition will give similar transitions as well as thermodynamic parameters for all probes.The authors concluded that hAR undergo a three-state transition during unfolding by chemical denaturants. The present data could not support this statement. The emission maximum wavelength of the Trp fluorescence revealed no populations of any intermediates. An intermediate seems to populate at low denaturant concentrations when observed by Trp fluorescence intensity. However, intrinsic fluorescence intensity is easily affected by numerous environmental factors and is not as reliable as the emission maximum wavelength. Yes, there is an obvious increase in ANS fluorescence. Previously Turoverov and his co-authors have shown that the appearance of off-pathway aggregates, in some cases, can enhance ANS fluorescence.
[1] It is necessary to avoid the misleading interpretation by identifying the potential intermediate(s) via more methods.It is unclear whether all the solid lines are fitted data or not. Moreover, it is better to provide the residuals to facilitate the readers to evaluate whether the fitting is appropriate.In Figure 3, the ANS profiles seem to comprise more than one state and are likely to be the sum of at least two states centered at around 1 and 2 M GdnHCl?CD signals are recommended to report in molar ellipticities. The quality of the far-UV CD spectra is not enough. High concentration of denaturants may affect the signals at low wavelengths, generally below 210 nm. However, the low signal-to-noise seems to be affected by other factors. The authors can try to adjust the measuring parameters such as the slit width and number of repetitions. Furthermore, the authors mentioned that all chemicals were Sigma reagent grade. To my experience, GdnHCl with an ultrapure grade will produce much better signals that the lower grades. Note that the impurities or contaminated solvents will greatly affect spectroscopic measurements. Before spectroscopic measurements, the samples should be filtered and degassed. Detailed guidelines of spectroscopic measurements can be found in the paper by Kelly et al.
[2]The units of the fluorescence intensities in Figures 2 and 3 are quite different. Are these data normalized by the same method? Furthermore, the spectra of N in panels E-2 should refer to the state in the absence of denaturant. Why the spectrum of N in Figure 2E-2 is different from that in Figure 3E-2? To me, the spectrum of N in Figure 2E-2 is strange. Please verify that the solutions are not affected by contaminants.I confirm that I have read this submission and believe that I have an appropriate level of expertise to confirm that it is of an acceptable scientific standard, however I have significant reservations, as outlined above.In this manuscript, Gurprit Sekhon and Ranvir Singh describe the results of the equilibrium unfolding of humanaldose reductase. The authors conclude that urea- and GuHCL-induced unfolding of this protein is characterized by the presence of partially folded intermediate. This is an interesting study with some potential. However, additional experiments are needed to provide more evidence for the existence of an intermediate state and also to provide better description of the structural properties of this protein.The authors are encouraged to add near-UV CD spectroscopy to the arsenal of techniques utilized in their study. This will give very important information on the effect of denaturants on the tertiary structure of a protein. In fact, combined use of near- and far-UV spectroscopy is the accepted practice in studies on the conformational stability of proteins.Quality of the reported far-UV CD spectra is very low. Spectra are very noisy and their utilization for the analysis of the effects of urea and GuHCl on secondary structure of humanaldose reductase is questionable. The authors have to change settings to generate more reliable and less noisy spectra. Probably, protein concentrations should be increased too.Far-UV CD spectra of completely unfolded forms induced by high concentrations of urea and GuHCl are very different, suggesting that these unfolded forms are not similar. Why?Compactness of an intermediate state should be evaluated. This can be done by a whole host of hydrodynamic techniques.The authors should make sure that the protein does not aggregate at the conditions promoting formation of a partially unfolded intermediate. This can be done by simple light scattering experiments.Data shown in Figure 3 suggest that very low GuHCl concentrations cause very noticeable changes in some of the analyzed parameters. Similar behavior was described previously
[1] and was attributed to the aggregation of partially unfolded species. This phenomenon should be discussed.The presence of partially folded intermediates can be visualized using "phase diagram" method described in previous work by Kuznetsova et al.
[2] The authors are encouraged to use this approach for the analysis of their data.The authors are encouraged to reconsider the use of terms "denatured/denaturation" and "unfolded/unfolding". Typically, unfolding is attributed to the process resulting in the formation of coil-like conformation, whereas denaturation is referred to the elimination of functional tertiary structure. Obviously, these terms are not equivalent - molten globule is denatured, but is not unfolded. Furthermore, although temperature increase typically causes melting of a protein tertiary structure, temperature-denatured species are often rather compact and preserve high levels of secondary structure.The manuscript contains some linguistic issues and errors and definitely needs editing.I confirm that I have read this submission and believe that I have an appropriate level of expertise to confirm that it is of an acceptable scientific standard, however I have significant reservations, as outlined above.Drs Sekhon and Singh have presented an in-vitro study "Humanaldose reductase unfolds through an intermediate" to F1000Research to be considered for indexing. Foremost, the authors should describe the newly described intermediate more with respect to its potential physiological relevance to stress the relevance of this study.Please remove "ultra" from any description for protein crystal structures, but add the actual resolution of that structure, 0.66 A from Ref. 5, and explain what additional features were derived from these, compared to average crystal structures."Understanding the capability of a polypeptide..." - this sentence only has limited linkage to the sentence beforehand and should be deleted. Instead, novel concepts of protein folding and stability should be briefly covered.
[1]
,
[2]Figure 1:A) it is not clear what the blue and red dots mean.B) Please do not use green and red as colors.Figure 2:A+B) Please label the different graphs within the diagrams.The manuscript needs careful copy-editing.I confirm that I have read this submission and believe that I have an appropriate level of expertise to confirm that it is of an acceptable scientific standard, however I have significant reservations, as outlined above.