Thermal transition of core particle is not a two-state process

Thermal transition of core particle which occurs before melting of DNA and can be followed by circular dichroism is not a two-state process; it is the result of two processes which cannot be dissociated in static experiments: unfolding of core particles is immediately followed by their aggregation. It is thus impossible to get thermodynamic parameters of core particle unfolding from its thermal transition monitored by circular dichroism. Thermal denaturation kinetics of core particles gives some information about their stability. Finally core particle structure is more stable in chromatin than in its isolated state.     

Reassessment of the electronic circular dichroism criteria for random coil conformations of poly(L-lysine) and the implications for protein folding and denaturation studies

The circular dichroism (CD) spectra of poly(L-lysine) in water and ethanediol/water (2:1) solutions in the temperature range -110 to 85 degrees C are presented. The results combined with vibrational CD data are interpreted in terms of a two-state conformational equilibrium with a left-handed trans polyproline II conformation being preferred at low temperatures. The relevance of these studies to the CD criteria for random-coil conformations, the study of helix-coil transitions and protein/peptide folding is pointed out.     

Effect of Base Stacking on the Relative Thermodynamic Stability of Oligonucleotide Complexes: A Spectroscopic Study

Abstract

Three-strand oligonucleotide complexes are employed to assess the effect of base stacking and base pair mismatch on the relative thermodynamic stabilities of oligonucleotide duplexes. The melting behavior of three-strand oligonucleotide complexes incorporating nicks and gaps as well as internal single base mismatches is monitored using temperature-dependent optical absorption spectroscopy. A sequential three-state equilibrium model is used to analyze the measured melting profiles and evaluate thermodynamic parameters associated with dissociation of the complexes. The free-energy of stabilization of a nick complex compared to a gap complex due to base stacking is determined to be ?1.9 kcal/mol. The influence of a mispaired base in these systems is shown to destabilize a nick complex by 3.1 kcal/mol and a gap complex by 2.8 kcal/mol, respectively.     

Effect of base stacking on the relative thermodynamic stability of oligonucleotide complexes: a spectroscopic study

Three-strand oligonucleotide complexes are employed to assess the effect of base stacking and base pair mismatch on the relative thermodynamic stabilities of oligonucleotide duplexes. The melting behavior of three-strand oligonucleotide complexes incorporating nicks and gaps as well as internal single base mismatches is monitored using temperature-dependent optical absorption spectroscopy. A sequential three-state equilibrium model is used to analyze the measured melting profiles and evaluate thermodynamic parameters associated with dissociation of the complexes. The free-energy of stabilization of a nick complex compared to a gap complex due to base stacking is determined to be -1.9 kcal/mol. The influence of a mispaired base in these systems is shown to destabilize a nick complex by 3.1 kcal/mol and a gap complex by 2.8 kcal/mol, respectively.     

Guanidinium chloride and urea denaturations of beta-lactoglobulin A at pH 2.0 and 25 degrees C: the equilibrium intermediate contains non-native structures (helix, tryptophan and hydrophobic patches)

We have carried out guanidinium chloride (GdmCl) and urea denaturations of bovine beta-lactoglobulin A (beta-lgA) at pH 2.0 and 25 degrees C, using far-UV and near-UV circular dichroism, near-UV absorption and tryptophan fluorescence spectroscopies. The stable intermediate state that occurs during GdmCl denaturation has been characterized by the far- and near-UV circular dichroism, tryptophan difference absorption, tryptophan fluorescence and 8-anilino-1-naphthalene sulphonic acid binding measurements. Following conclusions have been reached. (a) Urea-induced denaturation is not a two-state process. (b) GdmCl-induced denaturation is composed of two distinct two-state processes. (c) alpha-Helical content, burial of tryptophan residues and burial of hydrophobic surface area are more in the GdmCl-induced stable intermediate than those originally present in the native protein.     

Conformational changes in chemically modified Escherichia coli thioredoxin monitored by H/D exchange and electrospray ionization mass spectrometry

Hydrogen/deuterium (H/D) exchange in combination with electrospray ionization mass spectrometry and near-ultraviolet (UV) circular dichroism (CD) was used to study the conformational properties and thermal unfolding of Escherichia coli thioredoxin and its Cys32-alkylated derivatives in 1% acetic acid (pH 2.7). Thermal unfolding of oxidized (Oxi) and reduced (Red) -thioredoxin (TRX) and Cys-32-ethylglutathionyl (GS-ethyl-TRX) and Cys-32-ethylcysteinyl (Cys-ethyl-TRX), which are derivatives of Red-TRX, follow apparent EX1 kinetics as charge-state envelopes, H/D mass spectral exchange profiles, and near-UV CD appear to support a two-state folding/unfolding model. Minor mass peaks in the H/D exchange profiles and nonsuperimposable MS- and CD-derived melting curves, however, suggest the participation of unfolding intermediates leading to the conclusion that the two-state model is an oversimplification of the process. The relative stabilities as measured by melting temperatures by both CD and mass spectral charge states are, Oxi-TRX, GS-ethyl-TRX, Cys-ethyl-TRX, and Red-TRX. The introduction of the Cys-32-ethylglutathionyl group provides extra stabilization that results from additional hydrogen bonding interactions between the ethylglutathionyl group and the protein. Near-UV CD data show that the local environment near the active site is perturbed to almost an identical degree regardless of whether alkylation at Cys-32 is by the ethylglutathionyl group, or the smaller, nonhydrogen-bonding ethylcysteinyl group. Mass spectral data, however, indicate a tighter structure for GS-ethyl-TRX.     

Salt-dependent and protein-concentration-dependent changes in the solution structure of the DNA-binding histone-like protein, HBsu, from Bacillus subtilis.

The solution structure of the histone-like DNA-binding protein, HBsu, from Bacillus subtilis in 2 mM sodium cacodylate, pH 7.5, is sensitive to the ionic strength of the buffer. This was shown by circular dichroism measurements at different concentrations of sodium chloride and potassium fluoride. The stability of HBsu is also influenced; at HBsu concentrations of about 0.1 mg.ml-1, melting temperatures of 32 degrees C and 55 degrees C were found in the absence of potassium fluoride and in the presence of 0.5 M potassium fluoride, respectively, exhibiting drastic ionic-strength-dependent differences in the temperature-induced unfolding of HBsu. Furthermore, at low ionic strength, circular dichroism spectra vary markedly depending on the HBsu concentration in the approximate range 0.2-3 mg.ml-1. Such protein-concentration-dependent differences in the spectra were not observed in the presence of 0.5 M potassium fluoride. Very similar circular dichroism spectra of HBsu and the histone-like DNA-binding protein of Bacillus stearothermophilus (HBst) at high ionic strength, indicate comparable structures of both proteins under these conditions. Estimation of the secondary structure content from the circular dichroism spectra yields data which are in satisfactory agreement with the values obtained from the crystal structure of HBst. Transition temperatures of 45 degrees C and 61 degrees C were found in differential scanning calorimetric measurements performed with HBsu in potassium-fluoride-free buffer and in the presence of 0.5 M potassium fluoride, respectively. The thermodynamic data point to the melting of native HBsu dimers into two denatured monomers.     

Equilibrium and kinetic studies of protein cooperativity using urea-induced folding/unfolding of a Ubq-UIM fusion protein

Understanding the origins of cooperativity in proteins remains an important topic in protein folding. This study describes experimental folding/unfolding equilibrium and kinetic studies of the engineered protein Ubq-UIM, consisting of ubiquitin (Ubq) fused to the sequence of the ubiquitin interacting motif (UIM) via a short linker. Urea-induced folding/unfolding profiles of Ubq-UIM were monitored by far-UV circular dichroism and fluorescence spectroscopies and compared to those of the isolated Ubq domain. It was found that the equilibrium data for Ubq-UIM is inconsistent with a two-state model. Analysis of the kinetics of folding shows similarity in the folding transition state ensemble between Ubq and Ubq-UIM, suggesting that formation of Ubq domain is independent of UIM. The major contribution to the stabilization of Ubq-UIM, relative to Ubq, was found to be in the rates of unfolding. Moreover, it was found that the kinetic m-values for Ubq-UIM unfolding, monitored by different probes (far-UV circular dichroism and fluorescence spectroscopies), are different; thereby, further supporting deviations from a two-state behavior. A thermodynamic linkage model that involves four states was found to be applicable to the urea-induced unfolding of Ubq-UIM, which is in agreement with the previous temperature-induced unfolding study. The applicability of the model was further supported by site-directed variants of Ubq-UIM that have altered stabilities of Ubq/UIM interface and/or stabilities of individual Ubq- and UIM-domains. All variants show increased cooperativity and one variant, E43N_Ubq-UIM, appears to behave very close to an equilibrium two-state.     

Thermodynamic characterization of cytochrome c at low pH. Observation of the molten globule state and of the cold denaturation process.

Several reports have pointed out the existence of intermediate states (both kinetic and equilibrium intermediate) between the native and the denatured states. The molten globule state, a compact intermediate state in which the secondary structure is formed but the tertiary structure fluctuates considerably, is currently being studied intensively because of its possible implication in the folding process of several proteins. We have examined the thermal stability of horse cytochrome c at low pH between 2.0 and 3.2 and different potassium chloride concentrations by absorbance of the Soret band, far and near-ultraviolet circular dichroism (u.v. c.d.) and tryptophan fluorescence using a multidimensional spectrophotometer. The concentration of potassium chloride ranged from 0 M to 0.5 M. The experimental thermal denaturation curves show that: (1) the helical content of cytochrome c remains stable at higher temperature when the concentration of salt is increased; whereas (2) the extent of ordering of the tertiary structure is weakly dependent on salt concentration; and (3) for cytochrome c, the stabilization of the molten globule state is induced by the binding of anions. Other salts such as NaCl, LiCl, potassium ferricyanide (K3Fe(CN)6) and Na2SO4 may also be used to stabilize the molten globule state. The thermodynamic analysis of the denaturation curves of c.d. at 222 nm and c.d. at 282 nm shows that, whereas a two-state (native and denatured) transition is observed at low-salt concentration, the far and near-u.v. c.d. melting curves of cytochrome c do not coincide with each other at high-salt concentration, and a minimum of three different thermodynamic states (IIb, intermediate or IIc, and denatured) is necessary to achieve a sufficient analysis. The intermediate state (called IIc) is attributed to the molten globule state because of its high secondary structure content and the absence of tertiary structure. Therefore, at low pH, cytochrome c is present in at least four states (native, IIb, IIc and denatured) depending on the salt concentration and temperature. The thermodynamic parameters, i.e. the Gibbs free energy differences (delta G), the enthalpy differences (delta H), the midpoint temperatures (Tm) of the transition (IIb in equilibrium intermediate (IIc in equilibrium denatured) are determined. We also give estimates of the heat capacity differences (delta Cp) from the temperature dependence of the enthalpy differences. The enthalpy change and the heat capacity difference of the IIc in equilibrium denatured transition are non-zero. The number of charges (protons or chloride anions) released upon transitions are determined by analysing the pH and chloride anion concentration dependence of the Gibbs free energy.(ABSTRACT TRUNCATED AT 400 WORDS)     

Detection and characterization of partially unfolded oligomers of the SH3 domain of alpha-spectrin

For the purpose of equilibrium and kinetic folding-unfolding studies, the SH3 domain of alpha-spectrin (spc-SH3) has long been considered a classic two-state folding protein. In this work we have indeed observed that the thermal unfolding curves of spc-SH3 measured at pH 3.0 by differential scanning calorimetry, circular dichroism, and NMR follow apparently the two-state model when each unfolding profile is considered individually. Nevertheless, we have found that protein concentration has a marked effect upon the thermal unfolding profiles. This effect cannot be properly explained in terms of the two-state unfolding model and can only be interpreted in terms of the accumulation of intermediate associated states in equilibrium with the monomeric native and unfolded states. By chemical cross-linking and pulsed-field gradient NMR diffusion experiments we have been able to confirm the existence of associated states formed during spc-SH3 unfolding. A three-state model, in which a dimeric intermediate state is assumed to be significantly populated, provides the simplest interpretation of the whole set of thermal unfolding data and affords a satisfactory explanation for the concentration effects observed. Whereas at low concentrations the population of the associated intermediate state is negligible and the unfolding process consequently takes place in a two-state fashion, at concentrations above approximately 0.5 mM the population of the intermediate state becomes significant at temperatures between 45 degrees C and 80 degrees C and reaches up to 50% at the largest concentration investigated. The thermodynamic properties of the intermediate state implied by this analysis fall in between those of the unfolded state and the native ones, indicating a considerably disordered conformation, which appears to be stabilized by oligomerization.     

Unfolding of human lens recombinant betaB2- and gammaC-crystallins

betaB2- and gammaC-crystallins belong to the betagamma-crystallin superfamily and have very similar structures. Molecular spectroscopic techniques such as UV-visible absorption, circular dichroism, and fluorescence indicate they have similar biophysical properties. Their structures are characterized by the presence of two domains consisting of four Greek key motifs. The only difference is the connecting peptide of the two domains, which is flexible in gamma-crystallin but extended in beta-crystallin; thus, an intradomain association and a monomer are formed in gamma-crystallin and an interdomain association and a dimer are formed in beta-crystallin. The difference may be reflected in the thermodynamic stability. In the present study, we calculated the standard free-energy by equilibrium unfolding transition in guanidine hydrochloride using three spectroscopic parameters: absorbance at 235nm, Trp fluorescence intensity at 320nm, and far-UV circular dichroism at 223nm. Global analyses indicate that both dimeric betaB2- and monomeric gammaC-crystallins are a better fit to a three-state model than to a two-state model. In terms of standard free-energy, deltaG(0)(H(2)O,i) both betaB2-crystallin and gammaC-crystallin are stable proteins and dimeric betaB2-crystallin is more stable than the monomeric gammaC-crystallin. The significance of the thermodynamic stability for betaB2- and gammaC-crystallins may be related to their functions in the lens.     

Melting studies of short DNA hairpins containing the universal base 5-nitroindole.

Effects of the universal base 5-nitroindole on the thermodynamic stability of DNA hairpins having a 6 bp stem and four base loops were investigated by optical absorbance and differential scanning calorimetry techniques. Melting studies were conducted in buffer containing 115 mM Na(+). Five different modified versions of DNA hairpins containing a 5-nitroindole base or bases substituted at different positions in the stem and loop regions were examined. Thermo-dynamic parameters of the melting transitions estimated from a two-state analysis of optical melting curves and measured directly by calorimetry revealed that the presence of 5-nitroindole bases in the duplex stem or loop regions of short DNA hairpins significantly affects both their enthalpic and entropic melting components in a compensating manner, while the transition free energy varies linearly with the transition temperature. The calorimetrically determined enthalpy and entropy values of the modified hairpins were considerably smaller (43-53%) than the two-state optical parameters, suggesting that solvent effects may be significant in the melting processes of these hairpins. Results of circular dichroism measurements also revealed slight differences between the modified hairpins and the control in both the duplex and melted states, suggesting subtle structural differences between the control and DNA hairpins containing a 5-nitroindole base or bases.     

Studies of the optical properties of the proflavine–DNA complex

We report studies of the optical properties of the proflavine–DNA complex, using absorbance and circular dichroism spectroscopy. From comparison of the absorption spectra of proflavine complexed with calf thymus and T2 DNA, we conclude that stacking of the dyes external to the double helix is comparatively much weaker with T2 DXA, probably because of its glucosylation. Several sources are found for the circular dichroism induced in proflavine when it is complexed with DNA. There is a relatively weak circular dichroism induced when the dye is infinitely dilute on the DNA lattice; this presumably arises from the environmental asymmetry of the binding site. Stronger circular dichroism effects are induced by interaction of intercalated and stacked dyes; studies with T2 DNA, for which stacking seems to be blocked, permit a tentative resolution of effects due to the two modes of binding. One recurring theme of these studies is the observation that the optical properties are quite dependent on environment. The most dramatic example is a strong variation with salt concentration of the amplitude of the circular dichroism induced in the isolated (intercalated) monomer by the surrounding DNA. This suggests that the structure of the intercalated complex is quite sensitive to external conditions.     

Urea-induced denaturation of human serum albumin labeled with acrylodan

We induced the denaturation of unlabeled human serum albumin (HSA) and of similar albumin labeled with acrylodan (6-acryloyl-2-dimethylamino naphthalene) with urea and studied the transition profiles using circular dichroism and fluorescence spectroscopy. The circular dichroism spectra for both albumin preparations resulted in the same curves, thus indicating that labeling with acrylodan does not perturb the conformation of HSA. Our results indicate that the denaturation of both albumin preparations takes place at a single, two-state transition with midpoint at about 6 M urea, due to the unfolding of its domain II. It is important to point out that even at 8 M urea, some residual structure remains in the HSA. Great changes in the fluorescence of the dye bound to the protein were observed by addition of solid guanidine hydrochloride to the protein labeled with acrylodan dissolved in 8 M urea, indicating that domain I of this protein was not denatured by urea.     

Resolution of a structural competition involving dimeric G-quadruplex and its C-rich complementary strand

The resolution of the dimeric intermolecular G-quadruplex/duplex competition of the telomeric DNA sequence 5′-TAG GGT TAG GGT-3′ and of its complementary 5′ ACC CTA ACC CTA-3′ is reported. To achieve this goal, melting experiments of both sequences and of the mixtures of these sequences were monitored by molecular absorption, molecular fluorescence and circular dichroism spectroscopies. Molecular fluorescence measurements were carried out using molecular beacons technology, in which the 5′-TAG GGT TAG GGT-3′ sequence was labelled with a fluorophore and a quencher at the ends of the strand. Mathematical analysis of experimental spectroscopic data was performed by means of multivariate curve resolution, allowing the calculation of concentration profiles and pure spectra of all resolved structures (dimeric antiparallel and parallel G-quadruplexes, Watson–Crick duplex and single strands) present in solution. Our results show that parallel G-quadruplex is more stable than antiparallel G-quadruplex. When the complementary C-rich strand is present, a mixture of both G-quadruplex structures and Watson–Crick duplex is observed, the duplex being the major species. In addition to melting temperatures, equilibrium constants for the parallel/antiparallel G-quadruplex equilibrium and for the G-quadruplex/duplex equilibrium were determined from the concentration profiles.     

Evidence for identity between the equilibrium unfolding intermediate and a transient folding intermediate: a comparative study of the folding reactions of alpha-lactalbumin and lysozyme

The refolding kinetics of alpha-lactalbumin at different concentrations of guanidine hydrochloride have been investigated by means of kinetic circular dichroism and stopped-flow absorption measurements. The refolding reaction consists of at least two stages, the instantaneous accumulation of the transient intermediate that has peptide secondary structure and the subsequent slow process associated with formation of tertiary structure. The transient intermediate is compared with the well-characterized equilibrium intermediate observed during the denaturant-induced unfolding. Stabilities of the secondary structures against the denaturant, affinities for Ca2+, and tryptophan absorption properties of the transient and equilibrium intermediates were investigated. In all of these respects, the transient intermediate is identical with the equilibrium one, demonstrating the validity of the use of the equilibrium intermediate as a model of the folding intermediate. Essentially the same transient intermediate was also detected in the folding of lysozyme, the protein known to be homologous to alpha-lactalbumin but whose equilibrium unfolding is represented as a two-state reaction. The stability and cooperativity of the secondary structure of the intermediate of lysozyme are compared with those of alpha-lactalbumin. The results show that the protein folding occurring via the intermediate is not limited to the proteins that show equilibrium intermediates. Although the unfolding equilibria of most proteins are well approximated as a two-state reaction, the two-state hypothesis may not be applicable to the folding reaction under the native condition. Two models of protein folding, intermediate-controlled folding model and multiple-pathway folding model, which are different in view of the role of the intermediate in determining the pathway of folding, are also discussed.     

Conformational stability of the covalent complex between elastase and alpha 1-proteinase inhibitor

The equilibrium unfolding-refolding process of the elastase-alpha 1-proteinase inhibitor complex, induced by guanidinium chloride, was followed by spectroscopic methods. A reversible transition with a midpoint at 2.04 +/- 0.04 M guanidinium chloride was observed by fluorescence. This transition was attributed to elastase on the basis of circular dichroism and uv absorption difference data obtained for the covalent complex and for the free proteins. The conformational stability of elastase in the complex was analyzed considering the approximation of a two-state transition. The free energy of denaturation delta GH2O was 4.2 kcal.mol-1 for complexed elastase compared to 10.5 kcal.mol-1 for the free enzyme. Such a decrease in the stability of elastase suggests that, after forming the covalent complex with the inhibitor, the enzyme undergoes not only the expected local modifications of the active site, but also an extensive structural reorganization.     

Ca2+-induced alteration in the unfolding behavior of alpha-lactalbumin

Comparative studies of the unfolding equilibria of two homologous proteins, bovine alpha-lactalbumin and hen lysozyme, induced by treatment with guanidine hydrochloride have been made by analysis of the peptide and the aromatic circular dichroism spectra. The effect of the specific binding of Ca2+ ion by the former protein was taken into account in interpreting the unfolding equilibria of the protein. Proton nuclear magnetic resonance spectra of alpha-lactalbumin were also measured for the purpose of characterizing an intermediate structural state of the protein. In previous studies, alpha-lactalbumin was shown to be an exceptional protein whose equilibrium unfolding does not obey the two-state model of unfolding, although lysozyme is known to follow the two-state unfolding mechanism. The present results show that the apparent unfolding behavior of alpha-lactalbumin depends on Ca2+ concentration. At a low concentration of Ca2+, alpha-lactalbumin unfolds with a stable intermediate that has unfolded tertiary structure, as evidenced by the featureless nuclear magnetic resonance and aromatic circular dichroism spectra, but has folded secondary structure as evidenced by the peptide circular dichroism spectra. However, in the presence of a sufficiently high concentration of Ca2+, the unfolding transition of alpha-lactalbumin resembles that of lysozyme. The transition occurs between the two states, the native and the fully unfolded states, and the cooperativity of the unfolding is essentially the same as that of lysozyme. Such a change in the apparent unfolding behavior evidently results from an increase in the stability of the native state relative to that of the intermediate induced by the specific Ca2+ binding to native alpha-lactalbumin. The results are useful for understanding the relationship between the protein stability and the apparent unfolding behavior.     

Thermodynamic, spectroscopic, and equilibrium binding studies of DNA sequence context effects in six 22-base pair deoxyoligonucleotides.

Effects of different end sequences on stability, circular dichroism spectra (CD), and enzyme binding properties were investigated for six 22-base pair, non-self-complementary duplex DNA oligomers. The center sequences of these deoxyoligonucleotides have 8-14 base pairs in common and are flanked on both sides by sequences differing in context and A-T content. Temperature-induced melting transitions monitored by differential scanning calorimetry (DSC) and ultraviolet absorbance were measured for the six duplexes in buffered 115 mM Na(+) solutions. Values of the melting transition enthalpy, DeltaH(cal), and entropy, DeltaS(cal), were obtained directly from DSC experiments. Melting transition parameters, DeltaH(vH) and DeltaS(vH), were also estimated from van't Hoff analysis of optical melting curves collected as a function of DNA concentration, assuming a two-state melting transition. Melting free energies (20 degrees C) of the six DNAs evaluated from DSC experiments ranged from -18.7 to -32.7 kcal/mol. van't Hoff estimates of the free energies ranged from -18.5 to -48.0 kcal/mol. With either method, the trends in free energy as a function of sequence were identical. Equilibrium binding by BamHI restriction endonuclease to the 22-base pair DNAs was also investigated. The central eight base pairs of all six molecules, 5'-A-GGATCC-A-3', contained a BamHI recognition sequence bounded by A-T base pairs. Magnesium free binding assays were performed by titering BamHI against a constant concentration of each of the deoxyoligonucleotide substrates and analyzing reaction products by gel retardation. Binding isotherms of the total amount of bound DNA versus protein concentration were constructed which provided semiquantitative estimates of the equilibrium dissociation constants for dissociation of BamHI from the six DNA oligomers. Dissociation constants ranged from 0.5 x 10(-)(9) to 12.0 x 10(-)(9) M with corresponding binding free energies of -12.5 to -10.6 (+/-0. 1) kcal/mol. An inverse relationship is found when binding and stability are compared.     

Thermostability of the N-Terminal RNA-Binding Domain of the SARS-CoV Nucleocapsid Protein: Experiments and Numerical Simulations

Differential scanning calorimetry, circular dichroism spectroscopy, nuclear magnetic resonance spectroscopy, and numerical simulations were used to study the thermostability of the N-terminal RNA-binding domain (RBD) of the SARS-CoV nucleocapsid protein. The transition temperature of the RBD in a mixing buffer, composed of glycine, sodium acetate, and sodium phosphate with 100 mM sodium chloride, at pH 6.8, determined by differential scanning calorimetry and circular dichroism, is 48.74°C. Experimental results showed that the thermal-induced unfolding-folding transition of the RBD follows a two-state model with a reversibility >90%. Using a simple Gō-like model and Langevin dynamics we have shown that, in agreement with our experiments, the folding of the RBD is two-state. Theoretical estimates of thermodynamic quantities are in reasonable agreement with the experiments. Folding and thermal unfolding pathways of the RBD also were experimentally and numerically studied in detail. It was shown that the strand β₁ from the N-terminal folds last and unfolds first, while the remaining β-strands fold/unfold cooperatively.