Equilibrium and kinetic studies on folding of the authentic and recombinant forms of human alpha-lactalbumin by circular dichroism spectroscopy

The equilibrium and kinetics of the unfolding and refolding of authentic and recombinant human alpha-lactalbumin, the latter of which had an extra methionine residue at the N-terminus, were studied by circular dichroism spectroscopy, and the results were compared with the results for bovine and goat alpha-lactalbumins obtained in our previous studies. As observed in the bovine and goat proteins, the presence of the extra methionine residue in the recombinant protein remarkably destabilized the native state, and the destabilization was entirely ascribed to an increase in the rate of unfolding. The thermodynamic stability of the native state against the unfolded state was lower, and the thermodynamic stability of the molten globule state against the unfolded state was higher for the human protein than for the other alpha-lactalbumins previously studied. Thus, the population of the molten globule intermediate was higher during the equilibrium unfolding of human alpha-lactalbumin by guanidine hydrochloride. Unlike the molten globule states of the bovine and goat proteins, the human alpha-lactalbumin molten globule showed remarkably more intense circular dichroism ellipticity than the native state in the far-ultraviolet region below 225 nm. During refolding from the unfolded state, human alpha-lactalbumin thus exhibited overshoot kinetics, in which the alpha-helical peptide ellipticity exceeded the native value when the molten globule folding intermediate was formed in the burst phase. The subsequent folding involved reorganization of nonnative secondary structures. It should be noted that the rate constant of the major refolding phase was approximately the same among the three types of alpha-lactalbumin and that the rate constant of unfolding was accelerated 18-600 times in the human protein, and these results interpreted the lower thermodynamic stability of this protein.     

Unfolding pathways of goat alpha-lactalbumin as revealed in multiple alignment of molecular dynamics trajectories

Molecular dynamics simulations of protein unfolding were performed at an elevated temperature for the authentic and recombinant forms of goat alpha-lactalbumin. Despite very similar three-dimensional structures, the two forms have significantly different unfolding rates due to an extra N-terminal methionine in the recombinant protein. To identify subtle differences between the two forms in the highly stochastic kinetics of unfolding, we classified the unfolding trajectories using the multiple alignment method based on the analogy between the biological sequences and the molecular dynamics trajectories. A dendrogram derived from the multiple trajectory alignment revealed a clear difference in the unfolding pathways of the authentic and recombinant proteins, i.e. the former reached the transition state in an all-or-none manner while the latter unfolded less cooperatively. It was also found in the classification that the two forms of the protein shared a common transition state structure, which was in excellent agreement with the transition state structure observed experimentally in the Phi-value analysis.     

Retardation of the unfolding process by single N-glycosylation of ribonuclease A based on molecular dynamics simulations

The conformational characteristics of glycosylated- and unglycosylated bovine pancreatic ribonuclease A (RNaseA) were traced with unfolding molecular dynamics simulations using CHARMM program at 470 K. The glycosylated RNase (Glc_RNase) possesses nearly identical protein structure with RNaseA, differing only by presence of a single acetylglucosamine residue N-linked to Asn34 in the RNaseA. Attaching of monomeric N-acetylglucosamine residue to the Asn34 in RNaseA resulted in a change of denaturing process of Glc_RNase. Simulations showed that the unfolding of RNaseA involved significant weakening of nonlocal interactions whereas the glycosylation led Glc_RNase to preserve the nonlocal interactions even in its denatured form. Even in simulations over 8 ns at 470 K, Glc_RNase remained relatively stable as a less denatured conformation. However, conformation of RNaseA was changed to a fully unfolded state before 3 ns of the simulations at 470 K. This difference was due to fact that formation of hydrogen bond bridges and nonlocal contacts induced by the attached N-acetylglucosamine of Glc_RNase showing in the unfolding simulations. These high-temperature unfolding MD simulations provided a theoretical basis for the previous experimental work in which Glc_RNase showed slower unfolding kinetics compared with unglycosylated RNaseA, suggesting that single N-glycosylation induced retardation of unfolding process of the ribonuclease protein.     

Increasing temperature accelerates protein unfolding without changing the pathway of unfolding

We have traditionally relied on extremely elevated temperatures (498K, 225 degrees C) to investigate the unfolding process of proteins within the timescale available to molecular dynamics simulations with explicit solvent. However, recent advances in computer hardware have allowed us to extend our thermal denaturation studies to much lower temperatures. Here we describe the results of simulations of chymotrypsin inhibitor 2 at seven temperatures, ranging from 298K to 498K. The simulation lengths vary from 94ns to 20ns, for a total simulation time of 344ns, or 0.34 micros. At 298K, the protein is very stable over the full 50ns simulation. At 348K, corresponding to the experimentally observed melting temperature of CI2, the protein unfolds over the first 25ns, explores partially unfolded conformations for 20ns, and then refolds over the last 35ns. Above its melting temperature, complete thermal denaturation occurs in an activated process. Early unfolding is characterized by sliding or breathing motions in the protein core, leading to an unfolding transition state with a weakened core and some loss of secondary structure. After the unfolding transition, the core contacts are rapidly lost as the protein passes on to the fully denatured ensemble. While the overall character and order of events in the unfolding process are well conserved across temperatures, there are substantial differences in the timescales over which these events take place. We conclude that 498K simulations are suitable for elucidating the details of protein unfolding at a minimum of computational expense.     

Unfolding a linker between helical repeats

In many multi-repeat proteins, linkers between repeats have little secondary structure and place few constraints on folding or unfolding. However, the large family of spectrin-like proteins, including alpha-actinin, spectrin, and dystrophin, share three-helix bundle, spectrin repeats that appear in crystal structures to be linked by long helices. All of these proteins are regularly subjected to mechanical stress. Recent single molecule atomic force microscopy (AFM) experiments demonstrate not only forced unfolding but also simultaneous unfolding of tandem repeats at finite frequency, which suggests that the contiguous helix between spectrin repeats can propagate a cooperative helix-to-coil transition. Here, we address what happens atomistically to the linker under stress by steered molecular dynamics simulations of tandem spectrin repeats in explicit water. The results for alpha-actinin repeats reveal rate-dependent pathways, with one pathway showing that the linker between repeats unfolds, which may explain the single-repeat unfolding pathway observed in AFM experiments. A second pathway preserves the structural integrity of the linker, which explains the tandem-repeat unfolding event. Unfolding of the linker begins with a splay distortion of proximal loops away from hydrophobic contacts with the linker. This is followed by linker destabilization and unwinding with increased hydration of the backbone. The end result is an unfolded helix that mechanically decouples tandem repeats. Molecularly detailed insights obtained here aid in understanding the mechanical coupling of domain stability in spectrin family proteins.     

Effects of helix and fingertip mutations on the thermostability of xyn11A investigated by molecular dynamics simulations and enzyme activity assays

Local conformational changes and global unfolding pathways of wildtype xyn11A recombinant and its mutated structures were studied through a series of atomistic molecular dynamics (MD) simulations, along with enzyme activity assays at three incubation temperatures to investigate the effects of mutations at three different sites to the thermostability. The first mutation was to replace an unstable negatively charged residue at a surface beta turn near the active site (D32G) by a hydrophobic residue. The second mutation was to create a disulphide bond (S100C/N147C) establishing a strong connection between an alpha helix and a distal beta hairpin associated with the thermally sensitive Thumb loop, and the third mutation add an extra hydrogen bond (A155S) to the same alpha helix. From the MD simulations performed, MM/PBSA energy calculations of the unfolding energy were in a good agreement with the enzyme activities measured from the experiment, as all mutated structures demonstrated the improved thermostability, especially the S100C/N147C proved to be the most stable mutant both by the simulations and the experiment. Local conformational analysis at the catalytic sites and the xylan access region also suggested that mutated xyn11A structures could accommodate xylan binding. However, the analysis of global unfolding pathways showed that structural disruptions at the beta sheet regions near the N-terminal were still imminent. These findings could provide the insight on the molecular mechanisms underlying the enhanced thermostability due to mutagenesis and changes in the protein unfolding pathways for further protein engineering of the GH11 family xylanase enzymes.     

Structural and dynamic effects of alpha-helix deletion in Sso7d: implications for protein thermal stability

Sso7d is a 62-residue protein from the hyperthemophilic archaeon Sulfolobus solfataricus with a denaturation temperature close to 100 degrees C around neutral pH. An engineered form of Sso7d truncated at leucine 54 (L54Delta) is significantly less stable, with a denaturation temperature of 53 degrees C. Molecular dynamics (MD) studies of Sso7d and its truncated form at two different temperatures have been performed. The results of the MD simulations at 300 K indicate that: (1) the flexibility of Sso7d chain at 300 K agrees with that detected from X-ray and NMR structural studies; (2) L54Delta remains stable in the native folded conformation and possesses an overall dynamic behavior similar to that of the parent protein. MD simulations performed at 500 K, 10 ns long, indicate that, while Sso7d is in-silico resistant to high temperature, the truncated variant partially unfolds, revealing the early phases of the thermal unfolding pathway of the protein. Analysis of the trajectories of L54Delta suggests that the unzipping of the N-terminal and C-terminal beta-strands should be the first event of the unfolding pathway, and points out the regions more resistant to thermal unfolding. These findings allow one to understand the role played by specific interactions connecting the two ends of the chain for the high thermal stability of Sso7d, and support recent hypotheses on its folding mechanism emerged from site-directed mutagenesis studies.     

Unfolding of the cold shock protein studied with biased molecular dynamics

The cold shock protein from Bacillus caldolyticus is a small beta-barrel protein that folds in a two-state mechanism. For the native protein and for several mutants, a wealth of experimental data are available on stability and folding, so that it is an optimal system to study this process. We compare data from unfolding simulations (trajectories of 5 and up to 12 ns) obtained with a bias potential at room temperature and from unbiased thermal unfolding simulations with experimental data. The unfolding patterns derived from the trajectories starting from different native-like conformations and subject to different unfolding conditions agree. The transition state found in the simulations of unfolding is close to the native structure in agreement with experiment. Moreover, a lower value of the free energy barrier of unfolding was found for the mutant R3E than for the mutant E46A and the native protein, as indicated by experimental data. The first unfolding event involves the three-stranded beta-sheet whose decomposition corresponds to the transition state. In contrast to conclusions drawn from experiments, we found that the two-stranded beta-strand forms the most stable substructure, which decomposes very late in the unfolding process. However, assuming that this structure forms very early in the folding process, our findings would not contradict the experiments but require a different interpretation of them.     

Diffusing and colliding: the atomic level folding/unfolding pathway of a small helical protein

Proteins with ultra-fast folding/unfolding kinetics are excellent candidates for study by molecular dynamics. Here, we describe such simulations of a three helix bundle protein, the engrailed homeodomain (En-HD), which folds via the diffusion-collision model. The unfolding pathway of En-HD was characterized by seven simulations of the protein and 12 simulations of its helical fragments yielding over 1.1 micros of simulation time in water. Various conformational states along the unfolding pathway were identified. There is the compact native-like transition state, a U-shaped helical intermediate and an unfolded state with dynamic helical segments. Each of these states is in good agreement with experimental data. Examining these states as well as the transitions between them, we find the role of long-range tertiary contacts, specifically salt-bridges, important in the folding/unfolding pathway. In the folding direction, charged residues form long-range tertiary contacts before the hydrophobic core is formed. The formation of HII is assisted by a specific salt-bridge and by non-specific (fluctuating) tertiary contacts, which we call contact-assisted helix formation. Salt-bridges persist as the protein approaches the transition state, stabilizing HII until the hydrophobic core is formed. To complement this information, simulations of fragments of En-HD illustrate the helical propensities of the individual segments. By thermal denaturation, HII proved to be the least stable helix, unfolding in less than 450 ps at high temperature. We observed the low helical propensity of C-terminal residues from HIII in fragment simulations which, when compared to En-HD unfolding simulations, link the unraveling of HIII to the initial event that drives the unfolding of En-HD.     

Sequential four-state folding/unfolding of goat α-lactalbumin and its N-terminal variants

Equilibria and kinetics of folding/unfolding of α-lactalbumin and its two N-terminal variants were studied by circular dichroism spectroscopy. The two variants were wild-type recombinant and Glu1-deletion (E1M) variants expressed in Escherichia coli. The presence of an extra methionine at the N terminus in recombinant α-lactalbumin destabilized the protein by 2 kcal/mol, while the stability was recovered in the E1M variant in which Glu1 was replaced by Met1. Kinetic folding/unfolding reactions of the proteins, induced by stopped-flow concentration jumps of guanidine hydrochloride, indicated the presence of a burst-phase in refolding, and gave chevron plots with significant curvatures in both the folding and unfolding limbs. The folding-limb curvature was interpreted in terms of accumulation of the burst-phase intermediate. However, there was no burst phase observed in the unfolding kinetics to interpret the unfolding-limb curvature. We thus assumed a sequential four-state mechanism, in which the folding from the burst-phase intermediate takes place via two transition states separated by a high-energy intermediate. We estimated changes in the free energies of the burst-phase intermediate and two transition states, caused by the N-terminal variations and also by the presence of stabilizing calcium ions. The Φ values at the N terminus and at the Ca(2+)-binding site thus obtained increased successively during folding, demonstrating the validity of the sequential mechanism. The stability and the folding behavior of the E1M variant were essentially identical to those of the authentic protein, allowing us to use this variant as a pseudo-wild-type α-lactalbumin in future studies.     

Temperature-dependent downhill unfolding of ubiquitin. II. Modeling the free energy surface

To provide evidence for the interpretation of temperature‐dependent unfolding kinetics and the downhill unfolding scenario presented in the accompanying experimental article (Part I), the free energy surface of ubiquitin unfolding is calculated using statistical mechanical models of the Muñoz‐Eaton (ME) form. The models allow only two states for each amino acid residue, folded or unfolded, and permutations of these states generate an ensemble of microstates. One‐dimensional free energy curves are calculated using the number of folded residues as a reaction coordinate. The proposed sequential unfolding of ubiquitin's β‐sheet is tested by mapping the free energy onto two reaction coordinates inspired by the experiment as follows: the number of folded residues in ubiquitin's stable β‐strands I and II and those of the less stable strands III–V. Although the original ME model successfully captures folding features of zipper‐like one‐dimensional folders, it misses important tertiary interactions between residues that are far from each other in primary sequence. To take tertiary contacts into account, partially folded microstates based on a spherical growth model are included in the calculation and compared with the original model. By calculating the folding probability of each residue for a given point on the free energy surface, the unfolding pathway of ubiquitin is visualized. At low temperature, thermal unfolding occurs along a sequential unfolding pathway as follows: disruption of the β‐strands III–V followed by unfolding of the strands I and II. At high temperature, multiple unfolding routes are formed. The heterogeneity of the transition state explains the global nonexponential unfolding observed in the T‐jump experiment at high temperature. The calculation also reports a high stability for the α‐helix of ubiquitin. Proteins 2008. © 2008 Wiley‐Liss, Inc.     

A study of the structural properties and thermal stability of human Bcl-2 by molecular dynamics simulations

The anti-apoptotic B-cell lymphoma 2 (Bcl-2) protein interacts with several proteins that regulate the apoptotic properties of cells. In this research, we conduct several all-atom molecular dynamics (MD) simulations under high-temperature unfolding conditions, from 400 to 800?K, for 25?ns. These simulations were performed using a model of an engineered Bcl-2 human protein (Bcl-2-Δ22Σ3), which lacks 22 C-terminal residues of the transmembrane domain. The aim of this study is to gain insight into the structural behavior of Bcl-2-Δ22Σ3 by mapping the conformational movements involved in Bcl-2 stability and its biological function. To build a Bcl-2-Δ22Σ3 three-dimensional model, the protein core was built by homology modeling and the flexible loop domain (FLD, residues 33-91) by ab initio methods. Further, the entire protein model was refined by MD simulations. Afterwards, the production MD simulations showed that the FLD at 400 and 500?K has several conformations reaching into the protein core, whereas at 600?K some of the alpha-helices were lost. At 800?K, the Bcl-2 core is destabilized suggesting a possible mechanism for protein unfolding, where the alpha helices 1 and 6 were the most stable, and a reduction in the number of hydrogen bonds initially occurs. In conclusion, the structural changes and the internal protein interactions suggest that the core and the FLD are crucial components of Bcl-2 in its function of regulate ng access to the recognition sites of kinases and caspases.     

Temperature-dependent downhill unfolding of ubiquitin. I. Nanosecond-to-millisecond resolved nonlinear infrared spectroscopy

Transient thermal unfolding of ubiquitin is investigated using nonlinear infrared spectroscopy after a nanosecond laser temperature jump (T-jump). The abrupt change in the unfolding free energy surface and the ns time resolution allow us to observe a fast response on ns to micros time-scales, which we attribute to downhill unfolding, before a cross-over to ms kinetics. The downhill unfolding by a sub-population of folded proteins is induced through a shift of the barrier toward the native state. By adjusting the T-jump width, the effect of the initial (T(i)) and final (T(f)) temperature on the unfolding dynamics can be separated. From the amplitude of the fast downhill unfolding, the fractional population prepared at the unfolding transition state is obtained. This population increases with both T(i) and with T(f). A two-state kinetic analysis of the ms refolding provides thermodynamic information about the barrier height. By a combination of the fast and slow unfolding and folding parameters, a quasi-two-state kinetic analysis is performed to calculate the time-dependent population changes of the folded state. This calculation coincides with the experimentally obtained population changes at low temperature but deviations are found in the T-jump from 67 to 78 degrees C. Using temperature-dependent barrier height changes, a temperature Phi value analysis is performed. The result shows a decreasing trend of Phi(T) with temperature, which indicates an increase of the heterogeneity of the transition state. We conclude that ubiquitin unfolds along a well-defined pathway at low temperature which expands with increasing temperature to include multiple routes.     

Potentially amyloidogenic conformational intermediates populate the unfolding landscape of transthyretin: Insights from molecular dynamics simulations

Protein aggregation into insoluble fibrillar structures known as amyloid characterizes several neurodegenerative diseases, including Alzheimer's, Huntington's and Creutzfeldt‐Jakob. Transthyretin (TTR), a homotetrameric plasma protein, is known to be the causative agent of amyloid pathologies such as FAP (familial amyloid polyneuropathy), FAC (familial amyloid cardiomiopathy) and SSA (senile systemic amyloidosis). It is generally accepted that TTR tetramer dissociation and monomer partial unfolding precedes amyloid fibril formation. To explore the TTR unfolding landscape and to identify potential intermediate conformations with high tendency for amyloid formation, we have performed molecular dynamics unfolding simulations of WT‐TTR and L55P‐TTR, a highly amyloidogenic TTR variant. Our simulations in explicit water allow the identification of events that clearly discriminate the unfolding behavior of WT and L55P‐TTR. Analysis of the simulation trajectories show that (i) the L55P monomers unfold earlier and to a larger extent than the WT; (ii) the single α‐helix in the TTR monomer completely unfolds in most of the L55P simulations while remain folded in WT simulations; (iii) L55P forms, early in the simulations, aggregation‐prone conformations characterized by full displacement of strands C and D from the main β‐sandwich core of the monomer; (iv) L55P shows, late in the simulations, severe loss of the H‐bond network and consequent destabilization of the CBEF β‐sheet of the β‐sandwich; (v) WT forms aggregation‐compatible conformations only late in the simulations and upon extensive unfolding of the monomer. These results clearly show that, in comparison with WT, L55P‐TTR does present a much higher probability of forming transient conformations compatible with aggregation and amyloid formation.     

Molecular dynamics simulation and experimental validation of the effect of pH on protein desorption in hydrophobic charge induction chromatography

The ligands in hydrophobic charge induction chromatography (HCIC) are hydrophobic and ionisable. Thus, the pH is crucial for the separation performance in HCIC, especially for elution. However, it is difficult to obtain the microscopic information in HCIC through experimental means. In this work, molecular dynamics simulations are performed to examine the effect of pH on elution and protein conformational transition in HCIC, using a 46-bead β-barrel coarse-grained model protein and an HCIC adsorbent pore model constructed in an earlier work. Corresponding experiments are carried out for the validation of simulation results, using lysozyme and MEP Hypercel. Both the activities and fluorescence of lysozyme are examined to evaluate the conformational transition. The simulations indicate that the elution efficiency of protein increases with decreasing pH value in a non-linear manner. This is qualitatively consistent with the experimental results. MD simulations indicate that protein unfolding occurs in elution at all pH values. However, the experimental data show that the activity and conformation of lysozyme is independent of pH of the elution buffer. The microscopic information from simulation shows that protein unfolding is mainly observed on the adsorbent surface, but it cannot be detected in the experiments that only probe the proteins in the bulk liquid.     

Predicting Atomic Details of the Unfolding Pathway for YibK,a Knotted Protein from the SPOUT Superfamily

Abstract

Several protein structures have been reported to contain intricate knots of the polypeptide backbone but the mechanism of the (un)folding process of knotted proteins remains unknown. The members of the SPOUT superfamily of RNA methyltransferases are some of the most intensely studied systems for investigation of the knot formation and function. YibK (whose biochemical function remains unknown) is the representative protein of the SPOUT superfamily. This protein exhibits a deep trefoil knot at the C-terminus.

We conducted an extensive computational analysis of the unfolding process for the monomeric form of YibK. In order to predict the (un)folding pathway of YibK, we have calculated the order of secondary structure disassembly using UNFOLD, and performed thermal unfolding simulations using classical Molecular Dynamics (MD), as well as simulations employing reduced representation of the peptide chain using either MD with the UNRES method or the Monte Carlo (MC) unfolding with the REFINER method.

Results obtained from all methods used in this work are in qualitative agreement. We found that YibK unfolds through four intermediate states. The trefoil knot in YibK disappears at the end of the unfolding process, long after the protein loses its native topology. We observed that the C-terminus leaves the knotting loop folded into a hairpin-like structure, in agreement with the results of coarse-grained simulation reported earlier. We propose that the folding pathway of YibK corresponds to the reversed sequence of events observed in the unfolding pathway elucidated in this study. Thus, we predict that the knot formation is the slowest part of the YibK folding process.     

Thermal unfolding simulations of apo-calmodulin using leap-dynamics

The simulation method leap-dynamics (LD) has been applied to protein thermal unfolding simulations to investigate domain-specific unfolding behavior. Thermal unfolding simulations of the 148-residue protein apo-calmodulin with implicit solvent were performed at temperatures 290 K, 325 K, and 360 K and compared with the corresponding molecular dynamics trajectories in terms of a number of calculated conformational parameters. The main experimental results of unfolding are reproduced in showing the lower stability of the C-domain: at 290 K, both the N- and C-domains are essentially stable; at 325 K, the C-domain unfolds, whereas the N-domain remains folded; and at 360 K, both domains unfold extensively. This behavior could not be reproduced by molecular dynamics simulations alone under the same conditions. These results show an encouraging degree of convergence between experiment and LD simulation. The simulations are able to describe the overall plasticity of the apo-calmodulin structure and to reveal details such as reversible folding/unfolding events within single helices. The results show that by using the combined application of a fast and efficient sampling routine with a detailed molecular dynamics force field, unfolding simulations of proteins at atomic resolution are within the scope of current computational power.     

Structure and stability of recombinant protein depend on the extra N-terminal methionine residue: S6 permutein from direct and fusion expression systems.

Two permuted variants of S6 ribosomal protein were obtained in direct and fusion expression systems, respectively. The product of direct expression contained the extra N-terminal methionine residue. The structural properties and conformational stability of these permuteins were compared using 1-D (1)H-NMR, circular dichroism, intrinsic fluorescence, differential scanning calorimetry and resistance to urea-induced unfolding. A pronounced difference in all the parameters studied has been demonstrated. This means that the structure of recombinant protein can be sensitive to peculiarities of the expression and purification procedures, leading particularly to the presence or absence of the Met at the first position in the target protein sequence.     

Dynamics and unfolding pathways of a hyperthermophilic and a mesophilic rubredoxin.

Molecular dynamics simulations in solution are performed for a rubredoxin from the hyperthermophilic archaeon Pyrococcus furiosus (RdPf) and one from the mesophilic organism Desulfovibrio vulgaris (RdDv). The two proteins are simulated at four temperatures: 300 K, 373 K, 473 K (two sets), and 500 K; the various simulations extended from 200 ps to 1,020 ps. At room temperature, the two proteins are stable, remain close to the crystal structure, and exhibit similar dynamic behavior; the RMS residue fluctuations are slightly smaller in the hyperthermophilic protein. An analysis of the average energy contributions in the two proteins is made; the results suggest that the intraprotein energy stabilizes RdPf relative to RdDv. At 373 K, the mesophilic protein unfolds rapidly (it begins to unfold at 300 ps), whereas the hyperthermophilic does not unfold over the simulation of 600 ps. This is in accord with the expected stability of the two proteins. At 473 K, where both proteins are expected to be unstable, unfolding behavior is observed within 200 ps and the mesophilic protein unfolds faster than the hyperthermophilic one. At 500 K, both proteins unfold; the hyperthermophilic protein does so faster than the mesophilic protein. The unfolding behavior for the two proteins is found to be very similar. Although the exact order of events differs from one trajectory to another, both proteins unfold first by opening of the loop region to expose the hydrophobic core. This is followed by unzipping of the beta-sheet. The results obtained in the simulation are discussed in terms of the factors involved in flexibility and thermostability.